Descriptor-Driven de Novo Design Algorithms for DOCK6 Using RDKit.

Structure-based methods that employ principles of de novo design can be used to construct small organic molecules from scratch using pre-existing fragment libraries to sample chemical space and are an important class of computational algorithms for drug-lead discovery. Here, we present a powerful new design method for DOCK6 that employs a Descriptor-Driven De Novo strategy (termed D3N) in which user-defined cheminformatics descriptors (and their target ranges) are calculated at each layer of growth using the open-source toolkit RDKit. The objective is to tailor ligand growth toward desirable regions of chemical space. The approach was extensively validated through: (1) comparison of cheminformatics descriptors computed using the new DOCK6/RDKit interface versus the standard Python/RDKit installation, (2) examination of descriptor distributions generated using D3N growth under different conditions (target ranges and environments), and (3) construction of ligands with very tight (pinpoint) descriptor ranges using clinically relevant compounds as a reference. Our testing confirms that the new DOCK6/RDKit integration is robust, showcases how the new D3N routines can be used to direct sampling around user-defined chemical spaces, and highlights the utility of on-the-fly descriptor calculations for ligand design to important drug targets.

Neutral ceramidase-active site inhibitor chemotypes and binding modes.

Ceramides impact a diverse array of biological functions and have been implicated in disease pathogenesis. The enzyme neutral ceramidase (nCDase) is a zinc-containing hydrolase and mediates the metabolism of ceramide to sphingosine (Sph), both in cells and in the intestinal lumen. nCDase inhibitors based on substrate mimetics, for example C6-urea ceramide, have limited potency, aqueous solubility, and micelle-free fraction. To identify non-ceramide mimetic nCDase inhibitors, hit compounds from an HTS campaign were evaluated in biochemical, cell based and in silico modeling approaches. A majority of small molecule nCDase inhibitors contained pharmacophores capable of zinc interaction but retained specificity for nCDase over zinc-containing acid and alkaline ceramidases, as well as matrix metalloprotease-3 and histone deacetylase-1. nCDase inhibitors were refined by SAR, were shown to be substrate competitive and were active in cellular assays. nCDase inhibitor compounds were modeled by in silico DOCK screening and by molecular simulation. Modeling data supports zinc interaction and a similar compound binding pose with ceramide. nCDase inhibitors were identified with notably improved activity and solubility in comparison with the reference lipid-mimetic C6-urea ceramide.

A molecular evolution algorithm for ligand design in DOCK.

As a complement to virtual screening, de novo design of small molecules is an alternative approach for identifying potential drug candidates. Here, we present a new 3D genetic algorithm to evolve molecules through breeding, mutation, fitness pressure, and selection. The method, termed DOCK_GA, builds upon and leverages powerful sampling, scoring, and searching routines previously implemented into DOCK6. Three primary experiments were used during development: Single-molecule evolution evaluated three selection methods (elitism, tournament, and roulette), in four clinically relevant systems, in terms of mutation type and crossover success, chemical properties, ensemble diversity, and fitness convergence, among others. Large scale benchmarking assessed performance across 651 different protein-ligand systems. Ensemble-based evolution demonstrated using multiple inhibitors simultaneously to seed growth in a SARS-CoV-2 target. Key takeaways include: (1) The algorithm is robust as demonstrated by the successful evolution of molecules across a large diverse dataset. (2) Users have flexibility with regards to parent input, selection method, fitness function, and molecular descriptors. (3) The program is straightforward to run and only requires a single executable and input file at run-time. (4) The elitism selection method yields more tightly clustered molecules in terms of 2D/3D similarity, with more favorable fitness, followed by tournament and roulette.

Free Energy Landscapes from SARS-CoV-2 Spike Glycoprotein Simulations Suggest that RBD Opening Can Be Modulated via Interactions in an Allosteric Pocket.

The SARS-CoV-2 coronavirus is an enveloped, positive-sense single-stranded RNA virus that is responsible for the COVID-19 pandemic. The spike is a class I viral fusion glycoprotein that extends from the viral surface and is responsible for viral entry into the host cell and is the primary target of neutralizing antibodies. The receptor binding domain (RBD) of the spike samples multiple conformations in a compromise between evading immune recognition and searching for the host-cell surface receptor. Using atomistic simulations of the glycosylated wild-type spike in the closed and 1-up RBD conformations, we map the free energy landscape for RBD opening and identify interactions in an allosteric pocket that influence RBD dynamics. The results provide an explanation for experimental observation of increased antibody binding for a clinical variant with a substitution in this pocket. Our results also suggest the possibility of allosteric targeting of the RBD equilibrium to favor open states via binding of small molecules to the hinge pocket. In addition to potential value as experimental probes to quantify RBD conformational heterogeneity, small molecules that modulate the RBD equilibrium could help explore the relationship between RBD opening and S1 shedding.

Identification of Zika Virus Inhibitors Using Homology Modeling and Similarity-Based Screening to Target Glycoprotein E.

The World Health Organization has designated Zika virus (ZIKV) as a dangerous, mosquito-borne pathogen that can cause severe developmental defects. The primary goal of this work was identification of small molecules as potential ZIKV inhibitors that target the viral envelope glycoprotein (ZIKV E) involved in membrane fusion and viral entry. A homology model of ZIKV E containing the small molecule ß-octyl glucoside (BOG) was constructed, on the basis of an analogous X-ray structure from dengue virus, and >4 million commercially available compounds were computationally screened using the program DOCK6. A key feature of the screen involved the use of similarity-based scoring to identify inhibitor candidates that make similar interaction energy patterns (molecular footprints) as the BOG reference. Fifty-three prioritized compounds underwent experimental testing using cytotoxicity, cell viability, and tissue culture infectious dose 50% (TCID50) assays. Encouragingly, relative to a known control (NITD008), six compounds were active in both the cell viability assay and the TCID50 infectivity assay, and they showed activity in a third caspase activity assay. In particular, compounds 8 and 15 (tested at 25 µM) and compound 43 (tested at 10 µM) appeared to provide significant protection to infected cells, indicative of anti-ZIKV activity. Overall, the study highlights how similarity-based scoring can be leveraged to computationally identify potential ZIKV E inhibitors that mimic a known reference (in this case BOG), and the experimentally verified hits provide a strong starting point for further refinement and optimization efforts.

Identification of Ebola Virus Inhibitors Targeting GP2 Using Principles of Molecular Mimicry.

A key step in the Ebola virus (EBOV) replication cycle involves conformational changes in viral glycoprotein 2 (GP2) which facilitate host-viral membrane fusion and subsequent release of the viral genome. Ebola GP2 plays a critical role in virus entry and has similarities in mechanism and structure to the HIV gp41 protein for which inhibitors have been successfully developed. In this work, a putative binding pocket for the C-terminal heptad repeat in the N-terminal heptad repeat trimer was targeted for identification of small molecules that arrest EBOV-host membrane fusion. Two computational structure-based virtual screens of ∼1.7 M compounds were performed (DOCK program) against a GP2 five-helix bundle, resulting in 165 commercially available compounds purchased for experimental testing. Based on assessment of inhibitory activity, cytotoxicity, and target specificity, four promising candidates emerged with 50% inhibitory concentration values in the 3 to 26 µM range. Molecular dynamics simulations of the two most potent candidates in their DOCK-predicted binding poses indicate that the majority of favorable interactions involve seven highly conserved residues that can be used to guide further inhibitor development and refinement targeting EBOV.IMPORTANCE The most recent Ebola virus disease outbreak, from 2014 to 2016, resulted in approximately 28,000 individuals becoming infected, which led to over 12,000 causalities worldwide. The particularly high pathogenicity of the virus makes paramount the identification and development of promising lead compounds to serve as inhibitors of Ebola infection. To limit viral load, the virus-host membrane fusion event can be targeted through the inhibition of the class I fusion glycoprotein of Ebolavirus In the current work, several promising small-molecule inhibitors that target the glycoprotein GP2 were identified through systematic application of structure-based computational and experimental drug design procedures.

Identification of Fatty Acid Binding Protein 5 Inhibitors Through Similarity-Based Screening.

Fatty acid binding protein 5 (FABP5) is a promising target for development of inhibitors to help control pain and inflammation. In this work, computer-based docking (DOCK6 program) was employed to screen ∼2 M commercially available compounds to FABP5 based on an X-ray structure complexed with the small molecule inhibitor SBFI-26 previously identified by our group (also through virtual screening). The goal was discovery of additional chemotypes. The screen resulted in the purchase of 78 candidates, which led to the identification of a new inhibitor scaffold (STK-0) with micromolar affinity and apparent selectivity for FABP5 over FABP3. A second similarity-based screen resulted in three additional hits (STK-15, STK-21, STK-22) from which preliminary SAR could be derived. Notably, STK-15 showed comparable activity to the SBFI-26 reference under the same assay conditions (1.40 vs 0.86 µM). Additional molecular dynamics simulations, free energy calculations, and structural analysis (starting from DOCK-generated poses) revealed that R enantiomers (dihydropyrrole scaffold) of STK-15 and STK-22 have a more optimal composition of functional groups to facilitate additional H-bonds with Arg109 of FABP5. This observation suggests enantiomerically pure compounds could show enhanced activity. Overall, our study highlights the utility of using similarity-based screening methods to discover new inhibitor chemotypes, and the identified FABP5 hits provide a strong starting point for future efforts geared to improve activity.

Identification of a Water-Coordinating HER2 Inhibitor by Virtual Screening Using Similarity-Based Scoring.

Human epidermal growth factor receptor 2 (HER2) is a validated breast cancer drug target for small molecule inhibitors that target the ATP-binding pocket of the kinase domain. In this work, a large-scale virtual screen was performed to a novel homology model of HER2, in a hypothesized "fully active" state, that considered water-mediated interactions during the prioritization of compounds for experimental testing. This screen led to the identification of a new inhibitor with micro molar affinity and potency ( Kd = 7.0 µM, IC50 = 4.6 µM). Accompanying molecular dynamics simulations showed that inhibitor binding likely involves water coordination through an important water-mediated network previously identified in our laboratory. The predicted binding geometry also showed a remarkable overlap with the crystallographic poses for two previously reported inhibitors of the related Chk1 kinase. Concurrent with the HER2 studies, we developed formalized computational protocols that leverage solvated footprints (per-residue interaction maps that include bridging waters) to identify ligands that can "coordinate" or "displace" key binding site waters. Proof-of-concept screens targeting HIVPR and PARP1 demonstrate that molecules with high footprint overlap can be effectively identified in terms of their coordination or displacement patterns relative to a known reference. Overall, the procedures developed as a result of this study should be useful for researchers targeting HER2 and, more generally, for any protein in which the identification of compounds that exploit binding site waters is desirable.

The Antinociceptive Agent SBFI-26 Binds to Anandamide Transporters FABP5 and FABP7 at Two Different Sites.

Human FABP5 and FABP7 are intracellular endocannabinoid transporters. SBFI-26 is an α-truxillic acid 1-naphthyl monoester that competitively inhibits the activities of FABP5 and FABP7 and produces antinociceptive and anti-inflammatory effects in mice. The synthesis of SBFI-26 yields several stereoisomers, and it is not known how the inhibitor binds the transporters. Here we report co-crystal structures of SBFI-26 in complex with human FABP5 and FABP7 at 2.2 and 1.9 Å resolution, respectively. We found that only (S)-SBFI-26 was present in the crystal structures. The inhibitor largely mimics the fatty acid binding pattern, but it also has several unique interactions. Notably, the FABP7 complex corroborates key aspects of the ligand binding pose at the canonical site previously predicted by virtual screening. In FABP5, SBFI-26 was unexpectedly found to bind at the substrate entry portal region in addition to binding at the canonical ligand-binding pocket. Our structural and binding energy analyses indicate that both R and S forms appear to bind the transporter equally well. We suggest that the S enantiomer observed in the crystal structures may be a result of the crystallization process selectively incorporating the (S)-SBFI-26-FABP complexes into the growing lattice, or that the S enantiomer may bind to the portal site more rapidly than to the canonical site, leading to an increased local concentration of the S enantiomer for binding to the canonical site. Our work reveals two binding poses of SBFI-26 in its target transporters. This knowledge will guide the development of more potent FABP inhibitors based upon the SBFI-26 scaffold.