Which DFT Functional Performs Well in the Calculation of Methylcobalamin? Comparison of the B3LYP and BP86 Functionals and Evaluation of the Impact of Empirical Dispersion Correction

Dual External Field Strategy in Regulating the Superhalogen Characteristics of the Non-Noble Metal Constituted Tantalum Oxide Clusters

The identification of the non-noble metal constituted TaO cluster as a potential analogue to the noble metal Au is significant for the development of tailored materials. It leverages the superatom concept to engineer properties with precision. However, the impact of incrementally integrating TaO units on the electronic configurations and properties within larger TaO-based clusters remains to be elucidated. By employing the density functional theory calculations, the global minima and low-lying isomers of the TanOn (n = 2-5) clusters were determined, and their structural evolution was disclosed. In the cluster series, Ta5O5 was found to possess the highest electron affinity (EA) with a value of 2.14 eV, based on which a dual external field (DEF) strategy was applied to regulate the electronic property of the cluster. Initially, the electron-withdrawing CO ligand was affixed to Ta5O5, followed by the application of an oriented external electric field (OEEF). The CO ligation was found to be able to enhance the Ta5O5 cluster's electron capture capability by adjusting its electron energy levels, with the EA of Ta5O5(CO)4 peaking at 2.58 eV. Subsequently, the introduction of OEEF further elevated the EA of the CO-ligated cluster. Notably, OEEF, when applied along the +x axis, was observed to sharply increase the EA to 3.26 eV, meeting the criteria for superhalogens. The enhancement of EA in response to OEEF intensity can be quantified as a functional relationship. This finding highlights the advantage of OEEF over conventional methods, demonstrating its capacity for precise and continuous modulation of cluster EAs. Consequently, this research has adeptly transformed tantalum oxide clusters into superhalogen structures, underscoring the effectiveness of the DEF strategy in augmenting cluster EAs and its promise as a viable tool for the creation of superhalogens.

6-Bromo quinazoline derivatives as cytotoxic agents: design, synthesis, molecular docking and MD simulation

Based on unselectively, several side effects and drug resistance of available anticancer agents, the development and research for novel anticancer agents is necessary. In this study, a new series of quinazoline-4(3H)-one derivatives having a thiol group at position 2 of the quinazoline ring (8a-8 h) were designed and synthesized as potential anticancer agents. The Chemical structures of all compounds were characterized by 1H-NMR, 13C-NMR, and Mass spectroscopy. The antiproliferative activity of all derivatives were determined against two cancer cell lines (MCF-7 and SW480) and one normal cell lines (MRC-5) by the MTT method. Cisplatin, Erlotinib and Doxorubicin were used as positive controls. The results of in vitro screening showed that 8a with an aliphatic linker to SH group was the most potent compound with IC50 values of 15.85 ± 3.32 and 17.85 ± 0.92 µM against MCF-7 and SW480 cell lines, respectively. 8a indicated significantly better potency compared to Erlotinib in the MCF-7 cell line. The cytotoxic results obtained from testing compound 8a on the normal cell line, revealing an IC50 value of 84.20 ± 1.72 µM, provide compelling evidence of its selectivity in distinguishing between tumorigenic and non-tumorigenic cell lines. Structure-activity relationship indicated that the variation in the anticancer activities of quinazoline-4(3H)-one derivatives was affected by different substitutions on the SH position. Molecular docking and MD simulation were carried out for consideration of the binding affinity of compounds against EGFR and EGFR-mutated. The binding energy of compounds 8a and 8c were calculated at -6.7 and - 5.3 kcal.mol- 1, respectively. Compounds 8a and 8c were found to establish hydrogen bonds and some other important interactions with key residue. The DFT analysis was also performed at the B3LYP/6-31 + G(d, p) level for compounds 8a, 8c and Erlotinib. Compound 8a was thermodynamically more stable than 8c. Also, the calculated theoretical and experimental data for the IR spectrum were in agreement. The obtained results delineated that the 8a can be considered an appropriate pharmacophore to develop as an anti-proliferative agent.

Spodium bonding with noble gas atoms

The nature of the bonding between a neutral group 12 member (Zn3, Cd3 and Hg3) ring and a noble gas atom was explored using quantum chemical simulations. Natural bond orbital, quantum theory of atoms in molecules, symmetry-adapted perturbation theory, and molecular electrostatic potential surface analysis were also used to investigate the type of interaction between the noble gas atom and the metal rings (Zn3, Cd3 and Hg3). The Zn3, Cd3 and Hg3 rings are bonded to the noble gas through non-covalent interactions, which was revealed by the non-covalent interaction index. Additionally, energy decomposition analysis reveals that dispersion energy is the key factor in stabilizing these systems.

Simultaneous adsorption of ciprofloxacin drug and methyl violet dye on boron nitride nanosheets: experimental and theoretical insights

In this study, hexagonal boron nitride (BN) with a sheet-like morphology is successfully synthesized by reacting borax (Na2B4O7·10H2O) and urea (CO(NH2)2) powders in air via a facile microwave-assisted method within a short reaction time (15 min). The as-prepared product is structurally characterized via Fourier transformation infrared spectroscopy (FTIR), X-ray powder diffraction (XRD), Raman spectroscopy, field emission scanning electron microscopy (FESEM), energy dispersion X-ray analyzer (EDX), transmission electron microscopy (TEM), and Brunauer-Emmett-Teller (BET) surface area measurements. The adsorption process of methyl violet (MV) as a model of organic dyes and ciprofloxacin (CIP) as a model of antibiotics onto the boron nitride nanosheets has been experimentally and theoretically studied. The BN nanosheets exhibit the maximum adsorption capacity of 320.94 mg g-1 for MV dye and 266.29 mg g-1 for CIP antibiotic. The Freundlich isotherm model was suitable to describe the adsorption equilibrium isotherm data and the pseudo second-order model reflected the adsorption kinetics well. The calculated thermodynamic parameters show that the adsorption process is spontaneous under the measured conditions. The adsorption of CIP, MV and CIP + MV molecules on the surface of BN has been investigated through DFT calculations. The charge transfer and high adsorption capacity demonstrate the potential of BN nanosheets as an adsorbent for the simultaneous removal of MV dye and CIP drug from contaminated water.

Insight into the Mechanism Underlying Dehalococcoides mccartyi Strain CBDB1-Mediated B12-Dependent Aromatic Reductive Dehalogenation

Anaerobic bacteria transform aromatic halides through reductive dehalogenation. This dehalorespiration is catalyzed by the supernucleophilic coenzyme vitamin B₁₂, cob(I)alamin, in reductive dehalogenases. So far, the underlying inner-sphere electron transfer (ET) mechanism has been discussed controversially. In the present study, all 36 chloro-, bromo-, and fluorobenzenes and full-size cobalamin are analyzed at the quantum chemical density functional theory level with respect to a wide range of theoretically possible inner-sphere ET mechanisms. The calculated reaction free energies within the framework of Co^I···X (X = F, Cl, and Br) attack rule out most of the inner-sphere pathways. The only route with feasible energetics is a proton-coupled two-ET mechanism that involves a B₁₂ side-chain tyrosine (modeled by phenol) as a proton donor. For 12 chlorobenzenes and 9 bromobenzenes with experimental data from Dehalococcoides mccartyi strain CBDB1, the newly proposed PC-TET mechanism successfully discriminates 16 of 17 active from 4 inactive substrates and correctly predicts the observed regiospecificity to 100%. Moreover, fluorobenzenes are predicted to be recalcitrant in agreement with experimental findings. Conceptually, based on the Bell-Evans-Polanyi principle, the computational approach provides novel mechanistic insights and may serve as a tool for predicting the energetic feasibility of reductive aromatic dehalogenation.

The Perfluoro Cage Effect: A Search for Electron-Encapsulating Molecules

Quantum chemical calculations have for some time predicted that perfluorinated polyhedral organic molecules should exhibit a low-energy LUMO consisting of the overlapping inward-pointing lobes of the C-F σ* orbitals. Accordingly, these molecules should be able to encapsulate an electron within the interior of their cavities. Inspired by the recent confirmation of this prediction for perfluorocubane, we have sought to identify additional perfluorinated cage molecules capable of this remarkable behavior, which we refer to as the perfluoro cage effect (PCE). Using DFT calculations with multiple well-tested exchange-correlation functionals and large STO-QZ4P basis sets, we have identified several systems including [n]prismanes (n = 3-6), [n]asteranes (n = 3-5), twistane, and two norbornadiene dimer cages that clearly exhibit the PCE. In other words, they exhibit a low-energy LUMO belonging to the total symmetric irreducible representation of the point group in question and adiabatic electron affinities ranging from somewhat under 1 eV to over 2 eV. A pronounced size effect appears to hold, with larger cages exhibiting higher electron affinities (EAs). The largest adiabatic EAs, well over 3 eV, are predicted for perfluorinated dodecahedrane and C₆₀. In contrast, the PCE is barely discernible for perfluorinated tetrahedrane and bicyclo[1.1.1]pentane.

Beryllium bonding with noble gas atoms

Quantum chemical calculations were carried out to investigate the nature of the bonding between a neutral Be₃ ring and noble gas atom. Electronic structure calculation for these complexes was carried out at different computational levels in association with natural bond orbital, quantum theory of atoms in molecules, electron localization function, symmetry adapted perturbation theory, and molecular electrostatic potential surface analysis of Be₃ complexes. The Be atoms in the Be₃ moiety are chemically bonded to one another, with the BeBe bond dissociation energy being ~125 kJ mol^-1 . The Be₃ ring interacts with the noble gases through non-covalent interactions. The binding energies of the noble gas atoms with the Be₃ ring increases with increase in their atomic number. The non-covalent interaction index, density overlap region indicator and independent gradient model analyses reveal the presence of non-covalent inter-fragment interactions in the complexes. Energy decomposition analysis reveals that dispersion plays the major role towards stabilizing these systems.

A Theoretical Account of the Coupling between Metal- and Ligand-centred Spins

This study addresses the magnetic interaction between paramagnetic metal ions and the radical ligands taking the [Cu^II (hfac)₂ (imVDZ)] and [M^II (hfac)₂ (pyDTDA)] (imVDZ=1,5-dimethyl-3-(1-methyl-2-imidazolyl)-6-oxoverdazyl; hfac=(1,1,1,5,5,5)hexafluroacetylacetonate; pyDTDA=4-(2'-pyridyl)-1,2,3,5-dithiadiazolyl), (M=Cu, Ni, Co, Fe, Mn) compounds as reference systems. The coupling between the metal and ligand spins is quantified in terms of the exchange coupling constant (J) in the platform of density functional theory (DFT) and the wave function-based complete active space self-consistent field (CASSCF) method. Application of DFT and broken symmetry (BS) formalism results ferromagnetic coupling for all the transition metal complexes except the Mn(II) complex. This DFT-BS prediction of magnetic nature matches with the experimental finding for all the complexes other than the Fe(II)-pyDTDA complex, for which an antiferromagnetic coupling between high spin iron and the thiazyl ligand has been reported. However, evaluation of spin state energetics through the multiconfigurational wave function-based method produces the S=3/2 ground spin state for the iron-thiazyl in parity with experiment. Electronic structure analyses find the overlap between the metal- and ligand-based singly occupied molecular orbitals (SOMOs) to be one of the major reasons attributing to different extent of exchange coupling in the systems under investigation.

Porphyryne

Density functional theory calculations with the B3LYP*-D3 method with large STO-QZ4P basis sets unambiguously predict a singlet ground state for Zn-porphyryne. However, the calculations also predict a low singlet-triplet gap of about 0.4 eV and a high adiabatic electron affinity of 2.4 eV. Accordingly, the reactivity of porphyryne species may be dominated by electron transfer, hydrogen abstraction, and proton-coupled electron transfer processes.

Benchmark Density Functional Theory Approach for the Calculation of Bond Dissociation Energies of the M-O2 Bond: A Key Step in Water Splitting Reactions

A very fascinating aspect in quantum chemical research is to determine the accurate and cost-effective methods for the calculation of electronic and structural properties through a benchmark study. The current study focuses on the performance evaluation of density functional theory methods for the accurate measurement of bond dissociation energies (BDEs) of chemically important M-O₂ bonds in water splitting reactions. The BDE measurement has got noteworthy attention due to its importance in all areas of chemistry. For BDE measurements of M-O₂ bonds in five metal complexes with oxygen molecules, 14 density functionals (DFs) are chosen from seven classes of DFs with two series of mixed basis sets. A combination of pseudopotential and Pople basis sets [LANL2DZ & 6-31G(d) and SDD & 6-31+G(d)] are used as a series of mixed basis sets. The B3LYP-GD3BJ functional with LANL2DZ & 6-31G(d) gives outstanding results due to low deviations, error, and the best Pearson's correlation (R) between the experimental and theoretical data. Our study suggested an efficient, low-cost, precise, and accurate B3LYP-GD3BJ/LANL2DZ & 6-31G(d) level of theory for BDE of the M-O₂ bond, which may be useful for chemists working in the field of energy generation and utilization.

Fluoromethylsulfinyl radicals: spectroscopic characterization and photoisomerization via intramolecular hydrogen shift

Two new sulfinyl radicals, CHF₂SO˙ and CH₂FSO˙, have been generated in the gas phase through homolytic cleavage of the weak S-S bonds in disulfane oxides CHF₂S(O)SCF₃ and CH₂FS(O)SCF₃ by high-vacuum flash pyrolysis (HVFP) at ca. 500 °C. The IR spectroscopy characterization of the two fluoromethylsulfinyl radicals in solid N₂ (10 K), Ar (10 K), and Ne (3 K) matrices reveals the presence of two conformers for CHF₂SO˙ (gauche and cis) and one conformer for CH₂FSO˙ (gauche). Upon 266 nm laser irradiation, these radicals undergo both isomerization and decomposition in the matrices. In addition to the dominant formation of the elusive oxathiyl radicals CHF₂OS˙ (gauche and cis) and CH₂FOS˙ (gauche) via 1,2-alkyl migration, two higher-energy carbon-centered radicals ˙CF₂SOH and ˙CHFSOH bearing similar molecular structures to hydroperoxyalkyl radicals (˙QOOH) form via intramolecular 1,3-hydrogen shift in the two sulfinyl radicals. Additionally, the involvement of 1,3-hydrogen shift in CHF₂OS˙ and CH₂FOS˙ is also indicated by the observation of the fragmentation species. The identification of these radicals by matrix-isolation IR and UV-vis spectroscopy is aided by the quantum chemical calculations at the B3LYP/6-311++G(3df,3pd) level of theory. The stability of the isomers of the two sulfinyl radicals CHF₂SO˙ and CH₂FSO˙ has been discussed according to the experimental observations and also based on the CCSD(T)-F12a/aug-cc-pVTZ//B3LYP/6-311++G(3df,3pd) calculated energy profiles.

Fullerene binding effects in Al(III)/Zn(II) Porphyrin/Phthalocyanine photophysical properties and charge transport

We evaluate the fullerene C₆₀ binding effect; through the metal (Al) and through the ligand (Pc,TPP), on the photophysical and charge transport properties of M-porphyrin(TPP)/phthalocyanine(Pc) (M = Al(III), Zn(II)). We perform density functional theory (DFT) and time-dependent DFT calculations for the macrocycle-C₆₀ dyads, showing that all systems studied are thermodynamically favorable. The C₆₀ binding effect on the absorption spectrum is a red-shift of the Q and Soret (B) bands of TPPs and Pcs. The Pc-dyads show longer λ for Q bands (673 nm) than those with TPP (568 nm). AlTPP-C₆₀ and ZnTPP-C₆₀ show a more favorable electron injection to TiO₂ than the analogs Pcs, and the regeneration of the dye is preferred in AlTPP-C₆₀ and AlPc-C₆₀. Zero-bias conductance is computed (10^-4-10^-7 G₀) for the dyads using molecular junctions with Au(111)-based electrodes. When a bias voltage of around 0.6 V up to 1 V is applied, an increase in current is obtained for AlTPP-C₆₀ (10^-7 A), ZnTPP-C₆₀ (10^-7 A), and AlPc-C₆₀ (10^-8 A). Although there is not a unique trend in the behavior of the dyads, Pcs have better photophysical properties than TPPs and the latter are better in the charge transport. We conclude that AlTPP(ZnTPP)-C₆₀ dyads are an excellent alternative for designing new materials for dye-sensitized solar cells or optoelectronic devices.

Stabilization of Interstellar CSi2 Species by Donor Base Ligands: L-CSi2-L; L = cAACMe, NHCMe, and PMe3

The donor ligand bonded singlet (L)₂Si₂C containing a bent Si₂C unit in the middle has been studied by theoretical quantum mechanical calculations (NBO, QTAIM, EDA-NOCV analyses) [L = cAAC, NHC, Me₃P]. EDA-NOCV analysis suggests that this Si₂C is possible to stabilize by a pair of donor base ligands. The bond dissociation energy of the Si₂C fragment is endothermic (85-45 kcal/mol) with a sufficiently high intrinsic interaction energy (ΔE_int = -89 to -48 kcal/mol). Fifty percent of the total stabilization energy arises from electrostatic interactions, and nearly 45% is contributed by covalent orbital interaction between Si₂C and (L)₂ fragments in their singlet states. 75-80% of the orbital interaction energy is contributed by two sets of σ-donation L → SiCSi ← L. The π-back-donation is only 15-10%. The dispersion energy is not negligible (3-5%). The interaction energy is highest for 1 (L = cAAC) among three compounds. Additionally, (cAAC)₂Si₂C-Ni(CO)₃ (4) has been studied. The interaction energy between 1 and Ni(CO)₃ is nearly 61 kcal/mol with the major contribution coming from donation of electron cloud from electron rich Si₂C backbone to empty hybrid orbital of Ni(CO)₃ fragment. A sufficiently strong π-back-donation from (OC)₃Ni to Si₂C has also been identified.

Outer-sphere electron transfer does not underpin B12-dependent olefinic reductive dehalogenation in anaerobes

Anaerobic microbial B₁₂-dependent reductive dehalogenation may pave a way to remediate soil, sediment, and underground water contaminated with halogenated olefins. The chemical reaction is initiated by electron transfer (ET) from supernucleophilic cob(I)alamin (B_12s). However, the inherent mechanism as outer-sphere or inner-sphere route is still under debate. To clarify the possibility of an outer-sphere pathway, we calculated free energy barriers of the initial steps of all outer-sphere ET routes by Marcus theory employing density functional theory (DFT). For 18 fluorinated, chlorinated, and brominated ethenes as representative olefins, 164 of 165 reactions with free energy barriers larger than 20 kcal mol^-1 are not feasible under physiological dehalogenase conditions. Moreover, electronic structure analysis of perbromoethene with an outer-sphere free energy barrier of 18.2 kcal mol^-1 reveals no ET initiation down to Co⋯Br and Co⋯C distances of 3.15 Å. The results demonstrate that the B₁₂-catalyzed reductive dechlorination of olefins in microbes should proceed through an inner-sphere ET pathway.

Application of borophene as catechol sensor: a computational study

The efficacy of borophene (BP) as catechol (CC) sensor was explored using density functional theory (DFT) method. All calculations were performed at B3LYP level of theory and 6-31 + G(d) basis set employing the dispersion correction term of Grimme to consider dispersion interactions. The CC molecule is adsorbed on top of BP horizontally with the adsorption energy (E_ads) of about - 13.5 kcal·mol^-1. The HOMO and LUMO levels of nanosheet destabilize by about 0.36 and 0.14 eV, respectively, going from bare BP to BP-CC complex. Therefore, the E_g value decreases by about 10.5% upon adsorption process, which is a reasonable energy gap change for detection of CC. The negligible difference between the work function values (Φ, defined as the minimum amount of the energy needed to remove an electron from a solid to a point in the vacuum immediately outside the solid surface) of BP and its complex with CC indicates that the BP sheet is not an appropriate Φ-type sensor (in these sensors, adsorption of a chemical changes the gate voltage and produces an electrical signal that leads to the detection of chemical agent) for CC detection. The electrical conductivity of BP becomes 72 times higher after CC adsorption. The time needed for CC desorption from BP sheet is 7.6 ns, based on conventional transition state theory, showing that BP benefits from a short recovery time. The effect of CC concentration was explored by adsorption of 2 and 3 CC molecules on top of BP nanosheet and the results showed that the sensor response does not change by increasing the CC concentration. Also, the effect of lateral dimensions of BP on the adsorption energy was explored and it was shown that E_ads increases by enlargement of the nanosheet.

Surprising Homolytic Gas Phase Co-C Bond Dissociation Energies of Organometallic Aryl-Cobinamides Reveal Notable Non-Bonded Intramolecular Interactions

Aryl-cobalamins are a new class of organometallic structural mimics of vitamin B12 designed as potential 'antivitamins B12 '. Here, the first cationic aryl-cobinamides are described, which were synthesized using the newly developed diaryl-iodonium method. The aryl-cobinamides were obtained as pairs of organometallic coordination isomers, the stereo-structure of which was unambiguously assigned based on homo- and heteronuclear NMR spectra. The availability of isomers with axial attachment of the aryl group, either at the 'beta' or at the 'alpha' face of the cobalt-center allowed for an unprecedented comparison of the organometallic reactivity of such pairs. The homolytic gas-phase bond dissociation energies (BDEs) of the coordination-isomeric phenyl- and 4-ethylphenyl-cobinamides were determined by ESI-MS threshold CID experiments, furnishing (Co-Csp 2 )-BDEs of 38.4 and 40.6 kcal mol-1 , respectively, for the two β-isomers, and the larger BDEs of 46.6 and 43.8 kcal mol-1 for the corresponding α-isomers. Surprisingly, the observed (Co-Csp 2 )-BDEs of the Coβ -aryl-cobinamides were smaller than the (Co-Csp 3 )-BDE of Coβ -methyl-cobinamide. DFT studies and the magnitudes of the experimental (Co-Csp 2 )-BDEs revealed relevant contributions of non-bonded interactions in aryl-cobinamides, notably steric strain between the aryl and the cobalt-corrin moieties and non-bonded interactions with and among the peripheral sidechains.

Unsupported Donor-Acceptor Complexes of Noble Gases with Group 13 Elements

Unsupported donor-acceptor complexes of noble gases (Ng) with group 13 elements have been theoretically studied using density functional theory. Calculations reveal that heavier noble gases form thermodynamically stable compounds. The present study reveals that no rigid framework is necessary to stabilize the donor-acceptor complexes. Rather, prepyramidalization at the Lewis acid center may be an interesting alternative to stabilize these complexes. Detailed bonding analyses reveal the formation of two-center-two-electron dative bonding, where Ng atoms act as a donor.

The bimetallic and the anchoring group effects on both optical and charge transport properties of hexaphyrin amethyrin

Bimetallic Cu(ii)-hexaphyrin amethyrin proposed as a molecular switch operated by the application of an external magnetic field.

Dehalococcoides-Mediated B12-Dependent Reductive Dehalogenation of Aromatics Does Not Proceed through Outer-Sphere Electron Transfer

Several anaerobic bacteria can couple the reduction of aromatic halides to energy conservation. This organohalide respiration is catalyzed by enzymes containing cob(I)alamin, an activated supernucleophilic form of the coenzyme vitamin B₁₂. However, the mechanism underlying the electron transfer (inner-sphere vs outer-sphere ET) still remains elusive. To clarify this issue, we selected 36 fluoro-, chloro-, and bromobenzenes as representative substrates and calculated their free-energy barriers at the quantum chemical density functional theory level, considering a wide range of theoretically possible outer-sphere ET mechanisms. Across all 336 reaction routes addressed, 334 routes involve free-energy barriers larger than 20 kcal/mol. For two reaction routes with highly brominated benzenes, free-energy barriers below 20 kcal/mol imply abiotic reduction as observed in experiments. Thus, microbial B₁₂-dependent aromatic reductive dehalogenation does not proceed through an outer-sphere ET mechanism. Instead, the present study strongly suggests that microbe-catalyzed reductive dehalogenation of aromatic halides is governed by inner-sphere ET.

Fluorescence properties and quantum-chemical modeling of tert-butyl-substituted porphyrazines: Structural and ionization effect

Synthesis and identification of tetrakis-[5,6-bis(4-tert-butylphenyl)pyrazino] porphyrazine, tetra-(4-tert-butyl)phthalocyanine and octakis-(4-tert-butylphenyl)porphyrazine were carried out. Spectrophotometric method was used to study the spectral, acidic and fluorescence properties of the synthesized compounds. It was determined that the synthesized tert-butyl-substituted porphyrazines exhibit a high sensitivity of fluorescence to the molecule ionization. To understand the features of the spectral properties the geometry optimization and an analysis of energy levels and localization of highest occupied and lowest unoccupied molecular orbitals of the studied compounds were performed on the basis of density functional theory with the BP86 functional and the def2-TZVP basis set. The effect of substituents in molecular fragments of the macrocycle on the acidic and electro-optical properties of the studied compounds is revealed. Materials with pH-tunable fluorescence were designed.

Spectral-Fluorescence Properties of Zn(II)-Octaphenyltetraazaporphyrins

Zn(II)-octa-(4-chlorophenyl)- and Zn(II)-octa-(4-bromophenyl)tetraazaporphyrins were synthesized by the reaction of cyclotetramerization of di-(4-chlorophenyl)- and di-(4-bromophenyl)maleonitriles with zinc(II) chloride. The obtained compounds were identified by UV-vis, IR, NMR ¹H spectroscopy and mass spectrometry. Geometry optimization of the series of halogenated Zn(II)-octaaryltetraazaporphyrins was performed using the density functional method with the BP86 functional and the def2-TZVP basis set. An analysis of the distribution of molecular orbital energies in the neighborhood of highest occupied molecular orbitals (HOMO and HOMO-1) and lowest unoccupied molecular orbitals (LUMO and LUMO+1) and the width of the HOMO - LUMO energy gaps (E_H-L) was performed for the studied compounds. Fluorimetric measurements of the Zn(II)-octaphenyltetraazaporphyrins in toluene were carried out and fluorescence quantum yields of studied compounds were determined and analyzed. It has been shown that the halogen on the para-position of the phenyl groups significantly affects the value of the obtained quantum yields of fluorescence emission but does not significantly affect the Stokes shifts.

A reaction pathway to compound 0 intermediates in oxy-myoglobin through interactions with hydrogen sulfide and His64

Myoglobin (Mb) binds oxygen with high affinity as a low spin singlet complex and thus functions as an oxygen storage protein. Yet, hybrid Density Functional Theory/Molecular Mechanical (DFT/MM) calculations of oxy-Mb models predict that the O₂ bond is much less resistant to breaking in the presence of hydrogen sulfide (H₂S) compared with water. Specifically, a hydrogen atom from H₂S can be transferred to the distal oxygen atom through homolytic cleavage of the S-H bond to form the intermediate Compound (Cpd) 0 structure and a thiyl radical. In the presence of a neutral His64 (Nε protonation, His64-ε) and H₂S, only a metastable Cpd 0 would be formed as the active site is devoid of any additional proton donor to fully break the O₂ bond. In contrast, the calculations predict that the triplet state is significantly favored over the open shell singlet diradical state throughout the entire reaction coordinate in the presence of H₂S and a positively charged His64. Furthermore, a positively charged His64 can readily donate a proton to Cpd 0 to fully break the O₂ bond resulting in a configuration analogous to reported reaction models of a hemoglobin mutant bound to H₂O₂ with H₂S present. Typically, exotic techniques are required to generate Cpd 0 but under the conditions just described the intermediate is readily detected in UV-Vis spectra at room temperature. The effect is observed as a 2 nm red shift of the Soret band from 414 nm to 416 nm (pH 5.0, His64-εδ) and from 416 nm to 418 nm (pH 6.6, His64-ε).

Methodological Survey of Simplified TD-DFT Methods for Fast and Accurate Interpretation of UV-Vis-NIR Spectra of Phthalocyanines

A methodological survey of density functional theory (DFT) methods for the prediction of UV-visible (vis)-near-infrared (NIR) spectra of phthalocyanines is reported. Four methods, namely, full time-dependent (TD)-DFT and its Tamm-Dancoff approximation (TDA), together with their simplified modifications (sTD-DFT and sTDA, respectively), were tested by using the examples of unsubstituted and alkoxy-substituted metal-free ligands and zinc complexes. The theoretical results were compared with experimental data derived from UV-visible absorption and magnetic circular dichroism spectroscopy. Seven popular exchange-correlation functionals (BP86, B3LYP, TPSSh, M06, CAM-B3LYP, LC-BLYP, and ωB97X) were tested within these four approaches starting at a relatively modest level using 6-31G(d) basis sets and gas-phase BP86/def2-SVP optimized geometries. A gradual augmentation of the computational levels was used to identify the influence of starting geometry, solvation effects, and basis sets on the results of TD-DFT and sTD-DFT calculations. It was found that although these factors do influence the predicted energies of the vertical excitations, they do not affect the trends predicted in the spectral properties across series of structurally related substituted free bases and metallophthalocyanines. The best accuracy for the gas-phase vertical excitations was observed in the lower-energy Q-band region for calculations that made use of range-separated hybrids for both full and simplified TD-DFT approaches. The CAM-B3LYP functional provided particularly accurate results in the context of the sTD-DFT approach. The description of the higher-energy B-band region is considerably less accurate, and this demonstrates the need for further advances in the accuracy of theoretical calculations. Together with a general increase in accuracy, the application of simplified TD-DFT methods affords a 2-3 orders of magnitude speedup of the calculations in comparison to the full TD-DFT approach. It is anticipated that this approach will be widely used on desktop computers during the interpretation of UV-vis-NIR spectra of phthalocyanines and related macrocycles in the years ahead.

Exploring the Structure, Energetic, and Magnetic Properties of Neutral Small Lithium Clusters Doped with Yttrium: Supermagnetic Atom Research

Density functional theory calculations based on magnetic and energetic stability criteria were performed to study a series of yttrium-doped lithium neutral clusters. A relativistic approximation was employed to properly describe the energy and multiplicity of the given clusters' fundamental states. The interaction of the 4d-Y atomic orbitals with the sp-Li states had an important role in the magnetic and energetic behavior of the selected systems. The spin density was concentrated over the yttrium atom regardless of the size of the cluster. Li7Y is a new stable superatom due to its enhanced magnetic properties.

Revisiting the catalytic mechanism of Mo-Cu carbon monoxide dehydrogenase using QM/MM and DFT calculations

Previous density functional theory (DFT) studies have shown that the release of the produced carbon dioxide (CO2) from an active-site cluster is a thermodynamically or kinetically difficult step in the enzymatic carbon monoxide (CO) oxidation catalyzed by Mo-Cu carbon monoxide dehydrogenase (Mo-Cu CODH). To better understand the effect of the protein environment on this difficult CO2 release step as well as other reaction steps, we applied hybrid quantum mechanics and molecular mechanics (QM/MM) calculations to the Mo-Cu CODH enzyme. The results show that in the first step, the equatorial Mo[double bond, length as m-dash]O group in the active-site cluster attacks the nearby CO molecule bound to the Cu site. Afterward, a stable thiocarbonate intermediate is formed in which the CO2 molecule is embedded and the copper-S(μ-sulfido) bond is broken. A free CO2 molecule, i.e., the final product, is then released from the active-site cluster, not directly from the thiocarbonate intermediate but via a previously formed intermediate that also contains CO2 but retains the Cu-S(μ-sulfido) bond. In contrast to the previous DFT results, the calculated barrier for this process was low in our QM/MM calculations. An additional QM/MM analysis of the barrier height showed that the effect of the protein environment on this barrier lowering is not very large. We found that the reason for the low barrier obtained by QM/MM is that the barrier for CO2 release is already not high at the DFT level. These results allow us to conclude that the CO oxidation reaction passes through the formation of a thiocarbonate intermediate, and that the subsequent CO2 release is kinetically not difficult. Nevertheless, the protein environment has an important role to play in making the latter process thermodynamically favored. No low-barrier pathway for the product release could be obtained for the reaction of n-butylisocyanide, which is consistent with the experimental fact that n-butylisocyanide inhibits Mo-Cu CODH.

Computational Study on N-N Homolytic Bond Dissociation Enthalpies of Hydrazine Derivatives

The hydrazine derivatives have been regarded as the important building blocks in organic chemistry for the synthesis of organic N-containing compounds. It is important to understand the structure-activity relationship of the thermodynamics of N-N bonds, in particular, their strength as measured by using the homolytic bond dissociation enthalpies (BDEs). We calculated the N-N BDEs of 13 organonitrogen compounds by eight composite high-level ab initio methods including G3, G3B3, G4, G4MP2, CBS-QB3, ROCBS-QB3, CBS-Q, and CBS-APNO. Then 25 density functional theory (DFT) methods were selected for calculating the N-N BDEs of 58 organonitrogen compounds. The M05-2X method can provide the most accurate results with the smallest root-mean-square error (RMSE) of 8.9 kJ/mol. Subsequently, the N-N BDE predictions of different hydrazine derivatives including cycloalkylhydrazines, N-heterocyclic hydrazines, arylhydrazines, and hydrazides as well as the substituent effects were investigated in detail by using the M05-2X method. In addition, the analysis including the natural bond orbital (NBO) as well as the energies of frontier orbitals were performed in order to further understand the essence of the N-N BDE change patterns.

The Binding Mode of the Sonic Hedgehog Inhibitor Robotnikinin, a Combined Docking and QM/MM MD Study

Erroneous activation of the Hedgehog pathway has been linked to a great amount of cancerous diseases and therefore a large number of studies aiming at its inhibition have been carried out. One leverage point for novel therapeutic strategies targeting the proteins involved, is the prevention of complex formation between the extracellular signaling protein Sonic Hedgehog and the transmembrane protein Patched 1. In 2009 robotnikinin, a small molecule capable of binding to and inhibiting the activity of Sonic Hedgehog has been identified, however in the absence of X-ray structures of the Sonic Hedgehog-robotnikinin complex, the binding mode of this inhibitor remains unknown. In order to aid with the identification of novel Sonic Hedgehog inhibitors, the presented investigation elucidates the binding mode of robotnikinin by performing an extensive docking study, including subsequent molecular mechanical as well as quantum mechanical/molecular mechanical molecular dynamics simulations. The attained configurations enabled the identification of a number of key protein-ligand interactions, aiding complex formation and providing stabilizing contributions to the binding of the ligand. The predicted structure of the Sonic Hedgehog-robotnikinin complex is provided via a PDB file as Supplementary Material and can be used for further reference.

Co(iii) complexes of (1,3-selenazol-2-yl)hydrazones and their sulphur analogues

The first Co(iii) complexes with (1,3-selenazol-2-yl)hydrazones as an unexplored class of ligands were prepared and characterized by NMR spectroscopy and X-ray diffraction analysis. The novel ligands act as NNN tridentate chelators forming octahedral Co(iii) complexes. The impact of structural changes on ligands' periphery as well as that of isosteric replacement of sulphur with selenium on the electrochemical and electronic absorption features of complexes are explored. To support the experimental data, density functional theory (DFT) calculations were also conducted. Theoretical NMR chemical shifts, the relative energies and natural bond orbital (NBO) analysis are calculated within the DFT approach, while the singlet excited state energies and HOMO-LUMO energy gap were calculated with time-dependent density functional theory (TD-DFT). The electrophilic f^- and nucleophilic f⁺ Fukui functions are well adapted to find the electrophile and nucleophile centres in the molecules. Both (1,3-selenazol-2-yl)- and (1,3-thiazol-2-yl)hydrazone Co(iii) complexes showed potent antimicrobial and antioxidant activity. A significant difference among them was a smaller cytotoxicity of selenium compounds.

Computational study on C–B homolytic bond dissociation enthalpies of organoboron compounds

The theoretical study of three hybridized C–B BDEs with different substituents can provide corresponding guidance to experimental research studies.