Using NMR and molecular dynamics to link structure and dynamics effects of the universal base 8-aza, 7-deaza, N8 linked adenosine analog

MIM-ML: A Novel Quantum Chemical Fragment-Based Random Forest Model for Accurate Prediction of NMR Chemical Shifts of Nucleic Acids

We developed a random forest machine learning (ML) model for the prediction of 1H and 13C NMR chemical shifts of nucleic acids. Our ML model is trained entirely on reproducing computed chemical shifts obtained previously on 10 nucleic acids using a Molecules-in-Molecules (MIM) fragment-based density functional theory (DFT) protocol including microsolvation effects. Our ML model includes structural descriptors as well as electronic descriptors from an inexpensive low-level semiempirical calculation (GFN2-xTB) and trained on a relatively small number of DFT chemical shifts (2080 1H chemical shifts and 1780 13C chemical shifts on the 10 nucleic acids). The ML model is then used to make chemical shift predictions on 8 new nucleic acids ranging in size from 600 to 900 atoms and compared directly to experimental data. Though no experimental data was used in the training, the performance of our model is excellent (mean absolute deviation of 0.34 ppm for 1H chemical shifts and 2.52 ppm for 13C chemical shifts for the test set), despite having some nonstandard structures. A simple analysis suggests that both structural and electronic descriptors are critical for achieving reliable predictions. This is the first attempt to combine ML from fragment-based DFT calculations to predict experimental chemical shifts accurately, making the MIM-ML model a valuable tool for NMR predictions of nucleic acids.

Accurate and Cost-Effective NMR Chemical Shift Predictions for Nucleic Acids Using a Molecules-in-Molecules Fragmentation-Based Method

We have developed, implemented, and assessed an efficient protocol for the prediction of NMR chemical shifts of large nucleic acids using our molecules-in-molecules (MIM) fragment-based quantum chemical approach. To assess the performance of our approach, MIM-NMR calculations are calibrated on a test set of three nucleic acids, where the structure is derived from solution-phase NMR studies. For DNA systems with multiple conformers, the one-layer MIM method with trimer fragments (MIM1trimer) is benchmarked to get the lowest energy structure, with an average error of only 0.80 kcal/mol with respect to unfragmented full molecule calculations. The MIMI-NMRdimer calibration with respect to unfragmented full molecule calculations shows a mean absolute deviation (MAD) of 0.06 and 0.11 ppm, respectively, for 1H and 13C nuclei, but the performance with respect to experimental NMR chemical shifts is comparable to the more expensive MIM1-NMR and MIM2-NMR methods with trimer subsystems. To compare with the experimental chemical shifts, a standard protocol is derived using DNA systems with Protein Data Bank (PDB) IDs 1SY8, 1K2K, and 1KR8. The effect of structural minimizations is employed using a hybrid mechanics/semiempirical approach and used for computations in solution with implicit and explicit-implicit solvation models in our MIM1-NMRdimer methodology. To demonstrate the applicability of our protocol, we tested it on seven nucleic acids, including structures with nonstandard residues, heteroatom substitutions (F and B atoms), and side chain mutations with a size ranging from ∼300 to 1100 atoms. The major improvement for predicted MIM1-NMRdimer calculations is obtained from structural minimizations and implicit solvation effects. A significant improvement with the explicit-implicit solvation model is observed only for two smaller nucleic acid systems (1KR8 and 7NBK), where the expensive first solvation shell is replaced by the microsolvation model, in which a single water molecule is added for each solvent-exposed amino and imino protons, along with the implicit solvation. Overall, our target accuracy of ∼0.2-0.3 ppm for 1H and ∼2-3 ppm for 13C has been achieved for large nucleic acids. The proposed MIM-NMR approach is accurate and cost-effective (linear scaling with system size), and it can aid in the structural assignments of a wide range of complex biomolecules.

Isoguanine (2-Hydroxyadenine) and 2-Aminoadenine Nucleosides with an 8-Aza-7-deazapurine Skeleton: Synthesis, Functionalization with Fluorescent and Clickable Side Chains, and Impact of 7-Substituents on Physical Properties

7-Functionalized 8-aza-7-deaza-2'-deoxyisoguanine and 8-aza-7-deaza-2-aminoadenine 2'-deoxyribonucleosides decorated with fluorescent pyrene or benzofuran sensor tags or clickable side chains with terminal triple bonds were synthesized. 8-Aza-7-deaza-7-iodo-2-amino-2'-deoxyadenosine was used as the central intermediate and was accessible by an improved two-step glycosylation/amination protocol. Functionalization of position-7 was performed either on 8-aza-7-deaza-7-iodo-2-amino-2'-deoxyadenosine followed by selective deamination of the 2-amino group or on 7-iodinated 8-aza-7-deaza-2'-deoxyisoguanosine. Sonogashira and Suzuki-Miyaura cross-coupling reactions were employed for this purpose. Octadiynyl side chains were selected as linkers for click reactions with azido pyrenes. K_Taut values calculated from H₂O/dioxane mixtures revealed that side chains have a significant influence on the tautomeric equilibrium. Photophysical properties (fluorescence, solvatochromism, and quantum yields) of the new 8-aza-7-deazapurine nucleosides with fluorescent side chains were determined. Remarkably, a strong excimer fluorescence in H₂O was observed for pyrene dye conjugates of 8-aza-7-deazaisoguanine and 2-aminoadenine nucleosides with a long linker. In other solvents including methanol, excimer fluorescence was negligible. The 2-aminoadenine and isoguanine nucleosides with the 8-aza-7-deazapurine skeleton expand the class of nucleosides applicable to fluorescence detection with respect to diagnostic and therapeutic purposes.

NMR Structure Determination for Oligonucleotides

NMR spectroscopy is a versatile tool for determining the structure and dynamics of nucleic acids under solution conditions. In this unit, we provide an overview and detail of the experiments and methods used in our laboratory to determine the structure of oligonucleotides at natural abundance, thus limiting our approach to ¹ H, ¹³ C, and ³¹ P NMR techniques. Isotopic labeling is heavily used in RNA NMR studies, however, labeling of DNA is still less common and, if modified nucleotides are investigated, is exceptionally expensive or not feasible. Each method described here is extensively documented and annotated with tips and observations to facilitate their application. Sections are devoted to sample preparation, NMR experiments and setup, resonance assignment, structure generation protocols, evaluation, tips that may be useful, and software sources. © 2018 by John Wiley & Sons, Inc.

Molecular dynamics simulations and applications in computational toxicology and nanotoxicology

Nanotoxicology studies toxicity of nanomaterials and has been widely applied in biomedical researches to explore toxicity of various biological systems. Investigating biological systems through in vivo and in vitro methods is expensive and time taking. Therefore, computational toxicology, a multi-discipline field that utilizes computational power and algorithms to examine toxicology of biological systems, has gained attractions to scientists. Molecular dynamics (MD) simulations of biomolecules such as proteins and DNA are popular for understanding of interactions between biological systems and chemicals in computational toxicology. In this paper, we review MD simulation methods, protocol for running MD simulations and their applications in studies of toxicity and nanotechnology. We also briefly summarize some popular software tools for execution of MD simulations.

Silver-Mediated Base Pairs in DNA Incorporating Purines, 7-Deazapurines, and 8-Aza-7-deazapurines: Impact of Reduced Nucleobase Binding Sites and an Altered Glycosylation Position

Formation of silver-mediated DNA was studied with oligonucleotides incorporating 8-aza-7-deazapurine, 7-deazapurine, and purine nucleosides. The investigation was performed on non-self-complementary duplexes with one or two modifications and self-complementary duplexes with an alternating dA-dT motif. Homo base pairs as well as base pair mismatches of dA analogues with dC and Watson-Crick pairs with dT were studied by stoichiometric silver ion titration and T_m measurements. N⁸ -Glycosylated 8-aza-7-deazaadenine forms silver-ion-mediated base pairs capturing two silver ions (low silver content) whereas regularly glycosylated 8-aza-7-deazapurine, 7-deazapurine (c⁷ A_d ), and dA do not form comparable structures. Stable silver-mediated "dA-dC" base pair mismatches were detected for all nucleosides. Two silver ions per base pair are bound by 8-aza-7-deazapurine whereas c⁷ A_d binds only one silver ion. The situation is different when the equivalents of silver ions were increased to the number of total base pairs. Surprisingly, in 12-mer duplexes as well as in related 25-mer duplexes every base pair consumed one silver ion.

Structural Impact of Single Ribonucleotide Residues in DNA

Single ribonucleotide intrusions represent the most common nonstandard nucleotide type found incorporated in genomic DNA, yet little is known of their structural impact. This lesion incurs genomic instability in addition to affecting the physical properties of the DNA. To probe for structural and dynamic effects of single ribonucleotides in various sequence contexts-AxC, CxG, and GxC, where x=rG or dG-we report the structures of three single-ribonucleotide-containing DNA duplexes and the corresponding DNA controls. The lesion subtly and locally perturbs the structure asymmetrically on the 3' side of the lesion in both the riboguanosine-containing and the complementary strand of the duplex. The perturbations are mainly restricted to the sugar and phosphodiester backbone. The ribose and 3'-downstream deoxyribose units are predominately in N-type conformation; backbone torsion angles ϵ and/or ζ of the ribonucleotide or upstream deoxyribonucleotide are affected. Depending on the flanking sequences, the C2'-OH group forms hydrogen bonds with the backbone, 3'-neighboring base, and/or sugar. Interestingly, even in similar purine-rG-pyrimidine environments (A-rG-C and G-rG-C), a riboguanosine unit affects DNA in a distinct manner and manifests different hydrogen bonds, which makes generalizations difficult.