Folding propensities of synthetic peptide fragments covering the entire sequence of phage 434 Cro protein.

The phage 434 Cro protein, the N-terminal domain of its repressor (R1-69) and that of phage lambda (lambda6-85) constitute a group of small, monomeric, single-domain folding units consisting of five helices with striking structural similarity. The intrinsic helix stabilities in lambda6-85 have been correlated to its rapid folding behavior, and a residual hydrophobic cluster found in R1-69 in 7 M urea has been proposed as a folding initiation site. To understand the early events in the folding of 434 Cro, and for comparison with R1-69 and lambda6-85, we examined the conformational behavior of five peptides covering the entire 434 Cro sequence in water, 40% (by volume) TFE/water, and 7 M urea solutions using CD and NMR. Each peptide corresponds to a helix and adjacent residues as identified in the native 434 Cro NMR and crystal structures. All are soluble and monomeric in the solution conditions examined except for the peptide corresponding to the 434 Cro helix 4, which has low water solubility. Helix formation is observed for the 434 Cro helix 1 and helix 2 peptides in water, for all the peptides in 40% TFE and for none in 7 M urea. NMR data indicate that the helix limits in the peptides are similar to those in the native protein helices. The number of side-chain NOEs in water and TFE correlates with the helix content, and essentially none are observed in 7 M urea for any peptide, except that for helix 5, where a hydrophobic cluster may be present. The low intrinsic folding propensities of the five helices could account for the observed stability and folding behavior of 434 Cro and is, at least qualitatively, in accord with the results of the recently described diffusion-collision model incorporating intrinsic helix propensities.     

α-Helix-forming propensities in peptides and proteins

Much effort has been invested in seeking to understand the thermodynamic basis of helix stability in both peptides and proteins. Recently, several groups have measured the helix-forming propensities of individual residues (Lyu, P. C., Liff, M. I., Marky, L. A., Kallenbach, N. R. Science 250:669–673, 1990; O'Neil, K. T., DeGrado, W. F. Science 250:646–651, 1990; Padmanabhan, S., Marqusee, S., Ridgeway, T., Laue, T. M., Baldwin, R. L. Nature (London) 344:268–270, 1990). Using Monte Carlo computer simulations, we tested the hypothesis that these differences in measured helix-forming propensity are due primarily to loss of side chain conformational entropy upon helix formation (Creamer, T. P., Rose, G. D. Proc. Natl. Acad. Sci. U.S.A. 89:5937–5941, 1992). Our previous study employed a rigid helix backbone, which is here generalized to a completely flexible helix model in order to ensure that earlier results were not a methodological artifact. Using this flexible model, side chain rotamer distributions and entropy losses are calculated and shown to agree with those obtained earlier. We note that the side chain conformational entropy calculated for Trp in our previous study was in error; a corrected value is presented. Extending earlier work, calculated entropy losses are found to correlate strongly with recent helix propensity scales derived from substitutions made within protein helices (Horovitz, A., Matthews, J. M., Fersht, A. R. J. Mol. Biol. 227:560–568, 1992; Blaber, M., Zhang, X.-J., Matthews, B. M. Science 260:1637–1640, 1993). In contrast, little correlation is found between these helix propensity scales and the accessible surface area buried upon formation of a model polyalanyl α-helix. Taken in sum, our results indicate that loss of side chain entropy is a major determinant of the helix-forming tendency of residues in both peptide and protein helices. © 1994 Wiley-Liss, Inc.     

Conformational analysis of peptides corresponding to all the secondary structure elements of protein L B1 domain: secondary structure propensities are not conserved in proteins with the same fold.

The solution conformation of three peptides corresponding to the two beta-hairpins and the alpha-helix of the protein L B1 domain have been analyzed by circular dichroism (CD) and nuclear magnetic resonance spectroscopy (NMR). In aqueous solution, the three peptides show low populations of native and non-native locally folded structures, but no well-defined hairpin or helix structures are formed. In 30% aqueous trifluoroethanol (TFE), the peptide corresponding to the alpha-helix adopts a high populated helical conformation three residues longer than in the protein. The hairpin peptides aggregate in TFE, and no significant conformational change occurs in the NMR observable fraction of molecules. These results indicate that the helical peptide has a significant intrinsic tendency to adopt its native structure and that the hairpin sequences seem to be selected as non-helical. This suggests that these sequences favor the structure finally attained in the protein, but the contribution of the local interactions alone is not enough to drive the formation of a detectable population of native secondary structures. This pattern of secondary structure tendencies is different to those observed in two structurally related proteins: ubiquitin and the protein G B1 domain. The only common feature is a certain propensity of the helical segments to form the native structure. These results indicate that for a protein to fold, there is no need for large native-like secondary structure propensities, although a minimum tendency to avoid non-native structures and to favor native ones could be required.     

Reconstitution of a native-like SH2 domain from disordered peptide fragments examined by multidimensional heteronuclear NMR

The N-terminal SH2 domain from the p85alpha subunit of phosphatidylinositol 3' kinase is cleaved specifically into 9- and 5-kD fragments by limited proteolytic digestion with trypsin. The noncovalent SH2 domain complex and its constituent tryptic peptides have been investigated using high-resolution heteronuclear magnetic resonance (NMR). These studies have established the viability of the SH2 domain as a fragment complementation system. The individual peptide fragments are predominantly unstructured in solution. In contrast, the noncovalent 9-kD + 5-kD complex shows a native-like (1)H-(15)N HSQC spectrum, demonstrating that the two fragments fold into a native-like structure on binding. Chemical shift analysis of the noncovalent complex compared to the native SH2 domain reveals that the highest degree of perturbation in the structure occurs at the cleavage site within a flexible loop and along the hydrophobic interface between the two peptide fragments. Mapping of these chemical shift changes on the structure of the domain reveals changes consistent with the reduction in affinity for the target peptide ligand observed in the noncovalent complex relative to the intact protein. The 5-kD fragment of the homologous Src protein is incapable of structurally complementing the p85 9-kD fragment, either in complex formation or in the context of the full-length protein. These high-resolution structural studies of the SH2 domain fragment complementation features establish the suitability of the system for further protein-folding and design studies.     

Identification and characterization of folding inhibitors of hen egg lysozyme: an example of a new paradigm of drug design

Studies of protein folding indicate the presence of native contacts in the denatured state, giving rise to folding elements which contribute to the accomplishment of the native state. The possibility of finding molecules which can interact with specific folding elements of a target protein preventing it from reaching its native state, and hence from becoming biologically active, is particularly attractive. The notion that folding elements not only provide molecular recognition directing the folding process, but also have conserved sequence, implies that targeting such elements will make protein folding inhibitors less susceptible to mutations which, in many cases, abrogate drug effects. The folding-inhibition strategy can lead to a truly novel and rational approach to drug design, aside from providing new insight into folding. This is illustrated in the case of hen egg lysozyme.     

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.     

Stabilization of myoglobin by multiple alanine substitutions in helical positions.

We have carried out a series of multiple Xaa-->Ala changes at nonadjacent surface positions in the sequence of sperm whale myoglobin. Although the corresponding single substitutions do not increase the thermal stability of the protein, multiple substitutions enhance the stability of the resulting myoglobins. The effect observed is an increase in the observed Tm (midpoint unfolding temperature) relative to that predicted from assuming additivity of the free energy changes corresponding to single mutations. The stabilization occurs in the presence of urea, as measured by the dependence of the unfolding temperature on urea concentration. The sites that have been altered occur in different helices and are not close in sequence or in the native structure of myoglobin. The observed effect is consistent with a role of multiple alanines in residual interactions in the unfolded state of the mutant proteins.     

A comprehensive library of blocked dipeptides reveals intrinsic backbone conformational propensities of unfolded proteins

Despite prolonged scientific efforts to elucidate the intrinsic peptide backbone preferences of amino-acids based on understanding of intermolecular forces, many open questions remain, particularly concerning neighboring peptide interaction effects on the backbone conformational distribution of short peptides and unfolded proteins. Here, we show that spectroscopic studies of a complete library of 400 dipeptides reveal that, irrespective of side-chain properties, the backbone conformation distribution is narrow and they adopt polyproline II and β-strand, indicating the importance of backbone peptide solvation and electronic effects. By directly comparing the dipeptide circular dichroism and NMR results with those of unfolded proteins, the comprehensive dipeptides form a complete set of structural motifs of unfolded proteins. We thus anticipate that the present dipeptide library with spectroscopic data can serve as a useful database for understanding the nature of unfolded protein structures and for further refinements of molecular mechanical parameters.     

Local propensities and statistical potentials of backbone dihedral angles in proteins

The following three issues concerning the backbone dihedral angles of protein structures are presented. (1) How do the dihedral angles of the 20 amino acids depend on the identity and conformation of their nearest residues? (2) To what extent are the native dihedral angles determined by local (dihedral) potentials? (3) How to build a knowledge-based potential for a residue's dihedral angles, considering the identity and conformation of its nearest residues? We find that the dihedral angle distribution for a residue can significantly depend on the identity and conformation of its adjacent residues. These correlations are in sharp contrast to the Flory isolated-pair hypothesis. Statistical potentials are built for all combinations of residue triplets and depend on the dihedral angles between consecutive residues. First, a low-resolution potential is obtained, which only differentiates between the main populated basins in the dihedral angle density plots. Minimization of the dihedral potential for 125 test proteins reveals that most native alpha-helical residues (89%) and a large fraction of native beta-sheet residues (47%) adopt conformations close to their native one. For native loop residues, the percentage is 48%. It is also found that this fraction is higher for residues away from the ends of alpha or beta secondary structure elements. In addition, a higher resolution potential is built as a function of dihedral angles by a smoothing procedure and continuous functions interpolations. Monte Carlo energy minimization with this potential results in a lower fraction for native beta-sheet residues. Nevertheless, because of the higher flexibility and entropy of beta structures, they could be preferred under the influence of non-local interactions. In general, most alpha-helices and many beta-sheets are strongly determined by the local potential, while the conformations in loops and near the end of beta-sheets are more influenced by non-local interactions.     

Folding determinants of disulfide bond forming protein B explored by solution nuclear magnetic resonance spectroscopy

The two-stage model for membrane protein folding postulates that individual helices form first and are subsequently packed against each other. To probe the two-stage model, the structures of peptides representing individual transmembrane helices of the disulfide bond forming protein B have been studied in trifluoroethanol solution as well as in detergent micelles using nuclear magnetic resonance (NMR) and circular dichroism spectroscopy. In TFE solution, peptides showed well-defined α-helical structures. Peptide structures in TFE were compared to the structures of full-length protein obtained by X-ray crystallography and NMR. The extent of α-helical secondary structure coincided well, lending support for the two-stage model for membrane protein folding. However, the conformation of some amino acid side chains differs between the structures of peptide and full-length protein. In micellar solution, the peptides also adopted a helical structure, albeit of reduced definition. Using measurements of paramagnetic relaxation enhancement, peptides were confirmed to be embedded in micelles. These observations may indicate that in the native protein, tertiary interactions additionally stabilize the secondary structure of the individual transmembrane helices.     

Thermodynamics of the reconstitution of tuna cytochrome c from two peptide fragments

Two peptide fragments from tuna cytochrome c (cyt c), N-fragment (residues 1-44 containing the heme) and C-fragment (residues 45-103), combine to form a 1:1 fragment complex. This was clearly proved by ion-spray mass spectrometry. It was found from CD and NMR spectra that the structure of the fragment complex formed is similar to that of an intact cyt c, although each isolated fragment itself is unstructured. Binding constants and enthalpies upon the complex formation were directly observed by isothermal titration calorimetry. Thermodynamic parameters (deltaG(o)b, deltaHb, deltaS(o)b, and deltaC(b)p)) associated with the complex formation were determined at various pHs and temperatures. DeltaHb was found to be almost independent of pH values. The change in heat capacity accompanying the complex formation (deltaC(b)p) was directly determined from the temperature dependence of deltaHb. In addition, the change in heat capacity and enthalpy upon tuna cyt c unfolding were determined by differential scanning calorimetry. Thermodynamic parameters for the unfolding/dissociation process of the fragment complex were compared with those for cyt c unfolding at pH 3.9 and 303 K. In a comparison of two unfolding processes, the heat capacity change of each was very close to the other, while both the unfolding enthalpy and entropy of the fragment complex were larger than those of tuna cyt c. These thermodynamic data suggest that the internal interactions between polar groups (hydrogen bonding) and nonpolar groups (van der Waals interactions) are preserved in the fragment complex as well as in the native state of cyt c.     

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.     

Beyond basins: φ,ψ preferences of a residue depend heavily on the φ,ψ values of its neighbors

The Ramachandran plot distributions of nonglycine residues from experimentally determined structures are routinely described as grouping into one of six major basins: β, P_II, α, α_L, ξ and γ'. Recent work describing the most common conformations adopted by pairs of residues in folded proteins [i.e., (φ,ψ)₂‐motifs] showed that commonly described major basins are not true single thermodynamic basins, but are composed of distinct subregions that are associated with various conformations of either the preceding or following neighbor residue. Here, as documentation of the extent to which the conformational preferences of a central residue are influenced by the conformations of its two neighbors, we present a set of φ,ψ‐plots that are delimited simultaneously by the φ,ψ‐angles of its neighboring residues on both sides. The level of influence seen here is typically greater than the influence associated with considering the identities of neighboring residues, implying that the use of this heretofore untapped information can improve the accuracy of structure prediction algorithms and low resolution protein structure refinement.     

Identifying the structural boundaries of independent folding domains in the alpha subunit of tryptophan synthase, a beta/alpha barrel protein.

Two equilibrium intermediates have previously been observed in the urea denaturation of the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli, an eight-stranded beta/alpha barrel protein. In the current study, a series of amino-terminal fragments were characterized to probe the elementary folding units that may be in part responsible for this complex behavior. Stop-codon mutagenesis was used to produce eight fragments ranging in size from 105-214 residues and containing incremental elements of secondary structure. Equilibrium studies by circular dichroism indicate that all of these fragments are capable of adopting secondary structure. All except for the shortest fragment fold cooperatively. The addition of the fourth, sixth, and eighth beta-strands leads to distinct increases in structure, cooperativity, and/or stability, suggesting that folding involves the modular assembly of betaalphabeta supersecondary structural elements. One-dimensional NMR titrations at high concentrations of urea, probing the environment around His92, were also performed to test for the presence of residual structure in the fragments. All fragments that contained the first four betaalpha units of structure exhibited a cooperative unfolding transition at high concentrations of urea with significant but reduced stability relative to the full-length protein. These results suggest that the residual structure in alphaTS requires the participation of hydrophobic residues in multiple beta-strands that span the entire sequence.     

Effects of alanine substitutions in alpha-helices of sperm whale myoglobin on protein stability.

The peptide backbones in folded native proteins contain distinctive secondary structures, alpha-helices, beta-sheets, and turns, with significant frequency. One question that arises in folding is how the stability of this secondary structure relates to that of the protein as a whole. To address this question, we substituted the alpha-helix-stabilizing alanine side chain at 16 selected sites in the sequence of sperm whale myoglobin, 12 at helical sites on the surface of the protein, and 4 at obviously internal sites. Substitution of alanine for bulky side chains at internal sites destabilizes the protein, as expected if packing interactions are disrupted. Alanine substitutions do not uniformly stabilize the protein, either in capping positions near the ends of helices or at mid-helical sites near the surface of myoglobin. When corrected for the extent of exposure of each side chain replaced by alanine at a mid-helix position, alanine replacement still has no clear effect in stabilizing the native structure. Thus linkage between the stabilization of secondary structure and tertiary structure in myoglobin cannot be demonstrated, probably because of the relatively small free energy differences between side chains in stabilizing isolated helix. By contrast, about 80% of the variance in free energy observed can be accounted for by the loss in buried surface area of the native residue substituted by alanine. The differential free energy of helix stabilization does not account for any additional variation.     

Folding simulations of Trp‐cage mini protein in explicit solvent using biasing potential replica‐exchange molecular dynamics simulations

Replica exchange molecular dynamics (RexMD) simulations are frequently used for studying structure formation and dynamics of peptides and proteins. A significant drawback of standard temperature RexMD is, however, the rapid increase of the replica number with increasing system size to cover a desired temperature range. A recently developed Hamiltonian RexMD method has been used to study folding of the Trp‐cage protein. It employs a biasing potential that lowers the backbone dihedral barriers and promotes peptide backbone transitions along the replica coordinate. In two independent applications of the biasing potential RexMD method including explicit solvent and starting from a completely unfolded structure the formation of near‐native conformations was observed after 30–40 ns simulation time. The conformation representing the most populated cluster at the final simulation stage had a backbone root mean square deviation of ～1.3 Å from the experimental structure. This was achieved with a very modest number of five replicas making it well suited for peptide and protein folding and refinement studies including explicit solvent. In contrast, during five independent continuous 70 ns molecular dynamics simulations formation of collapsed states but no near native structure formation was observed. The simulations predict a largely collapsed state with a significant helical propensity for the helical domain of the Trp‐cage protein already in the unfolded state. Hydrogen bonded bridging water molecules were identified that could play an active role by stabilizing the arrangement of the helical domain with respect to the rest of the chain already in intermediate states of the protein. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

New Fourier transform infrared based computational method for peptide secondary structure determination. II. Application to study of peptide fragments reproducing processing site of ocytocin-neurophysin precursor

A new method for the quantitative determination of the percentage of intramolecular H-bonds, based on Fourier transform infrared techniques, is applied to the conformational analysis of a series of synthetic peptides spanning the processing site of the ocytocin and neurophysin precursor. Even though the method uses traditional tools such as Fourier self-deconvolution, the Nth derivative, and curve-fitting procedures for the analysis of the spectra, the assignment of the absorptions due to peptide groups participating into secondary structures is based on the direct observation and quantification of the isotopic effect induced on the groups participating in intramolecular H-bonds in the presence of organic solvents. This permits the quantification of the different populations of molecules containing intramolecular H-bonds involved in beta-turns and alpha-helices. The results are consistent with those previously obtained by NMR techniques in the same solvent systems.     

Describing protein structure: a general algorithm yielding complete helicoidal parameters and a unique overall axis

We present a general and mathematically rigorous algorithm which allows the helicoidal structure of a protein to be calculated starting from the atomic coordinates of its peptide backbone. This algorithm yields a unique curved axis which quantifies the folding of the backbone and a full set of helicoidal parameters describing the location of each peptide unit. The parameters obtained form a complete and independent set and can therefore be used for analyzing, comparing, or reconstructing protein backbone geometry. This algorithm has been implemented in a computer program named P-Curve. Several examples of its possible applications are discussed.     

蛋白质折迭的不可逆热力学理论及蛋白质动态热力学结构(英)

在机械力学及热力学基础上阐述蛋白质的各种物理性质是从物理学角度理解生物学的基础性工作.详细讨论了蛋白质折迭的不可逆热力学理论,蛋白质动态热力学结构理论.理论推断蛋白质动态热力学变化是一切生物学状态变化的基本热力学状态单位并作为分子生物学变化的分子开关.利用此理论解释了蛋白质的物理学性质及麻醉药的生物物理学机制.     

The Protein Coil Library: a structural database of nonhelix, nonstrand fragments derived from the PDB

Approximately half the structure of folded proteins is either alpha-helix or beta-strand. We have developed a convenient repository of all remaining structure after these two regular secondary structure elements are removed. The Protein Coil Library (http://roselab.jhu.edu/coil/) allows rapid and comprehensive access to non-alpha-helix and non-beta-strand fragments contained in the Protein Data Bank (PDB). The library contains both sequence and structure information together with calculated torsion angles for both the backbone and side chains. Several search options are implemented, including a query function that uses output from popular PDB-culling servers directly. Additionally, several popular searches are stored and updated for immediate access. The library is a useful tool for exploring conformational propensities, turn motifs, and a recent model of the unfolded state.