Thermodynamics and folding kinetics analysis of the SH3 domain form discrete molecular dynamics

Oligomer Formation by Physiologically Relevant C-Terminal Isoforms of Amyloid β-Protein

Alzheimer's disease (AD) is a neurological disorder associated with amyloid β-protein (Aβ) assembly into toxic oligomers. In addition to the two predominant alloforms, Aβ1-40 and Aβ1-42, other C-terminally truncated Aβ peptides, including Aβ1-38 and Aβ1-43, are produced in the brain. Here, we use discrete molecular dynamics (DMD) and a four-bead protein model with amino acid-specific hydropathic interactions, DMD4B-HYDRA, to examine oligomer formation of Aβ1-38, Aβ1-40, Aβ1-42, and Aβ1-43. Self-assembly of 32 unstructured monomer peptides into oligomers is examined using 32 replica DMD trajectories for each of the four peptides. In a quasi-steady state, Aβ1-38 and Aβ1-40 adopt similar unimodal oligomer size distributions with a maximum at trimers, whereas Aβ1-42 and Aβ1-43 oligomer size distributions are multimodal with the dominant maximum at trimers or tetramers, and additional maxima at hexamers and unidecamers (for Aβ1-42) or octamers and pentadecamers (for Aβ1-43). The free energy landscapes reveal isoform- and oligomer-order specific structural and morphological features of oligomer ensembles. Our results show that oligomers of each of the four isoforms have unique features, with Aβ1-42 alone resulting in oligomers with disordered and solvent-exposed N-termini. Our findings help unravel the structure-function paradigm governing oligomers formed by various Aβ isoforms.

A singularity-free torsion angle potential for coarse-grained molecular dynamics simulations

Conventional torsion angle potentials used in molecular dynamics (MD) have a singularity problem when three bonded particles are collinearly aligned. This problem is often encountered in coarse-grained (CG) simulations. Here, we propose a new form of the torsion angle potential, which introduces an angle-dependent modulating function. By carefully tuning the parameters for this modulating function, our method can eliminate the problematic angle-dependent singularity while being combined with existing models. As an example, we optimized the modulating function of the torsion angle potential for popular CG models of biomolecules based on the statistics over experimental structures deposited in the Protein Data Bank. By applying our method to designed and natural biomolecules, we show that the new torsion angle potential is able to eliminate the singularity problem while maintaining the structural features in the original models. Furthermore, by comparing our design with previous methods, we found that our new potential has advantages in computational efficiency and numerical stability. We strongly recommend the usage of our new potential in the CG simulations of flexible molecules.

Amyloid Self-Assembly of hIAPP8-20 via the Accumulation of Helical Oligomers, α-Helix to β-Sheet Transition, and Formation of β-Barrel Intermediates

The self-assembly of human islet amyloid polypeptide (hIAPP) into β-sheet-rich nanofibrils is associated with the pathogeny of type 2 diabetes. Soluble hIAPP is intrinsically disordered with N-terminal residues 8-17 as α-helices. To understand the contribution of the N-terminal helix to the aggregation of full-length hIAPP, here the oligomerization dynamics of the hIAPP fragment 8-20 (hIAPP8-20) are investigated with combined computational and experimental approaches. hIAPP8-20 forms cross-β nanofibrils in silico from isolated helical monomers via the helical oligomers and α-helices to β-sheets transition, as confirmed by transmission electron microscopy, atomic force microscopy, circular dichroism spectroscopy, Fourier transform infrared spectroscopy, and reversed-phase high performance liquid chromatography. The computational results also suggest that the critical nucleus of aggregation corresponds to hexamers, consistent with a recent mass-spectroscopy study of hIAPP8-20 aggregation. hIAPP8-20 oligomers smaller than hexamers are helical and unstable, while the α-to-β transition starts from the hexamers. Converted β-sheet-rich oligomers first form β-barrel structures as intermediates before aggregating into cross-β nanofibrils. This study uncovers a complete picture of hIAPP8-20 peptide oligomerization, aggregation nucleation via conformational conversion, formation of β-barrel intermediates, and assembly of cross-β protofibrils, thereby shedding light on the aggregation of full-length hIAPP, a hallmark of pancreatic beta-cell degeneration.

Microcanonical thermostatistics of coarse-grained proteins with amyloidogenic propensity

The formation of fibrillar aggregates seems to be a common characteristic of polypeptide chains, although the observation of these aggregates may depend on appropriate experimental conditions. Partially folded intermediates seem to have an important role in the generation of protein aggregates, and a mechanism for this fibril formation considers that these intermediates also correspond to metastable states with respect to the fibrillar ones. Here, using a coarse-grained (CG) off-lattice model, we carry out a comparative analysis of the thermodynamic aspects characterizing the folding transition with respect to the propensity for aggregation of four different systems: two isoforms of the amyloid β-protein, the Src SH3 domain, and the human prion proteins (hPrP). Microcanonical analysis of the data obtained from replica exchange method is conducted to evaluate the free-energy barrier and latent heat in these models. The simulations of the amyloid β isoforms and Src SH3 domain indicated that the folding process described by this CG model is related to a negative specific heat, a phenomenon that can only be verified in the microcanonical ensemble in first-order phase transitions. The CG simulation of the hPrP heteropolymer yielded a continuous folding transition. The absence of a free-energy barrier and latent heat favors the presence of partially unfolded conformations, and in this context, this thermodynamic aspect could explain the reason why the hPrP heteropolymer is more aggregation-prone than the other heteropolymers considered in this study. We introduced the hydrophobic radius of gyration as an order parameter and found that it can be used to obtain reliable information about the hydrophobic packing and the transition temperatures in the folding process.

An effective solvent theory connecting the underlying mechanisms of osmolytes and denaturants for protein stability

An all-atom Gō model of Trp-cage protein is simulated using discontinuous molecular dynamics in an explicit minimal solvent, using a single, contact-based interaction energy between protein and solvent particles. An effective denaturant or osmolyte solution can be constructed by making the interaction energy attractive or repulsive. A statistical mechanical equivalence is demonstrated between this effective solvent model and models in which proteins are immersed in solutions consisting of water and osmolytes or denaturants. Analysis of these studies yields the following conclusions: 1), Osmolytes impart extra stability to the protein by reducing the entropy of the unfolded state. 2), Unfolded states in the presence of osmolyte are more collapsed than in water. 3), The folding transition in osmolyte solutions tends to be less cooperative than in water, as determined by the ratio of van 't Hoff to calorimetric enthalpy changes. The decrease in cooperativity arises from an increase in native structure in the unfolded state, and thus a lower thermodynamic barrier at the transition midpoint. 4), Weak denaturants were observed to destabilize small proteins not by lowering the unfolded enthalpy, but primarily by swelling the unfolded state and raising its entropy. However, adding a strong denaturant destabilizes proteins enthalpically. 5), The folding transition in denaturant-containing solutions is more cooperative than in water. 6), Transfer to a concentrated osmolyte solution with purely hard-sphere steric repulsion significantly stabilizes the protein, due to excluded volume interactions not present in the canonical Tanford transfer model. 7), Although a solution with hard-sphere interactions adds a solvation barrier to native contacts, the folding is nevertheless less cooperative for reasons 1-3 above, because a hard-sphere solvent acts as a protecting osmolyte.

Turbulent phenomena in protein folding

Protein folding and hydrodynamic turbulence are two long-standing challenges, in molecular biophysics and fluid dynamics, respectively. The theories of these phenomena have been developed independently and used different formalisms. Here we show that the protein folding flows can be surprisingly similar to turbulent fluid flows. Studying a benchmark model protein (an SH3 domain), we have found that the flows for the slow folding trajectories of the protein, in which a partly formed N- and C-terminal β sheet hinders the RT loop from attaching to the protein core, have many properties of turbulent flows of a fluid. The flows are analyzed in a three-dimensional (3D) space of collective variables, which are the numbers of native contacts between the terminal β strands, between the RT loop and the protein core, and the rest of the native contacts. We have found that the flows have fractal nature and are filled with 3D eddies; the latter contain strange attractors, at which the tracer flow paths behave as saddle trajectories. Two regions of the space increment have been observed, in which the flux variations are self-similar with the scaling exponent h=1/3, in surprising agreement with the Kolmogorov inertial range theory of turbulence. In one region, the cascade of protein rearrangements is directed from larger to smaller scales (net folding), and in the other, it is oppositely directed (net unfolding). Folding flows for the fast trajectories are essentially "laminar" and do not have the property of self-similarity. Based on the results of our study, we infer, and support this inference by simulations, that the origin of the similarity between the protein folding and turbulent motion of a fluid is in a cascade mechanism of structural transformations in the systems that underlies these phenomena.

Nucleation process in the folding of a domain-swapped dimer

Nucleation processes are important for the understanding in protein dynamics. To evaluate the effect of nucleation mechanism in dimerization process, a domain-swapped dimer (Esp8) is simulated with the symmetrized Gō model and the classical Gō model. The pathways of the dimerization are analyzed with computational phi -analysis method. It is found out that some nuclei are observed in the kinetic steps of the dimeric association though the whole pathway is a process with multiple intermediate states. The key residues in the nuclei are rather similar to those observed in the monomeric folding. The differences with the monomeric cases are also discussed. These differences illustrate the effects of dimeric feature on the nucleation process. Besides, manual mutations are carried out to illustrate the importance of the interactions related to the nuclei. It is observed that the mutations in the nuclei-related interactions apparently change the dynamics while other mutations have little effect on the kinetics. All of these results outline a picture that the nucleation processes act as the fundamental steps of high-order organization of protein systems.

Analysis of the free-energy surface of proteins from reversible folding simulations

Computer generated trajectories can, in principle, reveal the folding pathways of a protein at atomic resolution and possibly suggest general and simple rules for predicting the folded structure of a given sequence. While such reversible folding trajectories can only be determined ab initio using all-atom transferable force-fields for a few small proteins, they can be determined for a large number of proteins using coarse-grained and structure-based force-fields, in which a known folded structure is by construction the absolute energy and free-energy minimum. Here we use a model of the fast folding helical λ-repressor protein to generate trajectories in which native and non-native states are in equilibrium and transitions are accurately sampled. Yet, representation of the free-energy surface, which underlies the thermodynamic and dynamic properties of the protein model, from such a trajectory remains a challenge. Projections over one or a small number of arbitrarily chosen progress variables often hide the most important features of such surfaces. The results unequivocally show that an unprojected representation of the free-energy surface provides important and unbiased information and allows a simple and meaningful description of many-dimensional, heterogeneous trajectories, providing new insight into the possible mechanisms of fast-folding proteins.

The process of protein folding is a complex transition from a disordered to an ordered state. Here, we simulate a specific fast-folding protein at the point at which the native and denatured states are at equilibrium and show that obtaining an accurate description of the mechanisms of folding and unfolding is far from trivial. Using simple quantities which quantify the degree of native order is, in the case of this protein, clearly misleading. We show that an unbiased representation of the free-energy surface can be obtained; using such a representation we are able to redesign the landscape and thus modify, upon site-specific “mutations”, the folding and unfolding rates. This leads us to formulate a hypothesis to explain the very fast folding of many proteins.

Multiscale approaches for studying energy transduction in dynein

Cytoplasmic dynein is an important motor that drives all minus-end directed movement along microtubules. Dynein is a complex motor whose processive motion is driven by ATP-hydrolysis. Dynein's run length has been measured to be several millimetres with typical velocities in the order of a few nanometres per second. Therefore, the average time between steps is a fraction of a second. When this time scale is compared with typical time scales for protein side chain and backbone movements (approximately 10(-9) s and approximately 10(-5) s, respectively), it becomes clear that a multi-timescale modelling approach is required to understand energy transduction in this protein. Here, we review recent efforts to use computational and mathematical modelling to understand various aspects of dynein's chemomechanical cycle. First, we describe a structural model of dynein's motor unit showing a heptameric organization of the motor subunits. Second, we describe our molecular dynamics simulations of the motor unit that are used to investigate the dynamics of the various motor domains. Third, we present a kinetic model of the coordination between the two dynein heads. Lastly, we investigate the various potential geometries of the dimer during its hydrolytic and stepping cycle.

Satisfaction of hydrogen-bonding potential influences the conservation of polar sidechains

Although polar amino acids tend to be found on the surface of proteins due to their hydrophilic nature, their important roles within the core of proteins are now becoming better recognized. It has long been understood that a significant number of mainchain functions will not achieve hydrogen bond satisfaction through the formation of secondary structures; in these circumstances, it is generally buried polar residues that provide hydrogen bond satisfaction. Here, we describe an analysis of the hydrogen-bonding of polar amino acids in a set of structurally aligned protein families. This allows us not only to calculate the conservation of each polar residue but also to assess whether conservation is correlated with the hydrogen-bonding potential of polar sidechains. We show that those polar sidechains whose hydrogen-bonding potential is satisfied tend to be more conserved than their unsatisfied or nonhydrogen-bonded counterparts, particularly when buried. Interestingly, these buried and satisfied polar residues are significantly more conserved than buried hydrophobic residues. Forming hydrogen bonds to mainchain amide atoms also influences conservation, with those satisfied buried polar residues that form two hydrogen bonds to mainchain amides being significantly more conserved than those that form only one or none. These results indicate that buried polar residues whose hydrogen-bonding potential is satisfied are likely to have important roles in maintaining protein structure.

Calculation of protein heat capacity from replica-exchange molecular dynamics simulations with different implicit solvent models

The heat capacity has played a major role in relating microscopic and macroscopic properties of proteins and their disorder-order phase transition of folding. Its calculation by atomistic simulation methods remains a significant challenge due to the complex and dynamic nature of protein structures, their solvent environment, and configurational averaging. To better understand these factors on calculating a protein heat capacity, we provide a comparative analysis of simulation models that differ in their implicit solvent description and force-field resolution. Our model protein system is the src Homology 3 (SH3) domain of alpha-spectrin, and we report a series of 10 ns replica-exchange molecular dynamics simulations performed at temperatures ranging from 298 to 550 K, starting from the SH3 native structure. We apply the all-atom CHARMM22 force field with different modified analytical generalized Born solvent models (GBSW and GBMV2) and compare these simulation models with the distance-dependent dielectric screening of charge-charge interactions. A further comparison is provided with the united-atom CHARMM19 plus a pairwise GB model. Unfolding-folding transition temperatures of SH3 were estimated from the temperature-dependent profiles of the heat capacity, root-mean-square distance from the native structure, and the fraction of native contacts, each calculated from the density of states by using the weighted histogram analysis method. We observed that, for CHARMM22, the unfolding transition and energy probability density were quite sensitive to the implicit solvent description, in particular, the treatment of the protein-solvent dielectric boundary in GB models and their surface-area-based hydrophobic term. Among the solvent models tested, the calculated melting temperature varied in the range 353-438 K and was higher than the experimental value near 340 K. A reformulated GBMV2 model of employing a smoother molecular-volume dielectric interface was the most accurate in reproducing the native conformation and a two-state folding landscape, although the melting transition temperature did not show the smallest deviation from experiment. For the lower-resolution CHARMM19/GB model, the simulations failed to yield a bimodal energy distribution, yet the melting temperature was observed to be a good estimate of higher-resolution simulation models. We also demonstrate that a careful analysis of a relatively long simulation is necessary to avoid trapping in local minima and to find a true thermodynamic transition temperature.

Simple continuous and discrete models for simulating replica exchange simulations of protein folding

The efficiency of temperature replica exchange (RE) simulations hinge on their ability to enhance conformational sampling at physiological temperatures by taking advantage of more rapid conformational interconversions at higher temperatures. While temperature RE is a parallel simulation technique that is relatively straightforward to implement, kinetics in the RE ensemble is complicated, and there is much to learn about how best to employ RE simulations in computational biophysics. Protein folding rates often slow down above a certain temperature due to entropic bottlenecks. This "anti-Arrhenius" behavior represents a challenge for RE. However, it is far from straightforward to systematically explore the impact of this on RE by brute force molecular simulations, since RE simulations of protein folding are very difficult to converge. To understand some of the basic mechanisms that determine the efficiency of RE, it is useful to study simplified low dimensionality systems that share some of the key characteristics of molecular systems. Results are presented concerning the efficiency of temperature RE on a continuous two-dimensional potential that contains an entropic bottleneck. Optimal efficiency was obtained when the temperatures of the replicas did not exceed the temperature at which the harmonic mean of the folding and unfolding rates is maximized. This confirms a result we previously obtained using a discrete network model of RE. Comparison of the efficiencies obtained using the continuous and discrete models makes it possible to identify non-Markovian effects, which slow down equilibration of the RE ensemble on the more complex continuous potential. In particular, the rate of temperature diffusion and also the efficiency of RE is limited by the time scale of conformational rearrangements within free energy basins.

Probing protein aggregation using discrete molecular dynamics

Understanding the role of biomolecular dynamics in cellular processes leading to human diseases and the ability to rationally manipulate these processes is of fundamental importance in scientific research. The last decade has witnessed significant progress in probing biophysical behavior of proteins. However, we are still limited in understanding how changes in protein dynamics and inter-protein interactions occurring in short length- and time-scales lead to aberrations in their biological function. Bridging this gap in biology probed using computer simulations marks a challenging frontier in computational biology. Here we examine hypothesis-driven simplified protein models in conjunction with discrete molecular dynamics in the study of protein aggregation, implicated in series of neurodegenerative diseases, such as Alzheimer's and Huntington's diseases. Discrete molecular dynamics simulations of simplified protein models have emerged as a powerful methodology with its ability to bridge the gap in time and length scales from protein dynamics to aggregation, and provide an indispensable tool for probing protein aggregation.

Coarse-grained lattice model simulations of sequence-structure fitness of a ribosome-inactivating protein

Many realistic protein-engineering design problems extend beyond the computational limits of what is considered practical when applying all-atom molecular-dynamics simulation methods. Lattice models provide computationally robust alternatives, yet most are regarded as too simplistic to accurately capture the details of complex designs. We revisit a coarse-grained lattice simulation model and demonstrate that a multiresolution modeling approach of reconstructing all-atom structures from lattice chains is of sufficient accuracy to resolve the comparability of sequence-structure modifications of the ricin A-chain (RTA) protein fold. For a modeled structure, the unfolding-folding transition temperature was calculated from the heat capacity using either the potential energy from the lattice model or the all-atom CHARMM19 force-field plus a generalized Born solvent approximation. We found, that despite the low-resolution modeling of conformational states, the potential energy functions were capable of detecting the relative change in the thermodynamic transition temperature that distinguishes between a protein design and the native RTA fold in excellent accord with reported experimental studies of thermal denaturation. A discussion is provided of different sequences fitted to the RTA fold and a possible unfolding model.

Protein folding: then and now

Over the past three decades the protein folding field has undergone monumental changes. Originally a purely academic question, how a protein folds has now become vital in understanding diseases and our abilities to rationally manipulate cellular life by engineering protein folding pathways. We review and contrast past and recent developments in the protein folding field. Specifically, we discuss the progress in our understanding of protein folding thermodynamics and kinetics, the properties of evasive intermediates, and unfolded states. We also discuss how some abnormalities in protein folding lead to protein aggregation and human diseases.

Parallel folding pathways in the SH3 domain protein

The transition-state ensemble (TSE) is the set of protein conformations with an equal probability to fold or unfold. Its characterization is crucial for an understanding of the folding process. We determined the TSE of the src-SH3 domain protein by using extensive molecular dynamics simulations of the Go model and computing the folding probability of a generated set of TSE candidate conformations. We found that the TSE possesses a well-defined hydrophobic core with variable enveloping structures resulting from the superposition of three parallel folding pathways. The most preferred pathway agrees with the experimentally determined TSE, while the two least preferred pathways differ significantly. The knowledge of the different pathways allows us to design the interactions between amino acids that guide the protein to fold through the least preferred pathway. This particular design is akin to a circular permutation of the protein. The finding motivates the hypothesis that the different experimentally observed TSEs in homologous proteins and circular permutants may represent potentially available pathways to the wild-type protein.

Folding behavior of ribosomal protein S6 studied by modified Gō-like model

Recent experimental and theoretical studies suggest that, although topology is the determinant factor in protein folding, especially for small single-domain proteins, energetic factors also play an important role in the folding process. The ribosomal protein S6 has been subjected to intensive studies. A radical change of the transition state in its circular permutants has been observed, which is believed to be caused by a biased distribution of contact energies. Since the simplistic topology-only Gō-like model is not able to reproduce such an observation, we modify the model by introducing variable contact energies between residues based on their physicochemical properties. The modified Gō-like model can successfully reproduce the Phi-value distributions, folding nucleus, and folding pathways of both the wild-type and circular permutants of S6. Furthermore, by comparing the results of the modified and the simplistic models, we find that the hydrophobic effect constructs the major force that balances the loop entropies. This may indicate that nature maintains the folding cooperativity of this protein by carefully arranging the location of hydrophobic residues in the sequence. Our study reveals a strategy or mechanism used by nature to get out of the dilemma when the native structure, possibly required by biological function, conflicts with folding cooperativity. Finally, the possible relationship between such a design of nature and amyloidosis is also discussed.

Sequence of events in folding mechanism: beyond the Gō model

Simplified Gō models, where only native contacts interact favorably, have proven useful to characterize some aspects of the folding of small proteins. The success of these models is limited by the fact that all residues interact in the same way so that the folding features of a protein are determined only by the geometry of its native conformation. We present an extended version of a Calpha-based Gō model where different residues interact with different energies. The model is used to calculate the thermodynamics of three small proteins (Protein G, Src-SH3, and CI2) and the effect of mutations (DeltaDeltaGU-N, DeltaDeltaGdouble dagger-N, DeltaDeltaGdouble dagger-U, and phi-values) on the wild-type sequence. The model allows us to investigate some of the most controversial areas in protein folding, such as its earliest stages and the nature of the unfolded state, subjects that have lately received particular attention.

Simple but predictive protein models

The traditional approach to computational biophysics studies of molecular systems is brute force molecular dynamics simulations under the conditions of interest. The disadvantages of this approach are that the time and length scales that are accessible to computer simulations often do not reach biologically relevant scales. An alternative approach, which we call intuitive modeling, is hypothesis-driven and based on tailoring simplified protein models to the systems of interest. Using intuitive modeling, the length and time scales that can be achieved using simplified protein models exceed those of traditional molecular-dynamic simulations. Here, we describe several recent studies that signify the predictive power of simplified protein models within the intuitive-modeling approach.

A Single Disulfide Bond Differentiates Aggregation Pathways of ß2-Microglobulin

Deposition of wild-type beta2-microglobulin (beta2m) into amyloid fibrils is a complication in patients undergoing long-term hemodialysis. The native beta-sandwich fold of beta2m has a highly conserved disulfide bond linking Cys25 and Cys80. Oxidized beta2m forms needle-like amyloid fibrils at pH 2.5 in vitro, whereas reduced beta2m, at acid pH, in which the intra-chain disulfide bond is disrupted, cannot form typical fibrils. Instead, reduced beta2m forms thinner and more flexible filaments. To uncover the difference in molecular mechanisms underlying the aggregation of the oxidized and reduced beta2m, we performed molecular dynamics simulations of beta2m oligomerization under oxidized and reduced conditions. We show that, consistent with experimental observations, the oxidized beta2m forms domain-swapped dimer, in which the two proteins exchange their N-terminal segments complementing each other. In contrast, both dimers and trimers, formed by reduced beta2m, are comprised of parallel beta-sheets between monomers and stabilized by the hydrogen bond network along the backbone. The oligomerized monomers are in extended conformations, capable of further aggregation. We find that both reduced and oxidized dimers are thermodynamically less stable than their corresponding monomers, indicating that beta2m oligomerization is not accompanied by the formation of a thermodynamically stable dimer. Our studies suggest that the different aggregation pathways of oxidized and reduced beta2m are dictated by the formation of distinct precursor oligomeric species that are modulated by Cys25-Cys80 disulfide-bonds. We propose that the propagation of domain swapping is the aggregation mechanism for the oxidized beta2m, while "parallel stacking" of partially unfolded beta2m is the aggregation mechanism for the reduced beta2m.

Effects of disulfide bonds on folding behavior and mechanism of the beta-sheet protein tendamistat

Tendamistat, a small disulfide-bonded beta-sheet protein, and its three single/double-disulfide mutants are investigated by using a modified Gō-like model, aiming to understand the folding mechanism of disulfide-bonded protein as well as the effects of removal of disulfide bond on the folding process. Our simulations show that tendamistat and its two single-disulfide mutants are all two-state folders, consistent with the experimental observations. It is found that the disulfide bonds as well as three hydrogen bonds between the N-terminal loop-0 and strand-6 are of significant importance for the folding of tendamistat. Without these interactions, their two-state behaviors become unstable and the predictions of the model are inconsistent with experiments. In addition, the effect of disulfide bonds on the folding process are studied by comparing the wild-type tendamistat and its two mutants; it is found that the removal of either of the C11-C27 or C45-C73 disulfide bond leads to a large decrease in the thermodynamical stability and loss of structure in the unfolded state, and the effect of the former is stronger than that of the later. These simulation results are in good agreement with experiments and, thus, validate our model. Based on the same model, the detailed folding pathways of the wild-type tendamistat and two mutants are studied, and the effect of disulfide bonds on the folding kinetics are discussed. The obtained results provide a detailed folding picture of these proteins and complement experimental findings. Finally, the folding nuclei predicted to be existent in this protein tendamistat as well as its mutants are firstly identified in this work. The positions of the nucleus are consistent with those argued in experimental studies. Therefore, a nucleation/growth folding mechanism that can explain the two-state folding manner is clearly characterized. Moreover, the effect by the removal of each disulfide bond on the folding thermodynamics and dynamics can also be well interpreted from their influence on the folding nucleus. The implementation of this work indicates that the modified Gō-like model really describes the folding behavior of protein tendamistat and could be used to study the folding of other disulfide-bonded proteins.

Structure and stability of a model three-helix-bundle protein on tailored surfaces

The interaction of protein molecules with surfaces is important in numerous applications. Theoretical work on protein adsorption has been limited. In particular, it is difficult to obtain quantitative predictions about the structure and stability of proteins on surfaces. In this study, density-of-states-based simulations were performed on a Gō-like model of a three-helix-bundle fragment from protein A (PDB ID: 1bdd). Both mechanical and thermal stability were investigated on neutral and attractive surfaces and compared to that in the absence of a surface. It was found that attaching the peptide to any type of surface decreases its melting temperature by as much as 9 K, depending upon orientation. Calorimetric cooperativity, as measured by van't Hoff to calorimetric enthalpy ratios, similarly decreased. It was also found that the mechanical strength of the peptide attached to surfaces is degraded to varying extents, depending upon the surface type and protein orientation. A comparison of mechanical and thermal stability showed that the two are not synonymous, but occur through different pathways, and that system configurations that are more thermally stable are not always so mechanically.

Acquisition of Fyn-Selective SH3 Domain Ligands via a Combinatorial Library Strategy

A stepwise library-based strategy has been employed to acquire a potent ligand for the SH3 domain of Fyn, a Src kinase family member that plays a key role in T cell activation. The easily automated methodology is designed to identify potential interaction sites that circumscribe the protein/peptide binding region on the SH3 domain. The library protocol creates peptide/nonpeptide chimeras that are able to bind to these interaction sites that are otherwise inaccessible to natural amino acid residues. The peptide-derived lead and the Fyn-SH3 domain form a complex that exhibits a K(D) of 25 +/- 5 nM, approximately 1000-fold more potent than that displayed by the corresponding conventional peptide ligand. Furthermore, the lead ligand exhibits selectivity against SH3 domains derived from other Src kinases, in spite of a sequence identity of approximately 80%.

How well can simulation predict protein folding kinetics and thermodynamics?

Simulation of protein folding has come a long way in five years. Notably, new quantitative comparisons with experiments for small, rapidly folding proteins have become possible. As the only way to validate simulation methodology, this achievement marks a significant advance. Here, we detail these recent achievements and ask whether simulations have indeed rendered quantitative predictions in several areas, including protein folding kinetics, thermodynamics, and physics-based methods for structure prediction. We conclude by looking to the future of such comparisons between simulations and experiments.

Scaling Behavior and Structure of Denatured Proteins

An ensemble of random-coil conformations with no persistent structures has long been accepted as the classical model of denatured proteins due to its consistency with the experimentally determined scaling of protein sizes. However, recent NMR spectroscopy studies on proteins at high chemical denaturant concentrations suggest the presence of significant amounts of native-like structures, in contrast to the classical random-coil picture. To reconcile these seemingly controversial observations, we examine thermally denatured states of experimentally characterized proteins by using molecular dynamics simulations. For all studied proteins, we find that denatured states indeed have strong local conformational bias toward native states while a random-coil power law scaling of protein sizes is preserved. In addition, we explain why experimentally determined size of the protein creatine kinase does not follow general scaling. In simulations, we observe that this protein exhibits a stable intermediate state, the size of which is consistent with the reported experimental observation.

Folding events in the 21-30 region of amyloid beta-protein (Abeta) studied in silico

Oligomeric assemblies of the amyloid beta-protein (Abeta) have been implicated in the pathogenesis of Alzheimer's disease as a primary source of neurotoxicity. Recent in vitro studies have suggested that a 10-residue segment, Ala-21-Ala-30, forms a turn-like structure that nucleates the folding of the full-length Abeta. To gain a mechanistic insight, we simulated Abeta(21-30) folding by using a discrete molecular dynamics algorithm and a united-atom model incorporating implicit solvent and a variable electrostatic interaction strength (EIS). We found that Abeta(21-30) folds into a loop-like conformation driven by an effective hydrophobic attraction between Val-24 and the butyl portion of the Lys-28 side chain. At medium EIS [1.5 kcal/mol (1 cal = 4.18 J)], unfolded conformations almost disappear, in agreement with experimental observations. Under optimal conditions for folding, Glu-22 and Asp-23 form transient electrostatic interactions (EI) with Lys-28 that stabilize the loop conformations. Glu-22-Lys-28 is the most favored interaction. High EIS, as it occurs in the interior of proteins and aggregates, destabilizes the packing of Val-24 and Lys-28. Analysis of the unpacked structures reveals strong EI with predominance of the Asp-23-Lys-28 interaction, in agreement with studies of molecular modeling of full-length Abeta fibrils. The binary nature of the EI involving Lys-28 provides a mechanistic explanation for the linkage of amino acid substitutions at Glu-22 with Alzheimer's disease and cerebral amyloid angiopathy. Substitutions may alter the frequency of Glu-22 or Asp-23 involvement in contact formation and affect the stability of the folding nucleus formed in the Abeta(21-30) region.

Nucleation and the transition state of the SH3 domain

We present a verified computational model of the SH3 domain transition state (TS) ensemble. This model was built for three separate SH3 domains using experimental phi-values as structural constraints in all-atom protein folding simulations. While averaging over all conformations incorrectly considers non-TS conformations as transition states, quantifying structures as pre-TS, TS, and post-TS by measurement of their transmission coefficient ("probability to fold", or p(fold)) allows for rigorous conclusions regarding the structure of the folding nucleus and a full mechanistic analysis of the folding process. Through analysis of the TS, we observe a highly polarized nucleus in which many residues are solvent-exposed. Mechanistic analysis suggests the hydrophobic core forms largely after an early nucleation step. SH3 presents an ideal system for studying the nucleation-condensation mechanism and highlights the synergistic relationship between experiment and simulation in the study of protein folding.

Molecular dynamics simulation of amyloid beta dimer formation

Recent experiments with amyloid beta (Abeta) peptide indicate that formation of toxic oligomers may be an important contribution to the onset of Alzheimer's disease. The toxicity of Abeta oligomers depends on their structure, which is governed by assembly dynamics. Due to limitations of current experimental techniques, a detailed knowledge of oligomer structure at the atomic level is missing. We introduce a molecular dynamics approach to study Abeta dimer formation. 1), We use discrete molecular dynamics simulations of a coarse-grained model to identify a variety of dimer conformations; and 2), we employ all-atom molecular mechanics simulations to estimate thermodynamic stability of all dimer conformations. Our simulations of a coarse-grained Abeta peptide model predicts 10 different planar beta-strand dimer conformations. We then estimate the free energies of all dimer conformations in all-atom molecular mechanics simulations with explicit water. We compare the free energies of Abeta(1-42) and Abeta(1-40) dimers. We find that 1), dimer conformations have higher free energies compared to their corresponding monomeric states; and 2), the free-energy difference between the Abeta(1-42) and the corresponding Abeta(1-40) dimer conformation is not significant. Our results suggest that Abeta oligomerization is not accompanied by the formation of thermodynamically stable planar beta-strand dimers.

Energy landscapes for cooperative processes: nearly ideal glass transitions, liquid-liquid transitions and folding transitions

We describe basic phenomenology in the physics of supercooling liquids at constant volume (most simulations), and at constant pressure (most laboratory experiments) before focusing attention on the exceptional cases that exhibit liquid-liquid phase transitions on constant-pressure cooling. We give evidence for point defects in glasses and liquids near T(g). Models based on defects predict transitions with density gaps in constant-pressure systems. We describe the energy landscape representation of such systems. Water, in these terms, is post-critical, and its nearly ideal glass formation can be related to nucleation-free protein 'funnel-folding'. For nucleated folding of proteins, a pseudo-gap should be present. Experimental methods of distinguishing between alternative folding scenarios are described.

Multiple folding pathways of the SH3 domain

Experimental observations suggest that proteins follow different folding pathways under different environmental conditions. We perform molecular dynamics simulations of a model of the c-Crk SH3 domain over a broad range of temperatures, and identify distinct pathways in the folding transition. We determine the kinetic partition temperature-the temperature for which the c-Crk SH3 domain undergoes a rapid folding transition with minimal kinetic barriers-and observe that below this temperature the model protein may undergo a folding transition by multiple folding pathways via only one or two intermediates. Our findings suggest the hypothesis that the SH3 domain, a protein fold for which only two-state folding kinetics was observed in previous experiments, may exhibit intermediate states under conditions that strongly stabilize the native state.

Statistical mechanics of protein folding by cluster distance geometry

This is our second type of model for protein folding where the configurational parameters and the effective potential energy function are chosen in such a way that all conformations are described and the canonical partition function can be evaluated analytically. Structure is described in terms of distances between pairs of sequentially contiguous blocks of eight residues, and all possible conformations are grouped into 71 subsets in terms of bounds on these distances. The energy is taken to be a sum of pairwise interactions between such blocks. The 210 energy parameters were adjusted so that the native folds of 32 small proteins are favored in free energy over the denatured state. We then found 146 proteins having negligible sequence similarity to any of the training proteins, yet the free energy of the respective correct native states were favored over the denatured state.

In silico study of amyloid beta-protein folding and oligomerization

Experimental findings suggest that oligomeric forms of the amyloid beta protein (Abeta) play a critical role in Alzheimer's disease. Thus, elucidating their structure and the mechanisms of their formation is critical for developing therapeutic agents. We use discrete molecular dynamics simulations and a four-bead protein model to study oligomerization of two predominant alloforms, Abeta40 and Abeta42, at the atomic level. The four-bead model incorporates backbone hydrogen-bond interactions and amino acid-specific interactions mediated through hydrophobic and hydrophilic elements of the side chains. During the simulations we observe monomer folding and aggregation of monomers into oligomers of variable sizes. Abeta40 forms significantly more dimers than Abeta42, whereas pentamers are significantly more abundant in Abeta42 relative to Abeta40. Structure analysis reveals a turn centered at Gly-37-Gly-38 that is present in a folded Abeta42 monomer but not in a folded Abeta40 monomer and is associated with the first contacts that form during monomer folding. Our results suggest that this turn plays an important role in Abeta42 pentamer formation. Abeta pentamers have a globular structure comprising hydrophobic residues within the pentamer's core and hydrophilic N-terminal residues at the surface of the pentamer. The N termini of Abeta40 pentamers are more spatially restricted than Abeta42 pentamers. Abeta40 pentamers form a beta-strand structure involving Ala-2-Phe-4, which is absent in Abeta42 pentamers. These structural differences imply a different degree of hydrophobic core exposure between pentamers of the two alloforms, with the hydrophobic core of the Abeta42 pentamer being more exposed and thus more prone to form larger oligomers.

New Insights into FAK Signaling and Localization Based on Detection of a FAT Domain Folding Intermediate

Mounting evidence suggests that the focal adhesion targeting (FAT) domain, an antiparallel four-helix bundle, exists in alternative conformations that may modulate phosphorylation, ligand binding, and the subcellular localization of focal adhesion kinase (FAK). In order to characterize the conformational dynamics of the FAT domain, we have developed a novel method for reconstructing the folding pathway of the FAT domain by using discrete molecular dynamics (DMD) simulations, with free energy constraints derived from NMR hydrogen exchange data. The DMD simulations detect a folding intermediate, in which a cooperative unfolding event causes helix 1 to lose helical character while separating from the helix bundle. The conformational dynamic features of helix 1 in the intermediate state of the FAT domain are likely to facilitate Y926 phosphorylation, yet interfere with paxillin binding. The presence of this intermediate state in vivo may promote FAK signaling via the ERK/MAPK pathway and by release of FAK from focal adhesions.

Formation of the folding nucleus of an SH3 domain investigated by loosely coupled molecular dynamics simulations

The experimentally well-established folding mechanism of the src-SH3 domain, and in particular the phi-value analysis of its transition state, represents a sort of testing table for computational investigations of protein folding. Here, parallel molecular dynamics simulations of the src-SH3 domain have been performed starting from denatured conformations. By rescuing and restarting only trajectories approaching the folding transition state, an ensemble of conformations was obtained with a completely structured central beta-sheet and a native-like packing of residues Ile-110, Ala-121, and Ile-132. An analysis of the trajectories shows that there are several pathways leading to the formation of the central beta-sheet whereas its two hairpins form in a different but consistent way.

Temperature dependence of the free energy landscape of the src-SH3 protein domain

The temperature dependence of the free energy landscape of the src-SH3 protein domain is investigated through fully atomic simulations in explicit solvent. Simulations are performed above and below the folding transition temperature, enabling an analysis of both protein folding and unfolding. The transition state for folding and unfolding, identified from the free energy surfaces, is found to be very similar, with structure in the central hydrophobic sheet and little structure throughout the rest of the protein. This is a result of a polarized folding (unfolding) mechanism involving early formation (late loss) of the central hydrophobic sheet at the transition state. Unfolding simulations map qualitatively well onto low-temperature free energy surfaces but appear, however, to miss important features observed in folding simulations. In particular, details of the folding mechanism involving the opening and closing of the hydrophobic core are not captured by unfolding simulations performed under strongly denaturing conditions. In addition, free energy surfaces at high temperatures do not display a desolvation barrier found at lower temperatures, involving the expulsion of water molecules from the hydrophobic core.

Statistical mechanics of protein folding with separable energy functions

We have initiated an entirely new approach to statistical mechanical models of strongly interacting systems where the configurational parameters and the potential energy function are both constructed so that the canonical partition function can be evaluated analytically. For a simplified model of proteins consisting of a single, fairly short polypeptide chain without cross-links, we can adjust the energy parameters to favor the experimentally determined native state of seven proteins having diverse types of folds. Then 497 test proteins are predicted to have stable native folds, even though they are also structurally diverse, and 480 of them have no significant sequence similarity to any of the training proteins.

Roles of physical interactions in determining protein-folding mechanisms: molecular simulation of protein G and alpha spectrin SH3

Protein-folding mechanisms of two small globular proteins, IgG binding domain of protein G and alpha spectrin SH3 domain are investigated via Brownian dynamics simulations with a model made of coarse-grained physical energy functions responsible for sequence-specific interactions and weak Gō-like energies. The folding pathways of alpha spectrin SH3 are known to be mainly controlled by the native topology, while protein G folding is anticipated to be more sensitive to the sequence-specific effects than native topology. We found in the folding of protein G that the C terminal beta hairpin is formed earlier and is rigid, once ordered, in the presence of an intact C terminal turn. The alpha helix is found to exhibit repeated partial formations/deformations during folding and to be stabilized via the tertiary contact with preformed beta sheets. This predicted scenario is fully consistent with experimental phi value data. Moreover, we found that the folding route is critically affected when the hydrophobic interaction is excluded from physical energy terms, suggesting that the hydrophobicity critically contributes to the folding propensity of protein G. For the folding of alpha spectrin SH3, we found that the distal beta hairpin and diverging turn are parts formed early, fully in harmony with previous results of simple Gō-like and experimental analysis, supporting that the folding route of SH3 domain is robust and coded by the native topology. The hybrid method provides useful tools for analyzing roles of physical interactions in determining folding mechanisms.

Discrete molecular dynamics simulations of peptide aggregation

We study the aggregation of peptides using the discrete molecular dynamics simulations. Specifically, at temperatures above the alpha-helix melting temperature of a single peptide, the model peptides aggregate into a multilayer parallel beta-sheet structure. This structure has an interstrand distance of 4.8 A and an intersheet distance of 10 A, which agree with experimental observations. Our model explains these results as follows: hydrogen-bond interactions give rise to the interstrand spacing in beta sheets, while Gō interactions between side chains make beta strands parallel to each other and allow beta sheets to pack into layers. An important feature of our results is that the aggregates contain free edges, which may allow for further aggregation of model peptides to form elongated fibrils.

Mechanism for the alpha-helix to beta-hairpin transition

The aggregation of alpha-helix-rich proteins into beta-sheet-rich amyloid fibrils is associated with fatal diseases, such as Alzheimer's disease and prion disease. During an aggregation process, protein secondary structure elements-alpha-helices-undergo conformational changes to beta-sheets. The fact that proteins with different sequences and structures undergo a similar transition on aggregation suggests that the sequence nonspecific hydrogen bond interaction among protein backbones is an important factor. We perform molecular dynamics simulations of a polyalanine model, which is an alpha-helix in its native state and observe a metastable beta-hairpin intermediate. Although a beta-hairpin has larger potential energy than an alpha-helix, the entropy of a beta-hairpin is larger because of fewer constraints imposed by the hydrogen bonds. In the vicinity of the transition temperature, we observe the interconversion of the alpha-helix and beta-sheet states via a random coil state. We also study the effect of the environment by varying the relative strength of side-chain interactions for a designed peptide-an alpha-helix in its native state. For a certain range of side-chain interaction strengths, we find that the intermediate beta-hairpin state is destabilized and even disappears, suggesting an important role of the environment in the aggregation propensity of a peptide.

A binding event converted into a folding event

We have designed a chimeric protein by connecting a circular permutant of the alpha-spectrin SH3 domain to the proline-rich decapeptide APSYSPPPPP with a three-residue link. Our aim was to obtain a single-chain protein with a tertiary fold that would mimic the binding between SH3 domains and proline-rich peptides. A comparison of the circular-dichroism and fluorescence spectra of the purified chimera and the SH3 circular permutant showed that the proline-rich sequence occupies the putative SH3 binding site in a similar conformation and with comparable interactions to those found in complexes between SH3 and proline-rich peptides. Differential scanning calorimetry indicated that the interactions in the binding motif interface are highly cooperative with the rest of the structure and thus the protein unfolds in a two-state process. The chimera is more stable than the circular permutant SH3 by 6-8 kJ mol(-1) at 25 degrees C and the difference in their unfolding enthalpy is approximately 32 kJ mol(-1), which coincides with the values found for the binding of proline-rich peptides to SH3 domains. This type of chimeric protein may be useful in designing SH3 peptide ligands with improved affinity and specificity.

Folding of Cu, Zn Superoxide Dismutase and Familial Amyotrophic Lateral Sclerosis

Cu, Zn superoxide dismutase (SOD1) has been implicated in the familial form of the neurodegenerative disease amyotrophic lateral sclerosis (ALS). It has been suggested that mutant mediated SOD1 misfolding/aggregation is an integral part of the pathology of ALS. We study the folding thermodynamics and kinetics of SOD1 using a hybrid molecular dynamics approach. We reproduce the experimentally observed SOD1 folding thermodynamics and find that the residues which contribute the most to SOD1 thermal stability are also crucial for apparent two-state folding kinetics. Surprisingly, we find that these residues are located on the surface of the protein and not in the hydrophobic core. Mutations in some of the identified residues are found in patients with the disease. We argue that the identified residues may play an important role in aggregation. To further characterize the folding of SOD1, we study the role of cysteine residues in folding and find that non-native disulfide bond formation may significantly alter SOD1 folding dynamics and aggregation propensity.

Posttransition state desolvation of the hydrophobic core of the src-SH3 protein domain

The folding thermodynamics of the src-SH3 protein domain were characterized under refolding conditions through biased fully atomic molecular dynamics simulations with explicit solvent. The calculated free energy surfaces along several reaction coordinates revealed two barriers. The first, larger barrier was identified as the transition state barrier for folding, associated with the formation of the first hydrophobic sheet of the protein. phi values calculated from structures residing at the transition state barrier agree well with experimental phi values. The microscopic information obtained from our simulations allowed us to unambiguously assign intermediate phi values as the result of multiple folding pathways. The second, smaller barrier occurs later in the folding process and is associated with the cooperative expulsion of water molecules between the hydrophobic sheets of the protein. This posttransition state desolvation barrier cannot be observed through traditional folding experiments, but is found to be critical to the correct packing of the hydrophobic core in the final stages of folding. Hydrogen exchange and NMR experiments are suggested to probe this barrier.

Molecular dynamics simulation of the SH3 domain aggregation suggests a generic amyloidogenesis mechanism

We use molecular dynamics simulation to study the aggregation of Src SH3 domain proteins. For the case of two proteins, we observe two possible aggregation conformations: the closed form dimer and the open aggregation state. The closed dimer is formed by "domain swapping"-the two proteins exchange their RT-loops. All the hydrophobic residues are buried inside the dimer so proteins cannot further aggregate into elongated amyloid fibrils. We find that the open structure-stabilized by backbone hydrogen bond interactions-packs the RT-loops together by swapping the two strands of the RT-loop. The packed RT-loops form a beta-sheet structure and expose the backbone to promote further aggregation. We also simulate more than two proteins, and find that the aggregate adopts a fibrillar double beta-sheet structure, which is formed by packing the RT-loops from different proteins. Our simulations are consistent with a possible generic amyloidogenesis scenario.