Next Article in Journal
In Search of Small Molecules That Selectively Inhibit MBOAT4
Next Article in Special Issue
A Comprehensive Review on Metal Catalysts for the Production of Cyclopentanone Derivatives from Furfural and HMF
Previous Article in Journal
Scientometric Overview of Coffee By-Products and Their Applications
Previous Article in Special Issue
Organic Phase Change Materials for Thermal Energy Storage: Influence of Molecular Structure on Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Computational Studies of Au(I) and Au(III) Anticancer MetalLodrugs: A Survey

1
Institut de Chimie Moleculaire de l’Université de Bourgogne (ICMUB), Université de Bourgogne Franche-Comté (UBFC), Avenue Alain Savary 9, 21078 Dijon, France
2
Dipartimento di Farmacia, Università degli Studi “G. D’Annunzio” Chieti-Pescara, Via dei Vestini, 66100 Chieti, Italy
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(24), 7600; https://doi.org/10.3390/molecules26247600
Submission received: 19 October 2021 / Revised: 29 November 2021 / Accepted: 12 December 2021 / Published: 15 December 2021
(This article belongs to the Special Issue Feature Review Papers in Physical Chemistry)

Abstract

:
Owing to the growing hardware capabilities and the enhancing efficacy of computational methodologies, computational chemistry approaches have constantly become more important in the development of novel anticancer metallodrugs. Besides traditional Pt-based drugs, inorganic and organometallic complexes of other transition metals are showing increasing potential in the treatment of cancer. Among them, Au(I)- and Au(III)-based compounds are promising candidates due to the strong affinity of Au(I) cations to cysteine and selenocysteine side chains of the protein residues and to Au(III) complexes being more labile and prone to the reduction to either Au(I) or Au(0) in the physiological milieu. A correct prediction of metal complexes’ properties and of their bonding interactions with potential ligands requires QM computations, usually at the ab initio or DFT level. However, MM, MD, and docking approaches can also give useful information on their binding site on large biomolecular targets, such as proteins or DNA, provided a careful parametrization of the metal force field is employed. In this review, we provide an overview of the recent computational studies of Au(I) and Au(III) antitumor compounds and of their interactions with biomolecular targets, such as sulfur- and selenium-containing enzymes, like glutathione reductases, glutathione peroxidase, glutathione-S-transferase, cysteine protease, thioredoxin reductase and poly (ADP-ribose) polymerase 1.

1. Introduction

The earliest recorded use of gold in medicine is traced back to 2500 BC when it was used in ancient China; much later, in seventeenth-century Europe, it was used against melancholy, fainting, fevers, and against syphilis in the late nineteenth century [1]. Gold complexes are ubiquitously utilized for the treatment of severe rheumatoid arthritis [2]. Moreover, the interest in gold(I) and gold(III) complexes in anticancer therapy has been recently growing owing to the strong antitumor activity in vitro and in vivo as well as to the relative chemical stability easily modulated by the addition of various ligands [3,4].
Computational chemistry nowadays plays an important role in the search for new drug candidates due to the increasing efficacy of methods based on quantitative structure-activity relationships [5]. These methods yield reliable statistical models which permit establishing a correlation between chemical properties of complexes and their biological activities, thus giving insights into their therapeutic mode of action [6]. This methodology is indispensable for the development of chemical compounds for medical use since it allows to find relationships between their affinities toward binding to their expected biotargets and their structural properties. This approach is based on scanning the databases of chemical complexes versus the intended targets, such a tactic being predominant for organic compounds. However, it is much less efficient in the case of metal complexes due to a specific activation path usual for all metallodrugs, relying on the ligand substitution by the biotarget and their ability to operate against a plethora of biotargets owing to the promiscuity of metallodrugs [7,8,9]. This results in the widespread usage of ab initio and DFT methods for studies of metal-based scaffolds intended for medical use [10,11,12,13,14].
In this work, we review recent computational studies of Au(I) and Au(III) antitumor compounds, including joint experimental/computational investigations. We consider this topic to be crucial for the design of novel gold-based metallodrugs and for the further development of existing metalloscaffolds for clinical use.

2. Computational Studies of Au(I) Metallodrugs

The mechanism of action of Au(I) complexes is not yet elucidated but several investigations have suggested that they may trigger apoptosis via inhibition of selenium- and sulfur-containing enzymes such as glutathione reductases, glutathione peroxidase, glutathione-S-transferase, cysteine protease, thioredoxin reductase (TrxR) and poly (ADP-ribose) polymerase 1 (PARP-1) [2,3,4].
Indeed, metal ion Au(I) is a soft metal center, demonstrating a strongly pronounced affinity toward soft ligands such as thiols of cysteines and thioethers of methionines [15] and even more with selenols of selenocysteine residues [16]. Nevertheless, several X-ray crystallographic studies proved that Au(I) ions are able to bind solvent-exposed His [17,18] even when free thiols are available. Moreover, Au(I) complexes may bind the Arg and Lys side chains as well as the N-terminal of Ala [19], albeit in the absence of available cysteines, methionines, and histidines. For instance, binding of Au(NHC)Cl (with NHC = 1-butyl-3-methyl-imidazole-2-ylidene) to the model protein thaumatin occurs at lysine side chains and at the N-terminal tail; the metal binds the protein after releasing Cl ligand, but retaining NHC fragment (Figure 1) [19]. Linear Au(I) complexes are robust inhibitors of the Se-free enzyme glutathione reductase (GR) due to the high selectivity toward thiols [20,21].
Among the most common Au(I) compounds synthesized and tested as potential anticancer drugs are those bearing phosphine [22], thiosugar [23], N-heterocyclic carbenes (NHC)[24], alkynyl [25], and other sulfur-based ligands such as thiosemicarbazone [26] (Figure 2).

2.1. Auranofin and Auranofin Analogs

The most prominent gold-based metallodrug is auranofin (2,3,4,6-tetra-O-acetyl-L-thio-β-D-glyco-pyranosato-S-(triethyl-phosphine)-gold(I) (AF) (Figure 3, structure 1) [23]. Initially known as an anti-arthritic agent, it was later revealed as an antiviral, antibacterial, antiparasitic, and antifungal therapeutics [27,28,29]. Additionally, the antitumoral activity of this agent makes it efficacious in inducing apoptosis in various types of human cancer cells: prostate, ovarian, lung, breast, blood, and bone [30,31]. Its mode of action greatly differs from the platinum-based complexes, based on DNA binding, since it selectively targets sulfur- and selenium-containing proteins. For example, AF extensively and quickly binds albumin [32], proteasome system [33], the NF-κB protein complex [34], and thioredoxin reductase (TrxR) [35], all of which are incorporated in defining the anticancer activity. However, TrxR is believed to be the main target of AF since the gold fragment of AF coordinates tightly to the redox-active site of TrxR, therefore causing oxidative stress and ultimately resulting in cellular apoptosis [36]. However, despite numerous experimental and computational studies, the modus operandi of AF is still elusive at the biomolecular level. There is a ubiquitous agreement that the anticancer effect is produced by the [Au(PEt3)]+ cation whilst the thiosugar operates principally as a carrier ligand [37]. In the last decade, mechanistic consequences of the replacement of thiosugar moiety with different ligands have been studied, and it has been noted that the biological anticancer activity of the [Au(PEt3)]+ fragment is not affected by this alteration [38].
The reactivity of the [Au(PMe3)]+ fragment of auranofin with the cysteine, selenocysteine was studied at DFT level of theory, as well as its interaction with the solvent-exposed subunit of TrxR, a tetrapeptide H2NGlyCysXGlyCOOH (X = Cys, Sec) which is responsible for the catalysis of the NADPH-dependent reduction of thioredoxins [39]. These calculations displayed a higher acidity of Se-H with respect to S-H, causing a stronger binding of gold at the selenium site compared to sulfur. Moreover, it was shown that the reducing capability of H2NGlyCysSecGlyCOOH increases after deprotonation of Sec or diminishes in case of gold coordination at Sec.
To shed light on the mechanism behind the inhibition of TrxR, a computational investigation at DFT level analyzed the ligand exchange reactions of auranofin with the potential target amino acids Cys, Sec, His, and Lys [40]. Both the cleavage of thioglucose and triethylphosphine were studied in aqueous solution at physiological conditions. The thioglucose substitution was established to be kinetically more favorable with the activation barriers in the order Lys < His < Sec < Cys yet the reaction free energies following the trend Sec << Cys<< Lys < His, as expected from the soft nature of gold(I) metal center and the softness of Sec and Cys ligands. The structural analysis of transition states permitted to correlate the reduction of activation enthalpy with the increase in reactive angle, thus envisioning the possibility to modulate auranofin reactivity by adding bulky ligands on the thiosugar moiety.
Another combined NMR and computational investigation at DFT level of the reactions between auranofin and models of thiol and selenol nucleophiles present in TrxR, conducted in chloroform and methanol, revealed important differences related to the polarity of the reaction milieu [41]. It was found that the thiosugar moiety was reversibly substituted by both Cys and Sec, albeit with large differences in the equilibrium constants, being ca. 1 and more than 1000 for the binding at S and Se, respectively. In the polar methanol, the reaction was demonstrated to be more complex, with the phosphine group being able to participate in the ligand exchange. Similar conclusions were made on the reactivity of Et3PAuCl (Figure 3, structure 2).
Computational studies on auranofin and the analogous chlorotriethylphosphine gold complex by means of hybrid density functional theory at the B3LYP/DZVP level [42] revealed metal binding to cysteine and selenocysteine and showed that auranofin could bind selenium in glutathione peroxidase enzymes. These results are corroborated by a mass spectrometry investigation [42], demonstrating that substitution of thiosugar in auranofin with chloride does not alter its reactivity and proving the gold triethylphosphine moiety to be the active ligand.
The iodide analog of auranofin, Au(PEt3)I (AF-I), was found to have an exceptionally advantageous biochemical profile demonstrating powerful cytotoxic activity in vitro against several cancer types and an almost full and prompt remission in an orthotopic in vivo mouse model of ovarian cancer (Figure 3, structure 3) [38]. Furthermore, the lack of severe toxic side effects leads to excellent toleration of AF-I even when high doses are administered. This was the motivation behind a recent combined experimental/theoretical investigation of the mechanistic aspects of the antitumor activity of AF-I, where its reactivity toward molecular models of relevant aminoacidic residues, namely histidine, cysteine, methionine, and selenocysteine, was analyzed, and the enhanced antitumoral activity of AF-I compared to AF was shown [43]. The thermodynamic and kinetic computational results corroborate the conclusion that only Cys and Sec in their anionic forms can favorably substitute thiosugar and iodide in AF and AF-I, respectively, with activation energy for iodide substitution in AF-I being substantially lower [43] than the thiosugar one in AF. This unequivocal preference toward Cys and Sec residues reinforces the understanding that thioredoxin reductase is the most probable target of these gold complexes. Moreover, the marked reactivity of AF-I toward these nucleophile sites makes this Au(I) complex a promising scaffold for the further exploration of its anticancer properties in vivo.
Although auranofin is currently widely studied as an anticancer drug, its original use is the treatment of rheumatoid arthritis [44]. A joint experimental and computational study on the effect of different decoration of the auranofin phosphine ligand towards its in vitro inhibition of cathepsin B, a highly homologous lysosomal cysteine protease participating in the processing of antigens, revealed that sequential replacement of ethyl by phenyl moieties causes a drastic increase in efficacy [45]. It was shown that while AF approaches at cathepsin B by fitting only one hydrophobic site, the resulting Au(I) complexes, endowed with bulky phenyl substituents, may fit two hydrophobic pockets of the enzyme, thus enhancing the binding of this metallodrug at its protein target (Figure 3, structure 4).
There is a body of experimental evidence on the efficacy of auranofin as an antiparasitic and anti-Schistosoma agent, yet the molecular mechanisms behind its inhibition of Schistosoma mansoni, the pathogenic parasite that causes intestinal schistosomiasis, is largely unknown [46]. A possible drug target is the selenoprotein thioredoxin-glutathione reductase (TGR), a key enzyme in the pathway of the parasite for detoxification of reactive oxygen species. The X-ray structure of wild type TGR incubated with AF has shown that the inhibitor is indeed the gold ion released from AF, which is bound at three different sites––mainly Cys residues––not directly involving the C-terminal Sec residue, see Figure 4; further kinetic studies have suggested that the terminal Sec residue mediates the transfer of gold from its ligands in AF to the redox-active Cys pairs of TGR [46].
The usage of in silico molecular docking allowed to infer that AF targets the Sec in the active site of thioredoxin (Trx) glutathione reductase (SmTGR) which shuttles electrons from NADPH to Trx and to the oxidized form of glutathione (GSSG). The tight binding of the [Et3PAu(I)]+ cation to the SmTGR, forming transient interaction with the neighboring protein residues (in particular with Cys596), was found to eventuate with the transfer of gold atoms to SmTGR active site [47].

2.2. Au(I) Complexes with N-heterocyclic Carbene (NHC) Ligands

Metal complexes with N-heterocyclic carbene (NHC) ligands are ubiquitously utilized in chemistry due to their excellent catalytic properties [48]. In the last decade, a lot of studies unveiled the potential of metal NHC compounds in medicinal chemistry since these scaffolds can be employed for the development of powerful and efficacious metallodrugs, active against infectious diseases and cancer [21,49]. Indeed, gold complexes with N-heterocyclic carbene (NHC) ligands were found to constitute a class of promising anticancer metallodrugs with vigorous in vitro and in vivo activities. Nevertheless, their mode of action was not fully elucidated [24,50]. Possible mechanisms for initiating the cell apoptosis process are (i) direct DNA damage, (ii) mitochondrial damage via thioredoxin reductase (TrxR) inhibition and subsequent mitochondrial impairment, (iii) alteration of specific kinases, and (iv) inhibition of proteasome [51]. Nowadays, there is a consensus that the inhibition of TrxR has a crucial part in the therapeutic effect of gold complexes due to the high affinity of gold toward thiol and selenol groups [52]. Moreover, the reduction of Au(III) cations to Au(I) in the cytoplasm augments their preference to favorably bind with the thiol side chains of cysteine residues, including those involved in the coordination of other metals. In this frame, Au(I) breaks the metabolic pathways via the displacement of functional metal atoms, impeding the reduction-oxidation balance and increasing the permeability of cells, and, consequently, inducing the apoptotic response [53].
Computational investigations at the DFT level on the possible usage of Au(I) bis-N-heterocyclic carbene complexes (Figure 5, structure 5) for the attack on Cys and Sec protein residues focused on the initial steps of the carbene ligand substitution by these aminoacidic nucleophiles [54,55,56]. In ref. [55] and [56], the authors employed the capped models of neutral and deprotonated Cys and Sec and explicitly considered the presence of acidic moiety in solution, such as a neighboring acidic residue or the acidic component of the buffer milieu. It was inferred that, although able to afford the concerted protonation of the released carbene, the targeting of neutral cysteine or selenocysteine is kinetically unfavorable in physiological conditions. In contrast, the deprotonated thiolate and selenolate forms of Cys and Sec, respectively, albeit necessitating an external proton source for the protonation of the leaving carbene, tend to be more nucleophilic and, thus, more reactive. The crucial role played by the buffer in the protonation of the leaving carbene spotlights the significance assumed by the biological milieu in the mechanism of action of Au(I) bis-N-heterocyclic carbene complexes.
Neutral N-heterocyclic carbene gold(I) compounds such as IMeAuCl with IMe = (Me2Imy) were recently reported as promising antitumor agents in medicinal chemistry [57]. A detailed density functional theory study was performed by aiming at the thermodynamics and kinetics of the reactions where chloride was substituted with a variety of nucleophilic protein side chains such as arginine, aspartic acid, asparagine, cysteine, glutamic acid, glutamine, histidine, lysine, methionine, selenocysteine, and the N-terminal group (Figure 6, structure 8) [58,59]. It was shown that despite the 23.0 kcal/mol activation-free energy of chloride substitution by water, a rather accessible barrier at physiological conditions, the barriers for the analogous reactions with other nucleophiles are lower (12–20 kcal/mol), thus suggesting that the [Au(I)(NHC)Cl] complex directly attacks its biomolecular targets before the occurrence of hydrolysis. The calculated kinetic reactivity order has been established to be Sec > Lys > activated-Gln/Asn > activated-Cys > Cys > Met > Arg > N-terminal > His > Gln/Asn ≈ Glu/Asp >, well corroborated by the experimentally observed prevalent binding of a [IMeAu]+ fragment to Sec and Cys in thioredoxin. The overall promiscuity of this gold-based scaffold, displayed in the favorable thermodynamic and kinetic values in its reaction with the selected nucleophiles, suggests that its binding preference to a protein side-chain might be rather controlled by the accessibility (bulk exposure) of the targetable sites.
The interaction of a set of gold-based complexes (Figure 7, structures 10a, 10b) with the mammalian and bacterial TrxR was computationally studied by means of docking studies focused on the binding of [(NHC)Au]+ and [(Et3P)Au]+ fragments with both types of TrxR isoforms [60]. Calculations displayed an exceptionally good embedding of the [(NHC)Au]+ fragment of compound 3 in both crystal structures, interacting more actively with the protein residues. Moreover, it was concluded that the selenocysteine in human TrxR favors the breaking up of 3, yielding the [(NHC)Au]+ fragment release.
An accurate computational investigation reveals a detailed kinetic analysis of the reaction of six alkyl-substituted NHC with cysteine (Cys), which is selected for the study as an important bionucleophilic molecule (Figure 5, structures 7a7f) [54]. In this study, the first and second ligand substitution was investigated with the full characterization of their mechanism. For the first reaction step, which turned out to be rate-limiting, the calculated trend of activation enthalpies resulted in being 7a/Me2 < 7f/Me, Et < 7c/n-Bu2 < 7b/i-Pr2 < 7e/Cy2 < 7d/t-Bu2; this order correlates well with the Au-S distances observed in transition states: longer distances relate to higher values of activation enthalpies, consistently with the steric hindrance exerted by the bulky alkyl-substituted NHC ligand on the gold center. On the other hand, the higher activation barrier can be ascribed to electronic modulations originating from the enhanced electron-richness in the substituted NHC ligands, which translates in the strengthening of metal-NHC coordination. Additionally, it was demonstrated that a similar, although somewhat faded, structure––reactivity relationship characterizes the substitution of the second ligand.
Another gold(I)-carbene complex, i.e., [Au(IPr)(Seu)]PF6 with Seu = selenourea and IPr = 1,3-Bis(2,6-diisopropylphenyl)imidazol-2-ylidene, was found to be less potent than cisplatin (cis-diamminedichloroplatinum) in inhibiting the cellular growth in lung carcinoma A549, colon cancer HCT15, and breast cancer MCF7 lines (Figure 6, structure 9) [61], and was also investigated at DFT level of theory. The usage of DFT allowed to optimize its structure and to analyze the charge distribution [61].
The NHC*-Au-thiosugar complexes based on 1,3-dibenzyl-4,5-diphenyl-imidazol-2-ylidene were investigated by experiment and DFT computations and were found to have similar strong and redox-active Au-S bond as in the structurally related auranofin (Figure 7, structures 11a11c) [62]. These complexes are soluble in the biological milieu and show good cytotoxic activity in the medium to low micromolar range, and complex 11a shows higher activity in the low micromolar to nanomolar range against the tested cell lines [62].
Another experimental/theoretical study focused on the NHC*-Gold(I)-X complexes, where NHC* is here represented by 1,3-dibenzyl-4,5-diphenylimidazol-2-ylidene and X = chloride, cyanide, dithiocarbamates, p-mercaptobenzoate or N-acetyl-l-cysteine (Figure 7, structures 12a12e) [63]. The evaluation of the biological activities and the calculated stabilities of these complexes led to the conclusion that NHC*-Au(I)-thiolates 12c, 12d, 12e are biologically more active than the complexes 12a, 12b, which do not possess the Au-S bond.
The joint use of experimental and computational methods allowed the evaluation of a set of gold(I) complexes with a 1,3-diethylbenzimidazol-2-ylidene N-heterocyclic carbene (NHC) ligand of the type NHC-Au-L (L = -Cl, -NHC, or -PPh3) as potential TrxR inhibitors (Figure 8, structures 13a13c) [64]. These metal scaffolds displayed exceptional antiproliferative properties, a strong induction of apoptosis, and enhancement of reactive oxygen species (ROS) formation, thus demonstrating their potential as new antitumor agents. The computational analyses of bond dissociation energies and charge distribution allows an explanation of the stability of these complexes. The employment of DFT calculations disclosed the differences in these three Au(I) complexes: the bond dissociation energy is the lowest in the case of chloride, thus indicating its capability to inhibit TrxR selectively and its reactivity comparable to that of auranofin. On the other hand, it was shown that the most stable among the three studied complexes are less effective in inhibiting TrxR. Indeed, the extensively delocalized positive charge in NHC-Au-NHC and NHC-Au-PPh3, together with the lipophilic structure of the bound ligands, make these Au(I) complexes perfect examples of delocalized lipophilic cations, which can easily permeate the hydrophobic barriers of cellular membranes, in force of the electric gradient between the cytosolic and extracellular compartments. The uptake of these two complexes, therefore, results in three times higher than for the neutral NHC-Au-Cl. Complexes 13b and 13c are well-balanced candidates for the role of antitumor agents.
Gold(I) NHC complexes with phosphines (Figure 8, structures 14a14d) strongly inhibit the proliferation of cancer cells, thanks to an effective intracellular uptake, as shown in a joint experimental-computational study [65]. Their activity toward TrxR and PARP-1 inhibition seems to be strongly related to the size of the alkyl and aryl groups decorating the phosphine ligand. In fact, DFT calculations demonstrated that the dissociation energy of Au-PPh3 bond in [Au(I)(NHC)(PPh3)] is lower than the Au-PR3 bond dissociation energy in [Au(I) (NHC)(PR3)] with R = alkyl, indicating a higher kinetic reactivity of this compound.
A different non-proteic target of Au(I)-NHC complexes are represented by the human telomeric DNA G-quadruplex (Figure 9) [66,67], representing an exciting new possibility in anticancer treatments, and its more complex computational modeling has become an intriguing challenge.
The targeting of the human telomeric (hTelo) and of a promoter sequence (C-KIT1) by Au(I)-NHC complexes was investigated via metadynamics simulations and revealed that these metal fragments might adopt several binding modes as corroborated by FRET DNA melting assays (Figure 10, structures 15a, 15b) [67]. The binding free energies for these Au(I)-DNA adducts were consistent with the corresponding experimental estimates. This study shows that metadynamics is a viable and promising approach for drug design, although its employment in the study of metallodrugs’ pharmacodynamics is still limited.
The biscarbene gold(I) complex [Au(1-butyl-3-methyl-2-ylidene)2]PF6 ([Au(NHC)2]PF6 hereafter) started to receive attention as a perspective anticancer metallodrug (Figure 5, structure 6) [68] due to its remarkable stability under physiological conditions. For instance, it preserves its original structure even in the presence of glutathione (GSH). A combination of high-resolution mass spectrometry and computational DFT and QM/MM methods was used to investigate its interaction [Au(NHC)2]PF6 with the human telomeric DNA G quadruplex, see Figure 9, and it was concluded that [Au(NHC)2]PF6 binding affects the G quadruplex melting temperature and its conformation, although only marginally [66].

2.3. Au(I) Complexes with Other Structures

The growing attention for the entering of thiosemicarbazones and their metal complexes in the field of medicinal chemistry stimulated the development of new gold(I) complexes bearing aryl-thiosemicarbazone moieties (Figure 11, structures 16a, 16b) [26]. These complexes display a satisfactory cytotoxic activity against either tumor (B16-F10 and CT26.WT) or non-tumor cell lines (BHK-21), being even more cytotoxic in tumor cells compared to cisplatin and showing a high selectivity with the phosphine derivatives being more selective than their corresponding non-phosphine gold complexes. Docking investigations gave insights into the binding mode of these gold-based complexes at the TrxR enzyme. Both docking and experimental investigations inferred that either 16a or 16b complexes bearing the p-hydroxyphenyl substituent (Figure 11) turn out to be the most effective inhibitors of TrxR due to their size and peculiar shape that well fit the residues surrounding the anchoring site of TrxR and favor the stabilization of the metal-protein adduct.
Gold(I) complexes with alkynyl and phosphine ligands are perspective scaffolds for metallodrug design and have recently been found to exert an inhibiting activity against TrxR in a joint experimental and computational study [69]. All the six synthesized complexes (Figure 12, structures 17a17f) displayed similar cytotoxicity; however, they were not selective for tumor cells. It was noted that the cytotoxic activity was only marginally affected by the residues at phosphine; a slightly stronger cytotoxic effect was produced when the ethyl moiety (17b) was present. Moreover, it was indicated that the inhibition of TrxR is stronger when alkyl/phenyl residues are available (17a, 17b, 17c), while the complexes incorporating ligands with N/O heteroatoms (17d, 17e, 17f) were less effective.
Novel gold(I) complexes containing tertiary phosphine and the new ligands 5-adamantyl-1,3-thiazolidine-2-thione, 3-methyladamantane–1,3,4-oxadiazole-2-thione (Figure 13, structure 18a18d) were recently synthesized [70]. Coordination of gold to the exocyclic sulfur atom was observed spectroscopically in all four complexes and was corroborated by X-ray results and quantum chemical calculations. All the four compounds were found to be cytotoxic in four different tumor cell lines, colon cancer (CT26.WT), metastatic skin melanoma (B16F10), mammary adenocarcinoma (4T1), and kidney normal cell (BHK-21). Moreover, molecular docking permitted to characterize the interaction between novel antitumor adamantane––azole gold(I) complexes and their potential target, again, thioredoxin reductase [70]. It was demonstrated that the adamantane ring is crucial for the stabilization of the metallodrug-TrxR complex before the formation of a covalent bond between gold and Sec.

3. Computational Studies of Au(III) Metallodrugs

Au(III) complexes have the metal center with d8 configuration, thus they are isostructural and isoelectronic with Pt(II) complexes, assume square planar geometry, and exhibit cytotoxic activity. It was first supposed they could show the same mode of action of platinum complexes, binding to DNA. However, contrarily to platinum(II) compounds, gold (III) analogs are more labile and, in the aqueous and reducing intracellular environment, are particularly suitable to hydrolysis and to reduction to either Au(I) or Au(0). The choice of ligand is therefore particularly important and, indeed, the stability of Au(III) compounds may be substantially enhanced by the appropriate selection of inert ligands, i.e., polydentate ligands with sulfur, oxygen, or nitrogen as donor atoms, including porphyrins, or organometallic cyclometalated scaffolds [71,72,73]. In addition, gold(III) complexes may be more easily modulated to either increase stability or decrease toxicity [74,75,76]. At variance with initial assumptions, later studies have shown that their mode of action nalogously to gold(I) compounds—is based on inhibition of sulfur and selenium-containing enzymes, such as thioredoxin reductase (TrxR), poly (ADP-ribose) polymerase 1 (PARP-1), cathepsins, aquaporins.
A number of gold(III)-based scaffolds were investigated as perspective anticancer drugs both experimentally and computationally, demonstrating effectiveness against cisplatin-resistant tumors and some of them displaying excellent cytotoxicity in vitro and in vivo against solid cancers together with low systemic toxicity [20,74]. Most of these gold(III) antitumor complexes are based on multidentate ligands, including N^N, N^N^N, C^N, C^N^N, C^N^C, porphyrins, and dithiocarbamate, see Figure 14, Figure 15, Figure 16, Figure 17, Figure 18 and Figure 19.
Several computational studies have been carried out in the last years to unveil the interaction between Au(III) center and possible biomolecular targets.

3.1. Interaction with Biologically Relevant Targets

3.1.1. Au(III) Complexes with Chelating N Donor Ligands

The kinetics and the mechanism of the substitution reactions between the mono-functional Au(III) complexes, [Au(dien)Cl]2+ and [Au(terpy)Cl]2+ (dien = 3-azapentane-1,5-diamine, terpy = 2,2′;6′,2″-terpyridine) (Figure 14, structure 19 and 33) and bi-functional Au(III) complexes, [Au(bipy)Cl2]+ and [Au(dach)Cl2]+ (bipy = 2.2′-bipyridine, dach = (1R,2R)-1,2-diaminocyclohexane) (Figure 14, structures 22a and 20) and biologically relevant targets such as L-histidine (L-His), inosine (Ino), inosine-5′-monophosphate (5′-IMP) and guanosine-5′-monophosphate (5′-GMP), were comprehensively investigated [77], revealing an intriguing result in that the reactions of mono-functional complexes were shown to be more rapid than the reactions of bi-functional complexes. DFT calculations corroborated the higher reactivity of [Au(terpy)Cl]2+ compared to [Au(dien)Cl]2+, with activation energy for the chloride exchange in the terpy complex a half lower than the dien complex. On the other hand, the bi-functional [Au(bipy)Cl2]+ complex disclosed an augmented reactivity, compared to [Au(dach)Cl2]+. The reactivity order with the biomolecular targets is similar for all the above Au(III) complexes, hence L-His > 5′-GMP > 5′-IMP > Ino. Quantum chemical calculations [77] indicated that the ligand exchange in most of these complexes takes place via an associative mechanism, whereas the imidazole attack on [Au(terpy)Cl]2+ and [Au(dien)Cl]2+ is better described with an associative interchange mechanism. There is a pronounced correspondence between the reactivity of the studied metal scaffolds toward biological targets and their structural and electronic properties. Moreover, this study elucidated the coordinative binding of gold(III) complexes to 5′-GMP, a DNA nucleophile site of recognized importance since its binding is presumed to account for the anticancer activity.
A recent DFT study focused on the interaction of pyridine gold (III) complexes Au(Hpm)Cl3 and Au(pm)Cl2 (Hpm = 2-pyridylmethanol) with cysteine and purine bases (Figure 14, structure 21a, 21b) [78]. It was demonstrated that hydrogen bonding and proton transfer are indispensable for stabilizing the produced complexes and for diminishing the activation energy. Interestingly, guanine displayed a more marked reactivity with gold (III) complexes in comparison with both adenine and cysteine. Moreover, the biological milieu in which the reaction takes place greatly affects the reaction due to the possibility of the formation of zwitterions. It was shown that the binding of studied gold(III)-based complexes to oxygens of Cys was kinetically more selective than the analogous metalation of cysteine sulfur. In contrast, both guanine and adenine demonstrate complete independence on the substrate types because their reactivity was found to be only marginally affected by the solvent.
Interactions of novel gold(III) complexes with essential biomolecular targets and DNA/BSA were characterized by various experimental methods and by theoretical calculations (Figure 14, structure 24a24c) [79]. Complexes 24a24c showed a good affinity toward different biomolecules, such as Guo, 5′-GMP, and DNA, indicating the capability of these gold(III) complexes to target DNA/BSA. Indeed, the UV-Vis and fluorescence measurements corroborated their interactions with DNA characterized by high binding constants in the range 103–104 (M−1), whereas the calculated constants for the binding of BSA were found to vary between 104 and 105 (M−1). Calculations confirmed the experimental data, and molecular docking studies were used to detail the anchoring of these metal complexes at DNA/BSA. The anticancer potential of the studied Au(III) complexes was supported by cytotoxic activity assays against various cancer cell lines.
Aquaporins (AQPs) are the cellular membrane channels that transfer water, glycerol, and other smaller chemical species, such as hydrogen peroxide, under the influence of osmotic gradients [80]. The AQPs are found in a myriad of human cell types and participate in the control of the urine concentration, metabolism of fat, wound healing, and hydration of skin [81]. Moreover, various isoforms of these transmembrane channels are overexpressed in human tumors, being connected to various types, grades, and stages of cancers [82]. This variability serves as the basis for the identification of AQPs as amenable targets of anticancer drugs. Up to date, no organic inhibitor of AQP was reported as a satisfactory choice for the clinical development of antitumor therapeutic agents. Nevertheless, a set of recent studies on Au(III) metal complexes indicates their unique abilities in providing a selective targeting of the aquaglyceroporins AQP3 [83].
QM/MM and DFT calculations were successfully used for the investigation of the most likely molecular scenario through which the gold(III) complexes with nitrogen donor ligands, such as 1,10-phenatroline, 2,2′-bipyridine, 4,4′-dimethyl-2,2′-bipyridine, 4,4′-diamino-2,2′-bipyridine, and 2,2′;6′,2″-terpyridine (Figure 14, structures 22a22c, Figure 15, structure 29, Figure 16, structure 33), interact with aquaporin AQP3 [84]. This study supports the previously developed model of inhibition of AQP3 by gold(III)-complexes, in which the metal tightly coordinates the thiol side chain of Cys40, which is located near the selectivity filter of this protein channel and induces the channel blockage. It was found that the considered gold(III) complexes manifest selectivity toward AQP3 over AQP1. The identification of the Cys40 of AQP3, proximal to the SF domain, as the most plausible Au(III) coordination site, was initially gained via homology modeling calculations [84] and, subsequently, corroborated via site-directed mutagenesis experiments [85] (vide infra).
The bipyridyl gold(III) complexes Auphen ([Au(phen)Cl2]Cl phen = phenantroline) and Aubipy ([Au(bpy)Cl2]PF6 bpy = bipyridine) have a unique ability to selectively cause the blockage of AQP3. Classical MD simulations produced a comprehensive picture of Aubipy interacting with AQP3 by assessing the structural response of AQP3 to either covalent or non-covalent binding of Aubipy (Figure 14, structure 22a) [86]. The comparison of the unbound, and either covalently—at Cys40 or non-covalently bound Au(III)-AQP3 adducts in terms of their structure and dynamics led to conclude that coordination of gold(III)-based complex induces conformational changes within the AQP3 channel that eventually lead to the impairment of its functionality. Such structural modifications of the AQP3 extracellular pore resulted in a decline of the water molecules flux, as assessed by the analysis of the channel radius and water density along the transmembrane. The most evident structural alteration of the pore ascribed to the binding of the gold fragment is represented by its shrinkage that prevents or limits the fluxing of water along this protein channel. The computational analyses allow gaining insights into the design of AQP3 inhibitors with ameliorated characteristics, in particular, via the insertion of substituents of the aromatic ligands to stabilize and/or reinforce the hydrophobic and π-stacking interactions involved in the Au(III)-AQP3 binding.
Figure 14. Au(III) complexes with azapentanediamine (19), diaminocyclohexane (20), pyridines (21a, 21b), bipyridyl (22a22c), dipyridin-2-ylamine (23), bispyrazolpyridine (24a, 24b) and bisindazolpyridine (24c).
Figure 14. Au(III) complexes with azapentanediamine (19), diaminocyclohexane (20), pyridines (21a, 21b), bipyridyl (22a22c), dipyridin-2-ylamine (23), bispyrazolpyridine (24a, 24b) and bisindazolpyridine (24c).
Molecules 26 07600 g014
Other gold(III)-based complexes (Figure 15, structures 2629, Figure 17, structure 34) [85] were studied by stopped-flow spectroscopy and molecular dynamics combined with DFT. They support the Cys40 as the main target of gold metallodrugs in AQP3 and its blockage induction causing disruption of glycerol and water permeation, as the main modus operandi of the studied metal complexes. Interestingly, the cationic compound 28 was found to be the most effective inhibitor of glycerol permeation, whereas complex 27 with a similar ligand was only marginally active. The computational study allowed an explanation of this aspect with the facile complexation of the cationic complex 28 [Au(PbImMe)Cl2]+ by cysteinato residues, whereas the neutral complexes [Au(pyb)Cl2] 34 and [Au(PbIm)Cl2] 27 are less prone to be complexed due to less favorable electrostatics.
Auphen was found to be an exceptionally effective inhibitor of AQP3 while being a rather weak inhibitor of AQP1, thus displaying a remarkable degree of selectivity (Figure 15, structure 29) [87]. Computational studies joint with experimental insights provided further evidence that its inhibitive properties originate from the Au(III) propensity toward interaction with sulfhydryl moieties of proteins, such as the thiol group of Cys40 of AQP3. Additionally, molecular docking corroborated the selectivity of Auphen to AQP3 with respect to AQP1, which, together with its high solubility in water and a proven inhibitory effect even at low concentrations (in the nanomolar range), makes this metallodrug a perfect candidate for future in vivo studies.
A set of six Au(III) complexes with substituted 1,10-phenanthroline ligands (Figure 15, structures 30a30c, 31, Figure 16, structures 32a, 32b), which inhibit the water and glycerol channel AQP3, were synthesized and studied in a joint experimental-computational study [88]. Both DFT computations and UV-Vis spectrophotometry data showed that these complexes do not undergo hydrolysis via the exchange of chloride ligands with water or hydroxyl anion, which makes them different from Au(III) bipyridyl complexes. Their cytotoxic activity against A549 human lung cancer cells is a sign of their ability to inhibit cell proliferation in vitro, whereas these metal-based scaffolds strongly disrupt the glycerol penetration in human red blood cells (hRBC) by means of AQP3 inhibition. Their affinity for thiols is more pronounced compared to Aubipy ([Au(bipy)Cl2]PF6, bipy = 2,2′-bipyridine), as shown by a QM study of the ligand exchange with methanethiol modeling the Cys40 side chain of AQP3. Moreover, a multilayered computational approach incorporating QM, molecular dynamics MD, hybrid QM/MM methodologies, as well as the analyses of atoms in molecules (AIM) and natural bond orbitals (NBO), was employed for investigating both coordinative and non-covalent binding of these Au(III) complexes to the AQP3 channel, characterizing it as result of complex protein conformational changes taking place concomitantly with the coordinative gold binding, leading to the pore blockage. These kinds of computations are an indispensable basis for the future design of isoform-selective AQP inhibitors since they shed light on the significance of non-coordinative adducts in fine-tuning the AQP3 inhibition.
Figure 15. Au(III) complexes with bipyridyl (25a25c), dipyridin-2-ylamine (26), pyridylbenzimidazole (27, 28), and phenanthroline (29, 30a30c, 31) groups.
Figure 15. Au(III) complexes with bipyridyl (25a25c), dipyridin-2-ylamine (26), pyridylbenzimidazole (27, 28), and phenanthroline (29, 30a30c, 31) groups.
Molecules 26 07600 g015
Figure 16. Au(III) complexes with substituted phenanthroline (32a, 32b) and terpyridine (33).
Figure 16. Au(III) complexes with substituted phenanthroline (32a, 32b) and terpyridine (33).
Molecules 26 07600 g016

3.1.2. Au(III) Complexes with Cyclometalated Ligands

The high redox stability is one of the crucial features of cyclometalated Au(III) compounds characterized by the presence of at least one metal-carbon σ bond [89,90]. However, the decrease in the formal positive charge of in Au(III) compounds gained by cyclometalation substantially alters its biological properties, in particular its reactivity with serum proteins and its uptake in cells [91]. The majority of cyclometalated Au(III) complexes with significant cytotoxicity include a bipyridine or terpyridine ligand and monodentate ligand(s) [92,93]. Cyclometalated Au(III) complexes with tetradentate ligands are instead uncommon in a biological milieu [91,94]. The cyclometalated complexes 35a35c, 36a36f, 37a37c (Figure 17) feature a remarkable redox activity, giving rise to peculiar pathways of protein targeting. Theoretical investigations on the structure and reactivity of such complexes [95,96,97,98] are reported below (vide infra, Section 3.2. Redox stability).
Recently, the anticancer properties of two bis-cyclometalated gold(III) complexes by N6,N6′-di(quinolin-2-yl)-[2,2′-bipyridine]-6,6′-diamine (Figure 17, structure 38) and N6,N6′-di(isoquinolin-3-yl)-[2,2′-bipyridine]-6,6′-diamine (structure 39) were studied in a joint experimental and computational study [99]. Complex 39 was found to be stable in the presence of biological thiols, whereas 38 produced the metastable Au(I) species after reduction in a millimolar concentration of glutathione, thus releasing the TrxR-inhibiting Au+ ions. The redox stability of 39 decreases significantly its ability to inhibit thioredoxin TrxR; however, unlike complex 38, it aggregates into nanoparticles in a biological medium, thus causing a more efficacious gold uptake by cells. Moreover, the lesser promiscuity of 39 makes it ten times more cytotoxic (with respect to cisplatin and 38) to human cancer cells (A549, A431, A375, and MCF7) than to noncancerous cells (MRC5). DFT calculations elucidated the contrasting mechanisms of action of both complexes: 38 was found to produce the TrxR-inhibiting Au+ ions and toxic tetrapyridyl moieties, which were found to target the hERG potassium channel responsible for the generation of cardiac toxicity in vivo, whereas 39 is characterized by a high gold cellular uptake, nuclear DNA damage, and interaction with hERG.
Figure 17. Au(III) complexes with benzylpyridine (34), diphenylpyridine (35a, 35b), terpyridine (35c), diphenylpyridine (36a36f) and naphtalene-2-yl-pyridine (37a37c), tetrapyridyl (38, 39).
Figure 17. Au(III) complexes with benzylpyridine (34), diphenylpyridine (35a, 35b), terpyridine (35c), diphenylpyridine (36a36f) and naphtalene-2-yl-pyridine (37a37c), tetrapyridyl (38, 39).
Molecules 26 07600 g017

3.1.3. Au(III) Complexes with Cyclometalated Ligands

The reactivity and selectivity of gold(III) complexes with dithiocarbamate ligands were investigated computationally (Figure 18, structures 40a, 40b, 41a, 41b) [100], revealing the importance of interatomic sulfur-halogen coordination in the stabilization of the metallic scaffold. Reactivity and selectivity toward TrxR were found to be related to the presence of charge transfer effects––either interatomic, within the gold(III) complex, or intermolecular, between the metal complex and TrxR enzyme. Both in vitro analyses and binding free energy estimates have identified complex 40a as the most reactive towards the TrxR targeting.
Figure 18. Au(III) complexes with dithiocarbamates (40a, 40b, 41a, 41b).
Figure 18. Au(III) complexes with dithiocarbamates (40a, 40b, 41a, 41b).
Molecules 26 07600 g018
Another way to stabilize gold(III) ions in the development of physiologically stable antitumor gold(III) complexes is represented by the tight chelation with porphyrinato ligands [101,102]. A set of gold(III) tetraarylporphyrins with porphyrinato ligands containing different peripheral substituents on the meso-aryl rings were recently synthesized, and their biological activity was investigated both in silico and in vitro (Figure 19, structure 42) [103]. A combination of docking and protein expression analysis experiments indicated the anti-apoptotic protein bcl-2 as the biomolecular target for these gold(III)-porphyrin complexes, the inhibition of which leads to cellular apoptosis. The versatility of the porphyrin structure functionalization together with the highly stable chelation of gold(III) are the features that make this scaffold promising for the further development of physiologically stable anticancer gold(III) complexes.
It was recently shown that the oxo-bridged binuclear gold(III) compounds, [Au2(μ-O)2(N^N)2](PF6)2, where N^N is 2,2′-bipyridine or a substituted 2,2′-bipyridine, possess substantial stability in physiological conditions and reveal cytotoxic activity against several cell lines of human tumors [104]. Selected complexes, such as [Au2(μ-O)2(bipy)2](PF6)2, cis-[Au2(μ-O)2(6-Mebipy)2](PF6)2, trans-[Au2(μ-O)2(6-oXylbipy)2](PF6)2, and [Au2(μ-O)2(6,6′-Me2bipy)2](PF6)2, were crystallographically and computationally characterized. Interestingly, these complexes disclose a correlation between their oxidizing power and their antiproliferative activity (Figure 19, structure 43a43c, 44, 45a, 45b) [105]. Indeed, the most oxidant species [Au2(μ-O)2(6,6′-Me2bipy)2](PF6)2 was found to be also the most reactive towards protein models as well as the most cytotoxic.
Figure 19. Au(III) complexes with porphyrinate (42) and oxo-bridged binuclear gold(III) compounds (43a43c, 44, 45a, 45b).
Figure 19. Au(III) complexes with porphyrinate (42) and oxo-bridged binuclear gold(III) compounds (43a43c, 44, 45a, 45b).
Molecules 26 07600 g019

3.2. Redox Stability

The mechanism of action of most Au(III) complexes is a multistep process, in which the exchange of ligands between the metal scaffold and the cellular nucleophiles takes place after the permeation of the complex in the cytoplasm [106]. Additionally, Au(III) may reduce to Au(I) under physiological conditions, thus opening a route for further ligand substitution reactions of Au(I) compound [73]. In the design of effective Au(III) metallodrugs, the formulation of computational protocols for an accurate prediction of reduction potentials may be particularly insightful, although such estimates, expectedly quite sensitive to the energy contributing terms (mostly nuclear repulsion, electronic energy, thermal correction, and solvation energy) can only be afforded at a high level of theory. The DFT studies reported in the present Section demonstrate the effectiveness of quantum chemical approaches for the determination of redox potential for bioinorganic systems incorporating gold-based scaffolds.
Distinctive features of Au(III) anticancer complexes, which make them unique, are their redox instability and the high affinity of gold towards cellular nucleophilic targets. An exceptionally detailed computational study focused on the interaction of [Au(C^N^C)Cl] probe complex (C^N^C = 2,6-diphenylpyridine) with water and the biomolecular targets represented by simplified molecular models CH3SH/CH3S, CH3Se, and 4-methylimidazole (Figure 17, structure 35a, Scheme 1) [95]. DFT calculations permitted to conclude that the lowest energy reaction path is composed of two consecutive processes: (a) the substitution of chloride by the nucleophile and (b) the reduction of the resulting Au(III) complex to the corresponding Au(I) derivative with the opening of the chelate ring. This study disentangles the reaction mechanism of Au(III) complex and its biomolecular targets and elucidates the processes behind the therapeutic effects.
The biological milieu also has a crucial effect on the reactivity of the AuIII-C^N^C antitumor complexes (Figure 17, structure 35a), as shown in a recent computational study [96]. As mentioned above, the activation of Au(III) complexes can be described as occurring in two steps. In this case, the reduction of Au(III) to Au(I) takes place only after the initial ligand exchange reaction has led the metal fragment to bind at the cellular cysteines in the active site of TrxR. Simplified models were employed in order to imitate the structure of the enzyme characterized by the presence of the C-terminus near the active site of TrxR and in close proximity to the gold fragment bound to Cys. The reduction potential increases with the ligand substitution by Cys in the active site of TrxR, thus, demonstrating that the redox stability of this prototypical Au(III) complex is substantially affected by the nature of the auxiliary ligand, giving a further understanding of how to further fine-tune Au(III) metal scaffolds.
One of the strategies adopted to control the reactivity of Au(III)-based antitumor complexes towards their biomolecular targets is the use of appropriate ligands. Due to the exceptional stability and lipophilicity, tertiary phosphines are widely utilized in the design of Au(III)-based metallodrugs [97]. In order to investigate the effects of various phosphines onto the redox stability of [Au(III)(C^N^C)PR3]+ (Figure 17, structure 35b) (C^N^C = 2,6-diphenylpyridine, PR3 = tertiary phosphine) and [Au(III)(N^N^N)PR3]3+ (N^N^N = 2,2′:6′,2″-terpyridine) complexes (structure 35c), a computational study was carried out to estimate the standard reduction potential εo for Au3+/Au+ [97]. It was found that εo values were negative for structures 35b and positive for structures 35c with varying spreads of 829 and 507 mV, respectively. Interestingly, it was inferred that the buried volume of phosphines, a descriptor proposed in ref. [107], is the major parameter for predicting the stability of complex. The electron-donating ability of phosphines has a greater impact on the redox stability of complexes 35c than complexes 35b since the phosphine affects both Au(III) and Au(I) forms of 35c where the steric effects turn out to be stronger in Au(I) reduced complexes. These considerations are of utmost importance for the prediction of the redox stability of Au(III) complexes, characterizing their reactivity with biomolecular targets.
A source of instability for anticancer Au(III) complexes under physiological conditions originate from their high Au3+/Au+ reduction potential. It is therefore essential to take into account the redox potential in the design of the Au(III) complexes with overall higher stability and improved biological activity. The standard reduction potential of nine Au(III) organometallic complexes of the type [Au3+(R-C^N^C)L]n (Figure 17, structures 36a36f, 37a37c) was computed by the density functional theory approach by producing estimates with the absolute error lower than 170 mV [98]. It was found that the range-corrected functional CAM-B3LYP with the SMD solvation model performed better than B3LYP and ωB97XD and B3LYP, with an absolute error of 82 mV in comparison to the ferrocene/ferrocenium reference electrode potential. Moreover, the employment of variable-temperature H-atom addition/abstraction method for the redox couples with non-equivalent charges resulted in a crucial enhancement of the accuracy in the computation of reduction potential, yielding the mean absolute error of only 87 mV without any scaling procedure, being 1.65-fold lower than 144 mV produced by standard procedure.

4. Challenges, Strengths, and Limitations of Current Computational Approaches

Various computational tools utilized in the studies of the reactivity of gold-based metallodrugs with their biomolecular targets––ab initio quantum chemistry, density functional theory, hybrid quantum mechanics/molecular mechanics (QM/MM), classical molecular dynamics (MD), and molecular docking approaches have their advantages and limitations. Quantum chemistry methods are employed for describing the complex chemistry behind the binding of metal scaffolds with biological targets, whereas the conformational variability and the steric effects are unaccounted for by these methods due to their poor scalability. On the other hand, classical approaches are used for the description of macromolecules and of their large conformational space, yet they lack the accuracy of quantum chemistry. The hybrid QM/MM methods share the strengths and weaknesses of both accurate quantum chemical and scalable classical approaches. The key point is the appropriate selection of the methodology, depending on the size of the studied system and the computational goals.
In most cases contemplated in this review, the computational investigation of the anticancer properties of either Au(I) or Au(III) stems from the identification of the chemical events culminating in the metalation of relevant biological targets. For instance, a myriad of the mononuclear Au(I) complexes [39,40,41,42,43,52,54,55,56,57] were mainly investigated for their capability to exchange labile ligands with suitable protein nucleophilic sites, such as Sec, Cys, His side chains, or other accessible coordinative groups. The support of accurate geometry optimization, typically achieved by means of DFT approach, has allowed in several instances to parallel the experimental data and to further rationalize the trend of activity for Au(I)-based metallodrugs [26,43,45,62,63,64,65,67,69]. To this purpose, the possibility to perform an accurate benchmarking of the level of theory by a thorough examination of basis set/XC-functional combinations is particularly valuable. However, to the best of our knowledge, most computational studies devoted to the structural characterization of Au(I)-metallodrugs reside on the choice of one DFT level of theory without the support of preliminary benchmarking calculations. Among the targets recognized by Au(I)-metallodrugs, the TrxR selenocysteine is undoubtedly identified as the most promising one for its potential role in cancer development [92,93]. Many computational studies [26,39,46,60,70] have therefore been specifically addressed to the metalation of this protein. DFT approaches have been valuably applied to the investigation of the thermodynamics and kinetics of the reaction between Au(I)-metallodrugs and Sec [39,40,41,42,43,52,54,55,56,57]. A challenging aspect that we believe important in the computational study of Au(I)-metallodrugs, and their TrxR targeting, is represented by the proton exchange events concomitant to the metalation process. For example, (NHC)2Au(I) complexes react by exchanging one NHC ligand with the target nucleophilic site; however, the concomitant protonation of the leaving ligand is the thermodynamic driving force. The correct modeling of such a parallel prototropic step is fundamental to correctly assess the thermodynamics of such an NHC/Sec exchange. These considerations have driven the DFT investigation performed by Re et al. [55,56] that hypothesized the bulk assistance to the NHC ligand exchange in the reaction between (NHC)2Au(I) and plausible protein targets, and that led to a significant improvement of the computational prediction of thermodynamics compared to previously reported studies [54]. The correct treatment of the bulk effects is a crucial aspect in the DFT investigation of Au(I)- or Au(III)-metallodrugs. In the latter case, based on the harder character of the Au(III) metal center, the computational study of the possible hydrolytic processes is fundamental to identifying the active species of the metallodrug. For instance, the Aubipy targeting of the Cys40 of AQP3 hypothesized [87] and subsequently confirmed by Casini et al. [84,85] has been computationally analyzed by considering the viability of preliminary aquation steps [86,88]. From a computational point of view, the assessment of the thermodynamics of these bulk processes is particularly challenging because affected by substantial entropy contributions that cannot be easily estimated by routine DFT approaches. In these instances, the use of thermodynamic cycles, allowing the admixture of experimental and theoretical estimates, has allowed establishing that Aubipy and its analogs react with the Cys40 of AQP3 in their chloro-hydroxo form, thus behaving as monofunctional metalation agents [86,88]. Another relevant aspect in the computational investigation of Au(I) and Au(III) metallodrugs is represented by the viability of either Au(III)→Au(I) or Au(I)→Au(0) reduction processes. Indeed, due to the possible reductive elimination affecting either Au(I) or Au(III) scaffolds, many metallodrugs are coordinated by ligands that provide for redox stabilization [95,99,100,101,102,103]. DFT and, more generally, ab initio approaches can be potentially employed to probe the redox stability of gold metallodrugs, even before and after the target metalation has occurred. However, the computational estimate of the reduction potential of these complexes is particularly challenging and requires a careful choice of the exchange-correlation potential and solvation model. On the one hand, a fully theoretical assessment of reduction potentials would generally involve the investigation of reduction half-reactions requiring the explicit handling of elusive species like free electrons. On the other, significant relativistic effects in the core electrons of gold are expected to be a major determinant in the theoretical estimate of Au(III)→Au(I) or Au(I)→Au(0) reduction potentials. In this respect, the routine employment of pseudopotentials for the treatment of relativistic effects may not be adequate in the investigation of the redox stability of gold metallodrugs.
The most common limitation to the tout court application of full DFT approaches in the investigation of gold metallodrugs is ascribed to the effectiveness in the use of reduced models to investigate the metal target. Indeed, when a metallodrug is expected to bind at biological macromolecules, such as proteins or nucleic acids, the structure of the target itself may have a marked influence on the metalation process, for example, by facilitating/limiting the accessibility of the metal scaffold, or by exerting a chemical control on the metalation [9,108,109]. In these cases, QM/MM approaches may be preferred to the full DFT treatment because they allow treating the steric and electrostatic influence of the target structure on the metalation process. For example, the formation of non-covalent adducts between Aubipy-based metallodrugs and AQP3 has been ascertained by QM/MM calculations, which have evidenced the importance of the metal center in the formation of weak interactions with AQP3 residues [86,88].
Classical mechanics-based, MD and docking approaches can also give useful information on the binding site in large biomolecular targets, such as proteins or DNA, provided a careful parametrization of the metal force field is employed. Indeed, the crucial drawback of classical MD simulation of bioinorganic systems is the predominant unavailability of parametrization of the metal ion in the most commonly used empirical force fields. For this purpose, the use of MD and docking approaches to investigate the interaction between gold complexes and biomolecular target requires precursory quantum chemical studies to parametrize the gold scaffold and obtain appropriate parameters for their successive employment in the empirical force field. For instance, as discussed above, MD simulations allowed a comprehensive picture of Aubipy interacting with AQP3 by assessing the structural response of AQP3 to either covalent or non-covalent binding of Aubipy [86]. In a recent study, stopped-flow spectroscopy and molecular dynamics combined with DFT allowed to identify Cys40 as the main target of gold metallodrugs 2629 and 34 in AQP3 and the induced mechanism of disruption of glycerol and water permeation, shedding light on the modus operandi of these metal complexes [85].

5. Conclusions

The success of cisplatin as an anticancer drug has led to numerous studies focused on the design of novel transition metal compounds. Gold-based complexes are one of the main directions in this research due to the strong affinity of Au(I) to cysteine and selenocysteine side chains of the protein residues. This targeting allows the inhibition of sulfur- and selenium-containing enzymes such as glutathione reductases, glutathione peroxidase, glutathione-S-transferase, cysteine protease, thioredoxin reductase (TrxR), and poly (ADP-ribose) polymerase 1 (PARP-1).
For the delivery of Au(I) into the tumor cells, it is necessary to decorate it with ligands that fine-tune its reactivity and convey these gold compounds to their targets, preventing them from participating in chemical side reactions beforehand. One of the possible routes is the use of specific ligands at the linear di-coordinate Au(I) metal center, such as phosphine, thiosugar, N-heterocyclic carbene, alkynyl, and thiosemicarbazone. These drawbacks are even more pronounced for Au(III) complexes which are labile and prone to reduction to either Au(I) or Au(0) in the aqueous or physiological milieu, thus requiring an appropriate selection of inert ligands which produce the stabilizing effect on the metal core. Usually, polydentate ligands with sulfur, oxygen, or nitrogen as donor atoms, including porphyrins or organometallic cyclometalated scaffolds, are employed. A balanced choice allows for an accurate adjustment of stability and toxicity, making gold(III)-based scaffolds perspective anticancer drugs.
The role of computational chemistry in the quest for new metallodrugs has grown due to the increasing efficacy of the methods and to the parallel increase in the hardware capabilities. Ab initio and DFT methods are nowadays the best choices for describing the mode of action of the metal-based complexes characterized by the promiscuity of metal center, its multitargeting modus operandi, the widespread necessity for activation via hydrolysis or solvolysis, and the substitution character of the metallodrug reactions with their targets.
This review is focused on the recent computational studies of Au(I) and Au(III) antitumor compounds, this topic being of utmost importance for the design of novel gold-based metallodrugs. Considering the described successes of computational modeling of gold-based complexes, we highlight the importance of in silico studies of developed metal compounds as a crucial phase even before in vivo and in vitro testing. The comprehension of the mechanistic aspects underlying the metallodrug activation, target selection, target binding is important for the design of metal-based compounds and their further structural improvement, as well as for the reduction in their toxicity effects and the enhancing of targeted drug delivery, allowing a rational drug design.

Author Contributions

Conceptualization, C.C. and N.R.; formal analysis, I.T.; writing—original draft preparation, I.T.; writing—review and editing, A.M, C.C. and N.R.; visualization, I.T. and A.M.; supervision, N.R; funding, CC. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

I.T. gratefully acknowledges the hosting of Institut de Chimie Moleculaire de l’Université de Bourgogne (ICMUB).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pricker, S.P. Medical uses of gold compounds: Past, present and future. Gold Bull. 1996, 29, 53–60. [Google Scholar] [CrossRef] [Green Version]
  2. Messori, L.; Marcon, G. Gold complexes in the treatment of rheumatoid arthritis. Met. Ions Biol. Syst. 2004, 41, 279–304. [Google Scholar]
  3. Ott, I. On the medicinal chemistry of gold complexes as anticancer drugs. Coord. Chem. Rev. 2009, 253, 1670–1681. [Google Scholar] [CrossRef]
  4. Zou, T.; Lum, C.T.; Lok, C.N.; Zhang, J.J.; Che, C.M. Chemical biology of anticancer gold (III) and gold (I) complexes. Chem. Soc. Rev. 2015, 44, 8786–8801. [Google Scholar] [CrossRef] [PubMed]
  5. Verma, J.; Khedkar, V.M.; Coutinho, E.C. 3D-QSAR in drug design-a review. Curr. Top. Med. Chem. 2010, 10, 95–115. [Google Scholar] [CrossRef] [PubMed]
  6. Palermo, G.; Spinello, A.; Saha, A.; Magistrato, A. Frontiers of metal-coordinating drug design. Expert Opin. Drug Dis. 2021, 16, 497–511. [Google Scholar] [CrossRef]
  7. Lee, R.F.; Menin, L.; Patiny, L.; Ortiz, D.; Dyson, P.J. Versatile tool for the analysis of metal–protein interactions reveals the promiscuity of metallodrug–protein interactions. Anal. Chem. 2017, 89, 11985–11989. [Google Scholar] [CrossRef]
  8. Sullivan, M.P.; Cziferszky, M.; Tolbatov, I.; Truong, D.; Mercadante, D.; Re, N.; Gust, R.; Goldstone, D.C.; Hartinger, C.G. Probing the paradigm of promiscuity for N-heterocyclic carbene complexes and their protein adduct formation. Angew. Chem. Int. Ed. 2021, 60, 19928. [Google Scholar] [CrossRef] [PubMed]
  9. Tolbatov, I.; Marrone, A.; Paciotti, R.; Re, N.; Coletti, C. Multilayered Modelling of the Metallation of Biological Targets. In International Conference on Computational Science and Its Applications; Springer: Cham, Swizerland, 2021; pp. 398–412. [Google Scholar]
  10. Paciotti, R.; Tolbatov, I.; Graziani, V.; Marrone, A.; Re, N.; Coletti, C. Insights on the activity of platinum-based anticancer complexes through computational methods. In AIP Conference Proceedings; AIP Publishing LLC: Melville, NY, USA, 2018; Volume 2040, No. 1; p. 020019. [Google Scholar]
  11. Barresi, E.; Tolbatov, I.; Marzo, T.; Zappelli, E.; Marrone, A.; Re, N.; Pratesi, A.; Martini, C.; Taliani, S.; Da Settimo, F.; et al. Two mixed valence diruthenium (II, III) isomeric complexes show different anticancer properties. Dalton Trans. 2021, 50, 9643–9647. [Google Scholar] [CrossRef] [PubMed]
  12. Tolbatov, I.; Coletti, C.; Marrone, A.; Re, N. Reactivity of arsenoplatin complex versus water and thiocyanate: A DFT benchmark study. Theor. Chem. Acc. 2020, 139, 1–11. [Google Scholar] [CrossRef]
  13. Barresi, E.; Tolbatov, I.; Pratesi, A.; Notarstefano, V.; Baglini, E.; Daniele, S.; Taliani, S.; Re, N.; Giorgini, E.; Martini, C.; et al. A mixed-valence diruthenium (II, III) complex endowed with high stability: From experimental evidence to theoretical interpretation. Dalton Trans. 2020, 49, 14520–14527. [Google Scholar] [CrossRef] [PubMed]
  14. Tolbatov, I.; Marzo, T.; Cirri, D.; Gabbiani, C.; Coletti, C.; Marrone, A.; Paciotti, R.; Messori, L.; Re, N. Reactions of cisplatin and cis-[PtI2(NH3)2] with molecular models of relevant protein side-chains: A comparative analysis. J. Inorg. Biochem. 2020, 209, 111096. [Google Scholar] [CrossRef] [PubMed]
  15. 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]
  16. 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]
  17. 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. Comm. 2013, 49, 10100–10102. [Google Scholar] [CrossRef] [PubMed]
  18. 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] [Green Version]
  19. Ferraro, G.; Gabbiani, C.; Merlino, A. First crystal structure for a gold carbene–protein adduct. Bioconjug. Chem. 2016, 27, 1584–1587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. 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] [PubMed]
  21. Oehninger, L.; Rubbiani, R.; Ott, I. N-Heterocyclic carbene metal complexes in medicinal chemistry. Dalton Trans. 2013, 42, 3269–3284. [Google Scholar] [CrossRef]
  22. Carlos Lima, J.; Rodriguez, L. Phosphine-gold (I) compounds as anticancer agents: General description and mechanisms of action. Anticancer Agents Med. Chem. 2011, 11, 921–928. [Google Scholar] [CrossRef] [PubMed]
  23. Onodera, T.; Momose, I.; Kawada, M. Potential anticancer activity of auranofin. Chem. Pharmaceut. Bull. 2019, 67, 186–191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhang, C.; Maddelein, M.L.; Sun, R.W.Y.; Gornitzka, H.; Cuvillier, O.; Hemmert, C. Pharmacomodulation on gold-NHC complexes for anticancer applications–is lipophilicity the key point? Eur. J. Med. Chem. 2018, 157, 320–332. [Google Scholar] [CrossRef] [PubMed]
  25. Yang, Z.; Jiang, G.; Xu, Z.; Zhao, S.; Liu, W. Advances in alkynyl gold complexes for use as potential anticancer agents. Coord. Chem. Rev. 2020, 423, 213492. [Google Scholar] [CrossRef]
  26. Tavares, T.T.; Azevedo, G.C.; Garcia, A.; Carpanez, A.G.; Lewer, P.M.; Paschoal, D.; Müller, B.L.; Dos Santos, H.F.; Matos, R.C.; Silva, H.; et al. Gold (I) complexes with aryl-thiosemicarbazones: Molecular modeling, synthesis, cytotoxicity and TrxR inhibition. Polyhedron 2017, 132, 95–104. [Google Scholar] [CrossRef]
  27. AbdelKhalek, A.; Abutaleb, N.S.; Mohammad, H.; Seleem, M.N. Antibacterial and antivirulence activities of auranofin against Clostridium difficile. Int. J. Antimicrob. 2019, 53, 54–62. [Google Scholar] [CrossRef]
  28. Mbouaka, A.L.; Leitsch, D.; Koehsler, M.; Walochnik, J. Antimicrobial effect of auranofin against Acanthamoeba spp. Int. J. Antimicrob. 2021, 58, 106425. [Google Scholar] [CrossRef] [PubMed]
  29. Wiederhold, N.P.; Patterson, T.F.; Srinivasan, A.; Chaturvedi, A.K.; Fothergill, A.W.; Wormley, F.L.; Ramasubramanian, A.K.; Lopez-Ribot, J.L. Repurposing auranofin as an antifungal: In vitro activity against a variety of medically important fungi. Virulence 2017, 8, 138–142. [Google Scholar] [CrossRef] [Green Version]
  30. Li, H.; Hu, J.; Wu, S.; Wang, L.; Cao, X.; Zhang, X.; Dai, B.; Cao, M.; Shao, R.; Zhang, R.; et al. Auranofin-mediated inhibition of PI3K/AKT/mTOR axis and anticancer activity in non-small cell lung cancer cells. Oncotarget 2016, 7, 3548. [Google Scholar] [CrossRef] [Green Version]
  31. Yue, S.; Luo, M.; Liu, H.; Wei, S. Recent advances of gold compounds in anticancer immunity. Front. Chem. 2020, 8, 543. [Google Scholar] [CrossRef]
  32. Pratesi, A.; Cirri, D.; Ciofi, L.; Messori, L. Reactions of auranofin and its pseudohalide derivatives with serum albumin investigated through ESI-Q-TOF MS. Inorg. Chem. 2018, 57, 10507–10510. [Google Scholar] [CrossRef]
  33. Liu, N.; Guo, Z.; Xia, X.; Liao, Y.; Zhang, F.; Huang, C.; Liu, Y.; Deng, X.; Jiang, L.; Wang, X.; et al. Auranofin lethality to prostate cancer includes inhibition of proteasomal deubiquitinases and disrupted androgen receptor signaling. Eur. J. Pharmacol. 2019, 846, 1–11. [Google Scholar] [CrossRef] [PubMed]
  34. Hwangbo, H.; Ji, S.Y.; Kim, M.Y.; Kim, S.Y.; Lee, H.; Kim, G.Y.; Kim, S.; Cheong, J.; Choi, Y.H. Anti-inflammatory effect of auranofin on palmitic acid and LPS-induced inflammatory response by modulating TLR4 and NOX4-mediated NF-κB signaling pathway in RAW264.7 Macrophages. Int. J. Mol. Sci. 2021, 22, 5920. [Google Scholar] [CrossRef]
  35. Zhang, X.; Selvaraju, K.; Saei, A.A.; D’Arcy, P.; Zubarev, R.A.; Arnér, E.S.; Linder, S. Repurposing of auranofin: Thioredoxin reductase remains a primary target of the drug. Biochimie 2019, 162, 46–54. [Google Scholar] [CrossRef]
  36. Landini, I.; Massai, L.; Cirri, D.; Gamberi, T.; Paoli, P.; Messori, L.; Mini, E.; Nobili, S. Structure-activity relationships in a series of auranofin analogues showing remarkable antiproliferative properties. J. Inorg. Biochem. 2020, 208, 111079. [Google Scholar] [CrossRef]
  37. Marzo, T.; Cirri, D.; Pollini, S.; Prato, M.; Fallani, S.; Cassetta, M.I.; Novelli, A.; Rossolini, G.M.; Messori, L. Auranofin and its analogues show potent antimicrobial activity against multidrug-resistant pathogens: Structure-activity relationships. ChemMedChem 2018, 13, 2448–2454. [Google Scholar] [CrossRef] [PubMed]
  38. Marzo, T.; Cirri, D.; Gabbiani, C.; Gamberi, T.; Magherini, F.; Pratesi, A.; Guerri, A.; Biver, T.; Binacchi, F.; Stefanini, M.; et al. Auranofin, Et3PauCl, and Et3PAuI are highly cytotoxic on colorectal cancer cells: A chemical and biological study. ACS Med. Chem. Lett. 2017, 8, 997–1001. [Google Scholar] [CrossRef]
  39. Howell, J.A. DFT investigation of the interaction between gold (I) complexes and the active site of thioredoxin reductase. J. Organomet. Chem. 2009, 694, 868–873. [Google Scholar] [CrossRef]
  40. Dos Santos, H.F. Reactivity of auranofin with S-, Se-and N-containing amino acids. Comput. Theor. Chem. 2014, 1048, 95–101. [Google Scholar] [CrossRef]
  41. Di Sarra, F.; Fresch, B.; Bini, R.; Saielli, G.; Bagno, A. Reactivity of auranofin with selenols and thiols–implications for the anticancer activity of gold (I) compounds. Eur. J. Inorg. Chem. 2013, 2013, 2718–2727. [Google Scholar] [CrossRef]
  42. Shoeib, T.; Atkinson, D.W.; Sharp, B.L. Structural analysis of the anti-arthritic drug Auranofin: Its complexes with cysteine, selenocysteine and their fragmentation products. Inorg. Chim. Acta 2010, 363, 184–192. [Google Scholar] [CrossRef]
  43. Tolbatov, I.; Cirri, D.; Marchetti, L.; Marrone, A.; Coletti, C.; Re, N.; La Mendola, D.; Messori, L.; Marzo, T.; Gabbiani, C.; et al. Mechanistic insights into the anticancer properties of the auranofin analog Au(PEt3)I: A theoretical and experimental study. Front. Chem. 2020, 8, 812. [Google Scholar] [CrossRef]
  44. Madeira, J.M.; Gibson, D.L.; Kean, W.F.; Klegeris, A. The biological activity of auranofin: Implications for novel treatment of diseases. Inflammopharmacology 2012, 20, 297–306. [Google Scholar] [CrossRef]
  45. Gunatilleke, S.S.; De Oliveira, C.A.F.; McCammon, J.A.; Barrios, A.M. Inhibition of cathepsin B by Au (I) complexes: A kinetic and computational study. J. Biol. Inorg. Chem. 2008, 13, 555–561. [Google Scholar] [CrossRef] [PubMed]
  46. Angelucci, F.; Sayed, A.A.; Williams, D.L.; Boumis, G.; Brunori, M.; Dimastrogiovanni, D.; Miele, A.E.; Pauly, F.; Bellelli, A. Inhibition of Schistosoma mansoni thioredoxin-glutathione reductase by auranofin: Structural and kinetic aspects. J. Biol. Chem. 2009, 284, 28977–28985. [Google Scholar] [CrossRef] [Green Version]
  47. Caroli, A.; Simeoni, S.; Lepore, R.; Tramontano, A.; Via, A. Investigation of a potential mechanism for the inhibition of SmTGR by Auranofin and its implications for Plasmodium falciparum inhibition. Biochem. Biophys. Res. Comm. 2012, 417, 576–581. [Google Scholar] [CrossRef] [PubMed]
  48. Velazquez, H.D.; Verpoort, F. N-heterocyclic carbene transition metal complexes for catalysis in aqueous media. Chem. Soc. Rev. 2012, 41, 7032–7060. [Google Scholar] [CrossRef] [PubMed]
  49. Gautier, A.; Cisnetti, F. Advances in metal–carbene complexes as potent anti-cancer agents. Metallomics 2012, 4, 23–32. [Google Scholar] [CrossRef]
  50. Porchia, M.; Pellei, M.; Marinelli, M.; Tisato, F.; Del Bello, F.; Santini, C. New insights in Au-NHCs complexes as anticancer agents. Eur. J. Med. Chem. 2018, 146, 709–746. [Google Scholar] [CrossRef] [PubMed]
  51. Mora, M.; Gimeno, M.C.; Visbal, R. Recent advances in gold–NHC complexes with biological properties. Chem. Soc. Rev. 2019, 48, 447–462. [Google Scholar] [CrossRef]
  52. Rubbiani, R.; Schuh, E.; Meyer, A.; Lemke, J.; Wimberg, J.; Metzler-Nolte, N.; Meyer, F.; Mohr, F.; Ott, I. TrxR inhibition and antiproliferative activities of structurally diverse gold N-heterocyclic carbene complexes. MedChemComm 2013, 4, 942–948. [Google Scholar] [CrossRef]
  53. Magherini, F.; Fiaschi, T.; Valocchia, E.; Becatti, M.; Pratesi, A.; Marzo, T.; Massai, L.; Gabbiani, C.; Landini, I.; Nobili, S.; et al. Antiproliferative effects of two gold(I)-N-heterocyclic carbene complexes in A2780 human ovarian cancer cells: A comparative proteomic study. Oncotarget 2018, 9, 28042–28068. [Google Scholar] [CrossRef] [PubMed]
  54. Dos Santos, H.F.; Vieira, M.A.; Sánchez Delgado, G.Y.; Paschoal, D. Ligand exchange reaction of Au (I) RN-heterocyclic carbene complexes with cysteine. J. Phys. Chem. A 2016, 120, 2250–2259. [Google Scholar] [CrossRef]
  55. Tolbatov, I.; Coletti, C.; Marrone, A.; Re, N. Insight into the substitution mechanism of antitumor Au (I) N-heterocyclic carbene complexes by cysteine and selenocysteine. Inorg. Chem. 2020, 59, 3312–3320. [Google Scholar] [CrossRef] [PubMed]
  56. Tolbatov, I.; Marzo, T.; Coletti, C.; La Mendola, D.; Storchi, L.; Re, N.; Marrone, A. Reactivity of antitumor coinage metal-based N-heterocyclic carbene complexes with cysteine and selenocysteine protein sites. J. Inorg. Biochem. 2021, 223, 111533. [Google Scholar] [CrossRef]
  57. Liu, W.; Gust, R. Metal N-heterocyclic carbene complexes as potential antitumor metallodrugs. Chem. Soc. Rev. 2013, 42, 755–773. [Google Scholar] [CrossRef]
  58. 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. [Google Scholar] [CrossRef] [Green Version]
  59. Paciotti, R.; Tolbatov, I.; Marrone, A.; Storchi, L.; Re, N.; Coletti, C. Computational investigations of bioinorganic complexes: The case of calcium, gold and platinum ions. In AIP Conference Proceedings; AIP Publishing LLC: Melville, NY, USA, 2019; Volume 2186, No. 1; p. 030011. [Google Scholar]
  60. Schmidt, C.; Albrecht, L.; Balasupramaniam, S.; Misgeld, R.; Karge, B.; Brönstrup, M.; Prokop, A.; Baumann, K.; Reichl, S.; Ott, I. A gold (I) biscarbene complex with improved activity as a TrxR inhibitor and cytotoxic drug: Comparative studies with different gold metallodrugs. Metallomics 2019, 11, 533–545. [Google Scholar] [CrossRef]
  61. Seliman, A.A.; Altaf, M.; Odewunmi, N.A.; Kawde, A.N.; Zierkiewicz, W.; Ahmad, S.; Altuwaijri, S.; Isab, A.A. Synthesis, X-ray structure, DFT calculations and anticancer activity of a selenourea coordinated gold (I)-carbene complex. Polyhedron 2017, 137, 197–206. [Google Scholar] [CrossRef]
  62. Dada, O.; Sánchez-Sanz, G.; Tacke, M.; Zhu, X. Synthesis and anticancer activity of novel NHC-gold (I)-sugar complexes. Tetrahedron Lett. 2018, 59, 2904–2908. [Google Scholar] [CrossRef]
  63. Curran, D.; Dada, O.; Müller-Bunz, H.; Rothemund, M.; Sánchez-Sanz, G.; Schobert, R.; Zhu, X.; Tacke, M. Synthesis and cytotoxicity studies of novel NHC*-gold (I) complexes derived from lepidiline A. Molecules 2018, 23, 2031. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Rubbiani, R.; Can, S.; Kitanovic, I.; Alborzinia, H.; Stefanopoulou, M.; Kokoschka, M.; Mönchgesang, S.; Sheldrick, W.S.; Wölfl, S.; Ott, I. Comparative in vitro evaluation of N-heterocyclic carbene gold (I) complexes of the benzimidazolylidene type. J. Med. Chem. 2011, 54, 8646–8657. [Google Scholar] [CrossRef] [PubMed]
  65. Rubbiani, R.; Salassa, L.; De Almeida, A.; Casini, A.; Ott, I. Cytotoxic gold (I) N-heterocyclic carbene complexes with phosphane ligands as potent enzyme inhibitors. ChemMedChem 2014, 9, 1205–1210. [Google Scholar] [CrossRef] [PubMed]
  66. Guarra, F.; Marzo, T.; Ferraroni, M.; Papi, F.; Bazzicalupi, C.; Gratteri, P.; Pescitelli, G.; Messori, L.; Biver, T.; Gabbiani, C. Interaction of a gold (I) dicarbene anticancer drug with human telomeric DNA G-quadruplex: Solution and computationally aided X-ray diffraction analysis. Dalton Trans. 2018, 47, 16132–16138. [Google Scholar] [CrossRef] [PubMed]
  67. Wragg, D.; De Almeida, A.; Bonsignore, R.; Kühn, F.E.; Leoni, S.; Casini, A. On the mechanism of Gold/NHC compounds binding to DNA G-quadruplexes: Combined metadynamics and biophysical methods. Angew. Chem. 2018, 130, 14732–14736. [Google Scholar] [CrossRef]
  68. Maji, B.; Bhattacharya, S. Advances in the molecular design of potential anticancer agents via targeting of human telomeric DNA. ChemComm 2014, 50, 6422–6438. [Google Scholar]
  69. Andermark, V.; Göke, K.; Kokoschka, M.; El Maaty, M.A.A.; Lum, C.T.; Zou, T.; Sun, R.W.-Y.; Aguiló, E.; Oehninger, L.; Rodríguez, L.; et al. Alkynyl gold (I) phosphane complexes: Evaluation of structure–activity-relationships for the phosphane ligands, effects on key signaling proteins and preliminary in-vivo studies with a nanoformulated complex. J. Inorg. Biochem. 2016, 160, 140–148. [Google Scholar] [CrossRef] [Green Version]
  70. Garcia, A.; Machado, R.C.; Grazul, R.M.; Lopes, M.T.P.; Corrêa, C.C.; Dos Santos, H.F.; de Almeida, M.V.; Silva, H. Novel antitumor adamantane–azole gold (I) complexes as potential inhibitors of thioredoxin reductase. J. Biol. Inorg. Chem. 2016, 21, 275–292. [Google Scholar] [CrossRef]
  71. 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]
  72. Gabbiani, C.; Casini, A.; Messori, L. Gold(III) compounds as anticancer drugs. Gold Bull. 2007, 40, 73–81. [Google Scholar] [CrossRef] [Green Version]
  73. Bertrand, B.; Williams, M.R.; Bochmann, M. Gold (III) complexes for antitumor applications: An overview. Chem. Eur. J. 2018, 24, 11840–11851. [Google Scholar] [CrossRef]
  74. Casini, A.; Hartinger, C.; Gabbiani, C.; Mini, E.; Dyson, P.J.; Keppler, B.K.; Messori, L. Gold(III) compounds as anticancer agents: Relevance of gold–protein interactions for their mechanism of action. J. Inorg. Biochem. 2008, 102, 564–575. [Google Scholar] [CrossRef]
  75. Tiekink, E.R. Gold derivatives for the treatment of cancer. Crit. Rev. Oncol. 2002, 42, 225–248. [Google Scholar] [CrossRef]
  76. 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]
  77. Djeković, A.; Petrović, B.; Bugarčić, Ž.D.; Puchta, R.; van Eldik, R. Kinetics and mechanism of the reactions of Au(III) complexes with some biologically relevant molecules. Dalton Trans. 2012, 41, 3633–3641. [Google Scholar] [CrossRef]
  78. He, Y.; Zhou, L. A theoretical study on pyridine gold (III) complexes AuCl3 (Hpm) and AuCl2 (pm) targeting purine bases and cysteine. Comput. Theor. Chem. 2016, 1093, 20–28. [Google Scholar] [CrossRef]
  79. Radisavljević, S.; Petrović, B. Gold(III) complexes: An overview on their kinetics, interactions with DNA/BSA, cytotoxic activity, and computational calculations. Front. Chem. 2020, 8, 379. [Google Scholar] [CrossRef] [PubMed]
  80. Verkman, A.S.; Anderson, M.O.; Papadopoulos, M.C. Aquaporins: Important but elusive drug targets. Nat. Rev. Drug Discov. 2014, 13, 259–277. [Google Scholar] [CrossRef] [Green Version]
  81. Verkman, A.S. Aquaporins in clinical medicine. Ann. Rev. Med. 2012, 63, 303–316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Ribatti, D.; Ranieri, G.; Annese, T.; Nico, B. Aquaporins in cancer. Biochim. Biophys. Acta (BBA)-Gen. Subjects 2014, 1840, 1550–1553. [Google Scholar] [CrossRef] [PubMed]
  83. Aikman, B.; De Almeida, A.; Meier-Menches, S.M.; Casini, A. Aquaporins in cancer development: Opportunities for bioinorganic chemistry to contribute novel chemical probes and therapeutic agents. Metallomics 2018, 10, 696–712. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Martins, A.P.; Ciancetta, A.; De Almeida, A.; Marrone, A.; Re, N.; Soveral, G.; Casini, A. Aquaporin inhibition by gold (III) compounds: New insights. ChemMedChem 2013, 8, 1086–1092. [Google Scholar] [CrossRef]
  85. De Almeida, A.; Mósca, A.F.; Wragg, D.; Wenzel, M.; Kavanagh, P.; Barone, G.; Leoni, S.; Soveral, G.; Casini, A. The mechanism of aquaporin inhibition by gold compounds elucidated by biophysical and computational methods. ChemComm 2017, 53, 3830–3833. [Google Scholar] [CrossRef]
  86. Graziani, V.; Marrone, A.; Re, N.; Coletti, C.; Platts, J.A.; Casini, A. A multilevel 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] [PubMed]
  87. Martins, A.P.; Marrone, A.; Ciancetta, A.; Galán Cobo, A.; Echevarría, M.; Moura, T.F.; Re, N.; Casini, A.; Soveral, G. Targeting aquaporin function: Potent inhibition of aquaglyceroporin-3 by a gold-based compound. PLoS ONE 2012, 7, e37435. [Google Scholar]
  88. Wenzel, M.N.; Mósca, A.F.; Graziani, V.; Aikman, B.; Thomas, S.R.; De Almeida, A.; Platts, J.A.; 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] [PubMed]
  89. Cutillas, N.; Yellol, G.S.; de Haro, C.; Vicente, C.; Rodríguez, V.; Ruiz, J. Anticancer cyclometalated complexes of platinum group metals and gold. Coord. Chem. Rev. 2013, 257, 2784–2797. [Google Scholar] [CrossRef]
  90. Li, C.K.L.; Sun, R.W.Y.; Kui, S.C.F.; Zhu, N.; Che, C.M. Anticancer cyclometalated [AuIIIm (C^N^C)mL]n+ compounds: Synthesis and cytotoxic properties. Chem. Eur. J. 2006, 12, 5253–5266. [Google Scholar] [CrossRef]
  91. Kumar, R.; Nevado, C. Cyclometalated gold (III) complexes: Synthesis, reactivity, and physicochemical properties. Angew. Chem. Int. Ed. 2017, 56, 1994–2015. [Google Scholar] [CrossRef]
  92. 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]
  93. Jia, J.J.; Geng, W.S.; Wang, Z.Q.; Chen, L.; Zeng, X.S. The role of thioredoxin system in cancer: Strategy for cancer therapy. Cancer Chemother. Pharmacol. 2019, 84, 453–470. [Google Scholar] [CrossRef]
  94. Tang, M.C.; Li, L.K.; Lai, S.L.; Cheung, W.L.; Ng, M.; Wong, C.Y.; Chan, M.-Y.; Yam, V.W.W. Design strategy towards horizontally oriented luminescent tetradentate-ligand-containing gold(III) systems. Angew. Chem. Int. Ed. 2020, 59, 21023–21031. [Google Scholar] [CrossRef]
  95. Delgado, G.S.; Paschoal, D.; Dos Santos, H.F. Reactivity of the [Au(C^N^C)Cl] complex in the presence of H2O and N-, S- and Se-containing nucleophiles: A DFT study. J. Biol. Inorg. Chem. 2018, 23, 1283–1293. [Google Scholar] [CrossRef] [PubMed]
  96. Delgado, G.Y.S.; Arvellos, J.F.; Paschoal, D.F.; Dos Santos, H.F. Role of the enzymatic environment in the reactivity of the AuIII-C^N^C anticancer complexes. Inorg. Chem. 2021, 60, 3181–3195. [Google Scholar] [CrossRef]
  97. Delgado, G.Y.S.; Paschoal, D.; de Oliveira, M.A.; Dos Santos, H.F. Structure and redox stability of [Au(III)(X^N^X)PR3] complexes (X = C or N) in aqueous solution: The role of phosphine auxiliary ligand. J. Inorg. Biochem. 2019, 200, 110804. [Google Scholar] [CrossRef]
  98. Delgado, G.S.; Paschoal, D.; Dos Santos, H.F. Predicting standard reduction potential for anticancer Au (III)-complexes: A DFT study. Comput. Theor. Chem. 2017, 1108, 86–92. [Google Scholar] [CrossRef]
  99. Zhou, X.Q.; Carbo-Bague, I.; Siegler, M.A.; Hilgendorf, J.; Basu, U.; Ott, I.; Liu, R.; Zhang, L.; Ramu, V.; IJzerman, A.P.; et al. Rollover cyclometalation vs nitrogen coordination in tetrapyridyl anticancer gold (III) complexes: Effect on protein interaction and toxicity. JACS Au 2021, 1, 380–395. [Google Scholar] [CrossRef]
  100. Lawal, M.M.; Lawal, I.A.; Klink, M.J.; Tolufashe, G.F.; Ndagi, U.; Kumalo, H.M. Density functional theory study of gold (III)-dithiocarbamate complexes with characteristic anticancer potentials. J. Inorg. Biochem. 2020, 206, 111044. [Google Scholar] [CrossRef] [PubMed]
  101. Che, C.M.; Sun, R.W.Y.; Yu, W.Y.; Ko, C.B.; Zhu, N.; Sun, H. Gold (III) porphyrins as a new class of anticancer drugs: Cytotoxicity, DNA binding and induction of apoptosis in human cervix epitheloid cancer cells. Chem. Comm. 2003, 1718–1719. [Google Scholar] [CrossRef]
  102. Sun, R.W.Y.; Che, C.M. The anti-cancer properties of gold (III) compounds with dianionic porphyrin and tetradentate ligands. Coord. Chem. Rev. 2009, 253, 1682–1691. [Google Scholar] [CrossRef]
  103. 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] [PubMed]
  104. Casini, A.; Cinellu, M.A.; Minghetti, G.; Gabbiani, C.; Coronnello, M.; Mini, E.; Messori, L. Structural and solution chemistry, antiproliferative effects, and DNA and protein binding properties of a series of dinuclear gold (III) compounds with bipyridyl ligands. J. Med. Chem. 2006, 49, 5524–5531. [Google Scholar] [CrossRef]
  105. Gabbiani, C.; Casini, A.; Messori, L.; Guerri, A.; Cinellu, M.A.; Minghetti, G.; Corsini, M.; Rosani, C.; Zanello, P.; Arca, M. Structural characterization, solution studies, and DFT calculations on a series of binuclear gold (III) oxo complexes: Relationships to biological properties. Inorg. Chem. 2008, 47, 2368–2379. [Google Scholar] [CrossRef] [PubMed]
  106. Messori, L.; Abbate, F.; Marcon, G.; Orioli, P.; Fontani, M.; Mini, E.; Mazzei, T.; Carotti, S.; O’Connell, T.; Zanello, P. Gold (III) complexes as potential antitumor agents: Solution chemistry and cytotoxic properties of some selected gold (III) compounds. J. Med. Chem. 2000, 43, 3541–3548. [Google Scholar] [CrossRef]
  107. Clavier, H.; Nolan, S.P. Percent buried volume for phosphine and N-heterocyclic carbene ligands: Steric properties in organometallic chemistry. ChemComm 2010, 46, 841–861. [Google Scholar]
  108. Tolbatov, I.; Marrone, A. Reaction of dirhodium and diruthenium paddlewheel tetraacetate complexes with nucleophilic protein sites: A computational study. Inorg. Chim. Acta 2021, 530, 120684. [Google Scholar] [CrossRef]
  109. Tolbatov, I.; Marrone, A. Computational strategies to model the interaction and the reactivity of biologically-relevant transition metal complexes. Inorg. Chim. Acta 2021, 530, 120686. [Google Scholar] [CrossRef]
Figure 1. Overall structure of thaumatin-Au(NHC)Cl (with NHC = 1-butyl-3-methyl-imidazole-2-ylidene) (see ref. [19]).
Figure 1. Overall structure of thaumatin-Au(NHC)Cl (with NHC = 1-butyl-3-methyl-imidazole-2-ylidene) (see ref. [19]).
Molecules 26 07600 g001
Figure 2. Most common ligands in Au(I) antitumor complexes: phosphine with R = Ph, Me, Et (a); thiosugar (b); NHC complexes imidazol-2-ylidene (c) and benzimidazol-2-ylidene with R = Me, Et (d); alkynyl (e); thiosemicarbazone (f) with R = phenyl, p-methoxyphenyl, thiophene, furan, pyridinyl, and p-hydroxyphenyl. Atoms directly bound to Au(I) are in bold.
Figure 2. Most common ligands in Au(I) antitumor complexes: phosphine with R = Ph, Me, Et (a); thiosugar (b); NHC complexes imidazol-2-ylidene (c) and benzimidazol-2-ylidene with R = Me, Et (d); alkynyl (e); thiosemicarbazone (f) with R = phenyl, p-methoxyphenyl, thiophene, furan, pyridinyl, and p-hydroxyphenyl. Atoms directly bound to Au(I) are in bold.
Molecules 26 07600 g002
Figure 3. Au(I)-based structures, (1) auranofin, (2) Et3PAuCl, (3) auranofin iodide, and (4) auranofin with phenyl substituents.
Figure 3. Au(I)-based structures, (1) auranofin, (2) Et3PAuCl, (3) auranofin iodide, and (4) auranofin with phenyl substituents.
Molecules 26 07600 g003
Figure 4. A three-dimensional model of one subunit of Schistosoma mansoni thioredoxin-glutathione reductase, SmTGR, is shown with the three gold-binding sites. The glutaredoxin domain of TGR is represented in green, whereas the thioredoxin domain is shown in red, and the bound flavin is also highlighted. (see ref. [46]).
Figure 4. A three-dimensional model of one subunit of Schistosoma mansoni thioredoxin-glutathione reductase, SmTGR, is shown with the three gold-binding sites. The glutaredoxin domain of TGR is represented in green, whereas the thioredoxin domain is shown in red, and the bound flavin is also highlighted. (see ref. [46]).
Molecules 26 07600 g004
Figure 5. Gold(I) bis-NHC-carbene complexes (5, 6, 7a7f).
Figure 5. Gold(I) bis-NHC-carbene complexes (5, 6, 7a7f).
Molecules 26 07600 g005
Figure 6. Gold(I) complexes with NHC carbenes (8, 9).
Figure 6. Gold(I) complexes with NHC carbenes (8, 9).
Molecules 26 07600 g006
Figure 7. Gold(I) complexes with NHC carbenes (10a10b, 11a11c, 12a12e).
Figure 7. Gold(I) complexes with NHC carbenes (10a10b, 11a11c, 12a12e).
Molecules 26 07600 g007
Figure 8. Gold(I) complexes with a 1,3-diethylbenzimidazol-2-ylidene N-heterocyclic carbene (13a13c, 14a14d).
Figure 8. Gold(I) complexes with a 1,3-diethylbenzimidazol-2-ylidene N-heterocyclic carbene (13a13c, 14a14d).
Molecules 26 07600 g008
Figure 9. The crystal structure of a G-quadruplex unit bound to one gold complex in the [Au(NHC)2]+/Tel24 adduct [66]. The following color code has been used to identify each single atom: magenta (Au+),yellow (K+), green (Sr2+) and orange (C of carbene ligands). (from PDB, code 6H5R).
Figure 9. The crystal structure of a G-quadruplex unit bound to one gold complex in the [Au(NHC)2]+/Tel24 adduct [66]. The following color code has been used to identify each single atom: magenta (Au+),yellow (K+), green (Sr2+) and orange (C of carbene ligands). (from PDB, code 6H5R).
Molecules 26 07600 g009
Figure 10. Gold(I) complexes with two (15a) and one (15b) methylcaffeine-ylidene group.
Figure 10. Gold(I) complexes with two (15a) and one (15b) methylcaffeine-ylidene group.
Molecules 26 07600 g010
Figure 11. Gold(I) complexes with two (16a) and one (16b) aryl-thiosemicarbazone moieties.
Figure 11. Gold(I) complexes with two (16a) and one (16b) aryl-thiosemicarbazone moieties.
Molecules 26 07600 g011
Figure 12. Gold(I) phosphine complexes (17a17f).
Figure 12. Gold(I) phosphine complexes (17a17f).
Molecules 26 07600 g012
Figure 13. Gold(I) complexes with tertiary phosphine and 3-methyladamantane–1,3,4-oxadiazole-2-thione (18a, 18b) and 5-adamantyl-1,3-thiazolidine-2-thione (18c, 18d).
Figure 13. Gold(I) complexes with tertiary phosphine and 3-methyladamantane–1,3,4-oxadiazole-2-thione (18a, 18b) and 5-adamantyl-1,3-thiazolidine-2-thione (18c, 18d).
Molecules 26 07600 g013
Scheme 1. The lowest energy reaction path for a typical Au(III)-based complex is composed of two consecutive steps: (1) substitution of chloride by the nucleophile and (2) reduction of the resulting Au(III) complex to the corresponding Au(I) derivative with the opening of the chelate ring.
Scheme 1. The lowest energy reaction path for a typical Au(III)-based complex is composed of two consecutive steps: (1) substitution of chloride by the nucleophile and (2) reduction of the resulting Au(III) complex to the corresponding Au(I) derivative with the opening of the chelate ring.
Molecules 26 07600 sch001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tolbatov, I.; Marrone, A.; Coletti, C.; Re, N. Computational Studies of Au(I) and Au(III) Anticancer MetalLodrugs: A Survey. Molecules 2021, 26, 7600. https://doi.org/10.3390/molecules26247600

AMA Style

Tolbatov I, Marrone A, Coletti C, Re N. Computational Studies of Au(I) and Au(III) Anticancer MetalLodrugs: A Survey. Molecules. 2021; 26(24):7600. https://doi.org/10.3390/molecules26247600

Chicago/Turabian Style

Tolbatov, Iogann, Alessandro Marrone, Cecilia Coletti, and Nazzareno Re. 2021. "Computational Studies of Au(I) and Au(III) Anticancer MetalLodrugs: A Survey" Molecules 26, no. 24: 7600. https://doi.org/10.3390/molecules26247600

APA Style

Tolbatov, I., Marrone, A., Coletti, C., & Re, N. (2021). Computational Studies of Au(I) and Au(III) Anticancer MetalLodrugs: A Survey. Molecules, 26(24), 7600. https://doi.org/10.3390/molecules26247600

Article Metrics

Back to TopTop