1
|
Dhayalan A, Thillainathan N, Velramar B, Athiyappagounder P, Sundaramoorthy D, Pachiappan P. Pectinase from a Fish Gut Bacterium, Aeromonas guangheii (SS6): Production, Cloning and Characterization. Protein J 2022; 41:572-590. [PMID: 36208356 DOI: 10.1007/s10930-022-10077-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/20/2022] [Indexed: 11/28/2022]
Abstract
During the present research, 11 gut bacteria were isolated from the freshwater fish, Systomus sarana (General name: olive barb) and upon screening, the strains produced extracellular pectinase enzyme. Among them, the SS6 strain was found to produce a high quantity of 208.731 U/ml pectinase and through molecular characterization the SS6 strain was identified as Aeromonas guangheii. During the culture of SS6 strain, a set of parameters were optimized through the response surface methodology with a Box-Behnken design, for the production of the enzyme. The optimal conditions were found to be 2.11% of maltose, 2.20% of yeast extract, 6.5 of pH, and a temperature of 27.3 °C at 32-h incubation. Under the above conditions, the activity of pectinase production was enhanced to 371 U/ml. The purified pectinase's molecular weight was determined to be ~ 50 kDa (by 10% 2-D PAGE). Totally, nine peptides were identified from the purified pectinase enzyme through the MALDI-TOF-MS analysis and MASCOT tool was used to get the mass spectrum of the peak at 2211 of peptide that indicated the reference pectinase protein. The referenced gene primer (pectate lyases) was PCR amplified and its nucleotide sequence was analyzed. The exo-pelA gene was cloned in pREST vector, which was found to be over expressed in Escherichia coli BL21. The ORF encoded for a mature protein comprising of 425 amino acids (1236 nucleotides) with a predicted molecular weight of ~ 48.7 kDa. The present findings underline the potential of the fish-gut microbes as a source of biotechnologically important enzymes.
Collapse
Affiliation(s)
- Arul Dhayalan
- Department of Biotechnology, School of Biosciences, Periyar University, Salem, Tamil Nadu, 636011, India.,SRS of ICAR- National Dairy Research Institute, Adugodi, Bengaluru, Karnataka, 560030, India
| | - Natarajan Thillainathan
- Department of Biotechnology, School of Biosciences, Periyar University, Salem, Tamil Nadu, 636011, India.,Department of Biomedical Engineering, Central University of Rajasthan, Ajmer, Rajasthan, 305817, India
| | - Balasubramanian Velramar
- Department of Biotechnology, School of Biosciences, Periyar University, Salem, Tamil Nadu, 636011, India.,Amity Institute of Biotechnology, Amity University, Raipur, Chhattisgarh, 493225, India
| | - Palanisammi Athiyappagounder
- Veterinary College & Research Institute, Tamil Nadu Veterinary & Animal Science University, Tirunelveli, Tamil Nadu, 627358, India
| | - Dhanasundaram Sundaramoorthy
- Department of Marine Science, School of Marine Sciences, Bharathidasan University, Tiruchirappalli, Tamil Nadu, 620024, India
| | - Perumal Pachiappan
- Department of Biotechnology, School of Biosciences, Periyar University, Salem, Tamil Nadu, 636011, India. .,Department of Marine Science, School of Marine Sciences, Bharathidasan University, Tiruchirappalli, Tamil Nadu, 620024, India.
| |
Collapse
|
2
|
Bavo F, Pallavicini M, Appiani R, Bolchi C. Determinants for α4β2 vs. α3β4 Subtype Selectivity of Pyrrolidine-Based nAChRs Ligands: A Computational Perspective with Focus on Recent cryo-EM Receptor Structures. Molecules 2021; 26:molecules26123603. [PMID: 34204637 PMCID: PMC8231201 DOI: 10.3390/molecules26123603] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Revised: 06/08/2021] [Accepted: 06/10/2021] [Indexed: 12/26/2022] Open
Abstract
The selectivity of α4β2 nAChR agonists over the α3β4 nicotinic receptor subtype, predominant in ganglia, primarily conditions their therapeutic range and it is still a complex and challenging issue for medicinal chemists and pharmacologists. Here, we investigate the determinants for such subtype selectivity in a series of more than forty α4β2 ligands we have previously reported, docking them into the structures of the two human subtypes, recently determined by cryo-electron microscopy. They are all pyrrolidine based analogues of the well-known α4β2 agonist N-methylprolinol pyridyl ether A-84543 and differ in the flexibility and pattern substitution of their aromatic portion. Indeed, the direct or water mediated interaction with hydrophilic residues of the relatively narrower β2 minus side through the elements decorating the aromatic ring and the stabilization of the latter by facing to the not conserved β2-Phe119 result as key distinctive features for the α4β2 affinity. Consistently, these compounds show, despite the structural similarity, very different α4β2 vs. α3β4 selectivities, from modest to very high, which relate to rigidity/extensibility degree of the portion containing the aromatic ring and to substitutions at the latter. Furthermore, the structural rationalization of the rat vs. human differences of α4β2 vs. α3β4 selectivity ratios is here proposed.
Collapse
Affiliation(s)
- Francesco Bavo
- Dipartimento di Scienze Farmaceutiche, Universita’degli Studi di Milano, I-20133 Milano, Italy; (F.B.); (M.P.); (R.A.)
- Department of Drug Design and Pharmacology, University of Copenhagen, DK-2100 Copenhagen, Denmark
| | - Marco Pallavicini
- Dipartimento di Scienze Farmaceutiche, Universita’degli Studi di Milano, I-20133 Milano, Italy; (F.B.); (M.P.); (R.A.)
| | - Rebecca Appiani
- Dipartimento di Scienze Farmaceutiche, Universita’degli Studi di Milano, I-20133 Milano, Italy; (F.B.); (M.P.); (R.A.)
| | - Cristiano Bolchi
- Dipartimento di Scienze Farmaceutiche, Universita’degli Studi di Milano, I-20133 Milano, Italy; (F.B.); (M.P.); (R.A.)
- Correspondence:
| |
Collapse
|
3
|
Chen Z, Yao X, Dong F, Duan H, Shao X, Chen X, Yang T, Wang G, Zheng Y. Ecological toxicity reduction of dinotefuran to honeybee: New perspective from an enantiomeric level. ENVIRONMENT INTERNATIONAL 2019; 130:104854. [PMID: 31200156 DOI: 10.1016/j.envint.2019.05.048] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Revised: 05/15/2019] [Accepted: 05/17/2019] [Indexed: 06/09/2023]
Abstract
In last decade, there has been a concerted effort to reduce the potential threats of honeybees' population due to exposure to neonicotinoid pesticides. A new perspective was put forward to reduce the potential ecological toxicity of neonicotinoid dinotefuran to honeybee in terms of an enantiomeric level in the study. Toxicity of dinotefuran was enantioselective, and S-dinotefuran was 41.1- to 128.4-fold more toxic than R-dinotefuran to honeybee Apis mellifera (Apis mellifera Linnaeus), whereas R-dinotefuran exhibited comparative insecticidal activities (1.7-2.4 times) to typical sucking pests Aphis gossypii and Apolygus lucorum compared to racemic mixtures. Our data suggested that use of R-dinotefuran could have a good efficacy in controlling target pests while minimizing hazard to honeybees. The mechanism for chiral specific toxicity to honeybee was further characterized by electrophysiological studies and molecular docking. S-dinotefuran appears to be more toxic by binding to α8 subunit of nAChR of Apis mellifera. The α8 also have a more stable, functional binding cavity to S-dinotefuran with a higher binding score of 7.15, primarily due to an extensive hydrogen bond network. Therefore, new chiral products with a high proportion of or an enantiomeric pure R-dinotefuran are recommended to achieve effective pests control reducing hazard to honeybee populations.
Collapse
Affiliation(s)
- Zenglong Chen
- State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, PR China; State Key Laboratory of Integrated Management of Pest Insects and Rodents, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, PR China
| | - Xiangmei Yao
- State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, PR China
| | - Fengshou Dong
- State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, PR China.
| | - Hongxia Duan
- Department of Applied Chemistry, College of Science, China Agricultural University, Beijing 100193, PR China
| | - Xusheng Shao
- School of Pharmacy, East China University of Science and Technology, Shanghai 200237, PR China
| | - Xiu Chen
- State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, PR China
| | - Ting Yang
- State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, PR China
| | - Guirong Wang
- State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, PR China.
| | - Yongquan Zheng
- State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, PR China.
| |
Collapse
|
4
|
Molecular docking and in vitro studies of soap nut trypsin inhibitor (SNTI) against phospholipase A2 isoforms in therapeutic intervention of inflammatory diseases. Int J Biol Macromol 2018; 114:556-564. [DOI: 10.1016/j.ijbiomac.2018.03.139] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2017] [Revised: 02/26/2018] [Accepted: 03/22/2018] [Indexed: 11/20/2022]
|
5
|
Inhibition of Mycobacterium-RmlA by Molecular Modeling, Dynamics Simulation, and Docking. Adv Bioinformatics 2016; 2016:9841250. [PMID: 26981117 PMCID: PMC4769735 DOI: 10.1155/2016/9841250] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2015] [Revised: 12/12/2015] [Accepted: 12/24/2015] [Indexed: 11/17/2022] Open
Abstract
The increasing resistance to anti-tb drugs has enforced strategies for finding new drug targets against Mycobacterium tuberculosis (Mtb). In recent years enzymes associated with the rhamnose pathway in Mtb have attracted attention as drug targets. The present work is on α-D-glucose-1-phosphate thymidylyltransferase (RmlA), the first enzyme involved in the biosynthesis of L-rhamnose, of Mtb cell wall. This study aims to derive a 3D structure of RmlA by using a comparative modeling approach. Structural refinement and energy minimization of the built model have been done with molecular dynamics. The reliability assessment of the built model was carried out with various protein checking tools such as Procheck, Whatif, ProsA, Errat, and Verify 3D. The obtained model investigates the relation between the structure and function. Molecular docking interactions of Mtb-RmlA with modified EMB (ethambutol) ligands and natural substrate have revealed specific key residues Arg13, Lys23, Asn109, and Thr223 which play an important role in ligand binding and selection. Compared to all EMB ligands, EMB-1 has shown better interaction with Mtb-RmlA model. The information thus discussed above will be useful for the rational design of safe and effective inhibitors specific to RmlA enzyme pertaining to the treatment of tuberculosis.
Collapse
|
6
|
Soleimani M, Mahnam K, Mirmohammad-Sadeghi H, Sadeghi-Aliabadi H, Jahanian-Najafabadi A. Theoretical design of a new chimeric protein for the treatment of breast cancer. Res Pharm Sci 2016; 11:187-99. [PMID: 27499788 PMCID: PMC4962299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022] Open
Abstract
p28 and NRC peptides are two anticancer peptides with various mechanisms have shown to be effective against breast cancer. Therefore, it seems that construction of a chimeric protein containing the two peptides might cause synergistic cytotoxic effects. However, since the two peptides bear opposite charges, production of a chimeric protein in which the two moieties do not intervene each other is difficult. In this study, our goal was to find a suitable peptide linker for the new chimeric protein in a manner that none of the peptides intervene the other's function. We selected some linkers with different characteristics and lengths and created a small library of the chimeric proteins harboring these linkers. Homology modeling and molecular dynamic simulation revealed that (PA)5P and (EAAAK)3 linkers can separate the p28 and NRC peptides effectively. Thus, the chimeric protein linked with (PA)5P or (EAAAK)3 linkers might show synergistic and stronger anticancer effects than the separate peptide moieties because they could exert their cytotoxic effects freely which is not influenced by the other part.
Collapse
Affiliation(s)
- Meysam Soleimani
- Department of Pharmaceutical Biotechnology, School of Pharmacy and Pharmaceutical Sciences, Isfahan University of Medical Sciences, Isfahan, I.R. Iran
| | - Karim Mahnam
- Biology Department, Faculty of Sciences, Shahrekord University, Shahrekord, I.R. Iran
| | - Hamid Mirmohammad-Sadeghi
- Department of Pharmaceutical Biotechnology, School of Pharmacy and Pharmaceutical Sciences, Isfahan University of Medical Sciences, Isfahan, I.R. Iran,Isfahan Pharmaceutical Sciences Research Center, School of Pharmacy and Pharmaceutical Sciences, Isfahan University of Medical Sciences, Isfahan, I.R. Iran
| | - Hojjat Sadeghi-Aliabadi
- Isfahan Pharmaceutical Sciences Research Center, School of Pharmacy and Pharmaceutical Sciences, Isfahan University of Medical Sciences, Isfahan, I.R. Iran,Department of Medicinal Chemistry, School of Pharmacy and Pharmaceutical Sciences, Isfahan University of Medical Sciences, Isfahan, I.R. Iran
| | - Ali Jahanian-Najafabadi
- Department of Pharmaceutical Biotechnology, School of Pharmacy and Pharmaceutical Sciences, Isfahan University of Medical Sciences, Isfahan, I.R. Iran,Isfahan Pharmaceutical Sciences Research Center, School of Pharmacy and Pharmaceutical Sciences, Isfahan University of Medical Sciences, Isfahan, I.R. Iran,Corresponding author: A. Jahanian-Najafabadi Tel: 0098 31 37927056, Fax: 0098 31 36680011
| |
Collapse
|
7
|
Salari R, Murlidaran S, Brannigan G. Pentameric Ligand-gated Ion Channels : Insights from Computation. MOLECULAR SIMULATION 2014; 40:821-829. [PMID: 25931676 PMCID: PMC4412168 DOI: 10.1080/08927022.2014.896462] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
Pentameric ligand-gated ion channels (pLGICs) conduct upon the binding of an agonist and are fundamental to neurotransmission. New insights into the complex mechanisms underlying pLGIC gating, ion selectivity, and modulation have recently been gained via a series of crystal structures in prokaryotes and C .elegans, as well as computational studies relying on these structures. Here we review contributions from a variety of computational approaches, including normal mode analysis, automated docking, and fully atomistic molecular dynamics simulation. Examples from our own research, particularly concerning interactions with general anesthetics and lipids, are used to illustrate predictive results complementary to crystallographic studies.
Collapse
Affiliation(s)
- Reza Salari
- Center for Computational and Integrative Biology, Rutgers University, Camden, NJ
- Department of Physics, Rutgers University, Camden, NJ
| | - Sruthi Murlidaran
- Center for Computational and Integrative Biology, Rutgers University, Camden, NJ
| | - Grace Brannigan
- Center for Computational and Integrative Biology, Rutgers University, Camden, NJ
- Department of Physics, Rutgers University, Camden, NJ
| |
Collapse
|
8
|
Structural analysis and ion translocation mechanisms of the muscle-type acetylcholine receptor channel. J Appl Biomater Funct Mater 2013; 11:e53-60. [PMID: 23728540 DOI: 10.5301/jabfm.5000148] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/31/2012] [Indexed: 12/13/2022] Open
Abstract
PURPOSE The aim of this work is to analyze the conformational changes in the acetylcholine receptor caused by channel opening and to investigate the electrostatic profile during ion translocation through the channel. METHODS A computational model of the human muscle-type acetylcholine receptor (AChR) was built and used to analyze channel structure and its interactions with different ions. Using the Torpedo AChR crystal structure as a homologous template, the 3D structure of the human muscle-type AChR was reconstructed. RESULTS This first model is optimized and an open structure of the channel is generated using Normal Mode Analysis in order to assess morphologic and energetic differences between open and closed structures. In addition, the issue of ion translocation is investigated in further detail. Results elucidate different aspects of the channel: channel gate structure, channel interactions with translocating ions, differences between muscle-type AChR and previous neuronal-type AChR models. CONCLUSIONS The model constructed here is ideal for further computational studies on muscle-type AChR and its pathologic mutations.
Collapse
|
9
|
Ashby JA, McGonigle IV, Price KL, Cohen N, Comitani F, Dougherty DA, Molteni C, Lummis SCR. GABA binding to an insect GABA receptor: a molecular dynamics and mutagenesis study. Biophys J 2012. [PMID: 23200041 DOI: 10.1016/j.bpj.2012.10.016] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
RDL receptors are GABA-activated inhibitory Cys-loop receptors found throughout the insect CNS. They are a key target for insecticides. Here, we characterize the GABA binding site in RDL receptors using computational and electrophysiological techniques. A homology model of the extracellular domain of RDL was generated and GABA docked into the binding site. Molecular dynamics simulations predicted critical GABA binding interactions with aromatic residues F206, Y254, and Y109 and hydrophilic residues E204, S176, R111, R166, S176, and T251. These residues were mutated, expressed in Xenopus oocytes, and their functions assessed using electrophysiology. The data support the binding mechanism provided by the simulations, which predict that GABA forms many interactions with binding site residues, the most significant of which are cation-π interactions with F206 and Y254, H-bonds with E204, S205, R111, S176, T251, and ionic interactions with R111 and E204. These findings clarify the roles of a range of residues in binding GABA in the RDL receptor, and also show that molecular dynamics simulations are a useful tool to identify specific interactions in Cys-loop receptors.
Collapse
|
10
|
Li Q, Fang Y, Li X, Zhang H, Liu M, Yang H, Kang Z, Li Y, Wang Y. Mechanism of the plant cytochrome P450 for herbicide resistance: a modelling study. J Enzyme Inhib Med Chem 2012; 28:1182-91. [PMID: 23057845 DOI: 10.3109/14756366.2012.719505] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Plant cytochrome P450 is a key enzyme responsible for the herbicide resistance but the molecular basis of the mechanism is unclear. To understand this, four typical plant P450s and a widely resistant herbicide chlortoluron were analysed by carrying out homology modelling, molecular docking, molecular dynamics simulations and binding free energy analysis. Our results demonstrate that: (i) the putative hydrophobic residues located in the F-helix and polar residues in I-helix are critical in the herbicide resistance; (ii) the binding mode analysis and binding free energy calculation indicate that the distance between catalytic site of chlortoluron and heme of P450, as well as the binding affinity are key elements affecting the resistance for plants. In conclusion, this work provides a new insight into the interactions of plant P450s with herbicide from a molecular level, offering valuable information for the future design of novel effective herbicides which also escape from the P450 metabolism.
Collapse
Affiliation(s)
- Qinfan Li
- College of Veterinary Medicine, Northwest A&F University , Yangling, Shaanxi , China
| | | | | | | | | | | | | | | | | |
Collapse
|
11
|
Yu R, Kaas Q, Craik DJ. Delineation of the unbinding pathway of α-conotoxin ImI from the α7 nicotinic acetylcholine receptor. J Phys Chem B 2012; 116:6097-105. [PMID: 22571488 DOI: 10.1021/jp301352d] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
α-Conotoxins potently and specifically inhibit isoforms of nicotinic acetylcholine receptors (nAChRs) and are used as molecular probes and as drugs or drug leads. Interactions occurring during binding and unbinding events are linked to binding kinetics, and knowledge of these interactions could help in the development of α-conotoxins as drugs. Here, the unbinding process for the prototypical α-conotoxin ImI/α7-nAChR system was investigated theoretically, and three exit routes were identified using random accelerated molecular dynamics simulations. The route involving the smallest conformation perturbation was further divided into three subpathways, which were studied using steered molecular dynamics simulations. Of the three subpathways, two had better experimental support and lower potential of mean force, indicating that they might be sampled more frequently. Additionally, these subpathways were supported by previous experimental studies. Several pairwise interactions, including a cation-π interaction and charge and hydrogen bond interactions, were identified as potentially playing important roles in the unbinding event.
Collapse
Affiliation(s)
- Rilei Yu
- Division of Chemistry and Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Australia
| | | | | |
Collapse
|
12
|
NMR structures of the transmembrane domains of the α4β2 nAChR. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2012; 1818:1261-8. [PMID: 22361591 DOI: 10.1016/j.bbamem.2012.02.008] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/20/2011] [Revised: 01/30/2012] [Accepted: 02/07/2012] [Indexed: 11/21/2022]
Abstract
The α4β2 nicotinic acetylcholine receptor (nAChR) is the predominant heteromeric subtype of nAChRs in the brain, which has been implicated in numerous neurological conditions. The structural information specifically for the α4β2 and other neuronal nAChRs is presently limited. In this study, we determined structures of the transmembrane (TM) domains of the α4 and β2 subunits in lauryldimethylamine-oxide (LDAO) micelles using solution NMR spectroscopy. NMR experiments and size exclusion chromatography-multi-angle light scattering (SEC-MALS) analysis demonstrated that the TM domains of α4 and β2 interacted with each other and spontaneously formed pentameric assemblies in the LDAO micelles. The Na(+) flux assay revealed that α4β2 formed Na(+) permeable channels in lipid vesicles. Efflux of Na(+) through the α4β2 channels reduced intra-vesicle Sodium Green™ fluorescence in a time-dependent manner that was not observed in vesicles without incorporating α4β2. The study provides structural insight into the TM domains of the α4β2 nAChR. It offers a valuable structural framework for rationalizing extensive biochemical data collected previously on the α4β2 nAChR and for designing new therapeutic modulators.
Collapse
|
13
|
Mohan S, Sheena A, Poulose N, Anilkumar G. Molecular dynamics simulation studies of GLUT4: substrate-free and substrate-induced dynamics and ATP-mediated glucose transport inhibition. PLoS One 2010; 5:e14217. [PMID: 21151967 PMCID: PMC2997047 DOI: 10.1371/journal.pone.0014217] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2010] [Accepted: 11/15/2010] [Indexed: 02/03/2023] Open
Abstract
Background Glucose transporter 4 (GLUT4) is an insulin facilitated glucose transporter that plays an important role in maintaining blood glucose homeostasis. GLUT4 is sequestered into intracellular vesicles in unstimulated cells and translocated to the plasma membrane by various stimuli. Understanding the structural details of GLUT4 will provide insights into the mechanism of glucose transport and its regulation. To date, a crystal structure for GLUT4 is not available. However, earlier work from our laboratory proposed a well validated homology model for GLUT4 based on the experimental data available on GLUT1 and the crystal structure data obtained from the glycerol 3-phosphate transporter. Methodology/Principal Findings In the present study, the dynamic behavior of GLUT4 in a membrane environment was analyzed using three forms of GLUT4 (apo, substrate and ATP-substrate bound states). Apo form simulation analysis revealed an extracellular open conformation of GLUT4 in the membrane favoring easy exofacial binding of substrate. Simulation studies with the substrate bound form proposed a stable state of GLUT4 with glucose, which can be a substrate-occluded state of the transporter. Principal component analysis suggested a clockwise movement for the domains in the apo form, whereas ATP substrate-bound form induced an anti-clockwise rotation. Simulation studies suggested distinct conformational changes for the GLUT4 domains in the ATP substrate-bound form and favor a constricted behavior for the transport channel. Various inter-domain hydrogen bonds and switching of a salt-bridge network from E345-R350-E409 to E345-R169-E409 contributed to this ATP-mediated channel constriction favoring substrate occlusion and prevention of its release into cytoplasm. These data are consistent with the biochemical studies, suggesting an inhibitory role for ATP in GLUT-mediated glucose transport. Conclusions/Significance In the absence of a crystal structure for any glucose transporter, this study provides mechanistic details of the conformational changes in GLUT4 induced by substrate and its regulator.
Collapse
Affiliation(s)
- Suma Mohan
- Amrita School of Biotechnology, Amrita Vishwa Vidyapeetham, Kollam, Kerala, India
| | - Aswathy Sheena
- Amrita School of Biotechnology, Amrita Vishwa Vidyapeetham, Kollam, Kerala, India
| | - Ninu Poulose
- Amrita School of Biotechnology, Amrita Vishwa Vidyapeetham, Kollam, Kerala, India
| | | |
Collapse
|
14
|
Willenbring D, Xu Y, Tang P. The role of structured water in mediating general anesthetic action on alpha4beta2 nAChR. Phys Chem Chem Phys 2010; 12:10263-9. [PMID: 20661501 PMCID: PMC3265171 DOI: 10.1039/c003573d] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Water is an essential component for many biological processes. Pauling proposed that water might play a critical role in general anesthesia by forming water clathrates around anesthetic molecules. To examine potential involvement of water in general anesthesia, we analyzed water within alpha4beta2 nAChR, a putative protein target hypersensitive to volatile anesthetics. Experimental structure-derived closed- and open-channel nAChR systems in a fully hydrated lipid bilayer were examined using all-atom molecular dynamics simulations. At the majority of binding sites in alpha4beta2 nAChR, halothane replaced the slow-exchanging water molecules and caused a regional water population decrease. Only two binding sites had an increased quantity of water in the presence of halothane, where water arrangements resemble clathrate-like structures. The small number of such clathrate-like water clusters suggests that the formation of water clathrates is unlikely to be a primary cause for anesthesia. Despite the decrease in water population at most of the halothane binding sites, the number of sites that can be occupied transiently by water is actually increased in the presence of halothane. Many of these water sites were located between two subunits or in regions containing agonist binding sites or critical structural elements for transducing agonist binding to channel gating. Changes in water sites in the presence of halothane affected water-mediated protein-protein interactions and the protein dynamics, which can have direct impact on protein function. Collectively, water contributes to the action of anesthetics in proteins by mediating interactions between protein subunits and altering protein dynamics, instead of forming water clathrates around anesthetics.
Collapse
Affiliation(s)
- Dan Willenbring
- Department of Anesthesiology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261
| | - Yan Xu
- Department of Anesthesiology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261
- Department of Pharmacology and Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261
| | - Pei Tang
- Department of Anesthesiology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261
- Department of Pharmacology and Chemical Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261
- Department of Computational Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261
| |
Collapse
|
15
|
NMR structure of the transmembrane domain of the n-acetylcholine receptor beta2 subunit. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2010; 1798:1608-14. [PMID: 20441771 DOI: 10.1016/j.bbamem.2010.04.014] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/15/2010] [Revised: 04/21/2010] [Accepted: 04/26/2010] [Indexed: 11/20/2022]
Abstract
Nicotinic acetylcholine receptors (nAChRs) are involved in fast synaptic transmission in the central and peripheral nervous system. Among the many different types of subunits in nAChRs, the beta2 subunit often combines with the alpha4 subunit to form alpha4beta2 pentameric channels, the most abundant subtype of nAChRs in the brain. Besides computational predictions, there is limited experimental data available on the structure of the beta2 subunit. Using high-resolution NMR spectroscopy, we solved the structure of the entire transmembrane domain (TM1234) of the beta2 subunit. We found that TM1234 formed a four-helix bundle in the absence of the extracellular and intracellular domains. The structure exhibited many similarities to those previously determined for the Torpedo nAChR and the bacterial ion channel GLIC. We also assessed the influence of the fourth transmembrane helix (TM4) on the rest of the domain. Although secondary structures and tertiary arrangements were similar, the addition of TM4 caused dramatic changes in TM3 dynamics and subtle changes in TM1 and TM2. Taken together, this study suggests that the structures of the transmembrane domains of these proteins are largely shaped by determinants inherent in their sequence, but their dynamics may be sensitive to modulation by tertiary and quaternary contacts.
Collapse
|
16
|
Exploring the molecular basis for selective binding of Mycobacterium tuberculosis Asp kinase toward its natural substrates and feedback inhibitors: a docking and molecular dynamics study. J Mol Model 2010; 16:1357-67. [PMID: 20140471 DOI: 10.1007/s00894-010-0653-4] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2009] [Accepted: 01/05/2010] [Indexed: 10/19/2022]
Abstract
Tuberculosis (TB) is still a major public health problem, compounded by the human immunodeficiency virus (HIV)-TB co-infection and recent emergence of multidrug-resistant (MDR) and extensively drug resistant (XDR)-TB. In this context, aspartokinase of mycobacterium tuberculosis has drawn attention for designing novel anti-TB drugs. Asp kinase is an enzyme responsible for the synthesis of 4-phospho-L-aspartate from L-aspartate and involved in the branched biosynthetic pathway leading to the synthesis of amino acids lysine, threonine, methionine and isoleucine. An intermediate of lysine biosynthetic branch, mesodiaminopimelate is also a component of the peptidoglycan which is a component of bacterial cell wall. To interfere with the production of all these amino acids and cell wall, it is possible to inhibit Asp kinase activity. This can be achieved using Asp kinase inhibitors. In order to design novel Asp kinase inhibitors as effective anti-TB drugs, it is necessary to have an understanding of the binding sites of Asp kinase. As no crystal structure of the enzyme has yet been published, we built a homology model of Asp kinase using the crystallized Asp kinase from M. Jannaschii, as template structures (2HMF and 3C1M). After the molecular dynamics refinement, the optimized homology model was assessed as a reliable structure by PROCHECK, ERRAT, WHAT-IF, PROSA2003 and VERIFY-3D. The results of molecular docking studies with natural substrates, products and feedback inhibitors are in agreement with the published data and showed that ACT domain plays an important role in binding to ligands. Based on the docking conformations, pharmacophore model can be developed by probing the common features of ligands. By analyzing the results, ACT domain architecture, certain key residues that are responsible for binding to feedback inhibitors and natural substrates were identified. This would be very helpful in understanding the blockade mechanism of Asp kinase and providing insights into rational design of novel Asp kinase inhibitors for M.tuberculosis.
Collapse
|
17
|
Liu LT, Willenbring D, Xu Y, Tang P. General anesthetic binding to neuronal alpha4beta2 nicotinic acetylcholine receptor and its effects on global dynamics. J Phys Chem B 2009; 113:12581-9. [PMID: 19697903 DOI: 10.1021/jp9039513] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is a target for general anesthetics. Currently available experimental structural information is inadequate to understand where anesthetics bind and how they modulate the receptor motions essential to function. Using our published open-channel structure model of alpha4beta2 nAChR, we identified and evaluated six amphiphilic interaction sites for the volatile anesthetic halothane via flexible ligand docking and subsequent 20-ns molecular dynamics simulations. Halothane binding energies ranged from -6.8 to -2.4 kcal/mol. The primary binding sites were located at the interface of extracellular and transmembrane domains, where halothane perturbed conformations of, and widened the gap among, the Cys loop, the beta1-beta2 loop, and the TM2-TM3 linker. The halothane with the highest binding affinity at the interface between the alpha4 and beta2 subunits altered interactions between the protein and nearby lipids by competing for hydrogen bonds. Gaussian network model analyses of the alpha4beta2 nAChR structures at the end of 20-ns simulations in the absence or presence of halothane revealed profound changes in protein residue mobility. The concerted motions critical to protein function were also perturbed considerably. Halothane's effect on protein dynamics was not confined to the residues adjacent to the binding sites, indicating an action on a more global scale.
Collapse
Affiliation(s)
- Lu Tian Liu
- Department of Anesthesiology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | | | | | | |
Collapse
|
18
|
Cheng X, Ivanov I, Wang H, Sine SM, McCammon JA. Molecular-dynamics simulations of ELIC-a prokaryotic homologue of the nicotinic acetylcholine receptor. Biophys J 2009; 96:4502-13. [PMID: 19486673 DOI: 10.1016/j.bpj.2009.03.018] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2008] [Revised: 03/05/2009] [Accepted: 03/11/2009] [Indexed: 11/19/2022] Open
Abstract
The ligand-gated ion channel from Erwinia chrysanthemi (ELIC) is a prokaryotic homolog of the eukaryotic nicotinic acetylcholine receptor (nAChR) that responds to the binding of neurotransmitter acetylcholine and mediates fast signal transmission. ELIC is similar to the nAChR in its primary sequence and overall subunit organization, but despite their structural similarity, it is not clear whether these two ligand-gated ion channels operate in a similar manner. Further, it is not known to what extent mechanistic insights gleaned from the ELIC structure translate to eukaryotic counterparts such as the nAChR. Here we use molecular-dynamics simulations to probe the conformational dynamics and hydration of the transmembrane pore of ELIC. The results are compared with those from our previous simulation of the human alpha7 nAChR. Overall, ELIC displays increased stability compared to the nAChR, whereas the two proteins exhibit remarkable similarity in their global motion and flexibility patterns. The majority of the increased stability of ELIC does not stem from the deficiency of the models used in the simulations, and but rather seems to have a structural basis. Slightly altered dynamical correlation features are also observed among several loops within the membrane region. In sharp contrast to the nAChR, ELIC is completely dehydrated from the pore center to the extracellular end throughout the simulation. Finally, the simulation of an ELIC mutant substantiates the important role of F246 on the stability, hydration and possibly function of the ELIC channel.
Collapse
Affiliation(s)
- Xiaolin Cheng
- Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA.
| | | | | | | | | |
Collapse
|
19
|
Pauli I, Macedo Timmers LFS, Andrade Caceres R, Augusto Basso L, Santiago Santos D, Filgueira de Azevedo Jr. W. Molecular modeling and dynamics studies of purine nucleoside phosphorylase from Bacteroides fragilis. J Mol Model 2009; 15:913-22. [DOI: 10.1007/s00894-008-0445-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2008] [Accepted: 12/08/2008] [Indexed: 10/21/2022]
|
20
|
Parthiban M, Rajasekaran MB, Ramakumar S, Shanmughavel P. Molecular Modeling of Human Pentameric α7Neuronal Nicotinic Acetylcholine Receptor and Its Interaction with its Agonist and Competitive Antagonist. J Biomol Struct Dyn 2009; 26:535-47. [DOI: 10.1080/07391102.2009.10507269] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
|
21
|
Perez PC, Caceres RA, Canduri F, de Azevedo WF. Molecular modeling and dynamics simulation of human cyclin-dependent kinase 3 complexed with inhibitors. Comput Biol Med 2009; 39:130-40. [PMID: 19152876 DOI: 10.1016/j.compbiomed.2008.11.004] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2007] [Revised: 08/15/2008] [Accepted: 11/20/2008] [Indexed: 11/27/2022]
Abstract
The complex CDK3-cyclin is involved in the control of the progression of G0. While the mechanisms governing early and late G1 progression are well understood, very little is known about the G0-G1 transition. Human CDK3 is closely related to cyclin-dependent kinase 2 (CDK2). Since there is no crystallographic structure of human CDK3, this work describes for the first time a molecular model of human CDK3 complexed with several inhibitors. Comparison of the binary complexes with different inhibitors strongly indicates that those inhibitors should inhibit CDK3 as well as CDK2.
Collapse
Affiliation(s)
- Patrícia Cardoso Perez
- Faculdade de Biociências-Laboratório de Bioquímica Estrutural-PUCRS, Porto Alegre, RS, Brazil
| | | | | | | |
Collapse
|
22
|
Timmers LFSM, Caceres RA, Vivan AL, Gava LM, Dias R, Ducati RG, Basso LA, Santos DS, de Azevedo WF. Structural studies of human purine nucleoside phosphorylase: Towards a new specific empirical scoring function. Arch Biochem Biophys 2008; 479:28-38. [DOI: 10.1016/j.abb.2008.08.015] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2008] [Revised: 08/20/2008] [Accepted: 08/20/2008] [Indexed: 10/21/2022]
|
23
|
Haddadian EJ, Cheng MH, Coalson RD, Xu Y, Tang P. In silico models for the human alpha4beta2 nicotinic acetylcholine receptor. J Phys Chem B 2008; 112:13981-90. [PMID: 18847252 DOI: 10.1021/jp804868s] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is one of the most widely expressed nAChR subtypes in the brain. Its subunits have high sequence identity (54 and 46% for alpha4 and beta2, respectively) with alpha and beta subunits in Torpedo nAChR. Using the known structure of the Torpedo nAChR as a template, the closed-channel structure of the alpha4beta2 nAChR was constructed through homology modeling. Normal-mode analysis was performed on this closed structure and the resulting lowest frequency mode was applied to it for a "twist-to-open" motion, which increased the minimum pore radius from 2.7 to 3.4 A and generated an open-channel model. Nicotine could bind to the predicted agonist binding sites in the open-channel model but not in the closed one. Both models were subsequently equilibrated in a ternary lipid mixture via extensive molecular dynamics (MD) simulations. Over the course of 11 ns MD simulations, the open channel remained open with filled water, but the closed channel showed a much lower water density at its hydrophobic gate comprised of residues alpha4-V259 and alpha4-L263 and their homologous residues in the beta2 subunits. Brownian dynamics simulations of Na+ permeation through the open channel demonstrated a current-voltage relationship that was consistent with experimental data on the conducting state of alpha4beta2 nAChR. Besides establishment of the well-equilibrated closed- and open-channel alpha4beta2 structural models, the MD simulations on these models provided valuable insights into critical factors that potentially modulate channel gating. Rotation and tilting of TM2 helices led to changes in orientations of pore-lining residue side chains. Without concerted movement, the reorientation of one or two hydrophobic side chains could be enough for channel opening. The closed- and open-channel structures exhibited distinct patterns of electrostatic interactions at the interface of extracellular and transmembrane domains that might regulate the signal propagation of agonist binding to channel opening. A potential prominent role of the beta2 subunit in channel gating was also elucidated in the study.
Collapse
Affiliation(s)
- Esmael J Haddadian
- Department of Anesthesiology, University of Pittsburgh School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | | | | | | |
Collapse
|
24
|
Behavioral differences between phosphatidic acid and phosphatidylcholine in the presence of the nicotinic acetylcholine receptor. Biophys J 2008; 95:5637-47. [PMID: 18835908 DOI: 10.1529/biophysj.108.136895] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
It has been found experimentally that negatively charged phosphatidic acid (PA) lipids and cholesterol molecules stabilize the nicotinic acetylcholine receptor (nAChR) in a functional resting state that can participate in an agonist-induced conformational change. In this study, we compare phosphatidylcholine (PC) and PA lipid behavior in the presence of the nAChR to determine why PC lipids do not support a functional nAChR. For lipids that are located within 1.0 nm of the protein, both PC and PA lipids have very similar order-parameter and bilayer-thickness values, which indicate that the annular lipid properties are protein-dependent. The most significant difference between the PC and PA bilayers is the formation of a lipid domain around the protein, which is visible in the PA bilayer but not the PC bilayer. This suggests that the PA domain may help stabilize the nAChR resting state. The PA lipids in the microdomain have a decreased order compared to a homogeneous PA bilayer and the lipids near the protein attempt to increase the free space in their vicinity by residing in multiple lateral planes.
Collapse
|
25
|
Caceres RA, Saraiva Timmers LF, Dias R, Basso LA, Santos DS, de Azevedo WF. Molecular modeling and dynamics simulations of PNP from Streptococcus agalactiae. Bioorg Med Chem 2008; 16:4984-93. [DOI: 10.1016/j.bmc.2008.03.044] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2007] [Revised: 03/12/2008] [Accepted: 03/14/2008] [Indexed: 11/15/2022]
|
26
|
Wells GB. Structural answers and persistent questions about how nicotinic receptors work. FRONTIERS IN BIOSCIENCE : A JOURNAL AND VIRTUAL LIBRARY 2008; 13:5479-510. [PMID: 18508600 PMCID: PMC2430769 DOI: 10.2741/3094] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
The electron diffraction structure of nicotinic acetylcholine receptor (nAChR) from Torpedo marmorata and the X-ray crystallographic structure of acetylcholine binding protein (AChBP) are providing new answers to persistent questions about how nAChRs function as biophysical machines and as participants in cellular and systems physiology. New high-resolution information about nAChR structures might come from advances in crystallography and NMR, from extracellular domain nAChRs as high fidelity models, and from prokaryotic nicotinoid proteins. At the level of biophysics, structures of different nAChRs with different pharmacological profiles and kinetics will help describe how agonists and antagonists bind to orthosteric binding sites, how allosteric modulators affect function by binding outside these sites, how nAChRs control ion flow, and how large cytoplasmic domains affect function. At the level of cellular and systems physiology, structures of nAChRs will help characterize interactions with other cellular components, including lipids and trafficking and signaling proteins, and contribute to understanding the roles of nAChRs in addiction, neurodegeneration, and mental illness. Understanding nAChRs at an atomic level will be important for designing interventions for these pathologies.
Collapse
Affiliation(s)
- Gregg B Wells
- Department of Molecular and Cellular Medicine, College of Medicine, Texas A&M Health Science Center, College Station, TX 77843-1114, USA.
| |
Collapse
|
27
|
Liu X, Xu Y, Li H, Wang X, Jiang H, Barrantes FJ. Mechanics of channel gating of the nicotinic acetylcholine receptor. PLoS Comput Biol 2008; 4:e19. [PMID: 18225945 PMCID: PMC2211534 DOI: 10.1371/journal.pcbi.0040019] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2007] [Accepted: 12/05/2007] [Indexed: 12/31/2022] Open
Abstract
The nicotinic acetylcholine receptor (nAChR) is a key molecule involved in the propagation of signals in the central nervous system and peripheral synapses. Although numerous computational and experimental studies have been performed on this receptor, the structural dynamics of the receptor underlying the gating mechanism is still unclear. To address the mechanical fundamentals of nAChR gating, both conventional molecular dynamics (CMD) and steered rotation molecular dynamics (SRMD) simulations have been conducted on the cryo-electron microscopy (cryo-EM) structure of nAChR embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer and water molecules. A 30-ns CMD simulation revealed a collective motion amongst C-loops, M1, and M2 helices. The inward movement of C-loops accompanying the shrinking of acetylcholine (ACh) binding pockets induced an inward and upward motion of the outer beta-sheet composed of beta9 and beta10 strands, which in turn causes M1 and M2 to undergo anticlockwise motions around the pore axis. Rotational motion of the entire receptor around the pore axis and twisting motions among extracellular (EC), transmembrane (TM), and intracellular MA domains were also detected by the CMD simulation. Moreover, M2 helices undergo a local twisting motion synthesized by their bending vibration and rotation. The hinge of either twisting motion or bending vibration is located at the middle of M2, possibly the gate of the receptor. A complementary twisting-to-open motion throughout the receptor was detected by a normal mode analysis (NMA). To mimic the pulsive action of ACh binding, nonequilibrium MD simulations were performed by using the SRMD method developed in one of our laboratories. The result confirmed all the motions derived from the CMD simulation and NMA. In addition, the SRMD simulation indicated that the channel may undergo an open-close (O <--> C) motion. The present MD simulations explore the structural dynamics of the receptor under its gating process and provide a new insight into the gating mechanism of nAChR at the atomic level.
Collapse
Affiliation(s)
- Xinli Liu
- Department of Engineering Mechanics, State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian, Liaoning, China
| | - Yechun Xu
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Honglin Li
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Xicheng Wang
- Department of Engineering Mechanics, State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian, Liaoning, China
- * To whom correspondence should be addressed. E-mail: (XW); (HJ)
| | - Hualiang Jiang
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- School of Pharmacy, East China University of Science and Technology, Shanghai, China
- * To whom correspondence should be addressed. E-mail: (XW); (HJ)
| | - Francisco J Barrantes
- UNESCO Chair of Biophysics & Molecular Neurobiology and Instituto de Investigaciones Bioquímicas de Bahía Blanca, Bahía Blanca, Argentina
| |
Collapse
|
28
|
Cheng MH, Coalson RD, Cascio M. Molecular dynamics simulations of ethanol binding to the transmembrane domain of the glycine receptor: Implications for the channel potentiation mechanism. Proteins 2007; 71:972-81. [DOI: 10.1002/prot.21784] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
29
|
Szarecka A, Xu Y, Tang P. Dynamics of heteropentameric nicotinic acetylcholine receptor: implications of the gating mechanism. Proteins 2007; 68:948-60. [PMID: 17546671 DOI: 10.1002/prot.21462] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The dynamics characteristics of the currently available structure of Torpedo nicotinic acetylcholine receptor (nAChR), including the extracellular, transmembrane, and intracellular domains (ICDs), were analyzed using the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM). We found that a symmetric quaternary twist motion, reported previously in the literature in a homopentameric receptor (Cheng et al. J Mol Biol 2006;355:310-324; Taly et al. Biophys J 2005;88:3954-3965), occurred also in the heteropentameric Torpedo nAChR. We believe, however, that the symmetric twist alone is not sufficient to explain a large body of experimental data indicating asymmetry and subunit nonequivalence during gating. Here we report our results supporting the hypothesis that a combination of symmetric and asymmetric motions opens the gate. We show that the asymmetric motion involves tilting of the TM2 helices. Furthermore, our study reveals three additional aspects of channel dynamics: (1) loop A serves as an allosteric mediator between the ligand binding loops and those at the domain interface, particularly the linker between TM2 and TM3; (2) the ICD can modulate the pore dynamics and thus should not be neglected in gating studies; and (3) the F loops, which are peculiarly longer and poorly-conserved in non-alpha-subunits, have important dynamical implications.
Collapse
Affiliation(s)
- Agnieszka Szarecka
- Department of Anesthesiology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania 15261, USA
| | | | | |
Collapse
|
30
|
Cheng MH, Cascio M, Coalson RD. Homology modeling and molecular dynamics simulations of the alpha1 glycine receptor reveals different states of the channel. Proteins 2007; 68:581-93. [PMID: 17469203 DOI: 10.1002/prot.21435] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Homology modeling is used to build initial models of the transmembrane domain of the human alpha1 glycine receptor (GlyR) based on the most recently published refined structure of nAChR (PDB ID: 2BG9). Six preliminary GlyR models are constructed using two different approaches. In one approach, five different homopentamers are built by symmetric assembly of alpha1 GlyR subunits using only one of the five unique chains of nAChR as a template. In a second approach, each nAChR subunit serves as a template for an alpha1 GlyR subunit. All six initial GlyR constructs are then embedded into a hydrated POPC lipid bilayer and subjected to molecular dynamics simulation for at least six nanoseconds. Each model is stable throughout the simulation, and the final models fall into three distinct categories. Homopentameric GlyR bundles using a single alpha nAChR subunit as a template appear to be in an open conformation. Under an applied external potential, permeation of Cl(-) ions is observed within several ns in a channel built on an alpha chain. Model channels built on non-alpha chains have a constriction either near the intracellular mouth or more centrally located in the pore domain, both of which may be narrow enough to close the channel and whose locations correspond to putative gates observed in nicotinicoid receptors. The differences between these three general models suggest that channel closure may be effected by either rotation or tangential tilting of TM2.
Collapse
Affiliation(s)
- Mary Hongying Cheng
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | | |
Collapse
|
31
|
Taly A. Opened by a twist: a gating mechanism for the nicotinic acetylcholine receptor. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2007; 36:911-8. [PMID: 17609938 DOI: 10.1007/s00249-007-0189-4] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2007] [Revised: 04/27/2007] [Accepted: 05/07/2007] [Indexed: 12/25/2022]
Affiliation(s)
- Antoine Taly
- Unité de Recherche Associée Centre National de la Recherche Scientifique 2182 "Récepteurs et Cognition", Département de Neuroscience, Institut Pasteur, 25 Rue du Dr Roux, 75724, Paris Cedex 15, France.
| |
Collapse
|
32
|
Cheng X, Ivanov I, Wang H, Sine SM, McCammon JA. Nanosecond-timescale conformational dynamics of the human alpha7 nicotinic acetylcholine receptor. Biophys J 2007; 93:2622-34. [PMID: 17573436 PMCID: PMC1989720 DOI: 10.1529/biophysj.107.109843] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
We explore the conformational dynamics of a homology model of the human alpha7 nicotinic acetylcholine receptor using molecular dynamics simulation and analyses of root mean-square fluctuations, block partitioning of segmental motion, and principal component analysis. The results reveal flexible regions and concerted global motions of the subunits encompassing extracellular and transmembrane domains of the subunits. The most relevant motions comprise a bending, hinged at the beta10-M1 region, accompanied by concerted tilting of the M2 helices that widens the intracellular end of the channel. Despite the nanosecond timescale, the observations suggest that tilting of the M2 helices may initiate opening of the pore. The results also reveal direct coupling between a twisting motion of the extracellular domain and dynamic changes of M2. Covariance analysis of interresidue motions shows that this coupling arises through a network of residues within the Cys and M2-M3 loops where Phe135 is stabilized within a hydrophobic pocket formed by Leu270 and Ile271. The resulting concerted motion causes a downward shift of the M2 helices that disrupts a hydrophobic girdle formed by 9' and 13' residues.
Collapse
Affiliation(s)
- Xiaolin Cheng
- Howard Hughes Medical Institute, National Science Foundation Center for Theoretical Biophysics, Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA.
| | | | | | | | | |
Collapse
|
33
|
Zhou H, Singh NJ, Kim KS. Homology modeling and molecular dynamics study of chorismate synthase from Shigella flexneri. J Mol Graph Model 2006; 25:434-41. [PMID: 16616535 DOI: 10.1016/j.jmgm.2006.02.013] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2005] [Revised: 02/17/2006] [Accepted: 02/17/2006] [Indexed: 11/20/2022]
Abstract
Shigellosis is a major public health problem in many developing countries. Antibiotic therapy can reduce the severity of the dysentery and prevent potentially lethal complication. However, owing to the increased resistance to most of the widely used and inexpensive antibiotics, there is an urgent need for new antibacterial agents, particularly those that act on novel targets. Chorismate synthase (CS) is a key enzyme in the shikimic acid pathway, which is essential for the synthesis of aromatic amino acids in bacteria. As an anti-bacterial drug target, CS has been well validated. A homology model of Shigella-CS with the flavin mononucleotide (FMN) binding was constructed using the crystal structure of CS from other species. The substrate 5-enolpyruvylshikimate 3-phosphate (EPSP) was subsequently docked into the active site based on previous theoretical studies. Molecular dynamics (MD) was used to refine the starting ternary model. The model was well conserved during the 1.8 ns MD simulation with the equilibrium root mean square deviation (RMSD) value of 3.5 angstrom. The substrate binding energy was calculated and the electrostatic energy was found to be the most important term for binding. Decomposition of binding energies revealed that R129, R125, R327, R134 and R48 are important residues involved in substrate binding, which is useful for further site-directed mutagenesis experiments. In the absence of crystal structure, our study provides an early insight into the structure of CS from Shigella flexneri and its binding to the substrate and cofactor, thus facilitating the inhibitor design.
Collapse
Affiliation(s)
- Hong Zhou
- Department of Chemistry, Division of Molecular and Life Sciences, Pohang University of Science and Technology, San 31, Hyojadong, Namgu, Pohang 790-784, Republic of Korea
| | | | | |
Collapse
|
34
|
Zhou H, Singh NJ, Kim KS. Homology modeling and molecular dynamics study of West Nile virus NS3 protease: A molecular basis for the catalytic activity increased by the NS2B cofactor. Proteins 2006; 65:692-701. [PMID: 16972281 DOI: 10.1002/prot.21129] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The West Nile virus (WNV) NS3 serine protease, which plays an important role in assembly of infective virion, is an attractive target for anti-WNV drug development. Cofactors NS2B and NS4A increase the catalytic activity of NS3 in dengue virus and Hepatitis C virus, respectively. Recent studies on the WNV-NS3 characterize the catalytically active form of NS3 by tethering the 40-residue cofactor NS2B. It is suggested that NS2B is essential for the NS3 activity in WNV, while there is no information of the WNV-NS3-related crystal structure. To understand the role of NS2B/substrate in the NS3 catalytic activity, we built a series of models: WNV-NS3 and WNV-NS3-NS2B and WNV-NS3-NS2B-substrate using homology modeling and molecular modeling techniques. Molecular dynamics (MD) simulations were performed for 2.75 ns on each model, to investigate the structural stabilization and catalytic triad motion of the WNV NS3 protease with and without NS2B/substrate. The simulations show that the NS3 rearrangement occurs upon the NS2B binding, resulting in the stable D75-OD1...H51-NH hydrogen bonding. After the substrate binds to the NS3-NS2B active site, the NS3 protease becomes more stable, and the catalytic triad is formed. These results provide a structural basis for the activation and stabilization of the enzyme by its cofactor and substrate.
Collapse
Affiliation(s)
- Hong Zhou
- Department of Chemistry, Division of Molecular and Life Sciences, Pohang University of Science and Technology, Hyojadong, Namgu, Pohang 790-784, Korea
| | | | | |
Collapse
|
35
|
Cheng X, Wang H, Grant B, Sine SM, McCammon JA. Targeted molecular dynamics study of C-loop closure and channel gating in nicotinic receptors. PLoS Comput Biol 2006; 2:e134. [PMID: 17009865 PMCID: PMC1584325 DOI: 10.1371/journal.pcbi.0020134] [Citation(s) in RCA: 105] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2006] [Accepted: 08/23/2006] [Indexed: 11/24/2022] Open
Abstract
The initial coupling between ligand binding and channel gating in the human α7 nicotinic acetylcholine receptor (nAChR) has been investigated with targeted molecular dynamics (TMD) simulation. During the simulation, eight residues at the tip of the C-loop in two alternating subunits were forced to move toward a ligand-bound conformation as captured in the crystallographic structure of acetylcholine binding protein (AChBP) in complex with carbamoylcholine. Comparison of apo- and ligand-bound AChBP structures shows only minor rearrangements distal from the ligand-binding site. In contrast, comparison of apo and TMD simulation structures of the nAChR reveals significant changes toward the bottom of the ligand-binding domain. These structural rearrangements are subsequently translated to the pore domain, leading to a partly open channel within 4 ns of TMD simulation. Furthermore, we confirmed that two highly conserved residue pairs, one located near the ligand-binding pocket (Lys145 and Tyr188), and the other located toward the bottom of the ligand-binding domain (Arg206 and Glu45), are likely to play important roles in coupling agonist binding to channel gating. Overall, our simulations suggest that gating movements of the α7 receptor may involve relatively small structural changes within the ligand-binding domain, implying that the gating transition is energy-efficient and can be easily modulated by agonist binding/unbinding. Nicotinic acetylcholine receptors are ligand-gated ion channels responsible for neurotransmitter-mediated signal transduction at synapses throughout the central and peripheral nervous systems. Binding of neurotransmitter molecules to subunit interfaces in the N-terminal extracellular domain induces structural rearrangements of the membrane-spanning domain permitting the influx of cations. A full understanding of how the conformational changes propagate from the ligand-binding site to the pore domain is of great interest to biologists, yet remains to be established. Using a special simulation technique known as targeted molecular dynamics, Cheng and colleagues probed the early stages of ligand-induced conformational rearrangements that may lead to channel opening. During the simulation, Cheng et al. observed a sequence of conformational changes that stem from the ligand-binding site to the transmembrane domain resulting in a wider channel. From these results, they suggest that gating movements may entail only small structural changes in the ligand-binding domain, implying that channel gating is energy-efficient and can readily be modulated by the binding/unbinding of agonist molecules.
Collapse
Affiliation(s)
- Xiaolin Cheng
- Howard Hughes Medical Institute, University of California San Diego, La Jolla, California, United States of America.
| | | | | | | | | |
Collapse
|
36
|
Khalid S, Bond PJ, Deol SS, Sansom MSP. Modeling and simulations of a bacterial outer membrane protein: OprF from Pseudomonas aeruginosa. Proteins 2006; 63:6-15. [PMID: 16397890 DOI: 10.1002/prot.20845] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
OprF is a major outer membrane protein from Pseudomonas aeruginosa, a homolog of OmpA from Escherichia coli. The N-terminal domains of both proteins have been demonstrated to form low conductance channels in lipid bilayers. Homology models, consisting of an eight-stranded beta-barrel, of the N-terminal domain OprF have been constructed based on the crystal structure of the corresponding domain from E. coli OmpA. OprF homology models have been evaluated via a set (6 x 10 ns) of simulations of the beta-barrel embedded within a solvated dimyristoyl-phosphatidylcholine (DMPC) bilayer. The conformational stability of the models is similar to that of the crystal structure of OmpA in comparable simulations. There is a degree of water penetration into the pore-like center of the OprF barrel. The presence of an acidic/basic (E8/K121) side-chain interaction within the OprF barrel may form a "gate" able to close/open a central pore. Lipid-protein interactions within the simulations were analyzed and revealed that aromatic side-chains (Trp, Tyr) of OprF interact with lipid headgroups. Overall, the behavior of the OprF model in simulations supports the suggestion that this molecule is comparable to OmpA. The simulations help to explain the mechanism of formation of low conductance pores within the outer membrane.
Collapse
Affiliation(s)
- Syma Khalid
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | | | | | | |
Collapse
|
37
|
Corry B. An energy-efficient gating mechanism in the acetylcholine receptor channel suggested by molecular and Brownian dynamics. Biophys J 2005; 90:799-810. [PMID: 16284265 PMCID: PMC1367105 DOI: 10.1529/biophysj.105.067868] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Acetylcholine receptors mediate electrical signaling between nerve and muscle by opening and closing a transmembrane ion conductive pore. Molecular and Brownian dynamics simulations are used to shed light on the location and mechanism of the channel gate. Four separate 5 ns molecular dynamics simulations are carried out on the imaged structure of the channel, a hypothetical open structure with a slightly wider pore and a mutant structure in which a central ring of hydrophobic residues is replaced by polar groups. Water is found to partially evacuate the pore during molecular simulations of the imaged structure, whereas ions face a large energy barrier and do not conduct through the channel in Brownian dynamics simulations. The pore appears to be in a closed configuration despite containing an unobstructed pathway across the membrane as a series of hydrophobic residues in the center of the channel provide an unfavorable home to water and ions. When the channel is widened slightly, water floods into the channel and ions conduct at a rate comparable to the currents measured experimentally in open channels. The pore remains permeable to ions provided the extracellular end of the pore-lining helix is restrained near the putative open configuration to mimic the presence of the ligand binding domain. Replacing some of the hydrophobic residues with polar ones decreases the barrier for ion permeation but does not result in significant currents. The channel is posited to utilize an energy efficient gating mechanism in which only minor conformational changes of the hydrophobic region of the pore are required to create macroscopic changes in conductance.
Collapse
Affiliation(s)
- Ben Corry
- Chemistry, School of Biomedical, Biomolecular and Chemical Sciences, The University of Western Australia, Crawley, WA 6009, Australia.
| |
Collapse
|