t-Distributed Stochastic Neighbor Embedding Method with the Least Information Loss for Macromolecular Simulations

TRIM44, a Novel Prognostic Marker, Supports the Survival of Proteasome-Resistant Multiple Myeloma Cells

TRIM44, a tripartite motif (TRIM) family member, is pivotal in linking the ubiquitin-proteasome system (UPS) to autophagy in multiple myeloma (MM). However, its prognostic impact and therapeutic potential remain underexplored. Here, we report that TRIM44 overexpression is associated with poor prognosis in a Multiple Myeloma Research Foundation (MMRF) cohort of 858 patients, persisting across primary and recurrent MM cases. TRIM44 expression notably increases in advanced MM stages, indicating its potential role in disease progression. Single-cell RNA sequencing across MM stages showed significant TRIM44 upregulation in smoldering MM (SMM) and MM compared to normal bone marrow, especially in patients with t(4;14) cytogenetic abnormalities. This analysis further identified high TRIM44 expression as predictive of lower responsiveness to proteasome inhibitor (PI) treatments, underscoring its critical function in the unfolded protein response (UPR) in TRIM44-high MM cells. Our findings also demonstrate that TRIM44 facilitates SQSTM1 oligomerization under oxidative stress, essential for its phosphorylation and subsequent autophagic degradation. This process supports the survival of PI-resistant MM cells by activating the NRF2 pathway, which is crucial for oxidative stress response and, potentially, other chemotherapy-induced stressors. Additionally, TRIM44 counters the TRIM21-mediated suppression of the antioxidant response, enhancing MM cell survival under oxidative stress. Collectively, our discoveries highlight TRIM44's significant role in MM progression and resistance to therapy, suggesting its potential value as a therapeutic target.

An isotopically-labelled temporal mass spectrometry imaging data analysis workflow to reveal glucose spatial metabolism patterns in bovine lens tissue

Application of stable isotopically labelled (SIL) molecules in Matrix-Assisted Laser Desorption/Ionization Mass Spectrometry Imaging (MALDI-MSI) over a series of time points allows the temporal and spatial dynamics of biochemical reactions to be tracked in a biological system. However, these large kinetic MSI datasets and the inherent variability of biological replicates presents significant challenges to the rapid analysis of the data. In addition, manual annotation of downstream SIL metabolites involves human input to carefully analyse the data based on prior knowledge and personal expertise. To overcome these challenges to the analysis of spatiotemporal MALDI-MSI data and improve the efficiency of SIL metabolite identification, a bioinformatics pipeline has been developed and demonstrated by analysing normal bovine lens glucose metabolism as a model system. The pipeline consists of spatial alignment to mitigate the impact of sample variability and ensure spatial comparability of the temporal data, dimensionality reduction to rapidly map regional metabolic distinctions within the tissue, and metabolite annotation coupled with pathway enrichment modules to summarise and display the metabolic pathways induced by the treatment. This pipeline will be valuable for the spatial metabolomics community to analyse kinetic MALDI-MSI datasets, enabling rapid characterisation of spatio-temporal metabolic patterns from tissues of interest.

Artificial intelligence in endodontics: Data preparation, clinical applications, ethical considerations, limitations, and future directions

Artificial intelligence (AI) is emerging as a transformative technology in healthcare, including endodontics. A gap in knowledge exists in understanding AI's applications and limitations among endodontic experts. This comprehensive review aims to (A) elaborate on technical and ethical aspects of using data to implement AI models in endodontics; (B) elaborate on evaluation metrics; (C) review the current applications of AI in endodontics; and (D) review the limitations and barriers to real-world implementation of AI in the field of endodontics and its future potentials/directions. The article shows that AI techniques have been applied in endodontics for critical tasks such as detection of radiolucent lesions, analysis of root canal morphology, prediction of treatment outcome and post-operative pain and more. Deep learning models like convolutional neural networks demonstrate high accuracy in these applications. However, challenges remain regarding model interpretability, generalizability, and adoption into clinical practice. When thoughtfully implemented, AI has great potential to aid with diagnostics, treatment planning, clinical interventions, and education in the field of endodontics. However, concerted efforts are still needed to address limitations and to facilitate integration into clinical workflows.

CENsible: Interpretable Insights into Small-Molecule Binding with Context Explanation Networks

We present a novel and interpretable approach for assessing small-molecule binding using context explanation networks. Given the specific structure of a protein/ligand complex, our CENsible scoring function uses a deep convolutional neural network to predict the contributions of precalculated terms to the overall binding affinity. We show that CENsible can effectively distinguish active vs inactive compounds for many systems. Its primary benefit over related machine-learning scoring functions, however, is that it retains interpretability, allowing researchers to identify the contribution of each precalculated term to the final affinity prediction, with implications for subsequent lead optimization.

Melting entropy of crystals determined by electron-beam-induced configurational disordering

Upon melting, the molecules in a crystal explore numerous configurations, reflecting an increase in disorder. The molar entropy of disorder can be defined by Boltzmann's formula ΔSd = Rln(Wd), where Wd is the increase in the number of microscopic states, so far inaccessible experimentally. We found that the Arrhenius frequency factor A of the electron diffraction signal decay provides Wd through an experimental equation A = AINTWd, where AINT is an inelastic scattering cross section. The method connects Clausius and Boltzmann experimentally and supplements the Clausius approach, being applicable to a femtogram quantity of thermally unstable and biomolecular crystals. The data also showed that crystal disordering and crystallization of melt are reciprocal, both governed by the entropy change but manifesting in opposite directions.

Machine Learning Deciphered Molecular Mechanistics with Accurate Kinetic and Thermodynamic Prediction

Time-lagged independent component analysis (tICA) and the Markov state model (MSM) have been extensively employed for extracting conformational dynamics and kinetic community networks from unbiased trajectory ensembles. However, these techniques may not be the optimal choice for elucidating transition mechanisms within low-dimensional representations, especially for intricate biosystems. Unraveling the association mechanism in such complex systems always necessitates permutations of several essential independent components or collective variables, a process that is inherently obscure and may require empirical knowledge for selection. To address these challenges, we have implemented an integrated unsupervised dimension reduction model: uniform manifold approximation and projection (UMAP) with hierarchy density-based spatial clustering of applications with noise (HDBSCAN). This approach effectively generates low-dimensional configurational embeddings. The hierarchical application of this architecture, in conjunction with MSM, reveals global kinetic connectivity while identifying local conformational states. Consequently, our methodology establishes a multiscale mechanistic elucidation framework. Leveraging the benefits of the uniform sample distribution and a denoising approach, our model demonstrates robustness in preserving global and local data structures compared to traditional dimension reduction methods in the field of MD analysis area. The interpretability of hyperparameter selection and compatibility with downstream tasks are cross-validated across various simulation data sets, utilizing both computational evaluation metrics and experimental kinetic observables. Furthermore, the predicted Mcl1-BH3 association kinetics (0.76 s-1) is in close agreement with surface plasmon resonance experiments (0.12 s-1), affirming the plausibility of the identified pathway composed of representative conformations. We anticipate that the devised workflow will serve as a foundational framework for studying recognition patterns in complex biological systems. Its contributions extend to the exploration of protein functional dynamics and rational drug design, offering a potent avenue for advancing research in these domains.

Dynamic changes of immunocyte subpopulations in thermogenic activation of adipose tissues

Objectives

The effects of cold exposure on whole-body metabolism in humans have gained increasing attention. Brown or beige adipose tissues are crucial in cold-induced thermogenesis to dissipate energy and thus have the potential to combat metabolic disorders. Despite the immune regulation of thermogenic adipose tissues, the overall changes in vital immune cells during distinct cold periods remain elusive. This study aimed to discuss the overall changes in immune cells under different cold exposure periods and to screen several potential immune cell subpopulations on thermogenic regulation.

Methods

Cibersort and mMCP-counter algorithms were employed to analyze immune infiltration in two (brown and beige) thermogenic adipose tissues under distinct cold periods. Changes in some crucial immune cell populations were validated by reanalyzing the single-cell sequencing dataset (GSE207706). Flow cytometry, immunofluorescence, and quantitative real-time PCR assays were performed to detect the proportion or expression changes in mouse immune cells of thermogenic adipose tissues under cold challenge.

Results

The proportion of monocytes, naïve, and memory T cells increased, while the proportion of NK cells decreased under cold exposure in brown adipose tissues.

Conclusion

Our study revealed dynamic changes in immune cell profiles in thermogenic adipose tissues and identified several novel immune cell subpopulations, which may contribute to thermogenic activation of adipose tissues under cold exposure.

Research progress of DDR1 inhibitors in the treatment of multiple human diseases

Discoidin domain receptor 1 (DDR1) is a collagen-activated receptor tyrosine kinase (RTK) and plays pivotal roles in regulating cellular functions such as proliferation, differentiation, invasion, migration, and matrix remodeling. DDR1 is involved in the occurrence and progression of many human diseases, including cancer, fibrosis, and inflammation. Therefore, DDR1 represents a highly promising therapeutic target. Although no selective small-molecule inhibitors have reached clinical trials to date, many molecules have shown therapeutic effects in preclinical studies. For example, BK40143 has demonstrated significant promise in the therapy of neurodegenerative diseases. In this context, our perspective aims to provide an in-depth exploration of DDR1, encompassing its structure characteristics, biological functions, and disease relevance. Furthermore, we emphasize the importance of understanding the structure-activity relationship of DDR1 inhibitors and highlight the unique advantages of dual-target or multitarget inhibitors. We anticipate offering valuable insights into the development of more efficacious DDR1-targeted drugs.

Acceleration of Molecular Simulations by Parametric Time-Lagged tSNE Metadynamics

The potential of molecular simulations is limited by their computational costs. There is often a need to accelerate simulations using some of the enhanced sampling methods. Metadynamics applies a history-dependent bias potential that disfavors previously visited states. To apply metadynamics, it is necessary to select a few properties of the system─collective variables (CVs) that can be used to define the bias potential. Over the past few years, there have been emerging opportunities for machine learning and, in particular, artificial neural networks within this domain. In this broad context, a specific unsupervised machine learning method was utilized, namely, parametric time-lagged t-distributed stochastic neighbor embedding (ptltSNE) to design CVs. The approach was tested on a Trp-cage trajectory (tryptophan cage) from the literature. The trajectory was used to generate a map of conformations, distinguish fast conformational changes from slow ones, and design CVs. Then, metadynamic simulations were performed. To accelerate the formation of the α-helix, we added the α-RMSD collective variable. This simulation led to one folding event in a 350 ns metadynamics simulation. To accelerate degrees of freedom not addressed by CVs, we performed parallel tempering metadynamics. This simulation led to 10 folding events in a 200 ns simulation with 32 replicas.

Analyzing Molecular Dynamics Trajectories Thermodynamically through Artificial Intelligence

Molecular dynamics simulations produce trajectories that correspond to vast amounts of structure when exploring biochemical processes. Extracting valuable information, e.g., important intermediate states and collective variables (CVs) that describe the major movement modes, from molecular trajectories to understand the underlying mechanisms of biological processes presents a significant challenge. To achieve this goal, we introduce a deep learning approach, coined DIKI (deep identification of key intermediates), to determine low-dimensional CVs distinguishing key intermediate conformations without a-priori assumptions. DIKI dynamically plans the distribution of latent space and groups together similar conformations within the same cluster. Moreover, by incorporating two user-defined parameters, namely, coarse focus knob and fine focus knob, to help identify conformations with low free energy and differentiate the subtle distinctions among these conformations, resolution-tunable clustering was achieved. Furthermore, the integration of DIKI with a path-finding algorithm contributes to the identification of crucial intermediates along the lowest free-energy pathway. We postulate that DIKI is a robust and flexible tool that can find widespread applications in the analysis of complex biochemical processes.

Construction and validation of a prognostic signature based on seven endoplasmic reticulum stress-related lncRNAs for patients with head and neck squamous cell carcinoma

Endoplasmic reticulum stress (ERS) occurs when misfolded or unfolded proteins accumulate in the endoplasmic reticulum (ER), and it is often observed in tumors, including head and neck squamous cell carcinoma (HNSCC). Relevant studies have demonstrated the prognostic significance of ERS-related long non-coding RNAs (lncRNAs) in various cancers. However, the relationship between ERS and lncRNAs in HNSCC has received limited attention in previous studies. In this study, we aimed to develop an ERS-related lncRNAs prognostic model using correlation analysis, Cox regression analysis, least absolute shrinkage, and selection operator (LASSO) regression analysis based on data from The Cancer Genome Atlas (TCGA) database. The survival and predictive ability of this model were evaluated using Kaplan-Meier analysis and time-dependent receiver operating characteristics (ROC), while nomograms and calibration curves were constructed. Then, functional enrichment analyses, tumor mutation burden (TMB), tumor infiltration of immune cells, single sample Gene Set Enrichment Analysis (ssGSEA), and drug sensitivity analysis were performed. Additionally, we conducted a consensus cluster analysis to compare differences between subtypes of tumors. Finally, we validated the expression of the ERS-related lncRNAs that constructed prognostic risk score model in HNSCC tissues through quantitative real-time PCR (qRT-PCR). We developed a prognostic signature based on seven ERS-related lncRNAs, which showed better predictive performance than other clinicopathological features. The high-risk poor prognosis group had a poorer prognosis in comparison to the low-risk good prognosis. The area under the ROC curve (AUC) predicted by this model for 3-year survival rates of HNSCC patients was 0.805. Enrichment analysis revealed that the differentially expressed genes were primarily enriched in pathways related to immune responses and signal transduction. Low-risk patients had lower TMB, more immune cell infiltrations, and enhanced anti-tumor immunity. Cluster analysis indicated that cluster 3 may have a better prognosis and immunotherapy effect. In addition, the result of qRT-PCR was consistent with our analysis. This prognostic model based on seven ERS-related lncRNAs is a promising tool for risk stratification, survival prediction, and immune cell infiltration status assessment.

Application of a Two-Dimensional Mapping-Based Visualization Technique: Nutrient-Value-Based Food Grouping

Worldwide, several food-based dietary guidelines, with diverse food-grouping methods in various countries, have been developed to maintain and promote public health. However, standardized international food-grouping methods are scarce. In this study, we used two-dimensional mapping to classify foods based on their nutrient composition. The Standard Tables of Food Composition in Japan were used for mapping with a novel technique-t-distributed stochastic neighbor embedding-to visualize high-dimensional data. The mapping results showed that most foods formed food group-based clusters in the Standard Tables of Food Composition in Japan. However, the beverages did not form large clusters and demonstrated scattered distribution on the map. Green tea, black tea, and coffee are located within or near the vegetable cluster whereas cocoa is near the pulse cluster. These results were ensured by the k-nearest neighbors. Thus, beverages made from natural materials can be categorized based on their origin. Visualization of food composition could enable an enhanced comprehensive understanding of the nutrients in foods, which could lead to novel aspects of nutrient-value-based food classifications.

CENsible: Interpretable Insights into Small-Molecule Binding with Context Explanation Networks

We present a novel and interpretable approach for predicting small-molecule binding affinities using context explanation networks (CENs). Given the specific structure of a protein/ligand complex, our CENsible scoring function uses a deep convolutional neural network to predict the contributions of pre-calculated terms to the overall binding affinity. We show that CENsible can effectively distinguish active vs. inactive compounds for many systems. Its primary benefit over related machine-learning scoring functions, however, is that it retains interpretability, allowing researchers to identify the contribution of each pre-calculated term to the final affinity prediction, with implications for subsequent lead optimization.

Clustering Heterogeneous Conformational Ensembles of Intrinsically Disordered Proteins with t-Distributed Stochastic Neighbor Embedding

Intrinsically disordered proteins (IDPs) populate a range of conformations that are best described by a heterogeneous ensemble. Grouping an IDP ensemble into "structurally similar" clusters for visualization, interpretation, and analysis purposes is a much-desired but formidable task, as the conformational space of IDPs is inherently high-dimensional and reduction techniques often result in ambiguous classifications. Here, we employ the t-distributed stochastic neighbor embedding (t-SNE) technique to generate homogeneous clusters of IDP conformations from the full heterogeneous ensemble. We illustrate the utility of t-SNE by clustering conformations of two disordered proteins, Aβ42, and α-synuclein, in their APO states and when bound to small molecule ligands. Our results shed light on ordered substates within disordered ensembles and provide structural and mechanistic insights into binding modes that confer specificity and affinity in IDP ligand binding. t-SNE projections preserve the local neighborhood information, provide interpretable visualizations of the conformational heterogeneity within each ensemble, and enable the quantification of cluster populations and their relative shifts upon ligand binding. Our approach provides a new framework for detailed investigations of the thermodynamics and kinetics of IDP ligand binding and will aid rational drug design for IDPs.

Density Functional Theory Computation and Machine Learning Studies of Interaction between Au3 Clusters and 20 Natural Amino Acid Molecules

The optimal adsorption sites and the binding energies of neutral Au₃ clusters with 20 natural amino acids under the gas phase and water solvation were systematically investigated based on density functional theory (DFT) calculations. The calculation results showed that in the gas phase Au₃ tends to bind with N atoms of amino groups in amino acids, except methionine, which tends to bind with Au₃ through S atoms. Under water solvation, Au₃ clusters tended to bind to N atoms of amino groups and N atoms of side chain amino groups in amino acids. However, methionine and cysteine bind more strongly to the gold atom through the S atom. Based on the binding energy data of Au₃ clusters and 20 natural amino acids under water solvation calculated by DFT, a machine learning model (gradient boosted decision tree) was proposed to predict the optimal binding Gibbs free energy (ΔG) of the interaction between Au₃ clusters and amino acids. The main factors affecting the strength of the interaction between Au₃ and amino acids were uncovered by the feature importance analysis.

Immune-related gene signature associates with immune landscape and predicts prognosis accurately in patients with skin cutaneous melanoma

Skin cutaneous melanoma (SKCM) is the skin cancer that causes the highest number of deaths worldwide. There is growing evidence that the tumour immune microenvironment is associated with cancer prognosis, however, there is little research on the role of immune status in melanoma prognosis. In this study, data on patients with Skin cutaneous melanoma were downloaded from the GEO, TCGA, and GTEx databases. Genes associated with the immune pathway were screened from published papers and lncRNAs associated with them were identified. We performed immune microenvironment and functional enrichment analyses. The analysis was followed by applying univariate/multivariate Cox regression algorithms to finally identify three lncRNAs associated with the immune pathway for the construction of prognostic prediction models (CXCL10, RXRG, and SCG2). This stepwise downscaling method, which finally screens out prognostic factors and key genes and then uses them to build a risk model, has excellent predictive power. According to analyses of the model's reliability, it was able to differentiate the prognostic value and continued existence of Skin cutaneous melanoma patient populations more effectively. This study is an analysis of the immune pathway that leads lncRNAs in Skin cutaneous melanoma in an effort to open up new treatment avenues for Skin cutaneous melanoma.

A study of autoencoders as a feature extraction technique for spike sorting

Spike sorting is the process of grouping spikes of distinct neurons into their respective clusters. Most frequently, this grouping is performed by relying on the similarity of features extracted from spike shapes. In spite of recent developments, current methods have yet to achieve satisfactory performance and many investigators favour sorting manually, even though it is an intensive undertaking that requires prolonged allotments of time. To automate the process, a diverse array of machine learning techniques has been applied. The performance of these techniques depends however critically on the feature extraction step. Here, we propose deep learning using autoencoders as a feature extraction method and evaluate extensively the performance of multiple designs. The models presented are evaluated on publicly available synthetic and real "in vivo" datasets, with various numbers of clusters. The proposed methods indicate a higher performance for the process of spike sorting when compared to other state-of-the-art techniques.

Unsupervised machine learning framework for discriminating major variants of concern during COVID-19

Due to the high mutation rate of the virus, the COVID-19 pandemic evolved rapidly. Certain variants of the virus, such as Delta and Omicron emerged with altered viral properties leading to severe transmission and death rates. These variants burdened the medical systems worldwide with a major impact to travel, productivity, and the world economy. Unsupervised machine learning methods have the ability to compress, characterize, and visualize unlabelled data. This paper presents a framework that utilizes unsupervised machine learning methods to discriminate and visualize the associations between major COVID-19 variants based on their genome sequences. These methods comprise a combination of selected dimensionality reduction and clustering techniques. The framework processes the RNA sequences by performing a k-mer analysis on the data and further visualises and compares the results using selected dimensionality reduction methods that include principal component analysis (PCA), t-distributed stochastic neighbour embedding (t-SNE), and uniform manifold approximation projection (UMAP). Our framework also employs agglomerative hierarchical clustering to visualize the mutational differences among major variants of concern and country-wise mutational differences for selected variants (Delta and Omicron) using dendrograms. We also provide country-wise mutational differences for selected variants via dendrograms. We find that the proposed framework can effectively distinguish between the major variants and has the potential to identify emerging variants in the future.

CNN6mA: Interpretable neural network model based on position-specific CNN and cross-interactive network for 6mA site prediction

N6-methyladenine (6mA) plays a critical role in various epigenetic processing including DNA replication, DNA repair, silencing, transcription, and diseases such as cancer. To understand such epigenetic mechanisms, 6 mA has been detected by high-throughput technologies on a genome-wide scale at single-base resolution, together with conventional methods such as immunoprecipitation, mass spectrometry and capillary electrophoresis, but these experimental approaches are time-consuming and laborious. To complement these problems, we have developed a CNN-based 6 mA site predictor, named CNN6mA, which proposed two new architectures: a position-specific 1-D convolutional layer and a cross-interactive network. In the position-specific 1-D convolutional layer, position-specific filters with different window sizes were applied to an inquiry sequence instead of sharing the same filters over all positions in order to extract the position-specific features at different levels. The cross-interactive network explored the relationships between all the nucleotide patterns within the inquiry sequence. Consequently, CNN6mA outperformed the existing state-of-the-art models in many species and created the contribution score vector that intelligibly interpret the prediction mechanism. The source codes and web application in CNN6mA are freely accessible at https://github.com/kuratahiroyuki/CNN6mA.git and http://kurata35.bio.kyutech.ac.jp/CNN6mA/, respectively.

pH Induced Switch in the Conformational Ensemble of Intrinsically Disordered Protein Prothymosin-α and Its Implications for Amyloid Fibril Formation

Aggregation of intrinsically disordered proteins (IDPs) can lead to neurodegenerative diseases. Although there is experimental evidence that acidic pH promotes IDP monomer compaction leading to aggregation, the general mechanism is unclear. We studied the pH effect on the conformational ensemble of prothymosin-α (proTα), which is involved in multiple essential functions, and probed its role in aggregation using computer simulations. We show that compaction in the proTα dimension at low pH is due to the protein's collapse in the intermediate region (E41-D80) rich in glutamic acid residues, enhancing its β-sheet content. We observed by performing dimer simulations that the conformations with high β-sheet content could act as aggregation-prone (N*) states and nucleate the aggregation process. The simulations initiated using N* states form dimers within a microsecond time scale, whereas the non-N* states do not form dimers within this time scale. This study contributes to understanding the general principles of pH-induced IDP aggregation.

An eleven autophagy-related genes-based prognostic signature for endometrial carcinoma

Background

Endometrial cancer (EC) is a common malignant tumor in women with increasing mortality. The prognosis of EC is highly heterogeneous which needs more effective biomarkers for clinical decision. Here, we reported the effect of autophagy-related genes (ARGs) on the prognosis of EC.

Methods

The expression data of EC tissues and adjacent non-tumor samples were available from the TCGA dataset and 232 autophagy-related genes were from The Human Autophagy Database. A prognostic ARGs risk model was further constructed by using LASSO-Cox regression, and its prognostic and predictive value were evaluated by nomogram. Further functional analysis was conducted to reveal a significant signaling pathway.

Results

A total of 45 differentially expressed ARGs were obtained, including 18 upregulated and 27 downregulated genes. Eleven ARGs (BID, CAPN2, CDKN2A, DLC1, GRID2, IFNG, MYC, NRG3, P4HB, PTK6, and TP73) were finally selected to build ARGs risk. This signature could well distinguish between the high- and low-risk patients (survival analysis: P = 1.18E-10; AUC: 0.733 at 1 year, 0.795 at 3 years, and 0.823 at 5 years). Furthermore, a nomogram was plotting to predict the possibility of overall survival and suggested good value for clinical utility.

Conclusion

We established an eleven-ARG signature, which was probably effective in the prognostic prediction of patients with EC.

Prediction of Aptamer-Small-Molecule Interactions Using Metastable States from Multiple Independent Molecular Dynamics Simulations

Understanding aptamer-ligand interactions is necessary to rationally design aptamer-based systems. Commonly used in silico tools have proven to be accurate to predict RNA and DNA oligonucleotide tertiary structures. However, given the complexity of nucleic acids, the most thermodynamically stable conformation is not necessarily the one with the highest affinity for a specific ligand. Because many metastable states may coexist, it remains challenging to predict binding sites through molecular docking simulations using available computational pipelines. In this study, we used independent simulations to broaden the conformational diversity sampled from DNA initial models of distinct stability and assessed the binding affinity of selected metastable representative structures. In our results, utilizing multiple metastable conformations for molecular docking analysis helped identify structures favorable for ligand binding and accurately predict the binding sites. Our workflow was able to correctly identify the binding sites of the characterized adenosine monophosphate and l-argininamide aptamers. Additionally, we demonstrated that our pipeline can be used to aid the design of competition assays that are conducive to aptasensing strategies using an uncharacterized aflatoxin B1 aptamer. We foresee that this approach may help rationally design effective and truncated aptamer sequences interacting with protein biomarkers or small molecules of interest for drug design and sensor applications.

Sparse group selection and analysis of function-related residue for protein-state recognition

Machine learning methods have helped to advance wide range of scientific and technological field in recent years, including computational chemistry. As the chemical systems could become complex with high dimension, feature selection could be critical but challenging to develop reliable machine learning based prediction models, especially for proteins as bio-macromolecules. In this study, we applied sparse group lasso (SGL) method as a general feature selection method to develop classification model for an allosteric protein in different functional states. This results into a much improved model with comparable accuracy (Acc) and only 28 selected features comparing to 289 selected features from a previous study. The Acc achieves 91.50% with 1936 selected feature, which is far higher than that of baseline methods. In addition, grouping protein amino acids into secondary structures provides additional interpretability of the selected features. The selected features are verified as associated with key allosteric residues through comparison with both experimental and computational works about the model protein, and demonstrate the effectiveness and necessity of applying rigorous feature selection and evaluation methods on complex chemical systems.

From Data to Knowledge: Systematic Review of Tools for Automatic Analysis of Molecular Dynamics Output

An increasing number of crystal structures available on one side, and the boost of computational power available for computer-aided drug design tasks on the other, have caused that the structure-based drug design tools are intensively used in the drug development pipelines. Docking and molecular dynamics simulations, key representatives of the structure-based approaches, provide detailed information about the potential interaction of a ligand with a target receptor. However, at the same time, they require a three-dimensional structure of a protein and a relatively high amount of computational resources. Nowadays, as both docking and molecular dynamics are much more extensively used, the amount of data output from these procedures is also growing. Therefore, there are also more and more approaches that facilitate the analysis and interpretation of the results of structure-based tools. In this review, we will comprehensively summarize approaches for handling molecular dynamics simulations output. It will cover both statistical and machine-learning-based tools, as well as various forms of depiction of molecular dynamics output.

A hybrid deep learning model for forecasting lymphocyte depletion during radiation therapy

PURPOSE

Recent studies have shown that severe depletion of the absolute lymphocyte count (ALC) induced by radiation therapy (RT) has been associated with poor overall survival of patients with many solid tumors. In this paper, we aimed to predict radiation-induced lymphocyte depletion in esophageal cancer patients during the course of RT based on patient characteristics and dosimetric features.

METHODS

We proposed a hybrid deep learning model in a stacked structure to predict a trend toward ALC depletion based on the clinical information before or at the early stages of RT treatment. The proposed model consisted of four channels, one channel based on long short-term memory (LSTM) network and three channels based on neural networks, to process four categories of features followed by a dense layer to integrate the outputs of four channels and predict the weekly ALC values. Moreover, a discriminative kernel was developed to extract temporal features and assign different weights to each part of the input sequence which enabled the model to focus on the most relevant parts. The proposed model was trained and tested on a dataset of 860 esophageal cancer patients who received concurrent chemoradiotherapy.

RESULTS

The performance of the proposed model was evaluated based on several important prediction metrics and compared to other commonly used prediction models. The results showed that the proposed model outperformed off-the-shelf prediction methods with at least a 30% reduction in the mean squared error (MSE) of weekly ALC predictions based on pretreatment data.Moreover, using an extended model based on augmented first-week treatment data reduced the MSE of predictions by 70% compared to the model based on the pretreatment data.

CONCLUSIONS

In conclusion, our model performed well in predicting radiation-induced lymphocyte depletion for RT treatment planning. The ability to predict ALC will enable physicians to evaluate individual RT treatment plans for lymphopenia risk and to identify patients at high risk who would benefit from modified treatment approaches. This article is protected by copyright. All rights reserved.

A Hypoxia-Related Long Non-Coding RNAs Signature Associated With Prognosis in Lower-Grade Glioma

Accumulating evidence suggests that hypoxia microenvironment and long non-coding lncRNAs (lncRNAs) exert critical roles in tumor development. Herein, we aim to develop a hypoxia-related lncRNA (HRL) model to predict the survival outcomes of patient with lower-grade glioma (LGG). The RNA-sequencing data of 505 LGG samples were acquired from The Cancer Genome Atlas (TCGA). Using consensus clustering based on the expression of hypoxia-related mRNAs, these samples were divided into three subsets that exhibit distinct hypoxia content, clinicopathologic features, and survival status. The differentially expressed lncRNAs across the subgroups were documented as candidate HRLs. With LASSO regression analysis, eight informative lncRNAs were selected for constructing the prognostic HRL model. This signature had a good performance in predicting LGG patients’ overall survival in the TCGA cohort, and similar results could be achieved in two validation cohorts from the Chinese Glioma Genome Atlas. The HRL model also showed correlations with important clinicopathologic characteristics such as patients’ age, tumor grade, IDH mutation, 1p/19q codeletion, MGMT methylation, and tumor progression risk. Functional enrichment analysis indicated that the HLR signature was mainly involved in regulation of inflammatory response, complement, hypoxia, Kras signaling, and apical junction. More importantly, the signature was related to immune cell infiltration, estimated immune score, tumor mutation burden, neoantigen load, and expressions of immune checkpoints and immunosuppressive cytokines. Finally, a nomogram was developed by integrating the HRL signature and clinicopathologic features, with a concordance index of 0.852 to estimate the survival probability of LGG patients. In conclusion, our study established an effective HRL model for prognosis assessment of LGG patients, which may provide insights for future research and facilitate the designing of individualized treatment.

Exploration and augmentation of pharmacological space via adversarial auto-encoder model for facilitating kinase-centric drug development

Predicting compound-protein interactions (CPIs) is of great importance for drug discovery and repositioning, yet still challenging mainly due to the sparse nature of CPI matrixes, resulting in poor generalization performance. Hence, unlike typical CPI prediction models focused on representation learning or model selection, we propose a deep neural network-based strategy, PCM-AAE, that re-explores and augments the pharmacological space of kinase inhibitors by introducing the adversarial auto-encoder model (AAE) to improve the generalization of the prediction model. To complete the data space, we constructed Ensemble of PCM-AAE (EPA), an ensemble model that quickly and accurately yields quantitative predictions of binding affinity between any human kinase and inhibitor. In rigorous internal validation, EPA showed excellent performance, consistently outperforming the model trained with the imbalanced set, especially for targets with relatively fewer training data points. Improved prediction accuracy of EPA for external datasets enhances its generalization ability, making it possible to gracefully handle previously unseen kinases and inhibitors. EPA showed promising potential when directly applied to virtual screening and off-target prediction, exhibiting its practicality in hit prediction. Our strategy is expected to facilitate kinase-centric drug development, as well as to solve more challenging prediction problems with insufficient data points.

Single-cell transcriptome identifies molecular subtype of autism spectrum disorder impacted by de novo loss-of-function variants regulating glial cells

BACKGROUND

In recent years, several hundred autism spectrum disorder (ASD) implicated genes have been discovered impacting a wide range of molecular pathways. However, the molecular underpinning of ASD, particularly from the point of view of 'brain to behaviour' pathogenic mechanisms, remains largely unknown.

METHODS

We undertook a study to investigate patterns of spatiotemporal and cell type expression of ASD-implicated genes by integrating large-scale brain single-cell transcriptomes (> million cells) and de novo loss-of-function (LOF) ASD variants (impacting 852 genes from 40,122 cases).

RESULTS

We identified multiple single-cell clusters from three distinct developmental human brain regions (anterior cingulate cortex, middle temporal gyrus and primary visual cortex) that evidenced high evolutionary constraint through enrichment for brain critical exons and high pLI genes. These clusters also showed significant enrichment with ASD loss-of-function variant genes (p < 5.23 × 10^-11) that are transcriptionally highly active in prenatal brain regions (visual cortex and dorsolateral prefrontal cortex). Mapping ASD de novo LOF variant genes into large-scale human and mouse brain single-cell transcriptome analysis demonstrate enrichment of such genes into neuronal subtypes and are also enriched for subtype of non-neuronal glial cell types (astrocyte, p < 6.40 × 10^-11, oligodendrocyte, p < 1.31 × 10^-09).

CONCLUSION

Among the ASD genes enriched with pathogenic de novo LOF variants (i.e. KANK1, PLXNB1), a subgroup has restricted transcriptional regulation in non-neuronal cell types that are evolutionarily conserved. This association strongly suggests the involvement of subtype of non-neuronal glial cells in the pathogenesis of ASD and the need to explore other biological pathways for this disorder.

Single-Cell Proteomic Profiling Identifies Nanoparticle Enhanced Therapy for Triple Negative Breast Cancer Stem Cells

Breast cancer remains a major cause of cancer-related deaths in women worldwide. Chemotherapy-promoted stemness and enhanced stem cell plasticity in breast cancer is a cause for great concern. The discovery of drugs targeting BCSCs was suggested to be an important advancement in the establishment of therapy that improves the efficacy of chemotherapy. In this work, by using single-cell mass cytometry, we observed that stemness in spheroid-forming cells derived from MDA-MB-231 cells was significantly increased after doxorubicin administration and up-regulated integrin αvβ3 expression was also observed. An RGD-included nanoparticle (CS-V) was designed, and it was found that it could promote doxorubicin’s efficacy against MDA-MB-231 spheroid cells. The above observations suggested that the combination of RGD-included nanoparticles (CS-V) with the chemo-drug doxorubicin could be developed as a potential therapy for breast cancer.

Collective variable-based enhanced sampling and machine learning

ABSTRACT

Collective variable-based enhanced sampling methods have been widely used to study thermodynamic properties of complex systems. Efficiency and accuracy of these enhanced sampling methods are affected by two factors: constructing appropriate collective variables for enhanced sampling and generating accurate free energy surfaces. Recently, many machine learning techniques have been developed to improve the quality of collective variables and the accuracy of free energy surfaces. Although machine learning has achieved great successes in improving enhanced sampling methods, there are still many challenges and open questions. In this perspective, we shall review recent developments on integrating machine learning techniques and collective variable-based enhanced sampling approaches. We also discuss challenges and future research directions including generating kinetic information, exploring high-dimensional free energy surfaces, and efficiently sampling all-atom configurations.

GRAPHIC ABSTRACT

AKT inhibitor AZD5363 suppresses stemness and promotes anti-cancer activity of 3,3'-diindolylmethane in human breast cancer cells

3,3'-diindolylmethane (DIM) is a dimer compound converted from Indoly-3-carbinol that had been studied as promising chemo-preventive agent against breast cancer. In this study, we observed that proportion of CD133⁺Nanog⁺ subpopulation in MCF-7 cells was significantly increased after DIM administration with up-regulated AKT activity by using CyTOF assay. SPADE analysis revealed this stem-like subpopulation exhibited apoptosis-resistance property against DIM treatment. By combining with AKT inhibitor AZD5363, DIM induced CD133 expression could be suppressed. In addition, a combination treatment of MCF-7 and MDA-MB-231 breast cancer cells with DIM and AZD5363 showed synergistic decreases in cell proliferation and induced apoptosis. Furthermore, results from imaging flow cytometry suggested that FoxO3a nuclear localization and PUMA expression could be improved by combination of AZD5363 with DIM. Taken together, the above observations suggested that the combination of AZD5363 with DIM could be developed as potential therapy for breast cancer.

Unsupervised Learning Methods for Molecular Simulation Data

Unsupervised learning is becoming an essential tool to analyze the increasingly large amounts of data produced by atomistic and molecular simulations, in material science, solid state physics, biophysics, and biochemistry. In this Review, we provide a comprehensive overview of the methods of unsupervised learning that have been most commonly used to investigate simulation data and indicate likely directions for further developments in the field. In particular, we discuss feature representation of molecular systems and present state-of-the-art algorithms of dimensionality reduction, density estimation, and clustering, and kinetic models. We divide our discussion into self-contained sections, each discussing a specific method. In each section, we briefly touch upon the mathematical and algorithmic foundations of the method, highlight its strengths and limitations, and describe the specific ways in which it has been used-or can be used-to analyze molecular simulation data.

Multiscale Reweighted Stochastic Embedding: Deep Learning of Collective Variables for Enhanced Sampling

Machine learning methods provide a general framework for automatically finding and representing the essential characteristics of simulation data. This task is particularly crucial in enhanced sampling simulations. There we seek a few generalized degrees of freedom, referred to as collective variables (CVs), to represent and drive the sampling of the free energy landscape. In theory, these CVs should separate different metastable states and correspond to the slow degrees of freedom of the studied physical process. To this aim, we propose a new method that we call multiscale reweighted stochastic embedding (MRSE). Our work builds upon a parametric version of stochastic neighbor embedding. The technique automatically learns CVs that map a high-dimensional feature space to a low-dimensional latent space via a deep neural network. We introduce several new advancements to stochastic neighbor embedding methods that make MRSE especially suitable for enhanced sampling simulations: (1) weight-tempered random sampling as a landmark selection scheme to obtain training data sets that strike a balance between equilibrium representation and capturing important metastable states lying higher in free energy; (2) a multiscale representation of the high-dimensional feature space via a Gaussian mixture probability model; and (3) a reweighting procedure to account for training data from a biased probability distribution. We show that MRSE constructs low-dimensional CVs that can correctly characterize the different metastable states in three model systems: the Müller-Brown potential, alanine dipeptide, and alanine tetrapeptide.

Dimeric allostery mechanism of the plant circadian clock photoreceptor ZEITLUPE

In Arabidopsis thaliana, the Light-Oxygen-Voltage (LOV) domain containing protein ZEITLUPE (ZTL) integrates light quality, intensity, and duration into regulation of the circadian clock. Recent structural and biochemical studies of ZTL indicate that the protein diverges from other members of the LOV superfamily in its allosteric mechanism, and that the divergent allosteric mechanism hinges upon conservation of two signaling residues G46 and V48 that alter dynamic motions of a Gln residue implicated in signal transduction in all LOV proteins. Here, we delineate the allosteric mechanism of ZTL via an integrated computational approach that employs atomistic simulations of wild type and allosteric variants of ZTL in the functional dark and light states, together with Markov state and supervised machine learning classification models. This approach has unveiled key factors of the ZTL allosteric mechanisms, and identified specific interactions and residues implicated in functional allosteric changes. The final results reveal atomic level insights into allosteric mechanisms of ZTL function that operate via a non-trivial combination of population-shift and dynamics-driven allosteric pathways.

CyAnno: A semi-automated approach for cell type annotation of mass cytometry datasets

MOTIVATION

For immune system monitoring in large-scale studies at the single-cell resolution using CyTOF, (semi-)automated computational methods are applied for annotating live cells of mixed cell types. Here, we show that the live cell pool can be highly enriched with undefined heterogeneous cells, i.e., 'ungated' cells, and that current semi-automated approaches ignore their modeling resulting in misclassified annotations.

RESULT

We introduce 'CyAnno', a novel semi-automated approach for deconvoluting the unlabeled cytometry dataset based on a machine learning framework utilizing manually gated training data that allows the integrative modeling of 'gated' cell types and the 'ungated' cells. By applying this framework on several CyTOF datasets, we demonstrated that including the 'ungated' cells can lead to a significant increase in the precision of the 'gated' cell types prediction. CyAnno can be used to identify even a single cell type, including rare cells, with higher efficacy than current state-of-the-art semi-automated approaches.

AVAILABILITY

The CyAnno is available as a python script with a user-manual and sample dataset at https://github.com/abbioinfo/CyAnno.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

UMAP as a Dimensionality Reduction Tool for Molecular Dynamics Simulations of Biomacromolecules: A Comparison Study

Proteins are the molecular machines of life. The multitude of possible conformations that proteins can adopt determines their free-energy landscapes. However, the inherently high dimensionality of a protein free-energy landscape poses a challenge to deciphering how proteins perform their functions. For this reason, dimensionality reduction is an active field of research for molecular biologists. The uniform manifold approximation and projection (UMAP) is a dimensionality reduction method based on a fuzzy topological analysis of data. In the present study, the performance of UMAP is compared with that of other popular dimensionality reduction methods such as t-distributed stochastic neighbor embedding (t-SNE), principal component analysis (PCA), and time-structure independent components analysis (tICA) in the context of analyzing molecular dynamics simulations of the circadian clock protein VIVID. A good dimensionality reduction method should accurately represent the data structure on the projected components. The comparison of the raw high-dimensional data with the projections obtained using different dimensionality reduction methods based on various metrics showed that UMAP has superior performance when compared with linear reduction methods (PCA and tICA) and has competitive performance and scalable computational cost.

Pepsin-like aspartic proteases (PAPs) as model systems for combining biomolecular simulation with biophysical experiments

Pepsin-like aspartic proteases (PAPs) are a class of aspartic proteases which shares tremendous structural similarity with human pepsin. One of the key structural features of PAPs is the presence of a β-hairpin motif otherwise known as flap. The biological function of the PAPs is highly dependent on the conformational dynamics of the flap region. In apo PAPs, the conformational dynamics of the flap is dominated by the rotational degrees of freedom associated with χ1 and χ2 angles of conserved Tyr (or Phe in some cases). However it is plausible that dihedral order parameters associated with several other residues might play crucial roles in the conformational dynamics of apo PAPs. Due to their size, complexities associated with conformational dynamics and clinical significance (drug targets for malaria, Alzheimer's disease etc.), PAPs provide a challenging testing ground for computational and experimental methods focusing on understanding conformational dynamics and molecular recognition in biomolecules. The opening of the flap region is necessary to accommodate substrate/ligand in the active site of the PAPs. The BIG challenge is to gain atomistic details into how reversible ligand binding/unbinding (molecular recognition) affects the conformational dynamics. Recent reports of kinetics (K i, K d) and thermodynamic parameters (ΔH, TΔS, and ΔG) associated with macro-cyclic ligands bound to BACE1 (belongs to PAP family) provide a perfect challenge (how to deal with big ligands with multiple torsional angles and select optimum order parameters to study reversible ligand binding/unbinding) for computational methods to predict binding free energies and kinetics beyond typical test systems e.g. benzamide-trypsin. In this work, i reviewed several order parameters which were proposed to capture the conformational dynamics and molecular recognition in PAPs. I further highlighted how machine learning methods can be used as order parameters in the context of PAPs. I then proposed some open ideas and challenges in the context of molecular simulation and put forward my case on how biophysical experiments e.g. NMR, time-resolved FRET etc. can be used in conjunction with biomolecular simulation to gain complete atomistic insights into the conformational dynamics of PAPs.

Three Immune-Associated Subtypes of Diffuse Glioma Differ in Immune Infiltration, Immune Checkpoint Molecules, and Prognosis

Diffuse glioma is one of the most prevalent malignancies of the brain, with high heterogeneity of tumor-infiltrating immune cells. However, immune-associated subtypes of diffuse glioma have not been determined, nor has the effect of different immune-associated subtypes on disease prognosis and immune infiltration of diffuse glioma patients. We retrieved the expression profiles of immune-related genes from The Cancer Genome Atlas (TCGA) (n = 672) and GSE16011 (n = 268) cohorts and used them to identify subtypes of diffuse glioma via Consensus Cluster Plus analysis. We used the limma, clusterProfiler, ESTIMATE, and survival packages of R for differential analysis, functional enrichment, immune and stromal score evaluation respectively in three subtypes, and performed log-rank tests in immune subtypes of diffuse glioma. The immune-associated features of diffuse glioma in the two cohorts were characterized via bioinformatic analyses of the mRNA expression data of immune-related genes. Three subtypes (C1–3) of diffuse glioma were identified from TCGA data, and were verified using the GSE16011 cohort. We then evaluated their immune characteristics and clinical features. Our mRNA profiling analyses indicated that the different subtypes of diffuse glioma presented differential expression profile of specific genes and signal pathways in the TCGA cohort. Patients with subtype C1, who were mostly diagnosed with grade IV glioma, had poorer outcomes than patients with subtype C2 or C3. Subtype C1 was characterized by a higher degree of immune cell infiltration as estimated by GSVA, and more frequent wildtype IDH1. By contrast, subtype C3 included more grade II and IDH1-mutated glioma, and was associated with more infiltration of CD4⁺T cells. Most subtype C2 had the features between subtypes C1 and C3. Meanwhile, immune checkpoints and their ligand molecules, including PD1/(PD-L1/PDL2), CTLA4/(CD80/CD86), and B7H3/TLT2, were significantly upregulated in subtype C1 and downregulated in subtype C3. In addition, patients with subtype C1 exhibited more frequent gene mutations. Univariate and multivariate Cox regression analyses revealed that diffuse glioma subtype was an effective, independent, and better prognostic factor. Therefore, we established a novel immune-related classification of diffuse glioma, which provides potential immunotherapy targets for diffuse glioma.

Conformational Landscapes of Halohydrin Dehalogenases and Their Accessible Active Site Tunnels

Halohydrin dehalogenases (HHDH) are industrially relevant biocatalysts exhibiting a promiscuous epoxide-ring opening reactivity in the presence of small nucleophiles, thus giving access to novel carbon–carbon, carbon–oxygen, carbon–nitrogen, and carbon–sulfur bonds. Recently, the repertoire of HHDH has been expanded, providing access to some novel HHDH subclasses exhibiting a broader epoxide substrate scope. In this work, we develop a computational approach based on the application of linear and non-linear dimensionality reduction techniques to long time-scale Molecular Dynamics (MD) simulations to study the HHDH conformational landscapes. We couple the analysis of the conformational landscapes to CAVER calculations to assess their impact on the active site tunnels and potential ability towards bulky epoxide ring opening reaction. Our study indicates that the analyzed HHDHs subclasses share a common breathing motion of the halide binding pocket, but present large deviations in the loops adjacent to the active site pocket and N-terminal regions. Such conformational differences affect the available tunnels for epoxide binding to the active site. The superior activity of the HHDH G subclass towards bulkier substrates is explained by the additional structural elements delimiting the active site region, its rich conformational heterogeneity, and the substantially wider and frequently observed active site tunnels. This study therefore provides key information for HHDH promiscuity and engineering.

Utilizing a multi-class classification approach to detect therapeutic and recreational misuse of opioids on Twitter

BACKGROUND

Opioid misuse (OM) is a major health problem in the United States, and can lead to addiction and fatal overdose. We sought to employ natural language processing (NLP) and machine learning to categorize Twitter chatter based on the motive of OM.

MATERIALS AND METHODS

We collected data from Twitter using opioid-related keywords, and manually annotated 6988 tweets into three classes-No-OM, Pain-related-OM, and Recreational-OM-with the No-OM class representing tweets indicating no use/misuse, and the Pain-related misuse and Recreational-misuse classes representing misuse for pain or recreation/addiction. We trained and evaluated multi-class classifiers, and performed term-level k-means clustering to assess whether there were terms closely associated with the three classes.

RESULTS

On a held-out test set of 1677 tweets, a transformer-based classifier (XLNet) achieved the best performance with F₁-score of 0.71 for the Pain-misuse class, and 0.79 for the Recreational-misuse class. Macro- and micro-averaged F₁-scores over all classes were 0.82 and 0.92, respectively. Content-analysis using clustering revealed distinct clusters of terms associated with each class.

DISCUSSION

While some past studies have attempted to automatically detect opioid misuse, none have further characterized the motive for misuse. Our multi-class classification approach using XLNet showed promising performance, including in detecting the subtle differences between pain-related and recreation-related misuse. The distinct clustering of class-specific keywords may help conduct targeted data collection, overcoming under-representation of minority classes.

CONCLUSION

Machine learning can help identify pain-related and recreational-related OM contents on Twitter to potentially enable the study of the characteristics of individuals exhibiting such behavior.

Unraveling the energetic significance of chemical events in enzyme catalysis via machine-learning based regression approach

The bacterial enzyme class of β-lactamases are involved in benzylpenicillin acylation reactions, which are currently being revisited using hybrid quantum mechanical molecular mechanical (QM/MM) chain-of-states pathway optimizations. Minimum energy pathways are sampled by reoptimizing pathway geometry under different representative protein environments obtained through constrained molecular dynamics simulations. Predictive potential energy surface models in the reaction space are trained with machine-learning regression techniques. Herein, using TEM-1/benzylpenicillin acylation reaction as the model system, we introduce two model-independent criteria for delineating the energetic contributions and correlations in the predicted reaction space. Both methods are demonstrated to effectively quantify the energetic contribution of each chemical process and identify the rate limiting step of enzymatic reaction with high degrees of freedom. The consistency of the current workflow is tested under seven levels of quantum chemistry theory and three non-linear machine-learning regression models. The proposed approaches are validated to provide qualitative compliance with experimental mutagenesis studies.

Sensitivity analysis of FDG PET tumor voxel cluster radiomics and dosimetry for predicting mid-chemoradiation regional response of locally advanced lung cancer

We investigated the sensitivity of regional tumor response prediction to variability in voxel clustering techniques, imaging features, and machine learning algorithms in 25 patients with locally advanced non-small cell lung cancer (LA-NSCLC) enrolled on the FLARE-RT clinical trial. Metabolic tumor volumes (MTV) from pre-chemoradiation (PETpre) and mid-chemoradiation fluorodeoxyglucose-positron emission tomography (FDG PET) images (PETmid) were subdivided into K-means or hierarchical voxel clusters by standardized uptake values (SUV) and 3D-positions. MTV cluster separability was evaluated by CH index, and morphologic changes were captured by Dice similarity and centroid Euclidean distance. PETpre conventional features included SUVmean, MTV/MTV cluster size, and mean radiation dose. PETpre radiomics consisted of 41 intensity histogram and 3D texture features (PET Oncology Radiomics Test Suite) extracted from MTV or MTV clusters. Machine learning models (multiple linear regression, support vector regression, logistic regression, support vector machines) of conventional features or radiomic features were constructed to predict PETmid response. Leave-one-out-cross-validated root-mean-squared-error (RMSE) for continuous response regression (ΔSUVmean) and area-under-receiver-operating-characteristic-curve (AUC) for binary response classification were calculated. K-means MTV 2-clusters (MTVhi, MTVlo) achieved maximum CH index separability (Friedman p < 0.001). Between PETpre and PETmid, MTV cluster pairs overlapped (Dice 0.70-0.87) and migrated 0.6-1.1 cm. PETmid ΔSUVmean response prediction was superior in MTV and MTVlo (RMSE = 0.17-0.21) compared to MTVhi (RMSE = 0.42-0.52, Friedman p < 0.001). PETmid ΔSUVmean response class prediction performance trended higher in MTVlo (AUC = 0.83-0.88) compared to MTVhi (AUC = 0.44-0.58, Friedman p = 0.052). Models were more sensitive to MTV/MTV cluster regions (Friedman p = 0.026) than feature sets/algorithms (Wilcoxon signed-rank p = 0.36). Top-ranked radiomic features included GLZSM-LZHGE (large-zone-high-SUV), GTSDM-CP (cluster-prominence), GTSDM-CS (cluster-shade) and NGTDM-CNT (contrast). Top-ranked features were consistent between MTVhi and MTVlo cluster pairs but varied between MTVhi-MTVlo clusters, reflecting distinct regional radiomic phenotypes. Variability in tumor voxel cluster response prediction can inform robust radiomic target definition for risk-adaptive chemoradiation in patients with LA-NSCLC. FLARE-RT trial: NCT02773238.

Deciphering the Allosteric Process of the Phaeodactylum tricornutum Aureochrome 1a LOV Domain

The conformational-driven allosteric protein diatom Phaeodactylum tricornutum aureochrome 1a (PtAu1a) differs from other light-oxygen-voltage (LOV) proteins for its uncommon structural topology. The mechanism of signaling transduction in the PtAu1a LOV domain (AuLOV) including flanking helices remains unclear because of this dissimilarity, which hinders the study of PtAu1a as an optogenetic tool. To clarify this mechanism, we employed a combination of tree-based machine learning models, Markov state models, machine-learning-based community analysis, and transition path theory to quantitatively analyze the allosteric process. Our results are in good agreement with the reported experimental findings and reveal a previously overlooked Cα helix and protein linkers as important in promoting the protein conformational changes. This integrated approach can be considered as a general workflow and applied on other allosteric proteins to provide detailed information about their allosteric mechanisms.

ivis Dimensionality Reduction Framework for Biomacromolecular Simulations

Molecular dynamics (MD) simulations have been widely applied to study macromolecules including proteins. However, the high dimensionality of the data sets produced by simulations makes thorough analysis difficult and further hinders a deeper understanding of biomacromolecules. To gain more insights into the protein structure-function relations, appropriate dimensionality reduction methods are needed to project simulations onto low-dimensional spaces. Linear dimensionality reduction methods, such as principal component analysis (PCA) and time-structure-based independent component analysis (t-ICA), could not preserve sufficient structural information. Though better than linear methods, nonlinear methods, such as t-distributed stochastic neighbor embedding (t-SNE), still suffer from the limitations in avoiding system noise and keeping inter-cluster relations. ivis is a novel deep learning-based dimensionality reduction method originally developed for single-cell data sets. Here, we applied this framework for the study of light, oxygen, and voltage (LOV) domains of diatom Phaeodactylum tricornutum aureochrome 1a (PtAu1a). Compared with other methods, ivis is shown to be superior in constructing a Markov state model (MSM), preserving information of both local and global distances, and maintaining similarity between high and low dimensions with the least information loss. Moreover, the ivis framework is capable of providing new perspectives for deciphering residue-level protein allostery through the feature weights in the neural network. Overall, ivis is a promising member of the analysis toolbox for proteins.

Allosteric Regulation at the Crossroads of New Technologies: Multiscale Modeling, Networks, and Machine Learning

Allosteric regulation is a common mechanism employed by complex biomolecular systems for regulation of activity and adaptability in the cellular environment, serving as an effective molecular tool for cellular communication. As an intrinsic but elusive property, allostery is a ubiquitous phenomenon where binding or disturbing of a distal site in a protein can functionally control its activity and is considered as the "second secret of life." The fundamental biological importance and complexity of these processes require a multi-faceted platform of synergistically integrated approaches for prediction and characterization of allosteric functional states, atomistic reconstruction of allosteric regulatory mechanisms and discovery of allosteric modulators. The unifying theme and overarching goal of allosteric regulation studies in recent years have been integration between emerging experiment and computational approaches and technologies to advance quantitative characterization of allosteric mechanisms in proteins. Despite significant advances, the quantitative characterization and reliable prediction of functional allosteric states, interactions, and mechanisms continue to present highly challenging problems in the field. In this review, we discuss simulation-based multiscale approaches, experiment-informed Markovian models, and network modeling of allostery and information-theoretical approaches that can describe the thermodynamics and hierarchy allosteric states and the molecular basis of allosteric mechanisms. The wealth of structural and functional information along with diversity and complexity of allosteric mechanisms in therapeutically important protein families have provided a well-suited platform for development of data-driven research strategies. Data-centric integration of chemistry, biology and computer science using artificial intelligence technologies has gained a significant momentum and at the forefront of many cross-disciplinary efforts. We discuss new developments in the machine learning field and the emergence of deep learning and deep reinforcement learning applications in modeling of molecular mechanisms and allosteric proteins. The experiment-guided integrated approaches empowered by recent advances in multiscale modeling, network science, and machine learning can lead to more reliable prediction of allosteric regulatory mechanisms and discovery of allosteric modulators for therapeutically important protein targets.