Simulation study of a rectifying bipolar ion channel: Detailed model versus reduced model

Abstract

We study a rectifying mutant of the OmpF porin ion channel using both all-atom and reduced models. The mutant was created by Miedema et al. [Nano Lett., 2007, 7, 2886] on the basis of the NP semiconductor diode, in which an NP junction is formed. The mutant contains a pore region with positive amino acids on the left-hand side and negative amino acids on the right-hand side. Experiments show that this mutant rectifies. Although we do not know the structure of this mutant, we can build an all-atom model for it on the basis of the structure of the wild type channel. Interestingly, molecular dynamics simulations for this all-atom model do not produce rectification. A reduced model that contains only the important degrees of freedom (the positive and negative amino acids and free ions in an implicit solvent), on the other hand, exhibits rectification. Our calculations for the reduced model (using the Nernst-Planck equation coupled to Local Equilibrium Monte Carlo simulations) reveal a rectification mechanism that is different from that seen for semiconductor diodes. The basic reason is that the ions are different in nature from electrons and holes (they do not recombine). We provide explanations for the failure of the all-atom model including the effect of all the other atoms in the system as a noise that inhibits the response of ions (that would be necessary for rectification) to the polarizing external field.

Ми вивчаємо випростовуючий мутант iонного каналу OmpF порiну, використовуючи повну атомiстичну i спрощену моделi. Даний мутант був створений Мiедемою та iн. [Nano Lett., 2007, 7, 2886] на основi напiвпровiдникового дiода, в якому сформувався N-P-перехiд. Мутант мiстить пористу зону з позитивно зарядженими амiнокислотами злiва вiд бiксторону та негативно зарядженi амiнокислоти справа. Дослiди показали, що цей мутант має випростовуючi властивостi. Хоча структура цього мутанта невiдома, можна побудувати його повнiстю атомарну модель. Моделювання молекулярної динамiки для цiєї атомарної моделi не забезпечує ефекту випростовування. Водночас, спрощена модель, яка мiстить лише важливi ступенi вiльностi (додатнi та вiд’ємнi амiнокислоти i вiльнi iони у неявному розчиннику), забезпечує випростовування. Дослiдження, виконанi для спрощеної моделi (з використанням рiвняння Нернста-Планка у поєднаннi з моделюванням Монте Карло iз локальною рiвновагою), показали механiзм випростовування, який суттєво вiдрiзняється вiд напiвпровiдникових дiодiв. Головна причина полягає в тому, що iони за своєю природою вiдрiзняються вiд електронiв та дiрок (iони не рекомбiнуються). Ми пояснюємо неспроможнiсть повної атомарної моделi, включно з ефектом решти атомiв як шуму, який блокує вiдгук iонiв на зовнiшнє поляризуюче поле необхiдний для появи ефекту випростовування.