The elusive atomic rationale for DNA base pair stability

Multiple hydrogen-bonded dimers: are only the frontier atoms relevant?

Non-frontier atom exchanges in hydrogen-bonded aromatic dimers can induce significant interaction energy changes (up to 6.5 kcal mol-1). Our quantum-chemical analyses reveal that the relative hydrogen-bond strengths of N-edited guanine-cytosine base pair isosteres, which cannot be explained from the frontier atoms, follow from the charge accumulation in the monomers.

Stability of Quadruple Hydrogen Bonds in an Ionic Liquid Environment

Hydrogen bonds (H-bonds) are highly sensitive to the surrounding environments owing to their dipolar nature, with polar solvents kown to significantly weaken H-bonds. Herein, the stability of the H-bonding motif ureidopyrimidinone (UPy) is investigated, embedded into a highly polar polymeric ionic liquid (PIL) consisting of pendant pyrrolidinium bis(trifluoromethylsulfonyl)imide (IL) moieties, to study the influence of such ionic environments on the UPy H-bonds. The content of the surrounding IL is changed by addition of an additional low molecular weight IL to further boost the IL content around the UPy moieties in molar ratios of UPy/IL ranging from 1/4 up to 1/113, thereby promoting the polar microenvironment around the UPy-H-bonds. Variable-temperature solid-state MAS NMR spectroscopy and FT-IR spectroscopy demonstrate that the UPy H-bonds are largely present as (UPy-) dimers, but sensitive to elevated temperatures (>70 °C). Subsequent rheology and DSC studies reveal that the ILs only solvate the polymeric chains but do not interfere with the UPy-dimer H-bonds, thus accounting for their high stability and applicability in many material systems.

Quantification of secondary electrostatic interactions in H-bonded complexes

The H-bonding properties of compounds that contain multiple functional groups are difficult to predict, because there are through-bond polarisation effects and long-range secondary electrostatic interactions that have significant effects on the interactions with solvents and other molecules. Here we use experimental measurements of association constants for formation of 1 : 1 H-bonded complexes that contain a single well-defined H-bond and a single well-defined secondary electrostatic interaction to quantify the magnitude of this effect. The results were used to develop a computational method for calculating functional group H-bond parameters that accurately reproduce the magnitudes of both primary H-bonding interaction and secondary electrostatic interactions. The effects of secondary electrostatic interactions are observed in calculations of ab initio Molecular Electrostatic Potential (MEP) values, but at the van der Waals surface, the magnitude of the effect is highly overestimated. MEP values calculated on electron density isosurfaces that lie closer to the nuclei provide a more accurate description of the experimental observations. H-bond parameters calculated using this approach successfully account for the properties of arrays of multiple H-bond donor and acceptor groups in different configurations. The results provide insight into the factors that govern the interaction properties of molecules that contain multiple functional groups and provide an accurate method for prediction of solution phase complexation free energies based on gas phase calculations of individual molecules.

Rare Tautomers of Artificially Expanded Genetic Letters and their Effects on the Base pair Stabilities

To expand the existing genetic letters, it is necessary to design robust nucleotides that can function naturally in living cells. Therefore, it is desirable to examine the roles of recently proposed second-generation artificially expanded genetic letters in producing stable duplex DNA. Here, a reliable dispersion-corrected density functional theory method is used to understand the electronic structures and properties of different rare tautomers of proposed expanded genetic letters and their effects on the base pair stabilities in the duplex DNA. It is found that the rare tautomers are not only stable in the aqueous medium but can also base pair with natural bases to produce stable mispairs. Except for J and V, all the artificial genetic letters are found to produce mispairs that are about 1-7 kcal/mol more stable than their complementary counterparts. They are also appreciably more stable than the naturally occurring G:C, A:T, and G:T pairs. The higher base pair stabilities are found to be mainly because of the polarity of monomers and attractive electrostatic interactions.

Hydrogen bond design principles

Hydrogen bonding principles are at the core of supramolecular design. This overview features a discussion relating molecular structure to hydrogen bond strengths, highlighting the following electronic effects on hydrogen bonding: electronegativity, steric effects, electrostatic effects, π-conjugation, and network cooperativity. Historical developments, along with experimental and computational efforts, leading up to the birth of the hydrogen bond concept, the discovery of nonclassical hydrogen bonds (C-H…O, O-H…π, dihydrogen bonding), and the proposal of hydrogen bond design principles (e.g., secondary electrostatic interactions, resonance-assisted hydrogen bonding, and aromaticity effects) are outlined. Applications of hydrogen bond design principles are presented.

The Nature of Hydrogen Bonds: A Delineation of the Role of Different Energy Components on Hydrogen Bond Strengths and Lengths

Hydrogen bonds are a complex interplay between different energy components, and their nature is still subject of an ongoing debate. In this minireview, we therefore provide an overview of the different perspectives on hydrogen bonding. This will be done by discussing the following individual energy components: 1) electrostatic interactions, 2) charge-transfer interactions, 3) π-resonance assistance, 4) steric repulsion, 5) cooperative effects, 6) dispersion interactions and 7) secondary electrostatic interactions. We demonstrate how these energetic factors are essential in a correct description of the hydrogen bond, and discuss several examples of systems whose energetic and geometrical features are not captured by easy-to-use predictive models.

Secondary Electrostatic Interaction Model Revised: Prediction Comes Mainly from Measuring Charge Accumulation in Hydrogen-Bonded Monomers

The secondary electrostatic interaction (SEI) model is often used to predict and explain relative hydrogen bond strengths of self-assembled systems. The SEI model oversimplifies the hydrogen-bonding mechanisms by viewing them as interacting point charges, but nevertheless experimental binding strengths are often in line with the model’s predictions. To understand how this rudimentary model can be predictive, we computationally studied two tautomeric quadruple hydrogen-bonded systems, DDAA-AADD and DADA-ADAD. Our results reveal that when the proton donors D (which are electron-donating) and the proton acceptors A (which are electron-withdrawing) are grouped together as in DDAA, there is a larger accumulation of charge around the frontier atoms than when the proton donor and acceptor groups are alternating as in DADA. This accumulation of charge makes the proton donors more positive and the proton acceptors more negative, which enhances both the electrostatic and covalent interactions in the DDAA dimer. The SEI model is thus predictive because it provides a measure for the charge accumulation in hydrogen-bonded monomers. Our findings can be understood from simple physical organic chemistry principles and provide supramolecular chemists with meaningful understanding for tuning hydrogen bond strengths and thus for controlling the properties of self-assembled systems.

Stability of doubly and triply H-bonded complexes governed by acidity-basicity relationships

We used our recently proposed acidity-basicity interplay (ABI) model (Chem. Sci., 2018, 9, 4402) and the Jorgensen secondary interactions hypothesis (JSIH) to rationalise the experimentally observed trends in the formation constants of doubly and triply H-bonded systems with -NHO[double bond, length as m-dash]C- and -NHN- interactions. Unlike the JSIH, the ABI interpretation can explain the trends in the complexation of amide/imide homo- and heterodimers as well as ADA-DAD clusters. We found that the strongest H-bonds play a very important role, a condition which offers an alternative to the well established JSIH to modulate the stability of these relevant systems.

A relative energy gradient (REG) study of the planar and perpendicular torsional energy barriers in biphenyl

Biphenyl is a prototype molecule, the study of which is important for a proper understanding of stereo-electronic effects. In the gas phase it has an equilibrium central torsion angle of ~ 45° and shows both a planar (0°) and a perpendicular (90°) torsional energy barrier. The latter is analysed for the first time. We use the newly proposed REG method, which is an exhaustive procedure that automatically ranks atomic energy contributions according to their importance in explaining the energy profile of a total system. Here, the REG method operates on energy contributions computed by the interacting quantum atoms method. This method is minimal in architecture and provides a crisp picture of well-defined and well-separated electrostatic, steric and exchange (covalent) energies at atomistic level. It is shown that the bond critical point occurring between the ortho-hydrogens in the planar geometry has been wrongly interpreted as a sign of repulsive interaction. A convenient metaphor of analysing football matches is introduced to clarify the role of a REG analysis.

Using the Relative Energy Gradient Method with Interacting Quantum Atoms to Determine the Reaction Mechanism and Catalytic Effects in the Peptide Hydrolysis in HIV-1 Protease

The reaction mechanism in an active site is of the utmost importance when trying to understand the role that an enzyme plays in biological processes. In a recently published paper [Theor. Chem. Acc. 2017, 136, 86], we formalised the Relative Energy Gradient (REG) method for automating an Interacting Quantum Atoms (IQA) analysis. Here, the REG method is utilised to determine the mechanism of peptide hydrolysis in the aspartic active site of the enzyme HIV-1 Protease. Using the REG method along with the IQA approach we determine the mechanism of peptide hydrolysis without employing any arbitrary parameters and with remarkable ease (albeit at large computational cost: the system contains 133 atoms, which means that there are 17 689 individual IQA terms to be calculated). When REG and IQA work together it is possible to determine a reaction mechanism at atomistic resolution from data directly derived from quantum calculations, without arbitrary parameters. Moreover, the mechanism determined by this novel method gives concrete insight into how the active site residues catalyse peptide hydrolysis.

Acidity and basicity interplay in amide and imide self-association

Simple acid–base properties explain the differences in amide and imide dimerisation, and represent an alternative to the secondary interactions hypothesis.

Amides dimerise more strongly than imides despite their lower acidity. Such an unexpected result has been rationalised in terms of the Jorgensen Secondary Interactions Hypothesis (JSIH) that involves the spectator (C Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O_S) and H-bonded (C Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O_HB) carbonyl groups in imides. Notwithstanding the considerable body of experimental and theoretical evidence supporting the JSIH, there are some computational studies which suggest that there might be other relevant intermolecular interactions than those considered in this model. We conjectured that the spectator carbonyl moieties could disrupt the resonance-assisted hydrogen bonds in imide dimers, but our results showed that this was not the case. Intrigued by this phenomenon, we studied the self-association of a set of amides and imides via¹H-NMR, ¹H-DOSY experiments, DFT calculations, QTAIM topological analyses of the electron density and IQA partitions of the electronic energy. These analyses revealed that there are indeed repulsions of the type O_S···O_HB in accordance with the JSIH but our data also indicate that the C Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O_S group has an overall attraction with the interacting molecule. Instead, we found correlations between self-association strength and simple Brønsted–Lowry acid/base properties, namely, N–H acidities and C Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O basicities. The results in CDCl₃ and CCl₄ indicate that imides dimerise less strongly than structurally related amides because of the lower basicity of their carbonyl fragments, a frequently overlooked aspect in the study of H-bonding. Overall, the model proposed herein could provide important insights in diverse areas of supramolecular chemistry such as the study of multiple hydrogen-bonded adducts which involve amide or imide functional groups.

Aromaticity gain increases the inherent association strengths of multipoint hydrogen-bonded arrays

Besides the textbook checklist for ways to control the association strengths of multipoint hydrogen-bonded arrays, “aromaticity” also makes the list!

The ANANKE relative energy gradient (REG) method to automate IQA analysis over configurational change

Much chemical insight ultimately comes down to finding out which fragment of a total system behaves like the total system, in terms of an energy profile. A simple example is that of the water dimer, where this system is regarded as held together by a hydrogen bond. The hydrogen bond consists of two atoms (H···O), which energetically behave similarly to the total system (H₂O)₂. However, from a quantum mechanical point of view, each atom in the total system interacts with any other atom. Thus, the view that the hydrogen bond by itself governs the energetic stability of the water dimer needs rigorous justification. In this work, we propose a method that provides such a justification, in general, but only illustrated on the water dimer here. This method is based on the topological energy partitioning method called interacting quantum atoms (IQA). The method is implemented in the program ANANKE, which calculates correlations between the energy profile of the total system and those of subsystems (or fragments). ANANKE acts on the IQA energy contributions obtained for a sequence of full-system geometries controlled by a coordinate of interest (e.g. the O···H distance in the water dimer). Although applied only for the water dimer in this work, the method is general and able to explain the gauche effect, the torsional barrier in biphenyl, the arrow-pushing scheme of an enzymatic reaction (peptide hydrolysis in the HIV-1 Protease active site), and halogen-alkane nucleophilic substitution (S_N2) reactions. Those applications will appear elsewhere as separate and elaborated case studies; here we focus on the details of the ANANKE method and its justification, using the water dimer as a concrete case.

Exploiting directional long range secondary forces for regulating electrostatics-dominated noncovalent interactions

It has been well established that long range secondary electrostatic interactions (SEIs) have a significant effect on the stability of supramolecular complexes. However, general rules for exploiting SEIs in the rational design of diverse supramolecular complexes have been difficult to obtain. In this work, we outline a quantum chemical approach for understanding the strength of electrostatic interactions. This approach is seen to provide excellent correlation between the electrostatic force and the binding energy between two partners in hydrogen-bonded complexes, as well as that between two ions in ion-pair complexes. Furthermore, we illustrate how the understanding of the binding allows for the rational design of new complexes where the association constant between the two partners can be increased or decreased, as desired, by several orders of magnitude. Hence, the current work showcases a general, simple and powerful method of understanding and exploiting long range secondary electrostatic interactions.

Roles of the amino group of purine bases in the thermodynamic stability of DNA base pairing

The energetic aspects of hydrogen-bonded base-pair interactions are important for the design of functional nucleotide analogs and for practical applications of oligonucleotides. The present study investigated the contribution of the 2-amino group of DNA purine bases to the thermodynamic stability of oligonucleotide duplexes under different salt and solvent conditions, using 2'-deoxyriboinosine (I) and 2'-deoxyribo-2,6-diaminopurine (D) as non-canonical nucleotides. The stability of DNA duplexes was changed by substitution of a single base pair in the following order: G•C > D•T ≈ I•C > A•T > G•T > I•T. The apparent stabilization energy due to the presence of the 2-amino group of G and D varied depending on the salt concentration, and decreased in the water-ethanol mixed solvent. The effects of salt concentration on the thermodynamics of DNA duplexes were found to be partially sequence-dependent, and the 2-amino group of the purine bases might have an influence on the binding of ions to DNA through the formation of a stable base-paired structure. Our results also showed that physiological salt conditions were energetically favorable for complementary base recognition, and conversely, low salt concentration media and ethanol-containing solvents were effective for low stringency oligonucleotide hybridization, in the context of conditions employed in this study.

On the physical role of exchange in the formation of an intramolecular bond path between two electronegative atoms

In this paper, we present a detailed energetic decomposition of intramolecular O···X interactions (X being O, S, or a halogen atom) based on the interacting quantum atoms approach of Pendás and co-workers. The nature of these interactions (repulsive or attractive, more or less electrostatic) is discussed in the framework of Bader's atoms in molecules theory, a particular emphasis being put on delocalization (measured by delocalization indexes and in terms of the source function) and on the exchange contributions. Notably, the concept of exchange channels introduced by Pendás and collaborators provides means of rationalizing and predicting the presence of bond critical points, enhancing the physical meaning of bond paths.

An AAAA–DDDD quadruple hydrogen-bond array

Secondary electrostatic interactions between adjacent hydrogen bonds can have a significant effect on the stability of a supramolecular complex. In theory, the binding strength should be maximized if all the hydrogen-bond donors (D) are on one component and all the hydrogen-bond acceptors (A) are on the other. Here, we describe a readily accessible AAAA–DDDD quadruple hydrogen-bonding array that exhibits exceptionally strong binding for a small-molecule hydrogen-bonded complex in a range of different solvents (K(a) > 3 × 10(12) M(-1) in CH2Cl2, 1.5 × 10(6) M(-1) in CH3CN and 3.4 × 10(5) M(-1) in 10% v/v DMSO/CHCl3). The association constant in CH2Cl2 corresponds to a binding free energy (ΔG) in excess of –71 kJ mol(-1) (more than 20% of the thermodynamic stability of a carbon–carbon covalent bond), which is remarkable for a supramolecular complex held together by just four intercomponent hydrogen bonds.

AAA-DDD triple hydrogen bond complexes

Experiment and theory both suggest that the AAA-DDD pattern of hydrogen bond acceptors (A) and donors (D) is the arrangement of three contiguous hydrogen bonding centers that results in the strongest association between two species. Murray and Zimmerman prepared the first example of such a system (complex 3*2) and determined the lower limit of its association constant (K(a)) in CDCl(3) to be 10(5) M(-1) by (1)H NMR spectroscopy (Murray, T. J. and Zimmerman, S. C. J. Am. Chem. Soc. 1992, 114, 4010-4011). The first cationic AAA-DDD pair (3*4(+)) was described by Bell and Anslyn (Bell, D. A. and Anslyn, E. A. Tetrahedron 1995, 51, 7161-7172), with a K(a) > 5 x 10(5) M(-1) in CH(2)Cl(2) as determined by UV-vis spectroscopy. We were recently able to quantify the strength of a neutral AAA-DDD arrangement using a more chemically stable AAA-DDD system, 6*2, which has an association constant of 2 x 10(7) M(-1) in CH(2)Cl(2) (Djurdjevic, S., Leigh, D. A., McNab, H., Parsons, S., Teobaldi, G. and Zerbetto, F. J. Am. Chem. Soc. 2007, 129, 476-477). Here we report on further AA(A) and DDD partners, together with the first precise measurement of the association constant of a cationic AAA-DDD species. Complex 6*10(+)[B(3,5-(CF(3))(2)C(6)H(3))(4)(-)] has a K(a) = 3 x 10(10) M(-1) at RT in CH(2)Cl(2), by far the most strongly bound triple hydrogen bonded system measured to date. The X-ray crystal structure of 6*10(+) with a BPh(4)(-) counteranion shows a planar array of three short (NH...N distances 1.95-2.15 A), parallel (but staggered rather than strictly linear; N-H...N angles 165.4-168.8 degrees), primary hydrogen bonds. These are apparently reinforced, as theory predicts, by close electrostatic interactions (NH-*-N distances 2.78-3.29 A) between each proton and the acceptor atoms of the adjacent primary hydrogen bonds.

Static and dynamic descriptions of bond breaking/formation: A complementary view?

Ab initio molecular-dynamic simulations using density-functional theory and the recent atom-centered density-matrix propagation (ADMP) method were used to study the bond breaking and formation for a case-study substitution nucleophilic bimolecular reaction, namely, the Walden inversion. Using the atoms-in-molecule approach, we have performed a detailed analysis to investigate intra- and intermolecular charge transfer along the ADMP trajectory. These results were compared to those obtained considering a static approach, such as the intrinsic reaction path. In particular, the topological properties computed along the dynamic trajectory well evidence a stronger electron exchange tending to spontaneously maximize the rising covalent interaction. Furthermore, their analysis suggests that the bond formation mechanism involves a reactive intermediate with a bonding interaction stronger than in the final product.

Hydrogen Bonding and Covalent Effects in Binding of Cisplatin to Purine Bases: Ab Initio and Atoms in Molecules Studies

Ab initio and density functional calculations are employed to investigate the role of hydrogen bonding in the binding of cisplatin to the purine bases guanine and adenine. Through the use of the theory of atoms in molecules (AIM), it is shown that hydrogen bonds are ubiquitous in such systems, with N-H...N and N-H...Cl interactions present in addition to the expected N-H...O. This in turn means that the known stability of cisplatin-guanine complexes cannot be ascribed solely to hydrogen bonding and allows decomposition of total binding energy into contributions from covalent and hydrogen bonds. To do so, a new method for predicting hydrogen bond energies from bond critical point properties is proposed, employing partial least-squares analysis to remove the family dependence of simple models. Still more hydrogen bond motifs are found in bifunctional complexes of the general type purine-[Pt(NH(3))(2)](2+)-purine, including purine...purine contacts, though again the energetics of these are insufficient to explain the observed trends in stability. Finally, the effect of platination on the pairing of guanine with cytosine is studied in a similar manner, revealing large redistributions of hydrogen bonding but surprisingly small overall changes in pairing energy.