1
|
Jernigan RL, Khade P, Kumar A, Kloczkowski A. Using Surface Hydrophobicity Together with Empirical Potentials to Identify Protein-Protein Binding Sites: Application to the Interactions of E-cadherins. Methods Mol Biol 2022; 2340:41-50. [PMID: 35167069 PMCID: PMC9131873 DOI: 10.1007/978-1-0716-1546-1_3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Studying the interactions within protein structures can inform about the details of how proteins of various types interact and aggregate. Empirical contact potentials have proven to be extremely important in the evaluation of individual modeled protein structures, but have found few applications to protein-protein interactions. In part, this is caused by a lack of properly formulated potentials with a proper reference state. Since the comparisons are made between different bound structures, the proper reference state should take into account other contacts. Therefore, a preferred reference state should be defined with respect to a given residue type interacting with an average residue instead of interacting with solvent as typically is used in derivation of statistical contact potentials. Here, a two-stage procedure for generating and evaluating interacting protein pairs is described, and an example of E-cadherin interactions is shown.
Collapse
Affiliation(s)
- Robert L Jernigan
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, IA, USA.
| | - Pranav Khade
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, IA, USA
| | - Ambuj Kumar
- Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, IA, USA
| | - Andrzej Kloczkowski
- Battelle Center for Mathematical Medicine, Nationwide Children's Hospital, Columbus, OH, USA
- Department of Pediatrics, The Ohio State University College of Medicine, Columbus, OH, USA
| |
Collapse
|
2
|
Liu Y, Tu G, Lai X, Kuang B, Li S. Exploring ligand dissociation pathways from aminopeptidase N using random acceleration molecular dynamics simulation. J Mol Model 2016; 22:236. [PMID: 27624165 DOI: 10.1007/s00894-016-3105-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2016] [Accepted: 09/01/2016] [Indexed: 11/30/2022]
Abstract
Aminopeptidase N (APN) is a zinc-dependent ectopeptidase involved in cell proliferation, secretion, invasion, and angiogenesis, and is widely recognized as an important cancer target. However, the mechanisms whereby ligands leave the active site of APN remain unknown. Investigating ligand dissociation processes is quite difficult, both in classical simulation methods and in experimental approaches. In this study, random acceleration molecular dynamics (RAMD) simulation was used to investigate the potential dissociation pathways of ligand from APN. The results revealed three pathways (channels A, B and C) for ligand release. Channel A, which matches the hypothetical channel region, was the most preferred region for bestatin to dissociate from the enzyme, and is probably the major channel for the inner bound ligand. In addition, two alternative channels (channels B and C) were shown to be possible pathways for ligand egression. Meanwhile, we identified key residues controlling the dynamic features of APN channels. Identification of the dissociation routes will provide further mechanistic insights into APN, which will benefit the development of more promising APN inhibitors. Graphical Abstract The release pathways of bestatin inside active site of aminopeptidase N were simulated using RAMD simulation.
Collapse
Affiliation(s)
- Ya Liu
- Department of Medicinal Chemistry, School of Pharmaceutical Science, NanChang University, 461, BaYi Road, NanChang, 330006, People's Republic of China
| | - GuoGang Tu
- Department of Medicinal Chemistry, School of Pharmaceutical Science, NanChang University, 461, BaYi Road, NanChang, 330006, People's Republic of China.
| | - XiaoPing Lai
- Department of Medicinal Chemistry, School of Pharmaceutical Science, NanChang University, 461, BaYi Road, NanChang, 330006, People's Republic of China
| | - BinHai Kuang
- Department of Medicinal Chemistry, School of Pharmaceutical Science, NanChang University, 461, BaYi Road, NanChang, 330006, People's Republic of China
| | - ShaoHua Li
- Department of Medicinal Chemistry, School of Pharmaceutical Science, NanChang University, 461, BaYi Road, NanChang, 330006, People's Republic of China
| |
Collapse
|
3
|
Shadrina MS, English AM, Peslherbe GH. Benchmarking Rapid TLES Simulations of Gas Diffusion in Proteins: Mapping O2 Migration and Escape in Myoglobin as a Case Study. J Chem Theory Comput 2016; 12:2038-46. [PMID: 26938707 DOI: 10.1021/acs.jctc.5b01132] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Standard molecular dynamics (MD) simulations of gas diffusion consume considerable computational time and resources even for small proteins. To combat this, temperature-controlled locally enhanced sampling (TLES) examines multiple diffusion trajectories per simulation by accommodating multiple noninteracting copies of a gas molecule that diffuse independently, while the protein and water molecules experience an average interaction from all copies. Furthermore, gas migration within a protein matrix can be accelerated without altering protein dynamics by increasing the effective temperature of the TLES copies. These features of TLES enable rapid simulations of gas diffusion within a protein matrix at significantly reduced (∼98%) computational cost. However, the results of TLES and standard MD simulations have not been systematically compared, which limits the adoption of the TLES approach. We address this drawback here by benchmarking TLES against standard MD in the simulation of O2 diffusion in myoglobin (Mb) as a case study since this model system has been extensively characterized. We find that 2 ns TLES and 108 ns standard simulations map the same network of diffusion tunnels in Mb and uncover the same docking sites, barriers, and escape portals. We further discuss the influence of simulation time as well as the number of independent simulations on the O2 population density within the diffusion tunnels and on the sampling of Mb's conformational space as revealed by principal component analysis. Overall, our comprehensive benchmarking reveals that TLES is an appropriate and robust tool for the rapid mapping of gas diffusion in proteins when the kinetic data provided by standard MD are not required. Furthermore, TLES provides explicit ligand diffusion pathways, unlike most rapid methods.
Collapse
Affiliation(s)
- Maria S Shadrina
- Centre for Research in Molecular Modeling (CERMM) and Department of Chemistry and Biochemistry, Concordia University , 7141 Sherbrooke Street West, Montréal, Québec, Canada H4B 1R6
| | - Ann M English
- Centre for Research in Molecular Modeling (CERMM) and Department of Chemistry and Biochemistry, Concordia University , 7141 Sherbrooke Street West, Montréal, Québec, Canada H4B 1R6
| | - Gilles H Peslherbe
- Centre for Research in Molecular Modeling (CERMM) and Department of Chemistry and Biochemistry, Concordia University , 7141 Sherbrooke Street West, Montréal, Québec, Canada H4B 1R6
| |
Collapse
|
4
|
Gee LB, Leontyev I, Stuchebrukhov A, Scott AD, Pelmenschikov V, Cramer SP. Docking and migration of carbon monoxide in nitrogenase: the case for gated pockets from infrared spectroscopy and molecular dynamics. Biochemistry 2015; 54:3314-9. [PMID: 25919807 DOI: 10.1021/acs.biochem.5b00216] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Evidence of a CO docking site near the FeMo cofactor in nitrogenase has been obtained by Fourier transform infrared spectroscopy-monitored low-temperature photolysis. We investigated the possible migration paths for CO from this docking site using molecular dynamics calculations. The simulations support the notion of a gas channel with multiple internal pockets from the active site to the protein exterior. Travel between pockets is gated by the motion of protein residues. Implications for the mechanism of nitrogenase reactions with CO and N2 are discussed.
Collapse
Affiliation(s)
- Leland B Gee
- †Department of Chemistry, University of California, Davis, California 95616, United States
| | - Igor Leontyev
- §InterX Inc., Berkeley, California 94710, United States
| | - Alexei Stuchebrukhov
- †Department of Chemistry, University of California, Davis, California 95616, United States
| | - Aubrey D Scott
- †Department of Chemistry, University of California, Davis, California 95616, United States
| | | | - Stephen P Cramer
- †Department of Chemistry, University of California, Davis, California 95616, United States.,‡Physical Biosciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
5
|
Na H, Song G. Quantitative delineation of how breathing motions open ligand migration channels in myoglobin and its mutants. Proteins 2015; 83:757-70. [PMID: 25645487 DOI: 10.1002/prot.24770] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2014] [Revised: 01/01/2015] [Accepted: 01/14/2015] [Indexed: 11/09/2022]
Abstract
Ligand migration and binding are central to the biological functions of many proteins such as myoglobin (Mb) and it is widely thought that protein breathing motions open up ligand channels dynamically. However, how a protein exerts its control over the opening and closing of these channels through its intrinsic dynamics is not fully understood. Specifically, a quantitative delineation of the breathing motions that are needed to open ligand channels is lacking. In this work, we present and apply a novel normal mode-based method to quantitatively delineate what and how breathing motions open ligand migration channels in Mb and its mutants. The motivation behind this work springs from the observation that normal mode motions are closely linked to the breathing motions that are thought to open ligand migration channels. In addition, the method provides a direct and detailed depiction of the motions of each and every residue that lines a channel and can identify key residues that play a dominating role in regulating the channel. The all-atom model and the full force-field employed in the method provide a realistic energetics on the work cost required to open a channel, and as a result, the method can be used to efficiently study the effects of mutations on ligand migration channels and on ligand entry rates. Our results on Mb and its mutants are in excellent agreement with MD simulation results and experimentally determined ligand entry rates.
Collapse
Affiliation(s)
- Hyuntae Na
- Department of Computer Science, Iowa State University, Ames, Iowa, 50011
| | | |
Collapse
|
6
|
Hosseini A, Brouk M, Lucas MF, Glaser F, Fishman A, Guallar V. Atomic picture of ligand migration in toluene 4-monooxygenase. J Phys Chem B 2014; 119:671-8. [PMID: 24798294 DOI: 10.1021/jp502509a] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Computational modeling combined with mutational and activity assays was used to underline the substrate migration pathways in toluene 4-monooxygenase, a member of the important family of bacterial multicomponent monooxygenases (BMMs). In all structurally defined BMM hydroxylases, several hydrophobic cavities in the α-subunit map a preserved path from the protein surface to the diiron active site. Our results confirm the presence of two pathways by which different aromatic molecules can enter/escape the active site. While the substrate is observed to enter from both channels, the more hydrophilic product is withdrawn mainly from the shorter channel ending at residues D285 and E214. The long channel ends in the vicinity of S395, whose variants have been seen to affect activity and specificity. These mutational effects are clearly reproduced and rationalized by the in silico studies. Furthermore, the combined computational and experimental results highlight the importance of residue F269, which is located at the intersection of the two channels.
Collapse
Affiliation(s)
- Ali Hosseini
- Department of Life Sciences, Barcelona Supercomputing Center , Nexus II Building, 08034 Barcelona, Spain
| | | | | | | | | | | |
Collapse
|
7
|
Abstract
![]()
Myoglobin
(Mb) binds diatomic ligands, like O2, CO,
and NO, in a cavity that is only transiently accessible. Crystallography
and molecular simulations show that the ligands can migrate through
an extensive network of transiently connected cavities but disagree
on the locations and occupancy of internal hydration sites. Here,
we use water 2H and 17O magnetic relaxation
dispersion (MRD) to characterize the internal water molecules in Mb
under physiological conditions. We find that equine carbonmonoxy Mb
contains 4.5 ± 1.0 ordered internal water molecules with a mean
survival time of 5.6 ± 0.5 μs at 25 °C. The likely
locations of these water molecules are the four polar hydration sites,
including one of the xenon-binding cavities, that are fully occupied
in all high-resolution crystal structures of equine Mb. The finding
that water escapes from these sites, located 17–31 Å apart
in the protein, on the same μs time scale suggests a global
exchange mechanism. We propose that this mechanism involves transient
penetration of the protein by H-bonded water chains. Such a mechanism
could play a functional role by eliminating trapped ligands. In addition,
the MRD results indicate that 2 or 3 of the 11 histidine residues
of equine Mb undergo intramolecular hydrogen exchange on a μs
time scale.
Collapse
Affiliation(s)
- Shuji Kaieda
- Department of Biophysical Chemistry, Lund University , P.O. Box 124, SE-22100 Lund, Sweden
| | | |
Collapse
|
8
|
Kalyaanamoorthy S, Chen YPP. Ligand release mechanisms and channels in histone deacetylases. J Comput Chem 2013; 34:2270-83. [PMID: 23893931 DOI: 10.1002/jcc.23390] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2013] [Revised: 05/31/2013] [Accepted: 06/28/2013] [Indexed: 01/08/2023]
Abstract
Exploring the molecular channels of class I histone deacetylases (HDACs) with buried active sites are important to understand their structures and functionalities. In this work, we perform hybrid classical molecular dynamics and random acceleration molecular dynamics simulations to explore the B3N [i.e., (4-(dimethylamino)N-[7(hydroxyamino)-7-oxoheptyle] benzamide)] exit channels in the x-ray crystal structures of HDAC3 and HDAC8 enzymes. Our simulations identify B3N release through four different channels in HDAC3 (denoted as A1, A2, B1, and B2) and HDAC8 (referred as A1, B1, B2, and B3) enzymes, among which egression through channel A1 is more predominant in both the enzymes. This mechanism is similar to ligand release in HDAC1 and HDAC2 described in our previous study and can be the fingerprint ligand release mechanisms in class I HDACs. Ligand release events through B channels, on the other hand, are different among HDAC3 and HDAC8, highlighting the significances of substituted residues in controlling the access to these channels This study reveals a novel aromatic gating mechanism elicited by TYR154-TRP141-TYR111 that controls the B3N access to all the B channels in HDAC8. The TRP141 in HDAC8 is substituted by LEU133 in HDAC3, which do not hinder the access to B channels in HDAC3. However, two hydrogen bonded barricades formed as ARG28-GLY297-GLY295-GLY131 and TRP129-ARG28-ALA130-LEU29-TRP129 obstruct the B3N from exploring the B channels in HDAC3. The structural and dynamical characterizations of molecular channels and ligand unbinding mechanisms reported in this study provide novel structural insights and atomic level perspectives on HDAC3 and HDAC8 enzymes, thereby potentially aiding in the design of more specific HDAC inhibitors.
Collapse
Affiliation(s)
- Subha Kalyaanamoorthy
- Department of Computer Science and Computer Engineering, La Trobe University, Melbourne, Australia; CSIRO Ecosystem Sciences, Canberra, Australia
| | | |
Collapse
|
9
|
Chen HY, Chang SS, Chan YC, Chen CYC. Discovery of novel insomnia leads from screening traditional Chinese medicine database. J Biomol Struct Dyn 2013; 32:776-91. [DOI: 10.1080/07391102.2013.790849] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
10
|
Abbruzzetti S, Spyrakis F, Bidon-Chanal A, Luque FJ, Viappiani C. Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations. Phys Chem Chem Phys 2013; 15:10686-701. [PMID: 23733145 DOI: 10.1039/c3cp51149a] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.
Collapse
Affiliation(s)
- Stefania Abbruzzetti
- Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Parma, viale delle Scienze 7A, 43124, Parma, Italy
| | | | | | | | | |
Collapse
|
11
|
IeongTou W, Chang SS, Wu D, Lai TW, Wang YT, Hsu CY, Yu-ChianChen C. Molecular level activation insights from a NR2A/NR2B agonist. J Biomol Struct Dyn 2013; 32:683-93. [DOI: 10.1080/07391102.2013.787371] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
12
|
Hung IC, Chang SS, Chang PC, Lee CC, Chen CYC. Memory enhancement by traditional Chinese medicine? J Biomol Struct Dyn 2012; 31:1411-39. [PMID: 23249175 DOI: 10.1080/07391102.2012.741052] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
Cognitive repair by insulin-like growth factor-I (IGF-I) through activation of insulin-like growth factor-I receptor (IGF-IR) is well established, but not used for clinical therapy due to its link to cancer. We hypothesize that IGF-IR activation rather than IGF-I per se may be essential for cognitive repair and attempted to identify ligands from traditional Chinese medicine (TCM) with drug-like potential towards IGF-IR. TCM ligands, 3-(2-carboxyphenyl)-4(3H)-quinazolinone from Isatisin digotica, (+)-N-methyllaurotetanine from Lindera aggregate, and (+)-1(R)-Coclaurine from Nelumbonucifera Gaertn, exhibited high binding affinities and good blood brain barrier (BBB) penetration crucial for accessing IGF-IR. Stable complex formation of the candidates was observed during molecular dynamics (MD) simulation. Interactions with Leu975 and Gly1055 or Asp1056 were important for ligand binding. Amino acid distance analysis revealed residues 974/975, 984-986, 996-1006, 1040-1056, and 1122-1135 as "hotspots" for ligand binding in IGF-IR. Versatile entry pathways for the TCM candidates suggest high accessibility to the binding site. Blockage of the binding site opening by the TCM candidates limits binding site access by other compounds. Multiple linear regression (R² = 0.9715), support vector machine (R² = 0.9084), Bayesian network (R² =0.8233) comparative molecular field analysis (CoMFA, R² = 0.9941), and comparative molecular similarity indices analysis (CoMSIA, R² = 0.9877) models consistently suggest that the TCM candidates might exert bioactivity on IGF-IR. Contour of representative MD conformations to CoMFA and CoMSIA maps exhibits similar results. Properties including BBB passage, evidence of ability to form stable complexes with IGF-IR by MD simulation, and predicted bioactivity suggest that the TCM candidates have drug-like properties and might have potential as cognitive-enhancing drugs.
Collapse
Affiliation(s)
- I-Chi Hung
- a School of Pharmacy, China Medical University , Taichung , 40402 , Taiwan
| | | | | | | | | |
Collapse
|
13
|
Tsou YA, Chen KC, Lin HC, Chang SS, Chen CYC. Uroporphyrinogen decarboxylase as a potential target for specific components of traditional Chinese medicine: a virtual screening and molecular dynamics study. PLoS One 2012; 7:e50087. [PMID: 23209648 PMCID: PMC3510221 DOI: 10.1371/journal.pone.0050087] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2012] [Accepted: 10/15/2012] [Indexed: 11/23/2022] Open
Abstract
Uroporphyrinogen decarboxylase (UROD) has been suggested as a protectant against radiation for head and neck cancer (HNC). In this study, we employed traditional Chinese medicine (TCM) compounds from TCM Database@Taiwan (http://tcm.cmu.edu.tw/) to screen for drug-like candidates with potential UROD inhibition characteristics using virtual screening techniques. Isopraeroside IV, scopolin, and nodakenin exhibited the highest Dock Scores, and were predicted to have good Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) properties. Two common moieties, 2H-chromen-2-one and glucoside, were observed among the top TCM candidates. Cross comparison of the docking poses indicated that candidates formed stable interactions with key binding and catalytic residues of UROD through these two moieties. The 2H-chromen-2-one moiety enabled pi-cation interactions with Arg37 and H-bonds with Tyr164. The glucoside moiety was involved in forming H-bonds with Arg37 and Asp86. From our computational results, we propose isopraeroside IV, scopolin, and nodakenin as ligands that might exhibit drug-like inhibitory effects on UROD. The glucoside and 2H-chromen-2-one moieties may potentially be used for designing inhibitors of UROD.
Collapse
Affiliation(s)
- Yung-An Tsou
- Department of Otolaryngology, China Medical University Hospital, Taichung, Taiwan
- School of Medicine, College of Medicine, China Medical University, Taichung, Taiwan
| | - Kuan-Chung Chen
- Graduate Institute of Pharmaceutical Chemistry, China Medical University, Taichung, Taiwan
| | - Hung-Che Lin
- School of Pharmacy, China Medical University, Taichung, Taiwan
| | - Su-Sen Chang
- Department of Medical Research, China Medical University Hospital, Taichung, Taiwan
| | - Calvin Yu-Chian Chen
- Department of Medical Research, China Medical University Hospital, Taichung, Taiwan
- Laboratory of Computational and Systems Biology, China Medical University, Taichung, Taiwan
- Department of Biotechnology, Asia University, Taichung, Taiwan
- Department of Biomedical Informatics, Asia University, Taichung, Taiwan
- China Medical University Beigang Hospital, Yunlin, Taiwan
- * E-mail:
| |
Collapse
|
14
|
Chen KC, Jian YR, Sun MF, Chang TT, Lee CC, Chen CYC. Investigation of silent information regulator 1 (Sirt1) agonists from Traditional Chinese Medicine. J Biomol Struct Dyn 2012; 31:1207-18. [PMID: 23075283 DOI: 10.1080/07391102.2012.726191] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
Silent information regulator 1 (Sirt1), a class III nicotinamide adenine dinucleotide dependent histone deacetylases, is important in cardioprotection, neuroprotection, metabolic disease, calorie restriction, and diseases associated with aging. Traditional Chinese Medicine (TCM) compounds from TCM Database@Taiwan ( http://tcm.cmu.edu.tw/ ) were employed for screening potent Sirt1 agonists, and molecular dynamics (MD) simulation was implemented to simulate ligand optimum docking poses and protein structure under dynamic conditions. TCM compounds such as (S)-tryptophan-betaxanthin, 5-O-feruloylquinic acid, and RosA exhibited good binding affinity across different computational methods, and their drug-like potential were validated by MD simulation. Docking poses indicate that the carboxylic group of the three candidates generated H-bonds with residues in the protein chain from Ser441 to Lys444 and formed H-bond, π-cation interactions, or hydrophobic contacts with Phe297 and key active residue, His363. During MD, stable π-cation interactions with residues Phe273 or Arg274 were formed by (S)-tryptophan-betaxanthin and RosA. All candidates were anchored to His363 by stable π- or H-bonds. Hence, we propose (S)-tryptophan-betaxanthin, 5-O-feruloylquinic acid, and RosA as potential lead compounds that can be further tested in drug development process for diseases associated with aging An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:28.
Collapse
Affiliation(s)
- Kuan-Chung Chen
- a Graduate Institute of Pharmaceutical Chemistry, China Medical University , Taichung , 40402 , Taiwan
| | | | | | | | | | | |
Collapse
|
15
|
Chen KC, Chang SS, Huang HJ, Lin TL, Wu YJ, Chen CYC. Three-in-one agonists for PPAR-α, PPAR-γ, and PPAR-δ from traditional Chinese medicine. J Biomol Struct Dyn 2012; 30:662-83. [DOI: 10.1080/07391102.2012.689699] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
16
|
Kalyaanamoorthy S, Chen YPP. Exploring inhibitor release pathways in histone deacetylases using random acceleration molecular dynamics simulations. J Chem Inf Model 2012; 52:589-603. [PMID: 22263580 DOI: 10.1021/ci200584f] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Molecular channel exploration perseveres to be the prominent solution for eliciting structure and accessibility of active site and other internal spaces of macromolecules. The volume and silhouette characterization of these channels provides answers for the issues of substrate access and ligand swapping between the obscured active site and the exterior of the protein. Histone deacetylases (HDACs) are metal-dependent enzymes that are involved in the cell growth, cell cycle regulation, and progression, and their deregulations have been linked with different types of cancers. Hence HDACs, especially the class I family, are widely recognized as the important cancer targets, and the characterizations of their structures and functions have been of special interest in cancer drug discovery. The class I HDACs are known to possess two different protein channels, an 11 Å and a 14 Å (named channels A and B1, respectively), of which the former is a ligand or substrate occupying tunnel that leads to the buried active site zinc ion and the latter is speculated to be involved in product release. In this work, we have carried out random acceleration molecular dynamics (RAMD) simulations coupled with the classical molecular dynamics to explore the release of the ligand, N-(2-aminophenyl) benzamide (LLX) from the active sites of the recently solved X-ray crystal structure of HDAC2 and the computationally modeled HDAC1 proteins. The RAMD simulations identified significant structural and dynamic features of the HDAC channels, especially the key 'gate-keeping' amino acid residues that control these channels and the ligand release events. Further, this study identified a novel and unique channel B2, a subchannel from channel B1, in the HDAC1 protein structure. The roles of water molecules in the LLX release from the HDAC1 and HDAC2 enzymes are also discussed. Such structural and dynamic properties of the HDAC protein channels that govern the ligand escape reactions will provide further mechanistic insights into the HDAC enzymes, which, in the long run, have a potential to bring new ideas for developing more promising HDAC inhibitors as well as extend our atomic level understandings on their mechanisms of action.
Collapse
Affiliation(s)
- Subha Kalyaanamoorthy
- Department of Computer Science and Computer Engineering, Faculty of Science, Technology and Engineering, La Trobe University, Melbourne, Australia
| | | |
Collapse
|
17
|
Schmidtke P, Bidon-Chanal A, Luque FJ, Barril X. MDpocket: open-source cavity detection and characterization on molecular dynamics trajectories. ACTA ACUST UNITED AC 2011; 27:3276-85. [PMID: 21967761 DOI: 10.1093/bioinformatics/btr550] [Citation(s) in RCA: 212] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
MOTIVATION A variety of pocket detection algorithms are now freely or commercially available to the scientific community for the analysis of static protein structures. However, since proteins are dynamic entities, enhancing the capabilities of these programs for the straightforward detection and characterization of cavities taking into account protein conformational ensembles should be valuable for capturing the plasticity of pockets, and therefore allow gaining insight into structure-function relationships. RESULTS This article describes a new method, called MDpocket, providing a fast, free and open-source tool for tracking small molecule binding sites and gas migration pathways on molecular dynamics (MDs) trajectories or other conformational ensembles. MDpocket is based on the fpocket cavity detection algorithm and a valuable contribution to existing analysis tools. The capabilities of MDpocket are illustrated for three relevant cases: (i) the detection of transient subpockets using an ensemble of crystal structures of HSP90; (ii) the detection of known xenon binding sites and migration pathways in myoglobin; and (iii) the identification of suitable pockets for molecular docking in P38 Map kinase. AVAILABILITY MDpocket is free and open-source software and can be downloaded at http://fpocket.sourceforge.net. CONTACT pschmidtke@ub.edu SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Peter Schmidtke
- Departament de Fisicoquímica and Institut de Biomedicina, Universitat de Barcelona, 08028 Barcelona, Spain.
| | | | | | | |
Collapse
|