Elucidating the intercalation of methylated 1,10-phenanthroline with DNA: the important weight of the CH/H interactions and the selectivity of CH/π and CH/n interactions

Flat molecules like phenanthroline derivatives intercalate between base pairs of deoxyribonucleic acid and produce cytotoxic effects against tumoral cells. Elucidating the way of intercalation and its modulation on their efficiency by substitution still remains a challenging topic of research. In this work we analysed the intercalation via the major groove of methylated derivatives of phenanthroline, in different number and position, between guanine–cytosine base pairs. We studied our systems by using semi-empirical methods and density functional theory including dispersion corrections with the PM6-DH2 Hamiltonian and the B3LYP-D3 functional. We explored the geometry and electronic structure by means of the quantum theory of atoms in molecules and non-covalent interactions index analyses, whereas the interaction energy was estimated by means of two different approaches: one taking into account the results from the quantum theory of atoms in molecules analysis and the other based on the so-called energy decomposition analysis. The effect of solvation was also taken into consideration. Our studies show that CH/π and CH/n interactions by means of the –CH3 groups of methylated phen follow a clear pattern for any number of –CH3 groups and their position in the methylated phen ligand. That is, they try to produce the CH/π and CH/n interactions with the O and N heteroatoms of the base pairs and with the O atoms of the sugar and phosphate backbone. These findings suggest that the modulation of the intercalation of ligands that are able to form CH/π and CH/n weak interactions with the deoxyribonucleic acid is ruled not only by the number and position of the substitutions of the ligands but also by some key sites, which are the O and N atoms of the deoxyribonucleic acid in our analysed systems. It suggests some key and lock mechanism in which the interacting fragments fit like puzzle pieces in order to achieve the optimal interaction for the stabilization of the system. Interaction energies were calculated by using different approaches which converged to similar trends about the number and position of the –CH3 groups. The important weight of the CH/H interactions in the total interaction energy must be highlighted.


Introduction
Nowadays the intercalation of metal complexes in duplex DNA structures arises as an alternative and less aggressive option to the use of cisplatin in processes of chemotherapy to face cancer diseases. 1,2 The intercalation between DNA base pairs (bps) of any metal complex needs a at ligand in its structure and several studies have appeared in the literature considering this kind of structure as intercalators of DNA during the last few years. 1,3,4 Some theoretical studies have been carried out considering at ligands like 1,10-phenanthroline (phen) derivatives, 2,2 0 -bipyridine (bipy) and dipyridophenazine (dppz) by several authors [5][6][7][8][9] and also in our group. [10][11][12][13][14] On the other hand, the cytotoxic activity of metal complexes containing phen has been also proven experimentally by our collaborators. 15 In addition, the interest of phen derivatives for medical applications has been increased with its potential applications as alternative treatments against antibiotic resistant bacteria. 16 Because of the structural characteristics of duplex DNA, which is ruled by the hydrogen bonds of the bps and the stacking of these bps forming the double helix, the sugars of the nucleotides remain always in the same side, which produces some kind of narrow groove, which is called the minor groove. On the other hand, the opposite side of DNA without the sugars of the nucleotides is called the major groove because the width of such a groove has a higher size. Some experimental studies based on the determination of crystal structures pointed out that intercalation of metal complexes including at ligands between bps may take place via minor groove opposite to the classic opinion, which says that intercalation through the minor groove would produce more strain to the system. [17][18][19][20] Nevertheless, NMR structure analyses and several solution spectroscopic techniques seem to propose that the intercalation of these metal complexes including at ligands may occur via major groove according to the classical point of view and it is tempting to speculate that a delicate balance between major and minor groove intercalation is found in solution. 18,19 Due to the wide range of medical applications where the at intercalators can be used 21 an understanding at fundamental molecular level of the processes involved in the intercalation becomes very useful. For this reason, it is crucial to gain insight on the modulation of these processes of intercalation by means of different substituents and to stablish differences in the toxicity, which is very important in drug design to obtain efficient biological activity. Regard to the modulation by means of substituents, an inspiring experimental work was published by Brodie et al., 22 where a family of methylated phen derivatives of [Pt(en)(phen)]Cl 2 was studied to observe their trends in cytotoxicity against L1210 mouse leukemia. The authors pointed out that only the derivatives 5,6-Me 2 phen and 5-Mephen presented remarkable cytotoxicity, which suggested that substitution at position 5 of phen could have and important role when analysing the interactions with DNA at molecular level. In order to get insight on the relation between structure and cytotoxicity we already performed a previous computational study 13 in which we studied all the intercalating systems considered by Brodie et al. 22 in their experiments by means of model systems. Our results based on interaction energies veried our hypothesis on the important role of the position of the methylation and the systems 5-Mephen and specially 5,6-Me 2 phen had the most negative interaction energy, which was in agreement with our hypothesis on the importance of the position of the substitution.
In order to explore the whole scheme for the intercalation, we present in this work a complete analysis of the interactions between methylated phen derivatives and GC/CG bps. We chose GC/CG bps because of the observed better exothermicity for other intercalators when comparing to ATTA bps, 23 whereas we took into account the intercalation via major groove (see Scheme 1) because this orientation could be preferred in an aqueous environment, which is the general environment in biological systems. 18,19 We consider a complete geometrical study attending to the roll distortion 24,25 produced in the DNA bps aer the intercalation. Bond properties and non-covalent interactions like p-p, CH/p, CH/n and CH/H are also analysed with the topological tools of Quantum Theory of Atoms in Molecules (QTAIM) 26 and Non-Covalent Interaction (NCI) index. 27 Specially interest arise from the non-negligible role that we have found for the CH/H interactions. Actually, these CH/H interactions may have a similar weight as for the usual p-p interactions found in at small molecules intercalating between bps. On the other hand, we found that the CH/p and CH/n interactions coming from the methyl groups, which we called Me/p and Me/n, have preferences for specic O and N atoms localized in some key sites when interacting with duplex DNA. It suggests some key and lock mechanism where the intercalator try to t like a puzzle piece to achieve the optimal interaction with the intercalation pocket of the DNA to stabilize better the system. Interaction energies are discussed not only by means of the approach of Molins-Espinosa-Lecomte (EML) 28 following the results of the electronic density coming from the QTAIM analysis but also with the so-called Energy Decomposition Analysis (EDA). 29,30 We expect that this work will help to shed light on the behaviour of these important processes of intercalation, which are of biological and medical interest.

Computational details
In order to study the intercalation of methylated phen derivatives the so-called ring models 31 have been considered. We built these ring models from the 2ROU 32 crystal structure taken from the Protein Data Bank (PDB). Starting from this structure, the (1R-trans-antibenzo[c]phenanthrene) ligand between the bps of DNA was removed and the N2 of the guanine was saturated. From the whole DNA structure, we only kept the intercalation pocket (bps, phosphates and sugars) and the rest of the sequence was removed. The remaining O atoms that limited the sugar and phosphate backbone were saturated with hydrogen atoms. The intercalator, that is, phen and methylated derivatives, was placed manually between bps considering the maximum overlap and an equidistant arrangement. Finally, in order to keep the neutral charge of the whole system a Na + cation was put for each phosphate group present in the nal ring model.
The semiempirical Hamiltonian PM6-DH2, 33 which includes dispersion corrections, was employed to carry out full geometry optimizations of the ring models for each studied system. Due to the nature of the studied structures, stacked aromatic moieties with noncovalent interactions, such dispersion corrections are required. 33,34 This PM6-DH2 semi-empirical method was already tested with success optimizing the 1BNA structure from PDB, 35 and comparing the geometrical results with the original PDB structure. 13 The MOPAC2016 soware 36 was employed to carry out such optimizations without constrains.
In order to understand the nature of the interaction between the DNA structure and the intercalators at fundamental molecular level, from the optimized geometries, wave functions were generated at M06-2X/6-31+G(d,p)//PM6-DH2 level of theory for all the studied systems with the Gaussian soware 37 to carry out the analysis of the electron density (r) with QTAIM. 26 QTAIM describes accurately the concepts of atom and chemical structure by exploring the topology of r. From the QTAIM point of view, two atoms are bonded when they share a common interatomic surface (zero-ux surface) through which they can interact, and where there is a point (contained in the zero-ux surface) where r is a minimum in an specic direction but a maximum in a plane perpendicular to it. These points are known as bond critical points (BPCs) and the pair of gradient paths that connects the BCP with each nucleus is referred as atomic interaction line or bond path. From the QTAIM studies and r values we also estimated the energy of interaction for the studied systems by means of the EML approach. 28 Moreover, the non-covalent interactions, which play an important role in many processes related to biological systems, have been analysed with the NCI index approach developed by Johnson et al. 27 This approach provides a rich and 3D representation of noncovalent interactions (p-p and different kind of hydrogen bonds), with surfaces based on the peaks that appear in the reduced density gradient at low values of r. These surfaces are mapped according to the second Hessian matrix eigenvalue, negative values (stabilizing interactions) are depicted in blue and pale green, while positive values (destabilizing interactions) are represented in yellow and red. This approach provides useful information in systems where p-p interactions are present between aromatic moieties (DNA bps and at intercalators). The AIMAll soware 38 has been used to perform both, the QTAIM topologies and NCI index isosurfaces.
The interaction energy (DE int ) between the intercalator and the DNA structure has been studied with the EDA that is based on the Morokuma-type decomposition method. 29,30 Such methodology allows the calculation of the contributions of the different terms of the interaction energy between two fragments: (1) the intercalator and (2) the rest of DNA structure (bps, sugars and phosphates). With this analysis DE int is split into different contributions: orbital energy related to polarization and charge transfer terms (DE orb ), Pauli repulsion contribution related to the destabilizing interactions between occupied orbitals (DE Pauli ), electrostatic contribution related to the classical electrostatic interaction between the unperturbed charge distributions of the rigid fragments (DE elstat ) and dispersion energy associated to van der Waals forces (DE disp ). In order to perform the EDA, ADF soware [39][40][41] has been employed with single-point calculations on the PM6-DH2 optimized geometries, using the uncontracted polarized triple-z basis set of Slater-type orbitals (TZP) and the B3LYP-D3 functional with the Grimme's D3 correction to include dispersion forces. [42][43][44][45] The solvent effects have been included with the COSMO solvent approach 46 because the relevance of solvation in these biological systems has been reported previously. 12 The Basis Set Superposition Error (BSSE) was not considered due to the fact that in previous works, 10,11 where the intercalation of phen derivatives between DNA bps was already studied, the error was only 8-12%.

Results and discussion
First of all the geometrical results will be presented. In Fig. 1 the coplanar arrangements for the optimized structures of the corresponding methylated derivatives of phen are shown along with the values obtained for the roll distortion 24,25 (see also ESI †) produced in the process of intercalation between the bps. It is observed that the roll distortion presents relative low values for all the intercalators between 1.6 and 6.2 . Surprisingly, for the most substituted intercalator, 3,4,7,8-Me 4 phen, the intercalation via major groove presents a lower value in comparison to the intercalation via minor groove checked in a previous work. 13 This difference can be explained because of the location of the -CH 3 groups in positions 7 and 8, which remain far away from the sugar and phosphate backbone in the case of the intercalation via major groove. On the contrary, in the case of the intercalation via minor groove the -CH 3 groups in positions 7 and 8 were located close to the sugar and phosphate backbone and produced destabilizing interactions with it. Another geometrical characteristic to take into account is the relative position of the aromatic moieties and how the molecular planes are facing each other. The main conclusion that arise from the geometrical arrangement is that the -CH 3 groups tend to be located near and between the O and N atoms of the GC/CG bps. This behaviour allows the non-substituted phen, 5-Mephen and 5,6-Me 2 phen to present the best overlapped arrangement between the ligand and DNA bps planes. On the other hand, 4-Mephen, 4,7-Me 2 phen and 3,4,7,8-Me4phen structures present the most distorted arrangements between the aromatic moieties, leaving part of the aromatic plane of the ligand outside of the intercalation pocked and in consequence, showing less overlap. Such result will be corroborated below with the QTAIM and NCI analyses, where interactions between the -CH 3 groups and O and N atoms of the bps were found.
To gain insight into the role played by the electronic density in the intercalation of phen derivatives in DNA structure NCI analysis is shown, such analysis provides a 3D representation of non-covalent interactions. Due to the presence of the weak interactions between the intercalator and the DNA (p-p stacking and several kinds of hydrogen bonds), the NCI analysis will give us a complete view to describe the bonding scheme between the intercalator and DNA structure towards the noncovalent interactions. In Fig. 2 the NCI isosurfaces are plotted. Blue and pale green indicate negative values, which corresponds to stabilizing interactions, while yellow and red indicate positive values and therefore destabilizing interactions. In order to show a clearer and cleaner picture, the NCI isosurfaces have been computed only for regions surrounding the intercalator. In this way the weak interactions between phosphates and sugars in the backbone and the isosurfaces corresponding to destabilizing interactions in the inner regions of sugar rings, that will overload the picture, are not shown. Moreover, the connectivity towards the bond paths according to the QTAIM methodology is also plotted, showing a good agreement between NCI and QTAIM analyses as expected. The characteristics of the corresponding BCPs will be described below.
It is observed that for all intercalators a large green region between the aromatic moieties of DNA structure and intercalator is predominant along with a more stable negative value for the NCI index in the surrounding zones of the BCPs (depicted in cyan). This is a common result in systems where pp stacking is produced. 13,14,47,48 For the non-substituted phen it is observed that in addition to the large NCI isosurface corresponding to the interplanar region, the presence of CH/H interactions is shown with H belonging to the sugar structures for the H atoms located in positions 4 and 8 of the intercalator. Similar kind of interactions have been already depicted previously with different computational approaches. [49][50][51][52] Moreover, the H atom in position 7 interacts with the O atom of sugar moiety. All these interactions present high and negative value for the NCI index, which indicates a strong and stabilizing interaction and thus, the important role of this kind of interactions in the studied systems.
When -CH 3 groups are present in phen structure, as we have observed in geometrical results depicted above, -CH 3 groups tend to be located in the proximity of the O and N atoms belonging to GC/CG bps. It is observed in Fig. 2 that the -CH 3 groups yield interactions with a considerable negative value of the NCI index with the O atom of the cytosine and the N atom of the guanine. Moreover, another interaction with the O atom belonging to the sugar appears. This is the preferred arrangement for the -CH 3 groups and it is present in all kind of substitutions at least for one -CH 3 group. According to this, the most favourable scenario corresponds to the 5,6-Me 2 phen where both -CH 3 groups present interactions with N and O atoms of the bps and at the same time with the O atoms belonging to the sugars of the sugar and phosphate backbone. All the mentioned interactions tend to present considerable negative values and they show cyan or dark cyan isosurfaces (see ESI †). In the following lines we are going to focus on the discussion on the topological analysis of r provided by the QTAIM analysis. The value of r at the BCPs indicates the strength of the bond formed between two atoms, while the electron-energy density (E d ) describes the stability of the chemical bond. In this kind of systems where the weak interactions prevail between the intercalator and the DNA structure, the values of E d are close to zero as characteristic of weak interactions. We observe the general agreement between the QTAIM and the NCI analyses where the regions of negative values in the NCI isosurfaces are surrounding the BCPs. In addition, the geometrical results, where the -CH 3 groups tend to be located close to the heteroatoms of bps, can be explained with the interactions depicted with the analysis of the r. Fig. 3 shows the topological analysis of r for all the systems intercalating via major groove. In order to show a clearer picture, only the interactions between the intercalator and the DNA structure are presented (interactions between sugar and phosphates are omitted). For all the structures several p-p interactions were found (blue boxes), and to avoid overloaded representations these p-p interactions have been only described for the non-methylated phen (see ESI † for the whole scheme). Aside from the p-p BCPs, interactions between the DNA structure and H atoms of phen also appear: CH/n sugar , CH/p bps and CH/H sugar , which are presented in green boxes.
Aer the implementation of -CH 3 groups in the phen structure we can see that several interactions appear due to the presence of -CH 3 groups in the surrounding environment of the pocket formed by the DNA bps. It is observed for 4-Mephen that four interactions appear because of the presence of the -CH 3 group. Three of these interactions are produced with the heteroatoms belonging to the plane of bps with values of 0.005 and 0.008 a.u. for r, and one interaction with an O atom of the nearby sugar with a high value of r (0.014 a.u.). We have previously observed in the geometrical results that the preferred position of the -CH 3 is located in the proximity of the O and N atoms of the interacting bps. Thus, the localization of the O and N atoms rule the global arrangement of phen derivatives intercalated between the bps. In this case where -CH 3 is in position 4, the intercalator plane is slightly turned in order to satisfy such position and interactions for the -CH 3 . A very similar scheme is presented for 5-Mephen, showing two interactions with heteroatoms of the bps and one with the O atom of the nearby sugar, being again the interaction with the O atom of the sugar the most strong with a high value of r in such BCP (0.011 a.u.). A different behaviour is presented for 4,7-Me 2 phen and 3,4,7,8-Me 4 phen. When different -CH 3 groups are present in positions far away from each other. That is, one of the -CH 3 groups is located close to the O and N atoms of the interacting bps, the preferred area, while the other -CH 3 groups are le out, which yields most distorted arrangements between the intercalator and the DNA structure, as we have observed previously in the geometrical analysis. For these structures the -CH 3 groups that are le out can also interact with the sugar and phosphate backbone and bps at the same time. Finally, it is noteworthy the bonding scheme presented by 5,6-Me 2 phen. In this case both -CH 3 groups are located in the centre of the structure of the bps, the preferred area, close to N and O. In this position both -CH 3 groups may form interactions with the heteroatoms of the bps and at the same time with the O atoms of the sugars. In this case both -CH 3 groups yield eight interactions, ve with the heteroatoms of the bps and three with the sugars, being the interactions formed with the O atoms of the sugar those that present higher values of r (0.015 a.u.). Moreover, another characteristic that should be highlighted is that when the -CH 3 groups are located at positions 5 and 6, the whole aromatic plane of phen is completely overlapped with the aromatic environment of the bps. It can be appreciated in the NCI results where for this system large isosurfaces between the aromatic moieties are observed. All these contributions place the 5,6-Me 2 phen as the preferred phen methylated derivative for an effective intercalation between CG/GC bps via major groove. Table 1 summarizes all the BCPs presented between the corresponding intercalator and the DNA structure. As expected, the increasing number of -CH 3 groups yields more BCPs, which increase the number of weak interactions. The increasing is remarkable for 5,6-Me 2 phen and 3,4,7,8-Me 4 phen due to the interactions of the -CH 3 groups with the sugars and the p environment of the bps.
Attending to the results obtained in Table 1 and Fig. 3 the question that arises is of how much the different kind of interactions contribute to the global stability of the intercalation. Thus, based on the properties of r at the BCPs, the interaction energy (DE int ) between two fragments can be calculated with the EML approach 28 (eqn (1)). Such approach has been previously employed in DNA systems where the role of the backbone, in the interaction of minor groove ligand binding, has been described. 53 The EML approximation was developed to consider H bonds only with O atoms, while when other kind of interactions are present, Díaz-Gómez et al. 53 implemented a correction for the BCPs different from the H/O hydrogen bonds (eqn (2)). Attending to such approaches Fig. 4 and Table  S1 of ESI † summarise all the stabilization contributions depending on the different weak interactions present between the phen derivatives and the GC/CG bps along with the total stabilization for each studied system according to the corrected EML approach.
where DE wi is the energy associated to any specic kind of weak interaction (XH/n, CH/n, CH/p, CH/H, etc.), whereas V(r cp ) is the potential energy density at the same BCP associated to such specic weak interaction localized in the studied system. This V(r cp ) is the average effective potential eld experienced by r and is proportional to the strength of this kind of weak interaction, DE wi . The summation of all the values of DE wi for each kind of weak interaction leads to an estimation of the total interaction energy of the studied system, DE int .
It is observed in Fig. 4 that the general trend is the increasing of the strengthening of the interaction with more negative values of DE int when the number of methyl groups increases. For all the structures, the main stabilization in general is due to all the BCPs corresponding to p-p interactions. However, such p-p contribution is very similar for all the studied systems (from À7.3 to À9.9, which is a difference of only 2.6 kcal mol À1 ). Therefore, the differences in the total interaction energy, DE int , will be given by the other weak interactions localized in the studied systems. Moreover, the values obtained for CH/H interactions must be also highlighted since they present considerable high values (from À5.9 to À9.9 kcal mol À1 ). The addition of such CH/H contributions lead to estimations of the total interaction energies, DE int , which are very similar to the DE int values obtained with the EDA discussed below. It indicates not only the correct consideration of these non-conventional interactions in the extrapolation of the DE int with the EML approach but also the relevance of these CH/H interactions that have been reported previously in the bibliography. [49][50][51][52] This behaviour can be understood with the high values for the r obtained from the QTAIM analysis in such interactions and in addition with the high negative values for the isosurfaces in the NCI analysis.
Another considerable contribution to the strengthening of the interaction energy and more negative values of DE int is obtained with the Me-p interactions, being À5.0, À4.5 and À8.0 kcal mol À1 for 4,7-Me 2 phen, 5,6-Me 2 phen and 3,4,7,8-Me 4 phen, respectively. Due to this contribution these three derivatives show a total DE int of À32.7, À28.4 and À31.1 kcal mol À1 , respectively, being the most preferred intercalators for the intercalation between GC/CG bps through the major groove. Considering the different kinds of interactions with different nature it can be assumed that the values for the constants in EML approach are still too general to describe properly the interaction energies and in addition, such approach is very local and only focus on the properties of the BCPs, whereas it does not take into account the whole effect of the fragments and the feasible cooperative effects that could lead to a more strengthened interaction and therefore to more negative interaction energies.
At this point we have to say that the EDA provides an accepted and accurate way to determine the interaction energy (DE int ) between fragments that can be decomposed into different contributions: electrostatic (DE elstat ), charge transfer and polarization terms in the orbital contribution (DE orb ), dispersion (DE disp ), and the destabilizing interaction between occupied orbitals (DE Pauli ): Fig . 5 shows the EDA results for the intercalation of phen and the studied methylated phen derivatives when intercalating between GC/CG bps via major groove. In general, the attractive terms are slightly increased with the presence of -CH 3 groups. That is, the DE disp is gradually increased with the presence of -CH 3 groups from À50.4 to À59.4 kcal mol À1 . This is also the general trend for the repulsive Pauli contribution. In the graph of Fig. 5 it is observed that such repulsion requires the small contributions of the electrostatic and orbital terms to stabilize the whole systems. Considering this fact and attending to the values of the electrostatic forces, the 5,6-Me 2 phen and 3,4,7,8-Me 4 phen are the most stabilized intercalators due to the high value of the electrostatic contribution À39.1 and À37.8 kcal mol À1 , respectively.
From the values of the EDA it is observed that the implementation of -CH 3 groups stabilized the intercalation, but the most effective intercalation in terms of number of -CH 3 groups vs. DE int corresponds to 5,6-Me 2 phen. With only two -CH 3 groups this compound reaches a very similar value of the DE int when comparing to the DE int value of the 3,4,7,8-Me 4 phen that presents different kinds of weak interactions in the four -CH 3 groups. Summarising, we can suggest that the most effective positions for the methylation are 5 and 6. This is in agreement with the results obtained from the topology of r, where we  showed that in such positions -CH 3 groups yield very effective weak interactions with the N and O atoms of the upper and down bps, and at the same time with O atoms of the sugars present in the sugar and phosphate backbone. In the same way the stabilization presented by 5-Mephen, À31.9 kcal mol À1 , can be understood due to the position of the -CH 3 that forms 3 effective weak interactions with O atoms of the bps and the sugar.
To gain more insight on the role of the solvent in the process of the intercalation, the desolvation penalty has been computed and added to the total DE int of the EDA and the obtained values are presented in Table 2.
Looking at the DE Aq energies, we observe that aer including the solvent penalty, DE Solv , to the interaction energy, DE int , 5-Mephen, 5,6-Me 2 phen and 3,4,7,8-Me 4 phen still remain as the intercalators that lead to more negative DE Aq energies being À20.7, À19.9 and À20.2 kcal mol À1 , respectively, and therefore achieving better intercalation via major groove with GC/CG bps.

Conclusions
Our present work offers important and detailed comprehension to understand the chemical reasons of the effect of methylation, in number and position, on the process of the intercalation of phen and methylated derivatives and how the strength of the interaction with DNA may be modulated. The number of -CH 3 groups present in the phen aromatic rings favour the process of intercalation between GC/CG bps through the major groove, and in consequence its cytotoxicity. Nevertheless, we may affirm that the position of -CH 3 groups is more important than the number. In the intercalation of methylated phen derivatives between GC/CG bps via major groove the number of -CH 3 groups does not increase signicantly the steric repulsion. However, from the topological analysis of the r it can be observed that the 5 and 6 positions yields effective weak interactions with O and N heteroatoms of the environment, O and N belonging to the bps and also O atoms belonging to sugars of the sugar and phosphate backbone. At the same time the methylation in 5 and 6 positions favours a better overlap and coplanarity between phen and bps aromatic moieties yielding a more effective p-p stacking and thus to a better inclusion and stabilization of the ligand between bps of DNA. Thus, we may affirm that the process of the intercalation of methylated phen derivatives with DNA will follow the key and lock mechanism in which each fragment, intercalator and intercalation pocket between bps, t like puzzle pieces in order to achieve the most effective interaction. On the other hand, the importance of CH/ H interactions are highlighted when calculating the interaction energy with the EML approach in which the weight of the electronic density corresponding to these CH/H interactions is very important in the total interaction energy. We compared the DE int obtained from the EML approximation with the DE int coming from the EDA, which is another different procedure to obtain the DE int , and the values and trends converge and are very similar, which corroborates the important role that CH/H interactions play in this kind of processes of biological and medical interest.

Conflicts of interest
There are no conicts to declare.