Treatment with a Cholecystokinin Receptor Antagonist, Proglumide, Improves Efficacy of Immune Checkpoint Antibodies in Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-associated deaths worldwide. Treatment with immune checkpoint antibodies has shown promise in advanced HCC, but the response is only 15-20%. We discovered a potential target for the treatment of HCC, the cholecystokinin-B receptor (CCK-BR). This receptor is overexpressed in murine and human HCC and not in normal liver tissue. Mice bearing syngeneic RIL-175 HCC tumors were treated with phosphate buffer saline (PBS; control), proglumide (a CCK-receptor antagonist), an antibody to programmed cell death protein 1 (PD-1Ab), or the combination of proglumide and the PD-1Ab. In vitro, RNA was extracted from untreated or proglumide-treated murine Dt81Hepa1-6 HCC cells and analyzed for expression of fibrosis-associated genes. RNA was also extracted from human HepG2 HCC cells or HepG2 cells treated with proglumide and subjected to RNA sequencing. Results showed that proglumide decreased fibrosis in the tumor microenvironment and increased the number of intratumoral CD8+ T cells in RIL-175 tumors. When proglumide was given in combination with the PD-1Ab, there was a further significant increase in intratumoral CD8+ T cells, improved survival, and alterations in genes regulating tumoral fibrosis and epithelial-to-mesenchymal transition. RNAseq results from human HepG2 HCC cells treated with proglumide showed significant changes in differentially expressed genes involved in tumorigenesis, fibrosis, and the tumor microenvironment. The use of the CCK receptor antagonist may improve efficacy of immune checkpoint antibodies and survival in those with advanced HCC.

Understanding COVID-19 Pathogenesis: A Drug-Repurposing Effort to Disrupt Nsp-1 Binding to Export Machinery Receptor Complex.

Non-structural protein 1 (Nsp1) is a virulence factor found in all beta coronaviruses (b-CoVs). Recent studies have shown that Nsp1 of SARS-CoV-2 virus interacts with the nuclear export receptor complex, which includes nuclear RNA export factor 1 (NXF1) and nuclear transport factor 2-like export factor 1 (NXT1). The NXF1-NXT1 complex plays a crucial role in the transport of host messenger RNA (mRNA). Nsp1 interferes with the proper binding of NXF1 to mRNA export adaptors and its docking to the nuclear pore complex. We propose that drugs targeting the binding surface between Nsp1 and NXF1-NXT1 may be a useful strategy to restore host antiviral gene expression. Exploring this strategy forms the main goals of this paper. Crystal structures of Nsp1 and the heterodimer of NXF1-NXT1 have been determined. We modeled the docking of Nsp1 to the NXF1-NXT1 complex, and discovered repurposed drugs that may interfere with this binding. To our knowledge, this is the first attempt at drug-repurposing of this complex. We used structural analysis to screen 1993 FDA-approved drugs for docking to the NXF1-NXT1 complex. The top hit was ganirelix, with a docking score of -14.49. Ganirelix competitively antagonizes the gonadotropin releasing hormone receptor (GNRHR) on pituitary gonadotrophs, and induces rapid, reversible suppression of gonadotropin secretion. The conformations of Nsp1 and GNRHR make it unlikely that they interact with each other. Additional drug leads were inferred from the structural analysis of this complex, which are discussed in the paper. These drugs offer several options for therapeutically blocking Nsp1 binding to NFX1-NXT1, which may normalize nuclear export in COVID-19 infection.

Host gene expression modulated by Zika virus infection of human-293 cells.

The HEK-293 cell line was created in 1977 by transformation of primary human embryonic kidney cells with sheared adenovirus type 5 DNA. A previous study determined that the HEK-293 cells have neuronal markers rather than kidney markers. In this study, we tested the hypothesis whether Zika virus (ZIKV), a neurotropic virus, is able to infect and replicate in the HEK-293 cells. We show that the HEK-293 cells infected with ZIKV support viral replication as shown by indirect immunofluorescence (IFA) and quantitative reverse transcriptase-PCR (qRT-PCR). We performed RNA-seq analysis on the ZIKV-infected and the control uninfected HEK-293 cells and find 659 genes that are differentially transcribed in ZIKV-infected HEK-293 cells as compared to uninfected cells. The results show that the top 10 differentially transcribed and upregulated genes are involved in antiviral and inflammatory responses. Seven upregulated genes, IFNL1, DDX58, CXCL10, ISG15, KCNJ15, IFNIH1, and IFIT2, were validated by qRT-PCR. Altogether, our findings show that ZIKV infection alters host gene expression by affecting their antiviral and inflammatory responses.

Network and Systems Medicine: Position Paper of the European Collaboration on Science and Technology Action on Open Multiscale Systems Medicine.

Introduction: Network and systems medicine has rapidly evolved over the past decade, thanks to computational and integrative tools, which stem in part from systems biology. However, major challenges and hurdles are still present regarding validation and translation into clinical application and decision making for precision medicine. Methods: In this context, the Collaboration on Science and Technology Action on Open Multiscale Systems Medicine (OpenMultiMed) reviewed the available advanced technologies for multidimensional data generation and integration in an open-science approach as well as key clinical applications of network and systems medicine and the main issues and opportunities for the future. Results: The development of multi-omic approaches as well as new digital tools provides a unique opportunity to explore complex biological systems and networks at different scales. Moreover, the application of findable, applicable, interoperable, and reusable principles and the adoption of standards increases data availability and sharing for multiscale integration and interpretation. These innovations have led to the first clinical applications of network and systems medicine, particularly in the field of personalized therapy and drug dosing. Enlarging network and systems medicine application would now imply to increase patient engagement and health care providers as well as to educate the novel generations of medical doctors and biomedical researchers to shift the current organ- and symptom-based medical concepts toward network- and systems-based ones for more precise diagnoses, interventions, and ideally prevention. Conclusion: In this dynamic setting, the health care system will also have to evolve, if not revolutionize, in terms of organization and management.

Recent Advances in Systems and Network Medicine: Meeting Report from the First International Conference in Systems and Network Medicine.

The First International Conference in Systems and Network Medicine gathered together 200 global thought leaders, scientists, clinicians, academicians, industry and government experts, medical and graduate students, postdoctoral scholars and policymakers. Held at Georgetown University Conference Center in Washington D.C. on September 11-13, 2019, the event featured a day of pre-conference lectures and hands-on bioinformatic computational workshops followed by two days of deep and diverse scientific talks, panel discussions with eminent thought leaders, and scientific poster presentations. Topics ranged from: Systems and Network Medicine in Clinical Practice; the role of -omics technologies in Health Care; the role of Education and Ethics in Clinical Practice, Systems Thinking, and Rare Diseases; and the role of Artificial Intelligence in Medicine. The conference served as a unique nexus for interdisciplinary discovery and dialogue and fostered formation of new insights and possibilities for health care systems advances.

Ridge regression estimated linear probability model predictions of O-glycosylation in proteins with structural and sequence data.

BACKGROUND: To-date, no claim regarding finding a consensus sequon for O-glycosylation has been made. Thus, predicting the likelihood of O-glycosylation with sequence and structural information using classical regression analysis is quite difficult. In particular, if a binary response is used to distinguish between O-glycosylated and non-O-glycosylated sequences, an appropriate set of non-O-glycosylatable sequences is hard to find. RESULTS: Three sequences from similar post-translational modifications (PTMs) of proteins occurring at, or very near, the S/T-site are analyzed: N-glycosylation, O-mucin type (O-GalNAc) glycosylation, and phosphorylation. Results found include: 1) The consensus composite sequon for O-glycosylation is: ~(W-S/T-W), where "~" denotes the "not" operator. 2) The consensus sequon for phosphorylation is ~(W-S/T/Y/H-W); although W-S/T/Y/H-W is not an absolute inhibitor of phosphorylation. 3) For linear probability model (LPM) estimation, N-glycosylated sequences are good approximations to non-O-glycosylatable sequences; although N - ~P - S/T is not an absolute inhibitor of O-glycosylation. 4) The selective positioning of an amino acid along the sequence, differentiates the PTMs of proteins. 5) Some N-glycosylated sequences are also phosphorylated at the S/T-site in the N - ~P - S/T sequon. 6) ASA values for N-glycosylated sequences are stochastically larger than those for O-GlcNAc glycosylated sequences. 7) Structural attributes (beta turn II, II´, helix, beta bridges, beta hairpin, and the phi angle) are significant LPM predictors of O-GlcNAc glycosylation. The LPM with sequence and structural data as explanatory variables yields a Kolmogorov-Smirnov (KS) statistic of 99%. 8) With only sequence data, the KS statistic erodes to 80%, and 21% of out-of-sample O-GlcNAc glycosylated sequences are mispredicted as not being glycosylated. The 95% confidence interval around this mispredictions rate is 16% to 26%. CONCLUSIONS: The data indicates the existence of a consensus sequon for O-glycosylation; and underscores the germaneness of structural information for predicting the likelihood of O-glycosylation.

Nucleotide-Induced Conformational Changes in Escherichia coli DnaA Protein Are Required for Bacterial ORC to Pre-RC Conversion at the Chromosomal Origin.

DnaA oligomerizes when bound to origins of chromosomal replication. Structural analysis of a truncated form of DnaA from Aquifex aeolicus has provided insight into crucial conformational differences within the AAA+ domain that are specific to the ATP- versus ADP- bound form of DnaA. In this study molecular docking of ATP and ADP onto Escherichia coli DnaA, modeled on the crystal structure of Aquifex aeolicus DnaA, reveals changes in the orientation of amino acid residues within or near the vicinity of the nucleotide-binding pocket. Upon limited proteolysis with trypsin or chymotrypsin ADP-DnaA, but not ATP-DnaA generated relatively stable proteolytic fragments of various sizes. Examined sites of limited protease susceptibility that differ between ATP-DnaA and ADP-DnaA largely reside in the amino terminal half of DnaA. The concentration of adenine nucleotide needed to induce conformational changes, as detected by these protease susceptibilities of DnaA, coincides with the conversion of an inactive bacterial origin recognition complex (bORC) to a replication efficient pre-replication complex (pre-RC) at the E. coli chromosomal origin of replication (oriC).

Naturally occurring polyphenol, morin hydrate, inhibits enzymatic activity of N-methylpurine DNA glycosylase, a DNA repair enzyme with various roles in human disease.

Interest in the mechanisms of DNA repair pathways, including the base excision repair (BER) pathway specifically, has heightened since these pathways have been shown to modulate important aspects of human disease. Modulation of the expression or activity of a particular BER enzyme, N-methylpurine DNA glycosylase (MPG), has been demonstrated to play a role in carcinogenesis and resistance to chemotherapy as well as neurodegenerative diseases, which has intensified the focus on studying MPG-related mechanisms of repair. A specific small molecule inhibitor for MPG activity would be a valuable biochemical tool for understanding these repair mechanisms. By screening several small molecule chemical libraries, we identified a natural polyphenolic compound, morin hydrate, which inhibits MPG activity specifically (IC50=2.6µM). Detailed mechanism analysis showed that morin hydrate inhibited substrate DNA binding of MPG, and eventually the enzymatic activity of MPG. Computational docking studies with an x-ray derived MPG structure as well as comparison studies with other structurally-related flavonoids offer a rationale for the inhibitory activity of morin hydrate observed. The results of this study suggest that the morin hydrate could be an effective tool for studying MPG function and it is possible that morin hydrate and its derivatives could be utilized in future studies focused on the role of MPG in human disease.

Structural and functional studies of S-adenosyl-L-methionine binding proteins: a ligand-centric approach.

BACKGROUND: The post-genomic era poses several challenges. The biggest is the identification of biochemical function for protein sequences and structures resulting from genomic initiatives. Most sequences lack a characterized function and are annotated as hypothetical or uncharacterized. While homology-based methods are useful, and work well for sequences with sequence identities above 50%, they fail for sequences in the twilight zone (<30%) of sequence identity. For cases where sequence methods fail, structural approaches are often used, based on the premise that structure preserves function for longer evolutionary time-frames than sequence alone. It is now clear that no single method can be used successfully for functional inference. Given the growing need for functional assignments, we describe here a systematic new approach, designated ligand-centric, which is primarily based on analysis of ligand-bound/unbound structures in the PDB. Results of applying our approach to S-adenosyl-L-methionine (SAM) binding proteins are presented. RESULTS: Our analysis included 1,224 structures that belong to 172 unique families of the Protein Information Resource Superfamily system. Our ligand-centric approach was divided into four levels: residue, protein/domain, ligand, and family levels. The residue level included the identification of conserved binding site residues based on structure-guided sequence alignments of representative members of a family, and the identification of conserved structural motifs. The protein/domain level included structural classification of proteins, Pfam domains, domain architectures, and protein topologies. The ligand level included ligand conformations, ribose sugar puckering, and the identification of conserved ligand-atom interactions. The family level included phylogenetic analysis. CONCLUSION: We found that SAM bound to a total of 18 different fold types (I-XVIII). We identified 4 new fold types and 11 additional topological arrangements of strands within the well-studied Rossmann fold Methyltransferases (MTases). This extends the existing structural classification of SAM binding proteins. A striking correlation between fold type and the conformation of the bound SAM (classified as types) was found across the 18 fold types. Several site-specific rules were created for the assignment of functional residues to families and proteins that do not have a bound SAM or a solved structure.

Recent advances in biocuration: meeting report from the Fifth International Biocuration Conference.

The 5th International Biocuration Conference brought together over 300 scientists to exchange on their work, as well as discuss issues relevant to the International Society for Biocuration's (ISB) mission. Recurring themes this year included the creation and promotion of gold standards, the need for more ontologies, and more formal interactions with journals. The conference is an essential part of the ISB's goal to support exchanges among members of the biocuration community. Next year's conference will be held in Cambridge, UK, from 7 to 10 April 2013. In the meanwhile, the ISB website provides information about the society's activities (http://biocurator.org), as well as related events of interest.

Structure-function analysis of human TYW2 enzyme required for the biosynthesis of a highly modified Wybutosine (yW) base in phenylalanine-tRNA.

Posttranscriptional modifications are critical for structure and function of tRNAs. Wybutosine (yW) and its derivatives are hyper-modified guanosines found at the position 37 of eukaryotic and archaeal tRNA(Phe). TYW2 is an enzyme that catalyzes α-amino-α-carboxypropyl transfer activity at the third step of yW biogenesis. Using complementation of a ΔTYW2 strain, we demonstrate here that human TYW2 (hTYW2) is active in yeast and can synthesize the yW of yeast tRNA(Phe). Structure-guided analysis identified several conserved residues in hTYW2 that interact with S-adenosyl-methionine (AdoMet), and mutation studies revealed that K225 and E265 are critical residues for the enzymatic activity. We previously reported that the human TYW2 is overexpressed in breast cancer. However, no difference in the tRNA(Phe) modification status was observed in either normal mouse tissue or a mouse tumor model that overexpresses Tyw2, indicating that hTYW2 may have a role in tumorigenesis unrelated to yW biogenesis.

Proteome-wide analysis of single-nucleotide variations in the N-glycosylation sequon of human genes.

N-linked glycosylation is one of the most frequent post-translational modifications of proteins with a profound impact on their biological function. Besides other functions, N-linked glycosylation assists in protein folding, determines protein orientation at the cell surface, or protects proteins from proteases. The N-linked glycans attach to asparagines in the sequence context Asn-X-Ser/Thr, where X is any amino acid except proline. Any variation (e.g. non-synonymous single nucleotide polymorphism or mutation) that abolishes the N-glycosylation sequence motif will lead to the loss of a glycosylation site. On the other hand, variations causing a substitution that creates a new N-glycosylation sequence motif can result in the gain of glycosylation. Although the general importance of glycosylation is well known and acknowledged, the effect of variation on the actual glycoproteome of an organism is still mostly unknown. In this study, we focus on a comprehensive analysis of non-synonymous single nucleotide variations (nsSNV) that lead to either loss or gain of the N-glycosylation motif. We find that 1091 proteins have modified N-glycosylation sequons due to nsSNVs in the genome. Based on analysis of proteins that have a solved 3D structure at the site of variation, we find that 48% of the variations that lead to changes in glycosylation sites occur at the loop and bend regions of the proteins. Pathway and function enrichment analysis show that a significant number of proteins that gained or lost the glycosylation motif are involved in kinase activity, immune response, and blood coagulation. A structure-function analysis of a blood coagulation protein, antithrombin III and a protease, cathepsin D, showcases how a comprehensive study followed by structural analysis can help better understand the functional impact of the nsSNVs.

Structure-guided rule-based annotation of protein functional sites in UniProt knowledgebase.

The rapid growth of protein sequence databases has necessitated the development of methods to computationally derive annotation for uncharacterized entries. Most such methods focus on "global" annotation, such as molecular function or biological process. Methods to supply high-accuracy "local" annotation to functional sites based on structural information at the level of individual amino acids are relatively rare. In this chapter we will describe a method we have developed for annotation of functional residues within experimentally-uncharacterized proteins that relies on position-specific site annotation rules (PIR Site Rules) derived from structural and experimental information. These PIR Site Rules are manually defined to allow for conditional propagation of annotation. Each rule specifies a tripartite set of conditions whereby candidates for annotation must pass a whole-protein classification test (that is, have end-to-end match to a whole-protein-based HMM), match a site-specific profile HMM and, finally, match functionally and structurally characterized residues of a template. Positive matches trigger the appropriate annotation for active site residues, binding site residues, modified residues, or other functionally important amino acids. The strict criteria used in this process have rendered high-confidence annotation suitable for UniProtKB/Swiss-Prot features.

Sequence signatures in envelope protein may determine whether flaviviruses produce hemorrhagic or encephalitic syndromes.

We analyzed the envelope proteins in pathogenic flaviviruses to determine whether there are sequence signatures associated with the tendency of viruses to produce hemorrhagic disease (H-viruses) or encephalitis (E-viruses). We found that, at the position corresponding to the glycosylated Asn-67 in dengue virus, asparagine (Asn) occurs in all seven viral species that cause hemorrhagic disease in humans. Furthermore, Asn was extremely rare at position 67 in six flaviviruses that cause encephalitis, being replaced by Asp in four of them. Of the 3,246 sequences from H- and E-viruses, we found that 2,916 sequences (90%) contained Asn in position 67 for H-viruses or Asp in position 67 for E-viruses. The change from Asn-67 that is prevalent in H-viruses to Asp-67 (common in E-viruses) contributes to a stronger electrostatically negative surface in the E-viruses as compared to the H-viruses. These findings should help predicting the disease potential of emerging and re-emerging flaviviruses and understanding the relationship between protein structure and disease outcome.

Protein functional annotation by homology.

Genome sequencing projects have resulted in a rapid accumulation of predicted protein sequences. With experimentally verified information on protein function lagging far behind, computational methods are used for functional annotation of proteins. Here we describe a number of protocols for protein sequence and structure analysis that can be used to infer function of uncharacterized proteins. These protocols rely on publicly available computational resources and tools and can be utilized by anyone with an Internet access.

Intestinal epithelial CD98: an oligomeric and multifunctional protein.

The intestinal epithelial cell-surface molecule, CD98 is a type II membrane glycoprotein. Molecular orientation studies have demonstrated that the C-terminal tail of human CD98 (hCD98), which contains a PDZ-binding domain, is extracellular. In intestinal epithelial cells, CD98 is covalently linked to an amino-acid transporter with which it forms a heterodimer. This heterodimer associates with beta(1)-integrin and intercellular adhesion molecular 1 (ICAM-1) to form a macromolecular complex in the basolateral membranes of polarized intestinal epithelial cells. This review focuses on the multifunctional roles of CD98, including involvement in extracellular signaling, adhesion/polarity, and amino-acid transporter expression in intestinal epithelia. A role for CD98 in intestinal inflammation, such as Intestinal Bowel Disease (IBD), is also proposed.

Extracellular interaction between hCD98 and the PDZ class II domain of hCASK in intestinal epithelia.

The extracellular domain of the glycoprotein-associated integrin hCD98 protrudes into the basolateral extracellular space of the intestine and contains a PDZ class II-binding domain (GLLLRFPYAA, amino acids 520-529). Protein-protein interaction studies in vitro as well as in human colonic sections and Caco2-BBE cells have revealed that hCD98 coimmunoprecipitated with the basolateral membrane-associated guanylate kinase hCASK and that this interaction occurred in a PDZ domain-dependent manner. These novel results, which provide the first evidence for a PDZ domain-dependent interaction between a membrane protein and an extracellular protein, open a new field of investigation related to extracellular signaling in cell biology.

A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes.

BACKGROUND: Sequencing the genomes of multiple, taxonomically diverse eukaryotes enables in-depth comparative-genomic analysis which is expected to help in reconstructing ancestral eukaryotic genomes and major events in eukaryotic evolution and in making functional predictions for currently uncharacterized conserved genes. RESULTS: We examined functional and evolutionary patterns in the recently constructed set of 5,873 clusters of predicted orthologs (eukaryotic orthologous groups or KOGs) from seven eukaryotic genomes: Caenorhabditis elegans, Drosophila melanogaster, Homo sapiens, Arabidopsis thaliana, Saccharomyces cerevisiae, Schizosaccharomyces pombe and Encephalitozoon cuniculi. Conservation of KOGs through the phyletic range of eukaryotes strongly correlates with their functions and with the effect of gene knockout on the organism's viability. The approximately 40% of KOGs that are represented in six or seven species are enriched in proteins responsible for housekeeping functions, particularly translation and RNA processing. These conserved KOGs are often essential for survival and might approximate the minimal set of essential eukaryotic genes. The 131 single-member, pan-eukaryotic KOGs we identified were examined in detail. For around 20 that remained uncharacterized, functions were predicted by in-depth sequence analysis and examination of genomic context. Nearly all these proteins are subunits of known or predicted multiprotein complexes, in agreement with the balance hypothesis of evolution of gene copy number. Other KOGs show a variety of phyletic patterns, which points to major contributions of lineage-specific gene loss and the 'invention' of genes new to eukaryotic evolution. Examination of the sets of KOGs lost in individual lineages reveals co-elimination of functionally connected genes. Parsimonious scenarios of eukaryotic genome evolution and gene sets for ancestral eukaryotic forms were reconstructed. The gene set of the last common ancestor of the crown group consists of 3,413 KOGs and largely includes proteins involved in genome replication and expression, and central metabolism. Only 44% of the KOGs, mostly from the reconstructed gene set of the last common ancestor of the crown group, have detectable homologs in prokaryotes; the remainder apparently evolved via duplication with divergence and invention of new genes. CONCLUSIONS: The KOG analysis reveals a conserved core of largely essential eukaryotic genes as well as major diversification and innovation associated with evolution of eukaryotic genomes. The results provide quantitative support for major trends of eukaryotic evolution noticed previously at the qualitative level and a basis for detailed reconstruction of evolution of eukaryotic genomes and biology of ancestral forms.