In Silico Prediction of Diffusion Interaction Parameter (kD), a Key Indicator of Antibody Solution Behaviors.

PURPOSE: To develop resource-sparing in silico approaches that aim to reduce experimental effort and material required by developability assessments (DA) of monoclonal antibody (mAb) drug candidates. METHODS: A battery of standardized biophysical experiments was performed on high concentration formulations of 16 drug product development stage mAbs using a platform buffer. Full-length molecular models of these mAbs were also generated via molecular modeling. These models were used to computationally estimate molecular descriptors of these 16 mAbs. Pairwise and multi-parameter correlations among experimentally measured biophysical attributes and calculated molecular descriptors were obtained via statistical analyses. RESULTS: Diffusion interaction parameter (kD) showed statistically significant pairwise correlations (p-values <0.005) with thermal stability, viscosity, isoelectric point, and apparent solubility of the antibodies in our dataset. kD also showed statistically significant pairwise correlations (p-values <0.005) with several computationally calculated molecular descriptors (pI, net charge, charge on the Fv region, and zeta potential.) These pairwise correlations were further refined by multivariate analyses. These analyses yielded several useful equations for prediction of kD from antibody sequences, structural models, and experimentally measured biophysical attributes. CONCLUSIONS: Diffusion interaction parameter (kD) was found to be a key biophysical property for the mAbs in our dataset. It connects conformational heterogeneity of an antibody with its colloidal and rheological behaviors. The equations derived in this work shall enable rapid, resource-sparing, and cost-effective DAs of biologic drug candidates.

Solvophobic and solvophilic contributions in the water-to-aqueous guanidinium chloride transfer free energy of model peptides.

We study the solvation free energy of two different conformations (helix and extended) of two different peptides (deca-alanine and deca-glycine) in two different solvents (water and aqueous guanidinium chloride, GdmCl). The free energies are obtained using the quasichemical organization of the potential distribution theorem, an approach that naturally provides the repulsive (solvophobic or cavity) and attractive (solvophilic) contributions to solvation. The solvophilic contribution is further parsed into a chemistry contribution arising from solute interaction with the solvent in the first solvation shell and a long-range contribution arising from non-specific interactions between the solute and the solvent beyond the first solvation shell. The cavity contribution is obtained for two different envelopes, ΣSE, which theory helps identify as the solvent excluded volume, and ΣG, a larger envelope beyond which solute-solvent interactions are Gaussian. The ΣSE envelope is independent of the solvent, as expected on the basis of the insensitivity to the solvent type of the distance of closest approach between protein heavy atoms and solvent heavy atoms, but contrary to the intuition based on treating solvent constituents as spheres of some effective radii. For both envelopes, the cavity contribution in water is proportional to the surface area of the envelope. The same does not hold for GdmCl(aq), revealing the limitation of using molecular area to assess solvation energetics. The ΣG-cavity contribution predicts that GdmCl(aq) should favor the more compact state, contrary to the role of GdmCl in unfolding proteins. The chemistry contribution attenuates this effect, but still the net local (chemistry plus ΣG-packing) contribution is inadequate in capturing the role of GdmCl. With the inclusion of the long-range contribution, which is dominated by van der Waals interaction, aqueous GdmCl favors the extended conformation over the compact conformation. Our finding emphasizes the importance of weak, but attractive, long-range dispersion interactions in protein solution thermodynamics.

Characterization and Higher-Order Structure Assessment of an Interchain Cysteine-Based ADC: Impact of Drug Loading and Distribution on the Mechanism of Aggregation.

The impact of drug loading and distribution on higher order structure and physical stability of an interchain cysteine-based antibody drug conjugate (ADC) has been studied. An IgG1 mAb was conjugated with a cytotoxic auristatin payload following the reduction of interchain disulfides. The 2-D LC-MS analysis shows that there is a preference for certain isomers within the various drug to antibody ratios (DARs). The physical stability of the unconjugated monoclonal antibody, the ADC, and isolated conjugated species with specific DAR, were compared using calorimetric, thermal, chemical denaturation and molecular modeling techniques, as well as techniques to assess hydrophobicity. The DAR was determined to have a significant impact on the biophysical properties and stability of the ADC. The CH2 domain was significantly perturbed in the DAR6 species, which was attributable to quaternary structural changes as assessed by molecular modeling. At accelerated storage temperatures, the DAR6 rapidly forms higher molecular mass species, whereas the DAR2 and the unconjugated mAb were largely stable. Chemical denaturation study indicates that DAR6 may form multimers while DAR2 and DAR4 primarily exist in monomeric forms in solution at ambient conditions. The physical state differences were correlated with a dramatic increase in the hydrophobicity and a reduction in the surface tension of the DAR6 compared to lower DAR species. Molecular modeling of the various DAR species and their conformers demonstrates that the auristatin-based linker payload directly contributes to the hydrophobicity of the ADC molecule. Higher order structural characterization provides insight into the impact of conjugation on the conformational and colloidal factors that determine the physical stability of cysteine-based ADCs, with implications for process and formulation development.

XTX301, a Tumor-Activated Interleukin-12 Has the Potential to Widen the Therapeutic Index of IL12 Treatment for Solid Tumors as Evidenced by Preclinical Studies.

IL12 is a proinflammatory cytokine, that has shown promising antitumor activity in humans by promoting the recruitment and activation of immune cells in tumors. However, the systemic administration of IL12 has been accompanied by considerable toxicity, prompting interest in researching alternatives to drive preferential IL12 bioactivity in the tumor. Here, we have generated XTX301, a tumor-activated IL12 linked to the human Fc protein via a protease cleavable linker that is pharmacologically inactivated by an IL12 receptor subunit beta 2 masking domain. In vitro characterization demonstrates multiple matrix metalloproteases, as well as human primary tumors cultured as cell suspensions, can effectively activate XTX301. Intravenous administration of a mouse surrogate mXTX301 demonstrated significant tumor growth inhibition (TGI) in inflamed and non-inflamed mouse models without causing systemic toxicities. The superiority of mXTX301 in mediating TGI compared with non-activatable control molecules and the greater percentage of active mXTX301 in tumors versus other organs further confirms activation by the tumor microenvironment-associated proteases in vivo. Pharmacodynamic characterization shows tumor selective increases in inflammation and upregulation of immune-related genes involved in IFNγ cell signaling, antigen processing, presentation, and adaptive immune response. XTX301 was tolerated following four repeat doses up to 2.0 mg/kg in a nonhuman primate study; XTX301 exposures were substantially higher than those at the minimally efficacious dose in mice. Thus, XTX301 has the potential to achieve potent antitumor activity while widening the therapeutic index of IL12 treatment and is currently being evaluated in a phase I clinical trial.

Solvation free energy of the peptide group: its model dependence and implications for the additive-transfer free-energy model of protein stability.

The group-additive decomposition of the unfolding free energy of a protein in an osmolyte solution relative to that in water poses a fundamental paradox: whereas the decomposition describes the experimental results rather well, theory suggests that a group-additive decomposition of free energies is, in general, not valid. In a step toward resolving this paradox, here we study the peptide-group transfer free energy. We calculate the vacuum-to-solvent (solvation) free energies of (Gly)n and cyclic diglycine (cGG) and analyze the data according to experimental protocol. The solvation free energies of (Gly)n are linear in n, suggesting group additivity. However, the slope interpreted as the free energy of a peptide unit differs from that for cGG scaled by a factor of half, emphasizing the context dependence of solvation. However, the water-to-osmolyte transfer free energies of the peptide unit are relatively independent of the peptide model, as observed experimentally. To understand these observations, a way to assess the contribution to the solvation free energy of solvent-mediated correlation between distinct groups is developed. We show that linearity of solvation free energy with n is a consequence of uniformity of the correlation contributions, with apparent group-additive behavior in the water-to-osmolyte transfer arising due to their cancellation. Implications for inferring molecular mechanisms of solvent effects on protein stability on the basis of the group-additive transfer model are suggested.

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, storage, and shipping. In this work, we have studied the conformational dynamics of human IgG1κ 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 transient 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.

Hydrophilic Interactions Dominate the Inverse Temperature Dependence of Polypeptide Hydration Free Energies Attributed to Hydrophobicity.

We address the association of the hydrophobic driving forces in protein folding with the inverse temperature dependence of protein hydration, wherein stabilizing hydration effects strengthen with increasing temperature in a physiological range. All-atom calculations of the free energy of hydration of aqueous deca-alanine conformers, holistically including backbone and side-chain interactions together, show that attractive peptide-solvent interactions and the thermal expansion of the solvent dominate the inverse temperature signatures that have been interpreted traditionally as the hydrophobic stabilization of proteins in aqueous solution. Equivalent calculations on a methane solute are also presented as a benchmark for comparison. The present study calls for a reassessment of the forces that stabilize folded protein conformations in aqueous solutions and of the additivity of hydrophobic/hydrophilic contributions.

Rational optimization of a monoclonal antibody for simultaneous improvements in its solution properties and biological activity.

Developability considerations should be integrated with lead engineering of antibody drug candidates in interest of their cost effective translations into medicines. To explore feasibility of this imperative, we have performed rational mutagenesis studies on a monoclonal antibody (MAB1) whose development was discontinued owing to manufacturability hurdles. Seven computationally designed variants of MAB1 containing single point (V44K, E59S, E59T and E59Y) and double (V44KE59S, V44KE59T and V44KE59Y) mutations in its light chain were produced in Chinese Hamster Ovary (CHO) cells and purified by using platform processes employed during commercial scale production of monoclonal antibodies. MAB1 and its variants were formulated in the same platform buffer and subjected to a battery of experiments to assess their solution behaviors, and biological activities. Five of the seven (71%) variants of MAB1 demonstrated improved biophysical attributes in multiple experimental testings. Contrary to the commonly expressed reservations about potential biological activity loss upon developability optimizations, the improvements in solution behavior of MAB1 also increased its biological activity up to ~180%. In particular, concentrate-ability and apparent solubility of V44KE59S improved to ~150% and ~160%, respectively. Its diffusion interaction parameter (kD) reduced to 28% and viscosity at ~100 mg/ml decreased to less than half of the corresponding values for MAB1. V44KE59S is also slightly more active and its transfections in CHO cells were more productive. It also degraded slower than MAB1 in three month long 25°C and 40°C formulation stability studies. These results open doors to an exciting realm of structure-based biologic drug design where developability and biological activity can be simultaneously optimized at the molecular engineering stages.

Intramolecular Interactions Overcome Hydration to Drive the Collapse Transition of Gly15.

Simulations and experiments show oligo-glycines, polypeptides lacking any side chains, can collapse in water. We assess the hydration thermodynamics of this collapse by calculating the hydration free energy at each of the end points of the reaction coordinate, here taken as the end-to-end distance (r) in the chain. To examine the role of the various conformations for a given r, we study the conditional distribution, P(Rg|r), of the radius of gyration for a given value of r. The free energy change versus Rg, -kBT ln P(Rg|r), is found to vary more gently compared to the corresponding variation in the excess hydration free energy. Using this observation within a multistate generalization of the potential distribution theorem, we calculate a tight upper bound for the hydration free energy of the peptide for a given r. On this basis, we find that peptide hydration greatly favors the expanded state of the chain, despite primitive hydrophobic effects favoring chain collapse. The net free energy of collapse is seen to be a delicate balance between opposing intrapeptide and hydration effects, with intrapeptide contributions favoring collapse.

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 VH, VL, 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.

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

Effective translation of breakthrough discoveries into innovative products in the clinic requires proactive mitigation or elimination of several drug development challenges. These challenges can vary depending upon the type of drug molecule. In the case of therapeutic antibody candidates, a commonly encountered challenge is high viscosity of the concentrated antibody solutions. Concentration-dependent viscosity behaviors of mAbs and other biologic entities may depend on pairwise and higher-order intermolecular interactions, non-native aggregation, and concentration-dependent fluctuations of various antibody regions. This article reviews our current understanding of molecular origins of viscosity behaviors of antibody solutions. We discuss general strategies and guidelines to select low viscosity candidates or optimize lead candidates for lower viscosity at early drug discovery stages. Moreover, strategies for formulation optimization and excipient design are also presented for candidates already in advanced product development stages. Potential future directions for research in this field are also explored.

Importance of Hydrophilic Hydration and Intramolecular Interactions in the Thermodynamics of Helix-Coil Transition and Helix-Helix Assembly in a Deca-Alanine Peptide.

For a model deca-alanine peptide the cavity (ideal hydrophobic) contribution to hydration favors the helix state over extended states and the paired helix bundle in the assembly of two helices. The energetic contributions of attractive protein-solvent interactions are separated into quasi-chemical components consisting of a short-range part arising from interactions with solvent in the first hydration shell and the remaining long-range part that is well described by a Gaussian. In the helix-coil transition, short-range attractive protein-solvent interactions outweigh hydrophobic hydration and favor the extended coil states. Analysis of enthalpic effects shows that it is the favorable hydration of the peptide backbone that favors the unfolded state. Protein intramolecular interactions favor the helix state and are decisive in favoring folding. In the pairing of two helices, the cavity contribution outweighs the short-range attractive protein-water interactions. However, long-range, protein-solvent attractive interactions can either enhance or reverse this trend depending on the mutual orientation of the helices. In helix-helix assembly, change in enthalpy arising from change in attractive protein-solvent interactions favors disassembly. In helix pairing as well, favorable protein intramolecular interactions are found to be as important as hydration effects. Overall, hydrophilic protein-solvent interactions and protein intramolecular interactions are found to play a significant role in the thermodynamics of folding and assembly in the system studied.

Conditional solvation thermodynamics of isoleucine in model peptides and the limitations of the group-transfer model.

The hydration thermodynamics of the amino acid X relative to the reference G (glycine) or the hydration thermodynamics of a small-molecule analog of the side chain of X is often used to model the contribution of X to protein stability and solution thermodynamics. We consider the reasons for successes and limitations of this approach by calculating and comparing the conditional excess free energy, enthalpy, and entropy of hydration of the isoleucine side chain in zwitterionic isoleucine, in extended penta-peptides, and in helical deca-peptides. Butane in gauche conformation serves as a small-molecule analog for the isoleucine side chain. Parsing the hydrophobic and hydrophilic contributions to hydration for the side chain shows that both of these aspects of hydration are context-sensitive. Furthermore, analyzing the solute-solvent interaction contribution to the conditional excess enthalpy of the side chain shows that what is nominally considered a property of the side chain includes entirely nonobvious contributions of the background. The context-sensitivity of hydrophobic and hydrophilic hydration and the conflation of background contributions with energetics attributed to the side chain limit the ability of a single scaling factor, such as the fractional solvent exposure of the group in the protein, to map the component energetic contributions of the model-compound data to their value in the protein. But ignoring the origin of cancellations in the underlying components the group-transfer model may appear to provide a reasonable estimate of the free energy for a given error tolerance.