Role of magnesium ions in the reaction mechanism at the interface between Tm1631 protein and its DNA ligand

Competitive Influence of Alkali Metals in the Ion Atmosphere on Nucleic Acid Duplex Stability

The nonspecific atmosphere around nucleic acids, often termed the ion atmosphere, encompasses a collection of weak ion-nucleic acid interactions. Although nonspecific, the ion atmosphere has been shown to influence nucleic acid folding and structural stability. Studies investigating the composition of the ion atmosphere have shown competitive occupancy of the atmosphere between metal ions in the same solution. Many studies have investigated single ion effects on nucleic acid secondary structure stability; however, no comprehensive studies have investigated how the competitive occupancy of mixed ions in the ion atmosphere influences nucleic acid secondary structure stability. Here, six oligonucleotides were optically melted in buffers containing molar quantities, or mixtures, of either XCl (X = Li, K, Rb, or Cs) or NaCl. A correction factor was developed to better predict RNA duplex stability in solutions containing mixed XCl/NaCl. For solutions containing a 1:1 mixture of XCl/NaCl, one alkali metal chloride contributed more to duplex stability than the other. Overall, there was a 54% improvement in predictive capabilities with the correction factor compared with the standard 1.0 M NaCl nearest-neighbor models. This correction factor can be used in models to better predict RNA secondary structure in solutions containing mixed XCl/NaCl.

Ligand Binding Mechanism and Its Relationship with Conformational Changes in Adenine Riboswitch

Riboswitches are naturally occurring RNA aptamers that control the expression of essential bacterial genes by binding to specific small molecules. The binding with both high affinity and specificity induces conformational changes. Thus, riboswitches were proposed as a possible molecular target for developing antibiotics and chemical tools. The adenine riboswitch can bind not only to purine analogues but also to pyrimidine analogues. Here, long molecular dynamics (MD) simulations and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) computational methodologies were carried out to show the differences in the binding model and the conformational changes upon five ligands binding. The binding free energies of the guanine riboswitch aptamer with C74U mutation complexes were compared to the binding free energies of the adenine riboswitch (AR) aptamer complexes. The calculated results are in agreement with the experimental data. The differences for the same ligand binding to two different aptamers are related to the electrostatic contribution. Binding dynamical analysis suggests a flexible binding pocket for the pyrimidine ligand in comparison with the purine ligand. The 18 μs of MD simulations in total indicate that both ligand-unbound and ligand-bound aptamers transfer their conformation between open and closed states. The ligand binding obviously affects the conformational change. The conformational states of the aptamer are associated with the distance between the mass center of two key nucleotides (U51 and A52) and the mass center of the other two key nucleotides (C74 and C75). The results suggest that the dynamical character of the binding pocket would affect its biofunction. To design new ligands of the adenine riboswitch, it is recommended to consider the binding affinities of the ligand and the conformational change of the ligand binding pocket.

Cations in motion: QM/MM studies of the dynamic and electrostatic roles of H+ and Mg2+ ions in enzyme reactions

Here we discuss current trends in the simulations of enzymatic reactions focusing on phosphate catalysis. The mechanistic details of the proton transfers coupled to the phosphate cleavage is one of the key challenges in QM/MM calculations of these and other enzyme catalyzed reactions. The lack of experimental information offers both an opportunity for computations as well as often unresolved controversies. We discuss the example of small GTPases including the important human Ras protein. The high dimensionality and chemical complexity of these reactions demand carefully chosen computational techniques both in terms of the underlying quantum chemical theory and the sampling of the conformational ensemble. We also point out the important role of Mg²⁺ ions, and recent advances in their transient involvement in the catalytic mechanisms.

Ligand Selectivity Mechanism and Conformational Changes in Guanine Riboswitch by Molecular Dynamics Simulations and Free Energy Calculations

Riboswitches regulate gene expression through direct and specific interactions with small metabolite molecules. Binding of a ligand to its RNA target is high selectivity and affinity and induces conformational changes of the RNA's secondary and tertiary structure. The structural difference of two purine riboswitches aptamers is caused by only one single mutation, where cytosine 74 in the guanine riboswitch is corresponding to a uracil 74 in adenine riboswitch. Here we employed molecular dynamics (MD) simulation, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and thermodynamic integration computational methodologies to evaluate the energetic and conformational changes of ligands binding to purine riboswitches. The snapshots used in MM-PBSA calculation were extracted from ten 50 ns MD simulation trajectories for each complex. These free energy results are in consistent with the experimental data and rationalize the selectivity of the riboswitches for different ligands. In particular, it is found that the loss in binding free energy upon mutation is mainly electrostatic in guanine (GUA) and riboswitch complex. Furthermore, new hydrogen bonds are found in mutated complexes. To reveal the conformational properties of guanine riboswitch, we performed a total of 6 μs MD simulations in both the presence and the absence of the ligand GUA. The MD simulations suggest that the conformation of guanine riboswitch depends on the distance of two groups in the binding pocket of ligand. The conformation is in a close conformation when U51-A52 is close to C74-U75.