Deep-GB: A novel deep learning model for globular protein prediction using CNN-BiLSTM architecture and enhanced PSSM with trisection strategy
Globular proteins (GPs) play vital roles in a wide range of biological processes, encompassing enzymatic catalysis and immune responses. Enzymes, among these globular proteins, facilitate biochemical reactions, while others, such as haemoglobin, contribute to essential physiological functions such as oxygen transport. Given the importance of these considerations, accurately identifying Globular proteins is essential. To address the need for precise GP identification, this research introduces an innovative approach that employs a hybrid-based deep learning model called Deep-GP. We generated two datasets based on primary sequences and developed a novel feature descriptor called, Consensus Sequence-based Trisection-Position Specific Scoring Matrix (CST-PSSM). The model training phase involved the application of deep learning techniques, including the bidirectional long short-term memory network (BiLSTM), gated recurrent unit (GRU), and convolutional neural network (CNN). The BiLSTM and CNN were hybridised for ensemble learning. The CST-PSSM-based ensemble model achieved the most accurate predictive outcomes, outperforming other competitive predictors across both training and testing datasets. This demonstrates the potential of harnessing deep learning for precise GB prediction as a robust tool to expedite research, streamline drug discovery, and unveil novel therapeutic targets.
siRNAEfficacyDB: An experimentally supported small interfering RNA efficacy database
Small interfering RNA (siRNA) has revolutionised biomedical research and drug development through precise post-transcriptional gene silencing technology. Despite its immense potential, siRNA therapy still faces technical challenges, such as delivery efficiency, targeting specificity, and molecular stability. To address these challenges and facilitate siRNA drug development, we have developed siRNAEfficacyDB, a comprehensive database that integrates experimentally validated siRNA efficacy data. This database contains 3544 siRNA records, covering 42 target genes and 5 cell lines. It provides detailed information on siRNA sequences, target genes, inhibition efficiencies, experimental techniques, cell lines, siRNA concentrations, and incubation times. siRNAEfficacyDB offers a user-friendly web interface that makes it easy to query, browse and analyse data, enabling efficient access to siRNA-related information. In summary, siRNAEfficacyDB provides a useful data foundation for siRNA drug design and optimisation, serving as a valuable resource for advancing computer-aided siRNA design research and nucleic acid drug development. siRNAEfficacyDB is freely available at https://cellknowledge.com.cn/siRNAEfficacy for non-commercial use.
A tumour-associated macrophage-based signature for deciphering prognosis and immunotherapy response in prostate cancer
For the multistage progression of prostate cancer (PCa) and resistance to immunotherapy, tumour-associated macrophage is an essential contributor. Although immunotherapy is an important and promising treatment modality for cancer, most patients with PCa are not responsive towards it. In addition to exploring new therapeutic targets, it is imperative to identify highly immunotherapy-sensitive individuals. This research aimed to establish a signature risk model, which derived from the macrophage, to assess immunotherapeutic responses and predict prognosis. Data from the UCSC-XENA, GEO and TISCH databases were extracted for analysis. Based on both single-cell datasets and bulk transcriptome profiles, a macrophage-related score (MRS) consisting of the 10-gene panel was constructed using the gene set variation analysis. MRS was highly correlated with hypoxia, angiogenesis, and epithelial-mesenchymal transition, suggesting its potential as a risk indicator. Moreover, poor immunotherapy responses and worse prognostic performance were observed in the high-MRS group of various immunotherapy cohorts. Additionally, APOE, one of the constituent genes of the MRS, affected the polarisation of macrophages. In particular, the reduced level of M2 macrophage and tumour progression suppression were observed in PCa xenografts which implanted in Apolipoprotein E-knockout mice. The constructed MRS has the potential as a robust prognostic prediction tool, and can aid in the treatment selection of PCa, especially immunotherapy options.
iGATTLDA: Integrative graph attention and transformer-based model for predicting lncRNA-Disease associations
Long non-coding RNAs (lncRNAs) have emerged as significant contributors to the regulation of various biological processes, and their dysregulation has been linked to a variety of human disorders. Accurate prediction of potential correlations between lncRNAs and diseases is crucial for advancing disease diagnostics and treatment procedures. The authors introduced a novel computational method, iGATTLDA, for the prediction of lncRNA-disease associations. The model utilised lncRNA and disease similarity matrices, with known associations represented in an adjacency matrix. A heterogeneous network was constructed, dissecting lncRNAs and diseases as nodes and their associations as edges. The Graph Attention Network (GAT) is employed to process initial features and corresponding adjacency information. GAT identified significant neighbouring nodes in the network, capturing intricate relationships between lncRNAs and diseases, and generating new feature representations. Subsequently, the transformer captures global dependencies and interactions across the entire sequence of features produced by the GAT. Consequently, iGATTLDA successfully captures complex relationships and interactions that conventional approaches may overlook. In evaluating iGATTLDA, it attained an area under the receiver operating characteristic (ROC) curve (AUC) of 0.95 and an area under the precision recall curve (AUPRC) of 0.96 with a two-layer multilayer perceptron (MLP) classifier. These results were notably higher compared to the majority of previously proposed models, further substantiating the model's efficiency in predicting potential lncRNA-disease associations by incorporating both local and global interactions. The implementation details can be obtained from https://github.com/momanyibiffon/iGATTLDA.
Comprehensive transcriptome and scRNA-seq analyses uncover the expression and underlying mechanism of SYNJ2 in papillary thyroid carcinoma
Synaptojanin 2 (SYNJ2) has crucial role in various tumors, but its role in papillary thyroid carcinoma (PTC) remains unexplored. This study first detected SYNJ2 protein expression in PTC using immunohistochemistry method and further assessed SYNJ2 mRNA expression through mRNA chip and RNA sequencing data and its association with clinical characteristics. Additionally, KEGG, GSVA, and GSEA analyses were conducted to investigate potential biological functions, while single-cell RNA sequencing data were used to explore SYNJ2's underlying mechanisms in PTC. Meanwhile, immune infiltration status in different SYNJ2 expression groups were analyzed. Besides, we investigated the immune checkpoint gene expression and implemented drug sensitivity analysis. Results indicated that SYNJ2 is highly expressed in PTC (SMD = 0.66 [95% CI: 0.17-1.15]) and could distinguish between PTC and non-PTC tissues (AUC = 0.74 [0.70-0.78]). Furthermore, the study identified 134 intersecting genes of DEGs and CEGs, mainly enriched in the angiogenesis and epithelial-mesenchymal transition (EMT) pathways. Subsequent analysis showed the above pathways were activated in PTC epithelial cells. PTC patients with high SYNJ2 expression showed higher sensitivity to the six common drugs. Summarily, SYNJ2 may promote PTC progression through angiogenesis and EMT pathways. High SYNJ2 expression is associated with better response to immunotherapy and chemotherapy.
Mechanism of action of "cistanche deserticola-Polygala" in treating Alzheimer's disease based on network pharmacology methods and molecular docking analysis
This article used network pharmacology, molecular docking, GEO analysis, and Gene Set Enrichment Analysis to obtain 38 main chemical components and 66 corresponding targets involved in Alzheimer's disease (AD) treatment in "Cistanche deserticola-Polygala". Through further Gene Ontology and Kyoto Encyclopaedia of Genes and Genomes enrichment analysis, we obtained AD signalling pathways, calcium signalling pathways, and other signalling pathways related to the treatment of AD with "Cistanche deserticola-Polygala". Molecular docking showed that most of the core chemical components had good binding ability with the core targets. This article aims to reveal the mechanism of "Cistanche deserticola-Polygala" in treating AD and provide a basis for the treatment of AD with "Cistanche deserticola-Polygala".
Developing a machine learning model with enhanced performance for predicting COVID-19 from patients presenting to the emergency room with acute respiratory symptoms
Artificial Intelligence is playing a crucial role in healthcare by enhancing decision-making and data analysis, particularly during the COVID-19 pandemic. This virus affects individuals across all age groups, but its impact is more severe on the elderly and those with underlying health issues like chronic diseases. This study aimed to develop a machine learning model to improve the prediction of COVID-19 in patients with acute respiratory symptoms. Data from 915 patients in two hospitals in Saudi Arabia were used, categorized into four groups based on chronic lung conditions and COVID-19 status. Four supervised machine learning algorithms-Random Forest, Bagging classifier, Decision Tree, and Logistic Regression-were employed to predict COVID-19. Feature selection identified 12 key variables for prediction, including CXR abnormalities, smoking status, and WBC count. The Random Forest model showed the highest accuracy at 99.07%, followed by Decision Tree, Bagging classifier, and Logistic Regression. The study concluded that machine learning algorithms, particularly Random Forest, can effectively predict and classify COVID-19 cases, supporting the development of computer-assisted diagnostic tools in healthcare.
Revealing the potential role of hub metabolism-related genes and their correlation with immune cells in acute ischemic stroke
Acute ischemic stroke (AIS) is caused by cerebral ischemia due to thrombosis in the blood vessel. The purpose of this study is to identify key genes related to metabolism to aid in the mechanism research and management of AIS.
Identification and analysis of epithelial-mesenchymal transition-related key long non-coding RNAs in hypospadias
EMT dysfunction is a dominant mechanisms of hypospadias. Thus, identification of EMT-related lncRNAs based on transcriptome sequencing data of hypospadias might provide novel molecular markers and therapeutic targets for hypospadias. First, the microarray data related to hypospadias were downloaded from Gene Expression Omnibus (GEO). Besides, the differentially expressed lncRNAs and messenger RNAs (mRNAs) related to EMT were screened to construct lncRNA-mRNA co-expression interaction pairs. In addition, the microRNA (miRNA) prediction analysis was performed through bioinformatics methods to construct a ceRNA network. Moreover, function prediction and function enrichment and pathway analyses were also performed. Finally, the core EMT-related lncRNAs were verified based on mRNA expression changes and cell functions. A total of 6 EMT-related lncRNAs were identified and 123 mRNA-lncRNA co-expression interaction pairs were screened in this study. Additionally, a ceRNA regulatory network comprising 17 mRNAs, 4 lncRNAs, and 28 miRNAs was constructed based on the prediction of hypospadias-related miRNAs. The validation results of the dataset GSE121712 revealed that only BEX1 was positively correlated with the expression of the lncRNA GNAS-AS1 (r = 0.874, P < 0.01), both of which had high expression. The cell experiment results demonstrated that interfering with the expression of GNAS-AS1 significantly promoted the proliferation, migration, and EMT of cells. Importantly, it was confirmed that GNAS-AS1 can serve as a ceRNA and play an important role in the EMT of hypospadias. Hence, it may be considered as a potential target in the treatment of this disease.
Adjusting laser power to control the heat generated by nanoparticles at the site of a patient's cells
Cancer treatment often involves heat therapy, commonly administered alongside chemotherapy and radiation therapy. The authors address the challenges posed by heat treatment methods and introduce effective control techniques. These approaches enable the precise adjustment of laser radiation over time, ensuring the tumour's core temperature attains an acceptable level with a well-defined transient response. In these control strategies, the input is the actual tumour temperature compared to the desired value, while the output governs laser radiation power. Efficient control methods are explored for regulating tumour temperature in the presence of nanoparticles and laser radiation, validated through simulations on a relevant physiological model. Initially, a Proportional-Integral-Derivative (PID) controller serves as the foundational compensator. The PID controller parameters are optimised using a combination of trial and error and the Imperialist Competitive Algorithm (ICA). ICA, known for its swift convergence and reduced computational complexity, proves instrumental in parameter determination. Furthermore, an intelligent controller based on an artificial neural network is integrated with the PID controller and compared against alternative methods. Simulation results underscore the efficacy of the combined neural network-PID controller in achieving precise temperature control. This comprehensive study illuminates promising avenues for enhancing heat therapy's effectiveness in cancer treatment.
Integrated bioinformatics analysis identified leucine rich repeat containing 15 and secreted phosphoprotein 1 as hub genes for calcific aortic valve disease and osteoarthritis
Calcific aortic valve disease (CAVD) and osteoarthritis (OA) are common diseases in the ageing population and share similar pathogenesis, especially in inflammation. This study aims to discover potential diagnostic and therapeutic targets in patients with CAVD and OA. Three CAVD datasets and one OA dataset were obtained from the Gene Expression Omnibus database. We used bioinformatics methods to search for key genes and immune infiltration, and established a ceRNA network. Immunohistochemical staining was performed to verify the expression of candidate genes in human and mice aortic valve tissues. Two key genes obtained, leucine rich repeat containing 15 (LRRC15) and secreted phosphoprotein 1 (SPP1), were further screened using machine learning and verified in human and mice aortic valve tissues. Compared to normal tissues, the infiltration of immune cells in CAVD tissues was significantly higher, and the expressions of LRRC15 and SPP1 were positively correlated with immune cells infiltration. Moreover, the ceRNA network showed extensive regulatory interactions based on LRRC15 and SPP1. The authors' findings identified LRRC15 and SPP1 as hub genes in immunological mechanisms during CAVD and OA initiation and progression, as well as potential targets for drug development.
DNA methylation plays important roles in lifestyle transition of Arthrobotrys oligospora
Trap formation is the key indicator of carnivorous lifestyle transition of nematode-trapping fungi (NTF). Here, the DNA methylation profile was explored during trap induction of Arthrobotrys oligospora, a typical NTF that captures nematodes by developing adhesive networks. Whole-genome bisulfite sequencing identified 871 methylation sites and 1979 differentially methylated regions (DMRs). This first-of-its-kind investigation unveiled the widespread presence of methylation systems in NTF, and suggested potential regulation of ribosomal RNAs through DNA methylation. Functional analysis indicated DNA methylation's involvement in complex gene regulations during trap induction, impacting multiple biological processes like response to stimulus, transporter activity, cell reproduction and molecular function regulator. These findings provide a glimpse into the important roles of DNA methylation in trap induction and offer new insights for understanding the molecular mechanisms driving carnivorous lifestyle transition of NTF.
Gene signatures of endoplasmic reticulum stress and mitophagy for prognostic risk prediction in lung adenocarcinoma
Genes associated with endoplasmic reticulum stress (ERS) and mitophagy can be conducive to predicting solid tumour prognosis. The authors aimed to develop a prognosis prediction model for these genes in lung adenocarcinoma (LUAD). Relevant gene expression and clinical information were collected from public databases including Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA). A total of 265 differentially expressed genes was finally selected (71 up-regulated and 194 downregulated) in the LUAD dataset. Among these, 15 candidate ERS and mitophagy genes (ATG12, CSNK2A1, MAP1LC3A, MAP1LC3B, MFN2, PGAM5, PINK1, RPS27A, SQSTM1, SRC, UBA52, UBB, UBC, ULK1, and VDAC1) might be critical to LUAD based on the expression analysis after crossing with the ERS and mitochondrial autophagy genes. The prediction model demonstrated the ability to effectively predict the 5-, 3-, and 1-year prognoses of LUAD patients in both GEO and TCGA databases. Moreover, high VDAC1 expression was associated with poor overall survival in LUAD (p < 0.001), suggesting it might be a critical gene for LUAD prognosis prediction. Overall, the prognosis model based on ERS and mitophagy genes in LUAD can be useful for evaluating the prognosis of patients with LUAD, and VDAC1 may serve as a promising biomarker for LUAD prognosis.
Anoikis-related genes as potential prognostic biomarkers in gastric cancer: A multilevel integrative analysis and predictive therapeutic value
Gastric cancer (GC) is a frequent malignancy of the gastrointestinal tract. Exploring the potential anoikis mechanisms and pathways might facilitate GC research.
Analysis of disulfidptosis- and cuproptosis-related LncRNAs in modulating the immune microenvironment and chemosensitivity in colon adenocarcinoma
The main objective was to establish a prognostic model utilising long non-coding RNAs associated with disulfidptosis and cuproptosis. The data for RNA-Sequence and clinicopathological information of Colon adenocarcinoma (COAD) were acquired from The Cancer Genome Atlas. A prognostic model was constructed using Cox regression and the Least Absolute Shrinkage and Selection Operator method. The model's predictive ability was assessed through principal component analysis, Kaplan-Meier analysis, nomogram etc. The ability of identifying the rates of overall survival, infiltration of immune cells, and chemosensitivity was also explored. In vitro experiments were conducted for the validation of differential expression and function of lncRNAs. A disulfidptosis and cuproptosis-related lncRNA prognostic model was constructed. The prognostic model exhibits excellent independent predictive capability for patient outcomes. Based on the authors' model, the high-risk group exhibited higher tumour mutation burdened worse survival. Besides, differences in immune cell infiltration and responsiveness to chemotherapeutic medications exist among patients with different risk scores. Furthermore, aberrant expressions in certain lncRNAs have been validated in HCT116 cells. In particular, FENDRR and SNHG7 could affect the proliferation and migration of colorectal cancer cells. Our study developed a novel prognostic signature, providing valuable insights into prognosis, immune infiltration, and chemosensitivity in COAD patients.
Excavation of gene markers associated with pancreatic ductal adenocarcinoma based on interrelationships of gene expression
Pancreatic ductal adenocarcinoma (PDAC) accounts for 95% of all pancreatic cancer cases, posing grave challenges to its diagnosis and treatment. Timely diagnosis is pivotal for improving patient survival, necessitating the discovery of precise biomarkers. An innovative approach was introduced to identify gene markers for precision PDAC detection. The core idea of our method is to discover gene pairs that display consistent opposite relative expression and differential co-expression patterns between PDAC and normal samples. Reversal gene pair analysis and differential partial correlation analysis were performed to determine reversal differential partial correlation (RDC) gene pairs. Using incremental feature selection, the authors refined the selected gene set and constructed a machine-learning model for PDAC recognition. As a result, the approach identified 10 RDC gene pairs. And the model could achieve a remarkable accuracy of 96.1% during cross-validation, surpassing gene expression-based models. The experiment on independent validation data confirmed the model's performance. Enrichment analysis revealed the involvement of these genes in essential biological processes and shed light on their potential roles in PDAC pathogenesis. Overall, the findings highlight the potential of these 10 RDC gene pairs as effective diagnostic markers for early PDAC detection, bringing hope for improving patient prognosis and survival.
Research on symptoms composition, time series evolution, and network visualisation of interstitial cystitis based on complex network community discovery algorithm
We analyzed the symptoms composition of Interstitial Cystitis (IC), the regularity of the evolution of symptoms before and after treatment, and the visualization of the community network, to provide a reference for clinical diagnosis and treatment of Interstitial Cystitis. Based on the outpatient electronic case data of 552 patients with Interstitial Cystitis, we used a complex network community discovery algorithm, directed weighted complex network, and Sankey map to mine the data of the symptoms composition of Interstitial Cystitis, the evolution of symptoms before and after treatment and the visualization of the community network, to analyze the epidemiological characteristics of interstitial cystitis symptoms in the real world. By the community division of the complex network of interstitial cystitis symptoms, We finally obtained three core symptom communities. Among them, symptom community A (bladder-related symptoms) is the symptom community with the highest proportion of nodes (60.00%) in the complex network of Interstitial Cystitis, symptom community B (non-bladder-related symptoms 1) ranks second (32.00%) in a complex network of Interstitial Cystitis, and symptom community C (non-bladder-related symptoms 2) has the lowest proportion (8.00%). There is a complex evolutionary relationship between the symptoms of Interstitial Cystitis before and after treatment. Among the single symptoms before and after treatment, the decreased rate of Day frequency is 93.22%, and the reduced urgency rate is 93.07%. The decline rate of Nocturia was 82.33%. From the perspective of different communities, the overall symptoms of symptom community A decreased by 34.39% after treatment, the general symptoms of symptom community B decreased by 35.37%, and the prevalent symptoms of symptom community C decreased by 71.43%. In the case of using diet regulation treatment to treat bladder pain, the cure rate of bladder pain can reach 22.67%. The cure rate of burning in bladders can get 15.38% with Percutaneous Sacral neuromodulation, 96.95% with diet regulation treatment, and 100% with Percutaneous Sacral neuromodulation. When using behavioral physiotherapy to treat bladder pain, 3.57% of the patient's symptoms change to bladder discomfort; 4% of the patient's symptoms change to bladder discomfort when using oral medicine to treat bladder pain.Symptom research methods based on community findings can effectively explore the rule of symptom outcome of Interstitial Cystitis before and after treatment, and the results are highly interpretable by professionals. The cover image is based on the Original Article Research on symptoms composition, time series evolution, and network visualisation of interstitial cystitis based on complex network community discovery algorithm by Lei Pang et al., https://doi.org/10.1049/syb2.12083.
Machine learning unveils RNA polymerase II binding as a predictor for SMAD2-dependent transcription dynamics in response to Actvin signalling
The transforming growth factor-β (TGF-β) superfamily, including Nodal and Activin, plays a critical role in various cellular processes. Understanding the intricate regulation and gene expression dynamics of TGF-β signalling is of interest due to its diverse biological roles. A machine learning approach is used to predict gene expression patterns induced by Activin using features, such as histone modifications, RNA polymerase II binding, SMAD2-binding, and mRNA half-life. RNA sequencing and ChIP sequencing datasets were analysed and differentially expressed SMAD2-binding genes were identified. These genes were classified into activated and repressed categories based on their expression patterns. The predictive power of different features and combinations was evaluated using logistic regression models and their performances were assessed. Results showed that RNA polymerase II binding was the most informative feature for predicting the expression patterns of SMAD2-binding genes. The authors provide insights into the interplay between transcriptional regulation and Activin signalling and offers a computational framework for predicting gene expression patterns in response to cell signalling.
The 'Other' subfamily of HECT E3 ubiquitin ligases evaluate the tumour immune microenvironment and prognosis in patients with hepatocellular carcinoma
Primary liver cancer is the sixth most common cancer and the third leading cause of cancer-related death worldwide. The role of the 'Other' subfamily of HECT E3 ligases (E3s) in hepatocellular carcinoma (HCC) remains unknown. The expression of the 'Other' HECT E3s was performed using The Cancer Genome Atlas (TCGA) data, and the authors found that the 'Other' HECT E3s were differentially expressed in HCC. Prognostic values were assessed using the Kaplan-Meier method and indicated that the high expressions of HECTD2, HECTD3, and HACE1 were associated with a worse clinical prognosis of HCC patients. The expression of HECTD2 was significantly correlated with the infiltration of CD4 T cells and neutrophils. The levels of HECTD3 and HACE1 were notably related to the dendritic cells and memory B cells infiltrated in HCC. In addition, the three previously mentioned genes have shown to be associated with immune checkpoint genes, such as FOXP3, CCR8, STAT5B, TGFB1 and TIM-3. Moreover, HECTD2 could promote the proliferative activity, cell migration and invasive ability of HCC cells. Collectively, the authors' study demonstrated that HECTD2 was a novel immune-related prognostic biomarker for HCC, providing new insight into the treatment and prognosis of HCC.
Bioinformatics identification of a T-cell-related signature for predicting prognosis and drug sensitivity in hepatocellular carcinoma
Hepatocellular carcinoma (HCC) is a fatal disease with poor clinical outcomes. T cells play a vital role in the crosstalk between the tumour microenvironment and HCC. Single-cell RNA sequencing data were downloaded from the GSE149614 dataset. The T-cell-related prognostic signature (TRPS) was developed with the integrative procedure including 10 machine learning algorithms. The TRPS was established using 7 T-cell-related markers in the Cancer Genome Atlas cohort with 1-, 2- and 3-year area under curve values of 0.820, 0.725 and 0.678, respectively. TRPS acted as an independent risk factor for HCC patients. HCC patients with a high TRPS-based risk score had a higher Tumour Immune Dysfunction and Exclusion score, lower PD1 and CTLA4 immunophenoscore and lower level of immunoactivated cells, including CD8 T cells and NK cells. The response rate was significantly higher in patients with low-risk scores in immunotherapy cohorts, including IMigor210 and GSE91061. The TRPS-based nomogram had a relatively good predictive value in evaluating the mortality risk at 1, 3 and 5 years in HCC. Overall, this study develops a TRPS by integrated bioinformatics analysis. This TRPS acted as an independent risk factor for the OS rate of HCC patients. It can screen for HCC patients who might benefit from immunotherapy, chemotherapy and targeted therapy.
Kinetic modelling of β-cell metabolism reveals control points in the insulin-regulating pyruvate cycling pathways
Insulin, a key hormone in the regulation of glucose homoeostasis, is secreted by pancreatic β-cells in response to elevated glucose levels. Insulin is released in a biphasic manner in response to glucose metabolism in β-cells. The first phase of insulin secretion is triggered by an increase in the ATP:ADP ratio; the second phase occurs in response to both a rise in ATP:ADP and other key metabolic signals, including a rise in the NADPH:NADP ratio. Experimental evidence indicates that pyruvate-cycling pathways play an important role in the elevation of the NADPH:NADP ratio in response to glucose. The authors developed a kinetic model for the tricarboxylic acid cycle and pyruvate cycling pathways. The authors successfully validated the model against experimental observations and performed a sensitivity analysis to identify key regulatory interactions in the system. The model predicts that the dicarboxylate carrier and the pyruvate transporter are the most important regulators of pyruvate cycling and NADPH production. In contrast, the analysis showed that variation in the pyruvate carboxylase flux was compensated by a response in the activity of mitochondrial isocitrate dehydrogenase (ICD ) resulting in minimal effect on overall pyruvate cycling flux. The model predictions suggest starting points for further experimental investigation, as well as potential drug targets for the treatment of type 2 diabetes.