Discovery and Mode-of-Action Characterization of a New Class of Acetolactate Synthase-Inhibiting Herbicides.

Herbicides are effective tools to manage weeds and enable food production and sustainable agriculture. Corteva Agriscience R&D has recently discovered new diphenyl-ether compounds displaying excellent postemergent efficacy on important weed species along with corn safety. Here, we describe the chemistry, biology, biochemistry, and computational modeling research that led to the discovery and elucidation of the primary mode of action for these compounds. The target protein was found to be acetolactate synthase (ALS), a key enzyme in the biosynthesis of branched chain amino acids (valine, leucine, and isoleucine). While weed resistance evolution to ALS herbicides is widespread, the molecular interaction of the diphenyl-ether compounds at the active site of the ALS enzyme differs significantly from that of some commercial ALS inhibitors. The unique biochemical profile of these molecules along with their excellent herbicidal activity and corn selectivity make them a noteworthy development in the pursuit of novel, safe, and sustainable weed control solutions.

A New Class of Diaryl Ether Herbicides: Structure-Activity Relationship Studies Enabled by a Rapid Scaffold Hopping Approach.

We report on the development of a novel class of diaryl ether herbicides. After the discovery of a phenoxybenzoic acid with modest herbicidal activity, optimization led to several molecules with improved control of broadleaf and grass weeds. To facilitate this process, we first employed a three-step combinatorial approach, then pivoted to a one-step Ullmann-type coupling that provided faster access to new analogs. After determining that the primary target site of our benchmark diaryl ethers was acetolactate synthase (ALS), we further leveraged this copper-catalyzed methodology to conduct a scaffold hopping campaign in the hope of uncovering an additional mode of action with fewer documented cases of resistance. Our comprehensive and systematic investigation revealed that while the herbicidal activity of this area seems to be exclusively linked to ALS inhibition, our molecules represent a structurally distinct class of Group 2 herbicides. The structure-activity relationships that led us to this conclusion are described herein.

Peptide Inhibitors Targeting the Neisseria gonorrhoeae Pivotal Anaerobic Respiration Factor AniA.

Neisseria gonorrhoeae causes the sexually transmitted infection gonorrhea, which is highly prevalent worldwide and has a major impact on reproductive and neonatal health. The superbug status of N. gonorrhoeae necessitates the development of drugs with different mechanisms of action. Here, we focused on targeting the nitrite reductase AniA, which is a pivotal component of N. gonorrhoeae anaerobic respiration and biofilm formation. Our studies showed that gonococci expressing AniA containing the altered catalytic residues D137A and H280A failed to grow under anaerobic conditions, demonstrating that the nitrite reductase function is essential. To facilitate the pharmacological targeting of AniA, new crystal structures of AniA were refined to 1.90-Å and 2.35-Å resolutions, and a phage display approach with libraries expressing randomized linear dodecameric peptides or heptameric peptides flanked by a pair of cysteine residues was utilized. Biopanning experiments led to the identification of 29 unique peptides, with 1 of them, C7-3, being identified multiple times. Evaluation of their ability to interact with AniA using enzyme-linked immunosorbent assay and computational docking studies revealed that C7-3 was the most promising inhibitor, binding near the type 2 copper site of the enzyme, which is responsible for interaction with nitrite. Subsequent enzymatic assays and biolayer interferometry with a synthetic C7-3 and its derivatives, C7-3m1 and C7-3m2, demonstrated potent inhibition of AniA. Finally, the MIC50 value of C7-3 and C7-3m2 against anaerobically grown N. gonorrhoeae was 0.6 mM. We present the first peptide inhibitors of AniA, an enzyme that should be further exploited for antigonococcal drug development.

Application of the 4D fingerprint method with a robust scoring function for scaffold-hopping and drug repurposing strategies.

Two factors contribute to the inefficiency associated with screening pharmaceutical library collections as a means of identifying new drugs: [1] the limited success of virtual screening (VS) methods in identifying new scaffolds; [2] the limited accuracy of computational methods in predicting off-target effects. We recently introduced a 3D shape-based similarity algorithm of the SABRE program, which encodes a consensus molecular shape pattern of a set of active ligands into a 4D fingerprint descriptor. Here, we report a mathematical model for shape similarity comparisons and ligand database filtering using this 4D fingerprint method and benchmarked the scoring function HWK (Hamza-Wei-Korotkov), using the 81 targets of the DEKOIS database. Subsequently, we applied our combined 4D fingerprint and HWK scoring function VS approach in scaffold-hopping and drug repurposing using the National Cancer Institute (NCI) and Food and Drug Administration (FDA) databases, and we identified new inhibitors with different scaffolds of MycP1 protease from the mycobacterial ESX-1 secretion system. Experimental evaluation of nine compounds from the NCI database and three from the FDA database displayed IC50 values ranging from 70 to 100 µM against MycP1 and possessed high structural diversity, which provides departure points for further structure-activity relationship (SAR) optimization. In addition, this study demonstrates that the combination of our 4D fingerprint algorithm and the HWK scoring function may provide a means for identifying repurposed drugs for the treatment of infectious diseases and may be used in the drug-target profile strategy.

Arylquins target vimentin to trigger Par-4 secretion for tumor cell apoptosis.

The tumor suppressor protein prostate apoptosis response-4 (Par-4), which is secreted by normal cells, selectively induces apoptosis in cancer cells. We identified a 3-arylquinoline derivative, designated Arylquin 1, as a potent Par-4 secretagogue in cell cultures and mice. Mechanistically, Arylquin 1 binds vimentin, displaces Par-4 from vimentin for secretion and triggers the efficient paracrine apoptosis of diverse cancer cells. Thus, targeting vimentin with Par-4 secretagogues efficiently induces paracrine apoptosis of tumor cells.

Novel mycosin protease MycP1 inhibitors identified by virtual screening and 4D fingerprints.

The rise of drug-resistant Mycobacterium tuberculosis lends urgency to the need for new drugs for the treatment of tuberculosis (TB). The identification of a serine protease, mycosin protease-1 (MycP1), as the crucial agent in hydrolyzing the virulence factor, ESX-secretion-associated protein B (EspB), potentially opens the door to new tuberculosis treatment options. Using the crystal structure of mycobacterial MycP1 in the apo form, we performed an iterative ligand- and structure-based virtual screening (VS) strategy to identify novel, nonpeptide, small-molecule inhibitors against MycP1 protease. Screening of ∼485,000 ligands from databases at the Genomics Research Institute (GRI) at the University of Cincinnati and the National Cancer Institute (NCI) using our VS approach, which integrated a pharmacophore model and consensus molecular shape patterns of active ligands (4D fingerprints), identified 81 putative inhibitors, and in vitro testing subsequently confirmed two of them as active inhibitors. Thereafter, the lead structures of each VS round were used to generate a new 4D fingerprint that enabled virtual rescreening of the chemical libraries. Finally, the iterative process identified a number of diverse scaffolds as lead compounds that were tested and found to have micromolar IC50 values against the MycP1 target. This study validated the efficiency of the SABRE 4D fingerprints as a means of identifying novel lead compounds in each screening round of the databases. Together, these results underscored the value of using a combination of in silico iterative ligand- and structure-based virtual screening of chemical libraries with experimental validation for the identification of promising structural scaffolds, such as the MycP1 inhibitors.

SABRE: ligand/structure-based virtual screening approach using consensus molecular-shape pattern recognition.

We present an efficient and rational ligand/structure shape-based virtual screening approach combining our previous ligand shape-based similarity SABRE (shape-approach-based routines enhanced) and the 3D shape of the receptor binding site. Our approach exploits the pharmacological preferences of a number of known active ligands to take advantage of the structural diversities and chemical similarities, using a linear combination of weighted molecular shape density. Furthermore, the algorithm generates a consensus molecular-shape pattern recognition that is used to filter and place the candidate structure into the binding pocket. The descriptor pool used to construct the consensus molecular-shape pattern consists of four dimensional (4D) fingerprints generated from the distribution of conformer states available to a molecule and the 3D shapes of a set of active ligands computed using SABRE software. The virtual screening efficiency of SABRE was validated using the Database of Useful Decoys (DUD) and the filtered version (WOMBAT) of 10 DUD targets. The ligand/structure shape-based similarity SABRE algorithm outperforms several other widely used virtual screening methods which uses the data fusion of multiscreening tools (2D and 3D fingerprints) and demonstrates a superior early retrieval rate of active compounds (EF(0.1%) = 69.0% and EF(1%) = 98.7%) from a large size of ligand database (∼95,000 structures). Therefore, our developed similarity approach can be of particular use for identifying active compounds that are similar to reference molecules and predicting activity against other targets (chemogenomics). An academic license of the SABRE program is available on request.

Structural insight into a novel human CCR5-V130I variant associated with resistance to HIV-1 infection.

We report the identification of a novel CC chemokine receptor 5 (CCR5) variant that seems associated with resistance to HIV-1 infection. The V130I mutation of the CCR5 receptor is located in the intracellular loop ICL2 known as DRY box and described in the literature as a nonsynonymous mutation present in nonhuman primates group. Extensive molecular modeling and dynamics simulations were performed to elucidate the mechanism by which the V130I mutation may induce conformational change of the CCR5 folding protein and prevent the interaction with the ß-arrestin protein. Our study provides new mechanistic insight into how a specific mutation in the regulatory domain of CCR5 might alter the structural folding of the DRY box and the possible ICL2 loop binding with the ß-arrestin protein, as described in our previous computational study. The results from our large-scale simulations complement recent experimental results and clinical features and offer useful insights into the mechanism behind CCR5 protein folding and signal transduction. In order for HIV, the entry of the virus to the cells must fuse with the CCR5 receptor that sits on the surface of T-helper immune cells. The described V130I mutation in the gene encoding the CCR5 protein may results in a defective CCR5-Arrestin binding complex that blocks entry of the virus.

Microscopic Modes and Free Energies for Topoisomerase I-DNA Covalent Complex Binding with Non-campothecin Inhibitors by Molecular Docking and Dynamics Simulations.

Topoisomerase I (Topo1) has been identified as an attractive target for anticancer drug development due to its central role in facilitating the nuclear process of the DNA. It is essential for rational design of novel Topo1 inhibitors to reliably predict the binding structures of the Topo1 inhibitors interacting with the Topo1-DNA complex. The detailed binding structures and binding free energies for the Topo1-DNA complex interacting with typical non-camptothecin (CPT) Topo1 inhibitors have been examined by performing molecular docking, molecular dynamic (MD) simulations, and binding free energy calculations. The computational results provide valuable insights into the binding modes of the inhibitors binding with the Topo1-DNA complex and the key factors affecting the binding affinity. It has been demonstrated that the - stacking interaction with the DNA base pairs and the hydrogen bonding with Topo1 have the pivotal contributions to the binding structures and binding free energies, although the van der Waals and electrostatic interactions also significantly contribute to the stabilization of the binding structures. The calculated binding free energies are in good agreement with the available experiment activity data. The detailed binding modes and the crucial factors affecting the binding free energies obtained from the present computational studies may provide valuable insights for future rational design of novel, more potent Topo1 inhibitors.

Fluorinated N,N-dialkylaminostilbenes repress colon cancer by targeting methionine S-adenosyltransferase 2A.

Methionine S-adenosyltransferase 2A (MAT2A) is the catalytic subunit for synthesis of S-adenosylmethionine (SAM), the principal methyl donor in many biological processes. MAT2A is up-regulated in many cancers, including liver cancer and colorectal cancer (CRC) and is a potentially important drug target. We developed a family of fluorinated N,N-dialkylaminostilbene agents, called FIDAS agents, that inhibit the proliferation of CRC cells in vitro and in vivo. Using a biotinylated FIDAS analogue, we identified the catalytic subunit of MAT2A as the direct and exclusive binding target of these FIDAS agents. MAT2B, an associated regulatory subunit of MAT2A, binds indirectly to FIDAS agents through its association with MAT2A. FIDAS agents inhibited MAT2A activity in SAM synthesis, and depletion of MAT2A by shRNAs inhibited CRC cell growth. A novel FIDAS agent delivered orally repressed CRC xenografts in athymic nude mice. These findings suggest that FIDAS analogues targeting MAT2A represent a family of novel and potentially useful agents for cancer treatment.

Assessing the regioselectivity of OleD-catalyzed glycosylation with a diverse set of acceptors.

To explore the acceptor regioselectivity of OleD-catalyzed glucosylation, the products of OleD-catalyzed reactions with six structurally diverse acceptors flavones- (daidzein), isoflavones (flavopiridol), stilbenes (resveratrol), indole alkaloids (10-hydroxycamptothecin), and steroids (2-methoxyestradiol)-were determined. This study highlights the first synthesis of flavopiridol and 2-methoxyestradiol glucosides and confirms the ability of OleD to glucosylate both aromatic and aliphatic nucleophiles. In all cases, molecular dynamics simulations were consistent with the determined product distribution and suggest the potential to develop a virtual screening model to identify additional OleD substrates.

Protein flexibility and conformational states of Leishmania antigen eIF-4A: identification of a novel plausible protein adjuvant using comparative genomics and molecular modeling.

Recent homology modeling studies have identified specific residues (epitope) of the Leishmania RNA helicase protein (LmeIF) that stimulates production of IL-12 cytokine. However, question remains concerning how LmeIF's N-terminal moiety initiates adjuvant effects. Extensive molecular modeling combining the normal mode analysis (NMA) and molecular dynamics simulations, in the present study, has demonstrated that the LmeIF structure may exist in two different forms corresponding to the extended and collapsed (closed) states of the entire structure. The computational results showed that the two domains of the LmeIF structure tend to undergo large fluctuations in a concerted fashion and have strong effect on the solvent accessible surface of the epitope situated on the N-terminal structure. The conformational freedom of the C-terminal domains may explain why the entire LmeIF protein is not as active as the N-terminal moiety. Thereafter, a comparative genome analysis with subsequent homology modeling and molecular electrostatic potential (MEP) techniques allowed us to predict a novel and plausible RNA helicase (LI-helicase) from the Listeria source with adjuvant property as observed for the Leishmania eIF-4A protein. The structural folding and MEP maps revealed similar topologies of the epitope of both LmeIF and LI-helicase proteins and striking identity in the local disposition of the charged groups. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:7.

A novel and efficient ligand-based virtual screening approach using the HWZ scoring function and an enhanced shape-density model.

In this work, we extend our previous ligand shape-based virtual screening approach by using the scoring function Hamza-Wei-Zhan (HWZ) score and an enhanced molecular shape-density model for the ligands. The performance of the method has been tested against the 40 targets in the Database of Useful Decoys and compared with the performance of our previous HWZ score method. The virtual screening results using the novel ligand shape-based approach demonstrated a favorable improvement (area under the receiver operator characteristics curve AUC = .89 ± .02) and effectiveness (hit rate HR(1%) = 53.0% ± 6.3 and HR(10%) = 71.1% ± 4.9). The comparison of the overall performance of our ligand shape-based method with the highest ligand shape-based virtual screening approach using the data fusion of multi queries showed that our strategy takes into account deeper the chemical information of the set of active ligands. Furthermore, the results indicated that our method are suitable for virtual screening and yields superior prediction accuracy than the other study derived from the data fusion using five queries. Therefore, our novel ligand shape-based screening method constitutes a robust and efficient approach to the 3D similarity screening of small compounds and open the door to a whole new approach to drug design by implementing the method in the structure-based virtual screening.

Probing the regiospecificity of enzyme-catalyzed steroid glycosylation.

The potential of a uniquely permissive engineered glycosyltransferase (OleD ASP) as a catalyst for steroid glycosylation is highlighted. The ability of OleD ASP to glucosylate a range of cardenolides and bufadienolides was assessed using a rapid LC-UV/MS-SPE-NMR analytical platform. While a bias toward OleD-catalyzed C3 monoglucosylation was observed, subtle alterations of the steroidal architecture, in some cases, invoked diglucosylation or, in one case (digoxigenin), C12 glucosylation. This latter case represents the first, and highly efficient, synthesis of digoxigenin 12-O-ß-D-glucoside.

Ligand-based virtual screening approach using a new scoring function.

In this study, we aimed to develop a new ligand-based virtual screening approach using an effective shape-overlapping procedure and a more robust scoring function (denoted by the HWZ score for convenience). The HWZ score-based virtual screening approach was tested against the compounds for 40 protein targets available in the Database of Useful Decoys (DUD; dud.docking.org/jahn/ ), and the virtual screening performance was evaluated in terms of the area under the receiver operator characteristic (ROC) curve (AUC), enrichment factor (EF), and hit rate (HR), demonstrating an improved overall performance compared to other popularly used approaches examined. In particular, the HWZ score-based virtual screening led to an average AUC value of 0.84 ± 0.02 (95% confidence interval) for the 40 targets. The average HR values at the top 1% and 10% of the active compounds for the 40 targets were 46.3% ± 6.7% and 59.2% ± 4.7%, respectively. In addition, the performance of the HWZ score-based virtual screening approach is less sensitive to the choice of the target.

Unveiling the unfolding pathway of F5F8D disorder-associated D81H/V100D mutant of MCFD2 via multiple molecular dynamics simulations.

Combined factor deficiency (F5F8D) is a rare autosomal recessive disorder caused by mutations in the LMAN1 or MCFD2 genes. It has been proposed that this pathogenic process occurs via a multi-step pathway involving metal loss, EF-hand-Ca21 dissociation and assembly of misfolded MCFD2-LMAN1 complex. Here, we have investigated the solution conformations of the MCFD2((D81H,V100D)) protein mutant through extensive molecular dynamics (MD) simulations. The V100D, one of the many MCFD2 mutations known to be associated to F5F8D, is difficult to be reconciled with the pathway model because it is located far from the metal sites and the MCFD2/LMAN1 interface. Consequently, an inspection of all the steps involved in D81H/V100D MCFD2 misfolding is expected to provide hints in the understanding of the molecular basis of the disease. A comparison with parallel studies carried out for the Wild-Type (WT) MCFD2 pointed out that the mutation decreases the affinity of the protein for the Ca21 ion. Multiple explicit solvents MD simulations (50 ns) performed on the two proteins revealed that in the WT protein, stable H-bond network and compact hydrophobic core region are created thus confirming a pivotal role of this region in driving the biophysical properties of the entire protein. In fact it is shown that the V100D mutation, although located far away the EF-hand domain, may induce subtle modification in the structural core of MCFD2 leading to the loosening of metal binding and to the formation of metastable intermediate states along the unfolding pathway. The native-like hydrophobic cluster formed near the V100 residue in the wild-type protein is disrupted by the negatively charged Asparagine residue. Furthermore, the presence of the D81H mutation in the EF-1 hand domain may also increase the protein unfolding rate and consequently prevent the formation of the MCFD2-LMAN1 complex. The detailed structural insights obtained from our large-scale simulations complement the clinical features and offer useful insights into the mechanism behind MCFD2 protein misfolding.

New inhibitor of 3-phosphoinositide dependent protein kinase-1 identified from virtual screening.

3-Phosphoinositide-dependent protein kinase-1 (PDK1) has been recognized as a promising anticancer target. Thus, it is interesting to identify new inhibitors of PDK1 for anticancer drug discovery. Through a combined use of virtual screening and wet experimental activity assays, we have identified a new PDK1 inhibitor with IC(50)=~200 nM. The anticancer activities of this compound have been confirmed by the anticancer activity assays using 60 cancer cell lines. The obtained new PDK1 inhibitor and its PDK1-inhibitor binding mode should be valuable in future de novo design of novel, more potent and selective PDK1 inhibitors for future development of anticancer therapeutics.

Corneal antifibrotic switch identified in genetic and pharmacological deficiency of vimentin.

The type III intermediate filaments (IFs) are essential cytoskeletal elements of mechanosignal transduction and serve critical roles in tissue repair. Mice genetically deficient for the IF protein vimentin (Vim(-/-)) have impaired wound healing from deficits in myofibroblast development. We report a surprising finding made in Vim(-/-) mice that corneas are protected from fibrosis and instead promote regenerative healing after traumatic alkali injury. This reparative phenotype in Vim(-/-) corneas is strikingly recapitulated by the pharmacological agent withaferin A (WFA), a small molecule that binds to vimentin and down-regulates its injury-induced expression. Attenuation of corneal fibrosis by WFA is mediated by down-regulation of ubiquitin-conjugating E3 ligase Skp2 and up-regulation of cyclin-dependent kinase inhibitors p27(Kip1) and p21(Cip1). In cell culture models, WFA exerts G(2)/M cell cycle arrest in a p27(Kip1)- and Skp2-dependent manner. Finally, by developing a highly sensitive imaging method to measure corneal opacity, we identify a novel role for desmin overexpression in corneal haze. We demonstrate that desmin down-regulation by WFA via targeting the conserved WFA-ligand binding site shared among type III IFs promotes further improvement of corneal transparency without affecting cyclin-dependent kinase inhibitor levels in Vim(-/-) mice. This dissociates a direct role for desmin in corneal cell proliferation. Taken together, our findings illuminate a previously unappreciated pathogenic role for type III IF overexpression in corneal fibrotic conditions and also validate WFA as a powerful drug lead toward anti-fibrosis therapeutic development.

Insight into the binding of the wild type and mutated alginate lyase (AlyVI) with its substrate: a computational and experimental study.

The homology model of the wild type alginate lyase (AlyVI) marine bacterium Vibrio sp. protein, was built using the crystal structure of the Family 7 alginate lyase from Sphingomonas sp. A1. To rationalize the observed structure-affinity relationships of aliginate lyase alyVI with its (GGG) substrate, molecular docking, MD imulations and binding free energy calculations followed by site-directed mutagenesis and alyVI activity assays were carried out. Per-residue decomposition of the (GGG) binding energy revealed that the most important contributions were from polar and charged residues, such as Asn138, Arg143, Asn217, and Lys308, while van der Waals interactions were responsible for binding with the catalytic His200 and Tyr312 residues. The mutants H200A, K308A, Y312A, Y312F, and W165A were found to be inactive or almost inactive. However, the catalytic efficiency (k(cat)/K(m)) of the double mutant L224V/D226G increased by two-fold compared to the wild type enzyme. This first structural model with its substrate binding mode and the agreement with experimental results provide a suitable base for the future rational design of new mutated alyVI structures with improved catalytic activity.

Novel human mPGES-1 inhibitors identified through structure-based virtual screening.

Microsomal prostaglandin E synthase-1 (mPGES-1) is an inducible prostaglandin E synthase after exposure to pro-inflammatory stimuli and, therefore, represents a novel target for therapeutic treatment of acute and chronic inflammatory disorders. It is essential to identify mPGES-1 inhibitors with novel scaffolds as new leads or hits for the purpose of drug design and discovery that aim to develop the next-generation anti-inflammatory drugs. Herein we report novel mPGES-1 inhibitors identified through a combination of large-scale structure-based virtual screening, flexible docking, molecular dynamics simulations, binding free energy calculations, and in vitro assays on the actual inhibitory activity of the computationally selected compounds. The computational studies are based on our recently developed three-dimensional (3D) structural model of mPGES-1 in its open state. The combined computational and experimental studies have led to identification of new mPGES-1 inhibitors with new scaffolds. In particular, (Z)-5-benzylidene-2-iminothiazolidin-4-one is a promising novel scaffold for the further rational design and discovery of new mPGES-1 inhibitors. To our best knowledge, this is the first time a 3D structural model of the open state mPGES-1 is used in structure-based virtual screening of a large library of available compounds for the mPGES-1 inhibitor identification. The positive experimental results suggest that our recently modeled trimeric structure of mPGES-1 in its open state is ready for the structure-based drug design and discovery.