Cis–trans isomerization of omega dihedrals in proteins

Peptide bonds in protein structures are mainly found in trans conformation with a torsion angle ω close to 180°. Only a very low proportion is observed in cis conformation with ω angle around 0°. Cis–trans isomerization leads to local conformation changes which play an important role in many biological processes. In this paper, we reviewed the recent discoveries and research achievements in this field. First, we presented some interesting cases of biological processes in which cis–trans isomerization is directly implicated. It is involved in protein folding and various aspect of protein function like dimerization interfaces, autoinhibition control, channel gating, membrane binding. Then we reviewed conservation studies of cis peptide bonds which emphasized evolution constraints in term of sequence and local conformation. Finally we made an overview of the numerous molecular dynamics studies and prediction methodologies already developed to take into account this structural feature in the research area of protein modeling. Many cis peptide bonds have not been recognized as such due to the limited resolution of the data and to the refinement protocol used. Cis–trans proline isomerization reactions represents a vast and promising research area that still needs to be further explored for a better understanding of isomerization mechanism and improvement of cis peptide bond predictions.     

Ligand Binding Turns Moth Pheromone-binding Protein into a pH Sensor: EFFECT ON THE ANTHERAEA POLYPHEMUS PBP1 CONFORMATION*

In moths, pheromone-binding proteins (PBPs) are responsible for the transport of the hydrophobic pheromones to the membrane-bound receptors across the aqueous sensillar lymph. We report here that recombinant Antheraea polyphemus PBP1 (ApolPBP1) picks up hydrophobic molecule(s) endogenous to the Escherichia coli expression host that keeps the protein in the “open” (bound) conformation at high pH but switches to the “closed” (free) conformation at low pH. This finding has bearing on the solution structures of undelipidated lepidopteran moth PBPs determined thus far. Picking up a hydrophobic molecule from the host expression system could be a common feature for lipid-binding proteins. Thus, delipidation is critical for bacterially expressed lipid-binding proteins. We have shown for the first time that the delipidated ApolPBP1 exists primarily in the closed form at all pH levels. Thus, current views on the pH-induced conformational switch of PBPs hold true only for the ligand-bound open conformation of the protein. Binding of various ligands to delipidated ApolPBP1 studied by solution NMR revealed that the protein in the closed conformation switches to the open conformation only at or above pH 6.0 with a protein to ligand stoichiometry of ∼1:1. Mutation of His⁷⁰ and His⁹⁵ to alanine drives the equilibrium toward the open conformation even at low pH for the ligand-bound protein by eliminating the histidine-dependent pH-induced conformational switch. Thus, the delipidated double mutant can bind ligand even at low pH in contrast to the wild type protein as revealed by fluorescence competitive displacement assay using 1-aminoanthracene and solution NMR.     

Oligomeric States along the Folding Pathways of β2-Microglobulin: Kinetics,Thermodynamics, and Structure

The transition of proteins from their soluble functional state to amyloid fibrils and aggregates is associated with the onset of several human diseases. Protein aggregation often requires some structural reshaping and the subsequent formation of intermolecular contacts. Therefore, the study of the conformation of excited protein states and their ability to form oligomers is of primary importance for understanding the molecular basis of amyloid fibril formation. Here, we investigated the oligomerization processes that occur along the folding of the amyloidogenic human protein β2-microglobulin. The combination of real-time two-dimensional NMR data with real-time small-angle X-ray scattering measurements allowed us to derive thermodynamic and kinetic information on protein oligomerization of different conformational states populated along the folding pathways. In particular, we could demonstrate that a long-lived folding intermediate (I-state) has a higher propensity to oligomerize compared to the native state. Our data agree well with a simple five-state kinetic model that involves only monomeric and dimeric species. The dimers have an elongated shape with the dimerization interface located at the apical side of β2-microglobulin close to Pro32, the residue that has a trans conformation in the I-state and a cis conformation in the native (N) state. Our experimental data suggest that partial unfolding in the apical half of the protein close to Pro32 leads to an excited state conformation with enhanced propensity for oligomerization. This excited state becomes more populated in the transient I-state due to the destabilization of the native conformation by the trans-Pro32 configuration.     

Conformational changes induced by the transforming amino acid substitution in the transmembrane domain of theneu oncogene-encoded p185 protein

Theneu oncogene is frequently found in certain types of human carcinomas and has been shown to be activated in animal models by nitrosourea-induced mutation. The activating mutation in theneu oncogene results in the substitution of a glutamic acid for a valine at position 664 in the transmembrane domain of the encoded protein product of 185 kda (designated p185), which, on the basis of homology studies, is presumed to be a receptor for an as yet unidentified growth factor. It has been proposed that activating amino acid substitutions in this region of p185 lead to a conformational change in the protein which causes signal transduction via an increase in tyrosine kinase activity in the absence of any external signal. Using conformational energy analysis, we have determined the preferred three-dimensional structures for the transmembrane decapeptide (residues 658–667) of the p185 protein with valine and glutamic acid at the critical position 664. The results indicate that the global minimum energy conformation of the decapeptide from the normal protein with Val at position 664 is an α-helix with a sharp bend (CD* conformation at residues 664 and 665) in this region, whereas the global minimum conformation for the decapeptide from the mutant transforming protein with Glu at position 664 assumes an all α-helical configuration. Furthermore, the second highest energy conformation for the decapeptide from the normal protein is identical to the global minimum energy conformation for the decapeptide from the transforming protein, providing a possible explanation why overexpression of the normal protein also has a transforming effect. These results suggest there may be a normal and a transforming conformation for theneu-encoded p185 proteins which may explain their differences in transforming activity.     

FRET-based detection of different conformations of MK2

MAP kinase-activated protein kinase 2 (MK2 or MAPKAP K2) is a stress-activated enzyme downstream to p38 MAPK. By fusion of green fluorescent protein variants to the N- and C-terminus we analysed conformational changes in the kinase molecule in vitro and in vivo. Activation of MK2 is accompanied by a decrease in fluorescence resonance energy transfer, indicating a transition from an inactive/closed to an active/open conformation with an increase in the apparent distance between the fluorophores of ~9 Å. The closed conformation exists exclusively in the nucleus. Upon stress, the open conformation of MK2 rapidly becomes detectable in the cytoplasm and accumulates in the nucleus only when Crm1-dependent nuclear export is blocked. Hence, in living cells activation of MK2 and its nuclear export are coupled by a phosphorylation-dependent conformational switch.     

Two-dimensional infrared correlation spectroscopy study of the interaction of oxidized and reduced cytochrome c with phospholipid model membranes

We have used two-dimensional infrared correlation spectroscopy (2D-IR) to study the interaction and conformation of cytochrome c in the presence of a binary phospholipid mixture composed of a zwitterionic perdeuterated phospholipid and a negatively-charged one. The influence of the main temperature phase transition of the phospholipid model membranes on the conformation of cytochrome c has been evaluated by monitoring both the Amide I′ band of the protein and the CH₂ and CD₂ stretching bands of the phospholipids. Synchronous 2D-IR analysis has been used to determine the different secondary structure components of cytochrome c which are involved in the specific interaction with the phospholipids, revealing the existence of a specific interaction between the protein with cardiolipin-containing vesicles but not with phosphatidic acid-containing ones. Interestingly, 2D-IR is capable of showing the existence of significant changes in the protein conformation at the same time that the phospholipid transition occurs. In summary, 2D-IR revealed an important effect of the phospholipid phase transition of cardiolipin on the secondary structure of oxidized cytochrome c but not to either reduced cytochrome c or in the presence of phosphatidic acid, demonstrating the existence of specific intermolecular interactions between cardiolipin and cytochrome c.     

Fluid Shear Induces Conformation Change in Human Blood Protein von Willebrand Factor in Solution

Many of the physiological functions of von Willebrand Factor (VWF), including its binding interaction with blood platelets, are regulated by the magnitude of applied fluid/hydrodynamic stress. We applied two complementary strategies to study the effect of fluid forces on the solution structure of VWF. First, small-angle neutron scattering was used to measure protein conformation changes in response to laminar shear rates (G) up to 3000/s. Here, purified VWF was sheared in a quartz Couette cell and protein conformation was measured in real time over length scales from 2-140 nm. Second, changes in VWF structure up to 9600/s were quantified by measuring the binding of a fluorescent probe 1,1′-bis(anilino)-4-,4′-bis(naphtalene)-8,8′-disulfonate (bis-ANS) to hydrophobic pockets exposed in the sheared protein. Small angle neutron scattering studies, coupled with quantitative modeling, showed that VWF undergoes structural changes at G < 3000/s. These changes were most prominent at length scales <10 nm (scattering vector (q) range >0.6/nm). A mathematical model attributes these changes to the rearrangement of domain level features within the globular section of the protein. Studies with bis-ANS demonstrated marked increase in bis-ANS binding at G > 2300/s. Together, the data suggest that local rearrangements at the domain level may precede changes at larger-length scales that accompany exposure of protein hydrophobic pockets. Changes in VWF conformation reported here likely regulate protein function in response to fluid shear.     

Effect of proline residues on protein folding

Conformational energy calculations have been used to study the role of the proline residues in the folding of bovine pancreatic trypsin inhibitor. In the calculation, each of the four proline residues of this small protein is forced from the trans to cis peptide isomer while still part of the native folded structure. The cis proline residue can always be accommodated by small changes of the native conformation (< 1 Å root-mean-square deviation). For three of the four proline residues, Pro2, Pro9 and Pro 13, being in the cis form is calculated to destabilize the folded conformation by less than 11 kcal/mol, suggesting that rapid folding to a stable native-like conformation can occur with either isomeric form. For one of these three, Pro13, the destabilization is only 1 kcal/mol, suggesting the existence of an alternative folded native conformation with Pro13 cis. The fourth proline residue, Pro8, is calculated to destabilize the native conformation by so much (33 kcal/mol) that it will block folding in the manner proposed by Brandts et al. (1975).     

The GPI anchor signal sequence dictates the folding and functionality of the Als5 adhesin from Candida albicans

Background

Proteins destined to be Glycosylphosphatidylinositol (GPI) anchored are translocated into the ER lumen completely before the C-terminal GPI anchor attachment signal sequence (SS) is removed by the GPI-transamidase and replaced by a pre-formed GPI anchor precursor. Does the SS have a role in dictating the conformation and function of the protein as well?

Methodology/Principal Findings

We generated two variants of the Als5 protein without and with the SS in order to address the above question. Using a combination of biochemical and biophysical techniques, we show that in the case of Als5, an adhesin of C. albicans, the C-terminal deletion of 20 amino acids (SS) results in a significant alteration in conformation and function of the mature protein.

Conclusions/Significance

We propose that the locking of the conformation of the precursor protein in an alternate conformation from that of the mature protein is one probable strategy employed by the cell to control the behaviour and function of proteins intended to be GPI anchored during their transit through the ER.     

Solution of Levinthal’s paradox is possible at the level of the formation and assembly of protein secondary structures

The complete volume of a protein’s conformation space is smaller by many orders of magnitude at the level of secondary-structure elements as compared with the conformation of amino-acid residues. According to Levinthal’s estimate, the latter is ~10^2L, with L being the number of residues in the chain, while the former, at the level of secondary structures, increases no faster than ~L^N with N being the number of the secondary-structure elements. N is approximately L/15 according to the statistics of protein structures. This drastic decrease in the exponent (L/15 instead of 2L) considerably reduces the sampling space and explains the reason that sampling of the conformation space at the level of secondary-structure elements does not prevent the protein chain from finding its most stable structure.     

In vitro protein folding by E. coli ribosome: unfolded protein splitting 70S to interact with 50S subunit

Folding of unfolded protein on Escherichia coli 70S ribosome is accompanied by rapid dissociation of the ribosome into 50S and 30S subunits. The dissociation rate of 70S ribosome with unfolded protein is much faster than that caused by combined effect of translation and polypeptide release factors known to be involved in the dissociation of ribosome into subunits. The protein then reaches a “folding competent” state on 50S and is released to take up native conformation by itself. Release before attaining the folding competent state or prevention of release by cross-linking it with ribosome, would not allow the protein to get back to its native conformation.     

Fusion of barnase to antiferritin antibody F11 VH domain results in a partially folded functionally active protein

A chimeric protein, VH-barnase, was obtained by fusing the VH domain of anti-human ferritin monoclonal antibody F11 to barnase, a bacterial RNase from Bacillus amyloliquefaciens. After refolding from inclusion bodies, the fusion protein formed insoluble aggregates. Off-pathway aggregation was significantly reduced by adding either purified GroEL/GroES chaperones or arginine, with 10–12-fold increase in the yield of the soluble protein. The final protein conformation was identical by calorimetric criteria and CD and fluorescence spectroscopy to that obtained without additives, thus suggesting that VH-barnase structure does not depend on folding conditions. Folding of VH-barnase resulted in a single calorimetrically revealed folding unit, the so-called “calorimetric domain”, with conformation consistent with a molten globule that possessed well-defined secondary structure and compact tertiary conformation with partial exposure of hydrophobic patches and low thermodynamic stability. The unique feature of VH-barnase is that, despite the partially unfolded conformation and coupling into a single “calorimetric domain”, this immunofusion retained both the antigen-binding and RNase activities that belong to the two heterologous domains.     

Modeling of the three-dimensional structure of luffin-α and its simulated reaction with the substrate oligoribonucleotide GAGA

A fundamental problem in biochemistry and molecular biology is understanding the spatial structure of macromolecules and then analyzing their functions. In this study, the three-dimensional structure of a ribosome-inactivating protein luffin-α was predicted using a neural network method and molecular dynamics simulation. A feedforward neural network with the backpropagation learning algorithm were trained on model class of homologous proteins including trichosanthin andα-momorcharin. The distance constraints for the Cα atoms in the protein backbone were utilized to generate a folded crude conformation of luffin-α by model building and the steepest descent minimization approach. The crude conformation was refined by molecular dynamics techniques and a simulated annealing procedure. The interaction between luffin-α and its analogous substrate GAGA was also simulated to understand its action mechanism.     

Single-Stranded DNA within Nanopores: Conformational Dynamics and Implications for Sequencing; a Molecular Dynamics Simulation Study

Engineered protein nanopores, such as those based on α-hemolysin from Staphylococcus aureus have shown great promise as components of next-generation DNA sequencing devices. However, before such protein nanopores can be used to their full potential, the conformational dynamics and translocation pathway of the DNA within them must be characterized at the individual molecule level. Here, we employ atomistic molecular dynamics simulations of single-stranded DNA movement through a model α-hemolysin pore under an applied electric field. The simulations enable characterization of the conformations adopted by single-stranded DNA, and allow exploration of how the conformations may impact on translocation within the wild-type model pore and a number of mutants. Our results show that specific interactions between the protein nanopore and the DNA can have a significant impact on the DNA conformation often leading to localized coiling, which in turn, can alter the order in which the DNA bases exit the nanopore. Thus, our simulations show that strategies to control the conformation of DNA within a protein nanopore would be a distinct advantage for the purposes of DNA sequencing.     

Role of Protein Stabilizers on the Conformation of the Unfolded State of Cytochrome c and Its Early Folding Kinetics: INVESTIGATION AT SINGLE MOLECULAR RESOLUTION*

An insight into the conformation and dynamics of unfolded and early intermediate states of a protein is essential to understand the mechanism of its aggregation and to design potent inhibitor molecules. Fluorescence correlation spectroscopy has been used to study the effects of several model protein stabilizers on the conformation of the unfolded state and early folding dynamics of tetramethyl rhodamine-labeled cytochrome c from Saccharomyces cerevisiae at single molecular resolution. Special attention has been given to arginine, which is a widely used stabilizer for improving refolding yield of different proteins. The value of the hydrodynamic radius (r_H) obtained by analyzing the intensity fluctuations of the diffusing molecules has been found to increase in a two-state manner as the protein is unfolded by urea. The results further show that the presence of arginine and other protein stabilizers favors a relatively structured conformation of the unfolded states (r_H of 29 Å) over an extended one (r_H of 40 Å), which forms in their absence. Also, the time constant of a kinetic component (τ_R) of about 30 μs has been observed by analyzing the correlation functions, which represents formation of a collapsed state. This time constant varies with urea concentration representing an inverted Chevron plot that shows a roll-over and behavior in the absence of arginine. To the best of our knowledge, this is one of the first applications of fluorescence correlation spectroscopy to study direct folding kinetics of a protein.     

Protein L20 from the large subunit of Escherichia coli ribosomes is an assembly protein

When 50 S subunits from Escherichia coli are incubated in the presence of 4.3 m-LiCl, the resulting 4.3c core particle quantitatively lacks L20 in addition to other proteins. The 4.3c core can be reconstituted to an active 50 S subunit in the presence of total 50 S proteins by means of the second step incubation of the two-step reconstitution procedure. This finding indicates that the conformation of the 4.3c core is at least equivalent to the conformation of the reconstitution intermediate RI₅₀¹(1) particle, which is exclusively formed in the first-step incubation. It follows that L20 is not necessary for the maintenance of the 4.3c core conformation. In contrast, the total reconstitution of an active 50 S particle from (23 S + 5 S) RNA and a protein preparation lacking L20 was fully dependent on the addition of L20. However, when the 4.3c core, which does not contain L20, is reconstituted with the same protein fraction, the activity of the resulting particle did not depend on the presence of L20. Thus, L20 is essential for the early assembly (occurring in the first-step incubation) but plays no role either in the late assembly steps, or the functions of the mature 50 S particle.Heat treatment of the 4.3c core distorts the 4.3c core conformation and leads to particles with lower s values. The degradation of the 4.3c core conformation is reduced when L20 is added. A further stabilization is obtained by the addition of (L20 + L24). Thus, L20 is dispensable for the maintenance of the 4.3c core conformation, but stabilizes this conformation.     

Reinforced Polyproline II Conformation in a Hydroxyproline-Rich Cell Wall Glycoprotein from Carrot Root

The salt-extractable hydroxyproline-rich cell wall glycoprotein from carrot (Daucus carota L.) roots is composed of 35% (w/w) protein, 3% (w/w) galactose, and 62% (w/w) arabinose. The arabinose is attached to hydroxyproline as tetra- and trisaccharides. The circular dichroism of the glycoprotein shows that it is completely in the polyproline II conformation. After deglycosylation of the glycoprotein, the polyproline II conformation of the peptide backbone was lost. This indicates that the carbohydrate reinforces the polyproline II conformation.     

Structure of H/ACA RNP protein Nhp2p reveals cis/trans isomerization of a conserved proline at the RNA and Nop10 binding interface

H/ACA small nucleolar and Cajal body ribonucleoproteins (RNPs) function in site-specific pseudouridylation of eukaryotic rRNA and snRNA, rRNA processing, and vertebrate telomerase biogenesis. Nhp2, one of four essential protein components of eukaryotic H/ACA RNPs, forms a core trimer with the pseudouridylase Cbf5 and Nop10 that binds to H/ACA RNAs specifically. Crystal structures of archaeal H/ACA RNPs have revealed how the protein components interact with each other and with the H/ACA RNA. However, in place of Nhp2p, archaeal H/ACA RNPs contain L7Ae, which binds specifically to an RNA K-loop motif absent from eukaryotic H/ACA RNPs, while Nhp2 binds a broader range of RNA structures. We report solution NMR studies of Saccharomyces cerevisiae Nhp2 (Nhp2p), which reveal that Nhp2p exhibits two major conformations in solution due to cis/trans isomerization of the evolutionarily conserved Pro83. The equivalent proline is in the cis conformation in all reported structures of L7Ae and other homologous proteins. Nhp2p has the expected α-β-α fold, but the solution structures of the major conformation of Nhp2p with trans Pro83 and of Nhp2p-S82W with cis Pro83 reveal that Pro83 cis/trans isomerization affects the positions of numerous residues at the Nop10 and RNA binding interface. An S82W substitution, which stabilizes the cis conformation, also stabilizes the association of Nhp2p with H/ACA snoRNPs expressed in vivo. We propose that Pro83 plays a key role in the assembly of the eukaryotic H/ACA RNP, with the cis conformation locking in a stable Cbf5-Nop10-Nhp2 ternary complex and positioning the protein backbone to interact with the H/ACA RNA.