Molecular dynamics simulation of thionated hen egg white lysozyme

Understanding of the driving forces of protein folding is a complex challenge because different types of interactions play a varying role. To investigate the role of hydrogen bonding involving the backbone, the effect of thio substitutions in a protein, hen egg white lysozyme (HEWL), was investigated through molecular dynamics simulations of native as well as partly (only residues in loops) and fully thionated HEWL using the GROMOS 54A7 force field. The results of the three simulations show that the structural properties of fully thionated HEWL clearly differ from those of the native protein, while for partly thionated HEWL they only changed slightly compared with native HEWL. The analysis of the torsional-angle distributions and hydrogen bonds in the backbone suggests that the α-helical segments of native HEWL tend to show a propensity to convert to 3(10)-helical geometry in fully thionated HEWL. A comparison of the simulated quantities with experimental NMR data such as nuclear overhauser effect (NOE) atom-atom distance bounds and (3)J((H)(N)(H)(α))-couplings measured for native HEWL illustrates that the information content of these quantities with respect to the structural changes induced by thionation of the protein backbone is rather limited.     

Study of the compact denatured state of a protein by molecular dynamics simulation]

The native state can be considered as a unique conformation of the protein molecule with the lowest free energy of residue contacts. In this case, all other conformations correspond to the denatured state. The degree of their compactness varies significantly. Under folding conditions, the compact denatured state rather than the random coil is in equilibrium with native protein. The balance between the main forces of protein folding, the solvophobic interactions and conformational entropy, suggests that some properties of the compact denatured state are close to those of native protein, whereas other properties are close to those of the random coil. To investigate the molecular structure of the compact denatured state, the method of molecular dynamics simulation seems to be very useful.     

Cooperativity of a protein folding reaction probed at multiple chain positions by real-time 2D NMR spectroscopy

The refolding reaction of S54G/P55N ribonuclease T1 is a two-step process, where fast formation of a partly folded intermediate is followed by the slow reaction to the native state, limited by a trans --> cis isomerization of Pro39. The hydrodynamic radius of this kinetic folding intermediate was determined by real-time diffusion NMR spectroscopy. Its folding to the native state was monitored by a series of 128 very fast 2D (15)N-HMQC spectra, to observe the kinetics of 66 individual backbone amide probes. We find that the intermediate is as compact as the native protein with many native chemical shifts. All 66 analyzed amide probes follow the rate-limiting prolyl isomerization, which indicates that this cooperative refolding reaction is fully synchronized. The stability of the folding intermediate was determined from the protection factors of 45 amide protons derived from a competition between refolding and H/D exchange. The intermediate has already gained 40% of the Gibbs free energy of refolding with many protected amides in not-yet-native regions.     

Structure of hen egg-white lysozyme solvated in TFE/water: a molecular dynamics simulation study based on NMR data

Various experimental studies of hen egg white lysozyme (HEWL) in water and TFE/water clearly indicate structural differences between the native state and TFE state of HEWL, e.g. the helical content of the protein in the TFE state is much higher than in the native state. However, the available detailed NMR studies were not sufficient to determine fully a structure of HEWL in the TFE state. Different molecular dynamics (MD) simulations, i.e. at room temperature, at increased temperature and using proton–proton distance restraints derived from NMR NOE data, have been used to generate configurational ensembles corresponding to the TFE state of HEWL. The configurational ensemble obtained at room temperature using atom-atom distance restraints measured for HEWL in TFE/water solution satisfies the experimental data and has the lowest protein energy. In this ensemble residues 50–58, which are part of the β-sheet in native HEWL, adopt fluctuating α-helical secondary structure.     

Engineering a compact non-native state of intestinal fatty acid-binding protein

The last three C-terminal residues (129-131) of intestinal fatty acid-binding protein (IFABP) participate in four main-chain hydrogen bonds and two electrostatic interactions to sequentially distant backbone and side-chain atoms. To assess if these interactions are involved in the final adjustment of the tertiary structure during folding, we engineered an IFABP variant truncated at residue 128. An additional mutation, Trp-6-->Phe, was introduced to simplify the conformational analysis by optical methods. Although the changes were limited to a small region of the protein surface, they resulted in an IFABP with altered secondary and tertiary structure. Truncated IFABP retains some cooperativity, is monomeric, highly compact, and has the molecular dimensions and shape of the native protein. Our results indicated that residues 129-131 are part of a crucial conformational determinant in which several long-range interactions, essential for the acquisition of the native state, are established. This work suggests that carefully controlled truncation can populate equilibrium non-native states under physiological conditions. These non-native states hold a great promise as experimental models for protein folding.     

Structural characterization of partially folded intermediates of apomyoglobin H64F

We present a detailed investigation of unfolded and partially folded states of a mutant apomyoglobin (apoMb) where the distal histidine has been replaced by phenylalanine (H64F). Previous studies have shown that substitution of His64, located in the E helix of the native protein, stabilizes the equilibrium molten globule and native states and leads to an increase in folding rate and a change in the folding pathway. Analysis of changes in chemical shift and in backbone flexibility, detected via [1H]-15N heteronuclear nuclear Overhauser effect measurements, indicates that the phenylalanine substitution has only minor effects on the conformational ensemble in the acid- and urea-unfolded states, but has a substantial effect on the structure, dynamics, and stability of the equilibrium molten globule intermediate formed near pH 4. In H64F apomyoglobin, additional regions of the polypeptide chain are recruited into the compact core of the molten globule. Since the phenylalanine substitution has negligible effect on the unfolded ensemble, its influence on folding rate and stability comes entirely from interactions within the compact folded or partly folded states. Replacement of His64 with Phe leads to favorable hydrophobic packing between the helix E region and the molten globule core and leads to stabilization of helix E secondary structure and overall thermodynamic stabilization of the molten globule. The secondary structure of the equilibrium molten globule parallels that of the burst phase kinetic intermediate; both intermediates contain significant helical structure in regions of the polypeptide that comprise the A, B, E, G, and H helices of the fully folded protein.     

Mapping long-range contacts in a highly unfolded protein

Insights into the earliest events in protein folding can be obtained by analysis of the conformational propensities of unfolded or partly folded states. The structure of the acid-unfolded state of apomyoglobin has been characterized using paramagnetic spin labeling and NMR. Nitroxide side-chains, introduced by coupling to mutant cysteine residues at positions 18, 77, and 133, were used as probes of chain compaction and long-range tertiary contacts. Significant interactions are observed within and between the N and C termini, while the central region of the polypeptide chain behaves as a random polymer. Even in this highly denatured form, the protein samples transient compact states in which there are native-like contacts between the N and C-terminal regions.     

Use of a quantitative structure-property relationship to design larger model proteins that fold rapidly

A quantitative structure-property relationship (QSPR) was used to design model protein sequences that fold repeatedly and relatively rapidly to stable target structures. The specific model was a 125-residue heteropolymer chain subject to Monte Carlo dynamics on a simple cubic lattice. The QSPR was derived from an analysis of a database of 200 sequences by a statistical method that uses a genetic algorithm to select the sequence attributes that are most important for folding and a neural network to determine the corresponding functional dependence of folding ability on the chosen attributes. The QSPR depends on the number of anti-parallel sheet contacts, the energy gap between the native state and quasi-continuous part of the spectrum and the total energy of the contacts between surface residues. Two Monte Carlo procedures were used in series to optimize both the target structures and the sequences. We generated 20 fully optimized sequences and 60 partially optimized control sequences and tested each for its ability to fold in dynamic MC simulations. Although sequences in which either the number of anti-parallel sheet contacts or the energy of the surface residues is non-optimal are capable of folding almost as well as fully optimized ones, sequences in which only the energy gap is optimized fold markedly more slowly. Implications of the results for the design of proteins are discussed.     

The dynamical response of hen egg white lysozyme to the binding of a carbohydrate ligand

It has become clear that the binding of small and large ligands to proteins can invoke significant changes in side chain and main chain motion in the fast picosecond to nanosecond timescale. Recently, the use of a "dynamical proxy" has indicated that changes in these motions often reflect significant changes in conformational entropy. These entropic contributions are sometimes of the same order as the total entropy of binding. Thus, it is important to understand the connections amongst motion between the manifold of states accessible to the native state of proteins, the corresponding entropy, and how this impacts the overall energetics of protein function. The interaction of proteins with carbohydrate ligands is central to a range of biological functions. Here, we examine a classic carbohydrate interaction with an enzyme: the binding of wild-type hen egg white lysozyme (HEWL) to the natural, competitive inhibitor chitotriose. Using NMR relaxation experiments, backbone amide and side chain methyl axial order parameters were obtained across apo and chitotriose-bound HEWL. Upon binding, changes in the apparent amplitude of picosecond to nanosecond main chain and side chain motions are seen across the protein. Indeed, binding of chitotriose renders a large contiguous fraction of HEWL effectively completely rigid. Changes in methyl flexibility are most pronounced closest to the binding site, but average to only a small overall change in the dynamics across the protein. The corresponding change in conformational entropy is unfavorable and estimated to be a significant fraction of the total binding entropy.     

Directionality in protein fold prediction

Background  

Ever since the ground-breaking work of Anfinsen et al. in which a denatured protein was found to refold to its native state, it has been frequently stated by the protein fold prediction community that all the information required for protein folding lies in the amino acid sequence. Recent in vitro experiments and in silico computational studies, however, have shown that cotranslation may affect the folding pathway of some proteins, especially those of ancient folds. In this paper aspects of cotranslational folding have been incorporated into a protein structure prediction algorithm by adapting the Rosetta program to fold proteins as the nascent chain elongates. This makes it possible to conduct a pairwise comparison of folding accuracy, by comparing folds created sequentially from each end of the protein.     

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.     

Effect of solvent conditions upon refolding pathways and intermediates for a simple lattice protein

Folding pathways and intermediates for a two-dimensional lattice protein have been investigated via computer simulation at various denaturant concentrations. The protein is represented as a chain of 8 hydrophobic (H) and 12 polar (P) beads on a square lattice sequenced in such a way that the native state is a compact hydrophobic core surrounded by a shell of polar beads. Two nonbonded H beads are said to attract each other with a potential of mean force of strength ϵ. Increasing |ϵ/kT| mimics decreasing the denaturant concentration in the solution. Dynamic Monte Carlo simulations have been performed in order to investigate the folding transition and the folding pathways. Sharp folding—unfolding transitions are observed and the folding process proceeds along well-defined pathways that are populated by partially folded intermediates. The folding pathways as well as the populations of the intermediates are strongly dependent upon the denaturant concentration. Generally, intermediates containing long open stretches of H beads are more populated at high denaturant concentration, whereas compact intermediates containing a substantial number of hydrophobic contacts are more populated at low denaturant concentrations. The folding process is also observed to be cooperative in nature in that the chain does not start folding until a key fold in the middle section of the chain is formed correctly. © 1997 John Wiley & Sons, Inc. Biopoly 42: 399–409, 1997     

Tumour necrosis factor is a compact trimer

Recombinant produced human tumour necrosis factor (TNF) has been studied to characterise the subunit structure of the protein. TNF is shown to be a trimer Mr 52000 in which the subunits are associated in a compact, triangular form. In secondary structure it belongs to the all-beta class of proteins. It has high thermodynamic stability and the unfolded subunits can fold and associate spontaneously to form native, biologically active TNF.     

Trehalose favors a cutinase compact intermediate off-folding pathway

The folding of cutinase, an enzyme displaying lipolytic activity, has been studied in the presence of trehalose. Equilibrium unfolding data show that trehalose increases the free energy change between folded and unfolded states. Unfolding kinetics reveal the presence of an intermediate which is ca. 60% folded in terms of solvent exposure. Trehalose stabilizes this intermediate relative to the folded state. In contrast, the intermediate revealed by folding kinetics is more compact than the transition state, as shown by the positive slope observed at low denaturant concentration in the chevron plot, as well as the decrease in the observable rate constant for folding with the increase in trehalose concentration. This intermediate displays more than 50% of area buried from the solvent (relative to the native state) compared to around 40% for the transition state for folding and therefore appears to be off the folding pathway. Trehalose stabilizes and guanidine hydrochloride destabilizes this compact intermediate. Both unfolding and folding kinetics show that compact conformational states are stabilized by trehalose, in agreement with current models on the effect of compatible solutes. This effect occurs even for compact states that decelerate the folding as in the case of the intermediate revealed by folding kinetics.     

Protein topology affects the appearance of intermediates during the folding of proteins with a flavodoxin-like fold

The topology of a native protein influences the rate with which it is formed, but does topology affect the appearance of folding intermediates and their specific role in kinetic folding as well? This question is addressed by comparing the folding data recently obtained on apoflavodoxin from Azotobacter vinelandii with those available on all three other alpha-beta parallel proteins the kinetic folding mechanism of which has been studied, i.e. Anabaena apoflavodoxin, Fusarium solani pisi cutinase and CheY. Two kinetic folding intermediates, one on-pathway and the other off-pathway, seem to be present during the folding of proteins with an alpha-beta parallel, also called flavodoxin-like, topology. The on-pathway intermediate lies on a direct route from the unfolded to the native state of the protein involved. The off-pathway intermediate needs to unfold to allow the production of native protein. Available simulation data of the folding of CheY show the involvement of two intermediates with characteristics that resemble those of the two intermediates experimentally observed. Apparently, protein topology governs the appearance and kinetic roles of protein folding intermediates during the folding of proteins that have a flavodoxin-like fold.     

Correlation between changes in nuclear magnetic resonance order parameters and conformational entropy: molecular dynamics simulations of native and denatured staphylococcal nuclease

Recent work has suggested that changes in NMR order parameters may quantitatively reflect changes in the conformational entropy of a protein ensemble. The extent of the mathematical relationship between local entropy changes as seen by NMR order parameters and the full protein entropy change is a complex issue. As a step towards a fuller understanding of this problem, molecular dynamics calculations of both native and denatured staphylococcal nuclease were performed. The N-H bond vector motion, in both explicit and implicit solvent, was analyzed to estimate local and global entropy changes. The calculated N-H bond vector order parameters from simulation agreed on average with experimental values for both native and denatured structures. However, the inverted-U profile of order parameters versus residue number observed experimentally for denatured nuclease was only partially reproduced by simulation of compact denatured structures. Comparisons made across the full set of simulations revealed a correlation between the N-H order parameter-based conformational entropy change and the total quasiharmonic-based conformational entropy change between the native and denatured structures. The calculations showed that about 25% of the total entropy change was reflected by changes in simulated S2 values. This result suggests that NMR-derived order parameters may be used to provide a reasonable estimate of the total conformational entropy change on protein folding.     

Loop fold structure of proteins: resolution of Levinthas paradox

According to Levinthal a protein chain of ordinary size would require enormous time to sort its conformational states before the final fold is reached. Experimentally observed time of folding suggests an estimate of the chain length for which the time would be sufficient. This estimate by order of magnitude fits to experimentally observed universal closed loop elements of globular proteins - 25-30 residues.     

Factors that affect the folding ability of proteins

The folding ability of a heteropolymer model for proteins subject to Monte Carlo dynamics on a simple cubic lattice is shown to be strongly correlated with the stability of the native state. We consider a number of estimates of the stability that can be determined without simulation, including the energy gap between the native state and the structurally dissimilar part of the spectrum (Z score) and, for sequences with fully compact native states, the gap in energy between the native and first excited fully compact states. These estimates are found to be more robust predictors of folding ability than a parameter sigma that requires simulation for its evaluation: sigma = 1 - Tf/Ttheta, where Tf is the temperature at which the fluctuation of an order parameter is at its maximum and Ttheta is the temperature at which the specific heat is at its maximum. We show that the interpretation of Ttheta as the collapse transition temperature is not correct in general and that the correlation between sigma and the folding ability arises from the fact that sigma is related to the energy gap (Z score).     

Alpha-lactalbumin forms a compact molten globule in the absence of disulfide bonds.

Human alpha-lactalbumin (alpha-LA) is a four disulfide-bonded protein that adopts partially structured conformations under a variety of mildly denaturing conditions. At low pH, the protein is denatured but compact, with a high degree of secondary structure and a native-like fold. This is commonly referred to as a molten globule. A variant of alpha-LA, in which all eight cysteines have been mutated to alanine (all-Ala alpha-LA), has been studied using NMR spectroscopy. At low pH all-Ala alpha-LA is nearly as compact as wild type alpha-LA. Urea-induced unfolding experiments reveal that the residues that remain compact in the absence of disulfide bonds are those that are most resistant to unfolding in the wild-type alpha-LA molten globule. This is particularly remarkable because this stable core is formed by segments of the polypeptide chain from both the N- and C-termini. These results show that the overall architecture of the protein fold of alpha-LA is determined by the polypeptide sequence itself, and not as the result of cross-linking by disulfide bonds, and provide insight into the way in which the sequence codes for the fold.     

All-atom ab initio folding of a diverse set of proteins

Natural proteins fold to a unique, thermodynamically dominant state. Modeling of the folding process and prediction of the native fold of proteins are two major unsolved problems in biophysics. Here, we show successful all-atom ab initio folding of a representative diverse set of proteins by using a minimalist transferable-energy model that consists of two-body atom-atom interactions, hydrogen bonding, and a local sequence-energy term that models sequence-specific chain stiffness. Starting from a random coil, the native-like structure was observed during replica exchange Monte Carlo (REMC) simulation for most proteins regardless of their structural classes; the lowest energy structure was close to native-in the range of 2-6 A root-mean-square deviation (rmsd). Our results demonstrate that the successful folding of a protein chain to its native state is governed by only a few crucial energetic terms.