Computational Perspectives into Plasmepsins Structure—Function Relationship: Implications to Inhibitors Design

The development of efficient and selective antimalariais remains a challenge for the pharmaceutical industry. The aspartic proteases plasmepsins, whose inhibition leads to parasite death, are classified as targets for the design of potent drugs. Combinatorial synthesis is currently being used to generate inhibitor libraries for these enzymes, and together with computational methodologies have been demonstrated capable for the selection of lead compounds. The high structural flexibility of plasmepsins, revealed by their X-ray structures and molecular dynamics simulations, made even more complicated the prediction of putative binding modes, and therefore, the use of common computational tools, like docking and free-energy calculations. In this review, we revised the computational strategies utilized so far, for the structure-function relationship studies concerning the plasmepsin family, with special focus on the recent advances in the improvement of the linear interaction estimation (LIE) method, which is one of the most successful methodologies in the evaluation of plasmepsin-inhibitor binding affinity

Functional analysis of archaeal MBF1 by complementation studies in yeast

Contains fulltext : 95916.pdf (publisher's version ) (Open Access)BACKGROUND: Multiprotein-bridging factor 1 (MBF1) is a transcriptional co-activator that bridges a sequence-specific activator (basic-leucine zipper (bZIP) like proteins (e.g. Gcn4 in yeast) or steroid/nuclear-hormone receptor family (e.g. FTZ-F1 in insect)) and the TATA-box binding protein (TBP) in Eukaryotes. MBF1 is absent in Bacteria, but is well- conserved in Eukaryotes and Archaea and harbors a C-terminal Cro-like Helix Turn Helix (HTH) domain, which is the only highly conserved, classical HTH domain that is vertically inherited in all Eukaryotes and Archaea. The main structural difference between archaeal MBF1 (aMBF1) and eukaryotic MBF1 is the presence of a Zn ribbon motif in aMBF1. In addition MBF1 interacting activators are absent in the archaeal domain. To study the function and therefore the evolutionary conservation of MBF1 and its single domains complementation studies in yeast (mbf1Delta) as well as domain swap experiments between aMBF1 and yMbf1 were performed. RESULTS: In contrast to previous reports for eukaryotic MBF1 (i.e. Arabidopsis thaliana, insect and human) the two archaeal MBF1 orthologs, TMBF1 from the hyperthermophile Thermoproteus tenax and MMBF1 from the mesophile Methanosarcina mazei were not functional for complementation of an Saccharomyces cerevisiae mutant lacking Mbf1 (mbf1Delta). Of twelve chimeric proteins representing different combinations of the N-terminal, core domain, and the C-terminal extension from yeast and aMBF1, only the chimeric MBF1 comprising the yeast N-terminal and core domain fused to the archaeal C-terminal part was able to restore full wild-type activity of MBF1.However, as reported previously for Bombyx mori, the C-terminal part of yeast Mbf1 was shown to be not essential for function. In addition phylogenetic analyses revealed a common distribution of MBF1 in all Archaea with available genome sequence, except of two of the three Thaumarchaeota; Cenarchaeum symbiosum A and Nitrosopumilus maritimus SCM1. CONCLUSIONS: The absence of MBF1-interacting activators in the archaeal domain, the presence of a Zn ribbon motif in the divergent N-terminal domain of aMBF1 and the complementation experiments using archaeal- yeast chimeric proteins presented here suggests that archaeal MBF1 is not able to functionally interact with the transcription machinery and/or Gcn4 of S. cerevisiae. Based on modeling and structural prediction it is tempting to speculate that aMBF1 might act as a single regulator or non-essential transcription factor, which directly interacts with DNA via the positive charged linker or the basal transcription machinery via its Zn ribbon motif and the HTH domain. However, also alternative functions in ribosome biosynthesis and/or functionality have been discussed and therefore further experiments are required to unravel the function of MBF1 in Archaea

Structures of a bi-functional Kunitz-type STI family inhibitor of serine and aspartic proteases: could the aspartic protease inhibition have evolved from a canonical serine protease-binding loop?

Bi-functional inhibitors from the Kunitz-type soybean trypsin inhibitor (STI) family are glycosylated proteins able to inhibit serine and aspartic proteases. Here we report six crystal structures of the wild-type and a non-glycosylated mutant of the bifunctional inhibitor E3Ad obtained at different pH values and space groups. The crystal structures show that E3Ad adopts the typical β-trefoil fold of the STI family exhibiting some conformational changes due to pH variations and crystal packing. Despite the high sequence identity with a recently reported potato cathepsin D inhibitor (PDI), three-dimensional structures obtained in this work show a significant conformational change in the protease-binding loop proposed for aspartic protease inhibition. The E3Ad binding loop for serine protease inhibition is also proposed, based on structural similarity with a novel non-canonical conformation described for the double-headed inhibitor API-A from the Kunitz-type STI family. In addition, structural and sequence analyses suggest that bifunctional inhibitors of serine and aspartic proteases from the Kunitz-type STI family are more similar to double-headed inhibitor API-A than other inhibitors with a canonical protease-binding loop

KinMutRF: a random forest classifier of sequence variants in the human protein kinase superfamily

Background: The association between aberrant signal processing by protein kinases and human diseases such as cancer was established long time ago. However, understanding the link between sequence variants in the protein kinase superfamily and the mechanistic complex traits at the molecular level remains challenging: cells tolerate most genomic alterations and only a minor fraction disrupt molecular function sufficiently and drive disease. Results: KinMutRF is a novel random-forest method to automatically identify pathogenic variants in human kinases. Twenty six decision trees implemented as a random forest ponder a battery of features that characterize the variants: a) at the gene level, including membership to a Kinbase group and Gene Ontology terms; b) at the PFAM domain level; and c) at the residue level, the types of amino acids involved, changes in biochemical properties, functional annotations from UniProt, Phospho.ELM and FireDB. KinMutRF identifies disease-associated variants satisfactorily (Acc: 0.88, Prec:0.82, Rec:0.75, F-score:0.78, MCC:0.68) when trained and cross-validated with the 3689 human kinase variants from UniProt that have been annotated as neutral or pathogenic. All unclassified variants were excluded from the training set. Furthermore, KinMutRF is discussed with respect to two independent kinase-specific sets of mutations no included in the training and testing, Kin-Driver (643 variants) and Pon-BTK (1495 variants). Moreover, we provide predictions for the 848 protein kinase variants in UniProt that remained unclassified. A public implementation of KinMutRF, including documentation and examples, is available online (http://kinmut2.bioinfo.cnio.es). The source code for local installation is released under a GPL version 3 license, and can be downloaded from https://github.com/Rbbt-Workflows/KinMut2. Conclusions: KinMutRF is capable of classifying kinase variation with good performance. Predictions by KinMutRF compare favorably in a benchmark with other state-of-the-art methods (i.e. SIFT, Polyphen-2, MutationAssesor, MutationTaster, LRT, CADD, FATHMM, and VEST). Kinase-specific features rank as the most elucidatory in terms of information gain and are likely the improvement in prediction performance. This advocates for the development of family-specific classifiers able to exploit the discriminatory power of features unique to individual protein families

Recommendations for the classification of germline variants in the exonuclease domain of POLE and POLD1

BackgroundGermline variants affecting the proofreading activity of polymerases epsilon and delta cause a hereditary cancer and adenomatous polyposis syndrome characterized by tumors with a high mutational burden and a specific mutational spectrum. In addition to the implementation of multiple pieces of evidence for the classification of gene variants, POLE and POLD1 variant classification is particularly challenging given that non-disruptive variants affecting the proofreading activity of the corresponding polymerase are the ones associated with cancer. In response to an evident need in the field, we have developed gene-specific variant classification recommendations, based on the ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology) criteria, for the assessment of non-disruptive variants located in the sequence coding for the exonuclease domain of the polymerases.MethodsA training set of 23 variants considered pathogenic or benign was used to define the usability and strength of the ACMG/AMP criteria. Population frequencies, computational predictions, co-segregation data, phenotypic and tumor data, and functional results, among other features, were considered.ResultsGene-specific variant classification recommendations for non-disruptive variants located in the exonuclease domain of POLE and POLD1 were defined. The resulting recommendations were applied to 128 exonuclease domain variants reported in the literature and/or public databases. A total of 17 variants were classified as pathogenic or likely pathogenic, and 17 as benign or likely benign.ConclusionsOur recommendations, with room for improvement in the coming years as more information become available on carrier families, tumor molecular characteristics and functional assays, are intended to serve the clinical and scientific communities and help improve diagnostic performance, avoiding variant misclassifications

KSR induces RAS-independent MAPK pathway activation and modulates the efficacy of KRAS inhibitors.

The kinase suppressor of rat sarcoma (RAS) proteins (KSR1 and KSR2) have long been considered as scaffolding proteins required for optimal mitogen-activated protein kinase (MAPK) pathway signalling. However, recent evidence suggests that they play a more complex role within this pathway. Here, we demonstrate that ectopic expression of KSR1 or KSR2 is sufficient to activate the MAPK pathway and to induce cell proliferation in the absence of RAS proteins. In contrast, the ectopic expression of KSR proteins is not sufficient to induce cell proliferation in the absence of either rapidly accelerated fibrosarcoma (RAF) or MAPK-ERK kinase proteins, indicating that they act upstream of RAF. Indeed, KSR1 requires dimerization with at least one member of the RAF family to stimulate proliferation, an event that results in the translocation of the heterodimerized RAF protein to the cell membrane. Mutations in the conserved aspartic acid-phenylalanine-glycine motif of KSR1 that affect ATP binding impair the induction of cell proliferation. We also show that increased expression levels of KSR1 decrease the responsiveness to the KRASG12C inhibitor sotorasib in human cancer cell lines, thus suggesting that increased levels of expression of KSR may make tumour cells less dependent on KRAS oncogenic signalling.We thank M. San Roman and R. Villar for technical assistance. This work was supported by grants from the European Research Council (ERC-AG/695566, THERACAN), the Spanish Ministry of Science, Innovation and Universities (RTI2018-094664-B-I00 and RTC2017-6576-1), the Autonomous Community of Madrid (B2017/BMD-3884 iLUNG-CM), the CRIS Cancer Foundation and the Asociacion Espanola contra el Cancer (AECC) (GC166173694BARB). MB is a recipient of an Endowed Chair from the AXA Research Fund. GP has been supported by a fellowship from the Programa de Atraccion de Talento of the Autonomous Community of Madrid. SGA is a recipient of a postdoctoral fellowship from the AECC. OB is a recipient of a fellowship from the Formacion de Personal Investigador (FPI) program of the Spanish Ministry of Science, Innovation and Universities.S

Germline variation in O6-methylguanine-DNA methyltransferase (MGMT) as cause of hereditary colorectal cancer

Somatic epigenetic inactivation of the DNA repair protein O6-methylguanine DNA methyltransferase (MGMT) is frequent in colorectal cancer (CRC); however, its involvement in CRC predisposition remains unexplored. We assessed the role and relevance of MGMT germline mutations and epimutations in familial and early-onset CRC. Mutation and promoter methylation screenings were performed in 473 familial and/or early-onset mismatch repair-proficient nonpolyposis CRC cases. No constitutional MGMT inactivation by promoter methylation was observed. Of six rare heterozygous germline variants identified, c.346C > T (p.H116Y) and c.476G > A (p.R159Q), detected in three and one families respectively, affected highly conserved residues and showed segregation with cancer in available family members. In vitro, neither p.H116Y nor p.R159Q caused statistically significant reduction of MGMT repair activity. No evidence of somatic second hits was found in the studied tumors. Case-control data showed over-representation of c.346C > T (p.H116Y) in familial CRC compared to controls, but no overall association of MGMT mutations with CRC predisposition. In conclusion, germline mutations and constitutional epimutations in MGMT are not major players in hereditary CRC. Nevertheless, the over-representation of c.346C > T (p.H116Y) in our familial CRC cohort warrants further research

Germline mutations in the spindle assembly checkpoint genes BUB1 and BUB3 are infrequent in familial colorectal cancer and polyposis

Germline mutations in BUB1 and BUB3 have been reported to increase the risk of developing colorectal cancer (CRC) at young age, in presence of variegated aneuploidy and reminiscent dysmorphic traits of mosaic variegated aneuploidy syndrome. We performed a mutational analysis of BUB1 and BUB3 in 456 uncharacterized mismatch repair-proficient hereditary non-polyposis CRC families and 88 polyposis cases. Four novel or rare germline variants, one splice-site and three missense, were identified in four families. Neither variegated aneuploidy nor dysmorphic traits were observed in carriers. Evident functional effects in the heterozygous form were observed for c.1965-1G>A, but not for c.2296G>A (p.E766K), in spite of the positive co-segregation in the family. BUB1 c.2473C>T (p.P825S) and BUB3 c.77C>T (p.T26I) remained as variants of uncertain significance. As of today, the rarity of functionally relevant mutations identified in familial and/or early onset series does not support the inclusion of BUB1 and BUB3 testing in routine genetic diagnostics of familial CRC

Germline Mutations in FAN1 Cause Hereditary Colorectal Cancer by Impairing DNA Repair

Identification of genes associated with hereditary cancers facilitates management of patients with family histories of cancer. We performed exome sequencing of DNA from 3 individuals from a family with colorectal cancer who met the Amsterdam criteria for risk of hereditary nonpolyposis colorectal cancer. These individuals had mismatch repair-proficient tumors and each carried nonsense variant in the FANCD2/FANCI-associated nuclease 1 gene (FAN1), which encodes a nuclease involved in DNA inter-strand cross-link repair. We sequenced FAN1 in 176 additional families with histories of colorectal cancer and performed in vitro functional analyses of the mutant forms of FAN1 identified. We detected FAN1 mutations in approximately 3% of families who met the Amsterdam criteria and had mismatch repair-proficient cancers with no previously associated mutations. These findings link colorectal cancer predisposition to the Fanconi anemia DNA repair pathway, supporting the connection between genome integrity and cancer risk