Concentration dependent viscosity of monoclonal antibody solutions: explaining experimental behavior in terms of molecular properties

A framework for the biophysical screening of antibody mutations targeting solvent-accessible hydrophobic and electrostatic patches for enhanced viscosity profiles

The formulation of high-concentration monoclonal antibody (mAb) solutions in low dose volumes for autoinjector devices poses challenges in manufacturability and patient administration due to elevated solution viscosity. Often many therapeutically potent mAbs are discovered, but their commercial development is stalled by unfavourable developability challenges. In this work, we present a systematic experimental framework for the computational screening of molecular descriptors to guide the design of 24 mutants with modified viscosity profiles accompanied by experimental evaluation. Our experimental observations using a model anti-IL8 mAb and eight engineered mutant variants reveal that viscosity reduction is influenced by the location of hydrophobic interactions, while targeting positively charged patches significantly increases viscosity in comparison to wild-type anti-IL-8 mAb. We conclude that most predicted in silico physicochemical properties exhibit poor correlation with measured experimental parameters for antibodies with suboptimal developability characteristics, emphasizing the need for comprehensive case-by-case evaluation of mAbs. This framework combining molecular design and triage via computational predictions with experimental evaluation aids the agile and rational design of mAbs with tailored solution viscosities, ensuring improved manufacturability and patient convenience in self-administration scenarios.

Current state of implementation of in silico tools in the biopharmaceutical industry-Proceedings of the 5th modeling workshop

The fifth modeling workshop (5MW) was held in June 2023 at Favrholm, Denmark and sponsored by Recovery of Biological Products Conference Series. The goal of the workshop was to assemble modeling practitioners to review and discuss the current state, progress since the last fourth mini modeling workshop (4MMW), gaps and opportunities for development, deployment and maintenance of models in bioprocess applications. Areas of focus were four categories: biophysics and molecular modeling, mechanistic modeling, computational fluid dynamics (CFD) and plant modeling. Highlights of the workshop included significant advancements in biophysical/molecular modeling to novel protein constructs, mechanistic models for filtration and initial forays into modeling of multiphase systems using CFD for a bioreactor and mapped strategically to cell line selection/facility fit. A significant impediment to more fully quantitative and calibrated models for biophysics is the lack of large, anonymized datasets. A potential solution would be the use of specific descriptors in a database that would allow for detailed analyzes without sharing proprietary information. Another gap identified was the lack of a consistent framework for use of models that are included or support a regulatory filing beyond the high-level guidance in ICH Q8-Q11. One perspective is that modeling can be viewed as a component or precursor of machine learning (ML) and artificial intelligence (AI). Another outcome was alignment on a key definition for "mechanistic modeling." Feedback from participants was that there was progression in all of the fields of modeling within scope of the conference. Some areas (e.g., biophysics and molecular modeling) have opportunities for significant research investment to realize full impact. However, the need for ongoing research and development for all model types does not preclude the application to support process development, manufacturing and use in regulatory filings. Analogous to ML and AI, given the current state of the four modeling types, a prospective investment in educating inter-disciplinary subject matter experts (e.g., data science, chromatography) is essential to advancing the modeling community.

Reduction of monoclonal antibody viscosity using interpretable machine learning

Early identification of antibody candidates with drug-like properties is essential for simplifying the development of safe and effective antibody therapeutics. For subcutaneous administration, it is important to identify candidates with low self-association to enable their formulation at high concentration while maintaining low viscosity, opalescence, and aggregation. Here, we report an interpretable machine learning model for predicting antibody (IgG1) variants with low viscosity using only the sequences of their variable (Fv) regions. Our model was trained on antibody viscosity data (>100 mg/mL mAb concentration) obtained at a common formulation pH (pH 5.2), and it identifies three key Fv features of antibodies linked to viscosity, namely their isoelectric points, hydrophobic patch sizes, and numbers of negatively charged patches. Of the three features, most predicted antibodies at risk for high viscosity, including antibodies with diverse antibody germlines in our study (79 mAbs) as well as clinical-stage IgG1s (94 mAbs), are those with low Fv isoelectric points (Fv pIs < 6.3). Our model identifies viscous antibodies with relatively high accuracy not only in our training and test sets, but also for previously reported data. Importantly, we show that the interpretable nature of the model enables the design of mutations that significantly reduce antibody viscosity, which we confirmed experimentally. We expect that this approach can be readily integrated into the drug development process to reduce the need for experimental viscosity screening and improve the identification of antibody candidates with drug-like properties.

Prediction of Antibody Viscosity from Dilute Solution Measurements

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

Application of Formulation Principles to Stability Issues Encountered During Processing, Manufacturing, and Storage of Drug Substance and Drug Product Protein Therapeutics

The field of formulation and stabilization of protein therapeutics has become rather extensive. However, most of the focus has been on stabilization of the final drug product. Yet, proteins experience stress and degradation through the manufacturing process, starting with fermentaition. This review describes how formulation principles can be applied to stabilize biopharmaceutical proteins during bioprocessing and manufacturing, considering each unit operation involved in prepration of the drug substance. In addition, the impact of the container on stabilty is discussed as well.

Biophysical Determinants for the Viscosity of Concentrated Monoclonal Antibody Solutions

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

Low-data interpretable deep learning prediction of antibody viscosity using a biophysically meaningful representation

Deep learning, aided by the availability of big data sets, has led to substantial advances across many disciplines. However, many scientific problems of practical interest lack sufficiently large datasets amenable to deep learning. Prediction of antibody viscosity is one such problem where deep learning methods have not yet been explored due to the relative scarcity of relevant training data. In this work, we overcome this limitation using a biophysically meaningful representation that enables us to develop generalizable models even under limited training data. We present, PfAbNet-viscosity, a 3D convolutional neural network architecture, to predict high-concentration viscosity of therapeutic antibodies. We show that with the electrostatic potential surface of the antibody variable region as the only input to the network, the models trained on as few as couple dozen datapoints can generalize with high accuracy. Our feature attribution analysis shows that PfAbNet-viscosity has learned key biophysical drivers of viscosity. The applicability of our approach to other biological systems is discussed.

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

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

Relationship between aggregation of therapeutic proteins and agitation parameters: Acceleration and frequency

An increase in protein aggregates during transportation should be suppressed in therapeutic protein products because the aggregates have a potential risk of immunogenicity. In this study, three protein solutions in vials were exposed to tri-axial vibration with various combinations of frequency and acceleration using a transportation test system to investigate the relationship between low g-force stresses and protein aggregate generation. The number concentration of micron aggregates detected by flow imaging analysis increased markedly when the acceleration and frequency of agitation were within a specific range, in other words, above a threshold. This threshold was common among the three protein solutions. The suppression of micron aggregate formation by adding a surfactant suggested that agitation above the threshold increased micron aggregates mainly via interface-mediated routes. Notably, agitation, including agitation below the threshold, accelerated spontaneous oligomerization (nanometer aggregate generation) of proteins in bulk solution even in the presence of the surfactant. Studies of stability against mechanical stresses (e.g., a random vibration test to simulate actual shipment, with a time-compressed setting by increasing acceleration) need to be performed and discussed with careful consideration of the threshold for generating micron aggregates.

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

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

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

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

High concentration formulation developability approaches and considerations

The growing need for biologics to be administered subcutaneously and ocularly, coupled with certain indications requiring high doses, has resulted in an increase in drug substance (DS) and drug product (DP) protein concentrations. With this increase, more emphasis must be placed on identifying critical physico-chemical liabilities during drug development, including protein aggregation, precipitation, opalescence, particle formation, and high viscosity. Depending on the molecule, liabilities, and administration route, different formulation strategies can be used to overcome these challenges. However, due to the high material requirements, identifying optimal conditions can be slow, costly, and often prevent therapeutics from moving rapidly into the clinic/market. In order to accelerate and derisk development, new experimental and in-silico methods have emerged that can predict high concentration liabilities. Here, we review the challenges in developing high concentration formulations, the advances that have been made in establishing low mass and high-throughput predictive analytics, and advances in in-silico tools and algorithms aimed at identifying risks and understanding high concentration protein behavior.

Embedding Dynamics in Intrinsic Physicochemical Profiles of Market-Stage Antibody-Based Biotherapeutics

Adequate stability, manufacturability, and safety are crucial to bringing an antibody-based biotherapeutic to the market. Following the concept of holistic in silico developability, we introduce a physicochemical description of 91 market-stage antibody-based biotherapeutics based on orthogonal molecular properties of variable regions (Fvs) embedded in different simulation environments, mimicking conditions experienced by antibodies during manufacturing, formulation, and in vivo. In this work, the evaluation of molecular properties includes conformational flexibility of the Fvs using molecular dynamics (MD) simulations. The comparison between static homology models and simulations shows that MD significantly affects certain molecular descriptors like surface molecular patches. Moreover, the structural stability of a subset of Fv regions is linked to changes in their specific molecular interactions with ions in different experimental conditions. This is supported by the observation of differences in protein melting temperatures upon addition of NaCl. A DEvelopability Navigator In Silico (DENIS) is proposed to compare mAb candidates for their similarity with market-stage biotherapeutics in terms of physicochemical properties and conformational stability. Expanding on our previous developability guidelines (Ahmed et al. Proc. Natl. Acad. Sci. 2021, 118 (37), e2020577118), the hydrodynamic radius and the protein strand ratio are introduced as two additional descriptors that enable a more comprehensive in silico characterization of biotherapeutic drug candidates. Test cases show how this approach can facilitate identification and optimization of intrinsically developable lead candidates. DENIS represents an advanced computational tool to progress biotherapeutic drug candidates from discovery into early development by predicting drug properties in different aqueous environments.

Ionic Liquid-Based Strategy for Predicting Protein Aggregation Propensity and Thermodynamic Stability

Novel drug candidates are continuously being developed to combat the most life-threatening diseases; however, many promising protein therapeutics are dropped from the pipeline. During biological and industrial processes, protein therapeutics are exposed to various stresses such as fluctuations in temperature, solvent pH, and ionic strength. These can lead to enhanced protein aggregation propensity, one of the greatest challenges in drug development. Recently, ionic liquids (ILs), in particular, biocompatible choline chloride ([Cho]Cl)-based ILs, have been used to hinder stress-induced protein conformational changes. Herein, we develop an IL-based strategy to predict protein aggregation propensity and thermodynamic stability. We examine three key variables influencing protein misfolding: pH, ionic strength, and temperature. Using dynamic light scattering, zeta potential, and variable temperature circular dichroism measurements, we systematically evaluate the structural, thermal, and thermodynamic stability of fresh immunoglobin G4 (IgG4) antibody in water and 10, 30, and 50 wt % [Cho]Cl. Additionally, we conduct molecular dynamics simulations to examine IgG4 aggregation propensity in each system and the relative favorability of different [Cho]Cl-IgG4 packing interactions. We re-evaluate each system following 365 days of storage at 4 °C and demonstrate how to predict the thermodynamic properties and protein aggregation propensity over extended storage, even under stress conditions. We find that increasing [Cho]Cl concentration reduced IgG4 aggregation propensity both fresh and following 365 days of storage and demonstrate the potential of using our predictive IL-based strategy and formulations to radically increase protein stability and storage.

Utilizing cross-product prior knowledge to rapidly de-risk chemical liabilities in therapeutic antibody candidates

There is considerable pressure in the pharmaceutical industry to advance better molecules faster. One pervasive concern for protein-based therapeutics is the presence of potential chemical liabilities. We have developed a simple methodology for rapidly de-risking specific chemical concerns in antibody-based molecules using prior knowledge of each individual liability at a specific position in the molecule’s sequence. Our methodology hinges on the development of sequence-aligned chemical liability databases of molecules from different stages of commercialization and on sequence-aligned experimental data from prior molecules that have been developed at Amgen. This approach goes beyond the standard practice of simply flagging all instances of each motif that fall in a CDR. Instead, we de-risk motifs that are common at a specific site in commercial mAb-based molecules (and therefore did not previously pose an insurmountable barrier to commercialization) and motifs at specific sites for which we have prior experimental data indicating acceptably low levels of modification. We have used this approach successfully to identify candidates in a discovery phase program with exclusively very low risk potential chemical liabilities. Identifying these candidates in the discovery phase allowed us to bypass protein engineering and accelerate the program’s timeline by 6 months.

Computational models for studying physical instabilities in high concentration biotherapeutic formulations

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

Antibodies with Weakly Basic Isoelectric Points Minimize Trade-offs between Formulation and Physiological Colloidal Properties

The widespread interest in antibody therapeutics has led to much focus on identifying antibody candidates with favorable developability properties. In particular, there is broad interest in identifying antibody candidates with highly repulsive self-interactions in standard formulations (e.g., low ionic strength buffers at pH 5-6) for high solubility and low viscosity. Likewise, there is also broad interest in identifying antibody candidates with low levels of non-specific interactions in physiological solution conditions (PBS, pH 7.4) to promote favorable pharmacokinetic properties. To what extent antibodies that possess both highly repulsive self-interactions in standard formulations and weak non-specific interactions in physiological solution conditions can be systematically identified remains unclear and is a potential impediment to successful therapeutic drug development. Here, we evaluate these two properties for 42 IgG1 variants based on the variable fragments (Fvs) from four clinical-stage antibodies and complementarity-determining regions from 10 clinical-stage antibodies. Interestingly, we find that antibodies with the strongest repulsive self-interactions in a standard formulation (pH 6 and 10 mM histidine) display the strongest non-specific interactions in physiological solution conditions. Conversely, antibodies with the weakest non-specific interactions under physiological conditions display the least repulsive self-interactions in standard formulations. This behavior can be largely explained by the antibody isoelectric point, as highly basic antibodies that are highly positively charged under standard formulation conditions (pH 5-6) promote repulsive self-interactions that mediate high colloidal stability but also mediate strong non-specific interactions with negatively charged biomolecules at physiological pH and vice versa for antibodies with negatively charged Fv regions. Therefore, IgG1s with weakly basic isoelectric points between 8 and 8.5 and Fv isoelectric points between 7.5 and 9 typically display the best combinations of strong repulsive self-interactions and weak non-specific interactions. We expect that these findings will improve the identification and engineering of antibody candidates with drug-like biophysical properties.

Utility of High Resolution 2D NMR Fingerprinting in Assessing Viscosity of Therapeutic Monoclonal Antibodies

Purpose

The viscosity of highly concentrated therapeutic monoclonal antibody (mAb) formulations at concentrations ≥ 100 mg/mL can significantly affect the stability, processing, and drug product development for subcutaneous delivery. An early identification of a viscosity prone mAb during candidate selection stages are often beneficial for downstream processes. Higher order structure of mAbs may often dictate their viscosity behavior at high concentration. Thus it is beneficial to gauge or rank-order their viscosity behavior using noninvasive structural fingerprinting methods and to potentially screen for suitable viscosity lowering excipients.

Methods

In this study, Dynamic Light Scattering (DLS) and 2D NMR based methyl fingerprinting were used to correlate viscosity behavior of a set of Pfizer mAbs. The viscosities of mAbs were determined. Respective Fab and Fc domains were generated for studies.

Result

Methyl fingerprinting of intact mAbs allows for differentiation of viscosity prone mAbs from well behaved ones even at 30–40 mg/ml, where bulk viscosity of the solutions are near identical. For viscosity prone mAbs, peak broadening and or distinct chemical shift changes were noted in intact and fragment fingerprints, unlike the well-behaved mAbs, indicative of protein protein interactions (PPI).

Conclusion

Fab-Fab or Fab-Fc interactions may lead to formation of protein networks at high concentration. The early transients to these network formation may be manifested through peak broadening or peak shift in the 2D NMR spectrum of mAb/mAb fragments. Such insights go beyond rank ordering mAbs based on viscosity behavior, which can be obtained by other methods as well..

Supplementary Information

The online version contains supplementary material available at 10.1007/s11095-022-03200-6.

Current advances in biopharmaceutical informatics: guidelines, impact and challenges in the computational developability assessment of antibody therapeutics

Therapeutic monoclonal antibodies and their derivatives are key components of clinical pipelines in the global biopharmaceutical industry. The availability of large datasets of antibody sequences, structures, and biophysical properties is increasingly enabling the development of predictive models and computational tools for the "developability assessment" of antibody drug candidates. Here, we provide an overview of the antibody informatics tools applicable to the prediction of developability issues such as stability, aggregation, immunogenicity, and chemical degradation. We further evaluate the opportunities and challenges of using biopharmaceutical informatics for drug discovery and optimization. Finally, we discuss the potential of developability guidelines based on in silico metrics that can be used for the assessment of antibody stability and manufacturability.

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

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

Differences in human IgG1 and IgG4 S228P monoclonal antibodies viscosity and self-interactions: Experimental assessment and computational predictions of domain interactions

Human/humanized IgG4 antibodies have reduced effector function relative to IgG1 antibodies, which is desirable for certain therapeutic purposes. However, the developability and biophysical properties for IgG4 antibodies are not well understood. This work focuses on the head-to-head comparison of key biophysical properties, such as self-interaction and viscosity, for 14 human/humanized, and chimeric IgG1 and IgG4 S228P monoclonal antibody pairs that contain the identical variable regions. Experimental measurements showed that the IgG4 S228P antibodies have similar or higher self-interaction and viscosity than that of IgG1 antibodies in 20 mM sodium acetate, pH 5.5. We report sequence and structural drivers for the increased viscosity and self-interaction detected in IgG4 S228P antibodies through a combination of experimental data and computational models. Further, we applied and extended a previously established computational model for IgG1 antibodies to predict the self-interaction and viscosity behavior for each antibody pair, providing insight into the structural characteristics and differences of these two isotypes. Interestingly, we observed that the IgG4 S228P swapped variants, where the CH3 domain was swapped for that of an IgG1, showed reduced self-interaction behavior. These domain swapped IgG4 S228P molecules also showed reduced viscosity from experiment and coarse-grained simulations. We also observed that experimental diffusion interaction parameter (kD) values have a high correlation with computational diffusivity prediction for both IgG1 and IgG4 S228P isotypes.

Abbreviations: AHc, constant region Hamaker constant; AHv, variable region Hamaker constant; CDRs, Complementarity-determining regions; CG, Coarse-grained model; CH1, Constant heavy chain 1; CH2 Constant heavy chain 2; CH3 Constant heavy chain 3; chgCH3 Effective charge on the CH3 region; CL Constant light chain; cP, Centipoise; DLS, Dynamic light scattering; Fab, Fragment antigen-binding; Fc, Fragment crystallizable; Fv, Variable domaing; (r) Radial distribution function; H1 CDR1 of Heavy Chain; H2 CDR2 of Heavy Chain; H3 CDR3 of Heavy Chain; HVI, High viscosity index; IgG1 human immunoglobulin of IgG1 subclass; IgG4 human immunoglobulin of IgG4 subclass; kD, Diffusion interaction parameter; L1 CDR1 of Light Chain; L2 CDR2 of Light Chain; L3 CDR3 of Light Chain; mAb, Monoclonal antibody; MD, Molecular dynamics; PPI Protein–protein interactions; SCM, Spatial charge map; UP-SEC, Ultra-high-performance size-exclusion chromatography; VH, Variable domain of Heavy Chain; VL, Variable domain of Light Chain

Stress-dependent flexibility of a full-length human monoclonal antibody: Insights from molecular dynamics to support biopharmaceutical development

After several decades of advancements in drug discovery, product development of biopharmaceuticals remains a time- and resource-consuming endeavor. One of the main reasons is associated to the lack of fundamental understanding of conformational dynamics of such biologic entities, and how they respond to various stresses encountered during manufacturing. In this work, we have studied the conformational dynamics of human IgG₁κ b12 monoclonal antibody (mAb) using molecular dynamics simulations. The hundreds of nanoseconds long trajectories reveal that b12 mAb is highly flexible. Its variable domains show greater conformational fluctuations than the constant domains. Additionally, it collapses towards a more globular shape in response to thermal stress, leading to decrease in the total solvent exposed surface area and radius of gyration. This behavior is more pronounced for the deglycosylated b12 mAb, and it appears to correlate with increase in inter-domain contacts between specific regions of the antibody. Conformational fluctuations also cause temporary formation and disruption of hydrophobic and charged patches on the antibody surface, which is particularly important for the prediction of CMC properties during development phases of antibody-based biotherapeutics. The insights gained through these simulations may help the development of biologic drugs, especially with regards to manufacturing processes where antibodies may undergo significant thermal stress.

Nanoencapsulation as a General Solution for Lyophilization of Labile Substrates

Protein macromolecules occur naturally at the nanoscale. The use of a dedicated nanoparticle as a lyophilization excipient, however, has not been reported. Because biopolymeric and lipid nanoparticles often denature protein macromolecules and commonly lack the structural rigidity to survive the freeze-drying process, we hypothesized that surrounding an individual protein substrate with a nanoscale, thermostable exoshell (tES) would prevent aggregation and protect the substrate from denaturation during freezing, sublimation, and storage. We systematically investigated the properties of tES, including secondary structure and its homogeneity, throughout the process of lyophilization and found that tES have a near 100% recovery following aqueous reconstitution. We then tested the hypothesis that tES could encapsulate a model substrate, horseradish peroxidase (HRP), using charge complementation and pH-mediated controlled assembly. HRP were encapsulated within the 8 nm internal tES aqueous cavity using a simplified loading procedure. Time-course experiments demonstrated that unprotected HRP loses 95% of activity after 1 month of lyophilized storage. After encapsulation within tES nanoparticles, 70% of HRP activity was recovered, representing a 14-fold improvement and this effect was reproducible across a range of storage temperatures. To our knowledge, these results represent the first reported use of nanoparticle encapsulation to stabilize a functional macromolecule during lyophilization. Thermostable nanoencapsulation may be a useful method for the long-term storage of labile proteins.

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.

Intrinsic physicochemical profile of marketed antibody-based biotherapeutics

Successful biologic drug discovery and development involves finding functional as well as developable candidates. Once a candidate has been demonstrated to be functional, the next step is to determine whether it can be translated into a drug product. This requires that the candidate can withstand stresses encountered during manufacturing, shipping, and storage. Additionally, it must be safe, efficacious, and possess good pharmacology. In silico analyses of the variable regions of 77 marketed antibody-based biotherapeutics have revealed five nonredundant physicochemical descriptors. Distributions of these descriptors, observed for marketed biotherapeutics, can help prioritize a drug candidate for experimental testing at early discovery stages, guide engineering efforts to further optimize it, and help increase the productivity of biologic drug discovery and development.

Feeding biopharma pipelines with biotherapeutic candidates that possess desirable developability profiles can help improve the productivity of biologic drug discovery and development. Here, we have derived an in silico profile by analyzing computed physicochemical descriptors for the variable regions (Fv) found in 77 marketed antibody-based biotherapeutics. Fv regions of these biotherapeutics demonstrate significant diversities in their germlines, complementarity determining region loop lengths, hydrophobicity, and charge distributions. Furthermore, an analysis of 24 physicochemical descriptors, calculated using homology-based molecular models, has yielded five nonredundant descriptors whose distributions represent stability, isoelectric point, and molecular surface characteristics of their Fv regions. Fv regions of candidates from our internal discovery campaigns, human next-generation sequencing repertoires, and those in clinical-stages (CST) were assessed for similarity with the physicochemical profile derived here. The Fv regions in 33% of CST antibodies show physicochemical properties that are dissimilar to currently marketed biotherapeutics. In comparison, physicochemical characteristics of ∼29% of the Fv regions in human antibodies and ∼27% of our internal hits deviated significantly from those of marketed biotherapeutics. The early availability of this information can help guide hit selection, lead identification, and optimization of biotherapeutic candidates. Insights from this work can also help support portfolio risk assessment, in-licensing, and biopharma collaborations.

Molecular computations of preferential interactions of proline, arginine.HCl, and NaCl with IgG1 antibodies and their impact on aggregation and viscosity

Preferential interactions of excipients with the antibody surface govern their effect on the stability of antibodies in solution. We probed the preferential interactions of proline, arginine.HCl (Arg.HCl), and NaCl with three therapeutically relevant IgG1 antibodies via experiment and simulation. With simulations, we examined how excipients interacted with different types of surface patches in the variable region (Fv). For example, proline interacted most strongly with aromatic surfaces, Arg.HCl was included near negative residues, and NaCl was excluded from negative residues and certain hydrophobic regions. The differences in interaction of different excipients with the same surface patch on an antibody may be responsible for variations in the antibody's aggregation, viscosity, and self-association behaviors in each excipient. Proline reduced self-association for all three antibodies and reduced aggregation for the antibody with an association-limited aggregation mechanism. The effects of Arg.HCl and NaCl on aggregation and viscosity were highly dependent on the surface charge distribution and the extent of exclusion from highly hydrophobic patches. At pH 5.5, both tended to increase the aggregation of an antibody with a strongly positive charge on the Fv, while only NaCl reduced the aggregation of the antibody with a large negative charge patch on the Fv. Arg.HCl reduced the viscosities of antibodies with either a hydrophobicity-driven mechanism or a charge-driven mechanism. Analysis of this data presents a framework for understanding how amino acid and ionic excipients interact with different protein surfaces, and how these interactions translate to the observed stability behavior.

Protein-protein interactions at high concentrations: Effects of ArgHCl and NaCl on the stability, viscosity and aggregation mechanisms of protein solution

The aim of this work was to use the diffusion coefficient ration (D_m/D_line) as a parameter to characterize the stability of protein at high concentration, to compare the effects of ArgHCl and NaCl on the interaction of highly concentrated proteins under different pH conditions, and to explore the correlation with protein stability. For this purpose, a high-concentration bovine serum albumin solution (BSA) was selected as the model system, and the diffusion coefficient, aggregation degree, conformational stability, and solution viscosity of the protein were studied by dynamic light scattering (DLS) and spectral detection techniques. The result showed that there was a significant correlation between the D_m/D_line and the protein aggregation. The D_m/D_line of the protein was minimum at pH 7.4, which corresponded to the maximum degree of aggregation and the highest solution viscosity. At pH 7.4, the hydrophobic interactions and the increased conformational stability of ArgHCl maximized the stability of the protein and reduced the viscosity of the solution by 69.3%. At pH 3.0, the strong charge shielding effect of ArgHCl and NaCl and the decreased conformational stability induced protein aggregation and the gel formation. These findings provided valuable insights into the mechanism of protein aggregation and the diffusion coefficient ration (D_m/D_line) could be a potential tool for the pre-formulation studies.

Characterization and Modeling of Reversible Antibody Self-Association Provide Insights into Behavior, Prediction, and Correction

Reversible antibody self-association, while having major developability and therapeutic implications, is not fully understood or readily predictable and correctable. For a strongly self-associating humanized mAb variant, resulting in unacceptable viscosity, the monovalent affinity of self-interaction was measured in the low μM range, typical of many specific and biologically relevant protein-protein interactions. A face-to-face interaction model extending across both the heavy-chain (HC) and light-chain (LC) Complementary Determining Regions (CDRs) was apparent from biochemical and mutagenesis approaches as well as computational modeling. Light scattering experiments involving individual mAb, Fc, Fab, and Fab'2 domains revealed that Fabs self-interact to form dimers, while bivalent mAb/Fab'2 forms lead to significant oligomerization. Site-directed mutagenesis of aromatic residues identified by homology model patch analysis and self-docking dramatically affected self-association, demonstrating the utility of these predictive approaches, while revealing a highly specific and tunable nature of self-binding modulated by single point mutations. Mutagenesis at these same key HC/LC CDR positions that affect self-interaction also typically abolished target binding with notable exceptions, clearly demonstrating the difficulties yet possibility of correcting self-association through engineering. Clear correlations were also observed between different methods used to assess self-interaction, such as Dynamic Light Scattering (DLS) and Affinity-Capture Self-Interaction Nanoparticle Spectroscopy (AC-SINS). Our findings advance our understanding of therapeutic protein and antibody self-association and offer insights into its prediction, evaluation and corrective mitigation to aid therapeutic development.

Machine Learning Applied to Determine the Molecular Descriptors Responsible for the Viscosity Behavior of Concentrated Therapeutic Antibodies

Predicting the solution viscosity of monoclonal antibody (mAb) drug products remains as one of the main challenges in antibody drug design, manufacturing, and delivery. In this work, the concentration-dependent solution viscosity of 27 FDA-approved mAbs was measured at pH 6.0 in 10 mM histidine-HCl. Six mAbs exhibited high viscosity (>30 cP) in solutions at 150 mg/mL mAb concentration. Combining molecular modeling and machine learning feature selection, we found that the net charge in the mAbs and the amino acid composition in the Fv region are key features which govern the viscosity behavior. For mAbs whose behavior was not dominated by charge effects, we observed that high viscosity is correlated with more hydrophilic and fewer hydrophobic residues in the Fv region. A predictive model based on the net charges of mAbs and a high viscosity index is presented as a fast screening tool for classifying low- and high-viscosity mAbs.

Generation of mAb Variants with Less Attractive Self-Interaction but Preserved Target Binding by Well-Directed Mutation

Strongly attractive self-interaction of therapeutic protein candidates can impose challenges for manufacturing, filling, stability, and administration due to elevated viscosity or aggregation propensity. Suitable formulations can mitigate these issues to a certain extent. Understanding the self-interaction mechanism on a molecular basis and rational protein engineering provides a more fundamental approach, and it can save costs and efforts as well as alleviate risks at later stages of development. In this study, we used computational methods for the identification of aggregation-prone regions in a mAb and generated mutants based on these findings. We applied hydrogen-deuterium exchange mass spectrometry to identify distinct self-interaction hot spots. Ultimately, we generated mAb variants based on a combination of both approaches and identified mutants with low attractive self-interaction propensity, minimal off-target binding, and even improved target binding. Our data show that the introduction of arginine in spatial proximity to hydrophobic patches is highly beneficial on all these levels. For our mAb, variants that contain more than one aspartate residue flanking to the hydrophobic HCDR3 show decreased attractive self-interaction at unaffected off-target and target binding. The combined engineering strategy described here underlines the high potential of understanding self-interaction in the early stages of development to predict and reduce the risk of failure in subsequent development.

Prediction Machines: Applied Machine Learning for Therapeutic Protein Design and Development

The rapid growth in technological advances and quantity of scientific data over the past decade has led to several challenges including data storage and analysis. Accurate models of complex datasets were previously difficult to develop and interpret. However, improvements in machine learning algorithms have since enabled unparalleled classification and prediction capabilities. The application of machine learning can be seen throughout diverse industries due to their ease of use and interpretability. In this review, we describe popular machine learning algorithms and highlight their application in pharmaceutical protein development. Machine learning models have now been applied to better understand the nonlinear concentration dependent viscosity of protein solutions, predict protein oxidation and deamidation rates, classify sub-visible particles and compare the physical stability of proteins. We also applied several machine learning algorithms using previously published data and describe models with improved predictions and classification. The authors hope that this review can be used as a resource to others and encourage continued application of machine learning algorithms to problems in pharmaceutical protein development.

Analysis of antibody self-interaction by bio-layer interferometry as tool to support lead candidate selection during preformulation and developability assessments

Developability assessment of therapeutic mAb candidates before entering CMC development mitigates the risk of later failure because of manufacturing and stability issues. For mAbs derived from library based screenings, such evaluation starts with the first panning and ends with the selection of a lead candidate. This candidate should show, amongst others, high affine target binding and beneficial conformational as well as chemical stability. In addition, colloidal stability, reflected by the self-interaction propensity, should be superior in order to reduce aggregate formation and unacceptably high viscosity at elevated protein concentrations. Here, we present a study demonstrating the application of self-interaction bio-layer interferometry (SI-BLI) in a developability assessment, including the evaluation of preformulations. We reveal that the formulation rankings based on SI-BLI, DLS and viscosity measurements correlate. SI-BLI provides a deeper understanding of influencing factors on mAb self-interaction such as ionic strength or cation species. The attractive mAb self-interaction propensity was significantly more suppressed by Mg²⁺ compared to Na⁺. SI-BLI can be performed in high throughput with minimal material and sample preparation needs. Therefore, it can be applied in early stages of developability assessment going beyond the use of a platform formulation and a small number of analysis, to screen more parameters before proceeding with candidate selection and further extensive development.

Toward Drug-Like Multispecific Antibodies by Design

The success of antibody therapeutics is strongly influenced by their multifunctional nature that couples antigen recognition mediated by their variable regions with effector functions and half-life extension mediated by a subset of their constant regions. Nevertheless, the monospecific IgG format is not optimal for many therapeutic applications, and this has led to the design of a vast number of unique multispecific antibody formats that enable targeting of multiple antigens or multiple epitopes on the same antigen. Despite the diversity of these formats, a common challenge in generating multispecific antibodies is that they display suboptimal physical and chemical properties relative to conventional IgGs and are more difficult to develop into therapeutics. Here we review advances in the design and engineering of multispecific antibodies with drug-like properties, including favorable stability, solubility, viscosity, specificity and pharmacokinetic properties. We also highlight emerging experimental and computational methods for improving the next generation of multispecific antibodies, as well as their constituent antibody fragments, with natural IgG-like properties. Finally, we identify several outstanding challenges that need to be addressed to increase the success of multispecific antibodies in the clinic.

Physicochemical Rules for Identifying Monoclonal Antibodies with Drug-like Specificity

The ability of antibodies to recognize their target antigens with high specificity is fundamental to their natural function. Nevertheless, therapeutic antibodies display variable and difficult-to-predict levels of nonspecific and self-interactions that can lead to various drug development challenges, including antibody aggregation, abnormally high viscosity, and rapid antibody clearance. Here we report a method for predicting the overall specificity of antibodies in terms of their relative risk for displaying high levels of nonspecific or self-interactions at physiological conditions. We find that individual and combined sets of chemical rules that limit the maximum and minimum numbers of certain solvent-exposed amino acids in antibody variable regions are strong predictors of specificity for large panels of preclinical and clinical-stage antibodies. We also demonstrate how the chemical rules can be used to identify sites that mediate nonspecific interactions in suboptimal antibodies and guide the design of targeted sublibraries that yield variants with high antibody specificity. These findings can be readily used to improve the selection and engineering of antibodies with drug-like specificity.

Highland games: A benchmarking exercise in predicting biophysical and drug properties of monoclonal antibodies from amino acid sequences

Biopharmaceutical product and process development do not yet take advantage of predictive computational modeling to nearly the degree seen in industries based on smaller molecules. To assess and advance progress in this area, spirited coopetition (mutually beneficial collaboration between competitors) was successfully used to motivate industrial scientists to develop, share, and compare data and methods which would normally have remained confidential. The first "Highland Games" competition was held in conjunction with the October 2018 Recovery of Biological Products Conference in Ashville, NC, with the goal of benchmarking and assessment of the ability to predict development-related properties of six antibodies from their amino acid sequences alone. Predictions included purification-influencing properties such as isoelectric point and protein A elution pH, and biophysical properties such as stability and viscosity at very high concentrations. Essential contributions were made by a large variety of individuals, including companies which consented to provide antibody amino acid sequences and test materials, volunteers who undertook the preparation and experimental characterization of these materials, and prediction teams who attempted to predict antibody properties from sequence alone. Best practices were identified and shared, and areas in which the community excels at making predictions were identified, as well as areas presenting opportunities for considerable improvement. Predictions of isoelectric point and protein A elution pH were especially good with all-prediction average errors of 0.2 and 1.6 pH unit, respectively, while predictions of some other properties were notably less good. This manuscript presents the events, methods, and results of the competition, and can serve as a tutorial and as a reference for in-house benchmarking by others. Organizations vary in their policies concerning disclosure of methods, but most managements were very cooperative with the Highland Games exercise, and considerable insight into common and best practices is available from the contributed methods. The accumulated data set will serve as a benchmarking tool for further development of in silico prediction tools.

Modeling and mitigation of high-concentration antibody viscosity through structure-based computer-aided protein design

For an antibody to be a successful therapeutic many competing factors require optimization, including binding affinity, biophysical characteristics, and immunogenicity risk. Additional constraints may arise from the need to formulate antibodies at high concentrations (>150 mg/ml) to enable subcutaneous dosing with reasonable volume (ideally <1.0 mL). Unfortunately, antibodies at high concentrations may exhibit high viscosities that place impractical constraints (such as multiple injections or large needle diameters) on delivery and impede efficient manufacturing. Here we describe the optimization of an anti-PDGF-BB antibody to reduce viscosity, enabling an increase in the formulated concentration from 80 mg/ml to greater than 160 mg/ml, while maintaining the binding affinity. We performed two rounds of structure guided rational design to optimize the surface electrostatic properties. Analysis of this set demonstrated that a net-positive charge change, and disruption of negative charge patches were associated with decreased viscosity, but the effect was greatly dependent on the local surface environment. Our work here provides a comprehensive study exploring a wide sampling of charge-changes in the Fv and CDR regions along with targeting multiple negative charge patches. In total, we generated viscosity measurements for 40 unique antibody variants with full sequence information which provides a significantly larger and more complete dataset than has previously been reported.

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.

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.

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.