Magnesium-binding architectures in RNA crystal structures: validation, binding preferences, classification and motif detection

Structural flexibility in the ordered domain of the dengue virus strain 2 capsid protein is critical for chaperoning viral RNA replication

Viral replication necessitates intricate nucleic acid rearrangements, including annealing and strand displacement to achieve the viral RNA functional structure. Often a single RNA chaperone performs these seemingly incompatible functions. This raises the question of what structural and dynamic features of such chaperones govern distinct RNA rearrangements. While cationic intrinsically disordered regions promote annealing by playing a charge-screening role, how the same chaperone mediates strand displacement remains elusive. Here, we investigate the annealing and strand displacement of the 5' upstream AUG region (5UAR) as chaperoned by the Dengue virus strain 2 capsid protein (Denv2C) as a model RNA chaperone. Through single molecule analysis and molecular simulations, we demonstrate that Denv2C regulates nucleic acid melting, folding, annealing, and strand displacement via flexibility in its ordered region. A mutation that renders the Denv2C ordered region rigid, converts Denv2C into a mere annealer. Our findings underscore the role of Denv2C's disordered region as a "macromolecular counterion" during RNA annealing, while a flexible ordered region is crucial for effective strand displacement.

PinMyMetal: a hybrid learning system to accurately model transition metal binding sites in macromolecules

Metal ions are vital components in many proteins for the inference and engineering of protein function, with coordination complexity linked to structural, catalytic, or regulatory roles. Modeling transition metal ions, especially in transient, reversible, and concentration-dependent regulatory sites, remains challenging. We present PinMyMetal (PMM), a hybrid machine learning system designed to accurately predict transition metal localization and environment in macromolecules, tailored to tetrahedral and octahedral geometries. PMM outperforms other predictors, achieving high accuracy in ligand and coordinate predictions. It excels in predicting regulatory sites (median deviation 0.36 Å), demonstrating superior accuracy in locating catalytic sites (0.33 Å) and structural sites (0.19 Å). Each predicted site is assigned a certainty score based on local structural and physicochemical features, independent of homologs. Interactive validation through our server, CheckMyMetal, expands PMM's scope, enabling it to pinpoint and validate diverse functional metal sites from different structure sources (predicted structures, cryo-EM, and crystallography). This facilitates residue-wise assessment and robust metal binding site design. The lightweight PMM system demands minimal computing resources and is available at https://PMM.biocloud.top . The PMM workflow can interrogate with protein sequence to characterize the localization of the most probable transition metals, which is often interchangeable and hard to differentiate by nature.

Complex water networks visualized by cryogenic electron microscopy of RNA

The stability and function of biomolecules are directly influenced by their myriad interactions with water1-16. Here we investigated water through cryogenic electron microscopy (cryo-EM) on a highly solvated molecule: the Tetrahymena ribozyme. By using segmentation-guided water and ion modelling (SWIM)17,18, an approach combining resolvability and chemical parameters, we automatically modelled and cross-validated water molecules and Mg2+ ions in the ribozyme core, revealing the extensive involvement of water in mediating RNA non-canonical interactions. Unexpectedly, in regions where SWIM does not model ordered water, we observed highly similar densities in both cryo-EM maps. In many of these regions, the cryo-EM densities superimpose with complex water networks predicted by molecular dynamics, supporting their assignment as water and suggesting a biophysical explanation for their elusiveness to conventional atomic coordinate modelling. Our study demonstrates an approach to unveil both rigid and flexible waters that surround biomolecules through cryo-EM map densities, statistical and chemical metrics, and molecular dynamics simulations.

Predicting RNA structure and dynamics with deep learning and solution scattering

Advanced deep learning and statistical methods can predict structural models for RNA molecules. However, RNAs are flexible, and it remains difficult to describe their macromolecular conformations in solutions where varying conditions can induce conformational changes. Small-angle x-ray scattering (SAXS) in solution is an efficient technique to validate structural predictions by comparing the experimental SAXS profile with those calculated from predicted structures. There are two main challenges in comparing SAXS profiles to RNA structures: the absence of cations essential for stability and charge neutralization in predicted structures and the inadequacy of a single structure to represent RNA's conformational plasticity. We introduce a solution conformation predictor for RNA (SCOPER) to address these challenges. This pipeline integrates kinematics-based conformational sampling with the innovative deep learning model, IonNet, designed for predicting Mg2+ ion binding sites. Validated through benchmarking against 14 experimental data sets, SCOPER significantly improved the quality of SAXS profile fits by including Mg2+ ions and sampling of conformational plasticity. We observe that an increased content of monovalent and bivalent ions leads to decreased RNA plasticity. Therefore, carefully adjusting the plasticity and ion density is crucial to avoid overfitting experimental SAXS data. SCOPER is an efficient tool for accurately validating the solution state of RNAs given an initial, sufficiently accurate structure and provides the corrected atomistic model, including ions.

Complex Water Networks Visualized through 2.2-2.3 Å Cryogenic Electron Microscopy of RNA

The stability and function of biomolecules are directly influenced by their myriad interactions with water. In this study, we investigated water through cryogenic electron microscopy (cryo-EM) on a highly solvated molecule, the Tetrahymena ribozyme, determined at 2.2 and 2.3 Å resolutions. By employing segmentation-guided water and ion modeling (SWIM), an approach combining resolvability and chemical parameters, we automatically modeled and cross-validated water molecules and Mg2+ ions in the ribozyme core, revealing the extensive involvement of water in mediating RNA non-canonical interactions. Unexpectedly, in regions where SWIM does not model ordered water, we observed highly similar densities in both cryo-EM maps. In many of these regions, the cryo-EM densities superimpose with complex water networks predicted by molecular dynamics (MD), supporting their assignment as water and suggesting a biophysical explanation for their elusiveness to conventional atomic coordinate modeling. Our study demonstrates an approach to unveil both rigid and flexible waters that surround biomolecules through cryo-EM map densities, statistical and chemical metrics, and MD simulations.

Metals in Motion: Understanding Labile Metal Pools in Bacteria

Metal ions are essential for all life. In microbial cells, potassium (K+) is the most abundant cation and plays a key role in maintaining osmotic balance. Magnesium (Mg2+) is the dominant divalent cation and is required for nucleic acid structure and as an enzyme cofactor. Microbes typically require the transition metals manganese (Mn), iron (Fe), copper (Cu), and zinc (Zn), although the precise set of metal ions needed to sustain life is variable. Intracellular metal pools can be conceptualized as a chemically complex mixture of rapidly exchanging (labile) ions, complemented by those reservoirs that exchange slowly relative to cell metabolism (sequestered). Labile metal pools are buffered by transient interactions with anionic metabolites and macromolecules, with the ribosome playing a major role. Sequestered metal pools include many metalloproteins, cofactors, and storage depots, with some pools redeployed upon metal depletion. Here, I review the size, composition, and dynamics of intracellular metal pools and highlight the major gaps in understanding.

Predicting RNA Structure and Dynamics with Deep Learning and Solution Scattering

Advanced deep learning and statistical methods can predict structural models for RNA molecules. However, RNAs are flexible, and it remains difficult to describe their macromolecular conformations in solutions where varying conditions can induce conformational changes. Small-angle X-ray scattering (SAXS) in solution is an efficient technique to validate structural predictions by comparing the experimental SAXS profile with those calculated from predicted structures. There are two main challenges in comparing SAXS profiles to RNA structures: the absence of cations essential for stability and charge neutralization in predicted structures and the inadequacy of a single structure to represent RNA's conformational plasticity. We introduce Solution Conformation Predictor for RNA (SCOPER) to address these challenges. This pipeline integrates kinematics-based conformational sampling with the innovative deep-learning model, IonNet, designed for predicting Mg2+ ion binding sites. Validated through benchmarking against fourteen experimental datasets, SCOPER significantly improved the quality of SAXS profile fits by including Mg2+ ions and sampling of conformational plasticity. We observe that an increased content of monovalent and bivalent ions leads to decreased RNA plasticity. Therefore, carefully adjusting the plasticity and ion density is crucial to avoid overfitting experimental SAXS data. SCOPER is an efficient tool for accurately validating the solution state of RNAs given an initial, sufficiently accurate structure and provides the corrected atomistic model, including ions.

Principles of ion binding to RNA inferred from the analysis of a 1.55 Å resolution bacterial ribosome structure - Part I: Mg2

The importance of Mg2+ ions for RNA structure and function cannot be overstated. Several attempts were made to establish a comprehensive Mg2+ binding site classification. However, such descriptions were hampered by poorly modelled ion binding sites as observed in a recent cryo-EM 1.55 Å Escherichia coli ribosome structure where incomplete ion assignments blurred our understanding of their binding patterns. We revisited this model to establish general binding principles applicable to any RNA of sufficient resolution. These principles rely on the 2.9 Å distance separating two water molecules bound in cis to Mg2+. By applying these rules, we could assign all Mg2+ ions bound with 2-4 non-water oxygens. We also uncovered unanticipated motifs where up to five adjacent nucleotides wrap around a single ion. The formation of such motifs involves a hierarchical Mg2+ ion dehydration process that plays a significant role in ribosome biogenesis and in the folding of large RNAs. Besides, we established a classification of the Mg2+…Mg2+ and Mg2+…K+ ion pairs observed in this ribosome. Overall, the uncovered binding principles enhance our understanding of the roles of ions in RNA structure and will help refining the solvation shell of other RNA systems.

High-throughput identification of calcium-regulated proteins across diverse proteomes

Calcium ions play important roles in nearly every biological process, yet whole-proteome analysis of calcium effectors has been hindered by a lack of high-throughput, unbiased, and quantitative methods to identify protein-calcium engagement. To address this, we adapted protein thermostability assays in budding yeast, human cells, and mouse mitochondria. Based on calcium-dependent thermostability, we identified 2,884 putative calcium-regulated proteins across human, mouse, and yeast proteomes. These data revealed calcium engagement of signaling hubs and cellular processes, including metabolic enzymes and the spliceosome. Cross-species comparison of calcium-protein engagement and mutagenesis experiments identified residue-specific cation engagement, even within well-known EF-hand domains. Additionally, we found that the dienoyl-coenzyme A (CoA) reductase DECR1 binds calcium at physiologically relevant concentrations with substrate-specific affinity, suggesting direct calcium regulation of mitochondrial fatty acid oxidation. These discovery-based proteomic analyses of calcium effectors establish a key resource to dissect cation engagement and its mechanistic effects across multiple species and diverse biological processes.

Metal Ion Sensing by Tetraloop-like RNA Fragment: Role of Compact Intermediates with Non-Native Metal Ion-RNA Inner-Shell Contacts

Divalent metal ions influence the folding and function of ribonucleic acid (RNA) in the cells. The mechanism of how RNA structural elements in riboswitches sense specific metal ions is unclear. RNA interacts with ions through two distinct binding modes: direct interaction between the ion and RNA (inner-shell (IS) coordination) and indirect interaction between the ion and RNA mediated through water molecules (outer-shell (OS) coordination). To understand how RNA senses metal ions such as Mg2+ and Ca2+, we studied the folding of a small RNA segment from the Mg2+ sensing M-Box riboswitch using computer simulations. This RNA segment has the characteristics of a GNRA tetraloop motif and interestingly requires the binding of a single Mg2+ ion. The folding free energy surface of this simple tetraloop system is multidimensional, with a population of multiple intermediates where the tetraloop and cation interact through IS and OS coordination. The partially folded compact tetraloop intermediates form multiple non-native IS contacts with the metal ion. Thermal fluctuations should break these strong non-native IS contacts so that the tetraloop can fold to the native state, resulting in higher folding free energy barriers. Ca2+ undergoes rapid OS to IS transitions and vice versa due to its lower charge density than Mg2+. However, the ability of Ca2+ to stabilize the native tetraloop state is weaker, as it could not hold the loop-closing nucleotides together due to its weaker interactions with the nucleotides. These insights are critical to understanding the specific ion sensing mechanisms in riboswitches, and the predictions are amenable for verification by nuclear magnetic resonance (NMR) experiments.

CheckMyMetal (CMM): validating metal-binding sites in X-ray and cryo-EM data

Identifying and characterizing metal-binding sites (MBS) within macromolecular structures is imperative for elucidating their biological functions. CheckMyMetal (CMM) is a web based tool that facilitates the interactive validation of MBS in structures determined through X-ray crystallography and cryo-electron microscopy (cryo-EM). Recent updates to CMM have significantly enhanced its capability to efficiently handle large datasets generated from cryo-EM structural analyses. In this study, we address various challenges inherent in validating MBS within both X-ray and cryo-EM structures. Specifically, we examine the difficulties associated with accurately identifying metals and modeling their coordination environments by considering the ongoing reproducibility challenges in structural biology and the critical importance of well annotated, high-quality experimental data. CMM employs a sophisticated framework of rules rooted in the valence bond theory for MBS validation. We explore how CMM validation parameters correlate with the resolution of experimentally derived structures of macromolecules and their complexes. Additionally, we showcase the practical utility of CMM by analyzing a representative cryo-EM structure. Through a comprehensive examination of experimental data, we demonstrate the capability of CMM to advance MBS characterization and identify potential instances of metal misassignment.

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.

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.

Understanding the Thermodynamics of Magnesium Binding to RNA Structural Motifs

Divalent magnesium ions (Mg2+) serve a vital role in defining the structural and catalytic chemistry of a wide array of RNA molecules. The body of structural information on RNA motifs continues to expand and, in turn, the functional importance of Mg2+ is revealed. A combination of prior work on the structural characterization of magnesium binding ligands with inner- and outer-sphere coordination modes, with recorded experimental binding energies for inner- and outer-sphere contacts, demonstrates the relative affinity and thermodynamic hierarchy for these sites. In turn, these can be correlated with cellular concentrations of free available magnesium ions, allowing the prioritization of populating important functional sites and a correlation with physiological function. This paper summarizes some of the key results of that analysis and provides predictive rules for the affinity and role of newly identified Mg binding sites on complex RNA structures. The influence of crystal packing on magnesium binding to RNA motifs, relative to their solution form, is addressed and caveats made.

Challenges and Pitfalls of Research Designs Involving Magnesium-Based Biomaterials: An Overview

Magnesium-based biomaterials hold remarkable promise for various clinical applications, offering advantages such as reduced stress-shielding and enhanced bone strengthening and vascular remodeling compared to traditional materials. However, ensuring the quality of preclinical research is crucial for the development of these implants. To achieve implant success, an understanding of the cellular responses post-implantation, proper model selection, and good study design are crucial. There are several challenges to reaching a safe and effective translation of laboratory findings into clinical practice. The utilization of Mg-based biomedical devices eliminates the need for biomaterial removal surgery post-healing and mitigates adverse effects associated with permanent biomaterial implantation. However, the high corrosion rate of Mg-based implants poses challenges such as unexpected degradation, structural failure, hydrogen evolution, alkalization, and cytotoxicity. The biocompatibility and degradability of materials based on magnesium have been studied by many researchers in vitro; however, evaluations addressing the impact of the material in vivo still need to be improved. Several animal models, including rats, rabbits, dogs, and pigs, have been explored to assess the potential of magnesium-based materials. Moreover, strategies such as alloying and coating have been identified to enhance the degradation rate of magnesium-based materials in vivo to transform these challenges into opportunities. This review aims to explore the utilization of Mg implants across various biomedical applications within cellular (in vitro) and animal (in vivo) models.

CytoSIP: an annotated structural atlas for interactions involving cytokines or cytokine receptors

Therapeutic agents targeting cytokine-cytokine receptor (CK-CKR) interactions lead to the disruption in cellular signaling and are effective in treating many diseases including tumors. However, a lack of universal and quick access to annotated structural surface regions on CK/CKR has limited the progress of a structure-driven approach in developing targeted macromolecular drugs and precision medicine therapeutics. Herein we develop CytoSIP (Single nucleotide polymorphisms (SNPs), Interface, and Phenotype), a rich internet application based on a database of atomic interactions around hotspots in experimentally determined CK/CKR structural complexes. CytoSIP contains: (1) SNPs on CK/CKR; (2) interactions involving CK/CKR domains, including CK/CKR interfaces, oligomeric interfaces, epitopes, or other drug targeting surfaces; and (3) diseases and phenotypes associated with CK/CKR or SNPs. The database framework introduces a unique tri-level SIP data model to bridge genetic variants (atomic level) to disease phenotypes (organism level) using protein structure (complexes) as an underlying framework (molecule level). Customized screening tools are implemented to retrieve relevant CK/CKR subset, which reduces the time and resources needed to interrogate large datasets involving CK/CKR surface hotspots and associated pathologies. CytoSIP portal is publicly accessible at https://CytoSIP.biocloud.top , facilitating the panoramic investigation of the context-dependent crosstalk between CK/CKR and the development of targeted therapeutic agents.

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.

PinMyMetal: A hybrid learning system to accurately model metal binding sites in macromolecules

Metal ions are vital components in many proteins for the inference and engineering of protein function, with coordination complexity linked to structural (4-residue predominate), catalytic (3-residue predominate), or regulatory (2-residue predominate) roles. Computational tools for modeling metal ions in protein structures, especially for transient, reversible, and concentration-dependent regulatory sites, remain immature. We present PinMyMetal (PMM), a sophisticated hybrid machine learning system for predicting zinc ion localization and environment in macromolecular structures. Compared to other predictors, PMM excels in predicting regulatory sites (median deviation of 0.34 Å), demonstrating superior accuracy in locating catalytic sites (median deviation of 0.27 Å) and structural sites (median deviation of 0.14 Å). PMM assigns a certainty score to each predicted site based on local structural and physicochemical features independent of homolog presence. Interactive validation through our server, CheckMyMetal, expands PMM's scope, enabling it to pinpoint and validates diverse functional zinc sites from different structure sources (predicted structures, cryo-EM and crystallography). This facilitates residue-wise assessment and robust metal binding site design. The lightweight PMM system demands minimal computing resources and is available at https://PMM.biocloud.top. While currently trained on zinc, the PMM workflow can easily adapt to other metals through expanded training data.

High-Throughput Identification of Calcium Regulated Proteins Across Diverse Proteomes

Calcium ions play important roles in nearly every biological process, yet whole-proteome analysis of calcium effectors has been hindered by lack of high-throughput, unbiased, and quantitative methods to identify proteins-calcium engagement. To address this, we adapted protein thermostability assays in the budding yeast, human cells, and mouse mitochondria. Based on calcium-dependent thermostability, we identified 2884 putative calcium-regulated proteins across human, mouse, and yeast proteomes. These data revealed calcium engagement of novel signaling hubs and cellular processes, including metabolic enzymes and the spliceosome. Cross-species comparison of calcium-protein engagement and mutagenesis experiments identified residue-specific cation engagement, even within well-known EF-hand domains. Additionally, we found that the dienoyl-CoA reductase DECR1 binds calcium at physiologically-relevant concentrations with substrate-specific affinity, suggesting direct calcium regulation of mitochondrial fatty acid oxidation. These unbiased, proteomic analyses of calcium effectors establish a key resource to dissect cation engagement and its mechanistic effects across multiple species and diverse biological processes.

Origin of ribonucleotide recognition motifs through ligand mimicry at early earth

In an RNA world, the emergence of template-specific self-replication and catalysis necessitated the presence of motifs facilitating reliable recognition between RNA molecules. What did these motifs entail, and how did they evolve into the proteinaceous RNA recognition entities observed today? Direct observation of these primordial entities is hindered by rapid degradation over geological time scales. To overcome this challenge, researchers employ diverse approaches, including scrutiny of conserved sequences and structural motifs across extant organisms and employing directed evolution experiments to generate RNA molecules with specific catalytic abilities. In this review, we delve into the theme of ribonucleotide recognition across key periods of early Earth's evolution. We explore scenarios of RNA interacting with small molecules and examine hypotheses regarding the role of minerals and metal ions in enabling structured ribonucleotide recognition and catalysis. Additionally, we highlight instances of RNA-protein mimicry in interactions with other RNA molecules. We propose a hypothesis where RNA initially recognizes small molecules and metal ions/minerals, with subsequent mimicry by proteins leading to the emergence of proteinaceous RNA binding domains.

Simultaneous Localization of Two High Affinity Divalent Metal Ion Binding Sites in the Tetracycline RNA Aptamer with Mn2+-Based Pulsed Dipolar EPR Spectroscopy

Mg2+ ions play an essential part in stabilizing the tertiary structure of nucleic acids. While the importance of these ions is well documented, their localization and elucidation of their role in the structure and dynamics of nucleic acids are often challenging. In this work, pulsed electron-electron double resonance spectroscopy (PELDOR, also known as DEER) was used to localize two high affinity divalent metal ion binding sites in the tetracycline RNA aptamer with high accuracy. For this purpose, the aptamer was labeled at different positions with a semirigid nitroxide spin label and diamagnetic Mg2+ was replaced with paramagnetic Mn2+, which did not alter the folding process or ligand binding. Out of the several divalent metal ion binding sites that are known from the crystal structure, two binding sites with high affinity were detected: one that is located at the ligand binding center and another at the J1/2 junction of the RNA.

Magnesium ions mediate ligand binding and conformational transition of the SAM/SAH riboswitch

The SAM/SAH riboswitch binds S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH) with similar affinities. Mg2+ is generally known to stabilize RNA structures by neutralizing phosphates, but how it contributes to ligand binding and conformational transition is understudied. Here, extensive molecular dynamics simulations (totaling 120 μs) predicted over 10 inner-shell Mg2+ ions in the SAM/SAH riboswitch. Six of them line the two sides of a groove to widen it and thereby pre-organize the riboswitch for ligand entry. They also form outer-shell coordination with the ligands and stabilize an RNA-ligand hydrogen bond, which effectively diminishes the selectivity between SAM and SAH. One Mg2+ ion unique to the apo form maintains the Shine-Dalgarno sequence in an autonomous mode and thereby facilitates its release for ribosome binding. Mg2+ thus plays vital roles in SAM/SAH riboswitch function.

Identification of metal ion-binding sites in RNA structures using deep learning method

Metal ion is an indispensable factor for the proper folding, structural stability and functioning of RNA molecules. However, it is very difficult for experimental methods to detect them in RNAs. With the increase of experimentally resolved RNA structures, it becomes possible to identify the metal ion-binding sites in RNA structures through in-silico methods. Here, we propose an approach called Metal3DRNA to identify the binding sites of the most common metal ions (Mg2+, Na+ and K+) in RNA structures by using a three-dimensional convolutional neural network model. The negative samples, screened out based on the analysis for binding surroundings of metal ions, are more like positive ones than the randomly selected ones, which are beneficial to a powerful predictor construction. The microenvironments of the spatial distributions of C, O, N and P atoms around a sample are extracted as features. Metal3DRNA shows a promising prediction power, generally surpassing the state-of-the-art methods FEATURE and MetalionRNA. Finally, utilizing the visualization method, we inspect the contributions of nucleotide atoms to the classification in several cases, which provides a visualization that helps to comprehend the model. The method will be helpful for RNA structure prediction and dynamics simulation study. Availability and implementation: The source code is available at https://github.com/ChunhuaLiLab/Metal3DRNA.

Magnesium ions mediate ligand binding and conformational transition of the SAM/SAH riboswitch

The SAM/SAH riboswitch binds S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH) with similar affinities. Mg 2+ is generally known to stabilize RNA structures by neutralizing phosphates, but how it contributes to ligand binding and conformational transition is understudied. Here, extensive molecular dynamics simulations (totaling 120 μs) identified over 10 inner-shell Mg 2+ ions in the SAM/SAH riboswitch. Six of them line the two sides of a groove to widen it and thereby pre-organize the riboswitch for ligand entry. They also form outer-shell coordination with the ligands and stabilize an RNA-ligand hydrogen bond, which effectively diminish the selectivity between SAM and SAH. One Mg 2+ ion unique to the apo form maintains the Shine-Dalgarno sequence in an autonomous mode and thereby facilitates its release for ribosome binding. Mg 2+ thus plays vital roles in SAM/SAH riboswitch function.

Developing Connections Between LINC00298 RNA and Alzheimer's Disease Through Mapping Its Interactome and Through Biochemical Characterization

BACKGROUND

Long non-coding RNAs are ubiquitous throughout the human system, yet many of their biological functions remain unknown. LINC00298 RNA, a long intergenic non-coding RNA, has been shown to have preferential expression in the central nervous system where it contributes to neuronal differentiation and development. Furthermore, previous research has indicated that LINC00298 RNA is known to be a genetic risk factor for the development of Alzheimer's disease.

OBJECTIVE

To biochemically characterize LINC00298 RNA and to elucidate its biological function within hippocampal neuronal cells, thereby providing a greater understanding of its role in Alzheimer's disease pathogenesis.

METHODS

LINC00298 RNA was in vitro transcribed and then subjected to structural analysis using circular dichroism, and UV-Vis spectroscopy. Additionally, affinity column chromatography was used to capture LINC00298 RNA's protein binding partners from hippocampal neuronal cells, which were then identified using liquid chromatography and mass spectrometry (LC/MS).

RESULTS

LINC00298 RNA is comprised of stem-loop secondary structural elements, with a cylindrical tertiary structure that has highly dynamic regions, which result in high positional entropy. LC/MS identified 24 proteins within the interactome of LINC00298 RNA.

CONCLUSION

Through analysis of LINC00298 RNA's 24 protein binding partners, it was determined that LINC00298 RNA may play significant roles in neuronal development, proliferation, and cellular organization. Furthermore, analysis of LINC00298 RNA's interactome indicated that LINC00298 RNA is capable of intracellular motility with dual localization in the nucleus and the cytosol. This biochemical characterization of LINC00298 RNA has shed light on its role in Alzheimer's disease pathogenesis.

Potential effects of metal ion induced two-state allostery on the regulatory mechanism of add adenine riboswitch

Riboswitches normally regulate gene expression through structural changes in response to the specific binding of cellular metabolites or metal ions. Taking add adenine riboswitch as an example, we explore the influences of metal ions (especially for K⁺ and Mg²⁺ ions) on the structure and dynamics of riboswitch aptamer (with and without ligand) by using molecular dynamic (MD) simulations. Our results show that a two-state transition marked by the structural deformation at the connection of J12 and P1 (C_J12-P1) is not only related to the binding of cognate ligands, but also strongly coupled with the change of metal ion environments. Moreover, the deformation of the structure at C_J12-P1 can be transmitted to P1 directly connected to the expression platform in multiple ways, which will affect the structure and stability of P1 to varying degrees, and finally change the regulation state of this riboswitch.

Molecular dynamic simulations are employed to assess the influence of metal ions on riboswitch structure and dynamics, suggesting a conformational control of riboswitch aptamers by metal ions before ligand binding.

Measurement of the specific and non-specific binding energies of Mg2+ to RNA

Determining the non-specific and specific electrostatic contributions of magnesium binding to RNA is a challenging problem. We introduce a single-molecule method based on measuring the folding energy of a native RNA in magnesium and at its equivalent sodium concentration. The latter is defined so that the folding energy in sodium equals the non-specific electrostatic contribution in magnesium. The sodium equivalent can be estimated according to the empirical 100/1 rule (1 M NaCl is equivalent to 10 mM MgCl2), which is a good approximation for most RNAs. The method is applied to an RNA three-way junction (3WJ) that contains specific Mg2+ binding sites and misfolds into a double hairpin structure without binding sites. We mechanically pull the RNA with optical tweezers and use fluctuation theorems to determine the folding energies of the native and misfolded structures in magnesium (10 mM MgCl2) and at the equivalent sodium condition (1 M NaCl). While the free energies of the misfolded structure are equal in magnesium and sodium, they are not for the native structure, the difference being due to the specific binding energy of magnesium to the 3WJ, which equals ΔG≃ 10 kcal/mol. Besides stabilizing the 3WJ, Mg2+ also kinetically rescues it from the misfolded structure over timescales of tens of seconds in a force-dependent manner. The method should generally be applicable to determine the specific binding energies of divalent cations to other tertiary RNAs.

Structural Bioinformatics and Deep Learning of Metalloproteins: Recent Advances and Applications

All living organisms require metal ions for their energy production and metabolic and biosynthetic processes. Within cells, the metal ions involved in the formation of adducts interact with metabolites and macromolecules (proteins and nucleic acids). The proteins that require binding to one or more metal ions in order to be able to carry out their physiological function are called metalloproteins. About one third of all protein structures in the Protein Data Bank involve metalloproteins. Over the past few years there has been tremendous progress in the number of computational tools and techniques making use of 3D structural information to support the investigation of metalloproteins. This trend has been boosted by the successful applications of neural networks and machine/deep learning approaches in molecular and structural biology at large. In this review, we discuss recent advances in the development and availability of resources dealing with metalloproteins from a structure-based perspective. We start by addressing tools for the prediction of metal-binding sites (MBSs) using structural information on apo-proteins. Then, we provide an overview of the methods for and lessons learned from the structural comparison of MBSs in a fold-independent manner. We then move to describing databases of metalloprotein/MBS structures. Finally, we summarizing recent ML/DL applications enhancing the functional interpretation of metalloprotein structures.

Petkov P, Milenov TI, Vayssilov GN. Sodium and Magnesium Ion Location at the Backbone and at the Nucleobase of RNA: Ab Initio Molecular Dynamics in Water Solution

The interactions between Na⁺ or Mg²⁺ ions with different parts of single-stranded RNA molecules, namely, the oxygen atoms from the phosphate groups or the guanine base, in water solution have been studied using first-principles molecular dynamics. Sodium ions were found to be much more mobile than Mg²⁺ ions and readily underwent transitions between a state directly bonded to RNA oxygen atoms and a completely solvated state. The inner solvation shell of Na⁺ ions fluctuated stochastically at a femtosecond timescale coordinating on average 5 oxygen atoms for bonded Na⁺ ions and 5.5 oxygen atoms for solvated Na⁺ ions. In contrast, the inner solvation shell of Mg²⁺ ions was stable in both RNA-bonded and completely solvated states. In both cases, Mg²⁺ ions coordinated 6 oxygen atoms from the inner solvation shell. Consistent with their stable solvation shells, Mg²⁺ ions were more effective than Na⁺ ions in stabilizing the RNA backbone conformation. The exclusion zones between the first and second solvation shells, solvation shell widths, and angles for binding to carbonyl oxygen of guanine for solvated Na⁺ or Mg²⁺ ions exhibited a number of quantitative differences when compared with RNA crystallographic data. The presented results support the distinct capacity of Mg²⁺ ions to support the RNA structure not only in the crystal phase but also in the dynamic water environment both on the side of the phosphate moiety and on the side of the nucleobase.

Interchangeable utilization of metals: New perspectives on the impacts of metal ions employed in ancient and extant biomolecules

Metal ions provide considerable functionality across biological systems, and their utilization within biomolecules has adapted through changes in the chemical environment to maintain the activity they facilitate. While ancient earth's atmosphere was rich in iron and manganese and low in oxygen, periods of atmospheric oxygenation significantly altered the availability of certain metal ions, resulting in ion replacement within biomolecules. This adaptation mechanism has given rise to the phenomenon of metal cofactor interchangeability, whereby contemporary proteins and nucleic acids interact with multiple metal ions interchangeably, with different coordinated metals influencing biological activity, stability, and toxic potential. The ability of extant organisms to adapt to fluctuating metal availability remains relevant in a number of crucial biomolecules, including the superoxide dismutases of the antioxidant defense systems and ribonucleotide reductases. These well-studied and ancient enzymes illustrate the potential for metal interchangeability and adaptive utilization. More recently, the ribosome has also been demonstrated to exhibit interchangeable interactions with metal ions with impacts on function, stability, and stress adaptation. Using these and other examples, here we review the biological significance of interchangeable metal ions from a new angle that combines both biochemical and evolutionary viewpoints. The geochemical pressures and chemical properties that underlie biological metal utilization are discussed in the context of their impact on modern disease states and treatments.

Structure and Mechanical Stabilities of the Three-Way Junction Motifs in Prohead RNA

The core structure of phi29 prohead RNA (pRNA) is composed of three major helices organized into three-way junction pRNA (3WJ-pRNA) and has stout structural rigidity along the coaxial helices. Prohead RNAs of the other Bacillus subtilis bacteriophages such as GA1 and SF5 share similar secondary structure and function with phi29; whether these pRNAs have similar mechanical rigidity remains to be elucidated. In this study, we constructed the tertiary structures of GA1 and SF5 3WJ-pRNAs by comparative modeling. Both GA1 and SF5 3WJ-pRNAs adopt a similar structure, in which three helices are organized as the three-way junction and two of the three helices are stacked coaxially. Moreover, detailed structural features of GA1 and SF5 3WJ-pRNAs are also similar to those of phi29 3WJ-pRNA: all of the bases of the coaxial helices are paired, and all of the adenines in the junction region are paired, which eliminates the interference of A-minor tertiary interactions. Hence, the coaxial helices tightly join to each other, and the major groove between them is very narrow. Two Mg²⁺ ions can thus fit into this major groove and form double Mg clamps. A steered molecular dynamics simulation was used to study the mechanical properties of these 3WJ-pRNAs. Both GA1 and SF5 3WJ-pRNAs show strong resistance to applied force in the direction of their coaxial helices. Such mechanical stability can be attributed to the Mg clamps.

virusMED: an atlas of hotspots of viral proteins

Metal binding sites, antigen epitopes and drug binding sites are the hotspots in viral proteins that control how viruses interact with their hosts. virusMED (virus Metal binding sites, Epitopes and Drug binding sites) is a rich internet application based on a database of atomic interactions around hotspots in 7041 experimentally determined viral protein structures. 25306 hotspots from 805 virus strains from 75 virus families were characterized, including influenza, HIV-1 and SARS-CoV-2 viruses. Just as Google Maps organizes and annotates points of interest, virusMED presents the positions of individual hotspots on each viral protein and creates an atlas upon which newly characterized functional sites can be placed as they are being discovered. virusMED contains an extensive set of annotation tags about the virus species and strains, viral hosts, viral proteins, metal ions, specific antibodies and FDA-approved drugs, which permits rapid screening of hotspots on viral proteins tailored to a particular research problem. The virusMED portal (https://virusmed.biocloud.top) can serve as a window to a valuable resource for many areas of virus research and play a critical role in the rational design of new preventative and therapeutic agents targeting viral infections.

IonoBiology: The functional dynamics of the intracellular metallome, with lessons from bacteria

Metal ions are essential for life and represent the second most abundant constituent (after water) of any living cell. While the biological importance of inorganic ions has been appreciated for over a century, we are far from a comprehensive understanding of the functional roles that ions play in cells and organisms. In particular, recent advances are challenging the traditional view that cells maintain constant levels of ion concentrations (ion homeostasis). In fact, the ionic composition (metallome) of cells appears to be purposefully dynamic. The scientific journey that started over 60 years ago with the seminal work by Hodgkin and Huxley on action potentials in neurons is far from reaching its end. New evidence is uncovering how changes in ionic composition regulate unexpected cellular functions and physiology, especially in bacteria, thereby hinting at the evolutionary origins of the dynamic metallome. It is an exciting time for this field of biology, which we discuss and refer to here as IonoBiology.

Oxidative Stress in Bacteria and the Central Dogma of Molecular Biology

Ever since the "great oxidation event," Earth's cellular life forms had to cope with the danger of reactive oxygen species (ROS) affecting the integrity of biomolecules and hampering cellular metabolism circuits. Consequently, increasing ROS levels in the biosphere represented growing stress levels and thus shaped the evolution of species. Whether the ROS were produced endogenously or exogenously, different systems evolved to remove the ROS and repair the damage they inflicted. If ROS outweigh the cell's capacity to remove the threat, we speak of oxidative stress. The injuries through oxidative stress in cells are diverse. This article reviews the damage oxidative stress imposes on the different steps of the central dogma of molecular biology in bacteria, focusing in particular on the RNA machines involved in transcription and translation.

Covid-19.bioreproducibility.org: A web resource for SARS-CoV-2-related structural models

The COVID-19 pandemic has triggered numerous scientific activities aimed at understanding the SARS-CoV-2 virus and ultimately developing treatments. Structural biologists have already determined hundreds of experimental X-ray, cryo-EM, and NMR structures of proteins and nucleic acids related to this coronavirus, and this number is still growing. To help biomedical researchers, who may not necessarily be experts in structural biology, navigate through the flood of structural models, we have created an online resource, covid19.bioreproducibility.org, that aggregates expert-verified information about SARS-CoV-2-related macromolecular models. In this article, we describe this web resource along with the suite of tools and methodologies used for assessing the structures presented therein.

Understanding the characteristics of nonspecific binding of drug-like compounds to canonical stem-loop RNAs and their implications for functional cellular assays

Identifying small molecules that selectively bind an RNA target while discriminating against all other cellular RNAs is an important challenge in RNA-targeted drug discovery. Much effort has been directed toward identifying drug-like small molecules that minimize electrostatic and stacking interactions that lead to nonspecific binding of aminoglycosides and intercalators to many stem-loop RNAs. Many such compounds have been reported to bind RNAs and inhibit their cellular activities. However, target engagement and cellular selectivity assays are not routinely performed, and it is often unclear whether functional activity directly results from specific binding to the target RNA. Here, we examined the propensities of three drug-like compounds, previously shown to bind and inhibit the cellular activities of distinct stem-loop RNAs, to bind and inhibit the cellular activities of two unrelated HIV-1 stem-loop RNAs: the transactivation response element (TAR) and the rev response element stem IIB (RREIIB). All compounds bound TAR and RREIIB in vitro, and two inhibited TAR-dependent transactivation and RRE-dependent viral export in cell-based assays while also exhibiting off-target interactions consistent with nonspecific activity. A survey of X-ray and NMR structures of RNA-small molecule complexes revealed that aminoglycosides and drug-like molecules form hydrogen bonds with functional groups commonly accessible in canonical stem-loop RNA motifs, in contrast to ligands that specifically bind riboswitches. Our results demonstrate that drug-like molecules can nonspecifically bind stem-loop RNAs most likely through hydrogen bonding and electrostatic interactions and reinforce the importance of assaying for off-target interactions and RNA selectivity in vitro and in cells when assessing novel RNA-binders.

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.

Molecular mechanisms underlying the extreme mechanical anisotropy of the flaviviral exoribonuclease-resistant RNAs (xrRNAs)

Mechanical anisotropy is an essential property for many biomolecules to assume their structures, functions and applications, however, the mechanisms for their direction-dependent mechanical responses remain elusive. Herein, by using a single-molecule nanopore sensing technique, we explore the mechanisms of directional mechanical stability of the xrRNA1 RNA from ZIKA virus (ZIKV), which forms a complex ring-like architecture. We reveal extreme mechanical anisotropy in ZIKV xrRNA1 which highly depends on Mg²⁺ and the key tertiary interactions. The absence of Mg²⁺ and disruption of the key tertiary interactions strongly affect the structural integrity and attenuate mechanical anisotropy. The significance of ring structures in RNA mechanical anisotropy is further supported by steered molecular dynamics simulations in combination with force distribution analysis. We anticipate the ring structures can be used as key elements to build RNA-based nanostructures with controllable mechanical anisotropy for biomaterial and biomedical applications.

Mg2+ Sensing by an RNA Fragment: Role of Mg2+-Coordinated Water Molecules

RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg²⁺, Ca²⁺, and K⁺) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg²⁺ solution. We show that the rigid network of Mg²⁺-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K⁺ and Ca²⁺ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.

A Distinct Motif in a Prokaryotic Small Ras-Like GTPase Highlights Unifying Features of Walker B Motifs in P-Loop NTPases

A hallmark of the catalytically essential Walker B motif of P-loop NTPases is the presence of an acidic residue (aspartate/glutamate) for efficient Mg²⁺ coordination. Although the Walker B motif has been identified in well-studied examples of P-loop NTPases, its identity is ambiguous in many families, for example, in the prokaryotic small Ras-like GTPase family of MglA. MglA, belonging to TRAFAC class of P-loop NTPases, possesses a threonine at the position equivalent to Walker B aspartate in eukaryotic Ras-like GTPases. To resolve the identity of the Walker B residue in MglA, we carried out a comprehensive analysis of Mg²⁺ coordination on P-loop NTPase structures. Atoms in the octahedral coordination of Mg²⁺ and their interactions comprise a network including water molecules, Walker A, Walker B and switch motifs of P-loop NTPases. Based on the conserved geometry of Mg²⁺ coordination, we confirm that a conserved aspartate functions as the Walker B residue of MglA, and validate it through mutagenesis and biochemical characterization. Location of the newly identified aspartate is spatially equivalent to the Walker B residue of the ASCE division of P-loop NTPases. Furthermore, similar to the allosteric regulation of the Walker B aspartate conformation in MglA, we identify protein families in which large conformational changes involving Walker B motif potentially function as allosteric regulators. The study unravels conserved features of Mg²⁺ coordination among divergent families of P-loop NTPases, especially between ancient Ras-like GTPases and ASCE family of ATPases. The conserved geometric features provide a foundation for design of nucleotide-hydrolyzing enzymes.

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.

Structural studies of geranylgeranylglyceryl phosphate synthase, a prenyltransferase found in thermophilic Euryarchaeota

Archaea are uniquely adapted to thrive in harsh environments, and one of these adaptations involves the archaeal membrane lipids, which are characterized by their isoprenoid alkyl chains connected via ether linkages to glycerol 1-phosphate. The membrane lipids of the thermophilic and acidophilic euryarchaeota Thermoplasma volcanium are exclusively glycerol dibiphytanyl glycerol tetraethers. The first committed step in the biosynthetic pathway of these archaeal lipids is the formation of the ether linkage between glycerol 1-phosphate and geranylgeranyl diphosphate, and is catalyzed by the enzyme geranylgeranylglyceryl phosphate synthase (GGGPS). The 1.72 Å resolution crystal structure of GGGPS from T. volcanium (TvGGGPS) in complex with glycerol and sulfate is reported here. The crystal structure reveals TvGGGPS to be a dimer, which is consistent with the absence of the aromatic anchor residue in helix α5a that is required for hexamerization in other GGGPS homologs; the hexameric quaternary structure in GGGPS is thought to provide thermostability. A phylogenetic analysis of the Euryarchaeota and a parallel ancestral state reconstruction investigated the relationship between optimal growth temperature and the ancestral sequences. The presence of an aromatic anchor residue is not explained by temperature as an ecological parameter. An examination of the active site of the TvGGGPS dimer revealed that it may be able to accommodate longer isoprenoid substrates, supporting an alternative pathway of isoprenoid membrane-lipid synthesis.

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.

In silico design and validation of high-affinity RNA aptamers targeting epithelial cellular adhesion molecule dimers

Nucleic acid aptamers hold great promise for therapeutic applications due to their favorable intrinsic properties, as well as high-throughput experimental selection techniques. Despite the utility of the systematic evolution of ligands by the exponential enrichment (SELEX) method for aptamer determination, complementary in silico aptamer design is highly sought after to facilitate virtual screening and increased understanding of important nucleic acid-protein interactions. Here, with a combined experimental and theoretical approach, we have developed two optimal epithelial cellular adhesion molecule (EpCAM) aptamers. Our structure-based in silico method first predicts their binding modes and then optimizes them for EpCAM with molecular dynamics simulations, docking, and free energy calculations. Our isothermal titration calorimetry experiments further confirm that the EpCAM aptamers indeed exhibit enhanced affinity over a previously patented nanomolar aptamer, EP23. Moreover, our study suggests that EP23 and the de novo designed aptamers primarily bind to EpCAM dimers (and not monomers, as hypothesized in previous published works), suggesting a paradigm for developing EpCAM-targeted therapies.

Visualizing Disordered Single-Stranded RNA: Connecting Sequence, Structure, and Electrostatics

Disordered homopolymeric regions of single-stranded RNA, such as U or A tracts, are found within functional RNAs where they play distinct roles in defining molecular structure and facilitating recognition by partners. Despite this prominence, details of conformational and biophysical properties of these regions have not yet been resolved. We apply a number of experimental techniques to investigate the conformations of these biologically important motifs and provide quantitative measurements of their ion atmospheres. Single strands of RNA display pronounced sequence-dependent conformations that relate to the unique ion atmospheres each attracts. Chains of rU bases are relatively unstructured under all conditions, while chains of rA bases display distinct ordering through stacking or clustering motifs, depending on the composition of the surrounding solution. These dramatic structural differences are consistent with the measured disparity in ion composition and atmospheres around each homopolymer, revealing a complex interplay of base, ion, and single-strand ordering. The unique structural and ionic signatures of homopolymer ssRNAs explains their role(s) in folding structured RNAs and may explain their distinct recognition by protein partners.

Metal Coordination Is Crucial for Geranylgeranyl Diphosphate Synthase-Bisphosphonate Interactions: A Crystallographic and Computational Analysis

Geranylgeranyl diphosphate synthase (GGPPS) is a central metalloenzyme in the mevalonate pathway, crucial for the prenylation of small GTPases. As small GTPases are pivotal for cellular survival, GGPPS was highlighted as a potential target for treating human diseases, including solid and hematologic malignancies and parasitic infections. Most available GGPPS inhibitors are bisphosphonates, but the clinically available compounds demonstrate poor pharmacokinetic properties. Although the design of novel bisphosphonates with improved physicochemical properties is highly desirable, the structure of wild-type human GGPPS (hGGPPS) bound to a bisphosphonate has not been resolved. Moreover, various metal-bisphosphonate-binding stoichiometries were previously reported in structures of yeast GGPPS (yGGPPS), hampering computational drug design with metal-binding pharmacophores (MBP). In this study, we report the 2.2 Å crystal structure of hGGPPS in complex with ibandronate, clearly depicting the involvement of three Mg2+ ions in bisphosphonate-protein interactions. Using drug-binding assays and computational docking, we show that the assignment of three Mg2+ ions to the binding site of both hGGPPS and yGGPPS greatly improves the correlation between calculated binding energies and experimentally measured affinities. This work provides a structural basis for future rational design of additional MBP-harboring drugs targeting hGGPPS. SIGNIFICANCE STATEMENT: Bisphosphonates are inhibitors of geranylgeranyl diphosphate synthase (GGPPS), a metalloenzyme crucial for cell survival. Bisphosphonate binding depends on coordination by Mg2+ ions, but various Mg2+-bisphosphonate-binding stoichiometries were previously reported. In this study, we show that three Mg2+ ions are vital for drug binding and provide a structural basis for future computational design of GGPPS inhibitors.

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.

Binding of metal ions and water molecules to nucleic acid bases: the influence of water molecule coordination to a metal ion on water-nucleic acid base hydrogen bonds

The interactions of nucleic acid bases with non-coordinated and coordinated water molecules were studied by analyzing data in the Protein Data Bank (PDB) and by quantum chemical calculations. The analysis of the data in the crystal structures from the PDB indicates that hydrogen bonds involving oxygen or nitrogen atoms of nucleic acid bases and water molecules are shorter when water is bonded to a metal ion. These results are in agreement with the quantum chemical calculations on geometries and interaction energies of hydrogen bonds; the calculations on model systems show that hydrogen bonds of nucleic acid bases with water bonded to a metal ion are stronger than hydrogen bonds with non-coordinated water. These calculated values are similar to the strength of hydrogen bonds between nucleic acid bases. The results presented in this paper may be relevant to understand the role of water molecules and metal ions in the process of replication and stabilization of nucleic acids and also to understand the possible toxicity of metal ion interactions with nucleic acids.

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.