Probing the physical determinants of thermal expansion of folded proteins

Thermal Lens Measurements of Thermal Expansivity in Thermosensitive Polymer Solutions

The weak absorption of a laser beam generates in a fluid an inhomogeneous refractive index profile acting as a negative lens. This self-effect on beam propagation, known as Thermal Lensing (TL), is extensively exploited in sensitive spectroscopic techniques, and in several all-optical methods for the assessment of thermo-optical properties of simple and complex fluids. Using the Lorentz-Lorenz equation, we show that the TL signal is directly proportional to the sample thermal expansivity α, a feature allowing minute density changes to be detected with high sensitivity in a tiny sample volume, using a simple optical scheme. We took advantage of this key result to investigate the compaction of PniPAM microgels occurring around their volume phase transition temperature, and the temperature-driven formation of poloxamer micelles. For both these different kinds of structural transitions, we observed a significant peak in the solute contribution to α, indicating a decrease in the overall solution density-rather counterintuitive evidence that can nevertheless be attributed to the dehydration of the polymer chains. Finally, we compare the novel method we propose with other techniques currently used to obtain specific volume changes.

Protein-Ligand Binding Volume Determined from a Single 2D NMR Spectrum with Increasing Pressure

Proteins undergo changes in their partial volumes in numerous biological processes such as enzymatic catalysis, unfolding–refolding, and ligand binding. The change in the protein volume upon ligand binding—a parameter termed the protein–ligand binding volume—can be extensively studied by high-pressure NMR spectroscopy. In this study, we developed a method to determine the protein–ligand binding volume from a single two-dimensional (2D) ¹H–¹⁵N heteronuclear single quantum coherence (HSQC) spectrum at different pressures, if the exchange between ligand-free and ligand-bound states of a protein is slow in the NMR time-scale. This approach required a significantly lower amount of protein and NMR time to determine the protein–ligand binding volume of two carbonic anhydrase isozymes upon binding their ligands. The proposed method can be used in other protein–ligand systems and expand the knowledge about protein volume changes upon small-molecule binding.

Single molecule secondary structure determination of proteins through infrared absorption nanospectroscopy

The chemical and structural properties of biomolecules determine their interactions, and thus their functions, in a wide variety of biochemical processes. Innovative imaging methods have been developed to characterise biomolecular structures down to the angstrom level. However, acquiring vibrational absorption spectra at the single molecule level, a benchmark for bulk sample characterization, has remained elusive. Here, we introduce off-resonance, low power and short pulse infrared nanospectroscopy (ORS-nanoIR) to allow the acquisition of infrared absorption spectra and chemical maps at the single molecule level, at high throughput on a second timescale and with a high signal-to-noise ratio (~10-20). This high sensitivity enables the accurate determination of the secondary structure of single protein molecules with over a million-fold lower mass than conventional bulk vibrational spectroscopy. These results pave the way to probe directly the chemical and structural properties of individual biomolecules, as well as their interactions, in a broad range of chemical and biological systems.

q-Canonical Monte Carlo Sampling for Modeling the Linkage between Charge Regulation and Conformational Equilibria of Peptides

The overall charge content and the patterning of charged residues have a profound impact on the conformational ensembles adopted by intrinsically disordered proteins. These parameters can be altered by charge regulation, which refers to the effects of post-translational modifications, pH-dependent changes to charge, and conformational fluctuations that modify the pKa values of ionizable residues. Although atomistic simulations have played a prominent role in uncovering the major sequence-ensemble relationships of IDPs, most simulations assume fixed charge states for ionizable residues. This may lead to erroneous estimates for conformational equilibria if they are linked to charge regulation. Here, we report the development of a new method we term q-canonical Monte Carlo sampling for modeling the linkage between charge regulation and conformational equilibria. The method, which is designed to be interoperable with the ABSINTH implicit solvation model, operates as follows: For a protein sequence with n ionizable residues, we start with all 2n charge microstates and use a criterion based on model compound pKa values to prune down to a subset of thermodynamically relevant charge microstates. This subset is then grouped into mesostates, where all microstates that belong to a mesostate have the same net charge. Conformational distributions, drawn from a canonical ensemble, are generated for each of the charge microstates that make up a mesostate using a method we designate as proton walk sampling. This method combines Metropolis Monte Carlo sampling in conformational space with an auxiliary Markov process that enables interconversions between charge microstates along a mesostate. Proton walk sampling helps identify the most likely charge microstate per mesostate. We then use thermodynamic integration aided by the multistate Bennett acceptance ratio method to estimate the free energies for converting between mesostates. These free energies are then combined with the per-microstate weights along each mesostate to estimate standard state free energies and pH-dependent free energies for all thermodynamically relevant charge microstates. The results provide quantitative estimates of the probabilities and preferred conformations associated with every thermodynamically accessible charge microstate. We showcase the application of q-canonical sampling using two model systems. The results establish the soundness of the method and the importance of charge regulation in systems characterized by conformational heterogeneity.

Lessons from pressure denaturation of proteins

Although it is now relatively well understood how sequence defines and impacts global protein stability in specific structural contexts, the question of how sequence modulates the configurational landscape of proteins remains to be defined. Protein configurational equilibria are generally characterized by using various chemical denaturants or by changing temperature or pH. Another thermodynamic parameter which is less often used in such studies is high hydrostatic pressure. This review discusses the basis for pressure effects on protein structure and stability, and describes how the unique mechanisms of pressure-induced unfolding can provide unique insights into protein conformational landscapes.

Practical aspects of high-pressure NMR spectroscopy and its applications in protein biophysics and structural biology

Pressure and temperature are the two fundamental variables of thermodynamics. Temperature and chemical perturbation are central experimental tools for the exploration of macromolecular structure and dynamics. Though it has long been recognized that hydrostatic pressure offers a complementary and often unique view of macromolecular structure, stability and dynamics, it has not been employed nearly as much. For solution NMR applications the limited use of high-pressure is undoubtedly traced to difficulties of employing pressure in the context of modern multinuclear and multidimensional NMR. Limitations in pressure tolerant NMR sample cells have been overcome and enable detailed studies of macromolecular energy landscapes, hydration, dynamics and function. Here we review the practical considerations for studies of biological macromolecules at elevated pressure, with a particular emphasis on applications in protein biophysics and structural biology.

Molecular Determinants of Temperature Dependence of Protein Volume Change upon Unfolding

Pressure is a well-known environmental stressor that can either stabilize or destabilize proteins. The volumetric change upon protein unfolding determines the effect of pressure on protein stability, where negative volume changes destabilized proteins at high pressures. High temperature often accompanies high pressure, for example, in the ocean depths near hydrothermal vents or near faults, so it is important to study the effect of temperature on the volumetric change upon unfolding. We previously detailed the magnitude and sign of the molecular determinants of volumetric change, allowing us to quantitatively predict the volumetric change upon protein unfolding. Here, we present a comprehensive analysis of the temperature dependence of the volumetric components of proteins, showing that hydration volume is the primary component that defines expansivities of the native and unfolded states and void volume only contributes slightly to the folded state expansivity.

A complete thermodynamic analysis of enzyme turnover links the free energy landscape to enzyme catalysis

Our understanding of how enzymes work is coloured by static structure depictions where the enzyme scaffold is presented as either immobile, or in equilibrium between well-defined static conformations. Proteins, however, exhibit a large degree of motion over a broad range of timescales and magnitudes and this is defined thermodynamically by the enzyme free energy landscape (FEL). The role and importance of enzyme motion is extremely contentious. Much of the challenge is in the experimental detection of so called 'conformational sampling' involved in enzyme turnover. Herein we apply combined pressure and temperature kinetics studies to elucidate the full suite of thermodynamic parameters defining an enzyme FEL as it relates to enzyme turnover. We find that the key thermodynamic parameters governing vibrational modes related to enzyme turnover are the isobaric expansivity term and the change in heat capacity for enzyme catalysis. Variation in the enzyme FEL affects these terms. Our analysis is supported by a range of biophysical and computational approaches that specifically capture information on protein vibrational modes and the FEL (all atom flexibility calculations, red edge excitation shift spectroscopy and viscosity studies) that provide independent evidence for our findings. Our data suggest that restricting the enzyme FEL may be a powerful strategy when attempting to rationally engineer enzymes, particularly to alter thermal activity. Moreover, we demonstrate how rational predictions can be made with a rapid computational approach.

High-Resolution Mapping of a Repeat Protein Folding Free Energy Landscape

A complete description of the pathways and mechanisms of protein folding requires a detailed structural and energetic characterization of the conformational ensemble along the entire folding reaction coordinate. Simulations can provide this level of insight for small proteins. In contrast, with the exception of hydrogen exchange, which does not monitor folding directly, experimental studies of protein folding have not yielded such structural and energetic detail. NMR can provide residue specific atomic level structural information, but its implementation in protein folding studies using chemical or temperature perturbation is problematic. Here we present a highly detailed structural and energetic map of the entire folding landscape of the leucine-rich repeat protein, pp32 (Anp32), obtained by combining pressure-dependent site-specific ¹H-¹⁵N HSQC data with coarse-grained molecular dynamics simulations. The results obtained using this equilibrium approach demonstrate that the main barrier to folding of pp32 is quite broad and lies near the unfolded state, with structure apparent only in the C-terminal region. Significant deviation from two-state unfolding under pressure reveals an intermediate on the folded side of the main barrier in which the N-terminal region is disordered. A nonlinear temperature dependence of the population of this intermediate suggests a large heat capacity change associated with its formation. The combination of pressure, which favors the population of folding intermediates relative to chemical denaturants; NMR, which allows their observation; and constrained structure-based simulations yield unparalleled insight into protein folding mechanisms.

α-chymotrypsin in water-acetone and water-dimethyl sulfoxide mixtures: Effect of preferential solvation and hydration

We investigated water/organic solvent sorption and residual enzyme activity to simultaneously monitor preferential solvation/hydration of protein macromolecules in the entire range of water content at 25°C. We applied this approach to estimate protein destabilization/stabilization due to the preferential interactions of bovine pancreatic α-chymotrypsin with water-acetone (moderate-strength H-bond acceptor) and water-DMSO (strong H-bond acceptor) mixtures. There are three concentration regimes for the dried α-chymotrypsin. α-Chymotrypsin is preferentially hydrated at high water content. The residual enzyme activity values are close to 100%. At intermediate water content, the dehydrated α-chymotrypsin has a higher affinity for acetone/DMSO than for water. Residual enzyme activity is minimal in this concentration range. The acetone/DMSO molecules are preferentially excluded from the protein surface at the lowest water content, resulting in preferential hydration. The residual catalytic activity in the water-poor acetone is ∼80%, compared with that observed after incubation in pure water. This effect is very small for the water-poor DMSO. Two different schemes are operative for the hydrated enzyme. At high and intermediate water content, α-chymotrypsin exhibits preferential hydration. However, at intermediate water content, in contrast to the dried enzyme, the initially hydrated α-chymotrypsin possesses increased preferential hydration parameters. At low water content, no residual enzyme activity was observed. Preferential binding of DMSO/acetone to α-chymotrypsin was detected. Our data clearly demonstrate that the hydrogen bond accepting ability of organic solvents and the protein hydration level constitute key factors in determining the stability of protein-water-organic solvent systems.

Structural and Thermodynamic Properties of Septin 3 Investigated by Small-Angle X-Ray Scattering

Septins comprise a family of proteins involved in a variety of cellular processes and related to several human pathologies. They are constituted by three structural domains: the N- and C-terminal domains, highly variable in length and composition, and the central domain, involved in the guanine nucleotide (GTP) binding. Thirteen different human septins are known to form heterogeneous complexes or homofilaments, which are stabilized by specific interactions between the different interfaces present in the domains. In this work, we have investigated by in-solution small-angle x-ray scattering the structural and thermodynamic properties of a human septin 3 construct, SEPT3-GC, which contains both of both interfaces (G and NC) responsible for septin-septin interactions. In order to shed light on the role of these interactions, small-angle x-ray scattering measurements were performed in a wide range of temperatures, from 2 up to 56°C, both with and without a nonhydrolysable form of GTP (GTPγS). The acquired data show a temperature-dependent coexistence of monomers, dimers, and higher-order aggregates that were analyzed using a global fitting approach, taking into account the crystallographic structure of the recently reported SEPT3 dimer, PDB:3SOP. As a result, the enthalpy, entropy, and heat capacity variations that control the dimer-monomer dissociation equilibrium in solution were derived and GTPγS was detected to increase the enthalpic stability of the dimeric species. Moreover, a temperature increase was observed to induce dissociation of SEPT3-GC dimers into monomers just preceding their reassembling into amyloid aggregates, as revealed by the Thioflavin-T fluorescence assays.

The Principles of Buoyancy in Marine Fish Eggs and Their Vertical Distributions across the World Oceans

Buoyancy acting on plankton, i.e. the difference in specific gravity between plankton and the ambient water, is a function of salinity and temperature. From specific gravity measurements of marine fish eggs salinity appears to be the only determinant of the buoyancy indicating that the thermal expansions of the fish egg and the ambient seawater are equal. We analyze the mechanisms behind thermal expansion in fish eggs in order to determine to what extent it can be justified to neglect the effects of temperature on buoyancy. Our results confirm the earlier assumptions that salinity is the basic determinant on buoyancy in marine fish eggs that, in turn, influence the vertical distributions and, consequently, the dispersal of fish eggs from the spawning areas. Fish populations have adapted accordingly by producing egg specific gravities that tune the egg buoyancy to create specific vertical distributions for each local population. A wide variety of buoyancy adaptations are found among fish populations. The ambient physical conditions at the spawning sites form a basic constraint for adaptation. In coastal regions where salinity increases with depth, and where the major fraction of the fish stocks spawns, pelagic and mesopelagic egg distributions dominate. However, in the larger part of worlds’ oceans salinity decreases with depth resulting in different egg distributions. Here, the principles of vertical distributions of fish eggs in the world oceans are presented in an overarching framework presenting the basic differences between regions, mainly coastal, where salinity increases with depth and the major part of the world oceans where salinity decreases with depth. We show that under these latter conditions, steady-state vertical distribution of mesopelagic fish eggs cannot exist as it does in most coastal regions. In fact, a critical spawning depth must exist where spawning below this depth threshold results in eggs sinking out of the water column and become lost for recruitment to the population. An example of adaptation to such conditions is Cape hake spawning above the critical layer in the Northern Benguela upwelling ecosystem. The eggs rise slowly in the onshore subsurface current below the Ekman layer, hence being advected inshore where the hatched larvae concentrate with optimal feeding conditions.

Physicochemical origin of high correlation between thermal stability of a protein and its packing efficiency: a theoretical study for staphylococcal nuclease mutants

There is an empirical rule that the thermal stability of a protein is related to the packing efficiency or core volume of the folded state and the protein tends to exhibit higher stability as the backbone and side chains are more closely packed. Previously, the wild type and its nine mutants of staphylococcal nuclease were compared by examining their folded structures. The results obtained were as follows: The stability was not correlated with the number of intramolecular hydrogen bonds, intramolecular electrostatic interaction energy, or degree of burial of the hydrophobic surface; though the empirical rule mentioned above held, it was not the proximate cause of higher stability; and the number of van der Waals contacts N vdW, or equivalently, the intramolecular van der Waals interaction energy was an important factor governing the stability. Here we revisit the wild type and its nine mutants of staphylococcal nuclease using our statistical-mechanical theory for hydration of a protein. A molecular model is employed for water. We show that the pivotal factor is the magnitude of the water-entropy gain upon folding. The gain originates from an increase in the total volume available to the translational displacement of water molecules coexisting with the protein in the system. The magnitude is highly correlated with the denaturation temperature T m. Moreover, the apparent correlation between N vdW and T m as well as the empirical rule is interpretable (i.e., their physicochemical meanings can be clarified) on the basis of the water-entropy effect.

Evolutionarily Conserved Pattern of Interactions in a Protein Revealed by Local Thermal Expansion Properties

The way in which the network of intramolecular interactions determines the cooperative folding and conformational dynamics of a protein remains poorly understood. High-pressure NMR spectroscopy is uniquely suited to examine this problem because it combines the site-specific resolution of the NMR experiments with the local character of pressure perturbations. Here we report on the temperature dependence of the site-specific volumetric properties of various forms of staphylococcal nuclease (SNase), including three variants with engineered internal cavities, as measured with high-pressure NMR spectroscopy. The strong temperature dependence of pressure-induced unfolding arises from poorly understood differences in thermal expansion between the folded and unfolded states. A significant inverse correlation was observed between the global thermal expansion of the folded proteins and the number of strong intramolecular hydrogen bonds, as determined by the temperature coefficient of the backbone amide chemical shifts. Comparison of the identity of these strong H-bonds with the co-evolution of pairs of residues in the SNase protein family suggests that the architecture of the interactions detected in the NMR experiments could be linked to a functional aspect of the protein. Moreover, the temperature dependence of the residue-specific volume changes of unfolding yielded residue-specific differences in expansivity and revealed how mutations impact intramolecular interaction patterns. These results show that intramolecular interactions in the folded states of proteins impose constraints against thermal expansion and that, hence, knowledge of site-specific thermal expansivity offers insight into the patterns of strong intramolecular interactions and other local determinants of protein stability, cooperativity, and potentially also of function.

Wanted: scalable tracers for diffusion measurements

Scalable tracers are potentially a useful tool to examine diffusion mechanisms and to predict diffusion coefficients, particularly for hindered diffusion in complex, heterogeneous, or crowded systems. Scalable tracers are defined as a series of tracers varying in size but with the same shape, structure, surface chemistry, deformability, and diffusion mechanism. Both chemical homology and constant dynamics are required. In particular, branching must not vary with size, and there must be no transition between ordinary diffusion and reptation. Measurements using scalable tracers yield the mean diffusion coefficient as a function of size alone; measurements using nonscalable tracers yield the variation due to differences in the other properties. Candidate scalable tracers are discussed for two-dimensional (2D) diffusion in membranes and three-dimensional diffusion in aqueous solutions. Correlations to predict the mean diffusion coefficient of globular biomolecules from molecular mass are reviewed briefly. Specific suggestions for the 3D case include the use of synthetic dendrimers or random hyperbranched polymers instead of dextran and the use of core-shell quantum dots. Another useful tool would be a series of scalable tracers varying in deformability alone, prepared by varying the density of crosslinking in a polymer to make say "reinforced Ficoll" or "reinforced hyperbranched polyglycerol."

Probing volumetric properties of biomolecular systems by pressure perturbation calorimetry (PPC)--the effects of hydration, cosolvents and crowding

Pressure perturbation calorimetry (PPC) is an efficient technique to study the volumetric properties of biomolecules in solution. In PPC, the coefficient of thermal expansion of the partial volume of the biomolecule is deduced from the heat consumed or produced after small isothermal pressure-jumps. The expansion coefficient strongly depends on the interaction of the biomolecule with the solvent or cosolvent as well as on its packing and internal dynamic properties. This technique, complemented with molecular acoustics and densimetry, provides valuable insights into the basic thermodynamic properties of solvation and volume effects accompanying interactions, reactions and phase transitions of biomolecular systems. After outlining the principles of the technique, we present representative examples on protein folding, including effects of cosolvents and crowding, together with a discussion of the interpretation, and further applications.