Unfolding of the cold shock protein studied with biased molecular dynamics

Terahertz Frequency Spectroscopy to Determine Cold Shock Protein Stability upon Solvation and Evaporation - A Molecular Dynamics Study

Infrared (IR) and Terahertz (THz) spectroscopy simulations were carried out using CHARMM35b2 to determine protein stability. The stabilities of three bacterial cold shock proteins (Csps) originating from mesophiles, thermophiles and hyper- thermophiles respectively were investigated in this study. The three different Csps were investigated by Normal-Mode analysis and Molecular Dynamics simulation of THz spectra using the Hessian matrix for solvated systems, interpreted in the harmonic approximation at optimum near-melting temperatures of each homologue, by incorporating differences in the hydrous and anhydrous states of the Csps. The results show slight variations in the large scale protein motion. However, the IR spectra of Csps observed at the low frequency saddle surface region, clearly distinguishes the thermophilic and mesophilic proteins based on their stability. Further studies on protein stability employing low-frequency collective modes have the potential to reveal functionally important conformational changes that are biologically significant.

Probing the protein-folding mechanism using denaturant and temperature effects on rate constants

Protein folding has been extensively studied, but many questions remain regarding the mechanism. Characterizing early unstable intermediates and the high-free-energy transition state (TS) will help answer some of these. Here, we use effects of denaturants (urea, guanidinium chloride) and temperature on folding and unfolding rate constants and the overall equilibrium constant as probes of surface area changes in protein folding. We interpret denaturant kinetic m-values and activation heat capacity changes for 13 proteins to determine amounts of hydrocarbon and amide surface buried in folding to and from TS, and for complete folding. Predicted accessible surface area changes for complete folding agree in most cases with structurally determined values. We find that TS is advanced (50-90% of overall surface burial) and that the surface buried is disproportionately amide, demonstrating extensive formation of secondary structure in early intermediates. Models of possible pre-TS intermediates with all elements of the native secondary structure, created for several of these proteins, bury less amide and hydrocarbon surface than predicted for TS. Therefore, we propose that TS generally has both the native secondary structure and sufficient organization of other regions of the backbone to nucleate subsequent (post-TS) formation of tertiary interactions. The approach developed here provides proof of concept for the use of denaturants and other solutes as probes of amount and composition of the surface buried in coupled folding and other large conformational changes in TS and intermediates in protein processes.

Discover binding pathways using the sliding binding-box docking approach: application to binding pathways of oseltamivir to avian influenza H5N1 neuraminidase

Drug binding and unbinding are transient processes which are hardly observed by experiment and difficult to analyze by computational techniques. In this paper, we employed a cost-effective method called "pathway docking" in which molecular docking was used to screen ligand-receptor binding free energy surface to reveal possible paths of ligand approaching protein binding pocket. A case study was applied on oseltamivir, the key drug against influenza a virus. The equilibrium pathways identified by this method are found to be similar to those identified in prior studies using highly expensive computational approaches.

Is there an en route folding intermediate for Cold shock proteins?

Cold shock proteins (Csps) play an important role in cold shock response of a diverse number of organisms ranging from bacteria to humans. Numerous studies of the Csp from various species showed that a two-state folding mechanism is conserved and the transition state (TS) appears to be very compact. However, the atomic details of the folding mechanism of Csp remain unclear. This study presents the folding mechanism of Csp in atomic detail using an all-atom Go model-based simulations. Our simulations predict that there may exist an en route intermediate, in which β strands 1-2-3 are well ordered and the contacts between β1 and β4 are almost developed. Such an intermediate might be too unstable to be detected in the previous fluorescence energy transfer experiments. The transition state ensemble has been determined from the P(fold) analysis and the TS appears even more compact than the intermediate state.

A Comparison of Three Perturbation Molecular Dynamics Methods for Modeling Conformational Transitions

Targeted, steered, and biased molecular dynamics (MD) are widely used methods for studying transition processes of biomolecules. They share the common feature of adding external perturbations along a conformational progress variable to guide the transition in a predefined direction in conformational space, yet differ in how these perturbations are applied. In the present paper, we report a comparison of these three methods on generating transition paths for two different processes: the unfolding of the B domain of protein A and a conformational transition of the catalytic domain of a Src kinase Lyn. Transition pathways were calculated with different simulation parameters including the choice of progress variable and the simulation length or biasing force constant. A comparison of the generated paths based on structural similarity finds that the three perturbation MD methods generate similar transition paths for a given progress variable in most cases. On the other hand, the path depends more strongly on the choice of progress variable used to move the system between the initial and final states. Potentials of mean force (PMF) were calculated starting from unfolding trajectories to estimate the relative probabilities of the paths. A lower PMF was found for the lowest biasing force constant with BMD.

Prediction of flexible/rigid regions from protein sequences using k-spaced amino acid pairs

Background

Traditionally, it is believed that the native structure of a protein corresponds to a global minimum of its free energy. However, with the growing number of known tertiary (3D) protein structures, researchers have discovered that some proteins can alter their structures in response to a change in their surroundings or with the help of other proteins or ligands. Such structural shifts play a crucial role with respect to the protein function. To this end, we propose a machine learning method for the prediction of the flexible/rigid regions of proteins (referred to as FlexRP); the method is based on a novel sequence representation and feature selection. Knowledge of the flexible/rigid regions may provide insights into the protein folding process and the 3D structure prediction.

Results

The flexible/rigid regions were defined based on a dataset, which includes protein sequences that have multiple experimental structures, and which was previously used to study the structural conservation of proteins. Sequences drawn from this dataset were represented based on feature sets that were proposed in prior research, such as PSI-BLAST profiles, composition vector and binary sequence encoding, and a newly proposed representation based on frequencies of k-spaced amino acid pairs. These representations were processed by feature selection to reduce the dimensionality. Several machine learning methods for the prediction of flexible/rigid regions and two recently proposed methods for the prediction of conformational changes and unstructured regions were compared with the proposed method. The FlexRP method, which applies Logistic Regression and collocation-based representation with 95 features, obtained 79.5% accuracy. The two runner-up methods, which apply the same sequence representation and Support Vector Machines (SVM) and Naïve Bayes classifiers, obtained 79.2% and 78.4% accuracy, respectively. The remaining considered methods are characterized by accuracies below 70%. Finally, the Naïve Bayes method is shown to provide the highest sensitivity for the prediction of flexible regions, while FlexRP and SVM give the highest sensitivity for rigid regions.

Conclusion

A new sequence representation that uses k-spaced amino acid pairs is shown to be the most efficient in the prediction of the flexible/rigid regions of protein sequences. The proposed FlexRP method provides the highest prediction accuracy of about 80%. The experimental tests show that the FlexRP and SVM methods achieved high overall accuracy and the highest sensitivity for rigid regions, while the best quality of the predictions for flexible regions is achieved by the Naïve Bayes method.

The folding mechanism of collagen-like model peptides explored through detailed molecular simulations

Collagen has a unique folding mechanism that begins with the formation of a triple-helical structure near its C terminus followed by propagation of this structure to the N terminus. To elucidate factors that affect the folding of collagen, we explored the folding pathway of collagen-like model peptides using detailed molecular simulations with explicit solvent. Using biased molecular dynamics we examined the latter stages of folding of a peptide model of native collagen, (Pro-Hyp-Gly)10, and a peptide that models a Gly --> Ser mutation found in several forms of osteogenesis imperfecta, (Pro-Hyp-Gly)3-Pro-Hyp-Ser-(Pro-Hyp-Gly)6. Starting from an unfolded state that contains a C-terminal nucleated trimer, (Pro-Hyp-Gly)10 folds to a structure where two of the three chains associate through water-mediated hydrogen bonds and the third is relatively separated from this dimer. Calculated free-energy profiles for folding from this intermediate to the final triple-helical structure suggest that further folding occurs at a rate of approximately one Pro-Hyp-Gly triplet per msec. In contrast, after 6 nsec of biased dynamics, the region N-terminal to the Ser residue in (Pro-Hyp-Gly)3-Pro-Hyp-Ser-(Pro-Hyp-Gly)6 folds to a structure where the three chains form close contacts near the N terminus, away from the mutation site. Further folding to an ideal triple-helical structure at the site of the mutation is unfavorable as the free energy of a triple-helical conformation at this position is more than 20 kcal/mol higher than that of a structure with unassociated chains. These data provide insights into the folding pathway of native collagen and the events underlying the formation of misfolded structures.

Substructural cooperativity and parallel versus sequential events during protein unfolding

According to the "old view," proteins fold along well-defined sequential pathways, whereas the "new view" sees protein folding as a highly parallel stochastic process on funnel-shaped energy landscapes. We have analyzed parallel and sequential processes on a large number of molecular dynamics unfolding trajectories of the protein CI2 at high temperatures. Using rigorous statistical measures, we quantify the degree of sequentiality on two structural levels. The unfolding process is highly parallel on the microstructural level of individual contacts. On a coarser, macrostructural level of contact clusters, characteristic parallel and sequential events emerge. These characteristic events can be understood from loop-closure dependencies between the contact clusters. A correlation analysis of the unfolding times of the contacts reveals a high degree of substructural cooperativity within the contact clusters.

Similarity and difference in the unfolding of thermophilic and mesophilic cold shock proteins studied by molecular dynamics simulations

Molecular dynamics simulations were performed to unfold a homologous pair of thermophilic and mesophilic cold shock proteins at high temperatures. The two proteins differ in just 11 of 66 residues and have very similar structures with a closed five-stranded antiparallel beta-barrel. A long flexible loop connects the N-terminal side of the barrel, formed by three strands (beta1-beta3), with the C-terminal side, formed by two strands (beta4-beta5). The two proteins were found to follow the same unfolding pathway, but with the thermophilic protein showing much slower unfolding. Unfolding started with the melting of C-terminal strands, leading to exposure of the hydrophobic core. Subsequent melting of beta3 and the beta-hairpin formed by the first two strands then resulted in unfolding of the whole protein. The slower unfolding of the thermophilic protein could be attributed to ion pair formation of Arg-3 with Glu-46, Glu-21, and the C-terminal. These ion pairs were also found to be important for the difference in folding stability between the pair of proteins. Thus electrostatic interactions appear to play similar roles in the difference in folding stability and kinetics between the pair of proteins.

Tryptophan rotamers as evidenced by X-ray, fluorescence lifetimes, and molecular dynamics modeling

We investigated the native-state dynamics of the Bacillus caldolyticus cold-shock protein mutant Bc-Csp L66E, using fluorescence and appropriate molecular dynamics methods. Two fluorescence lifetimes were found, the amplitudes of which agree very well with tryptophan rotamer populations, obtained from parallel tempering calculations. Rotamer lifetimes were predicted by transition-state theory from high-temperature simulations. Transition pathways were extracted from the transition rates between individual rotameric states. The molecular dynamics also reveal the loop fluctuations in the native state.

Src kinase activation: A switched electrostatic network

Src tyrosine kinases are essential in numerous cell signaling pathways, and improper functioning of these enzymes has been implicated in many diseases. The activity of Src kinases is regulated by conformational activation, which involves several structural changes within the catalytic domain (CD): the orientation of two lobes of CD; rearrangement of the activation loop (A-loop); and movement of an alpha-helix (alphaC), which is located at the interface between the two lobes, into or away from the catalytic cleft. Conformational activation was investigated using biased molecular dynamics to explore the transition pathway between the active and the down-regulated conformation of CD for the Src-kinase family member Lyn kinase, and to gain insight into the interdependence of these changes. Lobe opening is observed to be a facile motion, whereas movement of the A-loop motion is more complex requiring secondary structure changes as well as communication with alphaC. A key result is that the conformational transition involves a switch in an electrostatic network of six polar residues between the active and the down-regulated conformations. The exchange between interactions links the three main motions of the CD. Kinetic experiments that would demonstrate the contribution of the switched electrostatic network to the enzyme mechanism are proposed. Possible implications for regulation conferred by interdomain interactions are also discussed.

Loop-closure events during protein folding: rationalizing the shape of Phi-value distributions

In the past years, the folding kinetics of many small single-domain proteins has been characterized by mutational Phi-value analysis. In this article, a simple, essentially parameter-free model is introduced which derives folding routes from native structures by minimizing the entropic loop-closure cost during folding. The model predicts characteristic folding sequences of structural elements such as helices and beta-strand pairings. Based on few simple rules, the kinetic impact of these structural elements is estimated from the routes and compared to average experimental Phi-values for the helices and strands of 15 small, well-characterized proteins. The comparison leads on average to a correlation coefficient of 0.62 for all proteins with polarized Phi-value distributions, and 0.74 if distributions with negative average Phi-values are excluded. The diffuse Phi-value distributions of the remaining proteins are reproduced correctly. The model shows that Phi-value distributions, averaged over secondary structural elements, can often be traced back to entropic loop-closure events, but also indicates energetic preferences in the case of a few proteins governed by parallel folding processes.

Stability and fluctuations of amide hydrogen bonds in a bacterial cytochrome c: a molecular dynamics study

Molecular dynamics (MD) simulations on a bacterial cytochrome c were performed to investigate the lifetime and fluctuations of backbone hydrogen bonds and to correlate these data with protection factors for hydrogen exchange measured by NMR spectroscopy (Bartalesi et al. in Biochemistry, 42:10923-10930, 2003). The MD simulations provide a consistent pattern in that long lifetimes of hydrogen bonds go along with small amplitude fluctuations. In agreement with experiments, differences in stability were found with a rather flexible N-terminal segment as compared with a more rigid C-terminal part. Protection factors of backbone hydrogen exchange correlate strongly with the number of contacts but also with hydrogen-bond occupancy, hydrogen-bond survival times, as well as the inverse of fluctuations of backbone atoms and hydrogen-bond lengths derived from MD simulation data. We observed a conformational transition in the C-terminal loop, and significant motion in the N-terminal loop, which can be interpreted as being the structural units involved in the onset of the protein unfolding process in agreement with experimental evidence on mitochondrial cytochrome c.

The Folding Transition State of the Cold Shock Protein is Strongly Polarized

The cold shock protein CspB from Bacillus subtilis consists of a three-stranded (beta1-beta3) and a two stranded (beta4-beta5) sheet, which form a closed beta barrel structure. CspB folds and unfolds rapidly in a two-state reaction, and the unfolded and the folded molecules interconvert with a time constant of 30 ms at the midpoint of the urea-induced transition (at 25 degrees C). The transition state of folding is native-like, as judged by the Tanford betaT value of > or =0.9. By using a mutational approach and Phi value analysis, we find that the folding transition state of CspB is energetically polarized. Despite the high betaT value, most Phi values are low. Values close to 1 were found for only a few residues, particularly in strand beta1 (Lys5, Val6, Lys7, Asn10). The interactions of the Asn10 side-chain with the backbone at positions 12 and 13 define the turn that connects the strands beta1 and beta2. Lys5 and Val6 in beta1 interact with residues in beta4, and their high Phi values indicate that an energetic linkage between beta1 and beta4 and thus between the two sheets exists already in the transition state. We compared our experimental Phi values with theoretical predictions of the folding pathway of cold shock proteins. Several of them suggest that the entire first sheet is formed in the transition state, and some identify the beta1-beta4 pairing as a crucial step in folding. Alternative paths that involve formation of the second sheet and beta3-beta5 pairing reactions were, however, suggested as well. The calculations gave coarse-grained pictures that are limited in resolution to the two sheets of CspB or to the elements of secondary structure. They did not identify the key residues with the high Phi values within these structural elements.

Structural and dynamic effects of α-Helix deletion in Sso7d: Implications for protein thermal stability

Sso7d is a 62-residue protein from the hyperthemophilic archaeon Sulfolobus solfataricus with a denaturation temperature close to 100 degrees C around neutral pH. An engineered form of Sso7d truncated at leucine 54 (L54Delta) is significantly less stable, with a denaturation temperature of 53 degrees C. Molecular dynamics (MD) studies of Sso7d and its truncated form at two different temperatures have been performed. The results of the MD simulations at 300 K indicate that: (1) the flexibility of Sso7d chain at 300 K agrees with that detected from X-ray and NMR structural studies; (2) L54Delta remains stable in the native folded conformation and possesses an overall dynamic behavior similar to that of the parent protein. MD simulations performed at 500 K, 10 ns long, indicate that, while Sso7d is in-silico resistant to high temperature, the truncated variant partially unfolds, revealing the early phases of the thermal unfolding pathway of the protein. Analysis of the trajectories of L54Delta suggests that the unzipping of the N-terminal and C-terminal beta-strands should be the first event of the unfolding pathway, and points out the regions more resistant to thermal unfolding. These findings allow one to understand the role played by specific interactions connecting the two ends of the chain for the high thermal stability of Sso7d, and support recent hypotheses on its folding mechanism emerged from site-directed mutagenesis studies.