One NR molecule, which is centered in the simulation box, comprises the quantum subsystem. It was described at the same level of theory as in the BLG protein, i.e., fully long-range corrected CAM-B3LYP/6-31+G*. The lipids were modeled using parameters from the averaged lipid embedding potential (ALEP), parameters for ions and water molecules were from SEP, and the remaining NR molecules were modeled by LoProp 69 multipoles (up to and including quadrupoles) and dipole-dipole polarizabilities calculated at B3LYP/6-31+G* level of theory (using an ANO-type recontraction of 6-31+G*).