Long- and short-range electrostatic interactions affect the rheology of highly concentrated antibody solutions

FMO-guided design of darunavir analogs as HIV-1 protease inhibitors

The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.

Poly(glutamic acid)-Based Viscosity Reducers for Concentrated Formulations of a Monoclonal IgG Antibody

Above a concentration threshold, the viscosity of solutions of proteins increases abruptly, which hampers the injectability of therapeutic formulations. Concentrations above 200 g/L are an ideal goal for subcutaneous application of antibodies. Molecular additives, such as amino acids (e.g., arginine) help decrease the viscosity, but they are used at concentrations as high as about 200 mmol/L. We addressed the question of whether poly(amino acids) could be more efficient than small molecular additives. We observed marked fluidification of a model therapeutic monoclonal antibody (mAb) solution by poly(d,l-glutamic acid) and poly(l-glutamic acid) derivatives added at concentrations of <6.5 g/L (i.e., a mAb/polymer chain molar ratio between 4:1 and 1:1 mol/mol). The bare poly(glutamate) parent chains were compared with polyethylene glycol-grafted chains as PEGylation is a common way to enhance stability. Viscosity could be decreased to ∼20 mPa s as compared to values of ∼100 mPa s in the absence of polymers at 200 g/L mAb. Formation of complexes between the mAb and the polyglutamates was characterized by capillary electrophoresis analysis in dilute solutions (1 g/L mAb) and by observation of phase separation at higher concentrations, suggesting tight association at about 2:1 mol/mol mAb/polymer. Altogether, these results show that polyglutamate derivatives hold an untapped potential as an excipient for fluidification of concentrated protein solutions.

Molecular surface descriptors to predict antibody developability: sensitivity to parameters, structure models, and conformational sampling

In silico assessment of antibody developability during early lead candidate selection and optimization is of paramount importance, offering a rapid and material-free screening approach. However, the predictive power and reproducibility of such methods depend heavily on the selection of molecular descriptors, model parameters, accuracy of predicted structure models, and conformational sampling techniques. Here, we present a set of molecular surface descriptors specifically designed for predicting antibody developability. We assess the performance of these descriptors by benchmarking their correlations with an extensive array of experimentally determined biophysical properties, including viscosity, aggregation, hydrophobic interaction chromatography, human pharmacokinetic clearance, heparin retention time, and polyspecificity. Further, we investigate the sensitivity of these surface descriptors to methodological nuances, such as the choice of interior dielectric constant, hydrophobicity scales, structure prediction methods, and the impact of conformational sampling. Notably, we observe systematic shifts in the distribution of surface descriptors depending on the structure prediction method used, driving weak correlations of surface descriptors across structure models. Averaging the descriptor values over conformational distributions from molecular dynamics mitigates the systematic shifts and improves the consistency across different structure prediction methods, albeit with inconsistent improvements in correlations with biophysical data. Based on our benchmarking analysis, we propose six in silico developability risk flags and assess their effectiveness in predicting potential developability issues for a set of case study molecules.

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.

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.

A comprehensive evaluation of arginine and its derivatives as protein formulation stabilizers

Arginine and its derivatives (such as arginine ethyl ester and acetyl arginine) have varying degrees of protein aggregation suppressor effect across different protein solutions. To understand this performance ambiguity, we evaluated the activity of arginine, acetyl arginine, and arginine ethyl ester for aggregation suppressor effect against human intravenous immunoglobulin G (IgG) solution at pH 4.8. Both arginine and its cationic derivative arginine ethyl ester in their hydrochloride salt forms significantly reduced the colloidal and conformational stability (reduced kd and Tm) of IgG. Consequently, the monomer content was decreased with an increase in subvisible particulates after agitation or thermal stress. Furthermore, compared to arginine, arginine ethyl ester with one more cationic charge and hydrochloride salt form readily precipitated IgG at temperatures higher than 25 °C. On the contrary, acetyl arginine, which mostly exists in a neutral state at pH 4.8, efficiently suppressed the formation of subvisible particles retaining a high amount of monomer owing to its higher colloidal and conformational stability. Concisely, the charged state of additives significantly impacts protein stability. This study demonstrated that contrary to popular belief, arginine and its derivatives may either enhance or suppress protein aggregation depending on their net charge and concentration.

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.

Viscosity Prediction of High-Concentration Antibody Solutions with Atomistic Simulations

The computational prediction of the viscosity of dense protein solutions is highly desirable, for example, in the early development phase of high-concentration biopharmaceutical formulations where the material needed for experimental determination is typically limited. Here, we use large-scale atomistic molecular dynamics (MD) simulations with explicit solvation to de novo predict the dynamic viscosities of solutions of a monoclonal IgG1 antibody (mAb) from the pressure fluctuations using a Green-Kubo approach. The viscosities at simulated mAb concentrations of 200 and 250 mg/mL are compared to the experimental values, which we measured with rotational rheometry. The computational viscosity of 24 mPa·s at the mAb concentration of 250 mg/mL matches the experimental value of 23 mPa·s obtained at a concentration of 213 mg/mL, indicating slightly different effective concentrations (or activities) in the MD simulations and in the experiments. This difference is assigned to a slight underestimation of the effective mAb-mAb interactions in the simulations, leading to a too loose dynamic mAb network that governs the viscosity. Taken together, this study demonstrates the feasibility of all-atom MD simulations for predicting the properties of dense mAb solutions and provides detailed microscopic insights into the underlying molecular interactions. At the same time, it also shows that there is room for further improvements and highlights challenges, such as the massive sampling required for computing collective properties of dense biomolecular solutions in the high-viscosity regime with reasonable statistical precision.

A workflow for the development of template-assisted membrane crystallization downstream processing for monoclonal antibody purification

Monoclonal antibodies (mAbs) are commonly used biologic drugs for the treatment of diseases such as rheumatoid arthritis, multiple sclerosis, COVID-19 and various cancers. They are produced in Chinese hamster ovary cell lines and are purified via a number of complex and expensive chromatography-based steps, operated in batch mode, that rely heavily on protein A resin. The major drawback of conventional procedures is the high cost of the adsorption media and the extensive use of chemicals for the regeneration of the chromatographic columns, with an environmental cost. We have shown that conventional protein A chromatography can be replaced with a single crystallization step and gram-scale production can be achieved in continuous flow using the template-assisted membrane crystallization process. The templates are embedded in a membrane (e.g., porous polyvinylidene fluoride with a layer of polymerized polyvinyl alcohol) and serve as nucleants for crystallization. mAbs are flexible proteins that are difficult to crystallize, so it can be challenging to determine the optimal conditions for crystallization. The objective of this protocol is to establish a systematic and flexible approach for the design of a robust, economic and sustainable mAb purification platform to replace at least the protein A affinity stage in traditional chromatography-based purification platforms. The procedure provides details on how to establish the optimal parameters for separation (crystallization conditions, choice of templates, choice of membrane) and advice on analytical and characterization methods.

Discovery of Novel Naphthoquinone-Chalcone Hybrids as Potent FGFR1 Tyrosine Kinase Inhibitors: Synthesis, Biological Evaluation, and Molecular Modeling

This work presents a flexible synthesis of 10 novel naphthoquinone-chalcone derivatives (1-10) by nucleophilic substitution of readily accessible aminochalcones and 2,3-dichloro-1,4-naphthoquinone. All compounds displayed broad-spectrum cytotoxic activities against all the tested cancer cell lines (i.e., HuCCA-1, HepG2, A549, MOLT-3, T47D, and MDA-MB-231) with IC50 values in the range of 0.81-62.06 μM, especially the four most potent compounds 1, 3, 8, and 9. The in vitro investigation on the fibroblast growth factor receptor 1 (FGFR1) inhibitory effect indicated that eight derivatives (1-2, 4-5, and 7-10) were active FGFR1 inhibitors (IC50 = 0.33-3.13 nM) with more potency than that of the known FGFR1 inhibitor, AZD4547 (IC50 = 12.17 nM). Promisingly, compounds 5 (IC50 = 0.33 ± 0.01 nM), 9 (IC50 = 0.50 ± 0.04 nM), and 7 (IC50 = 0.85 ± 0.08 nM) were the three most potent FGFR1 inhibitors. Molecular docking, molecular dynamics simulations, and MM/GBSA-based free energy calculation revealed that the key amino acid residues involved in the binding of the compounds 5, 7, and 9 and the target FGFR1 protein were similar with those of the AZD4547 (i.e., Val492, Lys514, Ile545, Val561, Ala640, and Asp641). These findings revealed that the newly synthesized naphthoquinone-chalcone scaffold is a promising structural feature for an efficient inhibition of FGFR1.

Discovery of amphotericin B, an antifungal drug as tyrosinase inhibitor with potent anti-melanogenic activity

Tyrosinase, a rate-limiting enzyme for melanin production, has been the most efficient target for the development of depigmenting agents. Although hydroquinone, kojic acid, and arbutin are the most well-known tyrosinase inhibitors, their adverse effects are inevitable. In the present study, an in silico drug repositioning combined with experimental validation was performed to search for novel potent tyrosinase inhibitors. Docking-based virtual screening results revealed that, among the 3210 FDA-approved drugs available in the ZINC database, amphotericin B, an antifungal drug exhibited the highest binding efficiency against human tyrosinase. Results from tyrosinase inhibition assay demonstrated that amphotericin B could inhibit the activity of mushroom and cellular tyrosinases, especially from MNT-1 human melanoma cells. Molecular modeling results revealed that amphotericin B/human tyrosinase complex exhibited high stability in an aqueous environment. Melanin assay results demonstrated that amphotericin B significantly suppressed melanin production in α-MSH-induced B16F10 murine melanoma and MNT-1 human melanoma cell lines better than the known inhibitor, kojic acid. Mechanistically, amphotericin B treatment significantly activated ERK and Akt signaling pathways, resulting in the decreased expression of MITF and tyrosinase. The obtained results may pursue pre-clinical and clinical studies to examine the possibility of using amphotericin B as an alternative treatment for hyperpigmentation disorders.

Raman Spectroscopic Analysis of Highly-Concentrated Antibodies under the Acid-Treated Conditions

PURPOSE

Antibody drugs are usually formulated as highly-concentrated solutions, which would easily generate aggregates, resulting in loss of efficacy. Although low pH increases the colloidal dispersion of antibodies, acid denaturation can be an issue. Therefore, knowing the physical properties at low pH under high concentration conditions is important.

METHODS

Raman spectroscopy was used to investigate pH-induced conformational changes of antibodies at 50 mg/ml. Experiments in pH 3 to 7 were performed for human serum IgG and recombinant rituximab.

RESULTS

We detected the evident changes at pH 3 in Tyr and Trp bands, which are the sensitive markers of intermolecular interactions. Thermal transition analysis over the pH range demonstrated that the thermal transition temperature (Tm) was highest at pH 3. Acid-treated and neutralized one showed higher Tm than that of pH 7, indicating that their extent of intermolecular interactions correlated with the Tm values. Onset temperature was clearly different between concentrated and diluted samples. Colloidal analyses confirmed the findings of the Raman analysis.

CONCLUSION

Our studies demonstrated the positive correlation between Raman analysis and colloidal information, validating as a method for evaluating antibody conformation associated with aggregation propensities.

Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis

Herein, an Energy Decomposition Analysis (EDA) scheme extended to the framework of QM/MM calculations in the context of electrostatic embeddings (QM/MM-EDA) including atomic charges and dipoles is applied to assess the effect of the QM region size on the convergence of the different interaction energy components, namely, electrostatic, Pauli, and polarization, for cationic, anionic, and neutral systems interacting with a strong polar environment (water). Significant improvements are found when the bulk solvent environment is described by a MM potential in the EDA scheme as compared to pure QM calculations that neglect bulk solvation. The predominant electrostatic interaction requires sizable QM regions. The results reported here show that it is necessary to include a surprisingly large number of water molecules in the QM region to obtain converged values for this energy term, contrary to most cluster models often employed in the literature. Both the improvement of the QM wave function by means of a larger basis set and the introduction of polarization into the MM region through a polarizable force field do not translate to a faster convergence with the QM region size, but they lead to better results for the different interaction energy components. The results obtained in this work provide insight into the effect of each energy component on the convergence of the solute-solvent interaction energy with the QM region size. This information can be used to improve the MM FFs and embedding schemes employed in QM/MM calculations of solvated systems.

Structural dynamics and susceptibility of isobutylamido thiazolyl resorcinol (ThiamidolTM) against human and mushroom tyrosinases

Tyrosinase, a key enzyme catalyzing a rate-limiting step of the melanin production, has been the most promising target for suppressing hyperpigmentation. Although a number of tyrosinase inhibitors have been developed, most of those lack clinical efficacy as they were identified from using mushroom tyrosinase (mTyr) as the target. Previous study revealed that the inhibitory effect of isobutylamido thiazolyl resorcinol (ThiamidolTM) on human tyrosinase (hTyr) is ∼100 times higher than that on mTyr. In the present study, we aimed to investigate the structural dynamics and susceptibility of ThiamidolTM against hTyr and mTyr at the atomic level using molecular docking, molecular dynamics simulation, and free energy calculation based on the molecular mechanics/Poisson-Boltzmann surface area method. The obtained results revealed that the resorcinol moiety of ThiamidolTM was found to be embedded in the catalytic copper center, interacting with H180, H202, H211, F386, and H390 residues of hTyr as well as with F264 residue of mTyr, mostly through van der Waals interactions. However, the number of destabilizing residues was found to be more pronounced in the ThiamidolTM/mTyr complex than the ThiamidolTM/hTyr system, supported by the lower binding affinity of ThiamidolTM/mTyr complex as well as the higher water accessibility and the lower number of atomic contacts at the active site of mTyr. Altogether, the structural and energetic information from this work would be useful for further optimization of more potent human tyrosinase inhibitors based on ThiamidolTM scaffold.Communicated by Ramaswamy H. Sarma.

Raman Marker Bands for Secondary Structure Changes of Frozen Therapeutic Monoclonal Antibody Formulations During Thawing

In this work we use Raman spectroscopy for protein characterization in the frozen state. We investigate the behavior of frozen therapeutic monoclonal antibody IgG1 formulation upon thawing by Raman spectroscopy. Secondary and tertiary structure of the protein in three different mab formulations in the frozen state are followed through observation of marker bands for α-helix, β-sheet and random coil. We identify the tyrosine intensity ratio I₈₅₆/I₈₃₀ as a marker for mab aggregation. Upon fast cooling (40 °C/min) to -80 °C we observe a significant increase of random coil and α -helical structures, while this is not the case for slower cooling (20 °C/min) to -80 °C. Most changes in the protein's secondary structure are observed in the course of thawing in the range up to -20 °C, when passing through the glass transitions and cold-crystallization of the two types of freeze-concentrated solutions formed through macro- and microcryoconcentration. An increase of protein concentration and the addition of mannitol suppress secondary structural changes but do no impact on aggregation.

Synergistic Effects of Multiple Excipients on Controlling Viscosity of Concentrated Protein Dispersions

Viscosity control is essential for the manufacturing and delivery of concentrated therapeutic proteins. Limited availability of the precious protein-based drugs hinders the characterization and screening of the formulation conditions with new types or different combinations of excipients. In this work, a droplet-based microfluidic device with incorporated multiple particle tracking microrheology (MPT) is developed to quantify the effects of two excipients, arginine hydrochloride (ArgHCl) and caffeine, on the viscosity of concentrated bovine gamma globulin (BGG) dispersions at two different values of pH. The effectiveness of both ArgHCl and caffeine show dependence on the BGG concentration and solution pH. The data set with high compositional resolution provides useful information to guide formulation with multiple viscosity-reducing excipients and quantification appropriate to start elucidating the connection to protein-protein interaction mechanisms. Overall, this work has demonstrated that the developed microfluidic approach has the potential to effectively assess the impact of multiple excipients on the viscosity and provide data for computational methods to predict viscosity for high concentration protein formulations.

Discovery of Anilino-1,4-naphthoquinones as Potent EGFR Tyrosine Kinase Inhibitors: Synthesis, Biological Evaluation, and Comprehensive Molecular Modeling

Epidermal growth factor receptor (EGFR) has been recognized as one of the attractive targets for anticancer drug development. Herein, a set of anilino-1,4-naphthoquinone derivatives (3-18) was synthesized and investigated for their anticancer and EGFR inhibitory potentials. Among all tested compounds, three derivatives (3, 8, and 10) were selected for studying EGFR inhibitory activity (in vitro and in silico) due to their most potent cytotoxic activities against six tested cancer cell lines (i.e., HuCCA-1, HepG2, A549, MOLT-3, MDA-MB-231, and T47D; IC50 values = 1.75-27.91 μM), high selectivity index (>20), and good predicted drug-like properties. The experimental results showed that these three promising compounds are potent EGFR inhibitors with nanomolar IC50 values (3.96-18.64 nM). Interestingly, the most potent compound 3 bearing 4-methyl substituent on the phenyl ring displayed 4-fold higher potency than the known EGFR inhibitor, erlotinib. Molecular docking, molecular dynamics simulation, and MM/GBSA-based free energy calculation revealed that van der Waals force played a major role in the accommodations of compound 3 within the ATP-binding pocket of EGFR. Additionally, the 4-CH3 moiety of the compound was noted to be a key chemical feature contributing to the highly potent EGFR inhibitory activity via its formations of alkyl interactions with A743, K745, M766, and L788 residues as well as additional interactions with M766 and T790.

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.

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.

Discovery of compounds with viscosity-reducing effects on biopharmaceutical formulations with monoclonal antibodies

•

Computational screening yielded 44 new viscosity-reducing agents on two model mAbs.

•

Dual excipients for viscosity reduction and solution buffering were discovered.

•

Compounds with three or more charges reduce the viscosity of model mAb formulations.

•

Filtering based on physicochemical properties can be applied to other mAb formulations.

For the development of concentrated monoclonal antibody formulations for subcutaneous administration, the main challenge is the high viscosity of the solutions. To compensate for this, viscosity reducing agents are commonly used as excipients. Here, we applied two computational chemistry approaches to discover new viscosity-reducing agents: fingerprint similarity searching, and physicochemical property filtering. In total, 94 compounds were selected and experimentally evaluated on two model monoclonal antibodies, which led to the discovery of 44 new viscosity-reducing agents. Analysis of the results showed that using a simple filter that selects only compounds with three or more charge groups is a good ‘rule of thumb’ for selecting potential viscosity-reducing agents for two model monoclonal antibody formulations.

Next Generation Strategy for Tuning the Thermoresponsive Properties of Micellar and Hydrogel Drug Delivery Vehicles Using Ionic Liquids

Amongst the greatest challenges in developing injectable controlled thermoresponsive micellar and hydrogel drug delivery vehicles include tuning the cloud point (CP) and reducing the gelation temperature (Tgel), below 37 °C,...

Antibody coated conductive polymer for the electrochemical immunosensing of Human Cardiac Troponin I in blood plasma

Cardiac troponin I (cTnI) is a sensitive biomarker for cardiovascular disease (CVD). Rapid determination of cTnI concentration in blood can greatly reduce the potential of significant heart damage and heart failure. Herein, we demonstrate a new electrochemical immunosensor for selective affinity binding and rapid detection of cTnI in blood plasma by an electrochemical method. A conductive film of "poly 2,5-bis(2-thienyl)3,4-diamine-terthiophene (PDATT)" was deposited onto an Indium Tin Oxide (ITO) electrode using chronoamperometry. Anti-cardiac troponin I antibody was then attached to the two amine (NH₂) groups substituted on the central thiophene of terthiophene repeating unit of the polymer chain via amide bond formation. The gaps on the surface of the antibody coated immunosensor were backfilled with bovine serum albumin (BSA) to prevent nonspecific binding of interfering molecules. Differential pulse voltammetry (DPV) was used to determine cTnI upon the formation of cTnI immunocomplex on the sensing surface, appearing a peak at 0.27 V. The response range was 0.01-100 ng mL^-1 with limit of quantification down to 0.01 ng mL^-1. The developed immunosensor was used to determine cTnI in spiked blood plasma without interference from cardiac troponin T (cTnT). Therefore, this new sensor can be utilised for the detection of cTnI biomarker in pathological laboratories and points of care in less than 15 min.

Towards an improved prediction of concentrated antibody solution viscosity using the Huggins coefficient

The viscosity of a monoclonal antibody solution must be monitored and controlled as it can adversely affect product processing, packaging and administration. Engineering low viscosity mAb formulations is challenging as prohibitive amounts of material are required for concentrated solution analysis, and it is difficult to predict viscosity from parameters obtained through low-volume, high-throughput measurements such as the interaction parameter, k_D, and the second osmotic virial coefficient, B₂₂. As a measure encompassing the effect of intermolecular interactions on dilute solution viscosity, the Huggins coefficient, k_h, is a promising candidate as a parameter measureable at low concentrations, but indicative of concentrated solution viscosity. In this study, a differential viscometry technique is developed to measure the intrinsic viscosity, [η], and the Huggins coefficient, k_h, of protein solutions. To understand the effect of colloidal protein-protein interactions on the viscosity of concentrated protein formulations, the viscometric parameters are compared to k_D and B₂₂ of two mAbs, tuning the contributions of repulsive and attractive forces to the net protein-protein interaction by adjusting solution pH and ionic strength. We find a strong correlation between the concentrated protein solution viscosity and the k_h but this was not observed for the k_D or the b₂₂, which have been previously used as indicators of high concentration viscosity. Trends observed in [η] and k_h values as a function of pH and ionic strength are rationalised in terms of protein-protein interactions.

Insights into the Conformation and Self-Association of a Concentrated Monoclonal Antibody using Isothermal Chemical Denaturation and Nuclear Magnetic Resonance

The purpose of this investigation was to highlight the utility of nuclear magnetic resonance (NMR) as a multi-attribute method for the characterization of therapeutic antibodies. In this case study, we compared results from isothermal chemical denaturation (ICD) and NMR with standard methods to relate conformational states of a model monoclonal antibody (mAb1) with protein-protein interactions (PPI) that lead to self - association in concentrated solutions. The increase in aggregation rate and relative viscosity for mAb1 was found to be both concentration and pH dependent. The free energy of unfolding (∆G⁰) from ICD and thermal analysis in dilute solutions indicated that although the native state predominated between pH 4 - pH 7, it was disrupted at the CH₂ and unfolded noncooperatively under acidic conditions. One-dimensional (1D) ¹H NMR and two-dimensional (2D) ¹³C-¹H NMR performed, in concentrated solutions, confirmed that PPI between pH 4-7 occurred while mAb1 was in the native state. NMR corroborated that mAb1 maintained a dominant native state at formulation-relevant conditions at the tested pH range, had increased global molecular tumbling dynamics at lower pH and confirmed increased PPI at higher pH conditions. This report aligns and compares typical characterization of an IgG1 with assessment of structure by NMR and provided a more precise assessment and deeper insight into the conformation of an IgG1 in concentrated solutions.

Self-assembly of stimuli-responsive coiled-coil fibrous hydrogels

Owing to their tunable properties, hydrogels comprised of stimuli-sensitive polymers are one of the most appealing scaffolds with applications in tissue engineering, drug delivery and other biomedical fields. We previously reported a thermoresponsive hydrogel formed using a coiled-coil protein, Q. Here, we expand our studies to identify the gelation of Q protein at distinct pH conditions, creating a protein hydrogel system that is sensitive to temperature and pH. Through secondary structure analysis, transmission electron microscopy, and rheology, we observed that Q self-assembles and forms fiber-based hydrogels exhibiting upper critical solution temperature behavior with increased elastic properties at pH 7.4 and pH 10. At pH 6, however, Q forms polydisperse nanoparticles, which do not further self-assemble and undergo gelation. The high net positive charge of Q at pH 6 creates significant electrostatic repulsion, preventing its gelation. This study will potentially guide the development of novel scaffolds and functional biomaterials that are sensitive towards biologically relevant stimuli.

Ultradilute Measurements of Self-Association for the Identification of Antibodies with Favorable High-Concentration Solution Properties

There is significant interest in formulating antibody therapeutics as concentrated liquid solutions, but early identification of developable antibodies with optimal manufacturability, stability, and delivery attributes remains challenging. Traditional methods of identifying developable mAbs with low self-association in common antibody formulations require relatively concentrated protein solutions (>1 mg/mL), and this single challenge has frustrated early-stage and large-scale identification of antibody candidates with drug-like colloidal properties. Here, we describe charge-stabilized self-interaction nanoparticle spectroscopy (CS-SINS), an affinity-capture nanoparticle assay that measures colloidal self-interactions at ultradilute antibody concentrations (0.01 mg/mL), and is predictive of antibody developability issues of high viscosity and opalescence that manifest at four orders of magnitude higher concentrations (>100 mg/mL). CS-SINS enables large-scale, high-throughput selection of developable antibodies during early discovery.

Identification of Vinyl Sulfone Derivatives as EGFR Tyrosine Kinase Inhibitor: In Vitro and In Silico Studies

Epidermal growth factor receptor (EGFR), overexpressed in many types of cancer, has been proved as a high potential target for targeted cancer therapy due to its role in regulating proliferation and survival of cancer cells. In the present study, a series of designed vinyl sulfone derivatives was screened against EGFR tyrosine kinase (EGFR-TK) using in silico and in vitro studies. The molecular docking results suggested that, among 78 vinyl sulfones, there were eight compounds that could interact well with the EGFR-TK at the ATP-binding site. Afterwards, these screened compounds were tested for the inhibitory activity towards EGFR-TK using ADP-Glo™ kinase assay, and we found that only VF16 compound exhibited promising inhibitory activity against EGFR-TK with the IC₅₀ value of 7.85 ± 0.88 nM. In addition, VF16 showed a high cytotoxicity with IC₅₀ values of 33.52 ± 2.57, 54.63 ± 0.09, and 30.38 ± 1.37 µM against the A431, A549, and H1975 cancer cell lines, respectively. From 500-ns MD simulation, the structural stability of VF16 in complex with EGFR-TK was quite stable, suggesting that this compound could be a novel small molecule inhibitor targeting EGFR-TK.

Application of Affinity-Capture Self-Interaction Nanoparticle Spectroscopy in Predicting Protein Stability, Especially for Co-Formulated Antibodies

PURPOSE

From traditional monoclonal antibodies to more and more complex mAb-based formulations, biopharmaceutical faces one challenge after another. To avoid these issues, identification of therapeutic proteins in the initial discovery process that has high stability and low self-interaction would simplify the development of safe and effective antibody therapeutics.

METHOD

Affinity-capture self-interaction nanoparticle spectroscopy (AC-SINS) is a new prediction method capable of identifying mAbs with different self-association propensity. In this study, 10 formulated monoclonal antibody (mAb) therapeutics include different mAb isotypes and co-formulated antibodies were measured by AC-SINS and some biophysical methods to predict protein stability. The prediction results of all 10 mAbs were compared to their stability data (Δ%monomer and Δ%HMWs) at accelerated (25°C and 40°C) and long-term storage conditions (4°C) as measured by size exclusion chromatography.

RESULT

AC-SINS method has a good predictive correlation with each mAbs and co-formulated antibodies. There were no physicochemical, intermolecular, or biological interactions that occurred between the two components of co-formulated antibodies which confirmed by Analytical ultracentrifugation (AUC).

CONCLUSION

Here we discuss the correlation between each method and protein stability, and also use AC-SINS assay to predict the stability of co-formulated antibodies for the first time. This may be an effective way to predict the stability of these complex mAb-based formulations such as co-formulated mAbs.

Comparison of Huggins Coefficients and Osmotic Second Virial Coefficients of Buffered Solutions of Monoclonal Antibodies

The Huggins coefficient kH is a well-known metric for quantifying the increase in solution viscosity arising from intermolecular interactions in relatively dilute macromolecular solutions, and there has been much interest in this solution property in connection with developing improved antibody therapeutics. While numerous kH measurements have been reported for select monoclonal antibodies (mAbs) solutions, there has been limited study of kH in terms of the fundamental molecular interactions that determine this property. In this paper, we compare measurements of the osmotic second virial coefficient B22, a common metric of intermolecular and interparticle interaction strength, to measurements of kH for model antibody solutions. This comparison is motivated by the seminal work of Russel for hard sphere particles having a short-range "sticky" interparticle interaction, and we also compare our data with known results for uncharged flexible polymers having variable excluded volume interactions because proteins are polypeptide chains. Our observations indicate that neither the adhesive hard sphere model, a common colloidal model of globular proteins, nor the familiar uncharged flexible polymer model, an excellent model of intrinsically disordered proteins, describes the dependence of kH of these antibodies on B22. Clearly, an improved understanding of protein and ion solvation by water as well as dipole-dipole and charge-dipole effects is required to understand the significance of kH from the standpoint of fundamental protein-protein interactions. Despite shortcomings in our theoretical understanding of kH for antibody solutions, this quantity provides a useful practical measure of the strength of interprotein interactions at elevated protein concentrations that is of direct significance for the development of antibody formulations that minimize the solution viscosity.

Dissecting the molecular basis of high viscosity of monospecific and bispecific IgG antibodies

Some antibodies exhibit elevated viscosity at high concentrations, making them poorly suited for therapeutic applications requiring administration by injection such as subcutaneous or ocular delivery. Here we studied an anti-IL-13/IL-17 bispecific IgG₄ antibody, which has anomalously high viscosity compared to its parent monospecific antibodies. The viscosity of the bispecific IgG₄ in solution was decreased by only ~30% in the presence of NaCl, suggesting electrostatic interactions are insufficient to fully explain the drivers of viscosity. Intriguingly, addition of arginine-HCl reduced the viscosity of the bispecific IgG₄ by ~50% to its parent IgG level. These data suggest that beyond electrostatics, additional types of interactions such as cation-π and/or π-π may contribute to high viscosity more significantly than previously understood. Molecular dynamics simulations of antibody fragments in the mixed solution of free arginine and explicit water were conducted to identify hotspots involved in self-interactions. Exposed surface aromatic amino acids displayed an increased number of contacts with arginine. Mutagenesis of the majority of aromatic residues pinpointed by molecular dynamics simulations effectively decreased the solution's viscosity when tested experimentally. This mutational method to reduce the viscosity of a bispecific antibody was extended to a monospecific anti-GCGR IgG₁ antibody with elevated viscosity. In all cases, point mutants were readily identified that both reduced viscosity and retained antigen-binding affinity. These studies demonstrate a new approach to mitigate high viscosity of some antibodies by mutagenesis of surface-exposed aromatic residues on complementarity-determining regions that may facilitate some clinical applications.

Product-Specific Impact of Viscosity Modulating Formulation Excipients During Ultra-High Concentration Biotherapeutics Drug Product Development

Developing ultra-high concentration biotherapeutics drug products can be challenging due to increased viscosity, processing, and stability issues. Excipients used to alleviate these concerns are traditionally evaluated at lower protein concentrations. This study investigates whether classically known modulators of stability and viscosity at low (<50 mg/mL) to high (>50 - 150 mg/mL) protein concentrations are beneficial in ultra-high (>150 mg/mL) concentration protein formulations and drug products. This study evaluates the effect of arginine monohydrochloride, proline, and lysine monohydrochloride on viscosity and concentratability at different high and ultra-high protein concentrations using a monoclonal antibody, mAbN, formulation as a candidate protein system. The effect of excipients on the viscosity and concentratability (rate and extent) was different at high versus ultra-high protein concentrations. These results highlight that classical excipients in literature known to modulate protein interactions at low protein concentrations and reduce viscosity at high protein concentrations may need to be evaluated at target protein concentrations in a product-specific manner while developing ultra-high concentration biologics drug products.

On the Aggregation and Nucleation Mechanism of the Monoclonal Antibody Anti-CD20 Near Liquid-Liquid Phase Separation (LLPS)

The crystallization of Anti-CD20, a full-length monoclonal antibody, has been studied in the PEG400/Na2SO4/Water system near Liquid-Liquid Phase Separation (LLPS) conditions by both sitting-drop vapour diffusion and batch methods. In order to understand the Anti-CD20 crystallization propensity in the solvent system of different compositions, we investigated some measurable parameters, normally used to assess protein conformational and colloidal stability in solution, with the aim to understand the aggregation mechanism of this complex biomacromolecule. We propose that under crystallization conditions a minor population of specifically aggregated protein molecules are present. While this minor species hardly contributes to the measured average solution behaviour, it induces and promotes crystal formation. The existence of this minor species is the result of the LLPS occurring concomitantly under crystallization conditions.

Multiscale Coarse-Grained Approach to Investigate Self-Association of Antibodies

Self-association of therapeutic monoclonal antibodies (mabs) are thought to modulate the undesirably high viscosity observed in their concentrated solutions. Computational prediction of such a self-association behavior is advantageous early during mab drug candidate selection when material availability is limited. Here, we present a coarse-grained (CG) simulation method that enables microsecond molecular dynamics simulations of full-length antibodies at high concentrations. The proposed approach differs from others in two ways: first, charges are assigned to CG beads in an effort to reproduce molecular multipole moments and charge asymmetry of full-length antibodies instead of only localized charges. This leads to great improvements in the agreement between CG and all-atom electrostatic fields. Second, the distinctive hydrophobic character of each antibody is incorporated through empirical adjustments to the short-range van der Waals terms dictated by cosolvent all-atom molecular dynamics simulations of antibody variable regions. CG simulations performed on a set of 15 different mabs reveal that diffusion coefficients in crowded environments are markedly impacted by intermolecular interactions. Diffusion coefficients computed from the simulations are in correlation with experimentally measured observables, including viscosities at a high concentration. Further, we show that the evaluation of electrostatic and hydrophobic characters of the mabs is useful in predicting the nonuniform effect of salt on the viscosity of mab solutions. This CG modeling approach is particularly applicable as a material-free screening tool for selecting antibody candidates with desirable viscosity properties.

Coarse-Grained Molecular Dynamics Simulations for Understanding the Impact of Short-Range Anisotropic Attractions on Structure and Viscosity of Concentrated Monoclonal Antibody Solutions

Understanding protein-protein interactions in concentrated therapeutic monoclonal antibody (mAb) solutions is desirable for improved drug discovery, processing, and administration. Here, we deduce both the net protein charge and the magnitude and geometry of short-ranged, anisotropic attractions of a mAb across multiple concentrations and cosolute conditions by comparing structure factors S(q) obtained from small-angle X-ray scattering experiments with those from molecular dynamics (MD) simulations. The simulations, which utilize coarse-grained 12-bead models exhibiting a uniform van der Waals attraction, uniform electrostatic repulsion, and short-range attractions between specific beads, are versatile enough to fit S(q) of a wide range of protein concentrations and ionic strength with the same charge on each bead and a single anisotropic short-range attraction strength. Cluster size distributions (CSDs) obtained from best fit simulations reveal that the experimental structure is consistent with small reversible oligomers in even low viscosity systems and help quantify the impact of these clusters on viscosity. The ability to systematically use experimental S(q) data together with MD simulations to discriminate between different possible protein-protein interactions, as well as to predict viscosities from protein CSDs, is beneficial for designing mAbs and developing formulation strategies that avoid high viscosities and aggregation at high concentration.

An Empirical Quantitative Model Describing Simultaneously Temperature and Concentration Effects on Protein Solution Viscosity

The viscosity of high-concentration protein solutions can lead to a range of challenges in drug product manufacturing and administration. Accurately modeling the viscosity of biologics solutions in response to changes in the formulation and surrounding environment is of significant interest and remains a challenge. Here, we show a practical method of modeling the viscosity of a therapeutic solution in response to changes in temperature and protein concentration. Our viscosity model consists of a Ross-Minton model of concentration dependence and a modified Arrhenius temperature dependence. We measured the viscosity as a function of concentration and temperature of 4 therapeutic antibodies in a range of potential clinical formulations. With these data, our model shows surprising generality, proving effective with different types of antibodies, formulations, and a range of more than 2 orders of magnitude in viscosity. Our approach is built on existing theory but provides a practical approach to modeling the viscosity of formulated drug product over the range of process-relevant concentrations and temperatures to better mitigate challenges in the drug manufacturing process.

Suck-Back Impact on Fluid Behavior at Filling Needle Tip

Needle clogging induces several issues during the filling step of injectable drugs, which makes essential to avoiding it to ensure a favorable outcome for the process. The suck-back function, present in peristaltic pumps, is often used empirically to that end. This study aims at describing and understanding the fluid behavior after suck-back application, which provides some quantitative specifications to prevent needle clogging.

Protein aggregation and mitigation strategy in low pH viral inactivation for monoclonal antibody purification

Significant amounts of soluble product aggregates were observed during low-pH viral inactivation (VI) scale-up for an IgG4 monoclonal antibody (mAb IgG4-N1), while small-scale experiments in the same condition showed negligible aggregation. Poor mixing and product exposure to low pH were identified as the root cause. To gain a mechanistic understanding of the problem, protein aggregation properties were studied by varying critical parameters including pH, hold time and protein concentration. Comprehensive biophysical characterization of product monomers and aggregates was performed using fluorescence-size-exclusion chromatography, differential scanning fluorimetry, fluorescence spectroscopy, and dynamic light scattering. Results showed IgG4-N1 partially unfolds at about pH 3.3 where the product molecules still exist largely as monomers owing to strong inter-molecular repulsions and favorable colloidal stability. In the subsequent neutralization step, however, the conformationally changed monomers are prone to aggregation due to weaker inter-molecular repulsions following the pH transition from 3.3 to 5.5. Surface charge calculations using homology modeling suggested that intra-molecular repulsions, especially between CH2 domains, may contribute to the IgG4-N1 unfolding at ≤ pH 3.3. Computational fluid dynamics (CFD) modeling was employed to simulate the conditions of pH titration to reduce the risk of aggregate formation. The low-pH zones during acid addition were characterized using CFD modeling and correlated to the condition causing severe product aggregation. The CFD tool integrated with the mAb solution properties was used to optimize the VI operating parameters for successful scale-up demonstration. Our research revealed the governing aggregation mechanism for IgG4-N1 under acidic conditions by linking its molecular properties and various process-related parameters to macroscopic aggregation phenomena. This study also provides useful insights into the cause and mitigation of low-pH-induced IgG4 aggregation in downstream VI operation.

Self-diffusion of a highly concentrated monoclonal antibody by fluorescence correlation spectroscopy: insight into protein-protein interactions and self-association

The dynamic behavior of monoclonal antibodies (mAbs) at high concentration provides insight into protein microstructure and protein-protein interactions (PPI) that influence solution viscosity and protein stability. At high concentration, interpretation of the collective-diffusion coefficient Dc, as determined by dynamic light scattering (DLS), is highly challenging given the complex hydrodynamics and PPI at close spacings. In contrast, self-diffusion of a tracer particle by Brownian motion is simpler to understand. Herein, we develop fluorescence correlation spectroscopy (FCS) for the measurement of the long-time self-diffusion of mAb2 over a wide range of concentrations and viscosities in multiple co-solute formulations with varying PPI. The normalized self-diffusion coefficient D0/Ds (equal to the microscopic relative viscosity ηeff/η0) was found to be smaller than η/η0. Smaller ratios of the microscopic to macroscopic viscosity (ηeff/η) are attributed to a combination of weaker PPI and less self-association. The interaction parameters extracted from fits of D0/Ds with a length scale dependent viscosity model agree with previous measurements of PPI by SLS and SAXS. Trends in the degree of self-association, estimated from ηeff/η with a microviscosity model, are consistent with oligomer sizes measured by SLS. Finally, measurements of collective diffusion and osmotic compressibility were combined with FCS data to demonstrate that the changes in self-diffusion between formulations are due primarily to changes in the protein-protein friction in these systems, and not to protein-solvent friction. Thus, FCS is a robust and accessible technique for measuring mAb self-diffusion, and, by extension, microviscosity, PPI and self-association that govern mAb solution dynamics.

Recovery, overloading, and protein interactions in asymmetrical flow field-flow fractionation

In asymmetrical flow field-flow fractionation (AF4), similar to other separation techniques, mass recovery and overloading require special attention in order to obtain quantitative results. We conducted a systematic study with five globular proteins of different molecular weight (36.7–669 kDa) and isoelectric point (4.0–6.5), and ultrafiltration membranes that are commonly used in aqueous AF4, regenerated cellulose (RC) and polyethersulfone (PES). Phosphate-buffered saline (PBS) with ionic strength 0.15 M and pH 7.2 was used as the carrier liquid in this study. The actual molecular weight cutoff (MWCO) was found to be higher than the nominal value and varied between membranes of different chemistry but the same nominal MWCO. Adsorption on the membrane was found to be dependent on the membrane chemistry (RC had lower adsorption compared to PES), and independent of the protein standard for the examined proteins. On the other hand, the mass overloading effects (i.e., higher retention times, peak broadening, and fronting peaks) were significantly more pronounced for γ-globulin than for the other proteins. The overloading effects could be rationalized with the increase of the local viscosity close to the membrane, depending on the properties of the proteins, and we derived theoretical equations that related the dependency of the migration velocity on the protein concentration through this non-ideal viscosity effect.

Web-based display of protein surface and pH-dependent properties for assessing the developability of biotherapeutics

Protein instability leads to reversible self-association and irreversible aggregation which is a major concern for developing new biopharmaceutical leads. Protein solution behaviour is dictated by the physicochemical properties of the protein and the solution. Optimising protein solutions through experimental screens and targeted protein engineering can be a difficult and time consuming process. Here, we describe development of the protein-sol web server, which was previously restricted to protein solubility prediction from amino acid sequence. Tools are presented for calculating and mapping patches of hydrophobicity and charge on the protein surface. In addition, predictions of folded state stability and net charge are displayed as a heatmap for a range of pH and ionic strength conditions. Tools are evaluated in the context of antibodies, their fragments and interactions. Surprisingly, antibody-antigen interfaces are, on average, at least as polar as Fab surfaces. This benchmarking process provides the user with thresholds with which to assess non-polar surface patches, and possible solubility implications, in proteins of interest. Stability heatmaps compare favourably with experimental data for CH2 and CH3 domains. Display and quantification of surface polarity and pH/ionic strength dependence will be useful generally for investigation of protein biophysics.

Correlating multidimensional short-term empirical protein properties to long-term protein physical stability data via empirical phase diagrams

Identification of long-term stable biopharmaceutical formulations is essential for biopharmaceutical product development. Reduction of the number of long-term storage experiments and a well-defined formulation search space requires knowledge-based formulation screenings and a detailed protein phase behavior understanding. To achieve this, short-term analytical techniques can serve as predictors for long-term protein phase behavior. Protein phase behavior studies that investigate this concept commonly display shortcomings such as limited and small datasets, sample adjustments, or simplistic data analysis. To overcome these shortcomings, 150 unique lysozyme solutions were analyzed using six different short-term analytical techniques. Lysozyme's structural properties, conformational stability, colloidal stability, surface charge, and surface hydrophobicity were obtained directly after formulation preparation. Employing the empirical phase diagram method, this short-term data was correlated to long-term physical stability data obtained during 40 days of storage. Short-term protein properties showed partial correlation to long-term phase behavior. Structural differences, changing surface properties, colloidal stability, and conformation stability as a function of formulation conditions were observed. This study contributes to long-term protein phase behavior research by presenting a systematic, data-dependent, and multidimensional data evaluation workflow to create a comprehensive overview of short-term protein analytics in relation to long-term protein phase behavior.

Protein-Protein Interactions of Highly Concentrated Monoclonal Antibody Solutions via Static Light Scattering and Influence on the Viscosity

The ability to design and formulate mAbs to minimize attractive interactions at high concentrations is important for protein processing, stability, and administration, particularly in subcutaneous delivery, where high viscosities are often challenging. The strength of protein-protein interactions (PPIs) of an IgG1 and IgG4 monoclonal antibody (mAb) from low to high concentration was determined by static light scattering (SLS) and used to understand viscosity data. The PPI were tuned using NaCl and five organic ionic co-solutes. The PPI strength was quantified by the normalized structure factor S(0)/ S(0)_HS and Kirkwood-Buff integral G₂₂/ G_22,HS (HS = hard sphere) determined from the SLS data and also by fits with (1) a spherical Yukawa potential and (2) an interacting hard sphere (IHS) model, which describes attraction in terms of hypothetical oligomers. The IHS model was better able to capture the scattering behavior of the more strongly interacting systems (mAb and/or co-solute) than the spherical Yukawa potential. For each descriptor of PPI, linear correlations were obtained between the viscosity at high concentration (200 mg/mL) and the interaction strengths evaluated both at low (20 mg/mL) and high concentrations (200 mg/mL) for a given mAb. However, the only parameter that provided a correlation across both mAbs was the oligomer mass ratio ( m_oligomer/ m_{monomer+dimer}) from the IHS model, indicating the importance of self-association (in addition to the direct influence of the attractive PPI) on the viscosity.

Spectroscopic methods for assessing the molecular origins of macroscopic solution properties of highly concentrated liquid protein solutions

In cases of subcutaneous injection of therapeutic monoclonal antibodies, high protein concentrations (>50 mg/ml) are often required. During the development of these high concentration liquid formulations (HCLF), challenges such as aggregation, gelation, opalescence, phase separation, and high solution viscosities are more prone compared to low concentrated protein formulations. These properties can impair manufacturing processes, as well as protein stability and shelf life. To avoid such unfavourable solution properties, a detailed understanding about the nature of these properties and their driving forces are required. However, the fundamental mechanisms that lead to macroscopic solution properties, as above mentioned, are complex and not fully understood, yet. Established analytical methods for assessing the colloidal stability, i.e. the ability of a native protein to remain dispersed in solution, are restricted to dilute conditions and provide parameters such as the second osmotic virial coefficient, B₂₂, and the diffusion interaction coefficient, k_D. These parameters are routinely applied for qualitative estimations and identifications of proteins with challenging solution behaviours, such as high viscosities and aggregation, although the assays are prepared for low protein concentration conditions, typically between 0.1 and 20 mg/ml ("ideal" solution conditions). Quantitative analysis of samples of high protein concentration is difficult and it is hard to obtain information about the driving forces of such solution properties and corresponding protein-protein self-interactions. An advantage of using specific spectroscopic methods is the potential of directly analysing highly concentrated protein solutions at different solution conditions. This allows for collecting/gaining valuable information about the fundamental mechanisms of solution properties of the high protein concentration regime. In addition, the derived parameters might be more predictive as compared to the parameters originating from assays which are optimized for the low protein concentration range. The provided information includes structural data, molecular dynamics at various timescales and protein-solvent interactions, which can be obtained at molecular resolution. Herein, we provide an overview about spectroscopic techniques for analysing the origins of macroscopic solution behaviours in general, with a specific focus on pharmaceutically relevant high protein concentration and formulation conditions.