Potential Energy Surface of the Cytosine Dimer: MP2 Complete Basis Set Limit Interaction Energies, CCSD(T) Correction Term, and Comparison with the AMBER Force Field

Benchmarking two-body contributions to crystal lattice energies and a range-dependent assessment of approximate methods

Using the many-body expansion to predict crystal lattice energies (CLEs), a pleasantly parallel process, allows for flexibility in the choice of theoretical methods. Benchmark-level two-body contributions to CLEs of 23 molecular crystals have been computed using interaction energies of dimers with minimum inter-monomer separations (i.e., closest contact distances) up to 30 Å. In a search for ways to reduce the computational expense of calculating accurate CLEs, we have computed these two-body contributions with 15 different quantum chemical levels of theory and compared these energies to those computed with coupled-cluster in the complete basis set (CBS) limit. Interaction energies of the more distant dimers are easier to compute accurately and several of the methods tested are suitable as replacements for coupled-cluster through perturbative triples for all but the closest dimers. For our dataset, sub-kJ mol^-1 accuracy can be obtained when calculating two-body interaction energies of dimers with separations shorter than 4 Å with coupled-cluster with single, double, and perturbative triple excitations/CBS and dimers with separations longer than 4 Å with MP2.5/aug-cc-pVDZ, among other schemes, reducing the number of dimers to be computed with coupled-cluster by as much as 98%.

Perspectives on weak interactions in complex materials at different length scales

Nanocomposite materials consist of nanometer-sized quantum objects such as atoms, molecules, voids or nanoparticles embedded in a host material. These quantum objects can be exploited as a super-structure, which can be designed to create material properties targeted for specific applications. For electromagnetism, such targeted properties include field enhancements around the bandgap of a semiconductor used for solar cells, directional decay in topological insulators, high kinetic inductance in superconducting circuits, and many more. Despite very different application areas, all of these properties are united by the common aim of exploiting collective interaction effects between quantum objects. The literature on the topic spreads over very many different disciplines and scientific communities. In this review, we present a cross-disciplinary overview of different approaches for the creation, analysis and theoretical description of nanocomposites with applications related to electromagnetic properties.

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.

Probing the Nature of Noncovalent Interactions in Dimers of Linear Tyrosine-Based Dipeptides

Tyrosine-based dipeptides self-assemble to form higher order structures. To gain insights into the nature of intermolecular interactions contributing to the early stages of the self-assembly of aromatic dipeptides, we study the dimers of linear dityrosine (YY) and tryptophan-tyrosine (WY) using quantum-chemical methods with dispersion corrections and universal solvation model based on density in combination with energy decomposition and natural bond orbital (NBO) analyses. We find that hydrogen bonding is a dominant stabilizing force. The lowest energy structure for the linear YY dimer is characterized by Ocarboxyl···H(O)tyr. In contrast, the lowest energy dimer of linear WY is stabilized by Ocarboxyl···H(N)trp and πtyr···πtyr. The solvent plays a critical role as it can change the strength and nature of interactions. The lowest energy for linear WY dimer in acetone is stabilized by Ocarboxyl···H(O)tyr, πtrp···H(C), and πtrp···H(N). The ΔG of dimerization and stabilization energies of solvated dipeptides reveal that the dipeptide systems are more stable in the solvent phase than in gas phase. NBO confirms increased magnitudes for donor-acceptor interaction for the solvated dipeptides.

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.

Density Functionals for Noncovalent Interaction Energies of Biological Importance

Forty density functionals and one wavefunction method are assessed against a recently published database of accurate noncovalent interaction energies of biological importance. The comparison shows that two newly developed density functional theory (DFT) methods, PWB6K and M05-2X, give the best performance for this benchmark database of 22 noncovalent complexes, including both hydrogen-bonding and dispersion-dominated complexes. In contrast, the more popular B3LYP and PBEh functionals fail to describe the interactions in the dispersion-dominated complexes. The local spin density approximation and BHandH functionals give good performance for dispersion-dominated interactions at the expense of a large error for hydrogen bonding. PWB6K and M05-2X constitute a new generation of DFT methods based on simultaneously optimized exchange and correlation functionals that include kinetic energy density in both the exchange and correlation functional, and the present study confirms that they have greatly improved performance for noncovalent interactions as compared to previous DFT methods. We interpret this as being due to an improved treatment of medium-range correlation effects by the exchange-correlation functional. We recommend the PWB6K and M05-2X methods for investigating large biological systems and soft materials.

First principles potential for the cytosine dimer

A new first principles potential for the cytosine dimer.

DNA-protein π-interactions in nature: abundance, structure, composition and strength of contacts between aromatic amino acids and DNA nucleobases or deoxyribose sugar

Four hundred twenty-eight high-resolution DNA-protein complexes were chosen for a bioinformatics study. Although 164 crystal structures (38% of those searched) contained no interactions, 574 discrete π-contacts between the aromatic amino acids and the DNA nucleobases or deoxyribose were identified using strict criteria, including visual inspection. The abundance and structure of the interactions were determined by unequivocally classifying the contacts as either π-π stacking, π-π T-shaped or sugar-π contacts. Three hundred forty-four nucleobase-amino acid π-π contacts (60% of all interactions identified) were identified in 175 of the crystal structures searched. Unprecedented in the literature, 230 DNA-protein sugar-π contacts (40% of all interactions identified) were identified in 137 crystal structures, which involve C-H···π and/or lone-pair···π interactions, contain any amino acid and can be classified according to sugar atoms involved. Both π-π and sugar-π interactions display a range of relative monomer orientations and therefore interaction energies (up to -50 (-70) kJ mol(-1) for neutral (charged) interactions as determined using quantum chemical calculations). In general, DNA-protein π-interactions are more prevalent than perhaps currently accepted and the role of such interactions in many biological processes may yet to be uncovered.

Significant strength of charged DNA-protein π-π interactions: a preliminary study of cytosine

The present work characterized the preferred gas-phase structure and optimum interaction energy of both parallel stacked and perpendicular T-shaped dimers between cytosine (C), as a representative nucleobase, and aspartic/glutamic acid (DE), aspartate/glutamate (DE(-)) or arginine (R(+)), using detailed M06-2X/6-31+G(d,p) potential energy surface scans as a function of the relative monomer orientation. Through comparison to previous literature on the π-π interactions between the DNA nucleobases and the aromatic amino acid residues, this work will allow for comparisons between DNA-protein interactions involving aromatic and acyclic R-side chains, as well as comparisons of the relative geometric dependence and magnitude of π-π (C:DE), πcation-π (C:R(+)), and πanion-π (C:DE(-)) interactions. Our results show that the preferred relative monomer orientation is highly dependent on the monomer composition and charge, and is dictated by electrostatic-driven interactions. More importantly, for the first time, we report that the π-π interactions between cytosine and (neutral) aspartic/glutamic acid are up to approximately -40 kJ mol(-1), while the πcation-π or πanion-π interactions between cytosine and arginine or aspartate/glutamate are up to approximately -90 and -99 kJ mol(-1), respectively. An extensive investigation of the effects of the computational methodology implemented, including comparisons to detailed CCSD(T)/CBS potential energy surfaces and interaction energies, supports the use of M06-2X, as well as ωB97X-D, to study DNA-protein π-π interactions of varying composition and charge. Most importantly, the CCSD(T)/CBS results verify the strong nature of these DNA-protein π-π interactions, as well as the unique nature of the πcation-π and πanion-π counterparts. Therefore, our results emphasize that a wide variety of different types of noncovalent interactions between both cyclic and acyclic π-containing components can significantly contribute to the stability of DNA-protein complexes and likely play a larger role in biology than currently accepted.

Configuration of wobble base pairs having pyrimidines as anticodon wobble bases: significance for codon degeneracy

Degeneracy of the genetic code was attributed by Crick to imprecise hydrogen-bonded base-pairing at the wobble position during codon-anticodon pairing. The Crick wobble rules define but do not explain the RNA base pair combinations allowed at this position. We select six pyrimidine bases functioning as anticodon wobble bases (AWBs) to study their H-bonded pairing properties with the four major RNA bases using density functional theory at the B3LYP/6-31G(d,p) level. This is done to assess the extent to which the configuration of a solitary RNA wobble base pair may in itself determine specificity and degeneracy of the genetic code by allowing or disallowing the given base pair during codon-anticodon pairing. Calculated values of select configuration markers for the base pairs screen well between allowed and disallowed base pairs for most cases examined here, where the base pair width emerges as an important factor. A few allowed wobble pairs invoke the involvement of RNA nucleoside conformation, as well as involvement of the exocyclic substituent in H-bonding. This study, however, cannot explain the disallowed status of the Ura⋯Gua wobble pair on the basis of configuration alone. Explanation of the allowed status of the V⋯Ura pair requires further study on the mediatory role of water molecules. Apart from these two cases, these computational results are sufficient, on the basis of base pair configuration alone, to account for the specificity and degeneracy of the genetic code for all known cases of codon-anticodon pairing which involve the pyrimidine AWBs studied here.

yDNA versus yyDNA pyrimidines: computational analysis of the effects of unidirectional ring expansion on the preferred sugar-base orientation, hydrogen-bonding interactions and stacking abilities

The properties of natural, y- and yy-pyrimidines are compared using computational (B3LYP, MP2) methods. Ring expansion upon incorporation of benzene or naphthalene into the natural pyrimidines affects the preferred orientation of the base about the glycosidic bond in the corresponding nucleoside to a similar extent. Specifically, although the natural pyrimidines preferentially adopt the anti orientation with respect to the 2'-deoxyribose moiety, the expanded analogues will likely display (anti/syn) conformational heterogeneity, which may lead to alternate hydrogen-bonding modes in double-stranded duplexes. Nevertheless, the A:T Watson-Crick hydrogen-bond strengths do not significantly change upon base expansion, while the G:C interaction energy is slightly strengthened upon incorporation of either expanded pyrimidine. The largest effect of base expansion occurs in the stacking energies. Specifically, the maximum (most negative) stacking energies in isolated dimers formed by aligning the nucleobase centers of mass can be increased up to 45% by inclusion of a single y-pyrimidine and up to 55% by consideration of a yy-pyrimidine. Similar increases in the stacking interactions are found when a simplified duplex model composed of two stacked (hydrogen-bonded) base pairs is considered, where both the intrastrand and interstrand stacking interactions can be increased and the effects are more pronounced for the yy-pyrimidines. Moreover, the total stability (sum of all hydrogen-bonding and stacking interactions) is greater for duplexes containing expanded yy-pyrimidines compared to y-pyrimidines, which is mainly due to enhanced stacking interactions. Thus, our calculations suggest that multiple unidirectional increases in the size of the nucleobase spacer can continuously enhance the stability of expanded duplexes.

Can We Accurately Describe the Structure of Adenine Tracts in B-DNA? Reference Quantum-Chemical Computations Reveal Overstabilization of Stacking by Molecular Mechanics

Sequence-dependent local variations of helical parameters, structure, and flexibility are crucial for molecular recognition processes involving B-DNA. A-tracts, i.e., stretches of several consecutive adenines in one strand that are in phase with the DNA helical repeat, mediate significant DNA bending. During the past few decades, there have been intense efforts to understand the sequence dependence of helical parameters in DNA. Molecular dynamics (MD) simulations can provide valuable insights into the molecular mechanism behind the relationship between sequence and structure. However, although recent improvements in empirical force fields have helped to capture many sequence-dependent B-DNA properties, several problems remain, such as underestimation of the helical twist and suspected underestimation of the propeller twist in A-tracts. Here, we employ reference quantum mechanical (QM) calculations, explicit solvent MD, and bioinformatics to analyze the underestimation of propeller twisting of A-tracts in simulations. Although we did not identify a straightforward explanation, we discovered two imbalances in the empirical force fields. The first was overestimation of stacking interactions accompanied by underestimation of base-pairing energy, which we attribute to anisotropic polarizabilities that are not reflected by the isotropic force fields. This may lead to overstacking with potentially important consequences for MD simulations of nucleic acids. The second observed imbalance was steric clash between A(N1) and T(N3) nitrogens of AT base pairs in force-field descriptions, resulting in overestimation of the AT pair stretch in MD simulations. We also substantially extend the available set of benchmark estimated CCSD(T)/CBS data for B-DNA base stacking and provide a code that allows the generation of diverse base-stacking geometries suitable for QM computations with predefined intra- and interbase pair parameters.

Role of wobble base pair geometry for codon degeneracy: purine-type bases at the anticodon wobble position

Codon degeneracy is a key feature of the genetic code, explained by Crick (J Mol Biol 19:548-555, 1966) in terms of imprecision of base pairing at the codon third position (the wobble position) of the codon-anticodon duplex. The Crick wobble rules define, but do not explain, which base pairs are allowed/disallowed at the wobble position of this duplex. This work examines whether the H-bonded configurations of solitary RNA base pairs can in themselves help decide which base pairs are allowed at the wobble position during codon-anticodon pairing. Taking the purine-type bases guanine, hypoxanthine, queuine and adenine as anticodon wobble bases, H-bonded pairing energies and optimized configurations of numerous RNA base pairs are calculated in gas and modeled aqueous phase at the B3LYP/6-31 G(d,p) level. Calculated descriptors of alignment of these solitary base pairs are able to screen between allowed and disallowed base pairs for all cases studied here, except two cases which invoke base-sugar interactions in the codon wobble nucleoside. The exclusion of adenine from the anticodon wobble position cannot be explained on the basis of pairing facility or base pair geometry. These DFT results thus account for the specificity and degeneracy of the genetic code for all cases involving guanine, hypoxanthine and queuine as anticodon wobble bases.

COMPUTATION OF LARGE SYSTEMS WITH AN ECONOMIC BASIS SET: AB INITIO CALCULATIONS OF BIOLOGICAL NUCLEIC ACID BASE PAIRS

We show that an economic basis set, in which the polarization functions are considered only for oxygen and nitrogen atoms of strong electronegativity, can be used to determine reliable structures of nucleic acid base pairs. Mulliken charge analysis and the HOMO-LUMO gap in single-point energy calculations using standard basis sets on the geometric structures optimized with the economic basis set found reasonable agreements with the ones of standard calculations. This study is expected to provide a general guideline for basis set selection in the computation of large biological systems being performed with considerable high accuracy, using a low cost computation resource.

Counterpoise correction is not useful for short and Van der Waals distances but may be useful at long range

This article investigates the errors in supermolecule calculations for the helium dimer. In a full CI calculation, there are two errors. One is the basis set superposition error (BSSE), the other is the basis set convergence error (BSCE). Both of the errors arise from the incompleteness of the basis set. These two errors make opposite contributions to the interaction energies. The BSCE is by far the largest error in the short range and larger than (but much closer to) BSSE around the Van der Waals minimum. Only at the long range, the BSSE becomes the larger error. The BSCE and BSSE largely cancel each other over the Van der Waals well. Accordingly, it may be recommended to not include the BSSE for the calculation of the potential energy curve from short distance till well beyond the Van der Waals minimum, but it may be recommended to include the BSSE correction if an accurate tail behavior is required. Only if the calculation has used a very large basis set, one can refrain from including the counterpoise correction in the full potential range. These results are based on full CI calculations with the aug-cc-pVXZ (X = D, T, Q, 5) basis sets.

Effects of the biological backbone on stacking interactions at DNA-protein interfaces: the interplay between the backbone···π and π···π components

The (gas-phase) MP2/6-31G*(0.25) π···π stacking interactions between the five natural bases and the aromatic amino acids calculated using (truncated) monomers composed of conjugated rings and/or (extended) monomers containing the biological backbone (either the protein backbone or deoxyribose sugar) were previously compared. Although preliminary energetic results indicated that the protein backbone strengthens, while the deoxyribose sugar either strengthens or weakens, the interaction calculated using truncated models, the reasons for these effects were unknown. The present work explains these observations by dissecting the interaction energy of the extended complexes into individual backbone···π and π···π components. Our calculations reveal that the total interaction energy of the extended complex can be predicted as a sum of the backbone···π and π···π components, which indicates that the biological backbone does not significantly affect the ring system through π-polarization. Instead, we find that the backbone can indirectly affect the magnitude of the π···π contribution by changing the relative ring orientations in extended dimers compared with truncated dimers. Furthermore, the strengths of the individual backbone···π contributions are determined to be significant (up to 18 kJ mol(-1)). Therefore, the origin of the energetic change upon model extension is found to result from a balance between an additional (attractive) backbone···π component and differences in the strength of the π···π interaction. In addition, to understand the effects of the biological backbone on the stacking interactions at DNA-protein interfaces in nature, we analyzed the stacking interactions found in select DNA-protein crystal structures, and verified that an additive approach can be used to examine the strength of these interactions in biological complexes. Interestingly, although the presence of attractive backbone···π contacts is qualitatively confirmed using the quantum theory of atoms in molecules (QTAIM), QTAIM electron density analysis is unable to quantitatively predict the additive relationship of these interactions. Most importantly, this work reveals that both the backbone···π and π···π components must be carefully considered to accurately determine the overall stability of DNA-protein assemblies.

The assessment of density functionals for DNA–protein stacked and T-shaped complexes

The present work uses 129 nucleobase – amino acid CCSD(T)/CBS stacking and T-shaped interaction energies as reference data to test the ability of various density functionals with double-zeta quality basis sets, as well as some semi-empirical and molecular mechanics methods, to accurately describe noncovalent DNA–protein π–π and π+–π interactions. The goal of this work is to identify methods that can be used in hybrid approaches (QM/MM, ONIOM) for large-scale modeling of enzymatic systems involving active-site (substrate) π–π contacts. Our results indicate that AMBER is a more appropriate choice for the lower-level method in hybrid techniques than popular semi-empirical methods (AM1, PM3), and suggest that AMBER accurately describes the π–π interactions found throughout DNA–protein complexes. The M06–2X and PBE-D density functionals were found to provide very promising descriptions of the 129 nucleobase – amino acid interaction energies, which suggests that these may be the most suitable methods for describing high-level regions. Therefore, M06–2X and PBE-D with both the 6–31G(d) and 6–31+G(d,p) basis sets were further examined through potential-energy surface scans to better understand how these techniques describe DNA–protein π–π interactions in both minimum and nonminimum regions of the potential-energy surfaces, which is critical information when modeling enzymatic reaction pathways. Our results suggest that studies of stacked nucleobase – amino acid systems should implement the PBE-D/6–31+G(d,p) method. However, if T-shaped contacts are involved and (or) smaller basis sets must be considered due to limitations in computational resources, then M06–2X/6–31G(d) provides an overall excellent description of both nucleobase – amino acid stacking and T-shaped interactions for a range of DNA–protein π–π and π+–π interactions.

Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations?

High level ab initio methods have been used to study stacking interactions in ten unique base pair steps both in A-RNA and in B-DNA duplexes. The protocol for selection of geometries based on molecular dynamics (MD) simulations is proposed, and its suitability is demonstrated by comparison with stacking in steps at fiber diffraction geometries. It is shown that fiber diffraction geometries are not sufficiently accurate for interaction energy calculations. In addition, the protocol for selection of geometries based on MD simulations allows for the evaluation of the variability of the intrinsic stacking energies along the MD trajectories. The uncertainty in stacking energies (difference between the most and least stable geometry) due to the dynamical nature of systems can be, in some cases, as large as 3.0 kcal x mol(-1), which is almost 50% of the actual sequence dependence of base stacking energies (the energy difference between the most and least stable sequences). Thus, assessing the relative magnitude of the gas phase stacking energy using a single geometry for each sequence is insufficient to obtain an unambiguous order of gas phase stacking energies in canonical double helices. Though the ordering of ten unique dinucleotide steps cannot be definitive, some general conclusions were drawn. The stacking energies of base pair steps in A-RNA are more evenly separated compared to B-DNA, and their ordering is less sensitive to the dynamics of the system compared to be B-DNA. The most stable step both in B-DNA and A-RNA is the GC/GC [corrected] step that is well separated from the second most stable step CG/CG. [corrected] Also the least stable step (the CC/GG step) is well separated from the rest of the structures. The calculations further show that B-DNA stacking is favorable only marginally (on average by 1.14 kcal x mol(-1) per base pair step) over A-RNA stacking, and this difference vanishes after subtracting the stabilizing van der Waals effect of the thymine 5-methyl group that is absent in RNA. Basically, no correlation between the sequence dependence of gas phase stacking energies and the sequence dependence of DeltaG degrees(37) free energies used in nearest-neighbor models was found either for B-DNA or for A-RNA. This reflects the complexity of the balance of forces that are responsible for the sequence dependence of thermodynamics stability of nucleic acids, which masks the effect of the intrinsic interactions between the stacked base pairs.

Assessment of standard force field models against high-quality ab initio potential curves for prototypes of pi-pi, CH/pi, and SH/pi interactions

Several popular force fields, namely, CHARMM, AMBER, OPLS-AA, and MM3, have been tested for their ability to reproduce highly accurate quantum mechanical potential energy curves for noncovalent interactions in the benzene dimer, the benzene-CH(4) complex, and the benzene-H(2)S complex. All of the force fields are semi-quantitatively correct, but none of them is consistently reliable quantitatively. Re-optimization of Lennard-Jones parameters and symmetry-adapted perturbation theory analysis for the benzene dimer suggests that better agreement cannot be expected unless more flexible functional forms (particularly for the electrostatic contributions) are employed for the empirical force fields.

Effects of the biological backbone on DNA-protein stacking interactions

The pi-pi stacking (face-to-face) interactions between the five natural DNA or RNA nucleobases and the four aromatic amino acids were compared using three different types of dimers: (1) a truncated nucleoside (nucleobase) stacked with a truncated amino acid; (2) a truncated nucleoside (nucleobase) stacked with an extended amino acid; and (3) a nucleoside (extended nucleobase) stacked with a truncated amino acid. Systematic (MP2/6-31G*(0.25)) potential energy surface scans reveal important information about the effects of the deoxyribose sugar and protein backbone on the structure and binding energy between truncated nucleobase and amino acid models that are typically implemented in the literature. Most notably, electrostatic and steric interactions arising from the bulkiness of the biological backbones can change the preferred relative orientations of DNA and protein pi-systems. More importantly, the protein backbone can strengthen the stacking energy (by up to 10 kJ mol(-1)), while the deoxyribose moiety can strengthen or weaken the stacking interaction depending on the positioning of the amino acid relative to the sugar residue. These effects are likely due to additional interactions between the amino acid or nucleobase ring and the backbone in the extended monomer rather than significant changes in the properties of the biological pi-systems upon model extension. Since the present work reveals that all calculated DNA-protein stacking interactions are significant and approach the strength of other noncovalent interactions between biomolecules, both pi-pi and backbone-pi interactions must be considered when attempting to gain a complete picture of DNA-protein binding.

Hydrogen bonding and pi-stacking: how reliable are force fields? A critical evaluation of force field descriptions of nonbonded interactions

We have evaluated the performance of a set of widely used force fields by calculating the geometries and stabilization energies for a large collection of intermolecular complexes. These complexes are representative of a range of chemical and biological systems for which hydrogen bonding, electrostatic, and van der Waals interactions play important roles. Benchmark energies are taken from the high-level ab initio values in the JSCH-2005 and S22 data sets. All of the force fields underestimate stabilization resulting from hydrogen bonding, but the energetics of electrostatic and van der Waals interactions are described more accurately. OPLSAA gave a mean unsigned error of 2 kcal mol(-1) for all 165 complexes studied, and outperforms DFT calculations employing very large basis sets for the S22 complexes. The magnitude of hydrogen bonding interactions are severely underestimated by all of the force fields tested, which contributes significantly to the overall mean error; if complexes which are predominantly bound by hydrogen bonding interactions are discounted, the mean unsigned error of OPLSAA is reduced to 1 kcal mol(-1). For added clarity, web-based interactive displays of the results have been developed which allow comparisons of force field and ab initio geometries to be performed and the structures viewed and rotated in three dimensions.

Theoretical insights into the nature of intermolecular interactions in cytosine dimer

In this study we discuss stacking interactions in cytosine dimer in conformations appearing in B-DNA crystals. The variational-perturbational scheme was applied for decomposition of the intermolecular interaction energy at the MP2 level of theory. The significant influence of the mutual orientation of cytosine monomers was observed not only on the total intermolecular interaction energy but also on its components: Different components of intermolecular interaction energy depend in different manner on parameters describing mutual orientation of cytosine monomers.

Probing phenylalanine/adenine pi-stacking interactions in protein complexes with explicitly correlated and CCSD(T) computations

To examine the effects of pi-stacking interactions between aromatic amino acid side chains and adenine bearing ligands in crystalline protein structures, 26 toluene/(N9-methyl)adenine model configurations have been constructed from protein/ligand crystal structures. Full geometry optimizations with the MP2 method cause the 26 crystal structures to collapse to six unique structures. The complete basis set (CBS) limit of the CCSD(T) interaction energies has been determined for all 32 structures by combining explicitly correlated MP2-R12 computations with a correction for higher-order correlation effects from CCSD(T) calculations. The CCSD(T) CBS limit interaction energies of the 26 crystal structures range from -3.19 to -6.77 kcal mol (-1) and average -5.01 kcal mol (-1). The CCSD(T) CBS limit interaction energies of the optimized complexes increase by roughly 1.5 kcal mol (-1) on average to -6.54 kcal mol (-1) (ranging from -5.93 to -7.05 kcal mol (-1)). Corrections for higher-order correlation effects are extremely important for both sets of structures and are responsible for the modest increase in the interaction energy after optimization. The MP2 method overbinds the crystal structures by 2.31 kcal mol (-1) on average compared to 4.50 kcal mol (-1) for the optimized structures.

Geometrical and electronic structure variability of the sugar-phosphate backbone in nucleic acids

The anionic sugar-phosphate backbone of nucleic acids substantially contributes to their structural flexibility. To model nucleic acid structure and dynamics correctly, the potentially sampled substates of the sugar-phosphate backbone must be properly described. However, because of the complexity of the electronic distribution in the nucleic acid backbone, its representation by classical force fields is very challenging. In this work, the three-dimensional potential energy surfaces with two independent variables corresponding to rotations around the alpha and gamma backbone torsions are studied by means of high-level ab initio methods (B3LYP/6-31+G*, MP2/6-31+G*, and MP2 complete basis set limit levels). The ability of the AMBER ff99 [Wang, J. M.; Cieplak, P.; Kollman, P. A. J. Comput. Chem. 2000, 21, 1049-1074] and parmbsc0 [Perez, A.; Marchan, I.; Svozil, D.; Sponer, J.; Cheatham, T. E.; Laughten, C. A.; Orozco, M. Biophys. J. 2007, 92, 3817-3829] force fields to describe the various alpha/gamma conformations of the DNA backbone accurately is assessed by comparing the results with those of ab initio quantum chemical calculations. Two model systems differing in structural complexity were used to describe the alpha/gamma energetics. The simpler one, SPM, consisting of a sugar and methyl group linked through a phosphodiester bond was used to determine higher-order correlation effects covered by the CCSD(T) method. The second, more complex model system, SPSOM, includes two deoxyribose residues (without the bases) connected via a phosphodiester bond. It has been shown by means of a natural bond orbital analysis that the SPSOM model provides a more realistic representation of the hyperconjugation network along the C5'-O5'-P-O3'-C3' linkage. However, we have also shown that quantum mechanical investigations of this model system are nontrivial because of the complexity of the SPSOM conformational space. A comparison of the ab initio data with the ff99 potential energy surface clearly reveals an incorrect ff99 force-field description in the regions where the gamma torsion is in the trans conformation. An explanation is proposed for why the alpha/gamma flips are eliminated so successfully when the parmbsc0 force-field modification is used.