Computational Studies of Gramicidin Permeation: an Entry Way Sulfonate Enhances Cation Occupancy at Entry Sites

Overview

Journal Biochim Biophys Acta

Specialties Biochemistry
Biophysics

Date 2009 Apr 14

PMID 19361485

Authors

Morad Mustafa

Douglas J Henderson

David D Busath

Affiliations

Soon will be listed here.

Abstract

The impact on the cation-transport free-energy profile of replacing the C-terminal ethanolamine in the gramicidin A channel with a taurine residue is studied using molecular dynamics simulations of gramicidin A (1JNO) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl saline solution. The potential of mean force for ion transport is obtained by umbrella sampling. The presence of a negatively charged sulfonate group at the entrance of the gramicidin channel affects the depth and the location of the binding sites, producing a strong attraction for the cations in the bulk. The potential of mean force by the sulfonate acting directly through electrostatics and van der Waals interactions on the test ion is highly modulated by indirect effects (i.e., sulfonate effects on other components of the system that, in turn, affect the ion free-energy profile in the channel). Because the "entry" sites are located symmetrically at both entry and exit of the channel, the deeper free-energy wells should inhibit exit. Given that the channel has increased conductance experimentally, the simulation results suggest that the channel conductance is normally entry limited.

References

Capone R, Blake S, Restrepo M, Yang J, Mayer M . Designing nanosensors based on charged derivatives of gramicidin A. J Am Chem Soc. 2007; 129(31):9737-45. DOI: 10.1021/ja0711819. View

Roeske R, Pottorf R, Bridal T, Jin X, Busath D . Synthesis and channel properties of [Tau 16]gramicidin A. Biochim Biophys Acta. 1989; 982(2):223-7. DOI: 10.1016/0005-2736(89)90058-8. View

Allen T, Andersen O, Roux B . Molecular dynamics - potential of mean force calculations as a tool for understanding ion permeation and selectivity in narrow channels. Biophys Chem. 2006; 124(3):251-67. DOI: 10.1016/j.bpc.2006.04.015. View

Sansom M, Kerr I . Influenza virus M2 protein: a molecular modelling study of the ion channel. Protein Eng. 1993; 6(1):65-74. DOI: 10.1093/protein/6.1.65. View

Urry D, GOODALL M, Glickson J, Mayers D . The gramicidin A transmembrane channel: characteristics of head-to-head dimerized (L,D) helices. Proc Natl Acad Sci U S A. 1971; 68(8):1907-11. PMC: 389319. DOI: 10.1073/pnas.68.8.1907. View

Bastug T, Kuyucak S . Free energy simulations of single and double ion occupancy in gramicidin A. J Chem Phys. 2007; 126(10):105103. DOI: 10.1063/1.2710267. View

Wu Y, Voth G . A computational study of the closed and open states of the influenza a M2 proton channel. Biophys J. 2005; 89(4):2402-11. PMC: 1366740. DOI: 10.1529/biophysj.105.066647. View

Errington J, Kofke D . Calculation of surface tension via area sampling. J Chem Phys. 2007; 127(17):174709. DOI: 10.1063/1.2795698. View

Bastug T, Kuyucak S . Energetics of ion permeation, rejection, binding, and block in gramicidin A from free energy simulations. Biophys J. 2006; 90(11):3941-50. PMC: 1459526. DOI: 10.1529/biophysj.105.074633. View

10.

Boda D, Valisko M, Eisenberg B, Nonner W, Henderson D, Gillespie D . The effect of protein dielectric coefficient on the ionic selectivity of a calcium channel. J Chem Phys. 2006; 125(3):34901. DOI: 10.1063/1.2212423. View

11.

Chung S, Corry B . Conduction properties of KcsA measured using brownian dynamics with flexible carbonyl groups in the selectivity filter. Biophys J. 2007; 93(1):44-53. PMC: 1914447. DOI: 10.1529/biophysj.106.098954. View

12.

Kucerka N, Liu Y, Chu N, Petrache H, Tristram-Nagle S, Nagle J . Structure of fully hydrated fluid phase DMPC and DLPC lipid bilayers using X-ray scattering from oriented multilamellar arrays and from unilamellar vesicles. Biophys J. 2005; 88(4):2626-37. PMC: 1305359. DOI: 10.1529/biophysj.104.056606. View

13.

Bastug T, Patra S, Kuyucak S . Molecular dynamics simulations of gramicidin A in a lipid bilayer: from structure-function relations to force fields. Chem Phys Lipids. 2006; 141(1-2):197-204. DOI: 10.1016/j.chemphyslip.2006.02.012. View

14.

Mustafa M, Busath D . The gramicidin channel ion permeation free-energy profile: direct and indirect effects of CHARMM force field improvements. Interdiscip Sci. 2010; 1(2):113-27. PMC: 2806686. DOI: 10.1007/s12539-009-0025-3. View

15.

Harroun T, Heller W, Weiss T, Yang L, Huang H . Experimental evidence for hydrophobic matching and membrane-mediated interactions in lipid bilayers containing gramicidin. Biophys J. 1999; 76(2):937-45. PMC: 1300095. DOI: 10.1016/S0006-3495(99)77257-7. View

16.

Trzesniak D, Kunz A, van Gunsteren W . A comparison of methods to compute the potential of mean force. Chemphyschem. 2006; 8(1):162-9. DOI: 10.1002/cphc.200600527. View

17.

SARGES R, WITKOP B . GRAMICIDIN A. V. THE STRUCTURE OF VALINE- AND ISOLEUCINE-GRAMICIDIN A. J Am Chem Soc. 1965; 87:2011-20. DOI: 10.1021/ja01087a027. View

18.

Humphrey W, Dalke A, Schulten K . VMD: visual molecular dynamics. J Mol Graph. 1996; 14(1):33-8, 27-8. DOI: 10.1016/0263-7855(96)00018-5. View

19.

MacKerell Jr A, Feig M, Brooks 3rd C . Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations. J Comput Chem. 2004; 25(11):1400-15. DOI: 10.1002/jcc.20065. View

20.

Nonner W, Catacuzzeno L, Eisenberg B . Binding and selectivity in L-type calcium channels: a mean spherical approximation. Biophys J. 2000; 79(4):1976-92. PMC: 1301088. DOI: 10.1016/S0006-3495(00)76446-0. View