Crystal structure of the ternary complex of Leishmania major pteridine reductase 1 with the cofactor NADP+/NADPH and the substrate folic acid

Experimental structure based drug design (SBDD) applications for anti-leishmanial drugs: A paradigm shift?

Leishmaniasis is a group of neglected tropical diseases caused by at least 20 species of Leishmania protozoa, which are spread by the bite of infected sandflies. There are three main forms of the disease: cutaneous leishmaniasis (CL, the most common), visceral leishmaniasis (VL, also known as kala-azar, the most serious), and mucocutaneous leishmaniasis. One billion people live in areas endemic to leishmaniasis, with an annual estimation of 30,000 new cases of VL and more than 1 million of CL. New treatments for leishmaniasis are an urgent need, as the existing ones are inefficient, toxic, and/or expensive. We have revised the experimental structure-based drug design (SBDD) efforts applied to the discovery of new drugs against leishmaniasis. We have grouped the explored targets according to the metabolic pathways they belong to, and the key achieved advances are highlighted and evaluated. In most cases, SBDD studies follow high-throughput screening campaigns and are secondary to pharmacokinetic optimization, due to the majoritarian belief that there are few validated targets for SBDD in leishmaniasis. However, some SBDD strategies have significantly contributed to new drug candidates against leishmaniasis and a bigger number holds promise for future development.

Binding profile of quinonoid-dihydrobiopterin to quinonoid-dihydropteridine reductase examined by in silico and in vitro analyses

Quinonoid dihydropteridine reductase (QDPR) catalyses the reduction of quinonoid-form dihydrobiopterin (qBH2) to tetrahydrobiopterin (BH4). BH4 metabolism is a drug target for neglected tropical disorders because trypanosomatid protozoans, including Leishmania and Trypanosoma, require exogenous sources of biopterin for growth. Although QDPR is a key enzyme for maintaining intracellular BH4 levels, the precise catalytic properties and reaction mechanisms of QDPR are poorly understood due to the instability of quinonoid-form substrates. In this study, we analysed the binding profile of qBH2 to human QDPR in combination with in silico and in vitro methods. First, we performed docking simulation of qBH2 to QDPR to obtain possible binding modes of qBH2 at the active site of QDPR. Then, among them, we determined the most plausible binding mode using molecular dynamics simulations revealing its atomic-level interactions and confirmed it with the in vitro assay of mutant enzymes. Moreover, it was found that not only qBH2 but also quinonoid-form dihydrofolate (qDHF) could be potential physiological substrates for QDPR, suggesting that QDPR may be a bifunctional enzyme. These findings in this study provide important insights into biopterin and folate metabolism and would be useful for developing drugs for neglected tropical diseases.

Structure and dynamics of pteridine reductase 1: the key phenomena relevant to enzyme function and drug design

Pteridine reductase 1 (PTR1) is a folate and pterin pathway enzyme unique for pathogenic trypanosomatids. As a validated drug target, PTR1 has been the focus of recent research efforts aimed at finding more effective treatments against human parasitic diseases such as leishmaniasis or sleeping sickness. Previous PTR1-centered structural studies highlighted the enzyme characteristics, such as flexible regions around the active site, highly conserved structural waters, and species-specific differences in pocket properties and dynamics, which likely impacts the binding of natural substrates and inhibitors. Furthermore, several aspects of the PTR1 function, such as the substrate inhibition phenomenon and the level of ligand binding cooperativity in the enzyme homotetramer, likely related to the global enzyme dynamics, are poorly known at the molecular level. We postulate that future drug design efforts could greatly benefit from a better understanding of these phenomena through studying both the local and global PTR1 dynamics. This review highlights the key aspects of the PTR1 structure and dynamics relevant to structure-based drug design that could be effectively investigated by modeling approaches. Particular emphasis is given to the perspective of molecular dynamics, what has been accomplished in this area to date, and how modeling could impact the PTR1-targeted drug design in the future.

Tetrahydrobenzo[h]quinoline derivatives as a novel chemotype for dual antileishmanial-antimalarial activity graced with antitubercular activity: Design, synthesis and biological evaluation

Derivatives with tetrahydrobenzo[h]quinoline chemotype were synthesized via one-pot reactions and evaluated for their antileishmanial, antimalarial and antitubercular activities. Based on a structure-guided approach, they were designed to possess antileishmanial activity through antifolate mechanism, via targeting Leishmania major pteridine reductase 1 (Lm-PTR1). The in vitro antipromastigote and antiamastigote activity are promising for all candidates and superior to the reference miltefosine, in a low or sub micromolar range of activity. Their antifolate mechanism was confirmed via the ability of folic and folinic acids to reverse the antileishmanial activity of these compounds, comparably to Lm-PTR1 inhibitor trimethoprim. Molecular dynamics simulations confirmed a stable and high potential binding of the most active candidates against leishmanial PTR1. For the antimalarial activity, most of the compounds exhibited promising antiplasmodial effect against P. berghei with suppression percentage of up to 97.78%. The most active compounds were further screened in vitro against the chloroquine resistant strain P. falciparum, (RKL9) and showed IC50 value range of 0.0198-0.096 μM, compared to IC50 value of 0.19420 μM for chloroquine sulphate. Molecular docking of the most active compounds against the wild-type and quadruple mutant pf DHFR-TS structures rationalized the in vitro antimalarial activity. Some candidates showed good antitubercular activity against sensitive Mycobacterium tuberculosis in a low micromolar range of MIC, compared to 0.875 μM of isoniazid. The top active ones were further tested against a multidrug-resistant strain (MDR) and extensively drug-resistant strain (XDR) of Mycobacterium tuberculosis. Interestingly, the in vitro cytotoxicity test of the best candidates displayed high selectivity indices emphasizing their safety on mammalian cells. Generally, this work introduces a fruitful matrix for new dual acting antileishmanial-antimalarial chemotype graced with antitubercular activity. This would help in tackling drug-resistance issues in treating some Neglected Tropical Diseases.

Delving in folate metabolism in the parasite Leishmania major through a chemogenomic screen and methotrexate selection

Most of our understanding of folate metabolism in the parasite Leishmania is derived from studies of resistance to the antifolate methotrexate (MTX). A chemical mutagenesis screen of L. major Friedlin and selection for resistance to MTX led to twenty mutants with a 2- to 400-fold decrease in MTX susceptibility in comparison to wild-type cells. The genome sequence of the twenty mutants highlighted recurrent mutations (SNPs, gene deletion) in genes known to be involved in folate metabolism but also in novel genes. The most frequent events occurred at the level of the locus coding for the folate transporter FT1 and included gene deletion and gene conversion events, as well as single nucleotide changes. The role of some of these FT1 point mutations in MTX resistance was validated by gene editing. The gene DHFR-TS coding for the dihydrofolate reductase-thymidylate synthase was the second locus with the most mutations and gene editing confirmed a role in resistance for some of these. The pteridine reductase gene PTR1 was mutated in two mutants. The episomal overexpression of the mutated versions of this gene, but also of DHFR-TS, led to parasites several fold more resistant to MTX than those overexpressing the wild-type versions. Genes with no known link with folate metabolism and coding for a L-galactolactone oxidase or for a methyltransferase were mutated in specific mutants. Overexpression of the wild-type versions of these genes in the appropriate mutants reverted their resistance. Our Mut-seq approach provided a holistic view and a long list of candidate genes potentially involved in folate and antifolate metabolism in Leishmania.

New antileishmanial quinoline linked isatin derivatives targeting DHFR-TS and PTR1: Design, synthesis, and molecular modeling studies

In a search for new drug candidates for one of the neglected tropical diseases, leishmaniasis, twenty quinoline-isatin hybrids were synthesized and tested for their in vitro antileishmanial activity against Leishmaniamajor strain. All the synthesized compounds showed promising in vitro activity against the promastigote form in a low micromolar range (IC₅₀ = 0.5084-5.9486 μM) superior to the reference miltefosine (IC₅₀ = 7.8976 μM). All the target compounds were then tested against the intracellular amastigote form and showed promising inhibition effects (IC₅₀ = 0.60442-8.2948 μM versus 8.08 μM for miltefosine). Compounds 4e, 4b and 4f were shown to possess the highest antileishmanial activity against both promastigote and amastigote forms. The most active compounds were proven to exhibit their significant antileishmanial effects through antifolate mechanism, targeting DHFR-TS and PTR1. To evaluate the safety profile of the most active derivatives 4e, 4b and 4f, the in vitro cytotoxicity test was carried out and displayed higher selectivity indices than the reference miltefosine. Molecular docking within putative target protein PTR1 confirmed the high potentiality of the most active compounds 4e, 4b and 4f to block the catalytic activity of Lm-PTR1.

Selene-Ethylenelacticamides and N-Aryl-Propanamides as Broad-Spectrum Leishmanicidal Agents

The World Health Organization classifies Leishmania as one of the 17 “neglected diseases” that burden tropical and sub-tropical climate regions with over half a million diagnosed cases each year. Despite this, currently available anti-leishmania drugs have high toxicity and the potential to be made obsolete by parasite drug resistance. We chose to analyze organoselenides for leishmanicidal potential given the reduced toxicity inherent to selenium and the displayed biological activity of organoselenides against Leishmania. Thus, the biological activities of 77 selenoesters and their N-aryl-propanamide derivatives were predicted using robust in silico models of Leishmania infantum, Leishmania amazonensis, Leishmania major, and Leishmania (Viannia) braziliensis. The models identified 28 compounds with >60% probability of demonstrating leishmanicidal activity against L. infantum, and likewise, 26 for L. amazonesis, 25 for L. braziliensis, and 23 for L. major. The in silico prediction of ADMET properties suggests high rates of oral absorption and good bioavailability for these compounds. In the in silico toxicity evaluation, only seven compounds showed signs of toxicity in up to one or two parameters. The methodology was corroborated with the ensuing experimental validation, which evaluated the inhibition of the Promastigote form of the Leishmania species under study. The activity of the molecules was determined by the IC50 value (µM); IC50 values < 20 µM indicated better inhibition profiles. Sixteen compounds were synthesized and tested for their activity. Eight molecules presented IC50 values < 20 µM for at least one of the Leishmania species under study, with compound NC34 presenting the strongest parasite inhibition profile. Furthermore, the methodology used was effective, as many of the compounds with the highest probability of activity were confirmed by the in vitro tests performed.