1
|
Hayakawa D, Watanabe Y, Gouda H. Molecular Interaction Fields Describing Halogen Bond Formable Areas on Protein Surfaces. J Chem Inf Model 2024; 64:6003-6013. [PMID: 39012240 PMCID: PMC11323840 DOI: 10.1021/acs.jcim.4c00896] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2024] [Revised: 06/22/2024] [Accepted: 07/08/2024] [Indexed: 07/17/2024]
Abstract
Molecular interaction fields (MIFs) are three-dimensional interaction maps that describe the intermolecular interactions expected to be formed around target molecules. In this paper, a method for the fast computation of MIFs using the approximation functions of quantum mechanics-level MIFs of small model molecules is proposed. MIF functions of N-methylacetamide with chlorobenzene, bromobenzene, and iodobenzene probes were precisely approximated and used to calculate the MIFs on protein surfaces. This method appropriately reproduced halogen-bond-formable areas around the ligand-binding sites of proteins, where halogen bond formation was suggested in a previous study.
Collapse
Affiliation(s)
- Daichi Hayakawa
- Division of Biophysical
Chemistry,
Department of Pharmaceutical Sciences, Graduate School of Pharmacy, Showa University, 1-5-8 Hatanodai, Shinagawa-ku, Tokyo 142-8555, Japan
| | - Yurie Watanabe
- Division of Biophysical
Chemistry,
Department of Pharmaceutical Sciences, Graduate School of Pharmacy, Showa University, 1-5-8 Hatanodai, Shinagawa-ku, Tokyo 142-8555, Japan
| | - Hiroaki Gouda
- Division of Biophysical
Chemistry,
Department of Pharmaceutical Sciences, Graduate School of Pharmacy, Showa University, 1-5-8 Hatanodai, Shinagawa-ku, Tokyo 142-8555, Japan
| |
Collapse
|
2
|
Khan MA, Al Mamun Khan MA, Mahfuz AMUB, Sanjana JM, Ahsan A, Gupta DR, Hoque MN, Islam T. Highly potent natural fungicides identified in silico against the cereal killer fungus Magnaporthe oryzae. Sci Rep 2022; 12:20232. [PMID: 36418863 PMCID: PMC9684433 DOI: 10.1038/s41598-022-22217-w] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2022] [Accepted: 10/11/2022] [Indexed: 11/25/2022] Open
Abstract
Magnaporthe oryzae is one of the most notorious fungal pathogens that causes blast disease in cereals, and results in enormous loss of grain production. Many chemical fungicides are being used to control the pathogen but none of them are fully effective in controlling blast disease. Therefore, there is a demand for the discovery of a new natural biofungicide to manage the blast disease efficiently. A large number of new natural products showed inhibitory activities against M. oryzae in vitro. To find out effective biofungicides, we performed in silico molecular docking analysis of some of the potent natural compounds targeting four enzymes namely, scytalone dehydratase, SDH1 (PDB ID:1STD), trihydroxynaphthalene reductase, 3HNR (PDB ID:1YBV), trehalose-6-phosphate synthase, Tps1 (PDB ID:6JBI) and isocitrate lyase, ICL1 (PDB ID:5E9G) of M. oryzae fungus that regulate melanin biosynthesis and/or appresorium formation. Thirty-nine natural compounds that were previously reported to inhibit the growth of M. oryzae were subjected to rigid and flexible molecular docking against aforementioned enzymes followed by molecular dynamic simulation. The results of virtual screening showed that out of 39, eight compounds showed good binding energy with any one of the target enzymes as compared to reference commercial fungicides, azoxystrobin and strobilurin. Among the compounds, camptothecin, GKK1032A2 and chaetoviridin-A bind with more than one target enzymes of M. oryzae. All of the compounds except tricyclazole showed good bioactivity score. Taken together, our results suggest that all of the eight compounds have the potential to develop new fungicides, and remarkably, camptothecin, GKK1032A2 and chaetoviridin-A could act as multi-site mode of action fungicides against the blast fungus M. oryzae.
Collapse
Affiliation(s)
- Md. Arif Khan
- grid.443057.10000 0004 4683 7084Department of Biotechnology and Genetic Engineering, University of Development Alternative, Dhaka, 1209 Bangladesh
| | - Md. Abdullah Al Mamun Khan
- grid.443019.b0000 0004 0479 1356Department of Biotechnology and Genetic Engineering, Mawlana Bhashani Science and Technology University, Tangail, 1902 Bangladesh
| | - A. M. U. B. Mahfuz
- grid.443057.10000 0004 4683 7084Department of Biotechnology and Genetic Engineering, University of Development Alternative, Dhaka, 1209 Bangladesh
| | - Jannatul Maowa Sanjana
- grid.443019.b0000 0004 0479 1356Department of Biotechnology and Genetic Engineering, Mawlana Bhashani Science and Technology University, Tangail, 1902 Bangladesh
| | - Asif Ahsan
- grid.411511.10000 0001 2179 3896Department of Biotechnology, Bangladesh Agricultural University, Mymensingh, 2202 Bangladesh
| | - Dipali Rani Gupta
- grid.443108.a0000 0000 8550 5526Institute of Biotechnology and Genetic Engineering (IBGE), Bangabandhu Sheikh Mujibur Rahman Agricultural University (BSMRAU), Gazipur, 1706 Bangladesh
| | - M. Nazmul Hoque
- grid.443108.a0000 0000 8550 5526Department of Gynecology, Obstetrics and Reproductive Health, BSMRAU, Gazipur, 1706 Bangladesh
| | - Tofazzal Islam
- grid.443108.a0000 0000 8550 5526Institute of Biotechnology and Genetic Engineering (IBGE), Bangabandhu Sheikh Mujibur Rahman Agricultural University (BSMRAU), Gazipur, 1706 Bangladesh
| |
Collapse
|
3
|
Kim KH. Outliers in SAR and QSAR: 3. Importance of considering the role of water molecules in protein-ligand interactions and quantitative structure-activity relationship studies. J Comput Aided Mol Des 2021; 35:371-396. [PMID: 33712973 DOI: 10.1007/s10822-021-00377-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2020] [Accepted: 02/05/2021] [Indexed: 11/30/2022]
Abstract
It is frequently mentioned that QSARs have not generally lived up to expectations, especially in cases where high predictability is expected yet failed to deliver satisfactory results. Even though outliers can provide an increased opportunity in drug discovery research, outliers in SAR and QSAR can contort predictions and affect the accuracy if proper attention is not given. The percentages of outliers in QSARs have not changed appreciably over the last decade. In our previous studies, we suggested two possible sources of outliers in SAR and QSAR. In this paper, we suggest an additional possible source of outliers in QSAR. We presented several literature examples that show one or more water molecules that play a critical role in protein-ligand binding interactions as observed in their crystal structures. These examples illustrate that failing to account for the effects of water molecules in protein-ligand interactions could mislead interpretation and possibly yield outliers in SAR and QSAR. Examples include cases where QSAR, considering the role of water molecules in protein-ligand crystal structures, provided deeper insight into the understanding and interpretation of the developed QSAR.
Collapse
|
4
|
Crystal structures of non-oxidative decarboxylases reveal a new mechanism of action with a catalytic dyad and structural twists. Sci Rep 2021; 11:3056. [PMID: 33542397 PMCID: PMC7862292 DOI: 10.1038/s41598-021-82660-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2020] [Accepted: 01/22/2021] [Indexed: 12/01/2022] Open
Abstract
Hydroxybenzoic acids, like gallic acid and protocatechuic acid, are highly abundant natural compounds. In biotechnology, they serve as critical precursors for various molecules in heterologous production pathways, but a major bottleneck is these acids’ non-oxidative decarboxylation to hydroxybenzenes. Optimizing this step by pathway and enzyme engineering is tedious, partly because of the complicating cofactor dependencies of the commonly used prFMN-dependent decarboxylases. Here, we report the crystal structures (1.5–1.9 Å) of two homologous fungal decarboxylases, AGDC1 from Arxula adenivorans, and PPP2 from Madurella mycetomatis. Remarkably, both decarboxylases are cofactor independent and are superior to prFMN-dependent decarboxylases when heterologously expressed in Saccharomyces cerevisiae. The organization of their active site, together with mutational studies, suggests a novel decarboxylation mechanism that combines acid–base catalysis and transition state stabilization. Both enzymes are trimers, with a central potassium binding site. In each monomer, potassium introduces a local twist in a β-sheet close to the active site, which primes the critical H86-D40 dyad for catalysis. A conserved pair of tryptophans, W35 and W61, acts like a clamp that destabilizes the substrate by twisting its carboxyl group relative to the phenol moiety. These findings reveal AGDC1 and PPP2 as founding members of a so far overlooked group of cofactor independent decarboxylases and suggest strategies to engineer their unique chemistry for a wide variety of biotechnological applications.
Collapse
|
5
|
Guo J, Collins S, Miller WT, Rizzo RC. Identification of a Water-Coordinating HER2 Inhibitor by Virtual Screening Using Similarity-Based Scoring. Biochemistry 2018; 57:4934-4951. [PMID: 29975516 PMCID: PMC6110523 DOI: 10.1021/acs.biochem.8b00524] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
![]()
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.
Collapse
|
6
|
Forli S. Charting a Path to Success in Virtual Screening. Molecules 2015; 20:18732-58. [PMID: 26501243 PMCID: PMC4630810 DOI: 10.3390/molecules201018732] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2015] [Revised: 10/07/2015] [Accepted: 10/12/2015] [Indexed: 12/27/2022] Open
Abstract
Docking is commonly applied to drug design efforts, especially high-throughput virtual screenings of small molecules, to identify new compounds that bind to a given target. Despite great advances and successful applications in recent years, a number of issues remain unsolved. Most of the challenges and problems faced when running docking experiments are independent of the specific software used, and can be ascribed to either improper input preparation or to the simplified approaches applied to achieve high-throughput speed. Being aware of approximations and limitations of such methods is essential to prevent errors, deal with misleading results, and increase the success rate of virtual screening campaigns. In this review, best practices and most common issues of docking and virtual screening will be discussed, covering the journey from the design of the virtual experiment to the hit identification.
Collapse
Affiliation(s)
- Stefano Forli
- Molecular Graphics Laboratory, Department of Integrative Structural and Computational Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA.
| |
Collapse
|
7
|
Bayden AS, Moustakas DT, Joseph-McCarthy D, Lamb ML. Evaluating Free Energies of Binding and Conservation of Crystallographic Waters Using SZMAP. J Chem Inf Model 2015; 55:1552-65. [DOI: 10.1021/ci500746d] [Citation(s) in RCA: 56] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Alexander S. Bayden
- Oncology and Infection Innovative Medicines Units, AstraZeneca R&D Boston, 35 Gatehouse Drive, Waltham, Massachusetts 02451, United States
| | - Demetri T. Moustakas
- Oncology and Infection Innovative Medicines Units, AstraZeneca R&D Boston, 35 Gatehouse Drive, Waltham, Massachusetts 02451, United States
| | - Diane Joseph-McCarthy
- Oncology and Infection Innovative Medicines Units, AstraZeneca R&D Boston, 35 Gatehouse Drive, Waltham, Massachusetts 02451, United States
| | - Michelle L. Lamb
- Oncology and Infection Innovative Medicines Units, AstraZeneca R&D Boston, 35 Gatehouse Drive, Waltham, Massachusetts 02451, United States
| |
Collapse
|
8
|
Positive evolutionary selection of an HD motif on Alzheimer precursor protein orthologues suggests a functional role. PLoS Comput Biol 2012; 8:e1002356. [PMID: 22319430 PMCID: PMC3271017 DOI: 10.1371/journal.pcbi.1002356] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2011] [Accepted: 12/07/2011] [Indexed: 12/31/2022] Open
Abstract
HD amino acid duplex has been found in the active center of many different enzymes. The dyad plays remarkably different roles in their catalytic processes that usually involve metal coordination. An HD motif is positioned directly on the amyloid beta fragment (Aβ) and on the carboxy-terminal region of the extracellular domain (CAED) of the human amyloid precursor protein (APP) and a taxonomically well defined group of APP orthologues (APPOs). In human Aβ HD is part of a presumed, RGD-like integrin-binding motif RHD; however, neither RHD nor RXD demonstrates reasonable conservation in APPOs. The sequences of CAEDs and the position of the HD are not particularly conserved either, yet we show with a novel statistical method using evolutionary modeling that the presence of HD on CAEDs cannot be the result of neutral evolutionary forces (p<0.0001). The motif is positively selected along the evolutionary process in the majority of APPOs, despite the fact that HD motif is underrepresented in the proteomes of all species of the animal kingdom. Position migration can be explained by high probability occurrence of multiple copies of HD on intermediate sequences, from which only one is kept by selective evolutionary forces, in a similar way as in the case of the “transcription binding site turnover.” CAED of all APP orthologues and homologues are predicted to bind metal ions including Amyloid-like protein 1 (APLP1) and Amyloid-like protein 2 (APLP2). Our results suggest that HDs on the CAEDs are most probably key components of metal-binding domains, which facilitate and/or regulate inter- or intra-molecular interactions in a metal ion-dependent or metal ion concentration-dependent manner. The involvement of naturally occurring mutations of HD (Tottori (D7N) and English (H6R) mutations) in early onset Alzheimer's disease gives additional support to our finding that HD has an evolutionary preserved function on APPOs. HD amino acid duplex can be found in the active center of different metallo-enzymes. An HD motif is positioned directly on the amyloid beta (Aβ) fragment and on the carboxy-terminal region of the extracellular domain of the human amyloid precursor protein (APP) and a taxonomically well defined group of APP orthologues (APPOs). The conservation of the HD dyad is not position specific and it cannot be seen in a multiple alignment. Yet we show with a novel statistical method using evolutionary modeling that HD motif is positively selected by evolution on APPOs, despite the fact that HD dyad is underrepresented in the proteomes of all species of the animal kingdom. CAED of all APP orthologues and homologues are predicted to bind metal ions including Amyloid-like protein 1 (APLP1) and Amyloid-like protein 2 (APLP2). Our results suggest that HDs on the APPOs are most probably key components of metal-binding domains, which facilitate and/or regulate inter- or intra-molecular interactions in a metal ion-dependent or metal ion concentration-dependent manner. The involvement of naturally occurring mutations of HD (Tottori (D7N) and English (H6R)) in early onset Alzheimer's disease gives additional support to our finding that HD has an evolutionary preserved function on APPOs.
Collapse
|
9
|
Cottrell SJ, Gillet VJ, Taylor R, Wilton DJ. Generation of multiple pharmacophore hypotheses using multiobjective optimisation techniques. J Comput Aided Mol Des 2005; 18:665-82. [PMID: 15865060 DOI: 10.1007/s10822-004-5523-7] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
Pharmacophore methods provide a way of establishing a structure activity relationship for a series of known active ligands. Often, there are several plausible hypotheses that could explain the same set of ligands and, in such cases, it is important that the chemist is presented with alternatives that can be tested with different synthetic compounds. Existing pharmacophore methods involve either generating an ensemble of conformers and considering each conformer of each ligand in turn or exploring conformational space on-the-fly. The ensemble methods tend to produce a large number of hypotheses and require considerable effort to analyse the results, whereas methods that vary conformation on-the-fly typically generate a single solution that represents one possible hypothesis, even though several might exist. We describe a new method for generating multiple pharmacophore hypotheses with full conformational flexibility being explored on-the-fly. The method is based on multiobjective evolutionary algorithm techniques and is designed to search for an ensemble of diverse yet plausible overlays which can then be presented to the chemist for further investigation.
Collapse
Affiliation(s)
- Simon J Cottrell
- Department of Information Studies, University of Sheffield, Regent Court, 211 Portobello Street, Sheffield S1 4DP, UK
| | | | | | | |
Collapse
|
10
|
Paulini R, Müller K, Diederich F. Orthogonale multipolare Wechselwirkungen in chemischen und biologischen Strukturen. Angew Chem Int Ed Engl 2005. [DOI: 10.1002/ange.200462213] [Citation(s) in RCA: 99] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
11
|
Paulini R, Müller K, Diederich F. Orthogonal Multipolar Interactions in Structural Chemistry and Biology. Angew Chem Int Ed Engl 2005; 44:1788-805. [PMID: 15706577 DOI: 10.1002/anie.200462213] [Citation(s) in RCA: 378] [Impact Index Per Article: 19.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
The past few decades of molecular recognition studies have greatly enhanced our knowledge on apolar, ion-dipole, and hydrogen-bonding interactions. However, much less attention has been given to the role that multipolar interactions, in particular those with orthogonal dipolar alignment, play in organizing a crystal lattice or stabilizing complexes involving biological receptors. By using results from database mining, this review attempts to give an overview of types and structural features of these previously rather overlooked interactions. A number of illustrative examples of these interactions found in X-ray crystal structures of small molecules and protein-ligand complexes demonstrate their propensity and thus potential importance for both, chemical and biological molecular recognition processes.
Collapse
Affiliation(s)
- Ralph Paulini
- Laboratorium für Organische Chemie, ETH-Hönggerberg, HCI, 8093 Zürich, Switzerland
| | | | | |
Collapse
|
12
|
Okimoto N, Nakamura T, Suenaga A, Futatsugi N, Hirano Y, Yamaguchi I, Ebisuzaki T. Cooperative Motions of Protein and Hydration Water Molecules: Molecular Dynamics Study of Scytalone Dehydratase. J Am Chem Soc 2004; 126:13132-9. [PMID: 15469312 DOI: 10.1021/ja048053u] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Two molecular dynamics (MD) simulations totaling 25 ns of simulation time of monomeric scytalone dehydratase (SD) were performed. The enzyme has a ligand-binding pocket containing a cone-shaped alpha+beta barrel, and the C-terminal region covers the binding pocket. Our simulations clarified the difference in protein dynamics and conformation between the liganded protein and the unliganded protein. The liganded protein held the ligand molecule tightly and the initial structure was maintained during the simulation. The unliganded protein, on the other hand, fluctuated dynamically and its structure changed largely from the initial structure. In the equilibrium state, the binding pocket was fully solvated by opening of the C-terminal region, and the protein dynamics was connected with hydration water molecules entry into and release from the binding pocket. In addition, the cooperative motions of the unliganded protein and the hydration water molecules produced the path through the protein interior for ligand binding.
Collapse
Affiliation(s)
- Noriaki Okimoto
- Contribution from the Computational Astrophysics Laboratory, Institute of Physical and Chemical Research (RIKEN), 2-1 Hirosawa, Wako, Saitama 351-0198, Japan.
| | | | | | | | | | | | | |
Collapse
|
13
|
Langfelder K, Streibel M, Jahn B, Haase G, Brakhage AA. Biosynthesis of fungal melanins and their importance for human pathogenic fungi. Fungal Genet Biol 2003; 38:143-58. [PMID: 12620252 DOI: 10.1016/s1087-1845(02)00526-1] [Citation(s) in RCA: 409] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
For more than 40 years fungi have been known to produce pigments known as melanins. Predominantly these have been dihydroxyphenylalanine (DOPA)-melanin and dihydroxynaphthalene (DHN)-melanin. The biochemical and genetical analysis of the biosynthesis pathways have led to the identification of the genes and corresponding enzymes of the pathways. Only recently have both these types of melanin been linked to virulence in some human pathogenic and phytopathogenic fungi. The absence of melanin in human pathogenic and phytopathogenic fungi often leads to a decrease in virulence. In phytopathogenic fungi such as Magnaporthe grisea and Colletotrichum lagenarium, besides other possible functions in pathogenicity, DHN-melanin plays an essential role in generating turgor for plant appressoria to penetrate plant leaves. While the function of melanin in human pathogenic fungi such as Cryptococcus neoformans, Wangiella dermatitidis, Sporothrix schenckii, and Aspergillus fumigatus is less well defined, its role in protecting fungal cells has clearly been shown. Specifically, the ability of both DOPA- and DHN-melanins to quench free radicals is thought to be an important factor in virulence. In addition, in several fungi the production of fungal virulence factors, such as melanin, has been linked to a cAMP-dependent signaling pathway. Many of the components involved in the signaling pathway have been identified.
Collapse
Affiliation(s)
- Kim Langfelder
- Institut für Mikrobiologie, Universität Hannover, Schneiderberg 50, 30167, Hannover, Germany
| | | | | | | | | |
Collapse
|
14
|
Basarab GS, Jordan DB, Gehret TC, Schwartz RS. Design of inhibitors of scytalone dehydratase: probing interactions with an asparagine carboxamide. Bioorg Med Chem 2002; 10:4143-54. [PMID: 12413868 DOI: 10.1016/s0968-0896(02)00272-9] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
Among the active-site residues of scytalone dehydratase, the side-chain carboxamide of asparagine 131 has the greatest potential for strong electrostatic interactions. Structure-based inhibitor design aimed at enhancing interactions with this residue led to the synthesis of a series of highly potent inhibitors that have a five- or six-membered ring containing a carbonyl functionality for hydrogen bonding. To achieve a good orientation for hydrogen bonding, the inhibitors incorporate a phenyl substituent that displaces a phenylalanine residue away from the five- or six-membered rings. Without the phenyl substituent, inhibitor binding potency is diminished by three orders of magnitude. Larger K(i) values of a site-directed mutant (Asn131Ala) of scytalone dehydratase in comparison to those of wild-type enzyme validate the design concept. The most potent inhibitor (K(i)=15 pM) contains a tetrahydrothiophenone that can form a single hydrogen bond with the asparagine carboxamide. Inhibitors with a butyrolactam that can form two hydrogen bonds with the asparagine carboxamide demonstrate excellent in vivo fungicidal activity.
Collapse
Affiliation(s)
- Gregory S Basarab
- DuPont Central Research & Development, Experimental Station, Wilmington, DE 19880, USA.
| | | | | | | |
Collapse
|
15
|
Jordan DB, Basarab GS, Liao DI, Johnson WM, Winzenberg KN, Winkler DA. Structure-based design of inhibitors of the rice blast fungal enzyme trihydroxynaphthalene reductase. J Mol Graph Model 2002; 19:434-47, 470-1. [PMID: 11552692 DOI: 10.1016/s1093-3263(00)00103-0] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Rice Blast Disease, caused by the fungus Pyricularia oryzae, is one of the most important diseases of rice. Several enzymes in the melanin biosynthetic pathway have proven to be valuable targets for development of rice blast fungicides. In particular, inhibitors of trihydroxynaphthalene reductase (3HNR), which catalyzes the conversion of trihydroxynaphthalene to vermelone, have yielded commercially useful rice fungicides. The X-ray structure of 3HNR has been published recently, presenting an opportunity to use this information in the de novo design of novel 3HNR inhibitors that may exhibit useful rice blast activity. We used the LeapFrog program to develop a docking model for interaction of ligands with the active site of THNR. The final model gave a good correlation between calculated binding energy and log Ki and was used to design novel ligands and score compounds for synthesis. Using this as a tool, we synthesized inhibitors in the nanomolar range and also developed several inhibitors that did not conform to the properties of the THNR active site. Leapfrog was able to locate a previously unrecognized binding pocket that could accommodate these otherwise anomalous regions of structure.
Collapse
Affiliation(s)
- D B Jordan
- E. I. DuPont de Nemours and Co., Experimental Station, Wilmington, Delaware, USA
| | | | | | | | | | | |
Collapse
|
16
|
Liao D, Basarab GS, Gatenby AA, Valent B, Jordan DB. Structures of trihydroxynaphthalene reductase-fungicide complexes: implications for structure-based design and catalysis. Structure 2001; 9:19-27. [PMID: 11342131 DOI: 10.1016/s0969-2126(00)00548-7] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
BACKGROUND Trihydroxynaphthalene reductase catalyzes two intermediate steps in the fungal melanin biosynthetic pathway. The enzyme, a typical short-chain dehydrogenase, is the biochemical target of three commercial fungicides. The fungicides bind preferentially to the NADPH form of the enzyme. RESULTS Three X-ray structures of the Magnaporthe grisea enzyme complexed with NADPH and two commercial and one experimental fungicide were determined at 1.7 A (pyroquilon), 2.0 A (2,3-dihydro-4-nitro-1H-inden-1-one, 1), and 2.1 A (phthalide) resolutions. The chemically distinct inhibitors occupy similar space within the enzyme's active site. The three inhibitors share hydrogen bonds with the side chain hydroxyls of Ser-164 and Tyr-178 via a carbonyl oxygen (pyroquilon and 1) or via a carbonyl oxygen and a ring oxygen (phthalide). Active site residues occupy similar positions among the three structures. A buried water molecule that is hydrogen bonded to the NZ nitrogen of Lys-182 in each of the three structures likely serves to stabilize the cationic form of the residue for participation in catalysis. CONCLUSIONS The pro S hydrogen of NADPH (which is transferred as a hydride to the enzyme's naphthol substrates) is directed toward the carbonyl carbon of the inhibitors that mimic an intermediate along the reaction coordinate. Modeling tetrahydroxynaphthalene and trihydroxynaphthalene in the active site shows steric and electrostatic repulsion between the extra hydroxyl oxygen of the former substrate and the sulfur atom of Met-283 (the C-terminal residue), which accounts, in part, for the 4-fold greater substrate specificity for trihydroxynaphthalene over tetrahydroxynaphthalene.
Collapse
Affiliation(s)
- D Liao
- DuPont Central Research and Development Experimental Station, Wilmington, DE 19880, USA.
| | | | | | | | | |
Collapse
|
17
|
Thompson JE, Fahnestock S, Farrall L, Liao DI, Valent B, Jordan DB. The second naphthol reductase of fungal melanin biosynthesis in Magnaporthe grisea: tetrahydroxynaphthalene reductase. J Biol Chem 2000; 275:34867-72. [PMID: 10956664 DOI: 10.1074/jbc.m006659200] [Citation(s) in RCA: 104] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Mutants of Magnaporthe grisea harboring a defective gene for 1,3, 8-trihydroxynaphthalene reductase retain the capability to produce scytalone, thus suggesting the existence of a second naphthol reductase that can catalyze the reduction of 1,3,6, 8-tetrahydroxynaphthalene to scytalone within the fungal melanin biosynthetic pathway. The second naphthol reductase gene was cloned from M. grisea by identification of cDNA fragments with weak homology to the cDNA of trihydroxynaphthalene reductase. The amino acid sequence for the second naphthol reductase is 46% identical to that of trihydroxynaphthalene reductase. The second naphthol reductase was produced in Esherichia coli and purified to homogeneity. Substrate competition experiments indicate that the second reductase prefers tetrahydroxynaphthalene over trihydroxynaphthalene by a factor of 310; trihydroxynaphthalene reductase prefers trihydroxynaphthalene over tetrahydroxynaphthalene by a factor of 4.2. On the basis of the 1300-fold difference in substrate specificities between the two reductases, the second reductase is designated tetrahydroxynaphthalene reductase. Tetrahydroxynaphthalene reductase has a 200-fold larger K(i) for the fungicide tricyclazole than that of trihydroxynaphthalene reductase, and this accounts for the latter enzyme being the primary physiological target of the fungicide. M. grisea mutants lacking activities for both trihydroxynaphthalene and tetrahydroxynaphthalene reductases do not produce scytalone, indicating that there are no other metabolic routes to scytalone.
Collapse
Affiliation(s)
- J E Thompson
- DuPont Agricultural Products, Stine-Haskell Research Center, Newark, Delaware 19714, USA
| | | | | | | | | | | |
Collapse
|
18
|
Abstract
[reaction--see text] The major conformation of scytalone has an envelope shape with C3 forming the flap and the C3 hydroxyl in the equatorial position as determined by quantum mechanical calculations and corroborated by NMR. The C2 axial pro-R is slower to exchange with solvent than the equatorial pro-S hydrogen. Modeling the transition state for enolate formation points to a deprotonation through the flipped envelope conformation in which the C3-hydroxyl and the C2 pro-S hydrogen are axial.
Collapse
Affiliation(s)
- G S Basarab
- E. I. DuPont de Nemours Central Research and Development, Wilmington, Delaware 19880-0328, USA.
| | | | | |
Collapse
|
19
|
Jennings LD, Rayner DR, Jordan DB, Okonya JF, Basarab GS, Amorose DK, Anaclerio BM, Lee JK, Schwartz RS, Whitmore KA. Cyclobutane carboxamide inhibitors of fungal melanin: biosynthesis and their evaluation as fungicides. Bioorg Med Chem 2000; 8:897-907. [PMID: 10882002 DOI: 10.1016/s0968-0896(00)00034-1] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
A new fungicide lead has been identified by enzyme screening of a focused combinatorial library. The lead compound 4, a potent inhibitor of scytalone dehydratase (SD), exhibits fungicidal activity upon foliar application but does not show systemic activity. The X-ray crystal structure of the enzyme-inhibitor complex and an appreciation for the relationship between physical properties and systemic activity enabled us to rapidly improve upon this initial lead. The geminal halogen-methyl group combination was found to be optimal for interaction with the bounding serine and asparagine side-chain residues. Replacement of CF3 with methyl was a key discovery, giving inhibitors with slightly diminished enzyme inhibition potency while significantly increasing systemic activity. Amides prepared from amines with 2,4-dichloro substitution on the phenyl ring gave the most potent enzyme inhibitors. Two compounds from this series showed systemic activity comparable to the commercial standard and were selected for outdoor testing in flooded plots which simulate rice paddies.
Collapse
Affiliation(s)
- L D Jennings
- DuPont Agricultural Products, Stine Haskell Research Center, Newark, DE 19714, USA
| | | | | | | | | | | | | | | | | | | |
Collapse
|
20
|
Firestine SM, Salinas F, Nixon AE, Baker SJ, Benkovic SJ. Using an AraC-based three-hybrid system to detect biocatalysts in vivo. Nat Biotechnol 2000; 18:544-7. [PMID: 10802623 DOI: 10.1038/75414] [Citation(s) in RCA: 45] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Recent methods to create large libraries of proteins have greatly advanced the discovery of proteins with novel functions. However, one limitation in the discovery of new biocatalysts is the screening or selection methods employed to find enzymes from these libraries. We have developed a potentially general method termed QUEST (QUerying for EnzymeS using the Three-hybrid system), which allows the construction of an easily screened or selected phenotype for, in theory, any type of enzymatic reaction. The method couples the in vivo concentration of an enzyme's substrate to changes in the transcriptional level of a reporter operon. Using the arabinose operon activator AraC, we constructed a system capable of detecting the fungal enzyme scytalone dehydratase (SD) in bacteria, and demonstrated its sensitivity and usefulness in library screening.
Collapse
Affiliation(s)
- S M Firestine
- Department of Chemistry, The Pennsylvania State University, University Park, PA 16802, USA
| | | | | | | | | |
Collapse
|
21
|
Liao DI, Basarab GS, Gatenby AA, Jordan DB. Selection of a potent inhibitor of trihydroxynaphthalene reductase by sorting disease control data. Bioorg Med Chem Lett 2000; 10:491-4. [PMID: 10743955 DOI: 10.1016/s0960-894x(00)00037-8] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
Compounds that control rice blast, but not other crop diseases, were selected for testing as inhibitors of trihydroxynaphthalene reductase of the fungal melanin biosynthetic pathway. A potent inhibitor of the enzyme (2) (Ki = 25 nM) was identified. An X-ray structure of the enzyme-NADPH-2 complex was determined at 2.1 A resolution.
Collapse
Affiliation(s)
- D I Liao
- DuPont Central Research and Development, Experimental Station, Wilmington DE, 19880, USA
| | | | | | | |
Collapse
|
22
|
Jordan DB, Basarab GS. Binding dynamics of two water molecules constrained within the scytalone dehydratase binding pocket. Bioorg Med Chem Lett 2000; 10:23-6. [PMID: 10636235 DOI: 10.1016/s0960-894x(99)00586-7] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
Two water molecules reside between inhibitors and active site residues of scytalone dehydratase. A molecular dynamics study is consistent with one water molecule binding less tightly than the other. Inhibitor binding studies with site-directed mutants indicate that the hydrogen bonding network around the less mobile water molecule contributes much greater binding energy than that around the more mobile one.
Collapse
Affiliation(s)
- D B Jordan
- E. I. DuPont de Nemours Agricultural Products, Stine-Haskell Research Center, Newark, DE 19714, USA
| | | |
Collapse
|
23
|
Basarab GS, Jordan DB. Wild-type enzyme as a reporter of inhibitor binding by catalytically impaired mutant enzymes. Biochem Biophys Res Commun 1999; 263:617-20. [PMID: 10512727 DOI: 10.1006/bbrc.1999.1423] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
A method for the determination of inhibition constants for catalytically-debilitated mutant enzymes is described. The inhibitor is partitioned between the mutant and wild-type enzymes. Catalytic rates of the wild-type enzyme are used as the signal of inhibitor binding to the mutant enzyme. The method is validated with scytalone dehydratase, the Y50F mutant, and a potent inhibitor. The K(i) value for Y50F determined by this method is 0.49 +/- 0.10 nM. The K(i) value determined using the Y50F catalytic report for inhibitor binding in the absence of wild-type enzyme is 0.20 +/- 0.030 nM. The wild-type enzyme binds the inhibitor ten-fold less tightly, thus indicating that the hydrogen-bonding interaction between the Y50 hydroxyl group and the inhibitor (suggested by X-ray crystallography) is weak. The method is most useful when the catalytic activity of the wild-type enzyme is the most sensitive report of inhibitor binding and the mutant enzyme is greatly crippled in catalytic activity.
Collapse
Affiliation(s)
- G S Basarab
- DuPont Agricultural Products, Stine-Haskell Research Center, 1094 Elkton Road, Newark, Delaware, 19714-0030, USA
| | | |
Collapse
|
24
|
Jennings LD, Wawrzak Z, Amorose D, Schwartz RS, Jordan DB. A new potent inhibitor of fungal melanin biosynthesis identified through combinatorial chemistry. Bioorg Med Chem Lett 1999; 9:2509-14. [PMID: 10498198 DOI: 10.1016/s0960-894x(99)00424-2] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
A new fungicide lead has been identified by in vitro screening of a focused combinatorial library. Amides (768) were synthesized in pools of four and assayed as inhibitors of scytalone dehydratase. Deconvolution of one of the most active pools led to the discovery of a potent inhibitor of the enzyme 3b (K(i) = 26 pM), which has fungicidal properties.
Collapse
Affiliation(s)
- L D Jennings
- Dupont Agricultural Products, Stine Haskell Research Center, Newark, DE 19714, USA
| | | | | | | | | |
Collapse
|