Solution structure of the functional domain of Paracoccus denitrificans cytochrome c552 in the reduced state.

In order to determine the solution structure of Paracoccus denitrificans cytochrome c552 by NMR, we cloned and isotopically labeled a 10.5-kDa soluble fragment (100 residues) containing the functional domain of the 18.2-kDa membrane-bound protein. Using uniformly 15N-enriched samples of cytochrome c552 in the reduced state, a variety of two-dimensional and three-dimensional heteronuclear double-resonance NMR experiments was employed to achieve complete 1H and 15N assignments. A total of 1893 distance restraints was derived from homonuclear 2D-NOESY and heteronuclear 3D-NOESY spectra; 1486 meaningful restraints were used in the structure calculations. After restrained energy minimization a family of 20 structures was obtained with rmsd values of 0.56 +/- 0. 10 A and 1.09 +/- 0.09 A for the backbone and heavy atoms, respectively. The overall topology is similar to that seen in previously reported models of this class of proteins. The global fold consists of two long helices at the N-terminus and C-terminus and three shorter helices surrounding the heme moiety; the helices are connected by well-defined loops. Comparison with the X-ray structure shows some minor differences in the positions of the Trp57 and Phe65 side-chain rings as well as the heme propionate groups.     

1H and 15N NMR assignment and solution structure of the SH3 domain of spectrin: Comparison of unrefined and refined structure sets with the crystal structure

The assignment of the ¹H and ¹⁵Nnuclear magnetic resonance spectra of the Src-homology region 3 domain ofchicken brain -spectrin has been obtained. A set of solutionstructures has been determined from distance and dihedral angle restraints,which provide a reasonable representation of the protein structure insolution, as evaluated by a principal component analysis of the globalpairwise root-mean-square deviation (rmsd) in a large set of structuresconsisting of the refined and unrefined solution structures and the crystalstructure. The solution structure is well defined, with a lower degree ofconvergence between the structures in the loop regions than in the secondarystructure elements. The average pairwise rmsd between the 15 refinedsolution structures is 0.71 ± 0.13 Å for the backbone atoms and1.43 ± 0.14 Å for all heavy atoms. The solution structure isbasically the same as the crystal structure. The average rmsd between the 15refined solution structures and the crystal structure is 0.76 Å forthe backbone atoms and 1.45 ± 0.09 Å for all heavy atoms. Thereare, however, small differences probably caused by intermolecular contactsin the crystal structure.     

Subunit asymmetry in the three-dimensional structure of a human CuZnSOD mutant found in familial amyotrophic lateral sclerosis.

The X-ray crystal structure of a human copper/zinc superoxide dismutase mutant (G37R CuZnSOD) found in some patients with the inherited form of Lou Gehrig's disease (FALS) has been determined to 1.9 angstroms resolution. The two SOD subunits have distinct environments in the crystal and are different in structure at their copper binding sites. One subunit (subunit[intact]) shows a four-coordinate ligand geometry of the copper ion, whereas the other subunit (subunit[broken]) shows a three-coordinate geometry of the copper ion. Also, subunit(intact) displays higher atomic displacement parameters for backbone atoms ((B) = 30 +/- 10 angstroms2) than subunit(broken) ((B) = 24 +/- 11 angstroms2). This structure is the first CuZnSOD to show large differences between the two subunits. Factors that may contribute to these differences are discussed and a possible link of a looser structure to FALS is suggested.     

A proton-NMR investigation of the fully reduced cytochrome c7 from Desulfuromonas acetoxidans. Comparison between the reduced and the oxidized forms.

The solution structure via 1H NMR of the fully reduced form of cytochrome c7 has been obtained. The protein sample was kept reduced by addition of catalytic amounts of Desulfovibrio gigas iron hydrogenase in H2 atmosphere after it had been checked that the presence of the hydrogenase did not affect the NMR spectrum. A final family of 35 conformers with rmsd values with respect to the mean structure of 8.7 +/- 1.5 nm and 12.4 +/- 1.3 nm for the backbone and heavy atoms, respectively, was obtained. A highly disordered loop involving residues 54-61 is present. If this loop is ignored, the rmsd values are 6.2 +/- 1.1 nm and 10.2 +/- 1.0 nm for the backbone and heavy atoms, respectively, which represent a reasonable resolution. The structure was analyzed and compared with the already available structure of the fully oxidized protein. Within the indetermination of the two solution structures, the result for the two redox forms is quite similar, confirming the special structural features of the three-heme cluster. A useful comparison can be made with the available crystal structures of cytochromes c3, which appear to be highly homologous except for the presence of a further heme. Finally, an analysis of the factors affecting the reduction potentials of the heme irons was performed, revealing the importance of net charges in differentiating the reduction potential when the other parameters are kept constant.     

Solution structure and backbone dynamics of the holo form of the frenolicin acyl carrier protein

During polyketide biosynthesis, acyl carrier proteins (ACPs) perform the central role of transferring polyketide intermediates between active sites of polyketide synthase. The 4'-phosphopantetheine prosthetic group of a holo-ACP is a long and flexible arm that can reach into different active sites and provide a terminal sulfhydryl group for the attachment of acyl groups through a thioester linkage. We have determined the solution structure and characterized backbone dynamics of the holo form of the frenolicin acyl carrier protein (fren holo-ACP) by nuclear magnetic resonance (NMR). Unambiguous assignments were made for 433 hydrogen atoms, 333 carbon atoms, and 84 nitrogen atoms, representing a total of 94.6% of the assignable atoms in this protein. From 879 meaningful NOEs and 45 angle constraints, a family of 24 structures has been calculated. The solution structure is composed of three major alpha-helices packed in a bundle with three additional short helices in intervening loops; one of the short helices slowly exchanges between two conformations. Superposition of the major helical regions on the mean structure yields average atomic rmsd values of 0.49 +/- 0.09 and 0.91 +/- 0.08 A for backbone and non-hydrogen atoms, respectively. Although the three-helix bundle fold is conserved among acyl carrier proteins involved in fatty acid synthases and polyketide synthases, a detailed comparison revealed that ACPs from polyketide biosynthetic pathways are more related to each other in tertiary fold than to their homologues from fatty acid biosynthetic pathways. Comparison of the free form of ACPs (NMR structures of fren ACP and the Bacillus subtilis ACP) with the substrate-bound form of ACP (crystal structure of butyryl-ACP from Escherichia coli) suggests that conformational exchange plays a role in substrate binding.     

Solution structures of conformationally equilibrium forms of holo-acyl carrier protein (PfACP) from Plasmodium falciparum provides insight into the mechanism of activation of ACPs

Acyl Carrier Protein (ACP) from the malaria parasite, Plasmodium falciparum (PfACP) in its holo form is found to exist in two conformational states in solution. Unique 3D solution structures of holo-PfACP have been determined for both equilibrium conformations, using high-resolution NMR methods. Twenty high-resolution solution structures for each of the two forms of holo-PfACP have been determined on the basis of 1226 and 1218 unambiguously assigned NOEs (including NOEs between 4'-phosphopantetheine prosthetic group (4'-PP) and protein), 55 backbone dihedral angles and 26 hydrogen bonds. The atomic rmsd values of the determined structures of two equilibrium forms, about the mean coordinates of the backbone and heavy atoms, are 0.48 +/- 0.09 and 0.92 +/- 0.10 and 0.49 +/- 0.08 and 0.97 +/- 0.11 A, respectively. The interaction of 4'-PP with the polypeptide backbone is reported here for the first time for any of the ACPs. The structures of holo-PfACP consist of three well-defined helices that are tightly packed. The structured regions of the molecule are stabilized by extensive hydrophobic interactions. The difference between the two forms arises from a reorientation of the 4'-PP group. The enthalpy difference between the two forms, although small, implies that a conformational switch is essential for the activation of holo-ACP. Sequence and structures of holo-PfACP have been compared with those of the ACPs from type I and type II fatty acid biosynthesis pathways (FAS), in particular with the ACP from rat and the butyryl-ACP from E. coli. The PfACP structure, thus determined has several novel features hitherto not seen in other ACPs.     

Structure and dynamics of copper-free SOD: The protein before binding copper

The solution structure of the copper-free state of a monomeric form of superoxide dismutase (153 amino acids) was determined through (13)C and (15)N labeling. The protein contained two mutations at the native subunit-subunit interface (F50E and G51E) to obtain a soluble monomeric species and a mutation in the active site channel (E133Q). About 93% of carbon atoms, 95% of nitrogen atoms, and 96% of the protons were assigned. A total of 2467 meaningful NOEs and 170 dihedral angles provided a family of 35 conformers with RMSD values of 0.76 +/- 0.09 A for the backbone and 1.22 +/- 0.13 A for all heavy atoms. The secondary structure elements, connected by loops, produce the typical superoxide dismutase Greek key fold, formed by an eight-stranded beta-barrel. The comparison with the copper-bound monomeric and dimeric structures shows that the metal ligands have a conformation very close to that of the copper-bound forms. This feature indicates that the copper-binding site is preorganized and well ordered also in the absence of the copper ion. The active-site channel shows a sizable increase in width, achieving a suitable conformation to receive the copper ion. The histidines ring NH resonances that bind the copper ion and the region around the active-site channel experience, as found from (15)N relaxation studies, conformational exchange processes. The increased width of the channel and the higher mobility of the histidine rings of the copper site in the copper-free form with respect to the holoprotein is discussed in terms of the process of copper insertion.     

The polypeptide fold of the globular domain of histone H5 in solution. A study using nuclear magnetic resonance, distance geometry and restrained molecular dynamics.

The polypeptide fold of the 79-residue globular domain of chicken histone H5 (GH5) in solution has been determined by the combined use of distance geometry and restrained molecular dynamics calculations. The structure determination is based on 307 approximate interproton distance restraints derived from n.m.r. measurements. The structure is composed of a core made up of residues 3-18, 23-34, 37-60 and 71-79, and two loops comprising residues 19-22 and 61-70. The structure of the core is well defined with an average backbone atomic r.m.s. difference of 2.3 +/- 0.3 A between the final eight converged restrained dynamics structures and the mean structure obtained by averaging their coordinates best fitted to the core residues. The two loops are also well defined locally but their orientation with respect to the core could not be determined as no long range ([i-j[ greater than 5) proton-proton contacts could be observed between the loop and core residues in the two-dimensional nuclear Overhauser enhancement spectra. The structure of the core is dominated by three helices and has a similar fold to the C-terminal DNA binding domain of the cAMP receptor protein.     

Biosynthesis of D-alanyl-lipoteichoic acid: the tertiary structure of apo-D-alanyl carrier protein.

The D-alanylation of lipoteichoic acid (LTA) allows the Gram-positive organism to modulate its surface charge, regulate ligand binding, and control the electromechanical properties of the cell wall. The incorporation of D-alanine into LTA requires the D-alanine:D-alanyl carrier protein ligase (AMP-forming) (Dcl) and the carrier protein (Dcp). The high-resolution solution structure of the 81-residue (8.9 kDa) Dcp has been determined by multidimensional heteronuclear NMR. An ensemble of 30 structures was calculated using the torsion angle dynamics approach of DYANA. These calculations utilized 3288 NOEs containing 1582 unique nontrivial NOE distance constraints. Superposition of residues 4-81 on the mean structure yields average atomic rmsd values of 0.43 +/- 0.08 and 0.86 +/- 0.09 A for backbone and non-hydrogen atoms, respectively. The solution structure is composed of three alpha-helices in a bundle with additional short 3(10)- and alpha-helices in intervening loops. Comparisons of the three-dimensional structure with the acyl carrier proteins involved in fatty acid, polyketide, and nonribosomal peptide syntheses support the conclusion that Dcp is a homologue in this family. While there is conservation of the three-helix bundle fold, Dcp has a higher enthalpy of unfolding and no apparent divalent metal binding site(s), features that distinguish it from the fatty acid synthase acyl carrier protein of Escherichia coli. This three-dimensional structure also provides insights into the D-alanine ligation site recognized by Dcl, as well as the site which may bind the poly(glycerophosphate) acceptor moiety of membrane-associated LTA.     

High resolution solution structure of the protein part of Cu7 metallothionein.

The three-dimensional solution structure of the protein part of Cu7 metallothionein (Cu7MT) of Saccharomyces cerevisiae has been attempted by 1H two-dimensional NMR spectroscopy at 800 MHz. The protein part constitutes 53 amino acids. A total of 1192 NOEs, of which 1048 are meaningful, were used to determine the solution structure of the first 40 residues, the last 13 residues being disordered. A family of 30 structures was generated. Root-mean-square deviation (rmsd) values from the average structure of 0.32 +/- 0.13 A and 0.61 +/- 0.15 A for backbone and all heavy atoms, respectively, were obtained for the residues 2-40. The ten copper-coordinating cysteine sulfurs and the empty spaces around them are well defined. The structure of the protein part is similar but not identical to the available ones of the same holoprotein and of the Ag7 metallothionein, and is qualitatively superior. If the same metal-sulfur connectivities reported in the literature from 1H-109Ag heteronuclear multiple quantum coherence spectroscopy are assumed to hold for the present copper derivative, a peptide structure is obtained which is again similar, but still not identical, within indetermination, to that available. The structure of the copper polymetallic center may well be different from that proposed for the silver derivative, and indeed a number of different arrangements of the seven copper ions are consistent with the present highly refined structure of the protein part.     

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.     

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.     

NMR solution structure of plastocyanin from the photosynthetic prokaryote, Prochlorothrix hollandica

The solution structure of a divergent plastocyanin (PC) from the photosynthetic prokaryote Prochlorothrix hollandica was determined by homonuclear 1H NMR spectroscopy. Nineteen structures were calculated from 1222 distance restraints, yielding a family of structures having an average rmsd of 0.42 +/- 0.08 A for backbone atoms and 0.71 +/- 0.07 A for heavy atoms to the mean structure. No distance constraint was violated by more than 0.26 A in the structure family. Despite the low number of conserved residues shared with other PC homologues, the overall folding pattern of P. hollandica PC is similar to other PCs, in that the protein forms a two-sheet beta-barrel tertiary structure. The greatest variability among the backbone structures is seen in the loop region from residues 47-60. The differences seen in the P. hollandica PC homologue likely arise due to a small deletion of 2-4 residues compared to the PC consensus; this yields a less extended loop containing a short alpha-helix from residues Ala52-Leu55. Additionally, the protein has an altered hydrophobic patch thought to be important in binding reaction partners. Whereas the backbone structure is very similar within the loops of the hydrophobic region, the presence of two unique residues (Tyr12 and Pro14) yields a structurally different hydrophobic surface likely important in binding P. hollandica Photosystem I.     

Solution structure of human beta-parvalbumin and structural comparison with its paralog alpha-parvalbumin and with their rat orthologs

The aim of this research was to determine the structure of human beta-parvalbumin (109 amino acids) and to compare it with its paralog and ortholog proteins. The structure was determined in solution using multinuclear and multidimensional NMR methods and refined using substitution of the EF-hand Ca(2+) ion with a paramagnetic lanthanide. The resulting family of structures had a backbone rmsd of 0.50 A. Comparison with rat oncomodulin (X-ray, 1.3 A resolution) as well as with human (NMR, backbone rmsd of 0.49 A) and rat (X-ray, 2.0 A resolution) parvalbumins reveals small but reliable local differences, often but not always related to amino acid variability. The analysis of these structures has led us to propose an explanation for the different affinity for Ca(2+) between alpha- and beta-parvalbumins and between parvalbumins and calmodulins.     

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.     

High-resolution solution structure of gurmarin, a sweet-taste-suppressing plant polypeptide.

Gurmarin is a 35-residue polypeptide from the Asclepiad vine Gymnema sylvestre. It has been utilised as a pharmacological tool in the study of sweet-taste transduction because of its ability to selectively inhibit the neural response to sweet tastants in rats. We have chemically synthesised and folded gurmarin and determined its three-dimensional solution structure to high resolution using two-dimensional NMR spectroscopy. Structure calculations utilised 612 interproton-distance, 19 dihedral-angle, and 18 hydrogen-bond restraints. The structure is well defined for residues 3-34, with backbone and heavy atom rms differences of 0.27 +/- 0.09 A and 0.73 +/- 0.09 A, respectively. Gurmarin adopts a compact structure containing an antiparallel beta-hairpin (residues 22-34), several well-defined beta-turns, and a cystine-knot motif commonly observed in toxic and inhibitory polypeptides. Despite striking structural homology with delta-atracotoxin, a spider neurotoxin known to slow the inactivation of voltage-gated Na+ channels, we show that gurmarin has no effect on a variety of voltage-sensitive channels.     

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.     

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.     

New efficient statistical sequence-dependent structure prediction of short to medium-sized protein loops based on an exhaustive loop classification.

A bank of 13,563 loops from three to eight amino acid residues long, representing motifs between two consecutive regular secondary structures, has been derived from protein structures presenting less than 95 % sequence identity. Statistical analyses of occurrences of conformations and residues revealed length-dependent over-representations of particular amino acids (glycine, proline, asparagine, serine, and aspartate) and conformations (alphaL, epsilon, betaPregions of the Ramachandran plot). A position-dependent distribution of these occurrences was observed for N and C-terminal residues, which are correlated to the nature of the flanking regions. Loops of the same length were clustered into statistically meaningful families on the basis of their backbone structures when placed in a common reference frame, independent of the flanks. These clusters present significantly different distributions of sequence, conformations, and endpoint residue Calphadistances. On the basis of the sequence-structure correlation of this clustering, an automatic loop modeling algorithm was developed. Based on the knowledge of its sequence and of its flank backbone structures each query loop is assigned to a family and target loop supports are selected in this family. The support backbones of these target loops are then adjusted on flanking structures by partial exploration of the conformational space. Loop closure is performed by energy minimization for each support and the final model is chosen among connected supports based upon energy criteria. The quality of the prediction is evaluated by the root-mean-square deviation (rmsd) between the final model and the native loops when the whole bank is re-attributed on itself with a Jackknife test. This average rmsd ranges from 1.1 A for three-residue loops to 3.8 A for eight-residue loops. A few poorly predicted loops are inescapable, considering the high level of diversity in loops and the lack of environment data. To overcome such modeling problems, a statistical reliability score was assigned for each prediction. This score is correlated to the quality of the prediction, in terms of rmsd, and thus improves the selection accuracy of the model. The algorithm efficiency was compared to CASP3 target loop predictions. Moreover, when tested on a test loop bank, this algorithm was shown to be robust when the loops are not precisely delimited, therefore proving to be a useful tool in practice for protein modeling.     

The "second stalk" of Escherichia coli ATP synthase: structure of the isolated dimerization domain

The b subunit of E. coli F(0)F(1)-ATPase links the peripheral F(1) subunits to the membrane-integral F(0) portion and functions as a "stator", preventing rotation of F(1). The b subunit is present as a dimer in ATP synthase, and residues 62-122 are required to mediate dimerization. To understand how the b subunit dimer is formed, we have studied the structure of the isolated dimerization domain, b(62-122). Analytical ultracentrifugation and solution small-angle X-ray scattering (SAXS) indicate that the b(62-122) dimer is extremely elongated, with a frictional ratio of 1.60, a maximal dimension of 95 A, and a radius of gyration of 27 A, values that are consistent with an alpha-helical coiled-coil structure. The crystal structure of b(62-122) has been solved and refined to 1.55 A. The protein crystallized as an isolated, monomeric alpha helix with a length of 90 A. Combining the crystal structure of monomeric b(62-122) with SAXS data from the dimer in solution, we have constructed a model for the b(62-122) dimer in which the two helices form a coiled coil with a right-handed superhelical twist. Analysis of b sequences from E. coli and other prokaryotes indicates conservation of an undecad repeat, which is characteristic of a right-handed coiled coil and consistent with our structural model. Mutation of residue Arg-83, which interrupts the undecad pattern, to alanine markedly stabilized the dimer, as expected for the proposed two-stranded, right-handed coiled-coil structure.