Molecular dynamics simulations of CBD and THC with a DPPC / DPPG monolayer membrane: Comparative study with morphine and lidocaine


  • Gonzalo Rafael Vargas Castañeda Universidad Michoacana de San Nicolás de Hidalgo
  • J. A. González Universidad Michoacana de San Nicolás de Hidalgo
  • G. Espinosa Universidad Michoacana de San Nicolás de Hidalgo



Molecular dynamics, CBD, THC, Monolayer membrane


In recent years, the discussion about using marijuana to treat some specific diseases has been the reason for several studies to determine how Cannabis could affect or benefit the human body, such as chronic pain. However, the legal status of marijuana in different parts of the world remains a problem when conducting experiments for academic purposes. In this way, algorithms based on molecular dynamics are of great help to study the change in the physical properties of membranes, which is our first defense barrier, in the presence of external agents. Here we used equilibrium molecular dynamics to study the interactions of a DPPC/DPPG monolayer membrane with CBD and THC, these two molecules being one of the main components of Cannabis, and compare the interactions with morphine and lidocaine which are drugs with well-known anesthetic properties. We also used the umbrella sampling method to determine the variations of the potential of mean forces as every single drug crosses the membrane obtaining an approximation of the hydration free energy, we obtained a value of 57 kJ mol-1 for CBD and 5 kJ mol-1 for THC; values reported for the first time using molecular dynamics.


W. B. O’Shaughnessy, On the preparations of the Indian hemp, or Gunjah: Cannabis indica their effects on the animal system in health, and their utility in the treatment of tetanus and other convulsive diseases. Provincial Medical Journal and Retrospect of the Medical Sciences 5 (1843) 123.

V. Di Marzo, A brief history of cannabinoid and endocannabinoid pharmacology as inspired by the work of British scientists. Trends in pharmacological sciences 27.3 (2006) 134.

M. Quiroga, Cannabis: efectos nocivos sobre la salud f´ısica, Adicciones 12 (2000) 117

H. V. Curran, et al., Which biological and self-report measures of cannabis use predict cannabis dependency and acute psychotic-like effects?, Psychological medicine 49 (2019) 1574

H. Meng et al., Cannabis and cannabinoids in cancer pain management, Current Opinion in Supportive and Palliative Care 14 (2020) 87

J. Ahrens et al., The nonpsychotropic cannabinoid cannabidiol modulates and directly activates alpha-1 and alpha-1-Beta glycine receptor function, Pharmacology 83 (2009) 217

D. G. Boychuk et al., The effectiveness of cannabinoids in the management of chronic nonmalignant neuropathic pain: a systematic review, Joral facial pain headache 29 (2015) 7

A. Aran et al., Brief report: cannabidiol-rich cannabis in children with autism spectrum disorder and severe behavioral problems-a retrospective feasibility study, Journal of Autism and Developmental Disorders 49 (2019) 1284

C. Carroll et al., Cannabis for dyskinesia in Parkinson disease: a randomized double-blind crossover study, Neurology 63 (2004) 1245

M. W. Elmes et al., Fatty acid-binding proteins (FABPs) are intracellular carriers for ∆9-tetrahydrocannabinol (THC) and cannabidiol (CBD), Journal of Biological Chemistry 290 (2015) 8711

P. Morales, D. P. Hurst, and P. H. Reggio, Molecular targets of the phytocannabinoids: a complex picture, Phytocannabinoids 2 (2017) 103

Y.-m. Bian et al., Computational systems pharmacology analysis of cannabidiol: a combination of chemogenomicsknowledgebase network analysis and integrated in silico modeling and simulation, Acta Pharmacologica Sinica 40 (2019) 374

W. Kopec, J. Telenius, and H. Khandelia, Molecular dynam- ´ ics simulations of the interactions of medicinal plant extracts and drugs with lipid bilayer membranes, The FEBS journal 280 (2013) 2785

J. A. Campagna, K. W. Miller, and S. A. Forman, Mechanisms of actions of inhaled anesthetics, New England Journal of Medicine 348 (2003) 2110

M. J. Arcario, C. G. Mayne, and E. Tajkhorshid, A membraneembedded pathway delivers general anesthetics to two interacting binding sites in the Gloeobacter violaceus ion channel, Journal of Biological Chemistry 292 (2017) 9480

S. M. Saad et al., Mixed DPPC/DPPG monolayers at very high film compression, Langmuir 25 (2009) 10907

A. K. Sachan et al., High-resolution investigation of nanoparticle interaction with a model pulmonary surfactant monolayer, ACS nano 6 (2012) 1677

J. Hu et al., Investigation of drug for pulmonary administrationmodel pulmonary surfactant monolayer interactions using Langmuir-Blodgett monolayer and molecular dynamics simulation: a case study of Ketoprofen, Langmuir 35 (2019) 13452

H. Bekker et al., Gromacs-a parallel computer for moleculardynamics simulations, In 4th International Conference on Computational Physics (PC 92), World Scientific 4 (1993) 252-256.

M. J. Abraham et al., GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers, SoftwareX 1 (2015) 19

S. Jo et al., CHARMM-GUI Membrane Builder for mixed bilayers and its application to yeast membranes, Biophysical journal 97 (2009) 50

K. Vanommeslaeghe et al., CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields, Journal of computational chemistry 31 (2010) 671

D. J. Evans and B. L. Holian, The nose-hoover thermostat, The Journal of chemical physics 83 (1985) 4069

W. G. Hoover, Canonical dynamics: Equilibrium phase-space distributions, Physical review A 31 (1985) 1695

D. Rose et al., Molecular dynamics simulations of lung surfactant lipid monolayers, Biophysical chemistry 138 (2008) 67

A. Olzynska et al., Mixed DPPC/POPC monolayers: allatom molecular dynamics simulations and Langmuir monolayer experiments, Biochimica et Biophysica Acta (BBA)- Biomembranes 1858 (2016) 3120

C.J. Högberg and A. P. Lyubartsev, Effect of local anesthetic lidocaine on electrostatic properties of a lipid bilayer, Biophysical journal 94 (2008) 525

U. Essmann et al., A smooth particle mesh Ewald method, The Journal of chemical physics 103 (1995) 8577

B. Hess, P-LINCS: A parallel linear constraint solver for molecular simulation, Journal of chemical theory and computation 4(2008) 116

T. Broemstrup and N. Reuter, Molecular dynamics simulations of mixed acidic/zwitterionic phospholipid bilayers, Biophysical journal 99 (2010) 825

G. Patey and J. Valleau, The free energy of spheres with dipoles: Monte Carlo with multistage sampling, Chemical Physics Letters 21 (1973) 297

J. V. N. Ferreira et al., Mechanism of action of thymol on cell membranes investigated through lipid Langmuir monolayers at the air-water interface and molecular simulation, Langmuir 32 (2016) 3234

J. S. Hub, B. L. De Groot, and D. Van Der Spoel, g wham-A Free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates, Journal of chemical theory and computation 6 (2010) 3713

J. A. Lemkul and D. R. Bevan, Assessing the stability of Alzheimer’s amyloid protofibrils using molecular dynamics, The Journal of Physical Chemistry B 114 (2010) 1652

J. G. Kirkwood and F. P. Buff, The statistical mechanical theory of solutions. I, The Journal of chemical physics 19 (1951) 774

D. A. Dougherty, The cation-π interaction, Accounts of chemical research 46 (2013) 885

J. P. Gallivan and D. A. Dougherty, Cation-π interactions in structural biology, Proceedings of the National Academy of Sciences 96 (1999) 9459

B. G. Levine, J. E. Stone, and A. Kohlmeyer, Fast analysis of molecular dynamics trajectories with graphics processing unitsRadial distribution function histogramming, Journal of computational physics 230 (2011) 3556

H. Martinez-Seara and T. Rog, Molecular dynamics simulations of lipid bilayers: simple recipe of how to do it, In Biomolecular Simulations, pp. 407-429 (Springer, 2013).

L. S. Vermeer et al., Acyl chain order parameter profiles in phospho lipid bilayers: computation from molecular dynamics simulations and comparison with 2H NMR experiments, European Biophysics Journal 36 (2007) 919

G. Findenegg, JN Israelachvili: Intermolecular and Surface Forces (With Applications to Colloidal and Biological Systems). Academic Press, London, Orlando, San Diego, New York, Toronto, Montreal, Sydney, Tokyo 1985. 296 Seiten (1986).

D. E. Braun et al., Insights into hydrate formation and stability of morphinanes from a combination of experimental and computational approaches, Molecular pharmaceutics 11 (2014) 3145

M. Paloncyova et al., Amphiphilic drug-like molecules accumulate in a membrane below the head group region, The Journal of Physical Chemistry B 118 (2014) 1030

I. Shehatta, Transfer Thermodynamic Parameters of Local Anesthetic Lidocaine in Ethanol-Water Mixtures. Effect of Preferential Solvation, International journal of research in physical chemistry and chemical physics (2002)

M. Saeedi, A. P. Lyubartsev, and S. Jalili, Anesthetics mechanism on a DMPC lipid membrane model: Insights from molecular dynamics simulations, Biophysical Chemistry 226 (2017) 1

A. Kabedev et al., Molecular dynamics simulations reveal membrane interactions for poorly water-soluble drugs: impact of bile solubilization and drug aggregation, Journal of Pharmaceutical Sciences 110 (2021) 176

D. Lopes et al., Shedding light on the puzzle of drug-membrane interactions: Experimental techniques and molecular dynamics simulations, Progress in lipid research 65 (2017) 24




How to Cite

G. R. Vargas Castañeda, J. A. González, and G. Espinosa, “Molecular dynamics simulations of CBD and THC with a DPPC / DPPG monolayer membrane: Comparative study with morphine and lidocaine”, Rev. Mex. Fís., vol. 69, no. 1 Jan-Feb, pp. 011401 1–, Jan. 2023.



14 Other areas in Physics