NMR studies of the Escherichia coli Trp repressor.trpRs operator complex

To understand the specificity of the Escherichia coli Trp repressor for its operators, we have begun to study complexes of the protein with alternative DNA sequences, using 1H-NMR spectroscopy. We report here the 1H-NMR chemical shifts of a 20-bp oligodeoxynucleotide containing the sequence of a symmetrised form of the trpR operator in the presence and absence of the holorepressor. Deuterated protein was used to assign the spectrum of the oligodeoxynucleotide in a 37-kDa complex with the Trp holorepressor. Many of the resonances of the DNA shift on binding to the protein, which suggests changes in conformation throughout the sequence. The largest changes in shifts for the aromatic protons in the major groove are for A15 and G16, which are thought to hydrogen bond to the protein, possibly via water molecules. We have also examined the effect of DNA binding on the corepressor, tryptophan, in this complex. The indole proton resonance of the tryptophan undergoes a downfield shift of 1.2 ppm upon binding of DNA. This large shift is consistent with hydrogen bonding of the tryptophan to the phosphate backbone of the trpR operator DNA, as in the crystal structure of the holoprotein with the trp operator.     

RNA hydration: a detailed look

The crystal structure of the RNA duplex [r(CCCCGGGG)]2 has been refined to 1.46 A resolution with room temperature synchrotron diffraction data. This represents the highest resolution reported to date for an all-RNA oligonucleotide and is well beyond the best resolution ever achieved with an A-form DNA duplex. The analysis of the ordered hydration around the octamer duplex reveals conserved regular arrangements of water molecules in both grooves. In the major groove, all located first shell water molecules can be fitted into a pattern that is repeated through all eight base pairs, involves half the phosphate oxygens, and joins the two strands. In the minor groove, roughly across its narrowest dimension, tandem water molecules link the 2'-hydroxyl groups of adjacent nucleotides in base-pair steps in a similarly regular fashion. The structure provides evidence for an important role of the 2'-hydroxyl groups in the thermodynamic stabilization of RNA, beyond their known functions of locking the sugar pucker and mediating 3' --> 5' intrastrand O2'...O4' hydrogen bonds. The ribose 2'-hydroxyls lay the foundation for the enthalpic stability of the RNA relative to the DNA duplex, both as a scaffold for the water network in the minor groove and through their extensive individual hydration.     

Determination of the residence time of water molecules hydrating B'- DNA and B-DNA, by one-dimensional zero-enhancement nuclear Overhauser effect spectroscopy

The residence time of water in the minor groove of the d(CGCGAATTCGCG) duplex has been determined by a recent measurement combining nuclear Overhauser enhancements (NOE, ROE) and 17O relaxation dispersion. The time is in the range of nanoseconds, so that it may be measured by a rather simple method proposed here, namely the choice of conditions such that the NOE between the observed DNA proton and a nearby water proton is zero. This condition is realized when the residence time of the water molecule is 0.178 times the nuclear magnetic resonance period (e.g. 0.297 ns at 600 MHz). It may be achieved by varying the magnetic field and/or the temperature. The zero-NOE measurement may be performed by one-dimensional NMR, and has therefore good sensitivity. We have developed excitation sequences which suppress two spurious contributions to the NOE: from neighboring exchangeable protons and from H3' protons whose chemical shift is close to that of water. The method is applied here to the comparison of residence times of water next to B-DNA and next to B'-DNA, the latter corresponding to better stacked, propeller-twisted base-pairs and a correspondingly narrower minor groove. In the minor groove of [d(CGCGAATTCGCG)]2, a B'-DNA duplex, the residence time of the water molecule next to H2 of adenine(6) (underlined), is 0.6 ns at 10 degreesC, in good agreement with the value obtained previously. The residence time is slightly but distinctly shorter for the water next to A5, suggesting non-cooperative departure of these two molecules which are presumed to be part of the hydration spine. Near A5 and A4 of [d(AAAAATTTTT)]2, another B'-DNA duplex, the residence times are approximately twice as long, but the activation enthalpies are about the same, ca. 38 kJ/mol. The residence time in the minor groove of the regular B-DNA sequence d(CGCGATCGCG) was 0.3 ns at 10 degreesC, shorter than in the case of the B'-DNA sequences by factors of 2 and 4, respectively. The temperature dependence is less, with an activation enthalpy of 27 kJ/mol. The major groove residence times are comparable for the three sequences, and a few times shorter than those of minor groove water. A value of 0.36 ns, or even more in case of rotation of water, is obtained around -8 degreesC. The most striking aspect of these results is the relatively small difference in the residence times of reputedly fast and slow-exchanging water molecules bound to DNA in biological conditions. This suggests that the spine of hydration is perhaps not a major stabilizer of the B'-DNA structure as compared with B-DNA.     

Improved targeting of an anti-epidermal growth factor receptor variant III monoclonal antibody in tumor xenografts after labeling using N-succinimidyl 5-iodo-3-pyridinecarboxylate

The crystal structure is reported of a complex between the dodecanucleotide sequence d(CGCGAATTCGCG)2and an analogue of the DNA binding drug Hoechst 33258, in which the piperazine ring has been replaced by an amidinium group and the phenol ring by a phenylamidinium group. The structure has been refined to an R factor of 19.5% at 2.2 A resolution. The drug is held in the minor groove by five strong hydrogen bonds, together with bridging water molecules at both ends. There are few other contacts with the floor of the groove, indicating a lack of isohelicity with the groove and suggesting (i) that the observed high DNA affinity of this drug is primarily due to the array of hydrogen bonds and (ii) that these more than compensate for its poor isohelicity.     

Synergistic effects in the melting of DNA hydration shell: melting of the minor groove hydration spine in poly(dA).poly(dT) and its effect on base pair stability

We propose that water of hydration in contact with the double helix can exist in several states. One state, found in the narrow groove of poly(dA).poly(dT), should be considered as frozen to the helix, i.e., an integral part of the double helix. We find that this enhanced helix greatly effects the stability of that helix against base separation melting. Most water surrounding the helix is, however, melted or disassociated with respect to being an integral part of helix and plays a much less significant role in stabilizing the helix dynamically, although these water molecules play an important role in stabilizing the helix conformation statically. We study the temperature dependence of the melting of the hydration spine and find that narrow groove nonbonded interactions are necessary to stabilize the spine above room temperature and to show the broad transition observed experimentally. This calculation requires that synergistic effects of nonbonded interactions between DNA and its hydration shell affect the state of water-base atom hydrogen bonds. The attraction of waters into narrow groove tends to retain waters in the groove and compress or strain these hydrogen bonds.     

A model of water structure inside the HLA-A2 peptide binding groove

Based on molecular dynamics simulations, it is proposed that water within the binding groove of the human MHC class I molecule HLA-A2 plays a role in the formation of its complex with the influenza matrix protein (residues 58-66; GILGFVFTL) peptide. In these simulations, a loosely structured network of water molecules is present in the binding groove between the peptide and the MHC molecule, and may be important in completing the peptide-MHC interface. In two independent 400 ps simulations where groove-based water molecules were included, the peptide remained essentially in the conformation observed in the crystal structure. In contrast, in a 400 ps simulation in which no water molecules were placed between the peptide and the MHC molecule, the crystal structure conformation was rapidly lost. The basis for this behavior appears to be that the groove-based water molecules help to maintain the appropriate orientation of the Arg-97 side chain of HLA-A2 and, in turn, the conformation of the central part of the peptide.     

Crystal structure of a B-DNA dodecamer containing inosine, d(CGCIAATTCGCG), at 2.4 A resolution and its comparison with other B-DNA dodecamers

The crystal structure of the dodecamer, d(CGCIAATTCGCG), has been determined at 2.4 A resolution by molecular replacement, and refined to an R-factor of 0.174. The structure is isomorphous with that of the B-DNA dodecamer, d(CGCGAATTCGCG), in space group P2(1)2(1)2(1) with cell dimensions of a = 24.9, b = 40.4, and c = 66.4 A. The initial difference Fourier maps clearly indicated the presence of inosine instead of guanine. The structure was refined with 44 water molecules, and compared to the parent dodecamer. Overall the two structures are very similar, and the I:C forms Watson-Crick base pairs with similar hydrogen bond geometry to the G:C base pairs. The propeller twist angle is low for I4:C21 and relatively high for the I16:C9 base pair (-3.2 degrees compared to -23.0 degrees), and the buckle angles alter, probably due to differences in the contacts with symmetry related molecules in the crystal lattice. The central base pairs of d(CGCIAATTCGCG) show the large propeller twist angles, and the narrow minor groove that characterize A-tract DNA, although I:C base pairs cannot form the major groove bifurcated hydrogen bonds that are possible for A:T base pairs.     

The use of diaminopurine to investigate structural properties of nucleic acids and molecular recognition between ligands and DNA

2,6-Diaminopurine (DAP) is an analogue of adenine which can be converted to nucleotides that serve as substrates for incorporation into nucleic acids by polymerases in place of (d)AMP. It pairs with thymidine (or uracil), engaging in three hydrogen bonds of the Watson-Crick type. The result of DAP incorporation is to add considerable stability to the double helix and to impart other structural features, such as an altered groove width and disruption of the normal spine of hydration. DNA containing DAP may or may not be recognized by restriction endonucleases; RNA containing DAP may not engage in normal splicing. The DAP.T pair affects the local flexibility of DNA and impedes the interaction with helix bending proteins. By providing a non-canonical hydrogen bond donor in the minor groove and/or blocking access to the floor of that groove it strongly affects interactions with small molecules such as antibiotics and anticancer drugs. Examples which illustrate altered recognition of nucleotide sequences in DAP-containing DNA are presented: changed sites of cutting by bleomycin, photocleavage by uranyl nitrate and footprinting with mithramycin. Using DNA in which both A-->DAP and G-->Inosine substitutions have been made it is possible to assess precisely the role of the purine 2-amino group in ligand-DNA recognition.     

On the importance of van der Waals interaction in the groove binding of DNA with ligands: restrained molecular dynamics study

Comparison of interaction energy between an oligonucleotide and a DNA-binding ligand in the minor and major groove modes was made by use of restrained molecular dynamics. Distortion in DNA was found for the major groove mode whereas less significant changes for both ligand and DNA were detected for the minor groove binding after molecular dynamics simulation. The conformation of the ligand obtained from the major groove modes resembles that computed with the ligand soaked in water. The van der Waals contact energy was found to be as significant as electrostatic energy and more important for difference in binding energy between these two binding modes. The importance of van der Waals force in groove binding was supported by computations on the complex formed by the repressor peptide fragment from the bacteriophage 434 and its operator oligonucleotide.     

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.     

Structure of poly d(AI).poly d(CT) in two different packing arrangements

X-ray diffraction analysis of poly d(AI).poly d(CT) in oriented and polycrystalline fibers has revealed the DNA structure to be a 10-fold, right-handed, antiparallel, Watson-Crick base paired double helix in two distinct packing arrangements corresponding to one and two helices, respectively, in the unit cell. The helix pitch is 32.1 A and 32.4 A in the two cases, almost 1.5 A shorter than in classical B-DNA. The resulting B'-DNA geometry, described in terms of a tetranucleotide repeat which is conformationally similar to B-DNA, has its minor groove closely shut and major groove correspondingly widened, thus striking a sharp morphological contrast to B-DNA. According to difference electron density maps, a spine of hydration along the minor groove connects both strands and provides structural stability; ordered sodium ions and water molecules are actively involved in bridging the phosphate groups of neighboring helices. The crystallographic R-values for these two allomorphs are 0.26 and 0.20, respectively, for data up to 3.0 A resolution.     

The high resolution crystal structure of the DNA decamer d(AGGCATGCCT)

The crystal structure of the DNA decamer d(AGGCATGCCT) has been determined to a resolution of 1.3 A and R factor of 13.9%. The structure has a unique conformation with each of the decamer single strands forming base-pairing interactions with two symmetry-related strands. The central eight bases of the decamer form an A-DNA octamer duplex with one symmetry-related strand whilst the terminal 5'-A and T-3' bases are flipped out and away from the octamer helix axis to form base-pairing interactions with a second symmetry-related strand. These A.T base-pairs lie perpendicular to the crystallographic c axis and pack within the unit cell in conjunction with a symmetry-related A.T base-pair displaced by 3.4 A degrees along the c axis. A novel base triplet interaction of the type A*(G.C) is present in the structure with interaction from the major groove side of the terminal 5'-A base to the minor groove of the central A-DNA octamer. This structure reports the first example of cobalt hexammine binding to a right-handed DNA duplex. The crystallographic asymmetric unit contains two cobalt hexammine ligands with one site in the major groove coordinating via hydrogen bonds to the 5'-AGG bases, and the second site located between DNA molecules and interacting with the oxygen atoms of phosphate groups.     

Base excision repair initiation revealed by crystal structures and binding kinetics of human uracil-DNA glycosylase with DNA

Three high-resolution crystal structures of DNA complexes with wild-type and mutant human uracil-DNA glycosylase (UDG), coupled kinetic characterizations and comparisons with the refined unbound UDG structure help resolve fundamental issues in the initiation of DNA base excision repair (BER): damage detection, nucleotide flipping versus extrahelical nucleotide capture, avoidance of apurinic/apyrimidinic (AP) site toxicity and coupling of damage-specific and damage-general BER steps. Structural and kinetic results suggest that UDG binds, kinks and compresses the DNA backbone with a 'Ser-Pro pinch' and scans the minor groove for damage. Concerted shifts in UDG simultaneously form the catalytically competent active site and induce further compression and kinking of the double-stranded DNA backbone only at uracil and AP sites, where these nucleotides can flip at the phosphate-sugar junction into a complementary specificity pocket. Unexpectedly, UDG binds to AP sites more tightly and more rapidly than to uracil-containing DNA, and thus may protect cells sterically from AP site toxicity. Furthermore, AP-endonuclease, which catalyzes the first damage-general step of BER, enhances UDG activity, most likely by inducing UDG release via shared minor groove contacts and flipped AP site binding. Thus, AP site binding may couple damage-specific and damage-general steps of BER without requiring direct protein-protein interactions.     

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.     

Molecular dynamics investigations of DNA triple helical models: unique features of the Watson-Crick duplex

We have built computer models of triple helical structures with a third poly(dT) strand Hoogsteen base paired to the major groove of a poly(dA).poly(dT) Watson-Crick (WC) base-paired duplex in the canonical A-DNA as well as B-DNA. For the A-DNA form, the sugar-phosphate backbone of the third strand intertwines and clashes with the poly(dA) strand requiring a radical alteration of the duplex to access the hydrogen bonding sites in the major groove. In contrast, when the duplex was in the canonical B-DNA form, the third strand was readily accommodated in the major groove without perturbing the duplex. The triple helical model, with the duplex in the B-DNA form, was equilibrated for 400ps using molecular dynamics simulations including water molecules and counter-ions. During the entire simulations, the deoxyriboses of the adenine strand oscillate between the S-type and E-type conformations. However, 30% of the sugars of the thymine strands-II & III switch to the N-type conformation early in the simulations but return to the S-type conformation after 200ps. In the equilibrium structure, the WC duplex portion of the triplex is unique and its geometry differs from both the A- or B-DNA. the deoxyriboses of the three strands predominantly exhibit S-type conformation. Besides the sugar pucker, the major groove width and the base-tilt are analogous to B-DNA, while the X-displacement and helical twist resemble A-DNA, giving a unique structure to the triplex and the Watson & Crick and Hoogsteen duplexes.     

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.     

Study of a noncovalent trp repressor: DNA operator complex by electrospray ionization time-of-flight mass spectrometry

Electrospray ionization time-of-flight mass spectrometry (ESI-TOF MS) has been used to study noncovalent interactions between the trp apo-repressor (TrpR), its co-repressor tryptophan and its specific operator DNA. In 5 mM ammonium acetate, TrpR was detected as a partially unfolded monomer. In the presence of a 21-base-pair DNA possessing the two symmetrically arranged CTAG consensus sequences required for specific TrpR binding, a homodimer-dsDNA complex with a 1:1 stoichiometry was observed. Co-repressor was not needed for the complex to form under our experimental conditions. Collision induced dissociation (CID-MS) revealed that this complex was very stable in the gas phase since dissociation was achieved only at energies that also broke covalent bonds. We saw no evidence for the presence of the six water molecules that mediate the interaction between the protein and the DNA in the crystal structure. To check the binding specificity of the TrpR for its target DNA, a competitive experiment was undertaken: the protein was mixed with an equimolar amount of three different DNAs in which the two CTAG sequences were separated by 2, 4, and 6 bp, respectively. Only the DNA with the correct consensus spacing of 4 bp was able to form stable interactions with TrpR. This experiment demonstrates the potential of ESI-MS to test the sequence-specificity of protein-DNA complexes. The interactions between the TrpR-DNA complex and 5-methyl-, L- and D-tryptophan were also investigated. Two molecules of 5-methyl- or L-tryptophan were bound with high affinity to the TrpR-DNA complex. On the other hand, D-tryptophan appeared to bind to the complex with poor specificity and poor affinity.     

Sequence-based diagnosis of hemoglobinopathies in the clinical laboratory

A high-angle neutron fibre diffraction study of the hydration of A-DNA has been performed using the single-crystal diffractometer D19 at the Institut Laue-Langevin (Grenoble, France). The sample was prepared using deuterated DNA extracted from E. Coli cells cultured on deuterated nutrients. In common with our previous neutron fibre diffraction studies of DNA, this work exploits the ability to isotopically replace H2O around the DNA by D2O. However this study benefitted additionally from the fact that the hydrogen atoms which are covalently bonded to carbon atoms in the DNA sugars and bases were replaced by deuterium so that incoherent scattering and absorption effects were minimised. Successive cycles of Fourier synthesis and Fourier difference synthesis allowed water peaks to be identified and their positional and occupancy parameters to be refined against the observed diffraction data. The results confirm the main hydration features noted in our earlier studies with a clear network of water running along the inside edge of the major groove linking successive OI phosphate oxygen atoms. The central core of water running along the axis of the double helix is very much clearer in this work. Additionally this study shows chains of ordered water lying in the centre of the major groove.     

Structural studies of binding site tryptophan mutants in the high-affinity streptavidin-biotin complex

Previous thermodynamic and computational studies have pointed to the important energetic role of aromatic contacts in generating the exceptional binding free energy of streptavidin-biotin association. We report here the crystallographic characterization of single site tryptophan mutants in investigating structural consequences of alterations in these aromatic contacts. Four tryptophan residues, Trp79, Trp92, Trp108 and Trp120, play an important role in the hydrophobic binding contributions, which along with a hydrogen bonding network and a flexible binding loop give rise to tight ligand binding (Ka approximately 10(13) M-1). The crystal structures of ligand-free and biotin-bound mutants, W79F, W108F, W120F and W120A, in the resolution range from 1.9 to 2.3 A were determined. Nine data sets for these four different mutants were collected, and structural models were refined to R-values ranging from 0.15 to 0.20. The major question addressed here is how these mutations influence the streptavidin binding site and in particular how they affect the binding mode of biotin in the complex. The overall folding of streptavidin was not significantly altered in any of the tryptophan mutants. With one exception, only minor deviations in the unbound structures were observed. In one crystal form of unbound W79F, there is a coupled shift in the side-chains of Phe29 and Tyr43 toward the mutation site, although in a different crystal form these shifts are not observed. In the bound structures, the orientation of biotin in the binding pocket was not significantly altered in the mutant complex. Compared with the wild-type streptavidin-biotin complex, there were no additional crystallographic water molecules observed for any of the mutants in the binding pocket. These structural studies thus suggest that the thermodynamic alterations can be attributed to the local alterations in binding residue composition, rather than a rearrangement of binding site architectures.