Predicting proton titration in cationic micelle and bilayer environments

Knowledge of the protonation behavior of pH-sensitive molecules in micelles and bilayers has significant implications in consumer product development and biomedical applications. However, the calculation of pKa's in such environments proves challenging using traditional structure-based calculations. Here we apply all-atom constant pH molecular dynamics with explicit ions and titratable water to calculate the pKa of a fatty acid molecule in a micelle of dodecyl trimethylammonium chloride and liquid as well as gel-phase bilayers of diethyl ester dimethylammonium chloride. Interestingly, the pKa of the fatty acid in the gel bilayer is 5.4, 0.4 units lower than that in the analogous liquid bilayer or micelle, despite the fact that the protonated carboxylic group is significantly more desolvated in the gel bilayer. This work illustrates the capability of all-atom constant pH molecular dynamics in capturing the delicate balance in the free energies of desolvation and Coulombic interactions. It also shows the importance of the explicit treatment of ions in sampling the protonation states. The ability to model dynamics of pH-responsive substrates in a bilayer environment is useful for improving fabric care products as well as our understanding of the side effects of anti-inflammatory drugs.