Chen Ada Y, Brooks Bernard R, Damjanovic Ana
Department of Physics & Astronomy, Johns Hopkins University, Baltimore, Maryland; Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, Maryland.
Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, Maryland.
Biophys J. 2021 Aug 3;120(15):3050-3069. doi: 10.1016/j.bpj.2021.06.013. Epub 2021 Jun 30.
Through molecular dynamics (MD) and free energy simulations in electric fields, we examine the factors influencing conductance of bacterial voltage-gated sodium channel NaMs. The channel utilizes four glutamic acid residues in the selectivity filter (SF). Previously, we have shown, through constant pH and free energy calculations of pKa values, that fully deprotonated, singly protonated, and doubly protonated states are all feasible at physiological pH, depending on how many ions are bound in the SF. With 173 MD simulations of 450 or 500 ns and additional free energy simulations, we determine that the conductance is highest for the deprotonated state and decreases with each additional proton bound. We also determine that the pKa value of the four glutamic residues for the transition between deprotonated and singly protonated states is close to the physiological pH and that there is a small voltage dependence. The pKa value and conductance trends are in agreement with experimental work on bacterial Na channels, which show a decrease in maximal conductance with lowering of pH, with pKa in the physiological range. We examine binding sites for Na in the SF, compare with previous work, and note a dependence on starting structures. We find that narrowing of the gate backbone to values lower than the crystal structure's backbone radius reduces the conductance, whereas increasing the gate radius further does not affect the conductance. Simulations with some amount of negatively charged lipids as opposed to purely neutral lipids increases the conductance, as do simulations at higher voltages.
通过电场中的分子动力学(MD)和自由能模拟,我们研究了影响细菌电压门控钠通道NaMs电导率的因素。该通道在选择性过滤器(SF)中利用四个谷氨酸残基。此前,我们通过对pKa值进行恒定pH和自由能计算表明,在生理pH条件下,完全去质子化、单质子化和双质子化状态都是可行的,这取决于SF中结合的离子数量。通过对450或500 ns的173次MD模拟以及额外的自由能模拟,我们确定去质子化状态的电导率最高,并且随着每个额外质子的结合而降低。我们还确定,四个谷氨酸残基从去质子化状态转变为单质子化状态的pKa值接近生理pH,并且存在较小的电压依赖性。pKa值和电导率趋势与细菌钠通道的实验结果一致,实验表明随着pH降低,最大电导率降低,pKa值在生理范围内。我们研究了SF中Na的结合位点,与之前的工作进行比较,并注意到对起始结构的依赖性。我们发现,将通道主干变窄到低于晶体结构主干半径的值会降低电导率,而进一步增加通道半径不会影响电导率。与纯中性脂质相比,含有一定量带负电荷脂质的模拟会增加电导率,在更高电压下的模拟也是如此。