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

[Developability assessment]

The therapeutic antibodies and their by-products (antibody fragments, conjugated, etc.) establish one of the most dynamic biopharmaceutical market segments today. Due to their intrinsic properties of specificity towards their target, towards their flexible affinity and due to their stability, antibodies became therapeutic agents of the very first choice. One of the challenges of this sector is to create antibodies of very good quality, more and more quickly, while having less and less consequent development costs in fine.

Advancing Therapeutic Protein Discovery and Development through Comprehensive Computational and Biophysical Characterization

Therapeutic protein candidates should exhibit favorable properties that render them suitable to become drugs. Nevertheless, there are no well-established guidelines for the efficient selection of proteinaceous molecules with desired features during early stage development. Such guidelines can emerge only from a large body of published research that employs orthogonal techniques to characterize therapeutic proteins in different formulations. In this work, we share a study on a diverse group of proteins, including their primary sequences, purity data, and computational and biophysical characterization at different pH and ionic strength. We report weak linear correlations between many of the biophysical parameters. We suggest that a stability comparison of diverse therapeutic protein candidates should be based on a computational and biophysical characterization in multiple formulation conditions, as the latter can largely determine whether a protein is above or below a certain stability threshold. We use the presented data set to calculate several stability risk scores obtained with an increasing level of analytical effort and show how they correlate with protein aggregation during storage. Our work highlights the importance of developing combined risk scores that can be used for early stage developability assessment. We suggest that such scores can have high prediction accuracy only when they are based on protein stability characterization in different solution conditions.

Potential and limits of a colloid approach to protein solutions

Looking at globular proteins with the eyes of a colloid scientist has a long tradition, in fact a significant part of the early colloid literature was focused on protein solutions. However, it has also been recognized that proteins are much more complex than the typical hard sphere-like synthetic model colloids. Proteins are not perfect spheres, their interaction potentials are in general not isotropic, and using theories developed for such particles are thus clearly inadequate in many cases. In this perspective article, we now take a closer look at the field. In particular, we reflect on the fact that modern colloid science has been undergoing a tremendous development, where a multitude of novel systems have been developed in the lab and in silico. During the last decade we have seen a rapidly increasing number of reports on the synthesis of anisotropic, patchy and/or responsive synthetic colloids, that start to resemble their complex biological counterparts. This experimental development is also reflected in a corresponding theoretical and simulation effort. The experimental and theoretical toolbox of colloid science has thus rapidly expanded, and there is obviously an enormous potential for an application of these new concepts to protein solutions, which has already been realized and harvested in recent years. In this perspective article we make an attempt to critically discuss the exploitation of colloid science concepts to better understand protein solutions. We not only consider classical applications such as the attempt to understand and predict solution stability and phase behaviour, but also discuss new challenges related to the dynamics, flow behaviour and liquid-solid transitions found in concentrated or crowded protein solutions. It not only aims to provide an overview on the progress in experimental and theoretical (bio)colloid science, but also discusses current shortcomings in our ability to correctly reproduce and predict the structural and dynamic properties of protein solutions based on such a colloid approach.

Binding symmetry and surface flexibility mediate antibody self-association

Solution stability is an important factor in the optimization of engineered biotherapeutic candidates such as monoclonal antibodies because of its possible effects on manufacturability, pharmacology, efficacy and safety. A detailed atomic understanding of the mechanisms governing self-association of natively folded protein monomers is required to devise predictive tools to guide screening and re-engineering along the drug development pipeline. We investigated pairs of affinity-matured full-size antibodies and observed drastically different propensities to aggregate from variants differing by a single amino-acid. Biophysical testing showed that antigen-binding fragments (Fabs) from the aggregating antibodies also reversibly associated with equilibrium dissociation constants in the low-micromolar range. Crystal structures (PDB accession codes 6MXR, 6MXS, 6MY4, 6MY5) and bottom-up hydrogen-exchange mass spectrometry revealed that Fab self-association occurs in a symmetric mode that involves the antigen complementarity-determining regions. Subtle local conformational changes incurred upon point mutation of monomeric variants foster formation of complementary polar interactions and hydrophobic contacts to generate a dimeric Fab interface. Testing of popular in silico tools generally indicated low reliabilities for predicting the aggregation propensities observed. A structure-aggregation data set is provided here in order to stimulate further improvements of in silico tools for prediction of native aggregation. Incorporation of intermolecular docking, conformational flexibility, and short-range packing interactions may all be necessary features of the ideal algorithm.

Molecular Computations of Preferential Interaction Coefficients of IgG1 Monoclonal Antibodies with Sorbitol, Sucrose, and Trehalose and the Impact of These Excipients on Aggregation and Viscosity

Preferential interactions of formulation excipients govern their overall interactions with protein molecules, and molecular dynamics simulations allow for the examination of the interactions at the molecular level. We used molecular dynamics simulations to examine the interactions of sorbitol, sucrose, and trehalose with three different IgG1 antibodies to gain insight into how these excipients impact aggregation and viscosity. We found that sucrose and trehalose reduce aggregation more than sorbitol because of their larger size and their stronger interactions with high-spatial aggregation propensity residues compared to sorbitol. Two of the antibodies had high viscosity in sodium acetate buffer, and for these, we found that sucrose and trehalose tended to have opposite effects on viscosity. The data presented here provide further insight into the mechanisms of interactions of these three carbohydrate excipients with the antibody surface and thus their impact on excipient stabilization of antibody formulations.

Science and art of protein formulation development

Protein pharmaceuticals have become a significant class of marketed drug products and are expected to grow steadily over the next decade. Development of a commercial protein product is, however, a rather complex process. A critical step in this process is formulation development, enabling the final product configuration. A number of challenges still exist in the formulation development process. This review is intended to discuss these challenges, to illustrate the basic formulation development processes, and to compare the options and strategies in practical formulation development.

Concentration Effects in the Interaction of Monoclonal Antibodies (mAbs) with their Immediate Environment Characterized by EPR Spectroscopy

Monoclonal antibodies (mAbs) are often needed and applied in high concentration solutions, >100 mg/mL. Due to close intermolecular distances between mAbs at high concentrations (~10–20 nm at 200 mg/mL), intermolecular interactions between mAbs and mAbs and solvent/co-solute molecules become non-negligible. Here, EPR spectroscopy is used to study the high-concentration solutions of mAbs and their effect on co-solvated small molecules, using EPR “spin probing” assay in aqueous and buffered solutions. Such, information regarding the surrounding environments of mAbs at high concentrations were obtained and comparisons between EPR-obtained micro-viscosities (rotational correlation times) and macroscopic viscosities measured by rheology were possible. In comparison with highly viscous systems like glycerol-water mixtures, it was found that up to concentrations of 50 mg/mL, the mAb-spin probe systems have similar trends in their macro- (rheology) and micro-viscosities (EPR), whereas at very high concentrations they deviate strongly. The charged spin probes sense an almost unchanged aqueous solution even at very high concentrations, which in turn indicates the existence of large solvent regions that despite their proximity to large mAbs essentially offer pure water reservoirs for co-solvated charged molecules. In contrast, in buffered solutions, amphiphilic spin probes like TEMPO interact with the mAb network, due to slight charge screening. The application of EPR spectroscopy in the present work has enabled us to observe and discriminate between electrostatic and hydrophobic kinds of interactions and depict the potential underlying mechanisms of network formation at high concentrations of mAbs. These findings could be of importance as well for the development of liquid-liquid phase separations often observed in highly concentrated protein solutions.

IgG Charge: Practical and Biological Implications

Practically, IgG charge can contribute significantly to thermodynamic nonideality, and hence to solubility and viscosity. Biologically, IgG charge isomers exhibit differences in clearance and potency. It has been known since the 1930s that all immunoglobulins carry a weak negative charge in physiological solvents. However, there has been no systematic exploration of this fundamental property. Accurate charge measurements have been made using membrane confined electrophoresis in two solvents (pH 5.0 and pH 7.4) on a panel of twelve mAb IgGs, as well as their F(ab’)₂ and Fc fragments. The following observations were made at pH 5.0: (1) the measured charge differs from the calculated charge by ~40 for the intact IgGs, and by ~20 for the Fcs; (2) the intact IgG charge depends on both Fv and Fc sequences, but does not equal the sum of the F(ab)’₂ and Fc charge; (3) the Fc charge is consistent within a class. In phosphate buffered saline, pH 7.4: (1) the intact IgG charges ranged from 0 to −13; (2) the F(ab’)₂ fragments are nearly neutral for IgG1s and IgG2s, and about −5 for some of the IgG4s; (3) all Fc fragments are weakly anionic, with IgG1 < IgG2 < IgG4; (4) the charge on the intact IgGs does not equal the sum of the F(ab’)₂ and Fc charge. In no case is the calculated charge, based solely on H⁺ binding, remotely close to the measured charge. Some mAbs carried a charge in physiological salt that was outside the range observed for serum-purified human poly IgG. To best match physiological properties, a therapeutic mAb should have a measured charge that falls within the range observed for serum-derived human IgGs. A thermodynamically rigorous, concentration-dependent protein–protein interaction parameter is introduced. Based on readily measured properties, interaction curves may be generated to aid in the selection of proteins and solvent conditions. Example curves are provided.

Selecting and engineering monoclonal antibodies with drug-like specificity

Despite the recent explosion in the use of monoclonal antibodies (mAbs) as drugs, it remains a significant challenge to generate antibodies with a combination of physicochemical properties that are optimal for therapeutic applications. We argue that one of the most important and underappreciated drug-like antibody properties is high specificity - defined here as low levels of antibody non-specific and self-interactions - which is linked to low off-target binding and slow antibody clearance in vivo and high solubility and low viscosity in vitro. Here, we review the latest advances in characterizing antibody specificity and elucidating its molecular determinants as well as using these findings to improve the selection and engineering of antibodies with extremely high, drug-like specificity.

Impact of electroviscous effect on viscosity in developing highly concentrated protein formulations: Lessons from non-protein charged colloids

Subcutaneous delivery of highly concentrated protein formulations is paramount for reducing healthcare cost and improving patient compliance, where reducing the solution viscosity of formulations is critical for drug delivery. The objective of this paper is to provide some mechanistic understanding about the contribution of electrostatic repulsion to the viscosity of protein solutions at high concentrations, along with the effect of excipients such as salts on relative viscosity. Proteins are treated as charged colloids in this paper. At high concentrations, the electrical double layer starts to overlap, and secondary electroviscous effect becomes significant in addition to primary electroviscous effect. In other words, the hydrodynamic volume of proteins plays a great role in influencing their solution viscosity because of the excluded volume effect. Currently, it is hypothesized that the high viscosity of concentrated protein solutions is attributed to formation of clusters due to either electrostatic attraction or hydrophobic interactions, especially for monoclonal antibodies, in which anybody molecules in high concentration formulations may form networks. Consequently, viscosity reduction in the presence of inorganic or organic salts in these formulations is due to breaking up of these networks. In this review, authors hope to provide another point of view based on the effect of the electrostatic repulsion on the excluded volume-hydrodynamic volume. Finally, authors hope the proposed theoretical framework can be used to guide excipient selection in the product development of highly concentrated proteins.

Modeling the depletion effect caused by an addition of polymer to monoclonal antibody solutions

We present a theoretical study of colloidal stability of the model mixtures of monoclonal antibody molecules and non-adsorbing (no polymer-protein attraction) polymers. The antibodies are pictured as an assembly of seven hard spheres assuming a Y-like shape. Polymers present in the mixture are modeled as chain-like molecules having from 32 up to 128 monomers represented as hard spheres. We use Wertheim's thermodynamic perturbation theory to construct the two molecular species and to calculate measurable properties. The calculations are performed in the osmotic ensemble. In view that no direct attractive interaction is present in the model Hamiltonian, we only account for the entropic contribution to the phase equilibrium. We calculate chemical potentials and the equation of state for the model mixture to determine the liquid-liquid part of the phase diagram. We investigate how the critical antibody number density depends on the degree of polymerization and the bead size ratio of the polymer and protein components. The model mixture qualitatively correctly predicts some basic features of real systems. The effects of the model 'protein' geometry, that is the difference in results for the flexible Y-shaped protein versus the rigid spherical one, are also examined.

Controlling the viscosities of antibody solutions through control of their binding sites

For biotechnological drugs, it is desirable to formulate antibody solutions with low viscosities. We go beyond previous colloid theories in treating protein-protein self-association of molecules that are antibody-shaped and flexible and have spatially specific binding sites. We consider interactions either through fragment antigen (Fab-Fab) or fragment crystalizable (Fab-Fc) binding. Wertheim's theory is adapted to compute the cluster-size distributions, viscosities, second virial coefficients, and Huggins coefficients, as functions of antibody concentration. We find that the aggregation properties of concentrated solutions can be anticipated from simpler-to-measure dilute solutions. A principal finding is that aggregation is controllable, in principle, through modifying the antibody itself, and not just the solution it is dissolved in. In particular: (i) monospecific antibodies having two identical Fab arms can form linear chains with intermediate viscosities. (ii) Bispecific antibodies having different Fab arms can, in some cases, only dimerize, having low viscosities. (iii) Arm-to-Fc binding allows for three binding partners, leading to networks and high viscosities.

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.

Weak IgG self- and hetero-association characterized by fluorescence analytical ultracentrifugation

Weak protein-protein interactions may be important to binding cooperativity. A panel of seven fluorescently labeled tracer monoclonal IgG antibodies, differing in variable (V) and constant (C) region sequences, were sedimented in increasing concentrations of unlabeled IgGs of identical, similar, and different backgrounds. Weak IgG::IgG attractive interactions were detected and characterized by global analysis of the hydrodynamic nonideality coefficient, ks . The effects of salt concentration and temperature on ks suggest the interactions are predominantly enthalpic in origin. The interactions were found to be variable in strength, affected by both the variable and constant regions, but indiscriminate with respect to IgG subclass. Furthermore, weak attractive interactions were observed for all the mAbs with freshly purified human poly-IgG. The universality of the weak interactions suggest that they may contribute to effector function cooperativity in the normal immune response, and we postulate that the generality of the interactions allows for a broader range of epitope spacing for complement activation. These studies demonstrate the utility of analytical ultracentrifuge fluorescence detection in measuring weak protein-protein interactions. It also shows the strength of global analysis of sedimentation velocity data by SEDANAL to extract hydrodynamic nonideality ks to characterize weak macromolecular interactions.

Theory for the Liquid-Liquid Phase Separation in Aqueous Antibody Solutions

This study presents the theory for liquid-liquid phase separation for systems of molecules modeling monoclonal antibodies. Individual molecule is depicted as an assembly of seven hard spheres, organized to mimic the Y-shaped antibody. We consider the antibody-antibody interactions either through Fab, Fab' (two Fab fragments may be different), or Fc domain. Interaction between these three domains of the molecule (hereafter denoted as A, B, and C, respectively) is modeled by a short-range square-well attraction. To obtain numerical results for the model under study, we adapt Wertheim's thermodynamic perturbation theory. We use this model to calculate the liquid-liquid phase separation curve and the second virial coefficient B2. Various interaction scenarios are examined to see how the strength of the site-site interactions and their range shape the coexistence curve. In the asymmetric case, where an attraction between two sites is favored and the interaction energies for the other sites kept constant, critical temperature first increases and than strongly decreases. Some more microscopic information, for example, the probability for the particular two sites to be connected, has been calculated. Analysis of the experimental liquid-liquid phase diagrams, obtained from literature, is presented. In addition, we calculate the second virial coefficient under conditions leading to the liquid-liquid phase separation and present this quantity on the graph B2 versus protein concentration.

Predicting Protein-Protein Interactions of Concentrated Antibody Solutions Using Dilute Solution Data and Coarse-Grained Molecular Models

Protein-protein interactions for solutions of an IgG1 molecule were quantified using static light scattering (SLS) measurements from low to high protein concentrations (c₂). SLS was used to determine second osmotic virial coefficients (B₂₂) at low c₂, and excess Rayleigh profiles (R^ex/K vs. c₂) and zero-q structure factors (S_q=0) as a function of c₂ at higher c₂ for a series of conditions (pH, sucrose concentration, and total ionic strength [TIS]). Repulsive (attractive) interactions were observed at low TIS (high TIS) for pH 5 and 6.5, with increasing repulsions when 5% w/w sucrose was also present. Previously developed and refined coarse-grained antibody models were used to fit model parameters from B₂₂ versus TIS data. The resulting parameters from low-c₂ conditions were used as the sole input to multiprotein Monte Carlo simulations to predict high-c₂R^ex/K and S_q=0 behavior up to 150 g/L. Experimental results at high-c₂ conditions were quantitatively predicted by the simulations for the coarse-grained models that treated antibody molecules as either 6 or 12 (sub) domains, which preserved the basic shape of a monoclonal antibody. Finally, preferential accumulation of sucrose around the protein surface was identified via high-precision density measurements, which self-consistently explained the simulation and experimental SLS results.

Specific optical rotation is a versatile tool for the identification of critical micelle concentration and micellar growth of tartaric acid-based diastereomeric amphiphiles

Four novel tartaric acid-based diastereomeric chiral amphiphiles, two being enantiomers of the other two, have been synthesized and investigated using chiroptical spectroscopic methods, along with tensiometry and dynamic light scattering experiments. We found that an inflection point in specific optical rotation (SOR) values at ~0.32 mM corresponds to the critical micelle concentration (CMC). The increase in magnitude of SOR values beyond CMC corresponds to the growth of aggregates. For enantiomers, oppositely signed SOR values were observed, ruling out the possibility for the presence of aggregation size mediated artefacts. SOR values did not exhibit concentration dependence for a chiral tartaric acid based non-aggregating analogue further establishing the absence of artefacts or anomalous interaction of tartaric acid based head group with solvent. Electronic circular dichroism spectra showed no significant changes in band positions or intensities with concentration. Due to the requirement for higher concentrations (~200 mM) needed to obtain vibrational circular dichroism spectra, these measurements are not found to be useful for studying concentration dependent properties of chiral amphiphiles.

Rapid and accurate in silico solubility screening of a monoclonal antibody library

Antibodies represent essential tools in research and diagnostics and are rapidly growing in importance as therapeutics. Commonly used methods to obtain novel antibodies typically yield several candidates capable of engaging a given target. The development steps that follow, however, are usually performed with only one or few candidates since they can be resource demanding, thereby increasing the risk of failure of the overall antibody discovery program. In particular, insufficient solubility, which may lead to aggregation under typical storage conditions, often hinders the ability of a candidate antibody to be developed and manufactured. Here we show that the selection of soluble lead antibodies from an initial library screening can be greatly facilitated by a fast computational prediction of solubility that requires only the amino acid sequence as input. We quantitatively validate this approach on a panel of nine distinct monoclonal antibodies targeting nerve growth factor (NGF), for which we compare the predicted and measured solubilities finding a very close match, and we further benchmark our predictions with published experimental data on aggregation hotspots and solubility of mutational variants of one of these antibodies.

Preferential interactions of trehalose, L-arginine.HCl and sodium chloride with therapeutically relevant IgG1 monoclonal antibodies

Preferential interactions of weakly interacting formulation excipients govern their effect on the equilibrium and kinetics of several reactions of protein molecules in solution. Using vapor pressure osmometry, we characterized the preferential interactions of commonly used excipients trehalose, L-arginine.HCl and NaCl with three therapeutically-relevant, IgG1 monoclonal antibodies that have similar size and shape, but differ in their surface hydrophobicity and net charge. We further characterized the effect of these excipients on the reversible self-association, aggregation and viscosity behavior of these antibody molecules. We report that trehalose, L-arginine.HCl and NaCl are all excluded from the surface of the three IgG1 monoclonal antibodies, and that the exclusion behavior is linearly related to the excipient molality in the case of trehalose and NaCl, whereas a non-linear behavior is observed for L-arginine.HCl. Interestingly, we find that the magnitude of trehalose exclusion depends upon the nature of the protein surface. Such behavior is not observed in case of NaCl and L-arginine.HCl as they are excluded to the same extent from the surface of all three antibody molecules tested in this study. Analysis of data presented in this study provides further insight into the mechanisms governing excipient-mediated stabilization of mAb formulations.

Challenges in Predicting Protein-Protein Interactions from Measurements of Molecular Diffusivity

Dynamic light scattering can be used to measure the diffusivity of a protein within a formulation. The dependence of molecular diffusivity on protein concentration (traditionally expressed in terms of the interaction parameter k_D) is often used to infer whether protein-protein interactions are repulsive or attractive, resulting in solutions that are colloidally stable or unstable, respectively. However, a number of factors unrelated to intermolecular forces can also impact protein diffusion, complicating this interpretation. Here, we investigate the influence of multicomponent diffusion in a ternary protein-salt-water system on protein diffusion and k_D in the context of Nernst-Planck theory. This analysis demonstrates that large changes in protein diffusivity with protein concentration can result even for hard-sphere systems in the absence of protein-protein interactions. In addition, we show that dynamic light scattering measurements of diffusivity made at low ionic strength cannot be reliably used to detect protein conformational changes. We recommend comparing experimentally determined k_D values to theoretically predicted excluded-volume contributions, which will allow a more accurate assessment of protein-protein interactions.

Challenges for the pharmaceutical technical development of protein coformulations

Objectives

This review discusses challenges to stability, analytics and manufacturing of protein coformulations. Furthermore, general considerations to be taken into account for the pharmaceutical development of coformulated protein drug products are highlighted.

Key findings

Coformulation of two or more active substances in one single dosage form has recently seen increasing use offering several advantages, such as increased efficacy and/or the overall reduction of adverse event incidents in patients. Most marketed coformulated drug products are composed of small molecules. As proteins are not only comparatively large but also complex molecules, the maintenance of their physicochemical integrity within a formulation throughout pharmaceutical processing, storage, transport, handling and patient administration to ensure proper pharmacokinetics and pharmacodynamics in vivo already represents various challenges for single-entity products. Thus, nowadays, only sparse biologics-based coformulations can be found, as additional complexity during development is given for these products.

Summary

The complexity of the dosage form and the protein molecules results into additional challenges to formulation, manufacture, storage, transport, handling and patient administration, stability and analytics during the pharmaceutical development of protein coformulations. Various points have to be considered during different stages of development in order to obtain a safe and efficacious product.

Biopharmaceutical Informatics: supporting biologic drug development via molecular modelling and informatics

Objectives

The purpose of this article is to introduce an emerging field called ‘Biopharmaceutical Informatics’. It describes how tools from Information technology and Molecular Biophysics can be adapted, developed and gainfully employed in discovery and development of biologic drugs.

Key Findings

The findings described here are based on literature surveys and the authors’ collective experiences in the field of biologic drug product development. A strategic framework to forecast early the hurdles faced during drug product development is weaved together and elucidated using chemical degradation as an example. Efficiency of translating biologic drug discoveries into drug products can be significantly improved by combining learnings from experimental biophysical and analytical data on the drug candidates with molecular properties computed from their sequences and structures via molecular modeling and simulations.

Summary

Biopharmaceutical Informatics seeks to promote applications of computational tools towards discovery and development of biologic drugs. When fully implemented, industry-wide, it will enable rapid materials-free developability assessments of biologic drug candidates at early stages as well as streamline drug product development activities such as commercial scale production, purification, formulation, analytical characterization, safety and in vivo performance.

In-silico prediction of concentration-dependent viscosity curves for monoclonal antibody solutions

Early stage developability assessments of monoclonal antibody (mAb) candidates can help reduce risks and costs associated with their product development. Forecasting viscosity of highly concentrated mAb solutions is an important aspect of such developability assessments. Reliable predictions of concentration-dependent viscosity behaviors for mAb solutions in platform formulations can help screen or optimize drug candidates for flexible manufacturing and drug delivery options. Here, we present a computational method to predict concentration-dependent viscosity curves for mAbs solely from their sequence-structural attributes. This method was developed using experimental data on 16 different mAbs whose concentration-dependent viscosity curves were experimentally obtained under standardized conditions. Each concentration-dependent viscosity curve was fitted with a straight line, via logarithmic manipulations, and the values for intercept and slope were obtained. Intercept, which relates to antibody diffusivity, was found to be nearly constant. In contrast, slope, the rate of increase in solution viscosity with solute concentration, varied significantly across different mAbs, demonstrating the importance of intermolecular interactions toward viscosity. Next, several molecular descriptors for electrostatic and hydrophobic properties of the 16 mAbs derived using their full-length homology models were examined for potential correlations with the slope. An equation consisting of hydrophobic surface area of full-length antibody and charges on V_H, V_L, and hinge regions was found to be capable of predicting the concentration-dependent viscosity curves of the antibody solutions. Availability of this computational tool may facilitate material-free high-throughput screening of antibody candidates during early stages of drug discovery and development.

Mitigation of reversible self-association and viscosity in a human IgG1 monoclonal antibody by rational, structure-guided Fv engineering

Undesired solution behaviors such as reversible self-association (RSA), high viscosity, and liquid-liquid phase separation can introduce substantial challenges during development of monoclonal antibody formulations. Although a global mechanistic understanding of RSA (i.e., native and reversible protein-protein interactions) is sufficient to develop robust formulation controls, its mitigation via protein engineering requires knowledge of the sites of protein-protein interactions. In the study reported here, we coupled our previous hydrogen-deuterium exchange mass spectrometry findings with structural modeling and in vitro screening to identify the residues responsible for RSA of a model IgG1 monoclonal antibody (mAb-C), and rationally engineered variants with improved solution properties (i.e., reduced RSA and viscosity). Our data show that mutation of either solvent-exposed aromatic residues within the heavy and light chain variable regions or buried residues within the heavy chain/light chain interface can significantly mitigate RSA and viscosity by reducing the IgG's surface hydrophobicity. The engineering strategy described here highlights the utility of integrating complementary experimental and in silico methods to identify mutations that can improve developability, in particular, high concentration solution properties, of candidate therapeutic antibodies.