Electronic structures, bonding and energetics of non-heme mono and dinuclear iron-TPA complexes: a computational exploration

DFT and TDDFT exploration on electronic transitions and bonding aspect of DPA and PTDC ligated transition metal complexes

CONTEXT

In this study, we have investigated the structure, reactivity, bonding, and electronic transitions of DPA and PDTC along with their Ni-Zn complexes using DFT/TD-DFT methods. The energy gap between the frontier orbitals was computed to understand the reactivity pattern of the ligands and metal complexes. From the energies of FMO's, the global reactivity descriptors such as electron affinity, ionization potential, hardness (η), softness (S), chemical potential (μ), electronegativity (χ), and electrophilicity index (ω) have been calculated. The complexes show a strong NLO properties due to easily polarization as indicated by the narrow HOMO-LUMO gap. The polarizability and hyperpolarizabilities of the complexes indicate that they are good candidates for NLO materials. Molecular electrostatic potential (MEP) maps identified electrophilic and nucleophilic sites on the surfaces of the complexes. TDDFT and NBO analyses provided insights into electronic transitions, bonding, and stabilizing interactions within the studied complexes. DPA and PDTC exhibited larger HOMO-LUMO gaps and more negative electrostatic potentials compared to their metal complexes suggesting the higher reactivity. Ligands (DPA and PDTC) had absorption spectra in the range of 250 nm to 285 nm while their complexes spanned 250 nm to 870 nm. These bands offer valuable information on electronic transitions, charge transfer and optical behavior. This work enhances our understanding of the electronic structure and optical properties of these complexes.

METHODS

Gaussian16 program was used for the optimization of all the compounds. B3LYP functional in combination with basis sets, such as LanL2DZ for Zn, Ni and Cu while 6-311G** for other atoms like C, H, O, N, and S was used. Natural bond orbital (NBO) analysis is carried out to find out how the filled orbital of one sub-system interacts with the empty orbital of another sub-system. The ORCA software is used for computing spectral features along with the zeroth order regular approximation method (ZORA) to observe its relativistic effects. TD-DFT study is carried out to calculate the excitation energy by using B3LYP functional.

DFT and TDDFT exploration on the role of pyridyl ligands with copper toward bonding aspects and light harvesting

CONTEXT

Schiff base-containing metal complexes have been the subject of extensive research. In this work, a coordination polymer-derived complex called [Cu(L)] that is solution-stable (L = 2-(2-hydroxybenzylidene-amino)phenol) has been explored theoretically with five different pyridyl-based ligands using DFT/TDDFT in order to understand the structural-functional and electronic transitions of these five complexes. Frontier molecular orbital (FMO) analysis was carried out to assess the reactivity behavior of all five complexes. For the purpose of studying the charge energy distribution over complexes, electrostatic potential maps were also drawn. Furthermore, in order to identify any stabilizing interactions that may be present in the given complexes, an NBO analysis was studied. To learn more about any potential correlations between the properties of these five complexes, a comparative analysis was explored. Our calculations demonstrate that complex 3 having pyridine-4-carboxamide as a ligand has a lower energy gap and a higher negative electrostatic potential which may indicate its higher reactivity and this may be due to the electron withdrawing group (carboxamide). TDDFT results show that the highest light harvesting efficiency (LHE) of all the studied complexes is found in the range of 440-448 nm. Complexes 1, 2, and 4 show the higher light harvesting efficiency as compared to complexes 3 and 5. Our findings are in good accordance with the available experimental data.

METHODS

All DFT computations were performed using the Gaussian16 with unrestricted B3LYP-D2 functional with the basis sets 6-31G(d,p) for O, N, C, and H while LanL2DZ for Cu. The polarized continuum model (PCM) was used for the solvation. The software GaussView6.1 was utilized for the modeling of initial geometries and the plotting of MEP maps. The NBO6.0 program which is incorporated in Gaussian16 was utilized to investigate the bonding nature and stabilization energies of the complexes. The ORCA program was used to simulate the absorption spectra.

Electronic structures, bonding aspects and spectroscopic parameters of homo/hetero valent bridged dinuclear transition metal complexes

Bridged dinuclear metal complexes have fascinated scientists worldwide, and remarkable success has been achieved to unravel the electronic structures, structure-function relationship, coordination environments, and fine mechanistic details of the enzymes owing to the repercussion of biomimetic studies carried out on dinuclear model systems. Molecular level study of these systems integrated with spectroscopic study helps in gaining deep insights about structural and electronic aspects of natural enzymatic systems. Considering the same, here first time we report DFT study on bridged non-heme metal complexes based on N-Et-HPTB ligand system containing homovalent (M^IIM^II); {[(Mn^II)₂(O₂CCH₃)(N-Et-HPTB)]²⁺; Species I), [(Fe^II)₂(O₂CCH₃)(N-Et-HPTB)]²⁺; Species II), [(Co^II)₂(O₂CCH₃)(N-Et-HPTB)]²⁺; Species III)} and heterovalent (M^IIIM^II): {[(Mn^III)(Mn^II)(O₂)(N-Et-HPTB)]²⁺; Species Ia) [(Fe^III)(Fe^II)(O₂)(N-Et-HPTB)]²⁺; Species IIa) and [(Co^III)(Co^II)(O₂)(N-Et-HPTB)]²⁺; Species IIIa)} dinuclear metal centres. Bridging oxygen bears a significant spin density which may prompt important chemical reactions involving activation of bonds like C-H/O-H/N-H etc. TD-DFT calculations for UV-Visible absorption have been carried out to further shed light on structural-functional and electronic structures of these dinuclear species. Studying these dinuclear species may be a good starting point for the study of active sites of the bimetallic centre of dinuclear enzymes and thus may serve as fascinating spectroscopic models. Further, FMO analysis, MEP mapping, and NBO calculations were employed to analyze bonding aspects predict theoretical reactivity behaviour and any kind of stabilizing interactions present in the reported species.

Theoretical exploration on structures, bonding aspects and molecular docking of α-aminophosphonate ligated copper complexes against SARS-CoV-2 proteases

Recent years have witnessed a growing interest in the biological activity of metal complexes of α-aminophosphonates. Here for the first time, a detailed DFT study on five α-aminophosphonate ligated mononuclear/dinuclear CuII complexes is reported using the dispersion corrected density functional (B3LYP-D2) method. The electronic structures spin densities, FMO analysis, energetic description of spin states, and theoretical reactivity behaviour using molecular electrostatic potential (MEP) maps of all five species are reported. All possible spin states of the dinuclear species were computed and their ground state S values were determined along with the computation of their magnetic coupling constants. NBO analysis was also performed to provide details on stabilization energies. A molecular docking study was performed for the five complexes against two SARS-CoV-2 coronavirus protein targets (PDB ID: 6LU7 and 7T9K). The docking results indicated that the mononuclear species had a higher binding affinity for the targets compared to the dinuclear species. Among the species investigated, species I showed the highest binding affinity with the SARS-CoV-2 Omicron protease. NPA charge analysis showed that the heteroatoms of model species III had a more nucleophilic nature. A comparative study was performed to observe any variations and/or correlations in properties among all species.

Computational studies on potential new anti-Covid-19 agents with a multi-target mode of action

A compound that could inhibit multiple targets associated with SARS-CoV-2 infection would prove to be a drug of choice against the virus. Human receptor-ACE2, receptor binding domain (RBD) of SARS-CoV-2 S-protein, Papain-like protein of SARS-CoV-2 (PLpro), reverse transcriptase of SARS-CoV-2 (RdRp) were chosen for in silico study. A set of previously synthesized compounds (1-5) were docked into the active sites of the targets. Based on the docking score, ligand efficiency, binding free energy, and dissociation constants for a definite conformational position of the ligand, inhibitory potentials of the compounds were measured. The stability of the protein-ligand (P-L) complex was validated in silico by using molecular dynamics simulations using the YASARA suit. Moreover, the pharmacokinetic properties, FMO and NBO analysis were performed for ranking the potentiality of the compounds as drug. The geometry optimizations and electronic structures were investigated using DFT. As per the study, compound-5 has the best binding affinity against all four targets. Moreover, compounds 1, 3 and 5 are less toxic and can be considered for oral consumption.