Determination of Conformational Entropy of Fully and Partially Folded Conformations of Holo- and Apomyoglobin

Protein Stability─Analysis of Heat and Cold Denaturation without and with Unfolding Models

Protein stability is important in many areas of life sciences. Thermal protein unfolding is investigated extensively with various spectroscopic techniques. The extraction of thermodynamic properties from these measurements requires the application of models. Differential scanning calorimetry (DSC) is less common, but is unique as it measures directly a thermodynamic property, that is, the heat capacity C_p(T). The analysis of C_p(T) is usually performed with the chemical equilibrium two-state model. This is not necessary and leads to incorrect thermodynamic consequences. Here we demonstrate a straightforward model-independent evaluation of heat capacity experiments in terms of protein unfolding enthalpy ΔH(T), entropy ΔS(T), and free energy ΔG(T)). This now allows the comparison of the experimental thermodynamic data with the predictions of different models. We critically examined the standard chemical equilibrium two-state model, which predicts a positive free energy for the native protein, and diverges distinctly from the experimental temperature profiles. We propose two new models which are equally applicable to spectroscopy and calorimetry. The Θ_U(T)-weighted chemical equilibrium model and the statistical-mechanical two-state model provide excellent fits of the experimental data. They predict sigmoidal temperature profiles for enthalpy and entropy, and a trapezoidal temperature profile for the free energy. This is illustrated with experimental examples for heat and cold denaturation of lysozyme and β-lactoglobulin. We then show that the free energy is not a good criterion to judge protein stability. More useful parameters are discussed, including protein cooperativity. The new parameters are embedded in a well-defined thermodynamic context and are amenable to molecular dynamics calculations.

Complementary Experimental Methods to Obtain Thermodynamic Parameters of Protein Ligand Systems

In recent years, thermophoresis has emerged as a promising tool for quantifying biomolecular interactions. The underlying microscopic physical effect is still not understood, but often attributed to changes in the hydration layer once the binding occurs. To gain deeper insight, we investigate whether non-equilibrium coefficients can be related to equilibrium properties. Therefore, we compare thermophoretic data measured by thermal diffusion forced Rayleigh scattering (TDFRS) (which is a non-equilibrium process) with thermodynamic data obtained by isothermal titration calorimetry (ITC) (which is an equilibrium process). As a reference system, we studied the chelation reaction between ethylenediaminetetraacetic acid (EDTA) and calcium chloride (CaCl2) to relate the thermophoretic behavior quantified by the Soret coefficient ST to the Gibb's free energy ΔG determined in the ITC experiment using an expression proposed by Eastman. Finally, we have studied the binding of the protein Bovine Carbonic Anhydrase I (BCA I) to two different benzenesulfonamide derivatives: 4-fluorobenzenesulfonamide (4FBS) and pentafluorobenzenesulfonamide (PFBS). For all three systems, we find that the Gibb's free energies calculated from ST agree with ΔG from the ITC experiment. In addition, we also investigate the influence of fluorescent labeling, which allows measurements in a thermophoretic microfluidic cell. Re-examination of the fluorescently labeled system using ITC showed a strong influence of the dye on the binding behavior.

Probing terahertz dynamics of multidomain protein in cell-like confinement

The development of meaningful descriptions of multidomain proteins exhibiting complex inter-domain dynamics modes is a key challenge for understanding their roles in molecular recognition and signalling processes. Here we developed a generally applicable approach for probing the low frequency collective hydration dynamics of multidomain proteins that uses terahertz spectroscopy of a protein molecule confined in a phospholipid reverse micelles environment (named Droplet THz). With the combination of normal mode analysis, we demonstrated the binding of calcium ions modulates the local inter-domain motion of the human coagulant factor VIII protein in a concentration-dependent manner. These findings highlight the Droplet THz as a valuable tool for dissecting the ultrafast dynamics of domain motion in the multidomain proteins and suggest a modulating mechanism of calcium ions on the structural flexibility and function of human coagulant proteins.

Variation of Structural and Dynamical Flexibility of Myelin Basic Protein in Response to Guanidinium Chloride

Myelin basic protein (MBP) is intrinsically disordered in solution and is considered as a conformationally flexible biomacromolecule. Here, we present a study on perturbation of MBP structure and dynamics by the denaturant guanidinium chloride (GndCl) using small-angle scattering and neutron spin–echo spectroscopy (NSE). A concentration of 0.2 M GndCl causes charge screening in MBP resulting in a compact, but still disordered protein conformation, while GndCl concentrations above 1 M lead to structural expansion and swelling of MBP. NSE data of MBP were analyzed using the Zimm model with internal friction (ZIF) and normal mode (NM) analysis. A significant contribution of internal friction was found in compact states of MBP that approaches a non-vanishing internal friction relaxation time of approximately 40 ns at high GndCl concentrations. NM analysis demonstrates that the relaxation rates of internal modes of MBP remain unaffected by GndCl, while structural expansion due to GndCl results in increased amplitudes of internal motions. Within the model of the Brownian oscillator our observations can be rationalized by a loss of friction within the protein due to structural expansion. Our study highlights the intimate coupling of structural and dynamical plasticity of MBP, and its fundamental difference to the behavior of ideal polymers in solution.

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.

Self-Diffusive Properties of the Intrinsically Disordered Protein Histatin 5 and the Impact of Crowding Thereon: A Combined Neutron Spectroscopy and Molecular Dynamics Simulation Study

Intrinsically disordered proteins (IDPs) are proteins that, in comparison with globular/structured proteins, lack a distinct tertiary structure. Here, we use the model IDP, Histatin 5, for studying its dynamical properties under self-crowding conditions with quasi-elastic neutron scattering in combination with full atomistic molecular dynamics (MD) simulations. The aim is to determine the effects of crowding on the center-of-mass diffusion as well as the internal diffusive behavior. The diffusion was found to decrease significantly, which we hypothesize can be attributed to some degree of aggregation at higher protein concentrations, (≥100 mg/mL), as indicated by recent small-angle X-ray scattering studies. Temperature effects are also considered and found to, largely, follow Stokes–Einstein behavior. Simple geometric considerations fail to accurately predict the rates of diffusion, while simulations show semiquantitative agreement with experiments, dependent on assumptions of the ratio between translational and rotational diffusion. A scaling law that previously was found to successfully describe the behavior of globular proteins was found to be inadequate for the IDP, Histatin 5. Analysis of the MD simulations show that the width of the distribution with respect to diffusion is not a simplistic mirroring of the distribution of radius of gyration, hence, displaying the particular features of IDPs that need to be accounted for.

Molecular Dynamics of Lysozyme Amyloid Polymorphs Studied by Incoherent Neutron Scattering

Lysozyme amyloidosis is a hereditary disease, which is characterized by the deposition of lysozyme amyloid fibrils in various internal organs. It is known that lysozyme fibrils show polymorphism and that polymorphs formed at near-neutral pH have the ability to promote more monomer binding than those formed at acidic pH, indicating that only specific polymorphs become dominant species in a given environment. This is likely due to the polymorph-specific configurational diffusion. Understanding the possible differences in dynamical behavior between the polymorphs is thus crucial to deepen our knowledge of amyloid polymorphism and eventually elucidate the molecular mechanism of lysozyme amyloidosis. In this study, molecular dynamics at sub-nanosecond timescale of two kinds of polymorphic fibrils of hen egg white lysozyme, which has long been used as a model of human lysozyme, formed at pH 2.7 (LP27) and pH 6.0 (LP60) was investigated using elastic incoherent neutron scattering (EINS) and quasi-elastic neutron scattering (QENS). Analysis of the EINS data showed that whereas the mean square displacement of atomic motions is similar for both LP27 and LP60, LP60 contains a larger fraction of atoms moving with larger amplitudes than LP27, indicating that the dynamical difference between the two polymorphs lies not in the averaged amplitude, but in the distribution of the amplitudes. Furthermore, analysis of the QENS data showed that the jump diffusion coefficient of atoms is larger for LP60, suggesting that the atoms of LP60 undergo faster diffusive motions than those of LP27. This study thus characterizes the dynamics of the two lysozyme polymorphs and reveals that the molecular dynamics of LP60 is enhanced compared with that of LP27. The higher molecular flexibility of the polymorph would permit to adjust its conformation more quickly than its counterpart, facilitating monomer binding.

Complementary approaches to obtaining thermodynamic parameters from protein ligand systems-challenges and opportunities

Protein ligand interactions play an important role in biology. Increasingly the aim is to understand and influence protein ligand binding. The binding process is heavily influenced by its thermodynamic parameters. In order to understand how the whole system thermodynamics work it is important to characterise the individual contribution of each of the systems components. While the change in conformational entropy of the protein can be determined using QENS complementary methods are necessary in order to characterise all components. This paper will describe the challenges that can occur when combining the different methods, as well as how they can be overcome.

Light Scattering and Absorption Complementarities to Neutron Scattering: In Situ FTIR and DLS Techniques at the High-Intensity and Extended Q-Range SANS Diffractometer KWS-2

Understanding soft and biological materials requires global knowledge of their microstructural features from elementary units at the nm scale up to larger complex aggregates in the micrometer range. Such a wide range of scale can be explored using the KWS-2 small-angle neutron (SANS) diffractometer. Additional information obtained by in situ complementary techniques sometimes supports the SANS analysis of systems undergoing structural modifications under external stimuli or which are stable only for short times. Observations at the local molecular level structure and conformation assists with an unambiguous interpretation of the SANS data using appropriate structural models, while monitoring of the sample condition during the SANS investigation ensures the sample stability and desired composition and chemical conditions. Thus, we equipped the KWS-2 with complementary light absorption and scattering capabilities: Fourier transform infrared (FTIR) spectroscopy can now be performed simultaneously with standard and time-resolved SANS, while in situ dynamic light scattering (DLS) became available for routine experiments, which enables the observation of either changes in the sample composition, due to sedimentation effects, or in size of morphologies, due to aggregation processes. The performance of each setup is demonstrated here using systems representative of those typically investigated on this beamline and benchmarked to studies performed offline.

Mobile Protons Limit the Stability of Salt Bridges in the Gas Phase: Implications for the Structures of Electrosprayed Protein Ions

Electrosprayed protein ions can retain native-like conformations. The intramolecular contacts that stabilize these compact gas-phase structures remain poorly understood. Recent work has uncovered abundant salt bridges in electrosprayed proteins. Salt bridges are zwitterionic BH⁺/A^- contacts. The low dielectric constant in the vacuum strengthens electrostatic interactions, suggesting that salt bridges could be a key contributor to the retention of compact protein structures. A problem with this assertion is that H⁺ are mobile, such that H⁺ transfer can convert salt bridges into neutral B⁰/HA⁰ contacts. This possible salt bridge annihilation puts into question the role of zwitterionic motifs in the gas phase, and it calls for a detailed analysis of BH⁺/A^- versus B⁰/HA⁰ interactions. Here, we investigate this issue using molecular dynamics (MD) simulations and electrospray experiments. MD data for short model peptides revealed that salt bridges with static H⁺ have dissociation energies around 700 kJ mol^-1. The corresponding B⁰/HA⁰ contacts are 1 order of magnitude weaker. When considering the effects of mobile H⁺, BH⁺/A^- bond energies were found to be between these two extremes, confirming that H⁺ migration can significantly weaken salt bridges. Next, we examined the protein ubiquitin under collision-induced unfolding (CIU) conditions. CIU simulations were conducted using three different MD models: (i) Positive-only runs with static H⁺ did not allow for salt bridge formation and produced highly expanded CIU structures. (ii) Zwitterionic runs with static H⁺ resulted in abundant salt bridges, culminating in much more compact CIU structures. (iii) Mobile H⁺ simulations allowed for the dynamic formation/annihilation of salt bridges, generating CIU structures intermediate between scenarios (i) and (ii). Our results uncover that mobile H⁺ limit the stabilizing effects of salt bridges in the gas phase. Failure to consider the effects of mobile H⁺ in MD simulations will result in unrealistic outcomes under CIU conditions.

Effect of red blood cell shape changes on haemoglobin interactions and dynamics: a neutron scattering study

By using a combination of experimental neutron scattering techniques, it is possible to obtain a statistical perspective on red blood cell (RBC) shape in suspensions, and the inter-relationship with protein interactions and dynamics inside the confinement of the cell membrane. In this study, we examined the ultrastructure of RBC and protein-protein interactions of haemoglobin (Hb) in them using ultra-small-angle neutron scattering and small-angle neutron scattering (SANS). In addition, we used the neutron backscattering method to access Hb motion on the ns time scale and Å length scale. Quasi-elastic neutron scattering (QENS) experiments were performed to measure diffusive motion of Hb in RBCs and in an RBC lysate. By using QENS, we probed both internal Hb dynamics and global protein diffusion, on the accessible time scale and length scale by QENS. Shape changes of RBCs and variation of intracellular Hb concentration were induced by addition of the Na⁺-selective ionophore monensin and the K⁺-selective one, valinomycin. The experimental SANS and QENS results are discussed within the framework of crowded protein solutions, where free motion of Hb is obstructed by mutual interactions.

Thermophoresis: The Case of Streptavidin and Biotin

Thermophoretic behavior of a free protein changes upon ligand binding and gives access to information on the binding constants. The Soret effect has also been proven to be a promising tool to gain information on the hydration layer, as the temperature dependence of the thermodiffusion behavior is sensitive to solute–solvent interactions. In this work, we perform systematic thermophoretic measurements of the protein streptavidin (STV) and of the complex STV with biotin (B) using thermal diffusion forced Rayleigh scattering (TDFRS). Our experiments show that the temperature sensitivity of the Soret coefficient is reduced for the complex compared to the free protein. We discuss our data in comparison with recent quasi-elastic neutron scattering (QENS) measurements. As the QENS measurement has been performed in heavy water, we perform additional measurements in water/heavy water mixtures. Finally, we also elucidate the challenges arising from the quantiative thermophoretic study of complex multicomponent systems such as protein solutions.

Transition between protein-like and polymer-like dynamic behavior: Internal friction in unfolded apomyoglobin depends on denaturing conditions

Equilibrium dynamics of different folding intermediates and denatured states is strongly connected to the exploration of the conformational space on the nanosecond time scale and might have implications in understanding protein folding. For the first time, the same protein system apomyoglobin has been investigated using neutron spin-echo spectroscopy in different states: native-like, partially folded (molten globule) and completely unfolded, following two different unfolding paths: using acid or guanidinium chloride (GdmCl). While the internal dynamics of the native-like state can be understood using normal mode analysis based on high resolution structural information of myoglobin, for the unfolded and even for the molten globule states, models from polymer science are employed. The Zimm model accurately describes the slowly-relaxing, expanded GdmCl-denaturated state, ignoring the individuality of the different aminoacid side chain. The dynamics of the acid unfolded and molten globule state are similar in the framework of the Zimm model with internal friction, where the chains still interact and hinder each other: the first Zimm relaxation time is as large as the internal friction time. Transient formation of secondary structure elements in the acid unfolded and presence of α-helices in the molten globule state lead to internal friction to a similar extent.

Strong Adverse Contribution of Conformational Dynamics to Streptavidin-Biotin Binding

Molecular dynamics plays an important role for the biological function of proteins. For protein ligand interactions, changes of conformational entropy of protein and hydration layer are relevant for the binding process. Quasielastic neutron scattering (QENS) was used to investigate differences in protein dynamics and conformational entropy of ligand-bound and ligand-free streptavidin. Protein dynamics were probed both on the fast picosecond time scale using neutron time-of-flight spectroscopy and on the slower nanosecond time scale using high-resolution neutron backscattering spectroscopy. We found the internal equilibrium motions of streptavidin and the corresponding mean square displacements (MSDs) to be greatly reduced upon biotin binding. On the basis of the observed MSDs, we calculated the difference of conformational entropy ΔS_conf of the protein component between ligand-bound and ligand-free streptavidin. The rather large negative ΔS_conf value (-2 kJ mol^-1 K^-1 on the nanosecond time scale) obtained for the streptavidin tetramer seems to be counterintuitive, given the exceptionally high affinity of streptavidin-biotin binding. Literature data on the total entropy change ΔS observed upon biotin binding to streptavidin, which includes contributions from both the protein and the hydration water, suggest partial compensation of the unfavorable ΔS_conf by a large positive entropy gain of the surrounding hydration layer and water molecules that are displaced during ligand binding.

The first study on the impact of osmolytes in whole cells of high temperature-adapted microorganisms

The hyperthermophilic piezophile, Thermococcus barophilus displays a strong stress response characterized by the accumulation of the organic osmolyte, mannosylglycerate during growth under sub-optimal pressure conditions (0.1 MPa). Taking advantage of this known effect, the impact of osmolytes in piezophiles in an otherwise identical cellular context was investigated, by comparing T. barophilus cells grown under low or optimal pressures (40 MPa). Using neutron scattering techniques, we studied the molecular dynamics of live cells of T. barophilus at different pressures and temperatures. We show that in the presence of osmolytes, cells present a higher diffusion coefficient of hydration water and an increase of bulk water motions at a high temperature. In the absence of osmolytes, the T. barophilus cellular dynamics is more responsive to high temperature and high hydrostatic pressure. These results therefore give clear evidence for a protecting effect of osmolytes on proteins.

Conformational selection of the intrinsically disordered plant stress protein COR15A in response to solution osmolarity - an X-ray and light scattering study

The plant stress protein COR15A stabilizes chloroplast membranes during freezing. COR15A is an intrinsically disordered protein (IDP) in aqueous solution, but acquires an α-helical structure during dehydration or the increase of solution osmolarity. We have used small- and wide-angle X-ray scattering (SAXS/WAXS) combined with static and dynamic light scattering (SLS/DLS) to investigate the structural and hydrodynamic properties of COR15A in response to increasing solution osmolarity. Coarse-grained ensemble modelling allowed a structure-based interpretation of the SAXS data. Our results demonstrate that COR15A behaves as a biomacromolecule with polymer-like properties which strongly depend on solution osmolarity. Biomacromolecular self-assembly occurring at high solvent osmolarity is initiated by the occurrence of two specific structural subpopulations of the COR15A monomer. The osmolarity dependent structural selection mechanism is an elegant way for conformational regulation and assembly of COR15A. It highlights the importance of the polymer-like properties of IDPs for their associated biological function.

Role of Hydrogen Bonding of Cyclodextrin-Drug Complexes Probed by Thermodiffusion

Temperature gradient-induced migration of biomolecules, known as thermophoresis or thermodiffusion, changes upon ligand binding. In recent years, this effect has been used to determine protein-ligand binding constants. The mechanism through which thermodiffusive properties change when complexes are formed, however, is not understood. An important contribution to thermodiffusive properties originates from the thermal response of hydrogen bonds. Because there is a considerable difference between the degree of solvation of the protein-ligand complex and its isolated components, ligand-binding is accompanied by a significant change in hydration. The aim of the present work is therefore to investigate the role played by hydrogen bonding on the change in thermodiffusive behavior upon ligand-binding. As a model system, we use cyclodextrins (CDs) and acetylsalicylic acid (ASA), where quite a significant change in hydration is expected and where no conformational changes occur when a CD/ASA complex is formed in aqueous solution. Thermophoresis was investigated in the temperature range of 10-50 °C by infrared thermal diffusion forced Rayleigh scattering. Nuclear magnetic resonance measurements were performed at 25 °C to obtain information about the structure of the complexes. All CD/ASA complexes show a stronger affinity toward regions of lower temperature compared to the free CDs. We found that the temperature sensitivity of thermophoresis correlates with the 1-octanol/water partition coefficient. This observation not only establishes the relation between thermodiffusion and degree of hydrogen bonding but also opens the possibility to relate thermodiffusive properties of complexes to their partition coefficient, which cannot be determined otherwise. This concept is especially interesting for protein-ligand complexes where the protein undergoes a conformational change, different from the CD/ASA complexes, giving rise to additional changes in their hydrophilicity.

Photoactivation Reduces Side-Chain Dynamics of a LOV Photoreceptor

We used neutron-scattering experiments to probe the conformational dynamics of the light, oxygen, voltage (LOV) photoreceptor PpSB1-LOV from Pseudomonas putida in both the dark and light states. Global protein diffusion and internal macromolecular dynamics were measured using incoherent neutron time-of-flight and backscattering spectroscopy on the picosecond to nanosecond timescales. Global protein diffusion of PpSB1-LOV is not influenced by photoactivation. Observation-time-dependent global diffusion coefficients were found, which converge on the nanosecond timescale toward diffusion coefficients determined by dynamic light scattering. Mean-square displacements of localized internal motions and effective force constants, <k'>, describing the resilience of the proteins were determined on the respective timescales. Photoactivation significantly modifies the flexibility and the resilience of PpSB1-LOV. On the fast, picosecond timescale, small changes in the mean-square displacement and <k'> are observed, which are enhanced on the slower, nanosecond timescale. Photoactivation results in a slightly larger resilience of the photoreceptor on the fast, picosecond timescale, whereas in the nanosecond range, a significantly less resilient structure of the light-state protein is observed. For a residue-resolved interpretation of the experimental neutron-scattering data, we analyzed molecular dynamics simulations of the PpSB1-LOV X-ray structure. Based on these data, it is tempting to speculate that light-induced changes in the protein result in altered side-chain mobility mostly for residues on the protruding Jα helix and on the LOV-LOV dimer interface. Our results provide strong experimental evidence that side-chain dynamics play a crucial role in photoactivation and signaling of PpSB1-LOV via modulation of conformational entropy.

Dynamical Behavior of Human α-Synuclein Studied by Quasielastic Neutron Scattering

α-synuclein (αSyn) is a protein consisting of 140 amino acid residues and is abundant in the presynaptic nerve terminals in the brain. Although its precise function is unknown, the filamentous aggregates (amyloid fibrils) of αSyn have been shown to be involved in the pathogenesis of Parkinson's disease, which is a progressive neurodegenerative disorder. To understand the pathogenesis mechanism of this disease, the mechanism of the amyloid fibril formation of αSyn must be elucidated. Purified αSyn from bacterial expression is monomeric but intrinsically disordered in solution and forms amyloid fibrils under various conditions. As a first step toward elucidating the mechanism of the fibril formation of αSyn, we investigated dynamical behavior of the purified αSyn in the monomeric state and the fibril state using quasielastic neutron scattering (QENS). We prepared the solution sample of 9.5 mg/ml purified αSyn, and that of 46 mg/ml αSyn in the fibril state, both at pD 7.4 in D2O. The QENS experiments on these samples were performed using the near-backscattering spectrometer, BL02 (DNA), at the Materials and Life Science Facility at the Japan Accelerator Research Complex, Japan. Analysis of the QENS spectra obtained shows that diffusive global motions are observed in the monomeric state but largely suppressed in the fibril state. However, the amplitude of the side chain motion is shown to be larger in the fibril state than in the monomeric state. This implies that significant solvent space exists within the fibrils, which is attributed to the αSyn molecules within the fibrils having a distribution of conformations. The larger amplitude of the side chain motion in the fibril state than in the monomeric state implies that the fibril state is entropically favorable.

Picosecond to nanosecond dynamics provide a source of conformational entropy for protein folding

Myoglobin can be trapped in fully folded structures, partially folded molten globules, and unfolded states under stable equilibrium conditions.