Biophysical Determinants for the Viscosity of Concentrated Monoclonal Antibody Solutions.

Monoclonal antibodies (mAbs) are particularly relevant for therapeutics due to their high specificity and versatility, and mAb-based drugs are hence used to treat numerous diseases. The increased patient compliance of self-administration motivates the formulation of products for subcutaneous (SC) administration. The associated challenge is to formulate highly concentrated antibody solutions to achieve a significant therapeutic effect, while limiting their viscosity and preserving their physicochemical stability. Protein-protein interactions (PPIs) are in fact the root cause of several potential problems concerning the stability, manufacturability, and delivery of a drug product. The understanding of macroscopic viscosity requires an in-depth knowledge on protein diffusion, PPIs, and self-association/aggregation. Here, we study the self-diffusion of different mAbs of the IgG1 subtype in aqueous solution as a function of the concentration and temperature by quasi-elastic neutron scattering (QENS). QENS allows us to probe the short-time self-diffusion of the molecules and therefore to determine the hydrodynamic mAb cluster size and to gain information on the internal mAb dynamics. Small-angle neutron scattering (SANS) is jointly employed to probe structural details and to understand the nature and intensity of PPIs. Complementary information is provided by molecular dynamics (MD) simulations and viscometry, thus obtaining a comprehensive picture of mAb diffusion.

Effects of flexibility in coarse-grained models for bovine serum albumin and immunoglobulin G.

We construct a coarse-grained, structure-based, low-resolution, 6-bead flexible model of bovine serum albumin (BSA, PDB: 4F5S), which is a popular example of a globular protein in biophysical research. The model is obtained via direct Boltzmann inversion using all-atom simulations of a single molecule, and its particular form is selected from a large pool of 6-bead coarse-grained models using two suitable metrics that quantify the agreement in the distribution of collective coordinates between all-atom and coarse-grained Brownian dynamics simulations of solutions in the dilute limit. For immunoglobulin G (IgG), a similar structure-based 12-bead model has been introduced in the literature [Chaudhri et al., J. Phys. Chem. B 116, 8045 (2012)] and is employed here to compare findings for the compact BSA molecule and the more anisotropic IgG molecule. We define several modified coarse-grained models of BSA and IgG, which differ in their internal constraints and thus account for a variation of flexibility. We study denser solutions of the coarse-grained models with purely repulsive molecules (achievable by suitable salt conditions) and address the effect of packing and flexibility on dynamic and static behavior. Translational and rotational self-diffusivity is enhanced for more elastic models. Finally, we discuss a number of effective sphere sizes for the BSA molecule, which can be defined from its static and dynamic properties. Here, it is found that the effective sphere diameters lie between 4.9 and 6.1 nm, corresponding to a relative spread of about ±10% around a mean of 5.5 nm.

Molecular structural dynamics in water-ethanol mixtures: Spectroscopy with polarized neutrons simultaneously accessing collective and self-diffusion.

Binary mixtures of water with lower alcohols display non-linear phase behaviors upon mixing, which are attributed to potential cluster formation at the molecular level. Unravelling such elusive structures requires investigation of hydrogen-bonding sub-nanosecond dynamics. We employ high-resolution neutron time-of-flight spectroscopy with polarization analysis in combination with selective deuteration to study the concentration-dependent structural dynamics in the water rich part of the phase diagram of water-ethanol mixtures. This method enables simultaneous access to atomic correlations in space and time and allows us to separate spatially incoherent scattering probing self-diffusion of the ethanol fraction from the coherent scattering probing collective diffusion of the water network as a whole. Our observations indicate an enhanced rigidity of the hydrogen bond network at the mesoscopic length scale compared to the molecular scale as the ethanol fraction increases, which is consistent with the hypothesis of clusters.

Signature of functional enzyme dynamics in quasielastic neutron scattering spectra: The case of phosphoglycerate kinase.

We present an analysis of high-resolution quasi-elastic neutron scattering spectra of phosphoglycerate kinase which elucidates the influence of the enzymatic activity on the dynamics of the protein. We show that in the active state the inter-domain motions are amplified and the intra-domain asymptotic power-law relaxation ∝t-α is accelerated, with a reduced coefficient α. Employing an energy landscape picture of protein dynamics, this observation can be translated into a widening of the distribution of energy barriers separating conformational substates of the protein.

Strikingly Different Roles of SARS-CoV-2 Fusion Peptides Uncovered by Neutron Scattering.

Coronavirus disease-2019 (COVID-19), a potentially lethal respiratory illness caused by the coronavirus SARS-CoV-2, emerged in the end of 2019 and has since spread aggressively across the globe. A thorough understanding of the molecular mechanisms of cellular infection by coronaviruses is therefore of utmost importance. A critical stage in infection is the fusion between viral and host membranes. Here, we present a detailed investigation of the role of selected SARS-CoV-2 Spike fusion peptides, and the influence of calcium and cholesterol, in this fusion process. Structural information from specular neutron reflectometry and small angle neutron scattering, complemented by dynamics information from quasi-elastic and spin-echo neutron spectroscopy, revealed strikingly different functions encoded in the Spike fusion domain. Calcium drives the N-terminal of the Spike fusion domain to fully cross the host plasma membrane. Removing calcium, however, reorients the peptide back to the lipid leaflet closest to the virus, leading to significant changes in lipid fluidity and rigidity. In conjunction with other regions of the fusion domain, which are also positioned to bridge and dehydrate viral and host membranes, the molecular events leading to cell entry by SARS-CoV-2 are proposed.

Structure and diffusive dynamics of aspartate α-decarboxylase (ADC) liganded with D-serine in aqueous solution.

Incoherent neutron spectroscopy, in combination with dynamic light scattering, was used to investigate the effect of ligand binding on the center-of-mass self-diffusion and internal diffusive dynamics of Escherichia coli aspartate α-decarboxylase (ADC). The X-ray crystal structure of ADC in complex with the D-serine inhibitor was also determined, and molecular dynamics simulations were used to further probe the structural rearrangements that occur as a result of ligand binding. These experiments reveal that D-serine forms hydrogen bonds with some of the active site residues, that higher order oligomers of the ADC tetramer exist on ns-ms time-scales, and also show that ligand binding both affects the ADC internal diffusive dynamics and appears to further increase the size of the higher order oligomers.

Role of fast dynamics in the complexation of G-quadruplexes with small molecules.

G-quadruplexes (G4s) formed by the human telomeric sequence AG3 (TTAG3)3 (Tel22) play a key role in cancer and aging. We combined elastic incoherent neutron scattering (EINS) and quasielastic incoherent neutron scattering (QENS) to characterize the internal dynamics of Tel22 G4s and to assess how it is affected by complexation with two standard ligands, Berberine and BRACO19. We show that the interaction with the two ligands induces an increase of the overall mobility of Tel22 as quantified by the mean squared displacements (MSD) of hydrogen atoms. At the same time, the complexes display a lower stiffness than G4 alone. Two different types of motion characterize the G4 nanosecond timescale dynamics. Upon complexation, an increasing fraction of G4 atomic groups participate in this fast dynamics, along with an increase in the relevant characteristic length scales. We suggest that the entropic contribution to the conformational free energy of these motions might be crucial for the complexation mechanisms.

Multiscale relaxation dynamics and diffusion of myelin basic protein in solution studied by quasielastic neutron scattering.

We report an analysis of high-resolution quasielastic neutron scattering spectra from Myelin Basic Protein (MBP) in solution, comparing the spectra at three different temperatures (283, 303, and 323 K) for a pure D2O buffer and a mixture of D2O buffer with 30% of deuterated trifluoroethanol (TFE). Accompanying experiments with dynamic light scattering and Circular Dichroism (CD) spectroscopy have been performed to obtain, respectively, the global diffusion constant and the secondary structure content of the molecule for both buffers as a function of temperature. Modeling the decay of the neutron intermediate scattering function by the Mittag-Leffler relaxation function, ϕ(t) = Eα(-(t/τ)α) (0 < α < 1), we find that trifluoroethanol slows down the relaxation dynamics of the protein at 283 K and leads to a broader relaxation rate spectrum. This effect vanishes with increasing temperature, and at 323 K, its relaxation dynamics is identical in both solvents. These results are coherent with the data from dynamic light scattering, which show that the hydrodynamic radius of MBP in TFE-enriched solutions does not depend on temperature and is only slightly smaller compared to the pure D2O buffer, except for 283 K, where it is much reduced. In accordance with these observations, the CD spectra reveal that TFE induces essentially a partial transition from ß-strands to α-helices, but only a weak increase in the total secondary structure content, leaving about 50% of the protein unfolded. The results show that MBP is for all temperatures and in both buffers an intrinsically disordered protein and that TFE essentially induces a reduction in its hydrodynamic radius and its relaxation dynamics at low temperatures.

Zinc determines dynamical properties and aggregation kinetics of human insulin.

Protein aggregation is a widespread process leading to deleterious consequences in the organism, with amyloid aggregates being important not only in biology but also for drug design and biomaterial production. Insulin is a protein largely used in diabetes treatment, and its amyloid aggregation is at the basis of the so-called insulin-derived amyloidosis. Here, we uncover the major role of zinc in both insulin dynamics and aggregation kinetics at low pH, in which the formation of different amyloid superstructures (fibrils and spherulites) can be thermally induced. Amyloid aggregation is accompanied by zinc release and the suppression of water-sustained insulin dynamics, as shown by particle-induced x-ray emission and x-ray absorption spectroscopy and by neutron spectroscopy, respectively. Our study shows that zinc binding stabilizes the native form of insulin by facilitating hydration of this hydrophobic protein and suggests that introducing new binding sites for zinc can improve insulin stability and tune its aggregation propensity.

Molecular Flexibility of Antibodies Preserved Even in the Dense Phase after Macroscopic Phase Separation.

Antibody therapies are typically based on high-concentration formulations that need to be administered subcutaneously. These conditions induce several challenges, inter alia a viscosity suitable for injection, sufficient solution stability, and preservation of molecular function. To obtain systematic insights into the molecular factors, we study the dynamics on the molecular level under strongly varying solution conditions. In particular, we use solutions of antibodies with poly(ethylene glycol), in which simple cooling from room temperature to freezing temperatures induces a transition from a well-dispersed solution into a phase-separated and macroscopically arrested system. Using quasi-elastic neutron scattering during in situ cooling ramps and in prethermalized measurements, we observe a strong decrease in antibody diffusion, while internal flexibility persists to a significant degree, thus ensuring the movement necessary for the preservation of molecular function. These results are relevant for a more dynamic understanding of antibodies in high-concentration formulations, which affects the formation of transient clusters governing the solution viscosity.

Temperature and salt controlled tuning of protein clusters.

The formation of molecular assemblies in protein solutions is of strong interest both from a fundamental viewpoint and for biomedical applications. While ordered and desired protein assemblies are indispensable for some biological functions, undesired protein condensation can induce serious diseases. As a common cofactor, the presence of salt ions is essential for some biological processes involving proteins, and in aqueous suspensions of proteins can also give rise to complex phase diagrams including homogeneous solutions, large aggregates, and dissolution regimes. Here, we systematically study the cluster formation approaching the phase separation in aqueous solutions of the globular protein BSA as a function of temperature (T), the protein concentration (cp) and the concentrations of the trivalent salts YCl3 and LaCl3 (cs). As an important complement to structural, i.e. time-averaged, techniques we employ a dynamical technique that can detect clusters even when they are transient on the order of a few nanoseconds. By employing incoherent neutron spectroscopy, we unambiguously determine the short-time self-diffusion of the protein clusters depending on cp, cs and T. We determine the cluster size in terms of effective hydrodynamic radii as manifested by the cluster center-of-mass diffusion coefficients D. For both salts, we find a simple functional form D(cp, cs, T) in the parameter range explored. The calculated inter-particle attraction strength, determined from the microscopic and short-time diffusive properties of the samples, increases with salt concentration and temperature in the regime investigated and can be linked to the macroscopic behavior of the samples.

The modifying effect of supramolecular gel fibres on the diffusion of paracetamol and ibuprofen sodium on the picosecond timescale.

Employing neutron spectroscopy, we follow the tracer diffusion of two non-steroidal anti-inflammatory drug molecules, paracetamol (PCM) and ibuprofen sodium (IBU), in a supramolecular gel and the corresponding bulk solution. Both solutes show altered diffusion behaviour in the gel phase, deviating from each other and their bulk solution. Whilst picosecond diffusion of IBU is slightly quicker in the gel, this effect is significantly increased for PCM, which is up to 70% quicker in the gel than in solution. This effect is independent of changes in the solvent diffusion reported previously. An increased residence time of PCM in solution at lower temperatures points towards the onset of nucleation and crystallisation. This work reports one of the first experiments on the novel Backscattering and Time-of-Flight option (BATS) on the IN16B spectrometer at the Institut Laue-Langevin, France, which with its range and resolution in neutron energy and momentum transfer is ideally suited to observe this type of diffusion.

Evolution of the structure and dynamics of bovine serum albumin induced by thermal denaturation.

Protein denaturation in concentrated solutions consists of the unfolding of the native protein structure, and subsequent cross-linking into clusters or gel networks. While the kinetic evolution of structure has been studied for some cases, the underlying microscopic dynamics of proteins has so far been neglected. However, protein dynamics is essential to understand the specific nature of assembly processes, such as diffusion-limited growth, or vitrification of dense liquids. Here, we present a study on thermal denaturation of concentrated solutions of bovine serum albumin (BSA) in D2O with and without NaCl. Using small-angle scattering, we provide information on structure before, during and after denaturation. Using quasi-elastic neutron scattering, we monitor in real-time the microscopic dynamics and dynamical confinement throughout the entire denaturation process covering protein unfolding and cross-linking. After denaturation, the protein dynamics is slowed down in salty solutions compared to those in pure water, while the stability and dynamics of the native solution appears unaffected by salt. The approach presented here opens opportunities to link microscopic dynamics to emerging structural properties, with implications for assembly processes in soft and biological matter.

Picosecond self-diffusion in ethanol-water mixtures.

We report the self-diffusion in ethanol-water mixtures as a function of the water-ethanol ratio measured at different temperatures using quasi-elastic neutron spectroscopy (QENS). For our protiated samples, QENS is mainly sensitive to the dominant ensemble-averaged incoherent scattering from the hydrogen atoms of the liquid mixtures. The energy range and resolution render our experiment sensitive to the picosecond time scale and nanometer length scale. These observation scales complement different scales accessible by nuclear magnetic resonance techniques. Subsequent to testing different models, we find that a simple jump-diffusion model averaging over both types of molecules, water and ethanol, best fits our data.

Osmolytes modify protein dynamics and function of tetrameric lactate dehydrogenase upon pressurization.

We present a study of the combined effects of natural cosolvents (TMAO, glycine, urea) and pressure on the activity of the tetrameric enzyme lactate dehydrogenase (LDH). To this end, high-pressure stopped-flow methodology in concert with fast UV/Vis spectroscopic detection of product formation was applied. To reveal possible pressure effects on the stability and dynamics of the enzyme, FTIR spectroscopic and neutron scattering measurements were carried out. In neat buffer solution, the catalytic turnover number of the enzyme, kcat, increases up to 1000 bar, the pressure range where dissociation of the tetrameric species to dimers sets in. Accordingly, we obtain a negative activation volume, ΔV# = -45.3 mL mol-1. Further, the enzyme substrate complex has a larger volume compared to the enzyme and substrate in the unbound state. The neutron scattering data show that changes in the fast internal dynamics of the enzyme are not responsible for the increase of kcat upon compression. Whereas the magnitude of kcat is similar in the presence of the osmolytes, the pressure of deactivation is modulated by the addition of cosolvents. TMAO and glycine increase the pressure of deactivation, and in accordance with the observed stabilizing effect both cosolvents exhibit against denaturation and/or dissociation of proteins. While urea does not markedly affect the magnitude of the Michaelis constant, KM, both 1 M TMAO and 1 M glycine exhibit smaller KM values of about 0.07 mM and 0.05 mM below about 1 kbar. Such positive effect on the substrate affinity could be rationalized by the effect the two cosolutes impose on the thermodynamic activities of the reactants, which reflect changes in water-mediated intermolecular interactions. Our data show that the intracellular milieu, i.e., the solution conditions that have evolved, may be sufficient to maintain enzymatic activity under extreme environmental conditions, including the whole pressure range encountered on Earth.

Effect of Phosphorylation on a Human-like Osteopontin Peptide.

The last decade established that the dynamic properties of the phosphoproteome are central to function and its modulation. The temporal dimension of phosphorylation effects remains nonetheless poorly understood, particularly for intrinsically disordered proteins. Osteopontin, selected for this study due to its key role in biomineralization, is expressed in many species and tissues to play a range of distinct roles. A notable property of highly phosphorylated isoforms of osteopontin is their ability to sequester nanoclusters of calcium phosphate to form a core-shell structure, in a fluid that is supersaturated but stable. In Biology, this process enables soft and hard tissues to coexist in the same organism with relative ease. Here, we extend our understanding of the effect of phosphorylation on a disordered protein, the recombinant human-like osteopontin rOPN. The solution structures of the phosphorylated and unphosphorylated rOPN were investigated by small-angle x-ray scattering and no significant changes were detected on the radius of gyration or maximum interatomic distance. The picosecond-to-nanosecond dynamics of the hydrated powders of the two rOPN forms were further compared by elastic and quasi-elastic incoherent neutron scattering. Phosphorylation was found to block some nanosecond side-chain motions while increasing the flexibility of other side chains on the faster timescale. Phosphorylation can thus selectively change the dynamic behavior of even a highly disordered protein such as osteopontin. Through such an effect on rOPN, phosphorylation can direct allosteric mechanisms, interactions with substrates, cofactors and, in this case, amorphous or crystalline biominerals.

Alzheimer's peptide amyloid-ß, fragment 22-40, perturbs lipid dynamics.

The peptide amyloid-ß (Aß) interacts with membranes of cells in the human brain and is associated with Alzheimer's disease (AD). The intercalation of Aß in membranes alters membrane properties, including the structure and lipid dynamics. Any change in the membrane lipid dynamics will affect essential membrane processes, such as energy conversion, signal transduction and amyloid precursor protein (APP) processing, and may result in the observed neurotoxicity associated with the disease. The influence of this peptide on membrane dynamics was studied with quasi-elastic neutron scattering, a technique which allows a wide range of observation times from picoseconds to nanoseconds, over nanometer length scales. The effect of the membrane integral neurotoxic peptide amyloid-ß, residues 22-40, on the in- and out-of-plane lipid dynamics was observed in an oriented DMPC/DMPS bilayer at 15 °C, in its gel phase, and at 30 °C, near the phase transition temperature of the lipids. Near the phase-transition temperature, a 1.5 mol% of peptide causes up to a twofold decrease in the lipid diffusion coefficients. In the gel-phase, this effect is reversed, with amyloid-ß(22-40) increasing the lipid diffusion coefficients. The observed changes in lipid diffusion are relevant to protein-protein interactions, which are strongly influenced by the diffusion of membrane components. The effect of the amyloid-ß peptide fragment on the diffusion of membrane lipids will provide insight into the membrane's role in AD.

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. Here, we report an experimental study on the conformational dynamics of different folded conformational states of apo- and holomyoglobin in solution. Global protein diffusion and internal molecular motions were probed by neutron time-of-flight and neutron backscattering spectroscopy on the picosecond and nanosecond time scales. Global protein diffusion was found to depend on the α-helical content of the protein suggesting that charges on the macromolecule increase the short-time diffusion of protein. With regard to the molten globules, a gel-like phase due to protein entanglement and interactions with neighbouring macromolecules was visible due to a reduction of the global diffusion coefficients on the nanosecond time scale. Diffusion coefficients, residence and relaxation times of internal protein dynamics and root mean square displacements of localised internal motions were determined for the investigated structural states. The difference in conformational entropy ΔSconf of the protein between the unfolded and the partially or fully folded conformations was extracted from the measured root mean square displacements. Using thermodynamic parameters from the literature and the experimentally determined ΔSconf values we could identify the entropic contribution of the hydration shell ΔShydr of the different folded states. Our results point out the relevance of conformational entropy of the protein and the hydration shell for stability and folding of myoglobin.

Anomalous and anisotropic nanoscale diffusion of hydration water molecules in fluid lipid membranes.

We have studied nanoscale diffusion of membrane hydration water in fluid-phase lipid bilayers made of 1,2-dimyristoyl-3-phosphocholine (DMPC) using incoherent quasi-elastic neutron scattering. Dynamics were fit directly in the energy domain using the Fourier transform of a stretched exponential. By using large, 2-dimensional detectors, lateral motions of water molecules and motions perpendicular to the membranes could be studied simultaneously, resulting in 2-dimensional maps of relaxation time, τ, and stretching exponent, ß. We present experimental evidence for anomalous (sub-diffusive) and anisotropic diffusion of membrane hydration water molecules over nanometer distances. By combining molecular dynamics and Brownian dynamics simulations, the potential microscopic origins for the anomaly and anisotropy of hydration water were investigated. Bulk water was found to show intrinsic sub-diffusive motion at time scales of several picoseconds, likely related to caging effects. In membrane hydration water, however, the anisotropy of confinement and local dynamical environments leads to an anisotropy of relaxation times and stretched exponents, indicative of anomalous dynamics.

Hierarchical molecular dynamics of bovine serum albumin in concentrated aqueous solution below and above thermal denaturation.

The dynamics of proteins in solution is a complex and hierarchical process, affected by the aqueous environment as well as temperature. We present a comprehensive study on nanosecond time and nanometer length scales below, at, and above the denaturation temperature Td. Our experimental data evidence dynamical processes in protein solutions on three distinct time scales. We suggest a consistent physical picture of hierarchical protein dynamics: (i) self-diffusion of the entire protein molecule is confirmed to agree with colloid theory for all temperatures where the protein is in its native conformational state. At higher temperatures T > Td, the self-diffusion is strongly obstructed by cross-linking or entanglement. (ii) The amplitude of backbone fluctuations grows with increasing T, and a transition in its dynamics is observed above Td. (iii) The number of mobile side-chains increases sharply at Td while their average dynamics exhibits only little variations. The combination of quasi-elastic neutron scattering and the presented analytical framework provides a detailed microscopic picture of the protein molecular dynamics in solution, thereby reflecting the changes of macroscopic properties such as cluster formation and gelation.