Molecular dynamics of a vasopressin V2 receptor in a phospholipid bilayer membrane

Molecular dynamics simulations were carried out for a V2 receptor (V2R) model embedded in a dimyristoylphosphatidylcholine (DMPC) bilayer. Both free and ligand-bound states of V2R were modeled. Our initial V2R model was obtained using a rule-based automated method for GPCR modeling and refined using constrained simulated annealing in vacuo. The docking site of the native vasopressin ligand was selected and justified upon consideration of ligand-receptor interactions and structure-activity data. The primary purpose of this work was to investigate the usefulness of MD simulation of an integral membrane protein like a GPCR receptor, upon inclusion of a carefully parameterized surrounding lipid membrane and water. Physical properties of the system were evaluated and compared with the fully hydrated pure DMPC bilayer membrane. The solvation interactions, individual lipid-protein interaction and fluctuations of the protein, the lipid, and water were analyzed in detail. As expected, the membrane-spanning helices of the protein fluctuate less than the peripheral loops do. The protein appears to disturb the local lipid structure. Simulations were carried out using AMBER 4.1 package upon constant number-pressure-temperature (NPT) conditions on massively parallel computers Cray T3E and IBM SP2.