Machine learning-based modulation of Ca2+-binding affinity in EF-hand proteins and comparative structural insights into site-specific cooperative binding.

Ca2+-binding proteins are present in almost all living organisms and different types display different levels of binding affinities for the cation. Here, we report two new scoring schemes enabling the user to estimate and manipulate the calcium binding affinities in EF hand containing proteins. To validate this, we designed a unique EF-hand loop capable of binding calcium with high affinity by altering five residues. The N-terminal domain of Entamoeba histolytica calcium-binding protein1 (NtEhCaBP1) is used for site-directed mutagenesis to incorporate the designed loop sequence into the second EF-hand motif of this protein, referred as Nt-EhCaBP1-EF2 mutant. The binding isotherms calculated using ITC calorimetry showed that Nt-EhCaBP1-EF2 mutant site binds Ca2+ with higher affinity than Wt-Nt-EhCaBP1, by ∼600 times. The crystal structure of the mutant displayed more compact Ca2+-coordination spheres in both of its EF loops than the structure of the wildtype protein. The compact coordination sphere of EF-2 causes the bend in the helix-3, which leads to the formation of unexpected hexamer of NtEhCaBP1-EF2 mutant structure. Further dynamic correlation analysis revealed that the mutation in the second EF loop changed the entire residue network of the monomer, resulting in stronger coordination of Ca2+ even in another EF-hand loop.

Analysis of gene expression profiles to study malaria vaccine dose efficacy and immune response modulation.

Malaria is a life-threatening disease, and Africa is still one of the most affected endemic regions despite years of policy to limit infection and transmission rates. Further, studies into the variable efficacy of the vaccine are needed to provide a better understanding of protective immunity. Thus, the current study is designed to delineate the effect of each dose of vaccine on the transcriptional profiles of subjects to determine its efficacy and understand the molecular mechanisms underlying the protection this vaccine provides. Here, we used gene expression profiles of pre and post-vaccination patients after various doses of RTS,S based on samples collected from the Gene Expression Omnibus datasets. Subsequently, differential gene expression analysis using edgeR revealed the significantly (false discovery rate < 0.005) 158 downregulated and 61 upregulated genes between control vs. controlled human malaria infection samples. Further, enrichment analysis of significant genes delineated the involvement of CCL8, CXCL10, CXCL11, XCR1, CSF3, IFNB1, IFNE, IL12B, IL22, IL6, IL27, etc., genes which found to be upregulated after earlier doses but downregulated after the 3rd dose in cytokine-chemokine pathways. Notably, we identified 13 cytokine genes whose expression significantly varied during three doses. Eventually, these findings give insight into the dual role of cytokine responses in malaria pathogenesis. The variations in their expression patterns after various doses of vaccination are linked to the protection as it decreases the severe inflammatory effects in malaria patients. This study will be helpful in designing a better vaccine against malaria and understanding the functions of cytokine response as well.

Global Analysis of Plasmodium falciparum Dihydropteroate Synthase Variants Associated with Sulfadoxine Resistance Reveals Variant Distribution and Mechanisms of Resistance: A Computational-Based Study.

The continual rise in sulfadoxine (SDX) resistance affects the therapeutic efficacy of sulfadoxine-pyrimethamine; therefore, careful monitoring will help guide its prolonged usage. Mutations in Plasmodium falciparum dihydropteroate synthase (Pfdhps) are being surveilled, based on their link with SDX resistance. However, there is a lack of continuous analyses and data on the potential effect of molecular markers on the Pfdhps structure and function. This study explored single-nucleotide polymorphisms (SNPs) in Pfdhps that were isolated in Africa and other countries, highlighting the regional distribution and its link with structure. In total, 6336 genomic sequences from 13 countries were subjected to SNPs, haplotypes, and structure-based analyses. The SNP analysis revealed that the key SDX resistance marker, A437G, was nearing fixation in all countries, peaking in Malawi. The mutation A613S was rare except in isolates from the Democratic Republic of Congo and Malawi. Molecular docking revealed a general loss of interactions when comparing mutant proteins to the wild-type protein. During MD simulations, SDX was released from the active site in mutants A581G and A613S before the end of run-time, whereas an unstable binding of SDX to mutant A613S and haplotype A437A/A581G/A613S was observed. Conformational changes in mutant A581G and the haplotypes A581G/A613S, A437G/A581G, and A437G/A581G/A613S were seen. The radius of gyration revealed an unfolding behavior for the A613S, K540E/A581G, and A437G/A581G systems. Overall, tracking such mutations by the continuous analysis of Pfdhps SNPs is encouraged. SNPs on the Pfdhps structure may cause protein-drug function loss, which could affect the applicability of SDX in preventing malaria in pregnant women and children.

Functional implications of residues of the B' motif of geminivirus replication initiator protein in its helicase activity.

Geminivirus replication initiator protein (Rep) is a multifunctional viral protein required for replication. During the process of viral replication, Rep acts as a site- and strand-specific endonuclease, ligase, ATPase, and helicase. B' motif and ß-hairpin loop of the geminivirus Rep are conserved and important for Rep-mediated helicase activity required for viral replication. To dissect the roles of various amino acid residues of the B' motif and ß-hairpin loop of the geminivirus Rep helicase in its process of unwinding DNA, we investigated eight conserved residues near the ATP active site or the ssDNA contact channel. Our strategy was to mutate these residues to alanines and investigate the effects of these mutations on various biochemical activities associated with DNA unwinding. We looked into the ATP binding, ATP hydrolysis, DNA binding, and DNA unwinding activities of the wild-type and mutant Rep proteins. These investigations showed four residues (Arg279, Asp280, Tyr287, and Pro290) affecting the DNA unwinding activity. A structural model analysis confirmed the B' loop and ssDNA binding loop to be connected through a ß-hairpin structure, suggesting that changes on one loop might affect the other and that these residues function by acting in concert. Viral genomes containing Rep proteins having these mutations in the B' motif did not replicate in planta. Taken together, these results indicated all four residues to be implicated in helicase activity mediated by Rep and demonstrated the significance, for viral replication, of the B' motif and ß-hairpin loop of the C-terminal region of the Rep protein.

Key drug-targeting genes in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal type of cancer. In this study, we undertook a pairwise comparison of gene expression pattern between tumor tissue and its matching adjacent normal tissue for 45 PDAC patients and identified 22 upregulated and 32 downregulated genes. PPI network revealed that fibronectin 1 and serpin peptidase inhibitor B5 were the most interconnected upregulated-nodes. Virtual screening identified bleomycin exhibited reasonably strong binding to both proteins. Effect of bleomycin on cell viability was examined against two PDAC cell lines, AsPC-1 and MIA PaCa-2. AsPC-1 did not respond to bleomycin, however, MIA PaCa-2 responded to bleomycin with an IC50 of 2.6 µM. This implicates that bleomycin could be repurposed for the treatment of PDAC, especially in combination with other chemotherapy agents. In vivo mouse xenograft studies and patient clinical trials are warranted to understand the functional mechanism of bleomycin towards PDAC and optimize its therapeutic efficacy. Furthermore, we will evaluate the antitumor activity of the other identified drugs in our future studies.

An evolutionary non-conserved motif in Helicobacter pylori arginase mediates positioning of the loop containing the catalytic residue for catalysis.

The binuclear metalloenzyme Helicobacter pylori arginase is important for pathogenesis of the bacterium in the human stomach. Despite conservation of the catalytic residues, this single Trp enzyme has an insertion sequence (-153ESEEKAWQKLCSL165-) that is extremely crucial to function. This sequence contains the critical residues, which are conserved in the homolog of other Helicobacter gastric pathogens. However, the underlying basis for the role of this motif in catalytic function is not completely understood. Here, we used biochemical, biophysical and molecular dynamics simulations studies to determine that Glu155 of this stretch interacts with both Lys57 and Ser152. These interactions are essential for positioning of the motif through Trp159, which is located near Glu155 (His122-Trp159-Tyr125 contact is essential to tertiary structural integrity). The individual or double mutation of Lys57 and Ser152 to Ala considerably reduces catalytic activity with Lys57 to Ala being more significant, indicating they are crucial to function. Our data suggest that the Lys57-Glu155-Ser152 interaction influences the positioning of the loop containing the catalytic His133 so that this His can participate in catalysis, thereby providing a mechanistic understanding into the role of this motif in catalytic function. Lys57 was also found only in the arginases of other Helicobacter gastric pathogens. Based on the non-conserved motif, we found a new molecule, which specifically inhibits this enzyme. Thus, the present study not only provides a molecular basis into the role of this motif in function, but also offers an opportunity for the design of inhibitors with greater efficacy.

Accelerated identification of serine racemase inhibitor from Centella asiatica.

Serine racemase (SR) converts the free form of L-serine into D-serine (DS) in the mammalian brain. The DS functions as a co-agonist of N-methyl D-aspartate (NMDA) receptor. The over- activation of NMDA receptor leads to many neurological disorders like stroke, amyotrophic lateral sclerosis, Alzheimer's disease and an effective inhibitor of SR could be a corrective method for the receptor over-activation. We report for the first time here a rapid way of purifying and identifying an inhibitor from medicinal plants known to have the neuro-protective effect. We have purified SR inhibitor from the methanolic extract of Centella asiatica by affinity method. High resolution mass spectrometry and infrared spectroscopy were used to identify the ligand to be madecassoside. We have shown the madecassoside binding in silico and its inhibition of recombinant human serine racemase in vitro and ex vivo.

Metal-induced change in catalytic loop positioning in Helicobacter pylori arginase alters catalytic function.

Arginase is a bimetallic enzyme that utilizes mainly Mn2+ or Co2+ for catalytic function. In human homolog, the substitution of Mn2+ with Co2+ significantly reduces the Km value without affecting the kcat. However, in the Helicobacter pylori counterpart (important for pathogenesis), the kcat increases nearly 4-fold with Co2+ ions both in the recombinant holoenzyme and arginase isolated from H. pylori grown with Co2+ or Mn2+. This suggests that the active site of arginase in the two homologs is modulated differently by these two metal ions. To investigate the underlying mechanism for metal-induced difference in catalytic activity in the H. pylori enzyme, we used biochemical, biophysical and microsecond molecular dynamics simulations studies. The study shows that the difference in binding affinity of Co2+ and Mn2+ ions with the protein is linked to a different positioning of a loop (-122HTAYDSDSKHIHG134-) that contains a conserved catalytic His133. Consequently, the proximity of His133 and conserved Glu281 is varied. We found that the Glu281-His133 interaction is crucial for catalytic function and was previously unexplored in other homologs. We suggest that the proximity difference between these two residues in the Co2+- and Mn2+-proteins alters the proportion of protonated His133 via variation in its pKa. This affects the efficiency of proton transfer - an essential step of l-arginine hydrolysis reaction catalyzed by arginase and thus activity. Unlike in human arginase, the flexibility of the above segment observed in H. pylori homolog suggests that this region in the H. pylori enzyme may be explored to design its specific inhibitors.

Presence of actin binding motif in VgrG-1 toxin of Vibrio cholerae reveals the molecular mechanism of actin cross-linking.

Type VI secretion systems (T6SS) plays a crucial role in Vibrio cholerae mediated pathogenicity. Tip of T6SS is homologous to gp27/gp5 complex or tail spike of T4 bacteriophage. VgrG-1 of V. cholerae T6SS is unusual among other VgrG because its effector domain is trans-located into the cytosol of eukaryotic cells with an additional actin cross-linking domain (ACD) at its C terminal end. ACD of VgrG-1 (VgrG-1-ACD) causes T6SS dependent host cell cytotoxicity through actin cytoskeleton disruption to prevent bacterial engulfment by macrophages. ACD mediated actin cross-linking promotes survival of the bacteria in the small intestine of humans, along with other virulence factors; establishes successful infection with the onset of diarrhoea in humans. Our studies demonstrated VgrG-1-ACD can bind to actin besides actin cross-linking activity. Computational analysis of ACD revealed the presence of actin binding motif (ABM). Mutations in ABM lead to loss of actin binding in vitro. VgrG-1-ACD having the mutated ABM cannot cross-link actin efficiently in vitro and manifests less actin cytoskeleton disruption when transfected in HeLa cells.

Investigations on Binding Pattern of Kinase Inhibitors with PPARγ: Molecular Docking, Molecular Dynamic Simulations, and Free Energy Calculation Studies.

Peroxisome proliferator-activated receptor gamma (PPARγ) is a potential target for the treatment of several disorders. In view of several FDA approved kinase inhibitors, in the current study, we have investigated the interaction of selected kinase inhibitors with PPARγ using computational modeling, docking, and molecular dynamics simulations (MDS). The docked conformations and MDS studies suggest that the selected KIs interact with PPARγ in the ligand binding domain (LBD) with high positive predictive values. Hence, we have for the first time shown the plausible binding of KIs in the PPARγ ligand binding site. The results obtained from these in silico investigations warrant further evaluation of kinase inhibitors as PPARγ ligands in vitro and in vivo.

N-acetyl ornithine deacetylase is a moonlighting protein and is involved in the adaptation of Entamoeba histolytica to nitrosative stress.

Adaptation of the Entamoeba histolytica parasite to toxic levels of nitric oxide (NO) that are produced by phagocytes may be essential for the establishment of chronic amebiasis and the parasite's survival in its host. In order to obtain insight into the mechanism of E. histolytica's adaptation to NO, E. histolytica trophozoites were progressively adapted to increasing concentrations of the NO donor drug, S-nitrosoglutathione (GSNO) up to a concentration of 110 µM. The transcriptome of NO adapted trophozoites (NAT) was investigated by RNA sequencing (RNA-seq). N-acetyl ornithine deacetylase (NAOD) was among the 208 genes that were upregulated in NAT. NAOD catalyzes the deacetylation of N-acetyl-L-ornithine to yield ornithine and acetate. Here, we report that NAOD contributes to the better adaptation of the parasite to nitrosative stress (NS) and that this function does not depend on NAOD catalytic activity. We also demonstrated that glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is detrimental to E. histolytica exposed to NS and that this detrimental effect is neutralized by NAOD or by a catalytically inactive NAOD (mNAOD). These results establish NAOD as a moonlighting protein, and highlight the unexpected role of this metabolic enzyme in the adaptation of the parasite to NS.

Crystal structure of Arabidopsis thaliana calmodulin7 and insight into its mode of DNA binding.

Calmodulin (CaM) is a Ca(2+) sensor that participates in several cellular signaling cascades by interacting with various targets, including DNA. It has been shown that Arabidopsis thaliana CaM7 (AtCaM7) interacts with Z-box DNA and functions as a transcription factor [Kushwaha R et al. (2008) Plant Cell 20, 1747-1759; Abbas N et al. (2014) Plant Cell 26, 1036-1052]. The crystal structure of AtCaM7, and a model of the AtCAM7-Z-box complex suggest that Arg-127 determines the DNA-binding ability by forming crucial interactions with the guanine base. We validated the model using biolayer interferometry, which confirmed that AtCaM7 interacts with Z-box DNA with high affinity. In contrast, the AtCaM2/3/5 isoform does not show any binding, although it differs from AtCaM7 by only a single residue.

Structural insight into ß-Clamp and its interaction with DNA Ligase in Helicobacter pylori.

Helicobacter pylori, a gram-negative and microaerophilic bacterium, is the major cause of chronic gastritis, gastric ulcers and gastric cancer. Owing to its central role, DNA replication machinery has emerged as a prime target for the development of antimicrobial drugs. Here, we report 2Å structure of ß-clamp from H. pylori (Hpß-clamp), which is one of the critical components of DNA polymerase III. Despite of similarity in the overall fold of eubacterial ß-clamp structures, some distinct features in DNA interacting loops exists that have not been reported previously. The in silico prediction identified the potential binders of ß-clamp such as alpha subunit of DNA pol III and DNA ligase with identification of ß-clamp binding regions in them and validated by SPR studies. Hpß-clamp interacts with DNA ligase in micromolar binding affinity. Moreover, we have successfully determined the co-crystal structure of ß-clamp with peptide from DNA ligase (not reported earlier in prokaryotes) revealing the region from ligase that interacts with ß-clamp.

Structural investigation and inhibitory response of halide on phosphoserine aminotransferase from Trichomonas vaginalis.

BACKGROUND: Phosphoserine aminotransferase (PSAT) catalyses the second reversible step of the phosphoserine biosynthetic pathway in Trichomonas vaginalis, which is crucial for the synthesis of serine and cysteine. METHODS: PSAT from T. vaginalis (TvPSAT) was analysed using X-ray crystallography, enzyme kinetics, and molecular dynamics simulations. RESULTS: The crystal structure of TvPSAT was determined to 2.15Å resolution, and is the first protozoan PSAT structure to be reported. The active site of TvPSAT structure was found to be in a closed conformation, and at the active site PLP formed an internal aldimine linkage to Lys 202. In TvPSAT, Val 340 near the active site while it is Arg in most other members of the PSAT family, might be responsible in closing the active site. Kinetic studies yielded Km values of 54 µM and 202 µM for TvPSAT with OPLS and AKG, respectively. Only iodine inhibited the TvPSAT activity while smaller halides could not inhibit. CONCLUSION: Results from the structure, comparative molecular dynamics simulations, and the inhibition studies suggest that iodine is the only halide that can bind TvPSAT strongly and may thus inhibit the activity of TvPSAT. The long loop between ß8 and α8 at the opening of the TvPSAT active site cleft compared to other PSATs, suggests that this loop may help control the access of substrates to the TvPSAT active site and thus influences the enzyme kinetics. GENERAL SIGNIFICANCE: Our structural and functional studies have improved our understanding of how PSAT helps this organism persists in the environment.

Structure-Based Design of Inhibitors of the Crucial Cysteine Biosynthetic Pathway Enzyme O-Acetyl Serine Sulfhydrylase.

The cysteine biosynthetic pathway is of fundamental importance for the growth, survival, and pathogenicity of the many pathogens. This pathway is present in many species but is absent in mammals. The ability of pathogens to counteract the oxidative defences of a host is critical for the survival of these pathogens during their long latent phases, especially in anaerobic pathogens such as Entamoeba histolytica, Leishmania donovani, Trichomonas vaginalis, and Salmonella typhimurium. All of these organisms rely on the de novo cysteine biosynthetic pathway to assimilate sulphur and maintain a ready supply of cysteine. The de novo cysteine biosynthetic pathway, on account of its being important for the survival of pathogens and at the same time being absent in mammals, is an important drug target for diseases such as amoebiasis, trichomoniasis & tuberculosis. Cysteine biosynthesis is catalysed by two enzymes: serine acetyl transferase (SAT) followed by O-acetylserine sulfhydrylase (OASS). OASS is well studied, and with the availability of crystal structures of this enzyme in different conformations, it is a suitable template for structure-based inhibitor development. Moreover, OASS is highly conserved, both structurally and sequence-wise, among the above-mentioned organisms. There have been several reports of inhibitor screening and development against this enzyme from different organisms such as Salmonella typhimurium, Mycobacterium tuberculosis and Entamoeba histolytica. All of these inhibitors have been reported to display micromolar to nanomolar binding affinities for the open conformation of the enzyme. In this review, we highlight the structural similarities of this enzyme in different organisms and the attempts for inhibitor development so far. We also propose that the intermediate state of the enzyme may be the ideal target for the design of effective highaffinity inhibitors.

Ligand-induced conformation changes drive ATP hydrolysis and function in SMARCAL1.

Mutations and deletions in SMARCAL1, an SWI2/SNF2 protein, cause Schimke immuno-osseous dysplasia (SIOD). SMARCAL1 preferentially binds to DNA molecules possessing double-stranded to single-stranded transition regions and mediates annealing helicase activity. The protein is critical for alleviating replication stress and maintaining genome integrity. In this study, we have analysed the ATPase activity of three mutations ­ A468P, I548N and S579L ­ present in SIOD patients. These mutations are present in RecA-like domain I of the protein. Analysis using active DNA-dependent ATPase A domain (ADAAD), an N-terminal deleted construct of bovine SMARCAL1, showed that all three mutants were unable to hydrolyse ATP. Conformational studies indicated that the α-helix and ß-sheet content of the mutant proteins was altered compared to the wild-type protein. Molecular simulation studies confirmed that major structural changes had occurred in the mutant proteins. These changes included alteration of a loop region connecting motif Ia and II. As motif Ia has been implicated in DNA binding, ligand binding studies were done using fluorescence spectroscopy. These studies revealed that the Kd for protein-DNA interaction in the presence of ATP was indeed altered in the case of mutant proteins compared to the wild-type. Finally, in vivo studies were done to complement the in vitro and in silico studies. The results from these experiments demonstrate that mutations in human SMARCAL1 that result in loss in ATPase activity lead to increased replication stress and therefore possibly manifestation of SIOD.

Crystal structure of calcium binding protein-5 from Entamoeba histolytica and its involvement in initiation of phagocytosis of human erythrocytes.

Entamoeba histolytica is the etiological agent of human amoebic colitis and liver abscess, and causes a high level of morbidity and mortality worldwide, particularly in developing countries. There are a number of studies that have shown a crucial role for Ca2+ and its binding protein in amoebic biology. EhCaBP5 is one of the EF hand calcium-binding proteins of E. histolytica. We have determined the crystal structure of EhCaBP5 at 1.9 Šresolution in the Ca2+-bound state, which shows an unconventional mode of Ca2+ binding involving coordination to a closed yet canonical EF-hand motif. Structurally, EhCaBP5 is more similar to the essential light chain of myosin than to Calmodulin despite its somewhat greater sequence identity with Calmodulin. This structure-based analysis suggests that EhCaBP5 could be a light chain of myosin. Surface plasmon resonance studies confirmed this hypothesis, and in particular showed that EhCaBP5 interacts with the IQ motif of myosin 1B in calcium independent manner. It also appears from modelling of the EhCaBP5-IQ motif complex that EhCaBP5 undergoes a structural change in order to bind the IQ motif of myosin. This specific interaction was further confirmed by the observation that EhCaBP5 and myosin 1B are colocalized in E. histolytica during phagocytic cup formation. Immunoprecipitation of EhCaBP5 from total E. histolytica cellular extract also pulls out myosin 1B and this interaction was confirmed to be Ca2+ independent. Confocal imaging of E. histolytica showed that EhCaBP5 and myosin 1B are part of phagosomes. Overexpression of EhCaBP5 increases slight rate (∼20%) of phagosome formation, while suppression reduces the rate drastically (∼55%). Taken together, these experiments indicate that EhCaBP5 is likely to be the light chain of myosin 1B. Interestingly, EhCaBP5 is not present in the phagosome after its formation suggesting EhCaBP5 may be playing a regulatory role.

EhCoactosin stabilizes actin filaments in the protist parasite Entamoeba histolytica.

Entamoeba histolytica is a protist parasite that is the causative agent of amoebiasis, and is a highly motile organism. The motility is essential for its survival and pathogenesis, and a dynamic actin cytoskeleton is required for this process. EhCoactosin, an actin-binding protein of the ADF/cofilin family, participates in actin dynamics, and here we report our studies of this protein using both structural and functional approaches. The X-ray crystal structure of EhCoactosin resembles that of human coactosin-like protein, with major differences in the distribution of surface charges and the orientation of terminal regions. According to in vitro binding assays, full-length EhCoactosin binds both F- and G-actin. Instead of acting to depolymerize or severe F-actin, EhCoactosin directly stabilizes the polymer. When EhCoactosin was visualized in E. histolytica cells using either confocal imaging or total internal reflectance microscopy, it was found to colocalize with F-actin at phagocytic cups. Over-expression of this protein stabilized F-actin and inhibited the phagocytic process. EhCoactosin appears to be an unusual type of coactosin involved in E. histolytica actin dynamics.