Search NASA⌕ Search

SEARCH · Search NASA

Results for “biochemical simulations”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Structural Characterization and Dynamics of AdhE Ultrastructures from Clostridium thermocellum Show a Containment Strategy for Toxic Intermediates

Clostridium thermocellum, a cellulolytic thermophilic anaerobe, is considered by many to be a prime candidate for the realization of consolidated bioprocessing (CBP) and is known as an industry standard for biofuel production. C. thermocellum is among the best biomass degraders identified to date in nature and produces ethanol as one of its main products. Many studies have helped increase ethanol titers in this microbe; however, ethanol production using C. thermocellum is still not economically viable. Therefore, a better understanding of its ethanol synthesis pathway is required. The main pathway for ethanol production in C. thermocellum involves the bifunctional aldehyde-alcohol dehydrogenase (AdhE). To better understand the function of the C. thermocellum AdhE, we used cryo-electron microscopy (cryo-EM) to obtain a 3.28 A structure of the AdhE complex. This high-resolution structure, in combination with molecular dynamics simulations, provides insight into the substrate channeling of the toxic intermediate acetaldehyde, indicates the potential role of C. thermocellum AdhE to regulate activity and cofactor pools, and establishes a basis for future engineering studies. The containment strategy found in this enzyme offers a template that could be replicated in other systems where toxic intermediates need to be sequestered to increase the production of valuable biochemicals.

09 BIOMASS FUELS↗
DOE OSTI2024

Identifying Bayesian optimal experiments for uncertain biochemical pathway models

Abstract Pharmacodynamic (PD) models are mathematical models of cellular reaction networks that include drug mechanisms of action. These models are useful for studying predictive therapeutic outcomes of novel drug therapies in silico. However, PD models are known to possess significant uncertainty with respect to constituent parameter data, leading to uncertainty in the model predictions. Furthermore, experimental data to calibrate these models is often limited or unavailable for novel pathways. In this study, we present a Bayesian optimal experimental design approach for improving PD model prediction accuracy. We then apply our method using simulated experimental data to account for uncertainty in hypothetical laboratory measurements. This leads to a probabilistic prediction of drug performance and a quantitative measure of which prospective laboratory experiment will optimally reduce prediction uncertainty in the PD model. The methods proposed here provide a way forward for uncertainty quantification and guided experimental design for models of novel biological pathways.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Spatial access and resource limitations control carbon mineralization in soils

Core-scale soil carbon fluxes are ultimately regulated by pore-scale dynamics of substrate availability and microbial access. These are constrained by physicochemical and biochemical phenomena (e.g. spatial access and hydrologic connectivity, physical occlusion, adsorption-desorption with mineral surfaces, nutrient and resource limitations). We conducted an experiment to determine how spatial access and resource limitations influence core-scale SOM mineralization, and how these are regulated by antecedent moisture conditions. Intact soil cores were incubated at field-moist vs. drought conditions, after which they were saturated from above (to simulate precipitation) or below (to simulate groundwater recharge). Soluble C (acetate) and N (nitrate) forms were added to some cores during the rewetting process to alleviate potential nutrient limitations. Respiration showed evidence of C-limitation, as C amendments increased the cumulative CO 2 evolved. Drought and rewetting increased soil respiration, and these soils also exhibited increased complex aromatic molecules in porewater. This newly available substrate appeared to alleviate nutrient limitations on respiration, because there were no further respiration increases with subsequent C and N amendments. We expected that respiration would be influenced by wetting direction, as simulated precipitation would mobilize C from the surface. However, as a main effect, this response was seen only in the C-amended soils, indicating that surface-C may not have been bioavailable. At the pore scale (porewater samples), compared to the baseline soils, drought and the C, N amendments caused a net loss of identified molecules when the soils were rewet from below, whereas wetting from above caused a net increase in identified molecules. This indicates that as soils were wet from below, the fresh inputs simulated the C-and N-limited microbial populations present deeper in the soil profile. Furthermore, our experiment highlights the complex and interactive role of antecedent moisture conditions, wetting direction, and resource limitations in driving core-scale C fluxes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A subcellular biochemical model for T6SS dynamics reveals winning competitive strategies

Abstract The type VI secretion system (T6SS) is a broadly distributed interbacterial weapon that can be used to eliminate competing bacterial populations. Although unarmed target populations are typically used to study T6SS function in vitro, bacteria most likely encounter other T6SS-armed competitors in nature. However, the connection between subcellular details of the T6SS and the outcomes of such mutually lethal battles is not well understood. Here, we incorporate biological data derived from natural competitors of Vibrio fischeri light organ symbionts to build a biochemical model for T6SS at the single-cell level, which we then integrate into an agent-based model (ABM). Using the ABM, we isolate and experiment with strain-specific physiological differences between competitors in ways not possible with biological samples to identify winning strategies for T6SS-armed populations. Through in vitro experiments, we discover that strain-specific differences exist in T6SS activation speed. ABM simulations corroborate that faster activation is dominant in determining survival during competition. Once competitors are fully activated, the energy required for T6SS creates a tipping point where increased weapon building and firing becomes too costly to be advantageous. Through ABM simulations, we identify the threshold where this transition occurs in the T6SS parameter space. We also find that competitive outcomes depend on the geometry of the battlefield: unarmed target cells survive at the edges of a range expansion where unlimited territory can be claimed. Alternatively, competitions within a confined space, much like the light organ crypts where natural V. fischeri compete, result in the rapid elimination of the unarmed population.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Sensitivity of Riparian Buffer Designs to Climate Change—Nutrient and Sediment Loading to Streams: A Case Study in the Albemarle-Pamlico River Basins (USA) Using HAWQS

Riparian buffer zones (RBZs) provide multiple benefits to watershed ecosystems. We aimed to conduct an extensive sensitivity analysis of the RBZ designs to climate change nutrient and sediment loadings to streams. We designed 135 simulation scenarios starting with the six baselines RBZs (grass, urban, two-zone forest, three-zone forest, wildlife, and naturalized) in three 12-digit Hydrologic Unit Code watersheds within the Albemarle-Pamlico river basin (USA). Using the hydrologic and water quality system (HAWQS), we assessed the sensitivity of the designs to five water quality indicator (WQI) parameters: dissolved oxygen (DO), total phosphorous (TP), total nitrogen (TN), sediment (SD), and biochemical oxygen demand (BD). To understand the climate mitigation potential of RBZs, we identified a subset of future climate change projection models of air temperature and precipitation using EPA’s Locating and Selecting Scenarios Online tool. Analyses revealed optimal RBZ designs for the three watersheds. In terms of watershed ecosystem services sustainability, the optimal Urban RBZ in contemporary climate (1983–2018) reduced SD from 61–96%, TN from 34–55%, TP from 9–48%, and BD from 53–99%, and raised DO from 4–10% with respect to No-RBZ in the three watersheds. The late century’s (2070–2099) extreme mean annual climate changes significantly increased the projected SD and BD; however, the addition of urban RBZs was projected to offset the climate change reducing SD from 28–94% and BD from 69–93% in the watersheds. All other types of RBZs are also projected to fully mitigate the climate change impacts on WQI parameters except three-zone RBZ.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Virtual Engineering: Python framework for engineering process design

Virtual Engineering (VE) is a Python software framework designed to accelerate the research and development of engineering processes that are fundamentally defined by multiple unit operations executed in series. VE supports a wide variety of different multi-physics models and integrates them to simulate a complete end-to-end process. To automate the execution of this model sequence, VE provides (i) a robust method to communicate between models, (ii) a high-level, user-friendly interface to set model parameters and enable optimization, and (iii) an overall model-agnostic approach that allows new computational units to be swapped in and out of workflows. Although the VE framework was developed to support the biochemical conversion of biomass to fuel, we have designed each component to easily accommodate new domains and unit models.

09 BIOMASS FUELS↗
DOE OSTI2021

Data for: Spatial access and resource limitations control carbon mineralization in soils

This dataset contains data and code used for the paper "Spatial access and resource limitations control carbon mineralization in soils", https://doi.org/10.1016/j.soilbio.2021.108427. Core-scale soil carbon fluxes are ultimately regulated by pore-scale dynamics of substrate availability and microbial access. These are constrained by physicochemical and biochemical phenomena (e.g. spatial access and hydrologic connectivity, physical occlusion, adsorption-desorption with mineral surfaces, nutrient and resource limitations). We conducted an experiment to determine how spatial access and resource limitations influence core-scale water-soluble soil organic matter (SOM) mineralization, and how these are regulated by antecedent moisture conditions. Intact soil cores were incubated at field-moist vs. drought conditions, after which they were saturated from above (to simulate precipitation) or below (to simulate groundwater recharge). Soluble carbon (acetate) and nitrogen (nitrate) forms were added to some cores during the rewetting process to alleviate potential nutrient limitations. Soil respiration was measured during the incubation, after which pore water was extracted from the saturated soils and analyzed for water soluble organic carbon concentrations and characterization. Our results showed that carbon (C) amendments increased the cumulative carbon dioxide (CO2) evolved from the soil cores, suggesting that the soils were C-limited. Drought and rewetting increased soil respiration, and there was a greater abundance of complex aromatic molecules in pore waters sampled from these soils. This newly available substrate appeared to alleviate nutrient limitations on respiration, because there were no further respiration increases with subsequent C and N amendments. We had hypothesized that respiration would be influenced by wetting direction, as simulated precipitation would mobilize C from the surface. However, as a main effect, this response was seen only in the C-amended soils, indicating that surface-C may not have been bioavailable. At the pore scale (pore water samples), drought and the C, N amendments caused a net loss of identified molecules when the soils were rewet from below, whereas wetting from above caused a net increase in identified molecules, suggesting that fresh inputs stimulated the C-and N-limited microbial populations present deeper in the soil profile. Our experiment highlights the complex and interactive role of antecedent moisture conditions, wetting direction, and resource limitations in driving core-scale C fluxes.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/kaizadp/TES_spatial_access_2021 (DOI: 10.5281/zenodo.5522938).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems

The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗
DOE OSTI2021

Characterization of the structural forces governing the reversibility of the thermal unfolding of the human acidic fibroblast growth factor

Human acidic fibroblast growth factor (hFGF1) is an all beta-sheet protein that is involved in the regulation of key cellular processes including cell proliferation and wound healing. hFGF1 is known to aggregate when subjected to thermal unfolding. In this study, we investigate the equilibrium unfolding of hFGF1 using a wide array of biophysical and biochemical techniques. Systematic analyses of the thermal and chemical denaturation data on hFGF1 variants (Q54P, K126N, R136E, K126N/R136E, Q54P/K126N, Q54P/R136E, and Q54P/K126N/R136E) indicate that nullification of charges in the heparin-binding pocket can significantly increase the stability of wtFGF1. Triple variant (Q54P/K126N/R136E) was found to be the most stable of all the hFGF1 variants studied. With the exception of triple variant, thermal unfolding of wtFGF1 and the other variants is irreversible. Thermally unfolded triple variant refolds completely to its biologically native conformation. Microsecond-level molecular dynamic simulations reveal that a network of hydrogen bonds and salt bridges linked to Q54P, K126N, and R136E mutations, are responsible for the high stability and reversibility of thermal unfolding of the triple variant. In our opinion, the findings of the study provide valuable clues for the rational design of a stable hFGF1 variant that exhibits potent wound healing properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Structure and assembly of the diiron cofactor in the heme-oxygenase–like domain of the N-nitrosourea–producing enzyme SznF

In biosynthesis of the pancreatic cancer drug streptozotocin, the tridomain nonheme-iron oxygenase SznF hydroxylates N δ and N ω ' of N ω -methyl- L -arginine before oxidatively rearranging the triply modified guanidine to the N-methyl-N-nitrosourea pharmacophore. A previously published structure visualized the monoiron cofactor in the enzyme’s C-terminal cupin domain, which promotes the final rearrangement, but exhibited disorder and minimal metal occupancy in the site of the proposed diiron cofactor in the N-hydroxylating heme-oxygenase–like (HO-like) central domain. We leveraged our recent observation that the N-oxygenating µ-peroxodiiron(III/III) intermediate can form in the HO-like domain after the apo protein self-assembles its diiron(II/II) cofactor to solve structures of SznF with both of its iron cofactors bound. These structures of a biochemically validated member of the emerging heme-oxygenase–like diiron oxidase and oxygenase (HDO) superfamily with intact diiron cofactor reveal both the large-scale conformational change required to assemble the O 2 -reactive Fe 2 (II/II) complex and the structural basis for cofactor instability—a trait shared by the other validated HDOs. During cofactor (dis)assembly, a ligand-harboring core helix dynamically (un)folds. The diiron cofactor also coordinates an unanticipated Glu ligand contributed by an auxiliary helix implicated in substrate binding by docking and molecular dynamics simulations. The additional carboxylate ligand is conserved in another N-oxygenating HDO but not in two HDOs that cleave carbon–hydrogen and carbon–carbon bonds to install olefins. Among ~9,600 sequences identified bioinformatically as members of the emerging HDO superfamily, ~25% conserve this additional carboxylate residue and are thus tentatively assigned as N-oxygenases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ 54 , and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ 54 binding loops. Based on the DnaB analysis, we hypothesize that DnaB ‘grabs’ ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Elucidating interprotein energy transfer dynamics within the antenna network from purple bacteria

In photosynthesis, absorbed light energy transfers through a network of antenna proteins with near-unity quantum efficiency to reach the reaction center, which initiates the downstream biochemical reactions. While the energy transfer dynamics within individual antenna proteins have been extensively studied over the past decades, the dynamics between the proteins are poorly understood due to the heterogeneous organization of the network. Previously reported timescales averaged over such heterogeneity, obscuring individual interprotein energy transfer steps. Here, we isolated and interrogated interprotein energy transfer by embedding two variants of the primary antenna protein from purple bacteria, light-harvesting complex 2 (LH2), together into a near-native membrane disc, known as a nanodisc. We integrated ultrafast transient absorption spectroscopy, quantum dynamics simulations, and cryogenic electron microscopy to determine interprotein energy transfer timescales. By varying the diameter of the nanodiscs, we replicated a range of distances between the proteins. The closest distance possible between neighboring LH2, which is the most common in native membranes, is 25 Å and resulted in a timescale of 5.7 ps. Larger distances of 28 to 31 Å resulted in timescales of 10 to 14 ps. Corresponding simulations showed that the fast energy transfer steps between closely spaced LH2 increase transport distances by ∼15%. Overall, our results introduce a framework for well-controlled studies of interprotein energy transfer dynamics and suggest that protein pairs serve as the primary pathway for the efficient transport of solar energy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. Overall, this protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Integrative teaching of metabolic modeling and flux analysis with interactive python modules

Abstract The modeling of rates of biochemical reactions—fluxes—in metabolic networks is widely used for both basic biological research and biotechnological applications. A number of different modeling methods have been developed to estimate and predict fluxes, including kinetic and constraint‐based (Metabolic Flux Analysis and flux balance analysis) approaches. Although different resources exist for teaching these methods individually, to‐date no resources have been developed to teach these approaches in an integrative way that equips learners with an understanding of each modeling paradigm, how they relate to one another, and the information that can be gleaned from each. We have developed a series of modeling simulations in Python to teach kinetic modeling, metabolic control analysis, 13C‐metabolic flux analysis, and flux balance analysis. These simulations are presented in a series of interactive notebooks with guided lesson plans and associated lecture notes. Learners assimilate key principles using models of simple metabolic networks by running simulations, generating and using data, and making and validating predictions about the effects of modifying model parameters. We used these simulations as the hands‐on computer laboratory component of a four‐day metabolic modeling workshop and participant survey results showed improvements in learners' self‐assessed competence and confidence in understanding and applying metabolic modeling techniques after having attended the workshop. The resources provided can be incorporated in their entirety or individually into courses and workshops on bioengineering and metabolic modeling at the undergraduate, graduate, or postgraduate level.

Kaste, Joshua A. M.↗
DOE OSTI2025

Technology Case Study: Economic, Sustainability, and Deployment Considerations for Sustainable Aviation Fuels Produced via Lignocellulosic Sugar Catalysis

This report presents a technology case study reflecting one exemplary representative pathway for the conversion of lignocellulosic sugars to sustainable aviation fuels (SAF) via aqueous phase reforming (APR) catalysis, considered within a broader integrated biorefinery framework based on biochemical processing operations. While far from the only option for converting sugars to SAF, this pathway was selected as a case study here based on its relatively high technology maturity and simplistic processing approach (avoiding complex separations or other equipment scalability challenges), coupled with the potential for high fuel yields and favorable costs/carbon intensities with opportunities for further near-term optimization. The report considers key process integration and engineering design considerations for a modeled hypothetical, nth-plant commercial biorefinery, reflecting a number of processing options and parameters envisioned to be achievable as future goals. Resultant outputs from Aspen Plus process simulations are evaluated through techno-economic and life cycle analyses (TEA and LCA), including implications for marginal cost of CO2 abatement and inclusion of currently-applicable policy incentives. Moving beyond base case configurations, a number of alternative scenarios are also evaluated for their ability to further improve economics, greenhouse gas (GHG) emissions, and marginal cost of abatement, highlighting a path to achieve deep decarbonization goals of more than 70% GHG reduction for SAF (with the potential to reach net-negative carbon intensities in some cases) under reasonable fuel production costs. The report also highlights future opportunities and gaps for further research on this technology pathway.

09 BIOMASS FUELS↗
DOE OSTI2020

Base pairing, structural and functional insights into N 4 -methylcytidine (m 4 C) and N 4 ,N 4 -dimethylcytidine (m 4 2 C) modified RNA

The N 4 -methylation of cytidine (m 4 C and m 4 2 C) in RNA plays important roles in both bacterial and eukaryotic cells. In this work, we synthesized a series of m 4 C and m 4 2 C modified RNA oligonucleotides, conducted their base pairing and bioactivity studies, and solved three new crystal structures of the RNA duplexes containing these two modifications. Our thermostability and X-ray crystallography studies, together with the molecular dynamic simulation studies, demonstrated that m 4 C retains a regular C:G base pairing pattern in RNA duplex and has a relatively small effect on its base pairing stability and specificity. By contrast, the m 4 2 C modification disrupts the C:G pair and significantly decreases the duplex stability through a conformational shift of native Watson-Crick pair to a wobble-like pattern with the formation of two hydrogen bonds. This double-methylated m 4 2 C also results in the loss of base pairing discrimination between C:G and other mismatched pairs like C:A, C:T and C:C. The biochemical investigation of these two modified residues in the reverse transcription model shows that both mono- or di-methylated cytosine bases could specify the C:T pair and induce the G to T mutation using HIV-1 RT. In the presence of other reverse transcriptases with higher fidelity like AMV-RT, the methylation could either retain the normal nucleotide incorporation or completely inhibit the DNA synthesis. These results indicate the methylation at N 4 -position of cytidine is a molecular mechanism to fine tune base pairing specificity and affect the coding efficiency and fidelity during gene replication.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Proton Traffic Jam: Effect of Nanoconfinement and Acid Concentration on Proton Hopping Mechanism

Abstract The properties of the water network in concentrated HCl acid pools in nanometer‐sized reverse nonionic micelles were probed with TeraHertz absorption, dielectric relaxation spectroscopy, and reactive force field simulations capable of describing proton hopping mechanisms. We identify that only at a critical micelle size of W 0 =9 do solvated proton complexes form in the water pool, accompanied by a change in mechanism from Grotthuss forward shuttling to one that favors local oscillatory hopping. This is due to a preference for H + and Cl − ions to adsorb to the micelle interface, together with an acid concentration effect that causes a “traffic jam” in which the short‐circuiting of the hydrogen‐bonding motif of the hydronium ion decreases the forward hopping rate throughout the water interior even as the micelle size increases. These findings have implications for atmospheric chemistry, biochemical and biophysical environments, and energy materials, as transport of protons vital to these processes can be suppressed due to confinement, aggregation, and/or concentration.

Adams, Ellen M.↗