Folding and association of the antibody domain CH3: prolyl isomerization preceeds dimerization.

The simplest naturally occurring model system for studying immunoglobulin folding and assembly is the non-covalent homodimer formed by the C-terminal domains (CH3) of the heavy chains of IgG. Here, we describe the structure of recombinant CH3 dimer as determined by X-ray crystallography and an analysis of the folding pathway of this protein.Under conditions where prolyl isomerization does not contribute to the folding kinetics, formation of the beta-sandwich structure is the rate-limiting step. beta-Sheet formation of CH3 is a slow process, even compared to other antibody domains, while the subsequent association of the folded monomers is fast. After long-time denaturation, the majority of the unfolded CH3 molecules reaches the native state in two serial reactions, involving the re-isomerization of the Pro35-peptide bond to the cis configuration. The species with the wrong isomer accumulate as a monomeric intermediate. Importantly, the isomerization to the correct cis configuration is the prerequisite for dimerization of the CH3 domain. In contrast, in the Fab fragment of the same antibody, prolyl isomerization occurs after dimerization demonstrating that within one protein, comprised of highly homologous domains, both the kinetics of beta-sandwich formation and the stage at which prolyl isomerization occurs during the folding process can be completely different.     

Reversible and irreversible unfolding of multi-domain proteins

In contrast to single-domain proteins unfolding of larger multi-domain proteins is often irreversible. In a comparative case study on three different multi-domain proteins (phosphoglycerate kinase: PGK and two homologous alpha-amylases: TAKA and BLA) we investigated properties of unfolded states and their ability to fold back into the native state. For this purpose guanidine hydrochloride, alkaline pH, and thermal unfolded states were characterized. Structural alterations upon unfolding and refolding transitions were monitored using fluorescence and CD spectroscopy. Static and dynamic light scattering was employed to follow aggregation processes. Furthermore, proper refolding was also investigated by enzyme activity measurements. While for PGK at least partial reversible unfolding transitions were observed in most cases, we found reversible unfolding for TAKA in the case of alkaline pH and GndHCl induced unfolding. BLA exhibits reversible unfolding only under conditions with high concentrations of protecting osmolytes (glycerol), indicating that aggregation of the unfolded state is the main obstacle to achieve proper refolding for this protein. Structural properties, such as number and size of domains, secondary structure contents and compositions within domains, and domain topology were analyzed and considered in the interpretation of differences in refolding behavior of the investigated proteins.     

Folding mechanism of the CH2 antibody domain

The immunoglobulin C(H)2 domain is a simple model system suitable for the study of the folding of all-beta-proteins. Its structure consists of two beta-sheets forming a greek-key beta-barrel, which is stabilized by an internal disulfide bridge located in the hydrophobic core. Crystal structures of various antibodies suggest that the C(H)2 domains of the two heavy chains interact with their sugar moieties and form a homodimer. Here, we show that the isolated, unglycosylated C(H)2 domain is a monomeric protein. Equilibrium unfolding was a two-state process, and the conformational stability is remarkably low compared to other antibody domains. Folding kinetics of C(H)2 were found to consist of several phases. The reactions could be mapped to three parallel pathways, two of which are generated by prolyl isomerizations in the unfolded state. The slowest folding reaction, which was observed only after long-term denaturation, could be catalyzed by a prolyl isomerase. The majority of the unfolded molecules, however, folded more rapidly, on a time-scale of minutes. Presumably, these molecules also have to undergo prolyl isomerization before reaching the native state. In addition, we detected a small number of fast-folding molecules in which all proline residues appear to be in the correct conformation. On both prolyl isomerization limited pathways, the formation of partly structured intermediates could be observed.     

Stabilization of an immunoglobulin fold domain by an engineered disulfide bond at the buried hydrophobic region

We report for the first time the stabilization of an immunoglobulin fold domain by an engineered disulfide bond. In the llama single-domain antibody, which has human chorionic gonadotropin as its specific antigen, Ala49 and Ile70 are buried in the structure. A mutant with an artificial disulfide bond at this position showed a 10 degrees C higher midpoint temperature of thermal unfolding than that without the extra disulfide bond. The modified domains exhibited an antigen binding affinity comparable with that of the wild-type domain. Ala49 and Ile70 are conserved in camel and llama single-domain antibody frameworks. Therefore, domains against different antigens are expected to be stabilized by the engineered disulfide bond examined here. In addition to the effect of the loop constraints in the unfolded state, thermodynamic analysis indicated that internal interaction and hydration also control the stability of domains with disulfide bonds. The change in physical properties resulting from mutation often causes unpredictable and destabilizing effects on these interactions. The introduction of a hydrophobic cystine into the hydrophobic region maintains the hydrophobicity of the protein and is expected to minimize the unfavorable mutational effects.     

The thermal stability of immunoglobulin: unfolding and aggregation of a multi-domain protein

The denaturation of immunoglobulin G was studied by different calorimetric methods and circular dichroism spectroscopy. The thermogram of the immunoglobulin showed two main transitions that are a superimposition of distinct denaturation steps. It was shown that the two transitions have different sensitivities to changes in temperature and pH. The two peaks represent the F(ab) and F(c) fragments of the IgG molecule. The F(ab) fragment is most sensitive to heat treatment, whereas the F(c) fragment is most sensitive to decreasing pH. The transitions were independent, and the unfolding was immediately followed by an irreversible aggregation step. Below the unfolding temperature, the unfolding is the rate-determining step in the overall denaturation process. At higher temperatures where a relatively high concentration of (partially) unfolded IgG molecules is present, the rate of aggregation is so fast that IgG molecules become locked in aggregates before they are completely denatured. Furthermore, the structure of the aggregates formed depends on the denaturation method. The circular dichroism spectrum of the IgG is also strongly affected by both heat treatment and low pH treatment. It was shown that a strong correlation exists between the denaturation transitions as observed by calorimetry and the changes in secondary structure derived from circular dichroism. After both heat- and low-pH-induced denaturation, a significant fraction of the secondary structure remains.     

Folding and aggregation of a multi-domain engineered immunotoxin

Engineered immunotoxins with specific targeting mechanisms have potential applications for the treatment of cancer and other diseases; however, their folding behavior is often poorly understood and this presents challenges during process development, manufacturing, and formulation. Folding thermodynamics of an antibody variable domain (V_H/V_L) genetically fused to a biological toxin payload were characterized at pH 6.0 and pH 8.0 in order to assess the relative domain stabilities, along with time scales on which they fold, and the competition between aggregation and folding. The toxin and V_H/V_L domains had considerably different unfolding free energies (ΔG_UNF), leading to a thermodynamically-distinct intermediate species, with the toxin domain unfolded and the V_H/V_L folded. The intermediate is the majority species over a range of denaturant concentrations (∼4–6 M urea; ∼2–4 M guanidine HCl). Thermal unfolding resulted in reversible unfolding of the toxin domain at pH 8, but at pH 6 thermal unfolding was convoluted with aggregation due to irreversible unfolding and aggregation for the V_H/V_L domain. Chemical unfolding of both domains was more easily reversible, provided that the refold was done stepwise, allowing the antibody domain to fold first at intermediate denaturant concentration, as folding of the V_H/V_L domain played a key role in aggregation of this antibody fusion protein.     

Cold denaturation of alpha-lactalbumin

We have investigated the thermal unfolding of bovine alpha-lactalbumin by means of circular dichroism spectroscopy in the far- and near-ultraviolet regions, and shown that the native alpha-lactalbumin undergoes heat and cold denaturation. The guanidine hydrochloride-induced unfolding of alpha-lactalbumin was also investigated by circular dichroism spectroscopy at various temperatures from 261 to 318 K. It is shown that the population of the molten globule state is strongly dependent on temperature and that the molten globule state does not accumulate during the guanidine hydrochloride-induced unfolding transition at 261 K. Our results indicate that the molten globule state of alpha-lactalbumin undergoes cold denaturation as the native alpha-lactalbumin does, and that the heat capacity change of unfolding from the molten globule to the unfolded state is positive and significant. The present results further support the idea that the molten globule and the unfolded states do not belong to the same thermodynamic state, and that the native, molten globule and unfolded states are sufficient for interpreting the guanidine hydrochloride-induced unfolding behavior of alpha-lactalbumin.     

Linked thermal and solute perturbation analysis of cooperative domain interactions in proteins. Structural stability of diphtheria toxin

The temperature and guanidine hydrochloride (GuHCl) dependence of the structural stability of diphtheria toxin has been investigated by high-sensitivity differential scanning calorimetry. In 50 mM phosphate buffer at pH 8.0 and in the absence of GuHCl, the thermal unfolding of diphtheria toxin is characterized by a transition temperature (Tm) of 54.9 degrees C, a calorimetric enthalpy change (delta H) of 295 kcal/mol, and a van't Hoff to calorimetric enthalpy ratio of 0.57. Increasing the GuHCl concentration lowers the transition temperature and the calorimetric enthalpy change. At the same time, the van't Hoff to calorimetric enthalpy ratio increases until it reaches a value of 1 at 0.3 M GuHCl and remains constant thereafter. At low GuHCl concentrations (0-0.3 M), the thermal unfolding of diphtheria toxin is characterized by the presence of two transitions corresponding to the A and B domains of the protein. At higher GuHCl concentrations (0.3-1 M), the A domain is unfolded at all temperatures, and only one transition corresponding to the B domain is observed. Under these conditions, the most stable protein conformation at low temperatures is a partially folded state in which the A domain is unfolded and the B domain folded. A general model that explicitly considers the energetics of domain interactions has been developed in order to account for the stability and cooperative behavior of diphtheria toxin. It is shown that this cooperative domain interaction model correctly accounts for the temperature location as well as the shape and area of the calorimetric curves. Under physiological conditions, domain-domain interactions account for most of the structural stability of the A domain.(ABSTRACT TRUNCATED AT 250 WORDS)     

Kinetic stability plays a dominant role in the denaturant-induced unfolding of Erythrina indica lectin

The urea-induced denaturation of dimeric Erythrina indica lectin (EIL) has been studied at pH 7.2 under equilibrium and kinetic conditions in the temperature range of 40-55 degrees C. The structure of EIL is largely unaffected in this temperature range in absence of denaturant, and also in 8 M urea after incubation for 24 h at ambient temperature. The equilibrium denaturation of EIL exhibits a monophasic unfolding transition from the native dimer to the unfolded monomer as monitored by fluorescence, far-UV CD, and size-exclusion FPLC. The thermodynamic parameters determined for the two-state unfolding equilibrium show that the free energy of unfolding (DeltaGu, aq) remains practically same between 40 and 55 degrees C, with a value of 11.8 +/- 0.6 kcal mol(-1) (monomer units). The unfolding kinetics of EIL describes a single exponential decay pattern, and the apparent rate constants determined at different temperatures indicate that the rate of the unfolding reaction increases several fold with increase in temperature. The presence of probe like external metal ions (Mn2+, Ca2+) does not influence the unfolding reaction thermodynamically or kinetically; however, the presence of EDTA affects only kinetics. The present results suggest that the ability of EIL to preserve the structural integrity against the highly denaturing conditions is linked primarily to its kinetic stability, and the synergic action of heat and denaturant is involved in the unfolding of the protein.     

Exploring the temperature-pressure phase diagram of staphylococcal nuclease

The temperature dependence of the pressure-induced equilibrium unfolding of staphylococcal nuclease (Snase) was determined by fluorescence of the single tryptophan residue, FTIR absorption for the amide I' and tyrosine O-H bands, and small-angle X-ray scattering (SAXS). The results from these three techniques were similar, although the stability as measured by fluorescence was slightly lower than that measured by FTIR and SAXS. The resulting phase diagram exhibits the well-known curvature for heat and cold denaturation of proteins, due to the large decrease in heat capacity upon folding. The volume change for unfolding became less negative with increasing temperatures, consistent with a larger thermal expansivity for the unfolded state than for the folded state. Fluorescence-detected pressure-jump kinetics measurements revealed that the curvature in the phase diagram is due primarily to the rate constant for folding, indicating a loss in heat capacity for the transition state relative to the unfolded state. The similar temperature dependence of the equilibrium and activation volume changes for folding indicates that the thermal expansivities of the folded and transition states are similar. This, along with the fact that the activation volume for folding is positive over the temperature range examined, the nonlinear dependence of the folding rate constant upon temperature implicates significant dehydration in the rate-limiting step for folding of Snase.     

Thermodynamic analysis of the folding of the streptococcal protein G IgG-binding domains B1 and B2: why small proteins tend to have high denaturation temperatures.

We have cloned, expressed, and characterized two naturally occurring variations of the IgG-binding domain of streptococcal protein G. The domain is a stable cooperative folding unit of 56 amino acids, which maintains a unique folded structure without disulfide cross-links or tight ligand binding. We have studied the thermodynamics of the unfolding reaction for the two versions of this domain, designated B1 and B2, which differ by six amino acids. They have denaturation temperatures of 87.5 degrees C and 79.4 degrees C, respectively at pH 5.4, as determined by differential scanning calorimetry. Thermodynamic state functions for the unfolding reaction (delta G, delta H, delta S, and delta Cp) have been determined and reveal several interesting insights into the behavior of very small proteins. First, though the B1 domain has a heat denaturation point close to 90 degrees C, it is not unusually stable at physiologically relevant temperatures (delta G = 25 kJ/mol at 37 degrees C). This behavior occurs because the stability profile (delta G vs temperature) is flat and shallow due to the small delta S and delta Cp for unfolding. Related to this point is the second observation that small changes in the free energy of unfolding of the B-domain due to mutation or change in solvent conditions lead to large shifts in the heat denaturation temperature. Third, the magnitude and relative contributions of hydrophobic vs nonhydrophobic forces (per amino acid residue) to the total free energy of folding of the B-domain are remarkably typical of other globular proteins of much larger size.     

Calorimetric study of the heat and cold denaturation of beta-lactoglobulin.

Temperature-induced changes of the states of beta-lactoglobulin have been studied calorimetrically. In the presence of a high concentration of urea this protein shows not only heat but also cold denaturation. Its heat denaturation is approximated very closely by a two-state transition, while the cold denaturation deviates considerably from the two-state transition and this deviation increases as the temperature decreases. The heat effect of cold denaturation is opposite in sign to that of heat denaturation and is noticeably larger in magnitude. This difference in magnitude is caused by the temperature-dependent negative heat effect of additional binding of urea to the polypeptide chain of the protein upon its unfolding, which decreases the positive enthalpy of heat denaturation and increases the negative enthalpy of cold denaturation. The binding of urea considerably increases the partial heat capacity of the protein, especially in the denatured state. However, when corrected for the heat capacity effect of urea binding, the partial heat capacity of the denatured protein is close in magnitude to that expected for the unfolded polypeptide chain in aqueous solution without urea but only for temperatures below 10 degrees C. At higher temperatures, the heat capacity of the denatured protein is lower than that expected for the unfolded polypeptide chain. It appears that at temperatures above 10 degrees C not all the surface of the beta-lactoglobulin polypeptide chain is exposed to the solvent, even in the presence of 6 M urea; i.e., the denatured protein is not completely unfolded and unfolds only at temperatures lower than 10 degrees C.(ABSTRACT TRUNCATED AT 250 WORDS)     

Conformational stability and domain unfolding of the Von Willebrand factor A domains

Von Willebrand factor (VWF), a multimeric multidomain glycoprotein secreted into the blood from vascular endothelial cells, initiates platelet adhesion at sites of vascular injury. This process requires the binding of platelet glycoprotein Ib-IX-V to the A1 domain of VWF monomeric subunits under fluid shear stress. The A2 domain of VWF monomers contains a proteolytic site specific for a circulating plasma VWF metalloprotease, A Disintegrin and Metalloprotease with Thrombospondin motifs, member #13 of the ADAMTS enzyme family (ADAMTS-13), that functions to reduce adhesiveness of newly released, unusually large (UL), hyperactive forms of VWF. Shear stress assists ADAMTS-13 proteolysis of ULVWF multimers allowing ADAMTS-13 cleavage of an exposed peptide bond in the A2 domain. Shear stress may induce conformational changes in VWF, and even unfold regions of VWF monomeric subunits. We used urea as a surrogate for shear to study denaturation of the individual VWF recombinant A domains, A1, A2, and A3, and the domain triplet, A1-A2-A3. Denaturation was evaluated as a function of the urea concentration, and the intrinsic thermodynamic stability of the domains against unfolding was determined. The A1 domain unfolded in a 3-state manner through a stable intermediate. Domains A2 and A3 unfolded in a 2-state manner from native to denatured. The A1-A2-A3 triple domain unfolded in a 6-state manner through four partially folded intermediates between the native and denatured states. Urea denaturation of A1-A2-A3 was characterized by two major unfolding transitions: the first corresponding to the simultaneous complete unfolding of A2 and partial unfolding of A1 to the intermediate state; and the second corresponding to the complete unfolding of A3 followed by gradual unfolding of the intermediate state of A1 at high urea concentration. The A2 domain containing the cleavage site for ADAMTS-13 was the least stable of the three domains and was the most susceptible to unfolding. The low stability of the A2 domain is likely to be important in regulating the exposure of the A2 domain cleavage site in response to shear stress, ULVWF platelet adherence, and the attachment of ADAMTS-13 to ULVWF.     

Biophysical properties of camelid V(HH) domains compared to those of human V(H)3 domains

Camelidae possess an unusual form of antibodies lacking the light chains. The variable domain of these heavy chain antibodies (V(HH)) is not paired, while the V(H) domain of all other antibodies forms a heterodimer with the variable domain of the light chain (V(L)), held together by a hydrophobic interface. Here, we analyzed the biophysical properties of four camelid V(HH) fragments (H14, AMD9, RN05, and CA05) and two human consensus V(H)3 domains with different CDR3 loops to gain insight into factors determining stability and aggregation of immunoglobulin domains. We show by denaturant-induced unfolding equilibria that the free energies of unfolding of V(HH) fragments are characterized by Delta G(N-U) values between 21.1 and 35.0 kJ/mol and thus lie in the upper range of values for V(H) fragments from murine and human antibodies. Nevertheless, the V(HH) fragments studied here did not reach the high values between 39.7 and 52.7 kJ/mol of the human consensus V(H)3 domains with which they share the highest degree of sequence similarity. Temperature-induced unfolding of the V(HH) fragments that were studied proved to be reversible, and the binding affinity after cooling was fully retained. The melting temperatures were determined to be between 60.1 and 66.7 degrees C. In contrast, the studied V(H)3 domains aggregated during temperature-induced denaturation at 63-65 degrees C. In summary, the camelid V(HH) fragments are characterized by a favorable but not unusually high stability. Their hallmark is the ability to reversibly melt without aggregation, probably mediated by the surface mutations characterizing the V(HH) domains, which allow them to regain binding activity after heat renaturation.     

Engineered Human Antibody Constant Domains with Increased Stability

The immunoglobulin (Ig) constant CH2 domain is critical for antibody effector functions. Isolated CH2 domains are promising as scaffolds for construction of libraries containing diverse binders that could also confer some effector functions. However, previous work has shown that an isolated murine CH2 domain is relatively unstable to thermally induced unfolding. To explore unfolding mechanisms of isolated human CH2 and increase its stability γ1 CH2 was cloned and a panel of cysteine mutants was constructed. Human γ1 CH2 unfolded at a higher temperature (T_m = 54.1 °C, as measured by circular dichroism) than that previously reported for a mouse CH2 (41 °C). One mutant (m01) was remarkably stable (T_m = 73.8 °C). Similar results were obtained by differential scanning calorimetry. This mutant was also significantly more stable than the wild-type CH2 against urea induced unfolding (50% unfolding at urea concentration of 6.8 m versus 4.2 m). The m01 was highly soluble and monomeric. The existence of the second disulfide bond in m01 and its correct position were demonstrated by mass spectrometry and nuclear magnetic resonance spectroscopy, respectively. The loops were on average more flexible than the framework in both CH2 and m01, and the overall secondary structure was not affected by the additional disulfide bond. These data suggest that a human CH2 domain is relatively stable to unfolding at physiological temperature, and that both CH2 and the highly stable mutant m01 are promising new scaffolds for the development of therapeutics against human diseases.Monoclonal antibodies (mAbs)² with high affinity and specificity are now well established therapeutics and invaluable tools for biological research. It appears that their use will continue to expand in both targets and disease indications. However, a fundamental problem for full-size mAbs that limits their applications is their poor penetration into tissues (e.g. solid tumors) and poor or absent binding to regions on the surface of some molecules (e.g. on the HIV envelope glycoprotein) that are accessible by molecules of smaller size. Antibody fragments, e.g. Fabs (∼60 kDa) or single chain Fv fragments (scFvs) (20∼30 kDa), are significantly smaller than full-size antibodies (∼150 kDa), and have been used as imaging reagents and candidate therapeutics. Even smaller fragments of antibodies are of great interest and advantageous for pharmaceutical applications including cancer targeting and imaging.During the last decade a large amount of work has been aimed at developing of small size binders with scaffolds based on various highly stable human and non-human molecules (1–8). A promising direction is the development of binders based on the heavy or light chain variable region of an antibody; these fragments ranging in size from 11 kDa to 15 kDa were called “domain antibodies” or “dAbs” (7, 9). A unique kind of antibodies composed only of heavy chains are naturally formed in camels, dromedaries, and llamas, and their variable regions can also recognize antigens as single domain fragments (10). Not only is the overall size of the dAbs much smaller than that of full-size antibodies but also their paratopes are concentrated over a smaller area so that the dAbs provide the capability of interacting with novel epitopes that are inaccessible to conventional antibodies or antibody fragments with paired light and heavy chain variable domains.The structure of the constant antibody domains is similar to that of the variable domains consisting of β-strands connected mostly with loops or short helices. The second domain of the α, δ, and γ heavy chain constant regions, CH2, is unique in that it exhibits very weak carbohydrate-mediated interchain protein-protein interactions in contrast to the extensive interchain interactions that occur between the other domains. The expression of murine CH2 in bacteria, which does not support glycosylation, results in a monomeric domain (11). It has been hypothesized that the CH2 domain (CH2 of IgG, IgA, and IgD, and CH3 of IgE and IgM) could be used as a scaffold and could offer additional advantages compared with those of dAbs because it contains binding sites or portions of binding sites conferring effector and stability functions (12).It was found previously that an isolated murine CH2 is relatively unstable at physiological temperature with a temperature of 50% unfolding (T_m) slightly higher than 37 °C (11). We have hypothesized that human CH2 would exhibit different stability because of significant differences in the sequence compared with its murine counterpart. Therefore, we have extensively characterized the stability of an isolated unglycosylated single CH2 domain. We found that its stability is significantly higher than the previously reported for the murine CH2. We further increased the stability of the human CH2 by engineering an additional disulfide bond between the A and G strands. One of the newly developed mutants, denoted as m01, exhibited significantly higher stability (T_m = 73.8 °C) than that of wild type CH2. We suggest that both the wild type CH2 and the newly developed mutant, m01, could be used as scaffolds for binders. These results also demonstrate for the first time that the stability of constant antibody domains can be dramatically increased by engineering of an additional disulfide bond. The increase in stability of isolated domains may result in an increase in stability of larger antibody fragments, e.g. Fc, and therefore could have implications as a general method for increasing antibody stability. Thus, it appears that further development of CH2 and its more stable variants as scaffolds could provide new opportunities for identification of potentially useful therapeutics.     

Heat,Acid and Chemically Induced Unfolding Pathways,Conformational Stability and Structure-Function Relationship in Wheat α-Amylase

Wheat α-amylase, a multi-domain protein with immense industrial applications, belongs to α+β class of proteins with native molecular mass of 32 kDa. In the present study, the pathways leading to denaturation and the relevant unfolded states of this multi-domain, robust enzyme from wheat were discerned under the influence of temperature, pH and chemical denaturants. The structural and functional aspects along with thermodynamic parameters for α-amylase unfolding were probed and analyzed using fluorescence, circular dichroism and enzyme assay methods. The enzyme exhibited remarkable stability up to 70°C with tendency to aggregate at higher temperature. Acid induced unfolding was also incomplete with respect to the structural content of the enzyme. Strong ANS binding at pH 2.0 suggested the existence of a partially unfolded intermediate state. The enzyme was structurally and functionally stable in the pH range 4.0–9.0 with 88% recovery of hydrolytic activity. Careful examination of biophysical properties of intermediate states populated in urea and GdHCl induced denaturation suggests that α-amylase unfolding undergoes irreversible and non-coincidental cooperative transitions, as opposed to previous reports of two-state unfolding. Our investigation highlights several structural features of the enzyme in relation to its catalytic activity. Since, α-amylase has been comprehensively exploited for use in a range of starch-based industries, in addition to its physiological significance in plants and animals, knowledge regarding its stability and folding aspects will promote its biotechnological applications.     

Replica exchange simulation of reversible folding/unfolding of the Trp-cage miniprotein in explicit solvent: on the structure and possible role of internal water

We simulate the folding/unfolding equilibrium of the 20-residue miniprotein Trp-cage. We use replica exchange molecular dynamics simulations of the AMBER94 atomic detail model of the protein explicitly solvated by water, starting from a completely unfolded configuration. We employ a total of 40 replicas, covering the temperature range between 280 and 538 K. Individual simulation lengths of 100 ns sum up to a total simulation time of about 4 micros. Without any bias, we observe the folding of the protein into the native state with an unfolding-transition temperature of about 440 K. The native state is characterized by a distribution of root mean square distances (RMSD) from the NMR data that peaks at 1.8A, and is as low as 0.4A. We show that equilibration times of about 40 ns are required to yield convergence. A folded configuration in the entire extended ensemble is found to have a lifetime of about 31 ns. In a clamp-like motion, the Trp-cage opens up during thermal denaturation. In line with fluorescence quenching experiments, the Trp-residue sidechain gets hydrated when the protein opens up, roughly doubling the number of water molecules in the first solvation shell. We find the helical propensity of the helical domain of Trp-cage rather well preserved even at very high temperatures. In the folded state, we can identify states with one and two buried internal water molecules interconnecting parts of the Trp-cage molecule by hydrogen bonds. The loss of hydrogen bonds of these buried water molecules in the folded state with increasing temperature is likely to destabilize the folded state at elevated temperatures.     

Folding intermediates of the prion protein stabilized by hydrostatic pressure and low temperature

Prion diseases are associated with conformational conversion of the cellular prion protein, PrPC, into a misfolded form, PrPSc. We have investigated the equilibrium unfolding of the structured domain of recombinant murine prion protein, comprising residues 121-231 (mPrP-(121-231)). The equilibrium unfolding of mPrP-(121-231) by urea monitored by intrinsic fluorescence and circular dichroism (CD) spectroscopies indicated a two-state transition, without detectable folding intermediates. The fluorescent probe 4,4'-dianilino-1,1'-binaphthyl-5,5-disulfonic acid (bis-ANS) binds to native mPrP-(121-231), indicating exposure of hydrophobic domains on the protein surface. Increasing concentrations of urea (up to 4 M) caused the release of bound bis-ANS, whereas changes in intrinsic fluorescence and CD of mPrP took place only above 4 M urea. This indicates the existence of a partially unfolded conformation of mPrP, characterized by loss of bis-ANS binding and preservation of the overall structure of the protein, stabilized at low concentrations of urea. Hydrostatic pressure and low temperatures were also used to stabilize partially folded intermediates that are not detectable in the presence of chemical denaturants. Compression of mPrP to 3.5 kbar at 25 degrees C and pH 7 caused a slight decrease in intrinsic fluorescence emission and an 8-fold increase in bis-ANS fluorescence. Lowering the temperature to -9 degrees C under pressure reversed the decrease in intrinsic fluorescence and caused a marked (approximately 40-fold) increase in bis-ANS fluorescence. The increase in bis-ANS fluorescence at low temperatures was similar to that observed for mPrP at 1 atm at pH 4. These results suggest that pressure-assisted cold denaturation of mPrP stabilizes a partially folded intermediate that is qualitatively similar to the state obtained at acidic pH. Compression of mPrP in the presence of a subdenaturing concentration of urea stabilized another partially folded intermediate, and cold denaturation under these conditions led to complete unfolding of the protein. Possible implications of the existence of such partially folded intermediates in the folding of the prion protein and in the conversion to the PrPSc conformer are discussed.     

Microsecond simulations of the folding/unfolding thermodynamics of the Trp‐cage miniprotein

We study the unbiased folding/unfolding thermodynamics of the Trp‐cage miniprotein using detailed molecular dynamics simulations of an all‐atom model of the protein in explicit solvent using the Amberff99SB force field. Replica‐exchange molecular dynamics simulations are used to sample the protein ensembles over a broad range of temperatures covering the folded and unfolded states at two densities. The obtained ensembles are shown to reach equilibrium in the 1 μs/replica timescale. The total simulation time used in the calculations exceeds 100 μs. Ensemble averages of the fraction folded, pressure, and energy differences between the folded and unfolded states as a function of temperature are used to model the free energy of the folding transition, ΔG(P, T), over the whole region of temperatures and pressures sampled in the simulations. The ΔG(P, T) diagram describes an ellipse over the range of temperatures and pressures sampled, predicting that the system can undergo pressure‐induced unfolding and cold denaturation at low temperatures and high pressures, and unfolding at low pressures and high temperatures. The calculated free energy function exhibits remarkably good agreement with the experimental folding transition temperature (T_f = 321 K), free energy, and specific heat changes. However, changes in enthalpy and entropy are significantly different than the experimental values. We speculate that these differences may be due to the simplicity of the semiempirical force field used in the simulations and that more elaborate force fields may be required to describe appropriately the thermodynamics of proteins. Proteins 2010. © 2010 Wiley‐Liss, Inc.