Investigations of Stacked DNA Base-Pair Steps: Highly Accurate Stacking Interaction Energies, Energy Decomposition, and Many-Body Stacking Effects

Intra-strand phosphate-mediated pathways in microsolvated double-stranded DNA

We argue that dry DNA charge transport in molecular junctions, over distances of tens of nanometers, can take place via independent intra-strand pathways involving the phosphate groups. Such pathways explain recent single-molecule experiments that compare currents in intact and nicked 100 base-pair double-stranded DNA. We explore the conditions that favor independent intra-strand transport channels with the participation of the phosphate groups, as opposed to purely base-mediated transport involving the pi-stacked bases and inter-strand transitions. Our computations demonstrate how long-distance transport pathways in DNA are tuned by the degree of solvation, which affects the level of dynamic disorder in the pi-stacking, and the energies of phosphate-group molecular orbitals.

Building an ab initio solvated DNA model using Euclidean neural networks

Accurately modeling large biomolecules such as DNA from first principles is fundamentally challenging due to the steep computational scaling of ab initio quantum chemistry methods. This limitation becomes even more prominent when modeling biomolecules in solution due to the need to include large numbers of solvent molecules. We present a machine-learned electron density model based on a Euclidean neural network framework that includes a built-in understanding of equivariance to model explicitly solvated double-stranded DNA. By training the machine learning model using molecular fragments that sample the key DNA and solvent interactions, we show that the model predicts electron densities of arbitrary systems of solvated DNA accurately, resolves polarization effects that are neglected by classical force fields, and captures the physics of the DNA-solvent interaction at the ab initio level.

A novel fluorescent sensor with an overtone peak reference for highly sensitive detection of mercury (II) ions and hydrogen sulfide: Mechanisms and applications in environmental monitoring and bioanalysis

The present study introduces a novel fluorescent sensor with an overtone peak reference designed for the detection of mercury (Ⅱ) ions (Hg2+) and hydrogen sulfide (H2S). The study proposes two novel response mechanisms that hinges on the synergistic effect of cation exchange dissociation (CED) and photo-induced electron transfer (PET). This sensor exhibits a remarkable detection limit of 2.9 nM for Hg2+. Additionally, the sensor reacts with H2S to generate nickel sulfide (NiS) semiconductor nanoparticles, which amplify the fluorescence signal and enable a detection limit of 3.1 nM for H2S. The detection limit for H2S is further improved to 29.1 pM through the surface functionalization of the nanomaterial with pyridine groups (increasing reactivity) and chelation of gold nanoparticles (AuNPs), which enhances the sensor's specificity. This improvement is primarily due to the surface plasmon resonance (SPR) of AuNPs and their affinity for H2S. The single-emission strategy can yield skewed results due to environmental changes, whereas the overtone peak reference strategy enhances result accuracy and reliability by detecting environmental interference through reference emission peaks. In another observation, the low-toxicity dihydropyrene-bipyridine nanorods (TPP-BPY) has been successfully utilized for both endogenous and exogenous H2S detection in vivo using a mouse model. The successful development of TPP-BPY is expected to provide an effective tool for studying the role of H2S in biomedical systems.

Binding Mechanism of Riboswitch to Natural Ligand Elucidated by McMD-Based Dynamic Docking Simulations

Flavin mononucleotide riboswitches are common among many pathogenic bacteria and are therefore considered to be an attractive target for antibiotics development. The riboswitch binds riboflavin (RBF, also known as vitamin B2), and although an experimental structure of their complex has been solved with the ligand bound deep inside the RNA molecule in a seemingly unreachable state, the binding mechanism between these molecules is not yet known. We have therefore used our Multicanonical Molecular Dynamics (McMD)-based dynamic docking protocol to analyze their binding mechanism by simulating the binding process between the riboswitch aptamer domain and the RBF, starting from the apo state of the riboswitch. Here, the refinement stage was crucial to identify the native binding configuration, as several other binding configurations were also found by McMD-based docking simulations. RBF initially binds the interface between P4 and P6 including U61 and G62, which forms a gateway where the ligand lingers until this gateway opens sufficiently to allow the ligand to pass through and slip into the hidden binding site including A48, A49, and A85.

Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA

Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.

Atomistic Picture of Opening-Closing Dynamics of DNA Holliday Junction Obtained by Molecular Simulations

Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.

Polymorphism of Butyl Ester of Oleanolic Acid—The Dominance of Dispersive Interactions over Electrostatic

Oleanolic (OA) and glycyrrhetinic acids (GE), as well as their derivatives, show a variety of pharmacological properties. Their crystal structures provide valuable information related to the assembly modes of these biologically active compounds. In the known-to-date crystals of OA esters, their 11-oxo derivatives, and GE ester crystals, triterpenes associate, forming different types of ribbons and layers whose construction is based mainly on van der Waals forces and weak C-H···O interactions. New crystal structures of 11-oxo OA methyl ester and the polymorph of OA butyl ester reveal an alternative aggregation mode. Supramolecular architectures consist of helical chains which are stabilized by hydrogen bonds of O-H···O type. It was found that two polymorphic forms of butyl OA ester (layered and helical) are related monotropically. In a structure of metastable form, O-H···O hydrogen bonds occur, while the thermodynamically preferred phase is governed mainly by van der Waals interactions. The intermolecular interaction energies calculated using CrystalExplorer, PIXEL, and Psi4 programs showed that even in motifs formed through O-H···O hydrogen bonds, the dispersive forces have a significant impact.

Strong External Electric Fields Reduce Explosive Sensitivity: A Theoretical Investigation into the Reaction Selectivity in NH2NO2∙∙∙NH3

Controlling the selectivity of a detonation initiation reaction of explosive is essential to reduce sensitivity, and it seems impossible to reduce it by strengthening the external electric field. To verify this, the effects of external electric fields on the initiation reactions in NH2NO2∙∙∙NH3, a model system of the nitroamine explosive with alkaline additive, were investigated at the MP2/6-311++G(2d,p) and CCSD(T)/6-311++G(2d,p) levels. The concerted effect in the intermolecular hydrogen exchange is characterized by an index of the imaginary vibrations. Due to the weakened concerted effects by the electric field along the −x-direction opposite to the “reaction axis”, the dominant reaction changes from the intermolecular hydrogen exchange to 1,3-intramolecular hydrogen transference with the increase in the field strengths. Furthermore, the stronger the field strengths, the higher the barrier heights become, indicating the lower sensitivities. Therefore, by increasing the field strength and adjusting the orientation between the field and “reaction axis”, not only can the reaction selectivity be controlled, but the sensitivity can also be reduced, in particular under a super-strong field. Thus, a traditional concept, in which the explosive is dangerous under the super-strong external electric field, is theoretically broken. Compared to the neutral medium, a low sensitivity of the explosive with alkaline can be achieved under the stronger field. Employing atoms in molecules, reduced density gradient, and surface electrostatic potentials, the origin of the reaction selectivity and sensitivity change is revealed. This work provides a new idea for the technical improvement regarding adding the external electric field into the explosive system.

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the "gold standard" coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller-Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

A comprehensive benchmark investigation of quantum chemical methods for carbocations

The application of various density functional approximations (DFAs) and an emphasis on popular methods without any consensus have prevailed in computational studies dedicated to carbocations. More importantly, an extensive and rigorous benchmark investigation on density functionals for the class is still lacking. To close this gap, we present a comprehensive benchmark study of quantum chemical methods on a series of classical and nonclassical carbocations, the CARBO33 dataset. We evaluate a total of 107 DFT methods from all rungs giving particular attention to double hybrid density functionals as the potential of the class has been largely undermined in the context of carbocations. To support our findings, DLPNO-CCSD(T) at the complete basis set (CBS) limit and W1-F12 are used as reference methods. Our results indicate that the composite CBS-QB3 method performs poorly and should not be adopted for target energies. Oftentimes, the tested DFAs of a lower rung perform better than several DFAs in a higher rung of Perdew's "Jacob's ladder". Nonetheless, double hybrids DSD-PBEP86-NL and ωB97X-2-D3(BJ) stand out by showing the overall best performance. Among the hybrids evaluated, about half of them show mean absolute deviation (MAD) below 1.1 kcal mol^-1, including the popular hybrids M06-2X and mPW1PW91. In this family, MN15-D3(BJ) performs particularly well (MAD = 0.77 kcal mol^-1) displaying reliable results across various tests. Highly popular B3LYP exhibited one of the worst performances (MAD = 4.74 kcal mol^-1), and we do not recommend its application to carbocations. We also assess the 24 general-purpose basis sets of single- up to quadruple-ζ quality. The best compromise between accuracy and computational cost is achieved with cc-pVTZ followed by def2-TZVP. Computations on larger structures of general interest, including terpene carbocations, are also presented for selected DFT methods confirming general trends in the results.

Polarizable Molecular Block Model: Toward the Development of an Induced Dipole Force Field for DNA

For flexible and highly ionized macromolecules such as DNA, it is important to correctly evaluate the intramolecular polarization in an induced dipole force field. In a proposed polarizable molecular block (PMB) model, a large molecule is divided into several molecular blocks. The atomic charges of the blocks are optimized by using the respective electrostatic potentials (ESPs) on the molecular surface. By using the capped hydrogen removal operation, the total charge of the blocks is controlled exactly to have an integer charge. The atomic polarizabilities of the blocks are optimized by using the respective polarized one-electron potentials that are the differences between ESPs with and without an external test charge. Induced dipole-charge interactions between the blocks are all included, but those interactions within the blocks are strictly excluded. All dipole-dipole interactions are included, but the damping functions are applied to the close dipole-dipole pairs. Several types of damping (simple scaling, exponential, linear, and Gaussian) are evaluated. The validity of the PMB model was verified by using trinucleotide duplexes which have A-, B-, and Z-DNA forms. The reference energies of trinucleotide duplexes including counterions (GGT3Na-ACC3Na, GAC3Na-GTC3Na, and GCG3Na-CGC3Na) are calculated using ωB97XD/aug-cc-pVDZ. All damping types reproduced well the reference interaction energies, dipole moments, and ESPs. Among them, the simple scaling with strong attenuation to 1-2 atomic pairs showed the highest stability against the polarization catastrophe. This study shows that it is possible to develop a high-quality polarizable force field by treating the intramolecular polarization on a block-by-block basis.

Computer Aided Development of Nucleic Acid Applications in Nanotechnologies

Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.

Predicting accurate ab initio DNA electron densities with equivariant neural networks

One of the fundamental limitations of accurately modeling biomolecules like DNA is the inability to perform quantum chemistry calculations on large molecular structures. We present a machine learning model based on an equivariant Euclidean neural network framework to obtain accurate ab initio electron densities for arbitrary DNA structures that are much too large for conventional quantum methods. The model is trained on representative B-DNA basepair steps that capture both base pairing and base stacking interactions. The model produces accurate electron densities for arbitrary B-DNA structures with typical errors of less than 1%. Crucially, the error does not increase with system size, which suggests that the model can extrapolate to large DNA structures with negligible loss of accuracy. The model also generalizes reasonably to other DNA structural motifs such as the A- and Z-DNA forms, despite being trained on only B-DNA configurations. The model is used to calculate electron densities of several large-scale DNA structures, and we show that the computational scaling for this model is essentially linear. We also show that this machine learning electron density model can be used to calculate accurate electrostatic potentials for DNA. These electrostatic potentials produce more accurate results compared with classical force fields and do not show the usual deficiencies at short range.

Slow Rotational Dynamics of Cytosine NH2 Groups in Double-Stranded DNA

Aromatic NH₂ groups are essential as hydrogen-bond donors in secondary structures of DNA and RNA. Although rapid rotations of NH₂ groups of adenine and guanine bases were previously characterized, there has been a lack of quantitative information about slow rotations of cytosine NH₂ groups in Watson-Crick base pairs. In this study, using an NMR method we had recently developed, we determined the kinetic rate constants and energy barriers for cytosine NH₂ rotations in a 15-base-pair DNA duplex. Our data show that the rotational dynamics of cytosine NH₂ groups depend on local environments. Qualitative correlation between the ranges of ¹⁵N chemical shifts and rotational time scales for various NH₂ groups of nucleic acids and proteins illuminates a relationship between the partial double-bond character of the C-N bond and the time scale for NH₂ rotations.

Nonadiabatic Molecular Dynamics with Extended Density Functional Tight-Binding: Application to Nanocrystals and Periodic Solids

In this work, we report a new methodology for nonadiabatic molecular dynamics calculations within the extended tight-binding (xTB) framework. We demonstrate the applicability of the developed approach to finite and periodic systems with thousands of atoms by modeling "hot" electron relaxation dynamics in silicon nanocrystals and electron-hole recombination in both a graphitic carbon nitride monolayer and a titanium-based metal-organic framework (MOF). This work reports the nonadiabatic dynamic simulations in the largest Si nanocrystals studied so far by the xTB framework, with diameters up to 3.5 nm. For silicon nanocrystals, we find a non-monotonic dependence of "hot" electron relaxation rates on the nanocrystal size, in agreement with available experimental reports. We rationalize this relationship by a combination of decreasing nonadiabatic couplings related to system size and the increase of available coherent transfer pathways in systems with higher densities of states. We emphasize the importance of proper treatment of coherences for obtaining such non-monotonic dependences. We characterize the electron-hole recombination dynamics in the graphitic carbon nitride monolayer and the Ti-containing MOF. We demonstrate the importance of spin-adaptation and proper sampling of surface hopping trajectories in modeling such processes. We also assess several trajectory surface hopping schemes and highlight their distinct qualitative behavior in modeling the excited-state dynamics in superexchange-like models depending on how they handle coherences between nearly parallel states.

Affinity and Correlation in DNA

A statistical analysis of important DNA sequences and related proteins has been performed to study the relationships between monomers, and some general considerations about these macromolecules can be provided from the results. First, the most important relationship between sites in all the DNA sequences examined is that between two consecutive base pairs. This is an indication of an energetic stabilization due to the stacking interaction of these couples of base pairs. Secondly, the difference between human chromosome sequences and their coding parts is relevant both in the relationships between sites and in some specific compositional rules, such as the second Chargaff rule. Third, the evidence of the relationship in two successive triplets of DNA coding sequences generates a relationship between two successive amino acids in the proteins. This is obviously impossible if all the relationships between the sites are statistical evidence and do not involve causes; therefore, in this article, due to stacking interactions and this relationship in coding sequences, we will divide the concept of the relationship between sites into two concepts: affinity and correlation, the first with physical causes and the second without. Finally, from the statistical analyses carried out, it will emerge that the human genome is uniform, with the only significant exception being the Y chromosome.

Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials

Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.

Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials

There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.

Orientation Dependence of DNA Blunt-End Stacking Studied by Free-Energy Simulations

DNA blunt ends can associate mediated by stacking interactions between the terminal base pairs that form blunt ends. The blunt end association plays a role in DNA repair and recombination processes and can also be of importance for the design of DNA-based nano-materials. Its function depends on the sequence and on the geometric arrangement that leads to stable interaction. For a stacked state, the relative orientation (twisting) of the base pairs is important. Molecular dynamics and advanced sampling simulations were used to calculate free energy change associated with twist changes of the stacked blunt-end base pairs. The calculations reproduce blunt stacking arrangements found in crystal structures of DNA oligonucleotides as free energy minima. To elucidate the physical origin of the stabilization of certain angular arrangements, the interactions between backbone atoms in the blunt-end stack were switched off in additional free energy calculations. It allows us to decipher the contributions to stacking stabilization due to the nucleobases and the backbone and to analyze the sequence dependence of the angular stacking preferences. Good qualitative agreement was also found for the comparison with quantum mechanical calculations. The results may help in the design of novel DNA-based materials.

Short-Range Imbalances in the AMBER Lennard-Jones Potential for (Deoxy)Ribose···Nucleobase Lone-Pair···π Contacts in Nucleic Acids

The lone-pair···π (lp···π) (deoxy)ribose···nucleobase stacking is a recurring interaction in Z-DNA and RNAs that is characterized by sub-van der Waals lp···π contacts (<3.0 Å). It is a part of the structural signature of CpG Z-step motifs in Z-DNA and r(UNCG) tetraloops that are known to behave poorly in molecular dynamics (MD) simulations. Although the exact origin of the MD simulation issues remains unclear, a significant part of the problem might be due to an imbalanced description of nonbonded interactions, including the characteristic lp···π stacking. To gain insights into the links between lp···π stacking and MD, we present an in-depth comparison between accurate large-basis-set double-hybrid Kohn-Sham density functional theory calculations DSD-BLYP-D3/ma-def2-QZVPP (DHDF-D3) and data obtained with the nonbonded potential of the AMBER force field (AFF) for NpN Z-steps (N = G, A, C, and U). Among other differences, we found that the AFF overestimates the DHDF-D3 lp···π distances by ∼0.1-0.2 Å, while the deviation between the DHDF-D3 and AFF descriptions sharply increases in the short-range region of the interaction. Based on atom-in-molecule polarizabilities and symmetry-adapted perturbation theory analysis, we inferred that the DHDF-D3 versus AFF differences partly originate in identical nucleobase carbon atom Lennard-Jones (LJ) parameters despite the presence/absence of connected electron-withdrawing groups that lead to different effective volumes or vdW radii. Thus, to precisely model the very short CpG lp···π contact distances, we recommend revision of the nucleobase atom LJ parameters. Additionally, we suggest that the large discrepancy between DHDF-D3 and AFF short-range repulsive part of the interaction energy potential may significantly contribute to the poor performances of MD simulations of nucleic acid systems containing Z-steps. Understanding where, and if possible why, the point-charge-type effective potentials reach their limits is vital for developing next-generation FFs and for addressing specific issues in contemporary MD simulations.

Unveiling the complex pattern of intermolecular interactions responsible for the stability of the DNA duplex

Herein, we provide new insights into the intermolecular interactions responsible for the intrinsic stability of the duplex structure of a large portion of human B-DNA by using advanced quantum mechanical methods. Our results indicate that (i) the effect of non-neighboring bases on the inter-strand interaction is negligibly small, (ii) London dispersion effects are essential for the stability of the duplex structure, (iii) the largest contribution to the stability of the duplex structure is the Watson-Crick base pairing - consistent with previous computational investigations, (iv) the effect of stacking between adjacent bases is relatively small but still essential for the duplex structure stability and (v) there are no cooperativity effects between intra-strand stacking and inter-strand base pairing interactions. These results are consistent with atomic force microscope measurements and provide the first theoretical validation of nearest neighbor approaches for predicting thermodynamic data of arbitrary DNA sequences.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

Theoretical investigation of conformational deviation of the human parallel telomeric G-quadruplex DNA in the presence of different salt concentrations and temperatures under confinement

Various experimental reports address the stability of G-quadruplex DNA inside a close confinement such as α-hemolysin, nanocavity water pool and different metal-organic-frameworks (MOFs). To understand the conformational change of G-quadruplex DNA at the atomistic level, we have carried out a total of 40 μs simulation run under both non-polar and polar confinement conditions. To investigate the dynamics, we have considered two different KCl salt concentrations, i.e., 0.47 M (minimal salt concentration) and higher than 2 M (higher salt concentration), at two distinct temperatures, 300 K and 350 K. Here, we have observed that the human telomeric G-quadruplex DNA deviates more from its crystal structure at minimal salt concentration under both non-polar and polar confinement conditions. Besides, the loop regions deviate and fluctuate more compared to the other regions, i.e., sugar-phosphate backbone and tetrad regions. The presence of K+ ions is found to be primarily responsible for this phenomenon. From the spatial density function (SDF) plots, a higher density of K+ ions is observed in the backbone region. Furthermore, from the residue-wise first solvation shell estimation, we have noticed that the K+ ions mainly accumulate in the tetrad region under both non-polar and polar confinement conditions due to which the tetrad regions are more rigid than the loop regions. Higher salt concentration results in increased rigidity of the G-quadruplex DNA. Our study provides valuable insight into the conformational deviation of the G-quadruplex DNA under nanoconfinement conditions.

Analysis of DNA interactions and GC content with energy decomposition in large-scale quantum mechanical calculations

GC content is a contributing factor to the stability of nucleic acids due to hydrogen bonding. More hydrogen bonding generally results in greater stability. Empirical evidence, however, has suggested that the GC content of a nucleic acid is a poor predictor of its stability, implying that there are sequence-dependent interactions besides what its GC content indicates. To examine how much such sequence-dependent interactions affect the interaction energies of double-stranded DNA (dsDNA) molecules, dsDNA molecules of different sequences are generated and examined in silico for variabilities in the interaction energies within each group of dsDNA molecules of the same GC content. Since the amount of hydrogen bonding depends on the GC content, holding the GC content fixed when examining the differences in interaction energies allows sequence-dependent interactions to be isolated. The nature of sequence-dependent interactions is then dissected using energy decomposition analysis (EDA). By using EDA, the components of the interactions that depend on the neighboring base pairs help explain some of the variability in the interaction energies of the dsDNA molecules despite having the same GC content. This work provides a new paradigm and tool for the study and analysis of the distributions of interaction components in dsDNA with the same GC content using EDA within large-scale quantum chemistry calculations.

Replacing thymine with a strongly pairing fifth Base: A combined quantum mechanics and molecular dynamics study

The non-natural ethynylmethylpyridone C-nucleoside (W), a thymidine (T) analogue that can be incorporated in oligonucleotides by automated synthesis, has recently been reported to form a high fidelity base pair with adenosine (A) and to be well accommodated in B-DNA duplexes. The enhanced binding affinity for A of W, as compared to T, makes it an ideal modification for biotechnological applications, such as efficient probe hybridization for the parallel detection of multiple DNA strands. In order to complement the experimental study and rationalize the impact of the non-natural W nucleoside on the structure, stability and dynamics of DNA structures, we performed quantum mechanics (QM) calculations along with molecular dynamics (MD) simulations. Consistently with the experimental study, our QM calculations show that the A:W base pair has an increased stability as compared to the natural A:T pair, due to an additional CH-π interaction. Furthermore, we show that mispairing between W and guanine (G) causes a distortion in the planarity of the base pair, thus explaining the destabilization of DNA duplexes featuring a G:W pair. MD simulations show that incorporation of single or multiple consecutive A:W pairs in DNA duplexes causes minor changes to the intra- and inter-base geometrical parameters, while a moderate widening/shrinking of the major/minor groove of the duplexes is observed. QM calculations applied to selected stacks from the MD simulations also show an increased stacking energy for W, over T, with the neighboring bases.

An in silico investigation of the binding modes and pathway of APTO-253 on c-KIT G-quadruplex DNA

The stability of c-KIT G-quadruplex DNA via ligands has been a significant concern in the growing field of cancer therapy. Thus, it is very important to understand the mechanism behind the high binding affinity of the small drug molecules on the c-KIT G-quadruplex DNA. In this study, we have investigated the binding mode and pathway of the APTO-253 ligand on the c-KIT G-quadruplex DNA employing a total of 10 μs all atom molecular dynamics simulations and further 8.82 μs simulations via the umbrella sampling method using both OL15 and BSC1 latest force fields for DNA structures. From the cluster structure analysis, mainly three binding pathways i.e., top, bottom and side loop stacking modes are identified. Moreover, RMSD, RMSF and 2D-RMSD values indicate that the c-KIT G-quadruplex DNA and APTO-253 molecules are stable throughout the simulation run. Furthermore, the number of hydrogen bonds in each tetrad and the distance between the two central K⁺ cations confirm that the c-KIT G-quadruplex DNA maintains its conformation in the process of complex formation with the APTO-253 ligand. The binding free energies and the minimum values in the potential of mean forces suggest that the binding processes are energetically favorable. Furthermore, we have found that the bottom stacking mode is the most favorable binding mode among all the three modes for the OL15 force field. However, for the BSC1 force field, both the top and bottom binding modes of the APTO-253 ligand in c-KIT G-quadruplex DNA are comparable to each other. To investigate the driving force for the complex formation, we have noticed that the van der Waals (vdW) and π-π stacking interactions are mainly responsible. Our detailed studies provide useful information for the discovery of novel drugs in the field of stabilization of G-quadruplex DNAs.

The influence of base pair tautomerism on single point mutations in aqueous DNA

The relationship between base pair hydrogen bond proton transfer and the rate of spontaneous single point mutations at ambient temperatures and pressures in aqueous DNA is investigated. By using an ensemble-based multiscale computational modelling method, statistically robust rates of proton transfer for the A:T and G:C base pairs within a solvated DNA dodecamer are calculated. Several different proton transfer pathways are observed within the same base pair. It is shown that, in G:C, the double proton transfer tautomer is preferred, while the single proton transfer process is favoured in A:T. The reported range of rate coefficients for double proton transfer is consistent with recent experimental data. Notwithstanding the approximately 1000 times more common presence of single proton transfer products from A:T, observationally there is bias towards G:C to A:T mutations in a wide range of living organisms. We infer that the double proton transfer reactions between G:C base pairs have a negligible contribution towards this bias for the following reasons: (i) the maximum half-life of the G*:C* tautomer is in the range of picoseconds, which is significantly smaller than the milliseconds it takes for DNA to unwind during replication, (ii) statistically, the majority of G*:C* tautomers revert back to their canonical forms through a barrierless process, and (iii) the thermodynamic instability of the tautomers with respect to the canonical base pairs. Through similar reasoning, we also deduce that proton transfer in the A:T base pair does not contribute to single point mutations in DNA.

Importance of base-pair opening for mismatch recognition

Mismatch repair is a highly conserved cellular pathway responsible for repairing mismatched dsDNA. Errors are detected by the MutS enzyme, which most likely senses altered mechanical property of damaged dsDNA rather than a specific molecular pattern. While the curved shape of dsDNA in crystallographic MutS/DNA structures suggests the role of DNA bending, the theoretical support is not fully convincing. Here, we present a computational study focused on a base-pair opening into the minor groove, a specific base-pair motion observed upon interaction with MutS. Propensities for the opening were evaluated in terms of two base-pair parameters: Opening and Shear. We tested all possible base pairs in anti/anti, anti/syn and syn/anti orientations and found clear discrimination between mismatches and canonical base-pairs only for the opening into the minor groove. Besides, the discrimination gap was also confirmed in hotspot and coldspot sequences, indicating that the opening could play a more significant role in the mismatch recognition than previously recognized. Our findings can be helpful for a better understanding of sequence-dependent mutability. Further, detailed structural characterization of mismatches can serve for designing anti-cancer drugs targeting mismatched base pairs.

Biological activity of quinazoline analogues and molecular modeling of their interactions with G-quadruplexes

BACKGROUND

Quinazolines 1 to 6, with an aromatic or aryl-vinyl substituent in position 2 are selected with the aim to compare their structures and biological activity. The selection includes a natural alkaloid, schizocommunin, and the synthetic 2-(2'-quinolyl)-3H-quinazolin-4-one, known to interact with guanine-quadruplex dependent enzymes, respectively telomerase and topoisomerase.

METHODS

Breast cancer cells of the MDA cell line have been used to study the bioactivity of the tested compounds by the method of Comet Assay and FACS analyses. We model observed effects assuming stacking interactions of studied heterocycles with a naked skeleton of G-quadruplex, consisting of guanine quartet layers and potassium ions. Interaction energies are computed using a dispersion corrected density functional theory method, and an electron-correlated molecular orbital theory method.

RESULTS

Selected compounds do not remarkably delay nor change the dynamics of cellular progression through the cell cycle phases, while changing significantly cell morphology. Our computational models quantify structural effects on heterocyclic G4-complex stabilization energies, which directly correlate with observed biological activity.

CONCLUSION

Our computational model of G-quadruplexes is an acceptable tool for the study of interaction energies of G-quadruplexes and heterocyclic ligands, predicting, and allowing design of novel structures.

GENERAL SIGNIFICANCE

Genotoxicity of quinazolin-4-one analogues on human breast cancer cells is not related to molecular metabolism but rather to their interference with G-quadruplex regulatory mechanisms. Computed stabilization energies of heterocyclic ligand complexes of G-quadruplexes might be useful in the prediction of novel telomerase / helicase, topoisomerase and NA polymerase dependent drugs.

Do weak interactions affect the biological behavior of DNA? A DFT study of CpG island-like chains

The origin, stability, and contribution to the formation of noncovalent interactions, such as hydrogen bonds and π - π stacking, have been already widely discussed. However, there are few discussions about the relevance of these weak interactions in DNA performance. In this work, we seek to shed light on the effect of hydrogen bonds and π - π stacking interactions on the biological behavior of DNA through the description of these intermolecular forces in CpG island-like (GC-rich) chains. Furthermore, we made some comparisons with TATA box-like (TA-rich) chains in order to describe hydrogen bond and π - π stacking interactions as a function of the DNA sequence. For hydrogen bonds, we found that there is not a significant effect related to the number of base pairs. Whereas for π - π stacking interactions, the energy tended to decrease as the number of base pairs increased. We observed anticooperative effects for both hydrogen bonds and π - π stacking interactions. These results are in contrast with those of TATA box-like chains since cooperative and additive effects were found for both hydrogen bonds and π - π stacking, respectively. Based on the chemical hardness and density of states, we can conclude that proteins may interact easier with GC-rich chains. We conclude that regardless of the chain length, a protein could interact more easily with these genomics regions because the π - π stacking energies did not increase as a function of the number of base pairs, making, for the first time, a first approximation of the influence of noncovalent interaction on DNA behavior. We did all this work by means of DFT framework included in the DMol³ code (M06-L/DNP). Graphical Abstract Cartoon representation of how nocovalent interactions affect the interaction of DNA with a protein, i.e., how hydrogen bond and π - π stacking interactions influence the biological behavior of DNA.

Surprisingly broad applicability of the cc-pVnZ-F12 basis set for ground and excited states

Excellent convergence properties for the (aug-)cc-pVnZ-F12 basis set family, purpose-made for explicitly correlated calculations, are demonstrated with conventional wave function methods and Kohn-Sham density functional theory for various ground and excited-state calculations. Among the ground-state properties studied are dipole moments, covalent bond lengths, and interaction and reaction energies. For excited states, we looked at vertical excitation energies, UV absorption, and excited-state absorption spectra. Convergence is compared against the basis sets cc-pVnZ, def2-nVD, aug-pcseg-n, and nZaPa-NR. It is established that the cc-pVnZ-F12 family consistently yields results of n + 1 quality and better. Especially, the cc-pVDZ-F12 basis set is found to be a basis set of good cost vs performance trade-off.

Fine-Tuning of the AMBER RNA Force Field with a New Term Adjusting Interactions of Terminal Nucleotides

Determination of RNA structural-dynamic properties is challenging for experimental methods. Thus, atomistic molecular dynamics (MD) simulations represent a helpful technique complementary to experiments. However, contemporary MD methods still suffer from limitations of force fields (ffs), including imbalances in the nonbonded ff terms. We have recently demonstrated that some improvement of state-of-the-art AMBER RNA ff can be achieved by adding a new term for H-bonding called gHBfix, which increases tuning flexibility and reduces risk of side-effects. Still, the first gHBfix version did not fully correct simulations of short RNA tetranucleotides (TNs). TNs are key benchmark systems due to availability of unique NMR data, although giving too much weight on improving TN simulations can easily lead to overfitting to A-form RNA. Here we combine the gHBfix version with another term called tHBfix, which separately treats H-bond interactions formed by terminal nucleotides. This allows to refine simulations of RNA TNs without affecting simulations of other RNAs. The approach is in line with adopted strategy of current RNA ffs, where the terminal nucleotides possess different parameters for terminal atoms than the internal nucleotides. Combination of gHBfix with tHBfix significantly improves the behavior of RNA TNs during well-converged enhanced-sampling simulations using replica exchange with solute tempering. TNs mostly populate canonical A-form like states while spurious intercalated structures are largely suppressed. Still, simulations of r(AAAA) and r(UUUU) TNs show some residual discrepancies with primary NMR data which suggests that future tuning of some other ff terms might be useful. Nevertheless, the tHBfix has a clear potential to improve modeling of key biochemical processes, where interactions of RNA single stranded ends are involved.

Understanding The Role of Reline, a Natural DES, on Temperature-Induced Conformational Changes of C-Kit G-Quadruplex DNA: A Molecular Dynamics Study

The noncanonical guanine-rich DNAs have drawn particular attention to the scientific world due to their controllable diverse and polymorphic structures. Apart from biological and medical significance, G-quadruplex DNAs are widely used in various fields such as nanotechnology, nanomachine, biosensors, and biocatalyst. So far, the applications of the G-quadruplex DNA are mainly limited in the water medium. Recently, a new generation of solvent named deep eutectic solvent (DES) has become very popular and has been widely used as a reaction medium of biocatalytic reactions and long-term storage medium for nucleic acids, even at high temperature. Hence, it is essential to understand the role of DES on temperature-induced conformational changes of a G-quadruplex DNA. In this research work, we have explored the temperature-mediated conformational dynamics of c-kit oncogene promoter G-quadruplex DNA in reline medium in the temperature range of 300-500 K, using a total of 10 μs unbiased all-atom molecular dynamics simulation. Here, from RMSD, RMSF, R_g and principal component analyses, we notice that the c-kit G-quadruplex DNA is stable up to 450 K in reline medium. However, it unfolds in water medium at 450 K. It is found that the hydrogen bonding interactions between c-kit G-quadruplex DNA and reline play a key role in the stabilization of the G-quadruplex DNA even at high temperature. Furthermore, in this work we have observed a very interesting and distinctive phenomenon of the central cation of the G-quadruplex DNA. Its position was seen to fluctuate between the two tetrad cores, that is, the region between tetrad-1 and tetrad-2 and that between tetrad-2 and tetrad-3 and vice versa at 450 and 500 K in reline medium which is absent in water medium at 450 K. Moreover, the rate of its oscillation is increased when temperature is increased.

Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction

BACKGROUND

Molecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs.

METHODS

We employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair.

RESULTS

We noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear.

CONCLUSIONS

The quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs "tug of war" situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity.

GENERAL SIGNIFICANCE

We have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.

Toward a less costly but accurate calculation of the CCSD(T)/CBS noncovalent interaction energy

The popular method of calculating the noncovalent interaction energies at the coupled-cluster single-, double-, and perturbative triple-excitations [CCSD(T)] theory level in the complete basis set (CBS) limit was to add a CCSD(T) correction term to the CBS second-order Møller-Plesset perturbation theory (MP2). The CCSD(T) correction term is the difference between the CCSD(T) and MP2 interaction energies evaluated in a medium basis set. However, the CCSD(T) calculations with the medium basis sets are still very expensive for systems with more than 30 atoms. Comparatively, the domain-based local pair natural orbital coupled-cluster method [DLPNO-CCSD(T)] can be applied to large systems with over 1,000 atoms. Considering both the computational accuracy and efficiency, in this work, we propose a new scheme to calculate the CCSD(T)/CBS interaction energies. In this scheme, the MP2/CBS term keeps intact and the CCSD(T) correction term is replaced by a DLPNO-CCSD(T) correction term which is the difference between the DLPNO-CCSD(T) and DLPNO-MP2 interaction energies evaluated in a medium basis set. The interaction energies of the noncovalent systems in the S22, HSG, HBC6, NBC10, and S66 databases were recalculated employing this new scheme. The consistent and tight settings of the truncation parameters for DLPNO-CCSD(T) and DLPNO-MP2 in this noncanonical CCSD(T)/CBS calculations lead to the maximum absolute deviation and root-mean-square deviation from the canonical CCSD(T)/CBS interaction energies of less than or equal to 0.28 kcal/mol and 0.09 kcal/mol, respectively. The high accuracy and low cost of this new computational scheme make it an excellent candidate for the study of large noncovalent systems.

Importance of model size in quantum mechanical studies of DNA intercalation

The convergence of DFT-computed interaction energies with increasing binding site model size was assessed. The data show that while accurate intercalator interaction energies can be derived from binding site models featuring only the flanking nucleotides for uncharged intercalators that bind parallel to the DNA base pairs, errors remain significant even when including distant nucleotides for intercalators that are charged, exhibit groove-binding tails that engage in noncovalent interactions with distant nucleotides, or that bind perpendicular to the DNA base pairs. Consequently, binding site models that include at least three adjacent nucleotides are required to consistently predict converged binding energies. The computationally inexpensive HF-3c method is shown to provide reliable interaction energies and can be routinely applied to such large models.

Establishing the accuracy of density functional approaches for the description of noncovalent interactions in biomolecules

Biomolecules have complex structures, and noncovalent interactions are crucial to determine their conformations and functionalities.

Revisiting the Potential Energy Surface of the Stacked Cytosine Dimer: FNO-CCSD(T) Interaction Energies, SAPT Decompositions, and Benchmarking

Nucleobase stacking interactions are crucial for the stability of nucleic acids. This study investigates base stacking energies of the cytosine homodimer in different configurations, including intermolecular separation plots, detailed twist dependence, and displaced structures. Highly accurate ab initio quantum chemical single point energies using an energy function based on MP2 complete basis set extrapolation ([6 → 7]ZaPa-NR) and a CCSD(T)/cc-pVTZ-F12 high-level correction are presented as new reference data, providing the most accurate stacking energies of nucleobase dimers currently available. Accurate SAPT2+(3)δMP2 energy decomposition is used to obtain detailed insights into the nature of base stacking interactions at varying vertical distances and twist values. The ab initio symmetry adapted perturbation theory (SAPT) energy decomposition suggests that the base stacking originates from an intricate interplay between dispersion attraction, short-range exchange-repulsion, and Coulomb interaction. The interpretation of the SAPT data is a complex issue as key energy terms vary substantially in the region of optimal (low energy) base stacking geometries. Thus, attempts to highlight one leading stabilizing SAPT base stacking term may be misleading and the outcome strongly depends on the used geometries within the range of geometries sampled in nucleic acids upon thermal fluctuations. Modern dispersion-corrected density functional theory (among them DSD-BLYP-D3, ωB97M-V, and ωB97M-D3BJ) is benchmarked and often reaches up to spectroscopic accuracy (below 1 kJ/mol). The classical AMBER force field is benchmarked with multiple different sets of point-charges (e.g. HF, DFT, and MP2-based) and is found to produce reasonable agreement with the benchmark data.

Exploring Nature and Predicting Strength of Hydrogen Bonds: A Correlation Analysis Between Atoms-in-Molecules Descriptors, Binding Energies, and Energy Components of Symmetry-Adapted Perturbation Theory

This work studies the underlying nature of H-bonds (HBs) of different types and strengths and tries to predict binding energies (BEs) based on the properties derived from wave function analysis. A total of 42 HB complexes constructed from 28 neutral and 14 charged monomers were considered. This set was designed to sample a wide range of HB strengths to obtain a complete view about HBs. BEs were derived with the accurate coupled cluster singles and doubles with perturbative triples correction (CCSD(T))(T) method and the physical components of the BE were investigated by symmetry-adapted perturbation theory (SAPT). Quantum theory of atoms-in-molecules (QTAIM) descriptors and other HB indices were calculated based on high-quality density functional theory wave functions. We propose a new and rigorous classification of H-bonds (HBs) based on the SAPT decomposition. Neutral complexes are either classified as "very weak" HBs with a BE ≥ -2.5 kcal/mol that are mainly dominated by both dispersion and electrostatic interactions or as "weak-to-medium" HBs with a BE varying between -2.5 and -14.0 kcal/mol that are only dominated by electrostatic interactions. On the other hand, charged complexes are divided into "medium" HBs with a BE in the range of -11.0 to -15.0 kcal/mol, which are mainly dominated by electrostatic interactions, or into "strong" HBs whose BE is more negative than -15.0 kcal/mol, which are mainly dominated by electrostatic together with induction interactions. Among various explored correlations between BEs and wave function-based HB descriptors, a fairly satisfactory correlation was found for the electron density at the bond critical point (BCP; ρ_BCP ) of HBs. The fitted equation for neutral complexes is BE/kcal/mol =  - 223.08 × ρ_BCP /a. u. + 0.7423 with a mean absolute percentage error (MAPE) of 14.7%, while that for charged complexes is BE/kcal/mol =  - 332.34 × ρ_BCP /a. u. - 1.0661 with a MAPE of 10.0%. In practice, these equations may be used for a quick estimation of HB BEs, for example, for intramolecular HBs or large HB networks in biomolecules. © 2019 Wiley Periodicals, Inc.

Structural and Energetic Impact of Non-natural 7-Deaza-8-azaguanine, 7-Deaza-8-azaisoguanine, and Their 7-Substituted Derivatives on Hydrogen-Bond Pairing with Cytosine and Isocytosine

The impact of 7-deaza-8-azaguanine (DAG) and 7-deaza-8-azaisoguanine (DAiG) modifications on the geometry and stability of the G:C Watson-Crick (cWW) base pair and the G:iC and iG:C reverse Watson-Crick (tWW) base pairs has been characterized theoretically. In addition, the effect on the same base pairs of seven C7-substituted DAG and DAiG derivatives, some of which have been previously experimentally characterized, has been investigated. Calculations indicate that all of these modifications have a negligible impact on the geometry of the above base pairs, and that modification of the heterocycle skeleton has a small impact on the base-pair interaction energies. Instead, base-pair interaction energies are dependent on the nature of the C7 substituent. For the 7-substituted DAG-C cWW systems, a linear correlation between the base-pair interaction energy and the Hammett constant of the 7-substituent is found, with higher interaction energies corresponding to more electron-withdrawing substituents. Therefore, the explored modifications are expected to be accommodated in both parallel and antiparallel nucleic acid duplexes without perturbing their geometry, while the strength of a base pair (and duplex) featuring a DAG modification can, in principle, be tuned by incorporating different substituents at the C7 position.