The hydration problem in solution biophysics: an introduction

The molecular basis for hydrodynamic properties of PEGylated human serum albumin

Polyethylene glycol (PEG) conjugation provides a protective modification that enhances the pharmacokinetics and solubility of proteins for therapeutic use. A knowledge of the structural ensemble of these PEGylated proteins is necessary to understand the molecular details that contribute to their hydrodynamic and colligative properties. Because of the large size and dynamic flexibility of pharmaceutically important PEGylated proteins, the determination of structure is challenging. In addition, the hydration of these conjugates that contain large polymers is difficult to determine with traditional methods that identify only first shell hydration water, which does not account for the complete hydrodynamic volume of a macromolecule. Here, we demonstrate that structural ensembles, generated by coarse-grained simulations, can be analyzed with HullRad and used to predict sedimentation coefficients and concentration-dependent hydrodynamic and diffusion nonideality coefficients of PEGylated proteins. A knowledge of these concentration-dependent properties enhances the ability to design and analyze new modified protein therapeutics. HullRad accomplishes this analysis by effectively accounting for the complete hydration of a macromolecule, including that of flexible polymers.

Hydrodynamic and thermodynamic analysis of PEGylated human serum albumin

Covalent labeling of therapeutic drugs and proteins with polyethylene glycol (PEGylation) is an important modification for improving stability, solubility, and half-life. PEGylation alters protein solution behavior through its impact on thermodynamic nonideality by increasing the excluded volume, and on hydrodynamic nonideality by increasing the frictional drag. To understand PEGylation's impact, we investigated the thermodynamic and hydrodynamic properties of a model system consisting of PEGylated human serum albumin derivatives using analytical ultracentrifugation (AUC) and dynamic light scattering (DLS). We constructed PEGylated human serum albumin derivatives of single, linear 5K, 10K, 20K, and 40K PEG chains and a single branched-chain PEG of 40K (2 × 20K). Sedimentation velocity (SV) experiments were analyzed using SEDANAL direct boundary fitting to extract ideal sedimentation coefficients so, hydrodynamic nonideality ks, and thermodynamic nonideality 2BM1SV terms. These quantities allow the determination of the Stokes radius Rs, the frictional ratio f/fo, and the swollen or entrained volume Vs/v, which measure size, shape, and solvent interaction. We performed sedimentation equilibrium experiments to obtain independent measurements of thermodynamic nonideality 2BM1SE. From DLS measurements, we determined the interaction parameter, kD, the concentration dependence of the apparent diffusion coefficient, D, and from extrapolation of D to c = 0 a second estimate of Rs. Rs values derived from SV and DLS measurements and ensemble model calculations (see complementary study) are then used to show that ks + kD = theoretical 2B22M1. In contrast, experimental BM1 values from SV and sedimentation equilibrium data collectively allow for similar analysis for protein-PEG conjugates and show that ks + kD = 1.02-1.07∗BM1, rather than the widely used ks + kD = 2BM1 developed for hard spheres. The random coil behavior of PEG dominates the colloidal properties of PEG-protein conjugates and exceeds the sum of a random coil and hard-sphere volume due to excess entrained water.

PEG-Lipids: Quantitative Study of Unimers and Aggregates Thereof by the Methods of Molecular Hydrodynamics

Understanding the polymorphism of lipids in solution is the key to the development of intracellular delivery systems. Here, we study the dynamics of poly(ethylene glycol)-lipid (PEG-Lipid) conjugates aiming at a better understanding of their molecular properties and aggregation behavior in solution. Those PEG-Lipids are used as components of lipid nanoparticles (LNPs). LNPs are gaining increased popularity, e.g., by their utilization in modern vaccination strategies against SARS-CoV-2. Characterization of the systems is conducted by the classical methods of hydrodynamics in different solvents, such as ethanol and water, which are also commonly used for LNP formulation. We were able to elucidate the structurally associated hydrodynamic properties of isolated PEG-Lipids in ethanol, revealing the typically expected values of the hydrodynamic invariant for random coil polymers. By virtue of the same experimental setting, the PEG-Lipids' behavior in water was as well studied, which is a less good solvent than ethanol for the PEG-Lipids. Our experiments demonstrate that PEG-Lipids dissolved in water form well-defined micelles that can quantitatively be characterized in terms of their degree of aggregation of PEG-Lipid polymer unimers, their hydrodynamic size, and solvation, i.e., the quantitative determination of water contained or associated to the identified micelles. Quantitative results obtained from classical hydrodynamic analyses are fully supported by studies with standard dynamic light scattering (DLS). The obtained diffusion coefficients and hydrodynamic sizes are in excellent agreement with numerical results derived from analytical ultracentrifugation (AUC) data. Cryo-transmission electron microscopy (cryo-TEM) supports the structural insight from hydrodynamic studies, particularly, in terms of the observed spherical structure of the formed micelles. We demonstrate experimentally that the micelle systems can be considered as solvent-permeable, hydrated spheres.

Revisiting macromolecular hydration with HullRadSAS

Hydration of biological macromolecules is important for their stability and function. Historically, attempts have been made to describe the degree of macromolecular hydration using a single parameter over a narrow range of values. Here, we describe a method to calculate two types of hydration: surface shell water and entrained water. A consideration of these two types of hydration helps to explain the "hydration problem" in hydrodynamics. The combination of these two types of hydration allows accurate calculation of hydrodynamic volume and related macromolecular properties such as sedimentation and diffusion coefficients, intrinsic viscosities, and the concentration-dependent non-ideality identified with sedimentation velocity experiments.

Intrinsic viscosity of bovine serum albumin conformers

Intrinsic viscosity of bovine serum albumin (BSA) at different pH values (2.7, 4.3, 7.4, 8 and 10) has been determined, as well as the Mark-Houwink constant and expansion factor. The traditional technique for data analysis using extrapolation to obtain intrinsic viscosity values shows an unusual behavior regarding concentration that can be observed in the values obtained for Huggins' and Kraemer's constants.

Interfacial energetics of globular-blood protein adsorption to a hydrophobic interface from aqueous-buffer solution

Adsorption isotherms of nine globular proteins with molecular weight (MW) spanning 10-1000 kDa confirm that interfacial energetics of protein adsorption to a hydrophobic solid/aqueous-buffer (solid-liquid, SL) interface are not fundamentally different than adsorption to the water-air (liquid-vapour, LV) interface. Adsorption dynamics dampen to a steady-state (equilibrium) within a 1 h observation time and protein adsorption appears to be reversible, following expectations of Gibbs' adsorption isotherm. Adsorption isotherms constructed from concentration-dependent advancing contact angles theta(a) of buffered-protein solutions on methyl-terminated, self-assembled monolayer surfaces show that maximum advancing spreading pressure, Pi(a)max, falls within a relatively narrow 10 < Pi(a)max < 20 mN m(-1) band characteristic of all proteins studied, mirroring results obtained at the LV surface. Furthermore, Pi(a) isotherms exhibited a 'Traube-rule-like' progression in MW similar to the ordering observed at the LV surface wherein molar concentrations required to reach a specified spreading pressure Pi(a) decreased with increasing MW. Finally, neither Gibbs' surface excess quantities [Gamma(sl)-Gamma(sv)] nor Gamma(lv) varied significantly with protein MW. The ratio {[Gamma(sl)-Gamma(sv)]/Gamma(lv)} approximately 1, implying both that Gamma(sv) approximately 0 and chemical activity of protein at SL and LV surfaces was identical. These results are collectively interpreted to mean that water controls protein adsorption to hydrophobic surfaces and that the mechanism of protein adsorption can be understood from this perspective for a diverse set of proteins with very different composition.

Fluidity of bound hydration layers

We have measured the shear forces between solid surfaces sliding past each other across aqueous salt solutions, at pressures and concentrations typical of naturally occurring systems. In such systems the surface-attached hydration layers keep the compressed surfaces apart as a result of strongly repulsive hydration forces. We find, however, that the bound water molecules retain a shear fluidity characteristic of the bulk liquid, even when compressed down to films 1.0 +/- 0.3 nanometer thick. We attribute this to the ready exchange (as opposed to loss) of water molecules within the hydration layers as they rub past each other under strong compression.