Comparing Fast Pressure Jump and Temperature Jump Protein Folding Experiments and Simulations

Kinetics and dynamics of oligonucleotide hybridization

The hybridization of short nucleic acid strands is a remarkable spontaneous process that is foundational to biotechnology and nanotechnology and plays a crucial role in gene expression, editing and DNA repair. Decades of research into the mechanism of hybridization have resulted in a deep understanding of its thermodynamics, but many questions remain regarding its kinetics and dynamics. Recent advances in experiments and molecular dynamics simulations of nucleic acids are enabling more direct insight into the structural dynamics of hybridization, which can test long-standing assumptions regarding its mechanism. In this Review, we summarize the current state of knowledge of hybridization kinetics, discuss the barriers to a molecular description of hybridization dynamics, and highlight the new approaches that have begun uncovering the dynamics of hybridization and the duplex ensemble. The kinetics and dynamics of hybridization are highly sensitive to the composition of nucleic acids, and we emphasize recent discoveries and open questions on the role of nucleobase sequence and chemical modifications.

A Fluorescence-Based Temperature-Jump Apparatus for Illustrating Protein Dynamics on the Millisecond Time Scale

A fluorescence-based temperature jump (T-jump) module was constructed to illustrate the large-domain motion of a given protein upon thermal stimulus on the millisecond time scale. The aqueous sample was readily heated by 5.0 °C in ca. 2 ms with a lasting high temperature plateau (>1 s) upon irradiation with the "optical Riemann sum" of the discrete infrared pulses of different energy sequences from a 1467 nm diode laser operated at 1k Hz. The temperature evolution was revealed by the time-evolved fluorescence intensity change of the dissolved tryptophan. Bovine serum albumin (BSA) and human serum albumin (HSA) were chosen as model proteins, and their fluorescence intensity evolutions were recorded at 36.6-39.9 °C upon T-jump from 35.0 °C, within the range of physiological temperatures. The observed protein dynamics of BSA was characterized with an apparent activation energy of 276 ± 23 kJ mol-1, whereas HSA did not manifest the dynamic component. In this measurement, only a tiny amount of sample, ca. 1 μL, was required due to the conjugation of the microspot objective, and the initial temperature was readily controlled by a homemade thermostatic pad. This millisecond-resolution technique is advantageous for illustrating the large-domain dynamics of the targeted protein, bridging the characterizations of the localized protein dynamics on nanosecond to microsecond time scales using the fast techniques and the steady-state protein conformational features by conventional methods, such as Fourier-transform infrared and circular dichroism spectroscopies.

Effects of altered backbone composition on the folding kinetics and mechanism of an ultrafast-folding protein

Sequence-encoded protein folding is a ubiquitous biological process that has been successfully engineered in a range of oligomeric molecules with artificial backbone chemical connectivity. A remarkable aspect of protein folding is the contrast between the rapid rates at which most sequences in nature fold and the vast number of conformational states possible in an unfolded chain with hundreds of rotatable bonds. Research efforts spanning several decades have sought to elucidate the fundamental chemical principles that dictate the speed and mechanism of natural protein folding. In contrast, little is known about how protein mimetic entities transition between an unfolded and folded state. Here, we report effects of altered backbone connectivity on the folding kinetics and mechanism of the B domain of Staphylococcal protein A (BdpA), an ultrafast-folding sequence. A combination of experimental biophysical analysis and atomistic molecular dynamics simulations performed on the prototype protein and several heterogeneous-backbone variants reveal the interplay among backbone flexibility, folding rates, and structural details of the transition state ensemble. Collectively, these findings suggest a significant degree of plasticity in the mechanisms that can give rise to ultrafast folding in the BdpA sequence and provide atomic level insights into how protein mimetic chains adopt an ordered folded state.

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.

Force-Dependent Folding Kinetics of Single Molecules with Multiple Intermediates and Pathways

Most single-molecule studies derive the kinetic rates of native, intermediate, and unfolded states from equilibrium hopping experiments. Here, we apply the Kramers kinetic diffusive model to derive the force-dependent kinetic rates of intermediate states from nonequilibrium pulling experiments. From the kinetic rates, we also extract the force-dependent kinetic barriers and the equilibrium folding energies. We apply our method to DNA hairpins with multiple folding pathways and intermediates. The experimental results agree with theoretical predictions. Furthermore, the proposed nonequilibrium single-molecule approach permits us to characterize kinetic and thermodynamic properties of native, unfolded, and intermediate states that cannot be derived from equilibrium hopping experiments.

Temperature and Guanidine Hydrochloride Effects on the Folding Thermodynamics of WW Domain and Variants

We used simulations based on an all-atom Go model to calculate the folding temperatures (T_fs) and free energies (ΔGs) of two variants of the WW domain, which is a small all-β-sheet protein. The results, without adjusting any parameter, are in good agreement with experiments, thus validating the simulations. We then used the molecular transfer model to predict the changes in the ΔGs and T_fs as the guanidine hydrochloride concentration is varied. The predictions can be readily tested in experiments.

Tandem surface-induced dissociation of protein complexes on an ultrahigh resolution platform

We describe instrumentation for conducting tandem surface-induced dissociation (tSID) of native protein complexes on an ultrahigh resolution Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. The two stages of SID are accomplished with split lenses replacing the entrance lenses of the quadrupole mass filter (stage 1, referred to herein as SID-Q) and the collision cell (stage 2, Q-SID). After SID-Q, the scattered projectile ions and subcomplexes formed in transit traverse the 20 mm pre-filter prior to the mass-selecting quadrupole, providing preliminary insights into the SID fragmentation kinetics of noncovalent protein complexes. The isolated SID fragments (subcomplexes) are then fragmented by SID in the collision cell entrance lens (Q-SID), generating subcomplexes of subcomplexes. We show that the ultrahigh resolution of the FT-ICR can be used for deconvolving species overlapping in m/z, which are particularly prominent in tandem SID spectra due to the combination of symmetric charge partitioning and narrow product ion charge state distributions. Various protein complex topologies are explored, including homotetramers, homopentamers, a homohexamer, and a heterohexamer.

The Molecular Basis for Life in Extreme Environments

Sampling and genomic efforts over the past decade have revealed an enormous quantity and diversity of life in Earth's extreme environments. This new knowledge of life on Earth poses the challenge of understandingits molecular basis in such inhospitable conditions, given that such conditions lead to loss of structure and of function in biomolecules from mesophiles. In this review, we discuss the physicochemical properties of extreme environments. We present the state of recent progress in extreme environmental genomics. We then present an overview of our current understanding of the biomolecular adaptation to extreme conditions. As our current and future understanding of biomolecular structure-function relationships in extremophiles requires methodologies adapted to extremes of pressure, temperature, and chemical composition, advances in instrumentation for probing biophysical properties under extreme conditions are presented. Finally, we briefly discuss possible future directions in extreme biophysics.

The road less traveled in protein folding: evidence for multiple pathways

Free Energy Landscape theory of Protein Folding, introduced over 20 years ago, implies that a protein has many paths to the folded conformation with the lowest free energy. Despite the knowledge in principle, it has been remarkably hard to detect such pathways. The lack of such observations is primarily due to the fact that no one experimental technique can detect many parts of the protein simultaneously with the time resolution necessary to see such differences in paths. However, recent technical developments and employment of multiple experimental probes and folding prompts have illuminated multiple folding pathways in a number of proteins that had all previously been described with a single path.

Sampling of the conformational landscape of small proteins with Monte Carlo methods

Computer simulation provides an increasingly realistic picture of large-scale conformational change of proteins, but investigations remain fundamentally constrained by the femtosecond timestep of molecular dynamics simulations. For this reason, many biologically interesting questions cannot be addressed using accessible state-of-the-art computational resources. Here, we report the development of an all-atom Monte Carlo approach that permits the modelling of the large-scale conformational change of proteins using standard off-the-shelf computational hardware and standard all-atom force fields. We demonstrate extensive thermodynamic characterization of the folding process of the α-helical Trp-cage, the Villin headpiece and the β-sheet WW-domain. We fully characterize the free energy landscape, transition states, energy barriers between different states, and the per-residue stability of individual amino acids over a wide temperature range. We demonstrate that a state-of-the-art intramolecular force field can be combined with an implicit solvent model to obtain a high quality of the folded structures and also discuss limitations that still remain.

N-Terminal Decarboxylation as a Probe for Intramolecular Contact Formation in γ-Glu-(Pro)n-Met Peptides

The kinetics of intramolecular-contact formation between remote functional groups in peptides with restricted conformational flexibility were examined using designed peptides with variable-length proline bridges. As probes for this motion, free radicals were produced using the ^•OH-induced oxidation at the C-terminal methionine residue of γ-Glu-(Pro)_n-Met peptides (n = 0-3). The progress of the radicals' motion along the proline bridges was monitored as the radicals underwent reactions along the peptides' backbones. Of particular interest was the reaction between the sulfur atom located in the side chain of the oxidized Met residue and the unprotonated amino group of the glutamic acid moiety. Interactions between them were probed by the radiation-chemical yields (expressed as G values) of the formation of C-centered, α-aminoalkyl radicals (αN) on the Glu residue. These radicals were monitored directly or via their reaction with p-nitroacetophenone (PNAP) to generate the optically detected PNAP^•- radical anions. The yields of these αN radicals were found to be linearly dependent on the number of Pro residues. A constant decrease by 0.09 μM J^-1 per spacing Pro residue of the radiation-chemical yields of G(αN) was observed. Previous reports support the conclusion that the αN radicals in these cases would have to result from (S∴N)⁺-bonded cyclic radical cations that arose as a result from direct contact between the ends of the peptides. Furthermore, by analogy with the rate constants for the formation of intramolecularly (S∴S)⁺-bonded radical cations in Met-(Pro)_n-Met peptides ( J. Phys. Chem. B 2016, 120, 9732), the rate constants for the formation of intramolecularly (S∴N)⁺-bonded radical cations are activated to the same extent for all of the γ-Glu-(Pro)_n-Met peptides. Thus, the continuous decrease of G(αN) with the number of Pro residues (from 0 to 3) suggests that the formation of a contact between the S-atom in the C-terminal Met residue and the N-atom of a deprotonated N-terminal amino group of Glu is controlled in peptides with 0 to 3 Pro residues by the relative diffusion of the S^•+ and unoxidized N-atom. The overall rate constants of cyclization to form the (S∴N)-bonded radical cations were estimated to be 3.8 × 10⁶, 1.8 × 10⁶, and 8.1 × 10⁵ s^-1 for peptides with n = 0, 1, and 2 Pro residues, respectively. If activation is the same for all of the peptides, then these rate constants are a direct indication for the end-to-end dynamics along the chain.

Studying biomolecular folding and binding using temperature-jump mass spectrometry

Characterizing folding and complex formation of biomolecules provides a view into their thermodynamics, kinetics and folding pathways. Deciphering kinetic intermediates is particularly important because they can often be targeted by drugs. The key advantage of native mass spectrometry over conventional methods that monitor a single observable is its ability to identify and quantify coexisting species. Here, we show the design of a temperature-jump electrospray source for mass spectrometry that allows one to perform fast kinetics experiments (0.16-32 s) at different temperatures (10-90 °C). The setup allows recording of both folding and unfolding kinetics by using temperature jumps from high to low, and low to high, temperatures. Six biological systems, ranging from peptides to proteins to DNA complexes, exemplify the use of this device. Using temperature-dependent experiments, the folding and unfolding of a DNA triplex are studied, providing detailed information on its thermodynamics and kinetics.

Competition of individual domain folding with inter-domain interaction in WW domain engineered repeat proteins

Engineered repeat proteins have proven to be a fertile ground for studying the competition between folding, misfolding and transient aggregation of tethered protein domains. We examine the interplay between folding and inter-domain interactions of engineered FiP35 WW domain repeat proteins with n = 1 through 5 repeats. We characterize protein expression, thermal and guanidium melts, as well as laser T-jump kinetics. All experimental data is fitted by a global fitting model with two states per domain (U, N), plus a third state M to account for non-native states due to domain interactions present in all but the monomer. A detailed structural model is provided by coarse-grained simulated annealing using the AWSEM Hamiltonian. Tethered FiP35 WW domains with n = 2 and 3 domains are just slightly less stable than the monomer. The n = 4 oligomer is yet less stable, its expression yield is much lower than the monomer's, and depends on the purification tag used. The n = 5 plasmid did not express at all, indicating the sudden onset of aggregation past n = 4. Thus, tethered FiP35 has a critical nucleus size for inter-domain aggregation of n ≈ 4. According to our simulations, misfolded structures become increasingly prevalent as one proceeds from monomer to pentamer, with extended inter-domain beta sheets appearing first, then multi-sheet 'intramolecular amyloid' structures, and finally novel motifs containing alpha helices. We discuss the implications of our results for oligomeric aggregate formation and structure, transient aggregation of proteins whilst folding, as well as for protein evolution that starts with repeat proteins.

Proteins: "Boil 'Em, Mash 'Em, Stick 'Em in a Stew"

Cells of the vast majority of organisms are subject to temperature, pressure, pH, ionic strength, and other stresses. We discuss these effects in the light of protein folding and protein interactions in vitro, in complex environments, in cells, and in vivo. Protein phase diagrams provide a way of organizing different structural ensembles that occur under stress and how one can move among ensembles. Experiments that perturb biomolecules in vitro or in cells by stressing them have revealed much about the underlying forces that are competing to control protein stability, folding, and function. Two phenomena that emerge and serve to broadly classify effects of the cellular environment are crowding (mainly due to repulsive forces) and sticking (mainly due to attractive forces). The interior of cells is closely balanced between these emergent effects, and stress can tip the balance one way or the other. The free energy scale involved is small but significant on the scale of the "on/off switches" that control signaling in cells or of protein-protein association with a favorable function such as increased enzyme processivity. Quantitative tools from biophysical chemistry will play an important role in elucidating the world of crowding and sticking under stress.

Fast pressure-jump all-atom simulations and experiments reveal site-specific protein dehydration-folding dynamics

As theory and experiment have shown, protein dehydration is a major contributor to protein folding. Dehydration upon folding can be characterized directly by all-atom simulations of fast pressure drops, which create desolvated pockets inside the nascent hydrophobic core. Here, we study pressure-drop refolding of three λ-repressor fragment (λ_6-85) mutants computationally and experimentally. The three mutants report on tertiary structure formation via different fluorescent helix-helix contact pairs. All-atom simulations of pressure drops capture refolding and unfolding of all three mutants by a similar mechanism, thus validating the nonperturbative nature of the fluorescent contact probes. Analysis of simulated interprobe distances shows that the α-helix 1-3 pair distance displays a slower characteristic time scale than the 1-2 or 3-2 pair distance. To see whether slow packing of α-helices 1 and 3 is reflected in the rate-limiting folding step, fast pressure-drop relaxation experiments captured refolding on a millisecond time scale. These experiments reveal that refolding monitored by 1-3 contact formation indeed is much slower than when monitored by 1-2 or 3-2 contact formation. Unlike the case of the two-state folder [three-α-helix bundle (α₃D)], whose drying and core formation proceed in concert, λ_6-85 repeatedly dries and rewets different local tertiary contacts before finally forming a solvent-excluded core, explaining the non-two-state behavior observed during refolding in molecular dynamics simulations. This work demonstrates that proteins can explore desolvated pockets and dry globular states numerous times before reaching the native conformation.

A quantitative connection of experimental and simulated folding landscapes by vibrational spectroscopy

We break the barrier between simulation and experiment by comparing identical computed and experimental infrared observables.

For small molecule reaction kinetics, computed reaction coordinates often mimic experimentally measured observables quite accurately. Although nowadays simulated and measured biomolecule kinetics can be compared on the same time scale, a gap between computed and experimental observables remains. Here we directly compared temperature-jump experiments and molecular dynamics simulations of protein folding dynamics using the same observable: the time-dependent infrared spectrum. We first measured the stability and folding kinetics of the fastest-folding β-protein, the GTT35 WW domain, using its structurally specific infrared spectrum. The relaxation dynamics of the peptide backbone, β-sheets, turn, and random coil were measured independently by probing the amide I′ region at different frequencies. Next, the amide I′ spectra along folding/unfolding molecular dynamics trajectories were simulated by accurate mixed quantum/classical calculations. The simulated time dependence and spectral amplitudes at the exact experimental probe frequencies provided relaxation and folding rates in agreement with experimental observations. The calculations validated by experiment yield direct structural evidence for a rate-limiting reaction step where an intermediate state with either the first or second hairpin is formed. We show how folding switches from a more homogeneous (apparent two-state) process at high temperature to a more heterogeneous process at low temperature, where different parts of the WW domain fold at different rates.

Practical aspects of high-pressure NMR spectroscopy and its applications in protein biophysics and structural biology

Pressure and temperature are the two fundamental variables of thermodynamics. Temperature and chemical perturbation are central experimental tools for the exploration of macromolecular structure and dynamics. Though it has long been recognized that hydrostatic pressure offers a complementary and often unique view of macromolecular structure, stability and dynamics, it has not been employed nearly as much. For solution NMR applications the limited use of high-pressure is undoubtedly traced to difficulties of employing pressure in the context of modern multinuclear and multidimensional NMR. Limitations in pressure tolerant NMR sample cells have been overcome and enable detailed studies of macromolecular energy landscapes, hydration, dynamics and function. Here we review the practical considerations for studies of biological macromolecules at elevated pressure, with a particular emphasis on applications in protein biophysics and structural biology.

Protein Folding Free Energy Landscape along the Committor - the Optimal Folding Coordinate

Recent advances in simulation and experiment have led to dramatic increases in the quantity and complexity of produced data, which makes the development of automated analysis tools very important. A powerful approach to analyze dynamics contained in such data sets is to describe/approximate it by diffusion on a free energy landscape - free energy as a function of reaction coordinates (RC). For the description to be quantitatively accurate, RCs should be chosen in an optimal way. Recent theoretical results show that such an optimal RC exists; however, determining it for practical systems is a very difficult unsolved problem. Here we describe a solution to this problem. We describe an adaptive nonparametric approach to accurately determine the optimal RC (the committor) for an equilibrium trajectory of a realistic system. In contrast to alternative approaches, which require a functional form with many parameters to approximate an RC and thus extensive expertise with the system, the suggested approach is nonparametric and can approximate any RC with high accuracy without system specific information. To avoid overfitting for a realistically sampled system, the approach performs RC optimization in an adaptive manner by focusing optimization on less optimized spatiotemporal regions of the RC. The power of the approach is illustrated on a long equilibrium atomistic folding simulation of HP35 protein. We have determined the optimal folding RC - the committor, which was confirmed by passing a stringent committor validation test. It allowed us to determine a first quantitatively accurate protein folding free energy landscape. We have confirmed the recent theoretical results that diffusion on such a free energy profile can be used to compute exactly the equilibrium flux, the mean first passage times, and the mean transition path times between any two points on the profile. We have shown that the mean squared displacement along the optimal RC grows linear with time as for simple diffusion. The free energy profile allowed us to obtain a direct rigorous estimate of the pre-exponential factor for the folding dynamics.

Automated method to differentiate between native and mirror protein models obtained from contact maps

Mirror protein structures are often considered as artifacts in modeling protein structures. However, they may soon become a new branch of biochemistry. Moreover, methods of protein structure reconstruction, based on their residue-residue contact maps, need methodology to differentiate between models of native and mirror orientation, especially regarding the reconstructed backbones. We analyzed 130 500 structural protein models obtained from contact maps of 1 305 SCOP domains belonging to all 7 structural classes. On average, the same numbers of native and mirror models were obtained among 100 models generated for each domain. Since their structural features are often not sufficient for differentiating between the two types of model orientations, we proposed to apply various energy terms (ETs) from PyRosetta to separate native and mirror models. To automate the procedure for differentiating these models, the k-means clustering algorithm was applied. Using total energy did not allow to obtain appropriate clusters–the accuracy of the clustering for class A (all helices) was no more than 0.52. Therefore, we tested a series of different k-means clusterings based on various combinations of ETs. Finally, applying two most differentiating ETs for each class allowed to obtain satisfying results. To unify the method for differentiating between native and mirror models, independent of their structural class, the two best ETs for each class were considered. Finally, the k-means clustering algorithm used three common ETs: probability of amino acid assuming certain values of dihedral angles Φ and Ψ, Ramachandran preferences and Coulomb interactions. The accuracies of clustering with these ETs were in the range between 0.68 and 0.76, with sensitivity and selectivity in the range between 0.68 and 0.87, depending on the structural class. The method can be applied to all fully-automated tools for protein structure reconstruction based on contact maps, especially those analyzing big sets of models.

Study of protein folding under native conditions by rapidly switching the hydrostatic pressure inside an NMR sample cell

Development of specialized instrumentation enables rapid switching of the hydrostatic pressure inside an operating NMR spectrometer. This technology allows observation of protein signals during the repeated folding process. Applied to ubiquitin, a previously extensively studied model of protein folding, the methodology reveals an initially highly dynamic state that deviates relatively little from random coil behavior but also provides evidence for numerous repeatedly failed folding events, previously only observed in computer simulations. Above room temperature, direct NMR evidence shows a ∼50% fraction of proteins folding through an on-pathway kinetic intermediate, thereby revealing two equally efficient parallel folding pathways.

In general, small proteins rapidly fold on the timescale of milliseconds or less. For proteins with a substantial volume difference between the folded and unfolded states, their thermodynamic equilibrium can be altered by varying the hydrostatic pressure. Using a pressure-sensitized mutant of ubiquitin, we demonstrate that rapidly switching the pressure within an NMR sample cell enables study of the unfolded protein under native conditions and, vice versa, study of the native protein under denaturing conditions. This approach makes it possible to record 2D and 3D NMR spectra of the unfolded protein at atmospheric pressure, providing residue-specific information on the folding process. ¹⁵N and ¹³C chemical shifts measured immediately after dropping the pressure from 2.5 kbar (favoring unfolding) to 1 bar (native) are close to the random-coil chemical shifts observed for a large, disordered peptide fragment of the protein. However, ¹⁵N relaxation data show evidence for rapid exchange, on a ∼100-μs timescale, between the unfolded state and unstable, structured states that can be considered as failed folding events. The NMR data also provide direct evidence for parallel folding pathways, with approximately one-half of the protein molecules efficiently folding through an on-pathway kinetic intermediate, whereas the other half fold in a single step. At protein concentrations above ∼300 μM, oligomeric off-pathway intermediates compete with folding of the native state.

Extracting the Hidden Distributions Underlying the Mean Transition State Structures in Protein Folding

The inherent conflict between noncovalent interactions and the large conformational entropy of the polypeptide chain forces folding reactions and their mechanisms to deviate significantly from chemical reactions. Accordingly, measures of structure in the transition state ensemble (TSE) are strongly influenced by the underlying distributions of microscopic folding pathways that are challenging to discern experimentally. Here, we present a detailed analysis of 150,000 folding transition paths of five proteins at three different thermodynamic conditions from an experimentally consistent statistical mechanical model. We find that the underlying TSE structural distributions are rarely unimodal, and the average experimental measures arise from complex underlying distributions. Unfolding pathways also exhibit subtle differences from folding counterparts due to a combination of Hammond behavior and native-state movements. Local interactions and topological complexity, to a lesser extent, are found to determine pathway heterogeneity, underscoring the importance of the balance between local and nonlocal energetics in protein folding.

Protein unfolding mechanisms and their effects on folding experiments

In this review, I discuss the various methods researchers use to unfold proteins in the lab in order to understand protein folding both in vitro and in vivo. The four main techniques, chemical-, heat-, pressure- and force-denaturation, produce distinctly different unfolded conformational ensembles. Recent measurements have revealed different folding kinetics from different unfolding mechanisms. Thus, comparing these distinct unfolded ensembles sheds light on the underlying free energy landscape of folding.

Parallel folding pathways of Fip35 WW domain explained by infrared spectra and their computer simulation

We present a calculation of the amide I' infrared (IR) spectra of the folded, unfolded, and intermediate states of the WW domain Fip35, a model system for β-sheet folding. Using an all-atom molecular dynamics simulation in which multiple folding and unfolding events take place we identify six conformational states and then apply perturbed matrix method quantum-mechanical calculations to determine their amide I' IR spectra. Our analysis focuses on two states previously identified as Fip35 folding intermediates and suggests that a three-stranded core similar to the folded state core is the main source of the spectroscopic differences between the two intermediates. In particular, we propose a hypothesis for why folding via one of these intermediates was not experimentally observed by IR T-jump.

Globular Protein Folding In Vitro and In Vivo

In vitro, computational, and theoretical studies of protein folding have converged to paint a rich and complex energy landscape. This landscape is sensitively modulated by environmental conditions and subject to evolutionary pressure on protein function. Of these environments, none is more complex than the cell itself, where proteins function in the cytosol, in membranes, and in different compartments. A wide variety of kinetic and thermodynamics experiments, ranging from single-molecule studies to jump kinetics and from nuclear magnetic resonance to imaging on the microscope, have elucidated how protein energy landscapes facilitate folding and how they are subject to evolutionary constraints and environmental perturbation. Here we review some recent developments in the field and refer the reader to some original work and additional reviews that cover this broad topic in protein science.

Toward a quantitative description of microscopic pathway heterogeneity in protein folding

Experimentally consistent statistical modeling of protein folding thermodynamics reveals unprecedented complexity with numerous parallel folding routes in five different proteins.

Refined Parameterization of Nonbonded Interactions Improves Conformational Sampling and Kinetics of Protein Folding Simulations

Recent advances in computational technology have enabled brute-force molecular dynamics (MD) simulations of protein folding using physics-based molecular force fields. The extensive sampling of protein conformations afforded by such simulations revealed, however, considerable compaction of the protein conformations in the unfolded state, which is inconsistent with experiment. Here, we show that a set of surgical corrections to nonbonded interactions between amine nitrogen-carboxylate oxygen and aliphatic carbon-carbon atom pairs can considerably improve the realism of protein folding simulations. Specifically, we show that employing our corrections in ∼500 μs all-atom replica-exchange MD simulations of the WW domain and villin head piece proteins increases the size of the denatured proteins' conformations and does not destabilize the native conformations of the proteins. In addition to making the folded conformations a global minimum of the respective free energy landscapes at room temperature, our corrections also make the free energy landscape smoother, considerably accelerating the folding kinetics and, hence, reducing the computational expense of a protein folding simulation.

A chimeric protein-based malaria vaccine candidate induces robust T cell responses against Plasmodium vivax MSP119

The most widespread Plasmodium species, Plasmodium vivax, poses a significant public health threat. An effective vaccine is needed to reduce global malaria burden. Of the erythrocytic stage vaccine candidates, the 19 kDa fragment of the P. vivax Merozoite Surface Protein 1 (PvMSP1₁₉) is one of the most promising. Our group has previously defined several promiscuous T helper epitopes within the PvMSP1 protein, with features that allow them to bind multiple MHC class II alleles. We describe here a P. vivax recombinant modular chimera based on MSP1 (PvRMC-MSP1) that includes defined T cell epitopes genetically fused to PvMSP1₁₉. This vaccine candidate preserved structural elements of the native PvMSP1₁₉ and elicited cytophilic antibody responses, and CD4+ and CD8+ T cells capable of recognizing PvMSP1₁₉. Although CD8+ T cells that recognize blood stage antigens have been reported to control blood infection, CD8+ T cell responses induced by P. falciparum or P. vivax vaccine candidates based on MSP1₁₉ have not been reported. To our knowledge, this is the first time a protein based subunit vaccine has been able to induce CD8+ T cell against PvMSP1₁₉. The PvRMC-MSP1 protein was also recognized by naturally acquired antibodies from individuals living in malaria endemic areas with an antibody profile associated with protection from infection. These features make PvRMC-MSP1 a promising vaccine candidate.

Activation Barrier-Limited Folding and Conformational Sampling of a Dynamic Protein Domain

Folding reaction mechanisms of globular protein domains have been extensively studied by both experiment and simulation and found to be highly concerted chemical reactions in which numerous noncovalent bonds form in an apparent two-state fashion. However, less is known regarding intrinsically disordered proteins because their folding can usually be studied only in conjunction with binding to a ligand. We have investigated by kinetics the folding mechanism of such a disordered protein domain, the nuclear coactivator-binding domain (NCBD) from CREB-binding protein. While a previous computational study suggested that NCBD folds without an activation free energy barrier, our experimental data demonstrate that NCBD, despite its highly dynamic structure, displays relatively slow folding (∼10 ms at 277 K) consistent with a barrier-limited process. Furthermore, the folding kinetics corroborate previous nuclear magnetic resonance data showing that NCBD exists in two folded conformations and one more denatured conformation at equilibrium and, thus, that the folding mechanism is a three-state mechanism. The refolding kinetics is limited by unfolding of the less populated folded conformation, suggesting that the major route for interconversion between the two folded states is via the denatured state. Because the two folded conformations have been suggested to bind distinct ligands, our results have mechanistic implications for conformational sampling in protein-protein interactions.

Energetic and topological determinants of a phosphorylation-induced disorder-to-order protein conformational switch

We show that the phosphorylation of 4E-BP2 acts as a triggering event to shape its folding-function landscape that is delicately balanced between conflicting favorable energetics and intrinsically unfavorable topological connectivity. We further provide first evidence that the fitness landscapes of proteins at the threshold of disorder can differ considerably from ordered domains.

Quantifying Protein Disorder through Measures of Excess Conformational Entropy

Intrinsically disordered proteins (IDPs) and proteins with a large degree of disorder are abundant in the proteomes of eukaryotes and viruses, and play a vital role in cellular homeostasis and disease. One fundamental question that has been raised on IDPs is the process by which they offset the entropic penalty involved in transitioning from a heterogeneous ensemble of conformations to a much smaller collection of binding-competent states. However, this has been a difficult problem to address, as the effective entropic cost of fixing residues in a folded-like conformation from disordered amino acid neighborhoods is itself not known. Moreover, there are several examples where the sequence complexity of disordered regions is as high as well-folded regions. Disorder in such cases therefore arises from excess conformational entropy determined entirely by correlated sequence effects, an entropic code that is yet to be identified. Here, we explore these issues by exploiting the order-disorder transitions of a helix in Pbx-Homeodomain together with a dual entropy statistical mechanical model to estimate the magnitude and sign of the excess conformational entropy of residues in disordered regions. We find that a mere 2.1-fold increase in the number of allowed conformations per residue (∼0.7kBT favoring the unfolded state) relative to a well-folded sequence, or ∼2(N) additional conformations for a N-residue sequence, is sufficient to promote disorder under physiological conditions. We show that this estimate is quite robust and helps in rationalizing the thermodynamic signatures of disordered regions in important regulatory proteins, modeling the conformational folding-binding landscapes of IDPs, quantifying the stability effects characteristic of disordered protein loops and their subtle roles in determining the partitioning of folding flux in ordered domains. In effect, the dual entropy model we propose provides a statistical thermodynamic basis for the relative conformational propensities of amino acids in folded and disordered environments in proteins. Our work thus lays the foundation for understanding and quantifying protein disorder through measures of excess conformational entropy.

The Effect of Fluorescent Protein Tags on Phosphoglycerate Kinase Stability Is Nonadditive

It is frequently assumed that fluorescent protein tags used in biological imaging experiments are minimally perturbing to their host protein. As in-cell experiments become more quantitative and measure rates and equilibrium constants, rather than just "on-off" activity or the presence of a protein, it becomes more important to understand such perturbations. One criterion for a protein modification to be a perturbation is additivity of two perturbations (a linear effect on the protein free energy). Here we show that adding fluorescent protein tags to a host protein in vitro has a large nonadditive effect on its folding free energy. We compare an unlabeled, three singly labeled, and a doubly labeled enzyme (phosphoglycerate kinase). We propose two mechanisms for nonadditivity. In the "quinary interaction" mechanism, two tags interact transiently with one another, relieving the host protein from unfavorable tag-protein interactions. In the "crowding" mechanism, adding two tags provides the minimal crowding necessary to overcome destabilizing interactions of individual tags with the host protein. Both of these mechanisms affect protein stability in cells; we show here that they must also be considered for tagged proteins used for reference in vitro.

High-Resolution Mapping of the Folding Transition State of a WW Domain

Fast-folding WW domains are among the best-characterized systems for comparing experiments and simulations of protein folding. Recent microsecond-resolution experiments and long duration (totaling milliseconds) single-trajectory modeling have shown that even mechanistic changes in folding kinetics due to mutation can now be analyzed. Thus, a comprehensive set of experimental data would be helpful to benchmark the predictions made by simulations. Here, we use T-jump relaxation in conjunction with protein engineering and report mutational Φ-values (Φ(M)) as indicators for folding transition-state structure of 65 side chain, 7 backbone hydrogen bond, and 6 deletion and /or insertion mutants within loop 1 of the 34-residue hPin1 WW domain. Forty-five cross-validated consensus mutants could be identified that provide structural constraints for transition-state structure within all substructures of the WW domain fold (hydrophobic core, loop 1, loop 2, β-sheet). We probe the robustness of the two hydrophobic clusters in the folding transition state, discuss how local backbone disorder in the native-state can lead to non-classical Φ(M)-values (Φ(M) > 1) in the rate-determining loop 1 substructure, and conclusively identify mutations and positions along the sequence that perturb the folding mechanism from loop 1-limited toward loop 2-limited folding.

The Role of Electrostatic Interactions in Folding of β-Proteins

Atomic-level molecular dynamic simulations are capable of fully folding structurally diverse proteins; however, they are limited in their ability to accurately represent electrostatic interactions. Here we have experimentally tested the role of charged residues on stability and folding kinetics of one of the most widely simulated β-proteins, the WW domain. The folding of wild type Pin1 WW domain, which has two positively charged residues in the first turn, was compared to the fast folding mutant FiP35 Pin1, which introduces a negative charge into the first turn. A combination of FTIR spectroscopy and laser-induced temperature-jump coupled with infrared spectroscopy was used to probe changes in the amide I region. The relaxation dynamics of the peptide backbone, β-sheets and β-turns, and negatively charged aspartic acid side chain of FiP35 were measured independently by probing the corresponding bands assigned in the amide I region. Folding is initiated in the turns and the β-sheets form last. While the global folding mechanism is in good agreement with simulation predictions, we observe changes in the protonation state of aspartic acid during folding that have not been captured by simulation methods. The protonation state of aspartic acid is coupled to protein folding; the apparent pKa of aspartic acid in the folded protein is 6.4. The dynamics of the aspartic acid follow the dynamics of the intermediate phase, supporting assignment of this phase to formation of the first hairpin. These results demonstrate the importance of electrostatic interactions in turn stability and formation of extended β-sheet structures.

Fast-folding proteins under stress

Proteins are subject to a variety of stresses in biological organisms, including pressure and temperature, which are the easiest stresses to simulate by molecular dynamics. We discuss the effect of pressure and thermal stress on very-fast-folding model proteins, whose in vitro folding can be fully simulated on computers and compared with experiments. We then discuss experiments that can be used to subject proteins to low- and high-temperature unfolding, as well as low- and high-pressure unfolding. Pressure and temperature are prototypical perturbations that illustrate how close many proteins are to instability, a property that cells can exploit to control protein function. We conclude by reviewing some recent in-cell experiments, and progress being made in simulating and measuring protein stability and function inside live cells.