Mechanistic investigation of the uncatalyzed esterification reaction of acetic acid and acid halides with methanol: a DFT study

Linear and Nonlinear Optical Responses of Nitrobenzofurazan-Sulfide Derivatives: DFT-QTAIM Investigation on Twisted Intramolecular Charge Transfer

In this study, we report on the green fluorescence exhibited by nitrobenzofurazan-sulfide derivatives (NBD-Si, i = 1-4). The optical responses of these studied compounds in a polar methanol solvent were simulated by the use of time-dependent density functional theory (TD-DFT) employing the Becke-3-Parameter-Lee-Yang-Parr (B3LYP) functional along with the 6-31G(d,p) basis set. The computed energy and oscillator strength (f) results complement the experimental results. The band gap was calculated as the difference between the lowest unoccupied molecular orbital (LUMO) and the highest occupied molecular orbital (HOMO). Additionally, the density of states (DOS) was computed, providing a comprehensive understanding of the fundamental properties of these materials and further corroborating the experimental data. When the experimental data derived from ultraviolet/visible (UV/visible) and fluorescence spectroscopic techniques and those from simulated spectra are analyzed, the extracted values match up adequately. In addition, the NBD-sulfide compounds exhibit a large Stokes shift up to 85 nm in a polar methanol solvent. They are hypothesized to represent a novel paradigm of excited-state intramolecular charge transfer (ICT). To understand the intrinsic optical properties of NBD-Si materials, an ICT was identified, and its direction within the molecule was evaluated using the ratio of βvect and βtotal, values extracted from the computed nonlinear optical (NLO) properties. Moreover, the reduced density gradient (RDG)-based noncovalent interactions (NCIs) were employed to characterize the strength and type of NBD-Si interactions. Furthermore, noncovalent interactions were identified and categorized using the Quantum Theory of Atoms in Molecules (QTAIM) analysis. Ultimately, the combination of Hirshfeld surface analysis and DFT calculations was utilized to enhance the characterization and rationalization of these NCIs.

Formic acid-methanol complexation vs. esterification: elusive pre-reactive species identified by vibrational spectroscopy

Vibrational spectra of the mixed dimer and the two mixed trimers of methanol and formic acid as well as some of their isotopologues are presented. Out of the eight expected OH stretching fundamentals of these three pre-reactive hydrogen-bonded complexes, the three modes mainly involving an acid OH group bound to the alcohol appear to be missing in the jet-cooled spectra despite a combination of infrared and Raman probing. A possibility of spectral overlap is discussed in the mixed dimer case, but largely discarded. The missing modes correspond to (fractional) concerted elongation of all engaged OH bonds, promoting synchronous degenerate proton transfer between the molecules. One other trimer mode is very tentatively attributed to a broad spectral feature, whereas all OH bonds contacting carbonyl groups can be unambiguously identified by four relatively narrow infrared absorptions. The spectral features are confirmed by vibrational perturbation theory and deviate in a subtle but systematic way from scaled harmonic predictions which were previously validated for the formic acid complex with a more acidic alcohol. Despite being exothermic and exergonic, ester formation can only be detected in the rarefied gas expansions after extended pre-mixture of the gases, which somewhat contrasts the recent microwave spectroscopic evidence of in situ ester formation and in particular the lack of pre-reactive complex signals.

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

CONTEXT

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

METHODS

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

Ni2+ and Cu2+ complexes of N-(2,6-dichlorophenyl)-N-mesityl formamidine dithiocarbamate structural and functional properties as CYP3A4 potential substrates

Metal compounds continued to attract diverse applications due to their malleability in several capacities. In this study, we present our findings on the crystal structures and functional properties of Ni2+ and Cu2+ complexes of N'-(2,6-dichlorophenyl)-N-mesitylformamidine dithiocarbamate (L) comprising [Ni-(L)2] (1) and [Cu-(L)2] (2) with a four-coordinate metal center. We established the two complex structures through 1H and 13C nuclear magnetic resonance (NMR), elemental, and single-crystal X-ray analysis. The analyses showed that the two complexes are isomorphous, having P21/c as a space group and a unit-cell similarity index (π) of 0.002. The two complexes conform to a distorted square planar geometry around the metal centers. The calculated and experimental data, including bond lengths, angles, and NMR values, are similar. Hirshfeld surface analysis revealed the variational contribution of the different types of intermolecular contacts driven by the crystal lattice of the two solvated complexes. Our knowledge of the potential biological implication of these structures enabled us to probe the compounds as prospective CYP3A4 inhibitors. This approach mimics current trends in pharmaceutical design and biomedicine by incorporating potentially active molecules into various media to predict their biological efficacies. The simulations show appreciable binding of compounds 1 and 2 to CYP3A4 with average interaction energies of -97 and -87 kcal/mol, respectively. The protein attains at least five conformational states in the three studied models using a Gaussian Mixture Model-based clustering and free energy prediction. Electric field analysis shows the crucial residues to substrate binding at the active site, enabling CYP3A4 structure to function prediction. The predicted inhibition with these Ni2+ and Cu2+ complexes indicates that CYP3A4 overexpression in a diseased state like cancer would reduce, thereby increasing the chemotherapeutic compounds' shelf-lives for adsorption. This multidimensional study addresses various aspects of molecular metal electronics, including their application as substrate-mimicking inhibitors. The outcome would enable further research on bio-metal compounds of critical potential.

Quantitative In Silico Prediction of the Rate of Protodeboronation by a Mechanistic Density Functional Theory-Aided Algorithm

Computational reaction prediction has become a ubiquitous task in chemistry due to the potential value accurate predictions can bring to chemists. Boronic acids are widely used in industry; however, understanding how to avoid the protodeboronation side reaction remains a challenge. We have developed an algorithm for in silico prediction of the rate of protodeboronation of boronic acids. A general mechanistic model devised through kinetic studies of protodeboronation was found in the literature and forms the foundation on which the algorithm presented in this work is built. Protodeboronation proceeds through 7 distinct pathways, though for any particular boronic acid, only a subset of mechanistic pathways are active. The rate of each active mechanistic pathway is linearly correlated with its characteristic energy difference, which in turn can be determined using Density Functional Theory. We validated the algorithm using leave-one-out cross-validation on a data set of 50 boronic acids and made a further 50 rate predictions on academically and industrially important boronic acids out of sample. We believe this work will provide great assistance to chemists performing reactions that feature boronic acids, such as Suzuki-Miyaura and Chan-Evans-Lam couplings.

Coupled proton vibrations between two weak acids: the hinge complex between formic acid and trifluoroethanol

Raman and FTIR spectra of an acid–alcohol complex show complementary signatures from acidic and alcoholic OH stretching, proving its existence.

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.

A Theoretical and Experimental Study on Esterification of Citric Acid with the Primary Alcohols and the Hydroxyl Groups of Cellulose Chain (n = 1-2) in Parched Condition

Esterification of citric acid (CA) with the primary alcohols and the hydroxyl groups of cellulose chain (n = 1-2) in parched condition were investigated by using density functional theory (DFT) method and a two-layer ONIOM approach. Geometry and energy of reactants, products, and transition state (TS) structures were optimized at B3LYP/6-311g (d, p) level and ONIOM (B3LYP/6-311g (d, p):PM3MM) level. The computational results show that the esterification occurs in the two main steps: the first step is the dehydration reaction of CA to form anhydrides of 5-membered ring and 6-membered ring and the second step is the ring opening reaction with the hydroxyl –OH groups to form the ester products. The energy barrier of dehydration reaction step is much higher than that of ring opening reaction step. Effect of substituent R in primary alcohol R-CH2OH (R: CH=CH2, CH2NHCH3, CH2OCH3, CH2Cl) and cellulose chain (1G, 2G) on the reactivity, which has negative inductive effect –I, is significant. The combination of calculation data and experiment data were applied to make findings more rigorous. The activation energy of CA was determined by using differential scanning calorimetry (DSC) and thermal gravimetric (TG) analysis to be $E_a^{\exp }$  = 47.8 kcal/mol; the experimental data favoured the dehydration reaction step of CA.

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.

Does the Reaction of Cyclopropyl Acid Chlorides and Imines To Form 1,3-Oxazin-4-enone Heterocycles Proceed via a Ketene or an N-Acyl-iminium Mechanism?

The mechanism of the reaction between cyclopropyl acid chlorides and imines to form 1,3-oxazin-4-enones was probed through physical and computational experiments. The data gathered strongly support the reaction proceeding through an N-acyl iminium intermediate mechanism rather than a ketene intermediate mechanism.

Understanding the Hsp90 N-terminal Dynamics: Structural and Molecular Insights into the Therapeutic Activities of Anticancer Inhibitors Radicicol (RD) and Radicicol Derivative (NVP-YUA922)

Heat shock protein 90 (Hsp90) is a crucial component in carcinogenesis and serves as a molecular chaperone that facilitates protein maturation whilst protecting cells against temperature-induced stress. The function of Hsp90 is highly dependent on adenosine triphosphate (ATP) binding to the N-terminal domain of the protein. Thus, inhibition through displacement of ATP by means of competitive binding with a suitable organic molecule is considered an attractive topic in cancer research. Radicicol (RD) and its derivative, resorcinylic isoxazole amine NVP-AUY922 (NVP), have shown promising pharmacodynamics against Hsp90 activity. To date, the underlying binding mechanism of RD and NVP has not yet been investigated. In this study, we provide a comprehensive understanding of the binding mechanism of RD and NVP, from an atomistic perspective. Density functional theory (DFT) calculations enabled the analyses of the compounds' electronic properties and results obtained proved to be significant in which NVP was predicted to be more favorable with solvation free energy value of -23.3 kcal/mol and highest stability energy of 75.5 kcal/mol for a major atomic delocalization. Molecular dynamic (MD) analysis revealed NVP bound to Hsp90 (NT-NVP) is more stable in comparison to RD (NT-RD). The Hsp90 protein exhibited a greater binding affinity for NT-NVP (-49.4 ± 3.9 kcal/mol) relative to NT-RD (-28.9 ± 4.5 kcal/mol). The key residues influential in this interaction are Gly 97, Asp 93 and Thr 184. These findings provide valuable insights into the Hsp90 dynamics and will serve as a guide for the design of potent novel inhibitors for cancer treatment.

Concerted hydrolysis mechanism of HIV-1 natural substrate against subtypes B and C-SA PR: insight through molecular dynamics and hybrid QM/MM studies

It is well known that understanding the catalytic mechanism of HIV-1 PR is the rationale on which its inhibitors were developed; therefore, a better understanding of the mechanism of natural substrate hydrolysis is important. Herein, the reaction mechanism of HIV-1 natural substrates with subtypes B and common mutant in South Africa (subtype C-SA) protease were studied through transition state modelling, using a general acid-general base (GA-GB) one-step concerted process. The activation free energies of enzyme-substrate complexes were compared based on their rate of hydrolysis using a two-layered ONIOM (B3LYP/6-31++G(d,p):AMBER) method. We expanded our computational model to obtain a better understanding of the mechanism of hydrolysis as well as how the enzyme recognises or chooses the cleavage site of the scissile bonds. Using this model, a potential substrate-based inhibitor could be developed with better potency. The calculated activation energies of natural substrates in our previous study correlated well with experimental data. A similar trend was observed for the Gag and Gag-Pol natural substrates in the present work for both enzyme complexes except for the PR-RT substrate. Natural bond orbital (NBO) analysis was also applied to determine the extent of charge transfer within the QM part of both enzymes considered and the PR-RT natural substrate. The result of this study shows that the method can be utilized as a dependable computational technique to rationalize lead compounds against specific targets.

A comparison between observed and DFT calculations on structure of 5-(4-chlorophenyl)-2-amino-1,3,4-thiadiazole

The crystal and molecular structure of 5-(4-chlorophenyl)-2-amino-1,3,4-thiadiazole 3 was reported, which was characterized by various spectroscopic techniques (FT-IR, NMR and HRMS) and single-crystal X-ray diffraction. The crystal structure 3 (C₈H₆ClN₃S) crystallized in the orthorhombic space group Pna2₁ and the unit cell consisted of 8 asymmetric molecules. The unit cell parameters were a = 11.2027(2) Å, b = 7.6705(2) Å, c = 21.2166(6) Å, α = β = γ = 90°, V = 1823.15(8) ų, Z = 8. In addition, the structural geometry (bond lengths, bond angles, and torsion angles), the electronic properties of mono and dimeric forms of compound 3 were calculated by using the density functional theory (DFT) method at B3LYP level 6-31+ G(d,p), 6-31++ G(d,p) and 6-311+ G(d,p) basis sets in ground state. A good correlation was found (R² = 0.998) between the observed and theoretical vibrational frequencies. Frontier molecular orbitals (HOMO and LUMO) and Molecular Electrostatic Potential map of the compound was produced by using the optimized structures. The NBO analysis was suggested that the molecular system contains N-H…N hydrogen bonding, strong conjugative interactions and the molecule become more polarized owing to the movement of π-electron cloud from donor to acceptor. The calculated structural and geometrical results were in good rational agreement with the experimental X-ray crystal structure data of 1,3,4-thiadiazol-2-amine, 3. The compound 3 exhibited n→π* UV absorption peak of UV cutoff edge, and great magnitude of the first-order hyperpolarizability was observed. The obtained results suggest that compound 3 could have potential application as NLO material. Therefore, this study provides valuable insight experimentally and theoretically, for designing new chemical entities to meet the demands of specific applications.

The catalytic role of water in the binding site of l,d-transpeptidase 2 within acylation mechanism: A QM/MM (ONIOM) modelling

Mycobacterium tuberculosis is the causative agent of Tuberculosis. Formation of 3 → 3 crosslinks in the peptidoglycan layer of M. tuberculosis is catalyzed by l,d-transpeptidases. These enzymes can confer resistance against classical β-lactams that inhibit enzymes that generate 4 → 3 peptidoglycan crosslinks. The focus of this study is to investigate the catalytic role of water molecules in the acylation mechanism of the β-lactam ring within two models; 4- and 6-membered ring systems using two-layered our Own N-layer integrated Molecular Mechanics ONIOM (B3LYP/6-311++G(2d,2p): AMBER) model. The obtained thermochemical parameters revealed that the 6-membered ring model best describes the inhibition mechanism of acylation which indicates the role of water in the preference of 6-membered ring reaction pathway. This finding is in accordance with experimental data for the rate-limiting step of cysteine protease with the same class of inhibitor and binding affinity for both inhibitors. As expected, the ΔG^# results also reveal that the 6-membered ring reaction pathway is the most favourable. The electrostatic potential (ESP) and the natural bond orbital analysis (NBO) showed stronger interactions in 6-membered ring transition state (TS-6) mechanism involving water in the active site of the enzyme. This study could be helpful in the development of novel antibiotics against l,d-transpeptidase.

Applicability of DFT model in reactive distillation

The density functional theory (DFT) applicability to reactive distillation is discussed. Brief modeling techniques description of distillation and rectification with chemical reaction is provided as a background for quantum method usage description. The equilibrium and nonequilibrium distillation models are described for that purpose. The DFT quantum theory is concisely described. The usage of DFT in the modeling of reactive distillation is described in two parts. One of the fundamental and very important component of distillation modeling is vapor-liquid equilibrium description for which the DFT quantum approach can be used. The representative DFT models, namely COSMO-RS (Conductor like Screening Model for Real Solvents), COSMOSPACE (COSMO Surface Pair Activity Coefficient) and COSMO-SAC (SAC – segment activity coefficient) approaches are described. The second part treats the way in which the chemical reaction is described by means of quantum DFT method. The intrinsic reaction coordinate (IRC) method is described which is used to find minimum energy path of substrates to products transition. The DFT is one of the methods which can be used for that purpose. The literature data examples are provided which proves that IRC method is applicable for chemical reaction kinetics description.

Computational Investigation of the Aminolysis of RAFT Macromolecules

This work presents a detailed computational study and kinetic analysis of the aminolysis of dithioates, dithiobenzoates, trithiocarbonates, xanthates, and thiocarbamates, which are frequently used as chain-transfer agents for reversible addition-fragmentation chain-transfer (RAFT) polymerization. Rate coefficients were obtained from ab initio calculations, taking into account a diffusional contribution according to the encounter pair model. A kinetic model was constructed and reveals a reaction mechanism of four elementary steps: (i) formation of a zwitterionic intermediate, (ii) formation of a complex intermediate in which an assisting amine molecule takes over the proton from the zwitterionic intermediate, (iii) breakdown of the complex into a neutral tetrahedral intermediate with release of the assisting amine molecule, and (iv) amine-assisted breakdown of the neutral intermediate to the products. Furthermore, a comparative analysis indicates that the alkanedithioates and dithiobenzoates react the fastest, followed, respectively, by xanthates and trithiocarbonates, which react almost equally fast, and dithiocarbamates, which are not reactive at typical experimentally used conditions.