Substituted adenine quartets: interplay between substituent effect, hydrogen bonding, and aromaticity

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.

Impact of covalent modifications on the hydrogen bond strengths in diaminotriazine supramolecules

Melamine (M) is a popular triamine triazine compound in the field of supramolecular materials. In this work, we have computationally investigated how substituents can be exploited to improve the binding strength of M supramolecules. Two types of covalent modifications were studied: the substitution of an H atom within an amine group -NHR, and the replacement of the whole -NH 2 group (R = H, F, CH 3 and COCH 3 ). Through our dispersion-corrected density functional theory computations, we explain which covalent modification will show the best self-assembling capabilities, and why the binding energy is enhanced. Our charge density and molecular orbital analyses indicate that the best substituents are those that generate a charge accumulation on the endocyclic N atom, providing an improvement of the electrostatic attraction. At the same time the substituent assists the main N-H⋅⋅⋅N hydrogen bonds by interacting with the amino group of the other monomer. We also show how the selected group notably boosts the strength of hexameric rosettes. This research, therefore, provides molecular tools for the rational design of emerging materials based on uneven hydrogen-bonded arrangements.

Analysis of Conformational, Structural, Magnetic, and Electronic Properties Related to Antioxidant Activity: Revisiting Flavan, Anthocyanidin, Flavanone, Flavonol, Isoflavone, Flavone, and Flavan-3-ol

Understanding the antioxidant activity of flavonoids is important to investigate their biological activities as well as to design novel molecules with low toxicity and high activity. Aromaticity is a chemical property found in cyclic structures that plays an important role in their stability and reactivity, and its investigation can help us to understand the antioxidant activity of some heterocyclic compounds. In the present study, we applied the density functional theory (DFT) to investigate the properties of seven flavonoid structures with well-reported antioxidant activity: flavan, anthocyanidin, flavanone, flavonol, isoflavone, flavone, and flavan-3-ol. Conformational, structural, magnetic, and electronic analyses were performed using nuclear magnetic resonance, ionization potentials, electron affinity, bond dissociation energy, proton affinity, frontier molecular orbitals (highest occupied molecular orbital (HOMO)/lowest unoccupied molecular orbital (LUMO)), and aromaticity through nucleus-independent chemical shifts to analyze these seven flavonoid structures. We revised the influence of hydroxyl groups on the properties of flavonoids and also investigated the influence of the aromaticity of these seven flavonoids on the antioxidant activity.

Mutual Relations between Substituent Effect, Hydrogen Bonding, and Aromaticity in Adenine-Uracil and Adenine-Adenine Base Pairs

The electronic structure of substituted molecules is governed, to a significant extent, by the substituent effect (SE). In this paper, SEs in selected nucleic acid base pairs (Watson-Crick, Hoogsteen, adenine-adenine) are analyzed, with special emphasis on their influence on intramolecular interactions, aromaticity, and base pair hydrogen bonding. Quantum chemistry methods-DFT calculations, the natural bond orbital (NBO) approach, the Harmonic Oscillator Model of Aromaticity (HOMA) index, the charge of the substituent active region (cSAR) model, and the quantum theory of atoms in molecules (QTAIM)-are used to compare SEs acting on adenine moiety and H-bonds from various substitution positions. Comparisons of classical SEs in adenine with those observed in para- and meta-substituted benzenes allow for the better interpretation of the obtained results. Hydrogen bond stability and its other characteristics (e.g., covalency) can be significantly changed as a result of the SE, and its consequences are dependent on the substitution position. These changes allow us to investigate specific relations between H-bond parameters, leading to conclusions concerning the nature of hydrogen bonding in adenine dimers-e.g., H-bonds formed by five-membered ring nitrogen acceptor atoms have an inferior, less pronounced covalent nature as compared to those formed by six-membered ring nitrogen. The energies of individual H-bonds (obtained by the NBO method) are analyzed and compared to those predicted by the Espinosa-Molins-Lecomte (EML) model. Moreover, both SE and H-bonds can significantly affect the aromaticity of adenine rings; long-distance SEs on π-electron delocalization are also documented.