260 ps molecular dynamics simulation of substance P with hydrated dimyristoyl phosphatidyl choline bilayer

We present here results on 260 pico seconds (ps) molecular dynamics (MD) simulation of substance P (SP) in hydrated bilayer of dimyristoyl phosphatidyl choline (DMPC) (39 molecules of DMPC with 776 water molecules). 260 ps MD simulation has been carried out in 0.001 ps time interval with united atom force field, using AMBER 4.0 package. Non bonded pair list was updated every 20 cycles using 12.5 Angstrom cut off distance. Analysis of MD data is done using our package ANALMD. The obtained models are presented using graphics package RASMOL. All simulations, analysis of MD data and graphics is done on INDIGO-2, R-4400 extreme graphics work station. Our results show no systematic change in order parameter, but reduction in transfraction of the chain torsional angles, compared to our earlier results on MD simulation on hydrated DMPC bilayer without SP. C-terminal and central peptide residues adopt partial helical conformation. Helix type as classified on the basis of H-bonds is between alpha and 3(10). The peptide backbone shows flexibility during heating runs. Later, it is stabilized and there was not much change in the spatial position of the backbone. Lipid matrix serves the role of immobilization of the peptide backbone in a preferred conformation.