Molecular Dynamics Simulations and Computer-Aided Drug Discovery

MotivationThe ligand/protein binding interaction is typically investigated by docking and molecular dynamics (MD) simulations. In particular, docking-based virtual screening (VS) is used to select the best ligands from database of thousands of compounds, while MD calculations assess the energy stability of the ligand/protein binding complexes. Considering the broad use of these techniques, it is of great demand to have one single software that allows a combined and fast analysis of VS and MD results. With this in mind, we have developed the Drug Discovery Tool (DDT) that is an intuitive graphics user interface able to provide structural data and physico-chemical information on the ligand/protein interaction.ResultsDDT is designed as a plugin for the Visual Molecular Dynamics (VMD) software and is able to manage a large number of ligand/protein complexes obtained from AutoDock4 (AD4) docking calculations and MD simulations. DDT delivers four main outcomes: i) ligands ranking based on an energy score; ii) ligand ranking based on a ligands’ conformation cluster analysis; iii) identification of the aminoacids forming the most occurrent interactions with the ligands; iv) plot of the ligands’ center-of-mass coordinates in the Cartesian space. The flexibility of the software allows saving the best ligand/protein complexes using a number of user-defined options.Availability and implementationDDT_site_1 (alternative DDT_site_2); the DDT tutorial movie is available here.Supplementary informationSupplementary data are available at Bioinformatics online.

Inhibition of dipeptidyl peptidase-4 (DPP4) activity has emerged as a promising therapeutic approach for the treatment of type 2 diabetes mellitus (T2DM). Bioinformatics-driven approaches have emerged as crucial tools in drug discovery. Molecular docking and molecular dynamics (MD) simulations are effective tools in drug discovery, as they reduce the time and cost associated with experimental screening. In this study, we employed structure-assisted in-silico methods, including molecular docking and MD simulations, to identify SRT2183, a small molecule that may potentially inhibit the activity of DPP4 enzyme. The interaction between the small molecule "SRT2183" and DPP4 exhibited a binding affinity of −9.9 Kcal/Mol, leading to the formation of hydrogen bonds with the amino acid residues MET348, SER376, and THR351 of DPP4. The MD simulations over a period of 100 ns indicated stable protein-ligand interactions, with no significant conformational rearrangements observed within the simulated timeframe. In conclusion, our results suggest that the small molecule SRT2183 may have the potential to inhibit the DPP4 enzyme and pave the way for the therapeutics of T2DM.

BackgroundThe virtual screening (VS) of lead compounds using molecular docking and pharmacophore detection is now an important tool in drug discovery. VS tasks typically require a combination of several software tools and a molecular graphics system. Thus, the integration of all the requisite tools in a single operating environment could reduce the complexity of running VS experiments. However, only a few freely available integrated software platforms have been developed.ResultsA free open-source platform, IVSPlat 1.0, was developed in this study for the management and automation of VS tasks. We integrated several VS-related programs into a molecular graphics system to provide a comprehensive platform for the solution of VS tasks based on molecular docking, pharmacophore detection, and a combination of both methods. This tool can be used to visualize intermediate and final results of the VS execution, while also providing a clustering tool for the analysis of VS results. A case study was conducted to demonstrate the applicability of this platform.ConclusionsIVSPlat 1.0 provides a plug-in-based solution for the management, automation, and visualization of VS tasks. IVSPlat 1.0 is an open framework that allows the integration of extra software to extend its functionality and modified versions can be freely distributed. The open source code and documentation are available at http://kyc.nenu.edu.cn/IVSPlat/.

The ability to perform long, accurate molecular dynamics (MD) simulations involving proteins and other biological macromolecules could in principle lead to important scientific advances and provide a powerful new tool for drug discovery. A wide range of biologically interesting phenomena, however, occur over time scales on the order of a millisecond - several orders of magnitude beyond the duration of the longest current MD simulations. Our research group is currently building a specialized, massively parallel machine, called Anton, which should soon be capable of executing millisecond-scale MD simulations of proteins at an atomic level of detail. Antonpsilas highly accelerated execution of such simulations is attributable in large part to specialized logic for the high-speed calculation of pairwise interactions between particles and/or gridpoints separated by no more than some specified cutoff radius.In particular, each of Anton's 512 ASICs, which are implemented using 90-nm technology,includes a "high-throughput interaction subsystem" incorporating 32 highly specialized pipelines running at 800 MHz. During every cycle, each of these pipelines produces a pairwise-interaction result that would require approximately 50 arithmetic operations to calculate on a general-purpose processor. Novel algorithms and architectural features are used to greatly reduce the requirements for inter- and intra-chip communication, allowing Anton to feed these pipelines and collect their results at a speed sufficient to take advantage of the machine's computational power. The ASIC also includes a "flexible subsystem" based on eight programmable "geometry cores", each containing eight arithmetic pipelines. This talk will provide an overview of our work on parallel algorithms and machine architectures for high-speed MD simulation, with special attention to the respective roles of specialized vs. general-purpose hardware, and to the techniques used to minimize communication at various levels within the system.

Microscale thermophoresis (MST) and spectral shift (SpS) in drug discovery

IntroductionDengue Fever continues to pose a global threat due to the widespread distribution of its vector mosquitoes, Aedes aegypti and Aedes albopictus. While the WHO-approved vaccine, Dengvaxia, and antiviral treatments like Balapiravir and Celgosivir are available, challenges such as drug resistance, reduced efficacy, and high treatment costs persist. This study aims to identify novel potential inhibitors of the Dengue virus (DENV) using an integrative drug discovery approach encompassing machine learning and molecular docking techniques.MethodUtilizing a dataset of 21,250 bioactive compounds from PubChem (AID: 651640), alongside a total of 1,444 descriptors generated using PaDEL, we trained various models such as Support Vector Machine, Random Forest, k-nearest neighbors, Logistic Regression, and Gaussian Naïve Bayes. The top-performing model was used to predict active compounds, followed by molecular docking performed using AutoDock Vina. The detailed interactions, toxicity, stability, and conformational changes of selected compounds were assessed through protein-ligand interaction studies, molecular dynamics (MD) simulations, and binding free energy calculations.ResultsWe implemented a robust three-dataset splitting strategy, employing the Logistic Regression algorithm, which achieved an accuracy of 94%. The model successfully predicted 18 known DENV inhibitors, with 11 identified as active, paving the way for further exploration of 2683 new compounds from the ZINC and EANPDB databases. Subsequent molecular docking studies were performed on the NS2B/NS3 protease, an enzyme essential in viral replication. ZINC95485940, ZINC38628344, 2′,4′-dihydroxychalcone and ZINC14441502 demonstrated a high binding affinity of −8.1, −8.5, −8.6, and −8.0 kcal/mol, respectively, exhibiting stable interactions with His51, Ser135, Leu128, Pro132, Ser131, Tyr161, and Asp75 within the active site, which are critical residues involved in inhibition. Molecular dynamics simulations coupled with MMPBSA further elucidated the stability, making it a promising candidate for drug development.ConclusionOverall, this integrative approach, combining machine learning, molecular docking, and dynamics simulations, highlights the strength and utility of computational tools in drug discovery. It suggests a promising pathway for the rapid identification and development of novel antiviral drugs against DENV. These in silico findings provide a strong foundation for future experimental validations and in-vitro studies aimed at fighting DENV.

Beta 2 glycoprotein I (β2GPI) is a major antigen for autoantibodies present in antiphospholipid antibody syndrome (APS). β2GPI is a single polypeptide with five repeated domains and different conformations. The activated J-shaped conformation of β2GPI binds to negatively charged phospholipids in the membrane via the fifth domain and causes blood clotting reactions. We applied a drug repurposing strategy using virtual screening and molecular dynamics to find the best FDA drugs against the fifth domain of β2GPI. In the first phase, FDA drugs that had the most favorable ΔG with the fifth domain of β2GPI were selected by virtual screening. Among these drugs that had the most favorable ΔG, Vorapaxar and Antrafenine were selected for molecular dynamics (MD) simulation studies. MD simulation was performed to evaluate the stability of Vorapaxar and Antrafenine complexes and the effect of the two drugs on protein conformation. Also, MD simulation was done to investigate the effect of Antrafenine and Vorapaxar on the binding of β2GPI to the platelet model membrane. According to the results, Vorapaxar and Antrafenine were bound to the protein with the favorable binding energy (Vorapaxar and Antrafenine binding energies are - 49.641 and - 38.803kcal/mol, respectively). In this study, it was shown that unlike protein alone and protein in the Antrafenine complex, the protein in the Vorapaxar complex was completely separated from the model membrane after 350ns. Moreover, Vorapaxar led to more changes in the activated J-shape of β2GPI. Thus, Vorapaxar can be a suitable candidate for further investigations on the treatment of APS.

The role of computational tools in drug discovery and development is becoming increasingly important due to the rapid development of computing power and advancements in computational chemistry and biology, improving research efficiency and reducing the costs and potential risks of preclinical and clinical trials. Machine learning, especially deep learning, a subfield of artificial intelligence (AI), has demonstrated significant advantages in drug discovery and development, including high-throughput and virtual screening, ab initio design of drug molecules, and solving difficult organic syntheses. This review summarizes AI technologies used in drug discovery and development, including their roles in drug screening, design, and solving the challenges of clinical trials. Finally, it discusses the challenges of drug discovery and development based on AI technologies, as well as potential future directions.

Positron emission tomography neuroreceptor imaging as a tool in drug discovery, research and development

Infectious diseases, caused by bacteria, viruses, fungi, protozoa, or parasites are among the leading causes of death worldwide. The efficacy of all current anti-infectives is threatened by the spread of drug resistance factors, with some already made ineffective, and, as a result, there is a pressing need for a new pipeline of robust anti-infective drugs. In silico approaches, such as virtual high throughput screening and de novo structure-based rational drug design, have been established as powerful tools in drug discovery. In this review, we explore the exciting opportunities for antimalarial, antiviral, and antibacterial drug discovery arising from the new paradigm of structure-based drug design.

Proteins are essential units of life that govern several functions. Understanding their behavior is closely related to their conformations, native folds, and change in conformations. Thus, the dynamic information of protein becomes essential to understand its properties at the molecular level. The molecular dynamics (MD) simulation approach provides atomistic-level dynamic information about proteins. However, more extended or complex MD simulations of protein are challenging to analyze and to gather meaningful confirmation from several snapshots of the dynamic system. To achieve it, i.e., analyzing MD simulation data, Markov State Model (MSM) is a powerful tool that has a statistical background. It represents the MD simulation system as a combination of finite memoryless states, i.e., states that are not dependent on prior states and transition probability among such states. MSM applications have grown from peptides to membrane protein simulations. The present book chapter sheds light on MD simulation’s role in protein dynamics and why MSM is required. The brief theoretical aspects of MSM techniques are demonstrated. Lastly, the chapter discusses the application of MSM in different protein folding and dynamics.

In silico technology is an essential tool in parasitology, enabling rapid and cost-effective drug discovery, vaccine design, diagnosis, and prediction of host-parasite interactions. The computational tools, include molecular docking and dynamics simulations, virtual screening, homology modeling, and machine learning, allow the identification of novel therapeutic targets and potential drug candidates before laboratory validation.

Structure-based virtual screening is a key tool in early drug discovery, with growing interest in the screening of multi-billion chemical compound libraries. However, the success of virtual screening crucially depends on the accuracy of the binding pose and binding affinity predicted by computational docking. Here we develop a highly accurate structure-based virtual screen method, RosettaVS, for predicting docking poses and binding affinities. Our approach outperforms other state-of-the-art methods on a wide range of benchmarks, partially due to our ability to model receptor flexibility. We incorporate this into a new open-source artificial intelligence accelerated virtual screening platform for drug discovery. Using this platform, we screen multi-billion compound libraries against two unrelated targets, a ubiquitin ligase target KLHDC2 and the human voltage-gated sodium channel NaV1.7. For both targets, we discover hit compounds, including seven hits (14% hit rate) to KLHDC2 and four hits (44% hit rate) to NaV1.7, all with single digit micromolar binding affinities. Screening in both cases is completed in less than seven days. Finally, a high resolution X-ray crystallographic structure validates the predicted docking pose for the KLHDC2 ligand complex, demonstrating the effectiveness of our method in lead discovery.

Design, synthesis, and computational evaluation of ester prodrugs of isoguvacine as potential antiseizure medications.

Drug action can be rationalized as interaction of a molecule with proteins in a regulatory network of targets from a specific biological system. Both drug and side effects are often governed by interaction of the drug molecule with many, often unrelated, targets. Accordingly, arrays of protein-ligand interaction data from numerous in vitro profiling assays today provide growing evidence of polypharmacological drug interactions, even for marketed drugs. In vitro off-target profiling has therefore become an important tool in early drug discovery to learn about potential off-target liabilities, which are sometimes beneficial, but more often safety relevant. The rapidly developing field of in silico profiling approaches is complementing in vitro profiling. These approaches capitalize from large amounts of biochemical data from multiple sources to be exploited for optimizing undesirable side effects in pharmaceutical research. Therefore, current in silico profiling models are nowadays perceived as valuable tools in drug discovery, and promise a platform to support optimally informed decisions.