The Effect of Arginine on the Phase Stability of Aqueous Hen Egg-White Lysozyme Solutions.

The effect of arginine on the phase stability of the hen egg-white lysozyme (HEWL) has been studied via molecular dynamics computer simulations, as well as experimentally via cloud-point temperature determination. The experiments show that the addition of arginine increases the stability of the HEWL solutions. The computer simulation results indicate that arginine molecules tend to self-associate. If arginine residues are located on the protein surface, the free arginine molecules stay in their vicinity and prevent the way protein molecules "connect" through them to form clusters. The results are not sensitive to a particular force field and suggest a possible microscopic mechanism of the stabilizing role of arginine as an excipient.

The mechanism of self-association of human γ-D crystallin from molecular dynamics simulations.

The aggregation of human γ-D crystallin is associated with the age-onset cataract formation. Here, we extensively investigated the self-association mechanism of human γ-D crystallin through molecular dynamics computer simulations. By mutating the protein surface we found that electrostatic interactions between charged amino acids play a crucial role in its self-association. We have confirmed the two-fold role of arginine molecules. If they are located as residues on the protein surface they can initiate protein contacts and contribute to their stickiness with noteworthy hydrophobic interactions through stacking of their methylene groups. But if they are added as free arginine in the protein solution they can also stabilize it, by associating with the protein surface and also with themselves to form effective inter-protein spacers that obstruct protein aggregation.

BioMThermDB 1.0: Thermophysical Database of Proteins in Solutions.

We present here a freely available web-based database, called BioMThermDB 1.0, of thermophysical and dynamic properties of various proteins and their aqueous solutions. It contains the hydrodynamic radius, electrophoretic mobility, zeta potential, self-diffusion coefficient, solution viscosity, and cloud-point temperature, as well as the conditions for those determinations and details of the experimental method. It can facilitate the meta-analysis and visualization of data, can enable comparisons, and may be useful for comparing theoretical model predictions with experiments.

Influence of Low Molecular Weight Salts on the Viscosity of Aqueous-Buffer Bovine Serum Albumin Solutions.

Pharmaceutical design of protein formulations aims at maximum efficiency (protein concentration) and minimum viscosity. Therefore, it is important to know the nature of protein-protein interactions and their influence on viscosity. In this work, we investigated the dependence of the viscosity of BSA in an aqueous 20 mM acetate buffer at pH = 4.3 on protein concentration and on temperature (5-45 °C). The viscosity of the solution increased with protein concentration and was 230% higher than the viscosity of the protein-free formulation at 160 mg/mL. The viscosity decreased by almost 60% in the temperature range from 5 to 45 °C. The agreement of the modified Arrhenius theory with experiment was quantitative, whereas a hard-sphere model provided only a qualitative description of the experimental results. We also investigated the viscosity of a 100 mg/mL BSA solution as a function of the concentration of added low molecular weight salts (LiCl, NaCl, KCl, RbCl, CsCl, NaBr, NaI) in the range of salt concentrations up to 1.75 mol/L. In addition, the particle size and zeta potential of BSA-salt mixtures were determined for solutions containing 0.5 mol/L salt. The trends with respect to the different anions followed a direct Hofmeister series (Cl- > Br- > I-), whereas for cations in the case of viscosity the indirect Hofmeister series was observed (Li+ > Na+ > K+ > Rb+ > Cs+), but the values of particle sizes and zeta potential did not show cation-specific effects. Since the protein is positively charged at pH = 4.3, anions are more attracted to the protein surface and shield its charge, while the interaction with cations is less pronounced. We hypothesize that salt surface charge shielding reduces protein colloidal stability and promotes protein aggregate formation.

The influence of excipients on the viscosity of monoclonal antibody solutions.

The aggregation propensity of monoclonal antibodies can be modified by adding different cosolutes into the solution. A simple coarse-grained model in the combination with the thermodynamic perturbation theory was used to predict cluster distribution and viscosity of the solutions of IgG4 monoclonal anibody in the presence of L-Arginine Hydrochloride. The data were analysed using binding polynomial to describe the binding of cosolute (Arginine) to the antibody molecule. The results show that by binding to the antibody molecule the cosolute occupies some of the binding sites of the antibody, and in this way reduces the amount of binding sites available to other antibody molecules. The aggregation propensity of the antibody molecules is therefore reduced.

Studying the mechanism of phase separation in aqueous solutions of globular proteins via molecular dynamics computer simulations.

Proteins are the most abundant biomacromolecules in living cells, where they perform vital roles in virtually every biological process. To maintain their function, proteins need to remain in a stable (native) state. Inter- and intramolecular interactions in aqueous protein solutions govern the fate of proteins, as they can provoke their unfolding or association into aggregates. The initial steps of protein aggregation are difficult to capture experimentally, therefore we used molecular dynamics simulations in this study. We investigated the initial phase of aggregation of two different lysozymes, hen egg-white (HEWL) and T4 WT* lysozyme and also human lens γ-D crystallin by using atomistic simulations. We monitored the phase stability of their aqueous solutions by calculating time-dependent density fluctuations. We found that all proteins remained in their compact form despite aggregation. With an extensive analysis of intermolecular residue-residue interactions we discovered that arginine is of paramount importance in the initial stage of aggregation of HEWL and γ-D crystallin, meanwhile lysine was found to be the most involved amino acid in forming initial contacts between T4 WT* molecules.

How Water's Properties Are Encoded in Its Molecular Structure and Energies.

How are water's material properties encoded within the structure of the water molecule? This is pertinent to understanding Earth's living systems, its materials, its geochemistry and geophysics, and a broad spectrum of its industrial chemistry. Water has distinctive liquid and solid properties: It is highly cohesive. It has volumetric anomalies-water's solid (ice) floats on its liquid; pressure can melt the solid rather than freezing the liquid; heating can shrink the liquid. It has more solid phases than other materials. Its supercooled liquid has divergent thermodynamic response functions. Its glassy state is neither fragile nor strong. Its component ions-hydroxide and protons-diffuse much faster than other ions. Aqueous solvation of ions or oils entails large entropies and heat capacities. We review how these properties are encoded within water's molecular structure and energies, as understood from theories, simulations, and experiments. Like simpler liquids, water molecules are nearly spherical and interact with each other through van der Waals forces. Unlike simpler liquids, water's orientation-dependent hydrogen bonding leads to open tetrahedral cage-like structuring that contributes to its remarkable volumetric and thermal properties.

Protein aggregation in salt solutions.

Protein aggregation is broadly important in diseases and in formulations of biological drugs. Here, we develop a theoretical model for reversible protein-protein aggregation in salt solutions. We treat proteins as hard spheres having square-well-energy binding sites, using Wertheim's thermodynamic perturbation theory. The necessary condition required for such modeling to be realistic is that proteins in solution during the experiment remain in their compact form. Within this limitation our model gives accurate liquid-liquid coexistence curves for lysozyme and γ IIIa-crystallin solutions in respective buffers. It provides good fits to the cloud-point curves of lysozyme in buffer-salt mixtures as a function of the type and concentration of salt. It than predicts full coexistence curves, osmotic compressibilities, and second virial coefficients under such conditions. This treatment may also be relevant to protein crystallization.

Volumetric and calorimetric properties of aqueous ionene solutions.

The volumetric (partial and apparent molar volumes) and calorimetric properties (apparent heat capacities) of aqueous cationic polyelectrolyte solutions - ionenes - were studied using the oscillating tube densitometer and differential scanning calorimeter. The polyion's charge density and the counterion properties were considered as variables. The special attention was put to evaluate the contribution of electrostatic and hydrophobic effects to the properties studied. The contribution of the CH2 group of the polyion's backbone to molar volumes and heat capacities was estimated. Synergistic effect between polyion and counterions was found.

Phase behaviour of a continuous shouldered well model fluid. A grand canonical Monte Carlo study.

The phase behavior of the continuous shouldered well model fluid proposed by Franzese [J. Mol. Liq. 136 (2007) 267] was examined using the Monte Carlo computer simulations in the grand canonical ensemble. The essential parts of the vapour-liquid and liquid-liquid coexistence envelopes were obtained. The Widom lines departing from coexistence envelopes were calculated using maxima of the fluctuations of the number of particles as a function of chemical potential along various isotherms. The region embracing anomalies in the properties of the model was located using the approximate criterion that involves the excess pair entropy.. The temperature of maximum density line was built by performing canonical Monte Carlo simulations. Our results are consistent with previous results from molecular dynamics constant pressure-constant temperature simulations and provide wider insight into the phase behavior of the model by using the chemical potential as the external parameter.

Modelling the Correlation Between Molecular Electrostatic Potential and pKa on Sets of Carboxylic Acids, Phenols and Anilines.

Calculations of molecular electrostatic potential were correlated with experimental pKa values for different sets of acidic molecules (carboxylic acids, phenols, and anilines) to obtain linear relationships of variable quality. A single tri-parameter model function was constructed to describe the pKa dependence on MEP maxima together with two automatically generated molecular descriptors, namely the counts of carboxylic acid and amine functional groups.

Effects of the translational and rotational degrees of freedom on the hydration of simple solutes.

Molecular dynamics simulations with separate thermostats for rotational and translational motion were used to study the effect of these degrees of freedom on the structure of water around model solutes. To describe water molecules we used the SPC/E model. The simplest solute studied here, the hydrophobe, was represented as a Lennard-Jones particle. Since direct interaction between the hydrophobe and water molecules has no angular dependence the influence of the increase of the rotational temperature on the solvation of a hydrophobe is only indirect. In the next step the central solute was assumed to be charged with either a positive or a negative charge to mimic an ion in water. Hence, depending on the charge of the ion, the neighboring water molecules assumed different angular distributions. The principal conclusions of this work are: (i) an increase of the translational temperature always decreases the height of the first peak in the solute-water radial distribution function; (ii) an increase of the rotational temperature yields an increase in the first peak in the solute-water radial distribution function for hydrophobes and cations; (iii) in contrast to this, the solvation peak decreases around ions with sufficiently large negative charge; and (iv) an increase of the rotational temperature affects cations in an opposite way to anions. For this reason complex molecules with a small net charge may not be very sensitive to variation of the rotational temperature.

The application of the integral equation theory to study the hydrophobic interaction.

The Wertheim's integral equation theory was tested against newly obtained Monte Carlo computer simulations to describe the potential of mean force between two hydrophobic particles. An excellent agreement was obtained between the theoretical and simulation results. Further, the Wertheim's integral equation theory with polymer Percus-Yevick closure qualitatively correctly (with respect to the experimental data) describes the solvation structure under conditions where the simulation results are difficult to obtain with good enough accuracy.

An Investigation of Ion-Pairing of Alkali Metal Halides in Aqueous Solutions Using the Electrical Conductivity and the Monte Carlo Computer Simulation Methods.

The ion pairing is, in very dilute aqueous solutions, of rather small importance for solutions' properties, which renders its precise quantification quite a laborious task. Here we studied the ion pairing of alkali halides in water by using the precise electric conductivity measurements in dilute solutions, and in a wide temperature range. The low-concentration chemical model was used to analyze the results, and to estimate the association constant of different alkali halide salts. It has been shown that the association constant is related to the solubility of salts in water and produces a 'volcano relationship', when plotted against the difference between the free energy of hydration of the corresponding individual ions. The computer simulation, using the simple MB+dipole water model, were used to interprete the results, to find a microscopic basis for Collins' law of matching water affinities.

Antibody association in solution: cluster distributions and mechanisms.

Understanding factors that affect the clustering and association of antibodies molecules in solution is critical to their development as therapeutics. For 19 different monoclonal antibody (mAb) solutions, we measured the viscosities, the second virial coefficients, the Kirkwood-Buff integrals, and the cluster distributions of the antibody molecules as functions of protein concentration. Solutions were modeled using the statistical-physics Wertheim liquid-solution theory, representing antibodies as Y-shaped molecular structures of seven beads each. We found that high-viscosity solutions result from more antibody molecules per cluster. Multi-body properties such as viscosity are well predicted experimentally by the 2-body Kirkwood-Buff quantity, G22, but not by the second virial coefficient, B22, and well-predicted theoretically from the Wertheim protein-protein sticking energy. Weakly interacting antibodies are rate-limited by nucleation; strongly interacting ones by propagation. This approach gives a way to relate micro to macro properties of solutions of associating proteins.

The application of the thermodynamic perturbation theory to study the hydrophobic hydration.

The thermodynamic perturbation theory was tested against newly obtained Monte Carlo computer simulations to describe the major features of the hydrophobic effect in a simple 3D-Mercedes-Benz water model: the temperature and hydrophobe size dependence on entropy, enthalpy, and free energy of transfer of a simple hydrophobic solute into water. An excellent agreement was obtained between the theoretical and simulation results. Further, the thermodynamic perturbation theory qualitatively correctly (with respect to the experimental data) describes the solvation thermodynamics under conditions where the simulation results are difficult to obtain with good enough accuracy, e.g., at high pressures.

The Role of Buffers in Wild-Type HEWL Amyloid Fibril Formation Mechanism: A Methodological Approach.

The amyloidophilic dyes thioflavin T and Congo red are small, yet powerful, molecules that allow the in vitro and in vivo detection of amyloid fibrils in protein solutions. Even though Congo red and thioflavin T binding assays are widespread techniques for unveiling amyloid fibers and are gradually replacing the more demanding X-ray diffraction method, handling samples containing amyloid fibrils is still challenging and can lead to false-positive/negative results. Here we describe a relatively straightforward procedure of preparing hen egg-white lysozyme amyloid fibrils in different buffer solutions and their detection with thioflavin T and Congo red, supported by an indispensable method for determining the secondary structure of proteins - circular dichroism.

Biophysical Principles Emerging from Experiments on Protein-Protein Association and Aggregation.

Protein-protein association and aggregation are fundamental processes that play critical roles in various biological phenomena, from cellular signaling to disease progression. Understanding the underlying biophysical principles governing these processes is crucial for elucidating their mechanisms and developing strategies for therapeutic intervention. In this review, we provide an overview of recent experimental studies focused on protein-protein association and aggregation. We explore the key biophysical factors that influence these processes, including protein structure, conformational dynamics, and intermolecular interactions. We discuss the effects of environmental conditions such as temperature, pH and related buffer-specific effects, and ionic strength and related ion-specific effects on protein aggregation. The effects of polymer crowders and sugars are also addressed. We list the techniques used to study aggregation. We analyze emerging trends and challenges in the field, including the development of computational models and the integration of multidisciplinary approaches for a comprehensive understanding of protein-protein association and aggregation. Expected final online publication date for the Annual Review of Biophysics, Volume 53 is May 2024. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Protein Association in Solution: Statistical Mechanical Modeling.

Protein molecules associate in solution, often in clusters beyond pairwise, leading to liquid phase separations and high viscosities. It is often impractical to study these multi-protein systems by atomistic computer simulations, particularly in multi-component solvents. Instead, their forces and states can be studied by liquid state statistical mechanics. However, past such approaches, such as the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, were limited to modeling proteins as spheres, and contained no microscopic structure-property relations. Recently, this limitation has been partly overcome by bringing the powerful Wertheim theory of associating molecules to bear on protein association equilibria. Here, we review these developments.

Structural and thermodynamical properties of charged hard spheres in a mixture with core-softened model solvent.

The canonical Monte Carlo computer simulations and integral equation theory were applied to examine the structural and thermodynamic properties of a mixture of ions and a core-softened fluid molecules. The positive and negative ions forming a +1:-1 salt were modeled as charged hard spheres, immersed in the dielectric medium. It was shown previously that the core-softened fluid under study is characterized by a set of structural, thermodynamic, and dynamic anomalies. The principal objective of this work was to elucidate how the presence of ions alters this behavior. The structural properties of the mixtures are discussed in terms of the pair distribution functions; in addition, the pair contribution to the excess entropy was calculated. Thermodynamic properties are investigated by using the dependencies of energy and compressibility factor on density, composition of the mixture, and reduced temperature. The heat capacity was also evaluated. Our principal findings concern the description of structural anomalies in the mixture, the dependence of the temperature of maximum density on the ionic concentration, and establishing the regions delimiting the structural and thermodynamic anomalies of the model mixture.