Refinement of protein structures in explicit solvent

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

Simulated annealing with restrained molecular dynamics using a flexible restraint potential: theory and evaluation with simulated NMR constraints.

A new functional representation of NMR-derived distance constraints, the flexible restraint potential, has been implemented in the program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168) for molecular structure generation. In addition, flat-bottomed restraint potentials for representing dihedral angle and vicinal scalar coupling constraints have been introduced into CONGEN. An effective simulated annealing (SA) protocol that combines both weight annealing and temperature annealing is described. Calculations have been performed using ideal simulated NMR constraints, in order to evaluate the use of restrained molecular dynamics (MD) with these target functions as implemented in CONGEN. In this benchmark study, internuclear distance, dihedral angle, and vicinal coupling constant constraints were calculated from the energy-minimized X-ray crystal structure of the 46-amino acid polypeptide crambin (ICRN). Three-dimensional structures of crambin that satisfy these simulated NMR constraints were generated using restrained MD and SA. Polypeptide structures with extended backbone and side-chain conformations were used as starting conformations. Dynamical annealing calculations using extended starting conformations and assignments of initial velocities taken randomly from a Maxwellian distribution were found to adequately sample the conformational space consistent with the constraints. These calculations also show that loosened internuclear constraints can allow molecules to overcome local minima in the search for a global minimum with respect to both the NMR-derived constraints and conformational energy. This protocol and the modified version of the CONGEN program described here are shown to be reliable and robust, and are applicable generally for protein structure determination by dynamical simulated annealing using NMR data.     

Test of circular dichroism (CD) methods for crambin and CD-assisted secondary structure prediction of its homologous toxins

Methods that analyze protein circular dichroism (CD) spectra for fractions of secondary structure are evaluated for the plant protein crambin, which has a known high-resolution crystal structure. In addition, a two-step secondary structure prediction scheme is presented and used for the toxins homologous to crambin, shown by others to have secondary structures similar to crambin. The test of CD spectral analysis methods with the protein crambin employed two computer programs and several CD basis sets. Crambin's crystal structure, known to 0.945A resolution (Hendrickson, W.A., Teeter, M.M. Nature 290:107-113, 1981), allows accurate evaluation of results. Analysis with the protein spectra basis sets (Provencher, S.W., Gl?ckner, J. Biochemistry 20:33-37, 1981) as modified (Manavalan, P., Johnson, W.C., Jr. Anal. Biochem. 167:76-85, 1987) agreed most closely with crambin's crystal structure. This method was then applied to the CD spectra of the membrane-active toxins homologous to crambin (alpha 1- and beta-purothionin, phoratoxin A and B, and viscotoxin A3 and B). The new program SEQ (pronounced "seek") was developed to assign the secondary structure along the protein chain in a hierarchical fashion and applied to the plant toxins. The method constrained the secondary structure fractions to those from CD analysis and combined standard statistical methods with amphipathic helix location. Both CD-arrived secondary structure percentages and sequence assignment indicate that the viscotoxins are structurally most similar to crambin. Purothionin's secondary structure was predicted to be fundamentally similar to crambin's with a difference at the start of the first helix. This assignment agreed with Raman and NMR analyses of purothionin and lends validity to the method presented here. Differences from the NMR in the CD secondary structure fraction analysis for phoratoxin suggest interference in the CD from tryptophan residues.     

Structure of Paramecium tetraurelia calmodulin at 1.8 A resolution.

The crystal structure of calmodulin (CaM; M(r) 16,700, 148 residues) from the ciliated protozoan Paramecium tetraurelia (PCaM) has been determined and refined using 1.8 A resolution area detector data. The crystals are triclinic, space group P1, a = 29.66, b = 53.79, c = 25.49 A, alpha = 92.84, beta = 97.02, and gamma = 88.54 degrees with one molecule in the unit cell. Crystals of the mammalian CaM (MCaM; Babu et al., 1988) and Drosophila CaM (DCaM; Taylor et al., 1991) also belong to the same space group with very similar cell dimensions. All three CaMs have 148 residues, but there are 17 sequence changes between PCaM and MCaM and 16 changes between PCaM and DCaM. The initial difference in the molecular orientation between the PCaM and MCaM crystals was approximately 7 degrees as determined by the rotation function. The reoriented Paramecium model was extensively refitted using omit maps and refined using XPLOR. The R-value for 11,458 reflections with F > 3 sigma is 0.21, and the model consists of protein atoms for residues 4-147, 4 calcium ions, and 71 solvent molecules. The root mean square (rms) deviations in the bond lengths and bond angles in the model from ideal values are 0.016 A and 3 degrees, respectively. The molecular orientation of the final PCaM model differs from MCaM by only 1.7 degrees. The overall Paramecium CaM structure is very similar to the other calmodulin structures with a seven-turn long central helix connecting the two terminal domains, each containing two Ca-binding EF-hand motifs. The rms deviation in the backbone N, Ca, C, and O atoms between PCaM and MCaM is 0.52 A and between PCaM and DCaM is 0.85 A. The long central helix regions differ, where the B-factors are also high, particularly in PCaM and MCaM. Unlike the MCaM structure, with one kink at D80 in the middle of the linker region, and the DCaM structure, with two kinks at K75 and I85, in our PCaM structure there are no kinks in the helix; the distortion appears to be more gradually distributed over the entire helical region, which is bent with an apparent radius of curvature of 74.5(2) A. The different distortions in the central helical region probably arise from its inherent mobility.     

Crystal structure of a complex of HIV-1 protease with a dihydroxyethylene-containing inhibitor: comparisons with molecular modeling.

The structure of a crystal complex of recombinant human immunodeficiency virus type 1 (HIV-1) protease with a peptide-mimetic inhibitor containing a dihydroxyethylene isostere insert replacing the scissile bond has been determined. The inhibitor is Noa-His-Hch psi [CH(OH)CH(OH)]Vam-Ile-Amp (U-75875), and its Ki for inhibition of the HIV-1 protease is < 1.0 nM (Noa = 1-naphthoxyacetyl, Hch = a hydroxy-modified form of cyclohexylalanine, Vam = a hydroxy-modified form of valine, Amp = 2-pyridylmethylamine). The structure of the complex has been refined to a crystallographic R factor of 0.169 at 2.0 A resolution by using restrained least-squares procedures. Root mean square deviations from ideality are 0.02 A and 2.4 degrees, for bond lengths and angles, respectively. The bound inhibitor diastereomer has the R configurations at both of the hydroxyl chiral carbon atoms. One of the diol hydroxyl groups is positioned such that it forms hydrogen bonds with both the active site aspartates, whereas the other interacts with only one of them. Comparison of this X-ray structure with a model-built structure of the inhibitor, published earlier, reveals similar positioning of the backbone atoms and of the side-chain atoms in the P2-P2' region, where the interaction with the protein is strongest. However, the X-ray structure and the model differ considerably in the location of the P3 and P3' end groups, and also in the positioning of the second of the two central hydroxyl groups. Reconstruction of the central portion of the model revealed the source of the hydroxyl discrepancy, which, when corrected, provided a P1-P1' geometry very close to that seen in the X-ray structure.     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

Crystal structure of human annexin I at 2.5 A resolution.

cDNA coding for N-terminally truncated human annexin I, a member of the family of Ca(2+)-dependent phospholipid binding proteins, has been cloned and expressed in Escherichia coli. The expressed protein is biologically active, and has been purified and crystallized in space group P2(1)2(1)2(1) with cell dimensions a = 139.36 A, b = 67.50 A, and c = 42.11 A. The crystal structure has been determined by molecular replacement at 3.0 A resolution using the annexin V core structure as the search model. The average backbone deviation between these two structures is 2.34 A. The structure has been refined to an R-factor of 17.7% at 2.5 A resolution. Six calcium sites have been identified in the annexin I structure. Each is located in the loop region of the helix-loop-helix motif. Two of the six calcium sites in annexin I are not occupied in the annexin V structure. The superpositions of the corresponding loop regions in the four domains show that the calcium binding loops in annexin I can be divided into two classes: type II and type III. Both classes are different from the well-known EF-hand motif (type I).     

Crystal structure of IscS,a cysteine desulfurase from Escherichia coli

IscS is a widely distributed cysteine desulfurase that catalyzes the pyridoxal phosphate-dependent desulfuration of L-cysteine and plays a central role in the delivery of sulfur to a variety of metabolic pathways. We report the crystal structure of Escherichia coli IscS to a resolution of 2.1A. The crystals belong to the space group P2(1)2(1)2(1) and have unit cell dimensions a=73.70A, b=101.97A, c=108.62A (alpha=beta=gamma=90 degrees ). Molecular replacement with the Thermotoga maritima NifS model was used to determine phasing, and the IscS model was refined to an R=20.6% (R(free)=23.6%) with two molecules per asymmetric unit. The structure of E.coli IscS is similar to that of T.maritima NifS with nearly identical secondary structure and an overall backbone r.m.s. difference of 1.4A. However, in contrast to NifS a peptide segment containing the catalytic cysteine residue (Cys328) is partially ordered in the IscS structure. This segment of IscS (residues 323-335) forms a surface loop directed away from the active site pocket. Cys328 is positioned greater than 17A from the pyridoxal phosphate cofactor, suggesting that a large conformational change must occur during catalysis in order for Cys328 to participate in nucleophilic attack of a pyridoxal phosphate-bound cysteine substrate. Modeling suggests that rotation of this loop may allow movement of Cys328 to within approximately 3A of the pyridoxal phosphate cofactor.     

Solution structure of apamin determined by nuclear magnetic resonance and distance geometry

The solution structure of the bee venom neurotoxin apamin has been determined with a distance geometry program using distance constraints derived from NMR. Twenty embedded structures were generated and refined by using the program DSPACE. After error minimization using both conjugate gradient and dynamics algorithms, six structures had very low residual error. Comparisons of these show that the backbone of the peptide is quite well-defined with the largest rms difference between backbone atoms in these structures of 1.34 A. The side chains have far fewer constraints and show greater variability in their positions. The structure derived here is generally consistent with the qualitative model previously described, with most differences occurring in the loop between the beta-turn (residues 2-5) and the C-terminal alpha-helix (residues 9-17). Comparisons are made with previously derived models from NMR data and other methods.     

Crystal structure of alpha-hordothionin at 1.9 Angstrom resolution

Crystal structure of ubiquitous toxin from barley alpha-hordothionin (alpha-HT) has been determined at 1.9A resolution by X-ray crystallography. The primary sequence as well as the NMR solution structure of alpha-HT firmly established that alpha-HT belongs to a family of membrane active plant toxins-thionins. Since alpha-HT crystallized in a space group (P4(1)2(1)2) that is different from the space group (I422) of previously determined alpha(1)- and beta-purothionins, and visocotoxin A3, therefore, it provided independent information on protein-protein interactions that may be relevant to the toxin mechanism. The structure of alpha-HT not only confirms overall architectural features (crambin fold) but also provides an additional confirmation of the role for crucial solute molecules, that were postulated to be directly involved in the mechanism of toxicity for thionins.     

Exploration of disorder in protein structures by X-ray restrained molecular dynamics

Conformational disorder in crystal structures of ribonuclease-A and crambin is studied by including two independent structures in least-squares optimizations against X-ray data. The optimizations are carried out by X-ray restrained molecular dynamics (simulated annealing refinement) and by conventional least-squares optimization. Starting from two identical structures, the optimizations against X-ray data lead to significant deviations between the two, with rms backbone displacements of 0.45 A for refinement of ribonuclease at 1.53 A resolution, and 0.31 A for crambin at 0.945 A. More than 15 independent X-ray restrained molecular dynamics runs have been carried out for ribonuclease, and the displacements between the resulting structures are highly reproducible for most atoms. These include residues with two or more conformations with significant dihedral angle differences and alternative hydrogen bonding, as well as groups of residues that undergo displacements that are suggestive of rigid-body librations. The crystallographic R-values obtained are approximately 13%, as compared to 15.3% for a comparable refinement with a single structure. Least-squares optimization without an intervening restrained molecular dynamics stage is sufficient to reproduce most of the observed displacements. Similar results are obtained for crambin, where the higher resolution of the X-ray data allows for refinement of unconstrained individual anisotropic temperature factors. These are shown to be correlated with the displacements in the two-structure refinements.     

NMR structure of free RGS4 reveals an induced conformational change upon binding Galpha

Heterotrimeric guanine nucleotide-binding proteins (G-proteins) are transducers in many cellular transmembrane signaling systems where regulators of G-protein signaling (RGS) act as attenuators of the G-protein signal cascade by binding to the Galpha subunit of G-proteins (G(i)(alpha)(1)) and increasing the rate of GTP hydrolysis. The high-resolution solution structure of free RGS4 has been determined using two-dimensional and three-dimensional heteronuclear NMR spectroscopy. A total of 30 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 2871 experimental NMR restraints. The atomic rms distribution about the mean coordinate positions for residues 5-134 for the 30 structures is 0.47 +/- 0.05 A for the backbone atoms, 0. 86 +/- 0.05 A for all atoms, and 0.56 +/- 0.04 A for all atoms excluding disordered side chains. The NMR solution structure of free RGS4 suggests a significant conformational change upon binding G(i)(alpha)(1) as evident by the backbone atomic rms difference of 1. 94 A between the free and bound forms of RGS4. The underlying cause of this structural change is a perturbation in the secondary structure elements in the vicinity of the G(i)(alpha)(1) binding site. A kink in the helix between residues K116-Y119 is more pronounced in the RGS4-G(i)(alpha)(1) X-ray structure relative to the free RGS4 NMR structure, resulting in a reorganization of the packing of the N-terminal and C-terminal helices. The presence of the helical disruption in the RGS4-G(i)(alpha)(1) X-ray structure allows for the formation of a hydrogen-bonding network within the binding pocket for G(i)(alpha)(1) on RGS4, where RGS4 residues D117, S118, and R121 interact with residue T182 from G(i)(alpha)(1). The binding pocket for G(i)(alpha)(1) on RGS4 is larger and more accessible in the free RGS4 NMR structure and does not present the preformed binding site observed in the RGS4-G(i)(alpha)(1) X-ray structure. This observation implies that the successful complex formation between RGS4 and G(i)(alpha)(1) is dependent on both the formation of the bound RGS4 conformation and the proper orientation of T182 from G(i)(alpha)(1). The observed changes for the free RGS4 NMR structure suggest a mechanism for its selectivity for the Galpha-GTP-Mg(2+) complex and a means to facilitate the GTPase cycle.     

The NMR solution structure of the pheromone Er-1 from the ciliated protozoan Euplotes raikovi.

The 3-dimensional structure of the pheromone Er-1 isolated from the ciliated protozoan Euplotes raikovi has been determined in aqueous solution by 1H NMR spectroscopy. The structure of this 40-residue protein was calculated with the distance geometry program DIANA on the basis of 503 upper distance constraints derived from nuclear Overhauser effects and 77 dihedral angle constraints derived from spin-spin coupling constants, and refined by restrained energy minimization with the program OPAL. The Er-1 solution structure is represented by a group of 20 conformers with an average RMS deviation relative to the mean structure of 0.55 A for the backbone atoms N, C alpha, and C', and 0.93 A for all heavy atoms of the complete polypeptide chain, residues 1-40. The molecular architecture is dominated by an up-down-up bundle of 3 alpha-helices formed by residues 2-9, 12-19, and 24-33. Although this core part coincides closely with the previously determined structure of the homologous pheromone Er-10, the C-terminal peptide segment adopts a novel conformation. This is of interest in view of previous suggestions, based on sequence comparisons, that this molecular region may be important for the different specificity of receptor recognition by different pheromones.     

Energy minimization for tertiary structure prediction of homologous proteins: alpha 1-purothionin and viscotoxin A3 models from crambin

Homologous proteins may fold into similar three-dimensional structures. Spectroscopic evidence suggests this is true for the cereal grain thionins, the mistletoe toxins, and for crambin, three classes of plant proteins. We have combined primary sequence homology and energy minimization to predict the structures alpha 1-purothionin (from Durum wheat) and viscotoxin A3 (from Viscum album, European mistletoe) from the high resolution (0.945 A) crystal structure of crambin (from Crambe abyssinica). Our predictions will be verifiable because we have diffraction-quality crystals of alpha 1-purothionin whose structure we are have predicted. The potential energy minimizations for each protein were performed both with and without harmonic constraints to its initial backbone to explore the existence of local minima for the predicted proteins. Crambin was run as a control to examine the effects of the potential energy minimization on a protein with a well-known structure. Only alpha 1-purothionin which has one fewer residue in a turn region shows a significant difference for the two minimization paths. The results of these predictions suggest that alpha 1-purothionin and viscotoxin are amphipathic proteins, and this character may relate to the mechanism of action for these proteins. Both are mildly membrane-active and their amphipatic character is well suited for interaction with a lipid bilayer.     

Solution structure of alpha t alpha, a helical hairpin peptide of de novo design.

alpha t alpha is a 38-residue peptide designed to adopt a helical hairpin conformation in solution (Fezoui Y, Weaver DL Osterhout JJ, 1995, Protein Sci 4:286-295). A previous study of the carboxylate form of alpha t alpha by CD and two-dimensional NMR indicated that the peptide was highly helical and that the helices associated in approximately the intended orientation (Fezoui Y, Weaver DL, Osterhout JJ, 1994, Proc Natl Acad Sci USA 91:3675-3679). Here, the solution structure of alpha t alpha as determined by two-dimensional NMR is reported. A total of 266 experimentally derived distance restraints and 20 dihedral angle restraints derived from J-couplings were used. One-hundred initial structures were generated by distance geometry and refined by dynamical simulated annealing. Twenty-three of the lowest-energy structures consistent with the experimental restraints were analyzed. The results presented here show that alpha t alpha is comprised of two associating helices connected by a turn region.     

Monte carlo simulations of protein folding. II. Application to protein A,ROP, and crambin

The hierarchy of lattice Monte Carlo models described in the accompanying paper (Kolinski, A., Skolnick, J. Monte Carlo simulations of protein folding. I. Lattice model and interaction scheme. Proteins 18:338–352, 1994) is applied to the simulation of protein folding and the prediction of 3-dimensional structure. Using sequence information alone, three proteins have been successfully folded: the B domain of staphylococcal protein A, a 120 residue, monomeric version of ROP dimer, and crambin. Starting from a random expanded conformation, the model proteins fold along relatively well-defined folding pathways. These involve a collection of early intermediates, which are followed by the final (and rate-determining) transition from compact intermediates closely resembling the molten globule state to the native-like state. The predicted structures are rather unique, with native-like packing of the side chains. The accuracy of the predicted native conformations is better than those obtained in previous folding simulations. The best (but by no means atypical) folds of protein A have a coordinate rms of 2.25 Å from the native Cα trace, and the best coordinate rms from crambin is 3.18 Å. For ROP monomer, the lowest coordinate rms from equivalent Cαs of ROP dimer is 3.65 Å. Thus, for two simple helical proteins and a small α/β protein, the ability to predict protein structure from sequence has been demonstrated. © 1994 John Wiley & Sons, Inc.     

Revised 2.3 A structure of porcine pepsin: evidence for a flexible subdomain

A revised three-dimensional crystal structure of ethanol-inhibited porcine pepsin refined to an R-factor of 0.171 at 2.3 A resolution is presented and compared to the refined structures of the fungal aspartic proteinases: penicillopepsin, rhizopuspepsin, and endothiapepsin. Pepsin is composed of two nearly equal N and C domains related by an intra dyad. The overall polypeptide fold and active site structures are homologous for pepsin and the fungal enzymes. The weak inhibition of pepsin by ethanol can be explained by the presence of one or more ethanol molecules, in the vicinity of the active site carboxylates, which slightly alter the hydrogen-bonding network and which may compete with substrate binding in the active site. Structural superposition analysis showed that the N domains aligned better than the C-domains for pepsin and the fungal aspartic proteinases: 107-140 C alpha pairs aligned to 0.72-0.85 A rms for the N domains; 64-95 C alpha pairs aligned to 0.78-1.03 A rms for the C domains. The major structural difference between pepsin and the fungal enzymes concerns a newly described subdomain whose conformation varies markedly among these enzyme structures. The subdomain in pepsin comprises nearly 100 residues and is composed of two contiguous segments within the C domain (residues 192-212 and 223-299). the subdomain is connected, or "hinged," to a mixed beta-sheet that forms one of the structurally invariant, active site psi-loops. Relative subdomain displacements as large as a 21.0 degrees rotation and a 5.9 A translation were observed among the different enzymes. There is some suggestion in pepsin that the subdomain may be flexible and perhaps plays a structural role in mediating substrate binding, determining the substrate specificity, or in the activation of the zymogen.     

Three-dimensional solution structure of a curaremimetic toxin from Naja nigricollis venom: a proton NMR and molecular modeling study.

The solution conformation of toxin alpha from Naja nigricollis (61 amino acids and four disulfides), a snake toxin which specifically blocks the activity of the nicotinic acetylcholine receptor (AcChoR), has been determined using nuclear magnetic resonance spectroscopy and molecular modeling. The solution structures were calculated using 409 distance and 73 dihedral angle restraints. The average atomic rms deviation between the eight refined structures and the mean structure is approximately 0.5 A for the backbone atoms. The overall folding of toxin alpha consists of three major loops which are stabilized by three disulfide bridges and one short C terminal loop stabilized by a fourth disulfide bridge. All the disulfides are grouped in the same region of the molecule, forming a highly constrained structure from which the loops protrude. As predicted, this structure appears to be very similar to the 1.4-A resolution crystal structure of another snake neurotoxin, namely, erabutoxin b from Laticauda semifasciata. The atomic rms deviation for the backbone atoms between the solution and crystal structures is approximately 1.7 A. The minor differences which are observed between the two structures are partly related to the deletion of one residue from the chain of toxin alpha. It is notable that, although the two toxins differ from each other by 16 amino acid substitutions, their side chains have an essentially similar spatial organization. However, most of the side chains which constitute the presumed AcChoR binding site for the curaremimetic toxins are poorly resolved in toxin alpha.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structural analysis of zinc-substituted cytochrome c

Zinc-substituted cytochrome c has been widely used in studies of protein-protein interactions and photo-induced electron transfer reactions between proteins. However, the coordination geometry of zinc in zinc-substituted cyt c has not yet been determined; two different opinions about the coordination have been reached. Here the solution structures of zinc-substituted cytochrome c that might be five-coordinated and six-coordinated have been refined separately by using (1)H NMR spectroscopy, and the zinc coordination geometry was determined just by NOE distance constraints. Structural analysis of the energy-minimized average solution structures of both the pentacoordinated and hexacoordinated geometries indicate that that zinc in zinc-substituted cyt c should be bound to both His18 and Met80, which means that the zinc is six-coordinated. RMSD values of the family of 25 six-coordinated structures from the average structure are 0.66+/-0.13 A and 1.09+/-0.16 A for the backbone and all heavy atoms, respectively. A statistical analysis of the structure indicates its satisfactory quality. Comparison of the solution structure of the six-coordinated energy-minimized average structure of zinc-substituted cytochrome c with the solution structure of reduced cytochrome c reveals that for the overall folding the secondary structure elements are very close. The availability of the structure provides for a better understanding of the protein-protein complex and for electron transfer processes between Zn cyt c and other metalloproteins.     

Direct nuclear Overhauser effect refinement of crambin from two-dimensional nmr data using a slow-cooling annealing protocol

The solution structure of crambin has been refined using a direct nuclear Overhauser effect (NOE) simulation approach (DINOSAUR) following a slow-cooling simulated annealing protocol starting from eight previously obtained nmr and the x-ray structures of crambin. Theoretical NOE intensities calculated with inclusion of local motions were directly compared to the experimental nmr data and forces were derived using a simple first-order approximation for the calculation of the NOE gradient. A dynamic assignment procedure was applied for the peaks involving unassigned diastereotopic proton pairs or equivalent aromatic protons. With this approach, R factors could be minimized in a reasonable simulation time to low values (around 0.26) while deviations from ideal bond lengths and angles are still acceptable. The improvement in R factors is accompanied by an improvement of the precision of the structures, the rms deviations (rmsd; from the average) calculated on the ensemble of nine structures decreasing from 0.65 to 0.55 Å for backbone atoms and from 1.0 to 0.85 Å for all heavy atoms. The solution structure is significantly different from the x-ray structure with rmsd for all atoms of 1.35 Å compared to 0.85 Å between solution structures. The largest differences are found for residues Thr-21 and Pro-22 in the loop region between the two α-helices and for the side chain of Tyr-29. © 1994 John Wiley & Sons, Inc.