Kinetics and thermodynamics of ligand binding to a molten globular enzyme and its native counterpart

An engineered monomeric chorismate mutase (mMjCM) has been found to combine high catalytic activity with the characteristics of a molten globule. To gain insight into the dramatic structural changes that accompany binding of a transition-state analog, we examined mMjCM by isothermal calorimetry and compared it with its dimeric parent protein, MjCM (CM from Methanococcus jannaschii), a thermostable and conventionally folded enzyme. As expected for a ligand-induced ordering process, there is a large entropic penalty for binding to the monomer relative to the dimer (− TΔΔS = 5.1 ± 0.5 kcal/mol, at 20 °C). However, this unfavorable entropy term is largely offset by enthalpic gains (ΔΔH = − 3.5 ± 0.4 kcal/mol), presumably arising from tightening of non-covalent interactions throughout the monomeric complex. Stopped-flow kinetic measurements further reveal that the catalytic molten globule binds and releases ligands significantly faster than its natural counterpart, demonstrating that partial structural disorder can speed up molecular recognition. These results illustrate how structural plasticity may strongly perturb the thermodynamics and kinetics of transition-state recognition while negligibly affecting catalytic efficiency.     

Comparison of the molten globule states of thermophilic and mesophilic alpha-amylases

In recent years great interest has been generated in the process of protein folding, and the formation of intermediates during the folding process has been proven with new experimental strategies. In the present work, we have examined the molten globule state of Bacillus licheniformis alpha-amylase (BLA) by intrinsic fluorescence and circular dichroism spectra, 1-anilino naphthalene-8-sulfonate (ANS) binding and proteolytic digestion by pepsin, for comparison to its mesophilic counterpart, Bacillus amyloliquefaciens alpha-amylase (BAA). At pH 4.0, both enzymes acquire partially folded state which show characteristics of molten globule state. They unfold in such a way that their hydrophobic surfaces are exposed to a greater extent compared to the native forms. Chemical denaturation studies by guanidine hydrochloride and proteolytic digestion with pepsin show that molten globule state of BLA is more stable than from BAA. Results from gel filtration indicate that BAA has the same compactness at pH 4.0 and 7.5. However, molten globule state of BLA is less compact than its native state. The effects of polyols such as trehalose, sorbitol and glycerol on refolding of enzymes from molten globule to native state were also studied. These polyols are effective on refolding of mesophilic alpha-amylase but only slightly effect on BLA refolding. In addition, the folding pathway and stability of intermediate state of the thermophilic and the mesophilic alpha-amylases are discussed.     

Hydrophobic effect on the stability and folding of a hyperthermophilic protein

Ribonuclease HII from hyperthermophile Thermococcus kodakaraensis (Tk-RNase HII) is a kinetically robust monomeric protein. The conformational stability and folding kinetics of Tk-RNase HII were measured for nine mutant proteins in which a buried larger hydrophobic side chain is replaced by a smaller one (Leu/Ile to Ala). The mutant proteins were destabilized by 8.9 to 22.0 kJ mol^− 1 as compared with the wild-type protein. The removal of each -CH₂- group burial decreased the stability by 5.1 kJ mol^− 1 on average in the mutant proteins of Tk-RNase HII examined. This is comparable with the value of 5.3 kJ mol^− 1 obtained from experiments for proteins from organisms growing at moderate temperature. We conclude that the hydrophobic residues buried inside protein molecules contribute to the stabilization of hyperthermophilic proteins to a similar extent as proteins at normal temperature. In the folding experiments, the mutant proteins of Tk-RNase HII examined exhibited faster unfolding compared with the wild-type protein. These results indicate that the buried hydrophobic residues strongly contribute to the kinetic robustness of Tk-RNase HII. This is the first report that provides a practical cause of slow unfolding of hyperthermostable proteins.     

Increased folding stability of TEM-1 beta-lactamase by in vitro selection

In vitro selections of stabilized proteins lead to more robust enzymes and, at the same time, yield novel insights into the principles of protein stability. We employed Proside, a method of in vitro selection, to find stabilized variants of TEM-1 β-lactamase from Escherichia coli. Proside links the increased protease resistance of stabilized proteins to the infectivity of a filamentous phage. Several libraries of TEM-1 β-lactamase variants were generated by error-prone PCR, and variants with increased protease resistance were obtained by raising temperature or guanidinium chloride concentration during proteolytic selections. Despite the small size of phage libraries, several strongly stabilizing mutations could be obtained, and a manual combination of the best shifted the profiles for thermal unfolding and temperature-dependent inactivation of β-lactamase by almost 20 °C to a higher temperature. The wild-type protein unfolds in two stages: from the native state via an intermediate of the molten-globule type to the unfolded form. In the course of the selections, the native protein was stabilized by 27 kJ mol^− 1 relative to the intermediate and the cooperativity of unfolding was strongly increased. Three of our stabilizing replacements (M182T, A224V, and R275L) had been identified independently in naturally occurring β-lactamase variants with extended substrate spectrum. In these variants, they acted as global suppressors of destabilizations caused by the mutations in the active site. The comparison between the crystal structure of our best variant and the crystal structure of the wild-type protein indicates that most of the selected mutations optimize helices and their packing. The stabilization by the E147G substitution is remarkable. It removes steric strain that originates from an overly tight packing of two helices in the wild-type protein. Such unfavorable van der Waals repulsions are not easily identified in crystal structures or by computational approaches, but they strongly reduce the conformational stability of a protein.     

Trifluoroethanol stabilizes the molten globule state and induces non-amyloidic turbidity in stem bromelain near its isoelectric point

Stem bromelain (SBM) is a therapeutic protein that has been studied for alkaline denaturation in the intestines, the principal site of its absorption. In this study, we investigated fluorinated alcohol 2,2,2-trifluoroethanol (TFE)-induced conformational changes in the specific/pre-molten globule (SMG) state of SBM observed at pH 10 by spectroscopic methods. Far-UV circular dichroism (CD) spectra showed that the protein retained its native-like secondary structure at TFE concentrations of up to 30% with a pronounced minimum at 222 nm, characteristic of a helix. However, addition of slightly higher TFE concentrations (≥40%) resulted in an ∼2.5-fold induction of this helical feature and a time-dependent increase in non-amyloidic turbidity as evidenced by turbidometric, Congo red-binding, and Thioflavin T (ThT)-binding studies. Near-UV CD spectra suggested a gradual but significant loss of tertiary structure at 10-30% TFE. Tryptophan studies showed blue-shifted fluorescence, although the number of accessible tryptophans remained the same up to 30% TFE. The SMG showed enhanced binding of the fluorescent probe 1-anilino-8-naphthalene sulfonic acid (ANS) up to 30% TFE, beyond which binding plateaued. Thermal and guanidine hydrochloride (GdnHCl) transition studies in the near-UV range indicated a single cooperative transition for the SMG state in the presence of 30% TFE, similar to that observed for native SBM at pH 7.0 (although with different T_ms), unlike the SMG state. TFE (30%) appeared to induce native-like stability to the original SMG. These observations suggest a transformation of the SMG to a characteristic molten globule (MG) conformation at 30% TFE, possibly due to TFE-induced rearrangement of hydrophobic interactions at the protein's isoelectric point.     

Characterization of functional intermediates of endoglucanase from Aspergillus aculeatus during urea and guanidine hydrochloride unfolding

Low concentrations of urea and GuHCl (2 M) enhanced the activity of endoglucanase (EC 3.1.2.4) from Aspergillus aculeatus by 2.3- and 1.9-fold, respectively. The K_m values for controls, in the presence of 2 M urea and GuHCl, were found to be 2.4 ± 0.2 × 10⁻⁸ mol L⁻¹, 1.4 ± 0.2 × 10⁻⁸ mol L⁻¹, and 1.6 ± 0.2 × 10⁻⁸ mol L⁻¹, respectively. The dissociation constant (K_d) showed changes in the affinity of the enzyme for the substrate with increases in the K_cat suggesting an increased turnover number in the presence of urea and GuHCl. Fluorescence studies showed changes in the microenvironment of the protein. The increase in the activity of this intermediate state was due to conformational changes accompanied by increased flexibility at the active site.     

Electrostatic interactions in the denatured state and in the transition state for protein folding: effects of denatured state interactions on the analysis of transition state structure

The development of electrostatic interactions during the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) is investigated by pH-dependent rate equilibrium free energy relationships. We show that Asp8, among six acidic residues, is involved in non-native, electrostatic interactions with K12 in the transition state for folding as well as in the denatured state. The perturbed native state pK(a) of D8 (pK(a) = 3.0) appears to be maintained through non-native interactions in both the transition state and the denatured state. Mutational effects on the stability of the transition state for protein (un)folding are often analyzed in respect to change in ground states. Thus, the interpretation of transition state analysis critically depends on an understanding of mutational effects on both the native and denatured state. Increasing evidence for structurally biased denatured states under physiological conditions raises concerns about possible denatured state effects on folding studies. We show that the structural interpretation of transition state analysis can be altered dramatically by denatured state effects.     

Thermodynamic properties of BPTI variants with highly simplified amino acid sequences

We report the first detailed thermodynamic analysis of simplified proteins by differential scanning calorimetry (DSC). The experiments were carried out with five simplified BPTI variants, whose structures and activities have been reported, in which several residues not essential for specifying the tertiary structure were replaced by alanine. In most aspects, the thermodynamics of simplified proteins were very similar to, if not essentially identical with, those of natural proteins. In particular, they undergo a highly cooperative two-state thermal unfolding process with a large enthalpy change, which is a thermodynamic hallmark of the native state of natural globular proteins. Furthermore, the specific enthalpy and entropy changes upon unfolding at 110 degrees C were close to values invariably observed for small natural globular proteins (55 J g(-1) and ~16 J K(-1) g(-1), respectively). On the other hand, two simplified BPTI variants, BPTI-21 and BPTI-22 (containing 21 and 22 alanine residues), were enthalpically stabilized while entropically destabilized with respect to the reference BPTI-[5,55] molecule. This peculiar type of entropy-enthalpy compensation is in sharp contrast to the usual enthalpy destabilization/entropy stabilization observed in mutational studies of natural proteins. Overall, we conclude that a thermodynamic native state can be achieved by proteins encoded with extensively simplified sequences.     

The role of hERG1 K channels and a functional link between hERG1 K channels and SDF-1 in acute leukemic cell migration

Stromal cell-derived factor-1 (SDF-1) and its unique receptor, CXCR4, regulate stem/progenitor cell migration and retention in the bone marrow and are required for hematopoiesis. Recent studies found that hERG1 K⁺ channels were important regulators of tumor cell migration. In this study, we investigated whether SDF-1 induced acute leukemic cell migration associated with hERG1 K⁺ channels. Our results showed that E-4031, a specific hERG1 K⁺ channels inhibitor, significantly blocked SDF-1-induced migration of leukemic cell lines, primary acute leukemic cells, leukemic stem cells and HEK293T cells transfected with herg-pEGFP. The migration of phenotypically recognizable subsets gave the indication that lymphoblastic leukemic cells were inhibited more than myeloid cells while in the presence of E-4031 which maybe associated with herg expression. SDF-1 increased hERG1 K⁺ current expressed in oocytes and HEK293T cells transfected with herg-pEGFP. There were no significant changes of CXCR4 expression on both HL-60 cells and primary leukemic cells regardless if untreated or treated with E-4031 for 24 h (P > 0.05). The hERG1 K⁺ current increased by SDF-1 might contribute to the mechanism of SDF-1-induced leukemic cell migration. The data suggested that hERG1 K⁺ channels functionally linked to cell migration induced by SDF-1.     

Energetic coupling between native-state prolyl isomerization and conformational protein folding

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.     

Kinetic isotope effects reveal the presence of significant secondary structure in the transition state for the folding of the N-terminal domain of L9

Our present understanding of the nature of the transition state for protein folding depends predominantly on studies where individual side-chain contributions are mapped out by mutational analysis (phi value analysis). This approach, although extremely powerful, does not in general provide direct information about the formation of backbone hydrogen bonds. Here, we report the results of amide H/D isotope effect studies that probe the development of hydrogen bonded interactions in the transition state for the folding of a small alpha-beta protein, the N-terminal domain of L9. Replacement of amide protons by deuterons in a solvent of constant isotopic composition destabilized the domain, decreasing both its T(m) and Delta G(0) of unfolding. The folding rate also decreased. The parameter Phi(H/D), defined as the ratio of the effect of isotopic substitution upon the activation free energy to the equilibrium free energy was determined to be 0.6 in a D(2)O background and 0.75 in a H(2)O background, indicating that significant intraprotein hydrogen bond interactions are developed in the transition state for the folding of NTL9. The value is in remarkably good agreement with more traditional measures of the position of the transition state, which report on the relative burial of surface area. The results provide a picture of a compact folding transition state containing significant secondary structure. Indirect analysis argues that the bulk of the kinetic isotope effect arises from the beta-sheet-rich region of the protein, and suggests that the development of intraprotein hydrogen bonds in this region plays a critical role in the folding of NTL9.     

Order of steps in the cytochrome C folding pathway: evidence for a sequential stabilization mechanism

Previous work used hydrogen exchange (HX) experiments in kinetic and equilibrium modes to study the reversible unfolding and refolding of cytochrome c (Cyt c) under native conditions. Accumulated results now show that Cyt c is composed of five individually cooperative folding units, called foldons, which unfold and refold as concerted units in a stepwise pathway sequence. The first three steps of the folding pathway are linear and sequential. The ordering of the last two steps has been unclear because the fast HX of the amino acid residues in these foldons has made measurement difficult. New HX experiments done under slower exchange conditions show that the final two foldons do not unfold and refold in an obligatory sequence. They unfold separately and neither unfolding obligately contains the other, as indicated by their similar unfolding surface exposure and the specific effects of destabilizing and stabilizing mutations, pH change, and oxidation state. These results taken together support a sequential stabilization mechanism in which folding occurs in the native context with prior native-like structure serving to template the stepwise formation of subsequent native-like foldon units. Where the native structure of Cyt c requires sequential folding, in the first three steps, this is found. Where structural determination is ambiguous, in the final two steps, alternative parallel folding is found.     

Consequences of Domain Insertion on the Stability and Folding Mechanism of a Protein

SlyD, the sensitive-to-lysis protein from Escherichia coli, consists of two domains. They are not arranged successively along the protein chain, but one domain, the “insert-in-flap” (IF) domain, is inserted internally as a guest into a surface loop of the host domain, which is a prolyl isomerase of the FK506 binding protein (FKBP) type. We used SlyD as a model to elucidate how such a domain insertion affects the stability and folding mechanism of the host and the guest domain. For these studies, the two-domain protein was compared with a single-domain variant SlyDΔIF, SlyD* without the chaperone domain (residues 1-69 and 130-165) in which the IF domain was removed and replaced by a short loop, as present in human FKBP12. Equilibrium unfolding and folding kinetics followed an apparent two-state mechanism in the absence and in the presence of the IF domain. The inserted domain decreased, however, the stability of the host domain in the transition region and decelerated its refolding reaction by about 10-fold. This originates from the interruption of the chain connectivity by the IF domain and its inherent instability. To monitor folding processes in this domain selectively, a Trp residue was introduced as fluorescent probe. Kinetic double-mixing experiments revealed that, in intact SlyD, the IF domain folds and unfolds about 1000-fold more rapidly than the FKBP domain, and that it is strongly stabilized when linked with the folded FKBP domain. The unfolding limbs of the kinetic chevrons of SlyD show a strong downward curvature. This deviation from linearity is not caused by a transition-state movement, as often assumed, but by the accumulation of a silent unfolding intermediate at high denaturant concentrations. In this kinetic intermediate, the FKBP domain is still folded, whereas the IF domain is already unfolded.     

Tracing changes in protonation: a prerequisite to factorize thermodynamic data of inhibitor binding to aldose reductase

To prevent diabetic complications derived from enhanced glucose flux via the polyol pathway the development of aldose reductase inhibitors (ARIs) has been established as a promising therapeutic concept. Here, we study the binding process of inhibitors to aldose reductase (ALR2) with respect to changes of the protonation inventory upon complex formation. Knowledge of such processes is a prerequisite to factorize the binding free energy into enthalpic and entropic contributions on an absolute scale. Our isothermal titration calorimetry (ITC) measurements suggest a proton uptake upon complex formation with carboxylate-type inhibitors. As the protonation event will contribute strongly to the enthalpic signal recorded during ITC experiments, knowledge about the proton-accepting and releasing functional groups of the system is of utmost importance. However, this is intricate to retrieve, if, as in the present case, both, binding site and ligand possess several titratable groups. Here, we present pKa calculations complemented by mutagenesis and thermodynamic measurements suggesting a tyrosine residue located in the catalytic site (Tyr48) as a likely candidate to act as proton acceptor upon inhibitor binding, as it occurs deprotonated to a remarkable extent if only the cofactor NADP+ is bound. We furthermore provide evidence that the protonation state and binding thermodynamics depend strongly on the oxidation state of the cofactor;s nicotinamide moiety. Binding thermodynamics of IDD 388, IDD 393, tolrestat, sorbinil, and fidarestat are discussed in the context of substituent effects.     

Functional and structural roles of the Cys14-Cys38 disulfide of bovine pancreatic trypsin inhibitor

The disulfide bond between Cys14 and Cys38 of bovine pancreatic trypsin inhibitor lies on the surface of the inhibitor and forms part of the protease-binding region. The functional properties of three variants lacking this disulfide, with one or both of the Cys residues replaced with Ser, were examined, and X-ray crystal structures of the complexes with bovine trypsin were determined and refined to the 1.58-Å resolution limit. The crystal structure of the complex formed with the mutant with both Cys residues replaced was nearly identical with that of the complex containing the wild-type protein, with the Ser oxygen atoms positioned to replace the disulfide bond with a hydrogen bond. The two structures of the complexes with single replacements displayed small local perturbations with alternate conformations of the Ser side chains. Despite the absence of the disulfide bond, the crystallographic temperature factors show no evidence of increased flexibility in the complexes with the mutant inhibitors. All three of the variants were cleaved by trypsin more rapidly than the wild-type inhibitor, by as much as 10,000-fold, indicating that the covalent constraint normally imposed by the disulfide contributes to the remarkable resistance to hydrolysis displayed by the wild-type protein. The rates of hydrolysis display an unusual dependence on pH over the range of 3.5-8.0, decreasing at the more alkaline values, as compared with the increased hydrolysis rates for normal substrates under these conditions. These observations can be accounted for by a model for inhibition in which an acyl-enzyme intermediate forms at a significant rate but is rapidly converted back to the enzyme-inhibitor complex by nucleophilic attack by the newly created amino group. The model suggests that a lack of flexibility in the acyl-enzyme intermediate, rather than the enzyme-inhibitor complex, may be a key factor in the ability of bovine pancreatic trypsin inhibitor and similar inhibitors to resist hydrolysis.     

The impact of membrane lipid composition on antimicrobial function of an alpha-helical peptide

VP1, a putative alpha-helical antimicrobial peptide (alpha-AMP) inhibited growth of Bacillus subtilis and Escherichia coli at 500microM. The peptide induced stable surface pressure changes in monolayers formed from B. subtilis native lipid extract (circa 4.5mNm(-1)) but transient pressure changes in corresponding E. coli monolayers (circa 1.0mNm(-1)), which led to monolayer disintegration. Synthetic lipid monolayers mimetic of the extracts were used to generate compression isotherms. Thermodynamic analysis of B. subtilis isotherms indicated membrane stabilisation by VP1 (DeltaG(Mix)<0), via a mechanism dependent upon the phosphatidylglycerol to cardiolipin ratio. Corresponding analysis of E. coli isotherms indicated membrane destabilisation by the peptide (DeltaG(Mix)>0). Destabilisation correlated with PE levels present and appeared to involve a mechanism resembling those used by tilted peptides. These data emphasise that structure/function analysis of alpha-AMPs must consider not only their structural characteristics but also the lipid make-up of the target microbial membrane.