Symmetrical stabilization of bound Ca2+ ions in a cooperative pair of EF-hands through hydrogen bonding of coordinating water molecules in calbindin D(9k).

Water molecules are found to complete the Ca2+ coordination sphere when a protein fails to provide enough ligating oxygens. Hydrogen bonding of these water molecules to the protein backbone or side chains may contribute favorably to the Ca2+ affinity, as suggested in an earlier study of two calbindin D(9k) mutants [E60D and E60Q; Linse et al. (1994) Biochemistry 33, 12478-12486]. To investigate the generality of this conclusion, another side chain, Gln 22, which hydrogen bonds to a Ca2+-coordinating water molecule in calbindin D(9k), was mutated. Two calbindin D(9k) mutants, (Q22E+P43M) and (Q22N+P43M), were constructed to examine the interaction between Gln 22 and the water molecule in the C-terminal calcium binding site II. Shortening of the side chain, as in (Q22N+P43M), reduces the affinity of binding two calcium ions by a factor of 18 at low ionic strength, whereas introduction of a negative charge, as in (Q22E+P43M), leads to a 12-fold reduction. In 0.15 M KCl, a 7-fold reduction in affinity was observed for both mutants. The cooperativity of Ca2+ binding increases for (Q22E+P43M), while it decreases for (Q22N+P43M). The rates of Ca2+ dissociation are 5.5-fold higher for the double mutants than for P43M at low ionic strength. For both mutants, reduced strength of hydrogen bonding to calcium-coordinating water molecules is a likely explanation for the observed effects on Ca2+ affinity and dissociation. In the apo forms, the (Q22E+P43M) mutant has lower stability toward urea denaturation than (Q22N+P43M) and P43M. 2D (1)H NMR and crystallographic experiments suggest that the structure of (Q22E+P43M) and (Q22N+P43M) is unchanged relative to P43M, except for local perturbations in the loop regions.

Prediction of protein secondary structure at 80% accuracy.

Secondary structure prediction involving up to 800 neural network predictions has been developed, by use of novel methods such as output expansion and a unique balloting procedure. An overall performance of 77.2%-80.2% (77.9%-80.6% mean per-chain) for three-state (helix, strand, coil) prediction was obtained when evaluated on a commonly used set of 126 protein chains. The method uses profiles made by position-specific scoring matrices as input, while at the output level it predicts on three consecutive residues simultaneously. The predictions arise from tenfold, cross validated training and testing of 1032 protein sequences, using a scheme with primary structure neural networks followed by structure filtering neural networks. With respect to blind prediction, this work is preliminary and awaits evaluation by CASP4.

When size is important. Accommodation of magnesium in a calcium binding regulatory domain.

The accommodation of Mg2+ in the N-terminal domain of calmodulin was followed through amide 1H and 15N chemical shifts and line widths in heteronuclear single-quantum coherence spectroscopy NMR spectra. Mg2+ binds sequentially to the two Ca2+-binding loops in this domain, with affinities such that nearly half of the loops would be occupied by Mg2+ in resting eukaryotic cells. Mg2+ binding seems to occur without ligation to the residue in the 12th loop position, previously proven largely responsible for the major rearrangements induced by binding of the larger Ca2+. Consequently, smaller Mg2+-induced structural changes are indicated throughout the protein. The two Ca2+-binding loops have different Mg2+ binding characteristics. Ligands in the N-terminal loop I are better positioned for cation binding, resulting in higher affinity and slower binding kinetics compared with the C-terminal loop II (koff = 380 +/- 40 s-1 compared with approximately 10,000 s-1 at 25 degreesC). The Mg2+-saturated loop II undergoes conformational exchange on the 100-microseconds time scale. Available data suggest that this exchange occurs between a conformation providing a ligand geometry optimized for Mg2+ binding and a conformation more similar to that of the empty loop.

Distributed torsion angle grid search in high dimensions: a systematic approach to NMR structure determination.

Two complementary approaches for systematic search in torsion angle space are described for the generation of all conformations of polypeptides which satisfy experimental NMR restraints, hard-sphere van der Waals radii, and rigid covalent geometry. The first procedure is based on a recursive, tree search algorithm for the examination of linear chains of torsion angles, and uses a novel treatment to propagate the search results to neighboring regions so that the structural consequences of the restraints are fully realized. The second procedure is based on a binary combination of torsion vector spaces for connected submolecules, and produces intermediate results in Cartesian space for a more robust restraint analysis. Restraints for NMR applications include bounds on torsion angles and internuclear distances, including relational and degenerate restraints involving equivalent and nonstereoassigned protons. To illustrate these methods, conformation search results are given for the tetrapeptide APGA restrained to an idealized beta-turn conformation, an alanine octapeptide restrained to a right-handled helical conformation, and the structured region of the peptide SYPFDV.

Solution structure of the N-terminal EGF-like domain from human factor VII.

Blood coagulation is initiated by Ca(2+)-dependent binding of coagulation factor VIIa (FVIIa) to its cofactor, tissue factor (TF). The TF:FVIIa complex activates factors IX and X, ultimately leading to the formation of thrombin and the coagulation of blood. FVII consists of an N-terminal gamma-carboxyglutamic-acid-containing (Gla) domain followed by two epidermal growth factor (EGF) like domains, the first of which can bind one Ca2+ ion (Kd approximately 150 microM) and a C-terminal serine protease domain. Using 1H nuclear magnetic resonance spectroscopy, we have determined the solution structure of a synthetic N-terminal EGF-like domain (EGF1) of human FVII (residues 45-85) in the absence of Ca2+. A comparison of this structure of apo EGF1 with the Ca(2+)-bound EGF1 in the complex of FVIIa and TF [Banner, D. W., et al. (1996) Nature 380, 41-46] suggests that the structural changes in the EGF1 domain upon Ca2+ binding are minor and are concentrated near the Ca(2+)-binding site, which is facing away from the TF interaction surface. Amino acid side chains that are crucial for the binding of FVII to TF show a similar conformation in both structures and are therefore unlikely to directly influence the Ca(2+)-dependent binding of FVII to TF. As Ca2+ binding to EGF1 does not lead to a conformational change in the residues constituting the interaction surface for binding to TF, our results are consistent with the idea that the altered orientation between the Gla and EGF1 domains that result from Ca2+ binding is responsible for the increased affinity of FVII/FVIIa for TF.

Thermodynamics of a reverse turn motif. Solvent effects and side-chain packing.

The linear pentapeptide, Ala-Tyr-cis-Pro-Tyr-Asp-NMA (AYPYD) is known to have a significant population of type VI turn conformers in aqueous solvent. We have carried out theoretical studies of the conformational energetics of this peptide using a potential of mean force (PMF) consisting of the AMBER/OPLS empirical potential energy function, a macroscopic electrostatic model of polar solvation, and a surface area-based model of non-polar solvation. Conformers were taken from molecular dynamics simulations reported elsewhere, or generated by a random search method reported here. The chain entropy of folding was calculated by a systematic search of accessible dihedral angle space. The intra-peptide component was found to strongly favor folding and was nearly cancelled by the polar solvation term which disfavored folding. The non-polar solvation term had little effect. Fluctuations about the average value of the PMF were small and in accord with estimates from a simple harmonic model. When applied to conformers generated by a random search, the PMF selected a conformer close to the NMR-determined structure as the lowest energy conformer. The conformer with the second-lowest energy was extended, but was found to fold rapidly to the turn state in a subsequent molecular dynamics study, and may be an important state on the folding-unfolding pathway. Averages of the PMF were combined with the entropy estimates to provide an estimate of the free energy of folding that is in reasonable agreement with experimental results. In terms of the interplay between backbone electrostatic interactions and the packing of apolar side-chains, this peptide provides a model for the energetics of protein folding, and therefore makes a useful test case for calculations.

Low resolution solution structure of the Bacillus subtilis glucose permease IIA domain derived from heteronuclear three-dimensional NMR spectroscopy.

A low resolution solution structure of the IIA domain of the Bacillus subtilis phosphoenolpyruvate-sugar phosphotransferase system (PTS) glucose permease has been determined using 945 inter-residue and 724 intra-residue distance constraints derived from three-dimensional 15N and 13C edited NOESY spectra. A total of 15 structures was generated using distance geometry calculations. The protein is comprised of 13 beta-strands forming an antiparallel beta-barrel. The average backbone atomic RMS deviation about the average distance geometry structure for the beta-sheet residues is 1.1 A. The conformations of the loop regions between the beta-strands are less well determined. Backbone distance constraints obtained during the process of sequential assignment were insufficient to correctly calculate the polypeptide fold.

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.

Computational methods for determining protein structures from NMR data.

The general procedures by which solution structures of proteins may be deduced from distance and angular constraints derived from NMR are reviewed, with an emphasis on practical aspects of the calculations. In addition, novel methods based on chemical shift calculations and on quantitative fits to nuclear Overhauser effect intensities are presented; these should provide improved understanding of the limits of our ability to simulate complex spectra, and may permit higher precision structures to be determined.

Three-dimensional solution structure of the reduced form of Escherichia coli thioredoxin determined by nuclear magnetic resonance spectroscopy.

The three-dimensional solution structure of reduced (dithiol) thioredoxin from Escherichia coli has been determined with distance and dihedral angle constraints obtained from 1H NMR spectroscopy. Reduced thioredoxin has a well-defined global fold consisting of a central five-strand beta-sheet and three long helices. The beta-strands are packed in the sheet in the order beta 1 beta 3 beta 2 beta 4 beta 5, with beta 1, beta 3, and beta 2 parallel and beta 2, beta 4, and beta 5 arranged in an antiparallel fashion. Two of the helices connect strands of the beta-sheet: alpha 1 between beta 1 and beta 2 and alpha 2 between beta 2 and beta 3. Strands beta 4 and beta 5 are connected by a short loop that contains a beta-bulge. Strands beta 3 and beta 4 are connected by a long loop that contains a series of turn-like or 3(10) helical structures. The active site Cys-Gly-Pro-Cys sequence forms a protruding loop between strand beta 2 and helix alpha 2. The structure is very similar overall to that of oxidized (disulfide) thioredoxin obtained from X-ray crystal structure analysis but differs in the local conformation of the active site loop. The distance between the sulfurs of Cys 32 and Cys 35 increases from 2.05 A in the disulfide bridge to 6.8 +/- 0.6 A in the dithiol of reduced thioredoxin, as a result of a rotation of the side chain of Cys 35 and a significant change in the position of Pro 34. This conformational change has important implications for the mechanism of thioredoxin as a protein disulfide oxidoreductase.

Three-dimensional solution structure of a single zinc finger DNA-binding domain.

The three-dimensional solution structure of a zinc finger nucleic acid binding motif has been determined by nuclear magnetic resonance (NMR) spectroscopy. Spectra of a synthetic peptide corresponding to a single zinc finger from the Xenopus protein Xfin yielded distance and dihedral angle constraints that were used to generate structures from distance geometry and restrained molecular dynamics calculations. The zinc finger is an independently folded domain with a compact globular structure in which the zinc atom is bound by two cysteine and two histidine ligands. The polypeptide backbone fold consists of a well-defined helix, starting as alpha and ending as 3(10) helix, packed against two beta strands that are arranged in a hairpin structure. A high density of basic and polar amino acid side chains on the exposed face of the helix are probably involved in DNA binding.

Three-dimensional solution structure of plastocyanin from the green alga Scenedesmus obliquus.

The solution conformation of plastocyanin from the green alga Scenedesmus obliquus has been determined from distance and dihedral angle constraints derived by nuclear magnetic resonance (NMR) spectroscopy. Structures were generated with distance geometry and restrained molecular dynamics calculations. A novel molecular replacement method was also used with the same NMR constraints to generate solution structures of S. obliquus plastocyanin from the x-ray structure of the homologous poplar protein. Scenedesmus obliquus plastocyanin in solution adopts a beta-barrel structure. The backbone conformation is well defined and is similar overall to that of poplar plastocyanin in the crystalline state. The distinctive acidic region of the higher plant plastocyanins, which functions as a binding site for electron transfer proteins and inorganic complexes, differs in both shape and charge in S. obliquus plastocyanin.