Networks of Water Molecules in a Proflavine Deoxydinucleoside Phosphate Complex

Abstract

Using previously reported ab-initio atom-atom potentials for the interactions of a water molecule with phosphates, sugars and bases and newly computed ab-initio atom-atom potentials for the interaction between a proflavine ion and water, we have analyzed with the Monte-Carlo Metropolis method networks of water molecules hydrating a 2:2 complex of proflavine and deoxycytidylyl-3′,5′-guanosine, recently studied with X-ray crystallography. From our simulations we have i) verified the quality of our atom-atom potentials by obtaining patterns of oxygen atoms in very good agreement with the X-ray patterns for the minor groove and in reasonable agreement in the major groove, ii) predicted the water's hydrogen atoms positions and iii) preliminarily predicted the number of water molecules not reported in the X-ray study but present in the major groove. The above data, even if preliminary, and the analyses on the energetics of the water-water, water-proflavine and water-dCpG interactions indicate that very detailed accounts on the water filaments in the above crystal can be obtained optimally by merging computer and X-ray experiments.     

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.     

Hydration of transfer RNA molecules: a crystallographic study

Four crystal structures of transfer RNA molecules were refined at 3 A resolution with the inclusion of the solvent molecules found in the difference maps: yeast tRNA-phe in the orthorhombic form, yeast tRNA-phe in the monoclinic form and yeast tRNA-asp in the A and B forms. Over 100 solvent molecules were located in each tRNA crystal. Several hydration schemes are found repeatedly in the 4 crystals. The tertiary interactions in the corner of the L-shaped molecule attract numerous solvent molecules which bridge the ribose hydroxyl O(2') atoms, base exocyclic atoms and phosphate anionic oxygen atoms. Conservation of bases leads to conservative localized hydration patterns. Several solvent molecules are found stabilizing unusual base pairs like the G-U pairs and those involving the pseudouridine base. Water bridges between the O(2') and the exocyclic atom O2 of pyrimidines or the N3 atom of purines are common. Water bridges occur frequently between successive anionic oxygen atoms of each strand as well as between N7 or other exocyclic atoms of successive bases in the major groove. Magnesium ions or spermine molecules are found to bind in the major groove of tRNA helices without specific interactions.     

Hydration of nucleic acid bases studied using novel atom-atom potential functions

New simple atom-atom potential functions for simulating behavior of nucleic acids and their fragments in aqueous solutions are suggested. These functions contains terms which are inversely proportional to the first (electrostatics), sixth (or tenth for the atoms, forming hydrogen bonds) and twelfth (repulsion of all the atoms) powers of interatomic distance. For the refinement of the potential function parameters calculations of ice lattice energy, potential energy and configuration of small clusters consisting of water and nucleic acid base molecules as well as Monte Carlo simulation of liquid water were performed. Calculations using new potential functions give rise to more linear hydrogen bonds between water and base molecules than using other potentials. Sites of preferential hydration of five nucleic bases - uracil, thymine, cytosine, guanine and adenine as well as of 6,6,9-trimethyladenine were found. In the most energetically favourable sites water molecular interacts with two adjacent hydrophilic centres of the base. Studies of interaction of the bases with several water molecules showed that water-water interactions play an important role in the arrangement of the nearest to the base water molecules. Hydrophilic centres are connected by "bridges" formed by hydrogen bonded water molecules. The results obtained are consistent with crystallographic and mass-spectrometric data.     

Crystal and molecular structure of the ammonium salt of the dinucleoside monophosphate d(CpG)

The crystal and molecular structure of the ammonium salt of deoxycytidylyl-(3'-5')-deoxyguanosine has been determined from 0.85 A resolution single crystal X-ray diffraction data. The crystals obtained by acetone diffusion technique at -20 degrees C, are orthorhombic, P212121, a = 12.880(2), b = 17444(2) and c = 27.642(2) A. The structure was solved by high resolution Patterson and Fourier methods and refined to R = 0.136. There are two d(CpG) molecules in the asymmetric unit forming a mini left handed Z-DNA helix. This is in contrast to the earlier reported forms of d(CpG) where the molecules form self base paired duplexes. There are two ammonium ions in the asymmetric unit. The major groove NH+4 ion interacts with N7 of guanines through water bridges besides making H-bonded interactions directly with the phosphate oxygen atoms. A second NH+4 ion is found in the minor groove interacting directly with the phosphate oxygen atoms. Symmetry related molecules pack in such a way that the cytosine base stacks on cytosine and guanine base on guanine. Our structure demonstrates that alternating d(CpG) sequences have the ability to adopt the left handed Z-DNA structure even at the dimer level i.e., in a sequence which is only two base pairs long.     

Quantum-chemical and empirical calculations on phospholipids IV. Glycero-phosphorylethanolamine in a quasihexagonal planar lattice

The energetically most stable head group conformations of a racemic mixture of diacyl-glycero-phosphorylethanolamine in a planar quasihexagonal lattice were calculated using empirical 1-6-12 atom-atom potential functions for intra- and intermolecular interactions. The results demonstrate that the conformation of phospholipid head groups in bilayer systems is determined by intramolecular interactions as well as by intermolecular interactions with neighbouring phospholipid molecules and with solvent molecules. The most stable conformers are that with a φ2 = guache^? conformation of the phosphodiester group. All conformers with a φ2 = gauche⁽⁺⁾ or trans conformation have energies more than 15 kcal ☆ mol^?1 above that of the global minimum. The calculated torsional angles ?1 and φ1 are in very good agreement with the results of the X-ray diffraction analysis of 1,2-dilauroyl-DL-phosphatidylethanolamine (DLPE) acetic acid single crystals.     

Application of noninteracting probes in the protein structure space for the construction of statistical potentials for atom-atom interactions

Some effects hindering the construction of knowledge-based potentials of atom-atom interactions in problems of biopolymer modeling have been considered. It was shown that some limitations are overcome by considering the distribution of distances between noninteracting probes in the protein structure space. It was shown that knowledge-based potentials thus constructed can be effectively used to analyze the hydration of protein atoms. Using this approach, it is possible to predict the structure water location in a protein globule and recognize the correctly folded protein structure among decoys.     

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.     

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.     

Visualisation of extensive water ribbons and networks in a DNA minor-groove drug complex.

The crystal structure is reported of a complex between an ethyl derivative of the minor-groove drug furamidine and the dodecanucleotide duplex d(CGCGAATTCGCG)2, which has been refined to 1.85 A resolution and an R factor of 16.6% for data collected at -173 degreesC. An exceptionally large number (220) of water molecules have been located. The majority of these occur in the first coordination shell of solvation. There are extensive networks of connected waters, both in the major and minor grooves. In particular, there are 21 water molecules associated with the minor-groove drug, via hydrogen bonds from the four charged nitrogen atoms. One cluster of four waters is situated in the groove itself; the majority are on the outer edge of the groove, and serve to bridge between the outward-directed drug nitrogen atoms and backbone phosphate oxygen atoms. These bridges are both intra- and inter-strand, with the net effect that the outer edge of the drug molecule is covered by ribbons of water molecules.     

Testing water-mediated DNA recognition by the Hin recombinase

The Hin recombinase specifically recognizes its DNA-binding site by means of both major and minor groove interactions. A previous X-ray structure, together with new structures of the Hin DNA-binding domain bound to a recombination half-site that were solved as part of the present study, have revealed that two ordered water molecules are present within the major groove interface. In this report, we test the importance of these waters directly by X-ray crystal structure analysis of complexes with four mutant DNA sequences. These structures, combined with their Hin-binding properties, provide strong support for the critical importance of one of the intermediate waters. A lesser but demonstrable role is ascribed to the second water molecule. The mutant structures also illustrate the prominent roles of thymine methyls both in stabilizing intermediate waters and in interfering with water or amino acid side chain interactions with DNA.     

Stabilization of Z-DNA by demethylation of thymine bases: 1.3-A single-crystal structure of d(m5CGUAm5CG)

Methylation of cytosine bases at the C5 position has been known to stabilize Z-DNA. We had previously predicted from calculations of solvent-accessible surfaces that the methyl group at the same position of thymine has a destabilizing effect on Z-DNA. In the current studies, the sequence d(m5CGUAm5CG) has been crystallized and its structure solved as Z-DNA to 1.3-A resolution. A well-defined octahedral hexaaquomagnesium complex was observed to bridge the O4 oxygens of the adjacent uridine bases at the major groove surface, and four well-structured water molecules were found in the minor groove crevice at the d(UA) dinucleotide. These solvent interactions were not observed in the previously published Z-DNA structure of the analogous d(m5CGTAm5CG) sequence. A comparison of the thymine and uridine structures supports our prediction that demethylation of thymine bases helps to stabilize Z-DNA. A comparison of this d(UA)-containing Z-DNA structure with the analogous d(TA) structure shows that access of the O4 position is hindered by the C5 methyl of thymine due to steric and hydrophobic inhibition. In the absence of the methyl group, a magnesium-water complex binds to and slightly affects the structure of the Z-DNA major groove surface. This perturbation of the solvent structure at the major groove surface is translated into a much larger 1.41-A widening of the minor groove crevice, thereby allowing the specific binding of two water molecules at well-defined sites of each internal d(UA) base pair. Possible mechanisms by which modifications at the major groove surface of Z-DNA can affect the solvent properties of the minor groove crevice are discussed.     

Molecular dynamics simulations of dinucleoside and dinucleoside-drug crystal hydrates.

Molecular dynamics simulations have been performed on the dinucleoside monophosphates rGpC and dCpG, the latter in its intercalation complex with the acridine drug proflavine. The simulations were performed on the crystal structures, with crystallographically-located solvent molecules. It was found that satisfactory results were best obtained with restraints placed on the movements of the water molecules. Motions of individual atoms have been examined in terms of rms fluctuations and anisotropy and correlation functions. Relative motions of groups (phosphates, sugars, bases and proflavine molecules) have been analysed.     

Water structure around a left-handed Z-DNA fragment analyzed by cryo neutron crystallography

Even in high-quality X-ray crystal structures of oligonucleotides determined at a resolution of 1 Å or higher, the orientations of first-shell water molecules remain unclear. We used cryo neutron crystallography to gain insight into the H-bonding patterns of water molecules around the left-handed Z-DNA duplex [d(CGCGCG)]₂. The neutron density visualized at 1.5 Å resolution for the first time allows us to pinpoint the orientations of most of the water molecules directly contacting the DNA and of many second-shell waters. In particular, H-bond acceptor and donor patterns for water participating in prominent hydration motifs inside the minor groove, on the convex surface or bridging nucleobase and phosphate oxygen atoms are finally revealed. Several water molecules display entirely unexpected orientations. For example, a water molecule located at H-bonding distance from O6 keto oxygen atoms of two adjacent guanines directs both its deuterium atoms away from the keto groups. Exocyclic amino groups of guanine (N2) and cytosine (N4) unexpectedly stabilize waters H-bonded to O2 keto oxygens from adjacent cytosines and O6 keto oxygens from adjacent guanines, respectively. Our structure offers the most detailed view to date of DNA solvation in the solid-state undistorted by metal ions or polyamines.     

Structural insights into the effect of hydration and ions on A-tract DNA: a molecular dynamics study

DNA structure is known to be sensitive to hydration and ionic environment. To explore the dynamics, hydration, and ion binding features of A-tract sequences, a 7-ns Molecular dynamics (MD) study has been performed on the dodecamer d(CGCAAATTTGCG)(2). The results suggest that the intrusion of Na(+) ion into the minor groove is a rare event and the structure of this dodecamer is not very sensitive to the location of the sodium ions. The prolonged MD simulation successfully leads to the formation of sequence dependent hydration patterns in the minor groove, often called spine of hydration near the A-rich region and ribbon of hydration near the GC regions. Such sequence dependent differences in the hydration patterns have been seen earlier in the high resolution crystal structure of the Drew-Dickerson sequence, but not reported for the medium resolution structures (2.0 approximately 3.0 A). Several water molecules are also seen in the major groove of the MD simulated structure, though they are not highly ordered over the extended MD. The characteristic narrowing of the minor groove in the A-tract region is seen to precede the formation of the spine of hydration. Finally, the occurrence of cross-strand C2-H2.O2 hydrogen bonds in the minor groove of A-tract sequences is confirmed. These are found to occur even before the narrowing of the minor groove, indicating that such interactions are an intrinsic feature of A-tract sequences.     

Hydration of Nucleic Acid Bases Studied Using Novel Atom-Atom Potential Functions

Abstract

New simple atom-atom potential functions for simulating behavior of nucleic acids and their fragments in aqueous solutions are suggested. These functions contain terms which are inversely proportional to the first (electrostatics), sixth (or tenth for the atoms, forming hydrogen bonds) and twelfth (repulsion of all the atoms) powers of interatomic distance. For the refinement of the potential function parameters calculations of ice lattice energy, potential energy and configuration of small clusters consisting of water and nucleic acid base molecules as well as Monte Carlo simulation of liquid water were performed. Calculations using new potential functions give rise to more linear hydrogen bonds between water and base molecules than using other potentials. Sites of preferential hydration of five nucleic bases—uracil, thymine, cytosine, guanine and adenine as well as of 6,6,9-trimethyladenine were found. In the most energetically favourable sites water molecule interacts with two adjacent hydrophilic centres of the base. Studies of interaction of the bases with several water molecules showed that water-water interaction play an important role in the arrangement of the nearest to the base water molecules. Hydrophilic centres are connected by “bridges” formed by hydrogen bonded water molecules. The results obtained are consistent with crystallographic and mass-spectrometric data.     

A Monte Carlo simulation study of the aqueous hydration of r(GpC)2: comparison with crystallographic ordered water sites

Monte Carlo computer simulation is described for the dinucleotide duplex rGpC together with 562 water molecules at an environmental density of 1 g/cc in a cubic cell under periodic boundary conditions. Water-water interactions were treated using the TIP4P potential and the solute water interactions by TIP4P spliced with the nonbonded interactions from the AMBER 3.0 force field. The simulation was subjected to proximity analysis to obtain solute coordinate numbers and pair interaction energies for each solute atom. Hydration density distributions partitioned into contributions from the major groove side, the minor groove side and the sugar-phosphate backbone were examined, and the probabilities of occurrence for one- and two-water bridges in the simulation were enumerated. The results were compared with observations of crystallographic ordered water sites from x-ray diffraction studies on G and C containing small molecules, and in crystal structure determinations of the sodium, calcium, and ammonium salts of rGpC. The calculated results are generally consistent with the observed sites, except for cytosine N4, where a hydration site is predicted yet none observed in rGpC salts, and for guanine N3, which appears in this calculation to compete unfavorably with the adjacent donor site at guanine N2. There is, however, a significant probability of finding a one-water G-N3-W-G-N2 bridge indicated in the simulation. An explanation for the guanine N3 discrepancy in terms of electrostatic potentials is also offered. The calculated one- and two-water bridges in the rGpC hydration complex coincide in a number of cases to those observed in the ordered water structure of the sodium rGpC crystal hydrate.     

Atomic-resolution STM structure of DNA and localization of the retinoic acid binding site

Single-molecule imaging by scanning tunnelling microscopy (STM) yields the atomic-resolution (0.6A) structure of individual B-type DNA molecules. The strong correlation between these STM structures and those predicted from the known base sequence indicates that sequencing of single DNA molecules using STM may be feasible. There is excellent agreement between the STM and X-ray structures, but subtle differences exist due to radial distortions. We show that the interactions of other molecules with DNA, their binding configurations, and the structure of these complexes can be studied at the single-molecule level. The anti-cancer drug retinoic acid (RA) binds selectively to the minor groove of DNA with up to 6 RA molecules per DNA turn and with the plane of the RA molecule approximately parallel to the DNA symmetry axis. Similar studies for other drug molecules will be valuable in the a priori evaluation of the effectiveness of anti-cancer drugs.     

1 A crystal structures of B-DNA reveal sequence-specific binding and groove-specific bending of DNA by magnesium and calcium

The 1 A resolution X-ray crystal structures of Mg(2+) and Ca(2+) salts of the B-DNA decamers CCAACGTTGG and CCAGCGCTGG reveal sequence-specific binding of Mg(2+) and Ca(2+) to the major and minor grooves of DNA, as well as non-specific binding to backbone phosphate oxygen atoms. Minor groove binding involves H-bond interactions between cross-strand DNA base atoms of adjacent base-pairs and the cations' water ligands. In the major groove the cations' water ligands can interact through H-bonds with O and N atoms from either one base or adjacent bases, and in addition the softer Ca(2+) can form polar covalent bonds bridging adjacent N7 and O6 atoms at GG bases. For reasons outlined earlier, localized monovalent cations are neither expected nor found.Ultra-high atomic resolution gives an unprecedented view of hydration in both grooves of DNA, permits an analysis of individual anisotropic displacement parameters, and reveals up to 22 divalent cations per DNA duplex. Each DNA helix is quite anisotropic, and alternate conformations, with motion in the direction of opening and closing the minor groove, are observed for the sugar-phosphate backbone. Taking into consideration the variability of experimental parameters and crystal packing environments among these four helices, and 24 other Mg(2+) and Ca(2+) bound B-DNA structures, we conclude that sequence-specific and strand-specific binding of Mg(2+) and Ca(2+) to the major groove causes DNA bending by base-roll compression towards the major groove, while sequence-specific binding of Mg(2+) and Ca(2+) in the minor groove has a negligible effect on helix curvature. The minor groove opens and closes to accommodate Mg(2+) and Ca(2+) without the necessity for significant bending of the overall helix.The program Shelxdna was written to facilitate refinement and analysis of X-ray crystal structures by Shelxl-97 and to plot and analyze one or more Curves and Freehelix output files.