1
|
Wang L, Lu D, Wang Y, Xu X, Zhong P, Yang Z. Binding selectivity-dependent molecular mechanism of inhibitors towards CDK2 and CDK6 investigated by multiple short molecular dynamics and free energy landscapes. J Enzyme Inhib Med Chem 2023; 38:84-99. [DOI: 10.1080/14756366.2022.2135511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Affiliation(s)
- Lifei Wang
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Dan Lu
- Department of Physics, Jiangxi Agricultural University, Nanchang, PR China
| | - Yan Wang
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Xiaoyan Xu
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Peihua Zhong
- College of Computer Information and Engineering, Jiangxi Agriculture University, Nanchang, PR China
| | - Zhiyong Yang
- Department of Physics, Jiangxi Agricultural University, Nanchang, PR China
| |
Collapse
|
2
|
Ghorbanpour M, Soltani B, Mota A, Jahanbin Sardroodi J, Mehdizadeh Aghdam E, Shayanfar A, Molavi O, Mohammad-Rezaei R, Ebadi-Nahari M, Ziegler CJ. Copper (II) complexes with N, S donor pyrazole-based ligands as anticancer agents. Biometals 2022; 35:1095-1111. [PMID: 36001216 DOI: 10.1007/s10534-022-00426-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2021] [Accepted: 07/19/2022] [Indexed: 11/25/2022]
Abstract
A group of bidentate nitrogen and sulfur donor pyrazole derivative ligands abbreviated as Na[RNCS(Pz)], Na[RNCS(PzMe2)], Na[RNCS(PzMe3)], Na[RNCS(PzPhMe)], Na[RNCS(PzPh2)], where (R = Et, Ph), and their Cu (II) complexes were synthesized and characterized by spectroscopic and physicochemical methods. The crystal structure of [Cu(PhNCSPzMe3)2] was determined by X-ray crystallography analysis and the results described a distorted square planar coordination geometry for this complex. Also, the cyclic voltammetry investigations indicated that the synthesized copper complex is an electrochemically active species. Moreover, the cytotoxic activity of all of the twenty synthesized compounds was evaluated using MTT assay against the MCF-7 (human breast carcinoma) cell lines, in vitro. Cu (II) complexes indicate significant cytotoxicity against the MCF-7 cell lines as compared with the free ligands. The docking studies showed that the copper complexes have better interactions with EGFR and CDK2 proteins, compared to the free ligands, and most of the studied compounds have a higher value of binding energy relative to the studied controls. The results of QSAR analysis suggest that dipole moment is in direct correlation with the obtained IC50 values, and it strongly impact the anticancer effects generated by the compounds. Our findings suggest that the developed copper complexes can be good candidates for further evaluations as chemotherapeutic agents in the treatment of cancer.
Collapse
Affiliation(s)
- Monireh Ghorbanpour
- Department of Chemistry, Faculty of Basic Science, Azarbaijan Shahid Madani University, P. O. Box 53714-161, Tabriz, Iran
| | - Behzad Soltani
- Department of Chemistry, Faculty of Basic Science, Azarbaijan Shahid Madani University, P. O. Box 53714-161, Tabriz, Iran.
| | - Ali Mota
- Department of Clinical Biochemistry and Laboratory Medicine, Faculty of Medicine, Tabriz University of Medical Sciences, Tabriz, Iran
| | - Jaber Jahanbin Sardroodi
- Department of Chemistry, Faculty of Basic Science, Azarbaijan Shahid Madani University, P. O. Box 53714-161, Tabriz, Iran
| | - Elnaz Mehdizadeh Aghdam
- Department of Pharmaceutical Biotechnology, Faculty of Pharmacy, Tabriz University of Medical Sciences, Tabriz, Iran.,Molecular Medicine Research Center, Biomedicine Institute, Tabriz University of Medical Sciences, Tabriz, Iran
| | - Ali Shayanfar
- Pharmaceutical Analysis Research Center, Tabriz University of Medical Sciences, Tabriz, Iran
| | - Ommoleila Molavi
- Department of Pharmaceutical Biotechnology, Faculty of Pharmacy, Tabriz University of Medical Sciences, Tabriz, Iran
| | - Rahim Mohammad-Rezaei
- Department of Chemistry, Faculty of Basic Science, Azarbaijan Shahid Madani University, P. O. Box 53714-161, Tabriz, Iran
| | - Mostafa Ebadi-Nahari
- Department of Biology, Faculty of Basic Sciences, Azarbaijan Shahid Madani University, Tabriz, Iran
| | | |
Collapse
|
3
|
Wang T, Wu F, Luo L, Zhang Y, Ma J, Hu Y. Efficient synthesis and cytotoxic activity of polysubstituted thieno[2,3-d]pyrimidine derivatives. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.132497] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
4
|
Tang X, Zheng A, Wu F, Liao C, Hu Y, Luo C. Synthesis and anticancer activities of diverse furo[2,3- d]pyrimidine and benzofuro[3,2- d]pyrimidine derivatives. SYNTHETIC COMMUN 2022. [DOI: 10.1080/00397911.2022.2060117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Xiaoyu Tang
- School of Pharmaceutical Sciences and Institute of Medicinal Chemistry, Hubei University of Medicine, Shiyan, China
- Hubei Key Laboratory of Wudang Local Chinese Medicine Research, Hubei University of Medicine, Shiyan, China
| | - Aihua Zheng
- School of Pharmaceutical Sciences and Institute of Medicinal Chemistry, Hubei University of Medicine, Shiyan, China
- Hubei Key Laboratory of Wudang Local Chinese Medicine Research, Hubei University of Medicine, Shiyan, China
| | - Fengxu Wu
- School of Pharmaceutical Sciences and Institute of Medicinal Chemistry, Hubei University of Medicine, Shiyan, China
- Hubei Key Laboratory of Wudang Local Chinese Medicine Research, Hubei University of Medicine, Shiyan, China
| | - Chujie Liao
- School of Pharmaceutical Sciences and Institute of Medicinal Chemistry, Hubei University of Medicine, Shiyan, China
- Hubei Key Laboratory of Wudang Local Chinese Medicine Research, Hubei University of Medicine, Shiyan, China
| | - Yanggen Hu
- School of Pharmaceutical Sciences and Institute of Medicinal Chemistry, Hubei University of Medicine, Shiyan, China
- Hubei Key Laboratory of Wudang Local Chinese Medicine Research, Hubei University of Medicine, Shiyan, China
| | - Chao Luo
- School of Pharmaceutical Sciences and Institute of Medicinal Chemistry, Hubei University of Medicine, Shiyan, China
- School of Basic Medical Sciences, Hubei University of Medicine, Shiyan, China
| |
Collapse
|
5
|
Kumari G, Nigam VK, Pandey DM. The molecular docking and molecular dynamics study of flavonol synthase and flavonoid 3'-monooxygenase enzymes involved for the enrichment of kaempferol. J Biomol Struct Dyn 2022; 41:2478-2491. [PMID: 35105279 DOI: 10.1080/07391102.2022.2033324] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
Kaempferol is a natural flavonol that shows many pharmacological properties including anti-inflammatory, antioxidant, anticancer, antidiabetic activities etc. It has been reported in many vegetables, fruits, herbs and medicinal plants. The enzyme flavonol synthase (FLS, EC 1.14.20.6) catalyses the conversion of dihydroflavonols to flavonols. Whereas flavonoid 3'-monooxygenase (F3'H, EC 1.14.14.82) catalyses the hydroxylation of dihydroflavonol, and flavonol. FLS is involved in the synthesis of the kaempferol whereas F3'H causes degradation of kaempferol. The present study aimed to analyse the binding affinity, stability and activating activity of enzyme FLS as well as inhibitory activity of enzyme F3'H involved in the enrichment of the kaempferol using the in-silico approaches. Computational study for physico-chemical properties, conserved domain identification, 3-D structure prediction and its validation, conservation analysis, molecular docking followed by molecular dynamics analysis of FLS and F3'H, protein-activator (FLS-LIG Complex) and protein-inhibitor (F3'H-LIG Complex) complexes have been performed. Other structural analyses like root mean square fluctuation (RMSF), root mean square deviation (RMSD), surface area solvent accessibility (SASA), radius of gyration (Rg), hydrogen bond analysis, principal component analysis (PCA), Poisson-Boltzmann analysis (MM_PBSA) and the dynamic cross correlation map (DCCM) analysis to explore the structural, functional and thermodynamic stability of the proteins and the complexes were also studied. The molecular docking result showed that FLS binds strongly with the activator ascorbate (CID _54670067) while F3'H binds with the inhibitor ketoconazole (CID_456201). The most powerful inhibitor (ketoconazole for F3'H) and activator (ascorbate for FLS) is determined by computing the thermodynamic binding free energy through MM_PBSA analysis. The current work provides wide-ranging structural and functional information about FLS and F3'H enzymes showing detailed molecular mechanism of kaempferol biosynthesis and its degradation and hence kaempferol enrichment. Finding of the present work opens up new possibilities for future research towards enrichment of kaempferol by using activator (ascorbate) for FLS and inhibitor (ketoconazole) for F3'H as well as for its large-scale production using in vitro approaches.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Garima Kumari
- Department of Bioengineering and Biotechnology, Birla Institute of Technology, Mesra, Jharkhand, India
| | - Vinod Kumar Nigam
- Department of Bioengineering and Biotechnology, Birla Institute of Technology, Mesra, Jharkhand, India
| | - Dev Mani Pandey
- Department of Bioengineering and Biotechnology, Birla Institute of Technology, Mesra, Jharkhand, India
| |
Collapse
|
6
|
Liang S, Liu X, Zhang S, Li M, Zhang Q, Chen J. Binding mechanism of inhibitors to SARS-CoV-2 main protease deciphered by multiple replica molecular dynamics simulations. Phys Chem Chem Phys 2022; 24:1743-1759. [PMID: 34985081 DOI: 10.1039/d1cp04361g] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The outbreak caused by SARS-CoV-2 has received extensive worldwide attention. As the main protease (Mpro) in SARS-CoV-2 has no human homologues, it is feasible to reduce the possibility of targeting the host protein by accidental drugs. Thus, Mpro has been an attractive target of efficient drug design for anti-SARS-CoV-2 treatment. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), and the molecular mechanics-generalized Born surface area (MM-GBSA) method were integrated together to decipher the binding mechanism of four inhibitors masitinib, O6K, FJC and GQU to Mpro. The results indicate that the binding of four inhibitors clearly affects the structural flexibility and internal dynamics of Mpro along with dihedral angle changes of key residues. The analysis of FELs unveils that the stability in the relative orientation and geometric position of inhibitors to Mpro is favorable for inhibitor binding. Residue-based free energy decomposition reveals that the inhibitor-Mpro interaction networks involving hydrogen bonding interactions and hydrophobic interactions provide significant information for the design of potent inhibitors against Mpro. The hot spot residues including H41, M49, F140, N142, G143, C145, H163, H164, M165, E166 and Q189 identified by computational alanine scanning are considered as reliable targets of clinically available inhibitors inhibiting the activities of Mpro.
Collapse
Affiliation(s)
- Shanshan Liang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Meng Li
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250357, China.
| |
Collapse
|
7
|
Yu YX, Liu WT, Li HY, Wang W, Sun HB, Zhang LL, Wu SL. Decoding molecular mechanism underlying binding of drugs to HIV-1 protease with molecular dynamics simulations and MM-GBSA calculations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2021; 32:889-915. [PMID: 34551634 DOI: 10.1080/1062936x.2021.1979647] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Accepted: 09/08/2021] [Indexed: 06/13/2023]
Abstract
HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25', Gly27/Gly27', Ala28/Ala28', Asp29, Ile47/Ile47', Gly49/Gly49', Ile50/Ile50', Val82/Val82' and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.
Collapse
Affiliation(s)
- Y X Yu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W T Liu
- Shuifa Qilu Cultural Tourism Development Co., Ltd, Shuifa Ecological Industry Group, Jinan, China
| | - H Y Li
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L L Zhang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - S L Wu
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
8
|
Zhao J, Wang L, Bao H, Chen J. Cluster analysis on conformational changes of the GDP/KRAS complex induced by A59G and D33E. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138995] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
9
|
Yu YX, Wang W, Sun HB, Zhang LL, Wu SL, Liu WT. Insights into effect of the Asp25/Asp25' protonation states on binding of inhibitors Amprenavir and MKP97 to HIV-1 protease using molecular dynamics simulations and MM-GBSA calculations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2021; 32:615-641. [PMID: 34157882 DOI: 10.1080/1062936x.2021.1939149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/25/2021] [Accepted: 06/02/2021] [Indexed: 06/13/2023]
Abstract
The protonation states of two aspartic acids in the catalytic strands of HIV-1 protease (PR) remarkably affect bindings of inhibitors to PR. It is requisite for the design of potent inhibitors towards PR to investigate the influences of Asp25/Asp25' protonated states on dynamics behaviour of PR and binding mechanism of inhibitors to PR. In this work, molecular dynamics (MD) simulations, MM-GBSA method and principal component (PC) analysis were coupled to explore the effect of Asp25/Asp25' protonation states on conformational changes of PR and bindings of Amprenavir and MKP97 to PR. The results show that the Asp25/Asp25' protonation states exert different impacts on structural fluctuations, flexibility and motion modes of PR. Dynamics analysis verifies that Asp25/Asp25' protonated states highly affect conformational dynamics of two flaps in PR. The binding free energy calculations results suggest that the Asp25/Asp25' protonated states obviously strengthen bindings of inhibitors to PR compared to the non-protonation state. Calculations of residue-based free energy decomposition indicate that the Asp25/Asp25' protonation not only disturbs the interaction network of inhibitors with PR but also stabilizes bindings of inhibitors to PR by cancelling the electrostatic repulsive interaction. Therefore, special attentions should be paid to the Asp25/Asp25' protonation in the design of potent inhibitors towards PR.
Collapse
Affiliation(s)
- Y X Yu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L L Zhang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - S L Wu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W T Liu
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
10
|
Chen J, Zhang S, Wang W, Sun H, Zhang Q, Liu X. Binding of Inhibitors to BACE1 Affected by pH-Dependent Protonation: An Exploration from Multiple Replica Gaussian Accelerated Molecular Dynamics and MM-GBSA Calculations. ACS Chem Neurosci 2021; 12:2591-2607. [PMID: 34185514 DOI: 10.1021/acschemneuro.0c00813] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
To date, inhibiting the activity of β-amyloid cleaving enzyme 1 (BACE1) has been considered an efficient approach for treating Alzheimer's disease (AD). In the current work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations and the molecular mechanics general Born surface area (MM-GBSA) method were combined to investigate the effect of pH-dependent protonation on the binding of the inhibitors CS9, C6U, and 6WE to BACE1. Dynamic analyses based on the MR-GaMD trajectory show that pH-dependent protonation strongly affects the structural flexibility, correlated motions, and dynamic behavior of inhibitor-bound BACE1. According to the constructed free energy profiles, in the protonated state at low pH, inhibitor-bound BACE1 tends to populate at more conformations than in high pH. The binding free energies calculated by MM-GBSA suggest that inhibitors possess stronger binding abilities under the protonation conditions at high pH than under the protonation conditions at low pH. Moreover, pH-dependent protonation exerts a significant effect on the hydrogen bonding interactions of CS9, C6U, and 6WE to BACE1, which correspondingly alters the binding abilities of the three inhibitors to BACE1. Furthermore, in different protonated environments, three inhibitors share common interaction clusters and similar binding sites in BACE1, which are reliably used as efficient targets for the design of potent inhibitors of BACE1.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China
| |
Collapse
|
11
|
Chen J, Wang L, Wang W, Sun H, Pang L, Bao H. Conformational transformation of switch domains in GDP/K-Ras induced by G13 mutants: An investigation through Gaussian accelerated molecular dynamics simulations and principal component analysis. Comput Biol Med 2021; 135:104639. [PMID: 34247129 DOI: 10.1016/j.compbiomed.2021.104639] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Revised: 07/05/2021] [Accepted: 07/05/2021] [Indexed: 10/20/2022]
Abstract
Mutations in K-Ras are involved in a large number of all human cancers, thus, K-Ras is regarded as a promising target for anticancer drug design. Understanding the target roles of K-Ras is important for providing insights on the molecular mechanism underlying the conformational transformation of the switch domains in K-Ras due to mutations. In this study, multiple replica Gaussian accelerated molecular (MR-GaMD) simulations and principal component analysis (PCA) were applied to probe the effect of G13A, G13D and G13I mutations on conformational transformations of the switch domains in GDP-associated K-Ras. The results suggest that G13A, G13D and G13I enhance the structural flexibility of the switch domains, change the correlated motion modes of the switch domains and strengthen the total motion strength of K-Ras compared with the wild-type (WT) K-Ras. Free energy landscape analyses not only show that the switch domains of the GDP-bound inactive K-Ras mainly exist as a closed state but also indicate that mutations evidently alter the free energy profile of K-Ras and affect the conformational transformation of the switch domains between the closed and open states. Analyses of hydrophobic interaction contacts and hydrogen bonding interactions show that the mutations scarcely change the interaction network of GDP with K-Ras and only disturb the interaction of GDP with the switch (SW1). In summary, two newly introduced mutations, G13A and G13I, play similar adjustment roles in the conformational transformations of two switch domains to G13D and are possibly utilized to tune the activity of K-Ras and the binding of guanine nucleotide exchange factors.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250357, China.
| | - Lifei Wang
- School of Science, Shandong Jiaotong University, Jinan, 250357, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan, 250357, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan, 250357, China
| | - Laixue Pang
- School of Science, Shandong Jiaotong University, Jinan, 250357, China
| | - Huayin Bao
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, 250355, China.
| |
Collapse
|
12
|
Wang L, Wang Y, Yang Z, Xu S, Li H. Binding Selectivity of Inhibitors toward Bromodomains BAZ2A and BAZ2B Uncovered by Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations. ACS OMEGA 2021; 6:12036-12049. [PMID: 34056358 PMCID: PMC8154142 DOI: 10.1021/acsomega.1c00687] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Accepted: 04/19/2021] [Indexed: 06/12/2023]
Abstract
Two Bromodomain-Containing proteins BAZ2A and BAZ2B are responsible for remodeling chromatin and regulating noncoding RNAs. As for our current studies, integration of multiple short molecular dynamics simulations (MSMDSs) with molecular mechanics generalized Born surface area (MM-GBSA) method is adopted for insights into binding selectivity of three small molecules D8Q, D9T and UO1 to BAZ2A against BAZ2B. The calculations of MM-GBSA unveil that selectivity of inhibitors toward BAZ2A and BAZ2B highly depends on the enthalpy changes and the details uncover that D8Q has better selectivity toward BAZ2A than BAZ2B, D9T more favorably bind to BAZ2B than BAZ2A, and UO1 does not show obvious selectivity toward these two proteins. The analysis of interaction network between residues and inhibitors indicates that seven residues are mainly responsible for the selectivity of D8Q, six residues for D9T and four residues provide significant contributions to associations of UO1 with two proteins. Moreover the analysis of interaction network not only reveals warm spots of inhibitor bindings to BAZ2A and BAZ2B but also unveils that common residue pairs, including (W1816, W1887), (P1817, P1888), (F1818, F1889), (V1822, V1893), (N1823, N1894),(L1826, L1897), (V1827, V1898), (F1872, F1943), (N1873, N1944) and (V1879, I1950) belonging to (BAZ2A, BAZ2B), induce mainly binding differences of inhibitors to BAZ2A and BAZ2B. Hence, insights from our current studies offer useful dynamics information relating with conformational alterations and structure-affinity relationship at atomistic levels for novel therapeutic strategies toward BAZ2A and BAZ2B.
Collapse
Affiliation(s)
- Lifei Wang
- School
of Science, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| | - Yan Wang
- School
of Science, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| | - Zhiyong Yang
- Department
of Physics, Jiangxi Agricultural University, 1101 Zhimin Road, Economic and Technological
Development Zone, Nanchang, Jiangxi Province 330045, China
| | - Shuobo Xu
- School
of Information Science and Electrical Engineering, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| | - Hongyun Li
- School
of Science, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| |
Collapse
|
13
|
Wu S, Wang L, Zhang L, Xu X, Zhao J. Molecular dynamics insights into binding selectivity of inhibitors toward BRD4 and CBP. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
14
|
Chen J, Zhang S, Wang W, Pang L, Zhang Q, Liu X. Mutation-Induced Impacts on the Switch Transformations of the GDP- and GTP-Bound K-Ras: Insights from Multiple Replica Gaussian Accelerated Molecular Dynamics and Free Energy Analysis. J Chem Inf Model 2021; 61:1954-1969. [PMID: 33739090 DOI: 10.1021/acs.jcim.0c01470] [Citation(s) in RCA: 81] [Impact Index Per Article: 27.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Mutations yield significant effect on the structural flexibility of two switch domains, SW1 and SW2, in K-Ras, which is considered as an important target of anticancer drug design. To unveil a molecular mechanism with regard to mutation-mediated tuning on the activity of K-Ras, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations followed by analysis of free energy landscapes (FELs) are performed on the GDP- and GTP-bound wild-type (WT), G12V, and D33E K-Ras. The results suggest that G12V and D33E not only evidently change the flexibility of SW1 and SW2 but also greatly affect correlated motions of SW1 and SW2 separately relative to the P-loop and SW1, which exerts a certain tuning on the activity of K-Ras. The information stemming from the analyses of FELs reveals that the conformations of SW1 and SW2 are in high disorders in the GDP- and GTP-associated WT and mutated K-Ras, possibly producing significant effect on binding of guanine nucleotide exchange factors or effectors to K-Ras. The interaction networks of GDP and GTP with K-Ras are identified and the results uncover that the instability in hydrogen-bonding interactions of SW1 with GDP and GTP is mostly responsible for conformational disorder of SW1 and SW2 as well as tunes the activity of oncogenic K-Ras.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Laixue Pang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China
| |
Collapse
|
15
|
Wang Y, Wu S, Wang L, Yang Z, Zhao J, Zhang L. Binding selectivity of inhibitors toward the first over the second bromodomain of BRD4: theoretical insights from free energy calculations and multiple short molecular dynamics simulations. RSC Adv 2020; 11:745-759. [PMID: 35423696 PMCID: PMC8693360 DOI: 10.1039/d0ra09469b] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2020] [Accepted: 12/17/2020] [Indexed: 12/12/2022] Open
Abstract
Bromodomain-containing protein 4 (BRD4) plays an important role in mediating gene transcription involved in cancers and non-cancer diseases such as acute heart failure and inflammatory diseases. In this work, multiple short molecular dynamics (MSMD) simulations are integrated with a molecular mechanics generalized Born surface area (MM-GBSA) approach to decipher binding selectivity of three inhibitors 8NS, 82Y, and 837 toward two domains BD1 and BD2 of BRD4. The results demonstrate that the enthalpy effects play critical roles in selectivity identification of inhibitors toward BD1 and BD2, determining that 8NS has better selectivity toward BD2 than BD1, while 82Y and 837 more favorably bind to BD1 than BD2. A residue-based free-energy decomposition method was used to calculate an inhibitor-residue interaction spectrum and unveil contributions of separate residues to binding selectivity. The results identify six common residues, containing (P82, P375), (V87, V380), (L92, L385), (L94, L387), (N140, N433), and (I146, V439) individually belonging to (BD1, BD2) of BRD4, and yield a considerable binding difference of inhibitors to BD1 and BD2, suggesting that these residues play key roles in binding selectivity of inhibitors toward BD1 and BD2 of BRD4. Therefore, these results provide useful dynamics information and a structure affinity relationship for the development of highly selective inhibitors targeting BD1 and BD2 of BRD4.
Collapse
Affiliation(s)
- Yan Wang
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Shiliang Wu
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Lifei Wang
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Zhiyong Yang
- Department of Physics, Jiangxi Agricultural University Nanchang 330045 China
| | - Juan Zhao
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Lulu Zhang
- School of Science, Shandong Jiaotong University Jinan 250357 China
| |
Collapse
|
16
|
Zhao J, Sun H, Wang W, Zhang L, Chen J. Theoretical insights into mutation-mediated conformational changes of the GNP-bound H-RAS. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.138042] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
17
|
Chen J, Wang W, Sun H, Pang L, Yin B. Mutation-mediated influences on binding of anaplastic lymphoma kinase to crizotinib decoded by multiple replica Gaussian accelerated molecular dynamics. J Comput Aided Mol Des 2020; 34:1289-1305. [DOI: 10.1007/s10822-020-00355-5] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2020] [Accepted: 10/14/2020] [Indexed: 12/19/2022]
|
18
|
Chen J, Yin B, Wang W, Sun H. Effects of Disulfide Bonds on Binding of Inhibitors to β-Amyloid Cleaving Enzyme 1 Decoded by Multiple Replica Accelerated Molecular Dynamics Simulations. ACS Chem Neurosci 2020; 11:1811-1826. [PMID: 32459964 DOI: 10.1021/acschemneuro.0c00234] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The β-amyloid cleaving enzyme 1 (BACE1) has been thought to be an efficient target for treatment of Alzheimer's disease (AD). Deep insight into inhibitor-BACE1 binding mechanism is of significance for design of potent drugs toward BACE1. In this work, multiple replica accelerated molecular dynamics (MR-aMD) simulations, principal component (PC) analysis, and free energy landscapes were integrated to decode the effect of disulfide bonds (SSBs) in BACE1 on bindings of three inhibitors 3KO, 3KT, and 779 to BACE1. The results from cross-correlation analysis suggest that the breaking of SSBs exerts significant influence on structural flexibility and internal dynamics of inhibitor-bound BACE1. PC analysis and free energy landscapes reveal that the breaking of SSBs not only evidently induces the conformational rearrangement of BACE1 but also highly changes binding poses of three inhibitors in BACE1 and leads to more disordered binding of three inhibitors to BACE1, which is further supported by the increase in binding entropy of inhibitors to BACE1 due to the breaking of SSBs. Residue-based free energy decomposition method was utilized to evaluate contributions of separate residues to inhibitor-BACE1 binding. The results suggest that although the breaking of SSBs in BACE1 does not destroy the interaction network of inhibitors with BACE1 it changes interaction strength of some residues with inhibitors. Meanwhile, the information from residue-based free energy decomposition indicates that residues L91, S96, V130, Y132, Q134, W137, F169, I171, and I179 can be used as efficient targets of drug design toward BACE1.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Baohua Yin
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan 250355, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| |
Collapse
|
19
|
|
20
|
Ngo K, Collins-Kautz C, Gerstenecker S, Wagner B, Heine A, Klebe G. Protein-Induced Change in Ligand Protonation during Trypsin and Thrombin Binding: Hint on Differences in Selectivity Determinants of Both Proteins? J Med Chem 2020; 63:3274-3289. [DOI: 10.1021/acs.jmedchem.9b02061] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Khang Ngo
- Institute of Pharmaceutical Chemistry, Philipps-University Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| | - Chelsey Collins-Kautz
- Institute of Pharmaceutical Chemistry, Philipps-University Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| | - Stefan Gerstenecker
- Institute of Pharmaceutical Chemistry, Philipps-University Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| | - Björn Wagner
- Pharma Research Non-Clinical Safety, F. Hoffmann-La Roche AG, 4070 Basel, Switzerland
| | - Andreas Heine
- Institute of Pharmaceutical Chemistry, Philipps-University Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| | - Gerhard Klebe
- Institute of Pharmaceutical Chemistry, Philipps-University Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| |
Collapse
|
21
|
Rasafar N, Barzegar A, Mehdizadeh Aghdam E. Design and development of high affinity dual anticancer peptide-inhibitors against p53-MDM2/X interaction. Life Sci 2020; 245:117358. [PMID: 32001262 DOI: 10.1016/j.lfs.2020.117358] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2019] [Revised: 01/22/2020] [Accepted: 01/24/2020] [Indexed: 12/16/2022]
Abstract
AIMS Inhibition of P53-MDM2/X interaction is known as an effective cancer therapy strategy. In this regard, pDI peptide was introduced previously with the potential of targeting MDM2. In this research, the large-scale peptide mutation screening was used to achieve the best sequence of pDI with the highest affinity for inhibition activity against MDM2/X. MAIN METHODS Three mutant peptides of pDI as dual inhibitor peptides including single mutations of pDIm/4W, pDIm/11M and double mutations of pDIdm/4W11M were presented with the high affinities to inhibit both MDM2/X. The selected mutants were then evaluated comprehensively to confirm their ability as potent MDM2/X inhibitors, using a theoretical simulation approach. KEY FINDINGS MD simulations analyses confirmed their dual inhibition potential against both MDM2/X interactions with p53 protein. The developed pDIm and mainly pDIdm peptides showed stable conformations over the simulation time with conserved secondary structure and effective interaction with MDM2/X by physical binding such as hydrogen bonding. Besides, umbrella sampling free energy calculation indicated higher binding energy, ΔGbinding, of pDIm-MDM2/X and pDIdm-MDM2/X compared to pDI-MDM2/X. SIGNIFICANCE The optimized and improved mutant pDI, pDIdm, with more effective ΔGbinding values of -30 and -25 kcal/mol to MDMX and MDM2, respectively, is recommended as a promising anticancer agent and suitable candidate for experimental evaluations.
Collapse
Affiliation(s)
- Nasim Rasafar
- Research Institute of Bioscience and Biotechnology, University of Tabriz, Tabriz, Iran
| | - Abolfazl Barzegar
- Research Institute of Bioscience and Biotechnology, University of Tabriz, Tabriz, Iran.
| | - Elnaz Mehdizadeh Aghdam
- Molecular Medicine Research Center, Biomedicine Institute, Tabriz University of Medical Sciences, Tabriz, Iran; Department of Pharmaceutical Biotechnology, Faculty of Pharmacy, Tabriz University of Medical Sciences, Tabriz, Iran.
| |
Collapse
|
22
|
Chen J, Liu X, Zhang S, Chen J, Sun H, Zhang L, Zhang Q. Molecular mechanism with regard to the binding selectivity of inhibitors toward FABP5 and FABP7 explored by multiple short molecular dynamics simulations and free energy analyses. Phys Chem Chem Phys 2020; 22:2262-2275. [DOI: 10.1039/c9cp05704h] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Recently, fatty acid binding proteins 5 and 7 (FABP5 and FABP7) have been regarded as the prospective targets for clinically treating multiple diseases related to FABPs.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science
- Shandong Jiaotong University
- Jinan 250357
- People's Republic of China
| | - Xinguo Liu
- School of Physics and Electronics
- Shandong Normal University
- Jinan
- People's Republic of China
| | - Shaolong Zhang
- School of Physics and Electronics
- Shandong Normal University
- Jinan
- People's Republic of China
| | - Junxiao Chen
- School of Chemistry and Pharmaceutical Engineering
- Qilu University of Technology
- Jinan
- People's Republic of China
| | - Haibo Sun
- School of Science
- Shandong Jiaotong University
- Jinan 250357
- People's Republic of China
| | - Lin Zhang
- School of Construction Machinery
- Shandong Jiaotong University
- Jinan 250357
- People's Republic of China
| | - Qinggang Zhang
- School of Physics and Electronics
- Shandong Normal University
- Jinan
- People's Republic of China
| |
Collapse
|
23
|
Yan F, Liu X, Zhang S, Zhang Q, Chen J. Understanding conformational diversity of heat shock protein 90 (HSP90) and binding features of inhibitors to HSP90 via molecular dynamics simulations. Chem Biol Drug Des 2019; 95:87-103. [PMID: 31560152 DOI: 10.1111/cbdd.13623] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2019] [Revised: 08/06/2019] [Accepted: 08/26/2019] [Indexed: 12/23/2022]
Abstract
Heat shock protein 90 (HSP90) is a promising target for treatment of cancer, and inhibitor bindings can generate efficient suppression on tumor in multiple ways. In this work, 140-ns molecular dynamics simulations were performed on six systems. Principal component analysis was subsequently carried out to explore the conformational diversity of HSP90. The results suggest that inhibitor bindings induce large conformational changes of HSP90, which tends to enlarge the volume of the binding pocket to facilitate the entrance of inhibitors. Hierarchical clustering analyses, the calculation of the energy contribution of each atom, and the analyses of hydrogen-bonding interactions were performed. The results indicate that 20 residues in group A of the hierarchical tree are responsible for major contributions, and van der Waals interactions as well as hydrogen-bonding interactions between important residues in HSP90 and key regions of inhibitors are the main force for promoting inhibitor bindings. We expect that this work can provide useful theoretical information for development of efficient inhibitors targeting HSP90.
Collapse
Affiliation(s)
- Fangfang Yan
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
24
|
Chen J, Wang J, Yin B, Pang L, Wang W, Zhu W. Molecular Mechanism of Binding Selectivity of Inhibitors toward BACE1 and BACE2 Revealed by Multiple Short Molecular Dynamics Simulations and Free-Energy Predictions. ACS Chem Neurosci 2019; 10:4303-4318. [PMID: 31545898 DOI: 10.1021/acschemneuro.9b00348] [Citation(s) in RCA: 58] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
The β-amyloid cleaving enzymes 1 and 2 (BACE1 and BACE2) have been regarded as the prospective targets for clinically treating Alzheimer's disease (AD) in the last two decades. Thus, insight into the binding differences of inhibitors to BACE1 and BACE2 is of significance for designing highly selective inhibitors toward the two proteins. In this work, multiple short molecular dynamics (MSMD) simulations are coupled with the molecular mechanics generalized Born surface area (MM-GBSA) method to probe the binding selectivity of three inhibitors DBO, CS9, and SC7 on BACE1 over BACE2. The results show that the entropy effect plays a key role in selectivity identification of inhibitors toward BACE1 and BACE2, which determines that DBO has better selectivity toward BACE2 over BACE1, while CS9 and CS7 can more favorably bind to BACE1 than BACE2. The hierarchical clustering analysis based on energetic contributions of residues suggests that BACE1 and BACE2 share the common hot interaction spots. The residue-based free-energy decomposition method was applied to compute the inhibitor-residue interaction spectrum, and the results recognize four common binding subpockets corresponding to the different groups of inhibitors, which can be used as efficient targets for designing highly selective inhibitors toward BACE1 and BACE2. Therefore, these results provide a useful molecular basis and dynamics information for development of highly selective inhibitors targeting BACE1 and BACE2.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357 China
| | - Jinan Wang
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Baohua Yin
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan 250355, China
| | - Laixue Pang
- School of Science, Shandong Jiaotong University, Jinan 250357 China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357 China
| | - Weiliang Zhu
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| |
Collapse
|
25
|
Zhong S, Huang K, Xiao Z, Sheng X, Li Y, Duan L. Binding Mechanism of Thrombin–Ligand Systems Investigated by a Polarized Protein-Specific Charge Force Field and Interaction Entropy Method. J Phys Chem B 2019; 123:8704-8716. [DOI: 10.1021/acs.jpcb.9b08064] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Susu Zhong
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Kaifang Huang
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Zhengrong Xiao
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Xiehuang Sheng
- College of Chemistry, Chemical Engineering and Materials Science, Shandong Normal University, Jinan 250014, China
| | - Yuchen Li
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Lili Duan
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| |
Collapse
|
26
|
chen J, Yin B, Pang L, Wang W, Zhang JZH, Zhu T. Binding modes and conformational changes of FK506-binding protein 51 induced by inhibitor bindings: insight into molecular mechanisms based on multiple simulation technologies. J Biomol Struct Dyn 2019; 38:2141-2155. [DOI: 10.1080/07391102.2019.1624616] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Jianzhong chen
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Baohua Yin
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, China
| | - Laixue Pang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - John Z. H. Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
| | - Tong Zhu
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
| |
Collapse
|
27
|
Yan F, Liu X, Zhang S, Su J, Zhang Q, Chen J. Electrostatic interaction-mediated conformational changes of adipocyte fatty acid binding protein probed by molecular dynamics simulation. J Biomol Struct Dyn 2018; 37:3583-3595. [DOI: 10.1080/07391102.2018.1520648] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Fangfang Yan
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Jing Su
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
28
|
Chen J, Peng C, Wang J, Zhu W. Exploring molecular mechanism of allosteric inhibitor to relieve drug resistance of multiple mutations in HIV-1 protease by enhanced conformational sampling. Proteins 2018; 86:1294-1305. [PMID: 30260044 DOI: 10.1002/prot.25610] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2018] [Revised: 09/02/2018] [Accepted: 09/23/2018] [Indexed: 12/12/2022]
Abstract
Recently, allosteric regulations of HIV-1 protease (PR) are suggested as a promising approach to relieve drug resistance of mutations toward inhibitors targeting the active site of PR. Replica-exchange molecular dynamics (REMD) simulations and normal mode analysis (NMA) are integrated to enhance conformational sampling of PR. Molecular mechanics generalized Born surface area (MM-GBSA) method was applied to calculate binding free energies of three inhibitors APV, DRV, and NIT to the wild-type (WT) and multidrug resistance (MDR) PRs. The results suggest that binding free energies of APV and DRV are decreased in the MDR PR relative to the WT PR, suggesting drug resistance of mutations on these two inhibitors. However, the binding ability of the allosteric inhibitor NIT is not impaired in the MDR PR. In addition, internal dynamics analysis based on REMD simulations proves that mutations hardly produce obvious effect on the conformation of the MDR PR in comparison to the WT PR. Scanning of hydrophobic contacts and hydrogen bond contacts of inhibitors with residues of PRs on the concatenated trajectories of REMD demonstrates that mutations change the symmetric interaction networks of APV and DRV with PR, but do not generate obvious influence on the asymmetric interaction network of NIT with PR. In summary, allosteric inhibitor NIT can adapt the MDR PR better than those inhibitors toward the active site of PR, thus allosteric inhibitors of PR may be a possible channel to overcome drug resistance of PR.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China.,Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Cheng Peng
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China.,School of Chemistry, University of Chinese Academy of Sciences, Beijing, China
| | - Jinan Wang
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Weiliang Zhu
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| |
Collapse
|
29
|
Su J, Liu X, Zhang S, Yan F, Zhang Q, Chen J. Insight into selective mechanism of class of I-BRD9 inhibitors toward BRD9 based on molecular dynamics simulations. Chem Biol Drug Des 2018; 93:163-176. [PMID: 30225973 DOI: 10.1111/cbdd.13398] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2018] [Revised: 08/10/2018] [Accepted: 09/02/2018] [Indexed: 02/06/2023]
Abstract
Recently, bromodomain-containing protein 9 (BRD9), 7 (BRD7), and 4 (BRD4) have been potential targets of anticancer drug design. Molecular dynamic simulations followed by molecular mechanics Poisson-Boltzmann surface area calculation were performed to study the selective mechanism of I-BRD9 inhibitor H1B and its derivatives N1D, TVU, and 5V2 toward BRD9 and BRD4. The rank of our calculated binding free energies agrees with that of the experimental data. The results show that binding free energy of H1B to BRD7 is slightly lower than that of H1B to BRD9, and all four inhibitors bind more tightly to BRD9 than to BRD4. Decomposition of binding free energies into individual residues implies that Ile164 and Asn211 in BRD7 and Ile53 and Asn100 in BRD9 play a significant role in the selectivity of H1B toward BRD7 and BRD9. Besides, several key residues Phe44, Ile53, Asn100, Thr104 in BRD9 and Pro82, Lys91, Asn140, Asp144 in BRD4 that are located in the ZA-loop and BC-loop provide significant contributions to binding selectivity of inhibitors to BRD9 and BRD4. This study is expected to provide important theoretical guidance for rational designs of highly selective inhibitors targeting BRD9 and BRD4.
Collapse
Affiliation(s)
- Jing Su
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Fangfang Yan
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
30
|
Yan F, Liu X, Zhang S, Su J, Zhang Q, Chen J. Molecular Dynamics Exploration of Selectivity of Dual Inhibitors 5M7, 65X, and 65Z toward Fatty Acid Binding Proteins 4 and 5. Int J Mol Sci 2018; 19:ijms19092496. [PMID: 30142969 PMCID: PMC6164837 DOI: 10.3390/ijms19092496] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2018] [Revised: 08/14/2018] [Accepted: 08/18/2018] [Indexed: 12/11/2022] Open
Abstract
Designing highly selective inhibitors of fatty acid binding proteins 4 and 5 (FABP4 and FABP5) is of importance for treatment of some diseases related with inflammation, metabolism, and tumor growth. In this study, molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were performed to probe binding selectivity of three inhibitors (5M7, 65X, and 65Z) to FABP4/FABP5 with Ki values of 0.022/0.50 μM, 0.011/0.086 μM, and 0.016/0.12 μM, respectively. The results not only suggest that all inhibitors associate more tightly with FABP4 than FABP5, but also prove that the main forces driving the selective bindings of inhibitors to FABP4 and FABP5 stem from the difference in the van der Waals interactions and polar interactions of inhibitors with two proteins. Meanwhile, a residue-based free energy decomposition method was applied to reveal molecular basis that inhibitors selectively interact with individual residues of two different proteins. The calculated results show that the binding difference of inhibitors to the residues (Phe16, Phe19), (Ala33, Gly36), (Phe57, Leu60), (Ala75, Ala78), (Arg126, Arg129), and (Tyr128, Tyr131) in (FABP4, FABP5) drive the selectivity of inhibitors toward FABP4 and FABP5. This study will provide great help for further design of effective drugs to protect against a series of metabolic diseases, arteriosclerosis, and inflammation.
Collapse
Affiliation(s)
- Fangfang Yan
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China.
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China.
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China.
| | - Jing Su
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China.
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China.
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China.
| |
Collapse
|
31
|
Chen J, Wang J, Pang L, Zhu W. Inhibiting mechanism of small molecule toward the p53-MDM2 interaction: A molecular dynamic exploration. Chem Biol Drug Des 2018; 92:1763-1777. [DOI: 10.1111/cbdd.13345] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2018] [Revised: 05/01/2018] [Accepted: 05/28/2018] [Indexed: 12/23/2022]
Affiliation(s)
- Jianzhong Chen
- School of Science; Shandong Jiaotong University; Jinan China
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center; Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| | - Jinan Wang
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center; Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| | - Laixue Pang
- School of Science; Shandong Jiaotong University; Jinan China
| | - Weiliang Zhu
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center; Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai China
| |
Collapse
|
32
|
Trypsin-Ligand binding affinities calculated using an effective interaction entropy method under polarized force field. Sci Rep 2017; 7:17708. [PMID: 29255159 PMCID: PMC5735144 DOI: 10.1038/s41598-017-17868-z] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2017] [Accepted: 12/01/2017] [Indexed: 12/22/2022] Open
Abstract
Molecular dynamics (MD) simulation in the explicit water is performed to study the interaction mechanism of trypsin-ligand binding under the AMBER force field and polarized protein-specific charge (PPC) force field combined the new developed highly efficient interaction entropy (IE) method for calculation of entropy change. And the detailed analysis and comparison of the results of MD simulation for two trypsin-ligand systems show that the root-mean-square deviation (RMSD) of backbone atoms, B-factor, intra-protein and protein-ligand hydrogen bonds are more stable under PPC force field than AMBER force field. Our results demonstrate that the IE method is superior than the traditional normal mode (Nmode) method in the calculation of entropy change and the calculated binding free energy under the PPC force field combined with the IE method is more close to the experimental value than other three combinations (AMBER-Nmode, AMBER-IE and PPC-Nmode). And three critical hydrogen bonds between trypsin and ligand are broken under AMBER force field. However, they are well preserved under PPC force field. Detailed binding interactions of ligands with trypsin are further analyzed. The present work demonstrates that the polarized force field combined the highly efficient IE method is critical in MD simulation and free energy calculation.
Collapse
|
33
|
Shi S, Zhang S, Zhang Q. Insight into binding mechanisms of inhibitors MKP56, MKP73, MKP86, and MKP97 to HIV-1 protease by using molecular dynamics simulation. J Biomol Struct Dyn 2017; 36:981-992. [PMID: 28279118 DOI: 10.1080/07391102.2017.1305296] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
HIV-1 protease (PR) has been a significant target for design of potent inhibitors curing acquired immunodeficiency syndrome. Molecular dynamics simulations coupled with molecular mechanics Poisson-Boltzmann surface area method were performed to study interaction modes of four inhibitors MKP56, MKP73, MKP86, and MKP97 with PR. The results suggest that the main force controlling interactions of inhibitors with PR should be contributed by van der Waals interactions between inhibitors and PR. The cross-correlation analyses based on MD trajectories show that inhibitor binding produces significant effect on the flap dynamics of PR. Hydrogen bond analyses indicate that inhibitors can form stable hydrogen bonding interactions with the residues from the catalytic strands of PR. The contributions of separate residues to inhibitor bindings are evaluated by using residue-based free energy decomposition method and the results demonstrate that the CH-π and CH-CH interactions between the hydrophobic groups of inhibitors with residues drive the associations of inhibitors with PR. We expect that this study can provide a significant theoretical aid for design of potent inhibitors targeting PR.
Collapse
Affiliation(s)
- Shuhua Shi
- a School of Science , Shandong Jianzhu University , Jinan 250101 , China
| | - Shaolong Zhang
- b College of Physics and Electronics , Shandong Normal University , Jinan 250014 , China
| | - Qinggang Zhang
- b College of Physics and Electronics , Shandong Normal University , Jinan 250014 , China
| |
Collapse
|
34
|
Chen J, Wang J, Zhu W. Mutation L1196M-induced conformational changes and the drug resistant mechanism of anaplastic lymphoma kinase studied by free energy perturbation and umbrella sampling. Phys Chem Chem Phys 2017; 19:30239-30248. [DOI: 10.1039/c7cp05418a] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Anaplastic lymphoma kinase (ALK) has been regarded as a promising drug target in the treatment of tumors and the mutation L1196M induces different levels of drug resistance toward the existing inhibitors.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University
- Jinan
- China
| | - Jinan Wang
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences
- Shanghai
- China
| | - Weiliang Zhu
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences
- Shanghai
- China
| |
Collapse
|
35
|
Hu G, Ma A, Dou X, Zhao L, Wang J. Computational Studies of a Mechanism for Binding and Drug Resistance in the Wild Type and Four Mutations of HIV-1 Protease with a GRL-0519 Inhibitor. Int J Mol Sci 2016; 17:E819. [PMID: 27240358 PMCID: PMC4926353 DOI: 10.3390/ijms17060819] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2016] [Revised: 05/16/2016] [Accepted: 05/16/2016] [Indexed: 12/22/2022] Open
Abstract
Drug resistance of mutations in HIV-1 protease (PR) is the most severe challenge to the long-term efficacy of HIV-1 PR inhibitor in highly active antiretroviral therapy. To elucidate the molecular mechanism of drug resistance associated with mutations (D30N, I50V, I54M, and V82A) and inhibitor (GRL-0519) complexes, we have performed five molecular dynamics (MD) simulations and calculated the binding free energies using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. The ranking of calculated binding free energies is in accordance with the experimental data. The free energy spectra of each residue and inhibitor interaction for all complexes show a similar binding model. Analysis based on the MD trajectories and contribution of each residues show that groups R2 and R3 mainly contribute van der Waals energies, while groups R1 and R4 contribute electrostatic interaction by hydrogen bonds. The drug resistance of D30N can be attributed to the decline in binding affinity of residues 28 and 29. The size of Val50 is smaller than Ile50 causes the residue to move, especially in chain A. The stable hydrophobic core, including the side chain of Ile54 in the wild type (WT) complex, became unstable in I54M because the side chain of Met54 is flexible with two alternative conformations. The binding affinity of Ala82 in V82A decreases relative to Val82 in WT. The present study could provide important guidance for the design of a potent new drug resisting the mutation inhibitors.
Collapse
Affiliation(s)
- Guodong Hu
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University, Dezhou 253023, China.
| | - Aijing Ma
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University, Dezhou 253023, China.
| | - Xianghua Dou
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University, Dezhou 253023, China.
| | - Liling Zhao
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University, Dezhou 253023, China.
| | - Jihua Wang
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University, Dezhou 253023, China.
| |
Collapse
|
36
|
Chen J. Drug resistance mechanisms of three mutations V32I, I47V and V82I in HIV-1 protease toward inhibitors probed by molecular dynamics simulations and binding free energy predictions. RSC Adv 2016. [DOI: 10.1039/c6ra09201b] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
Molecular dynamics simulation and binding free energy calculations were used to probe drug resistance of HIV-1 protease mutations toward inhibitors.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science
- Shandong Jiaotong University
- Jinan 250357
- China
| |
Collapse
|
37
|
Chen J, Wang J, Zhang Q, Chen K, Zhu W. Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation. Sci Rep 2015; 5:17421. [PMID: 26616018 PMCID: PMC4663504 DOI: 10.1038/srep17421] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2015] [Accepted: 10/29/2015] [Indexed: 12/13/2022] Open
Abstract
Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250014, China
| | - Jinan Wang
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| | - Qinggang Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Kaixian Chen
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| | - Weiliang Zhu
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| |
Collapse
|
38
|
Shi S, Zhang S, Zhang Q. Probing Difference in Binding Modes of Inhibitors to MDMX by Molecular Dynamics Simulations and Different Free Energy Methods. PLoS One 2015; 10:e0141409. [PMID: 26513747 PMCID: PMC4625964 DOI: 10.1371/journal.pone.0141409] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2015] [Accepted: 10/08/2015] [Indexed: 01/20/2023] Open
Abstract
The p53-MDMX interaction has attracted extensive attention of anti-cancer drug development in recent years. This current work adopted molecular dynamics (MD) simulations and cross-correlation analysis to investigate conformation changes of MDMX caused by inhibitor bindings. The obtained information indicates that the binding cleft of MDMX undergoes a large conformational change and the dynamic behavior of residues obviously change by the presence of different structural inhibitors. Two different methods of binding free energy predictions were employed to carry out a comparable insight into binding mechanisms of four inhibitors PMI, pDI, WK23 and WW8 to MDMX. The data show that the main factor controlling the inhibitor bindings to MDMX arises from van der Waals interactions. The binding free energies were further divided into contribution of each residue and the derived information gives a conclusion that the hydrophobic interactions, such as CH-CH, CH-π and π-π interactions, are responsible for the inhibitor associations with MDMX.
Collapse
Affiliation(s)
- Shuhua Shi
- School of Science, Shandong Jianzhu University, Jinan, China
- * E-mail: ;
| | - Shaolong Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Qinggang Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, China
| |
Collapse
|
39
|
Cheng W, Liang Z, Wang W, Yi C, Li H, Zhang S, Zhang Q. Insight into binding modes of p53 and inhibitors to MDM2 based on molecular dynamic simulations and principal component analysis. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1087598] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
40
|
Chen J, Wang X, Zhu T, Zhang Q, Zhang JZH. A Comparative Insight into Amprenavir Resistance of Mutations V32I, G48V, I50V, I54V, and I84V in HIV-1 Protease Based on Thermodynamic Integration and MM-PBSA Methods. J Chem Inf Model 2015; 55:1903-13. [DOI: 10.1021/acs.jcim.5b00173] [Citation(s) in RCA: 98] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Jianzhong Chen
- School
of Science, Shandong Jiaotong University, Jinan 250357 China
| | - Xingyu Wang
- NYU−ECNU
Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Tong Zhu
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
| | - Qinggang Zhang
- Collage
of Physics and Electronic Science, Shandong Normal University, Jinan 250014, China
| | - John Z. H. Zhang
- NYU−ECNU
Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
| |
Collapse
|
41
|
Chen J, Liang Z, Wang W, Yi C, Zhang S, Zhang Q. Revealing origin of decrease in potency of darunavir and amprenavir against HIV-2 relative to HIV-1 protease by molecular dynamics simulations. Sci Rep 2014; 4:6872. [PMID: 25362963 PMCID: PMC4217091 DOI: 10.1038/srep06872] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2014] [Accepted: 10/13/2014] [Indexed: 12/30/2022] Open
Abstract
Clinical inhibitors Darunavir (DRV) and Amprenavir (APV) are less effective on HIV-2 protease (PR2) than on HIV-1 protease (PR1). To identify molecular basis associated with the lower inhibition, molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations were performed to investigate the effectiveness of the PR1 inhibitors DRV and APV against PR1/PR2. The rank of predicted binding free energies agrees with the experimental determined one. Moreover, our results show that two inhibitors bind less strongly to PR2 than to PR1, again in agreement with the experimental findings. The decrease in binding free energies for PR2 relative to PR1 is found to arise from the reduction of the van der Waals interactions induced by the structural adjustment of the triple mutant V32I, I47V and V82I. This result is further supported by the difference between the van der Waals interactions of inhibitors with each residue in PR2 and in PR1. The results from the principle component analysis suggest that inhibitor binding tends to make the flaps of PR2 close and the one of PR1 open. We expect that this study can theoretically provide significant guidance and dynamics information for the design of potent dual inhibitors targeting PR1/PR2.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Zhiqiang Liang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Changhong Yi
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Shaolong Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Qinggang Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| |
Collapse
|
42
|
Shi S, Zhang S, Zhang Q. Insight into the interaction mechanism of inhibitors P4 and WK23 with MDM2 based on molecular dynamics simulation and different free energy methods. COMPUT THEOR CHEM 2014. [DOI: 10.1016/j.comptc.2014.07.002] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
43
|
Kudryavtsev KV, Shulga DA, Chupakhin VI, Sinauridze EI, Ataullakhanov FI, Vatsadze SZ. Synthesis of novel bridged dinitrogen heterocycles and their evaluation as potential fragments for the design of biologically active compounds. Tetrahedron 2014. [DOI: 10.1016/j.tet.2014.09.009] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
44
|
A study of the interaction between HIV-1 protease and C 2-symmetric inhibitors by computational methods. J Mol Model 2014; 20:2369. [DOI: 10.1007/s00894-014-2369-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2014] [Accepted: 06/29/2014] [Indexed: 12/14/2022]
|
45
|
Chen J, Wang J, Zhu W. Binding modes of three inhibitors 8CA, F8A and I4A to A-FABP studied based on molecular dynamics simulation. PLoS One 2014; 9:e99862. [PMID: 24918907 PMCID: PMC4053400 DOI: 10.1371/journal.pone.0099862] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2014] [Accepted: 05/16/2014] [Indexed: 12/21/2022] Open
Abstract
Adipocyte fatty-acid binding protein (A-FABP) is an important target of drug designs treating some diseases related to lipid-mediated biology. Molecular dynamics (MD) simulations coupled with solvated interaction energy method (SIE) were carried out to study the binding modes of three inhibitors 8CA, F8A and I4A to A-FABP. The rank of our predicted binding affinities is in accordance with experimental data. The results show that the substitution in the position 5 of N-benzyl and the seven-membered ring of N-benzyl-indole carboxylic acids strengthen the I4A binding, while the substitution in the position 2 of N-benzyl weakens the F8A binding. Computational alanine scanning and dynamics analyses were performed and the results suggest that the polar interactions of the positively charged residue R126 with the three inhibitors provide a significant contribution to inhibitor bindings. This polar interaction induces the disappearance of the correlated motion of the C terminus of A-FABP relative to the N terminus and favors the stability of the binding complex. This study is helpful for the rational design of potent inhibitors within the fields of metabolic disease, inflammation and atherosclerosis.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
- * E-mail: (JC); (WZ)
| | - Jinan Wang
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Weiliang Zhu
- Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- * E-mail: (JC); (WZ)
| |
Collapse
|
46
|
Ermakova E, Kurbanov R. Effect of ligand binding on the dynamics of trypsin. Comparison of different approaches. J Mol Graph Model 2014; 49:99-109. [PMID: 24642055 DOI: 10.1016/j.jmgm.2014.02.001] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2014] [Revised: 02/07/2014] [Accepted: 02/08/2014] [Indexed: 11/17/2022]
Abstract
The intramolecular signal transduction induced by the binding of ligands to trypsin was investigated by molecular dynamics simulations. Ligand binding changes the residue-residue interaction energies and suppresses the mobility of loops that are in direct contact with the ligand. The reduced mobility of these loops results in the altered flexibility of the nearby loops and thereby transmits the information from ligand binding site to the remote sites. The analysis of the flexibility of all residues confirmed the coupling between loops L1 (185-188) and L2 (221-224) and the residues in the active center. The significance of S1 pocket residues for the signal transduction from the active center to the substrate-binding site was confirmed by the dynamical network and covariance matrix analyses. Gaussian network model and principal component analysis demonstrated that the active center residues had zero amplitude in the slowest fluctuations acting as hinges or anchors. Overall, our results provide a new insight into protein-ligand interactions and show how the allosteric signaling may occur.
Collapse
Affiliation(s)
- Elena Ermakova
- Kazan Institute of Biochemistry and Biophysics RAS, P.O. Box 30, Kazan 420111, Russia.
| | - Rauf Kurbanov
- Kazan Institute of Biochemistry and Biophysics RAS, P.O. Box 30, Kazan 420111, Russia
| |
Collapse
|
47
|
Gu Y, Wang W, Zhu X, Dong K. Molecular dynamic simulations reveal the mechanism of binding between xanthine inhibitors and DPP-4. J Mol Model 2014; 20:2075. [PMID: 24481594 DOI: 10.1007/s00894-014-2075-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2013] [Accepted: 11/12/2013] [Indexed: 11/24/2022]
Abstract
We apply molecular docking, molecular dynamics (MD) simulation, and binding free energy calculation to investigate and reveal the binding mechanism between five xanthine inhibitors and DPP-4. The electrostatic and van der Waals interactions of the five inhibitors with DPP-4 are analyzed and discussed. The computed binding free energies using MM-PBSA method are in qualitatively agreement with experimental inhibitory potency of five inhibitors. The hydrogen bonds of inhibitors with Ser630 and Asp663 can stabilize the inhibitors in binding sites. The van der Waals interactions, especially the key contacts with His740, Asn710, Trp629, and Tyr666 have larger contributions to the binding free energy and play important roles in distinguishing the variant bioactivity of five inhibitors.
Collapse
Affiliation(s)
- Yongliang Gu
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemistry and Chemical Engineering, Nanjing University of Technology, Nanjing, 210009, China
| | | | | | | |
Collapse
|
48
|
Chen J, Wang J, Zhu W, Li G. A computational analysis of binding modes and conformation changes of MDM2 induced by p53 and inhibitor bindings. J Comput Aided Mol Des 2013; 27:965-74. [PMID: 24264557 DOI: 10.1007/s10822-013-9693-z] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2013] [Accepted: 11/17/2013] [Indexed: 10/26/2022]
Abstract
Molecular dynamics (MD) simulations followed by principal component analysis were performed to study the conformational change of MDM2 induced by p53 and two inhibitor (P4 and MI63a) bindings. The results show that the hydrophobic cleft of MDM2 is very flexible and adaptive to different structural binding partners. The cleft tends to become wider and more stable as MDM2 binds to the three binding partners, while unbound MDM2 shows a narrower and pretty flexible cleft, which agrees with recent experimental data and theoretical studies. It was also found that the binding of P4 and p53 stabilizes the motion of the loop L2 linking the helix α2 and β strand (β3), but the presence of MI63a makes the motion of L2 disordered. In addition, the binding free energies of the three partners to MDM2 were calculated using molecular mechanics generalized Born surface area to explain the binding modes of these three partners to MDM2. This study will be helpful not only for better understanding the functional, concerted motion of MDM2, but also for the rational design of potent anticancer drugs targeting the p53-MDM2 interaction.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250014, China,
| | | | | | | |
Collapse
|
49
|
Ayan S, Dogan Ö, Ivantcova PM, Datsuk NG, Shulga DA, Chupakhin VI, Zabolotnev DV, Kudryavtsev KV. Asymmetric synthesis and molecular docking study of enantiomerically pure pyrrolidine derivatives with potential antithrombin activity. ACTA ACUST UNITED AC 2013. [DOI: 10.1016/j.tetasy.2013.05.023] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
50
|
Hu G, Wang J. Ligand selectivity of estrogen receptors by a molecular dynamics study. Eur J Med Chem 2013; 74:726-35. [PMID: 23694906 DOI: 10.1016/j.ejmech.2013.04.049] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2012] [Revised: 03/09/2013] [Accepted: 04/22/2013] [Indexed: 02/07/2023]
Abstract
Estrogen receptors α (ERα) and β (ERβ) have different physiological functions and expression levels in different tissues. ERα and ERβ are highly homologous and have only two residue substitutions in the binding pocket. This high similarity at the active site stimulates the interests for discovering subtype selective ligands. In this study, molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method have been carried out to analyze the basis of selectivity of three ligands (659, 818 and 041). The calculated binding free energies show that all the ligands bind more tightly to ERβ than to ERα. The dominant free energy components of selectivity for 659 are similar to that for 041, but different from that for 818. The decompositions of free energy contributions and structural analysis imply that there are eight residues primarily contributing to the selectivity for 659, five residues for 041, as well as two residues for 818. The structural analysis implies that two residue substitutions in binding packet cause the position of 659 in ERβ-659 complex to shift relative to that in ERα-659 complex and also cause the conformational changes of other residues in the binding pocket. The higher selectivity for 041 is mainly caused by three residues, Ile373 (Met421), His475 (His524) and Leu476 (Leu525).
Collapse
Affiliation(s)
- Guodong Hu
- Shandong Provincial Key Laboratory of Functional Macromolecular Biophysics, Dezhou University, Dezhou, Shandong 253023, China; Department of Physics, Dezhou University, Dezhou, Shandong 253023, China.
| | - Jihua Wang
- Shandong Provincial Key Laboratory of Functional Macromolecular Biophysics, Dezhou University, Dezhou, Shandong 253023, China; Department of Physics, Dezhou University, Dezhou, Shandong 253023, China.
| |
Collapse
|