(CaGly(1)), 55.2 (CaPhe), 118.7 (CbDEPhe), 119.0 + 124.8
(aromatic ring of pNA), 129.1 + 126.4 (aromatic rings), 130.8
(CaDEPhe), 137.2 + 134.0 (aromatic rings), 144.7 + 142.4
(aromatic ring of pNA), 162.54 (CO DEPhe), 164.9 (CO
Gly(3)), 168.4 (CO Gly(1)), 170.9 (CO Phe).
the next steps the side chain atoms from active side were
unfrozen and all atoms of the dehydropeptide were analysed
using the conjugate gradient algorithm and tether constraints
on the unfrozen heavy atoms, which was decreased from a
value of 200 to 20 kcal AÀ2, also the maximum derivative was
decreased from a value of 30 to 1 kcal AÀ1. In the last step we
optimized the complex without any force constant using a
conjugate gradient algorithm with the maximum derivative
equal to 0.1 and then 0.01 kcal AÀ1. To decide which arrange-
ments were the best for each dehydropeptide, we calculated the
Gly-Gly-DEPhe-Phe-pNA. This was obtained in identical
manner as described above. Yield: 0.323 g (98%). Elemental
analysis: calc. for C28H37O6N6 (553.52): C 60.7, H 6.73; found:
C 50.98, H 6.62%.
1H NMR (DMSO): 3.15 + 2.9 (m, 2H, CbH2 Phe), 3.67 (s,
2H, Ca H2Gly(1)), 4.04 (d, 2H, Ca H2Gly(2)), 4.77 (q, 1H, CaH
Phe), 6.55 (s, 1H, CbH DEPhe), 7.05 + 7.30 (m, 10H, aromatic
rings), 7.97 + 8.25 (2 Â d, 4H, aromatic ring of pNA), 8.14 (s,
2H, NH2 Gly(1)), 8.72 (d, 1H, NH Phe), 8.86 (t, 1H, NH
energy of interaction Eint, which was defined as Eint = Ecomp
À
(Een + Elig), where Ecomp was the energy of the enzyme–ligand
complex excluding the water layer; Een was the energy of the
enzyme; Elig was the energy of the ligand.
Gly(2)) 10.24 (s, 1H, NH DEPhe); 10.26 (s, 1H, NH pNA); 13
C
NMR (DMSO): 36.6 (CbPhe), 40.3 (CaGly(1)), 42.2
(CaGly(2)), 55.3 (CaPhe), 118.6 (CbDEPhe), 119.3 + 124.9
(aromatic ring of pNA), 129.2 + 126.5 (aromatic rings), 131.4
(CaDEPhe), 134.2 (aromatic ring of DEPhe), 137.4 (CgPhe),
144.9 + 142.4 (aromatic ring of pNA), 164.5 CO (DEPhe),
166.6 (CO Gly(1)), 167.9 (CO Gly(2)), 170.4 (CO Phe).
Acknowledgements
The molecular modeling were carried out using hardware and
software resources of The Supercomputing and Networking
Centre in Wroclaw.
References
Assay of the enzymatic activity
1 D. E. Palmer, C. Pattaroni, K. Nunami, R. K. Chadha, M.
Goodman, T. Wakamiya, K. Fukase, S. Horimoto, M. Kitazawa,
H. Fujita, A. Kubo and T. Shiba, J. Am. Chem. Soc., 1992, 114,
5634.
2 H. Allgaier, G. Jung, R. G. Werner, U. Schneider and H. Zamer,
Angew. Chem., 1985, 24, 1051.
3 O. Pieroni, A. Fissi, R. M. Jain and V. S. Chauhan, Biopolymers,
1996, 38, 141.
4 C. M. Venkatachalam, Biopolymers, 1968, 6, 1425.
5 S. M. Patel, J. O. Currie, Jr. and R. K. Olsen, J. Org. Chem., 1973,
38, 126.
Cathepsin C was isolated according to the described proce-
dure.22 Enzymatic reaction was assayed at 37 1C in acetate
buffer (pH 5) containing NaCl (10 mM final concentration)
and 2-mercaptoethanol (5 mM final concentration). The assay
mixture contained dehydrotetrapeptide (1–8 mM final con-
centration) and the course of reaction was monitored by
following the change in absorbance at 400 nm. Michaelis
constants (KM) and maximal velocities (Vmax) were obtained
using the computer program kindly provided by Dr J. Hurek
(University of Opole).
6 K. R. Rajashankar, S. Ramakumar and V. S. Chauhan, J. Am.
Chem. Soc., 1992, 114, 9225.
7 M. R. Ciajolo, A. Tuzi, C. R. Pratesi, A. Fissi and O. Pieroni,
Biopolymers, 1990, 30, 911.
Molecular modelling
8 Y. Inai, S. Kurashima, T. Hirabayashi and K. Yokota, Biopoly-
mers, 2000, 53, 482.
9 B. Turk, I. Dolenc and V. Turk, Handbook of Proteolytic Enzymes,
ed. A. J. Barrett, Academic Press, New York, 1998, pp. 631–634.
10 R. Butler, A. Michel, W. Kunz and M. O. Klinkert, Protein Peptide
Lett., 1995, 2, 313.
11 C. Toomes, J. James, A. J. Wood, C. L. Wu, D. McCormick, N.
Lench, C. Hewitt, L. Moynihan, E. Roberts, C. G. Woods, A.
Markham, M. Wong, R. Widmer, K. A. Ghaffar, M. Pemberton, I.
R. Hussein, S. A. Temtamy, R. Davies, A. P. Sloan, M. J. Dixon
and N. S. Thakker, Nat. Genet., 1999, 23, 421.
12 A. O. Breitholle and C. H. Stammer, Tetrahedron Lett., 1975, 28,
2381.
13 M. Makowski, M. Pawe"czak, R. Latajka, K. Nowak and P.
Kafarski, J. Pept. Sci., 2001, 7, 141.
14 M. Lisowski, R. Latajka, M. Makowski, T. Lis, B. Picur, J. Panek
and P. Kafarski, J. Am. Chem. Soc., submitted.
15 M. Makowski, B. Rzeszotarska, Z. Kubica, G. Pietrzynski and J.
Hepter, Liebigs Ann. Chem., 1986, 980.
16 M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A.
Robb, J. R. Cheeseman, J. A. Montgomery, Jr., T. Vreven, K. N.
Kudin, J. C. Burant, J. M. Millam, S. S. Iyengar, J. Tomasi, V.
Barone, B. Mennucci, M. Cossi, G. Scalmani, N. Rega, G. A.
Petersson, H. Nakatsuji, M. Hada, M. Ehara, K. Toyota, R.
Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O.
Kitao, H. Nakai, M. Klene, X. Li, J. E. Knox, H. P. Hratchian, J.
B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E.
Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J.
Ochterski, P. Y. Ayala, K. Morokuma, G. A. Voth, P. Salvador, J.
J. Dannenberg, V. G. Zakrzewski, S. Dapprich, A. D. Daniels,
The structures of the studied dehydropeptides were optimized
in program Gaussian03 at the HF/6-31g (d,p) level16 in the gas
phase using the Merz–Singh–Kollman scheme20 in the deter-
mination of the atomic charges. The calculation of the docking
process were performed using the AutoDock program.18 The
starting geometry of the dehydropeptides was taken from the
ab initio calculations, and also we assigned the charges using
charges from the ab initio calculation. The structure of cathe-
psin C was taken from the structure of human dipeptidyl
peptidase I deposited EC 3.4.14. in the Protein Data Bank.17
During the docking process the main chain of the dehydro-
peptide was fixed, whereas side chains and the terminal groups
(–NH2, –pNA, –Boc) were left as flexible. The coordinates of
the SH proton from the Cys234 were taken as a grid center in
the docking process. Several possible structures ligand–
enzyme complexes for each dehydropeptide were obtained in
this manner, which were grouped in clusters. For the next
stage we chose structures from the most possible clusters. In
the next step the selected complex compounds were optimized
using Accelrys’s DISCOVER program with the cff97 force
field, at neutral pH and considering a 10 A water layer. In the
first step all heavy atoms were frozen and a steep descent
algorithm with maximum derivative equal to 0.1 was used. In
ꢀc
This journal is the Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2006
New J. Chem., 2006, 30, 1009–1018 | 1017