SEMA 2.0: web-platform for B-cell conformational epitopes prediction using artificial intelligence.

Prediction of conformational B-cell epitopes is a crucial task in vaccine design and development. In this work, we have developed SEMA 2.0, a user-friendly web platform that enables the research community to tackle the B-cell epitopes prediction problem using state-of-the-art protein language models. SEMA 2.0 offers comprehensive research tools for sequence- and structure-based conformational B-cell epitopes prediction, accurate identification of N-glycosylation sites, and a distinctive module for comparing the structures of antigen B-cell epitopes enhancing our ability to analyze and understand its immunogenic properties. SEMA 2.0 website https://sema.airi.net is free and open to all users and there is no login requirement. Source code is available at https://github.com/AIRI-Institute/SEMAi.

Impact of human CD95 mutations on cell death and autoimmunity: a model.

CD95/Fas/APO-1 can trigger apoptotic as well as nonapoptotic pathways in immune cells. CD95 signaling in humans can be inhibited by several mechanisms, including mutations in the gene encoding CD95. CD95 mutations lead to autoimmune disorders, such as autoimmune lymphoproliferative syndrome (ALPS). Gaining further insight into the reported mutations of CD95 and resulting alterations of its signaling networks may provide further understanding of their presumed role in certain autoimmune diseases. For illustrative purposes and to better understand the potential outcomes of CD95 mutations, here we assign their positions to the recently determined 3D structures of human CD95. Based on this, we make certain predictions and speculate on the putative role of CD95 mutation defects in CD95-mediated signaling for certain autoimmune diseases.

PROSTATA: a framework for protein stability assessment using transformers.

MOTIVATION: Accurate prediction of change in protein stability due to point mutations is an attractive goal that remains unachieved. Despite the high interest in this area, little consideration has been given to the transformer architecture, which is dominant in many fields of machine learning. RESULTS: In this work, we introduce PROSTATA, a predictive model built in a knowledge-transfer fashion on a new curated dataset. PROSTATA demonstrates advantage over existing solutions based on neural networks. We show that the large improvement margin is due to both the architecture of the model and the quality of the new training dataset. This work opens up opportunities to develop new lightweight and accurate models for protein stability assessment. AVAILABILITY AND IMPLEMENTATION: PROSTATA is available at https://github.com/AIRI-Institute/PROSTATA and https://prostata.airi.net.

ANDDigest: a new web-based module of ANDSystem for the search of knowledge in the scientific literature.

BACKGROUND: The rapid growth of scientific literature has rendered the task of finding relevant information one of the critical problems in almost any research. Search engines, like Google Scholar, Web of Knowledge, PubMed, Scopus, and others, are highly effective in document search; however, they do not allow knowledge extraction. In contrast to the search engines, text-mining systems provide extraction of knowledge with representations in the form of semantic networks. Of particular interest are tools performing a full cycle of knowledge management and engineering, including automated retrieval, integration, and representation of knowledge in the form of semantic networks, their visualization, and analysis. STRING, Pathway Studio, MetaCore, and others are well-known examples of such products. Previously, we developed the Associative Network Discovery System (ANDSystem), which also implements such a cycle. However, the drawback of these systems is dependence on the employed ontologies describing the subject area, which limits their functionality in searching information based on user-specified queries. RESULTS: The ANDDigest system is a new web-based module of the ANDSystem tool, permitting searching within PubMed by using dictionaries from the ANDSystem tool and sets of user-defined keywords. ANDDigest allows performing the search based on complex queries simultaneously, taking into account many types of objects from the ANDSystem's ontology. The system has a user-friendly interface, providing sorting, visualization, and filtering of the found information, including mapping of mentioned objects in text, linking to external databases, sorting of data by publication date, citations number, journal H-indices, etc. The system provides data on trends for identified entities based on dynamics of interest according to the frequency of their mentions in PubMed by years. CONCLUSIONS: The main feature of ANDDigest is its functionality, serving as a specialized search for information about multiple associative relationships of objects from the ANDSystem's ontology vocabularies, taking into account user-specified keywords. The tool can be applied to the interpretation of experimental genetics data, the search for associations between molecular genetics objects, and the preparation of scientific and analytical reviews. It is presently available at https://anddigest.sysbio.ru/ .

Delineating the role of c-FLIP/NEMO interaction in the CD95 network via rational design of molecular probes.

BACKGROUND: Structural homology modeling supported by bioinformatics analysis plays a key role in uncovering new molecular interactions within gene regulatory networks. Here, we have applied this powerful approach to analyze the molecular interactions orchestrating death receptor signaling networks. In particular, we focused on the molecular mechanisms of CD95-mediated NF-κB activation and the role of c-FLIP/NEMO interaction in the induction of this pathway. RESULTS: To this end, we have created the homology model of the c-FLIP/NEMO complex using the reported structure of the v-FLIP/NEMO complex, and rationally designed peptides targeting this complex. The designed peptides were based on the NEMO structure. Strikingly, the experimental in vitro validation demonstrated that the best inhibitory effects on CD95-mediated NF-κB activation are exhibited by the NEMO-derived peptides with the substitution D242Y of NEMO. Furthermore, we have assumed that the c-FLIP/NEMO complex is recruited to the DED filaments formed upon CD95 activation and validated this assumption in silico. Further insight into the function of c-FLIP/NEMO complex was provided by the analysis of evolutionary conservation of interacting regions which demonstrated that this interaction is common in distinct mammalian species. CONCLUSIONS: Taken together, using a combination of bioinformatics and experimental approaches we obtained new insights into CD95-mediated NF-κB activation, providing manifold possibilities for targeting the death receptor network.

DNase and RNase activities of fresh cow milk lactoferrin.

Lactoferrin (LF) is an Fe3+ -binding glycoprotein first recognized in milk and then in other epithelial secretions and barrier body fluids to which many different functions have been attributed to LF, including protection from iron-induced lipid peroxidation, immunomodulation, cell growth regulation, DNA and RNA binding, as well as transcriptional activation, еtс. The polyfunctional physiological role of LF is still unclear, but it has been suggested to be responsible for primary defense against microbial and viral infections. Here, we present the first evidence that LF preparations isolated from milk of 18 cows of different breeds possess various levels of metal-dependent DNase and metal-independent RNase activities. For univocal assignment of DNase and RNase activities to cow LF, it was subjected to SDS-PAGE using gels with copolymerized calf thymus DNA or polymeric yeast RNA. In situ analysis was revealed DNase and RNase activities only in the gel zones corresponding to homogeneous LF. In contrast to human LF, cow LF possesses a relatively low cytotoxicity towards human tumor cells. The discovery that cow LF has these activities may contribute to understanding the multiple physiological functions of this extremely polyfunctional protein, including its protective role against microbial and viral infections. The computational spatial model of cow LF complex with DNA was obtained: according to the model positively charged residues of LF contact with DNA.

Organization and evolution of the chalcone synthase gene family in bread wheat and relative species.

BACKGROUND: Flavonoid compounds are secondary plant metabolites, having a functional importance in plant development, protection from pathogens and unfavorable environmental factors. Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoids; it is involved in biosynthesis of all classes of flavonoid compounds. Nevertheless, the Chs gene family in bread wheat (Triticum aestivum L.) has been not characterized yet. The aim of the current study was to investigate structural and functional organization of the Chs genes and evolution of this gene family in bread wheat and relative species. RESULTS: The nucleotide sequences of the eight Chs copies in T. aestivum were identified. Among them, two homoeologous sets of the Chs genes were located on the short (Chs-A1, -B1, -D1) and the long (Chs-A4, -B4, -D4) arms of homoeologous group 2 chromosomes. Two paralogous gene copies in the B-genome (Chs-B2, -B3) were located in the distal regions of 2BS chromosome. To clarify the origin of Chs duplications in the B-genome the phylogenetic analysis with the Chs sequences of Triticum and Aegilops species carrying ancestral genomes was conducted. It was estimated that the first duplication event occurred in the genome of the common ancestor of Triticum and Aegilops genera about 10-12 million years ago (MYA), then another copy was formed in the ancestor of the B-genome about 6-7 MYA. A homology modeling revealed high sequence similarity of bread wheat CHS enzymes. A number of short deletions in coding regions of some Chs sequences are not expected to have any significant functional effects. Estimation of transcriptional activity of the Chs copies along with a comparative analysis of their promoters structure suggested their functional specialization, which likely contributed to the maintaining of the duplicated Chs genes in wheat genome. CONCLUSIONS: From possible ten Chs copies in bread wheat genome, eight members of this family retained their intact structure and activity, while two copies appear to be lost at the level of diploid and tetraploid ancestors. Transcriptional assay along with a comparative analysis of the cis-regulatory elements revealed their functional diversification. The multiple functions supported by the Chs family are assumed to be a driving force for duplications of the Chs gene and their retention in plant genome.

Design, Synthesis and Molecular Modeling Study of Conjugates of ADP and Morpholino Nucleosides as A Novel Class of Inhibitors of PARP-1, PARP-2 and PARP-3.

We report on the design, synthesis and molecular modeling study of conjugates of adenosine diphosphate (ADP) and morpholino nucleosides as potential selective inhibitors of poly(ADP-ribose)polymerases-1, 2 and 3. Sixteen dinucleoside pyrophosphates containing natural heterocyclic bases as well as 5-haloganeted pyrimidines, and mimicking a main substrate of these enzymes, nicotinamide adenine dinucleotide (NAD+)-molecule, have been synthesized in a high yield. Morpholino nucleosides have been tethered to the ß-phosphate of ADP via a phosphoester or phosphoramide bond. Screening of the inhibiting properties of these derivatives on the autopoly(ADP-ribosyl)ation of PARP-1 and PARP-2 has shown that the effect depends upon the type of nucleobase as well as on the linkage between ADP and morpholino nucleoside. The 5-iodination of uracil and the introduction of the P-N bond in NAD+-mimetics have shown to increase inhibition properties. Structural modeling suggested that the P-N bond can stabilize the pyrophosphate group in active conformation due to the formation of an intramolecular hydrogen bond. The most active NAD+ analog against PARP-1 contained 5-iodouracil 2'-aminomethylmorpholino nucleoside with IC50 126 ± 6 µM, while in the case of PARP-2 it was adenine 2'-aminomethylmorpholino nucleoside (IC50 63 ± 10 µM). In silico analysis revealed that thymine and uracil-based NAD+ analogs were recognized as the NAD+-analog that targets the nicotinamide binding site. On the contrary, the adenine 2'-aminomethylmorpholino nucleoside-based NAD+ analogs were predicted to identify as PAR-analogs that target the acceptor binding site of PARP-2, representing a novel molecular mechanism for selective PARP inhibition. This discovery opens a new avenue for the rational design of PARP-1/2 specific inhibitors.

How human serum albumin recognizes DNA and RNA.

We show here for the first time that HSA possesses two nucleic acid-(NA) binding sites and we estimated the relative contributions of the nucleotide links of (pN)n to their total affinity for these binding sites with higher and lower affinity for NAs. The minimal ligands of these binding sites are orthophosphate (Kd=3.0 and 20.0 mm), various dNMPs (5.6-400 µm and 0.063-18 mm) and different rNMPs (4.9-30 µm and 14-250 µm). Maximal contribution to the total affinity of all NAs to the first and second sites was observed for one nucleotide and was remarkably lower for three additional nucleotide units of (pN)n (n=1-4) with a significant decrease in the contribution at n=5-6, and at n≥7-8 all dependencies reached plateaus. For d(pA)n and r(pA)n a relatively gradual decrease in the contribution to the affinity at n=1-6 was observed, while several d(pN)n, demonstrated a sharp increase in the contribution at n=2-4. Finally, all (pN)n>10 demonstrated high affinity for the first (1.4-150 nm) and the second (80-2400 nm) sites of HSA. Double-stranded NAs showed significantly lower affinity comparing with single-stranded ligands. The thermodynamic parameters characterizing the specific contribution of every nucleotide link of all (pN)1-9 (ΔG°) to their total affinity for HSA were estimated.

Antibodies against H3 and H4 histones from the sera of HIV-infected patients catalyze site-specific degradation of these histones.

Histones and their posttranslational modified forms play pivotal roles in chromatin functioning and gene transcription. Also, histones are harmful when they enter the intercellular space; their administration to animals results in systemic inflammatory and toxic responses. Autoantibodies having enzymatic activities (abzymes) are the specific feature of several autoimmune and viral diseases. Electrophoretically homogeneous IgGs containing no canonical proteases were purified from sera of HIV-infected patients by using several affinity chromatographies. In contrast to known canonical proteases, Abs from HIV-infected patients hydrolyzed exclusively only histones but no other control globular proteins. The H3 and H4 histone cleavage sites by antihistone IgGs were determined by matrix-assisted laser desorption/ionization mass spectrometry for the first time. Two clusters of H3 hydrolysis contain major (↕) and minor (*) cleavage sites: 18-K*Q*LA↕TK*A↕AR*KS↕A*P-30 and 34-G*VK*KPHR*YRPGTVA*L*R-50. H4 histone has only 1 cluster of cleavage sites containing additionally moderate (↓) cleavage sites: 15-A↕KR↕HR↕KVLR↓D*NIQ↓GIT*K-31. Sites of these histones cleavage correspond mainly to their known epitopes. It was surprising that most of the cleavage sites of histones are involved in the interaction with DNA of nucleosome core. Because histones act as damage-associated molecules, abzymes against H3 and H4 can play important role in pathogenesis of AIDs and probably other viral and immune diseases.

Molecular mechanisms underlying the impact of mutations in SOD1 on its conformational properties associated with amyotrophic lateral sclerosis as revealed with molecular modelling.

BACKGROUND: So far, little is known about the molecular mechanisms of amyotrophic lateral sclerosis onset and progression caused by SOD1 mutations. One of the hypotheses is based on SOD1 misfolding resulting from mutations and subsequent deposition of its cytotoxic aggregates. This hypothesis is complicated by the fact that known SOD1 mutations of similar clinical effect could be distributed over the whole protein structure. RESULTS: In this work, a measure of hydrogen bond stability in conformational states was studied with elastic network analysis of 35 SOD1 mutants. Twenty-eight hydrogen bonds were detected in nine of 35 mutants with their stability being significantly different from that with the wild-type. These hydrogen bonds were formed by the amino acid residues known from the literature to be located in contact between SOD1 aggregates. Additionally, residues disposed between copper binding sites of both protein subunits were found from the models to form a stiff core, which can be involved in mechanical impulse transduction between these active centres. CONCLUSIONS: The modelling highlights that both stability of the copper binding site and stability of the dimer can play an important role in ALS progression.

Correction to: selected articles from Belyaev Conference 2017: structural biology.

After publication of the article [1], it has been brought to our attention that there is a discrepancy between the publication date on the pdf and online formats. The date on the pdf is 6th February 2018 and online is 5th February 2018. The correct publication date is the one on the pdf, 6th February 2018.

Antibodies to H1 histone from the sera of HIV-infected patients recognize and catalyze site-specific degradation of this histone.

Histones and their posttranslational modifications have key roles in chromatin remodeling and gene transcription. Besides intranuclear functions, histones act as damage-associated molecules when they are released into the extracellular space. Administration of histones to animals leads to systemic inflammatory and toxic responses. Autoantibodies with enzymatic activities (abzymes) are distinctive feature of some autoimmune and viral diseases. Electrophoretically and immunologically homogeneous IgGs containing no canonical enzymes were isolated from sera of human immunodeficiency virus-infected patients by chromatography on several affinity sorbents. In contrast to canonical proteases (trypsin, chymotrypsin, and proteinase K), IgGs from human immunodeficiency virus-infected patients purified by affinity chromatography on Sepharose containing immobilized histones specifically recognized and hydrolyzed only histones but not many other tested globular proteins. Using matrix-assisted laser desorption/ionization mass spectrometry, the sites of H1 histone (193 amino acids [AAs]) cleavage by anti-H1 histone IgGs were determined for the first time. It was shown that 1 cluster of 2 major and 4 moderate sites of cleavage is located at the beginning (106-112 AAs) of the known antigenic determinants disposed at the long C-terminal sequence of H1. Two clusters of minor and very weak sites of the protein cleavage correspond to middle (8 sites, 138-158 AAs) and terminal (5 sites, 166-176 AAs) parts of the antigenic determinants. It was shown that in contrast to canonical proteases, N-terminal part of H1 histone (1-136 AAs) containing no antigenic determinants is an unpredictably very resistant against hydrolysis by abzymes, while it can be easily cleavage by canonical proteases. Because histones act as damage-associated molecules, abzymes against H1 and other histones can play important role in pathogenesis of acquired immune deficiency syndrome and probably other different diseases.

Systemic lupus erythematosus: molecular cloning and analysis of recombinant monoclonal kappa light chain NGTA2-Me-pro-ChTr possessing two different activities-trypsin-like and metalloprotease.

Polyclonal antibodies hydrolyzing myelin basic protein (MBP) can play an important role in the pathogenesis of multiple sclerosis and systemic lupus erythematosus (SLE). An immunoglobulin light chain phagemid library derived from peripheral blood lymphocytes of patients with SLE was used. The small pools of phage particles displaying light chains with different affinity for MBP were isolated by affinity chromatography on MBP-Sepharose. The fraction eluted with 0.5M NaCl was used for preparation of individual monoclonal light chains (MLChs, 26-27kDa). The clones were expressed in Escherichia coli in a soluble form; MLChs were purified by metal-chelating chromatography followed by gel filtration. In mammalians, there are serine proteases and metalloproteases. These and many other enzymes usually have only one active site and catalyze only one chemical reaction. In contrast to canonical proteases, one MLCh (NGTA2-Me-pro-ChTr) efficiently hydrolyzed MBP (but not other proteins) and four different oligopeptides corresponding to four immunodominant sequences containing cleavage sites of MBP. The proteolytic activity of MLCh was efficiently inhibited only by specific inhibitors of serine-like (phenylmethanesulfonylfluoride, PMSF) and metalloproteases (EDTA). It was shown that MLCh possess independent serine-like and metal-dependent activities. The principal existence of monoclonal antibodies with two different proteolytic activities is unexpected but very important for the further understanding of at present unknown biological functions of human antibodies.

Modulation of extrinsic apoptotic pathway by intracellular glycosylation.

The importance of post-translational modifications (PTMs), particularly O-GlcNAcylation, of cytoplasmic proteins in apoptosis has been neglected for quite a while. Modification of cytoplasmic proteins by a single N-acetylglucosamine sugar is a dynamic and reversible PTM exhibiting properties more like phosphorylation than classical O- and N-linked glycosylation. Due to the sparse information existing, we have only limited understanding of how GlcNAcylation affects cell death. Deciphering the role of GlcNAcylation in cell fate may provide further understanding of cell fate decisions. This review focus on the modulation of extrinsic apoptotic pathway via GlcNAcylation carried out by O-GlcNAc transferase (OGT) or by other bacterial effector proteins.

Arginine methylation of caspase-8 controls life/death decisions in extrinsic apoptotic networks.

Procaspase-8 is a key mediator of death receptor (DR)-mediated pathways. Recently, the role of post-translational modifications (PTMs) of procaspase-8 in controlling cell death has received increasing attention. Here, using mass spectrometry screening, pharmacological inhibition and biochemical assays, we show that procaspase-8 can be targeted by the PRMT5/RIOK1/WD45 methylosome complex. Furthermore, two potential methylation sites of PRMT5 on procaspase-8, R233 and R435, were identified in silico. R233 and R435 are highly conserved in mammals and their point mutations are among the most common mutations of caspase-8 in cancer. The introduction of mutations at these positions resulted in inhibitory effects on CD95L-induced caspase-8 activity, effector caspase activation and apoptosis. In addition, we show that procaspase-8 can undergo symmetric di-methylation. Finally, the pharmacological inhibition of PRMT5 resulted in the inhibitory effects on caspase activity and apoptotic cell death. Taken together, we have unraveled the additional control checkpoint in procaspase-8 activation and the arginine methylation network in the extrinsic apoptosis pathway.

Targeting type I DED interactions at the DED filament serves as a sensitive switch for cell fate decisions.

Activation of procaspase-8 in the death effector domain (DED) filaments of the death-inducing signaling complex (DISC) is a key step in apoptosis. In this study, a rationally designed cell-penetrating peptide, DEDid, was engineered to mimic the h2b helical region of procaspase-8-DED2 containing a highly conservative FL motif. Furthermore, mutations were introduced into the DEDid binding site of the procaspase-8 type I interface. Additionally, our data suggest that DEDid targets other type I DED interactions such as those of FADD. Both approaches of blocking type I DED interactions inhibited CD95L-induced DISC assembly, caspase activation and apoptosis. We showed that inhibition of procaspase-8 type I interactions by mutations not only diminished procaspase-8 recruitment to the DISC but also destabilized the FADD core of DED filaments. Taken together, this study offers insights to develop strategies to target DED proteins, which may be considered in diseases associated with cell death and inflammation.

Reconstruction of the regulatory hypermethylation network controlling hepatocellular carcinoma development during hepatitis C viral infection.

Hepatocellular carcinoma (HCC) has been associated with hepatitis C viral (HCV) infection as a potential risk factor. Nonetheless, the precise genetic regulatory mechanisms triggered by the virus, leading to virus-induced hepatocarcinogenesis, remain unclear. We hypothesized that HCV proteins might modulate the activity of aberrantly methylated HCC genes through regulatory pathways. Virus-host regulatory pathways, interactions between proteins, gene expression, transport, and stability regulation, were reconstructed using the ANDSystem. Gene expression regulation was statistically significant. Gene network analysis identified four out of 70 HCC marker genes whose expression regulation by viral proteins may be associated with HCC: DNA-binding protein inhibitor ID - 1 (ID1), flap endonuclease 1 (FEN1), cyclin-dependent kinase inhibitor 2A (CDKN2A), and telomerase reverse transcriptase (TERT). It suggested the following viral protein effects in HCV/human protein heterocomplexes: HCV NS3(p70) protein activates human STAT3 and NOTC1; NS2-3(p23), NS5B(p68), NS1(E2), and core(p21) activate SETD2; NS5A inhibits SMYD3; and NS3 inhibits CCN2. Interestingly, NS3 and E1(gp32) activate c-Jun when it positively regulates CDKN2A and inhibit it when it represses TERT. The discovered regulatory mechanisms might be key areas of focus for creating medications and preventative therapies to decrease the likelihood of HCC development during HCV infection.

Computer analysis of the relation between hydrogen bond stability in SOD1 mutants and the survival time of amyotrophic lateral sclerosis patients.

BACKGROUND AND OBJECTIVE: Mutations in the SOD1 protein can lead to the death of motor neurons, which, in turn, causes an incurable disease called amyotrophic lateral sclerosis (ALS). At the same time, the mechanism of the onset and development of this disease is not fully understood and is often contradictory. METHODS: Accelerated Molecular Dynamics as implemented in the OpenMM library, principal component analysis, regression analysis, random forest method. RESULTS: The stability of hydrogen bonds in 72 mutants of the SOD1 protein was calculated. Principal component analysis was carried out. Based on ten principal components acting as predictors, a multiple linear regression model was constructed. An analysis of the correlation of these ten principal components with the initial values of the stability of hydrogen bonds made it possible to characterize the contribution of known structurally and functionally important sites in the SOD1 to the scatter of ALS patients' survival time. CONCLUSION: Such an analysis made it possible to put forward hypotheses about the relationship between the stabilizing and destabilizing effects of mutations in different structurally and functionally important regions of SOD1 with the patients's survival time.

Regulation of extrinsic apoptotic signaling by c-FLIP: towards targeting cancer networks.

The extrinsic pathway is mediated by death receptors (DRs), including CD95 (APO-1/Fas) or TRAILR-1/2. Defects in apoptosis regulation lead to cancer and other malignancies. The master regulator of the DR networks is the cellular FLICE inhibitory protein (c-FLIP). In addition to its key role in apoptosis, c-FLIP may exert other cellular functions, including control of necroptosis, pyroptosis, nuclear factor κB (NF-κB) activation, and tumorigenesis. To gain further insight into the molecular mechanisms of c-FLIP action in cancer networks, we focus on the structure, isoforms, interactions, and post-translational modifications of c-FLIP. We also discuss various avenues to target c-FLIP in cancer cells for therapeutic benefit.