Measuring macromolecular size distributions and interactions at high concentrations by sedimentation velocity

A F420-dependent Single Domain Chemogenetic Tool for Protein De-dimerization

Protein-protein interactions (PPIs) mediate many fundamental cellular processes. Control of PPIs through optically or chemically responsive protein domains has had a profound impact on basic research and some clinical applications. Most chemogenetic methods induce the association, i.e., dimerization or oligomerization, of target proteins, whilst the few available dissociation approaches either break large oligomeric protein clusters or heteromeric complexes. Here, we have exploited the controlled dissociation of a homodimeric oxidoreductase from mycobacteria (MSMEG_2027) by its native cofactor, F420, which is not present in mammals, as a bioorthogonal monomerization switch. Using X-ray crystallography, we found that in the absence of F420 MSMEG_2027 forms a unique domain-swapped dimer that occludes the cofactor binding site. Rearrangement of the N-terminal helix upon F420 binding results in the dissolution of the dimer. We then showed that MSMEG_2027 can be fused to proteins of interest in human cells and applied it as a tool to induce and release MAPK/ERK signalling downstream of a chimeric fibroblast growth factor receptor 1 (FGFR1) tyrosine kinase. This F420-dependent chemogenetic de-homodimerization tool is stoichiometric and based on a single domain and thus represents a novel mechanism to investigate protein complexes in situ.

Analytical ultracentrifugation as a tool for exploring COSAN assemblies

The self-assembly of the cobaltabis(dicarbollide) (COSAN) anionic boron clusters into micelles above a critical micelle concentration (cmc) of 10-20 mM and its behavior as "sticky nano-ions" facilitating controlled protein aggregation have been previously investigated using scattering techniques. These techniques effectively provide average structural parameters but, when applied to colloidal systems, often rely on models assuming polydispersity or anisotropic shapes. Here, we employed sedimentation velocity analytical ultracentrifugation (SV-AUC), which offers the ability to resolve discrete species. We revisited two key questions: (1) the aggregation behavior of COSAN into micelles, a topic still under debate, and (2) the nature of the protein assemblies induced by COSAN, specifically their size/shape distribution and aggregation number. SV-AUC confirms the cmc of COSAN of 16 mM and reveals that COSAN micelles exhibit low aggregation numbers (8 in water and 14 in dilute salt), consistent with recent hypotheses. It shows that COSAN promotes myoglobin aggregation into discrete oligomeric species with well-defined aggregation numbers, such as dimers, tetramers, and higher-order assemblies, depending on the COSAN-to-protein ratio. COSAN binding could be quantified at the lower COSAN/myoglobin ratios. For example, at ratio 5, myoglobin monomer (25%) binds about two COSANs, dimer (45%) about 14 COSANs, and there are ≈ 30% very large aggregates. These results provide clarity on the discrete nature of COSAN micelle aggregation and protein assembly. This study highlights the complementary role of SV-AUC in understanding supramolecular assemblies, offering useful insights into the behavior of COSAN nano-ions and their interactions with biomacromolecules.

Is there reversible dimerization of albumin in blood plasma? And does it matter?

Most albumin in blood plasma is thought to be monomeric with some 5% covalently dimerized. However, many reports in the recent biophysics literature find that albumin is reversibly dimerized or even oligomerized. We review data on this from X-ray crystallography and diverse biophysical techniques. The number-average molecular weight of albumin would be increased by dimerization, affecting size-dependent filtration processes of albumin such as at the glycocalyx of the capillary endothelium and the podocyte slit-diaphragm of the renal glomerulus. If correct, and depending on characteristics of the process, such as Kd, reversible dimerization of albumin in plasma would have major implications for normal physiology and medicine. We present quantitative models of the impact of dimerization on albumin molecular forms, on the number-average molecular weight of albumin, and estimate the effect on the colloid osmotic pressure of albumin. Dimerization reduces colloid osmotic pressure as total albumin concentration increases below that expected in the absence of dimerization. Current models of albumin filtration by the renal glomerulus would need revision to account for the dynamic size of albumin molecules filtered. More robust biophysical data are needed to give a definitive answer to the questions posed and we suggest possible approaches to this.

Challenges in Peptide Solubilization - Amyloids Case Study

Peptide science has been a rapidly growing research field because of the enormous potential application of these biocompatible and bioactive molecules. However, many factors limit the widespread use of peptides in medicine, and low solubility is among the most common problems that hamper drug development in the early stages of research. Solubility is a crucial, albeit poorly understood, feature that determines peptide behavior. Several different solubility predictors have been proposed, and many strategies and protocols have been reported to dissolve peptides, but none of them is a one-size-fits-all method for solubilization of even the same peptide. In this review, we look for the reasons behind the difficulties in dissolving peptides, analyze the factors influencing peptide aggregation, conduct a critical analysis of solubilization strategies and protocols available in the literature, and give some tips on how to deal with the so-called difficult sequences. We focus on amyloids, which are particularly difficult to dissolve and handle such as amyloid beta (Aβ), insulin, and phenol-soluble modulins (PSMs).

BASIS: BioAnalysis SEDFIT integrated software for cGMP analysis of SV-AUC data

Sedimentation velocity analytical ultracentrifugation (SV-AUC) has long been an important method for characterization of antibody therapeutics. Recently, SV-AUC has experienced a wave of new interest and usage from the gene and cell therapy industry, where SV-AUC has proven itself to be the "gold standard" analytical approach for determining capsid loading ratios for adeno-associated virus (AAV) and other viral vectors. While other more common approaches have existed in the realm of cGMP-compliant techniques for years, SV-AUC has long been used strictly for characterization, but not for release testing. This manuscript describes the challenges faced in bringing SV-AUC to a cGMP environment and describes a new program, "BASIS", which allows for 21 CFR Part 11-compliant data handling and data analysis using the well-known and frequently cited SEDFIT analysis software.

An automated interface for sedimentation velocity analysis in SEDFIT

Sedimentation velocity analytical ultracentrifugation (SV-AUC) is an indispensable tool for the study of particle size distributions in biopharmaceutical industry, for example, to characterize protein therapeutics and vaccine products. In particular, the diffusion-deconvoluted sedimentation coefficient distribution analysis, in the software SEDFIT, has found widespread applications due to its relatively high resolution and sensitivity. However, a lack of suitable software compatible with Good Manufacturing Practices (GMP) has hampered the use of SV-AUC in this regulatory environment. To address this, we have created an interface for SEDFIT so that it can serve as an automatically spawned module with controlled data input through command line parameters and output of key results in files. The interface can be integrated in custom GMP compatible software, and in scripts that provide documentation and meta-analyses for replicate or related samples, for example, to streamline analysis of large families of experimental data, such as binding isotherm analyses in the study of protein interactions. To test and demonstrate this approach we provide a MATLAB script mlSEDFIT.

Retrospective rationalization of disparities between the concentration dependence of diffusion coefficients obtained by boundary spreading and dynamic light scattering

This study establishes the existence of substantial agreement between published results from traditional boundary spreading measurements (including synthetic boundary measurements in the analytical ultracenrifuge) on two globular proteins (bovine serum albumin, ovalbumin) and the concentration dependence of diffusion coefficient predicted for experiments conducted under the operative thermodynamic constraints of constant temperature and solvent chemical potential. Although slight negative concentration dependence of the translational diffusion coefficient is the experimentally observed as well as theoretically predicted, the extent of the concentration dependence is within the limits of experimental uncertainty inherent in diffusion coefficient measurement. Attention is then directed toward the ionic strength dependence of the concentration dependence coefficient ([Formula: see text]) describing diffusion coefficients obtained by dynamic light scattering, where, in principle, the operative thermodynamic constraints of constant temperature and pressure preclude consideration of results in terms of single-solute theory. Nevertheless, good agreement between predicted and published experimental ionic strength dependencies of [Formula: see text] for lysozyme and an immunoglobulin is observed by a minor adaptation of the theoretical treatment to accommodate the fact that thermodynamic activity is monitored on the molal concentration scale because of the constraint of constant pressure that pertains in dynamic light scattering experiments.

Simulation of Gilbert theory for self-association in sedimentation velocity experiments: a guide to evaluate best fitting models

There is a long tradition in the Biophysics community of using simulations as a means to understand macromolecular behavior in various physicochemical methods. This allows a rigorous means to interpret observations in terms of fundamental principles, including chemical equilibrium, reaction kinetics, transport processes and thermodynamics. Here we simulate data for the Gilbert Theory for self-association, a fundamental analytical ultracentrifuge (AUC) technique to understand the shape of sedimentation velocity reaction boundaries that involve reversible monomer-Nmer interactions. Simulating monomer-dimer through monomer-hexamer systems as a function of concentration about the equilibrium constant allows a visual means to differentiate reaction stoichiometry by determining end points and inflection positions. Including intermediates (eg A1-A2-A3-A4-A5-A6) in the simulations reveals the smoothing of the reaction boundary and the removal of sharp inflections between monomers and polymers. The addition of cooperativity restores sharp boundaries or peaks to the observation and allows more discrimination in the selection of possible fitting models. Thermodynamic nonideality adds additional features when applied across wide ranges of concentration that might be appropriate for high-concentration therapeutic monoclonal antibody (mAb) solutions. This presentation serves as a tutorial for using modern AUC analysis software like SEDANAL for selecting potential fitting models.

An automated interface for sedimentation velocity analysis in SEDFIT

Sedimentation velocity analytical ultracentrifugation (SV-AUC) is an indispensable tool for the study of particle size distributions in biopharmaceutical industry, for example, to characterize protein therapeutics and vaccine products. In particular, the diffusion-deconvoluted sedimentation coefficient distribution analysis, in the software SEDFIT, has found widespread applications due to its relatively high resolution and sensitivity. However, a lack of available software compatible with Good Manufacturing Practices (GMP) has hampered the use of SV-AUC in this regulatory environment. To address this, we have created an interface for SEDFIT so that it can serve as an automatically spawned module with controlled data input through command line parameters and output of key results in files. The interface can be integrated in custom GMP compatible software, and in scripts that provide documentation and meta-analyses for replicate or related samples, for example, to streamline analysis of large families of experimental data, such as binding isotherm analyses in the study of protein interactions. To test and demonstrate this approach we provide a MATLAB script mlSEDFIT.

Insights into the biochemical and biophysical mechanisms mediating the longevity of the transparent optics of the eye lens

In the human eye, a transparent cornea and lens combine to form the "refracton" to focus images on the retina. This requires the refracton to have a high refractive index "n," mediated largely by extracellular collagen fibrils in the corneal stroma and the highly concentrated crystallin proteins in the cytoplasm of the lens fiber cells. Transparency is a result of short-range order in the spatial arrangement of corneal collagen fibrils and lens crystallins, generated in part by post-translational modifications (PTMs). However, while corneal collagen is remodeled continuously and replaced, lens crystallins are very long-lived and are not replaced and so accumulate PTMs over a lifetime. Eventually, a tipping point is reached when protein aggregation results in increased light scatter, inevitably leading to the iconic protein condensation-based disease, age-related cataract (ARC). Cataracts account for 50% of vision impairment worldwide, affecting far more people than other well-known protein aggregation-based diseases. However, because accumulation of crystallin PTMs begins before birth and long before ARC presents, we postulate that the lens protein PTMs contribute to a "cataractogenic load" that not only increases with age but also has protective effects on optical function by stabilizing lens crystallins until a tipping point is reached. In this review, we highlight decades of experimental findings that support the potential for PTMs to be protective during normal development. We hypothesize that ARC is preventable by protecting the biochemical and biophysical properties of lens proteins needed to maintain transparency, refraction, and optical function.

Determination of specific and non-specific protein-protein interactions for beta-lactoglobulin by analytical ultracentrifugation and membrane osmometry experiments

Protein-protein interactions are essential for the understanding of biological processes. Specific protein aggregation is an important aspect for many biological systems. In particular, electrostatic interactions play the key role for protein-protein interactions, as many amino acids have pH-dependent charge states. Moreover, protein dissociation is directly related to the solution pH, ionic strength, temperature and protein concentration. The subtle interplay between different specific and non-specific interactions is demonstrated for beta-lactoglobulin (BLG) with a focus on low salt concentrations, thus mimicking technically relevant processing conditions. BLG is a well-characterized model system, proven to attain its monomer-dimer equilibrium strongly dependent upon the pH of the solution. In this manuscript, we present a unique combination of analytical ultracentrifugation and membrane osmometry experiments, which quantifies specific and non-specific interactions, i.e. in terms of the dimer dissociation constants and the second osmotic virial coefficient, at pH 3 and 7 and sodium chloride concentrations of 10 mM and 100 mM. This provides direct insight to protein-protein interactions for a system with a concentration-dependent monomer-dimer equilibrium. Moreover, using a coarse-grained extended DLVO model in combination with molecular dynamics simulations, we quantify non-specific monomer-monomer, monomer-dimer and dimer-dimer interactions as well as the binding free energy of BLG dimerization from theoretical calculations. The experimentally determined interactions are shown to be mainly governed by electrostatic interactions and further agree with free energy calculations. Our experimental protocol aims to determine non-specific and specific interactions for a dynamically interacting system and provides an understanding of protein-protein interactions for BLG at low salt concentrations.

Best Practices for Aggregate Quantitation of Antibody Therapeutics by Sedimentation Velocity Analytical Ultracentrifugation

Analytical ultracentrifugation (AUC) is a critical analytical tool supporting the development and manufacture of protein therapeutics. AUC is routinely used as an assay orthogonal to size exclusion chromatography for aggregate quantitation. This article distills the experimental and analysis procedures used by the authors for sedimentation velocity AUC into a series of best-practices considerations. The goal of this distillation is to help harmonize aggregate quantitation approaches across the biopharmaceutical industry. We review key considerations for sample and instrument suitability, experimental design, and data analysis best practices and conversely, highlight potential pitfalls to accurate aggregate analysis. Our goal is to provide experienced users benchmarks against which they can standardize their analyses and to provide guidance for new AUC analysts that will aid them to become proficient in this fundamental technique.

Label-free methods for optical in vitro characterization of protein-protein interactions

Protein-protein interactions are involved in the regulation and function of the majority of cellular processes. As a result, much effort has been aimed at the development of methodologies capable of quantifying protein-protein interactions, with label-free methods being of particular interest due to the associated simplified workflows and minimisation of label-induced perturbations. Here, we review recent advances in optical technologies providing label-free in vitro measurements of affinities and kinetics. We provide an overview and comparison of existing techniques and their principles, discussing advantages, limitations, and recent applications.

Probing interactions of therapeutic antibodies with serum via second virial coefficient measurements

Antibody-based therapeutics are the fastest-growing drug class on the market, used to treat aggressive forms of cancer, chronic autoimmune conditions, and numerous other disease states. Although the specificity, affinity, and versatility of therapeutic antibodies can provide an advantage over traditional small-molecule drugs, their development and optimization can be much more challenging and time-consuming. This is, in part, because the ideal formulation buffer systems used for in vitro characterization inadequately reflect the crowded biological environments (serum, endosomal lumen, etc.) that these drugs experience once administered to a patient. Such environments can perturb the binding of antibodies to their antigens and receptors, as well as homo- and hetero-aggregation, thereby altering therapeutic effect and disposition in ways that are incompletely understood. Although excluded volume effects are classically thought to favor binding, weak interactions with co-solutes in crowded conditions can inhibit binding. The second virial coefficient (B₂) parameter quantifies such weak interactions and can be determined by a variety of techniques in dilute solution, but analogous methods in complex biological fluids are not well established. Here, we demonstrate that fluorescence correlation spectroscopy is able to measure diffusive B₂-values directly in undiluted serum. Apparent second virial coefficient (B_2,app) measurements of antibodies in serum reveal that changes in the balance between attractive and repulsive interactions can dramatically impact global nonideality. Furthermore, our findings suggest that the approach of isolating specific components and completing independent cross-term virial coefficient measurements may not be an effective approach to characterizing nonideality in serum. The approach presented here could enrich our understanding of the effects of biological environments on proteins in general and advance the development of therapeutic antibodies and other protein-based therapeutics.

A multi-step nucleation process determines the kinetics of prion-like domain phase separation

Compartmentalization by liquid-liquid phase separation (LLPS) has emerged as a ubiquitous mechanism underlying the organization of biomolecules in space and time. Here, we combine rapid-mixing time-resolved small-angle X-ray scattering (SAXS) approaches to characterize the assembly kinetics of a prototypical prion-like domain with equilibrium techniques that characterize its phase boundaries and the size distribution of clusters prior to phase separation. We find two kinetic regimes on the micro- to millisecond timescale that are distinguished by the size distribution of clusters. At the nanoscale, small complexes are formed with low affinity. After initial unfavorable complex assembly, additional monomers are added with higher affinity. At the mesoscale, assembly resembles classical homogeneous nucleation. Careful multi-pronged characterization is required for the understanding of condensate assembly mechanisms and will promote understanding of how the kinetics of biological phase separation is encoded in biomolecules.

Boundary convection during velocity sedimentation in the Optima analytical ultracentrifuge

Analytical ultracentrifugation (AUC) provides the most widely applicable, precise, and accurate means for characterizing solution hydrodynamic and thermodynamic properties. While generally useful, boundary sedimentation velocity AUC (SV-AUC) analysis has become particularly important in assessing protein aggregation, fragmentation and conformational variants in the same solvents used during drug development and production. In early 2017 the only manufacturer of the analytical ultracentrifuge released its newest analytical ultracentrifuge, the Optima, to replace the aging second-generation XLA/I series ultracentrifuges. However, SV-AUC data from four Optima units used in the characterization of adeno-associated virus (AAV) have shown evidence of sample convection. Further investigation reveals this problem arises from the design of the temperature control system, which makes it prone to producing destabilizing temperature-induced density gradients that can lead to density inversions. The problem is intermittent and variable in severity within a given Optima unit and between Optima units. This convection appears to be associated mainly with low rotor speeds and dilute concentration of solvent components, i.e., AAV analysis conditions. Data features diagnostic for this problem and strategies for its elimination or minimization are provided.

Versatile Route toward Hydrophobically Polymer-Grafted Gold Nanoparticles from Aqueous Dispersions

Stabilization of gold nanoparticles in organic solvents is a key challenge in making them available for a wider range of material applications. Polymers are often used as stabilizing ligands because they also allow for the introduction of new properties and functionalities. Many of the established synthesis protocols for gold nanoparticles are water-based. However, the insolubility of many synthetic polymers in water renders the direct functionalization of aqueous particle dispersions with these ligands difficult. Here, we report on an approach for the functionalization of gold nanoparticles, which were prepared by aqueous synthesis, with hydrophobic polymer ligands and their characterization in nonpolar, organic dispersions. Our method employs an auxiliary ligand to first transfer gold nanoparticles from an aqueous to an organic medium. In the organic phase, the auxiliary ligand is then displaced by thiolated polystyrene ligands to form a dense polymer brush on the particle surface. We characterize the structure of the ligand shell using electron microscopy, scattering techniques, and ultracentrifugation and analyze the influence of the molecular weight of the polystyrene ligands on the structure of the polymer brush. We further investigate the colloidal stability of polystyrene-functionalized gold nanoparticles in various organic solvents. Finally, we extend the use of our protocol from small, spherical gold nanoparticles to larger gold nanorods and nanocubes.

Measuring aggregates, self-association, and weak interactions in concentrated therapeutic antibody solutions

Monoclonal antibodies are a class of biotherapeutics used for an increasing variety of disorders, including cancer, autoimmune, neurodegenerative, and viral diseases. Besides their antigen specificity, therapeutic use also mandates control of their solution interactions and colloidal properties in order to achieve a stable, efficacious, non-immunogenic, and low viscosity antibody solution at concentrations in the range of 50–150 mg/mL. This requires characterization of their reversible self-association, aggregation, and weak attractive and repulsive interactions governing macromolecular distance distributions in solution. Simultaneous measurement of these properties, however, has been hampered by solution nonideality. Based on a recently introduced sedimentation velocity method for measuring macromolecular size distributions in a mean-field approximation for hydrodynamic interactions, we demonstrate simultaneous measurement of polydispersity and weak and strong solution interactions in a panel of antibodies with concentrations up to 45 mg/mL. By allowing approximately an order of magnitude higher concentrations than previously possible in sedimentation velocity size distribution analysis, this approach can substantially improve efficiency and sensitivity for characterizing polydispersity and interactions of therapeutic antibodies at or close to formulation conditions.

Modification of cell wall polysaccharide guides cell division in Streptococcus mutans

In ovoid-shaped, Gram-positive bacteria, MapZ guides FtsZ-ring positioning at cell equators. The cell wall of the ovococcus Streptococcus mutans contains peptidoglycan decorated with serotype c carbohydrates (SCCs). In the present study, we identify the major cell separation autolysin AtlA as an SCC-binding protein. AtlA binding to SCC is attenuated by the glycerol phosphate (GroP) modification. Using fluorescently labeled AtlA constructs, we mapped SCC distribution on the streptococcal surface, revealing enrichment of GroP-deficient immature SCCs at the cell poles and equators. The immature SCCs co-localize with MapZ at the equatorial rings throughout the cell cycle. In GroP-deficient mutants, AtlA is mislocalized, resulting in dysregulated cellular autolysis. These mutants display morphological abnormalities associated with MapZ mislocalization, leading to FtsZ-ring misplacement. Altogether, our data support a model in which maturation of a cell wall polysaccharide provides the molecular cues for the recruitment of cell division machinery, ensuring proper daughter cell separation and FtsZ-ring positioning.

Quantifying the concentration dependence of sedimentation coefficients for globular macromolecules: a continuing age-old problem

This retrospective investigation has established that the early theoretical attempts to directly incorporate the consequences of radial dilution into expressions for variation of the sedimentation coefficient as a function of the loading concentration in sedimentation velocity experiments require concentration distributions exhibiting far greater precision than that achieved by the optical systems of past and current analytical ultracentrifuges. In terms of current methods of sedimentation coefficient measurement, until such improvement is made, the simplest procedure for quantifying linear s-c dependence (or linear concentration dependence of 1/s) for dilute systems therefore entails consideration of the sedimentation coefficient obtained by standard c(s), g*(s) or G(s) analysis) as an average parameter (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{s} $$\end{document}s¯) that pertains to the corresponding mean plateau concentration (following radial dilution) (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{c} $$\end{document}c¯) over the range of sedimentation velocity distributions used for the determination of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{s} $$\end{document}s¯. The relation of this with current descriptions of the concentration dependence of the sedimentation and translational diffusion coefficients is considered, together with a suggestion for the necessary improvement in the optical system.

Probing sedimentation non-ideality of particulate systems using analytical centrifugation

Analytical centrifugation is a versatile technique for the quantitative characterization of colloidal systems including colloidal stability. The recent developments in data acquisition and evaluation allow the accurate determination of particle size, shape anisotropy and particle density. High precision analytical centrifugation is in particular suited for the study of particle interactions and concentration-dependent sedimentation coefficients. We present a holistic approach for the quantitative determination of sedimentation non-ideality via analytical centrifugation for polydisperse, plain and amino-functionalized silica particles spanning over one order of magnitude in particle size between 100 nm and 1200 nm. These systems typically behave as neutral hard spheres as predicted by auxiliary lattice Boltzmann simulations. The extent of electrostatic interactions and their impact on sedimentation non-ideality can be quantified by the repulsion range, which is the ratio of the Debye length and the average interparticle distance. Experimental access to the repulsion range is provided through conductivity measurements. With the experimental repulsion range at hand, we estimate the effect of polydispersity on concentration-dependent sedimentation properties through a combination of lattice Boltzmann and Brownian dynamics simulations. Finally, we determine the concentration-dependent sedimentation properties of charge-stabilized, fluorescently-labeled silica particles with a nominal particle size of 30 nm and reduced interparticle distance, hence an elevated repulsion range. Overall, our results demonstrate how the influence of hard-sphere type and electrostatic interactions can be quantified when probing sedimentation non-ideality of particulate systems using analytical centrifugation even for systems exhibiting moderate sample heterogeneity and complex interactions.

Global multi-method analysis of interaction parameters for reversibly self-associating macromolecules at high concentrations

Weak macromolecular interactions assume a dominant role in the behavior of highly concentrated solutions, and are at the center of a variety of fields ranging from colloidal chemistry to cell biology, neurodegenerative diseases, and manufacturing of protein drugs. They are frequently measured in different biophysical techniques in the form of second virial coefficients, and nonideality coefficients of sedimentation and diffusion, which may be related mechanistically to macromolecular distance distributions in solution and interparticle potentials. A problem arises for proteins where reversible self-association often complicates the concentration-dependent behavior, such that grossly inconsistent coefficients are measured in experiments based on different techniques, confounding quantitative conclusions. Here we present a global multi-method analysis that synergistically bridges gaps in resolution and sensitivity of orthogonal techniques. We demonstrate the method with a panel of monoclonal antibodies exhibiting different degrees of self-association. We show how their concentration-dependent behavior, examined by static and dynamic light scattering and sedimentation velocity, can be jointly described in a self-consistent framework that separates nonideality coefficients from self-association properties, and thereby extends the quantitative interpretation of nonideality coefficients to probe dynamics in highly concentrated protein solutions.

Conformation of the von Willebrand factor/factor VIII complex in quasi-static flow

Von Willebrand factor (VWF) is a plasma glycoprotein that circulates noncovalently bound to blood coagulation factor VIII (fVIII). VWF is a population of multimers composed of a variable number of ∼280 kDa monomers that is activated in shear flow to bind collagen and platelet glycoprotein Ibα. Electron microscopy, atomic force microscopy, small-angle neutron scattering, and theoretical studies have produced a model in which the conformation of VWF under static conditions is a compact, globular “ball-of-yarn,” implying strong, attractive forces between monomers. We performed sedimentation velocity (SV) analytical ultracentrifugation measurements on unfractionated VWF/fVIII complexes. There was a 20% per mg/ml decrease in the weight-average sedimentation coefficient, s_w, in contrast to the ∼1% per mg/ml decrease observed for compact globular proteins. SV and dynamic light scattering measurements were performed on VWF/fVIII complexes fractionated by size-exclusion chromatography to obtain s_w values and z-average diffusion coefficients, D_z. Molecular weights estimated using these values in the Svedberg equation ranged from 1.7 to 4.1 MDa. Frictional ratios calculated from D_z and molecular weights ranged from 2.9 to 3.4, in contrast to values of 1.1–1.3 observed for globular proteins. The Mark–Houwink–Kuhn–Sakurada scaling relationships between s_w, D_z and molecular weight, s=k′Mas and D=k″MaD, yielded estimates of 0.51 and –0.49 for a_s and a_D, respectively, consistent with a random coil, in contrast to the a_s value of 0.65 observed for globular proteins. These results indicate that interactions between monomers are weak or nonexistent and that activation of VWF is intramonomeric.

Elucidation of the mechanism of subunit exchange in αB crystallin oligomers

AlphaB crystallin (αB-crystallin) is a key protein for maintaining the long-term transparency of the eye lens. In the eye lens, αB-crystallin is a "dynamical" oligomer regulated by subunit exchange between the oligomers. To elucidate the unsettled mechanism of subunit exchange in αB-crystallin oligomers, the study was carried out at two different protein concentrations, 28.5 mg/mL (dense sample) and 0.45 mg/mL (dilute sample), through inverse contrast matching small-angle neutron scattering. Interestingly, the exchange rate of the dense sample was the same as that of the dilute sample. From analytical ultracentrifuge measurements, the coexistence of small molecular weight components and oligomers was detected, regardless of the protein concentration. The model proposed that subunit exchange could proceed through the assistance of monomers and other small oligomers; the key mechanism is attaching/detaching monomers and other small oligomers to/from oligomers. Moreover, this model successfully reproduced the experimental results for both dense and dilute solutions. It is concluded that the monomer and other small oligomers attaching/detaching mainly regulates the subunit exchange in αB-crystallin oligomer.

Quantitative Analysis of Protein Self-Association by Sedimentation Velocity

Sedimentation velocity analytical ultracentrifugation is a powerful classical method to study protein self-association processes in solution based on the size-dependent macromolecular migration in the centrifugal field. This technique can elucidate the assembly scheme, measure affinities ranging from picomolar to millimolar K_d , and in favorable cases provide information on oligomer lifetimes and hydrodynamic shape. The present step-by-step protocols detail the essential steps of instrument calibration, experimental setup, and data analysis. Using a widely available commercial protein as a model system, the protocols invite replication and comparison with our results. A commentary discusses principles for modifications in the protocols that may be necessary to optimize application of sedimentation velocity analysis to other self-associating proteins. ©2020 Wiley Periodicals LLC. Basic Protocol 1: Measurement of external calibration factors Basic Protocol 2: Sedimentation velocity experiment for protein self-association Basic Protocol 3: Sedimentation coefficient distribution analysis in SEDFIT and isotherm analysis in SEDPHAT.

Sedimentation Velocity Methods for the Characterization of Protein Heterogeneity and Protein Affinity Interactions

Sedimentation velocity analytical ultracentrifugation is a powerful and versatile tool for the characterization of proteins and macromolecular complexes in solution. The direct modeling of the sedimentation process using modern computational strategies allows among others to assess the homogeneity/heterogeneity state of protein samples and to characterize protein associations. In this chapter, we will provide theoretical backgrounds and protocols to analyze the size distribution of protein samples and to determine the affinity of protein-protein hetero-associations.

Analysis of nonideality: insights from high concentration simulations of sedimentation velocity data

The Aviv fluorescence detection system (Aviv-FDS) has allowed the performance of sedimentation velocity experiments on therapeutic antibodies in highly concentrated environments like formulation buffers and serum. Methods were implemented in the software package SEDANAL for the analysis of nonideal, weakly associating AUC data acquired on therapeutic antibodies and proteins (Wright et al. Eur Biophys J 47:709–722, 2018, Anal Biochem 550:72–83, 2018). This involved fitting both hydrodynamic, k_s, and thermodynamic, BM₁, nonideality where concentration dependence is expressed as s = s^o/(1 + k_sc) and D = D^o(1 + 2BM₁c)/(1 + k_sc) and s^o and D^o are values extrapolated to c = 0 (mg/ml). To gain insight into the consequences of these phenomenological parameters, we performed simulations with SEDANAL of a monoclonal antibody as a function of k_s (0–100 ml/g) and BM₁ (0–100 ml/g). This provides a visual understanding of the separate and joint impact of k_s and BM₁ on the shape of high-concentration sedimentation velocity boundaries and the challenge of their unique determination by finite element methods. In addition, mAbs undergo weak self- and hetero-association (Yang et al. Prot Sci 27:1334–1348, 2018) and thus we have simulated examples of nonideal weak association over a wide range of concentrations (1–120 mg/ml). Here we demonstrate these data are best analyzed by direct boundary global fitting to models that account for k_s, BM₁ and weak association. Because a typical clinical dose of mAb is 50–200 mg/ml, these results have relevance for biophysical understanding of concentrated therapeutic proteins.

Sedimentation Velocity Analytical Ultracentrifugation Analysis of Marketed Rituximab Drug Product Size Distribution

PURPOSE

Analytical methods suitable for intact drug products are often necessary to evaluate the equivalence in physicochemical properties between two drug products (DP) containing the same drug substance (DS), e.g., an innovator biologic drug and its proposed biosimilar. Analytical Ultracentrifugation (AUC) is a biophysics technique applied to the analysis of size and shape of biomolecules. However, the application of AUC to formulated monoclonal antibody (mAb) DP at high concentration has not been reported.

METHODS

A sedimentation velocity (SV) AUC procedure with a short-pathlength centerpiece was applied to two marketed rituximab DPs, Rituxan® (US) and Reditux® (India), without any buffer exchange or dilution. Detailed precision analysis was performed.

RESULTS

Highly reproducible sedimentation coefficient values (S) and peak areas were obtained for the dominant (> 84%) monomeric rituximab peak. The minor mAb fragment peaks had large variation in both S values and peak areas (3-12%). The identification of oligomer peaks was only reproducible once the abundance was higher than 2%.

CONCLUSIONS

SV-AUC provides an orthogonal characterization tool for protein size distribution, composition and assay, which could be informative for biosimilar drug developers who mostly only have access to formulated mAb. However, AUC needs thorough validation on its accuracy, precision and sensitivity.

Structurally plastic NEMO and oligomerization prone IKK2 subunits define the behavior of human IKK2:NEMO complexes in solution

The human IκB Kinase (IKK) is a multisubunit protein complex of two kinases and one scaffolding subunit that controls induction of transcription factor NF-κB activity. IKK behaves as an entity of aberrantly high apparent molecular weight in solution. Recent X-ray crystallographic and cryo-electron microscopy structures of individual catalytic subunits (IKK1/IKKα and IKK2/IKKβ) reveal that they are both stably folded dimeric proteins that engage in extensive homo-oligomerization through unique surfaces that are required for activation of their respective catalytic activities. The NEMO/IKKγ subunit is a predominantly coiled coil protein that is required for activation of IKK through the canonical NF-κB signaling pathway. Here we report size-exclusion chromatography, multi-angle light scattering, analytical centrifugation, and thermal denaturation analyses of full-length human recombinant NEMO as well as deletion and disease-linked variants. We observe that NEMO is predominantly a dimer in solution, although by virtue of its modular coiled coil regions NEMO exhibits complicated solution dynamics involving portions that are mutually antagonistic toward homodimerization. This behavior causes NEMO to behave as a significantly larger sized particle in solution. Analyses of NEMO in complex with IKK2 indicate that NEMO preserves this structurally dynamic character within the multisubuit complex and provides the complex-bound IKK2 further propensity toward homo-oligomerization. These observations provide critical information on the structural plasticity of NEMO subunit dimers which helps clarify its role in diseases and in IKK regulation through oligomerization-dependent phosphorylation of catalytic IKK2 subunit dimers.

Modelling and Experimental Investigation on the Settling Rate of Kaolinite Particles in Non-Ideal Sedimentation Stage under Constant Gravity

We compared the catalytic effects of two polymers (soluble starch and apple pectin) on the flocculation of kaolinite suspension. Moreover, the relationship between the zeta potential value and the time when kaolin particle sedimentation occurred was verified, and the mechanism of flocculation was analyzed. Additionally, a constitutive model was proposed to simulate the non-ideal sedimentation of clay particles in an aqueous system under constant gravity. This model not only considers the inhomogeneity of the solute but also simulates the change in clay concentration with time during the deposition process. This model proposes a decay constant (α) and sedimentation coefficient (s). The model can also be used to calculate the instantaneous sedimentation rate of the clay suspensions at any time and any depth for the settling cylinder. These sedimentary characteristics were simulated by adopting the established deposition model. The results show that the model is capable of predicting the time required for the complete sedimentation of particles in the aqueous system, suggesting the feasibility of engineering wastewater treatment, site dredging, etc.

High Performance Size Exclusion Chromatography and High-Throughput Dynamic Light Scattering as Orthogonal Methods to Screen for Aggregation and Stability of Monoclonal Antibody Drug Products

The presence of aggregates in monoclonal antibody (mAb) drug product (DP) formulations can present product quality challenges. Here we show that use of High Performance Size Exclusion Chromatography (HP-SEC), in conjunction with high-throughput dynamic light scattering (HT-DLS) analyses of mAb DPs can be a useful strategy to determine monomer content and the presence of aggregates under simulated stress conditions. This analytical approach was used to evaluate four commercially available mAb DPs under different conditions i.e.; original formulations, diluted, and thermo-mechanical stressed condition. Due to particle size limitations of HP-SEC columns, resulting in particles accumulating in the column frits prior to reaching the detector for analysis, there is a possibility that large mAb aggregates may not be detected. Both HP-SEC and HT-DLS were able to detect and resolve the mAb monomer (~10-12 nm) of the DPs in their recommended storage conditions. However, the ability to detect large aggregates (>40 nm) by both analytical methods differed, and HT-DLS was able to detect aggregates between 60 nm and 1400 nm under stress conditions. Our data indicates that HP-SEC, in conjunction with HT-DLS, may be beneficial to detect both mAb DP monomer content and multiple aggregate species (1-1000 nm) in the submicron size range.

Natural Killer Cell Activation Receptor NKp30 Oligomerization Depends on Its N-Glycosylation

NKp30 is one of the main human natural killer (NK) cell activating receptors used in directed immunotherapy. The oligomerization of the NKp30 ligand binding domain depends on the length of the C-terminal stalk region, but our structural knowledge of NKp30 oligomerization and its role in signal transduction remains limited. Moreover, ligand binding of NKp30 is affected by the presence and type of N-glycosylation. In this study, we assessed whether NKp30 oligomerization depends on its N-glycosylation. Our results show that NKp30 forms oligomers when expressed in HEK293S GnTI^- cell lines with simple N-glycans. However, NKp30 was detected only as monomers after enzymatic deglycosylation. Furthermore, we characterized the interaction between NKp30 and its best-studied cognate ligand, B7-H6, with respect to glycosylation and oligomerization, and we solved the crystal structure of this complex with glycosylated NKp30, revealing a new glycosylation-induced mode of NKp30 dimerization. Overall, this study provides new insights into the structural basis of NKp30 oligomerization and explains how the stalk region and glycosylation of NKp30 affect its ligand affinity. This furthers our understanding of the molecular mechanisms involved in NK cell activation, which is crucial for the successful design of novel NK cell-based targeted immunotherapeutics.

Resuspending samples in analytical ultracentrifugation

In analytical ultracentrifugation it is often very useful to resuspend samples in situ after sedimentation experiments for further investigation. This can be achieved by manually subjecting the entire sample cell assembly to gentle motion that causes the air bubble in the sample compartment to repeatedly move through the solution and thereby cause convection. Here we describe a cell mixing device that can accomplish the same through axial rotation and slow rocking motion. This cell mixer is low-cost, open-source, and can be easily assembled from readily available components. It can efficiently mix multiple sample cells side-by-side and may be used with various centerpiece designs.

Revisiting very disperse macromolecule populations in hydrodynamic and light scattering studies of sodium carboxymethyl celluloses

One of the most abundant natural macromolecule, cellulose, is of high importance in technological research including medicine, energy application platforms, and many more. One of its most important ionic derivatives, sodium carboxymethyl cellulose, is known to be very disperse and heterogeneous. The experimental robustness of the methods of hydrodynamics and light scattering are put to test by studying these highly disperse, charged, and heterogeneous macromolecule populations. The following opportunities for molar mass estimations from experimental data were taken into consideration: (i) from the classical Svedberg equation, (ii) from size exclusion chromatography coupled to multi-angle laser light scattering, (iii) from the hydrodynamic invariant, and (iv) the sedimentation parameter. The orthogonality of such approach demonstrates a statistically robust assessment of chain conformational and chain dimensional characteristics of macromolecule populations. Quantitative comparison between the absolute techniques indicates that those have to be checked for accuracy of the obtained and derived characteristics.

Efficient data acquisition with three-channel centerpieces in sedimentation velocity

Three-channel 3D printed centerpieces with two sample sectors next to a joint solvent reference sector were recently described as a strategy to double the throughput of sedimentation velocity analytical ultracentrifugation experiments [Anal. Chem. 91 (2019) 5866-5873]. They are compatible with Rayleigh interference optical detection in commercial analytical ultracentrifuges, but require the rotor angles of data acquisition to be repeatedly adjusted during the experiment to record data from the two sample sectors. Here we present an approach to automate this data acquisition mode through the use of a secondary, general-purpose automation software, and an accompanying data pre-processing software for scan sorting.

A Comprehensive Brownian Dynamics Approach for the Determination of Non-ideality Parameters from Analytical Ultracentrifugation

Brownian dynamics (BD) has been applied as a comprehensive tool to model sedimentation and diffusion of nanoparticles in analytical ultracentrifugation (AUC) experiments. In this article, we extend the BD algorithm by considering space-dependent diffusion and solvent compressibility. With this, the changes in the sedimentation and diffusion coefficient from altered solvent properties at increased pressures are accurately taken into account. Moreover, it is demonstrated how the concept of space-dependent diffusion is employed to describe concentration-dependent sedimentation and diffusion coefficients, in particular, through the Gralen coefficient and the second virial coefficient. The influence of thermodynamic nonideality on diffusional properties can be accurately simulated and agree with well-known evaluation tools. BD simulations for sedimentation equilibrium and sedimentation velocity (SV) AUC experiments including effects of hydrodynamic and thermodynamic nonideality are validated by global evaluation in SEDANAL. The interplay of solvent compressibility and retrieved nonideality parameters can be studied utilizing BD. Finally, the second virial coefficient is determined for lysozyme from SV AUC experiments and BD simulations and compared to membrane osmometry. These results are in line with DLVO theory. In summary, BD simulations are established for the validation of nonideal sedimentation in AUC providing a sound basis for the evaluation of complex interactions even in polydisperse systems.

Impact of soft protein interactions on the excretion, extent of receptor occupancy and tumor accumulation of ultrasmall metal nanoparticles: a compartmental model simulation

Ultrasmall metal nanoparticles (NPs) are next-generation nano-based platforms for in vivo disease diagnosis and treatment. Due to their small size below the kidney filtration threshold and marked resistance to nonspecific serum protein adsorption, ultrasmall NPs can be rapidly excreted through the kidneys and escape liver uptake. However, although ultrasmall particles may be deemed highly resistant to protein adsorption, the real extent of this resistance is not known. Here, a simple compartmental model simulation was therefore implemented to understand how NP behavior in vivo could be modulated by soft, transient NP-plasma protein interactions characterized by dissociation constants in the millimolar range. In Model 1, ultrasmall NPs functionalized with a targeting probe, plasma proteins and target receptors were assumed to co-exist within a single compartment. Simulations were performed to understand the synergistic effect of soft interactions, systemic clearance and NP size on receptor occupancy in the single compartment. The results revealed the existence of a narrow range of ultraweak affinities and optimal particle sizes leading to greater target occupancy. In Model 2, simulations were performed to understand the impact of soft interactions on NP accumulation into a peripheral (tumor) compartment. The results revealed that soft interactions - but not active targeting - enhanced tumor uptake levels when tumor accumulation was limited by 'fast' plasma clearance and 'slow' vascular extravasation. The simple model presented here provides a basic framework to quantitatively understand the blood and tumor pharmacokinetics of ultrasmall NPs under the influence of transient protein interactions.

Enhanced Sample Handling for Analytical Ultracentrifugation with 3D-Printed Centerpieces

The centerpiece of the sample cell assembly in analytical ultracentrifugation holds the sample solution between windows, sealed against high vacuum, and is shaped such that macromolecular migration in centrifugal fields exceeding 200 000g can proceed undisturbed by walls or convection while concentration profiles are imaged with optical detection systems aligned perpendicular to the plane of rotation. We have recently shown that 3D printing using various materials allows inexpensive and rapid manufacturing of centerpieces. In the present work, we expand this endeavor to examine the accuracy of the measured sedimentation process, as well as short-term durability of the centerpieces. We find that 3D-printed centerpieces can be used many times and can provide data equivalent in quality to commonly used commercial epoxy resin centerpieces. Furthermore, 3D printing enables novel designs adapted to particular experimental objectives because they offer unique opportunities, for example, to create well-defined curved surfaces, narrow channels, and embossed features. We present examples of centerpiece designs exploiting these capabilities for improved AUC experiments. This includes narrow sector centerpieces that substantially reduce the required sample volume while maintaining the standard optical path length; thin centerpieces with integrated window holders to provide very short optical pathlengths that reduce optical aberrations at high macromolecular concentrations; long-column centerpieces that increase the observable distance of macromolecular migration for higher-precision sedimentation coefficients; and three-sector centerpieces that allow doubling the number of samples in a single run while reducing the sample volumes. We find each of these designs allows unimpeded macromolecular sedimentation and can provide high-quality sedimentation data.

A Reappraisal of Sedimentation Nonideality Coefficients for the Analysis of Weak Interactions of Therapeutic Proteins

The study of weak or colloidal interactions of therapeutic proteins in different formulations allows prediction and optimization of protein stability. Various biophysical techniques have been applied to determine the second osmotic virial coefficient B₂ as it reflects on the macromolecular distance distribution that governs solution behavior at high concentration. In the present work, we exploit a direct link predicted by hydrodynamic theory between B₂ and the nonideality of sedimentation, commonly measured in sedimentation velocity analytical ultracentrifugation through the nonideality coefficient of sedimentation, k_S. Using sedimentation equilibrium analytical ultracentrifugation for independent measurement of B₂, we have examined the dependence of k_S on B₂ for model proteins in different buffers. The data exhibit the expected linear relationship and highlight the impact of protein shape on the magnitude of the nonideality coefficient k_S. Recently, measurements of k_S have been considerably simplified allowing higher throughput and simultaneous polydispersity assessment at higher protein concentrations. Thus, sedimentation velocity may offer a useful approach to compare the impact of formulation conditions on weak interactions and simultaneously on higher-order structure of therapeutic proteins.

Measuring Ultra-Weak Protein Self-Association by Non-ideal Sedimentation Velocity

Ultra-weak self-association can govern the macroscopic solution behavior of concentrated macromolecular solutions ranging from food products to pharmaceutical formulations and the cytosol. For example, it can promote dynamic assembly of multi-protein signaling complexes, lead to intracellular liquid–liquid phase transitions, and seed crystallization or pathological aggregates. Unfortunately, weak self-association is technically extremely difficult to study, as it requires very high protein concentrations where short intermolecular distances cause strongly correlated particle motion. Additionally, protein samples near their solubility limit in vitro frequently show some degree of polydispersity. Here we exploit the strong mass-dependent separation of assemblies in the centrifugal field to study ultra-weak binding, using a sedimentation velocity technique that allows us to determine particle size distributions while accounting for colloidal hydrodynamic interactions and thermodynamic non-ideality (Chaturvedi, S. K.; et al. Nat. Commun.2018, 9, 4415; DOI: 10.1038/s41467-018-06902-x). We show that this approach, applied to self-associating proteins, can reveal a time-average association state for rapidly reversible self-associations from which the free energy of binding can be derived. The method is label-free and allows studying mid-sized proteins at millimolar protein concentrations in a wide range of solution conditions. We examine the performance of this method with hen egg lysozyme as a model system, reproducing its well-known ionic-strength-dependent weak self-association. The application to chicken γS-crystallin reveals weak monomer–dimer self-association with K_D = 24 mM, corresponding to a standard free energy change of approximately −9 kJ/mol, which is a large contribution to the delicate balance of forces ensuring eye lens transparency.