Protein phase diagrams II: nonideal behavior of biochemical reactions in the presence of osmolytes

Protein Function Analysis through Machine Learning

Machine learning (ML) has been an important arsenal in computational biology used to elucidate protein function for decades. With the recent burgeoning of novel ML methods and applications, new ML approaches have been incorporated into many areas of computational biology dealing with protein function. We examine how ML has been integrated into a wide range of computational models to improve prediction accuracy and gain a better understanding of protein function. The applications discussed are protein structure prediction, protein engineering using sequence modifications to achieve stability and druggability characteristics, molecular docking in terms of protein–ligand binding, including allosteric effects, protein–protein interactions and protein-centric drug discovery. To quantify the mechanisms underlying protein function, a holistic approach that takes structure, flexibility, stability, and dynamics into account is required, as these aspects become inseparable through their interdependence. Another key component of protein function is conformational dynamics, which often manifest as protein kinetics. Computational methods that use ML to generate representative conformational ensembles and quantify differences in conformational ensembles important for function are included in this review. Future opportunities are highlighted for each of these topics.

"Cooperative collapse" of the denatured state revealed through Clausius-Clapeyron analysis of protein denaturation phase diagrams

Protein phase diagrams have a unique potential to identify the presence of additional thermodynamic states even when non-2-state character is not readily apparent from the experimental observables used to follow protein unfolding transitions. Two-state analysis of the von Willebrand factor A3 domain has previously revealed a discrepancy in the calorimetric enthalpy obtained from thermal unfolding transitions as compared with Gibbs-Helmholtz analysis of free energies obtained from the Linear Extrapolation Method (Tischer and Auton, Prot Sci 2013; 22(9):1147-60). We resolve this thermodynamic conundrum using a Clausius-Clapeyron analysis of the urea-temperature phase diagram that defines how Δ H and the urea m-value interconvert through the slope of c_m versus T, ( ∂ c m / ∂ T ) = Δ H / ( m T ) . This relationship permits the calculation of Δ H at low temperature from m-values obtained through iso-thermal urea denaturation and high temperature m-values from Δ H obtained through iso-urea thermal denaturation. Application of this equation uncovers sigmoid transitions in both cooperativity parameters as temperature is increased. Such residual thermal cooperativity of Δ H and the m-value confirms the presence of an additional state which is verified to result from a cooperative phase transition between urea-expanded and thermally-compact denatured states. Comparison of the equilibria between expanded and compact denatured ensembles of disulfide-intact and carboxyamidated A3 domains reveals that introducing a single disulfide crosslink does not affect the presence of the additional denatured state. It does, however, make a small thermodynamically favorable free energy (∼-13 ± 1 kJ/mol) contribution to the cooperative denatured state collapse transition as temperature is raised and urea concentration is lowered. The thermodynamics of this "cooperative collapse" of the denatured state retain significant compensations between the enthalpy and entropy contributions to the overall free energy.

The N-Terminal Domain of ALS-Linked TDP-43 Assembles without Misfolding

Transactivation response element (TAR) DNA-binding protein 43 (TDP-43) misfolding is implicated in several neurodegenerative diseases characterized by aggregated protein inclusions. Misfolding is believed to be mediated by both the N- and C-terminus of TDP-43; however, the mechanistic basis of the contribution of individual domains in the process remained elusive. Here, using single-molecule fluorescence and ensemble biophysical techniques, and a wide range of pH and temperature conditions, we show that TDP-43^NTD is thermodynamically stable, well-folded and undergoes reversible oligomerization. We propose that, in full-length TDP-43, association between folded N-terminal domains enhances the propensity of the intrinsically unfolded C-terminal domains to drive pathological aggregation.

Synergy in protein-osmolyte mixtures

Virtually all taxa use osmolytes to protect cells against biochemical stress. Osmolytes often occur in mixtures, such as the classical combination of urea with TMAO (trimethylamine N-oxide) in cartilaginous fish or the cocktail of at least six different osmolytes in the kidney. The concentration patterns of osmolyte mixtures found in vivo make it likely that synergy between them plays an important role. Using statistical mechanical n-component Kirkwood-Buff theory, we show from first principles that synergy in protein-osmolyte systems can arise from two separable sources: (1) mutual alteration of protein surface solvation and (2) effects mediated through bulk osmolyte chemical activities. We illustrate both effects in a four-component system with the experimental example of the unfolding of a notch ankyrin domain in urea-TMAO mixtures, which make urea a less effective denaturant and TMAO a more effective stabilizer. Protein surface effects are primarily responsible for this synergy. The specific patterns of surface solvation point to denatured state expansion as the main factor, as opposed to direct competition.

Modulation of allostery by protein intrinsic disorder

Allostery is an intrinsic property of many globular proteins and enzymes that is indispensable for cellular regulatory and feedback mechanisms. Recent theoretical and empirical observations indicate that allostery is also manifest in intrinsically disordered proteins, which account for a substantial proportion of the proteome. Many intrinsically disordered proteins are promiscuous binders that interact with multiple partners and frequently function as molecular hubs in protein interaction networks. The adenovirus early region 1A (E1A) oncoprotein is a prime example of a molecular hub intrinsically disordered protein. E1A can induce marked epigenetic reprogramming of the cell within hours after infection, through interactions with a diverse set of partners that include key host regulators such as the general transcriptional coactivator CREB binding protein (CBP), its paralogue p300, and the retinoblastoma protein (pRb; also called RB1). Little is known about the allosteric effects at play in E1A-CBP-pRb interactions, or more generally in hub intrinsically disordered protein interaction networks. Here we used single-molecule fluorescence resonance energy transfer (smFRET) to study coupled binding and folding processes in the ternary E1A system. The low concentrations used in these high-sensitivity experiments proved to be essential for these studies, which are challenging owing to a combination of E1A aggregation propensity and high-affinity binding interactions. Our data revealed that E1A-CBP-pRb interactions have either positive or negative cooperativity, depending on the available E1A interaction sites. This striking cooperativity switch enables fine-tuning of the thermodynamic accessibility of the ternary versus binary E1A complexes, and may permit a context-specific tuning of associated downstream signalling outputs. Such a modulation of allosteric interactions is probably a common mechanism in molecular hub intrinsically disordered protein function.

Osmolyte-, binding-, and temperature-induced transitions of intrinsically disordered proteins

Structural studies of intrinsically disordered proteins (IDPs) entail unique experimental challenges due in part to the lack of well-defined three-dimensional structures exhibited by this class of proteins. Although IDPs can be studied in their native disordered conformations using a variety of ensemble and single-molecule biophysical techniques, one particularly informative experimental strategy is to probe protein disordered states as part of folding-unfolding transitions. In this chapter, we describe solution methods for probing conformational properties of IDPs (and unfolded proteins, in general), including the use of naturally occurring osmolytes to force protein folding, the quantification of coupled folding and ligand binding of IDPs, and the structural interrogation of solvent- and/or binding-induced folded conformations by thermal perturbations.

Osmolyte effects on protein stability and solubility: a balancing act between backbone and side-chains

In adaptation biology the discovery of intracellular osmolyte molecules that in some cases reach molar levels, raises questions of how they influence protein thermodynamics. We've addressed such questions using the premise that from atomic coordinates, the transfer free energy of a native protein (ΔG(tr,N)) can be predicted by summing measured water-to-osmolyte transfer free energies of the protein's solvent exposed side chain and backbone component parts. ΔG(tr,D) is predicted using a self avoiding random coil model for the protein, and ΔG(tr,D)-ΔG(tr,N), predicts the m-value, a quantity that measures the osmolyte effect on the N⇌D transition. Using literature and newly measured m-values we show 1:1 correspondence between predicted and measured m-values covering a range of 12 kcal/mol/M in protein stability for 46 proteins and 9 different osmolytes. Osmolytes present a range of side chain and backbone effects on N and D solubility and protein stability key to their biological roles.

Protein stability in the presence of cosolutes

Protein scientists have long used cosolutes to study protein stability. While denaturants, such as urea, have been employed for a long time, the attention became focused more recently on protein stabilizers, including osmolytes. Here, we provide practical experimental instructions for the use of both stabilizing and denaturing osmolytes with proteins, as well as data evaluation strategies. We focus on protein stability in the presence of cosolutes and their mixtures at constant and variable temperature.

Single-molecule fluorescence studies of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) (also referred to as natively unfolded proteins) play critical roles in a variety of cellular processes such as transcription and translation and also are linked to several human diseases. Biophysical studies of IDPs present unusual experimental challenges due in part to their broad conformational heterogeneity and potentially complex binding-induced folding behavior. By minimizing the averaging over an ensemble (which is typical of most conventional experiments), single-molecule fluorescence (SMF) techniques have recently begun to add advanced capabilities for structural studies to the experimental arsenal of IDP investigators. Here, we briefly discuss a few common SMF methods that are particularly useful for IDP studies, including SMF resonance energy transfer and fluorescence correlation spectroscopy, along with site-specific protein-labeling methods that are essential for application of these methods to IDPs. We then present an overview of a few studies in this area, highlighting how SMF methods are being used to gain valuable information about two amyloidogenic IDPs, the Parkinson's disease-linked alpha-synuclein and the NM domain of the yeast prion protein Sup 35. SMF experiments provided new information about the proteins' rapidly fluctuating IDP forms, and the complex alpha-synuclein folding behavior upon its binding to lipid and membrane mimics. We anticipate that SMF and single-molecule methods, in general, will find broad application for structural and mechanistic studies of a wide variety of IDPs, both of their disordered conformations, and their ordered ensembles relevant for function and disease.

Changes in thermodynamic stability of von Willebrand factor differentially affect the force-dependent binding to platelet GPIbalpha

In circulation, plasma glycoprotein von Willebrand Factor plays an important role in hemostasis and in pathological thrombosis under hydrodynamic forces. Mutations in the A1 domain of von Willebrand factor cause the hereditary types 2B and 2M von Willebrand disease that either enhance (2B) or inhibit (2M) the interaction of von Willebrand factor with the platelet receptor glycoprotein Ibalpha. To understand how type 2B and 2M mutations cause clinically opposite phenotypes, we use a combination of protein unfolding thermodynamics and atomic force microscopy to assess the effects of two type 2B mutations (R1306Q and I1309V) and a type 2M mutation (G1324S) on the conformational stability of the A1 domain and the single bond dissociation kinetics of the A1-GPIbalpha interaction. At physiological temperature, the type 2B mutations destabilize the structure of the A1 domain and shift the A1-GPIbalpha catch to slip bonding to lower forces. Conversely, the type 2M mutation stabilizes the structure of the A1 domain and shifts the A1-GPIbalpha catch to slip bonding to higher forces. As a function of increasing A1 domain stability, the bond lifetime at low force decreases and the critical force required for maximal bond lifetime increases. Our results are able to distinguish the clinical phenotypes of these naturally occurring mutations from a thermodynamic and biophysical perspective that provides a quantitative description of the allosteric coupling of A1 conformational stability with the force dependent catch to slip bonding between A1 and GPIbalpha.

Investigating protein unfolding kinetics by pulse proteolysis

Investigation of protein unfolding kinetics of proteins in crude samples may provide many exciting opportunities to study protein energetics under unconventional conditions. As an effort to develop a method with this capability, we employed "pulse proteolysis" to investigate protein unfolding kinetics. Pulse proteolysis has been shown to be an effective and facile method to determine global stability of proteins by exploiting the difference in proteolytic susceptibilities between folded and unfolded proteins. Electrophoretic separation after proteolysis allows monitoring protein unfolding without protein purification. We employed pulse proteolysis to determine unfolding kinetics of E. coli maltose binding protein (MBP) and E. coli ribonuclease H (RNase H). The unfolding kinetic constants determined by pulse proteolysis are in good agreement with those determined by circular dichroism. We then determined an unfolding kinetic constant of overexpressed MBP in a cell lysate. An accurate unfolding kinetic constant was successfully determined with the unpurified MBP. Also, we investigated the effect of ligand binding on unfolding kinetics of MBP using pulse proteolysis. On the basis of a kinetic model for unfolding of MBP*maltose complex, we have determined the dissociation equilibrium constant (K(d)) of the complex from unfolding kinetic constants, which is also in good agreement with known K(d) values of the complex. These results clearly demonstrate the feasibility and the accuracy of pulse proteolysis as a quantitative probe to investigate protein unfolding kinetics.

Molecular crowding and solvation: direct and indirect impact on protein reactions

The typical environment for biomolecules in vivo is highly crowded. Under such conditions chemical activities, rather than simply concentrations, govern the behavior of the molecules. In this chapter we discuss the underlying solvation principles that give rise to the chemical activities. We focus on simple experimentally accessible examples, macromolecular crowding, protein folding, and ligand binding under crowded conditions. We discuss effects of high concentrations of both macromolecules and small molecules in terms of the Kirkwood-Buff theory, which couples solution structure to thermodynamics.

Interplay of alpha-synuclein binding and conformational switching probed by single-molecule fluorescence

We studied the coupled binding and folding of alpha-synuclein, an intrinsically disordered protein linked with Parkinson's disease. Using single-molecule fluorescence resonance energy transfer and correlation methods, we directly probed protein membrane association, structural distributions, and dynamics. Results revealed an intricate energy landscape on which binding of alpha-synuclein to amphiphilic small molecules or membrane-like partners modulates conformational transitions between a natively unfolded state and multiple alpha-helical structures. Alpha-synuclein conformation is not continuously tunable, but instead partitions into 2 main classes of folding landscape structural minima. The switch between a broken and an extended helical structure can be triggered by changing the concentration of binding partners or by varying the curvature of the binding surfaces presented by micelles or bilayers composed of the lipid-mimetic SDS. Single-molecule experiments with lipid vesicles of various composition showed that a low fraction of negatively charged lipids, similar to that found in biological membranes, was sufficient to drive alpha-synuclein binding and folding, resulting here in the induction of an extended helical structure. Overall, our results imply that the 2 folded structures are preencoded by the alpha-synuclein amino acid sequence, and are tunable by small-molecule supramolecular states and differing membrane properties, suggesting novel control elements for biological and amyloid regulation of alpha-synuclein.

Effects of cell volume regulating osmolytes on glycerol 3-phosphate binding to triosephosphate isomerase

During cell volume regulation, intracellular concentration changes occur in both inorganic and organic osmolytes in order to balance the extracellular osmotic stress and maintain cell volume homeostasis. Generally, salt and urea increase the Km's of enzymes and trimethylamine N-oxide (TMAO) counteracts these effects by decreasing Km's. The hypothesis to account for these effects is that urea and salt shift the native state ensemble of the enzyme toward conformers that are substrate-binding incompetent (BI), while TMAO shifts the ensemble toward binding competent (BC) species. Km's are often complex assemblies of rate constants involving several elementary steps in catalysis, so to better understand osmolyte effects we have focused on a single elementary event, substrate binding. We test the conformational shift hypothesis by evaluating the effects of salt, urea, and TMAO on the mechanism of binding glycerol 3-phosphate, a substrate analogue, to yeast triosephosphate isomerase. Temperature-jump kinetic measurements promote a mechanism consistent with osmolyte-induced shifts in the [BI]/[BC] ratio of enzyme conformers. Importantly, salt significantly affects the binding constant through its effect on the activity coefficients of substrate, enzyme, and enzyme-substrate complex, and it is likely that TMAO and urea affect activity coefficients as well. Results indicate that the conformational shift hypothesis alone does not account for the effects of osmolytes on Km's.

α-Synuclein Multistate Folding Thermodynamics: Implications for Protein Misfolding and Aggregation

Alpha-synuclein aggregation has been tightly linked with the pathogenesis of Parkinson's disease and other neurodegenerative disorders. Despite the protein's putative function in presynaptic vesicle regulation, the roles of lipid binding in modulating alpha-synuclein conformations and the aggregation process remain to be fully understood. This study focuses on a detailed thermodynamic characterization of monomeric alpha-synuclein folding in the presence of SDS, a well-studied lipid mimetic. Far-UV CD spectroscopy was employed for detection of conformational transitions induced by SDS, temperature, and pH. The data we present here clearly demonstrate the multistate nature of alpha-synuclein folding, which involves two predominantly alpha-helical partially folded thermodynamic intermediates that we designate as F (most folded) and I (intermediately folded) states. Likely structures of these alpha-synuclein conformational states are also discussed. These partially folded forms can exist in the presence of either monomeric or micellar forms of SDS, which suggests that alpha-synuclein has an intrinsic propensity for adopting multiple alpha-helical structures even in the absence of micelle or membrane binding, a feature that may have implications for its biological activity and toxicity. Additionally, we discuss the relation between alpha-synuclein three-state folding and its aggregation, within the context of isothermal titration calorimetry and transmission electron microscopy measurements of SDS-initiated oligomer formation.

Molecular Basis of Osmolyte Effects on Protein and Metabolites

Osmolytes can have strong effects on biochemical reactions, such as protein folding or protein-ligand interaction. These effects are mediated through solvation-the nonspecific interaction between the solution components. Therefore, understanding the impact of osmolytes on cellular biochemistry requires an understanding of the underlying solvation processes. This chapter discusses the thermodynamic effects of osmolytes on proteins and small organic molecules in terms of the solvation of these molecules, as derived from Kirkwood-Buff theory. This approach allows experimental determination of solvation properties from thermodynamic data. Knowledge of solvation at this level provides insight into the observed behavior of proteins and small molecules in osmolyte solution on a microscopic level. As examples, we provide solvation effects on protein folding, ligand binding, and osmolyte thermodynamics.