Mg2+ ions: do they bind to nucleobase nitrogens?

Rapid DNA unwinding accelerates genome editing by engineered CRISPR-Cas9

Thermostable clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas9) enzymes could improve genome-editing efficiency and delivery due to extended protein lifetimes. However, initial experimentation demonstrated Geobacillus stearothermophilus Cas9 (GeoCas9) to be virtually inactive when used in cultured human cells. Laboratory-evolved variants of GeoCas9 overcome this natural limitation by acquiring mutations in the wedge (WED) domain that produce >100-fold-higher genome-editing levels. Cryoelectron microscopy (cryo-EM) structures of the wild-type and improved GeoCas9 (iGeoCas9) enzymes reveal extended contacts between the WED domain of iGeoCas9 and DNA substrates. Biochemical analysis shows that iGeoCas9 accelerates DNA unwinding to capture substrates under the magnesium-restricted conditions typical of mammalian but not bacterial cells. These findings enabled rational engineering of other Cas9 orthologs to enhance genome-editing levels, pointing to a general strategy for editing enzyme improvement. Together, these results uncover a new role for the Cas9 WED domain in DNA unwinding and demonstrate how accelerated target unwinding dramatically improves Cas9-induced genome-editing activity.

Geometric deep learning for the prediction of magnesium-binding sites in RNA structures

Magnesium ions (Mg2+) are essential for the folding, functional expression, and structural stability of RNA molecules. However, predicting Mg2+-binding sites in RNA molecules based solely on RNA structures is still challenging. The molecular surface, characterized by a continuous shape with geometric and chemical properties, is important for RNA modelling and carries essential information for understanding the interactions between RNAs and Mg2+ ions. Here, we propose an approach named RNA-magnesium ion surface interaction fingerprinting (RMSIF), a geometric deep learning-based conceptual framework to predict magnesium ion binding sites in RNA structures. To evaluate the performance of RMSIF, we systematically enumerated decoy Mg2+ ions across a full-space grid within the range of 2 to 10 Å from the RNA molecule and made predictions accordingly. Visualization techniques were used to validate the prediction results and calculate success rates. Comparative assessments against state-of-the-art methods like MetalionRNA, MgNet, and Metal3DRNA revealed that RMSIF achieved superior success rates and accuracy in predicting Mg2+-binding sites. Additionally, in terms of the spatial distribution of Mg2+ ions within the RNA structures, a majority were situated in the deep grooves, while a minority occupied the shallow grooves. Collectively, the conceptual framework developed in this study holds promise for advancing insights into drug design, RNA co-transcriptional folding, and structure prediction.

Mechanism of Fluoride Ion Encapsulation by Magnesium Ions in a Bacterial Riboswitch

Riboswitches sense various ions in bacteria and activate gene expression to synthesize proteins that help maintain ion homeostasis. The crystal structure of the aptamer domain (AD) of the fluoride riboswitch shows that the F- ion is encapsulated by three Mg2+ ions bound to the ligand-binding domain (LBD) located at the core of the AD. The assembly mechanism of this intricate structure is unknown. To this end, we performed computer simulations using coarse-grained and all-atom RNA models to bridge multiple time scales involved in riboswitch folding and ion binding. We show that F- encapsulation by the Mg2+ ions bound to the riboswitch involves multiple sequential steps. Broadly, two Mg2+ ions initially interact with the phosphate groups of the LBD using water-mediated outer-shell coordination and transition to a direct inner-shell interaction through dehydration to strengthen their interaction with the LBD. We propose that the efficient binding mode of the third Mg2+ and F- is that they form a water-mediated ion pair and bind to the LBD simultaneously to minimize the electrostatic repulsion between three Mg2+ bound to the LBD. The tertiary stacking interactions among the LBD nucleobases alone are insufficient to stabilize the alignment of the phosphate groups to facilitate Mg2+ binding. We show that the stability of the whole assembly is an intricate balance of the interactions among the five phosphate groups, three Mg2+, and the encapsulated F- ion aided by the Mg2+ solvated water. These insights are helpful in the rational design of RNA-based ion sensors and fast-switching logic gates.

Reanalysis of a μ opioid receptor crystal structure reveals a covalent adduct with BU72

BACKGROUND

The first crystal structure of the active μ opioid receptor (μOR) exhibited several unexplained features. The ligand BU72 exhibited many extreme deviations from ideal geometry, along with unexplained electron density. I previously showed that inverting the benzylic configuration resolved these problems, establishing revised stereochemistry of BU72 and its analog BU74. However, another problem remains unresolved: additional unexplained electron density contacts both BU72 and a histidine residue in the N-terminus, revealing the presence of an as-yet unidentified atom.

RESULTS

These short contacts and uninterrupted density are inconsistent with non-covalent interactions. Therefore, BU72 and μOR form a covalent adduct, rather than representing two separate entities as in the original model. A subsequently proposed magnesium complex is inconsistent with multiple lines of evidence. However, oxygen fits the unexplained density well. While the structure I propose is tentative, similar adducts have been reported previously in the presence of reactive oxygen species. Moreover, known sources of reactive oxygen species were present: HEPES buffer, nickel ions, and a sequence motif that forms redox-active nickel complexes. This motif contacts the unexplained density. The adduct exhibits severe strain, and the tethered N-terminus forms contacts with adjacent residues. These forces, along with the nanobody used as a G protein substitute, would be expected to influence the receptor conformation. Consistent with this, the intracellular end of the structure differs markedly from subsequent structures of active μOR bound to Gi protein.

CONCLUSIONS

Later Gi-bound structures are likely to be more accurate templates for ligand docking and modelling of active G protein-bound μOR. The possibility of reactions like this should be considered in the choice of protein truncation sites and purification conditions, and in the interpretation of excess or unexplained density.

Structure and function analysis of a type III preQ1-I riboswitch from Escherichia coli reveals direct metabolite sensing by the Shine-Dalgarno sequence

Riboswitches are small noncoding RNAs found primarily in the 5' leader regions of bacterial messenger RNAs where they regulate expression of downstream genes in response to binding one or more cellular metabolites. Such noncoding RNAs are often regulated at the translation level, which is thought to be mediated by the accessibility of the Shine-Dalgarno sequence (SDS) ribosome-binding site. Three classes (I-III) of prequeuosine1 (preQ1)-sensing riboswitches are known that control translation. Class I is divided into three subtypes (types I-III) that have diverse mechanisms of sensing preQ1, which is involved in queuosine biosynthesis. To provide insight into translation control, we determined a 2.30 Å-resolution cocrystal structure of a class I type III preQ1-sensing riboswitch identified in Escherichia coli (Eco) by bioinformatic searches. The Eco riboswitch structure differs from previous preQ1 riboswitch structures because it has the smallest naturally occurring aptamer and the SDS directly contacts the preQ1 metabolite. We validated structural observations using surface plasmon resonance and in vivo gene-expression assays, which showed strong switching in live E. coli. Our results demonstrate that the Eco riboswitch is relatively sensitive to mutations that disrupt noncanonical interactions that form the pseudoknot. In contrast to type II preQ1 riboswitches, a kinetic analysis showed that the type III Eco riboswitch strongly prefers preQ1 over the chemically similar metabolic precursor preQ0. Our results reveal the importance of noncanonical interactions in riboswitch-driven gene regulation and the versatility of the class I preQ1 riboswitch pseudoknot as a metabolite-sensing platform that supports SDS sequestration.

Structure of the actively translating plant 80S ribosome at 2.2 Å resolution

In plant cells, translation occurs in three compartments: the cytosol, the plastids and the mitochondria. While the structures of the (prokaryotic-type) ribosomes in plastids and mitochondria are well characterized, high-resolution structures of the eukaryotic 80S ribosomes in the cytosol have been lacking. Here the structure of translating tobacco (Nicotiana tabacum) 80S ribosomes was solved by cryo-electron microscopy with a global resolution of 2.2 Å. The ribosome structure includes two tRNAs, decoded mRNA and the nascent peptide chain, thus providing insights into the molecular underpinnings of the cytosolic translation process in plants. The map displays conserved and plant-specific rRNA modifications and the positions of numerous ionic cofactors, and it uncovers the role of monovalent ions in the decoding centre. The model of the plant 80S ribosome enables broad phylogenetic comparisons that reveal commonalities and differences in the ribosomes of plants and those of other eukaryotes, thus putting our knowledge about eukaryotic translation on a firmer footing.

Discovery of small molecules that target a tertiary-structured RNA

There is growing interest in therapeutic intervention that targets disease-relevant RNAs using small molecules. While there have been some successes in RNA-targeted small-molecule discovery, a deeper understanding of structure-activity relationships in pursuing these targets has remained elusive. One of the best-studied tertiary-structured RNAs is the theophylline aptamer, which binds theophylline with high affinity and selectivity. Although not a drug target, this aptamer has had many applications, especially pertaining to genetic control circuits. Heretofore, no compound has been shown to bind the theophylline aptamer with greater affinity than theophylline itself. However, by carrying out a high-throughput screen of low-molecular-weight compounds, several unique hits were identified that are chemically distinct from theophylline and bind with up to 340-fold greater affinity. Multiple atomic-resolution X-ray crystal structures were determined to investigate the binding mode of theophylline and four of the best hits. These structures reveal both the rigidity of the theophylline aptamer binding pocket and the opportunity for other ligands to bind more tightly in this pocket by forming additional hydrogen-bonding interactions. These results give encouragement that the same approaches to drug discovery that have been applied so successfully to proteins can also be applied to RNAs.

RNA as a Major-Groove Ligand: RNA-RNA and RNA-DNA Triplexes Formed by GAA and UUC or TTC Sequences

Friedreich's ataxia is associated with noncanonical nucleic acid structures that emerge when GAA:TTC repeats in the first intron of the FXN gene expand beyond a critical number of repeats. Specifically, the noncanonical repeats are associated with both triplexes and R-loops. Here, we present an in silico investigation of all possible triplexes that form by attaching a third RNA strand to an RNA:RNA or DNA:DNA duplex, complementing previous DNA-based triplex studies. For both new triplexes results are similar. For a pyridimine UUC⁺ third strand, the parallel orientation is stable while its antiparallel counterpart is unstable. For a neutral GAA third strand, the parallel conformation is stable. A protonated GA⁺A third strand is stable in both parallel and antiparallel orientations. We have also investigated Na⁺ and Mg²⁺ ion distributions around the triplexes. The presence of Mg²⁺ ions helps stabilize neutral, antiparallel GAA triplexes. These results (along with previous DNA-based studies) allow for the emergence of a complete picture of the stability and structural characteristics of triplexes based on the GAA and TTC/UUC sequences, thereby contributing to the field of trinucleotide repeats and the associated unusual structures that trigger expansion.

Structure and thermodynamics of a UGG motif interacting with Ba2+ and other metal ions: accommodating changes in the RNA structure and the presence of a G(syn)-G(syn) pair

The self-complementary triplet 5'UGG3'/5'UGG3' is a particular structural motif containing noncanonical G-G pair and two U·G wobble pairs. It constitutes a specific structural and electrostatic environment attracting metal ions, particularly Ba2+ ions. Crystallographic research has shown that two Ba2+ cations are located in the major groove of the helix and interact directly with the UGG triplet. A comparison with the unliganded structure has revealed global changes in the RNA structure in the presence of metal ions, whereas thermodynamic measurements have shown increased stability. Moreover, in the structure with Ba2+, an unusual noncanonical G(syn)-G(syn) pair is observed instead of the common G(anti)-G(syn). We further elucidate the metal binding properties of the UGG/UGG triplet by performing crystallographic and thermodynamic studies using DSC and UV melting with other metal ions. The results explain the preferences of the UGG sequence for Ba2+ cations and point to possible applications of this metal-binding propensity.

An electrochemical aptasensor based on target triggered multiple-channel DNAzymes cycling amplification strategy with PtFe@Co-MOF as signal amplifier

A novel electrochemical aptasensor for the detection of Aflatoxin B1 (AFB1) was developed for the first time by using the target-triggered multiple-channel deoxyribozymes (DNAzymes) cycling amplified assay with Pt Fe doped NH₂-Co-MOF (PtFe@Co-MOF) as a signal amplifier. In the presence of AFB1, a self-assembling cross-over nucleic structure could be triggered by AFB1 via two aptamers' structure switching for strand displacement, resulting in four channels of Mg²⁺-dependent DNAzyme recycling simultaneously to multiply the detection signals. These DNAzymes cyclically split the substrate sequence to release the PtFe@Co-MOF labeled detection probe (DP), which is subsequently hybridized with the capture probes on the Au-deposited glassy carbon electrode. The fabrication procedure was characterized by differential pulse voltammetry, and the results of the morphological and element composition characteristics methods were analyzed to determine the successful preparation of PtFe@Co-MOF. The limit of detection (LOD) for AFB1 detection was 2 pg mL^-1 with a linear range from 5 pg mL^-1 to 80 ng mL^-1. By comparison, the enhanced detection sensitivity has been found to originate from the efficient shearing of DNAzymes, enhanced peroxidase-like capability, and multiple active sites of PtFe@Co-MOF. Besides, this aptasensor showed high specificity for AFB1 compared with similar mycotoxins and exhibited high accuracy with low experimental cost and easy operation. Furthermore, the unique design of electrochemical aptasensors could provide a promising platform for the onsite determination of AFB1, as well as other targets by replacing the aptamer and other core recognition sequences.

DNA Antiadhesive Layer for Reusable Plasmonic Sensors: Nanostructure Pitch Effect

A long-term reusable sensor that provides the opportunity to easily regenerate the active surface and minimize the occurrence of undesired absorption events is an appealing solution that helps to cut down the costs and improve the device performances. Impressive advances have been made in the past years concerning the development of novel cutting-edge sensors, but the reusability can currently represent a challenge. Direct shielding of the sensor surface is not always applicable, because it can impact the device performance. This study reports an antiadhesive layer (AAL) made of 90 mg/mL DNA sodium salt from salmon testes (ssstDNA) for passivating gold plasmonic sensor surfaces. Our gold two-dimensional (2D) nanostructured plasmonic metasurfaces modified with AAL were used for DNA quantification. AAL is thin enough that the plasmonic sensor remains sensitive to subsequent deposition of DNA, which serves as an analyte. AAL protects the gold surface from unwanted nonspecific adsorption by enabling wash-off of the deposited analyte after analysis and thus recovery of the LSPR peak position (rLSPR). The calibration curve obtained on a single nanostructure (Achiral Octupolar, 100 nm pitch) gave an LOD = 105 ng/mL and an extraordinary dynamic range, performances comparable or superior to those of commercial UV-vis spectrometers for acid nucleic dosage. Two different analytes were tested: ssstDNA (∼2000 bp) in deionized water and double-strand DNA (dsDNA) of 546-1614 bp in 100 mM Tris buffer and 10 mM MgCl₂. The two nanostructures (Achiral Octupolar 25 and 100) were found to have the same sensitivity to DNA in deionized water but different sensitivity to DNA in a salt/buffer solution, opening a potential for solute discrimination. To the best of our knowledge, this is the first report on the use of AAL made of several kilobase-pairs-long dsDNA to produce a reusable plasmonic sensor. The working principle and limitations are drawn based on the LSPR and SERS study.

Structural Basis for Fluorescence Activation by Pepper RNA

Pepper is a fluorogenic RNA aptamer tag that binds to a variety of benzylidene-cyanophenyl (HBC) derivatives with tight affinity and activates their fluorescence. To investigate how Pepper RNA folds to create a binding site for HBC, we used antibody-assisted crystallography to determine the structures of Pepper bound to HBC530 and HBC599 to 2.3 and 2.7 Å resolutions, respectively. The structural data show that Pepper folds into an elongated structure and organizes nucleotides within an internal bulge to create the ligand binding site, assisted by an out-of-plane platform created by tertiary interactions with an adjacent bulge. As predicted from a lack of K⁺ dependence, Pepper does not use a G-quadruplex to form a binding pocket for HBC. Instead, Pepper uses a unique base-quadruple·base-triple stack to sandwich the ligand with a U·G wobble pair. Site-bound Mg²⁺ ions support ligand binding structurally and energetically. This research provides insight into the structural features that allow the Pepper aptamer to bind HBC and show how Pepper's function may expand to allow the in vivo detection of other small molecules and metals.

Lone Pair…π Contacts and Structure Signatures of r(UNCG) Tetraloops, Z-Turns, and Z-Steps: A WebFR3D Survey

Z-DNA and Z-RNA have long appeared as oddities to nucleic acid scientists. However, their Z-step constituents are recurrently observed in all types of nucleic acid systems including ribosomes. Z-steps are NpN steps that are isostructural to Z-DNA CpG steps. Among their structural features, Z-steps are characterized by the presence of a lone pair…π contact that involves the stacking of the ribose O4′ atom of the first nucleotide with the 3′-face of the second nucleotide. Recently, it has been documented that the CpG step of the ubiquitous r(UNCG) tetraloops is a Z-step. Accordingly, such r(UNCG) conformations were called Z-turns. It has also been recognized that an r(GAAA) tetraloop in appropriate conditions can shapeshift to an unusual Z-turn conformation embedding an ApA Z-step. In this report, we explore the multiplicity of RNA motifs based on Z-steps by using the WebFR3D tool to which we added functionalities to be able to retrieve motifs containing lone pair…π contacts. Many examples that underscore the diversity and universality of these motifs are provided as well as tutorial guidance on using WebFR3D. In addition, this study provides an extensive survey of crystallographic, cryo-EM, NMR, and molecular dynamics studies on r(UNCG) tetraloops with a critical view on how to conduct database searches and exploit their results.

Magnesium Force Fields for OPC Water with Accurate Solvation, Ion-Binding, and Water-Exchange Properties: Successful Transfer from SPC/E

Magnesium plays a vital role in a large variety of biological processes. To model such processes by molecular dynamics simulations, researchers rely on accurate force field parameters for Mg2+ and water. OPC is one of the most promising water models yielding an improved description of biomolecules in water. The aim of this work is to provide force field parameters for Mg2+ that lead to accurate simulation results in combination with OPC water. Using twelve different Mg2+ parameter sets, that were previously optimized with different water models, we systematically assess the transferability to OPC based on a large variety of experimental properties. The results show that the Mg2+ parameters for SPC/E are transferable to OPC and closely reproduce the experimental solvation free energy, radius of the first hydration shell, coordination number, activity derivative, and binding affinity toward the phosphate oxygens on RNA. Two optimal parameter sets are presented: MicroMg yields water exchange in OPC on the microsecond timescale in agreement with experiments. NanoMg yields accelerated exchange on the nanosecond timescale and facilitates the direct observation of ion binding events for enhanced sampling purposes.

Artificial Intelligence Resolves Kinetic Pathways of Magnesium Binding to RNA

Magnesium is an indispensable cofactor in countless vital processes. In order to understand its functional role, the characterization of the binding pathways to biomolecules such as RNA is crucial. Despite the importance, a molecular description is still lacking since the transition from the water-mediated outer-sphere to the direct inner-sphere coordination is on the millisecond time scale and therefore out of reach for conventional simulation techniques. To fill this gap, we use transition path sampling to resolve the binding pathways and to elucidate the role of the solvent in the binding process. The results reveal that the molecular void provoked by the leaving phosphate oxygen of the RNA is immediately filled by an entering water molecule. In addition, water molecules from the first and second hydration shell couple to the concerted exchange. To capture the intimate solute-solvent coupling, we perform a committor analysis as the basis for a machine learning algorithm that derives the optimal deep learning model from thousands of scanned architectures using hyperparameter tuning. The results reveal that the properly optimized deep network architecture recognizes the important solvent structures, extracts the relevant information, and predicts the commitment probability with high accuracy. Our results provide detailed insights into the solute-solvent coupling which is ubiquitous for kosmotropic ions and governs a large variety of biochemical reactions in aqueous solutions.

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D

Magnesium is essential in many vital processes. To correctly describe Mg2+ in physiological processes by molecular dynamics simulations, accurate force fields are fundamental. Despite the importance, force fields based on the commonly used 12-6 Lennard-Jones potential showed significant shortcomings. Recently progress was made by an optimization procedure that implicitly accounts for polarizability. The resulting microMg and nanoMg force fields (J. Chem. Theory Comput. 2021, 17, 2530-2540) accurately reproduce a broad range of experimental solution properties and the binding affinity to nucleic acids in TIP3P water. Since countless simulation studies rely on available water models and ion force fields, we here extend the optimization and provide Mg2+ parameters in combination with the SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D water models. For each water model, the Mg2+ force fields reproduce the solvation free energy, the distance to oxygens in the first hydration shell, the hydration number, the activity coefficient derivative in MgCl2 solutions, and the binding affinity and distance to the phosphate oxygens on nucleic acids. We present two parameter sets: MicroMg yields water exchange on the microsecond time scale and matches the experimental exchange rate. Depending on the water model, nanoMg yields accelerated water exchange in the range of 106 to 108 exchanges per second. The nanoMg parameters can be used to enhance the sampling of binding events, to obtain converged distributions of Mg2+, or to predict ion binding sites in biomolecular simulations. The parameter files are freely available at https://github.com/bio-phys/optimizedMgFFs.

Twisting and swiveling domain motions in Cas9 to recognize target DNA duplexes, make double-strand breaks, and release cleaved duplexes

The CRISPR-associated protein 9 (Cas9) has been engineered as a precise gene editing tool to make double-strand breaks. CRISPR-associated protein 9 binds the folded guide RNA (gRNA) that serves as a binding scaffold to guide it to the target DNA duplex via a RecA-like strand-displacement mechanism but without ATP binding or hydrolysis. The target search begins with the protospacer adjacent motif or PAM-interacting domain, recognizing it at the major groove of the duplex and melting its downstream duplex where an RNA-DNA heteroduplex is formed at nanomolar affinity. The rate-limiting step is the formation of an R-loop structure where the HNH domain inserts between the target heteroduplex and the displaced non-target DNA strand. Once the R-loop structure is formed, the non-target strand is rapidly cleaved by RuvC and ejected from the active site. This event is immediately followed by cleavage of the target DNA strand by the HNH domain and product release. Within CRISPR-associated protein 9, the HNH domain is inserted into the RuvC domain near the RuvC active site via two linker loops that provide allosteric communication between the two active sites. Due to the high flexibility of these loops and active sites, biophysical techniques have been instrumental in characterizing the dynamics and mechanism of the CRISPR-associated protein 9 nucleases, aiding structural studies in the visualization of the complete active sites and relevant linker structures. Here, we review biochemical, structural, and biophysical studies on the underlying mechanism with emphasis on how CRISPR-associated protein 9 selects the target DNA duplex and rejects non-target sequences.

The mutual interactions of RNA, counterions and water - quantifying the electrostatics at the phosphate-water interface

The structure and dynamics of polyanionic biomolecules, like RNA, are decisively determined by their electric interactions with the water molecules and the counterions in the environment. The solvation dynamics of the biomolecules involves a subtle balance of non-covalent and many-body interactions with structural fluctuations due to thermal motion occurring in a femto- to subnanosecond time range. This complex fluctuating many particle scenario is crucial in defining the properties of biological interfaces with far reaching significance for the folding of RNA structures and for facilitating RNA-protein interactions. Given the inherent complexity, suited model systems, carefully calibrated and benchmarked by experiments, are required to quantify the relevant interactions of RNA with the aqueous environment. In this feature article we summarize our recent progress in the understanding of the electrostatics at the biological interface of double stranded RNA (dsRNA) and transfer RNA (tRNA). Dimethyl phosphate (DMP) is introduced as a viable and rigorously accessible model system allowing the interaction strength with water molecules and counterions, their relevant fluctuation timescales and the spatial reach of interactions to be established. We find strong (up to ≈90 MV cm-1) interfacial electric fields with fluctuations extending up to ≈20 THz and demonstrate how the asymmetric stretching vibration νAS(PO2)- of the polarizable phosphate group can serve as the most sensitive probe for interfacial interactions, establishing a rigorous link between simulations and experiment. The approach allows for the direct interfacial observation of interactions of biologically relevant Mg2+ counterions with phosphate groups in contact pair geometries via the rise of a new absorption band imposed by exchange repulsion interactions at short interatomic distances. The systematic extension to RNA provides microscopic insights into the changes of the hydration structure that accompany the temperature induced melting of the dsRNA double helix and quantify the ionic interactions in the folded tRNA. The results show that pairs of negatively charged phosphate groups and Mg2+ ions represent a key structural feature of RNA embedded in water. They highlight the importance of binding motifs made of contact pairs in the electrostatic stabilization of RNA structures that have a strong impact on the surface potential and enable the fine tuning of the local electrostatic properties which are expected to be relevant for mediating the interactions between biomolecules.

Unique anticodon loop conformation with the flipped-out wobble nucleotide in the crystal structure of unbound tRNAVal

During protein synthesis on ribosome, tRNA recognizes its cognate codon of mRNA through base-pairing with the anticodon. The 5'-end nucleotide of the anticodon is capable of wobble base-pairing, offering a molecular basis for codon degeneracy. The wobble nucleotide is often targeted for post-transcriptional modification, which affects the specificity and fidelity of the decoding process. Flipping-out of a wobble nucleotide in the anticodon loop has been proposed to be necessary for modifying enzymes to access the target nucleotide, which has been captured in selective structures of protein-bound complexes. Meanwhile, all other structures of free or ribosome-bound tRNA display anticodon bases arranged in stacked conformation. We report the X-ray crystal structure of unbound tRNA^Val1 to a 2.04 Å resolution showing two different conformational states of wobble uridine in the anticodon loop, one stacked on the neighboring base and the other swiveled out toward solvent. In addition, the structure reveals a rare magnesium ion coordination to the nitrogen atom of a nucleobase, which has been sampled very rarely among known structures of nucleic acids.

Theoretical Studies of the Self Cleavage Pistol Ribozyme Mechanism

Ribozymes are huge complex biological catalysts composed of a combination of RNA and proteins. Nevertheless, there is a reduced number of small ribozymes, the self-cleavage ribozymes, that are formed just by RNA and, apparently, they existed in cells of primitive biological systems. Unveiling the details of these “fossils” enzymes can contribute not only to the understanding of the origins of life but also to the development of new simplified artificial enzymes. A computational study of the reactivity of the pistol ribozyme carried out by means of classical MD simulations and QM/MM hybrid calculations is herein presented to clarify its catalytic mechanism. Analysis of the geometries along independent MD simulations with different protonation states of the active site basic species reveals that only the canonical system, with no additional protonation changes, renders reactive conformations. A change in the coordination sphere of the Mg2+ ion has been observed during the simulations, which allows proposing a mechanism to explain the unique mode of action of the pistol ribozyme by comparison with other ribozymes. The present results are at the center of the debate originated from recent experimental and theoretical studies on pistol ribozyme.

Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations

Magnesium and calcium play an essential role in the folding and function of nucleic acids. To correctly describe their interactions with DNA and RNA in biomolecular simulations, an accurate parameterization is crucial. In most cases, the ion parameters are optimized based on a set of experimental solution properties such as solvation free energies, radial distribution functions, water exchange rates, and activity coefficient derivatives. However, the transferability of such bulk-optimized ion parameters to quantitatively describe biomolecular systems is limited. Here, we extend the applicability of our previous bulk-optimized parameters by including experimental binding affinities toward the phosphate oxygen on nucleic acids. In particular, we systematically adjust the combination rules that are an integral part of the pairwise interaction potentials of classical force fields. This allows us to quantitatively describe specific ion binding to nucleic acids without changing the solution properties in the most simple and efficient way. We show the advancement of the optimized Lorentz combination rule for two representative nucleic acid systems. For double-stranded DNA, the optimized combination rule for Ca²⁺ significantly improves the agreement with experiments, while the standard combination rule leads to unrealistically distorted DNA structures. For the add A-riboswitch, the optimized combination rule for Mg²⁺ improves the structure of two specifically bound Mg²⁺ ions as judged by the experimental distance to the binding site. Including experimental binding affinities toward specific ion binding sites on biomolecules, therefore, provides a promising perspective to develop a more accurate description of metal cations for biomolecular simulations.

Diverse relationships between metal ions and the ribosome

The ribosome requires metal ions for structural stability and translational activity. These metal ions are important for stabilizing the secondary structure of ribosomal RNA, binding of ribosomal proteins to the ribosome, and for interaction of ribosomal subunits. In this review, various relationships between ribosomes and metal ions, especially Mg2+ and Zn2+, are presented. Mg2+ regulates gene expression by modulating the translational stability and synthesis of ribosomes, which in turn contribute to the cellular homeostasis of Mg2+. In addition, Mg2+ can partly complement the function of ribosomal proteins. Conversely, a reduction in the cellular concentration of Zn2+ induces replacement of ribosomal proteins, which mobilizes free-Zn2+ in the cell and represses translation activity. Evolutional relationships between these metal ions and the ribosome are also discussed.

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties

Magnesium ions play an essential role in many vital processes. To correctly describe their interactions in molecular dynamics simulations, an accurate parametrization is crucial. Despite the importance and considerable scientific effort, current force fields based on the commonly used 12-6 Lennard-Jones interaction potential fail to reproduce a variety of experimental solution properties. In particular, no parametrization exists so far that simultaneously reproduces the solvation free energy and the distance to the water oxygens in the first hydration shell. Moreover, current Mg2+ force fields significantly underestimate the rate of water exchange leading to unrealistically slow exchange kinetics. In order to make progress in the development of improved models, we systematically optimize the Mg2+ parameters in combination with the TIP3P water model in a much larger parameter space than previously done. The results show that a long-ranged interaction potential and modified Lorentz-Berthelot combination rules allow us to accurately reproduce multiple experimental properties including the solvation free energy, the distances to the oxygens of the first hydration shell, the hydration number, the activity coefficient derivative in MgCl2 solutions, the self-diffusion coefficient, and the binding affinity to the phosphate oxygen of RNA. Matching this broad range of thermodynamic properties, we present two sets of optimal parameters: MicroMg yields water exchange on the microsecond timescale in agreement with experiments. NanoMg yields water exchange on the nanosecond timescale facilitating the direct observation of ion-binding events. As shown for the example of the add A-riboswitch, the optimized parameters correctly reproduce the structure of specifically bound ions and permit the de novo prediction of Mg2+-binding sites in biomolecular simulations.

Different tertiary interactions create the same important 3D features in a distinct flavivirus xrRNA

During infection by a flavivirus (FV), cells accumulate noncoding subgenomic flavivirus RNAs (sfRNAs) that interfere with several antiviral pathways. These sfRNAs are formed by structured RNA elements in the 3' untranslated region (UTR) of the viral genomic RNA, which block the progression of host cell exoribonucleases that have targeted the viral RNA. Previous work on these exoribonuclease-resistant RNAs (xrRNAs) from mosquito-borne FVs revealed a specific three-dimensional fold with a unique topology in which a ring-like structure protectively encircles the 5' end of the xrRNA. Conserved nucleotides make specific tertiary interactions that support this fold. Examination of more divergent FVs reveals differences in their 3' UTR sequences, raising the question of whether they contain xrRNAs and if so, how they fold. To answer this, we demonstrated the presence of an authentic xrRNA in the 3' UTR of the Tamana bat virus (TABV) and solved its structure by X-ray crystallography. The structure reveals conserved features from previously characterized xrRNAs, but in the TABV version these features are created through a novel set of tertiary interactions not previously seen in xrRNAs. This includes two important A-C interactions, four distinct backbone kinks, several ordered Mg²⁺ ions, and a C⁺-G-C base triple. The discovery that the same overall architecture can be achieved by very different sequences and interactions in distantly related flaviviruses provides insight into the diversity of this type of RNA and will inform searches for undiscovered xrRNAs in viruses and beyond.

Deflating the RNA Mg2+ bubble. Stereochemistry to the rescue!

Proper evaluation of the ionic structure of biomolecular systems through X ray and cryo-EM techniques remains challenging but is essential for advancing our understanding of the underlying structure/activity/solvent relationships. However, numerous studies overestimate the number of Mg2+ in deposited structures due to assignment errors finding their origin in improper consideration of stereochemical rules. Herein, to tackle such issues, we re-evaluate the PDBid 6QNR and 6SJ6 models of the ribosome ionic structure. We establish that stereochemical principles need to be carefully pondered when evaluating ion binding features, even when K+ anomalous signals are available as it is the case for the 6QNR PDB entry. For ribosomes, assignment errors can result in misleading conceptions of their solvent structure. For instance, present stereochemical analysis result in a significant decrease of the number of assigned Mg2+ in 6QNR, suggesting that K+ and not Mg2+ is the prevalent ion in the ribosome 1st solvation shell. We stress that the use of proper stereochemical guidelines in combination or not with other identification techniques, such as those pertaining to the detection of transition metals, of some anions and of K+ anomalous signals, is critical for deflating the current Mg2+ bubble witnessed in many ribosome and other RNA structures. We also stress that for the identification of lighter ions such as Mg2+, Na+, …, for which no anomalous signals can be detected, stereochemistry coupled with high resolution structures (<2.4 Å) remain the best currently available option.

Catalytic Mechanism of Non-Target DNA Cleavage in CRISPR-Cas9 Revealed by Ab Initio Molecular Dynamics

CRISPR-Cas9 is a cutting-edge genome editing technology, which uses the endonuclease Cas9 to introduce mutations at desired sites of the genome. This revolutionary tool is promising to treat a myriad of human genetic diseases. Nevertheless, the molecular basis of DNA cleavage, which is a fundamental step for genome editing, has not been established. Here, quantum-classical molecular dynamics (MD) and free energy methods are used to disclose the two-metal-dependent mechanism of phosphodiester bond cleavage in CRISPR-Cas9. Ab initio MD reveals a conformational rearrangement of the Mg2+-bound RuvC active site, which entails the relocation of H983 to act as a general base. Then, the DNA cleavage proceeds through a concerted associative pathway fundamentally assisted by the joint dynamics of the two Mg2+ ions. This clarifies previous controversial experimental evidence, which could not fully establish the catalytic role of the conserved H983 and the metal cluster conformation. The comparison with other two-metal-dependent enzymes supports the identified mechanism and suggests a common catalytic strategy for genome editing and recombination. Overall, the non-target DNA cleavage catalysis described here resolves a fundamental open question in the CRISPR-Cas9 biology and provides valuable insights for improving the catalytic efficiency and the metal-dependent function of the Cas9 enzyme, which are at the basis of the development of genome editing tools.

Structural Analysis of (p)ppGpp Reveals Its Versatile Binding Pattern for Diverse Types of Target Proteins

(p)ppGpp, highly phosphorylated guanosine, are global regulatory nucleotides that modulate several biochemical events in bacterial physiology ranging from core central dogma to various metabolic pathways. Conventionally, (p)ppGpp collectively refers to two nucleotides, ppGpp, and pppGpp in the literature. Initially, (p)ppGpp has been discovered as a transcription regulatory molecule as it binds to RNA polymerase and regulates transcriptional gene regulation. During the past decade, several other target proteins of (p)ppGpp have been discovered and as of now, more than 30 proteins have been reported to be regulated by the binding of these two signaling nucleotides. The regulation of diverse biochemical activities by (p)ppGpp requires fine-tuned molecular interactions with various classes of proteins so that it can moderate varied functions. Here we report a structural dynamics of (p)ppGpp in the unbound state using well-defined computational tools and its interactions with target proteins to understand the differential regulation by (p)ppGpp at the molecular level. We carried out replica exchange molecular dynamics simulation studies to enhance sampling of conformations during (p)ppGpp simulation. The detailed comparative analysis of torsion angle conformation of ribose sugar of unbound (p)ppGpp and bound states of (p)ppGpp was carried out. The structural dynamics shows that two linear phosphate chains provide plasticity to (p)ppGpp nucleotides for the binding to diverse proteins. Moreover, the intermolecular interactions between (p)ppGpp and target proteins were characterized through various physicochemical parameters including, hydrogen bonds, van der Waal’s interactions, aromatic stacking, and side chains of interacting residues of proteins. Surprisingly, we observed that interactions of (p)ppGpp to target protein have a consensus binding pattern for a particular functional class of enzymes. For example, the binding of (p)ppGpp to RNA polymerase is significantly different from the binding of (p)ppGpp to the proteins involved in the ribosome biogenesis pathway. Whereas, (p)ppGpp binding to enzymes involved in nucleotide metabolism facilitates the functional regulation through oligomerization. Analysis of these datasets revealed that guanine base-specific contacts are key determinants to discriminate functional class of protein. Altogether, our studies provide significant information to understand the differential interaction pattern of (p)ppGpp to its target and this information may be useful to design antibacterial compounds based on (p)ppGpp analogs.

Analysis of a preQ1-I riboswitch in effector-free and bound states reveals a metabolite-programmed nucleobase-stacking spine that controls gene regulation

Riboswitches are structured RNA motifs that recognize metabolites to alter the conformations of downstream sequences, leading to gene regulation. To investigate this molecular framework, we determined crystal structures of a preQ1-I riboswitch in effector-free and bound states at 2.00 Å and 2.65 Å-resolution. Both pseudoknots exhibited the elusive L2 loop, which displayed distinct conformations. Conversely, the Shine-Dalgarno sequence (SDS) in the S2 helix of each structure remained unbroken. The expectation that the effector-free state should expose the SDS prompted us to conduct solution experiments to delineate environmental changes to specific nucleobases in response to preQ1. We then used nudged elastic band computational methods to derive conformational-change pathways linking the crystallographically-determined effector-free and bound-state structures. Pathways featured: (i) unstacking and unpairing of L2 and S2 nucleobases without preQ1-exposing the SDS for translation and (ii) stacking and pairing L2 and S2 nucleobases with preQ1-sequestering the SDS. Our results reveal how preQ1 binding reorganizes L2 into a nucleobase-stacking spine that sequesters the SDS, linking effector recognition to biological function. The generality of stacking spines as conduits for effector-dependent, interdomain communication is discussed in light of their existence in adenine riboswitches, as well as the turnip yellow mosaic virus ribosome sensor.

Mg2+-Dependent Methyl Transfer by a Knotted Protein: A Molecular Dynamics Simulation and Quantum Mechanics Study

Mg²⁺ is required for the catalytic activity of TrmD, a bacteria-specific methyltransferase that is made up of a protein topological knot-fold, to synthesize methylated m¹G37-tRNA to support life. However, neither the location of Mg²⁺ in the structure of TrmD nor its role in the catalytic mechanism is known. Using molecular dynamics (MD) simulations, we identify a plausible Mg²⁺ binding pocket within the active site of the enzyme, wherein the ion is coordinated by two aspartates and a glutamate. In this position, Mg²⁺ additionally interacts with the carboxylate of a methyl donor cofactor S-adenosylmethionine (SAM). The computational results are validated by experimental mutation studies, which demonstrate the importance of the Mg²⁺-binding residues for the catalytic activity. The presence of Mg²⁺ in the binding pocket induces SAM to adopt a unique bent shape required for the methyl transfer activity and causes a structural reorganization of the active site. Quantum mechanical calculations show that the methyl transfer is energetically feasible only when Mg²⁺ is bound in the position revealed by the MD simulations, demonstrating that its function is to align the active site residues within the topological knot-fold in a geometry optimal for catalysis. The obtained insights provide the opportunity for developing a strategy of antibacterial drug discovery based on targeting of Mg²⁺-binding to TrmD.

VeriNA3d: an R package for nucleic acids data mining

SUMMARY

veriNA3d is an R package for the analysis of nucleic acids structural data, with an emphasis in complex RNA structures. In addition to single-structure analyses, veriNA3d also implements functions to handle whole datasets of mmCIF/PDB structures that could be retrieved from public/local repositories. Our package aims to fill a gap in the data mining of nucleic acids structures to produce flexible and high throughput analysis of structural databases.

AVAILABILITY AND IMPLEMENTATION

http://mmb.irbbarcelona.org/gitlab/dgallego/veriNA3d.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Hofmeister Series for Metal-Cation-RNA Interactions: The Interplay of Binding Affinity and Exchange Kinetics

A large variety of physicochemical properties involving RNA depends on the type of metal cation present in solution. In order to gain microscopic insight into the origin of these ion specific effects, we apply molecular dynamics simulations to describe the interactions of metal cations and RNA. For the three most common ion binding sites on RNA, we calculate the binding affinities and exchange rates of eight different mono- and divalent metal cations. Our results reveal that binding sites involving phosphate groups preferentially bind metal cations with high charge density (such as Mg²⁺) in inner-sphere conformations while binding sites involving N7 or O6 atoms preferentially bind cations with low charge density (such as K⁺). The binding affinity therefore follows a direct Hofmeister series at the backbone but is reversed at the nucleobases leading to a high selectivity of ion binding sites on RNA. In addition, the exchange rates for cation binding cover almost 5 orders of magnitude, leading to a vastly different time scale for the lifetimes of contact pairs. Taken together, the site-specific binding affinities and the specific lifetime of contact pairs provide the microscopic explanation of ion specific effects observed in a wide variety of macroscopic RNA properties. Finally, combining the results from atomistic simulations with extended Poisson-Boltzmann theory allows us to predict the distribution of metal cations around double-stranded RNA at finite concentrations and to reproduce the results of ion counting experiments with good accuracy.

Charge Density of Cation Determines Inner versus Outer Shell Coordination to Phosphate in RNA

Divalent cations are often required to fold RNA, which is a highly charged polyanion. Condensation of ions, such as Mg²⁺ or Ca²⁺, in the vicinity of RNA renormalizes the effective charges on the phosphate groups, thus minimizing the intra RNA electrostatic repulsion. The prevailing view is that divalent ions bind diffusively in a nonspecific manner. In sharp contrast, we arrive at the exact opposite conclusion using a theory for the interaction of ions with the phosphate groups using RISM theory in conjunction with simulations based on an accurate three-interaction-site RNA model. The divalent ions bind in a nucleotide-specific manner using either the inner (partially dehydrated) or outer (fully hydrated) shell coordination. The high charge density Mg²⁺ ion has a preference to bind to the outer shell, whereas the opposite is the case for Ca²⁺. Surprisingly, we find that bridging interactions, involving ions that are coordinated to two or more phosphate groups, play a crucial role in maintaining the integrity of the folded state. Their importance could become increasingly prominent as the size of the RNA increases. Because the modes of interaction of divalent ions with DNA are likely to be similar, we propose that specific inner and outer shell coordination could play a role in DNA condensation, and perhaps genome organization as well.

The L-platform/L-scaffold framework: a blueprint for RNA-cleaving nucleic acid enzyme design

We develop an L-platform/L-scaffold framework we hypothesize may serve as a blueprint to facilitate site-specific RNA-cleaving nucleic acid enzyme design. Building on the L-platform motif originally described by Suslov and coworkers, we identify new critical scaffolding elements required to anchor a conserved general base guanine ("L-anchor") and bind functionally important metal ions at the active site ("L-pocket"). Molecular simulations, together with a broad range of experimental structural and functional data, connect the L-platform/L-scaffold elements to necessary and sufficient conditions for catalytic activity. We demonstrate that the L-platform/L-scaffold framework is common to five of the nine currently known naturally occurring ribozyme classes (Twr, HPr, VSr, HHr, Psr), and intriguingly from a design perspective, the framework also appears in an artificially engineered DNAzyme (8-17dz). The flexibility of the L-platform/L-scaffold framework is illustrated on these systems, highlighting modularity and trends in the variety of known general acid moieties that are supported. These trends give rise to two distinct catalytic paradigms, building on the classifications proposed by Wilson and coworkers and named for the implicated general base and acid. The "G + A" paradigm (Twr, HPr, VSr) exclusively utilizes nucleobase residues for chemistry, and the "G + M + " paradigm (HHr, 8-17dz, Psr) involves structuring of the "L-pocket" metal ion binding site for recruitment of a divalent metal ion that plays an active role in the chemical steps of the reaction. Finally, the modularity of the L-platform/L-scaffold framework is illustrated in the VS ribozyme where the "L-pocket" assumes the functional role of the "L-anchor" element, highlighting a distinct mechanism, but one that is functionally linked with the hammerhead ribozyme.

Crucial Roles of Two Hydrated Mg2+ Ions in Reaction Catalysis of the Pistol Ribozyme

Pistol ribozymes constitute a new class of small self‐cleaving RNAs. Crystal structures have been solved, providing three‐dimensional snapshots along the reaction coordinate of pistol phosphodiester cleavage, corresponding to the pre‐catalytic state, a vanadate mimic of the transition state, and the product. The results led to the proposed underlying chemical mechanism. Importantly, a hydrated Mg²⁺ ion remains innersphere‐coordinated to N7 of G33 in all three states, and is consistent with its likely role as acid in general acid base catalysis (δ and β catalysis). Strikingly, the new structures shed light on a second hydrated Mg²⁺ ion that approaches the scissile phosphate from its binding site in the pre‐cleavage state to reach out for water‐mediated hydrogen bonding in the cyclophosphate product. The major role of the second Mg²⁺ ion appears to be the stabilization of product conformation. This study delivers a mechanistic understanding of ribozyme‐catalyzed backbone cleavage.

Dynamical ensemble of the active state and transition state mimic for the RNA-cleaving 8-17 DNAzyme in solution

We perform molecular dynamics simulations, based on recent crystallographic data, on the 8-17 DNAzyme at four states along the reaction pathway to determine the dynamical ensemble for the active state and transition state mimic in solution. A striking finding is the diverse roles played by Na+ and Pb2+ ions in the electrostatically strained active site that impact all four fundamental catalytic strategies, and share commonality with some features recently inferred for naturally occurring hammerhead and pistol ribozymes. The active site Pb2+ ion helps to stabilize in-line nucleophilic attack, provides direct electrostatic transition state stabilization, and facilitates leaving group departure. A conserved guanine residue is positioned to act as the general base, and is assisted by a bridging Na+ ion that tunes the pKa and facilitates in-line fitness. The present work provides insight into how DNA molecules are able to solve the RNA-cleavage problem, and establishes functional relationships between the mechanism of these engineered DNA enzymes with their naturally evolved RNA counterparts. This adds valuable information to our growing body of knowledge on general mechanisms of phosphoryl transfer reactions catalyzed by RNA, proteins and DNA.

Molecular simulations of the pistol ribozyme: unifying the interpretation of experimental data and establishing functional links with the hammerhead ribozyme

The pistol ribozyme (Psr) is among the most recently discovered RNA enzymes and has been the subject of experiments aimed at elucidating the mechanism. Recent biochemical studies have revealed exciting clues about catalytic interactions in the active site not apparent from available crystallographic data. The present work unifies the interpretation of the existing body of structural and functional data on Psr by providing a dynamical model for the catalytically active state in solution from molecular simulation. Our results suggest that a catalytic Mg2+ ion makes inner-sphere contact with G33:N7 and outer-sphere coordination to the pro-RP of the scissile phosphate, promoting electrostatic stabilization of the dianionic transition state and neutralization of the developing charge of the leaving group through a metal-coordinated water molecule that is made more acidic by a hydrogen bond donated from the 2'OH of P32. This model is consistent with experimental activity-pH and mutagenesis data, including sensitivity to G33(7cG) and phosphorothioate substitution/metal ion rescue. The model suggests several experimentally testable predictions, including the response of cleavage activity to mutations at G42 and P32 positions in the ribozyme, and thio substitutions of the substrate in the presence of different divalent metal ions. Further, the model identifies striking similarities of Psr to the hammerhead ribozyme (HHr), including similar global fold, organization of secondary structure around an active site three-way junction, catalytic metal ion binding mode, and guanine general base. However, the specific binding mode and role of the Mg2+ ion, as well as a conserved 2'-OH in the active site, are interrelated but subtly different between the ribozymes.

Theory and simulations for RNA folding in mixtures of monovalent and divalent cations

RNA molecules cannot fold in the absence of counterions. Experiments are typically performed in the presence of monovalent and divalent cations. How to treat the impact of a solution containing a mixture of both ion types on RNA folding has remained a challenging problem for decades. By exploiting the large concentration difference between divalent and monovalent ions used in experiments, we develop a theory based on the reference interaction site model (RISM), which allows us to treat divalent cations explicitly while keeping the implicit screening effect due to monovalent ions. Our theory captures both the inner shell and outer shell coordination of divalent cations to phosphate groups, which we demonstrate is crucial for an accurate calculation of RNA folding thermodynamics. The RISM theory for ion-phosphate interactions when combined with simulations based on a transferable coarse-grained model allows us to predict accurately the folding of several RNA molecules in a mixture containing monovalent and divalent ions. The calculated folding free energies and ion-preferential coefficients for RNA molecules (pseudoknots, a fragment of the rRNA, and the aptamer domain of the adenine riboswitch) are in excellent agreement with experiments over a wide range of monovalent and divalent ion concentrations. Because the theory is general, it can be readily used to investigate ion and sequence effects on DNA properties.

Deducing putative ancestral forms of GNRA/receptor interactions from the ribosome

Stable RNAs rely on a vast repertoire of long-range interactions to assist in the folding of complex cellular machineries such as the ribosome. The universally conserved L39/H89 interaction is a long-range GNRA-like/receptor interaction localized in proximity to the peptidyl transferase center of the large subunit of the ribosome. Because of its central location, L39/H89 likely originated at an early evolutionary stage of the ribosome and played a significant role in its early function. However, L39/H89 self-assembly is impaired outside the ribosomal context. Herein, we demonstrate that structural modularity principles can be used to re-engineer L39/H89 to self-assemble in vitro. The new versions of L39/H89 improve affinity and loop selectivity by several orders of magnitude and retain the structural and functional features of their natural counterparts. These versions of L39/H89 are proposed to be ancestral forms of L39/H89 that were capable of assembling and folding independently from proteins and post-transcriptional modifications. This work demonstrates that novel RNA modules can be rationally designed by taking advantage of the modular syntax of RNA. It offers the prospect of creating new biochemical models of the ancestral ribosome and increases the tool kit for RNA nanotechnology and synthetic biology.

No magnesium is needed for binding of the stimulator of interferon genes to cyclic dinucleotides

Stimulator of interferon genes (STING) binds cyclic dinucleotides (CDNs), which induce a large conformational change of the protein. The structural basis of activation of STING by CDNs is rather well understood. Unliganded STING forms an open dimer that undergoes a large conformational change (∼10 Å) to a closed conformation upon the binding of a CDN molecule. This event activates downstream effectors of STING and subsequently leads to activation of the type 1 interferon response. However, a previously solved structure of STING with 3',3'-c-di-GMP shows Mg atoms mediating the interaction of STING with this CDN. Here, it is shown that no Mg atoms are needed for this interaction; in fact, magnesium can in some cases obstruct the binding of a CDN to STING.

Importance of potassium ions for ribosome structure and function revealed by long-wavelength X-ray diffraction

The ribosome, the largest RNA-containing macromolecular machinery in cells, requires metal ions not only to maintain its three-dimensional fold but also to perform protein synthesis. Despite the vast biochemical data regarding the importance of metal ions for efficient protein synthesis and the increasing number of ribosome structures solved by X-ray crystallography or cryo-electron microscopy, the assignment of metal ions within the ribosome remains elusive due to methodological limitations. Here we present extensive experimental data on the potassium composition and environment in two structures of functional ribosome complexes obtained by measurement of the potassium anomalous signal at the K-edge, derived from long-wavelength X-ray diffraction data. We elucidate the role of potassium ions in protein synthesis at the three-dimensional level, most notably, in the environment of the ribosome functional decoding and peptidyl transferase centers. Our data expand the fundamental knowledge of the mechanism of ribosome function and structural integrity.

Comparison of the Structures and Mechanisms of the Pistol and Hammerhead Ribozymes

Comparison of the secondary and three-dimensional structures of the hammerhead and pistol ribozymes reveals many close similarities, so in this work we have asked if they are mechanistically identical. We have determined a new crystal structure of the pistol ribozyme and have shown that G40 acts as general base in the cleavage reaction. The conformation in the active site ensures an in-line attack of the O2' nucleophile, and the conformation at the scissile phosphate and the position of the general base are closely similar to those in the hammerhead ribozyme. However, the two ribozymes differ in the nature of the general acid. 2'-Amino substitution experiments indicate that the general acid of the hammerhead ribozyme is the O2' of G8, while that of the pistol ribozyme is a hydrated metal ion. The two ribozymes are related but mechanistically distinct.

Structure and functional reselection of the Mango-III fluorogenic RNA aptamer

Several turn-on RNA aptamers that activate small-molecule fluorophores have been selected in vitro. Among these, the ~30 nucleotide Mango-III is notable because it binds the thiazole orange derivative TO1-Biotin with high affinity and fluoresces brightly (quantum yield 0.55). Uniquely among related aptamers, Mango-III exhibits biphasic thermal melting, characteristic of molecules with tertiary structure. We report crystal structures of TO1-Biotin complexes of Mango-III, a structure-guided mutant Mango-III(A10U), and a functionally reselected mutant iMango-III. The structures reveal a globular architecture arising from an unprecedented pseudoknot-like connectivity between a G-quadruplex and an embedded non-canonical duplex. The fluorophore is restrained into a planar conformation by the G-quadruplex, a lone, long-range trans Watson-Crick pair (whose A10U mutation increases quantum yield to 0.66), and a pyrimidine perpendicular to the nucleobase planes of those motifs. The improved iMango-III and Mango-III(A10U) fluoresce ~50% brighter than enhanced green fluorescent protein, making them suitable tags for live cell RNA visualization.

Interactions between cyclic nucleotides and common cations: an ab initio molecular dynamics study

We present the first, to the best of our knowledge, ab initio molecular dynamics (AIMD) investigation on three aqueous solutions where an abasic cyclic nucleotide model is solvated in the presence of distinct cations (i.e., Na⁺, K⁺ and Mg²⁺). We elucidate the typical modalities of interaction between those ionic species and the nucleotide moiety by first-principles numerical simulations, starting from an inner-shell binding configuration on a time scale of 100 ps (total simulation time of ∼600 ps). Whereas the strong "structure-maker" Mg²⁺ is permanently bound to one of the two oxygen atoms of the phosphate group of the nucleotide model, Na⁺ and K⁺ show binding times τ_b of 65 ps and 10-15 ps, respectively, thus reflecting their chemical nature in aqueous solutions. Furthermore, we qualitatively relate these findings to approximate free-energy barriers of the cations' unbinding obtained by means of exploratory well-tempered metadynamics. With the aim of shedding light on the features of commonly employed force-fields (FFs), classical MD simulations (almost 200 trajectories with a total simulation time of ∼18 μs) using the biomolecular AMBER FF are also reported. By choosing several combinations of the parametrization for the water environment (i.e., TIP3P, SPC/E and OPC) and cations (i.e., Joung-Cheatham, Li-Merz 12-6 and Li-Merz 12-6-4), we found significant differences in the radial distribution functions and residence times compared to the ab initio results. The Na⁺ and K⁺ ions wrongly show quasi-identical radial distribution functions and the Li & Merz 12-6-4 Lennard-Jones parameters for Mg²⁺ were found to be essential in quickly reaching the binding state consistent with AIMD.