Molecular simulation of protein dynamics in nanopores. I. Stability and folding

Discontinuous molecular dynamics simulations, together with the protein intermediate resolution model, an intermediate-resolution model of proteins, are used to carry out several microsecond-long simulations and study folding transition and stability of alpha-de novo-designed proteins in slit nanopores. Both attractive and repulsive interaction potentials between the proteins and the pore walls are considered. Near the folding temperature T(f) and in the presence of the attractive potential, the proteins undergo a repeating sequence of folding/partially folding/unfolding transitions, with T(f) decreasing with decreasing pore sizes. The unfolded states may even be completely adsorbed on the pore's walls with a negative potential energy. In such pores the energetic effects dominate the entropic effects. As a result, the unfolded state is stabilized, with a folding temperature T(f) which is lower than its value in the bulk and that, compared with the bulk, the folding rate decreases. The opposite is true in the presence of a repulsive interaction potential between the proteins and the walls. Moreover, for short proteins in very tight pores with attractive walls, there exists an unfolded state with only one alpha-helical hydrogen bond and an energy nearly equal to that of the folded state. The proteins have, however, high entropies, implying that they cannot fold onto their native structure, whereas in the presence of repulsive walls the proteins do attain their native structure. There is a pronounced asymmetry between the two termini of the protein with respect to their interaction with the pore walls. The effect of a variety of factors, including the pore size and the proteins' length, as well as the temperature, is studied in detail.     

A localized specific interaction alters the unfolding pathways of structural homologues

Reductive unfolding studies of proteins are designed to provide information about intramolecular interactions that govern the formation (and stabilization) of the native state and about folding/unfolding pathways. By mutating Tyr92 to G, A, or L in the model protein, bovine pancreatic ribonuclease A, and through analysis of temperature factors and molecular dynamics simulations of the crystal structures of these mutants, it is demonstrated that the markedly different reductive unfolding rates and pathways of ribonuclease A and its structural homologue onconase can be attributed to a single, localized, ring-stacking interaction between Tyr92 and Pro93 in the bovine variant. The fortuitous location of this specific stabilizing interaction in a disulfide-bond-containing loop region of ribonuclease A results in the localized modulation of protein dynamics that, in turn, enhances the susceptibility of the disulfide bond to reduction leading to an alteration in the reductive unfolding behavior of the homologues. These results have important implications for folding studies involving topological determinants to obtain folding/unfolding rates and pathways, for protein structure-function prediction through fold recognition, and for predicting proteolytic cleavage sites.     

Influence of temperature, friction, and random forces on folding of the B-domain of staphylococcal protein A: all-atom molecular dynamics in implicit solvent

The influences of temperature, friction, and random forces on the folding of protein A have been analyzed. A series of all-atom molecular dynamics folding simulations with the Amber ff99 potential and Generalized Born solvation, starting from the fully extended chain, were carried out for temperatures from 300 to 500 K, using (a) the Berendsen thermostat (with no explicit friction or random forces) and (b) Langevin dynamics (with friction and stochastic forces explicitly present in the system). The simulation temperature influences the relative time scale of the major events on the folding pathways of protein A. At lower temperatures, helix 2 folds significantly later than helices 1 and 3. However, with increasing temperature, the folding time of helix 2 approaches the folding times of helices 1 and 3. At lower temperatures, the complete formation of secondary and tertiary structure is significantly separated in time whereas, at higher temperatures, they occur simultaneously. These results suggest that some earlier experimental and theoretical observations of folding events, e.g., the order of helix formation, could depend on the temperature used in those studies. Therefore, the differences in temperature used could be one of the reasons for the discrepancies among published experimental and computational studies of the folding of protein A. Friction and random forces do not change the folding pathway that was observed in the simulations with the Berendsen thermostat, but their explicit presence in the system extends the folding time of protein A.     

Model for the nucleation mechanism of protein folding

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

Monte Carlo vs molecular dynamics for all-atom polypeptide folding simulations

An efficient Monte Carlo (MC) algorithm including concerted rotations is directly compared to molecular dynamics (MD) in all-atom statistical mechanics folding simulations of small polypeptides. The previously reported algorithm "concerted rotations with flexible bond angles" (CRA) has been shown to successfully locate the native state of small polypeptides. In this study, the folding of three small polypeptides (trpzip2/H1/Trp-cage) is investigated using MC and MD, for a combined sampling time of approximately 10(11) MC configurations and 8 micros, respectively. Both methods successfully locate the experimentally determined native states of the three systems, but they do so at different speed, with 2-2.5 times faster folding of the MC runs. The comparison reveals that thermodynamic and dynamic properties can reliably be obtained by both and that results from folding simulations do not depend on the algorithm used. Similar to previous comparisons of MC and MD, it is found that one MD integration step of 2 fs corresponds to one MC scan, revealing the good sampling of MC. The simplicity and efficiency of the MC method will enable its future use in folding studies involving larger systems and the combination with replica exchange algorithms.     

A semi-analytical description of protein folding that incorporates detailed geometrical information

Much has been done to study the interplay between geometric and energetic effects on the protein folding energy landscape. Numerical techniques such as molecular dynamics simulations are able to maintain a precise geometrical representation of the protein. Analytical approaches, however, often focus on the energetic aspects of folding, including geometrical information only in an average way. Here, we investigate a semi-analytical expression of folding that explicitly includes geometrical effects. We consider a Hamiltonian corresponding to a Gaussian filament with structure-based interactions. The model captures local features of protein folding often averaged over by mean-field theories, for example, loop contact formation and excluded volume. We explore the thermodynamics and folding mechanisms of beta-hairpin and alpha-helical structures as functions of temperature and Q, the fraction of native contacts formed. Excluded volume is shown to be an important component of a protein Hamiltonian, since it both dominates the cooperativity of the folding transition and alters folding mechanisms. Understanding geometrical effects in analytical formulae will help illuminate the consequences of the approximations required for the study of larger proteins.     

Hierarchy of folding and unfolding events of protein G, CI2, and ACBP from explicit-solvent simulations

The study of the mechanism which is at the basis of the phenomenon of protein folding requires the knowledge of multiple folding trajectories under biological conditions. Using a biasing molecular-dynamics algorithm based on the physics of the ratchet-and-pawl system, we carry out all-atom, explicit solvent simulations of the sequence of folding events which proteins G, CI2, and ACBP undergo in evolving from the denatured to the folded state. Starting from highly disordered conformations, the algorithm allows the proteins to reach, at the price of a modest computational effort, nativelike conformations, within a root mean square deviation (RMSD) of approximately 1 A?. A scheme is developed to extract, from the myriad of events, information concerning the sequence of native contact formation and of their eventual correlation. Such an analysis indicates that all the studied proteins fold hierarchically, through pathways which, although not deterministic, are well-defined with respect to the order of contact formation. The algorithm also allows one to study unfolding, a process which looks, to a large extent, like the reverse of the major folding pathway. This is also true in situations in which many pathways contribute to the folding process, like in the case of protein G.     

Evidence of an intermediate and parallel pathways in protein unfolding from single-molecule fluorescence

Determining how proteins fold into their native structures is a subject of great importance, since ultimately it will allow protein structure and function to be predicted from primary sequence data. In addition, there is now a clear link between protein unfolding and misfolding events and many disease states. However, since proteins fold over rugged, multidimensional energy landscapes, this is a challenging experimental and theoretical problem. Single-molecule fluorescence methods developed over the past decade have the potential to follow the unfolding/folding of individual molecules. Mapping out the landscape without ensemble averaging will enable the identification of intermediate states which may not be significantly populated, in addition to the presence of multiple pathways. To date, there have been only a limited number of single-molecule folding/unfolding studies under nonequilibrium conditions and no intermediates have been observed. Here, for the first time, we present a single-molecule study of the unfolding of a large autofluorescent protein, Citrine, a variant of green fluorescent protein. Single-molecule fluorescence techniques are used to directly detect an intermediate on the unfolding/folding pathway and the existence of parallel unfolding pathways. This work, and the novel methods used, shows that single-molecule fluorescence can now provide new, hitherto experimentally inaccessible, insights into the folding/unfolding of proteins.     

Thermodynamics of Go-type models for protein folding

Go-type potentials, based on the inter-residue contacts present in the native structure of a protein, are frequently used to predict dynamic and structural features of the folding pathways through computer simulations. However, the mathematical form used to define the model interactions includes several arbitrary choices, whose consequences are not usually analyzed. In this work, we use a simple off-lattice protein model and a parallel tempering Monte Carlo simulation technique to carry out such analysis, centered in the thermodynamic characteristics of the folding transition. We show how the definition of a native contact has a deep impact on the presence of simple or complex transitions, with or without thermodynamic intermediates. In addition, we have checked that the width of the attractive wells has a profound effect on the free-energy barrier between the folded and unfolded states, mainly through its influence on the entropy of the denatured state.     

A ternary nucleation model for the nucleation pathway of protein folding

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

Infrared study of the stability and folding kinetics of a 15-residue beta-hairpin

The thermal stability and folding kinetics of a 15-residue beta-hairpin (SESYINPDGTWTVTE) have been studied by using infrared (IR) spectroscopy coupled with laser-induced temperature-jump (T-jump) technique for rapid folding-unfolding initiation. An alternative method based on analyzing IR difference spectra was also introduced to obtain thermodynamic properties of beta-sheets, which complements the commonly used circular dichroism (CD) and fluorescence techniques. Equilibrium IR measurements indicate that the thermal unfolding of this beta-hairpin is fairly broad. However, it can be described by a two-state transition with a thermal melting temperature of approximately 29 degrees C. Time-resolved IR measurements following a T-jump, probed at 1634 cm(-1), indicate that the folding of this beta-hairpin follows first-order kinetics and is amazingly fast. At 300 K, the folding time is approximately 0.8 micros, which is only 2-3 times slower than that of alpha-helix formation. Additionally, the energetic barrier for folding is small (approximately 2 kcal mol(-1)). These results, in conjunction with results from other studies, support a view that the details of native contacts play a dominant role in the kinetics of beta-hairpin folding.     

A new method for rapid characterization of the folding pathways of multidisulfide-containing proteins

Oxidative folding is a composite process that consists of both the conformational folding to the native three-dimensional structure and the regeneration of the native disulfide bonds of a protein, frequently involving over 100 disulfide intermediate species. Understanding the oxidative folding pathways of a multiple-disulfide-containing protein is a very difficult task that often requires years of devoted research due to the high complexity of the process and the very similar features of the large number of intermediates. Here we developed a method for rapidly delineating the major features of the oxidative folding pathways of a protein. The method examines the temperature dependence of the oxidative folding rate of the protein in combination with reduction pulses. Reduction pulses expose the presence of structured intermediates along the pathways. The correlation between the regeneration rate at different temperatures and the stability of the structured intermediates reveals the role that the intermediates play in determining the pathway. The method was first tested with bovine pancreatic ribonuclease A whose folding pathways were defined earlier. Then, it was explored to discern some of the major features of the folding pathways of its homologue, frog Onconase. The results suggest that the stability of the three-dimensional structure of the native protein is a major determinant of the folding rate in oxidative folding.     

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

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

Excluded volume entropic effects on protein unfolding times and intermediary stability

The dynamics of protein folding result from both enthalpic and entropic contributions to the free energy. In this paper we focus on entropic volume exclusion effects. We carry out computer simulations using a model that allows us to independently change the size or biochemical properties of amino acid residues. To determine the importance of excluded volume effects, we investigate the effects of changing the size of side chains on the unfolding dynamics of a model four-helix bundle protein. In addition, we also investigate the effects of changing the thickness of the chain's backbone. This has relevance to the behavior of synthetic polymers where the size of the constituent units can be varied. We find that entropic excluded volume effects are crucially important for stabilizing the organized native state relative to the molten globule.     

Balance of enthalpy and entropy in depletion forces

Solutes added to solutions often dramatically impact molecular processes ranging from the suspension or precipitation of colloids to biomolecular associations and protein folding. Here we revisit the origins of the effective attractive interactions that emerge between and within macromolecules immersed in solutions containing cosolutes that are preferentially excluded from the macromolecular interfaces. Until recently, these depletion forces were considered to be entropic in nature, resulting primarily from the tendency to increase the space available to the cosolute. However, recent experimental evidence indicates the existence of additional, energetically-dominated mechanisms. In this review we follow the emerging characteristics of these different mechanisms. By compiling a set of available thermodynamic data for processes ranging from protein folding to protein–protein interactions, we show that excluded cosolutes can act through two distinct mechanisms that correlate to a large extent with their molecular properties. For many polymers at low to moderate concentrations the steric interactions and molecular crowding effects dominate, and the mechanism is entropic. To contrast, for many small excluded solutes, such as naturally occurring osmolytes, the mechanism is dominated by favorable enthalpy, whereas the entropic contribution is typically unfavorable. We review the available models for these thermodynamic mechanisms, and comment on the need for new models that would be able to explain the full range of observed depletion forces.     

Expanding the realm of ultrafast protein folding: gpW, a midsize natural single-domain with alpha+beta topology that folds downhill

All ultrafast folding proteins known to date are either very small in size (less than 45 residues), have an alpha-helix bundle topology, or have been artificially engineered. In fact, many of them share two or even all three features. Here we show that gpW, a natural 62-residue alpha+beta protein expected to fold slowly in a two-state fashion, folds in microseconds (i.e., from tau = 33 micros at 310 K to tau = 1.7 micros at 355 K). Thermodynamic analyses of gpW reveal probe dependent thermal denaturation, complex coupling between two denaturing agents, and differential scanning calorimetry (DSC) thermogram characteristic of folding over a negligible thermodynamic folding barrier. The free energy surface analysis of gpW folding kinetics also produces a marginal folding barrier of about thermal energy ( RT) at the denaturation midpoint. From these results we conclude that gpW folds in the downhill regime and is close to the global downhill limit. This protein seems to be poised toward downhill folding by a loosely packed hydrophobic core with low aromatic content, large stabilizing contributions from local interactions, and abundance of positive charges on the native surface. These special features, together with a complex functional role in bacteriophage lambda assembly, suggest that gpW has been engineered to fold downhill by natural selection.     

Dissecting entropic coiling and poor solvent effects in protein collapse

The early events in protein collapse and folding are guided by the protein's elasticity. The contributions of entropic coiling and poor solvent effects like hydrophobic forces to the elastic response of proteins are currently unknown. Using molecular simulations of stretched ubiquitin in comparison with models of proteins as entropic chains, we find a surprisingly high stiffness of the protein backbone, reflected by a persistence length of 1.2 nm, which is significantly reduced by hydrophobic forces acting between protein side chains to an apparent persistence length of 0.3-0.6 nm frequently observed in single-molecule stretching experiments. Thus, the poor solvent conditions of a protein in water lead to a protein compaction much beyond the coiling of an entropic chain and thereby allow a protein to appear softer than when using good solvents.     

Protein folding kinetics: barrier effects in chemical and thermal denaturation experiments

Recent experimental work on fast protein folding brings about an intriguing paradox. Microsecond-folding proteins are supposed to fold near or at the folding speed limit (downhill folding), but yet their folding behavior seems to comply with classical two-state analyses, which imply the crossing of high free energy barriers. However, close inspection of chemical and thermal denaturation kinetic experiments in fast-folding proteins reveals systematic deviations from two-state behavior. Using a simple one-dimensional free energy surface approach we find that such deviations are indeed diagnostic of marginal folding barriers. Furthermore, the quantitative analysis of available fast-kinetic data indicates that many microsecond-folding proteins fold downhill in native conditions. All of these proteins are then promising candidates for an atom-by-atom analysis of protein folding using nuclear magnetic resonance.1 We also find that the diffusion coefficient for protein folding is strongly temperature dependent, corresponding to an activation energy of approximately 1 kJ.mol-1 per protein residue. As a consequence, the folding speed limit at room temperature is about an order of magnitude slower than the approximately 1 micros estimates from high-temperature T-jump experiments. Our analysis is quantitatively consistent with the available thermodynamic and kinetic data on slow two-state folding proteins and provides a straightforward explanation for the apparent fast-folding paradox.     

Constraint dynamics algorithm for simulation of semiflexible macromolecules

Semiflexible models are often used to study macromolecules containing stable structural elements. Based on rigid body dynamics, we developed a rigid fragment constraint dynamics algorithm for the simulation of semiflexible macromolecules. Stable structural elements are treated as rigid fragments. Rigid fragment constraints, defined as combinations of distance constraints and position constraints, are introduced to limit internal molecular motion to the required mode. The constraint forces are solved separately for each rigid fragment constraint and iteratively until all constraint conditions are satisfied within a given tolerance at each time step, as is done for the bond length constraint in the SHAKE algorithm. The orientation of a rigid fragment is represented by the quaternion parameters, and both translation and rotation are solved by the leap-frog formulation. We tested the algorithm with molecular dynamics simulations of a series of peptides and a small protein. The computation cost for the constraints is roughly proportional to the size of the molecule. In the microcanonical ensemble simulation of polyvalines, the total energy was conserved satisfactorily with time steps as large as 20 fs. A helix folding simulation of a synthetic peptide was carried out to show the efficiency of the algorithm in a conformational search. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1555–1566, 1998     

Predictive in silico all-atom folding of a four-helix protein with a free-energy model

We report the predictive all-atom folding of the 60 amino acid four-helix bacterial ribosomal protein (BRP) L20 with a stochastic evolutionary optimization method in a free-energy force field. The energetically best, as well as six of the 10 lowest conformations, converge to a near-native structure. All of the 10 best energy conformations share the secondary structure elements of the native conformation, but differ in their tertiary alignment. The best conformation has a backbone root-mean-square deviation of 4.6 A to the native conformation and reproduces most distance constraints of the NMR experiment to 1.5 A resolution. Starting from random initial conditions, the native content of the simulated population increases more than 60-fold in the course of the simulation. These data demonstrate the feasibility of predictive unbiased all-atom protein folding with present day computational resources for the BRP L20.