The quantum subsystem consisted of the solute described using CAM-B3LYP 64 together with the 6-31+G* basis set. The solvent environment was modeled using parameters from the solvent embedding potential (SEP). 65 The NR molecule was modeled using the fully long-range corrected CAM-B3LYP (α = 0.19, β = 0.81, µ = 0.33) functional and the 6-31G* basis set. Parameters from CP were used for the BLG protein, and SEP was used for the water molecules and ions. Thole damping was applied to dampen interactions between polarizable sites. The entire protein was always included for all models at all cutoff distances.