All QM single-point energies,
angle scans, and torsion scans were performed using Gaussian 09.
15 Angle and torsion scans on model hydrocarbon
molecules were performed at the MP2/cc-pVTZ level. Angles were scanned
from 90 to 150° in 1° increments. Torsions were scanned
over 360° in 10° increments. Single-point energy calculations
on model glyceride, ceramide, and phosphatidylcholine molecules were
performed at the MP2/cc-pVDZ level, with application of the polarizable
continuum model to create an implicit solvent environment for phosphatidylcholine
only. To obtain relevant glyceride, ceramide, and phosphatidylcholine
conformations, POPC and PSM simulations were performed for 100 ns
using initial Lipid21 parameters, and 1000 random lipid structures
were extracted. Prior to the QM single-point energy calculation, each
molecule was minimized using AMBER20
16 for
1000 steps, with the first 500 steps using the steepest descent and
the final 500 steps using the conjugate gradient method,
17 with restraints of 5000 kcal/mol/rad on each
of the torsions being fitted. The GBn generalized Born model (igb
= 7)
18 (link) was used during minimization of
model phosphatidylcholine molecules only. We repeated this process
on separate trajectories to create a test set, used for parameter
validation. Parameters for the model glyceride, ceramide, and phosphatidylcholine
molecules consisted of initial Lipid21 parameters and partial charges
derived at the MP2/cc-pVDZ level after optimization of a single molecular
conformation, allowing a two-stage RESP fit.
19 (link)
Dickson C.J., Walker R.C, & Gould I.R. (2022). Lipid21: Complex Lipid Membrane Simulations with AMBER. Journal of Chemical Theory and Computation, 18(3), 1726-1736.