The Application of MD Simulation to Lead Identification, Vaccine Design, and Structural Studies in Combat against Leishmaniasis - A Review

Drug discovery, vaccine design, and protein interaction studies are rapidly moving toward the routine use of molecular dynamics simulations (MDS) and related methods. As a result of MDS, it is possible to gain insights into the dynamics and function of identified drug targets, antibody-antigen interactions, potential vaccine candidates, intrinsically disordered proteins, and essential proteins. The MDS appears to be used in all possible ways in combating diseases such as cancer, however, it has not been well documented as to how effectively it is applied to infectious diseases such as Leishmaniasis. As a result, this review aims to survey the application of MDS in combating leishmaniasis. We have systematically collected articles that illustrate the implementation of MDS in drug discovery, vaccine development, and structural studies related to Leishmaniasis. Of all the articles reviewed, we identified that only a limited number of studies focused on the development of vaccines against Leishmaniasis through MDS. Also, the PCA and FEL studies were not carried out in most of the studies. These two were globally accepted utilities to understand the conformational changes and hence it is recommended that this analysis should be taken up in similar approaches in the future.

Global distribution of treatment resistance gene markers for leishmaniasis

BACKGROUND

Pentavalent antimonials (Sb(V)) such as meglumine antimoniate (Glucantime®) and sodium stibogluconate (Pentostam®) are used as first-line treatments for leishmaniasis, either alone or in combination with second-line drugs such as amphotericin B (Amp B), miltefosine (MIL), methotrexate (MTX), or cryotherapy. Therapeutic aspects of these drugs are now challenged because of clinical resistance worldwide.

METHODS

We reviewedthe recent original studies were assessed by searching in electronic databases such as Scopus, Pubmed, Embase, and Web of Science.

RESULTS

Studies on molecular biomarkers involved in drug resistance are essential for monitoring the disease. We reviewed genes and mechanisms of resistance to leishmaniasis, and the geographical distribution of these biomarkers in each country has also been thoroughly investigated.

CONCLUSION

Due to the emergence of resistant genes mainly in anthroponotic Leishmania species such as L. donovani and L. tropica, as the causative agents of ACL and AVL, respectively, selection of an appropriate treatment modality is essential. Physicians should be aware of the presence of such resistance for the selection of proper treatment modalities in endemic countries.

Evolutionary aspect of Miltefosine transporter proteins in Leishmania major

Transporter proteins, P-glycoprotein (P-gp) and P4ATPase-CDC50, are responsible for the transport of Miltefosine drug across cell membrane of a protozoan parasite Leishmania major. Mutations or change in activity of these proteins may lead to emergence of resistance in the parasite. Owing to the structural and functional importance of these transporter proteins, we have tried to decipher the evolutionary divergence of these Miltefosine transporter proteins across different forms of life including Protists, Fungi, Plants and Animals. We retrieved 96, 207, and 189 sequences of P-gp, P4ATPase and CDC50 proteins respectively, across diverse variety of organisms for the conserved analysis. Phylogenetic trees were constructed for these three transporter proteins based on Bayesian posterior probability inference. The evolutionary analysis concluded that these proteins remain highly conserved throughout the species diversity but still substantial differences in the proteins for host (Homo sapiens) and parasite (L. major) were observed which have led in targeting these Miltefosine transporter proteins in a parasite specific manner. The functional and structural components observed in terms of pattern resulting from the variability in the phylogenetic tree are outlined.

Computationally designed synthetic peptides for transporter proteins imparts allostericity in Miltefosine resistant L. major

The emergence of drug resistance is a major concern for combating against Cutaneous Leishmaniasis, a neglected tropical disease affecting 98 countries including India. Miltefosine is the only oral drug available for the disease and Miltefosine transporter proteins play a pivotal role in the emergence of drug-resistant Leishmania major. The cause of resistance is less accumulation of drug inside the parasite either by less uptake of the drug due to a decrease in the activity of P4ATPase-CDC50 complex or by increased efflux of the drug by P-glycoprotein (P-gp, an ABC transporter). In this paper, we are trying to allosterically modulate the behavior of resistant parasite (L. major) towards its sensitivity for the existing drug (Miltefosine, a phosphatidylcholine analog). We have used computational approaches to deal with the conservedness of the proteins and apparently its three-dimensional structure prediction through ab initio modeling. Long scale membrane-embedded molecular dynamics simulations were carried out to study the structural interaction and stability. Parasite-specific motifs of these proteins were identified based on the machine learning technique, against which a peptide library was designed. The protein-peptide docking shows good binding energy of peptides Pg5F, Pg8F and PC2 with specific binding to the motifs. These peptides were tested both in vitro and in vivo, where Pg5F in combination with PC2 showed 50-60% inhibition in resistant L. major's promastigote and amastigote forms and 80-90% decrease in parasite load in mice. We posit a model system wherein the data provide sufficient impetus for being novel therapeutics in order to counteract the drug resistance phenotype in Leishmania parasites.

Comparative transcriptomic analysis of antimony resistant and susceptible Leishmania infantum lines

Background

One of the major challenges to leishmaniasis treatment is the emergence of parasites resistant to antimony. To study differentially expressed genes associated with drug resistance, we performed a comparative transcriptomic analysis between wild-type and potassium antimonyl tartrate (Sb^III)-resistant Leishmania infantum lines using high-throughput RNA sequencing.

Methods

All the cDNA libraries were constructed from promastigote forms of each line, sequenced and analyzed using STAR for mapping the reads against the reference genome (L. infantum JPCM5) and DESeq2 for differential expression statistical analyses. All the genes were functionally annotated using sequence similarity search.

Results

The analytical pipeline considering an adjusted p-value < 0.05 and fold change > 2.0 identified 933 transcripts differentially expressed (DE) between wild-type and Sb^III-resistant L. infantum lines. Out of 933 DE transcripts, 504 presented functional annotation and 429 were assigned as hypothetical proteins. A total of 837 transcripts were upregulated and 96 were downregulated in the Sb^III-resistant L. infantum line. Using this DE dataset, the proteins were further grouped in functional classes according to the gene ontology database. The functional enrichment analysis for biological processes showed that the upregulated transcripts in the Sb^III-resistant line are associated with protein phosphorylation, microtubule-based movement, ubiquitination, host–parasite interaction, cellular process and other categories. The downregulated transcripts in the Sb^III-resistant line are assigned in the GO categories: ribonucleoprotein complex, ribosome biogenesis, rRNA processing, nucleosome assembly and translation.

Conclusions

The transcriptomic profile of L. infantum showed a robust set of genes from different metabolic pathways associated with the antimony resistance phenotype in this parasite. Our results address the complex and multifactorial antimony resistance mechanisms in Leishmania, identifying several candidate genes that may be further evaluated as molecular targets for chemotherapy of leishmaniasis.

Biophysical Approaches Facilitate Computational Drug Discovery for ATP-Binding Cassette Proteins

Although membrane proteins represent most therapeutically relevant drug targets, the availability of atomic resolution structures for this class of proteins has been limited. Structural characterization has been hampered by the biophysical nature of these polytopic transporters, receptors, and channels, and recent innovations to in vitro techniques aim to mitigate these challenges. One such class of membrane proteins, the ATP-binding cassette (ABC) superfamily, are broadly expressed throughout the human body, required for normal physiology and disease-causing when mutated, yet lacks sufficient structural representation in the Protein Data Bank. However, recent improvements to biophysical techniques (e.g., cryo-electron microscopy) have allowed for previously “hard-to-study” ABC proteins to be characterized at high resolution, providing insight into molecular mechanisms-of-action as well as revealing novel druggable sites for therapy design. These new advances provide ample opportunity for computational methods (e.g., virtual screening, molecular dynamics simulations, and structure-based drug design) to catalyze the discovery of novel small molecule therapeutics that can be easily translated from computer to bench and subsequently to the patient's bedside. In this review, we explore the utility of recent advances in biophysical methods coupled with well-established in silico techniques towards drug development for diseases caused by dysfunctional ABC proteins.

Computational modelling of efflux pumps and their inhibitors

Antimicrobial resistance is based on the multifarious strategies that bacteria adopt to face antibiotic therapies, making it a key public health concern of our era. Among these strategies, efflux pumps (EPs) contribute significantly to increase the levels and profiles of resistance by expelling a broad range of unrelated compounds - buying time for the organisms to develop specific resistance. In Gram-negative bacteria, many of these chromosomally encoded transporters form multicomponent 'pumps' that span both inner and outer membranes and are driven energetically by a primary or secondary transporter component.One of the strategies to reinvigorate the efficacy of antimicrobials is by joint administration with EP inhibitors (EPI), which either block the substrate binding and/or hinder any of the transport-dependent steps of the pump. In this review, we provide an overview of multidrug-resistance EPs, their inhibition strategies and the relevant findings from the various computational simulation studies reported to date with respect to deciphering the mechanism of action of inhibitors with the purpose of improving their rational design.

Specific mutations in mammalian P4-ATPase ATP8A2 catalytic subunit entail differential glycosylation of the accessory CDC50A subunit

P4-ATPases, or flippases, translocate phospholipids between the two leaflets of eukaryotic biological membranes. They are essential to the physiologically crucial phospholipid asymmetry and involved in severe diseases, but their molecular structure and mechanism are still unresolved. Here, we show that in an extensive mutational alanine screening of the mammalian flippase ATP8A2 catalytic subunit, five mutations stand out by leading to reduced glycosylation of the accessory subunit CDC50A. These mutations may disturb the interaction between the subunits.