The elusive atomic rationale for DNA base pair stability

A systematic analysis of the electrostatic interaction between 27 natural DNA base pairs was carried out, based on ab initio correlated wave functions and the topology of the electron density. Using high rank multipole moments we show that the atomic partitioning of the interaction energy contains many substantial contributions between distant atoms. Profiles of cumulative energy versus internuclear distance show large fluctuations and provide an electrostatic fingerprint of the partitioning of interaction energy in a complex. A quantified comparison between each pair of energy profiles, one for each base pair, makes clear that there is no correlation between the total base pair interaction energy and the shape of the profile. In other words, base pairs with similar interaction energy are not stable for the same reasons in terms of atomic partitioning. In summary, simple rules to rationalize the pattern of energetic stability of naturally occurring base pairs in terms of subsets of atoms are elusive. Our work cautions against inappropriate use of Jorgensen's secondary interaction hypothesis.     

Contact pair dynamics during folding of two small proteins: chicken villin head piece and the Alzheimer protein beta-amyloid

The folding of an extended protein to its unique native state requires establishment of specific, predetermined, often distant, contacts between amino acid residue pairs. The dynamics of contact pair formation between various hydrophobic residues during folding of two different small proteins, the chicken villin head piece (HP-36) and the Alzheimer protein beta-amyloid (betaA-40), are investigated by Brownian dynamics (BD) simulations. These two proteins represent two very different classes-HP-36 being globular while betaA-40 is nonglobular, stringlike. Hydropathy scale and nonlocal helix propensity of amino acids are used to model the complex interaction potential among the various amino acid residues. The minimalistic model we use here employs a connected backbone chain of atoms of equal size while an amino acid is attached to each backbone atom as an additional atom of differing sizes and interaction parameters, determined by the characteristics of each amino acid. Even for such simple models, we find that the low-energy structures obtained by BD simulations of both the model proteins mimic the native state of the real protein rather well, with a best root-mean-square deviation of 4.5 A for HP-36. For betaA-40 (where a single well-defined structure is not available), the simulated structures resemble the reported ensemble rather well, with the well-known beta-bend correctly reproduced. We introduce and calculate a contact pair distance time correlation function, C(P) (ij)(t), to quantify the dynamical evolution of the pair contact formation between the amino acid residue pairs i and j. The contact pair time correlation function exhibits multistage dynamics, including a two stage fast collapse, followed by a slow (microsecond long) late stage dynamics for several specific pairs. The slow late stage dynamics is in accordance with the findings of Sali et al. Analysis of the individual trajectories shows that the slow decay is due to the attempt of the protein to form energetically more favorable pair contacts to replace the less favorable ones. This late stage contact formation is a highly cooperative process, involving participation of several pairs and thus entropically unfavorable and expected to face a large free energy barrier. This is because any new pair contact formation among hydrophobic pairs will require breaking of several contacts, before the favorable ones can be formed. This aspect of protein folding dynamics is similar to relaxation in glassy liquids, where also alpha relaxation requires highly cooperative process of hopping. The present analysis suggests that waiting time for the necessary pair contact formation may obey the Poissonian distribution. We also study the dynamics of Forster energy transfer during folding between two tagged amino acid pairs. This dynamics can be studied by fluorescence resonance energy transfer (FRET). It is found that suitably placed donor-acceptor pairs can capture the slow dynamics during folding. The dynamics probed by FRET is predicted to be nonexponential.     

Multicanonical Monte Carlo calculation of the free‐energy map of the base–amino acid interaction

The specific interactions between base pairs and amino acids were studied by the multicanonical Monte Carlo method. We sampled numerous interaction configurations and side‐chain conformations of the amino acid by the multicanonical algorithm, and calculated the free energies of the interactions between an amino acid at given C_α positions and a fixed base pair. The contour maps of free energy derived from this calculation represent the preferred C_α position of the amino acid around the base, and these maps of various combinations of bases and amino acids can be used to quantify the specificity of intrinsic base–amino acid interactions. Similarly, enthalpy and entropy maps will provide further details of the specific interactions. We have also calculated the free‐energy map of the orientations of the C_α C_β bond vector, which indicates the preferential orientation of the amino acid against the base. We compared the results obtained by the multicanonical method with those of the exhaustive sampling and canonical Monte Carlo methods. The free‐energy map of the base–amino acid interaction obtained by the multicanonical simulation method was nearly identical to the accurate result derived from the exhaustive sampling method. This indicates that a single multicanonical Monte Carlo simulation can produce an accurate free‐energy map. Multicanonical Monte Carlo sampling produced free‐energy maps that were more accurate than those produced by canonical Monte Carlo sampling. Thus, the multicanonical Monte Carlo method can serve as a powerful tool for estimating the free‐energy landscape of base–amino acid interactions and for elucidating the mechanism by which amino acids of proteins recognize particular DNA base pairs. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 954–962, 2000     

The correspondence of detachment characteristics and acid-base interaction measures

A quantitative estimate for the intensity of interphase acid-base interactions is proposed based on the acid and base parameters of the free surface energy of polymer adhesives and substrates of a different nature. An interaction between acid-base and strength characteristics of adhesive compounds is established. An increase in the adhesive strength with an increase in the difference between the surface acid and base properties of connected materials is revealed for a number of compounds (polyolefine-metal, rubber mixtures-brass, and Thiokol compounds-glass).     

Intra- and intermolecular interactions between cyclic-AMP receptor protein and DNA: ab initio fragment molecular orbital study

The ab initio fragment molecular orbital (FMO) calculations were performed for the cAMP receptor protein (CRP) complexed with a cAMP and DNA duplex to elucidate their sequence-specific binding and the stability of the DNA duplex, as determined by analysis of their inter- and intramolecular interactions. Calculations were performed with the AMBER94 force field and at the HF and MP2 levels with several basis sets. The interfragment interaction energies (IFIEs) were analyzed for interactions of CRP-cAMP with each base pair, DNA duplex with each amino acid residue, and each base pair with each residue. In addition, base-base interactions were analyzed including hydrogen bonding and stacking of DNA. In the interaction between DNA and CRP-cAMP, there was a significant charge transfer (CT) from the DNA to CRP, and this CT interaction played an important role as well as the electrostatic interactions. It is necessary to apply a quantum mechanical approach beyond the "classical" force-field approach to describe the sequence specificity. In the DNA intramolecular interaction, the dispersion interactions dominated the stabilization of the base-pair stacking interactions. Strong, attractive 1,2-stacking interactions and weak, repulsive 1,3-stacking interactions were observed. Comparison of the intramolecular interactions of free and complexed DNA revealed that the base-pairing interactions were stronger, and the stacking interactions were weaker, in the complexed structure. Therefore, the DNA duplex stability appears to change due to both the electrostatic and the CT interactions that take place under conditions of DNA-CRP binding.     

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.     

Strength of Calpha-H...O=C hydrogen bonds in transmembrane proteins

A large number of Calpha-H...O contacts are present in transmembrane protein structures, but contribution of such interactions to protein stability is still not well understood. According to previous ab initio quantum calculations, the stabilization energy of a Calpha-H...O contact is about 2-3 kcal/mol. However, experimental studies on two different Calpha-H...O hydrogen bonds present in transmembrane proteins lead to conclusions that one contact is only weakly stabilizing and the other is not even stabilizing. We note that most previous computational studies were on optimized geometries of isolated molecules, but the experimental measurements were on those in the structural context of transmembrane proteins. In the present study, 263 Calpha-H...O=C contacts in alpha-helical transmembrane proteins were extracted from X-ray crystal structures, and interaction energies were calculated with quantum mechanical methods. The average stabilization energy of a Calpha-H...O=C interaction was computed to be 1.4 kcal/mol. About 13% of contacts were stabilizing by more than 3 kcal/mol, and about 11% were destabilizing. Analysis of the relationships between energy and structure revealed four interaction patterns: three types of attractive cases in which additional Calpha-H...O or N-H...O contact is present and a type of repulsive case in which repulsion between two carbonyl oxygen atoms occur. Contribution of Calpha-H...O=C contacts to protein stability is roughly estimated to be greater than 5 kcal/mol per helix pair for about 16% of transmembrane helices but for only 3% of soluble protein helices. The contribution would be larger if Calpha-H...O contacts involving side chain oxygen were also considered.     

Surface characterization of industrial fibers with inverse gas chromatography

Inverse gas chromatography (IGC) was applied for the determination of the surface characteristics of Tenax carbon fibers and Akzo Nobel Twaron fibers. Furthermore, IGC procedures for the determination of dispersive and acid-base interactions were validated. The data show that too high values for the dispersive component of the surface energy are obtained when the adsorption area occupied by a single adsorbed n-alkane molecule is estimated from parameters of the corresponding liquid. Comparable values are obtained when the Doris-Gray methodology (area per methylene unit) or measured probe areas are employed. For the fibers studied in this work meaningful Gibbs energy values of the acid-base interaction were only obtained with the polarizability approach. When the dispersive interaction of the polar probes with the fiber surface was scaled to the n-alkane interaction via surface tension, the boiling point, or the vapor pressure of the probes often negative acid-base interaction energies were found. From the temperature dependence of the Gibbs energy, the enthalpy of the acid-base interactions of various probes with the carbon and Twaron aramid fibers was determined. However, from these enthalpy values no meaningful acid-base surface parameters could be obtained. Generally, the limited accuracy with which these parameters can be obtained make the usefulness of this procedure questionable. Also the Gibbs energy data of acid-base interaction can provide a qualitative basis to classify the acidity-basicity of the fiber surface. This latter approach requires only a limited data set and is sufficiently rapid to enable the use of IGC as a screening tool for fibers at a production site. For several polar probes significant concentration effects on carbon fibers were observed. At very low probe loadings the interaction with the fiber surface suddenly increases. This effect is caused by the heterogeneity of the interaction energy of the active sites at the surface. A simple procedure to measure the adsorption isotherm at infinite dilution was developed. The determination of the concentration dependence of the interaction of an n-alkane, an acidic and a basic probe was incorporated in the IGC screening procedure of carbon fibers to monitor this heterogeneity.     

C–H...O and C–H...N interactions in RNA structures

Small molecule studies indicate that C–H...X interactions (X: O,N) constitute weak H-bonds. We have performed a comprehensive analysis of their occurrence and geometry in RNA structures. Here, we report on statistical properties of the total set of interactions identified and discuss selected motifs. The distance/angle distribution of all interactions exhibits an excluded region where the allowed C–H...X angle range increases with an increasing H...X distance. The preferred short C–H...X interactions in RNA are backbone-backbone contacts between neighbour nucleotides. Distance/angle distributions generated for various interaction types can be used for error recognition and modelling. The axial C2′(H)...O4′ and C5′(H)...O2′ interactions connect two backbone segments and form a seven-membered ring that is specific for RNA. An AA base pair with one standard H-bond and one C–H...N interaction has been identified in various structures. Despite the occurrence of short C–H...X contacts their free energy contribution to RNA stability remains to be assessed. Received: 17 May 1998 / Accepted: 4 August 1998 / Published online: 2 November 1998     

Ab initio molecular orbital study on the pairing and stacking interactions between nucleic acid bases in relation to the biological activities

Ab initio molecular orbital calculations were performed on the pairing and stacking interaction energies between bases in nucleic acids. Using these values we could explain the biologically important phenomena well. Thus the fact that O⁶-methylguanine (which is formed in small amounts) is more promutagenic than N₇-methylguanine (which is formed in larger amounts) could be explained by the difference in pairing interaction energies for these alkylguanines. To clarify the detailed mechanism of mutation induced by a base analogue (2-aminopurine) the interaction energy for the 2-aminopurine-cytosine pair was calculated by taking into account the tautomeric conversion of base and the base analogue. It was concluded that the base pair formed as an intermediate between the normal form of 2-aminopurine and the imino tautomer of cytosine has an important role in inducing the mutation by 2-aminopurine. The stacking interaction energy was found to be a principal factor in determining the conformation of nucleic acids, and it predicted the preference for the A-form or B-form of the deoxyoligonucleotides well. The stacking interaction energy was resolved into its components, and it was shown that electrostatic energy was base sequence-dependent, whereas the overall stability of the stacked polymers was largely dependent on the dispersion energy.     

Dual Recognition of Double-Stranded DNA by 2′-Aminoethoxy-Modified Oligonucleotides

Simultaneous interaction of the 2′-aminoethoxy-modified oligonucleotides with the phosphodiester backbone (shown on the right, A) and with the bases through Hoogsteen base contacts (B) is seen at each base-pair step of the duplex DNA target. The electrostatic interaction between the protonated amino group and the negatively charged phosphate group provides for a dramatic increase in the binding affinity and the association rate constant.     

6‐Aminopyridine‐3‐thiol

The title compound, C₅H₆N₂S, is a simple but novel pyridinethiol of pharmacological interest. The mol­ecule is planar. The crystal packing is dominated by hydro­phobic contacts and a pair of hydrogen‐bond interactions between the amino group of one pyridine mol­ecule and the ring N atom of a neighbouring base, stabilizing the structure.     

Hg(II) ion specifically binds with T:T mismatched base pair in duplex DNA

Metal-mediated base pair formation, resulting from the interaction between metal ions and artificial bases in oligonucleotides, has been developed for its potential application in nanotechnology. We have recently found that the T:T mismatched base pair binds with Hg(II) ions to generate a novel metal-mediated base pair in duplex DNA. The thermal stability of the duplex with the T-Hg-T base pair was comparable to that of the corresponding T:A or A:T. The novel T-Hg-T base pair involving the natural base thymine is more convenient than the metal-mediated base pairs involving artificial bases due to the lack of time-consuming synthesis. Here, we examine the specificity and thermodynamic properties of the binding between Hg(II) ions and the T:T mismatched base pair. Only the melting temperature of the duplex with T:T and not of the perfectly matched or other mismatched base pairs was found to specifically increase in the presence of Hg(II) ions. Hg(II) specifically bound with the T:T mismatched base pair at a molar ratio of 1:1 with a binding constant of 10(6) M(-1), which is significantly higher than that for nonspecific metal ion-DNA interactions. Furthermore, the higher-order structure of the duplex was not significantly distorted by the Hg(II) ion binding. Our results support the idea that the T-Hg-T base pair could eventually lead to progress in potential applications of metal-mediated base pairs in nanotechnology.     

Direct measurement of acid-base interaction energy at solid interfaces

We have studied acid-base interactions at solid-liquid and solid-solid interfaces using interface-sensitive sum frequency generation (SFG) spectroscopy. The shift of the sapphire hydroxyl peak in contact with several polar and nonpolar liquids and polymers was used to determine the interaction energy. The trend in the interaction energies cannot be explained by measuring only water contact angles. Molecular rearrangements at the sapphire interface, to maximize the interaction of the acid-base groups, play a dominant role, and these effects are not accounted for in the current theoretical models. These results provide important insights into understanding adhesion, friction, and wetting on solid interfaces.     

On the interaction of encapsulated pH indicator species within a silica matrix produced by three sol-gel routes

Solid acid-base sensors were prepared by encapsulating two pH indicators (brilliant yellow or acridine) within a silica matrix by the sol-gel method using three different routes: (1) non-hydrolytic, (2) acid catalyzed and (3) base catalyzed. The interactions of the silica-indicator with the resulting materials were then investigated by cyclic and differential pulse voltammetry. Complementary, ultraviolet-visible, photoacoustic spectroscopy was employed for the characterization of the interactions by monitoring the band shifts (bathochromic or hypsochromic, depending on the sol-gel route) between the neat pH indicators and those encapsulated within the silica network. Furthermore, X-ray photoelectron spectroscopy showed that the N 1s binding energy in brilliant yellow was shifted for the material resulting from the acid route. The electrochemical behavior and the pH indicator interactions with the silica network were dependent on the nature of the employed sol-gel route. For the sensors prepared with acridine, the interactions with the silica network took place through the nitrogen group from the pyridinic ring. For the brilliant yellow indicator, different behaviors were observed depending on the route, suggesting different processes during preparation or analysis. For the basic catalyzed and non-hydrolytic routes, it was not possible to assign a specific interaction. Nevertheless, it seemed that interactions might have taken place through the hydroxyl and/or sulphonic groups. Furthermore, for the brilliant yellow sensor prepared through the acid route, it was possible to show that the interaction probably or partially occurred through the azo groups.     

铂配合物与DNA碱基对间相互作用的理论研究

用量子化学方法研究一系列Pt(II)配合物作用于嘌呤碱基N7位点后对Watson-Crick碱基对AT、GC的影响. 计算结果显示铂配体与碱基对AT、GC的作用以静电作用为主，同时极化作用也是影响GC碱基对的重要因素. 静电作用极大地增强了铂化嘌呤碱基与嘧啶碱基间的相互作用，而嘌呤碱基与嘧啶碱基间作用与未铂化碱基对作用相近. Pd(II) 和 Ni(II)的相关研究得到类似的结果. 碱基对间氢键作用“二阶微扰能”分析结果与氢键强弱变化一致.     

A simple method for evaluation of the self-consistency of the radial distribution functions and the orientation of water molecules in the first coordination sphere from the results of molecular dynamics calculations

A simple method for determination of the angular orientation of water molecules in the first coordination sphere from the radial distribution functions is proposed. A comparative analysis of the ability of the model potentials of pair interaction to take into account the effects of manybody interactions (MBI) was performed. The responses of the model pair potentials to the MBI effects in the first and second coordination spheres were found to be poorly correlated with each other. It was concluded that it is necessary to derive a new analytical type of potential functions of pair interaction. Published inIzvestiya Akademii Nauk. Seriya Khimicheskaya. No. 11, pp. 1842–1846. November. 2000.     

Folding of a stable DNA motif involves a highly cooperative network of interactions

Hairpins are structural elements that play important roles in the folding and function of RNA and DNA. The extent of cooperativity in folding is an important aspect of the RNA folding problem. We reasoned that an investigation into the origin of cooperativity might be best carried out on a stable nucleic acid system with a limited number of interactions, such as a stable DNA hairpin loop. The stable d(cGNAg) hairpin loop motif (closing base pair in lower case; loop in upper case; N = A, C, G, or T) is stabilized through only three interactions: two loop-loop hydrogen bonds in a sheared GA base pair and a loop-closing base pair interaction. Herein, we investigate this network of interactions and test whether the loop-loop and loop-closing base pair interactions communicate. Thermodynamic measurements of nucleotide analogue substituted oligonucleotides were used to probe the additivity of the interactions. On the basis of double mutant cycles, all interactions were found to be nonadditive and interdependent, suggesting that loop-loop and loop-closing base pair interactions form in a highly cooperative manner. When double mutant cycles were repeated in the absence of the other interaction, nonadditivity was significantly reduced suggesting that coupling is indirect and requires all three interactions in order to be optimal. A cooperative network of interactions helps explain the structural and energetic bases of stability in certain DNA hairpins and paves the way for similar studies in more complex nucleic acid systems.     

Structure–affinity relationships for the binding of actinomycin D to DNA

Molecular models of the complexes between actinomycin D and 14 different DNA hexamers were built based on the X-ray crystal structure of the actinomycin–d(GAAGCTTC)2 complex. The DNA sequences included the canonical GpC binding step flanked by different base pairs, nonclassical binding sites such as GpG and GpT, and sites containing 2,6-diamino- purine. A good correlation was found between the intermolecular interaction energies calculated for the refined complexes and the relative preferences of actinomycin binding to standard and modified DNA. A detailed energy decomposition into van der Waals and electrostatic components for the interactions between the DNA base pairs and either the chromophore or the peptidic part of the antibiotic was performed for each complex. The resulting energy matrix was then subjected to principal component analysis, which showed that actinomycin D discriminates among different DNA sequences by an interplay of hydrogen bonding and stacking interactions. The structure–affinity relationships for this important antitumor drug are thus rationalized and may be used to advantage in the design of novel sequence-specific DNA-binding agents.