Arginine Citrullination at the C-Terminal Domain Controls RNA Polymerase II Transcription.

The post-translational modification of key residues at the C-terminal domain of RNA polymerase II (RNAP2-CTD) coordinates transcription, splicing, and RNA processing by modulating its capacity to act as a landing platform for a variety of protein complexes. Here, we identify a new modification at the CTD, the deimination of arginine and its conversion to citrulline by peptidyl arginine deiminase 2 (PADI2), an enzyme that has been associated with several diseases, including cancer. We show that, among PADI family members, only PADI2 citrullinates R1810 (Cit1810) at repeat 31 of the CTD. Depletion of PADI2 or loss of R1810 results in accumulation of RNAP2 at transcription start sites, reduced gene expression, and inhibition of cell proliferation. Cit1810 is needed for interaction with the P-TEFb (positive transcription elongation factor b) kinase complex and for its recruitment to chromatin. In this way, CTD-Cit1810 favors RNAP2 pause release and efficient transcription in breast cancer cells.

SBILib: a handle for protein modeling and engineering.

SUMMARY: The SBILib Python library provides an integrated platform for the analysis of macromolecular structures and interactions. It combines simple 3D file parsing and workup methods with more advanced analytical tools. SBILib includes modules for macromolecular interactions, loops, super-secondary structures, and biological sequences, as well as wrappers for external tools with which to integrate their results and facilitate the comparative analysis of protein structures and their complexes. The library can handle macromolecular complexes formed by proteins and/or nucleic acid molecules (i.e. DNA and RNA). It is uniquely capable of parsing and calculating protein super-secondary structure and loop geometry. We have compiled a list of example scenarios which SBILib may be applied to and provided access to these within the library. AVAILABILITY AND IMPLEMENTATION: SBILib is made available on Github at https://github.com/structuralbioinformatics/SBILib.

Novel SIX6 mutations cause recessively inherited congenital cataract, microcornea, and corneal opacification with or without coloboma and microphthalmia.

Purpose: To investigate the molecular basis of recessively inherited congenital cataract, microcornea, and corneal opacification with or without coloboma and microphthalmia in two consanguineous families. Methods: Conventional autozygosity mapping was performed using single nucleotide polymorphism (SNP) microarrays. Whole-exome sequencing was completed on genomic DNA from one affected member of each family. Exome sequence data were also used for homozygosity mapping and copy number variation analysis. PCR and Sanger sequencing were used to confirm the identification of mutations and to screen further patients. Evolutionary conservation of protein sequences was assessed using CLUSTALW, and protein structures were modeled using PyMol. Results: In family MEP68, a novel homozygous nucleotide substitution in SIX6 was found, c.547G>C, that converts the evolutionarily conserved aspartic acid residue at the 183rd amino acid in the protein to a histidine, p.(Asp183His). This residue mapped to the third helix of the DNA-binding homeobox domain in SIX6, which interacts with the major groove of double-stranded DNA. This interaction is likely to be disrupted by the mutation. In family F1332, a novel homozygous 1034 bp deletion that encompasses the first exon of SIX6 was identified, chr14:g.60975890_60976923del. Both mutations segregated with the disease phenotype as expected for a recessive condition and were absent from publicly available variant databases. Conclusions: Our findings expand the mutation spectrum in this form of inherited eye disease and confirm that homozygous human SIX6 mutations cause a developmental spectrum of ocular phenotypes that includes not only the previously described features of microphthalmia, coloboma, and congenital cataract but also corneal abnormalities.

SPServer: split-statistical potentials for the analysis of protein structures and protein-protein interactions.

BACKGROUND: Statistical potentials, also named knowledge-based potentials, are scoring functions derived from empirical data that can be used to evaluate the quality of protein folds and protein-protein interaction (PPI) structures. In previous works we decomposed the statistical potentials in different terms, named Split-Statistical Potentials, accounting for the type of amino acid pairs, their hydrophobicity, solvent accessibility and type of secondary structure. These potentials have been successfully used to identify near-native structures in protein structure prediction, rank protein docking poses, and predict PPI binding affinities. RESULTS: Here, we present the SPServer, a web server that applies the Split-Statistical Potentials to analyze protein folds and protein interfaces. SPServer provides global scores as well as residue/residue-pair profiles presented as score plots and maps. This level of detail allows users to: (1) identify potentially problematic regions on protein structures; (2) identify disrupting amino acid pairs in protein interfaces; and (3) compare and analyze the quality of tertiary and quaternary structural models. CONCLUSIONS: While there are many web servers that provide scoring functions to assess the quality of either protein folds or PPI structures, SPServer integrates both aspects in a unique easy-to-use web server. Moreover, the server permits to locally assess the quality of the structures and interfaces at a residue level and provides tools to compare the local assessment between structures. SERVER ADDRESS: https://sbi.upf.edu/spserver/ .

The rumen eukaryotome is a source of novel antimicrobial peptides with therapeutic potential.

BACKGROUND: The rise of microbial antibiotic resistance is a leading threat to the health of the human population. As such, finding new approaches to tackle these microbes, including development of novel antibiotics is vital. RESULTS: In this study, we mined a rumen eukaryotic metatranscriptomic library for novel Antimicrobial peptides (AMPs) using computational approaches and thereafter characterised the therapeutic potential of the AMPs. We identified a total of 208 potentially novel AMPs from the ruminal eukaryotome, and characterised one of those, namely Lubelisin. Lubelisin (GIVAWFWRLAR) is an α-helical peptide, 11 amino acid long with theoretical molecular weight of 1373.76 D. In the presence of Lubelisin, strains of methicillin-resistant Staphylococcus aureus (MRSA) USA300 and EMRSA-15 were killed within 30 min of exposure with ≥103 and 104 CFU/mL reduction in viable cells respectively. Cytotoxicity of Lubelisin against both human and sheep erythrocytes was low resulting in a therapeutic index of 0.43. Membrane permeabilisation assays using propidium iodide alongside transmission electron microscopy revealed that cytoplasmic membrane damage may contribute to the antimicrobial activities of Lubelisin. CONCLUSIONS: We demonstrate that the rumen eukaryotome is a viable source for the discovery of antimicrobial molecules for the treatment of bacterial infections and further development of these may provide part of the potential solution to the ongoing problem of antimicrobial resistance. The role of these AMPs in the ecological warfare within the rumen is also currently unknown.

A new genome allows the identification of genes associated with natural variation in aluminium tolerance in Brachiaria grasses.

Toxic concentrations of aluminium cations and low phosphorus availability are the main yield-limiting factors in acidic soils, which represent half of the potentially available arable land. Brachiaria grasses, which are commonly sown as forage in the tropics because of their resilience and low demand for nutrients, show greater tolerance to high concentrations of aluminium cations (Al3+) than most other grass crops. In this work, we explored the natural variation in tolerance to Al3+ between high and low tolerant Brachiaria species and characterized their transcriptional differences during stress. We identified three QTLs (quantitative trait loci) associated with root vigour during Al3+ stress in their hybrid progeny. By integrating these results with a new Brachiaria reference genome, we identified 30 genes putatively responsible for Al3+ tolerance in Brachiaria. We observed differential expression during stress of genes involved in RNA translation, response signalling, cell wall composition, and vesicle location homologous to aluminium-induced proteins involved in limiting uptake or localizing the toxin. However, there was limited regulation of malate transporters in Brachiaria, which suggests that exudation of organic acids and other external tolerance mechanisms, common in other grasses, might not be relevant in Brachiaria. The contrasting regulation of RNA translation and response signalling suggests that response timing is critical in high Al3+-tolerant Brachiaria.

Preclinical development of a humanized chimeric antigen receptor against B cell maturation antigen for multiple myeloma.

Multiple myeloma is a prevalent and incurable disease, despite the development of new and effective drugs. The recent development of chimeric antigen receptor (CAR)-T cell therapy has shown impressive results in the treatment of patients with relapsed or refractory hematological B cell malignancies. In the recent years, B-cell maturation antigen (BCMA) has appeared as a promising antigen to target using a variety of immuno-therapy treatments including CART cells, for MM patients. To this end, we generated clinical-grade murine CART cells directed against BCMA, named ARI2m cells. Having demonstrated its efficacy, and in an attempt to avoid the immune rejection of CART cells by the patient, the single chain variable fragment was humanized, creating ARI2h cells. ARI2h cells demonstrated comparable in vitro and in vivo efficacy to ARI2m cells, and superiority in cases of high tumor burden disease. In terms of inflammatory response, ARI2h cells showed a lower TNFα production and lower in vivo toxicity profile. Large-scale expansion of both ARI2m and ARI2h cells was efficiently conducted following Good Manufacturing Practice guidelines, obtaining the target CART cell dose required for treatment of multiple myeloma patients. Moreover, we demonstrate that soluble BCMA and BCMA released in vesicles impacts on CAR-BCMA activity. In summary, this study sets the bases for the implementation of a clinical trial (EudraCT code: 2019-001472-11) to study the efficacy of ARI2h cell treatment for multiple myeloma patients.

A Collection of Designed Peptides to Target SARS-CoV-2 Spike RBD-ACE2 Interaction.

The angiotensin-converting enzyme 2 (ACE2) is the receptor used by SARS-CoV and SARS-CoV-2 coronaviruses to attach to cells via the receptor-binding domain (RBD) of their viral spike protein. Since the start of the COVID-19 pandemic, several structures of protein complexes involving ACE2 and RBD as well as monoclonal antibodies and nanobodies have become available. We have leveraged the structural data to design peptides to target the interaction between the RBD of SARS-CoV-2 and ACE2 and SARS-CoV and ACE2, as contrasting exemplar, as well as the dimerization surface of ACE2 monomers. The peptides were modelled using our original method: PiPreD that uses native elements of the interaction between the targeted protein and cognate partner(s) that are subsequently included in the designed peptides. These peptides recapitulate stretches of residues present in the native interface plus novel and highly diverse conformations surrogating key interactions at the interface. To facilitate the access to this information we have created a freely available and dedicated web-based repository, PepI-Covid19 database, providing convenient access to this wealth of information to the scientific community with the view of maximizing its potential impact in the development of novel therapeutic and diagnostic agents.

Evaluation of Computationally Designed Peptides against TWEAK, a Cytokine of the Tumour Necrosis Factor Ligand Family.

The tumour necrosis factor-like weak inducer of apoptosis (TWEAK) is a member of the tumour necrosis factor ligand family and has been shown to be overexpressed in tumoral cells together with the fibroblast growth factor-inducible 14 (Fn14) receptor. TWEAK-Fn14 interaction triggers a set of intracellular pathways responsible for tumour cell invasion and migration, as well as proliferation and angiogenesis. Hence, modulation of the TWEAK-Fn14 interaction is an important therapeutic goal. The targeting of protein-protein interactions by external agents, e.g., drugs, remains a substantial challenge. Given their intrinsic features, as well as recent advances that improve their pharmacological profiles, peptides have arisen as promising agents in this regard. Here, we report, by in silico structural design validated by cell-based and in vitro assays, the discovery of four peptides able to target TWEAK. Our results show that, when added to TWEAK-dependent cellular cultures, peptides cause a down-regulation of genes that are part of TWEAK-Fn14 signalling pathway. The direct, physical interaction between the peptides and TWEAK was further elucidated in an in vitro assay which confirmed that the bioactivity shown in cell-based assays was due to the targeting of TWEAK. The results presented here are framed within early pre-clinical drug development and therefore these peptide hits represent a starting point for the development of novel therapeutic agents. Our approach exemplifies the powerful combination of in silico and experimental efforts to quickly identify peptides with desirable traits.

Methyl-CpG-binding (SmMBD2/3) and chromobox (SmCBX) proteins are required for neoblast proliferation and oviposition in the parasitic blood fluke Schistosoma mansoni.

While schistosomiasis remains a significant health problem in low to middle income countries, it also represents a recently recognised threat to more economically-developed regions. Until a vaccine is developed, this neglected infectious disease is primarily controlled by praziquantel, a drug with a currently unknown mechanism of action. By further elucidating how Schistosoma molecular components cooperate to regulate parasite developmental processes, next generation targets will be identified. Here, we continue our studies on schistosome epigenetic participants and characterise the function of a DNA methylation reader, the Schistosoma mansoni methyl-CpG-binding domain protein (SmMBD2/3). Firstly, we demonstrate that SmMBD2/3 contains amino acid features essential for 5-methyl cytosine (5mC) binding and illustrate that adult schistosome nuclear extracts (females > males) contain this activity. We subsequently show that SmMBD2/3 translocates into nuclear compartments of transfected murine NIH-3T3 fibroblasts and recombinant SmMBD2/3 exhibits 5mC binding activity. Secondly, using a yeast-two hybrid (Y2H) screen, we show that SmMBD2/3 interacts with the chromo shadow domain (CSD) of an epigenetic adaptor, S. mansoni chromobox protein (SmCBX). Moreover, fluorescent in situ hybridisation (FISH) mediated co-localisation of Smmbd2/3 and Smcbx to mesenchymal cells as well as somatic- and reproductive- stem cells confirms the Y2H results and demonstrates that these interacting partners are ubiquitously expressed and found within both differentiated as well as proliferating cells. Finally, using RNA interference, we reveal that depletion of Smmbd2/3 or Smcbx in adult females leads to significant reductions (46-58%) in the number of proliferating somatic stem cells (PSCs or neoblasts) as well as in the quantity of in vitro laid eggs. Collectively, these results further expand upon the schistosome components involved in epigenetic processes and suggest that pharmacological inhibition of SmMBD2/3 and/or SmCBX biology could prove useful in the development of future schistosomiasis control strategies.

On the mechanisms of protein interactions: predicting their affinity from unbound tertiary structures.

Motivation: The characterization of the protein-protein association mechanisms is crucial to understanding how biological processes occur. It has been previously shown that the early formation of non-specific encounters enhances the realization of the stereospecific (i.e. native) complex by reducing the dimensionality of the search process. The association rate for the formation of such complex plays a crucial role in the cell biology and depends on how the partners diffuse to be close to each other. Predicting the binding free energy of proteins provides new opportunities to modulate and control protein-protein interactions. However, existing methods require the 3D structure of the complex to predict its affinity, severely limiting their application to interactions with known structures. Results: We present a new approach that relies on the unbound protein structures and protein docking to predict protein-protein binding affinities. Through the study of the docking space (i.e. decoys), the method predicts the binding affinity of the query proteins when the actual structure of the complex itself is unknown. We tested our approach on a set of globular and soluble proteins of the newest affinity benchmark, obtaining accuracy values comparable to other state-of-art methods: a 0.4 correlation coefficient between the experimental and predicted values of ΔG and an error < 3 Kcal/mol. Availability and implementation: The binding affinity predictor is implemented and available at http://sbi.upf.edu/BADock and https://github.com/badocksbi/BADock. Contact: j.planas-iglesias@warwick.ac.uk or baldo.oliva@upf.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

MetaPred2CS: a sequence-based meta-predictor for protein-protein interactions of prokaryotic two-component system proteins.

MOTIVATION: Two-component systems (TCS) are the main signalling pathways of prokaryotes, and control a wide range of biological phenomena. Their functioning depends on interactions between TCS proteins, the specificity of which is poorly understood. RESULTS: The MetaPred2CS web-server interfaces a sequence-based meta-predictor specifically designed to predict pairing of the histidine kinase and response-regulator proteins forming TCSs. MetaPred2CS integrates six sequence-based methods using a support vector machine classifier and has been intensively tested under different benchmarking conditions: (i) species specific gene sets; (ii) neighbouring versus orphan pairs; and (iii) k-fold cross validation on experimentally validated datasets. AVAILABILITY AND IMPLEMENTATION: Web server at: http://metapred2cs.ibers.aber.ac.uk/, Source code: https://github.com/martinjvickers/MetaPred2CS or implemented as Virtual Machine at: http://metapred2cs.ibers.aber.ac.uk/download CONTACT: naf4@aber.ac.ukSupplementary information: Supplementary data are available at Bioinformatics online.

Prediction of homoprotein and heteroprotein complexes by protein docking and template-based modeling: A CASP-CAPRI experiment.

We present the results for CAPRI Round 30, the first joint CASP-CAPRI experiment, which brought together experts from the protein structure prediction and protein-protein docking communities. The Round comprised 25 targets from amongst those submitted for the CASP11 prediction experiment of 2014. The targets included mostly homodimers, a few homotetramers, and two heterodimers, and comprised protein chains that could readily be modeled using templates from the Protein Data Bank. On average 24 CAPRI groups and 7 CASP groups submitted docking predictions for each target, and 12 CAPRI groups per target participated in the CAPRI scoring experiment. In total more than 9500 models were assessed against the 3D structures of the corresponding target complexes. Results show that the prediction of homodimer assemblies by homology modeling techniques and docking calculations is quite successful for targets featuring large enough subunit interfaces to represent stable associations. Targets with ambiguous or inaccurate oligomeric state assignments, often featuring crystal contact-sized interfaces, represented a confounding factor. For those, a much poorer prediction performance was achieved, while nonetheless often providing helpful clues on the correct oligomeric state of the protein. The prediction performance was very poor for genuine tetrameric targets, where the inaccuracy of the homology-built subunit models and the smaller pair-wise interfaces severely limited the ability to derive the correct assembly mode. Our analysis also shows that docking procedures tend to perform better than standard homology modeling techniques and that highly accurate models of the protein components are not always required to identify their association modes with acceptable accuracy. Proteins 2016; 84(Suppl 1):323-348. © 2016 Wiley Periodicals, Inc.

Knowledge-based modeling of peptides at protein interfaces: PiPreD.

MOTIVATION: Protein-protein interactions (PPIs) underpin virtually all cellular processes both in health and disease. Modulating the interaction between proteins by means of small (chemical) agents is therefore a promising route for future novel therapeutic interventions. In this context, peptides are gaining momentum as emerging agents for the modulation of PPIs. RESULTS: We reported a novel computational, structure and knowledge-based approach to model orthosteric peptides to target PPIs: PiPreD. PiPreD relies on a precompiled and bespoken library of structural motifs, iMotifs, extracted from protein complexes and a fast structural modeling algorithm driven by the location of native chemical groups on the interface of the protein target named anchor residues. PiPreD comprehensive and systematically samples the entire interface deriving peptide conformations best suited for the given region on the protein interface. PiPreD complements the existing technologies and provides new solutions for the disruption of selected interactions. AVAILABILITY AND IMPLEMENTATION: Database and accessory scripts and programs are available upon request to the authors or at http://www.bioinsilico.org/PIPRED. CONTACT: narcis.fernandez@gmail.com.

InteractoMIX: a suite of computational tools to exploit interactomes in biological and clinical research.

Virtually all the biological processes that occur inside or outside cells are mediated by protein-protein interactions (PPIs). Hence, the charting and description of the PPI network, initially in organisms, the interactome, but more recently in specific tissues, is essential to fully understand cellular processes both in health and disease. The study of PPIs is also at the heart of renewed efforts in the medical and biotechnological arena in the quest of new therapeutic targets and drugs. Here, we present a mini review of 11 computational tools and resources tools developed by us to address different aspects of PPIs: from interactome level to their atomic 3D structural details. We provided details on each specific resource, aims and purpose and compare with equivalent tools in the literature. All the tools are presented in a centralized, one-stop, web site: InteractoMIX (http://interactomix.com).

ArchDB 2014: structural classification of loops in proteins.

The function of a protein is determined by its three-dimensional structure, which is formed by regular (i.e. ß-strands and α-helices) and non-periodic structural units such as loops. Compared to regular structural elements, non-periodic, non-repetitive conformational units enclose a much higher degree of variability--raising difficulties in the identification of regularities, and yet represent an important part of the structure of a protein. Indeed, loops often play a pivotal role in the function of a protein and different aspects of protein folding and dynamics. Therefore, the structural classification of protein loops is an important subject with clear applications in homology modelling, protein structure prediction, protein design (e.g. enzyme design and catalytic loops) and function prediction. ArchDB, the database presented here (freely available at http://sbi.imim.es/archdb), represents such a resource and has been an important asset for the scientific community throughout the years. In this article, we present a completely reworked and updated version of ArchDB. The new version of ArchDB features a novel, fast and user-friendly web-based interface, and a novel graph-based, computationally efficient, clustering algorithm. The current version of ArchDB classifies 149,134 loops in 5739 classes and 9608 subclasses.

Genome-wide prediction of prokaryotic two-component system networks using a sequence-based meta-predictor.

BACKGROUND: Two component systems (TCS) are signalling complexes manifested by a histidine kinase (receptor) and a response regulator (effector). They are the most abundant signalling pathways in prokaryotes and control a wide range of biological processes. The pairing of these two components is highly specific, often requiring costly and time-consuming experimental characterisation. Therefore, there is considerable interest in developing accurate prediction tools to lessen the burden of experimental work and cope with the ever-increasing amount of genomic information. RESULTS: We present a novel meta-predictor, MetaPred2CS, which is based on a support vector machine. MetaPred2CS integrates six sequence-based prediction methods: in-silico two-hybrid, mirror-tree, gene fusion, phylogenetic profiling, gene neighbourhood, and gene operon. To benchmark MetaPred2CS, we also compiled a novel high-quality training dataset of experimentally deduced TCS protein pairs for k-fold cross validation, to act as a gold standard for TCS partnership predictions. Combining individual predictions using MetaPred2CS improved performance when compared to the individual methods and in comparison with a current state-of-the-art meta-predictor. CONCLUSION: We have developed MetaPred2CS, a support vector machine-based metapredictor for prokaryotic TCS protein pairings. Central to the success of MetaPred2CS is a strategy of integrating individual predictors that improves the overall prediction accuracy, with the in-silico two-hybrid method contributing most to performance. MetaPred2CS outperformed other available systems in our benchmark tests, and is available online at http://metapred2cs.ibers.aber.ac.uk, along with our gold standard dataset of TCS interaction pairs.

Frag'r'Us: knowledge-based sampling of protein backbone conformations for de novo structure-based protein design.

MOTIVATION: The remodeling of short fragment(s) of the protein backbone to accommodate new function(s), fine-tune binding specificities or change/create novel protein interactions is a common task in structure-based computational design. Alternative backbone conformations can be generated de novo or by redeploying existing fragments extracted from protein structures i.e. knowledge-based. We present Frag'r'Us, a web server designed to sample alternative protein backbone conformations in loop regions. The method relies on a database of super secondary structural motifs called smotifs. Thus, sampling of conformations reflects structurally feasible fragments compiled from existing protein structures. Availability and implementation Frag'r'Us has been implemented as web application and is available at http://www.bioinsilico.org/FRAGRUS.

Next generation sequencing identifies mutations in Atonal homolog 7 (ATOH7) in families with global eye developmental defects.

The atonal homolog 7 (ATOH7) gene encodes a transcription factor involved in determining the fate of retinal progenitor cells and is particularly required for optic nerve and ganglion cell development. Using a combination of autozygosity mapping and next generation sequencing, we have identified homozygous mutations in this gene, p.E49V and p.P18RfsX69, in two consanguineous families diagnosed with multiple ocular developmental defects, including severe vitreoretinal dysplasia, optic nerve hypoplasia, persistent fetal vasculature, microphthalmia, congenital cataracts, microcornea, corneal opacity and nystagmus. Most of these clinical features overlap with defects in the Norrin/ß-catenin signalling pathway that is characterized by dysgenesis of the retinal and hyaloid vasculature. Our findings document Mendelian mutations within ATOH7 and imply a role for this molecule in the development of structures at the front as well as the back of the eye. This work also provides further insights into the function of ATOH7, especially its importance in retinal vascular development and hyaloid regression.

Homozygous mutations in PXDN cause congenital cataract, corneal opacity, and developmental glaucoma.

Anterior segment dysgenesis describes a group of heterogeneous developmental disorders that affect the anterior chamber of the eye and are associated with an increased risk of glaucoma. Here, we report homozygous mutations in peroxidasin (PXDN) in two consanguineous Pakistani families with congenital cataract-microcornea with mild to moderate corneal opacity and in a consanguineous Cambodian family with developmental glaucoma and severe corneal opacification. These results highlight the diverse ocular phenotypes caused by PXDN mutations, which are likely due to differences in genetic background and environmental factors. Peroxidasin is an extracellular matrix-associated protein with peroxidase catalytic activity, and we confirmed localization of the protein to the cornea and lens epithelial layers. Our findings imply that peroxidasin is essential for normal development of the anterior chamber of the eye, where it may have a structural role in supporting cornea and lens architecture as well as an enzymatic role as an antioxidant enzyme in protecting the lens, trabecular meshwork, and cornea against oxidative damage.