Molecular Informatics最新文献

The multi-step retrosynthesis problem can be solved by a search algorithm, such as Monte Carlo tree search (MCTS). The performance of multistep retrosynthesis, as measured by a trade-off in search time and route solvability, therefore depends on the hyperparameters of the search algorithm. In this paper, we demonstrated the effect of three MCTS hyperparameters (number of iterations, tree depth, and tree width) on metrics such as Linear integrated speed-accuracy score (LISAS) and Inverse efficiency score which consider both route solvability and search time. This exploration was conducted by employing three data-driven approaches, namely a systematic grid search, Bayesian optimization over an ensemble of molecules to obtain static MCTS hyperparameters, and a machine learning approach to dynamically predict optimal MCTS hyperparameters given an input target molecule. With the obtained results on the internal dataset, we demonstrated that it is possible to identify a hyperparameter set which outperforms the current AiZynthFinder default setting. It appeared optimal across a variety of target input molecules, both on proprietary and public datasets. The settings identified with the in-house dataset reached a solvability of 93 % and median search time of 151 s for the in-house dataset, and a 74 % solvability and 114 s for the ChEMBL dataset. These numbers can be compared to the current default settings which solved 85 % and 73 % during a median time of 110s and 84 s, for in-house and ChEMBL, respectively.

Phenothiazine derivatives can unselectively inhibit the trypanothione-dependent antioxidant system enzyme trypanothione reductase (TR). A virtual screening of 2163 phenothiazine derivatives from the ZINC15 and PubChem databases docked on the active site of T. cruzi TR showed that 285 compounds have higher affinity than the natural ligand trypanothione disulfide. 244 compounds showed higher affinity toward the parasite's enzyme than to its human homolog glutathione reductase. Protein-ligand interaction profiling predicted that the main interactions for the top scored compounds were with residues important for trypanothione disulfide binding: Phe396, Pro398, Leu399, His461, Glu466, and Glu467, particularly His461, which participates in catalysis. Two compounds with the desired profiles, ZINC1033681 (Zn_C687) and ZINC10213096 (Zn_C216), decreased parasite growth by 20 % and 50 %, respectively. They behaved as mixed-type inhibitors of recombinant TR, with Ki values of 59 and 47 μM, respectively. This study provides a further understanding of the potential of phenothiazine derivatives as TR inhibitors.

Conjugated QSPR models for reactions integrate fundamental chemical laws expressed by mathematical equations with machine learning algorithms. Herein we present a methodology for building conjugated QSPR models integrated with the Arrhenius equation. Conjugated QSPR models were used to predict kinetic characteristics of cycloaddition reactions related by the Arrhenius equation: rate constant $logk$ , pre-exponential factor $logA$ , and activation energy $E_a$ . They were benchmarked against single-task (individual and equation-based models) and multi-task models. In individual models, all characteristics were modeled separately, while in multi-task models $logk$ , $logA$ and $E_a$ were treated cooperatively. An equation-based model assessed $logk$ using the Arrhenius equation and $logA$ and $E_a$ values predicted by individual models. It has been demonstrated that the conjugated QSPR models can accurately predict the reaction rate constants at extreme temperatures, at which reaction rate constants hardly can be measured experimentally. Also, in the case of small training sets conjugated models are more robust than related single-task approaches.

SARS-CoV-2 gained crucial mutations at the receptor binding domain (RBD) that often changed the course of the pandemic leading to new waves with increased case fatality. Variants are observed with enhanced transmission and immune invasion abilities. Thus, predicting future variants with enhanced transmission ability is a problem of utmost research interest. Here, we have developed a multi-tier exhaustive SARS-CoV-2 mutation screening platform combining MM/GBSA, extensive molecular dynamics simulations, and steered molecular dynamics to identify RBD mutants with enhanced ACE2 binding capability. We have identified four RBM mutations (F490K, S494K, G504F, and the P499L) with significantly higher ACE2 binding abilities than wild-type RBD. Compared to wild-type RBD, they all form stable complexes with more hydrogen bonds and salt-bridge interactions with ACE2. Our simulation data suggest that these mutations allosterically alter the packing of the RBM interface of the RBD-ACE2 complex. As a result, the rupture force required to break the RBD-ACE2 contacts is significantly higher for these mutants.

Dimensionality reduction (DR) techniques are used for various purposes such as exploratory data analysis. A commonly employed linear DR technique is principal component analysis (PCA), which is one of the most popular methods for DR. Owing to its linear nature, PCA enables the determination of axes in a low-dimensional space and the calculation of corresponding loading vectors. However, PCA cannot necessarily extract important features of non-linearly distributed data. This study presents a technique aimed at aiding the interpretation of data reduced through non-linear DR methods. In the proposed method, non-linear dimensionally reduced data was clustered via a density-based clustering method. Thereafter, the obtained cluster labels were classified by random forest (RF) classifiers. Further, feature importance (FI) of RF classifiers and Spearman's rank correlation coefficients between predictive probabilities to obtained clusters and original feature values were utilized for characterizing the visualized dimensionally reduced data. The results revealed that the proposed method can provide the interpretable FI-based images of the handwritten digits dataset. Moreover, the proposed method was also applied to the polymer dataset. The study found that incorporating signed FI was advantageous in achieving a meaningful interpretation. Furthermore, Gaussian process regression was utilized to produce intuitive FI-based heatmaps on a 2-dimensional space for greater ease of understanding. Additionally, to enhance the interpretability of the obtained clusters, a feature selection technique called Boruta was applied. The Boruta feature selection method worked effectively to interpret the obtained clusters with limited and commonly important features. Additionally, the study suggested that computing FI solely from substructure-based descriptors could further enhance the interpretability of the results. Finally, the automation of the proposed method was investigated, and through maximizing the target score based on the quality of both the DR and clustering, indicative results were automatically obtained for both the handwritten digits and polymer datasets.

The new high selective mAChRs M3 inhibitors with IC₅₀ in nanomolecular ranges, which can be the prototypes for effective COPD and asthma treatment drugs, were discovered with computational approaches among trifluoromethyl containing hexahydropyrimidinones/thiones. Compounds [6-(4-ethoxy-3-methoxy-phenyl)-4-hydroxy-2-thioxo-4-(trifluoromethyl)hexahydropyrimidin-5-yl]-phenyl-methanone (THPT-1) and 5-benzoyl-6-(3,4-dimethoxyphenyl)-4-hydroxy-4-(trifluoromethyl)hexahydropyrimidin-2-one (THPO-4) have been proved to be a highly effective (with IC₅₀ values of 1.62 ⋅ 10^-7  M and 3.09 ⋅ 10^-9  M, respectively) at the same concentrations significantly competitive inhibit the signal conduction through mAChR3 in comparison with ipratropium bromide, without significant effect on mAChR2, nicotinic cholinergic and adrenergic receptors.

Androgen receptor (AR) inhibition remains the primary strategy to combat the progression of prostate cancer (PC). However, all clinically used AR inhibitors target the ligand-binding domain (LBD), which is highly susceptible to truncations through splicing or mutations that confer drug resistance. Thus, there exists an urgent need for AR inhibitors with novel modes of action. We thus launched a virtual screening of an ultra-large chemical library to find novel inhibitors of the AR DNA-binding domain (DBD) at two sites: protein-DNA interface (P-box) and dimerization site (D-box). The compounds selected through vigorous computational filtering were then experimentally validated. We identified several novel chemotypes that effectively suppress transcriptional activity of AR and its splice variant V7. The identified compounds represent previously unexplored chemical scaffolds with a mechanism of action that evades the conventional drug resistance manifested through LBD mutations. Additionally, we describe the binding features required to inhibit AR DBD at both P-box and D-box target sites.

In this study we present deimos, a computational methodology for optimal grouping, applied on the read-across prediction of engineered nanomaterials' (ENMs) toxicity-related properties. The method is based on the formulation and the solution of a mixed-integer optimization program (MILP) problem that automatically and simultaneously performs feature selection, defines the grouping boundaries according to the response variable and develops linear regression models in each group. For each group/region, the characteristic centroid is defined in order to allocate untested ENMs to the groups. The deimos MILP problem is integrated in a broader optimization workflow that selects the best performing methodology between the standard multiple linear regression (MLR), the least absolute shrinkage and selection operator (LASSO) models and the proposed deimos multiple-region model. The performance of the suggested methodology is demonstrated through the application to benchmark ENMs datasets and comparison with other predictive modelling approaches. However, the proposed method can be applied to property prediction of other than ENM chemical entities and it is not limited to ENMs toxicity prediction.

Computational de novo drug design is a challenging issue in medicine, and it is desirable to consider all of the relevant information of the biological systems in a disease state. Here, we propose a novel computational method to generate drug candidate molecular structures from patient gene expression profiles via deep learning, which we call DRAGONET. Our model can generate new molecules that are likely to counteract disease-specific gene expression patterns in patients, which is made possible by exploring the latent space constructed by a transformer-based variational autoencoder and integrating the substructures of disease-correlated molecules. We applied DRAGONET to generate drug candidate molecules for gastric cancer, atopic dermatitis, and Alzheimer's disease, and demonstrated that the newly generated molecules were chemically similar to registered drugs for each disease. This approach is applicable to diseases with unknown therapeutic target proteins and will make a significant contribution to the field of precision medicine.