We present a method of constructing the volume meshes of the membrane-channel protein system for finite element simulation of ion channels. The membrane channel system consists of the solvent region and the membrane-protein region. Our method focuses on labeling the tetrahedra in the solvent and membrane-protein regions and collecting the interface triangles between different regions. It contains two stages. Firstly, a volume mesh conforming the surface of the channel protein is generated by the surface and volume mesh generation tools: TMSmesh and TetGen. Then a walk-and-detect algorithm is used to identify the pore region to embed the membrane correctly. This method is shown to be robust because of its independence of the pore structure of the ion channels. In addition, we can also get the information of whether the ion channel is open or closed by the walk-and-detect algorithm. An on-line meshing procedure will be available at our website www.continuummodel.org.
2
Dostęp do pełnego tekstu na zewnętrznej witrynie WWW
Gramicidin A is a small and well characterized peptide that forms an ion channel in lipid membranes. An important feature of gramicidin A (gA) pore is that its conductance is affected by the electric charges near the its entrance. This property has led to the application of gramicidin A as a biochemical sensor for monitoring and quantifying a number of chemical and enzymatic reactions. Here, a mathematical model of conductance changes of gramicidin A pores in response to the presence of electrical charges near its entrance, either on membrane surface or attached to gramicidin A itself, is presented. In this numerical simulation, a two dimensional computational domain is set to mimic the structure of a gramicidin A channel in the bilayer surrounded by electrolyte. The transport of ions through the channel is modeled by the Poisson-Nernst-Planck (PNP) equations that are solved by Finite Element Method (FEM). Preliminary numerical simulations of this mathematical model are in qualitative agreement with the experimental results in the literature. In addition to the model and simulations, we also present the analysis of the stability of the solution to the boundary conditions and the convergence of FEM method for the two dimensional PNP equations in our model.
JavaScript jest wyłączony w Twojej przeglądarce internetowej. Włącz go, a następnie odśwież stronę, aby móc w pełni z niej korzystać.