Computational Study of the Interactions between Antimalarial Chemotherapies with Folate Pathway Receptors and Telomerase Reverse Transcriptase ()
1. Introduction
Malaria is one of the debilitating parasitic diseases caused by the protozoan parasite Plasmodium falciparum [1] [2]. It affects a large part of the world’s population and causes an average of half a million deaths each year [3]. The folate metabolic pathway is a critical source of metabolites and intermediates necessary for the synthesis of nucleotides and amino acids, paramount to the survival of the Plasmodium parasite [4] [5]. This pathway has been targeted for both treatment and prophylaxis of the disease [6]. Two key enzymes involved in this pathway, namely dihydrofolate reductase (DHFR) and dihydropteroate synthetase (DHPS) are highlighted in supplemental materials (Scheme 1).
Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate to tetrahydrofolate a key precursor in the synthesis of nucleotide bases (pyrimidines, and thymidylic acid) and amino acids, required by bacteria, Plasmodium as well as normal and cancerous human cells. Several antifolate drugs have been shown to prevent replication of parasitic DNA and therefore one effective way to fight malaria. This is primarily achieved by inhibiting DHFR enzymatic activity thereby blocking cell division of the parasite at the schizont stage [7] [8]. There is, however, evidence suggesting that the activity of antifolate drugs like chloroquine and pyrimethamine has been undermined by several cases of parasitic resistance caused by DHFR mutations [9]. The mutations affect the polarity and size of the DHFR binding domains. It is therefore vital to understand how such subtle changes affect the molecular details of the interactions. Dihydropteroate synthase (DHPS) is the other folate pathway enzyme that catalyzes the conversion of p-aminobenzoic acid to dihydrofolic acid. Targeting this enzyme is also crucial because it is only expressed in parasitic cells and therefore a good target for drug specificity. In addition to folate pathway enzymes, Plasmodium reverse transcriptase telomerase (pf-TERT) was recently identified as a target receptor for the development of antimalarial chemotherapies was discussed in a recent study by Parida and collaborators, Plasmodium falciparum telomerase reverse transcriptase (pf-TERT) also emerged as a plausible target receptor for antimalarial activity and his sequence provides valuable information for the design of specific anti-telomerase drugs [10].
In order to understand how these enzymes interact with antimalarial drugs at the molecular level requires the availability of good three-dimensional structures of the target receptors. The three-dimensional structures of some of these target receptors (DHFR) have been solved and available in the protein databank [11]. The absence of X-ray crystal structures for pf-DHPS and pf-TERT receptors pose a challenge.
In this study, homology modeling and molecular docking are employed, to characterize the interactions between key enzyme receptors in the critical pathways described above with current antimalarial medications. Specifically: 1) Homology modeling is used to build three-dimensional structures of Plasmodium enzymes (DHPS andTERT); 2) AutoDock Vina is used to estimating the binding affinity of antimalarial medications (chloroquine, pyrimethamine, sulfadoxine, pentamidine, primaquine and artemisinin) with human and parasitic receptors.
The main objective is to determine whether DHPS and pf-TERT can also serve as potential target receptors for antimalarial activity. In addition, we explore the difference in affinity between DHFR mutants and antimalarial drugs. The affinity and interaction patterns for multiple target enzyme receptors in different pathways can help explain the wide range of action of many antimalarial medications. Autodock Vina is used to rapidly access the stability and rank the drug-receptor complexes, with the goal to determine and explain any observed specificity and or selectivity of current antimalarial drugs like pyrimethamine. The drugs with enhanced selectivity will serve as pharmacophores for the development of novel sulfonamides, pyrimethamine or chloroquine-based derivatives.
2. Molecular Systems and Computational Methods
2.1. Molecular Systems
The interactions between 6 antimalarial drugs used in and target receptors from multiple metabolic pathways were characterized using molecular docking. The docked complexes were generated using a combination of the Chimera molecular visualization package, and docking software AutoDock Vina [12]. The three-dimensional structures of the antimalarial drugs (artemisinin, pyrimethamine, primaquine, chloroquine, pentamidine, and sulfadoxine) were downloaded from DrugBank database and represented in (Figure 1) [13].
2.2. Homology Modeling of Target Protein Receptors
Homology modeling is a computational technique used to approximate the tertiary structure of a protein whose primary amino acid sequence is available, based on the coordinates known 3D homologous structures. The atomic coordinates of the target protein is obtained based from the homologue structure based on the similarity scores obtained from the sequence and structural alignments. The four steps implemented to build the structures for pf-DHPS and pf-TERT include: target identifications, sequence alignment, model building and refinement [14] [15]. In this study, the blast sequence tool in Uniprot database was used to perform a sequence alignment between the target sequences (pf-TERT-uniprot ID 23479392 and pf-DHPS-uniprot ID Q5PZ24) and numerous template sequences
Figure 1. Chemical Structures of antimalarial ligands. (A) Pyrimethamine; (B) Chloroquine; (C) Primaquine; (D) Artemisinin; (E) Pentamidine and (F) Sulfadoxine.
TableSM1 [16] [17]. The crystal structure of Plasmodium falciparum Glutathione reductase, accession code (1ONF) was used to build pf-TERT model, while Plasmodium falciparum hydroxymethyldihydropterin pyrophosphokinase-dihydropteroate synthase accession code (6KCL) was used to build pf-DHPS). The models were built using Modeller with the Chimera graphics user interface [18]. The atomic overlaps in the resulting models were relaxed using the Chimera energy minimization tool. In each case, the AMBER force field was used with a dielectric constant of 78.3 (water). A strategy of steepest descents followed by conjugate gradients was adopted for a total of 10,000 iterations [19].
2.3. Docking Simulations
Two representative receptors from the folate pathway were used in these simulations. These include dihydrofolate reductase thymidylate synthase (DHFR-TS), and dihydropteroate synthase (DHPS). DHFR-TS mutations have been linked to resistance to antimalarial activity. The three-dimensional structures of DHFR have been solved and were downloaded from the RCSB Protein Data Bank are represented in (Figure 2). These include; human DHFR (h-DHFR and pdb-ID-3GYF), wild type Plasmodium DHFR (wt-pf-DHFR and pdb-ID-4DPD) double mutant Plasmodium (d-pf-DHFR and pdb-ID-3UM6) and quadruple mutant Plasmodium DHFR (q-DHFR and pdb-ID-4DP3) [20] [21] [22].
The target receptor enzymes with solved X-ray crystallized three-dimensional structures were refined prior to docking using Chimera. This involved adding deleted atoms (polar hydrogens) from crystal structures, removing all non-standard ligands and relaxing the enzyme structures. Other receptors used in the docking simulations include the homology models of pf-TERT and pf-DHPS (Figure 3).
The USCF Chimera Dock Prep tool was used to refine the receptor molecules and Gasteiger charges. A global search of the best bound complexes, ranked and
Figure 2. (A) Crystal Structure of human DHFR; (B) Wild type pf-DHFR, double mutant pf-DHFR with mutations A16V and S108N highlighted; and quadruple mutant pf-DHFR with mutations N51I, C59R, S108N, and I164L highlighted.
Figure 3. (A) Homology model of Plasmodium Teleromase; (B) Energy minimization profile; (C) homology model of plasmodium falciparum DHPS; and (D) Energy minimization profile.
scored based on favorable hydrophobic, hydrogen bond and van der Waals interactions was performed using the AutoDock Vina software [23]. Docking was performed using the protein and ligand coordinates along with docking grid box with size 50 Å × 50 Å × 50 Å center around active site of DHFR [24]. The docked complex with lowest binding energy is then extracted for further analysis.
3. Results and Discussion
3.1. Relative Binding Energies
The docking simulations screened the affinity and interactions of six current antimalarial medications with a multiple target receptors to understand what is driving these interactions at a molecular level, determine the key residues involved in binding and whether some receptor mutations significantly affect the binding process. The AutoDock Vina binding energies between the antimalarial drugs and folate pathway enzyme DHFR are represented in Figure 4. We also explored the interactions of pf-TERT and pf-DHPS using three dimensional structures obtained from homology modeling. All the antimalarial drugs modeled (pyrimethamine, primaquine, pentamidine, chloroquine, artemisinin and sulfadoxine) formed complexes with Plasmodium enzyme receptors.
The folate acid pathway enzymes were targeted because of their role in the synthesis of nucleotides and amino acids, critical for the survival of the plasmodium parasite. There are some indications that antimalarial drugs can interfere with the telomerase activities to induce shortening of telomere and chromosome loss [25]. This can be effective in fighting the malaria scourge caused by plasmodium parasite [26]. All ligands were successfully docked to the target DHFR enzyme receptors. Artemisinin, the benchmark drug present in many effective combination therapies against malaria interact more favorably with all the DHFR enzymes. There is some evidence indicating that the mutations in DHFR negatively impact the affinity of interactions involving artemisinin (Figure 4(A)). For example, the affinity becomes less negative from h-DHFR and wt-DHFR to quadruple mutant pf-DHFR (qm-pf-DHFR) (Figure 4(A)). The trend is slightly reversed for the interactions between pyrimethamine and the Plasmodium enzymes. For example, there is a bigger difference in energy between h-DHFR, or wt-pf-DHFR and dm-pf-DHFR which are respectively of −1.5 kcal/mol and −1.1 kcal/mol contrary to qm-pf-DHFR whose difference is only −0.4 kcal/mol with h-DHFR. The results of (Figure 4(A)) are different from those obtained by Forlemu and coworkers [5]. This difference could be explained not only by the active site chosen to perform the docking calculation but above all the nature of the amino acids (polar and non-polar) which involved these sites. The antimalarial drugs were also docked to thepf-DHPS and pf-TERT models built using
Figure 4. Relative binding energies between antimalarial drugs and folate enzymes. (A) DHFR receptors and (B) pf-TERT and pf-DHPS.
homology modeling (Figure 4(B)). The strength of these interactions does vary from receptor to receptor and depended on ligands applied. Specifically the affinity of the ligands with folate pathway enzyme DHFR-TS dimer is stronger compared to pf-DHPS and pf-TERT. The benchmark drug used to treat malaria formed the most stable complexes with all the receptors. The affinity of all the ligands with folate pathway enzyme DHFR is more favorable compared to DHPS and pf-TERT. For example, the relative binding energy between the antimalarial drugs and receptor were: −7.0 to −10.0 kcal/mol, for DHFR-TS, −6.00 to −7.00 kcal/mol for pf-TERT, and −6.00 to −9.00 kcal/mol for pf-DHPS.
3.2. Binding Site Residues and Complexes
In this section, we explore the residues and binding pockets in the receptors to identify key features that favor binding. The binding pockets of the receptors are populated by polar and nonpolar residues positioned to form hydrogen bonds, electrostatic, Van Der Waals and hydrophobic interactions with the ligands (Figure5). The combination of residues, their orientation in the active site of pf-DHFR, pf-DHPS and pf-TERT receptors is critical for variation in affinity (Figure5 and Figure6). The binding residues for rest of complexes are summarized in supplemental materials (TableSM2, FigureSM1, Figures SM7-SM12).
The complex between h-DHFR and pyrimethamine, is stabilized by hydrophobic interactions involving residues F31, F34, I60, L67 and L22 (Figure5). The wt-pf-DHFR-pyrimethamine complex is stabilized by hydrogen bonding involving Y170 and I164, Van Der Waals interactions involving I112, F58, L119, L46 and G41 (FigureSM2). The dm-pf-DHFR-pyrimethamine, complex is stabilized by hydrogen bonding with S111 and G44, hydrophobic interactions with F58, L40 and L46 (FigureSM3). The qm-pf-DHFR-pyrimethamine complex, is stabilized by three hydrogen bonding interactions involving L40, A16, N54 and hydrophobic interactions with V195, L164 and L46 (FigureSM4). The increase in hydrogen bonding contacts in the interactions between the mutant forms of DHFR and pyrimethamine is an important finding which may indicate some unique features of pyrimethamine. The pyrimethamine complexes formed are stabilized by and hydrophobic interactions involving non-polar residues. The binding site residues for pf-DHFR-TS receptor include: Isoleucine 14(I14), Leucine 46(L46), Phenylalanine 58(F58), Serine 111(S111), Isoleucine 164(I164), Threonine 185(T185), Phenylalanine 58(F58) is identified as a key residue involved in affinity with Plasmodium enzyme receptors. This residue has been identified as a key factor in the affinity of novel antimalarial ligands like P218 with DHFR [21]. Our calculation also identifies F58 presence in the wt-pf-DHFR-pyrimethamine complex and dm-pf-DHFR-pyrimethamine complex (TableSM2) as critical for complex formation.
The active site of pf-DHFR, pf-DHPS and pf-TERT receptors contains a mixture of polar and non-polar residues. The difference in affinity observed between these receptors towards the antimalarial ligands could be explained by the size
Figure 5. (A) h-DHFR-pyrimethamine complex and (B) 2D interaction diagram; (C) qm-pf-DHFR-pyrimethamine complex and (D) 2D interaction diagram.
Figure 6. Surface created by amino acid residues around ligands. (A) pf-DHFR-TS-artemisinin complex; (B) pf-DHPS-artemisinin complex; (C) pf-TERT-artemisinin complex.
and nature (polar or nonpolar, acid or basic) of amino acids residues around the ligand (Figure6). The binding pocket of the pf-DHFR-TS-artemisinin complex (FigureSM5), comprise multiple hydrophobic residues, involved in strong Van Der Waals interactions with the cycloalkane rings or artemisinin. The snapshot of the molecular surface of the artemisinin described in (Figure6(A)) also points to a strong contribution from non-polar interactions. The polar S111 and T108 are responsible for the strong polar contribution to the stability of the complex (Figure6). The binding pocket of the pf-DHPS-artemisinin complex comprises multiple basic amino acids (Arg and Lys), acidic amino acids (Asp, Asn) and polar amino acids (Ser, Thr, Tyr, Asn, and Asp). The reduced affinity for interactions here could be a result of the diminished impact of Van Der Waals interaction with the large hydrophobic cyclic rings of artemisinin. The energy obtained could be explained by the strong presence of polar amino acids. The favorable energy and stability of the complex involving pf-TERT-artemisinin (FigureSM6) is due to both polar and non-polar contacts as we an abundance of nonpolar amino acids mixed with polar ones like D427 and T481.
4. Conclusion
In this study, the 3D structures of pf-TERT and pf-DHPS enzymes were obtained through homology modelling using chimera and modeller software. The interactions of 6 antimalarial drugs with pf-TERT, pf-DHPS models and DHFR enzymes were screened using molecular docking with AutoDock Vina. The docking calculations show that all the ligands involved formed stable complexes with the receptors. DHFR receptors, however, seem to have a high affinity compared to pf-TERT and pf-DHPS. Their simulations show that the benchmark drug artemisinin shows the strongest affinity with all the receptors and this can be attributed to the size, shape and bipolar nature of the ligand. This bipolar nature ensures one side is mostly hydrophobic and the other polar. We observe that the stable complexes are formed when the ligands successfully bind using both favorable polar and non-polar contacts. The performed study provides us with a pharmacophore that can be used for the development of new antimalarial chemotherapies.
Funding Details
This work was supported with funds from the school of science and technology, Georgia Gwinnett College.
Acknowledgements
Djogang thanks the School of Science and Technology, Chemistry, Georgia Gwinnett College (USA) for providing the software used in this work.
Abbreviations
DHFR: dihydrofolate reductase.
DHPS: dihydropteroate synthetase.
pf: Plasmodium falciparum.
pf-TERT: Plasmodium telomerase reverse-transcriptase.
Appendix A. Supplementary Data
Supplementary data to this article can be found online at https://doi.org/..............
Supplemental Material
Scheme 1. Folate pathway and antifolate drug against malarial.
Table SM1.Sequence similarity scores between target protein and identified template structures.
Table SM2.Residues involved different complexes and their energies
Figure SM1. (A) h-DHFR-artemisinin complex; (B) 2D interaction diagram.
Figure SM2. (A) wt-pf-DHFR-pyrimethamine complex; (B) 2D interaction diagram.
Figure SM3. (A) dm-pf-DHFR-pyrimethamine complex; (B) 2D interaction diagram.
Figure SM4. (A) qm-pf-DHFR-pyrimethamine complex; (B) 2D interaction diagram.
Figure SM5. (A) pf-DHPS-artemisinin complex; (B) 2D interaction diagram.
Figure SM6. (A) pf-TERT-artemisinin complex; (B) 2D interaction diagram.
Figure SM7. (A) wt-pf-DHFR-artemisinin complex; (B) 2D interaction diagram.
Figure SM8. (A) dm-pf-DHFR-artemisinin complex; (B) 2D interaction diagram.
Figure SM9. (A) qm-pf-DHFR-artemisinin complex; (B) 2D interaction diagram.
Figure SM10. Binding domain amino acids residues involved pf-DHFR-TS-artemisinin complex.
Figure SM11. Binding domain amino acids residues involved pf-DHPS-artemisinin complex.
Figure SM12. Binding domain amino acids residues involved pf-TERT-artemisinin complex.