Thermodynamic effects of mutations on the denaturation of T4 lysozyme

Differential scanning calorimetry of proteins and Zimm-Bragg model in water

Differential Scanning Calorimetry (DSC) is a regular and powerful tool to measure the specific heat profile of various materials. Hydrogen bonds play a crucial role in stabilizing the three-dimensional structure of proteins. Naturally, information about the strength of hydrogen bonds is contained in the measured DSC profiles. Despite its obvious importance, there is no approach that would allow the extraction of such information from the heat capacity measurements. In order to connect the measured profile to microscopic properties of a polypeptide chain, a proper model is required to fit. Using recent advances in the Zimm-Bragg (ZB) theory of protein folding in water, we propose a new and efficient algorithm to process the DSC experimental data and to extract the H-bonding energy among other relevant constants. Thus, for the randomly picked set of 33 proteins, we have found a quite narrow distribution of hydrogen bonding energies from 1 to 8 kJ/mol with the average energy of intra-protein hydrogen bonds h¯=4.2±1.5 kJ/mol and the average energy of water-protein bonds as hps¯=3.8±1.5 kJ/mol. This is an important illustration of a tiny disbalance between the water-protein and intraprotein hydrogen bonds. Fitted values of the nucleation parameter σ belong to the range from 0.001 to 0.01, as expected. The reported method can be considered as complementary to the classical two-state approach and together with other parameters provides the protein-water and intraprotein H-bonding energies, not accessible within the two-state paradigm.

Dark nanodiscs for evaluating membrane protein thermostability by differential scanning fluorimetry

Measuring protein thermostability provides valuable information on the biophysical rules that govern the structure-energy relationships of proteins. However, such measurements remain a challenge for membrane proteins. Here, we introduce a new experimental system to evaluate membrane protein thermostability. This system leverages a recently developed nonfluorescent membrane scaffold protein to reconstitute proteins into nanodiscs and is coupled with a nano-format of differential scanning fluorimetry (nanoDSF). This approach offers a label-free and direct measurement of the intrinsic tryptophan fluorescence of the membrane protein as it unfolds in solution without signal interference from the "dark" nanodisc. In this work, we demonstrate the application of this method using the disulfide bond formation protein B (DsbB) as a test membrane protein. NanoDSF measurements of DsbB reconstituted in dark nanodiscs loaded with 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dimyristoyl-sn-glycero-3-phosphorylglycerol (DMPG) lipids show a complex biphasic thermal unfolding pattern with a minor unfolding transition followed by a major transition. The inflection points of the thermal denaturation curve reveal two distinct unfolding midpoint melting temperatures (Tm) of 70.5°C and 77.5°C, consistent with a three-state unfolding model. Further, we show that the catalytically conserved disulfide bond between residues C41 and C130 drives the intermediate state of the unfolding pathway for DsbB in a DMPC and DMPG nanodisc. To extend the utility of this method, we evaluate and compare the thermostability of DsbB in different lipid environments. We introduce this method as a new tool that can be used to understand how compositionally and biophysically complex lipid environments drive membrane protein stability.

Calculation of Protein Folding Thermodynamics Using Molecular Dynamics Simulations

Despite advances in artificial intelligence methods, protein folding remains in many ways an enigma to be solved. Accurate computation of protein folding energetics could help drive fields such as protein and drug design and genetic interpretation. However, the challenge of calculating the state functions governing protein folding from first-principles remains unaddressed. We present here a simple approach that allows us to accurately calculate the energetics of protein folding. It is based on computing the energy of the folded and unfolded states at different temperatures using molecular dynamics simulations. From this, two essential quantities (ΔH and ΔCp) are obtained and used to calculate the conformational stability of the protein (ΔG). With this approach, we have successfully calculated the energetics of two- and three-state proteins, representatives of the major structural classes, as well as small stability differences (ΔΔG) due to changes in solution conditions or variations in an amino acid residue.

Dark nanodiscs as a model membrane for evaluating membrane protein thermostability by differential scanning fluorimetry

Measuring protein thermostability provides valuable information on the biophysical rules that govern structure-energy relationships of proteins. However, such measurements remain a challenge for membrane proteins. Here, we introduce a new experimental system to evaluate membrane protein thermostability. This system leverages a recently-developed non-fluorescent membrane scaffold protein (MSP) to reconstitute proteins into nanodiscs and is coupled with a nano-format of differential scanning fluorimetry (nanoDSF). This approach offers a label-free and direct measurement of the intrinsic tryptophan fluorescence of the membrane protein as it unfolds in solution without signal interference from the "dark" nanodisc. In this work, we demonstrate the application of this method using the disulfide bond formation protein B (DsbB) as a test membrane protein. NanoDSF measurements of DsbB reconstituted in dark nanodiscs show a complex biphasic thermal unfolding pattern in the presence of lipids with a minor unfolding transition followed by a major transition. The inflection points of the thermal denaturation curve reveal two distinct unfolding midpoint melting temperatures (Tm) of 70.5 °C and 77.5 °C, consistent with a three-state unfolding model. Further, we show that the catalytically conserved disulfide bond between residues C41 and C130 drives the intermediate state of the unfolding pathway for DsbB in a nanodisc. We introduce this method as a new tool that can be used to understand how compositionally, and biophysically complex lipid environments drive membrane protein stability.

Effects of Ionic Liquids on Metalloproteins

In the past decade, innovative protein therapies and bio-similar industries have grown rapidly. Additionally, ionic liquids (ILs) have been an area of great interest and rapid development in industrial processes over a similar timeline. Therefore, there is a pressing need to understand the structure and function of proteins in novel environments with ILs. Understanding the short-term and long-term stability of protein molecules in IL formulations will be key to using ILs for protein technologies. Similarly, ILs have been investigated as part of therapeutic delivery systems and implicated in numerous studies in which ILs impact the activity and/or stability of protein molecules. Notably, many of the proteins used in industrial applications are involved in redox chemistry, and thus often contain metal ions or metal-associated cofactors. In this review article, we focus on the current understanding of protein structure-function relationship in the presence of ILs, specifically focusing on the effect of ILs on metal containing proteins.

A look back at the molten globule state of proteins: thermodynamic aspects

Interest in protein folding intermediates lies in their significance to protein folding pathways. The molten globule (MG) state is one such intermediate lying on the kinetic (and sometimes thermodynamic) pathway between native and unfolded states. Development of our qualitative and quantitative understanding of the MG state can provide deeper insight into the folding pathways and hence potentially facilitate solution of the protein folding problem. An extensive look at literature suggests that most studies into protein MG states have been largely qualitative. Attempts to obtain quantitative insights into MG states have involved application of high-sensitivity calorimetry (differential scanning calorimetry and isothermal titration calorimetry). This review addresses the progress made in this direction by discussing the knowledge gained to date, along with the future promise of calorimetry, in providing quantitative information on the structural features of MG states. Particular attention is paid to the question of whether such states share common structural features or not. The difference in the nature of the transition from the MG state to the unfolded state, in terms of cooperativity, has also been addressed and discussed.

Unfolding Pathways of Hen Egg-White Lysozyme in Ethanol

The aggregation of amyloid fibrils can lead to various diseases including Alzheimer's, Parkinson's disease, and transmissible spongiform encephalopathy. Amyloid fibrils can develop from a variety of proteins in the body as they misfold into a primarily β-sheet structure and aggregate. Human lysozyme has been shown to have far reaching effects in the human health-a homologous enzyme, hen egg-white lysozyme, has been shown to denature to a primarily β-sheet structure at low pH and high alcohol content solution. We have studied these systems in atomic-level detail with a combination of constant pH and microsecond long molecular dynamics simulation in explicit solvent, which cumulatively total over 10 μs of simulation time. These studies have allowed us to determine two potential unfolding pathways depending on the protonation state of a key glutamic acid residue as well as the effect of solution dynamics and pH on the unfolding process.

Thermal protein unfolding by differential scanning calorimetry and circular dichroism spectroscopy Two-state model versus sequential unfolding

Thermally-induced protein unfolding is commonly described with the two-state model. This model assumes only two types of protein molecules in solution, the native (N) and the denatured, unfolded (U) protein. In reality, protein unfolding is a multistep process, even if intermediate states are only sparsely populated. As an alternative approach we explore the Zimm-Bragg theory, originally developed for the α-helix-to-random coil transition of synthetic polypeptides. The theory includes intermediate structures with concentrations determined by the cooperativity of the unfolding reaction. We illustrate the differences between the two-state model and the Zimm-Bragg theory with measurements of apolipoprotein A-1 and lysozyme by differential scanning calorimetry (DSC) and CD spectroscopy. Nine further protein examples are taken from the literature. The Zimm-Bragg theory provides a perfect fit of the calorimetric unfolding transitions for all proteins investigated. In contrast, the transition curves and enthalpies predicted by the two-state model differ considerably from the experimental results. Apolipoprotein A-1 is ~50% α-helical at ambient temperature and its unfolding follows the classical α-helix-to-random coil equilibrium. The unfolding of proteins with little α-helix content, such as lysozyme, can also be analyzed with the Zimm-Bragg theory by introducing the concept of 'folded' and 'unfolded' peptide units assuming an average unfolding enthalpy per peptide unit. DSC is the method of choice to measure the unfolding enthalpy, , but CD spectroscopy in combination with the two-state model is often used to deduce the unfolding enthalpy. This can lead to erroneous result. Not only are different enthalpies required to describe the CD and DSC transition curves but these values deviate distinctly from the experimental result. In contrast, the Zimm-Bragg theory predicts the DSC and CD unfolding transitions with the same set of parameters.

The Compact and Biologically Relevant Structure of Inter-α-inhibitor Is Maintained by the Chondroitin Sulfate Chain and Divalent Cations

Inter-α-inhibitor is a proteoglycan of unique structure. The protein consists of three subunits, heavy chain 1, heavy chain 2, and bikunin covalently joined by a chondroitin sulfate chain originating at Ser-10 of bikunin. Inter-α-inhibitor interacts with an inflammation-associated protein, tumor necrosis factor-inducible gene 6 protein, in the extracellular matrix. This interaction leads to transfer of the heavy chains from the chondroitin sulfate of inter-α-inhibitor to hyaluronan and consequently to matrix stabilization. Divalent cations and heavy chain 2 are essential co-factors in this transfer reaction. In the present study, we have investigated how divalent cations in concert with the chondroitin sulfate chain influence the structure and stability of inter-α-inhibitor. The results showed that Mg(2+) or Mn(2+), but not Ca(2+), induced a conformational change in inter-α-inhibitor as evidenced by a decrease in the Stokes radius and a bikunin chondroitin sulfate-dependent increase of the thermodynamic stability. This structure was shown to be essential for the ability of inter-α-inhibitor to participate in extracellular matrix stabilization. In addition, the data revealed that bikunin was positioned adjacent to both heavy chains and that the two heavy chains also were in close proximity. The chondroitin sulfate chain interacted with all protein components and inter-α-inhibitor dissociated when it was degraded. Conventional purification protocols result in the removal of the Mg(2+) found in plasma and because divalent cations influence the conformation and affect function it is important to consider this when characterizing the biological activity of inter-α-inhibitor.

Existence of metastable intermediate lysozyme conformation highlights the role of alcohols in altering protein stability

Alcohols have a manifold effect on the conformational and thermodynamic stability of native proteins. Here, we study the effect of moderate concentrations of trifluoroethanol (TFE) on the thermal stability of hen egg-white lysozyme (HEWL), by far-UV circular dichroism and by steady-state and time-resolved photoluminescence of intrinsic tryptophans. Our results highlight that TFE affects lysozyme stability by preferential solvation of the protein molecule. Furthermore, we discovered the existence at 20% TFE of an equilibrium partially folded state of lysozyme, intermediate between the native and the unfolded state. A three-state model is therefore used to interpolate the thermal denaturation data. Our analysis explains how the stabilization of the intermediate conformation enhances the entropic contribution to unfolding, and thus decreases the unfolding temperature, while, at the same time, TFE enhances the conformational stability of the native fold at room temperature. Eventually, we challenged the ability of these intermediate structures to form supramolecular aggregates by heating experiments at different TFE concentrations.

Microcalorimetry of proteins and their complexes

Ultrasensitive microcalorimetric techniques for measuring the heat capacities of proteins in dilute solutions over a broad temperature range (DSC) and the heats of protein reactions at fixed temperatures (ITC) are described and the methods of working with these instruments are considered. Particular attention is paid to analyzing the thermal properties of individual proteins, their stability, the energetics of their folding, and their association with specific macromolecular partners. Use of these calorimetric methods is illustrated with examples of small compact globular proteins, small proteins having loose noncompact structure, multidomain proteins, and protein complexes, particularly with DNA.

Alanine-scanning mutagenesis of the beta-sheet region of phage T4 lysozyme suggests that tertiary context has a dominant effect on beta-sheet formation

In general, alpha-helical conformations in proteins depend in large part on the amino acid residues within the helix and their proximal interactions. For example, an alanine residue has a high propensity to adopt an alpha-helical conformation, whereas that of a glycine residue is low. The sequence preferences for beta-sheet formation are less obvious. To identify the factors that influence beta-sheet conformation, a series of scanning polyalanine mutations were made within the strands and associated turns of the beta-sheet region in T4 lysozyme. For each construct the stability of the folded protein was reduced substantially, consistent with removal of native packing interactions. However, the crystal structures showed that each of the mutants retained the beta-sheet conformation. These results suggest that the structure of the beta-sheet region of T4 lysozyme is maintained to a substantial extent by tertiary interactions with the surrounding parts of the protein. Such tertiary interactions may be important in determining the structures of beta-sheets in general.

Folding of beta-sheets in membranes: specificity and promiscuity in peptide model systems

The interactions that drive the folding of beta-barrel membrane proteins have not been well studied because there have been few available model systems for membrane beta-sheets. In this work, we expand on a recently described model system to explore the contributions of interstrand hydrogen bonds, side-chain/side-chain interactions and side-chain/membrane interactions to beta-sheet formation in membranes. These experiments are based on the observation that the hydrophobic hexapeptide acetyl-Trp-Leu-Leu-Leu-Leu-Leu-OH (AcWLLLLL) folds, cooperatively and reversibly, into oligomeric, antiparallel beta-sheets in phosphatidylcholine membranes. To systematically characterize the important interactions that drive beta-sheet formation in membranes, we have used circular dichroism spectroscopy to determine the membrane secondary structure of each member of a complete host-guest family of related peptides of the form AcWLL-X-LL, where X is one of the natural amino acids. Peptides with hydrophobic X-residues of any size or character (X=Ala, Val, Ile, Leu, Cys, Met, Phe and Trp) form similar beta-sheets in membranes, while peptides with any polar X-residue or Gly or Pro at the X-position are random-coils, even when bound to membranes at high concentrations. The observed membrane sheet preferences correlate poorly with intrinsic sheet propensity scales measured in soluble proteins, but they correlate well with several membrane hydrophobicity scales. These results support the idea that the predominant interactions of the side-chains in membrane-bound beta-sheets are with the membrane lipids, and that backbone hydrogen bonding is the major driving force for the stabilization of beta-sheets in membranes.

Contribution of disulfide bonds to the conformational stability and catalytic activity of ribonuclease A

Disulfide bonds between the side chains of cysteine residues are the only common crosslinks in proteins. Bovine pancreatic ribonuclease A (RNase A) is a 124-residue enzyme that contains four interweaving disulfide bonds (Cys26-Cys84, Cys40-Cys95, Cys58-Cys110, and Cys65-Cys72) and catalyzes the cleavage of RNA. The contribution of each disulfide bond to the conformational stability and catalytic activity of RNase A has been determined by using variants in which each cystine is replaced independently with a pair of alanine residues. Thermal unfolding experiments monitored by ultraviolet spectroscopy and differential scanning calorimetry reveal that wild-type RNase A and each disulfide variant unfold in a two-state process and that each disulfide bond contributes substantially to conformational stability. The two terminal disulfide bonds in the amino-acid sequence (Cys26-Cys84 and Cys58-Cys110) enhance stability more than do the two embedded ones (Cys40-Cys95 and Cys65-Cys72). Removing either one of the terminal disulfide bonds liberates a similar number of residues and has a similar effect on conformational stability, decreasing the midpoint of the thermal transition by almost 40 degrees C. The disulfide variants catalyze the cleavage of poly(cytidylic acid) with values of kcat/Km that are 2- to 40-fold less than that of wild-type RNase A. The two embedded disulfide bonds, which are least important to conformational stability, are most important to catalytic activity. These embedded disulfide bonds likely contribute to the proper alignment of residues (such as Lys41 and Lys66) that are necessary for efficient catalysis of RNA cleavage.

Generation of ligand binding sites in T4 lysozyme by deficiency-creating substitutions

Several variants of T4 lysozyme have been identified that sequester small organic ligands in cavities or clefts. To evaluate potential binding sites for non-polar molecules, we screened a number of hydrophobic large-to-small mutants for stabilization in the presence of benzene. In addition to Leu99-->Ala, binding was indicated for at least five other mutants. Variants Met102-->Ala and Leu133-->Gly, and a crevice mutant, Phe104-->Ala, were further characterized using X-ray crystallography and thermal denaturation. As predicted from the shape of the cavity in the benzene complex, mutant Leu133-->Gly also bound p-xylene. We attempted to enlarge the cavity of the Met102-->Ala mutant into a deep crevice through an additional substitution, but the double mutant failed to bind ligands because an adjacent helix rearranged into a non-helical structure, apparently due to the loss of packing interactions. In general, the protein structure contracted slightly to reduce the volume of the void created by truncating substitutions and expanded upon binding the non-polar ligand, with shifts similar to those resulting from the mutations.A polar molecule binding site was also created by truncating Arg95 to alanine. This creates a highly complementary buried polar environment that can be utilized as a specific "receptor" for a guanidinium ion. Our results suggest that creating a deficiency through truncating mutations of buried residues generates "binding potential" for ligands with characteristics similar to the deleted side-chain. Analysis of complex and apo crystal structures of binding and non-binding mutants suggests that ligand size and shape as well as protein flexibility and complementarity are all determinants of binding. Binding at non-polar sites is governed by hydrophobicity and steric interactions and is relatively permissive. Binding at a polar site is more restrictive and requires extensive complementarity between the ligand and the site.

Energetics of folding and DNA binding of the MAT alpha 2 homeodomain

Homeodomains are a class of DNA-binding protein domains which play an important role in genetic regulation in eukaryotes. We have characterized the thermodynamics of folding and sequence-specific association with DNA of the MAT alpha 2 homeodomain of yeast. Using differential scanning and isothermal titration calorimetry, we measured the enthalpy, heat capacity, and Gibbs free energy changes of these processes. The protein-DNA interaction is enthalpically driven at physiological temperatures. DSC data on the process of melting the protein-DNA complex at different salt concentrations were dissected into its endothermic components, yielding the enthalpy change and dissociation constant of binding. A comparison of the circular dichroism spectra of the free and DNA-bound protein species revealed the formation of additional alpha-helical structure upon binding to DNA. We propose that the latter half of helix 3, the recognition helix, is substantially unfolded in the free protein under the conditions used, as has been observed with other homeodomains [Tsao, D. H. H., et al. (1994) Biochemistry 33, 15053-15060: Cox, M., et al. (1995) J. Biomol. NMR 5, 23-32]. Formation of protein structure is induced by DNA binding, and the energies measured for association therefore include a component due to folding.