Structure, dynamics and hydration of the nogalamycin-d(ATGCAT)2Complex determined by NMR and molecular dynamics simulations in solution.

The structure of the 1:1 nogalamycin:d(ATGCAT)2 complex has been determined in solution from high-resolution NMR data and restrained molecular dynamics (rMD) simulations using an explicit solvation model. The antibiotic intercalates at the 5'-TpG step with the nogalose lying along the minor groove towards the centre of the duplex. Many drug-DNA nuclear Overhauser enhancements (NOEs) in the minor groove are indicative of hydrophobic interactions over the TGCA sequence. Steric occlusion prevents a second nogalamycin molecule from binding at the symmetry-related 5'-CpA site, leading to the conclusion that the observed binding orientation in this complex is the preferred orientation free of the complication of end-effects (drug molecules occupy terminal intercalation sites in all X-ray structures) or steric interactions between drug molecules (other NMR structures have two drug molecules bound in close proximity), as previously suggested. Fluctuations in key structural parameters such as rise, helical twist, slide, shift, buckle and sugar pucker have been examined from an analysis of the final 500 ps of a 1 ns rMD simulation, and reveal that many sequence-dependent structural features previously identified by comparison of different X-ray structures lie within the range of dynamic fluctuations observed in the MD simulations. Water density calculations on MD simulation data reveal a time-averaged pattern of hydration in both the major and minor groove, in good agreement with the extensive hydration observed in two related X-ray structures in which nogalamycin is bound at terminal 5'-TpG sites. However, the pattern of hydration determined from the sign and magnitude of NOE and ROE cross-peaks to water identified in 2D NOESY and ROESY experiments identifies only a few "bound" water molecules with long residence times. These solvate the charged bicycloaminoglucose sugar ring, suggesting an important role for water molecules in mediating drug-DNA electrostatic interactions within the major groove. The high density of water molecules found in the minor groove in X-ray structures and MD simulations is found to be associated with only weakly bound solvent in solution.     

NMR observation of individual molecules of hydration water bound to DNA duplexes: direct evidence for a spine of hydration water present in aqueous solution.

The residence times of individual hydration water molecules in the major and minor grooves of DNA were measured by nuclear magnetic resonance (NMR) spectroscopy in aqueous solutions of d-(CGCGAATTCGCG)2 and d-(AAAAATTTTT)2. The experimental observations were nuclear Overhauser effects (NOE) between water protons and the protons of the DNA. The positive sign of NOEs with the thymine methyl groups shows that the residence times of the hydration water molecules near these protons in the major groove of the DNA must be shorter than about 500 ps, which coincides with the behavior of surface hydration water in peptides and proteins. Negative NOEs were observed with the hydrogen atoms in position 2 of adenine in both duplexes studied. This indicates that a 'spine of hydration' in the minor groove, as observed by X-ray diffraction in DNA crystals, is present also in solution, with residence times significantly longer than 1 ns. Such residence times are reminiscent of 'interior' hydration water molecules in globular proteins, which are an integral part of the molecular architecture both in solution and in crystals.     

Water and ion binding around RNA and DNA (C,G) oligomers

The dynamics, hydration, and ion-binding features of two duplexes, the A(r(CG)(12)) and the B(d(CG)(12)), in a neutralizing aqueous environment with 0.25 M added KCl have been investigated by molecular dynamics (MD) simulations. The regular repeats of the same C=G base-pair motif have been exploited as a statistical alternative to long MD simulations in order to extend the sampling of the conformational space. The trajectories demonstrate the larger flexibility of DNA compared to RNA helices. This flexibility results in less well defined hydration patterns around the DNA than around the RNA backbone atoms. Yet, 22 hydration sites are clearly characterized around both nucleic acid structures. With additional results from MD simulations, the following hydration scale for C=G pairs can be deduced: A-DNA相似文献   

The role of flexibility and hydration on the sequence-specific DNA recognition by the Tn916 integrase protein: a molecular dynamics analysis

The N-terminal domain of the Tn916 integrase protein (INT-DBD) is responsible for DNA binding in the process of strand cleavage and joining reactions required for transposition of the Tn916 conjugative transposon. Site-specific association is facilitated by numerous protein-DNA contacts from the face of a three-stranded beta-sheet inserted into the major groove. The protein undergoes a subtle conformational transition and is slightly unfolded in the protein-DNA complex. The conformation of many charged residues is poorly defined by NMR data but mutational studies have indicated that removal of polar side chains decreases binding affinity, while non-polar contacts are malleable. Based on analysis of the binding enthalpy and binding heat capacity, we have reasoned that dehydration of the protein-DNA interface is incomplete. This study presents results from a molecular dynamics investigation of the INT-DBD-DNA complex aimed at a more detailed understanding of the role of conformational dynamics and hydration in site-specific binding. Comparison of simulations (total of 13 ns) of the free protein and of the bound protein conformation (in isolation or DNA-bound) reveals intrinsic flexibility in certain parts of the molecule. Conformational adaptation linked to partial unfolding appears to be induced by protein-DNA contacts. The protein-DNA hydrogen-bonding network is highly dynamic. The simulation identifies protein-DNA interactions that are poorly resolved or only surmised from the NMR ensemble. Single water molecules and water clusters dynamically optimize the complementarity of polar interactions at the 'wet' protein-DNA interface. The simulation results are useful to establish a qualitative link between experimental data on individual residue's contribution to binding affinity and thermodynamic properties of INT-DBD alone and in complex with DNA.     

Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.     

A unified picture of protein hydration: prediction of hydrodynamic properties from known structures.

Hydration is essential for the structural and functional integrity of globular proteins. How much hydration water is required for that integrity? A number of techniques such as X-ray diffraction, nuclear magnetic resonance (NMR) spectroscopy, calorimetry, infrared spectroscopy, and molecular dynamics (MD) simulations indicate that the hydration level is 0.3-0.5 g of water per gram of protein for medium sized proteins. Hydrodynamic properties, when accounted for by modeling proteins as ellipsoids, appear to give a wide range of hydration levels. In this paper we describe an alternative numerical technique for hydrodynamic calculations that takes account of the detailed protein structures. This is made possible by relating hydrodynamic properties (translational and rotational diffusion constants and intrinsic viscosity) to electrostatic properties (capacitance and polarizability). We show that the use of detailed protein structures in predicting hydrodynamic properties leads to hydration levels in agreement with other techniques. A unified picture of protein hydration emerges. There are preferred hydration sites around a protein surface. These sites are occupied nearly all the time, but by different water molecules at different times. Thus, though a given water molecule may have a very short residence time (approximately 100-500 ps from NMR spectroscopy and MD simulations) in a particular site, the site appears fully occupied in experiments in which time-averaged properties are measured.     

Dynamics of internal water in fatty acid binding protein: computer simulations and comparison with experiments

Multiple molecular dynamics (MD) simulations of fully solvated rat intestinal fatty acid binding protein (I-FABP) were conducted to investigate the dynamics of internal water molecules. Although the long time average of the number of internal water molecules in I-FABP is 22 as shown by the X-ray crystal structure, MD simulations predict large variations in the instantaneous number of internal water molecules on the nanosecond time scale. The computational model employed predicts that w135 (internal) and w217 (located on the protein surface) may be the water molecules with long residence times observed in previously reported magnetic relaxation dispersion studies. The average residence time of approximately 20 internal water molecules occupying the fatty acid binding cavity is estimated to be between 0.6 and 2.0 nanoseconds. Exchange of internal water in I-FABP appears to occur almost exclusively through the interface of beta-strands EF with the rest of the protein, which has significant implications for the pathways of the fatty acid entry and exit from the binding cavity. Proteins 2001;43:65-72.     

Rearrangement of side-chains in a Zif268 mutant highlights the complexities of zinc finger-DNA recognition.

Structural and biochemical studies of Cys(2)His(2) zinc finger proteins initially led several groups to propose a "recognition code" involving a simple set of rules relating key amino acid residues in the zinc finger protein to bases in its DNA site. One recent study from our group, involving geometric analysis of protein-DNA interactions, has discussed limitations of this idea and has shown how the spatial relationship between the polypeptide backbone and the DNA helps to determine what contacts are possible at any given position in a protein-DNA complex. Here we report a study of a zinc finger variant that highlights yet another source of complexity inherent in protein-DNA recognition. In particular, we find that mutations can cause key side-chains to rearrange at the protein-DNA interface without fundamental changes in the spatial relationship between the polypeptide backbone and the DNA. This is clear from a simple analysis of the binding site preferences and co-crystal structures for the Asp20-->Ala point mutant of Zif268. This point mutation in finger one changes the specificity of the protein from GCG TGG GCG to GCG TGG GC(G/T), and we have solved crystal structures of the D20A mutant bound to both types of sites. The structure of the D20A mutant bound to the GCG site reveals that contacts from key residues in the recognition helix are coupled in complex ways. The structure of the complex with the GCT site also shows an important new water molecule at the protein-DNA interface. These side-chain/side-chain interactions, and resultant changes in hydration at the interface, affect binding specificity in ways that cannot be predicted either from a simple recognition code or from analysis of spatial relationships at the protein-DNA interface. Accurate computer modeling of protein-DNA interfaces remains a challenging problem and will require systematic strategies for modeling side-chain rearrangements and change in hydration.     

Interfacial water as a "hydration fingerprint" in the noncognate complex of BamHI

The molecular code of specific DNA recognition by proteins as a paradigm in molecular biology remains an unsolved puzzle primarily because of the subtle interplay between direct protein-DNA interaction and the indirect contribution from water and ions. Transformation of the nonspecific, low affinity complex to a specific, high affinity complex is accompanied by the release of interfacial water molecules. To provide insight into the conversion from the loose to the tight form, we characterized the structure and energetics of water at the protein-DNA interface of the BamHI complex with a noncognate sequence and in the specific complex. The fully hydrated models were produced with Grand Canonical Monte Carlo simulations. Proximity analysis shows that water distributions exhibit sequence dependent variations in both complexes and, in particular, in the noncognate complex they discriminate between the correct and the star site. Variations in water distributions control the number of water molecules released from a given sequence upon transformation from the loose to the tight complex as well as the local entropy contribution to the binding free energy. We propose that interfacial waters can serve as a "hydration fingerprint" of a given DNA sequence.     

Long-residency hydration, cation binding, and dynamics of loop E/helix IV rRNA-L25 protein complex

Molecular dynamics simulations of RNA-protein complex between Escherichia coli loop E/helix IV (LE/HeIV) rRNA and L25 protein reveal a qualitative agreement between the experimental and simulated structures. The major groove of LE is a prominent rRNA cation-binding site. Divalent cations rigidify the LE major groove geometry whereas in the absence of divalent cations LE extensively interacts with monovalent cations via inner-shell binding. The HeIV region shows bistability of its major groove explaining the observed differences between x-ray and NMR structures. In agreement with the experiments, the simulations suggest that helix-alpha1 of L25 is the least stable part of the protein. Inclusion of Mg2+ cations into the simulations causes perturbation of basepairing at the LE/HeIV junction, which does not, however, affect the protein binding. The rRNA-protein complex is mediated by a number of highly specific hydration sites with long-residing water molecules and two of them are bound throughout the entire 24-ns simulation. Long-residing water molecules are seen also outside the RNA-protein contact areas with water-binding times substantially enhanced compared to simulations of free RNA. Long-residency hydration sites thus represent important elements of the three-dimensional structure of rRNA.     

Protein hydration studied with homonuclear 3D1H NMR experiments

Summary Homonuclear 3D¹H NOESY-TOCSY and 3D¹H ROESY-TOCSY experiments were used to resolve and assign nuclear Overhauser effect (NOE) cross peaks between the water signal and individual polypeptide proton resonances in H₂O solutions of the basic pancreatic trypsin inhibitor. Combined with a novel, robust water-suppression technique, positive and negative intermolecular NOEs were detected at 4°C. The observation of positive NOEs between water protons and protein protons enables more precise estimates of the very short residence times of the water molecules in the hydration sites on the protein surface.     

Dynamics of water molecules buried in cavities of apolipoprotein E studied by molecular dynamics simulations and continuum electrostatic calculations

Molecular dynamics (MD) simulations of several nanoseconds each were used to monitor the dynamic behavior of the five crystal water molecules buried in the interior of the N-terminal domain of apolipoprotein E. These crystal water molecules are fairly well conserved in several apolipoprotein E structures, suggesting that they are not an artifact of the crystal and that they may have a structural and/or functional role for the protein. All five buried crystal water molecules leave the protein interior in the course of the longest simulations and exchange with water molecules from the bulk. The free energies of binding evaluated from the electrostatic binding free energy computed using a continuum model and estimates of the binding entropy changes represent shallow minima. The corresponding calculated residence times of the buried water molecules range from tens of picoseconds to hundreds of nanoseconds, which denote rather short times as for buried water molecules. Several water exchanges monitored in the simulations show that water molecules along the exit/entrance pathway use a relay of H bonds primarily formed with charged residues which helps either the exit or the entrance from or into the buried site. The exit/entrance of water molecules from/into the sites is permitted essentially by local motions of, at most, two side chains, indicating that, in these cases, complex correlated atomic motions are not needed to open the buried site toward the surface of the protein. This provides a possible explanation for the short residence times.     

Binding of the estrogen receptor to DNA. The role of waters.

Molecular dynamics simulations are carried out to investigate the binding of the estrogen receptor, a member of the nuclear hormone receptor family, to specific and non-specific DNA. Two systems have been simulated, each based on the crystallographic structure of a complex of a dimer of the estrogen receptor DNA binding domain with DNA. One structure includes the dimer and a consensus segment of DNA, ds(CCAGGTCACAGTGACCTGG); the other structure includes the dimer and a nonconsensus segment of DNA, ds(CCAGAACACAGTGACCTGG). The simulations involve an atomic model of the protein-DNA complex, counterions, and a sphere of explicit water with a radius of 45 A. The molecular dynamics package NAMD was used to obtain 100 ps of dynamics for each system with complete long-range electrostatic interactions. Analysis of the simulations revealed differences in the protein-DNA interactions for consensus and nonconsensus sequences, a bending and unwinding of the DNA, a slight rearrangement of several amino acid side chains, and inclusion of water molecules at the protein-DNA interface region. Our results indicate that binding specificity and stability is conferred by a network of direct and water mediated protein-DNA hydrogen bonds. For the consensus sequence, the network involves three water molecules, residues Glu-25, Lys-28, Lys-32, Arg-33, and bases of the DNA. The binding differs for the nonconsensus DNA sequence in which case the fluctuating network of hydrogen bonds allows water molecules to enter the protein-DNA interface. We conclude that water plays a role in furnishing DNA binding specificity to nuclear hormone receptors.     

A "structural" water molecule in the family of fatty acid binding proteins

A single water molecule (w135), buried within the structure of rat intestinal fatty acid binding protein (I-FABP), is investigated by NMR, molecular dynamics simulations, and analysis of known crystal structures. An ordered water molecule was found in structurally analogous position in 24 crystal structures of nine different members of the family of fatty acid binding proteins. There is a remarkable conservation of the local structure near the w135 binding site among different proteins from this family. NMR cross-relaxation measurements imply that w135 is present in the I-FABP:ANS (1-sulfonato-8-(1')anilinonaphthalene) complex in solution with the residence time of >300 ps. Mean-square positional fluctuations of w135 oxygen observed in MD simulations (0.18 and 0.13 A2) are comparable in magnitude to fluctuations exhibited by the backbone atoms and result from highly constrained binding pocket as revealed by Voronoi volumes (averages of 27.0 +/- 1.8 A3 and 24.7 +/- 2.2 A3 for the two simulations). Escape of w135 from its binding pocket was observed only in one MD simulation. The escape process was initiated by interactions with external water molecules and was accompanied by large deformations in beta-strands D and E. Immediately before the release, w135 assumed three distinct states that differ in hydrogen bonding topology and persisted for about 15 ps each. Computer simulations suggest that escape of w135 from the I-FABP matrix is primarily determined by conformational fluctuations of the protein backbone and interactions with external water molecules.     

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480?K) and MM-PBSA at 300 and 360?K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.     

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson–Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480?K) and MM-PBSA at 300 and 360?K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.