基于分子对接的thanatin与NDM-1分子动力学模拟研究

耐碳青霉烯类抗生素的超级细菌给人类健康带来了严重威胁,其所携带的金属 β-内酰胺酶编码基因是耐药性的主要来源。NDM-1作为其中传播最广、活性最强的 β-内酰胺酶,其抑制剂的研发刻不容缓。具有广谱作用的抗菌肽thanatin对NDM-1展现出了较好的抑制效果,但抑制机理并不清楚。本文使用HPEPDOCK与Rosetta FlexPepDock服务器,将thanatin与NDM-1进行了分子对接,并使用Desmond软件包对对接模型进行了分子动力学模拟。结果表明,thanatin与NDM-1活性中心的Zn2+ 并无直接相互作用,而作为竞争性抑制剂结合于NDM-1的活性口袋,阻止抗生素分子进入活性口袋与Zn2+ 结合,从而抑制NDM-1的水解活性。本文为研发有效的NDM临床抑制剂探索了可行的方法。     

Molecular docking studies of a group of hydroxamate inhibitors with gelatinase-A by molecular dynamics

We have performed docking and molecular dynamics simulations of hydroxamates complexed with human gelatinase-A (MMP-2) to gain insight into the structural and energetic preferences of these inhibitors. The study was conducted on a selected set of eleven compounds with variation in structure and activity. Molecular dynamics simulations were performed at 300 K for 100 ps with equilibration for 50 ps. The structural analyses of the trajectories indicate that the coordinate bond interactions, the hydrogen bond interactions, the van der Waals interactions as well as the hydrophobic interactions between ligand and receptor are responsible simultaneously for the preference of inhibition and potency. The ligand hydroxamate group is coordinated to the catalytic zinc ion and form stable hydrogen bonds with the carbonyl oxygen of Gly 162. The P1 group makes extensive van der Waals and hydrophobic contacts with the nonpolar side chains of several residues in the S1 subsite, including Leu 197, Val 198, Leu 218 and Tyr 223. Moreover, four to eight hydrogen bonds between hydroxamates and MMP-2 are formed to stabilize the inhibitors in the active site. Compared with the P2 and P3 groups, the P1 groups of inhibitors are oriented regularly, which is produced by the restrain of the S1 subsite. From the relationship between the length of the nonpolar P1 group and the biological activity, we confirm that MMP-2 has a pocket-like S1 subsite, not a channel-like S1 subsite proposed by Kiyama (Kiyama, R. et al., J. Med. Chem. 42 (1999), 1723). The energetic analyses show that the experimental binding free energies can be well correlated with the interactions between the inhibitors and their environments, which could be used as a simple score function to evaluate the binding affinities for other similar hydroxamates. The validity of the force field parameters and the MD simulations can be fully testified by the satisfactory agreements between the experimental structure-activity relationship and the information from the structural and energetic analyses. The information generated from the predicted complexes should be useful for further work in the area of structure-based design of new compounds.     

Integrated 3D-QSAR,molecular docking,and molecular dynamics simulation studies on 1,2,3-triazole based derivatives for designing new acetylcholinesterase inhibitors

Alzheimer’s disease (AD) is a multifactorial and polygenic disease. It is the most prevalent reason for dementia in the aging population. A dataset of twenty-six 1,2,3-triazole-based derivatives previously synthetized and evaluated for acetylcholinesterase inhibitory activity were subjected to the three-dimensional quantitative structure-activity relationship (3D-QSAR) study. Good predictability was achieved for comparative molecular field analysis (CoMFA) (Q² = 0.604, R² = 0.863, r_ext² = 0.701) and comparative molecular similarity indices analysis (CoMSIA) (Q² = 0.606, R² = 0.854, r_ext² = 0.647). The molecular features characteristics provided by the 3D-QSAR contour plots were quite useful for designing and improving the activity of acetylcholinesterase of this class. Based on these findings, a new series of 1,2,3-triazole based derivatives were designed, among which compound A1 with the highest predictive activity was subjected to detailed molecular docking and compared to the most active compound. The selected compounds were further subjected to 20 ns molecular dynamics (MD) simulations to study the comparative conformation dynamics of the protein after ligand binding, revealing promising results for the designed molecule. Therefore, this study could provide worthy guidance for further experimental analysis of highly effective acetylcholinesterase inhibitors.     

A structure-activity relationship study of catechol-O-methyltransferase inhibitors combining molecular docking and 3D QSAR methods

A panel of 92 catechol-O-methyltransferase (COMT) inhibitors was used to examine the molecular interactions affecting their biological activity. COMT inhibitors are used as therapeutic agents in the treatment of Parkinson's disease, but there are limitations in the currently marketed compounds due to adverse side effects. This study combined molecular docking methods with three-dimensional structure-activity relationships (3D QSAR) to analyse possible interactions between COMT and its inhibitors, and to incite the design of new inhibitors. Comparative molecular field analysis (CoMFA) and GRID/GOLPE models were made by using bioactive conformations from docking experiments, which yielded q2 values of 0.594 and 0.636, respectively. The docking results, the COMT X-ray structure, and the 3D QSAR models are in agreement with each other. The models suggest that an interaction between the inhibitor's catechol oxygens and the Mg2+ ion in the COMT active site is important. Both hydrogen bonding with Lys144, Asn170 and Glu199, and hydrophobic contacts with Trp38, Pro174 and Leu198 influence inhibitor binding. Docking suggests that a large R1 substituent of the catechol ring can form hydrophobic contacts with side chains of Val173, Leu198, Met201 and Val203 on the COMT surface. Our models propose that increasing steric volume of e.g. the diethylamine tail of entacapone is favourable for COMT inhibitory activity.     

Synthesis and molecular docking study of novel COVID-19 inhibitors

In 2020, the world tried to combat the corona virus (COVID-19) pandemic. A proven treatment method specific to Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) is still not found. In this study, seven new antiviral compounds were designed for COVID-19 treatment. The ability of these compounds to inhibit COVID-19’s RNA processing was calculated by the molecular docking study. It has been observed that the compounds can have high binding affinities especially against NSP12 (between -9.06 and -8.00 kcal/mol). The molecular dynamics simulation of NSP12-ZG 7 complex proved the stability of interaction. The synthesis of two most active molecules was performed by one-pot reaction and characterized by FT-IR, ¹H-NMR, ¹³C-NMR, and mass spectroscopy. The compounds presented with their synthesis are inhibitory core structures against SARS-CoV-2 infection.     

On the Topological Sub-Structural Molecular Design (TOSS-MODE) in QSPR/QSAR and Drug Design Research

Abstract

A recently introduced graph-theoretical approach to the study of structure-property-activity relationships is presented. The theoretical approach and the computational strategy for the use of the TOSS-MODE approach are given with details. Several QSPR and QSAR applications are reviewed including the study of physical properties of organic compounds, diamagnetic susceptibilities, and biological properties. The applications of the TOSS-MODE approach to discrimination of active/inactive compounds, the virtual screening of compounds with a desired property from databases of chemical structures, identification of active/inactive fragments and its relationships with 2D/3D pharmacophores, and to the design of novel compounds with desired biological activities are also reviewed.     

Monte Carlo method based QSAR modelling of natural lipase inhibitors using hybrid optimal descriptors

Obesity is one of the most provoking health burdens in the developed countries. One of the strategies to prevent obesity is the inhibition of pancreatic lipase enzyme. The aim of this study was to build QSAR models for natural lipase inhibitors by using the Monte Carlo method. The molecular structures were represented by the simplified molecular input line entry system (SMILES) notation and molecular graphs. Three sets – training, calibration and test set of three splits – were examined and validated. Statistical quality of all the described models was very good. The best QSAR model showed the following statistical parameters: r² = 0.864 and Q² = 0.836 for the test set and r² = 0.824 and Q² = 0.819 for the validation set. Structural attributes for increasing and decreasing the activity (expressed as pIC₅₀) were also defined. Using defined structural attributes, the design of new potential lipase inhibitors is also presented. Additionally, a molecular docking study was performed for the determination of binding modes of designed molecules.     

Molecular description of pyrimidine-based inhibitors with activity against FAK combining 3D-QSAR analysis,molecular docking and molecular dynamics

Focal adhesion kinase (FAK) is a promising target for developing more effective anticancer drugs. To better understand the structure-activity relationships and mechanism of actions of FAK inhibitors, a molecular modeling study using 3D-QSAR, molecular docking, molecular dynamics simulations, and binding free energy analysis were conducted. Two types of satisfactory 3D-QSAR models were generated, comprising the CoMFA model (R²_cv = 0.528, R²_pred = 0.7557) and CoMSIA model (R²_cv = 0.757, R²_pred = 0.8362), for predicting the inhibitory activities of novel inhibitors. The derived contour maps indicate structural characteristics for substituents on the template. Molecular docking, molecular dynamic simulations and binding free energy calculations further reveal that the binding of inhibitors to FAK is mainly contributed from hydrophobic, electrostatic and hydrogen bonding interactions. In addition, some key residues (Arg14, Glu88, Cys90, Arg138, Asn139, Leu141, and Leu155) responsible for ligand-receptor binding are highlighted. All structural information obtained from 3D-QSAR models and molecular dynamics is consist with the available experimental activities. All the results will facilitate the optimization of this series of FAK inhibitors with higher inhibitory activities.     

Effects of mussel-inspired co-deposition of 2-hydroxymethyl methacrylate and poly (2-methoxyethyl acrylate) on the hydrophilicity and binding tendency of common hemodialysis membranes: Molecular dynamics simulations and molecular docking studies

Despite advances in the field, hemoincompatibility remains a critical issue for hemodialysis (HD) as interactions between various human blood constituents and the polymeric structure of HD membranes results in complications such as activation of immune system cascades. Adding hydrophilic polymer structures to the membranes is one modification approach that can decrease the extent of protein adsorption. This study conducted molecular dynamics (MD) simulations to understand the interactions between three human serum proteins (fibrinogen [FB], human serum albumin, and transferrin) and common HD membranes in untreated and modified forms. Poly(aryl ether sulfone) (PAES) and cellulose triacetate were used as the common dialyzer polymers, and membrane modifications were performed with 2-hydroxymethyl methacrylate (HEMA) and poly (2-methoxyethyl acrylate) (PMEA), using polydopamine-assisted co-deposition. The MD simulations were used as the framework for binding energy simulations, and molecular docking simulations were also performed to conduct molecular-level investigations between the two modifying polymers (HEMA and PMEA) and FB. Each of the three proteins acted differently with the membranes due to their unique nature and surface chemistry. The simulations show PMEA binds less intensively to FB with a higher number of hydrogen bonds, which reflects PMEA's superior performance compared to HEMA. The simulations suggest PAES membranes could be used in modified forms for blood-contact applications as they reflect the lowest binding energy to blood proteins.     

Protein farnesyltransferase: Flexible docking studies on inhibitors using computational modeling

Summary Using MacroModel, peptide, peptidomimetic and non-peptidomimetic inhibitors of the zinc metalloenzyme, farnesyltransferase (FTase), were docked into the enzyme binding site. Inhibitor flexibility, farnesyl pyrophosphate substrate flexibility, and partial protein flexibility were taken into account in these docking studies. In addition to CVFM and CVIM, as well as our own inhibitors FTI-276 and FTI-2148, we have docked other farnesyltransferase inhibitors (FTIs) including Zarnestra, which presently is in advanced clinical trials. The AMBER* force field was employed, augmented with parameters that were derived for zinc. A single binding site model that was derived from the crystal structure of CVFM complexed with farnesyltransferase and farnesylpyrophosphate was used for these studies. The docking results using the lowest energy structure from the simulation, or one of the lowest energy structures, were generally in excellent agreement with the X-ray structures. One of the most important findings of this study is that numerous alternative conformations for the methionine side chain can be accommodated by the enzyme suggesting that the methionine pocket can tolerate groups larger than methionine at the C-terminus of the tetrapeptide and suggesting alternative locations for the placement of side chains that may improve potency.     

The MM2QM tool for combining docking,molecular dynamics,molecular mechanics,and quantum mechanics†

The use of the MM2QM tool in a combined docking + molecular dynamics (MD) + molecular mechanics (MM) + quantum mechanical (QM) binding affinity prediction study is presented, and the tool itself is discussed. The system of interest is Mycobacterium tuberculosis (MTB) pantothenate synthetase in complexes with three highly similar sulfonamide inhibitors, for which crystal structures are available. Starting from the structure of MTB pantothenate synthetase in the “open” conformation and following the combined docking + MD + MM + QM procedure, we were able to capture the closing of the enzyme binding pocket and to reproduce the position of the ligands with an average root mean square deviation of 1.6 Å. Protein–ligand interaction energies were reproduced with an average error lower than 10%. The discussion on the MD part and a protein flexibility importance is carried out. The presented approach may be useful especially for finding analog inhibitors or improving drug candidates. © 2012 Wiley Periodicals, Inc.     

Pharmacophore modelling,virtual screening and molecular docking studies on PLD1 inhibitors

Lipid metabolism plays a significant role in influenza virus replication and subsequent infection. The regulatory mechanism governing lipid metabolism and viral replication is not properly understood to date, but both Phospholipase D (PLD1 and PLD2) activities are stimulated in viral infection. In vitro studies indicate that chemical inhibition of PLD1 delays viral entry and reduction of viral loads. The current study reports a three-dimensional pharmacophore model based on 35 known PLD1 inhibitors. A sub-set of 25 compounds was selected as the training set and the remaining 10 compounds were kept in the test set. One hundred and twelve pharmacophore models were generated; a six-featured pharmacophore model (AADDHR.57) with survival score (2.69) produced a statistically significant three-dimensional quantitative structure–activity relationship model with r² = 0.97 (internal training set), r² = 0.71 (internal test set) and Q² = 0.64. The predictive power of the pharmacophore model was validated with an external test set (r² = 0.73) and a systematic virtual screening work-flow was employed showing an enrichment factor of 23.68 at the top 2% of the dataset (active and decoys). Finally, the model was used for screening of the filtered PubChem database to fetch molecules which can be proposed as potential PLD1 inhibitors for blocking influenza infection.     

Exploring the effect of aplidin on low molecular weight protein tyrosine phosphatase by molecular docking and molecular dynamic simulation study

The low molecular weight protein tyrosine phosphatase (LMW-PTP) could regulate many signaling pathways, and it had drawn attention as a potential target for cancer. As previous report has indicated that the aplidin could inhibit the LMW-PTP, and thus, the relevant cancer caused by the abnormal regulation of the LMW-PTP could be remission. However, the molecular mechanism of inhibition of the LMW-PTP by the aplidin had not been fully understood. In this study, various computational approaches, namely molecular docking, MDs and post-dynamic analyses were utilized to explore the effect of the aplidin on the LMW-PTP. The results suggested that the intramolecular interactions of the residues in the two sides of the active site (Ser43-Ala55 and Pro121-Asn134) and the P-loop region (Leu13-Ser19) in the LMW-PTP was disturbed owing to the aplidin, meanwhile, the π-π interaction between Tyr131 and Tyr132 might be broken. The Asn15 might be the key residue to break the residues interactions. In a word, this study may provide more information for understanding the effect of inhibition of the aplidin on the LMW-PTP.     

表皮生长因子受体和抑制剂之间分子对接的研究

研究了表皮生长因子受体(EGFR)和4-苯胺喹唑啉类抑制剂之间的相互作用模式，表皮生长因子受体的三维结构通过同源蛋白模建的方法得到，而抑制剂和靶酶结合复合物结构则通过分子力学和分子动力学结合的方法计算得到。从模拟结果得到的抑制剂和靶酶之间的相互作用模式表明范德华相互作用、疏水相互作用以及氢键相互作用对抑制剂的活性都有重要的影响，抑制剂的苯胺部分位于活性口袋的底部，能够与受体残基的非极性侧链产生很强的范德华和疏水相互作用，抑制剂双环上的取代基团也能和活性口袋外部的部分残基形成一定的范德华和疏水性相互作用，而抑制剂喹唑啉环上的氮原子能和周围的残基形成较强的氢键相互作用，对抑制剂的活性有较大的影响，计算得到抑制剂和靶酶之间的非键相互作用能以及抑制剂和靶酶之间的相互作用信息能够很好地解释抑制剂活性和结构的关系，为全新抑制剂的设计提供了重要的结构信息。