Folding processes of the B domain of protein A to the native state observed in all-atom ab initio folding simulations

The B domain of protein A retains residual structures in 6 M guanidinium chloride as revealed by hydrogen/deuterium-exchange NMR spectroscopy

The characterization of residual structures persistent in unfolded proteins is an important issue in studies of protein folding, because the residual structures present, if any, may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the residual structures of the isolated B domain (BDPA) of staphylococcal protein A in 6 M guanidinium chloride. BDPA is a small three-helix-bundle protein, and until recently its folding/unfolding reaction has been treated as a simple two-state process between the native and the fully unfolded states. We employed a dimethylsulfoxide (DMSO)-quenched hydrogen/deuterium (H/D)-exchange 2D NMR techniques with the use of spin desalting columns, which allowed us to investigate the H/D-exchange behavior of individually identified peptide amide (NH) protons. We obtained H/D-exchange protection factors of the 21 NH protons that form an α-helical hydrogen bond in the native structure, and the majority of these NH protons were significantly protected with a protection factor of 2.0-5.2 in 6 M guanidinium chloride, strongly suggesting that these weakly protected NH protons form much stronger hydrogen bonds under native folding conditions. The results can be used to deduce the structure of an early folding intermediate, when such an intermediate is shown by other methods. Among three native helical regions, the third helix in the C-terminal side was highly protected and stabilized by side-chain salt bridges, probably acting as the folding initiation site of BDPA. The present results are discussed in relation to previous experimental and computational findings on the folding mechanisms of BDPA.

Evaluation of implicit solvent models in molecular dynamics simulation of α-Synuclein

We report conventional and accelerated molecular dynamics simulations of α-Synuclein, designed to assess performance of using different starting conformation, solvation environment and force field combination. Backbone and sidechain chemical shifts, radius of gyration, presence of β-hairpin structures in KTK(E/Q)GV repeats and secondary structure percentages were used to evaluate how variations in forcefield, solvation model and simulation protocol provide results that correlate with experimental findings. We show that with suitable choice of forcefield and solvent, ff03ws and OBC implicit model, respectively, acceptable reproduction of experimental data on size and secondary structure is obtained by both conventional and accelerated MD. In contrast to the implicit solvent model, simulations in explicit TIP4P/2005 solvent do not properly represent size or secondary structure of α-Synuclein.Communicated by Ramaswamy H. Sarma.

Implications of the unfolded state in the folding energetics of heterogeneous-backbone protein mimetics

Sequence-encoded folding is the foundation of protein structure and is also possible in synthetic chains of artificial chemical composition. In natural proteins, the characteristics of the unfolded state are as important as those of the folded state in determining folding energetics. While much is known about folded structures adopted by artificial protein-like chains, corresponding information about the unfolded states of these molecules is lacking. Here, we report the consequences of altered backbone composition on the structure, stability, and dynamics of the folded and unfolded states of a compact helix-rich protein. Characterization through a combination of biophysical experiments and atomistic simulation reveals effects of backbone modification that depend on both the type of artificial monomers employed and where they are applied in sequence. In general, introducing artificial connectivity in a way that reinforces characteristics of the unfolded state ensemble of the prototype natural protein minimizes the impact of chemical changes on folded stability. These findings have implications in the design of protein mimetics and provide an atomically detailed picture of the unfolded state of a natural protein and artificial analogues under non-denaturing conditions.

Biophysical experiments and atomistic simulation reveal impacts of protein backbone alteration on the ensemble that defines the unfolded state. These effects have implications on folded stability of protein mimetics.

Protein A superantigen: structure, engineering and molecular basis of antibody recognition

Staphylococcus aureus interacts with the human immune system through the production of secreted factors. Key among these is protein A, a B-cell superantigen capable of interacting with both antibody Fc and VH regions. Here, we review structural and molecular features of this important example of naturally occurring bacterial superantigens, as well as engineered variants and their application in biotechnology.

The effects of implicit modeling of nonpolar solvation on protein folding simulations

Implicit solvent models, in which the polar and nonpolar solvation free-energies of solute molecules are treated separately, have been widely adopted for molecular dynamics simulation of protein folding. While the development of the implicit models is mainly focused on the methodological improvement and key parameter optimization for polar solvation, nonpolar solvation has been either ignored or described by a simplistic surface area (SA) model. In this work, we performed the folding simulations of multiple β-hairpin and α-helical proteins with varied surface tension coefficients embedded in the SA model to clearly demonstrate the effects of nonpolar solvation treated by a popular SA model on protein folding. The results indicate that the change in the surface tension coefficient does not alter the ability of implicit solvent simulations to reproduce a protein native structure but indeed controls the components of the equilibrium conformational ensemble and modifies the energetic characterization of the folding transition pathway. The suitably set surface tension coefficient can yield explicit solvent simulations and/or experimentally suggested folding mechanism of protein. In addition, the implicit treatment of both polar and nonpolar components of solvation free-energy contributes to the overestimation of the secondary structure in implicit solvent simulations.

Simultaneous Determination of Two Subdomain Folding Rates Using the "Transfer-Quench" Method

The investigation of the mechanism of protein folding is complicated by the context dependence of the rates of intramolecular contact formation. Methods based on site-specific labeling and ultrafast spectroscopic detection of fluorescence signals were developed for monitoring the rates of individual subdomain folding transitions in situ, in the context of the whole molecule. However, each site-specific labeling modification might affect rates of folding of near-neighbor structural elements, and thus limit the ability to resolve fine differences in rates of folding of these elements. Therefore, it is highly desirable to be able to study the rates of folding of two or more neighboring subdomain structures using a single mutant to facilitate resolution of the order and interdependence of such steps. Here, we report the development of the "Transfer-Quench" method for measuring the rate of formation of two structural elements using a single triple-labeled mutant. This method is based on Förster resonance energy transfer combined with fluorescence quenching. We placed the donor and acceptor at the loop ends, and a quencher at an α-helical element involved in the node forming the loop. The folding of the triple-labeled mutant is monitored by the acceptor emission. The formation of nonlocal contact (loop closure) increases the time-dependent acceptor emission, while the closure of the labeled helix turn reduces this emission. The method was applied in a study of the folding mechanism of the common model protein, the B domain of staphylococcal protein A. Only natural amino acids were used as probes, and thus possible structural perturbations were minimized. Tyr and Trp residues served as donor and acceptor at the ends of a long loop between helices I and II, and a Cys residue as a quencher for the acceptor. We found that the closure of the loop (segment 14-33) occurs with the same rate constant as the nucleation of helix HII (segment 33-29), in line with the nucleation-condensation model.

Aggregation of Chameleon Peptides: Implications of α-Helicity in Fibril Formation

We investigate the relationship between the inherent secondary structure and aggregation propensity of peptides containing chameleon sequences (i.e., sequences that can adopt either α or β structure depending on context) using a combination of replica exchange molecular dynamics simulations, ion-mobility mass spectrometry, circular dichroism, and transmission electron microscopy. We focus on an eight-residue long chameleon sequence that can adopt an α-helical structure in the context of the iron-binding protein from Bacillus anthracis (PDB id 1JIG ) and a β-strand in the context of the baculovirus P35 protein (PDB id 1P35 ). We show that the isolated chameleon sequence is intrinsically disordered, interconverting between α-helical and β-rich conformations. The inherent conformational plasticity of the sequence can be constrained by addition of flanking residues with a given secondary structure propensity. Intriguingly, we show that the chameleon sequence with helical flanking residues aggregates rapidly into fibrils, whereas the chameleon sequence with flanking residues that favor β-conformations has weak aggregation propensity. This work sheds new insights into the possible role of α-helical intermediates in fibril formation.

Dissecting protein architecture with communication blocks and communicating segment pairs

BACKGROUND

Proteins adapt to environmental conditions by changing their shape and motions. Characterising protein conformational dynamics is increasingly recognised as necessary to understand how proteins function. Given a conformational ensemble, computational tools are needed to extract in a systematic way pertinent and comprehensive biological information.

RESULTS

Here, we present a method, Communication Mapping (COMMA), to decipher the dynamical architecture of a protein. The method first extracts residue-based dynamic properties from all-atom molecular dynamics simulations. Then, it integrates them in a graph theoretic framework, where it identifies groups of residues or protein regions that mediate short- and long-range communication. COMMA introduces original concepts to contrast the different roles played by these regions, namely communication blocks and communicating segment pairs, and evaluates the connections and communication strengths between them. We show the utility and capabilities of COMMA by applying it to three archetypal proteins, namely protein A, the tyrosine kinase KIT and the tumour suppressor p53.

CONCLUSION

Our method permits to compare in a direct way the dynamical behaviour either of proteins with different characteristics or of the same protein in different conditions. It is useful to identify residues playing a key role in protein allosteric regulation and to explain the effects of deleterious mutations in a mechanistic way. COMMA is a fully automated tool with broad applicability. It is freely available to the community at www.lcqb.upmc.fr/COMMA .

Partial unfolding and refolding for structure refinement: A unified approach of geometric simulations and molecular dynamics

The most successful protein structure prediction methods to date have been template-based modeling (TBM) or homology modeling, which predicts protein structure based on experimental structures. These high accuracy predictions sometimes retain structural errors due to incorrect templates or a lack of accurate templates in the case of low sequence similarity, making these structures inadequate in drug-design studies or molecular dynamics simulations. We have developed a new physics based approach to the protein refinement problem by mimicking the mechanism of chaperons that rehabilitate misfolded proteins. The template structure is unfolded by selectively (targeted) pulling on different portions of the protein using the geometric based technique FRODA, and then refolded using hierarchically restrained replica exchange molecular dynamics simulations (hr-REMD). FRODA unfolding is used to create a diverse set of topologies for surveying near native-like structures from a template and to provide a set of persistent contacts to be employed during re-folding. We have tested our approach on 13 previous CASP targets and observed that this method of folding an ensemble of partially unfolded structures, through the hierarchical addition of contact restraints (that is, first local and then nonlocal interactions), leads to a refolding of the structure along with refinement in most cases (12/13). Although this approach yields refined models through advancement in sampling, the task of blind selection of the best refined models still needs to be solved. Overall, the method can be useful for improved sampling for low resolution models where certain of the portions of the structure are incorrectly modeled.

Elucidation of the Aggregation Pathways of Helix-Turn-Helix Peptides: Stabilization at the Turn Region Is Critical for Fibril Formation

Aggregation of proteins to fiberlike aggregates often involves a transformation of native monomers to β-sheet-rich oligomers. This general observation underestimates the importance of α-helical segments in the aggregation cascade. Here, using a combination of experimental techniques and accelerated molecular dynamics simulations, we investigate the aggregation of a 43-residue, apolipoprotein A-I mimetic peptide and its E21Q and D26N mutants. Our study indicates a strong propensity of helical segments not to adopt cross-β-fibrils. The helix-turn-helix monomeric conformation of the peptides is preserved in the mature fibrils. Furthermore, we reveal opposite effects of mutations on and near the turn region in the self-assembly of these peptides. We show that the E21-R24 salt bridge is a major contributor to helix-turn-helix folding, subsequently leading to abundant fibril formation. On the other hand, the K19-D26 interaction is not required to fold the native helix-turn-helix peptide. However, removal of the charged D26 residue decreases the stability of the helix-turn-helix monomer and consequently reduces the level of aggregation. Finally, we provide a more refined assembly model for the helix-turn-helix peptides from apolipoprotein A-I based on the parallel stacking of helix-turn-helix dimers.

Fast helix formation in the B domain of protein A revealed by site-specific infrared probes

Comparison of experimental and computational protein folding studies can be difficult because of differences in structural resolution. Isotope-edited infrared spectroscopy offers a direct measure of structural changes involved in protein folding at the single-residue level. Here we demonstrate the increased resolution of site-specific infrared probes to the peptide backbone in the B domain of staphylococcal protein A (BdpA). (13)C═(18)O-labeled methionine was incorporated into each of the helices using recombinant protein expression. Laser-induced temperature jumps coupled with infrared spectroscopy were used to probe changes in the peptide backbone on the submillisecond time scale. The relaxation kinetics of the buried helices, solvated helices, and labeled positions were measured independently by probing the corresponding bands assigned in the amide I region. Using these wavelength-dependent measurements, we observe a fast nanosecond phase and slower microsecond phase at each position. We find at least partial formation of helices 1-3 in the fast intermediate state that precedes the transition state. These measurements provide direct, time-resolved experimental evidence of the early formation of partial helical structure in helices 1 and 3, supporting folding models proposed by computer simulations.

Important roles of hydrophobic interactions in folding and charge interactions in misfolding of α-helix bundle protein

An enhanced-sampling molecular dynamics simulation is presented to quantitatively demonstrate the important roles of hydrophobic and charge interactions in the folding and misfolding of α-helix bundle protein, respectively.

Quantifying hub-like behavior in protein folding networks

The free energy landscape of a protein is a function of many interdependent degrees of freedom. For this reason, conceptual constructs (e.g., funnels) have been useful to visualize these landscapes. One relatively new construct is the idea of a hub-like native state that is the final destination of many non-interconverting folding pathways. This is in contrast to the idea of a single predominant folding pathway connecting the native state to a rapidly interconverting ensemble of unfolded states. The key quantity to distinguish between these two ideas is the connectivity of the unfolded ensemble. We present a metric to determine this connectivity for a given network, which can be calculated either from continuous folding trajectories, or a Markov model. The metric determines how often a region of space is used as an intermediate on transition paths that connect two other regions of space, and we use it here to determine how often two parts of the unfolded ensemble are connected directly, versus how often these transitions are mediated by the native state.

Smoothing protein energy landscapes by integrating folding models with structure prediction

Decades of work has investigated the energy landscapes of simple protein models, but what do the landscapes of real, large, atomically detailed proteins look like? We explore an approach to this problem that systematically extracts simple funnel models of actual proteins using ensembles of structure predictions and physics-based atomic force fields and sampling. Central to our effort are calculations of a quantity called the relative entropy, which quantifies the extent to which a given set of structure decoys and a putative native structure can be projected onto a theoretical funnel description. We examine 86 structure prediction targets and one coupled folding-binding system, and find that in a majority of cases the relative entropy robustly signals which structures are nearest to native (i.e., which appear to lie closest to a funnel bottom). Importantly, the landscape model improves substantially upon purely energetic measures in scoring decoys. Our results suggest that physics-based models-including both folding theories and all-atom force fields-may be successfully integrated with structure prediction efforts. Conversely, detailed predictions of structures and the relative entropy approach enable one to extract coarse topographic features of protein landscapes that may enhance the development and application of simpler folding models.

Computational model for protein unfolding simulation

The protein folding problem is one of the fundamental and important questions in molecular biology. However, the all-atom molecular dynamics studies of protein folding and unfolding are still computationally expensive and severely limited by the time scale of simulation. In this paper, a simple and fast protein unfolding method is proposed based on the conformational stability analyses and structure modeling. In this method, two structure-based conditions are considered to identify the unstable regions of proteins during the unfolding processes. The protein unfolding trajectories are mimicked through iterative structure modeling according to conformational stability analyses. Two proteins, chymotrypsin inhibitor 2 (CI2) and α -spectrin SH3 domain (SH3) were simulated by this method. Their unfolding pathways are consistent with the previous molecular dynamics simulations. Furthermore, the transition states of the two proteins were identified in unfolding processes and the theoretical Φ values of these transition states showed significant correlations with the experimental data (the correlation coefficients are >0.8). The results indicate that this method is effective in studying protein unfolding. Moreover, we analyzed and discussed the influence of parameters on the unfolding simulation. This simple coarse-grained model may provide a general and fast approach for the mechanism studies of protein folding.

Folding network of villin headpiece subdomain

Protein folding is a complex multidimensional process that is difficult to illustrate by the traditional analyses based on one- or two-dimensional profiles. Analyses based on transition networks have become an alternative approach that has the potential to reveal detailed features of protein folding dynamics. However, due to the lack of successful reversible folding of proteins from conventional molecular-dynamics simulations, this approach has rarely been utilized. Here, we analyzed the folding network from several 10 μs conventional molecular-dynamics reversible folding trajectories of villin headpiece subdomain (HP35). The folding network revealed more complexity than the traditional two-dimensional map and demonstrated a variety of conformations in the unfolded state, intermediate states, and the native state. Of note, deep enthalpic traps at the unfolded state were observed on the folding landscape. Furthermore, in contrast to the clear separation of the native state and the primary intermediate state shown on the two-dimensional map, the two states were mingled on the folding network, and prevalent interstate transitions were observed between these two states. A more complete picture of the folding mechanism of HP35 emerged when the traditional and network analyses were considered together.

Trends in template/fragment-free protein structure prediction

Predicting the structure of a protein from its amino acid sequence is a long-standing unsolved problem in computational biology. Its solution would be of both fundamental and practical importance as the gap between the number of known sequences and the number of experimentally solved structures widens rapidly. Currently, the most successful approaches are based on fragment/template reassembly. Lacking progress in template-free structure prediction calls for novel ideas and approaches. This article reviews trends in the development of physical and specific knowledge-based energy functions as well as sampling techniques for fragment-free structure prediction. Recent physical- and knowledge-based studies demonstrated that it is possible to sample and predict highly accurate protein structures without borrowing native fragments from known protein structures. These emerging approaches with fully flexible sampling have the potential to move the field forward.

Folding simulations of a de novo designed protein with a betaalphabeta fold

betaalphabeta structural motifs are commonly used building blocks in protein structures containing parallel beta-sheets. However, to our knowledge, no stand-alone betaalphabeta structure has been observed in nature to date. Recently, for the first time that we know of, a small protein with an independent betaalphabeta structure (DS119) was successfully designed in our laboratory. To understand the folding mechanism of DS119, in the study described here, we carried out all-atom molecular dynamics and coarse-grained simulations to investigate its folding pathways and energy landscape. From all-atom simulations, we successfully observed the folding event and got a stable folded structure with a minimal root mean-square deviation of 2.6 A with respect to the NMR structure. The folding process can be described as a fast collapse phase followed by rapid formation of the central helix, and then slow formation of a parallel beta-sheet. By using a native-centric Gō-like model, the cooperativity of the system was characterized in terms of the calorimetric criterion, sigmoidal transitions, conformation distribution shifts, and free-energy profiles. DS119 was found to be an incipient downhill folder that folds more cooperatively than a downhill folder, but less cooperatively than a two-state folder. This may reflect the balance between the two structural elements of DS119: the rapidly formed alpha-helix and the slowly formed parallel beta-sheet. Folding times estimated from both the all-atom simulations and the coarse-grained model were at microsecond level, making DS119 another fast folder. Compared to fast folders reported previously, DS119 is, to the best of our knowledge, the first that exhibits a parallel beta-sheet.

Characterizing the denatured state of human prion 121-230

Misfolding and aggregation of the prion protein (PrP) are responsible for the development of fatal transmissible neurodegenerative diseases. PrP undergoes structural conversion from a natively folded state into a misfolded state, resulting in insoluble amyloid fibrils. Partial unfolding has been recognized as an essential step in fibrillation. The strong correlation of unfolding and fibrillation emphasizes the importance of denatured states. To gain insight into possible aggregation-prone denatured states, we characterized the denatured state of human prion (huPrP) 121-230 near extended conformation by self-guided Langevin dynamics simulations. Our results revealed that denatured huPrP is partially folded with alpha-helical structure.

Temperature-induced unfolding of epidermal growth factor (EGF): insight from molecular dynamics simulation

Thermal disruption of protein structure and function is a potentially powerful therapeutic vehicle. With the emerging nanoparticle-targeting and femtosecond laser technology, it is possible to deliver heating locally to specific molecules. It is therefore important to understand how fast a protein can unfold or lose its function at high temperatures, such as near the water boiling point. In this study, the thermal damage of EGF was investigated by combining the replica exchange (136 replicas) and conventional molecular dynamics simulations. The REMD simulation was employed to rigorously explore the free-energy landscape of EGF unfolding. Interestingly, besides the native and unfolded states, we also observed a distinct molten globule (MG) state that retained substantial amount of native contacts. Based on the understanding that which the unfolding of EGF is a three-state process, we have examined the unfolding kinetics of EGF (N-->MG and MG-->D) with multiple 20-ns conventional MD simulations. The Arrhenius prefactors and activation energy barriers determined from the simulation are within the range of previously studied proteins. In contrast to the thermal damage of cells and tissues which take place on the time scale of seconds to hours at relatively low temperatures, the denaturation of proteins occur in nanoseconds when the temperature of heat bath approaches the boiling point.

Nearly symmetrical proteins: folding pathways and transition states

The folding pathways of the B domain of protein A have been the subject of many experimental and computational studies. Based on a statistical mechanical model, it has been suggested that the native state symmetry leads to multiple pathways, highly dependent on temperature and denaturant concentration. Experiments, however, have not confirmed this scenario. By considering four nearly symmetrical proteins, one of them being the above molecule, here we show that, if contact energies are properly taken into account, a different picture emerges from kinetic simulations of the above-mentioned model. This is characterized by a dominant folding pathway, which is consistent with the most recent experimental results. Given the simplicity of the model, we also report on a direct sampling of the transition state.

Dual folding pathways of an alpha/beta protein from all-atom ab initio folding simulations

Successful ab initio folding of proteins with both alpha-helix and beta-sheet requires a delicate balance among a variety of forces in the simulation model, which may explain that the successful folding of any alpha/beta proteins to within experimental error has yet to be reported. Here we demonstrate that it is an achievable goal to fold alpha/beta proteins with a force field emphasizing the balance between the two major secondary structures. Using our newly developed force field, we conducted extensive ab initio folding simulations on an alpha/beta protein full sequence design (FSD) employing both conventional molecular dynamics and replica exchange molecular dynamics in combination with a generalized-Born solvation model. In these simulations, the folding of FSD to the native state with high population (>64.2%) and high fidelity (C(alpha)-Root Mean Square Deviation of 1.29 A for the most sampled conformation when compared to the experimental structure) was achieved. The folding of FSD was found to follow two pathways. In the major pathway, the folding started from the formation of the helix. In the minor pathway, however, folding of the beta-hairpin started first. Further examination revealed that the helix initiated from the C-terminus and propagated toward the N-terminus. The formation of the hydrophobic contacts coincided with the global folding. Therefore the hydrophobic force does not appear to be the driving force of the folding of this protein.

SPlitSSI-SVM: an algorithm to reduce the misleading and increase the strength of domain signal

Protein domains contain information about the prediction of protein structure, function, evolution and design since the protein sequence may contain several domains with different or the same copies of the protein domain. In this study, we proposed an algorithm named SplitSSI-SVM that works with the following steps. First, the training and testing datasets are generated to test the SplitSSI-SVM. Second, the protein sequence is split into subsequence based on order and disorder regions. The protein sequence that is more than 600 residues is split into subsequences to investigate the effectiveness of the protein domain prediction based on subsequence. Third, multiple sequence alignment is performed to predict the secondary structure using bidirectional recurrent neural networks (BRNN) where BRNN considers the interaction between amino acids. The information of about protein secondary structure is used to increase the protein domain boundaries signal. Lastly, support vector machines (SVM) are used to classify the protein domain into single-domain, two-domain and multiple-domain. The SplitSSI-SVM is developed to reduce misleading signal, lower protein domain signal caused by primary structure of protein sequence and to provide accurate classification of the protein domain. The performance of SplitSSI-SVM is evaluated using sensitivity and specificity on single-domain, two-domain and multiple-domain. The evaluation shows that the SplitSSI-SVM achieved better results compared with other protein domain predictors such as DOMpro, GlobPlot, Dompred-DPS, Mateo, Biozon, Armadillo, KemaDom, SBASE, HMMPfam and HMMSMART especially in two-domain and multiple-domain.

Structural stability of electrosprayed proteins: temperature and hydration effects

Electrospray ionization is a gentle method for sample delivery, routinely used in gas-phase studies of proteins. It is crucial for structural investigations that the protein structure is preserved, and a good understanding of how structure is affected by the transition to the gas phase is needed for the tuning of experiments to meet that requirement. Small amounts of residual solvent have been shown to protect the protein, but temperature is important too, although it is not well understood how the latter affects structural details. Using molecular dynamics we have simulated four sparingly hydrated globular proteins (Trp-cage; Ctf, a C-terminal fragment of a bacterial ribosomal protein; ubiquitin; and lysozyme) in vacuum starting at temperatures ranging from 225 K to 425 K. For three of the proteins, our simulations show that a water layer corresponding to 3 A preserves the protein structure in vacuum, up to starting temperatures of 425 K. Only Ctf shows minor secondary structural changes at lower starting temperatures. The structural conservation stems mainly from interactions with the surrounding water. Temperature scales in simulations are not directly translatable into experiments, but the wide temperature range in which we find the proteins to be stable is reassuring for the success of future single particle imaging experiments. The water molecules aggregate in clusters and form patterns on the protein surface, maintaining a reproducible hydrogen bonding network. The simulations were performed mainly using OPLS-AA/L, with cross checks using AMBER03 and GROMOS96 53a6. Only minor differences between the results from the three different force fields were observed.

The fast-folding HP35 double mutant has a substantially reduced primary folding free energy barrier

The LYS24/29NLE double mutant of villin headpiece subdomain (HP35) is the fastest folding protein known so far with a folding time constant of 0.6 micros. In this work, the folding mechanism of the mutant has been investigated by both conventional and replica exchange molecular dynamics (CMD and REMD) simulations with AMBER FF03 force field and a generalized-Born solvation model. Direct comparison to the ab initio folding of the wild type HP35 enabled a close examination on the mutational effect on the folding process. The mutant folded to the native state, as demonstrated by the 0.50 A C(alpha)-root mean square deviation (RMSD) sampled in both CMD and REMD simulations and the high population of the folded conformation compared with the denatured conformations. Consistent with experiments, the significantly reduced primary folding free energy barrier makes the mutant closer to a downhill folder than the wild type HP35 that directly leads to the faster transition and higher melting temperature. However, unlike the proposed downhill folding which envisages a smooth shift between unfolded and folded states without transition barrier, we observed a well-defined folding transition that was consistent with experiments. Further examination of the secondary structures revealed that the two mutated residues have higher intrinsic helical preference that facilitated the formation of both helix III and the intermediate state which contains the folded segment helix II/III. Other factors contributing to the faster folding include the more favorable electrostatic interactions in the transition state with the removal of the charged NH(3)(+) groups from LYS. In addition, both transition state ensemble and denatured state ensemble are shifted in the mutant.