In silico evolution of substrate selectivity: comparison of organometallic ruthenium complexes with the anticancer drug cisplatin

Sulfur Analogs of the Core Formose Cycle: A Free Energy Map

Using computational methods, we examine if the presence of H2S can tame the unruly formose reaction by generating a free energy map of the reaction thermodynamics and kinetics of sulfur analogs within the core cycle. With mercaptoaldehyde as the linchpin C2 species, and feeding the cycle with CH2O, selected aldol additions and enolizations are kinetically more favorable. Thione formation is thermodynamically less favored compared to aldehydes and ketones, but all these species can be connected by enolization reactions. In some sulfur analogs, the retroaldol transformation of a C4 species back into linchpin species is thermodynamically favorable, and we have found one route incorporating where incorporating sulfur selects for a specific pathway over others. However, as CH2O diminishes, the aldol addition of larger species is less favorable for the sulfur analogs. Our results also suggest that competing Cannizzaro side reactions are kinetically less favored and thermodynamically disfavored when H2S is abundant.

Computation of rate coefficients in solutions based on transition state theory combined with a heuristically corrected polarizable continuum model: intermolecular Diels-Alder reactions as case studies

Transition state theory (TST) based on activation parameters computed using quantum mechanics calculations combined with the polarizable continuum model (QM/PCM) is a fundamental tool for investigating reaction rates in the liquid phase. In conventional QM/PCM methods, thermodynamic data and partition functions for a solute are often derived from a quasi-ideal gas treatment (IGT) widely implemented in commercially available computation packages. This approach tends to overestimate entropy because calculations of thermodynamic parameters in the liquid phase ignore hindered translational and rotational modes in real solutions. The present work formulated partition functions for more realistic solutes hindered by surrounding solvent molecules in conjunction with the basic QM/PCM concept. In addition, a configuration partition function for solute molecules at a standard concentration of 1 mol dm-3 was incorporated using a simple lattice model. The canonical partition function and thermodynamic functions were derived based on statistical thermodynamics for localized systems. Expressions for rate coefficients within TST were also derived with a consistent formalism based on the standard state selected in partition function calculations. The performance of the proposed method was assessed by predicting rate coefficients for three different Diels-Alder reactions and comparing these with experimental results. QM/PCM calculations at the G4//ωB97X-D/6-311++G(d,p)/IEF-PCM level of theory with corrections for the dispersion and repulsion energies were performed to obtain the electronic structures of stationary points on potential energy surfaces as a means of finding activation enthalpy, entropy and Gibbs energy values based on revised partition functions as well as predicting rate coefficients. The activation Gibbs energies obtained from our proposed method were lower than those obtained from the IGT method due to reasonable entropy computations. The proposed method overestimated the rate coefficients by one to two orders of magnitude compared to the experimental values, whereas the IGT method underestimated them by the same amount. This discrepancy arises because the proposed method calculates the partition function from the viewpoint of a localized system, whereas the IGT method calculates it from the viewpoint of a non-localized system. Given that actual liquids exist in a state between non-localized and localized systems, it is essential to formulate the partition function in a way that more accurately represents the liquid state.

Exploring the Core Formose Cycle: Catalysis and Competition

The core autocatalytic cycle of the formose reaction may be enhanced or eroded by the presence of simple molecules at life's origin. Utilizing quantum chemistry, we calculate the thermodynamics and kinetics of reactions both within the core cycle and those that deplete the reactants and intermediates, such as the Cannizzaro reaction. We find that via disproportionation of aldehydes into carboxylic acids and alcohols, the Cannizzaro reaction furnishes simple catalysts for a variety of reactions. We also find that ammonia can catalyze both in-cycle and Cannizzaro reactions while hydrogen sulfide does not; both, however, play a role in sequestering reactants and intermediates in the web of potential reactions.

Quantum Chemical Investigation on Hydrolysis of Orally Active Organometallic Ruthenium(II) and Osmium(II) Anticancer Drugs and Their Interaction with Histidine

Influence of the metal center on hydrolysis of organometallic anticancer complexes containing an N-phenyl-2-pyridinecarbothioamide (PCA) ligand, [M(η⁶-p-cymene)(N-phenyl-2-pyridinecarbothioamide)Cl]⁺ (M = Ru^II, 1A, and Os^II, 2A), as well as their N-fluorophenyl derivatives [M(η⁶-p-cymene)(N-fluorophenyl-2-pyridinecarbothioamide)Cl]⁺ (M = Ru^II, 1B, and Os^II, 2B) have been investigated using the DFT method in aqueous medium. The activation energy barriers for the hydrolysis of 1A (21.5 kcal/mol) and 1B (20.7 kcal/mol) are found to be significantly lower than those of their corresponding osmium analogs 2A (28.6 kcal/mol) and 2B (27.5 kcal/mol). DFT evaluated results reveal the inertness of Os(II)-PCA complex toward the hydrolysis that rationalizes the experimental observations. However, the incorporation of fluoride substituent slightly decreases the activation energy for the hydrolysis of Ru(II)- and Os(II)-PCA. In addition, the interaction of hydrolyzed Ru(II)-PCAs (1AH and 1BH) and Os(II)-PCAs (2AH and 2BH) complexes with the histidine (Hist) have also been investigated. The aquated 1BH and 2BH show an enhanced propensity toward the interaction with histidine, and their activation Gibbs free energies are calculated to be 15.9 and 18.9 kcal/mol, respectively. ONIOM (QM/MM) study of the resulting aquated complexes inside histone protein shows the maximum stability of the 2BH complex having a binding energy of -43.6 kcal/mol.

Preliminary Free Energy Map of Prebiotic Compounds Formed from CO2, H2 and H2S

What kinds of CHOS compounds might be formed in a prebiotic milieu by reducing CO2 in the presence of H2 and H2S? How might the presence of sulfur influence the chemical composition of the mixture? We explore these questions by using first-principles quantum chemistry to calculate the free energies of CHOS compounds in aqueous solution, by first generating a thermodynamic map of one- and two-carbon species. We find that while thiols are thermodynamically favored, thioesters, thioacids, and thiones are less favorable than their non-sulfur counterparts. We then focus on the key role played by mercaptoacetaldehyde in sulfur analogs of the autocatalytic formose reaction, whereby the thiol group introduces asymmetry and potential thermodynamic selectivity of some compounds over others.

Thermodynamics of Potential CHO Metabolites in a Reducing Environment

How did metabolism arise and evolve? What chemical compounds might be suitable to support and sustain a proto-metabolism before the advent of more complex co-factors? We explore these questions by using first-principles quantum chemistry to calculate the free energies of CHO compounds in aqueous solution, allowing us to probe the thermodynamics of core extant cycles and their closely related chemical cousins. By framing our analysis in terms of the simplest feasible cycle and its permutations, we analyze potentially favorable thermodynamic cycles for CO₂ fixation with H₂ as a reductant. We find that paying attention to redox states illuminates which reactions are endergonic or exergonic. Our results highlight the role of acetate in proto-metabolic cycles, and its connection to other prebiotic molecules such as glyoxalate, glycolaldehyde, and glycolic acid.

Using general computational chemistry strategy to unravel the reactivity of emerging pollutants: An example of sulfonamide chlorination

Increasing number of emerging pollutants in environments requires an effective approach which can facilitate the prediction of reactivity and provide insights into the reaction mechanisms. Computational chemistry is exactly the tool to fulfill this demand with its good performance in theoretical investigation of chemical reactions at molecular level. In this study, chlorination of sulfonamide antibiotics is used as an illustration to present a systematic strategy demonstrating how computational chemistry can be applied to investigate the reaction behavior of emerging pollutants. Sulfonamides is a class of micropollutants that contain the common structure of 4-aminobenzenesulfonmaide while differ in their heterocycles. Based on the calculated conceptual DFT indices, the reactive sites of sulfonamide are successfully predicted, which locate on their common structure of 4-aminobenzenesulfonmaide. Therefore, all sulfonamides follow the similar reaction pathway. Product identification by LTQ-Orbitrap MS further verifies the in silico prediction. Three critical pathways are discovered, i.e., S-N bond cleavage, Cl-substitution onto aniline-N, and the following rearrangement to lose -SO₂- group, among which Cl-substitution is the key step due to its lowest free energy barrier. Heterocycles impact the reaction rate by affecting the electronic density of aniline group. In general, the more electron-donating the heterocycle is, the more readily sulfonamides to be chlorinated.

A detailed quantum chemical investigation on the hydrolysis mechanism of osmium(iii) anticancer drug, (ImH)[trans-OsCl4(DMSO)(Im)] (Os-NAMI-A; Im = imidazole)

Detailed hydrolysis mechanism of osmium(iii) anticancer drug, (ImH)[trans-OsCl₄(DMSO)(Im)] (Os-NAMI-A; Im = imidazole, DMSO = dimethyl sulfoxide) has been investigated using density functional theory (DFT) in combination with CPCM solvation model.

Early Steps of Glycolonitrile Oligomerization: A Free-Energy Map

Building on our previous free-energy map (J. Phys. Chem. A 2018, 122, 6769-6779) examining the reactions of monomeric glycolonitrile, we explore the formation of its dimers and trimers in aqueous solution under neutral conditions. While 5-membered rings are kinetically favored, open-chain oligomers with ester or amide linkages are thermodynamically favored. Accessing the 5-membered rings also provides a potential route to glyoxal that bypasses preforming glycolamide, the thermodynamic sink for monomers. However, finding a kinetically accessible route to glycine starting from glycolonitrile in the absence of added ammonia at room temperature proved challenging; the best case involved an intramolecular nucleophilic substitution reaction in a dimer containing neighboring imine and amide groups. Our free-energy map also examines routes to experimentally proposed moieties, explaining why some are observed in low yield or not at all.

Activation of the Ir-N(pyridine) Bond in Half-Sandwich Tethered Iridium(III) Complexes

We present four new organometallic half-sandwich iridium(III) complexes of formula [Ir(η⁵:κ¹-C₅Me₄CH₂py)(N,N)](PF₆)₂, bearing a N,N-chelating ligand [ethylenediamine (en), 1; 1,3-diaminopropane (dap), 2; 2,2'-bipyridine (bipy), 3; 1,10-phenanthroline (phen), 4]; and a derivatized cyclopentadienyl ligand, C₅Me₄CH₂C₅H₄N, which forms an additional five-membered chelate. The latter is hemilabile, and the Ir-N(py) bond can be reversibly cleaved by various stimuli. The four complexes are unreactive toward hydrolysis at pH 7. Interestingly, 1 and 2 react with hydrochloric acid and formate, and speciation between closed and open tether complexes can be followed by ¹H NMR spectroscopy. Complex 1 binds to nucleobase guanine (9-ethylguanine, 9-EtG), yet interaction to calf-thymus DNA was not observed. New X-ray structures of closed tether complexes 1-4 and open tether complexes [Ir(η⁵-C₅Me₄CH₂pyH)(en)Cl](PF₆)₂ (1·HCl) and [Ir(η⁵-C₅Me₄CH₂py)(en)H]PF₆ (1·hyd) have been determined. Hydride capture is efficient for 1 and 2. The kinetics of Ir-H bond formation and hydride transfer in a model organic molecule have been investigated, revealing a strong dependence on the temperature. Coincubation of complex 1 with nontoxic concentrations of sodium formate decreases the IC₅₀ value in MCF7 breast cancer cells, indicating the possibility of intracellular activation of the Ir-N(py) tether bond to generate cytotoxic activity via iridium-mediated transfer hydrogenation.

Revisiting the Active Sites at the MoS2/H2O Interface via Grand-Canonical DFT: The Role of Water Dissociation

MoS₂ is a promising low-cost catalyst for the hydrogen evolution reaction (HER). However, the nature of the active sites remains a subject of debate. By taking the electrochemcal potential explicitly into account using grand-canonical density functional theory (DFT) in combination with the linearized Poisson-Boltzmann equation, we herein revisit the active sites of 2H-MoS₂. In addition to the well-known catalytically active edge sites, also specific point defects on the otherwise inert basal plane provide highly active sites for HER. Given that HER takes place in water, we also assess the reactivity of these active sites with respect to H₂O. The thermodynamics of proton reduction as a function of the electrochemical potential reveals that four edge sites and three basal plane defects feature thermodynamic overpotentials below 0.2 V. In contrast to current proposals, many of these active sites involve adsorbed OH. The results demonstrate that even though H₂O and OH block "active" sites, HER can also occur on these "blocked" sites, reducing protons on surface OH/H₂O entities. As a consequence, our results revise the active sites, highlighting the so far overlooked need to take the liquid component (H₂O) of the functional interface into account when considering the stability and activity of the various active sites.

Imidazole and Benzimidazole Modified Half-Sandwich IridiumIII N-Heterocyclic Carbene Complexes: Synthesis, Anticancer Application, and Organelle Targeting

Herein, we report the synthesis, characterization and anticancer activity of a series of half-sandwich iridium^III imidazole and benzimidazole N-heterocyclic carbene (NHC) anticancer complexes, and the general formula of which can be expressed as [(η⁵-Cp^x)Ir(C^∧N)Cl]Cl (Cp^x: pentamethylcyclopentadienyl (Cp^*) or biphenyl derivatives (Cp^xbiph); C^∧N: imidazole and benzimidazole NHC chelating ligands). Compared with cis-platin, these complexes showed interesting antitumor activity against A549 cells. Complexes could bind to bovine serum albumin (BSA) by means of static quenching mode, catalyze the oxidation of nicotinamide adenine dinucleotide (NADH) and increase the levels of reactive oxygen species (ROS). Meanwhile, these complexes could arrest the cell cycles of A549 cells and influence the mitochondrial membrane potential significantly. Due to the inherent luminescence property, laser confocal test show that complexes could enter cells followed an energy-dependent mechanism and effectively accumulate in lysosome (the value of Pearson's co-localization coefficient is 0.70 after 1 h), further destroy lysosome integrity and induce apoptosis.

Role of Acid in the Co-oligomerization of Formaldehyde and Pyrrole

Building on previous work (J. Phys. Chem. A 2017, 121, 8154-8166) under neutral conditions, we examined the co-oligomerization of CH₂O and pyrrole to form porphryinogen under acidic conditions using density functional theory (B3LYP//6-311G**). Thermodynamically, we found that azafulvene intermediates were significantly stabilized under highly acidic conditions. Kinetically, energy barriers were lowered for C-C bond formation, discriminating in favor of reactions that lead to porphyrinogen. However, it was challenging to satisfactorily combine our thermodynamic and kinetic profiles into a unified free-energy profile because of difficulties in optimizing transition states of cationic species involving proton hops. Instead, we used neutral carboxylic acids as a proxy to study how energy barriers changed. By combining data from both neutral and acidic conditions, we estimate a free-energy profile for the initial steps of oligomerization under milder acidic conditions more relevant to prebiotic chemistry.

Theoretical analysis of glyoxal condensation with ammonia in aqueous solution

The reactions of glyoxal with ammonia, ammonium salts, and amines cause the formation of the secondary organic aerosol (SOA) components (imidazole and its derivatives) in the atmosphere. The interaction of glyoxal and ammonia in aqueous solution is a primary reaction for these processes, and the explanation of its mechanism will allow developing the methods to control the formation of the SOA components. A detailed mechanism for the formation of key intermediates, namely, ethanediimine, diaminoethanediol, and aminoethanetriol, required for the imidazole ring cyclization, is proposed, and its potential energy surface (PES) has been constructed. This mechanism includes the experimentally identified intermediate compounds and takes into account the conformational and hydration equilibria of glyoxal. The schemes are proposed for further conversion of the key intermediates to the products of condensation between glyoxal and ammonia in the aqueous solution, C-N cyclic oligomers, that were identified. The products are shown to correspond to low positions on the PES in terms of Gibbs free energy, from -30.8 to -68.3 kcal mol-1, which confirms the high probability of their formation. The preferable thermodynamic pathway for formation of the imidazole products does not comprise the conversion of the diimine intermediate with the participation of the proton, but rather the interaction of either the diaminoalcohol with glyoxal monohydrate or two monoamine derivatives between themeselves (aminoethantriol and aminohydroxyacetaldehyde).

Exploring Free Energy Profiles of Uracil and Cytosine Reactions with Formaldehyde

Simple polymers can be potentially formed by the co-oligomerization of pyrimidine nucleobases, uracil and cytosine, with the small molecule formaldehyde. Using density functional calculations, we have constructed a free energy map outlining the thermodynamics and kinetics for (1) the addition of formaldehyde to uracil and cytosine to form hydroxymethylated uracil (HMU) and hydroxymethylated cytosine (HMC), (2) the deamination of cytosine and HMC to uracil and HMU, respectively, and (3) the initial oligomerization of 5-HMU. For the initial formation of monomeric HMU, addition of formaldehyde to the C5 and C6 positions is thermodynamically favored over N1 and N3, but faces higher kinetic barriers, and explains why 5-HMU is the main product observed experimentally. Oligomerization of 5-HMU is thermodynamically favorable although decreasingly so at the tetramer stage. In addition, decreasing concentrations of initial monomer shifts the equilibrium disfavoring oligomer formation.

A simple heuristic approach to estimate the thermochemistry of condensed-phase molecules based on the polarizable continuum model

A simple model based on a quantum chemical approach with polarizable continuum models (PCMs) to provide reasonable translational and rotational entropies for liquid phase molecules was developed.

Reactions of Glycolonitrile with Ammonia and Water: A Free Energy Map

Glycolonitrile, the product of combining CH₂O and HCN, is an intermediate in the Strecker reaction leading to the synthesis of the amino acid glycine. However, besides glycine, a plethora of other compounds are also generated when CH₂O and HCN react in the presence of ammonia and water. As a starting point to analyze the possible components of this complex mixture, we have employed density functional theory to construct a free energy map of all two-carbon (C2) species that may be present when glycolonitrile participates in addition or elimination reactions with ammonia and water. By identifying thermodynamic sinks and kinetic barriers, we find that the myriad C2 species can be grouped into three broad regions across the free energy landscape. This allows us to trace possible routes to glycine and other molecules of interest in the reaction mixture. The present map also extends our previous work on one-carbon (C1) species. We had previously found one issue with our computational protocol in the C1 map; however, our present C2 map provides a larger data set that supports using an empirical correction to our original protocol for imidic acid to amide transformations, without increasing the computational cost, while retaining the original protocol for other classes of reactions.

Influence of the substituent on the phosphine ligand in novel rhenium(i) aldehydes. Synthesis, computational studies and first insights into the antiproliferative activity

Cyrhetrenyl phosphine derivatives were synthesized and evaluated as potential anticancer agents. Electrochemical and computational studies were carried out.

Electronic effects on reactivity and anticancer activity by half-sandwich N,N-chelated iridium(iii) complexes

This work demonstrated how the chemical reactivity and anticancer activity as well as the selectivity of these half-sandwich N,N-chelated iridium(iii) complexes can be controlled and fine-tuned by the modification of the ligand electronic perturbations.

Nucleophilic Addition of Ketones To Acetylenes and Allenes: A Quantum-Chemical Insight

A CBS-Q//B3 based study has been carried out to elucidate the mechanism of the KOH/DMSO superbase catalyzed ketones nucleophilic addition to alkyl propargyl and alkyl allenyl ethers yielding, along with (Z)-monoadducts, up to 26% of unexpected (E)-diadducts. The impact of different substrates (alkynes versus allenes) on the reaction mechanism has been discussed in detail. Along with the model reaction of acetone addition to propyne and allene, the addition of acetone and acetophenone to methyl propargyl and methyl allenyl ethers is considered. The limiting reaction stage of the starting ketone carbanion addition to propargyl and allenyl systems occurs with activation energies typical for vinylation of ketones. In contrast, the addition of intermediate α-carbanions to the terminal position of methyl allenyl ether is associated with unusually low activation barriers. The results obtained explain the composition of the reaction products and indicate the participation of mainly the allene form in the reaction.

Half-sandwich iridium N-heterocyclic carbene anticancer complexes

Half-sandwich iridium N-heterocyclic carbene complexes display potent anticancer activities and are attractive for development as new anticancer agents.

Theoretical estimation of kinetic parameters for nucleophilic substitution reactions in solution: an application of a solution translational entropy model

The kinetic parameters, such as activation entropy, activation enthalpy, activation free-energy, and reaction rate constant, for a series of nucleophilic substitution (SN) reactions in solution, are investigated using both a solution-phase translational entropy model and an ideal gas-phase translational entropy model. The results obtained from the solution translational entropy model are in excellent agreement with the experimental values, while the overestimation of activation free-energy from the ideal gas-phase translational entropy model is as large as 6.9 kcal mol(-1). For some of the reactions studied, such as and in methanol, and and in aqueous solution, the explicit + implicit model, namely, a cluster-continuum type model, should be employed to account for the strong solvent-solute interactions. In addition, the explicit + implicit models have also been applied to the DMSO-H2O mixtures, which would open up a door to investigate the reactions in a mixed solvent using density functional theory (DFT) methods.

Interactions of the "piano-stool" [ruthenium(II)(η(6) -arene)(quinolone)Cl](+) complexes with water; DFT computational study

Full optimizations of stationary points along the reaction coordinate for the hydration of several quinolone Ru(II) half-sandwich complexes were performed in water environment using the B3PW91/6-31+G(d)/PCM/UAKS method. The role of diffuse functions (especially on oxygen) was found crucial for correct geometries along the reaction coordinate. Single-point (SP) calculations were performed at the B3LYP/6-311++G(2df,2pd)/DPCM/saled-UAKS level. In the first part, two possible reaction mechanisms-associative and dissociative were compared. It was found that the dissociative mechanism of the hydration process is kinetically slightly preferred. Another important conclusion concerns the reaction channels. It was found that substitution of chloride ligand (abbreviated in the text as dechlorination reaction) represents energetically and kinetically the most feasible pathway. In the second part the same hydration reaction was explored for reactivity comparison of the Ru(II)-complexes with several derivatives of nalidixic acid: cinoxacin, ofloxacin, and (thio)nalidixic acid. The hydration process is about four orders of magnitude faster in a basic solution compared to neutral/acidic environment with cinoxacin and nalidixic acid as the most reactive complexes in the former and latter environments, respectively. The explored hydration reaction is in all cases endergonic; nevertheless the endergonicity is substantially lower (by ∼6 kcal/mol) in basic environment. © 2016 Wiley Periodicals, Inc.

A DFT kinetic study on 1,3-dipolar cycloaddition reactions in solution

Several popular density functional theory (DFT) methods have been employed to characterize a series of 1,3-dipolar cycloaddition reactions, including the exploration of reaction mechanisms and the calculations of kinetic parameters.

Quantum Chemical Studies on Detail Mechanism of Nitrosylation of NAMI-A-HSA Adduct

Hydrolysis of NAMI-A in NAMI-A-HSA (HSA = human serum albumin) and nitrosylation of hydrolyzed NAMI-A-HSA adduct have been studied in detail using density functional theory method. It has been observed that the chloride exchange reaction with water in the NAMI-A-HSA adduct follows an interchange dissociative mechanism passing through an unstable heptacoordinated activated complex. The computed free energy of activation (ΔG) and rate constant (k) for the hydrolysis process in aqueous medium are observed to be 24.85 kcal mol(-1) and 3.81 × 10(-6) s(-1), respectively. Nitrosylation of hydrolyzed NAMI-A-HSA adduct with nitric oxide is found to be thermodynamically more favorable with the incorporation of solvent effect and provides a detailed understanding related to the antimetastatic activity of the NAMI-A drug. This investigation shows that nitric oxide coordinates linearly to NAMI-A-HSA adduct leading to the reduction of ruthenium(III) to more active ruthenium(II), with the reduction potential of -2.32 V. Negative relative solvation and relative binding free energies suggest that the hydrolysis and nitrosylation reactions are found to be thermodynamically favorable and faster. Our computed results provide a detailed thermodynamics and kinetics which may be highly beneficial for understanding antimetastatic activity as well as the nitric oxide scavenging ability of NAMI-A.

A case study of the mechanism of alcohol-mediated Morita Baylis-Hillman reactions. The importance of experimental observations

The mechanism of the Morita Baylis-Hillman reaction has been heavily studied in the literature, and a long series of computational studies have defined complete theoretical energy profiles in these reactions. We employ here a combination of mechanistic probes, including the observation of intermediates, the independent generation and partitioning of intermediates, thermodynamic and kinetic measurements on the main reaction and side reactions, isotopic incorporation from solvent, and kinetic isotope effects, to define the mechanism and an experimental mechanistic free-energy profile for a prototypical Morita Baylis-Hillman reaction in methanol. The results are then used to critically evaluate the ability of computations to predict the mechanism. The most notable prediction of the many computational studies, that of a proton-shuttle pathway, is refuted in favor of a simple but computationally intractable acid-base mechanism. Computational predictions vary vastly, and it is not clear that any significant accurate information that was not already apparent from experiment could have been garnered from computations. With care, entropy calculations are only a minor contributor to the larger computational error, while literature entropy-correction processes lead to absurd free-energy predictions. The computations aid in interpreting observations but fail utterly as a replacement for experiment.

Water-soluble oxoglaucine-Y(iii), Dy(iii) complexes: in vitro and in vivo anticancer activities by triggering DNA damage, leading to S phase arrest and apoptosis

The complexes exhibited considerable in vitro and in vivo anticancer activity, and higher safety than ciplatin.

Potent Half-Sandwich Iridium(III) Anticancer Complexes Containing C∧N-Chelated and Pyridine Ligands

We report the synthesis and characterization of eight half-sandwich cyclopentadienyl IrIII pyridine complexes of the type [(η5-Cpxph)Ir(phpy)Z]PF6, in which Cpxph = C5Me4C6H5 (tetramethyl(phenyl)cyclopentadienyl), phpy = 2-phenylpyridine as C∧N-chelating ligand, and Z = pyridine (py) or a pyridine derivative. Three X-ray crystal structures have been determined. The monodentate py ligands blocked hydrolysis; however, antiproliferative studies showed that all the Ir compounds are highly active toward A2780, A549, and MCF-7 human cancer cells. In general the introduction of an electron-donating group (e.g., Me, NMe2) at specific positions on the pyridine ring resulted in increased antiproliferative activity, whereas electron-withdrawing groups (e.g., COMe, COOMe, CONEt2) decreased anticancer activity. Complex 5 displayed the highest anticancer activity, exhibiting submicromolar potency toward a range of cancer cell lines in the National Cancer Institute NCI-60 screen, ca. 5 times more potent than the clinical platinum(II) drug cisplatin. DNA binding appears not to be the major mechanism of action. Although complexes [(η5-Cpxph)Ir(phpy)(py)]+ (1) and [(η5-Cpxph)Ir(phpy)(4-NMe2-py)]+ (5) did not cause cell apoptosis or cell cycle arrest after 24 h drug exposure in A2780 human ovarian cancer cells at IC50 concentrations, they increased the level of reactive oxygen species (ROS) dramatically and led to a loss of mitochondrial membrane potential (ΔΨm), which appears to contribute to the anticancer activity. This class of organometallic Ir complexes has unusual features worthy of further exploration in the design of novel anticancer drugs.

Computational studies of the electronic, conductivities, and spectroscopic properties of hydrolysed Ru(II) anticancer complexes

The mechanism of activation of metal-based anticancer agents was reported to be through hydrolysis. In this study, computational method was used to gain insight to the correlation between the chemistry of the hydrolysis and the anticancer activities of selected Ru(II)-based complexes. Interestingly, we observed that the mechanism of activation by hydrolysis and their consequential anticancer activities is associated with favourable thermodynamic changes, higher hyperpolarizability (β), lower band-gap and higher first-order net current. The Fermi contact (FC) and spin dipole (SD) are found to be the two most significant Ramsey terms that determine the spin-spin couplings (J(HZ)) of most of the existing bonds in the complexes. Many of the computed properties give insights into the change in the chemistry of the complexes due to hydrolysis. Besides strong correlations of the computed properties to the anticancer activities of the complexes, using the quantum theory of atoms in a molecule (QTAIM) to analyse the spectroscopic properties shows a stronger correlation between the spectroscopic properties of Ru atom to the reported anticancer activities than the sum over of the spectroscopic properties of all atoms in the complexes.

Exploration of various electronic properties along the reaction coordinate for hydration of Pt(II) and Ru(II) complexes; the CCSD, MPx, and DFT computational study

In the study behavior of molecular electrostatic potential, averaged local ionization energy, and reaction electronic flux along the reaction coordinate of hydration process of three representative Ru(II) and Pt(II) complexes were explored using both post-HF and DFT quantum chemical approximations. Previously determined reaction mechanisms were explored by more detailed insight into changes of electronic properties using ωB97XD functional and MP2 method with 6-311++G(2df,2pd) basis set and CCSD/6-31(+)G(d,p) approach. The dependences of all examined properties on reaction coordinate give more detailed understanding of the hydration process.

A comparative DFT study on aquation and nucleobase binding of ruthenium (II) and osmium (II) arene complexes

The potential energy surfaces of the reactions of organometallic arene complexes of the type [(η (6)-arene)M(II)(pic)Cl] (where pic = 2-picolinic acid, M = Ru or Os) were examined by a DFT computational study. Among the seven density functional methods, hybrid exchange functional B3LYP outperforms the others to explain the aquation of the complexes. The reactions and binding energies of Ru(II) and Os(II) arene complexes with both 9EtG and 9EtA were studied to gain insight into the reactivity of these types of organometallic complexes with DNA. The obtained data rationalize experimental observation, contributing to partly understanding the potential biological and medical applications of organometallic complexes.