Abstract:We propose a theoretical approach to estimate the permeability coefficient 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 membranes. 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 \pm 0.01 ~\mathrm{cm~s^{-1}}$) than for methylamine ($0.069\pm 0.006~\mathrm{cm~s^{-1}}$), while the values of the permeability coefficient are satisfactorily close to those obtained from the brute-force MD simulations [$0.18\pm 0.03 ~\mathrm{cm~s^{-1}}$ and $0.052 \pm 0.005 ~\mathrm{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.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to quantify the permeability coefficient of substances through cell membranes by molecular dynamics (MD) simulations. Specifically, the author proposes a new method based on the Returning Probability (RP) theory, which can be used to estimate the permeability coefficient of substances passing through the membrane from the donor phase (D) to the receptor phase (A). The core of this method lies in reformulating the RP theory to describe the permeation phenomenon, thus theoretically providing a new framework for analyzing the process of substances passing through the membrane. ### Main contributions of the paper 1. **Theoretical framework**: By applying the RP theory to the membrane permeation process, the author establishes the relationship between the permeability coefficient and the thermodynamic and kinetic properties of substances in the membrane. This enables the calculation of these properties through MD simulations and further the estimation of the permeability coefficient. 2. **Application verification**: Ethanol and methylamine are used as research objects in the paper, and verification is carried out under the conditions of infinite dilution and 1 mol% concentration respectively. The results show that the permeability coefficients obtained by the new method are consistent with the traditional MD simulation results, proving the effectiveness of the method. 3. **Mechanism analysis**: Through the analysis of thermodynamic and kinetic contributions, the author further clarifies the differences in the permeability of ethanol and methylamine, especially that the sensitivity of the free - energy barrier of ethanol in the membrane to concentration leads to a higher concentration - dependence of its permeability coefficient. ### Formula analysis - **Definition of permeability coefficient**: \[ J_{\text{ss}} = c_D P_{\text{ss}} \] where $ J_{\text{ss}} $ is the flux at the steady state, $ c_D $ is the concentration in the donor phase, and $ P_{\text{ss}} $ is the permeability coefficient. - **Key equation in RP theory**: \[ s \hat{P}(s) = \frac{1}{\sigma} \int dr \int d\Lambda S(z, \Lambda) g_{\text{eq}}(r, \Lambda) \left(1 + k_{R \to A} \int_0^\infty dt P_{\text{RET}}(t)\right)^{-1} \] where $ \hat{P}(s) $ is the Laplace transform of the permeability coefficient, $ S(z, \Lambda) $ is the reaction sedimentation function, $ g_{\text{eq}}(r, \Lambda) $ is the equilibrium distribution function, and $ P_{\text{RET}}(t) $ is the return probability. - **Calculation of free - energy difference**: \[ \Delta G(z) = \Delta \mu_R(z) - \Delta \mu_D \] where $ \Delta \mu_R(z) $ and $ \Delta \mu_D $ are the solvation free energies of substances at the position $ z $ in the membrane and in the donor phase respectively. ### Conclusion By introducing the RP theory, this paper provides a new perspective for the theoretical analysis of the membrane permeation process. Through MD simulations, the author successfully estimates the permeability coefficients of ethanol and methylamine and reveals their permeation mechanisms under different concentration conditions. This method not only provides theoretical support for understanding the process of substances passing through the membrane but also provides valuable references for practical applications such as drug delivery.

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

Molecular dynamics simulation on permeation of acetone/nitrogen mixed gas through Al2O3 microporous membranes

[Development of Membrane Permeability Coefficient by Means of Novel Molecular Dynamics Methods]

Study on Permeation of Acetone/Nitrogen Mixed Gas Through Al2o3 Microporous Membranes: Molecular Dynamics

Molecular Transport Through Membranes: Accurate Permeability Coefficients from Multidimensional Potentials of Mean Force and Local Diffusion Constants

Statistical Mechanical Theories of Membrane Permeability

Calculation of Permeability Coefficients from Solute Equilibration Dynamics: An Assessment of Various Methods

Permeability of a Fluid Lipid Bilayer to Short-Chain Alcohols from First Principles

Mechanistic Insights into Passive Membrane Permeability of Drug-like Molecules from a Weighted Ensemble of Trajectories

Determining small-molecule permeation through lipid membranes

Assessing the Intestinal Permeability of Small Molecule Drugs Via Diffusion Motion on a Multidimensional Free Energy Surface.

Impact of the Unstirred Water Layer on the Permeation of Small-Molecule Drugs.

Free Energy Calculations of Membrane Permeation: Challenges due to Strong Headgroup-Solute Interactions

Data-driven equation for drug-membrane permeability across drugs and membranes

A machine-learning-assisted study of the permeability of small drug-like molecules across lipid membranes

Calculating Single-Channel Permeability and Conductance from Transition Paths.

Non-equilibrium Molecular Dynamics Simulation on Pure Gas Permeability Through Carbon Membranes

Effect of Molecular Association on Solubility, Diffusion, and Permeability in Polymeric Membranes

A mass transfer model for pure alcoholic permeation through the PDMS membrane

How Usable Are Published Permeability Data?

Interpreting effective energy barriers to membrane permeation in terms of a heterogeneous energy landscape