Protein-ligand binding site recognition using complementary binding-specific substructure comparison and sequence profile alignment

Prediction of Structure and Molecular Interaction with DNA of BvrR, a Virulence-Associated Regulatory Protein of Brucella

Brucellosis, also known as "undulant fever" is a zoonotic disease caused by Brucella, which is a facultative intracellular bacterium. Despite efforts to eradicate this disease, infection in uncontrolled domestic animals persists in several countries and therefore transmission to humans is common. Brucella evasion of the innate immune system depends on its ability to evade the mechanisms of intracellular death in phagocytic cells. The BvrR-BvrS two-component system allows the bacterium to detect adverse conditions in the environment. The BvrS protein has been associated with genes of virulence factors, metabolism, and membrane transport. In this study, we predicted the DNA sequence recognized by BvrR with Gibbs Recursive Sampling and identified the three-dimensional structure of BvrR using I-TASSER suite, and the interaction mechanism between BvrR and DNA with Protein-DNA docking and molecular dynamics (MD) simulation. Based on the Gibbs recursive Sampling analysis, we found the motif AAHTGC (H represents A, C, and T nucleotides) as a possible sequence recognized by BvrR. The docking and EMD simulation results showed that C-terminal effector domain of BvrR protein is likely to interact with AAHTGC sequence. In conclusion, we predicted the structure, recognition motif, and interaction of BvrR with DNA.

Characterization of TK1646, a highly thermostable 3'-5' single strand specific exonuclease from Thermococcus kodakarensis

Exonucleases catalyze the hydrolysis of terminal phosphodiester bond in nucleic acid. They play important role in maintaining the integrity of DNA in eukaryotes, prokaryotes and archaea. Limited studies have been done on archaeal exonucleases. Here we report molecular cloning of TK1646, a putative exonuclease from the hyperthermophilic archaeon Thermococcus kodakarensis, and expression of the gene in Escherichia coli. Recombinant TK1646, produced in soluble and active form, was purified to apparent homogeneity. Characterization of the recombinant enzyme indicated that it was single strand specific 3'-5' exonuclease which cleaved the substrate DNA after every two nucleotides. It exhibited highest activity at 85-100 °C and pH 9.0. Unique property of TK1646 was its thermostability as it maintained its activity even at 100 °C with a half-life of 180 min. Recombinant TK1646 followed Michaelis-Menten kinetics and exhibited apparent K_m and V_max values of 33 ± 4 μM and 812 ± 48 nmol/min/mg, respectively. To the best of our knowledge this is the most thermostable single strand specific 3'-5' exonuclease characterized to date.

Identification of lead molecules against potential drug target protein MAPK4 from L. donovani: An in-silico approach using docking, molecular dynamics and binding free energy calculation

Leishmaniasis caused by obligate intracellular parasites of genus Leishmania is one of the most neglected tropical diseases threatening 350 million people worldwide. Protein kinases have drawn much attention as potential drug targets due to their important role in various cellular processes. In Leishmania sp. mitogen-activated protein kinase 4 is essential for the parasite survival because of its involvement in various regulatory, apoptotic and developmental pathways. The current study reveals the identification of natural inhibitors of L. donovani mitogen-activated protein kinase-4 (LdMPK4). We have performed in silico docking of 110 natural inhibitors of Leishmania parasite that have been reported earlier and identified two compounds Genistein (GEN) and Chrysin (CHY). The homology model of LdMPK4 was developed, followed by binding affinity studies, and pharmacokinetic properties of the inhibitors were calculated by maintaining ATP as a standard molecule. The modelled structure was deposited in the protein model database with PMDB ID: PM0080988. Molecular dynamic simulation of the enzyme-inhibitor complex along with the free energy calculations over 50 ns showed that GEN and CHY are more stable in their binding. These two molecules, GEN and CHY, can be considered as lead molecules for targeting LdMPK4 enzyme and could emerge as potential LdMPK4 inhibitors.

Computational analysis of high-risk SNPs in human CHK2 gene responsible for hereditary breast cancer: A functional and structural impact

Nowadays CHK2 mutation is studied frequently in hereditary breast and ovarian cancer patients in addition to BRCA1/BRCA2. CHK2 is a tumor suppressor gene that encodes a serine/threonine kinase, also involved in pathways such as DNA repair, cell cycle regulation and apoptosis in response to DNA damage. CHK2 is a well-studied moderate penetrance gene that correlates with third high risk susceptibility gene with an increased risk for breast cancer. Hence before planning large population study, it is better to scrutinize putative functional SNPs of CHK2 using different computational tools. In this study, we have used various computational approaches to identify nsSNPs which are deleterious to the structure and/or function of CHK2 protein that might be causing this disease. Computational analysis was performed by different in silico tools including SIFT, Align GVGD, SNAP-2, PROVEAN, Poly-Phen-2, PANTHER, PhD-SNP, MUpro, iPTREE-STAB, Consurf, InterPro, NCBI Conserved Domain Search tool, ModPred, SPARKS-X, RAMPAGE, Verify-3D, FT Site, COACH and PyMol. Out of 78 nsSNP of human CHK2 gene, seven nsSNPs were predicted functionally most significant SNPs. Among these seven nsSNP, p.Arg160Gly, p.Gly210Arg and p.Ser415Phe are highly conserved residues with conservation score of 9 and three nsSNP were predicted to be involved in post translational modification. The p.Arg160Gly and p.Gly210Arg may interfere in phosphopeptide binding site on FHA conserved domain. The p.Ser415Phe may interfere in formation of activation loop of protein-kinase domain and might interfere in interactions of CHK2 with ligand. The study concludes that mutation of serine to phenylalanine at position 415 is a major mutation in native CHK2 protein which might contribute to its malfunction, ultimately causing disease. This is the first comprehensive study, where CHK2 gene variants are analyzed using in silico tools hence it will be of great help while considering large scale studies and also in developing precision medicines related to these polymorphisms in the era of personalized medicine.

mTM-align: a server for fast protein structure database search and multiple protein structure alignment

With the rapid increase of the number of protein structures in the Protein Data Bank, it becomes urgent to develop algorithms for efficient protein structure comparisons. In this article, we present the mTM-align server, which consists of two closely related modules: one for structure database search and the other for multiple structure alignment. The database search is speeded up based on a heuristic algorithm and a hierarchical organization of the structures in the database. The multiple structure alignment is performed using the recently developed algorithm mTM-align. Benchmark tests demonstrate that our algorithms outperform other peering methods for both modules, in terms of speed and accuracy. One of the unique features for the server is the interplay between database search and multiple structure alignment. The server provides service not only for performing fast database search, but also for making accurate multiple structure alignment with the structures found by the search. For the database search, it takes about 2-5 min for a structure of a medium size (∼300 residues). For the multiple structure alignment, it takes a few seconds for ∼10 structures of medium sizes. The server is freely available at: http://yanglab.nankai.edu.cn/mTM-align/.

PrankWeb: a web server for ligand binding site prediction and visualization

PrankWeb is an online resource providing an interface to P2Rank, a state-of-the-art method for ligand binding site prediction. P2Rank is a template-free machine learning method based on the prediction of local chemical neighborhood ligandability centered on points placed on a solvent-accessible protein surface. Points with a high ligandability score are then clustered to form the resulting ligand binding sites. In addition, PrankWeb provides a web interface enabling users to easily carry out the prediction and visually inspect the predicted binding sites via an integrated sequence-structure view. Moreover, PrankWeb can determine sequence conservation for the input molecule and use this in both the prediction and result visualization steps. Alongside its online visualization options, PrankWeb also offers the possibility of exporting the results as a PyMOL script for offline visualization. The web frontend communicates with the server side via a REST API. In high-throughput scenarios, therefore, users can utilize the server API directly, bypassing the need for a web-based frontend or installation of the P2Rank application. PrankWeb is available at http://prankweb.cz/, while the web application source code and the P2Rank method can be accessed at https://github.com/jendelel/PrankWebApp and https://github.com/rdk/p2rank, respectively.

Porphyromonas gingivalis PgFur Is a Member of a Novel Fur Subfamily With Non-canonical Function

Porphyromonas gingivalis, a keystone pathogen of chronic periodontitis, uses ferric uptake regulator homolog (PgFur) to regulate production of virulence factors. This study aimed to characterize PgFur protein in regard to its structure-function relationship. We experimentally identified the 5′ mRNA sequence encoding the 171-amino-acid-long PgFur protein in the A7436 strain and examined this PgFur version as a full-length protein. PgFur protein did not bind to the canonical Escherichia coli Fur box, but the wild-type phenotype of the mutant Δpgfur strain was restored partially when expression of the ecfur gene was induced from the native pgfur promoter. The full-length PgFur protein contained one zinc atom per protein monomer, but did not bind iron, manganese, or heme. Single cysteine substitutions of CXXC motifs resulted in phenotypes similar to the mutant Δpgfur strain. The modified proteins were produced in E. coli at significantly lower levels, were highly unstable, and did not bind zinc. The pgfur gene was expressed at the highest levels in bacteria cultured for 24 h in the absence of iron and heme or at higher levels in bacteria cultured for 10 h in the presence of protoporphyrin IX source. No influence of high availability of Fe²⁺, Zn²⁺, or Mn²⁺ on pgfur gene expression was observed. Two chromosomal mutant strains producing protein lacking 4 (pgfurΔ4aa) or 13 (pgfurΔ13aa) C-terminal amino acid residues were examined in regard to importance of the C-terminal lysine-rich region. The pgfurΔ13aa strain showed a phenotype typical for the mutant Δpgfur strain, but both the wild-type PgFur protein and its truncated version bound zinc with similar ability. The Δpgfur mutant strain produced higher amounts of HmuY protein compared with the wild-type strain, suggesting compromised regulation of its expression. Potential PgFur ligands, Fe²⁺, Mn²⁺, Zn²⁺, PPIX, or serum components, did not influence HmuY production in the Δpgfur mutant strain. The mutant pgfurΔ4aa and pgfurΔ13aa strains exhibited affected HmuY protein production. PgFur, regardless of the presence of the C-terminal lysine-rich region, bound to the hmu operon promoter. Our data suggest that cooperation of PgFur with partners/cofactors and/or protein/DNA modifications would be required to accomplish its role played in an in vivo multilayer regulatory network.

MPLs-Pred: Predicting Membrane Protein-Ligand Binding Sites Using Hybrid Sequence-Based Features and Ligand-Specific Models

Membrane proteins (MPs) are involved in many essential biomolecule mechanisms as a pivotal factor in enabling the small molecule and signal transport between the two sides of the biological membrane; this is the reason that a large portion of modern medicinal drugs target MPs. Therefore, accurately identifying the membrane protein-ligand binding sites (MPLs) will significantly improve drug discovery. In this paper, we propose a sequence-based MPLs predictor called MPLs-Pred, where evolutionary profiles, topology structure, physicochemical properties, and primary sequence segment descriptors are combined as features applied to a random forest classifier, and an under-sampling scheme is used to enhance the classification capability with imbalanced samples. Additional ligand-specific models were taken into consideration in refining the prediction. The corresponding experimental results based on our method achieved an appreciable performance, with 0.63 MCC (Matthews correlation coefficient) as the overall prediction precision, and those values were 0.604, 0.7, and 0.692, respectively, for the three main types of ligands: drugs, metal ions, and biomacromolecules. MPLs-Pred is freely accessible at http://icdtools.nenu.edu.cn/.

Semi-rational evolution of the 3-(3-hydroxyalkanoyloxy)alkanoate (HAA) synthase RhlA to improve rhamnolipid production in Pseudomonas aeruginosa and Burkholderia glumae

The 3-(3-hydroxyalkanoyloxy)alkanoate (HAA) synthase RhlA is an essential enzyme involved in the biosynthesis of HAAs in Pseudomonas and Burkholderia species. RhlA modulates the aliphatic chain length in rhamnolipids, conferring distinct physicochemical properties to these biosurfactants exhibiting promising industrial and pharmaceutical value. A detailed molecular understanding of substrate specificity and catalytic performance in RhlA could offer protein engineering tools to develop designer variants involved in the synthesis of novel rhamnolipid mixtures for tailored eco-friendly products. However, current directed evolution progress remains limited due to the absence of high-throughput screening methodologies and lack of an experimentally resolved RhlA structure. In the present work, we used comparative modeling and chimeric-based approaches to perform a comprehensive semi-rational mutagenesis of RhlA from Pseudomonas aeruginosa. Our extensive RhlA mutational variants and chimeric hybrids between the Pseudomonas and Burkholderia homologs illustrate selective modulation of rhamnolipid alkyl chain length in both Pseudomonas aeruginosa and Burkholderia glumae. Our results also demonstrate the implication of a putative cap-domain motif that covers the catalytic site of the enzyme and provides substrate specificity to RhlA. This semi-rational mutant-based survey reveals promising 'hot-spots' for the modulation of RL congener patterns and potential control of enzyme activity, in addition to uncovering residue positions that modulate substrate selectivity between the Pseudomonas and Burkholderia functional homologs. DATABASE: Model data are available in the PMDB database under the accession number PM0081867.

Molecular Informatics Studies of the Iron-Dependent Regulator (ideR) Reveal Potential Novel Anti-Mycobacterium ulcerans Natural Product-Derived Compounds

Buruli ulcer is a neglected tropical disease caused by the bacterium Mycobacterium ulcerans. Its virulence is attributed to the dermo-necrotic polyketide toxin mycolactone, whose synthesis is regressed when its iron acquisition system regulated by the iron-dependent regulator (ideR) is deactivated. Interfering with the activation mechanism of ideR to inhibit the toxin’s synthesis could serve as a possible cure for Buruli ulcer. The three-dimensional structure of the ideR for Mycobacterium ulcerans was generated using homology modeling. A library of 832 African natural products (AfroDB), as well as five known anti-mycobacterial compounds were docked against the metal binding site of the ideR. The area under the curve (AUC) values greater than 0.7 were obtained for the computed Receiver Operating Characteristics (ROC) curves, validating the docking protocol. The identified top hits were pharmacologically profiled using Absorption, Distribution, Metabolism, Elimination and Toxicity (ADMET) predictions and their binding mechanisms were characterized. Four compounds with ZINC IDs ZINC000018185774, ZINC000095485921, ZINC000014417338 and ZINC000005357841 emerged as leads with binding energies of −7.7 kcal/mol, −7.6 kcal/mol, −8.0 kcal/mol and −7.4 kcal/mol, respectively. Induced Fit Docking (IFD) was also performed to account for the protein’s flexibility upon ligand binding and to estimate the best plausible conformation of the complexes. Results obtained from the IFD were consistent with that of the molecular docking with the lead compounds forming interactions with known essential residues and some novel critical residues Thr14, Arg33 and Asp17. A hundred nanoseconds molecular dynamic simulations of the unbound ideR and its complexes with the respective lead compounds revealed changes in the ideR’s conformations induced by ZINC000018185774. Comparison of the lead compounds to reported potent inhibitors by docking them against the DNA-binding domain of the protein also showed the lead compounds to have very close binding affinities to those of the potent inhibitors. Interestingly, structurally similar compounds to ZINC000018185774 and ZINC000014417338, as well as analogues of ZINC000095485921, including quercetin are reported to possess anti-mycobacterial activity. Also, ZINC000005357841 was predicted to possess anti-inflammatory and anti-oxidative activities, which are relevant in Buruli ulcer and iron acquisition mechanisms, respectively. The leads are molecular templates which may serve as essential scaffolds for the design of future anti-mycobacterium ulcerans agents.

Comparative analysis of C-type lectin domain proteins in the ghost moth, Thitarodes xiaojinensis (Lepidoptera: Hepialidae)

Insects have a large family of C-type lectins involved in cell adhesion, pathogen recognition and activation of immune responses. In this study, 32 transcripts encoding C-type lectin domain proteins (CTLDPs) were identified from the Thitarodes xiaojinensis transcriptome. According to their domain structures, six CTLDPs with one carbohydrate-recognition domain (CRD) were classified into the CTL-S subfamily. The other 23 CTLDPs with two CRDs were grouped into the immulectin (IML) subfamily. The remaining three with extra regulatory domains were sorted into the CTL-X subfamily. Phylogenetic analysis showed that CTL-S and CTL-X members from different insects could form orthologous groups. In contrast, no T. xiaojinensis IML orthologues were found in other insects. Remarkable lineage-specific expansion in this subfamily was observed reflecting that these CTLDPs, as important receptors, have evolved diversified members in response to a variety of microbes. Prediction of binding ligands revealed that T. xiaojinensis, a cold-adapted species, conserved the ability of CRDs to combine with Ca²⁺ to keep its receptors from freezing. Comparative analysis of induction of CTLDP genes after different immune challenges indicated that IMLs might play critical roles in immune defenses. This study examined T. xiaojinensis CTLDPs and provides a basis for further studies of their characteristics.

Structure of the SARS-Unique Domain C From the Bat Coronavirus HKU4

Coronaviruses (CoVs) that cause infections such as severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome phylogenetically originate from bat CoVs. The coronaviral nonstructural protein 3 (nsp3) has been implicated in viral replication, polyprotein cleavage, and host immune interference. We report the structure of the C domain from the SARS-Unique Domain of bat CoV HKU4. The protein has a frataxin fold, consisting of 5 antiparallel β strands packed against 2 α helices. Bioinformatics analyses and nuclear magnetic resonance experiments were conducted to investigate the function of HKU4 C. The results showed that HKU4 C engages in protein-protein interactions with the nearby M domain of nsp3. The HKU4 C residues involved in protein-protein interactions are conserved in group 2c CoVs, indicating a conserved function.

Molecular cloning, expression and characterization of a novel feruloyl esterase from a soil metagenomic library with phthalate-degrading activity

OBJECTIVES

To discover novel feruloyl esterases (FAEs) by the function-driven screening procedure from soil metagenome.

RESULTS

A novel FAE gene bds4 was isolated from a soil metagenomic library and over-expressed in Escherichia coli. The recombinant enzyme BDS4 was purified to homogeneity with a predicted molecular weight of 38.8 kDa. BDS4 exhibited strong activity (57.05 U/mg) toward methyl ferulate under the optimum pH and temperature of 8.0 and 37°C. Based on its amino acid sequence and model substrates specificity, BDS4 was classified as a type-C FAE. The quantity of the releasing ferulic acid can be enhanced significantly in the presence of xylanase compared with BDS4 alone from de-starched wheat bran. In addition, BDS4 can also hydrolyze several phthalates such as diethyl phthalate, dimethyl phthalate and dibutyl phthalate.

CONCLUSION

The current investigation discovered a novel FAE with phthalate-degrading activity and highlighted the usefulness of metagenomic approaches as a powerful tool for discovery of novel FAEs.

In silico and in vitro studies of lupeol and iso-orientin as potential antidiabetic agents in a rat model

Background: In silico characterization can help to explain the interaction between molecules and predict three-dimensional structures. Various studies have confirmed the glucose-lowering effects of plant extracts, ie, lupeol and iso-orientin, which enable them to be used as antidiabetic agents. Purpose: Aims of the present study were to evaluate the hypoglycemic activities of lupeol and iso-orientin in a rat model. The study proposed the effects of alloxan on blood glucose level, body weight, and oxidative stress. Materials and Methods: Thirty (n=30) Wistar albino rats were divided into six groups and were subjected to different combinations of the compounds. Levels of different stress markers, ie, malondialdehyde, superoxide dismutase, catalase, nitric oxide, glutathione, glutathione peroxide, glutathione reductase, and blood glucose levels were estimated with their respective methods. Whereas, for their in silico analysis, identified target proteins, GPR40, glucose-6-phosphatase, UCP2, glycogen phosphorylase, aldose reductase, and glucose transporter-4 were docked with lupeol and iso-orientin. Three-dimensional structures were predicted by ERRAT, Rampage, Verify3D, threading and homology approaches. Results: Blood glucose levels were significantly increased in rats receiving intraperitoneal injection of alloxan (208±6.94 mg/dL) as compared to controls (90±7.38 mg/dL). Infected rats were administered plant extracts; combined treatment of both extracts (lupeol+iso-orientin) significantly reduced the levels of blood glucose (129.06±6.29 mg/dL) and improved the antioxidant status. Fifteen structures of each selected protein were evaluated using various techniques. Consequently, satisfactory quality factors [GPR40 (96.41%), glucose-6-phosphatase (96.56%), UCP2 (72.56%), glycogen phosphorylase (87.24%), aldose reductase (82.46%), and glucose transporter-4 (94.29%)] were selected. Molecular docking revealed interacting residues, effective drug properties and their binding affinities (ie, -8.9 to -12.6 Kcal/mol). Conclusion: Results of the study affirmed the antidiabetic activities of lupeol and iso-orientin. Administration of these extracts (either individually or in combination) significantly reduced blood glucose levels and oxidative stress. Hence, it may be considered beneficial in the treatment of diabetes.

Three dimensional structure prediction of panomycocin, a novel Exo-β-1,3-glucanase isolated from Wickerhamomyces anomalus NCYC 434 and the computational site-directed mutagenesis studies to enhance its thermal stability for therapeutic applications

Panomycocin is a naturally produced potent antimycotic/antifungal protein secreted by the yeast Wickerhamomyces anomalus NCYC 434 with an exo-β-1,3-glucanase activity. In this study the three dimensional structure of panomycocin was predicted and the computational site-directed mutagenesis was performed to enhance its thermal stability in liquid formulations over the body temperature for topical therapeutic applications. Homology modeling was performed with MODELLER and I-TASSER. Among the generated models, the model with the lowest energy and DOPE score was selected for further loop modeling. The loop model was optimized and the reliability of the model was confirmed with ERRAT, Verify 3D and Ramachandran plot values. Enhancement of the thermal stability of the model was done using contemporary servers and programs such as SPDBViewer, CNA, I-Mutant2.0, Eris, AUTO-MUTE and MUpro. In the region outside the binding site of the model Leu52 Arg, Phe223Arg and Gly254Arg were found to be the best thermostabilizing mutations with 6.26 K, 6.26 K and 8.27 K increases, respectively. In the binding site Glu186Arg was found to be the best thermostabilizer mutation with a 9.58 K temperature increase. The results obtained in this study led us to design a mutant panomycocin that can be used as a novel antimycotic/antifungal drug in a liquid formulation for topical applications over the normal body temperature.

In silico structural elucidation of RNA-dependent RNA polymerase towards the identification of potential Crimean-Congo Hemorrhagic Fever Virus inhibitors

The Crimean-Congo Hemorrhagic Fever virus (CCHFV) is a segmented negative single-stranded RNA virus (-ssRNA) which causes severe hemorrhagic fever in humans with a mortality rate of ~50%. To date, no vaccine has been approved. Treatment is limited to supportive care with few investigational drugs in practice. Previous studies have identified viral RNA dependent RNA Polymerase (RdRp) as a potential drug target due to its significant role in viral replication and transcription. Since no crystal structure is available yet, we report the structural elucidation of CCHFV-RdRp by in-depth homology modeling. Even with low sequence identity, the generated model suggests a similar overall structure as previously reported RdRps. More specifically, the model suggests the presence of structural/functional conserved RdRp motifs for polymerase function, the configuration of uniform spatial arrangement of core RdRp sub-domains, and predicted positively charged entry/exit tunnels, as seen in sNSV polymerases. Extensive pharmacophore modeling based on per-residue energy contribution with investigational drugs allowed the concise mapping of pharmacophoric features and identified potential hits. The combination of pharmacophoric features with interaction energy analysis revealed functionally important residues in the conserved motifs together with in silico predicted common inhibitory binding modes with highly potent reference compounds.

Two purified proteins from royal jelly with in vitro dual anti-hepatic damage potency: Major royal jelly protein 2 and its novel isoform X1

Liver diseases are serious life-threating conditions that should be controlled. Here, we identify a protein fraction from royal-jelly (RJ) that represents the most effective composite against CCl4-induced hepatotoxicity and HepG2 cell growth. Two closely related proteins were purified from this fraction by a new simple method and identified by MALDI-TOF MS as major RJ protein 2 (MRJP2) and its predicted isoform X1. The in silico assessment (3D structures and functions) of these proteins were performed using Iterative Threading ASSEmbly Refinement (I-TASSER) analysis and RAMPAGE program. These two purified proteins were able to relieve the necrotic hepatocytes (by 60.4%) via reducing tumor necrosis factor (TNF)-α, mixed lineage kinase domain-like protein (MLKL) and intracellular reactive species. The latter effects associated with improving hepatocyte functions. Furthermore, they revealed the potent anticancer effect via induction of caspase-dependent apoptosis and controlling the expression of both Bcl-2 and p53 in HepG2 cells. Thus, MRJP2 and its isoform X1 can be a promising dual strategy for fighting hepatic injury and cancer in future animal and human studies.

Novel interactions of Selenium Binding Protein family with the PICOT containing proteins AtGRXS14 and AtGRXS16 in Arabidopsis thaliana

During abiotic stress the primary symptom of phytotoxicity can be ROS production which is strictly regulated by ROS scavenging pathways involving enzymatic and non-enzymatic antioxidants. Furthermore, ROS are well-described secondary messengers of cellular processes, while during the course of evolution, plants have accomplished high degree of control over ROS and used them as signalling molecules. Glutaredoxins (GRXs) are small and ubiquitous glutathione (GSH) -or thioredoxin reductase (TR)-dependent oxidoreductases belonging to the thioredoxin (TRX) superfamily which are conserved in most eukaryotes and prokaryotes. In Arabidopsis thaliana GRXs are subdivided into four classes playing a central role in oxidative stress responses and physiological functions. In this work, we describe a novel interaction of AtGRXS14 with the Selenium Binding Protein 1 (AtSBP1), a protein proposed to be integrated in a regulatory network that senses alterations in cellular redox state and acts towards its restoration. We further show that SBP protein family interacts with AtGRXS16 that also contains a PICOT domain, like AtGRXS14.

Novel mutation in FTHL17 gene in pedigree with 46,XY pure gonadal dysgenesis

OBJECTIVE

To identify the genetic cause of a pedigree with four patients with 46,XY pure gonadal dysgenesis (PGD).

DESIGN

Genetic mutation study.

SETTING

Academic medical center.

PATIENT(S)

Four first cousins, from three households of a Chinese pedigree, affected by 46,XY PGD.

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

The patients were studied from clinical and genetic perspectives. Whole-genome sequencing was conducted in family members.

RESULT(S)

Four first cousins in the third generation were affected by 46,XY PGD. A specific familial characteristic was the prevalence of as high as 100% of gonadal tumors in patients. Whole-genome sequencing identified a new ferritin heavy chain-like 17 (FTHL17) mutation, c.GA442_443TT (p.E148L), which has the potential to interfere with protein function and cause 46,XY PGD. Moreover, the location (Xp21.2) of the FTHL17 gene proves that the family is X-linked recessive. In vitro functional study revealed that the perturbation of FTHL17 caused the decrease of protein expression and cell proliferation.

CONCLUSION(S)

We describe the first 46,XY PGD pedigree that may be attributed to mutations of the FTHL17 gene. We speculated that the FTHL17 gene is involved in the testis-determining pathway and tumorigenesis.

Structural studies on dihydrouridine synthase A (DusA) from Pseudomonas aeruginosa

Dihydrouridination is one of the abundant modifications in tRNA editing. The presence of dihydrouridine is attributed to tRNA stability desired for the efficient gene translation process. The conversion of uridine to dihydrouridine is catalyzed by flavine containing enzyme called dihydrouridine synthase (Dus). We report first ever information about DusA enzyme from Pseudomonas aeruginosa in form of structural and functional studies. The gene coding for DusA from P. aeruginosa (PADusA) was cloned, expressed and purified, using recombinant DNA technology methods. Thermal and chemical stability of PADusA was determined with respect to temperature and urea-induced equilibrium unfolding experiments, with monitoring the change of ellipticity at 200-260 nm by Circular Dichroism (CD) spectroscopy. Unfolding studies revealed that PADusA has acquired a stable tertiary structure fold with a T_m value of 46.2 °C and C_m of 2.7 M for urea. The enzyme contains 43% α-helices and 16% β-strands. The three dimensional structure of PADusA was modeled using insilico methods. In order to understand the mechanism of substrate recognition and catalysis, tRNA and puromycin were docked on PADusA structure and their binding was analyzed. The structural features suggested that PADusA may also form a novel target for structure based drug design of antimicrobial agents.

Detection and characterization of an albumin-like protein in Leishmania donovani

The protozoan parasite, Leishmania donovani, undergoes several molecular adaptations and secretes many effector molecules for host cell manipulation and successful parasitism. The current study identifies an albumin-like secretory protein, expressed in its extracellular promastigote forms. A leishmanial complementary DNA sequence of a partial gene has been cloned, and the encoded peptide (14 kD) is used for the production of polyclonal antibody. This targeted antibody identifies a large native protein (66.421 kD), expressed stage-specifically in promastigotes. Through electron microscopic studies, the native protein is found to be localized in the flagellar pocket and flagella and at the surface of the promastigotes. This native protein is purified with the same customized antibody for future characterization and sequencing. The sequence analysis reveals its homology with the mammalian serum albumin. It is evidenced from in silico studies that this albumin-like protein remains associated with long-chain fatty acids while in vitro studies indicate its close association with membrane cholesterol. Since antibody-mediated blocking compromises the parasite infectivity, these leishmanial albumin-like molecules are hereby proposed to play an instrumental role in the infectivity of L. donovani to peripheral blood monocyte cells. Thus, identification and characterization of an albumin-like protein in L. donovani promastigotes may be interpreted as a molecular adaptation candidate. It may be hypothesized that the parasite mimics the mammalian system for importing fatty acids into the intracellular amastigotes, facilitating its host cell infectivity.

Identification of deleterious SNPs and their effects on BCL11A, the master regulator of fetal hemoglobin expression

The B-cell lymphoma/leukemia 11A protein (encoded by BCL11A gene) is a key regulator of fetal-to-adult hemoglobin switching as seen in post-natal life. Although genetic polymorphisms like SNPs in BCL11A gene are expected to affect fetal hemoglobin (HbF) expression levels, yet their implications are poorly studied. This study utilizes a computational approach to identify the deleterious SNPs which may affect the structure and function of BCL11A protein. The study also generated a 3D structure of native and mutants. The analysis identified two SNPs in BCL11A as highly deleterious: N391K and C414S which are expected to affect structure and stability of the protein. According to conservation analysis, both residues N391 and C414 were identified as highly conserved. Additionally, post-translational modification sites were predicted at both sites. Ligand binding sites were also predicted in N391 and C414. Therefore, N391K and C414S in BCL11A can considered as important candidates to mediate HbF variation.

Comparative genome analysis provides novel insight into the interaction of Aquimarina sp. AD1, BL5 and AD10 with their macroalgal host

The Aquimarina genus is widely distributed throughout the marine environment, however little is understood regarding its ecological role, particularly when in association with eukaryotic hosts. Here, we examine the genomes of two opportunistic pathogens, Aquimarina sp. AD1 and BL5, and a non-pathogenic strain Aquimarina sp. AD10, that were isolated from diseased individuals of the red alga Delisea pulchra. Each strain encodes multiple genes for the degradation of marine carbohydrates and vitamin biosynthesis. These traits are hypothesised to promote nutrient exchange between the Aquimarina strains and their algal host, facilitating a close symbiotic relationship. Moreover, each strain harbours the necessary genes for the assembly of a Type 9 Secretion System (T9SS) and the associated gliding motility apparatus. In addition to these common features, pathogenic strains AD1 and BL5, encode genes for the production of flexirubin type pigments and a number of unique non-ribosomal peptide synthesis (NRPS) gene clusters, suggesting a role for these uncharacterised traits in virulence. This study provides valuable insight into the potential ecological role of Aquimarina in the marine environment and the complex factors driving pathogenesis and symbiosis in this genus.

Predicting protein-ligand binding residues with deep convolutional neural networks

Background

Ligand-binding proteins play key roles in many biological processes. Identification of protein-ligand binding residues is important in understanding the biological functions of proteins. Existing computational methods can be roughly categorized as sequence-based or 3D-structure-based methods. All these methods are based on traditional machine learning. In a series of binding residue prediction tasks, 3D-structure-based methods are widely superior to sequence-based methods. However, due to the great number of proteins with known amino acid sequences, sequence-based methods have considerable room for improvement with the development of deep learning. Therefore, prediction of protein-ligand binding residues with deep learning requires study.

Results

In this study, we propose a new sequence-based approach called DeepCSeqSite for ab initio protein-ligand binding residue prediction. DeepCSeqSite includes a standard edition and an enhanced edition. The classifier of DeepCSeqSite is based on a deep convolutional neural network. Several convolutional layers are stacked on top of each other to extract hierarchical features. The size of the effective context scope is expanded as the number of convolutional layers increases. The long-distance dependencies between residues can be captured by the large effective context scope, and stacking several layers enables the maximum length of dependencies to be precisely controlled. The extracted features are ultimately combined through one-by-one convolution kernels and softmax to predict whether the residues are binding residues. The state-of-the-art ligand-binding method COACH and some of its submethods are selected as baselines. The methods are tested on a set of 151 nonredundant proteins and three extended test sets. Experiments show that the improvement of the Matthews correlation coefficient (MCC) is no less than 0.05. In addition, a training data augmentation method that slightly improves the performance is discussed in this study.

Conclusions

Without using any templates that include 3D-structure data, DeepCSeqSite significantlyoutperforms existing sequence-based and 3D-structure-based methods, including COACH. Augmentation of the training sets slightly improves the performance. The model, code and datasets are available at https://github.com/yfCuiFaith/DeepCSeqSite.

Gene-Specific Intron Retention Serves as Molecular Signature that Distinguishes Melanoma from Non-Melanoma Cancer Cells in Greek Patients

BACKGROUND

Skin cancer represents the most common human malignancy, and it includes BCC, SCC, and melanoma. Since melanoma is one of the most aggressive types of cancer, we have herein attempted to develop a gene-specific intron retention signature that can distinguish BCC and SCC from melanoma biopsy tumors.

METHODS

Intron retention events were examined through RT-sqPCR protocols, using total RNA preparations derived from BCC, SCC, and melanoma Greek biopsy specimens. Intron-hosted miRNA species and their target transcripts were predicted via the miRbase and miRDB bioinformatics platforms, respectively. Ιntronic ORFs were recognized through the ORF Finder application. Generation and visualization of protein interactomes were achieved by the IntAct and Cytoscape softwares, while tertiary protein structures were produced by using the I-TASSER online server.

RESULTS

c-MYC and Sestrin-1 genes proved to undergo intron retention specifically in melanoma. Interaction maps of proteins encoded by genes being potentially targeted by retained intron-accommodated miRNAs were generated and SRPX2 was additionally delivered to our melanoma-specific signature. Novel ORFs were identified in MCT4 and Sestrin-1 introns, with potentially critical roles in melanoma development.

CONCLUSIONS

The property of c-MYC, Sestrin-1, and SRPX2 genes to retain specific introns could be clinically used to molecularly differentiate non-melanoma from melanoma tumors.

Molecular characterization of two novel proteins All1122 and Alr0750 of Anabaena PCC 7120 conferring tolerance to multiple abiotic stresses in Escherichia coli

In- silico and functional genomics approaches have been used to determine cellular functions of two hypothetical proteins All1122 and Alr0750 of Anabaena sp. PCC 7120. Motif analysis and multiple sequence alignment predicted them as typical α/β ATP binding universal stress family protein-A (UspA) with G-(2×)-G-(9×)-G(S/T) as conserved motif. qRT-PCR data under UV-B, NaCl, heat, As, CdCl_2, mannitol and methyl viologen registered approximately 1.4 to 4.3 fold induction of all1122 and alr0750 thus confirming their multiple abiotic stress tolerance potential. The recombinant E. coli (BL21) cells harboring All1122 and Alr0750 showed 12-41% and 23-41% better growth respectively over wild type control under said abiotic stresses thus revalidating their stress coping ability. Functional complementation on heterologous expression in UspA mutant E. coli strain LN29MG1655 (ΔuspA::Kan) attested their UspA family membership. This study tempted us to suggest that recombinant Anabaena PCC 7120 over expressing all1122 and alr0750 might contribute to the nitrogen economy in paddy fields experiencing array of abiotic stresses including drought and nutrient limitation.

A combined computational strategy of sequence and structural analysis predicts the existence of a functional eicosanoid pathway in Drosophila melanogaster

This study reports on a putative eicosanoid biosynthesis pathway in Drosophila melanogaster and challenges the currently held view that mechanistic routes to synthesize eicosanoid or eicosanoid-like biolipids do not exist in insects, since to date, putative fly homologs of most mammalian enzymes have not been identified. Here we use systematic and comprehensive bioinformatics approaches to identify most of the mammalian eicosanoid synthesis enzymes. Sensitive sequence analysis techniques identified candidate Drosophila enzymes that share low global sequence identities with their human counterparts. Twenty Drosophila candidates were selected based upon (a) sequence identity with human enzymes of the cyclooxygenase and lipoxygenase branches, (b) similar domain architecture and structural conservation of the catalytic domain, and (c) presence of potentially equivalent functional residues. Evaluation of full-length structural models for these 20 top-scoring Drosophila candidates revealed a surprising degree of conservation in their overall folds and potential analogs for functional residues in all 20 enzymes. Although we were unable to identify any suitable candidate for lipoxygenase enzymes, we report structural homology models of three fly cyclooxygenases. Our findings predict that the D. melanogaster genome likely codes for one or more pathways for eicosanoid or eicosanoid-like biolipid synthesis. Our study suggests that classical and/or novel eicosanoids mediators must regulate biological functions in insects–predictions that can be tested with the power of Drosophila genetics. Such experimental analysis of eicosanoid biology in a simple model organism will have high relevance to human development and health.

Identification of coenzyme-binding proteins with machine learning algorithms

The coenzyme-binding proteins play a vital role in the cellular metabolism processes, such as fatty acid biosynthesis, enzyme and gene regulation, lipid synthesis, particular vesicular traffic, and β-oxidation donation of acyl-CoA esters. Based on the theory of Star Graph Topological Indices (SGTIs) of protein primary sequences, we proposed a method to develop a first classification model for predicting protein with coenzyme-binding properties. To simulate the properties of coenzyme-binding proteins, we created a dataset containing 2897 proteins, among 456 proteins functioned as coenzyme-binding activity. The SGTIs of peptide sequence were calculated with Sequence to Star Network (S2SNet) application. We used the SGTIs as inputs to several classification techniques with a machine learning software - Weka. A Random Forest classifier based on 3 features of the embedded and non-embedded graphs was identified as the best predictive model for coenzyme-binding proteins. This model developed was with the true positive (TP) rate of 91.7%, false positive (FP) rate of 7.6%, and Area Under the Receiver Operating Characteristic Curve (AUROC) of 0.971. The prediction of new coenzyme-binding activity proteins using this model could be useful for further drug development or enzyme metabolism researches.

New Binding Sites, New Opportunities for GPCR Drug Discovery

Many central biological events rely on protein-ligand interactions. The identification and characterization of protein-binding sites for ligands are crucial for the understanding of functions of both endogenous ligands and synthetic drug molecules. G protein-coupled receptors (GPCRs) typically detect extracellular signal molecules on the cell surface and transfer these chemical signals across the membrane, inducing downstream cellular responses via G proteins or β-arrestin. GPCRs mediate many central physiological processes, making them important targets for modern drug discovery. Here, we focus on the most recent breakthroughs in finding new binding sites and binding modes of GPCRs and their potentials for the development of new medicines.

Phytochemical composition and the cholinesterase and xanthine oxidase inhibitory properties of seed extracts from the Washingtonia filifera palm fruit

The chemical composition and biological properties of palm Washingtonia filifera (Lindl.) H. Wendl. seeds are seldom studied.

Genome-wide identification and biochemical characterization of the UGT88F subfamily in Malus x domestica Borkh

The UDP-glycosyltransferase UGT88F subfamily has been described first in Malus x domestica with the characterization of UGT88F1. Up to now UGT88F1 was one of the most active UGT glycosylating dihydrochalcones in vitro. The involvement of UGT88F1 in phloridzin (phloretin 2'-O-glucoside) synthesis, the main apple tree dihydrochalcone, was further confirmed in planta. Since the characterization of UGT88F1, this new UGT subfamily has been poorly studied probably because it seemed restricted to Maloideae. In the present study, we investigate the apple tree genome to identify and biochemically characterize the whole UGT88F subfamily. The apple tree genome contains five full-length UGT88F genes out of which three newly identified members (UGT88F6, UGT88F7 and UGT88F8) and a pseudogene. These genes are organized into two genomic clusters resulting from the recent global genomic duplication event in the apple tree. We show that recombinant UGT88F8 protein specifically glycosylates phloretin in the 2'OH position to synthetize phloridzin in vitro and was therefore named UDP-glucose: phloretin 2'-O-glycosyltransferase. The Km values of UGT88F8 are 7.72 μM and 10.84 μM for phloretin and UDP-glucose respectively and are in the same range as UGT88F1 catalytic parameters thus constituting two isoforms. Co-expression patterns of both UGT88F1 and UGT88F8 argue for a redundant function in phloridzin biosynthesis in planta. Contrastingly, recombinant UGT88F6 protein is able to glycosylate in vitro a wide range of flavonoids including flavonols, flavones, flavanones, chalcones and dihydrochalcones, although flavonols are the preferred substrates, e.g. Km value for kaempferol is 2.1 μM. Depending on the flavonoid, glycosylation occurs at least on the 3-OH and 7-OH positions. Therefore UGT88F6 corresponds to an UDP-glucose: flavonoid 3/7-O-glycosyltransferase. Finally, a molecular modeling study highlights a very high substitution rate of residues in the acceptor binding pocket between UGT88F8 and UGT88F6 which is responsible for the enzymes divergence in substrate and regiospecificity, despite an overall high protein homology.

Molecular modelling and competitive inhibition of a Mycobacterium tuberculosis multidrug-resistance efflux pump

Tuberculosis is a major cause of mortality and morbidity in developing countries, and the emergency of multidrug and extensive drug resistance cases is an utmost issue on the control of the disease. Despite the efforts on the development of new antibiotics, eventually there will be strains resistant to them as well. Efflux plays an important role in the evolution of resistance in Mycobacterium tuberculosis. Tap is an important efflux pump associated with tuberculosis resistant to isoniazid, rifampicine and ofloxacin and with multidrug resistance. The development of efflux inhibitors for Tap could raise the effectiveness of second line drugs and reduce the duration of the current treatment. Therefore the objective of this study is to build a reliable molecular model of Tap efflux pump and test the possible competitive inhibition between efflux inhibitors and antibiotics in the optimized structure. We built twenty five Tap models with molecular modelling to elect the best according to the results of the validation analysis. The elected model went through to a 100 ns molecular dynamics simulation in a lipid bilayer, and the resulting optimized structure was used in docking studies to test if the used efflux inhibitors may act via competitive inhibition on antibiotics. The validation results pointed the model built by Phyre2 as the closest to a possible native Tap structure, and therefore it was the elected model. RSMD analysis revealed the model is stable, where the predicted binding site stabilized between 15 and 20 ns, maintaining the RMSD at around 0.35 Å throughout the molecular dynamics simulation in a lipid bilayer. Therefore this model is reliable and can also be used for further studies. The docking studies showed a possibility of competitive inhibition by NUNL02 on ofloxacin and bedaquiline, and by verapamil on ofloxacin and rifampicin. This presents the possibility that NUNL02 and verapamil are possible inhibitors of Tap efflux and highlights the importance of including efflux inhibitors as adjuvants to the tuberculosis therapy, as it indicates a possible extrusion of ofloxacin, rifampicin and bedaquilin by Tap.

In silico prediction, phylogenetic and bioinformatic analysis of SoPCS gene, survey of its protein characterization and gene expression in response to cadmium in Saccharum officinarum

Phytochelatin synthase isolated from microorganisms, yeasts, algae and plant, serve a fundamental role in reducing heavy metals. In this research the in silico PCS gene structure (SoPCS) of sugarcane, its secondary and 3D protein structure, physicochemical properties, cell localization and phylogenetic tree were predicted utilizing bioinformatics tools. SoPCS expression in the leaves and roots of sugarcane in tissue culture treated with cadmium was also studied utilizing real time PCR. The predicted SoPCS gene contains 1524 nucleotides, a protein encoded with 508 amino acids of which the molecular weight is 55953.3 Da, 6 exons and 5 introns. The subcellular position of the enzyme is mitochondrion or cytoplasmic. Two domains belonging to the phytochelatin synthase family with similar features was found in Pfam having more than 97% similarity with the predicted SoPCS protein. Phylogeny analyses of plant species were well isolated from other organisms. Ten disulfide-bonded cysteines were excluded from the structure of SoPCS. The predicted 3D structure of SoPCS showed that it is able to bind to L-gamma-glutamylcysteine as substrate. The binding site sequence of PCS included amino acids 52(Q),55(P),56(A),57(F), 58(C),103(G),104(I),151(S),163(G),165(F),206(D), 213(R). The common amino acid with conserved sequence in the binding site of the plant was 103Gly. Gene expression indicated that SoPCS has an important role in the response of sugarcane to cadmium with potential use in genetic engineering to remove metal contaminants in the environment. This is the first characterization of a PCS from sugarcane.

LIBRA-WA: a web application for ligand binding site detection and protein function recognition

Summary

Recently, LIBRA, a tool for active/ligand binding site prediction, was described. LIBRA's effectiveness was comparable to similar state-of-the-art tools; however, its scoring scheme, output presentation, dependence on local resources and overall convenience were amenable to improvements. To solve these issues, LIBRA-WA, a web application based on an improved LIBRA engine, has been developed, featuring a novel scoring scheme consistently improving LIBRA's performance, and a refined algorithm that can identify binding sites hosted at the interface between different subunits. LIBRA-WA also sports additional functionalities like ligand clustering and a completely redesigned interface for an easier analysis of the output. Extensive tests on 373 apoprotein structures indicate that LIBRA-WA is able to identify the biologically relevant ligand/ligand binding site in 357 cases (∼96%), with the correct prediction ranking first in 349 cases (∼98% of the latter, ∼94% of the total). The earlier stand-alone tool has also been updated and dubbed LIBRA+, by integrating LIBRA-WA's improved engine for cross-compatibility purposes.

Availability and implementation

LIBRA-WA and LIBRA+ are available at: http://www.computationalbiology.it/software.html.

Contact

polticel@uniroma3.it.

Supplementary information

Supplementary data are available at Bioinformatics online.

Crystal structure and functional analysis of human C1ORF123

Proteins of the DUF866 superfamily are exclusively found in eukaryotic cells. A member of the DUF866 superfamily, C1ORF123, is a human protein found in the open reading frame 123 of chromosome 1. The physiological role of C1ORF123 is yet to be determined. The only available protein structure of the DUF866 family shares just 26% sequence similarity and does not contain a zinc binding motif. Here, we present the crystal structure of the recombinant human C1ORF123 protein (rC1ORF123). The structure has a 2-fold internal symmetry dividing the monomeric protein into two mirrored halves that comprise of distinct electrostatic potential. The N-terminal half of rC1ORF123 includes a zinc-binding domain interacting with a zinc ion near to a potential ligand binding cavity. Functional studies of human C1ORF123 and its homologue in the fission yeast Schizosaccharomyces pombe (SpEss1) point to a role of DUF866 protein in mitochondrial oxidative phosphorylation.

Genomic screening of new putative antiviral lectins from Amazonian cyanobacteria based on a bioinformatics approach

Lectins are proteins of nonimmune origin, which are capable of recognizing and binding to glycoconjugate moieties. Some of them can block the interaction of viral glycoproteins to the host cell receptors acting as antiviral agents. Although cyanobacterial lectins have presented broad biotechnological potential, little research has been directed to Amazonian Cyanobacterial diversity. In order to identify new antiviral lectins, we performed genomic analysis in seven cyanobacterial strains from Coleção Amazônica de Cianobactérias e Microalgas (CACIAM). We found 75 unique CDS presenting one or more lectin domains. Since almost all were annotated as hypothetical proteins, we used homology modeling and molecular dynamics simulations to evaluate the structural and functional properties of three CDS that were more similar to known antiviral lectins. Nostoc sp. CACIAM 19 as well as Tolypothrix sp. CACIAM 22 strains presented cyanovirin‐N homologues whose function was confirmed by binding free energy calculations. Asn, Glu, Thr, Lys, Leu, and Gly, which were described as binding residues for cyanovirin, were also observed on those structures. As for other known cyanovirins, those residues in both our models also made favorable interactions with dimannose. Finally, Alkalinema sp. CACIAM 70d presented one CDS, which was identified as a seven‐bladed beta‐propeller structure with binding sites predicted for sialic acid and N‐acetylglucosamine. Despite its singular structure, our analysis suggested this molecule as a new putative antiviral lectin. Overall, the identification and the characterization of new lectins and their homologues are a promising area in antiviral research, and Amazonian cyanobacteria present biotechnological potential to be explored in this regard.

Design and implementation of structural bioinformatics projects for biological sciences undergraduate students

Contemporary biology is currently undergoing a revolution, driven by the availability of high-throughput technologies and a wide variety of bioinformatics tools. However, bioinformatics education and practice is still in its infancy in most of the African continent. Consequently, concerted efforts have been made in recent years to incorporate bioinformatics modules into biological sciences curriculum of African Universities. Despite this, one aspect of bioinformatics that is yet to be incorporated is structural bioinformatics. In this article, we report on a structural bioinformatics project carried out by final year project students in a Nigerian university. The target protein was the thermoacidophilic Sulfolobus islandicus rod-shaped virus 1 (SIRV1) Rep protein, which was further characterized using various free, user-friendly and online sequence-based and structure-based bioinformatics tools. This exercise gave students the opportunity to generate new data, interpret the data, and acquire collaborative research skills. In this report, emphasis is placed on analysis of the data generated to further encourage analytical skills. By sharing this experience, it is anticipated that other similar institutions would adopt parallel strategies to expose undergraduate students to structural biology, and increase awareness of freely available bioinformatics tools for tackling pertinent biological questions. © 2018 International Union of Biochemistry and Molecular Biology, 46(5):547-554, 2018.

The Functional Characterization of Podosphaera xanthii Candidate Effector Genes Reveals Novel Target Functions for Fungal Pathogenicity

Podosphaera xanthii is the main causal agent of powdery mildew disease in cucurbits. In a previous study, we determined that P. xanthii expresses approximately 50 Podosphaera effector candidates (PECs), identified based on the presence of a predicted signal peptide and the absence of functional annotation. In this work, we used host-induced gene silencing (HIGS), employing Agrobacterium tumefaciens as a vector for the delivery of the silencing constructs (ATM-HIGS), to identify genes involved in early plant-pathogen interaction. The analysis of seven selected PEC-encoding genes showed that six of them, PEC007, PEC009, PEC019, PEC032, PEC034, and PEC054, are required for P. xanthii pathogenesis, as revealed by reduced fungal growth and increased production of hydrogen peroxide by host cells. In addition, protein models and protein-ligand predictions allowed us to identify putative functions for these candidates. The biochemical activities of PEC019, PEC032, and PEC054 were elucidated using their corresponding proteins expressed in Escherichia coli. These proteins were confirmed as phospholipid-binding protein, α-mannosidase, and cellulose-binding protein. Further, BLAST searches showed that these three effectors are widely distributed in phytopathogenic fungi. These results suggest novel targets for fungal effectors, such as host-cell plasma membrane, host-cell glycosylation, and damage-associated molecular pattern-triggered immunity.

Screening and analysis of acetyl-cholinesterase (AChE) inhibitors in the context of Alzheimer's disease

Acetyl-cholinesterase enzyme (AChE) is a known target for identifying potential inhibitors against Alzheimer diseases (AD). Therefore, it is of interest to screen AChE with the CNS-BBB database. An AChE enzyme is a member of hydrolase family is activated by acetylcholine (ACh), so, targeting the AChE enzyme with the potential inhibitor may block the binding of the ACh. In this study we carried out virtual screening of drug-like molecules from Chemical Diversity Database particularly CNS-BBB compounds, to identify potential inhibitors using Glide docking program. Top ranking ten compounds, which have lower Glide Score when compared to known drugs (Tacrine and Galantamine) for AChE. For top three molecules MD simulation was carried out and calculated binding free energy. We report the best binding compounds with AChE compared to known drugs (Taine and Galantamine) for AD. We further document the salient features of their molecular interaction with the known target. Three molecules (1-benzyl-3-(2- hydroxyethyl)-N-[2-(3-pyridyl)ethyl]-3-pyrrolidinecarboxamide, N-{3[benzyl(methyl)amino]propyl}-1,5-dimethyl-4-oxo-4,5-dihydro- 1H-pyrrolo[3,2-c]quinoline-2-carboxamide, and 6-chloro-N-[2-(diethylamino)-2-phenylethyl]-4-oxo-4H-chromene-2-carboxamide) have -196.36, -204.27, -214.40 kJ/mol, binding free energy values respectively which are much lower than values calculated for the reference ligands Tacrine and Galantamine having -119.65 and -142.18 kJ/mol respectively. Thus these molecules can be very novel potential inhibitors against AChE involved in Alzheimer's disease.