Solvent bridging sites in A- and B-DNA helices

Nucleotide hydration is important for the understanding of the stability of and the transitions between the different helical conformations of DNA. We have used energy minimization and geometric criteria in order to look for possible sites for solvent which can bridge more than one polar or charged atomic group on a nucleotide. Such bridging sites between phosphate groups have been seen experimentally and used to explain the A to B transition. We show that these phosphate bridging sites occur at energy minima around A-DNA but do not occur around B-DNA. We also find that there are further low energy bridging sites which depend on sequence and which enable the more economical hydration of the A form.     

A systematic study of patterns of hydration in nucleic acids:(I) guanine and cytosine

The hydration sites of guanine and cytosine are defined by examination of the crystal structures of bases, nucleosides, nucleotides, and three dinucleoside phosphate salts. The patterns of hydration for two guanine and cytosine containing oligonucleotides are then predicted. The relationship between these structural motifs and thermodynamic parameters is discussed.     

A Systematic Study of Patterns of Hydration in Nucleic Acids:(I) Guanine and Cytosine

Abstract

The hydration sites of guanine and cytosine are defined by examination of the crystal structures of bases, nucleosides, nucleotides, and three dinucleoside phosphate salts. The patterns of hydration for two guanine and cytosine containing oligonucleotides are then predicted. The relationship between these structural motifs and thermodynamic parameters is discussed.     

Hydration of RNA base pairs

The hydration patterns around the RNA Watson-Crick and non-Watson-Crick base pairs in crystals are analyzed and described. The results indicate that (i) the base pair hydration is mostly "in-plane"; (ii) eight hydration sites surround the Watson-Crick G-C and A-U base pairs, with five in the deep and three in the shallow groove, an observation which extends the characteristic isostericity of Watson-Crick pairs; (iii) while the hydration around G-C base pairs is well defined, the hydration around A-U base pairs is more diffuse; (iv) the hydration sites close to the phosphate groups are the best defined and the most recurrent ones; (v) a string of water molecules links the two shallow groove 2'-hydroxyl groups, and (vi) the water molecules fit into notches, the size and accessibility of which are almost as important as the number and strength of the hydrophilic groups lining the cavity. Residence times of water molecules at specific hydration sites, inferred from molecular dynamics simulations, are discussed in the light of present data.     

Function of the F1-motor (F1-ATPase) of ATP synthase by apolar-polar repulsion through internal interfacial water

THESIS: Within the structurally-confined internal aqueous cavity of the F1-motor of ATP synthase, function results from free energy changes that shift the balance between interfacial charge hydration and interfacial hydrophobic hydration. TRANSITION STATE DESCRIPTION: At the beta-P end of ADP x Mg occurs an inorganic phosphate, P(i). This P(i) resides at the base of a water-filled cleft that functions like an aperture to focus, into an aqueous chamber, a competition for hydration (an apolar-polar repulsion) between charged phosphate and hydrophobic surface of the gamma-rotor. Two means available for the phosphate and the hydrophobic surface to improve their hydration free energies are physically to separate by rotation of the gamma-rotor or chemically to combine P(i) with ADP to form less charged ATP. This proposal derives from calculated changes in Gibbs free energy for hydrophobic association of amino acid side chains and chemical modifications thereof and from experimentally demonstrated water-mediated repulsion between hydrophobic and charged sites that resulted from extensive studies on designed elastic-contractile model proteins.     

Ordered water structure around a B-DNA dodecamer: A quantitative study

The crystal structure of the double-helical B-DNA dodecamer of sequence C-G-C-G-A-A-T-T-C-G-C-G has been solved and refined independently in three forms: (1) the parent sequence at room temperature; (2) the same sequence at 16 K; and (3) the 9-bromo variant C-G-C-G-A-A-T-T^BrC-G-C-G at 7 °C in 60% (v/v) 2-methyl-2.4-pentanediol. The latter two structures show extensive hydration along the phosphate backbone, a feature that was invisible in the native structure because of high temperature factors (indicating thermal or static disorder) of the backbone atoms. Sixty-five solvent peaks are associated with the phosphate backbone, or an average of three per phosphate group. Nineteen other molecules form a first shell of hydration to base edge N and O atoms within the major groove, and 36 more are found in upper hydration layers. The latter tend to occur in strings or clusters spanning the major groove from one phosphate group to another. A single spermine molecule also spans the major groove. In the minor groove, the zig-zag spine of hydration that we believe to be principally responsible for stabilizing the B form of DNA is found in all three structures. Upper level hydration in the minor groove is relatively sparse, and consists mainly of strings of water molecules extending across the groove, with few contacts to the spine below. Sugar O-1′ atoms are closely associated with water molecules, but these are chiefly molecules in the spine, so the association may reflect the geometry of the minor groove rather than any intrinsic attraction of O-1′ atoms for hydration. The phosphate O-3′ and O-5′ atoms within the backbone chain are least hydrated of all, although no physical or steric impediment seems to exist that would deny access to these oxygen atoms by water molecules.     

Solvent interactions stabilising nucleic acid conformers

The transition of oligonucleotides from the B to the A conformation has been studied by the use of simple geometric calculations aimed at finding possible hydration sites which could stabilize these conformations. The method involves the classification of equally spaced grid points, surrounding the oligonucleotide, into groups depending on whether a water molecule, so placed, could form single, multiple or nil contacts to polar oligonucleotide atoms. The occurrence of the multiple, and therefore bridging, sites is more extensive for the 'A' than the 'B' conformation. Thus, more general evidence is presented in support of the economy of hydration hypothesis in which phosphate groups, in the low humidity 'A' form, have been seen to be bridged by single water molecules. Similar calculations for the 'Z' DNA conformation show a different preference for multiple bridging sites.     

Study of conformation-dependent isotherms of DNA hydration

By the method of gravimetry and piezomicrobalance native and denaturated DNA hydration isotherms are obtained. Their adsorptional ability difference observed is interpreted according to IR-data, which show that if sugar-phosphate skeleton hydration levels are similar for these conformations, their nitrous base hydration is strongly different. BET-fitting of isotherms obtained discovered significant water binding energy heterogeneity of native DNA hydration-active sites.     

Water and ion binding around r(UpA)12 and d(TpA)12 oligomers--comparison with RNA and DNA (CpG)12 duplexes

The structural and dynamic properties of the water and ion first coordination shell of the r(A-U) and d(A-T) base-pairs embedded within the r(UpA)12 and d(TpA)12 duplexes are described on the basis of two 2.4 ns molecular dynamics simulations performed in a neutralizing aqueous environment with 0.25 M added KCl. The results are compared to previous molecular dynamics simulations of the r(CpG)12 and d(CpG)12 structures performed under similar conditions. It can be concluded that: (i) RNA helices are more rigid than DNA helices of identical sequence, as reflected by the fact that RNA duplexes keep their initial A-form shape while DNA duplexes adopt more sequence-specific shapes. (ii) Around these base-pairs, the water molecules occupy 21 to 22 well-defined hydration sites, some of which are partially occupied by potassium ions. (iii) These hydration sites are occupied by an average of 21.9, 21.0, 20.1, and 19.8 solvent molecules (water and ions) around the r(G=C), r(A-U), d(G=C), and d(A-T) pairs, respectively. (iv) From a dynamic point of view, the stability of the hydration shell is the strongest for the r(G=C) pairs and the weakest for the d(A-T) pairs. (v) For RNA, the observed long-lived hydration patterns are essentially non-sequence dependent and involve water bridges located in the deep groove and linking OR atoms of adjacent phosphate groups. Maximum lifetimes are close to 400 ps. (vi) In contrast, for DNA, long-lived hydration patterns are sequence dependent and located in the minor groove. For d(CpG)12, water bridges linking the (G)N3 and (C)O2 with the O4' atoms of adjacent nucleotides with 400 ps maximum lifetimes are characterized while no such bridges are observed for d(TpA)12. (vii) Potassium ions are observed to bind preferentially to deep/major groove atoms at RpY steps, essentially d(GpC), r(GpC), and r(ApU), by forming ion-bridges between electronegative atoms of adjacent base-pairs. On average, about half an ion is observed per base-pair. Positive ion-binding determinants are related to the proximity of two or more electronegative atoms. Negative binding determinants are associated with the electrostatic and steric hindrance due to the proximity of electropositive amino groups and neutral methyl groups. Potassium ions form only transient contacts with phosphate groups.     

Hydration study of homopolypeptides by (2)H NMR

Deuteron T(1) and T(2) was studied as a function of hydration in homopolyglycine (PG) and homopolyproline (PP). Water deuteron relaxation rates in PG conform to a hydration model involving two types of primary hydration sites where water is directly bonded to the polymer. Once these sites are filled, additional water only bonds to water molecules at the primary sites and in so doing affect their dynamics. PP exhibits an anomalous T(1) and T(2) hydration dependence which has been interpreted in terms of a cooperative water molecule-PP molecule helical conformational rearrangement which occurs once a certain hydration level is reached. The proposal of a water-PP structure is tested using molecular dynamics simulations.     

Proton polarizability of poly(L-tyrosine)–hydrogen phosphate hydrogen bonds as a function of alkali cations

We studied films of poly(L -tyrosine) with hydrogen phosphate (residue/phosphate, 1:1) by ir spectroscopy. The influences of the alkali cations (Li⁺, Na⁺, K⁺) and of the degree of hydration were clarified. If Li⁺ ions are present, the OH ?^?OP hydrogen bonds formed in the dried films between the tyrosine OH groups and hydrogen phosphate are asymmetrical. The formation of hydrogen phosphate–hydrogen phosphate hydrogen bonds is prevented by the presence of the Li⁺ ions. With an increase in the degree of hydration, the tyrosine–phosphate bonds are not broken but become slightly stronger. Completely different behaviour is found if K⁺ ions are present. In dry films, the OH ?^?OP ? O^? ?HOP hydrogen bonds formed between tyrosine and hydrogen phosphate show large proton polarizability. The tyrosine proton has a noticeable residence time at the acceptor O atom of the phosphate. The difference in the behaviour of the system with K⁺ ions when compared to the system with Li⁺ ions can be explained, since the hydrogen acceptor O atom of phosphate ions is more negatively charged due to the weaker influence of the K⁺ ions. Furthermore, POH ?^?OP hydrogen bonds between hydrogen phosphate molecules are formed. With an increase in the degree of hydration, the tyrosine–hydrogen phosphate hydrogen bonds are broken, all tyrosine protons are found at the tyrosine residues, and the -PO₃^? groupings are in a symmetrical environment, indicating that the K⁺ ions are removed from these groupings. If the degree of hydration increases further, hydrogen-bonded systems such as hydrogen phosphate–water–hydrogen phosphate are formed that show large proton polarizability due to collective proton motion. When Na⁺ ions are present, the OH ?^?OP ? O^? ?HOP hydrogen bonds formed in dry films still show proton polarizability, but the residence time of the tyrosine proton at the phosphate is very short.     

A molecular dynamics study on binding recognition between several 4,5 and 4,6‐linked aminoglycosides with A‐site RNA

A molecular dynamics (MD) simulation has been performed for two sets of aminoglycoside antibiotics bound with an RNA duplex corresponding to the aminoacyl‐tRNA decoding site of the 16S rRNA to characterize the energetics and dynamics of binding for several aminoglycosides. The binding free energy, essential dynamics and hydration analysis have been conducted to characterize the dynamics' properties associated with the binding recognition between each set of antibiotics and the RNA duplex. We have built several dynamic models with reasonable binding free energies showing good correlation with the experimental data. We have also conducted a hydration analysis on some long residency water molecules detected as W8 and W49 sites around the U1406 · U1495 pair and which are found to be important in binding recognition and in causing some apparent stretch variations of this pair during the dynamic studies. In addition, we also find that the hydration sites with long residence time identified between the ring III of two 4,6‐linked antibiotics (tobramycin and kanamycin) and phosphate oxygen atoms of G1405/U1406 may be worthy of further exploration for rational drug design. Copyright © 2009 John Wiley & Sons, Ltd.     

Study of the structure of N-acyldipalmitoylphosphatidylethanolamines in aqueous dispersion by infrared and Raman spectroscopies

The effect of the headgroup chain length on the structure and on the thermotropic behavior of N-acyldipalmitoylphosphatidylethanolamines (N-acyl-DPPEs) has been studied by infrared and Raman spectroscopies. The results show that the N-acyl-DPPEs can be divided in two classes depending on the N-acyl chain length. When the N-acyl chain contains 10 carbon atoms or more, it penetrates into the bilayer while it remains at the level of the glycerol backbone for shorter N-acyl chains. For both classes of N-acyl-DPPEs, the rotation of the lipid chains in the liquid-crystalline phase is hindered by the presence of the N-acyl group. In addition, the disruption of the hydrogen bonds between the amino and phosphate groups by N-acylation of the amino group results in an increase of the hydration of the phosphate group compared to that in DPPE. The hydration occurred at both the phosphate and amide group levels; the phosphate group is more hydrated for phospholipids with long N-acyl chains while in the case of short-chain derivatives both the phosphate and amide groups are hydrated. This higher degree of hydration coupled with the immobilization of the lipid molecule may contribute to the bilayer stabilizer role of N-acyl-PEs since hydration is an important factor in bilayer stability.     

Thermodynamic analysis of water molecules at the surface of proteins and applications to binding site prediction and characterization

Water plays an essential role in determining the structure and function of all biological systems. Recent methodological advances allow for an accurate and efficient estimation of the thermodynamic properties of water molecules at the surface of proteins. In this work, we characterize these thermodynamic properties and relate them to various structural and functional characteristics of the protein. We find that high-energy hydration sites often exist near protein motifs typically characterized as hydrophilic, such as backbone amide groups. We also find that waters around alpha helices and beta sheets tend to be less stable than waters around loops. Furthermore, we find no significant correlation between the hydration site-free energy and the solvent accessible surface area of the site. In addition, we find that the distribution of high-energy hydration sites on the protein surface can be used to identify the location of binding sites and that binding sites of druggable targets tend to have a greater density of thermodynamically unstable hydration sites. Using this information, we characterize the FKBP12 protein and show good agreement between fragment screening hit rates from NMR spectroscopy and hydration site energetics. Finally, we show that water molecules observed in crystal structures are less stable on average than bulk water as a consequence of the high degree of spatial localization, thereby resulting in a significant loss in entropy. These findings should help to better understand the characteristics of waters at the surface of proteins and are expected to lead to insights that can guide structure-based drug design efforts.     

Hydration behaviour of POPC/C(12)-Bet mixtures investigated by sorption gravimetry, (31)P NMR spectroscopy and X-ray diffraction

The hydration behaviour of mixtures of the zwitterionic phospholipid 1-palmitoyl-2-oleolyl-sn-glycero-3-phosphocholine (POPC) and the zwitterionic surfactant N,N-dimethyl-N-dodecyl-betain (C(12)-Bet) was investigated by sorption gravimetry, solid-state (31)P NMR-spectroscopy and small angle X-ray diffraction (SAXD). Negative excess hydration (dehydration) was found for almost all hydration degrees investigated. This behaviour is explained by the formation of an inner salt between the dipoles of phospholipid and surfactant headgroups that show a reverse sequence of partial charges with respect to the hydrocarbon backbone. The formation of an inner-salt most probably reduces potential water binding sites. Moreover, NMR data suggest that the incorporation of the zwitterionic surfactant into the phospholipid membrane is correlated with reorientation of the phosphate axis towards the membrane director as well as with reduced lateral and wobbling diffusion.     

Effect of calcium and magnesium on phosphatidylserine membranes: experiments and all-atomic simulations

It is known that phosphatidylserine (PS(-)) lipids have a very similar affinity for Ca(2+) and Mg(2+) cations, as revealed by electrokinetic and stability experiments. However, despite this similar affinity, experimental evidence shows that the presence of Ca(2+) or Mg(2+) induces very different aggregation behavior for PS(-) liposomes as characterized by their fractal dimensions. Also, turbidity measurements confirm substantial differences in aggregation behavior depending on the presence of Ca(2+) or Mg(2+) cations. These puzzling results suggest that although these two cations have a similar affinity for PS(-) lipids, they induce substantial structural differences in lipid bilayers containing each of these cations. In other words, these cations have strong ion-specific effects on the structure of PS(-) membranes. This interpretation is supported by all-atomic molecular-dynamics simulations showing that Ca(2+) and Mg(2+) cations have different binding sites and induce different membrane hydration. We show that although both ions are incorporated deep into the hydrophilic region of the membrane, they have different positions and configurations at the membrane. Absorbed Ca(2+) cations present a peak at a distance ~2 nm from the center of the lipid bilayer, and their most probable binding configuration involves two oxygen atoms from each of the charged moieties of the PS molecule (phosphate and carboxyl groups). In contrast, the distribution of absorbed Mg(2+) cations has two different peaks, located a few angstroms before and after the Ca(2+) peak. The most probable configurations (corresponding to these two peaks) involve binding to two oxygen atoms from carboxyl groups (the most superficial binding peak) or two oxygen atoms from phosphate groups (the most internal peak). Moreover, simulations also show differences in the hydration structure of the membrane: we obtained a hydration of 7.5 and 9 water molecules per lipid in simulations with Ca(2+) and Mg(2+), respectively.     

Solvent participation in Serratia marcescens endonuclease complexes

The monomer and dimer of the bacterium Serratia marcescens endonuclease (SMnase) are each catalytically active and the two subunits of the dimer function independently of each other. Specific interfacial waters may play a role in stability, complex formation, and functionality. We performed molecular dynamics simulations of both a subunit of SMnase and its model built complex with DNA and analyzed the relation of the hydration sites to the catalytic mechanism. It was found that the binding of DNA has little influence on the global hydration properties of the protein, including occupancy and water residence time distributions. DNA and protein recognition in our model mainly involves direct contacts by hydrogen bond and hydrophobic interactions. Water-mediated contacts exist, but are less common. Three interior water clusters were identified for SMnase. One cluster around the active site in the monomer-DNA complex shows relatively strong interactions between hydration sites as well as between the sites and the biomolecules. The simulated cluster properties agreed well with experimental data. The magnesium ion shows ligand exchange. Although Mg2+ keeps six ligands during the entire simulation, upon the binding of DNA, Asn119 loses its coordination with Mg2+, while one nonbridging oxygen of the phosphate of a DNA residue and two oxygen atoms of the side chain of Glu127 become the ligands. Waters in a nearby cluster exchange and participate in the resolvation of groups in the presence of DNA. Water thus not only participates in the cleavage of DNA but also can stabilize the transition state and the leaving groups in our model.     

Anion binding to nucleic acids

Nucleic acids are generally considered as efficient cation binders. Therefore, the likelihood that negatively charged ions might intrude their first hydration shell is rarely considered. Here, we show on the basis of (i) a survey of the Nucleic Acid Database, (ii) several structures extracted from the Cambridge Structural Database, and (iii) molecular dynamics simulations, that the nucleotide electropositive edges involving mainly amino, imino, and hydroxyl groups can cast specific anion binding sites. These binding sites constitute also good locations for the binding of the negatively charged groups of the Asp and Glu residues or the nucleic acid phosphate groups. Furthermore, it is observed in several instances that anions, like water molecules and cations, do mediate protein/nucleic acid interactions. Thus, anions as well as negatively charged groups are directly involved in specific recognition and folding phenomena involving polyanionic nucleic acids.     

A Connected-cluster of hydration around myoglobin: Correlation between molecular dynamics simulations and experiment

An analysis of a molecular dynamics simulation of metmyoglobin in an explicit solvent environment of 3,128 water molecules has been performed. Both statics and dynamics of the protein-solvent interface are addressed in a comparison with experiment. Three-dimensional density distributions, temperature factors, and occupancy weights are computed for the solvent by using the trajectory coordinates. Analysis of the hydration leads to the localization of more than 500 hydration sites distributed into multiple layers of solvation located between 2.6 and 6.8 Å from the atomic protein surface. After locating the local solvent density maxima or hydration sites we conclude that water molecules of hydration positions and hydration sites are distinct concepts. Both global and detailed properties of the hydration cluster around myoglobin are compared with recent neutron and X-ray data on myoglobin. Questions arising from differences between X-ray and neutron data concerning the locations of the protein-bound water are investigated. Analysis of water site differences found from X-ray and neutron experiments compared with our simulation shows that the simulation gives a way to unify the hydration picture given by the two experiments. © 1994 John Wiley & Sons, Inc.     

An analysis of the relationship between hydration and protein-DNA interactions.

Eleven protein-DNA crystal structures were analyzed to test the hypothesis that hydration sites predicted in the first hydration shell of DNA mark the positions where protein residues hydrogen-bond to DNA. For nine of those structures, protein atoms, which form hydrogen bonds to DNA bases, were found within 1.5 A of the predicted hydration positions in 86% of the interactions. The correspondence of the predicted hydration sites with the hydrogen-bonded protein side chains was significantly higher for bases inside the conserved DNA recognition sequences than outside those regions. In two CAP-DNA complexes, predicted base hydration sites correctly marked 71% (within 1.5 A) of protein atoms, which form hydrogen bonds to DNA bases. Phosphate hydration was compared to actual protein binding sites in one CAP-DNA complex with 78% marked contacts within 2.0 A. These data suggest that hydration sites mark the binding sites at protein-DNA interfaces.