In complex with the cosubstrate UDP-N-acetylglucosamine (UDP-GlcNAc),O-linked-GlcNAc transferase (OGT) catalyzes Ser/ThrO-GlcNAcylation of many cellular proteins and proteolysis of the transcriptional coregulator HCF-1. Such a dual glycosyltransferase-protease activity, which occurs in the same active site, is unprecedented and integrates both reversible and irreversible forms of protein post-translational modification within one enzyme. Although occurring within the same active site, we show here that glycosylation and proteolysis occur through separable mechanisms. OGT consists of tetratricopeptide repeat (TPR) and catalytic domains, which, together with UDP-GlcNAc, are required for both glycosylation and proteolysis. Nevertheless, a specific TPR domain contact with the HCF-1 substrate is critical for proteolysis but not Ser/Thr glycosylation. In contrast, key catalytic domain residues and even a UDP-GlcNAc oxygen important for Ser/Thr glycosylation are irrelevant for proteolysis. Thus, from a dual glycosyltransferase-protease, essentially single-activity enzymes can be engineered both in vitro and in vivo. Curiously, whereas OGT-mediated HCF-1 proteolysis is limited to vertebrate species, invertebrate OGTs can cleave human HCF-1. We present a model for the evolution of HCF-1 proteolysis by OGT.

The SwissBioisostere database (http://www.swissbioisostere.ch) contains information on molecular replacements and their performance in biochemical assays. It is meant to provide researchers in drug discovery projects with ideas for bioisosteric modifications of their current lead molecule, as well as to give interested scientists access to the details on particular molecular replacements. As of August 2012, the database contains 21,293,355 datapoints corresponding to 5,586,462 unique replacements that have been measured in 35,039 assays against 1948 molecular targets representing 30 target classes. The accessible data were created through detection of matched molecular pairs and mining bioactivity data in the ChEMBL database. The SwissBioisostere database is hosted by the Swiss Institute of Bioinformatics and available via a web-based interface.

Hepatoblastoma is the most common malignant pediatric liver cancer. Histological evaluation of tumor biopsies is used to distinguish among the different subtypes of hepatoblastoma, with fetal and embryonal representing the two main epithelial components. With frequent CTNNB1 mutations, hepatoblastoma is a Wnt/β-catenin-driven malignancy. Considering that Wnt activation has been associated with tumor metabolic reprogramming, we characterized the metabolic profile of cells from hepatoblastoma and compared it to cells from hepatocellular carcinoma. First, we demonstrated that glucose transporter GLUT3 is a direct TCF4/β-catenin target gene. RNA sequencing enabled to identify molecular and metabolic features specific to hepatoblastoma and revealed that several glycolytic enzymes are overexpressed in embryonal-like compared to fetal-like tumor cells. This led us to implement successfully three biomarkers to distinguish embryonal from fetal components by immunohistochemistry from a large panel of human hepatoblastoma samples. Functional analyses demonstrated that embryonal-like hepatoblastoma cells are highly glycolytic and sensitive to hexokinase-1 silencing. Altogether, our findings reveal a new, metabolic classification of human hepatoblastoma, with potential future implications for patients' diagnosis and treatment.

Recent clinical developments in antitumor immunotherapy involving T-cell related therapeutics have led to a renewed interest for human leukocyte antigen class I (HLA-I) binding peptides, given their potential use as peptide vaccines. Databases of HLA-I binding peptides hold therefore information on therapeutic targets essential for understanding immunity. In this work, we use in depth and accurate HLA-I peptidomics datasets determined by mass-spectrometry (MS) and analyze properties of the HLA-I binding peptides with structure-based computational approaches. HLA-I binding peptides are studied grouping all alleles together or in allotype-specific contexts. We capitalize on the increasing number of structurally determined proteins to (1) map the 3D structure of HLA-I binding peptides into the source proteins for analyzing their secondary structure and solvent accessibility in the protein context, and (2) search for potential differences between these properties in HLA-I binding peptides and in a reference dataset of HLA-I motif-like peptides. This is performed by an in-house developed heuristic search that considers peptides across all the human proteome and converges to a collection of peptides that exhibit exactly the same motif as the HLA-I peptides. Our results, based on 9-mers matched to protein 3D structures, clearly show enriched sampling for HLA-I presentation of helical fragments in the source proteins. This enrichment is significant, as compared to 9-mer HLA-I motif-like peptides, and is not entirely explained by the helical propensity of the preferred residues in the HLA-I motifs. We give possible hypothesis for the secondary structure biases observed in HLA-I peptides. This contribution is of potential interest for researchers working in the field of antigen presentation and proteolysis. This knowledge refines the understanding of the rules governing antigen presentation and could be added to the parameters of the current peptide-MHC class I binding predictors to increase their antigen predictive ability.

The CD8 T cell response to the HLA-A2-restricted epitope LLWNGPMAV (LLW) of the non-structural protein 4b of Yellow Fever Virus (YFV) is remarkably immunodominant, highly prevalent and powerful in YFV-vaccinated humans. Here we used a combinatorial peptide library screening in the context of an A2/LLW-specific CD8 T cell clone to identify a superagonist that features a methionine to isoleucine substitution at position 7. Based on in silico modeling, the functional enhancement of this LLW-7I mutation was associated with alterations in the structural dynamics of the peptide in the major histocompatibility complex (pMHC) binding with the T cell receptor (TCR). While the TCR off-rate of LLW-7I pMHC is comparable to the wild type peptide, the rigidity of the 7I peptide seems to confer less entropy loss upon TCR binding. This LLW-7I superagonist is an example of improved functionality in human CD8 T cells associated with optimized ligand rigidity for TCR binding and not with changes in TCR:pMHC off-rate kinetics.

Cargo transport along the cytoplasmic microtubular network is essential for neuronal function, and cytoplasmic dynein-1 is an established molecular motor that is critical for neurogenesis and homeostasis. We performed whole-exome sequencing, homozygosity mapping, and chromosomal microarray studies in five individuals from three independent pedigrees and identified likely-pathogenic variants in DYNC1I2 (Dynein Cytoplasmic 1 Intermediate Chain 2), encoding a component of the cytoplasmic dynein 1 complex. In a consanguineous Pakistani family with three affected individuals presenting with microcephaly, severe intellectual disability, simplification of cerebral gyration, corpus callosum hypoplasia, and dysmorphic facial features, we identified a homozygous splice donor site variant (GenBank: NM_001378.2:c.607+1G>A). We report two additional individuals who have similar neurodevelopmental deficits and craniofacial features and harbor deleterious variants; one individual bears a c.740A>G (p.Tyr247Cys) change in trans with a 374 kb deletion encompassing DYNC1I2, and an unrelated individual harbors the compound-heterozygous variants c.868C>T (p.Gln290∗) and c.740A>G (p.Tyr247Cys). Zebrafish larvae subjected to CRISPR-Cas9 gene disruption or transient suppression of dync1i2a displayed significantly altered craniofacial patterning with concomitant reduction in head size. We monitored cell death and cell cycle progression in dync1i2a zebrafish models and observed significantly increased apoptosis, likely due to prolonged mitosis caused by abnormal spindle morphology, and this finding offers initial insights into the cellular basis of microcephaly. Additionally, complementation studies in zebrafish demonstrate that p.Tyr247Cys attenuates gene function, consistent with protein structural analysis. Our genetic and functional data indicate that DYNC1I2 dysfunction probably causes an autosomal-recessive microcephaly syndrome and highlight further the critical roles of the dynein-1 complex in neurodevelopment.

The mechanisms regulating exhaustion of tumor-infiltrating lymphocytes (TIL) and responsiveness to PD-1 blockade remain partly unknown. In human ovarian cancer, we show that tumor-specific CD8+ TIL accumulate in tumor islets, where they engage antigen and upregulate PD-1, which restrains their functions. Intraepithelial PD-1+CD8+ TIL can be, however, polyfunctional. PD-1+ TIL indeed exhibit a continuum of exhaustion states, with variable levels of CD28 costimulation, which is provided by antigen-presenting cells (APC) in intraepithelial tumor myeloid niches. CD28 costimulation is associated with improved effector fitness of exhausted CD8+ TIL and is required for their activation upon PD-1 blockade, which also requires tumor myeloid APC. Exhausted TIL lacking proper CD28 costimulation in situ fail to respond to PD-1 blockade, and their response may be rescued by local CTLA-4 blockade and tumor APC stimulation via CD40L.

Chimeric antigen receptor (CAR) T cell therapy targeting B cell maturation antigen (BCMA) on multiple myeloma (MM) produces fast but not long-lasting responses. Reasons for treatment failure are poorly understood. CARs simultaneously targeting two antigens may represent an alternative. Here, we (1) designed and characterized novel A proliferation inducing ligand (APRIL) based dual-antigen targeting CARs, and (2) investigated mechanisms of resistance to CAR T cells with three different BCMA-binding moieties (APRIL, single-chain-variable-fragment, heavy-chain-only).

SARS-CoV-2 is the causative agent of COVID-19, the ongoing global pandemic. It has posed a worldwide challenge to human health as no effective treatment is currently available to combat the disease. Its severity has led to unprecedented collaborative initiatives for therapeutic solutions against COVID-19. Studies resorting to structure-based drug design for COVID-19 are plethoric and show good promise. Structural biology provides key insights into 3D structures, critical residues/mutations in SARS-CoV-2 proteins, implicated in infectivity, molecular recognition and susceptibility to a broad range of host species. The detailed understanding of viral proteins and their complexes with host receptors and candidate epitope/lead compounds is the key to developing a structure-guided therapeutic design. Since the discovery of SARS-CoV-2, several structures of its proteins have been determined experimentally at an unprecedented speed and deposited in the Protein Data Bank. Further, specialized structural bioinformatics tools and resources have been developed for theoretical models, data on protein dynamics from computer simulations, impact of variants/mutations and molecular therapeutics. Here, we provide an overview of ongoing efforts on developing structural bioinformatics tools and resources for COVID-19 research. We also discuss the impact of these resources and structure-based studies, to understand various aspects of SARS-CoV-2 infection and therapeutic development. These include (i) understanding differences between SARS-CoV-2 and SARS-CoV, leading to increased infectivity of SARS-CoV-2, (ii) deciphering key residues in the SARS-CoV-2 involved in receptor-antibody recognition, (iii) analysis of variants in host proteins that affect host susceptibility to infection and (iv) analyses facilitating structure-based drug and vaccine design against SARS-CoV-2.

Apart from efficacy and toxicity, many drug development failures are imputable to poor pharmacokinetics and bioavailability. Gastrointestinal absorption and brain access are two pharmacokinetic behaviors crucial to estimate at various stages of the drug discovery processes. To this end, the Brain Or IntestinaL EstimateD permeation method (BOILED-Egg) is proposed as an accurate predictive model that works by computing the lipophilicity and polarity of small molecules. Concomitant predictions for both brain and intestinal permeation are obtained from the same two physicochemical descriptors and straightforwardly translated into molecular design, owing to the speed, accuracy, conceptual simplicity and clear graphical output of the model. The BOILED-Egg can be applied in a variety of settings, from the filtering of chemical libraries at the early steps of drug discovery, to the evaluation of drug candidates for development.

Molecular docking is a computational approach for predicting the most probable position of ligands in the binding sites of macromolecules and constitutes the cornerstone of structure-based computer-aided drug design. Here, we present a new algorithm called Attracting Cavities that allows molecular docking to be performed by simple energy minimizations only. The approach consists in transiently replacing the rough potential energy hypersurface of the protein by a smooth attracting potential driving the ligands into protein cavities. The actual protein energy landscape is reintroduced in a second step to refine the ligand position. The scoring function of Attracting Cavities is based on the CHARMM force field and the FACTS solvation model. The approach was tested on the 85 experimental ligand-protein structures included in the Astex diverse set and achieved a success rate of 80% in reproducing the experimental binding mode starting from a completely randomized ligand conformer. The algorithm thus compares favorably with current state-of-the-art docking programs.

Specific metabolic pathways are activated by different nutrients to adapt the organism to available resources. Although essential, these mechanisms are incompletely defined. Here, we report that medium-chain fatty acids contained in coconut oil, a major source of dietary fat, induce the liver ω-oxidation genes Cyp4a10 and Cyp4a14 to increase the production of dicarboxylic fatty acids. Furthermore, these activate all ω- and β-oxidation pathways through peroxisome proliferator activated receptor (PPAR) α and PPARγ, an activation loop normally kept under control by dicarboxylic fatty acid degradation by the peroxisomal enzyme L-PBE. Indeed, L-pbe(-/-) mice fed coconut oil overaccumulate dicarboxylic fatty acids, which activate all fatty acid oxidation pathways and lead to liver inflammation, fibrosis, and death. Thus, the correct homeostasis of dicarboxylic fatty acids is a means to regulate the efficient utilization of ingested medium-chain fatty acids, and its deregulation exemplifies the intricate relationship between impaired metabolism and inflammation.

The mucosa-associated lymphoid tissue protein-1 (MALT1, also known as paracaspase) is a protease whose activity is essential for the activation of lymphocytes and the growth of cells derived from human diffuse large B-cell lymphomas of the activated B-cell subtype (ABC DLBCL). Crystallographic approaches have shown that MALT1 can form dimers via its protease domain, but why dimerization is relevant for the biological activity of MALT1 remains largely unknown. Using a molecular modeling approach, we predicted Glu 549 (E549) to be localized within the MALT1 dimer interface and thus potentially relevant. Experimental mutation of this residue into alanine (E549A) led to a complete impairment of MALT1 proteolytic activity. This correlated with an impaired capacity of the mutant to form dimers of the protease domain in vitro, and a reduced capacity to promote NF-κB activation and transcription of the growth-promoting cytokine interleukin-2 in antigen receptor-stimulated lymphocytes. Moreover, this mutant could not rescue the growth of ABC DLBCL cell lines upon MALT1 silencing. Interestingly, the MALT1 mutant E549A was unable to undergo monoubiquitination, which we identified previously as a critical step in MALT1 activation. Collectively, these findings suggest a model in which E549 at the dimerization interface is required for the formation of the enzymatically active, monoubiquitinated form of MALT1.

APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G) has antiretroviral activity associated with the hypermutation of viral DNA through cytosine deamination. APOBEC3G has two cytosine deaminase (CDA) domains; the catalytically inactive amino-terminal domain of APOBEC3G (N-CDA) carries the Vif interaction domain. There is no 3-D structure of APOBEC3G solved by X-ray or nuclear magnetic resonance.

To be effective as a drug, a potent molecule must reach its target in the body in sufficient concentration, and stay there in a bioactive form long enough for the expected biologic events to occur. Drug development involves assessment of absorption, distribution, metabolism and excretion (ADME) increasingly earlier in the discovery process, at a stage when considered compounds are numerous but access to the physical samples is limited. In that context, computer models constitute valid alternatives to experiments. Here, we present the new SwissADME web tool that gives free access to a pool of fast yet robust predictive models for physicochemical properties, pharmacokinetics, drug-likeness and medicinal chemistry friendliness, among which in-house proficient methods such as the BOILED-Egg, iLOGP and Bioavailability Radar. Easy efficient input and interpretation are ensured thanks to a user-friendly interface through the login-free website http://www.swissadme.ch. Specialists, but also nonexpert in cheminformatics or computational chemistry can predict rapidly key parameters for a collection of molecules to support their drug discovery endeavours.

Modulation by neuropeptides enhances several functions of acid-sensing ion channels (ASICs), such as pain sensation and acid-induced neuronal injury. The acid-induced opening of ASICs is transient, because of a rapid desensitization. Neuropeptides containing an Arg-Phe-amide motif affect ASIC desensitization and allow continuous activity of ASICs. In spite of the importance of the sustained ASIC activity during prolonged acidification, the molecular mechanisms of ASIC modulation by neuropeptides is only poorly understood. To identify the FRRFa (Phe-Arg-Arg-Phe-amide) binding site on ASIC1a, we carried out an in silico docking analysis and verified functionally the docking predictions. The docking experiments indicated three possible binding pockets, located (1) in the acidic pocket between the thumb, finger, β-ball and palm domains, (2) in a pocket at the bottom of the thumb domain, and (3) in the central vestibule along with the connected side cavities. Functional measurements of mutant ASIC1a confirmed the importance of residues of the lower palm, which encloses the central vestibule and its side cavities, for the FRRFa effects. The combined docking and functional experiments strongly suggest that FRRFa binds to the central vestibule and its side cavities to change ASIC desensitization.

The identification of patient-specific tumor antigens is complicated by the low frequency of T cells specific for each tumor antigen. Here we describe NeoScreen, a method that enables the sensitive identification of rare tumor (neo)antigens and of cognate T cell receptors (TCRs) expressed by tumor-infiltrating lymphocytes. T cells transduced with tumor antigen-specific TCRs identified by NeoScreen mediate regression of established tumors in patient-derived xenograft mice.

The haem enzyme indoleamine 2,3-dioxygenase 1 (IDO1) catalyses the rate-limiting step in the kynurenine pathway of tryptophan metabolism and plays an essential role in immunity, neuronal function, and ageing. Expression of IDO1 in cancer cells results in the suppression of an immune response, and therefore IDO1 inhibitors have been developed for use in anti-cancer immunotherapy. Here, we report an extension of our previously described highly efficient haem-binding 1,2,3-triazole and 1,2,4-triazole inhibitor series, the best compound having both enzymatic and cellular IC50 values of 34 nM. We provide enzymatic inhibition data for almost 100 new compounds and X-ray diffraction data for one compound in complex with IDO1. Structural and computational studies explain the dramatic drop in activity upon extension to pocket B, which has been observed in diverse haem-binding inhibitor scaffolds. Our data provides important insights for future IDO1 inhibitor design.

(5aR)-5a-C-pentyl-4-epi-isofagomine 1 is a powerful inhibitor of lysosomal β-galactosidase and a remarkable chaperone for mutations associated with GM1-gangliosidosis and Morquio disease type B. We report herein an improved synthesis of this compound and analogs (5a-C-methyl, pentyl, nonyl and phenylethyl derivatives), and a crystal structure of a synthetic intermediate that confirms its configuration resulting from the addition of a Grignard reagent. These compounds were evaluated as glycosidase inhibitors and their potential as chaperones for mutant lysosomal galactosidases determined. Based on these results and on docking studies, the 5-C-pentyl derivative 1 was selected as the optimal structure for further investigations: this compound induces the maturation of mutated β-galactosidase in fibroblasts of a GM1-gangliosidosis patient and promote the decrease of keratan sulfate and oligosaccharide load in patient cells. Compound 1 is clearly capable of restoring β-galactosidase activity and of promoting maturation of the protein, which should result in significant clinical benefit. These properties strongly support the development of compound 1 for the treatment of GM1-gangliosidosis and Morquio disease type B patients harboring β-galactosidase mutations sensitive to pharmacological chaperoning.

Investigate the genotype-phenotype correlations for five TGFBI (transforming growth factor, beta-induced) mutations including one novel pathogenic variant and one complex allele affecting the fourth FAS1 domain of keratoepithelin, and their potential effects on the protein's structure.