Experimental studies of protein folding and unfolding

MIA40 circumvents the folding constraints imposed by TRIAP1 function

The MIA40 relay system mediates the import of small cysteine-rich proteins into the intermembrane mitochondrial space (IMS). MIA40 substrates are synthesized in the cytosol and assumed to be disordered in their reduced state in this compartment. As they cross the outer mitochondrial membrane, MIA40 promotes the oxidation of critical native disulfides to facilitate folding, trapping functional species in the IMS. Here, we study the redox-controlled folding of TRIAP1, a small cysteine-rich protein with moonlighting function: regulating phospholipid trafficking between mitochondrial membranes in the IMS and preventing apoptosis in the cytosol. TRIAP1 dysregulation is connected to oncogenesis. Although TRIAP1 contains a canonical twin CX9C motif, its sequence characteristics and folding pathway deviate from typical MIA40 substrates. In its reduced state, TRIAP1 rapidly populates a hydrophobic collapsed, alpha-helical, and marginally stable molten globule. This intermediate, biases oxidative folding towards a non-native Cys37-Cys47 kinetic trap, slowing the reaction. MIA40 accelerates TRIAP1 folding rate by 30-fold, bypassing the formation of this folding trap. MIA40 drives the oxidation of the inner disulfide bond Cys18-Cys37, and subsequently, it can catalyze the formation of the outer disulfide bond Cys8-Cys47 to attain the native two-disulfide-bridged structure. We demonstrate that, unlike most MIA40 substrates, TRIAP1's folding pathway is strongly constrained by the structural requirements for its function in phospholipid traffic at the IMS. The obligatory population of a reduced, alpha-helical, metastable molten globule in the cytoplasm may explain TRIAP1's connection to the p53-dependent cell survival pathway, constituting a remarkable example of a functional molten globule state.

The Role of Force Fields and Water Models in Protein Folding and Unfolding Dynamics

Protein folding is a fascinating, not fully understood phenomenon in biology. Molecular dynamics (MD) simulations are an invaluable tool to study conformational changes in atomistic detail, including folding and unfolding processes of proteins. However, the accuracy of the conformational ensembles derived from MD simulations inevitably relies on the quality of the underlying force field in combination with the respective water model. Here, we investigate protein folding, unfolding, and misfolding of fast-folding proteins by examining different force fields with their recommended water models, i.e., ff14SB with the TIP3P model and ff19SB with the OPC model. To this end, we generated long conventional MD simulations highlighting the perks and pitfalls of these setups. Using Markov state models, we defined kinetically independent conformational substates and emphasized their distinct characteristics, as well as their corresponding state probabilities. Surprisingly, we found substantial differences in thermodynamics and kinetics of protein folding, depending on the combination of the protein force field and water model, originating primarily from the different water models. These results emphasize the importance of carefully choosing the force field and the respective water model as they determine the accuracy of the observed dynamics of folding events. Thus, the findings support the hypothesis that the water model is at least equally important as the force field and hence needs to be considered in future studies investigating protein dynamics and folding in all areas of biophysics.

Rethinking the protein folding problem from a new perspective

One of the main concerns of Anfinsen was to reveal the connection between the amino-acid sequence and their biologically active conformation. This search gave rise to two crucial questions in structural biology, namely, why the proteins fold and how a sequence encodes its folding. As to the why, he proposes a plausible answer, namely, the thermodynamic hypothesis. As to the how, this remains an unsolved challenge. Consequently, the protein folding problem is examined here from a new perspective, namely, as an 'analytic whole'. Conceiving the protein folding in this way enabled us to (i) examine in detail why the force-field-based approaches have failed, among other purposes, in their ability to predict the three-dimensional structure of a protein accurately; (ii) propose how to redefine them to prevent these shortcomings, and (iii) conjecture on the origin of the state-of-the-art numerical-methods success to predict the tridimensional structure of proteins accurately.

The Non-native Disulfide-Bond-Containing Landscape Orthogonal to the Oxidative Protein-Folding Trajectory: A Necessary Evil?

Oxidative protein folding describes the process by which disulfide-bond-containing proteins mature from their ribosomal, fully reduced and unfolded, origins. Over the past 40 years, a number of exemplar proteins including bovine pancreatic ribonuclease A (RNaseA), bovine pancreatic trypsin inhibitor (BPTI), and hen egg-white lysozyme (HEWL), among others, have provided rich insight into the nature of the intermolecular interactions that drive the formation of the native, biologically active fold. In this Review Article, we revisit the oxidative folding process of RNase A with a focus on reconciling the role of non-native disulfide-bond-containing species that populate the oxidative folding landscape. Toward gaining such an understanding, we project the regeneration pathway onto a Cartesian coordinate system. This helps not only to recognize the magnitude of the seemingly "fruitless", non-native disulfide-bond-containing species that lie orthogonal to the "native-protein-forming" reaction progress but also to reconcile a role for their existence in the regenerative trajectory. Finally, we superimpose the folding funnel onto the regeneration trajectory to draw parallels between oxidative folders and conformational folders (proteins that lack disulfide bonds). The overall objective is to provide the reader with a semi-quantitative description of oxidative protein folding and the barriers to successful regeneration while underscoring a role of seemingly fruitless intermediates.

Protein folding problem: enigma, paradox, solution

The ability of protein chains to spontaneously form their three-dimensional structures is a long-standing mystery in molecular biology. The most conceptual aspect of this mystery is how the protein chain can find its native, "working" spatial structure (which, for not too big protein chains, corresponds to the global free energy minimum) in a biologically reasonable time, without exhaustive enumeration of all possible conformations, which would take billions of years. This is the so-called "Levinthal's paradox." In this review, we discuss the key ideas and discoveries leading to the current understanding of protein folding kinetics, including folding landscapes and funnels, free energy barriers at the folding/unfolding pathways, and the solution of Levinthal's paradox. A special role here is played by the "all-or-none" phase transition occurring at protein folding and unfolding and by the point of thermodynamic (and kinetic) equilibrium between the "native" and the "unfolded" phases of the protein chain (where the theory obtains the simplest form). The modern theory provides an understanding of key features of protein folding and, in good agreement with experiments, it (i) outlines the chain length-dependent range of protein folding times, (ii) predicts the observed maximal size of "foldable" proteins and domains. Besides, it predicts the maximal size of proteins and domains that fold under solely thermodynamic (rather than kinetic) control. Complementarily, a theoretical analysis of the number of possible protein folding patterns, performed at the level of formation and assembly of secondary structures, correctly outlines the upper limit of protein folding times.

Conformational Stability and Denaturation Processes of Proteins Investigated by Electrophoresis under Extreme Conditions

The functional structure of proteins results from marginally stable folded conformations. Reversible unfolding, irreversible denaturation, and deterioration can be caused by chemical and physical agents due to changes in the physicochemical conditions of pH, ionic strength, temperature, pressure, and electric field or due to the presence of a cosolvent that perturbs the delicate balance between stabilizing and destabilizing interactions and eventually induces chemical modifications. For most proteins, denaturation is a complex process involving transient intermediates in several reversible and eventually irreversible steps. Knowledge of protein stability and denaturation processes is mandatory for the development of enzymes as industrial catalysts, biopharmaceuticals, analytical and medical bioreagents, and safe industrial food. Electrophoresis techniques operating under extreme conditions are convenient tools for analyzing unfolding transitions, trapping transient intermediates, and gaining insight into the mechanisms of denaturation processes. Moreover, quantitative analysis of electrophoretic mobility transition curves allows the estimation of the conformational stability of proteins. These approaches include polyacrylamide gel electrophoresis and capillary zone electrophoresis under cold, heat, and hydrostatic pressure and in the presence of non-ionic denaturing agents or stabilizers such as polyols and heavy water. Lastly, after exposure to extremes of physical conditions, electrophoresis under standard conditions provides information on irreversible processes, slow conformational drifts, and slow renaturation processes. The impressive developments of enzyme technology with multiple applications in fine chemistry, biopharmaceutics, and nanomedicine prompted us to revisit the potentialities of these electrophoretic approaches. This feature review is illustrated with published and unpublished results obtained by the authors on cholinesterases and paraoxonase, two physiologically and toxicologically important enzymes.

De novo protein folding on computers. Benefits and challenges

There has been recent success in prediction of the three-dimensional folded native structures of proteins, most famously by the AlphaFold Algorithm running on Google's/Alphabet's DeepMind computer. However, this largely involves machine learning of protein structures and is not a de novo protein structure prediction method for predicting three-dimensional structures from amino acid residue sequences. A de novo approach would be based almost entirely on general principles of energy and entropy that govern protein folding energetics, and importantly do so without the use of the amino acid sequences and structural features of other proteins. Most consider that problem as still unsolved even though it has occupied leading scientists for decades. Many consider that it remains one of the major outstanding issues in modern science. There is crucial continuing help from experimental findings on protein unfolding and refolding in the laboratory, but only to a limited extent because many researchers consider that the speed by which real proteins folds themselves, often from milliseconds to minutes, is itself still not fully understood. This is unfortunate, because a practical solution to the problem would probably have a major effect on personalized medicine, the pharmaceutical industry, biotechnology, and nanotechnology, including for example "smaller" tasks such as better modeling of flexible "unfolded" regions of the SARS-COV-2 spike glycoprotein when interacting with its cell receptor, antibodies, and therapeutic agents. Some important ideas from earlier studies are given before moving on to lessons from periodic and aperiodic crystals, and a possible role for quantum phenomena. The conclusion is that better computation of entropy should be the priority, though that is presented guardedly.

Structural Analysis of the Black-Legged Tick Saliva Protein Salp15

Salp15 is one of the proteins in the saliva of the tick Ixodes scapularis. Together with other biomolecules injected into the mammalian host at the biting site, it helps the tick to sustain its blood meal for days. Salp15 interferes with the cellular immune response of the mammalian host by inhibiting the activation of CD4⁺ T-lymphocytes. This function is co-opted by pathogens that use the tick as a vector and invade the host when the tick bites, such as Borrelia burgdorferi, the causative agent of Lyme borreliosis. Because of the immunity-suppressing role of Salp15, it has been proposed as a candidate for therapeutic applications in disorders of the immune system. The protein is produced as a 135-residue long polypeptide and secreted without its N-terminal signal 1–21 sequence. Detailed structural studies on Salp15 are lacking because of the difficulty in producing large amounts of the folded protein. We report the production of Salp15 and its structural analysis by NMR. The protein is monomeric and contains a flexible N-terminal region followed by a folded domain with mixed α + β secondary structures. Our results are consistent with a three-dimensional structural model derived from AlphaFold, which predicts the formation of three disulfide bridges and a free C-terminal cysteine.

The Protein Folding Problem: The Role of Theory

The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.

Diselenide crosslinks for enhanced and simplified oxidative protein folding

The in vitro oxidative folding of proteins has been studied for over sixty years, providing critical insight into protein folding mechanisms. Hirudin, the most potent natural inhibitor of thrombin, is a 65-residue protein with three disulfide bonds, and is viewed as a folding model for a wide range of disulfide-rich proteins. Hirudin's folding pathway is notorious for its highly heterogeneous intermediates and scrambled isomers, limiting its folding rate and yield in vitro. Aiming to overcome these limitations, we undertake systematic investigation of diselenide bridges at native and non-native positions and investigate their effect on hirudin's folding, structure and activity. Our studies demonstrate that, regardless of the specific positions of these substitutions, the diselenide crosslinks enhanced the folding rate and yield of the corresponding hirudin analogues, while reducing the complexity and heterogeneity of the process. Moreover, crystal structure analysis confirms that the diselenide substitutions maintained the overall three-dimensional structure of the protein and left its function virtually unchanged. The choice of hirudin as a study model has implications beyond its specific folding mechanism, demonstrating the high potential of diselenide substitutions in the design, preparation and characterization of disulfide-rich proteins.

Flexible Folding: Disulfide-Containing Peptides and Proteins Choose the Pathway Depending on the Environments

In the last few decades, development of novel experimental techniques, such as new types of disulfide (SS)-forming reagents and genetic and chemical technologies for synthesizing designed artificial proteins, is opening a new realm of the oxidative folding study where peptides and proteins can be folded under physiologically more relevant conditions. In this review, after a brief overview of the historical and physicochemical background of oxidative protein folding study, recently revealed folding pathways of several representative peptides and proteins are summarized, including those having two, three, or four SS bonds in the native state, as well as those with odd Cys residues or consisting of two peptide chains. Comparison of the updated pathways with those reported in the early years has revealed the flexible nature of the protein folding pathways. The significantly different pathways characterized for hen-egg white lysozyme and bovine milk α-lactalbumin, which belong to the same protein superfamily, suggest that the information of protein folding pathways, not only the native folded structure, is encoded in the amino acid sequence. The application of the flexible pathways of peptides and proteins to the engineering of folded three-dimensional structures is an interesting and important issue in the new realm of the current oxidative protein folding study.

Chemical Synthesis of Interleukin-2 and Disulfide Stabilizing Analogues

Chemical protein synthesis allows the construction of well-defined structural variations and facilitates the development of deeper understanding of protein structure-function relationships and new protein engineering strategies. Herein, we report the chemical synthesis of interleukin-2 (IL-2) variants on a multimilligram scale and the formation of non-natural disulfide mimetics that improve stability against reduction. The synthesis was accomplished by convergent KAHA ligations; the acidic conditions of KAHA ligation proved to be valuable for the solubilization of the hydrophobic segments of IL-2. The bioactivity of the synthetic IL-2 and its analogues were shown to be equipotent to recombinant IL-2 and exhibit improved stability against reducing agents.

Protein unfolding by SDS: the microscopic mechanisms and the properties of the SDS-protein assembly

The effects of detergent sodium dodecyl sulfate (SDS) on protein structure and dynamics are fundamental to the most common laboratory technique used to separate proteins and determine their molecular weights: polyacrylamide gel electrophoresis. However, the mechanism by which SDS induces protein unfolding and the microstructure of protein-SDS complexes remain largely unknown. Here, we report a detailed account of SDS-induced unfolding of two proteins-I27 domain of titin and β-amylase-obtained through all-atom molecular dynamics simulations. Both proteins were found to spontaneously unfold in the presence of SDS at boiling water temperature on the time scale of several microseconds. The protein unfolding was found to occur via two distinct mechanisms in which specific interactions of individual SDS molecules disrupt the protein's secondary structure. In the final state of the unfolding process, the proteins are found to wrap around SDS micelles in a fluid necklace-and-beads configuration, where the number and location of bound micelles changes dynamically. The global conformation of the protein was found to correlate with the number of SDS micelles bound to it, whereas the number of SDS molecules directly bound to the protein was found to define the relaxation time scale of the unfolded protein. Our microscopic characterization of SDS-protein interactions sets the stage for future refinement of SDS-enabled protein characterization methods, including protein fingerprinting and sequencing using a solid-state nanopore.

Miklós Bodanszky Award Lecture: Selective chalcogen chemistry to study protein science

In recent decades, chemical protein synthesis and the development of chemoselective reactions-including ligation reactions-have led to significant breakthroughs in protein science. Among them are a better understanding of protein structure-function relationships, the study of protein posttranslational modifications, exploration of protein design, unnatural amino acid incorporation, and the study of therapeutic proteins and protein folding. Chalcogen chemistry, especially that of sulfur and selenium, is quite rich, and we have witnessed continuous progress in this field in recent years. In this short review, we will instead summarize three stories that we have recently presented on chalcogen chemistry and its impact on protein science, which was presented in the Miklós Bodanszky Award Lecture at the 35th European Peptide Society Meeting in Dublin, Ireland, 26 August 2018.

BPTI folding revisited: switching a disulfide into methylene thioacetal reveals a previously hidden path

The folding mechanism of the model protein bovine pancreatic trypsin inhibitor was revisited. By switching the solvent exposed disulfide bond with methylene thioacetal we uncovered a hidden pathway in its folding mechanism. In addition, this moiety enhanced protein stability while fully maintaining the protein structure and biological function.

Bovine pancreatic trypsin inhibitor (BPTI) is a 58-residue protein that is stabilized by three disulfide bonds at positions 5–55, 14–38 and 30–51. Widely studied for about 50 years, BPTI represents a folding model for many disulfide-rich proteins. In the study described below, we replaced the solvent exposed 14–38 disulfide bond with a methylene thioacetal bridge in an attempt to arrest the folding pathway of the protein at its two well-known intermediates, N′ and N*. The modified protein was expected to be unable to undergo the rate-determining step in the widely accepted BPTI folding mechanism: the opening of the 14–38 disulfide bond followed by rearrangements that leads to the native state, N. Surprisingly, instead of halting BPTI folding at N′ and N*, we uncovered a hidden pathway involving a direct reaction between the N* intermediate and the oxidizing reagent glutathione (GSSG) to form the disulfide-mixed intermediate N*–SG, which spontaneously folds into N. On the other hand, N′ was unable to fold into N. In addition, we found that the methylene thioacetal bridge enhances BPTI stability while fully maintaining its structure and biological function. These findings suggest a general strategy for enhancing protein stability without compromising on function or structure, suggesting potential applications for future therapeutic protein production.

A component analysis of the free energies of folding of 35 proteins: A consensus view on the thermodynamics of folding at the molecular level

What factors favor protein folding? This is a textbook question. Parsing the experimental free energies of folding/unfolding into diverse enthalpic and entropic components of solute and solvent favoring or disfavoring folding is not an easy task. In this study, we present a computational protocol for estimating the free energy contributors to protein folding semi-quantitatively using ensembles of unfolded and native states generated via molecular dynamics simulations. We tested the methodology on 35 proteins with diverse structural motifs and sizes and found that the calculated free energies correlate well with experiment (correlation coefficient ∼ 0.85), enabling us to develop a consensus view of the energetics of folding. As a more sensitive test of the methodology, we also investigated the free energies of folding of an additional 33 single point mutants and obtained a correlation coefficient of 0.8. A notable observation is that the folding free energy components appear to carry signatures of the fold (SCOP classification) of the protein. © 2017 Wiley Periodicals, Inc.

Distant Phe345 mutation compromises the stability and activity of Mycobacterium tuberculosis isocitrate lyase by modulating its structural flexibility

Isocitrate lyase (ICL), a potential anti-tubercular drug target, catalyzes the first step of the glyoxylate shunt. In the present investigation, we studied the conformational flexibility of MtbICL to better understand its stability and catalytic activity. Our biochemical results showed that a point mutation at Phe345, which is topologically distant (>10 Å) to the active site signature sequence (¹⁸⁹KKCGH¹⁹³), completely abolishes the activity of the enzyme. In depth computational analyses were carried out for understanding the structural alterations using molecular dynamics, time-dependent secondary structure and principal component analysis. The results showed that the mutated residue increased the structural flexibility and induced conformational changes near the active site (residues 170–210) and in the C-terminal lid region (residues 411–428). Both these regions are involved in the catalytic activity of MtbICL. Upon mutation, the residual mobility of the enzyme increased, resulting in a decrease in the stability, which was confirmed by the lower free energy of stabilization in the mutant enzyme suggesting the destabilization in the structure. Our results have both biological importance and chemical novelty. It reveals internal dynamics of the enzyme structure and also suggests that regions other than the active site should be exploited for targeting MtbICL inhibition and development of novel anti-tuberculosis compounds.

There and back again: Two views on the protein folding puzzle

The ability of protein chains to spontaneously form their spatial structures is a long-standing puzzle in molecular biology. Experimentally measured folding times of single-domain globular proteins range from microseconds to hours: the difference (10-11 orders of magnitude) is the same as that between the life span of a mosquito and the age of the universe. This review describes physical theories of rates of overcoming the free-energy barrier separating the natively folded (N) and unfolded (U) states of protein chains in both directions: "U-to-N" and "N-to-U". In the theory of protein folding rates a special role is played by the point of thermodynamic (and kinetic) equilibrium between the native and unfolded state of the chain; here, the theory obtains the simplest form. Paradoxically, a theoretical estimate of the folding time is easier to get from consideration of protein unfolding (the "N-to-U" transition) rather than folding, because it is easier to outline a good unfolding pathway of any structure than a good folding pathway that leads to the stable fold, which is yet unknown to the folding protein chain. And since the rates of direct and reverse reactions are equal at the equilibrium point (as follows from the physical "detailed balance" principle), the estimated folding time can be derived from the estimated unfolding time. Theoretical analysis of the "N-to-U" transition outlines the range of protein folding rates in a good agreement with experiment. Theoretical analysis of folding (the "U-to-N" transition), performed at the level of formation and assembly of protein secondary structures, outlines the upper limit of protein folding times (i.e., of the time of search for the most stable fold). Both theories come to essentially the same results; this is not a surprise, because they describe overcoming one and the same free-energy barrier, although the way to the top of this barrier from the side of the unfolded state is very different from the way from the side of the native state; and both theories agree with experiment. In addition, they predict the maximal size of protein domains that fold under solely thermodynamic (rather than kinetic) control and explain the observed maximal size of the "foldable" protein domains.

Small molecule diselenide additives for in vitro oxidative protein folding

The in vitro oxidative folding of disulfide-rich proteins can be challenging. Here we show a new class of small molecule diselenides, which can be easily prepared from inexpensive starting materials, used to enhance oxidative protein folding. These compounds were tested on a model protein, bovine pancreatic trypsin inhibitor. Two of the tested diselenides showed considerable improvement over glutathione and were on par with the previously described selenoglutathione.

Harnessing selenocysteine reactivity for oxidative protein folding

Turbo-charged folding with selenium: targeted replacement of cysteines in proteins with selenocysteines is a valuable strategy for increasing the rates of oxidative protein folding, altering folding mechanisms, and rescuing kinetically trapped intermediates.

Although oxidative folding of disulfide-rich proteins is often sluggish, this process can be significantly enhanced by targeted replacement of cysteines with selenocysteines. In this study, we examined the effects of a selenosulfide and native versus nonnative diselenides on the folding rates and mechanism of bovine pancreatic trypsin inhibitor. Our results show that such sulfur-to-selenium substitutions alter the distribution of key folding intermediates and enhance their rates of interconversion in a context-dependent manner.

Biochemical and functional characterization of the klotho-VS polymorphism implicated in aging and disease risk

Klotho (KL) is an age-regulating protein named after the Greek goddess who spins the thread of life. Mice deficient in KL are normal throughout development, but rapidly degenerate and display a variety of aging-associated abnormalities that eventually lead to decreased life expectancy. While multiple genetic association studies have identified KL polymorphisms linked with changes in disease risk, there is a paucity of concrete mechanistic data to explain how these amino acid substitutions alter KL protein function. The KLVS polymorphism is suggested to lead to changes in protein trafficking although the mechanism is unclear. Our studies have sought to further investigate the functional differences in the KLVS variant that result in increased risk of many age-related diseases. Our findings suggest that the F352V and C370S substitutions lead to alterations in processing as seen by differences in shedding and half-life. Their co-expression in KLVS results in a phenotype resembling wild-type, but despite this intragenic complementation there are still changes in homodimerization and interactions with FGFR1c. Taken together, these studies suggest that KLVS leads to altered homodimerization that indirectly leads to changes in processing and FGFR1c interactions. These findings help elucidate the functional differences that result from the VS polymorphism, which will help clarify how alterations in KL function can lead to human disease and affect cognition and lifespan.

Unfolding thermodynamics of cysteine-rich proteins and molecular thermal-adaptation of marine ciliates

Euplotes nobilii and Euplotes raikovi are phylogenetically closely allied species of marine ciliates, living in polar and temperate waters, respectively. Their evolutional relation and the sharply different temperatures of their natural environments make them ideal organisms to investigate thermal-adaptation. We perform a comparative study of the thermal unfolding of disulfide-rich protein pheromones produced by these ciliates. Recent circular dichroism (CD) measurements have shown that the two psychrophilic (E. nobilii) and mesophilic (E. raikovi) protein families are characterized by very different melting temperatures, despite their close structural homology. The enhanced thermal stability of the E. raikovi pheromones is realized notwithstanding the fact that these proteins form, as a rule, a smaller number of disulfide bonds. We perform Monte Carlo (MC) simulations in a structure-based coarse-grained (CG) model to show that the higher stability of the E. raikovi pheromones is due to the lower locality of the disulfide bonds, which yields a lower entropy increase in the unfolding process. Our study suggests that the higher stability of the mesophilic E. raikovi phermones is not mainly due to the presence of a strongly hydrophobic core, as it was proposed in the literature. In addition, we argue that the molecular adaptation of these ciliates may have occurred from cold to warm, and not from warm to cold. To provide a testable prediction, we identify a point-mutation of an E. nobilii pheromone that should lead to an unfolding temperature typical of that of E. raikovi pheromones.

Diverse pathways of oxidative folding of disulfide proteins: underlying causes and folding models

The pathway of oxidative folding of disulfide proteins exhibits a high degree of diversity, which is manifested mainly by distinct structural heterogeneity and diverse rearrangement pathways of folding intermediates. During the past two decades, the scope of this diversity has widened through studies of more than 30 disulfide-rich proteins by various laboratories. A more comprehensive landscape of the mechanism of protein oxidative folding has emerged. This review will cover three themes. (1) Elaboration of the scope of diversity of disulfide folding pathways, including the two opposite extreme models, represented by bovine pancreatic trypsin inhibitor (BPTI) and hirudin. (2) Demonstration of experimental evidence accounting for the underlying mechanism of the folding diversity. (3) Discussion of the convergence between the extreme models of oxidative folding and models of conventional conformational folding (framework model, hydrophobic collapse model).

Structural heterogeneity of 6 M GdmCl-denatured proteins: implications for the mechanism of protein folding

An in vitro experiment with protein folding is typically initiated with 6 M GdmCl-denatured proteins, which are generally considered fully unfolded. However, studies conducted by various laboratories have shown that many 6 M GdmCl-denatured proteins are structurally heterogeneous and still retain nativelike residual structures. The extent of conformational heterogeneity of the 6 M GdmCl-denatured protein has significant implications for the folding landscape as well as the interpretation of the observed early stage folding mechanism. Using the method of disulfide scrambling, we are able to gain rough insight into the diverse structural properties of 6 M GdmCl-denatured proteins. It demonstrates that most 6 M GdmCl-denatured proteins are approximately fully denatured, but partially unfolded. Most of them comprise diverse conformational isomers. We review here the cumulative evidence obtained from various laboratories and also provide experimental data obtained in our laboratory.

Isotope-coded, iodoacetamide-based reagent to determine individual cysteine pK(a) values by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry

Cysteine reactivity in enzymes is imparted to a large extent by the stabilization of the deprotonated form of the reduced cysteine (i.e., the thiolate) within the active site. Although this is likely to be an important chemical attribute of many thiol-based enzymes, including cysteine-dependent peroxidases (peroxiredoxins) and proteases, only relatively few pK(a) values have been determined experimentally. Presented here is a new technique for determining the pK(a) value of cysteine residues through quantitative mass spectrometry following chemical modification with an iodoacetamide-based reagent over a range of pH buffers. This isotope-coded reagent, N-phenyl iodoacetamide (iodoacetanilide), is readily prepared in deuterated (d(5)) and protiated (d(0)) versions and is more reactive toward free cysteine than is iodoacetamide. Using this approach, the pK(a) values for the two cysteine residues in Escherichia coli thioredoxin were determined to be 6.5 and greater than 10.0, in good agreement with previous reports using chemical modification approaches. This technique allows the pK(a) of specific cysteine residues to be determined in a clear, fast, and simple manner and, because cysteine residues on separate tryptic peptides are measured separately, is not complicated by the presence of multiple cysteines within the protein of interest.

Self-assembling layers created by membrane proteins on gold

Membrane systems are based on several types of organization. First, amphiphilic lipids are able to create monolayer and bilayer structures which may be flat, vesicular or micellar. Into these structures membrane proteins can be inserted which use the membrane to provide signals for lateral and orientational organization. Furthermore, the proteins are the product of highly specific self-assembly otherwise known as folding, which mostly places individual atoms at precise places in three dimensions. These structures all have dimensions in the nanoscale, except for the size of membrane planes which may extend for millimetres in large liposomes or centimetres on planar surfaces such as monolayers at the air/water interface. Membrane systems can be assembled on to surfaces to create supported bilayers and these have uses in biosensors and in electrical measurements using modified ion channels. The supported systems also allow for measurements using spectroscopy, surface plasmon resonance and atomic force microscopy. By combining the roles of lipids and proteins, highly ordered and specific structures can be self-assembled in aqueous solution at the nanoscale.

Threading a peptide through a peptide: protein loops, rotaxanes, and knots

Proteins adopt complex folds in nature that typically avoid conformations that are knotted or "threaded" through closed loops. Is this the result of fundamental barriers to folding, or have proteins simply evolved to avoid threaded conformations? Organic synthesis has been used in supramolecular chemistry to install topological links in small molecules. By following these principles, we now show that it is possible to assemble a topologically linked protein complex by threading a linear protein through a cyclic protein to form a [2]pseudo-rotaxane. Subsequent ring closure using native chemical ligation cyclizes the linear protein, forming a [2]heterocatenane. Although the kinetics of protein threading are slower than the folding kinetics of the native protein, threading appears to be a highly efficient process.

Unfolding and Refolding of Bovine α-Crystallin in Urea and Its Chaperone Activity

We undertook an unfolding and refolding study of alpha(L)-crystallin in presence of urea to explore the breakdown and formation of various levels of structure and to find out whether the breakdown of various levels of structure occurs simultaneously or in a hierarchal manner. We used various techniques such as circular dichroism, fluorescence spectroscopy, light scattering, polarization to determine the changes in secondary, tertiary, and quaternary structure. Unfolding and refolding occurred through a number of intermediates. The results showed that all levels of structure in alpha(L)-crystallin collapsed or reformed simultaneously. The intermediates that occurred in the 2-4 M urea concentration range during unfolding and refolding differed from each other in terms of the polarity of the tryptophan environment. The ANS binding experiments revealed that refolded alpha(L)-crystallin had higher number of hydrophobic pockets compared to native one. On the other hand, polarity of these pockets remained same as that of the native protein. Both light scattering and polarization measurements showed smaller oligomeric size of refolded alpha(L)-crystallin. Thus, although the secondary structural changes were almost reversible, the tertiary and quaternary structural changes were not. The refolded alpha(L)-crystallin had more exposed hydrophobic sites with increased binding affinity. The refolded form also showed higher chaperone activity than native one. Since the refolded form was smaller in oligomeric size, some buried hydrophobic sites were available. The higher chaperone activity of lower sized oligomer of alpha(L)-crystallin again revealed that chaperone activity was dependent on hydrophobicity and not on oligomeric size.

A ternary nucleation model for the nucleation pathway of protein folding

Recently [Y. S. Djikaev and E. Ruckenstein, J. Phys. Chem. B 111, 886 (2007)], the authors proposed a kinetic model for the nucleation mechanism of protein folding where a protein was modeled as a heteropolymer consisting of hydrophobic and hydrophilic beads and the composition of the growing cluster of protein residues was assumed to be constant and equal to the overall protein composition. Here, they further develop the model by considering a protein as a three-component heteropolymer and by allowing the composition of the growing cluster of protein residues to vary independently of the overall one. All the bonds in the heteropolymer (now consisting of hydrophobic, hydrophilic, and neutral beads) have the same constant length, and all the bond angles are equal and fixed. As a crucial idea of the model, an overall potential around the cluster wherein a residue performs a chaotic motion is considered to be a combination of the average dihedral and average pairwise potentials assigned to the bead. The overall potential as a function of the distance from the cluster center has a double well shape which allows one to determine its emission and absorption rates by using a first passage time analysis. Knowing these rates as functions of three independent variables of a ternary cluster, one can develop a self-consistent kinetic theory for the nucleation mechanism of folding of a protein using a ternary nucleation formalism and evaluate the size and composition of the nucleus and the protein folding time. As an illustration, the model is applied to the folding of bovine pancreatic ribonuclease consisting of 124 amino acids whereof 40 are hydrophobic, 81 hydrophilic, and 3 neutral. With a reasonable choice of diffusion coefficients of the residues in the native state and potential parameters, the model predicts folding times in the range of 1-100 s.

Model for the Nucleation Mechanism of Protein Folding

A nucleation-like pathway of protein folding involves the formation of a cluster containing native residues that grows by including residues from the unfolded part of the protein. This pathway is examined by using a heteropolymer as a protein model. The model heteropolymer consists of hydrophobic and hydrophilic beads with fixed bond lengths and bond angles. The total energy of the heteropolymer is determined by the pairwise repulsive/attractive interactions between nonlinked beads and by the contribution from the dihedral angles involved. The parameters of these interactions can be rigorously defined, unlike the ill-defined surface tension of a cluster of protein residues that constitutes the basis of a previous nucleation model. The main idea underlying the new model consists of averaging the dihedral potential of a selected residue over all possible configurations of all neighboring residues along the protein chain. The resulting average dihedral potential depends on the distance between the selected residue and the cluster center. Its combination with the average pairwise potential of the selected residue and with a confining potential caused by the bonds between the residues leads to an overall potential around the cluster that has a double-well shape. Residues in the inner (closer to the cluster) well are considered as belonging to the folded cluster, whereas those in the outer well are treated as belonging to the unfolded part of the protein. Transitions of residues from the inner well into the outer one and vice versa are considered as elementary emission and absorption events, respectively. The double-well character of the potential well around the cluster allows one to determine the rates of both emission and absorption of residues by the cluster using a first passage time analysis. Once these rates are found as functions of the cluster size, one can develop a self-consistent kinetic theory for the nucleation mechanism of folding of a protein. The model allows one to evaluate the size of the nucleus and the protein folding time. The latter is evaluated as the sum of the times necessary for the first nucleation event to occur and for the nucleus to grow to the maximum size (of the folded protein). Depending on the diffusion coefficients of the native residues in the range from 10(-6) to 10(-8) cm2/s, numerical calculations for a protein of 2500 residues suggest that the folding time ranges from several seconds to several hundreds of seconds.

Time-dependent folding of immunological epitopes of the human chorionic gonadotropin beta-subunit

We have explored the possibility to use 14 different monoclonal antibodies in order to follow the formation of the respective epitopes during the biosynthesis of hCG subunits and their association in JEG-3 choriocarcinoma cells using pulse (30s to 5 min)-chase (0-180 min) experiments. We found central cystine knot epitope structures (epitope beta1) to be formed immediately and simultaneously with epitopes on the protruding hCG-beta loops 1 and 3. We found also differences in the time-dependent folding of beta2 and beta4 epitopes, which are highly overlapping structures on the loops 1+3. These differences were reinforced by decreasing the temperature during the pulse-chase experiments to 25 degrees C. Moreover, we describe for the first time an intracellular intact hCG beta-subunit form that showed the transient expression of the hCG-beta-core fragment epitope beta11 in the course of the maturation of this subunit which casts new light on the presence of hCG-beta-core fragment in Down's syndrome, tumors and pregnancy.

A Stable Disulfide-free Gene-3-protein of Phage fd Generated by In vitro Evolution

Disulfide bonds provide major contributions to the conformational stability of proteins, and their cleavage often leads to unfolding. The gene-3-protein of the filamentous phage fd contains two disulfides in its N1 domain and one in its N2 domain, and these three disulfide bonds are essential for the stability of this protein. Here, we employed in vitro evolution to generate a disulfide-free variant of the N1-N2 protein with a high conformational stability. The gene-3-protein is essential for the phage infectivity, and we exploited this requirement for a proteolytic selection of stabilized protein variants from phage libraries. First, optimal replacements for individual disulfide bonds were identified in libraries, in which the corresponding cysteine codons were randomized. Then stabilizing amino acid replacements at non-cysteine positions were selected from libraries that were created by error-prone PCR. This stepwise procedure led to variants of N1-N2 that are devoid of all three disulfide bonds but stable and functional. The best variant without disulfide bonds showed a much higher conformational stability than the disulfide-containing wild-type form of the gene-3-protein. Despite the loss of all three disulfide bonds, the midpoints of the thermal transitions were increased from 48.5 degrees C to 67.0 degrees C for the N2 domain and from 60.0 degrees C to 78.7 degrees C for the N1 domain. The major loss in conformational stability caused by the removal of the disulfides was thus over-compensated by strongly improved non-covalent interactions. The stabilized variants were less infectious than the wild-type protein, probably because the domain mobility was reduced. Only a small fraction of the sequence space could be accessed by using libraries created by error-prone PCR, but still many strongly stabilized variants could be identified. This is encouraging and indicates that proteins can be stabilized by mutations in many different ways.

The Study of Protein Folding and Dynamics by Determination of Intramolecular Distance Distributions and Their Fluctuations Using Ensemble and Single-Molecule FRET Measurements

The folding and dynamics of globular proteins is a multidimensional problem. The structures of the heterogeneous population of refolding protein molecules are characterized by multiple distances and time constants. Deciphering the mechanism of folding depends on studies of the processes rather than the folded structures alone. Spectroscopy is indispensable for these sorts of studies. Herein, it is shown that the determination of intramolecular distance distributions by ensemble and single-molecule FRET experiments enable the exploration of partially folded states of refolding protein molecules.