Molecular Dynamics of Double Stranded Xylo-Nucleic Acid

Structural insights into the morpholino nucleic acid/RNA duplex using the new XNA builder Ducque in a molecular modeling pipeline

The field of synthetic nucleic acids with novel backbone structures [xenobiotic nucleic acids (XNAs)] has flourished due to the increased importance of XNA antisense oligonucleotides and aptamers in medicine, as well as the development of XNA processing enzymes and new XNA genetic materials. Molecular modeling on XNA structures can accelerate rational design in the field of XNAs as it contributes in understanding and predicting how changes in the sugar-phosphate backbone impact on the complementation properties of the nucleic acids. To support the development of novel XNA polymers, we present a first-in-class open-source program (Ducque) to build duplexes of nucleic acid analogs with customizable chemistry. A detailed procedure is described to extend the Ducque library with new user-defined XNA fragments using quantum mechanics (QM) and to generate QM-based force field parameters for molecular dynamics simulations within standard packages such as AMBER. The tool was used within a molecular modeling workflow to accurately reproduce a selection of experimental structures for nucleic acid duplexes with ribose-based as well as non-ribose-based nucleosides. Additionally, it was challenged to build duplexes of morpholino nucleic acids bound to complementary RNA sequences.

Efficient and Accurate Potential Energy Surfaces of Puckering in Sugar-Modified Nucleosides

Puckering of the sugar unit in nucleosides and nucleotides is an important structural aspect that directly influences the helical structure of nucleic acids. The preference for specific puckering modes in nucleic acids can be analyzed via in silico conformational analysis, but the large amount of conformations and the accuracy of the analysis leads to an extensive amount of computational time. In this paper, we show that the combination of geometry optimizations with the HF-3c method with single point energies at the RI-MP2 level results in accurate results for the puckering potential energy surface (PES) of DNA and RNA nucleosides while significantly reducing the necessary computational time. Applying this method to a series of known xeno nucleic acids (XNAs) allowed us to rapidly explore the puckering PES of each of the respective nucleosides and to explore the puckering PES of six-membered modified XNA (HNA and β-homo-DNA) for the first time.

Stable Hairpin Structures Formed by Xylose-Based Nucleic Acids

Xenobiology explores synthetic nucleic acid polymers as alternative carriers of genetic information to expand the central dogma. The xylo- and deoxyxylo-nucleic acids (XyNA and dXyNA), containing 3' epimers of riboses and deoxyriboses, are considered to be potential candidates for an orthogonal system. In this study, thermal and spectroscopic analyses show that XyNA and dXyNA form stable hairpins. The dXyNA hairpin structure determined by NMR spectroscopy contains a flexible loop that locks the stem into a stable ladder-like duplex with marginal right-handed helicity. The reduced flexibility of the dXyNA duplex observed in the stem of the hairpin demonstrates that folding of dXyNA yields more stable structures described so far.

The proto-Nucleic Acid Builder: a software tool for constructing nucleic acid analogs

The helical structures of DNA and RNA were originally revealed by experimental data. Likewise, the development of programs for modeling these natural polymers was guided by known structures. These nucleic acid polymers represent only two members of a potentially vast class of polymers with similar structural features, but that differ from DNA and RNA in the backbone or nucleobases. Xeno nucleic acids (XNAs) incorporate alternative backbones that affect the conformational, chemical, and thermodynamic properties of XNAs. Given the vast chemical space of possible XNAs, computational modeling of alternative nucleic acids can accelerate the search for plausible nucleic acid analogs and guide their rational design. Additionally, a tool for the modeling of nucleic acids could help reveal what nucleic acid polymers may have existed before RNA in the early evolution of life. To aid the development of novel XNA polymers and the search for possible pre-RNA candidates, this article presents the proto-Nucleic Acid Builder (https://github.com/GT-NucleicAcids/pnab), an open-source program for modeling nucleic acid analogs with alternative backbones and nucleobases. The torsion-driven conformation search procedure implemented here predicts structures with good accuracy compared to experimental structures, and correctly demonstrates the correlation between the helical structure and the backbone conformation in DNA and RNA.

Energy Landscapes of Deoxyxylo- and Xylo-Nucleic Acid Octamers

Artificial analogues of the natural nucleic acids have attracted interest as a diverse class of information storage molecules capable of self-replication. In this study, we use the computational potential energy landscape framework to investigate the structural and dynamical properties of xylo- and deoxyxylo-nucleic acids (XyNA and dXyNA), which are derived from their respective RNA and DNA analogues by inversion of a single chiral center in the sugar moiety of the nucleotides. For an octameric XyNA sequence and the analogue dXyNA, we observe facile conformational transitions between a left-handed helix, which is the free energy global minimum, and a ladder-type structure with approximately zero helicity. The competing ensembles are better separated in the dXyNA, making it a more suitable candidate for a molecular switch, whereas the XyNA exhibits additional flexibility. Both energy landscapes exhibit greater frustration than we observe in RNA or DNA, in agreement with the higher degree of optimization expected from the principle of minimal frustration in evolved biomolecules.

The furanosidic scaffold of d-ribose: a milestone for cell life

The recruitment of the furanosidic scaffold of ribose as the crucial step for nucleotides and then for nucleic acids synthesis is presented. Based on the view that the selection of molecules to be used for relevant metabolic purposes must favor structurally well-defined molecules, the inadequacy of ribose as a preferential precursor for nucleotides synthesis is discussed. The low reliability of ribose in its furanosidic hemiacetal form must have played ab initio against the choice of d-ribose for the generation of d-ribose-5-phosphate, the fundamental precursor of the ribose moiety of nucleotides. The latter, which is instead generated through the 'pentose phosphate pathway' is strictly linked to the affordable and reliable pyranosidic structure of d-glucose.

Probing the role of dispersion energy on structural transformation of double-stranded xylo- and ribo-nucleic acids

The structural transformation of double-stranded octameric xyloNA and RNA were probed by modulating the dispersion energy. For the RNA, the increase and the decrease in dispersion energy lead to over-winding and unwinding of the helix. These structural transformations resemble the features observed due to the action of the topoisomerases and helicases enzymes, respectively. On the other hand, an increase in the dispersion energy has minimal effect on the structural transformation of double-strand xyloNA, whilst a decrease in the dispersion energy results in a structural transformation which happens due to the action of the helicases. The unresponsive behaviour of xyloNA to an increase in the dispersion energy is attributed to the presence of an Lpπ interaction between the oxygen atom of the xylose sugar and the adjacent nucleobase.

Intradomain Interactions in an NMDA Receptor Fragment Mediate N-Glycan Processing and Conformational Sampling

The structural and functional roles of highly conserved asparagine-linked (N)-glycans on the extracellular ligand-binding domain (LBD) of the N-methyl-D-aspartate receptors are poorly understood. We applied solution- and computation-based methods that identified N-glycan-mediated intradomain and interglycan interactions. Nuclear magnetic resonance (NMR) spectra of the GluN1 LBD showed clear signals corresponding to each of the three N-glycans and indicated the reducing end of glycans at N440 and N771 potentially contacted nearby amino acids. Molecular dynamics simulations identified contacts between nearby amino acids and the N440- and N771-glycans that were consistent with the NMR spectra. The distal portions of the N771-glycan also contacted the core residues of the nearby N471-glycan. This result was consistent with mass spectrometry data indicating the limited N471-glycan core fucosylation and reduced branch processing of the N771-glycan could be explained by interglycan contacts. We discuss a potential role for the GluN1 LBD N-glycans in interdomain contacts formed in NMDA receptors.