Molecular simulation studies of the structure of phosphorylcholine self-assembled monolayers

We report a study of the structure of phosphorylcholine self-assembled monolayers (PC-SAMs) on Au(111) surfaces using both molecular mechanics (MM) and molecular dynamics (MD) simulation techniques. The lattice structure (i.e., packing densities and patterns) of the PC chains was determined first, by examining the packing energies of different structures by MM simulations in an implicit solvent. The chain orientation (i.e., antiparallel and parallel arrangements of the PC head groups) was then evaluated. The initial azimuthal angles of the PC chains were also adjusted to ensure that the optimal lattice structure was found. Finally, the two most probable lattice structures were solvated with explicit water molecules and their energies were compared after 1.5 ns of MD simulations to verify the optimal structures obtained from MM. We found that the optimal lattice structure of the PC-SAM corresponds to a radical7 x radical7 R19degree lattice structure (i.e., surface coverage of 50.4 A(2)molecule) with a parallel arrangement of the head groups. The corresponding thickness of the optimal PC-SAM is 13.4 A which is in agreement with that from experiments. The head groups of the PC chains are aligned on the surface in such a way that their dipole components are minimized. The P-->N vector of the head groups forms an angle of 82 degrees with respect to the surface normal. The tilt direction of molecular chains was observed to be towards their next nearest neighbor.