Surface patches induce nonspecific binding and phase separation of antibodies.

Nonspecific interactions are a key challenge in the successful development of therapeutic antibodies. The tendency for nonspecific binding of antibodies is often difficult to reduce by rational design, and instead, it is necessary to rely on comprehensive screening campaigns. To address this issue, we performed a systematic analysis of the impact of surface patch properties on antibody nonspecificity using a designer antibody library as a model system and single-stranded DNA as a nonspecificity ligand. Using an in-solution microfluidic approach, we find that the antibodies tested bind to single-stranded DNA with affinities as high as KD = 1 µM. We show that DNA binding is driven primarily by a hydrophobic patch in the complementarity-determining regions. By quantifying the surface patches across the library, the nonspecific binding affinity is shown to correlate with a trade-off between the hydrophobic and total charged patch areas. Moreover, we show that a change in formulation conditions at low ionic strengths leads to DNA-induced antibody phase separation as a manifestation of nonspecific binding at low micromolar antibody concentrations. We highlight that phase separation is driven by a cooperative electrostatic network assembly mechanism of antibodies with DNA, which correlates with a balance between positive and negative charged patches. Importantly, our study demonstrates that both nonspecific binding and phase separation are controlled by the size of the surface patches. Taken together, these findings highlight the importance of surface patches and their role in conferring antibody nonspecificity and its macroscopic manifestation in phase separation.

Knowledge-Based Strategy to Improve Ligand Pose Prediction Accuracy for Lead Optimization.

Accurately predicting how a small molecule binds to its target protein is an essential requirement for structure-based drug design (SBDD) efforts. In structurally enabled medicinal chemistry programs, binding pose prediction is often applied to ligands after a related compound's crystal structure bound to the target protein has been solved. In this article, we present an automated pose prediction protocol that makes extensive use of existing X-ray ligand information. It uses spatial restraints during docking based on maximum common substructure (MCS) overlap between candidate molecule and existing X-ray coordinates of the related compound. For a validation data set of 8784 docking runs, our protocol's pose prediction accuracy (80-82%) is almost two times higher than that of one unbiased docking method software (43%). To demonstrate the utility of this protocol in a project setting, we show its application in a chronological manner for a number of internal drug discovery efforts. The accuracy and applicability of this algorithm (>70% of cases) to medicinal chemistry efforts make this the approach of choice for pose prediction in lead optimization programs.

Structure-Based Optimization of Antibody-Based Biotherapeutics for Improved Developability: A Practical Guide for Molecular Modelers.

A great effort to avoid known developability risks is now more often being made earlier during the lead candidate discovery and optimization phase of biotherapeutic drug development. Predictive computational strategies, used in the early stages of antibody discovery and development, to mitigate the risk of late-stage failure of antibody candidates, are highly valuable. Various structure-based methods exist for accurately predicting properties critical to developability, and, in this chapter, we discuss the history of their development and demonstrate how they can be used to filter large sets of candidates arising from target affinity screening and to optimize lead candidates for developability. Methods for modeling antibody structures from sequence and detecting post-translational modifications and chemical degradation liabilities are also discussed.

Structural Considerations in Affinity Maturation of Antibody-Based Biotherapeutic Candidates.

Affinity maturation is an important stage in biologic drug discovery as is the natural process of generating an immune response inside the human body. In this chapter, we describe in silico approaches to affinity maturation via a worked example. Both advantages and limitations of the computational methods used are critically examined. Furthermore, construction of affinity maturation libraries and how their outputs might be implemented in an experimental setting are also described. It should be noted that structure-based design of biologic drugs is an emerging field and the tools currently available require further development. Furthermore, there are no standardized structure-based strategies yet for antibody affinity maturation as this research relies heavily on scientific logic as well as creative intuition.

Antibody modeling assessment.

A blinded study to assess the state of the art in three-dimensional structure modeling of the variable region (Fv) of antibodies was conducted. Nine unpublished high-resolution x-ray Fab crystal structures covering a wide range of antigen-binding site conformations were used as benchmark to compare Fv models generated by four structure prediction methodologies. The methodologies included two homology modeling strategies independently developed by CCG (Chemical Computer Group) and Accerlys Inc, and two fully automated antibody modeling servers: PIGS (Prediction of ImmunoGlobulin Structure), based on the canonical structure model, and Rosetta Antibody Modeling, based on homology modeling and Rosetta structure prediction methodology. The benchmark structure sequences were submitted to Accelrys and CCG and a set of models for each of the nine antibody structures were generated. PIGS and Rosetta models were obtained using the default parameters of the servers. In most cases, we found good agreement between the models and x-ray structures. The average rmsd (root mean square deviation) values calculated over the backbone atoms between the models and structures were fairly consistent, around 1.2 Å. Average rmsd values of the framework and hypervariable loops with canonical structures (L1, L2, L3, H1, and H2) were close to 1.0 Å. H3 prediction yielded rmsd values around 3.0 Å for most of the models. Quality assessment of the models and the relative strengths and weaknesses of the methods are discussed. We hope this initiative will serve as a model of scientific partnership and look forward to future antibody modeling assessments.

Structure-based charge calculations for predicting isoelectric point, viscosity, clearance, and profiling antibody therapeutics.

The effect of hydrophobicity on antibody aggregation is well understood, and it has been shown that charge calculations can be useful for high-concentration viscosity and pharmacokinetic (PK) clearance predictions. In this work, structure-based charge descriptors are evaluated for their predictive performance on recently published antibody pI, viscosity, and clearance data. From this, we devised four rules for therapeutic antibody profiling which address developability issues arising from hydrophobicity and charged-based solution behavior, PK, and the ability to enrich for those that are approved by the U.S. Food and Drug Administration. Differences in strategy for optimizing the solution behavior of human IgG1 antibodies versus the IgG2 and IgG4 isotypes and the impact of pH alterations in formulation are discussed.

Identifying biophysical assays and in silico properties that enrich for slow clearance in clinical-stage therapeutic antibodies.

Understanding the pharmacokinetic (PK) properties of a drug, such as clearance, is a crucial step for evaluating efficacy. The PK of therapeutic antibodies can be complex and is influenced by interactions with the target, Fc-receptors, anti-drug antibodies, and antibody intrinsic factors. A growing body of literature has linked biophysical properties of antibodies, particularly nonspecific-binding propensity, hydrophobicity and charged regions to rapid clearance in preclinical species and selected human PK studies. A clear understanding of the connection between biophysical properties and their impact on PK would allow for early selection and optimization of antibodies and reduce costly attrition during clinical trials due to sub-optimal human clearance. Due to the difficulty in obtaining large and unbiased human PK data, previous studies have focused mostly on preclinical PK. For this study, we obtained and curated the most comprehensive clinical PK dataset to date and calculated accurate estimates of linear clearance for 64 monoclonal antibodies ranging from investigational candidates in Phase 2 trials to marketed products. This allows for the first time a deep analysis of the influence of biophysical and sequence-based in silico properties directly on human clearance. We use statistical analysis and a Random Forest classifier to identify properties that have the greatest influence in our dataset. Our findings indicate that in vitro poly-specificity assay and in silico estimated isoelectric point can discriminate fast and slow clearing antibodies, extending previous observations on preclinical clearance. This provides a simple yet powerful approach to select antibodies with desirable PK during early-stage screening.

Beyond affinity: selection of antibody variants with optimal biophysical properties and reduced immunogenicity from mammalian display libraries.

The early phase of protein drug development has traditionally focused on target binding properties leading to a desired mode of therapeutic action. As more protein therapeutics pass through the development pipeline; however, it is clear that non-optimal biophysical properties can emerge, particularly as proteins are formulated at high concentrations, causing aggregation or polyreactivity. Such late-stage "developability" problems can lead to delay or failure in traversing the development process. Aggregation propensity is also correlated with increased immunogenicity, resulting in expensive, late-stage clinical failures. Using nucleases-directed integration, we have constructed large mammalian display libraries where each cell contains a single antibody gene/cell inserted at a single locus, thereby achieving transcriptional normalization. We show a strong correlation between poor biophysical properties and display level achieved in mammalian cells, which is not replicated by yeast display. Using two well-documented examples of antibodies with poor biophysical characteristics (MEDI-1912 and bococizumab), a library of variants was created based on surface hydrophobic and positive charge patches. Mammalian display was used to select for antibodies that retained target binding and permitted increased display level. The resultant variants exhibited reduced polyreactivity and reduced aggregation propensity. Furthermore, we show in the case of bococizumab that biophysically improved variants are less immunogenic than the parental molecule. Thus, mammalian display helps to address multiple developability issues during the earliest stages of lead discovery, thereby significantly de-risking the future development of protein drugs.

In silico identification of anthropogenic chemicals as ligands of zebrafish sex hormone binding globulin.

Anthropogenic compounds with the capacity to interact with the steroid-binding site of sex hormone binding globulin (SHBG) pose health risks to humans and other vertebrates including fish. Building on studies of human SHBG, we have applied in silico drug discovery methods to identify potential binders for SHBG in zebrafish (Danio rerio) as a model aquatic organism. Computational methods, including; homology modeling, molecular dynamics simulations, virtual screening, and 3D QSAR analysis, successfully identified 6 non-steroidal substances from the ZINC chemical database that bind to zebrafish SHBG (zfSHBG) with low-micromolar to nanomolar affinities, as determined by a competitive ligand-binding assay. We also screened 80,000 commercial substances listed by the European Chemicals Bureau and Environment Canada, and 6 non-steroidal hits from this in silico screen were tested experimentally for zfSHBG binding. All 6 of these compounds displaced the [(3)H]5alpha-dihydrotestosterone used as labeled ligand in the zfSHBG screening assay when tested at a 33 microM concentration, and 3 of them (hexestrol, 4-tert-octylcatechol, and dihydrobenzo(a)pyren-7(8H)-one) bind to zfSHBG in the micromolar range. The study demonstrates the feasibility of large-scale in silico screening of anthropogenic compounds that may disrupt or highjack functionally important protein:ligand interactions. Such studies could increase the awareness of hazards posed by existing commercial chemicals at relatively low cost.

An updated steroid benchmark set and its application in the discovery of novel nanomolar ligands of sex hormone-binding globulin.

A benchmark data set of steroids with known affinity for sex hormone-binding globulin (SHBG) has been widely used to validate popular molecular field-based QSAR techniques. We have expanded the data set by adding a number of nonsteroidal SHBG ligands identified both from the literature and in our previous experimental studies. This updated molecular set has been used herein to develop 4D QSAR models based on "inductive" descriptors and to gain insight into the molecular basis of protein-ligand interactions. Molecular alignment was generated by means of docking active compounds into the active site of the SHBG. Surprisingly, the alignment of the benchmark steroids contradicted the classical ligand-based alignment utilized in previous CoMFA and CoMSIA models yet afforded models with higher statistical significance and predictive power. The resulting QSAR models combined with CoMFA and CoMSiA models as well as structure-based virtual screening allowed discovering several low-micromolar to nanomolar nonsteroidal inhibitors for human SHBG.

Homology modeling and structure-based design improve hydrophobic interaction chromatography behavior of integrin binding antibodies.

Monoclonal antibody (mAb) candidates from high-throughput screening or binding affinity optimization often contain mutations leading to liabilities for further development of the antibody, such as aggregation-prone regions and lack of solubility. In this work, we optimized a candidate integrin α11-binding mAb for developability using molecular modeling, rational design, and hydrophobic interaction chromatography (HIC). A homology model of the parental mAb Fv region was built, and this revealed hydrophobic patches on the surface of the complementarity-determining region loops. A series of 97 variants of the residues primarily responsible for the hydrophobic patches were expressed and their HIC retention times (RT) were measured. As intended, many of the computationally designed variants reduced the HIC RT compared to the parental mAb, and mutating residues that contributed most to hydrophobic patches had the greatest effect on HIC RT. A retrospective analysis was then performed where 3-dimentional protein property descriptors were evaluated for their ability to predict HIC RT using the current series of mAbs. The same descriptors were used to train a simple multi-parameter protein quantitative structure-property relationship model on this data, producing an improved correlation. We also extended this analysis to recently published HIC data for 137 clinical mAb candidates as well as 31 adnectin variants, and found that the surface area of hydrophobic patches averaged over a molecular dynamics sample consistently correlated to the experimental data across a diverse set of biotherapeutics.