Large-scale integrated super-computing platform for next generation virtual drug discovery

The New Era of Drug Discovery: The Power of Computer-aided Drug Design (CADD)

Abstract:

Drug design and discovery is a process that requires high financial costs and is timeconsuming. For many years, this process focused on empirical pharmacology. However, over the years, the target-based approach allowed a significant discovery in this field, initiating the rational design era. In view, to decrease the time and financial cost, rational drug design is benefited by increasing computer engineering and software development, and computer-aided drug design (CADD) emerges as a promising alternative. Since the 1970s, this approach has been able to identify many important and revolutionary compounds, like protease inhibitors, antibiotics, and others. Many anticancer compounds identified through this approach have shown their importance, being CADD essential in any drug discovery campaign. Thus, this perspective will present the prominent successful cases utilizing this approach and entering into the next stage of drug design. We believe that drug discovery will follow the progress in bioinformatics, using high-performance computing with molecular dynamics protocols faster and more effectively. In addition, artificial intelligence and machine learning will be the next process in the rational design of new drugs. Here, we hope that this paper generates new ideas and instigates research groups worldwide to use these methods and stimulate progress in drug design.

An Improved Method for P2X7R Antagonist Screening

ATP physiologically activates the P2X7 receptor (P2X7R), a member of the P2X ionotropic receptor family. When activated by high concentrations of ATP (i.e., at inflammation sites), this receptor is capable of forming a pore that allows molecules of up to 900 Da to pass through. This receptor is upregulated in several diseases, particularly leukemia, rheumatoid arthritis and Alzheimer's disease. A selective antagonist of this receptor could be useful in the treatment of P2X7R activation-related diseases. In the present study, we have evaluated several parameters using in vitro protocols to validate a high-throughput screening (HTS) method to identify P2X7R antagonists. We generated dose-response curves to determine the EC₅₀ value of the known agonist ATP and the IC_s50 values for the known antagonists Brilliant Blue G (BBG) and oxidized ATP (OATP). The values obtained were consistent with those found in the literature (0.7 ± 0.07 mM, 1.3-2.6 mM and 173-285 μM for ATP, BBG and OATP, respectively). The Z-factor, an important statistical tool that can be used to validate the robustness and suitability of an HTS assay, was 0.635 for PI uptake and 0.867 for LY uptake. No inter-operator variation was observed, and the results obtained using our improved method were reproducible. Our data indicate that our assay is suitable for the selective and reliable evaluation of P2X7 activity in multiwell plates using spectrophotometry-based methodology. This method might improve the high-throughput screening of conventional chemical or natural product libraries for possible candidate P2X7R antagonist or agonist

Prediction of the binding affinity of compounds with diverse scaffolds by MP-CAFEE

Accurate methods to predict the binding affinities of compounds for target molecules are powerful tools in structure-based drug design (SBDD). A recently developed method called massively parallel computation of absolute binding free energy with a well-equilibrated system (MP-CAFEE) successfully predicted the binding affinities of compounds with relatively similar scaffolds. We investigate the applicability of MP-CAFEE for predicting the affinity of compounds having more diverse scaffolds for the target p38α, a mitogen-activated protein kinase. The calculated and experimental binding affinities correlate well, showing that MP-CAFEE can accurately rank the compounds with diverse scaffolds. We propose a method to determine the optimal number of sampling runs with respect to a predefined level of accuracy, which is established according to the stage in the SBDD process being considered. The optimal number of sampling runs for two key stages-lead identification and lead optimization-is estimated to be five and eight or more, respectively, in our model system using Cochrans sample size formula.

New insights for drug design from the X-ray crystallographic structures of G-protein-coupled receptors

Methodological advances in X-ray crystallography have made possible the recent solution of X-ray structures of pharmaceutically important G protein-coupled receptors (GPCRs), including receptors for biogenic amines, peptides, a nucleoside, and a sphingolipid. These high-resolution structures have greatly increased our understanding of ligand recognition and receptor activation. Conformational changes associated with activation common to several receptors entail outward movements of the intracellular side of transmembrane helix 6 (TM6) and movements of TM5 toward TM6. Movements associated with specific agonists or receptors have also been described [e.g., extracellular loop (EL) 3 in the A(2A) adenosine receptor]. The binding sites of different receptors partly overlap but differ significantly in ligand orientation, depth, and breadth of contact areas in TM regions and the involvement of the ELs. A current challenge is how to use this structural information for the rational design of novel potent and selective ligands. For example, new chemotypes were discovered as antagonists of various GPCRs by subjecting chemical libraries to in silico docking in the X-ray structures. The vast majority of GPCR structures and their ligand complexes are still unsolved, and no structures are known outside of family A GPCRs. Molecular modeling, informed by supporting information from site-directed mutagenesis and structure-activity relationships, has been validated as a useful tool to extend structural insights to related GPCRs and to analyze docking of other ligands in already crystallized GPCRs.

Established and emerging trends in computational drug discovery in the structural genomics era

Bioinformatics and chemoinformatics approaches contribute to hit discovery, hit-to-lead optimization, safety profiling, and target identification and enhance our overall understanding of the health and disease states. A vast repertoire of computational methods has been reported and increasingly combined in order to address more and more challenging targets or complex molecular mechanisms in the context of large-scale integration of structure and bioactivity data produced by private and public drug research. This review explores some key computational methods directly linked to drug discovery and chemical biology with a special emphasis on compound collection preparation, virtual screening, protein docking, and systems pharmacology. A list of generally freely available software packages and online resources is provided, and examples of successful applications are briefly commented upon.