Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations

Comparing ANI-2x, ANI-1ccx neural networks, force field, and DFT methods for predicting conformational potential energy of organic molecules

In this study, the conformational potential energy surfaces of Amylmetacresol, Benzocaine, Dopamine, Betazole, and Betahistine molecules were scanned and analyzed using the neural network architecture ANI-2 × and ANI-1ccx, the force field method OPLS, and density functional theory with the exchange-correlation functional B3LYP and the basis set 6-31G(d). The ANI-1ccx and ANI-2 × methods demonstrated the highest accuracy in predicting torsional energy profiles, effectively capturing the minimum and maximum values of these profiles. Conformational potential energy values calculated by B3LYP and the OPLS force field method differ from those calculated by ANI-1ccx and ANI-2x, which account for non-bonded intramolecular interactions, since the B3LYP functional and OPLS force field weakly consider van der Waals and other intramolecular forces in torsional energy profiles. For a more comprehensive analysis, electronic parameters such as dipole moment, HOMO, and LUMO energies for different torsional angles were calculated at two levels of theory, B3LYP/6-31G(d) and ωB97X/6-31G(d). These calculations confirmed that ANI predictions are more accurate than density functional theory calculations with B3LYP functional and OPLS force field for determining potential energy surfaces. This research successfully addressed the challenges in determining conformational potential energy levels and shows how machine learning and deep neural networks offer a more accurate, cost-effective, and rapid alternative for predicting torsional energy profiles.

A density functional theory study on interactions in water-bridged dimeric complexes of lignin

Lignin is the main plant cell wall component responsible for recalcitrance in the process of lignocellulosic biomass conversion into biofuels. The recalcitrance and insolubility of lignin in different reaction media are due in part to the hydrogen bonds and π interactions that hold syringyl (S) and guaiacyl (G) units together and promote the formation of stable water-bridged dimeric complexes (WBDCs): S⋯G and S⋯S, in native lignin. The current understanding of how each type of interaction influences the stability of these complexes within lignin native cell walls is still limited. Here, we found by DFT calculations that hydrogen bonding is more dominant than π-stacking interaction between aromatic rings of WBDCs. Although there is a stronger interaction of hydrogen bonds between subunits and water and higher π-stacking interaction in the S⋯S complex compared to the S⋯G complex, the former complex is less thermodynamically stable than the latter due to the entropic contribution coming from the methoxy substituents in the S-unit. Our results demonstrate that the methoxylation degree of lignin units does not significantly influence the structural geometries of WBDCs; if anything, an enhanced dispersion interaction between ring aromatics results in quasi-sandwich geometries as found in "coiled" lignin structures in the xylem tissue of wood. In the same way as that with ionic liquids, polar solvents can dissolve S-lignin by favorable interactions with the aliphatic hydroxyl group in the α-position as the key site or the aromatic hydroxyl group as the secondary site.

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.

The Role of Bond Functions in Describing Intermolecular Electron Correlation for Van der Waals Dimers: A Study of (CH4)2 and Ne2

We present a study of the intermolecular interactions in van der Waals complexes of methane and neon dimers within the framework of the CCSD method. This approach was implemented and applied to calculate and examine the behavior of the contracted two-particle reduced density matrix (2-RDM). It was demonstrated that the region near the minimum of the two-particle density matrix correlation part, corresponding to the primary bulk of the Coulomb hole contribution, exerts a significant influence on the dispersion interaction energetics of the studied systems. As a result, the bond functions approach was applied to improve the convergence performance for the intermolecular correlation energy results with respect to the size of the atomic basis. For this, substantial acceleration was achieved by introducing an auxiliary basis of bond functions centered on the minima of the 2-RDM. For both methane and neon dimers, this general conclusion was confirmed with a series of CCSD calculations for the 2-RDM and the correlation energies.

NICE-FF: A non-empirical, intermolecular, consistent, and extensible force field for nucleic acids and beyond

A new non-empirical ab initio intermolecular force field (NICE-FF in buffered 14-7 potential form) has been developed for nucleic acids and beyond based on the dimer interaction energies (IEs) calculated at the spin component scaled-MI-second order Møller-Plesset perturbation theory. A fully automatic framework has been implemented for this purpose, capable of generating well-polished computational grids, performing the necessary ab initio calculations, conducting machine learning (ML) assisted force field (FF) parametrization, and extending existing FF parameters by incorporating new atom types. For the ML-assisted parametrization of NICE-FF, interaction energies of ∼18 000 dimer geometries (with IE < 0) were used, and the best fit gave a mean square deviation of about 0.46 kcal/mol. During this parametrization, atom types apparent in four deoxyribonucleic acid (DNA) bases have been first trained using the generated DNA base datasets. Both uracil and hypoxanthine, which contain the same atom types found in DNA bases, have been considered as test molecules. Three new atom types have been added to the DNA atom types by using IE datasets of both pyrazinamide and 9-methylhypoxanthine. Finally, the last test molecule, theophylline, has been selected, which contains already-fitted atom-type parameters. The performance of NICE-FF has been investigated on the S22 dataset, and it has been found that NICE-FF outperforms the well-known FFs by generating the most consistent IEs with the high-level ab initio ones. Moreover, NICE-FF has been integrated into our in-house developed crystal structure prediction (CSP) tool [called FFCASP (Fast and Flexible CrystAl Structure Predictor)], aiming to find the experimental crystal structures of all considered molecules. CSPs, which were performed up to 4 formula units (Z), resulted in NICE-FF being able to locate almost all the known experimental crystal structures with sufficiently low RMSD20 values to provide good starting points for density functional theory optimizations.

Molecular van der Waals Fluids in Cavity Quantum Electrodynamics

Intermolecular van der Waals interactions are central to chemical and physical phenomena ranging from biomolecule binding to soft-matter phase transitions. In this work, we demonstrate that strong light-matter coupling can be used to control the thermodynamic properties of many-molecule systems. Our analyses reveal orientation dependent single molecule energies and interaction energies for van der Waals molecules. For example, we find intermolecular interactions that depend on the distance between the molecules R as R-3 and R0. Moreover, we employ ab initio cavity quantum electrodynamics calculations to develop machine-learning-based interaction potentials for molecules inside optical cavities. By simulating systems ranging from 12 H2 to 144 H2 molecules, we observe varying degrees of orientational order because of cavity-modified interactions, and we explain how quantum nuclear effects, light-matter coupling strengths, number of cavity modes, molecular anisotropies, and system size all impact the extent of orientational order.

MP2-Based Correction Scheme to Approach the Limit of a Complete Pair Natural Orbitals Space in DLPNO-CCSD(T) Calculations

The domain-based local pair natural orbital (PNO) coupled-cluster DLPNO-CCSD(T) method has been proven to provide accurate single-point energies at a fraction of the cost of canonical CCSD(T) calculations. However, the desired "chemical accuracy" can only be obtained with a large PNO space and extended basis set. We present a simple yet accurate and efficient correction scheme based on a perturbative approach. Here, in addition to DLPNO-CCSD(T) energy, one calculates DLPNO-MP2 correlation energy with the same settings as in the preceding coupled-cluster calculation. In the next step, the canonical MP2 correlation energy is obtained in the same orbital basis. This can be efficiently performed for essentially all molecule sizes accessible with the DLPNO-CCSD(T) method. By taking the difference between the canonical MP2 and DLPNO-MP2 energies, we obtain a correction term that can be added to the DLPNO-CCSD(T) correlation energy. This way, one can obtain the total correlation energy close to the limit of the complete PNO space (cPNO). The presented approach allows us to significantly increase the accuracy of the DLPNO-CCSD(T) method for both closed- and open-shell systems. The latter are known to be especially challenging for locally correlated methods. Unlike the previously developed PNO extrapolation procedure by Altun, Neese, and Bistoni ( J. Chem. Theory Comput. 2020, 16, 6142-6149), this strategy allows us to get the DLPNO-CCSD(T) correlation energy at the cPNO limit in a cost-efficient way, resulting in a minimal overall increase in calculation time as compared to the uncorrected method.

High-Order Quantum-Mechanical Analysis of Hydrogen Bonding in Hachimoji and Natural DNA Base Pairs

High-order quantum chemistry is applied to hydrogen-bonded natural DNA nucleobase pairs [adenine:thymine (A:T) and guanine:cytosine (G:C)] and non-natural Hachimoji nucleobase pairs [isoguanine:1-methylcytosine (B:S) and 2-aminoimidazo[1,2a][1,3,5]triazin-4(1H)-one:6-amino-5-nitropyridin-2-one (P:Z)] to see how the intermolecular interaction energies and their energetic components (electrostatics, exchange-repulsion, induction/polarization, and London dispersion interactions) vary among the base pairs. We examined the Hoogsteen (HG) geometries in addition to the traditional Watson-Crick (WC) geometries. Coupled-cluster theory through perturbative triples [CCSD(T)] extrapolated to the complete basis set (CBS) limit and high-order symmetry-adapted perturbation theory (SAPT) at the SAPT2+(3)(CCD)δMP2/aug-cc-pVTZ level are used to estimate highly accurate noncovalent interaction energies. Electrostatic interactions are the most attractive component of the interaction energies, but the sum of induction/polarization and London dispersion is nearly as large, for all base pairs and geometries considered. Interestingly, the non-natural Hachimoji base pairs interact more strongly than the corresponding natural base pairs, by -21.8 (B:S) and -0.3 (P:Z) kcal mol^-1 in the WC geometries, according to CCSD(T)/CBS. This is consistent with the H-bond distances being generally shorter in the non-natural base pairs. The natural base pairs are energetically more stabilized in their Hoogsteen geometries than in their WC geometries. The Hoogsteen geometry makes the A:T base pair slightly more stable, by -0.8 kcal mol^-1, and it greatly stabilizes the G:C⁺ base pair, by -15.3 kcal mol^-1. The G:C⁺ stabilization is mainly due to the fact that C has typically added a proton when found in Hoogsteen geometries. By contrast, Hoogsteen geometries are substantially less favorable than WC geometries for non-natural Hachimoji base pairs, by 17.3 (B:S) and 13.8 (P:Z) kcal mol^-1.

A computational study of cellulose regeneration: Coarse-grained molecular dynamics simulations

Understanding the microscopic mechanisms of regeneration of cellulose is prerequisite for engineering and controlling its material properties. In this paper, we performed coarse-grained Martini 3 molecular dynamics simulations of cellulose regeneration at a scale comparable to the experiments. The X-ray diffraction (XRD) curves were monitored to follow the structural changes of regenerated cellulose and trace formation of cellulose sheets and crystallites. The calculated coarse-grained morphologies of regenerated cellulose were backmapped to atomistic ones. After the backmapping we find that the regenerated coarse-grained cellulose structures calculated for both topology parameters of cellulose I_β and cellulose II/III, are transformed to cellulose II, where the calculated XRD curves exhibit the main peak at approximately 20-21 degrees, corresponding to the (110)/(020) planes of cellulose II. This result is in good quantitative agreement with the available experimental observations.

Discovery, implications and initial use of semi-synthetic organisms with an expanded genetic alphabet/code

Much recent interest has focused on developing proteins for human use, such as in medicine. However, natural proteins are made up of only a limited number of canonical amino acids with limited functionalities, and this makes the discovery of variants with some functions difficult. The ability to recombinantly express proteins containing non-canonical amino acids (ncAAs) with properties selected to impart the protein with desired properties is expected to dramatically improve the discovery of proteins with different functions. Perhaps the most straightforward approach to such an expansion of the genetic code is through expansion of the genetic alphabet, so that new codon/anticodon pairs can be created to assign to ncAAs. In this review, I briefly summarize more than 20 years of effort leading ultimately to the discovery of synthetic nucleotides that pair to form an unnatural base pair, which when incorporated into DNA, is stably maintained, transcribed and used to translate proteins in Escherichia coli. In addition to discussing wide ranging conceptual implications, I also describe ongoing efforts at the pharmaceutical company Sanofi to employ the resulting 'semi-synthetic organisms' or SSOs, for the production of next-generation protein therapeutics. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

Transfer learning for chemically accurate interatomic neural network potentials

Developing machine learning-based interatomic potentials from ab initio electronic structure methods remains a challenging task for computational chemistry and materials science. This work studies the capability of transfer learning, in particular discriminative fine-tuning, for efficiently generating chemically accurate interatomic neural network potentials on organic molecules from the MD17 and ANI data sets. We show that pre-training the network parameters on data obtained from density functional calculations considerably improves the sample efficiency of models trained on more accurate ab initio data. Additionally, we show that fine-tuning with energy labels alone can suffice to obtain accurate atomic forces and run large-scale atomistic simulations, provided a well-designed fine-tuning data set. We also investigate possible limitations of transfer learning, especially regarding the design and size of the pre-training and fine-tuning data sets. Finally, we provide GM-NN potentials pre-trained and fine-tuned on the ANI-1x and ANI-1ccx data sets, which can easily be fine-tuned on and applied to organic molecules.

Accurate Hellmann-Feynman forces from density functional calculations with augmented Gaussian basis sets

The Hellmann-Feynman (HF) theorem provides a way to compute forces directly from the electron density, enabling efficient force calculations for large systems through machine learning (ML) models for the electron density. The main issue holding back the general acceptance of the HF approach for atom-centered basis sets is the well-known Pulay force which, if naively discarded, typically constitutes an error upward of 10 eV/Å in forces. In this work, we demonstrate that if a suitably augmented Gaussian basis set is used for density functional calculations, the Pulay force can be suppressed, and HF forces can be computed as accurately as analytical forces with state-of-the-art basis sets, allowing geometry optimization and molecular dynamics to be reliably performed with HF forces. Our results pave a clear path forward for the accurate and efficient simulation of large systems using ML densities and the HF theorem.

Force-Field-Dependent DNA Breathing Dynamics: A Case Study of Hoogsteen Base Pairing in A6-DNA

The Hoogsteen (HG) base pairing conformation, commonly observed in damaged and mutated DNA helices, facilitates DNA repair and DNA recognition. The free energy difference between HG and Watson-Crick (WC) base pairs has been computed in previous studies. However, the mechanism of the conformational transition is not well understood. A detailed understanding of the process of WC to HG base pair transition can provide a deeper understanding of DNA repair and recognition. In an earlier study, we explored the free energy landscape for this process using extensive computer simulation with the CHARMM36 force field. In this work, we study the impact of force field models in describing the WC to HG base pairing transition using meta-eABF enhanced sampling, quasi-harmonic entropy calculation, and nonbonded energy analysis. The secondary structures of both base pairing forms and the topology of the free energy landscapes were consistent over different force field models, although the relative free energy, entropy, and the interaction energies tend to vary. The relative stability of the WC and HG conformations is dictated by a delicate balance between the enthalpic stabilization and the reduced entropy of the structurally rigid HG structure. These findings highlight the impact that subtleties in force field models can have on accurately modeling DNA base pair dynamics and should stimulate further computational investigations into other dynamically important motions in DNA.

Strategies for Controlling the Spatial Orientation of Single Molecules Tethered on DNA Origami Templates Physisorbed on Glass Substrates: Intercalation and Stretching

Nanoarchitectural control of matter is crucial for next-generation technologies. DNA origami templates are harnessed to accurately position single molecules; however, direct single molecule evidence is lacking regarding how well DNA origami can control the orientation of such molecules in three-dimensional space, as well as the factors affecting control. Here, we present two strategies for controlling the polar (θ) and in-plane azimuthal (ϕ) angular orientations of cyanine Cy5 single molecules tethered on rationally-designed DNA origami templates that are physically adsorbed (physisorbed) on glass substrates. By using dipolar imaging to evaluate Cy5's orientation and super-resolution microscopy, the absolute spatial orientation of Cy5 is calculated relative to the DNA template. The sequence-dependent partial intercalation of Cy5 is discovered and supported theoretically using density functional theory and molecular dynamics simulations, and it is harnessed as our first strategy to achieve θ control for a full revolution with dispersion as small as ±4.5°. In our second strategy, ϕ control is achieved by mechanically stretching the Cy5 from its two tethers, being the dispersion ±10.3° for full stretching. These results can in principle be applied to any single molecule, expanding in this way the capabilities of DNA as a functional templating material for single-molecule orientation control. The experimental and modeling insights provided herein will help engineer similar self-assembling molecular systems based on polymers, such as RNA and proteins.

Collision-induced dissociation of homodimeric and heterodimeric radical cations of 9-methylguanine and 9-methyl-8-oxoguanine: correlation between intra-base pair proton transfer originating from the N1-H at a Watson-Crick edge and non-statistical dissociation

It has been shown previously in protonated, deprotonated and ionized guanine-cytosine base pairs that intra-base pair proton transfer from the N1-H at the Watson-Crick edge of guanine to the complementary nucleobase prompts non-statistical dissociation of the base-pair system, and the dissociation of a proton-transferred base-pair structure is kinetically more favored than that of the starting, conventional base-pair structure. However, the fundamental chemistry underlying this anomalous and intriguing kinetics has not been completely revealed, which warrants the examination of more base-pair systems in different structural contexts in order to derive a generalized base-pair structure-kinetics correlation. The purpose of the present work is to expand the investigation to the non-canonical homodimeric and heterodimeric radical cations of 9-methylguanine (9MG) and 9-methyl-8-oxoguanine (9MOG), i.e., [9MG·9MG]˙⁺, [9MOG·9MG]˙⁺ and [9MOG·9MOG]˙⁺. Experimentally, collision-induced dissociation tandem mass spectrometry coupled with an electrospray ionization (ESI) source was used for the formation of base-pair radical cations, followed by detection of dissociation product ions and cross sections in the collisions with Xe gas under single ion-molecule collision conditions and as a function of the center-of-mass collision energy. Computationally, density functional theory and coupled cluster theory were used to calculate and identify probable base-pair structures and intra-base pair proton transfer and hydrogen transfer reactions, followed by kinetics modeling to explore the properties of dissociation transition states and kinetic factors. The significance of this work is twofold: it provides insight into base-pair opening kinetics in three biologically-important, non-canonical systems upon oxidative and ionization damage; and it links non-statistical dissociation to intra-base pair proton-transfer originating from the N1-H at the Watson-Crick edge of 8-oxoguanine, enhancing understanding towards the base-pair fragmentation assisted by proton transfer.

Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA-ellipticine and buckycatcher-fullerene

In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA-ellipticine intercalation complex, and the buckycatcher-C₆₀ complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller-Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA-ellipticine, and buckycatcher-C₆₀, respectively. Such sizable discrepancy above "chemical accuracy" for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree-Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.

A unified computational view of DNA duplex, triplex, quadruplex and their donor-acceptor interactions

DNA can assume various structures as a result of interactions at atomic and molecular levels (e.g., hydrogen bonds, π–π stacking interactions, and electrostatic potentials), so understanding of the consequences of these interactions could guide development of ways to produce elaborate programmable DNA for applications in bio- and nanotechnology. We conducted advanced ab initio calculations to investigate nucleobase model structures by componentizing their donor-acceptor interactions. By unifying computational conditions, we compared the independent interactions of DNA duplexes, triplexes, and quadruplexes, which led us to evaluate a stability trend among Watson–Crick and Hoogsteen base pairing, stacking, and even ion binding. For a realistic solution-like environment, the influence of water molecules was carefully considered, and the potassium-ion preference of G-quadruplex was first analyzed at an ab initio level by considering both base-base and ion-water interactions. We devised new structure factors including hydrogen bond length, glycosidic vector angle, and twist angle, which were highly effective for comparison between computationally-predicted and experimentally-determined structures; we clarified the function of phosphate backbone during nucleobase ordering. The simulated tendency of net interaction energies agreed well with that of real world, and this agreement validates the potential of ab initio study to guide programming of complicated DNA constructs.

DNA Ring-Opening Polymerization Driven by Base Stacking

Supramolecular polymers, relying on reversible intermolecular interactions, promise a wide range of applications, including optoelectronic materials, self-healing materials, and biomedical delivery materials. Among potential molecular candidates, DNA strands act as an excellent platform. DNA has a well-established secondary structure (double helix), and its intermolecular interactions can be readily thermodynamically engineered and kinetically controlled. Extensive studies have demonstrated that various DNA motifs can polymerize/assemble into large polymers with different topology, geometry, and dimensionalities. Most of the reported polymerization is driven by hybridization of DNA strands. Herein, we report a novel system of DNA supramolecular polymerization that is driven by DNA base stacking. The polymerization has been confirmed by native polyacrylamide gel electrophoresis (PAGE) and atomic force microscopy (AFM). We believe that this work will expand the toolbox for DNA supramolecular polymerization and would, with further development, increase further control of DNA supramolecular polymerization.

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.

Frequency and hydrogen bonding of nucleobase homopairs in small molecule crystals

We used the high resolution and accuracy of the Cambridge Structural Database (CSD) to provide detailed information regarding base pairing interactions of selected nucleobases. We searched for base pairs in which nucleobases interact with each other through two or more hydrogen bonds and form more or less planar structures. The investigated compounds were either free forms or derivatives of adenine, guanine, hypoxanthine, thymine, uracil and cytosine. We divided our findings into categories including types of pairs, protonation patterns and whether they are formed by free bases or substituted ones. We found base pair types that are exclusive to small molecule crystal structures, some that can be found only in RNA containing crystal structures and many that are native to both environments. With a few exceptions, nucleobase protonation generally followed a standard pattern governed by pKa values. The lengths of hydrogen bonds did not depend on whether the nucleobases forming a base pair were charged or not. The reasons why particular nucleobases formed base pairs in a certain way varied significantly.

The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation

In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.

Ab Initio Prediction of the Phase Transition for Solid Ammonia at High Pressures

Ammonia is one of the most basic components on the planet and its high-pressure characteristics play an important role in planetary science. Solid ammonia crystals frequently adopt multiple distinct polymorphs exhibiting different properties. Predicting the crystal structure of these polymorphs and under what thermodynamic conditions these polymorphs are stable would be of great value to environmental industry and other fields. Theoretical calculations based on the classical force fields and density-functional theory (DFT) are versatile methods but lack of accurate description of weak intermolecular interactions for molecular crystals. In this study, we employ an ab initio computational study on the solid ammonia at high pressures, using the second-order Møller-Plesset perturbation (MP2) theory and the coupled cluster singles, doubles, and perturbative triples (CCSD(T)) theory along with the embedded fragmentation method. The proposed algorithm is capable of performing large-scale calculations using high-level wavefunction theories, and accurately describing covalent, ionic, hydrogen bonding, and dispersion interactions within molecular crystals, and therefore can predict the crystal structures, Raman spectra and phase transition of solid ammonia phases I and IV accurately. We confirm the crystal structures of solid ammonia phases I and IV that have been controversial for a long time and predict their phase transition that occurs at 1.17 GPa and 210 K with small temperature dependence, which is in line with experiment.

The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules

Maximum diversification of data is a central theme in building generalized and accurate machine learning (ML) models. In chemistry, ML has been used to develop models for predicting molecular properties, for example quantum mechanics (QM) calculated potential energy surfaces and atomic charge models. The ANI-1x and ANI-1ccx ML-based general-purpose potentials for organic molecules were developed through active learning; an automated data diversification process. Here, we describe the ANI-1x and ANI-1ccx data sets. To demonstrate data diversity, we visualize it with a dimensionality reduction scheme, and contrast against existing data sets. The ANI-1x data set contains multiple QM properties from 5 M density functional theory calculations, while the ANI-1ccx data set contains 500 k data points obtained with an accurate CCSD(T)/CBS extrapolation. Approximately 14 million CPU core-hours were expended to generate this data. Multiple QM calculated properties for the chemical elements C, H, N, and O are provided: energies, atomic forces, multipole moments, atomic charges, etc. We provide this data to the community to aid research and development of ML models for chemistry.

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.

Binding of Sulfoxaflor to Aplysia californica-AChBP: Computational Insights from Multiscale Approaches

Structural features and binding properties of sulfoxaflor (SFX) with Ac-AChBP, the surrogate of the insect nAChR ligand binding domain (LBD), are reported herein using various complementary molecular modeling approaches (QM, molecular docking, molecular dynamics, and QM/QM'). The different SFX stereoisomers show distinct behaviors in terms of binding and interactions with Ac-AChBP. Molecular docking and Molecular Dynamics (MD) simulations highlight the specific intermolecular contacts involved in the binding of the different SFX isomers and the relative contribution of the SFX functional groups. QM/QM' calculations provide further insights and a significant refinement of the geometric and energetic contributions of the various residues leading to a preference for the SS and RR stereoisomers. Notable differences in terms of binding interactions are pointed out for the four stereoisomers. The results point out the induced fit of the Ac-AChBP binding site according to the SFX stereoisomer. In this process, the water molecules-mediated contacts play a key role, their energetic contribution being among the most important for the various stereoisomers. In all cases, the interaction with Trp147 is the major binding component, through CH···π and π···π interactions. This study provides a rationale for the binding of SFX to insect nAChR, in particular with respect to the new class of sulfoximine-based insect nAChR competitive modulators, and points out the requirements of various levels of theory for an accurate description of ligand-receptor interactions.

Approaching coupled cluster accuracy with a general-purpose neural network potential through transfer learning

Computational modeling of chemical and biological systems at atomic resolution is a crucial tool in the chemist's toolset. The use of computer simulations requires a balance between cost and accuracy: quantum-mechanical methods provide high accuracy but are computationally expensive and scale poorly to large systems, while classical force fields are cheap and scalable, but lack transferability to new systems. Machine learning can be used to achieve the best of both approaches. Here we train a general-purpose neural network potential (ANI-1ccx) that approaches CCSD(T)/CBS accuracy on benchmarks for reaction thermochemistry, isomerization, and drug-like molecular torsions. This is achieved by training a network to DFT data then using transfer learning techniques to retrain on a dataset of gold standard QM calculations (CCSD(T)/CBS) that optimally spans chemical space. The resulting potential is broadly applicable to materials science, biology, and chemistry, and billions of times faster than CCSD(T)/CBS calculations.

Disruption of nucleobase stacking to restore reactivity

Strong intermolecular interaction can prevent an organic molecule from dissolving in a reaction solution, thereby jeopardizing its reactivity and usefulness. Nucleobases and nucleosides (especially many purines and their derivatives) are notoriously difficult to dissolve in most organic solvents, generally attributed to their strong intermolecular interactions caused by the aromaticity, polarity and hydrogen-bonding. Guided by our computational study and prediction, to address this challenge, we have found that by doping the reaction solution with toluene (an inert aromatic compound), the added solvent molecules are capable of generating the stacking interaction with the solute molecules (e.g., purine derivatives) and disrupting the intermolecular stacking of the solute molecules. Thus, this inert doping can successfully address the insoluble challenge, dissolve the poorly soluble reactants (such as purine phosphoramidites), and restore the amidite reactivity for oligonucleotide synthesis. Our research has offered a simple strategy to efficiently synthesize labile oligonucleotides, via disrupting stacking interaction with inert aromatic molecules.

Does the stability of the stacking motif surpass the planar motif in 2-amino-4-nitrophenol? - a CCSD(T) analysis

In this work we analyzed O-H...O, O-H...N, and N-H...O contacts existing in the 2-amino-4-nitrophenol structure engaged in ANP molecules through quantum chemical methods. Furthermore, the above contacts were favored to comprehensively understand the stability of noncovalent interactions, π stacking and hydrogen bonding, surviving in 2-amino-4-nitrophenol. The geometries of π stacking and hydrogen bond interactions between two 2-amino-4-nitrophenols were optimized at BLYP-D3/def2-QZVP with dispersion 3 and MP2/cc-pVTZ levels of theory, and their stability was compared using the CCSD(T) interaction energies. The analyses predicted a particularly strong π stacking interaction of 2-amino-4-nitrophenol with hydrogen bond due to the narrow equivalent configuration of NO2 interactions with the other 2-amino-4-nitrophenols. Furthermore, this work focused on analyzing the stability of the individual hydrogen bonds existing in planar and stacked arrangements. Graphical abstract Stacked and planar motif in 2-amino-4-nitrophenol.

Hydrogen-Rich Cation Radicals of DNA Dinucleotides: Generation and Structure Elucidation by UV-Vis Action Spectroscopy

Hydrogen-rich DNA dinucleotide cation radicals (dGG + 2H)^+•, (dCG + 2H)^+•, and (dGC + 2H)^+• represent transient species comprising protonated and hydrogen atom adducted nucleobase rings that serve as models for proton and radical migrations in ionized DNA. These DNA cation radicals were generated in the gas phase by electron-transfer dissociation of dinucleotide dication-crown-ether complexes and characterized by UV-vis photodissociation action spectra, ab initio calculations of structures and relative energies, and time-dependent density functional theory calculations of UV-vis absorption spectra. Theoretical calculations indicate that (dGG + 2H)^+• cation radicals formed by electron transfer underwent an exothermic conformational collapse that was accompanied by guanine ring stacking and facile internucleobase hydrogen atom transfer, forming 3'-guanine C-8-H radicals. In contrast, exothermic hydrogen transfer from the 5'-cytosine radical onto the guanine ring in (dCG + 2H)^+• was kinetically hampered, resulting in the formation of a mixture of 5'-cytosine and 3'-guanine radicals. Conformational folding and nucleobase stacking were energetically unfavorable in (dGC + 2H)^+• that retained its structure of a 3'-cytosine radical, as formed by one-electron reduction of the dication. Hydrogen-rich guanine (G + H)^• and cytosine (C + H)^• radicals were calculated to have vastly different basicities in water, as illustrated by the respective p K_a values of 20.0 and 4.6, which is pertinent to their different abilities to undergo proton-transfer reactions in solution.

Assessment of Density Functional Methods for Geometry Optimization of Bimolecular van der Waals Complexes

We explore the suitability of three popular density functionals (B97-D3, B3LYP-D3, M05-2X) for producing accurate equilibrium geometries of van der Waals (vdW) complexes with diverse binding motifs. For these functionals, optimizations using Dunning's aug-cc-pVDZ basis set best combine accuracy and a reasonable computational expense. Each DFT/aug-cc-pVDZ combination produces optimized equilibrium geometries for 21 small vdW complexes of organic molecules (up to four non-hydrogen atoms total) that agree with high-level CCSD(T)/CBS reference geometries to within ±0.1 Å for the averages of the center-of-mass displacement and the mean least root-mean-squared displacement. The DFT/aug-cc-pVDZ combinations are also able to reproduce the optimal center-of-mass displacements interpolated from CCSD(T)/CBS radial potential energy surfaces in both NBC7x and HBC6 test sets to within ±0.1 Å. We therefore conclude that each of these denisty functional methods, together with the aug-cc-pVDZ basis set, is suitable for producing equilibrium geometries of generic nonbonded complexes.

An Electron Density Source-Function Study of DNA Base Pairs in Their Neutral and Ionized Ground States†

The source function (SF) decomposes the electron density at any point into contributions from all other points in the molecule, complex, or crystal. The SF "illuminates" those regions in a molecule that most contribute to the electron density at a point of reference. When this point of reference is the bond critical point (BCP), a commonly used surrogate of chemical bonding, then the SF analysis at an atomic resolution within the framework of Bader's Quantum Theory of Atoms in Molecules returns the contribution of each atom in the system to the electron density at that BCP. The SF is used to locate the important regions that control the hydrogen bonds in both Watson-Crick (WC) DNA dimers (adenine:thymine (AT) and guanine:cytosine (GC)) which are studied in their neutral and their singly ionized (radical cationic and anionic) ground states. The atomic contributions to the electron density at the BCPs of the hydrogen bonds in the two dimers are found to be delocalized to various extents. Surprisingly, gaining or loosing an electron has similar net effects on some hydrogen bonds concealing subtle compensations traced to atomic sources contributions. Coarser levels of resolutions (groups, rings, and/or monomers-in-dimers) reveal that distant groups and rings often have non-negligible effects especially on the weaker hydrogen bonds such as the third weak CH⋅⋅⋅O hydrogen bond in AT. Interestingly, neither the purine nor the pyrimidine in the neutral or ionized forms dominate any given hydrogen bond despite that the former has more atoms that can act as source or sink for the density at its BCP. © 2018 Wiley Periodicals, Inc.

Combining Accuracy and Efficiency: An Incremental Focal-Point Method Based on Pair Natural Orbitals

In this work, we present a new pair natural orbitals (PNO)-based incremental scheme to calculate CCSD(T) and CCSD(T0) reaction, interaction, and binding energies. We perform an extensive analysis, which shows small incremental errors similar to previous non-PNO calculations. Furthermore, slight PNO errors are obtained by using T_PNO = T_TNO with appropriate values of 10^-7 to 10^-8 for reactions and 10^-8 for interaction or binding energies. The combination with the efficient MP2 focal-point approach yields chemical accuracy relative to the complete basis-set (CBS) limit. In this method, small basis sets (cc-pVDZ, def2-TZVP) for the CCSD(T) part are sufficient in case of reactions or interactions, while some larger ones (e.g., (aug)-cc-pVTZ) are necessary for molecular clusters. For these larger basis sets, we show the very high efficiency of our scheme. We obtain not only tremendous decreases of the wall times (i.e., factors >10²) due to the parallelization of the increment calculations as well as of the total times due to the application of PNOs (i.e., compared to the normal incremental scheme) but also smaller total times with respect to the standard PNO method. That way, our new method features a perfect applicability by combining an excellent accuracy with a very high efficiency as well as the accessibility to larger systems due to the separation of the full computation into several small increments.

A DFT study of 2-aminopurine-containing dinucleotides: prediction of stacked conformations with B-DNA structure

The fluorescence properties of dinucleotides incorporating 2-aminopurine (2AP) suggest that the simplest oligonucleotides adopt conformations similar to those found in duplex DNA. However, there is a lack of structural data for these systems. We report a density functional theory (DFT) study of the structures of 2AP-containing dinucleotides (deoxydinucleoside monophosphates), including full geometry optimisation of the sugar-phosphate backbone. Our DFT calculations employ the M06-2X functional for reliable treatment of dispersion interactions and include implicit aqueous solvation. Dinucleotides with 2AP in the 5'-position and each of the natural bases in the 3'-position are examined, together with the analogous 5'-adenine-containing systems. Computed structures are compared in detail with typical B-DNA base-step parameters, backbone torsional angles and sugar pucker, derived from crystallographic data. We find that 2AP-containing dinucleotides adopt structures that closely conform to B-DNA in all characteristic parameters. The structures of 2AP-containing dinucleotides closely resemble those of their adenine-containing counterparts, demonstrating the fidelity of 2AP as a mimic of the natural base. As a first step towards exploring the conformational heterogeneity of dinucleotides, we also characterise an imperfectly stacked conformation and one in which the bases are completely unstacked.

Noncovalent Interactions in the Catechol Dimer

Noncovalent interactions play a significant role in a wide variety of biological processes and bio-inspired species. It is, therefore, important to have at hand suitable computational methods for their investigation. In this paper, we report on the contribution of dispersion and hydrogen bonds in both stacked and T-shaped catechol dimers, with the aim of delineating the respective role of these classes of interactions in determining the most stable structure. By using second-order Møller⁻Plesset (MP2) calculations with a small basis set, specifically optimized for these species, we have explored a number of significant sections of the interaction potential energy surface and found the most stable structures for the dimer, in good agreement with the highly accurate, but computationally more expensive coupled cluster single and double excitation and the perturbative triples (CCSD(T))/CBS) method.

Intermolecular interactions of oligothienoacenes: Do S⋯S interactions positively contribute to crystal structures of sulfur-containing aromatic molecules?

Intermolecular interactions in the crystals of tetra- and penta-thienoacene were studied using ab initio molecular orbital calculations for evaluating the magnitude of characteristic S⋯S interactions with great attention paid to their origin. The interactions between the π-stacked neighboring molecules are significantly greater than those between the neighboring molecules exhibiting the S⋯S contact, although it has sometimes been claimed that the S⋯S interactions play important roles in adjusting the molecular arrangement of sulfur-containing polycyclic aromatic molecules in the crystals owing to short S⋯S contacts. The coupled cluster calculations with single and double substitutions with noniterative triple excitation interaction energies at the basis set limit estimated for the π-stacked and S⋯S contacted neighboring molecules in the tetrathienoacene crystal are -11.17 and -4.27 kcal/mol, respectively. Those for π-stacked molecules in the pentathienoacene crystal is -14.38 kcal/mol, while those for S⋯S contacted molecules are -7.02 and -6.74 kcal/mol. The dispersion interaction is the major source of the attraction between the π-stacked and S⋯S contacted molecules, while the orbital-orbital interactions are repulsive: The orbital-orbital interactions, which are significant for charge carrier transport properties, are not much more than the results of the short S⋯S contact caused by the strong dispersion interactions. Besides, the intermolecular interaction energy calculated for a trithienoacene dimer has strong orientation dependence.

Programmed Transfer of Sequence Information into a Molecularly Imprinted Polymer for Hexakis(2,2'-bithien-5-yl) DNA Analogue Formation toward Single-Nucleotide-Polymorphism Detection

A new strategy of simple, inexpensive, rapid, and label-free single-nucleotide-polymorphism (SNP) detection using robust chemosensors with piezomicrogravimetric, surface plasmon resonance, or capacitive impedimetry (CI) signal transduction is reported. Using these chemosensors, selective detection of a genetically relevant oligonucleotide under FIA conditions within 2 min is accomplished. An invulnerable-to-nonspecific interaction molecularly imprinted polymer (MIP) with electrochemically synthesized probes of hexameric 2,2'-bithien-5-yl DNA analogues discriminating single purine-nucleobase mismatch at room temperature was used. With density functional theory modeling, the synthetic procedures developed, and isothermal titration calorimetry quantification, adenine (A)- or thymine (T)-substituted 2,2'-bithien-5-yl functional monomers capable of Watson-Crick nucleobase pairing with the TATAAA oligodeoxyribonucleotide template or its peptide nucleic acid (PNA) analogue were designed. Characterized by spectroscopic techniques, molecular cavities exposed the ordered nucleobases on the 2,2'-bithien-5-yl polymeric backbone of the TTTATA hexamer probe designed to hybridize the complementary TATAAA template. In that way, an artificial TATAAA-promoter sequence was formed in the MIP. The purine nucleobases of this sequence are known to be recognized by RNA polymerase to initiate the transcription in eukaryotes. The hexamer strongly hybridized TATAAA with the complex stability constant K_s^{TTTATA-TATAAA} = k_a/k_d ≈ 10⁶ M^-1, as high as that characteristic for longer-chain DNA-PNA hybrids. The CI chemosensor revealed a 5 nM limit of detection, quite appreciable as for the hexadeoxyribonucleotide. Molecular imprinting increased the chemosensor sensitivity to the TATAAA analyte by over 4 times compared to that of the nonimprinted polymer. The herein-devised detection platform enabled the generation of a library of hexamer probes for typing the majority of SNP probes as well as studying a molecular mechanism of the complex transcription machinery, physics of single polymer molecules, and stable genetic nanomaterials.