The effect of phosphorylation on the conformation of oligo-peptides with Ser–Pro motif: a molecular dynamics simulation

Model tetrapeptide system was designed to investigate the cis/trans isomerization of peptidyl-prolyl imide bond of Ser–Pro motif. To establish the side-chain O-phosphorylation effect in regulating the peptides conformations, molecular dynamics (MD) simulations where carried out on the designed tetrapeptides and their corresponding phosphorylated forms by MD Insight II Discovery3 approach. The most stable configurations and the statistic cis/trans concentration distribution demonstrated that the phosphorylation evidently influences the peptidyl-prolyl imide bond isomerization and works as a key effect in regulating the peptide conformations. The charge state and the site provided for the charge of the phosphate moiety might be an important key. The results also demonstrated that phosphorylation changes the cis conformation ratio of the peptide and the maximum cis value is obtained when the phosphate group has no negative charge.     

Prolyl isomerization: How significant for in vivo protein folding?

Suggestive but not decisive evidence indicates that in vivo peptide chain folding is completed in a time not much longer than that required for covalent peptide synthesis. Extrapolation of model peptide rates of the cis–trans prolyl isomerization leads to the prediction tht protein folding should be much slower than the apparent in vivo rates. On the assumption that rapid protein folding in vivo is the rule, three routes are suggested by which a protein undergoing biosynthesis can avoid a strongly slowed folding rate: (1) by a peptide chain-elongation process that adds only trans peptide bonds, follwed by a rapid folding process that incorporates them into a three-dimensional structure, raising the energy barrier to isomerization; (2) by folding to produce three dimensional structures that position prolyl residues largely in chain turns on the protein surface, where the residue may be either cis or trans without large effects on the protein structure and function; (3) prolyl cis–trans isomerization may be speeded by the formation of peptide loops.     

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 β 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 β 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 β 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 β 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.     

Kinetic coupling of folding and prolyl isomerization of beta2-microglobulin studied by mutational analysis

β₂-Microglobulin (β2-m), a protein responsible for dialysis-related amyloidosis, adopts a typical immunoglobulin domain fold with the N-terminal peptide bond of Pro32 in a cis isomer. The refolding of β2-m is limited by the slow trans-to-cis isomerization of Pro32, implying that intermediates with a non-native trans-Pro32 isomer are precursors for the formation of amyloid fibrils. To obtain further insight into the Pro-limited folding of β2-m, we studied the Gdn-HCl-dependent unfolding/refolding kinetics using two mutants (W39 and P32V β2-ms) as well as the wild-type β2-m. W39 β2-m is a triple mutant in which both of the authentic Trp residues (Trp60 and Trp95) are replaced by Phe and a buried Trp common to other immunoglobulin domains is introduced at the position of Leu39 (i.e., L39W/W60F/W95F). W39 β2-m exhibits a dramatic quenching of fluorescence upon folding, enabling a detailed analysis of Pro-limited unfolding/refolding. On the other hand, P32V β2-m is a mutant in which Pro32 is replaced by Val, useful for probing the kinetic role of the trans-to-cis isomerization of Pro32. A comparative analysis of the unfolding/refolding kinetics of these mutants including three types of double-jump experiments revealed the prolyl isomerization to be coupled with the conformational transitions, leading to apparently unusual kinetics, particularly for the unfolding. We suggest that careful consideration of the kinetic coupling of unfolding/refolding and prolyl isomerization, which has tended to be neglected in recent studies, is essential for clarifying the mechanism of protein folding and, moreover, its biological significance.     

Bioactivity of folding intermediates studied by the recovery of enzymatic activity during refolding

The peptide bond preceding proline residues realizes a cis/trans conformational switch with high switching resistance in native proteins and folding intermediates. Therefore, individual isomers have the potential to differ in bioactivity. However, information about isomer-specific bioactivities is difficult to obtain because of the risk of affecting isomeric distribution by bioactivity assay components.Here we present an approach that allows for the measurement of the recovery of enzymatic activities of wild-type RNase T₁ and RNase T₁ variants during refolding under conditions where the population of enzyme-substrate or enzyme-product complexes is negligible. Recovery of enzymatic activity was continuously monitored within the visible range of the spectrum by addition of a fluorescence-labeled nucleotide substrate to the refolding sample. We found that a nonnative trans conformation at Pro39 renders the RNase T₁ almost completely inactive. A folding intermediate having a nonnative trans conformation at Pro55 shows about 46% of the enzymatic activity referred to the native state. Pro55, in contrast to the active site located Pro39, is situated in a solvent-exposed loop region remote from active-site residues. In both cases, peptidyl prolyl cis/trans isomerases accelerate the regain of nucleolytic activity. Our findings show that even if there is a considerable distance between the site of isomerization and the active site, conformational control of the bioactivity of proteins is likely to occur, and that the surface location of prolyl bonds suffices for the control of buried active sites mediated by peptidyl prolyl cis/trans isomerases.     

Structure of an early native‐like intermediate of β2‐microglobulin amyloidogenesis

To investigate early intermediates of β2‐microglobulin (β2m) amyloidogenesis, we solved the structure of β2m containing the amyloidogenic Pro32Gly mutation by X‐ray crystallography. One nanobody (Nb24) that efficiently blocks fibril elongation was used as a chaperone to co‐crystallize the Pro32Gly β2m monomer under physiological conditions. The complex of P32G β2m with Nb24 reveals a trans peptide bond at position 32 of this amyloidogenic variant, whereas Pro32 adopts the cis conformation in the wild‐type monomer, indicating that the cis to trans isomerization at Pro32 plays a critical role in the early onset of β2m amyloid formation.     

Steric effects of cis-trans isomerism on neighboring residues in proline oligopeptides: A 13C-nmr study of conformational heterogeneity in linear tripeptides

A series of proline-containing linear oligopeptides (4 dipeptides and 15 tripeptides) were synthesized and examined in aqueous and nonaqueous solutions using ¹³C-nmr spectroscopy. Spectra of linear tripeptides showing cis-trans isomerism about the X-Pro bond (X = Pro, Gly, and Ala) also show neighboring effects on the chemical shifts of residues both preceding and following the prolyl moiety. The extent of cis-trans isomerism observed about the X-Pro peptide bond correlates not only with the nature of X, but also depends on the size of the residue following proline; the larger substituents favor an increase in cis content about the X-Pro bond.     

Phosphorylation and Prolyl Isomerization Independently Regulate the Signal Adapter Function of CrkII

The signaling protein CrkII switches between forms with high or low binding affinity. Both phosphorylation and native-state prolyl isomerization were suggested to regulate the transition between these forms. Here we analyzed how phosphorylation at Tyr222 and Tyr252 and the Pro238Ala substitution affect signal transfer of human and chicken CrkII to a downstream target. Human CrkII is regulated by phosphorylation only, but chicken CrkII is regulated by Pro238 trans → cis isomerization and by Tyr222 phosphorylation. Surprisingly, they act in an independent fashion. Apparently, the allosteric transition to a low-activity form can be induced by phosphorylation or prolyl isomerization located at distant sites in CrkII.     

A nonessential role for Arg 55 in cyclophilin18 for catalysis of proline isomerization during protein folding

The protein folding process is often in vitro rate‐limited by slow cis‐trans proline isomerization steps. Importantly, the rate of this process in vivo is accelerated by prolyl isomerases (PPIases). The archetypal PPIase is the human cyclophilin 18 (Cyp18 or CypA), and Arg 55 has been demonstrated to play a crucial role when studying short peptide substrates in the catalytic action of Cyp18 by stabilizing the transition state of isomerization. However, in this study we show that a R55A mutant of Cyp18 is as efficient as the wild type to accelerate the refolding reaction of human carbonic anhydrase II (HCA II). Thus, it is evident that the active‐site located Arg 55 is not required for catalysis of the rate‐limiting prolyl cis‐trans isomerization steps during the folding of a protein substrate as HCA II. Nevertheless, catalysis of cis‐trans proline isomerization in HCA II occurs in the active‐site of Cyp18, since binding of the inhibitor cyclosporin A abolishes rate acceleration of the refolding reaction. Obviously, the catalytic mechanisms of Cyp18 can differ when acting upon a simple model peptide, four residues long, with easily accessible Pro residues compared with a large protein molecule undergoing folding with partly or completely buried Pro residues. In the latter case, the isomerization kinetics are significantly slower and simpler mechanistic factors such as desolvation and/or strain might operate during folding‐assisted catalysis, since binding to the hydrophobic active site is still a prerequisite for catalysis.     

Mechanistic insights into Pin1 peptidyl‐prolyl cis‐trans isomerization from umbrella sampling simulations

The peptidyl‐proyl isomerase Pin1 plays a key role in the regulation of phospho(p)‐Ser/Thr‐Pro proteins, acting as a molecular timer of the cell cycle. After recognition of these motifs, Pin1 catalyzes the rapid cis‐trans isomerization of proline amide bonds of substrates, contributing to maintain the equilibrium between the two conformations. Although a great interest has arisen on this enzyme, its catalytic mechanism has long been debated. Here, the cis‐trans isomerization of a model peptide system was investigated by means of umbrella sampling simulations in the Pin1‐bound and unbound states. We obtained free energy barriers consistent with experimental data, and identified several enzymatic features directly linked to the acceleration of the prolyl bond isomerization. In particular, an enhanced autocatalysis, the stabilization of perturbed ground state conformations, and the substrate binding in a procatalytic conformation were found as main contributions to explain the lowering of the isomerization free energy barrier. Proteins 2014; 82:2943–2956. © 2014 Wiley Periodicals, Inc.     

Solid-state geometry and conformation of linear,diastereoisomeric oligoprolines

We have carried out a systematic analysis of the solid-state conformational preferences of a number of linear homo-oligoprolines (to the tetramer) by ir absorption and x-ray diffraction. The peptides present different chiral sequences (tacticities), various types (urethane and amide) of N-protecting groups, and free and blocked C-termini (which imply different capabilities of forming H-bonds). The following conclusions can be drawn: (i) values for the geometry of the prolyl residue and the peptide bond in the cis and in the trans conformations are proposed; (ii) in general the conformational angles φ and ψ in the linear homo-oligoprolines have values appropriate for the polyproline II structure (conformation F); (iii) the pyrrolidine ring shows various types of puckering with no apparent relation to the backbone conformation; (iv) Pro-Pro peptide bonds generally take the trans conformation, the few cases of cis conformation being formed by Pro residues of different chirality; (v) the single H-bond donor — OH, when present, is always bonded to H-acceptors, which can be either the urethane or the amide or the peptide carbonyl but never the carbonyl group of the — COOH moiety.     

The slow step of folding of Staphylococcus aureus PC1 β-lactamase involves the collapse of a surface loop rate limited by the Trans to Cis isomerization of a non-proline peptide bond

We wished to test the hypothesis that the non proline cis to trans isomerization of the peptide bond at position 167 in the S. aureus β-lactamase PC1 exerts a significant controlling effect on the folding pathway of this enzyme. The previous data presented in support of this hypothesis could not rule out the effect of factors unrelated to non-proline cis/trans isomerization. We have used the plasmid pET9d to direct soluble overproduction of the S. aureus β-lactamase PC1 and a site-directed mutant (Ile 167 to Pro) in Escherichia coli. Following purification the proteins were subjected to a comparative analysis of the kinetics of unfolding and refolding using the techniques of near- and far-UV circular dichroism spectroscopy and fluorescence spectroscopy in conjunction with “double-jump” experiments. Results show that the fully-unfolded I167P mutant enzyme retains 20% of molecules in a fast-refolding form and that slower-refolding molecules fold faster than the recombinant wild-type enzyme. The final stage of folding involves folding of the Ω-loop into a conformation essential for enzymatic activity. In support of the original hypothesis, the folding of this Ω-loop is rate limited by the isomerization of the Glu 166-Ile 167 peptide bond. Proteins 33:550–557, 1998. © 1998 Wiley-Liss, Inc.     

Cis-trans isomerization of the Epstein-Barr virus determinant peptide EENLLDFVRF after the DM1 TCR recognition of the HLA-B*4405/peptide complex

The Epstein–Barr virus determinant peptide EENLLDFVRF shows high immunogenicity when presented by HLA-B*4405 allotype. This fact is accompanied by a cis–trans isomerization of the Leu5-Asp6 peptide bond upon TCR binding of the pMHC complex. Molecular dynamics simulations of pMHC/TCR structures, with the EENLLDFVRF peptide in cis and trans conformations have been employed in order to examine the structure and dynamics of the pMHC complex with such an unusual conformation. The results, based on MM-PBSA free energy computations as well as buried surface area analysis and interactions at the pMHC/TCR interface, indicate that the TCR binds preferably the pMHC complex with the Leu5-Asp6 peptide bond in cis conformation. It is the first time that this notable conformational feature of T-cell epitope is investigated.     

Solid Phase Synthesis of Phosphopeptides Incorporating 2,2-Dimethyloxazolidine Pseudoproline Analogs: Evidence for <Emphasis Type="Italic">trans</Emphasis> Leu-Pro Peptide Bonds in Stat3 Inhibitors

To answer the question of whether the conformation of the Leu-Pro bond is cis or trans in Ac-pTyr-Leu-Pro-Gln-Thr-Val-NH₂ when complexed with the SH2 domain of Stat3, we substituted 2,2-dimethyloxazolidines derived from serine (Ser(Ψ^Me,Mepro)) and threonine (Thr(Ψ^Me,Mepro)) for proline. The 2,2-dimethyloxazolidine and 2,2-dimethylthiazolidine pseudoproline (ΨPro) analogs induce predominantly cis Xxx-ΨPro peptide bonds. As these ΨPro analogs are acid-labile, the phosphopeptides were synthesized using Fmoc-based SPPS using unprotected phosphotyrosine and 4-hydroxybenzoate as the linker that allowed release from the support by alkaline ammonolysis, conditions that kept the oxazolidine rings intact. Incorporation of Ser(Ψ^Me,Mepro) resulted in 69% cis Leu-ΨPro bond content in aqueous solution whereas that for Thr(Ψ^Me,Mepro) analog was 63%. Affinities for Stat3 were 3–5 fold lower than the lead compound and were inversely correlated with cis content. Thus we conclude that the Leu-Pro peptide bond is trans when the peptide is bound to Stat3.     

The rate-limiting step in the folding of the cis-Pro167Thr mutant of TEM-1 β-lactamase is the trans to cis isomerization of a non-proline peptide bond

The stability and kinetics of unfolding and refolding of the P167T mutant of the TEM-1 β-lactamase have been investigated as a function of guanidine hydrochloride concentration. The activity of the mutant enzyme was not significantly modified, which strongly suggests that the Glu166–Thr167 peptide bond, like the Glu166–Pro167, is cis. The mutation, however, led to a significant decrease in the stability of the native state relative to both the thermodynamically stable intermediate and the fully unfolded state of the protein. In contrast to the two slower phases seen in the refolding of the wild-type enzyme, only one phase was detected in the refolding of the mutant, indicating a determining role of proline 167 in the kinetics of folding of the wild-type enzyme. The former phases are replaced by rapid refolding when the enzyme is unfolded for short periods of time, but the latter is independent of the time of unfolding. The monophasic refolding reaction of the mutant is proposed to reflect mainly the trans→cis isomerization of the Glu166–Thr167 peptide bond. © 1996 John Wiley & Sons, Inc.     

Cation binding cyclic peptides composed of imino acid residues

Cyclo(L -Pro-Sar)_n (n = 2–4) with moderate flexibility and hydrophobicity of molecular structure was synthesized, and the characteristics of these cyclic peptides and their metal complexes in acetonitrile were investigated in connection with the residual properties using ¹³C-nmr measurements. The cyclic tetrapeptide cyclo(L -Pro-Sar)₂ showed a sterically hindered phenomenon in acetonitrile in which the amide backbone adopted a cis-trans-cis-trans sequence. The cyclic hexapeptide cyclo(L -Pro-Sar)₃ existed as a mixture of several conformers whose interconversion is slow on the nmr time scale, including cis-cis-trans and/or cis-trans-trans arrangement of the Sar-Pro bond. Finally, it was demonstrated that the cyclic octapeptide cyclo(L -Pro-Sar)₄ behaved as a mixture of multiple conformers which allowed for cis-trans isomerism about the Pro-Sar peptide bond, of which 20–30% had the all-cis Sar-Pro bond isomer and the remaining 70–80% had one (or more) cis Sar-Pro bond isomer. ¹³C-nmr spectra also demonstrated that cyclo(L -Pro-Sar)_n (n = 3,4) formed a 1:1 ion complex whose conformation was characterized by an all-trans peptide bond in the presence of excess metal salt. Cation binding studies, using CD measurements, established that the ion selectivity of cyclo(L -Pro-Sar)₄ in acetonitrile decreased in the order, Ba²⁺ > Ca²⁺ > Na⁺ > Mg²⁺ > Li⁺.     

Prediction ofcis/transisomerization in proteins using PSI-BLAST profiles and secondary structure information

Background  

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

Loss of intramolecular electrostatic interactions and limited conformational ensemble may promote self‐association of cis–tau peptide

Self‐association of proteins can be triggered by a change in the distribution of the conformational ensemble. Posttranslational modification, such as phosphorylation, can induce a shift in the ensemble of conformations. In the brain of Alzheimer's disease patients, the formation of intra‐cellular neurofibrillary tangles deposition is a result of self‐aggregation of hyper‐phosphorylated tau protein. Biochemical and NMR studies suggest that the cis peptidyl prolyl conformation of a phosphorylated threonine‐proline motif in the tau protein renders tau more prone to aggregation than the trans isomer. However, little is known about the role of peptidyl prolyl cis/trans isomerization in tau aggregation. Here, we show that intra‐molecular electrostatic interactions are better formed in the trans isomer. We explore the conformational landscape of the tau segment containing the phosphorylated‐Thr²³¹‐Pro²³² motif using accelerated molecular dynamics and show that intra‐molecular electrostatic interactions are coupled to the isomeric state of the peptidyl prolyl bond. Our results suggest that the loss of intra‐molecular interactions and the more restricted conformational ensemble of the cis isomer could favor self‐aggregation. The results are consistent with experiments, providing valuable complementary atomistic insights and a hypothetical model for isomer specific aggregation of the tau protein. Proteins 2015; 83:436–444. © 2014 Wiley Periodicals, Inc.     

Biotransformation of Unsaturated Long-Chain Fatty Acids by Eubacterium lentum

Eubacterium lentum (33 strains) isomerized the 12-cis double bond of C₁₈ fatty acids with cis double bonds at C-9 and C-12 into an 11-trans double bond before reduction of the 9-cis double bond. The 14-cis double bond of homo-γ-linolenic acid was isomerized by 29 strains into a 13-trans double bond. The same strains isomerized the 14-cis double bond of arachidonic acid into a 13-trans double bond and then isomerized the 8-cis double bond into a 7-trans double bond; the 13-cis double bond of 10-cis, 13-cis-nonadecadienoic acid was isomerized into a 12-trans double bond. None of these isomerization products was further reduced. Studies with resting cells showed optimal isomerization velocity at a linoleic acid concentration of 37.5 μM; higher concentrations were inhibitory. The pH optimum for isomerization was 7.5 to 8.5. The isomerase was inhibited by the sulfhydryl reagents iodoacetamide, bromoacetate, and N-ethylmaleimide and by the chelators EDTA and 1,10-phenanthroline.     

Carbon isotope fractionation during <Emphasis Type="Italic">cis</Emphasis>–<Emphasis Type="Italic">trans</Emphasis> isomerization of unsaturated fatty acids in <Emphasis Type="Italic">Pseudomonas putida</Emphasis>

The molecular mechanism of the unique cis to trans isomerization of unsaturated fatty acids in the solvent-tolerant bacterium Pseudomonas putida S12 was studied. For this purpose, the carbon isotope fractionation of the cis–trans isomerase was estimated. In resting cell experiments, addition of 3-nitrotoluene for activation of the cis–trans isomerase resulted in the conversion of the cis-unsaturated fatty acids into the corresponding trans isomers. For the conversion of C16:1 cis to its corresponding trans isomer, a significant fractionation was measured. The intensity of this fractionation strongly depended on the rate of cis–trans isomerization and the added concentration of 3-nitrotoluene, respectively. The presence of a significant fractionation provides additional indication for a transition from the sp² carbon linkage of the cis-double bond to an intermediate sp³ within an enzyme–substrate complex. The sp² linkage is reconstituted after rotation to the trans configuration has occurred. As cytochrome c plays a major role in the catabolism of Cti polypeptide, these findings favour a mechanism for the enzyme in which electrophilic iron (Fe³⁺), provided by a heme domain, removes an electron of the cis double bond thereby transferring the sp² linkage into sp³.