Effect of Guanine to Inosine Substitution on Stability of Canonical DNA and RNA Duplexes: Molecular Dynamics Thermodynamics Integration Study

Structural basis of water-mediated cis Watson-Crick/Hoogsteen base-pair formation in non-CpG methylation

Non-CpG methylation is associated with several cellular processes, especially neuronal development and cancer, while its effect on DNA structure remains unclear. We have determined the crystal structures of DNA duplexes containing -CGCCG- regions as CCG repeat motifs that comprise a non-CpG site with or without cytosine methylation. Crystal structure analyses have revealed that the mC:G base-pair can simultaneously form two alternative conformations arising from non-CpG methylation, including a unique water-mediated cis Watson-Crick/Hoogsteen, (w)cWH, and Watson-Crick (WC) geometries, with partial occupancies of 0.1 and 0.9, respectively. NMR studies showed that an alternative conformation of methylated mC:G base-pair at non-CpG step exhibits characteristics of cWH with a syn-guanosine conformation in solution. DNA duplexes complexed with the DNA binding drug echinomycin result in increased occupancy of the (w)cWH geometry in the methylated base-pair (from 0.1 to 0.3). Our structural results demonstrated that cytosine methylation at a non-CpG step leads to an anti→syntransition of its complementary guanosine residue toward the (w)cWH geometry as a partial population of WC, in both drug-bound and naked mC:G base pairs. This particular geometry is specific to non-CpG methylated dinucleotide sites in B-form DNA. Overall, the current study provides new insights into DNA conformation during epigenetic regulation.

Biomolecular dynamics with machine-learned quantum-mechanical force fields trained on diverse chemical fragments

The GEMS method enables molecular dynamics simulations of large heterogeneous systems at ab initio quality.

A-to-I edited miR-411-5p targets MET and promotes TKI response in NSCLC-resistant cells

Non-small cell lung cancer (NSCLC) patients carrying an epidermal growth factor receptor (EGFR) mutation have an initial favorable clinical response to the tyrosine kinase inhibitors (TKIs). Unfortunately, rapid resistance occurs mainly because of genetic alterations, including amplification of the hepatocyte growth factor receptor (MET) and its abnormal activity. The RNA post-transcriptional modifications that contribute to aberrant expression of MET in cancer are largely under-investigated and among them is the adenosine-to-inosine (A-to-I) RNA editing of microRNAs. A reduction of A-to-I editing in position 5 of miR-411-5p has been identified in several cancers, including NSCLC. In this study, thanks to cancer-associated gene expression analysis, we assessed the effect of the edited miR-411-5p on NSCLC cell lines. We found that edited miR-411-5p directly targets MET and negatively affects the mitogen-activated protein kinases (MAPKs) pathway. Considering the predominant role of the MAPKs pathway on TKIs resistance, we generated NSCLC EGFR mutated cell lines resistant to TK inhibitors and evaluated the effect of edited miR-411-5p overexpression. We found that the edited miR-411-5p reduces proliferation and induces apoptosis, promoting EGFR TKIs response in NSCLC-resistant cells.

Systematically Modulating Aptamer Affinity and Specificity by Guanosine-to-Inosine Substitution

Aptamers are widely used in small molecule detection applications due to their specificity, stability, and cost effectiveness. One key challenge in utilizing aptamers in sensors is matching the binding affinity of the aptamer to the desired concentration range for analyte detection. The most common methods for modulating affinity have inherent limitations, such as the likelihood of drastic changes in aptamer folding. Here, we propose that substituting guanosine for inosine at specific locations in the aptamer sequence provides a less perturbative approach to modulating affinity. Inosine is a naturally occurring nucleotide that results from hydrolytic deamination of adenosine, and like guanine, it base pairs with cytosine. Using the well-studied cocaine binding aptamer, we systematically replaced guanosine with inosine and were able to generate sequences having a range of binding affinities from 230 nM to 80 μM. Interestingly, we found that these substitutions could also modulate the specificity of the aptamers, leading to a range of binding affinities for structurally related analytes. Analysis of folding stability via melting temperature shows that, as expected, aptamer structure is impacted by guanosine-to-inosine substitutions. The ability to tune binding affinity and specificity through guanosine-to-inosine substitution provides a convenient and reliable approach for rapidly generating aptamers for diverse biosensing applications.

Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA

Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.

Mutual Protein-Ligand Conformational Selection Drives cGMP vs. cAMP Selectivity in Protein Kinase G

Protein kinase G (PKG) is a major receptor of cGMP, and controls signaling pathways distinct from those regulated by cAMP. However, the contributions of the two substituents that differentiate cGMP from cAMP (i.e. 6-oxo and 2-NH₂) to the cGMP-versus-cAMP selectivity of PKG remain unclear. Here, using NMR to map how binding affinity and dynamics of the protein and ligand vary along a ligand double-substitution cycle, we show that the contributions of the two substituents to binding affinity are surprisingly non-additive. Such non-additivity stems primarily from mutual protein-ligand conformational selection, whereby not only does the ligand select for a preferred protein conformation upon binding, but also, the protein selects for a preferred ligand conformation. The 6-oxo substituent mainly controls the conformational equilibrium of the bound protein, while the 2-NH₂ substituent primarily controls the conformational equilibrium of the unbound ligand (i.e. syn versus anti). Therefore, understanding the conformational dynamics of both the protein and ligand is essential to explain the cGMP-versus-cAMP selectivity of PKG.

Recognition of N6-Methyladenosine by the YTHDC1 YTH Domain Studied by Molecular Dynamics and NMR Spectroscopy: The Role of Hydration

The YTH domain of YTHDC1 belongs to a class of protein "readers", recognizing the N6-methyladenosine (m⁶A) chemical modification in mRNA. Static ensemble-averaged structures revealed details of N6-methyl recognition via a conserved aromatic cage. Here, we performed molecular dynamics (MD) simulations along with nuclear magnetic resonance (NMR) and isothermal titration calorimetry (ITC) to examine how dynamics and solvent interactions contribute to the m⁶A recognition and negative selectivity toward an unmethylated substrate. The structured water molecules surrounding the bound RNA and the methylated substrate's ability to exclude bulk water molecules contribute to the YTH domain's preference for m⁶A. Intrusions of bulk water deep into the binding pocket disrupt binding of unmethylated adenosine. The YTHDC1's preference for the 5'-Gm⁶A-3' motif is partially facilitated by a network of water-mediated interactions between the 2-amino group of the guanosine and residues in the m⁶A binding pocket. The 5'-Im⁶A-3' (where I is inosine) motif can be recognized too, but disruption of the water network lowers affinity. The D479A mutant also disrupts the water network and destabilizes m⁶A binding. Our interdisciplinary study of the YTHDC1 protein-RNA complex reveals an unusual physical mechanism by which solvent interactions contribute toward m⁶A recognition.

Free-Energy Calculation of Ribonucleic Inosines and Its Application to Nearest-Neighbor Parameters

Can current simulations quantitatively predict the stability of ribonucleic acids (RNAs)? In this research, we apply a free-energy perturbation simulation of RNAs containing inosine, a modified ribonucleic base, to the derivation of RNA nearest-neighbor parameters. A parameter set derived solely from 30 simulations was used to predict the free-energy difference of the RNA duplex with a mean unbiased error of 0.70 kcal/mol, which is a level of accuracy comparable to that obtained with parameters derived from 25 experiments. We further show that the error can be lowered to 0.60 kcal/mol by combining the simulation-derived free-energy differences with experimentally measured differences. This protocol can be used as a versatile method for deriving nearest-neighbor parameters of RNAs with various modified bases.

Direct DNA Methylation Profiling with an Electric Biosensor

DNA methylation is one of the principal epigenetic mechanisms that control gene expression in humans, and its profiling provides critical information about health and disease. Current profiling methods require chemical modification of bases followed by sequencing, which is expensive and time-consuming. Here, we report a direct and rapid determination of DNA methylation using an electric biosensor. The device consists of a DNA-tweezer probe integrated on a graphene field-effect transistor for label-free, highly sensitive, and specific methylation profiling. The device performance was evaluated with a target DNA that harbors a sequence of the methylguanine-DNA methyltransferase, a promoter of glioblastoma multiforme, a lethal brain tumor. The results show that we successfully profiled the methylated and nonmethylated forms at picomolar concentrations. Further, fluorescence kinetics and molecular dynamics simulations revealed that the position of the methylation site(s), their proximity, and accessibility to the toe-hold region of the tweezer probe are the primary determinants of the device performance.

Inosine 5'-diphosphate, a molecular decoy rescues Nucleoside diphosphate kinase from c-MYC G-Quadruplex unfolding

BACKGROUND

The transcription-inhibitory G-Quadruplex(Pu27-GQ) at c-MYC promoter is challenging to target due to structural heterogeneity. Nucleoside diphosphate kinase (NM23-H2) specifically binds and unfolds Pu27-GQ to increase c-MYC transcription. Here, we used Inosine 5'-diphosphate (IDP) to disrupt NM23-H2-Pu27-GQ interactions and arrest c-MYC transcription without compromising NM23-H2-mediated kinase properties.

METHODS

Site-directed mutagenesis,³¹P-NMR and STD-NMR studies delineate the epitope of NM23-H2-IDP complex and characterize specific amino acids in NM23-H2 involved in Pu27-GQ and IDP interactions. Immunoprecipitations and phosphohistidine-immunoblots reveal how IDP blocks NM23-H2-Pu27 association to downregulate c-MYC transcription in MDAMB-231 cells exempting NM23-H2-mediated kinase properties.

RESULTS

NMR studies show that IDP binds to the Guanosine diphosphate-binding pocket of NM23-H2 (K_D = 5.0 ± 0.276 μM). Arg88-driven hydrogen bonds to the terminal phosphate of IDP restricts P-O-P bond-rotation increasing its pKa (∆pKa = 0.85 ± 0.0025).9-inosinyl moiety of IDP is stacked over Phe60 phenyl ring driving trans-conformation of inosine and axial geometry of pyrophosphates. Chromatin immunoprecipitations revealed that these interactions rescue NM23-H2-driven Pu27-GQ unfolding, which triggers Nucleolin recruitment and lowers Sp1 occupancy at c-MYC promoter stabilizing Pu27-GQ. This silences c-MYC transcription that reduces c-MYC-Sp1 association amplifying Sp1 recruitment across P21 promoter stimulating P21 transcription and G₂/M arrest.

CONCLUSIONS

IDP synergizes the effects of Pu27-GQ-interacting compounds to abrogate c-MYC transcription and induce apoptosis in MDAMB-231 cells by disrupting NM23-H2-Pu27-GQ interactions without affecting NM23-H2-mediated kinase properties.

GENERAL SIGNIFICANCE

Our study provides a pragmatic approach for developing NM23-H2-targeting regulators to rescue NM23-H2 binding at structurally ambiguous Pu27-GQ that synergizes the anti-tumorigenic effects of GQ-based therapeutics with minimized off-target effects.

Molecular Dynamics Study of the Hybridization between RNA and Modified Oligonucleotides

MicroRNAs (miRNAs) are attractive drug candidates for many diseases as they can modulate the expression of gene networks. Recently, we discovered that DNAs targeting microRNA-22-3p (miR-22-3p) hold the potential for treating obesity and related metabolic disorders (type 2 diabetes mellitus, hyperlipidemia, and nonalcoholic fatty liver disease (NAFLD)) by turning fat-storing white adipocytes into fat-burning adipocytes. In this work, we explored the effects of chemical modifications, including phosphorothioate (PS), locked nucleic acid (LNA), and peptide nucleic acid (PNA), on the structure and energy of DNA analogs by using molecular dynamics (MD) simulations. To achieve a reliable prediction of the hybridization free energy, the AMOEBA polarizable force field and the free energy perturbation technique were employed. The calculated hybridization free energies are generally compatible with previous experiments. For LNA and PNA, the enhanced duplex stability can be explained by the preorganization mechanism, i.e., the single strands adopt stable helical structures similar to those in the duplex. For PS, the S and R isomers (Sp and Rp) have preferences for C2'-endo and C3'-endo sugar puckering conformations, respectively, and therefore Sp is less stable than Rp in DNA/RNA hybrids. In addition, the solvation penalty of Rp accounts for its destabilization effect. PS-LNA is similar to LNA as the sugar puckering is dominated by the locked sugar ring. This work demonstrated that MD simulations with polarizable force fields are useful for the understanding and design of modified nucleic acids.

NMR solution structure of an asymmetric intermolecular leaped V-shape G-quadruplex: selective recognition of the d(G2NG3NG4) sequence motif by a short linear G-rich DNA probe

Aside from classical loops among G-quadruplexes, the unique leaped V-shape scaffold spans over three G-tetrads, without any intervening residues. This scaffold enables a sharp reversal of two adjacent strand directions and simultaneously participates in forming the G-tetrad core. These features make this scaffold itself distinctive and thus an essentially more accessible target. As an alternative to the conventional antisense method using a complementary chain, forming an intermolecular G-quadruplex from two different oligomers, in which the longer one as the target is captured by a short G-rich fragment, could be helpful for recognizing G-rich sequences and structural motifs. However, such an intermolecular leaped V-shape G-quadruplex consisting of DNA oligomers of quite different lengths has not been evaluated. Here, we present the first nuclear magnetic resonance (NMR) study of an asymmetric intermolecular leaped V-shape G-quadruplex assembled between an Oxytricha nova telomeric sequence d(G2T4G4T4G4) and a single G-tract fragment d(TG4A). Furthermore, we explored the selectivity of this short fragment as a potential probe, examined the kinetic discrimination for probing a specific mutant, and proposed the key sequence motif d(G2NG3NG4) essential for building the leaped V-shape G-quadruplexes.

DNA binding preferences of S. cerevisiae RNA polymerase I Core Factor reveal a preference for the GC-minor groove and a conserved binding mechanism

In Saccharomyces cerevisiae, Core Factor (CF) is a key evolutionarily conserved transcription initiation factor that helps recruit RNA polymerase I (Pol I) to the ribosomal DNA (rDNA) promoter. Upregulated Pol I transcription has been linked to many cancers, and targeting Pol I is an attractive and emerging anti-cancer strategy. Using yeast as a model system, we characterized how CF binds to the Pol I promoter by electrophoretic mobility shift assays (EMSA). Synthetic DNA competitors along with anti-tumor drugs and nucleic acid stains that act as DNA groove blockers were used to discover the binding preference of yeast CF. Our results show that CF employs a unique binding mechanism where it prefers the GC-rich minor groove within the rDNA promoter. In addition, we show that yeast CF is able to bind to the human rDNA promoter sequence that is divergent in DNA sequence and demonstrate CF sensitivity to the human specific Pol I inhibitor, CX-5461. Finally, we show that the human Core Promoter Element (CPE) can functionally replace the yeast Core Element (CE) in vivo when aligned by conserved DNA structural features rather than DNA sequence. Together, these findings suggest that the yeast CF and the human ortholog Selectivity Factor 1 (SL1) use an evolutionarily conserved, structure-based mechanism to target DNA. Their shared mechanism may offer a new avenue in using yeast to explore current and future Pol I anti-cancer compounds.

Substituting Inosine for Guanosine in DNA: Structural and Dynamic Consequences

Inosine differs from the guanosine nucleoside only by the absence of the N2 amino group. Both nucleosides also have similar electrostatic potentials. Therefore, substituting I for G has been used to probe various properties of nucleic acids and to facilitate the interpretation of binding studies. In particular, the absence of the amino group permits the assessment of its importance in the binding of ligands to the minor groove of duplex DNA. It has been known for some time that an I-C base pair is of lower stability than a regular G-C base pair, which needs to be considered when making DNA constructs containing inosine. However, it is generally assumed that both base pairs are structurally highly similar. To test this assumption in an identical sequence environment, we have determined the fine structure of two hairpin DNA substrates that differ only in the substitution of an I-C base pair for a G-C base pair. The structures have been solved using nuclear magnetic resonance (NMR) restraints in conjunction with Mardigras and molecular dynamics. The structural data are complemented with thermodynamic and dynamic data to get a comprehensive evaluation of the consequences of G-C vs I-C base pair substitutions. Our data show a strong similarity in the structures of the hairpins, but a significant difference in the melting temperatures, T m. This difference is also reflected in the drastically decreased base pair lifetime of 7.4 milliseconds compared to the G-C base pair lifetime of 155 milliseconds. The substitution of I-C for G-C is to probe for specific effect due to the amino group is satisfactory, as long as the lowered thermal stability and the drastically increased local dynamics are considered.

Chemically Accurate Relative Folding Stability of RNA Hairpins from Molecular Simulations

To benchmark RNA force fields, we compared the folding stabilities of three 12-nucleotide hairpin stem loops estimated by simulation to stabilities determined by experiment. We used umbrella sampling and a reaction coordinate of end-to-end (5' to 3' hydroxyl oxygen) distance to estimate the free energy change of the transition from the native conformation to a fully extended conformation with no hydrogen bonds between non-neighboring bases. Each simulation was performed four times using the AMBER FF99+bsc0+χ_OL3 force field, and each window, spaced at 1 Å intervals, was sampled for 1 μs, for a total of 552 μs of simulation. We compared differences in the simulated free energy changes to analogous differences in free energies from optical melting experiments using thermodynamic cycles where the free energy change between stretched and random coil sequences is assumed to be sequence-independent. The differences between experimental and simulated ΔΔ G° are, on average, 0.98 ± 0.66 kcal/mol, which is chemically accurate and suggests that analogous simulations could be used predictively. We also report a novel method to identify where replica free energies diverge along a reaction coordinate, thus indicating where additional sampling would most improve convergence. We conclude by discussing methods to more economically perform these simulations.

Translation of non-standard codon nucleotides reveals minimal requirements for codon-anticodon interactions

The precise interplay between the mRNA codon and the tRNA anticodon is crucial for ensuring efficient and accurate translation by the ribosome. The insertion of RNA nucleobase derivatives in the mRNA allowed us to modulate the stability of the codon-anticodon interaction in the decoding site of bacterial and eukaryotic ribosomes, allowing an in-depth analysis of codon recognition. We found the hydrogen bond between the N1 of purines and the N3 of pyrimidines to be sufficient for decoding of the first two codon nucleotides, whereas adequate stacking between the RNA bases is critical at the wobble position. Inosine, found in eukaryotic mRNAs, is an important example of destabilization of the codon-anticodon interaction. Whereas single inosines are efficiently translated, multiple inosines, e.g., in the serotonin receptor 5-HT2C mRNA, inhibit translation. Thus, our results indicate that despite the robustness of the decoding process, its tolerance toward the weakening of codon-anticodon interactions is limited.

γ-Hemolysin Nanopore Is Sensitive to Guanine-to-Inosine Substitutions in Double-Stranded DNA at the Single-Molecule Level

Biological nanopores provide a unique single-molecule sensing platform to detect target molecules based on their specific electrical signatures. The γ-hemolysin (γ-HL) protein produced by Staphylococcus aureus is able to assemble into an octamer nanopore with a ∼2.3 nm diameter β-barrel. Herein, we demonstrate the first application of γ-HL nanopore for DNA structural analysis. To optimize conditions for ion-channel recording, the properties of the γ-HL pore (e.g., conductance, voltage-dependent gating, and ion-selectivity) were characterized at different pH, temperature, and electrolyte concentrations. The optimal condition for DNA analysis using γ-HL corresponds to 3 M KCl, pH 5, and T = 20 °C. The γ-HL protein nanopore is able to translocate dsDNA at about ∼20 bp/ms, and the unique current-signature of captured dsDNA can directly distinguish guanine-to-inosine substitutions at the single-molecule level with ∼99% accuracy. The slow dsDNA threading and translocation processes indicate this wild-type γ-HL channel has potential to detect other base modifications in dsDNA.

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

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

Tetrahelical structural family adopted by AGCGA-rich regulatory DNA regions

Here we describe AGCGA-quadruplexes, an unexpected addition to the well-known tetrahelical families, G-quadruplexes and i-motifs, that have been a focus of intense research due to their potential biological impact in G- and C-rich DNA regions, respectively. High-resolution structures determined by solution-state nuclear magnetic resonance (NMR) spectroscopy demonstrate that AGCGA-quadruplexes comprise four 5′-AGCGA-3′ tracts and are stabilized by G-A and G-C base pairs forming GAGA- and GCGC-quartets, respectively. Residues in the core of the structure are connected with edge-type loops. Sequences of alternating 5′-AGCGA-3′ and 5′-GGG-3′ repeats could be expected to form G-quadruplexes, but are shown herein to form AGCGA-quadruplexes instead. Unique structural features of AGCGA-quadruplexes together with lower sensitivity to cation and pH variation imply their potential biological relevance in regulatory regions of genes responsible for basic cellular processes that are related to neurological disorders, cancer and abnormalities in bone and cartilage development.

DNA tetrahelical structures such as G-quadruplexes are known to play important roles in DNA replication and repair. Here the authors present the structure of 5′-AGCGA-3′-quadruplexes enriched in genetic regulatory regions.

Synergy between NMR measurements and MD simulations of protein/RNA complexes: application to the RRMs, the most common RNA recognition motifs

RNA recognition motif (RRM) proteins represent an abundant class of proteins playing key roles in RNA biology. We present a joint atomistic molecular dynamics (MD) and experimental study of two RRM-containing proteins bound with their single-stranded target RNAs, namely the Fox-1 and SRSF1 complexes. The simulations are used in conjunction with NMR spectroscopy to interpret and expand the available structural data. We accumulate more than 50 μs of simulations and show that the MD method is robust enough to reliably describe the structural dynamics of the RRM-RNA complexes. The simulations predict unanticipated specific participation of Arg142 at the protein-RNA interface of the SRFS1 complex, which is subsequently confirmed by NMR and ITC measurements. Several segments of the protein-RNA interface may involve competition between dynamical local substates rather than firmly formed interactions, which is indirectly consistent with the primary NMR data. We demonstrate that the simulations can be used to interpret the NMR atomistic models and can provide qualified predictions. Finally, we propose a protocol for 'MD-adapted structure ensemble' as a way to integrate the simulation predictions and expand upon the deposited NMR structures. Unbiased μs-scale atomistic MD could become a technique routinely complementing the NMR measurements of protein-RNA complexes.

Predicting RNA Duplex Dimerization Free-Energy Changes upon Mutations Using Molecular Dynamics Simulations

The dimerization free energies of RNA-RNA duplexes are fundamental values that represent the structural stability of RNA complexes. We report a comparative analysis of RNA-RNA duplex dimerization free-energy changes upon mutations, estimated from a molecular dynamics simulation and experiments. A linear regression for nine pairs of double-stranded RNA sequences, six base pairs each, yielded a mean absolute deviation of 0.55 kcal/mol and an R(2) value of 0.97, indicating quantitative agreement between simulations and experimental data. The observed accuracy indicates that the molecular dynamics simulation with the current molecular force field is capable of estimating the thermodynamic properties of RNA molecules.

A Computational Model for Predicting Experimental RNA Nearest-Neighbor Free Energy Rankings: Inosine•Uridine Pairs

A computational model for predicting RNA nearest neighbor free energy rankings has been expanded to include the nonstandard nucleotide inosine. The model uses average fiber diffraction data and molecular dynamic simulations to generate input geometries for Quantum mechanic calculations. This resulted in calculated intrastrand stacking, interstrand stacking, and hydrogen bonding energies that were combined to give total binding energies. Total binding energies for RNA dimer duplexes containing inosine were ranked and compared to experimentally determined free energy ranks for RNA duplexes containing inosine. Statistical analysis showed significant agreement between the computationally determined ranks and the experimentally determined ranks.

Quantum chemical benchmark study on 46 RNA backbone families using a dinucleotide unit

We have created a benchmark set of quantum chemical structure-energy data denoted as UpU46, which consists of 46 uracil dinucleotides (UpU), representing all known 46 RNA backbone conformational families. Penalty-function-based restrained optimizations with COSMO TPSS-D3/def2-TZVP ensure a balance between keeping the target conformation and geometry relaxation. The backbone geometries are close to the clustering-means of their respective RNA bioinformatics family classification. High-level wave function methods (DLPNO-CCSD(T) as reference) and a wide-range of dispersion-corrected or inclusive DFT methods (DFT-D3, VV10, LC-BOP-LRD, M06-2X, M11, and more) are used to evaluate the conformational energies. The results are compared to the Amber RNA bsc0χOL3 force field. Most dispersion-corrected DFT methods surpass the Amber force field significantly in accuracy and yield mean absolute deviations (MADs) for relative conformational energies of ∼0.4-0.6 kcal/mol. Double-hybrid density functionals represent the most accurate class of density functionals. Low-cost quantum chemical methods such as PM6-D3H+, HF-3c, DFTB3-D3, as well as small basis set calculations corrected for basis set superposition errors (BSSEs) by the gCP procedure are also tested. Unfortunately, the presently available low-cost methods are struggling to describe the UpU conformational energies with satisfactory accuracy. The UpU46 benchmark is an ideal test for benchmarking and development of fast methods to describe nucleic acids, including force fields.

Mechanical properties of base-modified DNA are not strictly determined by base stacking or electrostatic interactions

This work probes the mystery of what balance of forces creates the extraordinary mechanical stiffness of DNA to bending and twisting. Here we explore the relationship between base stacking, functional group occupancy of the DNA minor and major grooves, and DNA mechanical properties. We study double-helical DNA molecules substituting either inosine for guanosine or 2,6-diaminopurine for adenine. These DNA variants, respectively, remove or add an amino group from the DNA minor groove, with corresponding changes in hydrogen-bonding and base stacking energy. Using the techniques of ligase-catalyzed cyclization kinetics, atomic force microscopy, and force spectroscopy with optical tweezers, we show that these DNA variants have bending persistence lengths within the range of values reported for sequence-dependent variation of the natural DNA bases. Comparison with seven additional DNA variants that modify the DNA major groove reveals that DNA bending stiffness is not correlated with base stacking energy or groove occupancy. Data from circular dichroism spectroscopy indicate that base analog substitution can alter DNA helical geometry, suggesting a complex relationship among base stacking, groove occupancy, helical structure, and DNA bend stiffness.

Guanine to inosine substitution leads to large increases in the population of a transient G·C Hoogsteen base pair

We recently showed that Watson–Crick base pairs in canonical duplex DNA exist in dynamic equilibrium with G(syn)·C⁺ and A(syn)·T Hoogsteen base pairs that have minute populations of ∼1%. Here, using nuclear magnetic resonance R_1ρ relaxation dispersion, we show that substitution of guanine with the naturally occurring base inosine results in an ∼17-fold increase in the population of transient Hoogsteen base pairs, which can be rationalized by the loss of a Watson–Crick hydrogen bond. These results provide further support for transient Hoogsteen base pairs and demonstrate that their population can increase significantly upon damage or chemical modification of the base.

Comparative genomics of nucleotide metabolism: a tour to the past of the three cellular domains of life

BACKGROUND

Nucleotide metabolism is central to all biological systems, due to their essential role in genetic information and energy transfer, which in turn suggests its possible presence in the last common ancestor (LCA) of Bacteria, Archaea and Eukarya. In this context, elucidation of the contribution of the origin and diversification of de novo and salvage pathways of nucleotide metabolism will allow us to understand the links between the enzymatic steps associated with the LCA and the emergence of the first metabolic pathways.

RESULTS

In this work, the taxonomical distribution of the enzymes associated with nucleotide metabolism was evaluated in 1,606 complete genomes. 151 sequence profiles associated with 120 enzymatic reactions were used. The evaluation was based on profile comparisons, using RPS-Blast. Organisms were clustered based on their taxonomical classifications, in order to obtain a normalized measure of the taxonomical distribution of enzymes according to the average of presence/absence of enzymes per genus, which in turn was used for the second step, to calculate the average presence/absence of enzymes per Clade.

CONCLUSION

From these analyses, it was suggested that divergence at the enzymatic level correlates with environmental changes and related modifications of the cell wall and membranes that took place during cell evolution. Specifically, the divergence of the 5-(carboxyamino) imidazole ribonucleotide mutase to phosphoribosylaminoimidazole carboxylase could be related to the emergence of multicellularity in eukaryotic cells. In addition, segments of salvage and de novo pathways were probably complementary in the LCA to the synthesis of purines and pyrimidines. We also suggest that a large portion of the pathway to inosine 5'-monophosphate (IMP) in purines could have been involved in thiamine synthesis or its derivatives in early stages of cellular evolution, correlating with the fact that these molecules may have played an active role in the protein-RNA world. The analysis presented here provides general observations concerning the adaptation of the enzymatic steps in the early stages of the emergence of life and the LCA.

Roles of the amino group of purine bases in the thermodynamic stability of DNA base pairing

The energetic aspects of hydrogen-bonded base-pair interactions are important for the design of functional nucleotide analogs and for practical applications of oligonucleotides. The present study investigated the contribution of the 2-amino group of DNA purine bases to the thermodynamic stability of oligonucleotide duplexes under different salt and solvent conditions, using 2'-deoxyriboinosine (I) and 2'-deoxyribo-2,6-diaminopurine (D) as non-canonical nucleotides. The stability of DNA duplexes was changed by substitution of a single base pair in the following order: G•C > D•T ≈ I•C > A•T > G•T > I•T. The apparent stabilization energy due to the presence of the 2-amino group of G and D varied depending on the salt concentration, and decreased in the water-ethanol mixed solvent. The effects of salt concentration on the thermodynamics of DNA duplexes were found to be partially sequence-dependent, and the 2-amino group of the purine bases might have an influence on the binding of ions to DNA through the formation of a stable base-paired structure. Our results also showed that physiological salt conditions were energetically favorable for complementary base recognition, and conversely, low salt concentration media and ethanol-containing solvents were effective for low stringency oligonucleotide hybridization, in the context of conditions employed in this study.

Modified Amber Force Field Correctly Models the Conformational Preference for Tandem GA pairs in RNA

Molecular mechanics with all-atom models was used to understand the conformational preference of tandem guanine-adenine (GA) noncanonical pairs in RNA. These tandem GA pairs play important roles in determining stability, flexibility, and structural dynamics of RNA tertiary structures. Previous solution structures showed that these tandem GA pairs adopt either imino (cis Watson-Crick/Watson-Crick A-G) or sheared (trans Hoogsteen/sugar edge A-G) conformations depending on the sequence and orientation of the adjacent closing base pairs. The solution structures (GCGGACGC)2 [Biochemistry, 1996, 35, 9677-9689] and (GCGGAUGC)2 [Biochemistry, 2007, 46, 1511-1522] demonstrate imino and sheared conformations for the two central GA pairs, respectively. These systems were studied using molecular dynamics and free energy change calculations for conformational changes, using umbrella sampling. For the structures to maintain their native conformations during molecular dynamics simulations, a modification to the standard Amber ff10 force field was required, which allowed the amino group of guanine to leave the plane of the base [J. Chem. Theory Comput., 2009, 5, 2088-2100] and form out-of-plane hydrogen bonds with a cross-strand cytosine or uracil. The requirement for this modification suggests the importance of out-of-plane hydrogen bonds in stabilizing the native structures. Free energy change calculations for each sequence demonstrated the correct conformational preference when the force field modification was used, but the extent of the preference is underestimated.

Inosine octamer stabilized by alkali earth metal cations - as studied by electrospray ionization mass spectrometry

By using electrospray ionization mass spectrometry, inosine was found to be able to form an octamer stabilized by alkali earth metal cation, namely Ca(2+), Sr(2+) and Ba(2+), of which the most stable is that stabilized by Ca(2+) (ion [I8+Ca](2+)). It was established that 9-methylhypoxanthine (M) did not form an analogical octamer, since ion [M8+Ca](2+) was not detected. On the other hand, 9-methylhypoxanthine can form "mixed" octamers together with inosine (ions [InMm+Ca](2+), n + m = 8, were detected).

Higher order structural effects stabilizing the reverse Watson-Crick Guanine-Cytosine base pair in functional RNAs

The G:C reverse Watson-Crick (W:W trans) base pair, also known as Levitt base pair in the context of tRNAs, is a structurally and functionally important base pair that contributes to tertiary interactions joining distant domains in functional RNA molecules and also participates in metabolite binding in riboswitches. We previously indicated that the isolated G:C W:W trans base pair is a rather unstable geometry, and that dicationic metal binding to the Guanine base or posttranscriptional modification of the Guanine can increase its stability. Herein, we extend our survey and report on other H-bonding interactions that can increase the stability of this base pair. To this aim, we performed a bioinformatics search of the PDB to locate all the occurencies of G:C trans base pairs. Interestingly, 66% of the G:C trans base pairs in the PDB are engaged in additional H-bonding interactions with other bases, the RNA backbone or structured water molecules. High level quantum mechanical calculations on a data set of representative crystal structures were performed to shed light on the structural stability and energetics of the various crystallographic motifs. This analysis was extended to the binding of the preQ1 metabolite to a preQ1-II riboswitch.

Propensities for loop structures of RNA & DNA backbones

RNA oligonucleotides exhibit a large tendency to bend and form a loop conformation which is a major motif contributing to their complex three-dimensional structure. This is in contrast to DNA molecules that predominantly form the double-helix structure. In this paper we investigate by molecular dynamics simulation, as well as, by its combination with the replica-exchange method, the propensity of RNA chains containing the GCUAA pentaloop to form spontaneously a hairpin conformation. The results were then compared with those of analogous hybrid oligonucleotides in which the ribose groups in the loop-region were substituted by deoxyriboses. We find that the RNA oligomers exhibit a marginal excess stability to form loop structures. The equilibrium constant for opening the loop to an extended conformation is twice as large in the hybrid than it is in the RNA chain. Analyses of the hydrogen bonds indicate that the excess stability for forming a hairpin is a result of hydrogen bonds the 2'-hydroxyls in the loop region form with other groups in the loop. Of these hydrogen bonds, the most important is the hydrogen bond donated from the 2'-OH at the first position of the loop to N7 of adenine at the forth position. RNA and DNA backbones are characterized by different backbone dihedral angles and sugar puckering that can potentially facilitate or hamper the hydrogen bonds involving the 2'-OH. Nevertheless, the sugar puckerings of all the pentaloop nucleotides were not significantly different between the two chains displaying the C3'-endo conformation characteristic to the A-form double helix. All of the other backbone dihedrals also did not show any considerable difference in the loop-region except of the δ-dihedral. In this case, the RNA loop exhibited bimodal distributions corresponding to, both, the RNA and DNA backbones, whereas the loop of the hybrid chain behaved mostly as that of a DNA backbone. Thus, it is possible that the behavior of the δ-dihedrals in the loop-region of the RNA adopts conformations that facilitate the intra-nucleotide hydrogen bondings of the 2'-hydroxyls, and consequently renders loop structures in RNA more stable.

Ultrafast electronic deactivation dynamics of the inosine dimer--a model case for H-bonded purine bases

The structural properties and ultrafast electronic deactivation dynamics of the inosine dimer in CHCl3 have been investigated by two-dimensional (1)H NMR and static FTIR spectroscopy and by femtosecond time-resolved transient absorption spectroscopy, respectively. The (1)H NMR and IR spectra show the formation of a well-defined, symmetric dimer with an association equilibrium constant of KI·I = 690 ± 100 M(-1). The excited-state dynamics after photoexcitation at λpump = 260 nm monitored by ultrafast absorption spectroscopy show great similarity with those of the monomer inosine in an aqueous solution and are governed by a decay time of τ = 90 ± 10 fs, which is one of the shortest electronic lifetimes of all nucleobases and nucleobase dimers studied so far. On the basis of these observations, the inosine dimer is expected to follow a similar relaxation pathway as the monomer, involving an out-of-plane deformation of the six-membered ring. The importance of the C(2) position for the electronic deactivation of hypoxanthine and guanine is discussed. The obtained well-determined structure and straightforward dynamics qualify the inosine dimer as an excellent reference case for more complicated systems such as the G·G dimer and the G·C and A·T Watson-Crick pairs.

How to understand quantum chemical computations on DNA and RNA systems? A practical guide for non-specialists

In this review primarily written for non-experts we explain basic methodological aspects and interpretation of modern quantum chemical (QM) computations applied to nucleic acids. We introduce current reference QM computations on small model systems consisting of dozens of atoms. Then we comment on recent advance of fast and accurate dispersion-corrected density functional theory methods, which will allow computations of small but complete nucleic acids building blocks in the near future. The qualitative difference between QM and molecular mechanics (MM, force field) computations is discussed. We also explain relation of QM and molecular simulation computations to experiments.