Molecular Dynamics Simulations of Conformational Behavior of Linear RGD Peptidomimetics and Cyclic Prodrugs in Aqueous and Octane Solutions

Abstract

Conformations available to a class of cyclic prodrugs and corresponding linear RGD peptidomimetics were explored using 1 ns length molecular dynamics simulations performed with the program CHARMM. Water and octane, modeled explicitly, were used as solvents to mimic the change of the environment experienced by the solutes upon partition from water to membrane in the trans-cellular transport process. In water, the linear peptidomimetics tended to populate extended-like structures, characterized by strong favorable interactions with solvent and low intrinsic stability. In these extended conformations the charged termini are able to assume large distances, above 15 Å for the longest systems. These linear peptidomimetics have been found to exhibit the highest potency in experimental studies, in accord with the trends experimentally observed for RGD peptides. In contrast, in octane compact conformers of the linear peptidomimetics were favored, with all charged groups aggregated and shielded from solvent, exhibiting high intrinsic stability and weak solute-solvent interactions. Our calculations predict a large unfavorable energy change for transferring the linear systems from water to octane, in agreement with experimental findings that these compounds are not transported via the trans-cellular pathway. The cyclic pro- drugs did not exhibit major structural differences between the simulations in water and octane, adopting turn-like conformations in both solvents. The limited response of the cyclic structures to changes in the environment leads to energies of transfer from water to octane that are also unfavorable, but much less so than for the linear molecules. This effect is in accord with the observed enhanced passive trans-cellular transport of the cyclic prodrugs.