Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complementary”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

Programmable soft valves for digital and analog control

In soft devices, complex actuation sequences and precise force control typically require hard electronic valves and microcontrollers. Existing designs for entirely soft pneumatic control systems are capable of either digital or analog operation, but not both, and are limited by speed of actuation, range of pressure, time required for fabrication, or loss of power through pull-down resistors. Using the nonlinear mechanics intrinsic to structures composed of soft materials—in this case, by leveraging membrane inversion and tube kinking—two modular soft components are developed: a piston actuator and a bistable pneumatic switch. These two components combine to create valves capable of analog pressure regulation, simplified digital logic, controlled oscillation, nonvolatile memory storage, linear actuation, and interfacing with human users in both digital and analog formats. Three demonstrations showcase the capabilities of systems constructed from these valves: 1) a wearable glove capable of analog control of a soft artificial robotic hand based on input from a human user’s fingers, 2) a human-controlled cushion matrix designed for use in medical care, and 3) an untethered robot which travels a distance dynamically programmed at the time of operation to retrieve an object. This work illustrates pathways for complementary digital and analog control of soft robots using a unified valve design.

42 ENGINEERING↗

Deciphering the evolution of flavin-dependent monooxygenase stereoselectivity using ancestral sequence reconstruction

Controlling the selectivity of a reaction is critical for target-oriented synthesis. Accessing complementary selectivity profiles enables divergent synthetic strategies, but is challenging to achieve in biocatalytic reactions given enzymes’ innate preferences of a single selectivity. Thus, it is critical to understand the structural features that control selectivity in biocatalytic reactions to achieve tunable selectivity. Here, we investigate the structural features that control the stereoselectivity in an oxidative dearomatization reaction that is key to making azaphilone natural products. Crystal structures of enantiocomplementary biocatalysts guided the development of multiple hypotheses centered on the structural features that control the stereochemical outcome of the reaction; however, in many cases, direct substitutions of active site residues in natural proteins led to inactive enzymes. Ancestral sequence reconstruction (ASR) and resurrection were employed as an alternative strategy to probe the impact of each residue on the stereochemical outcome of the dearomatization reaction. These studies suggest that two mechanisms are active in controlling the stereochemical outcome of the oxidative dearomatization reaction: one involving multiple active site residues in AzaH and the other dominated by a single Phe to Tyr switch in TropB and AfoD. Moreover, this study suggests that the flavin-dependent monooxygenases (FDMOs) adopt simple and flexible strategies to control stereoselectivity, which has led to stereocomplementary azaphilone natural products produced by fungi. This paradigm of combining ASR and resurrection with mutational and computational studies showcases sets of tools for understanding enzyme mechanisms and provides a solid foundation for future protein engineering efforts.

Science & Technology - Other Topics↗

Decoupling of catalysis and transition state analog binding from mutations throughout a phosphatase revealed by high-throughput enzymology

Using high-throughput microfluidic enzyme kinetics (HT-MEK), we measured over 9,000 inhibition curves detailing impacts of 1,004 single-site mutations throughout the alkaline phosphatase PafA on binding affinity for two transition state analogs (TSAs), vanadate and tungstate. As predicted by catalytic models invoking transition state complementary, mutations to active site and active-site-contacting residues had highly similar impacts on catalysis and TSA binding. Unexpectedly, most mutations to more distal residues that reduced catalysis had little or no impact on TSA binding and many even increased tungstate affinity. These disparate effects can be accounted for by a model in which distal mutations alter the enzyme’s conformational landscape, increasing the occupancy of microstates that are catalytically less effective but better able to accommodate larger transition state analogs. In support of this ensemble model, glycine substitutions (rather than valine) were more likely to increase tungstate affinity (but not more likely to impact catalysis), presumably due to increased conformational flexibility that allows previously disfavored microstates to increase in occupancy. These results indicate that residues throughout an enzyme provide specificity for the transition state and discriminate against analogs that are larger only by tenths of an Ångström. Thus, engineering enzymes that rival the most powerful natural enzymes will likely require consideration of distal residues that shape the enzyme’s conformational landscape and fine-tune active-site residues. Biologically, the evolution of extensive communication between the active site and remote residues to aid catalysis may have provided the foundation for allostery to make it a highly evolvable trait.

59 BASIC BIOLOGICAL SCIENCES↗

Correlating chemistry and mass transport in sustainable iron production

Steelmaking contributes 8% to the total CO 2 emissions globally, primarily due to coal-based iron ore reduction. Clean hydrogen-based ironmaking has variable performance because the dominant gas–solid reduction mechanism is set by the defects and pores inside the mm- to nm-sized oxide particles that change significantly as the reaction progresses. While these governing dynamics are essential to establish continuous flow of iron and its ores through reactors, the direct link between agglomeration and chemistry is still contested due to missing measurements. In this work, we directly measure the connection between chemistry and agglomeration in the smallest iron oxides relevant to magnetite ores. Using synthesized spherical 10-nm magnetite particles reacting in H 2 , we resolve the formation and consumption of wüstite (Fe 1-x O)—the step most commonly attributed to whiskering. Using X-ray diffraction, we resolve crystallographic anisotropy in the rate of the initial reaction. Complementary imaging demonstrated how the particles self-assemble, subsequently react, and grow into elongated “whisker” structures. Our insights into how morphologically uniform iron oxide particles react and agglomerate in H 2 reduction enable future size-dependent models to effectively describe the multiscale aspects of iron ore reduction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The exceptional form and function of the giant bacterium Ca. Epulopiscium viviparus revolves around its sodium motive force

Epulopiscium spp. are the largest known heterotrophic bacteria; a large cigar-shaped individual is a million times the volume of Escherichia coli. To better understand the metabolic potential and relationship of Epulopiscium sp. type B with its host Naso tonganus, we generated a high-quality draft genome from a population of cells taken from a single fish. We propose the name Candidatus Epulopiscium viviparus to describe populations of this best-characterized Epulopiscium species. Metabolic reconstruction reveals more than 5% of the genome codes for carbohydrate active enzymes, which likely degrade recalcitrant host-diet algal polysaccharides into substrates that may be fermented to acetate, the most abundant short-chain fatty acid in the intestinal tract. Moreover, transcriptome analyses and the concentration of sodium ions in the host intestinal tract suggest that the use of a sodium motive force (SMF) to drive ATP synthesis and flagellar rotation is integral to symbiont metabolism and cellular biology. In natural populations, genes encoding both F-type and V-type ATPases and SMF generation via oxaloacetate decarboxylation are among the most highly expressed, suggesting that ATPases synthesize ATP and balance ion concentrations across the cell membrane. High expression of these and other integral membrane proteins may allow for the growth of its extensive intracellular membrane system. Further, complementary metabolism between microbe and host is implied with the potential provision of nitrogen and B vitamins to reinforce this nutritional symbiosis. The few features shared by all bacterial behemoths include extreme polyploidy, polyphosphate synthesis, and thus far, they have all resisted cultivation in the lab.

59 BASIC BIOLOGICAL SCIENCES↗

Spatially resolved land and grid model of carbon neutrality in China

China has committed to achieve net carbon neutrality by 2060 to combat global climate change, which will require unprecedented deployment of negative emissions technologies, renewable energies (RE), and complementary infrastructure. At terawatt-scale deployment, land use limitations interact with operational and economic features of power systems. To address this, we developed a spatially resolved resource assessment and power systems planning optimization that models a full year of power system operations, sub-provincial RE siting criteria, and transmission connections. Our modeling results show that wind and solar must be expanded to 2,000 to 3,900 GW each, with one plausible pathway leading to 300 GW/yr combined annual additions in 2046 to 2060, a three-fold increase from today. Over 80% of solar and 55% of wind is constructed within 100 km of major load centers when accounting for current policies regarding land use. Large-scale low-carbon systems must balance key trade-offs in land use, RE resource quality, grid integration, and costs. Under more restrictive RE siting policies, at least 740 GW of distributed solar would become economically feasible in regions with high demand, where utility-scale deployment is limited by competition with agricultural land. Effective planning and policy formulation are necessary to achieve China’s climate goals.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Prevention efficacy of the broadly neutralizing antibody VRC01 depends on HIV-1 envelope sequence features

In the Antibody Mediated Prevention (AMP) trials (HVTN 704/HPTN 085 and HVTN 703/HPTN 081), prevention efficacy (PE) of the monoclonal broadly neutralizing antibody (bnAb) VRC01 (vs. placebo) against HIV-1 acquisition diagnosis varied according to the HIV-1 Envelope (Env) neutralization sensitivity to VRC01, as measured by 80% inhibitory concentration (IC80). Here, we performed a genotypic sieve analysis, a complementary approach to gaining insight into correlates of protection that assesses how PE varies with HIV-1 sequence features. We analyzed HIV-1 Env amino acid (AA) sequences from the earliest available HIV-1 RNA-positive plasma samples from AMP participants diagnosed with HIV-1 and identified Env sequence features that associated with PE. The strongest Env AA sequence correlate in both trials was VRC01 epitope distance that quantifies the divergence of the VRC01 epitope in an acquired HIV-1 isolate from the VRC01 epitope of reference HIV-1 strains that were most sensitive to VRC01-mediated neutralization. In HVTN 704/HPTN 085, the Env sequence-based predicted probability that VRC01 IC80 against the acquired isolate exceeded 1 µg/mL also significantly associated with PE. In HVTN 703/HPTN 081, a physicochemical-weighted Hamming distance across 50 VRC01 binding-associated Env AA positions of the acquired isolate from the most VRC01-sensitive HIV-1 strain significantly associated with PE. These results suggest that incorporating mutation scoring by BLOSUM62 and weighting by the strength of interactions at AA positions in the epitope:VRC01 interface can optimize performance of an Env sequence-based biomarker of VRC01 prevention efficacy. Future work could determine whether these results extend to other bnAbs and bnAb combinations.

60 APPLIED LIFE SCIENCES↗

Reexpansion of charged nanoparticle assemblies in concentrated electrolytes

Electrostatic forces in solutions are highly relevant to a variety of fields, ranging from electrochemical energy storage to biology. However, their manifestation in concentrated electrolytes is not fully understood, as exemplified by counterintuitive observations of colloidal stability and long-ranged repulsions in molten salts. Highly charged biomolecules, such as DNA, respond sensitively to ions in dilute solutions. Here, we use non-base-pairing DNA-coated nanoparticles (DNA-NP) to analyze electrostatic interactions in concentrated salt solutions. Despite their negative charge, these conjugates form colloidal crystals in solutions of sufficient divalent cation concentration. We utilize small-angle X-ray scattering (SAXS) to study such DNA-NP assemblies across the full accessible concentration ranges of aqueous CaCl 2 , MgCl 2 , and SrCl 2 solutions. SAXS shows that the crystallinity and phases of the assembled structures vary with cation type. For all tested salts, the aggregates contract with added ions at low salinities and then begin expanding above a cation-dependent threshold salt concentration. Wide-angle X-ray scattering (WAXS) reveals enhanced positional correlations between ions in the solution at high salt concentrations. Complementary molecular dynamics simulations show that these ion–ion interactions reduce the favorability of dense ion configurations within the DNA brushes below that of the bulk solution. Measurements in solutions with lowered permittivity demonstrate a simultaneous increase in ion coupling and decrease in the concentration at which aggregate expansion begins, thus confirming the connection between these phenomena. Our work demonstrates that interactions between charged objects continue to evolve considerably into the high-concentration regime, where classical theories project electrostatics to be of negligible consequence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hopping and crawling DNA-coated colloids

Understanding the motion of particles with multivalent ligand-receptors is important for biomedical applications and material design. Yet, even among a single design, the prototypical DNA-coated colloids, seemingly similar micrometric particles hop or roll, depending on the study. We shed light on this problem by observing DNA-coated colloids diffusing near surfaces coated with complementary strands for a wide array of coating designs. We find colloids rapidly switch between 2 modes: They hop—with long and fast steps—and crawl—with short and slow steps. Both modes occur at all temperatures around the melting point and over various designs. The particles become increasingly subdiffusive as temperature decreases, in line with subsequent velocity steps becoming increasingly anticorrelated, corresponding to switchbacks in the trajectories. Overall, crawling (or hopping) phases are more predominant at low (or high) temperatures; crawling is also more efficient at low temperatures than hopping to cover large distances. We rationalize this behavior within a simple model: At lower temperatures, the number of bound strands increases, and detachment of all bonds is unlikely, hence, hopping is prevented and crawling favored. We thus reveal the mechanism behind a common design rule relying on increased strand density for long-range self-assembly: Dense strands on surfaces are required to enable crawling, possibly facilitating particle rearrangements.

Science & Technology - Other Topics↗

Atomic view of photosynthetic metabolite permeability pathways and confinement in synthetic carboxysome shells

Carboxysomes are protein microcompartments found in cyanobacteria, whose shell encapsulates rubisco at the heart of carbon fixation in the Calvin cycle. Carboxysomes are thought to locally concentrate CO 2 in the shell interior to improve rubisco efficiency through selective metabolite permeability, creating a concentrated catalytic center. However, permeability coefficients have not previously been determined for these gases, or for Calvin-cycle intermediates such as bicarbonate (HCO$^–_3$), 3-phosphoglycerate, or ribulose-1,5-bisphosphate. Starting from a high-resolution cryogenic electron microscopy structure of a synthetic β-carboxysome shell, we perform unbiased all-atom molecular dynamics to track metabolite permeability across the shell. The synthetic carboxysome shell structure, lacking the bacterial microcompartment trimer proteins and encapsulation peptides, is found to have similar permeability coefficients for multiple metabolites, and is not selectively permeable to HCO$^–_3$ relative to CO 2 . To resolve how these comparable permeabilities can be reconciled with the clear role of the carboxysome in the CO 2 -concentrating mechanism in cyanobacteria, complementary atomic-resolution Brownian Dynamics simulations estimate the mean first passage time for CO 2 assimilation in a crowded model carboxysome. Despite a relatively high CO 2 permeability of approximately 10 -2 cm/s across the carboxysome shell, the shell proteins reflect enough CO 2 back toward rubisco that 2,650 CO 2 molecules can be fixed by rubisco for every 1 CO 2 molecule that escapes under typical conditions. The permeabilities determined from all-atom molecular simulation are key inputs into flux modeling, and the insight gained into carbon fixation can facilitate the engineering of carboxysomes and other bacterial microcompartments for multiple applications.

59 BASIC BIOLOGICAL SCIENCES↗

Dissecting neurofilament tail sequence-phosphorylation-structure relationships with multicomponent reconstituted protein brushes

Neurofilaments (NFs) are multisubunit, bottlebrush-shaped intermediate filaments abundant in the axonal cytoskeleton. Each NF subunit contains a long intrinsically disordered tail domain, which protrudes from the NF core to form a “brush” surrounding each NF. Precisely how the tails’ variable charge patterns and repetitive phosphorylation sites mediate their conformation within the brush remains an open question in axonal biology. We address this problem by grafting recombinant NF tail protein constructs NF-Light, -Medium, and -Heavy (NFL, NFM, and NFH) to surfaces, yielding protein brushes of defined stoichiometry that can be phosphorylated in vitro. Atomic force microscopy measurements reveal that brush height depends on composition monotonically but not always linearly for binary NFL:NFM or NFL:NFH systems, and that NFM-based brushes are highly extended, while brushes incorporating the much larger NFH are surprisingly compact even after multisite phosphorylation. Complementary self-consistent field theory (SCFT) predicts multilayer brush morphologies for NFM and phosphorylated NFH brushes. Further experiments and SCFT analysis with designed mutants reveal that N-terminal negative charges in the NFH tail repel phosphorylated residues to generate the multilayer morphology, while the C-terminal charge-neutral region contributes to multilayer brush morphology but not total brush height. Charge-shuffled NFM variants show that charge segregation promotes brush collapse near physiological ionic strengths. Collectively, this study supports a role for NFM in establishing a dynamic range for NF brush conformation, lending insight into previous in vitro and in vivo findings. More broadly, this work establishes a platform for dissecting contributions of disordered protein sequence to conformation at interfaces.

Science & Technology - Other Topics↗

Prioritizing urban heat adaptation infrastructure based on multiple outcomes: Comfort, health, and energy

Globally, cities face increasing extreme heat, impacting comfort, health, and energy consumption. Infrastructure-based heat adaptation strategies can improve these outcomes, but each strategy has a unique mix of benefits and drawbacks. Here, we apply an urbanized meteorological model (WRF) with the newly integrated multilayer BEP-Tree street tree model to dynamically downscale Earth System Model projections and a 3-D microclimate model (TUF-Pedestrian) to simulate the street-scale radiation environment impacting pedestrians. We evaluate the performance of five heat adaptation strategies (street trees, cool roofs, green roofs, rooftop photovoltaics (PV), and reflective pavements) during extreme heat events in three cities with contrasting background climates (Toronto, Phoenix, and Miami), under contemporary and end-of-century projected climates, based on three metrics: outdoor heat stress, air conditioning (AC) energy use, and ventilation of vehicular air pollution. No single adaptation strategy improves all three outcomes. While street trees inhibit ventilation, they reduce outdoor heat stress four times more effectively than the next best strategy via shade provision, fully offsetting heat stress increases under a high-emissions end-of-century climate scenario in all cities studied. Cool roofs and green roofs moderately reduce heat stress and energy use. Alternatively, rooftop PV with energy storage can generate sufficient power for space cooling but have marginal effects on heat stress. Reflective pavements are the least effective across metrics. Where the ventilation of street-level emissions is of less concern, our results clearly support the combination of street trees and rooftop PV as a highly complementary and effective means of adaptive mitigation across different climates and neighborhood densities.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Conformation-specific synthetic intrabodies modulate mTOR signaling with subcellular spatial resolution

Subcellular compartmentalization is integral to the spatial regulation of mechanistic target of rapamycin (mTOR) signaling. However, the biological outputs associated with location-specific mTOR signaling events are poorly understood and challenging to decouple. Here, we engineered synthetic intracellular antibodies (intrabodies) that are capable of modulating mTOR signaling with genetically programmable spatial resolution. Epitope-directed phage display was exploited to generate high affinity synthetic antibody fragments (Fabs) against the FKBP12–Rapamycin binding site of mTOR (mTOR FRB ). We determined high-resolution crystal structures of two unique Fabs that discriminate distinct conformational states of mTOR FRB through recognition of its substrate recruitment interface. By leveraging these conformation-specific binders as intracellular probes, we uncovered the structural basis for an allosteric mechanism governing mTOR complex 1 (mTORC1) stability mediated by subtle structural adjustments within mTOR FRB . Furthermore, our results demonstrated that synthetic binders emulate natural substrates by employing divergent yet complementary hydrophobic residues at defined positions, underscoring the broad molecular recognition capability of mTOR FRB . Intracellular signaling studies showed differential time-dependent inhibition of S6 kinase 1 and Akt phosphorylation by genetically encoded intrabodies, thus supporting a mechanism of inhibition analogous to the natural product rapamycin. Finally, we implemented a feasible approach to selectively modulate mTOR signaling in the nucleus through spatially programmed intrabody expression. These findings establish intrabodies as versatile tools for dissecting the conformational regulation of mTORC1 and should be useful to explore how location-specific mTOR signaling influences disease progression.

Science & Technology - Other Topics↗

An enzyme-based approach for highly efficient self-replication of DNA origami dimers

Self-replication and exponential growth are essential to all living things, the driving force for Darwinian evolution, and potentially useful in nanotechnology for large-scale production of nanoscopic materials. An artificial (nonliving) self-replication system has been shown to exhibit exponential growth and selection using DNA monomer origami tiles templated on a dimer seed. That system purposefully avoided the use of enzymes to get a hint of how self-replication might have evolved in a prebiotic world by using CNV K and UV light to crosslink complementary DNA single strands. For further investigations into competition and extinction and for potential applications involving biocompatibility, we wanted to investigate enzymatic ligation to replace the chemical photo crosslinking step. Here, we present a system which uses thermotolerant T4 DNA ligase and no UV. This system has several additional advantages including a much faster cycling time, yielding 2,000,000 amplifications in 12 h. We also introduce competition to study the possibility of Darwinian-like evolution. Two pairs of DNA origami tiles compete for the same connection strands and show different growth rates under different connection strand concentrations. This system has the potential to combine with other enzymes, such as RNA polymerase to support feedback, allowing us to fine-tune replication dynamics and achieve sophisticated, life-like behaviors. The highly efficient self-replication and exponential growth of DNA origami dimers demonstrated in this work not only enhances our understanding of Darwinian evolution in nature but also opens the door to applications ranging from synthetic biology to smart materials.

Science & Technology - Other Topics↗

Designing for cooperative grain boundary segregation in multicomponent alloys

Tailoring the nanoscale distribution of chemical species at grain boundaries is a powerful method to dramatically influence the properties of polycrystalline materials. However, classical approaches to the problem have tacitly assumed that only competition is possible between solute species. In this paper, we show that solute elements can cooperate in the way they segregate to grain boundaries: In properly targeted alloys, the different chemical species cooperate to each fill complementary grain boundary sites disfavored by the other. By developing a theoretical “spectral” approach to this problem based on quantum-accurate grain boundary site distributions, we show how grain boundaries can be cooperatively alloyed, whether by depletion or enrichment. We provide machine-learned cosegregation information for over 700 ternary aluminum-based alloys and experimentally validate the concept in one ternary alloy where cosegregation is not expected by prior models but is expected based on the cooperative model.

grain boundaries↗

Hund’s flat band in a frustrated spinel oxide

Electronic flat bands associated with quenched kinetic energy and heavy electron mass have attracted great interest for promoting strong electronic correlations and emergent phenomena such as high-temperature charge fractionalization and superconductivity. Intense experimental and theoretical research has been devoted to establishing the rich nontrivial metallic and heavy fermion phases intertwined with such localized electronic states. Here, in this study, we investigate the transition metal oxide spinel LiV 2 O 4 , an enigmatic heavy fermion compound lacking localized f orbital states. We use angle-resolved photoemission spectroscopy and dynamical mean-field theory to reveal a kind of correlation-induced flat band with suppressed interatomic electron hopping arising from intra-atomic Hund’s coupling. The appearance of heavy quasiparticles is ascribed to a proximate orbital-selective Mott state characterized by fluctuating local moments as evidenced by complementary magnetotransport measurements. The spectroscopic fingerprints of long-lived quasiparticles and their disappearance with increasing temperature further support the emergence of a high-temperature “bad” metal state observed in transport data. This work resolves a long-standing puzzle on the origin of heavy fermion behavior and unconventional transport in LiV 2 O 4 . Simultaneously, it opens a path to achieving flat bands through electronic interactions in d-orbital systems with geometrical frustration, potentially enabling the realization of exotic phases of matter such as the fractionalized Fermi liquids.

36 MATERIALS SCIENCE↗

Tandem sialoglycan-binding modules in a Streptococcus sanguinis serine-rich repeat adhesin create target dependent avidity effects

Sialic acid–binding immunoglobulin-like lectins (Siglec)–like domains of streptococcal serine-rich repeat (SRR) adhesins recognize sialylated glycans on human salivary, platelet, and plasma glycoproteins via a YTRY sequence motif. The SRR adhesin from Streptococcus sanguinis strain SK1 has tandem sialoglycan-binding domains and has previously been shown to bind sialoglycans with high affinity. However, both domains contain substitutions within the canonical YTRY motif, making it unclear how they interact with host receptors. To identify how the S. sanguinis strain SK1 SRR adhesin affects interactions with sialylated glycans and glycoproteins, we determined high-resolution crystal structures of the binding domains alone and with purified trisaccharides. These structural studies determined that the ligands still bind at the noncanonical binding motif, but with fewer hydrogen-bonding interactions to the protein than is observed in structures of other Siglec-like adhesins. Complementary biochemical studies identified that each of the two binding domains has a different selectivity profile. Interestingly, the binding of SK1 to platelets and plasma glycoproteins identified that the interaction to some host targets is dominated by the contribution of one binding domain, whereas the binding to other host receptors is mediated by both binding domains. These results provide insight into outstanding questions concerning the roles of tandem domains in targeting host receptors and suggest mechanisms for how pathogens can adapt to the availability of a range of related but nonidentical host receptors. They further suggest that the definition of the YTRY motif should be changed to ΦTRX, a more rigorous description of this sialic acid–recognition motif given recent findings.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling bicarbonate formation in an alkaline solution with multi-level quantum mechanics/molecular dynamics simulations

Understanding carbonate speciation and how it may be modulated is essential for the advancement of carbon dioxide (CO 2 ) capture and storage technologies, which often rely on the transformation of CO 2 into carbonate, e.g. via the formation of carbonate minerals. To date, few atomic-level, quantum-mechanics-based simulations have been carried out to characterize how carbonic acid (H 2 CO 3 ) and bicarbonate ($HCO^{-}_{3}$) form in aqueous solution, and how pH affects this process. Recently, Martirez and Carter utilized rare-event sampling density functional theory molecular dynamics simulations in combination with multi-level embedded correlated wavefunction theory, thus accounting for both solvent dynamics and electron correlation accurately, to elucidate the mechanism of H 2 CO 3 formation in neutral solution (J. Am. Chem. Soc., 145, 12561, 2023). Here, we perform a complementary simulation using the same method to map out the energetics of $HCO^{-}_{3}$ formation from dissolved CO 2 in basic solution. We find that, as in H 2 CO 3 formation, including water dynamics is important to obtain an accurate prediction of the energetics for the aforementioned reaction. Furthermore, only with MD did we identify the correct pathway for the reaction, in which water – not hydroxide – acts as the initial nucleophile and only at the transition state does it lose a proton.

74 ATOMIC AND MOLECULAR PHYSICS↗