Evolutionary conservation of the folding nucleus

Here, we present statistical analysis of conservation profiles in families of homologous sequences for nine proteins whose folding nucleus was determined by protein engineering methods. We show that in all but one protein (AcP) folding nucleus residues are significantly more conserved than the rest of the protein. Two aspects of our study are especially important: (i) grouping of amino acid residues into classes according to their physical-chemical properties and (ii) proper normalization of amino acid probabilities that reflects the fact that evolutionary pressure to conserve some amino acid types may itself affect concentration of various amino acid types in protein families. Neglect of any of those two factors may make physical and biological signals from conservation profiles disappear.     

Evolutionary conservation in protein folding kinetics

The sequence and structural conservation of folding transition states have been predicted on theoretical grounds. Using homologous sequence alignments of proteins previously characterized via coupled mutagenesis/kinetics studies, we tested these predictions experimentally. Only one of the six appropriately characterized proteins exhibits a statistically significant correlation between residues' roles in transition state structure and their evolutionary conservation. However, a significant correlation is observed between the contributions of individual sequence positions to the transition state structure across a set of homologous proteins. Thus the structure of the folding transition state ensemble appears to be more highly conserved than the specific interactions that stabilize it.     

Evolutionary aspects of protein structure and folding

Four basic stages of evolution of protein structure are described based on recent work of the authors targeted specifically on reconstruction of the earliest events in the protein evolution. According to this reconstruction, the initial stage of short peptides of, probably, only few amino-acid residues had been followed by formation of closed loops of the size 25-30 residues, which corresponds to the polymer-statistically optimal ring closure size for mixed polypeptide chains. The next stage involved fusion of the respective small linear genes and formation of protein structures consisting of several closed loops of the nearly standard size, up to 4-6 loops (100-200 amino acid residues) in a typical protein fold. The last, modern stage began with combinatorial fusion of the presumably circular 300-600 bp DNA units and, accordingly, formation of multidomain proteins.     

Evolutionary pressures on apicoplast transit peptides

Malaria parasites (species of the genus Plasmodium) harbor a relict chloroplast (the apicoplast) that is the target of novel antimalarials. Numerous nuclear-encoded proteins are translocated into the apicoplast courtesy of a bipartite N-terminal extension. The first component of the bipartite leader resembles a standard signal peptide present at the N-terminus of secreted proteins that enter the endomembrane system. Analysis of the second portion of the bipartite leaders of P. falciparum, the so-called transit peptide, indicates similarities to plant transit peptides, although the amino acid composition of P. falciparum transit peptides shows a strong bias, which we rationalize by the extraordinarily high AT content of P. falciparum DNA. 786 plastid transit peptides were also examined from several other apicomplexan parasites, as well as from angiosperm plants. In each case, amino acid biases were correlated with nucleotide AT content. A comparison of a spectrum of organisms containing primary and secondary plastids also revealed features unique to secondary plastid transit peptides. These unusual features are explained in the context of secondary plastid trafficking via the endomembrane system.     

Evolutionary computer programming of protein folding and structure predictions

In order to understand the mechanism of protein folding and to assist the rational de-novo design of fast-folding, non-aggregating and stable artificial enzymes it is very helpful to be able to simulate protein folding reactions and to predict the structures of proteins and other biomacromolecules. Here, we use a method of computer programming called "evolutionary computer programming" in which a program evolves depending on the evolutionary pressure exerted on the program. In the case of the presented application of this method on a computer program for folding simulations, the evolutionary pressure exerted was towards faster finding deep minima in the energy landscape of protein folding. Already after 20 evolution steps, the evolved program was able to find deep minima in the energy landscape more than 10 times faster than the original program prior to the evolution process.     

Metal ion complexation and folding of linear peptides

A Linear peptide, GASYQDLG was synthesised and used as a model to evaluate the effects of nickel additions to increase the conformational stability. The NMR data obtained for the peptide and its histidyl derivative (H)(3)GASYQDLG(H)(3) suggest that in solution folded structures are present only for the H-tagged peptide-Ni(II) ion system. These results suggest that metal ions and additions of a double histidine tags of suitable length can be used as efficient tools to reduce peptide flexibility without other internal modifications. Synthesis of H-tagged analogs could offer a promising strategy for large-scale preparation of diagnostic tools and, in general, whenever more rigid molecular structures should be advisable.     

Synthesis and folding of native collagen III model peptides.

Solid-phase synthesis of triple-helical peptides, including native collagen III sequences, was started with a trimeric branch, based upon the lysine dipeptide [Fields, C. G., Mickelson, D. J., Drake, S. L., McCarthy, J. B., and Fields, G. B. (1993) J. Biol. Chem. 268, 14153-14160]. Branch synthesis was modified, using TentaGel R as resin, p-hydroxybenzyl alcohol (HMP) as linker, Dde as N(epsilon)-protective group, and HATU/HOAT as coupling reagent. Three homotrimeric sequences, each containing the Gly 606-Gly 618 portion of human type III collagen, were added to the amino functions of the branch. The final incorporation of GlyProHyp triplets, stabilizing the collagen III triple helix, was performed by using protected GlyProHyp tripeptides, each containing tert-butylated hydroxyproline [P(tBu)] instead of hydroxyproline (P). Among the protected tripeptides FmocP(tBu)PG, FmocPP(tBu)G, and FmocGPP(tBu), prepared manually on a chlorotrityl resin, incorporation of FmocPP(tBu)Gly was best suited for synthesis of large and stable peptides, such as PPG(8), containing 8 (PPG)(3) trimers (115 residues, 10 610 Da). The structures of five peptides, differing from each other by the type and number of the triplets incorporated, were verified by MALDI-TOF-MS. Their conformations and thermodynamic data were studied by circular dichroism and differential scanning calorimetry. Except for PPG(8), containing 8 (PPG)(3) trimers with hydroxyproline in the X position and adopting a polyproline II structure, all peptides were triple-helical in 0.1 M acetic acid and their thermal stabilities ranged from t(1/2) = 39. 4 to t(1/2) = 62.5 degrees C, depending on the identity and number of the triplets used. Similar values of the van't Hoff enthalpy, DeltaH(vH), derived from melting curves, and the calorimetric enthalpy, DeltaH(cal), obtained from heat capacity curves, indicate a cooperative ratio of CR = DeltaH(vH)/DeltaH(cal) = 1, establishing a two-state process for unfolding of THP(III) peptides. The independence of the transition temperatures t(1/2) on peptide concentration as well as equilibrium centrifugation data indicate monomolecular dimer(f) to dimer(u) (F(2) <--> U(2)) transitions and, in addition, bimolecular dimer(f) to monomer(u) transitions (F(2) <--> 2U). The dominance of the concentration-independent monomolecular reaction over the concentration-dependent bimolecular reaction makes thermal unfolding of THP(III) peptides appear to be monomolecular. If one designates the molecularity described by the term pseudomonomolecular, unfolding of the dimeric peptides PPG(6-8) follows a two-state, pseudomonomolecular reaction.     

Diversity in the disulfide folding pathways of cystine knot peptides

Summary The plant cyclotides are a fascinating family of circular proteins that contain a cyclic cystine knot motif (CCK). This unique family was discovered only recently but contains over 50 known sequences to date. Various biological activities are associated with these peptides including antimicrobial and insecticidal activity. The knotted topology and cyclic nature of the cyclotides poses interesting questions about the folding mechanisms and how the knotted arrangement of disulfide bonds is formed. Some studies have been performed on related inhibitor cystine knot (ICK) containing peptides, but little is known about the folding mechanisms of CCK molecules. We have examined the oxidative refolding and reductive unfolding of the prototypic member of the cyclotide family, kalata B1. Analysis of the rates of formation of the intermediates along the reductive unfolding pathway highlights the stability conferred by the cystine knot motif. Significant differences are observed between the folding of kalata B1 and an acyclic cystine knot protein, EETI-II, suggesting that the circular backbone has a significant influence in directing the folding pathway.     

Switch-peptides as folding precursors in self-assembling peptides and amyloid fibrillogenesis

The study of conformational transitions of peptides has obtained considerable attention recently because of their importance as a molecular key event in a variety of degenerative diseases. However, the study of peptide self-assembly into beta-sheets and amyloid beta (Abeta) fibrils is strongly hampered by their difficult synthetic access and low solubility. We have recently developed a new concept termed switch-peptides that allows the controlled onset of polypeptide folding and misfolding at physiologic conditions. As a major feature, the folding process is initiated by chemically or enzyme triggered O,N-acyl migration in flexible and soluble folding precursors containing Ser- or Thr-derived switch (S)-elements. The elaborated methodologies are exemplified for the in situ conversion of NPY- and Cyclosporine A-derived prodrugs, as well as for the onset and reversal of alpha and beta conformational transitions in Abeta peptides. In combining orthogonally addressable switch-elements, the consecutive switching on of S-elements gives new insights into the role of individual peptide segments (hot spots) in early processes of polypeptide self-assembly and fibrillogenesis. Finally, the well-known secondary structure disrupting effect of pseudoprolines (PsiPro) is explored for its use as a building block (S-element) in switch-peptides. To this end, synthetic strategies are described, allowing for the preparation of PsiPro-containing folding precursors, exhibiting flexible random-coil conformations devoid of fibril forming propensity. The onset of beta-sheet and fibril formation by restoring the native peptide chain in a single step classify PsiPro-units as the most powerful tool for inhibiting peptide self-assembly, and complement the present methodologies of the switch-concept for the study of fibrillogenesis.     

Molecular dynamics simulations of the folding of poly(alanine) peptides

The secondary structures and the shapes of long-chain polyalanine (PA) molecules were investigated by all-atom molecular dynamics simulations using a modified Amber force field. Homopolymers of polyaminoacids such as PA are convenient models to study the mechanism of protein folding. It was found that the conformational structures of PA peptides are highly sensitive to the chain length. In the absence of solvent, straight α-helices dominate in short (n ∼ 20) peptides at room temperature. A shape transition occurs at a chain length n of 40–45; the compact helix-turn-helix structure (the double-leg hairpin) becomes favored over a straight α-helix. For n = 60, double-leg and the triple-leg hairpins are the only structures present in PA molecules. An exploration of a chain organization in a cubic cavity revealed a clear predisposition of PA molecules for additional breaks in α-helices and the formation of multifolded hairpins. Furthermore, under confinement the hairpin structure becomes much looser, the antiparallel positions of helical stems are disturbed, and a sizeable proportion of the helical stems are transformed from α-helices into 3₁₀-helices.     

Evolutionary similarity among peptide segments is a basis for prediction of protein folding

Short segments of polypeptide, from a protein for which the primary sequence but not the three-dimensional structure is known, are compared to a library of known structures. The basis of comparison is the probability with which residues in the unknown segment might substitute through evolution for residues in segments of known structure. In test cases, segments from known structures that are similar in sequence to those from a protein treated as unknown are often found to be similar in three-dimensional structure to one another and to the true structure of the “unknown” segment. This provides a basis for prediction of the local configuration (secondary structure) of polypeptides.     

Stability of peptide-HLA-I complexes and tapasin folding facilitation--tools to define immunogenic peptides

Only a small fraction of the peptides generated inside the cell end up being presented by HLA-I on the cell surface. High stability of peptide–HLA-I complexes and a low HLA-I tapasin-facilitation have been proposed to predict immunogenicity. We here set out to investigate if these parameters correlated and defined immunogenic peptides. Both peptide–HLA–B¹08:01 and peptide–HLA–A¹02:01 complexes showed small differences in tapasin-facilitation and larger differences in stability. This suggests that the stability of immunogenic peptide–HLA-I complexes vary above an HLA-I allomorph dependent lower limit (e.g. >2 h for HLA–A¹02:01), immunogenicity predicted by tapasin-facilitation may be defined by an equally allomorph unique upper value (e.g. tapasin-facilitation <1.5 for HLA–A¹02:01), and variation above the stability-threshold does not directly reflect a variation in tapasin-facilitation.     

Membrane adsorption, folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

Evolutionary origin and divergence of PQRFamide peptides and LPXRFamide peptides in the RFamide peptide family. Insights from novel lamprey RFamide peptides

Among the RFamide peptide groups, PQRFamide peptides, such as neuropeptide FF (NPFF) and neuropeptide AF (NPAF), share a common C-terminal Pro-Gln-Arg-Phe-NH(2) motif. LPXRFamide (X = L or Q) peptides, such as gonadotropin-inhibitory hormone (GnIH), frog growth hormone-releasing peptide (fGRP), goldfish LPXRFamide peptide and mammalian RFamide-related peptides (RFRPs), also share a C-terminal Leu-Pro-Leu/Gln-Arg-Phe-NH(2) motif. Such a similar C-terminal structure suggests that these two groups may have diverged from a common ancestral gene. In this study, we sought to clarify the evolutionary origin and divergence of these two groups, by identifying novel RFamide peptides from the brain of sea lamprey, one of only two extant groups of the oldest lineage of vertebrates, Agnatha. A novel lamprey RFamide peptide was identified by immunoaffinity purification using the antiserum against LPXRFamide peptide. The lamprey RFamide peptide did not contain a C-terminal LPXRFamide motif, but had the sequence SWGAPAEKFWMRAMPQRFamide (lamprey PQRFa). A cDNA of the precursor encoded one lamprey PQRFa and two related peptides. These related peptides, which also had the C-terminal PQRFamide motif, were further identified as mature endogenous ligands. Phylogenetic analysis revealed that lamprey PQRFamide peptide precursor belongs to the PQRFamide peptide group. In situ hybridization demonstrated that lamprey PQRFamide peptide mRNA is expressed in the regions predicted to be involved in neuroendocrine and behavioral functions. This is the first demonstration of the presence of RFamide peptides in the agnathan brain. Lamprey PQRFamide peptides are considered to have retained the most ancestral features of PQRFamide peptides.     

Analysis of the kinetics of folding of proteins and peptides using circular dichroism

Circular dichroism (CD) is a useful spectroscopic technique for studying the secondary structure, folding and binding properties of proteins. This protocol covers how to use the intrinsic circular dichroic properties of proteins to follow their folding and unfolding as a function of time. Included are methods of obtaining data and for analyzing the folding and unfolding data to determine the rate constants and the order of the folding and unfolding reactions. The protocol focuses on the use of CD to follow folding when it is relatively slow, on the order of minutes to days. The methods for analyzing the data, however, can also be applied to data collected with a CD machine equipped with stopped-flow accessories in the range of milliseconds to seconds and folding analyzed by other spectroscopic methods including changes in absorption or fluorescence spectra as a function of time.     

Antimicrobial and cell-penetrating peptides induce lipid vesicle fusion by folding and aggregation

According to their distinct biological functions, membrane-active peptides are generally classified as antimicrobial (AMP), cell-penetrating (CPP), or fusion peptides (FP). The former two classes are known to have some structural and physicochemical similarities, but fusogenic peptides tend to have rather different features and sequences. Nevertheless, we found that many CPPs and some AMPs exhibit a pronounced fusogenic activity, as measured by a lipid mixing assay with vesicles composed of typical eukaryotic lipids. Compared to the HIV fusion peptide (FP23) as a representative standard, all designer-made peptides showed much higher lipid-mixing activities (MSI-103, MAP, transportan, penetratin, Pep1). Native sequences, on the other hand, were less fusogenic (magainin 2, PGLa, gramicidin S), and pre-aggregated ones were inactive (alamethicin, SAP). The peptide structures were characterized by circular dichroism before and after interacting with the lipid vesicles. A striking correlation between the extent of conformational change and the respective fusion activities was found for the series of peptides investigated here. At the same time, the CD data show that lipid mixing can be triggered by any type of conformation acquired upon binding, whether α-helical, β-stranded, or other. These observations suggest that lipid vesicle fusion can simply be driven by the energy released upon membrane binding, peptide folding, and possibly further aggregation. This comparative study of AMPs, CPPs, and FPs emphasizes the multifunctional aspects of membrane-active peptides, and it suggests that the origin of a peptide (native sequence or designer-made) may be more relevant to define its functional range than any given name.     

Albumin is a redox-active crowding agent that promotes oxidative folding of cysteine-rich peptides

Oxidative folding that occurs in a crowded cellular milieu is characterized by multifaceted interactions that occur among nascent polypeptides and resident components of the endoplasmic reticulum (ER) lumen. Macromolecular crowding has been considered an essential factor in the folding of polypeptides, but the excluded volume effect has not been evaluated for small, disulfide-rich peptides. In the research presented, we examined how macromolecular crowding agents, such as albumin, ovalbumin, and polysaccharides, influenced the kinetics and thermodynamics of forming disulfide bonds in four model peptides of varying molecular size from 13 residues (1.4 kDa) to 58-residues (6.5 kDa): conotoxins: GI, PVIIA, r11a, and bovine pancreatic trypsin inhibitor. Our results indicate that the excluded volume effect does not significantly alter the folding rates nor equilibria for these peptides. In stark contrast, folding reactions were dramatically accelerated, when protein-based crowding agents were present at concentrations lower than those predicted to provide the excluded volume effect. Submillimolar albumin alone was as effective as glutathione in promoting the oxidative folding of GI conotoxin at concentrations typically found in the ER. To the best of our knowledge, this is the first report and quantitative characterization of oxidative folding of peptides mediated by other than thioredoxin-based protein disulfide bonds. Our work raises a possibility that concurrent secretory and ER-resident proteins may influence the oxidative folding of small, cysteine-rich peptides not as crowding agents, but as redox-active factors.     

Electron spin resonance of TOAC labeled peptides: folding transitions and high frequency spectroscopy

The unnatural, conformationally constrained nitroxide amino acid TOAC (2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid) stabilizes helical structure and provides a means for studying rigidly spin labeled peptides by electron spin resonance (ESR). Two new directions in TOAC research are described. The first investigates intermediates formed during alpha-helix unfolding. Double TOAC labeled alpha-helical peptides were unfolded at low temperature in aqueous solution with increasing concentrations of guanidine hydrochloride. Comparison of ESR spectra from two doubly labeled peptides suggests that 3(10)-helix emerges as an intermediate. The second research direction involves the use of high frequency ESR (140 GHz) at low temperature to assess dipolar couplings and, hence, distances between TOAC pairs in a series of 3(10)-helical peptides. Preliminary simulations suggest that high frequency ESR is able to extract correct distances between 6 and 11 A. In addition, the spectra appear to be very sensitive to the relative orientation of the TOAC labels.