Double-Metal-Ion/Single-Metal-Ion Mechanisms of the Cleavage Reaction of Ribozymes: First-Principles Molecular Dynamics Simulations of a Fully Hydrated Model System

On the Recognition of Natural Substrate CTP and Endogenous Inhibitor ddhCTP of SARS-CoV-2 RNA-Dependent RNA Polymerase: A Molecular Dynamics Study

The novel coronavirus SARS-CoV-2 is the causative agent of the COVID-19 outbreak that is affecting the entire planet. As the pandemic is still spreading worldwide, with multiple mutations of the virus, it is of interest and of help to employ computational methods for identifying potential inhibitors of the enzymes responsible for viral replication. Attractive antiviral nucleotide analogue RNA-dependent RNA polymerase (RdRp) chain terminator inhibitors are investigated with this purpose. This study, based on molecular dynamics (MD) simulations, addresses the important aspects of the incorporation of an endogenously synthesized nucleoside triphosphate, ddhCTP, in comparison with the natural nucleobase cytidine triphosphate (CTP) in RdRp. The ddhCTP species is the product of the viperin antiviral protein as part of the innate immune response. The absence of the ribose 3'-OH in ddhCTP could have important implications in its inhibitory mechanism of RdRp. We built an in silico model of the RNA strand embedded in RdRp using experimental methods, starting from the cryo-electron microscopy structure and exploiting the information obtained by spectrometry on the RNA sequence. We determined that the model was stable during the MD simulation time. The obtained results provide deeper insights into the incorporation of nucleoside triphosphates, whose molecular mechanism by the RdRp active site still remains elusive.

Quantum Mechanical/Molecular Mechanical Analysis of the Catalytic Mechanism of Phosphoserine Phosphatase

Phosphoserine phosphatase (PSP), a member of the haloacid dehalogenase (HAD) superfamily that comprises the vast majority of phosphotransferases, is likely a steady-state regulator of the level of d-serine in the brain. The proposed catalytic cycle of PSP consists of a two-step mechanism: formation of a phospho-enzyme intermediate by phosphate transfer to Asp11 and its subsequent hydrolysis. Our combined quantum mechanical/molecular mechanical (QM/MM) calculations of the reaction pathways favour a dissociative mechanism of nucleophilic substitution via a trigonal-planar metaphosphate-like configuration for both steps, associated with proton transfer to the leaving group or from the nucleophile. This proton transfer is facilitated by active site residue Asp13 that acts as both a general base and a general acid. Free energy calculation on the reaction pathways further support the structural role of the enzymatic environment and the active site architecture. The choice of a proper reaction coordinate along which to bias the free energy calculations can be guided by a projection of the canonical reaction coordinate obtained from a chain-of-state optimisation onto important internal coordinates.

Development of Site-Specific Mg(2+)-RNA Force Field Parameters: A Dream or Reality? Guidelines from Combined Molecular Dynamics and Quantum Mechanics Simulations

The vital contribution of Mg²⁺ ions to RNA biology is challenging to dissect at the experimental level. This calls for the integrative support of atomistic simulations, which at the classical level are plagued by limited accuracy. Indeed, force fields intrinsically neglect nontrivial electronic effects that Mg²⁺ exerts on its surrounding ligands in varying RNA coordination environments. Here, we present a combined computational study based on classical molecular dynamics (MD) and Density Functional Theory (DFT) calculations, aimed at characterizing (i) the performance of five Mg²⁺ force field (FF) models in RNA systems and (ii) how charge transfer and polarization affect the binding of Mg²⁺ ions in different coordination motifs. As a result, a total of ∼2.5 μs MD simulations (100/200 ns for each run) for two prototypical Mg²⁺-dependent ribozymes showed remarkable differences in terms of populations of inner-sphere coordination site types. Most importantly, complementary DFT calculations unveiled that differences in charge transfer and polarization among recurrent Mg²⁺-RNA coordination motifs are surprisingly small. In particular, the charge of the Mg²⁺ ions substantially remains constant through different coordination sites, suggesting that the common philosophy of developing site-specific Mg²⁺ ion parameters is not in line with the physical origin of the Mg²⁺-RNA MD simulations inaccuracies. Overall, this study constitutes a guideline for an adept use of current Mg²⁺ models and provides novel insights for the rational development of next-generation Mg²⁺ FFs to be employed for atomistic simulations of RNA.

Catalytic metal ions and enzymatic processing of DNA and RNA

CONSPECTUS: Two-metal-ion-dependent nucleases cleave the phosphodiester bonds of nucleic acids via the two-metal-ion (2M) mechanism. Several high-resolution X-ray structures portraying the two-metal-aided catalytic site, together with mutagenesis and kinetics studies, have demonstrated a functional role of the ions for catalysis in numerous metallonucleases. Overall, the experimental data confirm the general mechanistic hypothesis for 2M-aided phosphoryl transfer originally reported by Steitz and Steitz ( Proc. Natl. Acad. Sci. U.S.A. 1993 , 90 ( 14 ), 6498 - 6502 ). This seminal paper proposed that one metal ion favors the formation of the nucleophile, while the nearby second metal ion facilitates leaving group departure during RNA hydrolysis. Both metals were suggested to stabilize the enzymatic transition state. Nevertheless, static X-ray structures alone cannot exhaustively unravel how the two ions execute their functional role along the enzymatic reaction during processing of DNA or RNA strands when moving from reactants to products, passing through metastable intermediates and high-energy transition states. In this Account, we discuss the role of multiscale molecular simulations in further disclosing mechanistic insights of 2M-aided catalysis for two prototypical enzymatic targets for drug discovery, namely, ribonuclease H (RNase H) and type II topoisomerase (topoII). In both examples, first-principles molecular simulations, integrated with structural data, emphasize a cooperative motion of the bimetal motif during catalysis. The coordinated motion of both ions is crucial for maintaining a flexible metal-centered structural architecture exquisitely tailored to accommodate the DNA or RNA sugar-phosphate backbone during phosphodiester bond cleavage. Furthermore, our analysis of RNase H and the N-terminal domain (PAN) of influenza polymerase shows that classical molecular dynamics simulations coupled with enhanced sampling techniques have contributed to describe the modulatory effect of metal ion concentration and metal uptake on the 2M mechanism and efficiency. These aspects all point to the emerging and intriguing role of additional adjacent ions potentially involved in the modulation of phosphoryl transfer reactions and enzymatic turnover in 2M-catalysis, as recently observed experimentally in polymerase η and homing endonuclease I-DmoI. These computational results, integrated with experimental findings, describe and reinforce the nascent concept of a functional and cooperative dynamics of the catalytic metal ions during the 2M-dependent enzymatic processing of DNA and RNA. Encouraged by the insights provided by computational approaches, we foresee further experiments that will feature the functional and joint dynamics of the catalytic metal ions for nucleic acid processing. This could impact the de novo design of artificial metallonucleases and the rational design of potent metal-chelating inhibitors of pharmaceutically relevant enzymes.

QM/MM molecular dynamics studies of metal binding proteins

Mixed quantum-classical (quantum mechanical/molecular mechanical (QM/MM)) simulations have strongly contributed to providing insights into the understanding of several structural and mechanistic aspects of biological molecules. They played a particularly important role in metal binding proteins, where the electronic effects of transition metals have to be explicitly taken into account for the correct representation of the underlying biochemical process. In this review, after a brief description of the basic concepts of the QM/MM method, we provide an overview of its capabilities using selected examples taken from our work. Specifically, we will focus on heme peroxidases, metallo-β-lactamases, α-synuclein and ligase ribozymes to show how this approach is capable of describing the catalytic and/or structural role played by transition (Fe, Zn or Cu) and main group (Mg) metals. Applications will reveal how metal ions influence the formation and reduction of high redox intermediates in catalytic cycles and enhance drug metabolism, amyloidogenic aggregate formation and nucleic acid synthesis. In turn, it will become manifest that the protein frame directs and modulates the properties and reactivity of the metal ions.

Aqueous solvation of methane from first principles

Structural, dynamical, bonding, and electronic properties of water molecules around a soluted methane molecule are studied from first principles. The results are compatible with experiments and qualitatively support the conclusions of recent classical molecular dynamics simulations concerning the controversial issue on the presence of "immobilized" water molecules around hydrophobic groups: the hydrophobic solute slightly reduces (by a less than 2 factor) the mobility of many surrounding water molecules rather than immobilizing just the few ones which are closest to methane, similarly to what was obtained by previous first-principles simulations of soluted methanol. Moreover, the rotational slowing down is compatible with the one predicted on the basis of the excluded volume fraction, which leads to a slower hydrogen bond exchange rate. The analysis of simulations performed at different temperatures suggests that the target temperature of the soluted system must be carefully chosen, in order to avoid artificial slowing-down effects. By generating maximally localized Wannier functions, a detailed description of the polarization effects in both solute and solvent molecules is obtained, which better characterizes the solvation process.

The structural role of Mg2+ ions in a class I RNA polymerase ribozyme: a molecular simulation study

According to the RNA world hypothesis, self-replicating ribozymes, storing the genetic information and being able to perform catalysis, were the constituents of the first living organisms. In particular, RNA polymerase ribozymes, similar to current proteinaceous enzymatic polymerases, may have been able to promote the synthesis of RNA strands in a primitive world. Polymerase catalysis is usually assisted by Mg(2+) ions, but it is not always trivial to find out experimentally the number of Mg(2+) ions placed in the active site as well as the identity and the number of their coordination ligands. Here, we addressed this issue in an artificial class I ligase ribozyme. On the basis of a recently solved crystal structure, we constructed computational models of reactant and product states of this ribozyme, considering monometallic and bimetallic species. Our models were relaxed by force field based molecular dynamics (MD) simulations and mixed quantum-classical (QM/MM) MD. The structural and dynamical properties of these models were consistent with experimental data and were validated by a comparison with the catalytic sites of proteinaceous DNA and RNA polymerases. Consistently with enzymatic polymerases, our results suggest that class I RNA ligases most probably contain two magnesium ions in the active site and they may, therefore, catalyze the junction of two RNA strands via "a two Mg(2+) ions" mechanism.

Modeling the RNA 2'OH activation: possible roles of metal ion and nucleobase as catalysts in self-cleaving ribozymes

The RNA 2'OH activation as taking place in the first chemical step of self-cleaving ribozymes is studied theoretically by DFT and MP2 methods using a continuum solvation model (CPCM). The reaction of proton transfer is studied in the presence of two kinds of catalysts: a fully hydrated metal ion (Mg(2+)) or partially hydrated nucleobase (guanine), taken separately or together leading to three different modes of activation. The metal ion is either directly bound (inner-sphere) or indirectly bound (outer-sphere) to the 2'OH group and a hydroxide ion acts as a general or specific base; the nucleobase is taken in anionic or in neutral enol-tautomeric forms playing itself the role of general base. The presence of a close metal ion (outer-sphere) lowers the pK(a) value of the 2'OH group by several log units in both metal-ion and nuleobase catalysis. The direct metal coordination to the 2'OH group (inner-sphere) further stabilizes the developing negative charge on the nucleophile. The switching from the inner-sphere to the outer-sphere coordination appears to be driven by the energy cost for reorganizing the first coordination shell rather than by the electrostatic repulsion between the ligands. The metal-ion catalysis is more effective with a specific base in the dianionic mechanism. On the other hand, the nucleobase catalysis is more effective in the monoanionic mechanism and in the presence of a metal ion acting as a cofactor through nonspecific electrostatic interactions. The results establish a baseline to study the possible roles of metal and nucleobase catalysts and their environment in more realistic models for self-cleaving ribozymes.

Protonation of a hydroxide anion bridging two divalent magnesium cations in water probed by first-principles metadynamics simulation

The protonation of a hydroxide anion (OH(-)) located between two magnesium cations (Mg(2+)) in aqueous solution has been investigated by first-principles metadynamics simulation. We observe that the complex Mg(2+)-OH(-)-Mg(2+) is stabilized by the coparticipation of the hydroxide anion to the first hydration shells of both the Mg(2+) cations. Contrary to the cases of OH(-) in pure water, the transfer of protons in the presence of the divalent metal ions turns out to be a slow chemical event. This can be ascribed to the decreased proton affinity of the bridging OH(-). Metadynamics simulation, used to overcome the difficulty of the long time scale required by the protonation of the bridging OH(-), has shown that the system possesses a great stability on the reactant state, characterized by a bioctahedral (6,6) solvation structure around the two Mg(2+) cations. The exploration of the free energy landscape shows that this stable bioctahedral configuration converts into a lower coordinated (5,6) structure, leading to a proton transfer from a water molecule belonging to the first solvation shell of the Mg(2+) ion having the lower coordination to the bridging OH(-); the free energy barrier for the protonation reaction is 11 kcal/mol, meaning that the bridging hydroxide is a weak base. During the proton transfer, the bridging OH(-) reverts to an H(2)O molecule, and this breaks the electrostatic coupling of the two Mg(2+) ions, which depart independently with their own hydration shells, one of which is entirely formed by water molecules. The second one carries the newly created OH(-). Our results show that the flexibility in the metal coordination plays a crucial role in both the protonation process of the bridging OH(-) and the separation of the metal cations, providing useful insight into the nature of proton transfer in binuclear divalent metal ions, with several biological implications, such as, for instance, transesterification of catalytic RNA.

Hammerhead ribozymes: true metal or nucleobase catalysis? Where is the catalytic power from?

The hammerhead ribozyme was first considered as a metalloenzyme despite persistent inconsistencies between structural and functional data. In the last decade, metal ions were confirmed as catalysts in self-splicing ribozymes but displaced by nucleobases in self-cleaving ribozymes. However, a model of catalysis just relying on nucleobases as catalysts does not fully fit some recent data. Gathering and comparing data on metal ions in self-cleaving and self-splicing ribozymes, the roles of divalent metal ions and nucleobases are revisited. Hypothetical models based on cooperation between metal ions and nucleobases are proposed for the catalysis and evolution of this prototype in RNA catalysis.

The effects of water on beta-D-xylose condensation reactions

Car-Parrinello-based ab initio molecular dynamics simulations (CPMD) combined with metadynamics (MTD) simulations were used to determine the reaction energetics for the beta-D-xylose condensation reaction to form beta-1,4-linked xylobiose in a dilute acid solution. Protonation of the hydroxyl group on the xylose molecule and the subsequent breaking of the C-O bond were found to be the rate-limiting step during the xylose condensation reaction. Water and water structure was found to play a critical role in these reactions due to the proton's high affinity for water molecules. The reaction free energy and reaction barrier were determined using CPMD-MTD. We found that solvent reorganization due to proton partial desolvation must be taken into account in order to obtain the correct reaction activation energy. Our calculated reaction free energy and reaction activation energy compare well with available experimental results.

Are there immobilized water molecules around hydrophobic groups? Aqueous solvation of methanol from first principles

Structural, dynamical, bonding, and electronic properties of water molecules around a soluted methanol molecule are studied from first principles. The results are compatible with experiments and qualitatively support the conclusions of recent classical molecular dynamics simulations concerning the controversial issue on the presence of "immobilized" water molecules around hydrophobic groups: the hydrophobic solute slightly reduces the mobility of many surrounding water molecules rather than immobilizing just the few ones which are closest to the methyl group. By generating maximally localized Wannier functions, a detailed description of the polarization effects in both solute and solvent molecules is obtained, which better elucidates the solvation process.

Roles of metal ions in nucleases

The hydrolysis of phosphodiester bonds by metallonucleases is crucial to most aspects of nucleic acid processing. In recent years, studies of the classical restriction endonucleases have given way to the characterization of metallonucleases with widely divergent active site motifs. These developments fuel debates regarding the roles of metal ions in these enzymes. It is fortuitous that the current literature also includes the increased application of a variety of computational techniques to test the roles of metal ions in nucleic acid hydrolysis by these systems. This includes recent proposals and indirect evidence that these enzymes utilize metal ion movement in these reactions.

First principles molecular dynamics study of catalytic reactions of biological macromolecular systems: toward analyses with QM/MM hybrid molecular simulations

First principles molecular dynamics simulations performed on a fully solvated RNA model structure allowed us to investigate the mechanism for enzymatic cleavage reactions, in vitro, of RNA enzymes (ribozymes). The concerted action of two metal catalysts turns out to be the most efficient way to promote, on the one hand, the proton abstraction from 2(')-OH that triggers the nucleophilic attack and, on the other hand, the cleavage of the P-O(5(')) bond. In fact, the elimination of one of the two metal cations leads to an increase in the activation energy of the reaction. The simulated pathway shows that an OH(-) in the coordination shell of the Mg(2+) close to O(2(')) promotes the initial proton abstraction and prevents its transfer to the ribozyme. This suggests that, in a real ribozyme, the double-metal-ion reaction mechanism is preferred with respect to single-metal-ion mechanisms either in the presence or in absence of the OH(-) anion. Finally, an insight into the importance of hybrid quantum mechanics/molecular mechanics (QM/MM) schemes is discussed in view of the modelling of a realistic system carrying all the features of a true ribozyme.

Hydration properties of magnesium and calcium ions from constrained first principles molecular dynamics

We studied the solvation structures of the divalent metal cations Mg(2+) and Ca(2+) in ambient water by applying a Car-Parrinello-based constrained molecular dynamics method. By employing the metal-water oxygen coordination number as a reaction coordinate, we could identify distinct aqua complexes characterized by structural variations of the first coordination shell. In particular, our estimated free-energy profile clearly shows that the global minimum for Mg(2+) is represented by a rather stable sixfold coordination in the octahedral arrangement, in agreement with experiments. Conversely, for Ca(2+) the free-energy curve shows several shallow local minima, suggesting that the hydration structure of Ca(2+) is highly variable. Implications for water exchange reactions are also discussed.

Possible Mechanism of Proton Transfer through Peptide Groups in the H-Pathway of the Bovine CytochromecOxidase

The peptide group connecting Tyr440 and Ser441 of the bovine cytochrome c oxidase is involved in a recently proposed proton-transfer path (H-path) where, at variance with other pathways (D- and K-paths), a usual hydrogen-bond network is interrupted, thus making this proton propagation rather unconventional. Our density-functional based molecular dynamics simulations show that, despite this anomaly and provided that a proton can reach a nearby water, a multistep proton-transfer pathway can become a viable pathway for such a reaction: a proton is initially transferred to the carbonyl oxygen of a keto form of the Tyr440-Ser441 peptide group [-CO-NH-], producing an imidic acid [-C(OH)-NH-] as a metastable state; the amide proton of the imidic acid is then transferred, spontaneously to the deprotonated carboxyl group of the Asp51 side chain, leading to the formation of an enol form [-C(OH)=N-] of the Tyr440-Ser441 peptide group. Then a subsequent enol-to-keto tautomerization occurs via a double proton-transfer path realized in the two adjacent Tyr440-Ser441 and Ser441-Asp442 peptide groups. An analysis of this multistep proton-transfer pathway shows that each elementary process occurs through the shortest distance, no permanent conformational changes are induced, thus preserving the X-ray crystal structure, and the reaction path is characterized by a reasonable activation barrier.

Two-metal-ion mechanism for hammerhead-ribozyme catalysis

The hammerhead ribozyme is one of the best studied ribozymes, but it still presents challenges for our understanding of RNA catalysis. It catalyzes a transesterification reaction that converts a 5',3' diester to a 2',3' cyclic phosphate diester via an S(N)2 mechanism. Thus, the overall reaction corresponds to that catalyzed by bovine pancreatic ribonuclease. However, an essential distinguishing aspect is that metal ions are not involved in RNase catalysis but appear to be important in ribozymes. Although various techniques have been used to assign specific functions to metals in the hammerhead ribozyme, their number and roles in catalysis is not clear. Two recent theoretical studies on RNA catalysis examined the reaction mechanism of a single-metal-ion model. A two-metal-ion model, which is supported by experiment and based on ab initio and density functional theory calculations, is described here. The proposed mechanism of the reaction has four chemical steps with three intermediates and four transition states along the reaction pathway. Reaction profiles are calculated in the gas phase and in solution. The early steps of the reaction are found to be fast (with low activation barriers), and the last step, corresponding to the departure of the leaving group, is rate limiting. This two-metal-ion model differs from the models proposed previously in that the two metal ions function not only as Lewis acids but also as general acids/bases. Comparison with experiment shows good agreement with thermodynamic and kinetic data. A detailed analysis based on natural bond orbitals (NBOs) and natural energy decomposition (NEDA) provides insights into the role of metal ions and other factors important for catalysis.

Hydration of alkali ions from first principles molecular dynamics revisited

Structural and dynamical properties of the hydration of Li(+), Na(+), and K(+) in liquid water at ambient conditions were studied by first principles molecular dynamics. Our simulations successfully captured the different hydration behavior shown by the three alkali ions as observed in experiments. The present analyses of the dependence of the self-diffusion coefficient and rotational correlation time of water on the ion concentration suggest that Li(+) (K(+)) is certainly categorized as a structure maker (breaker), whereas Na(+) acts as a weak structure breaker. An analysis of the relevant electronic structures, based on maximally localized Wannier functions, revealed that the dipole moment of H(2)O molecules in the first solvation shell of Na(+) and K(+) decreases by about 0.1 D compared to that in the bulk, due to a contraction of the oxygen lone pair orbital pointing toward the metal ion.

Investigating biological systems using first principles Car-Parrinello molecular dynamics simulations

Density functional theory (DFT)-based Car-Parrinello molecular dynamics (CPMD) simulations describe the time evolution of molecular systems without resorting to a predefined potential energy surface. CPMD and hybrid molecular mechanics/CPMD schemes have recently enabled the calculation of redox properties of electron transfer proteins in their complex biological environment. They provided structural and spectroscopic information on novel platinum-based anticancer drugs that target DNA, also setting the basis for the construction of force fields for the metal lesion. Molecular mechanics/CPMD also lead to mechanistic hypotheses for a variety of metalloenzymes. Recent advances that increase the accuracy of DFT and the efficiency of investigating rare events are further expanding the domain of CPMD applications to biomolecules.

Hsc70 ATPase: an insight into water dissociation and joint catalytic role of K+ and Mg2+ metal cations in the hydrolysis reaction

Hybrid quantum mechanics/molecular mechanics simulations, coupled to the recently introduced metadynamics method, performed on the adenosine triphosphate (ATP) of the bovine Hsc70 ATPase protein, show which specific water molecule of the solvation shell of the Mg2+ metal cation acts as a trigger in the initial phase of the ATP hydrolysis reaction in ATP synthase. Furthermore, we provide a detailed picture of the reaction mechanism, not accessible to experimental probes, that allows us to address two important issues not yet unraveled: (i) the pathway followed by a proton and a hydroxyl anion, produced upon dissociation of a putative catalytic H2O molecule, that is crucial in the selection of the reaction channel leading to the hydrolysis; (ii) the unique and cooperative role of K+ and Mg2+ metal ions in the reaction, acting as co-catalysts and promoting the release of the inorganic phosphate via an exchange of the OH- hydroxyl anion between their respective solvation shells. This is deeply different from the proton wire mechanism evidenced, for instance, in actin and lowers significantly the free energy barrier of the reaction.

Change in a protein's electronic structure induced by an explicit solvent: Anab initio fragment molecular orbital study of ubiquitin

The effect of solvation on the electronic structure of the ubiquitin protein was analyzed using the ab initio fragment molecular orbital (FMO) method. FMO calculations were performed for the protein in vacuo, and the protein was immersed in an explicit solvent shell as thick as 12 A at the HF or MP2 level by using the 6-31G* basis set. The protein's physical properties examined were the net charge, the dipole moment, the internal energy, and the solvent interaction energy. Comparison of the computational results revealed the following changes in the protein upon solvation. First, the positively charged amino acid residues on the protein surface drew electrons from the solvent, while the negatively charged ones transfer electrons to the solvent. Second, the dipole moment of the protein was enhanced as a result of the polarization. Third, the internal energy of the protein was destabilized, but the destabilization was more than compensated for by the generation of a favorable protein-solvent interaction. Finally, the energetic changes were elicited both by the electron correlation effect of the first solvent shell and by the electrostatic effect of more distant solvent molecules. These findings were consistent with the picture of the solvated protein being a polarizable molecule dissolved in a dielectric media.

Enol-to-keto Tautomerism of Peptide Groups

Density functional based simulations, performed on polyglycine containing an enol peptide group [-C(OH)N-] which is a structural isomer of a keto form [-CONH-], show that in the enol-to-keto tautomeric reaction, the enol peptide group is less stable than the keto form, and that the enol-to-keto tautomerism is characterized by a cis/trans isomerization of the C-N peptide bond. The rate-limiting step in the cis/trans isomerization is a hydrogen migration from O to N atoms in the peptide group with a transition state consisting of a four-membered ring in the cis configuration. An analysis of the cis/trans isomerization pathway shows that the mechanisms for the cis/trans isomerization are essentially different between the enol and keto forms.