Dataset Information

Computational Characterization of Antibody-Excipient Interactions for Rational Excipient Selection Using the Site Identification by Ligand Competitive Saturation-Biologics Approach.

ABSTRACT: Protein therapeutics typically require a concentrated protein formulation, which can lead to self-association and/or high viscosity due to protein-protein interaction (PPI). Excipients are often added to improve stability, bioavailability, and manufacturability of the protein therapeutics, but the selection of excipients often relies on trial and error. Therefore, understanding the excipient-protein interaction and its effect on non-specific PPI is important for rational selection of formulation development. In this study, we validate a general workflow based on the site identification by ligand competitive saturation (SILCS) technology, termed SILCS-Biologics, that can be applied to protein therapeutics for rational excipient selection. The National Institute of Standards and Technology monoclonal antibody (NISTmAb) reference along with the CNTO607 mAb is used as model antibody proteins to examine PPIs, and NISTmAb was used to further examine excipient-protein interactions, in silico. Metrics from SILCS include the distribution and predicted affinity of excipients, buffer interactions with the NISTmAb Fab, and the relation of the interactions to predicted PPI. Comparison with a range of experimental data showed multiple SILCS metrics to be predictive. Specifically, the number of favorable sites to which an excipient binds and the number of sites to which an excipient binds that are involved in predicted PPIs correlate with the experimentally determined viscosity. In addition, a combination of the number of binding sites and the predicted binding affinity is indicated to be predictive of relative protein stability. Comparison of arginine, trehalose, and sucrose, all of which give the highest viscosity in combination with analysis of B₂₂ and k_D and the SILCS metrics, indicates that higher viscosities are associated with a low number of predicted binding sites, with lower binding affinity of arginine leading to its anomalously high impact on viscosity. The present study indicates the potential for the SILCS-Biologics approach to be of utility in the rational design of excipients during biologics formulation.

SUBMITTER: Jo S

PROVIDER: S-EPMC7606568 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Predicting relative protein-ligand binding affinities is a central pillar of lead optimization efforts in structure-based drug design. The site identification by ligand competitive saturation (SILCS) methodology is based on functional group affinity patterns in the form of free energy maps that may be used to compute protein-ligand binding poses and affinities. Presented are results obtained from the SILCS methodology for a set of eight target proteins as reported originally in Wang et al. (J. Am. Chem. Soc., 2015, 137, 2695-2703) using free energy perturbation (FEP) methods in conjunction with enhanced sampling and cycle closure corrections. These eight targets have been subsequently studied by many other authors to compare the efficacy of their method while comparing with the outcomes of Wang et al. In this work, we present results for a total of 407 ligands on the eight targets and include specific analysis on the subset of 199 ligands considered previously. Using the SILCS methodology we can achieve an average accuracy of up to 77% and 74% when considering the eight targets with their 199 and 407 ligands, respectively, for rank-ordering ligand affinities as calculated by the percent correct metric. This accuracy increases to 82% and 80%, respectively, when the SILCS atomic free energy contributions are optimized using a Bayesian Markov-chain Monte Carlo approach. We also report other metrics including Pearson's correlation coefficient, Pearlman's predictive index, mean unsigned error, and root mean square error for both sets of ligands. The results obtained for the 199 ligands are compared with the outcomes of Wang et al. and other published works. Overall, the SILCS methodology yields similar or better-quality predictions without a priori need for known ligand orientations in terms of the different metrics when compared to current FEP approaches with significant computational savings while additionally offering quantitative estimates of individual atomic contributions to binding free energies. These results further validate the SILCS methodology as an accurate, computationally efficient tool to support lead optimization and drug discovery.

Project description:Drug-induced cardiotoxicity is a potentially lethal and yet one of the most common side effects with the drugs in clinical use. Most of the drug-induced cardiotoxicity is associated with an off-target pharmacological blockade of K+ currents carried out by the cardiac Human-Ether-a-go-go-Related (hERG1) potassium channel. There is a compulsory preclinical stage safety assessment for the hERG1 blockade for all classes of drugs, which adds substantially to the cost of drug development. The availability of a high-resolution cryogenic electron microscopy (cryo-EM) structure for the channel in its open/depolarized state solved in 2017 enabled the application of molecular modeling for rapid assessment of drug blockade by molecular docking and simulation techniques. More importantly, if successful, in silico methods may allow a path to lead-compound salvaging by mapping out key block determinants. Here, we report the blind application of the site identification by the ligand competitive saturation (SILCS) protocol to map out druggable/regulatory hotspots in the hERG1 channel available for blockers and activators. The SILCS simulations use small solutes representative of common functional groups to sample the chemical space for the entire protein and its environment using all-atom simulations. The resulting chemical maps, FragMaps, explicitly account for receptor flexibility, protein-fragment interactions, and fragment desolvation penalty allowing for rapid ranking of potential ligands as blockers or nonblockers of hERG1. To illustrate the power of the approach, SILCS was applied to a test set of 55 blockers with diverse chemical scaffolds and pIC50 values measured under uniform conditions. The original SILCS model was based on the all-atom modeling of the hERG1 channel in an explicit lipid bilayer and was further augmented with a Bayesian-optimization/machine-learning (BML) stage employing an independent literature-derived training set of 163 molecules. BML approach was used to determine weighting factors for the FragMaps contributions to the scoring function. pIC50 predictions from the combined SILCS/BML approach to the 55 blockers showed a Pearson correlation (PC) coefficient of >0.535 relative to the experimental data. SILCS/BML model was shown to yield substantially improved performance as compared to commonly used rigid and flexible molecular docking methods for a well-established cohort of hERG1 blockers, where no correlation with experimental data was recorded. SILCS/BML results also suggest that a proper weighting of protonation states of common blockers present at physiological pH is essential for accurate predictions of blocker potency. The precalculated and optimized SILCS FragMaps can now be used for the rapid screening of small molecules for their cardiotoxic potential as well as for exploring alternative binding pockets in the hERG1 channel with applications to the rational design of activators.

Dataset Information

Computational Characterization of Antibody-Excipient Interactions for Rational Excipient Selection Using the Site Identification by Ligand Competitive Saturation-Biologics Approach.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets