Calculating Kinetic Rates and Membrane Permeability from Biased Simulation

Magd Badaoui, Adam Kells, Carla Molteni, Callum J. Dickson, Viktor Hornak, Edina Rosta

Research output: Contribution to journalArticlepeer-review

20 Citations (Scopus)
284 Downloads (Pure)


We present a simple approach to calculate the kinetic properties of lipid membrane crossing processes from biased molecular dynamics simulations. We demonstrate that by using biased simulations, one can obtain highly accurate kinetic information with significantly reduced computational time with respect to unbiased simulations. We describe how to conveniently calculate the transition rates to enter, cross and exit the membrane in terms of mean first passage times. To obtain free energy barriers and relaxation times from biased simulations only, we constructed Markov models using the Dynamic Histogram Analysis Method (DHAM). The permeability coefficients that are calculated from the relaxation times are found to correlate highly with experimentally evaluated values. We show that more generally, certain calculated kinetic properties linked to the cross-ing of the membrane layer (e.g., barrier height and barrier crossing rates) are good indicators of ordering drugs by permeabil-ity. Extending the analysis to a 2D Markov model provides a physical description of the membrane crossing mechanism.
Original languageEnglish
Publication statusE-pub ahead of print - 24 Sept 2018


Dive into the research topics of 'Calculating Kinetic Rates and Membrane Permeability from Biased Simulation'. Together they form a unique fingerprint.

Cite this