Anesthetics mechanism on a DMPC lipid membrane model: Insights from molecular dynamics simulations.

To provide insight into the molecular mechanisms of local anesthetic action, we have carried out an extensive investigation of two amide type local anesthetics, lidocaine and articaine in both charged and uncharged forms, interacting with DMPC lipid membrane. We have applied both standard molecular dynamics simulations and metadynamics simulations to provide a detailed description of the free energy landscape of anesthetics embedded in the lipid bilayer. The global minimum of the free energy surface (equilibrium position of anesthetics in the lipid membrane) occurred around 1nm of the bilayer center. The uncharged anesthetics show more affinity to bind to this region compared to the charged drugs. The binding free energy of uncharged lidocaine in the membrane (-30.3kJ/mol) is higher than uncharged articaine (-24.0kJ/mol), which is in good agreement with higher lipid solubility of lidocaine relative to the articaine. The octanol/water partition coefficient of uncharged drugs was also investigated using expanded ensemble simulations. In addition, complementary standard MD simulations were carried out to study the partitioning behavior of multiple anesthetics inside the lipid bilayer. The results obtained here are in line with previously reported simulations and suggest that the different forms of anesthetics induce different structural modifications in the lipid bilayer, which can provide new insights into their complex membrane translocation phenomena.

Study of procaine and tetracaine in the lipid bilayer using molecular dynamics simulation.

Despite available experimental results, the molecular mechanism of action of local anesthetics upon the nervous system and contribution of the cell membrane to the process are still controversial. In this work, molecular dynamics simulations were performed to investigate the effect of two clinically used local anesthetics, procaine and tetracaine, on the structure and dynamics of a fully hydrated dimyristoylphosphatidylcholine lipid bilayer. We focused on comparing the main effects of uncharged and charged drugs on various properties of the lipid membrane: mass density distribution, diffusion coefficient, order parameter, radial distribution function, hydrogen bonding, electrostatic potential, headgroup angle, and water dipole orientation. To compare the diffusive nature of anesthetic through the lipid membrane quantitatively, we investigated the hexadecane/water partition coefficient using expanded ensemble simulation. We predicted the permeability coefficient of anesthetics in the following order: uncharged tetracaine > uncharged procaine > charged tetracaine > charged procaine. We also shown that the charged forms of drugs are more potent in hydrogen bonding, disturbing the lipid headgroups, changing the orientation of water dipoles, and increasing the headgroup electrostatic potential more than uncharged drugs, while the uncharged drugs make the lipid diffusion faster and increase the tail order parameter. The results of these simulation studies suggest that the different forms of anesthetics induce different structural modifications in the lipid bilayer, which provides new insights into their molecular mechanism.

Study of curcumin behavior in two different lipid bilayer models of liposomal curcumin using molecular dynamics simulation.

Liposomal formulation of curcumin is an important therapeutic agent for the treatment of various cancers. Despite extensive studies on the biological effects of this formulation in cancer treatment, much remains unknown about curcumin-liposome interactions. Understanding how different lipid bilayers respond to curcumin molecule may help us to design more effective liposomal curcumin. Here, we used molecular dynamics simulation method to investigate the behavior of curcumin in two lipid bilayers commonly used in preparation of liposomal curcumin, namely dipalmitoylphosphatidylcholine (DPPC) and dimyristoylphosphatidylglycerol (DMPG). First, the free energy barriers for translocation of one curcumin molecule from water to the lipid bilayer were determined by using the potential of mean force (PMF). The computed free energy profile exhibits a global minimum at the solvent-headgroup interface (LH region) for both lipid membranes. We also evaluated the free energy difference between the equilibrium position of curcumin in the lipid bilayer and bulk water as the excess chemical potential. Our results show that curcumin has the higher affinity in DMPG compared to DPPC lipid bilayer (-8.39 vs. -1.69 kBT) and this is related to more hydrogen bond possibility for curcumin in DMPG lipid membrane. Next, using an unconstrained molecular dynamic simulation with curcumin initially positioned at the center of lipid bilayer, we studied various properties of each lipid bilayer system in the presence of curcumin molecule that was in full agreement with PMF and experimental data. The results of these simulation studies suggest that membrane composition could have a large effect on interaction of curcumin-lipid bilayer.