Structure-based identification of inhibitors targeting obstruction of the HIVgp41 N-heptad repeat trimer

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.

Ebola and Marburg virus VP35 coiled-coil validated as antiviral target by tripartite split-GFP complementation

Ebola virus (EBOV) and Marburg virus (MARV) are highly pathogenic viruses in humans, against which approved antivirals are lacking. During EBOV and MARV infection, coiled-coil mediated oligomerization is essential for the virion protein 35 (VP35) polymerase co-factor function and type I interferon antagonism, making VP35 coiled-coil an elective drug target. We established a tripartite split-green fluorescent protein (GFP) fluorescence complementation (FC) system based on recombinant GFP-tagged EBOV and MARV VP35, which probes VP35 coiled-coil assembly by monitoring fluorescence on E. coli colonies, or in vitro in 96/384-multiwell. Oligomerization-defective VP35 mutants showed that correct coiled-coil knobs-into-holes pairing within VP35 oligomer is pre-requisite for GFP tags and GFP detector to reconstitute fluorescing full-length GFP. The method was validated by screening a small compound library, which identified Myricetin and 4,5,6,7-Tetrabromobenzotriazole as inhibitors of EBOV and MARV VP35 oligomerization-dependent FC with low-micromolar IC₅₀ values. These findings substantiate the VP35 coiled-coil value as antiviral target.

•

Ebola and Marburg virus VP35 oligomerize via trimeric and tetrameric coiled-coil

•

VP35 coiled-coil assembly triggers fluorescence of a tripartite split-GFP system

•

Mutations perturbing VP35 coiled-coil hamper split-GFP complementation

•

Myricetin and TBBT inhibit split-GFP complementation mediated by VP35 coiled-coil

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.

Suitable fusion of N-terminal heptad repeats to achieve covalently stabilized potent N-peptide inhibitors of HIV-1 infection

N-terminal heptad repeat (NHR)-derived peptide (N-peptide) fusion inhibitors, which are derived from human immunodeficiency virus (HIV) envelope glycoprotein 41 (gp41), are limited by aggregation and unstable trimer conformation. However, they could function as potent inhibitors of viral infection by forming a coiled-coil structure covalently stabilized by interchain disulfide bonds. We previously synthesized N-peptides with potent anti-HIV-1 activity and high stability by coiled-coil fusion and covalent stabilization. Here, we attempted to study the effects of NHRs of chimeric N-peptides by fusing de novo coiled-coil isopeptide bridge-tethered T21 peptides of different NHR lengths. Peptides (T21N23)₃ and (T21N36)₃ was a more potent HIV-1 fusion inhibitor than (T21N17)₃. The site of isopeptide bond formation was precisely controlled and had little influence on N-peptide properties. The N-peptide (T21N36)₃, which had a similar conformation as the NHR trimer and interacted well with the C34 peptide, may be useful for screening other C-peptides and small-molecule fusion inhibitors, and for studying the interactions between the NHR trimer and C-terminal heptad repeats.

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 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 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 (K_d = 7.0 μM, IC₅₀ = 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.