Molecular dynamics simulation of anesthetic-phospholipid bilayer interactions.

To probe the hypothesis of a lipid-mediated mechanism of general anesthetic action on a molecular level, and to help elucidate the nature of the interactions of bioactive compounds with membranes, the effects of trichloroethylene (TCE), an inhalational general anesthetic, on a dioleoylphosphatidylcholine (DOPC) lipid bilayer have been investigated by molecular dynamics (MD) simulations at 37 degrees C and 1 atm and the results compared with 31P and 2H NMR experimental studies (Ref 1). The model used included a single TCE molecule embedded in a lipid bilayer consisting of 24 DOPC molecules and an 8 A layer of explicit water of solvation in each polar head group region of the bilayer, together with constant-pressure periodic boundary conditions in three dimensions. A comparison of the bilayer properties calculated in the presence and absence of the anesthetic led to the detection of three major perturbations of the bilayer caused by the anesthetic at 1 atm: i) an increase in the ratio of the effective areas of hydrocarbon tails and the head group per lipid, predicting the tendency of lipids near the anesthetic site of action to form a hexagonal phase (HII); ii) a slight increase in the frequency of chain dihedral angles found in the gauche conformation; and iii) a significant increase in the lateral mean-square displacement of lipid molecules, an indication of increased lipid lateral diffusion and membrane fluidity. The pressure antagonism of these effects was also studied by MD simulations at pressures of 200 and 400 atm. The study of the pressure reversibility of these effects at 200 and 400 atm indicated that they were partially prevented at 200 atm and essentially blocked at 400 atm, suggesting their probable relevance to the pressure reversal effect seen with general anesthesia. These results may thus provide insights into the interaction between general anesthetics and similar small organic molecules with membranes.