Bending of DNA upon binding of ecteinascidin 743 and phthalascidin 650 studied by unrestrained molecular dynamics simulations

被引:37
作者
García-Nieto, R
Manzanares, I
Cuevas, C
Gago, F [1 ]
机构
[1] Univ Alcala de Henares, Dept Farmacol, E-28871 Alcala De Henares, Madrid, Spain
[2] Pharma Mar SA, Madrid 28760, Spain
关键词
D O I
10.1021/ja000964q
中图分类号
O6 [化学];
学科分类号
0703 ;
摘要
Recognition of DNA sequence information by the natural antitumor agent ecteinascidin 743 (ET743) has been proposed to operate through a direct readout mechanism involving specific hydrogen bonding interactions between the drug and the DNA minor groove prior to covalent alkylation of a guanine base. 5'-AGC and 5'-CGG are examples of high-reactivity target triplets whereas 5'-CGA is a low-reactivity sequence. Molecular dynamics computer simulations were first used to explore the stability and behavior of the pre-covalent complex between ET743 and a DNA nonamer containing the 5'-AGC binding site in the central region. The pre-covalent complex was stable, and some unreported distinctive features were observed that may not be amenable to direct experimental verification. A similar simulation with ET743 bound to another nonamer containing a central 5'-CGA triplet did not result in a stable association supporting the proposed role of a hydrogen bonding network in the stabilization of these complexes. The covalent complexes between ET743 and the nonamers containing the 5'-AGC and 5'-CGG target sites were then simulated. In each case the drug displayed the predicted binding mode and gave rise to a widening of the minor groove. In addition, we show that as a consequence of ET743 binding to the target sequences, positive roll is introduced that results in smooth bending of the helix toward the major groove, in agreement with results from gel electrophoresis experiments. Similar results were obtained with the synthetic compound phthalascidin, which presents a biological profile almost indistinguishable from that of ET743. The local bending elements in AGC and CGG were found to be different, and the distinct behavior of these sequences in the absence of bound drugs was in consonance with their intrinsic bending propensities.
引用
收藏
页码:7172 / 7182
页数:11
相关论文
共 56 条
[1]  
AQVIST J, 1990, J PHYS CHEM-US, V94, P8021, DOI 10.1021/j100384a009
[2]   OPTIMIZED PARAMETERS FOR A-DNA AND B-DNA [J].
ARNOTT, S ;
HUKINS, DWL .
BIOCHEMICAL AND BIOPHYSICAL RESEARCH COMMUNICATIONS, 1972, 47 (06) :1504-&
[3]   A WELL-BEHAVED ELECTROSTATIC POTENTIAL BASED METHOD USING CHARGE RESTRAINTS FOR DERIVING ATOMIC CHARGES - THE RESP MODEL [J].
BAYLY, CI ;
CIEPLAK, P ;
CORNELL, WD ;
KOLLMAN, PA .
JOURNAL OF PHYSICAL CHEMISTRY, 1993, 97 (40) :10269-10280
[4]   Flexibility of the B-DNA backbone: effects of local and neighbouring sequences on pyrimidine-purine steps [J].
Bertrand, HO ;
Ha-Duong, T ;
Fermandjian, S ;
Hartmann, B .
NUCLEIC ACIDS RESEARCH, 1998, 26 (05) :1261-1267
[5]   ATOMIC CHARGES DERIVED FROM SEMIEMPIRICAL METHODS [J].
BESLER, BH ;
MERZ, KM ;
KOLLMAN, PA .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1990, 11 (04) :431-439
[6]   MOLECULAR-DYNAMICS STUDIES OF DNA [J].
BEVERIDGE, DL ;
RAVISHANKER, G .
CURRENT OPINION IN STRUCTURAL BIOLOGY, 1994, 4 (02) :246-255
[7]  
Bonfanti M, 1999, ANTI-CANCER DRUG DES, V14, P179
[8]   A modified version of the Cornell et al. force field with improved sugar pucker phases and helical repeat [J].
Cheatham, TE ;
Cieplak, P ;
Kollman, PA .
JOURNAL OF BIOMOLECULAR STRUCTURE & DYNAMICS, 1999, 16 (04) :845-862
[9]  
Cheatham TE, 1996, J MOL BIOL, V259, P434
[10]   Molecular dynamics simulations highlight the structural differences among DNA:DNA, RNA:RNA, and DNA:RNA hybrid duplexes [J].
Cheatham, TE ;
Kollman, PA .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1997, 119 (21) :4805-4825