Thermodynamics of structural stability and cooperative folding behavior in proteins

Collapse and Protein Folding: Should We Be Surprised That Biothermodynamics Works So Well?

A complete understanding of protein function and dynamics requires the characterization of the multiple thermodynamic states, including the denatured state ensemble (DSE). Whereas residual structure in the DSE (as well as in partially folded states) is pertinent in many biological contexts, here we are interested in how such structure affects protein thermodynamics. We examine issues related to chain collapse in light of new developments, focusing on potential complications arising from differences in the DSE's properties under various conditions. Despite some variability in the degree of collapse and structure in the DSE, stability measurements are remarkably consistent between two standard methods, calorimetry and chemical denaturation, as well as with hydrogen-deuterium exchange. This robustness is due in part to the DSEs obtained with different perturbations being thermodynamically equivalent and hence able to serve as a common reference state. An examination of the properties of the DSE points to it as being a highly expanded ensemble with minimal amounts of stable hydrogen bonded structure. These two features are likely to be critical in the broad and successful application of thermodynamics to protein folding. Our review concludes with a discussion of the impact of these findings on folding mechanisms and pathways.

The homeodomain regulates stable DNA binding of prostate cancer target ONECUT2

The CUT and homeodomain are ubiquitous DNA binding elements often tandemly arranged in multiple transcription factor families. However, how the CUT and homeodomain work concertedly to bind DNA remains unknown. Using ONECUT2, a driver and therapeutic target of advanced prostate cancer, we show that while the CUT initiates DNA binding, the homeodomain thermodynamically stabilizes the ONECUT2-DNA complex through allosteric modulation of CUT. We identify an arginine pair in the ONECUT family homeodomain that can adapt to DNA sequence variations. Base interactions by this ONECUT family-specific arginine pair as well as the evolutionarily conserved residues are critical for optimal DNA binding and ONECUT2 transcriptional activity in a prostate cancer model. The evolutionarily conserved base interactions additionally determine the ONECUT2-DNA binding energetics. These findings provide insights into the cooperative DNA binding by CUT-homeodomain proteins.

A cell-free antigen processing system informs HIV-1 epitope selection and vaccine design

Distinct CD4+ T cell epitopes have been associated with spontaneous control of HIV-1 replication, but analysis of antigen-dependent factors that influence epitope selection is lacking. To examine these factors, we used a cell-free antigen processing system that incorporates soluble HLA-DR (DR1), HLA-DM (DM), cathepsins, and full-length protein antigens for epitope identification by LC-MS/MS. HIV-1 Gag, Pol, Env, Vif, Tat, Rev, and Nef were examined using this system. We identified 35 novel epitopes, including glycopeptides. Epitopes from smaller HIV-1 proteins mapped to regions of low protein stability and higher solvent accessibility. HIV-1 antigens associated with limited CD4+ T cell responses were processed efficiently, while some protective epitopes were inefficiently processed. 55% of epitopes obtained from cell-free processing induced memory CD4+ T cell responses in HIV-1+ donors, including eight of 19 novel epitopes tested. Thus, an in vitro processing system utilizing the components of Class II processing reveals factors influencing epitope selection of HIV-1 and represents an approach to understanding epitope selection from non-HIV-1 antigens.

MCL-1 promiscuity and the structural resilience of its binding partners

The allosteric protein MCL-1 and its natural inhibitors, the BH3-only proteins PUMA, BIM, and NOXA regulate apoptosis by interacting promiscuously within an entangled binding network. Little is known about the transient processes and dynamic conformational fluctuations that are the basis for the formation and stability of the MCL-1/BH3-only complex. In this study, we designed photoswitchable versions of MCL-1/PUMA and MCL-1/NOXA, and investigated the protein response after an ultrafast photo-perturbation with transient infrared spectroscopy. We observed partial α-helical unfolding in all cases, albeit on strongly varying timescales (1.6 ns for PUMA, 9.7 ns for the previously studied BIM, and 85 ns for NOXA). These differences are interpreted as a BH3-only-specific "structural resilience" to defy the perturbation while remaining in MCL-1's binding pocket. Thus, the presented insights could help to better understand the differences between PUMA, BIM, and NOXA, the promiscuity of MCL-1, in general, and the role of the proteins in the apoptotic network.

Recombinant pseudorabies virus (PRV) expressing stabilized E2 of classical swine fever virus (CSFV) protects against both PRV and CSFV

Pseudorabies (PR) and classical swine fever (CSF) are economically important infectious diseases of pigs. Most pig farms in China are immunized against these two diseases. Here, we describe a stabilized E2 protein as an immunogen inserted into the PRV genome as a bivalent live virus-vectored vaccine. The E2 protein has 48 variant sites, there are 2-5 candidate amino acids per variant site, and the relative energy contribution of each amino acid to E2 energy was calculated. Combined substitutions of amino acids at the neighbor variant site (neighbor substitution) were performed to obtain the E2 protein sequence with the lowest energy (stabilized E2). Multiple amino acid substitutions at 48 variant sites were performed, and the results were consistent with neighbor substitutions. The stabilized E2 sequence was obtained, and its energy decreased by 22 Rosetta Energy Units (REUs) compared with the original sequence. After the recombinant PRV expressing stabilized E2 of CSFV was constructed, the secretion efficiency of stabilized E2 was increased by 2.97 times, and the thermal stability was increased by 10.5 times. Immunization of mice resulted in a 2-fold increase in antibody production, and a balanced antibody level against subtype 1.1 and subtype 2.1d E2 was achieved. In rabbits immunized, the lethal challenge of PRV-ZJ and the fever response induced by CSFV could be prevented simultaneously. These findings suggest that rPRV-muta/287aaE2 is a promising bivalent vaccine against CSFV and PRV infections.

Pre-Molten, Wet, and Dry Molten Globules en Route to the Functional State of Proteins

Transitions between the unfolded and native states of the ordered globular proteins are accompanied by the accumulation of several intermediates, such as pre-molten globules, wet molten globules, and dry molten globules. Structurally equivalent conformations can serve as native functional states of intrinsically disordered proteins. This overview captures the characteristics and importance of these molten globules in both structured and intrinsically disordered proteins. It also discusses examples of engineered molten globules. The formation of these intermediates under conditions of macromolecular crowding and their interactions with nanomaterials are also reviewed.

Propensities of Some Amino Acid Pairings in α-Helices Vary with Length

The results of secondary structure prediction methods are widely used in applications in biotechnology and bioinformatics. However, the accuracy limit of these methods could be improved up to 92%. One approach to achieve this goal is to harvest information from the primary structure of the peptide. This study aims to contribute to this goal by investigating the variations in propensity of amino acid pairings to α-helices in globular proteins depending on helix length. (n):(n + 4) residue pairings were determined using a comprehensive peptide data set according to backbone hydrogen bond criterion which states that backbone hydrogen bond is the dominant driving force of protein folding. Helix length is limited to 13 to 26 residues. Findings of this study show that propensities of ALA:GLY and GLY:GLU pairings to α-helix in globular protein increase with increasing helix length but of ALA:ALA and ALA:VAL decrease. While the frequencies of ILE:ALA, LEU:ALA, LEU:GLN, LEU:GLU, LEU:LEU, MET:ILE and VAL:LEU pairings remain roughly constant with length, the 25 residue pairings have varying propensities in narrow helix lengths. The remaining pairings have no prominent propensity to α-helices.

Thermodynamics for the Self-Assembly of Alkylated Peptides

Self-assembling peptides form aggregates with various nanostructures such as spheres, sheets, and fibers and have potential applications in nanomedicine and drug delivery. The alkylation of peptides is a promising strategy for controlling the self-assembly of peptides. In this study, we investigated the thermodynamic properties associated with the aggregation of alkyl-chain-modified self-assembling peptides. The tripeptide sequence, KYF, which has been reported to form fibrous aggregates via self-assembly, was modified with various fatty acids at the N-terminus. The fibrous morphology of the aggregates was observed by transmission electron microscopy and atomic force microscopy. Thioflavin T fluorescence and circular dichroism spectroscopy revealed the formation of β-sheet structures. The critical micelle concentration and its temperature dependence were determined to obtain the thermodynamic parameters for aggregation. The results showed that the aggregation was an entropy-driven process at low temperatures, whereas it was enthalpy-driven at high temperatures. The negative heat capacity changes for aggregation suggested that hydrophobic interactions were the major driving force for self-assembly. Other entropic and enthalpic interactions were also contributed in part to the self-assembly. We individually identified the contributions of the peptide and alkyl chain moiety to the self-assembly. These contributions can be explained by the theoretical values for the self-assembly of each component. The results of this study provide fundamental insights into the design of self-associating peptides.

Intrinsic affinity of protein - ligand binding by differential scanning calorimetry

Differential scanning calorimetry (DSC) determines the enthalpy change upon protein unfolding and the melting temperature of the protein. Performing DSC of a protein in the presence of increasing concentrations of specifically-binding ligand yields a series of curves that can be fit to obtain the protein-ligand dissociation constant as done in the fluorescence-based thermal shift assay (FTSA, ThermoFluor, DSF). The enthalpy of unfolding, as directly determined by DSC, helps improving the precision of the fit. If the ligand binding is linked to protonation reactions, the intrinsic binding constant can be determined by performing the affinity determination at a series of pH values. Here, the intrinsic, pH-independent, affinity of acetazolamide binding to carbonic anhydrase (CA) II was determined. A series of high-affinity ligands binding to CAIX, an anticancer drug target, and CAII showed recognition and selectivity for the anticancer isozyme. Performing the DSC experiment in buffers of highly different enthalpies of protonation enabled to observe the ligand unbinding-linked protonation reactions and estimate the intrinsic enthalpy of binding. The heat capacity of combined unfolding and unbinding was determined by varying the ligand concentrations. Taken together, these parameters provided a detailed thermodynamic picture of the linked ligand binding and protein unfolding process.

Altered Conformational Landscape upon Sensing Guanine Nucleotides in a Disease Mutant of Elongation Factor-like 1 (EFL1) GTPase

The final maturation step of the 60S ribosomal subunit requires the release of eukaryotic translation initiation factor 6 (human eIF6, yeast Tif6) to enter the pool of mature ribosomes capable of engaging in translation. This process is mediated by the concerted action of the Elongation Factor-like 1 (human EFL1, yeast Efl1) GTPase and its effector, the Shwachman-Bodian-Diamond syndrome protein (human SBDS, yeast Sdo1). Mutations in these proteins prevent the release of eIF6 and cause a disease known as Shwachman–Diamond Syndrome (SDS). While some mutations in EFL1 or SBDS result in insufficient proteins to meet the cell production of mature large ribosomal subunits, others do not affect the expression levels with unclear molecular defects. We studied the functional consequences of one such mutation using Saccharomyces cerevisiae Efl1 R1086Q, equivalent to human EFL1 R1095Q described in SDS patients. We characterised the enzyme kinetics and energetic basis outlining the recognition of this mutant to guanine nucleotides and Sdo1, and their interplay in solution. From our data, we propose a model where the conformational change in Efl1 depends on a long-distance network of interactions that are disrupted in mutant R1086Q, whereby Sdo1 and the guanine nucleotides no longer elicit the conformational changes previously described in the wild-type protein. These findings point to the molecular malfunction of an EFL1 mutant and its possible impact on SDS pathology.

Molecular understanding of calorimetric protein unfolding experiments

Testing and predicting protein stability gained importance because proteins, including antibodies, became pharmacologically relevant in viral and cancer therapies. Isothermal scanning calorimetry is the principle method to study protein stability. Here, we use the excellent experimental heat capacity C_p(T) data from the literature for a critical inspection of protein unfolding as well as for the test of a new cooperative model. In the relevant literature, experimental temperature profiles of enthalpy, H_cal(T), entropy, S_cal(T), and free energy, G_cal(T) are missing. First, we therefore calculate the experimental H_cal(T), S_cal(T), and G_cal(T) from published C_p(T) thermograms. Considering only the unfolding transition proper, the heat capacity and all thermodynamic functions are zero in the region of the native protein. In particular, the free energy of the folded proteins is also zero and G_cal(T) displays a trapezoidal temperature profile when cold denaturation is included. Second, we simulate the DSC-measured thermodynamic properties with a new molecular model based on statistical-mechanical thermodynamics. The model quantifies the protein cooperativity and predicts the aggregate thermodynamic variables of the system with molecular parameters only. The new model provides a perfect simulation of all thermodynamic properties, including the observed trapezoidal G_cal(T) temperature profile. Importantly, the new cooperative model can be applied to a broad range of protein sizes, including antibodies. It predicts not only heat and cold denaturation but also provides estimates of the unfolding kinetics and allows a comparison with molecular dynamics calculations and quasielastic neutron scattering experiments.

Investigating single amino acid substitutions in PIM1 kinase: A structural genomics approach

PIM1, is a serine/threonine proto-oncogene kinase, involved in many biological functions, including cell survival, proliferation, and differentiation, thus play a key role in oncogenesis. It plays a crucial role in the onset and progression of various hematopoietic and non-hematopoietic malignancies, including acute myeloid leukemia and prostate cancer. Mutations in PIM1, especially in its kinase domain, can induce abnormal structural changes and thus alter functionalities that can lead to disease progression and other complexities. Herein, we have performed an extensive analysis of the PIM1 mutations at sequence and structure level while utilizing state-of-the-art computational approaches. Based on the impact on PIM1, numerous pathogenic and destabilizing mutations were identified and subsequently analyzed in detail. Finally, two amino acid substitutions (W109C and F147C) in the kinase domain of PIM1 were selected to explore their impact on the PIM1 structure in a time evolution manner using all-atom molecular dynamics (MD) simulations for 200 ns. MD results indicate significant conformational altercations in the structure of PIM1, especially upon F147C mutation. This study provides a significant insight into the PIM1 dysfunction upon single amino acid substitutions, which can be utilized to get insights into the molecular basis of PIM1-associated disease progression.

Quantification of Conformational Entropy Unravels Effect of Disordered Flanking Region in Coupled Folding and Binding

Intrinsic disorder (ID) constitutes a new dimension to the protein structure-function relationship. The ability to undergo conformational changes upon binding is a key property of intrinsically disordered proteins and remains challenging to study using conventional methods. A 1994 paper by R. S. Spolar and M. T. Record presented a thermodynamic approach for estimating changes in conformational entropy based on heat capacity changes, allowing quantification of residues folding upon binding. Here, we adapt the method for studies of intrinsically disordered proteins. We integrate additional data to provide a broader experimental foundation for the underlying relations and, based on >500 protein-protein complexes involving disordered proteins, reassess a key relation between polar and nonpolar surface area changes, previously determined using globular protein folding. We demonstrate the improved suitability of the adapted method to studies of the folded αα-hub domain RST from radical-induced cell death 1, whose interactome is characterized by ID. From extensive thermodynamic data, quantifying the conformational entropy changes upon binding, and comparison to the NMR structure, the adapted method improves accuracy for ID-based studies. Furthermore, we apply the method, in conjunction with NMR, to reveal hitherto undetected effects of interaction-motif context. Thus, inclusion of the disordered context of the DREB2A RST-binding motif induces structuring of the binding motif, resulting in major enthalpy-entropy compensation in the interaction interface. This study, also evaluating additional interactions, demonstrates the strength of the ID-adapted Spolar-Record thermodynamic approach for dissection of structural features of ID-based interactions, easily overlooked in traditional studies, and for translation of these into mechanistic knowledge.

Stability of the Retinoid X Receptor-α Homodimer in the Presence and Absence of Rexinoid and Coactivator Peptide

Differential scanning calorimetry and differential scanning fluorimetry were used to measure the thermal stability of human retinoid X receptor-α ligand binding domain (RXRα LBD) homodimer in the absence or presence of rexinoid and coactivator peptide, GRIP-1. The apo-RXRα LBD homodimer displayed a single thermal unfolding transition with a Tm of 58.7 °C and an unfolding enthalpy (ΔH) of 673 kJ/mol (12.5 J/g), much lower than average value (35 J/g) of small globular proteins. Using a heat capacity change (ΔCp) of 15 kJ/(mol K) determined by measurements at different pH values, the free energy of unfolding (ΔG) of the native state was 33 kJ/mol at 37 °C. Rexinoid binding to the apo-homodimer increased Tm by 5 to 9 °C and increased the ΔG of the native homodimer by 12 to 20 kJ/mol at 37 °C, consistent with the nanomolar dissociation constant (Kd) of the rexinoids. GRIP-1 binding to holo-homodimers containing rexinoid resulted in additional increases in ΔG of 14 kJ/mol, a value that was the same for all three rexinoids. Binding of rexinoid and GRIP-1 resulted in a combined 50% increase in unfolding enthalpy, consistent with reduced structural fluidity and more compact folding observed in other published structural studies. The complexes of UAB110 and UAB111 are each more stable than the UAB30 complex by 8 kJ/mol due to enhanced hydrophobic interactions in the binding pocket because of their larger end groups. This increase in thermodynamic stability positively correlates with their improved RXR activation potency. Thermodynamic measurements are thus valuable in predicting agonist potency.

DNA-binding function of c-Myb R2R3 around thermal denaturation temperature

The minimum DNA-binding domain of the transcrip­tional factor c-Myb R2R3 remarkably fluctuates in the solution. In the present study, we evaluated the protein fluctuation of R2R3 C130I mutant, R2R3*, on its DNA-binding and folding thermodynamics. DNA-binding analysis using isothermal titration calorimetry revealed that the heat capacity change determined from the correlation between temperature and binding enthalpy change is highly negative above 35°C, indicating that the fluctuation increases with increasing temperature and elevates the conformational change on DNA binding. The results were in accordance with those of differential scanning calorimetry, which revealed that the heat capacity corresponding to thermal denatu­ration gradually increased above 35°C, followed by the broad transition peak. In contrast, the transition peak of R2R3* in the DNA-bound state was sharper and larger than that in the DNA-unbound state. The fluctuating form could transform into lesser fluctuating form upon DNA binding, resulting in a larger enthalpy change for denaturation of R2R3* in the DNA-bound state. It should also be noted that R2R3* could specifi­cally bind to DNA around thermal denaturation temperature. This would be due to proteins with numerous fluctuations. Moreover, we discuss specific and non-specific DNA binding accompanied by the conformational change between well-ordered and disordered forms of R2R3* observed around the denaturation temperature.

Halogen Atoms in the Protein-Ligand System. Structural and Thermodynamic Studies of the Binding of Bromobenzotriazoles by the Catalytic Subunit of Human Protein Kinase CK2

Binding of a family of brominated benzotriazoles to the catalytic subunit of human protein kinase CK2 (hCK2α) was used as a model system to assess the contribution of halogen bonding to protein–ligand interaction. CK2 is a constitutively active pleiotropic serine/threonine protein kinase that belongs to the CMGC group of eukaryotic protein kinases (EPKs). Due to the addiction of some cancer cells, CK2 is an attractive and well-characterized drug target. Halogenated benzotriazoles act as ATP-competitive inhibitors with unexpectedly good selectivity for CK2 over other EPKs. We have characterized the interaction of bromobenzotriazoles with hCK2α by X-ray crystallography, low-volume differential scanning fluorimetry, and isothermal titration calorimetry. Properties of free ligands in solution were additionally characterized by volumetric and RT-HPLC measurements. Thermodynamic data indicate that the affinity increases with bromo substitution, with greater contributions from 5- and 6-substituents than 4- and 7-substituents. Except for 4,7-disubstituted compounds, the bromobenzotriazoles adopt a canonical pose with the triazole close to lysine 68, which precludes halogen bonding. More highly substituted benzotriazoles adopt many additional noncanonical poses, presumably driven by a large hydrophobic contribution to binding. Some noncanonical ligand orientations allow the formation of halogen bonds with the hinge region. Consistent with a predominantly hydrophobic interaction, the isobaric heat capacity decreases upon ligand binding, the more so the higher the substitution.

Structural and Energetic Characterization of the Denatured State from the Perspectives of Peptides, the Coil Library, and Intrinsically Disordered Proteins

The α and polyproline II (PPII) basins are the two most populated regions of the Ramachandran map when constructed from the protein coil library, a widely used denatured state model built from the segments of irregular structure found in the Protein Data Bank. This indicates the α and PPII conformations are dominant components of the ensembles of denatured structures that exist in solution for biological proteins, an observation supported in part by structural studies of short, and thus unfolded, peptides. Although intrinsic conformational propensities have been determined experimentally for the common amino acids in short peptides, and estimated from surveys of the protein coil library, the ability of these intrinsic conformational propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs), an increasingly important class of proteins in cell function, has thus far proven elusive to establish. Recently, we demonstrated that the sequence dependence of the mean hydrodynamic size of IDPs in water and the impact of heat on the coil dimensions, provide access to both the sequence dependence and thermodynamic energies that are associated with biases for the α and PPII backbone conformations. Here, we compare results from peptide-based studies of intrinsic conformational propensities and surveys of the protein coil library to those of the sequence-based analysis of heat effects on IDP hydrodynamic size, showing that a common structural and thermodynamic description of the protein denatured state is obtained.

Heat Capacity Changes (ΔCp) for Interconversions between Differentially-Ordered DNA States within Physiological Temperature Domains: Implications for Biological Regulatory Switches

Knowledge of differences in heat capacity changes (ΔC_p) between biopolymer states provides essential information about the temperature dependence of the thermodynamic properties of these states, while also revealing insights into the nature of the forces that drive the formation of functional and dysfunctional biopolymer "order." In contrast to proteins, for nucleic acids there is a dearth of direct experimental determination of this information-rich parameter, a deficiency that compromises interpretations of the ever-increasing thermodynamic analyses of nucleic acid properties; particularly as they relate to differential nucleic acid (meta)stability states and their potential biological functions. Here we demonstrate that such heat capacity differences, in fact, exist not only between traditionally measured native to fully unfolded (assumed "random coil") DNA states, but also between competing order-to-order transformations. We illustrate the experimental approach by measuring the heat capacity change between "native"/ordered, sequence homologous, "isomeric" DNA states that differ in conformation but not sequence. Importantly, these heat capacity differences occur within biologically relevant temperature ranges. In short, we describe a new and general method to measure the value of such heat capacity differences anywhere in experimentally accessible conformational and temperature space; in this case, between two metastable bulge loop states, implicated in DNA expansion diseases, and their competing, fully paired, thermodynamically more stable duplex states. This measurement reveals a ΔC_p of 61 ± 7 cal mol_bp ^-1 K ^-1. Such heat capacity differences between competing DNA "native" ensemble states must be considered when evaluating equilibria between different DNA "ordered" conformations, including the assessment of the differential stabilizing forces and potential biological functions of competing DNA "structured" motifs.

The structure of the extracellular domains of human interleukin 11α receptor reveals mechanisms of cytokine engagement

Interleukin (IL) 11 activates multiple intracellular signaling pathways by forming a complex with its cell surface α-receptor, IL-11Rα, and the β-subunit receptor, gp130. Dysregulated IL-11 signaling has been implicated in several diseases, including some cancers and fibrosis. Mutations in IL-11Rα that reduce signaling are also associated with hereditary cranial malformations. Here we present the first crystal structure of the extracellular domains of human IL-11Rα and a structure of human IL-11 that reveals previously unresolved detail. Disease-associated mutations in IL-11Rα are generally distal to putative ligand-binding sites. Molecular dynamics simulations showed that specific mutations destabilize IL-11Rα and may have indirect effects on the cytokine-binding region. We show that IL-11 and IL-11Rα form a 1:1 complex with nanomolar affinity and present a model of the complex. Our results suggest that the thermodynamic and structural mechanisms of complex formation between IL-11 and IL-11Rα differ substantially from those previously reported for similar cytokines. This work reveals key determinants of the engagement of IL-11 by IL-11Rα that may be exploited in the development of strategies to modulate formation of the IL-11-IL-11Rα complex.

Consensus protein engineering on the thermostable histone-like bacterial protein HUs significantly improves stability and DNA binding affinity

Consensus-based protein engineering strategy has been applied to various proteins and it can lead to the design of proteins with enhanced biological performance. Histone-like HUs comprise a protein family with sequence variety within a highly conserved 3D-fold. HU function includes compacting and regulating bacterial DNA in a wide range of biological conditions in bacteria. To explore the possible impact of consensus-based design in the thermodynamic stability of HU proteins, the approach was applied using a dataset of sequences derived from a group of 40 mesostable, thermostable, and hyperthermostable HUs. The consensus-derived HU protein was named HUBest, since it is expected to perform best. The synthetic HU gene was overexpressed in E. coli and the recombinant protein was purified. Subsequently, HUBest was characterized concerning its correct folding and thermodynamic stability, as well as its ability to interact with plasmid DNA. A substantial increase in HUBest stability at high temperatures is observed. HUBest has significantly improved biological performance at ambience temperature, presenting very low K_d values for binding plasmid DNA as indicated from the Gibbs energy profile of HUBest. This K_d may be associated to conformational changes leading to decreased thermodynamic stability and, therefore, higher flexibility at ambient temperature.

Impact of Heat on Coil Hydrodynamic Size Yields the Energetics of Denatured State Conformational Bias

Conformational equilibria in the protein denatured state have key roles regulating folding, stability, and function. The extent of conformational bias in the protein denatured state under folding conditions, however, has thus far proven elusive to quantify, particularly with regard to its sequence dependence and energetic character. To better understand the structural preferences of the denatured state, we analyzed both the sequence dependence to the mean hydrodynamic size of disordered proteins in water and the impact of heat on the coil dimensions, showing that the sequence dependence and thermodynamic energies associated with intrinsic biases for the α and polyproline II (PPII) backbone conformations can be obtained. Experiments that evaluate how the hydrodynamic size changes with compositional changes in the protein reveal amino acid specific preferences for PPII that are in good quantitative agreement with calorimetry-measured values from unfolded peptides and those inferred by survey of the protein coil library. At temperatures above 25 °C, the denatured state follows the predictions of a PPII-dominant ensemble. Heat effects on coil hydrodynamic size indicate the α bias is comparable to the PPII bias at cold temperatures. Though historically thought to give poor resolution to structural details, the hydrodynamic size of the unfolded state is found to be an effective reporter on the extent of the biases for the α and PPII backbone conformations.

Thermodynamics of DNA: heat capacity changes on duplex unfolding

The heat capacity change, ΔCp, accompanying the folding/unfolding of macromolecules reflects their changing state of hydration. Thermal denaturation of the DNA duplex is characterized by an increase in ΔCp but of much lower magnitude than observed for proteins. To understand this difference, the changes in solvent accessible surface area (ΔASA) have been determined for unfolding the B-form DNA duplex into disordered single strands. These showed that the polar component represents ~ 55% of the total increase in ASA, in contrast to globular proteins of similar molecular weight for which the polar component is only about 1/3rd of the total. As the exposure of polar surface results in a decrease of ΔCp, this explains the much reduced heat capacity increase observed for DNA and emphasizes the enhanced role of polar interactions in maintaining duplex structure. Appreciation of a non-zero ΔCp for DNA has important consequences for the calculation of duplex melting temperatures (T_m). A modified approach to T_m prediction is required and comparison is made of current methods with an alternative protocol.

Sequence Reversal Prevents Chain Collapse and Yields Heat-Sensitive Intrinsic Disorder

Sequence patterns of charge, hydrophobicity, hydrogen bonding, and other amino acid physicochemical properties contribute to mechanisms of protein folding, but how sequence composition and patterns influence the conformational dynamics of the denatured state ensemble is not fully understood. To investigate structure-sequence relationships in the denatured state, we reversed the sequence of staphylococcal nuclease and characterized its structure, thermodynamic character, and hydrodynamic radius using circular dichroism spectroscopy, dynamic light scattering, analytical ultracentrifugation, and size-exclusion chromatography as a function of temperature. The macromolecular size of "Retro-nuclease" is highly expanded in solution with characteristics similar to biological intrinsically disordered proteins. In contradistinction to a disordered state, Retro-nuclease exhibits a broad sigmoid transition of its hydrodynamic dimensions as temperature is increased, indicating a thermodynamically controlled compaction. Counterintuitively, the magnitude of these temperature-induced hydrodynamic changes exceed that observed from thermal denaturation of folded unaltered staphylococcal nuclease. Undetectable by calorimetry and intrinsic tryptophan fluorescence, the lack of heat capacity or fluorescence changes throughout the thermal transition indicate canonical hydrophobic collapse did not drive the Retro-nuclease structural transitions. Temperature-dependent circular dichroism spectroscopy performed on Retro-nuclease and computer simulations correlate to temperature sensitivity in the intrinsic sampling of backbone conformations for polyproline II and α-helix. The experimental results indicate a role for sequence direction in mediating the collapse of the polypeptide chain, whereas the simulation trends illustrate the generality of the observed heat effects on disordered protein structure.

Thermodynamic Stability of the Transcription Regulator PaaR2 from Escherichia coli O157:H7

PaaR2 is a putative transcription regulator encoded by a three-component parDE-like toxin-antitoxin module from Escherichia coli O157:H7. Although this module's toxin, antitoxin, and toxin-antitoxin complex have been more thoroughly investigated, little remains known about its transcription regulator PaaR2. Using a wide range of biophysical techniques (circular dichroism spectroscopy, size-exclusion chromatography-multiangle laser light scattering, dynamic light scattering, small-angle x-ray scattering, and native mass spectrometry), we demonstrate that PaaR2 mainly consists of α-helices and displays a concentration-dependent octameric build-up in solution and that this octamer contains a global shape that is significantly nonspherical. Thermal unfolding of PaaR2 is reversible and displays several transitions, suggesting a complex unfolding mechanism. The unfolding data obtained from spectroscopic and calorimetric methods were combined into a unifying thermodynamic model, which suggests a five-state unfolding trajectory. Furthermore, the model allows the calculation of a stability phase diagram, which shows that, under physiological conditions, PaaR2 mainly exists as a dimer that can swiftly oligomerize into an octamer depending on local protein concentrations. These findings, based on a thorough biophysical and thermodynamic analysis of PaaR2, may provide important insights into biological function such as DNA binding and transcriptional regulation.

Structural, kinetic and thermodynamic characterizations of SDS-induced molten globule state of a highly negatively charged cytochrome c

This study presents the structural, kinetic and thermodynamic characterizations of previously unknown submicellar concentrations of SDS-induced molten globule (MGSDS) state of a highly negatively charged base-denatured ferricytochrome c (UB-state) at pH ∼12.8 (±0.2). The far-UV CD, near-UV CD, ANS-fluorescence data of UB-state in the presence of different concentrations of SDS indicate that the submicellar concentrations of SDS (≤0.4 mM) transform the UB-state to MGSDS-state. The MGSDS-state has native-like α-helical secondary structure but lacks tertiary structure. The free energy change (ΔG°D) for UB→ MGSDS transition determined by far-UV CD (∼2.7 kcal mol-1) is slightly higher than those determined by fluorescence (∼2.0 kcal mol-1) at 25°C. At very low SDS and NaCl concentrations, the MGSDS-state undergoes cold denaturation. As SDS concentration is increased, the thermal denaturation temperature increases and the cold denaturation temperature decrease. Kinetic experiments involving the measurement of the CO-association rate to the base-denatured ferrocytochrome c at pH ≈12.8 (±0.2), 25°C indicate that the submicellar concentrations of SDS restrict the internal dynamics of base-denatured protein.

Cooperative energetic effects elicited by the yeast Shwachman-Diamond syndrome protein (Sdo1) and guanine nucleotides modulate the complex conformational landscape of the elongation factor-like 1 (Efl1) GTPase

One of the final maturation steps of the large ribosomal subunit requires the joint action of the elongation factor-like 1 (human EFL1, yeast Efl1) GTPase and the Shwachman-Diamond syndrome protein (human SBDS, yeast Sdo1) to release the eukaryotic translation initiation factor 6 (human eIF6, yeast Tif6) and allow the assembly of mature ribosomes. EFL1 function is driven by conformational changes. However, the nature of such conformational changes or the mechanism by which they are prompted are still largely unknown. In previous studies, it has been established that this GTPase interacts with its cofactor in solution in an inverted orientation with respect to the binding mode derived from 60S ribosome subunit cryo-EM data. To shed new light on this conundrum, we characterized calorimetrically the energetic basis describing the recognition of Efl1 to GT(D)P, Sdo1 and their intercommunication in solution. A structural-based analysis of the binding signatures indicates that Efl1 has a large structural flexibility. The mutual effects of Sdo1 and nucleotides on Efl1 modulate in a very specific and robust way the complex conformational landscape of Efl1, resembling the behavior observed with other GTPases and their cofactors.

Hydration differences between the major and minor grooves of DNA revealed from heat capacity measurements

The nature of water on the surface of a macromolecule is reflected in the temperature dependence of the heat effect, i.e., the heat capacity change, ΔCp, that accompanies its removal on forming a complex. The relationship between ΔCp and the nature of the surface dehydrated cannot be modeled for DNA by the use of small molecules, as previously done for proteins, since the contiguous surfaces of the grooves cannot be treated as the sum of small component molecules such as nucleotides. An alternative approach is used here in which ΔCp is measured for the formation of several protein/DNA complexes and the calculated contribution from protein dehydration subtracted to yield the heat capacity change attributable to dehydration of the DNA. The polar and apolar surface areas of the DNA dehydrated on complex formation were calculated from the known structures of the complexes, allowing heat capacity coefficients to be derived representing dehydration of unit surface area of polar and apolar surface in both grooves. Dehydration of apolar surfaces in both grooves is essentially identical and accompanied by a reduction in ΔCp by about 3 J K⁻¹ mol⁻¹ (Ų)⁻¹, a value of somewhat greater magnitude than observed for proteins {ΔCp = − 1.79 J K⁻¹ mol⁻¹ (Ų)⁻¹}. In contrast, dehydration of polar surfaces is very different in the two grooves: in the minor groove ΔCp increases by 2.7 J K⁻¹ mol⁻¹ (Ų)⁻¹, but in the major groove, although ΔCp is also positive, it is low in value: + 0.4 J K⁻¹ mol⁻¹ (Ų)⁻¹. Physical explanations for the magnitudes of ΔCp are discussed.

The Stabilities of the Soluble Ectodomain and Fusion Peptide Hairpins of the Influenza Virus Hemagglutinin Subunit II Protein Are Positively Correlated with Membrane Fusion

Cellular entry of influenza virus is mediated by the viral protein hemagglutinin (HA), which forms an initial complex of three HA1 and three HA2 subunits. Each HA2 includes a fusion peptide (FP), a soluble ectodomain (SE), and a transmembrane domain. HA1 binds to cellular sialic acids, followed by virus endocytosis, pH reduction, dissociation of HA1, and structural rearrangement of HA2 into a final trimer-of-SE hairpins. A decrease in pH also triggers HA2-mediated virus/endosome membrane fusion. SE hairpins have an interior parallel helical bundle and C-terminal strands in the grooves of the exterior of the bundle. FPs are separate helical hairpins. This study compares wild-type HA2 (WT-HA2) with G1E(FP) and I173E(SE strand) mutants. WT-HA2 induces vesicle fusion at pH 5.0, whereas the extent of fusion is greatly reduced for both mutants. Circular dichroism for HA2 and FHA2≡FP+SE constructs shows dramatic losses of stability for the mutants, including a T_m reduced by 40 °C for I173E-FHA2. This is evidence of destabilization of SE hairpins via dissociation of strands from the helical bundle, which is also supported by larger monomer fractions for mutant versus WT proteins. The G1E mutant may have disrupted FP hairpins, with consequent non-native FP binding to dissociated SE strands. It is commonly proposed that free energy released by the HA2 structural rearrangement catalyzes HA-mediated fusion. This study supports an alternate mechanistic model in which fusion is preceded by FP insertion in the target membrane and formation of the final SE hairpin. Less fusion by the mutants is due to the loss of hairpin stability and consequent reduced level of membrane apposition of the virus and target membranes.

Self-protection of cytosolic malate dehydrogenase against oxidative stress in Arabidopsis

Plant malate dehydrogenase (MDH) isoforms are found in different cell compartments and function in key metabolic pathways. It is well known that the chloroplastic NADP-dependent MDH activities are strictly redox regulated and controlled by light. However, redox dependence of other NAD-dependent MDH isoforms have been less studied. Here, we show by in vitro biochemical characterization that the major cytosolic MDH isoform (cytMDH1) is sensitive to H2O2 through sulfur oxidation of cysteines and methionines. CytMDH1 oxidation affects the kinetics, secondary structure, and thermodynamic stability of cytMDH1. Moreover, MS analyses and comparison of crystal structures between the reduced and H2O2-treated cytMDH1 further show that thioredoxin-reversible homodimerization of cytMDH1 through Cys330 disulfide formation protects the protein from overoxidation. Consistently, we found that cytosolic thioredoxins interact specifically with cytMDH in a yeast two-hybrid system. Importantly, we also show that cytosolic and chloroplastic, but not mitochondrial NAD-MDH activities are sensitive to H2O2 stress in Arabidopsis. NAD-MDH activities decreased both in a catalase2 mutant and in an NADP-thioredoxin reductase mutant, emphasizing the importance of the thioredoxin-reducing system to protect MDH from oxidation in vivo. We propose that the redox switch of the MDH activity contributes to adapt the cell metabolism to environmental constraints.

The forkhead domain hinge-loop plays a pivotal role in DNA binding and transcriptional activity of FOXP2

Forkhead box (FOX) proteins are a ubiquitously expressed family of transcription factors that regulate the development and differentiation of a wide range of tissues in animals. The FOXP subfamily members are the only known FOX proteins capable of forming domain-swapped forkhead domain (FHD) dimers. This is proposed to be due to an evolutionary mutation (P539A) that lies in the FHD hinge loop, a key region thought to fine-tune DNA sequence specificity in the FOX transcription factors. Considering the importance of the hinge loop in both the dimerisation mechanism of the FOXP FHD and its role in tuning DNA binding, a detailed investigation into the implications of mutations within this region could provide important insight into the evolution of the FOX family. Isothermal titration calorimetry and hydrogen exchange mass spectroscopy were used to study the thermodynamic binding signature and changes in backbone dynamics of FOXP2 FHD DNA binding. Dual luciferase reporter assays were performed to study the effect that the hinge-loop mutation has on FOXP2 transcriptional activity in vivo. We demonstrate that the change in dynamics of the hinge-loop region of FOXP2 alters the energetics and mechanism of DNA binding highlighting the critical role of hinge loop mutations in regulating DNA binding characteristics of the FOX proteins.

Localized conformational changes trigger the pH-induced fibrillogenesis of an amyloidogenic λ light chain protein

Solvent conditions modulate the expression of the amyloidogenic potential of proteins. In this work the effect of pH on the fibrillogenic behavior and the conformational properties of 6aJL2, a model protein of the highly amyloidogenic variable light chain λ6a gene segment, was examined. Ordered aggregates showing the ultrastructural and spectroscopic properties observed in amyloid fibrils were formed in the 2.0-8.0 pH range. At pH <3.0 a drastic decrease in lag time and an increase in fibril formation rate were found. In the 4.0-8.0 pH range there was no spectroscopic evidence for significant conformational changes in the native state. Likewise, heat capacity measurements showed no evidence for residual structure in the unfolded state. However, at pH <3.0 stability is severely decreased and the protein suffers conformational changes as detected by circular dichroism, tryptophan and ANS fluorescence, as well as by NMR spectroscopy. Molecular dynamics simulations indicate that acid-induced conformational changes involve the exposure of the loop connecting strands E and F. These results are compatible with pH-induced changes in the NMR spectra. Overall, the results indicate that the mechanism involved in the acid-induced increase in the fibrillogenic potential of 6aJL2 is profoundly different to that observed in κ light chains, and is promoted by localized conformational changes in a region of the protein that was previously not known to be involved in acid-induced light chain fibril formation. The identification of this region opens the potential for the design of specific inhibitors.

Phosphorylation decelerates conformational dynamics in bacterial translation elongation factors

Bacterial protein synthesis is intricately connected to metabolic rate. One of the ways in which bacteria respond to environmental stress is through posttranslational modifications of translation factors. Translation elongation factor Tu (EF-Tu) is methylated and phosphorylated in response to nutrient starvation upon entering stationary phase, and its phosphorylation is a crucial step in the pathway toward sporulation. We analyze how phosphorylation leads to inactivation of Escherichia coli EF-Tu. We provide structural and biophysical evidence that phosphorylation of EF-Tu at T382 acts as an efficient switch that turns off protein synthesis by decoupling nucleotide binding from the EF-Tu conformational cycle. Direct modifications of the EF-Tu switch I region or modifications in other regions stabilizing the β-hairpin state of switch I result in an effective allosteric trap that restricts the normal dynamics of EF-Tu and enables the evasion of the control exerted by nucleotides on G proteins. These results highlight stabilization of a phosphorylation-induced conformational trap as an essential mechanism for phosphoregulation of bacterial translation and metabolism. We propose that this mechanism may lead to the multisite phosphorylation state observed during dormancy and stationary phase.

Zinc Bioavailability from Phytate-Rich Foods and Zinc Supplements. Modeling the Effects of Food Components with Oxygen, Nitrogen, and Sulfur Donor Ligands

Aqueous solubility of zinc phytate (K_sp = (2.6 ± 0.2) × 10^-47 mol⁷/L⁷), essential for zinc bioavailability from plant foods, was found to decrease with increasing temperature corresponding to ΔH_dis of -301 ± 22 kJ/mol and ΔS_dis of -1901 ± 72 J/(mol K). Binding of zinc to phytate was found to be exothermic for the stronger binding site and endothermic for the weaker binding site. The solubility of the slightly soluble zinc citrate and insoluble zinc phytate was found to be considerably enhanced by the food components with oxygen donor, nitrogen donor, and sulfur donor ligands. The driving force for the enhanced solubility is mainly due to the complex formation between zinc and the investigated food components rather than ligand exchange and ternary complex formation as revealed by quantum mechanical calculations and isothermal titration calorimetry. Histidine and citrate are promising ligands for improving zinc absorption from phytate-rich foods.

Comparative analyses of the thermodynamic RNA binding signatures of different types of RNA recognition motifs

RNA recognition motifs (RRMs) are structurally versatile domains important in regulation of alternative splicing. Structural mechanisms of sequence-specific recognition of single-stranded RNAs (ssRNAs) by RRMs are well understood. The thermodynamic strategies are however unclear. Therefore, we utilized microcalorimetry and semi-empirical analyses to comparatively analyze the cognate ssRNA binding thermodynamics of four different RRM domains, each with a different RNA binding mode. The different binding modes are: canonical binding to the β-sheet surface; canonical binding with involvement of N- and C-termini; binding to conserved loops; and binding to an α-helix. Our results identify enthalpy as the sole and general force driving association at physiological temperatures. Also, networks of weak interactions are a general feature regulating stability of the different RRM–ssRNA complexes. In agreement, non-polyelectrolyte effects contributed between ∼75 and 90% of the overall free energy of binding in the considered complexes. The various RNA binding modes also displayed enormous heat capacity differences, that upon dissection revealed large differential changes in hydration, conformations and dynamics upon binding RNA. Altogether, different modes employed by RRMs to bind cognate ssRNAs utilize various thermodynamics strategies during the association process.

Ribosome-dependent Vibrio cholerae mRNAse HigB2 is regulated by a β-strand sliding mechanism

Toxin–antitoxin (TA) modules are small operons involved in bacterial stress response and persistence. higBA operons form a family of TA modules with an inverted gene organization and a toxin belonging to the RelE/ParE superfamily. Here, we present the crystal structures of chromosomally encoded Vibrio cholerae antitoxin (VcHigA2), toxin (VcHigB2) and their complex, which show significant differences in structure and mechanisms of function compared to the higBA module from plasmid Rts1, the defining member of the family. The VcHigB2 is more closely related to Escherichia coli RelE both in terms of overall structure and the organization of its active site. VcHigB2 is neutralized by VcHigA2, a modular protein with an N-terminal intrinsically disordered toxin-neutralizing segment followed by a C-terminal helix-turn-helix dimerization and DNA binding domain. VcHigA2 binds VcHigB2 with picomolar affinity, which is mainly a consequence of entropically favorable de-solvation of a large hydrophobic binding interface and enthalpically favorable folding of the N-terminal domain into an α-helix followed by a β-strand. This interaction displaces helix α3 of VcHigB2 and at the same time induces a one-residue shift in the register of β-strand β3, thereby flipping the catalytically important Arg64 out of the active site.

Onset of disorder and protein aggregation due to oxidation-induced intermolecular disulfide bonds: case study of RRM2 domain from TDP-43

We have investigated the behavior of second RNA-recognition motif (RRM2) of neuropathological protein TDP43 under the effect of oxidative stress as modeled in vitro. Toward this end we have used the specially adapted version of H/D exchange experiment, NMR relaxation and diffusion measurements, dynamic light scattering, controlled proteolysis, gel electrophoresis, site-directed mutagenesis and microsecond MD simulations. Under oxidizing conditions RRM2 forms disulfide-bonded dimers that experience unfolding and then assemble into aggregate particles (APs). These particles are strongly disordered, highly inhomogeneous and susceptible to proteolysis; some of them withstand the dithiothreitol treatment. They can recruit/release monomeric RRM2 through thiol-disulfide exchange reactions. By using a combination of dynamic light scattering and NMR diffusion data we were able to approximate the size distribution function for the APs. The key to the observed aggregation behavior is the diminished ability of disulfide-bonded RRM2 dimers to refold and their increased propensity to misfold, which makes them vulnerable to large thermal fluctuations. The emerging picture provides detailed insight on how oxidative stress can contribute to neurodegenerative disease, with unfolding, aggregation, and proteolytic cleavage as different facets of the process.

Prp40 Homolog A Is a Novel Centrin Target

Pre-mRNA processing protein 40 (Prp40) is a nuclear protein that has a role in pre-mRNA splicing. Prp40 possesses two leucine-rich nuclear export signals, but little is known about the function of Prp40 in the export process. Another protein that has a role in protein export is centrin, a member of the EF-hand superfamily of Ca²⁺-binding proteins. Prp40 was found to be a centrin target by yeast-two-hybrid screening using both Homo sapiens centrin 2 (Hscen2) and Chlamydomonas reinhardtii centrin (Crcen). We identified a centrin-binding site within H. sapiens Prp40 homolog A (HsPrp40A), which contains a hydrophobic triad W₁L₄L₈ that is known to be important in the interaction with centrin. This centrin-binding site is highly conserved within the first nuclear export signal consensus sequence identified in Saccharomyces cerevisiae Prp40. Here, we examine the interaction of HsPrp40A peptide (HsPrp40Ap) with both Hscen2 and Crcen by isothermal titration calorimetry. We employed the thermodynamic parameterization to estimate the polar and apolar surface area of the interface. In addition, we have defined the molecular mechanism of thermally induced unfolding and dissociation of the Crcen-HsPrp40Ap complex using two-dimensional infrared correlation spectroscopy. These complementary techniques showed for the first time, to our knowledge, that HsPrp40Ap interacts with centrin in vitro, supporting a coupled functional role for these proteins in pre-mRNA splicing.

Structural Stability of the Coiled-Coil Domain of Tumor Susceptibility Gene (TSG)-101

The tumor susceptibility gene-101 coiled coil domain (TSG101cc) is an integral component of the endosomal maturation machinery and cytokinesis, and also interacts with several transcription factors. The TSG101cc has been crystallized as a homotetramer but is known to interact with two of its binding partners as a heterotrimer. To investigate this apparent discrepancy, we examined the solution thermodynamics of the TSG101cc. Here, we use circular dichroism, differential scanning calorimetry, analytical ultracentrifugation, fluorescence, and structural thermodynamic analysis to investigate the structural stability and the unfolding of the TSG101cc. We demonstrate that TSG101cc exists in solution primarily as a tetramer, which unfolds in a two-state manner. Surprisingly, no homodimeric or homotrimeric species were detected. Structural thermodynamic analysis of the homotetrameric structure and comparison with known oligomeric coiled-coils suggests that the TSG101cc homotetramer is comparatively unstable on a per residue basis. Furthermore, the homotrimeric coiled-coil is predicted to be much less stable than the functional heterotrimeric coiled-coil in the endosomal sorting complex required for transport 1 (ESCRT1). These results support a model whereby the tetramer–monomer equilibrium of TSG101 serves as the cellular reservoir of TSG101, which is effectively outcompeted when its binding partners are present and the heteroternary complex can form.

Temperature stability of proteins: Analysis of irreversible denaturation using isothermal calorimetry

The structural stability of proteins has been traditionally studied under conditions in which the folding/unfolding reaction is reversible, since thermodynamic parameters can only be determined under these conditions. Achieving reversibility conditions in temperature stability experiments has often required performing the experiments at acidic pH or other nonphysiological solvent conditions. With the rapid development of protein drugs, the fastest growing segment in the pharmaceutical industry, the need to evaluate protein stability under formulation conditions has acquired renewed urgency. Under formulation conditions and the required high protein concentration (∼100 mg/mL), protein denaturation is irreversible and frequently coupled to aggregation and precipitation. In this article, we examine the thermal denaturation of hen egg white lysozyme (HEWL) under irreversible conditions and concentrations up to 100 mg/mL using several techniques, especially isothermal calorimetry which has been used to measure the enthalpy and kinetics of the unfolding and aggregation/precipitation at 12°C below the transition temperature measured by DSC. At those temperatures the rate of irreversible protein denaturation and aggregation of HEWL is measured to be on the order of 1 day^-1 . Isothermal calorimetry appears a suitable technique to identify buffer formulation conditions that maximize the long term stability of protein drugs.

Binding Mechanism of the N-Terminal SH3 Domain of CrkII and Proline-Rich Motifs in cAbl

The N-terminal Src homology 3 (nSH3) domain of a signaling adaptor protein, CT-10 regulator of kinase II (CrkII), recognizes proline-rich motifs (PRMs) of binding partners, such as cAbl kinase. The interaction between CrkII and cAbl kinase is involved in the regulation of cell spreading, microbial pathogenesis, and cancer metastasis. Here, we report the detailed biophysical characterizations of the interactions between the nSH3 domain of CrkII and PRMs in cAbl. We identified that the nSH3 domain of CrkII binds to three PRMs in cAbl with virtually identical affinities. Structural studies, by using x-ray crystallography and NMR spectroscopy, revealed that the binding modes of all three nSH3:PRM complexes are highly similar to each other. Van 't Hoff analysis revealed that nSH3:PRM interaction is associated with favorable enthalpy and unfavorable entropy change. The combination of experimentally determined thermodynamic parameters, structure-based calculations, and (15)N NMR relaxation analysis highlights the energetic contribution of conformational entropy change upon the complex formation, and water molecules structured in the binding interface of the nSH3:PRM complex. Understanding the molecular basis of nSH3:PRM interaction will provide, to our knowledge, new insights for the rational design of small molecules targeting the interaction between CrkII and cAbl.

Multiple helical conformations of the helix-turn-helix region revealed by NOE-restrained MD simulations of tryptophan aporepressor, TrpR

The nature of flexibility in the helix-turn-helix region of E. coli trp aporepressor has been unexplained for many years. The original ensemble of nuclear magnetic resonance (NMR structures showed apparent disorder, but chemical shift and relaxation measurements indicated a helical region. Nuclear Overhauser effect (NOE) data for a temperature-sensitive mutant showed more helical character in its helix-turn-helix region, but nevertheless also led to an apparently disordered ensemble. However, conventional NMR structure determination methods require all structures in the ensemble to be consistent with every NOE simultaneously. This work uses an alternative approach in which some structures of the ensemble are allowed to violate some NOEs to permit modeling of multiple conformational states that are in dynamic equilibrium. Newly measured NOE data for wild-type aporepressor are used as time-averaged distance restraints in molecular dynamics simulations to generate an ensemble of helical conformations that is more consistent with the observed NMR data than the apparent disorder in the previously reported NMR structures. The results indicate the presence of alternating helical conformations that provide a better explanation for the flexibility of the helix-turn-helix region of trp aporepressor. Structures representing these conformations have been deposited with PDB ID: 5TM0. Proteins 2017; 85:731-740. © 2016 Wiley Periodicals, Inc.

Intrinsic α helix propensities compact hydrodynamic radii in intrinsically disordered proteins

Proteins that lack tertiary stability under normal conditions, known as intrinsically disordered, exhibit a wide range of biological activities. Molecular descriptions for the biology of intrinsically disordered proteins (IDPs) consequently rely on disordered structural models, which in turn require experiments that assess the origins to structural features observed. For example, while hydrodynamic size is mostly insensitive to sequence composition in chemically denatured proteins, IDPs show strong sequence-specific effects in the hydrodynamic radius (R_h ) when measured under normal conditions. To investigate sequence-modulation of IDP R_h , disordered ensembles generated by a hard sphere collision model modified with a structure-based parameterization of the solution energetics were used to parse the contributions of net charge, main chain dihedral angle bias, and excluded volume on hydrodynamic size. Ensembles for polypeptides 10-35 residues in length were then used to establish power-law scaling relationships for comparison to experimental R_h from 26 IDPs. Results showed the expected outcomes of increased hydrodynamic size from increases in excluded volume and net charge, and compaction from chain-solvent interactions. Chain bias representing intrinsic preferences for α helix and polyproline II (PP_II ), however, modulated R_h with intricate dependence on the simulated propensities. PP_II propensities at levels expected in IDPs correlated with heightened R_h sensitivity to even weak α helix propensities, indicating bias for common (φ, ψ) are important determinants of hydrodynamic size. Moreover, data show that IDP R_h can be predicted from sequence with good accuracy from a small set of physicochemical properties, namely intrinsic conformational propensities and net charge. Proteins 2017; 85:296-311. © 2016 Wiley Periodicals, Inc.