Delineation of protein structure classes from multivariate analysis of protein Raman optical activity data

Vibrational Raman optical activity (ROA), measured as a small difference in the intensity of Raman scattering from chiral molecules in right and left-circularly polarized incident light, or as the intensity of a small circularly polarized component in the scattered light, is a powerful probe of the aqueous solution structure of proteins. On account of the large number of structure-sensitive bands in protein ROA spectra, multivariate analysis techniques such as non-linear mapping (NLM) are especially favourable for determining structural relationships between different proteins. Here NLM is used to map a dataset of 80 polypeptide, protein and virus ROA spectra, considered as points in a multidimensional space with axes representing the digitized wavenumbers, into readily visualizable two and three-dimensional spaces in which points close to or distant from each other, respectively, represent similar or dissimilar structures. Discrete clusters are observed which correspond to the seven structure classes all alpha, mainly alpha, alphabeta, mainly beta, all beta, mainly disordered/irregular and all disordered/irregular. The average standardised ROA spectra of the proteins falling within each structure class have distinct features characteristic of each class. A distinct cluster containing the wheat protein A-gliadin and the plant viruses potato virus X, narcissus mosaic virus, papaya mosaic virus and tobacco rattle virus, all of which appear in the mainly alpha cluster in the two-dimensional representation, becomes clearly separated in the direction of increasing disorder in the three-dimensional representation. This suggests that the corresponding five proteins, none of which to date has yielded high-resolution X-ray structures, consist mainly of alpha-helix and disordered structure with little or no beta-sheet. This combination of structural elements may have functional significance, such as facilitating disorder-to-order transitions (and vice versa) and suppressing aggregation, in these proteins and also in sequences within other proteins. The use of ROA to identify proteins containing significant amounts of disordered structure will, inter alia, be valuable in structural genomics/proteomics since disordered regions often inhibit crystallization.     

Changes in protein conformation and dynamics upon complex formation of brain-derived neurotrophic factor and its receptor: investigation by isotope-edited Fourier transform IR spectroscopy

The interactions of brain-derived neurotrophic factor (BDNF) with the extracellular domain of its receptor (trkB) are investigated by employing isotope-edited Fourier transform IR (FTIR) spectroscopy. The protein secondary structures of individual BDNF and trkB in solutions are compared with those in their complex. The temperature dependence of the secondary structures of BDNF, trkB, and their complex is also investigated. Consistent with the crystal structure, we observe by FTIR spectroscopy that BDNF in solution contains predominantly beta strands (approximately 53%) and relatively low contents of other secondary structures including beta turns (approximately 16%), disordered structures (approximately 12%), and loops (approximately 18%) and is deficient in alpha helix. We also observe that trkB in solution contains mostly beta strands (52%) and little alpha helix. Conformational changes in both BDNF and trkB are observed upon complex formation. Specifically, upon binding of BDNF, the conformational changes in trkB appear to involve mostly beta turns and disordered structures while the majority of the beta-strand conformation remains unchanged. The IR data indicate that some of the disordered structures in the loop regions are likely converted to beta strands upon complex formation. The FTIR spectral data of BDNF, trkB, and their complex indicate that more amide NH groups of trkB undergo H-D exchange within the complex than those of the ligand-free receptor and that the thermal stability of trkB is decreased slightly upon binding of BDNF. The FT-Raman spectra of BDNF, trkB, and their complex show that the six intramolecular disulfide bonds of trkB undergo significant conformational changes upon binding of BDNF as a result of changes in the tertiary structure of trkB. Taken together, the FTIR and Raman data are consistent with the loosening of the tertiary structure of trkB upon binding of BDNF, which leads to more solvent exposure of the amide NH group and decreased thermal stability of trkB. This finding reveals an intriguing structural property of the neurotrophin ligand-receptor complex that is in contrast to other ligand-receptor complexes such as a cytokine-receptor complex that usually shows protection of the amide NH group and increased thermal stability upon complex formation.     

Addition of missing loops and domains to protein models by x-ray solution scattering

Inherent flexibility and conformational heterogeneity in proteins can often result in the absence of loops and even entire domains in structures determined by x-ray crystallographic or NMR methods. X-ray solution scattering offers the possibility of obtaining complementary information regarding the structures of these disordered protein regions. Methods are presented for adding missing loops or domains by fixing a known structure and building the unknown regions to fit the experimental scattering data obtained from the entire particle. Simulated annealing was used to minimize a scoring function containing the discrepancy between the experimental and calculated patterns and the relevant penalty terms. In low-resolution models where interface location between known and unknown parts is not available, a gas of dummy residues represents the missing domain. In high-resolution models where the interface is known, loops or domains are represented as interconnected chains (or ensembles of residues with spring forces between the C(alpha) atoms), attached to known position(s) in the available structure. Native-like folds of missing fragments can be obtained by imposing residue-specific constraints. After validation in simulated examples, the methods have been applied to add missing loops or domains to several proteins where partial structures were available.     

Solution structure of the phosphocarrier protein HPr from Bacillus subtilis by two-dimensional NMR spectroscopy.

The solution structure of the phosphocarrier protein, HPr, from Bacillus subtilis has been determined by analysis of two-dimensional (2D) NMR spectra acquired for the unphosphorylated form of the protein. Inverse-detected 2D (1H-15N) heteronuclear multiple quantum correlation nuclear Overhauser effect (HMQC NOESY) and homonuclear Hartmann-Hahn (HOHAHA) spectra utilizing 15N assignments (reported here) as well as previously published 1H assignments were used to identify cross-peaks that are not resolved in 2D homonuclear 1H spectra. Distance constraints derived from NOESY cross-peaks, hydrogen-bonding patterns derived from 1H-2H exchange experiments, and dihedral angle constraints derived from analysis of coupling constants were used for structure calculations using the variable target function algorithm, DIANA. The calculated models were refined by dynamical simulated annealing using the program X-PLOR. The resulting family of structures has a mean backbone rmsd of 0.63 A (N, C alpha, C', O atoms), excluding the segments containing residues 45-59 and 84-88. The structure is comprised of a four-stranded antiparallel beta-sheet with two antiparallel alpha-helices on one side of the sheet. The active-site His 15 residue serves as the N-cap of alpha-helix A, with its N delta 1 atom pointed toward the solvent to accept the phosphoryl group during the phosphotransfer reaction with enzyme I. The existence of a hydrogen bond between the side-chain oxygen atom of Tyr 37 and the amide proton of Ala 56 is suggested, which may account for the observed stabilization of the region that includes the beta-turn comprised of residues 37-40. If the beta alpha beta beta alpha beta (alpha) folding topology of HPr is considered with the peptide chain polarity reversed, the protein fold is identical to that described for another group of beta alpha beta beta alpha beta proteins that include acylphosphatase and the RNA-binding domains of the U1 snRNP A and hnRNP C proteins.     

The HPV16 E7 viral oncoprotein self-assembles into defined spherical oligomers

Despite the fact that E7 is a major transforming oncoprotein in papillomavirus, its structure and precise molecular mechanism of action remain puzzling to date. E7 proteins share sequence homology and proteasome targeting properties of tumor suppressors with adenovirus E1A and SV40 T antigen, two other paradigmatic oncoproteins from DNA tumor viruses. High-risk HPV16 E7, a nonglobular dimer with some properties of intrinsically disordered proteins, is capable of undergoing pH-dependent conformational transitions that expose hydrophobic surfaces to the solvent. We found that treatment with a chelating agent produced a protein that can readily assemble into homogeneous spherical particles with an average molecular mass of 790 kDa and a diameter of 50 nm, as determined from dynamic light scattering and electron microscopy. The protein undergoes a substantial conformational transition from coil to beta-sheet structure, with concomitant consolidation of tertiary structure as judged by circular dichroism and fluorescence. The assembly process is very slow, in agreement with a substantial energy barrier caused by structural rearrangements. The resulting particles are highly stable, cooperatively folded, and capable of binding both Congo Red and thioflavin T, reporters of repetitive beta-sheet structures similar to those found in amyloids, although no fibrillar or insoluble material was observed under our experimental conditions.     

Differential dehydration effects on globular proteins and intrinsically disordered proteins during film formation

Globular proteins composed of different secondary structures and fold types were examined by synchrotron radiation circular dichroism spectroscopy to determine the effects of dehydration on their secondary structures. They exhibited only minor changes upon removal of bulk water during film formation, contrary to previously reported studies of proteins dehydrated by lyophilization (where substantial loss of helical structure and gain in sheet structure was detected). This near lack of conformational change observed for globular proteins contrasts with intrinsically disordered proteins (IDPs) dried in the same manner: the IDPs, which have almost completely unordered structures in solution, exhibited increased amounts of regular (mostly helical) secondary structures when dehydrated, suggesting formation of new intra‐protein hydrogen bonds replacing solvent‐protein hydrogen bonds, in a process which may mimic interactions that occur when IDPs bind to partner molecules. This study has thus shown that the secondary structures of globular and intrinsically disordered proteins behave very differently upon dehydration, and that films are a potentially useful format for examining dehydrated soluble proteins and assessing IDPs structures.     

Model systems for beta-hairpins and beta-sheets

beta-Sheets and alpha-helices are the two principal secondary structures in proteins. However, our understanding of beta-sheet structure lags behind that of alpha-helices, largely because, until recently, there was no model system to study the beta-sheet secondary structure in isolation. With the development of well-folded beta-hairpins, this is changing rapidly. Recent advances include: increased understanding of the relative contributions of turn, strand and sidechain interactions to beta-hairpin and beta-sheet stability, with the role of aromatic residues as a common subtheme; experimental and theoretical kinetic and thermodynamic studies of beta-hairpin and beta-sheet folding; de novo protein design, including all-beta structures, mixed alpha/beta motifs and switchable systems; and the creation of functional beta-hairpins.     

Standard conformations of a polypeptide chain in irregular protein regions

A detailed stereochemical analysis of known protein structures has been made which shows that: (1) irregular regions of proteins consist of a limited number of standard structures formed by three, four of more residues; (2) an amino acid residue of a protein can adopt one of the six sterically allowed conformations designated here as alpha, alpha L, beta, gamma, delta, and epsilon. It is shown that there are two allowed conformations of a polypeptide chain at the N-end of an alpha-helix, beta alpha n- and beta gamma alpha n-conformations, where n is a number of residues in the alpha-helix. At the C-end of the alpha-helix there are two conformations as well, alpha n gamma beta- and alpha n gamma alpha L beta-ones. Two beta-strands in a beta-hairpin can be joined, for example, by standard structures with beta beta alpha L beta-, beta alpha gamma alpha L beta-, beta alpha alpha gamma alpha L beta-conformations which are referred to as turns. In the regions where a polypeptide chain passes from one layer to another there are standard structures with beta gamma beta-, beta alpha beta beta-, beta alpha gamma beta-conformations etc., referred to as cross-overs. A structure of any protein irregular region can be represented as a combination of these and other standard turns and cross-overs considered in the paper. The major part of the turns and cross-overs has residues in alpha L- or epsilon-conformations which must be glycine or other residues with small or flexible side chains. Massive hydrophobic residues must not occupy the first beta-positions of the most standard structures. The results obtained can be successfully applied for prediction of the location of the turns and cross-overs in proteins from their amino acid sequences and for interpretation of electron density maps.     

NMR characterization of three-disulfide variants of lysozyme, C64A/C80A, C76A/C94A, and C30A/C115A--a marginally stable state in folded proteins

Our earlier NMR study showed that a two-disulfide variant of hen lysozyme containing intra-alpha-domain disulfide bridges, C6-C127 and C30-C115, is partially folded, with the alpha domain tightly folded to the nativelike conformation and the beta domain flexible or unfolded. With a view that the formation of a third disulfide bridge is a key for the accomplishment of the overall chain fold, three-dimensional structures of three-disulfide variants of hen lysozyme lacking one disulfide bridge (C64A/C80A, C76A/C94A, and C30A/C115A) were studied in detail using NMR spectroscopy. Amide hydrogen exchange rates were measured to estimate the degree of conformational fluctuation in a residue-specific manner. The structure of C76A/C94A was found to be quite similar to that of the wild type, except for the peptide segment of residues 74-78. The structure of C64A/C80A was considerably disordered in the entire region of the loop (residues 62-79). Further, it was found that a network of hydrogen bonds within the beta sheet and the 3(10) helix in the beta domain were disrupted and fluctuating. In C30A/C115A, the D helix was unstructured and the interface of the B helix with the D helix was significantly perturbed. However, the structural disorder generated in the hydrophobic core of the alpha domain was prevented by the C helix from propagating toward the beta domain. A marginally stable state in folded proteins is discussed based on the structures remaining in each three-disulfide variant.     

Structure-function mapping of interleukin 1 precursors. Cleavage leads to a conformational change in the mature protein

The two interleukin 1 (IL-1) genes (IL-1 alpha and beta) encode 31-kDa precursor molecules, which are cleaved upon secretion to generate the mature, active, carboxyl-terminal 17-kDa proteins. The IL-1 beta precursor is inactive, whereas the IL-1 alpha precursor is as active as the mature IL-1 alpha. In this report, we demonstrate that when either of the recombinant precursors is processed to the mature form, the mature region undergoes a conformational change from a proteinase K-sensitive structure to one that is proteinase K-insensitive. In addition, cysteine residues that are exposed to solvent in the IL-1 beta precursor become buried in the mature protein. Limited structure-activity mapping of the IL-1 beta precursor indicates that the amino-terminal 76 residues are responsible for the conformational change, whereas the most dramatic change in biological activity occurs after further removal of residues 77-94. These findings suggest that the altered structure of the mature region in precursor IL-1s has been conserved for some function. Denaturation/renaturation experiments implicate the precursor domain in protein folding, and by analogy with signal-directed secretory proteins, the unique conformation of the precursors may play a role in IL-1 secretion.     

Temperature-induced reversible conformational change in the first 100 residues of alpha-synuclein

Natively disordered proteins are a growing class of anomalies to the structure-function paradigm. The natively disordered protein alpha-synuclein is the primary component of Lewy bodies, the cellular hallmark of Parkinson's disease. We noticed a dramatic difference in dilute solution 1H-15N Heteronuclear Single Quantum Coherence (HSQC) spectra of wild-type alpha-synuclein and two disease-related mutants (A30P and A53T), with spectra collected at 35 degrees C showing fewer cross-peaks than spectra acquired at 10 degrees C. Here, we show the change to be the result of a reversible conformational exchange linked to an increase in hydrodynamic radius and secondary structure as the temperature is raised. Combined with analytical ultracentrifugation data showing alpha-synuclein to be monomeric at both temperatures, we conclude that the poor quality of the 1H-15N HSQC spectra obtained at 35 degrees C is due to conformational fluctuations that occur on the proton chemical shift time scale. Using a truncated variant of alpha-synuclein, we show the conformational exchange occurs in the first 100 amino acids of the protein. Our data illustrate a key difference between globular and natively disordered proteins. The properties of globular proteins change little with solution conditions until they denature cooperatively, but the properties of natively disordered proteins can vary dramatically with solution conditions.     

Energy Landscapes of Protein Aggregation and Conformation Switching in Intrinsically Disordered Proteins

The protein folding problem was apparently solved recently by the advent of a deep learning method for protein structure prediction called AlphaFold. However, this program is not able to make predictions about the protein folding pathways. Moreover, it only treats about half of the human proteome, as the remaining proteins are intrinsically disordered or contain disordered regions. By definition these proteins differ from natively folded proteins and do not adopt a properly folded structure in solution. However these intrinsically disordered proteins (IDPs) also systematically differ in amino acid composition and uniquely often become folded upon binding to an interaction partner. These factors preclude solving IDP structures by current machine-learning methods like AlphaFold, which also cannot solve the protein aggregation problem, since this meta-folding process can give rise to different aggregate sizes and structures. An alternative computational method is provided by molecular dynamics simulations that already successfully explored the energy landscapes of IDP conformational switching and protein aggregation in multiple cases. These energy landscapes are very different from those of ‘simple’ protein folding, where one energy funnel leads to a unique protein structure. Instead, the energy landscapes of IDP conformational switching and protein aggregation feature a number of minima for different competing low-energy structures. In this review, I discuss the characteristics of these multifunneled energy landscapes in detail, illustrated by molecular dynamics simulations that elucidated the underlying conformational transitions and aggregation processes.     

A comparative study of the relationship between protein structure and beta-aggregation in globular and intrinsically disordered proteins

A growing number of proteins are being identified that are biologically active though intrinsically disordered, in sharp contrast with the classic notion that proteins require a well-defined globular structure in order to be functional. At the same time recent work showed that aggregation and amyloidosis are initiated in amino acid sequences that have specific physico-chemical properties in terms of secondary structure propensities, hydrophobicity and charge. In intrinsically disordered proteins (IDPs) such sequences would be almost exclusively solvent-exposed and therefore cause serious solubility problems. Further, some IDPs such as the human prion protein, synuclein and Tau protein are related to major protein conformational diseases. However, this scenario contrasts with the large number of unstructured proteins identified, especially in higher eukaryotes, and the fact that the solubility of these proteins is often particularly good. We have used the algorithm TANGO to compare the beta aggregation tendency of a set of globular proteins derived from SCOP and a set of 296 experimentally verified, non-redundant IDPs but also with a set of IDPs predicted by the algorithms DisEMBL and GlobPlot. Our analysis shows that the beta-aggregation propensity of all-alpha, all-beta and mixed alpha/beta globular proteins as well as membrane-associated proteins is fairly similar. This illustrates firstly that globular structures possess an appreciable amount of structural frustration and secondly that beta-aggregation is not determined by hydrophobicity and beta-sheet propensity alone. We also show that globular proteins contain almost three times as much aggregation nucleating regions as IDPs and that the formation of highly structured globular proteins comes at the cost of a higher beta-aggregation propensity because both structure and aggregation obey very similar physico-chemical constraints. Finally, we discuss the fact that although IDPs have a much lower aggregation propensity than globular proteins, this does not necessarily mean that they have a lower potential for amyloidosis.     

Ions binding to S100 proteins. II. Conformational studies and calcium-induced conformational changes in S100 alpha alpha protein: the effect of acidic pH and calcium incubation on subunit exchange in S100a (alpha beta) protein

A rapid separation method for bovine brain S100 alpha alpha, S100a, and S100b protein using fast protein liquid chromatography on a Mono Q column and its application in preparation of a large amount of S100 alpha alpha protein are described. The conformation of S100 alpha alpha in the metal-free forms as well as in the presence of calcium were studied by UV absorption, circular dichroism, intrinsic fluorescence, sulfhydryl reactivity, and interaction with a hydrophobic fluorescent probe. The alpha-subunit appears to have nearly identical conformation in S100 alpha alpha and S100a protein dimers. We also confirmed that only the alpha-subunit exposes hydrophobic domains to solvent in the presence of calcium and that cysteine residues exposed upon Ca2+ binding to S100 proteins correspond to Cys 85 alpha and Cys 84 beta. Incubation of S100a with calcium and KCl proved that calcium binding to the putative calcium-binding sites (site I alpha, I beta) triggers a time- and temperature-dependent conformational change in the protein structure which decreases the antagonistic effect of KCl on calcium binding to sites II alpha and II beta and provokes subunit exchanges between protein dimers and the emergence of S100 alpha alpha and S100b (beta beta) proteins. Dynamic fluorescence measurements showed that incubating calcium at high S100a protein concentrations (greater than 10(-5) M) induces an apparent slow dimer-monomer equilibrium which might result in total subunit dissociation at lower protein concentrations. The effect of acidic pH on subunit dissociation in S100a protein (Morero, R. D., and Weber, G. (1982) Biochim. Biophys. Acta 703, 231-240) arises from conformational changes in the protein structure that are similar to those induced by Ca2+ incubation.     

Probing the activation-promoted structural rearrangements in preassembled receptor-G protein complexes

Activation of heterotrimeric G proteins by their cognate seven transmembrane domain receptors is believed to involve conformational changes propagated from the receptor to the G proteins. However, the nature of these changes remains unknown. We monitored the conformational rearrangements at the interfaces between receptors and G proteins and between G protein subunits by measuring bioluminescence resonance energy transfer between probes inserted at multiple sites in receptor-G protein complexes. Using the data obtained for the alpha(2A)AR-G alpha(i1) beta1gamma2 complex and the available crystal structures of G alpha(i1) beta1gamma2, we propose a model wherein agonist binding induces conformational reorganization of a preexisting receptor-G protein complex, leading the G alpha-G betagamma interface to open but not dissociate. This conformational change may represent the movement required to allow nucleotide exit from the G alpha subunit, thus reflecting the initial activation event.     

Intrinsically Disordered Proteins and Their Environment: Effects of Strong Denaturants, Temperature, pH, Counter Ions, Membranes, Binding Partners, Osmolytes, and Macromolecular Crowding

Intrinsically disordered proteins (IDPs) differ from “normal” ordered proteins at several levels, structural, functional and conformational. Amino acid biases characteristic for IDPs determine their structural variability and lack of rigid well-folded structure. This structural plasticity is necessary for the unique functional repertoire of IDPs, which is complementary to the catalytic activities of ordered proteins. Amino acid biases also drive atypical responses of IDPs to changes in their environment. The conformational behavior of IDPs is characterized by the low cooperativity (or the complete lack thereof) of the denaturant-induced unfolding, lack of the measurable excess heat absorption peak(s) characteristic for the melting of ordered proteins, “turned out” response to heat and changes in pH, the ability to gain structure in the presence of various counter ions, osmolytes, membranes and binding partners, and by the unique response to macromolecular crowding. This review describes some of the most characteristic features of the IDP conformational behavior and the unique response of IDPs to changes in their environment.     

3D structure and significance of the GPhiXXG helix packing motif in tetramers of the E1beta subunit of pyruvate dehydrogenase from the archeon Pyrobaculum aerophilum.

As part of a structural genomics project, we have determined the 2.0 A structure of the E1beta subunit of pyruvate dehydrogenase from Pyrobaculum aerophilum (PA), a thermophilic archaeon. The overall fold of E1beta from PA is closely similar to the previously determined E1beta structures from humans (HU) and P. putida (PP). However, unlike the HU and PP structures, the PA structure was determined in the absence of its partner subunit, E1alpha. Significant structural rearrangements occur in E1beta when its E1alpha partner is absent, including rearrangement of several secondary structure elements such as helix C. Helix C is buried by E1alpha in the HU and PP structures, but makes crystal contacts in the PA structure that lead to an apparent beta(4) tetramer. Static light scattering and sedimentation velocity data are consistent with the formation of PA E1beta tetramers in solution. The interaction of helix C with its symmetry-related counterpart stabilizes the tetrameric interface, where two glycine residues on the same face of one helix create a packing surface for the other helix. This GPhiXXG helix-helix interaction motif has previously been found in interacting transmembrane helices, and is found here at the E1alpha-E1beta interface for both the HU and PP alpha(2)beta(2) tetramers. As a case study in structural genomics, this work illustrates that comparative analysis of protein structures can identify the structural significance of a sequence motif.     

The conformation of glucagon: predictions and consequences.

It is proposed that glucagon, a polypeptide hormone, is delicately balanced between two major conformational states. Utilizing a new predictive model [Chou, P.Y., and Fasman, G.D. (1974), Biochemistry 13, 222] which considers all the conformational states in proteins (helix, beta sheet, random coil, and beta turns), the secondary structural regions of glucagon are computed herein. The conformational sensitivity of glucagon may be due to residues 19-27 which have both alpha-helical potential (mean value of Palpha = 1.19) as well as beta-sheet potential (mean value of Pbeta = 1.25). Two conformational states are predicted for glucagon. In predicted form (a), residues 5-10 form a beta-sheet region while residues 19-27 form an alpha-helical region (31% alpha, 21% beta) agreeing well with the circular dichroism (CD) spectra of glucagon. The similarity in the CD spectra of glucagon and insulin further suggests the presence of beta structure in glucagon, since X-ray analysis of insulin showed 24% beta sheet. In predicted form (b), both regions, residues 5-10 and residues 19-27, are beta sheets sheets (0% alpha, 52% beta) in agreement with the infrared spectral evidence that glucagon gels and fibrils have a predominant beta-sheet conformation. Since three reverse beta turns are predicted at residues 2-5, 10-13, and 15-18, glucagon may possess tertiary structure in agreement with viscosity and tritium-hydrogen exchange experiments. A proposal is offered concerning an induced alpha yields beta transition at residues 22-27 in glucagon during receptor site binding. Amino acid substitutions are proposed which should disrupt the beta sheets of glucagon with concomitant loss of biological activity. The experimental findings that glucagon aggregates to form dimers, trimers, and hexamers can be explained in terms of beta-sheet interactions as outlined in the present predictive model. Thus the conflicting conclusions of previous workers, concerning the conformation of glucagon in different environments, can be rationalized by the suggested conformational transition occurring within the molecule.     

Optimal relationship between average conformational entropy and average energy of interactions between residues for fast protein folding

Based on the known experimental data and using the theoretical modeling of protein folding, we demonstrate that there exists an optimal relationship between the average conformational entropy and the average energy of contacts per residue, that is an entropy capacity, for fast protein folding. Statistical analysis of conformational entropy and the number of contacts per residue for 5829 protein structures from four general structural classes (all-alpha, all-beta, +/-/beta, alpha+beta) demonstrates that each class of proteins has its own class-specific average number of contacts and average conformational entropy per residue. These class-specific features determine the folding rates: a proteins are the fastest folding proteins, then follow beta and alpha+beta proteins, and finally alpha/beta proteins are the slowest ones.     

Helicity of alpha(404-451) and beta(394-445) tubulin C-terminal recombinant peptides

We have investigated the solution conformation of the functionally relevant C-terminal extremes of alpha- and beta-tubulin, employing the model recombinant peptides RL52alpha3 and RL33beta6, which correspond to the amino acid sequences 404-451(end) and 394-445(end) of the main vertebrate isotypes of alpha- and beta-tubulin, respectively, and synthetic peptides with the alpha-tubulin(430-443) and beta-tubulin(412-431) internal sequences. Alpha(404-451) and beta(394-445) are monomeric in neutral aqueous solution (as indicated by sedimentation equilibrium), and have circular dichroism (CD) spectra characteristic of nearly disordered conformation, consistent with low scores in peptide helicity prediction. Limited proteolysis of beta(394-445) with subtilisin, instead of giving extensive degradation, resulted in main cleavages at positions Thr409-Glu410 and Tyr422-Gln423-Gln424, defining the proteolysis resistant segment 410-422, which corresponds to the central part of the predicted beta-tubulin C-terminal helix. Both recombinant peptides inhibited microtubule assembly, probably due to sequestration of the microtubule stabilizing associated proteins. Trifluoroethanol (TFE)-induced markedly helical CD spectra in alpha(404-451) and beta(394-445). A substantial part of the helicity of beta(394-445) was found to be in the CD spectrum of the shorter peptide beta(412-431) with TFE. Two-dimensional 1H-NMR parameters (nonsequential nuclear Overhauser effects (NOE) and conformational C alphaH shifts) in 30% TFE permitted to conclude that about 25% of alpha(404-451) and 40% of beta(394-451) form well-defined helices encompassing residues 418-432 and 408-431, respectively, flanked by disordered N- and C-segments. The side chains of beta(394-451) residues Leu418, Val419, Ser420, Tyr422, Tyr425, and Gln426 are well defined in structure calculations from the NOE distance constraints. The apolar faces of the helix in both alpha and beta chains share a characteristic sequence of conserved residues Ala,Met(+4),Leu(+7),Tyr(+11). The helical segment of alpha(404-451) is the same as that described in the electron crystallographic model structure of alphabeta-tubulin, while in beta(394-451) it extends for nine residues more, supporting the possibility of a functional coil --> helix transition at the C-terminus of beta-tubulin. These peptides may be employed to construct model complexes with microtubule associated protein binding sites.