Geometry of intercalation of psoralens in DNA approached by molecular mechanics

The results of molecular mechanical calculations on intercalation complexes of 3-carbethoxypsoralen, 5-methoxypsoralen, 8-methoxypsoralen, 7-methylpyrido[3,4-c]psoralen (MepyPs) and 7-methylpyrido[4,3-c]psoralen (2N-MePyPs) with the double stranded duodecanucleotide d(CGCGATATCGCG)2 are presented. In the energy-minimized structures, the psoralens are intercalated with their plane orthogonal to the helix axis. Stacking interactions between the furan ring of the psoralen and the adjacent bases are maximized in most derivatives studied, whereas the effect of the various substituents of the psoralen ring is to specifically push part of the molecule towards either the minor or the major groove, preventing a symmetrical intercalation (with respect to the two strands of the DNA). The relative position of the psoralen ring and of the adjacent thymine foreshadows the formation of furan-side monoadducts in 3-CPs, MePyPs and 2N-MePyPs, whereas the formation of a pyrone-side monoadduct appears as geometrically more favourable in 5-MOP and both furan- and pyrone-side monoadducts can be geometrically envisaged in 8-MOP. A good correlation therefore exists between the more or less favourable equilibrium geometries and the experimentally observed photoreactions. The present study is the first attempt to characterize the geometrical parameters as part of a complex set of geometrical, dynamical and excited state parameters governing the overall DNA-psoralen photoreaction.