Molecular simulation of polymer assisted protein refolding

Protein refolding in vitro, the formation of the tertiary structure that enables the protein to display its biological function, can be significantly enhanced by adding a polymer of an appropriate hydrophobicity and concentration into the refolding buffer. A molecular simulation of the refolding of a two-dimensional simple lattice protein was presented. A protein folding map recording the occurrence frequency of specified conformations was derived, from which the refolding thermodynamics and kinetics were interpreted. It is shown that, in the absence of polymer, the protein falls into the "energy trapped" conformations characterized by a high intramolecular hydrophobic interaction, denoted as HH contact, and a high magnitude of the structure overlap function, chi. This makes it difficult for the protein to fold to the native state. The polymer with a suitable chain length, concentration, and hydrophobicity has formed complex with partially folded protein and created diversified intermediates with low chi. This gives more pathways for the protein to fold to the native state. At a given hydrophobicity, the short chain polymer has a broader concentration range where it assists protein folding than those of long chains. The above simulation agrees well with the experimental results reported elsewhere [Cleland et al., J. Biol. Chem. 267, 13327 (1992); ibid., Bio/Technology 10, 1013 (1992); Chen et al., Enzyme Microb. Technol. 32, 120 (2003); Lu et al., Biochem. Eng. J. 24, 55 (2005); ibid., J. Chem. Phys. 122, 134902 (2005); ibid., Biochem. Eng. J. (to be published)] and is of fundamental importance for the design and application of polymers for protein refolding.     

Molecular dynamics for surfactant-assisted protein refolding

Surfactants are widely used to refold recombinant proteins that are produced as inclusion bodies in E. Coli. However, the microscopic details of the surfactant-assisted protein refolding processes are yet to be uncovered. In the present work, the authors aim to provide insights into the effect of hydrophobic interactions of a denatured protein with surfactant molecules on the refolding kinetics and equilibrium by using the Langevin dynamics for coarse-grained models. The authors have investigated the folding behavior of a beta-barrel protein in the presence of surfactants of different hydrophobicities and concentrations. It is shown that the protein folding process follows a "collapse-rearrangement" mechanism, i.e., the denatured protein first falls into a collapsed state before acquiring the native conformation. In comparison with the protein folding without surfactants, the protein-surfactant hydrophobic interactions promote the collapse of a denatured protein and, consequently, the formation of a hydrophobic core. However, the surfactants must be released from the hydrophobic core during the rearrangement step, in which the native conformation is formed. The simulation results can be qualitatively reproduced by experiments.     

Variationally determined free energy profiles for structural models of proteins: characteristic temperatures for folding and trapping

Characterizing the phase diagram for proteins is important both for laboratory studies and for the development of structure prediction algorithms. Using a variational scheme, we calculated the generic features of the protein thermostability over a large range of temperatures for a set of more than 50 different proteins using a model based on native structure alone. Focusing on a specific system, protein G, we further examined, using a more realistic model that includes the nonnative interaction, the thermostability of both the native state and a collection of trap structures. By surveying the native structures for many proteins and by paying closer attention to the various trap structures of protein G, we obtained an overall understanding of the folding dynamics far from the conditions usually focused on; namely, those near the folding temperature alone. Two characteristic temperatures (shown to scale with folding temperature in general) signal drastic changes in the folding mechanism. The variational calculations suggest that most proteins would, indeed, fold in a barrierless manner below a critical temperature analogous to a spinodal in crystallization. For fixed interaction strengths, this temperature, however, seems to be generally very low, approximately 50% of the equilibrium folding temperature. Likewise, native proteins, in general, would unfold in a completely barrierless way at a temperature 25% above folding temperature according to these variational calculations. We also studied the distribution of free energy profiles for escape from a set of trap structures generated by simulations.     

Computer simulations of the refolding of sperm whale apomyoglobin from high-temperature denaturated state

The refolding mechanism of apomyoglobin (apoMb) subsequent to high-temperature unfolding has been examined using computer simulations with atomic level detail. The folding of this protein has been extensively studied experimentally, providing a large database of folding parameters which can be probed using simulations. In the present study, 4-folding trajectories of apoMb were computed starting from coiled structures. A crystal structure of sperm whale myoglobin taken from the Protein Data Bank was used to construct the final native conformation by removal of the heme group followed by energy optimization. The initial unfolded conformations were obtained from high-temperature molecular dynamics simulations. Room-temperature refolding trajectories at neutral pH were obtained using the stochastic difference equation in length algorithm. The folding trajectories were compared with experimental results and two previous molecular dynamics studies at low pH. In contrast to the previous simulations, an extended intermediate with large helical content was not observed. In the present study, a structural collapse occurs without formation of helices or native contacts. Once the protein structure is more compact (radius of gyration<18 A) secondary and tertiary structures appear. These results suggest that apoMb follows a different folding pathway after high-temperature denaturation.     

Role of water mediated interactions in protein-protein recognition landscapes

The energy landscape picture of protein folding and binding is employed to optimize a number of pair potentials for direct and water-mediated interactions in protein complex interfaces. We find that water-mediated interactions greatly complement direct interactions in discriminating against various types of trap interactions that model those present in the cell. We highlight the context dependent nature of knowledge-based binding potentials, as contrasted with the situation for autonomous folding. By performing a Principal Component Analysis (PCA) of the corresponding interaction matrixes, we rationalize the strength of the recognition signal for each combination of the contact type and reference trap states using the differential in the idealized "canonical" amino acid compositions of native and trap layers. The comparison of direct and water-mediated contact potential matrixes emphasizes the importance of partial solvation in stabilizing charged groups in the protein interfaces. Specific water-mediated interresidue interactions are expected to influence significantly the kinetics as well as thermodynamics of protein association.     

Dynamic control of protein folding pathway with a polymer of tunable hydrophobicity

While the knowledge of protein folding in a dilute solution is now well-advanced, little is known of the influence of surrounding conditions on the folding kinetics, in particular when the protein is in a dynamically responsive environment. Here we report a new procedure to control the pathways of protein folding by using a thermally responsive polymer that varies its hydrophobicity concomitant with the protein structural changes. The advantages of folding in a dynamic environment have been demonstrated first by Langevin dynamics simulations on the basis of coarse-grained models for both the protein and polymer and then by experiments for lysozyme refolding in the presence of poly(N-isopropylacrylamide-co-N-tert-butylacrylamide), a thermal responsive polymer that varies its hydrophobicity in response to temperature. The simulation suggests that decreasing the polymer hydrophobicity during the folding process may result in an optimized free-energy landscape that enhances both the folding yield and kinetics. The experiments affirm that an optimal folding condition can be identified when structural transitions of the protein collaborate with the polymer hydrophobicity tuned by variation of temperature.     

Confinement in nanopores can destabilize α-helix folding proteins and stabilize the β structures

Protein folding in confined media has attracted wide attention over the past decade due to its importance in both in vivo and in vitro applications. Currently, it is generally believed that protein stability increases by decreasing the size of the confining medium, if its interaction with the confining walls is repulsive, and that the maximum folding temperature in confinement occurs for a pore size only slightly larger than the smallest dimension of the folded state of a protein. Protein stability in pore sizes, very close to the size of the folded state, has not however received the attention that it deserves. Using detailed, 0.3-ms-long molecular dynamics simulations, we show that proteins with an α-helix native state can have an optimal folding temperature in pore sizes that do not affect the folded-state structure. In contradiction to the current theoretical explanations, we find that the maximum folding temperature occurs in larger pores for smaller α-helices. In highly confined pores the free energy surface becomes rough, and a new barrier for protein folding may appear close to the unfolded state. In addition, in small nanopores the protein states that contain the β structures are entropically stabilized, in contrast to the bulk. As a consequence, folding rates decrease notably and the free energy surface becomes rougher. The results shed light on many recent experimental observations that cannot be explained by the current theories, and demonstrate the importance of entropic effects on proteins' misfolded states in highly confined environments. They also support the concept of passive effect of chaperonin GroEL on protein folding by preventing it from aggregation in crowded environment of biological cells, and provide deeper clues to the α → β conformational transition, believed to contribute to Alzheimer's and Parkinson's diseases. The strategy of protein and enzyme stabilization in confined media may also have to be revisited in the case of tight confinement. For in silico studies of protein folding in confined media, use of non-Go potentials may be more appropriate.     

Structural and thermodynamic investigations on the aggregation and folding of acylphosphatase by molecular dynamics simulations and solvation free energy analysis

Protein engineering method to study the mutation effects on muscle acylphosphatase (AcP) has been actively applied to describe kinetics and thermodynamics associated with AcP aggregation as well as folding processes. Despite the extensive mutation experiments, the molecular origin and the structural motifs for aggregation and folding kinetics as well as thermodynamics of AcP have not been rationalized at the atomic resolution. To this end, we have investigated the mutation effects on the structures and thermodynamics for the aggregation and folding of AcP by using the combination of fully atomistic, explicit-water molecular dynamics simulations, and three-dimensional reference interaction site model theory. The results indicate that the A30G mutant with the fastest experimental aggregation rate displays considerably decreased α1-helical contents as well as disrupted hydrophobic core compared to the wild-type AcP. Increased solvation free energy as well as hydrophobicity upon A30G mutation is achieved due to the dehydration of hydrophilic side chains in the disrupted α1-helix region of A30G. In contrast, the Y91Q mutant with the slowest aggregation rate shows a non-native H-bonding network spanning the mutation site to hydrophobic core and α1-helix region, which rigidifies the native state protein conformation with the enhanced α1-helicity. Furthermore, Y91Q exhibits decreased solvation free energy and hydrophobicity compared to wild type due to more exposed and solvated hydrophilic side chains in the α1-region. On the other hand, the experimentally observed slower folding rates in both mutants are accompanied by decreased helicity in α2-helix upon mutation. We here provide the atomic-level structures and thermodynamic quantities of AcP mutants and rationalize the structural origin for the changes that occur upon introduction of those mutations along the AcP aggregation and folding processes.     

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.     

Dynamic redox environment-intensified disulfide bond shuffling for protein refolding in vitro: molecular simulation and experimental validation

One challenge in protein refolding is to dissociate the non-native disulfide bonds and promote the formation of native ones. In this study, we present a coarse-grained off-lattice model protein containing disulfide bonds and simulate disulfide bond shuffling during the folding of this model protein. Introduction of disulfide bonds in the model protein led to enhanced conformational stability but reduced foldability in comparison to counterpart protein without disulfide bonds. The folding trajectory suggested that the model protein retained the two-step folding mechanism in terms of hydrophobic collapse and structural rearrangement. The disulfide bonds located in the hydrophobic core were formed before the collapsing step, while the bonds located on the protein surface were formed during the rearrangement step. While a reductive environment at the initial stage of folding favored the formation of native disulfide bonds in the hydrophobic core, an oxidative environment at a later stage of folding was required for the formation of disulfide bonds at protein surface. Appling a dynamic redox environment, that is, one that changes from reductive to oxidative, intensified disulfide bond shuffling and thus resulted in improved recovery of the native conformation. The above-mentioned simulation was experimentally validated by refolding hen-egg lysozyme at different urea concentrations and oxidized glutathione/reduced glutathione (GSSG/GSH) ratios, and an optimal redox environment, in terms of the GSSG to GSH ratio, was identified. The implementation of a dynamic redox environment by tuning the GSSG/GSH ratio further improved the refolding yield of lysozyme, as predicted by molecular simulation.     

Sampling of states for estimating the folding funnel entropy and energy landscape of a model alpha-helical hairpin peptide

Protein folding times are many orders of magnitude shorter than would occur if the peptide chain randomly sampled possible configurations, which implies that protein folding is a directed process. The detailed shape of protein's energy landscape determines the rate and reliability of folding to the native state, but the large number of structural degrees of freedom generates an energy landscape that is hard to visualize because of its high dimensionality. A commonly used picture is that of an energy funnel leading from high energy random coil state down to the low energy native state. As lattice computer models of protein dynamics become more realistic, the number of possible configurations becomes too large to count directly. Statistical mechanic and thermodynamic approaches allow us to count states in an approximate manner to quantify the entropy and energy of the energy landscape within a folding funnel for an alpha-helical protein. We also discuss the problems that arise in attempting to count the huge number of individual states of the random coil at the top of the funnel.     

Accelerated molecular dynamics simulations of protein folding

Folding of four fast‐folding proteins, including chignolin, Trp‐cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred‐of‐microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2–2.1 Å of the native NMR or X‐ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second‐order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein‐folding studies. © 2015 Wiley Periodicals, Inc.     

A component analysis of the free energies of folding of 35 proteins: A consensus view on the thermodynamics of folding at the molecular level

What factors favor protein folding? This is a textbook question. Parsing the experimental free energies of folding/unfolding into diverse enthalpic and entropic components of solute and solvent favoring or disfavoring folding is not an easy task. In this study, we present a computational protocol for estimating the free energy contributors to protein folding semi‐quantitatively using ensembles of unfolded and native states generated via molecular dynamics simulations. We tested the methodology on 35 proteins with diverse structural motifs and sizes and found that the calculated free energies correlate well with experiment (correlation coefficient ∼ 0.85), enabling us to develop a consensus view of the energetics of folding. As a more sensitive test of the methodology, we also investigated the free energies of folding of an additional 33 single point mutants and obtained a correlation coefficient of 0.8. A notable observation is that the folding free energy components appear to carry signatures of the fold (SCOP classification) of the protein. © 2017 Wiley Periodicals, Inc.     

Spectroscopic studies on the interaction of bovine (BSA) and human (HSA) serum albumins with ionic surfactants

Bovine (BSA) and human (HSA) serum albumins are frequently used in biophysical and biochemical studies since they have a similar folding, a well known primary structure, and they have been associated with the binding of many different categories of small molecules. One important difference of BSA and HSA is the fact that bovine albumin has two tryptophan residues while human albumin has a unique tryptophan. In this work results are presented for the interaction of BSA and HSA with several ionic surfactants, namely, anionic sodium dodecyl sulfate (SDS), cationic cethyltrimethylammonium chloride (CTAC) and zwitterionic N-hexadecyl-N,N-dimethyl-3-ammonium-1-propanesulfonate (HPS), as monitored by fluorescence spectroscopy of intrinsic tryptophans and circular dichroism spectroscopy. On the interaction of all three surfactants with BSA, at low concentrations, a quenching of fluorescence takes place and Stern-Volmer analysis allowed to estimate their 'effective' association constants to the protein: for SDS, CTAC and HPS at pH 7.0 these constants are, respectively, (1.4+/-0.1) x 10(5) M(-1), (8.9+/-0.1) x 10(3) M(-1) and (1.4+/-0.1) x 10(4) M(-1). A blue shift of maximum emission is observed from 345 to 330 nm upon surfactant binding. Analysis of fluorescence emission spectra allowed to separate three species in solution which were associated to native protein, a surfactant protein complex and partially denatured protein. The binding at low surfactant concentrations follows a Hill plot model displaying positive cooperativity and a number of surfactant binding sites very close to the number of cationic or anionic residues present in the protein. Circular dichroism data corroborated the partial loss of secondary structure upon surfactant addition showing the high stability of serum albumin. The interaction of the surfactants with HSA showed an enhancement of fluorescence at low concentrations, opposite to the effect on BSA, consistent with the existence of a unique buried tryptophan residue in this protein with considerable static quenching in the native state. The effects of surfactants at low concentrations were very similar to those of myristic acid suggesting a non specific binding through hydrophobic interaction modulated by eletrostatic interactions. The changes in the vicinity of the tryptophan residues are discussed based on the recently published crystallographic structure of HSA myristate complex (S. Curry et al., Nat. Struct. Biol. 5 (1998) 827).     

Enhanced Wang Landau sampling of adsorbed protein conformations

Using computer simulations to model the folding of proteins into their native states is computationally expensive due to the extraordinarily low degeneracy of the ground state. In this paper, we develop an efficient way to sample these folded conformations using Wang Landau sampling coupled with the configurational bias method (which uses an unphysical "temperature" that lies between the collapse and folding transition temperatures of the protein). This method speeds up the folding process by roughly an order of magnitude over existing algorithms for the sequences studied. We apply this method to study the adsorption of intrinsically disordered hydrophobic polar protein fragments on a hydrophobic surface. We find that these fragments, which are unstructured in the bulk, acquire secondary structure upon adsorption onto a strong hydrophobic surface. Apparently, the presence of a hydrophobic surface allows these random coil fragments to fold by providing hydrophobic contacts that were lost in protein fragmentation.     

KOH catalysed preparation of activated carbon aerogels for dye adsorption

Surfactants prevent the irreversible aggregation of partially refolded proteins, and they are also known to assist in protein refolding. A novel approach to protein refolding that utilizes a pair of low molecular weight folding assistants, a detergent and cyclodextrin, was proposed by Rozema and Gellman (D. Rozema, S.H. Gellman, J. Am. Chem. Soc. 117 (1995) 2373). We report the refolding of bovine serum albumin (BSA) assisted by these artificial chaperones, utilizing gemini surfactants for the first time. A combination of cationic gemini surfactants, bis(cetyldimethylammonium)pentane dibromide (C(16)H(33)(CH(3))(2)N(+)-(CH(2))(5)-N(+)(CH(3))(2)C(16)H(33)·2Br(-) designated as G5 and bis(cetyldimethylammonium)hexane dibromide (C(16)H(33)(CH(3))(2)N(+)-(CH(2))(6)-N(+)(CH(3))(2)C(16)H(33)·2Br(-) designated as G6 and cyclodextrins, was used to refold guanidinium chloride (GdCl) denatured BSA in the artificial chaperone assisted two step method. The single chain cationic surfactant cetyltrimethylammonium bromide (CTAB) was used for comparative studies. The studies were carried out in an aqueous medium at pH 7.0 using circular dichroism, dynamic light scattering and ANS binding studies. The denatured BSA was found to get refolded by very small concentrations of gemini surfactant at which the single chain counterpart was found to be ineffective. Different from the single chain surfactant, the gemini surfactants exhibit much stronger electrostatic and hydrophobic interactions with the protein and are thus effective at much lower concentrations. Based on the present study it is expected that gemini surfactants may prove useful in the protein refolding operations and may thus be effectively employed to circumvent the problem of misfolding and aggregation.     

Prediction of native-state hydrogen exchange from perfectly funneled energy landscapes

Simulations based on perfectly funneled energy landscapes often capture many of the kinetic features of protein folding. We examined whether simulations based on funneled energy functions can also describe fluctuations in native-state protein ensembles. We quantitatively compared the site-specific local stability determined from structure-based folding simulations, with hydrogen exchange protection factors measured experimentally for ubiquitin, chymotrypsin inhibitor 2, and staphylococcal nuclease. Different structural definitions for the open and closed states based on the number of native contacts for each residue, as well as the hydrogen-bonding state, or a combination of both criteria were evaluated. The predicted exchange patterns agree with the experiments under native conditions, indicating that protein topology indeed has a dominant effect on the exchange kinetics. Insights into the simplest mechanistic interpretation of the amide exchange process were thus obtained.