Molecular dynamics simulation of six β-blocker drugs passing across POPC bilayer

<div><p>Six selected β-blocker drugs (alprenolol, atenolol, metoprolol, nadolol, pindolol and propranolol) passing across 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer were studied using all-atom molecular dynamics simulation. The free energy profiles can be divided into two groups, according to their shapes: the free energy curve of group one (atenolol, nadolol and pindolol) has an obvious minimum while that of the other group (propranolol, metoprolol and alprenolol) is flat inside membrane. Energy analysis shows that electrostatic interaction plays an important role for the first group drugs. The hydrogen bond analysis results also certify that the first group drugs form more hydrogen bonds than the other β-blockers. The calculated permeability sequence agrees with the experimental ones. Our calculation suggests that the permeability model using potential of mean force (PMF) method can be also applied to chemically similar compounds besides chemically diverse compounds.</p>