Peptide model systems: Correlation between thermophilicity and hydrophilicity

Deciphering the guanidinium cation: Insights into thermal diffusion

Thermophoresis, or thermodiffusion, is becoming a more popular method for investigating the interactions between proteins and ligands due to its high sensitivity to the interactions between solutes and water. Despite its growing use, the intricate mechanisms behind thermodiffusion remain unclear. This gap in knowledge stems from the complexities of thermodiffusion in solvents that have specific interactions as well as the intricate nature of systems that include many components with both non-ionic and ionic groups. To deepen our understanding, we reduce complexity by conducting systematic studies on aqueous salt solutions. In this work, we focused on how guanidinium salt solutions behave in a temperature gradient, using thermal diffusion forced Rayleigh scattering experiments at temperatures ranging from 15 to 35 °C. We looked at the thermodiffusive behavior of four guanidinium salts (thiocyanate, iodide, chloride, and carbonate) in solutions with concentrations ranging from 1 to 3 mol/kg. The guanidinium cation is disk-shaped and is characterized by flat hydrophobic surfaces and three amine groups, which enable directional hydrogen bonding along the edges. We compare our results to the behavior of salts with spherical cations, such as sodium, potassium, and lithium. Our discussions are framed around how different salts are solvated, specifically in the context of the Hofmeister series, which ranks ions based on their effects on the solvation of proteins.

Non-monotonic Soret coefficients of aqueous LiCl solutions with varying concentrations

We investigate the thermodiffusive properties of aqueous solutions of lithium chloride, using thermal diffusion forced Rayleigh scattering in a concentration range of 0.5-2 mole per kg of solvent and a temperature range of 5 to 45 °C. All solutions exhibit non-monotonic variations of the Soret coefficient ST with a concentration exhibiting a minimum at about one mole per kg of solvent. The depth of the minimum decreases with increasing temperature and shifts slightly towards higher concentrations. We compare the experimental data with published data and apply a recent model based on overlapping hydration shells. Additionally, we calculate the ratio of the phenomenological Onsager coefficients using our experimental results and published data to calculate the thermodynamic factor. Simple linear, quadratic and exponential functions can be used to describe this ratio accurately, and together with the thermodynamic factors, the experimental Soret coefficients can be reproduced. The main conclusion from this analysis is that the minimum of the Soret coefficients results from a maximum in the thermodynamic factor, which appears itself at concentrations far below the experimental concentrations. Only after multiplication by the (negative) monotonous Onsager ratio does the minimum move into the experimental concentration window.

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.

Mass effects for thermodiffusion in dilute aqueous solutions

Thermodiffusion is the phenomenon by which molecules in a mixture present concentration gradients in response to an imposed temperature gradient. Despite decades of investigations, this effect remains poorly understood at a molecular level. A common, phenomenological approach is to individuate the molecular factors that influence the Soret coefficient, the parameter that quantifies the resulting concentration-gradient. Experimental studies, often performed on organic mixtures, as well as simulations of model particle systems have evidenced that the difference in masses between the mixture components has an important effect on the amplitude of the Soret coefficient. Here, we use molecular dynamics simulations of a thermophoretic setting to investigate the mass dependence of the Soret coefficient in dilute aqueous solutions. An advantage of simulation approaches is that they are not limited in the range of explored molecular masses, which is often limited to isotopic substitutions in the experiments. Our simulations reveal that the mass dependence of the Soret coefficient in these solutions is in agreement with previous experimental and simulation work on molecular-size systems. In particular, it is sensitive to the relative mass difference between the solute and the solvent, but not to their absolute mass. Adjusting the mass of the solvent and of the solute can turn a thermophobic solution into a thermophilic one, where solute accumulation is reversed. This demonstrates that the mass effect can indeed compensate for the other contributions to the Soret coefficient. Finally, we find that changing the molecular moments of inertia has a much more limited impact as compared to a change in the total molecular mass.

Towards understanding specific ion effects in aqueous media using thermodiffusion

Specific ion effects play an important role in scientific and technological processes. According to Hofmeister, the influence on the hydrogen bond network depends on the ion and leads to a specific order of the ions. Also thermodiffusion the mass transport caused by a temperature gradient is very sensitive to changes of the hydrogen bond network leading to a ranking according to hydrophilicity of the salt. Hence, we investigate various salt solutions in order to compare with the Hofmeister concept. We have studied three different sodium salts in water as a function of temperature (25-45[Formula: see text]C) and concentration (0.5-5 mol kg[Formula: see text]) using Thermal Diffusion Forced Rayleigh Scattering (TDFRS). The three anions studied, carbonate, acetate and thiocyanate, span the entire range of the Hofmeister series from hydrophilic to hydrophobic. We compare the results with the recent measurements of the corresponding potassium salts to see to what extent the cation changes the thermodiffusion of the salt.

Overlapping hydration shells in salt solutions causing non-monotonic Soret coefficients with varying concentration

We develop an intuitive picture that overlapping hydration shells in salt solutions cause non-monotonic Soret coefficients with varying concentration.

Molecular Dynamics Approach to Calculate the Thermodiffusion (Soret and Seebeck) Coefficients of Salts in Aqueous Solutions

An approach to investigate the physical parameters related to ion thermodiffusion in aqueous solutions is proposed herein by calculating the equilibrium hydration free energy and the self-diffusion coefficient as a function of temperature, ranging from 293 to 353 K, using molecular dynamics simulations of infinitely diluted ions in aqueous solutions. Several ion force field parameters are used in the simulations, and new parameters are proposed for some ions to better describe their hydration free energy. Such a theoretical framework enables the calculation of some single-ion properties, such as heat of transport, Soret coefficient, and mass current density, as well as properties of salts, such as effective mass and thermal diffusion, Soret and Seebeck, coefficients. These calculated properties are compared with experimental data available from optical measurements and showed good agreement revealing an excellent theoretical predictability of salt thermodiffusion properties. Differences in single-ion Soret and self-diffusion coefficients of anions and cations give rise to a thermoelectric field, which affects the system response that is quantified by the Seebeck coefficient. The fast and slow Seebeck coefficients are calculated and discussed, resulting in values with mV/K order of magnitude, as observed in experiments involving several salts, such as K⁺Cl^-, Na⁺Cl^-, H⁺Cl^-, Na⁺OH^-, TMA⁺OH^-, and TBA⁺OH^-. The present approach can be adopted for any ion or charged particle dispersed in water with the aim of predicting the thermoelectric field induced through the fluid. It has potential applications in designing electrolytes for ionic thermoelectric devices in order to harvest energy and thermoelectricity in biological nanofluids.

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.

Thermophoresis of biological and biocompatible compounds in aqueous solution

With rising popularity of microscale thermophoresis for the characterisation of protein-ligand binding reactions and possible applications in microfluidic devices, there is a growing interest in considering thermodiffusion in the context of life sciences. But although the understanding of thermodiffusion in non-polar mixtures has grown rapidly in recent years, predictions for associated mixtures like aqueous solutions remain challenging. This review aims to give an overview of the literature on thermodiffusion in aqueous systems, show the difficulties in theoretical description that arise from the non-ideal behaviour of water-mixtures, and highlight the relevance of thermodiffusion in a biological context. We find that the thermodiffusion in aqueous systems is dominated by contributions from heat of transfer, hydrogen bond interactions and charge effects. However, the separation of these effects is often difficult, especially in case of biological systems where a systematic exclusion of contributions may not be feasible.

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.

Configurational contribution to the Soret effect of a protein ligand system : An investigation with density-of-states simulations

Many of the biological functions of proteins are closely associated with their ability to bind ligands and change conformations in response to changing conditions. Since binding state and conformation of a protein affect its response to a temperature gradient, they may be probed with thermophoresis. In recent years, thermophoretic techniques to investigate biomolecular interactions, quantify ligand binding, and probe conformational changes have become established. To develop a better understanding of the mechanisms underlying the thermophoretic behavior of proteins and ligands, we employ a simple, off-lattice model for a protein and ligand in explicit solvent. To investigate the partitioning of the particles in a temperature gradient, we perform Wang-Landau-type simulations in a divided simulation box and construct the density of states over a two-dimensional state space. This method gives us access to the entropy and energy of the divided system and allows us to estimate the configurational contribution to the Soret coefficient. In this work, we focus on dilute solutions of hydrophobic proteins and investigate the effect of ligand binding on their thermophoretic behavior. We find that our simple model captures important aspects of protein-ligand interactions and allows us to relate the binding energy to the change in Soret coefficient upon ligand binding.

Correlation between thermophoretic behavior and hydrophilicity for various alcohols⋆

Recent experiments for various amides and sugars showed a clear correlation of the temperature dependence of the Soret coefficient with the hydrophilicity, quantitatively described by the logarithm of the 1-octanol/water partition coefficient log P . This coefficient is a measure for the hydrophilicity/hydrophobicity balance of a solute and is often used to model the transport of a compound in the environment or to screen for potential pharmaceutical compounds. In order to validate whether this concept works also for other water soluble molecules we investigated systematically the thermophoresis of mono- and polyhydric alcohols. As experimental method we use a holographic grating technique called infrared Thermal Diffusion Forced Rayleigh Scattering (IR-TDFRS). Experiments showed that the temperature dependence of the Soret coefficient of polyhydric alcohols also correlates with log P and lies on the same master plot as amides and sugars.

Thermophoresis of a Colloidal Rod: Contributions of Charge and Grafted Polymers

In this study, we investigated the thermodiffusion behavior of a colloidal model system as a function of the Debye length, λ_DH, which is controlled by the ionic strength. Our system consists of an fd-virus grafted with poly(ethylene glycol) (PEG) with a molecular mass of 5000 g mol^-1. The results are compared with recent measurements on a bare fd-virus and with results of PEG. The diffusion coefficients of both viruses are comparable and increase with the increasing Debye length. The thermal diffusion coefficient, D_T, of the bare virus increases strongly with the Debye length, whereas D_T of the grafted fd-virus shows only a very weak increase. The Debye length dependence of both systems can be described with an expression derived for charged rods using the surface charge density and an offset of D_T as adjustable parameters. It turns out that the ratio of the determined surface charges is inverse to the ratio of the surfaces of the two systems, which means that the total charge remains almost constant. The determined offset of the grafted fd-virus describing the chemical contributions is the sum of D_T of PEG and the offset of the bare fd-virus. At high λ_DH, corresponding to the low ionic strength, the S_T values of both colloidal model systems approach each other. This implies a contribution from the polymer layer, which is strong at short λ_DH and fades out for the longer Debye lengths, when the electric double layer reaches further than the polymer chains and therefore dominates interactions with the surrounding water.