Membrane proteins are involved in a large variety of functions. Most of these protein functions are regulated by ligand binding with diverse modes of action: agonists, partial agonists, antagonists, and allosteric modulators, potentiators and inhibitors. From the pharmacological point of view, membrane proteins are one if not the major target for drug development. However, experimental structure determination of membrane proteins in complex or in free form still represents a great challenge. Molecular dynamics (MD) simulations commonly reach the microsecond scale on membrane systems. This numerical tool is mature enough to predict and add molecular details on the different ligand-binding modes. In the present chapter, I will present the different steps to design, simulate, and analyze a MD simulation system containing a protein embedded in a membrane and surrounded by water and ligand. As an illustration, the simulation of the ligand-gated ion channel γ-aminobutyric acid type A receptor (GABAAR) surrounded by one of its allosteric potentiators, bromoform, will be presented and discussed.
Keywords: Bromoform; Cys-loop receptors; GABAAR; Ligand binding; Membrane dynamic; Membrane protein; Molecular dynamics; System preparation.