Sequence-specific resonance assignments in the 1H nuclear-magnetic-resonance spectrum of the lac repressor DNA-binding domain 1-51 from Escherichia coli by two-dimensional spectroscopy

Brownian motion, spin diffusion and protein structure determination in solution

This paper presents my recollections on the development of protein structure determination by NMR in solution from 1968 to 1992. The key to success was to identify NMR-accessible parameters that unambiguously determine the spatial arrangement of polypeptide chains. Inspired by work with cyclopeptides, model considerations showed that enforcing short non-bonding interatomic distances imposes «ring closure conditions» on polypeptide chains. Given that distances are scalar parameters, this indicated an avenue for studies of proteins in solution, i.e., under the regime of stochastic rotational and translational motions at frequencies in the nanosecond range (Brownian motion), where sharp pictures could not be obtained by photography-related methods. Later-on, we used distance geometry calculations with sets of inter-atomic distances derived from protein crystal structures to confirm that measurements of short proton-proton distances could provide atomic-resolution structures of globular proteins. During the years 1976-1984 the following four lines of research then led to protein structure determination by NMR in solution. First, the development of NMR experiments enabling the use of the nuclear Overhauser effect (NOE) for measurements of interatomic distances between pairs of hydrogen atoms in proteins. Second, obtaining sequence-specific resonance assignment solved the "phase problem" for protein structure determination by NMR. Third, generating and programming novel distance geometry algorithms enabled the calculation of atomic-resolution protein structures from limited sets of distance constraints measured by NMR. Fourth, the introduction of two-dimensional NMR provided greatly improved spectral resolution of the complex spectra of proteins as well as efficient delineation of scalar and dipole-dipole ¹H-¹H connectivities, thus making protein structure determination in solution viable and attractive.

Introduction to a special issue of Magnetic Resonance in honour of Robert Kaptein at the occasion of his 80th birthday

This publication, in honour of Robert Kaptein's 80th birthday, contains contributions from colleagues, many of whom have worked with him, and others who admire his work and have been stimulated by his research. The contributions show current research in biomolecular NMR, spin hyperpolarisation and spin chemistry, including CIDNP (chemically induced dynamic nuclear polarisation), topics to which he has contributed enormously. His proposal of the radical pair mechanism was the birth of the field of spin chemistry, and the laser CIDNP NMR experiment on a protein was a major breakthrough in hyperpolarisation research. He set milestones for biomolecular NMR by developing computational methods for protein structure determination, including restrained molecular dynamics and 3D NMR methodology. With a lac repressor headpiece, he determined one of the first protein structures determined by NMR. His studies of the lac repressor provided the first examples of detailed studies of protein nucleic acid complexes by NMR. This deepened our understanding of protein DNA recognition and led to a molecular model for protein sliding along the DNA. Furthermore, he played a leading role in establishing the cluster of NMR large-scale facilities in Europe. This editorial gives an introduction to the publication and is followed by a biography describing his contributions to magnetic resonance.

T1BT* structural study of an anti-plasmodial peptide through NMR and molecular dynamics

Background

T1BT* is a peptide construct containing the T1 and B epitopes located in the 5’ minor repeat and the 3’ major repeat of the central repeat region of the Plasmodium falciparum circumsporozoite protein (CSP), respectively, and the universal T* epitope located in the C-terminus of the same protein. This peptide construct, with B = (NANP)₃, has been found to elicit antisporozoite antibodies and gamma-interferon-screening T-cell responses in inbred strains of mice and in outbred nonhuman primates. On the other hand, NMR and CD spectroscopies have identified the peptide B’ = (NPNA)₃ as the structural unit of the major repeat in the CSP, rather than the more commonly quoted NANP. With the goal of assessing the structural impact of the NPNA cadence on a proven anti-plasmodial peptide, the solution structures of T1BT* and T1B’T* were determined in this work.

Methods

NMR spectroscopy and molecular dynamics calculations were used to determine the solution structures of T1BT* and T1B’T*. These structures were compared to determine the main differences and similarities between them.

Results

Both peptides exhibit radically different structures, with the T1B’T* showing strong helical tendencies. NMR and CD data, in conjunction with molecular modelling, provide additional information about the topologies of T1BT* and T1B’T*. Knowing the peptide structures required to elicit the proper immunogenic response can help in the design of more effective, conformationally defined malaria vaccine candidates. If peptides derived from the CSP are required to have helical structures to interact efficiently with their corresponding antibodies, a vaccine based on the T1B’T* construct should show higher efficiency as a pre-erythrocyte vaccine that would prevent infection of hepatocytes by sporozoites.

Effect of solvent on an NMR chemical shift difference between glycyl geminal alpha-protons as a probe of beta-turn formation of short peptides

Proton NMR spectra of short peptides with a glycyl (Gly) or N-methylglycyl (sarcosyl, Sar) residue were measured in various mixed solvents with a wide range of dielectric constants: 78.3-2.3. From analyses of the octet and quartet signals of the geminal alpha-protons of Gly and Sar residues, respectively, we have estimated chemical shift differences between the two alpha-protons, Delta delta(alpha/alpha'). It is found that the Delta delta(alpha/alpha') values increase with decreasing solvent polarity and the increasing rates depend significantly on amino acid sequences. By referring to infrared spectra and chemical shift of the terminal NH protons, deltaNH, of the peptides, the Delta delta(alpha/alpha') values were found to be a good probe of beta-turn formation. From solvent-dependent change of Delta delta(alpha/alpha'), we estimated the free energies for the beta-turn formation and compared the results with those estimated from deltaNH. Using the resulting free energies, we have discussed effects of solvent on the beta-turn formation.

Graph-theoretical assignment of secondary structure in multidimensional protein NMR spectra: application to the lac repressor headpiece

A novel procedure is presented for the automatic identification of secondary structures in proteins from their corresponding NOE data. The method uses a branch of mathematics known as graph theory to identify prescribed NOE connectivity patterns characteristic of the regular secondary structures. Resonance assignment is achieved by connecting these patterns of secondary structure together, thereby matching the connected spin systems to specific segments of the protein sequence. The method known as SERENDIPITY refers to a set of routines developed in a modular fashion, where each program has one or several well-defined tasks. NOE templates for several secondary structure motifs have been developed and the method has been successfully applied to data obtained from NOESY-type spectra. The present report describes the application of the SERENDIPITY protocol to a 3D NOESY-HMQC spectrum of the 15N-labelled lac repressor headpiece protein. The application demonstrates that, under favourable conditions, fully automated identification of secondary structures and semi-automated assignment are feasible.

Polyanion-induced alpha-helical structure of a synthetic 23-residue peptide representing the lysine-rich segment of the N-terminal extension of yeast cytoplasmic aspartyl-tRNA synthetase

Conformational studies were performed on the synthetic tricosapeptide N-acetyl-SKKALKKLQKEQEKQRKKEERAL-amide, representing the highly basic segment (residues 30-52) of the N-terminal extension of yeast cytoplasmic aspartyl-tRNA synthetase. Circular dichroism experiments show that, in aqueous solution at neutral pH, the peptide adopts a random conformation. The effects of pH, temperature, addition of trifluoroethanol (TFE), and titration with polyanions on the conformation of the peptide were studied. In TFE or in the presence of an equimolar concentration of (phosphate)18, the peptide adopts a 100% alpha-helical conformation. A partially alpha-helical conformation is induced by (phosphate)4 or d(pT)8 (respectively 40% and 35% helical content). Raising the pH in aqueous solution promotes 75% alpha-helicity, with a transition pK of 9.9 reflecting deprotonation of lysine residues. On the basis of these results, nuclear magnetic resonance studies were carried out in TFE as well as in aqueous solution in the presence of (phosphate)18, to determine the structure of the molecule. Complete 1H resonance assignments were obtained by conventional two-dimensional NMR techniques. A total of 138 interproton constraints derived from NOESY experiments were used to calculate the three-dimensional structure by a two-stage distance geometry/simulated annealing procedure. The two deduced structures were highly similar and show that nine cationic residues are segregated on one face of a helical structure, providing an ideal polycationic interface for binding to polyanionic surfaces.

On the interpretation of biochemical data by molecular dynamics computer simulation

The application of computer simulation to molecular systems of biochemical interest is reviewed. It is shown that computer simulation is a tool complementary to experimental methods, which can be used to access atomic details inaccessible to experimental probes. Examples are given in which computer simulation augments the experimental information by providing an atomic picture of high resolution with respect to space, energy or time. The usefulness of a computer simulation largely depends on its quality. The most important factors that limit the accuracy of simulated results are discussed. The accuracy of different simulation studies can differ by orders of magnitude. The accuracy will depend on the type of biomolecular system and process studied. It will also depend on the choice of force field, the simulation set-up and the protocol that is used. A list of quality-determining factors is given, which may be useful when interpreting simulation studies appearing in the literature.

The chemical shift index: a fast and simple method for the assignment of protein secondary structure through NMR spectroscopy

Previous studies by Wishart et al. [Wishart, D. S., Sykes, B. D., & Richards, F. M. (1991) J. Mol. Biol. (in press)] have demonstrated that 1H NMR chemical shifts are strongly dependent on the character and nature of protein secondary structure. In particular, it has been found that the 1H NMR chemical shift of the alpha-CH proton of all 20 naturally occurring amino acids experiences an upfield shift (with respect to the random coil value) when in a helical configuration and a comparable downfield shift when in a beta-strand extended configuration. On the basis of these observations, a technique is described for rapidly and quantitatively determining the identity, extent, and location of secondary structural elements in proteins based on the simple inspection of the alpha-CH 1H resonance assignments. A number of examples are provided to demonstrate both the simplicity and the accuracy of the technique. This new method is found to be almost as accurate as the more traditional NOE-based methods of determining secondary structure and could prove to be particularly useful in light of the recent development of sequential assignment techniques which are now almost NOE-independent [Ikura, M., Kay, L. E., & Bax, A. (1990) Biochemistry 29, 4659-4667]. We suggest that this new procedure should not necessarily be seen as a substitute to existing rigorous methods for secondary structure determination but, rather, should be viewed as a complement to these approaches.

Relationship between nuclear magnetic resonance chemical shift and protein secondary structure

An analysis of the 1H nuclear magnetic resonance chemical shift assignments and secondary structure designations for over 70 proteins has revealed some very strong and unexpected relationships. Similar studies, performed on smaller databases, for 13C and 15N chemical shifts reveal equally strong correlations to protein secondary structure. Among the more interesting results to emerge from this work is the finding that all 20 naturally occurring amino acids experience a mean alpha-1H upfield shift of 0.39 parts per million (from the random coil value) when placed in a helical configuration. In a like manner, the alpha-1H chemical shift is found to move downfield by an average of 0.37 parts per million when the residue is placed in a beta-strand or extended configuration. Similar changes are also found for amide 1H, carbonyl 13C, alpha-13C and amide 15N chemical shifts. Other relationships between chemical shift and protein conformation are also uncovered; in particular, a correlation between helix dipole effects and amide proton chemical shifts as well as a relationship between alpha-proton chemical shifts and main-chain flexibility. Additionally, useful relationships between alpha-proton chemical shifts and backbone dihedral angles as well as correlations between amide proton chemical shifts and hydrogen bond effects are demonstrated.

Simple techniques for the quantification of protein secondary structure by 1H NMR spectroscopy

Previous work by Wishart et al. (in press) and others [(1989) J. Magn. Reson. 83, 441-449; (1990) J. Magn. Reson. 90, 165-176] has shown a strong tendency for protein secondary structure to be manifested in 1H NMR chemical shifts. Based on these earlier results, two techniques have been developed for the quantification of secondary structure in proteins. Both methods allow for the rapid and accurate determination of the percent content of helix, coil, and beta-strand based on the integration (or peak enumeration) of selected portions of either 1-D or 2-D 1H NMR spectra. These new and very simple procedures have been found to compare quite favorably to other well established techniques for secondary structure determination such as CD, Raman and IR spectroscopy.

Three-dimensional correlated NMR study of Megasphaera elsdenii flavodoxin in the oxidized state

The value of a three-dimensional (3D) non-selective total correlation/nuclear Overhauser enhancement spectroscopy (TOCSY-NOESY) spectrum for making sequential resonance assignments in proteins is demonstrated using the relatively large Megasphaera elsdenii flavodoxin (molecular mass 15 kDa) in the oxidized state. An easy and concise method for the analysis of 3D-NMR spectra and a strategy for the resonance assignment of 3D-NMR protein spectra is introduced. In this context, non-selective TOCSY-NOESY is compared with selective TOCSY-NOESY and non-selective NOESY-TOCSY. Sequential assignments in various secondary structure elements of flavodoxin are made using the method of analysis introduced. NOEs not previously identified in 2D-NMR spectra due to resonance overlap are found in the 3D Clean-TOCSY-NOESY spectrum. Also additional side-chain assignments could be made.

Theoretical variation of the H alpha chemical shift in acetyl-glycyl-N-methylamide and oligoglycines with molecular conformation and environment

The sum of the magnetic anisotropy and polarization contributions to the magnetic shielding constants of the alpha protons is calculated as a function of the torsion angles about the NC alpha (phi) and C alpha C' (psi) bonds of the dipeptide. The results show that the polarization or electric field effect is several fold larger than the magnetic anisotropy contribution. The calculated variations are large enough to account for the spread of the values measured for these protons in peptides and proteins. The results obtained for polyglycine alpha helices and antiparallel beta sheets are discussed in relation with molecular conformation and environmental effects on the one hand and experimental data on the other.

A two-dimensional 1H NMR study on Megasphaera elsdenii flavodoxin in the reduced state. Sequential assignments

Assignments for the 137 amino acid residues of Megasphaera elsdenii flavodoxin in the reduced state have been made using the sequential resonance assignment procedure. Several hydroxyl and sulfhydryl protons were observed at 41 degrees C at pH 8.3. Spin systems were sequentially assigned using phase-sensitive two-dimensional-correlated spectroscopy and phase-sensitive nuclear Overhauser enhancement spectroscopy. Spectra of the protein in H2O and of protein preparations either completely or partly exchanged against 2H2O were obtained. Use of the fast electron shuttle between the paramagnetic semiquinone and the diamagnetic hydroquinone state greatly simplified the NMR spectra, making it possible to assign easily the 1H resonances of amino acid residues located in the immediate neighbourhood of the isoalloxazine ring. The majority of the nuclear Overhauser effect contracts between the flavin and the apoprotein correspond to the crystal structure of the flavin domain of Clostridium MP flavodoxin, but differences are also observed. The assignments provide the basis for the structure determination of M. elsdenii flavodoxin in the reduced state as well as for assigning the resonances of the oxidized flavodoxin.

Proton nuclear magnetic resonance assignments

The procedures outlined here have been used successfully for more than 30 proteins to date, and are nearly routine for molecules up to a molecular weight of 10,000. Some of the proteins assigned have a molecular weight greater than 10,000. For these larger proteins, relayed-COSY and TOCSY experiments have been essential for the identification of spin systems, although for thioredoxin these experiments could not be used. In this case, assignments were accomplished using nonspecific deuteration to the level of 75% and specific, nearly complete, deuteration of certain kinds of residues (see LeMaster [2], this volume). Nonspecific deuteration reduces the cross-relaxation rates of each proton to the rest of the molecule, thus reducing the linewidths. The cross-peak patterns were also narrowed due to simplification of the coupling patterns. Such a laborious procedure of nonspecific deuteration may not be necessary for complete proton assignments of proteins in this size range, as evidenced by the fact that this method was not used for the other two molecules mentioned above. It may prove, however, to be quite valuable in the study of larger molecules, where linewidths are expected to increase due to longer rotational correlation times. Overlap problems in the NH chemical shifts can be dealt with by making use of the differential temperature dependence of these shifts. Another technique is to take advantage of the wide range of exchange rates between these protons and the solvent. Spectra containing only the slowly exchanging NH protons can be obtained by acquiring spectra of the protein soon after dilution in D2O, and spectra of only the rapidly exchanging protons can be obtained by obtaining spectra in a freshly prepared H2O solution of the protein after having completely exchanged all the NH protons with deuterium. Variation of the pH will resolve problems of overlap in all regions of the spectrum, although many chemical shifts may be unaffected by pH. In some cases, pH variation may change the conformation of the molecule. This may, in fact, assist in the sequential assignment if the chemical shifts can be followed with pH. Finally, the relayed-NOESY experiments can resolve overlap problems with the alpha-proton chemical shifts. Thus, it is very likely that the assignment methods outlined here will be successful for the assignment of the proton spectra of even larger molecules if there is significant secondary structure and significant variety of residues to provide enough dispersion of the chemical shifts.(ABSTRACT TRUNCATED AT 400 WORDS)

Photo-CIDNP study of the interaction between lac repressor headpiece and lac operator DNA

Lac repressor headpiece (HP) and intact lac repressor have been studied using the photo-CIDNP method. At neutral pH histidine 29, tyrosines 7, 12 and 17 and methionine 1 are polarised. His-29 polarizations are weaker and broader in HP59 than in HP51 indicating that the C-terminal octapeptide in HP59 adopts a conformation that allows an interaction with His-29. The photo-CIDNP spectra of intact lac repressor and HP51 are very similar, showing that the same residues are accessible to the photo-excited flavin. An equimolar mixture of HP51 and a 14 base pair lac operator fragment strongly suppresses the photo-CIDNP effect of tyrosines 7 and 17 and abolishes the His-29 polarizations. The results are compared with earlier photo-CIDNP measurements on a complex of headpiece with poly[d(AT)] and with a model derived from a 2D NMR study on a lac headpiece-operator complex.

Protein NMR resonance assignment by isotropic mixing experiments on random fractionally deuterated samples

The 108-residue protein E. coli thioredoxin has been uniformly enriched to 50% with deuterium at all carbon-bound hydrogen positions. Isotropic mixing (i.e. TOCSY) experiments have been conducted for both the deuterated and natural-abundance samples. Using a 54 ms mixing time correlation peaks can be seen for all four protons on the benzenoid ring of tryptophan in both samples. The deuteration results in an average decrease in cross-sectional area of a factor of 2-3 for the TOCSY cross-peaks. The cross-peak intensities for the deuterated sample systematically decrease as a function of the number of protons involved in the transfer process thus overcoming a common ambiguity in the TOCSY experiment.

Flexibility of the DNA-binding domains of trp repressor

An orthorhombic crystal form of trp repressor (aporepressor plus L-tryptophan ligand) was solved by molecular replacement, refined to 1.65 A resolution, and compared to the structure of the repressor in trigonal crystals. Even though these two crystal forms of repressor were grown under identical conditions, the refined structures have distinctly different conformations of the DNA-binding domains. Unlike the repressor/aporepressor structural transition, the conformational shift is not caused by the binding or loss of the L-tryptophan ligand. We conclude that while L-tryptophan binding is essential for forming a specific complex with trp operator DNA, the corepressor ligand does not lock the repressor into a single conformation that is complementary to the operator. This flexibility may be required by the various binding modes proposed for trp repressor in its search for and adherence to its three different operator sites.

Combined procedure of distance geometry and restrained molecular dynamics techniques for protein structure determination from nuclear magnetic resonance data: application to the DNA binding domain of lac repressor from Escherichia coli

The technique of two-dimensional nuclear magnetic resonance (2D-NMR) has recently assumed an active role in obtaining information on structures of polypeptides, small proteins, sugars, and DNA fragments in solution. In order to generate spatial structures from the atom-atom distance information obtained by the NMR method, different procedures have been developed. Here we introduce a combined procedure of distance geometry (DG) and molecular dynamics (MD) calculations for generating 3D structures that are consistent with the NMR data set and have reasonable internal energies. We report the application of the combined procedure on the lac repressor DNA binding domain (headpiece) using a set of 169 NOE and 17 "hydrogen bond" distance constraints. Eight of ten structures generated by the distance geometry algorithm were refined within 10 ps MD simulation time to structures with low internal energies that satisfied the distance constraints. Although the combination of DG and MD was designed to combine the good sampling properties of the DG algorithm with an efficient method of lowering the internal energy of the molecule, we found that the MD algorithm contributes significantly to the sampling as well.

1H-NMR sequential assignments and cation-binding studies of spinach plastocyanin

The essentially complete assignment of the 1H-NMR spectrum of the Cu(i) form of spinach plastocyanin has been achieved using two-dimensional NMR techniques and sequence-specific resonance assignment procedures. A variety of pH and temperature conditions was utilised to overcome the problems of resonance overlap in the spectrum, degeneracy of C alpha H and solvent H2O chemical shifts, and cross-saturation of labile NH resonances. A qualitative analysis of the long-range nuclear Overhauser effects observed indicates that the backbone fold of spinach plastocyanin is very similar to that of poplar plastocyanin, whose structure has been solved by X-ray crystallography and differs in 22 of its 99 amino acid residues. The assignments provide a basis for further investigations into the structural and ion- and protein-binding properties of plastocyanin in solution.

Sequence-specific assignments and their use in NMR studies of DNA structure

There has been a surge of recent interest, reflected by a sharp increase in the number of publications, in the area of high-resolution nuclear magnetic resonance (NMR) studies of DNA. The goal of many of these studies is to monitor the structure of biologically important DNA sequences directly in solution; the impetus for such studies was the realization, from early single-crystal X-ray structures, that nearest-neighbor context effects are a major determinant of local structure in short double-helical DNAs (Dickerson & Drew, 1981; Dickerson, 1983). Thus, instead of the previously assumed regular averaged structure of the double helix derived from fibre diffraction analysis, the more interesting concept emerged that specific sequence-dependent distortions from ‘classical’ DNA structure might be responsible for the recognition of such sequences by a variety of ligands such as repressors, polymerases, drugs, etc.

1H nuclear magnetic resonance studies of an integral membrane protein: subunit c of the F1F0 ATP synthase

The membrane-traversing subunit c parallel from the F0 part of the ATP synthase molecule has been studied in chloroform/methanol by high-resolution 1H n.m.r. Various one-dimensional and two-dimensional techniques have been used for assignment purposes, some NOE connectivities were established and some 3JHN alpha coupling constants were measured from spin--echo experiments. The effects of varying pH, solvent composition, lanthanide concentration and temperature have been investigated. Evidence is presented that the molecule has extensive alpha-helical segments, and the hairpin structure suggested by other groups is supported by our n.m.r. data. Only one ionizable group, assigned to the C-terminal carboxyl, is observed to titrate in the pH range 2 to 10; so the conserved residue, Asp61, which binds dicyclohexylcarbodiimide, presumably has (at least in this solvent system) an abnormally high pK value.

Nuclear magnetic resonance study of the solution structure of alpha 1-purothionin. Sequential resonance assignment, secondary structure and low resolution tertiary structure

The solution structure of the 45-residue plant protein, alpha 1-purothionin, is investigated by nuclear magnetic resonance (n.m.r.) spectroscopy. Using a combination of two-dimensional n.m.r. techniques to demonstrate through-bond and through-space (less than 5 A) connectivities, the 1H n.m.r. spectrum of alpha 1-purothionin is assigned in a sequential manner. The secondary structure elements are then delineated on the basis of a qualitative interpretation of short-range nuclear Overhauser effects (NOE) involving the NH, C alpha H and C beta H protons. There are two helices extending from residues 10 to 19 and 23 to 28, two short beta-strands from residues 3 to 5 and 31 to 34 which form a mini anti-parallel beta-sheet, and five turns. In addition, a number of long-range NOE connectivities are assigned and a low resolution tertiary structure is proposed.

Complex of lac repressor headpiece with a 14 base-pair lac operator fragment studied by two-dimensional nuclear magnetic resonance

Two-dimensional nuclear Overhauser enhancement spectra are presented of the complex of lac repressor headpiece with a 14 base-pair lac operator fragment. Analysis of nuclear Overhauser enhancements observed between protein and DNA shows that the second helix of the headpiece ("the recognition helix") binds in the major groove of DNA as has been suggested, but that the orientation of this helix is approximately 180 degrees different from the proposed models.

Residue-specific assignments of resonances in the 1H nuclear magnetic resonance spectrum of ribosomal protein E-L30 by systematic application of two-dimensional Fourier transform nuclear magnetic resonance methods

A two-dimensional Fourier transform nuclear magnetic resonance study of the ribosomal protein E-L30 is reported. Five two-dimensional techniques, namely: nuclear magnetic resonance J-resolved spectroscopy, correlated spectroscopy, double quantum spectroscopy, relayed coherence transfer and nuclear Overhauser enhancement spectroscopy were used. Qualitative inspection of the spectra obtained by these techniques provided evidence that the E-L30 molecule has a well-defined structure in solution. This analysis indicated that, despite the fact that the protein is stable only at moderate temperatures and neutral pH, a structural analysis of the molecule would be feasible. A detailed analysis of the spectra permitted unambiguous discrimination between the spin systems of different amino acids, resulting in residue-specific resonance assignments. We were able to assign all resonances of all six threonine, four valine, five alanine, two histidine, two serine, one phenylalanine, one asparagine and one aspartic acid residue of E-L30. Complete resonance assignment was obtained for two glycine residues. Partial assignments became available for all six isoleucine, three glycine and one glutamine residue. These results form a sound basis for the structure determination of the protein described in the accompanying paper.

Application of molecular dynamics with interproton distance restraints to three-dimensional protein structure determination. A model study of crambin

The applicability of restrained molecular dynamics for the determination of three-dimensional protein structures on the basis of short interproton distances (less than 4 A) that can be realistically determined from nuclear magnetic resonance measurements in solution is assessed. The model system used is the 1.2 A resolution crystal structure of the 46 residue protein crambin, from which a set of 240 approximate distance restraints, divided into three ranges (2.5 +/- 0.5, 3.0+0.5(-1.0) and 4 +/- 1 A), is derived. This interproton distance set comprises 159 short-range ([i-j] less than or equal to 5) and 56 ([i-j] greater than 5) long-range inter-residue distances and 25 intra-residue distances. Restrained molecular dynamics are carried out using a number of different protocols starting from two initial structures: a completely extended beta-strand; and an extended structure with two alpha-helices in the same positions as in the crystal structure (residues 7 to 19, and 23 to 30) and all other residues in the form of extended beta-strands. The root-mean-square (r.m.s.) atomic differences between these two initial structures and the crystal structure are 43 A and 23 A, respectively. It is shown that, provided protocols are used that permit the secondary structure elements to form at least partially prior to folding into a tertiary structure, convergence to the correct final structure, both globally and locally, is achieved. The r.m.s. atomic differences between the converged restrained dynamics structures and the crystal structure range from 1.5 to 2.2 A for the backbone atoms and from 2.0 to 2.8 A for all atoms. The r.m.s. atomic difference between the X-ray structure and the structure obtained by first averaging the co-ordinates of the converged restrained dynamics structures is even smaller: 1.0 A for the backbone atoms and 1.6 A for all atoms. These results provide a measure with which to judge future experimental results on proteins whose crystal structures are unknown. In addition, from an examination of the dynamics trajectories, it is shown that the convergence pathways followed by the various simulations are different.

Solution structure of human growth hormone releasing factor. Combined use of circular dichroism and nuclear magnetic resonance spectroscopy

The solution structures of two human growth hormone releasing factor analogues, 27Leu45Gly-hGHRF(1-45)OH and 27Nle-hGHRF(1-29)NH2, are investigated by means of circular dichroism and nuclear magnetic resonance spectroscopy. Using circular dichroism spectroscopy, it is shown that both peptides adopt ordered structures at low concentrations of trifluoroethanol (approximately 30%). Quantitative analysis of the circular dichroism spectra indicates that the same number of residues, approximately 23 to 25, are in a helical state in both peptides. Using two-dimensional nuclear magnetic resonance methods all proton resonances of the 27Nle-hGHRF(1-29)NH2 fragment are assigned and its secondary structure is determined from a qualitative interpretation of the nuclear Overhauser enhancement data. Two distinctive regions of alpha-helix are present extending from residues 6 to 13 and 16 to 29.

Nuclear magnetic resonance study of the globular domain of chicken histone H5: resonance assignment and secondary structure

A 1H NMR study of the 79-residue globular domain of chicken erythrocyte histone H5 (GH5) is presented. Using a combination of two-dimensional NMR techniques to demonstrate through-bond and through-space (less than 5 A) connectivities, the resonances of GH5 are assigned in a sequential manner. From a qualitative interpretation of the short-range nuclear Overhauser effects (NOEs) involving the NH and C alpha H protons, it is shown that GH5 has four alpha-helices. The approximate spatial relationship of three of these four helices relative to each other is deduced from the observation of a number of long-range NOEs. The peptide chain outside the helices appears to have little regular secondary structure and no NOEs characteristic of beta-sheets are apparent.

Structure refinement of oligonucleotides by molecular dynamics with nuclear Overhauser effect interproton distance restraints: application to 5' d(C-G-T-A-C-G)2

The solution structure of the self-complementary DNA hexamer 5' d(C-G-T-A-C-G)2 is refined by restrained molecular dynamics in which 192 interproton distances, determined from pre-steady-state nuclear Overhauser enhancement measurements, are incorporated into the total energy of the system in the form of effective potentials. First the method is tested by applying an idealized set of distance restraints taken from classical B-DNA to a simulation starting off from A-DNA and vice versa. It is shown that in both cases the expected transition between A- and B-DNA occurs. Second, a set of restrained molecular dynamics calculations is carried out starting from both A- and B-DNA with the experimental interproton distances for 5' d(C-G-T-A-C-G)2 as restraints. Convergence to the same B-type structure is achieved with the interproton distances equal to the measured values within experimental error. The root-mean-square atomic difference between the two average restrained dynamics structures (less than 1 A) is approximately the same as the root-mean-square fluctuations of the atoms.

Specific assignment of resonances in the 1H nuclear magnetic resonance spectrum of the polypeptide cardiac stimulant anthopleurin-A

The specific assignment of resonances in the 300-MHz 1H nuclear magnetic resonance (NMR) spectrum of anthopleurin-A, a polypeptide cardiac stimulant from the sea anemone Anthopleura xanthogrammica, is described. Assignments have been made using two-dimensional NMR techniques, in particular the method of sequential assignments, where through-bond and through-space connectivities to the peptide backbone NH resonances are used to identify the spin systems of residues adjacent in the amino acid sequence. Complete assignments have been made of the resonances from 33 residues out of a total of 49, and partial assignments of a further 3. The resonances from several of the remaining residues have been identified but not yet specifically assigned. A complicating factor in making these assignments is the conformational heterogeneity exhibited by anthopleurin-A in solution. The resonances from a number of amino acid residues in the minor conformer have also been assigned. These assignments contribute towards identification of the origin of this heterogeneity, and permit some preliminary conclusions to be drawn regarding the secondary structure of the polypeptide.

Solution conformation of a heptadecapeptide comprising the DNA binding helix F of the cyclic AMP receptor protein of Escherichia coli. Combined use of 1H nuclear magnetic resonance and restrained molecular dynamics

A nuclear magnetic resonance study on a heptadecamer (17-mer) peptide comprising the DNA binding helix F of the cyclic AMP receptor protein of Escherichia coli is presented under solution conditions (viz. 40% (v/v) trifluorethanol) where it adopts an ordered helical structure as judged by circular dichroism. Using a combination of two-dimensional nuclear magnetic resonance techniques, complete resonance assignments are obtained in a sequential manner. From the two-dimensional nuclear Overhauser enhancement spectra, a set of 87 approximate distance restraints is derived and used as the basis for three-dimensional structure determination with a restrained molecular dynamics algorithm in which the interproton distances are incorporated into the total energy function of the system in the form of an additional effective potential term. The convergence properties of this approach are tested by starting from three different initial structures, namely an alpha-helix, a beta-strand and a 3-10 helix. In all three cases, convergence to an alpha-helical structure is achieved with a root mean square difference of less than 3 A for all atoms and less than 2 A for the backbone atoms.

Two-dimensional 1H NMR investigation of ribonuclease A and ribonuclease-A--pyrimidine-nucleotide complexes

Ribonuclease A was studied by two-dimensional 1H NMR spectroscopy. 10 out of 12 alanine and 9 out of 10 threonine spin systems as well as all valine [9], leucine [2] and isoleucine [3] spin systems were identified from the correlated spectroscopy (COSY) and relayed coherence transfer spectroscopy (RCT). Sequence-specific assignments were obtained from nuclear Overhauser effect spectra for proton resonances of 21 amino acid moieties. 2' and 3'-pyrimidine-nucleotide-RNase-A complexes were also investigated by two-dimensional NMR. We were able to monitor structural changes in the active center, the vicinity of the active center and in regions far from the catalytic region. Chemical shift changes of resonances of protons near Thr-45 reflected the binding of the same moiety. This in turn is also dependent on the position of the nucleotide phosphate group. Binding of 2' nucleotides led to characteristic changes in protein regions not affected by the binding of 3' nucleotides. These results are interpreted in terms of structural differences between the 2' and 3'-nucleotide-RNase-A complexes; the structure of the complex of the native 3' nucleotide inhibitor being more closely related to that of the free protein.

Preliminary identification of the secondary structure of the trp repressor from Escherichia coli

The probable secondary structure content of the trp repressor from Escherichia coli has been inferred from NMR and circular dichroic measurements; the results are compared with those of prediction algorithms. 70% of the amide protons have exchange rate constants orders of magnitude smaller than the intrinsic rate constants, identifying them as participating in hydrogen bonds. The exchange rate constants fall into two distinct classes, one having half-lives of 20 min and the other more than 24 h. The latter class, consisting of 50% of all amide protons, indicates a stable core. The exchange data are consistent with circular dichroism and predictions that suggest that about 55% of the peptides from alpha helix, and 20% form beta sheets and turns. The NMR spectrum further indicates that there is little beta sheet, suggesting that the secondary structure class is alpha.

NMR studies of the trp repressor from Escherichia coli. Characterisation and assignments of residue types

High-resolution proton nuclear magnetic resonance spectra of the trp repressor of Escherichia coli under various conditions are reported and analysed. The spectrum of the denatured state agrees with that predicted from the amino acid composition, with the exception of the two histidine residues, which have different chemical shifts although they titrate normally. The spectrum of the native protein shows the presence of extensive secondary and tertiary structure. Using information from chemical shifts, numbers of protons, titration behaviour, homonuclear chemical-shift-correlated spectroscopy and nuclear Overhauser enhancement correlated spectroscopy, most of the aromatic protons have been assigned to residue type. Further, about 30% of the aliphatic protons have been assigned to residue type by two-dimensional spectroscopy. Nuclear Overhauser enhancements establish that high-field methyl groups belonging to a valine residue lie directly over an aromatic ring.

Hydrogen exchange of individual amide protons in the E. coli lac repressor DNA-binding domain: a nuclear magnetic resonance study

Proton exchange in lac repressor headpiece was studied by COSY and 2D NOE spectroscopy. The exchange rates of amide protons, stabilized by the hydrogen bonds of the three alpha-helices of the headpiece, could be determined quantitatively. The exchange rates in these helices showed repetitive patterns of about three to four residues. A correlation with the position of the amide proton in the interior or the exterior of the alpha-helix of the protein was found. The exchange data strongly support the validity of the three-dimensional structure, as determined recently (Kaptein, R. et al., J. Mol. Biol. 182, 179-182 (1985)).

Two-dimensional 1H NMR studies of cytochrome c

Two-dimensional nuclear magnetic resonance techniques were used to assign the NH, C alpha H, and C beta H protons of over 60 of the 104 amino acid residues in the 1H NMR spectrum of horse ferrocytochrome c. The majority of these amino acids were completely assigned. Assignments were based on the analysis of two-dimensional J-correlated (COSY), nuclear Overhauser effect (NOESY), and relayed COSY spectra and on comparisons of the J-correlated spectra of various cytochrome c species. Spin diffusion is not a problem with monomeric proteins the size of cytochrome c. Here these advances are illustrated with data that lead to the assignment of the heme-associated residues cysteine-14 and tryptophan-59, the axial ligands methionine-80 and histidine-18, the entire N-terminal helix, and several other amino acid spin systems. With these approaches, structure, structure change, the internal dynamics of cytochrome c, and the interaction of these with function are being studied, especially by observation of the hydrogen exchange behavior of essentially all the H-bonded amides and some side chain protons in both the reduced and oxidized proteins.