Molecular dynamics study of free energy profiles for organic cations in gramicidin A channels

被引:18
作者
Hao, Y [1 ]
Pear, MR [1 ]
Busath, DD [1 ]
机构
[1] BROWN UNIV, DEPT PHYSIOL, PROVIDENCE, RI 02912 USA
关键词
D O I
10.1016/S0006-3495(97)78202-X
中图分类号
Q6 [生物物理学];
学科分类号
071011 ;
摘要
The free energy profiles for four organic cations in right-handed single-helix gramicidin A dimers were computed by using umbrella sampling molecular dynamics with CHARMM. ion-water column translocations were facilitated by using a novel ''water-tunnel'' approach. The overlapping pieces of free energy profile for adjacent windows were selected from three trajectories that differed in initial ion rotation and were aligned by the method of umbrella potential differences. Neglected long-range electrostatic energies from the bulk water and the bilayer were computed with DelPhi and added to the profile. The approach was corroborated for the formamidinium-guanidinium pair by using perturbation dynamics at axial positions 0, 6, 12, and 15 Angstrom from the channel center. The barrier to ethylammonium entry was prohibitive at 21 kcal/mol, whereas for methylammonium it was 5.5 kcal/mol, and the profile was quite flat through the channel, roughly consistent with conductance measurements. The profile for formamidinium was very similar to that of methylammonium. Guanidinium had a high entry barrier (Delta F = +8.6 kcal/mol) and a narrow deep central well (Delta F = -2.6 kcal/mol), qualitatively consistent with predictions from voltage-dependent potassium current block measurements. Its deep central well, contrasting with the flat profile for formamidinium, was verified with perturbation dynamics and was correlated with its high propensity to form hydrogen bonds with the channel at the dimer junction (not shared by the other three cations). Analysis of the ensemble average radial forces on the ions demonstrates that all four ions undergo compressive forces in the channel that are at maximum at the center of the monomer and relieved at the dimer junction, illustrating increased flexibility of the channel walls in the center of the channel.
引用
收藏
页码:1699 / 1716
页数:18
相关论文
共 56 条
[1]  
[Anonymous], REC CHEM PROG
[2]   ENERGETICS OF ION PERMEATION THROUGH MEMBRANE CHANNELS - SOLVATION OF NA+ BY GRAMICIDIN-A [J].
AQVIST, J ;
WARSHEL, A .
BIOPHYSICAL JOURNAL, 1989, 56 (01) :171-182
[3]  
ARSENIEV AS, 1986, BIOL MEMBRANY, V3, P1077
[4]   TEMPERATURE-DEPENDENT PROPERTIES OF GRAMICIDIN A CHANNELS [J].
BAMBERG, E ;
LAUGER, P .
BIOCHIMICA ET BIOPHYSICA ACTA, 1974, 367 (02) :127-133
[5]   AMINO-ACID SUBSTITUTIONS AND ION CHANNEL FUNCTION - MODEL-DEPENDENT CONCLUSIONS [J].
BECKER, MD ;
KOEPPE, RE ;
ANDERSEN, OS .
BIOPHYSICAL JOURNAL, 1992, 62 (01) :25-27
[6]  
BOGUSZ SJ, 1995, THESIS BROWN U PROVI
[7]   CHARMM - A PROGRAM FOR MACROMOLECULAR ENERGY, MINIMIZATION, AND DYNAMICS CALCULATIONS [J].
BROOKS, BR ;
BRUCCOLERI, RE ;
OLAFSON, BD ;
STATES, DJ ;
SWAMINATHAN, S ;
KARPLUS, M .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1983, 4 (02) :187-217
[8]  
BUSATH D, 1988, JERUS SYM Q, V21, P187
[9]  
BUSATH DD, 1993, ANNU REV PHYSIOL, V55, P473
[10]   A MOLECULAR-DYNAMICS STUDY OF GATING IN DIOXOLANE-LINKED GRAMICIDIN-A CHANNELS [J].
CROUZY, S ;
WOOLF, TB ;
ROUX, B .
BIOPHYSICAL JOURNAL, 1994, 67 (04) :1370-1386