Simplified Molecular Input Line Entry System Research Articles

Cheminformatic Identification of Tyrosyl-DNA Phosphodiesterase 1 (Tdp1) Inhibitors: A Comparative Study of SMILES-Based Supervised Machine Learning Models

Background: Tyrosyl-DNA phosphodiesterase 1 (Tdp1) repairs damages in DNA induced by abortive topoisomerase 1 activity; however, maintenance of genetic integrity may sustain cellular division of neoplastic cells. It follows that Tdp1-targeting chemical inhibitors could synergize well with existing chemotherapy drugs to deny cancer growth; therefore, identification of Tdp1 inhibitors may advance precision medicine in oncology. Objective: Current computational research efforts focus primarily on molecular docking simulations, though datasets involving three-dimensional molecular structures are often hard to curate and computationally expensive to store and process. We propose the use of simplified molecular input line entry system (SMILES) chemical representations to train supervised machine learning (ML) models, aiming to predict potential Tdp1 inhibitors. Methods: An open-sourced consensus dataset containing the inhibitory activity of numerous chemicals against Tdp1 was obtained from Kaggle. Various ML algorithms were trained, ranging from simple algorithms to ensemble methods and deep neural networks. For algorithms requiring numerical data, SMILES were converted to chemical descriptors using RDKit, an open-sourced Python cheminformatics library. Results: Out of 13 optimized ML models with rigorously tuned hyperparameters, the random forest model gave the best results, yielding a receiver operating characteristics-area under curve of 0.7421, testing accuracy of 0.6815, sensitivity of 0.6444, specificity of 0.7156, precision of 0.6753, and F1 score of 0.6595. Conclusions: Ensemble methods, especially the bootstrap aggregation mechanism adopted by random forest, outperformed other ML algorithms in classifying Tdp1 inhibitors from non-inhibitors using SMILES. The discovery of Tdp1 inhibitors could unlock more treatment regimens for cancer patients, allowing for therapies tailored to the patient’s condition.

Prediction of Glass Transition Temperature of Polymers Using Simple Machine Learning.

Polymer materials have garnered significant attention due to their exceptional mechanical properties and diverse industrial applications. Understanding the glass transition temperature (Tg) of polymers is critical to prevent operational failures at specific temperatures. Traditional methods for measuring Tg, such as differential scanning calorimetry (DSC) and dynamic mechanical analysis, while accurate, are often time-consuming, costly, and susceptible to inaccuracies due to random and uncertain factors. To address these limitations, the aim of the present study is to investigate the potential of Simplified Molecular Input Line Entry System (SMILES) as descriptors in simple machine learning models to predict Tg efficiently and reliably. Five models were utilized: k-nearest neighbors (KNNs), support vector regression (SVR), extreme gradient boosting (XGBoost), artificial neural network (ANN), and recurrent neural network (RNN). SMILES descriptors were converted into numerical data using either One Hot Encoding (OHE) or Natural Language Processing (NLP). The study found that SMILES inputs with fewer than 200 characters were inadequate for accurately describing compound structures, while inputs exceeding 200 characters diminished model performance due to the curse of dimensionality. The ANN model achieved the highest R2 value of 0.79; however, the XGB model, with an R2 value of 0.774, exhibited the highest stability and shorter training times compared to other models, making it the preferred choice for Tg prediction. The efficiency of the OHE method over NLP was demonstrated by faster training times across the KNN, SVR, XGB, and ANN models. Validation of new polymer data showed the XGB model's robustness, with an average prediction deviation of 9.76 from actual Tg values. These findings underscore the importance of optimizing SMILES conversion methods and model parameters to enhance prediction reliability. Future research should focus on improving model accuracy and generalizability by incorporating additional features and advanced techniques. This study contributes to the development of efficient and reliable predictive models for polymer properties, facilitating the design and application of new polymer materials.

Unlocking comprehensive molecular design across all scenarios with large language model and unordered chemical language.

Molecular generation stands at the forefront of AI-driven technologies, playing a crucial role in accelerating the development of small molecule drugs. The intricate nature of practical drug discovery necessitates the development of a versatile molecular generation framework that can tackle diverse drug design challenges. However, existing methodologies often struggle to encompass all aspects of small molecule drug design, particularly those rooted in language models, especially in tasks like linker design, due to the autoregressive nature of large language model-based approaches. To empower a language model for a wider range of molecular design tasks, we introduce an unordered simplified molecular-input line-entry system based on fragments (FU-SMILES). Building upon this foundation, we propose FragGPT, a universal fragment-based molecular generation model. Initially pretrained on extensive molecular datasets, FragGPT utilizes FU-SMILES to facilitate efficient generation across various practical applications, such as de novo molecule design, linker design, R-group exploration, scaffold hopping, and side chain optimization. Furthermore, we integrate conditional generation and reinforcement learning (RL) methodologies to ensure that the generated molecules possess multiple desired biological and physicochemical properties. Experimental results across diverse scenarios validate FragGPT's superiority in generating molecules with enhanced properties and novel structures, outperforming existing state-of-the-art models. Moreover, its robust drug design capability is further corroborated through real-world drug design cases.

Molecular interaction analysis of ferulic acid (4-hydroxy-3-methoxycinnamic acid) as main bioactive compound from palm oil waste against MCF-7 receptors: An in silico study.

Ferulic acid (4-hydroxy-3-methoxycinnamic acid) is a phytochemical compound that is commonly found in conjugated forms within mono-, di-, polysaccharides and other organic compounds in cell walls of grain, fruits, and vegetables. This compound is highly abundant in the palm oil waste. The aim of the study was to predict the anticancer activity of ferulic acid against the breast cancer cell lines (MCF-7) receptors through a computational analysis. MCF-7 receptors with PDB IDs of 1R5K, 2IOG, 4IV2, 4IW6, 5DUE, 5T92, and 5U2B were selected based on the Simplified Molecular Input Line Entry System (SMILES) similarity of the native ligand. Thereafter, the protein was prepared on Chimera 1.16 and docked with ferulic acid on Autodock Vina 1.2.5. The ligand-protein complex interaction was validated by computing the root mean square fluctuation (RMSF) and radius of gyration (Rg) through molecular dynamic simulation. In addition, an absorption, distribution, metabolism, excretion, and toxicity (ADMET) prediction was performed on ferulic acid using the pkCSM platform. The molecular docking revealed that the ferulic acid could interact with all receptors as indicated by the affinity energy <-5 kcal/mol. The compound had the most optimum interaction with receptor 2IOG (affinity energy=-6.96 kcal/mol), involving hydrophobic interaction (n=12) and polar hydrogen interaction (n=4). The molecular dynamic simulation revealed that the complex had an RMSF of 1.713 Å with a fluctuation of Rg value around 1.000 Å. The ADMET properties of ferulic acid suggested that the compound is an ideal drug candidate. In conclusion, this study suggested that ferulic acid, which can be isolated from palm oil waste, has the potential to interact with MCF-7 receptors.

The Block Copolymer Phase Behavior Database.

The Block Copolymer Database (BCDB) is a platform that allows users to search, submit, visualize, benchmark, and download experimental phase measurements and their associated characterization information for di- and multiblock copolymers. To the best of our knowledge, there is no widely accepted data model for publishing experimental and simulation data on block copolymer self-assembly. This proposed data schema with traceable information can accommodate any number of blocks and at the time of publication contains over 5400 block copolymer total melt phase measurements mined from the literature and manually curated and simulation data points of the phase diagram generated from self-consistent field theory that can rapidly be augmented. This database can be accessed via the Community Resource for Innovation in Polymer Technology (CRIPT) web application and the Materials Data Facility. The chemical structure of the polymer is encoded in BigSMILES, an extension of the Simplified Molecular-Input Line-Entry System (SMILES) into the macromolecular domain, and the user can search repeat units and functional groups using the SMARTS search syntax (SMILES Arbitrary Target Specification). The user can also query characterization and phase information using Structured Query Language (SQL) and download custom sets of block copolymer data to train machine learning models. Finally, a protocol is presented in which GPT-4, an AI-powered large language model, can be used to rapidly screen and identify block copolymer papers from the literature using only the abstract text and determine whether they have BCDB data, allowing the database to grow as the number of published papers on the World Wide Web increases. The F1 score for this model is 0.74. This platform is an important step in making polymer data more accessible to the broader community.

Comparison of Machine Learning Approaches for Prediction of the Equivalent Alkane Carbon Number for Microemulsions Based on Molecular Properties.

The chemical properties of oils are vital in the design of microemulsion systems. The hydrophilic-lipophilic difference equation used to predict microemulsions' phase behavior expresses the oils' physiochemical properties as the equivalent alkane carbon number (EACN). The experimental determination of EACN requires knowledge of the temperature dependence of the microemulsion system and the effects of different surfactant concentrations. Thus, the experimental determination is time-intensive and tedious, requiring days to months for proper separations. Furthermore, the experiments require high purity of chemicals because microemulsions are sensitive to impurities. Our work focuses on the quick and reliable predictions of the EACN with machine learning (ML) models. Due to the immaturity of ML chemical predictions, we compare three graph neural networks (GNNs) and a gradient-boosted tree algorithm, known as XGBoost. The GNNs use the molecular structures represented as simplified molecular-input line-entry system (SMILES) codes for the initial input, which allows us to assess whether geometry optimization is necessary for reliable results. The XGBoost model also begins with the SMILES representations of the molecules but uses molecular descriptors instead of geometry optimizations. The best model tested (crystal graph convolutional neural network with Merck molecular force field-94) has an error of 1.15 EACN units of the true EACN for unknown data with the errors skewed toward zero and an R2 score of 0.9.

Toxicity Prediction for Immune Thrombocytopenia Caused by Drugs Based on Logistic Regression with Feature Importance

Background: One of the problems in drug discovery that can be solved by artificial intelligence is toxicity prediction. In drug-induced immune thrombocytopenia, toxicity can arise in patients after five to ten days by significant bleeding caused by drugdependent antibodies. In clinical trials, when this condition occurs, all the drugs consumed by patients should be stopped, although sometimes this is not possible, especially for older patients who are dependent on their medication. Therefore, being able to predict toxicity in drug-induced immune thrombocytopenia is very important. Computational technologies, such as machine learning, can help predict toxicity better than empirical techniques owing to the lower cost and faster processing. Objective: Previous studies used the KNN method. However, the performance of these approaches needs to be enhanced. This study proposes a Logistic Regression to improve accuracy scores. Methods: In this study, we present a new model for drug-induced immune thrombocytopenia using a machine learning method. Our model extracts several features from the Simplified Molecular Input Line Entry System (SMILES). These features were fused and cleaned, and the important features were selected using the SelectKBest method. The model uses a Logistic Regression that is optimized and tuned by the Grid Search Cross Validation. Results: The highest accuracy occurred when using features from PADEL, CDK, RDKIT, MORDRED, BLUEDESC combinations, resulting in an accuracy of 80%. Conclusion: Our proposed model outperforms previous studies in accuracy categories. The information and source code is accessible online at Github: https://github.com/Osphanie/Thrombocytopenia

Multicellular target QSAR models for predicting of novel inhibitors against pancreatic cancer by Monte Carlo approach

Pancreatic cancer is one of the most fatal malignant tumors, with a poor prognosis. Nowadays, computational methods for designing and predicting anticancer agents have gained increased attention. In this work, a QSAR study is conducted to further investigate the activity of 5797 agents against pancreatic cancer cell lines. QSAR models are a system of classification approaches that are used to construct and validate the predictive potential of created models. This approach demonstrates the prediction in the active (1) and inactive (−1) format. Here, the classification QSAR technique is employed to construct and evaluate models. The classification QSAR model has been constructed by the Monte Carlo method based on the representation of the molecular structure by SMILES (simplified molecular input line entry systems) using the CORAL software. Four random splits of data by making four sets such as active training, passive training, calibration, and validation sets are assessed. The numerical values of sensitivity, specificity, and accuracy for the validation set of splits 1 to 4 are in the range of 0.7726–0.7901, 0.9007–0.9220, and 0.8548–0.8730 respectively. The Matthews correlation coefficient (MCC) of these models for external validation sets is found to be in the range of 0.6832 to 0.7215.

Predicting blood–brain barrier permeability of molecules with a large language model and machine learning

Predicting the blood–brain barrier (BBB) permeability of small-molecule compounds using a novel artificial intelligence platform is necessary for drug discovery. Machine learning and a large language model on artificial intelligence (AI) tools improve the accuracy and shorten the time for new drug development. The primary goal of this research is to develop artificial intelligence (AI) computing models and novel deep learning architectures capable of predicting whether molecules can permeate the human blood–brain barrier (BBB). The in silico (computational) and in vitro (experimental) results were validated by the Natural Products Research Laboratories (NPRL) at China Medical University Hospital (CMUH). The transformer-based MegaMolBART was used as the simplified molecular input line entry system (SMILES) encoder with an XGBoost classifier as an in silico method to check if a molecule could cross through the BBB. We used Morgan or Circular fingerprints to apply the Morgan algorithm to a set of atomic invariants as a baseline encoder also with an XGBoost classifier to compare the results. BBB permeability was assessed in vitro using three-dimensional (3D) human BBB spheroids (human brain microvascular endothelial cells, brain vascular pericytes, and astrocytes). Using multiple BBB databases, the results of the final in silico transformer and XGBoost model achieved an area under the receiver operating characteristic curve of 0.88 on the held-out test dataset. Temozolomide (TMZ) and 21 randomly selected BBB permeable compounds (Pred scores = 1, indicating BBB-permeable) from the NPRL penetrated human BBB spheroid cells. No evidence suggests that ferulic acid or five BBB-impermeable compounds (Pred scores < 1.29423E−05, which designate compounds that pass through the human BBB) can pass through the spheroid cells of the BBB. Our validation of in vitro experiments indicated that the in silico prediction of small-molecule permeation in the BBB model is accurate. Transformer-based models like MegaMolBART, leveraging the SMILES representations of molecules, show great promise for applications in new drug discovery. These models have the potential to accelerate the development of novel targeted treatments for disorders of the central nervous system.

Advancements in hand-drawn chemical structure recognition through an enhanced DECIMER architecture

Accurate recognition of hand-drawn chemical structures is crucial for digitising hand-written chemical information in traditional laboratory notebooks or facilitating stylus-based structure entry on tablets or smartphones. However, the inherent variability in hand-drawn structures poses challenges for existing Optical Chemical Structure Recognition (OCSR) software. To address this, we present an enhanced Deep lEarning for Chemical ImagE Recognition (DECIMER) architecture that leverages a combination of Convolutional Neural Networks (CNNs) and Transformers to improve the recognition of hand-drawn chemical structures. The model incorporates an EfficientNetV2 CNN encoder that extracts features from hand-drawn images, followed by a Transformer decoder that converts the extracted features into Simplified Molecular Input Line Entry System (SMILES) strings. Our models were trained using synthetic hand-drawn images generated by RanDepict, a tool for depicting chemical structures with different style elements. A benchmark was performed using a real-world dataset of hand-drawn chemical structures to evaluate the model's performance. The results indicate that our improved DECIMER architecture exhibits a significantly enhanced recognition accuracy compared to other approaches.Scientific contributionThe new DECIMER model presented here refines our previous research efforts and is currently the only open-source model tailored specifically for the recognition of hand-drawn chemical structures. The enhanced model performs better in handling variations in handwriting styles, line thicknesses, and background noise, making it suitable for real-world applications. The DECIMER hand-drawn structure recognition model and its source code have been made available as an open-source package under a permissive license.Graphical

Identifying High-Quality Leads among Screened Anticancerous Compounds Using SMILES Representations.

Cancer is a lethal disease that affects numerous people worldwide. Chemotherapy stands as one of the most effective treatment regimens to combat cancer. Nevertheless, anticancer drugs face a high failure rate due to safety and efficacy issues. Drug failure could be subdued by instigating drug leads with reduced toxicity and enhanced efficacy. Computer-aided drug discovery endorses drug leads in manoeuvring protein and ligand structures or representations. Simplified molecular input line entry system (SMILES) is a linear notation representing the three-dimensional structure of a molecule using symbols and alphanumeric characters. SMILES representation hoards rings and scaffold structures in its depiction. Mining ring and scaffold patterns from molecular SMILES would assist in ascertaining biological properties based on molecular patterns. Moreover, the emergence of artificial intelligence (AI) technologies would accelerate identification of efficient anticancer drug leads. AI algorithms proclaimed for their pattern recognition ability could be employed for identifying molecular patterns from SMILES representation, thereby enabling property prediction. Consequently, we developed a multilayer perceptron (MLP) model for the prediction of anticancer activity using SMILES of NCI-60 cancer growth inhibition data. Furthermore, the top 8 frequent scaffolds were identified on preliminary analysis of cancer growth inhibition data and ChEMBL drugs. The developed MLP model classified anticancer and nonanticancer compounds with a classification accuracy of 0.92. Also, benchmarking of the developed model with machine learning algorithms exhibited better performance of the MLP model.

Can large language models understand molecules?

PurposeLarge Language Models (LLMs) like Generative Pre-trained Transformer (GPT) from OpenAI and LLaMA (Large Language Model Meta AI) from Meta AI are increasingly recognized for their potential in the field of cheminformatics, particularly in understanding Simplified Molecular Input Line Entry System (SMILES), a standard method for representing chemical structures. These LLMs also have the ability to decode SMILES strings into vector representations.MethodWe investigate the performance of GPT and LLaMA compared to pre-trained models on SMILES in embedding SMILES strings on downstream tasks, focusing on two key applications: molecular property prediction and drug-drug interaction prediction.ResultsWe find that SMILES embeddings generated using LLaMA outperform those from GPT in both molecular property and DDI prediction tasks. Notably, LLaMA-based SMILES embeddings show results comparable to pre-trained models on SMILES in molecular prediction tasks and outperform the pre-trained models for the DDI prediction tasks.ConclusionThe performance of LLMs in generating SMILES embeddings shows great potential for further investigation of these models for molecular embedding. We hope our study bridges the gap between LLMs and molecular embedding, motivating additional research into the potential of LLMs in the molecular representation field. GitHub: https://github.com/sshaghayeghs/LLaMA-VS-GPT.

ImageDTA: A Simple Model for Drug-Target Binding Affinity Prediction.

Predicting the drug-target binding affinity (DTA) is crucial in drug discovery, and an increasing number of researchers are using artificial intelligence techniques to make such predictions. Many effective deep neural network prediction models have been proposed. However, current methods need improvement in accuracy, complexity, and efficiency. In this study, we propose a method based on a multiscale 2-dimensional convolutional neural network (CNN), namely ImageDTA. Many studies have shown that CNN achieves good learning effects with limited data. Therefore, we take a unique perspective by treating the word vector encoded with a simplified molecular input line entry system (SMILES) string as an "image" and processing it like handling images, fully leveraging the efficient processing capabilities of CNN for image data. Furthermore, we show that ImageDTA has higher training and inference efficiency than pretrained large models and outperforms attention-based graph neural network models in accuracy and interpretability. We also use visualization techniques to select appropriate convolutional kernel sizes, thereby increasing the network's interpretability.

Molecular sharing and molecular-specific representations for multimodal molecular property prediction

Molecular property prediction plays a crucial role in drug discovery and development. However, traditional experimental measurements and Quantitative Structure-Activity Relationship (QSAR) models are often expensive, time-consuming, and data acquisition is challenging. To overcome these limitations and challenges, this study innovatively proposes a fusion molecular property prediction method called molecular property prediction model (MSSP) to address the non-uniqueness of Simplified Molecular Input Line Entry System (SMILES) string representation and the difficulty of capturing global information in molecular graphs. This method extracts multiple fingerprint features and utilizes graph neural network encoding to map different modalities of molecules into molecular sharing and molecular-specific representation spaces, achieving modal alignment and fusion of molecules by combining molecular invariance and representation specificity. To enhance the interpretability and visualization capabilities of the model, graph attention mechanisms are introduced, enabling the identification and inference of important chemical fragments within molecules. Experimental results on publicly available cell line phenotype and kinase activity datasets demonstrate that MSSP outperforms the current state-of-the-art methods in molecular property prediction. Additionally, MSSP exhibits strong competitiveness across nine benchmark molecular property prediction datasets. Furthermore, in the task of predicting SRC kinase data properties, this study successfully screens promising therapeutic compounds from compound libraries by validating the predictions of the MSSP model and combining them with traditional methods such as molecular docking and molecular dynamics simulations. Multiple potential Lyn inhibitors have been discovered through this approach. The application of MSSP model is helpful to discover new molecules with new drug properties or functions, accelerate the process of drug discovery, save time and resources, and provide important guidance for drug discovery.

Predicting ADMET Properties from Molecule SMILE: A Bottom-Up Approach Using Attention-Based Graph Neural Networks.

Understanding the pharmacokinetics, safety and efficacy of candidate drugs is crucial for their success. One key aspect is the characterization of absorption, distribution, metabolism, excretion and toxicity (ADMET) properties, which require early assessment in the drug discovery and development process. This study aims to present an innovative approach for predicting ADMET properties using attention-based graph neural networks (GNNs). The model utilizes a graph-based representation of molecules directly derived from Simplified Molecular Input Line Entry System (SMILE) notation. Information is processed sequentially, from substructures to the whole molecule, employing a bottom-up approach. The developed GNN is tested and compared with existing approaches using six benchmark datasets and by encompassing regression (lipophilicity and aqueous solubility) and classification (CYP2C9, CYP2C19, CYP2D6 and CYP3A4 inhibition) tasks. Results show the effectiveness of our model, which bypasses the computationally expensive retrieval and selection of molecular descriptors. This approach provides a valuable tool for high-throughput screening, facilitating early assessment of ADMET properties and enhancing the likelihood of drug success in the development pipeline.

Graphormer-IR: Graph Transformers Predict Experimental IR Spectra Using Highly Specialized Attention.

Infrared (IR) spectroscopy is an important analytical tool in various chemical and forensic domains and a great deal of effort has gone into developing in silico methods for predicting experimental spectra. A key challenge in this regard is generating highly accurate spectra quickly to enable real-time feedback between computation and experiment. Here, we employ Graphormer, a graph neural network (GNN) transformer, to predict IR spectra using only simplified molecular-input line-entry system (SMILES) strings. Our data set includes 53,528 high-quality spectra, measured in five different experimental media (i.e., phases), for molecules containing the elements H, C, N, O, F, Si, S, P, Cl, Br, and I. When using only atomic numbers for node encodings, Graphormer-IR achieved a mean test spectral information similarity (SISμ) value of 0.8449 ± 0.0012 (n = 5), which surpasses that the current state-of-the-art model Chemprop-IR (SISμ = 0.8409 ± 0.0014, n = 5) with only 36% of the encoded information. Augmenting node embeddings with additional node-level descriptors in learned embeddings generated through a multilayer perceptron improves scores to SISμ = 0.8523 ± 0.0006, a total improvement of 19.7σ (t = 19). These improved scores show how Graphormer-IR excels in capturing long-range interactions like hydrogen bonding, anharmonic peak positions in experimental spectra, and stretching frequencies of uncommon functional groups. Scaling our architecture to 210 attention heads demonstrates specialist-like behavior for distinct IR frequencies that improves model performance. Our model utilizes novel architectures, including a global node for phase encoding, learned node feature embeddings, and a one-dimensional (1D) smoothing convolutional neural network (CNN). Graphormer-IR's innovations underscore its value over traditional message-passing neural networks (MPNNs) due to its expressive embeddings and ability to capture long-range intramolecular relationships.

Cholinesterase Inhibitory and In Silico Toxicity Assessment of Thirty-Four Isoquinoline Alkaloids - Berberine as the Lead Compound.

Cholinesterase (ChE) inhibitors used currently in clinics for the treatment of Alzheimer's disease (AD) are the most prescribed drug class with nitrogen-containing chemical formula. Galanthamine, the latest generation anti-ChE drug, contains an isoquinoline structure. The aim of the current study was to investigate the inhibitory potential of thirty-four isoquinoline alkaloids, e.g. (-)-adlumidine, β-allocryptopine, berberine, (+)-bicuculline, (-)-bicuculline, (+)-bulbocapnine, (-)-canadine, (±)-chelidimerine, corydaldine, (±)-corydalidzine, (-)-corydalmine, (+)-cularicine, dehydrocavidine, (+)-fumariline, (-)-fumarophycine, (+)-α-hydrastine, (+)-isoboldine, 13-methylcolumbamine, (-)-norjuziphine, norsanguinarine, (-)-ophiocarpine, (-)-ophiocarpine-Noxide, oxocularine, oxosarcocapnine, palmatine, (+)-parfumine, protopine, (+)-reticuline, sanguinarine, (+)-scoulerine, (±)-sibiricine, (±)-sibiricine acetate, (-)-sinactine, and (-)-stylopine isolated from several Fumaria (fumitory) and Corydalis species towards acetyl- (AChE) and butyrylcholinesterase (BChE) by microtiter plate assays. The alkaloids with strong ChE inhibition were proceeded to molecular docking simulations as well as in silico toxicity screening for their mutagenic capacity through VEGA QSAR (AMES test) consensus model and VEGA platform as statistical approaches. The inputs were evaluated in a simplified molecular input-line entry system (SMILES). ChE inhibition assays indicated that the highest AChE inhibition was caused by berberine (IC50: 0.72 ± 0.04 μg/mL), palmatine (IC50: 6.29 ± 0.61 μg/mL), β-allocryptopine (IC50: 10.62 ± 0.45 μg/mL), (-)-sinactine (IC50: 11.94 ± 0.44 μg/mL), and dehydrocavidine (IC50: 15.01 ± 1.87 μg/mL) as compared to that of galanthamine (IC50: 0.74 ± 0.01 μg/mL), the reference drug with isoquinoline skeleton. Less number of the tested alkaloids exhibited notable BChE inhibition. Among them, berberine (IC50: 7.67 ± 0.36 μg/mL) and (-)-corydalmine (IC50: 7.78 ± 0.38 μg/mL) displayed a stronger inhibition than that of galanthamine (IC50: 12.02 ± 0.25 μg/mL). The mutagenic activity was shown for β-allocryptopine, (+)- and (-)-bicuculline, (±)-corydalidzine, (-)-corydalmine, (+)-cularicine, (-)- fumarophycine, (-)-norjuziphine, (-)-ophiocarpine-N-oxide, (+)-scoulerine, (-)-sinactine, and (-)- stylopine by means of in silico experiments. The results obtained by molecular docking simulations of berberine, palmatine, and (-)- corydalmine suggested that the estimated free ligand-binding energies of these compounds inside the binding domains of their targets are reasonable to make them capable of establishing strong polar and nonpolar bonds with the atoms of the active site amino acids. Our findings revealed that berberine, palmatin, and (-)-corydalmine stand out as the most promising isoquinoline alkaloids in terms of ChE inhibition. Among them, berberine has displayed a robust dual inhibition against both ChEs and could be evaluated further as a lead compound for AD.

MvMRL: a multi-view molecular representation learning method for molecular property prediction.

Effective molecular representation learning is very important for Artificial Intelligence-driven Drug Design because it affects the accuracy and efficiency of molecular property prediction and other molecular modeling relevant tasks. However, previous molecular representation learning studies often suffer from limitations, such as over-reliance on a single molecular representation, failure to fully capture both local and global information in molecular structure, and ineffective integration of multiscale features from different molecular representations. These limitations restrict the complete and accurate representation of molecular structure and properties, ultimately impacting the accuracy of predicting molecular properties. To this end, we propose a novel multi-view molecular representation learning method called MvMRL, which can incorporate feature information from multiple molecular representations and capture both local and global information from different views well, thus improving molecular property prediction. Specifically, MvMRL consists of four parts: a multiscale CNN-SE Simplified Molecular Input Line Entry System (SMILES) learning component and a multiscale Graph Neural Network encoder to extract local feature information and global feature information from the SMILES view and the molecular graph view, respectively; a Multi-Layer Perceptron network to capture complex non-linear relationship features from the molecular fingerprint view; and a dual cross-attention component to fuse feature information on the multi-views deeply for predicting molecular properties. We evaluate the performance of MvMRL on 11 benchmark datasets, and experimental results show that MvMRL outperforms state-of-the-art methods, indicating its rationality and effectiveness in molecular property prediction. The source code of MvMRL was released in https://github.com/jedison-github/MvMRL.

In-silico green toxicology approach toward discovering safer ligands for development of safe-by-design metal–organic frameworks

A vast variety of chemical compounds have been fabricated and commercialized, they not only result in industrial exposure during manufacturing and usage, but also have environmental impacts throughout their whole life cycle. Consequently, attempts to assess the risk of chemicals in terms of toxicology have never ceased. In-silico toxicology, also known as predictive toxicology, has advanced significantly over the last decade as a result of the drawbacks of experimental investigations. In this study, ProTox-III was applied to predict the toxicity of the ligands used for metal–organic framework (MOF) design and synthesis. Initially, 35 ligands, that have been frequently utilized for MOF synthesis and fabrication, were selected. Subsequently, canonical simplified molecular-input line-entry system (SMILES) of ligands were extracted from the PUBCHEM database and inserted into the ProTox-III online server. Ultimately, webserver outputs including LD50 and the probability of toxicological endpoints (cytotoxicity, carcinogenicity, mutagenicity, immunotoxicity, and ecotoxicity) were obtained and organized. According to retrieved LD50 data, the safest ligand was 5-hydroxyisophthalic. In contrast, the most hazardous ligand was 5-chlorobenzimidazole, with an LD50 of 8 mg/kg. Among evaluated endpoints, ecotoxicity was the most active and was detected in several imidazolate ligands. This data can open new horizons in design and development of green MOFs.

HBCVTr: an end-to-end transformer with a deep neural network hybrid model for anti-HBV and HCV activity predictor from SMILES

Hepatitis B and C viruses (HBV and HCV) are significant causes of chronic liver diseases, with approximately 350 million infections globally. To accelerate the finding of effective treatment options, we introduce HBCVTr, a novel ligand-based drug design (LBDD) method for predicting the inhibitory activity of small molecules against HBV and HCV. HBCVTr employs a hybrid model consisting of double encoders of transformers and a deep neural network to learn the relationship between small molecules’ simplified molecular-input line-entry system (SMILES) and their antiviral activity against HBV or HCV. The prediction accuracy of HBCVTr has surpassed baseline machine learning models and existing methods, with R-squared values of 0.641 and 0.721 for the HBV and HCV test sets, respectively. The trained models were successfully applied to virtual screening against 10 million compounds within 240 h, leading to the discovery of the top novel inhibitor candidates, including IJN04 for HBV and IJN12 and IJN19 for HCV. Molecular docking and dynamics simulations identified IJN04, IJN12, and IJN19 target proteins as the HBV core antigen, HCV NS5B RNA-dependent RNA polymerase, and HCV NS3/4A serine protease, respectively. Overall, HBCVTr offers a new and rapid drug discovery and development screening method targeting HBV and HCV.