Next Article in Journal
Mutation S115T in IMP-Type Metallo-β-Lactamases Compensates for Decreased Expression Levels Caused by Mutation S119G
Next Article in Special Issue
Comparative Characterization of Plasmodium falciparum Hsp70-1 Relative to E. coli DnaK Reveals the Functional Specificity of the Parasite Chaperone
Previous Article in Journal
The Potential Involvement of an ATP-Dependent Potassium Channel-Opening Mechanism in the Smooth Muscle Relaxant Properties of Tamarix dioica Roxb.
Previous Article in Special Issue
Zika and Flavivirus Shell Disorder: Virulence and Fetal Morbidity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Novel Serine/Threonine Protein Kinase LmjF.22.0810 from Leishmania major May Be Involved in the Resistance to Drugs such as Paromomycin

by
Andrés Vacas
1,
Celia Fernández-Rubio
1,
Miriam Algarabel
1,
José Peña-Guerrero
1,
Esther Larrea
1,
Fabio Rocha Formiga
2,
Alfonso T. García-Sosa
3 and
Paul A. Nguewa
1,*
1
Department of Microbiology and Parasitology, ISTUN Institute of Tropical Health, University of Navarra, IdiSNA (Navarra Institute for Health Research), Pamplona, E-31008 Navarra, Spain
2
Aggeu Magalhães Institute, Oswaldo Cruz Foundation (FIOCRUZ). Graduate Program in Applied Cellular and Molecular Biology (BCMA), University of Pernambuco (UPE), Recife/PE 50.670-420, Brazil
3
Institute of Chemistry, University of Tartu, Ravila 14a, 54011 Tartu, Estonia
*
Author to whom correspondence should be addressed.
Biomolecules 2019, 9(11), 723; https://doi.org/10.3390/biom9110723
Submission received: 6 October 2019 / Revised: 1 November 2019 / Accepted: 5 November 2019 / Published: 11 November 2019

Abstract

:
The identification and clarification of the mechanisms of action of drugs used against leishmaniasis may improve their administration regimens and prevent the development of resistant strains. Herein, for the first time, we describe the structure of the putatively essential Ser/Thr kinase LmjF.22.0810 from Leishmania major. Molecular dynamics simulations were performed to assess the stability of the kinase model. The analysis of its sequence and structure revealed two druggable sites on the protein. Furthermore, in silico docking of small molecules showed that aminoglycosides preferentially bind to the phosphorylation site of the protein. Given that transgenic LmjF.22.0810-overexpressing parasites displayed less sensitivity to aminoglycosides such as paromomycin, our predicted models support the idea that the mechanism of drug resistance observed in those transgenic parasites is the tight binding of such compounds to LmjF.22.0810 associated with its overexpression. These results may be helpful to understand the complex machinery of drug response in Leishmania.

1. Introduction

In 2006, paromomycin (PMM) was introduced for the treatment of visceral leishmaniasis. Its effectiveness and low cost allowed its use as a first-line alternative in the setting of resistance to traditional antileishmanial treatments [1]. PMM belongs to the aminoglycoside (AG) antibiotic family that has been broadly used for treating Gram-negative bacterial infections [2,3]. However, the emergent resistance to these antibiotics, conferred principally by bacterial eukaryotic-like protein kinases (ePKs), has left AGs ineffective against several pathogenic strains [2,3,4]. Concerning PMM alternatives against leishmaniasis, the emergence of drug-resistant strains has already rendered antimony useless against Leishmania infections; furthermore, other second-line drugs such as pentamidine have reported failure rates of around 30% in India [5]. On the other hand, the cost concerns of miltefosine or lipid formulations of amphotericin B remain a considerable challenge for the control and reduction of the burden of leishmaniasis [5,6]. Therefore, PMM stands as a convincing treatment for a broader usage [1]. However, even after its approval and introduction into the clinic, its mechanisms of action still need to be further studied to understand the development of resistance in Leishmania [7].
Trypanosomatid signaling pathways are not likely to culminate at the level of transcription since Trypanosoma and Leishmania transcribe their genes as polycistronic units with little evidence of regulation [8]. Therefore, as a consequence of inference with other eukaryotic cells, it is expected that protein phosphorylation may be critical for parasite’s development [9]. In fact, changes in the phosphorylation pattern and expression of ePKs and phosphatases in Leishmania and other kinetoplastids during their life cycles are common phenomena and have been demonstrated [10,11,12,13]. Protein kinases (PKs) are key regulators of cellular signaling and have long been recognized as important drug targets in a variety of diseases (such as cancer and inflammatory and infectious diseases) [14,15,16,17]. The human kinome (518 kinases, ~2% of the human genome) is the largest “druggable” family of molecular targets in humans [18]. Kinases also represent a relatively large group of the predicted protein-coding genes from trypanosomatid parasites (approximately 2% of each genome). Protein essentiality, druggability, and the absence of orthologs in the human proteome are some of the key features that should be displayed by a potential therapeutic target [19]. Interestingly, PKs play vital roles in Leishmania and Trypanosoma such as the control of virulence and growth of parasites, allowing their consideration as potential therapeutic targets against neglected tropical diseases (NTDs) [20,21,22,23].
Several dozens of small-molecule kinase inhibitors have been either approved or studied for the treatment of various human diseases which include cancer, Alzheimer’s disease, as well as infections [20,24,25,26,27]. Nevertheless, due to the phylogenetic proximity with their hosts, targeting ePKs from eukaryotic microorganisms is more challenging than dealing with bacterial PKs [6]. Analyses conducted by Florez et al. (2010) determined phosphoproteins as the largest cluster among the 384 potential drug targets discovered in Leishmania major [28]. Amongst these proteins, 91 kinases were predicted to be essential for parasite’s survival and also shown to lack orthologous proteins in the human kinome. Thus, focusing on peculiarities found in PKs from eukaryotic microorganisms is essential to avoid off-target activity and cross-reactivity of inhibitors with kinases from their hosts [6,18].
Fang et al. (2013) noticed that understanding the structural details of kinase conformational changes eventually allows the rational design of innovative kinase inhibitors [29]. Moreover, the structure of the enzyme-bound antagonist complex is used in the classification of kinase inhibitors [30]. Most of the FDA-approved kinase inhibitors (type I) reversibly bind to the adenine site of the ATP pocket from kinases in their active state [14], therefore preventing the co-factor ATP from binding to the site. However, they have to compete with ATP towards binding to this region [18]. The structure of active kinases is highly conserved in all organisms; on the contrary, kinases in their inactive state are structurally more diverse and dynamic [31]. Type II and III inhibitors target kinases during their inactive state, whereas type IV inhibitors target regions outside of the ATP cleft and the phosphoacceptor site in active and inactive kinases [30,31,32]. Similarly, in recent years, increased attention has been placed on small molecules that target the protein substrate binding site [18]. Nevertheless, it is well known that designing molecules for such relatively flat surfaces is tricky. Consequently, such a challenging strategy has not yet been extensively studied [32].
Since calcium ion (Ca2+) is an important second messenger in the biology of trypanosomatid parasites and is implicated in their differentiation, host cell invasion, and flagellum motility [33], our study was focused on the characterization of a protein from the Ca2+/calmodulin-dependent group of kinases (CAMK) in Leishmania. This novel molecule has also been predicted to be a therapeutic target [28]. Our work represents the first report describing this gene (LmjF.22.0810). Considering the lack of orthologous sequences, we decided to name this gene “Jean3”. Moreover, due to the unavailability of an X-ray 3D structure of the protein LmjF.22.0810, we were prompted to develop for the first time a homology model for LmjF.22.0810. Various methods were then employed to validate the predicted structure before it was used to study in silico its ability to bind different pharmacological substrates, including several AGs. After cloning and overexpressing LmjF.22.0810 in L. major, we found that such parasites were less sensitive to PMM and geneticin (G418). These outcomes, together with the comparative analysis of the primary sequence of the protein, allow us to propose a novel mechanism of AG resistance in Leishmania mediated by the tight binding of AGs to the protein substrate site of LmjF.22.0810. Our results will be valuable not only to understand the complexity of drug resistance in Leishmania but also to design novel inhibitors against homologues of LmjF.22.0810 kinases.

2. Materials and Methods

2.1. Database Inspection and Nucleotide and Protein Sequence Alignment

Nucleotide and genomic sequence data, as well as Jean3 orthologous protein sequences, were obtained from TriTrypDB [34]. Multiple sequence alignments were prepared using the MUSCLE algorithm [35] implemented in Geneious v9.1.7 [36] with a maximum number of iterations of 8.

2.2. Evaluation of Secondary Structure, Conserved Domains, and Post-Translational Modifications

The secondary structure of LmJean3 was predicted using the EMBOSS v6.5.7 tool [37] implemented in Geneious v9.1.7. Conserved domain identification was performed using the Conserved Domains Database (CDD) [38,39]. LmJean3 disordered regions were inspected by the Meta-Disorder method from the PredictProtein server [40]. The kinase active site (Prosite pattern PS00108: [LIVMFYC]-x-[HY]-x-D-[LIVMFY]-K-x(2)-N-[LIVMFYCT](3)) and ATP-binding site (Prosite pattern PS00107: [LIV]-G-{P}min-G-{P}-[FYWMGSTNH]-[SGA]-{PW}-[LIVCAT]-{PD}-x-[GSTACLIVMFY]-x(5,18)-[LIVMFYWCSTAR]-[AIVP]-[LIVMFAGCKR]-K) were determined by InterPro: Protein sequence analysis and classification [41]. Kinase active residues were predicted using the Pfam database from the European Bioinformatics Institute [42] (Uniprot accession number: Q4QBR6). EF-hand motifs were scanned using the PATTINPROT tool from the PRABI-Lyon-Gerland Network Protein Sequence analysis (NPS@) Web server (https://npsa-prabi.ibcp.fr). Prediction of EF-hand motifs was performed by searching for the motif signature of the canonical EF loop with both flanking helices described by Zhou et al., (2006) (x-{DNQ}-x(2)-{GP}-{ENPQS}-x(2)-{DPQR}-[DNS]-x-[DNS]-{FLIVWY}-[DNESTG]-[DNQGHRK]-{GP}-[LIVMC]-[DENQSTAGC]-x(2)-[ED]-[FLYMVIW]-x(2)-{NPS}-{DENQ}-x(3)) [43].
The following tools were used to assess post-translational modifications: The NetNGlyc 1.0 server from the Technical University of Denmark was used for N-terminal glycosylation prediction (http://www.cbs.dtu.dk/services/NetNGlyc/). The MYR Predictor of the IMP Bioinformatics Group assessed Jean3 sequences for possible myristoylation patterns (http://mendel.imp.ac.at/myristate/SUPLpredictor.htm). The TMHMM Server v2.0 from the Technical University of Denmark was used to determine protein transmembrane regions (http://www.cbs.dtu.dk/services/TMHMM/). PredGPI [44] was used to investigate the presence of GPI anchor sites. ProtParam [45] was used to predict the molecular mass of the protein and its theoretical isoelectric point. Serine and threonine phosphorylation sites were assessed by neural network prediction using the NetPhos 2.0 server [46]. Prediction of the presence and location of signal peptide cleavage sites was performed using the SignalP 4.1 server [47].
Genomic sequence distribution, multiple alignments, and plasmid diagrams were computed and arranged using Geneious version 9.1.7 created by Biomatters. Schematic representation of the protein with its motifs, regions, and conserved domains was drawn using Inkscape software v0.91 [48].

2.3. Phylogenetic Analysis

The phylogenetic tree was constructed with the sequences retrieved from the NCBI protein database and TriTrypDB [34]. LmJean3 (XP_001683232) full deduced amino acid sequence and individual conserved domains were used to perform BLASTp searches within the non-redundant protein sequences database (nr), with the default search parameters. Sequences were then aligned using the MAFFT algorithm [49] implemented in Geneious software, with a JTT200 scoring matrix.
ProtTest 2.4 server [50] was used to determine the best-fit model for the phylogenetic reconstruction of Jean3. According to ProtTest, the best model for protein evolution of the 23 sequences included in the analysis, with the smallest Akaike Information Criterion (AIC), was Le and Gascuel (LG) [51] with the following parameters: gamma shape (4 rate categories) of 1.295, proportion of invariable sites of 0.035; the observed amino acid frequencies obtained from the alignment were: A: 0.066, C: 0.016, D: 0.053, E: 0.069, F: 0.037, G: 0.066, H: 0.031, I: 0.060, K: 0.059, L: 0.100, M: 0.023, N: 0.045, P: 0.048, Q: 0.030, R: 0.061, S: 0.083, T: 0.047, V: 0.068, W: 0.007, and Y: 0.030. Geneious was used to portray the Newick output from the best tree (with optimized topology and branch lengths).

2.4. Screening for Homology Modeling Templates

The 262-amino acid sequence corresponding to the Ser/Thr domain of LmJean3 was used to perform a BLAST search within the Swiss-Model Repository [52]. Sequence alignment of the kinase domains from CIPK23, CIPK24/SOS2, Snf1, LdBPK_220630.1.1 (LdJean3), LbrM.22.0750 (LbJean3), and LmjF.22.0810 (LmJean3) was performed using the ClustalW iterative algorithm [53].

2.5. Molecular Modeling and 3D-Structure Validation

For homology modeling, the prediction server Robetta [54] was employed. Robetta uses the ROSETTA software to predict structures under the Ginzu protocol [54] either by homology modeling or by ab initio modeling. The C-terminal 97 amino acids (disordered region) from LmJean3 sequence were removed to improve the template selection by the software and consequently the final model. Ginzu alignment from the resultant N-terminal 275 amino acids from LmJean3 showed that the best candidate (0.8112 confidence) for comparative modeling was the B-chain from the 4CZT entry from the Protein Data Bank (PDB) [55]. 4CZT corresponds to the crystal structure of the kinase domain of CIPK23, determined by Chaves-Sanjuan et al. (2014) using X-ray diffraction with 2.3 Å resolution [56].
The backbone Cα root-mean-square deviation (RMSD) between 4CZT and the LmJean3 structures obtained from Robetta was determined using PyMOL [57]. The model with the lowest RMSD value was chosen for energy minimization using the YASARA force field [58]. The final model solvation accessibility was assessed by SolvX (http://ekhidna.biocenter.helsinki.fi/solvx/start). The three-dimensional profile of LmJean3 was determined by Verify3D [59]. The absolute quality of the model was estimated with QMEAN6 [60,61,62,63], and the atomic empirical mean force potential was determined by atomic non-local environment assessment (ANOLEA) [64,65,66]. The Swiss-Model Protein Structure and Assessment tools server [52] was used to assess the stereochemical quality of the model. The correctness of residue prediction was measured using PROCHECK [67] with a minimum resolution of 2.5 Å; the empirical force field energy was evaluated by GROMOS [68]. DSSP was used to compare the predicted model with the structural motifs from the secondary structure [69,70]. The three-dimensional model and the template were also checked using PROSA [71,72]. Finally, the structure quality of the template and the predicted LmJean3 structure was estimated using ERRAT [73]. All ribbon figures were designed using PyMOL v1.8.4.0 (https://sourceforge.net/projects/pymol/files/pymol/1.8/).

2.6. MD Simulation

The OPLS-AA force field [74] for the protein and the SPC water model [75] were employed. The predicted structure of LmJean3 was solvated in a rhombic dodecahedron box of water. The distance between the structure and the edge of the box was set to 12 Å, resulting in a minimal distance between periodic images of at least 24 Å. The net charge of the protein was neutralized by adding three positive counter-ions. Periodic boundaries were implemented in all three dimensions. Energy minimization of 1373 steps using the steepest descent algorithm was followed by two equilibration phases for 200 ps (first NVT—constant system, volume and temperature and then NPT—constant system, pressure and temperature ensemble). The temperature and pressure were controlled by a V-rescale thermostat and a Parrinello-Rahman barostat, respectively [76,77]. The equilibration phase was followed by 20 ns MD simulations, carried out using GROMACS v5.1.1 [78]. The integration time-step was of 2 fs, and the atom positions of the solute were written to file every 10 ps (every 5000 steps). To check the stability of the simulations, the RMSDs of the Cα atoms with respect to the minimized starting structure were calculated and monitored.
Data evaluation from MD simulations was carried out using the software package VMD v1.9.3 [79]. Salt bridges were identified by the following criteria: distance between any of the oxygen atoms of acidic residues and the nitrogen atoms of basic residues found within 3.2 Å, and centroids of the sidechain charged groups found within 4.0 Å of each other. Then, the trajectories were analyzed, and the presence of the identified salt bridges was assessed in each frame of the simulation. Similarly, hydrogen bonds were calculated with a cutoff distance of 3.2 Å and a cutoff angle of 20 degrees. A list was generated with those hydrogen bonds found in at least 30% of the trajectory (Supplemental Table S1). Residue interactions previously identified as salt bridges were not considered for hydrogen bond evaluation. In addition, the average coordinates over the trajectory were calculated and used to write a PDB file using VMD. Then, the RMSDs of the structure of the initial frame and the average structure were calculated using PyMOL v1.8.4.0.

2.7. Protein Preparation, Binding Site Identification, and High-Throughput in Silico Docking

Binding site identification, protein preparation, ligand preparation, and Glide docking were performed using Schrödinger Suite (Release 2016-1). Using Schrödinger’s protein preparation wizard [80], hydrogen atoms were added after deleting any original ones, followed by adjustment of bond orders for amino acid residues and ligand. The protonation and tautomeric states of Asp, Glu, Arg, Lys, and His were adjusted to match the physiological pH of 7.4. Possible orientations of Asn and Gln residues were generated. Hydrogen bond optimization (ProtAssign) sampling was performed using PROPKA [81] at a pH of 7.4.
Binding sites were identified using Schrödinger SiteMap v3.8 [82,83]. The region corresponding to the peptide substrate binding site was selected for grid files generation with Glide v7.0 [84,85]. The generated grid was used for in silico docking of small molecular compounds from the Zim dataset from ZINC [86] using Schrödinger’s Virtual Screening Workflow (VSW) with Epik penalties turned on [87,88]. From the results obtained, the molecular structures with the highest predicted affinity and lowest docking scores to the active site were selected (within Glide’s reported error of 2 kcal/mol).

2.8. Parasite Culture Conditions

L. major (Lv39c5) promastigotes were grown at 26 °C in M199 medium (Sigma-Aldrich, St. Louis, MO, USA) supplemented with 25 mM HEPES (pH 7.2; Sigma-Aldrich), 0.1 mM adenine (Sigma-Aldrich), 0.0005% (w/v) hemin (Sigma-Aldrich), 2 mg/mL biopterin (Sigma-Aldrich), 0.0001% (w/v) biotin (Sigma-Aldrich), 10% (v/v) heat-inactivated fetal bovine serum (Gibco Laboratories, Grand Island, NY, USA), and an antibiotic cocktail (50 U/mL penicillin, 50 mg/mL streptomycin) (Sigma-Aldrich). L. major cultures used for qPCR analysis were grown in Schneider’s medium (Gibco Laboratories) supplemented with 10% (v/v) heat-inactivated fetal bovine serum (Gibco Laboratories) and 40 μg/mL gentamicin (Sigma) at 26 °C.

2.9. Genetic Manipulation of L. major

The plasmid pXG-mCherry12 was constructed by our group as previously reported [89]. L. major DNA was extracted following the protocol previously described by Medina-Acosta et al., (1993) [90]. Then, the coding DNA sequence (CDS) corresponding to LmJean3 was amplified by PCR from L. major genomic DNA using the primers J3XStop-Fw and J3XStop-Rv for the construct pXG-LmJean3-mCherry12 (Supplemental Figure S1A), J3NotI-Fw and J3NotI-Rv for pXG-GFP2+-LmJean3 (Supplemental Figure S1B), and J3SF and J3SR for pXG-Hyg-LmJean3 (Supplemental Figure S1C). Primer sequences are detailed in Supplemental Table S2. Each of LmJean3 PCR products was then ligated into pCR®2.1-TOPO® (ThermoFisher Scientific, Rockville, MD, USA) cloning vectors following the manufacturer’s protocol. The ligation products were used to transform DH5α Escherichia coli bacteria by the heat shock method. Plasmids were then extracted from the bacteria and digested with the following restriction enzymes: BstXI (Clonetech, Palo Alto, CA, USA) for the pXG-LmJean3-mCherry12 construct (Supplemental Figure S1A), NotI-HF (New England Biolabs, Massachusetts, USA) for the pXG-GFP2+-LmJean3 construct (Supplemental Figure S1B), or SmaI (Clonetech) for the pXG-Hyg-LmJean3 plasmid (Supplemental Figure S1C). Digestions were then gel-purified, and LmJean3 sequences were ligated into the cloning cage of pXG-mCherry12, pXG-GFP2+ or pXG-Hyg plasmids (Supplemental Figure S1). LmJean3 sequence and orientation within the constructed vectors, pXG-LmJean3-mCherry12, pXG-GFP2+-LmJean3, and pXG-Hyg-LmJean3, were assessed by DNA sequencing and PCR, respectively.
A total number of 108 log-phase L. major parasites were used to transfect the plasmids pXG-LmJean3-mCherry12, pXG-GFP2+-LmJean3, and pXG-Hyg-LmJean3 by electroporation. A BioRad Gene Pulser II machine (Bio-Rad Laboratories, Hercules, CA, USA) was used for electroporation, following a method previously described by Cruz et al., (1991) [91]. Recombinant colonies were isolated from M199 agar plates supplemented with 100 μg/mL of hygromycin B Gold (InvivoGen Europe, Toulouse, France) for pXG-LmJean3-mCherry12 and pXG-Hyg-LmJean3, or 30 μg/mL geneticin (Sigma Aldrich) for pXG-GFP2+-LmJean3.

2.10. Fluorescence Microscopy

A final number of 2 × 107 log-phase pXG-LmJean3-mCherry12 and pXG-GFP2+-LmJean3 L. major promastigotes were harvested, centrifuged at 7000× g for 10 min, and then fixed using a 1% paraformaldehyde/PBS solution. The pellets were then washed twice with PBS, and the cells were resuspended in a 1 mg/mL DAPI (Sigma Aldrich, USA) solution for 30 min at 4 °C for kinetoplast and nucleus staining. The parasites were washed twice with PBS before visualization. Images from the slides were acquired using a PerkinElmer ultraVIEW confocal microscope, employing 405/400–600 nm, 488/500–560 nm, and 561/570–700 nm excitation/emission wavelengths.

2.11. Metacyclic Forms Isolation

Metacyclic parasites (unagglutinated cells) were isolated using the peanut agglutinin (PNA) method (Sacks et al., 1985) from 5 × 108 stationary parasites [92]. The cells were then counted using a Z1 Coulter counter (Beckman Coulter, Fullerton, CA, USA) and used for the extraction of their RNA.

2.12. In Vitro Infections

A total number of 3 × 105 RAW 264.7 murine macrophages per well were seeded in 6-well plates and grown at 37 °C in DMEM medium (Gibco Laboratories) supplemented with 10% (v/v) heat-inactivated FBS (Gibco Laboratories) and an antibiotic cocktail (50 U/mL penicillin, 50 mg/mL streptomycin) (Sigma). Macrophages were infected by pXG-LmJean3 and pXG-Hyg metacyclic parasites isolated by the PNA method [92] using a 25:1 (parasites/macrophage) ratio. The plates were incubated at 37 °C in a 5% CO2 atmosphere. After 24 h, non-phagocytosed parasites were washed with PBS, and fresh medium was added to the wells. The cells were harvested after 48 h by incubation in a trypsin–EDTA solution (Gibco Laboratories) for 5 min at 37 °C. Finally, the cells were centrifuged at 300× g for 10 min, and the pellets were washed twice with PBS before RNA extraction.

2.13. RNA Expression Quantification

For wild-type (WT) gene quantification, the RNA from log phase and stationary phase procyclic and metacyclic parasites from the in vitro infections of RAW 264.7 macrophages was extracted using a Qiagen RNAeasy extraction mini kit (Qiagen, Hilden, Germany). For the quantification of gene expression from transgenic parasites, the RNA from log-phase procyclic parasites was also extracted with the same kit.
RNA samples were then treated with Ambion DNA-free Kit (Invitrogen, Vilnius, Lithuania) following the manufacturer’s instructions. Retrotranscription was performed using 800 ng of RNA and M-MLV Reverse Transcriptase (Promega, Madison, WI, USA), following the protocol of the manufacturer. The obtained cDNA was then used for the qPCR assay. qPCR assays were performed using a 7500 real-time PCR system (Applied Biosystems, Foster City, CA, USA), 96-well plates (Applied Biosystems), and SYBR Green PCR master mix (Applied Biosystems). The primers used for qPCR are listed in Supplemental Table S3. The expression levels of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were used to normalize gene expression.

2.14. Cell Cycle Analysis by Propidium Iodide Staining and Flow Cytometry

A total amount of 106 log-phase parasites were harvested, washed twice with cold PBS, and fixed with 70% ethanol overnight. Then, the cells were centrifuged and resuspended in a solution of PBS with RNase A (Sigma) at a final concentration of 50 µg/mL for 20 min (37 °C). The cell suspensions were centrifuged and washed twice with PBS; the parasites were then stained with a 5 µg/mL solution of propidium iodide (Sigma). Finally, the cell suspensions were subjected to flow cytometry using an Attune flow cytometer (Applied Biosystems, CA, USA). Ten thousand events were collected for each sample. The number of gated cells in the G1, G2/M, and S-phase is presented as a percentage.

2.15. Cytotoxicity Evaluation

Log-phase L. major parasites were seeded in 96-well plates and grown at 26 °C with increasing concentrations of amphotericin B (Sigma), miltefosine (Calbiochem, Darmstadt, Germany), paromomycin (Sigma), and geneticin (Sigma), diluted in M199 medium. After 48 and 72 h of incubation, parasites’ viability was measured using 3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyl-tetrazolium bromide (MTT), as previously described [93]. The half-maximal effective concentration (EC50) was determined by fitting a sigmoidal Emax model to dose–response curves. Promastigotes viability was evaluated by comparison with untreated control cells. The results were expressed as means ± standard deviation (SD).

2.16. Statistical Analysis

Statistical analyses were carried out with GraphPad Prism v7.0b [94]. Two-group comparisons were performed by employing unpaired, two-tailed Student’s t-tests; p values < 0.05 were considered statistically significant (* p < 0.05, ** p < 0.01, *** p < 0.001).

3. Results

3.1. LmjF.22.0810 (LmJean3), a Novel Predicted Trypanosomatid Protein Kinase

LmjF.22.0810 was identified and retrieved from GeneDB (www.genedb.org) as one of the L. major genes encoding proteins harboring a kinase catalytic domain of the CAMK family (determined by phylogenetic inference [22]). Interestingly and on the basis of previous studies, LmjF.22.0810 had been described as a possible drug target [28]. LmjF.22.0810 orthologs are located on chromosome 22 of all Leishmania species whose genome was sequenced and included in TriTrypDB (www.tritrypdb.org). Gene sequence alignments displayed a high identity among Leishmania spp. (>85%) (Supplemental Figure S2). This novel gene encoding a protein with a kinase domain was presently named ‘LmJean3′. Finally, Sanger sequencing of LmJean3 from L. major (Lv39c5) was performed and showed a single nucleotide mutation (C282T) when compared to the gene sequence from the database (Supplemental Figure S3). Nevertheless, this polymorphism (C282T) does not change the translated protein sequence, given that both codons (AUC and AUT) encode the same amino acid (isoleucine).
Genes encoding protein orthologs were also found in T. cruzi and T. brucei on chromosomes 13 and 7, respectively. Gene sequence alignments displayed more than 50% identity when compared to the Trypanosoma spp. genes (Supplemental Figure S2).

3.2. LmJean3 (LmjF.22.0810) Domains and Motifs

LmJean3 is a 372-amino acid protein with an estimated molecular mass of 40.87 kDa and a theoretical pI of 6.60. A putative protein kinase domain was found between residues 8 and 266 (Figure 1A), containing all the currently known essential motifs and residues required for protein kinase activity among its 11 kinase subdomains (Supplemental Figure S4). Moreover, we also identified several sites and regions involved in the phosphorylation activity of kinases, such as ATP-binding site, catalytic loop, and activation segment (Figure 1A). No putative glycosylation, myristoylation, signal peptide, or GPI anchor was found within the sequence of LmJean3. Nonetheless, sequence analysis revealed EF-hand patterns at two C-terminal regions of the protein (244–272 and 333–361, with a percentage of similarity to the query PROSITE pattern of 83% and 81%, respectively; Figure 1A).
Additionally, we identified 17 serine and three threonine residues predicted to be putative phosphorylation sites in LmJean3 (Supplemental Figure S5). Two of those residues were found within the activation loop (S157 and S162). Furthermore, the region containing these residues (157S–X–X–X–X–S162) is conserved within Leishmania and shows similarity to the motifs of MAPKKs of animal, yeast (S–X–X–X–S/T), and plants (S/T–X–X–X–X–X–S/T) that turn the protein into the active state after phosphorylation [95]. On the other hand, in T. cruzi and T. brucei, the C-terminal serine from the motif is different. In T. cruzi, a glutamine is found at this position (Q161), while in T. brucei, a glutamate is present (E161), which may mimic the required phosphoresidue, analogously to other kinases such as Phk [96].
The kinase domain is followed by a less conserved disordered region (residues 275–372) (Figure 1A) that did not display significant similarity to other proteins after BLAST analysis. Furthermore, most of the amino acid substitutions may occur at the C-terminal region of Jean3 sequences, while the N-terminal region where the kinase domain is located is highly conserved (Figure 1B). On the contrary, Leishmania braziliensis Jean3 protein displayed at its N-terminal region an additional extension of 60 amino acids compared to its homologs from other trypanosomatids (Figure 1B).

3.3. The Phylogeny of LmjF.22.0810 Homologues

The kinase domains of non-kinetoplastid organisms (Arabidopsis thaliana, Homo sapiens, and Saccharomyces cerevisiae) were retrieved from a BLASTp search and were used to reconstruct the phylogram of LmjF.22.0810 homologues. Such phylogenetic reconstruction not only supported the aforementioned homology among trypanosomatids but also confirmed that LmjF.22.0810 orthologs are restricted to kinetoplastid organisms (Figure 1C). The Ser/Thr kinase domains of CIPK23 (calcineurin B-like interacting protein kinase 23), SOS2 (salt overly sensitive 2), Snfp1 (sucrose non-fermenting kinase 1), AMPK (AMP-activated protein kinase), MARK2 (microtubule affinity-regulating kinase 2), and TSSK3 (testis-specific serine kinase 3) displayed sequence homology to the catalytic domain of LmJean3 as determined by BLAST analysis and multiple-sequence alignment. However, the phylogenetic reconstruction did not show a clear orthologous relationship with other specific ePKs, except for those in the kinetoplastids branch (Figure 1C). The phylogenetic reconstruction also included the protein CDPK1 (calcium-dependent protein kinase 1) from Plasmodium falciparum, an extensively characterized protein with EF-hand motifs within its sequence [98]. As shown in Figure 1C and previously observed [22], LmJean3 kinase and its homologues do not belong to the clustering of CDPK.

3.4. LmJean3 was Localized in the Cytoplasm, Nucleus, and Flagellum of Leishmania Promastigotes and Was Significantly Expressed in the Infective and Amastigote Forms

The red and green fluorescent fusion protein constructs pXG-LmJean3-mCherry12 and pXG-GFP2+-LmJean3 were used to locate LmJean3 inside transgenic L. major parasites. Both fluorescent approaches (LmJean3-mCherry12 and GFP2+-LmJean3) allowed the visualization of LmJean3 in the cytosol, nucleus, and flagellum of the parasites (Figure 2A). Despite the observed localization of this L. major protein, analysis of the previously mentioned additional N-terminal 60-residue fragment from L. braziliensis by the TMHMM Server predicted a transmembrane region exclusive to the L. braziliensis protein (Supplemental Figure S6).
The expression levels of LmJean3 were quantified during the in vitro growth of L. major promastigotes and also from amastigotes derived from in vitro infections. The gene expression levels of LmJean3 increased ~2-fold for log-phase parasites as they reached the stationary phase as procyclic or metacyclic promastigotes (Figure 2B). After that, these expression levels were maintained in the infective form (metacyclic promastigotes) and the amastigote parasites (Figure 2B).

3.5. The Kinase Domain of LmJean3 Is Conserved and Exhibits a High Similarity to CIPK24/SOS2, CIPK23, and Snf1 Kinase Domains

Since no crystallographic structure of LmJean3 kinase homologues was available, homology modeling was chosen as an approach to generate the protein structure of LmJean3. The most suitable templates found in the Swiss-Model Repository were two CIPK protein kinases from A. thaliana (CIPK23 and CIPK24/SOS2) and Snf1 from S. cerevisiae. SOS2 and Snf1 are two kinases that belong to the CAMK family. The multiple-sequence alignment of critical regions for kinase phosphotransferase activity is shown in Figure 3. The catalytic loop, activation loop, glycine-rich loop, and ATP-binding region signature of Jean3 kinases were found to share a high homology with the possible templates (Figure 3). Moreover, catalytically important residues for kinases, such as R128 and D129 from the H–X–D triad [99], as well as D147 from the DFG (Asp–Phe–Gly) triplet, K37 from the A–X–K motif, and E53 from the C-helix [100], were conserved in all the analyzed sequences (Figure 3). Sequence identity was 34.22% between LmJean3 and Snf1, 37.26% between LmJean3 and CIPK23, and 41.83% between LmJean3 and SOS2. Consequently, sequence analysis indicated the protein structures of CIPK23, CIPK24/SOS2, and Snf1 represented suitable templates for homology modeling. The complete sequence alignment can be found in the Supplementary data (Supplemental Figure S7).

3.6. Homology Modeling, Refinement, and Validation of LmJean3 Structure

The Robetta webserver was used to generate the tertiary structure of the protein [54]. After Ginzu alignment of the LmJean3 kinase domain, Robetta determined CIPK23 as the best template for homology modeling (0.81 confidence). Consequently, an intermediate model was generated based on chain B of the 4CZT entry in the Protein Data Bank (PDB) repository. 4CZT corresponds to the crystal structure of the kinase domain of CIPK23 determined by Chaves-Sanjuan et al. (2014) using X-ray diffraction with a resolution of 2.3 Å [56]. Finally, energy minimization was performed with the YASARA force field to alleviate steric clashes and other peculiarities in the intermediate model [58]. The optimized resulting model was then validated. ERRAT determined the overall quality factor of the final model to be 98.11% (Figure 4A), as opposed to 95.49% of the template (4CZT) (Figure 4B). ProSA displayed a similar overall quality (Z-score) for LmJean3 (−7.66; Figure 4C) and the template (−8.47; Figure 4D). Negative Z-scores are required to consider the validity of structures [102]. The Ramachandran plot showed 88.4% of the residues to be localized in the favored regions, while none was found in the disallowed regions (Figure 4E).
Local model quality estimations by GROMOS and ANOLEA, in addition to LmJean3 structural features determined by DSSP [64,68,69,70], were also evaluated (Supplemental Figure S8). These data confirmed the validity of the model. Finally, the observed backbone Cα RMSD between LmJean3 and the template (4CZT_B) was 1.87 Å, while superposition of the five active site relevant residues (K37, E53, R128, D129, and D147) displayed an RMSD value of 0.58 Å (Figure 5A). Thus, the catalytically relevant residues were aligned to their equivalents in the template. A summary of the global quality of the optimized predicted structure is shown in Supplemental Table S4.

3.7. Overall Description of the Predicted Structure of LmJean3 Catalytic Domain

The structure of the predicted catalytic domain of LmJean3 displayed a canonical Ser/Thr kinase fold and was in accordance with the template used for homology modeling. The protein folded into two lobes (Figure 5B), with the catalytic cleft located between them. The core of the N-lobe consisted of five-stranded antiparallel β-sheets (β1-β5) and the C-helix. The glycine-rich loop (residues 15–20) was found within the N-lobe between β-strands 1b and 2. On the other hand, the C-lobe comprised six well-defined helices (αD-αI), with the catalytic loop (residues 125–137) and the activation loop (residues 147–180) within them. Inside the C-lobe, β8 and αT1 were found flanking the DFG motif, where D147 was inferred to recognize one of the ATP-bound Mg2+ ions [103,104]. Residues R128 and D129 from the H–X–D triplet were found inside the catalytic cleft, within the peptide substrate binding region (Figure 5B). Finally, the activation loop, which is delimited by two α-helical turns (αT1 and αT2), was found placing the G-loop aside from the active site (Figure 5B). This disposition of the activation loop impedes the access of protein substrates to the catalytic cleft and is a major indicator of the inactive state of a kinase [104]. Since kinase conformations are determinant for the specificity of small-molecule inhibitors [14,30,31,105], we decided to assess the local spatial pattern (LSP) alignment of LmJean3 catalytic and regulatory spines. LSP showed the alignment between both spines to be broken and the C-helix directed outward from the active site (αC-out) (Figure 5C). Additionally, as shown in Figure 5A, D147 from the DFG triplet appeared directed towards the ATP-binding site in what is known as the “DFG-in” conformation. Therefore, the LmJean3 structure displays a DFG-in/αC-out inactive conformational state, in agreement with the conformation of the template.

3.8. Molecular Dynamics Simulation

Over the course of the simulation, the RMSD per residue was calculated against the initial frame. The protein RMSD increased to a maximum deviation of 3.89 Å after 15.34 ns, with an average value of 2.94 ± 0.52 Å (Supplemental Figure S9). Importantly, the observed RMSD for the catalytically relevant residues of LmJean3 (K37, E53, R128, D129, and D147) was rather constant during the simulation (Supplemental Figure S9), having a maximum of 1.47 Å at 14.67 ns and an average value of 0.89 ± 0.15 Å.
The structure obtained from the average atomic coordinates of the trajectory was aligned against the predicted catalytic residues of the initial LmJean3 structure (Figure 5D). The catalytically relevant residues of both structures displayed an RMSD value of 1.15 Å. Hence, their structural resemblance throughout the duration of the simulation provides confidence in the homology model structure.
Non-covalent residue interactions such as salt bridges and hydrogen bonds were also calculated during the simulation (Supplemental Tables S1 and S5). The main interactions found within the kinase domain and those involving the catalytically relevant residues, as well as their persistence during the trajectory, are described in Table 1. The observed non-covalent interactions for LmJean3 were similar to those reported for other kinases such as Src, PKA, and Twitching [106,107,108]. Figure 5E shows that R254 from subdomain XI formed a highly conserved stabilizing salt bridge with E180 from the APE motif [107]. Similarly, an intracatalytic salt bridge, known for stabilizing the catalytic loop, was observed in subdomain VIb between K131 and D129 (Table 1; Figure 5E) [107,108]. As a consequence of the αC-out conformation found in the structure, R152 from the activation segment and E53 from the αC-helix were found bound by a salt bridge [108]. Conversely, K37 from the A–X–K motif formed a salt bridge with D147, indicating an inward direction of DFG–Asp towards the active site, a feature known as one of the prerequisites for switching into the active state (“DFG-in” conformation) (Table 1; Figure 5E). Additionally, residues L137–G88, R128–D192, and N134–D129 displayed several stabilizing hydrogen bonds with a persistence higher than 40% (Table 1), equivalent to interactions reported for other protein kinases [106,107,108,109].

3.9. LmJean3 Binding Sites: Prediction and Analysis

The top-five ranked potential ligand-binding sites of LmJean3 (named site A, B, C, D, and E) were identified using SiteMap v3.8 [82,83] and are summarized in Supplemental Table S6. Site A (SiteScore = 1.05) corresponds to the solvent-exposed regions of the allosteric site, the ATP-binding pocket, and deeper cavities of the ATP site (Figure 6A). Site B (SiteScore = 0.92) is located in a smaller allosteric cleft found at the protein substrate binding site (Figure 6B). This site includes C-lobe residues of the catalytic loop, the activation loop, the P+1 loop, A188, A191, and S195 of the F-helix, and D49 and V46 of the αC-helix (Figure 6B). Interestingly, residues of the catalytic loop were found in Site B, not obstructed by the activation loop. The distribution of Site A appeared wider than that of Site B, similarly to what has been described for other DFG-in/αC-out conformations [30]. Both sites (A and B) displayed more hydrophilic than hydrophobic components (Figure 6 and Supplemental Table S6). On the other hand, sites C, D, and E did not surpass the recommended SiteScore cut-off of 0.80 [83]. The druggability (Dscore) of sites A and B was also evaluated using SiteMap. Site A displayed a Dscore of 1.03, and Site B a Dscore of 0.85, both being promising values for targeting by small-molecule, drug-like ligands [83]. Due to the low selectivity that ATP-binding-site inhibitors typically display [14,30,31,104] and to minimize any off-target activity [18], we decided to target the allosteric site found in the catalytic pocket of LmJean3 (Figure 6B).

3.10. Aminoglycosides Are Predicted Ligands of LmJean3.

More than 11,400 molecules from the Zim dataset (April 2016) from ZINC [86] were used for in silico docking to the Site B of LmJean3. Zim (ZINC in man) is a subset of experimental compounds, including drugs, which have been approved for use in humans. Eleven compounds were predicted to bind LmJean3, with predicted docking scores ranging from −9.06 to −11.46 kcal/mol (within GLIDE’s error margin of 2 kcal/mol) (Supplemental Table S7). Seven of the best solutions were AG antibiotics, three were flavonoids, and one was the anthracenedione mitoxantrone. The seven best docking solutions were ranked by their calculated docking score (Table 2). Limitations of the database, such as stereochemical ambiguity [86], were considered. In fact, only stereoisomers of neomycin B (NEO) and amikacin (AMK) were predicted to bind to Site B (Table 2). On the other hand, paromomycin (PMM) and tobramycin (TBR) were docked as their natural purified structures (Table 2).
Four of the docked compounds from Table 2 are either AGs (PMM and TBR) or isomers of AGs (NEO and AMK). In Leishmania, AGs are supposed to bind directly to the eukaryotic A-site of the small ribosomal subunit (40S) [110]; however, their mechanism of action is not completely clear. The four docked AGs contain a 2-deoxystreptamine core (2-DOS; Ring II; Figure 7) and can be classified on the basis of their substitution positions in the 2-DOS ring. PMM and NEO 2-DOS groups are substituted at positions 4 and 5 (Figure 7A), whereas AMK and TBR have a 4,6-disubstituted 2-DOS group (Figure 7B). Additionally, only PMM displays a 6′-OH substituent group at Ring I, while the other docked AGs are 6′-NH2 derivatives (green-colored circles in Figure 7). The nature of the substituent group at the 6′-position of Ring I has been described to be an important component of selectivity toward eukaryotic species [111]. The aforementioned statement is in agreement with previous works in Leishmania according to which NEO, AMK, and TBR (6′-NH2 derivatives) displayed either high EC50 values in vitro or incomplete healing in vivo, in contrast to PMM and geneticin (G418) which are 6′-OH derivatives [110,112,113].

3.11. Predicted Interactions of Paromomycin with Site B of LmJean3

The predicted interactions of PMM within Site B are presented in Figure 7C (DockScore of −11.46 kcal/mol). The 2-DOS ring from PMM was found hydrogen-bonded to residues from the H–X–D triplet, as well as to residues in the activation loop and the αC-helix. Three hydrogen bonds were found between R128 and the PMM 2-DOS ring (Ring II). Also, D129 formed a hydrogen bond and a salt bridge with the protonated amine moiety at position 1 from Ring II. R155 in the activation loop was found to form a hydrogen bond with the 6′-OH substituent group from Ring I, while PMM Ring IV was hydrogen-bonded to E156 and R152. Finally, D49 from the αC-helix displayed an H-bond and a salt bridge with position 6′′′ of Ring IV. Overall, these results showed that residues from the H–X–D triplet and C-terminal to the DFG motif, two key regions for substrate binding and γ-phosphate transfer [99,114,115], may play a crucial role in PMM binding to Site B.

3.12. Docking in Site A and Cross-Docking

The top-ranked molecules resulting from the Site B docking analysis as well as geneticin were also docked into Site A. Moreover, the detergent molecule CHAPS (C32H58N2O7S, HET-ID CPS, 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate), found in the ATP-binding site of the template CIPK23 during the crystallization process, and ATP were also analyzed to predict their docking score to Site A (self-dock for ATP) and to Site B.
As expected, ATP showed a stronger docking score for Site A than for Site B (Figure 8). Geneticin displayed a weaker docking score for both sites as compared to those of other dataset molecules. On the other hand, PMM displayed the strongest docking score (lowest free energy) for both sites (Figure 8). These data show a predicted strong putative interaction between PMM and LmJean3. Such predicted interaction was stronger than that of other cofactors and possible ligands.

3.13. Generation of LmJean3-Overexpressing Parasites (LmJ3OE)

We next assessed the activity of some of the docked drugs onto L. major highly expressing LmJean3. By using the constructed expression vector pXG-LmJean3, we developed a transgenic L. major strain overexpressing LmJean3 (LmJ3OE). qPCR revealed that LmJ3OE parasites showed a 30-fold higher LmJean3 gene expression level compared to the controls (pXG-Hyg parasites) (Supplemental Figure S10). In addition, the cell cycle distribution and growth rate of LmJ3OE parasites were similar to those observed for the controls (Supplemental Figure S10).
The transcriptional profiles of several genes related to cell morphology, proliferation, and treatment resistance were also evaluated in promastigotes.
Expression levels of ABC-transporter genes (ABCA3, ABCH1, ATPase 1-like, ABCC2, and MRPA) and α-tubulin (an important component of cell cytoskeleton) were measured in log-phase promastigotes from LmJ3OE and control parasites (Supplemental Figure S11). No significant alteration was detected between the two groups, with the exception of ABCA3 (LmjF.11.1240), which was significantly downregulated in LmJ3OE parasites (p < 0.001; Supplemental Figure S11).

3.14. LmJ3OE Parasites Exhibited Less Sensitivity to Paromomycin and Other Aminoglycosides

Since PMM, a first-line treatment against leishmaniasis, was previously proposed to be a putative ligand of LmJean3, we investigated this hypothesis.
Through the analysis of the growth inhibition of LmJ3OE and control parasites after PMM treatment, EC50 values were estimated. Our data showed that LmJ3OE promastigotes were significantly less sensitive to PMM after 48 h and 72 h of exposure (~2-fold increase of EC50) (Figure 9A). In addition to PMM (a 4,5-disubstituted aminoglycoside), we also tested geneticin (a 4,6-disubstituted aminoglycoside and another 6′-OH derivative). Our results confirmed LmJ3OE parasites exhibited less sensitivity to paromomycin and geneticin (Figure 9A,B; Supplemental Table S8). On the contrary, compared to the controls, LmJ3OE promastigotes displayed a significantly higher sensitivity to both leishmanicidal compounds amphotericin B and miltefosine (Figure 9C,D; Supplemental Table S8), predicted not to interact with our target.

4. Discussion

PKs are key regulators of cellular signaling and have long been recognized as important drug targets in a variety of diseases (such as cancer and inflammatory and infectious diseases) [14,15,16,17]. LmjF.22.0810 was identified to encode a protein kinase that we named “LmJean3”. This PK harbors EF-hand sequences, which are intrinsic motifs for calcium-dependent protein kinases (CDPKs). CDPKs are a prominent group of kinases found in plants and parasites, such as P. falciparum and Toxoplasma gondii, but absent in human and yeast kinomes [22]. However, the phylogenetic inference of LmJean3 kinase homologues did not support its clustering with CDPKs. Additionally, the known regulatory motifs of CAMKs were also absent in LmJean3. Thus, as previously mentioned by Parsons et al. (2005), it is likely that Jean3 genes encode a novel class of ePKs containing EF hands [22].
LmJean3 was found to be located in the flagellum, cytoplasm, and nucleus of L. major. On the basis of the high protein sequence identity found in all the Leishmania species analyzed, LmJean3 orthologs are expected to display a similar localization. Since trypanosomatids lack tyrosine kinases, it has been proposed that receptor Ser/Thr kinases could respond to host or parasite ligands and mediate intercellular communication [20,22]. Moreover, the predicted N-terminal transmembrane region of Jean3 from L. braziliensis may represent an additional feature that could be further evaluated in order to identify biological functions of LmJean3 kinase homologues. However, the primary focus for this research study was to characterize LmJean3 as a drug target.
In the first step, computations were performed to predict a 3D structure for LmJean3 suitable for molecular docking. The evaluation of the protein structure and molecular dynamics simulations showed that the LmJean3 model and its catalytic residues were stable, providing validation of the protein structure. We also showed in silico that PMM and TBR may bind to LmJean3 with multiple hydrogen bonds, which are major components of stable protein–ligand complexes [116]. Cross-docking indicated that the strength of the interaction between PMM and the protein is stronger than that of the interaction with other cofactors and possible ligands, including those co-crystallized in the original template. Finally, to complement the in silico studies, we evaluated the cytotoxicity of PMM and geneticin in transgenic LmJean3-overexpressing L. major parasites (LmJ3OE). Geneticin has an analog structure to TBR, with a 4,6-disubstituted 2-DOS ring; on the contrary, PMM bears a 4,5-disubstituted 2-DOS ring. Nonetheless, previous studies have shown that the leishmanicidal activity of both compounds is mainly due to their 6′-OH substituent group in Ring I [110]. In vitro tests determined that LmJ3OE parasites’ sensitivity is significantly lower to both drugs. In agreement with this finding, the overexpression of a specific target has been described as a putative mechanism for drug resistance in eukaryotic pathogens [6].
In bacteria, the enzymatic modification of AG antibiotics is the most prevalent type of resistance found in the clinical setting [3]. It may result from the addition of an acetyl, a nucleotydil, or a phosphate group by acetyltransferases, nucleotidyltransferases, or phosphotransferases, respectively [117]. It is remarkable that aminoglycoside phosphotransferases (APHs) are proteins that share structural and functional homology with ePKs but display very little primary sequence conservation (<5%) [2,118]. It is also interesting to highlight that some APHs from pathogenic bacteria confer resistance by sequestering AGs by tight binding, rather than by inactivation of the target [119,120]. On the basis of this strategy, and in agreement with the docking results and cytotoxicity evaluation, we postulate that Jean3 may be capable of sequestering PMM and geneticin.
Additionally, it has been proposed that PMM might not only interact with ribosomes but also alter membrane fluidity, interfere with the mitochondrial membrane potential, and inhibit respiration [121,122,123]. Hence, we measured the gene expression levels of PCNA, α-tubulin, and various ABC-transporters. In Leishmania, those molecules had previously been related to resistance against several drugs such as miltefosine, antimony, and PMM [122,124,125]. Interestingly, in the LmJ3OE strain, we found that ABCA3 (LmjF.11.1240) was highly downregulated. ABCA3 overexpression has been observed in antimony-resistant strains; however, transfection of the gene did not confer a drug-resistant phenotype, excluding it as the unique factor determining antimony resistance [126]. On the contrary to antimony-resistant strains, the expression of ABCA3 has been observed to be inversely correlated with miltefosine susceptibility in Leishmania (Viannia) [127]. Nevertheless, in LmJ3OE parasites, the low expression of ABCA3 did not correlate directly with their increased susceptibility to miltefosine. These data seem to indicate the existence of a broader mechanism underlying the observed drug response in LmJ3OE parasites.
Interestingly, the low number of cases of patients infected with PMM-resistant strains is linked to its occasional use as a clinical treatment for visceral leishmaniasis [7]. PMM has been used more widely for the treatment of the cutaneous disease, in which two cases caused by PMM-resistant L. aethiopica reported three- to five-fold drug sensitivity decrease [7,128]. A review of published research data from drug-resistant strains does not implicate Jean3 transcript or its coding protein in PMM resistance. Therefore, it would be insightful to further analyze LmJean3 homologues transcripts in drug-resistant strains.
The predicted complexes of LmJean3 with different AGs within its protein binding site rationalize the structural basis for this enzyme broad-spectrum AG resistance activity. Furthermore, the generated predicted model of LmJean3 could be exploited for structure-based drug design of compounds to combat leishmaniasis. From the perspective of drug resistance, the examination of AG binding to LmJean3 homologues may highlight avenues for the development of next-generation aminoglycosidic drugs.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-273X/9/11/723/s1, Figure S1: Schematic representation of the cloning approach to generate the expression vectors, Figure S2: LmjF.22.0810 orthologs are highly conserved, Figure S3: LmJean3 sequencing study, Figure S4: LmjF.22.0810 encodes a putative Ser/Thr kinase and contains all the essential motifs and residues required for protein kinase activity among its 11 kinase subdomains, Figure S5: LmJean3 phosphorylation sites prediction, Figure S6: L. braziliensis Jean3 transmembrane helices prediction, Figure S7: Multiple-sequence alignment of the kinase domain sequences from LmJean3, LdJean3, LbJean3, CIPK24/SOS22, and Snf1, Figure S8: ANOLEA and GROMOS quality estimation and DSSP defined secondary structure of LmJean3 model, Figure S9: Monitoring of MD trajectories, Figure S10: Generation of LmJean3-overexpressing parasites, Figure S11: Gene expression analysis from LmJ3OE and control parasites, Table S1: H-bonds found in LmJean3 kinase domain at least in 30% of the trajectory, Table S2: Primer sequences used for the construction of expression plasmids, Table S3: Primer sequences used for qPCR analyses, Table S4: LmJean3 global quality estimation, Table S5: Predicted salt bridges in LmJean3 kinase domain, Table S6: LmJean3 SiteMap property values and Dscorea, Table S7: LmJean3 site B docking results*, Table S8: Drug activity profile of L. major cell lines.

Author Contributions

Conceptualization, A.V. and P.A.N.; methodology, P.A.N., C.F.-R., E.L., A.V., A.T.G.-S., and F.R.F.; software, A.V. and A.T.G.-S.; validation, A.V., C.F.-R., M.A., E.L., and J.P.-G.; formal analysis, A.V.; investigation, A.V., C.F.-R., M.A., E.L., J.P.-G., and P.A.N.; resources, P.A.N.; writing—original draft preparation, A.V. and P.A.N.; writing—review and editing, P.A.N., A.V., and A.T.G.-S.; visualization, A.V.; supervision, P.A.N., F.R.F., and A.T.G.-S.; project administration, P.A.N.; funding acquisition, P.A.N.

Funding

This research was funded by Obra Social la Caixa and Fundación Caja Navarra, Gobierno de Navarra Salud (12/2017), Fundación Roviralta, Ubesol, Inversiones Garcilaso de la Vega S. L., by Government of Navarre and Laser Ebro. J.P.-G. was supported by a Ministerio de Educacion Cultura y Deporte fellowship (FPU17/03304). A.T.G.-S. thanks Haridus- ja teadusministeerium for grant IUT34-14.

Acknowledgments

We acknowledge Susana Murillo and Julia Villarreal from the University of Navarra for administrative support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wiwanitkit, V. Interest in paromomycin for the treatment of visceral leishmaniasis (Kala-azar). Ther. Clin. Risk Manag. 2012, 8, 323–328. [Google Scholar] [CrossRef] [PubMed]
  2. Daigle, D.M.; McKay, G.A.; Thompson, P.R.; Wright, G.D. Aminoglycoside antibiotic phosphotransferases are also serine protein kinases. Chem. Biol. 1999, 6, 11–18. [Google Scholar] [CrossRef]
  3. Ramirez, M.S.; Tolmasky, M.E. Aminoglycoside modifying enzymes. Drug Resist. Updates 2010, 13, 151–171. [Google Scholar] [CrossRef] [PubMed]
  4. Stogios, P.J.; Shakya, T.; Evdokimova, E.; Savchenko, A.; Wright, G.D. Structure and function of APH(4)-Ia, a hygromycin B resistance enzyme. J. Biol. Chem. 2011, 286, 42786. [Google Scholar] [CrossRef]
  5. Sundar, S. Drug resistance in Indian visceral leishmaniasis. Trop. Med. Int. Health 2001, 6, 849–854. [Google Scholar] [CrossRef]
  6. Fairlamb, A.H.; Gow, N.A.R.; Matthews, K.R.; Waters, A.P. Drug resistance in eukaryotic microorganisms. Nat. Microbiol. 2016, 1, 16092. [Google Scholar] [CrossRef]
  7. Chemother, J.A. Drug resistance in. Commun. Dis. Rep. Cdr. Rev. 2000, 19, 149–150. [Google Scholar]
  8. Clayton, C.E. Life without transcriptional control? From fly to man and back again. EMBO J. 2002, 21, 1881–1888. [Google Scholar] [CrossRef]
  9. Wiese, M. A mitogen-activated protein (MAP) kinase homologue of Leishmania mexicana is essential for parasite survival in the infected host. EMBO J. 1998, 17, 2619–2628. [Google Scholar] [CrossRef]
  10. Nett, I.R.E.; Martin, D.M.A.; Miranda-Saavedra, D.; Lamont, D.; Barber, J.D.; Mehlert, A.; Ferguson, M.A.J. The phosphoproteome of bloodstream form Trypanosoma brucei, causative agent of African sleeping sickness. Mol. Cell. Proteom. 2009, 8, 1527–1538. [Google Scholar] [CrossRef]
  11. Nakayasu, E.S.; Gaynor, M.R.; Sobreira, T.J.P.; Ross, J.A.; Almeida, I.C. Phosphoproteomic analysis of the human pathogen Trypanosoma cruzi at the epimastigote stage. Proteomics 2009, 9, 3489–3506. [Google Scholar] [CrossRef] [PubMed]
  12. Queiroz, R.M.L.; Charneau, S.; Mandacaru, S.C.; Schwämmle, V.; Lima, B.D.; Roepstorff, P.; Ricart, C.A.O. Quantitative proteomic and phosphoproteomic analysis of Trypanosoma cruzi amastigogenesis. Mol. Cell. Proteom. 2014, 13, 3457–3472. [Google Scholar] [CrossRef] [PubMed]
  13. Tsigankov, P.; Gherardini, P.F.; Helmer-Citterich, M.; Späth, G.F.; Zilberstein, D. Phosphoproteomic analysis of differentiating Leishmania parasites reveals a unique stage-specific phosphorylation motif. J. Proteome Res. 2013, 12, 3405–3412. [Google Scholar] [CrossRef] [PubMed]
  14. Wu, P.; Nielsen, T.E.; Clausen, M.H. Small-molecule kinase inhibitors: An analysis of FDA-approved drugs. Drug Discov. Today 2016, 21, 5–10. [Google Scholar] [CrossRef] [PubMed]
  15. Diaz, R.; Nguewa, P.A.; Redrado, M.; Manrique, I.; Calvo, A. Sunitinib reduces tumor hypoxia and angiogenesis, and radiosensitizes prostate cancer stem-like cells. Prostate 2015, 75, 1137–1149. [Google Scholar] [CrossRef]
  16. Ibanez, E.; Agliano, A.; Prior, C.; Nguewa, P.; Redrado, M.; Gonzalez-Zubeldia, I.; Plano, D.A.; Palop, J.; Sanmartin, C.; Calvo, A. The Quinoline Imidoselenocarbamate EI201 Blocks the AKT/mTOR Pathway and Targets Cancer Stem Cells Leading to a Strong Antitumor Activity. Curr. Med. Chem. 2012, 19, 3031–3043. [Google Scholar] [CrossRef]
  17. Nguewa, P.A.; Calvo, A.; Pullamsetti, S.S.; Banat, G.A.; Grimminger, F.; Savai, R. Tyrosine kinase inhibitors with antiangiogenic properties for the treatment of non-small cell lung cancer. Expert Opin. Investig. Drugs 2011, 20, 61–74. [Google Scholar] [CrossRef]
  18. Breen, M.E.; Soellner, M.B. Small molecule substrate phosphorylation site inhibitors of protein kinases: Approaches and challenges. ACS Chem. Biol. 2015, 10, 175–189. [Google Scholar] [CrossRef]
  19. Crowther, G.J.; Shanmugam, D.; Carmona, S.J.; Doyle, M.A.; Hertz-Fowler, C.; Berriman, M.; Nwaka, S.; Ralph, S.A.; Roos, D.S.; van Voorhis, W.C.; et al. Identification of attractive drug targets in neglected- disease pathogens using an in Silico approach. PLoS Negl. Trop. Dis. 2010, 4, e804. [Google Scholar] [CrossRef]
  20. Naula, C.; Parsons, M.; Mottram, J.C. Protein kinases as drug targets in trypanosomes and Leishmania. Biochim. Biophys. Acta-Proteins Proteom. 2005, 1754, 151–159. [Google Scholar] [CrossRef]
  21. Merritt, C.; Silva, L.E.; Tanner, A.L.; Stuart, K.; Pollastri, M.P. Kinases as druggable targets in trypanosomatid protozoan parasites. Chem. Rev. 2014, 114, 11280–11304. [Google Scholar] [CrossRef] [PubMed]
  22. Parsons, M.; Worthey, E.A.; Ward, P.N.; Mottram, J.C. Comparative analysis of the kinomes of three pathogenic trypanosomatids: Leishmania major, Trypanosoma brucei and Trypanosoma cruzi. BMC Genom. 2005, 6, 127. [Google Scholar] [CrossRef] [PubMed]
  23. Klug, D.M.; Gelb, M.H.; Pollastri, M.P. Repurposing strategies for tropical disease drug discovery. Bioorganic Med. Chem. Lett. 2016, 26, 2569–2576. [Google Scholar] [CrossRef] [PubMed]
  24. Kufareva, I.; Abagyan, R. Type-II kinase inhibitor docking, screening, and profiling using modified structures of active kinase states. J. Med. Chem. 2008, 51, 7921–7932. [Google Scholar] [CrossRef] [PubMed]
  25. Danilenko, V.N.; Simonov, A.Y.; Lakatosh, S.A.; Kubbutat, M.H.G.; Totzke, F.; Schächtele, C.; Elizarov, S.M.; Bekker, O.B.; Printsevskaya, S.S.; Luzikov, Y.N.; et al. Search for inhibitors of bacterial and human protein kinases among derivatives of diazepines[1,4] annelated with maleimide and indole cycles. J. Med. Chem. 2008, 51, 7731–7736. [Google Scholar] [CrossRef] [PubMed]
  26. Diaz, R.; Nguewa, P.A.; Parrondo, R.; Perez-Stable, C.; Manrique, I.; Redrado, M.; Catena, R.; Collantes, M.; Peñuelas, I.; Díaz-González, J.A.; et al. Antitumor and antiangiogenic effect of the dual EGFR and HER-2 tyrosine kinase inhibitor lapatinib in a lung cancer model. BMC Cancer 2010, 10, 188. [Google Scholar] [CrossRef] [PubMed]
  27. Diaz, R.; Nguewa, P.A.; Diaz-Gonzalez, J.A.; Hamel, E.; Gonzalez-Moreno, O.; Catena, R.; Serrano, D.; Redrado, M.; Sherris, D.; Calvo, A. The novel Akt inhibitor Palomid 529 (P529) enhances the effect of radiotherapy in prostate cancer. Br. J. Cancer 2009, 100, 932–940. [Google Scholar] [CrossRef]
  28. Flórez, A.F.; Park, D.; Bhak, J.; Kim, B.C.; Kuchinsky, A.; Morris, J.H.; Espinosa, J.; Muskus, C. Protein network prediction and topological analysis in Leishmania major as a tool for drug target selection. BMC Bioinform. 2010, 11, 484. [Google Scholar] [CrossRef]
  29. Fang, Z.; Grütter, C.; Rauh, D. Strategies for the selective regulation of kinases with allosteric modulators: Exploiting exclusive structural features. ACS Chem. Biol. 2013, 8, 58–70. [Google Scholar] [CrossRef]
  30. Roskoski, R. Classification of small molecule protein kinase inhibitors based upon the structures of their drug-enzyme complexes. Pharmacol. Res. 2016, 103, 26–48. [Google Scholar] [CrossRef]
  31. Müller, S.; Chaikuad, A.; Gray, N.S.; Knapp, S. The ins and outs of selective kinase inhibitor development. Nat. Chem. Biol. 2015, 11, 818–821. [Google Scholar] [CrossRef]
  32. Lamba, V.; Ghosh, I. New Directions in Targeting Protein Kinases: Focusing Upon True Allosteric and Bivalent Inhibitors. Curr. Pharm. Des. 2012, 18, 2936–2945. [Google Scholar] [CrossRef]
  33. Docampo, R.; Huang, G. Calcium signaling in trypanosomatid parasites. Cell Calcium 2015, 57, 194–202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Aslett, M.; Aurrecoechea, C.; Berriman, M.; Brestelli, J.; Brunk, B.P.; Carrington, M.; Depledge, D.P.; Fischer, S.; Gajria, B.; Gao, X.; et al. TriTrypDB: A functional genomic resource for the Trypanosomatidae. Nucleic Acids Res. 2009, 38, D457–D462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  36. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
  37. Rice, P.; Longden, L.; Bleasby, A. EMBOSS: The European Molecular Biology Open Software Suite. Trends Genet. 2000, 16, 276–277. [Google Scholar] [CrossRef]
  38. Marchler-Bauer, A.; Anderson, J.B.; Chitsaz, F.; Derbyshire, M.K.; Deweese-Scott, C.; Fong, J.H.; Geer, L.Y.; Geer, R.C.; Gonzales, N.R.; Gwadz, M.; et al. CDD: Specific functional annotation with the Conserved Domain Database. Nucleic Acids Res. 2009, 37, D205–D210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Marchler-Bauer, A.; Lu, S.; Anderson, J.B.; Chitsaz, F.; Derbyshire, M.K.; DeWeese-Scott, C.; Fong, J.H.; Geer, L.Y.; Geer, R.C.; Gonzales, N.R.; et al. CDD: A Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39, D225–D229. [Google Scholar] [CrossRef] [Green Version]
  40. Yachdav, G.; Kloppmann, E.; Kajan, L.; Hecht, M.; Goldberg, T.; Hamp, T.; Hönigschmid, P.; Schafferhans, A.; Roos, M.; Bernhofer, M.; et al. PredictProtein—An open resource for online prediction of protein structural and functional features. Nucleic Acids Res. 2014, 42, W337–W343. [Google Scholar] [CrossRef] [Green Version]
  41. Finn, R.D.; Attwood, T.K.; Babbitt, P.C.; Bateman, A.; Bork, P.; Bridge, A.J.; Chang, H.Y.; Dosztanyi, Z.; El-Gebali, S.; Fraser, M.; et al. InterPro in 2017-beyond protein family and domain annotations. Nucleic Acids Res. 2017, 45, D190–D199. [Google Scholar] [CrossRef] [PubMed]
  42. Finn, R.D.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Mistry, J.; Mitchell, A.L.; Potter, S.C.; Punta, M.; Qureshi, M.; Sangrador-Vegas, A.; et al. The Pfam protein families database: Towards a more sustainable future. Nucleic Acids Res. 2016, 44, D279–D285. [Google Scholar] [CrossRef] [PubMed]
  43. Zhou, Y.; Yang, W.; Kirberger, M.; Lee, H.W.; Ayalasomayajula, G.; Yang, J.J. Prediction of EF-hand calcium-binding proteins and analysis of bacterial EF-hand proteins. Proteins Struct. Funct. Genet. 2006, 65, 643–655. [Google Scholar] [CrossRef]
  44. Pierleoni, A.; Martelli, P.L.; Casadio, R. PredGPI: A GPI-anchor predictor. BMC Bioinform. 2008, 9, 392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Wilkins, M.R.; Gasteiger, E.; Bairoch, A.; Sanchez, J.C.; Williams, K.L.; Appel, R.D.; Hochstrasser, D.F. Protein identification and analysis tools in the ExPASy server. Methods Mol. Biol. 1999, 112, 531–552. [Google Scholar] [PubMed]
  46. Blom, N.; Gammeltoft, S.; Brunak, S. Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. J. Mol. Biol. 1999, 294, 1351–1362. [Google Scholar] [CrossRef]
  47. Petersen, T.N.; Brunak, S.; Von Heijne, G.; Nielsen, H. SignalP 4.0: Discriminating signal peptides from transmembrane regions. Nat. Methods 2011, 8, 785–786. [Google Scholar] [CrossRef]
  48. Harrington, B.; Engelen, J. Inkscape. Available online: http://www.inkscape.Org (accessed on 20 Ocotrber 2009).
  49. Katoh, K. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [Green Version]
  50. Abascal, F.; Zardoya, R.; Posada, D. ProtTest: Selection of best-fit models of protein evolution. Bioinformatics 2005, 21, 2104–2105. [Google Scholar] [CrossRef] [Green Version]
  51. Le, S.Q.; Gascuel, O. An improved general amino acid replacement matrix. Mol. Biol. Evol. 2008, 25, 1307–1320. [Google Scholar] [CrossRef] [Green Version]
  52. Arnold, K.; Bordoli, L.; Kopp, J.; Schwede, T. The SWISS-MODEL workspace: A web-based environment for protein structure homology modelling. Bioinformatics 2006, 22, 195–201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Chenna, R.; Sugawara, H.; Koike, T.; Lopez, R.; Gibson, T.J.; Higgins, D.G.; Thompson, J.D. Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31, 3497–3500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Kim, D.E.; Chivian, D.; Baker, D. Protein structure prediction and analysis using the Robetta server. Nucleic Acids Res. 2004, 32, W526–W531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Berman, H.M.; Battistuz, T.; Bhat, T.N.; Bluhm, W.F.; Bourne, P.E.; Burkhardt, K.; Feng, Z.; Gilliland, G.L.; Iype, L.; Jain, S.; et al. The protein data bank. Acta Crystallogr. Sect. D Biol. Crystallogr. 2002, 58, 899–907. [Google Scholar] [CrossRef]
  56. Chaves-Sanjuan, A.; Sanchez-Barrena, M.J.; Gonzalez-Rubio, J.M.; Moreno, M.; Ragel, P.; Jimenez, M.; Pardo, J.M.; Martinez-Ripoll, M.; Quintero, F.J.; Albert, A. Structural basis of the regulatory mechanism of the plant CIPK family of protein kinases controlling ion homeostasis and abiotic stress. Proc. Natl. Acad. Sci. USA 2014, 111, E4532–E4541. [Google Scholar] [CrossRef] [Green Version]
  57. Vinet, L.; Zhedanov, A. A “missing” family of classical orthogonal polynomials. J. Phys. A Math. Theor. 2011, 44, 085201. [Google Scholar] [CrossRef]
  58. Krieger, E.; Joo, K.; Lee, J.; Lee, J.; Raman, S.; Thompson, J.; Tyka, M.; Baker, D.; Karplus, K. Improving physical realism, stereochemistry, and side-chain accuracy in homology modeling: Four approaches that performed well in CASP8. Proteins Struct. Funct. Bioinform. 2009, 77, 114–122. [Google Scholar] [CrossRef] [Green Version]
  59. Eisenberg, D.; Lüthy, R.; Bowie, J.U. VERIFY3D: Assessment of protein models with three-dimensional profiles. Methods Enzymol. 1997, 277, 396–404. [Google Scholar]
  60. Benkert, P.; Biasini, M.; Schwede, T. Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics 2011, 27, 343–350. [Google Scholar] [CrossRef]
  61. Benkert, P.; Künzli, M.; Schwede, T. QMEAN server for protein model quality estimation. Nucleic Acids Res. 2009, 37, W510–W514. [Google Scholar] [CrossRef] [Green Version]
  62. Benkert, P.; Tosatto, S.C.E.; Schomburg, D. QMEAN: A comprehensive scoring function for model quality assessment. Proteins Struct. Funct. Genet. 2008, 71, 261–277. [Google Scholar] [CrossRef] [PubMed]
  63. Benkert, P.; Schwede, T.; Tosatto, S.C. QMEANclust: Estimation of protein model quality by combining a composite scoring function with structural density information. BMC Struct. Biol. 2009, 9, 35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Melo, F.; Devos, D.; Depiereux, E.; Feytmans, E. ANOLEA: A www server to assess protein structures. Proc. Int. Conf. Intell. Syst. Mol. Biol. 1997, 5, 187–190. [Google Scholar] [PubMed]
  65. Dong, Q.W.; Wang, X.L.; Lin, L. Novel knowledge-based mean force potential at the profile level. BMC Bioinform. 2006, 7, 207–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Melo, F.; Feytmans, E. Assessing protein structures with a non-local atomic interaction energy. J. Mol. Biol. 1998, 277, 1141–1152. [Google Scholar] [CrossRef] [PubMed]
  67. Laskowski, R.A.; MacArthur, M.W.; Moss, D.S.; Thornton, J.M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 1993, 26, 283–291. [Google Scholar] [CrossRef]
  68. van Gunsteren, W.F.; Billeter, S.R.; Eising, A.A.; Hünenberger, P.H.; Krüger, P.; Mark, A.E.; Scott, W.R.P.; Tironi, I.G. Biomolecular Simulation: The {GROMOS96} Manual and User Guide; Vdf Hochschulverlag AG an der ETH Zürich: Zürich, Switzerland, 1996. [Google Scholar]
  69. Kabsch, W.; Sander, C. Dictionary of protein secondary structure: Pattern recognition of hydrogen--bonded and geometrical features. Biopolymers 1983, 22, 2577–2637. [Google Scholar] [CrossRef]
  70. Hutchinson, E.G.; Thornton, J.M. PROMOTIF-A program to identify and analyze structural motifs in proteins. Protein Sci. 2008, 5, 212–220. [Google Scholar] [CrossRef] [Green Version]
  71. Wiederstein, M.; Sippl, M.J. ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007, 35, W407–W410. [Google Scholar] [CrossRef] [Green Version]
  72. Sippl, M.J. Recognition of errors in three--dimensional structures of proteins. Proteins Struct. Funct. Bioinform. 1993, 17, 355–362. [Google Scholar] [CrossRef]
  73. Colovos, C.; Yeates, T.O. Verification of protein structures: Patterns of nonbonded atomic interactions. Protein Sci. 1993, 2, 1511–1519. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Kaminski, G.A.; Friesner, R.A.; Tirado-Rives, J.; Jorgensen, W.L. Evaluation and reparametrization of the OPLS-AA force field for proteins via comparison with accurate quantum chemical calculations on peptides. J. Phys. Chem. B 2001, 105, 6474–6487. [Google Scholar] [CrossRef]
  75. Berendsen, H.J.C.; Postma, J.P.M.; van Gunsteren, W.F.; Hermans, J. Interaction Models for Water in Relation to Protein Hydration. In Intermolecular Forces; Springer, Dordrecht: Dordrecht, The Netherlands, 1981; pp. 331–342. [Google Scholar]
  76. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 14101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  78. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindah, E. Gromacs: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef] [Green Version]
  79. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  80. Madhavi Sastry, G.; Adzhigirey, M.; Day, T.; Annabhimoju, R.; Sherman, W. Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments. J. Comput. Aided Mol. Des. 2013, 27, 221–234. [Google Scholar] [CrossRef]
  81. Rostkowski, M.; Olsson, M.H.; Søndergaard, C.R.; Jensen, J.H. Graphical analysis of pH-dependent properties of proteins predicted using PROPKA. BMC Struct. Biol. 2011, 11, 6. [Google Scholar] [CrossRef] [Green Version]
  82. Halgren, T. New method for fast and accurate binding-site identification and analysis. Chem. Biol. Drug Des. 2007, 69, 146–148. [Google Scholar] [CrossRef]
  83. Halgren, T.A. Identifying and characterizing binding sites and assessing druggability. J. Chem. Inf. Model. 2009, 49, 377–389. [Google Scholar] [CrossRef]
  84. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Halgren, T.A.; Murphy, R.B.; Friesner, R.A.; Beard, H.S.; Frye, L.L.; Pollard, W.T.; Banks, J.L. Glide: A New Approach for Rapid, Accurate Docking and Scoring. 2. Enrichment Factors in Database Screening. J. Med. Chem. 2004, 47, 1750–1759. [Google Scholar] [CrossRef] [PubMed]
  86. Irwin, J.J.; Shoichet, B.K. ZINC—A free database of commercially available compounds for virtual screening. J. Chem. Inf. Model. 2005, 45, 177–182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Shelley, J.C.; Cholleti, A.; Frye, L.L.; Greenwood, J.R.; Timlin, M.R.; Uchimaya, M. Epik: A software program for pKa prediction and protonation state generation for drug-like molecules. J. Comput. Aided Mol. Des. 2007, 21, 681–691. [Google Scholar] [CrossRef] [PubMed]
  88. Greenwood, J.R.; Calkins, D.; Sullivan, A.P.; Shelley, J.C. Towards the comprehensive, rapid, and accurate prediction of the favorable tautomeric states of drug-like molecules in aqueous solution. J. Comput. Aided Mol. Des. 2010, 24, 591–604. [Google Scholar] [CrossRef]
  89. Vacas, A.; Sugden, C.; Velasco-Rodriguez, Ó.; Algarabel-Olona, M.; Peña-Guerrero, J.; Larrea, E.; Fernández-Rubio, C.; Nguewa, P.A. Construction of Two mCherry Plasmids (pXG-mCherry) for Transgenic Leishmania: Valuable Tools for Future Molecular Analysis. J. Parasitol. Res. 2017, 2017, 1964531. [Google Scholar] [CrossRef] [Green Version]
  90. Medina-Acosta, E.; Cross, G.A.M. Rapid isolation of DNA from trypanosomatid protozoa using a simple “mini-prep” procedure. Mol. Biochem. Parasitol. 1993, 59, 327–329. [Google Scholar] [CrossRef]
  91. Cruz, A.; Coburn, C.M.; Beverley, S.M. Double targeted gene replacement for creating null mutants. Proc. Natl. Acad. Sci. USA 1991, 88, 7170–7174. [Google Scholar] [CrossRef] [Green Version]
  92. Sacks, D.L.; Hieny, S.; Sher, A. Identification of cell surface carbohydrate and antigenic changes between noninfective and infective developmental stages of Leishmania major promastigotes. J. Immunol. 1985, 135, 564–569. [Google Scholar]
  93. Fernández-Rubio, C.; Campbell, D.; Vacas, A.; Ibañez, E.; Moreno, E.; Espuelas, S.; Calvo, A.; Palop, J.A.; Plano, D.; Sanmartin, C.; et al. Leishmanicidal activities of novel methylseleno-imidocarbamates. Antimicrob. Agents Chemother. 2015, 59, 5705–5713. [Google Scholar] [CrossRef] [Green Version]
  94. Dame, N. GraphPad Prism, Version 7.0, Software for Technical Computation; GraphPad Software, Inc.: San Diego, CA, USA, 2015. [Google Scholar]
  95. Wang, J.; Pan, C.; Wang, Y.; Ye, L.; Wu, J.; Chen, L.; Zou, T.; Lu, G. Genome-wide identification of MAPK, MAPKK, and MAPKKK gene families and transcriptional profiling analysis during development and stress response in cucumber. BMC Genomics 2015, 16, 386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Nolen, B.; Taylor, S.; Ghosh, G. Regulation of protein kinases: Controlling activity through activation segment conformation. Mol. Cell 2004, 15, 661–675. [Google Scholar] [CrossRef] [PubMed]
  97. Henikoff, S.; Henikoff, J.G. Amino acid substitution matrices. Adv. Protein Chem. 2000, 54, 73–97. [Google Scholar] [PubMed]
  98. Bansal, A.; Singh, S.; More, K.R.; Hans, D.; Nangalia, K.; Yogavel, M.; Sharma, A.; Chitnis, C.E. Erratum: Characterization of plasmodium falciparum calcium-dependent protein kinase 1 (PfCDPK1) and its role in microneme secretion during erythrocyte invasion. J. Biol. Chem. 2013, 288, 19643. [Google Scholar] [CrossRef] [Green Version]
  99. Zhang, L.; Wang, J.C.; Hou, L.; Cao, P.R.; Wu, L.; Zhang, Q.S.; Yang, H.Y.; Zang, Y.; Ding, J.P.; Li, J. Functional role of histidine in the conserved His-x-Asp motif in the catalytic core of protein kinases. Sci. Rep. 2015, 5, 10115. [Google Scholar] [CrossRef] [Green Version]
  100. Taylor, S.S.; Keshwani, M.M.; Steichen, J.M.; Kornev, A.P. Evolution of the eukaryotic protein kinases as dynamic molecular switches. Philos. Trans. R. Soc. B Biol. Sci. 2012, 367, 2517–2528. [Google Scholar] [CrossRef] [Green Version]
  101. Ren, J.; Wen, L.; Gao, X.; Jin, C.; Xue, Y.; Yao, X. DOG 1.0: Illustrator of protein domain structures. Cell Res. 2009, 19, 271–273. [Google Scholar] [CrossRef]
  102. Negi, A.; Bhushan, S.; Gupta, P.; Garg, P.; Kumar, R. Cystathionine β -Lyase-Like Protein with Pyridoxal Binding Domain Characterized in Leishmania major by Comparative Sequence Analysis and Homology Modelling. ISRN Comput. Biol. 2013, 2013, 1–9. [Google Scholar] [CrossRef]
  103. Taylor, S.S.; Kornev, A.P. Protein kinases: Evolution of dynamic regulatory proteins. Trends Biochem. Sci. 2011, 36, 65–77. [Google Scholar] [CrossRef] [Green Version]
  104. Fabbro, D.; Cowan-Jacob, S.W.; Moebitz, H. Ten things you should know about protein kinases: IUPHAR Review 14. Br. J. Pharmacol. 2015, 172, 2675–2700. [Google Scholar] [CrossRef] [Green Version]
  105. Xu, M.; Yu, L.; Wan, B.; Yu, L.; Huang, Q. Predicting inactive conformations of protein kinases using active structures: Conformational selection of Type-II inhibitors. PLoS ONE 2011, 6, e22644. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Gräter, F.; Shen, J.; Jiang, H.; Gautel, M.; Grubmüller, H. Mechanically induced titin kinase activation studied by force-probe molecular dynamics simulations. Biophys. J. 2005, 88, 790–804. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Roskoski, R. A historical overview of protein kinases and their targeted small molecule inhibitors. Pharmacol. Res. 2015, 100, 1–23. [Google Scholar] [CrossRef] [PubMed]
  108. Roskoski, R. Src protein-tyrosine kinase structure, mechanism, and small molecule inhibitors. Pharmacol. Res. 2015, 94, 9–25. [Google Scholar] [CrossRef] [PubMed]
  109. Oruganty, K.; Talathi, N.S.; Wood, Z.A.; Kannan, N. Identification of a hidden strain switch provides clues to an ancient structural mechanism in protein kinases. Proc. Natl. Acad. Sci. USA 2013, 110, 924–929. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  110. Shalev, M.; Kondo, J.; Kopelyanskiy, D.; Jaffe, C.L.; Adir, N.; Baasov, T. Identification of the molecular attributes required for aminoglycoside activity against Leishmania. Proc. Natl. Acad. Sci. USA 2013, 110, 13333–13338. [Google Scholar] [CrossRef] [Green Version]
  111. Pfister, P.; Hobbie, S.; Brüll, C.; Corti, N.; Vasella, A.; Westhof, E.; Böttger, E.C. Mutagenesis of 16S rRNA C1409-G1491 base-pair differentiates between 6′OH and 6′NH3+ aminoglycosides. J. Mol. Biol. 2005, 346, 467–475. [Google Scholar] [CrossRef]
  112. Davidson, R.N.; den Boer, M.; Ritmeijer, K. Paromomycin. Trans. R. Soc. Trop. Med. Hyg. 2009, 103, 653–660. [Google Scholar] [CrossRef]
  113. El-On, J.; Jacobs, G.P.; Witztum, E.; Greenblatt, C.L. Development of topical treatment for cutaneous leishmaniasis caused by Leishmania major in experimental animals. Antimicrob. Agents Chemother. 1984, 26, 745–751. [Google Scholar] [CrossRef] [Green Version]
  114. Oruganty, K.; Talevich, E.E.; Neuwald, A.F.; Kannan, N. Identification and classification of small molecule kinases: Insights into substrate recognition and specificity Genome evolution and evolutionary systems biology. BMC Evol. Biol. 2016, 16, 7. [Google Scholar] [CrossRef] [Green Version]
  115. Kannan, N.; Neuwald, A.F. Did protein kinase regulatory mechanisms evolve through elaboration of a simple structural component? J. Mol. Biol. 2005, 351, 956–972. [Google Scholar] [CrossRef] [PubMed]
  116. Wu, M.Y.; Dai, D.Q.; Yan, H. PRL-dock: Protein-ligand docking based on hydrogen bond matching and probabilistic relaxation labeling. Proteins Struct. Funct. Bioinform. 2012, 80, 2137–2153. [Google Scholar] [CrossRef] [PubMed]
  117. Shaw, K.J.; Rather, P.N.; Hare, R.S.; Miller, G.H. Molecular genetics of aminoglycoside resistance genes and familial relationships of the aminoglycoside-modifying enzymes. Microbiol. Rev. 1993, 57, 138–163. [Google Scholar] [PubMed]
  118. Thompson, P.R.; Boehr, D.D.; Berghuis, A.M.; Wright, G.D. Mechanism of aminoglycoside antibiotic kinase APH(3′)-IIIa: Role of the nucleotide positioning loop. Biochemistry 2002, 41, 7001–7007. [Google Scholar] [CrossRef] [PubMed]
  119. Menard, R.; Molinas, C.; Arthur, M.; Duval, J.; Courvalin, P.; Leclercq, R. Overproduction of 3′-aminoglycoside phosphotransferase type I confers resistance to tobramycin in Escherichia coli. Antimicrob. Agents Chemother. 1993, 37, 78–83. [Google Scholar] [CrossRef] [Green Version]
  120. Magnet, S.; Smith, T.A.; Zheng, R.; Nordmann, P.; Blanchard, J.S. Aminoglycoside resistance resulting from tight drug binding to an altered aminoglycoside acetyltransferase. Antimicrob. Agents Chemother. 2003, 47, 1577–1583. [Google Scholar] [CrossRef] [Green Version]
  121. Maarouf, M.; De Kouchkovsky, Y.; Brown, S.; Petit, P.X.; Robert-Gero, M. In vivo interference of paromomycin with mitochondrial activity of Leishmania. Exp. Cell Res. 1997, 232, 339–348. [Google Scholar] [CrossRef]
  122. Bhandari, V.; Sundar, S.; Dujardin, J.C.; Salotra, P. Elucidation of cellular mechanisms involved in experimental paromomycin resistance in leishmania donoVani. Antimicrob. Agents Chemother. 2014, 58, 2580–2585. [Google Scholar] [CrossRef] [Green Version]
  123. Chawla, B.; Jhingran, A.; Panigrahi, A.; Stuart, K.D.; Madhubala, R. Paromomycin affects translation and vesicle-mediated trafficking as revealed by proteomics of paromomycin -susceptible -resistant leishmania donovani. PLoS ONE 2011, 6, e26660. [Google Scholar] [CrossRef]
  124. El Fadili, K.; Messier, N.; Leprohon, P.; Roy, G.; Guimond, C.; Trudel, N.; Saravia, N.G.; Papadopoulou, B.; Légaré, D.; Ouellette, M. Role of the ABC transporter MRPA (PGPA) in antimony resistance in Leishmania infantum axenic and intracellular amastigotes. Antimicrob. Agents Chemother. 2005, 49, 1988–1993. [Google Scholar] [CrossRef] [Green Version]
  125. Khanra, S.; Sarraf, N.R.; Das, A.K.; Roy, S.; Manna, M. Miltefosine Resistant Field Isolate from Indian Kala-Azar Patient Shows Similar Phenotype in Experimental Infection. Sci. Rep. 2017, 7, 10330. [Google Scholar] [CrossRef] [PubMed]
  126. Leprohon, P.; Légaré, D.; Girard, I.; Papadopoulou, B.; Ouellette, M. Modulation of Leishmania ABC protein gene expression through life stages and among drug-resistant parasites. Eukaryot. Cell 2006, 5, 1713–1725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  127. Obonaga, R.; Fernández, O.L.; Valderrama, L.; Rubiano, L.C.; Castro, M.D.M.; Barrera, M.C.; Gomez, M.A.; Saravia, N.G. Treatment failure and miltefosine susceptibility in dermal leishmaniasis caused by leishmania subgenus viannia species. Antimicrob. Agents Chemother. 2014, 58, 144–152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  128. Teklemariam, S.; Hiwot, A.G.; Frommel, D.; Mike, T.L.; Ganlov, G.; Bryceson, A. Aminosidine and its combination with sodium stibogluconate in the treatment of diffuse cutaneous leishmaniasis caused by leishmania aethiopica. Trans. R. Soc. Trop. Med. Hyg. 1994, 88, 334–339. [Google Scholar] [CrossRef]
Figure 1. Homologues of LmjF.22.0810 (LmJean3) kinase are highly conserved in trypanosomatids. (A) Schematic representation of the predicted catalytically relevant regions and domains of LmJean3; (B) Leishmania major, Leishmania mexicana, Leishmania braziliensis, Leishmania infantum, Trypanosoma cruzi, and Trypanosoma brucei Jean3 orthologs sequences alignment. Residues are colored on the basis of their similarity under a BLOSUM62 score matrix [97]; (C) phylogenetic analysis of Jean3. Saccharomyces cerevisiae CDC28 was used to root the tree. Protein sequences from Jean3 trypanosomatid orthologues (accession numbers: XP_001465592, XP_003860813, XP_010699073, XP_001564994, XP_001683232, XP_015664142, KPI87780, EPY38836, CCC91397, XP_845898, CCD21423, XP_819355, XP_009309601, CUI14844, EPY19901) and kinases from Arabidopsis thaliana (CCH26589, NP_564353), Homo sapiens (AAA64745, NP_001156768, NP_443073), S. cerevisiae (NP_010765, NP_009718) and Plasmodium falciparum (P62343) were used to generate the phylogram.
Figure 1. Homologues of LmjF.22.0810 (LmJean3) kinase are highly conserved in trypanosomatids. (A) Schematic representation of the predicted catalytically relevant regions and domains of LmJean3; (B) Leishmania major, Leishmania mexicana, Leishmania braziliensis, Leishmania infantum, Trypanosoma cruzi, and Trypanosoma brucei Jean3 orthologs sequences alignment. Residues are colored on the basis of their similarity under a BLOSUM62 score matrix [97]; (C) phylogenetic analysis of Jean3. Saccharomyces cerevisiae CDC28 was used to root the tree. Protein sequences from Jean3 trypanosomatid orthologues (accession numbers: XP_001465592, XP_003860813, XP_010699073, XP_001564994, XP_001683232, XP_015664142, KPI87780, EPY38836, CCC91397, XP_845898, CCD21423, XP_819355, XP_009309601, CUI14844, EPY19901) and kinases from Arabidopsis thaliana (CCH26589, NP_564353), Homo sapiens (AAA64745, NP_001156768, NP_443073), S. cerevisiae (NP_010765, NP_009718) and Plasmodium falciparum (P62343) were used to generate the phylogram.
Biomolecules 09 00723 g001
Figure 2. LmJean3 was localized in the cytoplasm of L. major parasites, and its gene expression level s doubled in the stationary phase. (A) Visualization of L. major promastigotes expressing GFP2+-LmJean3 (left lane) or LmJean3-mCherry (right lane) fusion proteins; (B) relative expression quantification of LmJean3 from log-phase, procyclic, and metacyclic (infective) promastigotes and from intracellular amastigotes. Bars represent gene expression mean fold change (± SD) from three independent experiments (* p < 0.05; ** p < 0.01).
Figure 2. LmJean3 was localized in the cytoplasm of L. major parasites, and its gene expression level s doubled in the stationary phase. (A) Visualization of L. major promastigotes expressing GFP2+-LmJean3 (left lane) or LmJean3-mCherry (right lane) fusion proteins; (B) relative expression quantification of LmJean3 from log-phase, procyclic, and metacyclic (infective) promastigotes and from intracellular amastigotes. Bars represent gene expression mean fold change (± SD) from three independent experiments (* p < 0.05; ** p < 0.01).
Biomolecules 09 00723 g002
Figure 3. The kinase domain of LmJean3 homologues shares a high overall similarity to CIPK24/SOS2, CIPK23, and Snf1 kinase domains. On top: schematic representation of LmJean3 secondary structure, sequence, and regions. Underneath: LmJean3 sequence alignment (residues 88-211) with Leishmania donovani Jean3 (LdJean3), L. braziliensis Jean3 (LbJean3), A. thaliana CIPK24/SOS2, A. thaliana CIPK23, and S. cerevisiae Snf1. Within the ATP-binding region signature (yellow rectangle), the glycine-rich loop is highlighted as an orange rectangle (residues 15–20). The catalytic and activation loops are marked as blue and green rectangles, respectively. Residues S157 and S162 from the activation loop are highlighted. Cyan and red squares denote the putative catalytically relevant residues: K37, E53, R128, D129, and D147. Alignment was performed using the ClustalW iterative algorithm [53]. The secondary structure of LmJean3 was predicted using the EMBOSS v6.5.7 tool [37] implemented in Geneious v9.1.7 [36]. The figure was designed using IBS v1.02 illustrator [101].
Figure 3. The kinase domain of LmJean3 homologues shares a high overall similarity to CIPK24/SOS2, CIPK23, and Snf1 kinase domains. On top: schematic representation of LmJean3 secondary structure, sequence, and regions. Underneath: LmJean3 sequence alignment (residues 88-211) with Leishmania donovani Jean3 (LdJean3), L. braziliensis Jean3 (LbJean3), A. thaliana CIPK24/SOS2, A. thaliana CIPK23, and S. cerevisiae Snf1. Within the ATP-binding region signature (yellow rectangle), the glycine-rich loop is highlighted as an orange rectangle (residues 15–20). The catalytic and activation loops are marked as blue and green rectangles, respectively. Residues S157 and S162 from the activation loop are highlighted. Cyan and red squares denote the putative catalytically relevant residues: K37, E53, R128, D129, and D147. Alignment was performed using the ClustalW iterative algorithm [53]. The secondary structure of LmJean3 was predicted using the EMBOSS v6.5.7 tool [37] implemented in Geneious v9.1.7 [36]. The figure was designed using IBS v1.02 illustrator [101].
Biomolecules 09 00723 g003
Figure 4. LmJean3 protein structure analysis and evaluation. Statistics of non-bonded interactions analyzed by ERRAT for (A) LmJean3 and (B) the template 4CZT_B. * On the error axis, two lines are drawn to indicate the confidence with which it is possible to reject regions that exceed that error value. ** Expressed as the percentage of the protein for which the calculated error value falls below the 95% rejection limit. (C,D) Comparative protein structure analysis (ProSA) Z-score evaluation for (C) LmJean3 and (D) 4CZT_B; (E) Ramachandran plot from the optimized structure for LmJean3. The homology model displayed 88.4% of the residues in favorable positions, and none was found in the disallowed regions.
Figure 4. LmJean3 protein structure analysis and evaluation. Statistics of non-bonded interactions analyzed by ERRAT for (A) LmJean3 and (B) the template 4CZT_B. * On the error axis, two lines are drawn to indicate the confidence with which it is possible to reject regions that exceed that error value. ** Expressed as the percentage of the protein for which the calculated error value falls below the 95% rejection limit. (C,D) Comparative protein structure analysis (ProSA) Z-score evaluation for (C) LmJean3 and (D) 4CZT_B; (E) Ramachandran plot from the optimized structure for LmJean3. The homology model displayed 88.4% of the residues in favorable positions, and none was found in the disallowed regions.
Biomolecules 09 00723 g004
Figure 5. LmJean3 structure evaluation. (A) Superposition of the optimized LmJean3 homology model (pink N-lobe and cyan C-lobe) and 4CZT_B (purple) revealed LmJean3 catalytic residues K37, E53, R128, D129, and D147 aligned to their equivalents in the template (4CZT_B). The calculated RMSD was 1.87 Å; (B) predicted homology model for LmJean3 showing the ATP-binding region in yellow, glycine-rich loop in orange, activation loop in green, catalytic loop in blue, and catalytic residues R128 and D129 in red; (C) molecular surface representation of the catalytic and regulatory spines of LmJean3. The residues forming the spines are displayed as white sticks; (D) structure alignment of LmJean3 catalytically relevant residues (K37, E53, R128, D129, and D147) at the initial point of the simulation (magenta N-lobe and orange C-lobe) and the average calculated structure from the trajectory (pink N-lobe and cyan C-lobe). The catalytic residues from the initial frame are colored yellow, and their equivalents from the average structure are colored white. The calculated RMSD was 1.15 Å; (E) salt bridges (SB) and hydrogen bonds (HB), with persistence higher than 40% during the trajectory of the molecular dynamics’ simulation are shown in the average structure.
Figure 5. LmJean3 structure evaluation. (A) Superposition of the optimized LmJean3 homology model (pink N-lobe and cyan C-lobe) and 4CZT_B (purple) revealed LmJean3 catalytic residues K37, E53, R128, D129, and D147 aligned to their equivalents in the template (4CZT_B). The calculated RMSD was 1.87 Å; (B) predicted homology model for LmJean3 showing the ATP-binding region in yellow, glycine-rich loop in orange, activation loop in green, catalytic loop in blue, and catalytic residues R128 and D129 in red; (C) molecular surface representation of the catalytic and regulatory spines of LmJean3. The residues forming the spines are displayed as white sticks; (D) structure alignment of LmJean3 catalytically relevant residues (K37, E53, R128, D129, and D147) at the initial point of the simulation (magenta N-lobe and orange C-lobe) and the average calculated structure from the trajectory (pink N-lobe and cyan C-lobe). The catalytic residues from the initial frame are colored yellow, and their equivalents from the average structure are colored white. The calculated RMSD was 1.15 Å; (E) salt bridges (SB) and hydrogen bonds (HB), with persistence higher than 40% during the trajectory of the molecular dynamics’ simulation are shown in the average structure.
Biomolecules 09 00723 g005
Figure 6. LmJean3 displays two promising binding sites for small-molecule, drug-like ligands. Hydrophobic (yellow), donor (blue), and acceptor (red) maps for (A) the ATP-binding pocket and the adjacent allosteric site and (B) the protein substrate binding site. SiteMap version 3.8 [82,83] was used to identify the binding sites, and Maestro Schrödinger Release 2017-2 to compose the figures.
Figure 6. LmJean3 displays two promising binding sites for small-molecule, drug-like ligands. Hydrophobic (yellow), donor (blue), and acceptor (red) maps for (A) the ATP-binding pocket and the adjacent allosteric site and (B) the protein substrate binding site. SiteMap version 3.8 [82,83] was used to identify the binding sites, and Maestro Schrödinger Release 2017-2 to compose the figures.
Biomolecules 09 00723 g006
Figure 7. Representative aminoglycosides (AG) antibiotics of 4,5- and 4,6-disubstituted subclasses. Structure representation of (A) 4,5-disubsituted AGs (paromomycin and neomycin B) and (B) 4,6-disubstituted AGs (amikacin, tobramycin, and geneticin); HABA: α-hydroxy-γ-aminobutyric acid. The position 6′ from Ring I is highlighted by a green-colored circle; (C) predicted interactions of paromomycin with LmJean3 Site B. The 3D representation of the docking pose is shown on the left (distances up to 2.8 Å), and the ligand interaction diagram (distances up to 4 Å) on the right. Species in the ligand interaction diagram are colored by type: hydrophobic residues green, charged residues red (−) or violet (+), polar residues light blue and glycine light green. H-bond interactions with amino acids are highlighted as purple lines, salt bridges as blue- and red-colored lines. The figure was prepared using Maestro Schrödinger Release 2017-2.
Figure 7. Representative aminoglycosides (AG) antibiotics of 4,5- and 4,6-disubstituted subclasses. Structure representation of (A) 4,5-disubsituted AGs (paromomycin and neomycin B) and (B) 4,6-disubstituted AGs (amikacin, tobramycin, and geneticin); HABA: α-hydroxy-γ-aminobutyric acid. The position 6′ from Ring I is highlighted by a green-colored circle; (C) predicted interactions of paromomycin with LmJean3 Site B. The 3D representation of the docking pose is shown on the left (distances up to 2.8 Å), and the ligand interaction diagram (distances up to 4 Å) on the right. Species in the ligand interaction diagram are colored by type: hydrophobic residues green, charged residues red (−) or violet (+), polar residues light blue and glycine light green. H-bond interactions with amino acids are highlighted as purple lines, salt bridges as blue- and red-colored lines. The figure was prepared using Maestro Schrödinger Release 2017-2.
Biomolecules 09 00723 g007
Figure 8. Site A and Site B docking scores. The top-ranked molecules resulting from docking to Site B were also employed to predict their docking scores to Site A. Additionally, CHAPS (present in the co-crystal structure) and ATP (co-factor) were also included in the analysis. The predicted GlideScores (kcal/mol) for each molecule are represented on the X-axis for Site A and on the Y-axis for Site B.
Figure 8. Site A and Site B docking scores. The top-ranked molecules resulting from docking to Site B were also employed to predict their docking scores to Site A. Additionally, CHAPS (present in the co-crystal structure) and ATP (co-factor) were also included in the analysis. The predicted GlideScores (kcal/mol) for each molecule are represented on the X-axis for Site A and on the Y-axis for Site B.
Biomolecules 09 00723 g008
Figure 9. LmJean3-overexpression in Leishmania (LmJ3OE parasites) reduced the sensitivity to AG treatment after 48 h and 72 h. Leishmanicidal activity determination (EC50) by MTT assays of (A) paromomycin, (B) geneticin (G418), (C) amphotericin B, and (D) miltefosine in (mock) control and LmJ3OE parasites. Bars represent the normalized mean EC50 (± SD) from three independent experiments (* p < 0.05, ** p < 0.01, *** p < 0.001).
Figure 9. LmJean3-overexpression in Leishmania (LmJ3OE parasites) reduced the sensitivity to AG treatment after 48 h and 72 h. Leishmanicidal activity determination (EC50) by MTT assays of (A) paromomycin, (B) geneticin (G418), (C) amphotericin B, and (D) miltefosine in (mock) control and LmJ3OE parasites. Bars represent the normalized mean EC50 (± SD) from three independent experiments (* p < 0.05, ** p < 0.01, *** p < 0.001).
Biomolecules 09 00723 g009
Table 1. Persistence of the SB and HB found within the kinase domain of LmJean3 during the molecular dynamics’ simulation (20 ns). See Tables S2 and S3 for an exhaustive list of all the possible interactions found during the trajectory.
Table 1. Persistence of the SB and HB found within the kinase domain of LmJean3 during the molecular dynamics’ simulation (20 ns). See Tables S2 and S3 for an exhaustive list of all the possible interactions found during the trajectory.
DonorAcceptorBond TypePersistence (>40%)
ARG254GLU180SB100.00
LYS131ASP129SB98.45
ARG152GLU53SB85.66
LYS37ASP147SB82.46
ARG128-MainASP192-SideHB60.49
LEU137-MainGLY88-MainHB48.95
ASN134-SideASP129-MainHB40.61
Table 2. Pharmacokinetic parameters important for good oral bioavailability and docking results.
Table 2. Pharmacokinetic parameters important for good oral bioavailability and docking results.
RankIUPAC NameChemical NameMWLogPHBA (Lipinski)HBD (Lipinski)Violations (Lipinski)Docking Score (kcal/mol)
1(2S,3S,4R,5R,6R)-5-amino-2-(aminomethyl)-6-[(2R,3S,4R,5S)-5-[(1R,2R,3S,5R,6S)-3,5-diamino-2-[(2S,3R,4R,5S,6R)-3-amino-4,5-dihydroxy-6-(hydroxymethyl)oxan-2-yl]oxy-6-hydroxycyclohexyl]oxy-4-hydroxy-2-(hydroxymethyl)oxolan-3-yl]oxyoxane-3,4-diolParomomycin sulfate615.6−8.6719182−11.46
2(2R,3S,4R,5R,6R)-5-amino-2-(aminomethyl)-6-[(1R,2R,3S,4R,6S)-4,6-diamino-2-[(2S,3R,4S,5R)-4-[(2R,3R,4R,5S,6S)-3-amino-6-(aminomethyl)-4,5-dihydroxyoxan-2-yl]oxy-3-hydroxy-5-(hydroxymethyl)oxolan-2-yl]oxy-3-hydroxycyclohexyl]oxyoxane-3,4-diolNeomycin sulfate stereoisomer A614.6−8.9619192−10.93
32-(3,4-dihydroxyphenyl)-5,7-dihydroxy-3-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl]oxan-2-yl]oxychromen-4-oneRutin610.5−1.1616102−10.84
41,4-dihydroxy-5,8-bis[2-(2-hydroxyethylamino)ethylamino]anthracene-9,10-dioneMitoxantrone444.50.071081−9.14
5(2S,3R,4S,5S,6R)-4-amino-2-[(1S,2S,3R,4S,6R)-4,6-diamino-3-[(2R,3R,5S,6R)-3-amino-6-(aminomethyl)-5-hydroxyoxan-2-yl]oxy-2-hydroxycyclohexyl]oxy-6-(hydroxymethyl)oxane-3,5-diolTobramycin sulfate467.5−6.8614152−9.06
65-hydroxy-2-(3-hydroxy-4-methoxyphenyl)-7-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl]oxan-2-yl]oxychromen-4-oneDiosmin608.5−0.41582−9.05
7(2S)-4-amino-N-[(1R,2S,3S,4R,5S)-5-amino-2-[(2S,3R,4S,5S,6R)-4-amino-3,5-dihydroxy-6-(hydroxymethyl)oxan-2-yl]oxy-4-[(2R,3R,4S,5S,6R)-6-(aminomethyl)-3,4,5-trihydroxyoxan-2-yl]oxy-3-hydroxycyclohexyl]-2-hydroxybutanamideAmikacin stereoisomer585.6−8.4318172−9.04
MW, molecular weight; LogP, logarithm of compound partition coefficient between n-octanol and water; HBA, number of hydrogen bond acceptors; HBD, number of hydrogen bond donors.

Share and Cite

MDPI and ACS Style

Vacas, A.; Fernández-Rubio, C.; Algarabel, M.; Peña-Guerrero, J.; Larrea, E.; Rocha Formiga, F.; García-Sosa, A.T.; Nguewa, P.A. The Novel Serine/Threonine Protein Kinase LmjF.22.0810 from Leishmania major May Be Involved in the Resistance to Drugs such as Paromomycin. Biomolecules 2019, 9, 723. https://doi.org/10.3390/biom9110723

AMA Style

Vacas A, Fernández-Rubio C, Algarabel M, Peña-Guerrero J, Larrea E, Rocha Formiga F, García-Sosa AT, Nguewa PA. The Novel Serine/Threonine Protein Kinase LmjF.22.0810 from Leishmania major May Be Involved in the Resistance to Drugs such as Paromomycin. Biomolecules. 2019; 9(11):723. https://doi.org/10.3390/biom9110723

Chicago/Turabian Style

Vacas, Andrés, Celia Fernández-Rubio, Miriam Algarabel, José Peña-Guerrero, Esther Larrea, Fabio Rocha Formiga, Alfonso T. García-Sosa, and Paul A. Nguewa. 2019. "The Novel Serine/Threonine Protein Kinase LmjF.22.0810 from Leishmania major May Be Involved in the Resistance to Drugs such as Paromomycin" Biomolecules 9, no. 11: 723. https://doi.org/10.3390/biom9110723

APA Style

Vacas, A., Fernández-Rubio, C., Algarabel, M., Peña-Guerrero, J., Larrea, E., Rocha Formiga, F., García-Sosa, A. T., & Nguewa, P. A. (2019). The Novel Serine/Threonine Protein Kinase LmjF.22.0810 from Leishmania major May Be Involved in the Resistance to Drugs such as Paromomycin. Biomolecules, 9(11), 723. https://doi.org/10.3390/biom9110723

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop