1. Please check and comment entries here.
Table of Contents

    Topic review

    Duplex DNA and G-Quadruplexes

    View times: 5
    Submitted by: Adrià Gil-Mestres

    Definition

    Biomedicine may be considered the cornerstone of modern health care.  It includes key enabled technologies such as molecular biology, biotechnology, nanobiotechnology, biological engineering, etc., and concerns a wide range of scientific and technological approaches that range from the understanding of molecular interactions to the study of gene therapy.

    1. Introduction

    Biomedicine may be considered the cornerstone of modern health care. [1] It includes key enabled technologies such as molecular biology, biotechnology, nanobiotechnology, biological engineering, etc., and concerns a wide range of scientific and technological approaches that range from the understanding of molecular interactions to the study of gene therapy. According to the most recent statistics in Europe, there were 1.93 million deaths caused by cancer in Europe during 2018 (~36.3% of the total number of deaths), which means that cancer is still an important topic for biomedical research and that improvements in chemotherapy are still necessary. On the other hand, the emerging problem of bacterial antibiotic resistance (BAR) has caused 33,000 deaths per year in the EU (https://ec.europa.eu/health/antimicrobial-resistance/eu-action-on-antimicrobial-resistance_en accessed on 1 August 2021) and has become a current hot topic of research. Thus, research developed within these topics for the next years becomes important, timely, and within the priorities of the world concerning health, healthy ageing, and wellbeing. Cisplatin is the reference drug in chemotherapy treatments against cancer. Ref. [2] However, this drug may cause serious side effects due to the damage of normal tissues [3]. Research aimed at the development of more efficient alternatives is valuable.
    The use of flat ligands such as 1,10-phenanthroline (phen), both in their isolated form or included in metal complexes, was devised some years ago as an alternative to cisplatin for chemotherapy treatments [4][5][6]. Their use as an innovative method to fight against BAR has been considered more recently. Refs. [7][8][9][10] They intercalate between DNA base pairs (bps) and can either inhibit the replication of DNA or cleave the DNA chain causing the death of cancer cells or bacteria. Nevertheless, some competition between intercalation and groove binding modes of interaction between these flat ligands and DNA has been proposed in the literature. Refs. [11][12][13][14][15] Indeed, whereas the groove binding occurs very fast (i.e., tenths of a millisecond), the intercalation mode of interaction takes more time to occur (i.e., in the range of few milliseconds). Ref. [12] Moreover, the intercalation mode is usually more related to cytotoxic effects. Since the cytotoxic effect of any intercalator depends on the time of residence of the drug between bps, [16] the design of any efficient drug should aim at an increased drug-DNA interaction to stabilize the intercalated state, but it should do so at less stable groove binding states to make the kinetics faster. Such modulation for the binding sites, and thus for the cytotoxicity, can be achieved by substitution of phen in number and position. For this reason, the comprehension and rationalization of the interactions between flat ligands and duplex DNA (dDNA) becomes crucial. They will give us information about how to modulate the interactions in both modes by substitution in the flat ligands and how to optimize the drug design. In this sense, the effect of substitution in phen in the intercalation process is still not clear. On the other hand, the antibacterial efficiency of phen derivatives is higher when they are coordinated to any metal than when the ligand is alone [10]. However, no satisfactory explanation has been given for it. Thus, the influence of ancillary ligands in the intercalation process is an interesting topic for study. Our contribution to state-of-the-art studies has been the analysis of the nature of the interaction and the investigation concerning how the substitution of flat ligands, such as phen, in number and position (with different kind of functional groups as -CH3, -OH, -NH2, =O, -Ph, -Cl, -COOH, etc.) favors the intercalation or the groove binding competitive mode of interaction. Moreover, the inclusion of metal atoms has been also analyzed.
    Another alternative to overcome the problem with cisplatin is based on the stabilization of G-quadruplexes (GQ), which are alternative non-canonical quadruple-stranded helical DNA structures found in guanine rich sequences of DNA. The formation of GQ gives other singularity to the DNA, which may lead to more selective interactions. Moreover, by avoiding the abovementioned side effects, this new target could replace therapies based on cisplatin. Indeed, the formation and stabilization of GQ were shown to decrease the activity of telomerase [17], which is the enzyme responsible for the elongation of telomeres, a phenomenon that prevents cell apoptosis. Since high telomerase activity is involved in 85% of cancers [18], it is recognized as a potential cancer specific target. Thus, inhibition of telomerase becomes a key method for stopping tumoral cell growth. The presence of GQ in promoters also represents a subject of study. Ref. [19] In this case, the presence of stable GQ in oncogene promoters can alter the expression of the gen reducing some key processes in the growth of tumor cells. Ref. [20] Thus, stabilization of GQ may also be an innovative strategy to fight against BAR since the stabilization of GQ in bacteria by small molecules may also inhibit the expression of genes responsible for BAR. Ref. [21] Several organic ligands, either alone or included in metal complexes, were reported in the literature as GQ stabilizers, inhibiting telomerase activity or disrupting the transcriptional activity of some oncogenes. Ref. [22] There are three main sites in the GQ where these stabilizing small molecules may interact: end-stacking, grooves, and loops (with the first being the most common). In order to favor such interactions, small planar molecules may induce end-stacking binding, whereas the interactions with side loops and grooves are enhanced by the presence of side chains, which are positively charged or have affinity for protons that are attached to the planar aromatic cores. These side chains participate in electrostatic interactions with the negatively charged DNA phosphate backbones. Ref. [23] Moreover, the same organic ligands bound to different metal centers may still generate different GQ affinities and interactions.
    On the other hand, and trying to devise other strategies, it was only during the early-mid 2000s that systematic studies on the application of nanostructures based on polyoxometalates (POMs) for cancer treatment became a booming domain. Refs. [24][25][26] Recent advances in the use of POMs for medical applications and their antitumor activity have been reviewed by Bijelic et al. [27] In addition, Bijelic et al. and Kortz et al. also reviewed the latest developments on the use of POMs against BAR. Refs. [28][29] According to these reviews, the proposed modes of antitumoral action of POMs involve several possibilities: (1) the activation of cell death pathways; (2) inhibition of angiogenesis; (3) interaction with proteins; or (4) DNA interaction, among other mechanisms. In the case of POMs action against bacteria the proposed mechanisms considers: (a) inhibition of both PBP2a and β-lactamases; (b) targeting P-type ATPases; (c) impairment of the bacterial electron-transport chain (respiratory system); (d) POM-mediated increase of the reactive oxygen species level via oxidation; (e) interaction with important membrane-anchored proteins and enzymes; (f) disruption of the bacterial cytoskeleton dynamics by POM-interactions with cytoskeletal elements; (g) disruption of the bacterial cell wall leading to leakage of intracellular substances; or (h) interaction with cytoplasmic elements of proteins that are anion-sensitive like nucleic acid-binding proteins. In our studies, among all the processes described in the reviews of Bijelic et al. [27][28], we focused on the interaction of POMs with DNA to promote its phosphoester hydrolysis as artificial phosphoesterases. The seminal works found in the bibliography on this topic are the experimental studies of Parac-Vogt et al. based on [Mo7O24]6− ({Mo7}) species, [30][31][32][33] whereas other experimental works of Abrantes et al. [34][35] were based on the MoO2X2L species.

    2. Computational Techniques, Methods and Tools Used in Our Studies

    We have used different computational approaches to carry out the different calculations for the three topics developed in these studies. In the lines of research regarding the interaction of small molecules based on phen derivatives and Mo complexes, including phen with dDNA and GQ, we used several models along with different levels of calculation. That is, in the case of sandwich models [36], we used the M06-2X/6-31+G(d,p) level of calculation [37][38][39] with Gaussian09 in order to study the intercalation of phen derivatives with dDNA. Ref. [40] The M06-2X functional is recommended to study systems in which weak interactions such as dispersion are important to explain their behavior (as the intercalation of small molecules between bps of DNA). In the case of the ring models [36] including sugar, phosphates, and Na+ counterions, and for the study of the groove binding interactions with the d(GTCGAC)2 hexamer, we used the semi-empirical Hamiltonian PM6-DH2 [41] including dispersion effects with the MOPAC software. Ref. [42] In the case of the ring model 2 ps semi-empirical MD, simulations with the PM6-DH2 Hamiltonian were also performed with a time step of 1fs. Explicit water molecules were considered by means of the TIP3PBOX solvent model for the simulations by using a rectangular box with edges no closer than 5 Å to any atom of the solute. Finally, in the case of the studies on the interaction of the Mo[(η3-C3H5)Br(CO)2(phen)] complex with dDNA and GQ, we took into account the d(AGACGTCT)2 octamer for the dDNA, coming from the 1n37 PDB structure, and the GQ coming from the 2jwq PDB structure to study the interaction of the Mo[(η3-C3H5)Br(CO)2(phen)] with this non-canonical structure. These structures were studied at the LMKLL/DZDP level, [43][44] which includes van der Waals corrections. The core electrons were substituted by norm-conserving pseudopotentials. Refs. [45][46] Such LS-DFT computations with ~500 atoms for the dDNA and ~1000 atoms for the GQ were carried out with the Spanish Initiative for Electronic Simulations of Thousands of Atoms (SIESTA) method and associated software. Ref. [47] The cut-off radii for the atomic orbitals of each element were obtained for an energy shift [47] of 30 meV. The tolerances used for the optimizations were 10−5 eV for the energetics, whereas the tolerance for the forces was 0.02 eV/Å for the dDNA and 0.1 eV/Å for the GQ. In order to gain deeper insight into the interaction between the small molecules and DNA substrates, we performed the Energy Decomposition Analysis (EDA). Refs. [48][49] To carry out the EDA, we mainly used the B3LYP-D3/TZP level of calculation [50][51][52][53] since the B3LYP-D3 functional includes an explicit Grimme’s D3 correction for dispersion and, therefore, an additional ΔEdisp term appears in the EDA for this functional. For this reason, we thought that the discussion was simpler to visualize better the trends of the intercalated systems. It must be said that the M06-2X/TZP and M06-L/TZP levels were also checked for the EDA. In any case, the three functionals led to results comparable to the MP2/6-31G* (0.25) level of theory already used by Řeha et al. [54] which gave results comparable to the benchmark CCSD(T) data. Ref. [55] These EDA were carried out with the ADF software. Refs. [56][57][58] Another way to gain deeper insight into the interaction between the studied small molecules and DNA structures is the topological analysis of the structures. Two kinds of approaches were used, namely the classical QTAIM developed by Bader et al. [59] and the most recent approach developed by Johnson et al. [60] based on the NCI. The latter provides a rich, 3D representation for the non-covalent interactions with surfaces based on the peaks that appear in the reduced density gradient at low values of ρ. Such isosurfaces are mapped according to values of the sign of the second Hessian eigenvalue, (λ2), and while negative values (i.e., stabilizing interactions) are depicted in blue and pale green, positive values (i.e., destabilizing interactions) are represented in yellow and red. QTAIM and NCI computations were performed with the AIM2000 [61] and AIMALL [62] software with the wave functions generated at M06-2X/6-31+G(d,p) level with Gaussian09 in all cases with the exception of the interactions of [Pt(en)(phen)]2+ derivatives via groove binding with dDNA in which the wave function was obtained at B3LYP/6-31G(d,p) level with Gaussian09. Finally, for some of these systems in which small molecules are interacting with DNA, we studied the polarization/charge transfer for the different modes of interaction with several charge schemes arising from different approaches (Mulliken, [63] APT, [64] Hirshfeld [65] and NPA/NBO [66][67]).
    In the study of the activity of Mo-oxo species as promoters and catalysts for the hydrolysis of the phosphoester bond, we carried out the DFT calculations by using ADF [56][57][58] and the Gaussian 16 Revision A.03 software. Ref. [68] In the case of ADF, geometry optimizations were carried out with the BP86 functional, which uses the Vosko-Wilk-Nusair exchange-correlation potential [69] with the generalised gradient approximation exchange correction reported by the Becke (1998) exchange functional [70] and Perdew (1986) correlation correction [71], as well as Grimme dispersion corrections (BP86-D3). Ref. [53] Relativistic effects were treated with the zero order regular approximation (ZORA) Hamiltonian. Refs. [72][73] The frozen core approximation and triple-ζ Slater-type orbitals (STO) were used to describe the valence shells of C and N (2s and 2p). One polarization function was added to C, N, O, and Mo (single-ζ, 3d, 4f). Triple-ζ STOs were used to describe the valence shells of H (1s) augmented with one polarization function (single-ζ, 2s, 2p). Solvent effects were included with the COSMO [74] and standard parameters (Water, ε = 78.39). Analytical frequencies were calculated to characterise the obtained stationary points and calculate the Gibbs free energies (standard state T = 298.15 K, p = 1 atm). Transition states were followed after a fractional displacement of the imaginary vibrational mode to both the reactant(s) and product(s). On the other hand, Gaussian calculations were carried out with the B3LYP hybrid functional. Refs. [50][51][52] The LANL2DZ effective core potential with the associated double zeta basis set supplemented with f polarisation functions was used for Mo atoms [75][76][77][78][79] and the 6-31+G(d,p) basis set [39] for the rest of the atoms. Dispersion effects were included by using the third version of Grimme dispersion with the Becke–Johnson damping approach. Ref. [80] Stationary points were characterised by means of frequency calculations, and the intrinsic reaction coordinate was followed in order to obtain the geometries of the reactants and products followed by unconstrained optimisations. Refs. [81][82] For some mechanistically relevant stationary points the energies were refined with the 6-311++G(3df,2p) basis set. Ref. [83] For these calculations with Gaussian 16 Revision A.03 the Polarisable Continuum Model (PCM) was used to take into account the solvent in an implicit way [84][85].

    The entry is from 10.3390/molecules26164737

    References

    1. Scheffler, R.W.; Strasser, B.J. Biomedical Sciences, History and Sociology of. In International Encyclopedia of the Social & Behavioral Sciences, 2nd ed.; Elsevier: Orlando, FL, USA, 2015; pp. 663–669.
    2. Harrap, K.R. Preclinical studies identifying carboplatin as a viable cisplatin alternative. Cancer Treat. Rev. 1985, 12, 21–33.
    3. Kostova, I. Ruthenium complexes as anticancer agents. Curr. Med. Chem. 2006, 13, 1085–1107.
    4. Bandarra, D.; Lopes, M.; Lopes, T.; Almeida, J.; Saraiva, M.S.; Dias, M.V.; Nunes, C.D.; Félix, V.; Brandão, P.; Vaz, P.D.; et al. Mo(II) complexes: A new family of cytotoxic agents? J. Inorg. Biochem. 2010, 104, 1171–1177.
    5. Ikotun, O.F.; Higbee, E.M.; Ouelette, W.; Doyle, R.P. Pyrophosphate-bridged complexes with picomolar toxicity. J. Inorg. Biochem. 2009, 103, 1254–1264.
    6. Liang, X.; Zou, X.; Tan, L.; Zhy, W. Study on nucleic acid (CT-DNA and yeast tRNA) binding behaviors and cytotoxic properties of a heterodinuclear Ru(II)-Co(III) polypyridyl complex. J. Inorg. Biochem. 2010, 104, 1259–1266.
    7. Dasari, S.; Singh, S.; Sivakumar, S.; Patra, A.K. Dual-Sensitized Luiminiscent Europium(III) and Terbioum(III) Complexes as Biomaging and Light-Responsive Therapeutic Agents. Chem. Eur. J. 2016, 22, 1–11.
    8. Viganor, L.; Howe, O.; McCarron, P.; McCann, M.; Devereux, M. The Antibacterial Activity of Metal Complexes Containing 1,10-phenanthroline: Potential as Alternative Therapeutics in the Era of Antibiotic Resistance. Curr. Top. Med. Chem. 2017, 17, 1280–1302.
    9. Páez, P.L.; Bazán, C.M.; Bongiovanni, M.E.; Toneatto, J.; Albesa, I.; Becerra, M.C.; Argüello, G.A. Oxidative Stress and Antimicrobial Activity of Chromium(III) and Ruthenium(II) Complexes on Staphylococcus aureus and Escherichia coli. Biomed. Res. Int. 2013, 2013, 906912.
    10. Viganor, L.; Galdino, A.C.M.; Nunes, A.P.F.; Santos, K.R.N.; Branquinha, M.H.; Devereux, M.; Kellett, A.; McCann, M.; Santos, A.L.S. Anti-Pseudomonas aeruginosa activity of 1,10-phenanthroline-based drugs against both planktonic- and biofilm-growing cells. J. Antimicrob. Chemother. 2016, 71, 128–134.
    11. Mukherjee, A.; Lavery, R.; Bagchi, B.; Hynes, J.T. On the Molecular Mechanism of Drug Intercalation Study of the Intercalation Pathway, Free Energy, and DNA Structural Changes. J. Am. Chem. Soc. 2008, 130, 9747–9755.
    12. Sasikala, W.D.; Mukherjee, A. Intercalation and deintercalation pathway of proflavine through the minor and major grooves of DNA: Roles of water and entropy. Phys. Chem. Chem. Phys. 2013, 15, 6446–6455.
    13. Mukherjee, A.; Sasikala, W.D. Drug-DNA intercalation: From discovery to the molecular mechanism. Adv. Protein Chem. Struct. Biol. 2013, 92, 1–62.
    14. Franco, D.; Vargiu, A.V.; Magistrato, A. Ru[(bpy)2(dppz)]2+ and Rh[(bpy)2(chrysi)]3+ Targeting Double Strand DNA: The Shape of the Intercalating Ligand Tunes the Free Energy Landscape of Deintercalation. Inorg. Chem. 2014, 53, 7999–8008.
    15. Vargiu, A.V.; Magistrato, A. Atomistic-Level Portrayal of Drug-DNA Interplay: A History of Courtships and Meetings Revealed by Molecular Simulations. ChemMedChem 2014, 9, 1966–1981.
    16. Snyder, R.D. Assessment of atypical DNA intercalating agents in biological and in silico systems. Mutat. Res.-Fund. Mol. Mech. Mutagen. 2007, 623, 72–82.
    17. Fletcher, T.M.; Sun, D.; Salazar, M.; Hurley, L.H. Effect of dna secondary structure on human telomerase activity. Biochemistry 1998, 37, 5536–5541.
    18. Mergny, J.L.; Riou, J.F.; Mailliet, P.; Teulade-Fichou, M.P.; Gilson, E. Natural and pharmacological regulation of telomerase. Nucleic Acids Res. 2002, 30, 839–865.
    19. Balasubramanian, S.; Hurley, H.L.; Neidle, S. Targeting G-quadruplexes in gene promoters: A novel anticancer strategy? Nat. Rev. Drug Discov. 2011, 10, 261–275.
    20. Feng, Y.; Yang, D.; Chen, H.; Cheng, W.; Wang, L.; Sun, H.; Tang, Y. Stabilization of G-quadruplex DNA and inhibition of Bcl-2 expression by a pyrodostatin analog. Bioorg. Med. Chem. Lett. 2016, 26, 1660–1663.
    21. Saranathan, N.; Vivekanandan, P. G-Quadruplexes: More Than Just a Kink in Microbial Genomes. Trends Microbiol. 2019, 27, 148–163.
    22. Cao, Q.; Li, Y.; Freisinger, E.; Qin, P.Z.; Sigel, R.K.O.; Mao, Z.-W. G-quadruplex DNA targeted metal complexes acting as potential anticancer drugs. Inorg. Chem. 2017, 4, 10–32.
    23. Yang, D.; Okamoto, K. Structural insights into G-quadruplexes: Towards new anticancer drugs. Fut. Med. Chem. 2010, 2, 619–646.
    24. Hasenknopt, B. Polyoxometalates: Introduction to a class of inorganic compounds and their biomedical applications. Front. Biosci. 2005, 10, 275–287.
    25. Wang, X.; Liu, J.; Li, J.; Yang, Y.; Liu, J.; Li, B.; Pope, M.T. Synthesis and antitumor activity of cyclopentadienyltitanium substituted polyoxotungstate [CoW11O39(CpTi)]7− (Cp = η5-C5H5). J. Inorg. Biochem. 2003, 94, 279–284.
    26. Li, J.; Qi, Y.; Li, H.; Wang, X.; Wu, L.; Duan, E.; Wang, J. Heteropolymolybdate-amino acid complexes: Synthesis, characterization and biological activity. J. Coord. Chem. 2004, 57, 1309–1319.
    27. Bijelic, A.; Aureliano, M.; Rompel, A. Polyoxometalates as Potential Next-Generation Metallodrugs in the Combat Against Cancer. Angew. Chem. Int. Ed. 2018, 58, 2980–2999.
    28. Bijelic, A.; Aureliano, M.; Rompel, A. The antibacterial activity of polyoxometalates: Structures, antibiotic effects and future perspectives. Chem. Commun. 2018, 54, 1153–1169.
    29. Ma, T.; Yang, P.; Damman, I.; Lin, Z.; Mougharbel, S.; Li, M.-X.; Adǎscǎlitei, F.; Mitea, R.; Silvestru, C.; Thorstenson, C.; et al. Tetra-(p-tolyl)antimony(III)-Containing Heteropolytungstates, [{(p-tolyl)SbIII}4(A-α-XW9O34)2]n− (X = P, As, or Ge): Synthesis, Structure, and Study of Antibacterial and Antitumor Activity. Inorg. Chem. 2020, 59, 2978–2987.
    30. Cartuyvels, E.; Absillis, G.; Parac-Vogt, T.N. Questioning the paradigm of metal complex promoted phosphodiester hydrolysis: [Mo7O24]6− polyoxometalate cluster as an unlikely catalyst for the hydrolysis of a DNA model substrate. Chem. Commun. 2008, 2008, 85–87.
    31. Van Lokeren, L.; Cartuyvels, E.; Absillis, G.; Willem, R.; Parac-Vogt, T.N. Phosphoesterase activity of polyoxomolybdates: Diffusion ordered NMR spectroscopy as a tool for obtaining insights into the reactivity of polyoxometalate clusters. Chem. Commun. 2008, 2008, 2774–2776.
    32. Absillis, G.; Cartuyvels, E.; Van Deun, R.; Parac-Vogt, T.N. Hydrolytic Cleavage of an RNA-Model Phosphodiester Catalyzed by a Highly Negatively Charged Polyoxomolybdate [Mo7O24]6− Cluster. J. Am. Chem. Soc. 2008, 130, 17400–17408.
    33. Absillis, G.; Van Deun, R.; Parac-Vogt, T.N. Polyoxomolybdate Promoted Hydrolysis of a DNA-Model Phosphoester Studied by NMR and EXAFS Spectroscopy. Inorg. Chem. 2011, 50, 11552–11560.
    34. Tomé, C.M.; Oliveira, M.C.; Pillinger, M.; Gonçalves, I.S.; Abrantes, M. Use of MoO2Cl2(DMF)2 as a precursor for molybdate promoted hydrolysis of phosphoester bonds. Dalton Trans. 2013, 42, 3901–3907.
    35. Gomes, A.C.; Pillinger, M.; Nunes, P.; Gonçalves, I.S.; Abrantes, M. Promotion of phosphoester hydrolysis by MoO2Cl2L (L = bipyridine derivatives, H2O, no ligand), MoO2(CH3)2L (L = bipyridine derivatives) and related inorganic-organic hybrids in aqueous media. J. Organomet. Chem. 2014, 760, 42–47.
    36. Biancardi, A.; Biver, T.; Marini, A.; Mennucci, B.; Secco, F. Thiazole orange (TO) as a light-switch probe: A combined quantum-mechanical and spectroscopic study. Phys. Chem. Chem. Phys. 2011, 13, 12595–12602.
    37. Zhao, Y.; Truhlar, D.G. Density functionals with broad applicability in chemistry. Acc. Chem. Res. 2008, 41, 157–167.
    38. Zaho, Y.; Truhlar, D.G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: Two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215–241.
    39. Hehre, W.J.; Radom, L.; Schleyer, P.v.R.; People, J.A. Ab Initio Molecular Quantum Theory; Wiley: New York, NY, USA, 1986.
    40. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Sclmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09; Gaussian, Inc.: Wallingford, CT, USA, 2009.
    41. Korth, M.; Pitoňák, M.; Řezáč, J.; Hobza, P. A Transferable H-Bonding Correction for Semiempirical Quantum-Chemical Methods. J. Chem. Theory Comput. 2010, 6, 344–352.
    42. MOPAC; Stewart, J.J.P. Stewart Computational Chemistry, Colorado Springs, CO. 2016. Available online: http://OpenMOPAC.net (accessed on 24 June 2021).
    43. Lee, K.; Murray, É.D.; Kong, L.; Lundqvist, B.; Langreth, D.C. Higher-accuracy van der Waals density functional. Phys. Rev. B 2010, 82, 081101.
    44. Junquera, J.; Paz, O.; Sánchez-Portal, D.; Artacho, E. Numerical atomic orbitals for linear-scaling calculations. Phys. Rev. B 2001, 64, 235111.
    45. Troullier, N.; Martins, J.L. Efficient pseudopotentials for plane-wave calculations. II. Operators for fast iterative diagonalization. Phys. Rev. B 1991, 43, 8861–8869.
    46. Kleinman, L.; Bylander, D.M. Efficacious Form for Model Pseudopotentials. Phys. Rev. Lett. 1982, 48, 1425–1428.
    47. Soler, J.M.; Artacho, E.; Gale, J.D.; García, A.; Junquera, J.; Ordejón, P.; Sánchez-Portal, D. The SIESTA method for ab initio order-N materials simulation. J. Phys. Condens. Matter. 2002, 14, 2745–2779.
    48. Kitaura, K.; Morokuma, K. A new energy decomposition scheme for molecular interactions within the Hartree-Fock approximation. Int. J. Quantum Chem. 1975, 10, 325–340.
    49. von Hopffgarten, M.; Frenking, G. Energy decomposition analysis. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 43–62.
    50. Becke, A.D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648–5652.
    51. Miehlich, B.; Savin, A.; Stoll, H.; Preuss, H. Results Obtained with the Correlation Energy Density Functionals of Becke and Lee, Yang and Parr. Chem. Phys. Lett. 1989, 157, 200–206.
    52. Lee, C.; Yang, W.; Parr, G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B Condens. Matter. 1988, 37, 785–789.
    53. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H.A. Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104.
    54. Řeha, D.; Kabelač, M.; Ryjaček, F.; Šponer, J.; Šponer, J.E.; Elstner, M.; Suhai, S.; Hobza, P. Intercalators. 1. Nature of Stacking Interactions between Intercalators (Ethidium, Daunomycin, Ellipticine, and 4’,6-Diaminide-2-phenylindole) and DNA Base Pairs. Ab Inition Quantum Chemical, Density Functional Theory, and Empirical Potential Study. J. Am. Chem. Soc. 2002, 124, 3366–3376.
    55. Hobza, P.; Šponer, J. Structure, Energetics, and Dynamics of the Nucleic Acid Base Pair: Nonempirical Ab Initio Calculations. Chem. Rev. 1999, 99, 3247–3276.
    56. te Velde, G.; Bickelhaupt, F.M.; Baerends, E.J.; Fonseca Guerra, C.; van Gisbergen, S.J.A.; Snijders, J.G.; Ziegler, T. Chemistry with ADF. J. Comput. Chem. 2001, 22, 931–967.
    57. Fonseca Guerra, C.; Snijders, J.G.; Te Velde, G.; Baerends, E.J. Towards an Order-N DFT Method. Theor. Chem. Acc. 1998, 99, 391–403.
    58. Software Chemistry Materials. Powerful Computational Chemistry. ADF2013 SCM. Available online: http://www.scm.com (accessed on 24 June 2021).
    59. Bader, R.F.W. Atoms in Molecules: A Quantum Theory; Claredon: Oxford, UK, 1990.
    60. Johnson, E.R.; Keinan, S.; Mori-Sánchez, P.; Contreras-Garcia, J.; Cohen, A.J.; Yang, W. Revealing noncovalent interactions. J. Am. Chem. Soc. 2010, 132, 6498–6506.
    61. Friedrich Biegler-König, J.S. Update of the AIM2000-program for atoms in molecules. J. Comput. Chem. 2002, 23, 1489–1494.
    62. Keith, T.A. AIMALL (Version 17.11.14); TK Gristmill Software: Overland Park, KS, USA, 2017.
    63. Mulliken, R.S. Electronic Population Analysis on LCAO-MO Molecular Wave Functions. II. Overlap Populations, Bond Orders, and Covalent Energies. J. Chem. Phys. 1955, 23, 1841.
    64. Cioslowski, J. A new population analysis based on atomic polar tensors. J. Am. Chem. Soc. 1989, 111, 8333–8336.
    65. Hirshfeld, F.L. Bonded-atom fragments for describing molecular charge densities. Theor. Chim. Acta 1977, 44, 129–138.
    66. Reed, A.E.; Weistock, R.B.; Weinhold, F. Natural population analysis. J. Chem. Phys. 1985, 73, 735–746.
    67. Reed, A.E.; Curtiss, L.A.; Weinhold, F. Intermolecular interactions from a natural bons orbital, donor-acceptor viewpoint. Chem. Rev. 1988, 88, 899–926.
    68. Frish, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision A.03; Gaussian Inc.: Wallingford, CT, USA, 2016.
    69. Vosko, S.H.; Wilk, L.; Nusair, M. Accurate spin-dependent electron liquid correlation energies for local spin density calculations: A critical analysis. Can. J. Phys. 1980, 58, 1200–1211.
    70. Becke, A.D. Density-functional exchange energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098.
    71. Perdew, J.P. Density-functional approximation for the correlation energy of the inhomogeneous electron gas. Phys. Rev. B 1986, 33, 8822.
    72. Van Lenthe, E.; Baerends, E.J.; Snijders, J.G. Relativistic total energy using regular approximations. J. Chem Phys. 1994, 101, 9783–9792.
    73. Van Lenthe, E.; Baerends, E.J.; Snijders, J.G. Relativistic regular two-component Hamiltonians. J. Chem. Phys. 1993, 99, 4597–4610.
    74. Klamt, A.; Schuuman, G. COSMO: A new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradiens. J. Chem. Soc. Perkin Trans. 1993, 20, 799–805.
    75. Hay, P.J.; Wadt, W.R. Ab initio effective core potentials for molecular calculations. Potentials for the transition metal atoms Sc to Hg. J. Chem. Phys. 1985, 82, 270–283.
    76. Hay, P.J.; Wadt, W.R. Ab initio effective core potentials for molecular calculations. Potentials for main groups elements Na to Bi. J. Chem. Phys. 1985, 82, 284–298.
    77. Wadt, W.R.; Hay, P.J. Ab initio effective core potentials for molecular calculations. Potentials for K to Au including the outermost core orbitals. J. Chem. Phys. 1985, 82, 299–310.
    78. Check, C.E.; Faust, T.O.; Bailey, J.M.; Wright, B.J.; Gilbert, T.M.; Sunderlin, L.S. Addition of Polarization and Diffuse Functions to the LAN2L2DZ Basis Set for P-Block Elements. J. Phys. Chem. A 2001, 105, 8111–8116.
    79. Ehlers, A.W.; Böhme, M.; Dapprich, S.; Gobbi, A.; Höllwarth, A.; Jonas, V.; Köhler, K.F.; Stegmann, R.; Veldkamp, A.; Frenking, G.A. Set of f-polarization functions for pseudopotential basis sets of the transition metals Sc-Cu, Y-Ag and La-Au. Chem. Phys. Lett. 1993, 208, 111–114.
    80. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465.
    81. Hratchian, H.P.; Schlegel, H.B. Accurate reaction paths using a Hessian based predictor-corrector integrator. J. Chem. Phys. 2004, 120, 9918–9924.
    82. Hratchian, H.P.; Schlegel, H.B. Using Hessian Updating to Increase the Efficiency of a Hessian Based Predictor-Corrector Reaction Path Following Method. J. Chem. Theory Comput. 2005, 1, 61–69.
    83. Hehre, W.J.; Ditchfield, K.; People, J.A. Self-Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian-Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules. J. Chem. Phys. 1972, 56, 2257.
    84. Cances, E.; Menucci, B.; Tomasi, J. A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics. J. Chem. Phys. 1997, 107, 3032–3041.
    85. Cances, E.; Menucci, B. New applications of integral equations methods for solvation continuum models: Ionic solutions and liquid crystals. J. Math. Chem. 1998, 23, 309–326.
    More