chemalot and chemalot_knime: Command line programs as workflow tools for drug discovery

Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Utilizing Molecular Simulations to Examine Nanosuspension Stability

Drug nanosuspensions offer a promising approach to improve bioavailability for poorly soluble drug candidates. Such formulations often necessitate the inclusion of an excipient to stabilize the drug nanoparticles. However, the rationale for the choice of the correct excipient for a given drug candidate remains unclear. To gain molecular insight into formulation design, this work first utilizes a molecular dynamics simulation to computationally investigate drug-excipient interactions for a number of combinations that have been previously studied experimentally. We find that hydrophobic interactions drive excipient adsorption to drug nanoparticles and that the fraction of polar surface area serves as a predictor for experimental measurements of nanosuspension stability. To test these ideas prospectively, we applied our model to an uncharacterized drug compound, GDC-0810. Our simulations predicted that a salt form of GDC-0810 would lead to more stable nanosuspensions than the neutral form; therefore, we tested the stability of salt GDC-0810 nanosuspensions and found that the salt form readily formed nanosuspensions even without the excipient. To avoid computationally expensive simulations in the future, we extended our model by showing that simple, two-dimensional properties of single drug molecules can be used to rationalize nanosuspension designs without simulations. In all, our work demonstrates how computational tools can provide molecular insight into drug-excipient interactions and aid in rational formulation design.

Identification of Macrocyclic Peptide Families from Combinatorial Libraries Containing Noncanonical Amino Acids Using Cheminformatics and Bioinformatics Inspired Clustering

In the past decade, macrocyclic peptides gained increasing interest as a new therapeutic modality to tackle intracellular and extracellular therapeutic targets that had been previously classified as "undruggable". Several technological advances have made discovering macrocyclic peptides against these targets possible: 1) the inclusion of noncanonical amino acids (NCAAs) into mRNA display, 2) increased availability of next generation sequencing (NGS), and 3) improvements in rapid peptide synthesis platforms. This type of directed-evolution based screening can produce large numbers of potential hit sequences given that DNA sequencing is the functional output of this platform. The current standard for selecting hit peptides from these selections for downstream follow-up relies on the frequency counting and sorting of unique peptide sequences which can result in the generation of false negatives due to technical reasons including low translation efficiency or other experimental factors. To overcome our inability to detect weakly enriched peptide sequences among our large data sets, we wanted to develop a clustering method that would enable the identification of peptide families. Unfortunately, utilizing traditional clustering algorithms, such as ClustalW, is not possible for this technology due to the incorporation of NCAAs in these libraries. Therefore, we developed a new atomistic clustering method with a Pairwise Aligned Peptide (PAP) chemical similarity metric to perform sequence alignments and identify macrocyclic peptide families. With this method, low enriched peptides, including isolated sequences (singletons), can now be clustered into families providing a comprehensive analysis of NGS data resulting from macrocycle discovery selections. Additionally, upon identification of a hit peptide with the desired activity, this clustering algorithm can be used to identify derivatives from the initial data set for structure-activity relationship (SAR) analysis without requiring additional selection experiments.

Toward Quantitative Models in Safety Assessment: A Case Study to Show Impact of Dose-Response Inference on hERG Inhibition Models

Due to challenges with historical data and the diversity of assay formats, in silico models for safety-related endpoints are often based on discretized data instead of the data on a natural continuous scale. Models for discretized endpoints have limitations in usage and interpretation that can impact compound design. Here, we present a consistent data inference approach, exemplified on two data sets of Ether-à-go-go-Related Gene (hERG) K+ inhibition data, for dose-response and screening experiments that are generally applicable for in vitro assays. hERG inhibition has been associated with severe cardiac effects and is one of the more prominent safety targets assessed in drug development, using a wide array of in vitro and in silico screening methods. In this study, the IC₅₀ for hERG inhibition is estimated from diverse historical proprietary data. The IC₅₀ derived from a two-point proprietary screening data set demonstrated high correlation (R = 0.98, MAE = 0.08) with IC_50s derived from six-point dose-response curves. Similar IC₅₀ estimation accuracy was obtained on a public thallium flux assay data set (R = 0.90, MAE = 0.2). The IC₅₀ data were used to develop a robust quantitative model. The model's MAE (0.47) and R² (0.46) were on par with literature statistics and approached assay reproducibility. Using a continuous model has high value for pharmaceutical projects, as it enables rank ordering of compounds and evaluation of compounds against project-specific inhibition thresholds. This data inference approach can be widely applicable to assays with quantitative readouts and has the potential to impact experimental design and improve model performance, interpretation, and acceptance across many standard safety endpoints.

Chemical space docking enables large-scale structure-based virtual screening to discover ROCK1 kinase inhibitors

With the ever-increasing number of synthesis-on-demand compounds for drug lead discovery, there is a great need for efficient search technologies. We present the successful application of a virtual screening method that combines two advances: (1) it avoids full library enumeration (2) products are evaluated by molecular docking, leveraging protein structural information. Crucially, these advances enable a structure-based technique that can efficiently explore libraries with billions of molecules and beyond. We apply this method to identify inhibitors of ROCK1 from almost one billion commercially available compounds. Out of 69 purchased compounds, 27 (39%) have K_i values < 10 µM. X-ray structures of two leads confirm their docked poses. This approach to docking scales roughly with the number of reagents that span a chemical space and is therefore multiple orders of magnitude faster than traditional docking.

Comparison of logP and logD correction models trained with public and proprietary data sets

In drug discovery, partition and distribution coefficients, logP and logD for octanol/water, are widely used as metrics of the lipophilicity of molecules, which in turn have a strong influence on the bioactivity and bioavailability of potential drugs. There are a variety of established methods, mostly fragment or atom-based, to calculate logP while logD prediction generally relies on calculated logP and pKa for the estimation of neutral and ionized populations at a given pH. Algorithms such as ClogP have limitations generally leading to systematic errors for chemically related molecules while pKa estimation is generally more difficult due to the interplay of electronic, inductive and conjugation effects for ionizable moieties. We propose an integrated machine learning QSAR modeling approach to predict logD by training the model with experimental data while using ClogP and pKa predicted by commercial software as model descriptors. By optimizing the loss function for the ClogD calculated by the software, we build a correction model that incorporates both descriptors from the software and available experimental logD data. Additionally, we calculate logP from the logD model using the software predicted pKa's. Here, we have trained models using publicly or commercial available logD data to show that this approach can improve on commercial software predictions of lipophilicity. When applied to other logD data sets, this approach extends the domain of applicability of logD and logP predictions over commercial software. Performance of these models favorably compare with models built with a larger set of proprietary logD data.

Ligand Strain Energy in Large Library Docking

While small molecule internal strain is crucial to molecular docking, using it in evaluating ligand scores has remained elusive. Here, we investigate a technique that calculates strain using relative torsional populations in the Cambridge Structural Database, enabling fast precalculation of these energies. In retrospective studies of large docking screens of the dopamine D4 receptor and of AmpC β-lactamase, where close to 600 docking hits were tested experimentally, including such strain energies improved hit rates by preferentially reducing the ranks of strained high-scoring decoy molecules. In a 40-target subset of the DUD-E benchmark, we found two thresholds that usefully distinguished between ligands and decoys: one based on the total strain energy of the small molecules and another based on the maximum strain allowed for any given torsion within them. Using these criteria, about 75% of the benchmark targets had improved enrichment after strain filtering. Relying on precalculated population distributions, this approach is rapid, taking less than 0.04 s to evaluate a conformation on a standard core, making it pragmatic for precalculating strain in even ultralarge libraries. Since it is scoring function agnostic, it may be useful to multiple docking approaches; it is openly available at http://tldr.docking.org.

Novel Homodimer Metabolites of GDC-0994 via Cytochrome P450-Catalyzed Radical Coupling

Two novel homodimer metabolites were identified in rat samples collected during the in vivo study of GDC-0994. In this study, we investigated the mechanism of the formation of these metabolites. We generated and isolated the dimer metabolites using a biomimetic oxidation system for NMR structure elucidation to identify a symmetric dimer formed via carbon-carbon bond between two pyrazoles and an asymmetric dimer formed via an aminopyrazole-nitrogen to pyrazole-carbon bond. In vitro experiments demonstrated formation of these dimers was catalyzed by cytochrome P450 enzymes (P450s) with CYP3A4/5 being the most efficient. Using density functional theory, we determined these metabolites share a mechanism of formation, initiated by an N-H hydrogen atom abstraction by the catalytically active iron-oxo of P450s. Molecular modeling studies also show these dimer metabolites fit in the CYP3A4 binding site in low energy conformations with minimal protein rearrangement. Collectively, the results of these experiments suggest that formation of these two homodimer metabolites is mediated by CYP3A, likely involving activation of two GDC-0994 molecules by a single P450 enzyme and proceeding through a radical coupling mechanism. SIGNIFICANCE STATEMENT: These studies identified structures and enzymology for two distinct homodimer metabolites and indicate a novel biotransformation reaction mediated by CYP3A. In it, two molecules may bind within the active site and combine through radical coupling. The mechanism of dimerization was elucidated using density functional theory computations and supported by molecular modeling.

A path to next-generation reproducibility in cheminformatics

Currently, the submission guidelines for the Journal of Cheminformatics say it will "only publish research or software that is entirely reproducible by third parties." They go on to specify that being reproducible means that anything essential to the conclusion of the paper be freely accessible and states that source code must be provided. I submit that this definition of reproducibility is too narrow-that a cheminformatics method can only truly be replicated by reimplementing it from a detailed, step-by-step high-level description to determine how reliably the algorithm per se does what it is intended to do.

A chemoinformatics approach for the characterization of hybrid nanomaterials: safer and efficient design perspective

In this study, photocatalytic properties and in vitro cytotoxicity of 29 TiO2-based multi-component nanomaterials (i.e., hybrids of more than two composition types of nanoparticles) were evaluated using a combination of the experimental testing and supervised machine learning modeling. TiO2-based multi-component nanomaterials with metal clusters of silver, and their mixtures with gold, palladium, and platinum were successfully synthesized. Two activities, photocatalytic activity and cytotoxicity, were studied. A novel cheminformatic approach was developed and applied for the computational representation of the photocatalytic activity and cytotoxicity effect. In this approach, features of investigated TiO2-based hybrid nanomaterials were reflected by a series of novel additive descriptors for hybrid and hybrid nanostructures (denoted as "hybrid nanosctructure descriptors"). These descriptors are based on quantum chemical calculations and the Smoluchowski equation. The obtained experimental data and calculated hybrid-nanostructure descriptors were used to develop novel predictive Quantitative Structure-Activity Relationship computational models (called "nano-QSARmix"). The proposed modeling approach is an initial step in the understanding of the relationships between physicochemical properties of hybrid nanoparticles, their toxicity, and photochemical activity under UV-vis irradiation. Acquired knowledge supports the safe-by-design approaches relevant to the development of efficient hybrid nanomaterials with reduced hazardous effects.

Computational Support of Medicinal Chemistry in Industrial Settings

The practice of computational chemistry in an industrial setting poses unique opportunities and challenges. Industrial computational chemists must manage large amounts of data, master modeling software, write scripts to perform custom calculations, and stay abreast of scientific advances in the field. Just as importantly, because computational chemists are full partners in the drug discovery effort at companies, in order to influence and streamline the drug discovery process, they must communicate effectively with medicinal chemists and other scientists to deliver results of their calculations in a timely fashion. The skills necessary to play this role require education that emphasizes a combination of chemistry, programming, and communication skills. Professors are encouraged to incorporate such training in their curriculum.