Dynamic behavior of chemical reactivity indices in density functional theory: A Bohn-Oppenheimer quantum molecular dynamics study

Design of New Schiff-Base Copper(II) Complexes: Synthesis, Crystal Structures, DFT Study, and Binding Potency toward Cytochrome P450 3A4

We report the synthesis and crystal structures of three new copper(II) Schiff-base complexes. The complexes have been characterized by elemental analysis and Fourier transform infrared (FT-IR) and UV-visible spectroscopies. The X-ray diffraction (XRD) analysis reveals that complexes 1 and 3 crystallize in a monoclinic space group C2/c and 2 in a triclinic space group P1̅, each adopting a square planar geometry around the metal center. We use a density functional theory method to explore the quantum chemical properties of these complexes. The calculation proceeds with the three-dimensional (3D) crystal structure characterization of the complexes in which the calculated IR and UV-vis values are comparable to the experimental results. Charge distribution and molecular orbital analyses enabled quantum chemical property prediction of these complexes. We study the drug-likeness properties and binding potentials of the synthesized complexes. The in silico outcome showed that they could serve as permeability-glycoprotein (P-gp) and different cytochrome P450 substrates. Our calculations showed that the complexes significantly bind to cytochrome P450 3A4.

Reactivity Dynamics

The chemical reactivity of a molecule as a whole or of an atom in a molecule varies during a chemical reaction. A variation of global and local reactivity descriptors in the course of a physicochemical process was studied within a quantum fluid density functional theory framework. Effects of a physical confinement and the electronic excitation therein were studied. In this Perspective, we also highlight the direction of a spontaneous chemical reaction in the light of the dynamical variants of the conceptual density functional theory-based electronic structure principles. An exhaustive state-of-the-art dynamical study is warranted in order to understand a chemical reaction from a reactivity perspective augmenting the associated molecular reaction dynamics analysis.

Multi-target approach for Alzheimer's disease treatment: computational biomolecular modeling of cholinesterase enzymes with a novel 4-N-phenylaminoquinoline derivative reveal promising potentials

The identification of dual inhibitors targeting the active sites of the cholinesterase enzymes, acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE), have lately surfaced as a multi-approach towards Alzheimer treatment. More recently, a novel series of 4-N-phenylaminoquinolines was synthesized and evaluated against AChE and BuChE in which one of the compounds displayed appreciable inhibition compared to the standard compound, galantamine. To provide a clearer picture of the inhibition mechanism of this potent compound at the molecular level, computational biomolecular modeling was carried out. The investigation was initiated with the exploration of the chemical properties of the identified compound 11 b and reference drug, galantamine. Density functional theory (DFT) calculations reveal some conceptual parameters that provide information on the stability and reactivity of the compounds as potential inhibitors. To unveil the binding mechanism, energetics and enzyme-ligand interactions, molecular dynamics (MD) simulations of six different systems were executed over a period. Calculated binding free energy values are in the same order with experimental IC₅₀ data. Identification of the main residues driving optimum binding of the active compound 11 b to the binding region of both AChE and BuChE showed Trp81 and Trp110 as the most important, respectively. It was proposed that the studied compound could serve as a dual inhibitor for AChE and BuChE, therefore, would potentially be a promising moiety in a multi-target approach for the treatment of Alzheimer's disorder.Communicated by Ramaswamy H. Sarma.

Broad-band luminescence involving fluconazole antifungal drug in a lead-free bismuth iodide perovskite: Combined experimental and computational insights

The synthesis and characterization of a lead-free perovskite-type material, (C₁₃H₁₄N₆F₂O)₂ Bi₂I₁₀ is reported. It exhibits a zero-dimensional (0D) Bi₂I₁₀^4- octahedral unit, surrounded by a flexible tripodal antifungal ligand (H₂Fluconazole)²⁺. The several intermolecular interactions of the independent cation and the bismuth iodide octahedra were tested via the Hirshfeld surface analysis. The detailed interpretation of the vibrational modes was carried out. The band gap (Eg) of 2.10 eV agrees with the theoretical values. Upon photoexcitation, the crystals exhibit a broadband green emission peaked at 534 nm, which originates from electronic transitions within the inorganic cluster [Bi₂I₁₀]^4-. The theoretical calculations were carried out using DFT and TD-DFT methods to appraise the molecular geometry, vibrational spectra, electronic absorption spectra, frontier molecular orbitals (FOMs) and global reactivity descriptors. Calculations reveal that the energy gap (Eg) and other chemical reactivity descriptors are primarily linked to the inorganic anion and the triazolium rings (A and B) of the organic cation reflecting their importance in the activity and the antioxidant ability of the molecule.

Atom-Condensed Fukui Function in Condensed Phases and Biological Systems and Its Application to Enzymatic Fixation of Carbon Dioxide

Local reactivity descriptors such as atom-condensed Fukui functions are promising computational tools to study chemical reactivity at specific sites within a molecule. Their applications have been mainly focused on isolated molecules in their most stable conformation without considering the effects of the surroundings. Here we propose to combine quantum mechanics/molecular mechanics Born-Oppenheimer molecular dynamics simulations to obtain the microstates (configurations) of a molecular system using different representations of the molecular environment and calculate Boltzmann-weighted atom-condensed local reactivity descriptors based on conceptual density functional theory. Our approach takes the conformational fluctuations of the molecular system and the polarization of its electron density by the environment into account, allowing us to analyze the effect of the molecular environment on reactivity. In this contribution, we apply the method mentioned above to the catalytic fixation of carbon dioxide by crotonyl-CoA carboxylase/reductase and study if the enzyme alters the reactivity of its substrate compared with an aqueous solution. Our main result is that the protein environment activates the substrate by the elimination of solute-solvent hydrogen bonds from aqueous solution in the two elementary steps of the reaction mechanism: the nucleophilic attack of a hydride anion from NADPH on the α,β-unsaturated thioester and the electrophilic attack of carbon dioxide on the formed enolate species.

Quantum chemical study of a derivative of 3-substituted dithiocarbamic flavanone

The aim of this work is to characterize a quite novel 3-dithiocarbamic flavonoid by vibrational spectroscopy in conjunction with Density Functional Theory (DFT) calculations. Quantum mechanics calculations of energies, geometries and vibrational wavenumbers in the ground state were carried out by using hybrid functional B3LYP with 6-311G(d,p) as basis set. The results indicate a remarkable agreement between the calculated molecular geometries, as well as vibrational frequencies, and the corresponding experimental data. In addition, a complete assignment of all the absorption bands present in the vibrational spectrum has been performed. In order to assess its chemical potential, quantum molecular descriptors characterizing the interactions between the 3-dithiocarbamic flavonoid and its biological receptors have been computed. The frontier molecular orbitals and the HOMO-LUMO energy gap have been used in order to explain the way in which the new molecule can interact with other species and to characterize its molecular chemical stability/reactivity. The molecular electrostatic potential (MEP) map, computed in order to identify the sites of the studied flavonoid that are most likely to interact with electrophilic and nucleophilic species, is discussed.

Quantum chemical aspects of solvent effects on the Diels–Alder reaction of 2,3-dimethyl-1,3-butadiene and diethyl azodicarboxylate

A comprehensive theoretical study was performed regarding the solvent effect on the Diels–Alder reaction of 2,3-dimethyl-1,3-butadiene and the hetero-dienophile, diethyl azodicarboxylate. Reaction rates in the various solutions were studied using SCRF–DFT theory and good agreement with the experimental results was obtained. According to natural bond orbital analysis, it was confirmed that the stabilisation energy of the C–N bonds contributes to the strong interaction between the solvents and transition state (TS) structures. Topological analyses showed a linear correlation between the interaction energy and the sum of electron density at the bond critical point of C–N at the TS. Moreover, with an increase in the overall electron density of the C–N bond critical point, the vibrational frequencies were decreased. Finally, some correlations between the quantum reactivity indices and the rate constants in the presence of 23 solvents were investigated.

Directed electrostatic activation in enantioselective organocatalytic cyclopropanation reactions: a computational study

Cyclopropane rings are versatile building blocks in organic chemistry. Their synthesis, by the reaction of sulfur ylides with α,β-unsaturated carbonyl compounds, has recently aroused renewed interest after the discovery of efficient catalysis by using (S)-indoline-2-carboxylic acid. In order to rationalize the behavior of this catalyst, MacMillan proposed a directed electrostatic activation (DEA) mechanism, in which the negative carboxylate group interacts with the positive thionium moiety, thus reducing the activation energy and increasing the reaction rate. More recently, Mayr refuted some of MacMillan conclusions, but accepted the DEA mechanism as a justification for the experimental high reaction rates. In contrast, our results indicate that the selectivity obtained in the process seems to result from several strong hydrogen bond interactions between the two reacting species, while no strong evidence for a DEA mechanism was found. We also concluded that the hydrogen bonds don't improve the reaction rate by lowering the activation energy of the rate-determining step, but can do it by promoting efficient reaction trajectories due to long-range complexation of the reagents. Finally, our results confirm that the cyclopropanation reaction occurs by a two-step mechanism, and that the overall enantioselectivity depends on the relative energies of the two steps, averaged by the relative populations of the iminium intermediates that are initially formed in the reaction.

Complexes of (η6-benzene)ruthenium(ii) with 1,4-bis(phenylthio/seleno-methyl)-1,2,3-triazoles: synthesis, structure and applications in catalytic activation of oxidation and transfer hydrogenation

1,4-Bis(phenylthio/seleno methyl)-1,2,3-triazoles (L1–L4) synthesized by a ‘Click’ reaction form complexes [(η⁶-C₆H₆)RuClL]PF₆ (L = L1–L4) suitable as catalyst for oxidation/transfer hydrogenation.

Reactivity of aldehydes at the air–water interface. Insights from molecular dynamics simulations and ab initio calculations

Computer simulations show that solvation effects at the air–water interface significantly influence the chemistry of aldehydes, enhancing for instance the benzaldehyde photolysis rate constant by one order of magnitude.

A computational mechanistic study of pH-dependent alcohol dehydrogenation catalyzed by a novel [C,N] or [C,C] cyclometalated Cp*Ir complex in aqueous solution

Our results have elucidated the impressive performance of a catalyst in alcohol dehydrogenation in aqueous solution by tuning the pH value.

Toward understanding the nature of internal rotation barriers with a new energy partition scheme: ethane and n-butane

On the basis of an alternative energy partition scheme where density-based quantification of the steric effect was proposed [Liu, S. B. J. Chem. Phys. 2007, 126, 244103], the origin of the internal rotation barrier between the eclipsed and staggered conformers of ethane and n-butane is systematically investigated in this work. Within the new scheme, the total electronic energy is decomposed into three independent components, steric, electrostatic, and fermionic quantum. The steric energy defined in this way is repulsive, exclusive, and extensive and intrinsically linked to Bader's atoms in molecules approach. Two kinds of differences, adiabatic (with optimal structure) and vertical (with fixed geometry), are considered for the molecules in this work. We find that in the adiabatic case the eclipsed conformer possesses a larger steric repulsion than the staggered conformer for both molecules, but in the vertical cases the staggered conformer retains a larger steric repulsion. For ethane, a linear relationship between the total energy difference and the fermionic quantum energy difference is discovered. This linear relationship, however, does not hold for n-butane, whose behaviors in energy component differences are found to be more complicated. The impact of basis set and density functional choices on energy components from the new energy partition scheme has been investigated, as has its comparison with another definition of the steric effect in the literature in terms of the natural bond orbital analysis through the Pauli Exclusion Principle. In addition, profiles of conceptual density functional theory reactivity indices as a function of dihedral angle changes have been examined. Put together, these results suggest that the new energy partition scheme provides insights from a different perspective of internal rotation barriers.

Structural and Dynamic Properties of (SiO2)6 Silica Nanostructures: A Quantum Molecular Dynamics Study

Structural and dynamic properties of the building block of silica nanowires, (SiO2)6, are investigated by Born-Oppenheimer quantum molecular dynamics simulations. Thirteen conformers have been identified, seven of which have not been reported before. The energy component analysis shows that the lower electrostatic interaction differentiates the global minimum from the other structures. We also observe that the maximum hardness principle can be employed to justify the molecular stability for this system. Time profiles of a few density functional reactivity indices exhibit correlations of dynamic fluctuations between HOMO and LUMO and between chemical potential and hardness. Electrophilicity, nucleaofugality, and electrofugality indices are found to change concurrently and significantly, indicating that the nanostructures sampled during the dynamic process are exceedingly reactive and rich in chemistry.