서 론
백금 및 팔라듐착물의 항암활성에 관한 연구는 실험적으로나 이론적으로 많이 연구되고 있으나, 니켈착물에 대해서는 실험적으로나 이론적으로 아직까지 미흡한상태이다. Zakharova등1은 ATPase가결합된 membrane에 대한 백금(II) 및 팔라듐(II)착물들의 저해효과는 산성리간드의 주개-받개 성질에 의존한다고 보고하였다.
그러나 항종양제의 대부분은 유기화합물 또는 alkylating agents, antibiotic alkaloids, enzymes, hormones등을 포함하는 자연생성물로써 알려졌으나,2 금속을 포함하는 무기화합물에 있어서는 대부분의 금속들을 잠재적인 발암제라고 생각했을 뿐 그다지 체계적으로 연구되지는 못했다. 니켈,3 팔라듐4 및 백금착물5에서 니켈 화합물이 발암제로 알려졌으나, 몇 가지 니켈 착물들에 대해 실험한 결과단지 bis(diethyldithiophosphato) Nickel(II) 착물만이 Walker 256 carcinosarcoma에대항할 만큼 백금 및 팔라듐 착물 유도체에 대해서 보다 더 강한 항암활성임이 보고되었다.6 또한 팔라듐(II) 착물7을 동물속에 투입시켜 항암성 정도를 실험하였으나, 같은 구조를 가진 백금(II) 착물의 항암활성보다는 일반적으로 더 낮음이 보고되었다.8,9 한편 Tang등10은 cis-[Pt(am)2Cl2] (am=amines)형 착물에 대해 CNDO/2법으로 계산한 결과 전자배치와 항암활성과의 관계를 Pt-N간의 결합이 강할수록, Cl원자의 알짜전하가 음의 값으로 클수록 항암활성에 크게 관여하고 있다고 발표하였다. Qu등11은 백금(II) 아민 착물의 항암활성과 구조 사이의 관계에 대해서 CNDO/2법으로 계산된 양자화학적 양과 관련지워 설명하였는데, 그 결과 독성(toxicity)과 착물의활성은 아민과 라디칼리간드에 밀접한관련이 있다고 보고하였다. 그러나 Cu(II)와 Ni(II) 착물의 저해효과에 대해서는 Tatjanenko등12이 polyamines과 imidazole 유도체가배위된 Cu(II) 및 Ni(II) 착물의 Ca와 Mg에의존하는 ATPase 기질에 대해 저해활성을 조사한 결과 안정한 사각평면형, 팔면체형 및 5배위된 Cu(II) 착물과 사면체형, 6배위된 Ni(II) 착물에서는 저해활성이 약해지고 나머지 다소 불안한 Cu(II) 및 Ni(II) 착물에 대해서는 저해활성값이 크게 나타남을 보고하였다. 그러나 이론적 연구는 미비한 상태이다. 따라서본연구에서는 polyamines와 imidazole 유도체가 배위된 Ni(II)착물의항암활성을 양자화학적인 양인 형성엔탈피, 결합에너지, 중심금속인 Ni(II)이온의 알짜전하, 쌍극자힘등과 관련지워 착물의 기하학적 구조의 차이, 리간드 치환체의 수와 입체효과가 어느정도 영향을 미치는지에 대해 ZINDO/1 (Zerner’s Intermediate Neglect of Differential Overlap)법13으로 조사하고자 한다.
이론 및 계산
계산방법-반경험적 분자궤도함수 계산
자유리간드인 polyamine류와 imidazole 유도체에 대해서는 AM1(Austin Model No.1)법14과 PM3(Parametric Method No.3)법15을이용하여계산하였다. Polyamine류와 imidazole 유도체가 배위된 Ni(II)착물에 대해서는 우선 분자역학(MM+)법16으로 최적화시키고, 최적화된 Ni(II)착물에 대해서 ZINDO/1(Zerner’s Intermediate Neglect of Differential Overlap No.1)법13으로 다시 최적화시킨후 single point계산을 실행하였다. PM3는 반경험적인 SCF-MO법으로 AM1방법과 동일한 연산방식을 사용하지만 AM1방법에서 사용하는 변수를 재변수화한방법으로 AM1법과 함께다소정확한 계산방법 중의 하나로 알려져 있다. AM1, PM3 및 ZINDO/1 계산은 Hyperchem package(V 5.0)17를이용하여계산하였고, ab initio법은제한적 Hartree-Fock(RHF/6-31G*)법으로하였고, 최적화조건의연산은 Polak-Ribiere법을이용하여 RMS gradient를 0.001 kcal/Å·mol로하여계산하였다. ZINDO/1(Zerner’s INDO/1)13은 INDO(Intermediate Neglect of Differential Overlap)의 변형으로 전이금속에대해 EHMO(Extended Hückel Molecular Orbital)법보다 다소 정확한 계산법으로 일정한 궤도함수 지수를 사용함으로써 INDO와 INDO/1의가장 최근버전과 동일하다. 특히, ZINDO/1 계산법은 d-궤도함수를 포함한 전이금속착물에대해그에너지상태와기하학적구조를 최적화하는 계산뿐만 아니라 수용액 상태의 분자구조의 연구를위한 SCRF(Self-Consistent Reaction Field) 계산과 함께 atomic partial charge, bond order, UV-visible spectra, reactivity surface 연구를수행하는데 많이 이용된다.
모델 설정
자유 리간드(free ligands) 종류
자유리간드인 폴리아민류와 이미다졸유도체에 대한 화학적 구조를 Scheme 1에 나타내었다.
Scheme 1.Chemical structures for seven free ligands, polyamines and imidazole derivatives.
최적화된 자유 리간드 분자 구조
AM1 계산에 의한 최적화된 7개의 자유리간드인 폴리아민류와 이미다졸 유도체에 대한 분자 구조를 Scheme 2에 나타내었다.
Scheme 2.Optimized structures for seven free ligands, polyamines and imidazole derivatives at the AM1 calculation.
결과 및 고찰
Ni(II)착물의 저해활성 정도를 조사하기 위하여 우선 Ni(II)착물에 배위되는 자유 리간드인 폴리아민류와 이미다졸 유도체에대해서 반경험적인 AM1, PM3법으로 양자화학적양을 얻었다. 즉, 형성엔탈피(ΔHf), HOMO, LUMO 에너지및 쌍극자힘(μD)의 값을 Table 1에나타내었다.
Table 1.AM1, PM3 enthalpy of formation(ΔHf), HOMO and LUMO energy, and dipole moment(μD) of polyamines and imidazole derivatives
Table 1에 의하면 자유리간드인 폴리아민류와 이미다졸유도체들은 전자주개이므로 전자를 줄 수 있는 능력의 척도인 HOMO 에너지값을 살펴보면, AM1 및 PM3에서 폴리아민류와 이미다졸 유도체 모두에서 대체로 -9.30 eV 정도로나타났다. HOMO 에너지와전자를 받을 수 있는 능력의 척도인 LUMO 에너지와의 차이 값은 폴리아민류에서는 대략 12.20 eV 정도로 나타났고, 이미다졸유도체인 N-PropIm과 N-EtIm에서는 10.02 eV 정도로 폴리아민류보다 다소 작은 값을 보였다. 이는 자유 리간드 상태에서는 대체로 폴리아민류가 이미다졸 유도체 보다 더 안정하다는 것을 알 수 있다. 또한 HOMO-LUMO에너지 차이 값(ΔE)은 N-PropIm 또한 극성의 정도를 나타내는 쌍극자힘(μD)은 이미다졸유도체인 N-PropIm와 N-EtIm에서 4.10 Debye(D)로서 폴리아민류 보다 더 큰 값을 보여서 더 극성임을 알 수 있었다. 다음은 Ni(II) 이온에 이들 자유 리간드인 폴리아민류와 이미다졸 유도체들을 배위시켜 ZINDO/1법 계산에 의해 양자화학적 양인 각 착물의 전체에너지(Etot), 결합에너지(Ebind), 형성엔탈피(ΔHf), HOMO, LUMO에너지, 쌍극자힘(μD) 값을 Table 2에 나타내었다. Table 2.ZINDO/1 calculated quantum chemical quantities for Ni(II) complexes with polyamines and imidazole derivatives Table 3.a)Ref. 12. Table 2에서 Ni(II)착물의 각 원자 간의 σ- 및 π-결합에 의한 즉, 겹칩적분과 공명적분에 의해 얻어진 총결합에너지(Ebind)는 [Ni(trien)(en)]2+<[Ni(dpt)(tn)]2+ <[Ni(tn)3]2+<[Ni(dpt)(tn)(NCS)]+<[Ni(N-PropIm)2(NCS)2]<[Ni(Me6tren)(NCS)]+<[Ni(N-EtIm)6]2+<[Ni(N-PropIm)6]2+순으로 나타났다. 결합강도의 정도를 나타내는 총결합에너지는 폴리아민류가 배위된 Ni(II)착물에서 보다 이미다졸 유도체인 N-PropIm과 N-EtIm이 Ni(II) 이온에 배위되었을 때 각각 -1361.148와 -1169.069 eV로서 음의 값으로 크게 나타났다. 이는 Table 1에서 보여진 자유리간드 자체인 N-PropIm과 N-EtIm의 HOMO-LUMO 에너지 차가 다른폴리아민류의 값보다대체로 작은값이 나타나서 즉, 자유 리간드 자체의 HOMO-LUMO에너지 차이값(ΔE)이 작을수록 Ni(II)착물의 총결합에너지가 크게 나타났다. 이는 자유리간드 자체가 불안정할수록 Ni(II)이온에 더 강하게 결합함을 알 수가 있었다. 형성엔탈피(ΔHf) 값은 [Ni(N-PropIm)2(NCS)2]<[Ni(dpt)(tn)]2+<[Ni(tn)3]2+<[Ni(dpt)(tn)(NCS)]+<[Ni(trien)(en)]2+<[Ni(Me6tren)(NCS)]+<[Ni(N-EtIm)6]2+<[Ni(N-PropIm)6]2+순으로 나타났다. 또한 HOMO-LUMO 에너지 차이의 값은 [Ni(N-PropIm)2(NCS)2]<[Ni(N-PropIm)6]2+<[Ni(NEtIm)6]2+<[Ni(Me6tren)(NCS)]+<[Ni(dpt)(tn)(NCS)]+<[Ni(tn)3]2+<[Ni(trien)(en)]2+<[Ni(dpt)(tn)]2+순으로 증가하였다. 다음은 중심금속 Ni(II)의 알짜전하(QNi), Ni(II)착물의 형성엔탈피(ΔHf)와 Ni(II)착물에 대한 저해활성(IA)값 사이의 상관성을 알아보기 위하여 이들 값들을 Table 3에 나타내었다. Table 3을 살펴보면 폴리아민류와 이미다졸유도체가 배위된 Ni(II) 착물에서 중심금속 Ni(II)의 원자전하는 0.020에서 0.084까지의 값을 나타내었다. Ni(II) 착물에서 중심금속 Ni(II)의 원자전하 값은 [Ni(N-PropIm)2(NCS)2]<[Ni(dpt)(tn)]2+<[Ni(Me6tren)(NCS)]+< [Ni(trien)(en)]2+<[Ni(tn)3]2+<[Ni(N-PropIm)6]2+<[Ni(dpt) (tn)(NCS)]+<[Ni(N-EtIm)6]2+ 순으로 나타났다. 여기서 Ni(II)이온의 원자전하 값과 저해활성과의 관계를 살펴보면, 5배위 사각피라미드형인 [Ni(N-PropIm)2(NCS)2]와 일그러진 사면체 구조인 [Ni(dpt)(tn)]2+ 착물에서 Ni(II) 이온의 원자전하 값이 각각 0.020과 0.029로 가장낮게 나타나서 다른 Ni(II) 착물에 비해 저해활성값이 거의 없다는 것을 알 수가 있었다. 다음은 Ni(II)착물의 형성엔탈피(ΔHf)와저해활성(IA)과의 상관성을 살펴보면 형성엔탈피가 음의 값으로 적게 나타난 -263.606 eV인 [Ni(N-PropIm)2(NCS)2]와 -286.188 eV인 [Ni(dpt)(tn)]2+착물에서 저해활성 값이 거의 없는 것으로 나타났다. 반면에 형성엔탈피 값이 음의 값으로 -300 eV 이상인 착물에서 저해활성정도가 대체로 25% 이상인 값을 가진다는 것을 알 수 있다. Table 4.aGuanine + 2'-Deoxyribose
bGuanine + 2'-Deoxyribose + triphosphate 다음은 이들 Ni(II)착물들과 DNA 및 RNA 염기쌍들과의 전자전이 상관관계를 알아보기 위하여 DNA 및 RNA 염기쌍에 대해 AM1, PM3 및 RHF/6-31G*법에 의한 HOMO와 LUMO 에너지 값을 Table 4에 나타내었다. Table 4에서 DNA 및 RNA 염기쌍의 HOMO 에너지값은 AM1법계산에서는 -9.971 eV에서 -8.862 eV까지로 Uracil < Thymine < Cytosine < Adenine < Guanine 순으로에너지가 높아졌다. PM3법계산에서는 -9.710 eV 에서 -8.821 eV까지로서 Uracil < Thymine < Cytosine < Adenine < Guanine 순으로 에너지가 높아졌다. RHF/6-31G* 계산에서는 -9.852 eV에서 -7.912 eV까지 Uracil < Thymine < Cytosine < Adenine < Guanine 순으로 에너지가 높아졌다. AM1, PM3, RHF/6-31G* 계산 모두에서는 같은 양상을 보였다. 이들 계산 방법 모두에서 전자를 줄 수 있는 능력의 척도인 HOMO 에너지가 가장 낮은 염기쌍은 Uracil이고, 가장 높은 염기쌍은 Guanine으로 나타났다. 따라서 Uracil보다는 Guanine 염기쌍이 Ni(II)착물에 전자를 제공하기가 더 용이함을 알 수 있다. 따라서 DNA 및 RNA 염기쌍의 HOMO 에너지에서 Ni(II)착물의 LUMO 에너지로 전자 전이가 일어나서 분자착물을 형성하여 저해활성 및 항암활성에 관여할 것으로 생각된다. DNA 및 RNA의 염기쌍의 HOMO 에너지와 폴리아민류와 이미다졸 유도체가 배위된 Ni(II) 착물의 LUMO 에너지 사이의 전자전이의정도를알아보기위하여 Fig. 1에도시하였다. Fig. 1.The different between HOMO energy of DNA and RNA bases(by AM1) and Ni(II) complexes with polyamines and imidazole derivatives(by ZINDO/1). Fig. 1에서 보면 DNA 및 RNA 염기쌍의 중에 Uracil의 HOMO 에너지와 Ni(II) 착물의 LUMO 에너지 사이의 에너지 차이가 가장 크게 나타났으며, Guanine의 HOMO 에너지와 Ni(II) 착물의 LUMO 에너지 사이의 차이가 가장 작게 나타남을 볼 수 있다. 따라서 5개의 염기쌍중에 Guanine 염기쌍에서 Ni(II) 착물로전자전이가 가장 용이하게 일어날 수 있다는 알 수 있다. 폴리아민류와 이미다졸 유도체가 배위된 Ni(II) 착물의 중심금속의 알짜전하가 0.03이상에서저해활성값이 나타났고, 또한형성엔탈피(ΔHf)와 저해활성(IA)과의 상관성을 살펴보면 형성엔탈피값이음의값으로 -300 eV 이상인 착물에서 저해활성 정도가 대체로 25% 이상인 값을 가지는 것을 알 수 있었다. DNA 및 RNA 염기쌍과 Ni(II) 착물 사이의 관계에서는 우선, DNA 및 RNA 염기쌍 자체에서 전자를 줄 수 있는 능력의 척도인 HOMO 에너지가 가장 낮은 염기쌍은 Uracil이고, 가장 높은 염기쌍은 Guanine으로 나타났다. 따라서 Uracil보다는 Guanine 염기쌍이 Ni(II)착물에 전자를 제공하기가 더 용이함을 알수있었다. 따라서 DNA 및 RNA 염기쌍의 HOMO 에너지에서 Ni(II)착물의 LUMO 에너지로 전자 전이가 일어나서 분자착물을 형성하여 저해활성 및 항암활성에 관여할 것으로 생각된다. 본 연구는 2001년도 대구대학교 학술연구비 지원에 의하여 수행된 바 이에 감사드립니다.결 론
References
- Zakharova, I. A.; Tatjanenko, L. V.; Moshkovsky, Y. S.;Raykhman, L. M.; Kondratjeva, T. A. Biofisica 1977,22, 418.
- Haiduc, I.; Silverstru, C. Coord. Chem. Rev. 1990, 99, 253. https://doi.org/10.1016/0010-8545(90)80065-2
- Williams, R. D. Chem. Rev. 1972, 72, 203. https://doi.org/10.1021/cr60277a001
- Cleare, M. J. Coord. Chem. Rev. 1974, 12, 349. https://doi.org/10.1016/S0010-8545(00)82029-9
- Sttetsento, A. I.; Presnov, M. A.; Konovalova, A. L. Chem. Abstr.1981, 95, 46.
- Connors, T. A. In Platinum Coordination Complexesin Cancer Chemotherapy; Roborts, J. J., Ed.; Springer-Verlag: Berlin, Germamany, 1974; p 237.
- Scovell, W. M.; Kroos, L. R.; Capponi, V. J. In Platinum Gold and Biochemotherapeutic Agents; Chemistry and Biochemistry; Lippard, S. J., Ed.; ACS Symposium Series: Washington, D. C., 1983, 209, 101. https://doi.org/10.1021/bk-1983-0209.ch005
- Benesch, E.; Harly, H. A.; Benesch, R. J. Biol. Chem. 1965, 216, 663.
- Hacker, M. P.; Douple, E. B. In Platinum CoordinationComplexes in Cancer Chemotherapy Development inOncology Series; Krakoff, I. M., Ed.; Martinus Nijhoff:Boston, MA, U. S. A., 1984, 17, 127.
- Livingstone, S. E.; Mihkelson, A. E. Inorg. Chem. 1970, 9, 2545. https://doi.org/10.1021/ic50093a033
- Cleare, M. J.; Hydes, H. C. In Metal Ions in Biological Systems; Sigel, H., Ed.; Marcel Dekker: New York, U.S. A., 1980, 11, 1.
- Cleare, M. J.; Hoeschele, J. D. Bioinorg. Chem. 1973,2, 187. https://doi.org/10.1016/S0006-3061(00)80249-5
- Graham, R. D.; Williams, D. R. J. Inorg. Nucl. Chem.1979, 41, 1245. https://doi.org/10.1016/0022-1902(79)80496-0
- Tang, W. X.; Dong, Y.; Qu, Y.; Dai, A. B. 分子科學與化學硏究 第2期 1985, 8, 167.
- Qu, Y.; Tang, W. X.; Dai, A. B. 藥學學報, Acta PharmaceuticaSinica 1986, 21, 586.
- Tatzanenko, L. V.; Moshkovsky. Y. S.; Zakharova, I.A.; Ponticelli, G.; Massacesi, M.; Devoto, G. Polyhedron,1985, 4(1), 159. https://doi.org/10.1016/S0277-5387(00)84237-0
- Zerner, M. C.; Bacom, A. D. Theor. Chem. Acta.1979, 53, 21. https://doi.org/10.1007/BF00547605
- Angerson, W. P.; Edwards, W. D.;Zerner, M. C. Inorg. Chem. 1986, 28, 2728.
- Angerson, W. P.; Cundar, T. R.; Zerner, M. C. Int. J. Quantum.Chem. 1991, 39, 31. https://doi.org/10.1002/qua.560390106
- Dewar, M. J. S., Zeobisch, E. G., Healy, E. F., and Stewart, J. J. P., J. Am. Chem. Soc. 1985, 107, 3902. https://doi.org/10.1021/ja00299a024
- Dewar, M. J. S., and Dieter, K. M., J. Am. Chem. Soc. 1986, 108, 8075 https://doi.org/10.1021/ja00285a033
- Stewart, J. J. P. Compt. Chem. 1989, 10, 209. https://doi.org/10.1002/jcc.540100208
- Stewart, J. J. P. Compt. Chem. 1989, 10, 221. https://doi.org/10.1002/jcc.540100209
- Burket, U.; Allinger, N. L. Molecular Mechanics, ACSMonography, Am. Chem. Soc. Washington D. C. 1981,177.
- Hyperchem Package, Copyright HyperCube, Inc. 419Phillip Street, Waterloo, Ontario N2L 3X2, Canada, 1996.