The present view of the mechanism of protein folding

Heat-activated growth of metastable and length-defined DNA fibers expands traditional polymer assembly

Biopolymers such as nucleic acids and proteins exhibit dynamic backbone folding, wherein site-specific intramolecular interactions determine overall structure. Proteins then hierarchically assemble into supramolecular polymers such as microtubules, that are robust yet dynamic, constantly growing or shortening to adjust to cellular needs. The combination of dynamic, energy-driven folding and growth with structural stiffness and length control is difficult to achieve in synthetic polymer self-assembly. Here we show that highly charged, monodisperse DNA-oligomers assemble via seeded growth into length-controlled supramolecular fibers during heating; when the temperature is lowered, these metastable fibers slowly disassemble. Furthermore, the specific molecular structures of oligomers that promote fiber formation contradict the typical theory of block copolymer self-assembly. Efficient curling and packing of the oligomers - or 'curlamers' - determine morphology, rather than hydrophobic to hydrophilic ratio. Addition of a small molecule stabilises the DNA fibers, enabling temporal control of polymer lifetime and underscoring their potential use in nucleic-acid delivery, stimuli-responsive biomaterials, and soft robotics.

The thermodynamics of neurodegenerative disease

The formation of protein aggregates in the brain is a central aspect of the pathology of many neurodegenerative diseases. This self-assembly of specific proteins into filamentous aggregates, or fibrils, is a fundamental biophysical process that can easily be reproduced in the test tube. However, it has been difficult to obtain a clear picture of how the biophysical insights thus obtained can be applied to the complex, multi-factorial diseases and what this means for therapeutic strategies. While new, disease-modifying therapies are now emerging, for the most devastating disorders, such as Alzheimer's and Parkinson's disease, they still fall well short of offering a cure, and few drug design approaches fully exploit the wealth of mechanistic insights that has been obtained in biophysical studies. Here, I attempt to provide a new perspective on the role of protein aggregation in disease, by phrasing the problem in terms of a system that, under constant energy consumption, attempts to maintain a healthy, aggregate-free state against the thermodynamic driving forces that inexorably push it toward pathological aggregation.

Significance of the Disulfide Bridge in the Structure and Stability of Metalloprotein Azurin

Metalloproteins make up a class of proteins that incorporate metal ions into their structures, enabling them to perform essential functions in biological systems, such as catalysis and electron transport. Azurin is one such metalloprotein with copper cofactor, having a β-barrel structure with exceptional thermal stability. The copper metal ion is coordinated at one end of the β-barrel structure, and there is a disulfide bond at the opposite end. In this study, we explore the effect of this disulfide bond in the high thermal stability of azurin by analyzing both the native S-S bonded and S-S nonbonded (S-S open) forms using temperature replica exchange molecular dynamics (REMD). Similar to experimental observations, we find a 35 K decrease in denaturation temperature for S-S open azurin compared to that of the native holo form (420 K). As observed in the case of native holo azurin, the unfolding process of the S-S open form also started with disruptions of the α-helix. The free energy surfaces of the unfolding process revealed that the denaturation event of the S-S open form progresses through different sets of conformational ensembles. Subsequently, we compared the stabilities of individual β-sheet strands of both the S-S bonded and the S-S nonbonded forms of azurin. Further, we examined the contacts between individual residues for the central structures from the free energy surfaces of the S-S nonbonded form. The microscopic origin of the lowering in the denaturation temperature is further supplemented by thermodynamic analysis.

Resolving Hidden Solution Conformations of Hemoglobin Using IMS-IMS on a Cyclic Instrument

Ion mobility spectrometry-mass spectrometry (IMS-MS) experiments on a cyclic IMS instrument were used to examine heterogeneous distributions of structures found in the 15+ to 18+ charge states of the hemoglobin tetramer (Hb). The resolving power of IMS measurements is known to increase with increasing drift-region length. This effect is not significant for Hb charge states as peaks were shown to broaden with increasing drift-region length. This observation suggests that multiple structures with similar cross sections may be present. To examine this hypothesis, selections of drift time distributions were isolated and subsequently reinjected into the mobility region for additional separation. These IMS-IMS experiments demonstrate that selected regions separate further upon additional passes around the drift cell, consistent with the idea that initial resolving power was limited due to the presence of many closely related conformations. Additional variable temperature electrospray ionization (vT-ESI) experiments were conducted to study how changing the solution temperature affects solution conformations. Some features in these IMS-IMS studies were observed to change similarly with solution temperature compared to features in the single IMS distribution. Other features changed differently in the selected mobility data, indicating that solution structures that were obscured upon IMS analysis because of the complex heterogeneity of the original distribution are discernible after reducing the number of conformers that are analyzed by further IMS analysis. These results illustrate that the combination of vT-ESI with IMS-IMS is useful for resolving and exploring conformer distributions and stabilities in systems that exhibit a large degree of structural heterogeneity.

Stability of 20S Proteasome Configurations: Preopening the Axial Gate

Mass spectrometry studies of the stability of the S. cerevisiae 20S proteasome from 11 to 55 °C reveal a series of related configurations and coupled transitions that appear to be associated with opening of the proteolytic core. We find no evidence for dissociation, and all transitions are reversible. A thermodynamic analysis indicates that configurations fall into three general types of structures: enthalpically stabilized, tightly closed (observed as the +54 to +58 charge states) configurations; high-entropy (+60 to +66) states that are proposed as precursors to pore opening; and larger (+70 to +79) partially and fully open pore structures. In the absence of the 19S regulatory unit, the mechanism for opening the 20S pore appears to involve a charge-priming process that loosens the closed-pore configuration. Only a small fraction (≤2%) of these 20S precursor configurations appear to open and thus expose the catalytic cavity.

Role of Metal Cofactor in Enhanced Thermal Stability of Azurin

Metal cofactors are critical centers for different biochemical processes of metalloproteins, and often, this metal coordination renders additional structural stability. In this study, we explore the additional stability conferred by the copper ion on azurin by analyzing both the apo and holo forms using temperature replica exchange molecular dynamics (REMD) data. We find a 14 K decrease in denaturation temperature for apo (406 K) azurin relative to that of holo (420 K), indicating a copper ion-induced additional thermal stability for holo azurin. The unfolding of apo azurin begins with the melting of α-helix and β-sheet V, similar to that of holo form. β-Sheets IV, VII, and VIII are comparatively more stable than other β-strands and melt at higher temperatures. Similar to holo azurin, the strong hydrophobic interactions among the apolar residues in the protein core is the key factor that renders high stability to apo protein as well. We construct free energy surfaces at different temperatures to capture the major conformations along the unfolding basins of the protein. Using contact maps from different basins we show the changes in the interaction between different residues along the unfolding pathway. Furthermore, we compare the Cα root-mean-square fluctuations (Cα-RMSF) and B-factor of all residues of apo and holo forms to understand the flexibility of different regions. The concerted displacement of α-helix and β-sheets V and VI from the protein core is another distinction we observe for apo compared to the holo form, where β-sheet VI was relatively stable.

Disordered proteins mitigate the temperature dependence of site-specific binding free energies

Biophysical characterization of protein-protein interactions involving disordered proteins is challenging. A common simplification is to measure the thermodynamics and kinetics of disordered site binding using peptides containing only the minimum residues necessary. We should not assume, however, that these few residues tell the whole story. Son of sevenless, a multidomain signaling protein from Drosophila melanogaster, is critical to the mitogen-activated protein kinase pathway, passing an external signal to Ras, which leads to cellular responses. The disordered 55 kDa C-terminal domain of Son of sevenless is an autoinhibitor that blocks guanidine exchange factor activity. Activation requires another protein, Downstream of receptor kinase (Drk), which contains two Src homology 3 domains. Here, we utilized NMR spectroscopy and isothermal titration calorimetry to quantify the thermodynamics and kinetics of the N-terminal Src homology 3 domain binding to the strongest sites incorporated into the flanking disordered sequences. Comparing these results to those for isolated peptides provides information about how the larger domain affects binding. The affinities of sites on the disordered domain are like those of the peptides at low temperatures but less sensitive to temperature. Our results, combined with observations showing that intrinsically disordered proteins become more compact with increasing temperature, suggest a mechanism for this effect.

Direct Observation and Real-Time Tracking of an Extraordinarily Stable Folding Intermediate in Mitotic Arrest Deficient Protein 2 Folding by Single-Molecule Fluorescence Resonance Energy Transfer

Although ensemble experiments have suggested that mitotic arrest deficient protein 2 (Mad2), a metamorphic protein, has folding intermediates, direct evidence and characterization are not available. It remains an outstanding challenge to capture the folding intermediates in real time, which is crucial to elucidate the folding mechanism, but the folding intermediates are normally unstable and only exist transiently. By combining confocal-microscopy-based and total internal reflection fluorescence (TIRF)-microscopy-based single-molecule Förster resonance energy transfer (sm-FRET) techniques, we have investigated the folding/unfolding process of Mad2 and captured its folding intermediate at the single-molecule level. This provides direct evidence for the existence of an intermediate along the folding pathway of Mad2. The folding intermediate proved to be extraordinarily stable, with an extremely long average dwell time of 2.3 s under the conditions of 3 M GdmCl at ambient temperature. The folding trajectories obtained from TIRF experiments further suggest that the intermediate is on-pathway to native Mad2.

Amyloid-β Oligomers: Multiple Moving Targets

Alzheimer’s Disease (AD) is a neurodegenerative disorder that is characterized clinically by progressive cognitive decline and pathologically by the β-sheet rich fibril plaque deposition of the amyloid-β (Aβ) peptide in the brain. While plaques are a hallmark of AD, plaque burden is not correlated with cognitive impairment. Instead, Aβ oligomers formed during the aggregation process represent the main agents of neurotoxicity, which occurs 10–20 years before patients begin to show symptoms. These oligomers are dynamic in nature and represented by a heterogeneous distribution of aggregates ranging from low- to high-molecular weight, some of which are toxic while others are not. A major difficulty in determining the pathological mechanism(s) of Aβ, developing reliable diagnostic markers for early-stage detection, as well as effective therapeutics for AD are the differentiation and characterization of oligomers formed throughout disease propagation based on their molecular features, effects on biological function, and relevance to disease propagation and pathology. Thus, it is critical to methodically identify the mechanisms of Aβ aggregation and toxicity, as well as describe the roles of different oligomers and aggregates in disease progression and molecular pathology. Here, we describe a variety of biophysical techniques used to isolate and characterize a range of Aβ oligomer populations, as well as discuss proposed mechanisms of toxicity and therapeutic interventions aimed at specific assemblies formed during the aggregation process. The approaches being used to map the misfolding and aggregation of Aβ are like what was done during the fundamental early studies, mapping protein folding pathways using combinations of biophysical techniques in concert with protein engineering. Such information is critical to the design and molecular engineering of future diagnostics and therapeutics for AD.

Computational Insights of Unfolding of N-Terminal Domain of TDP-43 Reveal the Conformational Heterogeneity in the Unfolding Pathway

TDP-43 proteinopathies is a disease hallmark that characterizes amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). The N-terminal domain of TDP-43 (NTD) is important to both TDP-43 physiology and TDP-43 proteinopathy. However, its folding and dimerization process is still poorly characterized. In the present study, we have investigated the folding/unfolding of NTD employing all-atom molecular dynamics (MD) simulations in 8 M dimethylsulfoxide (DMSO) at high temperatures. The MD results showed that the unfolding of the NTD at high temperature evolves through the formation of a number of conformational states differing in their stability and free energy. The presence of structurally heterogeneous population of intermediate ensembles was further characterized by the different extents of solvent exposure of Trp80 during unfolding. We suggest that these non-natives unfolded intermediate ensembles may facilitate NTD oligomerization and subsequently TDP-43 oligomerization, which might lead to the formation of irreversible pathological aggregates, characteristics of disease pathogenesis.

Sp2- and sp3-C⋯O tetrel bonds in the 3-oxetanone homodimer

The structures and non-covalent interactions at play in the 3-oxetanone homodimer have been investigated using a pulsed jet Fourier transform microwave spectrometer supplemented with quantum chemical calculations. Two isomers were identified in the pulsed jet. With the analyses of non-covalent intermolecular interactions including the quantum theory of atoms, Johnson's non-covalent interactions and natural bond orbital, the observed global minimum is stabilized by a combination of one sp²-C⋯O tetrel bond and a network of multiple C-H⋯O weak hydrogen bonds. The second isomer is characterized by carbonyl-carbonyl interactions, with the formation of one sp²- and one sp³-C⋯O tetrel bond. The conformational population of the two observed isomers in the supersonic expansion was estimated to be N_CE1/N_CC1 ≈ 7/5.

Molecular Insight into the High Thermal Stability of Metalloprotein Azurin

We investigate the events characterizing the steps of the unfolding pathway of blue copper metalloprotein azurin using replica exchange molecular dynamics (REMD). Our studies show that the unfolding of azurin begins with the melting of α-helix and β-sheets II and V. This is followed by the melting of other β-sheets and the exposure of hydrophobic protein core to the solvent, resulting in disruptions of its tertiary structure. Free energy surfaces constructed at different temperatures portray different basins that signify the stability of different melted structures in the unfolding process. The contact maps at different temperatures reveal that the strong hydrophobic interaction within the core of the protein is the vital force that renders high stability to this protein. Analysis of the individual β-sheets by looking into their amino acid sequence shows that β-sheets with charged side chains on the surface melt fast compared to others. The β-barrel of azurin is able to dynamically rearrange, and it helps the protein to preserve its hydrophobic core, holding back the native topology from melting fast. B-factor analysis shows that residues of β-sheets III, IV, and VII deviate less from their initial structure at the transition temperature.

Variable-Temperature Native Mass Spectrometry for Studies of Protein Folding, Stabilities, Assembly, and Molecular Interactions

The structures and conformational dynamics of proteins, protein complexes, and their noncovalent interactions with other molecules are controlled specifically by the Gibbs free energy (entropy and enthalpy) of the system. For some organisms, temperature is highly regulated, but the majority of biophysical studies are carried out at room, nonphysiological temperature. In this review, we describe variable-temperature electrospray ionization (vT-ESI) mass spectrometry (MS)-based studies with unparalleled sensitivity, dynamic range, and selectivity for studies of both cold- and heat-induced chemical processes. Such studies provide direct determinations of stabilities, reactivities, and thermodynamic measurements for native and non-native structures of proteins and protein complexes and for protein-ligand interactions. Highlighted in this review are vT-ESI-MS studies that reveal 40 different conformers of chymotrypsin inhibitor 2, a classic two-state (native → unfolded) unfolder, and thermochemistry for a model membrane protein system binding lipid and its regulatory protein. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Reframing the Protein Folding Problem: Entropy as Organizer

It has been a long-standing conviction that a protein's native fold is selected from a vast number of conformers by the optimal constellation of enthalpically favorable interactions. In marked contrast, this Perspective introduces a different mechanism, one that emphasizes conformational entropy as the principal organizer in protein folding while proposing that the conventional view is incomplete. This mechanism stems from the realization that hydrogen bond satisfaction is a thermodynamic necessity. In particular, a backbone hydrogen bond may add little to the stability of the native state, but a completely unsatisfied backbone hydrogen bond would be dramatically destabilizing, shifting the U(nfolded) ⇌ N(ative) equilibrium far to the left. If even a single backbone polar group is satisfied by solvent when unfolded but buried and unsatisfied when folded, that energy penalty alone, approximately +5 kcal/mol, would rival almost the entire free energy of protein stabilization, typically between -5 and -15 kcal/mol under physiological conditions. Consequently, upon folding, buried backbone polar groups must form hydrogen bonds, and they do so by assembling scaffolds of α-helices and/or strands of β-sheet, the only conformers in which, with rare exception, hydrogen bond donors and acceptors are exactly balanced. In addition, only a few thousand viable scaffold topologies are possible for a typical protein domain. This thermodynamic imperative winnows the folding population by culling conformers with unsatisfied hydrogen bonds, thereby reducing the entropy cost of folding. Importantly, conformational restrictions imposed by backbone···backbone hydrogen bonding in the scaffold are sequence-independent, enabling mutation─and thus evolution─without sacrificing the structure.

Folding Intermediates, Heterogeneous Native Ensembles and Protein Function

Single domain proteins fold via diverse mechanisms emphasizing the intricate relationship between energetics and structure, which is a direct consequence of functional constraints and demands imposed at the level of sequence. On the other hand, elucidating the interplay between folding mechanisms and function is challenging in large proteins, given the inherent shortcomings in identifying metastable states experimentally and the sampling limitations associated with computational methods. Here, we show that free energy profiles and surfaces of large systems (>150 residues), as predicted by a statistical mechanical model, display a wide array of folding mechanisms with ubiquitous folding intermediates and heterogeneous native ensembles. Importantly, residues around the ligand binding or enzyme active site display a larger tendency to partially unfold and this manifests as intermediates or excited states along the folding coordinate in ligand binding domains, transcription repressors, and representative enzymes from all the six classes, including the SARS-CoV-2 receptor binding domain (RBD) of the spike protein and the protease M^pro. It thus appears that it is relatively easier to distill the imprints of function on the folding landscape of larger proteins as opposed to smaller systems. We discuss how an understanding of energetic-entropic features in ordered proteins can pinpoint specific avenues through which folding mechanisms, populations of partially structured states and function can be engineered.

Protein folding - seeing is deceiving

This Perspective is intended to raise questions about the conventional interpretation of protein folding. According to the conventional interpretation, developed over many decades, a protein population can visit a vast number of conformations under unfolding conditions, but a single dominant native population emerges under folding conditions. Accordingly, folding comes with a substantial loss of conformational entropy. How is this price paid? The conventional answer is that favorable interactions between and among the side chains can compensate for entropy loss, and moreover, these interactions are responsible for the structural particulars of the native conformation. Challenging this interpretation, the Perspective introduces a proposal that high energy (i.e., unfavorable) excluding interactions winnow the accessible population substantially under physical-chemical conditions that favor folding. Both steric clash and unsatisfied hydrogen bond donors and acceptors are classified as excluding interactions, so called because conformers with such disfavored interactions will be largely excluded from the thermodynamic population. Both excluding interactions and solvent factors that induce compactness are somewhat nonspecific, yet together they promote substantial chain organization. Moreover, proteins are built on a backbone scaffold consisting of α-helices and strands of β-sheet, where the number of hydrogen bond donors and acceptors is exactly balanced. These repetitive secondary structural elements are the only two conformers that can be both completely hydrogen-bond satisfied and extended indefinitely without encountering a steric clash. Consequently, the number of fundamental folds is limited to no more than ~10,000 for a protein domain. Once excluding interactions are taken into account, the issue of "frustration" is largely eliminated and the Levinthal paradox is resolved. Putting the "bottom line" at the top: it is likely that hydrogen-bond satisfaction represents a largely under-appreciated parameter in protein folding models.

Protons Are Fast and Smart; Proteins Are Slow and Dumb: On the Relationship of Electrospray Ionization Charge States and Conformations

We present simple considerations of how differences in time scales of motions of protons, the lightest and fastest chemical moiety, and the much longer time scales associated with the dynamics of proteins, among the heaviest and slowest analytes, may allow many protein conformations from solution to be kinetically trapped during the process of electrospraying protein solutions into the gas phase. In solution, the quantum nature of protons leads them to change locations by tunneling, an instantaneous process; moreover, the Grotthuss mechanism suggests that these small particles can respond nearly instantaneously to the dynamic motions of proteins that occur on much longer time scales. A conformational change is accompanied by favorable or unfavorable variations in the free energy of the system, providing the impetus for solvent ↔ protein proton exchange. Thus, as thermal distributions of protein conformations interconvert, protonation states rapidly respond, as specific acidic and basic sites are exposed or protected. In the vacuum of the mass spectrometer, protons become immobilized in locations that are specific to the protein conformations from which they were incorporated. In this way, conformational states from solution are preserved upon electrospraying them into the gas phase. These ideas are consistent with the exquisite sensitivity of electrospray mass spectra to small changes of the local environment that alter protein structure in solution. We might remember this approximation for the protonation of proteins in solution with the colloquial expression-protons are fast and smart; proteins are slow and dumb.

Cryptic intermediates and metastable states of proteins as predicted by OneG computational method

Understanding structural excursions of proteins under folding conditions is crucial to map energy landscapes of proteins. In the present study, OneG computational tool has been used for analyzing possible existence of cryptic intermediates and metastable states of 26 proteins for which three prerequisite inputs of the OneG such as atomic coordinates of proteins, free energy of unfolding (ΔG_U) and free energy of exchange (ΔG_HX) determined in the absence of denaturant were available during the course of the study. The veraciousness of the tool on predicting the partially folded states of the proteins has been comprehensively described using experimental data available for 15 of the 26 proteins. Meanwhile, possible existence of partially structured states in the folding pathways of 11 other proteins has also been delineated as predicted by the OneG. In addition to mapping the folding pathways of proteins, the salient merits of the tool on systematically addressing the discrepancy between the ΔG_U and the ΔG_HX of the proteins have also been dealt.Communicated by Ramaswamy H. Sarma.

Rationalizing the Role of Monosodium Glutamate in the Protein Aggregation Through Biophysical Approaches: Potential Impact on Neurodegeneration

Monosodium glutamate (MSG) is the world’s most extensively used food additive and is generally recognized as safe according to the FDA. However, it is well reported that MSG is associated with a number of neurological diseases, and in turn, neurological diseases are associated with protein aggregation. This study rationalized the role of MSG in protein aggregation using different biophysical techniques such as absorption, far-UV CD, DLS, and ITC. Kinetic measurements revealed that MSG causes significant enhancement of aggregation of BSA through a nucleation-dependent polymerization mechanism. Also, CTAB-BSA aggregation is enhanced by MSG significantly. MSG-induced BSA aggregation also exhibits the formation of irreversible aggregates, temperature dependence, non-Arrhenius behavior, and enhancement of hydrodynamic diameter. From the isothermal titration calorimetry measurement, the significant endothermic heat of the interaction of BSA-MSG indicates that protein aggregation may be due to the coupling of MSG with the protein. The determined enthalpy change (ΔH) is largely positive, also suggesting an endothermic nature, whereas entropy change (ΔS) is positive and Gibbs free energy change (ΔG) is largely negative, suggesting the spontaneous nature of the interaction. Furthermore, even a low concentration of MSG is involved in the unfolding of the secondary structure of protein with the disappearance of original peaks and the formation of a unique peak in the far-UV CD, which is an attention-grabbing observation. This is the first investigation which links the dietary MSG with protein aggregation and thus will be very instrumental in understanding the mechanism of various MSG-related human physiological as well as neurological diseases.

The concept of protein folding/unfolding and its impacts on human health

Proteins have evolved in specific 3D structures and play different functions in cells and determine various reactions and pathways. The newly synthesized amino acid chains once depart ribosome must crumple into three-dimensional structures so can be biologically active. This process of protein that makes a functional molecule is called protein folding. The protein folding is both a biological and a physicochemical process that depends on the sequence of it. In fact, this process occurs more complicated and in some cases and in exposure to some molecules like glucose (glycation), mistaken folding leads to amyloid structures and fatal disorders called conformational diseases. Such conditions are detected by the quality control system of the cell and these abnormal proteins undergo renovation or degradation. This scenario takes place by the chaperones, chaperonins, and Ubiquitin-proteasome complex. Understanding of protein folding mechanisms from different views including experimental and computational approaches has revealed some intermediate ensembles such as molten globule and has been subjected to biophysical and molecular biology attempts to know more about prevalent conformational diseases.

Molecular characterization of Hsf1 as a master regulator of heat shock response in the thermotolerant methylotrophic yeast Ogataea parapolymorpha

Ogataea parapolymorpha (Hansenula polymorpha DL-1) is a thermotolerant methylotrophic yeast with biotechnological applications. Here, O. parapolymorpha genes whose expression is induced in response to heat shock were identified by transcriptome analysis and shown to possess heat shock elements (HSEs) in their promoters. The function of O. parapolymorpha HSF1 encoding a putative heat shock transcription factor 1 (OpHsf1) was characterized in the context of heat stress response. Despite exhibiting low sequence identity (26%) to its Saccharomyces cerevisiae homolog, OpHsf1 harbors conserved domains including a DNA binding domain (DBD), domains involved in trimerization (TRI), transcriptional activation (AR1, AR2), transcriptional repression (CE2), and a C-terminal modulator (CTM) domain. OpHSF1 could complement the temperature sensitive (Ts) phenotype of a S. cerevisiae hsf1 mutant. An O. parapolymorpha strain with an H221R mutation in the DBD domain of OpHsf1 exhibited significantly retarded growth and a Ts phenotype. Intriguingly, the expression of heat-shock-protein-coding genes harboring HSEs was significantly decreased in the H221R mutant strain, even under non-stress conditions, indicating the importance of the DBD for the basal growth of O. parapolymorpha. Notably, even though the deletion of C-terminal domains (ΔCE2, ΔAR2, ΔCTM) of OpHsf1 destroyed complementation of the growth defect of the S. cerevisiae hsf1 strain, the C-terminal domains were shown to be dispensable in O. parapolymorpha. Overexpression of OpHsf1 in S. cerevisiae increased resistance to transient heat shock, supporting the idea that OpHsf1 could be useful in the development of heat-shock-resistant yeast host strains.

Evidence for Many Unique Solution Structures for Chymotrypsin Inhibitor 2: A Thermodynamic Perspective Derived from vT-ESI-IMS-MS Measurements

Chymotrypsin inhibitor 2 (CI-2) is a classic model for two-state cooperative protein folding and is one of the most extensively studied systems. Alan Fersht, a pioneer in the field of structural biology, has studied the wild-type (wt) and over 100 mutant forms of CI-2 with traditional analytical and biochemical techniques. Here, we examine wt CI-2 and three mutant forms (A16G, K11A, L32A) to demonstrate the utility of variable-temperature (vT) electrospray ionization (ESI) paired with ion mobility spectrometry (IMS) and mass spectrometry (MS) to map the free energy folding landscape. As the solution temperature is increased, the abundance of each of the six ESI charge states for wt CI-2 and each mutant is found to vary independently. These results require that at least six unique types of CI-2 solution conformers are present. Ion mobility analysis reveals that within each charge state there are additional conformers having distinct solution temperature profiles. A model of the data at ∼30 different temperatures for all four systems suggests the presence of 41 unique CI-2 solution conformations. A thermodynamic analysis of this system yields values of ΔC_p as well as ΔG, ΔH, and ΔS for each state at every temperature studied. Detailed energy landscapes derived from these data provide a rare glimpse into Anfinsen's thermodynamic hypothesis and the process of thermal denaturation, normally thought of as a cooperative two-state transition involving the native state and unstructured denatured species. Specifically, as the temperature is varied, the entropies and enthalpies of different conformers undergo dramatic changes in magnitude and relative order to maintain the delicate balance associated with equilibrium.

Real-time tracking of protein unfolding with time-resolved x-ray solution scattering

The correct folding of proteins is of paramount importance for their function, and protein misfolding is believed to be the primary cause of a wide range of diseases. Protein folding has been investigated with time-averaged methods and time-resolved spectroscopy, but observing the structural dynamics of the unfolding process in real-time is challenging. Here, we demonstrate an approach to directly reveal the structural changes in the unfolding reaction. We use nano- to millisecond time-resolved x-ray solution scattering to probe the unfolding of apomyoglobin. The unfolding reaction was triggered using a temperature jump, which was induced by a nanosecond laser pulse. We demonstrate a new strategy to interpret time-resolved x-ray solution scattering data, which evaluates ensembles of structures obtained from molecular dynamics simulations. We find that apomyoglobin passes three states when unfolding, which we characterize as native, molten globule, and unfolded. The molten globule dominates the population under the conditions investigated herein, whereas native and unfolded structures primarily contribute before the laser jump and 30 μs after it, respectively. The molten globule retains much of the native structure but shows a dynamic pattern of inter-residue contacts. Our study demonstrates a new strategy to directly observe structural changes over the cause of the unfolding reaction, providing time- and spatially resolved atomic details of the folding mechanism of globular proteins.

Exploring the role of hydrophilic amino acids in unfolding of protein in aqueous ethanol solution

Hydrophobic association is the key contributor behind the formation of well packed core of a protein which is often believed to be an important step for folding from an unfolded chain to its compact functional form. While most of the protein folding/unfolding studies have evaluated the changes in the hydrophobic interactions during chemical denaturation, the role of hydrophilic amino acids in such processes are not discussed in detail. Here we report the role of the hydrophilic amino acids behind ethanol induced unfolding of protein. Using free energy simulations, we show that chicken villin head piece (HP-36) protein unfolds gradually in presence of water-ethanol binary mixture with increasing composition of ethanol. However, upon mutation of hydrophilic amino acids by glycine while keeping the hydrophobic amino acids intact, the compact state of the protein is found to be stable at all compositions with gradual flattening of the free energy landscape upon increasing compositions. The local environment around the protein in terms of ethanol/water number significantly differs in wild type protein compared to the mutated protein. The calculated Wyman-Tanford preferential binding coefficient of ethanol for wild type protein reveals that a greater number of cosolutes (here ethanol) bind to the unfolded state compared to its folded state. However, no significant increase in binding coefficient of ethanol at the unfolded state is found for mutated protein. Local-bulk partition coefficient calculation also suggests similar scenarios. Our results reveal that the weakening of hydrophobic interactions in aqueous ethanol solution along with larger preferential binding of ethanol to the unfolded state mediated by hydrophilic amino acids combinedly helps unfolding of protein in aqueous ethanol solution.

London Dispersion and Hydrogen-Bonding Interactions in Bulky Molecules: The Case of Diadamantyl Ether Complexes

Diadamantyl ether (DAE, C20 H30 O) represents a good model to study the interplay between London dispersion and hydrogen-bond interactions. By using broadband rotational spectroscopy, an accurate experimental structure of the diadamantyl ether monomer is obtained and its aggregates with water and a variety of aliphatic alcohols of increasing size are analyzed. In the monomer, C-H⋅⋅⋅H-C London dispersion attractions between the two adamantyl subunits further stabilize its structure. Water and the alcohol partners bind to diadamantyl ether through hydrogen bonding and non-covalent Owater/alcohol ⋅⋅⋅H-CDAE and C-Halcohol ⋅⋅⋅H-CDAE interactions. Electrostatic contributions drive the stabilization of all the complexes, whereas London dispersion interactions become more pronounced with increasing size of the alcohol. Complexes with dominant dispersion contributions are significantly higher in energy and were not observed in the experiment. The results presented herein shed light on the first steps of microsolvation and aggregation of molecular complexes with London dispersion energy donor (DED) groups and the kind of interactions that control them.

Understanding the structure and conformation of bovine hemoglobin in presence of the drug hydroxyurea: multi-spectroscopic studies supported by docking and molecular dynamics simulation

Binding interaction between the small antitumor drug Hydroxyurea (HU) and Bovine Hemoglobin (BHb) has been explored in details from multi-spectroscopic and computational studies. The formation of ground state complex between BHb and HU has been suggested from the electronic UV-Vis and steady-state fluorescence spectroscopic studies. The quenching in fluorescence of BHb in presence of HU at varied concentrations has been analyzed from the SV plots. Static type of quenching has been suggested from time-resolved fluorescence spectroscopic studies. Binding parameters associated with the BHb-HU complex have also been estimated from the temperature dependent fluorescence spectroscopic studies. Alterations in the micro-environment of the Tyr and Trp residues of BHb in presence of HU have been observed from the synchronous fluorescence measurement. The result obtained from CD spectroscopic measurements signify partial unfolding in the secondary structure of BHb due to binding with HU molecule. The experimental observations are supported by theoretical studies. Molecular docking and molecular dynamics simulations have been performed to investigate the structural stability and compactness of BHb in the binding interaction between BHb and HU. The interaction of BHb with HU is expected to provide fundamental insights towards understanding the therapeutic effectiveness of HU upon interaction with BHb used in chemo-, radio therpeutic procedures and also in the treatment of SCD.

Mosaic Cooperativity in Slow Polypeptide Topological Isomerization Revealed by Residue-Specific NMR Thermodynamic Analysis

Slow polypeptide conformational changes on time scales of >1 s are generally assumed to be highly cooperative two-state transitions, reflecting the high energy barrier. However, few experimental characterizations have tested the validity of this assumption. We performed residue-specific NMR thermodynamic analysis of the 27-residue lantibiotic peptide, nukacin ISK-1, to characterize the isomerization between two topological states on the second time scale. Unexpectedly, the thermal transition behaviors were distinct among peptide regions, indicating that the topological isomerization process is a mosaic of different degrees of cooperativity. The conformational change path between the two NMR structures was deduced by a targeted molecular dynamics simulation. The unique side-chain threading motions through the monosulfide rings are the structural basis of the high energy barrier, and the nonlocal interactions in the hydrophobic core are the structural basis of the cooperativity. Taken together, we provide an energetic description of the topological isomerization of nukacin ISK-1.

Limits in accuracy and a strategy of RNA structure prediction using experimental information

RNA structural complexity and flexibility present a challenge for computational modeling efforts. Experimental information and bioinformatics data can be used as restraints to improve the accuracy of RNA tertiary structure prediction. Regarding utilization of restraints, the fundamental questions are: (i) What is the limit in prediction accuracy that one can achieve with arbitrary number of restraints? (ii) Is there a strategy for selection of the minimal number of restraints that would result in the best structural model? We address the first question by testing the limits in prediction accuracy using native contacts as restraints. To address the second question, we develop an algorithm based on the distance variation allowed by secondary structure (DVASS), which ranks restraints according to their importance to RNA tertiary structure prediction. We find that due to kinetic traps, the greatest improvement in the structure prediction accuracy is achieved when we utilize only 40-60% of the total number of native contacts as restraints. When the restraints are sorted by DVASS algorithm, using only the first 20% ranked restraints can greatly improve the prediction accuracy. Our findings suggest that only a limited number of strategically selected distance restraints can significantly assist in RNA structure modeling.

Nature's Shortcut to Protein Folding

This Feature Article presents a view of the protein folding transition based on the hypothesis that Nature has built features within the sequences that enable a Shortcut to efficient folding. Nature's Shortcut is proposed to be the early establishment of a set of nonlocal weak contacts, constituting protein loops that significantly constrain regions of the collapsed disordered protein into a native-like low-resolution fluctuating topology of major sections of the backbone. Nature's establishment of this scaffold of nonlocal contacts is claimed to bypass what would otherwise be a nearly hopeless unaided search for the final three-dimensional structure in proteins longer than ∼100 amino acids. To support this main contention of the Feature Article, the loop hypothesis (LH) description of early folding events is experimentally tested with time-resolved Förster resonance energy transfer techniques for adenylate kinase, and the data are shown to be consistent with theoretical predictions from the sequential collapse model (SCM). The experimentally based LH and the theoretically founded SCM are argued to provide a unified picture of the role of nonlocal contacts as constituting Nature's Shortcut to protein folding. Importantly, the SCM is shown to reliably predict key nonlocal contacts utilizing only primary sequence information. This view on Nature's Shortcut is open to the protein community for further detailed assessment, including its practical consequences, by suitable application of advanced experimental and computational techniques.

Visualizing Protein Folding and Unfolding

Protein folding/unfolding is a complicated process that defies high-resolution characterization by experimental methods. As an alternative, atomistic molecular dynamics simulations are now routinely employed to elucidate and magnify the accompanying conformational changes and the role of solvent in the folding process. However, the level of detail necessary to map the process at high spatial-temporal resolution provides an overwhelming amount of data. As more and better tools are developed for analysis of these large data sets and validation of the simulations, one is still left with the problem of visualizing the results in ways that provide insight into the folding/unfolding process. While viewing and interrogating static crystal structures has become commonplace, more and different approaches are required for dynamic, interconverting, unfolding, and refolding proteins. Here we review a variety of approaches, ranging from straightforward to complex and unintuitive for multiscale analysis and visualization of protein folding and unfolding.

Nonnative contact effects in protein folding

The folding simulations of three ββα-motifs and β-barrel structured proteins (NTL9, NuG2b, and CspA) were performed to determine the important roles of native and nonnative contacts in protein folding.

Solvent Mediation of Peptide Conformations: Polyproline Structures in Water, Methanol, Ethanol, and 1-Propanol as Determined by Ion Mobility Spectrometry-Mass Spectrometry

Ion mobility spectrometry and circular dichroism spectroscopy are used to examine the populations of the small model peptide, polyproline-13 in water, methanol, ethanol, and 1-propanol over a range of solution temperatures (from 288 to 318 K). At low temperatures, the less-polar solvents (1-propanol and ethanol) favor the all-cis polyproline I helix (PPI); as the temperature is increased, the trans-configured polyproline II helix (PPII) is formed. In polar solvents (methanol and water), PPII is favored at all temperatures. From the experimental data, we determine the relative stabilities of the eight structures in methanol, ethanol, and 1-propanol, as well as four in water, all with respect to PPII. Although these conformers show relatively small differences in free energies, substantial variability is observed in the enthalpies and entropies across the structures and solvents. This requires that enthalpies and entropies be highly correlated: in 1-propanol, cis-configured PPI conformations are energetically favorable but entropically disfavored. In more polar solvents, PPI is enthalpically less favorable and entropy favors trans-configured forms. While either ΔH0 or ΔS0 can favor different structures, no conformation in any solvent is simultaneously energetically and entropically stabilized. These data present a rare opportunity to examine the origin of conformational stability. Graphical Abstract ᅟ.

The Hsp70-Hsp90 Chaperone Cascade in Protein Folding

Conserved families of molecular chaperones assist protein folding in the cell. Here we review the conceptual advances on three major folding routes: (i) spontaneous, chaperone-independent folding; (ii) folding assisted by repetitive Hsp70 cycles; and (iii) folding by the Hsp70-Hsp90 cascades. These chaperones prepare their protein clients for folding on their own, without altering their folding path. A particularly interesting role is reserved for Hsp90. The function of Hsp90 in folding is its ancient function downstream of Hsp70, free of cochaperone regulation and present in all kingdoms of life. Eukaryotic signalling networks, however, embrace Hsp90 by a plethora of cochaperones, transforming the profolding machinery to a folding-on-demand factor. We discuss implications for biology and molecular medicine.

Vasoinhibin comprises a three-helix bundle and its antiangiogenic domain is located within the first 79 residues

Vasoinhibin belongs to a family of angiogenesis inhibitors generated when the fourth α-helix (H4) of the hormone prolactin (PRL) is removed by specific proteolytic cleavage. The antiangiogenic properties are absent in uncleaved PRL, indicating that conformational changes create a new bioactive domain. However, the solution structure of vasoinhibin and the location of its bioactive domain are unknown. Molecular dynamic simulation (MD) showed that the loss of H4 exposes the hydrophobic nucleus of PRL and leads to the compression of the molecule into a three-helix bundle that buries the hydrophobic nucleus again. Compression occurs by the movement of loop 1 (L1) and its interaction with α-helix 1 (H1) generating a new L1 conformation with electrostatic and hydrophobic surfaces distinct from those of PRL, that may correspond to a bioactive domain. Consistent with this model, a recombinant protein containing the first 79 amino acids comprising H1 and L1 of human PRL inhibited the proliferation and migration of endothelial cells and upregulated the vasoinhibin target genes, IL1A and ICAM1. This bioactivity was comparable to that of a conventional vasoinhibin having the 123 residues encompassing H1, L1, Η2, L2, and Η3 of human PRL. These findings extend the vasoinhibin family to smaller proteins and provide important structural information, which will aid in antiangiogenic drug development.

Geometrical Description of Protein Structural Motifs

We present a geometrical method that can identify secondary structural motifs in proteins via angular correlations. The method uses crystal structure coordinates to calculate angular and radial signatures of each residue relative to an external reference point as the number of nearest-neighbor residues increases. We apply our approach to the blue copper protein amicyanin using the copper cofactor as the external reference point. We define a signature termed Δβ which describes the change in angular correlation as the span of nearest neighbor residues increases. We find that three turn regions of amicyanin harbor residues with Δβ near zero, while residues in other secondary structures have Δβ greater than zero: for β-strands, Δβ changes gradually residue by residue along the strand. Extension of our analysis to other blue copper proteins demonstrated that the noted structural trends are general. Importantly, our geometrical description of the folded protein accounts for all forces holding the structure together. Through this analysis, we identified some of the turns in amicyanin as symmetrical anchor points.

Elucidation of the structural stability and dynamics of heterogeneous intermediate ensembles in unfolding pathway of the N-terminal domain of TDP-43

The N-terminal domain of the RNA binding protein TDP-43 (NTD) is essential to both physiology and proteinopathy; however, elucidation of its folding/unfolding still remains a major quest. In this study, we have investigated the biophysical behavior of intermediate ensembles employing all-atom molecular dynamics simulations in 8 M urea accelerated with high temperatures to achieve unfolded states in a confined computation time. The cumulative results of the 2.75 μs simulations show that unfolding of the NTD at 350 K evolves through different stable and meta-stable intermediate states. The free-energy landscape reveals two meta-stable intermediates (IN and IU) stabilized by non-native interactions, which are largely hydrophilic and highly energetically frustrated. A single buried tryptophan residue, W80, undergoes solvent exposure to different extents during unfolding; this suggests a structurally heterogeneous population of intermediate ensembles. Furthermore, the structure properties of the IN state show a resemblance to the molten globule (MG) state with most of the secondary structures intact. The unfolding of the NTD is initiated by the loss of β-strands, and the unfolded (U) states exhibit a population of non-native α-helices. These non-native unfolded intermediate ensembles may mediate protein oligomerization, leading to the formation of pathological, irreversible aggregates, characteristics of disease pathogenesis.

From a Highly Disordered to a Metastable State: Uncovering Insights of α-Synuclein

α-Synuclein (αS) is a major constituent of Lewy bodies, the insoluble aggregates that are the hallmark of one of the most prevalent neurodegenerative disorders, Parkinson's disease (PD). The vast majority of experiments in vitro and in vivo provide extensive evidence that a disordered monomeric form is the predominant state of αS in water solution, and it undergoes a large-scale disorder-to-helix transition upon binding to vesicles of different types. Recently, another form, tetrameric, of αS with a stable helical structure was identified experimentally. It has been shown that a dynamic intracellular population of metastable αS tetramers and monomers coexists normally; and the tetramer plays an essential role in maintaining αS homeostasis. Therefore, it is of interest to know whether the tetramer can serve as a means of preventing or delaying the start of PD. Before answering this very important question, it is, first, necessary to find out, on an atomistic level, a correlation between tetramers and monomers; what mediates tetramer formation and what makes a tetramer stable. We address these questions here by investigating both monomeric and tetrameric forms of αS. In particular, by examining correlations between the motions of the side chains and the main chain, steric parameters along the amino-acid sequence, and one- and two-dimensional free-energy landscapes along the coarse-grained dihedral angles γ and δ and principal components, respectively, in monomeric and tetrameric αS, we were able to shed light on a fundamental relationship between monomers and tetramers, and the key residues involved in mediating formation of a tetramer. Also, the reasons for the stability of tetrameric αS and inability of monomeric αS to fold are elucidated here.

Hsp90 Breaks the Deadlock of the Hsp70 Chaperone System

Protein folding in the cell requires ATP-driven chaperone machines such as the conserved Hsp70 and Hsp90. It is enigmatic how these machines fold proteins. Here, we show that Hsp90 takes a key role in protein folding by breaking an Hsp70-inflicted folding block, empowering protein clients to fold on their own. At physiological concentrations, Hsp70 stalls productive folding by binding hydrophobic, core-forming segments. Hsp90 breaks this deadlock and restarts folding. Remarkably, neither Hsp70 nor Hsp90 alters the folding rate despite ensuring high folding yields. In fact, ATP-dependent chaperoning is restricted to the early folding phase. Thus, the Hsp70-Hsp90 cascade does not fold proteins, but instead prepares them for spontaneous, productive folding. This stop-start mechanism is conserved from bacteria to man, assigning also a general function to bacterial Hsp90, HtpG. We speculate that the decreasing hydrophobicity along the Hsp70-Hsp90 cascade may be crucial for enabling spontaneous folding.

Finding the needle in the haystack: towards solving the protein-folding problem computationally

Prediction of protein tertiary structures from amino acid sequence and understanding the mechanisms of how proteins fold, collectively known as "the protein folding problem," has been a grand challenge in molecular biology for over half a century. Theories have been developed that provide us with an unprecedented understanding of protein folding mechanisms. However, computational simulation of protein folding is still difficult, and prediction of protein tertiary structure from amino acid sequence is an unsolved problem. Progress toward a satisfying solution has been slow due to challenges in sampling the vast conformational space and deriving sufficiently accurate energy functions. Nevertheless, several techniques and algorithms have been adopted to overcome these challenges, and the last two decades have seen exciting advances in enhanced sampling algorithms, computational power and tertiary structure prediction methodologies. This review aims at summarizing these computational techniques, specifically conformational sampling algorithms and energy approximations that have been frequently used to study protein-folding mechanisms or to de novo predict protein tertiary structures. We hope that this review can serve as an overview on how the protein-folding problem can be studied computationally and, in cases where experimental approaches are prohibitive, help the researcher choose the most relevant computational approach for the problem at hand. We conclude with a summary of current challenges faced and an outlook on potential future directions.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Conservation of folding and association within a family of spidroin N-terminal domains

Web spiders synthesize silk fibres, nature’s toughest biomaterial, through the controlled assembly of fibroin proteins, so-called spidroins. The highly conserved spidroin N-terminal domain (NTD) is a pH-driven self-assembly device that connects spidroins to super-molecules in fibres. The degree to which forces of self-assembly is conserved across spider glands and species is currently unknown because quantitative measures are missing. Here, we report the comparative investigation of spidroin NTDs originating from the major ampullate glands of the spider species Euprosthenops australis, Nephila clavipes, Latrodectus hesperus, and Latrodectus geometricus. We characterized equilibrium thermodynamics and kinetics of folding and self-association using dynamic light scattering, stopped-flow fluorescence and circular dichroism spectroscopy in combination with thermal and chemical denaturation experiments. We found cooperative two-state folding on a sub-millisecond time scale through a late transition state of all four domains. Stability was compromised by repulsive electrostatic forces originating from clustering of point charges on the NTD surface required for function. pH-driven dimerization proceeded with characteristic fast kinetics yielding high affinities. Results showed that energetics and kinetics of NTD self-assembly are highly conserved across spider species despite the different silk mechanical properties and web geometries they produce.

Thermostable exoshells fold and stabilize recombinant proteins

The expression and stabilization of recombinant proteins is fundamental to basic and applied biology. Here we have engineered a thermostable protein nanoparticle (tES) to improve both expression and stabilization of recombinant proteins using this technology. tES provides steric accommodation and charge complementation to green fluorescent protein (GFPuv), horseradish peroxidase (HRPc), and Renilla luciferase (rLuc), improving the yields of functional in vitro folding by ~100-fold. Encapsulated enzymes retain the ability to metabolize small-molecule substrates, presumably via four 4.5-nm pores present in the tES shell. GFPuv exhibits no spectral shifts in fluorescence compared to a nonencapsulated control. Thermolabile proteins internalized by tES are resistant to thermal, organic, chaotropic, and proteolytic denaturation and can be released from the tES assembly with mild pH titration followed by proteolysis.

Improving recombinant protein expression and stabilization remains a significant challenge. Here, the authors engineer Archaeoglobus fulgidus ferritin as a thermostable exoshell to provide steric accommodation and charge complementation for recombinant proteins, which can improve yields by 100 fold.

Specific and intrinsic sequence patterns extracted by deep learning from intra-protein binding and non-binding peptide fragments

The key finding in the DNA double helix model is the specific pairing or binding between nucleotides A-T and C-G, and the pairing rules are the molecule basis of genetic code. Unfortunately, no such rules have been discovered for proteins. Here we show that intrinsic sequence patterns between intra-protein binding peptide fragments exist, they can be extracted using a deep learning algorithm, and they bear an interesting semblance to the DNA double helix model. The intra-protein binding peptide fragments have specific and intrinsic sequence patterns, distinct from non-binding peptide fragments, and multi-millions of binding and non-binding peptide fragments from currently available protein X-ray structures are classified with an accuracy of up to 93%. The specific binding between short peptide fragments may provide an important driving force for protein folding and protein-protein interaction, two open and fundamental problems in molecular biology, and it may have significant potential in design, discovery, and development of peptide, protein, and antibody drugs.

Watching Proteins Wiggle: Mapping Structures with Two-Dimensional Infrared Spectroscopy

Proteins exhibit structural fluctuations over decades of time scales. From the picosecond side chain motions to aggregates that form over the course of minutes, characterizing protein structure over these vast lengths of time is important to understanding their function. In the past 15 years, two-dimensional infrared spectroscopy (2D IR) has been established as a versatile tool that can uniquely probe proteins structures on many time scales. In this review, we present some of the basic principles behind 2D IR and show how they have, and can, impact the field of protein biophysics. We highlight experiments in which 2D IR spectroscopy has provided structural and dynamical data that would be difficult to obtain with more standard structural biology techniques. We also highlight technological developments in 2D IR that continue to expand the scope of scientific problems that can be accessed in the biomedical sciences.

Plasticity in the Oxidative Folding Pathway of the High Affinity Nerita Versicolor Carboxypeptidase Inhibitor (NvCI)

Nerita Versicolor carboxypeptidase inhibitor (NvCI) is the strongest inhibitor reported so far for the M14A subfamily of carboxypeptidases. It comprises 53 residues and a protein fold composed of a two-stranded antiparallel β sheet connected by three loops and stabilized by three disulfide bridges. Here we report the oxidative folding and reductive unfolding pathways of NvCI. Much debate has gone on whether protein conformational folding guides disulfide bond formation or instead they are disulfide bonds that favour the arrangement of local or global structural elements. We show here that for NvCI both possibilities apply. Under physiological conditions, this protein folds trough a funnelled pathway involving a network of kinetically connected native-like intermediates, all sharing the disulfide bond connecting the two β-strands. In contrast, under denaturing conditions, the folding of NvCI is under thermodynamic control and follows a "trial and error" mechanism, in which an initial quasi-stochastic population of intermediates rearrange their disulfide bonds to attain the stable native topology. Despite their striking mechanistic differences, the efficiency of both folding routes is similar. The present study illustrates thus a surprising plasticity in the folding of this extremely stable small disulfide-rich inhibitor and provides the basis for its redesign for biomedical applications.

Structure-Based Prediction of Protein-Folding Transition Paths

We propose a general theory to describe the distribution of protein-folding transition paths. We show that transition paths follow a predictable sequence of high-free-energy transient states that are separated by free-energy barriers. Each transient state corresponds to the assembly of one or more discrete, cooperative units, which are determined directly from the native structure. We show that the transition state on a folding pathway is reached when a small number of critical contacts are formed between a specific set of substructures, after which folding proceeds downhill in free energy. This approach suggests a natural resolution for distinguishing parallel folding pathways and provides a simple means to predict the rate-limiting step in a folding reaction. Our theory identifies a common folding mechanism for proteins with diverse native structures and establishes general principles for the self-assembly of polymers with specific interactions.