Computational Modeling in Virus Infections and Virtual Screening, Docking, and Molecular Dynamics in Drug Design
Ontology highlight
ABSTRACT:
SUBMITTER: da Silva F
PROVIDER: S-EPMC7991479 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Computational docking can be used to predict bound conformations and free energies of binding for small-molecule ligands to macromolecular targets. Docking is widely used for the study of biomolecular interactions and mechanisms, and it is applied to structure-based drug design. The methods are fast enough to allow virtual screening of ligand libraries containing tens of thousands of compounds. This protocol covers the docking and virtual screening methods provided by the AutoDock suite of programs, including a basic docking of a drug molecule with an anticancer target, a virtual screen of this target with a small ligand library, docking with selective receptor flexibility, active site prediction and docking with explicit hydration. The entire protocol will require ?5 h.
| S-EPMC4868550 | biostudies-literature
Project description:Virtual compound screening using molecular docking is widely used in the discovery of new lead compounds for drug design. However, this method is not completely reliable and therefore unsatisfactory. In this study, we used massive molecular dynamics simulations of protein-ligand conformations obtained by molecular docking in order to improve the enrichment performance of molecular docking. Our screening approach employed the molecular mechanics/Poisson-Boltzmann and surface area method to estimate the binding free energies. For the top-ranking 1,000 compounds obtained by docking to a target protein, approximately 6,000 molecular dynamics simulations were performed using multiple docking poses in about a week. As a result, the enrichment performance of the top 100 compounds by our approach was improved by 1.6-4.0 times that of the enrichment performance of molecular dockings. This result indicates that the application of molecular dynamics simulations to virtual screening for lead discovery is both effective and practical. However, further optimization of the computational protocols is required for screening various target proteins.
| S-EPMC2746282 | biostudies-literature
Project description:AimTo construct a quantitative pharmacophore model of tubulin inhibitors and to discovery new leads with potent antitumor activities.MethodsLigand-based pharmacophore modeling was used to identify the chemical features responsible for inhibiting tubulin polymerization. A set of 26 training compounds was used to generate hypothetical pharmacophores using the HypoGen algorithm. The structures were further validated using the test set, Fischer randomization method, leave-one-out method and a decoy set, and the best model was chosen to screen the Specs database. Hit compounds were subjected to molecular docking study using a Molecular Operating Environment (MOE) software and to biological evaluation in vitro.ResultsHypo1 was demonstrated to be the best pharmacophore model that exhibited the highest correlation coefficient (0.9582), largest cost difference (70.905) and lowest RMSD value (0.6977). Hypo1 consisted of one hydrogen-bond acceptor, a hydrogen-bond donor, a hydrophobic feature, a ring aromatic feature and three excluded volumes. Hypo1 was validated with four different methods and had a goodness-of-hit score of 0.81. When Hypo1 was used in virtual screening of the Specs database, 952 drug-like compounds were revealed. After docking into the colchicine-binding site of tubulin, 5 drug-like compounds with the required interaction with the critical amino acid residues and the binding free energies < -4 kcal/mol were selected as representative leads. Compounds 1 and 3 exhibited inhibitory activity against MCF-7 human breast cancer cells in vitro.ConclusionHypo1 is a quantitative pharmacophore model for tubulin inhibitors, which not only provides a better understanding of their interaction with tubulin, but also assists in discovering new potential leads with antitumor activities.
| S-EPMC4088285 | biostudies-literature
Project description:Multicanonical molecular dynamics (McMD)-based dynamic docking is a powerful tool to not only predict the native binding configuration between two flexible molecules, but it can also be used to accurately simulate the binding/unbinding pathway. Furthermore, it can also predict alternative binding sites, including allosteric ones, by employing an exhaustive sampling approach. Since McMD-based dynamic docking accurately samples binding/unbinding events, it can thus be used to determine the molecular mechanism of binding between two molecules. We developed the McMD-based dynamic docking methodology based on the powerful, but woefully underutilized McMD algorithm, combined with a toolset to perform the docking and to analyze the results. Here, we showcase three of our recent works, where we have applied McMD-based dynamic docking to advance the field of computational drug design. In the first case, we applied our method to perform an exhaustive search between Hsp90 and one of its inhibitors to successfully predict the native binding configuration in its binding site, as we refined our analysis methods. For our second case, we performed an exhaustive search of two medium-sized ligands and Bcl-xL, which has a cryptic binding site that differs greatly between the apo and holo structures. Finally, we performed a dynamic docking simulation between a membrane-embedded GPCR molecule and a high affinity ligand that binds deep within its receptor's pocket. These advanced simulations showcase the power that the McMD-based dynamic docking method has, and provide a glimpse of the potential our methodology has to unravel and solve the medical and biophysical issues in the modern world.Supplementary informationThe online version contains supplementary material available at 10.1007/s12551-022-01010-z.
| S-EPMC9842809 | biostudies-literature
Project description:Among the different classes of enzymes involved in the ubiquitin pathway, E2 ubiquitin-conjugating enzymes occupy a central role in the ubiquitination cascade. Cdc34-like E2 enzymes are characterized by a 12-14 residue insertion in the proximity of the catalytic site, known as the acidic loop. Cdc34 ubiquitin-charging activity is regulated by CK2-dependent phosphorylation and the regulatory mechanism involves the acidic loop. Indeed, the phosphorylation stabilizes the loop in an open conformation that is competent for ubiquitin charging. Cdc34 is associated with a variety of diseases, such as hepatocellular carcinomas and prostatic adenocarcinomas. In light of its role, the discovery of potential inhibitory compounds would provide the mean to effectively modulate its activity. Here, we carried out a computational study based on molecular dynamics, virtual screening and docking to identify potential inhibitory compounds of Cdc34, modulating the acidic loop conformation. The molecules identified in this study have been designed to act as molecular hinges that can bind the acidic loop in its closed conformation, thus inhibiting the Cdc34-mediated ubiquitination cascade at the ubiquitin-charging step. In particular, we proposed a pharmacophore model featuring two amino groups in the central part of the model and two lateral aromatic chains, which respectively establish electrostatic interactions with the acidic loop (Asp 108 and Glu 109) and a hydrogen bond with Ser 139, which is one of the key residues for Cdc34 activity.
| S-EPMC4050183 | biostudies-literature
Project description:In case of fire in a building, how will people behave in the crowd? The behavior of each individual affects the behavior of others and, conversely, each one behaves considering the crowd as a whole and the individual others. In this article, I propose a three-step method to explore a brand new way to study behavior dynamics. The first step relies on the creation of specific situations with standard techniques (such as mental imagery, text, video, and audio) and an advanced technique [Virtual Reality (VR)] to manipulate experimental settings. The second step concerns the measurement of behavior in one, two, or many individuals focusing on parameters extractions to provide information about the behavior dynamics. Finally, the third step, which uses the parameters collected and measured in the previous two steps in order to simulate possible scenarios to forecast through computational models, understand, and explain behavior dynamics at the social level. An experimental study was also included to demonstrate the three-step method and a possible scenario.
| S-EPMC4635205 | biostudies-literature
Project description:Pseudomonas aeruginosa is an opportunistic Gram-negative bacterium responsible for acute and chronic infections in planktonic state or in biofilms. The sessile structures are known to confer physical stability, increase virulence, and work as a protective armor against antimicrobial compounds. P. aeruginosa can control the expression of genes, population density, and biofilm formation through a process called quorum sensing (QS), a rather complex and hierarchical system of communication. A recent strategy to try and overcome bacterial resistance is to target QS proteins. In this study, a combined multi-level computational approach was applied to find possible inhibitors against P. aeruginosa QS regulator protein MvfR, also known as PqsR, using a database of approved FDA drugs, as a repurposing strategy. Fifteen compounds were identified as highly promising putative MvfR inhibitors. On those 15 MvfR ligand complexes, molecular dynamic simulations and MM/GBSA free-energy calculations were performed to confirm the docking predictions and elucidate on the mode of interaction. Ultimately, the five compounds that presented better binding free energies of association than the reference molecules (a known antagonist, M64 and a natural inducer, 2-nonyl-4-hydroxyquinoline) were highlighted as very promising MvfR inhibitors.
| S-EPMC8868191 | biostudies-literature
Project description:The transient receptor potential vanilloid type 1 (TRPV1) is a heat-activated cation channel protein, which contributes to inflammation, acute and persistent pain. Antagonists of human TRPV1 (hTRPV1) represent a novel therapeutic approach for the treatment of pain. Developing various antagonists of hTRPV1, however, has been hindered by the unavailability of a 3D structure of hTRPV1. Recently, the 3D structures of rat TRPV1 (rTRPV1) in the presence and absence of ligand have been reported as determined by cryo-EM. rTRPV1 shares 85.7% sequence identity with hTRPV1. In the present work, we constructed and reported the 3D homology tetramer model of hTRPV1 based on the cryo-EM structures of rTRPV1. Molecular dynamics (MD) simulations, energy minimizations, and prescreen were applied to select and validate the best model of hTRPV1. The predicted binding pocket of hTRPV1 consists of two adjacent monomers subunits, which were congruent with the experimental rTRPV1 data and the cyro-EM structures of rTRPV1. The detailed interactions between hTRPV1 and its antagonists or agonists were characterized by molecular docking, which helped us to identify the important residues. Conformational changes of hTRPV1 upon antagonist/agonist binding were also explored by MD simulation. The different movements of compounds led to the different conformational changes of monomers in hTRPV1, indicating that TRPV1 works in a concerted way, resembling some other channel proteins such as aquaporins. We observed that the selective filter was open when hTRPV1 bound with an agonist during MD simulation. For the lower gate of hTRPV1, we observed large similarities between hTRPV1 bound with antagonist and with agonist. A five-point pharmacophore model based on several antagonists was established, and the structural model was used to screen in silico for new antagonists for hTRPV1. By using the 3D TRPV1 structural model above, the pilot in silico screening has begun to yield promising hits with activity as hTRPV1 antagonists, several of which showed substantial potency.
| S-EPMC4508124 | biostudies-literature
Project description:The pandemic of novel coronavirus disease 2019 (COVID-19) has rampaged the world, with more than 58.4 million confirmed cases and over 1.38 million deaths across the world by 23 November 2020. There is an urgent need to identify effective drugs and vaccines to fight against the virus. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) belongs to the family of coronaviruses consisting of four structural and 16 non-structural proteins (NSP). Three non-structural proteins, main protease (Mpro), papain-like protease (PLpro), and RNA-dependent RNA polymerase (RdRp), are believed to have a crucial role in replication of the virus. We applied computational ligand-receptor binding modeling and performed comprehensive virtual screening on FDA-approved drugs against these three SARS-CoV-2 proteins using AutoDock Vina, Glide, and rDock. Our computational studies identified six novel ligands as potential inhibitors against SARS-CoV-2, including antiemetics rolapitant and ondansetron for Mpro; labetalol and levomefolic acid for PLpro; and leucal and antifungal natamycin for RdRp. Molecular dynamics simulation confirmed the stability of the ligand-protein complexes. The results of our analysis with some other suggested drugs indicated that chloroquine and hydroxychloroquine had high binding energy (low inhibitory effect) with all three proteins-Mpro, PLpro, and RdRp. In summary, our computational molecular docking approach and virtual screening identified some promising candidate SARS-CoV-2 inhibitors that may be considered for further clinical studies.
| S-EPMC7928605 | biostudies-literature
Project description:Cyclophilin D (CypD) is a peptidyl prolyl isomerase F that resides in the mitochondrial matrix and associates with the inner mitochondrial membrane during the mitochondrial membrane permeability transition. CypD plays a central role in opening the mitochondrial membrane permeability transition pore (mPTP) leading to cell death and has been linked to Alzheimer's disease (AD). Because CypD interacts with amyloid beta (Aβ) to exacerbate mitochondrial and neuronal stress, it is a potential target for drugs to treat AD. Since appropriately designed small organic molecules might bind to CypD and block its interaction with Aβ, 20 trial compounds were designed using known procedures that started with fundamental pyrimidine and sulfonamide scaffolds know to have useful therapeutic effects. Two-dimensional (2D) quantitative structure-activity relationship (QSAR) methods were applied to 40 compounds with known IC50 values. These formed a training set and were followed by a trial set of 20 designed compounds. A correlation analysis was carried out comparing the statistics of the measured IC50 with predicted values for both sets. Selectivity-determining descriptors were interpreted graphically in terms of principle component analyses. These descriptors can be very useful for predicting activity enhancement for lead compounds. A 3D pharmacophore model was also created. Molecular dynamics simulations were carried out for the 20 trial compounds with known IC50 values, and molecular descriptors were determined by 2D QSAR studies using the Lipinski rule-of-five. Fifteen of the 20 molecules satisfied all 5 Lipinski rules, and the remaining 5 satisfied 4 of the 5 Lipinski criteria and nearly satisfied the fifth. Our previous use of 2D QSAR, 3D pharmacophore models, and molecular docking experiments to successfully predict activity indicates that this can be a very powerful technique for screening large numbers of new compounds as active drug candidates. These studies will hopefully provide a basis for efficiently designing and screening large numbers of more potent and selective inhibitors for CypD treatment of AD.
| S-EPMC3985759 | biostudies-literature