Conceptual and computational DFT in the study of aromaticity

Density functional theory and enzyme studies support interactions between angiotensin receptor blockers and angiotensin converting enzyme-2: Relevance to coronavirus 2019

The binding affinities and interactions between eight drug candidates, both commercially available (candesartan; losartan; losartan carboxylic acid; nirmatrelvir; telmisartan) and newly synthesized benzimidazole-N-biphenyltetrazole (ACC519T), benzimidazole bis-N,N'-biphenyltetrazole (ACC519T(2) and 4-butyl-N,N-bis([2-(2H-tetrazol-5-yl)biphenyl-4-yl]) methyl (BV6), and the active site of angiotensin-converting enzyme-2 (ACE2) were evaluated for their potential as inhibitors against SARS-CoV-2 and regulators of ACE2 function through Density Functional Theory methodology and enzyme activity assays, respectively. Notably, telmisartan and ACC519T(2) exhibited pronounced binding affinities, forming strong interactions with ACE2's active center, favorably accepting proton from the guanidinium group of arginine273. The ordering of candidates by binding affinity and reactivity descriptors, emerged as telmisartan > ACC519T(2) > candesartan > ACC519T > losartan carboxylic acid > BV6 > losartan > nirmatrelvir. Proton transfers among the active center amino acids revealed their interconnectedness, highlighting a chain-like proton transfer involving tyrosine, phenylalanine, and histidine. Furthermore, these candidates revealed their potential antiviral abilities by influencing proton transfer within the ACE2 active site. Furthermore, through an in vitro pharmacological assays we determined that candesartan and the BV6 derivative, 4-butyl-N,N0-bis[20-2Htetrazol-5-yl)bipheyl-4-yl]methyl)imidazolium bromide (BV6(K+)2) also contain the capacity to increase ACE2 functional activity. This comprehensive analysis collectively underscores the promise of these compounds as potential therapeutic agents against SARS-CoV-2 by targeting crucial protein interactions.

Aromaticity of six-membered nitro energetic compounds through molecular electrostatic potential, magnetic, electronic delocalization and reactivity-based indices

The electron density depletion associated with π-hole at the ring center typical of energetic compounds was clearly revealed by the molecular electrostatic potential (ESP). In addition, the spatial arrangement of NO2 groups appears to affect the ESP value in the ring center, and therefore sensitivity to detonation. Indeed, for monocyclic nitrobenzene compounds with the same number of NO2 groups, the ESP value in the ring center decreases as the NO2 groups are more closely spaced. As expected, the central rings become less aromatic as NO2 groups are added. The MCI, PDI, PLR, NICSzz(1), FLU indices are all strongly correlated with the ESP values observed in the ring center of the set of nitrobenzenes. Aromaticity indices based on electron delocalization criteria appear to be very sensitive to small variations in aromaticity. Among magnetic-based indices, only NICSzz(1) is capable to predict small changes in aromaticity. The PLR index derived from conceptual DFT is quite relevant for predicting small variations in aromaticity. According to our results, the most suitable aromaticity index is not based on a single criterion, and that selecting it is more subtle. Therefore, it is important to combine information from several criteria to obtain a more complete description of the aromaticity of the studied compounds.

Review and Theoretical Analysis of Fluorinated Radicals in Direct CAr-H Functionalization of (Hetero)arenes

We highlight key contributions in the field of direct radical CAr- H (hetero)aromatic functionalization involving fluorinated radicals. A compilation of Functional Group Transfer Reagents and their diverse activation mechanisms leading to the release of radicals are discussed. The substrate scope for each radical is analyzed and classified into three categories according to the electronic properties of the substrates. Density functional theory computational analysis provides insights into the chemical reactivity of several fluorinated radicals through their electrophilicity and nucleophilicity parameters. Theoretical analysis of their reduction potentials also highlights the remarkable correlation between electrophilicity and oxidizing ability. It is also established that highly fluorinated radicals (e.g. ⋅OCF3) are capable of engaging in single-electron transfer (SET) processes rather than radical addition, which is in good agreement with experimental literature data. A reactivity scale, based on activation barrier of addition of these radicals to benzene is also elaborated using the high accuracy DLPNO-(U)CCSD(T) method.

Advances for Triangular and Sandwich-Shaped All-Metal Aromatics

Much experimental work has been contributed to all-metal σ, π and δ-aromaticity among transition metals, semimetallics and other metals in the past two decades. Before our focused investigations on the properties of triangular and sandwich-shaped all-metal aromatics, A. I. Boldyrev presented general discussions on the concepts of all-metal σ-aromaticity and σ-antiaromaticity for metallo-clusters. Schleyer illustrated that Nucleus-Independent Chemical Shifts (NICS) were among the most authoritative criteria for aromaticity. Ugalde discussed the earlier developments of all-metal aromatic compounds with all possible shapes. Besides the theoretical predictions, many stable all-metal aromatic trinuclear clusters have been isolated as the metallic analogues of either the σ-aromatic molecule's [H3]+ ion or the π-aromatic molecule's [C3H3]+ ion. Different from Hoffman's opinion on all-metal aromaticity, triangular all-metal aromatics were found to hold great potential in applications in coordination chemistry, catalysis, and material science. Triangular all-metal aromatics, which were theoretically proved to conform to the Hückel (4n + 2) rule and possess the smallest aromatic ring, could also play roles as stable ligands during the formation of all-metal sandwiches. The triangular and sandwich-shaped all-metal aromatics have not yet been specifically summarized despite their diversity of existence, puissant developments and various interesting applications. These findings are different from the public opinion that all-metal aromatics would be limited to further applications due to their overstated difficulties in synthesis and uncertain stabilities. Our review will specifically focus on the summarization of theoretical predictions, feasible syntheses and isolations, and multiple applications of triangular and sandwich shaped all-metal aromatics. The appropriateness and necessities of this review will emphasize and disseminate their importance and applications forcefully and in a timely manner.

Theoretical investigation of substitution effect on the sixth and seventh positions of coumarin derivatives

The linear and non-linear optical properties of 6-donor and 7-donor substituted coumarins were compared using density functional theory (DFT) and time-dependent-DFT (TD-DFT). Charge transfer characteristics were investigated through natural bond order analysis, frontier molecular orbital, and molecular electrostatic potential plots. TD-DFT results suggested that the 6-donor substituted coumarins (PS1, PS3, and PS5) showed red-shifted absorption than the 7-donor substituted coumarins (PS2, PS4, and PS6). The chemical potential (μ) and electrophilicity index (ω) showed direct relation with the band gap and an inverse relation with chemical hardness (η) and hyperhardness (Γ). The global reactivity descriptors μ and ω showed direct and η and Γ showed an inverse correlation with first-order hyperpolarizability (β0) and second-order hyperpolarizability (γ). The β0 and γ for 7-donor substituted coumarin are higher than for 6-donor substituted coumarin.

In Silico and In Vitro Analyses to Repurpose Quercetin as a Human Pancreatic α-Amylase Inhibitor

Human pancreatic α-amylase (HPA), situated at the apex of the starch digestion hierarchy, is an attractive therapeutic approach to precisely regulate blood glucose levels, thereby efficiently managing diabetes. Polyphenols offer a natural and multifaceted approach to moderate postprandial sugar spikes, with their slight modulation in carbohydrate digestion and potential secondary benefits, such as antioxidant and anti-inflammatory effects. Taking into consideration the unfavorable side effects of currently available commercial medications, we aimed to study a library of polyphenols attributed to their remarkable antidiabetic properties and screened the most potent HPA inhibitor via a comprehensive in silico study encompassing molecular docking, molecular mechanics with generalized Born and surface area solvation (MM/GBSA) calculation, molecular dynamics (MD) simulation, density functional theory (DFT) study, and pharmacokinetic properties followed by an in vitro assay. Significant hydrogen bonding with the catalytic triad residues of HPA, prominent MM/GBSA binding energy of -27.03 kcal/mol, and the stable nature of the protein-ligand complex with regard to 100 ns MD simulation screened quercetin as the best HPA inhibitor. Additionally, quercetin showed strong reactivity in the substrate-binding pocket of HPA and exhibited favorable pharmacokinetic properties with a considerable inhibitory concentration (IC50) of 57.37 ± 0.9 μg/mL against α-amylase. This study holds prospects for HPA inhibition and suggests quercetin as an approach to therapy for diabetes; however, it is imperative to conduct further research.

Cyclo[n]carbons and catenanes from different perspectives: disentangling the molecular thread

All-carbon atomic rings, cyclo[n]carbons, have recently attracted vivid attention of experimentalists and theoreticians. Among them, cyclo[18]carbon is the most studied system. In this paper, we summarize and review various properties of cyclo[n]carbons, emphasising the aspects of their aromaticity/antiaromaticity. In the first part, the trends in bonding patterns and selected aromaticity indices with the increasing size of the rings are discussed. In the second part we explore the properties of catenane models based on interlocked cyclo[18]carbon rings from different perspectives and investigate their behaviour under the action of external force using computational experiments.

π-Aromaticity Dominating in a Saturated Ring: Neutral Aromatic Silicon Analogues of Cyclobutane-1,3-diyls

The synthesis, structures, and reactivity of the first neutral 2π-aromatic Si4 rings [LSiSiAr(X)]2 (3: X = Br; 4: X = Cl; L = PhC(NtBu)2, Ar = 2,4,6-Me3C6H2) were described. Compounds 3 and 4 were obtained by 1,3-halogenation of tetrasilacyclobutadiene (LSiSiAr)2 (2), which was prepared by the reductive cross-coupling of trisilane (ArSiCl2)2SiHAr with two equiv of chlorosilylene LSiCl. The reaction of 3 with two equiv of PhLi yielded the corresponding substitution Si4 ring [LSiSiAr(Ph)]2 (5). Single-crystal X-ray diffraction analysis of 3 disclosed that it adopts both puckered (3a) and planar (3b) structures in the solid state, whereas 4 and 5 exhibit only a puckered structure. DFT calculations suggested that the puckered 3a features almost the same electronic structure with fully delocalized 2π planar 3b. The dominant 2π-aromaticity of 3 in a σ-frame has been demonstrated by DFT calculations, providing the first example of aromatics featuring both planar and puckered structures.

Aromaticity: Quo Vadis

Aromaticity is one of the most deeply rooted concepts in chemistry. But why, if two-thirds of existing compounds can be classified as aromatic, is there no consensus on what aromaticity is? σ-, π-, δ-, spherical, Möbius, or all-metal aromaticity… why are so many attributes needed to specify a property? Is aromaticity a dubious concept? This perspective aims to reflect where the aromaticity community is and where it is going.

Single-molecule conductance studies on quasi- and metallaaromatic dibenzoylmethane coordination compounds and their aromatic analogs

The ability to predict the conductive behaviour of molecules, connected to macroscopic electrodes, represents a crucial prerequisite for the design of nanoscale electronic devices. In this work, we investigate whether the notion of a negative relation between conductance and aromaticity (the so-called NRCA rule) also pertains to quasi-aromatic and metallaaromatic chelates derived from dibenzoylmethane (DBM) and Lewis acids (LAs) that either do or do not contribute two extra d_π electrons to the central resonance-stabilised β-ketoenolate binding pocket. We therefore synthesised a family of methylthio-functionalised DBM coordination compounds and subjected them, along with their truly aromatic terphenyl and 4,6-diphenylpyrimidine congeners, to scanning tunneling microscope break-junction (STM-BJ) experiments on gold nanoelectrodes. All molecules share the common motif of three π-conjugated, six-membered, planar rings with a meta-configuration at the central ring. According to our results, their molecular conductances fall within a factor of ca. 9 in an ordering aromatic < metallaaromatic < quasi-aromatic. The experimental trends are rationalised by quantum transport calculations based on density functional theory (DFT).

Barrier Height Prediction by Machine Learning Correction of Semiempirical Calculations

Different machine learning (ML) models are proposed in the present work to predict density functional theory-quality barrier heights (BHs) from semiempirical quantum mechanical (SQM) calculations. The ML models include a multitask deep neural network, gradient-boosted trees by means of the XGBoost interface, and Gaussian process regression. The obtained mean absolute errors are similar to those of previous models considering the same number of data points. The ML corrections proposed in this paper could be useful for rapid screening of the large reaction networks that appear in combustion chemistry or in astrochemistry. Finally, our results show that 70% of the features with the highest impact on model output are bespoke predictors. This custom-made set of predictors could be employed by future Δ-ML models to improve the quantitative prediction of other reaction properties.

Quest for the Most Aromatic Pathway in Charged Expanded Porphyrins

Despite the central role of aromaticity in the chemistry of expanded porphyrins, the evaluation of aromaticity remains difficult for these extended macrocycles. The presence of multiple conjugation pathways and different planar and nonplanar π-conjugation topologies makes the quantification of global and local aromaticity even more challenging. In neutral expanded porphyrins, the predominance of the aromatic conjugation pathway passing through the imine-type nitrogens and circumventing the amino NH groups is established. However, for charged macrocycles, the question about the main conjugation circuit remains open. Accordingly, different conjugation pathways in a set of neutral, anionic, and cationic expanded porphyrins were investigated by means of several aromaticity indices rooted in the structural, magnetic, and electronic criteria. Overall, our results reveal the predominance of the conjugation pathway that passes through all nitrogen atoms to describe the aromaticity of deprotonated expanded porphyrins, while the outer pathway through the perimeter carbon atoms becomes the most aromatic in protonated macrocycles. In nonplanar and charged macrocycles, a discrepancy between electronic and magnetic descriptors is observed. Nevertheless, our work demonstrates AV_min remains the best tool to determine the main conjugation pathway of expanded porphyrins.

Acoplanarity, Aromaticity, Chirality, and Helical Twisting Power of Chlorin e6 13(N)-Methylamide-15,17-dimethyl Ester Complexes: Effect of a Metal

The experimental and theoretical study of the influence of metal complexing on geometry, aromaticity, chirality, and the ability to twist the nematic phase by complexes based on modified natural chlorin e6 was carried out. The geometry optimization of the chlorin e6 13(N)-methylamide-15,17-dimethyl ester (MADMECl) and its Zn, Cu, and Ni complexes by DFT (CAM-B3LYP/6–31 G(d,p) functional) method was performed. Based on these calculations, the acoplanarity degree of the macrocyclic ligand and the distortion energy of its dianion were estimated, which allowed the arrangement of the MADMECl complexes in the series Ni > Cu > Zn. Aromaticity was evaluated using the NICS criterion (nuclear independent chemical shift). An increase in the degree of aromaticity of the macrocycle upon complex formation was established. At the same time, the aromaticity of the inner conjugation contour corresponds to the same series as the acoplanarity, while the outer π-delocalization is characterized by the reverse sequence. An experimental evaluation of the electron circular dichroism of the Soret and the Q-bands, as well as the g-factor of dissymmetry, was carried out. The growth of these parameters with an increase in the degree of acoplanarity and aromaticity of the internal conjugation contour was determined. The induction of helical phases in mixtures of nematic liquid crystals (LCs) based on cyanobiphenyls and MADMECl macrocyclic metal complexes was studied by polarization microscopy, and the clearance temperatures and helix pitch of the mesophases were measured. A strong effect of the metal on the phase transition temperature and helical twisting power was established.

Divergent Synthesis of Alternant Bisanthenequinone and Nonalternant Heptalenodifluorenedione Ring Systems via a Concentration-Dependent Rearrangement

The synthesis of a diol containing a nonalternant aromatic core was investigated to access a nonalternant isomer of bisanthene with functional groups suitable for two-dimensional polymerization. An alternant diol and its nonalternant isomer were prepared in a short synthetic route from the same bifluorenylidene starting material. The bifluorenylidene reactant undergoes a Stone-Wales rearrangement in neat triflic acid, which unexpectedly provided both an alternant and nonalternant dione. The rearrangement was characterized by spectroscopy and single crystal X-ray diffraction of Grignard addition products of both isomers. The relative yield of the rearranged, alternant product increased along with the initial concentration of its polycyclic aromatic hydrocarbon (PAH) precursor, implicating a bimolecular rearrangement mechanism and enabling the divergent synthesis of both the nonalternant and alternant products. These findings offer convenient access to functional derivatives of two PAH classes of interest for their optoelectronic properties and serve as yet another warning about the importance of characterizing these materials with care, especially when insoluble products must be carried forward in a multistep synthetic route.

Atoms-In-Molecules' Faces of Chemical Hardness by Conceptual Density Functional Theory

The chemical hardness concept and its realization within the conceptual density functional theory is approached with innovative perspectives, such as the electronegativity and hardness equalization of atoms in molecules connected with the softness kernel, in order to examine the structure-reactivity equalization ansatz between the electronic sharing index and the charge transfer either in the additive or geometrical mean picture of bonding. On the other hand, the maximum hardness principle presents a relation with the chemical stability of the hardness concept. In light of the inverse relation between hardness and polarizability, the minimum polarizability principle has been proposed. Additionally, this review includes important applications of the chemical hardness concept to solid-state chemistry. The mentioned applications support the validity of the electronic structure principles regarding chemical hardness and polarizability in solid-state chemistry.

Experimental and Computational Studies on Bio-Inspired Flavylium Salts as Sensitizers for Dye-Sensitized Solar Cells

Six new bio-inspired flavylium salts were synthesized and investigated by a combined computational and experimental study for dye-sensitized solar cell applications. The compounds were characterized by FT-IR, UV-Vis, NMR spectroscopy, and LC-MS spectrometry techniques. The pH-dependent photochromic properties of the flavylium dyes were investigated through a UV-Vis spectroscopy study and revealed that they follow the same network of chemical reactions as anthocyanins upon pH changes. The structural and electronic properties of the dyes were investigated using density functional theory (DFT) and time-dependent density functional theory (TD-DFT). Geometry optimization calculation revealed that all dyes, regardless of the specie, flavylium cations or quinoidal bases, present a planar geometry. The photovoltaic performances of the dyes, in both flavylium and quinoidal base forms, were evaluated by the HOMO and LUMO energies and by calculating the light-harvesting efficiencies, the free energy change of electron injection, and the free energy change regeneration. The MO analysis showed that all dyes can inject electrons into the conduction band of the TiO₂ upon excitation and that the redox couple can regenerate the oxidized dyes. The results obtained for the free energy change of electron injection suggest that the quinoidal bases should inject electrons into the semiconductor more efficiently than the flavylium cations. The values for the free energy change regeneration showed that the redox electrolyte can easily regenerate all dyes. Dipole moment analysis was also performed. DSSCs based on the dyes, in both flavylium and quinoidal base forms, were assembled, and their photovoltaic performances were evaluated by measuring the open-circuit voltage, the short circuit current density, the fill factor, and the energy conversion efficiency. Results obtained by both experimental and computational studies showed that the overall performances of the DSSCs with the quinoidal forms were better than those obtained with the flavylium cations dyes.

Ambident reactivity of 5-aminopyrazoles towards donor-acceptor cyclopropanes

Lewis acid-catalysed reactions of donor-acceptor cyclopropanes with 1,3-disubstituted 5-aminopyrazoles were investigated. Under catalysis with gallium(III) chloride, products of the three-membered ring opening via a nucleophilic attack of the exocyclic amino group were obtained in a chemoselective manner. Oppositely, in the presence of scandium(III) triflate, products of either N-alkylation or C(4)-alkylation, or a mixture of both were formed. The products of the C(4) alkylation were transformed in one step into tetrahydropyrazolo[3,4-b]azepines that are attractive for medicinal chemistry and pharmacology.

Properties of the density functional response kernels and its implications on chemistry

An overview of mathematical properties of the non-local second order derivatives of the canonical, grand canonical, isomorphic, and grand isomorphic ensembles is given. The significance of their positive or negative semidefiniteness and the implications of these properties for atoms and molecules are discussed. Based on this property, many other interesting properties can be derived, such as the expansion in eigenfunctions, bounds on the diagonal and off-diagonal elements, and the eigenvalues of these kernels. We also prove Kato's theorem for the softness kernel and linear response and the dissociation limit of the linear responses as the sum of the linear responses of the individual fragments when dissociating a system into two non-interacting molecular fragments. Finally, strategies for the practical calculation of these kernels, their eigenfunctions, and their eigenvalues are discussed.

Application of DFT Calculations in Designing Polymer-Based Drug Delivery Systems: An Overview

Drug delivery systems transfer medications to target locations throughout the body. These systems are often made up of biodegradable and bioabsorbable polymers acting as delivery components. The introduction of density functional theory (DFT) has tremendously aided the application of computational material science in the design and development of drug delivery materials. The use of DFT and other computational approaches avoids time-consuming empirical processes. Therefore, this review explored how the DFT computation may be utilized to explain some of the features of polymer-based drug delivery systems. First, we went through the key aspects of DFT and provided some context. Then we looked at the essential characteristics of a polymer-based drug delivery system that DFT simulations could predict. We observed that the Gaussian software had been extensively employed by researchers, particularly with the B3LYP functional and 6-31G(d, p) basic sets for polymer-based drug delivery systems. However, to give researchers a choice of basis set for modelling complicated organic systems, such as polymer–drug complexes, we then offered possible resources and presented the future trend.

Synthesis, characterization and computational studies of 1,3-bis[(-furan-2-yl)methylene]urea and 1,3-bis[(-furan-2-yl)methylene]thiourea

Urea and thiourea derivatives: 1,3-bis[(E)-furan-2-yl)methylene]urea (BFMU) and 1,3-bis[(E)-furan-2-yl)methylene]thiourea (BFMT) were synthesized and characterized by spectrometry analyses (UV, IR, 1H NMR and 13C NMR). They were screened for antibacterial (Salmonella typhi, Staphylococcus aureus, Pseudomonas aeruginosa, Xanthomonas axonopodis and Streptococcus bovis) and antifungal (Fusarium oxysporum, Colletotrichum gloeosporioides and Cercospora zeae-maydis) activities. Quantum chemical calculations of frontier molecular orbital energies (EHOMO and ELUMO), and their associated global parameters were carried out by DFT levels of theory, with complete relaxation in the potential energy surface using 6-31G* basis set (DFT/B3LYP/6-31G*). Azomethine functional groups (C=N) appeared at δ 7.6 ppm and δ 7.0 ppm in the proton spectra, the peaks between δ 105 and δ 160 ppm of 13C spectra represent the methylene carbons (C=C). BFMU is a better inhibitor of P. aeruginosa and S. bovis, while BFMT is a better inhibitor of S. typhi, S. aureus and X. axonopodis and the fungi isolates (F. oxysporum, C. gloeosporioides and C. zeae-maydis) used. The global parameters agreed favorably with the experimental results, indicating the higher activity of BFMT.

Characterization of resonance structures in aromatic rings of benzene and its heavier-element analogues

We present the experimental visualization of the valence-electron-density distribution in benzene and its kinetically stabilized heavier-element analogues, i.e., 1,2-disilabenzene and 1,2-digermabenzene. The valence-electron-density-distribution (EDD) analysis on the 1,2-disila- and 1,2-digermabenzenes revealed that these contain incompletely delocalized π electrons on their cyclic conjugation systems, making them less aromatic compared to benzene. Based on the results of this EDD analysis in combination with anisotropy of the current-induced density (ACID) calculations, considerable contributions from the characteristic resonance structures of 1,2-disila- and 1,2-digermabenzenes with cleaved EE bonds can be expected.

Evolution of Triangular All-Metal Aromatic Complexes from Bonding Quandaries to Powerful Catalytic Platforms

This manuscript describes an overview on the literature detailing the observation of trinuclear complexes that present delocalized metal-metal bonds similar to those of regular aromatics, which are formed combining main group elements. A particular emphasis is given to the structural and electronic features of aromatic clusters that are sufficiently stable to allow their isolation. In parallel to the description of their key bonding properties, the work presents reported catalytic applications of these complexes, which already span from elaborated C-C-forming cascades to highly efficient cross-coupling methods. These examples present peculiar aspects of the unique reactivity exerted by all-metal aromatic complexes, which can often be superior to their established, popular mononuclear peers in terms of chemoselectivity and chemical robustness.

Ionic liquid treated leaves of Juglans regia as an adsorbent for the removal of methyl orange dye: experimental, computational, and statistical approach

The novel biosorbents prepared by surface modification from leaves of Juglans regia plant were exploited for removal of methyl orange dye from aqueous solution. The leaves in the form of dust and charcoal were separately impregnated with 1-butyl-3-methyl imidazolium bromide (I) to obtain adsorbents namely J. regia dust/charcoal impregnated with I (JRDI/JRCI) which were characterized using advanced analytical approaches. The impregnation of ionic liquid was confirmed by the appearance of new bands. Langmuir isotherm fitted well; the calculated adsorption capacity being 59.37 (JRDI) and 102.72 mg g^-1 (JRCI). The kinetic study revealed that sorption obeyed the pseudo-first order model; the experimental adsorption capacity being 53.53 (JRDI) and 86.82 mg g^-1 (JRCI) at selected conditions of pH 3, initial dye concentration 100 ppm, dosage of adsorbent 0.3 g and contact time 70 min. The mathematical models which predicted adsorption capacity as 51.5 (JRDI) and 82.1 mg g^-1 (JRCI) were found at par with experimental values. Fukui condensed functions revealed that adsorbents had electron deficient electrophilic reaction sites while dye had electron-rich nucleophilic reaction sites. The structural properties and good adsorption capability of adsorbents indicate that they could be used as potential, eco-friendly adsorbents for the treatment of negatively charged dye pollutants.

Doping P atom with a lone pair: an effective strategy to realize high HER catalytic activity and avoid deactivation under wide H* coverage on 2D silicene and germanene by increasing the structural rigidity

In view of the weak aromatic characteristic resulting from the weak π-bonding ability (different from the analogous graphene), employing two-dimensional (2D) silicene and germanene monolayers could be one of the most promising ways to realize a new type of highly efficient and nonprecious catalyst for the hydrogen evolution reaction (HER). However, the HER activity of pristine silicene and germanene has to be improved, although both of them can exhibit a good change trend. Particularly, the hydrogen phenomenon can occur under moderate or high H* coverage on 2D silicene and germanene. To overcome these bottlenecks, in this study we identify the most effective strategy through doping P with a lone pair to significantly improve the HER catalytic activity under a high H* coverage, by screening a series of IIIA (i.e., B, Al, Ga, In and Tl) and VA (i.e., N, P, As, Sb and Bi) heteroatoms with different electronegativity under detailed DFT calculations. It is revealed that the doped P atoms and almost all the Si/Ge atoms can uniformly serve as highly active sites. Especially, in view of the existence of the lone pair, doping P effectively prevents hydrogenation (even under full H* coverage) by increasing the structural rigidity. Moreover, the P-doping concentration also plays a crucial role in obtaining high HER activity. The relevant mechanisms have been analyzed in detail. Clearly, all these fascinating findings are beneficial for realizing new HER electrocatalysts based on the excellent silicene or germanene nanomaterials, and even other Si/Ge-related materials in the near future.

Conformational Properties and Putative Bioactive Targets for Novel Thiosemicarbazone Derivatives

The structure assignment and conformational analysis of the thiosemicarbazones, DKI21 and DKI24, were performed through homonuclear and heteronuclear 2D Nuclear Magnetic Resonance (NMR) spectroscopy (2D-COSY, 2D-NOESY, 2D-ROESY, 2D-HSQC, and 2D-HMBC) and quantum mechanics (QM) calculations, using Functional Density Theory (DFT). In addition, utilizing a combination of 2D-NOESY and 2D-ROESY spectra an exo structure was established for both of the analogs. This experimental results were confirmed by theoretical mechanistic studies, as the lowest minima conformations derived through DFT calculations were compatible with the spatial correlations observed in the 2D-NOESY and 2D-ROESY spectra. Finally, molecular binding experiments were performed to detect the potential targets for DKI21 and DKI24, derived from SwissAdme. In silico molecular binding experiments showed favorable binding energy values for the most of the enzymes studied. The ADMET calculations, using the preADMET and pKCSm software, showed that the two molecules appear as possible drug leads.

Structures, Energetics, and Spectra of (NH) and (OH) Tautomers of 2-(2-Hydroxyphenyl)-1-azaazulene: A Density Functional Theory/Time-Dependent Density Functional Theory Study

Tautomerization of 2-(2-hydroxyphenyl)-1-azaazulene (2OHPhAZ) in the gas phase and ethanol has been studied using B3LYP, M06-2X, and ωB97XD density functional theory (DFT) with different basis sets. For more accurate data, energies were refined at CCSD(T)/6-311++G(2d,2p) in the gas phase. Nuclear magnetic resonance (NMR), aromaticity, Fukui functions, acidity, and basicity were also calculated and compared with experimental data. Time-dependent density functional theory (TDDFT)-solvation model based on density (TDDFT-SMD) calculations in acetonitrile have been utilized for the simulation of UV-vis electronic spectra. In addition, electronic structures of the investigated system have been discussed. The results reveal that the enol form (2OHPhAZ) is thermodynamically and kinetically stable relative to the keto tautomer (2OPhAZ) and different rotamers (2OHPhAZ-R1:R3) in the gas phase and ethanol. A comparison with the experiment illustrates a good agreement and supports the computational findings.

Three types of noncovalent interactions studied between pyrazine and XF

Three types noncovalent interactions (type I, II and III) between pyrazine (C₄H₄N₂) and XF (X = F, Cl, Br, and I) have been discovered at the MP2/aug-cc-pVTZ level. TypeI is σ-hole interaction between the positive site on the halogen X of XF and the negative site on one of the pyrazine nitrogens. Type II is counterintuitive σ-hole interaction driven by polarization between the positive site on the halogen X of XF and a portion of the pyrazine ring. Type III is an interaction between the lateral regions of the halogen X of XF and the position of the pyrazine ring. Through comparing the calculated interaction energy, we can know that the type II and type III interactions are weaker than the corresponding type I interactions, and type III interactions are weaker than the corresponding type II interactions in C₄H₄N₂-XF complexes. SAPT analysis shows that the electrostatic energy are the major source of the attraction for the type I (σ-hole) interactions while the type III interactions are mainly dispersion energy. For the type II (counterintuitive σ-hole) interactions in C₄H₄N₂-XF (X = F and Cl) complexes, electrostatic energy are the major source of the attraction, while in C₄H₄N₂-XF (X = Br and I) complexes, the electrostatic term, induction and dispersion play equally important role in the total attractive interaction. NBO analysis, AIM theory, and conceptual DFT are also being utilized.

Novel Third-Order Nonlinear Optical Materials with Craig-Möbius Aromaticity

Electron delocalization in aromatic materials significantly impacts their third-order nonlinear optics (NLO). Despite organometallic complexes with Craig-Möbius aromaticity attracting great attention for their unusual physicochemical properties, their third-order NLO have been little studied to date. Herein, 12 Craig-Möbius aromatic organometallics with a stable structure similar to osmapentalyne, namely, carbolong complexes, are screened by DFT. They exhibit high third-order NLO responses because of the d and p electron delocalization in the organometallic ring. Furthermore, electron-hole distribution analyses draw a conclusion that extending the conjugated plane will increase the π-conjugation system to enhance the local excitation in the plane, and the introduction of typical aromatic ligands can result in the organometallic ring-to-ligand charge transfer (RLCT), which are effective methods to improve the third-order NLO response. This study opens a new window in the application of Craig-Möbius aromatic complexes and provides a new approach for third-order NLO materials design.

In silico studies on phytochemicals to combat the emerging COVID-19 infection

The current COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and its variants, remains a serious health hazard globally. The SARS-CoV-2 Mpro and spike proteins, as well as the human ACE2 receptor, have previously been reported as good targets for the development of new drug leads to combat COVID-19. Various ligands, including synthetic and plant-derived small molecules, can interact with the aforementioned proteins. In this study, we investigated the interaction of eight phytochemicals, from selected medicinal plants (Aegle marmelos, Azadirachta indica, and Ocimum sanctum) commonly used in Indian traditional medicine, with SARS-CoV-2 Mpro (PDBID: 6LU7), SARS-CoV-2S spike protein (PDB ID: 6M0J) and the human ACE2 receptor (PDB ID: 6M18). All compounds were subjected to density functional theory (DFT) and frontier molecular orbitals (FMO) analysis to determine their geometry, and key electronic and energetic properties. Upon examining the interactions of the phytochemicals with the human ACE2 receptor and the SARS-CoV-2 Mpro, spike protein targets, two compounds (C-5 and C-8) were identified as the best binding ligands. These were further examined in MD simulation studies to determine the stability of the ligand–protein interactions. QSAR, pharmacokinetic and drug-likeness properties studies revealed that C-5 may be the best candidate to serve as a template for the design and development of new drugs to combat COVID-19.

Is Aromaticity a Driving Force in Catalytic Cycles? A Case from the Cycloisomerization of Enynes Catalyzed by All-Metal Aromatic Pd3+ Clusters and Carboxylic Acids

The work details a mechanistic study based on density functional theory modeling on the cycloisomerization of polyunsaturated substrates catalyzed by all-metal aromatic tripalladium complexes and carboxylic acids. These clusters are an emerging class of catalysts for a variety of relevant transformations, including C-C forming processes that occur under mild conditions and display synthetic features complementary to those of established mononuclear complexes. This study is the first computational one devoted to the comprehension of the series of elementary steps involved in a synthetic transformation catalyzed by an all-metal aromatic complex. Present results confirm previous experimental hints on the striking mechanistic differences exerted by these clusters with respect to the usual cyclization pathways of related substrates. Moreover, the catalytic cycle involving present all-metal aromatic clusters closely parallels the mechanism of the aromatic substitution of regular arenes.

Spectral, NBO, NLO, NCI, aromaticity and charge transfer analyses of anthracene-9,10-dicarboxaldehyde by DFT

Anthracene-9,10-dicarboxaldehyde (ADCA) is a polynuclear aromatic compound that has a planar structure with double bonds which are in conjugation. The molecule is subjected to theoretical investigation with DFT/B3LYP/6-311++G(d,p) basis set to find out the electronic structural properties and stability. Theoretical and experimental vibrational analyses are carried out. NBO studies predict that the molecule has high stability. NCI interaction studies reveal that Van der Waals force and steric effect are seen in the molecule. A shaded surface map with a projection of LOL analysis pointed out that electron depletion area is seen in this molecule. The tunnelling current is more in the boundary rings than the central ring. It is docked with the protein 4COF and the binding energy is found to be -6.6 kcal/mol. Electrons excitation analysis is performed and found that local excitation takes place for the lowest five excitations. The aromaticity of the molecule is also thoroughly investigated.

Tuning the Computational Evaluation of Spectroscopic, ELF, LOL, NCI analysis and Molecular Docking of Novel Anti COVID-19 Molecule 4-Dimethylamino Pyridinium 3, 5-Dichlorosalicylate

In this work novel antiviral compound 4-(Dimethylamino) Pyridinium 3, 5-dichlorosalicylate was synthesized and characterized by UV-vis, FT-IR, FT-Raman, ¹H NMR and ¹³C NMR spectra. Quantum chemical computations were carried out by Density functional theory methods at B3LYP level. Electronic stability of the compound arising from hyper conjugative interactions and charge delocalization is investigated using natural bond orbital analysis. Assignments of vibrational spectra have been carried out with the aid of Normal coordinate analysis following the SQMFF methodology. TD-DFT approach was applied to assign the electronic transition observed in UV visible spectrum measured experimentally. Frontier molecular orbital energy gap affirms the bioactivity of the molecule and NCI analysis gives information about inter and intra non covalent interactions. ESP recognises the nucleophilic and electrophilic regions of molecule and the chemical implication of molecule was explained using ELF, LOL. The reactive sites of the compound were studied from the Fukui function calculations and chemical descriptors define the reactivity of the molecule. Molecular docking done with SARS and MERS proteins endorses the bioactivity of molecule and drug likeness factors were calculated to comprehend the biological assets of DADS.

From Macrocycles to Quantum Rings: Does Aromaticity Have a Size Limit?

ConspectusThe ring currents of aromatic and antiaromatic molecules are remarkable emergent phenomena. A ring current is a quantum-mechanical feature of the whole system, and its existence cannot be inferred from the properties of the individual components of the ring. Hückel's rule states that when an aromatic molecule with a circuit of [4n + 2] π electrons is placed in a magnetic field, the field induces a ring current that creates a magnetic field opposing the external field inside the ring. In contrast, antiaromatic rings with 4n π electrons exhibit ring currents in the opposite direction. This rule bears the name of Erich Hückel, and it grew from his molecular orbital theory, but modern formulations of Hückel's rule incorporate contributions from others, particularly William Doering and Ronald Breslow. It is often assumed that aromaticity is restricted to small molecular rings with up to about 22 π electrons. This Account outlines the discovery of global ring currents in large macrocycles with circuits of up to 162 π electrons. The largest aromatic rings yet investigated are cyclic porphyrin oligomers, which exhibit global ring currents after oxidation, reduction or optical excitation but not in the neutral ground state. The global aromaticity in these porphyrin nanorings leads to experimentally measurable aromatic stabilization energies in addition to magnetic effects that can be studied by NMR spectroscopy. Wheel-like templates can be bound inside these nanorings, providing excellent control over the molecular geometry and allowing the magnetic shielding to be probed inside the nanoring. The ring currents in these systems are well-reproduced by density functional theory (DFT), although the choice of DFT functional often turns out to be critical. Here we review recent contributions to this field and present a simple method for determining the ring current susceptibility (in nA/T) in any aromatic or antiaromatic ring from experimental NMR data by classical Biot-Savart calculations. We use this method to quantify the ring currents in a variety of aromatic rings. This survey confirms that Hückel's rule reliably predicts the direction of the ring current, and it reveals that the ring current susceptibility is surprisingly insensitive to the size of the ring. The investigation of aromaticity in even larger molecular rings is interesting because ring currents are also observed when mesoscopic metal rings are placed in a magnetic field at low temperatures. The striking similarity between the ring currents in molecules and mesoscopic metal rings arises because the effects have a common origin: a field-dependent phase shift in the electronic wave function. The main difference is that the magnetic flux through mesoscopic rings is much greater because of their larger areas, so their persistent currents are nonlinear and oscillatory with the applied field, whereas the flux through aromatic molecules is so small that their response is approximately linear in the applied field. We discuss how nonlinearity is expected to emerge in large molecular nanorings at high magnetic fields. The insights from this work are fundamentally important for understanding aromaticity and for bridging the gap between chemistry and mesoscopic physics, potentially leading to new functions in molecular electronics.

Electrostatic Potential Topology for Probing Molecular Structure, Bonding and Reactivity

Following the pioneering investigations of Bader on the topology of molecular electron density, the topology analysis of its sister field viz. molecular electrostatic potential (MESP) was taken up by the authors’ groups. Through these studies, MESP topology emerged as a powerful tool for exploring molecular bonding and reactivity patterns. The MESP topology features are mapped in terms of its critical points (CPs), such as bond critical points (BCPs), while the minima identify electron-rich locations, such as lone pairs and π-bonds. The gradient paths of MESP vividly bring out the atoms-in-molecule picture of neutral molecules and anions. The MESP-based characterization of a molecule in terms of electron-rich and -deficient regions provides a robust prediction about its interaction with other molecules. This leads to a clear picture of molecular aggregation, hydrogen bonding, lone pair–π interactions, π-conjugation, aromaticity and reaction mechanisms. This review summarizes the contributions of the authors’ groups over the last three decades and those of the other active groups towards understanding chemical bonding, molecular recognition, and reactivity through topology analysis of MESP.