Desolvation and Development of Specific Hydrophobic Core Packing during Im7 Folding

Development of a tightly packed hydrophobic core drives the folding of water-soluble globular proteins and is a key determinant of protein stability. Despite this, there remains much to be learnt about how and when the hydrophobic core becomes desolvated and tightly packed during protein folding. We have used the bacterial immunity protein Im7 to examine the specificity of hydrophobic core packing during folding. This small, four-helix protein has previously been shown to fold via a compact three-helical intermediate state. Here, overpacking substitutions, in which residue side-chain size is increased, were used to examine the specificity and malleability of core packing in the folding intermediate and rate-limiting transition state. In parallel, polar groups were introduced into the Im7 hydrophobic core via Val→Thr or Phe→Tyr substitutions and used to determine the solvation status of core residues at different stages of folding. Over 30 Im7 variants were created allowing both series of substitutions to cover all regions of the protein structure. Φ-value analysis demonstrated that the major changes in Im7 core solvation occur prior to the population of the folding intermediate, with key regions involved in docking of the short helix III remaining solvent-exposed until after the rate-limiting transition state has been traversed. In contrast, overpacking core residues revealed that some regions of the native Im7 core are remarkably malleable to increases in side-chain volume. Overpacking residues in other regions of the Im7 core result in substantial (> 2.5 kJ mol^− 1) destabilisation of the native structure or even prevents efficient folding to the native state. This study provides new insights into Im7 folding; demonstrating that whilst desolvation occurs early during folding, adoption of a specifically packed core is achieved only at the very last step in the folding mechanism.     

Ab initio folding simulation of the Trp-cage mini-protein approaches NMR resolution

Here, we report a 100 ns molecular dynamics simulation of the folding process of a recently designed autonomous-folding mini-protein designated as tc5b with a new AMBER force field parameter set developed based on condensed-phase quantum mechanical calculations and a Generalized Born continuum solvent model. Starting from its fully extended conformation, our simulation has produced a final structure resembling that of NMR native structure to within 1A main-chain root mean square deviation. Remarkably, the simulated structure stayed in the native state for most part of the simulation after it reached the state. Of greater significance is that our simulation has not only reached the correct main-chain conformation, but also a very high degree of accuracy in side-chain packing conformation. This feat has traditionally been a challenge for ab initio simulation studies. In addition to characterization of the trajectory, comparison of our results to experimental data is also presented. Analysis of the trajectory suggests that the rate-limiting step of folding of this mini-protein is the packing of the Trp side-chain.     

Side chain accessibility and dynamics in the molten globule state of alpha-lactalbumin: a (19)F-NMR study

The molten globule state of alpha-lactalbumin (alpha-LA) has been considered a prototype of partially folded proteins. Despite the importance of molten globules in understanding the mechanisms of protein folding and its relevance to some biological phenomena, site-specific information on the structure and dynamics of a molten globule is limited, largely because of the high conformational flexibility and heterogeneity. Here, we use selective isotope labeling and (19)F NMR to investigate the solvent accessibility and side-chain dynamics of aromatic residues in the molten globule of alpha-LA. Comparison of these properties with those of the native and unfolded protein indicates that the alpha-LA molten globule is highly heterogeneous; each residue has its unique solvent accessibility and motional environment. Many aromatic residues normally buried in the interior of native alpha-LA remain significantly buried in the molten globule and the side-chain dynamics of these residues are highly restricted. Our results suggest that hydrophobic and van der Waals interactions mediated by the inaccessible surface area could be sufficient to account for all the stability of the alpha-LA molten globule, which is approximately 50% of the value for the native protein.     

Steady-state and time-resolved fluorescence studies of the intestinal fatty acid binding protein

The intestinal fatty acid binding protein contains two tryptophan residues (Trp6 and Trp82) both of which have been shown by X-ray and NMR methods to be buried in hydrophobic clusters. By using a combination of steady-state and time-resolved fluorescence experiments, we have deconvoluted the lifetime weighted contribution of each of the tryptophans to the steady-state fluorescence quantum yield. While Trp82 has been implicated in an intermediate that appears at relatively high denaturant concentrations, the variation of the lifetime weighted contribution of Trp6 with urea or guanidium hydrochloride shows formation of an intermediate state at low concentrations of the denaturant before the actual unfolding starts. Trp82 did not show similar behavior. Fluorescence quenching experiments by acrylamide show that while Trp6 in the native protein is less solvent-exposed, its accessibility is increased significantly at low urea concentration indicating that the early intermediate state is partially unfolded. Time-resolved anisotropy experiments indicate that the volume of the partially unfolded intermediates is larger than the native protein and lead to the speculation that the last step of the protein folding might be the removal of solvent molecules from the protein.     

The slow folding reaction of barstar: the core tryptophan region attains tight packing before substantial secondary and tertiary structure formation and final compaction of the polypeptide chain

The slow folding of a single tryptophan-containing mutant of barstar has been studied in the presence of 2 M urea at 10 degrees C, using steady state and time-resolved fluorescence methods and far and near-UV CD measurements. The protein folds in two major phases: a fast phase, which is lost in the dead time of measurement during which the polypeptide collapses to a compact form, is followed by a slow observable phase. During the fast phase, the rotational correlation time of Trp53 increases from 2.2 ns to 7.2 ns, and its mean fluorescence lifetime increases from 2.3 ns to 3.4 ns. The fractional changes in steady-state fluorescence, far-UV CD, and near-UV CD signals, which are associated with the fast phase are, respectively, 36 %, 46 %, and 16 %. The product of the fast phase can bind the hydrophobic dye ANS. These observations together suggest that the folding intermediate accumulated at the end of the fast phase has: (a) about 20 % of the native-state secondary structure, (b) marginally formed or disordered tertiary structure, (c) a water-intruded and mobile protein interior; and (d) solvent-accessible patches of hydrophobic groups. Measurements of the anisotropy decay of Trp53 suggest that it undergoes two types of rotational motion in the intermediate: (i) fast (tau(r) approximately 1 ns) local motion of its indole side-chain, and (ii) a slower (tau(r) approximately 7.2 ns) motion corresponding to global tumbling of the entire protein molecule. The ability of the Trp53 side-chain to undergo fast local motion in the intermediate, but not in the fully folded protein where it is completely buried in the hydrophobic core, suggests that the core of the intermediate is still poorly packed. The global tumbling time of the fully folded protein is faster at 5.6 ns, suggesting that the volume of the intermediate is 25 % more than that of the fully folded protein. The rate of folding of this intermediate to the native state, measured by steady-state fluorescence, far-UV CD, and near-UV CD, is 0.07(+/-0.01) min(-1) This rate compares to a rate of folding of 0.03(+/-0.005) min(-1), determined by double-jump experiments which monitor directly formation of native protein; and to a rate of folding of 0.05 min(-1), when determined from time-resolved anisotropy measurements of the long rotational correlation time, which relaxes from an initial value of 7.2 ns to a final value of 5. 6 ns as the protein folds. On the other hand, the amplitude of the short correlation time decreases rapidly with a rate of 0.24(+/-0.06) min(-1). These results suggest that tight packing of residues in the hydrophobic core occurs relatively early during the observable slow folding reaction, before substantial secondary and tertiary structure formation and before final compaction of the protein.     

Equilibrium unfolding of mutant apomyoglobins with substitutions of conserved nonfunctional residues by alanine

The problems of protein aggregation and protein misfolding in the cell are connected with the appearance of many genetic diseases. Both processes can be a consequence of substitutions of certain amino acid residues in proteins. The substitutions can influence the protein stability and protein folding rates in both the intermediate and the native states. We have studied equilibrium urea unfolding of mutant forms of apomyoglobin with substitutions of conserved nonfunctional residues by Ala to estimate their influence on protein stability. These residues include Val10, Trp14, Ilel11, Leu115, Met131 and Leu135. Conformational transitions were monitored by intrinsic Trp fluorescence and by circular dichroism spectra in the far UV region. Free energy changes upon the transition from the native to intermediate state and from the intermediate to unfolded state were determined. It was shown that all substitutions used lead to an appreciable decrease of the apomyoglobin native state stability, whereas the stability of the intermediate state is affected substantially smaller.     

Water dynamics clue to key residues in protein folding

A computational method independent of experimental protein structure information is proposed to recognize key residues in protein folding, from the study of hydration water dynamics. Based on all-atom molecular dynamics simulation, two key residues are recognized with distinct water dynamical behavior in a folding process of the Trp-cage protein. The identified key residues are shown to play an essential role in both 3D structure and hydrophobic-induced collapse. With observations on hydration water dynamics around key residues, a dynamical pathway of folding can be interpreted.     

Unfolding transition state and intermediates of the tumor suppressor p16INK4a investigated by molecular dynamics simulations

The ankyrin repeat is one of the most common protein motifs and is involved in protein-protein interactions. It consists of 33 residues that assume a beta-hairpin helix-loop-helix fold. Mutagenesis and kinetic experiments (Phi-value analysis of the folding transition state) have shown that the tumor suppressor p16(INK4a), a four-repeat protein, unfolds sequentially starting from the two N-terminal repeats. Here, the flexibility of p16(INK4a) at room temperature and its unfolding mechanism at high temperature have been investigated by multiple molecular dynamics runs in explicit water for a total simulation time of 0.65 micros. The transition state ensemble (TSE) of p16(INK4a) was identified by monitoring both the deviation from the experimental Phi values and sudden conformational changes along the unfolding trajectories. Conformations in the TSE have a mainly unstructured second repeat whereas the other repeats are almost completely folded. A rigid-body displacement of the first repeat involving both a rotation and translation is observed in all molecular dynamics simulations at high temperature. The Trp(15), Pro(75), and Ala(76) side-chains are more buried in the TSE than the native state. The sequential unfolding starting at the second repeat is in agreement with the mutagenesis studies whereas the displacement of the first repeat and the presence of nonnative interactions at the TSE are simulation results which supplement the experimental data. Furthermore, the unfolding trajectories reveal the presence of two on-pathway intermediates with partial alpha-helical structure. Finally, on the basis of the available experimental and simulation results we suggest that in modular proteins the shift of the folding TSE toward the native structure upon reduction of the number of tandem repeats is consistent with the Hammond effect.     

Refolding simulations of an isolated fragment of barnase into a native-like β hairpin: Evidence for compactness and hydrogen bonding as concurrent stabilizing factors

Experimental evidence and theoretical models both suggest that protein folding is initiated within specific fragments intermittently adopting conformations close to that found in the protein native structure. These folding initiation sites encompassing short portions of the protein are ideally suited for study in isolation by computational methods aimed at peering into the very early events of folding. We have used Molecular Dynamics (MD) technique to investigate the behavior of an isolated protein fragment formed by residues 85 to 102 of barnase that folds into a β hairpin in the protein native structure. Three independent MD simulations of 1.3 to 1.8 ns starting from unfolded conformations of the peptide portrayed with an all-atom model in water were carried out at gradually decreasing temperature. A detailed analysis of the conformational preferences adopted by this peptide in the course of the simulations is presented. Two of the unfolded peptide conformations fold into a hairpin characterized by native and a larger bulk of nonnative interactions. Both refolding simulations substantiate the close relationship between interstrand compactness and hydrogen bonding network involving backbone atoms. Persistent compactness witnessed by side-chain interactions always occurs concomitantly with the formation of backbone hydrogen bonds. No highly populated conformations generated in a third simulation starting from the remotest unfolded conformer relative to the native structure are observed. However, nonnative long-range and medium-range contacts with the aromatic moiety of Trp94 are spotted, which are in fair agreement with a former nuclear magnetic resonance study of a denaturing solution of an isolated barnase fragment encompassing the β hairpin. All this lends reason to believe that the 85–102 barnase fragment is a strong initiation site for folding. Proteins 29:212–227, 1997. © 1997 Wiley-Liss, Inc.     

Folding mechanism of Pseudomonas aeruginosa cytochrome c551: role of electrostatic interactions on the hydrophobic collapse and transition state properties.

We report on the folding kinetics of the small 82 residue cytochrome c551from Pseudomonas aeruginosa. The presence of two Trp residues (Trp56 and Trp77) allows the monitoring of fluorescence quenching on refolding in two different regions of the protein. A single His residue (the iron-coordinating His16) permits the study of refolding in the absence of miscoordination events. After identification of the kinetic traps (Pro isomerization and aggregation of denatured protein), overall refolding kinetics is described by two processes: (i) a burstphase collapse (faster than milliseconds) which we show to be a global event leading to a state whose compactness depends on the overall net charge; at the isoeletric pH (4.7), it is maximally compact, while above and below it is more expanded; and (ii) an exponential phase (in the millisecond time range) leading to the native protein via a transition state(s) possibly involving the formation of a specific salt bridge between Lys10 and Glu70, at the contact between the N and C-terminal helices. Comparison with the widely studied horse cytochrome c allows the discussion of similarities and differences in the folding of two proteins which have the same "fold" despite a very low degree of sequence homology (<30 %).     

The folding thermodynamics and kinetics of crambin using an all-atom Monte Carlo simulation

We present a novel Monte Carlo simulation of protein folding, in which all heavy atoms are represented as interacting hard spheres. This model includes all degrees of freedom relevant to folding, all side-chain and backbone torsions, and uses a Go potential. In this study, we focus on the 46 residue alpha/beta protein crambin and two of its structural components, the helix and helix hairpin. For a wide range of temperatures, we recorded multiple folding events of these three structures from random coils to native conformations that differ by less than 1 A C(alpha) dRMS from their crystal structure coordinates. The thermodynamics and kinetic mechanism of the helix-coil transition obtained from our simulation shows excellent agreement with currently available experimental and molecular dynamics data. Based on insights obtained from folding its smaller structural components, a possible folding mechanism for crambin is proposed. We observed that the folding occurs via a cooperative, first order-like process, and that many folding pathways to the native state exist. One particular sequence of events constitutes a "fast-folding" pathway where kinetic traps are avoided. At very low temperatures, a kinetic trap arising from the incorrect packing of side-chains was observed. These results demonstrate that folding to the native state can be observed in a reasonable amount of time on desktop computers even when an all-atom representation is used, provided the energetics sufficiently stabilize the native state.     

Folding of Desulfovibrio desulfuricans flavodoxin is accelerated by cofactor fly-casting

Folding of cofactor-binding proteins involves ligand binding in addition to polypeptide folding. We here assess the kinetic folding/binding landscape for Desulfovibrio desulfuricans flavodoxin that coordinates an FMN cofactor. The apo-form folds in a two-step process involving a burst-phase intermediate. Studies on Tyr98Ala and Trp60Ala variants reveal that these aromatics-that stack with the FMN in the holo-form-are not participating in the apo-protein folding pathway. However, these residues are essential for FMN interactions with the unfolded protein during refolding of holo-flavodoxin. Unfolding of wild-type holo-flavodoxin is coupled to FMN dissociation whereas for Tyr98Ala and Trp60Ala holo-variants, FMN dissociates before polypeptide unfolding. Both variants refold as apo-proteins before FMN rebinds. In sharp contrast, refolding of unfolded wild-type holo-flavodoxin is over an order of magnitude faster than that of the apo-form, the pathway does not include a burst-phase intermediate, and the speed is independent of FMN excess ratio. These observations demonstrate that FMN binds rapidly to the unfolded polypeptide and guides folding straight to the native state. As this path to functional D. desulfuricans holo-flavodoxin is faster than if the cofactor binds to pre-folded apo-protein, this is one of few examples where molecular recognition via a "fly-casting" mechanism is kinetically favored.     

beta-hairpin stability and folding: molecular dynamics studies of the first beta-hairpin of tendamistat

The stability and (un)folding of the 19-residue peptide, SCVTLYQSWRYSQADNGCA, corresponding to the first beta-hairpin (residues 10 to 28) of the alpha-amylase inhibitor tendamistat (PDB entry 3AIT) has been studied by molecular dynamics simulations in explicit water under periodic boundary conditions at several temperatures (300 K, 360 K and 400 K), starting from various conformations for simulation lengths, ranging from 10 to 30 ns. Comparison of trajectories of the reduced and oxidized native peptides reveals the importance of the disulphide bridge closing the beta-hairpin in maintaining a proper turn conformation, thereby insuring a proper side-chain arrangement of the conserved turn residues. This allows rationalization of the conservation of those cysteine residues among the family of alpha-amylase inhibitors. High temperature simulations starting from widely different initial configurations (native beta-hairpin, alpha and left-handed helical and extended conformations) begin sampling similar regions of the conformational space within tens of nanoseconds, and both native and non-native beta-hairpin conformations are recovered. Transitions between conformational clusters are accompanied by an increase in energy fluctuations, which is consistent with the increase in heat capacity measured experimentally upon protein folding. The folding events observed in the various simulations support a model for beta-hairpin formation in which the turn is formed first, followed by hydrogen bond formation closing the hairpin, and subsequent stabilization by side-chain hydrophobic interactions.     

Heterogeneous Packing in the Folding/Unfolding Intermediate State of Bitter Gourd Trypsin Inhibitor

The conformation and dynamics of a protein are essential in characterizing the protein folding/unfolding intermediate state. They are closely involved in the packing and site-specific interactions of peptide elements to build and stabilize the tertiary structure of the protein. In this study, it was confirmed that trypsin inhibitor obtained from seeds of bitter gourd (BGTI) adopted a peculiar but plausible conformation and dynamics in the unfolding intermediate state. The fluorescence spectrum of one of two tryptophan residues of BGTI, Trp9, shifted to the blue side in the presence of 2–3 M guanidine hydrochloride, although the other, Trp54, did not show this spectral shift. At the same time, the motional freedom of Trp9 revealed by a time-resolved fluorescence study decreased, suggesting that the segmental motion of this residue was more restricted. These results indicate that BGTI takes such a conformation state that the hydrophobic core and loop domains arranging Trp9 and Trp54 respectively are heterogeneously packed in the unfolding intermediate state.     

Influence of Trp mutation on native, intermediate, and transition states of goat alpha-lactalbumin: an equilibrium and kinetic study

Equilibrium circular dichroism and kinetic stopped-flow fluorescence studies on the stability and the folding kinetics of a set of Trp to Phe mutants of goat alpha-lactalbumin (GLA) were used to characterize the native, intermediate, and transition states of these constructs. GLA contains four tryptophan residues, three of which, Trp26, Trp104, and Trp118, are located in the alpha-domain, while the fourth, Trp60, is located in the beta-domain. Trp26, Trp60, and Trp104 are part of a hydrophobic cluster, whereas Trp118 is situated in a more flexible region near the C-terminal end of the protein. In each case, the mutation leads to a reduction in the overall stability, but only for W26F and W60F is an equilibrium intermediate observed in guanidine hydrochloride-induced unfolding experiments. In kinetic refolding experiments, however, for all samples a burst phase is observed, the amplitude of which depends on the specific mutation. Refolding and unfolding kinetics can adequately be described by a sequential three-state mechanism. phi value analysis showed that the local structure around Trp26, Trp60, and Trp104 is formed in the intermediate and in the transition state of the folding reaction, while around Trp118 no persistent native contacts are observed. From these findings, we conclude that, although hydrophobicity is a major driving force for folding, minor steric changes induced by point mutation can considerably influence the overall stability and the folding process of the protein.     

Intrinsic tryptophans of CRABPI as probes of structure and folding.

The native state fluorescence and CD spectra of the predominantly beta-sheet cellular retinoic acid-binding protein I (CRABPI) include contributions from its three tryptophan residues and are influenced by the positions of these residues in the three-dimensional structure. Using a combination of spectroscopic approaches and single Trp-mutants of CRABPI, we have deconvoluted these spectra and uncovered several features that have aided in our analysis of the development of structure in the folding pathway of CRABPI. The emission spectrum of native CRABPI is dominated by Trp 7. Trp 109 is fluorescence-silent due to its interaction with the guanidino group of Arg 111. Although the far-UV CD spectrum of CRABPI is largely determined by the protein's secondary structure, aromatic clustering around Trp 87 and the aromatic-charge interaction between Arg 111 and Trp 109 give rise to a characteristic feature in the CD spectrum at 228 nm. The near-UV CD bands of CRABPI arise largely from additive contributions of the three tryptophan residues. Trp 7 and Trp 87 give a negative CD band at 275 nm. The near-UV CD band from Trp 109 is positive and shifted to longer wavelengths (to 302 nm) due to the charge-aromatic interaction between Arg 111 and Trp 109. Our deconvolution of the equilibrium spectra have been used to interpret kinetic folding experiments monitored by stopped-flow fluorescence. These dynamic experiments suggest the early evolution of a well-populated, hydrophobically collapsed intermediate, which undergoes global rearrangement to form the fully folded structure. The results presented here suggest several additional strategies for dissecting the folding pathway of CRABPI.     

Folding simulations of Trp‐cage mini protein in explicit solvent using biasing potential replica‐exchange molecular dynamics simulations

Replica exchange molecular dynamics (RexMD) simulations are frequently used for studying structure formation and dynamics of peptides and proteins. A significant drawback of standard temperature RexMD is, however, the rapid increase of the replica number with increasing system size to cover a desired temperature range. A recently developed Hamiltonian RexMD method has been used to study folding of the Trp‐cage protein. It employs a biasing potential that lowers the backbone dihedral barriers and promotes peptide backbone transitions along the replica coordinate. In two independent applications of the biasing potential RexMD method including explicit solvent and starting from a completely unfolded structure the formation of near‐native conformations was observed after 30–40 ns simulation time. The conformation representing the most populated cluster at the final simulation stage had a backbone root mean square deviation of ～1.3 Å from the experimental structure. This was achieved with a very modest number of five replicas making it well suited for peptide and protein folding and refinement studies including explicit solvent. In contrast, during five independent continuous 70 ns molecular dynamics simulations formation of collapsed states but no near native structure formation was observed. The simulations predict a largely collapsed state with a significant helical propensity for the helical domain of the Trp‐cage protein already in the unfolded state. Hydrogen bonded bridging water molecules were identified that could play an active role by stabilizing the arrangement of the helical domain with respect to the rest of the chain already in intermediate states of the protein. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

The role of carbon-donor hydrogen bonds in stabilizing tryptophan conformations

Although most side-chain torsion angles correspond to low-energy rotameric positions, deviations occur with significant frequency. One striking example arises in Trp residues, which have an important role in stabilizing protein structures because of their size and mixed hydrophobic/hydrophilic character. Ten percent of Trp side-chains have unexplained conformations with chi(2) near 0 degrees instead of the expected 90 degrees. The current work is a structural and energetic analysis of these conformations. It is shown that many Trp residues with these orientations are stabilized by three-center carbon-donor hydrogen bonds of the form C-H...X...H-C, where X is a polar hydrogen-bond acceptor in the environment of the side-chain. The bridging hydrogen bonds occur both within the Trp side-chain and between the side-chain and the local protein backbone. Free energy maps of an isolated Trp residue in an explicit water environment show a minimum corresponding to the off-rotamer peak observed in the crystallographic data. Bridging carbon-donor hydrogen bonds are also shown to stabilize on-rotamer Trp conformations, and similar bridging hydrogen bonds also stabilize some off-rotamer Asp conformations. The present results suggest a previously unrecognized role for three-center carbon-donor hydrogen bonds in protein structures and support the view that the off-rotamer Trp side-chain orientations are real rather than artifacts of crystallographic refinements. Certain of the off-rotamer Trp conformations appear to have a functional role.