High-resolution three-dimensional structure of reduced recombinant human thioredoxin in solution

The solution structure of recombinant human thioredoxin (105 residues) has been determined by nuclear magnetic resonance (NMR) spectroscopy combined with hybrid distance geometry-dynamical simulated annealing calculations. Approximate interproton distance restraints were derived from nuclear Overhauser effect (NOE) measurements. In addition, a large number of stereospecific assignments for beta-methylene protons and torsion angle restraints for phi, psi, and chi 1 were obtained by using a conformational grid search on the basis of the intraresidue and sequential NOE data in conjunction with 3JHN alpha and 3J alpha beta coupling constants. The structure calculations were based on 1983 approximate interproton distance restraints, 52 hydrogen-bonding restraints for 26 hydrogen bonds, and 98 phi, 71 psi, and 72 chi 1 torsion angle restraints. The 33 final simulated annealing structures obtained had an average atomic rms distribution of the individual structures about the mean coordinate positions of 0.40 +/- 0.06 A for the backbone atoms and 0.78 +/- 0.05 A for all atoms. The solution structure of human thioredoxin consists of a five-stranded beta-sheet surrounded by four alpha-helices, with an active site protrusion containing the two redox-active cysteines. The overall structure is similar to the crystal and NMR structures of oxidized [Katti, S. K., LeMaster, D. M., & Eklund, H. (1990) J. Mol. Biol. 212, 167-184] and reduced [Dyson, J. H., Gippert, G. P., Case, D. A., Holmgren, A., & Wright, P. (1990) Biochemistry 29, 4129-4136] Escherichia coli thioredoxin, respectively, despite the moderate 25% amino acid sequence homology. Several differences, however, can be noted. The human alpha 1 helix is a full turn longer than the corresponding helix in E. coli thioredoxin and is characterized by a more regular helical geometry. The helix labeled alpha 3 in human thioredoxin has its counterpart in the 3(10) helix of the E. coli protein and is also longer in the human protein. In contrast to these structural differences, the conformation of the active site loop in both proteins is very similar, reflecting the perfect sequence identity for a stretch of eight amino acid residues around the redox-active cysteines.     

Determination of the spatial structure of insectotoxin 15A from Buthus erpeus by (1)H-NMR spectroscopy data]

The solution structure of insectotoxin 15A (35 residues) from scorpion Buthus eupeus was determined on the basis of 386 interproton distance restraints 12 hydrogen-bonding restraints and 113 dihedral angle restraints derived from 1H NMR experiments. A group of 20 structures was calculated with the distance geometry program DIANA followed by the restrained energy minimization with the program CHARMM. The atomic RMS distribution about the mean coordinate position is 0.64 +/- 0.11 A for the backbone atoms and 1.35 +/- 0.20 A for all atoms. The structure contains an alpha-helix (residues 10-20) and a three-stranded antiparallel beta-sheet (residues 2-5, 24-28 and 29-33). A pairing of the eight cysteine residues of insectotoxin 15A was established basing on NMR data. Three disulfide bridges (residues 2-19, 16-31 and 20-33) connect the alpha-helix with the beta-sheet, and the fourth one (5-26) joins beta-strands together. The spatial fold of secondary structure elements (the alpha-helix and the beta-sheet) of the insectotoxin 15A is very similar to those of the other short and long scorpion toxins in spite of a low (about 20%) sequence homology.     

Determination of the three-dimensional solution structure of the antihypertensive and antiviral protein BDS-I from the sea anemone Anemonia sulcata: a study using nuclear magnetic resonance and hybrid distance geometry-dynamical simulated annealing

The three-dimensional solution structure of the antihypertensive and antiviral protein BDS-I from the sea anemone Anemonia sulcata has been determined on the basis of 489 interproton and 24 hydrogen-bonding distance restraints supplemented by 23 phi backbone and 21 chi 1 side-chain torsion angle restraints derived from nuclear magnetic resonance (NMR) measurements. A total of 42 structures is calculated by a hybrid metric matrix distance geometry-dynamical simulated annealing approach. Both the backbone and side-chain atom positions are well defined. The average atomic rms difference between the 42 individual SA structures and the mean structure obtained by averaging their coordinates is 0.67 +/- 0.12 A for the backbone atoms and 0.90 +/- 0.17 A for all atoms. The core of the protein is formed by a triple-stranded antiparallel beta-sheet composed of residues 14-16 (strand 1), 30-34 (strand 2), and 37-41 (strand 3) with an additional mini-antiparallel beta-sheet at the N-terminus (residues 6-9). The first and second strands of the triple-stranded antiparallel beta-sheet are connected by a long exposed loop (residues 17-30). A number of side-chain interactions are discussed in light of the structure.     

Solution structure and dynamics of a designed monomeric variant of the lambda Cro repressor.

The solution structure of a monomeric variant of the lambda Cro repressor has been determined by multidimensional NMR. Cro K56[DGEVK] differs from wild-type Cro by the insertion of five amino acids at the center of the dimer interface. 1H and 15N resonances for 70 of the 71 residues have been assigned. Thirty-two structures were calculated by hybrid distance geometry/simulated annealing methods using 463 NOE-distance restraints, 26 hydrogen-bond, and 39 dihedral-angle restraints. The root-mean-square deviation (RMSD) from the average structure for atoms in residues 3-60 is 1.03 +/- 0.44 A for the peptide backbone and 1.6 +/- 0.73 A for all nonhydrogen atoms. The overall structure conforms very well to the original design. Although the five inserted residues form a beta hairpin as expected, this engineered turn as well as other turns in the structure are not well defined by the NMR data. Dynamics studies of backbone amides reveal T1/T2 ratios of residues in the alpha2-alpha3, beta2-beta3, and engineered turn that are reflective of chemical exchange or internal motion. The solution structure and dynamics are discussed in light of the conformational variation that has been observed in other Cro structures, and the importance of flexibility in DNA recognition.     

Determination of a high precision structure of a novel protein, Linum usitatissimum trypsin inhibitor (LUTI), using computer-aided assignment of NOESY cross-peaks

The solution structure of a novel 69 residue proteinase inhibitor, Linum usitatissimum trypsin inhibitor (LUTI), was determined using a method based on computer aided assignment of nuclear Overhauser enhancement spectroscopy (NOESY) data. The approach applied uses the program NOAH/DYANA for automatic assignment of NOESY cross-peaks. Calculations were carried out using two unassigned NOESY peak lists and a set of determined dihedral angle restraints. In addition, hydrogen bonds involving amide protons were identified during calculations using geometrical criteria and values of HN temperature coefficients. Stereospecific assignment of beta-methylene protons was carried out using a standard procedure based on nuclear Overhauser enhancement intensities and 3J(alpha)(beta) coupling constants. Further stereospecific assignment of methylene protons and diastereotopic methyl groups were established upon structure-based method available in the program GLOMSA and chemical shift calculations. The applied algorithm allowed us to assign 1968 out of 2164 peaks (91%) derived from NOESY spectra recorded in H2O and 2H2O. The final experimental data input consisted of 1609 interproton distance restraints, 88 restraints for 44 hydrogen bonds, 63 torsion angle restraints and 32 stereospecifically assigned methylene proton pairs and methyl groups. The algorithm allowed the calculation of a high precision protein structure without the laborious manual assignment of NOESY cross-peaks. For the 20 best conformers selected out of 40 refined ones in the program CNS, the calculated average pairwise rmsd values for residues 3 to 69 were 0.38 A (backbone atoms) and 1.02 A (all heavy atoms). The three-dimensional LUTI structure consists of a mixed parallel and antiparallel beta-sheet, a single alpha-helix and shows the fold of the potato 1 family of proteinase inhibitors. Compared to known structures of the family, LUTI contains Arg and Trp residues at positions P6' and P8', respectively, instead of two Arg residues, involved in the proteinase binding loop stabilization. A consequence of the ArgTrp substitution at P8' is a slightly more compact conformation of the loop relative to the protein core.     

The determination of the three-dimensional structure of barley serine proteinase inhibitor 2 by nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The solution structure of the 64 residue structured domain (residues 20-83) of barley serine proteinase inhibitor 2 (BSPI-2) is determined on the basis of 403 interproton distance, 34 phi backbone torsion angle and 26 hydrogen bonding restraints derived from n.m.r. measurements. A total of 11 converged structures were computed using a metric matrix distance geometry algorithm and refined by restrained molecular dynamics. The average rms difference between the final 11 structures and the mean structure obtained by averaging their coordinates is 1.4 +/- 0.2 A for the backbone atoms and 2.1 +/- 0.1 A for all atoms. The overall structure, which is almost identical to that found by X-ray crystallography, is disc shaped and consists of a central four component mixed parallel and antiparallel beta-sheet flanked by a 13 residue alpha-helix on one side and the reactive site loop on the other.     

Determination of the three-dimensional solution structure of barnase using nuclear magnetic resonance spectroscopy

The solution conformation of the ribonuclease barnase has been determined by using 1H nuclear magnetic resonance (NMR) spectroscopy. The 20 structures were calculated by using 853 interproton distance restraints obtained from analyses of two-dimensional nuclear Overhauser spectra, 72 phi and 53 chi 1 torsion angle restraints, and 17 hydrogen-bond distance restraints. The calculated structures contain two alpha-helices (residues 6-18 and 26-34) and a five-stranded antiparallel beta-sheet (residues 50-55, 70-75, 85-91, 94-101, and 105-108). The core of the protein is formed by the packing of one of the alpha-helices (residues 6-18) onto the beta-sheet. The average RMS deviation between the calculated structures and the mean structure is 1.11 A for the backbone atoms and 1.75 A for all atoms. The protein is least well-defined in the N-terminal region and in three large loops. When these regions are excluded, the average RMS deviation between the calculated structures and the mean structure for residues 5-34, 50-56, 71-76, 85-109 is 0.62 A for the backbone atoms and 1.0 A for all atoms. The NMR-derived structure has been compared with the crystal structure of barnase [Mauguen et al. (1982) Nature (London) 297, 162-164].     

Three-dimensional solution structure of beta cryptogein, a beta elicitin secreted by a phytopathogenic fungus Phytophthora cryptogea.

Cryptogein belongs to a new family of 10-kDa proteins called elicitins. Elicitins are necrotic and signaling proteins secreted by Phytophthora spp. responsible for the incompatible reaction and systemic hypersensitive-like necroses of diverse plant species leading to resistance against fungal or bacterial plant pathogens. The solution structure of beta cryptogein from Phytophthora cryptogea fungus was determined by using multidimensional heteronuclear nuclear magnetic resonance spectroscopy. A set of 18 structures was calculated using 1360 NOE-derived distance restraints and 40 dihedral angle restraints obtained from 3JHNH alpha couplings. The RMS deviation from the mean structure is 0.87 +/- 0.14 A for backbone atoms and 1.34 +/- 0.14 A for all the non-hydrogen atoms of residues 2 to 98. The structure of beta cryptogein reveals a novel protein fold, with five helices and a double-stranded beta-sheet facing an omega-loop. One edge of the beta-sheet and the adjacent face of the omega-loop form a hydrophobic cavity. This cavity made of highly conserved residues represents a plausible binding site. Residue 13, which has been identified from directed mutagenesis and natural sequence comparison studies as a key amino acid involved in the differential control of necrosis, is surface exposed and could contribute to the binding to a ligand or a receptor. The solution structure is close to the X-ray structure, with slight differences lightly due to the crystal packing.     

Solution-state structure by NMR of zinc-substituted rubredoxin from the marine hyperthermophilic archaebacterium Pyrococcus furiosus.

The three-dimensional solution-state structure is reported for the zinc-substituted form of rubredoxin (Rd) from the marine hyperthermophilic archaebacterium Pyrococcus furiosus, an organism that grows optimally at 100 degrees C. Structures were generated with DSPACE by a hybrid distance geometry (DG)-based simulated annealing (SA) approach that employed 403 nuclear Overhauser effect (NOE)-derived interproton distance restraints, including 67 interresidue, 124 sequential (i-j = 1), 75 medium-range (i-j = 2-5), and 137 long-range (i-j > 5) restraints. All lower interproton distance bounds were set at the sum of the van Der Waals radii (1.8 A), and upper bounds of 2.7 A, 3.3 A, and 5.0 A were employed to represent qualitatively observed strong, medium, and weak NOE cross peak intensities, respectively. Twenty-three backbone-backbone, six backbone-sulfur (Cys), two backbone-side chain, and two side chain-side chain hydrogen bond restraints were include for structure refinement, yielding a total of 436 nonbonded restraints, which averages to > 16 restraints per residue. A total of 10 structures generated from random atom positions and 30 structures generated by molecular replacement using the backbone coordinates of Clostridium pasteurianum Rd converged to a common conformation, with the average penalty (= sum of the square of the distance bounds violations; +/- standard deviation) of 0.024 +/- 0.003 A2 and a maximum total penalty of 0.035 A2. Superposition of the backbone atoms (C, C alpha, N) of residues A1-L51 for all 40 structures afforded an average pairwise root mean square (rms) deviation value (+/- SD) of 0.42 +/- 0.07 A. Superposition of all heavy atoms for residues A1-L51, including those of structurally undefined external side chains, afforded an average pairwise rms deviation of 0.72 +/- 0.08 A. Qualitative comparison of back-calculated and experimental two-dimensional NOESY spectra indicate that the DG/SA structures are consistent with the experimental spectra. The global folding of P. furiosus Zn(Rd) is remarkably similar to the folding observed by X-ray crystallography for native Rd from the mesophilic organism C. pasteurianum, with the average rms deviation value for backbone atoms of residues A1-L51 of P. furiosus Zn(Rd) superposed with respect to residues K2-V52 of C. pasteurianum Rd of 0.77 +/- 0.06 A. The conformations of aromatic residues that compose the hydrophobic cores of the two proteins are also similar. However, P. furiosus Rd contains several unique structural elements, including at least four additional hydrogen bonds and three potential electrostatic interactions.(ABSTRACT TRUNCATED AT 400 WORDS)     

Determination of the three-dimensional solution structure of the C-terminal domain of cellobiohydrolase I from Trichoderma reesei. A study using nuclear magnetic resonance and hybrid distance geometry-dynamical simulated annealing

The solution structure of a synthetic 36-residue polypeptide comprising the C-terminal cellulose binding domain of cellobiohydrolase I (CT-CBH I) from Trichoderma reesei was investigated by nuclear magnetic resonance (NMR) spectroscopy. The 1H NMR spectrum was completely assigned in a sequential manner by two-dimensional NMR techniques. A large number of stereospecific assignments for beta-methylene protons, as well as ranges for the phi, psi, and chi 1 torsion angles, were obtained on the basis of sequential and intraresidue nuclear Overhauser enhancement (NOE) and coupling constant data in combination with a conformational data base search. The structure calculations were carried out in an iterative manner by using the hybrid distance geometry-dynamical simulated annealing method. This involved computing a series of initial structures from a subset of the experimental data in order to resolve ambiguities in the assignments of some NOE cross-peaks arising from chemical shift degeneracy. Additionally, this permitted us to extend the stereospecific assignments to the alpha-methylene protons of glycine using information on phi torsion angles derived from the initial structure calculations. The final experimental data set consisted of 554 interproton distance restraints, 24 restraints for 12 hydrogen bonds, and 33 phi, 24 psi, and 25 chi 1 torsion angle restraints. CT-CBH I has two disulfide bridges whose pairing was previously unknown. Analysis of structures calculated with all three possible combinations of disulfide bonds, as well as without disulfide bonds, indicated that the correct disulfide bridge pairing was 8-25 and 19-35. Forty-one structures were computed with the 8-25 and 19-35 disulfide bridges, and the average atomic rms difference between the individual structures and the mean structure obtained by averaging their coordinates was 0.33 +/- 0.04 A for the backbone atoms and 0.52 +/- 0.06 A for all atoms. The protein has a wedgelike shape with an amphiphilic character, one face being predominantly hydrophilic and the other mainly hydrophobic. The principal element of secondary structure is made up of an irregular triple-stranded antiparallel beta-sheet composed of residues 5-9 (beta 1), 24-28 (beta 2), and 33-36 (beta 3) in which strand beta 3 is hydrogen bonded to the other two strands.(ABSTRACT TRUNCATED AT 400 WORDS)     

Solution structure of a sweet protein single-chain monellin determined by nuclear magnetic resonance and dynamical simulated annealing calculations

Single-chain monellin (SCM), which is an engineered 94-residue polypeptide, has proven to be as sweet as native two-chain monellin. SCM is more stable than the native monellin for both heat and acidic environments. Data from gel filtration HPLC and NMR indicate that the SCM exists as a monomer in aqueous solution. The solution structure of SCM has been determined by nuclear magnetic resonance (NMR) spectroscopy and dynamical simulated annealing calculations. A stable alpha-helix spanning residues Phe11-Ile26 and an antiparallel beta-sheet formed by residues 2-5, 36-38, 41-47, 54-64, 69-75, and 83-88 have been identified. The sheet was well defined by backbone-backbone NOEs, and the corresponding beta-strands were further confirmed by hydrogen bond networks based on amide hydrogen exchange data. Strands beta2 and beta3 are connected by a small bulge comprising residues Ile38-Cys41. A total of 993 distance and 56 dihedral angle restraints were used for simulated annealing calculations. The final simulated annealing structures (k) converged well with a root-mean-square deviation (rmsd) between backbone atoms of 0.49 A for secondary structural regions and 0.70 A for backbone atoms excluding two loop regions. The average restraint energy-minimized (REM) structure exhibited root-mean-square deviations of 1.19 A for backbone atoms and 0.85 A for backbone atoms excluding two loop regions with respect to 20 k structures. The solution structure of SCM revealed that the long alpha-helix was folded into the concave side of a six-stranded antiparallel beta-sheet. The side chains of Tyr63 and Asp66 which are common to all sweet peptides showed an opposite orientation relative to H1 helix, and they were all solvent-exposed. Residues at the proposed dimeric interface in the X-ray structure were observed to be mostly solvent-exposed and demonstrated high degrees of flexibility.     

Solution structure of a lipid transfer protein extracted from rice seeds. Comparison with homologous proteins.

Nuclear magnetic resonance (NMR) spectroscopy was used to determine the three dimensional structure of rice nonspecific lipid transfer protein (ns-LTP), a 91 amino acid residue protein belonging to the broad family of plant ns-LTP. Sequence specific assignment was obtained for all but three HN backbone 1H resonances and for more than 95% of the 1H side-chain resonances using a combination of 1H 2D NOESY; TOCSY and COSY experiments at 293 K. The structure was calculated on the basis of four disulfide bridge restraints, 1259 distance constraints derived from 1H-1H Overhauser effects, 72 phi angle restraints and 32 hydrogen-bond restraints. The final solution structure involves four helices (H1: Cys3-Arg18, H2: Ala25-Ala37, H3: Thr41-Ala54 and H4: Ala66-Cys73) followed by a long C-terminal tail (T) with no observable regular structure. N-capping residues (Thr2, Ser24, Thr40), whose side-chain oxygen atoms are involved in hydrogen bonds with i + 3 amide proton additionally stabilize the N termini of the first three helices. The fourth helix involving Pro residues display a mixture of alpha and 3(10) conformation. The rms deviation of 14 final structures with respect to the average structure is 1.14 +/- 0.16 A for all heavy atoms (C, N, O and S) and 0.72 +/- 0.01 A for the backbone atoms. The global fold of rice ns-LTP is close to the previously published structures of wheat, barley and maize ns-LTPs exhibiting nearly identical pattern of the numerous sequence specific interactions. As reported previously for different four-helix topology proteins, hydrophobic, hydrogen bonding and electrostatic mechanisms of fold stabilization were found for the rice ns-LTP. The sequential alignment of 36 ns-LTP primary structures strongly suggests that there is a uniform pattern of specific long-range interactions (in terms of sequence), which stabilize the fold of all plant ns-LTPs.     

Three-dimensional structure of interleukin 8 in solution

The solution structure of the interleukin 8 (IL-8) dimer has been solved by nuclear magnetic resonance (NMR) spectroscopy and hybrid distance geometry-dynamical simulated annealing calculations. The structure determination is based on a total of 1880 experimental distance restraints (of which 82 are intersubunit) and 362 torsion angle restraints (comprising phi, psi, and chi 1 torsion angles). A total of 30 simulated annealing structures were calculated, and the atomic rms distribution about the mean coordinate positions (excluding residues 1-5 of each subunit) is 0.41 +/- 0.08 A for the backbone atoms and 0.90 +/- 0.08 A for all atoms. The three-dimensional solution structure of the IL-8 dimer reveals a structural motif in which two symmetry-related antiparallel alpha-helices, approximately 24 A long and separated by about 14 A, lie on top of a six-stranded antiparallel beta-sheet platform derived from two three-stranded Greek keys, one from each monomer unit. The general architecture is similar to that of the alpha 1/alpha 2 domains of the human class I histocompatibility antigen HLA-A2. It is suggested that the two alpha-helices form the binding site for the cellular receptor and that the specificity of IL-8, as well as that of a number of related proteins involved in cell-specific chemotaxis, mediation of cell growth, and the inflammatory response, is achieved by the distinct distribution of charged and polar residues at the surface of the helices.     

Solution structure and dynamics of an open beta-sheet, glycolytic enzyme, monomeric 23.7 kDa phosphoglycerate mutase from Schizosaccharomyces pombe

The structure and backbone dynamics of a double labelled (15N,13C) monomeric, 23.7 kD phosphoglycerate mutase (PGAM) from Schizosaccharomyces pombe have been investigated in solution using NMR spectroscopy. A set of 3125 NOE-derived distance restraints, 148 restraints representing inferred hydrogen bonds and 149 values of (3)J(HNHalpha) were used in the structure calculation. The mean rmsd from the average structure for all backbone atoms from residues 6-205 in the best 21 calculated structures was 0.59 A. The core of the enzyme includes an open, twisted, six-stranded beta-sheet flanked by four alpha-helices and a short 3(10)-helix. An additional smaller domain contains two short antiparallel beta-strands and a further pair of alpha-helices. The C(alpha) atoms of the S. pombe PGAM may be superimposed on their equivalents in one of the four identical subunits of Saccharomyces cerevisiae PGAM with an rmsd of 1.34 A (0.92 A if only the beta-sheet is considered). Small differences between the two structures are attributable partly to the deletion in the S. pombe sequence of a 25 residue loop involved in stabilising the S. cerevisiae tetramer. Analysis of 15N relaxation parameters indicates that PGAM tumbles isotropically with a rotational correlation time of 8.7 ns and displays a range of dynamic features. Of 178 residues analysed, only 77 could be fitted without invoking terms for fast internal motion or chemical exchange, and out of the remainder, 77 required a chemical exchange term. Significantly, 46 of the slowly exchanging (milli- to microsecond) residues lie in helices, and these account for two-thirds of all analysed helix residues. On the contrary, only one beta-sheet residue required an exchange term. In contrast to other analyses of backbone dynamics reported previously, residues in slow exchange appeared to correlate with architectural features of the enzyme rather than congregating close to ligand binding sites.     

Low resolution structure of interleukin-1 beta in solution derived from 1H-15N heteronuclear three-dimensional nuclear magnetic resonance spectroscopy

A low resolution solution structure of the cytokine interleukin-1 beta, a 153 residue protein of molecular weight 17,400, has been determined on the basis of 446 nuclear Overhauser effect (NOE) derived approximate interproton distance restraints involving solely NH, C alpha H and C beta H protons, supplemented by 90 distance restraints for 45 hydrogen bonds, and 79 phi torsion angle restraints. With the exception of 27 C alpha H-C alpha H NOEs, all the NOEs were assigned from a three-dimensional 1H-1H NOE 15N-1H heteronuclear multiple quantum coherence (HMQC) spectrum. The torsion angle restraints were obtained from accurate 3JHN alpha coupling constants measured from a HMQC-J spectrum, while the hydrogen bonds were derived from a qualitative analysis of the NOE, coupling constant and amide exchange data. A total of 20 simulated annealing (SA) structures was computed using the hybrid distance geometry-dynamical simulated annealing method. The solution structure of IL-1 beta comprises 12 beta-strands arranged in three pseudo-symmetrical topological units (each consisting of 5 anti-parallel beta-strands), joined by turns, short loops and long loops. The core of the structure, which is made up of the 12 beta-strands, together with the turns joining strands I and II, strands VIII and IX and strands X and XI, is well determined with a backbone atomic root-mean-square (r.m.s.) distribution about the mean co-ordinate positions of 1.2(+/- 0.1) A. The loop conformations, on the other hand, are poorly determined by the current data. A comparison of the core of the low resolution solution structure of IL-1 beta with that of the X-ray structure indicates that they are similar, with a backbone atomic r.m.s. difference of only 1.5 A between the co-ordinates of the restrained minimized mean of the SA structures and the X-ray structure.     

High-resolution solution structure of reduced French bean plastocyanin and comparison with the crystal structure of poplar plastocyanin

The three-dimensional solution structure of reduced (CuI) plastocyanin from French bean leaves has been determined by distance geometry and restrained molecular dynamics methods using constraints obtained from 1H n.m.r. (nuclear magnetic resonance) spectroscopy. A total of 1244 experimental constraints were used, including 1120 distance constraints, 103 dihedral angle constraints and 21 hydrogen bond constraints. Stereospecific assignments were made for 26 methylene groups and the methyls of 11 valines. Additional constraints on copper co-ordination were included in the restrained dynamics calculations. The structures are well defined with average atomic root-mean-square deviations from the mean of 0.45 A for all backbone heavy atoms and 1.08 A for side-chain heavy atoms. French bean plastocyanin adopts a beta-sandwich structure in solution that is similar to the X-ray structure of reduced poplar plastocyanin; the average atomic root-mean-square difference between 16 n.m.r. structures and the X-ray structure is 0.76 A for all backbone heavy atoms. The conformations of the side-chains that constitute the hydrophobic core of French bean plastocyanin are very well defined. Of 47 conserved residues that populate a single chi 1 angle in solution, 43 have the same rotamer in the X-ray structure. Many surface side-chains adopt highly preferred conformations in solution, although the 3J alpha beta coupling constants often indicate some degree of conformational averaging. Some surface side-chains are disordered in both the solution and crystal structures of plastocyanin. There is a striking correlation between measures of side-chain disorder in solution and side-chain temperature factors in the X-ray structure. Side-chains that form a distinctive acidic surface region, believed to be important in binding other electron transfer proteins, appear to be disordered. Fifty backbone amide protons form hydrogen bonds to carbonyls in more than 60% of the n.m.r. structures; 45 of these amide protons exchange slowly with solvent deuterons. Ten hydrogen bonds are formed between side-chain and backbone atoms, eight of which are correlated with decreased proton exchange. Of the 60 hydrogen bonds formed in French bean plastocyanin, 56 occur in the X-ray structure of the poplar protein; two of the missing hydrogen bonds are absent as a result of mutations. It appears that molecular dynamics refinement of highly constrained n.m.r. structures allows accurate prediction of the pattern of hydrogen bonding.     

NMR structure of cysteinyl-phosphorylated enzyme IIB of the N,N'-diacetylchitobiose-specific phosphoenolpyruvate-dependent phosphotransferase system of Escherichia coli

The determination by NMR of the solution structure of the phosphorylated enzyme IIB (P-IIB(Chb)) of the N,N'-diacetylchitobiose-specific phosphoenolpyruvate-dependent phosphotransferase system of Escherichia coli is presented. Most of the backbone and side-chain resonances were assigned using a variety of mostly heteronuclear NMR experiments. The remaining resonances were assigned with the help of the structure calculations.NOE-derived distance restraints were used in distance geometry calculations followed by molecular dynamics and simulated annealing protocols. In addition, combinations of ambiguous restraints were used to resolve ambiguities in the NOE assignments. By combining sets of ambiguous and unambiguous restraints into new ambiguous restraints, an error function was constructed that was less sensitive to information loss caused by assignment uncertainties. The final set of structures had a pairwise rmsd of 0.59 A and 1.16 A for the heavy atoms of the backbone and side-chains, respectively.Comparing the P-IIB(Chb) solution structure with the previously determined NMR and X-ray structures of the wild-type and the Cys10Ser mutant shows that significant differences between the structures are limited to the active-site region. The phosphoryl group at the active-site cysteine residue is surrounded by a loop formed by residues 10 through 16. NOE and chemical shift data suggest that the phosphoryl group makes hydrogen bonds with the backbone amide protons of residues 12 and 15. The binding mode of the phosphoryl group is very similar to that of the protein tyrosine phosphatases. The differences observed are in accordance with the presumption that IIB(Chb) has to be more resistant to hydrolysis than the protein tyrosine phosphatases. We propose a proton relay network by which a transfer occurs between the cysteine SH proton and the solvent via the hydroxyl group of Thr16.     

Solution structure of neurotoxin I from the sea anemone Stichodactyla helianthus. A nuclear magnetic resonance, distance geometry, and restrained molecular dynamics study

The three-dimensional structure of the sea anemone polypeptide Stichodactyla helianthus neurotoxin I in aqueous solution has been determined using distance geometry and restrained molecular dynamics simulations based on NMR data acquired at 500 MHz. A set of 470 nuclear Overhauser enhancement values was measured, of which 216 were used as distance restraints in the structure determination along with 15 dihedral angles derived from coupling constants. After restrained molecular dynamics refinement, the eight structures that best fit the input data form a closely related family. They describe a structure that consists of a core of twisted, four-stranded, antiparallel beta-sheet encompassing residues 1-3, 19-24, 29-34, and 40-47, joined by three loops, two of which are well defined by the NMR data. The third loop, encompassing residues 7-16, is poorly defined by the data and is assumed to undergo conformational averaging in solution. Pairwise root mean square displacement values for the backbone heavy atoms of the eight best structures are 1.3 +/- 0.2A when the poorly defined loop is excluded and 3.6 +/- 1.0A for all backbone atoms. Refinement using restrained molecular dynamics improved the quality of the structures generated by distance geometry calculations with respect to the number of nuclear Overhauser enhancements violated, the size of the total distance violations and the total potential energies of the structures. The family of structures for S. heliathus neurotoxin I is compared with structures of related sea anemone proteins that also bind to the voltage-gated sodium channel.     

The influence of stereospecific assignments on the determination of three-dimensional structures of proteins by nuclear magnetic resonance spectroscopy. Application to the sea anemone protein BDS-I

The influence of the stereospecific assignments of beta-methylene protons and the classification of chi 1 torsion angles on the definition of the three-dimensional structures of proteins determined from NMR data is investigated using the sea anemone protein BDS-I (43 residues) as a model system. Two sets of structures are computed. The first set comprises 42 converged structures (denoted STEREO structures) calculated on the basis of the complete list of restraints derived from the NMR data, consisting of 489 interproton and 24 hydrogen bonding distance restraints, supplemented by 23 phi backbone and 21 chi 1 side chain torsion angle restraints. The second set comprises 31 converged structures (denoted NOSTEREO structures) calculated from a reduced data set in which those restraints arising from stereospecific assignments, and the corresponding chi 1 torsion angle restraints, are explicitly omitted. The results show that the inclusion of the stereospecific restraints leads to a significant improvement in the definition of the structure of BDS-I, both with respect to the backbone and the detailed arrangement of the side chains. Average atomic rms differences between the individual structures and the mean structures for the backbone atoms are 0.67 +/- 0.12 A and 0.93 +/- 0.16 A for the STEREO and NOSTEREO structures, respectively; the corresponding values for all atoms are 0.90 +/- 0.17 A and 1.17 +/- 0.17 A, respectively. In addition, while the overall fold remains unchanged, there is a small but significant atomic displacement between the two sets of structures.     

Solution structure of the PDZ2 domain from human phosphatase hPTP1E and its interactions with C-terminal peptides from the Fas receptor

The solution structure of the second PDZ domain (PDZ2) from human phosphatase hPTP1E has been determined using 2D and 3D heteronuclear NMR experiments. The binding of peptides derived from the C-terminus of the Fas receptor to PDZ2 was studied via changes in backbone peptide and protein resonances. The structure is based on a total of 1387 nonredundant experimental NMR restraints including 1261 interproton distance restraints, 45 backbone hydrogen bonds, and 81 torsion angle restraints. Analysis of 30 lowest-energy structures resulted in rmsd values of 0.41 +/- 0.09 A for backbone atoms (N, Calpha, C') and 1.08 +/- 0.10 A for all heavy atoms, excluding the disordered N- and C-termini. The hPTP1E PDZ2 structure is similar to known PDZ domain structures but contains two unique structural features. In the peptide binding domain, the first glycine of the GLGF motif is replaced by a serine. This serine appears to replace a bound water observed in PDZ crystal structures that hydrogen bonds to the bound peptide's C-terminus. The hPTP1E PDZ2 structure also contains an unusually large loop following strand beta2 and proximal to the peptide binding site. This well-ordered loop folds back against the PDZ domain and contains several residues that undergo large amide chemical shift changes upon peptide binding. Direct observation of peptide resonances demonstrates that as many as six Fas peptide residues interact with the PDZ2 domain.