Spatio-temporal expression of polyphenol oxidase unveils the dynamics of L-DOPA accumulation in faba bean (Vicia faba L.).

Faba bean (Vicia faba L.) is a winter season grain legume and a rich source of the anti-parkinson drug, L-3,4-dihydroxyphenylalanine (L-DOPA). The biosynthesis of L-DOPA in plants is not uniform and remains largely unexplored. While the hydroxylase activities of Tyrosine Hydroxylase (TH), the Cytochrome P450 (CYP450) class of enzymes, and Polyphenol Oxidases (PPOs) on tyrosine substrate have been reported in plants, only the roles of PPOs in L-DOPA biosynthesis have been recently established in velvet bean (Mucuna pruriens). To understand the differential accumulation of L-DOPA in different tissues of faba bean, profiling of L-Tyrosine, L-DOPA, Tyramine, and Dopamine in different tissues was performed. Differential accumulation of L-DOPA depended on tissue type and maturity. Furthermore, dopamine biosynthesis through L-DOPA from L-Tyr was confirmed in faba bean. The expression analysis of PPOs in leaf and flower tissues revealed the selective induction of only four (HePPO-2, HePPO-7, HePPO-8b, and HePPO-10) out of ten genes encoding different PPOs mined from the faba bean genome. Higher accumulation of L-DOPA in young leaves and flower buds than in mature leaves and flowers was accompanied by significantly higher expression of HePPO-10 and HePPO-7, respectively. The role of various transcription factors contributing to such metabolite dynamics was also predicted. Further exploration of this mechanism using a multi-omics approach can provide meaningful insight and pave the way for enhancing L-DOPA content in crops. Supplementary Information: The online version contains supplementary material available at 10.1007/s12298-024-01449-2.

Unveiling the Molecular Mechanism of Trastuzumab Resistance in SKBR3 and BT474 Cell Lines for HER2 Positive Breast Cancer.

HER2-positive breast cancer is one of the most prevalent forms of cancer among women worldwide. Generally, the molecular characteristics of this breast cancer include activation of human epidermal growth factor receptor-2 (HER2) and hormone receptor activation. HER2-positive is associated with a higher death rate, which led to the development of a monoclonal antibody called trastuzumab, specifically targeting HER2. The success rate of HER2-positive breast cancer treatment has been increased; however, drug resistance remains a challenge. This fact motivated us to explore the underlying molecular mechanisms of trastuzumab resistance. For this purpose, a two-fold approach was taken by considering well-known breast cancer cell lines SKBR3 and BT474. In the first fold, trastuzumab treatment doses were optimized separately for both cell lines. This was done based on the proliferation rate of cells in response to a wide variety of medication dosages. Thereafter, each cell line was cultivated with a steady dosage of herceptin for several months. During this period, six time points were selected for further in vitro analysis, ranging from the untreated cell line at the beginning to a fully resistant cell line at the end of the experiment. In the second fold, nucleic acids were extracted for further high throughput-based microarray experiments of gene and microRNA expression. Such expression data were further analyzed in order to infer the molecular mechanisms involved in the underlying development of trastuzumab resistance. In the list of differentially expressed genes and miRNAs, multiple genes (e.g., BIRC5, E2F1, TFRC, and USP1) and miRNAs (e.g., hsa miR 574 3p, hsa miR 4530, and hsa miR 197 3p) responsible for trastuzumab resistance were found. Downstream analysis showed that TFRC, E2F1, and USP1 were also targeted by hsa-miR-8485. Moreover, it indicated that miR-4701-5p was highly expressed as compared to TFRC in the SKBR3 cell line. These results unveil key genes and miRNAs as molecular regulators for trastuzumab resistance.

3D-GNOME 3.0: a three-dimensional genome modelling engine for analysing changes of promoter-enhancer contacts in the human genome.

In the current update, we added a feature for analysing changes in spatial distances between promoters and enhancers in chromatin 3D model ensembles. We updated our datasets by the novel in situ CTCF and RNAPII ChIA-PET chromatin loops obtained from the GM12878 cell line mapped to the GRCh38 genome assembly and extended the 1000 Genomes SVs dataset. To handle the new datasets, we applied GPU acceleration for the modelling engine, which gives a speed-up of 30× versus the previous versions. To improve visualisation and data analysis, we embedded the IGV tool for viewing ChIA-PET arcs with additional genes and SVs annotations. For 3D model visualisation, we added a new viewer: NGL, where we provided colouring by gene and enhancer location. The models are downloadable in mmcif and xyz format. The web server is hosted and performs calculations on DGX A100 GPU servers that provide optimal performance with multitasking. 3D-GNOME 3.0 web server provides unique insights into the topological mechanism of human variations at the population scale with high speed-up and is freely available at https://3dgnome.mini.pw.edu.pl/.

PFP-GO: Integrating protein sequence, domain and protein-protein interaction information for protein function prediction using ranked GO terms.

Protein function prediction is gradually emerging as an essential field in biological and computational studies. Though the latter has clinched a significant footprint, it has been observed that the application of computational information gathered from multiple sources has more significant influence than the one derived from a single source. Considering this fact, a methodology, PFP-GO, is proposed where heterogeneous sources like Protein Sequence, Protein Domain, and Protein-Protein Interaction Network have been processed separately for ranking each individual functional GO term. Based on this ranking, GO terms are propagated to the target proteins. While Protein sequence enriches the sequence-based information, Protein Domain and Protein-Protein Interaction Networks embed structural/functional and topological based information, respectively, during the phase of GO ranking. Performance analysis of PFP-GO is also based on Precision, Recall, and F-Score. The same was found to perform reasonably better when compared to the other existing state-of-art. PFP-GO has achieved an overall Precision, Recall, and F-Score of 0.67, 0.58, and 0.62, respectively. Furthermore, we check some of the top-ranked GO terms predicted by PFP-GO through multilayer network propagation that affect the 3D structure of the genome. The complete source code of PFP-GO is freely available at https://sites.google.com/view/pfp-go/.

Multi-scale phase separation by explosive percolation with single-chromatin loop resolution.

The 2 m-long human DNA is tightly intertwined into the cell nucleus of the size of 10 µm. The DNA packing is explained by folding of chromatin fiber. This folding leads to the formation of such hierarchical structures as: chromosomal territories, compartments; densely-packed genomic regions known as Topologically Associating Domains (TADs), or Chromatin Contact Domains (CCDs), and loops. We propose models of dynamical human genome folding into hierarchical components in human lymphoblastoid, stem cell, and fibroblast cell lines. Our models are based on explosive percolation theory. The chromosomes are modeled as graphs where CTCF chromatin loops are represented as edges. The folding trajectory is simulated by gradually introducing loops to the graph following various edge addition strategies that are based on topological network properties, chromatin loop frequencies, compartmentalization, or epigenomic features. Finally, we propose the genome folding model - a biophysical pseudo-time process guided by a single scalar order parameter. The parameter is calculated by Linear Discriminant Analysis of chromatin features. We also include dynamics of loop formation by using Loop Extrusion Model (LEM) while adding them to the system. The chromatin phase separation, where fiber folds in 3D space into topological domains and compartments, is observed when the critical number of contacts is reached. We also observe that at least 80% of the loops are needed for chromatin fiber to condense in 3D space, and this is constant through various cell lines. Overall, our in-silico model integrates the high-throughput 3D genome interaction experimental data with the novel theoretical concept of phase separation, which allows us to model event-based time dynamics of chromatin loop formation and folding trajectories.

From DNA human sequence to the chromatin higher order organisation and its biological meaning: Using biomolecular interaction networks to understand the influence of structural variation on spatial genome organisation and its functional effect.

The three-dimensional structure of the human genome has been proven to have a significant functional impact on gene expression. The high-order spatial chromatin is organised first by looping mediated by multiple protein factors, and then it is further formed into larger structures of topologically associated domains (TADs) or chromatin contact domains (CCDs), followed by A/B compartments and finally the chromosomal territories (CTs). The genetic variation observed in human population influences the multi-scale structures, posing a question regarding the functional impact of structural variants reflected by the variability of the genes expression patterns. The current methods of evaluating the functional effect include eQTLs analysis which uses statistical testing of influence of variants on spatially close genes. Rarely, non-coding DNA sequence changes are evaluated by their impact on the biomolecular interaction network (BIN) reflecting the cellular interactome that can be analysed by the classical graph-theoretic algorithms. Therefore, in the second part of the review, we introduce the concept of BIN, i.e. a meta-network model of the complete molecular interactome developed by integrating various biological networks. The BIN meta-network model includes DNA-protein binding by the plethora of protein factors as well as chromatin interactions, therefore allowing connection of genomics with the downstream biomolecular processes present in a cell. As an illustration, we scrutinise the chromatin interactions mediated by the CTCF protein detected in a ChIA-PET experiment in the human lymphoblastoid cell line GM12878. In the corresponding BIN meta-network the DNA spatial proximity is represented as a graph model, combined with the Proteins-Interaction Network (PIN) of human proteome using the Gene Association Network (GAN). Furthermore, we enriched the BIN with the signalling and metabolic pathways and Gene Ontology (GO) terms to assert its functional context. Finally, we mapped the Single Nucleotide Polymorphisms (SNPs) from the GWAS studies and identified the chromatin mutational hot-spots associated with a significant enrichment of SNPs related to autoimmune diseases. Afterwards, we mapped Structural Variants (SVs) from healthy individuals of 1000 Genomes Project and identified an interesting example of the missing protein complex associated with protein Q6GYQ0 due to a deletion on chromosome 14. Such an analysis using the meta-network BIN model is therefore helpful in evaluating the influence of genetic variation on spatial organisation of the genome and its functional effect in a cell.

Aggregated network centrality shows non-random structure of genomic and proteomic networks.

Network analysis is a powerful tool for modelling biological systems. We propose a new approach that integrates the genomic interaction data at population level with the proteomic interaction data. In our approach we use chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) data from human genome to construct a set of genomic interaction networks, considering the natural partitioning of chromatin into chromatin contact domains (CCD). The genomic networks are then mapped onto proteomic interactions, to create protein-protein interaction (PPI) subnetworks. Furthermore, the network-based topological properties of these proteomic subnetworks are investigated, namely closeness centrality, betweenness centrality and clustering coefficient. We statistically confirm, that networks identified by our method significantly differ from random networks in these network properties. Additionally, we identify one of the regions, namely chr6:32014923-33217929, as having an above-random concentration of the single nucleotide polymorphisms (SNPs) related to autoimmune diseases. Then we present it in the form of a meta-network, which includes multi-omic data: genomic contact sites (anchors), genes, proteins and SNPs. Using this example we demonstrate, that the created networks provide a valid mapping of genes to SNPs, expanding on the raw SNP dataset used.

Analysis of protein targets in pathogen-host interaction in infectious diseases: a case study on Plasmodium falciparum and Homo sapiens interaction network.

Infection and disease progression is the outcome of protein interactions between pathogen and host. Pathogen, the role player of Infection, is becoming a severe threat to life as because of its adaptability toward drugs and evolutionary dynamism in nature. Identifying protein targets by analyzing protein interactions between host and pathogen is the key point. Proteins with higher degree and possessing some topologically significant graph theoretical measures are found to be drug targets. On the other hand, exceptional nodes may be involved in infection mechanism because of some pathway process and biologically unknown factors. In this article, we attempt to investigate characteristics of host-pathogen protein interactions by presenting a comprehensive review of computational approaches applied on different infectious diseases. As an illustration, we have analyzed a case study on infectious disease malaria, with its causative agent Plasmodium falciparum acting as 'Bait' and host, Homo sapiens/human acting as 'Prey'. In this pathogen-host interaction network based on some interconnectivity and centrality properties, proteins are viewed as central, peripheral, hub and non-hub nodes and their significance on infection process. Besides, it is observed that because of sparseness of the pathogen and host interaction network, there may be some topologically unimportant but biologically significant proteins, which can also act as Bait/Prey. So, functional similarity or gene ontology mapping can help us in this case to identify these proteins.