Noncovalent Interactions in Palladium(II)-Catalyzed Meta-Selective C-H Functionalization: Mechanistic Insights and Origins of Regioselectivity.

The use of noncovalent interactions to control the regioselectivity of transition metal-catalyzed C-H functionalization of arenes has received significant attention in recent years. Herein, we present a mechanistic study based on Density Functional Theory (DFT) of palladium(II)-catalyzed meta-selective C-H olefination employing a noncovalent directing template. We analyze the key steps of the mechanism and discuss the origins of reaction selectivity. The role of the directing template was elucidated, demonstrating its essential function in lowering reaction barriers and controlling selectivity. Our results reveal a competition in activation between ortho- and meta-C-H bonds. Contrary to the previous proposal in the literature, hydrogen bonds between the N-H bonds of the urea moiety and the carbonyl oxygen of the substrate predominantly favor ortho-selectivity over meta-selectivity. DFT results, alongside Quantum Theory of Atoms in Molecules (QTAIM) and Non-Covalent Interaction Index analysis, suggest that secondary interactions between the R group linked to the urea moiety and the catalyst exert a more pronounced influence compared to the aforementioned hydrogen bonds, directing the selectivity towards the meta C-H bond.

Exploration of Free Energy Surface of the Au10 Nanocluster at Finite Temperature.

The first step in comprehending the properties of Au10 clusters is understanding the lowest energy structure at low and high temperatures. Functional materials operate at finite temperatures; however, energy computations employing density functional theory (DFT) methodology are typically carried out at zero temperature, leaving many properties unexplored. This study explored the potential and free energy surface of the neutral Au10 nanocluster at a finite temperature, employing a genetic algorithm coupled with DFT and nanothermodynamics. Furthermore, we computed the thermal population and infrared Boltzmann spectrum at a finite temperature and compared it with the validated experimental data. Moreover, we performed the chemical bonding analysis using the quantum theory of atoms in molecules (QTAIM) approach and the adaptive natural density partitioning method (AdNDP) to shed light on the bonding of Au atoms in the low-energy structures. In the calculations, we take into consideration the relativistic effects through the zero-order regular approximation (ZORA), the dispersion through Grimme's dispersion with Becke-Johnson damping (D3BJ), and we employed nanothermodynamics to consider temperature contributions. Small Au clusters prefer the planar shape, and the transition from 2D to 3D could take place at atomic clusters consisting of ten atoms, which could be affected by temperature, relativistic effects, and dispersion. We analyzed the energetic ordering of structures calculated using DFT with ZORA and single-point energy calculation employing the DLPNO-CCSD(T) methodology. Our findings indicate that the planar lowest energy structure computed with DFT is not the lowest energy structure computed at the DLPN0-CCSD(T) level of theory. The computed thermal population indicates that the 2D elongated hexagon configuration strongly dominates at a temperature range of 50-800 K. Based on the thermal population, at a temperature of 100 K, the computed IR Boltzmann spectrum agrees with the experimental IR spectrum. The chemical bonding analysis on the lowest energy structure indicates that the cluster bond is due only to the electrons of the 6 s orbital, and the Au d orbitals do not participate in the bonding of this system.

The Separability Problem in Molecular Quantum Systems: Information-Theoretic Framework for Atoms in Molecules.

Even though molecules are fundamentally quantum entities, the concept of a molecule retains certain classical attributes concerning its constituents. This includes the empirical separability of a molecule into its three-dimensional, rigid structure in Euclidean space, a framework often obtained through experimental methods like X-Ray crystallography. In this work, we delve into the mathematical implications of partitioning a molecule into its constituent parts using the widely recognized Atoms-In-Molecules (AIM) schemes, aiming to establish their validity within the framework of Information Theory concepts. We have uncovered information-theoretical justifications for employing some of the most prevalent AIM schemes in the field of Chemistry, including Hirshfeld (stockholder partitioning), Bader's (topological dissection), and the quantum approach (Hilbert's space definition). In the first approach we have applied the generalized principle of minimum relative entropy derived from the Sharma-Mittal two-parameter functional, avoiding the need for an arbitrary selection of reference promolecular atoms. Within the ambit of topological-information partitioning, we have demonstrated that the Fisher information of Bader's atoms conform to a comprehensive theory based on the Principle of Extreme Physical Information avoiding the need of employing the Schwinger's principle, which has been proven to be problematic. For the quantum approach we have presented information-theoretic justifications for conducting Löwdin symmetric transformations on the density matrix to form atomic Hilbert spaces generating orthonormal atomic orbitals with maximum occupancy for a given wavefunction.

Bond analysis in meta- and para-substituted thiophenols: overlap descriptors, local mode analysis, and QTAIM.

CONTEXT: This study delves into the chemical nuances of thiophenols and their derivatives through a comprehensive computational analysis, moving beyond traditional energetic perspectives such as bond dissociation enthalpy and S-H dissociation dynamics. By employing the overlap model along with its topological descriptors (OP/TOP), quantum theory of atoms in molecules (QTAIM), and local vibrational mode (LVM) theories, the research provides a deeper understanding of the S-H and C-S bonding scenarios in substituted thiophenols. The investigation follows the electron-donating capacity of S-H substituent variation with the nature and positioning of other ring substituents. Energy profile analyses indicate distinct stability differences in the cis and trans conformations of meta- and para-PhSH systems, influenced by the electron-donating strength of these substituents. The study also uncovers significant variations in S-H bond distances and descriptor values, particularly in para-substituted PhSH, reflecting the influence of electron-donating or withdrawing substituents. In contrast, alterations at the meta-position show minimal effects on C-S bond descriptors, while para-substitutions markedly influence C-S bond characteristics, demonstrating a clear correlation with the electron-donating or withdrawing capabilities of the substituents. This research sheds light on the intricate bond dynamics in aromatic systems with diverse substituents, highlighting the complex interaction between electronic effects and molecular conformation. METHODS: The study employs the ω B97X-D/Def2TZVP level of theory for molecular geometries, ensuring accurate characterization of structures as true minima via analytical harmonic frequency determination. The electronic properties of S-H and C-S bonds in variously substituted thiophenols were analyzed using OP/TOP, QTAIM, and LVM methodologies. Computational processes, including conformational scans, geometry optimizations, and vibrational frequency calculations, were conducted using Gaussian 09, with ultra-fine integration grids and tight convergence criteria for the SCF procedure. Bond descriptors were computed utilizing ChemBOS, Multiwfn, and LModeA software, providing a robust and detailed examination of bond properties.

Second generation of pyrimidin-quinolone hybrids obtained from virtual screening acting as sphingosine kinase 1 inhibitors and potential anticancer agents.

We report here the virtual screening design, synthesis and activity of eight new inhibitors of SphK1. For this study we used a pre-trained Graph Convolutional Network (GCN) combined with docking calculations. This exploratory analysis proposed nine compounds from which eight displayed significant inhibitory effect against sphingosine kinase 1 (SphK1) demonstrating a high level of efficacy for this approach. Four of these compounds also displayed anticancer activity against different tumor cell lines, and three of them (5), (6) and (7) have shown a wide inhibitory action against many of the cancer cell line tested, with GI50 below 5 µM, being (5) the most promising with TGI below 10 µM for the half of cell lines. Our results suggest that the three most promising compounds reported here are the pyrimidine-quinolone hybrids (1) and (6) linked by p-aminophenylsulfanyl and o-aminophenol fragments respectively, and (8) without such aryl linker. We also performed an exhaustive study about the molecular interactions that stabilize the different ligands at the binding site of SphK1. This molecular modeling analysis was carried out by using combined techniques: docking calculations, MD simulations and QTAIM analysis. In this study we also included PF543, as reference compound, in order to better understand the molecular behavior of these ligands at the binding site of SphK1.These results provide useful information for the design of new inhibitors of SphK1 possessing these structural scaffolds.

A DFT study of ethanol interaction with the bimetallic clusters of PtSn and its implications on reactivity.

The comprehension of the factors affecting the adsorption of ethanol over metals and metal alloys is a crucial step for the rational development of new catalysts for hydrogen production through ethanol reforming. In this work, we analyze the effect of combining Pt and Sn on a metal cluster on the complexation energy and reactivity for OH dehydrogenation of ethanol. Metal clusters of Pt10, Sn10 and Pt5Sn5 had their putative minimum located with the help of the artificial bee colony algorithm. Whereas the isolated Pt cluster shows a high degree of polarization (ESP surface), the Sn cluster shows a quite uniform electron density surface. The PtSn cluster is strongly polarized, with Pt atoms withdrawing electron density of Sn atoms. Complexation occurs with the oxygen atom of ethanol directed towards the point of highest electron potential in the ESP surface. Pt presents the highest complexation energy, -20.90 kcal/mol, against only -7.83 kcal/mol (at the B97-3c level). For the PtSn cluster, the value is intermediate, namely -12.39 kcal/mol. The more malleable electron density of Pt and its electron affinity are responsible for its highest complexation energy. These characteristics are partially transferred to the PtSn cluster. QTAIM results show that, for the PtSn cluster, the O-H bond in ethanol is somewhat weaker than for pure Pt and Sn. As a consequence, the energy barrier for the O-H dehydrogenation has its lowest value for the PtSn cluster, which shows that the alloying of two metals can lead to quite quite unexpected results opening the perspective for a more rational fine tuning of catalysts properties.

Theoretical investigation of anion perfluorocubane.

CONTEXT: In this work, we did a theoretical exploration of C8F8 (Ib) and its anion radical analogue (IIb) in this work. By investigating the thermochemistry of electron capture, we find that the free energy associated with the conversion of C8H8 (Ia) into its anion radical analogue IIa is of the order of + 92.83 kcal.mol-1, while the conversion of Ib into IIb is - 6.42 kcal.mol-1. Therefore, species IIb is thermodynamically more stable than its neutral analogue. Natural bond orbitals (NBO) analyses revealed that compound Ib exhibits a relative electronic stability as a function of intramolecular delocalisations of the type [Formula: see text] of the order of 2.70 kcal.mol-1. Similar delocalizations for Ia are energetically lower (1.45 kcal.mol-1). Topological analyses of compounds Ib and IIb indicate that the addition of an electron to Ib enhances the covalency of the C-C bond, as can be seen by the reduction in the ellipticity of the C-C bond. The opposite is observed for Ia, whose addition of the electron (leading to IIa) reduces the covalency of the C-C bond. By comparing the free and packaged forms of the species, it is found that, in the crystalline form, the system will present greater relative stability due to the dispersive interactions involved, as evidenced by non-covalent interactions (NCI) analysis. Finally, it was possible to verify that the manifestation of the current density with a lower paratropic and less antiaromatic character in Ib and IIb point to C8F8 as a strong candidate for electron capture. METHODS: Geometry optimization calculations were carried out, for all monomer structures using the hybrid functional B3LYP-D3 and the 6-31+G(d,p) basis set. To determine the formation thermochemistry of the ions, electronic energy corrections was performed using the DLPNO-CCSD(T)/aug-cc-pVTZ/C method. Starting from the optimised forms, shielding, nuclear magnetic resonance (NMR) spectra employing gauge-independent atomic orbital (GIAO), and NBO calculations were performed for these monomers, using the PBE0 functional and the pCSseg-2 atomic basis set. The magnetochemical analysis of ring currents was performed using the GIMIC formalism. For the topological analysis, it was applied the combination DLPNO-CCSD(T)/aug-cc-pVTZ/C, previously used for correcting the electronic energy.

Weak noncovalent interactions in two positional isomers of acrylonitrile derivatives: inputs from PIXEL energy, Hirshfeld surface and QTAIM analyses.

A single crystal X-ray diffraction analysis was performed on two positional isomers (m-tolyl and p-tolyl) of acrylonitrile derivatives, namely, (Z)-3-(4-(pyridin-2-yl) phenyl)-2-(m-tolyl) acrylonitrile (1) and (Z)-3-(4-(pyridin-2-yl)phenyl)-2-(p-tolyl) acrylonitrile (2). Compound 1 crystallized in the monoclinic P21/n space group with two crystallographically independent molecules. Compound 2 also possesses two crystallographically independent molecules and crystallized in the triclinic P-1 space group. The Hirshfeld surface analysis revealed that, in both isomers, intermolecular H⋅⋅⋅H/C/N contacts contribute significantly to the crystal packing. More than 40% of the contribution arises from intermolecular C-H⋅⋅⋅C(π) contacts. In both compounds, the relative contribution of these contacts is comparable, indicating that the positional isomeric effects are marginal. The structures in which these isomers are arranged in the solid state are very similar, and the lattice energies are also comparable between the isomers. The Coulomb-London-Pauli-PIXEL (CLP-PIXEL) energy analysis identified the energetically significant dimers. The strength of the intra- and intermolecular interactions was evaluated using the quantum theory of atoms in molecules approach. The UV-Vis absorbance in three different solvents (chloroform, ethanol, and ethyl acetate) for isomers 1 and 2 are very similar. This result is in good agreement with the time-dependent density-functional theory (TD-DFT) calculations.

(Ro)vibrational Spectroscopic Constants, Lifetime and QTAIM Evaluation of Fullerene Dimers Stability.

The iconic caged shape of fullerenes gives rise to a series of unique chemical and physical properties; hence a deeper understanding of the attractive and repulsive forces between two buckyballs can bring detrimental information about the structural stability of such complexes, providing significant data applicable for several studies. The potential energy curves for the interaction of multiple van der Waals buckyball complexes with increasing mass were theoretically obtained within the DFT framework at ωB97xD/6-31G(d) compound model. These potential energy curves were employed to estimate the spectroscopic constants and the lifetime of the fullerene complexes with the Discrete Variable Representation and with the Dunham approaches. It was revealed that both methods are compatible in determining the rovibrational structure of the dimers and that they are genuinely stable, i.e., long-lived complexes. To further inquire into the nature of such interaction, Bader's QTAIM approach was applied. QTAIM descriptors indicate that the interactions of these closed-shell systems are dominated by weak van der Waals forces. This non-covalent interaction character was confirmed by the RDG analysis scheme. Indirectly, QTAIM also allowed us to confirm the stability of the non-covalent bonded fullerene dimers. Our lifetime calculations have shown that the studied dimers are stable for more than 1 ps, which increases accordingly with the number of carbon atoms.

Exploring the Non-Covalent Bonding in Water Clusters.

QTAIM and source function analysis were used to explore the non-covalent bonding in twelve different water clusters (H2O)n obtained by considering n = 2-7 and various geometrical arrangements. A total of seventy-seven O-H⋯O hydrogen bonds (HBs) were identified in the systems under consideration, and the examination of the electron density at the bond critical point (BCP) of these HBs revealed the existence of a great diversity of O-H⋯O interactions. Furthermore, the analysis of quantities, such as |V(r)|/G(r) and H(r), allowed a further description of the nature of analogous O-H⋯O interactions within each cluster. In the case of 2-D cyclic clusters, the HBs are nearly equivalent between them. However, significant differences among the O-H⋯O interactions were observed in 3-D clusters. The assessment of the source function (SF) confirmed these findings. Finally, the ability of SF to decompose the electron density (ρ) into atomic contributions allowed the evaluation of the localized or delocalized character of these contributions to ρ at the BCP associated to the different HBs, revealing that weak O-H⋯O interactions have a significant spread of the atomic contributions, whereas strong interactions have more localized atomic contributions. These observations suggest that the nature of the O-H⋯O hydrogen bond in water clusters is determined by the inductive effects originated by the different spatial arrangements of the water molecules in the studied clusters.

Fragmentation route of doubly ionized benzene, aniline, and nitroanilines monomers using a novel protocol from density functional theory and QTAIM.

The possibility of finding the fragmentation routes by theoretical methods led us to compare the molecular ions between neutral molecules of benzene, aniline, and o-, m-, and p-nitroaniline, using the density functional theory (DFT), under an aug-cc-pVDZ base set and a B3LYP exchange-correlation functional. After determining the structure and electronic energy of neutral and doubly ionized species, we used a new protocol based on analyzing Wiberg's binding indexes and the quantum theory of atoms in Bader molecules (QTAIM). The charge transfer and electronic distribution in aromatic monomers indicate the possibility of fragment formation in at least two pairs of carbon-carbon (CC) atoms. They show the possible loss of the -CNH2 and -NO2 groups in the aniline and nitroaniline molecules doubly ionized.

Nature and Strength of Weak O⋠⋠⋠O Interactions in Nitryl Halide Dimers.

The use of real space functions and molecular graphs has pushed some chemists to wonder: Are interactions between negatively charged oxygen atoms possible? In this contribution we analyze whether there is a real interaction between oxygen atoms in nitryl halide dimers (XNO2 )2 (X=F, Cl, Br and I) and in tetranitromethane and derivatives. Based on ab-initio and density functional theories (DFT) methods, we show these complexes are weakly stabilized. Energy decomposition analyses based on local molecular orbitals (LMOEDA) and interacting quantum atoms (IQA) reveal both dispersion and exchange play a crucial role in the stabilization of these complexes. Electron charge density and IQA analyses indicate that the oxygen atoms are connected by privileged exchange channels. In addition, electrostatic interactions between O and N atoms are also vital for the stabilization of the complexes. Finally, a reasonable explanation is given for the dynamic behavior of nitryl groups in tetranitromethane and derivatives.

Energetics and electronics of polar Diels-Alder reactions at the atomic level: QTAIM and IQA analyses of complete IRC paths.

The mechanism of Diels-Alder reactions between cyclopentadiene and several cyanoethylenes was studied by means of Density Functional Theory calculations using QTAIM and IQA (Interacting Quantum Atoms) analyses along complete IRC paths. Each geometry from the IRC had its wavefunction computed and the topology of the electronic density for it was then evaluated. By means of IQA, the global energetic profile was partitioned among the various atoms in the molecule, providing insight into what atoms are the main ones responsible for the magnitude of the energy barriers. The (a)synchronicity of the reaction mechanisms featuring non-symmetrically substituted dienophiles was characterized, from QTAIM, by the electron densities and Laplacians over the LCP's as well as by the different atomic energy barriers obtained from IQA. The magnitude of the atomic barrier nicely explains the (a)synchronicity of the reaction mechanisms, and the degree of (a)synchronicity is nicely revealed by the difference between the earlier and later bond breaking and bond formations. The main conclusion is that important energetic and electronic changes are occurring before and after the position of the transition state structure, mainly for those asynchronous mechanisms, and although these provide essential insight into the reaction mechanism, most studies cannot assess this kind of information because they are focusing solely on reactants, transition states, and products. We advocate that the additional computational effort required for such analyses is more than compensated by the great amount of useful information it provides.

Dissecting Bonding Interactions in Cysteine Dimers.

Neutral (n) and zwitterionic (z) forms of cysteine monomers are combined in this work to extensively explore the potential energy surfaces for the formation of cysteine dimers in aqueous environments represented by a continuum. A simulated annealing search followed by optimization and characterization of the candidate structures afforded a total of 746 structurally different dimers held together via 80 different types of intermolecular contacts in 2894 individual non-covalent interactions as concluded from Natural Bond Orbitals (NBO), Quantum Theory of Atoms in Molecules (QTAIM) and Non-Covalent Interactions (NCI) analyses. This large pool of interaction possibilities includes the traditional primary hydrogen bonds and salt bridges which actually dictate the structures of the dimers, as well as the less common secondary hydrogen bonds, exotic X⋯Y (X = C, N, O, S) contacts, and H⋯H dihydrogen bonds. These interactions are not homogeneous but have rather complex distributions of strengths, interfragment distances and overall stabilities. Judging by their Gibbs bonding energies, most of the structures located here are suitable for experimental detection at room conditions.

Conceptual DFT, QTAIM, and Molecular Docking Approaches to Characterize the T-Type Calcium Channel Blocker Anandamide.

The anandamide is a relevant ligand due to its capacity of interacting with several proteins, including the T-type calcium channels, which play an important role in neuropathic pain and depression disorders. Hence, a detailed characterization of the chemical properties and conformational stability of anandamide may provide valuable information to understand its behavior in a biological context. Herein, conceptual DFT and QTAIM analyses were performed to theoretically characterize the chemical reactivity properties and the structural stability of conformations of anandamide, using the BP86/cc-pVTZ level of theory. Global reactivity description, based on conceptual DFT, indicates that the hardness increases and the electrophilicity index decreases for both, the hairpin and U-shape conformers relative to the extended conformers. Also, an increase in the chemical potential value and a decrease in the electronegativity and the electrophilicity index is observed in the ethanolamide open ring conformers in comparison with the corresponding closed ring structures. In addition, regarding the characterization of local reactivity descriptors, the maximum values of the Fukui and Parr functions indicate that the most probable location for a nucleophilic attack is either the hydroxyl oxygen located in the ethanolamide closed ring conformers or the carbonyl oxygen present in the open ring conformers. The most probable location for an electrophilic attack is in the alkyl double bond region in all anandamide conformers. According to the QTAIM results, the intramolecular hydrogen bond formation stabilizing the structure of anandamide has interaction energy values for the closed ring conformations of 12.33-12.46 kcal mol-1, indicating a strong interaction. Lastly, molecular docking calculations determined that a region in the pore, denominate as pore-blocking, is a probable site for the interaction of anandamide with the human Cav3.2 isoform of the T-type calcium channel family. The pore-blocking site contains hydrophobic residues where the non-polar part in the final alkyl region of anandamide established mainly alkyl-alkyl interactions, while the polar part (the ethanolamide group) interacts with the polar residue S900. The information based on conceptual DFT presented may aid in the design of drugs with similar chemical characteristics as those identified in anandamide so as to bind anandamide-interacting proteins, including the T-type calcium channels.

Microsolvation versus Encapsulation in Mono, Di, and Trivalent Cations.

The effects of the formal charge in the stability and bonding of water cavities when solvating a cation are studied here using [X(H2 O)20 ]q+ clusters starting with the well known 512 isomer of (water)20 , placing a single mono, di, or trivalent Xq+ cation at the interior, and then optimizing and characterizing the resulting clusters. Highly correlated interaction and deformation energies are calculated using the CCSD(T)-DLPNO formalism. Bonding interactions are characterized using the tools provided by the quantum theory of atoms in molecules, natural bond orbitals, and non-covalent surfaces. Our results indicate that water to water hydrogen bonds are sensibly strengthened resulting in strong cooperative effects, which amount to ≈ 2 ${ \approx 2}$  kcal/mol per hydrogen bond in the bare cavity and to larger values for the systems including the cations. Approximate encapsulation, that is, surrounding the cation by a network of hydrogen bonds akin to the well known methane clathrate seems to be preferred by cations with smaller charge densities while microsolvation, that is, cluster structures having explicit X⋯O contacts seem to be preferred by cations with larger charge densities which severely deform the cavity.

A new quinolinone-chalcone hybrid with potential antibacterial and herbicidal properties using in silico approaches.

Quinolinone-chalcones are hybrid compounds consisting of chalcone and quinolone moieties with biological activity related to their hybrid structure. This work seeks to describe the structural and theoretical parameters related to the physicochemical properties and biological activity of a new quinolinone-chalcone. The synthesis, structural characterization by X-ray diffraction, molecular topology by Hirshfeld surfaces and QTAIM, molecular electronic calculations, and pharmacophore analysis were described. The weak interactions C-H…O, C-H…π, and C-H…Br were responsible for crystal growth and stabilized the crystalline state. The DFT analysis shows that the sulfonamide group region is susceptible to observed interactions, and the frontier molecular orbitals indicate high kinetic stability. Also, pharmacophore analysis revealed potential antibacterial and herbicidal activity; by docking within the active site of TtgR, a transcription regulator for the efflux pump TtgABC from the highly resistant Pseudomonas putida (P. putida) strain DOT-TIE, we showed that the activation of TtgR relies upon the binding of aromatic-harboring compounds, which plays a crucial role in bacterial evasion. In this context, a new quinolinone-chalcone has a higher binding affinity than tetracycline, which suggests it might be a better effector for TtgR.

Evaluating the conformational space of the active site of D2 dopamine receptor. Scope and limitations of the standard docking methods.

We report here for the first time the potential energy surfaces (PES) of phenyletilamine (PEA) and meta-tyramine (m-OH-PEA) at the D2 dopamine receptor (D2DR) binding site. PESs not only allow us to observe all the critical points of the surface (minimums, maximums, and transition states), but also to note the ease or difficulty that each local minima have for their conformational inter-conversions and therefore know the conformational flexibility that these ligands have in their active sites. Taking advantage of possessing this valuable information, we analyze how accurate a standard docking study is in these cases. Our results indicate that although we have to be careful in how to carry out this type of study and to consider performing some extra-simulations, docking calculations can be satisfactory. In order to analyze in detail the different molecular interactions that are stabilizing the different ligand-receptor (L-R) complexes, we carried out quantum theory of atoms in molecules (QTAIM) computations and NMR shielding calculations. Although some of these techniques are a bit tedious and require more computational time, our results demonstrate the importance of performing computational simulations using different types of combined techniques (docking/MD/hybrid QM-MM/QTAIM and NMR shielding calculations) in order to obtain more accurate results. Our results allow us to understand in details the molecular interactions stabilizing and destabilizing the different L-R complexes reported here. Thus, the different activities observed for dopamine (DA), m-OH-PEA, and PEA can be clearly explained at molecular level.

A fast approximate extension of the interacting quantum atoms energy decomposition to excited states.

An approach is developed for the fast calculation of the interacting quantum atoms energy decomposition (IQA) from the information contained in the first order reduced density matrix only. The proposed methodology utilizes an approximate exchange-correlation density from Density Matrix Functional Theory without the need to evaluate the correlation-exchange contribution directly. Instead, weight factors are estimated to decompose the exact Vxc into atomic and pairwise contributions. In this way, the sum of the IQA contributions recovers the energy obtained from the electronic structure calculation. This method can, hence, be applied to obtain atomic contributions in excited states on the same footing as in their ground states using any method that delivers the reduced first-order density matrix. In this way, one can locate chromophores from first principles quantum chemical calculations. Test calculations on the ground and excited states of a set of small molecules indicate that the scaled atomic contributions reproduce vertical electronic transition energies calculated exactly. This approach may be useful to extend the applicability of the IQA approach in the study of large photochemical systems especially when the calculations of the second order reduced density matrices is prohibitive or not possible.

Synthesis, Structural Characterization, and In Vitro and In Silico Antifungal Evaluation of Azo-Azomethine Pyrazoles (PhN2(PhOH)CHN(C3N2(CH3)3)PhR, R = H or NO2).

The azo-azomethine imines, R1-N=N-R2-CH=N-R3, are a class of active pharmacological ligands that have been prominent antifungal, antibacterial, and antitumor agents. In this study, four new azo-azomethines, R1 = Ph, R2 = phenol, and R3 = pyrazol-Ph-R' (R = H or NO2), have been synthesized, structurally characterized using X-ray, IR, NMR and UV-Vis techniques, and their antifungal activity evaluated against certified strains of Candida albicans and Cryptococcus neoformans. The antifungal tests revealed a high to moderate inhibitory activity towards both strains, which is regulated as a function of both the presence and the location of the nitro group in the aromatic ring of the series. These biological assays were further complemented with molecular docking studies against three different molecular targets from each fungus strain. Molecular dynamics simulations and binding free energy calculations were performed on the two best molecular docking results for each fungus strain. Better affinity for active sites for nitro compounds at the "meta" and "para" positions was found, making them promising building blocks for the development of new Schiff bases with high antifungal activity.