Davor Šakić

Prereactive complexes in chlorination of benzene, triazine, and tetrazine: a quantum chemical study.

In order to perform a complete search for prereactive complexes between arenes and chlorine, the stochastic search method was employed. Stationary points are optimized at B3LYP, M05-2X, and MP2 levels, while improved energetics are calculated using the B2PLYP-D method, which includes corrections important for accurate description of dispersion forces. New intermediates were located and their mechanistic relevance has been discussed. It has been suggested that, at least in the gas-phase, the T-shaped complex precedes the formation of classical benzene/chlorine π-complex. No σ-complex is found on the energy surface, unless explicit counterions are included in calculations. Neither π- nor σ-complexes were located on the reactant side of chlorination of triazine, but only linear and T-shaped complexes were identified as stable minima. These structures represent important prereactive complexes for chlorination of triazine. In the case of tetrazine, which is unlikely to undergo direct chlorination, only two complexes (resting and T-shaped) were located.

Chlorination of N-methylacetamide and amide-containing pharmaceuticals. Quantum-chemical study of the reaction mechanism.

Chlorination of amides is of utmost importance in biochemistry and environmental chemistry. Despite the huge body of data, the mechanism of reaction between amides and hypochlorous acid in aqueous environment remains unclear. In this work, the three different reaction pathways for chlorination of N-methylacetamide by HOCl have been considered: the one-step N-chlorination of the amide, the chlorination via O-chlorinated intermediate, and the N-chlorination of the iminol intermediate. The high-level quantum chemical G3B3 composite procedure, double-hybrid B2-PLYPD, B2K-PLYP methods, and global hybrid M06-2X and BMK methods have been employed. The calculated energy barriers have been compared to the experimental value of ΔG(#)298 ≈ 87 kJ/mol, which corresponds to reaction rate constant k(r) ≈ 0.0036 M(-1) s(-1). Only the mechanism in which the iminol form of N-methylacetamide reacts with HOCl is consistent (ΔG(#)298 = 87.3 kJ/mol at G3B3 level) with experimental results. The analogous reaction mechanism has been calculated as the most favorable pathway in the chlorination of small-sized amides and amide-containing pharmaceuticals: carbamazepine, acetaminophen, and phenytoin. We conclude that the formation of the iminol intermediate followed by its reaction with HOCl is the general mechanism of N-chlorination for a vast array of amides.

The chemical fate of paroxetine metabolites. Dehydration of radicals derived from 4-(4-fluorophenyl)-3-(hydroxymethyl)piperidine

Quantum chemical calculations have been used to model reactions which are important for understanding the chemical fate of paroxetine-derived radicals in the environment. In order to explain the experimental observation that the loss of water occurs along the (photo)degradation pathway, four different mechanisms of radical-induced dehydrations have been considered. The elimination of water from the N-centered radical cation, which results in the formation of an imine intermediate, has been calculated as the most feasible process. The predicted energy barrier (ΔG = 98.5 kJ mol(-1)) is within the barrier limits set by experimental measurements. All reaction intermediates and transition state structures have been calculated using the G3(MP2)-RAD composite procedure, and solvent effects have been determined using a mixed (cluster/continuum) solvation model. Several new products, which comply with the available experimental data, have been proposed. These structures could be relevant for the chemical fate of antidepressant paroxetine, but also for biologically and environmentally related substrates.

Structural and electronic determinants of flavonoid binding to human serum albumin: An extensive ligand-based study

Flavonoids are ubiquitous plant metabolites that interfere with different biological processes in the human organism. After absorption they bind to human serum albumin (HSA), the most abundant carrier protein in the blood which also binds various hormones and drugs. Binding of flavonoids to HSA may impact their distribution, influencing the active concentration in the blood. To determine the most prominent features responsible for binding of 20 different flavonoid aglycones to the IIA region of HSA, in vitro fluorescence measurements and density functional theory (DFT) calculations were conducted. These results were then integrated to elucidate structure–affinity relationships. The presented results reveal that flavones and flavonoles bind most strongly to the IIA region of HSA. There are several electronic and structural determinants associated with flavonoid binding to this HSA region: high C3 nucleophilicity and partial charge of O4, high HOMO and LUMO energies, and coplanarity of AC and B rings. Both steric and electronic characteristics of flavonoids have a great impact on their binding to HSA, with hydrogen donor and acceptor properties and coplanarity being the most prominent.

A quantum chemical study of HOCl-induced transformations of carbamazepine

The antiepileptic drug carbamazepine (CBZ) is one of the most persistent pharmaceuticals in the environment. Its chemical fate is influenced by the type of wastewater treatment. This study sets out to determine the degradation mechanism and products in the reaction between CBZ and hypochlorous acid (HOCl), which is the main chlorinating species in water. In the search for the most feasible pathways of HOCl-induced transformations of CBZ, a quantum chemical approach was employed. Chlorination and epoxidation of CBZ are two initial, competitive processes that result in two key intermediates: N-chloramide and 10,11-epoxide. The calculated free energy barriers (ΔG) for these reactions are 105.7 and 95.7 kJ mol-1 resp., which is in agreement with the experimental energy barrier of 98.2 kJ mol-1. All transformation products detected in chlorination experiments were located by computational models, and the reaction mechanism underlying their formation was described in detail. Different computational methods (density functional and ab initio theory) were applied, and the double hybrid B2-PLYPD functional was found to be superior in terms of efficiency and accuracy. Of special interest are oxoiminostilbene and formylacridine, which are the final products in the degradation cascade. Their exceptional thermodynamic stability, as predicted by quantum chemical methods, suggests that these structures should be considered as recalcitrants in chlorinated waters. Fruitful interplay between computational models and experimental data proves that the quantum chemical approach can be used as a predictive tool in environmental degradation studies.

Chlorination of 5-fluorouracil: Reaction mechanism and ecotoxicity assessment of chlorinated products

What happens to drugs in the chlorinating environment? Degradation products may vary in pharmacological profiles and in ecotoxicity potentials compared to the parent compound. This study combines synthesis, NMR spectroscopy, quantum chemical calculations, and toxicity experiments on Daphnia magna to investigate chemical fate of antineoplastic drug 5-fluorouracil (5-FU) in chlorinated environment, which is common in waste-water treatment procedures, but also endogenous in activated neutrophils. A reduction of toxicity (EC50 after 48 h is 50% higher than for the parent 5-FU) was observed after the first chlorination step, in which a chlorohydrin 5-chloro-5-fluoro-6-hydroxy-5,6-dihydrouracil was formed. Further chlorination leads to N-chlorinated intermediate, that undergoes the pyrimidine ring opening reaction. The final product, 2-chloro-2-fluoro-3,3-dihydroxypropanoic acid was obtained after the loss of the chlorinated urea fragment. This is the most potent compound in the reaction sequence, with toxicity parameter EC50, after 48 h, more than twice lower compared to the parent 5-FU. Clearly, the contact time between chlorinating species and degradation products provide different ecotoxicological properties of reaction mixtures. Interplay between experimental and theoretical procedures, to properly describe reaction pathways and provide more information on toxicity profiles, is a way forward in environmental science research.