Next Article in Journal
Bladder Cancer Diagnosis and Follow-Up: The Current Status and Possible Role of Extracellular Vesicles
Next Article in Special Issue
Platinum(II) O,S Complexes Inhibit the Aggregation of Amyloid Model Systems
Previous Article in Journal
High-Throughput Virtual Screening, Molecular Dynamics Simulation, and Enzyme Kinetics Identified ZINC84525623 as a Potential Inhibitor of NDM-1
Previous Article in Special Issue
Reaction with Proteins of a Five-Coordinate Platinum(II) Compound
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 4
10.3390/ijms20040820
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reactivity of Gold(I) Monocarbene Complexes with Protein Targets: A Theoretical Study

by
Iogann Tolbatov
,
Cecilia Coletti
,
Alessandro Marrone
and
Nazzareno Re
*
Dipartimento di Farmacia, Università degli Studi “G. D’Annunzio” Chieti-Pescara, Via dei Vestini, I-66100 Chieti, Italy
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(4), 820; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20040820
Submission received: 15 January 2019 / Revised: 6 February 2019 / Accepted: 12 February 2019 / Published: 14 February 2019
(This article belongs to the Special Issue Interaction between Metal Compounds and Proteins)
Browse Figures
Versions Notes

Abstract

:
Neutral N–heterocyclic carbene gold(I) compounds such as IMeAuCl are widely used both in homogeneous catalysis and, more recently, in medicinal chemistry as promising antitumor agents. In order to shed light on their reactivity with protein side chains, we have carried out density functional theory (DFT) calculations on the thermodynamics and kinetics of their reactions with water and various nucleophiles as a model of plausible protein binding sites such as arginine, aspartic acid, asparagine, cysteine, glutamic acid, glutamine, histidine, lysine, methionine, selenocysteine, and the N-terminal group. In agreement with recent experimental data, our results suggest that IMeAuCl easily interacts with all considered biological targets before being hydrated—unless sterically prevented—and allows the establishment of an order of thermodynamic stability and of kinetic reactivity for its binding to protein residues.

Graphical Abstract

1. Introduction

The administration of gold complexes has been assessed for various pharmaceutical aims, particularly in the treatment of cancer and arthritis [1,2,3,4,5,6,7]. A plethora of gold(I) and gold(III) complexes with different ligands display substantial cell-growth-inhibiting properties [4,8]. Nevertheless, taking into account the great structural variation of their ligands, it seems that there is no common mode of action that could describe all of them [4,6,8]. Gold compounds target various cellular processes that ultimately activate apoptosis: direct DNA damage, mitochondrial damage via thioredoxin reductase inhibition, proteasome suppression, the modification of the cell cycle, and the modulation of specific kinases [1,2,3,4,5,6,9,10,11,12].
Gold(I) N–heterocyclic carbenes (NHC) of the general formula [Au(I)(NHC)R], in which R is an anionic or a neutral ligand such as Cl−, PH3, or another NHC moiety, leading to a neutral or to a cationic complex, have been extensively utilized in homogeneous catalysis [13,14,15,16] and, more recently, in medicinal chemistry [4,17,18]. Their unique chemical and biological properties make them successful anti-arthritic [19], antibacterial [20], and anticancer [18,21,22,23,24,25,26,27] agents. At variance with platinum-based drugs, DNA is not the main pharmaceutical target for these gold-based cytotoxic agents that, conversely, mainly act through the modification of selected proteins with consequent loss of function. This implies that selective ‘‘protein metalation’’ is the key feature in the mechanism of action of anticancer gold drugs. In particular, several studies have recently highlighted that the crucial role in the mode of action of N–heterocyclic carbene gold(I) complexes emerges from an efficient selective inhibition of the enzyme thioredoxin reductase (TrxR) [1,7,25,28,29]. Electrospray ionization mass spectrometry (ESI-MS) measurements indicate that gold carbene complexes preferentially attach to free cysteine and selenocysteine side chains [19,30], with the thiol/selenol group being coordinated by [Au(NHC)]+ fragments or Au+ ions upon the release of carbene ligands. Linear Au(I) complexes have also been shown to be effective inhibitors of the Se-free enzyme glutathione reductase (GR), owing to the high thiol reactivity toward gold, as plainly manifested by the crystal structure of GR inhibited by a gold(I) complex displaying Au(I) attached to the active site of Cys thiol with linear S–Au–S coordination [23,27,31]. Detailed 1H NMR studies of the ligand exchange reactions of cationic Au(I)(NHC) complexes with cysteine and selenocysteine reveal a two-step reaction comprising the consecutive substitution of both NHC ligands to yield the [Au(Cys)2] compound [25]. Although as a soft metal ion Au(I) shows a marked preference for soft ligands such as thiols of cysteines and thioethers of methionines [32], multiple X-ray crystallography data have shown that Au(I) ions can coordinate solvent-exposed His [33,34,35], even in the presence of free thiols. In particular, a recent X-ray study [36] has shown that in the absence of free cysteines, methionines, and histidines, these Au(I)(NHC) gold compounds are able to bind the side chains of Lys and Arg residues or the N-terminal Ala amino group.
Even though the interaction of gold carbene compounds with many proteins (for example, TrxR, phosphatases, glutathione reductases, serum albumin, Atox-1, or thaumatin) has now been widely examined [1,7,17,19,25,28,29,30,36,37,38,39,40,41,42,43,44], their detailed mechanisms of reaction with protein side chains are still lacking. In this study, we compute the reactivity of the neutral Au(I)NHC complex (Me2Imy)AuCl, or IMeAuCl, with water, and with the main binding sites in a protein or a polypeptide, i.e., the side chains of arginine, aspartic acid, asparagine, cysteine, glutamic acid, glutamine, histidine, lysine, methionine, and selenocysteine residue, and the N-terminal group. We have used short models of these residues (Scheme 1). Each side chain was represented by the nucleophilic group, while the rest of the chain was represented by the ethyl group. A deeper insight into the binding mode of these Au(I) complexes with protein targets and their binding preference would be very useful to fully comprehend their mechanisms of action in vivo and may be useful to conceive novel and more effective gold(I) anticancer drugs.

2. Results and Discussion

A preliminary investigation on the isolated IMeAuCl complex has shown that at the selected level of theory, the minimized structure matches the crystallographic data [45] accurately with Au–C and Au–Cl bond distances within 0.05 Å, and Cl–Au–C angles within 2°, thus indicating that the B3LYP functional provides a good description of this molecular structure (Figure 1, Table S1).
We first considered a dissociative mechanism with detachment of a chloride anion. The reaction barrier of 26.5 kcal/mol and reaction energy of 13.1 kcal/mol make the dissociation unlikely. DFT calculations were therefore carried out on the thermodynamics and the kinetics of the ligand exchange of one chloride ligand by incoming nucleophile molecules. Simplified models of the residue’s side chains (Scheme 1) were preferred to either free or capped amino acids for several reasons. First, free amino acids are less representative of protein residues due to the presence of terminal amine and carboxylic acid groups, in ionized zwitterionic form, which do not exist in proteins. The capped forms, although chemically resembling protein residues, are connected to the side chain nucleophilic group through hydrocarbon chains of different lengths, thus leading to large variations in the size of the ligand reacting with the metal complex. This latter aspect may bias the calculation of solvation free energies, thus affecting the comparison of amino acid reactivity. All residues were assumed to be in their most stable protonation state at pH=7.2, i.e., positive for arginine, lysine, and N-terminal alanine; neutral for histidine, cysteine, methionine, and asparagine/glutamine; and anionic for selenocysteine and glutamic/aspartic acid. However, for arginine, lysine, and N-terminal alanine, we considered as nucleophiles the neutral form taking into account the deprotonation energy. For cysteine, we considered both the neutral form, expected to be the most stable at neutral pH and the anionic form, which is still present in low concentrations, or is even the most stable form, when cysteine is close to basic residues (such as histidine) which are known to stabilize the anionic form.
The reaction was investigated, taking into account that both reagents and products can form stable non-covalent adducts both prior to and after the reaction takes place. We have thus optimized the geometries of the reactant adducts (RA), transition states (TS), product adducts (PA) intermediates, and their energies with respect to the isolated species that have been evaluated and given below. The activation enthalpies and free energies have been calculated as the difference between TS and the lowest between reactants and reactant adducts, while the reaction enthalpies and free energies have been calculated as difference between reagents and products infinitely apart. The calculated values for the IMeAuCl reactions with all considered molecular models are reported in Table 1, and the results allow one to establish an order of thermodynamic stability and kinetic reactivity for the Au(I) binding to the considered protein residues side chains.
The results were compared to recent experimental X-ray and ESI-MS studies of the binding of this Au(I)NHC neutral complex to peptides and proteins [30,46] showing Au(I) binding mainly to cysteine and selenocysteine, and to the more recent results of the paper by Merlino et al. [36] reporting the first distinct crystal structure of a gold N–heterocyclic carbene compound bound to a protein—thaumatin—showing the Au(I) metal center binding lysine, arginine side chains, as well as the N-terminal group.
We first considered the hydrolysis reaction with substitution of the chlorido ligand by a water molecule to probe the stability of this chloro NHC complex in biological fluids. The transition state (imaginary frequency of −62.92 cm−1) (Figure 2 and Figure S2) for the reaction with water is characterized by the entering water molecule approaching the metal center with Au–O distance of 2.35 Å, whereas the leaving chlorine moves further at the distance of 2.73 Å.
An approximately planar trigonal geometry is observed, with an acute leaving ligand—metal—entering ligand angle of 70.4°. Table 1 and Figure 3 and Figure S3 show that this reaction is slightly endergonic, with a reaction-free energy of 1.3 kcal/mol.
The calculated activation-free energy of 23.0 kcal/mol predicts this reaction to be relatively fast at the physiological temperature, although its barrier is significantly higher than that of all other nucleophiles examined in the present study, indicating that the IMeAuCl complexes are likely to attack the biological targets before being hydrated. Although there are no experimental data for the barrier of the hydrolysis of any [Au(I)(NHC)Cl] complex, nor is there any clear indication of its stability in solution or neutral buffer, the calculated value of the activation-free energy is reasonably consistent with the sparse experimental evidence showing that chlorido 1-butyl-3-methyl-imidazol-2-ylidene gold(I) complex is stable for at least 24 hours at room temperature [47], while chlorido 1,3-methyl-benzimidazol-2-ylidene gold(I) complex undergoes hydrolysis processes within a few hours in buffer solution pH 7.0 incubated at 37 °C [48]. The B3LYP results for this prototypical process were benchmarked against latest-generation density functionals, such as M06-2X, wB97X, and CAM-B3LYP-D3, which have been recently shown to be among the best-performing density functionals in the calculations on main group molecules and metal complexes [49,50,51] and against more accurate ab initio calculations such as MP2 with larger basis sets. The results are reported in Table S2 and show that reaction and activation enthalpies and free energies within 2–3 kcal/mol form the B3LYP results, thus indicating that the B3LYP functional provides a good description of the energetics of this kind of process.
Transition states for methionine, cysteine, and selenocysteine in their neutral form (with imaginary frequencies of −50.43, −44.33, and −18.09 cm−1, respectively) have similar planar trigonal geometries, with gold–sulfur/selenium and gold–chlorine distances being 2.6–2.7 A and the S/Se–Au–Cl angle varying within 82–88°, and show analogous reaction profiles with reaction activation barriers of 18.2, 18.1, and 17.6 kcal/mol, respectively. When cysteine and selenocysteine are considered in their anionic form, lower activation energies have been calculated, respectively, 15.5 and 12.8 kcal/mol, as expected on the basis of the better nucleophilicity provided by the negative charge.
Nucleophiles attacking with their N-based moieties, His, Arg, Lys, and N-terminal amine, form planar trigonal transition states with slightly smaller angles within 77–79°, and Au–N distances of 2.3–2.4 Å, with imaginary frequencies of −59.6 −61.8, −60.8, and −52.9 cm−1 (Figure 2 and Figure S2). Reaction activation barriers for the N-based moieties are also quite low in this case: 20.6 kcal/mol for histidine, 18.3 kcal/mol for arginine, 16.2 kcal/mol for Lys, and 18.8 kcal/mol for N-terminal. When binding to glutamine/asparagine is considered, two different tautomers could in principle be taken into account, corresponding to the amide and imide forms (Scheme 2): for both tautomers, a TS was found with geometries similar to those of the other N-based nucleophiles and activation-free energies of, respectively, 24.8 and 16.8 kcal/mol.
The barrier for the Cl substitution by the imidic form is much lower than that of the amidic form, but it should be born in mind that the imidic form is significantly less stable, by 16.4 kcal/mol, and is stabilized only in few protein systems by suitable adjacent catalytic residues [52] (Figure S1).
Finally, binding to glutamic/aspartic acid was considered, and we found a TS with the usual distorted trigonal geometry and a relatively high activation free energies of 24.7 kcal/mol.
The results for chloride substitution from the other residue models are reported in Table 1 and show that all reactions are slightly to moderately exergonic with the following order of thermodynamic stability:
Sec > activated-Cys > Arg > Lys > His > activated-Gln/Asn > N-terminal > Met > Cys > Glu/Asp > Gln/Asn >
although the kinetic reactivity order is quite different:
Sec > Lys > activated-Gln/Asn > activated-Cys > Cys > Met > Arg > N-terminal > His > Gln/Asn ≈ Glu/Asp >
These results are consistent with the experimental evidence showing that the selenocysteine residue of thioredoxine, and, to a lesser extent, the histidine-activated cysteine residues of several enzymes such as glutathione reductase, react promptly with neutral Au(I)(NHC) complexes and are considered the main physiological targets of this class of compounds. In particular, the low barriers for Cys and especially Sec are in agreement with the results of a recent study on the thioredoxin reductase inhibition by [Au(NHC)Cl] and on the reaction of the same complex with the reduced linear dodecapeptide Ac-SGGDILQSGCUGNH2 investigated by ESI-MS, showing (i) a relevant inhibition of thioredoxin reductase activity soon after mixing (about 30 s) and (ii) the prevalent binding of a [Au(I)(NHC)]+ fragment to the selenocysteine residue and, to a much lesser extent, to cysteine [30].
On the other hand, glutamic/aspartic acid and glutamine/asparagine show high barriers around 25 kcal/mol, and slightly exergonic or even endergonic reaction free energies in agreement with the experimental evidence, indicating that such residues are not good binding sites for [Au(I)(NHC)]+ complexes. However, the interaction of some Au(III) complexes with hen egg white lysozyme has been observed to lead, upon metal reduction and ligand loss, to the binding of an [Au(H2O)]+ gold(I) fragment to a glutamine residue, probably preceded by a amide-imidic acid tautomerization of the glutamine [46].
Binding to histidine shows the highest barrier after glutamic/aspartic acid and glutamine/asparagine, 20.6 kcal/mol, and is therefore slower than the substitution by all other residues, although the barrier is low enough to make the reaction still relatively fast at physiological temperature. However, [IMeAu]+ binding to histidine is one of the most exergonic processes involving the neutral residues, so it can be considered to be plausible when metal binding to protein is under thermodynamic control. This result is consistent with previous X-ray studies showing that the reaction of an analogous gold(I) complex—[Au(PEt3)Cl]—with cyclophilin-3 forms an [AuPEt3]+ adduct with an histidine residue, despite the presence of four Cys thiol groups in the protein [53], and that the interaction of several Au(III) complexes with hen egg white lysozyme preferentially leads, upon metal reduction and ligand loss, to the binding of a AuCl gold(I) fragment to histidine, even though several methionine groups are available in the protein [34].
Relatively low barriers are calculated for the N-based nucleophile sites—lysine, arginine, and N-terminal amine (16–19 kcal/mol)—indicating little discrimination in terms of reactivity within either S/Se or N-based nucleophile sites and suggesting that these N-based sites could also bind an [Au(NHC)]+ moiety whenever cysteine or selenocysteine residues are absent or buried, and that the binding preference to a protein side chain could be determined by its accessibility rather than by its kinetic preference toward chloride substitution. These results are in agreement with a recent crystallographic study showing that the reaction of 1-butyl-3-methyl-imidazol-2-ylidene gold(I) complex with the model protein thaumatin leads to the binding of the [Au(NHC)]+ fragment to lysine and arginine side chains and to the N-terminal amine [36] and suggest that the accessibility of amino acids within the protein rather than reactivity could be the key point to comprehend the amino acid targeting of the [Au(NHC)Cl] complex displayed in the experiment [36].
To better investigate the latter point, solvent-accessible surface (SAS) analyses (Figure 4 caption) were carried out on the X-ray structure of the plant protein thaumatin-1 (pdb id 3qy5) showing that Lys, Arg, Glu, Asp, Gln, Asn, and the only Met residues in thaumatin are markedly more exposed than Cys and N-term, while His residues are absent.
As shown in Figure 4, the gold(I) centers are always bound to highly-exposed amino acids—mainly Arg and Lys—with the exception of N-terminal amine, which undergoes metal binding despite being less exposed. However, the binding of the [IMeAu]+ fragment to the N-terminal amine [36] could be related to the well-known flexibility of the terminal segment of the protein, whose SAS values may vary within a large range in solution despite being restricted to smaller values by the crystal packing. Notably, the preferential Au-binding on Lys106 compared to the sequentially close and similarly exposed Met112 is consistent with the lower values of either reaction or activation-free energies estimated for Lys.

3. Materials and Methods

All calculations were performed with the Jaguar 9.1 [54,55] and Gaussian 09 A.02 [56] quantum chemistry package. Optimizations were carried out in gas phase using the density functional theory (DFT) with the B3LYP hybrid functional [57,58], which is known to give a good description of geometries and reaction profiles for transition metal-containing compounds [59], and particularly antitumor compounds based on metals in different oxidation states [60,61,62,63,64,65,66,67,68] including gold(I) [28,69,70] and gold(III)-based [71,72] anticancer compounds. Other latest-generation functionals explicitly developed for systems with long-range electron correlations were used for benchmarking such as M06-2X [73], wB97X [74], and CAM-B3LYP-D3 [75], which have been recently shown to be among the best-performing density functionals in the calculations on main group molecules and metal complexes [49,50,51], as well as explicitly correlated wave function based methods, such as the Møller–Plesset perturbative approach MP2 [76,77,78].
All geometrical optimizations and single-point electronic energy calculations have been carried out with the LACVP** and LACV3P**++ basis sets, respectively [79,80,81,82]. These basis sets describe the 1s–4d core electrons of the gold atom with the Hay and Wadt core-valence relativistic (i.e., with an implicit treatment of scalar-relativistic effects) effective core-potential (ECP) leaving the outer electrons of gold as well as electrons of the remaining atoms to be treated explicitly by the 6-31G** and 6-311++G** basis sets of double- and triple-ξ quality, correspondingly. Additionally, the basis set aug–cc–pVTZ [83] was used for the MP2 single point energy calculations in order to assess the high-level ab initio methodology performance with a large basis set.
Frequency calculations were performed to verify the correct nature of the stationary points and to estimate zero-point energy (ZPE) and thermal corrections to thermodynamic properties. Intrinsic reaction coordinate (IRC) calculations were employed to locate reagents and products minima connected with the transition states for each considered reaction step.
Single-point electronic energy calculations were performed on the gas-phase geometries. The Poisson–Boltzmann (PB) continuum solvent method was used to describe the solvation. It represents the solvent as a layer of charges at the molecular surface, serving as a dielectric continuum boundary, thus accounting for the detailed molecular shape [84,85]. Solvation energies were calculated on gas-phase stationary points at the B3LYP/LACV3P**++ level of theory, taking a dielectric constant of 80.37 for water and the standard set of optimized radii in Jaguar. Experimental values have been taken for the solvation-free energies of small ions. Thermodynamic properties in solution were calculated as follows. Solvation-free energies were taken as the difference between the solution energies and the gas phase energies. The calculation of entropic contributions in solution is a rather delicate issue, and is particularly important when entropies (and free energies) need to be evaluated for bimolecular reactions involving separated reactants and products. Indeed, in this case, the translational degree of freedom in the reactants/products becomes a loose vibration for adducts or transition states, leading to a loss/gain of entropy, much larger in the gas phase than in the confined condensed phase. For this reason, the use of gas phase entropies in solution often leads to artificially too large results. A way to overcome this problem was proposed by Werz [86] and proved to lead to solution entropies and free energies in excellent agreement with experimental values [87,88,89,90,91]. Following this approach, a solute dissolved in a solvent loses a constant fraction of its entropy in vacuo. The amount of this loss, for water, can be approximated as 50% of the gas phase entropy changed of sign.
Enthalpies and Gibbs free energies were obtained at 298.15 K and 1 atm from total electronic energies including ZPE, thermal, and entropic corrections, and the addition of solvent corrections. A further correction of ΔG°→*=RT ln (24.6) was included to pass from ideal gas at 1 atm to ideal solution at 1 mol/L as standard state [92].
The solvent-accessible surface (SAS) of each considered amino acid was computed by using the sas option implemented in Gromacs software [93].

4. Conclusions

Detailed mechanisms of the reaction of gold carbene compounds with most amino acids and protein fragments remain largely undisclosed. In this study, we have computed the reactivity of IMeAuCl with the various nucleophiles water, side chains of arginine, cysteine, histidine, lysine, methionine, and selenocysteine, and the N-terminal of any amino acid.
The activation-free energy of 23.0 kcal/mol calculated for the chloride substitution by water makes this reaction relatively fast at a physiological temperature, although it shows a barrier significantly higher than all the other nucleophiles examined in the present study, indicating that the [Au(I)(NHC)Cl] complexes are likely to attack the biological targets before being hydrated.
The results for chloride substitution from the considered residues show that all reactions are slightly to moderately exergonic with the following order of thermodynamic stability:
Sec > activated-Cys > Arg > Lys > His > activated-Gln/Asn > N-terminal > Met > Cys > Glu/Asp > Gln/Asn >
but have a different kinetic reactivity order:
Sec > Lys > activated-Gln/Asn > activated- Cys > Cys > Met > Arg > N-terminal > His > Gln/Asn ≈ Glu/Asp >
The lowest barriers calculated for Cys and especially Sec are consistent with the experimental evidence showing the prevalent binding of a [IMeAu]+ fragment to the selenocysteine residue and, to a much lesser extent, the cysteine residue of thioredoxine.
Relatively low barriers are calculated also for the N-based nucleophile sites, lysine, arginine, and N-terminal amine (16–19 kcal/mol), indicating little discrimination in terms of reactivity within either S/Se or N-based nucleophile sites and suggesting that that these N-based sites could also bind an [IMeAu]+ moiety whenever cysteine or selenocysteine residues are absent or buried and that the binding preference to a protein side chain could be determined by its accessibility rather than by its kinetic preference toward chloride substitution.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/4/820/s1.

Author Contributions

Conceptualization, N.R.; methodology, C.C. and A.M.; investigation, I.T.; data curation, A.M.; writing—original draft preparation, I.T.; writing—review and editing, C.C., A.M., and N.R.; supervision, N.R.

Funding

This research received no external funding.

Acknowledgments

Luigi Messori is acknowledged for useful discussion about the experimental behavior of the [Au(NHC)Cl] complexes.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

DFTDensity Functional Theory
PBPoisson–Boltzmann
MP2Møller–Plesset 2nd order perturbation theory
NHCN-Heterocyclic Carbenes

References

  1. Bindoli, A.; Rigobello, M.P.; Scutari, G.; Gabbiani, C.; Casini, A.; Messori, L. Thioredoxin reductase: A target for gold compounds acting as potential anticancer drugs. Coord. Chem. Rev. 2009, 253, 1692–1707. [Google Scholar] [CrossRef]
  2. Milacic, V.; Dou, Q.P. The tumor proteasome as a novel target for gold (III) complexes: Implications for breast cancer therapy. Coord. Chem. Rev. 2009, 253, 1649–1660. [Google Scholar] [CrossRef] [PubMed]
  3. Molter, A.; Mohr, F. Gold complexes containing organoselenium and organotellurium ligands. Coord. Chem. Rev. 2010, 254, 19–45. [Google Scholar] [CrossRef]
  4. Ott, I. On the medicinal chemistry of gold complexes as anticancer drugs. Coord. Chem. Rev. 2009, 253, 1670–1681. [Google Scholar] [CrossRef]
  5. Berners-Price, S.J.; Filipovska, A. Gold compounds as therapeutic agents for human diseases. Metallomics 2011, 3, 863–873. [Google Scholar] [CrossRef] [PubMed]
  6. Nobili, S.; Mini, E.; Landini, I.; Gabbiani, C.; Casini, A.; Messori, L. Gold compounds as anticancer agents: Chemistry, cellular pharmacology, and preclinical studies. Med. Res. Rev. 2010, 30, 550–580. [Google Scholar] [CrossRef] [PubMed]
  7. Che, C.M.; Sun, R.W.Y. Therapeutic applications of gold complexes: Lipophilic gold (III) cations and gold (I) complexes for anti-cancer treatment. Chem. Commun. 2011, 47, 9554–9560. [Google Scholar] [CrossRef]
  8. Shaw, C.F. Gold-based therapeutic agents. Chem. Rev. 1999, 99, 2589–2600. [Google Scholar] [CrossRef]
  9. Gabbiani, C.; Messori, L. Protein targets for anticancer gold compounds: Mechanistic inferences. Anticancer Agents Med. Chem. 2011, 11, 929–939. [Google Scholar] [CrossRef]
  10. Sun, R.W.Y.; Li, C.K.L.; Ma, D.L.; Yan, J.J.; Lok, C.N.; Leung, C.H.; Zhu, N.; Che, C.M. Stable anticancer gold (III)–porphyrin complexes: Effects of porphyrin structure. Chem. Eur. J. 2010, 16, 3097–3113. [Google Scholar] [CrossRef]
  11. Casini, A.; Messori, L. Molecular mechanisms and proposed targets for selected anticancer gold compounds. Curr. Top. Med. Chem. 2011, 11, 2647–2660. [Google Scholar] [CrossRef] [PubMed]
  12. Bruijnincx, P.C.; Sadler, P.J. New trends for metal complexes with anticancer activity. Curr. Opin. Chem. Biol. 2008, 12, 197–206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Hashmi, A.S.K.; Hutchings, G.J. Gold Catalysis. Angew. Chem. Int. Ed. 2006, 45, 7896–7936. [Google Scholar] [CrossRef] [PubMed]
  14. Gorin, D.J.; Sherry, B.D.; Toste, F.D. Ligand effects in homogeneous Au catalysis. Chem. Rev. 2008, 108, 3351–3378. [Google Scholar] [CrossRef] [PubMed]
  15. Li, Z.; Brouwer, C.; He, C. Gold-catalyzed organic transformations. Chem. Rev. 2008, 108, 3239–3265. [Google Scholar] [CrossRef] [PubMed]
  16. Marion, N.; Nolan, S.P. N-heterocyclic carbenes in gold catalysis. Chem. Soc. Rev. 2008, 37, 1776–1782. [Google Scholar] [CrossRef]
  17. Barnard, P.J.; Berners-Price, S.J. Targeting the mitochondrial cell death pathway with gold compounds. Coord. Chem. Rev. 2007, 251, 1889–1902. [Google Scholar] [CrossRef]
  18. Navarro, M. Gold complexes as potential anti-parasitic agents. Coord. Chem. Rev. 2009, 253, 1619–1626. [Google Scholar] [CrossRef]
  19. Krishnamurthy, D.; Karver, M.R.; Fiorillo, E.; Orrú, V.; Stanford, S.M.; Bottini, N.; Barrios, M. Gold(I)-mediated inhibition of protein tyrosine phosphatases: A detailed in vitro and cellular study. J. Med. Chem. 2008, 51, 4790–4795. [Google Scholar] [CrossRef]
  20. Glišić, B.D.; Djuran, M.I. Gold complexes as antimicrobial agents: An overview of different biological activities in relation to the oxidation state of the gold ion and the ligand structure. Dalton Trans. 2014, 43, 5950–5969. [Google Scholar] [CrossRef]
  21. Liu, W.; Gust, R. Metal N-heterocyclic carbene complexes as potential antitumor metallodrugs. Chem. Soc. Rev. 2013, 42, 755–773. [Google Scholar] [CrossRef] [PubMed]
  22. Rubbiani, R.; Kitanovic, I.; Alborzinia, H.; Can, S.; Kitanovic, A.; Onambele, L.A.; Stefanopoulou, M.; Geldmacher, Y.; Sheldrick, W.S.; Wolber, G. Benzimidazol-2-ylidene gold(I) complexes are thioredoxin reductase inhibitors with multiple antitumor properties. J. Med. Chem. 2010, 53, 8608–8618. [Google Scholar] [CrossRef] [PubMed]
  23. Oehninger, L.; Rubbiani, R.; Ott, I. N-Heterocyclic carbene metal complexes in medicinal chemistry. Dalton Trans. 2013, 42, 3269–3284. [Google Scholar] [CrossRef] [PubMed]
  24. Weaver, J.; Gaillard, S.; Toye, C.; Macpherson, S.; Nolan, S.P.; Riches, A. Cytotoxicity of gold(I) N-heterocyclic carbene complexes assessed by using human tumor cell lines. Chem. Eur. J. 2011, 17, 6620–6624. [Google Scholar] [CrossRef] [PubMed]
  25. Hickey, J.J.; Ruhayel, R.A.; Barnard, P.J.; Baker, M.V.; Berners-Price, S.J.; Filipovska, A. Mitochondria-targeted chemotherapeutics: The rational design of gold(I) N-heterocyclic carbene complexes that are selectively toxic to cancer cells and target protein selenols in preference to thiols. J. Am. Chem. Soc. 2008, 130, 12570–12571. [Google Scholar] [CrossRef] [PubMed]
  26. Gasser, G.; Ott, I.; Metzler-Nolte, N. Organometallic anticancer compounds. J. Med. Chem. 2011, 54, 3–25. [Google Scholar] [CrossRef] [PubMed]
  27. Bertrand, B.; Casini, A. A golden future in medicinal inorganic chemistry: The promise of anticancer gold organometallic compounds. Dalton Trans. 2014, 43, 4209–4219. [Google Scholar] [CrossRef]
  28. Schuh, E.; Pflüger, C.; Citta, A.; Folda, A.; Rigobello, M.P.; Bindoli, A.; Casini, A.; Mohr, F. Gold (I) carbene complexes causing thioredoxin 1 and thioredoxin 2 oxidation as potential anticancer agents. J. Med. Chem. 2012, 55, 5518–5528. [Google Scholar] [CrossRef]
  29. Arambula, J.F.; McCall, R.; Sidoran, K.J.; Magda, D.; Mitchell, N.A.; Bielawski, C.W.; Lynch, V.M.; Sessler, J.L.; Arumugam, K. Gold(I) carbene complexes causing thioredoxin 1 and thioredoxin 2 oxidation as potential anticancer agents. Chem. Sci. 2016, 7, 1245–1256. [Google Scholar] [CrossRef]
  30. Pratesi, A.; Gabbiani, C.; Michelucci, E.; Ginanneschi, M.; Papini, A.M.; Rubbiani, R.; Ott, I.; Messori, L. Insights on the mechanism of thioredoxin reductase inhibition by gold N-heterocyclic carbene compounds using the synthetic linear selenocysteine containing C-terminal peptide hTrxR(488–499): An ESI-MS investigation. J. Inorg. Biochem. 2014, 136, 161–169. [Google Scholar] [CrossRef]
  31. Arcau, J.; Andermark, V.; Rodrigues, M.; Giannicchi, I.; PérezGarcia, L.; Ott, I.; Rodríguez, L. Synthesis and Biological Activity of Gold(I) N-Heterocyclic Carbene Complexes with Long Aliphatic Side Chains. Eur. J. Inorg. Chem. 2014, 35, 6117–6125. [Google Scholar] [CrossRef]
  32. Urig, S.; Fritz-Wolf, K.; Réau, R.; Herold-Mende, C.; Tóth, K.; Davioud-Charvet, E.; Becker, K. Undressing of Phosphine Gold(I) Complexes as Irreversible Inhibitors of Human Disulfide Reductases. Angew. Chem. Int. Ed. 2006, 45, 1881–1886. [Google Scholar] [CrossRef] [PubMed]
  33. Messori, L.; Scaletti, F.; Massai, L.; Cinellu, M.A.; Russo Krauss, I.; Di Martino, G.; Vergara, A.; Paduano, L.; Merlino, A. Interactions of gold-based drugs with proteins: Crystal structure of the adduct formed between ribonuclease A and a cytotoxic gold (III) compound. Metallomics 2014, 6, 233–236. [Google Scholar] [CrossRef] [PubMed]
  34. Messori, L.; Scaletti, F.; Massai, L.; Cinellu, M.A.; Gabbiani, C.; Vergara, A.; Merlino, A. The mode of action of anticancer gold-based drugs: A structural perspective. Chem. Commun. 2013, 49, 10100–10102. [Google Scholar]
  35. Ferraro, G.; Massai, L.; Messori, L.; Cinellu, M.A.; Merlino, A. Structural evidences for a secondary gold binding site in the hydrophobic box of lysozyme. Biometals 2015, 28, 745–754. [Google Scholar] [CrossRef] [PubMed]
  36. Ferraro, G.; Gabbiani, C.; Merlino, A. First crystal structure for a gold carbene–protein adduct. Bioconjugate Chem. 2016, 27, 1584–1587. [Google Scholar] [CrossRef] [PubMed]
  37. Gabbiani, C.; Scaletti, F.; Massai, L.; Michelucci, E.; Cinellu, M.A.; Messori, L. Medicinal gold compounds form tight adducts with the copper chaperone Atox-1: Biological and pharmacological implications. Chem. Commun. 2012, 48, 11623–11625. [Google Scholar] [CrossRef] [PubMed]
  38. Bhabak, K.P.; Bhuyan, B.J.; Mugesh, G. Bioinorganic and medicinal chemistry: Aspects of gold(I)-protein complexes. Dalton Trans. 2011, 40, 2099–2111. [Google Scholar] [CrossRef] [PubMed]
  39. Serebryanskaya, T.V.; Lyakhov, A.S.; Ivashkevich, L.S.; Schur, J.; Frias, C.; Prokop, A.; Ott, I. Gold (I) thiotetrazolates as thioredoxin reductase inhibitors and antiproliferative agents. Dalton Trans. 2015, 44, 1161–1169. [Google Scholar] [CrossRef]
  40. Vergara, E.; Casini, A.; Sorrentino, F.; Zava, O.; Cerrada, E.; Rigobello, M.P.; Bindoli, A.; Laguna, M.; Dyson, P.J. Anticancer therapeutics that target selenoenzymes: Synthesis, characterization, in vitro cytotoxicity, and thioredoxin reductase inhibition of a series of gold(I) complexes containing hydrophilic phosphine ligands. ChemMedChem 2010, 5, 96–102. [Google Scholar] [CrossRef]
  41. Ott, I.; Qian, X.; Xu, Y.; Vlecken, D.H.W.; Marques, I.J.; Kubutat, D.; Will, J.; Sheldrick, W.S.; Jesse, P.; Prokop, A. A gold(I) phosphine complex containing a naphthalimide ligand functions as a TrxR inhibiting antiproliferative agent and angiogenesis inhibitor. J. Med. Chem. 2009, 52, 763–770. [Google Scholar] [CrossRef] [PubMed]
  42. Baker, M.V.; Barnard, P.J.; Berners-Price, S.J.; Brayshaw, S.K.; Hickey, J.L.; Skelton, B.W.; White, A.H. Cationic, linear Au(I) N-heterocyclic carbene complexes: Synthesis, structure and anti-mitochondrial activity. Dalton Trans. 2006, 30, 3708–3715. [Google Scholar] [CrossRef] [PubMed]
  43. Aher, S.B.; Muskawar, P.N.; Thenmozhi, K.; Bhagat, P.R. Recent developments of metal N-heterocyclic carbenes as anticancer agents. Eur. J. Med. Chem. 2014, 81, 408–419. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, J.J.; Che, C.M.; Ott, I. Caffeine derived platinum(II) N-heterocyclic carbene complexes with multiple anti-cancer activities. J. Organomet. Chem. 2015, 782, 37–41. [Google Scholar] [CrossRef]
  45. Wang, H.M.; Vasam, C.S.; Tsai, T.Y.; Chen, S.H.; Chang AH, H.; Lin, I.J. Gold (I) N-heterocyclic carbene and carbazolate complexes. Organometallics 2005, 24, 486–493. [Google Scholar] [CrossRef]
  46. Messori, L.; Cinellu, M.A.; Merlino, A. Protein recognition of gold-based drugs: 3D structure of the complex formed when lysozyme reacts with Aubipyc. ACS Med. Chem. Lett. 2014, 5, 1110–1113. [Google Scholar] [CrossRef] [PubMed]
  47. Messori, L.; Marchetti, L.; Massai, L.; Scaletti, F.; Guerri, A.; Landini, I.; Nobili, S.; Perrone, G.; Mini, E.; Leoni, P.; et al. Chemistry and biology of two novel gold (I) carbene complexes as prospective anticancer agents. Inorg. Chem. 2014, 53, 2396–2403. [Google Scholar] [CrossRef]
  48. Rubbiani, R. Gold (I) N-Heterocyclic Carbene Complexes: A Chemical and Biological Study of Their Therapeutic Potential as New Anticancer Agents. Ph.D. Thesis, TU Braunschweig, Braunschweig, Germany, 2013. [Google Scholar]
  49. Minenkov, Y.; Singstad, Å.; Occhipinti, G.; Jensen, V.R. The accuracy of DFT-optimized geometries of functional transition metal compounds: A validation study of catalysts for olefin metathesis and other reactions in the homogeneous phase. Dalton Trans. 2012, 41, 5526–5541. [Google Scholar] [CrossRef]
  50. Qi, S.C.; Hayashi, J.I.; Zhang, L. Recent application of calculations of metal complexes based on density functional theory. RSC Adv. 2016, 6, 77375–77395. [Google Scholar] [CrossRef]
  51. Mardirossian, N.; Head-Gordon, M. How accurate are the Minnesota density functionals for noncovalent interactions, isomerization energies, thermochemistry, and barrier heights involving molecules composed of main-group elements? J. Chem. Theory Comput. 2016, 12, 4303–4325. [Google Scholar] [CrossRef]
  52. Grigorenko, B.R.; Khrenova, M.G.; Nemukin, A.V. Amide–imide tautomerization in the glutamine side chain in enzymatic and photochemical reactions in proteins. Phys. Chem. Chem. Phys. 2018, 20, 23827–23836. [Google Scholar] [PubMed]
  53. Zou, J.; Taylor, P.; Dornan, J.; Robinson, S.P.; Walkinshaw, M.D.; Sadler, P.J. First crystal structure of a medicinally relevant gold protein complex: Unexpected binding of [Au(PEt3)]+ to histidine. Angew. Chem. Int. Ed. Engl. 2000, 39, 2931–2934. [Google Scholar] [CrossRef]
  54. Bochevarov, A.D.; Harder, E.; Hughes, T.F.; Greenwood, J.R.; Braden, D.A.; Philipp, D.M.; Rinaldo, D.; Halls, M.D.; Zhang, J.; Friesner, R.A. Jaguar: A high-performance quantum chemistry software program with strengths in life and materials sciences. Int. J. Quantum Chem. 2013, 113, 2110–2142. [Google Scholar] [CrossRef] [Green Version]
  55. Jaguar Version 7.5; Schrodinger LLC: New York, NY, USA, 2008.
  56. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. (Eds.) Gaussian 09, Revision A.02; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  57. Becke, A.D. Becke’s three parameter hybrid method using the LYP correlation functional. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  58. Stephens, P.J.; Devlin, F.J.; Chabalowski, C.; Frisch, M.J. Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. J. Phys. Chem. 1994, 98, 11623–11627. [Google Scholar] [CrossRef]
  59. Niu, S.; Hall, B.M. Theoretical studies on reactions of transition-metal complexes. Chem. Rev. 2000, 100, 353–405. [Google Scholar] [CrossRef] [PubMed]
  60. Tolbatov, I.; Coletti, C.; Marrone, A.; Re, N. Insight into the electrochemical reduction mechanism of Pt (IV) anticancer complexes. Inorg. Chem. 2018, 57, 3411–3419. [Google Scholar] [CrossRef] [PubMed]
  61. Paciotti, R.; Corinti, D.; De Petris, A.; Ciavardini, A.; Piccirillo, S.; Coletti, C.; Re, N.; Maitre, P.; Bellina, B.; Barran, P.; et al. Cisplatin and transplatin interaction with methionine: Bonding motifs assayed by vibrational spectroscopy in the isolated ionic complexes. Phys. Chem. Chem. Phys. 2017, 19, 26697–26707. [Google Scholar] [PubMed]
  62. Corinti, D.; De Petris, A.; Coletti, C.; Re, N.; Chiavarino, B.; Crestoni, M.E.; Fornarini, S. Cisplatin primary complex with L-histidine target revealed by IR multiple photon dissociation (IRMPD) spectroscopy. Chem. Phys. Chem. 2017, 18, 318–325. [Google Scholar] [CrossRef] [PubMed]
  63. Corinti, D.; Coletti, C.; Re, N.; Chiavarino, B.; Crestoni, M.E.; Fornarini, S. Cisplatin binding to biological ligands revealed at the encounter complex level by IR action spectroscopy. Chem. Eur. J. 2016, 22, 3794–3803. [Google Scholar] [CrossRef] [PubMed]
  64. Graziani, V.; Coletti, C.; Marrone, A.; Re, N. Activation and reactivity of a bispidine analogue of cisplatin: A theoretical investigation. J. Phys. Chem. A 2016, 120, 5175–5186. [Google Scholar] [CrossRef] [PubMed]
  65. Ciancetta, A.; Coletti, C.; Marrone, A.; Re, N. Activation of carboplatin by carbonate: A theoretical investigation. Dalton Trans. 2012, 41, 12960–12969. [Google Scholar] [CrossRef] [PubMed]
  66. Ciancetta, A.; Coletti, C.; Marrone, A.; Re, N. Activation of carboplatin by chloride ions: A theoretical investigation. Theor. Chem. Acc. 2011, 129, 757–769. [Google Scholar] [CrossRef]
  67. Besker, N.; Coletti, C.; Marrone, A.; Re, N. Binding of antitumor ruthenium complexes to DNA and proteins: A theoretical approach. J. Phys. Chem. B 2007, 111, 9955–9964. [Google Scholar] [CrossRef] [PubMed]
  68. Paciotti, R.; Tolbatov, I.; Graziani, V.; Marrone, A.; Re, N.; Coletti, C. Insights on the Activity of Platinum-based Anticancer Complexes through Computational Methods. AIP Conf. Proc. 2018, 2040, 020019. [Google Scholar]
  69. Marino, T.; Russo, N.; Toscano, M.; Pavelka, M. Theoretical investigation on DNA/RNA base pairs mediated by copper, silver, and gold cations. Dalton Trans. 2012, 41, 1816–1823. [Google Scholar] [CrossRef] [PubMed]
  70. Nuevo, D.; Poyatos, M.; Peris, E. A dinuclear Au(I) complex with a pyrene-di-N-heterocyclic carbene linker: Supramolecular and catalytic studies. Organometallics 2018, 37, 3407–3411. [Google Scholar] [CrossRef]
  71. Graziani, V.; Marrone, A.; Re, N.; Coletti, C.; Platts, J.A.; Casini, A. A multi-level theoretical study to disclose the binding mechanisms of gold (III)–bipyridyl compounds as selective aquaglyceroporin inhibitors. Chem. Eur. J. 2017, 23, 13802–13813. [Google Scholar] [CrossRef]
  72. Wenzel, M.; Mósca, A.; Graziani, V.; Aikman, B.; Thomas, S.; de Almeida, A.; Platts, J.; Re, N.; Coletti, C.; Marrone, A.; et al. Insights into the Mechanisms of Aquaporin-3 Inhibition by Gold(III) Complexes: The Importance of Non-Coordinative Adduct Formation. Inorg. Chem. 2019, 58, 2140–2148. [Google Scholar] [CrossRef]
  73. Zhao, Y.; Truhlar, D. 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. [Google Scholar]
  74. Chai, J.D.; Head-Gordon, M. Systematic optimization of long-range corrected hybrid density functionals. J. Chem. Phys. 2008, 128, 084106. [Google Scholar] [CrossRef] [PubMed]
  75. Yanai, T.; Tew, D.P.; Handy, N.C. A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP). Chem. Phys. Lett. 2004, 393, 51–57. [Google Scholar] [CrossRef] [Green Version]
  76. Frisch, M.J.; Head-Gordon, M.; Pople, J.A. Direct MP2 gradient method. Chem. Phys. Lett. 1990, 166, 275–280. [Google Scholar] [CrossRef]
  77. Frisch, M.J.; Head-Gordon, M.; Pople, J.A. Semi-direct algorithms for the MP2 energy and gradient. Chem. Phys. Lett. 1990, 166, 281–289. [Google Scholar] [CrossRef]
  78. Head-Gordon, M.; Pople, J.A.; Frisch, M.J. MP2 energy evaluation by direct methods. Chem. Phys. Lett. 1988, 153, 503–506. [Google Scholar] [CrossRef]
  79. Hay, P.J.; Wadt, W.R. 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. [Google Scholar] [CrossRef]
  80. Ditchfield, R.; Hehre, W.J.; Pople, J.A. Self-consistent molecular orbital methods. 9. Extended Gaussian-type basis for molecular-orbital studies of organic molecules. J. Chem. Phys. 1971, 54, 724. [Google Scholar] [CrossRef]
  81. McLean, A.D.; Chandler, G.S. Contracted Gaussian basis sets for molecular calculations. I. Second row atoms, Z= 11–18. J. Chem. Phys. 1980, 72, 5639–5648. [Google Scholar] [CrossRef]
  82. Krishnan, R.B.; Binkley, J.S.; Seeger, R.; Pople, J.A. Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions. J. Chem. Phys. 1980, 72, 650–654. [Google Scholar] [CrossRef]
  83. Neto, A.C.; Muniz, E.P.; Centoducatte, R.; Jorge, F.E. Gaussian basis sets for correlated wave functions. Hydrogen, helium, first-and second-row atoms. J. Mol. Struct. 2005, 718, 219–224. [Google Scholar] [CrossRef]
  84. Tannor, D.J.; Marten, B.; Murphy, R.; Friesner, R.A.; Sitkoff, D.; Nicholls, A.; Ringnalda, M.; Goddard, W.A., III; Honig, B. Accurate first principles calculation of molecular charge distributions and solvation energies from ab initio quantum mechanics and continuum dielectric theory. J. Am. Chem. Soc. 1994, 116, 11875–11882. [Google Scholar] [CrossRef]
  85. Marten, B.; Kim, K.; Cortis, C.; Friesner, R.A.; Murphy, R.B.; Ringnalda, M.N.; Sitkoff, D.; Honig, B. New model for calculation of solvation free energies: Correction of self-consistent reaction field continuum dielectric theory for short-range hydrogen-bonding effects. J. Phys. Chem. 1996, 100, 11775–11788. [Google Scholar] [CrossRef]
  86. Wertz, D.H. Relationship between the gas-phase entropies of molecules and their entropies of solvation in water and 1-octanol. J. Am. Chem. Soc. 1980, 102, 5316–5322. [Google Scholar] [CrossRef]
  87. Cooper, J.; Ziegler, T. A density functional study of SN2 substitution at square-planar platinum (II) complexes. Inorg. Chem. 2002, 41, 6614–6622. [Google Scholar] [CrossRef]
  88. Lau, J.K.C.; Deubel, D.V. Hydrolysis of the anticancer drug cisplatin: Pitfalls in the interpretation of quantum chemical calculations. J. Chem. Theory Comput. 2006, 2, 103–106. [Google Scholar] [CrossRef] [PubMed]
  89. Ahlquist, M.; Nielsen, R.J.; Periana, R.A.; Goddard, W.A., III. Product protection, the key to developing high performance methane selective oxidation catalysts. J. Am. Chem. Soc. 2009, 131, 17110–17115. [Google Scholar] [CrossRef]
  90. Zhu, H.; Ziegler, T. Probing the influence of trans and leaving ligands on the ability of square-planar platinum (II) complexes to activate methane. A theoretical study. Organometallics 2009, 28, 2773–2777. [Google Scholar] [CrossRef]
  91. Bercaw, J.E.; Chen, G.S.; Labinger, J.A.; Lin, B.L. Protonolysis of platinum (II) and palladium (II) methyl complexes: A combined experimental and theoretical investigation. Organometallics 2010, 29, 4354–4359. [Google Scholar] [CrossRef]
  92. Ben-Naim, A.Y. Solvation Thermodynamics; Springer Science & Business Media: New York, NY, USA, 2013. [Google Scholar]
  93. Berendsen, H.J.C.; van der Spoel, D.; van Drunen, R. GROMACS: A message-passing parallel molecular dynamics implementation. Comp. Phys. Comm. 1995, 91, 43–56. [Google Scholar] [CrossRef] [Green Version]
Ijms 20 00820 sch001 550
Scheme 1. Reaction scheme and considered nucleophiles. Attacking atoms in bold.
Scheme 1. Reaction scheme and considered nucleophiles. Attacking atoms in bold.
Ijms 20 00820 sch001
Ijms 20 00820 g001 550
Figure 1. Calculated geometrical structure of the IMeAuCl complex. Bond lengths (purple in the figure) in Å, angles (green in the figure) in degrees.
Figure 1. Calculated geometrical structure of the IMeAuCl complex. Bond lengths (purple in the figure) in Å, angles (green in the figure) in degrees.
Ijms 20 00820 g001
Ijms 20 00820 g002 550
Figure 2. Transition states for selected nucleophiles. All distances in angstroms, all angles in degrees.
Figure 2. Transition states for selected nucleophiles. All distances in angstroms, all angles in degrees.
Ijms 20 00820 g002
Ijms 20 00820 g003 550
Figure 3. Reaction profiles for IMeAuCl interacting with selected nucleophiles. Values in kcal/mol are calculated at the B3LYP/LACV3P**++//B3LYP/LACVP** level. R, P, TS, RA, and PA stand for reactant, product, transition state, reactant-adduct, and product-adduct, respectively.
Figure 3. Reaction profiles for IMeAuCl interacting with selected nucleophiles. Values in kcal/mol are calculated at the B3LYP/LACV3P**++//B3LYP/LACVP** level. R, P, TS, RA, and PA stand for reactant, product, transition state, reactant-adduct, and product-adduct, respectively.
Ijms 20 00820 g003
Ijms 20 00820 sch002 550
Scheme 2. Amidic/imidic tautomer equilibrium of glutamine/asparagine.
Scheme 2. Amidic/imidic tautomer equilibrium of glutamine/asparagine.
Ijms 20 00820 sch002
Ijms 20 00820 g004 550
Figure 4. Solvent accessible surface (SAS) calculated on an X-ray structure of Thaumatin-1 (pdb entry 3qy5). Arrows indicate metal bound amino acids. Stars indicate metal centers with an attached carbene ligand.
Figure 4. Solvent accessible surface (SAS) calculated on an X-ray structure of Thaumatin-1 (pdb entry 3qy5). Arrows indicate metal bound amino acids. Stars indicate metal centers with an attached carbene ligand.
Ijms 20 00820 g004
Table
Table 1. Enthalpies/Gibbs free energies for reactions with various nucleophiles. Reported values are calculated at the B3LYP/LACV3P**++//B3LYP/LACVP** level and reported in kcal/mol in format Enthalpy (ΔH)/Gibbs free energy (ΔG) (in solution). R, P, TS, RA, and PA stand for reactant, product, transition state, reactant-adduct, and product-adduct, respectively.
Table 1. Enthalpies/Gibbs free energies for reactions with various nucleophiles. Reported values are calculated at the B3LYP/LACV3P**++//B3LYP/LACVP** level and reported in kcal/mol in format Enthalpy (ΔH)/Gibbs free energy (ΔG) (in solution). R, P, TS, RA, and PA stand for reactant, product, transition state, reactant-adduct, and product-adduct, respectively.
Reaction withR → PR → TS
ΔHΔGΔHΔG
H2O7.21.316.823.0
Ala N-terminal−6.6−11.211.918.8
Arg−11.0−16.111.318.3
Cys−3.5−9.110.918.1
Cys–−16.8−21.411.617.4
His−8.3−13.414.120.6
Lys−10.2−14.49.016.2
Met−5.4−9.811.518.2
Sec−2.2−6.711.217.6
Sec–−19.2−23.86.712.8
Gln/Asn (amidic form)6.02.417.524.8
Gln/Asn (imidic form) *−7.3−12.610.016.8
Glu/Asp acid−2.2−7.018.224.7
* Reactant for the imidic form of Gln/Asn is 16.4 kcal/mol higher than its amidic form, so all the values for the imidic form should be augmented by this value before its comparison with the amidic form. Additionally, see Figure S1.

Share and Cite

MDPI and ACS Style

Tolbatov, I.; Coletti, C.; Marrone, A.; Re, N. Reactivity of Gold(I) Monocarbene Complexes with Protein Targets: A Theoretical Study. Int. J. Mol. Sci. 2019, 20, 820. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20040820

AMA Style

Tolbatov I, Coletti C, Marrone A, Re N. Reactivity of Gold(I) Monocarbene Complexes with Protein Targets: A Theoretical Study. International Journal of Molecular Sciences. 2019; 20(4):820. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20040820

Chicago/Turabian Style

Tolbatov, Iogann, Cecilia Coletti, Alessandro Marrone, and Nazzareno Re. 2019. "Reactivity of Gold(I) Monocarbene Complexes with Protein Targets: A Theoretical Study" International Journal of Molecular Sciences 20, no. 4: 820. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20040820

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop