Multifunctional Organometallic Compounds Active against Infective Trypanosomes: Ru(II) Ferrocenyl Derivatives with Two Different Bioactive Ligands.

Human African trypanosomiasis (HAT, sleeping sickness) and American trypanosomiasis (Chagas disease) are endemic zoonotic diseases caused by genomically related trypanosomatid protozoan parasites (Trypanosoma brucei and Trypanosoma cruzi, respectively). Just a few old drugs are available for their treatment, with most of them sharing poor safety, efficacy, and pharmacokinetic profiles. Only fexinidazole has been recently incorporated into the arsenal for the treatment of HAT. In this work, new multifunctional Ru(II) ferrocenyl compounds were rationally designed as potential agents against these pathogens by including in a single molecule 1,1'-bis(diphenylphosphino)ferrocene (dppf) and two bioactive bidentate ligands: pyridine-2-thiolato-1-oxide ligand (mpo) and polypyridyl ligands (NN). Three [Ru(mpo)(dppf)(NN)](PF6) compounds and their derivatives with chloride as a counterion were synthesized and fully characterized in solid state and solution. They showed in vitro activity on bloodstream T. brucei (EC50 = 31-160 nM) and on T. cruzi trypomastigotes (EC50 = 190-410 nM). Compounds showed the lowest EC50 values on T. brucei when compared to the whole set of metal-based compounds previously developed by us. In addition, several of the Ru compounds showed good selectivity toward the parasites, particularly against the highly proliferative bloodstream form of T. brucei. Interaction with DNA and generation of reactive oxygen species (ROS) were ruled out as potential targets and modes of action of the Ru compounds. Biochemical assays and in silico analysis led to the insight that they are able to inhibit the NADH-dependent fumarate reductase from T. cruzi. One representative hit induced a mild oxidation of low molecular weight thiols in T. brucei. The compounds were stable for at least 72 h in two different media and more lipophilic than both bioactive ligands, mpo and NN. An initial assessment of the therapeutic efficacy of one of the most potent and selective candidates, [Ru(mpo)(dppf)(bipy)]Cl, was performed using a murine infection model of acute African trypanosomiasis. This hit compound lacks acute toxicity when applied to animals in the dose/regimen described, but was unable to control parasite proliferation in vivo, probably because of its rapid clearance or low biodistribution in the extracellular fluids. Future studies should investigate the pharmacokinetics of this compound in vivo and involve further research to gain deeper insight into the mechanism of action of the compounds.

Dismantling and Rebuilding the Trisulfide Cofactor Demonstrates Its Essential Role in Human Sulfide Quinone Oxidoreductase.

Sulfide quinone oxidoreductase (SQOR) catalyzes the first step in sulfide clearance, coupling H2S oxidation to coenzyme Q reduction. Recent structures of human SQOR revealed a sulfur atom bridging the SQOR active site cysteines in a trisulfide configuration. Here, we assessed the importance of this cofactor using kinetic, crystallographic, and computational modeling approaches. Cyanolysis of SQOR proceeds via formation of an intense charge transfer complex that subsequently decays to eliminate thiocyanate. We captured a disulfanyl-methanimido thioate intermediate in the SQOR crystal structure, revealing how cyanolysis leads to reversible loss of SQOR activity that is restored in the presence of sulfide. Computational modeling and MD simulations revealed an ∼105-fold rate enhancement for nucleophilic addition of sulfide into the trisulfide versus a disulfide cofactor. The cysteine trisulfide in SQOR is thus critical for activity and provides a significant catalytic advantage over a cysteine disulfide.

The Chemical Basis of Thiol Addition to Nitro-conjugated Linoleic Acid, a Protective Cell-signaling Lipid.

Nitroalkene fatty acids are formed in vivo and exert protective and anti-inflammatory effects via reversible Michael addition to thiol-containing proteins in key signaling pathways. Nitro-conjugated linoleic acid (NO2-CLA) is preferentially formed, constitutes the most abundant nitrated fatty acid in humans, and contains two carbons that could potentially react with thiols, modulating signaling actions and levels. In this work, we examined the reactions of NO2-CLA with low molecular weight thiols (glutathione, cysteine, homocysteine, cysteinylglycine, and ß-mercaptoethanol) and human serum albumin. Reactions followed reversible biphasic kinetics, consistent with the presence of two electrophilic centers in NO2-CLA located on the ß- and δ-carbons with respect to the nitro group. The differential reactivity was confirmed by computational modeling of the electronic structure. The rates (kon and koff) and equilibrium constants for both reactions were determined for different thiols. LC-UV-Visible and LC-MS analyses showed that the fast reaction corresponds to ß-adduct formation (the kinetic product), while the slow reaction corresponds to the formation of the δ-adduct (the thermodynamic product). The pH dependence of the rate constants, the correlation between intrinsic reactivity and thiol pKa, and the absence of deuterium solvent kinetic isotope effects suggested stepwise mechanisms with thiolate attack on NO2-CLA as rate-controlling step. Computational modeling supported the mechanism and revealed additional features of the transition states, anionic intermediates, and final neutral products. Importantly, the detection of cysteine-δ-adducts in human urine provided evidence for the biological relevance of this reaction. Finally, human serum albumin was found to bind NO2-CLA both non-covalently and to form covalent adducts at Cys-34, suggesting potential modes for systemic distribution. These results provide new insights into the chemical basis of NO2-CLA signaling actions.

Reaction of Hydrogen Sulfide with Disulfide and Sulfenic Acid to Form the Strongly Nucleophilic Persulfide.

Hydrogen sulfide (H2S) is increasingly recognized to modulate physiological processes in mammals through mechanisms that are currently under scrutiny. H2S is not able to react with reduced thiols (RSH). However, H2S, more precisely HS(-), is able to react with oxidized thiol derivatives. We performed a systematic study of the reactivity of HS(-) toward symmetric low molecular weight disulfides (RSSR) and mixed albumin (HSA) disulfides. Correlations with thiol acidity and computational modeling showed that the reaction occurs through a concerted mechanism. Comparison with analogous reactions of thiolates indicated that the intrinsic reactivity of HS(-) is 1 order of magnitude lower than that of thiolates. In addition, H2S is able to react with sulfenic acids (RSOH). The rate constant of the reaction of H2S with the sulfenic acid formed in HSA was determined. Both reactions of H2S with disulfides and sulfenic acids yield persulfides (RSSH), recently identified post-translational modifications. The formation of this derivative in HSA was determined, and the rate constants of its reactions with a reporter disulfide and with peroxynitrite revealed that persulfides are better nucleophiles than thiols, which is consistent with the α effect. Experiments with cells in culture showed that treatment with hydrogen peroxide enhanced the formation of persulfides. Biological implications are discussed. Our results give light on the mechanisms of persulfide formation and provide quantitative evidence for the high nucleophilicity of these novel derivatives, setting the stage for understanding the contribution of the reactions of H2S with oxidized thiol derivatives to H2S effector processes.

Deconstructing the catalytic efficiency of peroxiredoxin-5 peroxidatic cysteine.

Human peroxiredoxin-5 (PRDX5) is a thiol peroxidase that reduces H2O2 10(5) times faster than free cysteine. To assess the influence of two conserved residues on the reactivity of the critical cysteine (C47), we determined the reaction rate constants of PRDX5, wild type (WT), T44V and R127Q with one substrate electrophile (H2O2) and a nonspecific electrophile (monobromobimane). We also studied the corresponding reactions of low molecular weight (LMW) thiolates in order to construct a framework against which we could compare our proteins. To obtain a detailed analysis of the structural and energetic changes involved in the reaction between WT PRDX5 and H2O2, we performed ONIOM quantum mechanics/molecular mechanics (QM/MM) calculations with a QM region including 60 atoms of substrate and active site described by the B3LYP density functional and the 6-31+G(d,p) basis set; the rest of the protein was included in the MM region. Brønsted correlations reveal that the absence of T44 can increase the general nucleophilicity of the C47 but decreases the specific reactivity toward H2O2 by a factor of 10(3). The R127Q mutation causes C47 to behave like a LMW thiolate in the two studied reactions. QM/MM results with WT PRDX5 showed that hydrogen bonds in the active site are the cornerstone of two effects that make catalysis possible: the enhancement of thiolate nucleophilicity upon substrate ingress and the stabilization of the transition state. In both effects, T44 has a central role. These effects occur in a precise temporal sequence that ensures that the selective nucleophilicity of C47 is available only for peroxide substrates.

VHL-P138R and VHL-L163R Novel Variants: Mechanisms of VHL Pathogenicity Involving HIF-Dependent and HIF-Independent Actions.

The von Hippel-Lindau (VHL) disease is an autosomal dominant cancer syndrome caused by mutations in the VHL tumor suppressor gene. VHL protein (pVHL) forms a complex (VBC) with Elongins B-C, Cullin2, and Rbx1. Although other functions have been discovered, the most described function of pVHL is to recognize and target hypoxia-inducible factor (HIF) for degradation. This work comprises the functional characterization of two novel variants of the VHL gene (P138R and L163R) that have been described in our center in patients with VHL disease by in vitro, in vivo, and in silico approaches. In vitro, we found that these variants have a significantly shorter half-life compared to wild-type VHL but still form a functional VBC complex. Altered fibronectin deposition was evidenced for both variants using immunofluorescence. In vivo studies revealed that both variants failed to suppress tumor growth. By means of molecular dynamics simulations, we inspected in silico the nature of the changes introduced by each variant in the VBC complex. We have demonstrated the pathogenicity of P138R and L163R novel variants, involving HIF-dependent and HIF-independent mechanisms. These results provide the basis for future studies regarding the impact of structural alterations on posttranslational modifications that drive pVHL's fate and functions.

Nitrogen dioxide solubility and permeation in lipid membranes.

Nitrogen dioxide (()NO(2)) is an important oxidant molecule in biology that is produced by several biological processes, and it is also an important air pollutant. It can oxidize proteins and lipids with important consequences on their biological functions. Despite its relevance, the interaction of ()NO(2) with the cell barrier, the lipid membrane, is poorly understood. For instance, can lipid membranes limit ()NO(2) diffusion? To estimate the permeability of lipid membranes to ()NO(2) it is necessary to learn more about its solubility in the lipid phase. However, experimental data on ()NO(2) solubility is very limited. To improve our knowledge on this matter, we used a mixed approach consisting in calculating the solubility of ()NO(2) and related diatomic and triatomic gases (()NO, O(2), CO(2), etc.) in different solvents using quantum calculations and Tomasi's Polarizable Continuum Model and validating and correcting these results using experimental data available for the related gases. This approach led to an estimated partition coefficient for ()NO(2) of 2.7 between n-octanol and water, and 1.5 between lipid membranes and water, meaning that ()NO(2) is a moderately hydrophobic molecule (less than ()NO, more than CO(2)). Based on the solubility-diffusion permeability theory, the permeability coefficient was estimated to be 5 cms(-1), up to 4000 times higher than that of peroxynitrous acid. It is concluded that lipid membranes are not significant barriers to ()NO(2) transport.

Understanding the mechanism of H2S oxidation by flavin-dependent sulfide oxidases: a DFT/IEF-PCM study.

In the last years, H2S has been recognized as a signaling molecule in mammals, which can synthesize and catabolize (by oxidation) such species. The latter process is accelerated by a sulfide:quinone oxidoreductase (SQR, E.C. 1.8.5.4), a flavin-dependent sulfide oxidase (FDSO). FDSOs catalyze electron transfer from H2S to an acceptor in catalytic cycles involving two phases: (I) reduction of FAD by H2S (SH-) and (II) electron transfer from FADH- to the electron acceptor. The first step of FAD reduction consists on the reaction of SH- with a catalytic disulfide at the active site of the enzyme, to yield a thiolate and a persulfide in the protein. This step is ca. 106 times faster than the analogous reaction with low-molecular-weight disulfides (LMWDs) and the causes of such extraordinary acceleration remain unknown. Using the IEF-PCM(ε ≈ 10)/M06-2X-D3/6-31+G(d,p) level of theory, we have modeled the reaction of SH- with a disulfide as located in a representative model of the active site extracted from a prokaryotic SQR, assessing the effects of partial covalent interactions (PCIs) between the leaving sulfur atom and flavin ring on the activation Gibbs free-energy barrier at 298 K (∆‡G298K). To also evaluate the importance of entropic penalties on the first step, we have modeled at the same level of theory the reaction of (bis)hydroxyethyl disulfide in aqueous solution, a LMWD for which experimental data is available. Our results show that PCIs between the leaving sulfur atom and the flavin group only have a minor effect (∆‡G298K reduced by 1.6 kcal mol-1) while compensating entropic penalties could have a much larger effect (up to 8.3 kcal mol-1). Finally, we also present here a first model of some of further steps in the phase I of the catalytic cycle as in mammalian FDSOs, providing some light about their detailed mechanism. Graphical abstract .

The thiol of human serum albumin: Acidity, microenvironment and mechanistic insights on its oxidation to sulfenic acid.

Human serum albumin (HSA) has a single reduced cysteine residue, Cys34, whose acidity has been controversial. Three experimental approaches (pH-dependence of reactivity towards hydrogen peroxide, ultraviolet titration and infrared spectroscopy) are used to determine that the pKa value in delipidated HSA is 8.1±0.2 at 37°C and 0.1M ionic strength. Molecular dynamics simulations of HSA in the sub-microsecond timescale show that while sulfur exposure to solvent is limited and fluctuating in the thiol form, it increases in the thiolate, stabilized by a persistent hydrogen-bond (HB) network involving Tyr84 and bridging waters to Asp38 and Gln33 backbone. Insight into the mechanism of Cys34 oxidation by H2O2 is provided by ONIOM(QM:MM) modeling including quantum water molecules. The reaction proceeds through a slightly asynchronous SN2 transition state (TS) with calculated Δ‡G and Δ‡H barriers at 298K of respectively 59 and 54kJmol-1 (the latter within chemical accuracy from the experimental value). A post-TS proton transfer leads to HSA-SO- and water as products. The structured reaction site cages H2O2, which donates a strong HB to the thiolate. Loss of this HB before reaching the TS modulates Cys34 nucleophilicity and contributes to destabilize H2O2. The lack of reaction-site features required for differential stabilization of the TS (positive charges, H2O2 HB strengthening) explains the striking difference in kinetic efficiency for the same reaction in other proteins (e.g. peroxiredoxins). The structured HB network surrounding HSA-SH with sequestered waters carries an entropic penalty on the barrier height. These studies contribute to deepen the understanding of the reactivity of HSA-SH, the most abundant thiol in human plasma, and in a wider perspective, provide clues on the key aspects that modulate thiol reactivity against H2O2.

Connecting proline and γ-aminobutyric acid in stressed plants through non-enzymatic reactions.

The accumulation of proline (Pro) in plants exposed to biotic/abiotic stress is a well-documented and conserved response in most vegetal species. Stress conditions induce the overproduction of reactive oxygen species which can lead to cellular damage. In vitro assays have shown that enzyme inactivation by hydroxyl radicals (·OH) can be avoided in presence of Pro, suggesting that this amino acid could act as an ·OH scavenger. We applied Density Functional Theory coupled with a polarizable continuum model to elucidate how Pro reacts with ·OH. In this work we suggest that Pro reacts favourably with ·OH by H-abstraction on the amine group. This reaction produces the spontaneous decarboxylation of Pro leading to the formation of pyrrolidin-1-yl. In turn, pyrrolidin-1-yl can easily be converted to Δ1-pyrroline, the substrate of the enzyme Δ1-pyrroline dehydrogenase, which produces γ-aminobutyric acid (GABA). GABA and Pro are frequently accumulated in stressed plants and several protective roles have been assigned to these molecules. Thereby we present an alternative non-enzymatic way to synthetize GABA under oxidative stress. Finally this work sheds light on a new beneficial role of Pro accumulation in the maintenance of photosynthetic activity.

Aromatic amine N-oxide organometallic compounds: searching for prospective agents against infectious diseases.

In search of prospective agents against infectious diseases, 1,1'-bis(diphenylphosphino)ferrocene pyridine-2-thiolato-1-oxide M(ii) hexafluorophosphate compounds [M(mpo)(dppf)](PF6), where M = palladium or platinum, were synthesized and fully characterized in the solid state and in solution using experimental and DFT computational techniques. The compounds are isomorphous and the M(ii) transition metal ions are in a nearly planar trapezoidal cis-coordination bound to the pyridine-2-thiolato-1-oxide (mpo) and to the 1,1'-bis(diphenylphosphino)ferrocene molecules, both acting as bidentate ligands. Both compounds showed high cytotoxic activity on Trypanosoma cruzi and Mycobacterium tuberculosis (MTB) and acceptable selectivities towards MTB, but good to excellent selectivity index values as anti-T. cruzi compounds. The inclusion of the ferrocene moiety (dppf ligand) improved the selectivity towards the parasite when compared to the previously reported [M(mpo)2] complexes. Related to the probable mechanism of action of the complexes, molecular docking studies on modelled T. cruzi NADH-fumarate reductase (TcFR) predicted that both be very good inhibitors of the enzyme. The effect of the compounds on the enzyme activity was experimentally confirmed using T. cruzi protein extracts. According to all obtained results, both [M(mpo)(dppf)](PF6) compounds could be considered prospective anti-trypanosomal agents that deserve further research.

Expanding the family of heteroleptic oxidovanadium(IV) compounds with salicylaldehyde semicarbazones and polypyridyl ligands showing anti-Trypanosoma cruzi activity.

Searching for prospective vanadium-based drugs for the treatment of Chagas disease, a new series of heteroleptic [V(IV)O(L-2H)(NN)] compounds was developed by including the lipophilic 3,4,7,8-tetramethyl-1,10-phenanthroline (tmp) NN ligand and seven tridentate salicylaldehyde semicarbazone derivatives (L1-L7). The compounds were characterized in the solid state and in solution. EPR spectroscopy suggests that the NN ligand is bidentate bound through both nitrogen donor atoms in an axial-equatorial mode. The EPR and (51)V-NMR spectra of aerated solutions at room temperature indicate that the compounds are stable to hydrolysis and that no significant oxidation of V(IV) to V(V) takes place at least in 24h. The complexes are more active in vitro against Trypanosoma cruzi, the parasite responsible for Chagas disease, than the reference drug Nifurtimox and most of them are more active than previously reported [V(IV)O(L-2H)(NN)] complexes of other NN co-ligands. Selectivity towards the parasite was analyzed using J-774 murine macrophages as mammalian cell model. Due to both, high activity and high selectivity, L2, L4, L5 and L7 complexes could be considered new hits for further drug development. Lipophilicity probably plays a relevant role in the bioactivity of the new compounds. The [V(IV)O(L-2H)(NN)] compounds were designed aiming DNA as potential molecular target. Therefore, the novel L1-L7 tmp complexes were screened by computational modeling, comparing their DNA-binding features with those of previously reported [V(IV)O(L-2H)(NN)] compounds with different NN co-ligands. Whereas all the complexes interact well with DNA, with binding modes and strength tuned in different extents by the NN and semicarbazone co-ligands, molecular docking suggests that the observed anti-T. cruzi activity cannot be explained upon DNA intercalation as the sole mechanism of action.

Molecular mechanisms for the reaction between (˙)OH radicals and proline: insights on the role as reactive oxygen species scavenger in plant stress.

The accumulation of proline (Pro) and overproduction of reactive oxygen species (ROS) by plants exposed to stress is well-documented. In vitro assays show that enzyme inactivation by hydroxyl radicals ((•)OH) can be avoided in the presence of Pro, suggesting this amino acid might act as a (•)OH scavenger. Although production of hydroxyproline (Hyp) has been hypothesized in connection with such antioxidant activity, no evidence on the detailed mechanism of scavenging has been reported. To elucidate whether and how Hyp might be produced, we used density functional theory calculations coupled to a polarizable continuum model to explore 27 reaction channels including H-abstraction by (•)OH and (•)OH/H2O addition. The structure and energetics of stable species and transition states for each reaction channel were characterized at the PCM-(U)M06/6-31G(d,p) level in aqueous solution. Evidence is found for a main pathway in which Pro scavenges (•)OH by successive H-abstractions (ΔG(‡,298) = 4.1 and 7.5 kcal mol(-1)) to yield 3,4-Δ-Pro. A companion pathway with low barriers yielding Δ(1)-pyrroline-5-carboxylate (P5C) is also supported, linking with 5-Hyp through hydration. However, this connection remains unlikely in stressed plants because P5C would be efficiently recycled to Pro (contributing to its accumulation) by P5C reductase, hypothesis coined here as the "Pro-Pro cycle".

Homology modeling of T. cruzi and L. major NADH-dependent fumarate reductases: ligand docking, molecular dynamics validation, and insights on their binding modes.

Leishmania major and Trypanosoma cruzi are the main causes of leishmaniasis and Chagas disease, two endemic parasitosis identified as neglected diseases by the World Health Organization. Fumarate reductase (FR) is a central enzyme in the conversion of fumarate to succinate, an energy releasing path essential for the survival of these protozoans which is also absent in their mammalian hosts. FR can thus be considered as a good candidate for targeting specific inhibition by new drugs designed against L. major and T. cruzi. The lack of tertiary structures available for LmFR and TcFR has limited until now the possibility of performing structure-based drug design. Here we used homology modeling combined with enzyme-cofactor docking to propose tertiary structures for NADH-dependent LmFR and TcFR using an homologous X-ray crystallographic structure of flavine-adenine dinucleotide (FAD) dependent FR from Shewanella frigidimarina (PDB ID: 1QO8) as template. These models were refined and stabilized with/without substrate in the active site using classical molecular dynamics simulations under quasi-physiological conditions. Structural features relevant for understanding the mechanism of action of the enzyme were also analyzed, with special attention to the hydrogen bond network involving the cofactor and water molecules present at the binding sites. A small set of compounds previously synthesized and assayed for their inhibitory capacity against TcFR ([M(mpo)2] metal complexes with M=Pt(II), Pd(II) and V(IV)O and mpo=2-mercaptopyridine N-oxide) and LmFR (licochalcone A) were screened by protein-ligand docking using the NADH-LmFR and NADH-TcFR models here proposed and validated, gaining insight into their binding modes in each enzyme.

Density functional theory characterization and descriptive analysis of cisplatin and related compounds.

Quantum and nonquantum descriptors clearly related to physicochemical features and predictors of the trends to evolve along different stages of a known mechanism of action were determined for a set of square-planar compounds of general formula [M(II)A(1)A(2)L(1)L(2)] (M(II) = Pt(II)/Pd(II); A(i)/L(i) = carrier/labile ligands), structurally related to the anticancer agent Cisplatin. Selected compounds have been sorted and classified by Ward's Cluster Analysis and Principal Components Analysis data-mining techniques using seventeen 1D and two 3D of such theoretical descriptors calculated at the DFT level (PCM-B3LYP/LANL2DZ/6-31G*). A rationale emerging from the study is that whereas most significant differences come from substitution of Cisplatin ligands, cis/trans isomerism, and exchange of M(II) introduce minor alterations in the electronic/geometrical structure. This provides theoretical support to the assay of transplatinum compounds as potential anticancer drugs, a fact already pointed out by empirical evidence. Similarly, the little geometrical/electronic differences triggered by switching M(II) from Pt to Pd enable us to devise a rational path to propose new compounds with expected good anticancer profiles, tuning alterations introduced by simultaneously changing both metal and ligands. Current results serve thus to enlarge the Cleare-Hoeschele guides for Pt(II) square-planar anticancer potential drugs to Pd(II) compounds, both using cis/trans scaffolds.

Structural and Energetic Study of Cisplatin and Derivatives: Comparison of the Performance of Density Funtional Theory Implementations.

In this work, we compare the performance of different DFT implementations, using analytical and numerical basis sets for the expansion of the atomic wave function, in determining structural and energetic parameters of Cisplatin and some biorelevant derivatives. Characterization of the platinum-containing species was achieved at the HF, MP2, and DFT (PBE1PBE, mPW1PW91, B3LYP, B3PW91, and B3P86) levels of theory, using two relativistic effective core potentials to treat the Pt atom (LanL2DZ and SBK), together with analytical Gaussian-type basis sets as implemented in Gaussian03. These results were compared with those obtained with the SIESTA code that employs a pseudopotential derived from the Troullier-Martins procedure for the Pt atom and numerical pseudoatomic orbitals as basis set. All modeled properties were also compared with the experimental values when available or to the best theoretical calculations known to date. On the basis of the results, SIESTA is an excellent alternative to determine structure and energetics of platinum complexes derived from Cisplatin, with less computational efforts. This validates the use of the SIESTA code for this type of chemical systems and thus provides a computationally efficient quantum method (capable to linear scaling at large sizes and available in QM/MM implementations) for exploring larger and more complex chemical models which shall reproduce more faithfully the real chemistry of Cisplatin in physiological conditions.