1
|
Albano JMR, Mussini N, Toriano R, Facelli JC, Ferraro MB, Pickholz M. Calcium interactions with Cx26 hemmichannel: Spatial association between MD simulations biding sites and variant pathogenicity. Comput Biol Chem 2018; 77:331-342. [PMID: 30466042 DOI: 10.1016/j.compbiolchem.2018.11.004] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2018] [Revised: 07/08/2018] [Accepted: 11/08/2018] [Indexed: 01/23/2023]
Abstract
Connexinophaties are a collective of diseases related to connexin channels and hemichannels. In particular many Cx26 alterations are strongly associated to human deafness. Calcium plays an important role on this structures regulation. Here, using calcium as a probe, extensive atomistic Molecular Dynamics simulations were performed on the Cx26 hemichannel embedded in a lipid bilayer. Exploring different initial conditions and calcium concentration, simulation reached ∼4 μs. Several analysis were carried out in order to reveal the calcium distribution and localization, such as electron density profiles, density maps and distance time evolution, which is directly associated to the interaction energy. Specific amino acid interactions with calcium and their stability were capture within this context. Few of these sites such as, GLU42, GLU47, GLY45 and ASP50, were already suggested in the literature. Besides, we identified novel calcium biding sites: ASP2, ASP117, ASP159, GLU114, GLU119, GLU120 and VAL226. To the best of our knowledge, this is the first time that these sites are reported within this context. Furthermore, since various pathologies involving the Cx26 hemichannel are associated with pathogenic variants in the corresponding CJB2 gene, using ClinVar, we were able to spatially associate the 3D positions of the identified calcium binding sites within the framework of this work with reported pathogenic variants in the CJB2 gene. This study presents a first step on finding associations between molecular features and pathological variants of the Cx26 hemichannel.
Collapse
Affiliation(s)
- Juan M R Albano
- Facultad de Ciencias Exactas y Naturales, Departamento de Física, Universidad de Buenos Aires, Argentina; CONICET- Universidad de Buenos Aires, IFIBA, Buenos Aires, Argentina
| | - Nahuel Mussini
- Facultad de Ciencias Exactas y Naturales, Departamento de Física, Universidad de Buenos Aires, Argentina; CONICET- Universidad de Buenos Aires, IFIBA, Buenos Aires, Argentina
| | - Roxana Toriano
- Facultad de Medicina, Departamento de Ciencias Fisiológicas, Laboratorio de Biomembranas, Buenos Aires, Argentina; CONICET - Universidad de Buenos Aires, IFIBIO Houssay, Buenos Aires, Argentina
| | - Julio C Facelli
- Department of Biomedical Informatics, The University of Utah, 421 Wakara Way, Suite 140, Salt Lake City, UT 84108, USA.
| | - Marta B Ferraro
- Facultad de Ciencias Exactas y Naturales, Departamento de Física, Universidad de Buenos Aires, Argentina; CONICET- Universidad de Buenos Aires, IFIBA, Buenos Aires, Argentina
| | - Mónica Pickholz
- Facultad de Ciencias Exactas y Naturales, Departamento de Física, Universidad de Buenos Aires, Argentina; CONICET- Universidad de Buenos Aires, IFIBA, Buenos Aires, Argentina
| |
Collapse
|
2
|
Liu L, Patey GN. A molecular dynamics investigation of the influence of water structure on ion conduction through a carbon nanotube. J Chem Phys 2017; 146:074502. [DOI: 10.1063/1.4975690] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- L. Liu
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| | - G. N. Patey
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| |
Collapse
|
3
|
Ohba T, Yamamoto S, Kodaira T, Hata K. Changing water affinity from hydrophobic to hydrophilic in hydrophobic channels. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2015; 31:1058-1063. [PMID: 25586050 DOI: 10.1021/la504522x] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The behavior of water at hydrophobic interfaces can play a significant role in determining chemical reaction outcomes and physical properties. Carbon nanotubes and aluminophosphate materials have one-dimensional hydrophobic channels, which are entirely surrounded by hydrophobic interfaces. Unique water behavior was observed in such hydrophobic channels. In this article, changes in the water affinity in one-dimensional hydrophobic channels were assessed using water vapor adsorption isotherms at 303 K and grand canonical Monte Carlo simulations. Hydrophobic behavior of water adsorbed in channels wider than 3 nm was observed for both adsorption and desorption processes, owing to the hydrophobic environment. However, water showed hydrophilic properties in both adsorption and desorption processes in channels narrower than 1 nm. In intermediate-sized channels, the hydrophobic properties of water during the adsorption process were seen to transition to hydrophilic behavior during the desorption process. Hydrophilic properties in the narrow channels for both adsorption and desorption processes are a result of the relatively strong water-channel interactions (10-15 kJ mol(-1)). In the 2-3 nm channels, the water-channel interaction energy of 4-5 kJ mol(-1) was comparable to the thermal translational energy. The cohesive water interaction was approximately 35 kJ mol(-1), which was larger than the others. Thus, the water affinity change in the 2-3 nm channels for the adsorption and desorption processes was attributed to weak water-channel interactions and strong cohesive interactions. These results are inherently important to control the properties of water in hydrophobic environments.
Collapse
Affiliation(s)
- Tomonori Ohba
- Graduate School of Science, Chiba University , 1-33 Yayoi, Inage, Chiba 263-8522, Japan
| | | | | | | |
Collapse
|
4
|
Kulleperuma K, Smith SME, Morgan D, Musset B, Holyoake J, Chakrabarti N, Cherny VV, DeCoursey TE, Pomès R. Construction and validation of a homology model of the human voltage-gated proton channel hHV1. ACTA ACUST UNITED AC 2013; 141:445-65. [PMID: 23530137 PMCID: PMC3607825 DOI: 10.1085/jgp.201210856] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The topological similarity of voltage-gated proton channels (HV1s) to the voltage-sensing domain (VSD) of other voltage-gated ion channels raises the central question of whether HV1s have a similar structure. We present the construction and validation of a homology model of the human HV1 (hHV1). Multiple structural alignment was used to construct structural models of the open (proton-conducting) state of hHV1 by exploiting the homology of hHV1 with VSDs of K+ and Na+ channels of known three-dimensional structure. The comparative assessment of structural stability of the homology models and their VSD templates was performed using massively repeated molecular dynamics simulations in which the proteins were allowed to relax from their initial conformation in an explicit membrane mimetic. The analysis of structural deviations from the initial conformation based on up to 125 repeats of 100-ns simulations for each system reveals structural features consistently retained in the homology models and leads to a consensus structural model for hHV1 in which well-defined external and internal salt-bridge networks stabilize the open state. The structural and electrostatic properties of this open-state model are compatible with proton translocation and offer an explanation for the reversal of charge selectivity in neutral mutants of Asp112. Furthermore, these structural properties are consistent with experimental accessibility data, providing a valuable basis for further structural and functional studies of hHV1. Each Arg residue in the S4 helix of hHV1 was replaced by His to test accessibility using Zn2+ as a probe. The two outermost Arg residues in S4 were accessible to external solution, whereas the innermost one was accessible only to the internal solution. Both modeling and experimental data indicate that in the open state, Arg211, the third Arg residue in the S4 helix in hHV1, remains accessible to the internal solution and is located near the charge transfer center, Phe150.
Collapse
Affiliation(s)
- Kethika Kulleperuma
- Molecular Structure and Function, Hospital for Sick Children, Toronto, Ontario M5G 1X8, Canada
| | | | | | | | | | | | | | | | | |
Collapse
|
5
|
Maffeo C, Bhattacharya S, Yoo J, Wells D, Aksimentiev A. Modeling and simulation of ion channels. Chem Rev 2012; 112:6250-84. [PMID: 23035940 PMCID: PMC3633640 DOI: 10.1021/cr3002609] [Citation(s) in RCA: 148] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Christopher Maffeo
- Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL
| | - Swati Bhattacharya
- Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL
| | - Jejoong Yoo
- Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL
| | - David Wells
- Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL
| | - Aleksei Aksimentiev
- Department of Physics, University of Illinois, 1110 W. Green St., Urbana, IL
| |
Collapse
|
6
|
Jung YW, Lu B, Mascagni M. A computational study of ion conductance in the KcsA K+ channel using a Nernst–Planck model with explicit resident ions. J Chem Phys 2009; 131:215101. [DOI: 10.1063/1.3268774] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
7
|
Kraszewski S, Boiteux C, Ramseyer C, Girardet C. Determination of the charge profile in the KcsA selectivity filter using ab initio calculations and molecular dynamics simulations. Phys Chem Chem Phys 2009; 11:8606-13. [PMID: 19774294 DOI: 10.1039/b905991a] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The charge profile of K+ and Na+ ions moving in a single file across the filter of the KcsA channel is determined using both molecular dynamics simulations and ab initio calculations. We show a strong correlation between the charge variation and the ion location resulting in a saw-tooth profile, which provides additional information on the influence of charge transfer on the permeation and selectivity of the channel.
Collapse
Affiliation(s)
- Sebastian Kraszewski
- Laboratoire de Physique Moléculaire-UMR CNRS 6624, Faculté des Sciences Université de Franche-Comté, 16, Route de Gray, 25030, Besançon Cedex, France.
| | | | | | | |
Collapse
|
8
|
Shao Q, Zhou J, Lu L, Lu X, Zhu Y, Jiang S. Anomalous hydration shell order of Na+ and K+ inside carbon nanotubes. NANO LETTERS 2009; 9:989-994. [PMID: 19206198 DOI: 10.1021/nl803044k] [Citation(s) in RCA: 64] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
We performed molecular dynamics simulations of the hydration of Na+ and K+ in infinitely long single-walled armchair carbon nanotubes (CNTs) at 298 K. Simulation results indicate that the preferential orientation of water molecules in coordination shells of these two cations presents an anomalous change in the CNTs and causes a diameter-dependent variation for the interaction energy between the cation and water molecules in its coordination shell. In the five CNTs of this work, it is energetically favorable for confining a hydrated K+ inside the two narrow CNTs with diameters of 0.60 and 0.73 nm, whereas the situation is reverse inside the wide CNTs with diameters of 0.87, 1.0, and 1.28 nm. This finding is important for CNT applications in ionic systems that control the selectivity and the ionic flow.
Collapse
Affiliation(s)
- Qing Shao
- State Key Laboratory of Materials-oriented Chemical Engineering, Nanjing University of Technology, Nanjing 210009, People's Republic of China
| | | | | | | | | | | |
Collapse
|
9
|
Jørgensen AM, Tagmose L, Jørgensen AMM, Bøgesø KP, Peters GH. Molecular dynamics simulations of Na+/Cl(-)-dependent neurotransmitter transporters in a membrane-aqueous system. ChemMedChem 2008; 2:827-40. [PMID: 17436258 DOI: 10.1002/cmdc.200600243] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
We have performed molecular dynamics simulations of a homology model of the human serotonin transporter (hSERT) in a membrane environment and in complex with either the natural substrate 5-HT or the selective serotonin reuptake inhibitor escitalopram. We have also included a transporter homologue, the Aquifex aeolicus leucine transporter (LeuT), in our study to evaluate the applicability of a simple and computationally attractive membrane system. Fluctuations in LeuT extracted from simulations are in good agreement with crystallographic B factors. Furthermore, key interactions identified in the X-ray structure of LeuT are maintained throughout the simulations indicating that our simple membrane system is suitable for studying the transmembrane protein hSERT in complex with 5-HT or escitalopram. For these transporter complexes, only relatively small fluctuations are observed in the ligand-binding cleft. Specific interactions responsible for ligand recognition, are identified in the hSERT-5HT and hSERT-escitalopram complexes. Our findings are in good agreement with predictions from mutagenesis studies.
Collapse
Affiliation(s)
- Anne Marie Jørgensen
- MEMPHYS-Center for Biomembrane Physics, Department of Chemistry, Technical University of Denmark, Building 206, 2800 Kgs. Lyngby, Denmark
| | | | | | | | | |
Collapse
|
10
|
Kraszewski S, Yesylevskyy SO, Boiteux C, Ramseyer C, Kharkyanen VN. Is the mobility of the pore walls and water molecules in the selectivity filter of KcsA channel functionally important? Phys Chem Chem Phys 2008; 10:2249-55. [PMID: 18404233 DOI: 10.1039/b719512e] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We performed in-depth analysis of the forces which act on the K(+) ions in the selectivity filter of the KcsA channel in order to estimate the relative importance of static and dynamic influence of the filter wall and water molecules on ion permeation and selectivity. The forces were computed using the trajectories of all-atom molecular dynamics simulations. It is shown that the dynamics of the selectivity filter contributes about 3% to the net force acting on the ions and can be neglected in the studies focused on the macroscopic properties of the channel, such as the current. Among the filter atoms, only the pore-forming carbonyl groups can be considered as dynamic in the studies of microscopic events of conduction, while the dynamic effects from all other atoms are negligible. We also show that the dynamics of the water molecules in the filter can not be neglected. The fluctuating forces from the water molecules can be as strong as net forces from the pore walls and can effectively drive the ions through the local energy barriers in the filter.
Collapse
Affiliation(s)
- Sebastian Kraszewski
- Institut UTINAM, Laboratoire de Physique Moléculaire, UMR CNRS 6213, Faculté des Sciences et Techniques, Université de Franche-Comté, 16 route de Gray, 25030 Besançon Cedex, La Bouloie, France
| | | | | | | | | |
Collapse
|
11
|
Molecular dynamics simulations of asymmetric NaCl and KCl solutions separated by phosphatidylcholine bilayers: potential drops and structural changes induced by strong Na+-lipid interactions and finite size effects. Biophys J 2008; 94:3565-76. [PMID: 18222999 DOI: 10.1529/biophysj.107.116335] [Citation(s) in RCA: 99] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Differences of ionic concentrations across lipid bilayers are some of the primary energetic driving forces for cellular electrophysiology. While macroscopic models of asymmetric ionic solutions are well-developed, their connection to ion, water, and lipid interactions at the atomic scale are much more poorly understood. In this study, we used molecular dynamics to examine a system of two chambers of equal ionic strength, but differing amounts of NaCl and KCl, separated by a lipid bilayer. Our expectation was that the net electrostatic potential difference between the two chambers should be small or zero. Contrary to our expectation, a large potential difference (-70 mV) slowly evolved across the two water chambers over the course of our 172-ns simulation. This potential primarily originated from strong Na(+) binding to the carbonyls of the phosphatidylcholine lipids. This ion adsorption also led to significant structural and mechanical changes in the lipid bilayer. We discuss this surprising result in the context of indirect experimental evidence for Na(+) interaction with bilayers as well as potential caveats in current biomembrane simulation methodology, including force-field parameters and finite size effects.
Collapse
|
12
|
Garbuio V, Andreani C, Imberti S, Pietropaolo A, Reiter GF, Senesi R, Ricci MA. Proton quantum coherence observed in water confined in silica nanopores. J Chem Phys 2007; 127:154501. [DOI: 10.1063/1.2789436] [Citation(s) in RCA: 61] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Abstract
AbstractDespite the complexity of ion-channels, MD simulations based on realistic all-atom models have become a powerful technique for providing accurate descriptions of the structure and dynamics of these systems, complementing and reinforcing experimental work. Successful multidisciplinary collaborations, progress in the experimental determination of three-dimensional structures of membrane proteins together with new algorithms for molecular simulations and the increasing speed and availability of supercomputers, have made possible a considerable progress in this area of biophysics. This review aims at highlighting some of the work in the area of potassium channels and molecular dynamics simulations where numerous fundamental questions about the structure, function, folding and dynamics of these systems remain as yet unresolved challenges.
Collapse
|
14
|
Haider S, Khalid S, Tucker SJ, Ashcroft FM, Sansom MSP. Molecular dynamics simulations of inwardly rectifying (Kir) potassium channels: a comparative study. Biochemistry 2007; 46:3643-52. [PMID: 17326663 DOI: 10.1021/bi062210f] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Inward rectifier potassium (Kir) channels regulate cell excitability and transport K+ ions across membranes. Homotetrameric models of three mammalian Kir channels (Kir1.1, Kir3.1, and Kir6.2) have been generated, using the KirBac3.1 transmembrane and rat Kir3.1 intracellular domain structures as templates. All three models have been explored by 10 ns molecular dynamics simulations in phospholipid bilayers. Analysis of the initial structures revealed conservation of potential lipid interaction residues (Trp/Tyr and Arg/Lys side chains near the lipid headgroup-water interfaces). Examination of the intracellular domains revealed key structural differences between Kir1.1 and Kir6.2 which may explain the difference in channel inhibition by ATP. The behavior of all three models in the MD simulations revealed that they have conformational stability similar to that seen for comparable simulations of, for example, structures derived from cryoelectron microscopy data. Local distortions of the selectivity filter were seen during the simulations, as observed in previous simulations of KirBac and in simulations and structures of KcsA. These may be related to filter gating of the channel. The intracellular hydrophobic gate does not undergo any substantial changes during the simulations and thus remains functionally closed. Analysis of lipid-protein interactions of the Kir models emphasizes the key role of the M0 (or "slide") helix which lies approximately parallel to the bilayer-water interface and forms a link between the transmembrane and intracellular domains of the channel.
Collapse
Affiliation(s)
- Shozeb Haider
- Department of Biochemistry, University of Oxford, UK
| | | | | | | | | |
Collapse
|
15
|
Kóňa J, Minozzi M, Torre V, Carloni P. A gate mechanism indicated in the selectivity filter of the potassium channel KscA. Theor Chem Acc 2007. [DOI: 10.1007/s00214-006-0226-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
16
|
Faraldo-Gómez JD, Forrest LR, Baaden M, Bond PJ, Domene C, Patargias G, Cuthbertson J, Sansom MSP. Conformational sampling and dynamics of membrane proteins from 10-nanosecond computer simulations. Proteins 2006; 57:783-91. [PMID: 15317024 DOI: 10.1002/prot.20257] [Citation(s) in RCA: 81] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
In the current report, we provide a quantitative analysis of the convergence of the sampling of conformational space accomplished in molecular dynamics simulations of membrane proteins of duration in the order of 10 nanoseconds. A set of proteins of diverse size and topology is considered, ranging from helical pores such as gramicidin and small beta-barrels such as OmpT, to larger and more complex structures such as rhodopsin and FepA. Principal component analysis of the C(alpha)-atom trajectories was employed to assess the convergence of the conformational sampling in both the transmembrane domains and the whole proteins, while the time-dependence of the average structure was analyzed to obtain single-domain information. The membrane-embedded regions, particularly those of small or structurally simple proteins, were found to achieve reasonable convergence. By contrast, extra-membranous domains lacking secondary structure are often markedly under-sampled, exhibiting a continuous structural drift. This drift results in a significant imprecision in the calculated B-factors, which detracts from any quantitative comparison to experimental data. In view of such limitations, we suggest that similar analyses may be valuable in simulation studies of membrane protein dynamics, in order to attach a level of confidence to any biologically relevant observations.
Collapse
Affiliation(s)
- José D Faraldo-Gómez
- Department of Physiology and Biophysics, Weill Medical College of Cornell University, 1300 York Avenue, New York, New York 10021, USA.
| | | | | | | | | | | | | | | |
Collapse
|
17
|
Grottesi A, Domene C, Hall B, Sansom MSP. Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations. Biochemistry 2006; 44:14586-94. [PMID: 16262258 DOI: 10.1021/bi0510429] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.
Collapse
|
18
|
Law RJ, Capener C, Baaden M, Bond PJ, Campbell J, Patargias G, Arinaminpathy Y, Sansom MSP. Membrane protein structure quality in molecular dynamics simulation. J Mol Graph Model 2006; 24:157-65. [PMID: 16102990 DOI: 10.1016/j.jmgm.2005.05.006] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/04/2005] [Indexed: 10/25/2022]
Abstract
Our goal was to assess the relationship between membrane protein quality, output from protein quality checkers and output from molecular dynamics (MD) simulations. Membrane transport proteins are essential for a wide range of cellular processes. Structural features of integral membrane proteins are still under-explored due to experimental limitations in structure determination. Computational techniques can be used to exploit biochemical and medium resolution structural data, as well as sequence homology to known structures, and enable us to explore the structure-function relationships in several transmembrane proteins. The quality of the models produced is vitally important to obtain reliable predictions. An examination of the relationship between model stability in molecular dynamics (MD) simulations derived from RMSD (root mean squared deviation) and structure quality assessment from various protein quality checkers was undertaken. The results were compared to membrane protein structures, solved at various resolution, by either X-ray or electron diffraction techniques. The checking programs could predict the potential success of MD in making functional conclusions. MD stability was shown to be a good indicator for the quality of structures. The quality was also shown to be dependent on the resolution at which the structures were determined.
Collapse
Affiliation(s)
- Richard J Law
- Laboratory of Molecular Biophysics, Department of Biochemistry, The University of Oxford, The Rex Richards Building, South Parks Road, Oxford OX1 3QU, UK.
| | | | | | | | | | | | | | | |
Collapse
|
19
|
Candler A, Featherstone M, Ali R, Maloney L, Watts A, Fischer WB. Computational analysis of mutations in the transmembrane region of Vpu from HIV-1. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2005; 1716:1-10. [PMID: 16154109 DOI: 10.1016/j.bbamem.2005.07.012] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/16/2005] [Revised: 06/27/2005] [Accepted: 07/28/2005] [Indexed: 11/27/2022]
Abstract
Vpu is an 81 amino acid integral membrane protein encoded by HIV-1. Its alpha-helical transmembrane (TM) domain (residues approximately 6-28) enhances virion release by oligomerizing into bundles and forming ion-conducting channels across the plasma membrane. Its cytoplasmic domain (residues approximately 29-81) is also alpha-helical and binds to the transmembrane protein CD4, inducing its degradation. Mutations within the TM domain have been found to abrogate enhanced particle release from the infected cell (Tiganos et al. Virology (1998) 251 96-107). A series of computational models of monomeric, pentameric and hexameric Vpu(1-31) mutants have been constructed, embedded in fully hydrated lipid bilayers and subjected to a 3 ns molecular dynamics (MD) simulation. None of the mutations has any destabilizing effect on the secondary and tertiary structure. One of the mutants, in which the position of a tryptophan residue within the TM domain is altered, is known not to induce CD4 degradation; an extended kinked model of this mutant has been generated (Vpu(1-52)IVW-k) and during subsequent MD simulations, the bend between the TM and a part of the cytoplasmic domain is found to unwind and a complex salt bridge involving Lys-37 is formed.
Collapse
Affiliation(s)
- Andrew Candler
- Biomembrane Structure Unit, Department of Biochemistry, Oxford University, Oxford OX1 3QU, UK
| | | | | | | | | | | |
Collapse
|
20
|
Abstract
Organ function (the heart beat for example) can only be understood through knowledge of molecular and cellular processes within the constraints of structure-function relations at the tissue level. A quantitative modeling framework that can deal with these multiscale issues is described here under the banner of the International Union of Physiological Sciences Physiome Project.
Collapse
Affiliation(s)
- Peter Hunter
- Bioengineering Institute, University of Auckland, Auckland, New Zealand.
| | | |
Collapse
|
21
|
Grottesi A, Domene C, Haider S, Sansom MSP. Molecular dynamics simulation approaches to K channels: conformational flexibility and physiological function. IEEE Trans Nanobioscience 2005; 4:112-20. [PMID: 15816177 DOI: 10.1109/tnb.2004.842473] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Molecular modeling and simulations enable extrapolation for the structure of bacterial potassium channels to the function of their mammalian homologues. Molecular dynamics simulations have revealed the concerted single-file motion of potassium ions and water molecules through the selectivity filter of K channels and the role of filter flexibility in ion permeation and in "fast gating." Principal components analysis of extended K channel simulations suggests that hinge-bending of pore-lining M2 (or S6) helices plays a key role in K channel gating. Based on these and other simulations, a molecular model for gating of inward rectifier K channel gating is presented.
Collapse
Affiliation(s)
- Alessandro Grottesi
- Department of Biochemistry, University of Oxford, Oxford OX 3QU, United Kingdom.
| | | | | | | |
Collapse
|
22
|
Compoint M, Picaud F, Ramseyer C, Girardet C. Zip gating of the KcsA channel studied by targeted molecular dynamics. Chem Phys Lett 2005. [DOI: 10.1016/j.cplett.2005.03.087] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
23
|
Compoint M, Picaud F, Ramseyer C, Girardet C. Targeted molecular dynamics of an open-state KcsA channel. J Chem Phys 2005; 122:134707. [PMID: 15847489 DOI: 10.1063/1.1869413] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Pore opening of KcsA channel is studied using targeted molecular dynamics simulations. Conformational changes of the protein are determined, starting from the crystallized refined 2.0 A structure (pdb 1K4C) determined in x-ray experiments and arriving to the open-state structure constructed on the basis of electron paramagnetic resonance spectroscopic data (pdb 1JQ1). Our results corroborate the essential role played by the terminal residues located on the transmembrane helices M2 which were not taken into account at that time. The aperture mechanism of the channel appears to be ziplike. A small constraint (approximately equal to 5 x 10(-2) kcal mol(-1) A(-2) per C(alpha)) applied to the terminal residues located on the intracellular side is sufficient to initialize the pore opening at the innermost part of the gate, but additional constraint must be applied to definitely complete the pore aperture. The open structure is proved to be a metastable state since releasing the constraint leads to another relaxed open conformation which seems to reach stability.
Collapse
Affiliation(s)
- Mylène Compoint
- Laboratoire de Physique Moléculaire, UMR CNRS 6624-Faculté des Sciences, la Bouloie, Université de Franche-Comté, 25030 Besançon Cedex, France
| | | | | | | |
Collapse
|
24
|
Domene C, Grottesi A, Sansom MSP. Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1. Biophys J 2005; 87:256-67. [PMID: 15240462 PMCID: PMC1304348 DOI: 10.1529/biophysj.104.039917] [Citation(s) in RCA: 54] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.
Collapse
Affiliation(s)
- Carmen Domene
- Laboratory of Molecular Biophysics, Department of Biochemistry, University of Oxford, Oxford, OX1 3QU United Kingdom
| | | | | |
Collapse
|
25
|
Compoint M, Boiteux C, Huetz P, Ramseyer C, Girardet C. Role of water molecules in the KcsA protein channel by molecular dynamics calculations. Phys Chem Chem Phys 2005; 7:4138-45. [PMID: 16474879 DOI: 10.1039/b508281a] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Molecular dynamics simulations supported by electrostatic calculations have been conducted on the KcsA channel to determine the role of water molecules in the pore. Starting from the X-ray structure of the KcsA channel in its closed state at 2.0 angstroms resolution, the opening of the pore towards a conformation built on the basis of EPR results is studied. We show that water molecules act as a structural element for the K+ ions inside the filter and the hydrophobic cavity of the channel. In the filter, water tends to enhance the depth of the wells occupied by the K+ ions, while in the cavity there is a strong correlation between the water molecules and the cavity ion. As a consequence, the protein remains very stable in the presence of three K+ ions in the selectivity filter and one in the cavity. The analysis of the dynamics of water molecules in the cavity reveals preferred orientations of the dipoles along the pore axis, and a correlated behavior between this dipole orientation and the displacement of the K+ ion during the gating process.
Collapse
Affiliation(s)
- M Compoint
- Laboratoire de Physique Moléculaire, UMR CNRS 6624, Université de Franche-Comté, F-25030 Besançon Cedex, France
| | | | | | | | | |
Collapse
|
26
|
Park H, Lee S. Prediction of the mutation-induced change in thermodynamic stabilities of membrane proteins from free energy simulations. Biophys Chem 2004; 114:191-7. [PMID: 15829352 DOI: 10.1016/j.bpc.2004.12.006] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2004] [Revised: 12/02/2004] [Accepted: 12/03/2004] [Indexed: 11/17/2022]
Abstract
Comparative protein structure modeling and free energy perturbation simulation have been applied in a consecutive manner to investigate the mutation-induced stabilization of membrane proteins (MPs) in aqueous solution without knowledge of their three-dimensional structures. The calculated difference in protein solvation free energy between the wild type and a mutant compares well with their relative thermodynamic stabilities in solution. For monomeric MPs, a mutant reveals a higher stability than the wild type if the calculated solvation free energy indicates a favorable change. On the contrary, for oligomeric MPs the stability of a mutant increases as the solvation free energy of a mutated monomer becomes less favorable, indicating that the oligomeric MP mutant would be stabilized in solution due to the reduced desolvation cost for oligomerization. The present computational strategy is expected to find its way as a useful tool for assessing the relative stability of a mutant MP with respect to its wild type in solution.
Collapse
Affiliation(s)
- Hwangseo Park
- School of Chemistry and Molecular Engineering, Seoul National University, Korea
| | | |
Collapse
|
27
|
Gouldson PR, Kidley NJ, Bywater RP, Psaroudakis G, Brooks HD, Diaz C, Shire D, Reynolds CA. Toward the active conformations of rhodopsin and the beta2-adrenergic receptor. Proteins 2004; 56:67-84. [PMID: 15162487 DOI: 10.1002/prot.20108] [Citation(s) in RCA: 66] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Using sets of experimental distance restraints, which characterize active or inactive receptor conformations, and the X-ray crystal structure of the inactive form of bovine rhodopsin as a starting point, we have constructed models of both the active and inactive forms of rhodopsin and the beta2-adrenergic G-protein coupled receptors (GPCRs). The distance restraints were obtained from published data for site-directed crosslinking, engineered zinc binding, site-directed spin-labeling, IR spectroscopy, and cysteine accessibility studies conducted on class A GPCRs. Molecular dynamics simulations in the presence of either "active" or "inactive" restraints were used to generate two distinguishable receptor models. The process for generating the inactive and active models was validated by the hit rates, yields, and enrichment factors determined for the selection of antagonists in the inactive model and for the selection of agonists in the active model from a set of nonadrenergic GPCR drug-like ligands in a virtual screen using ligand docking software. The simulation results provide new insights into the relationships observed between selected biochemical data, the crystal structure of rhodopsin, and the structural rearrangements that occur during activation.
Collapse
|
28
|
Treptow W, Maigret B, Chipot C, Tarek M. Coupled motions between pore and voltage-sensor domains: a model for Shaker B, a voltage-gated potassium channel. Biophys J 2004; 87:2365-79. [PMID: 15454436 PMCID: PMC1304659 DOI: 10.1529/biophysj.104.039628] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2004] [Accepted: 06/23/2004] [Indexed: 11/18/2022] Open
Abstract
A high-resolution crystal structure of KvAP, an archeabacterial voltage-gated potassium (Kv) channel, complexed with a monoclonal Fab fragment has been recently determined. Based on this structure, a mechanism for the activation (opening) of Kv channels has been put forward. This mechanism has since been criticized, suggesting that the resolved structure is not representative of the family of voltage-gated potassium channels. Here, we propose a model of the transmembrane domain of Shaker B, a well-characterized Kv channel, built by homology modeling and docking calculations. In this model, the positively charged S4 helices are oriented perpendicular to the membrane and localized in the groove between segments S5 and S6 of adjacent subunits. The structure and the dynamics of the full atomistic model embedded in a hydrated lipid bilayer were investigated by means of two large-scale molecular dynamics simulations under transmembrane-voltage conditions known to induce, respectively, the resting state (closed) and the activation (opening) of voltage-gated channels. Upon activation, the model undergoes conformational changes that lead to an increase of the hydration of the charged S4 helices, correlated with an upward translation and a tilting of the latter, concurrently with movements of the S5 helices and the activation gate. Although small, these conformational changes ultimately result in an alteration of the ion-conduction pathway. Our findings support the transporter model devised by Bezanilla and collaborators, and further underline the crucial role played by internal hydration in the activation of the channel.
Collapse
Affiliation(s)
- Werner Treptow
- Equipe de Dynamique des Assemblages Membranaires, Unité Mixte de Recherche, Centre National de la Recherche Scientifique/Université Henri Poincaré 7565, Nancy, France
| | | | | | | |
Collapse
|
29
|
Compoint M, Ramseyer C, Huetz P. Ab initio investigation of the atomic charges in the KcsA channel selectivity filter. Chem Phys Lett 2004. [DOI: 10.1016/j.cplett.2004.09.033] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
30
|
Jiang FY, Bouret Y, Kindt JT. Molecular dynamics simulations of the lipid bilayer edge. Biophys J 2004; 87:182-92. [PMID: 15240456 PMCID: PMC1304341 DOI: 10.1529/biophysj.103.031054] [Citation(s) in RCA: 79] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2003] [Accepted: 03/29/2004] [Indexed: 11/18/2022] Open
Abstract
Phospholipid bilayers have been intensively studied by molecular dynamics (MD) simulation in recent years. The properties of bilayer edges are important in determining the structure and stability of pores formed in vesicles and biomembranes. In this work, we use molecular dynamics simulation to investigate the structure, dynamics, and line tension of the edges of bilayer ribbons composed of pure dimyristoylphosphatidylcholine (DMPC) or palmitoyl-oleoylphosphatidylethanolamine (POPE). As expected, we observe a significant reorganization of lipids at and near the edges. The treatment of electrostatic effects is shown to have a qualitative impact on the structure and stability of the edge, and significant differences are observed in the dynamics and structure of edges formed by DMPC and palmitoyl-oleoylphosphatidylethanolamine. From the pressure anisotropy in the simulation box, we calculate a line tension of approximately 10-30 pN for the DMPC edge, in qualitative agreement with experimental estimates for similar lipids.
Collapse
Affiliation(s)
- Frank Y Jiang
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | | | | |
Collapse
|
31
|
Law RJ, Sansom MSP. Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2004; 33:477-89. [PMID: 15071758 DOI: 10.1007/s00249-004-0398-z] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2003] [Revised: 02/10/2004] [Accepted: 02/12/2004] [Indexed: 10/26/2022]
Abstract
The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.
Collapse
Affiliation(s)
- Richard J Law
- Laboratory of Molecular Biophysics, Department of Biochemistry, The University of Oxford, South Parks Road, Oxford, OX1 3QU, UK
| | | |
Collapse
|
32
|
Compoint M, Carloni P, Ramseyer C, Girardet C. Molecular dynamics study of the KcsA channel at 2.0-A resolution: stability and concerted motions within the pore. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2004; 1661:26-39. [PMID: 14967472 DOI: 10.1016/j.bbamem.2003.11.019] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/21/2003] [Revised: 11/20/2003] [Accepted: 11/24/2003] [Indexed: 11/26/2022]
Abstract
The stability of the KcsA channel accommodating more than one ion in the pore has been studied with molecular dynamics. We have used the very last X-ray structure of the KcsA channel at 2.0-A resolution determined by Zhou et al. [Nature 414 (2001) 43]. In this channel, six of the seven experimentally evidenced sites have been considered. We show that the protein remains very stable in the presence of four K+ ions (three in the selectivity filter and one in the cavity). The locations and the respective distances of the different K+ ions and water molecules (W), calculated within our KWKWKK sequence, also fits well with the experimental observations. The analysis of the K+ ions and water molecules displacements shows concerted file motions on the simulated time scale (approximately 1 ns), which could act as precursor to the diffusion of K+ ions inside the channel. A simple one-dimensional dynamical model is used to interpret the concerted motions of the ions and water molecules in the pore leading ultimately to ion transfer.
Collapse
Affiliation(s)
- Mylène Compoint
- Laboratoire de Physique Moléculaire, UMR CNRS 6624, Faculté des Sciences, la Bouloie, Université de Franche, Comté, 25030, Besançon, France.
| | | | | | | |
Collapse
|
33
|
Nelson PH. Modeling the concentration-dependent permeation modes of the KcsA potassium ion channel. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2003; 68:061908. [PMID: 14754235 DOI: 10.1103/physreve.68.061908] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/16/2003] [Revised: 08/18/2003] [Indexed: 05/24/2023]
Abstract
The potassium channel from Streptomyces lividans (KcsA) is an integral membrane protein with sequence similarity to all known potassium channels, particularly in the selectivity filter region. A recently proposed model for ion channels containing either n or (n-1) single-file ions in their selectivity filters [P. H. Nelson, J. Chem. Phys. 177, 11396 (2002)] is applied to published KcsA channel K+ permeation data that exhibit a high-affinity process at low concentrations and a low-affinity process at high concentrations [M. LeMasurier et al., J. Gen. Physiol. 118, 303 (2001)]. The kinetic model is shown to provide a reasonable first-order explanation for both the high- and low-concentration permeation modes observed experimentally. The low-concentration mode ([K+]<200 mM) has a 200-mV dissociation constant of 56 mM and a conductance of 88 pS. The high-concentration mode ([K+]>200 mM) has a 200-mV dissociation constant of 1100 mM and a conductance of 500 pS. Based on the permeation model, and x-ray analysis [J. H. Morais-Cabral et al., Nature (London) 414, 37 (2001)], it is suggested that the experimentally observed K+ permeation modes correspond to an n=3 mechanism at high concentrations and an n=2 mechanism at low concentrations. The ratio of the electrical dissociation distances for the high- and low-concentration modes is 3:2, also consistent with the proposed n=3 and n=2 modes. Model predictions for K+ channels that exhibit asymmetric current-voltage (I-V) curves are presented, and further validation of the kinetic model via molecular simulation and experiment is discussed. The qualitatively distinct I-V characteristics exhibited experimentally by Tl+, NH+4, and Rb+ ions at 100 mM concentration can also be explained using the model, but more extensive experimental tests are required for quantitative validation of the model predictions.
Collapse
Affiliation(s)
- Peter Hugo Nelson
- Department of Physics, Benedictine University, Lisle, Illinois 60532, USA.
| |
Collapse
|
34
|
Domene C, Sansom MSP. Potassium channel, ions, and water: simulation studies based on the high resolution X-ray structure of KcsA. Biophys J 2003; 85:2787-800. [PMID: 14581184 PMCID: PMC1303560 DOI: 10.1016/s0006-3495(03)74702-x] [Citation(s) in RCA: 97] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2002] [Accepted: 04/15/2003] [Indexed: 11/25/2022] Open
Abstract
Interactions of Na(+), K(+), Rb(+), and Cs(+) ions within the selectivity filter of a potassium channel have been investigated via multiple molecular dynamics simulations (total simulation time, 48 ns) based on the high resolution structure of KcsA, embedded in a phospholipid bilayer. As in simulations based on a lower resolution structure of KcsA, concerted motions of ions and water within the filter are seen. Despite the use of a higher resolution structure and the inclusion of four buried water molecules thought to stabilize the filter, this region exhibits a significant degree of flexibility. In particular, pronounced distortion of filter occurs if no ions are present within it. The two most readily permeant ions, K(+) and Rb(+), are similar in their interactions with the selectivity filter. In contrast, Na(+) ions tend to distort the filter by binding to a ring of four carbonyl oxygens. The larger Cs(+) ions result in a small degree of expansion of the filter relative to the x-ray structure. Cs(+) ions also appear to interact differently with the gate region of the channel, showing some tendency to bind within a predominantly hydrophobic pocket. The four water molecules buried between the back of the selectivity filter and the remainder of the protein show comparable mobility to the surrounding protein and do not exchange with water molecules within the filter or the central cavity. A preliminary comparison of the use of particle mesh Ewald versus cutoff protocols for the treatment of long-range electrostatics suggests some difference in the kinetics of ion translocation within the filter.
Collapse
Affiliation(s)
- Carmen Domene
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | | |
Collapse
|
35
|
Arinaminpathy Y, Biggin PC, Shrivastava IH, Sansom MSP. A prokaryotic glutamate receptor: homology modelling and molecular dynamics simulations of GluR0. FEBS Lett 2003; 553:321-7. [PMID: 14572644 DOI: 10.1016/s0014-5793(03)01036-6] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
GluR0 is a prokaryotic homologue of mammalian glutamate receptors that forms glutamate-activated, potassium-selective ion channels. The topology of its transmembrane (TM) domain is similar to that of simple potassium channels such as KcsA. Two plausible alignments of the sequence of the TM domain of GluR0 with KcsA are possible, differing in the region of the P helix. We have constructed homology models based on both alignments and evaluated them using 6 ns duration molecular dynamics simulations in a membrane-mimetic environment. One model, in which an insertion in GluR0 relative to KcsA is located in the loop between the M1 and P helices, is preferred on the basis of lower structural drift and maintenance of the P helix conformation during simulation. This model also exhibits inter-subunit salt bridges that help to stabilise the TM domain tetramer. During the simulation, concerted K(+) ion-water movement along the selectivity filter is observed, as is the case in simulations of KcsA. K(+) ion exit from the central cavity is associated with opening of the hydrophobic gate formed by the C-termini of the M2 helices. In the intact receptor the opening of this gate will be controlled by interactions with the extramembranous ligand-binding domains.
Collapse
Affiliation(s)
- Yalini Arinaminpathy
- Department of Biochemistry, The University of Oxford, The Rex Richards Building, South Parks Road, OX1 3QU Oxford, UK
| | | | | | | |
Collapse
|
36
|
Rovere M, Gallo P. Effects of confinement on static and dynamical properties of water. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2003; 12:77-81. [PMID: 15007683 DOI: 10.1140/epje/i2003-10027-5] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
Molecular-dynamics results on water confined in a silica pore are reviewed and discussed in connection with experiments performed on water in Vycor and with studies of water in contact with proteins. The properties of confined water are studied as a function of both temperature and hydration level. The interaction of water in the film close to the substrate with the silica atoms induces a strong distortion of the hydrogen bond network. At high hydration levels a double dynamical regime is observed. At low hydration an anomalous diffusion is found upon supercooling with a transition from a Brownian to a non-Brownian regime on approaching the substrate in agreement with results found in studies of water in contact with globular proteins.
Collapse
Affiliation(s)
- M Rovere
- Dipartimento di Fisica E. Amaldi, Universitá Roma Tre, INFM Unitá Roma Tre, Italy.
| | | |
Collapse
|
37
|
Elmore DE, Dougherty DA. Investigating lipid composition effects on the mechanosensitive channel of large conductance (MscL) using molecular dynamics simulations. Biophys J 2003; 85:1512-24. [PMID: 12944269 PMCID: PMC1303328 DOI: 10.1016/s0006-3495(03)74584-6] [Citation(s) in RCA: 84] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Previous experimental work has shown that the functional properties of the mechanosensitive channel of large conductance (MscL) are affected by variations in lipid composition. Here, we utilize molecular dynamics simulations of Mycobacterium tuberculosis MscL to investigate such lipid composition effects on a molecular level. In particular, two sets of simulations were performed. In the first, trajectories using lipids with different headgroups (phosphatidylcholine and phosphatidylethanolamine) were compared. Protein-lipid interactions were clearly altered by the headgroup changes, leading to conformational differences in the C-terminal region of M. tuberculosis MscL. In the second set of simulations, lipid tails were gradually shortened, thinning the membrane over a molecular dynamics trajectory. These simulations showed evidence of hydrophobic matching between MscL and the lipid membrane, as previously proposed. For all simulations, protein-lipid interaction energies in the second transmembrane region were correlated to mutagenic data, emphasizing the importance of lipid interactions for proper MscL function.
Collapse
Affiliation(s)
- Donald E Elmore
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | | |
Collapse
|
38
|
Bond PJ, Sansom MSP. Membrane protein dynamics versus environment: simulations of OmpA in a micelle and in a bilayer. J Mol Biol 2003; 329:1035-53. [PMID: 12798692 DOI: 10.1016/s0022-2836(03)00408-x] [Citation(s) in RCA: 114] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The bacterial outer membrane protein OmpA is one of the few membrane proteins whose structure has been solved both by X-ray crystallography and by NMR. Crystals were obtained in the presence of detergent, and the NMR structure is of the protein in a detergent micelle. We have used 10 ns duration molecular dynamics simulations to compare the behaviour of OmpA in a detergent micelle and in a phospholipid bilayer. The dynamic fluctuations of the protein structure seem to be ca 1.5 times greater in the micelle environment than in the lipid bilayer. There are subtle differences between the nature of OmpA-detergent and OmpA-lipid interactions. As a consequence of the enhanced flexibility of the OmpA protein in the micellar environment, side-chain torsion angle changes are such as to lead to formation of a continuous pore through the centre of the OmpA molecule. This may explain the experimentally observed channel formation by OmpA.
Collapse
Affiliation(s)
- Peter J Bond
- Laboratory of Molecular Biophysics, Department of Biochemistry, The University of Oxford, The Rex Richards Building, South Parks Road, Oxford OX1 3QU, UK
| | | |
Collapse
|
39
|
Campbell JD, Biggin PC, Baaden M, Sansom MSP. Extending the structure of an ABC transporter to atomic resolution: modeling and simulation studies of MsbA. Biochemistry 2003; 42:3666-73. [PMID: 12667056 DOI: 10.1021/bi027337t] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular modeling and simulation approaches have been use to generate a complete model of the prokaryotic ABC transporter MsbA from Escherichia coli, starting from the low-resolution structure-based Calpha trace (PDB code 1JSQ). MsbA is of some biomedical interest as it is homologous to mammalian transporters such as P-glycoprotein and TAP. The quality of the MsbA model is assessed using a combination of molecular dynamics simulations and static structural analysis. These results suggest that the approach adopted for MsbA may be of general utility for generating all atom models from low-resolution crystal structures of membrane proteins. Molecular dynamics simulations of the MsbA model inserted in a fully solvated octane slab (a membrane mimetic environment) reveal that while the monomer is relatively stable, the dimer is unstable and undergoes significant conformational drift on a nanosecond time scale. This suggests that the MsbA crystal dimer may not correspond to the MsbA dimer in vivo. An alternative model of the dimer is discussed in the context of available experimental data.
Collapse
Affiliation(s)
- Jeff D Campbell
- Laboratory of Molecular Biophysics, Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, United Kingdom
| | | | | | | |
Collapse
|
40
|
Capener CE, Proks P, Ashcroft FM, Sansom MSP. Filter flexibility in a mammalian K channel: models and simulations of Kir6.2 mutants. Biophys J 2003; 84:2345-56. [PMID: 12668443 PMCID: PMC1302801 DOI: 10.1016/s0006-3495(03)75040-1] [Citation(s) in RCA: 52] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
The single-channel conductance varies significantly between different members of the inward rectifier (Kir) family of potassium channels. Mutations at three sites in Kir6.2 have been shown to produce channels with reduced single-channel conductance, the largest reduction (to 40% of wild-type) being for V127T. We have used homology modeling (based on a KcsA template) combined with molecular dynamics simulations in a phosphatidycholine bilayer to explore whether changes in structural dynamics of the filter were induced by three such mutations: V127T, M137C, and G135F. Overall, 12 simulations of Kir6.2 models, corresponding to a total simulation time of 27 ns, have been performed. In these simulations we focused on distortions of the selectivity filter, and on the presence/absence of water molecules lying behind the filter, which form interactions with the filter and the remainder of the protein. Relative to the wild-type simulation, the V127T mutant showed significant distortion of the filter such that approximately 50% of the simulation time was spent in a closed conformation. While in this conformation, translocation of K(+) ions between sites S1 and S2 was blocked. The distorted filter conformation resembles that of the bacterial channel KcsA when crystallized in the presence of a low [K(+)]. This suggests filter distortion may be a possible general model for determining the conductance of K channels.
Collapse
Affiliation(s)
- Charlotte E Capener
- Laboratory of Molecular Biophysics, Department of Biochemistry, Rex Richards Building, University of Oxford, South Parks Road, Oxford, OX1 3QU, United Kingdom
| | | | | | | |
Collapse
|
41
|
Domene C, Bond PJ, Sansom MS. Membrane protein simulations: ion channels and bacterial outer membrane proteins. ADVANCES IN PROTEIN CHEMISTRY 2003; 66:159-93. [PMID: 14631819 DOI: 10.1016/s0065-3233(03)66005-5] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Carmen Domene
- Laboratory of Molecular Biophysics (LMB), Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | | | | |
Collapse
|
42
|
MacCallum JL, Tieleman DP. Structures of neat and hydrated 1-octanol from computer simulations. J Am Chem Soc 2002; 124:15085-93. [PMID: 12475354 DOI: 10.1021/ja027422o] [Citation(s) in RCA: 95] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics computer simulations of 1-octanol and its mixtures with water have been performed. The liquid is composed of regions enriched in either hydrocarbons or hydroxyl groups. In neat octanol, the hydroxyl groups form clusters of long, thin chains. Upon the addition of water, the clusters become longer and more spherical, forming a structure that can be described as consisting of "overlapping elongated inverse micelles". The structures of the mixtures obtained at different hydration levels are consistent with those of experimental diffraction studies of water/octanol mixtures and previous computer simulations of neat and water-saturated octanol. The saturation point of the model has been calculated using the cavity-bias particle insertion method. The solubility of water in octanol is slightly too low compared to experimental results, and suggestions for possible improvements to the force field are made.
Collapse
Affiliation(s)
- Justin L MacCallum
- Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta, Canada T2N 1N4
| | | |
Collapse
|
43
|
Tieleman DP, Hess B, Sansom MSP. Analysis and evaluation of channel models: simulations of alamethicin. Biophys J 2002; 83:2393-407. [PMID: 12414676 PMCID: PMC1302328 DOI: 10.1016/s0006-3495(02)75253-3] [Citation(s) in RCA: 102] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Alamethicin is an antimicrobial peptide that forms stable channels with well-defined conductance levels. We have used extended molecular dynamics simulations of alamethicin bundles consisting of 4, 5, 6, 7, and 8 helices in a palmitoyl-oleolyl-phosphatidylcholine bilayer to evaluate and analyze channel models and to link the models to the experimentally measured conductance levels. Our results suggest that four helices do not form a stable water-filled channel and might not even form a stable intermediate. The lowest measurable conductance level is likely to correspond to the pentamer. At higher aggregation numbers the bundles become less symmetrical. Water properties inside the different-sized bundles are similar. The hexamer is the most stable model with a stability comparable with simulations based on crystal structures. The simulation was extended from 4 to 20 ns or several times the mean passage time of an ion. Essential dynamics analyses were used to test the hypothesis that correlated motions of the helical bundles account for high-frequency noise observed in open channel measurements. In a 20-ns simulation of a hexameric alamethicin bundle, the main motions are those of individual helices, not of the bundle as a whole. A detailed comparison of simulations using different methods to treat long-range electrostatic interactions (a twin range cutoff, Particle Mesh Ewald, and a twin range cutoff combined with a reaction field correction) shows that water orientation inside the alamethicin channels is sensitive to the algorithms used. In all cases, water ordering due to the protein structure is strong, although the exact profile changes somewhat. Adding an extra 4-nm layer of water only changes the water ordering slightly in the case of particle mesh Ewald, suggesting that periodicity artifacts for this system are not serious.
Collapse
Affiliation(s)
- D Peter Tieleman
- Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada.
| | | | | |
Collapse
|
44
|
Sansom MSP, Shrivastava IH, Bright JN, Tate J, Capener CE, Biggin PC. Potassium channels: structures, models, simulations. BIOCHIMICA ET BIOPHYSICA ACTA 2002; 1565:294-307. [PMID: 12409202 DOI: 10.1016/s0005-2736(02)00576-x] [Citation(s) in RCA: 92] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Potassium channels have been studied intensively in terms of the relationship between molecular structure and physiological function. They provide an opportunity to integrate structural and computational studies in order to arrive at an atomic resolution description of mechanism. We review recent progress in K channel structural studies, focussing on the bacterial channel KcsA. Structural studies can be extended via use of computational (i.e. molecular simulation) approaches in order to provide a perspective on aspects of channel function such as permeation, selectivity, block and gating. Results from molecular dynamics simulations are shown to be in good agreement with recent structural studies of KcsA in terms of the interactions of K(+) ions with binding sites within the selectivity filter of the channel, and in revealing the importance of filter flexibility in channel function. We discuss how the KcsA structure may be used as a template for developing structural models of other families of K channels. Progress in this area is explored via two examples: inward rectifier (Kir) and voltage-gated (Kv) potassium channels. A brief account of structural studies of ancillary domains and subunits of K channels is provided.
Collapse
Affiliation(s)
- Mark S P Sansom
- Laboratory of Molecular Biophysics, Department of Biochemistry, The University of Oxford, The Rex Richards Building, South Parks Road, Oxford, UK.
| | | | | | | | | | | |
Collapse
|
45
|
Shrivastava IH, Tieleman DP, Biggin PC, Sansom MSP. K(+) versus Na(+) ions in a K channel selectivity filter: a simulation study. Biophys J 2002; 83:633-45. [PMID: 12124253 PMCID: PMC1302175 DOI: 10.1016/s0006-3495(02)75197-7] [Citation(s) in RCA: 124] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022] Open
Abstract
Molecular dynamics simulations of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal significant differences in interactions of the selectivity filter with K(+) compared with Na(+) ions. K(+) ions and water molecules within the filter undergo concerted single-file motion in which they translocate between adjacent sites within the filter on a nanosecond timescale. In contrast, Na(+) ions remain bound to sites within the filter and do not exhibit translocation on a nanosecond timescale. Furthermore, entry of a K(+) ion into the filter from the extracellular mouth is observed, whereas this does not occur for a Na(+) ion. Whereas K(+) ions prefer to sit within a cage of eight oxygen atoms of the filter, Na(+) ions prefer to interact with a ring of four oxygen atoms plus two water molecules. These differences in interactions in the selectivity filter may contribute to the selectivity of KcsA for K(+) ions (in addition to the differences in dehydration energy between K(+) and Na(+)) and the block of KcsA by internal Na(+) ions. In our simulations the selectivity filter exhibits significant flexibility in response to changes in ion/protein interactions, with a somewhat greater distortion induced by Na(+) than by K(+) ions.
Collapse
Affiliation(s)
- Indira H Shrivastava
- Laboratory of Molecular Biophysics, Department of Biochemistry, University of Oxford, Oxford OX1 3QU, United Kingdom
| | | | | | | |
Collapse
|