Influence of modulated structural dynamics on the kinetics of alpha-chymotrypsin catalysis. Insights through chemical glycosylation, molecular dynamics and domain motion analysis

Although the chemical nature of the catalytic mechanism of the serine protease alpha-chymotrypsin (alpha-CT) is largely understood, the influence of the enzyme's structural dynamics on its catalysis remains uncertain. Here we investigate whether alpha-CT's structural dynamics directly influence the kinetics of enzyme catalysis. Chemical glycosylation [Solá RJ & Griebenow K (2006) FEBS Lett 580, 1685-1690] was used to generate a series of glycosylated alpha-CT conjugates with reduced structural dynamics, as determined from amide hydrogen/deuterium exchange kinetics (k(HX)). Determination of their catalytic behavior (K(S), k(2), and k(3)) for the hydrolysis of N-succinyl-Ala-Ala-Pro-Phe p-nitroanilide (Suc-Ala-Ala-Pro-Phe-pNA) revealed decreased kinetics for the catalytic steps (k(2) and k(3)) without affecting substrate binding (K(S)) at increasing glycosylation levels. Statistical correlation analysis between the catalytic (DeltaG( not equal)k(i)) and structurally dynamic (DeltaG(HX)) parameters determined revealed that the enzyme acylation and deacylation steps are directly influenced by the changes in protein structural dynamics. Molecular modelling of the alpha-CT glycoconjugates coupled with molecular dynamics simulations and domain motion analysis employing the Gaussian network model revealed structural insights into the relation between the protein's surface glycosylation, the resulting structural dynamic changes, and the influence of these on the enzyme's collective dynamics and catalytic residues. The experimental and theoretical results presented here not only provide fundamental insights concerning the influence of glycosylation on the protein biophysical properties but also support the hypothesis that for alpha-CT the global structural dynamics directly influence the kinetics of enzyme catalysis via mechanochemical coupling between domain motions and active site chemical groups.     

Engineering of protein thermodynamic, kinetic, and colloidal stability: Chemical Glycosylation with monofunctionally activated glycans

In this work we establish the relationship between chemical glycosylation and protein thermodynamic, kinetic, and colloidal stability. While there have been reports in the literature that chemical glycosylation modulates protein stability, mechanistic details still remain uncertain. To address this issue, we designed and coupled monofunctional activated glycans (lactose and dextran) to the model protein alpha-chymotrypsin (alpha-CT). This resulted in a series of glycoconjugates with variations in the glycan size and degree of glycosylation. Thermodynamic unfolding, thermal inactivation, and temperature-induced aggregation experiments revealed that chemical glycosylation increased protein thermodynamic (Delta G(25 degrees C)), kinetic (t(1/2)(45 degrees C)), and colloidal stability. These results highlight the potential of chemical glycosylation with monofunctional activated glycans as a technology for increasing the long-term stability of liquid protein formulations for industrial and biotherapeutic applications.     

The effect of glycosylation on the folding kinetics of erythropoietin

Glycosylation is a common posttranslational modification that generally increases protein solubility and thermodynamic stability. Less is known about how this modification influences protein folding, particularly folding processes involving intermediate species. In the present report, folding comparisons of a nonglycosylated erythropoietin (EPO) mutant are made with the fully glycosylated EPO, which was recently shown to fold by a three-state on-pathway mechanism. The absence of glycosylation did not alter the folding mechanism of EPO but did greatly decrease the stability of the intermediate species, change the rate-limiting step of the folding reaction, and accelerate the folding kinetics to both the intermediate state and the native state. Surprisingly, glycosylation stabilized the intermediate species to a greater extent than it increased the EPO equilibrium stability. These results suggest that glycosylation impedes the latter EPO folding steps rather than accelerating them by biasing particular folding pathways, as previously proposed for folding reactions initiated from unfolded ensembles with minimal residual structure. Due to the specific biological processes modulated by EPO glycosylation, however, there may be little evolutionary pressure to fold on a faster, more direct pathway at the expense of biological function, particularly given the protective role glycosylation has at preventing EPO aggregation. Lastly, evidence that is consistent with glycosylation destabilizing the unfolded state to some degree and contributing to the greater equilibrium stability of the glycosylated EPO is presented.     

Stabilization of alpha-chymotrypsin upon PEGylation correlates with reduced structural dynamics

Protein stability remains one of the main factors limiting the realization of the full potential of protein therapeutics. Poly(ethylene glycol) (PEG) conjugation to proteins has evolved into an important tool to overcome instability issues associated with proteins. The observed increase in thermodynamic stability of several proteins upon PEGylation has been hypothesized to arise from reduced protein structural dynamics, although experimental evidence for this hypothesis is currently missing. To test this hypothesis, the model protein alpha-chymotrypsin (alpha-CT) was covalently modified with PEGs with molecular weights (M(W)) of 700, 2,000 and 5,000 and the degree of modification was systematically varied. The procedure did not cause significant tertiary structure changes. Thermodynamic unfolding experiments revealed that PEGylation increased the thermal transition temperature (T(m)) of alpha-CT by up to 6 degrees C and the free energy of unfolding [DeltaG(U) (25 degrees C)] by up to 5 kcal/mol. The increase in stability was found to be independent of the PEG M(W) and it leveled off after an average of four PEG molecules were bound to alpha-CT. Fourier-transformed infrared (FTIR) H/D exchange experiments were conducted to characterize the conformational dynamics of the PEG-conjugates. It was found that the magnitude of thermodynamic stabilization correlates with a reduction in protein structural dynamics and was independent of the PEG M(W). Thus, the initial hypothesis proved positive. Similar to the thermodynamic stabilization of proteins by covalent modification with glycans, PEG thermodynamically stabilizes alpha-CT by reducing protein structural dynamics. These results provide guidance for the future development of stable protein formulations.     

Relation between stability, dynamics and enzyme activity in 3-phosphoglycerate kinases from yeast and Thermus thermophilus.

3-Phosphoglycerate kinases from yeast and the extreme thermophilic bacterium Thermus thermophilus HB8 have been used as models for investigating the relationship between stability, dynamics and activity. It was found that while at a given temperature the thermophilic protein is more stable, its conformational dynamics as measured by the ability of acrylamide to quench the fluorescence of a buried tryptophan as well as its specific activity, are both lower than for the mesophilic protein. As the temperature is increased, the thermodynamic stability of the thermophilic protein approaches that of the mesophilic protein at its working temperature. Its conformational dynamics and specific activity however were both shown to increase, until at the physiologically operational temperature, they become similar to those of the mesophilic enzyme at its operational temperature. These results confirm the proposal that a direct relationship and balance holds between thermodynamic stability, dynamics and specific activity in globular proteins. They demonstrate also the constraining effect of increased stability upon conformational dynamics and enzyme activity.     

Comparative NMR study on the impact of point mutations on protein stability of Pseudomonas mendocina lipase

In this work we compare the dynamics and conformational stability of Pseudomonas mendocina lipase enzyme and its F180P/S205G mutant that shows higher activity and stability for use in washing powders. Our NMR analyses indicate virtually identical structures but reveal remarkable differences in local dynamics, with striking correspondence between experimental data (i.e., (15)N relaxation and H/D exchange rates) and data from Molecular Dynamics simulations. While overall the cores of both proteins are very rigid on the pico- to nanosecond timescale and are largely protected from H/D exchange, the two point mutations stabilize helices alpha1, alpha4, and alpha5 and locally destabilize the H-bond network of the beta-sheet (beta7-beta9). In particular, it emerges that helix alpha5, undergoing some fast destabilizing motions (on the pico- to nanosecond timescale) in wild-type lipase, is substantially rigidified by the mutation of Phe180 for a proline at its N terminus. This observation could be explained by the release of some penalizing strain, as proline does not require any "N-capping" hydrogen bond acceptor in the i+3 position. The combined experimental and simulated data thus indicate that reduced molecular flexibility of the F180P/S205G mutant lipase underlies its increased stability, and thus reveals a correlation between microscopic dynamics and macroscopic thermodynamic properties. This could contribute to the observed altered enzyme activity, as may be inferred from recent studies linking enzyme kinetics to their local molecular dynamics.     

Stepwise adaptations to low temperature as revealed by multiple mutants of psychrophilic α-amylase from Antarctic Bacterium

The mutants Mut5 and Mut5CC from a psychrophilic α-amylase bear representative stabilizing interactions found in the heat-stable porcine pancreatic α-amylase but lacking in the cold-active enzyme from an Antarctic bacterium. From an evolutionary perspective, these mutants can be regarded as structural intermediates between the psychrophilic and the mesophilic enzymes. We found that these engineered interactions improve all the investigated parameters related to protein stability as follows: compactness; kinetically driven stability; thermodynamic stability; resistance toward chemical denaturation, and the kinetics of unfolding/refolding. Concomitantly to this improved stability, both mutants have lost the kinetic optimization to low temperature activity displayed by the parent psychrophilic enzyme. These results provide strong experimental support to the hypothesis assuming that the disappearance of stabilizing interactions in psychrophilic enzymes increases the amplitude of concerted motions required by catalysis and the dynamics of active site residues at low temperature, leading to a higher activity.     

A Genome-Wide Sequence–Structure Analysis Suggests Aggregation Gatekeepers Constitute an Evolutionary Constrained Functional Class

Protein aggregation is geared by aggregation-prone regions that self-associate by β-strand interactions. Charged residues and prolines are enriched at the flanks of aggregation-prone regions resulting in decreased aggregation. It is still unclear what drives the overrepresentation of these “aggregation gatekeepers”, that is, whether their presence results from structural constraints determining protein stability or whether they constitute a bona fide functional class selectively maintained to control protein aggregation. As functional residues are typically conserved regardless of their cost to protein stability, we compared sequence conservation and thermodynamic cost of these residues in 2659 protein families in Escherichia coli. Across protein families, we find gatekeepers to be under strong selective conservation while at the same time representing a significant thermodynamic cost to protein structure. This finding supports the notion that aggregation gatekeepers are not structurally determined but evolutionary selected to control protein aggregation.     

Effect of PEG modification on subtilisin Carlsberg activity, enantioselectivity, and structural dynamics in 1,4-dioxane

The employment of enzymes as catalysts within organic media has traditionally been hampered by the reduced enzymatic activities when compared to catalysis in aqueous solution. Although several complementary hypotheses have provided mechanistic insights into the causes of diminished activity, further development of biocatalysts would greatly benefit from effective chemical strategies (e.g., PEGylation) to ameliorate this event. Herein we explore the effects of altering the solvent composition from aqueous buffer to 1,4-dioxane on structural, dynamical, and catalytic properties of the model enzyme subtilisin Carlsberg (SBc). Furthermore, we also investigate the effects of dissolving the enzyme in 1,4-dioxane through chemical modification with poly(ethylene)-glycol (PEG, M(W) = 20 kDa) on these enzyme properties. In 1,4-dioxane a 10(4)-fold decrease in the enzyme's catalytic activity was observed for the hydrolysis reaction of vinyl butyrate with D(2)O and a 50% decrease in enzyme structural dynamics as evidenced by reduced amide H/D exchange kinetics occurred. Attaching increasing amounts of PEG to the enzyme reversed some of the activity loss. Evaluation of the structural dynamic behavior of the PEGylated enzyme within the organic solvent revealed an increase in structural dynamics at increased PEGylation. Correlation analysis between the catalytic and structural dynamic parameters revealed that the enzyme's catalytic activity and enantioselectivity depended on the changes in protein structural dynamics within 1,4-dioxane. These results demonstrate the importance of protein structural dynamics towards regulating the catalytic behavior of enzymes within organic media.     

Adenylation-dependent conformation and unfolding pathways of the NAD+-dependent DNA ligase from the thermophile Thermus scotoductus

In the last few years, an increased attention has been focused on NAD(+)-dependent DNA ligases. This is mostly due to their potential use as antibiotic targets, because effective inhibition of these essential enzymes would result in the death of the bacterium. However, development of an efficient drug requires that the conformational modifications involved in the catalysis of NAD(+)-dependent DNA ligases are understood. From this perspective, we have investigated the conformational changes occurring in the thermophilic Thermus scotoductus NAD(+)-DNA ligase upon adenylation, as well as the effect of cofactor binding on protein resistance to thermal and chemical (guanidine hydrochloride) denaturation. Our results indicate that cofactor binding induces conformational rearrangement within the active site and promotes a compaction of the enzyme. These data support an induced "open-closure" process upon adenylation, leading to the formation of the catalytically active enzyme that is able to bind DNA. These conformational changes are likely to be associated with the protein function, preventing the formation of nonproductive complexes between deadenylated ligases and DNA. In addition, enzyme adenylation significantly increases resistance of the protein to thermal denaturation and GdmCl-induced unfolding, establishing a thermodynamic link between ligand binding and increased conformational stability. Finally, chemical unfolding of deadenylated and adenylated enzyme is accompanied by accumulation of at least two equilibrium intermediates, the molten globule and premolten globule states. Maximal populations of these intermediates are shifted toward higher GdmCl concentrations in the case of the adenylated ligase. These data provide further insights into the properties of partially folded intermediates.     

Evolutionarily conserved linkage between enzyme fold, flexibility, and catalysis

Proteins are intrinsically flexible molecules. The role of internal motions in a protein''s designated function is widely debated. The role of protein structure in enzyme catalysis is well established, and conservation of structural features provides vital clues to their role in function. Recently, it has been proposed that the protein function may involve multiple conformations: the observed deviations are not random thermodynamic fluctuations; rather, flexibility may be closely linked to protein function, including enzyme catalysis. We hypothesize that the argument of conservation of important structural features can also be extended to identification of protein flexibility in interconnection with enzyme function. Three classes of enzymes (prolyl-peptidyl isomerase, oxidoreductase, and nuclease) that catalyze diverse chemical reactions have been examined using detailed computational modeling. For each class, the identification and characterization of the internal protein motions coupled to the chemical step in enzyme mechanisms in multiple species show identical enzyme conformational fluctuations. In addition to the active-site residues, motions of protein surface loop regions (>10 Å away) are observed to be identical across species, and networks of conserved interactions/residues connect these highly flexible surface regions to the active-site residues that make direct contact with substrates. More interestingly, examination of reaction-coupled motions in non-homologous enzyme systems (with no structural or sequence similarity) that catalyze the same biochemical reaction shows motions that induce remarkably similar changes in the enzyme–substrate interactions during catalysis. The results indicate that the reaction-coupled flexibility is a conserved aspect of the enzyme molecular architecture. Protein motions in distal areas of homologous and non-homologous enzyme systems mediate similar changes in the active-site enzyme–substrate interactions, thereby impacting the mechanism of catalyzed chemistry. These results have implications for understanding the mechanism of allostery, and for protein engineering and drug design.

Author''s Summary

Enzymes are nature''s molecular machines that catalyze biochemical reactions with remarkable efficiency. Recent evidence suggests that enzyme function may involve not only direct structural interactions between the enzyme and its substrate, but also internal motions of the enzyme itself. Here, we describe a computational investigation of three classes of enzymes that catalyze completely different biochemical reactions. Remarkably, the mobile enzyme regions and the nature of these motions are the same across species ranging from single-celled organisms to complex life-forms. Also surprisingly, non-homologous enzymes that catalyze the same chemical reaction but do not share sequence or structural similarity reveal a similar impact of enzyme motions on their reaction mechanisms. Flexible enzyme regions are found to be connected by conserved networks of coupled interactions that connect surface regions to active-site residues. These networks may provide a mechanism for the solvent on an enzyme''s surface to couple to the reaction catalyzed by the enzyme. These results have implications for understanding the mechanism of allostery (long-range effects), and for protein engineering and drug design.     

ATP binding enables broad antibiotic selectivity of aminoglycoside phosphotransferase(3')-IIIa: an elastic network analysis

The bacterial enzyme aminoglycoside phosphotransferase(3′)-IIIa (APH) confers resistance against a wide range of aminoglycoside antibiotics. In this study, we use the Gaussian network model to investigate how the binding of nucleotides and antibiotics influences the dynamics and thereby the ligand binding properties of APH. Interestingly, in NMR experiments, the dynamics differ significantly in various APH complexes, although crystallographic studies indicate that no larger conformational changes occur upon ligand binding. Isothermal titration calorimetry also shows different thermodynamic contributions to ligand binding. Formation of aminoglycoside-APH complexes is enthalpically driven, while the enthalpic change upon aminoglycoside binding to the nucleotide-APH complex is much smaller. The differential effects of nucleotide binding and antibiotic binding to APH can be explained theoretically by single-residue fluctuations and correlated motions of the enzyme. The surprising destabilization of β-sheet residues upon nucleotide binding, as seen in hydrogen/deuterium exchange experiments, shows that the number of closest neighbors does not fully explain residue flexibility. Additionally, we must consider correlated motions of dynamic protein domains, which show that not only connectivity but also the overall protein architecture is important for protein dynamics.     

A directed approach for engineering conditional protein stability using biologically silent small molecules

The ability to regulate the function of specific proteins using cell-permeable molecules can be a powerful method for interrogating biological systems. To bring this type of "chemical genetic" control to a wide range of proteins, we recently developed an experimental system in which the stability of a small protein domain expressed in mammalian cells depends on the presence of a high affinity ligand. This ligand-dependent stability is conferred to any fused partner protein. The FK506- and rapamycin-binding protein (FKBP12) has been the subject of extensive biophysical analyses, including both kinetic and thermodynamic studies of the wild-type protein as well as dozens of mutants. The goal of this study was to determine if the thermodynamic stabilities (DeltaDeltaG(U-F)) of various amino acid substitutions within a given protein are predictive for engineering additional ligand-dependent destabilizing domains. We used FKBP12 as a model system and found that in vitro thermodynamic stability correlates weakly with intracellular degradation rates of the mutants and that the ability of a given mutation to destabilize the protein is context-dependent. We evaluated several new FKBP12 ligands for their ability to stabilize these mutants and found that a cell-permeable molecule called Shield-1 is the most effective stabilizing ligand. We then performed an unbiased microarray analysis of NIH3T3 cells treated with various concentrations of Shield-1. These studies show that Shield-1 does not elicit appreciable cellular responses.     

Contribution of the multi-turn segment in the reversible thermal stability of hyperthermophile rubredoxin: NMR thermal chemical exchange analysis of sequence hybrids

Pyrococcus furiosus (Pf) rubredoxin is the most thermostable protein characterized to date. Reflecting the complications arising from irreversible denaturation of this protein, predictions of which structural regions confer differential thermal stability have utilized kinetic stability measurements, hydrogen exchange protection factors, long range hydrogen bond NMR spin couplings, and molecular dynamics simulations, and have primarily implicated the three-stranded beta-sheet and the adjacent metal binding site. Herein, NMR chemical exchange experiments demonstrate reversible two-state unfolding at the thermal transition temperature (T(m)) for hybrids of Pf and the mesophile Clostridium pasteurianum (Cp) rubredoxins which interchange residues 14-33, the so-called multi-turn segment. This complementary pair of hybrid rubredoxins exhibits largely additive incremental thermal stabilizations vs. the parental proteins. Both stabilization free energy measurements as well as incremental T(m) values indicate that a minimum of 37% of the total differential thermal stability resides in this multi-turn segment. Such a proportionality between DeltaDeltaG and incremental T(m) values is predicted for hybrid pairs exhibiting thermodynamic additivity in which the differential stability is predominantly enthalpic.     

Monte Carlo simulations of the peptide recognition at the consensus binding site of the constant fragment of human immunoglobulin G: the energy landscape analysis of a hot spot at the intermolecular interface

Monte Carlo simulations of molecular recognition at the consensus binding site of the constant fragment (Fc) of human immunoglobulin G (Ig) protein have been performed to analyze structural and thermodynamic aspects of binding for the 13-residue cyclic peptide DCAWHLGELVWCT. The energy landscape analysis of a hot spot at the intermolecular interface using alanine scanning and equilibrium-simulated tempering dynamics with the simplified, knowledge-based energy function has enabled the role of the protein hot spot residues in providing the thermodynamic stability of the native structure to be determined. We have found that hydrophobic interactions between the peptide and the Met-252, Ile-253, His-433, and His-435 protein residues are critical to guarantee the thermodynamic stability of the crystallographic binding mode of the complex. Binding free energy calculations, using a molecular mechanics force field and a solvation energy model, combined with alanine scanning have been conducted to determine the energetic contribution of the protein hot spot residues in binding affinity. The conserved Asn-434, Ser-254, and Tyr-436 protein residues contribute significantly to the binding affinity of the peptide-protein complex, serving as an energetic hot spot at the intermolecular interface. The results suggest that evolutionary conserved hot spot protein residues at the intermolecular interface may be partitioned in fulfilling thermodynamic stability of the native binding mode and contributing to the binding affinity of the complex.     

Architectural and thermodynamic principles underlying intramembrane protease function

Intramembrane proteases hydrolyze peptide bonds within the membrane as a signaling paradigm universal to all life forms and with implications in disease. Deciphering the architectural strategies supporting intramembrane proteolysis is an essential but unattained goal. We integrated new, quantitative and high-throughput thermal light-scattering technology, reversible equilibrium unfolding and refolding and quantitative protease assays to interrogate rhomboid architecture with 151 purified variants. Rhomboid proteases maintain low intrinsic thermodynamic stability (ΔG = 2.1-4.5 kcal mol(-1)) resulting from a multitude of generally weak transmembrane packing interactions, making them highly responsive to their environment. Stability is consolidated by two buried glycines and several packing leucines, with a few multifaceted hydrogen bonds strategically deployed to two peripheral regions. Opposite these regions lie transmembrane segment 5 and connected loops that are notably exempt of structural responsibility, suggesting intramembrane proteolysis involves considerable but localized protein dynamics. Our analyses provide a comprehensive 'heat map' of the physiochemical anatomy underlying membrane-immersed enzyme function at, what is to our knowledge, unprecedented resolution.     

Thermal stability of proteins does not correlate with the energy of intramolecular interactions

Small monomeric proteins from mesophilic and thermophilic organisms were studied. They have close structural and physical and chemical properties but vary in thermal stability. A thermodynamic analysis of heat unfolding was made and integral enthalpy of unfolding (DeltaH(unf)), heat capacity of hydration (DeltaC(p)(hyd)) and enthalpy of hydration (DeltaH(hyd)) and of the buried surface area (DeltaASA) of nonpolar and polar groups as well as the enthalpy of disruption of intramolecular interaction (DeltaH(int) in gas phase) at 298 K were determined. The absence of correlation between protein thermostability and energetic components suggests that regulatory mechanism of protein thermal stabilization has entropic nature.     

Alanine Scan of Core Positions in Ubiquitin Reveals Links between Dynamics,Stability, and Function

Mutations at solvent-inaccessible core positions in proteins can impact function through many biophysical mechanisms including alterations to thermodynamic stability and protein dynamics. As these properties of proteins are difficult to investigate, the impacts of core mutations on protein function are poorly understood for most systems. Here, we determined the effects of alanine mutations at all 15 core positions in ubiquitin on function in yeast. The majority (13 of 15) of alanine substitutions supported yeast growth as the sole ubiquitin. Both the two null mutants (I30A and L43A) were less stable to temperature-induced unfolding in vitro than wild type (WT) but were well folded at physiological temperatures. Heteronuclear NMR studies indicated that the L43A mutation reduces temperature stability while retaining a ground-state structure similar to WT. This structure enables L43A to bind to common ubiquitin receptors in vitro. Many of the core alanine ubiquitin mutants, including one of the null variants (I30A), exhibited an increased accumulation of high-molecular-weight species, suggesting that these mutants caused a defect in the processing of ubiquitin-substrate conjugates. In contrast, L43A exhibited a unique accumulation pattern with reduced levels of high-molecular-weight species and undetectable levels of free ubiquitin. When conjugation to other proteins was blocked, L43A ubiquitin accumulated as free ubiquitin in yeast. Based on these findings, we speculate that ubiquitin's stability to unfolding may be required for efficient recycling during proteasome-mediated substrate degradation.     

Analysis of side chain mobility among protein G B1 domain mutants with widely varying stabilities

"Host-guest" studies of the B1 domain from Streptococcal protein G have been used previously to establish a thermodynamic scale for the beta-sheet-forming propensities of the 20 common amino acids. To investigate the contribution of side chain conformational entropy to the relative stabilities of B1 domain mutants, we have determined the dynamics of side chain methyl groups in 10 of the 20 mutants used in a previous study. Deuterium relaxation rates were measured using two-dimensional NMR techniques for 13CH2D groups. Analysis of the relaxation data using the Lipari-Szabo model-free formalism showed that mutations introduced at the guest position caused small but statistically significant changes in the methyl group dynamics. In addition, there was a low level of covariation of the Lipari-Szabo order parameters among the 10 mutants. The variations in conformational free energy estimated from the order parameters were comparable in magnitude to the variations in global stability of the 10 mutants but did not correlate with the global stability of the domain or with the structural properties of the guest amino acids. The data support the view that conformational entropy in the folded state is one of many factors that can influence the folding thermodynamics of proteins.     

A miniaturized technique for assessing protein thermodynamics and function using fast determination of quantitative cysteine reactivity

Protein thermodynamic stability is a fundamental physical characteristic that determines biological function. Furthermore, alteration of thermodynamic stability by macromolecular interactions or biochemical modifications is a powerful tool for assessing the relationship between protein structure, stability, and biological function. High-throughput approaches for quantifying protein stability are beginning to emerge that enable thermodynamic measurements on small amounts of material, in short periods of time, and using readily accessible instrumentation. Here we present such a method, fast quantitative cysteine reactivity, which exploits the linkage between protein stability, sidechain protection by protein structure, and structural dynamics to characterize the thermodynamic and kinetic properties of proteins. In this approach, the reaction of a protected cysteine and thiol-reactive fluorogenic indicator is monitored over a gradient of temperatures after a short incubation time. These labeling data can be used to determine the midpoint of thermal unfolding, measure the temperature dependence of protein stability, quantify ligand-binding affinity, and, under certain conditions, estimate folding rate constants. Here, we demonstrate the fQCR method by characterizing these thermodynamic and kinetic properties for variants of Staphylococcal nuclease and E. coli ribose-binding protein engineered to contain single, protected cysteines. These straightforward, information-rich experiments are likely to find applications in protein engineering and functional genomics.