1712 J ournal of Medicinal Chemistry, 1999, Vol. 42, No. 10
D′Ursi et al.
and the residue chromatographed on a Sephadex LH 20
column (100 cm × 2.5 cm) in water as eluent. Fractions (20
mL) 49-50 afforded 90 mg of cyclic title compound which was
further purified by preparative HPLC using a C18 RP column
eluted by linear gradient H2O/0.1% TFA and CH3CN/0.1% TFA
obtaining 70 mg (36% yield) of pure 12. Fast atom bombard-
ment mass value for (M + H)+ was 539 as expected.
Cyclo[Th r -Lys(2-Cl-Z)-P r o-Ar g(NO2)-â-Ala ] (13). A total
of 0.9 g (1 mmol) of 10 was reacted with 0.136 g of iBCCl
adopting the identical cyclization conditions used for the
synthesis of cyclopentapeptide 11 as previously described.
After purification of the crude reaction material by subsequent
chromatographic procedures on Sephadex LH 20 in methanol
as eluent and on silica gel eluting with chloroform-methanol-
water (65:25:4), 195 mg (25% yield) of cyclic title compound
was recovered (Rf ) 0.45 in chloroform-methanol-water, 65:
25:4).
Cyclo(Th r -Lys-P r o-Ar g-â-Ala ) (14). A total of 195 mg of
cyclopentapeptide 13 in 40 mL of acetic acid-water (1:1) were
hydrogenated over 20 mg of Pd black following the same
procedure used for obtaining cyclopeptide 12. The residue was
purified by chromatography on Sephadex LH 20 in water as
eluent followed by a preparative HPLC purification on a C18
RP column eluted by a linear gradient H2O/0.1% TFA and CH3-
CN/0.1% TFA obtaining 55 mg (40% yield) of pure 14. FAB
(M + H)+ was 553 as expected.
starting structures of cyclic analogues, was modeled starting
from the conformation adopted by the corresponding residues
(T-K-P-R) in human IgG. Peptide coordinates were transferred
directly from the protein X-ray structure using InsightII
Homology module. The obtained molecule was subjected to an
energy minimization, in its zwitterionic form. To allow a
gradual structural relaxation, first a steepest descents mini-
mization using a larger convergence criterion (until maximum
rms derivative was less than 0.5 kcal mol-1 Å-1) was applied;
second the geometry optimization proceeded using conjugate
gradient until the maximum rms derivative was less than
0.001 kcal mol-1 Å-1. This protocol was adopted in each EM
run.
Using InsightII Biopolymer module, tuftsin has been cy-
clized with glycine by automatic addition of the new amino
acid to tuftsin NH terminus. In the case of â-Ala, it proved
necessary to introduce this nonstandard residue manually. The
cyclic models, obtained by these two methods, were subjected
to an unrestrained EM calculation to let the peptide conforma-
tion reach its local energy minimum and then to a restrained
EM calculation by using distance restraints from NOESY
spectra. All restrained calculations were made using a range
of distances according to the percent error calculated for NOE-
derived distances.
Restraints were applied by a “flat-bottomed” potential
function with a force constant value of 1000 kcal mol-1 Å-1
.
NMR. NMR samples were prepared by dissolving appropri-
ate amounts of tuftsin, ctuf2, and ctufâ in a DMSO-d6/water
(80:20, v:v) cryomixture to obtain 1 mM solutions. NMR
spectra were run at 600 MHz on a Bruker DRX-600 instru-
ment. The samples for amide proton exchange were prepared
dissolving the peptides in a DMSO-d6/D2O water (80:20, v:v)
cryomixture to obtain 1 mM solutions.
This function applies a further energy penalty on the atoms
which exceeds the limits of the distance restraints ranges.
Regarding this latter parameter, two different energy penalty
values were assigned, to force more and less strongly the
conformation to fit NOE distances. These values were set to 1
and 100 kcal mol-1 Å-1, respectively, leading to two different
sets of structures defined as ‘less restrained’ and ‘more
restrained’. The obtained structures were relaxed by steepest
descents minimization until maximum rms derivative was less
than 0.5 kcal mol-1 Å-1 to let them reach the nearest local
minimum. The observed and absent NOEs were fully checked.
The conformation with better respected NOEs distances has
been chosen as starting point for restrained molecular dynam-
ics calculations.
1D NMR spectra were recorded in the Fourier mode, with
quadrature detection, and the water signal was suppressed
by a low-power selective irradiation in the homogated mode.
DQF-COSY,21 TOCSY,22 and NOESY23 experiments were run
in the phase-sensitive mode using quadrature detection in ω1
by time-proportional phase incrementation of initial pulse.30
Data block sizes were 2048 addresses in t2 and 512 equidistant
t1 values. Before Fourier transformation, the time domain data
matrices were multiplied by shifted sin2 functions in both
dimensions. A mixing time of 70 ms was employed for the
TOCSY experiment. NOESY experiments were run at 275 K
with mixing times in the range 100-300 ms. Cross-peak
volumes were quantitated in NOESY spectra with 75, 150, 225,
and 300 ms mixing times, using the AURELIA software
(Bruker). Distances were calibrated using d(δ,δ′) Pro3 distances
of 0.180 nm and calculating other distances according to the
r-6 relationship, fitting the build-up curve to a second-order
polynomial function. The necessary pseudoatom corrections
were applied for nonstereospecifically assigned protons at
prochiral centers and for Me group of Thr. The resulting values
were used as input upper distance restraints, considering a
percent error limit of 10%.
A molecular dynamics procedure was performed on the
peptide side chains while keeping the backbone geometry fixed.
Molecular dynamics simulations were run at high temperature
for 300 000 ps so that thermal energy sufficient to cross energy
barriers was added. Local minima conformations were then
equilibrated by running dynamics at 310 K for 600 000 ps to
allow the system to escape from high-energy states.
An average conformation was extracted from molecular
dynamics runs. The whole molecule was then relaxed using a
steepest descents minimization until maximum rms derivative
was less than 0.5 kcal mol-1 Å-1 to let it reach the nearest
local minimum. Observed and absent NOEs were fully checked.
The conformation of “Suzuki’s peptide”, T-P-K-R, was built
using canonical values of the type I â-turn for the internal
rotation angles as suggested in the original paper29 and
subsequently performing a molecular dynamics procedure on
the peptide side chains while keeping the backbone geometry
fixed.
En er gy Ca lcu la tion s. Energy calculations were based on
the all-atom parametrization of the CVFF force field32 as
provided by the MSI InsightII package. Every atom parameter
was assigned explicitly since the CVFF force field is especially
parametrized to reproduce peptide and protein properties. EM
calculations have been performed in vacuo with no distance
cutoff for nonbonded interactions; charges were computed as
provided by the CVFF force field.
Refer en ces
(1) Najjar, V. A.; Nishioka, K. A natural phagocytosis stimulating
peptide. Nature 1970, 288, 672-673.
X-ray coordinates of human immunoglobulin G25 were
extracted from the Brookhaven Protein Data Bank. Model
building was performed according to two different methods
(vide infra). In method (a) starting models for ctuf2 and ctufâ
were built from standard geometries based on canonical
peptide structures containing a cis K-P peptide bond.27 All
molecules were subjected to extensive preliminary EM and
then compared to NMR data. Introduction of NMR-based
restraints led to a good agreement only in the case of the VIa
model, which was thus chosen for subsequent calculations. In
method (b) linear tuftsin, used as a guideline to generate
(2) Nishioka, K.; Amoscato, A. A.; Babcock, G. F.; Banks, R. A.;
Philip, J . H. An immunomodulating peptide hormone and its
clinical potential as a natural biological response modifier.
Cancer Invest. 1984, 2, 39-49.
(3) Fridkin, M.; Najjar, V. A. Tuftsin: Its chemistry, biology and
physiological function. Crit. Rev. Biochem. Mol. Biol. 1989, 24,
1-40.
(4) Wagle, J . R.; Anservin, A. T.; Dessens, S. E.; Nishioka, K. Specific
translocation of Tuftsin (Thr-Lys-Pro-Arg), a natural immuno-
modulating peptide, into the nuclei of human monocytes. Bio-
chem. Biophys. Res. Commun. 1989, 159, 1147-1153.
(5) Constantopoulos, A. Congenital tuftsin deficiency. Ann. N.Y.
Acad. Sci. 1983, 419, 214-219.