Electrofreezing of Liquid Ammonia

Here we prove that, in addition to temperature and pressure, another important thermodynamic variable permits the exploration of the phase diagram of ammonia: the electric field. By means of (path integral) ab initio molecular dynamics simulations, we predict that, upon applying intense electric fields on ammonia, the electrofreezing phenomenon occurs, leading the liquid toward a novel ferroelectric solid phase. This study proves that electric fields can generally be exploited as the access key to otherwise-unreachable regions in phase diagrams, unveiling the existence of new condensed-phase structures. Furthermore, the reported findings have manifold practical implications, from the safe storage and transportation of ammonia to the understanding of the solid structures this compound forms in planetary contexts.

A Computational Quantum-Based Perspective on the Molecular Origins of Life’s Building Blocks

The search for the chemical origins of life represents a long-standing and continuously debated enigma. Despite its exceptional complexity, in the last decades the field has experienced a revival, also owing to the exponential growth of the computing power allowing for efficiently simulating the behavior of matter—including its quantum nature—under disparate conditions found, e.g., on the primordial Earth and on Earth-like planetary systems (i.e., exoplanets). In this minireview, we focus on some advanced computational methods capable of efficiently solving the Schro¨dinger equation at different levels of approximation (i.e., density functional theory)—such as ab initio molecular dynamics—and which are capable to realistically simulate the behavior of matter under the action of energy sources available in prebiotic contexts. In addition, recently developed metadynamics methods coupled with first-principles simulations are here reviewed and exploited to answer to old enigmas and to propose novel scenarios in the exponentially growing research field embedding the study of the chemical origins of life.

Synthesis of (d)-erythrose from glycolaldehyde aqueous solutions under electric field

A Miller-like numerical experiment demonstrates that ubiquitarious molecules such as water and glycolaldehyde can synthesize (d)-erythrose, one of the direct precursors of ribose.

Channeling through Two Stacked Guanine Quartets of One and Two Alkali Cations in the Li+, Na+, K+, and Rb+ Series. Assessment of the Accuracy of the SIBFA Anisotropic Polarizable Molecular Mechanics Potential

Stacking of guanine quartets (GQs) can trigger the formation of DNA or RNA quadruple helices, which play numerous biochemical roles. The GQs are stabilized by alkali cations, mainly K⁺ and Na⁺, which can reside in, or channel through, the central axis of the GQ stems. Further, ion conduction through GQ wires can be leveraged for nanochemistry applications. G-quadruplex systems have been extensively studied by classical molecular dynamics (MD) simulations using pair-additive force fields or by quantum-chemical (QC) calculations. However, the non-polarizable force fields are very approximate, while QC calculations lack the necessary sampling. Thus, ultimate description of GQ systems would require long-enough simulations using advanced polarizable molecular mechanics (MM). However, to perform such calculations, it is first mandatory to evaluate the method's accuracy using benchmark QC. We report such an evaluation for SIBFA polarizable MM, bearing on the channeling (movement) of an alkali cation (Li⁺, Na⁺, K⁺, or Rb⁺) along the axis of two stacked G quartets interacting with either one or two ions. The QC energy profiles display markedly different features depending upon the cation but can be retrieved in the majority of cases by the SIBFA profiles. An appropriate balance of first-order (electrostatic and short-range repulsion) and second-order (polarization, charge-transfer, and dispersion) contributions within ΔE is mandatory. With two cations in the channel, the relative weights of the second-order contributions increase steadily upon increasing the ion size. In the G8 complexes with two K⁺ or two Rb⁺ cations, the sum of polarization and charge-transfer exceeds the first order terms for all ion positions.

Solvent effects on the photochemistry of 4-aminoimidazole-5-carbonitrile, a prebiotically plausible precursor of purines

4-Aminoimidazole-5-carbonitrile (AICN) was suggested as a prebiotically plausible precursor of purine nucleobases and nucleotides. Although it can be formed in a sequence of photoreactions, AICN is immune to further irradiation with UV-light. We present state-of-the-art multi-reference quantum-chemical calculations of potential energy surface cuts and conical intersection optimizations to explain the molecular mechanisms underlying the photostability of this compound. We have identified the N-H bond stretching and ring-puckering mechanisms that should be responsible for the photochemistry of AICN in the gas phase. We have further considered the photochemistry of AICN-water clusters, while including up to six explicit water molecules. The calculations reveal charge transfer to solvent followed by formation of an H3O(+) cation, both of which occur on the (1)πσ* hypersurface. Interestingly, a second proton transfer to an adjacent water molecule leads to a (1)πσ*/S0 conical intersection. We suggest that this electron-driven proton relay might be characteristic of low-lying (1)πσ* states in chromophore-water clusters. Owing to its nature, this mechanism might also be responsible for the photostability of analogous organic molecules in bulk water.

From formamide to RNA, the path is tenuous but continuous

Reactions of formamide (NH2COH) in the presence of catalysts of both terrestrial and meteoritic origin yield, in plausible and variegated conditions, a large panel of precursors of (pre)genetic and (pre)metabolic interest. Formamide chemistry potentially satisfies all of the steps from the very initial precursors to RNA. Water chemistry enters the scene in RNA non-enzymatic synthesis and recombination.

Molecular mechanism of diaminomaleonitrile to diaminofumaronitrile photoisomerization: an intermediate step in the prebiotic formation of purine nucleobases

The photoinduced isomerization of diaminomaleonitrile (DAMN) to diaminofumaronitrile (DAFN) was suggested to play a key role in the prebiotically plausible formation of purine nucleobases and nucleotides. In this work we analyze two competitive photoisomerization mechanisms on the basis of state-of-the-art quantum-chemical calculations. Even though it was suggested that this process might occur on the triplet potential-energy surface, our results indicate that the singlet reaction channel should not be disregarded either. In fact, the peaked topography of the S1 /S0 conical intersection suggests that the deexcitation should most likely occur on a sub-picosecond timescale and the singlet photoisomerization mechanism might effectively compete even with a very efficient intersystem crossing. Such a scenario is further supported by the relatively small spin-orbit coupling of the S1 and T2 states in the Franck-Condon region, which does not indicate a very effective triplet bypass for this photoreaction. Therefore, we conclude that the triplet reaction channel in DAMN might not be as prominent as was previously thought.

Structural and energetic factors controlling the enantioselectivity of dinucleotide formation under prebiotic conditions

Recently, it has been reported that the montmorillonite-catalyzed oligomerization of activated nucleotides exhibits remarkable enantioselectivity. In the current paper we investigate the structures and intrinsic energies of homochiral and heterochiral cyclic dinucleotides by means of accurate quantum chemical calculations in gas-phase and in bulk water. The steric effect of the clay is represented with geometrical constraints. Our computations reveal that the heterochiral dimer geometries are systematically less stable than their homochiral counterparts due to steric clashes inside the sugar-phosphate ring geometry. Thus we suggest that the homochiral selectivity observed in the cyclic dinucleotide formation in confined spaces may arise from the energetic destabilization of the heterochiral ring geometries as compared to their homochiral analogues. In the present paper we provide the first model of the 3D structure of d,l cyclic dinucleotides, which until now has eluded experimental observation.

How does hydroxyl introduction influence the double helical structure: the stabilization of an altritol nucleic acid:ribonucleic acid duplex

Altritol nucleic acids (ANAs) are a promising new tool in the development of artificial small interfering ribonucleic acids (siRNAs) for therapeutical applications. To mimic the siRNA:messenger RNA (mRNA) interactions, the crystal structure of the ANA:RNA construct a(CCGUAAUGCC-P):r(GGCAUUACGG) was determined to 1.96 Å resolution which revealed the hybrid to form an A-type helix. As this A-form is a major requirement in the RNAi process, this crystal structure confirms the potential of altritol-modified siRNAs. Moreover, in the ANA strands, a new type of intrastrand interactions was found between the O2' hydroxyl group of one residue and the sugar ring O4' atom of the next residue. These interactions were further investigated by quantum chemical methods. Besides hydration effects, these intrastrand hydrogen bonds may also contribute to the stability of ANA:RNA duplexes.

Noncanonical hydrogen bonding in nucleic acids. Benchmark evaluation of key base-phosphate interactions in folded RNA molecules using quantum-chemical calculations and molecular dynamics simulations

RNA molecules are stabilized by a wide range of noncanonical interactions that are not present in DNA. Among them, the recently classified base-phosphate (BPh) interactions belong to the most important ones. Twelve percent of nucleotides in the ribosomal crystal structures are involved in BPh interactions. BPh interactions are highly conserved and provide major constraints on RNA sequence evolution. Here we provide assessment of the energetics of BPh interactions using MP2 computations extrapolated to the complete basis set of atomic orbitals and corrected for higher-order electron correlation effects. The reference computations are compared with DFT-D and DFT-D3 approaches, the SAPT method, and the molecular mechanics force field. The computations, besides providing the basic benchmark for the BPh interactions, allow some refinements of the original classification, including identification of some potential doubly bonded BPh patterns. The reference computations are followed by analysis of some larger RNA fragments that consider the context of the BPh interactions. The computations demonstrate the complexity of interaction patterns utilizing the BPh interactions in real RNA structures. The BPh interactions are often involved in intricate interaction networks. We studied BPh interactions of protonated adenine that can contribute to catalysis of hairpin ribozyme, the key BPh interaction in the S-turn motif of the sarcin-ricin loop, which may predetermine the S-turn topology and complex BPh patterns from the glmS riboswitch. Finally, the structural stability of BPh interactions in explicit solvent molecular dynamics simulations is assessed. The simulations well preserve key BPh interactions and allow dissection of structurally/functionally important water-meditated BPh bridges, which could not be considered in earlier bioinformatics classification of BPh interactions.

Theoretical modeling on the kinetics of the arsenate-ester hydrolysis: implications to the stability of As-DNA

Quantum chemical calculations reveal that neither steric hindrance nor less polar solvent medium is able to reduce the otherwise high hydrolysis rate of arsenate-esters. These results question the stability of As-DNA not only in aqueous but also in non-aqueous environments.

On the role of the cis Hoogsteen:sugar-edge family of base pairs in platforms and triplets-quantum chemical insights into RNA structural biology

Base pairs belonging to the cis Hoogsteen:sugar-edge (H:S) family play important structural roles in folded RNA molecules. Several of these are present in internal loops, where they are involved in interactions leading to planar dinucleotide platforms which stabilize higher order structures such as base triplets and quartets. We report results of analysis of 30 representative examples spanning 16 possible base pair combinations, with several of them showing multimodality of base pairing geometry. The geometries of 23 of these base pairs were modeled directly from coordinates extracted from RNA crystal structures. The other seven were predicted structures which were modeled on the basis of observed isosteric analogues. After appropriate satisfaction of residual valencies, these structures were relaxed using the B3LYP/6-31G(d,p) method and interaction energies were derived at the RIMP2/aug-cc-pVDZ level of theory. The geometries for each of the studied base pairs have been characterized in terms of the number and nature of H-bonds, rmsd values observed on optimization, base pair geometrical parameters, and sugar pucker analysis. In addition to its evaluation, the nature of intermolecular interaction in these complexes was also analyzed using Morokuma decomposition. The gas phase interaction energies range between -5.2 and -20.6 kcal/mol and, in contrast to the H:S trans base pairs, show enhanced relative importance of the electron correlation component, indicative of the greater role of dispersion energy in stabilization of these base pairs. The rich variety of hydrogen bonding pattern, involving the flexible sugar edge, appears to hold the key to several features of structural motifs, such as planarity and propensity to participate in triplets, observed in this family of base pairs. This work explores these aspects by integrating database analysis, and detailed base pairing geometry analysis at the atomistic level, with ab initio computation of interaction energies. The study, involving alternative classification of base pairs and triplets, provides insights into intrinsic properties of these base pairs and their possible structural and functional roles.

Classification and energetics of the base-phosphate interactions in RNA

Structured RNA molecules form complex 3D architectures stabilized by multiple interactions involving the nucleotide base, sugar and phosphate moieties. A significant percentage of the bases in structured RNA molecules in the Protein Data Bank (PDB) hydrogen-bond with phosphates of other nucleotides. By extracting and superimposing base-phosphate (BPh) interactions from a reduced-redundancy subset of 3D structures from the PDB, we identified recurrent phosphate-binding sites on the RNA bases. Quantum chemical calculations were carried out on model systems representing each BPh interaction. The calculations show that the centers of each cluster obtained from the structure superpositions correspond to energy minima on the potential energy hypersurface. The calculations also show that the most stable phosphate-binding sites occur on the Watson-Crick edge of guanine and the Hoogsteen edge of cytosine. We modified the 'Find RNA 3D' (FR3D) software suite to automatically find and classify BPh interactions. Comparison of the 3D structures of the 16S and 23S rRNAs of Escherichia coli and Thermus thermophilus revealed that most BPh interactions are phylogenetically conserved and they occur primarily in hairpin, internal or junction loops or as part of tertiary interactions. Bases that form BPh interactions, which are conserved in the rRNA 3D structures are also conserved in homologous rRNA sequence alignments.

Trans Hoogsteen/sugar edge base pairing in RNA. Structures, energies, and stabilities from quantum chemical calculations

Trans Hoogsteen/sugar edge (H/SE) RNA base pairs form one of the six families of RNA base pairs that utilize the 2'-hydroxyl group of ribose for base pairing and play key roles in stabilizing folded RNA molecules. Here, we provide a detailed quantum chemical characterization of intrinsic structures and interaction energies of this base pair family, along with the evaluation of solvent screening effects by a continuum solvent approach. We report DFT-optimized geometries and MP2 interaction energies for all 10 crystallographically identified members of the family, for a representative set of them, using complete basis set extrapolation. For 6 of the 10 base pairs, we had to apply geometric constraints to keep the geometries relevant to RNA. We confirm that the remaining, hitherto undetected, possible members of this family do not have appropriate steric features required to establish stable base pairing in the trans H/SE fashion. The interaction patterns in the trans H/SE family are highly diverse, with gas-phase interaction energies in the range from -1 to -17 kcal/mol. Except for the C/rC and G/rG trans H/SE base pairs, the interaction energy is roughly evenly distributed between the HF and correlation components. Thus, in the trans H/SE base pairs, the relative importance of electron correlation is noticeably smaller than in the cis WC/SE or cis and trans SE/SE base pairs, but still larger than in canonical base pairs. The trans H/SE A/rG base pair is the intrinsically most stable member of this family. This base pair is also known as the sheared AG base pair and belongs to the most prominent set of RNA base pairs utilized in molecular building blocks of functional RNAs. For all trans H/SE base pairs that we identified, in addition to conventional base pairing, viable alternative structures were stabilized by amino-acceptor interactions. In the QM calculations, these amino-acceptor complexes appear to be equally as stable as those with common H-bonds, and more importantly, the switch to amino-acceptor interaction does not require any significant geometrical rearrangement of the base pairs. Such interactions are worthy of further investigations, as X-ray crystallography cannot unambiguously distinguish between conventional and amino-acceptor interactions involving the 2'-hydroxyl group, formation of such interactions is usually not considered, and molecular modeling force fields do not include such interactions properly as a result of neglect of aminogroup pyramidalization.

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.

Theoretical study of the scalar coupling constants across the noncovalent contacts in RNA base pairs: the cis- and trans-watson-crick/sugar edge base pair family

The structure and function of RNA molecules are substantially affected by non-Watson-Crick base pairs actively utilizing the 2'-hydroxyl group of ribose. Here we correlate scalar coupling constants across the noncovalent contacts calculated for the cis- and trans-WC/SE (Watson-Crick/sugar edge) RNA base pairs with the geometry of base to base and sugar to base hydrogen bond(s). 23 RNA base pairs from the 32 investigated were found in RNA crystal structures, and the calculated scalar couplings are therefore experimentally relevant with regard to the binding patterns occurring in this class of RNA base pairs. The intermolecular scalar couplings 1hJ(N,H), 2hJ(N,N), 2hJ(C,H), and 3hJ(C,N) were calculated for the N-H...N and N-H...O=C base to base contacts and various noncovalent links between the sugar hydroxyl and RNA base. Also, the intramolecular 1J(N,H) and 2J(C,H) couplings were calculated for the amino or imino group of RNA base and the ribose 2'-hydroxyl group involved in the noncovalent interactions. The calculated scalar couplings have implications for validation of local geometry, show specificity for the amino and imino groups of RNA base involved in the linkage, and can be used for discrimination between the cis- and trans-WC/SE base pairs. The RNA base pairs within an isosteric subclass of the WC/SE binding patterns can be further sorted according to the scalar couplings calculated across different local noncovalent contacts. The effect of explicit water inserted in the RNA base pairs on the magnitude of the scalar couplings was calculated, and the data for discrimination between the water-inserted and direct RNA base pairs are presented. The calculated NMR data are significant for structural interpretation of the scalar couplings in the noncanonical RNA base pairs.

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.

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.

Sugar edge/sugar edge base pairs in RNA: stabilities and structures from quantum chemical calculations

Cis and trans sugar edge/sugar edge (SE/SE) binding patterns are essential building units of RNAs. For example, SE/SE interactions form the A-minor motifs, the most important tertiary interaction type in functional RNAs. This study provides an in-depth structure and stability analysis for these two base pair families. Gas-phase-optimized geometries are reported for 12 cis and 7 trans SE/SE base pairs and contrasted to their X-ray counterparts. Interaction energies are computed at the RIMP2 level of theory using the density-functional-theory-optimized geometries. There is a good overall agreement between the optimized and X-ray geometries of the cis SE/SE base pairs. In contrast, only three of the seven trans SE/SE binding patterns could be optimized without a significant distortion of the X-ray geometry. Note, however, that many SE/SE base pairs participate in broader networks of interactions; thus it is not surprising to see some of them to deviate from the X-ray geometry in a complete isolation. Computed interaction energies reveal that all 12 known cis SE/SE binding patterns are very stable. Among the trans SE/SE binding patterns, only the rG/rG, rG/rC, and rA/rG base pairs are sufficiently stable in the crystal geometry. Prediction has been made for some structures not yet detected by crystallography, namely, cis rC/rC, rG/rC, rG/rU, and rU/rU and trans rG/rA base pairs. Interestingly, the new cis SE/SE binding patterns are not necessarily isosteric with the remaining 12 members of this family. The trans rG/rA base pair represents a viable option for base pairing in RNA to be identified by future X-ray studies. In a complete lack of structural information, prediction of other unknown members of the trans SE/SE family was not attempted. Analysis of the interaction energies shows a very large electron correlation component of the interaction energy, pointing at the elevated role of dispersion energy as compared to other types of base pairs. This likely is profitable for stabilization of SE/SE binding patterns in polar environments and could be one of the reasons why the A-minor motif is the leading type of tertiary interactions in RNAs.

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.

Theoretical Study on the Structure, Stability, and Electronic Properties of the Guanine−Zn−Cytosine Base Pair in M-DNA

M-DNA is a type of metalated DNA that forms at high pH and in the presence of Zn, Ni, and Co, with the metals placed in between each base pair, as in G-Zn-C. Experiments have found that M-DNA could be a promising candidate for a variety of nanotechnological applications, as it is speculated that the metal d-states enhance the conductivity, but controversy still clouds these findings. In this paper, we carry out a comprehensive ab initio study of eight G-Zn-C models in the gas phase to help discern the structure and electronic properties of Zn-DNA. Specifically, we study whether a model prefers to be planar and has electronic properties that correlate with Zn-DNA having a metallic-like conductivity. Out of all the studied models, there is only one which preserves its planarity upon full geometry optimization. Nevertheless, starting from this model, one can deduce a parallel Zn-DNA architecture only. This duplex would contain the imino proton, in contrast to what has been proposed experimentally. Among the nonplanar models, there is one that requires less than 8 kcal/mol to flatten (both in gas and solvent conditions), and we propose that it is a plausible model for building an antiparallel duplex. In this duplex, the imino proton would be replaced by Zn, in accordance with experimental models. Neither planar nor nonplanar models have electronic properties that correlate with Zn-DNA having a metallic-like conductivity due to Zn d-states. To understand whether density functional theory (DFT) can describe appropriately the electronic properties of M-DNAs, we have investigated the electronic properties of G-Co-C base pairs. We have found that when self-interaction corrections (SIC) are not included the HOMO state contains Co d-levels, whereas these levels are moved below the HOMO state when SIC are considered. This result indicates that caution should be exercised when studying the electronic properties of M-DNAs with functionals that do not account for strong electronic correlations.

Mechanism of Action of Anticancer Titanocene Derivatives: An Insight from Quantum Chemical Calculations

Titanocene derivatives exhibit high potential in the treatment of cisplatin-resistant tumor types. Density functional theory calculations were performed on the hydrated form of five drug candidates differing in the pendant arms attached to the aromatic rings. A qualitative correlation has been found between the experimentally measured anticancer activity of alkylammonium-functionalized titanocene derivatives and the computed free energy change of the proton-induced dissociation reaction of these compounds. The results indicate that differences in the cytotoxic activities could be related to the solvation properties of the protolysis products, whereas no correlation was found with gas-phase properties of these molecules. Contrary to the free energy change of the protolysis reaction, other molecular properties, such as the geometrical parameters or the binding energies of the cyclopentadienyl rings in solution, do not correlate with the in vitro cytotoxic activity of these drug candidates.

Molecular dynamics simulations of RNA kissing-loop motifs reveal structural dynamics and formation of cation-binding pockets

Explicit solvent molecular dynamics (MD) simulations were carried out for three RNA kissing-loop complexes. The theoretical structure of two base pairs (2 bp) complex of H3 stem-loop of Moloney murine leukemia virus agrees with the NMR structure with modest violations of few NMR restraints comparable to violations present in the NMR structure. In contrast to the NMR structure, however, MD shows relaxed intermolecular G-C base pairs. The core region of the kissing complex forms a cation-binding pocket with highly negative electrostatic potential. The pocket shows nanosecond-scale breathing motions coupled with oscillations of the whole molecule. Additional simulations were carried out for 6 bp kissing complexes of the DIS HIV-1 subtypes A and B. The simulated structures agree well with the X-ray data. The subtype B forms a novel four-base stack of bulged-out adenines. Both 6 bp kissing complexes have extended cation-binding pockets in their central parts. While the pocket of subtype A interacts with two hexacoordinated Mg2+ ions and one sodium ion, pocket of subtype B is filled with a string of three delocalized Na+ ions with residency times of individual cations 1-2 ns. The 6 bp complexes show breathing motions of the cation-binding pockets and loop major grooves.

Unique tertiary and neighbor interactions determine conservation patterns of Cis Watson-Crick A/G base-pairs

X-ray, phylogenetic and quantum chemical analysis of molecular interactions and conservation patterns of cis Watson-Crick (W.C.) A/G base-pairs in 16S rRNA, 23S rRNA and other molecules was carried out. In these base-pairs, the A and G nucleotides interact with their W.C. edges with glycosidic bonds oriented cis relative to each other. The base-pair is stabilised by two hydrogen bonds, the C1'-C1' distance is enlarged and the G(N2) amino group is left unpaired. Quantum chemical calculations show that, in the absence of other interactions, the unpaired amino group is substantially non-planar due to its partial sp(3) pyramidalization, while the whole base-pair is internally propeller twisted and very flexible. The unique molecular properties of the cis W.C. A/G base-pairs make them distinct from other base-pairs. They occur mostly at the ends of canonical helices, where they serve as interfaces between the helix and other motifs. The cis W.C. A/G base-pairs play crucial roles in natural RNA structures with salient sequence conservation patterns. The key contribution to conservation is provided by the unpaired G(N2) amino group that is involved in a wide range of tertiary and neighbor contacts in the crystal structures. Many of them are oriented out of the plane of the guanine base and utilize the partial sp(3) pyramidalization of the G(N2). There is a lack of A/G to G/A covariation, which, except for the G(N2) position, would be entirely isosteric. On the contrary, there is a rather frequent occurrence of G/A to G/U covariation, as the G/U wobble base-pair has an unpaired amino group in the same position as the cis W.C. G/A base-pair. The cis W.C. A/G base-pairs are not conserved when there is no tertiary or neighbor interaction. Obtaining the proper picture of the interactions and phylogenetic patterns of the cis W.C. A/G base-pairs requires a detailed analysis of the relation between the molecular structures and the energetics of interactions at a level of single H-bonds and contacts.

The influence of the thymine C5 methyl group on spontaneous base pair breathing in DNA

Sequences of four or more AT base pairs without a 5'-TA-3' step, so-called A-tracts, influence the global properties of DNA by causing curvature of the helix axis if phased with the helical repeat and also influence nucleosome packaging. Hence it is interesting to understand this phenomenon on the molecular level, and numerous studies have been devoted to investigations of dynamical and structural features of A-tract DNA. It was early observed that anomalously slow base pair-opening kinetics were a striking physical property unique to DNA A-tracts (Leroy, J. L., Charretier, E., Kochoyan, M., and Gueron, M. (1988) Biochemistry 27, 8894-8898). Furthermore, a strong correlation between DNA curvature and anomalously slow base pair-opening dynamics was found. In the present work it is shown, using imino proton exchange measurements by NMR spectroscopy that the main contribution to the dampening of the base pair-opening fluctuations in A-tracts comes from the C5 methylation of the thymine base. Because the methyl group has been shown to have a very limited effect on the DNA curvature as well as the structure of the DNA helix, the thymine C5 methyl group stabilizes the helix directly. Empirical potential energy calculations show that methylation of the tract improves the stacking energy of a base pair with its neighbors in the tract by 3-4 kcal/mol.

Loss of Hoogsteen pairing ability upon N1 adenine platinum binding

Chloroform- and Freon-soluble mixed thymine, adenine complexes trans-[Pt(MeNH(2))(2)(ChmT-N3)(ChmA-N1)]NO(3) (2) and trans-[Pt(MeNH(2))(2)(ChmT-N3)(TBDMS-ado-N1)]BF(4) (3) (ChmT = anion of 1-cyclohexylmethylthymine ChmTH, ChmA = 9-cyclohexylmethyladenine, TBDMS-ado = 2',3',5'-tri-tert-butyldimethylsilyladenosine) have been prepared and characterized to study their propensity to undergo Hoogsteen and/or reversed Hoogsteen pairing in solution with free ChmTH and free 3',5'-diacetyl-2'-deoxyuridine, respectively. No Hoogsteen or reversed Hoogsteen pairing between 2 and ChmT takes place in CDCl(3). In Freon, partial H bonding between N1 platinated TBDMS-ado and 3',5'-diacetyl-2'-deoxyuridine as well as its [3-(15)N] labeled analogue is unambiguously observed only below 150 K. Comparison of (1)J ((15)N-(1)H) coupling constants of 3',5'-diacetyl-2'-deoxyuridine involved in Hoogsteen pairing with free and N1 platinated adenine suggests that the interaction is inherently weaker in the case of platinated adenine. To better understand the complete absence of hydrogen bonding between the ChmA ligand in 2 and free ChmTH, ab initio calculations (gas phase, 0 K) have been carried out for Hoogsteen pairs involving adenine (A) and thymine (T), as well as simplified analogues of 2 and T, both in the presence and absence of counteranions. The data strongly suggest that reduction of the effective positive charge of the heavy metal ion Pt(2+) by counterions diminishes interaction energies. With regard to mixtures of 2 and ChmTH in chloroform, this implies that ion pair formation between the cation of 2 and NO(3)(-) may be responsible for the lack of any measurable Hoogsteen pairing in this solvent.

Intercalators. 1. Nature of stacking interactions between intercalators (ethidium, daunomycin, ellipticine, and 4',6-diaminide-2-phenylindole) and DNA base pairs. Ab initio quantum chemical, density functional theory, and empirical potential study

Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.

Protonation of platinated adenine nucleobases. Gas phase vs condensed phase picture

Protonation of adenine carrying a Pt(II) moiety either at N7, N3, or N1 is possible in solution, but the site of protonation is influenced by the location of the Pt(II) electrophile and to some extent also by the overall charge of the metal entity (+2, +1, 0, -1), hence the other ligands (NH(3), Cl(-), OH(-)) bound to Pt(II). Quantum chemical calculations based on density functional theory (DFT) have been carried out for intrinsic protonation energies of adenine complexes carrying the following Pt(II) species at either of the three ring N atoms: [Pt(NH(3))(3)](2+) (1), trans- [Pt(NH(3))(2)Cl](+) (2a), cis-[Pt(NH(3))(2)Cl](+) (2b), trans-[Pt(NH(3))(2)Cl(2)] (3a), cis-[Pt(NH(3))Cl(2)] (3b), [PtCl(3)](-) (4), trans-[Pt(NH(3))(2)OH](+) (5a), cis-[Pt(NH(3))(2)(OH)](+) (5b), trans-[Pt(NH(3))(OH)(2)] (6a), cis-[Pt(NH(3))(OH)(2)] (6b), and [Pt(OH)(3)](-) (7). The data have been compared with results derived from solution studies (water) and X-ray crystallography, whenever available. The electrostatic effects associated with the charge of the metal entity have the major influence on the calculated intrinsic (gas phase) proton affinities, unlike the condensed phase data. Nevertheless, the relative gas phase trends correlate surprisingly well with condensed phase data; i.e., variation of the pK(a) values measured in solution is consistent with the calculated gas phase protonation energies. In addition to a systematic study of the ring proton affinities, proton transfer processes within the platinated adenine species were often observed when investigating Pt adducts with OH(-) ligands, and they are discussed in more detail. To the best of our knowledge, this is the first study attempting to find a systematic correlation between gas phase and condensed phase data on protonation of metalated nucleobases. The gas phase data provide a very useful complement to the condensed phase and X-ray experiments, showing that the gas phase studies are capable of valuable predictions and contribute to our understanding of the solvent and counterion effects on metal-assisted proton shift processes.

Cation-pi and amino-acceptor interactions between hydrated metal cations and DNA bases. A quantum-chemical view

Cation-pi interactions between cytosine and hexahydrated cations have been characterized using ab initio method with inclusion of electron correlation effects, assuming idealized and crystal geometries of the interacting species. Hydrated metal cations can interact with nucleobases in a cation-pi manner. The stabilization energy of such complexes would be large and comparable to the one for cation-pi complex with benzene. Further, polarized water molecules belonging to the hydration shell of the cation are capable to form a strong hydrogen bond interaction with the nitrogen lone electron pair of the amino groups of bases and enforce a pronounced sp3 pyramidalization of the nucleobase amino groups. However, in contrast to the benzene-cation complexes, the cation-pi configurations are highly unstable for a nucleobase since the conventional in plane binding of hydrated cations to the acceptor sites on the nucleobase is strongly preferred. Thus, a cation-pi interaction with a nucleobase can occur only if the position of the cation is locked above the nucleobase plane by another strong interaction. This indeed can occur in biopolymers and may have an effect on the local DNA architecture. Nevertheless, nucleobases have no intrinsic propensity to form cation-pi interactions.