Is protein folding hierarchic? II. Folding intermediates and transition states

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

Intrinsically disordered proteins play diverse roles in cell signaling

Abstract

Signaling pathways allow cells to detect and respond to a wide variety of chemical (e.g. Ca²⁺ or chemokine proteins) and physical stimuli (e.g., sheer stress, light). Together, these pathways form an extensive communication network that regulates basic cell activities and coordinates the function of multiple cells or tissues. The process of cell signaling imposes many demands on the proteins that comprise these pathways, including the abilities to form active and inactive states, and to engage in multiple protein interactions. Furthermore, successful signaling often requires amplifying the signal, regulating or tuning the response to the signal, combining information sourced from multiple pathways, all while ensuring fidelity of the process. This sensitivity, adaptability, and tunability are possible, in part, due to the inclusion of intrinsically disordered regions in many proteins involved in cell signaling. The goal of this collection is to highlight the many roles of intrinsic disorder in cell signaling. Following an overview of resources that can be used to study intrinsically disordered proteins, this review highlights the critical role of intrinsically disordered proteins for signaling in widely diverse organisms (animals, plants, bacteria, fungi), in every category of cell signaling pathway (autocrine, juxtacrine, intracrine, paracrine, and endocrine) and at each stage (ligand, receptor, transducer, effector, terminator) in the cell signaling process. Thus, a cell signaling pathway cannot be fully described without understanding how intrinsically disordered protein regions contribute to its function. The ubiquitous presence of intrinsic disorder in different stages of diverse cell signaling pathways suggest that more mechanisms by which disorder modulates intra- and inter-cell signals remain to be discovered.

Graphical abstract

Navigating complex peptide structures using macrocycle conformational maps

Identification of turn motifs that are stabilized by intramolecular hydrogen bonds can be useful in describing the conformation of peptide systems. However, this approach is somewhat insufficient for cyclic peptides because peptide regions that are not positioned within a hydrogen bond can be left with no description. Furthermore, non-regular secondary structures and other rarely-observed conformations can be left without detailed evaluation. Herein, we describe “higher-order” ϕ/ψ plots termed macrocycle conformational maps (MCMs) as a tool for evaluating and comparing the conformations of a series of structurally related macrocyclic peptides.

Identification of turn motifs that are stabilized by hydrogen bonds can be useful in describing the conformation of peptides. Herein, we describe “higher-order” ϕ/ψ plots termed macrocycle conformational maps (MCMs) as a tool to evaluate and compare the conformations of related macrocycles.

Synthesis of the Prototypical Cyclopropyl Dipeptide Mimic and Evaluation of Its Turn-Inducing Capability

The (+) and (-) enantiomers of a new turn-inducing cyclopropyl dipeptide mimic have been synthesized and evaluated. The mimic derives its turn-inducing capabilities solely from the cyclopropyl group and without the conformational biasing that would be provided by side-chain substituents. The mimic and peptide-mimic hybrids prepared from it have been studied using a combination of spectroscopic techniques (NMR, IR, and CD). The dipeptide mimic itself displays intramolecular hydrogen bonding in organic solvents, which differs from that observed in natural peptide turns. In contrast, more elaborate peptide-mimic hybrids exhibit hydrogen bonding characteristics that vary with solvent but are consistent with structures found in the tetrapeptide portion (i → i + 3) of a native β-turn.

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.

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.

Effect of TMAO and Urea on Dimers and Tetramers of Amyloidogenic Heptapeptides (23FGAILSS29)

Human islet amyloid polypeptide (hIAPP) (1-37) is an intrinsically disordered protein that is released with insulin by β-cells found in the pancreas. Under certain environmental conditions, hIAPP can aggregate, which leads to β-cell death. FGAILSS (23-29) residues of the hIAPP protein form β sheets, which may be toxic species in type 2 diabetes (T2D) patients. All-atom molecular dynamics (MD) simulations have been used to analyze the effect of two distinct types of osmolytes trimethylamine N-oxide (TMAO) and urea on two and four FGAILSS heptapeptides. TMAO leads the individual peptide toward an extended conformation with a higher radius of gyration and favors the formation of antiparallel β-sheets with an increase in its concentration. However, urea mostly shows compaction of individual peptides except at 4.0 M in the case of a tetramer but does not show aggregation behavior as a whole. TMAO leads both the dimer and tetramer toward the native state with an increase in its concentration. Moreover, both the dimer and tetramer show irregular behavior in urea. The tetramer in 4.0 M urea shows the maximum fraction of native contacts due to the formation of antiparallel β-sheets. This formation of antiparallel β-sheets favors the aggregation of peptides. TMAO forms a smaller number of hydrogen bonds with peptides as compared to urea as the exclusion of TMAO and accumulation of urea around the peptides have occurred in the first solvation shell (FSS). Principal component analysis (PCA) results suggest that the minima in the free energy landscape (FEL) plot are homogeneous for a particular conformation in TMAO with smaller basins, while in urea, the dimer shows minima mostly for extended conformations. For a 4.0 M urea concentration, the tetramer shows the minimum for antiparallel β-sheets, which indicates the aggregation behavior of the tetramer, and for a higher concentration, it shows minima with wider basins of extended conformations.

Imaging and Force Spectroscopy of Single Transmembrane Proteins with the Atomic Force Microscope

The atomic force microscope (AFM) has opened avenues and provided opportunities to investigate biological soft matter and processes ranging from nanometer (nm) to millimeter (mm). The high temporal (millisecond) and spatial (nanometer) resolutions of the AFM are suited for studying many biological processes in their native conditions. The AFM cantilever-aptly termed as a "lab on a tip"-can be used as an imaging tool as well as a handle to manipulate single bonds and proteins. Recent examples have convincingly established AFM as a tool to study the mechanical properties and monitor processes of single proteins and cells with high sensitivity, thus affording insight into important mechanistic details. This chapter specifically focuses on practical and analytical protocols of single-molecule AFM methodologies related to high-resolution imaging and single-molecule force spectroscopy of transmembrane proteins in a lipid bilayer (reconstituted or native). Both these techniques are operator oriented, and require specialized working knowledge of the instrument, theory and practical skills.

Effects of Arginine Deimination and Citrulline Side-Chain Length on Peptide Secondary Structure Formation

Post-translational modifications expand the chemical functionality of peptides and proteins beyond that originating from the encoded amino acids, but studies on the structural effects of these modifications have been limited. Arginine undergoes deimination to give citrulline (Cit), converting the positively charged guanidinium moiety into a neutral urea group. Herein, we report the effect of Arg deimination on secondary structure formation. To understand the reason for the number of methylene units in Cit, the effect of Cit side-chain length on secondary structure formation was also studied. Ala-based peptides and β-hairpin peptides were used to study α-helix and β-sheet formation, respectively. Peptides containing Cit analogues were prepared by an orthogonal protecting group strategy coupled with solid-phase carbamylation. The CD data for the Ala-based peptides were analyzed by using modified Lifson-Roig theory, showing that the helix propensity of Arg decreased upon deimination and that either shortening or lengthening Cit also decreased the helix propensity. The β-hairpin peptides were analyzed by NMR methods, showing minimal change in strand formation energetics upon Arg deimination. Altering the Cit side-chain length did not affect strand formation energetics either. These results should be useful for the preparation of urea-bearing systems and the design of peptides incorporating urea-bearing residues with varying side-chain length.

Folding with a protein's native shortcut network

A complex network approach to protein folding is proposed, wherein a protein's contact map is reconceptualized as a network of shortcut edges, and folding is steered by a structural characteristic of this network. Shortcut networks are generated by a known message passing algorithm operating on protein residue networks. It is found that the shortcut networks of native structures (SCN0s) are relevant graph objects with which to study protein folding at a formal level. The logarithm form of their contact order (SCN0_lnCO) correlates significantly with folding rate of two-state and nontwo-state proteins. The clustering coefficient of SCN0s (C_SCN0 ) correlates significantly with folding rate, transition-state placement and stability of two-state folders. Reasonable folding pathways for several model proteins are produced when C_SCN0 is used to combine protein segments incrementally to form the native structure. The folding bias captured by C_SCN0 is detectable in non-native structures, as evidenced by Molecular Dynamics simulation generated configurations for the fast folding Villin-headpiece peptide. These results support the use of shortcut networks to investigate the role protein geometry plays in the folding of both small and large globular proteins, and have implications for the design of multibody interaction schemes in folding models. One facet of this geometry is the set of native shortcut triangles, whose attributes are found to be well-suited to identify dehydrated intraprotein areas in tight turns, or at the interface of different secondary structure elements.

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.

Application of an interpretable classification model on Early Folding Residues during protein folding

Background

Machine learning strategies are prominent tools for data analysis. Especially in life sciences, they have become increasingly important to handle the growing datasets collected by the scientific community. Meanwhile, algorithms improve in performance, but also gain complexity, and tend to neglect interpretability and comprehensiveness of the resulting models.

Results

Generalized Matrix Learning Vector Quantization (GMLVQ) is a supervised, prototype-based machine learning method and provides comprehensive visualization capabilities not present in other classifiers which allow for a fine-grained interpretation of the data. In contrast to commonly used machine learning strategies, GMLVQ is well-suited for imbalanced classification problems which are frequent in life sciences. We present a Weka plug-in implementing GMLVQ. The feasibility of GMLVQ is demonstrated on a dataset of Early Folding Residues (EFR) that have been shown to initiate and guide the protein folding process. Using 27 features, an area under the receiver operating characteristic of 76.6% was achieved which is comparable to other state-of-the-art classifiers. The obtained model is accessible at https://biosciences.hs-mittweida.de/efpred/.

Conclusions

The application on EFR prediction demonstrates how an easy interpretation of classification models can promote the comprehension of biological mechanisms. The results shed light on the special features of EFR which were reported as most influential for the classification: EFR are embedded in ordered secondary structure elements and they participate in networks of hydrophobic residues. Visualization capabilities of GMLVQ are presented as we demonstrate how to interpret the results.

Electronic supplementary material

The online version of this article (10.1186/s13040-018-0188-2) contains supplementary material, which is available to authorized users.

Characterizing the relation of functional and Early Folding Residues in protein structures using the example of aminoacyl-tRNA synthetases

Proteins are chains of amino acids which adopt a three-dimensional structure and are then able to catalyze chemical reactions or propagate signals in organisms. Without external influence, many proteins fold into their native structure, and a small number of Early Folding Residues (EFR) have previously been shown to initiate the formation of secondary structure elements and guide their respective assembly. Using the two diverse superfamilies of aminoacyl-tRNA synthetases (aaRS), it is shown that the position of EFR is preserved over the course of evolution even when the corresponding sequence conservation is small. Folding initiation sites are positioned in the center of secondary structure elements, independent of aaRS class. In class I, the predicted position of EFR resembles an ancient structural packing motif present in many seemingly unrelated proteins. Furthermore, it is shown that EFR and functionally relevant residues in aaRS are almost entirely disjoint sets of residues. The Start2Fold database is used to investigate whether this separation of EFR and functional residues can be observed for other proteins. EFR are found to constitute crucial connectors of protein regions which are distant at sequence level. Especially, these residues exhibit a high number of non-covalent residue-residue contacts such as hydrogen bonds and hydrophobic interactions. This tendency also manifests as energetically stable local regions, as substantiated by a knowledge-based potential. Despite profound differences regarding how EFR and functional residues are embedded in protein structures, a strict separation of structurally and functionally relevant residues cannot be observed for a more general collection of proteins.

Revisiting Chameleon Sequences in the Protein Data Bank

The steady growth of the Protein Data Bank (PDB) suggests the periodic repetition of searches for sequences that form different secondary structures in different protein structures; these are called chameleon sequences. This paper presents a fast (nlog(n)) algorithm for such searches and presents the results on all protein structures in the PDB. The longest such sequence found consists of 20 residues.

Compact Structure Patterns in Proteins

Globular proteins typically fold into tightly packed arrays of regular secondary structures. We developed a model to approximate the compact parallel and antiparallel arrangement of α-helices and β-strands, enumerated all possible topologies formed by up to five secondary structural elements (SSEs), searched for their occurrence in spatial structures of proteins, and documented their frequencies of occurrence in the PDB. The enumeration model grows larger super-secondary structure patterns (SSPs) by combining pairs of smaller patterns, a process that approximates a potential path of protein fold evolution. The most prevalent SSPs are typically present in superfolds such as the Rossmann-like fold, the ferredoxin-like fold, and the Greek key motif, whereas the less frequent SSPs often possess uncommon structure features such as split β-sheets, left-handed connections, and crossing loops. This complete SSP enumeration model, for the first time, allows us to investigate which theoretically possible SSPs are not observed in available protein structures. All SSPs with up to four SSEs occurred in proteins. However, among the SSPs with five SSEs, approximately 20% (218) are absent from existing folds. Of these unobserved SSPs, 80% contain two or more uncommon structure features. To facilitate future efforts in protein structure classification, engineering, and design, we provide the resulting patterns and their frequency of occurrence in proteins at: http://prodata.swmed.edu/ssps/.

Cytochrome unfolding pathways from computational analysis of crystal structures

We have developed a model to study the role of geometrical factors in influencing the early stages of unfolding in three cytochromes: cyt c', cyt c-b562 and cyt c. Each stage in unfolding is quantified by the spatial extension λ̂i of n-residue segments, and by their angular extension 〈βn〉. Similarities and differences between and among the three cytochromes in the unfolding of helical and non-helical regions can be determined by analyzing the data for each signature separately. Definite conclusions can be drawn when spatial and angular changes are considered in tandem. To facilitate comparisons, we present graphical portraits of the three cytochromes at the same stage of unfolding, and in relation to their native state structures. We also display specific segments at different stages of unfolding to illustrate differences in stability of defined domains thereby allowing us to make specific predictions on the unfolding of corresponding internal and terminal helices in cyt c' and cyt c-b562. Our work accords with an earlier experimental report on the presence and persistence of a hydrophobic core in cyt c.

279(Val→Phe) Polymorphism of lipoprotein-associated phospholipase A(2) resulted in changes of folding kinetics and recognition to substrate

INTRODUCTION

PLA2G7 encodes Lp-PLA2 having role in the formation of atherosclerotic plaques by catalyzing its substrate, phosphatydilcholine (PC), to be pro-inflammatory substances. The increased risk for coronary artery disease (CAD) in Asian population has been related with this enzyme. 279(Val→Phe) variant was reported to have a protective role against CAD due to, in part, secretion defect or loss of enzymatic function. Therefore, We study folding kinetics and enzyme-substrate interaction in 279(Val→Phe) by using clinical and computational biology approach.

METHODS

Polymorphisms were detected by genotyping among 103 acute myocardial infarction patients and 37 controls. Folding Lp-PLA2 was simulated using GROMACS software by assessing helicity, hydrogen bond formation and stability. The interactions of Lp-PLA2 and its substrate were simulated using Pyrx software followed by molecular dynamics simulation using YASARA software.

RESULT

Polymorphism of 279(Val→Phe) was represented by the change of nucleotide from G to T of 994th PLA2G7 gene. The folding simulation suggested a decreased percentage of α-helix, hydrogen bond formation, hydrogen bond stability and hydrophobicity in 279(Val→Phe). The PC did not interact with active site of 279(Val→Phe) as paradoxically observed in 279 valine. 279(Val→Phe) polymorphism is likely to cause unstable binding to the substrate and decrease the enzymatic activity as observed in molecular dynamics simulations. The results of our computational biology study supported a protected effect of 279(Val→Phe) Polymorphism showed by the odd ratio for MI of 0.22 (CI 95% 0.035-1.37) in this study.

CONCLUSION

279(Val→Phe) Polymorphism of Lp-PLA2 may lead to decrease the enzymatic activity via changes of folding kinetics and recognition to its substrate.

Native globular actin has a thermodynamically unstable quasi-stationary structure with elements of intrinsic disorder

The native form of globular actin, G-actin, is formed in vivo as a result of complex post-translational folding processes that require ATP energy expenditure and are assisted by the 70 kDa heat shock protein, prefoldin and chaperonin containing TCP-1. G-actin is stabilized by the binding of one ATP molecule and one Ca(2+) ion (or Mg(2+) in vivo). Chemical denaturants, heating or Ca(2+) removal transform native actin (N) into 'inactivated actin' (I), a compact oligomer comprising 14-16 subunits. Viscogenic and crowding agents slow this process but do not stop it. The lack of calcium in the solution accelerates the spontaneous N → I transition. Thus, native G-actin has a kinetically stable (as a result of the high free energy barrier between the N and I states) but thermodynamically unstable structure, which, in the absence of Ca(2+) or other bivalent metal ions, spontaneously converts to the thermodynamically stable I state. It was noted that native actin has much in common with intrinsically disordered proteins: it has functionally important disordered regions; it is constantly in complex with one of its numerous partners; and it plays key roles in many cellular processes, in a manner similar to disordered hub proteins. By analyzing actin folding in vivo and unfolding in vitro, we advanced the hypothesis that proteins in a native state may have a thermodynamically unstable quasi-stationary structure. The kinetically stable native state of these proteins appears forcibly under the influence of intracellular folding machinery. The denaturation of such proteins is always irreversible because the inactivated state, for which the structure is determined by the amino acid sequence of a protein, comprises the thermodynamically stable state under physiological conditions.

Effect of charged amino acid side chain length on lateral cross-strand interactions between carboxylate- and guanidinium-containing residues in a β-hairpin

β-Sheet is one of the major protein secondary structures. Oppositely charged residues are frequently observed across neighboring strands in antiparallel sheets, suggesting the importance of cross-strand ion pairing interactions. The charged amino acids Asp, Glu, Arg, and Lys have different numbers of hydrophobic methylenes linking the charged functionality to the backbone. To investigate the effect of side chain length of guanidinium- and carboxylate-containing residues on lateral cross-strand ion pairing interactions at non-hydrogen-bonded positions, β-hairpin peptides containing Zbb-Agx (Zbb = Asp, Glu, Aad in increasing length; Agx = Agh, Arg, Agb, Agp in decreasing length) sequence patterns were studied by NMR methods. The fraction folded population and folding energy were derived from the chemical shift deviation data. Peptides with high fraction folded populations involved charged residue side chain lengths that supported high strand propensity. Double mutant cycle analysis was used to determine the interaction energy for the potential lateral ion pairs. Minimal interaction was observed between residues with short side chains, most likely due to the diffused positive charge on the guanidinium group, which weakened cross-strand electrostatic interactions with the carboxylate side chain. Only the Aad-Arg/Agh interactions with long side chains clearly exhibited stabilizing energetics, possibly relying on hydrophobics. A survey of a non-redundant protein structure database revealed that the statistical sheet pair propensity followed the trend Asp-Arg < Glu-Arg, implying the need for matching long side chains. This suggested the need for long side chains on both guanidinium-bearing and carboxylate-bearing residues to stabilize the β-hairpin motif.

Real-time protein NMR spectroscopy and investigation of assisted protein folding

BACKGROUND

During protein-folding reactions toward the native structure, short-lived intermediate states can be populated. Such intermediates expose hydrophobic patches and can self-associate leading to non-productive protein misfolding. A major focus of current research is the characterization of short-lived intermediates and how molecular chaperones enable productive folding. Real-time NMR spectroscopy, together with the development of advanced methods, is reviewed here and the potential these methods have to characterize intermediate states as well as interactions with molecular chaperone proteins at single-residue resolution is highlighted.

SCOPE OF REVIEW

Various chaperone interactions can guide the protein-folding reaction and thus are important for protein structure formation, stability, and activity of their substrates. Chaperone-assisted protein folding, characterization of intermediates, and their molecular interactions using real-time NMR spectroscopy will be discussed. Additionally, recent advances in NMR methods employed for characterization of high-energy intermediates will be discussed.

MAJOR CONCLUSIONS

Real-time NMR combines high resolution with kinetic information of protein reactions, which can be employed not only for protein-folding studies and the characterization of folding intermediates but also to investigate the molecular mechanisms of assisted protein folding.

GENERAL SIGNIFICANCE

Real-time NMR spectroscopy remains an effective tool to reveal structural details about the interaction between chaperones and transient intermediates. Methodologically, it provides in-depth understanding of how kinetic intermediates and their thermodynamics contribute to the protein-folding reaction. This review summarizes the most recent advances in this field. This article is part of a Special Issue titled Proline-directed Foldases: Cell Signaling Catalysts and Drug Targets.

Actinous enigma or enigmatic actin: Folding, structure, and functions of the most abundant eukaryotic protein

Being the most abundant protein of the eukaryotic cell, actin continues to keep its secrets for more than 60 years. Everything about this protein, its structure, functions, and folding, is mysteriously counterintuitive, and this review represents an attempt to solve some of the riddles and conundrums commonly found in the field of actin research. In fact, actin is a promiscuous binder with a wide spectrum of biological activities. It can exist in at least three structural forms, globular, fibrillar, and inactive (G-, F-, and I-actin, respectively). G-actin represents a thermodynamically instable, quasi-stationary state, which is formed in vivo as a result of the energy-intensive, complex posttranslational folding events controlled and driven by cellular folding machinery. The G-actin structure is dependent on the ATP and Mg²⁺ binding (which in vitro is typically substituted by Ca²⁺) and protein is easily converted to the I-actin by the removal of metal ions and by action of various denaturing agents (pH, temperature, and chemical denaturants). I-actin cannot be converted back to the G-form. Foldable and “natively folded” forms of actin are always involved in interactions either with the specific protein partners, such as Hsp70 chaperone, prefoldin, and the CCT chaperonin during the actin folding in vivo or with Mg²⁺ and ATP as it takes place in the G-form. We emphasize that the solutions for the mysteries of actin multifunctionality, multistructurality, and trapped unfolding can be found in the quasi-stationary nature of this enigmatic protein, which clearly possesses many features attributed to both globular and intrinsically disordered proteins.

Thermodynamics of Deca-alanine Folding in Water

The determination of the folding dynamics of polypeptides and proteins is critical in characterizing their functions in biological systems. Numerous computational models and methods have been developed for studying structure formation at the atomic level. Due to its small size and simple structure, deca-alanine is used as a model system in molecular dynamics (MD) simulations. The free energy of unfolding in vacuum has been studied extensively using the end-to-end distance of the peptide as the reaction coordinate. However, few studies have been conducted in the presence of explicit solvent. Previous results show a significant decrease in the free energy of extended conformations in water, but the α-helical state is still notably favored over the extended state. Although sufficient in vacuum, we show that end-to-end distance is incapable of capturing the full complexity of deca-alanine folding in water. Using α-helical content as a second reaction coordinate, we deduce a more descriptive free-energy landscape, one which reveals a second energy minimum in the extended conformations that is of comparable free energy to the α-helical state. Equilibrium simulations demonstrate the relative stability of the extended and α-helical states in water as well as the transition between the two states. This work reveals both the necessity and challenge of determining a proper reaction coordinate to fully characterize a given process.

Redefining the dry molten globule state of proteins

Dynamics and function of proteins are governed by the structural and energetic properties of the different states they adopt and the barriers separating them. In earlier work, native-state triplet-triplet energy transfer (TTET) on the villin headpiece subdomain (HP35) revealed an equilibrium between a locked native state and an unlocked native state, which are structurally similar but have different dynamic properties. The locked state is restricted to low amplitude motions, whereas the unlocked state shows increased conformational flexibility and undergoes local unfolding reactions. This classified the unlocked state as a dry molten globule (DMG), which was proposed to represent an expanded native state with loosened side-chain interactions and a solvent-shielded core. To test whether the unlocked state of HP35 is actually expanded compared to the locked state, we performed high-pressure TTET measurements. Increasing pressure shifts the equilibrium from the locked toward the unlocked state, with a small negative reaction volume for unlocking (ΔV(0)=-1.6±0.5cm(3)/mol). Therefore, rather than being expanded, the unlocked state represents an alternatively packed, compact state, demonstrating that native proteins can exist in several compact folded states, an observation with implications for protein function. The transition state for unlocking/locking, in contrast, has a largely increased volume relative to the locked and unlocked state, with respective activation volumes of 7.1±0.4cm(3)/mol and 8.7±0.9cm(3)/mol, indicating an expansion of the protein during the locking/unlocking transition. The presented results demonstrate the existence of both compact, low-energy and expanded, high-energy DMGs, prompting a broader definition of this state.

Salt-bridge energetics in halophilic proteins

Halophilic proteins have greater abundance of acidic over basic and very low bulky hydrophobic residues. Classical electrostatic stabilization was suggested as the key determinant for halophilic adaptation of protein. However, contribution of specific electrostatic interactions (i.e. salt-bridges) to overall stability of halophilic proteins is yet to be understood. To understand this, we use Adaptive-Poison-Boltzmann-Solver Methods along with our home-built automation to workout net as well as associated component energy terms such as desolvation energy, bridge energy and background energy for 275 salt-bridges from 20 extremely halophilic proteins. We then perform extensive statistical analysis on general and energetic attributes on these salt-bridges. On average, 8 salt-bridges per 150 residues protein were observed which is almost twice than earlier report. Overall contributions of salt-bridges are −3.0 kcal mol⁻¹. Majority (78%) of salt-bridges in our dataset are stable and conserved in nature. Although, average contributions of component energy terms are equal, their individual details vary greatly from one another indicating their sensitivity to local micro-environment. Notably, 35% of salt-bridges in our database are buried and stable. Greater desolvation penalty of these buried salt-bridges are counteracted by stable network salt-bridges apart from favorable equal contributions of bridge and background terms. Recruitment of extensive network salt-bridges (46%) with a net contribution of −5.0 kcal mol⁻¹ per salt-bridge, seems to be a halophilic design wherein favorable average contribution of background term (−10 kcal mol⁻¹) exceeds than that of bridge term (−7 kcal mol⁻¹). Interiors of proteins from halophiles are seen to possess relatively higher abundance of charge and polar side chains than that of mesophiles which seems to be satisfied by cooperative network salt-bridges. Overall, our theoretical analyses provide insight into halophilic signature in its specific electrostatic interactions which we hope would help in protein engineering and bioinformatics studies.

Fast closure of long loops at the initiation of the folding transition of globular proteins studied by time-resolved FRET-based methods

The protein folding problem would be considered “solved” when it will be possible to “read genes”, i.e., to predict the native fold of proteins, their dynamics, and the mechanism of fast folding based solely on sequence data. The long-term goal should be the creation of an algorithm that would simulate the stepwise mechanism of folding, which constrains the conformational space and in which random search for stable interactions is possible. Here, we focus attention on the initial phases of the folding transition starting with the compact disordered collapsed ensemble, in search of the initial sub-domain structural biases that direct the otherwise stochastic dynamics of the backbone. Our studies are designed to test the “loop hypothesis”, which suggests that fast closure of long loop structures by non-local interactions between clusters of mainly non-polar residues is an essential conformational step at the initiation of the folding transition of globular proteins. We developed and applied experimental methods based on time-resolved resonance excitation energy transfer (trFRET) measurements combined with fast mixing methods and studied the initial phases of the folding of

Effect of charged amino acid side chain length on lateral cross-strand interactions between carboxylate-containing residues and lysine analogues in a β-hairpin

β-Sheets are one of the fundamental three-dimensional building blocks for protein structures. Oppositely charged amino acids are frequently observed directly across one another in antiparallel sheet structures, suggesting the importance of cross-strand ion pairing interactions. Despite the apparent electrostatic nature of ion pairing interactions, the charged amino acids Asp, Glu, Arg, Lys have different numbers of hydrophobic methylenes linking the charged functionality to the backbone. Accordingly, the effect of charged amino acid side chain length on cross-strand ion pairing interactions at lateral non-hydrogen bonded positions was investigated in a β-hairpin motif. The negatively charged residues with a carboxylate (Asp, Glu, Aad in increasing length) were incorporated at position 4, and the positively charged residues with an ammonium (Dap, Dab, Orn, Lys in increasing length) were incorporated at position 9. The fraction folded population and folding free energy were derived from the chemical shift deviation data. Double mutant cycle analysis was used to determine the interaction energy for the potential lateral ion pairs. Only the Asp/Glu-Dap interactions with shorter side chains and the Aad-Orn/Lys interactions with longer side chains exhibited stabilizing energetics, mostly relying on electrostatics and hydrophobics, respectively. This suggested the need for length matching of the interacting residues to stabilize the β-hairpin motif. A survey of a nonredundant protein structure database revealed that the statistical sheet pair propensity followed the trend Asp-Lys < Glu-Lys, also implying the need for length matching of the oppositely charged residues.

Computational and theoretical methods for protein folding

A computational approach is essential whenever the complexity of the process under study is such that direct theoretical or experimental approaches are not viable. This is the case for protein folding, for which a significant amount of data are being collected. This paper reports on the essential role of in silico methods and the unprecedented interplay of computational and theoretical approaches, which is a defining point of the interdisciplinary investigations of the protein folding process. Besides giving an overview of the available computational methods and tools, we argue that computation plays not merely an ancillary role but has a more constructive function in that computational work may precede theory and experiments. More precisely, computation can provide the primary conceptual clues to inspire subsequent theoretical and experimental work even in a case where no preexisting evidence or theoretical frameworks are available. This is cogently manifested in the application of machine learning methods to come to grips with the folding dynamics. These close relationships suggested complementing the review of computational methods within the appropriate theoretical context to provide a self-contained outlook of the basic concepts that have converged into a unified description of folding and have grown in a synergic relationship with their computational counterpart. Finally, the advantages and limitations of current computational methodologies are discussed to show how the smart analysis of large amounts of data and the development of more effective algorithms can improve our understanding of protein folding.

Effect of short- and long-range interactions on trp rotamer populations determined by site-directed tryptophan fluorescence of tear lipocalin

In the lipocalin family, the conserved interaction between the main α-helix and the β-strand H is an ideal model to study protein side chain dynamics. Site-directed tryptophan fluorescence (SDTF) has successfully elucidated tryptophan rotamers at positions along the main alpha helical segment of tear lipocalin (TL). The rotamers assigned by fluorescent lifetimes of Trp residues corroborate the restriction expected based on secondary structure. Steric conflict constrains Trp residues to two (t, g⁻) of three possible χ₁ (t, g⁻, g⁺) canonical rotamers. In this study, investigation focused on the interplay between rotamers for a single amino acid position, Trp 130 on the α-helix and amino acids Val 113 and Leu 115 on the H strand, i.e. long range interactions. Trp130 was substituted for Phe by point mutation (F130W). Mutations at positions 113 and 115 with combinations of Gly, Ala, Phe residues alter the rotamer distribution of Trp130. Mutations, which do not distort local structure, retain two rotamers (two lifetimes) populated in varying proportions. Replacement of either long range partner with a small amino acid, V113A or L115A, eliminates the dominance of the t rotamer. However, a mutation that distorts local structure around Trp130 adds a third fluorescence lifetime component. The results indicate that the energetics of long-range interactions with Trp 130 further tune rotamer populations. Diminished interactions, evident in W130G113A115, result in about a 22% increase of α-helix content. The data support a hierarchic model of protein folding. Initially the secondary structure is formed by short-range interactions. TL has non-native α-helix intermediates at this stage. Then, the long-range interactions produce the native fold, in which TL shows α-helix to β-sheet transitions. The SDTF method is a valuable tool to assess long-range interaction energies through rotamer distribution as well as the characterization of low-populated rotameric states of functionally important excited protein states.

Effect of charged amino acid side chain length at non-hydrogen bonded strand positions on β-hairpin stability

β-Sheets have been implicated in various neurological disorders, and ∼20% of protein residues adopt a sheet conformation. Therefore, studies on the structural origin of sheet stability can provide fundamental knowledge with potential biomedical applications. Oppositely charged amino acids are frequently observed across one another in antiparallel β-sheets. Interestingly, the side chains of natural charged amino acids Asp, Glu, Arg, Lys have different numbers of hydrophobic methylenes linking the backbone to the hydrophilic charged functionalities. To explore the inherent effect of charged amino acid side chain length on antiparallel sheets, the stability of a designed hairpin motif containing charged amino acids with varying side chain lengths at non-hydrogen bonded positions was studied. Peptides with the guest position on the N-terminal strand and the C-terminal strand were investigated by NMR methods. The charged amino acids (Xaa) included negatively charged residues with a carboxylate group (Asp, Glu, Aad in increasing length), positively charged residues with an ammonium group (Dap, Dab, Orn, Lys in increasing length), and positively charged residues with a guanidinium group (Agp, Agb, Arg, Agh in increasing length). The fraction folded and folding free energy for each peptide were derived from the chemical shift deviation data. The stability of the peptides with the charged residues at the N-terminal guest position followed the trends: Asp > Glu > Aad, Dap < Dab < Orn ∼ Lys, and Agb < Arg < Agh < Agp. The stability of the peptides with the charged residues at the C-terminal guest position followed the trends: Asp < Glu < Aad, Dap ∼ Dab < Orn ∼ Lys, and Agb < Arg ∼ Agp < Agh. These trends were rationalized by thermodynamic sheet propensity and cross-strand interactions.

Nonuniform chain collapse during early stages of staphylococcal nuclease folding detected by fluorescence resonance energy transfer and ultrarapid mixing methods

The development of tertiary structure during folding of staphylococcal nuclease (SNase) was studied by time-resolved fluorescence resonance energy transfer measured using continuous- and stopped-flow techniques. Variants of this two-domain protein containing intradomain and interdomain fluorescence donor/acceptor pairs (Trp and Cys-linked fluorophore or quencher) were prepared to probe the intradomain and interdomain structural evolution accompanying SNase folding. The intra-domain donor/acceptor pairs are within the β-barrel domain (Trp27/Cys64 and Trp27/Cys97) and the interdomain pair is between the α-helical domain and the β-barrel domain (Trp140/Cys64). Time-resolved energy transfer efficiency accompanying folding and unfolding at different urea concentrations was measured over a time range from 30 μs to ≈ 10 s. Information on average donor/acceptor distances at different stages of the folding process was obtained by using a quantitative kinetic modeling approach. The average distance for the donor/acceptor pairs in the β-barrel domain decreases to nearly native values whereas that of the interdomain donor/acceptor pairs remains unchanged in the earliest intermediate (<500 μs of refolding). This indicates a rapid nonuniform collapse resulting in an ensemble of heterogeneous conformations in which the central region of the β-barrel domain is well developed while the C-terminal α-helical domain remains disordered. The distance between Trp140 and Cys64 decreases to native values on the 100-ms time scale, indicating that the α-helical domain docks onto the preformed β-barrel at a late stage of the folding. In addition, the unfolded state is found to be more compact under native conditions, suggesting that changes in solvent conditions may induce a nonspecific hydrophobic collapse.

Misfolding and amyloid aggregation of apomyoglobin

Apomyoglobin is an excellent example of a monomeric all α-helical globular protein whose folding pathway has been extensively studied and well characterized. Structural perturbation induced by denaturants or high temperature as well as amino acid substitution have been described to induce misfolding and, in some cases, aggregation. In this article, we review the molecular mechanism of the aggregation process through which a misfolded form of a mutated apomyoglobin aggregates at physiological pH and room temperature forming an amyloid fibril. The results are compared with data showing that either amyloid or aggregate formation occurs under particular denaturing conditions or upon cleavage of the residues corresponding to the C-terminal helix of apomyoglobin. The results are discussed in terms of the sequence regions that are more important than others in determining the amyloid aggregation process.

The loop hypothesis: contribution of early formed specific non-local interactions to the determination of protein folding pathways

The extremely fast and efficient folding transition (in seconds) of globular proteins led to the search for some unifying principles embedded in the physics of the folding polypeptides. Most of the proposed mechanisms highlight the role of local interactions that stabilize secondary structure elements or a folding nucleus as the starting point of the folding pathways, i.e., a "bottom-up" mechanism. Non-local interactions were assumed either to stabilize the nucleus or lead to the later steps of coalescence of the secondary structure elements. An alternative mechanism was proposed, an "up-down" mechanism in which it was assumed that folding starts with the formation of very few non-local interactions which form closed long loops at the initiation of folding. The possible biological advantage of this mechanism, the "loop hypothesis", is that the hydrophobic collapse is associated with ordered compactization which reduces the chance for degradation and misfolding. In the present review the experiments, simulations and theoretical consideration that either directly or indirectly support this mechanism are summarized. It is argued that experiments monitoring the time-dependent development of the formation of specifically targeted early-formed sub-domain structural elements, either long loops or secondary structure elements, are necessary. This can be achieved by the time-resolved FRET-based "double kinetics" method in combination with mutational studies. Yet, attempts to improve the time resolution of the folding initiation should be extended down to the sub-microsecond time regime in order to design experiments that would resolve the classes of proteins which first fold by local or non-local interactions.

A Euclidean Perspective on the Unfolding of Azurin: Spatial Correlations

We investigate the stability to structural perturbation of Pseudomonas aeruginosa azurin using a previously developed geometric model. Our analysis considers Ru(2,2',6',2″-terpyridine)(1,10-phenanthroline)(His83)-labeled wild-type azurin and five variants with mutations to Cu-ligating residues. We find that in the early stages of unfolding, the β-strands exhibit the most structural stability. The conserved residues comprising the hydrophobic core are dislocated only after nearly complete unfolding of the β-barrel. Attachment of the Ru-complex at His83 does not destabilize the protein fold, despite causing some degree of structural rearrangement. Notably, replacing the Cys112 and/or Met121 Cu ligands does not affect the conformational integrity of the protein. Notably, these results are in accord with experimental evidence, as well as molecular dynamics simulations of the denaturation of azurin.

Atomic force microscopy and spectroscopy to probe single membrane proteins in lipid bilayers

The atomic force microscope (AFM) has opened vast avenues hitherto inaccessible to the biological scientist. The high temporal (millisecond) and spatial (nanometer) resolutions of the AFM are suited for studying many biological processes in their native conditions. The AFM cantilever stylus is aptly termed as a "lab on a tip" owing to its versatility as an imaging tool as well as a handle to manipulate single bonds and proteins. Recent examples assert that the AFM can be used to study the mechanical properties and monitor processes of single proteins and single cells, thus affording insight into important mechanistic details. This chapter specifically focuses on practical and analytical protocols of single-molecule AFM methodologies related to high-resolution imaging and single-molecule force spectroscopy of membrane proteins. Both these techniques are operator oriented, and require specialized working knowledge of the instrument, theoretical, and practical skills.

Kinetic Analysis of Guanidine Hydrochloride Inactivation of β-Galactosidase in the Presence of Galactose

Inactivation of purified β-Galactosidase was done with GdnHCl in the absence and presence of varying [galactose] at 50°C and at pH 4.5. Lineweaver-Burk plots of initial velocity data, in the presence and absence of guanidine hydrochloride (GdnHCl) and galactose, were used to determine the relevant K_m and V_max values, with p-nitrophenyl β-D-galactopyranoside (pNPG) as substrate, S. Plots of ln([P]_∞ − [P]_t) against time in the presence of GdnHCl yielded the inactivation rate constant, A. Plots of A versus [S] at different galactose concentrations were straight lines that became increasingly less steep as the [galactose] increased, showing that A was dependent on [S]. Slopes and intercepts of the 1/[P]_∞ versus 1/[S] yielded k₊₀ and k'₊₀, the microscopic rate constants for the free enzyme and the enzyme-substrate complex, respectively. Plots of k₊₀ and k'₊₀ versus [galactose] showed that galactose protected the free enzyme as well as the enzyme-substrate complex (only at the lowest and highest [galactose]) against GdnHCl inactivation. In the absence of galactose, GdnHCl exhibited some degree of non-competitive inhibition. In the presence of GdnHCl, galactose exhibited competitive inhibition at the lower [galactose] of 5 mM which changed to non-competitive as the [galactose] increased. The implications of our findings are further discussed.