Pharmacophore Modeling for Anti-Chagas Drug Design Using the Fragment Molecular Orbital Method

Abstract

Background

Chagas disease, caused by the parasite Trypanosoma cruzi, is a neglected tropical disease that causes severe human health problems. To develop a new chemotherapeutic agent for the treatment of Chagas disease, we predicted a pharmacophore model for T. cruzi dihydroorotate dehydrogenase (TcDHODH) by fragment molecular orbital (FMO) calculation for orotate, oxonate, and 43 orotate derivatives.

Methodology/Principal Findings

Intermolecular interactions in the complexes of TcDHODH with orotate, oxonate, and 43 orotate derivatives were analyzed by FMO calculation at the MP2/6-31G level. The results indicated that the orotate moiety, which is the base fragment of these compounds, interacts with the Lys43, Asn67, and Asn194 residues of TcDHODH and the cofactor flavin mononucleotide (FMN), whereas functional groups introduced at the orotate 5-position strongly interact with the Lys214 residue.

Conclusions/Significance

FMO-based interaction energy analyses revealed a pharmacophore model for TcDHODH inhibitor. Hydrogen bond acceptor pharmacophores correspond to Lys43 and Lys214, hydrogen bond donor and acceptor pharmacophores correspond to Asn67 and Asn194, and the aromatic ring pharmacophore corresponds to FMN, which shows important characteristics of compounds that inhibit TcDHODH. In addition, the Lys214 residue is not conserved between TcDHODH and human DHODH. Our analysis suggests that these orotate derivatives should preferentially bind to TcDHODH, increasing their selectivity. Our results obtained by pharmacophore modeling provides insight into the structural requirements for the design of TcDHODH inhibitors and their development as new anti-Chagas drugs.

About the Authors:

Ryunosuke Yoshino

Contributed equally to this work with: Ryunosuke Yoshino, Nobuaki Yasuo, Daniel Ken Inaoka, Yohsuke Hagiwara

Affiliations Global Scientific Information and Computing Center, Tokyo Institute of Technology, Meguro, Tokyo, 152–8550, Japan, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo, Tokyo, 113–8657, Japan

Nobuaki Yasuo

Contributed equally to this work with: Ryunosuke Yoshino, Nobuaki Yasuo, Daniel Ken Inaoka, Yohsuke Hagiwara

Affiliation: Department of Computer Science, Tokyo Institute of Technology, Meguro, Tokyo, 152–8550, Japan

Daniel Ken Inaoka

Contributed equally to this work with: Ryunosuke Yoshino, Nobuaki Yasuo, Daniel Ken Inaoka, Yohsuke Hagiwara

Affiliation: Graduate School of Medicine, The University of Tokyo, Bunkyo, Tokyo, 113–0033, Japan

Yohsuke Hagiwara

Contributed equally to this work with: Ryunosuke Yoshino, Nobuaki Yasuo, Daniel Ken Inaoka, Yohsuke Hagiwara

Affiliation: Chemistry Research Labs, Drug Discovery Research, Astellas Pharma Inc., Tsukuba, Ibaraki, 305–8585, Japan

Kazuki Ohno

¶‡ These authors also contributed equally to this work.

Affiliation: Chemistry Research Labs, Drug Discovery Research, Astellas Pharma Inc., Tsukuba, Ibaraki, 305–8585, Japan

Masaya Orita

Affiliations Global Scientific Information and Computing Center, Tokyo Institute of Technology, Meguro, Tokyo, 152–8550, Japan, Chemistry Research Labs, Drug Discovery Research, Astellas Pharma Inc., Tsukuba, Ibaraki, 305–8585, Japan

Masayuki Inoue

¶‡ These authors also contributed equally to this work.

Affiliation: Graduate School of Pharmaceutical Sciences, The University of Tokyo, Bunkyo, Tokyo, 113–0033, Japan

Tomoo Shiba

Affiliation: Graduate School of Science and Technology, Kyoto Institute of Technology, Sakyo, Kyoto, 606–8585, Japan

Shigeharu Harada

Affiliation: Graduate School of Science and Technology, Kyoto Institute of Technology, Sakyo, Kyoto, 606–8585, Japan

Teruki Honma

Affiliation: Graduate School of Medicine, The University of Tokyo, Bunkyo, Tokyo, 113–0033, Japan

Emmanuel Oluwadare Balogun

Affiliations Graduate School of Medicine, The University of Tokyo, Bunkyo, Tokyo, 113–0033, Japan, Department of Biochemistry, Ahmadu Bello University, Zaria, 2222, Nigeria

Josmar Rodrigues da Rocha

Affiliation: Instituto de Química de São Carlos, Universidade de São Paulo, São Carlos, 13566–590, Brazil

Carlos Alberto Montanari

Affiliation: Instituto de Química de São Carlos, Universidade de São Paulo, São Carlos, 13566–590, Brazil

Kiyoshi Kita

Affiliation: Graduate School of Medicine, The University of Tokyo, Bunkyo, Tokyo, 113–0033, Japan

Masakazu Sekijima

* E-mail: [email protected]

Affiliations Global Scientific Information and Computing Center, Tokyo Institute of Technology, Meguro, Tokyo, 152–8550, Japan, Department of Computer Science, Tokyo Institute of Technology, Meguro, Tokyo, 152–8550, Japan

Introduction

Chagas disease is an infectious disease caused by the parasitic protozoan Trypanosoma cruzi (T. cruzi) and is the third most common parasitic disease in the world [1]. It affects people from approximately 20 countries, particularly those living in the southern United States and Latin America [2–4], with 15 million people estimated to be infected [5]. T. cruzi is primarily transmitted by blood-sucking insects belonging to the subfamily Triatominae (family Reduviidae) or by infected blood transfusion [6]. Once infected, the host may experience influenza-like symptoms during the acute phase and gastrointestinal and cardiac disorders during the chronic phase [7–8]. Two drugs, nifurtimox and benznidazole, are currently available for the treatment of Chagas disease, but there are severe problems associated with their use, including adverse effects and limited effectiveness during the chronic phase of the disease [9–10]. Thus, developing new therapeutic agents against T. cruzi infection is desirable [11–12].

To develop a novel anti-Chagas drug, we focused on dihydroorotate dehydrogenase (DHODH) as the target protein. DHODH is an enzyme that takes part in the fourth step in the de novo biosynthesis of pyrimidines, which are heterocyclic compounds essential for RNA and DNA synthesis. This enzyme is an oxidoreductase that catalyzes the oxidation of dihydroorotate to orotate using flavin mononucleotide (FMN) as a cofactor. FMN can take either an oxidized form (FMN) or a reduced form (FMNH2), and the oxidized form serves as the oxidizing agent during orotate production. FMNH2 is re-oxidized by an electron acceptor that differs according to the cellular localization of DHODH [13]. In humans, DHODH is a mitochondrial inner-membrane protein that uses respiratory ubiquinone as the electron acceptor [14]. In contrast, T. cruzi DHODH (TcDHODH), a cytosolic protein, uses fumarate as the electron acceptor [15]. A previous study showed that a DHODH-knockout T. cruzi was not viable [16]. The differences in biochemical properties between human and T. cruzi DHODHs as well as its essentiality for the parasite make TcDHODH a promising target for developing novel therapeutic agents against Chagas disease.

DHODH is a validated drug target for humans [17–18], as an immunosuppressant and also against Plasmodium falciparum [17, 19, 20] and Helicobacter pylori [21]. Species-specific DHODH inhibitors have been developed and shown to be effective in vitro [22] and in vivo [23]. However, all inhibitors developed to date target the ubiquinone binding site and do not inhibit the cytosolic DHODHs, for which potent and selective inhibitors have never been reported.

The atomic resolution crystal structures of TcDHODH in complexes with its substrates and products have been determined [15]. Based on analysis of overall structure and the active site region of the TcDHODH-orotate complex (Fig 1), it is thought that a strong π-π interaction between orotate and the isoalloxazine ring of FMN occurs.

[Figure omitted. See PDF.]

Fig 1. TcDHODH overall structure.

Crystal structure of TcDHODH (PDB ID: 2E6A, A-chain).

https://doi.org/10.1371/journal.pone.0125829.g001

A pharmacophore is defined as “an ensemble of steric and electronic features that ensures optimal supramolecular interactions with a specific biological target and the trigger (or inhibit) of its biological function” [24]. Based on this definition, we define pharmacophore modeling as a process for predicting pharmacophores with common or specific characteristics among compounds. This definition is applied not only to molecular design but to protein—ligand docking simulation and quantitative structure-activity relationships (QSAR) as well [25]. However, pharmacophore modeling without ligand structural alignment information is difficult. Thus, knowledge of protein—ligand structure is useful for predicting pharmacophores.

The fragment molecular orbital (FMO) method [26] employs ab initio quantum mechanical calculations for large biomolecules such as protein—ligand complexes. Intermolecular interaction energies typically can be determined on the basis of molecular mechanics. However, this method is not universally applicable to all compounds, because there is a limit to the determination of molecular potentials based on atom type, especially of quantum chemical elements such as π electrons. For this reason, in this study, we used the FMO method to analyze the interaction energies between the target proteins and ligands with the aim of identifying important amino acid residues for ligand binding. Amino acids and ligands in the system of interest are divided into fragments, and molecular orbital calculations are performed for individual fragments. Because the effects of interfragment potentials are taken into account in these molecular orbital calculations, the FMO method can estimate the interaction energy between each pair of fragments. The method can clearly describe the detailed interactions between the ligand and each amino acid residue, and is frequently used in the design of new drugs [27–34]. Moreover, the method can extract specific interaction from a wide variety of derivatives. FMO calculation is thus suitable for obtaining pharmacophore models and is useful for guiding molecular design.

In the present study, we identified pharmacophores using the FMO method with the aim of designing anti-Chagas drugs, via analysis of interaction energy between TcDHODH and orotate, oxonate [35] as a competitive inhibitor of DHODH, and 43 orotate derivatives.

Materials and Methods

Co-crystallization of TcDHODH with orotate derivatives

Recombinant TcDHODH expression, purification, and crystallization in complex with 43 orotate derivatives were performed essentially as previously reported [15]. Briefly, TcDHODH was purified to homogeneity from BL21(DE3)pET3a/TcDHODH using DEAE Sepharose Fast Flow (GE Healthcare) followed by Phenyl Sepharose H.P. (GE Healthcare) and TSK G3000SW (Tosoh). By this method, a total of 10–15 mg of TcDHODH with specific activity ranging from 12 to 18 μmol/min/mg could be purified from 10 L of culture. TcDHODH crystals formed in the presence of oxonate [15] were used to soak overnight into crystallization buffer containing 1 mM of orotate derivatives that are listed in Fig 2. X-ray diffraction data were collected at SPring-8 (beam lines BL32XU, BL41XU and BL44XU) or Photon Factory (beam lines AR-NW12A, AR-NE3A, BL-5A and BL-17A). All the data were processed and scaled with HKL2000 [36]. The co-crystal structures of TcDHODH in complex with orotate derivatives were initially solved by molecular replacement using coordinates from TcDHODH-oxonate complexed structure (3W1Q) and later by one of the open form structures (e.g., 3W1R). Manual building and crystallographic refinement were performed with the programs COOT [37] and REFMAC5 [38], respectively. The PDB IDs, resolution, R-value, R-Free, average of the occupancy of the ligands and average of the B-factors are listed in Table 1. The detailed synthetic process of the 43 orotate derivatives and the X-ray analysis performed in this study will be described elsewhere (Inaoka et al., manuscript in preparation).

[Figure omitted. See PDF.]

Fig 2. Structures of orotate and R groups of synthesized derivatives.

Functional groups were introduced at the pyrimidine 5-position through an ethylene or butylene linker.

https://doi.org/10.1371/journal.pone.0125829.g002

[Figure omitted. See PDF.]

Table 1. Summary of co-crystal structures (UniProtKB: Q4D3W2).

https://doi.org/10.1371/journal.pone.0125829.t001

Interaction energy analysis

Interaction energy analysis was performed using the analytical tool Facio [39] based on pair interaction energy decomposition analysis, as proposed by Fedorov et al. [40]. This analysis using the FMO method provides a quantitative evaluation of hydrogen bonding and hydrophobic interactions that are important for ligand binding to a protein, as well as π-π, π–cation, and CH–π interactions, which require quantum chemical calculations [41]. These analyses are sometimes applied in structure-based drug design [42].

Calculation procedure

All structures of the TcDHODH—compound complexes were visualized and computations were performed with hydrogenation. Optimization of the structures was performed only for the added hydrogen atoms, with all heavy atoms fixed at the positions given in the PDB using the CHARMM force field [43] by Discovery Studio (Accelrys, San Diego, CA). FMO calculation job files were generated using FMOutil version 2.1, and calculations were performed for each A-chain monomer using GAMESS [44] at the MP2/6-31G level. Alternate positions at except active site adopted a type-A conformation. Alternate positions in the active site adopted both. All calculations were performed with the TSUBAME2.5 supercomputer at Tokyo Tech (a HP Proliant SL390s G7 server with an Intel Xeon X5670 2.93-GHz processor). The approximate calculation time of each structure was 3 h at the MP2/6-31G level.

Definition of pharmacophore

In this paper, we define a common pharmacophore as that having characteristics common to more than 80% conserved interactions between TcDHODH and compounds, and a particular pharmacophore as that having characteristics that have interactions stronger than −10 kcal mol-1 between TcDHODH and compounds.

Results

Analysis of interaction energies between TcDHODH and orotate and oxonate

Orotate binds to the active center of TcDHODH, which is positioned parallel to an aromatic ring of FMN. Certain amino acid residues such as Lys43, Asn67, and Asn194 form hydrogen bonds with orotate and are important for substrate binding (Fig 3).

[Figure omitted. See PDF.]

Fig 3. Orotate binding mode.

Crystal structure of the TcDHODH—orotate complex active site. Hydrogen bond distance with the Lys43 amine group: 2.93 Å, with the Asn67 amine group: 3.05 Å, with the Asn67 carbonyl group: 2.82, with the Asn194 amine group: 3.32 and with the Asn194 carbonyl group: 2.92.

https://doi.org/10.1371/journal.pone.0125829.g003

The side chain of Asn67 interacts with the 1-position hydrogen atom and 2-position carbonyl group of orotate, while the side chain of Asn194 side chain interacts with orotate’s 3-position hydrogen atom and 4-position carbonyl group through two hydrogen bonds. Moreover, the Lys43 side chain amide interacts with orotate’s 6-position carboxylate group.

The FMO-based interaction energy analysis indicates that the hydrogen bonds with Lys43, Asn67 and Asn194 residues are important in this interaction (-13.02, -33.71, and -24.79 kcal mol-1, respectively), with interactions with Asn67 being strongest (Fig 4).

[Figure omitted. See PDF.]

Fig 4. Interaction energy analysis of TcDHODH—orotate.

The vertical axis shows the interaction energy between the ligand and each fragment, and the horizontal axis shows the fragment numbers. Fragments are numbered in order from N to C termini, followed by ligands, such as substrate and cofactor, to preserve the order described in the PDB file.

https://doi.org/10.1371/journal.pone.0125829.g004

These findings indicate that Asn67 and Asn194 contribute to the specificity of orotate’s pyrimidine ring, and Lys43 contributes to the hydrogen bonding of orotate at the carboxylate moiety at the pyrimidine 6-position. The π-π interaction energy between orotate and FMN was −6.29 kcal mol-1 and contributed significantly to ligand binding (Fig 4). Given that Asn67 and Asn194 residues form two hydrogen bonds with orotate, large interaction energies would be predicted from this binding mode.

Oxonate, a well-known inhibitor of TcDHODH [15], also binds to the active center of TcDHODH, similarly to orotate. The FMO-based interaction energy analysis confirms the interaction with Lys43, Asn67, Asn194 and FMN as well as that of orotate (Fig 5: -7.41, -43.73, -16.21, and -7.26 kcal mol-1, respectively).

[Figure omitted. See PDF.]

Fig 5. Interaction energy analysis of TcDHODH—oxonate.

The vertical axis and horizontal axis are same as the Fig 4. (PDB ID: 2E6F, A-chain).

https://doi.org/10.1371/journal.pone.0125829.g005

These results suggest that inhibitors targeting the active site of TcDHODH should preferentially have an aromatic moiety containing hydrogen bond donors and acceptors.

Analysis of interaction energies between TcDHODH and orotate derivatives

Next, 43 orotate derivatives (Fig 2) predicted to form additional hydrogen bond and hydrophobic interactions were synthesized and their co-crystal structures were determined. These derivatives contain a functional group substituted at the pyrimidine 5-position through an ethylene or butylene linker. The binding mode between the active site of TcDHODH and derivative 27 is shown in Fig 6 as an example. Similarly to orotate, derivative 27 and other derivatives interact with Lys43, Asn67, Asn194 and FMN.

[Figure omitted. See PDF.]

Fig 6. Binding mode of derivative 27.

Crystal structures of the TcDHODH—derivative 27 complex active site. Hydrogen bond distance with the Lys43 amine group: 2.73 Å, with the Asn67 amine group: 2.89 Å, with the Asn67 carbonyl group: 2.96, with the Asn194 amine group: 2.86, with the Asn194 carbonyl group: 2.80, with the Lys214 residue: 2.91 Å.

https://doi.org/10.1371/journal.pone.0125829.g006

Furthermore, new interactions were predicted from the crystal structure, such as Lys214 by the introduction of a functional group at the pyrimidine 5-position.

These orotate derivatives were used to analyze the interaction energy by the FMO method. Table 2 shows the results of the interaction energy analyses for some amino acid residues and FMN. Lys43, Asn67, Asn194, and FMN showed strong interactions with the orotate moiety.

[Figure omitted. See PDF.]

Table 2. Interaction energies between TcDHODH and orotate derivatives (kcal mol-1).

https://doi.org/10.1371/journal.pone.0125829.t002

The interaction energies between three amino acid residues and all orotate derivatives were strongly similar to those of orotate. Our analyses indicated that interactions between the orotate moiety and these three amino acid residues as well as FMN were strongly conserved for all of the derivatives. In contrast, some derivatives such as the alternative conformer B of 21 showed the strongest interaction with Lys214, followed by 27, 5, 10, and 8 with interaction energies of -16.50, -13.55, -11.89, -7.48, and -7.13 kcal mol-1, respectively. The carboxylate group introduced in the phenyl moiety of 27 formed a hydrogen bond with the Lys214 residue (Fig 6). Derivatives 40 and 41 interacted with Lys214 with low energy despite the introduction of a carboxylate group on the aromatic ring moiety. Derivatives 40 and 41 were linked to a functional group at the pyrimidine 5-position through a butylene linker. In contrast, derivatives 21 and 27 interacted strongly (>10 kcal mol-1) when coupled with functional groups through an ethylene linker. These results suggest that the ethylene linker is appropriate for mediating interactions with Lys214.

Pharmacophore obtained from calculation

From these results, intermolecular interactions with the Lys43, Asn67 and Asn194 by hydrogen bond and FMN through π-π interaction were obtained as a common pharmacophore in the orotate moiety and intermolecular interactions with the Lys214 by hydrogen bond obtained as particular pharmacophore at derivatives 5, 21 and 27.

Discussion

This study provides a comprehensive FMO-based interaction analysis between TcDHODH residues and 43 orotate derivative inhibitors. Our analysis indicated that the orotate moiety formed hydrogen bonds with the Lys43, Asn67, and Asn194 residues of TcDHODH and a π-π interaction with FMN. These interactions were conserved in oxonate and orotate derivatives, and are also consistent with the magnitude of the interaction energy. From these results, the orotate moiety of orotate derivatives maintained these interactions whether or not functional groups were introduced at the pyrimidine ring 5-position. We accordingly expected a pharmacophore model containing molecules with hydrogen-bond donor/acceptor groups corresponding to Lys43, Asn67, and Asn194, such as amine, carbonyl, and carboxylate groups with an aromatic ring close to FMN.

In contrast, some derivatives containing an acceptor at a functional group form a hydrogen bond with the Lys214 residue. In particular, the energies of interaction between Lys214 and derivatives 5, 21, 27 were predicted to be stronger than -10 kcal mol-1. Furthermore, alignment analysis has shown that the TcDHODH Lys214 residue is replaced by Arg298 in human DHODH [15]. In their crystal structures, TcDHODH Lys214 and human DHODH Arg298 are each located in the loop connecting βF–βG and β6–βE [14–15]. These loop regions are poorly conserved within the family of DHODHs, an observation consistent with the difference in the loop structures (Fig 7); whereas, the amino group of TcDHODH Lys214 points towards the active site, the guanidium group of human DHODH Arg298 points to the opposite side, with a distance of 20 Å between the two.

[Figure omitted. See PDF.]

Fig 7. Comparison of loop structures.

Alignment structure based on orotate (green ribbon: 2E6A, blue ribbon: 1D3G, distance between TcDHODH Lys214 and human DHODH Arg298: 19.5 Å). Superimposition was performed by “Align Binding Sites” in Maestro based on orotate.

https://doi.org/10.1371/journal.pone.0125829.g007

Thus, even if TcDHODH Lys214 and human DHODH Arg298 are apparently conserved by alignment analysis, they are structurally not conserved. Consequently, orotate derivatives such as derivatives 5, 21, and 27 may specifically inhibit TcDHODH by interacting with the Lys214 residue. Thus, a hydrogen bond acceptor as a particular pharmacophore corresponding to Lys214 is necessary to inhibit TcDHODH selectively.

In conclusion, we propose a pharmacophore model inferred from FMO-based interaction analysis. Hydrogen bond acceptor pharmacophores correspond to Lys43 and Lys214, hydrogen bond donor and acceptor pharmacophores correspond to Asn67 and Asn197, and an aromatic ring pharmacophore corresponds to FMN, indicating important characteristics of compounds expected to inhibit Trypanosoma cruzi DHODH (Fig 8).

[Figure omitted. See PDF.]

Fig 8. Pharmacophore inferred from interaction energy analysis.

Acc.: Acceptor, Don.: Donor and Aro.: Aromatic.

https://doi.org/10.1371/journal.pone.0125829.g008

These characteristics inferred from single-point FMO calculations for a large number of structures provide insights into the ligand—amino acid residue interactions important for pharmacophore modeling and may facilitate the development of TcDHODH—targeted anti-Chagas drugs.

Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable comments and suggestions to improve the quality of the paper. We also thank all beamline staff members from SPring-8 as well as Photon Factory for their help with X-ray diffraction data collection. This research was partially supported by Japanese Society for the Promotion of Science (JSPS) KAKENHI Grant Numbers 25870222 and 24240044 (to M.S.). This work was supported in part by Creative Scientific Research Grant 18GS0314 (to K.K.), Grant-in-aid for Scientific Research 18073004 and 26253025 (to K.K.) and 26870119 (to D.K.I.) from the JSPS, and Targeted Proteins Research Program (to K.K. and S.H.) from the Japanese Ministry of Education, Science, Culture, Sports and Technology (MEXT). We also acknowledge the support from Science and Technology Research Promotion Program for Agriculture, Forestry, Fisheries and Food Industry (to K.K. and S.H.) and JST/JICA, SATREPS (Science and Technology Research Partnership for Sustainable Development) (10000284) to K.K. Astellas Pharma Inc., provided support in the form of salaries for authors YH, KO and MO.

Author Contributions

Conceived and designed the experiments: RY DKI KO MO SH TH KK MS. Performed the experiments: RY NY YH MI TS EOB JRR CAM. Analyzed the data: RY DKI YH EOB JRR CAM. Contributed reagents/materials/analysis tools: DKI EOB JRR CAM. Wrote the paper: RY NY DKI KO MO MI SH TH KK MS.

Citation: Yoshino R, Yasuo N, Inaoka DK, Hagiwara Y, Ohno K, Orita M, et al. (2015) Pharmacophore Modeling for Anti-Chagas Drug Design Using the Fragment Molecular Orbital Method. PLoS ONE 10(5): e0125829. https://doi.org/10.1371/journal.pone.0125829

References

1. Lockman JW, Hamilton AD. Recent developments in the identification of chemotherapeutics for Chagas disease. Curr Med Chem. 2005;12: 945–959. pmid:15853707

2. Schmunis GA, Yadon ZE. Chagas disease: a Latin American health problem becoming a world health problem. Acta Trop. 2010;115: 14–21. pmid:19932071

3. Roellig DM, Brown EL, Barnabe C, Tibayrenc M, Steurer FJ, Yabsley MJ. Molecular typing of Trypanosoma cruzi isolates. United States. Emerg Infect Dis. 2008;14: 1123–1125. pmid:18598637

4. Marcet PL, Duffy T, Cardinal MV, Burgos JM, Lauricella MA, Levin MJ, et al. PCR-based screening and lineage identification of Trypanosoma cruzi directly from faecal samples of triatomine bugs from northwestern Argentina. Parasitology. 2006;132: 57–65. pmid:16393354

5. Rassi A Jr, Rassi A, Marin-Neto JA. Chagas disease. Lancet. 2010;375: 1388–1402. pmid:20399979

6. Chagas C. Nova trypanozomiaze humana: estudos sobre a morfolojia e o ciclo evolutivo do Schizotrypanum cruzi n. gen., n. sp., ajente etiolojico de nova entidade morbida do homem. Mem Inst Oswaldo Cruz. 1909;1: 159–218.

7. Kirchhoff LV. American trypanosomiasis (Chagas’ disease). Gastroenterol Clin North Am. 1996;25: 517–533. pmid:8863038

8. Rossi MA, Bestetti RB. The challenge of chagasic cardiomyopathy. Cardiology. 1995;86: 1–7. pmid:7728781

9. Hugo VDB. Chemotherapy for parasitic infections. Nature. 1978;273: 626–630. pmid:351413

10. Linares GEG, Ravaschino EL, Rodriguez JB. Progresses in the field of drug design to combat tropical protozoan parasitic diseases. Curr Med Chem. 2006;13: 335–360. pmid:16475941

11. Maya JD, Cassels BK. Iturriaga-Vasquez P, Ferreira J, Faundez M, Galanti N, et al. Mode of action of natural and synthetic drugs against Trypanosoma cruzi and their interaction with the mammalian host. Comp Biochem Physiol A Mol Integr Physiol. 2007;146: 601–620. pmid:16626984

12. Buckner FS, Navabi N. Advances in Chagas disease drug development: 2009–2010. Curr Opin Infect Dis. 2010;23: 609–616. pmid:20885320

13. Jensen KF, Björnberg O. Evolutionary and functional families of dihydroorotate dehydrogenases. Paths Pyrimidines. 1998;6:20–28.

14. Liu S, Neidhardt EA, Grossman TH, Ocain T, Clardy J. Structures of human dihydroorotate dehydrogenase in complex with antiproliferative agents. Structure. 2000;8: 25–33. pmid:10673429

15. Inaoka DK, Sakamoto S, Shimizu H, Shiba T, Kurisu G, Nara T, et al. Structures of Trypanosoma cruzi dihydroorotate dehydrogenase complexed with substrates and products: atomic resolution insights into mechanisms of dihydroorotate oxidation and fumarate reduction. Biochemistry. 2008;47: 10881–10891. pmid:18808149

16. Annoura T, Nara T, Makiuchi T, Hashimoto T, Aoki T. The origin of dihydroorotate dehydrogenase genes of kinetoplastids, with special reference to their biological significance and adaptation to anaerobic, parasitic conditions. J Mol Evol. 2005;60: 113–127. pmid:15696374

17. Davies M, Heikkilä T, McConkey GA, Fishwick CW, Parsons MR, Johnson AP. Structure-based design, synthesis, and characterization of inhibitors of human and Plasmodium falciparum dihydroorotate dehydrogenases. J Med Chem. 2009;52: 2683–2693. pmid:19351152

18. Davis JP, Cain GA, Pitts WJ, Magolda RL, Copeland RA. The immunosuppressive metabolite of leflunomide is a potent inhibitor of human dihydroorotate dehydrogenase. Biochemistry. 1996;35: 1270–1273. pmid:8573583

19. Coteron JM, Marco M, Esquivias J, Deng X, White KL, White J, et al. Structure-guided lead optimization of triazolopyrimidine-ring substituents identifies potent Plasmodium falciparum dihydroorotate dehydrogenase inhibitors with clinical candidate potential. J Med Chem. 2011;54: 5540–5561. pmid:21696174

20. Phillips MA, Rathod PK. Plasmodium dihydroorotate dehydrogenase: a promising target for novel anti-malarial chemotherapy. Infect Disord Drug Targets. 2010;10: 226–239. pmid:20334617

21. Copeland RA, Marcinkeviciene J, Haque TS, Kopcho LM, Jiang W, Wang K, et al. Helicobacter pylori-selective antibacterials based on inhibition of pyrimidine biosynthesis. J Biol Chem. 2000;275: 33373–33378. pmid:10938275

22. Diao Y, Lu W, Jin H, Zhu J, Han L, Xu M, et al. Discovery of Diverse Human Dihydroorotate Dehydrogenase Inhibitors as Immunosuppressive Agents by Structure-Based Virtual Screening. J Med Chem. 2012;55: 8341–8349. pmid:22984987

23. Marwaha A, White J, El Mazouni F, Creason SA, Kokkonda S, Buckner FS, et al. Bioisosteric Transformations and Permutations in the Triazolopyrimidine Scaffold To Identify the Minimum Pharmacophore Required for Inhibitory Activity against Plasmodium falciparum Dihydroorotate Dehydrogenase. J Med Chem. 2012;55: 7425–7436 pmid:22877245

24. Wermuth CG, Ganellin CR, Lindberg P, Mitscher LA. Glossary of terms used in medicinal chemistry (IUPAC Recommendations 1998). Pure and Applied Chemistry. 1998;70: 1129–1143.

25. Lingling J, Robert CR. Pharmacophore-based similarity scoring method for DOCK. J Phys Chem B. 2015 Jan 22;119(3):1083–102. . Epub.2014.10.10. pmid:25229837

26. Kitaura K, Sawai T, Asadaa T, Nakanob T, Uebayasi M. Pair interaction molecular orbital method: an approximate computational method for molecular interactions. Chem Phys Lett. 1999;312: 319–324.

27. Tomonaga O, Kosuke O. CH/π hydrogen bonds determine the selectivity of the Src homology 2 domain to tyrosine phosphotyrosyl peptides: an ab initio fragment molecular orbital study. Journal of Computational Chemistry. 2008;29: 2656–2666. pmid:18484636

28. Tomonaga O, Eiichi T, Motoyasu O, Cchiaki H, Harunobu M, Tochihiro N,et al. The importance of CH/π hydrogen bonds in rational drug design: an ab initio fragment molecular orbital study to leukocyte-specific protein tyrosine (LCK) kinase. Bioorganic and Medicinal Chemistry. 2008;16: 10311–10318. pmid:18977146

29. Tomonaga O, Kosuke O, Kazuo K. CH/π hydrogen bonds play a role in ligand recognition and equilibrium between active and inactive states of the b2 adrenergic receptor: An ab initio fragment molecular orbital (FMO) study. Bioorganic and Medicinal Chemistry. 2011;19: 5231–5237. pmid:21820903

30. Tomonaga O, Kosuke O, Kazuo K. Importance of CH/π hydrogen bonds in recognition of the core motif in proline-recognition domains: an ab initio fragment molecular orbital study. Journal of Computational Chemistry. 2011;32: 2774–2782. pmid:21710635

31. Kenichi F, Yoshimasa S. The role of fluorine atoms in a fluorinated prostaglandin agonist. Chem Med Chem. 2010;5: 1254–1257. pmid:20486154

32. Ohno K., Mori K., Orita M., Takeuchi M.. Computational insights into binding of bisphosphates to farnesyl pyrophosphate synthase. Current Medicinal Chemistry. 2011;18: 220–233. pmid:21110804

33. Osamu I, John B, Richard JL, Mark W, Compound design by fragment linking. Molecular Informatics. 2011;30: 298–306.

34. Michael PM, Osamu I, Richard JL, Mark W, Prediction of cyclin dependent kinase 2 inhibitor potency using the fragment molecular orbital method. Journal of Cheminformatics. 2011;3: 2. pmid:21219630

35. Björnberg O., Jordan DB, Palfey BA, Jensen KF. Dihydrooxonate is a substrate of dihydroorotate dehydrogenase (DHOD) providing evidence for involvement of cysteine and serine residues in base catalysis. Arch. Biochem. Biophys. 2001;391: 286–94. pmid:11437361

36. Zbyszek O and Wladek M. Processing of X-ray Diffraction Data Collected in Oscillation Mode. Methods in Enzymology. 1997; 276: 307–326.

37. Emsley P and Cowtan K. Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr. 2004;60: 2126–2132. pmid:15572765

38. Murshudov GN, Vagin AA, Dodson EJ. Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr D Biol Crystallogr. 1997;53: 240–255. pmid:15299926

39. http://www1.bbiq.jp/zzzfelis/Facio.html

40. Fedorov DG, Kazuo K. Pair interaction energy decomposition analysis. J Comput Chem. 2007;28: 222–237. pmid:17109433

41. Ozawa T, Tsuji E, Ozawa M, Handa C, Mukaiyama H, Nishimura T, et al. The importance of CH/pi hydrogen bonds in rational drug design: An ab initio fragment molecular orbital study to leukocyte-specific protein tyrosine (LCK) kinase. Bioorg Med Chem. 2008;16: 10311–10318. pmid:18977146

42. Watanabe C, Fukuzawa K, Okiyama Y, Tsukamoto T, Kato A, Tanaka S, et al. Three- and four-body corrected fragment molecular orbital calculations with a novel subdividing fragmentation method applicable to structure-based drug design. J Mol Graph Model. 2013l;41: 31–42. pmid:23467020

43. David. AC, Tom D, Thomas EC, Carlos S, Junmei. W, Robert. ED, et al., Amber 11. San Francisco. 2010

44. Michael WS, Kim KB, Jerry AB, Steven TE, Mark SG, Jan HJ, et al. General atomic and molecular electronic structure system. J Comput Chem. 1993;14: 1347–1363.