Tyrosine kinase inhibitors (TKIs) for ovarian cancer treatment: from organic to inorganic chemotherapeutics towards selectivity-a perspective overview.

Ovarian cancer (OC) is a lethal gynecologic cancer in industrialized countries. Treatments for OC include the surgical removal and chemotherapy. In the last decades, improvements have been made in the surgery technologies, drug combinations and administration protocols, and in diagnosis. However, mortality from OC is still high owing to recurrences and insurgence of drug resistance. Accordingly, it is urgent the development of novel agents capable to effectively target OC. In this respect, tyrosine kinase inhibitors (TKIs) may play an important role. Most of TKIs developed and tested so far are organic. However, owing to their chemical versatility, also metals can be exploited to design selective and potent TKIs. We provide a short and easy-to-read overview on the main organic TKIs with a summary of those that entered clinical trials. Additionally, we describe the potential of metal-based TKIs, focusing on this overlooked family of compounds that may significantly contribute towards the concept of precision-medicine.

Mechanism of Action of Antitumor Au(I) N-Heterocyclic Carbene Complexes: A Computational Insight on the Targeting of TrxR Selenocysteine.

The targeting of human thioredoxin reductase is widely recognized to be crucially involved in the anticancer properties of several metallodrugs, including Au(I) complexes. In this study, the mechanism of reaction between a set of five N-heterocyclic carbene Au(I) complexes and models of the active Sec residue in human thioredoxin reductase was investigated by means of density functional theory approaches. The study was specifically addressed to the kinetics and thermodynamics of the tiled process by aiming at elucidating and explaining the differential inhibitory potency in this set of analogous Au(I) bis-carbene complexes. While the calculated free energy profile showed a substantially similar reactivity, we found that the binding of these Au(I) bis-carbene at the active CysSec dyad in the TrxR enzyme could be subjected to steric and orientational restraints, underlining both the approach of the bis-carbene scaffold and the attack of the selenol group at the metal center. A new and detailed mechanistic insight to the anticancer activity of these Au(I) organometallic complexes was thus provided by consolidating the TrxR targeting paradigm.

Inorganic Drugs as a Tool for Protein Structure Solving and Studies on Conformational Changes.

Inorganic drugs are capable of tight interactions with proteins through coordination towards aminoacidic residues, and this feature is recognized as a key aspect for their pharmacological action. However, the "protein metalation process" is exploitable for solving the phase problem and structural resolution. In fact, the use of inorganic drugs bearing specific metal centers and ligands capable to drive the binding towards the desired portions of the protein target could represent a very intriguing and fruitful strategy. In this context, a theoretical approach may further contribute to solve protein structures and their refinement. Here, we delineate the main features of a reliable experimental-theoretical integrated approach, based on the use of metallodrugs, for protein structure solving.

Mechanistic Evaluations of the Effects of Auranofin Triethylphosphine Replacement with a Trimethylphosphite Moiety.

Auranofin, a gold(I)-based complex, is under clinical trials for application as an anticancer agent for the treatment of nonsmall-cell lung cancer and ovarian cancer. In the past years, different derivatives have been developed, modifying gold linear ligands in the search for new gold complexes endowed with a better pharmacological profile. Recently, a panel of four gold(I) complexes, inspired by the clinically established compound auranofin, was reported by our research group. As described, all compounds possess an [Au{P(OMe)3}]+ cationic moiety, in which the triethylphosphine of the parent compound auranofin was replaced with an oxygen-rich trimethylphosphite ligand. The gold(I) linear coordination geometry was complemented by Cl-, Br-, I-, and the auranofin-like thioglucose tetraacetate ligand. As previously reported, despite their close similarity to auranofin, the panel compounds exhibited some peculiar and distinctive features, such as lower log P values which can induce relevant differences in the overall pharmacokinetic profiles. To get better insight into the P-Au strength and stability, an extensive study was carried out for relevant biological models, including three different vasopressin peptide analogues and cysteine, using 31P NMR and LC-ESI-MS. A DFT computational study was also carried out for a better understanding of the theoretical fundamentals of the disclosed differences with regard to triethylphosphine parent compounds.

Tetrazo[1,2-b]indazoles: Straightforward Access to Nitrogen-Rich Polyaromatics from s-Tetrazines.

The straightforward access to a new class of aza-polyaromatics is reported. Starting from readily available fluorinated s-tetrazine, a cyclization process with azide leads to the formation of an unprecedented tetrazo[1,2-b]indazole or a bis-tetrazo[1,2-b]indazole (cis and trans conformers). Based on the new nitrogen core, further N-directed palladium-catalyzed ortho-C-H bond functionalization allows the introduction of halides or acetates. The physicochemical properties of these compounds were studied by a joint experimental/theoretical approach. The tetrazo[1,2-b]indazoles display solid-state π-stacking, low reduction potential, absorption in the visible range up to the near-infrared, and intense fluorescence, depending on the molecular structure.

Kinetics of Reactions of Dirhodium and Diruthenium Paddlewheel Tetraacetate Complexes with Nucleophilic Protein Sites: Computational Insights.

Recently, dirhodium and diruthenium paddlewheel complexes have drawn attention as perspective anticancer drugs. In this study, the kinetics of reaction of typical paddlewheel scaffolds Rh2(µ-O2CCH3)4(H2O)2, Ru2(µ-O2CCH3)4(H2O)Cl, and [Ru2(µ-O2CCH3)4(HO)Cl]- with protein nucleophiles were investigated by means of the density functional theory. The substitution of axial ligands─water and chloride─by the models of protein residue side chains was analyzed, revealing the binding selectivity displayed by these paddlewheel metal scaffolds. The substitution of water is under a thermodynamic control, in which, although the Arg, Cys-, and Sec- residues are the most favorable, their binding is expected to be scarcely selective in a biological context. On the other hand, the replacement of the axial water with a more stable hydroxo ligand induces the chloride substitution in diRu complexes, which also targets Arg, Cys-, and Sec-, although with a moderately higher activation barrier for any examined protein residue. Additionally, the carried out characterization of the geometrical parameters of the transition states permitted determination of the impact of an increased steric hindrance of diRh and diRu complexes on their protein site selectivity. This study corroborates the idea of the substitution of the acetate ligands with biologically active, but more hindering, carboxylate ligands, in order to yield dual acting metallodrugs. This study allows us to assume that the delivery of diRu paddlewheel complexes in their monoanionic form [Ru2(µ-O2CR)4(OH)Cl]- decorated by the bulky substituents R may constitute an approach to augment the selectivity toward anticancer targets, such as TrxR in tumor cells, although under the condition that such a selectivity is operative only in high pH conditions.

Reactivity of N-Heterocyclic Carbene Half-Sandwich Ru-, Os-, Rh-, and Ir-Based Complexes with Cysteine and Selenocysteine: A Computational Study.

The structure and the reactivity of four half-sandwich metal complexes of RuII, OsII, RhIII, and IrIII were investigated by means of density functional theory approaches. These piano-stool complexes, grouped in cym-bound complexes, RuII(cym)(dmb)Cl2, 1, and OsII(cym)(dmb)Cl2, 2, and Cp*-bound complexes, RhIII(Cp*)(dmb)Cl2, 3, and IrIII(Cp*)(dmb)Cl2, 4, with cym = η6-p-cymene, Cp* = η5-pentamethylcyclopentadienyl, and dmb = 1,3-dimethylbenzimidazol-2-ylidene, were recently proposed as anticancer metallodrugs that preferably target Cys- or Sec-containing proteins. Thus, density functional theory calculations were performed here to characterize in detail the thermodynamics and the kinetics underlining the targeting of these metallodrugs at either neutral or anionic Cys and Sec side chains. Calculations evidenced that all these complexes preferably target at Cys or Sec via chloro exchange, although cym-bound and Cp*-bound complexes resulted to be more prone to bind at neutral or anionic forms, respectively, of these soft protein sites. Further decomposition analyses of the activation free energies for the reaction between 1-4 complexes and either Cys or Sec, paralleled with the comparison among the optimized transition-state structures, allowed us to spotlight the significant role played by solvation in determining the overall reactivity and selectivity expected for these prototypical metallodrugs.

Structural Reshaping of the Zinc-Finger Domain of the SARS-CoV-2 nsp13 Protein Using Bismuth(III) Ions: A Multilevel Computational Study.

The identification of novel therapeutics against the pandemic SARS-CoV-2 infection is an indispensable new address of current scientific research. In the search for anti-SARS-CoV-2 agents as alternatives to the vaccine or immune therapeutics whose efficacy naturally degrades with the occurrence of new variants, the salts of Bi3+ have been found to decrease the activity of the Zn2+-dependent non-structural protein 13 (nsp13) helicase, a key component of the SARS-CoV-2 molecular tool kit. Here, we present a multilevel computational investigation based on the articulation of DFT calculations, classical MD simulations, and MIF analyses, focused on the examination of the effects of Bi3+/Zn2+ exchange on the structure and molecular interaction features of the nsp13 protein. Our calculations confirmed that Bi3+ ions can replace Zn2+ in the zinc-finger metal centers and cause slight but appreciable structural modifications in the zinc-binding domain of nsp13. Nevertheless, by employing an in-house-developed ATOMIF tool, we evidenced that such a Bi3+/Zn2+ exchange may decrease the extension of a specific hydrophobic portion of nsp13, responsible for the interaction with the nsp12 protein. The present study provides for a detailed, atomistic insight into the potential anti-SARS-CoV-2 activity of Bi3+ and, more generally, evidences the hampering of the nsp13-nsp12 interaction as a plausible therapeutic strategy.

Reactions of Arsenoplatin-1 with Protein Targets: A Combined Experimental and Theoretical Study.

Arsenoplatin-1 (AP-1) is a dual-action anticancer metallodrug with a promising pharmacological profile that features the simultaneous presence of a cisplatin-like center and an arsenite center. We investigated its interactions with proteins through a joint experimental and theoretical approach. The reactivity of AP-1 with a variety of proteins, including carbonic anhydrase (CA), superoxide dismutase (SOD), myoglobin (Mb), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), and human serum albumin (HSA), was analyzed by means of electrospray ionization mass spectrometry (ESI MS) measurements. In accordance with previous observations, ESI MS experiments revealed that the obtained metallodrug-protein adducts originated from the binding of the [(AP-1)-Cl]+ fragment to accessible protein residues. Remarkably, in two cases, i.e., Mb and GAPDH, the formation of a bound metallic fragment that lacked the arsenic center was highlighted. The reactions of AP-1 with various nucleophiles side chains of neutral histidine, methionine, cysteine, and selenocysteine, in neutral form as well as cysteine and selenocysteine in anionic form, were subsequently analyzed through a computational approach. We found that the aquation of AP-1 is energetically disfavored, with a reaction free energy of +19.2 kcal/mol demonstrating that AP-1 presumably attacks its biological targets through the exchange of the chloride ligand. The theoretical analysis of thermodynamics and kinetics for the ligand-exchange processes of AP-1 with His, Met, Cys, Sec, Cys-, and Sec- side chain models unveils that only neutral histidine and deprotonated cysteine and selenocysteine are able to effectively replace the chloride ligand in AP-1.

Medicinal Hypervalent Tellurium Prodrugs Bearing Different Ligands: A Comparative Study of the Chemical Profiles of AS101 and Its Halido Replaced Analogues.

Ammonium trichloro (dioxoethylene-O,O') tellurate (AS101) is a potent immunomodulator prodrug that, in recent years, entered various clinical trials and was tested for a variety of potential therapeutic applications. It has been demonstrated that AS101 quickly activates in aqueous milieu, producing TeOCl3-, which likely represents the pharmacologically active species. Here we report on the study of the activation process of AS101 and of two its analogues. After the synthesis and characterization of AS101 and its derivatives, we have carried out a comparative study through a combined experimental and computational analysis. Based on the obtained results, we describe here, for the first time, the detailed reaction that AS101 and its bromido- and iodido-replaced analogues undergo in presence of water, allowing the conversion of the original molecule to the likely true pharmacophore. Interestingly, moving down in the halogens' group we observed a higher tendency to react, attributable to the ligands' effect. The chemical and mechanistic implications of these meaningful differences are discussed.

Anti-Staphylococcal Activity of the Auranofin Analogue Bearing Acetylcysteine in Place of the Thiosugar: An Experimental and Theoretical Investigation.

Auranofin (AF, hereafter) is an orally administered chrysotherapeutic agent approved for the treatment of rheumatoid arthritis that is being repurposed for various indications including bacterial infections. Its likely mode of action involves the impairment of the TrxR system through the binding of the pharmacophoric cation [AuPEt3]+. Accordingly, a reliable strategy to expand the medicinal profile of AF is the replacement of the thiosugar moiety with different ligands. Herein, we aimed to prepare the AF analogue bearing the acetylcysteine ligand (AF-AcCys, hereafter) and characterize its anti-staphylococcal activity. Biological studies revealed that AF-AcCys retains an antibacterial effect superimposable with that of AF against Staphylococcus aureus, whereas it is about 20 times less effective against Staphylococcus epidermidis. Bioinorganic studies confirmed that upon incubation with human serum albumin, AF-AcCys, similarly to AF, induced protein metalation through the [AuPEt3]+ fragment. Additionally, AF-AcCys appeared capable of binding the dodecapeptide Ac-SGGDILQSGCUG-NH2, corresponding to the tryptic C-terminal fragment (488-499) of hTrxR. To shed light on the pharmacological differences between AF and AF-AcCys, we carried out a comparative experimental stability study and a theoretical estimation of bond dissociation energies, unveiling the higher strength of the Au-S bond in AF-AcCys. From the results, it emerged that the lower lipophilicity of AF-AcCys with respect to AF could be a key feature for its different antibacterial activity. The differences and similarities between AF and AF-AcCys are discussed, alongside the opportunities and consequences that chemical structure modifications imply.

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

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.

Probing the Paradigm of Promiscuity for N-Heterocyclic Carbene Complexes and their Protein Adduct Formation.

Metal complexes can be considered a "paradigm of promiscuity" when it comes to their interactions with proteins. They often form adducts with a variety of donor atoms in an unselective manner. We have characterized the adducts formed between a series of isostructural N-heterocyclic carbene (NHC) complexes with Ru, Os, Rh, and Ir centers and the model protein hen egg white lysozyme by X-ray crystallography and mass spectrometry. Distinctive behavior for the metal compounds was observed with the more labile Ru and Rh complexes targeting mainly a surface l-histidine moiety through cleavage of p-cymene or NHC co-ligands, respectively. In contrast, the more inert Os and Ir derivatives were detected abundantly in an electronegative binding pocket after undergoing ligand exchange of a chlorido ligand for an amino acid side chain. Computational studies supported the binding profiles and hinted at the role of the protein microenvironment for metal complexes eliciting selectivity for specific binding sites on the protein.

Insight into the Substitution Mechanism of Antitumor Au(I) N-Heterocyclic Carbene Complexes by Cysteine and Selenocysteine.

We carried out a detailed theoretical study on the mechanism of the carbene ligand substitution by cysteine and selenocysteine residues in an Au(I) bis-N-heterocyclic carbene complex in order to model the initial stages of the mechanism of action of this promising class of antitumor metallodrug. Both neutral and deprotonated capped Cys and Sec species were considered as possible nucleophiles in the ligand exchange reaction on the metal center to model the corresponding protein side chains. Energies and geometric structures of the possible transition states and reactant- and product-adducts involved in the substitution process have been calculated using density functional theory and local MP2. Reaction and activation enthalpies and free energies have been evaluated and indicate a slightly exothermic and exergonic process with reasonably low barriers, 21.3 and 19.6 kcal mol-1, respectively, for capped Cys and Sec, in good agreement with the experimental data available for the reaction with free amino acids. The results suggest a mechanism for the ligand exchange reaction involving an anionic thiolate or selenothiolate species coupled to an explicit proton transfer to the leaving carbene from the acidic component of the buffer. The presence of a buffer is necessary both in in vitro experiments and under physiological conditions, and its proton reservoir behavior reveals the importance of the environmental effects in carbene substitution by biological nucleophiles.

Determinants of the Lead(II) Affinity in pbrR Protein: A Computational Study.

Investigation of the diverse evolutionary developed mechanisms enabling bacteria to maintain homeostasis and to be resistant to lead is crucial for the discovery of novel strategies for isolation of this highly toxic metal and its subsequent elimination from contaminated environments. The metalloregulatory protein pbrR and its homologues that were identified in the Cupriavidus metallidurans CH34 chromosome are the only characterized natural metalloproteins that have a special affinity toward Pb(II) and that bind it with at least a 1000-fold selectivity over other heavy metals. The X-ray structures of apo and Pb(II)-bound pbrR have been recently reported. In the present study, the binding of Pb(II) at pbrR was investigated by means of multiscale computational modeling. Molecular dynamics simulations substantiated how conformations amenable for the Pb(II) complexation through the tris-cysteine motif are formed from the antiparallel coiled-coil packing interaction of two dimerization helices of two pbrR monomers, and the phase space of apo-pbrR has been extensively sampled. Hybrid quantum mechanics/molecular mechanics (QM/MM) calculations on metal-bound structures of pbrR also allowed us to determine the most probable protonation state for the lead binding motif and evaluate the structural features mostly affecting the Pb(II) coordination in this protein. In agreement with available experimental data, we found that pbrR may control its Pb(II) affinity, probably, by conformational changes that affect the distance between Cys78' and Cys122 and their protonation states, thus being able to switch on the Pb(II) sequestration/release-prone states in response to external stimuli. The protein structure enveloping the metal binding motif favors the thiol-thiolate-thiolate protonation state of Pb(II)-pbrR, thus probably enhancing the binding selectivity for Pb(II), compared to other metal ions.

An Insight on the Gold(I) Affinity of golB Protein via Multilevel Computational Approaches.

Several bacterial species have evolutionary developed protein systems specialized in the control of intracellular gold ion concentration. In order to prevent the detrimental consequences that may be induced even at very low concentrations, bacteria such as Salmonella enterica and Cupriavidus metallidurans utilize Au-specific merR-type transcriptional regulators that detect these toxic ions and control the expression of specific resistance factors. Among these highly specialized proteins, golB has been investigated in depth, and X-ray structures of both apo and Au(I)-bound golB have been recently reported. Here, the binding of Au(I) at golB was investigated by means of multilevel computational approaches. Molecular dynamics simulations evidenced how conformations amenable for the Au(I) chelation through the Cys-XX-Cys motif on helix 1 are extensively sampled in the phase space of apo-golB. Hybrid QM/MM calculations on metal-bound structures of golB also allowed to characterize the most probable protonation state for gold binding motif and to assess the structural features mostly influencing the Au(I) coordination in this protein. Consistently with experimental evidence, we found that golB may control its Au(I) affinity by conformational changes that affect the distance between Cys10 and Cys13, thus being able to switch between the Au(I) sequestration/release-prone states in response to external stimuli. The protein structure enveloping the metal binding motif favors the thiol-thiolate protonation state of Au(I)-golB, thus probably enhancing the binding selectivity for Au(I) compared to other cations.

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

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

Insight into the Electrochemical Reduction Mechanism of Pt(IV) Anticancer Complexes.

We carried out a theoretical study on the mechanism of electrochemical reduction of the prototypical platinum(IV) anticancer complex [Pt(NH3)2(CH3COO)2Cl2] to the corresponding platinum(II) [Pt(NH3)2Cl2] derivative. Energies and geometric structures of the original Pt(IV) complex and all possible Pt(III) and Pt(II) intermediates and transition states involved in the reduction process have been calculated using density functional theory and Møller-Plesset perturbation theory. This study allowed us to formulate a detailed mechanism for the two-electron reduction of the [PtIV(NH3)2(CH3COO)2Cl2] complex. The results show that, in agreement with the experimental evidence from cyclic voltammetry, the initial one-electron reduction of the [PtIV(NH3)2(CH3COO)2Cl2] complex occurs through a stepwise mechanism via a metastable hexacoordinated platinum(III) [PtIII(NH3)2(CH3COO)2Cl2]- intermediate and a subsequent acetate ligand detachment with an activation free energy of 5.1 kcal mol-1. On the other hand, the second electron reduction of the resulting pentacoordinated [PtIII(NH3)2(CH3COO)Cl2] species occurs through a barrierless concerted process to the final [PtII(NH3)2Cl2] derivative.

Monocarbon cationic cluster yields from N2/CH4 mixtures embedded in He nanodroplets and their calculated binding energies.

The formation of monocarbon cluster ions has been investigated by electron ionization mass spectrometry of cold helium nanodroplets doped with nitrogen/methane mixtures. Ion yields for two groups of clusters, CHmN2(+) or CHmN4(+), were determined for mixtures with different molecular ratios of CH4. The possible geometrical structures of these clusters were analyzed using electronic structure computations. Little correlation between the ion yields and the associated binding energies has been observed indicating that in most cases kinetic control is more important than thermodynamic control for forming the clusters.

Diruthenium Paddlewheel Complexes Attacking Proteins: Axial versus Equatorial Coordination.

Metallodrugs are an important group of medicinal agents used for the treatment of various diseases ranging from cancers to viral, bacterial, and parasitic diseases. Their distinctive features include the availability of a metal centre, redox activity, as well as the ability to multitarget. Diruthenium paddlewheel complexes are an intensely developing group of metal scaffolds, which can securely coordinate bidentate xenobiotics and transport them to target tissues, releasing them by means of substitution reactions with biomolecular nucleophiles. It is of the utmost importance to gain a complete comprehension of which chemical reactions happen with them in physiological milieu to design novel drugs based on these bimetallic scaffolds. This review presents the data obtained in experiments and calculations, which clarify the chemistry these complexes undergo once administered in the proteic environment. This study demonstrates how diruthenium paddlewheel complexes may indeed embody a new paradigm in the design of metal-based drugs of dual-action by presenting and discussing the protein metalation by these complexes.