Protein language models are performant in structure-free virtual screening.

Hitherto virtual screening (VS) has been typically performed using a structure-based drug design paradigm. Such methods typically require the use of molecular docking on high-resolution three-dimensional structures of a target protein-a computationally-intensive and time-consuming exercise. This work demonstrates that by employing protein language models and molecular graphs as inputs to a novel graph-to-transformer cross-attention mechanism, a screening power comparable to state-of-the-art structure-based models can be achieved. The implications thereof include highly expedited VS due to the greatly reduced compute required to run this model, and the ability to perform early stages of computer-aided drug design in the complete absence of 3D protein structures.

TransFoxMol: predicting molecular property with focused attention.

Predicting the biological properties of molecules is crucial in computer-aided drug development, yet it's often impeded by data scarcity and imbalance in many practical applications. Existing approaches are based on self-supervised learning or 3D data and using an increasing number of parameters to improve performance. These approaches may not take full advantage of established chemical knowledge and could inadvertently introduce noise into the respective model. In this study, we introduce a more elegant transformer-based framework with focused attention for molecular representation (TransFoxMol) to improve the understanding of artificial intelligence (AI) of molecular structure property relationships. TransFoxMol incorporates a multi-scale 2D molecular environment into a graph neural network + Transformer module and uses prior chemical maps to obtain a more focused attention landscape compared to that obtained using existing approaches. Experimental results show that TransFoxMol achieves state-of-the-art performance on MoleculeNet benchmarks and surpasses the performance of baselines that use self-supervised learning or geometry-enhanced strategies on small-scale datasets. Subsequent analyses indicate that TransFoxMol's predictions are highly interpretable and the clever use of chemical knowledge enables AI to perceive molecules in a simple but rational way, enhancing performance.

Identification of Helicobacter pylori-carcinogenic TNF-alpha-inducing protein inhibitors via daidzein derivatives through computational approaches.

Gastric cancer is considered a class 1 carcinogen that is closely linked to infection with Helicobacter pylori (H. pylori), which affects over 1 million people each year. However, the major challenge to fight against H. pylori and its associated gastric cancer due to drug resistance. This research gap had led our research team to investigate a potential drug candidate targeting the Helicobacter pylori-carcinogenic TNF-alpha-inducing protein. In this study, a total of 45 daidzein derivatives were investigated and the best 10 molecules were comprehensively investigated using in silico approaches for drug development, namely pass prediction, quantum calculations, molecular docking, molecular dynamics simulations, Lipinski rule evaluation, and prediction of pharmacokinetics. The molecular docking study was performed to evaluate the binding affinity between the target protein and the ligands. In addition, the stability of ligand-protein complexes was investigated by molecular dynamics simulations. Various parameters were analysed, including root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), radius of gyration (Rg), hydrogen bond analysis, principal component analysis (PCA) and dynamic cross-correlation matrix (DCCM). The results has confirmed that the ligand-protein complex CID: 129661094 (07) and 129664277 (08) formed stable interactions with the target protein. It was also found that CID: 129661094 (07) has greater hydrogen bond occupancy and stability, while the ligand-protein complex CID 129664277 (08) has greater conformational flexibility. Principal component analysis revealed that the ligand-protein complex CID: 129661094 (07) is more compact and stable. Hydrogen bond analysis revealed favourable interactions with the reported amino acid residues. Overall, this study suggests that daidzein derivatives in particular show promise as potential inhibitors of H. pylori.

Using Multiscale Molecular Modeling to Analyze Possible NS2b-NS3 Protease Inhibitors from Philippine Medicinal Plants.

Within the field of Philippine folkloric medicine, the utilization of indigenous plants like Euphorbia hirta (tawa-tawa), Carica papaya (papaya), and Psidium guajava (guava) as potential dengue remedies has gained attention. Yet, limited research exists on their comprehensive effects, particularly their anti-dengue activity. This study screened 2944 phytochemicals from various Philippine plants for anti-dengue activity. Absorption, distribution, metabolism, excretion, and toxicity (ADMET) profiling provided 1265 compounds demonstrating pharmacokinetic profiles suitable for human use. Molecular docking targeting the dengue virus NS2b-NS3 protease's catalytic triad (Asp 75, Ser 135, and His 51) identified ten ligands with higher docking scores than reference compounds idelalisib and nintedanib. Molecular dynamics simulations confirmed the stability of eight of these ligand-protease complexes. Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) analysis highlighted six ligands, including veramiline (-80.682 kJ/mol), cyclobranol (-70.943 kJ/mol), chlorogenin (-63.279 kJ/mol), 25beta-Hydroxyverazine (-61.951 kJ/mol), etiolin (-59.923 kJ/mol), and ecliptalbine (-56.932 kJ/mol) with favorable binding energies, high oral bioavailability, and drug-like properties. This integration of traditional medical knowledge with advanced computational drug discovery methods paves new pathways for the development of treatments for dengue.

In Silico Screening of 1,3,4-Thiadiazole Derivatives as Inhibitors of Vascular Endothelial Growth Factor Receptor-2 (VEGFR-2).

Angiogenesis plays a pivotal role in the growth, survival, and metastasis of solid tumors, with Vascular Endothelial Growth Factor Receptor-2 (VEGFR-2) being overexpressed in many human solid tumors, making it an appealing target for anti-cancer therapies. This study aimed to identify potential lead compounds with azole moiety exhibiting VEGFR-2 inhibitory effects. A ligand-based pharmacophore model was constructed using the X-ray crystallographic structure of VEGFR-2 complexed with tivozanib (PDB ID: 4ASE) to screen the ZINC15 database. Following virtual screening, six compounds demonstrated promising docking scores and drug-likeness comparable to tivozanib. These hits underwent detailed pharmacokinetic analysis to assess their absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties. Furthermore, Density Functional Theory (DFT) analysis was employed to investigate the molecular orbital properties of the top hits from molecular docking. Molecular dynamics (MD) simulations were conducted to evaluate the conformational stability of the complexes over a 100 ns run. Results indicated that the compounds (ZINC8914312, ZINC8739578, ZINC8927502, and ZINC17138581) exhibited the most promising lead requirements for inhibiting VEGFR-2 and suppressing angiogenesis in cancer therapy. This integrated approach, combining pharmacophore modeling, molecular docking, ADMET studies, DFT analysis, and MD simulations, provides valuable insights into the identification of potential anti-cancer agents targeting VEGFR-2.

Cheminformatics-based analysis identified (Z)-2-(2,5-dimethoxy benzylidene)-6-(2-(4-methoxyphenyl)-2-oxoethoxy) benzofuran-3(2H)-one as an inhibitor of Marburg replication by interacting with NP.

The highly pathogenic Marburg virus (MARV) is a member of the Filoviridae family, a non-segmented negative-strand RNA virus. This article represents the computer-aided drug design (CADD) approach for identifying drug-like compounds that prevent the MARV virus disease by inhibiting nucleoprotein, which is responsible for their replication. This study used a wide range of in silico drug design techniques to identify potential drugs. Out of 368 natural compounds, 202 compounds passed ADMET, and molecular docking identified the top two molecules (CID: 1804018 and 5280520) with a high binding affinity of -6.77 and -6.672 kcal/mol, respectively. Both compounds showed interactions with the common amino acid residues SER_216, ARG_215, TYR_135, CYS_195, and ILE_108, which indicates that lead compounds and control ligands interact in the common active site/catalytic site of the protein. The negative binding free energies of CID: 1804018 and 5280520 were -66.01 and -31.29 kcal/mol, respectively. Two lead compounds were re-evaluated using MD modeling techniques, which confirmed CID: 1804018 as the most stable when complexed with the target protein. PC3 of the (Z)-2-(2,5-dimethoxybenzylidene)-6-(2-(4-methoxyphenyl)-2-oxoethoxy) benzofuran-3(2H)-one (CID: 1804018) was 8.74 %, whereas PC3 of the 2'-Hydroxydaidzein (CID: 5280520) was 11.25 %. In this study, (Z)-2-(2,5-dimethoxybenzylidene)-6-(2-(4-methoxyphenyl)-2-oxoethoxy) benzofuran-3(2H)-one (CID: 1804018) unveiled the significant stability of the proteins' binding site in ADMET, Molecular docking, MM-GBSA and MD simulation analysis studies, which also showed a high negative binding free energy value, confirming as the best drug candidate which is found in Angelica archangelica which may potentially inhibit the replication of MARV nucleoprotein.

Recent Advances in Automated Structure-Based De Novo Drug Design.

As the number of determined and predicted protein structures and the size of druglike 'make-on-demand' libraries soar, the time-consuming nature of structure-based computer-aided drug design calls for innovative computational algorithms. De novo drug design introduces in silico heuristics to accelerate searching in the vast chemical space. This review focuses on recent advances in structure-based de novo drug design, ranging from conventional fragment-based methods, evolutionary algorithms, and Metropolis Monte Carlo methods to deep generative models. Due to the historical limitation of de novo drug design generating readily available drug-like molecules, we highlight the synthetic accessibility efforts in each category and the benchmarking strategies taken to validate the proposed framework.

Development of a tetrahydroindazolone-based HDAC6 inhibitor with in-vivo anti-arthritic activity.

Histone deacetylase 6 (HDAC6) induces the expression of pro-inflammatory cytokines in macrophages; therefore, HDAC inhibitors may be beneficial for the treatment of macrophage-associated immune disorders and chronic inflammatory diseases, including atherosclerosis and rheumatoid arthritis. Structure-activity relationship studies were conducted on various phenyl hydroxamate HDAC6 inhibitors with indolone/indazolone-based bi- or tricyclic ring moieties as the cap group aiming to develop novel anti-arthritic drug candidates. Several compounds exhibited nanomolar activity and HDAC6 selectivity greater than 500-fold over HDAC1. Compound 21, a derivative with the tetrahydroindazolone cap group, is a potent HDAC6 inhibitor with an IC50 of 18 nM and 217-fold selectivity over HDAC1 and showed favorable oral bioavailability in animals. Compound 21 increases the acetylation level of tubulin without affecting histone acetylation in cutaneous T-cell lymphoma cells and inhibits TNF-α secretion in LPS-stimulated macrophage cells. The anti-arthritic effects of compound 21 were evaluated using a rat adjuvant-induced arthritis (AIA) model. Treatment with compound 21 significantly reduced the arthritis score, and combination treatment with methotrexate showed a synergistic effect in AIA models. We identified a novel HDAC6 inhibitor, compound 21, with excellent in vivo anti-arthritic efficacy, which can lead to the development of oral anti-arthritic drugs.

Design, in silico Evaluation, and Determination of Antitumor Activity of Potential Inhibitors Against Protein Kinases: Application to BCR-ABL Tyrosine Kinase.

Despite significant progress made over the past two decades in the treatment of chronic myeloid leukemia (CML), there is still an unmet need for effective and safe agents to treat patients with resistance and intolerance to the drugs used in clinic. In this work, we designed 2-arylaminopyrimidine amides of isoxazole-3-carboxylic acid, assessed in silico their inhibitory potential against Bcr-Abl tyrosine kinase, and determined their antitumor activity in K562 (CML), HL-60 (acute promyelocytic leukemia), and HeLa (cervical cancer) cells. Based on the analysis of computational and experimental data, three compounds with the antitumor activity against K562 and HL-60 cells were identified. The lead compound efficiently suppressed the growth of these cells, as evidenced by the low IC50 values of 2.8 ± 0.8 µM (K562) and 3.5 ± 0.2 µM (HL-60). The obtained compounds represent promising basic structures for the design of novel, effective, and safe anticancer drugs able to inhibit the catalytic activity of Bcr-Abl kinase by blocking the ATP-binding site of the enzyme.

Structure-guided discovery of novel dUTPase inhibitors with anti-Nocardia activity by computational design.

The zoonosis caused by Nocardia is increasing seriously. But commonly used antibiotic drugs often lead to resistance. N. seriolae dUTPase (NsdUTPase) plays a key role in the proliferation of Nocardia, and was regarded as a potent drug target. However, there was little report about the NsdUTPase inhibitors. In this study, we discovered a series of novel NsdUTPase inhibitors to fight against Nocardia. The first crystal structure of NsdUTPase was released, and a structure-based computational design was performed. Compounds 4b and 12b exhibited promising activities towards NsdUTPase (IC50 = 0.99 µM and 0.7 µM). In addition, they showed satisfied anti-Nocardia activity (MIC value ranges from 0.5 to 2 mg/L) and low cytotoxicity, which were better than approved drugs oxytetracycline and florfenicol. Molecular modelling study indicated that hydrophobic interaction might be the main contribution for ligand binding. Our results suggested that NsdUTPase inhibitors might be a useful way to repress Nocardia.

Rapid discovery of a new antifoulant: From in silico studies targeting barnacle chitin synthase to efficacy against barnacle settlement.

Due to the adverse environmental impacts of toxic heavy metal-based antifoulants, the screening of environmentally friendly antifoulants has become important for the development of marine antifouling technology. Compared with the traditional lengthy and costly screening method, computer-aided drug design (CADD) offers a promising and efficient solution that can accelerate the screening process of green antifoulants. In this study, we selected barnacle chitin synthase (CHS, an important enzyme for barnacle settlement and development) as the target protein for docking screening. Three CHS genes were identified in the barnacle Amphibalanus amphitrite, and their encoded proteins were found to share a conserved glycosyltransferase domain. Molecular docking of 31,561 marine natural products with AaCHSs revealed that zoanthamine alkaloids had the best binding affinity (-11.8 to -12.6 kcal/mol) to AaCHSs. Considering that the low abundance of zoanthamine alkaloids in marine organisms would limit their application as antifoulants, a marine fungal-derived natural product, mycoepoxydiene (MED), which has a similar chemical structure to zoanthamine alkaloids and the potential for large-scale production by fermentation, was selected and validated for stable binding to AaCHS2L2 using molecular docking and molecular dynamics simulations. Finally, the efficacy of MED in inhibiting cyprid settlement of A. amphitrite was confirmed by a bioassay that demonstrated an EC50 of 1.97 µg/mL, suggesting its potential as an antifoulant candidate. Our research confirmed the reliability of using AaCHSs as antifouling targets and has provided insights for the efficient discovery of green antifoulants by CADD.

Computer-Aided Drug Design Using the Fragment Molecular Orbital Method: Current Status and Future Applications for SBDD.

Owing to the increasing use of computers, computer-aided drug design (CADD) has become an essential component of drug discovery research. In structure-based drug design (SBDD), including inhibitor design and in silico screening of drug target molecules, concordance with wet experimental data is important to provide insights on unique perspectives derived from calculations. Fragment molecular orbital (FMO) method is a quantum chemical method that facilitates precise energy calculations. Fragmentation method makes it possible to apply the quantum chemical method to biological macromolecules for energy calculation based on the electron behavior. Furthermore, interaction energies calculated on a residue-by-residue basis via fragmentation aid in the analysis of interactions between the target and ligand molecule residues and molecular design. In this review, we outline the recent developments in SBDD and FMO methods and highlight the prospects of developing machine learning approaches for large computational data using the FMO method.

Design, synthesis, and anti-tumor activity of derivatives of ring A and C-28 of asiatic acid.

Based on computer-aided drug design (CADD), the active groups of the known active small molecule compounds that can bind to EGFR target protein were analyzed through the molecular docking method. Then, 12 novel asiatic acid derivatives were synthesized by introducing active groups at ring A and C-28 positions of asiatic acid. The structures of these novel compounds were determined by NMR and MS. Furthermore, the anti-tumor activities of these derivatives on human lung cancer cells (A549) and human breast cancer cells (MCF-7) were evaluated by MTT assay. In conclusion, compounds I4 and II3 have stronger anti-cancer activity than parent compounds, the activities were stronger than gefitinib and comparable to afatinib, which may be potential candidate compounds for tumor therapy.

Unveiling Novel Urease Inhibitors for Helicobacter pylori: A Multi-Methodological Approach from Virtual Screening and ADME to Molecular Dynamics Simulations.

Helicobacter pylori (Hp) infections pose a global health challenge demanding innovative therapeutic strategies by which to eradicate them. Urease, a key Hp virulence factor hydrolyzes urea, facilitating bacterial survival in the acidic gastric environment. In this study, a multi-methodological approach combining pharmacophore- and structure-based virtual screening, molecular dynamics simulations, and MM-GBSA calculations was employed to identify novel inhibitors for Hp urease (HpU). A refined dataset of 8,271,505 small molecules from the ZINC15 database underwent pharmacokinetic and physicochemical filtering, resulting in 16% of compounds for pharmacophore-based virtual screening. Molecular docking simulations were performed in successive stages, utilizing HTVS, SP, and XP algorithms. Subsequent energetic re-scoring with MM-GBSA identified promising candidates interacting with distinct urease variants. Lys219, a residue critical for urea catalysis at the urease binding site, can manifest in two forms, neutral (LYN) or carbamylated (KCX). Notably, the evaluated molecules demonstrated different interaction and energetic patterns in both protein variants. Further evaluation through ADMET predictions highlighted compounds with favorable pharmacological profiles, leading to the identification of 15 candidates. Molecular dynamics simulations revealed comparable structural stability to the control DJM, with candidates 5, 8 and 12 (CA5, CA8, and CA12, respectively) exhibiting the lowest binding free energies. These inhibitors suggest a chelating capacity that is crucial for urease inhibition. The analysis underscores the potential of CA5, CA8, and CA12 as novel HpU inhibitors. Finally, we compare our candidates with the chemical space of urease inhibitors finding physicochemical similarities with potent agents such as thiourea.

Structure-Aided Computational Design of Triazole-Based Targeted Covalent Inhibitors of Cruzipain.

Cruzipain (CZP), the major cysteine protease present in T. cruzi, the ethiological agent of Chagas disease, has attracted particular attention as a therapeutic target for the development of targeted covalent inhibitors (TCI). The vast chemical space associated with the enormous molecular diversity feasible to explore by means of modern synthetic approaches allows the design of CZP inhibitors capable of exhibiting not only an efficient enzyme inhibition but also an adequate translation to anti-T. cruzi activity. In this work, a computer-aided design strategy was developed to combinatorially construct and screen large libraries of 1,4-disubstituted 1,2,3-triazole analogues, further identifying a selected set of candidates for advancement towards synthetic and biological activity evaluation stages. In this way, a virtual molecular library comprising more than 75 thousand diverse and synthetically feasible analogues was studied by means of molecular docking and molecular dynamic simulations in the search of potential TCI of CZP, guiding the synthetic efforts towards a subset of 48 candidates. These were synthesized by applying a Cu(I)-catalyzed azide-alkyne cycloaddition (CuAAC) centered synthetic scheme, resulting in moderate to good yields and leading to the identification of 12 hits selectively inhibiting CZP activity with IC50 in the low micromolar range. Furthermore, four triazole derivatives showed good anti-T. cruzi inhibition when studied at 50 µM; and Ald-6 excelled for its high antitrypanocidal activity and low cytotoxicity, exhibiting complete in vitro biological activity translation from CZP to T. cruzi. Overall, not only Ald-6 merits further advancement to preclinical in vivo studies, but these findings also shed light on a valuable chemical space where molecular diversity might be explored in the search for efficient triazole-based antichagasic agents.

Lysine-Specific Demethylase 1 Inhibitors: A Comprehensive Review Utilizing Computer-Aided Drug Design Technologies.

Lysine-specific demethylase 1 (LSD1/KDM1A) has emerged as a promising therapeutic target for treating various cancers (such as breast cancer, liver cancer, etc.) and other diseases (blood diseases, cardiovascular diseases, etc.), owing to its observed overexpression, thereby presenting significant opportunities in drug development. Since its discovery in 2004, extensive research has been conducted on LSD1 inhibitors, with notable contributions from computational approaches. This review systematically summarizes LSD1 inhibitors investigated through computer-aided drug design (CADD) technologies since 2010, showcasing a diverse range of chemical scaffolds, including phenelzine derivatives, tranylcypromine (abbreviated as TCP or 2-PCPA) derivatives, nitrogen-containing heterocyclic (pyridine, pyrimidine, azole, thieno[3,2-b]pyrrole, indole, quinoline and benzoxazole) derivatives, natural products (including sanguinarine, phenolic compounds and resveratrol derivatives, flavonoids and other natural products) and others (including thiourea compounds, Fenoldopam and Raloxifene, (4-cyanophenyl)glycine derivatives, propargylamine and benzohydrazide derivatives and inhibitors discovered through AI techniques). Computational techniques, such as virtual screening, molecular docking and 3D-QSAR models, have played a pivotal role in elucidating the interactions between these inhibitors and LSD1. Moreover, the integration of cutting-edge technologies such as artificial intelligence holds promise in facilitating the discovery of novel LSD1 inhibitors. The comprehensive insights presented in this review aim to provide valuable information for advancing further research on LSD1 inhibitors.

Benchmarking applicability of medium-resolution cryo-EM protein structures for structure-based drug design.

Cryo-electron microscopy (cryo-EM) is gaining large attention for high-resolution protein structure determination in solutions. However, a very high percentage of cryo-EM structures correspond to resolutions of 3-5 Å, making the structures difficult to be used in in silico drug design. In this study, we analyze how useful cryo-EM protein structures are for in silico drug design by evaluating ligand docking accuracy. From realistic cross-docking scenarios using medium resolution (3-5 Å) cryo-EM structures and a popular docking tool Autodock-Vina, only 20% of docking succeeded, when the success rate doubles in the same kind of cross-docking but using high-resolution (<2 Å) crystal structures instead. We decipher the reason for failures by decomposing the contribution from resolution-dependent and independent factors. The heterogeneity in the protein side-chain and backbone conformations is identified as the major resolution-dependent factor causing docking difficulty from our analysis, while intrinsic receptor flexibility mainly comprises the resolution-independent factor. We demonstrate the flexibility implementation in current ligand docking tools is able to rescue only a portion of failures (10%), and the limited performance was majorly due to potential structural errors than conformational changes. Our work suggests the strong necessity of more robust method developments on ligand docking and EM modeling techniques in order to fully utilize cryo-EM structures for in silico drug design.

GPU-specific algorithms for improved solute sampling in grand canonical Monte Carlo simulations.

The Grand Canonical Monte Carlo (GCMC) ensemble defined by the excess chemical potential, µex , volume, and temperature, in the context of molecular simulations allows for variations in the number of particles in the system. In practice, GCMC simulations have been widely applied for the sampling of rare gasses and water, but limited in the context of larger molecules. To overcome this limitation, the oscillating µex GCMC method was introduced and shown to be of utility for sampling small solutes, such as formamide, propane, and benzene, as well as for ionic species such as monocations, acetate, and methylammonium. However, the acceptance of GCMC insertions is low, and the method is computationally demanding. In the present study, we improved the sampling efficiency of the GCMC method using known cavity-bias and configurational-bias algorithms in the context of GPU architecture. Specifically, for GCMC simulations of aqueous solution systems, the configurational-bias algorithm was extended by applying system partitioning in conjunction with a random interval extraction algorithm, thereby improving the efficiency in a highly parallel computing environment. The method is parallelized on the GPU using CUDA and OpenCL, allowing for the code to run on both Nvidia and AMD GPUs, respectively. Notably, the method is particularly well suited for GPU computing as the large number of threads allows for simultaneous sampling of a large number of configurations during insertion attempts without additional computational overhead. In addition, the partitioning scheme allows for simultaneous insertion attempts for large systems, offering considerable efficiency. Calculations on the BK Channel, a transporter, including a lipid bilayer with over 760,000 atoms, show a speed up of ~53-fold through the use of system partitioning. The improved algorithm is then combined with an enhanced µex oscillation protocol and shown to be of utility in the context of the site-identification by ligand competitive saturation (SILCS) co-solvent sampling approach as illustrated through application to the protein CDK2.

Design, synthesis, and biological testing of biphenylmethyloxazole inhibitors targeting HIV-1 reverse transcriptase.

We report non-nucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs) using a biphenylmethyloxazole pharmacophore. A crystal structure of benzyloxazole 1 was obtained and suggested the potential viability of biphenyl analogues. In particular, 6a, 6b, and 7 turned out to be potent NNRTIs with low-nanomolar activity in enzyme inhibition and infected T-cell assays, and with low cytotoxicity. Though modeling further suggested that analogues with fluorosulfate and epoxide warheads might provide covalent modification of Tyr188, synthesis and testing did not find evidence for this outcome.

MASSA Algorithm: an automated rational sampling of training and test subsets for QSAR modeling.

QSAR models capable of predicting biological, toxicity, and pharmacokinetic properties were widely used to search lead bioactive molecules in chemical databases. The dataset's preparation to build these models has a strong influence on the quality of the generated models, and sampling requires that the original dataset be divided into training (for model training) and test (for statistical evaluation) sets. This sampling can be done randomly or rationally, but the rational division is superior. In this paper, we present MASSA, a Python tool that can be used to automatically sample datasets by exploring the biological, physicochemical, and structural spaces of molecules using PCA, HCA, and K-modes. The proposed algorithm is very useful when the variables used for QSAR are not available or to construct multiple QSAR models with the same training and test sets, producing models with lower variability and better values for validation metrics. These results were obtained even when the descriptors used in the QSAR/QSPR were different from those used in the separation of training and test sets, indicating that this tool can be used to build models for more than one QSAR/QSPR technique. Finally, this tool also generates useful graphical representations that can provide insights into the data.