A methodology of quantifying membrane permeability based on returning probability theory and molecular dynamics simulation.

Autor:	Matsubara Y; Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan., Okabe R; Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan., Masayama R; Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan., Watanabe NM; Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan., Umakoshi H; Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan., Kasahara K; Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan., Matubayasi N; Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan.
Jazyk:	angličtina
Zdroj:	The Journal of chemical physics [J Chem Phys] 2024 Jul 14; Vol. 161 (2).
DOI:	10.1063/5.0214401
Abstrakt:	We propose a theoretical approach to estimate the permeability coefficients of substrates (permeants) for crossing membranes from donor (D) phase to acceptor (A) phase by means of molecular dynamics (MD) simulation. A fundamental aspect of our approach involves reformulating the returning probability (RP) theory, a rigorous bimolecular reaction theory, to describe permeation phenomena. This reformulation relies on the parallelism between permeation and bimolecular reaction processes. In the present method, the permeability coefficient is represented in terms of the thermodynamic and kinetic quantities for the reactive (R) phase that exists within the inner region of a membrane. One can evaluate these quantities using multiple MD trajectories starting from phase R. We apply the RP theory to the permeation of ethanol and methylamine at different concentrations (infinitely dilute and 1 mol % conditions of permeants). Under the 1 mol% condition, the present method yields a larger permeability coefficient for ethanol (0.12 ± 0.01 cm s-1) than for methylamine (0.069 ± 0.006 cm s-1), while the values of the permeability coefficient are satisfactorily close to those obtained from the brute-force MD simulations (0.18 ± 0.03 and 0.052 ± 0.005 cm s-1 for ethanol and methylamine, respectively). Moreover, upon analyzing the thermodynamic and kinetic contributions to the permeability, we clarify that a higher concentration dependency of permeability for ethanol, as compared to methylamine, arises from the sensitive nature of ethanol's free-energy barrier within the inner region of the membrane against ethanol concentration. (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)
Databáze:	MEDLINE
Externí odkaz:	Zobrazit plný text záznamu