Effect of Explicit Hydration on the Cisplatin Reaction Mechanism with Adenine and Guanine

1. Introduction

Cisplatin, cis-Pt(NH₃)₂Cl₂ (1), is one of the quintessential drugs in cancer treatment [1]. Since its discovery in 1965 [2], it has been be widely used in the fight against various types of cancer such as testicular [3], ovarian [4], lung [5], prostate [6], among many others [7]. For almost six decades, the elucidation of its mechanism of action has been a topic of interest from experimental [1] and theoretical [8] perspectives. Cisplatin is a $d^8$ platinum (II) metal complex with a square planar geometry. This complex, however, does not react directly with DNA [9]. It must first undergo hydration reactions, in which it loses one or both chloride anions (Figure 1). These transformations occur via a characteristic substitution reaction mechanism of square planar complexes [10], wherein the transition state has a pentacoordinate trigonal bipyramid geometry [11]. It is still not clear which of the activated complexes, cis-diamminechloroaquaPt(II) (2) or cis-diamminediaquaPt(II) (3) binds directly to the DNA bases[12]. This addition is the rate-limiting step of the whole process. Previous in vitro [13] and computational [14] studies indicated the propensity of cisplatin to bind to the N7 of adenine and guanine nitrogenous bases (Figure 2).

Despite the large body of research about cisplatin, there are still unsolved issues concerning its reaction mechanism with DNA. For example, the role of the solvation water molecules in this mechanism is still unclear. Given the biological relevance of the action pathway of drugs, the availability of a hydration model capable of reliably reproducing experimental values is often pivotal. Due to this circumstance, there are studies that compare distinct solvation methods for biologically relevant reactions. For example, Lakbaibi et. al. evaluated the role of explicit solvation in the Darzens reaction, one of the most useful synthetic procedures to form epoxides, as well as natural and pharmacological products [15]. These authors found that explicit solvation significantly reduces the values of Gibbs free activation energies values. This result was attributed to the capacity of the explicit solvation model to account for crucial hydrogen-bonding interactions, which are not properly described by implicit models. In another study, Wang and Cao [16] investigated the hydration of CO₂ with both explicit and implicit methods of microhydration. Their results show that explicit solvation is critical for a reasonable description of the proton relay mechanism which occurs in this system. Alberto et. al [17]. addressed the hydrolysis mechanism of a platinum derivative anticancer drug [18], Nedaplatin, with an implicit solvation model and the influence of an extra explicit water molecule. These authors concluded that the explicit water molecule incorporates the effects of hydrogen bonding and provides a more accurate description of the stabilisation of the leaving group. More specifically for the investigation reported herein, there are still unsolved issues regarding the role of microsolvation in the mechanism of cisplatin in spite of the considerable interest about this pharmaceutical. A related study on this matter [19] highlighted the significant influence of explicit water molecules for the underlying interactions governing the addition reaction mechanisms of cisplatin with DNA bases. In this investigation, Cossio et. al. concluded that the incorporation of explicit water molecules reduces the Gibbs free activation energy in every step of the action mechanism of cisplatin. The utilisation of explicit microhydration led to computational results which agree closely with experimental observations [19]. Unfortunately, the consideration of electronic energies without explicit H₂O solvation molecules in this same study yields questionable negative activation energies for the second hydration of cisplatin [19]. Besides, the examination of interatomic distances of some reactants and the corresponding transition states indicates that the water molecules surrounding the reactive system interact more strongly with the reactants than with the transition state. Hence, one would expect that the activation energy of such steps of the action mechanism of cisplatin increases as a result of the consideration of explicit microhydration, as opposed to the conclusions reported by Cossio et al. [19]. Therefore, we decided to revisit the effect of explicit microhydration on the biological mechanism of cisplatin. For this purpose, we also exploited state-of-the-art methods of wave function analyses to get further insights about the changes in the chemical bonding scenario in these reactive systems due to the incorporation of explicit H₂O solvation molecules. It is well known, indeed, that non-covalent interactions play a fundamental role in a drug mechanism of action [20]. Since the hydration of cisplatin is the key step for the activation of this drug, it is vital to understand the role of hydrogen bonding in this reaction. Indeed, we have previously investigated in detail the role of H-bonds in the catalysis of chemical reactions [21,22]. More specifically for the matter at hand, Robertazzi et al. [23], considered thermodynamical properties and effects of explicit solvation in the first hydration shell of cisplatin using electronic structure calculations and the Quantum Theory of Atoms in Molecules (QTAIM). As opposed to the conclusions of Cossio et. al. [19], these authors found that explicit solvation substantially distorts the cisplatin geometry and leads to higher activation energies while forming a solvation sphere with ten explicit water molecules.

In view of this background, we aim to elucidate the role of non-covalent interactions in the reactivity of cisplatin across its biological mechanism with the aid of wave function analyses. We also focused on the chemical bonds, which are formed and destroyed throughout the whole reaction mechanism of cisplatin, based on electronic structure theory and Quantum Chemical Topology (QCT) methods, namely QTAIM and the Interacting Quantum Atoms (IQA) energy partition. In short, we address in this study: (i) the steps of the biological mechanism of cisplatin and their corresponding activation energies and (ii) the influence of explicit water molecules on the chemical bonding of the reactants and transition states, for a better understanding of their effect on the activation energies along the biological mechanism of cisplatin. Electronic structure and QCT calculations indicate that the interaction of the explicit solvation water molecules stabilises to a greater extent (i) reactants or (ii) transition states, leading to an increase or decrease of the corresponding activation energy, respectively. Finally, we observe the changes in the chemical bonding scenario of leaving and entering groups throughout the biological mechanism of cisplatin via delocalisation indices as definided by QTAIM. Overall, this paper illustrates how the combination of suitable electronic structure calculations and QCT analyses provides (i) a well-suited assessment of the effects of explicit solvation water molecules on activation energies and (ii) valuable insights into such effects in terms of changes in the chemical bonding scenario within the actin pathway of cisplatin.

Figure 2. Numbering of atoms in guanine and adenine. We emphasise that cisplatin binds N7 when it interacts with either of these nitrogenous bases.

Figure 2. Numbering of atoms in guanine and adenine. We emphasise that cisplatin binds N7 when it interacts with either of these nitrogenous bases.

2. Theoretical Framework

In this section we briefly discuss the QCT methods used in this work, namely the QTAIM and IQA methods of wave function analyses. The QTAIM approach leads to a partition of the 3D-space into disjoints regions ${\Omega }_{\mathrm{A}}$, ${\Omega }_{\mathrm{B}}$, ...identified with the atoms of chemistry [24]. The QTAIM defines expectation values of quantum-mechanical observables within these atomic basins. These observables include energy and multipole moments. For example, the QTAIM atomic charge is defined as,

$$Q^{{\Omega }_{\mathrm{A}}}=Z^{{\Omega }_{\mathrm{A}}}-{\int }_{{\Omega }_{\mathrm{A}}}\rho \left(\mathit{r}\right)\mathrm{d}\mathit{r},$$

(1)

wherein $Z^{{\Omega }_{\mathrm{A}}}$ is the atomic number of the nucleus inside ${\Omega }_{\mathrm{A}}$, $\rho \left(\mathit{r}\right)$ is the electronic charge distribution and ${\int }_{{\Omega }_{\mathrm{A}}}\rho \left(\mathit{r}\right)\mathrm{d}\mathit{r}$ is the expectation value of the number of electrons within ${\Omega }_{\mathrm{A}}$. Topological analyses of the electronic density can be performed in terms of critical points of $\rho \left(\mathbf{r}\right)$, i.e., points where ∇$\rho$ vanishes. The QTAIM allows to compute the number of electrons shared between two atoms, denoted as Delocalization Index (DI),

$$\mathrm{DI}\left({\Omega }_{\mathrm{A}},{\Omega }_{\mathrm{B}}\right)=-2cov\left(N\left({\Omega }_{\mathrm{A}}\right),N\left({\Omega }_{\mathrm{B}}\right)\right),$$

(2)

in which $cov\left(N\left({\Omega }_{\mathrm{A}}\right),N\left({\Omega }_{\mathrm{B}}\right)\right)$ is the covariance between the expectation value of the number of electrons in ${\Omega }_{\mathrm{A}}$ and ${\Omega }_{\mathrm{B}}$. The value of $\mathrm{DI}\left({\Omega }_{\mathrm{A}},{\Omega }_{\mathrm{B}}\right)$ yields a measure of the exchange of electrons between ${\Omega }_{\mathrm{A}}$ and ${\Omega }_{\mathrm{B}}$, the distinctive feature of covalency.

On the other hand, the IQA energy partition can be used to divide the electronic energy of a system into net, $E_{\mathrm{net}}^{{\Omega }_{\mathrm{A}}}$, and interatomic, $E_{\mathrm{int}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}$, contributions,

$$E=\sum _{\mathrm{A}}{E}_{\mathrm{net}}^{{\Omega }_{\mathrm{A}}}+\sum _{\mathrm{A}}\sum _{\mathrm{B}>\mathrm{A}}{E}_{\mathrm{int}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}.$$

(3)

The quantity $E_{\mathrm{net}}^{{\Omega }_{\mathrm{A}}}$ is the net energy of atom A, while the value of $E_{\mathrm{int}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}$ is the interaction energy between atoms ${\Omega }_{\mathrm{A}}$ and ${\Omega }_{\mathrm{B}}$. In turn, we can split the IQA interaction energy into a classical, $E_{\mathrm{cl}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}$, and an exchange-correlation, $E_{\mathrm{xc}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}$ component, i.e.,

$$E_{\mathrm{int}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}=E_{\mathrm{cl}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}+E_{\mathrm{xc}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}},$$

(4)

wherein $E_{\mathrm{cl}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}$ and $E_{\mathrm{xc}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}$ are respectively related to ionic and covalent contributions to the chemical bond between basins ${\Omega }_{\mathrm{A}}$ and ${\Omega }_{\mathrm{B}}$. These quantities can be approximated using QTAIM charges and DIs [25] as follows,

$$\begin{array}{cc}\hfill E_{\mathrm{cl}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}& \approx {\displaystyle \frac{Q^{{\Omega }_{\mathrm{A}}}{Q}^{{\Omega }_{\mathrm{B}}}}{R}},\hfill \end{array}$$

(5)

$$\begin{array}{cc}\hfill E_{\mathrm{xc}}^{{\Omega }_{\mathrm{A}}{\Omega }_{\mathrm{B}}}& \approx -{\displaystyle \frac{\mathrm{DI}\left({\Omega }_{\mathrm{A}},{\Omega }_{\mathrm{B}}\right)}{2R}},\hfill \end{array}$$

(6)

in which R is the distance between the nuclei within ${\Omega }_{\mathrm{A}}$ and ${\Omega }_{\mathrm{B}}$.

3. Computational Details

We performed geometry optimisations of all Pt complexes involved in the cisplatin reaction mechanism using both the M06-2X [26] and PBE0 [27] exchange-correlation functionals as implemented in the Orca 5.0.3 program package [28,29,30]. For these electronic structure calculations, we used the def2-SVP basis set for all the main-group elements and def2-TZVP basis set for platinum (hereafter referred to as Basis Set 1), along with a Stuttgart-Dresden relativistic effective core potential for the Pt centre [31]. To identify the most suitable exchange-correlation functional to carry out the QCT computations, we benchmarked the first hydration step of cisplatin using the more extensive def2-QZVP basis set (Basis Set 2) developed by Ahlrichs [32]. The PBE0 functional was used together with the Grimme’s D3 dispersion correction [33,34], as this combination accurately accounts for non-covalent interactions such as hydrogen bonding between different molecular strands which are critical factors in the second addition of nitrogenous bases to cisplatin. Moreover, PBE0+D3 reliably captures the structure and energetics of heavy metal complexes [35], while M06-2X agrees very well with experimental results for the hydration steps of cisplatin [36,37]. Besides, the last mentioned exchange-correlation functional is particularly accurate to predict energy barriers across reaction mechanisms [38]. Notably, both functionals yielded similar trends throughout the reaction pathway involving (i) cisplatin and (ii) adenine or guanine. We emphasise that all calculations accounted for solvent effects via the conductor-like polarizable continuum model (C-PCM) [39]. We proceeded in this way, so that we could better assess the effect of the sorrunding water molecules on the reactive system.

In order to generate suitable electron densities for QTAIM analyses, we performed single-point calculations for every minimum and transition state using the Zero-Order Regular Approximation (ZORA) [40]. These electronic structure computations were carried out with relativistically contracted orbital and auxiliary basis sets [41]. Finally, QTAIM analyses [24] were performed using the Aimall program [42].

4. Results and Discussion

We discuss first the hydration of cisplatin as schematised in Figure 1. It is well known that the first hydration of cisplatin is the rate-limiting step for its activation in the cell [7]. Therefore, we set out to investigate under different approximations the effect of microhydration on the activation energy for this first hydration step. We considered two explicit water molecules along the entire reaction mechanism to have comparable results with previously reported activation energies. As Cossio et al. [19] wrote: “going from zero to two surrounding water molecules, the computed activation and reaction energies can vary substantially, particularly the Gibbs activation and reaction energies. However, our results also show that for a reliable representation of the system, only two molecules of water need to be considered ” [19]. The corresponding results are shown in Table 1, where the consideration of microhydration in the system increases the electronic and the Gibbs free activation energies in contrast with the results of reference [19]. We also observed that there are no large differences between $\Delta E^{‡}$ and $\Delta G^{‡}$ (>1.7 kcal/mol).The electronic activation energy computed with the M06-2X approximation along with Basis Set 1 with implicit solvation effects (20.57 kcal/mol) agrees with that experimentally determined by Repta and Long [37] (19–20 kcal/mol). However, by using the same approximate functional and basis set, but considering explicit microsolvation (22.04 kcal/mol), our results agree with the experimental findings of Brancroft et al. [36] (22.25 kcal/mol). Regardless of the consideration of microsolvation, the values of $\Delta G_{\mathrm{rxn}}$, i.e., the change in Gibbs free energy for the overall reaction, indicate that the reaction is endergonic. Nevertheless, microsolvation results in a higher value in $\Delta G_{\mathrm{rxn}}$ as represented in Table 1. In summary, implicit and explicit solvation represent experimental results adequately. Therefore, we will consider in the following section both sets of results with and without the effect of explicit water molecules to understand the role of microsolvation in the reaction mechanism of cisplatin with the nitrogenous bases adenine and guanine.

4.1. Monofunctional and Bifunctional Addition of DNA Bases to Cisplatin

Another significant step in the action pathway of cisplatin is the addition of nitrogenous bases. Table 2 presents the results concerning activation enthalpies for the first and second functionalisations of cisplatin with guanine as schematised in Figure 3. Relevant effects occur with the inclusion of explicit water molecules for each of the additions of the nitrogenous bases. We found for the monofunctionalisation of cisplatin that explicit microsolvation decreases the activation energy, while for the second functionalisation, it presents the opposite effect. The M06-2X functional closely aligns with the experimental activation enthalpy for the first functionalisation of cisplatin (18 ± 1 kcal/mol) [36] for both implicit and explicit solvation models. The experimental activation enthalpy for the subsequent addition of guanine is 21 ± 2 kcal/mol, a value which is accurately reproduced by the PBE0 functional with the Grimme’s correction for dispersion and explicit solvation water molecules. Still, the M06-2X functional with explicit solvation remains competitive in modelling the bifunctionalisation reaction of guanine failing just by 1 kcal/mol in the superior limit of the experimental data and it yields significantly better results regarding the first functionalisation of cisplatin. Thus, the remainder of the discussion will focus exclusively on the results obtained using the M06-2X functional.

4.2. Energetics for the Reaction of Cisplatin with Bases Adenine and Guanine

We will examine in this section the biological addition of adenine and guanine to cisplatin (Figure 4) with implicit and explicit solvation effects. More specifically, we will analyse in detail each of the reaction steps for the addition of these nitrogenous bases to cisplatin along with the corresponding changes in their electronic activation energies due to explicit microhydration. As shown in Figure 4, the mechanisms with guanine and adenine share certain steps, viz., the hydration of cisplatin (TS1, TS2). The mechanism involving the addition of adenine consists of transition states TS3A–TS7A. Ditto for the addition of guanine and the transition states TS3G–TS7G. The profiles corresponding to the Gibbs free activation energies exhibit the same trends as the activation electronic energies (see Figures S1–S24 in the Supporting Information). Given the lack of a direct relationship between (i) the respective quality and nucleophilicity of the leaving and entering groups and (ii) changes in the activation energies, we ruled out that the entering or leaving group determine the variations in activation energies. As mentioned above, a preliminary examination of interatomic distances led us to realise that in cases where the activation energy increases due to the contemplation of explicit water molecules (see e.g TS4A and TS6G in Table 3), the apparent role of the surrounding water molecules is a greater stabilisation in the reactants than in the corresponding transition states. The consideration of the differences in electronic activation energies,

$$\Delta \Delta E^{‡}=\Delta E^{‡}\left(\mathrm{explicit}\right)-\Delta E^{‡}\left(\mathrm{implicit}\right),$$

(7)

in Table 3 shows that the explicit consideration of solvation water molecules increases overall the activation energy of the reaction steps considered in Figure 4 with a few exceptions. These results are in direct opposition with the conclusions offered by Cossio et. al. [19]. The lower and upper limits of $\Delta \Delta E^{‡}$ are −4.58 and 15.07 kcal/mol. We considered for further wave function analyses the reactants and transition states of the steps for which $|\Delta \Delta E^{‡}|>$ 3 kcal/mol, i.e, those elementary reactions for which the incorporation of explicit water solvation molecules is most conspicuous.

4.3. QTAIM and IQA Wave Function Analyses

In order to better understand the effect of the surrounding water molecules on the activation energies of the above-mentioned reaction steps, we determined the IQA interaction energy $\left(E_{\mathrm{int}}\right)$ of the atoms with which the encompassing water molecules form a bond path with the reactive system, as defined in equations (4)−(6) of the Theoretical Framework section. Table 4 reports these differences for these IQA interaction energies for the transition state and the corresponding reactants,

$$\Delta E_{\mathrm{int}}^{H2O\cdots \mathrm{cisPt}}=E_{\mathrm{int}}^{H2O\cdots \mathrm{cisPt}}\left(TS\right)-E_{\mathrm{int}}^{H2O\cdots \mathrm{cisPt}}\left(\mathrm{reactants}\right).$$

(8)

When $\Delta E_{\mathrm{int}}^{{\mathrm{H}}_2\mathrm{O}\cdots \mathrm{cisPt}}>0$, then the surrounding water molecules stabilise the reactants more strongly than the transition state and therefore we would expect that the encompassing H₂O monomers increase the activation energy of the reaction under consideration, i.e., $\Delta \Delta E^{‡}>0$. Correspondingly, when $\Delta E_{\mathrm{int}}^{{\mathrm{H}}_2\mathrm{O}\cdots \mathrm{cisPt}}<0$, then the encircling water molecules stabilise the transition state more substantially than the reactants, and hence the solvation H₂O monomers conduce to a reduction of the activation energy, viz., $\Delta \Delta E^{‡}<0$. There is indeed a correspondence between the signs of $\Delta \Delta E^{‡}$ and $\Delta E_{\mathrm{int}}^{{\mathrm{H}}_2\mathrm{O}\cdots \mathrm{cisPt}}$ for all the systems shown in Table 4 apart from TS1, i.e., the transition state with the smallest value of $\Delta \Delta E^{‡}$ for which we performed the QTAIM and IQA analyses.

It is noteworthy that the interaction energies ($\Delta E_{\mathrm{int}}^{{\mathrm{H}}_2\mathrm{O}\cdots \mathrm{cisPt}}$) differ significantly between the TS3A and TS3G steps, despite their expected similarity. In TS3A, the interaction is highly attractive, partly due to a charge transfer of 0.03 a.u. from the water molecules to cisplatin. In contrast, TS3G shows strongly repulsive interactions, as there is no such charge transfer. Additionally, TS3A forms a closed coordination cycle, whereas in TS3G, one of the water molecules primarily acts as a spectator, accumulating a charge of $-0.02\mathrm{a}.\mathrm{u}.$. This leaves the other water molecule to engage in predominantly electrostatically repulsive interactions with the cisplatin molecule.

Finally, to get further insights into the adidtion of nitrogenous bases to cisplatin, we used QTAIM delocalisation indices as chemical bond descriptors. We proceeded in this way to describe the changes in the bond formation from reactants to transition states. There is a displacement of bonded species with the arrival of a nucleophilic ligand in substitution reactions of square planar complexes. The substitution takes place via an associative mechanism in which the nucleophilic ligand first binds to the metal centre, forming a trigonal bipyramidal structure in the transition state. The leaving group is subsequently replaced by the entering nucleophile. More concretely, the entering and leaving species are in equatorial positions throughout the geometric rearrangement of the transition state and the spectator ligands do not affect the direct substitution of the reaction. In the case of the reaction of cisplatin with nitrogenous bases of DNA, the entering group will always be the nitrogenous base and the leaving group will be either chloride, water, or a second nitrogenous base depending on the reaction step being analysed. Table 5 shows the variations in the values of DIs for the leaving and entering groups involved in the formation and breaking of chemical bonds of each transition state with explicit and implicit solvation models for the structures shown in Figure 5. The corresponding changes in delocalisation indices are,

$$\Delta {\mathrm{DI}}^{\mathrm{X}}=\mathrm{DI}{(\mathrm{Pt},\mathrm{X})}^{\mathrm{TS}}-\mathrm{DI}{(\mathrm{Pt},\mathrm{X})}^{\mathrm{reactant}},$$

(9)

wherein X is an entering or leaving group as shown in Figure 5. The greater the absolute value of $\Delta {\mathrm{DI}}^{\mathrm{X}}$, the greater the change in the chemical bonding of the Pt–X interaction. As expected, $\Delta {\mathrm{DI}}^{\mathrm{X}}$ for the entering and leaving groups are positive and negative respectively (Table 5). The addition of $\Delta \mathrm{DI}$, for the leaving and the entering groups,

$$\sum _{\mathrm{X}}{\mathrm{DI}}^{\mathrm{X}}=\Delta {\mathrm{DI}}^{\mathrm{entering}}+\Delta {\mathrm{DI}}^{\mathrm{leaving}},$$

(10)

represents a measure of the effects of the rupture and the formation of chemical bonds from the reactants to the transition state. Table 6 reports the values of ${\sum }_{\mathrm{X}}\Delta {\mathrm{DI}}^{\mathrm{X}}$ for the implicit and the explicit solvation models addressed herein. Because

• the breaking of chemical bonds involves an energy cost, while the formation of these interactions release energy to the surroundings and,
• the energy of the transition states is higher than those of the reactants in every case (all values of $\Delta E^{‡}>0$ in Table 3),

we expect the effects of the rupture of chemical bonds to be more important than those of the formation of these contacts. In other words, we anticipate that the value of ${\sum }_{\mathrm{X}}\Delta {\mathrm{DI}}^{\mathrm{X}}<0$ since $\Delta {\mathrm{DI}}^{\mathrm{X}}$ is positive for the entering group and negative for the leaving group, being the latter influence more relevant as mentioned above. Table 6 shows that this is indeed the case. The difference,

$$\Delta \Delta \mathrm{DI}\left(\mathrm{Pt}\right)=\sum _{\mathrm{X}}\Delta {\mathrm{DI}}^{\mathrm{X}}\left(\mathrm{expl}\right)-\sum _{\mathrm{X}}\Delta {\mathrm{DI}}^{\mathrm{X}}\left(\mathrm{impl}\right),$$

(11)

indicates whether the effects of rupture of chemical bonds are more important for the explicit ($\Delta \Delta \mathrm{DI}\left(\mathrm{Pt}\right)<0$) or the implicit ($\Delta \Delta \mathrm{DI}\left(\mathrm{Pt}\right)>0$) solvation models. Apart from TS1, which is the transition state with the smallest value of $|\Delta \Delta E^{‡}|$ for which we carried out QCT analyses, we note that for those transition states wherein $\Delta \Delta \mathrm{DI}\left(\mathrm{Pt}\right)<0$, it holds that $\Delta \Delta E^{‡}>0$. Indeed, the inclusion of explicit solvation water molecules makes the effect of rupture of chemical bonding from reactants to transition state more stringent and, hence, a higher cost in activation energy is observed. We note the opposite behaviour when $\Delta \Delta \mathrm{DI}\left(\mathrm{Pt}\right)>0$. These observations are consistent with Hammond’s postulate “structures close in energy that transform directly into each other are also similar in structure” [43]. In other words, when the inclusion of explicit solvation water molecules results in a more/less drastic change in the chemical bonding scenario, we expect a rise/diminution in activation energy with respect to implicit solvation models.

5. Conclusions

We examined the reaction mechanism of cisplatin with guanine and adenine, emphasising the impact of microsolvation on activation energies, which are shown to closely match reported experimental values. Our findings highlight that the presence of explicit water molecules can significantly influence activation energies by preferentially stabilising either the reactants or the transition state. State-of-the-art methods of wave function analyses, namely, the QTAIM and IQA approaches, have been instrumental in elucidating the role of microsolvation by providing detailed insights into the interaction energies and changes in the chemical bonding landscape of reactants and transition states. Overall, this research enhances our understanding of how explicit solvation shapes the biological mechanism of cisplatin, demonstrating that microsolvation can modulate activation barriers through a complex interplay of structural and energetic factors, consistent with Hammond’s postulate.