1,721,139 research outputs found
The Molecular dYnamics SHAred PharmacophorE (MYSHAPE) approacha new tool to arise docking and pharmacophore modeling performance: virtues and vices
In a recent paper, we presented a new virtual screening workflow that addresses the arising issues of molecular docking and pharmacophore modeling when using a single set of coordinates and a single active ligand [1]. MD simulations were carried out and ligand-protein interactions were analyzed and collected together with their appearance frequency. A pharmacophore model was then created using only the common feature patterns that all the ligands exhibited during MD simulations. This ‘Molecular dYnamics SHAred PharmacophorE’ was then used for virtual screening on active and inactive molecules library. MYSHAPE was also used as constraints for the creation of the docking grid. The application of the MYSHAPE model showed an interesting increase of the screening capability both in terms of sensitivity of the model and specificity when compared to the PDB models. This work [1] was a first essay for a workflow that should be applied to different proteins. In the present study we tried to apply the MYSHAPE approach to other three different ligand-protein systems (ERα; RXRα, and MAPKp38) with the aim to optimize the method to each different biological target taking in consideration the early recognition. The obtained results for these new targets confirmed that it is mandatory, to optimize the virtual screening campaign, the selection of dynamic features and constraints for docking. In particular, the addition of the constraints derived from MD simulation leads to an improvement in the model selectivity for RXRα and ERα in standard precision docking mode. For MAPKp38, validation metrics such as ROC, BEDROC, and AUAC are higher in extra precision mode. For the pharmacophore modeling, the addition of the features derived from the common interactions in MD simulations guarantee an improvement in the AUC for RXRα (37%), and ERα (77%), but light improvement for MAPKp38. MD simulation derived common interactions revealed fundamental for docking selectivity, while they are applied to pharmacophore modeling only when the number of final features in the common and dynamic pharmacophore is higher than the starting static pharmacophore. The strength behind the protocol is the ease of use related to the improvement of results. It also could represent a valid alternative to use very time-consuming techniques such as XP docking with constraints.
Reference: 1. Perricone, U., Wieder, M., Seidel, T., Langer, T., Padova, A., Almerico, A. M., & Tutone, M. (2017). A Molecular Dynamics–Shared Pharmacophore Approach to Boost Early-Enrichment Virtual Screening: A Case Study on Peroxisome Proliferator-Activated Receptor α. ChemMedChem, 2017. DOI: 10.1002/cmdc.20160052
Treatment of Complex Regional Pain Syndrome (Crps): New Perspectives in the Use of Sulfonamides as Modulators of P2x Receptors
Computational Approaches: Drug Discovery and Design in Medicinal Chemistry and Bioinformatics
To date, computational approaches have been recognized as a key component in drug design and discovery workflows. Developed to help researchers save time and reduce costs, several computational tools have been developed and implemented in the last twenty years. At present, they are routinely used to identify a therapeutic target, understand ligand–protein and protein–protein interactions, and identify orthosteric and allosteric binding sites, but their primary use remains the identification of hits through ligand-based and structure-based virtual screening and the optimization of lead compounds, followed by the estimation of the binding free energy. The repurposing of an old drug for the treatment of new diseases, helped by in silico tools, has also gained a prominent role in virtual screening campaigns.
Moreover, the availability and the decreasing cost of hardware and software, together with the development of several web servers on which to upload and download computational data, have contributed to the success of computer-assisted drug design. These improved, accurate, and reliable methods should help to add new and more potent molecules to the paraphernalia of approved drugs. Nevertheless, the ease of access of computational tools in drug design (software, databases, libraries, and web servers) should not encourage users with little or almost no knowledge of the underlying physical basis of the methods used, who could compromise the interpretation of the results. The figure of the computational (medicinal) chemist should be recognized and included in all research groups. These considerations led us to promote a volume collecting some original contributions regarding all aspects of the computational approaches, such as docking, induced-fit docking, molecular dynamics simulations, free energy calculations, and reverse modeling. We also include ligand-based approaches, such as molecular similarity fingerprints, shape methods, pharmacophore modeling, and QSAR. Drug design and the development process strive to predict the metabolic fate of a drug candidate to establish a relationship between the pharmacodynamics and pharmacokinetics and highlight the potential toxicity of the drug candidate. Even though the use of computational approaches is often combined, we tried to identify which of these play a central role in each manuscript.
In this Special Issue, the use of molecular dynamics simulations, both unbiased and biased, cover a major part of the contributions. Non-covalent inhibition of the immunoproteasome was investigated in-depth through MD-binding and binding pose metadynamics [1]. MD simulations provided insight into the structural features of hTSPO (Translocator Protein) and the previously unknown interplay between PK11195, a molecule routinely used in positron emission tomography (PET) for the imaging of neuroinflammatory sites, and cholesterol [2]. The interaction of certain endogen substrates, drug substrates, and inhibitors with wild-type MRP4 (WT-MRP4) and its variants G187W and Y556C were studied to determine differences in the intermolecular interactions and affinity related to SNPs using several approaches, but particularly all-atom, coarse-grained, and umbrella sampling molecular dynamics simulations (AA-MDS and CG-MDS, respectively) [3]. Natural sodium–glucose co-transporter 2 (SGLT2) inhibitors were selected to explore their potential against an emerging uropathogenic bacterial therapeutic target, i.e., FimH, which plays a critical role in the colonization of uropathogenic bacteria on the urinary tract surface, and molecular dynamics simulations were carried out to study the potential interactions [4]. Doxorubicin encapsulation in carbon nanotubes with haeckelite or Stone–Wales defects as drug carriers were investigated using a molecular dynamics approach [5]. The combined use of different approaches has been reported in a series of papers associated with the virtual screening of libraries. Almeelebia and co. screened 224,205 natural compounds from the ZINC database against the catalytic site of the Mtb proteasome [6]. Pharmacophore-based virtual screening and molecular docking were carried out to identify potential Src inhibitors starting from a total of 891 molecules. Finally, MD simulations identified two molecules as potential lead compounds against Src kinase [7]. An in silico study identified a methotrexate analog as a potential inhibitor of drug-resistant human dihydrofolate reductase for cancer therapeutics [8]. A structure-based method for high-throughput virtual screening aimed to identify new dual-target hit molecules for acetylcholinesterase, and the α7 nicotinic acetylcholine receptor was reported and confirmed in vitro [9]. A new complementary computational analysis called “dock binning” evaluates antibody–antigen docking models to identify why and where they might compete in terms of possible binding sites on the antigen [10]. Interesting drug repurposing strategies have been reported. Hudson and Samudrala presented a computational analysis of a novel drug opportunities (CANDO) platform for shotgun multitarget repurposing. It implements several pipelines for the large-scale modeling and simulation of interactions between comprehensive libraries of drugs/compounds and protein structures [11]. Qi and co. data-mined the crowd extracted expression of differential signatures (CREEDS) database to evaluate the similarities between gene expression signature (GES) profiles from drugs and their indicated diseases for GES-guided drug-repositioning approaches [12]. In late 2019, the SARS-CoV-2 pandemic focused the attention of many researchers intending to find not only vaccines but also new antiviral drugs. These reasons boosted the use of computational approaches to explore large libraries of natural compounds, already approved drugs, and in-house and commercial compounds [13,14]. In this issue, Baig and co. studied the efficacy of the Mpro inhibitor PF-00835231 against Mpro and its reported mutants in clinical trials. Several in silico approaches were used to investigate and compare the efficacy of PF-00835231 and five drugs previously documented to inhibit Mpro [15]. Li and co. computationally investigated the MPD3 phytochemical database along with the pool of reported natural antiviral compounds to be used against SARS-CoV-2 [16]. Pedretti and co., exploiting the availability of resolved structures, designed a structure-based computational approach. The innovative idea of their study was to exploit known inhibitors of SARS-CoV 3CL-Pro as a training set to perform and validate multiple virtual screening campaigns [17]. In the context of antiviral drugs, Regad and co. investigated the emergence of HIV-2 resistance. They proposed a structural analysis of 31 drug-resistant mutants of HIV-2 protease (PR2), an important target against HIV-2 infection [18]. A wide series of contributions regarding the use of QSAR, machine learning, and deep learning has reported interesting outcomes. A multiple-molecule drug design based on systems biology approaches and a deep neural network to mitigate human skin aging was developed by Yeh and co. With the proposed systems medicine design procedure, they not only shed light on the skin-aging molecular progression mechanisms, but they also suggested two multiple-molecule drugs to mitigate human skin aging [19]. The construction of quantitative structure–activity relationship (QSAR) models was used to predict the biological activities of the compounds obtained with virtual screening and identify new selective chemical entities for the COX-2 enzyme [20]. The three-dimensional QSAR model, employing a common-features pharmacophore as an alignment rule, was built on 20 highly active/selective HDAC1 inhibitors. The predictive power of the 3D QSAR model represents a useful filtering tool for screening large chemical databases, finding novel derivatives with improved HDAC1 inhibitory activity [21]. Different machine learning (ML) and deep learning (DL) algorithms using various integer and binary type fingerprints were evaluated to develop quantitative structure–activity relationship (QSAR) models, which are important for hERG potassium channel blocker prediction [22].
Throughout this Special Issue, all the recent aspects of the computational approaches applied to several research fields are reported. We express our deep gratitude to all the contributors to this Special Issue for their commitment, hard work, and outstanding papers. We also thank all the reviewers involved in the manuscript revisions for their unpaid contributions to improve any aspects of the submitted works. Last but not least, we deeply thank Mrs. Jessie Zhang for her assistance during the period in which we served as guest editors
Natural products as non-covalent and covalent modulators of the KEAP1/NRF2 pathway exerting antioxidant effects
By controlling several antioxidant and detoxifying genes at the transcriptional level, including NAD(P)H quinone oxidoreductase 1 (NQO1), multidrug resistance-associated proteins (MRPs), UDP-glucuronosyltransferase (UGT), glutamate-cysteine ligase catalytic (GCLC) and modifier (GCLM) subunits, glutathione S-transferase (GST), sulfiredoxin1 (SRXN1), and heme-oxygenase-1 (HMOX1), the KEAP1/NRF2 pathway plays a crucial role in the oxidative stress response. Accordingly, the discovery of modulators of this pathway, activating cellular signaling through NRF2, and targeting the antioxidant response element (ARE) genes is pivotal for the development of effective antioxidant agents. In this context, natural products could represent promising drug candidates for supplementation to provide antioxidant capacity to human cells. In recent decades, by coupling in silico and experimental methods, several natural products have been characterized to exert antioxidant effects by targeting the KEAP1/NRF2 pathway. In this review article, we analyze several natural products that were investigated experimentally and in silico for their ability to modulate KEAP1/NRF2 by non-covalent and covalent mechanisms. These latter represent the two main sections of this article. For each class of inhibitors, we reviewed their antioxidant effects and potential therapeutic applications, and where possible, we analyzed the structure-activity relationship (SAR). Moreover, the main computational techniques used for the most promising identified compounds are detailed in this survey, providing an updated view on the development of natural products as antioxidant agents
Exploring the non-covalent ligand-binding mechanism on immunoproteasome by enhanced Molecular Dynamics
Selective inhibition of immunoproteasome is a valuable strategy to treat autoimmune and
inflammatory diseases, and hematologic malignancies. In particular, non-covalent inhibition
is strongly desirable because it is free of the drawbacks and side effects associated with
covalent inhibition. Recently, a new series of amide derivatives with Ki values in the low/submicromolar ranges toward the β1i subunit have been identified as non-covalent inhibitors
1
. We
investigated the binding mechanism of the most potent and selective inhibitor (1) to elucidate
the steps from the ligand entrance into the binding pocket to the ligand-induced
conformational changes. We carried out a total of 400ns of MD-binding analysis, followed by
200ns of plain MD. The trajectories clustering allowed identifying three representative poses
evidencing new key interactions with Phe31 and Lys33 together to a flipped orientation of a
representative pose. Further, Binding pose metadynamics (BPMD) studies have been performed
to evaluate the binding affinity, comparing
(1) with other four inhibitors of β1i subunit (2, 3, 4, and 5). Results are consistent with
experimental values of inhibition, confirming (1) as a lead compound of this series. The
adopted methods provided a full dynamic description of the binding events and the information
obtained could be exploited for the rational design of new and more active inhibitor
Deciphering the Nonsense Readthrough Mechanism of Action of Ataluren: An in Silico Compared Study
Ataluren was reported to suppress nonsense mutations by promoting the readthrough of premature stop codons, although its mechanism of action (MOA) is still debated. The likely interaction of Ataluren with CFTR-mRNA has been previously studied by molecular dynamics. In this work we extended the modeling of Ataluren's MOA by complementary computational approaches such as induced fit docking (IFD), quantum polarized ligand docking (QPLD), MM-GBSA free-energy calculations, and computational mutagenesis. In addition to CFTR-mRNA, this study considered other model targets implicated in the translation process, such as eukaryotic rRNA 18S, prokaryotic rRNA 16S, and eukaryotic Release Factor 1 (eRF1), and we performed a comparison with a new promising Ataluren analogue (NV2445) and with a series of aminoglycosides, known to suppress the normal proofreading function of the ribosome. Results confirmed mRNA as the most likely candidate target for Ataluren and its analogue, and binding energies calculated after computational mutagenesis highlighted how Ataluren's interaction with the premature stop codon could be affected by ancillary nucleotides in the genetic context
Sulfonamide moiety as "molecular chimera" in the design of new drugs
The -SO2NH- group is of great significance in modern pharmaceutical use since in sulfa-drugs it is possible to introduce easily chemical modifications, and even small changes may lead to an improved version of an already existing drug
Reverse screening on indicaxanthin from Opuntia ficus-indica as natural chemoactive and chemopreventive agent
Indicaxanthin is a bioactive and bioavailable betalain pigment extracted from Opuntia ficus indica fruits. Indicaxanthin has pharmacokinetic proprieties, rarely found in other phytochemicals, and it has been demonstrated that it provides a broad-spectrum of pharmaceutical activity, exerting anti-proliferative, anti-inflammatory, and neuromodulator effects. The discovery of the Indicaxanthin physiological targets plays an important role in understanding the biochemical mechanism. In this study, combined reverse pharmacophore mapping, reverse docking, and text-based database search identified Inositol Trisphosphate 3-Kinase (ITP3K-A), Glutamate carboxypeptidase II (GCPII), Leukotriene-A4 hydrolase (LTA4H), Phosphoserine phosphatase (HPSP), Phosphodiesterase 4D (PDE4D), AMPA receptor (GluA3 and GluA2 subunits) and Kainate receptor (GluK1 isoform) as potential targets for Indicaxanthin. These targets are implicated in neuromodulation, and inflammatory regulation, normally expressed mostly in the CNS, and expressed (or overexpressed) in cancer tissues (i.e. breast, thyroid, and prostate cancer cells). Moreover, this study provides qualitative and quantitative information about dynamic interactions of Indicaxanthin at the binding site of target proteins, through molecular dynamics simulations and MM-GBSA
Indicaxanthin, a multi-target natural compound from Opuntia ficus-indica fruit: From its poly-pharmacological effects to biochemical mechanisms and molecular modelling studies
Over the latest years phytochemical consumption has been associated to a decreased risk of both theonset and the development of a number of pathological conditions. In this context indicaxanthin, abetalain pigment fromOpuntiaficus-indicafruit, has been the object of sound research. Explored, atfirst,for its mere antioxidant potential, Indicaxanthin is now regarded as a redox-active compound able toexert significant poly-pharmacological effects against several targets in a number of experimental con-ditions bothin vivoandin vitro. This paper aims to provide an overview on the therapeutical effects ofindicaxanthin, ranging from the anti-inflammatory to the neuro-modulatory and anti-tumoral ones andfavored by its high bioavailability. Moreover, biochemical and molecular modelling investigations areaimed to identify the pharmacological targets the compound is able to interact with and to address thechallenging development in the future researc
Immunoproteasome and Non-Covalent Inhibition: Exploration by Advanced Molecular Dynamics and Docking Methods
The selective inhibition of immunoproteasome is a valuable strategy to treat autoimmune, inflammatory diseases, and hematologic malignancies. Recently, a new series of amide derivatives as non-covalent inhibitors of the β1i subunit with Ki values in the low/submicromolar ranges have been identified. Here, we investigated the binding mechanism of the most potent and selective inhibitor, N-benzyl-2-(2-oxopyridin-1(2H)-yl)propanamide (1), to elucidate the steps from the ligand entrance into the binding pocket to the ligand-induced conformational changes. We carried out a total of 400 ns of MD-binding analyses, followed by 200 ns of plain MD. The trajectories clustering allowed identifying three representative poses evidencing new key interactions with Phe31 and Lys33 together in a flipped orientation of a representative pose. Further, Binding Pose MetaDynamics (BPMD) studies were performed to evaluate the binding stability, comparing 1 with four other inhibitors of the β1i subunit: N-benzyl-2-(2-oxopyridin-1(2H)-yl)acetamide (2), N-cyclohexyl-3-(2-oxopyridin-1(2H)-yl)propenamide (3), N-butyl-3-(2-oxopyridin-1(2H)-yl)propanamide (4), and (S)-2-(2-oxopyridin-1(2H)-yl)-N,4-diphenylbutanamide (5). The obtained results in terms of free binding energy were consistent with the experimental values of inhibition, confirming 1 as a lead compound of this series. The adopted methods provided a full dynamic description of the binding events, and the information obtained could be exploited for the rational design of new and more active inhibitors
- …
