Molecular insight into endosulfan degradation by Ese protein from Arthrobacter: Evidence-based structural bioinformatics and quantum mechanical calculations.

Endosulfan is an organochlorine insecticide widely used for agricultural pest control. Many nations worldwide have restricted or completely banned it due to its extreme toxicity to fish and aquatic invertebrates. Arthrobacter sp. strain KW has the ability to degrade α, ß endosulfan and its intermediate metabolite endosulfate; this degradation is associated with Ese protein, a two-component flavin-dependent monooxygenase (TC-FDM). Employing in silico tools, we obtained the 3D model of Ese protein, and our results suggest that it belongs to the Luciferase Like Monooxygenase family (LLM). Docking studies showed that the residues V59, V315, D316, and T335 interact with α-endosulfan. The residues: V59, T60, V315, D316, and T335 are implicated in the interacting site with ß-endosulfan, and the residues: H17, V315, D316, T335, N364, and Q363 participate in the interaction with endosulfate. Topological analysis of the electron density by means of the Quantum Theory of Atoms in Molecules (QTAIM) and the Non-Covalent Interaction (NCI) index reveals that the Ese-ligands complexes are formed mainly by dispersive forces, where Cl atoms have a predominant role. As Ese is a monooxygenase member, we predict the homodimer formation. However, enzymatic studies must be developed to investigate the Ese protein's enzymatic and catalytic activity.

Design, synthesis, in silico, and in vitro evaluation of benzylbenzimidazolone derivatives as potential drugs on α-glucosidase and glucokinase as pharmacological targets.

Benzimidazolones have shown biological activities, including antihyperglycemic and hypoglycemic, by inhibiting or activating of α-glu and GK. The aim of this study is the rational design of compounds using in silico assays to delimitate the selection of structures to synthesize and the in vitro evaluation of benzimidazolone derivatives in blood glucose control. A docking of 23 benzimidazolone derivatives was performed; selecting the compounds with better in silico profiles to synthesize by microwave-irradiation/conventional heat and evaluate in enzymatic in vitro evaluation. Compounds 2k, 2m, 2r, and 2s presented the best in silico profiles, showing good affinity energy (-10.9 to -8.6 kcal mol-1) and binding with catalytic-amino acids. They were synthesized at 70 °C and 24 h using DMF as the solvent and potassium carbonate (yield: 22-38%). The results with α-glu showed moderate inhibition of 2k (14 ± 1.23-29 ± 0.45), 2m (12 ± 2.21-36 ± 0.30), 2r (7 ± 2.21-13 ± 1.34), and 2s (11 ± 0.74-35 ± 2.95) at evaluated concentrations (0.1 to 100 µg mL-1). The GK activation assay showed an enzymatic activity increase; compound 2k increased 1.31 and 2.83 more than normal activity, 2m (2.13-fold), 2s (2.86 and 3.74-fold) at 100 and 200 µg mL-1 respectively. The present study showed that the 2s derivative presents moderate potential as an α-glu inhibitor and a good activator potential of GK, suggesting that this compound is a good candidate for blood glucose control through antihyperglycemic and hypoglycemic mechanisms.

Dynamic Association of ESCRT-II Proteins with ESCRT-I and ESCRT-III Complexes during Phagocytosis of Entamoeba histolytica.

By their active movement and voraux phagocytosis, the trophozoites of Entamoeba histolytica constitute an excellent system to investigate the dynamics of the Endosomal Sorting Complex Required for Transport (ESCRT) protein interactions through phagocytosis. Here, we studied the proteins forming the E. histolytica ESCRT-II complex and their relationship with other phagocytosis-involved molecules. Bioinformatics analysis predicted that EhVps22, EhVps25, and EhVps36 are E. histolytica bona fide orthologues of the ESCRT-II protein families. Recombinant proteins and specific antibodies revealed that ESCRT-II proteins interact with each other, with other ESCRT proteins, and phagocytosis-involved molecules, such as the adhesin (EhADH). Laser confocal microscopy, pull-down assays, and mass spectrometry analysis disclosed that during phagocytosis, ESCRT-II accompanies the red blood cells (RBCs) from their attachment to the trophozoites until their arrival to multivesicular bodies (MVBs), changing their interactive patterns according to the time and place of the process. Knocked-down trophozoites in the Ehvps25 gene presented a 50% lower rate of phagocytosis than the controls and lower efficiency to adhere RBCs. In conclusion, ESCRT-II interacts with other molecules during prey contact and conduction throughout the phagocytic channel and trophozoites membranous system. ESCRT-II proteins are members of the protein chain during vesicle trafficking and are fundamental for the continuity and efficiency of phagocytosis.

Alpha-terthienyl increases filamentous actin of Entamoeba histolytica.

This study aimed to know if alpha terthienyl (α-T) affects E. histolytica viability and to analyze its effect on the actin cytoskeleton. Trophozoites of E. histolytica HM1-IMSS were treated with α-T, then, cell viability and morphology were evaluated using tetrazolium salts and scanning electron microscopy, respectively; while actin filaments (F-actin) were stained with rhodamine-phalloidin, observed by confocal microscopy and quantified by fluorometry. Data showed that α-T inhibited cell viability of trophozoites (IC50, 19.43 µg / mL), affected the cell morphology, and increased the F-actin in a dose-dependent manner. Production of reactive oxygen species and RhoA-GTP levels remained normal in α-T-treated amebas. Two inhibitors that affect the organization of the trophozoites cytoskeleton, one that interacts directly with actin, Cytochalasin D (CD), and one that affects the Rho signaling pathway by inhibiting the downstream effector Rock, Y27632, were tested. Y27632 did not affect the increase of polymerized actin observed with α-T, this compound partially ameliorates the potent disrupting effects of CD on actin filaments. Docking results suggest that α-T could be an antagonist of CD for the same interaction zone in actin, however, more studies are needed to define the action mechanism of this compound.

Intra-host genetic population diversity: Role in emergence and persistence of drug resistance among Mycobacterium tuberculosis complex minor variants.

Drug resistant tuberculosis (DR-TB) is an important public health issue in different parts of the world. Mycobacterium tuberculosis complex variants (MTBC vars) preferentially infect certain hosts, limiting their distribution to different ecosystems. However, MTBC vars can infect other hosts beyond their preferred target potentially contributing to persistence of drug resistance (DR) in other niches. Here, we performed a comprehensive intra-host genetic analysis for the identification of DR-related mutations among all MTBC minor vars whole genome sequences (8,095 strains) publicly available worldwide. High confidence drug-resistance mutations in katG (isoniazid), rpsL (streptomycin), pncA (pyrazinamide), rpoB (rifampicin) and gyrA (fluoroquinolones) genes were identified among intrahost minor sub-populations in 197 different strains (2.43%) belonging to vars africanum, bovis, caprae, microti, orygis and pinnipedii. In addition, a three-dimensional structure modeling analysis to assess the role of novel mutations was also performed. Our findings highlight the importance of detecting discrete intra-host populations carrying DR mutations.

Molecular insights into the interaction of HPV-16 E6 variants against MAGI-1 PDZ1 domain.

Oncogenic protein E6 from Human Papilloma Virus 16 (HPV-16) mediates the degradation of Membrane-associated guanylate kinase with inverted domain structure-1 (MAGI-1), throughout the interaction of its protein binding motif (PBM) with the Discs-large homologous regions 1 (PDZ1) domain of MAG1-1. Generic variation in the E6 gene that translates to changes in the protein's amino acidic sequence modifies the interaction of E6 with the cellular protein MAGI-1. MAGI-1 is a scaffolding protein found at tight junctions of epithelial cells, where it interacts with a variety of proteins regulating signaling pathways. MAGI-1 is a multidomain protein containing two WW (rsp-domain-9), one guanylate kinase-like, and six PDZ domains. PDZ domains played an important role in the function of MAGI-1 and served as targets for several viral proteins including the HPV-16 E6. The aim of this work was to evaluate, with an in silico approach, employing molecular dynamics simulation and protein-protein docking, the interaction of the intragenic variants E-G350 (L83V), E-C188/G350 (E29Q/L83V), E-A176/G350 (D25N/L83V), E6-AAa (Q14H/H78Y/83V) y E6-AAc (Q14H/I27RH78Y/L83V) and E6-reference of HPV-16 with MAGI-1. We found that variants E-G350, E-C188/G350, E-A176/G350, AAa and AAc increase their affinity to our two models of MAGI-1 compared to E6-reference.

microRNA Profile Associated with Positive Lymph Node Metastasis in Early-Stage Cervical Cancer.

Lymph node metastasis (LNM) is an important prognostic factor in cervical cancer (CC). In early stages, the risk of LNM is approximately 3.7 to 21.7%, and the 5-year overall survival decreases from 80% to 53% when metastatic disease is identified in the lymph nodes. Few reports have analyzed the relationship between miRNA expression and the presence of LNM. The aim of this study was to identify a subset of miRNAs related to LNM in early-stage CC patients. Formalin-fixed paraffin-embedded tissue blocks were collected from patients with early-stage CC treated by radical hysterectomy with lymphadenectomy. We analyzed samples from two groups of patients-one group with LNM and the other without LNM. Global miRNA expression was identified by microarray analysis, and cluster analysis was used to determine a subset of miRNAs associated with LNM. Microarray expression profiling identified a subset of 36 differentially expressed miRNAs in the two groups (fold change (FC) ≥ 1.5 and p < 0.01). We validated the expression of seven miRNAs; miR-487b, miR-29b-2-5p, and miR-195 were underexpressed, and miR-92b-5p, miR-483-5p, miR-4534, and miR-548ac were overexpressed according to the microarray experiments. This signature exhibited prognostic value for identifying early-stage CC patients with LNM. These findings may help detect LNM that cannot be observed in imaging studies.

EhVps23: A Component of ESCRT-I That Participates in Vesicular Trafficking and Phagocytosis of Entamoeba histolytica.

The endosomal sorting complex required for transport (ESCRT) is formed by ESCRT-0, ESCRT-I, ESCRT-II, ESCRT-III complexes, and accessory proteins. It conducts vesicular trafficking in eukaryotes through the formation of vesicles and membrane fission and fusion events. The trophozoites of Entamoeba histolytica, the protozoan responsible for human amoebiasis, presents an active membrane movement in basal state that increases during phagocytosis and tissue invasion. ESCRT-III complex has a pivotal role during these events, but ESCRT-0, ESCRT-I and ESCRT-II have been poorly studied. Here, we unveiled the E. histolytica ESCRT-I complex and its implication in vesicular trafficking and phagocytosis, as well as the molecular relationships with other phagocytosis-involved molecules. We found a gene encoding for a putative EhVps23 protein with the ubiquitin-binding and Vps23 core domains. In basal state, it was in the plasma membrane, cytoplasmic vesicles and multivesicular bodies, whereas during phagocytosis it was extensively ubiquitinated and detected in phagosomes and connected vesicles. Docking analysis, immunoprecipitation assays and microscopy studies evidenced its interaction with EhUbiquitin, EhADH, EhVps32 proteins, and the lysobisphosphatidic acid phospholipid. The knocking down of the Ehvps23 gene resulted in lower rates of phagocytosis. Our results disclosed the concert of finely regulated molecules and vesicular structures participating in vesicular trafficking-related events with a pivotal role of EhVps23.

Protein Sumoylation Is Crucial for Phagocytosis in Entamoeba histolytica Trophozoites.

Posttranslational modifications provide Entamoeba histolytica proteins the timing and signaling to intervene during different processes, such as phagocytosis. However, SUMOylation has not been studied in E. histolytica yet. Here, we characterized the E. histolytica SUMO gene, its product (EhSUMO), and the relevance of SUMOylation in phagocytosis. Our results indicated that EhSUMO has an extended N-terminus that differentiates SUMO from ubiquitin. It also presents the GG residues at the C-terminus and the ΨKXE/D binding motif, both involved in target protein contact. Additionally, the E. histolytica genome possesses the enzymes belonging to the SUMOylation-deSUMOylation machinery. Confocal microscopy assays disclosed a remarkable EhSUMO membrane activity with convoluted and changing structures in trophozoites during erythrophagocytosis. SUMOylated proteins appeared in pseudopodia, phagocytic channels, and around the adhered and ingested erythrocytes. Docking analysis predicted interaction of EhSUMO with EhADH (an ALIX family protein), and immunoprecipitation and immunofluorescence assays revealed that the association increased during phagocytosis; whereas the EhVps32 (a protein of the ESCRT-III complex)-EhSUMO interaction appeared stronger since basal conditions. In EhSUMO knocked-down trophozoites, the bizarre membranous structures disappeared, and EhSUMO interaction with EhADH and EhVps32 diminished. Our results evidenced the presence of a SUMO gene in E. histolytica and the SUMOylation relevance during phagocytosis. This is supported by bioinformatics screening of many other proteins of E. histolytica involved in phagocytosis, which present putative SUMOylation sites and the ΨKXE/D binding motif.

Structural insights into SARS-CoV-2 spike protein and its natural mutants found in Mexican population.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a newly emerged coronavirus responsible for coronavirus disease 2019 (COVID-19); it become a pandemic since March 2020. To date, there have been described three lineages of SARS-CoV-2 circulating worldwide, two of them are found among Mexican population, within these, we observed three mutations of spike (S) protein located at amino acids H49Y, D614G, and T573I. To understand if these mutations could affect the structural behavior of S protein of SARS-CoV-2, as well as the binding with S protein inhibitors (cepharanthine, nelfinavir, and hydroxychloroquine), molecular dynamic simulations and molecular docking were employed. It was found that these punctual mutations affect considerably the structural behavior of the S protein compared to wild type, which also affect the binding of its inhibitors into their respective binding site. Thus, further experimental studies are needed to explore if these affectations have an impact on drug-S protein binding and its possible clinical effect.

Modeling and Molecular Dynamics of the 3D Structure of the HPV16 E7 Protein and Its Variants.

The oncogenic potential of high-risk human papillomavirus (HPV) is predicated on the production of the E6 and E7 oncoproteins, which are responsible for disrupting the control of the cell cycle. Epidemiological studies have proposed that the presence of the N29S and H51N variants of the HPV16 E7 protein is significantly associated with cervical cancer. It has been suggested that changes in the amino acid sequence of E7 variants may affect the oncoprotein 3D structure; however, this remains uncertain. An analysis of the structural differences of the HPV16 E7 protein and its variants (N29S and H51N) was performed through homology modeling and structural refinement by molecular dynamics simulation. We propose, for the first time, a 3D structure of the E7 reference protein and two of Its variants (N29S and H51N), and conclude that the mutations induced by the variants in N29S and H51N have a significant influence on the 3D structure of the E7 protein of HPV16, which could be related to the oncogenic capacity of this protein.

Analysis of the Relationship Between Metalloprotease-9 and Tau Protein in Alzheimer's Disease.

BACKGROUND: Neurofibrillary tangles (NFTs) and amyloid plaques are the neuropathological hallmarks in brains with Alzheimer's disease (AD). Post-translational modifications of tau, such as phosphorylation and truncation, have been proposed as initiators in the assembly of the abnormal paired helical filaments that constitute the NFTs. Neurons and NFTs are sites of matrix metalloproteinases (MMPs). OBJECTIVE: The aim of this study was to analyze the relationship of MMP-9 and tau protein in brain samples with AD. METHODS: This study was performed on brain tissue samples from patients with early, moderate, and late AD. MMPs and tau levels were analyzed by western blot and gelatin-substrate zymography. Immunofluorescence techniques and confocal microscopy were used to analyze the presence of both proteins in NFTs. Further, molecular dynamics simulations (MDS) and protein-protein docking were conducted to predict interaction between MMP-9 and tau protein. RESULTS: MMP-9 expression was greatest in moderate and late AD, whereas MMP-2 expression was only increased in late-stage AD. Interestingly, confocal microscopy revealed NFTs in which there was co-localization of MMP-9 and tau protein. MDS and protein-protein docking predictions indicate that a high-affinity complex can be formed between MMP-9 and full-length tau protein. CONCLUSION: These observations provide preliminary evidence of an interaction between these two proteins. Post-translational modifications of tau protein, such as C-terminal truncation or phosphorylation of amino acid residues in the MMP-9 recognition site and conformational changes in the protein, such as folding of the N-terminal sequence over the three-repeat domain, could preclude the interaction between MMP-9 and tau protein during stages of NFT development.

High prevalence of autosomal recessive congenital ichthyosis in a Mexican population caused by a new mutation in the TGM1 gene: epidemiological evidence of a founder effect.

BACKGROUND: Autosomal recessive congenital ichthyoses (ARCI) are inherited disorders produced by mutations in essential genes for the skin function. A low prevalence of this disease has been resported worldwide; however, in a recent study, we identified a large cluster of ARCI families who resided in the High Mountains Region from the Veracruz State, Mexico. Thus, we aimed to identify the causative mutation of ARCI and describe the high prevalence of this disease in this region. METHODS: We selected seven familiar trios and performed whole-exome sequencing to identify the mutation associated with ARCI. To validate the identified mutation, we performed Sanger sequencing in 62 patients, 30 unaffected relatives, and 100 healthy volunteers. Finally, we performed molecular modeling to investigate the possible functional consequences produced by the mutation. RESULTS: We identified a novel homozygous mutation (c.1054C>G [p.Pro352Ala]) in the exon 7 of the TGM1 gene in all the patients. We calculated a prevalence rate of ARCI of 74:100,000 (1:1,348) in the studied communities. Molecular modeling revealed that the mutation leads to a nonconservative amino acid substitution, which is very probably damaging to the protein structure/function. CONCLUSIONS: We report a novel mutation in the TGM1 gene in 62 Mexican patients. The unusually high frequency of this mutation suggests a founder effect; however, further haplotype analysis is necessary to corroborate this hypothesis. In this respect, to our knowledge, the prevalence of ARCI found in the studied communities is the highest observed worldwide.

Synthesis and molecular docking studies of imines as α-glucosidase and α-amylase inhibitors.

Imine functionality is found in many compounds with important biological activity. Thus, the development of novel synthetic approaches for imines is important. In this work, it is propose an easy, eco-friendly and straightforward synthesis pathway of aryl imines under microwave irradiation catalyzed by Alumina-sulfuric acid. In addition, the in vitro enzymatic inhibition, antioxidant activity and molecular docking studies were performed. The aryl imines were isolated with yields in the range of 37-94%. All aryl imines synthesized were evaluated for in vitro inhibitory potential against α-glucosidase and α-amylase enzymes and the results exhibited that the most of the compounds displayed inhibitory activity against both enzymes. The (E)-1-(4-nitrophenyl)-N-(pyridin-2-yl)methanimine (3d) was 1.15-fold more active than acarbose against α-amylase whilst the (E)-1-phenyl-N-(pyridin-2-yl)methanimine (3c) displayed similar activity that acarbose against α-glucosidase. The molecular docking studies in α-glucosidase and α-amylase reveal that aryl imines mainly establish an H-bond with the R2-subtituent and hydrophobic interactions with the R1-subtituent. The docking analysis reveals these synthetic aryl imines 3d-i interact in same active site than acarbose drug in both enzymes.

In silico prediction of structural changes in human papillomavirus type 16 (HPV16) E6 oncoprotein and its variants.

BACKGROUND: HPV16 infection is one of the main risk factors involved in the development of cervical cancer, mainly due to the high oncogenic potential of the viral proteins E6 and E7, which are involved in the different processes of malignant transformation. There is a broad spectrum of intratypical variation of E6, which is reflected in its high diversity, biological behavior, global distribution and risk of causing cervical cancer. Experimental studies have shown that the intratypical variants of the protein E6 from the European variants (E-G350, E-A176/G350, E-C188/G350) and Asian-American variants (AAa and AAc), are capable of inducing the differential expression of genes involved in the development of cervical cancer. RESULTS: An in silico analysis was performed to characterize the molecular effects of these variations using the structure of the HPV16 E6 oncoprotein (PDB: 4XR8; chain H) as a template. In particular, we evaluated the 3D structures of the intratypical variants by structural alignment, ERRAT, Ramachandran plots and prediction of protein disorder, which was further validated by molecular dynamics simulations. Our results, in general, showed no significant changes in the protein 3D structure. However, we observed subtle changes in protein physicochemical features and structural disorder in the N- and C-termini. CONCLUSIONS: Our results showed that mutations in the viral oncogene E6 of six high-risk HPV16 variants are effectively neutral and do not cause significant structural changes except slight variations of structural disorder. As structural disorder is involved in rewiring protein-protein interactions, these results suggest a differential pattern of interaction of E6 with the target protein P53 and possibly different patterns of tumor aggressiveness associated with certain types of variants of the E6 oncoprotein.

EhRabB mobilises the EhCPADH complex through the actin cytoskeleton during phagocytosis of Entamoeba histolytica.

Movement and phagocytosis are clue events in colonisation and invasion of tissues by Entamoeba histolytica, the protozoan causative of human amoebiasis. During phagocytosis, EhRab proteins interact with other functional molecules, conducting them to the precise cellular site. The gene encoding EhrabB is located in the complementary chain of the DNA fragment containing Ehcp112 and Ehadh genes, which encode for the proteins of the EhCPADH complex, involved in phagocytosis. This particular genetic organisation suggests that the three corresponding proteins may be functionally related. Here, we studied the relationship of EhRabB with EhCPADH and actin during phagocytosis. First, we obtained the EhRabB 3D structure to carry out docking analysis to predict the interaction sites involved in the EhRabB protein and the EhCPADH complex contact. By confocal microscopy, transmission electron microscopy, and immunoprecipitation assays, we revealed the interaction among these proteins when they move through different vesicles formed during phagocytosis. The role of the actin cytoskeleton in this event was also confirmed using Latrunculin A to interfere with actin polymerisation. This affected the movement of EhRabB and EhCPADH, as well as the rate of phagocytosis. Mutant trophozoites, silenced in EhrabB gene, evidenced the interaction of this molecule with EhCPADH and strengthened the role of actin during erythrophagocytosis.

Studying the collective motions of the adenosine A2A receptor as a result of ligand binding using principal component analysis.

Adenosine receptors (ARs) belong to family A of GPCRs that are involved in many diseases, including cerebral and cardiac ischemic diseases, immune and inflammatory disorders, etc. Thus, they represent important therapeutic targets to treat these conditions. Computational techniques such as molecular dynamics (MD) simulations permit researchers to obtain structural information about these proteins, and principal component analysis (PCA) allows for the identification of collective motions. There are available structures for the active form (3QAK) and the inactive form (3EML) of A2AR which permit us to gain insight about their activation/inactivation mechanism. In this work, we have proposed an inverse strategy using MD simulations where the active form was coupled to the antagonist caffeine and the inactive form was coupled to adenosine agonist. Moreover, we have included four reported thermostabilizing mutations in the inactive form to study A2AR structural differences under different conditions. Some observations stand out from the PCA studies. For instance, the apo structures showed remarkable similarities, and the principal components (PCs) were rearranged in a ligand-dependent manner. Additionally, the active conformation was less stable compared to the inactive one. Some PCs inverted their direction in the presence of a ligand, and comparison of the PCs between 3EML and 3EML_ADN showed that adenosine induced major changes in the structure of A2AR. Rearrangement of PCs precedes and drives conformational changes that occur after ligand binding. Knowledge about these conformational changes provides important insights about the activity of A2AR.

MicroRNAs in Tumor Cell Metabolism: Roles and Therapeutic Opportunities.

Dysregulated metabolism is a common feature of cancer cells and is considered a hallmark of cancer. Altered tumor-metabolism confers an adaptive advantage to cancer cells to fulfill the high energetic requirements for the maintenance of high proliferation rates, similarly, reprogramming metabolism confers the ability to grow at low oxygen concentrations and to use alternative carbon sources. These phenomena result from the dysregulated expression of diverse genes, including those encoding microRNAs (miRNAs) which are involved in several metabolic and tumorigenic pathways through its post-transcriptional-regulatory activity. Further, the identification of key actionable altered miRNA has allowed to propose novel targeted therapies to modulated tumor-metabolism. In this review, we discussed the different roles of miRNAs in cancer cell metabolism and novel miRNA-based strategies designed to target the metabolic machinery in human cancer.

Structural and functional characterization of the divergent Entamoeba Src using Src inhibitor-1.

BACKGROUND: The abundant number of kinases that Entamoeba histolytica possesses allows us to assume that the regulation of cellular functions by phosphorylation-dephosphorylation processes is very important. However, the kinases responsible for the phosphorylation in Entamoeba spp. vary in the structure of their domains and, therefore, could be responsible for the unusual biological characteristics of this parasite. In higher eukaryotes, Src kinases share conserved structural domains and are very important in the regulation of the actin cytoskeleton. In both Entamoeba histolytica and Entamoeba invadens, the major Src kinase homologue of higher eukaryotes lacks SH3 and SH2 domains, but does have KELCH domains; the latter are part of actin cross-linking proteins in higher eukaryotic cells. METHODS: The function of the EhSrc protein kinase of Entamoeba spp. was evaluated using Src inhibitor-1, microscopy assays, Src kinase activity and western blot. In addition, to define the potential inhibitory mechanism of Src-inhibitor-1 for the amoebic EhSrc protein kinase, molecular dynamic simulations using NAnoscale Molecular Dynamics (NAMD2) program and docking studies were performed with MOE software. RESULTS: We demonstrate that Src inhibitor-1 is able to prevent the activity of EhSrc protein kinase, most likely by binding to the catalytic domain, which affects cell morphology via the disruption of actin cytoskeleton remodeling and the formation of phagocytic structures without an effect on cell adhesion. Furthermore, in E. invadens, Src inhibitor-1 inhibited the encystment process by blocking RhoA GTPase activity, a small GTPase protein of Rho family. CONCLUSIONS: Even though the EhSrc molecule of Entamoeba is not a typical Src, because its divergent amino acid sequence, it is a critical factor in the biology of this parasite via the regulation of actin cytoskeleton remodeling via RhoA GTPase activation. Based on this, we conclude that EhSrc could become a target molecule for the future design of drugs that can prevent the transmission of the disease.

Heterodimerization of the Entamoeba histolytica EhCPADH virulence complex through molecular dynamics and protein-protein docking.

EhCPADH is a protein complex involved in the virulence of Entamoeba histolytica, the protozoan responsible for human amebiasis. It is formed by the EhCP112 cysteine protease and the EhADH adhesin. To explore the molecular basis of the complex formation, three-dimensional models were built for both proteins and molecular dynamics simulations (MDS) and docking calculations were performed. Results predicted that the pEhCP112 proenzyme and the mEhCP112 mature enzyme were globular and peripheral membrane proteins. Interestingly, in pEhCP112, the propeptide appeared hiding the catalytic site (C167, H329, N348); while in mEhCP112, this site was exposed and its residues were found structurally closer than in pEhCP112. EhADH emerged as an extended peripheral membrane protein with high fluctuation in Bro1 and V shape domains. 500 ns-long MDS and protein-protein docking predictions evidenced different heterodimeric complexes with the lowest free energy. pEhCP112 interacted with EhADH by the propeptide and C-terminal regions and mEhCP112 by the C-terminal through hydrogen bonds. In contrast, EhADH bound to mEhCP112 by 442-479 residues, adjacent to the target cell-adherence region (480-600 residues), and by the Bro1 domain (9-349 residues). Calculations of the effective binding free energy and per residue free energy decomposition showed that EhADH binds to mEhCP112 with a higher binding energy than to pEhCP112, mainly through van der Waals interactions and the nonpolar part of solvation energy. The EhADH and EhCP112 structural relationship was validated in trophozoites by immunofluorescence, TEM, and immunoprecipitation assays. Experimental findings fair agreed with in silico results.