Quantum Chemical Studies Of Nucleic Acids Can We Construct A Bridge To The Rna Structural Biology And Bioinformatics Communities?

In this feature article we provide a side-by-side introduction for two research fields quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics Our main aim is to demonstrate that these research areas while largely separated in contemporary literature have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution

Noncanonical Hydrogen Bonding In Nucleic Acids. Benchmark Evaluation Of Key Base-phosphate Interactions In Folded Rna Molecules Using Quantum-chemical Calculations And Molecular Dynamics Simulations

RNA molecules are stabilized by a wide range of non canonical interactions that are not present in DNA. Among them, the recently classified base phosphate (BPh) interactions belong to the most important ones. Twelve percent of nucleotides in the ribosomal crystal structures are involved in BPh interactions. BPh interactions are highly conserved and provide major constraints on RNA sequence evolution. Here we provide assessment of the energetics of BPh interactions using MP2 computations extrapolated to the complete basis set of atomic orbitals and corrected for higher-order electron correlation effects. The reference computations are compared with DFT-D and DFT-D3 approaches, the SAPT method, and the molecular mechanics force field. The computations, besides providing the basic benchmark for the BPh interactions, allow some refinements of the original classification, including identification of some potential doubly bonded BPh patterns. The reference computations are followed by analysis of some larger RNA fragments that consider the context of the BPh interactions. The computations demonstrate the complexity of interaction patterns utilizing the BPh interactions in real RNA structures. The BPh interactions are often involved in intricate interaction networks. We studied BPh interactions of protonated adenine that can contribute to catalysis of hairpin ribozyme, the key BPh interaction in the S-turn motif of the sarcin ricin loop, which may predetermine the S-turn topology and complex BPh patterns-from the glmS riboswitch. Finally, the structural stability of BPh interactions in explicit solvent molecular dynamics simulations is assessed. The simulations well preserve key BPh interactions and allow dissection of structurally/functionally important water-meditated BPh bridges, which could not be considered in earlier bioinformatics classification of BPh interactions

Exploring the Dynamics of Propeller Loops in Human Telomeric DNA Quadruplexes Using Atomistic Simulations

We have carried out a series of extended unbiased molecular dynamics (MD) simulations (up to 10 \u3bcs long, 3c162 \u3bcs in total) complemented by replica-exchange with the collective variable tempering (RECT) approach for several human telomeric DNA G-quadruplex (GQ) topologies with TTA propeller loops. We used different AMBER DNA force-field variants and also processed simulations by Markov State Model (MSM) analysis. The slow conformational transitions in the propeller loops took place on a scale of a few \u3bcs, emphasizing the need for long simulations in studies of GQ dynamics. The propeller loops sampled similar ensembles for all GQ topologies and for all force-field dihedral-potential variants. The outcomes of standard and RECT simulations were consistent and captured similar spectrum of loop conformations. However, the most common crystallographic loop conformation was very unstable with all force-field versions. Although the loss of canonical \u3b3-trans state of the first propeller loop nucleotide could be related to the indispensable bsc0 \u3b1/\u3b3 dihedral potential, even supporting this particular dihedral by a bias was insufficient to populate the experimentally dominant loop conformation. In conclusion, while our simulations were capable of providing a reasonable albeit not converged sampling of the TTA propeller loop conformational space, the force-field description still remained far from satisfactory

Revisiting the planarity of nucleic acid bases: Pyramidilization at glycosidic nitrogen in purine bases is modulated by orientation of glycosidic torsion

We describe a novel, fundamental property of nucleobase structure, namely, pyramidilization at the N1/9 sites of purine and pyrimidine bases. Through a combined analyses of ultra-high-resolution X-ray structures of both oligonucleotides extracted from the Nucleic Acid Database and isolated nucleotides and nucleosides from the Cambridge Structural Database, together with a series of quantum chemical calculations, molecular dynamics (MD) simulations, and published solution nuclear magnetic resonance (NMR) data, we show that pyramidilization at the glycosidic nitrogen is an intrinsic property. This property is common to isolated nucleosides and nucleotides as well as oligonucleotides—it is also common to both RNA and DNA. Our analysis suggests that pyramidilization at N1/9 sites depends in a systematic way on the local structure of the nucleoside. Of note, the pyramidilization undergoes stereo-inversion upon reorientation of the glycosidic bond. The extent of the pyramidilization is further modulated by the conformation of the sugar ring. The observed pyramidilization is more pronounced for purine bases, while for pyrimidines it is negligible. We discuss how the assumption of nucleic acid base planarity can lead to systematic errors in determining the conformation of nucleotides from experimental data and from unconstrained MD simulations

A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein–DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations. This article describes results providing information not only on all 10 unique base pair steps, but also on all possible nearest-neighbor effects on these steps. These results are derived from simulations of 50–100 ns on 39 different DNA oligomers in explicit solvent and using a physiological salt concentration. We demonstrate that the simulations are converged in terms of helical and backbone parameters. The results show that nearest-neighbor effects on base pair steps are very significant, implying that dinucleotide models are insufficient for predicting sequence-dependent behavior. Flanking base sequences can notably lead to base pair step parameters in dynamic equilibrium between two conformational sub-states. Although this study only provides limited data on next-nearest-neighbor effects, we suggest that such effects should be analyzed before attempting to predict the sequence-dependent behavior of DNA

Cation binding to 15-TBA quadruplex DNA is a multiple-pathway cation-dependent process

A combination of explicit solvent molecular dynamics simulation (30 simulations reaching 4 µs in total), hybrid quantum mechanics/molecular mechanics approach and isothermal titration calorimetry was used to investigate the atomistic picture of ion binding to 15-mer thrombin-binding quadruplex DNA (G-DNA) aptamer. Binding of ions to G-DNA is complex multiple pathway process, which is strongly affected by the type of the cation. The individual ion-binding events are substantially modulated by the connecting loops of the aptamer, which play several roles. They stabilize the molecule during time periods when the bound ions are not present, they modulate the route of the ion into the stem and they also stabilize the internal ions by closing the gates through which the ions enter the quadruplex. Using our extensive simulations, we for the first time observed full spontaneous exchange of internal cation between quadruplex molecule and bulk solvent at atomistic resolution. The simulation suggests that expulsion of the internally bound ion is correlated with initial binding of the incoming ion. The incoming ion then readily replaces the bound ion while minimizing any destabilization of the solute molecule during the exchange

Selectivity of major isoquinoline alkaloids from Chelidonium majus towards telomeric G-quadruplex: A study using a transition-FRET (t-FRET) assay

Background Natural bioproducts are invaluable resources in drug discovery. Isoquinoline alkaloids of Chelidonium majus constitute a structurally diverse family of natural products that are of great interest, one of them being their selectivity for human telomeric G-quadruplex structure and telomerase inhibition. Methods The study focuses on the mechanism of telomerase inhibition by stabilization of telomeric G-quadruplex structures by berberine, chelerythrine, chelidonine, sanguinarine and papaverine. Telomerase activity and mRNA levels of hTERT were estimated using quantitative telomere repeat amplification protocol (q-TRAP) and qPCR, in MCF-7 cells treated with different groups of alkaloids. The selectivity of the main isoquinoline alkaloids of Chelidonium majus towards telomeric G-quadruplex forming sequences were explored using a sensitive modified thermal FRET-melting measurement in the presence of the complementary oligonucleotide CT22. We assessed and monitored G-quadruplex topologies using circular dichroism (CD) methods, and compared spectra to previously well-characterized motifs, either alone or in the presence of the alkaloids. Molecular modeling was performed to rationalize ligand binding to the G-quadruplex structure. Results The results highlight strong inhibitory effects of chelerythrine, sanguinarine and berberine on telomerase activity, most likely through substrate sequestration. These isoquinoline alkaloids interacted strongly with telomeric sequence G-quadruplex. In comparison, chelidonine and papaverine had no significant interaction with the telomeric quadruplex, while they strongly inhibited telomerase at transcription level of hTERT. Altogether, all of the studied alkaloids showed various levels and mechanisms of telomerase inhibition. Conclusions We report on a comparative study of anti-telomerase activity of the isoquinoline alkaloids of Chelidonium majus. Chelerythrine was most effective in inhibiting telomerase activity by substrate sequesteration through G-quadruplex stabilization. General significance Understanding structural and molecular mechanisms of anti-cancer agents can help in developing new and more potent drugs with fewer side effects. Isoquinolines are the most biologically active agents from Chelidonium majus, which have shown to be telomeric G-quadruplex stabilizers and potent telomerase inhibitors

A Computational Quantum-Based Perspective on the Molecular Origins of Life’s Building Blocks

The search for the chemical origins of life represents a long-standing and continuously debated enigma. Despite its exceptional complexity, in the last decades the field has experienced a revival, also owing to the exponential growth of the computing power allowing for efficiently simulating the behavior of matter—including its quantum nature—under disparate conditions found, e.g., on the primordial Earth and on Earth-like planetary systems (i.e., exoplanets). In this minireview, we focus on some advanced computational methods capable of efficiently solving the Schrödinger equation at different levels of approximation (i.e., density functional theory)—such as ab initio molecular dynamics—and which are capable to realistically simulate the behavior of matter under the action of energy sources available in prebiotic contexts. In addition, recently developed metadynamics methods coupled with first-principles simulations are here reviewed and exploited to answer to old enigmas and to propose novel scenarios in the exponentially growing research field embedding the study of the chemical origins of life