Search

The entry of HIV-1 into a target cell requires gp120 and receptor CD4 as well as coreceptor CCR5/CXCR4 recognition events associated with conformational changes of the involved proteins. The binding of CD4 to gp120 is the initiation step of the whole process involving structural rearrangements that are crucial for subsequent pathways. Despite the wealth of knowledge about the gp120/CD4 interactions, details of the conformational changes occurring at this stage remain elusive. We have performed molecular dynamics simulations in explicit solvent based on the gp120/CD4/CD4i crystal structure in conjunction with modeled V3 and V4 loops to gain insight into the dynamics of the binding process. Three differentiated interaction modes between CD4 and gp120 were found, which involve electrostatics, hydrogen bond and van der Waals networks. A "binding funnel" model is proposed based on the dynamical nature of the binding interface together with a CD4-attraction gradient centered in gp120 at the CD4-Phe43-binding cavity. Distinct dynamical behaviors of free and CD4-bound gp120 were monitored, which likely represent the ground and pre-fusogenic states, respectively. The transition between these states revealed concerted motions in gp120 leading to: i) loop contractions around the CD4-Phe43-insertion cavity; ii) stabilization of the four-stranded "bridging sheet" structure; and iii) translocation and clustering of the V3 loop and the bridging sheet leading to the formation of the coreceptor binding site. Our results provide new insight into the dynamic of the underlying molecular recognition mechanism that complements the biochemical and structural studies., (Copyright 2004 Wiley-Liss, Inc.)

Several studies have shown that biomolecular NMR structures are often of lower quality when compared to crystal structures, and consequently they are often excluded from structural analyses. We present a publicly available database of re-refined NMR structures, exhibiting significantly improved quality. This database (available at http://www.cmbi.kun.nl/dress/) presents a uniformly refined and validated set of structural models that improves the value of these NMR structures as input for experimental and theoretical studies in many fields of research., (Copyright 2004 Wiley-Liss, Inc.)

The protein CNOT4 possesses an N-terminal RING finger domain that acts as an E3 ubiquitin ligase and specifically interacts with UbcH5B, a ubiquitin-conjugating enzyme. The structure of the CNOT4 RING domain has been solved and the amino acids important for the binding to UbcH5B have been mapped. Here, the residues of UbcH5B important for the binding to CNOT4 RING domain were identified by NMR chemical shift perturbation experiments, and these data were used to generate structural models of the complex with the program HADDOCK. Together with the NMR data, additional biochemical data were included in a second docking, and comparisons of the resulting model with the structure of the c-Cbl/UbcH7 complex reveal some significant differences, notably at specific residues, and give structural insights into the E2/E3 specificity.

A structural model of the transient complex between the yeast copper chaperone Atx1 and the first soluble domain of the copper transporting ATPase Ccc2 was obtained with HADDOCK, combining NMR chemical shift mapping information with in silico docking. These two proteins are involved in copper trafficking in yeast cells. Calculations were performed starting with the copper ion either bound to Atx1 or to Ccc2 and using the experimental structures of the copper-loaded and apo forms of each protein. The copper binding motifs of the two proteins are found in close proximity. Copper tends to move from Atx1 to Ccc2, consistent with the physiological direction of transfer, with concomitant structural rearrangements, in agreement with experimental observations. The interaction is mainly of an electrostatic nature with hydrogen bonds stabilizing the complex. The structural data are relevant for a number of proteins homologous to Atx1 and Ccc2 and conserved from bacteria to humans.

The high sensitivity of modern NMR instrumentation, in combination with full deuteration, enabled the measurement of long-range NOEs between amide protons in a fully deuterated protein corresponding to distances up to 8A. These are beyond the limit normally observed in protonated samples. Such long-distance NOEs could be observed using long mixing times, which became possible due to reduced spin diffusion and T1 relaxation of the amide protons in the fully deuterated sample. This information was used in combination with secondary structure restraints derived from secondary chemical shifts for structure calculations. With these backbone amide proton NOEs only, a unique fold could be obtained with positional root mean square deviations from the average of 1.30 and 2.25 A for backbone and heavy atoms, respectively. Despite the low density of restraints, no mirror image problems were observed. Addition of sidechain NOE information increased the precision of the ensemble and in particular of the core packing. The structures obtained in this way were close to the published crystal structure. NOE completeness analysis revealed that the cumulative completeness is still more than 80% for an 8.0 A cut-off distance.

We present an in silico, structure-based approach for design and evaluation of conformationally restricted peptide-vaccines. In particular, we designed four cyclic peptides of ten or 11 residues mimicking the crystallographically observed beta-turn conformation of a predicted immunodominant loop of PorA from Neisseria meningitidis. Conformational correctness and stability of the peptide designs, as evaluated by molecular dynamics simulations, correctly predicted the immunogenicity of the peptides. We observed a peptide-induced functional antibody response that, remarkably, exceeded the response induced by the native protein in outer membrane vesicles, without losing specificity for related strains. The presented approach offers tools for a priori design and selection of peptide-vaccine candidates with full biological activity. This approach could be widely applicable: to outer membrane proteins of Gram-negative bacteria, and to other epitopes in a large range of pathogens.

The glucocorticoid receptor (GR) is a hormone dependent nuclear receptor that regulates gene transcription when bound to the glucocorticoid response element (GRE). The GRE acts as an allosteric effector, inducing a structural change in the GR DNA-binding domain (GR DBD) upon binding, thereby switching the GR to an active conformation. A similar conformational change can be induced by two single point mutations: Ser459Ala and Pro493Arg. Structural and dynamical aspects of the conformational switch have been investigated by molecular dynamics simulations. Our results indicate that these two mutants, which share a similar phenotype, exert their action at a structural level through different mechanisms. In the Arg(493) mutant, the D-loop and the second helix are stabilized in a conformation that preforms the protein-protein dimer interface. In the Ala(459) mutant, the structurally important hydrogen bond between Arg(496) and Ser(459) is missing, which leads to a core reorganization and a reorientation of the second helical region. Although remote, both in sequence and three dimensional structure, these two mutations induce structural changes that are ultimately reflected in similar regions of the GR DBD structure, namely the D-loop and the short second helical region. These correspond to hot area of the GR DBD that are important both for DNA-binding and for the proper formation of the protein-protein interface. The conformational rearrangements in these area are proposed to decrease unfavorable protein-DNA and protein-protein contacts and allow unspecific DNA-binding leading to the squelching phenotype of the mutants. The GR DBD can thus exist in two states, a transcriptionally active and an inactive state. Switching between these states can be accomplished either by GRE binding or by the described mutations., (Copyright 2003 Elsevier Science Ltd.)

Mersacidin belongs to the type B lantibiotics (lanthionine-containing antibiotics) that contain post-translationally modified amino acids and cyclic ring structures. It targets the cell wall precursor lipid II and thereby inhibits cell wall synthesis. In light of the emerging antibiotics resistance problem, the understanding of the antibacterial activity on a structural basis provides a key to circumvent this issue. Here we present solution NMR studies of mersacidin-lipid II interaction in dodecylphosphocholine (DPC) micelles. Distinct solution structures of mersacidin were determined in three different states: in water/methanol solution and in DPC micelles with and without lipid II. The structures in various sample conditions reveal remarkable conformational changes in which the junction between Ala-12 and Abu-13 (where Abu is aminobutyric acid) effectively serves as the hinge for the opening and closure of the ring structures. The DPC micelle-bound form resembles the previously determined NMR and x-ray crystal structures of mersacidin in pure methanol but substantially deviates from the other two states in our current report. The structural changes delineate the large chemical shift perturbations observed during the course of a two-step (15)N-(1)H heteronuclear single quantum coherence titration. They also modulate the surface charge distribution of mersacidin suggesting that electrostatics play a central role in the mersacidin-lipid II interaction. The observed conformational adaptability of mersacidin might be a general feature of lipid II-interacting antibiotics/peptides.

Biomolecular structures provide the basis for many studies in research areas such as structure-based drug design and homology modeling. In order to use molecular coordinates it is important that they are reliable in terms of accurate description of the experimental data and in terms of the overall and local geometry. Besides these primary quality criteria an indication is needed for the uncertainty in the atomic coordinates that may arise from the dynamic behavior of the considered molecules as well as from experimental- and computational procedures. In contrast to the crystallographic B-factor, a good measure for the uncertainty in NMR-derived atomic coordinates is still not available. It has become clear in recent years that the widely used atomic Root Mean Square Deviation (RMSD), which is a measure for the precision of the data, overestimates the accuracy of NMR structure ensembles and therefore is a problematic measure for the uncertainty in the atomic coordinates. In this study we report a method that yields a more realistic estimate of the uncertainty in the atomic coordinates by maximizing the RMSD of an ensemble of structures, while maintaining the accordance with the experimentally derived data. The results indicate that the RMSD of most NMR structure ensembles can be significantly increased compromising neither geometric quality nor NMR data. This maximized RMSD therefore seems a better estimate of the true uncertainty in the atomic coordinates.

The high-resolution 3D structure of the octapeptide hormone angiotensin II (AII) in aqueous solution has been obtained by simulated annealing calculations, using high-resolution NMR-derived restraints. After final refinement in explicit water, a family of 13 structures was obtained with a backbone RMSD of 0.73 +/- 0.23 A. AII adopts a fairly compact folded structure, with its C-terminus and N-terminus approaching to within approximately 7.2 A of each other. The side chains of Arg2, Tyr4, Ile5 and His6 are oriented on one side of a plane defined by the peptide backbone, and the Val3 and Pro7 are pointing in opposite directions. The stabilization of the folded conformation can be explained by the stacking of the Val3 side chain with the Pro7 ring and by a hydrophobic cluster formed by the Tyr4, Ile5 and His6 side chains. Comparison between the NMR-derived structure of AII in aqueous solution and the refined crystal structure of the complex of AII with a high-affinity mAb (Fab131) [Garcia, K.C., Ronco, P.M., Verroust, P.J., Brunger, A.T., Amzel, L.M. (1992) Science257, 502-507] provides important quantitative information on two common structural features: (a) a U-shaped structure of the Tyr4-Ile5-His6-Pro7 sequence, which is the most immunogenic epitope of the peptide, with the Asp1 side chain oriented towards the interior of the turn approaching the C-terminus; (b) an Asx-turn-like motif with the side chain aspartate carboxyl group hydrogen-bonded to the main chain NH group of Arg2. It can be concluded that small rearrangements of the epitope 4-7 in the solution structure of AII are required by a mean value of 0.76 +/- 0.03 A for structure alignment and approximately 1.27 +/- 0.02 A for sequence alignment with the X-ray structure of AII bound to the mAb Fab131. These data are interpreted in terms of a biological "nucleus" conformation of the hormone in solution, which requires a limited number of structural rearrangements for receptor-antigen recognition and binding.

The structure determination of protein-protein complexes is a rather tedious and lengthy process, by both NMR and X-ray crystallography. Several methods based on docking to study protein complexes have also been well developed over the past few years. Most of these approaches are not driven by experimental data but are based on a combination of energetics and shape complementarity. Here, we present an approach called HADDOCK (High Ambiguity Driven protein-protein Docking) that makes use of biochemical and/or biophysical interaction data such as chemical shift perturbation data resulting from NMR titration experiments or mutagenesis data. This information is introduced as Ambiguous Interaction Restraints (AIRs) to drive the docking process. An AIR is defined as an ambiguous distance between all residues shown to be involved in the interaction. The accuracy of our approach is demonstrated with three molecular complexes. For two of these complexes, for which both the complex and the free protein structures have been solved, NMR titration data were available. Mutagenesis data were used in the last example. In all cases, the best structures generated by HADDOCK, that is, the structures with the lowest intermolecular energies, were the closest to the published structure of the respective complexes (within 2.0 A backbone RMSD).

We present a CPU efficient protocol for refinement of protein structures in a thin layer of explicit solvent and energy parameters with completely revised dihedral angle terms. Our approach is suitable for protein structures determined by theoretical (e.g., homology modeling or threading) or experimental methods (e.g., NMR). In contrast to other recently proposed refinement protocols, we put a strong emphasis on consistency with widely accepted covalent parameters and computational efficiency. We illustrate the method for NMR structure calculations of three proteins: interleukin-4, ubiquitin, and crambin. We show a comparison of their structure ensembles before and after refinement in water with and without a force field energy term for the dihedral angles; crambin was also refined in DMSO. Our results demonstrate the significant improvement of structure quality by a short refinement in a thin layer of solvent. Further, they show that a dihedral angle energy term in the force field is beneficial for structure calculation and refinement. We discuss the optimal weight for the energy constant for the backbone angle omega and include an extensive discussion of meaning and relevance of the calculated validation criteria, in particular root mean square Z scores for covalent parameters such as bond lengths., (Copyright 2003 Wiley-Liss, Inc.)

The C-terminal domain of the UvrC protein (UvrC CTD) is essential for 5' incision in the prokaryotic nucleotide excision repair process. We have determined the three-dimensional structure of the UvrC CTD using heteronuclear NMR techniques. The structure shows two helix-hairpin-helix (HhH) motifs connected by a small connector helix. The UvrC CTD is shown to mediate structure-specific DNA binding. The domain binds to a single-stranded-double-stranded junction DNA, with a strong specificity towards looped duplex DNA that contains at least six unpaired bases per loop ("bubble DNA"). Using chemical shift perturbation experiments, the DNA-binding surface is mapped to the first hairpin region encompassing the conserved glycine-valine-glycine residues followed by lysine-arginine-arginine, a positively charged surface patch and the second hairpin region consisting of glycine-isoleucine-serine. A model for the protein-DNA complex is proposed that accounts for this specificity.

Nisin is an example of type-A lantibiotics that contain cyclic lanthionine rings and unusual dehydrated amino acids. Among the numerous pore-forming antimicrobial peptides, type-A lantibiotics form an unique family of post-translationally modified peptides. Via the recognition of cell wall precursor lipid II, nisin has the capacity to form pores against Gram-positive bacteria with an extremely high activity in the nanomolar (nM) range. Here we report a high-resolution NMR spectroscopy study of nisin/lipid II interactions in SDS micelles as a model membrane system in order to elucidate the mechanism of molecular recognition at residue level. The binding to lipid II was studied through (15)N-(1)H HSQC titration, backbone amide proton temperature coefficient analysis, and heteronuclear (15)N[(1)H]-NOE relaxation dynamics experiments. Upon the addition of lipid II, significant changes were monitored in the N-terminal part of nisin. An extremely low amide proton temperature coefficient (Delta delta/Delta T) was found for the amide proton of Ala3 (> -0.1 ppb/K) in the complex form. This suggests tight hydrogen bonding and/or isolation from the bulk solvent for this residue. Large chemical shift perturbations were also observed in the first two rings. In contrast, the C-terminal part of nisin was almost unaffected. This part of the molecule remains flexible and solvent-exposed. On the basis of our results, a multistep pore-forming mechanism is proposed. The N-terminal part of nisin first binds to lipid II, and a subsequent structural rearrangement takes place. The C-terminal part of nisin is possibly responsible for the activation of the pore formation. In light of the emerging antibiotic resistance problems, an understanding of the specific recognition mechanism of nisin with lipid II at the residue specific level may therefore aid in the development of novel antibiotics.

The lac repressor-operator system is a model system for understanding protein-DNA interactions and allosteric mechanisms in gene regulation. Despite the wealth of biochemical data provided by extensive mutations of both repressor and operator, the specific recognition mechanism of the natural lac operators by lac repressor has remained elusive. Here we present the first high-resolution structure of a dimer of the DNA-binding domain of lac repressor bound to its natural operator 01. The global positioning of the dimer on the operator is dramatically asymmetric, which results in a different pattern of specific contacts between the two sites. Specific recognition is accomplished by a combination of elongation and twist by 48 degrees of the right lac subunit relative to the left one, significant rearrangement of many side chains as well as sequence-dependent deformability of the DNA. The set of recognition mechanisms involved in the lac repressor-operator system is unique among other protein-DNA complexes and presents a nice example of the adaptability that both proteins and DNA exhibit in the context of their mutual interaction.

The present work describes the dynamics of the apo form of cytochrome b(562), a small soluble protein consisting of 106 amino acid residues [Itagaki, E., and Hager, L. P. (1966) J. Biol. Chem. 241, 3687-3695]. The presence of exchange in the millisecond time scale is demonstrated for the last part of helix IV (residues 95-105 in the holo form). The chemical shift index analysis [Wishart, D. S., and Sykes, B. D. (1994) J. Biomol. NMR 4, 171-180] based on H(alpha), C(alpha), C(beta), and C' chemical shifts suggests a larger helical content than shown in the NMR structure based on NOEs. These results indicate the presence of helical-like conformations participating in the exchange process. This hypothesis is consistent with amide deuterium exchange rates and the presence of some hydrogen bonds identified from amide chemical shift temperature coefficients [Baxter, N. J., and Williamson, M. P. (1997) J. Biomol. NMR 9, 359-369]. (15)N relaxation indicates limited mobility for the amide protons of this part of the helix in the picosecond time scale. A 30 ns stochastic dynamics simulation shows small fluctuations around the helical conformation on this time scale. These fluctuations, however, do not result in a significant decrease of the calculated order parameters which are consistent with the experimental (15)N relaxation data. These results resolve an apparent discrepancy in the NMR structures between the disorder observed in helix IV due to a lack of NOEs and the secondary structure predictions based on H(alpha) chemical shifts [Feng, Y., Wand, A. J., and Sligar, S. G. (1994) Struct. Biol. 1, 30-35].

The possibility of generating protein folds at the stage of backbone assignment using structural restraints derived from experimentally measured cross-hydrogen bond scalar couplings and secondary chemical shift information is investigated using as a test case the small alpha/beta protein chymotrypsin inhibitor 2. Dihedral angle restraints for the phi and psi angles of 32 out of 64 residues could be obtained from secondary chemical shift analysis with the TALOS program (Corneliscu et al., 1999a). This information was supplemented by 18 hydrogen-bond restraints derived from experimentally measured cross-hydrogen bond 3hbJNC' coupling constants. These experimental data were sufficient to generate structures that are as close as 1.0 A backbone rmsd from the crystal structure. The fold is, however, not uniquely defined and several solutions are generated that cannot be distinguished on the basis of violations or energetic considerations. Correct folds could be identified by combining clustering methods with knowledge-based potentials derived from structural databases.

The retinoid X receptor (RXR) is a prominent member of the nuclear receptor family of ligand-inducible transcription factors. Many proteins of this family exert their function as heterodimers with RXR as a common upstream partner. Studies of the DNA-binding domains of several nuclear receptors reveal differences in structure and dynamics, both between the different proteins and between the free- and DNA-bound receptor DBDs. We investigated the differences in dynamics between RXR free in solution and in complex with a 14 base-pair oligonucleotide, using (1)H and (15)N relaxation studies. Nano- to picosecond dynamics were probed on (15)N, employing Lipari-Szabo analysis with an axially symmetric tumbling model to estimate the exchange contributions to the transverse relaxation rates. Furthermore, milli- to microsecond dynamics were estimated qualitatively for (1)H and (15)N, using CPMG-HSQC and CPMG-T(2) measurements with differential pulse spacing. RXR shows hardly any nano- to picosecond time-scale internal motion. Upon DNA binding, the order parameters show a tiny increase. Dynamics in the milli- to microsecond time scale is more prevalent. It is localized in the first and second zinc fingers of the free RXR. Upon DNA-binding, exchange associated with specific/aspecific DNA-binding of RXR is observed throughout the sequence, whereas conformational flexibility of the D-box and the second zinc finger of RXR is greatly reduced. Since this DNA-binding induced folding transition occurs remote from the DNA in a region which is involved in protein-protein interactions, it may very well be related to the cooperativity of dimeric DNA binding.

The hydration of the collagen-like Ac-(Gly-Pro-Hyp)(6)-NH(2) triple-helical peptide in solution was investigated using an integrated set of high-resolution NMR hydration experiments, including different recently developed exchange-network editing methods. This approach was designed to explore the hydration dynamics in the proximity of labile groups, such as the hydroxyproline hydroxyl group, and revealed that the first shell of hydration in collagen-like triple helices is kinetically labile with upper limits for water molecule residence times in the nanosecond to sub-nanosecond range. This result is consistent with a "hopping" hydration model in which solvent molecules are exchanged in and out of solvation sites at a rate that is not directly correlated to the degree of site localization. The hopping model thus reconciles the dynamic view of hydration revealed by NMR with the previously suggested partially ordered semi-clathrate-like cylinder of hydration. In addition, the nanosecond to sub-nanosecond upper limits for water molecule residence times imply that hydration-dehydration events are not likely to be the rate-limiting step for triple helix self-recognition, complementing previous investigations on water dynamics in collagen fibers. This study has also revealed labile proton features expected to facilitate the characterization of the structure and folding of triple helices in collagen peptides., (Copyright 2000 Academic Press.)

The stability and (un)folding of the 19-residue peptide, SCVTLYQSWRYSQADNGCA, corresponding to the first beta-hairpin (residues 10 to 28) of the alpha-amylase inhibitor tendamistat (PDB entry 3AIT) has been studied by molecular dynamics simulations in explicit water under periodic boundary conditions at several temperatures (300 K, 360 K and 400 K), starting from various conformations for simulation lengths, ranging from 10 to 30 ns. Comparison of trajectories of the reduced and oxidized native peptides reveals the importance of the disulphide bridge closing the beta-hairpin in maintaining a proper turn conformation, thereby insuring a proper side-chain arrangement of the conserved turn residues. This allows rationalization of the conservation of those cysteine residues among the family of alpha-amylase inhibitors. High temperature simulations starting from widely different initial configurations (native beta-hairpin, alpha and left-handed helical and extended conformations) begin sampling similar regions of the conformational space within tens of nanoseconds, and both native and non-native beta-hairpin conformations are recovered. Transitions between conformational clusters are accompanied by an increase in energy fluctuations, which is consistent with the increase in heat capacity measured experimentally upon protein folding. The folding events observed in the various simulations support a model for beta-hairpin formation in which the turn is formed first, followed by hydrogen bond formation closing the hairpin, and subsequent stabilization by side-chain hydrophobic interactions., (Copyright 2000 Academic Press.)

The localisation and dynamics of sodium counterions around the DNA duplex d(AGCGTACTAGTACGCT)2 corresponding to the trp operator fragment used in the crystal structure of the half site complex has been studied by a 1.4 ns molecular dynamics simulation in explicit solvent. A continuous and well-defined counterion density is shown to be present around the minor groove, while density patches are found in the major groove in regions where DNA bending is observed. A residence time analysis reveals the dynamic nature of these distributions. The resulting picture agrees with previous theoretical and experimental studies of A-tract DNA sequences, and is consistent with the polyelectrolyte condensation model.

Background: Lactose repressor protein (Lac) controls the expression of the lactose metabolic genes in Escherichia coli by binding to an operator sequence in the promoter of the lac operon. Binding of inducer molecules to the Lac core domain induces changes in tertiary structure that are propagated to the DNA-binding domain through the connecting hinge region, thereby reducing the affinity for the operator. Protein-protein and protein-DNA interactions involving the hinge region play a crucial role in the allosteric changes occurring upon induction, but have not, as yet, been analyzed in atomic detail., Results: We have used nuclear magnetic resonance (NMR) spectroscopy and restrained molecular dynamics (rMD) to determine the structure of the Lac repressor DNA-binding domain (headpeice 62; HP62) in complex with a symmetrized lac operator. Analysis of the structures reveals specific interactions between Lac repressor and DNA that were not found in previously investigated Lac repressor-DNA complexes. Important differences with the previously reported structures of the HP56-DNA complex were found in the loop following the helix-turn-helix (HTH) motif. The protein-protein and protein-DNA interactions involving the hinge region and the deformations in the DNA structure could be delineated in atomic detail. The structures were also used for comparison with the available crystallographic data on the Lac and Pur repressor-DNA complexes., Conclusions: The structures of the HP62-DNA complex provide the basis for a better understanding of the specific recognition in the Lac repressor-operator complex. In addition, the structural features of the hinge region provide detailed insight into the protein-protein and protein-DNA interactions responsible for the high affinity of the repressor for operator DNA.

The structure and hydration of the DNA duplex d-(AGCGTACTAGTACGCT)2 corresponding to the trp operator fragment used in the crystal structure of the half site complex (PDB entry 1TRR) was studied by a 1.4 ns molecular dynamics simulation in water. The simulation, starting from a B-DNA conformation, used a non-bonded cutoff of 1.4 nm with a reaction field correction and resulted in a stable trajectory. The average DNA conformation obtained was closer to the ones found in the crystal structures of the complexes (PDB entries 1TRO and 1TRR) than to the crystal structure of unbound trp operator (Nucleic Acid Database entry BDJ061). The DNA hydration was characterized in terms of hydrogen bond percentages and corresponding residence times. The residence times of water molecules within 0.35 nm of the DNA non-exchangeable protons were calculated for comparison with NMR measurements of intermolecular water-DNA NOEs and nuclear magnetic relaxation dispersion measurements. No significant difference was found between major and minor groove hydration. The DNA donors and acceptors were hydrogen bonded to water molecules for 77(+/-19)% of the time on average. The average residence time of the hydrogen bonded water molecules was 11(+/-11) ps with a maximum of 223 ps. When all water molecules within NOE distance (0.35 nm) of non-exchangeable protons were considered, the average residence times increased to an average of 100(+/-4) ps and a maximum of 608 ps. These results agree with the experimental NMR results of Sunnerhagen et al. which did not show any evidence for water molecules bound with more than 1 ns residence time on the DNA surface. The exchange of hydration water from the DNA occurred in the major groove primarily through direct exchange with the bulk solvent, while access to and from the minor groove frequently proceeded via pathways involving ribose O3' and O4' and phosphate O2P oxygen atoms. The most common water diffusion pathways in the minor groove were perpendicular to the groove direction. In general, water molecules visited only a limited number of sites in the DNA grooves before exiting. The hydrogen bonding sites, where hydrogen bonds could be formed with donor and acceptor groups of the DNA, were filled with water molecules with an average B-factor value of 0.58 mn2. No special values were observed at any of the sites, where water molecules were observed both in the trp repressor/operator co-crystals and in the crystal structure of unbound DNA., (Copyright 1998 Academic Press.)

The present NMR study investigates the residence times of the hydration water molecules associated with uncomplexed trp operator DNA in solution by measuring intermolecular nuclear Overhauser effects (NOE) between water and DNA protons, and the nuclear magnetic relaxation dispersion (NMRD) of the water 2H and 17O resonances. Both methods indicate that the hydration water molecules exchange with bulk water on the sub-nanosecond time scale at 4 degreesC. No evidence was obtained for water molecules bound with longer residence times. In particular, the water molecules at the sites of interfacial hydration in the trp repressor/operator complex do not seem kinetically stabilized in the uncomplexed DNA. Analysis of the crystal structures of two different trp repressor/operator complexes shows very similar structural environments for the water molecules mediating specific contacts between the protein and the DNA, whereas much larger variations are observed for the location of corresponding water molecules detected in the crystal structure of an uncomplexed trp operator DNA duplex. Therefore, it appears unlikely that the hydration characteristics of the uncomplexed DNA target would be a major determinant of trp repressor/operator recognition., (Copyright 1998 Academic Press.)

The solution structure of the DNA binding domain of lac repressor (headpiece 1-56; HP56) has been refined using data from 2D and 3D NMR spectroscopy. The structure was derived from 1546 restraints (giving an average of 27.6 per residue), comprising 389 intraresidual, 402 sequential, 385 medium range and 325 long range distance restraints and also 30 phi and 15 chi 1 dihedral angle restraints. The structures were determined by the method of direct refinement against nuclear Overhauser enhancement peak volumes with the program DINOSAUR. The final set of 32 selected structures displayed an r.m.s. deviation from the average of 0.43(+/-0.08) A angstroms (backbone) and 0.95(+/-0.08) angstroms (all heavy atoms) for the best defined region of the protein (residues 3 to 49). The ensemble R-factor was 0.35, which indicates close correspondence with the experimental data. The structures revealed good stereochemical qualities. The conformations of the NMR structures of free and DNA complexed lac repressor headpiece were compared. The regions comprising the secondary structure elements show close correspondence for both conformations. However, the conformation of the loop between helix II and III changes considerably upon complexation of the headpiece. This change in the conformation of the loop in lac HP56 is essential for binding of the side-chains of residues Asn25 and His29 to the lac operator DNA. Finally, the lac headpiece residues that are intolerant to mutations were analysed. Most of these mutation-sensitive residues are important for a correct folding of the headpiece region, and a number of these residues are also involved in contacting the operator DNA.

The feasibility of determining the relative populations of multi-conformer structures from NOE-derived distances alone is assessed. Without cross-validation of the NOE restraints, any population ratio can be refined to a similar quality of the fit. Complete cross-validation provides a less biased measure of fit and allows the estimation of the correct population ratio when used in conjunction with very tight distance restraints. With the qualitative distance restraints most commonly used in NMR structure determination, cross-validation is unsuccessful in providing the correct answer. Other experimental sources are therefore needed to determine relative populations of multi-conformer structures.

In structure determination by X-ray crystallography and solution NMR spectroscopy, experimental data are collected as time and ensemble-averages. Thus, in principle, appropriate time and ensemble-averaged models should be used. Refinement of an ensemble of conformers rather than one single structure against the experimental NMR data could, however, result in overfitting the data because of the significantly increased number of parameters. To avoid overfitting, complete cross-validation, which provides an unbiased measure of the fit, has been applied to nuclear Overhauser effect derived distance refinement. Using two synthetic test cases, a correlation was demonstrated between the cross-validated measure to the fit (defined in terms of root-mean-square deviations from the distance restraints and number of violations) and the number of models that best reproduce the conformational variability in solution. A new method, based on a probability map, has been used to generate good representations of the resulting ensembles of structures. The method has also been applied to observed NMR data for two proteins, interleukin 4 and interleukin 8. For interleukin 4, cross-validation indicates that a single-conformer model gives the most accurate representation of the structure, whereas conventional measures of fit between the experimental data and those calculated from the model decrease when increasing the number of conformers, indicating overfitting. For interleukin 8, complete cross-validation predicts a twin-conformer model to be the most faithful representation of the experimental data. Two distinct conformations for the loop formed by residues 16 to 22 emerge from the family of twin-conformer structures. The putative alternate conformation of the loop is not observed in the crystal structure of interleukin 8. However, because of crystal packing contacts in this region this does not necessarily exclude the presence of the alternate conformation in solution. The twin-conformer model is supported by observed chemical exchange line broadening for the amide of His18 obtained by 15N relaxation studies. This region has also been implied to be involved in receptor binding.

The solution structure of the glucocorticoid receptor (GR) DNA-binding domain (DBD), consisting of 93 residues, has been refined from two and three-dimensional NMR data using an ensemble iterative relaxation matrix approach followed by direct NOE refinement with DINOSAUR. A set of 47 structures of the rat GR fragment Cys440-Arg510 was generated with distance geometry and further refined with a combination of restrained energy minimization and restrained molecular dynamics in a parallel refinement protocol. Distance constraints were obtained from an extensive set of NOE build-up curves in H2O and 2H2O via relaxation matrix calculations (1186 distance constraints from NOE intensities, 10 phi and 22 chi 1 dihedral angle constraints from J- coupling data were used for the calculations). The root-mean-square deviation values of the 11 best structures on the well-determined part of the protein (Cys440 to Ser448, His451 to Glu469 and Pro493 to Glu508) are 0.60 A and 1.20 A from the average for backbone and all heavy atoms, respectively. The final structures have R-factors around 0.40 and good stereochemical qualities. The first zinc-coordinating domain of the GR DBD is very similar to the crystal structure with a root-mean-square difference of 1.4 A. The second zinc-coordinating domain is still disordered in solution. No secondary structure element is found in this domain in the free state. As suggested by crystallographic studies on the estrogen receptor DBD-DNA and GR DBD-DNA complexes, part of this region will form a distorted helix and the D-box will undergo a conformational change upon cooperative binding to DNA.

The Arc repressor of Salmonella bacteriophage P22 is a dimeric sequence-specific DNA-binding protein. The solution structure of Arc has been determined from 2D NMR data using an "ensemble" iterative relaxation matrix approach (IRMA) followed by direct NOE refinement with DINOSAUR. A set of 51 structures was generated with distance geometry and further refined with a combination of restrained energy minimization and restrained molecular dynamics in a parallel refinement protocol. Distance constraints were obtained from an extensive set of NOE build-ups in H2O and 2H2O via relaxation matrix calculations from the ensemble of structures. Methyl group rotation, aromatic ring flaps and internal mobility effects (via order parameters obtained from a free molecular dynamics run in water) were included in these calculations. The best structures were finally refined with direct NOE constraints following a slow-cooling simulated annealing protocol. In this final refinement stage, theoretical NOE intensities were directly compared with the experimental data and forces were derived using a simple two-spin approximation for the gradient of the NOE function. Dynamic assignment was applied to the peaks involving unassigned diastereotopic groups. The structure is determined to a precision (r.m.s.d. from the average excluding the ill defined C and N-terminal region) of 0.55 and 1.10 A for backbone and all atoms, respectively. The final structures, with R factor values around 0.35, have good stereochemical qualities, contain an extensive network of hydrogen bonds consistent with the secondary structure elements and structural features in concordance with genetic data. The overall folding of the solution and crystal structures is the same.

NMR data are collected as time- and ensemble-averaged quantities. Yet, in commonly used methods for structure determination of biomolecules, structures are required to satisfy simultaneously a large number of constrainsts. Recently, however, methods have been developed that allow a better fit of the experimental data by the use of time- or ensemble-averaged restraints. Thus far, these methods have been applied to structure refinement using distance and J-coupling restraints. In this paper, time and ensemble averaging is extended to the direct refinement with experimental NOE data. The implementation of time- and ensemble-averaged NOE restraints in DINOSAUR is described and illustrated with experimental NMR data for crambin, a 46-residue protein. Structure refinement with both time- and ensemble-averaged NOE restraints results in lower R-factors, indicating a better fit of the experimental NOE data.

The structure in solution of crambin, a small protein of 46 residues, has been determined from 2D NMR data using an iterative relaxation matrix approach (IRMA) together with distance geometry, distance bound driven dynamics, molecular dynamics, and energy minimization. A new protocol based on an "ensemble" approach is proposed and compared to the more standard initial rate analysis approach and a "single structure" relaxation matrix approach. The effects of fast local motions are included and R-factor calculations are performed on NOE build-ups to describe the quality of agreement between theory and experiment. A new method for stereospecific assignment of prochiral groups, based on a comparison of theoretical and experimental NOE intensities, has been applied. The solution structure of crambin could be determined with a precision (rmsd from the average structure) of 0.7 A on backbone atoms and 1.1 A on all heavy atoms and is largely similar to the crystal structure with a small difference observed in the position of the side chain of Tyr-29 which is determined in solution by both J-coupling and NOE data. Regions of higher structural variability (suggesting higher mobility) are found in the solution structure, in particular for the loop between the two helices (Gly-20 to Pro-22).

The three-dimensional structure of the DNA-binding domain of the human retinoic acid receptor-beta (hRAR-beta) has been determined by nuclear magnetic resonance spectroscopy in conjunction with distance geometry, restrained molecular dynamics and iterative relaxation matrix calculations. A total of 1244 distance restraints were obtained from NOE intensities, of which 448 were intra-residue and 796 inter-residue restraints. In addition 23 chi and 30 phi dihedral angle restraints were obtained from J-coupling data. The two 'zinc-finger' regions of the 80-amino acid residue protein are followed by two alpha-helices that cross each other perpendicularly. There is a short stretch of b-sheet near the N-terminus. The alpha-helical core of the protein is well determined with a backbone root-mean-square deviation (r.m.s.d.) with respect to the average of 0.18 A and 0.37 A when the side chains of residues 31, 32, 36, 61, 62, 65 and 69 are included. The r.m.s.d. for the backbone of residues 5-80 is 0.76 A. For the first finger (residues 8-28), the r.m.s.d. of the backbone is 0.79 A. For the second finger (residues 44-62) the r.m.s.d. is 0.64 A. The overall structure is similar to that of the corresponding domain of the glucocorticoid receptor, although the C-terminal part of the protein is different. The second alpha-helix is two residues shorter and is followed by a well-defined region of extended backbone structure.

A set of computer programs called DINOSAUR has been developed, which allows the refinement of biomolecular structures directly from 2D NOE intensities. The NOE restraining potential implemented emphasises the weak intensities corresponding to larger distances which are more likely to determine the three-dimensional structure. An approximation based on a two-spin approximation is proposed for the gradient of the NOE intensities instead of the exact solution which is extremely time-consuming. The DINOSAUR routines have been implemented in various refinement programs (Distance bound Driven Dynamics, Molecular Dynamics and Energy Minimisation) and tested on an eight-residue model peptide.