An official website of the United States government
Abstract Folding of ribozymes into well-defined tertiary structures usually requires divalent cations. How Mg2+ ions direct the folding kinetics has been a long-standing unsolved problem because experiments cannot detect the positions and dynamics of ions. To address this problem, we used molecular simulations to dissect the folding kinetics of the Azoarcus ribozyme by monitoring the path each molecule takes to reach the folded state. We quantitatively establish that Mg2+ binding to specific sites, coupled with counter-ion release of monovalent cations, stimulate the formation of secondary and tertiary structures, leading to diverse pathways that include direct rapid folding and trapping in misfolded structures. In some molecules, key tertiary structural elements form when Mg2+ ions bind to specific RNA sites at the earliest stages of the folding, leading to specific collapse and rapid folding. In others, the formation of non-native base pairs, whose rearrangement is needed to reach the folded state, is the rate-limiting step. Escape from energetic traps, driven by thermal fluctuations, occurs readily. In contrast, the transition to the native state from long-lived topologically trapped native-like metastable states is extremely slow. Specific collapse and formation of energetically or topologically frustrated states occur early in the assembly process.
more »
« less
Theory and simulations for RNA folding in mixtures of monovalent and divalent cations
RNA molecules cannot fold in the absence of counterions. Experiments are typically performed in the presence of monovalent and divalent cations. How to treat the impact of a solution containing a mixture of both ion types on RNA folding has remained a challenging problem for decades. By exploiting the large concentration difference between divalent and monovalent ions used in experiments, we develop a theory based on the reference interaction site model (RISM), which allows us to treat divalent cations explicitly while keeping the implicit screening effect due to monovalent ions. Our theory captures both the inner shell and outer shell coordination of divalent cations to phosphate groups, which we demonstrate is crucial for an accurate calculation of RNA folding thermodynamics. The RISM theory for ion–phosphate interactions when combined with simulations based on a transferable coarse-grained model allows us to predict accurately the folding of several RNA molecules in a mixture containing monovalent and divalent ions. The calculated folding free energies and ion-preferential coefficients for RNA molecules (pseudoknots, a fragment of the rRNA, and the aptamer domain of the adenine riboswitch) are in excellent agreement with experiments over a wide range of monovalent and divalent ion concentrations. Because the theory is general, it can be readily used to investigate ion and sequence effects on DNA properties.
more »
« less
- Award ID(s):
- 1900093
- PAR ID:
- 10188687
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 116
- Issue:
- 42
- ISSN:
- 0027-8424
- Page Range / eLocation ID:
- 21022 to 21030
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Rational design of molecular chelating agents requires a detailed understanding of physicochemical ligand–metal interactions in solvent phase. Computational quantum chemistry methods should be able to provide this, but computational reports have shown poor accuracy when determining absolute binding constants for many chelating molecules. To understand why, we compare and benchmark static- and dynamics-based computational procedures for a range of monovalent and divalent cations binding to a conventional cryptand molecule: 2.2.2-cryptand ([2.2.2]). The benchmarking comparison shows that dynamics simulations using standard OPLS-AA classical potentials can reasonably predict binding constants for monovalent cations, but these procedures fail for divalent cations. We also consider computationally efficient static procedure using Kohn–Sham density functional theory (DFT) and cluster-continuum modeling that accounts for local microsolvation and pH effects. This approach accurately predicts binding energies for monovalent and divalent cations with an average error of 3.2 kcal mol −1 compared to experiment. This static procedure thus should be useful for future molecular screening efforts, and high absolute errors in the literature may be due to inadequate modeling of local solvent and pH effects.more » « less
-
Electric double layers are crucial to energy storage and electrocatalytic device performance. While double layer formation originates in electrostatic interactions, electric double layer properties are governed by a balance of both electrostatic and entropic driving forces; favorable ion-surface electrostatic interactions attract counterions to charged surfaces to compensate, or "screen," potentials, but the confinement of these same ions from a bulk reservoir to the interface incurs an entropic penalty. Here, we use a dicationic imidazolium ionic liquid and its monovalent analogue to explore how cation valence and entropy influence double layer formation and electrochemical reactivity using CO2 electroreduction as a model reaction. We find that divalent and monovalent cations display similar CO2 reduction kinetics but differ vastly in steady-state reactivity due to rapid electrochemically induced precipitation of insulating dicationic (bi)carbonate films. Using in situ surface-enhanced Raman scattering spectroscopy, we find that potential-dependent reorientation occurs at similar potentials between the two ionic liquids, but the introduction of a covalent link in the divalent cation imparts a more ordered double layer structure that favors (bi)carbonate precipitation. In mixed monovalent-divalent electrolytes, we find that the divalent cations dominate interfacial properties by preferentially accumulating at surfaces even at very low relative concentrations. Our findings confirm that ion entropy plays a key role in modulating local electrochemical environments and highlight how double layer properties are very sensitive to the properties of counterions that pay the lowest entropic penalty to accumulate at interfaces. Overall, we illustrate that ion entropy provides a new knob to tune reaction microenvironments and unveil how entropy plays a major role in modulating electrochemical reactivity in mixed ion electrolytes.more » « less
-
Understanding the mechanisms of molecular transport in polyamide membranes is imperative to improve their solute-specific selectivity. We explored the partitioning behaviors of water and salts in polyamide membranes to elucidate the role of ion-membrane interactions in the transport. Quartz crystal microbalance (QCM) was employed to quantify the mass uptake at different temperatures and determine partition energies (Ek) for water and salts under two different pH values. Zeta potential and permeability tests were conducted to support the ionmembrane affinity trends observed with QCM and link these trends to ion-ion selectivity. Our results demonstrate a high affinity of water to the polyamide membrane (Ek < 0), with a significant swelling effect attributed to dipole interactions and hydrogen bonding. Ion partitioning revealed distinct differences between monovalent and divalent cations, as well as between kosmotropic and chaotropic anions. Specifically, divalent cations (Ca2+ and Mg2+) exhibited considerably lower partition energies (-0.99 and 0.29 kcal mol-1, respectively) and more efficient charge neutralization, indicating stronger interactions with the membrane compared to monovalent cations (~2.2 kcal mol-1). The partition energies of the chaotropic iodide and kosmotropic sulphate anions were substantially different (-5.5 and 4.0 kcal mol-1, respectively), likely due to the different tendency of these anions to shed their hydration shell and stick to the polymer. Last, our permeability tests indicate the potential existence of an intrinsic tradeoff between ion partitioning and intrapore diffusion, presumably due to the opposite effects that ion-membrane interactions have on these transport steps. Overall, our work underscores the role of ionspecific interactions in membrane transport and selectivity.more » « less
-
Abstract Detailed crystallographic characterization of a tri-aspartate metal-binding site previously identified on the three-fold symmetry axis of a hexameric enzyme, LarE from Lactobacillus plantarum, was conducted. By screening an array of monovalent, divalent, and trivalent metal ions, we demonstrated that this metal binding site stoichiometrically binds Ca2+, Mn2+, Fe2+/Fe3+, Co2+, Ni2+, Cu2+, Zn2+, and Cd2+, but not monovalent metal ions, Cr3+, Mg2+, Y3+, Sr2+ or Ba2+. Extensive database searches resulted in only 13 similar metal binding sites in other proteins, indicative of the rareness of tri-aspartate architectures, which allows for engineering such a selective multivalent metal ion binding site into target macromolecules for structural and biophysical characterization.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1073/pnas.1911632116
