Molecular modeling of hair keratin/peptide complex: Using MM-PBSA calculations to describe experimental binding results.

Molecular dynamics simulations of a keratin/peptide complex have been conducted to predict the binding affinity of four different peptides toward human hair. Free energy calculations on the peptides' interaction with the keratin model demonstrated that electrostatic interactions are believed to be the main driving force stabilizing the complex. The molecular mechanics-Poisson-Boltzmann surface area methodology used for the free energy calculations demonstrated that the dielectric constant in the protein's interior plays a major role in the free energy calculations, and the only way to obtain accordance between the free energy calculations and the experimental binding results was to use the average dielectric constant.