Watson-crick base pairs with thiocarbonyl groups: How sulfur changes the hydrogen bonds in DNA

Ab initio Modeling of Hydrogen Bonding of Remdesivir and Adenosine with Uridine

Remdesivir (RDV) emerged as an effective drug against the SARS-CoV-2 virus pandemic. One of the crucial steps in the mechanism of action of RDV is its incorporation into the growing RNA strand. RDV, an adenosine analogue, forms Watson-Crick (WC) type hydrogen bonds with uridine in the complementary strand and the strength of this interaction will control efficacy of RDV. While there is a plethora of structural and energetic information available about WC H-bonds in natural base pairs, the interaction of RDV with uridine has not been studied yet at the atomic level. In this article, we aim to bridge this gap, to understand RDV and its hydrogen bonding interactions, by employing density functional theory (DFT) at the M06-2X/cc-pVDZ level. The interaction energy, QTAIM analysis, NBO and SAPT2 are performed for RDV, adenosine, and their complex with uridine to gain insights into the nature of hydrogen bonding. The computations show that RDV has similar geometry, energetic, molecular orbitals, and aromaticity as adenosine, suggesting that RDV is an effective adenosine analogue. The important geometrical parameters, such as bond distances and red-shift in the stretching vibrational modes of adenosine, RDV and uridine identify two WC-type H-bonds. The relative strength of these two H-bonds is computed using QTAIM parameters and the computed hydrogen bond energy. Finally, the SAPT2 study is performed at the minima and at non-equilibrium base pair distances to understand the dominant intermolecular physical force. This study, based on a thorough analysis of a variety of computations, suggests that both adenosine and RDV have similar structure, energetic, and hydrogen bonding behaviour.

Sulfur-substitution-induced base flipping in the DNA duplex

Base flipping is widely observed in a number of important biological processes. The genetic codes deposited inside the DNA duplex become accessible to external agents upon base flipping. The sulfur substitution of guanine leads to thioguanine, which alters the thermodynamic stability of the GC base pairs and the GT mismatches. Experimental studies conclude that the sulfur substitution decreases the lifetime of the GC base pair. In this work, under three AMBER force fields for nucleotide systems, we firstly performed equilibrium and nonequilibrium free energy simulations to investigate the variation of the thermodynamic profiles in base flipping upon sulfur substitution. It is found that the bsc0 modification, the bsc1 modification and the OL15 modification of AMBER force fields are able to qualitatively describe the sulfur-substitution dependent behavior of the thermodynamics. However, only the two last-generation AMBER force fields are able to provide quantitatively correct predictions. The second computational study on the sulfur substitutions focused on the relative stability of the S6G-C base pair and the S6G-T mismatch. Two conflicting experimental observations were reported by the same authors. One suggested that the S6G-C base pair was more stable, while the other concludes that the S6G-T mismatch was more stable. We answered this question by constructing the free energy profiles along the base flipping pathway computationally.

tRNA structural and functional changes induced by oxidative stress

Oxidatively damaged biomolecules impair cellular functions and contribute to the pathology of a variety of diseases. RNA is also attacked by reactive oxygen species, and oxidized RNA is increasingly recognized as an important contributor to neurodegenerative complications in humans. Recently, evidence has accumulated supporting the notion that tRNA is involved in cellular responses to various stress conditions. This review focuses on the intriguing consequences of oxidative modification of tRNA at the structural and functional level.

Desulfurization of 2-thiouracil nucleosides: conformational studies of 4-pyrimidinone nucleosides

4-Pyrimidinone ribofuranoside (H(2)o(4)U) and 4-pyrimidinone 2'-deoxyribofuranoside (dH(2)o(4)U) were synthesized by the oxidative desulfurization of parent 2-thiouracil nucleosides with m-chloroperbenzoic acid. The crystal structures of H(2)o(4)U and dH(2)o(4)U and their conformations in solution were determined and compared with corresponding 2-thiouracil and uracil nucleosides. The absence of a large 2-thiocarbonyl/2-carbonyl group in the nucleobase moiety results in C2'-endo puckering of the ribofuranose ring (S conformer) in the crystal structure of H(2)o(4)U, which is not typical of RNA nucleosides. Interestingly, the hydrogen bonding network in the crystals of dH(2)o(4)U stabilizes the sugar moiety conformation in the C3'-endo form (N conformer), rarely found in DNA nucleosides. In aqueous solution, dH(2)o(4)U reveals a similar population of the C2'-endo conformation (65%) to that of 2'-deoxy-2-thiouridine (62%), while the 62% population of the S conformer for H(2)o(4)U is significantly different from that of the parent 2-thiouridine, for which the N conformer is dominant (71%). Such a difference may be of biological importance, as the desulfurization process of natural tRNA 2-thiouridines may occur under conditions of oxidative stress in the cell and may influence the decoding process.

The intermolecular interactions in the aminonitromethylbenzenes

The intermolecular non-covalent interactions in aminonitromethylbenzenes namely 2-methyl-4-nitroaniline, 4-methyl-3-nitroaniline, 2-methyl-6-nitroaniline, 4-amino-2,6-dinitrotoluene, 2-methyl-5-nitroaniline, 4-methyl-2-nitroaniline, 2,3-dimethyl-6-nitroaniline, 4,5-dimethyl-2-nitroaniline and 2-methyl-3,5-dinitroaniline were studied by quantum mechanical calculations at RHF/311++G(3df,2p) and B3LYP/311++G(3df,2p) level of theory. The calculations prove that solely geometrical study of hydrogen bonding can be very misleading because not all short distances (classified as hydrogen bonds on the basis of interaction geometry) are bonding in character. For studied compounds interaction energy ranges from 0.23 kcal mol−1 to 5.59 kcal mol−1. The creation of intermolecular hydrogen bonds leads to charge redistribution in donors and acceptors. The Natural Bonding Orbitals analysis shows that hydrogen bonds are created by transfer of electron density from the lone pair orbitals of the H-bond acceptor to the antibonding molecular orbitals of the H-bond donor and Rydberg orbitals of the hydrogen atom. The stacking interactions are the interactions of delocalized molecular π-orbitals of the one molecule with delocalized antibonding molecular π-orbitals and the antibonding molecular σ-orbital created between the carbon atoms of the second aromatic ring and vice versa.