The main proteinase (Mpro) of SARS-CoV-2 plays a critical role in cleaving viral polyproteins into functional proteins required for viral replication and assembly, making it a prime drug target for COVID-19. It is well known that noncompetitive inhibition offers potential therapeutic options for treating COVID-19, which can effectively reduce the likelihood of cross-reactivity with other proteins and increase the selectivity of the drug. Therefore, the discovery of allosteric sites of Mpro has both scientific and practical significance. In this study, we explored the binding characteristics and inhibiting process of Mpro activity by two recently reported allosteric inhibitors, pelitinib and AT7519 which were obtained by the X-ray screening experiments, to probe the allosteric mechanism via molecular dynamic (MD) simulations. We found that pelitinib and AT7519 can stably bind to Mpro far from the active site. The binding affinity is estimated to be -24.37 ± 4.14 and - 26.96 ± 4.05 kcal/mol for pelitinib and AT7519, respectively, which is considerably stable compared with orthosteric drugs. Furthermore, the strong binding caused clear changes in the catalytic site of Mpro, thus decreasing the substrate accessibility. The community network analysis also validated that pelitinib and AT7519 strengthened intra- and inter-domain communication of Mpro dimer, resulting in a rigid Mpro, which could negatively impact substrate binding. In summary, our findings provide the detailed working mechanism for the two experimentally observed allosteric sites of Mpro. These allosteric sites greatly enhance the \'druggability\' of Mpro and represent attractive targets for the development of new Mpro inhibitors.

Protein phosphorylation plays an important role in the signal transduction and is capable of regulation of cell activity. The death-associated protein kinase 1 (DAPK1), as a Ser/Thr kinase, interacts with calmodulin (CaM) to regulate apoptotic and autophagic signaling. Autophosphorylation of DAPK1 at Ser308 located at the autoregulatory domain (ARD) blocks CaM binding and inhibits kinase catalytic activity. However, the mechanism underlying the influence of Ser308 phosphorylation (pS308) on the DAPK1 activity remains unclear. Here, we performed multiple, microsecond length molecular dynamics (MD) simulations, the molecular mechanics generalized Born/surface area (MM-GBSA) binding free energy calculations, principal component analysis, and dynamic cross-correlation analysis to unravel the conformational dynamics and allostery of the DAPK1 - CaM interaction triggered by the pS308 at the ARD. MD simulations showed that pS308 affected the conformational stability of the DAPK1 - CaM complex. Further energetic and structural exploration revealed that pS308 weakened the association of the phosphorylated DAPK1 to CaM, which lowered the susceptibility of DAPK1 to be activated by CaM. This result can provide mechanistic insights into the molecular underpinning through which the DAPK1 kinase activity is modulated by the auto-phosphorylation.Communicated by Ramaswamy H. Sarma.

Major facilitator superfamily domain-containing 2a (Mfsd2a) is a sodium-dependent lysophosphatidylcholine cotransporter that plays an important role in maintaining the integrity of the blood-brain barrier and neurological function. Abnormal degradation of Mfsd2a often leads to dysfunction of the blood-brain barrier, while upregulation of Mfsd2a can retrieve neurological damage. It has been reported that Mfsd2a can be specifically recognized and ubiquitinated by neural precursor cell-expressed developmentally downregulated gene 4 type 2 (NEDD4-2) ubiquitin ligase and finally degraded through the proteasome pathway. However, the structural basis for the specific binding of Mfsd2a to NEDD4-2 is unclear. In this work, we combined deep learning and molecular dynamics simulations to obtain a Mfsd2a structure with high quality and a stable Mfsd2a/NEDD4-2-WW3 interaction model. Moreover, molecular mechanics generalized Born surface area (MM-GBSA) methods coupled with per-residue energy decomposition studies were carried out to analyze the key residues that dominate the binding interaction. Based on these results, we designed three peptides containing the key residues by truncating the Mfsd2a sequences. One of them was found to significantly inhibit Mfsd2a ubiquitination, which was further validated in an oxygen-glucose deprivation (OGD) model in a human microvascular endothelial cell line. This work provides some new insights into the understanding of Mfsd2a and NEDD4-2 interaction and might promote further development of drugs targeting Mfsd2a ubiquitination.

Succinate dehydrogenase inhibitors (SDHIs) are a class of fungicides targeting the pathogenic fungi mitochondrial SDH. Here, molecular docking, three-dimensional quantitative structure-activity relationship (3D-QSAR), and molecular dynamics (MD) simulations were used to guide SDHI innovation. Molecular docking was performed to explore the binding modes of SDH and its inhibitors. 3D-QSAR models were carried out on 33 compounds with activity against Rhizoctonia cerealis (R. cerealis); their structure-activity relationships were analyzed using comparative molecular field analysis and comparative molecular similarity indices analysis. MD simulations were used to assess the stability of the complexes under physiological conditions, and the results were consistent with molecular docking. Binding free energy was calculated through the molecular mechanics generalized born surface area method, and the binding free energy was decomposed. The results are consistent with the activity of bioassay and indicate that van der Waals and lipophilic interactions contribute the most in the molecular binding process. Afterward, we designed and synthesized 12 compounds under the guidance of the above-mentioned analyses, bioassay found that F9 was active against R. cerealis with the EC50 value of 9.43 μg/mL, and F4, F5, and F9 were active against Botrytis cinerea with an EC50 values of 5.80, 3.17, and 1.63 μg/mL, respectively. They all showed good activity between positive controls of pydiflumetofen and thifluzamide. Our study provides new considerations for effective SDHIs discovery.

BACKGROUND: ZAP-70 (zeta-chain-associated protein of 70 kDa), serving as a critical regulator for T cell antigen receptor signaling, represents an attractive therapeutic target for autoimmunity disease. How the mechanistical mechanism of ZAP-70 to a human autoimmune syndrome-associated R192W mutation remains unclear. The results indicated that the R192W mutation of ZAP-70 clearly affected the conformational flexibility of the N-terminal ITAM-Y2P. Structural analysis unveiled that the R192W mutation of ZAP-70 caused the exposure of the N-terminal ITAM-Y2P to the solvent. MM-GBSA binding free energy calculations exhibited that the R192W mutation decreased the binding affinity of ITAM-Y2P to the ZAP-70 mutant. Residue-based free energy decomposition further revealed that the protein-peptide interaction networks involving electrostatic interactions provide significant contributions for complex formation. The energy unfavorable residues include Arg43, Arg192, Tyr240, and Lys244 from ZAP-70 and Asn301, Leu303, pY304, and pY315 from ITAM-Y2P in the R192W mutant. Our obtained results may help the understanding of the deactivation mechanism of ZAP-70 induced by the R192W mutation.
METHODS: In the work, multiple replica molecular dynamics simulations and molecular mechanics-generalized Born surface area (MM-GBSA) method were performed to reveal the doubly phosphorylated ITAMs (ITAM-Y2P)-mediated deactivation mechanism of ZAP-70 induced by the R192W mutation.

Most of the breast cancers are estrogen receptor-positive recurring with a steady rate of up to 20 years dysregulating the normal cell cycle. Dinaciclib is still in clinical trials and considered as a research drug against such cancers targeting CDK2.The major goal of this study was to identify the potential inhibitors of CDK-2 present in Moringa oleifera for treating hormonal receptor positive breast cancers. For this purpose, in silico techniques; molecular docking, MM-GBSA and molecular dynamics simulations were employed to screen Moringa oleifera compounds and their anticancer potential was determined against CDK-2 protein targets. Among 36 compounds of Moringa oleifera reported in literature, chlorogenic acid (1), quercetin (2), ellagic acid (3), niazirin (4), and kaempferol (5) showed good affinity with the target. The interaction of the compounds was visualized using PYMOL software. The profiles of absorption, distribution, metabolism, excretion (ADME) and toxicity were determined using SWISS and ProTox II webservers. The MTT assay was performed in-vitro using MCF-7 cancer cell lines to validate the anticancer potential of Moringa oleifera leaf extract.MTT assay results revealed no significant change in proliferation of Mcf-7 cells following 24 h treatment with fraction A (petroleum ether). However, significant antiproliferative effect was observed at 200 µg/mL dose of fraction B (ethyl acetate) and cell viability was reduced to 40%.In conclusion, the data suggested that all the compounds with highest negative docking score than the reference could be the potential candidates for cyclin dependent kinase-2 (CDK-2) inhibition while ellagic acid, chlorogenic acid and quercetin being the most stable and potent inhibitors to treat estrogen receptor positive breast cancer targeting CDK-2. Moreover, the data suggested that further investigation is required to determine the optimum dose for significant antiproliferative effects using in-vivo models to validate our findings of in-silico analysis.

Computer-aided molecular modeling was applied to design a series of Spodoptera frugiperda RyR agonists. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to generate 3D-QSAR models. MD simulations in the complex with S. frugiperda native, mutant RyR, and mammalian RyR1 under physiological conditions were used to validate the detailed binding mechanism. Binding free energy calculation by molecular mechanics generalized surface area (MM-GBSA) explained the role of key amino acid residues in ligand-receptor binding. Therefore, 14 new compounds were effectively designed and synthesized, and a bioassay indicated that compounds A-2 and A-3 showed comparable activity to that of chloranthraniliprole with LC50 values of 0.27, 0.18, and 0.20 mg L-1, respectively, against S. frugiperda. Most target compounds also displayed good activity against Mythinma separata at 0.1 mg L-1. Molecular docking and MM-GBSA calculations demonstrated that A-3 had a better binding capacity with native and mutant S. frugiperda RyRs.

The heat shock protein (HSP90) has been an import target of drug design in the treatment of human disease. An exploration of the conformational changes in HSP90 can provide useful information for the development of efficient inhibitors targeting HSP90. In this work, multiple independent all-atom molecular dynamics (AAMD) simulations followed by calculations of the molecular mechanics generalized Born surface area (MM-GBSA) were performed to explore the binding mechanism of three inhibitors (W8Y, W8V, and W8S) to HSP90. The dynamics analyses verified that the presence of inhibitors impacts the structural flexibility, correlated movements, and dynamics behavior of HSP90. The results of the MM-GBSA calculations suggest that the selection of GB models and empirical parameters has important influences on the predicted results and verify that van der Waals interactions are the main forces that determine inhibitor-HSP90 binding. The contributions of separate residues to the inhibitor-HSP90 binding process indicate that hydrogen-bonding interactions (HBIs) and hydrophobic interactions play important roles in HSP90-inhibitor identifications. Moreover, residues L34, N37, D40, A41, D79, I82, G83, M84, F124, and T171 are recognized as hot spots of inhibitor-HSP90 binding and provide significant target sites of for the design of drugs related to HSP90. This study aims to contribute to the development of efficient inhibitors that target HSP90 by providing an energy-based and theoretical foundation.

β-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward the treatment of Alzheimer\'s disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively determine the identification mechanism of BACE1 for three inhibitors, 60W, 954 and 60X. The analyses of MD trajectories indicated that the presence of three inhibitors influences the structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor-BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor-BACE1 binding, which provides a direction for future drug design toward the treatment of AD.

BACKGROUND: In the replication of SARS-CoV-2, the main protease (Mpro/3CLpro) is significant. It is conserved in a number of novel coronavirus variations, and no known human proteases share its cleavage sites. Therefore, 3CLpro is an ideal target. In the report, we screened five potential inhibitors (1543, 2308, 3717, 5606, and 9000) of SARS-CoV-2 Mpro through a workflow. The calculation of MM-GBSA binding free energy showed that three of the five potential inhibitors (1543, 2308, 5606) had similar inhibitor effects to X77 against Mpro of SARS-CoV-2. In conclusion, the manuscript lays the groundwork for the design of Mpro inhibitors.
METHODS: In the virtual screening phase, we used structure-based virtual screening (Qvina2.1) and ligand-based virtual screening (AncPhore). In the molecular dynamic simulation part, we used the Amber14SB + GAFF force field to perform molecular dynamic simulation of the complex for 100 ns (Gromacs2021.5) and performed MM-GBSA binding free energy calculation according to the simulation trajectory.