Unbinding ligands from SARS-CoV-2 Mpro via umbrella sampling simulations
Ontology highlight
ABSTRACT:
SUBMITTER: Tam N
PROVIDER: S-EPMC8790385 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The main protease (3CLpro) is one of the essential components of the SARS-CoVs viral life cycle, which makes it an interesting target for overpowering these viruses. Although many covalent and noncovalent inhibitors have been designed to inhibit this molecular target, none have gained FDA approval as a drug. Because of the high rate of COVID-19 pandemic development, in addition to laboratory research, we require in silico methods to accelerate rational drug design. The unbinding pathways of two SARS-CoV and SARS-CoV-2 3CLpro noncovalent inhibitors with the PDB IDs: 3V3M, 4MDS, 6W63, 5RF7 were explored from a comparative perspective using unbiased molecular dynamics (UMD) simulations. We uncovered common weak points for selected inhibitors that could not interact significantly with a binding pocket at specific residues by all their fragments. So water molecules entered the free binding S regions and weakened protein-inhibitor fundamental interactions gradually. N142, G143, and H163 are the essential residues, which cause key protein-ligand interactions in the binding pocket. We believe that these results will help design new potent inhibitors against SARS-CoV-2.
| S-EPMC8827461 | biostudies-literature
Project description:To date, the COVID-19 pandemic has still been infectious around the world, continuously causing social and economic damage on a global scale. One of the most important therapeutic targets for the treatment of COVID-19 is the main protease (Mpro) of SARS-CoV-2. In this study, we combined machine-learning (ML) model with atomistic simulations to computationally search for highly promising SARS-CoV-2 Mpro inhibitors from the representative natural compounds of the National Cancer Institute (NCI) Database. First, the trained ML model was used to scan the library quickly and reliably for possible Mpro inhibitors. The ML output was then confirmed using atomistic simulations integrating molecular docking and molecular dynamic simulations with the linear interaction energy scheme. The results turned out to show that there was evidently good agreement between ML and atomistic simulations. Ten substances were proposed to be able to inhibit SARS-CoV-2 Mpro. Seven of them have high-nanomolar affinity and are very potential inhibitors. The strategy has been proven to be reliable and appropriate for fast prediction of SARS-CoV-2 Mpro inhibitors, benefiting for new emerging SARS-CoV-2 variants in the future accordingly.
| S-EPMC9950021 | biostudies-literature
Project description:Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a dreaded pandemic in lack of specific therapeutic agent. SARS-CoV-2 Mpro, an essential factor in viral pathogenesis, is recognized as a prospective therapeutic target in drug discovery against SARS-CoV-2. To tackle this pandemic, Food and Drug Administration-approved drugs are being screened against SARS-CoV-2 Mpro via in silico and in vitro methods to detect the best conceivable drug candidates. However, identification of natural compounds with anti-SARS-CoV-2 Mpro potential have been recommended as rapid and effective alternative for anti-SARS-CoV-2 therapeutic development. Thereof, a total of 653 natural compounds were identified against SARS-CoV-2 Mpro from NP-lib database at MTi-OpenScreen webserver using virtual screening approach. Subsequently, top four potential compounds, i.e. 2,3-Dihydroamentoflavone (ZINC000043552589), Podocarpusflavon-B (ZINC000003594862), Rutin (ZINC000003947429) and Quercimeritrin 6"-O-L-arabinopyranoside (ZINC000070691536), and co-crystallized N3 inhibitor as reference ligand were considered for stringent molecular docking after geometry optimization by DFT method. Each compound exhibited substantial docking energy >-12 kcal/mol and molecular contacts with essential residues, including catalytic dyad (His41 and Cys145) and substrate binding residues, in the active pocket of SARS-CoV-2 Mpro against N3 inhibitor. The screened compounds were further scrutinized via absorption, distribution, metabolism, and excretion - toxicity (ADMET), quantum chemical calculations, combinatorial molecular simulations and hybrid QM/MM approaches. Convincingly, collected results support the potent compounds for druglikeness and strong binding affinity with the catalytic pocket of SARS-CoV-2 Mpro. Hence, selected compounds are advocated as potential inhibitors of SARS-CoV-2 Mpro and can be utilized in drug development against SARS-CoV-2 infection.
| S-EPMC7929395 | biostudies-literature
Project description:We compute the absolute binding affinities for two ligands bound to the FKBP protein using nonequilibrium unbinding simulations. The methodology is straightforward requiring little or no modification to many modern molecular simulation packages. The approach makes use of a physical pathway, eliminating the need for complicated alchemical decoupling schemes. We compare our nonequilibrium results to those obtained via a fully equilibrium approach and to experiment. The results of this study suggest that to obtain accurate results using nonequilibrium approaches one should use the stiff-spring approximation with the second cumulant expansion. From this study we conclude that nonequilibrium simulation could provide a simple means to estimate protein-ligand binding affinities.
| S-EPMC2905451 | biostudies-literature
Project description:IntroductionRiboswitches are cis-acting regulatory RNA elements prevalently located in the leader sequences of bacterial mRNA. An adenine sensing riboswitch cis-regulates adeninosine deaminase gene (add) in Vibrio vulnificus. The structural mechanism regulating its conformational changes upon ligand binding mostly remains to be elucidated. In this open framework it has been suggested that the ligand stabilizes the interaction of the distal "kissing loop" complex. Using accurate full-atom molecular dynamics with explicit solvent in combination with enhanced sampling techniques and advanced analysis methods it could be possible to provide a more detailed perspective on the formation of these tertiary contacts.MethodsIn this work, we used umbrella sampling simulations to study the thermodynamics of the kissing loop complex in the presence and in the absence of the cognate ligand. We enforced the breaking/formation of the loop-loop interaction restraining the distance between the two loops. We also assessed the convergence of the results by using two alternative initialization protocols. A structural analysis was performed using a novel approach to analyze base contacts.ResultsContacts between the two loops were progressively lost when larger inter-loop distances were enforced. Inter-loop Watson-Crick contacts survived at larger separation when compared with non-canonical pairing and stacking interactions. Intra-loop stacking contacts remained formed upon loop undocking. Our simulations qualitatively indicated that the ligand could stabilize the kissing loop complex. We also compared with previously published simulation studies.Discussion and conclusionsKissing complex stabilization given by the ligand was compatible with available experimental data. However, the dependence of its value on the initialization protocol of the umbrella sampling simulations posed some questions on the quantitative interpretation of the results and called for better converged enhanced sampling simulations.
| S-EPMC4464220 | biostudies-literature
Project description:The first oral drug for the treatment of COVID-19, Paxlovid, has been authorized; however, nirmatrelvir, a major component of the drug, is reported to be associated with some side effects. Moreover, the appearance of many novel variants raises concerns about drug resistance, and designing new potent inhibitors to prevent viral replication is thus urgent. In this context, using a hybrid approach combining machine learning (ML) and free energy simulations, 6 compounds obtained by modifying nirmatrelvir were proposed to bind strongly to SARS-CoV-2 Mpro. The structural modification of nirmatrelvir significantly enhances the electrostatic interaction free energy between the protein and ligand and slightly decreases the vdW term. However, the vdW term is the most important factor in controlling the ligand-binding affinity. In addition, the modified nirmatrelvir might be less toxic to the human body than the original inhibitor.
| S-EPMC10233213 | biostudies-literature
Project description:Inhibition of the main protease (Mpro) of SARS-CoV-2 has been suggested to be vital in shutting down viral replication in a host. Most efforts aimed at inhibiting MPro activity have been channeled into competitive inhibition at the active site, but this strategy will require a high inhibitor concentration and impressive inhibitor-MPro binding affinity. Allosteric inhibition can potentially serve as an effective strategy for alleviating these limitations. In this study, the ability of antiviral natural products to inhibit MPro in an allosteric fashion was explored with in silico techniques. Molecular docking revealed a strong interaction between casticin, an antiviral flavonoid, and Mpro at a site distant from the active site. This site, characterized as a distal site, has been shown to have an interdependent dynamic effect with the active site region. Mpro only, Mpro-peptide (binary) and Mpro-peptide-casticin (ternary) complexes were subjected to molecular dynamics simulations for 50 ns to investigate the modulatory activity of casticin binding on Mpro. Molecular dynamic simulations revealed that binding of casticin at the distal site interferes with the proper orientation of the peptide substrate in the oxyanion hole of the active site, and this could lead to a halt or decrease in catalytic activity. This study therefore highlights casticin as a potential allosteric modulator of the SARS-CoV-2 main protease, which could be optimized and developed into a potential lead compound for anti-SARS-CoV-2 chemotherapy. Graphical Abstract Supplementary Information The online version contains supplementary material available at 10.1007/s42250-022-00411-7.
| S-EPMC9261893 | biostudies-literature
Project description:Drug discovery is expensive and high-risk. Its main reasons of failure are lack of efficacy and toxicity of a drug candidate. Binding affinity for the biological target has been usually considered one of the most relevant figures of merit to judge a drug candidate along with bioavailability, selectivity and metabolic properties, which could depend on off-target interactions. Nevertheless, affinity does not always satisfactorily correlate with in vivo drug efficacy. It is indeed becoming increasingly evident that the time a drug spends in contact with its target (aka residence time) can be a more reliable figure of merit. Experimental kinetic measurements are operatively limited by the cost and the time needed to synthesize compounds to be tested, to express and purify the target, and to setup the assays. We present here a simple and efficient molecular-dynamics-based computational approach to prioritize compounds according to their residence time. We devised a multiple-replica scaled molecular dynamics protocol with suitably defined harmonic restraints to accelerate the unbinding events while preserving the native fold. Ligands are ranked according to the mean observed scaled unbinding time. The approach, trivially parallel and easily implementable, was validated against experimental information available on biological systems of pharmacological relevance.
| S-EPMC4477625 | biostudies-literature
Project description:Corona Virus Disease 2019 (COVID-19) caused by Severe Acute Respiratory Syndrome coronavirus (SARS CoV-2) has been declared a worldwide pandemic by WHO recently. The complete understanding of the complex genomic structure of SARS CoV-2 has enabled the use of computational tools in search of SARS CoV-2 inhibitors against the multiple proteins responsible for its entry and multiplication in human cells. With this endeavor, 177 natural, anti-viral chemical entities and their derivatives, selected through the critical analysis of the literatures, were studied using pharmacophore screening followed by molecular docking against RNA dependent RNA polymerase and main protease. The identified hits have been subjected to molecular dynamic simulations to study the stability of ligand-protein complexes followed by ADMET analysis and Lipinski filters to confirm their drug likeliness. It has led to an important start point in the drug discovery and development of therapeutic agents against SARS CoV-2.
| S-EPMC8059878 | biostudies-literature
Project description:The SARS-CoV-2 virus has been identified as a causative agent for COVID-19 pandemic. About more than 6.3 million fatalities have been attributed to COVID-19 worldwide to date. Finding a viable cure for the illness is urgently needed in light of the present pandemic. The prominence of main protease in the life cycle of virus shapes the main protease as a viable target for design and development of antiviral agents to combat COVID-19. The current study presents the fragment linking strategy to design the novel Mpro inhibitors for COVID-19. A total of 293,451 fragments from diversified libraries have been screened for their binding affinity towards Mpro enzyme. The best 1600 fragment hits were subjected to fragment joining to achieve 100 new molecules using Schrödinger software. The resulting molecules were further screened for their Mpro binding affinity, ADMET, and drug-likeness features. The best 13 molecules were selected, and the first 6 compounds were investigated for their ligand-receptor complex stability through a molecular dynamics study using GROMACS software. The resulting molecules have the potential to be further evaluated for COVID-19 drug discovery.
| S-EPMC9399563 | biostudies-literature