Molecular dynamics simulations of a hydrated protein vectorially oriented on polar and nonpolar soft surfaces

We present a collection of molecular dynamics computer simulation studies on a model protein-membrane system, namely a cytochrome c monolayer attached to an organic self-assembled monolayer (SAM). Modifications of the system are explored, including the polarity of the SAM endgroups, the amount of water present for hydration, and the coordination number of the heme iron atom. Various structural parameters are measured, e.g., the protein radius of gyration and eccentricity, the deviation of the protein backbone from the x-ray crystal structure, the orientation of the protein relative to the SAM surface, and the profile structures of the SAM, protein, and water. The polar SAM appears to interact more strongly with the protein than does the nonpolar SAM. Increased hydration of the system tends to reduce the effects of other parameters. The choice of iron coordination model has a significant effect on the protein structure and the heme orientation. The overall protein structure is largely conserved, except at each end of the sequence and in one loop region. The SAM structure is only perturbed in the region of its direct contact with the protein. Our calculations are in reasonably good agreement with experimental measurements (polarized optical absorption/emission spectroscopy, x-ray interferometry, and neutron interferometry).