Identification of potential human beta-secretase 1 inhibitors by hierarchical virtual screening and molecular dynamics.

Alzheimer's disease (AD) is a neurodegenerative pathology responsible for 70% of dementia cases worldwide. Despite its relevance, the few drugs available for the treatment of this disease offer only symptomatic relief, with limited efficacy and serious adverse effects. The most accepted hypothesis about the pathogenesis involves the aggregation and deposition of ß-amyloid peptides, mainly in the cerebral cortex and hippocampus, through the catalytic action of beta-secretase 1 (BACE-1), making this enzyme a promising target for the development of new drugs. In order to prioritize candidates for BACE-1 inhibitors, a hierarchical virtual screening by pharmacophore model and molecular docking was performed against the 216,833 molecules contained in several databases. Our previously built pharmacophore model was used for the first filtering step, which resulted in the selection of 399 molecules. The remaining molecules were filtered through molecular docking with GOLD 5.4.0. In this step, molecules with scoring values ​​greater than the mean plus standard deviation were evaluated for commercial availability and absence of asymmetric centers. Four molecules were selected and evaluated for mutagenic potential by the AMES test with the help of the pkCSM server. Finally, they were tested against the descriptors on Lipinski and Veber rules, and ZINC01589617 (QFIT = 56.52/Score = 44.95) satisfied all requirements, being subjected to molecular dynamics simulations (t = 100 ns) in order to obtain robust data on the mode of bonding and profile of intermolecular interactions. Those in silico strategies demonstrated that ZINC01589617 is a potential candidate for biological tests.Communicated by Ramaswamy H. Sarma.

Evaluation of Docking Machine Learning and Molecular Dynamics Methodologies for DNA-Ligand Systems.

DNA is a molecular target for the treatment of several diseases, including cancer, but there are few docking methodologies exploring the interactions between nucleic acids with DNA intercalating agents. Different docking methodologies, such as AutoDock Vina, DOCK 6, and Consensus, implemented into Molecular Architect (MolAr), were evaluated for their ability to analyze those interactions, considering visual inspection, redocking, and ROC curve. Ligands were refined by Parametric Method 7 (PM7), and ligands and decoys were docked into the minor DNA groove (PDB code: 1VZK). As a result, the area under the ROC curve (AUC-ROC) was 0.98, 0.88, and 0.99 for AutoDock Vina, DOCK 6, and Consensus methodologies, respectively. In addition, we proposed a machine learning model to determine the experimental ∆Tm value, which found a 0.84 R2 score. Finally, the selected ligands mono imidazole lexitropsin (42), netropsin (45), and N,N'-(1H-pyrrole-2,5-diyldi-4,1-phenylene)dibenzenecarboximidamide (51) were submitted to Molecular Dynamic Simulations (MD) through NAMD software to evaluate their equilibrium binding pose into the groove. In conclusion, the use of MolAr improves the docking results obtained with other methodologies, is a suitable methodology to use in the DNA system and was proven to be a valuable tool to estimate the ∆Tm experimental values of DNA intercalating agents.

Structure-Based Virtual Screening: From Classical to Artificial Intelligence.

The drug development process is a major challenge in the pharmaceutical industry since it takes a substantial amount of time and money to move through all the phases of developing of a new drug. One extensively used method to minimize the cost and time for the drug development process is computer-aided drug design (CADD). CADD allows better focusing on experiments, which can reduce the time and cost involved in researching new drugs. In this context, structure-based virtual screening (SBVS) is robust and useful and is one of the most promising in silico techniques for drug design. SBVS attempts to predict the best interaction mode between two molecules to form a stable complex, and it uses scoring functions to estimate the force of non-covalent interactions between a ligand and molecular target. Thus, scoring functions are the main reason for the success or failure of SBVS software. Many software programs are used to perform SBVS, and since they use different algorithms, it is possible to obtain different results from different software using the same input. In the last decade, a new technique of SBVS called consensus virtual screening (CVS) has been used in some studies to increase the accuracy of SBVS and to reduce the false positives obtained in these experiments. An indispensable condition to be able to utilize SBVS is the availability of a 3D structure of the target protein. Some virtual databases, such as the Protein Data Bank, have been created to store the 3D structures of molecules. However, sometimes it is not possible to experimentally obtain the 3D structure. In this situation, the homology modeling methodology allows the prediction of the 3D structure of a protein from its amino acid sequence. This review presents an overview of the challenges involved in the use of CADD to perform SBVS, the areas where CADD tools support SBVS, a comparison between the most commonly used tools, and the techniques currently used in an attempt to reduce the time and cost in the drug development process. Finally, the final considerations demonstrate the importance of using SBVS in the drug development process.

Molecular Architect: A User-Friendly Workflow for Virtual Screening.

Computer-assisted drug design (CADD) methods have greatly contributed to the development of new drugs. Among CADD methodologies, virtual screening (VS) can enrich the compound collection with molecules that have the desired physicochemical and pharmacophoric characteristics that are needed to become drugs. Many free tools are available for this purpose, but they are difficult to use and do not have a graphical user interface. Furthermore, several free tools must be used to carry out the entire VS process, requiring the user to process the results of one software program so that they can be used in another program, adding a potential source of human error. Moreover, some software programs require knowledge of advanced computational skills, such as programming languages. This context has motivated us to develop Molecular Architect (MolAr). MolAr is a workflow with a simple and intuitive interface that acts in an integrated and automated form to perform the entire VS process, from protein preparation (homology modeling and protonation state) to virtual screening. MolAr carries out VS through AutoDock Vina, DOCK 6, or a consensus of the two. Two case studies were conducted to demonstrate the performance of MolAr. In the first study, the feasibility of using MolAr for DNA-ligand systems was assessed. Both AutoDock Vina and DOCK 6 showed good results in performing VS in DNA-ligand systems. However, the use of consensus virtual screening was able to enrich the results. According to the area under the ROC curve and the enrichment factors, consensus VS was better able to predict the positions of the active ligands. The second case study was performed on 8 targets from the DUD-E database and 10 active ligands for each target. The results demonstrated that using the final ligand conformation provided by AutoDock Vina as an input for DOCK 6 improved the DOCK 6 ROC curves by up to 42% in VS. These case studies demonstrated that MolAr is capable conducting the VS process and is an easy-to-use and effective tool. MolAr is available for download free of charge at http: //www.drugdiscovery.com.br/software/.

Docking-based virtual screening of Brazilian natural compounds using the OOMT as the pharmacological target database.

The demand for new therapies has encouraged the development of faster and cheaper methods of drug design. Considering the number of potential biological targets for new drugs, the docking-based virtual screening (DBVS) approach has occupied a prominent role among modern strategies for identifying new bioactive substances. Some tools have been developed to validate docking methodologies and identify false positives, such as the receiver operating characteristic (ROC) curve. In this context, a database with 31 molecular targets called the Our Own Molecular Targets Data Bank (OOMT) was validated using the root-mean-square deviation (RMSD) and the area under the ROC curve (AUC) with two different docking methodologies: AutoDock Vina and DOCK 6. Sixteen molecular targets showed AUC values of >0.8, and those targets were selected for molecular docking studies. The drug-likeness properties were then determined for 473 Brazilian natural compounds that were obtained from the ZINC database. Ninety-six compounds showed similar drug-likeness property values to the marked drugs (positive values). These compounds were submitted to DBVS for 16 molecular targets. Our results showed that AutoDock Vina was more appropriate than DOCK 6 for performing DBVS experiments. Furthermore, this work suggests that three compounds-ZINC13513540, ZINC06041137, and ZINC1342926-are inhibitors of the three molecular targets 1AGW, 2ZOQ, and 3EYG, respectively, which are associated with cancer. Finally, since ZINC and the PDB were solely created to store biomolecule structures, their utilization requires the application of filters to improve the first steps of the drug development process. Graphical Abstract Evaluation of docking methods used for virtual screening.

Octopus: a platform for the virtual high-throughput screening of a pool of compounds against a set of molecular targets.

Octopus is an automated workflow management tool that is scalable for virtual high-throughput screening (vHTS). It integrates MOPAC2016, MGLTools, PyMOL, and AutoDock Vina. In contrast to other platforms, Octopus can perform docking simulations of an unlimited number of compounds into a set of molecular targets. After generating the ligands in a drawing package in the Protein Data Bank (PDB) format, Octopus can carry out geometry refinement using the semi-empirical method PM7 implemented in MOPAC2016. Docking simulations can be performed using AutoDock Vina and can utilize the Our Own Molecular Targets (OOMT) databank. Finally, the proposed software compiles the best binding energies into a standard table. Here, we describe two successful case studies that were verified by biological assay. In the first case study, the vHTS process was carried out for 22 (phenylamino)urea derivatives. The vHTS process identified a metalloprotease with the PDB code 1GKC as a molecular target for derivative LE&007. In a biological assay, compound LE&007 was found to inhibit 80% of the activity of this enzyme. In the second case study, compound Tx001 was submitted to the Octopus routine, and the results suggested that Plasmodium falciparum ATP6 (PfATP6) as a molecular target for this compound. Following an antimalarial assay, Tx001 was found to have an inhibitory concentration (IC50) of 8.2 µM against PfATP6. These successful examples illustrate the utility of this software for finding appropriate molecular targets for compounds. Hits can then be identified and optimized as new antineoplastic and antimalarial drugs. Finally, Octopus has a friendly Linux-based user interface, and is available at www.drugdiscovery.com.br . Graphical Abstract Octopus: A platform for inverse virtual screening (IVS) to search new molecular targets for drugs.

A Ranking Approach for Probe Selection and Classification of Microarray Data with Artificial Neural Networks.

Acute leukemia classification into its myeloid and lymphoblastic subtypes is usually accomplished according to the morphology of the tumor. Nevertheless, the subtypes may have similar histopathological appearance, making screening procedures difficult. In addition, approximately one-third of acute myeloid leukemias are characterized by aberrant cytoplasmic localization of nucleophosmin (NPMc(+)), where the majority has a normal karyotype. This work is based on two DNA microarray datasets, available publicly, to differentiate leukemia subtypes. The datasets were split into training and test sets, and feature selection methods were applied. Artificial neural network classifiers were developed to compare the feature selection methods. For the first dataset, 50 genes selected using the best classifier was able to classify all patients in the test set. For the second dataset, five genes yielded 97.5% accuracy in the test set.