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

Third Metal Ion Dictates the Catalytic Activity of the Two-Metal-Ion Pre-Ribosomal RNA-Processing Machinery

Nucleic acid processing enzymes use a two-Mg2+-ion motif to promote the formation and cleavage of phosphodiester bonds. Yet, recent evidence demonstrates the presence of spatially conserved second-shell cations surrounding the catalytic architecture of proteinaceous and RNA-dependent enzymes. The RNase mitochondrial RNA processing (MRP) complex, which cleaves the ribosomal RNA (rRNA) precursor at the A3 cleavage site to yield mature 5'-end of 5.8S rRNA, hosts in the catalytic core one atypically-located Mg2+ ion, in addition to the ions forming the canonical catalytic motif. Here, we employ biased quantum classical molecular dynamics simulations of RNase MRP to discover that the third Mg2+ ion inhibits the catalytic process. Instead, its displacement in favour of a second-shell monovalent K+ ion propels phosphodiester bond cleavage by enabling the formation of a specific hydrogen bonding network that mediates the essential proton transfer step. This study points to a direct involvement of a transient K+ ion in the catalytic cleavage of the phosphodiester bond and implicates cation trafficking as a general mechanism in nucleic acid processing enzymes and ribozymes.

Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation

Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.

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.

Unveiling the "Template-Dependent" Inhibition on the Viral Transcription of SARS-CoV-2

Remdesivir is one nucleotide analogue prodrug capable to terminate RNA synthesis in SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) by two distinct mechanisms. Although the "delayed chain termination" mechanism has been extensively investigated, the "template-dependent" inhibitory mechanism remains elusive. In this study, we have demonstrated that remdesivir embedded in the template strand seldom directly disrupted the complementary NTP incorporation at the active site. Instead, the translocation of remdesivir from the +2 to the +1 site was hindered due to the steric clash with V557. Moreover, we have elucidated the molecular mechanism characterizing the drug resistance upon V557L mutation. Overall, our studies have provided valuable insight into the "template-dependent" inhibitory mechanism exerted by remdesivir on SARS-CoV-2 RdRp and paved venues for an alternative antiviral strategy for the COVID-19 pandemic. As the "template-dependent" inhibition occurs across diverse viral RdRps, our findings may also shed light on a common acting mechanism of inhibitors.

Unraveling the binding mechanism of the active form of Remdesivir to RdRp of SARS-CoV-2 and designing new potential analogues: Insights from molecular dynamics simulations

The binding of the active form of Remdesivir (RTP) to RNA-dependent RNA Polymerase (RdRp) of SARS-CoV-2 was studied using molecular dynamics simulation. The RTP maintained the interactions observed in the experimental cryo-EM structure. Next, we designed new analogues of RTP, which not only binds to the RNA primer strand in a similar pose as that of RTP, but also binds more strongly than RTP does as predicted by MM-PBSA binding energy. This suggest that these analogues might be able to covalently link to the primer strand as RTP, but their 3' modification would terminate the primer strand growth.

The mechanism of action of T-705 as a unique delayed chain terminator on influenza viral polymerase transcription

Favipiravir (T-705) has been developed as a potent anti-influenza drug and exhibited a strong inhibition effect against a broad spectrum of RNA viruses. Its active form, ribofuranosyl-triphosphate (T-705-RTP), functions as a competitive substrate for the RNA-dependent RNA polymerase (RdRp) of the influenza A virus (IAV). However, the exact inhibitory mechanisms of T-705 remain elusive and subject to a long-standing debate. Although T-705 has been proposed to inhibit transcription by acting as a chain terminator, it is also paradoxically suggested to be a mutagen towards IAV RdRp by inducing mutations due to its ambiguous base pairing of C and U. Here, we combined biochemical assay with molecular dynamics (MD) simulations to elucidate the molecular mechanism underlying the inhibitory functions exerted by T-705 in IAV RdRp. Our in vitro transcription assay illustrated that IAV RdRp could recognize T-705 as a purine analogue and incorporate it into the nascent RNA strand. Incorporating a single T-705 is incapable of inhibiting transcription as extra natural nucleotides can be progressively added. However, when two consecutive T-705 are incorporated, viral transcription is completely terminated. MD simulations reveal that the sequential appearance of two T-705 in the nascent strand destabilizes the active site and disrupts the base stacking of the nascent RNA. Altogether, our results provide a plausible explanation for the inhibitory roles of T-705 targeting IAV RdRp by integrating the computational and experimental methods. Our study also offers a comprehensive platform to investigate the inhibition effect of antivirals and a novel explanation for the designing of anti-flu drugs.

Incorporation efficiency and inhibition mechanism of 2'-substituted nucleotide analogs against SARS-CoV-2 RNA-dependent RNA polymerase

The ongoing pandemic caused by SARS-CoV-2 emphasizes the need for effective therapeutics. Inhibition of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) by nucleotide analogs provides a promising antiviral strategy. One common group of RdRp inhibitors, 2'-modified nucleotides, are reported to exhibit different behaviors in the SARS-CoV-2 RdRp transcription assay. Three of these analogs, 2'-O-methyl UTP, Sofosbuvir, and 2'-methyl CTP, act as effective inhibitors in previous biochemical experiments, while Gemcitabine and ara-UTP show no inhibitory activity. To understand the impact of the 2'-modification on their inhibitory effects, we conducted extensive molecular dynamics simulations and relative binding free energy calculations using the free energy perturbation method on SARS-CoV-2 replication-transcription complex (RTC) with these five nucleotide analogs. Our results reveal that the five nucleotide analogs display comparable binding affinities to SARS-CoV-2 RdRp and they can all be added to the nascent RNA chain. Moreover, we examine how the incorporation of these nucleotide triphosphate (NTP) analogs will impact the addition of the next nucleotide. Our results indicate that 2'-O-methyl UTP can weaken the binding of the subsequent NTP and consequently lead to partial chain termination. Additionally, Sofosbuvir and 2'-methyl CTP can cause immediate termination due to the strong steric hindrance introduced by their bulky 2'-methyl groups. In contrast, nucleotide analogs with smaller substitutions, such as the fluorine atoms and the ara-hydroxyl group in Gemcitabine and ara-UTP, have a marginal impact on the polymerization process. Our findings are consistent with experimental observations, and more importantly, shed light on the detailed molecular mechanism of SARS-CoV-2 RdRp inhibition by 2'-substituted nucleotide analogs, and may facilitate the rational design of antiviral agents to inhibit SARS-CoV-2 RdRp.

Computing Metal-Binding Proteins for Therapeutic Benefit

Over one third of biomolecules rely on metal ions to exert their cellular functions. Metal ions can play a structural role by stabilizing the structure of biomolecules, a functional role by promoting a wide variety of biochemical reactions, and a regulatory role by acting as messengers upon binding to proteins regulating cellular metal-homeostasis. These diverse roles in biology ascribe critical implications to metal-binding proteins in the onset of many diseases. Hence, it is of utmost importance to exhaustively unlock the different mechanistic facets of metal-binding proteins and to harness this knowledge to rationally devise novel therapeutic strategies to prevent or cure pathological states associated with metal-dependent cellular dysfunctions. In this compendium, we illustrate how the use of a computational arsenal based on docking, classical, and quantum-classical molecular dynamics simulations can contribute to extricate the minutiae of the catalytic, transport, and inhibition mechanisms of metal-binding proteins at the atomic level. This knowledge represents a fertile ground and an essential prerequisite for selectively targeting metal-binding proteins with small-molecule inhibitors aiming to (i) abrogate deregulated metal-dependent (mis)functions or (ii) leverage metal-dyshomeostasis to selectively trigger harmful cells death.

1'-Ribose cyano substitution allows Remdesivir to effectively inhibit nucleotide addition and proofreading during SARS-CoV-2 viral RNA replication

COVID-19 has recently caused a global health crisis and an effective interventional therapy is urgently needed. Remdesivir is one effective inhibitor for SARS-CoV-2 viral RNA replication. It supersedes other NTP analogues because it not only terminates the polymerization activity of RNA-dependent RNA polymerase (RdRp), but also inhibits the proofreading activity of intrinsic exoribonuclease (ExoN). Even though the static structure of Remdesivir binding to RdRp has been solved and biochemical experiments have suggested it to be a "delayed chain terminator", the underlying molecular mechanisms is not fully understood. Here, we performed all-atom molecular dynamics (MD) simulations with an accumulated simulation time of 24 microseconds to elucidate the inhibitory mechanism of Remdesivir on nucleotide addition and proofreading. We found that when Remdesivir locates at an upstream site in RdRp, the 1'-cyano group experiences electrostatic interactions with a salt bridge (Asp865-Lys593), which subsequently halts translocation. Our findings can supplement the current understanding of the delayed chain termination exerted by Remdesivir and provide an alternative molecular explanation about Remdesivir's inhibitory mechanism. Such inhibition also reduces the likelihood of Remdesivir to be cleaved by ExoN acting on 3'-terminal nucleotides. Furthermore, our study also suggests that Remdesivir's 1'-cyano group can disrupt the cleavage site of ExoN via steric interactions, leading to a further reduction in the cleavage efficiency. Our work provides plausible and novel mechanisms at the molecular level of how Remdesivir inhibits viral RNA replication, and our findings may guide rational design for new treatments of COVID-19 targeting viral replication.

Modeling the Binding Mechanism of Remdesivir, Favilavir, and Ribavirin to SARS-CoV-2 RNA-Dependent RNA Polymerase

Recent efforts to repurpose drugs to combat COVID-19 have identified Remdesivir as a candidate. It acts on the RNA-dependent, RNA polymerase (RdRp) of the SARS-CoV-2 virus, a protein complex responsible for mediating replication of the virus's genome. However, its exact action mechanism, and that of other nucleotide analogue inhibitors, is not known. In this study, we examine at the molecular level the interaction of this drug and that of similar nucleotide analogue inhibitors, ribavirin and favilavir, by relying on atomistic molecular simulations and advanced sampling. By analyzing the binding free energies of these different drugs, it is found that all of them bind strongly at the active site. Surprisingly, however, ribavirin and favilavir do not bind the nucleotide on the complementary strand as effectively and seem to act by a different mechanism than remdesivir. Remdesivir exhibits similar binding interactions to the natural base adenine. Moreover, by analyzing remdesivir at downstream positions of the RNA, we also find that, consistent with a "delayed" termination mechanism, additional nucleotides can be incorporated after remdesivir is added, and its highly polar 1'-cyano group induces a set of conformational changes that can affect the normal RdRp complex function. By analyzing the fluctuations of residues that are altered by remdesivir binding, and comparing them to those induced by lethal point mutations, we find a possible secondary mechanism in which remdesivir destabilizes the protein complex and its interactions with the RNA strands.

Atomic-Level Mechanism of Pre-mRNA Splicing in Health and Disease

Intron removal from premature-mRNA (pre-mRNA splicing) is an essential part of gene expression and regulation that is required for the production of mature, protein-coding mRNA. The spliceosome (SPL), a majestic machine composed of five small nuclear RNAs and hundreds of proteins, behaves as an eminent transcriptome tailor, efficiently performing splicing as a protein-directed metallo-ribozyme. To select and excise long and diverse intronic sequences with single-nucleotide precision, the SPL undergoes a continuous compositional and conformational remodeling, forming eight distinct complexes throughout each splicing cycle. Splicing fidelity is of paramount importance to preserve the integrity of the proteome. Mutations in splicing factors can severely compromise the accuracy of this machinery, leading to aberrant splicing and altered gene expression. Decades of biochemical and genetic studies have provided insights into the SPL's composition and function, but its complexity and plasticity have prevented an in-depth mechanistic understanding. Single-particle cryogenic electron microscopy techniques have ushered in a new era for comprehending eukaryotic gene regulation, providing several near-atomic resolution structures of the SPL from yeast and humans. Nevertheless, these structures represent isolated snapshots of the splicing process and are insufficient to exhaustively assess the function of each SPL component and to unravel particular facets of the splicing mechanism in a dynamic environment.In this Account, building upon our contributions in this field, we discuss the role of biomolecular simulations in uncovering the mechanistic intricacies of eukaryotic splicing in health and disease. Specifically, we showcase previous applications to illustrate the role of atomic-level simulations in elucidating the function of specific proteins involved in the architectural reorganization of the SPL along the splicing cycle. Moreover, molecular dynamics applications have uniquely contributed to decrypting the channels of communication required for critical functional transitions of the SPL assemblies. They have also shed light on the role of carcinogenic mutations in the faithful selection of key intronic regions and the molecular mechanism of splicing modulators. Additionally, we emphasize the role of quantum-classical molecular dynamics in unraveling the chemical details of pre-mRNA cleavage in the SPL and in its evolutionary ancestors, group II intron ribozymes. We discuss methodological pitfalls of multiscale calculations currently used to dissect the splicing mechanism, presenting future challenges in this field. The results highlight how atomic-level simulations can enrich the interpretation of experimental results. We envision that the synergy between computational and experimental approaches will aid in developing innovative therapeutic strategies and revolutionary gene modulation tools to fight the over 200 human diseases associated with splicing misregulation, including cancer and neurodegeneration.

Frontiers of metal-coordinating drug design

Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.

Structural Basis of the Potential Binding Mechanism of Remdesivir to SARS-CoV-2 RNA-Dependent RNA Polymerase

Starting from late 2019, the coronavirus disease 2019 (COVID-19) has emerged as a once-in-a-century pandemic with deadly consequences, which urgently calls for new treatments, cures, and supporting apparatuses. Recently, because of its positive results in clinical trials, remdesivir was approved by the Food and Drug Administration to treat COVID-19 through Emergency Use Authorization. Here, we used molecular dynamics simulations and free energy perturbation methods to study the inhibition mechanism of remdesivir to its target SARS-CoV-2 virus RNA-dependent RNA polymerase (RdRp). We first constructed the homology model of this polymerase based on a previously available structure of SARS-CoV NSP12 RdRp (with a sequence identity of 95.8%). We then built a putative preinsertion binding structure by aligning the remdesivir + RdRp complex to the ATP bound poliovirus RdRp without the RNA template. The putative binding structure was further optimized with molecular dynamics simulations. The resulting stable preinsertion state of remdesivir appeared to form hydrogen bonds with the RNA template when aligned with the newly solved cryo-EM structure of SARS-CoV-2 RdRp. The relative binding free energy between remdesivir and ATP was calculated to be -2.80 ± 0.84 kcal/mol, where remdesivir bound much stronger to SARS-CoV-2 RdRp than the natural substrate ATP. The ∼100-fold improvement in the K_d from remdesivir over ATP indicates an effective replacement of ATP in blocking of the RdRp preinsertion site. Key residues D618, S549, and R555 are found to be the contributors to the binding affinity of remdesivir. These findings suggest that remdesivir can potentially act as a SARS-CoV-2 RNA-chain terminator, effectively stopping its RNA replication, with key residues also identified for future lead optimization and/or drug resistance studies.

Structural Basis of the Potential Binding Mechanism of Remdesivir to SARS-CoV-2 RNA-Dependent RNA Polymerase

Starting from late 2019, the coronavirus disease 2019 (COVID-19) has emerged as a once-in-a-century pandemic with deadly consequences, which urgently calls for new treatments, cures, and supporting apparatuses. Recently, because of its positive results in clinical trials, remdesivir was approved by the Food and Drug Administration to treat COVID-19 through Emergency Use Authorization. Here, we used molecular dynamics simulations and free energy perturbation methods to study the inhibition mechanism of remdesivir to its target SARS-CoV-2 virus RNA-dependent RNA polymerase (RdRp). We first constructed the homology model of this polymerase based on a previously available structure of SARS-CoV NSP12 RdRp (with a sequence identity of 95.8%). We then built a putative preinsertion binding structure by aligning the remdesivir + RdRp complex to the ATP bound poliovirus RdRp without the RNA template. The putative binding structure was further optimized with molecular dynamics simulations. The resulting stable preinsertion state of remdesivir appeared to form hydrogen bonds with the RNA template when aligned with the newly solved cryo-EM structure of SARS-CoV-2 RdRp. The relative binding free energy between remdesivir and ATP was calculated to be -2.80 ± 0.84 kcal/mol, where remdesivir bound much stronger to SARS-CoV-2 RdRp than the natural substrate ATP. The ∼100-fold improvement in the Kd from remdesivir over ATP indicates an effective replacement of ATP in blocking of the RdRp preinsertion site. Key residues D618, S549, and R555 are found to be the contributors to the binding affinity of remdesivir. These findings suggest that remdesivir can potentially act as a SARS-CoV-2 RNA-chain terminator, effectively stopping its RNA replication, with key residues also identified for future lead optimization and/or drug resistance studies.

Regioselective formation of RNA strands in the absence of magnesium ions

The oligomerization of ribonucleotides can produce short RNA strands in the absence of enzymes. This reaction gives one of two regioisomeric phosphodiester linkages, a 2',5'- or a 3',5'-diester. The former, non-natural linkage is detrimental for duplex stability, and is known to form preferentially in oligomerizations occurring in homogeneous solution with preactivated nucleotides in the presence of magnesium cations. We have studied ribonucleotide oligomerization with in situ activation, using NMR as monitoring technique. Unexpectedly, the known preference for 2',5'-linkages in the oligomerization of AMP was reversed in the absence of magnesium ions at slightly basic pH. Further, oligomerization was surprisingly efficient in the absence of Mg2+ salts, producing oligomers long enough for duplex formation. A quantitative systems chemistry analysis then revealed that the absence of magnesium ions favors the activation of nucleotides, and that the high concentration of active species can compensate for slower coupling. Further, organocatalytic intermediates can help to overcome the unfavorable regioselectivity of the magnesium-catalyzed reactions. Our findings allay concerns that RNA may have been difficult to form in the absence of enzymes. They also show that there is an efficient path to genetic material that does not require mineral surfaces or cations known to catalyze RNA hydrolysis.

Effects of Salt Stress on Growth, Photosynthesis, and Mineral Nutrients of 18 Pomegranate (Punica granatum) Cultivars

Pomegranate (Punica granatum L.) is widely grown in arid and semiarid regions, where the salinization may have developed through irrigation. A greenhouse experiment was conducted to investigate NaCl stress on growth, photosynthesis, and nutrients of 18 pomegranate cultivars. One group was irrigated twice a week with a nutrient solution. The other group was watered twice a week with the same nutrient solution and 200 mM NaCl for five weeks. Dry weight, shoot length, new shoot number, root length and number, leaf area, leaf relative water content, and net photosynthesis of salt-treated plants were negatively impacted by salt stress, and there was a significant difference among cultivars. Few foliar damages were observed. Na content of plants significantly increased in all cultivars, while P, S, K, Ca, Mg, Si, Al, Zn content of plants decreased under salt stress. Fe, Mn, and Cu content increased in most cultivars. Pomegranate accumulated supraoptimal Na mostly in roots and transported more K and Ca to shoots, which was attributed to maintaining a higher ratio of K/Na and Ca/Na in the aerial part of plants. Ten of the 18 cultivars were considered salt-tolerant, which would offer a reference for pomegranate cultivation on saline lands.

Disclosing the Impact of Carcinogenic SF3b Mutations on Pre-mRNA Recognition Via All-Atom Simulations

The spliceosome accurately promotes precursor messenger-RNA splicing by recognizing specific noncoding intronic tracts including the branch point sequence (BPS) and the 3'-splice-site (3'SS). Mutations of Hsh155 (yeast)/SF3B1 (human), which is a protein of the SF3b factor involved in BPS recognition and induces altered BPS binding and 3'SS selection, lead to mis-spliced mRNA transcripts. Although these mutations recur in hematologic malignancies, the mechanism by which they change gene expression remains unclear. In this study, multi-microsecond-long molecular-dynamics simulations of eighth distinct ∼700,000 atom models of the spliceosome Bact complex, and gene sequencing of SF3B1, disclose that these carcinogenic isoforms destabilize intron binding and/or affect the functional dynamics of Hsh155/SF3B1 only when binding non-consensus BPSs, as opposed to the non-pathogenic variants newly annotated here. This pinpoints a cross-talk between the distal Hsh155 mutation and BPS recognition sites. Our outcomes unprecedentedly contribute to elucidating the principles of pre-mRNA recognition, which provides critical insights on the mechanism underlying constitutive/alternative/aberrant splicing.

The pivotal role of MBD4–ATP7B in the human Cu(i) excretion path as revealed by EPR experiments and all-atom simulations

Atox1–MBD4 interaction mediates the in-cell Cu(i) concentration.

Conformational Dynamics of thiM Riboswitch To Understand the Gene Regulation Mechanism Using Markov State Modeling and the Residual Fluctuation Network Approach

Thiamine pyrophosphate (TPP) riboswitch is a cis-regulatory element in the noncoding region of mRNA. The aptamer domain of TPP riboswitch detects the high abundance of coenzyme thiamine pyrophosphate (TPP) and modulates the gene expression for thiamine synthetic gene. The mechanistic understanding in recognition of TPP in aptamer domain and ligand-induced compactness for folding of expression platform are most important to designing novel modulators. To understand the dynamic behavior of TPP riboswitch upon TPP binding, molecular dynamics simulations were performed for 400 ns in both apo and TPP bound forms of thiM riboswitch from E. coli and analyzed in terms of eRMSD-based Markov state modeling and residual fluctuation network. Markov state models show good correlations in transition probability among metastable states from simulated trajectory and generated models. Structural compactness in TPP bound form is observed which is correlated with SAXS experiment. The importance of junction of P4 and P5 is evident during dynamics, which correlates with FRET analysis. The dynamic nature of two sensor forearms is due to the flexible P1 helix, which is its intrinsic property. The transient state in TPP-bound form was observed in the Markov state model, along with stable states. We believe that this transient state is responsible to assist the influx and outflux of ligand molecule by creating a solvent channel around the junction region of P4 and P5 and such a structure was anticipated in FRET analysis. The dynamic nature of riboswitch is dependent on the interaction between residues on distal loops L3 and L5/P3 and junction P4 and P5, J3/2 which stabilize the J2/4. It helps in the transfer of allosteric information between J2/4 and P3/L5 tertiary docking region through the active site residues. Understanding such information flow will benefit in highlighting crucial residues in highly dynamic and kinetic systems. Here, we report the residues and segments in riboswitch that play vital roles in providing stability and this can be exploited in designing inhibitors to regulate the functioning of riboswitches.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

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.

Molecular dynamic simulations of protein/RNA complexes: CRISPR/Csy4 endoribonuclease

BACKGROUND

Many prokaryotic genomes comprise Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) offering defense against foreign nucleic acids. These immune systems are conditioned by the production of small CRISPR-derived RNAs matured from long RNA precursors. This often requires a Csy4 endoribonuclease cleaving the RNA 3'-end.

METHODS

We report extended explicit solvent molecular dynamic (MD) simulations of Csy4/RNA complex in precursor and product states, based on X-ray structures of product and inactivated precursor (55 simulations; ~3.7μs in total).

RESULTS

The simulations identify double-protonated His29 and deprotonated terminal phosphate as the likely dominant protonation states consistent with the product structure. We revealed potential substates consistent with Ser148 and His29 acting as the general base and acid, respectively. The Ser148 could be straightforwardly deprotonated through solvent and could without further structural rearrangements deprotonate the nucleophile, contrasting similar studies investigating the general base role of nucleobases in ribozymes. We could not locate geometries consistent with His29 acting as general base. However, we caution that the X-ray structures do not always capture the catalytically active geometries and then the reactive structures may be unreachable by the simulation technique.

CONCLUSIONS

We identified potential catalytic arrangement of the Csy4/RNA complex but we also report limitations of the simulation technique. Even for the dominant protonation state we could not achieve full agreement between the simulations and the structural data.

GENERAL SIGNIFICANCE

Potential catalytic arrangement of the Csy4/RNA complex is found. Further, we provide unique insights into limitations of simulations of protein/RNA complexes, namely, the influence of the starting experimental structures and force field limitations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

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.

Choline ion interactions with DNA atoms explain unique stabilization of A-T base pairs in DNA duplexes: a microscopic view

Under physiological conditions, G-C base pairs are more stable than A-T base pairs. In a previous study, we showed that in the hydrated ionic liquid of choline dihydrogen phosphate, the stabilities of these base pairs are reversed. In the present study, we elucidated the unique binding interactions of choline ions with DNA atoms from a microscopic viewpoint using molecular dynamics simulations. Three times more choline ions bind to the DNA duplex than sodium ions. Sodium ions bind closely but not stably; in contrast, the choline ions bind through multiple hydrogen bonding networks with DNA atoms stably. The affinity of choline ion for the minor groove of A-T base pairs is more than 2 times that for other groove areas. In the narrow A-T minor groove, choline ion has high affinity for the ribose atoms of thymine. Choline ions also destabilize the formation of hydrogen bonds between G-C base pairs by binding to base atoms preferentially for both of duplex and single-strand DNA, which are associated with the bonds between G-C base pairs. Our new finding will not only lead to better control of DNA stability for use in DNA nanodevices, but also provide new insight into the stability of DNA duplexes under crowding conditions found in living cells.

Molecular modeling of nucleic acid structure: energy and sampling

An overview of computer simulation techniques as applied to nucleic acid systems is presented. This unit discusses methods used to treat energy and to sample representative configurations. Emphasis is placed on molecular mechanics and empirical force fields.

Substrate recognition by norovirus polymerase: microsecond molecular dynamics study

Molecular dynamics simulations of complexes between Norwalk virus RNA dependent RNA polymerase and its natural CTP and 2dCTP (both containing the O5'-C5'-C4'-O4' sequence of atoms bridging the triphosphate and sugar moiety) or modified coCTP (C5'-O5'-C4'-O4'), cocCTP (C5'-O5'-C4'-C4'') substrates were produced by means of CUDA programmable graphical processing units and the ACEMD software package. It enabled us to gain microsecond MD trajectories clearly showing that similar nucleoside triphosphates can bind surprisingly differently into the active site of the Norwalk virus RNA dependent RNA polymerase. It corresponds to their different modes of action (CTP-substrate, 2dCTP-poor substrate, coCTP-chain terminator, cocCTP-inhibitor). Moreover, extremely rare events-as repetitive pervasion of Arg182 into a potentially reaction promoting arrangement-were captured.

How calcium inhibits the magnesium-dependent kinase gsk3β: a molecular simulation study

Glycogen synthase kinase 3β (GSK3β) is a ubiquitous serine/threonine kinase that plays a pivotal role in many biological processes. GSK3β catalyzes the transfer of γ-phosphate of ATP to the unique substrate Ser/Thr residues with the assistance of two natural activating cofactors Mg(2+). Interestingly, the biological observation reveals that a non-native Ca(2+) ion can inhibit the GSK3β catalytic activity. Here, the inhibitory mechanism of GSK3β by the displacement of native Mg(2+) at site 1 by Ca(2+) was investigated by means of 80 ns comparative molecular dynamics (MD) simulations of the GSK3β···Mg(2+)-2/ATP/Mg(2+) -1 and GSK3β···Mg(2+)-2/ATP/Ca(2+)-1 systems. MD simulation results revealed that using the AMBER point charge model force field for Mg(2+) was more appropriate in the reproduction of the active site architectural characteristics of GSK3β than using the magnesium-cationic dummy atom model force field. Compared with the native Mg(2+) bound system, the misalignment of the critical triphosphate moiety of ATP, the erroneous coordination environments around the Mg(2+) ion at site 2, and the rupture of the key hydrogen bond between the invariant Lys85 and the ATP O(β2) atom in the Ca(2+) substituted system were observed in the MD simulation due to the Ca(2+) ion in active site in order to achieve its preferred sevenfold coordination geometry, which adequately abolish the enzymatic activity. The obtained results are valuable in understanding the possible mechanism by why Ca(2+) inhibits the GSK3β activity and also provide insights into the mechanism of Ca(2+) inhibition in other structurally related protein kinases.

Mapping L1 ligase ribozyme conformational switch

L1 ligase (L1L) molecular switch is an in vitro optimized synthetic allosteric ribozyme that catalyzes the regioselective formation of a 5'-to-3' phosphodiester bond, a reaction for which there is no known naturally occurring RNA catalyst. L1L serves as a proof of principle that RNA can catalyze a critical reaction for prebiotic RNA self-replication according to the RNA world hypothesis. L1L crystal structure captures two distinct conformations that differ by a reorientation of one of the stems by around 80Å and are presumed to correspond to the active and inactive state, respectively. It is of great interest to understand the nature of these two states in solution and the pathway for their interconversion. In this study, we use explicit solvent molecular simulation together with a novel enhanced sampling method that utilizes concepts from network theory to map out the conformational transition between active and inactive states of L1L. We find that the overall switching mechanism can be described as a three-state/two-step process. The first step involves a large-amplitude swing that reorients stem C. The second step involves the allosteric activation of the catalytic site through distant contacts with stem C. Using a conformational space network representation of the L1L switch transition, it is shown that the connection between the three states follows different topographical patterns: the stem C swing step passes through a narrow region of the conformational space network, whereas the allosteric activation step covers a much wider region and a more diverse set of pathways through the network.

On the active site of mononuclear B1 metallo β-lactamases: a computational study

Metallo-β-lactamases (MβLs) are Zn(II)-based bacterial enzymes that hydrolyze β-lactam antibiotics, hampering their beneficial effects. In the most relevant subclass (B1), X-ray crystallography studies on the enzyme from Bacillus Cereus point to either two zinc ions in two metal sites (the so-called '3H' and 'DCH' sites) or a single Zn(II) ion in the 3H site, where the ion is coordinated by Asp120, Cys221 and His263 residues. However, spectroscopic studies on the B1 enzyme from B. Cereus in the mono-zinc form suggested the presence of the Zn(II) ion also in the DCH site, where it is bound to an aspartate, a cysteine, a histidine and a water molecule. A structural model of this enzyme in its DCH mononuclear form, so far lacking, is therefore required for inhibitor design and mechanistic studies. By using force field based and mixed quantum-classical (QM/MM) molecular dynamics (MD) simulations of the protein in aqueous solution we constructed such structural model. The geometry and the H-bond network at the catalytic site of this model, in the free form and in complex with two common β-lactam drugs, is compared with experimental and theoretical findings of CphA and the recently solved crystal structure of new B2 MβL from Serratia fonticola (Sfh-I). These are MβLs from the B2 subclass, which features an experimentally well established mono-zinc form, in which the Zn(II) is located in the DCH site. From our simulations the εεδ and δεδ protomers emerge as possible DCH mono-zinc reactive species, giving a novel contribution to the discussion on the MβL reactivity and to the drug design process.