Principles of RNA base pairing: structures and energies of the trans Watson-Crick/sugar edge base pairs

Due to the presence of the 2'-OH hydroxyl group of ribose, RNA molecules utilize an astonishing variability of base pairing patterns to build up their structures and perform the biological functions. Many of the key RNA base pairing families have no counterparts in DNA. In this study, the trans Watson-Crick/sugar edge (trans WC/SE) RNA base pair family has been characterized using quantum chemical and molecular mechanics calculations. Gas-phase optimized geometries from density functional theory (DFT) calculations and RIMP2 interaction energies are reported for the 10 crystallographically identified trans WC/SE base pairing patterns. Further, stable structures are predicted for all of the remaining six possible members of this family not seen in RNAs so far. Among these novel six base pairs, the computations substantially refine two structures suggested earlier based on simple isosteric considerations. For two additional trans WC/SE base pairs predicted in this study, no arrangement was suggested before. Thus, our study brings a complete set of trans WC/SE base pairing patterns. The present results are also contrasted with calculations reported recently for the cis WC/SE base pair family. The computed base pair sizes are in sound correlation with the X-ray data for all WC/SE pairing patterns including both their cis and trans isomers. This confirms that the isostericity of RNA base pairs, which is one of the key factors determining the RNA sequence conservation patterns, originates in the properties of the isolated base pairs. In contrast to the cis structures, however, the isosteric subgroups of the trans WC/SE family differ not only in their H-bonding patterns and steric dimensions but also in the intrinsic strength of the intermolecular interactions. The distribution of the total interaction energy over the sugar-base and base-base contributions is controlled by the cis-trans isomerism.     

Non-Watson-Crick base pairing in RNA. quantum chemical analysis of the cis Watson-Crick/sugar edge base pair family

Large RNA molecules exhibit an astonishing variability of base-pairing patterns, while many of the RNA base-pairing families have no counterparts in DNA. The cis Watson-Crick/sugar edge (cis WC/SE) RNA base pairing is investigated by ab initio quantum chemical calculations. A detailed structural and energetic characterization of all 13 crystallographically detected members of this family is provided by means of B3LYP/6-31G and RIMP2/aug-cc-pVDZ calculations. Further, a prediction is made for the remaining 3 cis WC/SE base pairs which are yet to be seen in the experiments. The interaction energy calculations point at the key role of the 2'-OH group in stabilizing the sugar-base contact and predict all 16 cis WC/SE base-pairing patterns to be nearly isoenergetic. The perfect correlation of the main geometrical parameters in the gas-phase optimized and X-ray structures shows that the principle of isosteric substitutions in RNA is rooted from the intrinsic structural similarity of the isolated base pairs. The present quantum chemical calculations for the first time analyze base pairs involving the ribose 2'-OH group and unambiguously correlate the structural information known from experiments with the energetics of interactions. The calculations further show that the relative importance and absolute value of the dispersion energy in the cis WC/SE base pairs are enhanced compared to the standard base pairs. This may by an important factor contributing to the strength of such interactions when RNA folds in its polar environment. The calculations further demonstrate that the Cornell et al. force field commonly used in molecular modeling and simulations provides satisfactory performance for this type of RNA interactions.     

A benchmark for non-covalent interactions in solids

A benchmark for non-covalent interactions in solids (C21) based on the experimental sublimation enthalpies and geometries of 21 molecular crystals is presented. Thermal and zero-point effects are carefully accounted for and reference lattice energies and thermal pressures are provided, which allow dispersion-corrected density functionals to be assessed in a straightforward way. Other thermal corrections to the sublimation enthalpy (the 2RT term) are reexamined. We compare the recently implemented exchange-hole dipole moment (XDM) model with other approaches in the literature to find that XDM roughly doubles the accuracy of DFT-D2 and non-local functionals in computed lattice energies (4.8 kJ/mol mean absolute error) while, at the same time, predicting cell geometries within less than 2% of the experimental result on average. The XDM model of dispersion interactions is confirmed as a very promising approach in solid-state applications.     

The interaction of carbohydrates and amino acids with aromatic systems studied by density functional and semi-empirical molecular orbital calculations with dispersion corrections

Density functional theory (DFT-D) and semi-empirical (PM3-D) methods having an added dispersion correction have been used to study stabilising carbohydrate-aromatic and amino acid-aromatic interactions. The interaction energy for three simple sugars in different conformations with benzene, all give interaction energies close to 5 kcal mol(-1). Our original parameterization of PM3 (PM3-D) seriously overestimates this value, and has prompted a reparametrization which includes a modified core-core interaction term. With two additional parameters, the carbohydrate complexes, as well as the S22 data set, are well reproduced. The new PM3 scheme (PM3-D*) is found to describe the peptide bond-aromatic ring interactions accurately and, together with the DFT-D method, it is used to investigate the interaction of six amino acids with pyrene. Whilst the peptide backbone can adopt both stacked and T-shaped structures in the complexes with similar interaction energies, there is a preference for the unsaturated ring to adopt a stacked structure. Thus, peptides in which the latter interactions are maximised are likely to be the most effective for the functionalisation of carbon nanotubes.     

Van der Waals interactions in aromatic systems: structure and energetics of dimers and trimers of pyridine.

Full geometry optimizations at the dispersion-corrected DFT-BLYP level of theory were carried out for dimers and trimers of pyridine. The DFT-D interaction energies were checked against results from single-point SCS-MP2/aug-cc-pVTZ calculations. Three stacked structures and a planar H-bonded dimer were found to be very close in energy (interaction energies in the range from -3.4 to -4.0 kcal mol(-1)). Two T-shaped geometries are higher lying, by about 1 kcal mol(-1), which is explained by the more favorable electrostatic interactions in the stacked and H-bonded arrangements. The DFT-D approach has proved to be a reliable and efficient tool to explore the conformational space of aromatic van der Waals complexes and furthermore provides interaction energies with errors of less than 10-20 % of DeltaE. Comparisons with previous results obtained by using only partially optimized model geometries strongly indicate that unconstrained optimizations are mandatory in such weakly bonded low-symmetry systems.     

Intermolecular interaction energies in molecular crystals: comparison and agreement of localized Møller-Plesset 2, dispersion-corrected density functional, and classical empirical two-body calculations

A comparative analysis of the intermolecular energy for a data set including 60 molecular crystals with a large variety of functional groups has been carried out using three different computational approaches: (i) a method based on a physically meaningful empirical partition of the interaction energy (PIXEL), (ii) density functional methods with a posteriori empirical correction for the dispersion interactions (DFT-D), and (iii) a full periodic ab initio quantum mechanical method based on M?ller-Plesset perturbation theory for the electron correlation using localized crystal orbitals (LMP2). Due to the large computational cost, LMP2 calculations have been restricted to a subset of seven molecular crystal comprising benzene, formic acid, formamide, succinic anhydride, urea, oxalic acid, and nitroguanidine, and the results compared with PIXEL and DFT-D data as well as with the experimental data show excellent agreement among all adopted methods. This shows that both DFT-D and PIXEL approaches are robust predictive tools for studying molecular crystals. A detailed analysis shows a very similar dispersion contribution of the two methods across the 60 considered molecular crystals. The study also confirms that pure DFT shows serious deficiencies in properly handling molecular crystals in which the dispersive contribution is large. Due to the negligible requested computational resources, PIXEL is the method of choice in screening of a large number of molecular crystals, an essential step to predict crystal polymorphism or to study crystal growth processes. DFT-D can then be used to refine the ranking emerged from PIXEL calculations due to its general applicability and robustness in properly handling short-range interactions.     

Effect of the damping function in dispersion corrected density functional theory

It is shown by an extensive benchmark on molecular energy data that the mathematical form of the damping function in DFT-D methods has only a minor impact on the quality of the results. For 12 different functionals, a standard "zero-damping" formula and rational damping to finite values for small interatomic distances according to Becke and Johnson (BJ-damping) has been tested. The same (DFT-D3) scheme for the computation of the dispersion coefficients is used. The BJ-damping requires one fit parameter more for each functional (three instead of two) but has the advantage of avoiding repulsive interatomic forces at shorter distances. With BJ-damping better results for nonbonded distances and more clear effects of intramolecular dispersion in four representative molecular structures are found. For the noncovalently-bonded structures in the S22 set, both schemes lead to very similar intermolecular distances. For noncovalent interaction energies BJ-damping performs slightly better but both variants can be recommended in general. The exception to this is Hartree-Fock that can be recommended only in the BJ-variant and which is then close to the accuracy of corrected GGAs for non-covalent interactions. According to the thermodynamic benchmarks BJ-damping is more accurate especially for medium-range electron correlation problems and only small and practically insignificant double-counting effects are observed. It seems to provide a physically correct short-range behavior of correlation/dispersion even with unmodified standard functionals. In any case, the differences between the two methods are much smaller than the overall dispersion effect and often also smaller than the influence of the underlying density functional.     

π-Hole Interactions Involving Nitro Aromatic Ligands in Protein Structures

Studying noncanonical intermolecular interactions between a ligand and a protein constitutes an emerging research field. Identifying synthetically accessible molecular fragments that can engage in intermolecular interactions is a key objective in this area. Here, it is shown that so-called “π-hole interactions” are present between the nitro moiety in nitro aromatic ligands and lone pairs within protein structures (water and protein carbonyls and sulfurs). Ample structural evidence was found in a PDB analysis and computations reveal interaction energies of about −5 kcal mol⁻¹ for ligand–protein π-hole interactions. Several examples are highlighted for which a π-hole interaction is implicated in the superior binding affinity or inhibition of a nitro aromatic ligand versus a similar non-nitro analogue. The discovery that π-hole interactions with nitro aromatics are significant within protein structures parallels the finding that halogen bonds are biologically relevant. This has implications for the interpretation of ligand–protein complexation phenomena, for example, involving the more than 50 approved drugs that contain a nitro aromatic moiety.     

Potential-energy and free-energy surfaces of glycyl-phenylalanyl-alanine (GFA) tripeptide: experiment and theory

The free-energy surface (FES) of glycyl-phenylalanyl-alanine (GFA) tripeptide was explored by molecular dynamics (MD) simulations in combination with high-level correlated ab initio quantum chemical calculations and metadynamics. Both the MD and metadynamics employed the tight-binding DFT-D method instead of the AMBER force field, which yielded inaccurate results. We classified the minima localised in the FESs as follows: a) the backbone-conformational arrangement; and b) the existence of a COOH...OC intramolecular H-bond (families CO(2)H(free) and CO(2)H(bonded)). Comparison with experimental results showed that the most stable minima in the FES correspond to the experimentally observed structures. Remarkably, however, we did not observe experimentally the CO(2)H(bonded) family (also predicted by metadynamics), although its stability is comparable to that of the CO(2)H(free) structures. This fact was explained by the former's short excited-state lifetime. We also carried out ab initio calculations using DFT-D and the M06-2X functional. The importance of the dispersion energy in stabilising peptide conformers is well reflected by our pioneer analysis using the DFT-SAPT method to explore the nature of the backbone/side-chain interactions.     

Theoretical analysis of antisense duplexes: determinants of the RNase H susceptibility

The structure and dynamic properties of different antisense related duplexes (DNA x RNA, 2'O-Me-DNA x RNA, 2'F-ANA x RNA, C5(Y)-propynyl-DNA x RNA, ANA x RNA, and control duplexes DNA x DNA and RNA x RNA) have been determined by means of long molecular dynamics simulations (covering more than 0.5 micros of fully solvated unrestrained MD simulation). The massive analysis presented here allows us to determine the subtle differences between the different duplexes, which in all cases pertain to the same structural family. This analysis provides information on the molecular determinants that allow RNase H to recognize and degrade some of these duplexes, whereas others with apparently similar conformations are not affected. Subtle structural and deformability features define the key properties used by RNase H to discriminate between duplexes.     

Effects of hydrophobic and dipole-dipole interactions on the conformational transitions of a model polypeptide

We studied the effects of hydrophobicity and dipole-dipole interactions between the nearest-neighbor amide planes on the secondary structures of a model polypeptide by calculating the free energy differences between different peptide structures. The free energy calculations were performed with low computational costs using the accelerated Monte Carlo simulation (umbrella sampling) method, with a bias-potential method used earlier in our accelerated molecular dynamics simulations. It was found that the hydrophobic interaction enhances the stability of alpha helices at both low and high temperatures but stabilizes beta structures only at high temperatures at which alpha helices are not stable. The nearest-neighbor dipole-dipole interaction stabilizes beta structures under all conditions, especially in the low temperature region where alpha helices are the stable structures. Our results indicate clearly that the dipole-dipole interaction between the nearest neighboring amide planes plays an important role in determining the peptide structures. Current research provides a more unified and quantitative picture for understanding the effects of different forms of interactions on polypeptide structures. In addition, the present model can be extended to describe DNA/RNA, polymer, copolymer, and other chain systems.     

Stacking effects on local structure in RNA: changes in the structure of tandem GA pairs when flanking GC pairs are replaced by isoG-isoC pairs

The Watson-Crick-like isoG-isoC (iGiC) pair, with the amino and carbonyl groups transposed relative to the Watson-Crick GC pair, provides an expanded alphabet for understanding interactions that shape nucleic acid structure. Here, thermodynamic stabilities of tandem GA pairs flanked by iGiC pairs are reported along with the NMR structures of the RNA self-complementary duplexes (GCiGGAiCGCA)2 and (GGiCGAiGCCA)2. A sheared GA pairing forms in (GCiGGAiCGCA)2, and an imino GA pairing forms in (GGiCGAiGCCA)2. The structures contrast with the formation of tandem imino and sheared GA pairs flanked by GC pairs in the RNA self-complementary duplexes (GCGGACGC)2 and (GGCGAGCC)2, respectively. In both iGiC duplexes, Watson-Crick-like hydrogen bonds are formed between iG and iC, and iGiC substitutions result in less favorable loop stability. The results provide benchmarks for testing computations of molecular interactions that shape RNA three-dimensional structure.     

Anisotropic NMR interaction tensors in the decamethylaluminocenium cation

Solid-state NMR experiments, analytical and numerical simulations of solid NMR powder patterns, ab initio self-consistent field and hybrid density functional theory calculations, and single-crystal X-ray diffraction are used to characterize the molecular structure and anisotropic NMR interaction tensors in the bis(pentamethylcyclopentadienyl)aluminum cation, [Cp(2)Al](+). This highly symmetric main group metallocene has a structure analogous to that of ferrocene and the cobaltocenium cation. The single-crystal X-ray diffraction structure is reported for [Cp(2)Al][AlCl(4)]. Solid-state (27)Al[(1)H] magic-angle spinning and static NMR experiments are used to study the aluminum chemical shielding and electric field gradient tensors, revealing axial symmetry in both cases with a large chemical shielding span of Omega = 83(3) ppm and a small nuclear quadrupole coupling constant, C(Q)((27)Al) = 0.86(10) MHz. Carbon-13 CPMAS NMR experiments in combination with ab initio calculations and simulations of the effects of chemical exchange on (13)C static powder patterns reveal dynamic rotation of rings and suggest a low internal rotational barrier for this process. Theoretical computations of interaction tensors using the Gaussian 98 and Amsterdam Density Functional software packages are in good agreement with experiment and lend insight into the molecular origin of these NMR interactions. Orientations of the NMR tensors determined from theory, the large chemical shielding span, and the very small value of C(Q)((27)Al) can all be rationalized in terms of the high molecular symmetry.     

Prediction of protein–protein complexes using replica exchange with repulsive scaling

The realistic prediction of protein–protein complex structures is import to ultimately model the interaction of all proteins in a cell and for the design of new protein–protein interactions. In principle, molecular dynamics (MD) simulations allow one to follow the association process under realistic conditions including full partner flexibility and surrounding solvent. However, due to the many local binding energy minima at the surface of protein partners, MD simulations are frequently trapped for long times in transient association states. We have designed a replica-exchange based scheme employing different levels of a repulsive biasing between partners in each replica simulation. The bias acts only on intermolecular interactions based on an increase in effective pairwise van der Waals radii (repulsive scaling (RS)-REMD) without affecting interactions within each protein or with the solvent. For a set of five protein test cases (out of six) the RS-REMD technique allowed the sampling of near-native complex structures even when starting from the opposide site with respect to the native binding site for one partner. Using the same start structures and same computational demand regular MD simulations sampled near native complex structures only for one case. The method showed also improved results for the refinement of docked structures in the vicinity of the native binding geometry compared to regular MD refinement.     

Dispersion interactions of carbohydrates with condensate aromatic moieties: theoretical study on the CH-π interaction additive properties

In this article we present the first systematic study of the additive properties (i.e. degree of additivity) of the carbohydrate-aromatic moiety CH-π dispersion interaction. The additive properties were studied on the β-D-glucopyranose, β-D-mannopyranose and α-L-fucopyranose complexes with the naphthalene molecule by comparing the monodentate (single CH-π) and bidentate (two CH-π) complexes. All model complexes were optimized using the DFT-D approach, at the BP/def2-TZVPP level of theory. The interaction energies were refined using single point calculations at highly correlated ab initio methods at the CCSD(T)/CBS level, calculated as E + (E(CCSD(T))-E(MP2))(Small Basis). Bidentate complexes show very strong interactions in the range from -10.79 up to -7.15 and -8.20 up to -6.14 kcal mol(-1) for the DFT-D and CCSD(T)/CBS level, respectively. These values were compared with the sum of interaction energies of the appropriate monodentate carbohydrate-naphthalene complexes. The comparison reveals that the bidentate complex interaction energy is higher (interaction is weaker) than the sum of monodentate complex interaction energies. Bidentate complex interaction energy corresponds to 2/3 of the sum of the appropriate monodentate complex interaction energies (averaging over all modeled carbohydrate complexes). The observed interaction energies were also compared with the sum of interaction energies of the corresponding previously published carbohydrate-benzene complexes. Also in this case the interaction energy of the bidentate complex was higher (i.e. weaker interaction) than the sum of interaction energies of the corresponding benzene complexes. However, the obtained difference is lower than before, while the bidentate complex interaction energy corresponds to 4/5 of the sum of interaction energy of the benzene complexes, averaged over all structures. The mentioned comparison might aid protein engineering efforts where amino acid residues phenylalanine or tyrosine are to be replaced by a tryptophan and can help to predict the changes in the interactions. The observed results also show that DFT-D correctly describes the CH-π interaction energy and their additive properties in comparison to CCSD(T)/CBS calculated interaction energies. Thus, the DFT-D approach might be used for calculation of larger complexes of biological interest, where dispersion interaction plays an important role.     

Structural reorganizations control intermolecular conductance and charge trapping in paraquat-tetraphenylborate inverse photochemical cell

Miniaturization of electronic devices to the level of single molecules requires detailed understanding of the mechanisms of their operation. One of the questions here is the identification of the role of structural alterations in charge separation and stabilization in photoactive complexes. To address this question, we calculate optimized molecular and electronic structures, and optical and vibrational spectra of l,l'-dimethyl 4,4'-bipyridinium-bis tetraphenylborate PQ(BPh(4))(2) complex ab initio using density functional theory approach and compare them with the experimentally observed UV-Vis and Raman spectra of the molecules in solid-state films. The results indicate that the association of PQ and BPh(4) leads to the formation of an internally ionized structure that is accompanied by the structural reorganization of both PQ (the twisting of pyridinium rings) and BPh(4) (phenyl rings rotation) moieties. The quanta of light do not seem to be directly involved in the formation of this ionized structure, but provide energy for fast recombination of the separated charges between BPh(4)(-) and PQ(2+). The high efficiency of the dark charge separation and the stabilization of separated charges in the complex permit the using of PQ(BPh(4))(2) in various charge-transfer devices like molecular probes, photovoltaic devices or chemical memory units.     

Accurate description of van der Waals complexes by density functional theory including empirical corrections

An empirical method to account for van der Waals interactions in practical calculations with the density functional theory (termed DFT-D) is tested for a wide variety of molecular complexes. As in previous schemes, the dispersive energy is described by damped interatomic potentials of the form C6R(-6). The use of pure, gradient-corrected density functionals (BLYP and PBE), together with the resolution-of-the-identity (RI) approximation for the Coulomb operator, allows very efficient computations for large systems. Opposed to previous work, extended AO basis sets of polarized TZV or QZV quality are employed, which reduces the basis set superposition error to a negligible extend. By using a global scaling factor for the atomic C6 coefficients, the functional dependence of the results could be strongly reduced. The "double counting" of correlation effects for strongly bound complexes is found to be insignificant if steep damping functions are employed. The method is applied to a total of 29 complexes of atoms and small molecules (Ne, CH4, NH3, H2O, CH3F, N2, F2, formic acid, ethene, and ethine) with each other and with benzene, to benzene, naphthalene, pyrene, and coronene dimers, the naphthalene trimer, coronene. H2O and four H-bonded and stacked DNA base pairs (AT and GC). In almost all cases, very good agreement with reliable theoretical or experimental results for binding energies and intermolecular distances is obtained. For stacked aromatic systems and the important base pairs, the DFT-D-BLYP model seems to be even superior to standard MP2 treatments that systematically overbind. The good results obtained suggest the approach as a practical tool to describe the properties of many important van der Waals systems in chemistry. Furthermore, the DFT-D data may either be used to calibrate much simpler (e.g., force-field) potentials or the optimized structures can be used as input for more accurate ab initio calculations of the interaction energies.     

Leading RNA tertiary interactions: structures, energies, and water insertion of A-minor and P-interactions. A quantum chemical view

Complex molecular shapes of ribosomal RNA molecules are stabilized by recurrent types of tertiary interactions involving highly specific and conserved non-Watson-Crick base pairs, triplets, and quartets. We analyzed the intrinsic structure and stability of the P-motif and the four basic types of A-minor interactions (types I, II, III, and 0), which represent the most prominent RNA tertiary interaction patterns refined in the course of evolution. In the studied interactions, the electron correlation component of the stabilization usually exceeds the Hartree-Fock (HF) term, leading to a strikingly different balance of forces as compared to standard base pairing stabilized primarily by the HF term. In other words, the A-minor and P-interactions are considerably more influenced by the dispersion energy as compared to canonical base pairs, which makes them particularly suitable to zip the folded RNA structures that are substantially hydrated even in their interior. Continuum solvent COSMO calculations confirm that the stability of the canonical GC base pair is affected (reduced) by the continuum solvent screening considerably more than the stability of the A-minor interaction. Among the studied systems, the strong A-minor II and weak A-minor III interactions require water molecules to stabilize the experimental geometry. Gas-phase optimization of the canonical A-minor II A/CG triplet without water results in a geometry that is clearly inconsistent with the RNA structure. The gas-phase structure of the P-interaction and the most stable A-minor I interaction nicely agrees with the geometries occurring in the ribosome. A-minor I can also adopt an alternative water-mediated substate rather often observed in X-ray and molecular dynamics studies. The A-minor I water bridge, however, does not appear to stabilize the tertiary contact, and its role is to provide structural flexibility to this binding pattern within the context of the RNA structure. Interestingly, the insertion of a polar water molecule in the A-minor I A/CG tertiary contact occurring in the A/C tertiary pair is stabilized primarily by the HF (electrostatic) interaction energy, while the dispersion-controlled A/G contact remains firmly bound. Thus, the intrinsic balance of forces as revealed by quantum mechanics (QM) calculations nicely correlates with many behavioral aspects of the studied interactions inside RNA. The comparison of interaction energies computed using quantum chemistry and an AMBER force field reveals that common molecular mechanics calculations perform rather well, except that the strength of the P-interaction is modestly overestimated. We also briefly discuss the non-negligible methodological differences when evaluating simple base-base nucleic acids base pairs and the complex RNA tertiary base pairing patterns using QM procedures.