The Role of Frontier Orbitals in Chemical Reactions (Nobel Lecture)

Differences in the Reaction Mechanisms of Chlorine Atom and Hydroxyl Radical with Organic Compounds: From Thermodynamics to Kinetics

Hydroxyl radical (HO•) and chlorine atom (Cl•) are common reactive species in aqueous environments. However, the intrinsic difference in their reactions with organic compounds has not been revealed. This study compared the reaction mechanisms of HO• and Cl• with 13 aromatic and 11 aliphatic compounds by quantum chemical calculation and laser flash photolysis. Both HO• and Cl• can spontaneously react with aromatic compounds via radical adduct formation (RAF), hydrogen atom transfer (HAT), and single electron transfer (SET) pathways. The SET reactions of Cl• were more thermodynamically favorable than HO•, but contrary results were obtained for HAT reactions. According to the free energy of activation (ΔGaq‡), the dominant oxidation mechanisms of aromatic compounds were RAF and SET by HO• and SET by Cl•. The important role of SET in the HO• reactions with aromatic compounds was further verified by accurately calculating the solvation free energy of HO•/HO- and experimentally tracking the radical cations, which were generally neglected in previous studies. Meanwhile, the ΔGaq‡ value of each reaction pathway of Cl• was lower than that of HO•, resulting in higher rate constants of Cl• with aromatic compounds than HO•. For saturated aliphatic compounds, HAT was found to be the only mechanism accounting for their transformation by HO• and Cl•. This study proposed general rules for the reaction mechanisms of HO• and Cl• and unraveled their differences in the aspects of thermodynamics and kinetics, providing fundamental information for understanding contaminant transformation in processes involving HO• and Cl•.

Unveiling the Unique Reactivity of Anionic Mn(I) Complexes via Metal-Ligand Cooperation: Nucleophilic Attack on C(sp3)-X Bonds

Metal-ligand cooperation (MLC) has emerged as a pivotal strategy for the catalytic activation of small molecules within both synthetic and biological arenas. Leveraging this approach, a suite of potent catalytic reactions─encompassing hydrogenation, hydroelementation, and dehydrogenative processes─have been realized, with notable advances in manganese catalysis in recent years. However, the activation of alkyl halides by Mn complexes, which typically requires strong reductants to form Mn(-I) complexes that are incompatible with standard cross-coupling conditions, remains a significant challenge. This limitation underscores the urgent need to investigate alternative methods for activating C(sp3)-X bonds using higher valence state Mn complexes. In response to this challenge, we present the synthesis, characterization, and reactivity of a new anionic Mn(I) complex featuring a redox-active dianionic ligand that induces multiple MLC functionalities. We have discovered an innovative mechanism of MLC, characterized by a single ligand transferring two electrons to the metal center. This novel process facilitates an orbital-symmetry-allowed nucleophilic attack on C(sp3)-X bonds, preserving manganese's oxidative state at +1. To the best of our knowledge, this is the first instance where the MLC strategy via a two-electron transfer process has been utilized to execute an SN2 nucleophilic attack at a C(sp3)-X bond by a relatively electron-deficient metal center like Mn(I). Additionally, the dianionic ligand of the anionic Mn(I) complex exhibits ambident nucleophilicity by reacting with different electrophiles, further highlighting its versatile reactivity.

Targeting cyclin-dependent kinase 2 CDK2: Insights from molecular docking and dynamics simulation - A systematic computational approach to discover novel cancer therapeutics

Global public health is confronted with significant challenges due to the prevalence of cancer and the emergence of treatment resistance. This work focuses on the identification of cyclin-dependent kinase 2 (CDK2) through a systematic computational approach to discover novel cancer therapeutics. A ligand-based pharmacophore model was initially developed using a training set of seven potent CDK2 inhibitors. The obtained most robust model was characterized by three features: one donor (|Don|) and two acceptors (|Acc|). Screening this model against the ZINC database resulted in identifying 108 hits, which underwent further molecular docking studies. The docking results indicated binding affinity, with energy values ranging from -6.59 kcal mol⁻¹ to -7.40 kcal mol⁻¹ compared to the standard Roscovitine. The top 10 compounds (Z1-Z10) selected from the docking data were further screened for ADMET profiling, ensuring their compliance with pharmacokinetic and toxicological criteria. The top 3 compounds (Z1-Z3) chosen from the docking were subjected to Density Functional Theory (DFT) studies. They revealed significant variations in electronic properties, providing insights into the reactivity, stability, and polarity of these compounds. Molecular dynamics simulations confirmed the stability of the ligand-protein complexes, with acceptable RMSD and RMSF values. Specifically, compound Z1 demonstrated stability, around 2.4 Å, and maintained throughout the 100 ns simulation period with minimal conformational changes, stable RMSD, and consistent protein-ligand interactions.

Synthesis, structural characterization, DFT and molecular dynamics simulations of dinuclear (μ-hydroxo)-bridged triethanolamine copper(II) complexes: efficient candidates towards visible light-mediated photo-Fenton degradation of organic dyes

Multinuclear (di/tri) copper(II) complexes bridged through hydroxyl groups are very interesting coordination complexes owing to their potential applications in various fields. In this work, three novel dinuclear (μ-hydroxo)-bridged copper(II) complexes in the crystal form, namely, [Cu2(3,5-DIFLB)2(H2tea)2](H2O) (1), [Cu2(4-ClB)2(H2tea)2](H2O) (2), and [Cu2(4-ETHB)2(H2tea)2](H2O)2 (3) (where DIFLB = difluorobenzoate, CLB = chlorobenzoate, ETHB = ethoxybenzoate, and H3tea = triethanolamine), were isolated at room temperature using methanol and water in a 4 : 1 v/v ratio as a solvent. Furthermore, all three complexes (1-3) were characterised using spectroscopic (UV-vis, DRS, and FT-IR), electrochemical (CV) and single-crystal X-ray diffraction techniques. Structural insights gained by packing analysis revealed the role of steric constraints of substituents and various non-covalent interactions in lattice stabilization, which were indeed supported by theoretical and molecular electrostatic potential illustrations. Hirshfeld surface analysis provided quantitative verification about various non-covalent interactions (interatomic contacts) involved in the packing of molecules. Interestingly, as a potential application, complexes 1-3 all exhibited remarkable visible light-mediated photo-Fenton degradation of approximately 98% for 50 ppm concentration of organic dyes (fuchsin basic (FB) and methyl orange (MO)) in 90 minutes with the optimized conditions of 1 mg mL-1 of dye solution. In all the cases, dye degradation by these materials was ascribed to the symbiotic relations among the molecular structures of complexes 1-3, which were endowed with various electron-withdrawing and electron-releasing substituents and ionic strength, with respect to the structure, shape and interacting patterns of dye molecules. The adsorption mechanism indicates that various weak interactions between the donor and acceptor groups of complexes and dyes, such as electrostatic, hydrogen bonding, and direct coordination to metal sites, play a crucial role, which is confirmed by molecular dynamics (MD) simulations. Theoretical studies by DFT-based descriptors, molecular electrostatic potentials, and band gaps provided deep insights into various electronic and reactivity parameters. For subsequent processes of dye degradation, complexes 1-3 were stable and recoverable. The successful integration of experimental and theoretical approaches sheds light on copper-based dinuclear stable coordination complexes, showcasing a significant step towards the development of novel heterogeneous photo-Fenton catalysts.

(E)-2-Benzylidenecyclanones: Part XIX. Reaction of (E)-2-(4'-X-Benzylidene)-1-tetralones with Cellular Thiols: Comparison of Thiol Reactivities of Open-Chain Chalcones and Their Six- and Seven-Membered Cyclic Analogs

Non-enzyme-catalyzed thiol addition onto the α,β-unsaturated carbonyl system is associated with several biological effects. Kinetics and diastereoselectivity of non-enzyme catalyzed nucleophilic addition of reduced glutathione (GSH) and N-acetylcysteine (NAC) to the six-membered cyclic chalcone analogs 2a and 2b were investigated at different pH values (pH 3.2, 7.4 and 8.0). The selected compounds displayed in vitro cancer cell cytotoxicity (IC50) of different orders of magnitude. The chalcones intrinsically reacted with both thiols under all incubation conditions. The initial rates and compositions of the final mixtures depended both on the substitution and the pH. The stereochemical outcome of the reactions was evaluated using high-pressure liquid chromatography with UV detection (HPLC-UV). The structures of the formed thiol-conjugates and the retro-Michael products (Z)-2a and (Z)-2b were confirmed by high-pressure liquid chromatography-mass spectrometry (HPLC-MS). Frontier molecular orbitals and the Fukui function calculations were carried out to investigate their effects on the six-membered cyclic analogs. Data were compared with those obtained with the open-chain (1) and the seven-membered (3) analogs. The observed reactivities do not directly relate to the difference in in vitro cancer cell cytotoxicity of the compounds.

Adsorption studies of isoxazole derivatives as corrosion inhibitors for mild steel in 1M HCl solution: DFT studies and molecular dynamics simulation

CONTEXT

The corrosion of mild steel is a significant issue in various industries, prompting the need for effective corrosion inhibitors. This study focuses on understanding the corrosion inhibition properties of organic compounds derived from isoxazole, namely series Iso(a), Iso(b), Iso(c), Iso(d), Iso(e), Iso(f), Iso(g), and Iso(h), which could have implications for materials science and industrial applications. By investigating the influence of different substitutions on these compounds, valuable insights can be gained into designing better corrosion inhibitors for practical use.

METHODS

Theoretical studies were conducted using density functional theory (DFT) with the B3LYP functional and the 6-31G (d,p) basis set. These calculations enabled the evaluation of various parameters including frontier orbital energies (EHOMO, ELUMO), energy gap (∆E), electronegativity (χ), absolute hardness (η), softness (σ), fraction of transferred electrons (∆N), as well as local properties such as natural atomic populations and Fukui indices. Additionally, molecular dynamics simulations were performed to study the adsorption behavior of the inhibitors on the surface of Fe (110). The simulations were conducted using Materials Studio version 8.0 software package using COMPASS force field to understand the impact of different functional groups on the inhibitors before and after adsorption on the iron surface.

First insights into the bioaccumulation, biotransformation and trophic transfer of typical tetrabromobisphenol A (TBBPA) analogues along a simulated aquatic food chain

Tetrabromobisphenol A (TBBPA) analogues have been investigated for their prevalent occurrence in environments and potential hazardous effects to humans and wildlife; however, there is still limited knowledge regarding their toxicokinetics and trophic transfer in aquatic food chains. Using a developed toxicokinetic model framework, we quantified the bioaccumulation, biotransformation and trophic transfer of tetrabromobisphenol S (TBBPS) and tetrabromobisphenol A di(allyl ether) (TBBPA-DAE) during trophic transfer from brine shrimp (Artemia salina) to zebrafish (Danio rerio). The results showed that the two TBBPA analogues could be readily accumulated by brine shrimp, and the estimated bioconcentration factor (BCF) value of TBBPS (5.68 L kg-1 ww) was higher than that of TBBPA-DAE (1.04 L kg-1 ww). The assimilation efficiency (AE) of TBBPA-DAE in zebrafish fed brine shrimp was calculated to be 16.3%, resulting in a low whole-body biomagnification factor (BMF) in fish (0.684 g g-1 ww). Based on the transformation products screened using ultra-high-performance liquid chromatograph-high resolution mass spectrometry (UPLC-HRMS), oxidative debromination and hydrolysis were identified as the major transformation pathways of TBBPS, while the biotransformation of TBBPA-DAE mainly took place through ether bond breaking and phase-II metabolism. Lower accumulation of TBBPA as a metabolite than its parent chemical was observed in both brine shrimp and zebrafish, with metabolite parent concentration factors (MPCFs) < 1. The investigated BCFs for shrimp of the two TBBPA analogues were only 3.77 × 10-10 - 5.59 × 10-3 times of the theoretical Kshrimp-water based on the polyparameter linear free energy relationships (pp-LFERs) model, and the BMF of TBBPA-DAE for fish was 0.299 times of the predicted Kshrimp-fish. Overall, these results indicated the potential of the trophic transfer in bioaccumulation of specific TBBPA analogues in higher trophic-level aquatic organisms and pointed out biotransformation as an important mechanism in regulating their bioaccumulation processes. ENVIRONMENTAL IMPLICATION: The internal concentration of a pollutant in the body determines its toxicity to organisms, while bioaccumulation and trophic transfer play important roles in elucidating its risks to ecosystems. Tetrabromobisphenol A (TBBPA) analogues have been extensively investigated for their adverse effects on humans and wildlife; however, there is still limited knowledge regarding their toxicokinetics and trophic transfer in aquatic food chains. This study investigated the bioaccumulation, biotransformation and trophic transfer of TBBPS and TBBPA-DAE in a simulated di-trophic food chain. This state-of-art study will provide a reference for further research on this kind of emerging pollutant in aquatic environments.

Synthesis, characterization, biological activity and computation-based efficacy of cobalt(II) complexes of biphenyl-2-ol against SARS-CoV-2 virus

Cobalt(II) complexes of biphenyl-2-ol of composition, CoCl2-n(OC6H4C6H5-2)n(H2O)4 (where n = 1 or 2), were prepared by reacting cobaltous(II) chloride with equi- and bimolar ratios of sodium salt of biphenyl-2-ol. The structural characterization of the synthesized complexes was accomplished by NMR, FTIR, thermogravimetry (TGA), high resolution mass spectroscopy (HRMS), electronic spectroscopic techniques coupled with density functional theory (DFT). The stability of the complexes in different pH media of solvent was studied. Chemical reactivity parameters of the newly synthesized complexes, computed using DFT, indicated greater reactivity of complex 2 over complex 1 and free ligand as indicated by its low HOMO-LUMO energy gap corresponding to 1.71 eV. Molecular docking (MD) studies were carried out in order to study the binding affinities between amino acid residues of DNA duplex (PDB ID: 1BNA) and SARS-CoV-2 (PDB ID: 7T9K) with newly synthesized complexes. Complex 2 has shown promising antivirus behaviour with an inhibition constant value of 0.0423 µmol-1 with amino acid residues of SARS-CoV-2 virus. Toxicity of the complexes was predicted using ProTox-II online server. Antibacterial studies have indicated the complexes to exhibit greater efficacy than the free ligand, while the antioxidant activities have suggested them to display enhanced antioxidant behaviour as compared to reference compounds.Communicated by Ramaswamy H. Sarma.

Theoretical elucidation of the structure, bonding, and reactivity of the CaMn4Ox clusters in the whole Kok cycle for water oxidation embedded in the oxygen evolving center of photosystem II. New molecular and quantum insights into the mechanism of the O-O bond formation

This paper reviews our historical developments of broken-symmetry (BS) and beyond BS methods that are applicable for theoretical investigations of metalloenzymes such as OEC in PSII. The BS hybrid DFT (HDFT) calculations starting from high-resolution (HR) XRD structure in the most stable S1 state have been performed to elucidate structure and bonding of whole possible intermediates of the CaMn4Ox cluster (1) in the Si (i = 0 ~ 4) states of the Kok cycle. The large-scale HDFT/MM computations starting from HR XRD have been performed to elucidate biomolecular system structures which are crucial for examination of possible water inlet and proton release pathways for water oxidation in OEC of PSII. DLPNO CCSD(T0) computations have been performed for elucidation of scope and reliability of relative energies among the intermediates by HDFT. These computations combined with EXAFS, XRD, XFEL, and EPR experimental results have elucidated the structure, bonding, and reactivity of the key intermediates, which are indispensable for understanding and explanation of the mechanism of water oxidation in OEC of PSII. Interplay between theory and experiments have elucidated important roles of four degrees of freedom, spin, charge, orbital, and nuclear motion for understanding and explanation of the chemical reactivity of 1 embedded in protein matrix, indicating the participations of the Ca(H2O)n ion and tyrosine(Yz)-O radical as a one-electron acceptor for the O-O bond formation. The Ca-assisted Yz-coupled O-O bond formation mechanisms for water oxidation are consistent with recent XES and very recent time-resolved SFX XFEL and FTIR results.

The Nature of the Chemical Bonds of High-Valent Transition-Metal Oxo (M=O) and Peroxo (MOO) Compounds: A Historical Perspective of the Metal Oxyl-Radical Character by the Classical to Quantum Computations

This review article describes a historical perspective of elucidation of the nature of the chemical bonds of the high-valent transition metal oxo (M=O) and peroxo (M-O-O) compounds in chemistry and biology. The basic concepts and theoretical backgrounds of the broken-symmetry (BS) method are revisited to explain orbital symmetry conservation and orbital symmetry breaking for the theoretical characterization of four different mechanisms of chemical reactions. Beyond BS methods using the natural orbitals (UNO) of the BS solutions, such as UNO CI (CC), are also revisited for the elucidation of the scope and applicability of the BS methods. Several chemical indices have been derived as the conceptual bridges between the BS and beyond BS methods. The BS molecular orbital models have been employed to explain the metal oxyl-radical character of the M=O and M-O-O bonds, which respond to their radical reactivity. The isolobal and isospin analogy between carbonyl oxide R2C-O-O and metal peroxide LFe-O-O has been applied to understand and explain the chameleonic chemical reactivity of these compounds. The isolobal and isospin analogy among Fe=O, O=O, and O have also provided the triplet atomic oxygen (3O) model for non-heme Fe(IV)=O species with strong radical reactivity. The chameleonic reactivity of the compounds I (Cpd I) and II (Cpd II) is also explained by this analogy. The early proposals obtained by these theoretical models have been examined based on recent computational results by hybrid DFT (UHDFT), DLPNO CCSD(T0), CASPT2, and UNO CI (CC) methods and quantum computing (QC).

Understanding the reactivity and selectivity of Diels-Alder reactions involving furans

The reactivity and endo/exo selectivity of the Diels-Alder cycloaddition reactions involving furan and substituted furans as dienes have been computationally explored. In comparison to cyclopentadiene, it is found that furan is comparatively less reactive and also less endo-selective in the reaction with maleic anhydride as the dienophile. Despite that, both the reactivity and the selectivity can be successfully modified by the presence of substituents at either 2- or 3-positions of the heterocycle. In this sense, it is found that the presence of strong electron-donor groups significantly increases the reactivity of the system while the opposite is found in the presence of electron-withdrawing groups. The observed trends in both the reactivity and selectivity are analyzed quantitatively in detail by means of the activation strain model of reactivity in combination with the energy decomposition analysis methods.

Structural basis of antibacterial photodynamic action of curcumin against S. aureus

Antimicrobial photodynamic therapy (aPDT) is an alternative tool to commercial antibiotics for the inactivation of pathogenic bacteria (e.g., S. aureus). However, there is still a lack of understanding of the molecular modeling of the photosensitizers and their mechanism of action through oxidative pathways. Herein, a combined experimental and computational evaluation of curcumin as a photosensitizer against S. aureus was performed. The radical forms of keto-enol tautomers and the energies of curcumin's frontier molecular orbitals were evaluated by density functional theory (DFT) to point out the photodynamic action as well as the photobleaching process. Furthermore, the electronic transitions of curcumin keto-enol tautomers were undertaken to predict the transitions as a photosensitizer during the antibacterial photodynamic process. Moreover, molecular docking was used to evaluate the binding affinity with the S. aureus tyrosyl-tRNA synthetase as the proposed a target for curcumin. In this regard, the molecular orbital energies show that the curcumin enol form has a character of 4.5% more basic than the keto form - the enol form is a more promising electron donor than its tautomer. Curcumin is a strong electrophile, with the enol form being 4.6% more electrophilic than its keto form. In addition, the regions susceptible to nucleophilic attack and photobleaching were evaluated by the Fukui function. Regarding the docking analysis, the model suggested that four hydrogen bonds contribute to the binding energy of curcumin's interaction with the ligand binding site of S. aureus tyrosyl-tRNA synthetase. Finally, residues Tyr36, Asp40, and Asp177 contact curcumin and may contribute to orienting the curcumin in the active area. Moreover, curcumin presented a photoinactivation of 4.5 log unit corroborating the necessity of the combined action of curcumin, light, and O2 to promote the photooxidation damage of S. aureus. These computational and experimental data suggest insights regarding the mechanism of action of curcumin as a photosensitizer to inactivate S. aureus bacteria.

Origin of the Bürgi-Dunitz Angle

The Bürgi-Dunitz (BD) angle plays a pivotal role in organic chemistry to rationalize the nucleophilic addition to carbonyl groups. Yet, the origin of the obtuse trajectory of the nucleophile remains incompletely understood. Herein, we quantify the importance of the underlying physical factors quantum chemically. The obtuse BD angle appears to originate from the concerted action of a reduced Pauli repulsion between the nucleophile HOMO and carbonyl π bond, a more stabilizing HOMO-π*-LUMO(C=O) interaction, as well as a more favorable electrostatic attraction.

Modulating Escape Channels of Cycloheptatrienyl Rhodium Carbenes To Form Semibullvalene

We describe the various escape channels available to dirhodium carbene intermediates from cycloheptatrienyl diazo compounds located with density functional theory. An intramolecular cyclopropanation would, in principle, provide a new route to semibullvalenes (SBVs). A detailed exploration of the potential energy surface reveals that methylating carbon-7 suppresses a competing β-hydride migration pathway to heptafulvene products, giving SBV formation a reasonable chance. During our explorations, we additionally discovered unusual spirononatriene, spironorcaradiene, and metal-stabilized 9-barbaralyl cation structures as local minima.

Is DFT Accurate Enough to Calculate Regioselectivity? The Case of 1,3-Dipolar Cycloaddition of Azide to Alkynes and Alkenes

The importance of regioselectivity in 1,3-dipolar cycloadditions (DCs) makes it surprising that no benchmarking study on this problem has appeared. We investigated whether DFT calculations are an accurate tool to predict the regioselectivity of uncatalyzed thermal azide 1,3-DCs. We considered the reaction between HN3 and 12 dipolarophiles, comprising ethynes HC≡C-R and ethenes H2 C=CH-R (R=F, OH, NH2 , Me, CN, CHO), which cover a broad range of electron demand and conjugation ability. We established benchmark data by the W3X protocol [complete-basis-set-extrapolated CCSD(T)-F12 energy with T-(T) and (Q) corrections and MP2-calculated core/valence and relativistic effects] and showed that core/valence effects and high-order excitations are important for accurate regioselectivity. Regioselectivities calculated using an extensive set of density functional approximations (DFAs) were compared with benchmark data. Range-separated and meta-GGA hybrids gave the best results. Good treatment of self-interaction and electron exchange are the key features for accurate regioselectivity. Dispersion correction slightly improves agreement with W3X results. The best DFAs provide the isomeric TS energy difference with an expected error ≈0.7 mh and errors ≈2 mh can occur. The isomer yield provided by the best DFA has an expected error of ±5 %, though errors up to 20 % are not rare. At present, an accuracy of 1-2 % is unfeasible but it seems that we are not far from achieving this goal.

Synthesis, crystal structure, Hirshfeld surface analysis and DFT study of the 1,1'-(buta-1,3-diyne-1,4-di-yl)bis-(cyclo-hexan-1-ol)

The title compound, C16H22O2, was synthesized in order to obtain its guest-free form because 'wheel-and-axle'-shaped mol-ecules tend to crystallize from solutions as solvates or host-guest mol-ecules. It crystallizes in the monoclinic space group P2/c with two crystallographically non-equivalent mol-ecules, one situated on an inversion center and the other on a twofold axis. The rod-like 1,3-diyne fragments have the usual linear geometry. In the crystal, O-H ⋯ O bonds form eight-membered rings of the R 4 4(8) type, linking mol-ecules into layers. The Hirshfeld surface analysis indicates that the largest con-tributions are from inter-molecular H⋯H (ca 71%) and H⋯C/C⋯H (ca 19%) contacts. The energies of the frontier mol-ecular orbitals were determined by DFT calculations at the B3LYP/def2-TZVP level of theory.

(E)-2-Benzylidenecyclanones: Part XVIII Study the Possible Link between Glutathione Reactivity and Cancer Cell Cytotoxic Effects of Some Cyclic Chalcone Analogs A Comparison of the Reactivity of the Open-Chain and the Seven-Membered Homologs

Non-enzymatic thiol addition into the α,β-unsaturated carbonyl system is associated with several biological effects. In vivo, the reactions can form small-molecule thiol (e.g., glutathione) or protein thiol adducts. The reaction of two synthetic (4'-methyl- and 4'-methoxy substituted) cyclic chalcone analogs with reduced glutathione (GSH) and N-acetylcysteine (NAC) was studied by (high-pressure liquid chromatography-ultraviolet spectroscopy) HPLC-UV method. The selected compounds displayed in vitro cancer cell cytotoxicity (IC₅₀) of different orders of magnitude. The structure of the formed adducts was confirmed by (high-pressure liquid chromatography-mass spectrometry) HPLC-MS. The incubations were performed under three different pH conditions (pH 3.2/3.7, 6.3/6.8, and 8.0/7.4). The chalcones intrinsically reacted with both thiols under all incubation conditions. The initial rates and compositions of the final mixtures depended on the substitution and the pH. The frontier molecular orbitals and the Fukui function were carried out to investigate the effects on open-chain and seven-membered cyclic analogs. Furthermore, machine learning protocols were used to provide more insights into physicochemical properties and to support the different thiol-reactivity. HPLC analysis indicated diastereoselectivity of the reactions. The observed reactivities do not directly relate to the different in vitro cancer cell cytotoxicity of the compounds.

Predictive models for the aqueous phase reactivity of inorganic radicals with organic micropollutants

Single-electron transfer (SET) is one of the most common reaction mechanisms for degrading organic micropollutants (OMPs) in advanced oxidation processes. We collected 300 SET reactions (CO₃^•-, SO₄^•-, Cl₂^•-, and Br₂^•--mediated) and calculated three key parameters for understanding the SET mechanism: aqueous phase free energies of activation (ΔG^‡), free energies of reactions (ΔG), and orbital energy gaps of reactants (E_OMPsHOMO-E_RadiLUMO). Then, we classified the OMPs according to their structure, developed and evaluated linear energy relationships of the second-order rate constants (k) with ΔG^‡, ΔG, or E_OMPsHOMO-E_RadiLUMO in each class. Considering that a single descriptor cannot capture all the chemical diversity, we combined ΔG^‡, ΔG, and E_OMPsHOMO-E_RadiLUMO as inputs to develop multiple linear regression (MLR) models. Chemical classification is critical to the linear model described above. However, OMPs usually have multiple functional groups, making the classification challenging and uncertain. Therefore, we tried machine learning algorithms to predict k values without chemical classification. We found that decision trees (R² = 0.88-0.95) and random forest (R² = 0.90-0.94) algorithms show better performance on the prediction of the k values, whereas boosted tree algorithm cannot make an accurate prediction (R² = 0.19-0.36). Overall, our study provides a powerful tool to predict the aqueous phase reactivity of OMP to certain radicals without the need for chemical classification.

Unraveling Bonding Mechanisms and Electronic Structure of Pyridine Oximes on Fe(110) Surface: Deeper Insights from DFT, Molecular Dynamics and SCC-DFT Tight Binding Simulations

The development of corrosion inhibitors with outstanding performance is a never-ending and complex process engaged in by researchers, engineers and practitioners. The computational assessment of organic corrosion inhibitors' performance is a crucial step towards the design of new task-specific materials. Herein, the electronic features, adsorption characteristics and bonding mechanisms of two pyridine oximes, namely 2-pyridylaldoxime (2POH) and 3-pyridylaldoxime (3POH), with the iron surface were investigated using molecular dynamics (MD), and self-consistent-charge density-functional tight-binding (SCC-DFTB) simulations. SCC-DFTB simulations revealed that the 3POH molecule can form covalent bonds with iron atoms in its neutral and protonated states, while the 2POH molecule can only bond with iron through its protonated form, resulting in interaction energies of -2.534, -2.007, -1.897, and -0.007 eV for 3POH, 3POH⁺, 2POH⁺, and 2POH, respectively. Projected density of states (PDOSs) analysis of pyridines-Fe(110) interactions indicated that pyridine molecules were chemically adsorbed on the iron surface. Quantum chemical calculations (QCCs) revealed that the energy gap and Hard and Soft Acids and Bases (HSAB) principles were efficient in predicting the bonding trend of the molecules investigated with an iron surface. 3POH had the lowest energy gap of 1.706 eV, followed by 3POH⁺ (2.806 eV), 2POH⁺ (3.121 eV), and 2POH (3.431 eV). In the presence of a simulated solution, MD simulation showed that the neutral and protonated forms of molecules exhibited a parallel adsorption mode on an iron surface. The excellent adsorption properties and corrosion inhibition performance of 3POH may be attributed to its low stability compared to 2POH molecules.

New Insights on Glutathione's Supramolecular Arrangement and Its In Silico Analysis as an Angiotensin-Converting Enzyme Inhibitor

Angiotensin-converting enzyme (ACE) inhibitors are one of the most active classes for cardiovascular diseases and hypertension treatment. In this regard, developing active and non-toxic ACE inhibitors is still a continuous challenge. Furthermore, the literature survey shows that oxidative stress plays a significant role in the development of hypertension. Herein, glutathione's molecular structure and supramolecular arrangements are evaluated as a potential ACE inhibitor. The tripeptide molecular modeling by density functional theory, the electronic structure by the frontier molecular orbitals, and the molecular electrostatic potential map to understand the biochemical processes inside the cell were analyzed. The supramolecular arrangements were studied by Hirshfeld surfaces, quantum theory of atoms in molecules, and natural bond orbital analyses. They showed distinct patterns of intermolecular interactions in each polymorph, as well as distinct stabilizations of these. Additionally, the molecular docking study presented the interactions between the active site residues of the ACE and glutathione via seven hydrogen bonds. The pharmacophore design indicated that the hydrogen bond acceptors are necessary for the interaction of this ligand with the binding site. The results provide useful information for the development of GSH analogs with higher ACE inhibitor activity.

DFT study of therapeutic potential of graphitic carbon nitride as a carrier for controlled release of melphalan: an anticancer drug

In the present research, the drug-delivery efficiency of graphitic carbon nitride (g-CN) for melphalan (an anti-cancer drug) was evaluated. To investigate the efficacy of g-CN as a drug-delivery system, the electronic properties of melphalan drug, g-CN, and g-CN-melphalan were calculated at the ground and excited states. The adsorption energy calculated for g-CN-melphalan complex in the water phase is - 1.51 eV. The interactions between g-CN and melphalan were investigated by a non-covalent interactions (NCl) analysis, which showed that there were weak interactions between g-CN and melphalan drug. These low intermolecular forces will allow for easy off-loading of the melphalan at the targeted site. Frontier molecular-orbitals (FMOs) analysis showed that the charge was transferred from melphalan to g-CN during the excitation process. Charge transfer was studied by charge decomposition analysis. Calculations at the excited state revealed that the g-CN-melphalan complex's λ_max showed a redshift of 15 nm and 39 nm in the gas and water phase, respectively. The photoinduced electron transfer (PET) process was studied for 1-2 excited state by using electron hole theory. PET process suggests that fluorescence quenching may take place. The findings demonstrated that g-CN can be used as a drug-delivery system for melphalan drug to treat cancer. This investigation may also encourage more consideration of different 2D substances for drug delivery.

Polymyxin B Activity Rescue by (−)-Camphene-Based Thiosemicarbazide Against Carbapenem-Resistant Enterobacterales

Due to the significant shortage of therapeutic options for carbapenem-resistant Enterobacterales (CRE) infections, new drugs or therapeutic combinations are urgently required. We show in this study that (-)-camphene-based thiosemicarbazide (TSC) may act synergistically with polymyxin B (PMB) against CRE, rescuing the activity of this antimicrobial. With the specific aim of a better molecular understanding of this effect caused by the presence of TSC, theoretical calculations were also performed in this study. Based on these findings, it is concluded that the presence of TSC moieties contributes to significant changes in the hydrogen atom charge of PMB structure, which trend more positives for the PMB/TSC system studied. This could lead to the formation of stronger hydrogen bonds in the Enterobacterales active site and, thus contribute to a molecular understanding of the PMB rescue of activity promoted by the presence of TSC moiety. As such, the clinical potential of these drug combinations requires further evaluation.

Reactivity of 4,5-Dichlorophthalic Anhydride towards Thiosemicarbazide and Amines: Synthesis, Spectroscopic Analysis, and DFT Study

The cyclic anhydrides are broadly employed in several fields, such as the chemical, plastic, agrochemical, and pharmaceutical industries. This study describes the chemical reactivity of 4,5-dichlorophthalic anhydride towards several nucleophiles, including thiosemicarbazide and different amines, to produce the carboxylic acid derivatives resulting from anhydride’s opening, namely, phthalimide and dicarboxylic acid (1–12) products. Their chemical structures are confirmed by NMR, IR and MS spectra analyses. Density–functional theory (DFT) studies are performed using (DFT/B3LYP) with the 6-311G(d, p) basis sets to recognize different chemical and physical features of the target compounds.

Gamma radiolytic stability of the novel modified diglycolamide 2,2'-oxybis(N,N-didecylpropanamide) (mTDDGA) for grouped actinide extraction

Reprocessing of spent nuclear fuel aims at improving resource efficiency and reducing its radiotoxicity and heat production in the long term. The necessary separation of certain metal ions from the spent fuel solutions can be achieved using different solvent extraction processes. For the scenario of the EURO-GANEX process, the use of the new, modified diglycolamide 2,2′-oxybis(N,N-didecylpropanamide) (mTDDGA) was recently proposed to simplify the current solvent composition and reduce extraction of fission products. Before further developing the process based on this new ligand, its stability under ionizing radiation conditions needs to be studied. For this reason, gamma irradiation experiments were conducted followed by analyses with high performance liquid chromatography coupled to a mass spectrometer (HPLC-MS). The determined degradation rate of mTDDGA was found to be lower than that of the reference molecule N,N,N′,N′-tetra-n-octyl-diglycolamide (TODGA). Many identified degradation compounds of both molecules are analogues showing the same bond breaking, although also unreported de-methylation, double/triple de-alkylation and n-dodecane addition products were observed.

The radiolysis behavior of a new diglycolamide for solvent extraction of actinides and lanthanides was studied. The observed degradation rate was lower than for the reference molecule and 22 degradation compounds were identified.

Exploring the maximum Fukui function sites with the frontier-controlled soft-soft reactions using 1,3-dipolar cycloaddition reactions of nitrilium betaines

According to P.K. Chattaraj. J Phys Chem A 2001, 105, 511-513 "the maximum Fukui function site is the best for the frontier-controlled soft-soft reactions whereas for the charge-controlled hard-hard interactions the preferred site is associated with the maximum net charge and not necessarily the minimum Fukui function". Taking into account these outcomes in this research is explored this reactivity scheme using in first case the reaction between fulminic acid with ethylene (reference reaction), after is varying the dipolarophile in the reaction between fulminic acid with acetylene, and finally is varying the dipole in the reaction between formonitrile imine with ethylene. These results allow study parameter such as charge transfer, polarizability, covalent character on bonding, among other; also shown the preference by the sf^- / sf⁺ interactions in the transition state on the sf^- / sf^- interactions. On the other hand, these results also were justified using net electrophilicity which is defined as the electrophilic power of a system relative to its own nucleophilic power.

Kenichi Fukui, Frontier Molecular Orbital Theory, and the Woodward-Hoffmann Rules. Part I. The Person†

Kenichi Fukui shared the 1981 Nobel Prize in Chemistry with Roald Hoffmann for "their theories, developed independently, concerning the course of chemical reactions." This is Paper 4 - Part I, of a three-part trilogy within a 27-paper series on the history of the development of the Woodward-Hoffmann rules. The personality and professional style of Fukui is discussed as well as his style of teaching and his mottos and sayings. A brief chronological history of his academic and professional life is presented along with a listing of Fukui's great flexibility in behaviors and in science. Lastly, a short list of Fukui's key awards is provided.

Molecular Modelling of Optical Biosensor Phosphorene-Thioguanine for Optimal Drug Delivery in Leukemia Treatment

Thioguanine is an anti-cancer drug used for the treatment of leukemia. However, thioguanine has weak aqueous solubility and low biocompatibility, which limits its performance in the treatment of cancer. In the present work, these inadequacies were targeted using density functional theory-based simulations. Three stable configurations were obtained for the adsorption of thioguanine molecules on the phosphorene surface, with adsorption energies in the range of -76.99 to -38.69 kJ/mol, indicating physisorption of the drug on the phosphorene surface. The calculated bandgap energies of the individual and combined geometries of phosphorene and thioguanine were 0.97 eV, 2.81 eV and 0.91 eV, respectively. Owing to the physisorption of the drug molecule on the phosphorene surface, the bandgap energy of the material had a direct impact on optical conductivity, which was significantly altered. All parameters that determine the potential ability for drug delivery were calculated, such as the dipole moment, chemical hardness, chemical softness, chemical potential, and electrophilicity index. The higher dipole moment (1.74 D) of the phosphorene-thioguanine complex reflects its higher biodegradability, with no adverse physiological effects.

Resonant Electron Spectroscopy: Identification of Atomic Contributions to Valence States

Valence electronic structure is crucial for understanding and predicting reactivity. Valence non-resonant X-ray photoelectron spectroscopy (NRXPS) provides a direct method for probing the overall valence electronic structure. However, it is...

Effect of Metal Complexation on Diglycolamides Radiolysis: A Comparison between Ex-Situ Gamma and In-Situ Alpha Irradiation

Radiolytic degradation is an important aspect to consider when developping a ligand or a complexant for radionucleides. Diglycolamide extractants (DGAs) have been playing an important role in many partition processes...

Synthesis, crystal structure, computational study and anti-virus effect of mixed ligand copper (II) complex with ONS donor Schiff base and 1, 10-phenanthroline

This work deals with the synthesis, crystal structure, computational study and antiviral potential of mixed ligand copper(II) complex [Cu(L)(phen)](1), (where, H₂L = (Z)-N'-((E)-2-hydroxy-3,5-diiodobenzylidene)-N,N-dimethylcarbamohydrazonothioic acid, phen = 1,10-phenanthroline). The Schiff base ligand (H₂L) is coordinated with Cu(II) ion in O, N, S-tridentate mode. The copper complex (1) crystallized in the monoclinic system of the space group P21/c with eight molecules in the unit cell and reveals a square pyramidal geometry. Furthermore, we also perform quantum chemical calculations to get insights into the structure-property relationship and functional properties of ligand (H₂L) and its copper (II) complex [Cu(L)(phen)](1). Complex [Cu(L)(phen)](1) was also virtually designed in-silico evaluation by Swiss-ADME. Additionally, inspiring by recent developments to find a potential inhibitor for the COVID-19 virus, we have also performed molecular docking study of ligand and its copper complex (1) to see if our compounds shows an affinity for the main protease (M^pro) of COVID-19 spike protein (PDB ID: 7C8U). Interestingly, the results are found quite encouraging where the binding affinity and inhibition constant were found to be -7.14 kcal/mol and 5.82 μM for ligand (H₂L) and -6.18 kcal/mol and 0.76 μM for complex [Cu(L)(phen)](1) with M^pro protein. This binding affinity is reasonably well as compared to recently known antiviral drugs. For instance, the binding affinity of ligand and complex was found to be better than docking results of chloroquine (-6.293 kcal/mol), hydroxychloroquine (-5.573 kcal/mol) and remdesivir (-6.352 kcal/mol) with M^pro protein. The present study may offer the technological solutions and potential inhibition to the COVID-19 virus in the ongoing and future challenges of the global community. In the framework of synthesis and characterization of mixed ligand copper (II) complex; the major conclusions can be drawn as follow.

Exploring Dyson's Orbitals and Their Electron Binding Energies for Conceptualizing Excited States from Response Methodology

The molecular orbital (MO) concept is a useful tool, which relates the molecular ground-state energy with the energies (and occupations) of the individual orbitals. However, analysis of the excited states from linear response computations is performed in terms of the initial state MOs or some other forms of orbitals, e.g., natural or natural transition orbitals. Because these orbitals lack the respective energies, they do not allow developing a consistent orbital picture of the excited states. Herein, we argue that Dyson's orbitals enable description of the response states compatible with the concepts of molecular orbital theory. The Dyson orbitals and their energies obtained by mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) for the response ground state are remarkably similar to the canonical MOs obtained by the usual DFT calculation. For excited states, the Dyson orbitals provide a chemically sensible picture of the electronic transitions, thus bridging the chasm between orbital theory and response computations.

Experimental measurement and prediction of ionic liquid ionisation energies

Ionic liquid (IL) valence electronic structure provides key descriptors for understanding and predicting IL properties. The ionisation energies of 60 ILs are measured and the most readily ionised valence state of each IL (the highest occupied molecular orbital, HOMO) is identified using a combination of X-ray photoelectron spectroscopy (XPS) and synchrotron resonant XPS. A structurally diverse range of cations and anions were studied. The cation gave rise to the HOMO for nine of the 60 ILs presented here, meaning it is energetically more favourable to remove an electron from the cation than the anion. The influence of the cation on the anion electronic structure (and vice versa) were established; the electrostatic effects are well understood and demonstrated to be consistently predictable. We used this knowledge to make predictions of both ionisation energy and HOMO identity for a further 516 ILs, providing a very valuable dataset for benchmarking electronic structure calculations and enabling the development of models linking experimental valence electronic structure descriptors to other IL properties, e.g. electrochemical stability. Furthermore, we provide design rules for the prediction of the electronic structure of ILs.

Characterization of Phytochemicals in Ulva intestinalis L. and Their Action Against SARS-CoV-2 Spike Glycoprotein Receptor-Binding Domain

Coronavirus disease-2019 (COVID-19) has caused a severe impact on almost all aspects of human life and economic development. Numerous studies are being conducted to find novel therapeutic strategies to overcome COVID-19 pandemic in a much effective way. Ulva intestinalis L. (Ui), a marine microalga, known for its antiviral property, was considered for this study to determine the antiviral efficacy against severe acute respiratory syndrome-associated Coronavirus-2 (SARS-CoV-2). The algal sample was dried and subjected to ethanolic extraction, followed by purification and analysis using gas chromatography-coupled mass spectrometry (GC-MS). Forty-three known compounds were identified and docked against the S1 receptor binding domain (RBD) of the spike (S) glycoprotein. The compounds that exhibited high binding affinity to the RBD of S1 protein were further analyzed for their chemical behaviour using conceptual density-functional theory (C-DFT). Finally, pharmacokinetic properties and drug-likeliness studies were carried out to test if the compounds qualified as potential leads. The results indicated that mainly phenols, polyenes, phytosteroids, and aliphatic compounds from the extract, such as 2,4-di-tert-butylphenol (2,4-DtBP), doconexent, 4,8,13-duvatriene-1,3-diol (DTD), retinoyl-β-glucuronide 6',3'-lactone (RBGUL), and retinal, showed better binding affinity to the target. Pharmacokinetic validation narrowed the list to 2,4-DtBP, retinal and RBGUL as the possible antiviral candidates that could inhibit the viral spike protein effectively.

Proposal of a Fermi-Dirac-Derived Reactivity Descriptor: Beyond the Frontier MO Model

In this paper, we derive a reactivity descriptor stemming from the Fermi-Dirac population scheme, applied to density functional calculations on molecular systems. Assuming that molecular orbitals only marginally change when temperature is slightly increased from 0 K, we study the response of electron density to a change in temperature. Connection with usual conceptual density functional theory descriptors is made, and the T-variation of electron density for some representative examples is given and discussed.

Molecular modelling of compounds used for corrosion inhibition studies: a review

Molecular modelling of organic compounds using computational software has emerged as a powerful approach for theoretical determination of the corrosion inhibition potential of organic compounds. Some of the common techniques involved in the theoretical studies of corrosion inhibition potential and mechanisms include density functional theory (DFT), molecular dynamics (MD) and Monte Carlo (MC) simulations, and artificial neural network (ANN) and quantitative structure-activity relationship (QSAR) modeling. Using computational modelling, the chemical reactivity and corrosion inhibition activities of organic compounds can be explained. The modelling can be regarded as a time-saving and eco-friendly approach for screening organic compounds for corrosion inhibition potential before their wet laboratory synthesis would be carried out. Another advantage of computational modelling is that molecular sites responsible for interactions with metallic surfaces (active sites or adsorption sites) and the orientation of organic compounds can be easily predicted. Using different theoretical descriptors/parameters, the inhibition effectiveness and nature of the metal-inhibitor interactions can also be predicted. The present review article is a collection of major advancements in the field of computational modelling for the design and testing of the corrosion inhibition effectiveness of organic corrosion inhibitors.