Three 1 ns length molecular dynamics simulations of an RGD peptide (Ac-Pen-Arg-Gly-Asp-Cys-NH2, with Pen denoting penicillamine) have been performed in aqueous solution, one for the disulfide bridged, and two for the unbridged form. The trajectories were analyzed to identify conformations explored by the two forms and to calculate several properties: NMR vicinal coupling constants, order parameters, dipole moments and diffusion coefficients, in an effort to describe the physical role of the disulfide bond. The cyclic peptide was able to explore several distinct backbone conformations centered around a turn-extended-turn structure. However, its flexibility was limited and it appeared to be 'locked in' into a a family of structures characterized by a high dipole moment and a well-defined conformation of the pharmacophore, which has been previously identified as biologically active. Excellent agreement between the simulated and observed NMR vicinal coupling constants indicates that realistic structures were sampled in the cyclic peptide simulation. The linear form of the peptide was much more flexible than the cyclic one. In the two independent 1 ns simulations of the linear form the explored conformations could be roughly grouped into two classes, of cyclic-like and extended type. Within each simulation the peptide switched between the two classes of structures several times. Exact matches between conformations in the two linear peptide simulations were not found; several conformational regions with backbone rms deviations below 1A were identified, suggesting that representative structures of the linear form have also been identified. In the linear peptide simulations the RGD pharmacophore is able to adopt a wide range of conformations, including the one preferred by the cyclic form. The lower biological activity of the linear peptide compared to the cyclic one may be correlated with the lower population of this structure in the absence of the disulfide bond.