Molecular modeling of B-DNA site recognition by Ru intercalators: molecular shape selection

In this work, molecular modeling methods have been applied to the interaction characterization of polypyridyl transitional-metal complexes with the oligonucleotide (B-DNA fragment). In order to explore the factors governing the groove recognition and intercalative depth, we establish a simple and practical docking method (step-by-step docking operation) to obtain potential curves while making complexes inset into B-DNA along an assigned path. Energy values in the potential curve are obtained from energy minimization of binding geometries. Modeling results clearly show that the optimum binding conformation corresponding to the global minimum in the potential curve for each complex is found to correlate well with the experimental results. Our results also confirm that minor changes of the ligand structure can lead to profound influences on binding geometries, so the molecular shape of the complexes is a predominant factor in governing the binding mode. Moreover, we find that the vdW force and “water molecular effect” are strongly associated with molecular-shape selection in our model. These results complement and extend the knowledge of the nature of these complexes binding to B-DNA.