MOTIPS: Automated Motif Analysis for Predicting Targets of Modular Protein Domains

Background  

Many protein interactions, especially those involved in signaling, involve short linear motifs consisting of 5-10 amino acid residues that interact with modular protein domains such as the SH3 binding domains and the kinase catalytic domains. One straightforward way of identifying these interactions is by scanning for matches to the motif against all the sequences in a target proteome. However, predicting domain targets by motif sequence alone without considering other genomic and structural information has been shown to be lacking in accuracy.     

Differential responsiveness to immunoablative therapy in refractory rheumatoid arthritis is associated with level and avidity of anti-cyclic citrullinated protein autoantibodies: a case study

In order to identify pathogenic correlates of refractory rheumatoid arthritis (RA), antibodies against anti-cyclic citrullinated protein (ACPAs) were investigated in RA patients in whom the dysregulated immune system had been ablated by high-dose chemotherapy (HDC) and autologous haematopoietic stem cell transplantation (HSCT). Six patients with refractory RA were extensively characterized in terms of levels of total immunoglobulins, RA-specific autoantibodies (ACPAs and rheumatoid factor) and antibodies against rubella, tetanus toxoid (TT) and phosphorylcholine before and after HDC plus HSCT. Additionally, the avidity of ACPAs was measured before and after treatment and compared with the avidity of TT antibodies following repeated immunizations. Synovial biopsies were obtained by arthroscopy before HDC plus HSCT, and analyzed by immunohistochemistry. In the three patients with clinically long-lasting responses to HDC plus HSCT (median 423 days), significant reductions in ACPA-IgG levels after therapy were observed (median level dropped from 215 to 34 arbitrary units/ml; P = 0.05). In contrast, stable ACPA-IgG levels were observed in three patients who relapsed shortly after HDC plus HSCT (median of 67 days). Clinical responders had ACPA-IgG of lower avidity (r = 0.75; P = 0.08) and higher degree of inflammation histologically (r = 0.73; P = 0.09). Relapse (after 38 to 530 days) in all patients was preceded by rising levels of low avidity ACPA-IgG (after 30 to 388 days), in contrast to the stable titres of high avidity TT antibodies. In conclusion, humoral autoimmune responses were differentially modulated by immunoablative therapy in patients with synovial inflammation and low avidity ACPA-IgG autoantibodies as compared with patients with high levels of high avidity ACPA-IgG. The distinct clinical disease course after immunoablative therapy based on levels and avidity of ACPA-IgG indicates that refractory RA is not a single disease entity.     

Downregulation of the Host Gene jigr1 by miR-92 Is Essential for Neuroblast Self-Renewal in Drosophila

Intragenic microRNAs (miRNAs), located mostly in the introns of protein-coding genes, are often co-expressed with their host mRNAs. However, their functional interaction in development is largely unknown. Here we show that in Drosophila, miR-92a and miR-92b are embedded in the intron and 3’UTR of jigr1, respectively, and co-expressed with some jigr1 isoforms. miR-92a and miR-92b are highly expressed in neuroblasts of larval brain where Jigr1 expression is low. Genetic deletion of both miR-92a and miR-92b demonstrates an essential cell-autonomous role for these miRNAs in maintaining neuroblast self-renewal through inhibiting premature differentiation. We also show that miR-92a and miR-92b directly target jigr1 in vivo and that some phenotypes due to the absence of these miRNAs are partially rescued by reducing the level of jigr1. These results reveal a novel function of the miR-92 family in Drosophila neuroblasts and provide another example that local negative feedback regulation of host genes by intragenic miRNAs is essential for animal development.     

Electrochemical based detection of microRNA, mir21 in breast cancer cells

Hemagglutinin (HA) is a trimeric glycoprotein expressed on the influenza virus membrane. HA of influenza viruses binds to the host's cell surface complex glycans via a terminal sialic acid (Sia), as the first key step in the process of infection, transmission and virulence of influenza viruses. It is important to monitor and evaluate the receptor (glycan) binding preferences of the HAs derived from influenza A viruses, especially those originating from birds and swine, to understand their potential ability for interspecies transmission. From this viewpoint, in the present study, we have developed a protocol for analyzing the glycan-HA interactions efficiently and kinetically, based on surface plasmon resonance (SPR). Our results showed that glycan-HA binding analyses can be performed reliably and efficiently on Biacore-chips in the SPR system, using chemically synthesized biotinylated multivalent-glycans. Using the CAP-chip, we were able to regenerate the surface for multiple analyses, allowing us to derive, for the first time, the precise kinetic parameters for different HA-glycan complexes of newly emerging influenza viruses. These studies suggested that this SPR-based method is suitable for influenza surveillance to define the pandemic scenario as well as to screen of synthetic glycans and other compounds that may interfere with glycan-HA interactions.     

Ghrelin and obestatin expression in oral squamous cell carcinoma: an immunohistochemical and biochemical study

The underlying molecular mechanism of carcinogenesis in oral squamous cell carcinoma (OSCC) is poorly understood and appears to be controlled on many genetic, environmental, and hormonal factors. Obestatin and ghrelin, two recently discovered hormones, are co-expressed in endocrine cells. The purpose of this investigation was to examine the immunohistochemical features of OSCCs in relation to the tissue concentration of ghrelin and obestatin. The association between OSCC and Epstein Barr Virus (EBV) status was also explored. The expression of ghrelin and obestatin was examined by immunohistochemistry and immunoassay in oral biopsy specimens: 10 benign squamous epithelial cell samples, 10 microinvasive squamous cell carcinomas, and seven well-differentiated and seven poorly differentiated OSCCs. The presence of EBV was evaluated in these samples using immunohistochemistry. The concentrations of ghrelin and obestatin in tissue homogenates were measured by RIA and ELISA, respectively. Squamous cell carcinomas and benign tissue samples were positive for anti-EBV antibody, and obestatin and ghrelin were shown to be co-expressed in all stratified squamous epithelium samples. Expression of ghrelin and obestatin was decreased or absent in OSCCs in relation to the invasiveness of the carcinoma; ghrelin and obestatin levels in cancerous tissue homogenates were lower than in benign tissue homogenates. These results indicate that the concentrations and distribution of immunoreactive obestatin and ghrelin might be helpful in distinguishing OSCC from benign tumors. Maintaining normal levels of these hormones might be required for regulation of normal cell division. However, detailed studies will be required for better understanding of the complex mechanism of carcinogenesis relating to OSCCs.     

A flexible docking scheme to explore the binding selectivity of PDZ domains

Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTA LIGAND , we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density‐95/Dlg/ZO‐1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 Å. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately.     

ProbFAST: Probabilistic Functional Analysis System Tool

Background  

The post-genomic era has brought new challenges regarding the understanding of the organization and function of the human genome. Many of these challenges are centered on the meaning of differential gene regulation under distinct biological conditions and can be performed by analyzing the Multiple Differential Expression (MDE) of genes associated with normal and abnormal biological processes. Currently MDE analyses are limited to usual methods of differential expression initially designed for paired analysis.     

Blind Test of Physics-Based Prediction of Protein Structures

We report here a multiprotein blind test of a computer method to predict native protein structures based solely on an all-atom physics-based force field. We use the AMBER 96 potential function with an implicit (GB/SA) model of solvation, combined with replica-exchange molecular-dynamics simulations. Coarse conformational sampling is performed using the zipping and assembly method (ZAM), an approach that is designed to mimic the putative physical routes of protein folding. ZAM was applied to the folding of six proteins, from 76 to 112 monomers in length, in CASP7, a community-wide blind test of protein structure prediction. Because these predictions have about the same level of accuracy as typical bioinformatics methods, and do not utilize information from databases of known native structures, this work opens up the possibility of predicting the structures of membrane proteins, synthetic peptides, or other foldable polymers, for which there is little prior knowledge of native structures. This approach may also be useful for predicting physical protein folding routes, non-native conformations, and other physical properties from amino acid sequences.     

Describing the current status of Plasmodium falciparum population structure and drug resistance within mainland Tanzania using molecular inversion probes

High‐throughput Plasmodium genomic data is increasingly useful in assessing prevalence of clinically important mutations and malaria transmission patterns. Understanding parasite diversity is important for identification of specific human or parasite populations that can be targeted by control programmes, and to monitor the spread of mutations associated with drug resistance. An up‐to‐date understanding of regional parasite population dynamics is also critical to monitor the impact of control efforts. However, this data is largely absent from high‐burden nations in Africa, and to date, no such analysis has been conducted for malaria parasites in Tanzania countrywide. To this end, over 1,000 P. falciparum clinical isolates were collected in 2017 from 13 sites in seven administrative regions across Tanzania, and parasites were genotyped at 1,800 variable positions genome‐wide using molecular inversion probes. Population structure was detectable among Tanzanian P. falciparum parasites, approximately separating parasites from the northern and southern districts and identifying genetically admixed populations in the north. Isolates from nearby districts were more likely to be genetically related compared to parasites sampled from more distant districts. Known drug resistance mutations were seen at increased frequency in northern districts (including two infections carrying pfk13‐R561H), and additional variants with undetermined significance for antimalarial resistance also varied by geography. Malaria Indicator Survey (2017) data corresponded with genetic findings, including average region‐level complexity‐of‐infection and malaria prevalence estimates. The parasite populations identified here provide important information on extant spatial patterns of genetic diversity of Tanzanian parasites, to which future surveys of genetic relatedness can be compared.     

Cloning and expression of the Clostridium thermocellum L-lactate dehydrogenase gene in Escherichia coli and enzyme characterization

The structural gene for L-lactate dehydrogenase (LDH) (EC.1.1.1.27) from Clostridium thermocellum 27405 was cloned in Escherichia coli by screening the Lambda Zap II phage library of C. thermocellum genomic DNA. In one positive clone, an open reading frame of 948 base pairs corresponded to C. thermocellum ldh gene encoding for the predicted 315-residue protein. The ldh gene was successfully expressed in E. coli FMJ39 (ldh mutant) under the lac promoter. The recombinant enzyme was partially purified from E. coli cell extracts and its kinetic properties were determined. Clostridium thermocellum LDH was shown to catalyze a highly reversible reaction and to be an allosteric enzyme that is activated by fructose-1,6-diphosphate (FDP). For pyruvate, partially purified LDH had Km and Vmax values of 7.3 mmol/L and 87 micromol/min, respectively, and in the presence of FDP, a 24-fold decrease in Km and a 5.7-fold increase in Vmax were recorded. The enzyme exhibited no marked catalytic activity for lactate in the absence of FDP, whereas Km and Vmax values were 59.5 mmol/L and 52 micromol/min, respectively, in its presence. The enzyme did not lose activity when incubated at 65 degrees C for 5 min.