Do bioactive 8-hydroxyquinolines oxidovanadium(IV) and (V) complexes inhibit the growth of M. smegmatis?

The antiproliferative effects of four series of V^IVO- and V^VO-based compounds containing 8-hydroxyquinoline ligands on the bacterium Mycolicibacterium smegmatis (M. smeg) were investigated. The effects on M. smeg were compared to the antiproliferative effects on the protozoan parasite Trypanosoma cruzi (T. cruzi), the causative agent for Chagas disease. In this study, we investigate the speciation of these compounds under physiological conditions as well as the antiproliferative effects on the bacterium M. smeg. We find that the complexes are more stable the less H₂O is present, and that the stability increases in lipid-like environments. Only one heteroleptic complex and two homoleptic complexes were found to show similar antiproliferative effects on M. smeg as reported for T. cruzi so the responses generally observed by M.smeg. is less than observed by the pathogen. In summary, we find that M. smeg is more sensitive to the detailed structure of the V-complex but overall these complexes are less effective against M. smeg compared to T. cruzi.

Multiple and Variable Binding of Pharmacologically Active Bis(maltolato)oxidovanadium(IV) to Lysozyme

The interaction with proteins of metal-based drugs plays a crucial role in their transport, mechanism, and activity. For an active ML_n complex, where L is the organic carrier, various binding modes (covalent and non-covalent, single or multiple) may occur and several metal moieties (M, ML, ML₂, etc.) may interact with proteins. In this study, we have evaluated the interaction of [V^IVO(malt)₂] (bis(maltolato)oxidovanadium(IV) or BMOV, where malt = maltolato, i.e., the common name for 3-hydroxy-2-methyl-4H-pyran-4-onato) with the model protein hen egg white lysozyme (HEWL) by electrospray ionization mass spectrometry, electron paramagnetic resonance, and X-ray crystallography. The multiple binding of different V-containing isomers and enantiomers to different sites of HEWL is observed. The data indicate both non-covalent binding of cis-[VO(malt)₂(H₂O)] and [VO(malt)(H₂O)₃]⁺ and covalent binding of [VO(H₂O)_3–4]²⁺ and cis-[VO(malt)₂] and other V-containing fragments to the side chains of Glu35, Asp48, Asn65, Asp87, and Asp119 and to the C-terminal carboxylate. Our results suggest that the multiple and variable interactions of potential V^IVOL₂ drugs with proteins can help to better understand their solution chemistry and contribute to define the molecular basis of the mechanism of action of these intriguing molecules.

The interaction of [V^IVO(malt)₂] (BMOV, malt = maltolato) with hen egg white lysozyme (HEWL) reveals the multiple binding of different V-containing isomers and enantiomers to different sites and non-covalent and covalent binding of cis-[VO(malt)₂(H₂O)], [VO(malt)(H₂O)₃]⁺, [VO(H₂O)₃₋₄]²⁺, and cis-[VO(malt)₂] to Glu, Asp, and Asn residues.

Mo(VI) Potential Metallodrugs: Explaining the Transport and Cytotoxicity by Chemical Transformations

The transport and cytotoxicity of molybdenum-based drugs have been explained with the concept of chemical transformation, a very important idea in inorganic medicinal chemistry that is often overlooked in the interpretation of the biological activity of metal-containing systems. Two monomeric, [MoO₂(L¹)(MeOH)] (1) and [MoO₂(L²)(EtOH)] (2), and two mixed-ligand dimeric Mo^VIO₂ species, [{MoO₂(L^1-2)}₂(μ-4,4'-bipy)] (3-4), were synthesized and characterized. The structures of the solid complexes were solved through SC-XRD, while their transformation in water was clarified by UV-vis, ESI-MS, and DFT. In aqueous solution, 1-4 lead to the penta-coordinated [MoO₂(L^1-2)] active species after the release of the solvent molecule (1 and 2) or removal of the 4,4'-bipy bridge (3 and 4). [MoO₂(L^1-2)] are stable in solution and react with neither serum bioligand nor cellular reductants. The binding affinity of 1-4 toward HSA and DNA were evaluated through analytical and computational methods and in both cases a non-covalent interaction is expected. Furthermore, the in vitro cytotoxicity of the complexes was also determined and flow cytometry analysis showed the apoptotic death of the cancer cells. Interestingly, μ-4,4'-bipy bridged complexes 3 and 4 were found to be more active than monomeric 1 and 2, due to the mixture of species generated, that is [MoO₂(L^1-2)] and the cytotoxic 4,4'-bipy released after their dissociation. Since in the cytosol neither the reduction of Mo^VI to Mo^V/IV takes place nor the production of reactive oxygen species (ROS) through Fenton-like reactions of 1-4 with H₂O₂ occurs, the mechanism of cytotoxicity should be attributable to the direct interaction with DNA that happens with a minor-groove binding which results in cell death through an apoptotic mechanism.

Spectroscopic/Computational Characterization and the X-ray Structure of the Adduct of the VIVO-Picolinato Complex with RNase A

The structure, stability, and enzymatic activity of the adduct formed upon the reaction of the V-picolinato (pic) complex [VIVO(pic)2(H2O)], with an octahedral geometry and the water ligand in cis to the V═O group, with the bovine pancreatic ribonuclease (RNase A) were studied. While electrospray ionization-mass spectrometry, circular dichroism, and ultraviolet-visible absorption spectroscopy substantiate the interaction between the metal moiety and RNase A, electron paramagnetic resonance (EPR) allows us to determine that a carboxylate group, stemming from Asp or Glu residues, and imidazole nitrogen from His residues are involved in the V binding at acidic and physiological pH, respectively. Crystallographic data demonstrate that the VIVO(pic)2 moiety coordinates the side chain of Glu111 of RNase A, by substituting the equatorial water molecule at acidic pH. Computational methods confirm that Glu111 is the most affine residue and interacts favorably with the OC-6-23-Δ enantiomer establishing an extended network of hydrogen bonds and van der Waals stabilizations. By increasing the pH around neutrality, with the deprotonation of histidine side chains, the binding of the V complex to His105 and His119 could occur, with that to His105 which should be preferred when compared to that to the catalytically important His119. The binding of the V compound affects the enzymatic activity of RNase A, but it does not alter its overall structure and stability.

Influence of temperature on the equilibria of oxidovanadium(IV) complexes in solution

The equilibria in the solution of three different oxidovanadium(IV) complexes, VO(dhp)₂ (dhp = 1,2-dimethyl-3-hydroxy-4(1H)-pyridinonato), VO(ma)₂ (ma = maltolato) and VO(pic)₂(H₂O) (pic = picolinato), were examined in the temperature range of 120-352 K through a combination of instrumental (EPR spectroscopy) and computational techniques (DFT methods). The results revealed that a general equilibrium exists: VOL₂ + H₂O ⇄ cis-VOL₂(H₂O) ⇄ trans-VOL₂(H₂O), where cis and trans refer to the relative position of H₂O and the oxido ligand. The equilibrium is more or less shifted to the right depending on the ligand, the temperature, the ionic strength and the coordinating properties of the solvent. With VO(dhp)₂, only the square pyramidal species exists at 298 K in aqueous solution, while at 120 K the cis- and trans-VO(dhp)₂(H₂O) species are also present. The complex of maltol exists almost exclusively in the form cis-VO(ma)₂(H₂O) in aqueous solution at 298 K, while the trans species can be revealed only at higher temperatures, where the EPR linewidth significantly decreases. The equilibria involving 1-methylimidazole (MeIm), a model for the side chain His coordination, are also influenced by temperature, with its coordination being favored by decreasing the temperature. The implications of these results in the study of the (vanadium complex)-protein systems are discussed and the interaction with myoglobin (Mb) is examined as a representative example.

Integration of Mass Spectrometry Data for Structural Biology

Mass spectrometry (MS) is increasingly being used to probe the structure and dynamics of proteins and the complexes they form with other macromolecules. There are now several specialized MS methods, each with unique sample preparation, data acquisition, and data processing protocols. Collectively, these methods are referred to as structural MS and include cross-linking, hydrogen-deuterium exchange, hydroxyl radical footprinting, native, ion mobility, and top-down MS. Each of these provides a unique type of structural information, ranging from composition and stoichiometry through to residue level proximity and solvent accessibility. Structural MS has proved particularly beneficial in studying protein classes for which analysis by classic structural biology techniques proves challenging such as glycosylated or intrinsically disordered proteins. To capture the structural details for a particular system, especially larger multiprotein complexes, more than one structural MS method with other structural and biophysical techniques is often required. Key to integrating these diverse data are computational strategies and software solutions to facilitate this process. We provide a background to the structural MS methods and briefly summarize other structural methods and how these are combined with MS. We then describe current state of the art approaches for the integration of structural MS data for structural biology. We quantify how often these methods are used together and provide examples where such combinations have been fruitful. To illustrate the power of integrative approaches, we discuss progress in solving the structures of the proteasome and the nuclear pore complex. We also discuss how information from structural MS, particularly pertaining to protein dynamics, is not currently utilized in integrative workflows and how such information can provide a more accurate picture of the systems studied. We conclude by discussing new developments in the MS and computational fields that will further enable in-cell structural studies.

Insights into metalloproteins and metallodrugs from electron paramagnetic resonance spectroscopy

Metal ions play an important role in diverse biological processes, and much of the basic knowledge derived from studying native bioinorganic systems are applied in the synthesis of new molecules with the aim of diagnosing and treating diseases. At first glance, metalloproteins and metallodrugs are very different systems, but metal ion coordination, redox chemistry and substrate binding play essential roles in advancing both of these research fields. In this article, we discuss recent metalloprotein and metallodrug studies where electron paramagnetic resonance spectroscopy served as a major tool to gain a better understanding of metal-based structures and their function.

Rationalizing the Decavanadate(V) and Oxidovanadium(IV) Binding to G-Actin and the Competition with Decaniobate(V) and ATP

The experimental data collected over the past 15 years on the interaction of decavanadate(V) (V₁₀O₂₈^6-; V₁₀), a polyoxometalate (POM) with promising anticancer and antibacterial action, with G-actin, were rationalized by using several computational approaches (docking, density functional theory (DFT), and molecular dynamics (MD)). Moreover, a comparison with the isostructural and more stable decaniobate(V) (Nb₁₀O₂₈^6-; Nb₁₀) was carried out. Four binding sites were identified, named α, β, γ, and δ, the site α being the catalytic nucleotide site located in the cleft of the enzyme at the interface of the subdomains II and IV. It was observed that the site α is preferred by V₁₀, whereas Nb₁₀ is more stable at the site β; this indicates that, differently from other proteins, G-actin could contemporaneously bind the two POMs, whose action would be synergistic. Both decavanadate and decaniobate induce conformational rearrangements in G-actin, larger for V₁₀ than Nb₁₀. Moreover, the binding mode of oxidovanadium(IV) ion, V^IVO²⁺, formed upon the reduction of decavanadate(V) by the -SH groups of accessible cysteine residues, is also found in the catalytic site α with (His161, Asp154) coordination; this adduct overlaps significantly with the region where ATP is bound, accounting for the competition between V₁₀ and its reduction product V^IVO²⁺ with ATP, as previously observed by EPR spectroscopy. Finally, the competition with ATP was rationalized: since decavanadate prefers the nucleotide site α, Ca²⁺-ATP displaces V₁₀ from this site, while the competition is less important for Nb₁₀ because this POM shows a higher affinity for β than for site α. A relevant consequence of this paper is that other metallodrug-protein systems, in the absence or presence of eventual inhibitors and/or competition with molecules of the organism, could be studied with the same approach, suggesting important elements for an explanation of the biological data and a rational drug design.

Integrated experimental/computational approaches to characterize the systems formed by vanadium with proteins and enzymes

An integrated instrumental/computational approach to characterize metallodrug–protein adducts at the molecular level is reviewed. A series of applications are described, focusing on potential vanadium drugs with a generalization to other metals.

Photoactive PtII and PdII complexes of N,N-diethyl-N′-3,4,5-trimethoxybenzoylthiourea: synthesis, crystal structures, DFT and cytotoxicity studies

New trans-Pd^II/Pt^II acylthiourea complexes prepared by photo-induced isomerism are found to be more cytotoxic than cis-Pd^II/Pt^II against human prostate cancer and normal embryonic kidney cell lines.

Covalent and non-covalent binding in vanadium–protein adducts

An integrated method, generalizable to any metals and proteins, based on ESI-MS, EPR and molecular modelling was applied to study the covalent and non-covalent binding of the potential drug [V^IVO(nalidixato)₂(H₂O)] to lysozyme and cytochrome c.

BioMetAll: Identifying Metal-Binding Sites in Proteins from Backbone Preorganization

With a large amount of research dedicated to decoding how metallic species bind to proteins, in silico methods are interesting allies for experimental procedures. To date, computational predictors mostly work by identifying the best possible sequence or structural match of the target protein with metal-binding templates. These approaches are fundamentally focused on the first coordination sphere of the metal. Here, we present the BioMetAll predictor that is based on a different postulate: the formation of a potential metal-binding site is related to the geometric organization of the protein backbone. We first report the set of convenient geometric descriptors of the backbone needed for the algorithm and their parameterization from a statistical analysis. Then, the successful benchmark of BioMetAll on a set of more than 90 metal-binding X-ray structures is presented. Because BioMetAll allows structural predictions regardless of the exact geometry of the side chains, it appears extremely valuable for systems whose structures (either experimental or theoretical) are not optimal for metal-binding sites. We report here its application on three different challenging cases: (i) the modulation of metal-binding sites during conformational transition in human serum albumin, (ii) the identification of possible routes of metal migration in hemocyanins, and (iii) the prediction of mutations to generate convenient metal-binding sites for de novo biocatalysts. This study shows that BioMetAll offers a versatile platform for numerous fields of research at the interface between inorganic chemistry and biology and allows to highlight the role of the preorganization of the protein backbone as a marker for metal binding. BioMetAll is an open-source application available at https://github.com/insilichem/biometall.

New VIV, VIVO, VVO, and VVO2 Systems: Exploring their Interconversion in Solution, Protein Interactions, and Cytotoxicity

The synthesis and characterization of one oxidoethoxidovanadium(V) [V^VO(L¹)(OEt)] (1) and two nonoxidovanadium(IV) complexes, [V^IV(L^2-3)₂] (2 and 3), with aroylhydrazone ligands incorporating naphthalene moieties, are reported. The synthesized oxido and nonoxido vanadium complexes are characterized by various physicochemical techniques, and their molecular structures are solved by single crystal X-ray diffraction (SC-XRD). This revealed that in 1 the geometry around the vanadium atom corresponds to a distorted square pyramid, with a O₄N coordination sphere, whereas that of the two nonoxido V^IV complexes 2 and 3 corresponds to a distorted trigonal prismatic arrangement with a O₄N₂ coordination sphere around each "bare" vanadium center. In aqueous solution, the V^VO moiety of 1 undergoes a change to V^VO₂ species_, yielding [V^VO₂(L¹)]^- (1'), while the nonoxido V^IV-compounds 2 and 3 are partly converted into their corresponding V^IVO complexes, [V^IVO(L^2-3)(H₂O)] (2' and 3'). Interaction of these V^VO₂, V^IVO, and V^IV systems with two model proteins, ubiquitin (Ub) and lysozyme (Lyz), is investigated through docking approaches, which suggest the potential binding sites: the interaction is covalent for species 2' and 3', with the binding to Glu16, Glu18, and Asp21 for Ub, and His15 for Lyz, and it is noncovalent for species 1', 2, and 3, with the surface residues of the proteins. The ligand precursors and complexes are also evaluated for their in vitro antiproliferative activity against ovarian (A2780) and prostate (PC3) human cancer cells and in normal fibroblasts (V79) to check the selectivity of the compounds for cancer cells.

ESI-MS Study of the Interaction of Potential Oxidovanadium(IV) Drugs and Amavadin with Model Proteins

In this study, the binding to lysozyme (Lyz) of four important VIV compounds with antidiabetic and/or anticancer activity, [VIVO(pic)2(H2O)], [VIVO(ma)2], [VIVO(dhp)2], and [VIVO(acac)2], where pic-, ma-, dhp-, and acac- are picolinate, maltolate, 1,2-dimethyl-3-hydroxy-4(1H)-pyridinonate, and acetylacetonate anions, and of the vanadium-containing natural product amavadin ([VIV(hidpa)2]2-, with hidpa3- N-hydroxyimino-2,2'-diisopropionate) was investigated by ElectroSpray Ionization-Mass Spectrometry (ESI-MS). Moreover, the interaction of [VIVO(pic)2(H2O)], chosen as a representative VIVO2+ complex, was examined with two additional proteins, myoglobin (Mb) and ubiquitin (Ub), to compare the data. The examined vanadium concentration was in the range 15-150 μM, i.e., very close to that found under physiological conditions. With pic-, dhp-, and hidpa3-, the formation of adducts n[VIVOL2]-Lyz or n[VIVL2]-Lyz is favored, while with ma- and acac- the species n[VIVOL]-Lyz are detected, with n dependent on the experimental VIV/protein ratio. The behavior of the systems with [VIVO(pic)2(H2O)] and Mb or Ub is very similar to that of Lyz. The results suggested that under physiological conditions, the moiety cis-VIVOL2 (L = pic-, dhp-) is bound by only one accessible side-chain protein residue that can be Asp, Glu, or His, while VIVOL+ (L = ma-, acac-) can interact with the two equatorial and axial sites. If the VIV complex is thermodynamically stable and does not have available coordination positions, such as amavadin, the protein cannot interact with it through the formation of coordination bonds and, in such cases, noncovalent interactions are predicted. The formation of the adducts is dependent on the thermodynamic stability and geometry in aqueous solution of the VIVO2+ complex and affects the transport, uptake, and mechanism of action of potential V drugs.

Molecular Modeling for Artificial Metalloenzyme Design and Optimization

Artificial metalloenzymes (ArMs) are obtained by inserting homogeneous catalysts into biological scaffolds and are among the most promising strategies in the quest for new-to-nature biocatalysts. The quality of their design strongly depends on how three partners interact: the biological host, the "artificial cofactor," and the substrate. However, structural characterization of functional artificial metalloenzymes by X-ray or NMR is often partial, elusive, or absent. How the cofactor binds to the protein, how the receptor reorganizes upon the binding of the cofactor and the substrate, and which are the binding mode(s) of the substrate for the reaction to proceed are key questions that are frequently unresolved yet crucial for ArM design. Such questions may eventually be solved by molecular modeling but require a step change beyond the current state-of-the-art methodologies.Here, we summarize our efforts in the study of ArMs, presenting both the development of computational strategies and their application. We first focus on our integrative computational framework that incorporates a variety of methods such as protein-ligand docking, classical molecular dynamics (MD), and pure quantum mechanical (QM) methods, which, when properly combined, are able to depict questions that range from host-cofactor binding predictions to simulations of entire catalytic mechanisms. We also pay particular attention to the protein-ligand docking strategies that we have developed to accurately predict the binding of transition metal-containing molecules to proteins. While this aspect is fundamental to many bioinorganic fields beyond ArMs, it has been disregarded from the molecular modeling landscape until very recently.Next we describe how to apply this computational framework to particular ArMs including systems previously characterized experimentally as well as others where computation served to guide the design. We start with the prediction of the interactions between homogeneous catalysts and biological hosts. Protein-ligand docking is pivotal at that stage, but it needs to be combined with QM/MM or MD approaches when the binding of the cofactor implies significant conformational changes of the protein or involve changes of the electronic state of the metal.Then, we summarize molecular modeling studies aimed at identifying cofactor-substrate arrangements inside the ArM active pocket that are consistent with its reactivity. These calculations stand on "Theozyme"-like dockings, MD-refined or not, which provide molecular rationale of the catalytic profiles of the artificial systems.In the third section, we present case studies to decode the entire catalytic mechanism of two ArMs: (1) an iridium based asymmetric transfer hydrogenase obtained by insertion of Noyori's catalyst into streptavidin and (2) a metallohydrolase achieved by including a receptor. Transition states, second coordination sphere effects, as well as motions of the cofactors are identified as drivers of the enantiomeric profiles.Finally, we report computer-aided designs of ArMs to guide experiments toward chemical and mutational changes that improve their activity and/or enantioselective profiles and expand toward future directions.

Chemistry of mixed-ligand oxidovanadium(IV) complexes of aroylhydrazones incorporating quinoline derivatives: Study of solution behavior, theoretical evaluation and protein/DNA interaction

A series of eight hexacoordinated mixed-ligand oxidovanadium(IV) complexes [VO(L^x)(L^N-N)] (1-8), where L^x = L¹ - L⁴ are four differently substituted ONO donor aroylhydrazone ligands and L^N-N are N,N-donor bases like 2,2'-bipyridine (bipy) (1, 3, 5 and 7) and 1,10-phenanthroline (phen) (2, 4, 6 and 8), have been reported. All synthesized complexes have been characterized by various physicochemical techniques and molecular structures of 1 and 6 were determined by X-ray crystallography. With a view to evaluate the biological activity of the V^IVO species, the behavior of the systems V^IVO²⁺/L^x, V^IVO²⁺/L^x/bipy and V^IVO²⁺/L^x/phen was studied as a function of pH in a mixture of H₂O/DMSO 50/50 (v/v). DFT calculations allowed finding out the relative stability of the tautomeric forms of the ligands, and predicting the structure of vanadium complexes and their EPR parameters. To study their interaction with proteins, firstly the ternary systems V^IVO²⁺/L^1,2 with 1-methylimidazole, which is a good model for histidine binding, were examined. Subsequently the interaction of the complexes with lysozyme (Lyz), cytochrome c (Cyt) and bovine serum albumin (BSA) was studied. The results indicate that the complexes showed moderate binding affinity towards BSA, while no interaction takes place with lysozyme and cytochrome c. This could be explained with the higher number of accessible coordinating and polar residues for BSA than for Lyz and Cyt. Further, the complexes were also evaluated for their DNA binding propensity through UV-vis absorption titration and fluorescence spectral studies. These results were consistent with BSA binding affinity and showed moderate binding affinity towards CT-DNA.

Effect of secondary interactions, steric hindrance and electric charge on the interaction of VIVO species with proteins

The effect of secondary interactions (hydrogen bonds and van der Waals contacts), steric hindrance and electric charge, on the binding of V^IV complexes formed by pipemidic and 8-hydroxyquinoline-5-sulphonic acids with ubiquitin and lysozyme is studied.