Molecular basis of high viscosity in concentrated antibody solutions: Strategies for high concentration drug product development

Portable Microfluidic Viscometer for Formulation Development and in Situ Quality Control of Protein and Antibody Solutions

Viscosity of protein solutions is a critical product quality attribute for protein therapeutics such as monoclonal antibodies. Here we introduce a portable single-use analytical chip-based viscometer for determining the viscosity of protein solutions using low sample volumes of 10 μL. Through the combined use of a microfluidic viscometer, a smartphone camera for image capture, and an automated data processing algorithm for the calculation of the viscosity of fluids, we enable measurement of viscosity of multiple samples in parallel. We first validate the viscometer using glycerol-water mixtures and subsequently demonstrate the ability to perform rapid characterization of viscosity in four different monoclonal antibody formulations in a broad concentration (1 to 320 mg/mL) and viscosity (1 to 600 cP) range, showing excellent agreement with values obtained by a conventional cone-plate rheometer. Not only does the platform offer benefits of viscosity measurements using minimal sample volumes, but enables higher throughput compared to gold-standard methodologies owing to multiplexing of the measurement and single-use characteristics of the viscometer, thus showing great promise in developability studies. Additionally, as our platform has the capability of performing viscosity measurements at the point of sample collection, it offers the opportunity to employ viscosity measurement as an in situ quality control of therapeutic proteins and antibodies.

Novel EDTA mediated ethanol protein precipitation method and the application for polysorbate quantification in high protein concentration biopharmaceuticals

Non-ionic surfactants such as Polysorbate 20/ 80 (PS20/ PS80), are commonly used in protein drug formulations to increase protein stability by protecting against interfacial stress and surface absorption. Polysorbate is susceptible to degradation which can impact product stability, leading to the formation of sub-visible and/or visible particles in the drug product during its shelf-life, affecting patient safety and efficacy. Therefore, it is important to monitor polysorbate concentration in drug product formulations of biotherapeutic drugs. The common method for measuring polysorbate concentration in drug product formulations uses mixed mode ion exchange reversed phase HPLC (MAX) coupled to evaporative light scattering detection (ELSD). However, high protein concentration can adversely impact method performance due to high sample viscosity, gel formation, column clogging, interfering peaks and loss of accuracy. To overcome this, a new method was developed based on EDTA mediated ethanol protein precipitation (EDTA/EtOH). This method was successfully implemented for the analysis of polysorbate in antibody formulations with wide range of protein concentration (10-250 mg/mL).

Stability of Protein Pharmaceuticals: Recent Advances

There have been significant advances in the formulation and stabilization of proteins in the liquid state over the past years since our previous review. Our mechanistic understanding of protein-excipient interactions has increased, allowing one to develop formulations in a more rational fashion. The field has moved towards more complex and challenging formulations, such as high concentration formulations to allow for subcutaneous administration and co-formulation. While much of the published work has focused on mAbs, the principles appear to apply to any therapeutic protein, although mAbs clearly have some distinctive features. In this review, we first discuss chemical degradation reactions. This is followed by a section on physical instability issues. Then, more specific topics are addressed: instability induced by interactions with interfaces, predictive methods for physical stability and interplay between chemical and physical instability. The final parts are devoted to discussions how all the above impacts (co-)formulation strategies, in particular for high protein concentration solutions.'

Simultaneous Protein Adsorption and Viscosity Measurement using Micropillar-Enhanced Acoustic Wave (μPAW) Device for Pharmaceutical Applications

At the early stages of drug development, the amount of drug materials is rather limited. In this case, viscosity measurement is often postponed to the later stages, where grams of proteins can be produced. Therefore, it is necessary to develop a viscometer capable of measuring the viscosity with high accuracy while requiring low sample volume. This study presents a novel viscosity measurement technique based on measuring the resonance frequency and motional resistance of a micropillar-enhanced acoustic wave (μPAW) device. The μPAW was developed by fabricating micropillars on the quartz crystal microbalance substrate in order to achieve ultra-high sensitivity, thanks to a unique coupling between the micropillar and the resonator. The experimental measurements demonstrated a nonlinear relationship between the density and viscosity of the fluid and the response of μPAW. A calibration correlation was developed using the response of μPAW in aqueous glycerol and sucrose solutions. The measurements were then extended using high-concentration BSA solutions as the model of protein solution. The main advantage of the μPAW device in this work over other viscometers is the ability to simultaneously measure solution viscosity and protein adsorption on the surface. This is a huge step forward in the development of sensing systems for the pharmaceutical industry, where real-time sensing of target biological proteins and measuring the viscosity of a solution is required.

The search for novel proline analogs for viscosity reduction and stabilization of highly concentrated monoclonal antibody solutions

Administration of monoclonal antibodies (mAbs) is currently focused on subcutaneous injection associated with increased patient adherence and reduced treatment cost, leading to sustainable healthcare. The main bottleneck is low volume that can be injected, requiring highly concentrated mAb solutions. The latter results in increased solution viscosity with pronounced mAb aggregation propensity because of intensive protein-protein interactions. Small molecule excipients have been proposed to restrict the protein-protein interactions, contributing to reduced viscosity. The aim of the study was to discover novel compounds that reduce the viscosity of highly concentrated mAb solution. First, the chemical space of proline analogs was explored and 35 compounds were determined. Viscosity measurements revealed that 18 proline analogs reduced the mAb solution viscosity similar to or more than proline. The compounds forming both electrostatic and hydrophobic interactions with mAb reduced the viscosity of the formulation more efficiently without detrimentally effecting mAb physical stability. A correlation between the level of interaction and viscosity-reducing effect was confirmed with molecular dynamic simulations. Structure rigidity of the compounds and aromaticity contributed to their viscosity-reducing effect, dependent on molecule size. The study results highlight the novel proline analogs as an effective approach in viscosity reduction in development of biopharmaceuticals for subcutaneous administration.

Effect of Annealing on Visible-Bubble Formation and Stability Profiles of Freeze-Dried High Concentration Omalizumab Formulations

In the clinical application of freeze-dried highly concentrated omalizumab formulations, extensive visible bubbles (VBs) can be generated and remain for a long period of time in the reconstitution process, which greatly reduces the clinical use efficiency. It is necessary to understand the forming and breaking mechanism of VBs in the reconstitution process, which is a key factor for efficient and safe administration of biopharmaceutical injection. The effects of different thermal treatments on the volume of VBs and stability of omalizumab, mAb-1, and mAb-2 were investigated. The internal microvoids of the cake were characterized by scanning electron microscopy and mercury intrusion porosimetry. Electron paramagnetic resonance was applied to obtain the molecular mobility of the protein during annealing. A large number of VBs were generated in the reconstitution process of unannealed omalizumab and remained for a long period of time. When annealing steps were added, the volume of VBs was dramatically reduced. When annealed at an aggressive temperature (i.e., -6 °C), although the volume of VBs decreased, the aggregation and acidic species increased significantly. Thus, our observations highlight the importance of setting an additional annealing step with a suitable temperature, which contributes to reducing the VBs while maintaining the stability of the high concentration freeze-dried protein formulation.

Quantifying Protein Shape to Elucidate Its Influence on Solution Viscosity in High-Concentration Electrolyte Solutions

Therapeutic proteins with a high concentration and low viscosity are highly desirable for subcutaneous and certain local injections. The shape of a protein is known to influence solution viscosity; however, the precise quantification of protein shape and its relative impact compared to other factors like charge-charge interactions remains unclear. In this study, we utilized seven model proteins of varying shapes and experimentally determined their shape factors (v) based on Einstein's viscosity theory, which correlate strongly with the ratios of the proteins' surface area to the 2/3 power of their respective volumes, based on protein crystal structures resolved experimentally or predicted by AlphaFold. This finding confirms the feasibility of computationally estimating protein shape factors from amino acid sequences alone. Furthermore, our results demonstrated that, in high-concentration electrolyte solutions, a more spherical protein shape increases the protein's critical concentration (C*), the transition concentration beyond which protein viscosity increases exponentially relative to concentration increases. In summary, our work elucidates protein shape as a key determinant of solution viscosity through quantitative analysis and comparison with other contributing factors. This provides insights into molecular engineering strategies to optimize the molecular design of therapeutic proteins, thus optimizing their viscosity.

Modern day breeding approaches for improvement of castor

Castor (Ricinus communis L.) is an industrially important oil producing crop belongs to Euphorbiaceae family. Castor oil has unique chemical properties make it industrially important crop. It is a member of monotypic genus even though it has ample amount of variability. Using this variability, conventionally many varieties and hybrids have been developed. But, like other crops, the modern and unconventional methods of crop improvement has not fully explored in castor. This article discusses the use of polyploidy induction, distant/wide hybridization and mutation breeding as tools for generating variety. Modern approaches accelerate the speed of crop breeding as an alternative tool. To achieve this goal, molecular markers are employed in breeding to capture the genetic variability through molecular analysis and population structuring. Allele mining is used to trace the evolution of alleles, identify new haplotypes and produce allele specific markers for use in marker aided selection using Genome wide association studies (GWAS) and quantitative trait loci (QTL) mapping. Plant genetic transformation is a rapid and effective mode of castor improvement is also discussed here. The efforts towards developing stable regeneration protocol provide a wide range of utility like embryo rescue in distant crosses, development of somaclonal variation, haploid development using anther culture and callus development for stable genetic transformation has reviewed in this article. Omics has provided intuitions to the molecular mechanisms of (a)biotic stress management in castor along with dissected out the possible genes for improving the yield. Relating genes to traits offers additional scientific inevitability leading to enhancement and sympathetic mechanisms of yield improvement and several stress tolerance.

ProtT5 and random forests-based viscosity prediction method for therapeutic mAbs

Viscosity is a key characteristic of therapeutic antibodies for subcutaneous administration which requires low volume and high concentration formulations. It would be highly beneficial to accurately predict the viscosity of newly developed therapeutic antibodies in the early stages of development. In this work, a ProtT5-XL-UniRef50 (ProtT5) and Random Forests (RF)-based prediction method was proposed for accurately predicting the viscosity of monoclonal antibodies, with only corresponding sequences needed. Starting from the given heavy and light chain V-region sequences, corresponding features were first extracted from the ProtT5 pretrained model. Kernel principal analysis (Kernel-PCA) was then used for reducing the extracted 2048-D (1024-D for each sequence) feature vector to a reasonable level for efficient training of the RF-regressor. Then, the RF model was constructed on 40 commercially available therapeutic antibodies and tested with 3-folds cross-validation. Test results show that the model could reproduce the viscosity value at a high level (Pearson correlation coefficient (PCC) = 0.928). Performance on classifying high (>30 cP) and low (<30 cP) viscosity is much more satisfactory, the Accuracy (ACC) and the area under precision-recall curve (AUC) of the classification model from validation tests are 0.975 and 1.000, respectively. Compared to 5 existing state-of-the-art viscosity prediction methods, the proposed method performs best which would facilitate high concentration antibody viscosity high-throughput screening.

Recent advances in generative biology for biotherapeutic discovery

Generative biology combines artificial intelligence (AI), advanced life sciences technologies, and automation to revolutionize the process of designing novel biomolecules with prescribed properties, giving drug discoverers the ability to escape the limitations of biology during the design of next-generation protein therapeutics. Significant hurdles remain, namely: (i) the inherently complex nature of drug discovery, (ii) the bewildering number of promising computational and experimental techniques that have emerged in the past several years, and (iii) the limited availability of relevant protein sequence-function data for drug-like molecules. There is a need to focus on computational methods that will be most practically effective for protein drug discovery and on building experimental platforms to generate the data most appropriate for these methods. Here, we discuss recent advances in computational and experimental life sciences that are most crucial for impacting the pace and success of protein drug discovery.

Computational Procedure for Predicting Excipient Effects on Protein-Protein Affinities

Protein-protein interactions lie at the center of many biological processes and are a challenge in formulating biological drugs, such as antibodies. A key to mitigating protein association is to use small-molecule additives, i.e., excipients that can weaken protein-protein interactions. Here, we develop a computationally efficient model for predicting the viscosity-reducing effect of different excipient molecules by combining atomic-resolution MD simulations, binding polynomials, and a thermodynamic perturbation theory. In a proof of principle, this method successfully ranks the order of four types of excipients known to reduce the viscosity of solutions of a particular monoclonal antibody. This approach appears useful for predicting the effects of excipients on protein association and phase separation, as well as the effects of buffers on protein solutions.

The Osmolality and Hemolysis of High-Concentration Monoclonal Antibody Formulations

PURPOSES

Highly concentrated monoclonal antibody (mAb) formulations for subcutaneous administration are becoming increasingly preferred within the biopharmaceutical industry for ease of use and improved patient compliance. A common phenomenon observed in the industry is that osmolality detected via freezing-point depression (FPD) in high-concentration mAb formulations is much higher than the theoretical concentrations, yet the occurrence of this phenomenon and its possible safety issues have been rarely reported.

METHODS

The current study summarized theoretical osmolality of U.S. Food and Drug Administration approved high-concentration mAb formulations and evaluated effects of high osmolality on safety using hemolysis experiments for the first time. Two mAbs formulated at 150 mg/mL were used as models and configured into two isotonic solutions: a, a theoretically calculated molarity in the isotonic range (H) and b, an osmolality value measured via the FPD in the isotonic range (I). The H and I formulations of each mAb were individually subjected to hemolysis experiments, and the hemolysis rates of the two formulations of the same mAb were compared. Besides, the effect of mAb concentration on osmolality detected by FPD was explored as well.

RESULTS

The results indicated that the hemolysis rates were similar between the H and I formulations of mAbs at the same sample addition volume, and the osmolality values increased approximately linearly with the increase in mAb concentration.

CONCLUSIONS

High osmolality for high-concentration mAb formulations would not affect product safety and the excipients could be added at relatively high levels to maintain product stability, especially for labile products.

Antibody association in solution: cluster distributions and mechanisms

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

Modulation of the high concentration viscosity of IgG1 antibodies using clinically validated Fc mutations

The self-association of therapeutic antibodies can result in elevated viscosity and create problems in manufacturing and formulation, as well as limit delivery by subcutaneous injection. The high concentration viscosity of some antibodies has been reduced by variable domain mutations or by the addition of formulation excipients. In contrast, the impact of Fc mutations on antibody viscosity has been minimally explored. Here, we studied the effect of a panel of common and clinically validated Fc mutations on the viscosity of two closely related humanized IgG1, κ antibodies, omalizumab (anti-IgE) and trastuzumab (anti-HER2). Data presented here suggest that both Fab-Fab and Fab-Fc interactions contribute to the high viscosity of omalizumab, in a four-contact model of self-association. Most strikingly, the high viscosity of omalizumab (176 cP) was reduced 10.7- and 2.2-fold by Fc modifications for half-life extension (M252Y:S254T:T256E) and aglycosylation (N297G), respectively. Related single mutations (S254T and T256E) each reduced the viscosity of omalizumab by ~6-fold. An alternative half-life extension Fc mutant (M428L:N434S) had the opposite effect in increasing the viscosity of omalizumab by 1.5-fold. The low viscosity of trastuzumab (8.6 cP) was unchanged or increased by ≤ 2-fold by the different Fc variants. Molecular dynamics simulations provided mechanistic insight into the impact of Fc mutations in modulating electrostatic and hydrophobic surface properties as well as conformational stability of the Fc. This study demonstrates that high viscosity of some IgG1 antibodies can be mitigated by Fc mutations, and thereby offers an additional tool to help design future antibody therapeutics potentially suitable for subcutaneous delivery.

Variable domain mutational analysis to probe the molecular mechanisms of high viscosity of an IgG1 antibody

Subcutaneous injection is the preferred route of administration for many antibody therapeutics for reasons that include its speed and convenience. However, the small volume limit (typically ≤ 2 mL) for subcutaneous delivery often necessitates antibody formulations at high concentrations (commonly ≥100 mg/mL), which may lead to physicochemical problems. For example, antibodies with large hydrophobic or charged patches can be prone to self-interaction giving rise to high viscosity. Here, we combined X-ray crystallography with computational modeling to predict regions of an anti-glucagon receptor (GCGR) IgG1 antibody prone to self-interaction. An extensive mutational analysis was undertaken of the complementarity-determining region residues residing in hydrophobic surface patches predicted by spatial aggregation propensity, in conjunction with residue-level solvent accessibility, averaged over conformational ensembles from molecular dynamics simulations. Dynamic light scattering (DLS) was used as a medium throughput screen for self-interaction of ~ 200 anti-GCGR IgG1 variants. A negative correlation was found between the viscosity determined at high concentration (180 mg/mL) and the DLS interaction parameter measured at low concentration (2-10 mg/mL). Additionally, anti-GCGR variants were readily identified with reduced viscosity and antigen-binding affinity within a few fold of the parent antibody, with no identified impact on overall developability. The methods described here may be useful in the optimization of other antibodies to facilitate their therapeutic administration at high concentration.

A computational procedure for predicting excipient effects on protein-protein affinities

Protein-protein interactions lie at the center of much biology and are a challenge in formulating biological drugs such as antibodies. A key to mitigating protein association is to use small molecule additives, i.e. excipients that can weaken protein-protein interactions. Here, we develop a computationally efficient model for predicting the viscosity-reducing effect of different excipient molecules by combining atomic-resolution MD simulations, binding polynomials and a thermodynamic perturbation theory. In a proof of principle, this method successfully rank orders four types of excipients known to reduce the viscosity of solutions of a particular monoclonal antibody. This approach appears useful for predicting effects of excipients on protein association and phase separation, as well as the effects of buffers on protein solutions.

Improved Diffusion Interaction Parameter Measurement to Predict the Viscosity of Concentrated mAb Solutions

During the developability assessment of therapeutic monoclonal antibody (mAb) candidates, utilization of robust high-throughput predictive assays enables rapid selection of top candidates with low risks for late-stage development. Predicting the viscosities of highly concentrated mAbs using limited materials is an important aspect of developability assessment because high viscosity can complicate manufacturability, stability, and administration. Here, we report a high-throughput assay measuring protein-protein interactions to predict mAb viscosity. The diffusion interaction parameter (kD) measures colloidal self-association in dilute solutions and has been reported to be predictive of the mAb viscosity at high concentrations. However, kD of Amgen early stage IgG1 mAb candidates measured in 10 mM acetate at pH 5.2 containing sucrose and polysorbate (denoted A52SuT) shows only weak correlation to their viscosities at 140 mg/mL in A52SuT. We hypothesize that kD measured in A52SuT reflects primarily long-range electrostatic repulsions because most of these mAb candidates carry strong net positive charges in this low ionic strength formulation with pH (5.2) well below pI values of mAb candidates. However, the viscosities of high concentration mAbs depend heavily on short-range molecular interactions. We propose an improved kD method in which salt is added to suppress charge repulsions and to allow for detection of key short-range interactions in dilute solutions. Salt types and salt concentrations were screened, and an optimal salt condition was identified. This optimized method was further validated using two test mAb sets. Overall, the method improves the Pearson R2 between kD and viscosity (6-230 cP) from 0.24 to 0.80 for a data set consisting of 37 mAbs.

Prediction of Antibody Viscosity from Dilute Solution Measurements

The high antibody doses required to achieve a therapeutic effect often necessitate high-concentration products that can lead to challenging viscosity issues in production and delivery. Predicting antibody viscosity in early development can play a pivotal role in reducing late-stage development costs. In recent years, numerous efforts have been made to predict antibody viscosity through dilute solution measurements. A key finding is that the entanglement of long, flexible complexes contributes to the sharp rise in antibody viscosity at the required dosing. This entanglement model establishes a connection between the two-body binding affinity and the many-body viscosity. Exploiting this insight, this study connects dilute solution measurements of self-association to high-concentration viscosity profiles to quantify the relationship between these regimes. The resulting model has exhibited success in predicting viscosity at high concentrations (around 150 mg/mL) from dilute solution measurements, with only a few outliers remaining. Our physics-based approach provides an understanding of fundamental physics, interpretable connections to experimental data, the potential to extrapolate beyond training conditions, and the capacity to effectively explain the physical mechanics behind these outliers. Conducting hypothesis-driven experiments that specifically target the viscosity and relaxation mechanisms of outlier molecules may allow us to unravel the intricacies of their behavior and, in turn, enhance the performance of our model.

Identification, quantification, and structural role of N-glycans in two highly purified isoforms of sheep testicular hyaluronidase

Animal-derived hyaluronidase, which hydrolyzes the polysaccharide hyaluronic acid, has been used in medical applications despite its limited purity. Additionally, the N-glycan characterization of sheep testicular hyaluronidase (STH) and its structural role remain poorly understood. In this study, STH was purified from the commercially available STH preparation (containing at least 14 impurity proteins) using heparin-affinity chromatography followed by size exclusion chromatography. The structure and quantity of N-glycans of STH were investigated using liquid chromatography-electrospray ionization-high energy collision dissociation-tandem mass spectrometry. Two isoforms, H3S1 and H3S2, of STH were obtained (purity >98 %) with a yield of 3.4 % and 5.1 %, respectively. Fourteen N-glycans, including nine core-fucosylated N-glycans (important for the stability and function of glycoproteins), were identified in both H3S1 and H3S2, with similar quantities of each N-glycan. The amino acid sequences of the proteolytic peptides of H3S1 and H3S2 were compared with those reported in STH. The hyaluronic acid-degrading activity of deglycosylated H3S1 and H3S2 was reduced to 70.8 % and 71.1 % compared to that (100 %) of H3S1 and H3S2, respectively. This is the first report of N-glycan characterization of two highly purified isoforms of STH. These H3S1 and H3S2 will be useful for medical use without unwanted effects of partially purified STH.

Predicting Colloidal Stability of High-Concentration Monoclonal Antibody Formulations in Common Pharmaceutical Buffers Using Improved Polyethylene Glycol Induced Protein Precipitation Assay

Colloidal stability is an important consideration when developing high concentration mAb formulations. PEG-induced protein precipitation is a commonly used assay to assess the colloidal stability of protein solutions. However, the practical usefulness and the current theoretical model for this assay have yet to be verified over a large formulation space across multiple mAbs and mAb-based modalities. In the present study, we used PEG-induced protein precipitation assays to evaluate colloidal stability of 3 mAbs in 24 common formulation buffers at 20 and 5 °C. These prediction assays were conducted at low protein concentration (1 mg/mL). We also directly characterized high concentration (100 mg/mL) formulations for cold-induced phase separation, turbidity, and concentratibility by ultrafiltration. This systematic study allowed analysis of the correlation between the results of low concentration assays and the high concentration attributes. The key findings of this study include the following: (1) verification of the usefulness of three different parameters (Cmid, μB, and Tcloud) from PEG-induced protein precipitation assays for ranking colloidal stability of high concentration mAb formulations; (2) a new method to implement PEG-induced protein precipitation assay suitable for high throughput screening with low sample consumption; (3) improvement in the theoretical model for calculating robust thermodynamic parameters of colloidal stability (μB and εB) that are independent of specific experimental settings; (4) systematic evaluation of the effects of pH and buffer salts on colloidal stability of mAbs in common formulation buffers. These findings provide improved theoretical and practical tools for assessing the colloidal stability of mAbs and mAb-based modalities during formulation development.

Characterizing Protein-Protein Interactions and Viscosity of a Monoclonal Antibody from Low to High Concentration Using Small-Angle X-ray Scattering and Molecular Dynamics Simulations

Understanding protein-protein interactions and formation of reversible oligomers (clusters) in concentrated monoclonal antibody (mAb) solutions is necessary for designing stable, low viscosity (η) concentrated formulations for processing and subcutaneous injection. Here we characterize the strength (K) of short-range anisotropic attractions (SRA) for 75-200 mg/mL mAb2 solutions at different pH and cosolute conditions by analyzing structure factors (Seff(q)) from small-angle X-ray scattering (SAXS) using coarse-grained molecular dynamics simulations. Best fit simulations additionally provide cluster size distributions, fractal dimensions, cluster occluded volume, and mAb coordination numbers. These equilibrium properties are utilized in a model to account for increases in viscosity caused by occluded volume in the clusters (packing effects) and dissipation of stress across lubricated fractal clusters. Seff(q) is highly sensitive to K at 75 mg/mL where mAbs can mutually align to form SRA contacts but becomes less sensitive at 200 mg/mL as steric repulsion due to packing becomes dominant. In contrast, η at 200 mg/mL is highly sensitive to SRA and the average cluster size from SAXS/simulation, which is observed to track the cluster relaxation time from shear thinning. By analyzing the distribution of sub-bead hot spots on the 3D mAb surface, we identify a strongly attractive hydrophobic patch in the complementarity determining region (CDR) at pH 4.5 that contributes to the high K and consequently large cluster sizes and high η. Adding NaCl screens electrostatic interactions and increases the impact of hydrophobic attraction on cluster size and raises η, whereas nonspecific binding of Arg attenuates all SRA, reducing η. The hydrophobic patch is absent at higher pH values, leading to smaller K, smaller clusters, and lower η. This work constitutes a first attempt to use SAXS and CG modeling to link both structural and rheological properties of concentrated mAb solutions to the energetics of specific hydrophobic patches on mAb surfaces. As such, our work opens an avenue for future research, including the possibility of designing coarse-grained models with physically meaningful interacting hot spots.

Stimuli-Responsive Hydrogels for Protein Delivery

Proteins and peptides are potential therapeutic agents, but their physiochemical properties make their use as drug substances challenging. Hydrogels are hydrophilic polymeric networks that can swell and retain high amounts of water or biological fluids without being dissolved. Due to their biocompatibility, their porous structure, which enables the transport of various peptides and proteins, and their protective effect against degradation, hydrogels have gained prominence as ideal carriers for these molecules' delivery. Particularly, stimuli-responsive hydrogels exhibit physicochemical transitions in response to subtle modifications in the surrounding environment, leading to the controlled release of entrapped proteins or peptides. This review is focused on the application of these hydrogels in protein and peptide delivery, including a brief overview of therapeutic proteins and types of stimuli-responsive polymers.

A Research Journey: Over a Decade of Denaturing and Native-MS Analyses of Hydrophobic and Membrane Proteins in Amgen Therapeutic Discovery

Membrane proteins and associated complexes currently comprise the majority of therapeutic targets and remain among the most challenging classes of proteins for analytical characterization. Through long-term strategic collaborations forged between industrial and academic research groups, there has been tremendous progress in advancing membrane protein mass spectrometry (MS) analytical methods and their concomitant application to Amgen therapeutic project progression. Herein, I will describe a detailed and personal account of how electrospray ionization (ESI) native mass spectrometry (nMS), ion mobility-MS (IM-MS), reversed phase liquid chromatographic mass spectrometry (RPLC-MS), high-throughput solid phase extraction mass spectrometry, and matrix-assisted laser desorption ionization mass spectrometry methods were developed, optimized, and validated within Amgen Research, and importantly, how these analytical methods were applied for membrane and hydrophobic protein analyses and ultimately therapeutic project support and progression. Additionally, I will discuss all the highly important and productive collaborative efforts, both internal Amgen and external academic, which were key in generating the samples, methods, and associated data described herein. I will also describe some early and previously unpublished nano-ESI (nESI) native-MS data from Amgen Research and the highly productive University of California Los Angeles (UCLA) collaboration. I will also present previously unpublished examples of real-life Amgen biotherapeutic membrane protein projects that were supported by all the MS (and IM) analytical techniques described herein. I will start by describing the initial nESI nMS experiments performed at Amgen in 2011 on empty nanodisc molecules, using a quadrupole time-of-flight MS, and how these experiments progressed on to the 15 Tesla Fourier transform ion cyclotron resonance MS at UCLA. Then described are monomeric and multimeric membrane protein data acquired in both nESI nMS and tandem-MS modes, using multiple methods of ion activation, resulting in dramatic spectral simplification. Also described is how we investigated the far less established and less published subject, that is denaturing RPLC-MS analysis of membrane proteins, and how we developed a highly robust and reproducible RPLC-MS method capable of effective separation of membrane proteins differing in only the presence or absence of an N-terminal post translational modification. Also described is the evolution of the aforementioned RPLC-MS method into a high-throughput solid phase extraction MS method. Finally, I will give my opinion on key developments and how the area of nMS of membrane proteins needs to evolve to a state where it can be applied within the biopharmaceutical research environment for routine therapeutic project support.

How can we discover developable antibody-based biotherapeutics?

Antibody-based biotherapeutics have emerged as a successful class of pharmaceuticals despite significant challenges and risks to their discovery and development. This review discusses the most frequently encountered hurdles in the research and development (R&D) of antibody-based biotherapeutics and proposes a conceptual framework called biopharmaceutical informatics. Our vision advocates for the syncretic use of computation and experimentation at every stage of biologic drug discovery, considering developability (manufacturability, safety, efficacy, and pharmacology) of potential drug candidates from the earliest stages of the drug discovery phase. The computational advances in recent years allow for more precise formulation of disease concepts, rapid identification, and validation of targets suitable for therapeutic intervention and discovery of potential biotherapeutics that can agonize or antagonize them. Furthermore, computational methods for de novo and epitope-specific antibody design are increasingly being developed, opening novel computationally driven opportunities for biologic drug discovery. Here, we review the opportunities and limitations of emerging computational approaches for optimizing antigens to generate robust immune responses, in silico generation of antibody sequences, discovery of potential antibody binders through virtual screening, assessment of hits, identification of lead drug candidates and their affinity maturation, and optimization for developability. The adoption of biopharmaceutical informatics across all aspects of drug discovery and development cycles should help bring affordable and effective biotherapeutics to patients more quickly.

Predictive Nature of High-Throughput Assays in ADC Formulation Screening

Utilization of high-throughput biophysical screening techniques during early screening studies is warranted due to the limited amount of material and large number of samples. But the predictability of the data to longer-term storage stability is critical as the high-throughput methods assist in defining the design space for the longer-term studies. In this study, the biophysical properties of two ADCs in 16 formulation conditions were evaluated using high-throughput techniques. Conformational stability and colloidal stability were evaluated by determining Tm values, kD, B22, and Tagg. In addition, the samples were placed on stability and the extent of aggregate formation over the 8-week interval was determined. The rank order of the 16 different formulations in the high-throughput assays was compared to the rank order observed during the stability studies to assess the predictive capabilities of the screening methods. It was demonstrated that similar rank orders can be expected between high-throughput physical stability indicating assays such as Tagg and B22 and traditional aggregation by SEC data, whereas conformational stability read-outs (Tm) are less predictive. In addition, the high-throughput assays appropriately identified the poor performing formulation conditions, which is ultimately what is desired of screening assays.

Development and Biophysical Characterization of a Humanized FSH-Blocking Monoclonal Antibody Therapeutic Formulated at an Ultra-High Concentration

Highly concentrated antibody formulations are oftentimes required for subcutaneous, self-administered biologics. Here, we report the creation of a unique formulation for our first-in- class FSH-blocking humanized antibody, MS-Hu6, which we propose to move to the clinic for osteoporosis, obesity, and Alzheimer's disease. The studies were carried out using our Good Laboratory Practice (GLP) platform, compliant with the Code of Federal Regulations (Title 21, Part 58). We first used protein thermal shift, size exclusion chromatography, and dynamic light scattering to examine MS-Hu6 concentrations between 1 and 100 mg/mL. We found that thermal, monomeric, and colloidal stability of formulated MS-Hu6 was maintained at a concentration of 100 mg/mL. The addition of the antioxidant L-methionine and chelating agent disodium EDTA improved the formulation's long-term colloidal and thermal stability. Thermal stability was further confirmed by Nano differential scanning calorimetry (DSC). Physiochemical properties of formulated MS-Hu6, including viscosity, turbidity, and clarity, conformed with acceptable industry standards. That the structural integrity of MS-Hu6 in formulation was maintained was proven through Circular Dichroism (CD) and Fourier Transform Infrared (FTIR) spectroscopy. Three rapid freeze-thaw cycles at -80°C/25°C or -80°C/37°C further revealed excellent thermal and colloidal stability. Furthermore, formulated MS-Hu6, particularly its Fab domain, displayed thermal and monomeric storage stability for more than 90 days at 4°C and 25°C. Finally, the unfolding temperature (T _m ) for formulated MS-Hu6 increased by >4.80°C upon binding to recombinant FSH, indicating highly specific ligand binding. Overall, we document the feasibility of developing a stable, manufacturable and transportable MS-Hu6 formulation at a ultra-high concentration at industry standards. The study should become a resource for developing biologic formulations in academic medical centers.

Inverse design of triblock Janus spheres for self-assembly of complex structures in the crystallization slot via digital alchemy

The digital alchemy framework is an extended ensemble simulation technique that incorporates particle attributes as thermodynamic variables, enabling the inverse design of colloidal particles for desired behavior. Here, we extend the digital alchemy framework for the inverse design of patchy spheres that self-assemble into target crystal structures. To constrain the potentials to non-trivial solutions, we conduct digital alchemy simulations with constant second virial coefficient. We optimize the size, range, and strength of patchy interactions in model triblock Janus spheres to self-assemble the 2D kagome and snub square lattices and the 3D pyrochlore lattice, and demonstrate self-assembly of all three target structures with the designed models. The particles designed for the kagome and snub square lattices assemble into high quality clusters of their target structures, while competition from similar polymorphs lower the yield of the pyrochlore assemblies. We find that the alchemically designed potentials do not always match physical intuition, illustrating the ability of the method to find nontrivial solutions to the optimization problem. We identify a window of second virial coefficients that result in self-assembly of the target structures, analogous to the crystallization slot in protein crystallization.

Enhanced Thermal Stability and Reduced Aggregation in an Antibody Fab Fragment at Elevated Concentrations

The aggregation of protein therapeutics such as antibodies remains a major challenge in the biopharmaceutical industry. The present study aimed to characterize the impact of the protein concentration on the mechanisms and potential pathways for aggregation, using the antibody Fab fragment A33 as the model protein. Aggregation kinetics were determined for 0.05 to 100 mg/mL Fab A33, at 65 °C. A surprising trend was observed whereby increasing the concentration decreased the relative aggregation rate, ln(v) (% day^-1), from 8.5 at 0.05 mg/mL to 4.4 at 100 mg/mL. The absolute aggregation rate (mol L^-1 h^-1) increased with the concentration following a rate order of approximately 1 up to a concentration of 25 mg/mL. Above this concentration, there was a transition to an apparently negative rate order of -1.1 up to 100 mg/mL. Several potential mechanisms were examined as possible explanations. A greater apparent conformational stability at 100 mg/mL was observed from an increase in the thermal transition midpoint (T_m) by 7-9 °C, relative to those at 1-4 mg/mL. The associated change in unfolding entropy (△S_vh) also increased by 14-18% at 25-100 mg/mL, relative to those at 1-4 mg/mL, indicating reduced conformational flexibility in the native ensemble. Addition of Tween or the crowding agents Ficoll and dextran, showed that neither surface adsorption, diffusion limitations nor simple volume crowding affected the aggregation rate. Fitting of kinetic data to a wide range of mechanistic models implied a reversible two-state conformational switch mechanism from aggregation-prone monomers (N*) into non-aggregating native forms (N) at higher concentrations. k_D measurements from DLS data also suggested a weak self-attraction while remaining colloidally stable, consistent with macromolecular self-crowding within weakly associated reversible oligomers. Such a model is also consistent with compaction of the native ensemble observed through changes in T_m and △S_vh.

Patient-Centric Product Development: A Summary of Select Regulatory CMC and Device Considerations

Patient-centric drug development describes the systematic approach to incorporating the patient's perspectives and preferences into the design, assessment, and production of a therapeutic product. While a patient centric approach can be applied at any stage of the drug development lifecycle, an integrated end-to-end strategy is often most effective to create an optimized product for the patient at the earliest possible timepoint. The importance of patient centricity is well recognized by health authorities and biopharmaceutical organizations which have established toolsets, guidances, and methodologies for incorporating patient input during the clinical stage of development. However, in addition to clinical research, there are other significant aspects of product development that profoundly impact the patient experience. Specifically, chemistry, manufacturing, and control (CMC) and device aspects must also be acknowledged and addressed as part of a cohesive patient-centric development strategy. This review explores current applications and regulatory considerations for patient-centric approaches across the product lifecycle, including R&D, early product development, clinical development, device and combination product development, and post-approval change management. Specific topics of discussion include the contributions of product modality, formulation, and devices to the patient experience; usage of the Quality Target Product Profile (QTPP) as a patient-centered design tool; and post-approval product optimization. Future advancements in regulatory data management and information exchange are also explored as potential enablers of patient engagement which support enhanced communication and interconnectivity between stakeholders. Multidisciplinary collaboration between patients, health authorities, health care providers, and the biopharmaceutical industry is ultimately necessary for ensuring that medicinal products, and their corresponding regulatory processes, take on a patient-first mindset that prioritizes patient needs, values, and preferences.

Characterizing Experimental Monoclonal Antibody Interactions and Clustering Using a Coarse-Grained Simulation Library and a Viscosity Model

Attractive protein-protein interactions in concentrated monoclonal antibody (mAb) solutions may lead to the formation of clusters that increase viscosity. Here, we propose an analytical model that relates mAb solution viscosity to clustering by accounting for the contributions of suboptimal mAb packing within a cluster and cluster fractal dimension. The influence of short-range, anisotropic attractions and long-range Coulombic repulsion on cluster properties is investigated by analyzing the cluster-size distributions, cluster fractal dimensions, radial distribution functions, and static structure factors from a library of coarse-grained molecular dynamics simulations. The library spans a vast range of mAb charges and attractive interactions in solutions of varying ionic strength. We present a framework for combining the viscosity model and simulation library to successfully characterize the attraction, repulsion, and clustering of an experimental mAb in three different pH and cosolute conditions by fitting the measured viscosity or structure factor from small-angle X-ray scattering. At low ionic strength, the cluster-size distribution is impacted by strong charges, and both the viscosity and net charge or structure factor and net charge must be considered to deconvolute the effects of short-range attraction and long-range repulsion.

AF4 and PEG Precipitation as Predictive Assays for Antibody Self-Association

Monoclonal antibodies (mAbs) are often formulated as high-protein-concentration solutions, which in some cases can exhibit physical stability issues such as high viscosity and opalescence. To ensure that mAb-based drugs can meet their manufacturing, stability, and delivery requirements, it is advantageous to screen for and select mAbs during discovery that are not prone to such behaviors. It has been recently shown that both these macroscopic properties can be predicted to a certain extent from the diffusion interaction parameter (k_D), which is a measure of self-association under dilute conditions.¹ However, k_D can be challenging to measure at the early stage of discovery, where a relatively large amount of a high-purity material, which is required by traditional methods, is often not available. In this study, we demonstrate asymmetric field-flow fractionation (AF4) as a tool to measure self-association and therefore identify antibodies with problematic issues at high concentrations. The principle lies on the ability to concentrate the sample close to the membrane during the injection mode, which can reach formulation-relevant concentrations (>100 mg/mL).² By analyzing a well-characterized library of commercial antibodies, we show that the measured retention time of the antibodies allows us to pinpoint molecules that exhibit issues at high concentrations. Remarkably, our AF4 assay requires very little (30 μg) sample under dilute conditions and does not need extensive sample purification. Furthermore, we show that a polyethylene glycol (PEG) precipitation assay provides results consistent with AF4 and moreover can further differentiate molecules with issues of opalescence or high viscosity. Overall, our results delineate a two-step strategy for the identification of problematic variants at high concentrations, with AF4 for early developability screening, followed by a PEG assay to validate the problematic molecules and further discriminate between opalescence or high-viscosity issues. This two-step antibody selection strategy enables us to select antibodies early in the discovery process, which are compatible with high-concentration formulations.

An Intercompany Perspective on Practical Experiences of Predicting, Optimizing and Analyzing High Concentration Biologic Therapeutic Formulations

Developing high-dose biologic drugs for subcutaneous injection often requires high-concentration formulations and optimizing viscosity, solubility, and stability while overcoming analytical, manufacturing, and administration challenges. To understand industry approaches for developing high-concentration formulations, the Formulation Workstream of the BioPhorum Development Group, an industry-wide consortium, conducted an inter-company collaborative exercise which included several surveys. This collaboration provided an industry perspective, experience, and insight into the practicalities for developing high-concentration biologics. To understand solubility and viscosity, companies desire predictive tools, but experience indicates that these are not reliable and experimental strategies are best. Similarly, most companies prefer accelerated and stress stability studies to in-silico or biophysical-based prediction methods to assess aggregation. In addition, optimization of primary container-closure and devices are pursued to mitigate challenges associated with high viscosity of the formulation. Formulation strategies including excipient selection and application of studies at low concentration to high-concentration formulations are reported. Finally, analytical approaches to high concentration formulations are presented. The survey suggests that although prediction of viscosity, solubility, and long-term stability is desirable, the outcome can be inconsistent and molecule dependent. Significant experimental studies are required to confirm robust product definition as modeling at low protein concentrations will not necessarily extrapolate to high concentration formulations.

Predictive modeling of concentration-dependent viscosity behavior of monoclonal antibody solutions using artificial neural networks

Solutions of monoclonal antibodies (mAbs) can show increased viscosity at high concentration, which can be a disadvantage during protein purification, filling, and administration. The viscosity is determined by protein-protein-interactions, which are influenced by the antibody's sequence as well as solution conditions, like pH, buffer type, or the presence of salts and other excipients. To predict viscosity, experimental parameters, like the diffusion interaction parameter (kD), or computational tools harnessing information derived from primary sequence, are often used, but a reliable predictive tool is still missing. We present a modeling approach employing artificial neural networks (ANNs) using experimental factors combined with simulation-derived parameters plus viscosity data from 27 highly concentrated (180 mg/mL) mAbs. These ANNs can be used to predict if mAbs exhibit problematic viscosity at distinct concentrations or to model viscosity-concentration-curves.

Assessing developability early in the discovery process for novel biologics

Beyond potency, a good developability profile is a key attribute of a biological drug. Selecting and screening for such attributes early in the drug development process can save resources and avoid costly late-stage failures. Here, we review some of the most important developability properties that can be assessed early on for biologics. These include the influence of the source of the biologic, its biophysical and pharmacokinetic properties, and how well it can be expressed recombinantly. We furthermore present in silico, in vitro, and in vivo methods and techniques that can be exploited at different stages of the discovery process to identify molecules with liabilities and thereby facilitate the selection of the most optimal drug leads. Finally, we reflect on the most relevant developability parameters for injectable versus orally delivered biologics and provide an outlook toward what general trends are expected to rise in the development of biologics.

Development of in silico models to predict viscosity and mouse clearance using a comprehensive analytical data set collected on 83 scaffold-consistent monoclonal antibodies

Biologic drug discovery pipelines are designed to deliver protein therapeutics that have exquisite functional potency and selectivity while also manifesting biophysical characteristics suitable for manufacturing, storage, and convenient administration to patients. The ability to use computational methods to predict biophysical properties from protein sequence, potentially in combination with high throughput assays, could decrease timelines and increase the success rates for therapeutic developability engineering by eliminating lengthy and expensive cycles of recombinant protein production and testing. To support development of high-quality predictive models for antibody developability, we designed a sequence-diverse panel of 83 effector functionless IgG1 antibodies displaying a range of biophysical properties, produced and formulated each protein under standard platform conditions, and collected a comprehensive package of analytical data, including in vitro assays and in vivo mouse pharmacokinetics. We used this robust training data set to build machine learning classifier models that can predict complex protein behavior from these data and features derived from predicted and/or experimental structures. Our models predict with 87% accuracy whether viscosity at 150 mg/mL is above or below a threshold of 15 centipoise (cP) and with 75% accuracy whether the area under the plasma drug concentration-time curve (AUC0-672 h) in normal mouse is above or below a threshold of 3.9 × 106 h x ng/mL.

Blueprint for antibody biologics developability

Large-molecule antibody biologics have revolutionized medicine owing to their superior target specificity, pharmacokinetic and pharmacodynamic properties, safety and toxicity profiles, and amenability to versatile engineering. In this review, we focus on preclinical antibody developability, including its definition, scope, and key activities from hit to lead optimization and selection. This includes generation, computational and in silico approaches, molecular engineering, production, analytical and biophysical characterization, stability and forced degradation studies, and process and formulation assessments. More recently, it is apparent these activities not only affect lead selection and manufacturability, but ultimately correlate with clinical progression and success. Emerging developability workflows and strategies are explored as part of a blueprint for developability success that includes an overview of the four major molecular properties that affect all developability outcomes: 1) conformational, 2) chemical, 3) colloidal, and 4) other interactions. We also examine risk assessment and mitigation strategies that increase the likelihood of success for moving the right candidate into the clinic.

A Novel MATLAB®-Algorithm-Based Video Analysis to Quantitatively Determine Solution Creeping in Intact Pharmaceutical Glass Vials

During the filling process of a biopharmaceutical drug product (DP), a liquid DP film might creep up the inner vial wall which is barely discernible, appears as milky-white haze after lyophilisation and is known as fogging. Creeping and fogging are mainly dependent on the primary packaging material surface and its hydration, vial preparation process as well as DP composition. The occurrence of both can impede visual inspection and might lead to DP rejection. Hence, our studies focused on the early detection of liquid solution and glass vial surface interaction directly after filling. For a fast and highly sensitive evaluation a novel video-based analysis was used. To our knowledge, this is the first time a MATLAB®-algorithm-based video analysis was applied to quantitatively determine creeping time-resolved. Furthermore, creeping in dependence of vial processing sites, surfactant type and concentration, filling temperature, and vial format were investigated. The results were verified using orthogonal conventional methods such as surface tension, wetting behaviour, and contact angle measurements, as well as ToF-SIMS, ICP-MS, and SEM. Additionally, the methods applied were assessed regarding their cross-validation capability. The observations indicate that the vial preparation process can have a pronounced impact on alteration of the glass vial surface and related creeping behaviour of the filled solution.

Effects of Monovalent Salt on Protein-Protein Interactions of Dilute and Concentrated Monoclonal Antibody Formulations

In this study, we used sodium chloride (NaCl) to extensively modulate non-specific protein-protein interactions (PPI) of a humanized anti-streptavidin monoclonal antibody class 2 molecule (ASA-IgG2). The changes in PPI with varying NaCl (C_NaCl) and monoclonal antibody (mAb) concentration (C_mAb) were assessed using the diffusion interaction parameter k_D and second virial coefficient B₂₂ measured from solutions with low to moderate C_mAb. The effective structure factor S(q)_eff measured from concentrated mAb solutions using small-angle X-ray and neutron scattering (SAXS/SANS) was also used to characterize the PPI. Our results found that the nature of net PPI changed not only with C_NaCl, but also with increasing C_mAb. As a result, parameters measured from dilute and concentrated mAb samples could lead to different predictions on the stability of mAb formulations. We also compared experimentally determined viscosity results with those predicted from interaction parameters, including k_D and S(q)_eff. The lack of a clear correlation between interaction parameters and measured viscosity values indicates that the relationship between viscosity and PPI is concentration-dependent. Collectively, the behavior of flexible mAb molecules in concentrated solutions may not be correctly predicted using models where proteins are considered to be uniform colloid particles defined by parameters derived from low C_mAb.

Computational models for studying physical instabilities in high concentration biotherapeutic formulations

Computational prediction of the behavior of concentrated protein solutions is particularly advantageous in early development stages of biotherapeutics when material availability is limited and a large set of formulation conditions needs to be explored. This review provides an overview of the different computational paradigms that have been successfully used in modeling undesirable physical behaviors of protein solutions with a particular emphasis on high-concentration drug formulations. This includes models ranging from all-atom simulations, coarse-grained representations to macro-scale mathematical descriptions used to study physical instability phenomena of protein solutions such as aggregation, elevated viscosity, and phase separation. These models are compared and summarized in the context of the physical processes and their underlying assumptions and limitations. A detailed analysis is also given for identifying protein interaction processes that are explicitly or implicitly considered in the different modeling approaches and particularly their relations to various formulation parameters. Lastly, many of the shortcomings of existing computational models are discussed, providing perspectives and possible directions toward an efficient computational framework for designing effective protein formulations.

Design of a Novel Fab-Like Antibody Fragment with Enhanced Stability and Affinity for Clinical use

With increasing interest in applying recombinant monoclonal antibodies (mAbs) in human medicine, engineered mAb fragments with reduced size and improved stability are in demand to overcome current limitations in clinical use. Herein, a novel Fab-like antibody fragment generated via an in silico-based engineering approach where the CH1 and CL domains of Fab are replaced by the IgG1 CH3 domains is described. This construct, designated as FabCH3, maintains the natural N-terminus and C-terminus of IgG antibody, can be expressed at a high level in bacterial cells and, importantly, exhibits much higher stability and affinity than the parental Fab when tested in a mesothelin-specific Fab m912, as well as a vascular endothelial growth factor A (VEGFA)-specific Fab Ranibizumab (in vivo). The high-resolution crystal structures of m912 FabCH3 and m912 Fab are determined, and the comparative analysis reveals more rigid structures in both constant domains and complementarity-determining regions of FabCH3, explaining its enhanced stability and affinity. Overall, the stabilized FabCH3 described in this report provides a versatile platform for engineering Fab-like antibody fragments with higher stability and antigen-binding affinity that can be used as a distinct class of antibody therapeutics.

Machine learning prediction of antibody aggregation and viscosity for high concentration formulation development of protein therapeutics

Machine learning has been recently used to predict therapeutic antibody aggregation rates and viscosity at high concentrations (150 mg/ml). These works focused on commercially available antibodies, which may have been optimized for stability. In this study, we measured accelerated aggregation rates at 45°C and viscosity at 150 mg/ml for 20 preclinical and clinical-stage antibodies. Features obtained from molecular dynamics simulations of the full-length antibody and sequences were used for machine learning model construction. We found a k-nearest neighbors regression model with two features, spatial positive charge map on the CDRH2 and solvent-accessible surface area of hydrophobic residues on the variable fragment, gives the best performance for predicting antibody aggregation rates (r = 0.89). For the viscosity classification model, the model with the highest accuracy is a logistic regression model with two features, spatial negative charge map on the heavy chain variable region and spatial negative charge map on the light chain variable region. The accuracy and the area under precision recall curve of the classification model from validation tests are 0.86 and 0.70, respectively. In addition, we combined data from another 27 commercial mAbs to develop a viscosity predictive model. The best model is a logistic regression model with two features, number of hydrophobic residues on the light chain variable region and net charges on the light chain variable region. The accuracy and the area under precision recall curve of the classification model are 0.85 and 0.6, respectively. The aggregation rates and viscosity models can be used to predict antibody stability to facilitate pharmaceutical development.

Dipolar Interactions and Protein Hydration in Highly Concentrated Antibody Formulations

Molecular interaction mechanisms in high-concentrated protein systems are of fundamental importance for the rational development of biopharmaceuticals such as monoclonal antibody (mAb) formulations. In such high-concentrated protein systems, the intermolecular distances between mAb molecules are reduced to the size of the protein diameter (approx. 10 nm). Thus, protein-protein interactions are more pronounced at high concentrations; so a direct extrapolation of physicochemical properties obtained from measurements at a low protein concentration of the corresponding properties at a high protein concentration is highly questionable. Besides the charge-charge interaction, the effects of molecular crowding, dipolar interaction, changes in protein hydration, and self-assembling tendency become more relevant. Here, protein hydration, protein dipole moment, and protein-protein interactions were studied in protein concentrations up to 200 mg/mL (= 1.3 mM) in different formulations for selected mAbs using dielectric relaxation spectroscopy (DRS). These data are correlated with the second virial coefficient, A₂, the diffusion interaction parameter, k_D, the elastic shear modulus, G', and the dynamic viscosity, η. When large contributions of dipolar protein-protein interactions were observed, the tendency of self-assembling and an increase in solution viscosity were detected. These effects were examined using specific buffer conditions. Furthermore, different types of protein-water interactions were identified via DRS, whereby the effect of high protein concentration on protein hydration was investigated for different high-concentrated liquid formulations (HCLFs).

Assessment of Therapeutic Antibody Developability by Combinations of In Vitro and In Silico Methods

Although antibodies have become the fastest-growing class of therapeutics on the market, it is still challenging to develop them for therapeutic applications, which often require these molecules to withstand stresses that are not present in vivo. We define developability as the likelihood of an antibody candidate with suitable functionality to be developed into a manufacturable, stable, safe, and effective drug that can be formulated to high concentrations while retaining a long shelf life. The implementation of reliable developability assessments from the early stages of antibody discovery enables flagging and deselection of potentially problematic candidates, while focussing available resources on the development of the most promising ones. Currently, however, thorough developability assessment requires multiple in vitro assays, which makes it labor intensive and time consuming to implement at early stages. Furthermore, accurate in vitro analysis at the early stage is compromised by the high number of potential candidates that are often prepared at low quantities and purity. Recent improvements in the performance of computational predictors of developability potential are beginning to change this scenario. Many computational methods only require the knowledge of the amino acid sequences and can be used to identify possible developability issues or to rank available candidates according to a range of biophysical properties. Here, we describe how the implementation of in silico tools into antibody discovery pipelines is increasingly offering time- and cost-effective alternatives to in vitro experimental screening, thus streamlining the drug development process. We discuss in particular the biophysical and biochemical properties that underpin developability potential and their trade-offs, review various in vitro assays to measure such properties or parameters that are predictive of developability, and give an overview of the growing number of in silico tools available to predict properties important for antibody development, including the CamSol method developed in our laboratory.

Deciphering the high viscosity of a therapeutic monoclonal antibody in high concentration formulations by microdialysis-hydrogen/deuterium exchange mass spectrometry

High concentration formulations of therapeutic monoclonal antibodies (mAbs) are highly desired for subcutaneous injection. However, high concentration formulations can exhibit unusual molecular behaviors, such as high viscosity or aggregation, that present challenges for manufacturing and administration. To understand the molecular mechanism of the high viscosity exhibited by high concentration protein formulations, we analyzed a human IgG4 (mAb1) at high protein concentrations using sedimentation velocity analytical ultracentrifugation (SV-AUC), X-ray crystallography, hydrogen/deuterium exchange mass spectrometry (HDX-MS), and protein surface patches analysis. Particularly, we developed a microdialysis HDX-MS method to determine intermolecular interactions at different protein concentrations. SV-AUC revealed that mAb1 displayed a propensity for self-association of Fab-Fab, Fab-Fc, and Fc-Fc. mAb1 crystal structure and HDX-MS results demonstrated self-association between complementarity-determining regions (CDRs) and Fc through electrostatic interactions. HDX-MS also indicated Fab-Fab interactions through hydrophobic surface patches constructed by mAb1 CDRs. Our multi-method approach, including fast screening of SV-AUC as well as interface analysis by X-ray crystallography and HDX-MS, helped to elucidate the high viscosity of mAb1 at high concentrations as induced by self-associations of Fab-Fc and Fab-Fab.

Progress and challenges for the machine learning-based design of fit-for-purpose monoclonal antibodies

Although the therapeutic efficacy and commercial success of monoclonal antibodies (mAbs) are tremendous, the design and discovery of new candidates remain a time and cost-intensive endeavor. In this regard, progress in the generation of data describing antigen binding and developability, computational methodology, and artificial intelligence may pave the way for a new era of in silico on-demand immunotherapeutics design and discovery. Here, we argue that the main necessary machine learning (ML) components for an in silico mAb sequence generator are: understanding of the rules of mAb-antigen binding, capacity to modularly combine mAb design parameters, and algorithms for unconstrained parameter-driven in silico mAb sequence synthesis. We review the current progress toward the realization of these necessary components and discuss the challenges that must be overcome to allow the on-demand ML-based discovery and design of fit-for-purpose mAb therapeutic candidates.

New insight into the importance of formulation variables on parenteral growth hormone preparations: potential effect on the injection-site pain

Reducing injection-site pain (ISP) in patients with chronic conditions such as growth hormone deficiency is a valuable strategy to improve patient compliance and therapeutic efficiency. Thus understanding different aspects of pain induction following subcutaneous injection of biotherapeutics and identifying the responsible factors are vital. Here we have discussed the effects of formulation's viscosity, concentration, osmolality, buffering agents, pH, and temperature as well as injection volume, dosing frequency, and different excipients on ISP following subcutaneous injection of commercially available recombinant human growth hormone products. Our literature review found limited available data on the effects of different components of parenteral rhGH products on ISP. This may be due to high cost associated with conducting various clinical trials to assess each excipient in the formulation or to determine the complex interactions of different components and its impact on ISP. Recently, conducting molecular dynamics simulation studies before formulation design has been recommended as an alternative and less-expensive approach. On the other hand, the observed inconsistencies in the available data is mainly due to different pain measurement approaches used in each study. Moreover, it is difficult to translate data obtained from animal studies to human subjects. Despite all these limitations, our investigation showed that components of parenteral rhGH products can significantly contribute to ISP. We suggest further investigation is required for development of long acting, buffer-free, preservative-free formulations. Besides, various excipients are currently being investigated for reducing ISP which can be used as alternatives for common buffers, surfactants or preservatives in designing future rhGH formulations.

Separating clinical antibodies from repertoire antibodies, a path to in silico developability assessment

Approaches for antibody discovery have seen substantial improvement and success in recent years. Yet, advancing antibodies into the clinic remains difficult because therapeutic developability concerns are challenging to predict. We developed a computational model to simplify antibody developability assessment and enable accelerated early-stage screening. To this end, we quantified the ability of hundreds of sequence- and structure-based descriptors to differentiate clinical antibodies that have undergone rigorous screening and characterization for drug-like properties from antibodies in the human repertoire that are not natively paired. This analysis identified 144 descriptors capable of distinguishing clinical from repertoire antibodies. Five descriptors were selected and combined based on performance and orthogonality into a single model referred to as the Therapeutic Antibody Developability Analysis (TA-DA). On a hold-out test set, this tool separated clinical antibodies from repertoire antibodies with an AUC = 0.8, demonstrating the ability to identify developability attributes unique to clinical antibodies. Based on our results, the TA-DA score may serve as an approach for selecting lead antibodies for further development.

Abbreviations: Affinity-Capture Self-Interaction Nanoparticle Spectroscopy (AC-SINS), Area Under the Curve (AUC), Complementary-Determining Region (CDR), Clinical-Stage Therapeutics (CST), Framework (FR), Monoclonal Antibodies (mAbs), Observed Antibody Space (OAS), Receiver Operating Characteristic (ROC), Size-Exclusion Chromatography (SEC), Structural Aggregation Propensity (SAP), Therapeutic Antibody Developability Analysis (TA-DA), Therapeutic Antibody Profiler (TAP), Therapeutic Structural Antibody Database (Thera-SAbDab), Variable Heavy (VH), Variable Light (VL).