Recent experimental studies have shown unexpected chiroptical response from some chiral surfactant molecules, where the specific rotations changed significantly as a function of concentration. To establish a theoretical understanding of this experimentally observed phenomena, a novel methodology for studying chiral surfactants via combined molecular dynamics (MD) and quantum mechanical (QM) calculations is presented. MD simulations on the +10 000 atom surfactant systems have been performed using MD and QM/molecular mechanics (MM) approaches. QM calculations performed on MD snapshots coupled with extensive analysis on lauryl ester of phenylalanine (LEP) surfactant system indicate that the experimentally observed variation of specific rotation with concentration may be due to the conformational differences of the surfactant monomers in the aggregates. Though traditional MM simulations did not show significant differences in the conformer populations, QM/MM simulations using the forces derived from the PM6 method did predict conformational differences between aggregated and nonaggregated LEP molecules, which is consistent with experimental data. Additionally the electrostatic environment of charged surfactants may also be important, since dramatic changes in the Boltzmann populations of surfactant monomers can be noted in the presence of an electric field generated by the chiral ionic aggregates.