Distinguishing unfolding and functional conformational transitions of calmodulin using ultraviolet resonance Raman spectroscopy

Calmodulin (CaM) is a ubiquitous moderator protein for calcium signaling in all eukaryotic cells. This small calcium‐binding protein exhibits a broad range of structural transitions, including domain opening and folding–unfolding, that allow it to recognize a wide variety of binding partners in vivo. While the static structures of CaM associated with its various binding activities are fairly well‐known, it has been challenging to examine the dynamics of transition between these structures in real‐time, due to a lack of suitable spectroscopic probes of CaM structure. In this article, we examine the potential of ultraviolet resonance Raman (UVRR) spectroscopy for clarifying the nature of structural transitions in CaM. We find that the UVRR spectral change (with 229 nm excitation) due to thermal unfolding of CaM is qualitatively different from that associated with opening of the C‐terminal domain in response to Ca²⁺ binding. This spectral difference is entirely due to differences in tertiary contacts at the interdomain tyrosine residue Tyr138, toward which other spectroscopic methods are not sensitive. We conclude that UVRR is ideally suited to identifying the different types of structural transitions in CaM and other proteins with conformation‐sensitive tyrosine residues, opening a path to time‐resolved studies of CaM dynamics using Raman spectroscopy.     

Nucleation‐based prediction of the protein folding rate and its correlation with the folding nucleus size

The problem of protein self‐organization is in the focus of current molecular biology studies. Although the general principles are understood, many details remain unclear. Specifically, protein folding rates are of interest because they dictate the rate of protein aggregation which underlies many human diseases. Here we offer predictions of protein folding rates and their correlation with folding nucleus sizes. We calculated free energies of the transition state and sizes of folding nuclei for 84 proteins and peptides whose other parameters were measured at the point of thermodynamic equilibrium between their unfolded and native states. We used the dynamic programming method where each residue was considered to be either as folded as in its native state or completely disordered. The calculated and measured folding rates showed a good correlation at the temperature mid‐transition point (the correlation coefficient was 0.75). Also, we pioneered in demonstrating a moderate (‐0.57) correlation coefficient between the calculated sizes of folding nuclei and the folding rates. Predictions made by different methods were compared. The established good correlation between the estimated free energy barrier and the experimentally found folding rate of each studied protein/peptide indicates that our model gives reliable results for the considered data set. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

A novel method to determine thermal transition curves of disulfide-containing proteins and their structured folding intermediates

The stability of a protein or of its folding intermediates is frequently characterized by its resistance to chemical and/or thermal denaturation. The folding/unfolding process is generally followed by spectroscopic methods such as absorbance, fluorescence, circular dichroism spectroscopy, etc. Here, we demonstrate a new method, by using HPLC, for determining the thermal unfolding transitions of disulfide-containing proteins and their structured folding intermediates. The thermal transitions of a model protein, ribonuclease A (RNase A), and a recently found unfolding intermediate of onconase (ONC), des [30-75], have been estimated by this method. Finally, the advantages of this method over traditional techniques are discussed by providing specific examples.     

Shortening a loop can increase protein native state entropy

Protein loops are essential structural elements that influence not only function but also protein stability and folding rates. It was recently reported that shortening a loop in the AcP protein may increase its native state conformational entropy. This effect on the entropy of the folded state can be much larger than the lower entropic penalty of ordering a shorter loop upon folding, and can therefore result in a more pronounced stabilization than predicted by polymer model for loop closure entropy. In this study, which aims at generalizing the effect of loop length shortening on native state dynamics, we use all‐atom molecular dynamics simulations to study how gradual shortening a very long or solvent‐exposed loop region in four different proteins can affect their stability. For two proteins, AcP and Ubc7, we show an increase in native state entropy in addition to the known effect of the loop length on the unfolded state entropy. However, for two permutants of SH3 domain, shortening a loop results only with the expected change in the entropy of the unfolded state, which nicely reproduces the observed experimental stabilization. Here, we show that an increase in the native state entropy following loop shortening is not unique to the AcP protein, yet nor is it a general rule that applies to all proteins following the truncation of any loop. This modification of the loop length on the folded state and on the unfolded state may result with a greater effect on protein stability. Proteins 2015; 83:2137–2146. © 2015 Wiley Periodicals, Inc.     

Protein folding by motion planning

We investigate a novel approach for studying protein folding that has evolved from robotics motion planning techniques called probabilistic roadmap methods (PRMs). Our focus is to study issues related to the folding process, such as the formation of secondary and tertiary structures, assuming we know the native fold. A feature of our PRM-based framework is that the large sets of folding pathways in the roadmaps it produces, in just a few hours on a desktop PC, provide global information about the protein's energy landscape. This is an advantage over other simulation methods such as molecular dynamics or Monte Carlo methods which require more computation and produce only a single trajectory in each run. In our initial studies, we obtained encouraging results for several small proteins. In this paper, we investigate more sophisticated techniques for analyzing the folding pathways in our roadmaps. In addition to more formally revalidating our previous results, we present a case study showing that our technique captures known folding differences between the structurally similar proteins G and L.     

Protein folding at single-molecule resolution

The protein folding reaction carries great significance for cellular function and hence continues to be the research focus of a large interdisciplinary protein science community. Single-molecule methods are providing new and powerful tools for dissecting the mechanisms of this complex process by virtue of their ability to provide views of protein structure and dynamics without associated ensemble averaging. This review briefly introduces common FRET and force methods, and then explores several areas of protein folding where single-molecule experiments have yielded insights. These include exciting new information about folding landscapes, dynamics, intermediates, unfolded ensembles, intrinsically disordered proteins, assisted folding and biomechanical unfolding. Emerging and future work is expected to include advances in single-molecule techniques aimed at such investigations, and increasing work on more complex systems from both the physics and biology standpoints, including folding and dynamics of systems of interacting proteins and of proteins in cells and organisms. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.     

Slow motion analysis of protein folding intermediates within wet silica gels

Resolving the complete folding pathway of a protein is a major challenge to conventional experimental methods because of the rapidity and complexity of folding. Here, we show that entrapment of the protein cytochrome c in wet, optically transparent, porous silica gel matrices has enabled a dramatic expansion, to days or weeks, of the folding time, allowing direct observation of the entire folding pathway using a combination of three spectroscopic techniques, far-ultraviolet circular dichroism, tryptophan fluorescence, and Soret absorption spectroscopy. During refolding in silica gels, collapse and helix formation occur in a stepwise manner, as observed in aqueous solution. Analysis of kinetics and transient spectra indeed reveals a sequence of four distinct intermediates with progressively increasing degrees of folding, two of which closely resemble those previously characterized in solution, namely, the early collapsed and the molten globule intermediates. The other two are the precollapsed and pre-molten globule intermediates that may escape detection by conventional kinetic methods. Interestingly, varying the strength of the gel network has a dramatic effect on the folding time, but no significant effect on the structural features of each folding intermediate, indicating that the interaction between the protein and gel matrix has no measurable effect on the folding pathway. These results better define the major pathway of cytochrome c folding. In addition, in this paper we present the results of the application of this method to a simple, apparent two-state folder ubiquitin, helping to interpret the results for cytochrome c.     

Methods for molecular dynamics simulations of protein folding/unfolding in solution

All atom molecular dynamics simulations have become a standard method for mapping equilibrium protein dynamics and non-equilibrium events like folding and unfolding. Here, we present detailed methods for performing such simulations. Generic protocols for minimization, solvation, simulation, and analysis derived from previous studies are also presented. As a measure of validation, our water model is compared with experiment. An example of current applications of these methods, simulations of the ultrafast folding protein Engrailed Homeodomain are presented including the experimental evidence used to verify their results. Ultrafast folders are an invaluable tool for studying protein behavior as folding and unfolding events measured by experiment occur on timescales accessible with the high-resolution molecular dynamics methods we describe. Finally, to demonstrate the prospect of these methods for folding proteins, a temperature quench simulation of a thermal unfolding intermediate of the Engrailed Homeodomain is described.     

Computational investigation of the effect of thermal perturbation on the mechanical unfolding of titin I27

The emergence of single-molecule force measurement experiments has facilitated a better understanding of protein folding pathways and the thermodynamics involved. Computational methods such as steered molecular dynamics (SMD) simulations are helpful in providing atomistic level information on the unfolding pathways. Recent experimental studies have showed that combinations of single-molecule experiments with traditional methods such as chemical and/or thermal denaturation yield additional insights into the folding phenomenon. In this study, we report results from extensive computations (a total of about 60 SMD simulations with a total length of about 0.4 μs) that address the effect of thermal perturbation on the mechanical stability of the I27 domain of the protein titin. A wide range of temperatures (280-340 K) were considered for the pulling, which was done at both constant velocity and constant force using SMD simulations. Good agreement with experimental data, such as for the trends in changes in average force and the maximum force with respect to the temperature, was obtained. This study identifies two competing pathways for the mechanical unfolding of I27, and illustrates the significance of combining various techniques to examine protein folding.     

Calculations of protein circular dichroism from first principles

Understanding the relationship between the amino acid sequence of a protein and its unique, compact 3D structure is one of the grand challenges in molecular biophysics. One particularly exciting approach is time-resolved electronic circular dichroism (CD) spectroscopy, which offers resolution on a nanosecond (or faster) time scale, although it does not provide the spatial resolution of techniques like X-ray crystallography or NMR. The thrust of our work is to underpin fast time scale spectroscopic studies of protein folding with a stronger theoretical foundation. Ultimately, we seek to use molecular dynamics simulations to study the influence of conformational dynamics and conformational transitions on the electronic CD spectra of proteins. We discuss how improved quantum chemical models of individual chromophores, including aromatic sidechains, can be incorporated into calculations of the electronic structure of proteins and their CD.     

A dynamic model of long-range conformational adaptations triggered by nucleotide binding in GroEL-GroES

The molecular chaperone, GroEL, essential for correct protein folding in E. coli, is composed of 14 identical subunits organized in two interacting rings, each providing a folding chamber for non‐native substrate proteins. The oligomeric assembly shows positive cooperativity within each ring and negative cooperativity between the rings. Although it is well known that ATP and long‐range allosteric interactions drive the functional cycle of GroEL, an atomic resolution view of how ligand binding modulates conformational adaptations over long distances remains a major challenge. Moreover, little is known on the relation between equilibrium dynamics at physiological temperatures and the allosteric transitions in GroEL. Here we present multiple all‐atom molecular dynamics simulations of the GroEL‐GroES assemblies at different stages of the functional cycle. Combined with an extensive analysis of the complete set of experimentally available structures, principal component analysis and conformer plots, we provide an explicit evaluation of the accessible conformational space of unliganded GroEL. Our results suggest the presence of pre‐existing conformers at the equatorial domain level, and a shift of the conformational ensemble upon ATP‐binding. At the inter‐ring interface the simulations capture a remarkable offset motion of helix D triggered by ATP‐binding to the folding active ring. The reorientation of helix D, previously only observed upon GroES association, correlates with a change of the internal dynamics in the opposite ring. This work contributes to the understanding of the molecular mechanisms in GroEL and highlights the ability of all‐atom MD simulations to model long‐range structural changes and allosteric events in large systems. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Early events in protein folding

Recent advances have significantly increased the time and spectroscopic resolution of protein folding experiments. We can now study the timescale and nature of polypeptide collapse, and how this correlates with secondary and tertiary structure formation. Studies on ultrafast folding proteins and peptides provide experimental benchmarks on a timescale that overlaps directly with that of molecular dynamics simulations. This makes possible direct tests of both simulations and current models of protein folding.     

Reducing the dimensionality of the protein-folding search problem

How does a folding protein negotiate a vast, featureless conformational landscape and adopt its native structure in biological real time? Motivated by this search problem, we developed a novel algorithm to compare protein structures. Procedures to identify structural analogs are typically conducted in three-dimensional space: the tertiary structure of a target protein is matched against each candidate in a database of structures, and goodness of fit is evaluated by a distance-based measure, such as the root-mean-square distance between target and candidate. This is an expensive approach because three-dimensional space is complex. Here, we transform the problem into a simpler one-dimensional procedure. Specifically, we identify and label the 11 most populated residue basins in a database of high-resolution protein structures. Using this 11-letter alphabet, any protein''s three-dimensional structure can be transformed into a one-dimensional string by mapping each residue onto its corresponding basin. Similarity between the resultant basin strings can then be evaluated by conventional sequence-based comparison. The disorder → order folding transition is abridged on both sides. At the onset, folding conditions necessitate formation of hydrogen-bonded scaffold elements on which proteins are assembled, severely restricting the magnitude of accessible conformational space. Near the end, chain topology is established prior to emergence of the close-packed native state. At this latter stage of folding, the chain remains molten, and residues populate natural basins that are approximated by the 11 basins derived here. In essence, our algorithm reduces the protein-folding search problem to mapping the amino acid sequence onto a restricted basin string.     

Fluorescence characterization of denatured proteins

Characterization of unfolded states, while critical to a complete understanding of protein folding, is inherently difficult due to structural heterogeneity and dynamic interchange between states. The growing body of work focusing on single molecule fluorescence techniques for the study of protein folding, also highlights their potential for studies of unfolded proteins. These methods can obtain conformational information about individual subpopulations of molecules in an ensemble, and measure dynamics without the need for synchronization. The studies highlighted here demonstrate the promise of these techniques for obtaining novel information about unfolded states in vitro and in more physiologically relevant milieu.     

Complex folding pathways in a simple beta-hairpin

The determination of the folding mechanisms of proteins is critical to understand the topological change that can propagate Alzheimer and Creutzfeld-Jakobs diseases, among others. The computational community has paid considerable attention to this problem; however, the associated time scale, typically on the order of milliseconds or more, represents a formidable challenge. Ab initio protein folding from long molecular dynamics simulations or ensemble dynamics is not feasible with ordinary computing facilities and new techniques must be introduced. Here we present a detailed study of the folding of a 16-residue beta-hairpin, described by a generic energy model and using the activation-relaxation technique. From a total of 90 trajectories at 300 K, three folding pathways emerge. All involve a simultaneous optimization of the complete hydrophobic and hydrogen bonding interactions. The first two pathways follow closely those observed by previous theoretical studies (folding starting at the turn or by interactions between the termini). The third pathway, never observed by previous all-atom folding, unfolding, and equilibrium simulations, can be described as a reptation move of one strand of the beta-sheet with respect to the other. This reptation move indicates that non-native interactions can play a dominant role in the folding of secondary structures. Furthermore, such a mechanism mediated by non-native hydrogen bonds is not available for study by unfolding and Gō model simulations. The exact folding path followed by a given beta-hairpin is likely to be influenced by its sequence and the solvent conditions. Taken together, these results point to a more complex folding picture than expected for a simple beta-hairpin.     

Probing the Folding-Unfolding Transition of a Thermophilic Protein,MTH1880

The folding mechanism of typical proteins has been studied widely, while our understanding of the origin of the high stability of thermophilic proteins is still elusive. Of particular interest is how an atypical thermophilic protein with a novel fold maintains its structure and stability under extreme conditions. Folding-unfolding transitions of MTH1880, a thermophilic protein from Methanobacterium thermoautotrophicum, induced by heat, urea, and GdnHCl, were investigated using spectroscopic techniques including circular dichorism, fluorescence, NMR combined with molecular dynamics (MD) simulations. Our results suggest that MTH1880 undergoes a two-state N to D transition and it is extremely stable against temperature and denaturants. The reversibility of refolding was confirmed by spectroscopic methods and size exclusion chromatography. We found that the hyper-stability of the thermophilic MTH1880 protein originates from an extensive network of both electrostatic and hydrophobic interactions coordinated by the central β-sheet. Spectroscopic measurements, in combination with computational simulations, have helped to clarify the thermodynamic and structural basis for hyper-stability of the novel thermophilic protein MTH1880.     

First-principles calculations of protein circular dichroism in the far-ultraviolet and beyond

Understanding the relationship between the amino acid sequence of a protein and its unique, compact three-dimensional structure is one of the grand challenges in molecular biophysics. One exciting approach to the protein-folding problem is fast time-resolved spectroscopy in the ultra-violet (UV). Time-resolved electronic circular dichroism (CD) spectroscopy offers resolution on a nanosecond (or faster) timescale, but does not provide the spatial resolution of techniques like X-ray crystallography or NMR. There is a need to underpin fast timescale spectroscopic studies of protein folding with a stronger theoretical foundation. We review some recent studies in this regard and briefly highlight how modern quantum chemical models of aromatic groups have improved the accuracy of calculations of protein CD spectra near-UV. On the other side of the far-UV, we describe calculations indicating that charge-transfer transitions are likely to be responsible for bands observed in the vacuum UV in protein CD.