Refinement of the Cornell et al. Nucleic Acids Force Field Based on Reference Quantum Chemical Calculations of Glycosidic Torsion Profiles

Conformational plasticity of SpyCas9 induced by AcrIIA4 and AcrIIA2: Insights from molecular dynamics simulation

CRISPR-Cas9 systems constitute bacterial adaptive immune systems that protect against phage infections. Bacteriophages encode anti-CRISPR proteins (Acrs) that mitigate the bacterial immune response. However, the structural basis for their inhibitory actions from a molecular perspective remains elusive. In this study, through microsecond atomistic molecular dynamics simulations, we demonstrated the remarkable flexibility of Streptococcus pyogenes Cas9 (SpyCas9) and its conformational adaptability during interactions with AcrIIA4 and AcrIIA2. Specifically, we demonstrated that the binding of AcrIIA4 and AcrIIA2 to SpyCas9 induces a conformational rearrangement that causes spatial separation between the nuclease and cleavage sites, thus making the endonuclease inactive. This separation disrupts the transmission of signals between the protospacer adjacent motif recognition and nuclease domains, thereby impeding the efficient processing of double-stranded DNA. The simulation also reveals that AcrIIA4 and AcrIIA2 cause different structural variations of SpyCas9. Our research illuminates the precise mechanisms underlying the suppression of SpyCas9 by AcrIIA4 and AcrIIA2, thus presenting new possibilities for controlling genome editing with higher accuracy.

In silico analysis of aptamer-RNA conjugate interactions with human transferrin receptor

The human transmembrane protein Transferrin Receptor-1 is regarded as a promising target for the systemic delivery of therapeutic agents, particularly of nucleic acid therapeutics, such as short double stranded RNAs. This ubiquitous receptor is involved in cellular iron uptake, keeping intracellular homeostasis. It is overexpressed in multiple cancer cell types and is internalized via clathrin-mediated endocytosis. In previous studies, a human transferrin receptor-1 RNA aptamer, identified as TR14 ST1-3, was shown to be capable of effectively internalizing into cells in culture and to deliver small, double stranded RNAs in vitro and in vivo, via systemic administration. To understand, at the molecular level, the aptamer binding to the receptor and the impact of conjugation with the therapeutic RNA, a multi-level in silico protocol was employed, including protein-aptamer docking, molecular dynamics simulations and free energy calculations. The competition for the binding pocket, between the aptamer and the natural ligand human Transferrin, was also evaluated. The results show that the aptamer binds to the same region as Transferrin, with residues from the helical domain showing a critical role. Moreover, the conjugation to the therapeutic RNA, was shown not to affect aptamer binding. Overall, this study provides an atomic-level understanding of aptamer association to human Transferrin Receptor-1 and of its conjugation with a short model-therapeutic RNA, providing also important clues for futures studies aiming to deliver other oligonucleotide-based therapeutics via Transferrin Receptor.

Atomistic Insights into Sequence-Mediated Spontaneous Association of Short RNA Chains

RNA-RNA association and phase separation appear to be essential for the assembly of stress granules and underlie RNA foci formation in repeat expansion disorders. RNA molecules are found to play a significant role in gene-regulatory functions via condensate formation among themselves or with RNA-binding proteins. The interplay between driven versus spontaneous processes is likely to be an important factor for controlling the formation of RNA-mediated biomolecular condensate. However, the sequence-specific interactions and molecular mechanisms that drive the spontaneous RNA-RNA association and help to form RNA-mediated phase-separated condensate remain unclear. With microseconds-long atomistic molecular simulations here, we report how essential aspects of RNA chains, namely, base composition, metal ion binding, and hydration properties, contribute to the association of the series of simplest biologically relevant homopolymeric and heteropolymeric short RNA chains. We show that spontaneous processes make the key contributions governed by the sequence-intrinsic properties of RNA chains, where the definite roles of base-specific hydrogen bonding and stacking interactions are prominent in the association of the RNA chains. Purine versus pyrimidine contents of RNA chains can directly influence the association properties of RNA chains by modulating hydrogen bonding and base stacking interactions. This study determines the impact of ionic environment in sequence-specific spontaneous association of short RNA chains, hydration features, and base-specific interactions of Na+, K+, and Mg2+ ions with RNA chains.

Self-assembly and condensation of intermolecular poly(UG) RNA quadruplexes

Poly(UG) or 'pUG' dinucleotide repeats are highly abundant sequences in eukaryotic RNAs. In Caenorhabditis elegans, pUGs are added to RNA 3' ends to direct gene silencing within Mutator foci, a germ granule condensate. Here, we show that pUG RNAs efficiently self-assemble into gel condensates through quadruplex (G4) interactions. Short pUG sequences form right-handed intermolecular G4s (pUG G4s), while longer pUGs form left-handed intramolecular G4s (pUG folds). We determined a 1.05 Å crystal structure of an intermolecular pUG G4, which reveals an eight stranded G4 dimer involving 48 nucleotides, 7 different G and U quartet conformations, 7 coordinated potassium ions, 8 sodium ions and a buried water molecule. A comparison of the intermolecular pUG G4 and intramolecular pUG fold structures provides insights into the molecular basis for G4 handedness and illustrates how a simple dinucleotide repeat sequence can form complex structures with diverse topologies.

Mechanism of Nucleic Acid Phosphodiester Bond Cleavage by Human Endonuclease V: MD and QM/MM Calculations Reveal a Versatile Metal Dependence

Human endonuclease V (EndoV) catalytically removes deaminated nucleobases by cleaving the phosphodiester bond as part of RNA metabolism. Despite being implicated in several diseases (cancers, cardiovascular diseases, and neurological disorders) and potentially being a useful tool in biotechnology, details of the human EndoV catalytic pathway remain unclear due to limited experimental information beyond a crystal structure of the apoenzyme and select mutational data. Since a mechanistic understanding is critical for further deciphering the central roles and expanding applications of human EndoV in medicine and biotechnology, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations were used to unveil the atomistic details of the catalytic pathway. Due to controversies surrounding the number of metals required for nuclease activity, enzyme-substrate models with different numbers of active site metals and various metal-substrate binding configurations were built based on structural data for other nucleases. Subsequent MD simulations revealed the structure and stability of the human EndoV-substrate complex for a range of active site metal binding architectures. Four unique pathways were then characterized using QM/MM that vary in metal number (one versus two) and modes of substrate coordination [direct versus indirect (water-mediated)], with several mechanisms being fully consistent with experimental structural, kinetic, and mutational data for related nucleases, including members of the EndoV family. Beyond uncovering key roles for several active site amino acids (D240 and K155), our calculations highlight that while one metal is essential for human EndoV activity, the enzyme can benefit from using two metals due to the presence of two suitable metal binding sites. By directly comparing one- versus two-metal-mediated P-O bond cleavage reactions within the confines of the same active site, our work brings a fresh perspective to the "number of metals" controversy.

Magnesium Regulates RNA Ring Dynamics and Folding in Subgenomic Flaviviral RNA

Mosquito-borne flaviviruses including dengue, Zika, yellow fever, and regional encephalitis produce a large amount of short subgenomic flaviviral RNAs during infection. A segment of these RNAs named as xrRNA1 features a multi-pseudoknot (PK)-associated structure, which resists the host cell enzyme (XRN1) from degrading the viral RNA. We investigate how this long-range RNA PK folds in the presence of counterions, specifically in a mix of monovalent (K+) and divalent (Mg2+) salts at physiological concentrations. In this study, we use extensive explicit solvent molecular dynamics (MD) simulations to characterize the RNA ion environment of the folded RNA conformation, as determined by the crystal structure. This allowed us to identify the precise locations of various coordinated RNA-Mg2+ interactions, including inner-sphere/chelated and outer-sphere coordinated Mg2+. Given that RNA folding involves large-scale conformational changes, making it challenging to explore through classical MD simulations, we investigate the folding mechanism of xrRNA1 using an all-atom structure-based RNA model with a hybrid implicit-explicit treatment of the ion environment via the dynamic counterion condensation model, both with and without physiological Mg2+ concentration. The study reveals potential folding pathways for this xrRNA1, which is consistent with the results obtained from optical tweezer experiments. The equilibrium and free energy simulations both capture a dynamic equilibrium between the ring-open and ring-close states of the RNA, driven by a long-range PK interaction. Free energy calculations reveal that with the addition of Mg2+ ions, the equilibrium shifts more toward the ring-close state. A detailed analysis of the free energy pathways and ion-mediated contact probability map highlights the critical role of Mg2+ in bridging G50 and A33. This Mg2+-mediated connection helps form the long-range PK which in turn controls the transition between the ring-open and ring-close states. The study underscores the critical role of Mg2+ in the RNA folding transition, highlighting specific locations of Mg2+ contributing to the stabilization of long-range PK connections likely to enhance the robustness of Xrn1 resistance of flaviviral xrRNAs.

Unveiling Putative Excited State and Transmission of Binding Information in the Fluoride Riboswitch

Riboswitches regulate downstream gene expression by binding to specific small molecules or ions with multiple mechanisms to transfer the binding information. In the case of the fluoride riboswitch, the transcription termination signal is conveyed through a transient excited state (ES). In this work, we performed conventional molecular dynamics (MD) simulations, totaling 180 μs, to obtain the ES structure and investigate the mechanism underlying information transmission in Mg2+/F- binding within the fluoride riboswitch aptamer. The Mg2+/F- binding pocket exhibits various conformations in its apo form. A series of ES structures were extracted from the MD trajectories of the apo form. The dynamics of the Mg2+/F- binding pocket influenced key pair A40-U48 in ES structures. The pathway connecting the binding pocket to the pair involves interactions between the phosphate groups of U7 and G8 and the nucleobases of G8-C47-U48. Our work presents a structural ensemble of the ES and elucidates a pathway for transferring Mg2+/F- binding information, thereby facilitating the understanding of how the holo-like apo state achieves transcriptional repression.

The molecular mechanism for TERRA recruitment and annealing to telomeres

Telomeric repeat containing RNA (TERRA) is a noncoding RNA that is transcribed from telomeres. Previous study showed that TERRA trans anneals by invading into the telomeric duplex to form an R-loop in mammalian cells. Here, we elucidate the molecular mechanism underlying TERRA recruitment and invasion into telomeres in the context of shelterin proteins, RAD51 and RNase H using single molecule (sm) assays. We demonstrate that TERRA trans annealing into telomeric DNA exhibits dynamic movement that is stabilized by TRF2. TERRA annealing to the telomeric duplex results in the formation of a stable triplex structure which differs from a conventional R-loop. We identified that the presence of a sub-telomeric DNA and a telomeric overhang in the form of a G-quadruplex significantly enhances TERRA annealing to telomeric duplex. We also demonstrate that RAD51-TERRA complex invades telomere duplex more efficiently than TERRA alone. Additionally, TRF2 increases TERRA affinity to telomeric duplex and protects it from RNase H digestion. In contrast, TRF1 represses TERRA annealing to telomeric duplex and fails to provide protection against RNase H digestion. Our findings provide an in-depth molecular mechanism underpinning TERRA recruitment and annealing to the telomere.

RudS: bacterial desulfidase responsible for tRNA 4-thiouridine de-modification

In this study, we present an extensive analysis of a widespread group of bacterial tRNA de-modifying enzymes, dubbed RudS, which consist of a TudS desulfidase fused to a Domain of Unknown Function 1722 (DUF1722). RudS enzymes exhibit specific de-modification activity towards the 4-thiouridine modification (s4U) in tRNA molecules, as indicated by our experimental findings. The heterologous overexpression of RudS genes in Escherichia coli significantly reduces the tRNA 4-thiouridine content and diminishes UVA-induced growth delay, indicating the enzyme's role in regulating photosensitive tRNA s4U modification. Through a combination of protein modeling, docking studies, and molecular dynamics simulations, we have identified amino acid residues involved in catalysis and tRNA binding. Experimental validation through targeted mutagenesis confirms the TudS domain as the catalytic core of RudS, with the DUF1722 domain facilitating tRNA binding in the anticodon region. Our results suggest that RudS tRNA modification eraser proteins may play a role in regulating tRNA during prokaryotic stress responses.

N-terminal domain of polypyrimidine-tract binding protein is a dynamic folding platform for adaptive RNA recognition

The N-terminal RNA recognition motif domain (RRM1) of polypyrimidine tract binding protein (PTB) forms an additional C-terminal helix α3, which docks to one edge of the β-sheet upon binding to a stem-loop RNA containing a UCUUU pentaloop. Importantly, α3 does not contact the RNA. The α3 helix therefore represents an allosteric means to regulate the conformation of adjacent domains in PTB upon binding structured RNAs. Here we investigate the process of dynamic adaptation by stem-loop RNA and RRM1 using NMR and MD in order to obtain mechanistic insights on how this allostery is achieved. Relaxation data and NMR structure determination of the free protein show that α3 is partially ordered and interacts with the domain transiently. Stem-loop RNA binding quenches fast time scale dynamics and α3 becomes ordered, however microsecond dynamics at the protein-RNA interface is observed. MD shows how RRM1 binding to the stem-loop RNA is coupled to the stabilization of the C-terminal helix and helps to transduce differences in RNA loop sequence into changes in α3 length and order. IRES assays of full length PTB and a mutant with altered dynamics in the α3 region show that this dynamic allostery influences PTB function in cultured HEK293T cells.

Physiological significance of the two isoforms of initiator tRNAs in Escherichia coli

Escherichia coli possesses four initiator tRNA (i-tRNA) genes, three of which are present together as metZWV and the fourth one as metY. In E. coli B, all four genes (metZWV and metY) encode i-tRNAfMet1, in which the G at position 46 is modified to m7G46 by TrmB (m7G methyltransferase). However, in E. coli K, because of a single-nucleotide polymorphism, metY encodes a variant, i-tRNAfMet2, having an A in place of m7G46. We generated E. coli strains to explore the importance of this polymorphism in i-tRNAs. The strains were sustained either on metYA46 (metY of E. coli K origin encoding i-tRNAfMet2) or its derivative metYG46 (encoding i-tRNAfMet1) in single (chromosomal) or plasmid-borne copies. We show that the strains sustained on i-tRNAfMet1 have a growth fitness advantage over those sustained on i-tRNAfMet2. The growth fitness advantages are more pronounced for the strains sustained on i-tRNAfMet1 in nutrient-rich media than in nutrient-poor media. The growth fitness of the strains correlates well with the relative stabilities of the i-tRNAs in vivo. Furthermore, the atomistic molecular dynamics simulations support the higher stability of i-tRNAfMet1 than that of i-tRNAfMet2. The stability of i-tRNAfMet1 remains unaffected upon the deletion of TrmB. These studies highlight how metYG46 and metYA46 alleles might influence the growth fitness of E. coli under certain nutrient-limiting conditions.

IMPORTANCE

Escherichia coli harbors four initiator tRNA (i-tRNA) genes: three of these at metZWV and the fourth one at metY loci. In E. coli B, all four genes encode i-tRNAfMet1. In E. coli K, because of a single-nucleotide polymorphism, metY encodes a variant, i-tRNAfMet2, having an A in place of G at position 46 of i-tRNA sequence in metY. We show that G46 confers stability to i-tRNAfMet1. The strains sustained on i-tRNAfMet1 have a growth fitness advantage over those sustained on i-tRNAfMet2. Strains harboring metYG46 (B mimic) or metYA46 (K mimic) show that while in the nutrient-rich media, the K mimic is outcompeted rapidly; in the nutrient-poor medium, the K mimic is outcompeted less rapidly.

Molecular Simulations to Investigate the Impact of N6-Methylation in RNA Recognition: Improving Accuracy and Precision of Binding Free Energy Prediction

N6-Methyladenosine (m6A) is a prevalent RNA post-transcriptional modification that plays crucial roles in RNA stability, structural dynamics, and interactions with proteins. The YT521-B (YTH) family of proteins, which are notable m6A readers, functions through its highly conserved YTH domain. Recent structural investigations and molecular dynamics (MD) simulations have shed light on the mechanism of recognition of m6A by the YTHDC1 protein. Despite advancements, using MD to predict the stabilization induced by m6A on the free energy of binding between RNA and YTH proteins remains challenging due to inaccuracy of the employed force field and limited sampling. For instance, simulations often fail to sufficiently capture the hydration dynamics of the binding pocket. This study addresses these challenges through an innovative methodology that integrates metadynamics, alchemical simulations, and force-field refinement. Importantly, our research identifies hydration of the binding pocket as giving only a minor contribution to the binding free energy and emphasizes the critical importance of precisely tuning force-field parameters to experimental data. By employing a fitting strategy built on alchemical calculations, we refine the m6A partial charge parameters, thereby enabling the simultaneous reproduction of N6 methylation on both the protein binding free energy and the thermodynamic stability of nine RNA duplexes. Our findings underscore the sensitivity of binding free energies to partial charges, highlighting the necessity for thorough parametrization and validation against experimental observations across a range of structural contexts.

Activation of automethylated PRC2 by dimerization on chromatin

Polycomb repressive complex 2 (PRC2) is an epigenetic regulator that trimethylates lysine 27 of histone 3 (H3K27me3) and is essential for embryonic development and cellular differentiation. H3K27me3 is associated with transcriptionally repressed chromatin and is established when PRC2 is allosterically activated upon methyl-lysine binding by the regulatory subunit EED. Automethylation of the catalytic subunit enhancer of zeste homolog 2 (EZH2) stimulates its activity by an unknown mechanism. Here, we show that human PRC2 forms a dimer on chromatin in which an inactive, automethylated PRC2 protomer is the allosteric activator of a second PRC2 that is poised to methylate H3 of a substrate nucleosome. Functional assays support our model of allosteric trans-autoactivation via EED, suggesting a previously unknown mechanism mediating context-dependent activation of PRC2. Our work showcases the molecular mechanism of auto-modification-coupled dimerization in the regulation of chromatin-modifying complexes.

DMDA-PatA mediates RNA sequence-selective translation repression by anchoring eIF4A and DDX3 to GNG motifs

Small-molecule compounds that elicit mRNA-selective translation repression have attracted interest due to their potential for expansion of druggable space. However, only a limited number of examples have been reported to date. Here, we show that desmethyl desamino pateamine A (DMDA-PatA) represses translation in an mRNA-selective manner by clamping eIF4A, a DEAD-box RNA-binding protein, onto GNG motifs. By systematically comparing multiple eIF4A inhibitors by ribosome profiling, we found that DMDA-PatA has unique mRNA selectivity for translation repression. Unbiased Bind-n-Seq reveals that DMDA-PatA-targeted eIF4A exhibits a preference for GNG motifs in an ATP-independent manner. This unusual RNA binding sterically hinders scanning by 40S ribosomes. A combination of classical molecular dynamics simulations and quantum chemical calculations, and the subsequent development of an inactive DMDA-PatA derivative reveals that the positive charge of the tertiary amine on the trienyl arm induces G selectivity. Moreover, we identified that DDX3, another DEAD-box protein, is an alternative DMDA-PatA target with the same effects on eIF4A. Our results provide an example of the sequence-selective anchoring of RNA-binding proteins and the mRNA-selective inhibition of protein synthesis by small-molecule compounds.

Engineering salt-tolerant Cas12f1 variants for gene-editing applications

CRISPR has revolutionized the field of genome editing in life sciences by serving as a versatile and state-of-the-art tool. Cas12f1 is a small nuclease of the bacterial immunity CRISPR system with an ideal size for cellular delivery, in contrast to CRISPR-associated (Cas) proteins like Cas9 or Cas12. However, Cas12f1 works best at low salt concentrations. In this study, we find that the plasticity of certain Cas12f1 regions (K196-Y202 and I452-L515) is negatively affected by increased salt concentrations. On this line, key protein domains (REC1, WED, Nuc, lid) that are involved in the DNA-target recognition and the activation of the catalytic RuvC domain are in turn also affected. We suggest that salt concentration should be taken in to consideration for activity assessments of Cas engineered variants, especially if the mutations are on the protospacer adjacent motif interacting domain. The results can be exploited for the engineering of Cas variants and the assessment of their activity at varying salt concentrations. We propose that the K198Q mutation can restore at great degree the compromised plasticity and could potentially lead to salt-tolerant Cas12f1 variants. The methodology can be also employed for the study of biomolecules in terms of their salinity tolerance.Communicated by Ramaswamy H. Sarma.

The effect of the loop on the thermodynamic and kinetic of single base pair in pseudoknot

RNA pseudoknots are RNA molecules with specialized three-dimensional structures that play important roles in various biological processes. To understand the functions and mechanisms of pseudoknots, it is essential to elucidate their structures and folding pathways. The most fundamental step in RNA folding is the opening and closing of a base pair. The effect of flexible loops on the base pair in pseudoknots remains unclear. In this work, we use molecular dynamics simulations and Markov state model to study the configurations, thermodynamic and kinetic of single base pair in pseudoknots. We find that the presence of the loop leads to a trap state. In addition, the rate-limiting step for the formation of base pair is the disruption of the trap state, rather than the open state to the closed state, which is quite different from the previous studies on non-pseudoknot RNA. For the thermodynamic parameters in pseudoknots, we find that the entropy difference upon opening the base pair between this simulation and the nearest-neighbor model results from the different entropy of different lengths of loop in solution. The thermodynamic parameters of the stack in pseudoknot are close to the nearest-neighbor parameters. The bases on the loop have different distribution patterns in different states, and the slow transition states of the loop are determined by the orientation of the bases.

Adjusting the Energy Profile for CH-O Interactions Leads to Improved Stability of RNA Stem-Loop Structures in MD Simulations

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules

RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.

SPRank─A Knowledge-Based Scoring Function for RNA-Ligand Pose Prediction and Virtual Screening

The growing interest in RNA-targeted drugs underscores the need for computational modeling of interactions between RNA molecules and small compounds. Having a reliable scoring function for RNA-ligand interactions is essential for effective computational drug screening. An ideal scoring function should not only predict the native pose for ligand binding but also rank the affinity of the binding for different ligands. However, existing scoring functions are primarily designed to predict the native binding modes for a given RNA-ligand pair and have not been thoroughly assessed for virtual screening purposes. In this paper, we introduce SPRank, a combination of machine-learning and knowledge-based scoring functions developed through a weighted iterative approach, specifically designed to tackle both binding mode prediction and virtual screening challenges. Our approach incorporates third-party docking software, such as rDock and AutoDock Vina, to sample flexible ligands against an ensemble of RNA structures, capturing the conformational flexibility of both the RNA and the ligand. Through rigorous testing, SPRank demonstrates improved performance compared to the tested scoring functions across four test sets comprising 122, 42, 55, and 71 nucleic acid-ligand complexes. Furthermore, SPRank exhibits improved performance in virtual screening tests targeting the HIV-1 TAR ensemble, which highlights its advantage in drug discovery. These results underscore the advantages of SPRank as a potentially promising tool for the RNA-targeted drug design. The source code of SPRank and the data sets are freely accessible at https://github.com/Vfold-RNA/SPRank.

Machine-learned molecular mechanics force fields from large-scale quantum chemical data

The development of reliable and extensible molecular mechanics (MM) force fields-fast, empirical models characterizing the potential energy surface of molecular systems-is indispensable for biomolecular simulation and computer-aided drug design. Here, we introduce a generalized and extensible machine-learned MM force field, espaloma-0.3, and an end-to-end differentiable framework using graph neural networks to overcome the limitations of traditional rule-based methods. Trained in a single GPU-day to fit a large and diverse quantum chemical dataset of over 1.1 M energy and force calculations, espaloma-0.3 reproduces quantum chemical energetic properties of chemical domains highly relevant to drug discovery, including small molecules, peptides, and nucleic acids. Moreover, this force field maintains the quantum chemical energy-minimized geometries of small molecules and preserves the condensed phase properties of peptides and folded proteins, self-consistently parametrizing proteins and ligands to produce stable simulations leading to highly accurate predictions of binding free energies. This methodology demonstrates significant promise as a path forward for systematically building more accurate force fields that are easily extensible to new chemical domains of interest.

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes

Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.

Accurately Modeling RNA Stem-Loops in an Implicit Solvent Environment

Ribonucleic acid (RNA) molecules can adopt a variety of secondary and tertiary structures in solution, with stem-loops being one of the more common motifs. Here, we present a systematic analysis of 15 RNA stem-loop sequences simulated with molecular dynamics simulations in an implicit solvent environment. Analysis of RNA cluster ensembles showed that the stem-loop structures can generally adopt the A-form RNA in the stem region. Loop structures are more sensitive, and experimental structures could only be reproduced with modification of CH···O interactions in the force field, combined with an implicit solvent nonpolar correction to better model base stacking interactions. Accurately modeling RNA with current atomistic physics-based models remains challenging, but the RNA systems studied herein may provide a useful benchmark set for testing other RNA modeling methods in the future.

Novel Inhibitors of SARS-CoV-2 RNA Identified through Virtual Screening

We currently lack antivirals for most human viruses. In a quest for new molecules, focusing on viral RNA, instead of viral proteins, can represent a promising strategy. In this study, new inhibitors were identified starting from a published crystal structure of the tertiary SARS-CoV-2 RNA involved in the -1 programmed ribosomal frameshift. The pseudoknot structure was refined, and a virtual screening was performed using the repository of binders to the nucleic acid library, taking into consideration RNA flexibility. Hit compounds were validated against the wild-type virus and with a dual-luciferase assay measuring the frameshift efficiency. Several active molecules were identified. Our study reveals new inhibitors of SARS-CoV-2 but also highlights the feasibility of targeting RNA starting from virtual screening, a strategy that could be broadly applied to drug development.

Expanding the Available RNA Labeling Toolbox With CutA Nucleotidyltransferase for Efficient Transcript Labeling with Purine and Pyrimidine Nucleotide Analogs

RNA labeling is an invaluable tool for investigation of the function and localization of nucleic acids. Labels are commonly incorporated into 3' end of RNA and the primary enzyme used for this purpose is RNA poly(A) polymerase (PAP), which belongs to the class of terminal nucleotidyltransferases (NTases). However, PAP preferentially adds ATP analogs, thus limiting the number of available substrates. Here, we report the use of another NTase, CutA from the fungus Thielavia terrestris. Using this enzyme, we were able to incorporate into the 3' end of RNA not only purine analogs, but also pyrimidine analogs. We engaged strain-promoted azide-alkyl cycloaddition (SPAAC) to obtain fluorescently labeled or biotinylated transcripts from RNAs extended with azide analogs by CutA. Importantly, modified transcripts retained their biological properties. Furthermore, fluorescently labeled mRNAs were suitable for visualization in cultured mammalian cells. Finally, we demonstrate that either affinity studies or molecular dynamic (MD) simulations allow for rapid screening of NTase substrates, what opens up new avenues in the search for the optimal substrates for this class of enzymes.

Deciphering the Role of Trehalose in Chroococcidiopsis sp. 029's High-Desiccation Resistance: Sequence Determination, Structural Modelling and Simulative Analysis of the 30S Ribosomal Subunit

Desert strains of the genus Chroococcidiopsis are among the most desiccation-resistant cyanobacteria capable of anhydrobiosis. The accumulation of two sugars, sucrose and trehalose, facilitates the entrance of anhydrobiotes into a reversible state of dormancy by stabilizing cellular components upon water removal. This study aimed to evaluate, at the atomistic level, the role of trehalose in desiccation resistance by using as a model system the 30S ribosomal subunit of the desert cyanobacterium Chroococcidiopsis sp. 029. Molecular dynamic simulations provided atomistic evidence regarding its protective role on the 30S molecular structure. Trehalose forms an enveloping shell around the ribosomal subunit and stabilizes the structures through a network of direct interactions. The simulation confirmed that trehalose actively interacts with the 30S ribosomal subunit and that, by replacing water molecules, it ensures ribosomal structural integrity during desiccation, thus enabling protein synthesis to be carried out upon rehydration.

Cas12a domain flexibility guides R-loop formation and forces RuvC resetting

The specific nature of CRISPR-Cas12a makes it a desirable RNA-guided endonuclease for biotechnology and therapeutic applications. To understand how R-loop formation within the compact Cas12a enables target recognition and nuclease activation, we used cryo-electron microscopy to capture wild-type Acidaminococcus sp. Cas12a R-loop intermediates and DNA delivery into the RuvC active site. Stages of Cas12a R-loop formation-starting from a 5-bp seed-are marked by distinct REC domain arrangements. Dramatic domain flexibility limits contacts until nearly complete R-loop formation, when the non-target strand is pulled across the RuvC nuclease and coordinated domain docking promotes efficient cleavage. Next, substantial domain movements enable target strand repositioning into the RuvC active site. Between cleavage events, the RuvC lid conformationally resets to occlude the active site, requiring re-activation. These snapshots build a structural model depicting Cas12a DNA targeting that rationalizes observed specificity and highlights mechanistic comparisons to other class 2 effectors.

Simulation-Guided Conformational Space Exploration to Assess Reactive Conformations of a Ribozyme

Self-splicing ribozymes are small ribonucleic acid (RNA) enzymes that catalyze their own cleavage through a transphosphoesterification reaction. While this process is involved in some specific steps of viral RNA replication and splicing, it is also of importance in the context of the (putative) first autocatalytic RNA-based systems that could have preceded the emergence of modern life. The uncatalyzed phosphoester bond formation is thermodynamically very unfavorable, and many experimental studies have focused on understanding the molecular features of catalysis in these ribozymes. However, chemical reaction paths are short-lived and not easily characterized by experimental approaches, so molecular simulation approaches appear as an ideal tool to unveil the molecular details of the reaction. Here, we focus on the model hairpin ribozyme. We show that identifying a relevant initial conformation for reactivity studies, which is frequently overlooked in mixed quantum-classical studies that predominantly concentrate on the chemical reaction itself, can be highly challenging. These challenges stem from limitations in both available experimental structures (which are chemically altered to prevent self-cleavage) and the accuracy of force fields, together with the necessity for comprehensive sampling. We show that molecular dynamics simulations, combined with extensive conformational phase space exploration with Hamiltonian replica-exchange simulations, enable us to characterize the relevant conformational basins of the minimal hairpin ribozyme in the ligated state prior to self-cleavage. We find that what is usually considered a canonical reactive conformation with active site geometries and hydrogen-bond patterns that are optimal for the addition-elimination reaction with general acid/general base catalysis is metastable and only marginally populated. The thermodynamically stable conformation appears to be consistent with the expectations of a mechanism that does not require the direct participation of ribozyme residues in the reaction. While these observations may suffer from forcefield inaccuracies, all investigated forcefields lead to the same conclusions upon proper sampling, contrasting with previous investigations on shorter timescales suggesting that at least one reparametrization of the Amber99 forcefield allowed to stabilize aligned active site conformations. Our study demonstrates that identifying the most pertinent reactant state conformation holds equal importance alongside the accurate determination of the thermodynamics and kinetics of the chemical steps of the reaction.

Excited-Ground-State Transition of the RNA Strand Slippage Mechanism Captured by the Base-Specific Force Field

Excited-ground-state transition and strand slippage of RNA play key roles in transcription and translation of central dogma. Due to limitation of current experimental techniques, the dynamic structure ensembles of RNA remain inadequately understood. Molecular dynamics simulations offer a promising complementary approach, whose accuracy depends on the force field. Here, we develop the new version of RNA base-specific force field (BSFF2) to address underestimation of base pairing stability and artificial backbone conformations. Extensive evaluations on typical RNA systems have comprehensively confirmed the accuracy of BSFF2. Furthermore, BSFF2 demonstrates exceptional efficiency in de novo folding of tetraloops and reproducing base pair reshuffling transition between RNA excited and ground states. Then, we explored the RNA strand slippage mechanism with BSFF2. We conducted a comprehensive three-dimensional structural investigation into the strand slippage of the most complex r(G4C2)9 repeat element and presented the molecular details in the dynamic transition along with the underlying mechanism. Our results of capturing the strand slippage, excited-ground transition, de novo folding, and simulations for various typical RNA motifs indicate that BSFF2 should be one of valuable tools for dynamic conformation research and structure prediction of RNA, and a future contribution to RNA-targeted drug design as well as RNA therapy development.

Motif-VI loop acts as a nucleotide valve in the West Nile Virus NS3 Helicase

The Orthoflavivirus NS3 helicase (NS3h) is crucial in virus replication, representing a potential drug target for pathogenesis. NS3h utilizes nucleotide triphosphate (ATP) for hydrolysis energy to translocate on single-stranded nucleic acids, which is an important step in the unwinding of double-stranded nucleic acids. Intermediate states along the ATP hydrolysis cycle and conformational changes between these states, represent important yet difficult-to-identify targets for potential inhibitors. Extensive molecular dynamics simulations of West Nile virus NS3h+ssRNA in the apo, ATP, ADP+Pi and ADP bound states were used to model the conformational ensembles along this cycle. Energetic and structural clustering analyses depict a clear trend of differential enthalpic affinity of NS3h with ADP, demonstrating a probable mechanism of hydrolysis turnover regulated by the motif-VI loop (MVIL). Based on these results, MVIL mutants (D471L, D471N and D471E) were found to have a substantial reduction in ATPase activity and RNA replication compared to the wild-type. Simulations of the mutants in the apo state indicate a shift in MVIL populations favoring either a closed or open 'valve' conformation, affecting ATP entry or stabilization, respectively. Combining our molecular modeling with experimental evidence highlights a conformation-dependent role for MVIL as a 'valve' for the ATP-pocket, presenting a promising target for antiviral development.

Carbazole Derivatives Binding to Bcl-2 Promoter Sequence G-quadruplex

In this study, we used ultraviolet-visible (UV-Vis), fluorescence, and circular dichroism (CD) techniques, as well as molecular modeling, to probe the interactions between carbazole derivatives and the G-quadruplex structure formed in the promoter region of gene Bcl-2. This gene is a rational target for anticancer therapy due to its high expression in a variety of tumors as well as resistance to chemotherapy-induced apoptosis. We employed a sequence with a specific dual G-to-T mutation that may form a mixed-type hybrid G-quadruplex structure in the Bcl-2 P1 promoter region. The three tested carbazole compounds differing in substitution on the nitrogen atom of carbazole interact with the Bcl-2 G-quadruplex by the same binding mode with the very comparable binding affinities in the order of 105 M-1. During absorption and fluorescence measurements, large changes in the ligand spectra were observed at higher G4 concentrations. The spectrophotometric titration results showed a two-step complex formation between the ligands and the G-quadruplex in the form of initial hypochromicity followed by hyperchromicity with a bathochromic shift. The strong fluorescence enhancement of ligands was observed after binding to the DNA. All of the used analytical techniques, as well as molecular modeling, suggested the π-π interaction between carbazole ligands and a guanine tetrad of the Bcl-2 G-quadruplex. Molecular modeling has shown differences in the interaction between each of the ligands and the tested G-quadruplex, which potentially had an impact on the binding strength.

Identifying small-molecules binding sites in RNA conformational ensembles with SHAMAN

The rational targeting of RNA with small molecules is hampered by our still limited understanding of RNA structural and dynamic properties. Most in silico tools for binding site identification rely on static structures and therefore cannot face the challenges posed by the dynamic nature of RNA molecules. Here, we present SHAMAN, a computational technique to identify potential small-molecule binding sites in RNA structural ensembles. SHAMAN enables exploring the conformational landscape of RNA with atomistic molecular dynamics simulations and at the same time identifying RNA pockets in an efficient way with the aid of probes and enhanced-sampling techniques. In our benchmark composed of large, structured riboswitches as well as small, flexible viral RNAs, SHAMAN successfully identifies all the experimentally resolved pockets and ranks them among the most favorite probe hotspots. Overall, SHAMAN sets a solid foundation for future drug design efforts targeting RNA with small molecules, effectively addressing the long-standing challenges in the field.

Software Infrastructure for Next-Generation QM/MM-ΔMLP Force Fields

We present software infrastructure for the design and testing of new quantum mechanical/molecular mechanical and machine-learning potential (QM/MM-ΔMLP) force fields for a wide range of applications. The software integrates Amber's molecular dynamics simulation capabilities with fast, approximate quantum models in the xtb package and machine-learning potential corrections in DeePMD-kit. The xtb package implements the recently developed density-functional tight-binding QM models with multipolar electrostatics and density-dependent dispersion (GFN2-xTB), and the interface with Amber enables their use in periodic boundary QM/MM simulations with linear-scaling QM/MM particle-mesh Ewald electrostatics. The accuracy of the semiempirical models is enhanced by including machine-learning correction potentials (ΔMLPs) enabled through an interface with the DeePMD-kit software. The goal of this paper is to present and validate the implementation of this software infrastructure in molecular dynamics and free energy simulations. The utility of the new infrastructure is demonstrated in proof-of-concept example applications. The software elements presented here are open source and freely available. Their interface provides a powerful enabling technology for the design of new QM/MM-ΔMLP models for studying a wide range of problems, including biomolecular reactivity and protein-ligand binding.

Structure of an internal loop motif with three consecutive U•U mismatches from stem-loop 1 in the 3'-UTR of the SARS-CoV-2 genomic RNA

The single-stranded RNA genome of SARS-CoV-2 is highly structured. Numerous helical stem-loop structures interrupted by mismatch motifs are present in the functionally important 5'- and 3'-UTRs. These mismatches modulate local helical geometries and feature unusual arrays of hydrogen bonding donor and acceptor groups. However, their conformational and dynamical properties cannot be directly inferred from chemical probing and are difficult to predict theoretically. A mismatch motif (SL1-motif) consisting of three consecutive U•U base pairs is located in stem-loop 1 of the 3'-UTR. We combined NMR-spectroscopy and MD-simulations to investigate its structure and dynamics. All three U•U base pairs feature two direct hydrogen bonds and are as stable as Watson-Crick A:U base pairs. Plasmodium falciparum 25S rRNA contains a triple U•U mismatch motif (Pf-motif) differing from SL1-motif only with respect to the orientation of the two closing base pairs. Interestingly, while the geometry of the outer two U•U mismatches was identical in both motifs the preferred orientation of the central U•U mismatch was different. MD simulations and potassium ion titrations revealed that the potassium ion-binding mode to the major groove is connected to the different preferred geometries of the central base pair in the two motifs.

Targeting the conserved active site of splicing machines with specific and selective small molecule modulators

The self-splicing group II introns are bacterial and organellar ancestors of the nuclear spliceosome and retro-transposable elements of pharmacological and biotechnological importance. Integrating enzymatic, crystallographic, and simulation studies, we demonstrate how these introns recognize small molecules through their conserved active site. These RNA-binding small molecules selectively inhibit the two steps of splicing by adopting distinctive poses at different stages of catalysis, and by preventing crucial active site conformational changes that are essential for splicing progression. Our data exemplify the enormous power of RNA binders to mechanistically probe vital cellular pathways. Most importantly, by proving that the evolutionarily-conserved RNA core of splicing machines can recognize small molecules specifically, our work provides a solid basis for the rational design of splicing modulators not only against bacterial and organellar introns, but also against the human spliceosome, which is a validated drug target for the treatment of congenital diseases and cancers.

Simulations predict preferred Mg2+ coordination in a nonenzymatic primer-extension reaction center

The mechanism by which genetic information was copied prior to the evolution of ribozymes is of great interest because of its importance to the origin of life. The most effective known process for the nonenzymatic copying of an RNA template is primer extension by a two-step pathway in which 2-aminoimidazole-activated nucleotides first react with each other to form an imidazolium-bridged intermediate that subsequently reacts with the primer. Reaction kinetics, structure-activity relationships, and X-ray crystallography have provided insight into the overall reaction mechanism, but many puzzles remain. In particular, high concentrations of Mg2+ are required for efficient primer extension, but the mechanism by which Mg2+ accelerates primer extension remains unknown. By analogy with the mechanism of DNA and RNA polymerases, a role for Mg2+ in facilitating the deprotonation of the primer 3'-hydroxyl is often assumed, but no catalytic metal ion is seen in crystal structures of the primer-extension complex. To explore the potential effects of Mg2+ binding in the reaction center, we performed atomistic molecular dynamics simulations of a series of modeled complexes in which a Mg2+ ion was placed in the reaction center with inner-sphere coordination with different sets of functional groups. Our simulations suggest that coordination of a Mg2+ ion with both O3' of the terminal primer nucleotide and the pro-Sp nonbridging oxygen of the reactive phosphate of an imidazolium-bridged dinucleotide would help to pre-organize the structure of the primer/template substrate complex to favor the primer-extension reaction. Our results suggest that the catalytic metal ion may play an important role in overcoming electrostatic repulsion between a deprotonated O3' and the reactive phosphate of the bridged dinucleotide and lead to testable predictions of the mode of Mg2+ binding that is most relevant to catalysis of primer extension.

Similarities and Differences in Ligand Binding to Protein and RNA Targets: The Case of Riboflavin

It is nowadays clear that RNA molecules can play active roles in several biological processes. As a result, an increasing number of RNAs are gradually being identified as potentially druggable targets. In particular, noncoding RNAs can adopt highly organized conformations that are suitable for drug binding. However, RNAs are still considered challenging targets due to their complex structural dynamics and high charge density. Thus, elucidating relevant features of drug-RNA binding is fundamental for advancing drug discovery. Here, by using Molecular Dynamics simulations, we compare key features of ligand binding to proteins with those observed in RNA. Specifically, we explore similarities and differences in terms of (i) conformational flexibility of the target, (ii) electrostatic contribution to binding free energy, and (iii) water and ligand dynamics. As a test case, we examine binding of the same ligand, namely riboflavin, to protein and RNA targets, specifically the riboflavin (RF) kinase and flavin mononucleotide (FMN) riboswitch. The FMN riboswitch exhibited enhanced fluctuations and explored a wider conformational space, compared to the protein target, underscoring the importance of RNA flexibility in ligand binding. Conversely, a similar electrostatic contribution to the binding free energy of riboflavin was found. Finally, greater stability of water molecules was observed in the FMN riboswitch compared to the RF kinase, possibly due to the different shape and polarity of the pockets.

Splicing analysis of STAT3 tandem donor suggests non-canonical binding registers for U1 and U6 snRNAs

Tandem donor splice sites (5'ss) are unique regions with at least two GU dinucleotides serving as splicing cleavage sites. The Δ3 tandem 5'ss are a specific subclass of 5'ss separated by 3 nucleotides which can affect protein function by inserting/deleting a single amino acid. One 5'ss is typically preferred, yet factors governing particular 5'ss choice are not fully understood. A highly conserved exon 21 of the STAT3 gene was chosen as a model to study Δ3 tandem 5'ss splicing mechanisms. Based on multiple lines of experimental evidence, endogenous U1 snRNA most likely binds only to the upstream 5'ss. However, the downstream 5'ss is used preferentially, and the splice site choice is not dependent on the exact U1 snRNA binding position. Downstream 5'ss usage was sensitive to exact nucleotide composition and dependent on the presence of downstream regulatory region. The downstream 5'ss usage could be best explained by two novel interactions with endogenous U6 snRNA. U6 snRNA enables the downstream 5'ss usage in STAT3 exon 21 by two mechanisms: (i) binding in a novel non-canonical register and (ii) establishing extended Watson-Crick base pairing with the downstream regulatory region. This study suggests that U6:5'ss interaction is more flexible than previously thought.

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

Structure and Internal Dynamics of Short RNA Duplexes Determined by a Combination of Pulsed EPR Methods and MD Simulations

We used EPR spectroscopy to characterize the structure of RNA duplexes and their internal twist, stretch and bending motions. We prepared eight 20-base-pair-long RNA duplexes containing the rigid spin-label Çm, a cytidine analogue, at two positions and acquired orientation-selective PELDOR/DEER data. By using different frequency bands (X-, Q-, G-band), detailed information about the distance and orientation of the labels was obtained and provided insights into the global conformational dynamics of the RNA duplex. We used 19F Mims ENDOR experiments on three singly Çm- and singly fluorine-labeled RNA duplexes to determine the exact position of the Çm spin label in the helix. In a quantitative comparison to MD simulations of RNA with and without Çm spin labels, we found that state-of-the-art force fields with explicit parameterization of the spin label were able to describe the conformational ensemble present in our experiments. The MD simulations further confirmed that the Çm spin labels are excellent mimics of cytidine inducing only small local changes in the RNA structure. Çm spin labels are thus ideally suited for high-precision EPR experiments to probe the structure and, in conjunction with MD simulations, motions of RNA.

Assessing RNA atomistic force fields via energy landscape explorations in implicit solvent

Predicting the structure and dynamics of RNA molecules still proves challenging because of the relative scarcity of experimental RNA structures on which to train models and the very sensitive nature of RNA towards its environment. In the last decade, several atomistic force fields specifically designed for RNA have been proposed and are commonly used for simulations. However, it is not necessarily clear which force field is the most suitable for a given RNA molecule. In this contribution, we propose the use of the computational energy landscape framework to explore the energy landscape of RNA systems as it can bring complementary information to the more standard approaches of enhanced sampling simulations based on molecular dynamics. We apply the EL framework to the study of a small RNA pseudoknot, the Aquifex aeolicus tmRNA pseudoknot PK1, and we compare the results of five different RNA force fields currently available in the AMBER simulation software, in implicit solvent. With this computational approach, we can not only compare the predicted 'native' states for the different force fields, but the method enables us to study metastable states as well. As a result, our comparison not only looks at structural features of low energy folded structures, but provides insight into folding pathways and higher energy excited states, opening to the possibility of assessing the validity of force fields also based on kinetics and experiments providing information on metastable and unfolded states.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-024-01202-9.

Structural dynamics of therapeutic nucleic acids with phosphorothioate backbone modifications

Antisense oligonucleotides (ASOs) offer ground-breaking possibilities for selective pharmacological intervention for any gene product-related disease. Therapeutic ASOs contain extensive chemical modifications that improve stability to enzymatic cleavage and modulate binding affinity relative to natural RNA/DNA. Molecular dynamics (MD) simulation can provide valuable insights into how such modifications affect ASO conformational sampling and target binding. However, force field parameters for chemically modified nucleic acids (NAs) are still underdeveloped. To bridge this gap, we developed parameters to allow simulations of ASOs with the widely applied phosphorothioate (PS) backbone modification, and validated these in extensive all-atom MD simulations of relevant PS-modified NA systems representing B-DNA, RNA, and DNA/RNA hybrid duplex structures. Compared to the corresponding natural NAs, single PS substitutions had marginal effects on the ordered DNA/RNA duplex, whereas substantial effects of phosphorothioation were observed in single-stranded RNA and B-DNA, corroborated by the experimentally derived structure data. We find that PS-modified NAs shift between high and low twist states, which could affect target recognition and protein interactions for phosphorothioated oligonucleotides. Furthermore, conformational sampling was markedly altered in the PS-modified ssRNA system compared to that of the natural oligonucleotide, indicating sequence-dependent effects on conformational preference that may in turn influence duplex formation.

Dynamics of Protein-RNA Interfaces Using All-Atom Molecular Dynamics Simulations

Facing the current challenges posed by human health diseases requires the understanding of cell machinery at a molecular level. The interplay between proteins and RNA is key for any physiological phenomenon, as well protein-RNA interactions. To understand these interactions, many experimental techniques have been developed, spanning a very wide range of spatial and temporal resolutions. In particular, the knowledge of tridimensional structures of protein-RNA complexes provides structural, mechanical, and dynamical pieces of information essential to understand their functions. To get insights into the dynamics of protein-RNA complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on nine different protein-RNA complexes with different functions and interface size by taking into account the bound and unbound forms. First, we characterized structural changes upon binding and, for the RNA part, the change in the puckering. Second, we extensively analyzed the interfaces, their dynamics and structural properties, and the structural waters involved in the binding, as well as the contacts mediated by them. Based on our analysis, the interfaces rearranged during the simulation time showing alternative and stable residue-residue contacts with respect to the experimental structure.

Integrating cryo-OrbiSIMS with computational modelling and metadynamics simulations enhances RNA structure prediction at atomic resolution

The 3D architecture of RNAs governs their molecular interactions, chemical reactions, and biological functions. However, a large number of RNAs and their protein complexes remain poorly understood due to the limitations of conventional structural biology techniques in deciphering their complex structures and dynamic interactions. To address this limitation, we have benchmarked an integrated approach that combines cryogenic OrbiSIMS, a state-of-the-art solid-state mass spectrometry technique, with computational methods for modelling RNA structures at atomic resolution with enhanced precision. Furthermore, using 7SK RNP as a test case, we have successfully determined the full 3D structure of a native RNA in its apo, native and disease-remodelled states, which offers insights into the structural interactions and plasticity of the 7SK complex within these states. Overall, our study establishes cryo-OrbiSIMS as a valuable tool in the field of RNA structural biology as it enables the study of challenging, native RNA systems.

Near-Infrared Light Activated Formulation for the Spatially Controlled Release of CRISPR-Cas9 Ribonucleoprotein for Brain Gene Editing

The CRISPR/Cas9 system has emerged as a promising platform for gene editing; however, the lack of an efficient and safe delivery system to introduce it into cells continues to hinder clinical translation. Here, we report a rationally designed gene-editing nanoparticle (NP) formulation for brain applications: an sgRNA:Cas9 ribonucleoprotein complex is immobilized on the NP surface by oligonucleotides that are complementary to the sgRNA. Irradiation of the formulation with a near-infrared (NIR) laser generates heat in the NP, leading to the release of the ribonucleoprotein complex. The gene-editing potential of the formulation was demonstrated in vitro at the single-cell level. The safety and gene editing of the formulation were also demonstrated in the brains of reporter mice, specifically in the subventricular zone after intracerebral administration and in the olfactory bulb after intranasal administration. The formulation presented here offers a new strategy for the spatially controlled delivery of the CRISPR system to the brain.

The 8-17 DNAzyme can operate in a single active structure regardless of metal ion cofactor

DNAzymes - synthetic enzymes made of DNA - have long attracted attention as RNA-targeting therapeutic agents. Yet, as of now, no DNAzyme-based drug has been approved, partially due to our lacking understanding of their molecular mode of action. In this work we report the solution structure of 8-17 DNAzyme bound to a Zn2+ ion solved through NMR spectroscopy. Surprisingly, it turned out to be very similar to the previously solved Pb2+-bound form (catalytic domain RMSD = 1.28 Å), despite a long-standing literature consensus that Pb2+ recruits a different DNAzyme fold than other metal ion cofactors. Our follow-up NMR investigations in the presence of other ions - Mg2+, Na+, and Pb2+ - suggest that at DNAzyme concentrations used in NMR all these ions induce a similar tertiary fold. Based on these findings, we propose a model for 8-17 DNAzyme interactions with metal ions postulating the existence of only a single catalytically-active structure, yet populated to a different extent depending on the metal ion cofactor. Our results provide structural information on the 8-17 DNAzyme in presence of non-Pb2+ cofactors, including the biologically relevant Mg2+ ion.

Interaction of 3,9-disubstituted acridine with single stranded poly(rA), double stranded poly(rAU) and triple stranded poly(rUAU): molecular docking - A spectroscopic tandem study

RNA plays an important role in many biological processes which are crucial for cell survival, and it has been suggested that it may be possible to inhibit individual processes involved in many diseases by targeting specific sequences of RNA. The aim of this work is to determine the affinity of novel 3,9-disubstited acridine derivative 1 with three different RNA molecules, namely single stranded poly(rA), double stranded homopolymer poly(rAU) and triple stranded poly(rUAU). The results of the absorption titration assays show that the binding constant of the novel derivative to the RNA molecules was in the range of 1.7-6.2 × 104 mol dm-3. The fluorescence and circular dichroism titration assays revealed considerable changes. The most significant results in terms of interpreting the nature of the interactions were the melting temperatures of the RNA samples in complexes with the 1. In the case of poly(rA), denaturation resulted in a self-structure formation; increased stabilization was observed for poly(rAU), while the melting points of the ligand-poly(rUAU) complex showed significant destabilization as a result of the interaction. The principles of molecular mechanics were applied to propose the non-bonded interactions within the binding complex, pentariboadenylic acid and acridine ligand as the study model. Initial molecular docking provided the input structure for advanced simulation techniques. Molecular dynamics simulation and cluster analysis reveal π - π stacking and the hydrogen bonds formation as the main forces that can stabilize the binding complex. Subsequent MM-GBSA calculations showed negative binding enthalpy accompanied the complex formation and proposed the most preferred conformation of the interaction complex.

Origin of Enantioselectivity in Engineered Cytochrome c-Catalyzed Carbon-Radical FePP Hydrolysis Revealed Using QM/MM (ABEEM Polarizable Force Field) and MD Simulations

The origin of highly efficient asymmetric aminohydroxylation of styrene catalyzed by engineered cytochrome c is investigated by the developed Atom-Bond Electronegativity Equalization Method polarizable force field (ABEEM PFF), which is a combined outcome of electronic and steric effects. Model molecules were used to establish the charge parameters of the ABEEM PFF, for which the bond-stretching and angle-bending parameters were obtained by using a combination of modified Seminario and scan methods. The interactions between carbon-radical Fe-porphyrin (FePP) and waters are simulated by molecular dynamics, which shows a clear preference for the pre-R over the pre-S. This preference is attributed to the hydrogen-bond between the mutated 100S and 101P residues as well as van der Waals interactions, enforcing a specific conformation of the carbon-radical FePP complex within the binding pocket. Meanwhile, the hydrogen-bond between water and the nitrogen atom in the active intermediate dictates the stereochemical outcome. Quantum mechanics/molecular mechanics (QM/MM (ABEEM PFF)) and free-energy perturbation calculations elucidate that the 3RTS is characterized by sandwich-like structure among adjacent amino acid residues, which exhibits greater stability than crowed arrangement in 3STS and enables the R enantiomer to form more favorably. Thus, this study provides mechanistic insight into the catalytic reaction of hemoproteins.

Molecular Origin of Distinct Hydration Dynamics in Double Helical DNA and RNA Sequences

Water molecules are essential to determine the structure of nucleic acids and mediate their interactions with other biomolecules. Here, we characterize the hydration dynamics of analogous DNA and RNA double helices with unprecedented resolution and elucidate the molecular origin of their differences: first, the localization of the slowest hydration water molecules─in the minor groove in DNA, next to phosphates in RNA─and second, the markedly distinct hydration dynamics of the two phosphate oxygen atoms OR and OS in RNA. Using our Extended Jump Model for water reorientation, we assess the relative importance of previously proposed factors, including the local topography, water bridges, and the presence of ions. We show that the slow hydration dynamics at RNA OR sites is not due to bridging water molecules but is caused by both the larger excluded volume and the stronger initial H-bond next to OR, due to the different phosphate orientations in A-form double helical RNA.

Sequence- and Structure-Dependent Cytotoxicity of Phosphorothioate and 2'-O-Methyl Modified Single-Stranded Oligonucleotides

Single-stranded oligonucleotides (SSOs) are a rapidly expanding class of therapeutics that comprises antisense oligonucleotides, microRNAs, and aptamers, with ten clinically approved molecules. Chemical modifications such as the phosphorothioate backbone and the 2'-O-methyl ribose can improve the stability and pharmacokinetic properties of therapeutic SSOs, but they can also lead to toxicity in vitro and in vivo through nonspecific interactions with cellular proteins, gene expression changes, disturbed RNA processing, and changes in nuclear structures and protein distribution. In this study, we screened a mini library of 277 phosphorothioate and 2'-O-methyl-modified SSOs, with or without mRNA complementarity, for cytotoxic properties in two cancer cell lines. Using circular dichroism, nucleic magnetic resonance, and molecular dynamics simulations, we show that phosphorothioate- and 2'-O-methyl-modified SSOs that form stable hairpin structures through Watson-Crick base pairing are more likely to be cytotoxic than those that exist in an extended conformation. In addition, moderate and highly cytotoxic SSOs in our dataset have a higher mean purine composition than pyrimidine. Overall, our study demonstrates a structure-cytotoxicity relationship and indicates that the formation of stable hairpins should be a consideration when designing SSOs toward optimal therapeutic profiles.

Three-metal ion mechanism of cross-linked and uncross-linked DNA polymerase β: A theoretical study

In our recent publication, we have proposed a revised base excision repair pathway in which DNA polymerase β (Polβ) catalyzes Schiff base formation prior to the gap-filling DNA synthesis followed by β-elimination. In addition, the polymerase activity of Polβ employs the "three-metal ion mechanism" instead of the long-standing "two-metal ion mechanism" to catalyze phosphodiester bond formation based on the fact derived from time-resolved x-ray crystallography that a third Mg2+ was captured in the polymerase active site after the chemical reaction was initiated. In this study, we develop the models of the uncross-linked and cross-linked Polβ complexes and investigate the "three-metal ion mechanism" vs the "two-metal ion mechanism" by using the quantum mechanics/molecular mechanics molecular dynamics simulations. Our results suggest that the presence of the third Mg2+ ion stabilizes the reaction-state structures, strengthens correct nucleotide binding, and accelerates phosphodiester bond formation. The improved understanding of Polβ's catalytic mechanism provides valuable insights into DNA replication and damage repair.