Side-chain effects on peptidyl-prolyl cis/trans isomerisation

A library of fluorescent peptides for exploring the substrate specificities of prolyl isomerases

To fully explore the substrate specificities of prolyl isomerases, we synthesized a library of 20 tetrapeptides that are labeled with a 2-aminobenzoyl (Abz) group at the amino terminus and a p-nitroanilide (pNA) group at the carboxy terminus. In this peptide library of the general formula Abz-Ala-Xaa-Pro-Phe-pNA, the position Xaa before the proline is occupied by all 20 proteinogenic amino acids. A conformational analysis of the peptide by molecular dynamics simulations and by NMR spectroscopy showed that the mutual distance between the Abz and pNA moieties in the peptides depends on the isomeric state of the Xaa-Pro bond. In the cis, but not in the trans form, there are significant chemical shift changes of the Abz and pNA moieties, because their aromatic rings are close to each other. This proximity also leads to a strong quenching of Abz fluorescence, which, in combination with a solvent jump, was used to devise a sensitive assay for prolyl isomerases. Unlike the traditional assay, it is not coupled with peptide proteolysis and thus can be employed for protease-sensitive prolyl isomerases as well. The peptide library was used to provide a complete set of P1-site specificities for prototypic human members of the three prolyl isomerase families, FKBP12, cyclophilin 18, and parvulin 14. In a second application, the substrate specificity of SlyD, a protease-sensitive prolyl isomerase from Escherichia coli, was characterized and compared with that of human FKBP12 as well as with homologues from other bacteria.

Chaperone domains convert prolyl isomerases into generic catalysts of protein folding

The cis/trans isomerization of peptide bonds before proline (prolyl bonds) is a rate-limiting step in many protein folding reactions, and it is used to switch between alternate functional states of folded proteins. Several prolyl isomerases of the FK506-binding protein family, such as trigger factor, SlyD, and FkpA, contain chaperone domains and are assumed to assist protein folding in vivo. The prolyl isomerase activity of FK506-binding proteins strongly depends on the nature of residue Xaa of the Xaa-Pro bond. We confirmed this in assays with a library of tetrapeptides in which position Xaa was occupied by all 20 aa. A high sequence specificity seems inconsistent with a generic function of prolyl isomerases in protein folding. Accordingly, we constructed a library of protein variants with all 20 aa at position Xaa before a rate-limiting cis proline and used it to investigate the performance of trigger factor and SlyD as catalysts of proline-limited folding. The efficiencies of both prolyl isomerases were higher than in the tetrapeptide assays, and, intriguingly, this high activity was almost independent of the nature of the residue before the proline. Apparently, the almost indiscriminate binding of the chaperone domain to the refolding protein chain overrides the inherently high sequence specificity of the prolyl isomerase site. The catalytic performance of these folding enzymes is thus determined by generic substrate recognition at the chaperone domain and efficient transfer to the active site in the prolyl isomerase domain.

A cis-proline in alpha-hemoglobin stabilizing protein directs the structural reorganization of alpha-hemoglobin

alpha-Hemoglobin (alphaHb) stabilizing protein (AHSP) is expressed in erythropoietic tissues as an accessory factor in hemoglobin synthesis. AHSP forms a specific complex with alphaHb and suppresses the heme-catalyzed evolution of reactive oxygen species by converting alphaHb to a conformation in which the heme is coordinated at both axial positions by histidine side chains (bis-histidyl coordination). Currently, the detailed mechanism by which AHSP induces structural changes in alphaHb has not been determined. Here, we present x-ray crystallography, NMR spectroscopy, and mutagenesis data that identify, for the first time, the importance of an evolutionarily conserved proline, Pro(30), in loop 1 of AHSP. Mutation of Pro(30) to a variety of residue types results in reduced ability to convert alphaHb. In complex with alphaHb, AHSP Pro(30) adopts a cis-peptidyl conformation and makes contact with the N terminus of helix G in alphaHb. Mutations that stabilize the cis-peptidyl conformation of free AHSP, also enhance the alphaHb conversion activity. These findings suggest that AHSP loop 1 can transmit structural changes to the heme pocket of alphaHb, and, more generally, highlight the importance of cis-peptidyl prolyl residues in defining the conformation of regulatory protein loops.

The folding pathway of onconase is directed by a conserved intermediate

A promising approach to unravel the relationship between sequence information, tertiary structure, and folding mechanism of proteins is the analysis of the folding behavior of proteins with low sequence identity but comparable tertiary structures. Ribonuclease A (RNase A) and its homologues, forming the RNase A superfamily, provide an excellent model system for respective studies. RNase A has been used extensively as a model protein for folding studies. However, little is known about the folding of homologous RNases. Here, we analyze the folding pathway of onconase, a homologous protein from the Northern leopard frog with great potential as a tumor therapeutic, by high-resolution techniques. Although onconase and RNase A significantly differ in the primary structure (28% sequence identity) and in thermodynamic stability (DeltaDeltaG = 20 kJ mol(-1)), both enzymes possess very similar tertiary structures. The present folding studies on onconase by rapid mixing techniques in combination with fluorescence and NMR spectroscopy allow the structural assignment of the three kinetic phases observed in stopped-flow fluorescence spectroscopy. After a slow peptidyl-prolyl cis-to-trans isomerization reaction in the unfolded state, ONC folds via an on-pathway intermediate to the native state. By quenched-flow hydrogen/deuterium exchange experiments coupled with 2D NMR spectroscopy, 31 amino acid residues were identified to be involved in the structure formation of the intermediate. Twelve of these residues are identical in the RNase A sequence, which is a significantly higher percentage (39%) than the overall 28% sequence identity. Moreover, the structure of this intermediate closely resembles two of the intermediates that occur early during the refolding of RNase A. Obviously, in spite of considerable differences in their amino acid sequence the initial folding events of both proteins are comparable, guided by a limited number of conserved residues.

Green fluorescent protein: structure, folding and chromophore maturation

The prolific use of green fluorescent protein and its variants throughout cellular biology relies on the post-translational formation of the chromophore, which proceeds without the need for any additional enzymes or cofactors, except molecular oxygen. In order to form the mature chromophore, the polypeptide backbone must undergo four distinct processes: folding, cyclisation, oxidation and dehydration. This tutorial review looks in detail at the proposed mechanisms for chromophore formation arising out of experimental and computational studies. The folding process is discussed, and the role that the native state plays in catalysing the initial cyclisation and subsequent chemistry is analysed. The specific functions of four conserved residues (Y66, G67, R96 and E222) in the maturation process are also presented. A greater understanding of the maturation process of fluorescent proteins from both jellyfish and coral species will profit the ongoing quest for brighter, faster maturing, genetically-encodable fluorescent probes of all colours, thus increasing their utility throughout the biomedical sciences.

Dependence of the AmII'p proline Raman band on peptide conformation

We utilized UV resonance Raman (UVRR) measurements and density functional theory (DFT) calculations to relate the AmII'p frequency to the psi angle. The AmII'p frequency shifts by approximately 25 cm(-1) as the psi angle is varied over allowed angles of the Pro peptide bond. The AmII'p frequency does not show any significant dependence on the phi dihedral angle. The conformation sensitivity of the AmII'p frequency derives from conformation-induced changes in the planarity of the Pro peptide bond; psi angle changes push the amide nitrogen out of the peptide bond plane. We use this AmII'p frequency dependence on the psi angle to track temperature-induced conformation changes in a polyproline peptide. The temperature-induced 7 cm(-1) downshift in the AmII'p frequency of the polyproline peptide results from an approximately 45 degrees rotation of the psi dihedral angle from psi = 145 degrees (ideal PPII conformation) to psi = 100 degrees (collapsed PPII conformation).

Intramolecular hydrogen bond-controlled prolyl amide isomerization in glucosyl 3'(S)-hydroxy-5'-hydroxymethylproline hybrids: influence of a C-5'-hydroxymethyl substituent on the thermodynamics and kinetics of prolyl amide cis/trans isomerization

Peptide mimics containing spirocyclic glucosyl-(3'-hydroxy-5'-hydroxymethyl)proline hybrids (Glc3'(S)-5'(CH(2)OH)HypHs) with a polar hydroxymethyl substituent at the C-5' position, such as C-terminal ester Ac-Glc3'(S)-5'(CH(2)OH)Hyp-OMe and C-terminal amide Ac-Glc3'(S)-5'(CH(2)OH)Hyp-N'-CH(3), were synthesized. C-Terminal esters exhibit increased cis population (23-53%) relative to Ac-3(S)HyPro-OMe (17%) or Ac-Pro-OMe (14%) in D(2)O. The prolyl amide cis population is further increased to 38-74% in the C-terminal amide form in D(2)O. Our study shows that the stereochemistry of the hydroxymethyl substituent at the C-5' position of proline permits tuning of the prolyl amide cis/trans isomer ratio. Inversion-magnetization transfer NMR experiments indicate that the stereochemistry of the hydroxymethyl substituent has a dramatic effect on the kinetics of prolyl amide cis/trans isomerization. A 200-fold difference in the trans-to-cis (k(tc)) isomerization and a 90-fold rate difference in the cis-to-trans (k(ct)) isomerization is observed between epimeric C-5' 3 and 4. When compared to reference peptide mimics Ac-Pro-OMe and Ac-3(S)Hyp-OMe, our study demonstrates that a (13-16)-fold decrease in k(tc) and k(ct) is observed for the C-5'(S), while a (5-24)-fold acceleration is observed for the C-5'(R) epimer. DFT calculations indicate that the pyrrolidine ring prefers a C(beta) exo pucker in both Ac-Glc3'(S)-5'(CH(2)OH)Hyp-OMe diastereoisomers. Computational calculations and chemical shift temperature coefficient (Delta delta/Delta T) experiments indicate that the hydroxymethyl group at C-5' in Ac-Glc3'(S)-5'(CH(2)OH)Hyp-OMe forms a stabilizing intramolecular hydrogen bond to the carbonyl of the N-acetyl group in both epimeric cis isomers. However, a competing intramolecular hydrogen bond between the hydroxymethyl groups in the pyrrolidine ring and pyran ring stabilizes the trans isomer in the C-5'(S) diastereoisomer. The dramatic differences in the kinetic properties of the diastereoisomeric peptide mimics are rationalized by the presence or absence of an intramolecular hydrogen bond between the hydroxymethyl substituent located at C-5' and the developing lone pair on the nitrogen atom of the N-acetyl group in the transition state.

Noncooperative folding of subdomains in adenylate kinase

Conformational change is regulating the biological activity of a large number of proteins and enzymes. Efforts in structural biology have provided molecular descriptions of the interactions that stabilize the stable ground states on the reaction trajectories during conformational change. Less is known about equilibrium thermodynamic stabilities of the polypeptide segments that participate in structural changes and whether the stabilities are relevant for the reaction pathway. Adenylate kinase (Adk) is composed of three subdomains: CORE, ATPlid, and AMPbd. ATPlid and AMPbd are flexible nucleotide binding subdomains where large-scale conformational changes are directly coupled to catalytic activity. In this report, the equilibrium thermodynamic stabilities of Adk from both mesophilic and hyperthermophilic bacteria were investigated using solution state NMR spectroscopy together with protein engineering experiments. Equilibrium hydrogen to deuterium exchange experiments indicate that the flexible subdomains are of significantly lower thermodynamic stability compared to the CORE subdomain. Using site-directed mutagenesis, parts of ATPlid and AMPbd could be selectively unfolded as a result of perturbation of hydrophobic clusters located in these respective subdomains. Analysis of the perturbed Adk variants using NMR spin relaxation and C(alpha) chemical shifts shows that the CORE subdomain can fold independently of ATPlid and AMPbd; consequently, folding of the two flexible subdomains occurs independently of each other. Based on the experimental results it is apparent that the flexible subdomains fold into their native structure in a noncooperative manner with respect to the CORE subdomain. These results are discussed in light of the catalytically relevant conformational change of ATPlid and AMPbd.

3S-fluoroproline as a probe to monitor proline isomerization during protein folding by 19F-NMR

Variable-temperature inversion transfer NMR is used to determine the kinetic and thermodynamic parameters of cis-trans isomerization of N-Ac-(3R) and (3S)-fluoroproline-OMe.

Molecular determinants of a native-state prolyl isomerization

Prolyl cis/trans isomerizations determine the rates of many protein-folding reactions, and they can serve as molecular switches and timers. The energy required to shift the prolyl cis/trans equilibrium during these processes originates from conformational reactions that are linked structurally and energetically with prolyl isomerization. We used the N2 domain of the gene-3-protein of phage fd to elucidate how such an energetic linkage develops in the course of folding. The Asp160-Pro161 bond at the tip of a beta hairpin of N2 is cis in the crystal structure, but in fact, it exists as a mixture of conformers in folded N2. During refolding, about 10 kJ mol(-1) of conformational energy becomes available for a 75-fold shift of the cis/trans equilibrium constant at Pro161, from 7/93 in the unfolded to 90/10 in the folded form. We combined single- and double-mixing kinetic experiments with a mutational analysis to identify the structural origin of this proline shift energy and to elucidate the molecular path for the transfer of this energy to Pro161. It originates largely, if not entirely, from the two-stranded beta sheet at the base of the Pro161 hairpin. The two strands improve their stabilizing interactions when Pro161 is cis, and this stabilization is propagated to Pro161, because the connector peptides between the beta strands and Pro161 are native-like folded when Pro161 is cis. In the presence of a trans-Pro161, the connector peptides are locally unfolded, and thus, Pro161 is structurally and energetically uncoupled from the beta sheet. Such interrelations between local folding and prolyl isomerization and the potential modulation by prolyl isomerases might also be used to break and reestablish slow communication pathways in proteins.

Support vector machines for prediction of peptidyl prolyl cis/trans isomerization

A new method for peptidyl prolyl cis/trans isomerization prediction based on the theory of support vector machines (SVM) was introduced. The SVM represents a new approach to supervised pattern classification and has been successfully applied to a wide range of pattern recognition problems. In this study, six training datasets consisting of different length local sequence respectively were used. The polynomial kernel functions with different parameter d were chosen. The test for the independent testing dataset and the jackknife test were both carried out. When the local sequence length was 20-residue and the parameter d = 8, the SVM method archived the best performance with the correct rate for the cis and trans forms reaching 70.4 and 69.7% for the independent testing dataset, 76.7 and 76.6% for the jackknife test, respectively. Matthew's correlation coefficients for the jackknife test could reach about 0.5. The results obtained through this study indicated that the SVM method would become a powerful tool for predicting peptidyl prolyl cis/trans isomerization.

Thermodynamics, kinetics, and salt dependence of folding of YopM, a large leucine-rich repeat protein

Small globular proteins have many contacts between residues that are distant in primary sequence. These contacts create a complex network between sequence-distant segments of secondary structure, which may be expected to promote the cooperative folding of globular proteins. Although repeat proteins, which are composed of tandem modular units, lack sequence-distant contacts, several of considerable length have been shown to undergo cooperative two-state folding. To explore the limits of cooperativity in repeat proteins, we have studied the unfolding of YopM, a leucine-rich repeat (LRR) protein of over 400 residues. Despite its large size and modular architecture (15 repeats), YopM equilibrium unfolding is highly cooperative, and shows a very strong dependence on the concentration of urea. In contrast, kinetic studies of YopM folding indicate a mechanism that includes one or more transient intermediates. The urea dependence of the folding and unfolding rates suggests a relatively small transition state ensemble. As with the urea dependence, we have found an extreme dependence of the free energy of unfolding on the concentration of salt. This salt dependence likely results from general screening of a large number of unfavorable columbic interactions in the folded state, rather than from specific cation binding.

The dual-basin landscape in GFP folding

Recent experimental studies suggest that the mature GFP has an unconventional landscape composed of an early folding event with a typical funneled landscape, followed by a very slow search and rearrangement step into the locked, active chromophore-containing structure. As we have shown previously, the substantial difference in time scales is what generates the observed hysteresis in thermodynamic folding. The interconversion between locked and the soft folding structures at intermediate denaturant concentrations is so slow that it is not observed under the typical experimental observation time. Simulations of a coarse-grained model were used to describe the fast folding event as well as identify native-like intermediates on energy landscapes enroute to the fluorescent native fold. Interestingly, these simulations reveal structural features of the slow dynamic transition to chromophore activation. Experimental evidence presented here shows that the trapped, native-like intermediate has structural heterogeneity in residues previously linked to chromophore formation. We propose that the final step of GFP folding is a "locking" mechanism leading to chromophore formation and high stability. The combination of previous experimental work and current simulation work is explained in the context of a dual-basin folding mechanism described above.

Kinetic traps in the folding of beta alpha-repeat proteins: CheY initially misfolds before accessing the native conformation

The beta alpha-repeat class of proteins, represented by the (beta alpha)(8) barrel and the alpha/beta/alpha sandwich, are among the most common structural platforms in biology. Previous studies on the folding mechanisms of these motifs have revealed or suggested that the initial event involves the submillisecond formation of a kinetically trapped species that must at least partially unfold before productive folding to the respective native conformation can occur. To test the generality of these observations, CheY, a bacterial response regulator, was subjected to an extensive analysis of its folding reactions. Although earlier studies had proposed the formation of an off-pathway intermediate, the data available were not sufficient to rule out an alternative on-pathway mechanism. A global analysis of single- and double-jump kinetic data, combined with equilibrium unfolding data, was used to show that CheY folds and unfolds through two parallel channels defined by the state of isomerization of a prolyl peptide bond in the active site. Each channel involves a stable, highly structured folding intermediate whose kinetic properties are better described as the properties of an off-pathway species. Both intermediates subsequently flow through the unfolded state ensemble and adopt the native cis-prolyl isomer prior to forming the native state. Initial collapse to off-pathway folding intermediates is a common feature of the folding mechanisms of beta alpha-repeat proteins, perhaps reflecting the favored partitioning to locally determined substructures that cannot directly access the native conformation. Productive folding requires the dissipation of these prematurely folded substructures as a prelude to forming the larger-scale transition state that leads to the native conformation. Results from Gō-modeling studies in the accompanying paper elaborate on the topological frustration in the folding free-energy landscape of CheY.

The periplasmic peptidyl prolyl cis-trans isomerases PpiD and SurA have partially overlapping substrate specificities

One of the rate-limiting steps in protein folding has been shown to be the cis-trans isomerization of proline residues, catalysed by a range of peptidyl prolyl cis-trans isomerases (PPIases). In the periplasmic space of Escherichia coli and other Gram-negative bacteria, two PPIases, SurA and PpiD, have been identified, which show high sequence similarity to the catalytic domain of the small PPIase parvulin. This observation raises a question regarding the biological significance of two apparently similar enzymes present in the same cellular compartment: do they interact with different substrates or do they catalyse different reactions? The substrate-binding motif of PpiD has not been characterized so far, and no biochemical data were available on how this folding catalyst recognizes and interacts with substrates. To characterize the interaction between model peptides and the periplasmic PPIase PpiD from E. coli, we employed a chemical crosslinking strategy that has been used previously to elucidate the interaction of substrates with SurA. We found that PpiD interacted with a range of model peptides independently of whether they contained proline residues or not. We further demonstrate here that PpiD and SurA interact with similar model peptides, and therefore must have partially overlapping substrate specificities. However, the binding motif of PpiD appears to be less specific than that of SurA, indicating that the two PPIases might interact with different substrates. We therefore propose that, although PpiD and SurA have partially overlapping substrate specificities, they fulfil different functions in the cell.

Modulation of 14-3-3 interaction with phosphorylated histone H3 by combinatorial modification patterns

Post-translational modifications of histones are determining factors in the global and local regulation of genome activity. Phosphorylation of histone H3 is globally associated with mitotic chromatin compaction but occurs in a much more restricted manner during interphase transcriptional regulation of a limited subset of genes. In the course of gene regulation, serine 10 phosphorylation at histone H3 is targeted to a very small fraction of nucleosomes that is highly susceptible to additional acetylation events. Recently, we and others have identified 14-3-3 as a binding protein that recognizes both phosphorylated serine 10 and phosphorylated serine 28 on histone H3. In vitro, the affinity of 14-3-3 for phosphoserine 10 is weak but becomes significantly increased by additional acetylation of either lysine 9 or lysine 14 on the same histone tail. In contrast, the histone H3S28 site matches elements of 14-3-3 high affinity consensus motifs. This region mediates an initial stronger interaction that is less susceptible to modulation by "auxiliary" modifications. Here we discuss the binding of 14-3-3 proteins to histone H3 in detail and putative biological implications of these interactions.

Universal kinetics of helix-coil transition in gelatin

By covering a much wider concentration range than previous studies we find a very unusual exponential dependence of the rate of helix formation on concentration of gelatin in water and ethylene glycol solutions. By applying a procedure of concentration-temperature superposition we build a master curve describing the initial renaturation rates in both solvents. The growth of the normalized helical fraction chi(t) is a first-order process, with a rate constant consistent with cis-trans isomerization, in most situations. We propose that association of three separate chains to form a triple helical nucleus occurs rapidly and contributes less to the helical onset than previously thought. The measured helix content is a result of lengthening of the triple helix after nucleation, by zipping from the associated nuclei.

Protected aminooxyprolines for expedited library synthesis: application to Tsg101-directed proline-oxime containing peptides

The stereoselective synthesis of aminooxy-containing proline analogues bearing Fmoc/Boc or Fmoc/Mtt protection that renders them suitable for incorporation into peptides using Fmoc protocols is reported. Acid-catalyzed unmasking at the completion of peptide synthesis yields free aminooxy-functionalities for oxime formation through reaction with libraries of aldehydes. This allows post solid-phase diversification strategies that may facilitate structure-activity relationship studies.

Probing polyproline structure and dynamics by photoinduced electron transfer provides evidence for deviations from a regular polyproline type II helix

Polyprolines are well known for adopting a regular polyproline type II helix in aqueous solution, rendering them a popular standard as molecular ruler in structural molecular biology. However, single-molecule spectroscopy studies based on Förster resonance energy transfer (FRET) have revealed deviations of experimentally observed end-to-end distances of polyprolines from theoretical predictions, and it was proposed that the discrepancy resulted from dynamic flexibility of the polyproline helix. Here, we probe end-to-end distances and conformational dynamics of poly-l-prolines with 1-10 residues using fluorescence quenching by photoinduced-electron transfer (PET). A single fluorophore and a tryptophan residue, introduced at the termini of polyproline peptides, serve as sensitive probes for distance changes on the subnanometer length scale. Using a combination of ensemble fluorescence and fluorescence correlation spectroscopy, we demonstrate that polyproline samples exhibit static structural heterogeneity with subpopulations of distinct end-to-end distances that do not interconvert on time scales from nano- to milliseconds. By observing prolyl isomerization through changes in PET quenching interactions, we provide experimental evidence that the observed heterogeneity can be explained by interspersed cis isomers. Computer simulations elucidate the influence of trans/cis isomerization on polyproline structures in terms of end-to-end distance and provide a structural justification for the experimentally observed effects. Our results demonstrate that structural heterogeneity inherent in polyprolines, which to date are commonly applied as a molecular ruler, disqualifies them as appropriate tool for an accurate determination of absolute distances at a molecular scale.

Topology and Sequence in the Folding of a TIM Barrel Protein: Global Analysis Highlights Partitioning between Transient Off-pathway and Stable On-pathway Folding Intermediates in the Complex Folding Mechanism of a (βα)8 Barrel of Unknown Function from B. subtilis

The relative contributions of chain topology and amino acid sequence in directing the folding of a (betaalpha)(8) TIM barrel protein of unknown function encoded by the Bacillus subtilis iolI gene (IOLI) were assessed by reversible urea denaturation and a combination of circular dichroism, fluorescence and time-resolved fluorescence anisotropy spectroscopy. The equilibrium reaction for IOLI involves, in addition to the native and unfolded species, a stable intermediate with significant secondary structure and stability and self-associated forms of both the native and intermediate states. Global kinetic analysis revealed that the unfolded state partitions between an off-pathway refolding intermediate and the on-pathway equilibrium intermediate early in folding. Comparisons with the folding mechanisms of two other TIM barrel proteins, indole-3-glycerol phosphate synthase from the thermophile Sulfolobus solfataricus (sIGPS) and the alpha subunit of Escherichia coli tryptophan synthase (alphaTS), reveal striking similarities that argue for a dominant role of the topology in both early and late events in folding. Sequence-specific effects are apparent in the magnitudes of the relaxation times and relative stabilities, in the presence of additional monomeric folding intermediates for alphaTS and sIGPS and in rate-limiting proline isomerization reactions for alphaTS.

Enhanced stability of cis Pro-Pro peptide bond in Pro-Pro-Phe sequence motif

Identification of sequence motifs that favor cis peptide bonds in proteins is important for understanding and designing proteins containing turns mediated by cis peptide conformations. From (1)H NMR solution studies on short peptides, we show that the Pro-Pro peptide bond in Pro-Pro-Phe almost equally populates the cis and trans isomers, with the cis isomer stabilized by a CHc...pi interaction involving the terminal Pro and Phe. We also show that Phe is over-represented at sequence positions immediately following cis Pro-Pro motifs in known protein structures. Our results demonstrate that the Pro-Pro cis conformer in Pro-Pro-Phe sequence motifs is as important as the trans conformer, both in short peptides as well as in natively folded proteins.

Dynamics, Hydration, and Motional Averaging of a Loop-Gated Artificial Protein Cavity: The W191G Mutant of Cytochrome c Peroxidase in Water as Revealed by Molecular Dynamics Simulations

Five molecular dynamics simulations of the W191G cavity mutant of cytochrome c peroxidase in explicit water reveal distinct dynamic and hydration behavior depending on the closed or open state of the flexible loop gating the cavity, the binding of (K+ or small molecule) cations, and the system temperature. The conformational spaces sampled by the loop region and by the cavity significantly reduce upon binding. The largest ordering factor on water dynamics is the presence of the K+ ion occupying the gated cavity. Considerable water exchange occurs for the open-gate cavity when no ligand or cation is bound. In all cases, good correspondence is found between the calculated (ensemble-averaged) location of water molecules and the water sites determined by X-ray crystallography experiments. However, our simulations suggest that these sites do not necessarily correspond to the presence of bound water molecules. In fact, individual water molecules may repeatedly exchange within the cavity volume yet occupy on average these water sites. Four major conclusions emerge. First, it seems misleading to interpret the conformation of protein loop regions in terms of single dominant structures. Second, our simulations support the general picture of Pro 190 cis-trans isomerization as a determinant of the loop-opening mechanism. Third, receptor flexibility is fundamental for ligand binding and molecular recognition, and our results suggest its importance for the docking of small compounds to the artificial cavity. Fourth, after validation against the available experimental data, molecular dynamics simulations can be used to characterize the dynamics and exchange of water molecules and ions, providing atomic level and time-dependent information otherwise inaccessible to experiments.

The rough energy landscape of superfolder GFP is linked to the chromophore

Many green fluorescent protein (GFP) variants have been developed for use as fluorescent tags, and recently a superfolder GFP (sfGFP) has been developed as a robust folding reporter. This new variant shows increased stability and improved folding kinetics, as well as 100% recovery of native protein after denaturation. Here, we characterize sfGFP, and find that this variant exhibits hysteresis as unfolding and refolding equilibrium titration curves are non-coincident even after equilibration for more than eight half-lives as estimated from kinetic unfolding and refolding studies. This hysteresis is attributed to trapping in a native-like intermediate state. Mutational studies directed towards inhibiting chromophore formation indicate that the novel backbone cyclization is responsible for the hysteresis observed in equilibrium titrations of sfGFP. Slow equilibration and the presence of intermediates imply a rough landscape. However, de novo folding in the absence of the chromophore is dominated by a smoother energy landscape than that sampled during unfolding and refolding of the post-translationally modified polypeptide.

Insight into the role of hydration on protein dynamics

The potential energy surface of a protein is rough. This intrinsic energetic roughness affects diffusion, and hence the kinetics. The dynamics of a system undergoing Brownian motion on this surface in an implicit continuum solvent simulation can be tuned via the frictional drag or collision frequency to be comparable to that of experiments or explicit solvent simulations. We show that the kinetic rate constant for a local rotational isomerization in stochastic simulations with continuum solvent and a collision frequency of 2 ps(-1) is about 10(4) times faster than that in explicit water and experiments. A further increase in the collision frequency to 60 ps(-1) slows down the dynamics, but does not fully compensate for the lack of explicit water. We also show that the addition of explicit water does not only slow down the dynamics by increasing the frictional drag, but also increases the local energetic roughness of the energy landscape by as much as 1.0 kcal/mol.

A 10-Å Spectroscopic Ruler Applied to Short Polyprolines

Fluorescence resonance energy transfer (FRET) from the amino acid tryptophan (Trp) as donor and a 2,3-diazabicyclo[2.2.2]oct-2-ene-labeled asparagine (Dbo) as acceptor in peptides of the general structure Trp-(Pro)n-Dbo-NH2 (n = 1-6) was investigated by steady-state and time-resolved fluorescence, CD, and NMR spectroscopy as well as by molecular dynamics (MD) simulations (GROMOS96 force field). The Trp/Dbo FRET pair is characterized by a very short Förster radius (R0 ca. 9 A), which allowed distance determinations in such short peptides. Water and propylene glycol were investigated as solvents. The peptides were designed to show an early nucleation of the poly(Pro)II (PPII) secondary helix structure for n > or = 2, which was confirmed by their CD spectra. The shortest peptide (n = 1) adopts preferentially the trans conformation about the Trp-Pro bond, as confirmed by NMR spectra. The FRET efficiencies ranged 2-72% and were found to depend sensitively on the peptide length, i.e., the number of intervening proline residues. The analysis of the FRET data at different levels of theory (assuming either a fixed distance or distance distributions according to a wormlike chain or Gaussian model) afforded donor-acceptor distances between ca. 8 A (n = 1) and ca. 16 A (n = 6) in water, which were found to be similar or slightly higher in propylene glycol. The distances afforded by the Trp/Dbo FRET pair were found to be reasonable in comparison to literature data, expectations from the PPII helix structure, and the results from MD simulations. The persistence lengths for the longer peptides were found to lie at 30-70 A in water and 220 +/- 40 A in propylene glycol, suggesting a more rigid PPII helical structure in propylene glycol. A detailed comparison with literature data on FRET in polyprolines demonstrates that the donor-acceptor distances extracted by FRET are correlated with the Förster radii of the employed FRET pairs. This demonstrates the limitations of using FRET as a spectroscopic ruler for short polyprolines, which is presumably due to the breakdown of the point dipole approximation in Förster theory, when the size of the chromophores becomes comparable or larger than the distances under investigation.

Geometry of nonbonded interactions involving planar groups in proteins

Although hydrophobic interaction is the main contributing factor to the stability of the protein fold, the specificity of the folding process depends on many directional interactions. An analysis has been carried out on the geometry of interaction between planar moieties of ten side chains (Phe, Tyr, Trp, His, Arg, Pro, Asp, Glu, Asn and Gln), the aromatic residues and the sulfide planes (of Met and cystine), and the aromatic residues and the peptide planes within the protein tertiary structures available in the Protein Data Bank. The occurrence of hydrogen bonds and other nonconventional interactions such as C-H...pi, C-H...O, electrophile-nucleophile interactions involving the planar moieties has been elucidated. The specific nature of the interactions constraints many of the residue pairs to occur with a fixed sequence difference, maintaining a sequential order, when located in secondary structural elements, such as alpha-helices and beta-turns. The importance of many of these interactions (for example, aromatic residues interacting with Pro or cystine sulfur atom) is revealed by the higher degree of conservation observed for them in protein structures and binding regions. The planar residues are well represented in the active sites, and the geometry of their interactions does not deviate from the general distribution. The geometrical relationship between interacting residues provides valuable insights into the process of protein folding and would be useful for the design of protein molecules and modulation of their binding properties.

Loop anchor modification causes the population of an alternative native state in an SH3-like domain

Many stably folded proteins are proposed to contain long, unstructured loops. A series of hybrid proteins (EbE1-4) containing the folded scaffold of photosystem I accessory protein E (PsaE), an SH3-like protein, and the 40-residue heme-binding loop of cytochrome b(5) was created to inspect the dependence of thermodynamic and kinetic parameters on the residues at the interface of folded and flexible regions. Compared to the simplest hybrid (EbE1), the chimeras differed by Gly insertions (EbE2, EbE3) or an asymmetric four-residue restructuring of loop termini (EbE4). NMR spectroscopy indicated that the chimeras retained the PsaE topology; native and unfolded state solubilities, however, were affected to varying degrees. Thermal and chemical denaturation experiments revealed that the EbE2 and EbE1 constructs resulted in a modest destabilization of the PsaE core, whereas apparent stability was increased by >5 kJ/mol in EbE4. EbE3 aggregated at microM concentrations and was not studied in detail. EbE4 populated two native states (N1 and N2), which differed by hydrophobic core packing and C-terminal interactions. At room temperature, the population ratio ( approximately 3-4:1) favored the state whose spectroscopic properties most resembled those of PsaE (N1). EbE4 also demonstrated altered folding kinetics, displaying multiple slow phases related to the population of intermediates and possibly N2. It was concluded that loop anchors can affect protein properties, including stability, via short-range effects on local structure and long-range communication with the packed hydrophobic core. Modification of the attachment points appears to be a possible stepping stone in the transition from one three-dimensional structure to another.

Peptidyl-prolyl cis/trans isomerases and transcription: is there a twist in the tail?

Eukaryotic transcription is regulated predominantly by the post-translational modification of the participating components. One such modification is the cis-trans isomerization of peptidyl-prolyl bonds, which results in a conformational change in the protein involved. Enzymes that carry out this reaction include the yeast peptidyl-prolyl cis/trans isomerase Ess1 and its human counterpart Pin1, both of which recognize phosphorylated target motifs exclusively. Consequently, they operate together with proline-directed serine-threonine kinases and phosphatases. High-profile client proteins involved in transcription include steroid hormone receptors, cell-cycle regulators and immune mediators. Other key targets are elements of the transcription machinery, including the multiply phosphorylated carboxy-terminal domain of RNA polymerase II. Changes in isomerase activity have been shown to alter the transactivation potential, protein stability or intracellular localization of these client proteins. The resulting disruption to developmental processes and cell proliferation has been linked, in some cases, to human cancers.

Thermodynamic effects of proline introduction on protein stability

The amino acid Pro is more rigid than other naturally occurring amino acids and, in proteins, lacks an amide hydrogen. To understand the structural and thermodynamic effects of Pro substitutions, it was introduced at 13 different positions in four different proteins, leucine-isoleucine-valine binding protein, maltose binding protein, ribose binding protein, and thioredoxin. Three of the maltose binding protein mutants were characterized by X-ray crystallography to confirm that no structural changes had occurred upon mutation. In the remaining cases, fluorescence and CD spectroscopy were used to show the absence of structural change. Stabilities of wild type and mutant proteins were characterized by chemical denaturation at neutral pH and by differential scanning calorimetry as a function of pH. The mutants did not show enhanced stability with respect to chemical denaturation at room temperature. However, 6 of the 13 single mutants showed a small but significant increase in the free energy of thermal unfolding in the range of 0.3-2.4 kcal/mol, 2 mutants showed no change, and 5 were destabilized. In five of the six cases, the stabilization was because of reduced entropy of unfolding. However, the magnitude of the reduction in entropy of unfolding was typically several fold larger than the theoretical estimate of -4 cal K(-1) mol(-1) derived from the relative areas in the Ramachandran map accessible to Pro and Ala residues, respectively. Two double mutants were constructed. In both cases, the effects of the single mutations on the free energy of thermal unfolding were nonadditive.

Sequence-structure and structure-function analysis in cysteine-rich domains forming the ultrastable nematocyst wall

The nematocyst wall of cnidarians is a unique biomaterial that withstands extreme osmotic pressures, allowing an ultrafast discharge of the nematocyst capsules. Assembly of the highly robust nematocyst wall is achieved by covalent linkage of cysteine-rich domains (CRDs) from two main protein components, minicollagens and nematocyst outer wall antigen (NOWA). The bipolar minicollagens have different disulfide patterns and topologies in their N and C-terminal CRDs. The functional significance of this polarity has been elusive. Here, we show by NMR structural analysis that all representative cysteine-rich domains of NOWA are structurally related to N-terminal minicollagen domains. Natural sequence insertions in NOWA CRDs have very little effect on the tightly knit domain structures, nor do they preclude the efficient folding to a single native conformation. The different folds in NOWA CRDs and the atypical C-terminal minicollagen domain on the other hand can be directly related to different conformational preferences in the reduced states. Ultrastructural analysis in conjunction with aggregation studies argues for an association between the similar NOWA and N-terminal minicollagen domains in early stages of the nematocyst wall assembly, which is followed by the controlled association between the unusual structures of C-terminal minicollagen domains.

Influence of the Internal Disulfide Bridge on the Folding Pathway of the CL Antibody Domain

Disulfide bridges are one of the most important factors stabilizing the native structure of a protein. Whereas the basis for their stabilizing effect is well understood, their role in a protein folding reaction still seems to require further attention. We used the constant domain of the antibody light chain (C(L)), a representative of the ubiquitous immunoglobulin (Ig)-superfamily, to delineate the kinetic role of its single buried disulfide bridge. Independent of its redox state, the monomeric C(L) domain adopts a typical Ig-fold under native conditions and does not retain significant structural elements when unfolded. Interestingly, its folding pathway is strongly influenced by the disulfide bridge. The more stable oxidized protein folds via a highly structured on-pathway intermediate, whereas the destabilized reduced protein populates a misfolded off-pathway species on its way to the native state. In both cases, the formation of the intermediate species is shown to be independent of the isomerization state of the Tyr(141)-Pro(142) bond. Our results demonstrate that the internal disulfide bridge in an antibody domain restricts the folding pathway by bringing residues of the folding nucleus into proximity thus facilitating the way to the native state.

Structure and mechanism of the phosphotyrosyl phosphatase activator

Phosphotyrosyl phosphatase activator (PTPA), also known as PP2A phosphatase activator, is a conserved protein from yeast to human. Here we report the 1.9 A crystal structure of human PTPA, which reveals a previously unreported fold consisting of three subdomains: core, lid, and linker. Structural analysis uncovers a highly conserved surface patch, which borders the three subdomains, and an associated deep pocket located between the core and the linker subdomains. The conserved surface patch and the deep pocket are responsible for binding to PP2A and ATP, respectively. PTPA and PP2A A-C dimer together constitute a composite ATPase. PTPA binding to PP2A results in a dramatic alteration of substrate specificity, with enhanced phosphotyrosine phosphatase activity and decreased phosphoserine phosphatase activity. This function of PTPA strictly depends on the composite ATPase activity. These observations reveal significant insights into the function and mechanism of PTPA and have important ramifications for understanding PP2A function.

Population of On-pathway Intermediates in the Folding of Ubiquitin

The role that intermediate states play in protein folding is the subject of intense investigation and in the case of ubiquitin has been controversial. We present fluorescence-detected kinetic data derived from single and double mixing stopped-flow experiments to show that the F45W mutant of ubiquitin (WT*), a well-studied single-domain protein and most recently regarded as a simple two-state system, folds via on-pathway intermediates. To account for the discrepancy we observe between equilibrium and kinetic stabilities and m-values, we show that the polypeptide chain undergoes rapid collapse to an intermediate whose presence we infer from a fast lag phase in interrupted refolding experiments. Double-jump kinetic experiments identify two direct folding phases that are not associated with slow isomerisation reactions in the unfolded state. These two phases are explained by kinetic partitioning which allows molecules to reach the native state from the collapsed state via two possible competing routes, which we further examine using two destabilised ubiquitin mutants. Interrupted refolding experiments allow us to observe the formation and decay of an intermediate along one of these pathways. A plausible model for the folding pathway of ubiquitin is presented that demonstrates that obligatory intermediates and/or chain collapse are important events in restricting the conformational search for the native state of ubiquitin.

Impact of cis-proline analogs on peptide conformation

The beta-turn is a common motif in both proteins and peptides and often a recognition site in protein interactions. A beta-turn of four sequential residues reverses the direction of the peptide chain and is classified by the phi and psi backbone torsional angles of residues i + 1 and i + 2. The type VI turn usually contains a proline with a cis-amide bond at residue i + 2. Cis-proline analogs that constrain the peptide to adopt a type VI turn led to peptidomimetics with enhanced activity or metabolic stability. To compare the impact of different analogs on amide cis-trans isomerism and peptide conformation, the conformational preference for the cis-amide bond and the type VI turn was investigated at the MP2/6-31+G** level of theory in water (polarizable continuum water model). Analogs stabilize the cis-amide conformations through different mechanisms: (1) 5-alkylproline, with bulky hydrocarbon substituent on the C(delta) of proline, increases the cis-amide population through steric hindrance between the alkyl substituent and the N-terminal residues; (2) oxaproline or thioproline, the oxazolidine- or thiazolidine-derived proline analog, favors interactions between the dipole of the heterocyclic ring and the preceding carbonyl oxygen; and (3) azaproline, containing a nitrogen atom in place of the C(alpha) of proline, prefers the cis-amide bond by lone-pair repulsion between the alpha-nitrogen and the preceding carbonyl oxygen. Preference for the cis conformation was augmented by combining different modifications within a single proline. Azaproline and its derivatives are most effective in stabilizing cis-amide bonds without introducing additional steric bulk to compromise receptor interactions.

Effects of i and i+3 residue identity on cis-trans isomerism of the aromatic(i+1)-prolyl(i+2) amide bond: implications for type VI beta-turn formation

Cis-trans isomerization of amide bonds plays critical roles in protein molecular recognition, protein folding, protein misfolding, and disease. Aromatic-proline sequences are particularly prone to exhibit cis amide bonds. The roles of residues adjacent to a tyrosine-proline residue pair on cis-trans isomerism were examined. A short series of peptides XYPZ was synthesized and cis-trans isomerism was analyzed. Based on these initial studies, a series of peptides XYPN, X = all 20 canonical amino acids, was synthesized and analyzed by NMR for i residue effects on cis-trans isomerization. The following effects were observed: (a) aromatic residues immediately preceding Tyr-Pro disfavor cis amide bonds, with K(trans/cis)= 5.7-8.0, W > Y > F; (b) proline residues preceding Tyr-Pro lead to multiple species, exhibiting cis-trans isomerization of either or both X-Pro amide bonds; and (c) other residues exhibit similar values of K(trans/cis) (= 2.9-4.2), with Thr and protonated His exhibiting the highest fraction cis. beta-Branched and short polar residues were somewhat more favorable in stabilizing the cis conformation. Phosphorylation of serine at the i position modestly increases the stability of the cis conformer. In addition, the effect of the i+3 residue was examined in a limited series of peptides TYPZ. NMR data indicated that aromatic residues, Pro, Asn, Ala, and Val at the i+3 residue all favor cis amide bonds, with aromatic residues and Asn favoring more compact phi at Tyr(cis) and Ala and Pro favoring more extended phi at Tyr(cis). D-Alanine at the i+3 position particularly disfavors cis amide bonds.

Probing Nature’s Knots: The Folding Pathway of a Knotted Homodimeric Protein

The homodimeric protein YibK from Haemophilus influenzae belongs to a recently discovered superfamily of knotted proteins that has brought about a new protein-folding conundrum. Members of the alpha/beta-knot clan form deep trefoil knots in their native backbone structure, a topological feature that is currently unexplained in the protein-folding field. To help solve the puzzle of how a polypeptide chain can efficiently knot itself, the folding kinetics of YibK have been studied extensively and the results are reported here. Folding was monitored using probes for changes in both secondary and tertiary structure, and the monomer-dimer equilibrium was perturbed with a variety of solution conditions to allow characterisation of otherwise inaccessible states. Multiphasic kinetics were observed in the unfolding and refolding reactions of YibK, and under conditions where the dimer is favoured, dissociation and association were rate-limiting, respectively. A folding model consistent with all kinetic data is proposed: YibK appears to fold via two parallel pathways, partitioned by proline isomerisation events, to two distinct monomeric intermediates. These form a common third intermediate that is able to fold to native dimer. Kinetic simulations suggest that all intermediates are on-pathway. These results provide the valuable groundwork required to further understand how Nature codes for knot formation.

Diproline Templates as Folding Nuclei in Designed Peptides. Conformational Analysis of Synthetic Peptide Helices Containing Amino Terminal Pro-Pro Segments

The effect of N-terminal diproline segments in nucleating helical folding in designed peptides has been studied in two model sequences Piv-Pro-Pro-Aib-Leu-Aib-Phe-OMe (1) and Boc-Aib-Pro-Pro-Aib-Val-Ala-Phe-OMe (2). The structure of 1 in crystals, determined by X-ray diffraction, reveals a helical (alphaR) conformation for the segment residues 2 to 5, stabilized by one 4-->1 hydrogen bond and two 5-->1 interactions. The N-terminus residue, Pro(1) adopts a polyproline II (P(II)) conformation. NMR studies in three different solvent systems support a conformation similar to that observed in crystals. In the apolar solvent CDCl3, NOE data favor the population of both completely helical and partially unfolded structures. In the former, the Pro-Pro segment adopts an alphaR-alphaR conformation, whereas in the latter, a P(II)-alphaR structure is established. The conformational equilibrium shifts in favor of the P(II)-alphaR structure in solvents like methanol and DMSO. A significant population of the Pro(1)-Pro(2) cis conformer is also observed. The NMR results are consistent with the population of at least three conformational states about Pro-Pro segment: trans alphaR-alphaR, trans P(II)-alphaR and cis P(II)-alphaR. Of these, the two trans conformers are in rapid dynamic exchange on the NMR time scale, whereas the interconversion between cis and trans form is slow. Similar results are obtained with peptide 2. Analysis of 462 diproline segments in protein crystal structures reveals 25 examples of the alphaR-alphaR conformation followed by a helix. Modeling and energy minimization studies suggest that both P(II)-alphaR and alphaR-alphaR conformations have very similar energies in the model hexapeptide 1.

Conservation of cis prolyl bonds in proteins during evolution

In proteins and peptides, the vast majority of peptide bonds occurs in trans conformation, but a considerable fraction (about 5%) of X-Pro bonds adopts the cis conformation. Here we study the conservation of cis prolyl residues in evolutionary related proteins. We find that overall, in contrast to local, protein sequence similarity is a clear indicator for the conformation of prolyl residues. We observe that cis prolyl residues are more often conserved than trans prolyl residues, and both are more conserved than the surrounding amino acids, which show the same extent of conservation as the whole protein. The pattern of amino acid exchanges differs between cis and trans prolyl residues. Also, the cis prolyl bond is maintained in proteins with sequence identity as low as 20%. This finding emphasizes the importance of cis peptide bonds in protein structure and function.

Electronic control of amide cis-trans isomerism via the aromatic-prolyl interaction

The cis-trans isomerization of prolyl amide bonds results in large structural and functional changes in proteins and is a rate-determining step in protein folding. We describe a novel electronic strategy to control cis-trans isomerization, based on the demonstration that interactions between aromatic residues and proline are tunable by aromatic electronics. A series of peptides of sequence TXPN, X = Trp, pyridylalanine, pentafluorophenylalanine, or 4-Z-phenylalanine derivatives (Z = electron-donating, electron-withdrawing, or electron-neutral substituents), was synthesized and Ktrans/cis analyzed by NMR. Electron-rich aromatic residues stabilized cis amide bond formation, while electron-poor aromatics relatively favored trans amide bond formation. A Hammett correlation between aromatic electronics and cis-trans isomerization was observed. These results indicate that the interaction between aromatic residues and proline, which is observed to stabilize cis amide bonds and is also a general stabilizing interaction ubiquitous in proteins and protein-protein complexes, is not stabilized exclusively by a classical hydrophobic effect. To a large extent, the aromatic-prolyl interaction is driven and controllable by an electronic effect between the aromatic ring pi-electrons and the proline ring, consistent with a C-H-pi interaction as the key stabilizing force. The aromatic-prolyl interaction is electronically tunable by 0.9 kcal/mol and is enthalpic in nature. In addition, by combining aromatic ring electronics and stereoelectronic effects using 4-fluoroprolines, we demonstrate broad tuning (2.0 kcal/mol) of cis-trans isomerism in tetrapeptides. We demonstrate a simple tetrapeptide, TWflpN, that exhibits 60% cis amide bond and adopts a type VIa1 beta-turn conformation.

Laminated morphology of nontwisting beta-sheet fibrils constructed via peptide self-assembly

A synthetic peptide has been de novo designed that self-assembles into beta-sheet fibrils exhibiting a nontwisted, stacked morphology. The stacked morphology is constituted by 2.5 nm wide filaments that laterally associate to form flat fibril laminates exceeding 50 nm in width and micrometers in length. The height of each fibril is limited to the length of exactly one peptide monomer in an extended beta-strand conformation, approximately 7 nm. Once assembled, these highly ordered, 2-D structures are stable over a wide range of pH and temperature and exhibit characteristics similar to those of amyloid fibrils. Furthermore, the rate of assembly and degree of fibril lamination can be controlled with kinetic parameters of pH and temperature. Finally, the presence of a diproline peptide between two beta-sheet-forming strands in the peptide sequence is demonstrated to be an important factor in promoting the nontwisting, laminated fibril morphology.

Prediction of cis/trans isomerization in proteins using PSI-BLAST profiles and secondary structure information

Background

The majority of peptide bonds in proteins are found to occur in the trans conformation. However, for proline residues, a considerable fraction of Prolyl peptide bonds adopt the cis form. Proline cis/trans isomerization is known to play a critical role in protein folding, splicing, cell signaling and transmembrane active transport. Accurate prediction of proline cis/trans isomerization in proteins would have many important applications towards the understanding of protein structure and function.

Results

In this paper, we propose a new approach to predict the proline cis/trans isomerization in proteins using support vector machine (SVM). The preliminary results indicated that using Radial Basis Function (RBF) kernels could lead to better prediction performance than that of polynomial and linear kernel functions. We used single sequence information of different local window sizes, amino acid compositions of different local sequences, multiple sequence alignment obtained from PSI-BLAST and the secondary structure information predicted by PSIPRED. We explored these different sequence encoding schemes in order to investigate their effects on the prediction performance. The training and testing of this approach was performed on a newly enlarged dataset of 2424 non-homologous proteins determined by X-Ray diffraction method using 5-fold cross-validation. Selecting the window size 11 provided the best performance for determining the proline cis/trans isomerization based on the single amino acid sequence. It was found that using multiple sequence alignments in the form of PSI-BLAST profiles could significantly improve the prediction performance, the prediction accuracy increased from 62.8% with single sequence to 69.8% and Matthews Correlation Coefficient (MCC) improved from 0.26 with single local sequence to 0.40. Furthermore, if coupled with the predicted secondary structure information by PSIPRED, our method yielded a prediction accuracy of 71.5% and MCC of 0.43, 9% and 0.17 higher than the accuracy achieved based on the singe sequence information, respectively.

Conclusion

A new method has been developed to predict the proline cis/trans isomerization in proteins based on support vector machine, which used the single amino acid sequence with different local window sizes, the amino acid compositions of local sequence flanking centered proline residues, the position-specific scoring matrices (PSSMs) extracted by PSI-BLAST and the predicted secondary structures generated by PSIPRED. The successful application of SVM approach in this study reinforced that SVM is a powerful tool in predicting proline cis/trans isomerization in proteins and biological sequence analysis.

The folding energy landscape of apoflavodoxin is rugged: hydrogen exchange reveals nonproductive misfolded intermediates

Many native proteins occasionally form partially unfolded forms (PUFs), which can be detected by hydrogen/deuterium exchange and NMR spectroscopy. Knowledge about these metastable states is required to better understand the onset of folding-related diseases. So far, not much is known about where PUFs reside within the energy landscape for protein folding. Here, four PUFs of the relatively large apoflavodoxin (179 aa) are identified. Remarkably, at least three of them are partially misfolded conformations. The misfolding involves side-chain contacts as well as the protein backbone. The rates at which the PUFs interconvert with native protein have been determined. Comparison of these rates with stopped-flow data positions the PUFs in apoflavodoxin's complex folding energy landscape. PUF1 and PUF2 are unfolding excursions that start from native apoflavodoxin but do not continue to the unfolded state. PUF3 and PUF4 could be similar excursions, but their rates of formation suggest that they are on a dead-end folding route that starts from unfolded apoflavodoxin and does not continue all of the way to native protein. All PUFs detected thus are off the protein's productive folding route.