Efflux dynamics of the antiseizure drug, levetiracetam, through the P-glycoprotein channel revealed by advanced comparative molecular simulations

Understanding the precise mechanistic details of the possible binding and transport of antiseizure medications (ASMs) through the P-glycoprotein (P-gp) efflux pump is essential to find strategies for the treatment of patients with epilepsy resistant to ASMs. In the present work, conventional molecular dynamics, binding free energy calculations, steered molecular dynamics and umbrella sampling were applied to study the interactions of levetiracetam and brivaracetam with P-gp and their possible egress path from the binding site. Comparative results for the control drugs, zosuquidar and verapamil, confirmed their established P-gp inhibitory activity. Brivaracetam, a non-substrate of P-gp, demonstrated stronger static and dynamic interactions with the exporter protein, than levetiracetam. The potential of mean force calculations indicated that the energy barriers through the ligand export were the lowest for levetiracetam, suggesting the drug as a P-gp substrate with facile passage through the transporter channel. Our findings also stressed the contribution of nonpolar interactions with P-gp channel lining as well as with membrane lipid molecules to hamper the ASM efflux by the transmembrane exporter. Appropriate structural engineering of the ASMs is thus recommended to address drug-resistant epilepsy.