Primary and Secondary micro-RNA Modulation the Extrinsic Pathway of Apoptosis in Hepatocellular Carcinoma.

Abstract-One of the most common malignant liver diseases is hepatocellular carcinoma, which has a high recurrence rate and a low five-year survival rate. It is very heterogeneous both in structure and between patients, which complicates the diagnosis, prognosis and response to treatment. In this regard, an individualized, patient-centered approach becomes important, in which the use of mimetics and hsa-miRNA inhibitors involved in the pathogenesis of the disease may be determinative. From this point of view hsa-miRNAs are of interest, their aberrant expression is associated with poor prognosis for patients and is associated with tumor progression due to dysregulation of programmed cell death (apoptosis). However, the effect of hsa-miRNA on tumor development depends not only on its direct effect on expression of genes, the primary targets, but also on secondary targets mediated by regulatory pathways. While the former are actively studied, the role of secondary targets of these hsa-miRNAs in modulating apoptosis is still unclear. The present work summarizes data on hsa-miRNAs whose primary targets are key genes of the extrinsic pathway of apoptosis. Their aberrant expression is associated with early disease relapse and poor patient outcome. For these hsa-miRNAs, using the software package ANDSystem, we reconstructed the regulation of the expression of secondary targets and analyzed their impact on the activity of the extrinsic pathway of apoptosis. The potential effect of hsa-miRNAs mediated by action on secondary targets is shown to negatively correlate with the number of primary targets. It is also shown that hsa-miR-373, hsa-miR-106b and hsa-miR-96 have the highest priority as markers of hepatocellular carcinoma, whose action on secondary targets enhances their anti-apoptotic effect.

[Primary and Secondary micro-RNA Modulation the Extrinsic Pathway of Apoptosis in Hepatocellular Carcinoma].

One of the most common malignant liver diseases is hepatocellular carcinoma, which has a high recurrence rate and a low five-year survival rate. It is very heterogeneous both in structure and between patients, which complicates the diagnosis, prognosis and response to treatment. In this regard, an individualized, patient-centered approach becomes important, in which the use of mimetics and hsa-miRNA inhibitors involved in the pathogenesis of the disease may be determinative. From this point of view hsa-miRNAs are of interest, their aberrant expression is associated with poor prognosis for patients and is associated with tumor progression due to dysregulation of programmed cell death (apoptosis). However, the effect of hsa-miRNA on tumor development depends not only on its direct effect on expression of genes, the primary targets, but also on secondary targets mediated by regulatory pathways. While the former are actively studied, the role of secondary targets of these hsa-miRNAs in modulating apoptosis is still unclear. The present work summarizes data on hsa-miRNAs whose primary targets are key genes of the extrinsic pathway of apoptosis. Their aberrant expression is associated with early disease relapse and poor patient outcome. For these hsa-miRNAs, using the software package ANDSystem, we reconstructed the regulation of the expression of secondary targets and analyzed their impact on the activity of the extrinsic pathway of apoptosis. The potential effect of hsa-miRNAs mediated by action on secondary targets is shown to negatively correlate with the number of primary targets. It is also shown that hsa-miR-373, hsa-miR-106b and hsa-miR-96 have the highest priority as markers of hepatocellular carcinoma, whose action on secondary targets enhances their anti-apoptotic effect.

Prioritization of potential pharmacological targets for the development of anti-hepatocarcinoma drugs modulating the extrinsic apoptosis pathway: the reconstruction and analysis of associative gene networks help.

Hepatocellular carcinoma (HCC) is a common severe type of liver cancer characterized by an extremely aggressive course and low survival rates. It is known that disruptions in the regulation of apoptosis activation are some of the key features inherent in most cancer cells, which determines the pharmacological induction of apoptosis as an important strategy for cancer therapy. The computer design of chemical compounds capable of specifically regulating the external signaling pathway of apoptosis induction represents a promising approach for creating new effective ways of therapy for liver cancer and other oncological diseases. However, at present, most of the studies are devoted to pharmacological effects on the internal (mitochondrial) apoptosis pathway. In contrast, the external pathway induced via cell death receptors remains out of focus. Aberrant gene methylation, along with hepatitis C virus (HCV) infection, are important risk factors for the development of hepatocellular carcinoma. The reconstruction of gene networks describing the molecular mechanisms of interaction of aberrantly methylated genes with key participants of the extrinsic apoptosis pathway and their regulation by HCV proteins can provide important information when searching for pharmacological targets. In the present study, 13 criteria were proposed for prioritizing potential pharmacological targets for developing anti-hepatocarcinoma drugs modulating the extrinsic apoptosis pathway. The criteria are based on indicators of the structural and functional organization of reconstructed gene networks of hepatocarcinoma, the extrinsic apoptosis pathway, and regulatory pathways of virus-extrinsic apoptosis pathway interaction and aberrant gene methylation-extrinsic apoptosis pathway interaction using ANDSystem. The list of the top 100 gene targets ranked according to the prioritization rating was statistically significantly (p-value = 0.0002) enriched for known pharmacological targets approved by the FDA, indicating the correctness of the prioritization method. Among the promising potential pharmacological targets, six highly ranked genes (JUN, IL10, STAT3, MYC, TLR4, and KHDRBS1) are likely to deserve close attention.

Molecular-genetic pathways of hepatitis C virus regulation of the expression of cellular factors PREB and PLA2G4C, which play an important role in virus replication.

The participants of Hepatitis C virus (HCV) replication are both viral and host proteins. Therapeutic approaches based on activity inhibition of viral non-structural proteins NS3, NS5A, and NS5B are undergoing clinical trials. However, rapid mutation processes in the viral genome and acquisition of drug resistance to the existing drugs remain the main obstacles to fighting HCV. Identifying the host factors, exploring their role in HCV RNA replication, and studying viral effects on their expression is essential for understanding the mechanisms of viral replication and developing novel, effective curative approaches. It is known that the host factors PREB (prolactin regulatory element binding) and PLA2G4C (cytosolic phospholipase A2 gamma) are important for the functioning of the viral replicase complex and the formation of the platforms of HCV genome replication. The expression of PREB and PLA2G4C was significantly elevated in the presence of the HCV genome. However, the mechanisms of its regulation by HCV remain unknown. In this paper, using a text-mining technology provided by ANDSystem, we reconstructed and analyzed gene networks describing regulatory effects on the expression of PREB and PLA2G4C by HCV proteins. On the basis of the gene network analysis performed, we put forward hypotheses about the modulation of the host factors functions resulting from protein-protein interaction with HCV proteins. Among the viral proteins, NS3 showed the greatest number of regulatory linkages. We assumed that NS3 could inhibit the function of host transcription factor (TF) NOTCH1 by protein-protein interaction, leading to upregulation of PREB and PLA2G4C. Analysis of the gene networks and data on differential gene expression in HCV-infected cells allowed us to hypothesize further how HCV could regulate the expression of TFs, the binding sites of which are localized within PREB and PLA2G4C gene regions. The results obtained can be used for planning studies of the molecular-genetic mechanisms of viral-host interaction and searching for potential targets for anti-HCV therapy.

Reconstruction and analysis of the gene regulatory network for cell wall function in Arabidopsis thaliana L. leaves in response to water deficit.

The plant cell wall represents the outer compartment of the plant cell, which provides a physical barrier and triggers signaling cascades under the influence of biotic and abiotic stressors. Drought is a factor that negatively affects both plant growth and development. Cell wall proteins (CWP) play an important role in the plant response to water deficit. The adaptation mechanisms of the cell wall to water loss are of interest for identifying important genetic factors determining plant drought resistance and provide valuable information on biomarkers for further selection aimed at increasing the yield of crop plants. Using ANDSystem, a gene network describing the regulation of CWPs under water restriction conditions was reconstructed. The analysis of the gene network and the transcriptome data analysis allowed prioritizing transcription factors (TF) based on their enrichment of differentially expressed genes regulated by them. As a result, scores were calculated, acting as indicators of the association of TFs with water deficit. On the basis of the score values, eight most significant TFs were selected. The highest priority was given to the TF GBF3. CWPs were prioritized according to the criterion of summing up the scores of transcription factors regulating these genes. Among the most prioritized CWPs were the AT5G03350 gene encoding a lectin-like protein, AT4G20860 encoding BBE-like 22 required for the oxidation of cellulose degradation products, and AT4G37800 encoding xyloglucan endotransglucosylase/ hydrolase 7. Overall, the implemented algorithm could be used for prediction of regulatory interactions between transcription factors and target genes encoding cell wall proteins in plants.

Gene networks for use in metabolomic data analysis of blood plasma from patients with postoperative delirium.

Postoperative delirium (POD) is considered one of the most severe complications, resulting in impaired cognitive function, extended hospitalization, and higher treatment costs. The challenge of early POD diagnosis becomes particularly significant in cardiac surgery cases, as the incidence of this complication exceeds 50 % in certain patient categories. While it is known that neuroinflammation, neurotransmitter imbalances, disruptions in neuroendocrine regulation, and interneuronal connections contribute significantly to the development of POD, the molecular, genetic mechanisms of POD in cardiac surgery patients, along with potential metabolomic diagnostic markers, remain inadequately understood. In this study, blood plasma was collected from a group of patients over 65 years old after cardiac surgery involving artificial circulation. The collected samples were analyzed for sphingomyelin content and quantity using high-performance liquid chromatography coupled with mass spectrometry (HPLC-MS/MS) methods. The analysis revealed four significantly different sphingomyelin contents in patients with POD compared to those who did not develop POD (control group). Employing gene network reconstruction, we perceived a set of 82 regulatory enzymes affiliated with the genetic coordination of the sphingolipid metabolism pathway. Within this set, 47 are assumed to be regulators of gene expression, governing the transcription of enzymes pivotal to the metabolic cascade. Complementing this, an additional assembly of 35 regulators are considered to be regulators of activity, degradation, and translocation dynamics of enzymes integral to the aforementioned pathway. Analysis of the overrepresentation of diseases with which these regulatory proteins are associated showed that the regulators can be categorized into two groups, associated with cardiovascular pathologies (CVP) and neuropsychiatric diseases (NPD), respectively. The regulators associated with CVP are expectedly related to the effects on myocardial tissue during surgery. It is hypothesized that dysfunction of NPD-associated regulators may specifically account for the development of POD after cardiac surgery. Thus, the identified regulatory genes may provide a basis for planning further experiments, in order to study disorders at the level of expression of these genes, as well as impaired function of proteins encoded by them in patients with POD. The identified significant sphingolipids can be considered as potential markers of POD.

Plasma metabolomics and gene regulatory networks analysis reveal the role of nonstructural SARS-CoV-2 viral proteins in metabolic dysregulation in COVID-19 patients.

Metabolomic analysis of blood plasma samples from COVID-19 patients is a promising approach allowing for the evaluation of disease progression. We performed the metabolomic analysis of plasma samples of 30 COVID-19 patients and the 19 controls using the high-performance liquid chromatography (HPLC) coupled with tandem mass spectrometric detection (LC-MS/MS). In our analysis, we identified 103 metabolites enriched in KEGG metabolic pathways such as amino acid metabolism and the biosynthesis of aminoacyl-tRNAs, which differed significantly between the COVID-19 patients and the controls. Using ANDSystem software, we performed the reconstruction of gene networks describing the potential genetic regulation of metabolic pathways perturbed in COVID-19 patients by SARS-CoV-2 proteins. The nonstructural proteins of SARS-CoV-2 (orf8 and nsp5) and structural protein E were involved in the greater number of regulatory pathways. The reconstructed gene networks suggest the hypotheses on the molecular mechanisms of virus-host interactions in COVID-19 pathology and provide a basis for the further experimental and computer studies of the regulation of metabolic pathways by SARS-CoV-2 proteins. Our metabolomic analysis suggests the need for nonstructural protein-based vaccines and the control strategy to reduce the disease progression of COVID-19.

Computer analysis of regulation of hepatocarcinoma marker genes hypermethylated by HCV proteins.

Hepatitis C virus (HCV) is a risk factor that leads to hepatocellular carcinoma (HCC) development. Epigenetic changes are known to play an important role in the molecular genetic mechanisms of virus-induced oncogenesis. Aberrant DNA methylation is a mediator of epigenetic changes that are closely associated with the HCC pathogenesis and considered a biomarker for its early diagnosis. The ANDSystem software package was used to reconstruct and evaluate the statistical significance of the pathways HCV could potentially use to regulate 32 hypermethylated genes in HCC, including both oncosuppressor and protumorigenic ones identified by genome-wide analysis of DNA methylation. The reconstructed pathways included those affecting protein-protein interactions (PPI), gene expression, protein activity, stability, and transport regulations, the expression regulation pathways being statistically significant. It has been shown that 8 out of 10 HCV proteins were involved in these pathways, the HCV NS3 protein being implicated in the largest number of regulatory pathways. NS3 was associated with the regulation of 5 tumor-suppressor genes, which may be the evidence of its central role in HCC pathogenesis. Analysis of the reconstructed pathways has demonstrated that following the transcription factor inhibition caused by binding to viral proteins, the expression of a number of oncosuppressors (WT1, MGMT, SOCS1, P53) was suppressed, while the expression of others (RASF1, RUNX3, WIF1, DAPK1) was activated. Thus, the performed gene-network reconstruction has shown that HCV proteins can influence not only the methylation status of oncosuppressor genes, but also their transcriptional regulation. The results obtained can be used in the search for pharmacological targets to develop new drugs against HCV-induced HCC.

Prioritization of biological processes based on the reconstruction and analysis of associative gene networks describing the response of plants to adverse environmental factors.

Methods for prioritizing or ranking candidate genes according to their importance based on specif ic criteria via the analysis of gene networks are widely used in biomedicine to search for genes associated with diseases and to predict biomarkers, pharmacological targets and other clinically relevant molecules. These methods have also been used in other f ields, particularly in crop production. This is largely due to the development of technologies to solve problems in marker-oriented and genomic selection, which requires knowledge of the molecular genetic mechanisms underlying the formation of agriculturally valuable traits. A new direction for the study of molecular genetic mechanisms is the prioritization of biological processes based on the analysis of associative gene networks. Associative gene networks are heterogeneous networks whose vertices can depict both molecular genetic objects (genes, proteins, metabolites, etc.) and the higher-level factors (biological processes, diseases, external environmental factors, etc.) related to regulatory, physicochemical or associative interactions. Using a previously developed method, biological processes involved in plant responses to increased cadmium content, saline stress and drought conditions were prioritized according to their degree of connection with the gene networks in the SOLANUM TUBEROSUM knowledge base. The prioritization results indicate that fundamental processes, such as gene expression, post-translational modif ications, protein degradation, programmed cell death, photosynthesis, signal transmission and stress response play important roles in the common molecular genetic mechanisms for plant response to various adverse factors. On the other hand, a group of processes related to the development of seeds ("seeding development") was revealed to be drought specif ic, while processes associated with ion transport ("ion transport") were included in the list of responses specif ic to salt stress and processes associated with the metabolism of lipids were found to be involved specif ically in the response to cadmium.

[Prediction of Bacterial and Archaeal Allergenicity with AllPred Program].

Nowadays, allergic disorders have become one of the most important social problems in the world. This can be related to the advent of new allergenic agents in the environment, as well as an increasing density of human contact with known allergens, including various proteins. Thus, the development of computer programs designed for the prediction of allergenic properties of proteins becomes one of the urgent tasks of mo dern bioinformatics. Previously we developed a web accessible Allpred Program (http://www-bionet.sscc.ru/ psd/cgi-bin/programs/Allpred/allpred.cgi) that allows users to assess the allergenicity of proteins by taking into account the characteristics of their spatial structure. In this paper, using AllPred, we predicted the allergenicity of proteins from 462 archaea and bacteria species for which a complete genome was available. The segregation of considered proteins on archaea and bacteria has shown that allergens are predicted more often among archaea than among bacteria. The division of these proteins into groups according to their intracellular localization has revealed that the majority of allergenic proteins were among the secreted proteins. The application of methods for predicting the level of gene expression of microorganisms based on DNA sequence analysis showed a statistically significant relationship between the expression level of the proteins and their allergenicity. This analysis has revealed that potentially allergenic proteins were more common among highly expressed proteins. Sorting microorganisms into the pathogenic and nonpathogenic groups has shown that pathogens can potentially be more allergenic because of a statistically significant greater number of allergens predicted among their proteins.

The substitutions G245C and G245D in the Zn(2+)-binding pocket of the p53 protein result in differences of conformational flexibility of the DNA-binding domain.

Transcription activation of the proapoptotic target genes is a means by which the p53 protein implements its function of tumor suppression. Zn(2+) is a known regulator of p53 binding to the target genes. We have previously obtained an evidence that amino acid substitutions in the p53 Zn(2+)-binding pocket can presumably exert an influence on Zn(2+) position in the Zn(2+)-p53 complex and thereby affect p53 binding to DNA. With these background considerations, our aim was to estimate the effect of the putative changes in the Zn(2+) position in its binding pocket due to the G245C and G245D substitutions on the conformation of the p53 DNA-binding motif. Statistical analysis of the molecular dynamics (MD) trajectories of the mutant p53-Zn(2+) complexes was used to detect significant deviations in conformation of the mutant p53 forms. MD simulations demonstrated that (1) the two substitutions in the Zn(2+)-binding pocket caused changes in the conformation of the p53 DNA-binding motif, as compared with the wild-type (WT) p53; (2) binding of Zn(2+) to the p53 mutant forms reduced the effect of the substitutions on conformational change; and (3) Zn(2+) binding in the normal position compensated the effect of the mutations on the conformation in comparison to the altered Zn(2+) position.

[Protein Structure Discovery: software package to perform computational proteomics tasks].

Software-information system Protein Structure Discovery was developed. The system can be used for the wide range of tasks in the field of computer proteomics including prediction of function, structure and immunological properties of proteins. A specially created section of the system allows evaluating the quantitative and qualitative effects of mutations on the structural and functional properties of proteins. There are 19 of different programs integrated into the system, including the database of protein functional sites PDBSite, a PDBSiteScan program for the prediction of functional sites in three-dimensional structures of proteins, and WebProAnalyst program for the quantitative analysis of the structure-activity relationship of proteins. Protein Structure Discovery program has a Web interface and is available for users through the Internet (http://www-bionet.sscc.ru/psd/). For example, binding sites of zinc ion and ADP showed high stability of the method to errors PDBSiteScan reconstruction of spatial structures of proteins in the recognition of functional sites in model structures.

ANDVisio: a new tool for graphic visualization and analysis of literature mined associative gene networks in the ANDSystem.

The ANDVisio tool is designed to reconstruct and analyze associative gene networks in the earlier developed Associative Network Discovery System (ANDSystem) software package. The ANDSystem incorporates utilities for automated extraction of knowledge from Pubmed published scientific texts, analysis of factographic databases, also the ANDCell database containing information on molecular-genetic events retrieved from texts and databases. ANDVisio is a new user's interface to the ANDCell database stored in a remote server. ANDVisio provides graphic visualization, editing, search, also saving of associative gene networks in different formats resulting from user's request. The associative gene networks describe semantic relationships between molecular-genetic objects (proteins, genes, metabolites and others), biological processes, and diseases. ANDVisio is provided with various tools to support filtering by object types, relationships between objects and information sources; graph layout; search of the shortest pathway; cycles in graphs.