Effect of solvation on the condensed Fukui function and the generalized philicity index

DFT, ADME studies and evaluation of the binding with HSA and MAO-B inhibitory potential of protoberberine alkaloids from Guatteria friesiana: theoretical insights of promising candidates for the treatment of Parkinson's disease

CONTEXT

Parkinson's disease is a chronic neurodegenerative condition that has no cure, characterized by the progressive degeneration of specific brain cells responsible for producing dopamine, a crucial neurotransmitter for controlling movement and muscle coordination. Parkinson's disease is estimated to affect around 1% of the world's population over the age of 60, but it can be diagnosed at younger ages. One of the treatment strategies for Parkinson's disease involves the use of drugs that aim to increase dopamine levels or simulate the action of dopamine in the brain. A class of commonly prescribed drugs are the so-called monoamine oxidase B (MAO-B) inhibitors due to the fact that this enzyme is responsible for metabolizing dopamine, thus reducing its levels in the brain. Studies have shown that berberine-derived alkaloids have the ability to selectively inhibit MAO-B activity, resulting in increased dopamine availability in the brain. In this context, berberine derivatives 13-hydroxy-discretinine and 7,8-dihydro-8-hydroxypalmatine, isolated from Guatteria friesiana, were evaluated via density functional theory followed by ADME studies, docking and molecular dynamic simulations with MAO-B, aiming to evaluate their anti-Parkinson potential, which have not been reported yet. Docking simulations with HSA were carried out aiming to evaluate the transport of these molecules through the circulatory system.

METHODS

The 3D structures of the berberine-derived alkaloids were modeled via the DFT approach at B3LYP-D3(BJ)/6-311 + + G(2df, 2pd) theory level using Gaussian 09 software. Solvation free energies were determined through Truhlar's solvation model. MEP and ALIE maps were generated with Multiwfn software. Autodock Vina software was used for molecular docking simulations and analysis of the interactions in the binding sites. The 3D structure of MAO-B was obtained from the Protein Data Bank website under PDB code 2V5Z. For the interaction of studied alkaloids with human serum albumin (HSA) drug sites, 3D structures with PDB codes 2BXD, 2BXG, and 4L9K were used. Molecular dynamics simulations were carried out using GROMACS 2019.4 software, with the GROMOS 53A6 force field at 100 ns simulation time. The estimation of the ligand's binding free energies was obtained via molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method.

Electrophilicity index revisited

This review aims to be a comprehensive, authoritative, critical, and accessible review of general interest to the chemistry community; because the electrophilicity index is a very useful global reactivity descriptor defined within a conceptual density functional theory framework. Our group has also introduced electrophilicity based new global and local reactivity descriptors and also new associated electronic structure principles, which are important indicators of structure, stability, bonding, reactivity, interactions, and dynamics in a wide variety of physico-chemical systems and processes. This index along with its local counterpart augmented by the associated electronic structure principles could properly explain molecular vibrations, internal rotations and various types of chemical reactions. The concept of the electrophilicity index has been extended to dynamical processes, excited states, confined environment, spin-dependent and temperature-dependent situations, biological activity, site selectivity, aromaticity, charge removal and acceptance, presence of external perturbation through solvents, external electric and magnetic fields, and so forth. Although electrophilicity and its local variant can adequately interpret the behavior of a wide variety of systems and different physico-chemical processes involving them, their predictive potential remains to be explored. An exhaustive review on all these aspects will set the tone of the future research in that direction.

Regioselective Cycloaddition of Nitrile Imines to 5-Methylidene-3-phenyl-hydantoin: Synthesis and DFT Calculations

Nitrile imine cycloaddition to hydantoins containing an exocyclic C=C double bond has been previously described in a very limited number of examples. In this work, regioselective synthesis of spiro-pyrazoline-imidazolidine-2,4-diones based on a 1,3-dipolar cycloaddition reaction of nitrile imines to 5-methylidene-3-phenyl-hydantoin have been proposed. It was found that, regardless of the nature of the aryl substituents at the terminal C and N atoms of the C-N-N fragment of nitrile imine (electron donor or electron acceptor), cycloaddition to the 5-methylidenhydantoin exocyclic C=C bond proceeds regioselectively, and the terminal nitrogen atom of the nitrile imine connects to the more sterically hindered carbon atom of the double bond, which leads to the formation of a 5-disubstituted pyrazoline ring. The observed cycloaddition regioselectivity was rationalized using DFT calculations of frontier molecular orbital interactions, global CDFT reactivity indices, and minimum energy paths.

Is It Possible to Restrain OER on Simple Carbon Electrodes to Efficiently Electrooxidize Organic Pollutants?

This paper presents a comparative analysis of three carbon-based electrodes: bare multiwalled carbon nanotubes (MWCNT), SnO2/MWCNT, and PbO2/graphene-nanoribbons (PbO2/GNR) composites, as anodes for the electrooxidative degradation of Rhodamine B as a model organic pollutant. Anodic electrooxidation of Rhodamine B was performed on all three electrodes, and the decolorization efficiency was found to increase in the order MWCNT < PbO2/GNR < SnO2/MWCNT. The electrodes were characterized by X-ray photoelectron spectroscopy (XPS) and linear sweep voltammetry (LSV). It was proposed that, in the 0.1 M Na2SO4 applied as electrolyte, observed decolorization mainly occurs in the interaction of Rhodamine B with OH radical adsorbed on the anode. Finally, the obtained results were complemented with Density Functional Theory (DFT) calculations of OH-radical interaction with appropriate model surfaces: graphene(0001), SnO2(001), and PbO2(001). It was found that the stabilization of adsorbed OH-radical on metal oxide spots (SnO2 or PbO2) compared to carbon is responsible for the improved efficiency of composites in the degradation of Rhodamine B. The observed ability of metal oxides to improve the electrooxidative potential of carbon towards organic compounds can be useful in the future design of appropriate anodes.

XNgNSi (X = HCC, F; Ng = Kr, Xe, Rn): A New Class of Metastable Insertion Compounds Containing Ng-C/F and Ng-N Bonds and Possible Isomerization therein

Recently, astronomically important silaisocyanoacetylene (HCCNSi) possessing a large dipole moment has been detected for the first time with the help of crossed molecular beam experiments. Quantum chemical computations at higher levels of theory have also been performed to characterize the transient species. In this study, we have analyzed the equilibrium geometry, stability, reactivity, and energetics as well as the nature of bonding in the noble gas (Ng) inserted HCCNSi compound. We have also considered its F analogue to understand the influence of the most electronegative atom in the compound. Metastable behavior of the XNgNSi compounds (X = HCC, F; Ng = Kr-Rn) is examined by calculating thermochemical parameters like free energy change (ΔG) and zero-point-energy-corrected dissociation energy (D₀) at 298 K for all possible two-body (2B) and three-body (3B) (both neutral as well as ionic) dissociation channels using coupled-cluster theory [CCSD(T)] in addition to density functional theory (DFT) as well as second order Møller-Plesset perturbation theory (MP2). The set of predicted compounds is found to be endergonic in nature, having high positive free energy change suggesting the thermochemical stability of the compounds except for the 2B Ng-release paths. Though thermodynamically feasible, they are kinetically protected with very high activation free energy barriers. Interestingly, the release of Ng from the parent moiety XNgNSi produces the XSiN isomer, by 180° flipping of the NSi moiety. This can also be seen in the dynamical simulation carried out with the help of atom-centered density matrix propagation (ADMP) technique at 2000K for 1 ps. The bonding in Ng-C, Ng-F, and Ng-N bonds of the studied compounds is analyzed and described with the aid of natural bond orbital (NBO), topological parameters computed using atoms-in-molecules theory (AIM), energy decomposition analysis (EDA), and adaptive natural density partitioning (AdNDP) methods. The natural charge distribution on the constituent atoms suggests that the compounds can be partitioned into both ways of representations, viz., neutral radical as well as ionic fragments. Lastly, the reactivity of the compounds is scrutinized using certain reactivity descriptors calculated within the domain of conceptual density functional theory (CDFT).

Synthesis and Assessment of Two Malonyl Dihydrazide Derivatives as Corrosion Inhibitors for Carbon Steel in Acidic Media: Experimental and Theoretical Studies

Despite the extensive use of carbon steel in all industrial sectors, particularly in the petroleum industry, its low corrosion resistance is an ongoing problem for these industries. In the current work, two malonyl dihydrazide derivatives, namely 2,2'-malonylbis (N-phenylhydrazine-1-carbothiamide (MBC) and N'1, N'3-bis(-2-hydroxybenzylidene) malonohydrazide (HBM), were examined as inhibitors for the carbon steel corrosion in 1.0 M HCl. Both MBC and HBM were characterised using thin-layer chromatography, elemental analysis, infrared spectroscopy, and nuclear magnetic resonance techniques. The corrosion tests were performed using mass loss measurements, polarisation curves, and electrochemical impedance spectroscopy. It is obtained from the mass loss studies that the optimal concentration for both inhibitors is 2.0 × 10^-5 mol/L, and the inhibition efficiencies reached up to 90.7% and 84.5% for MBC and HBM, respectively. Electrochemical impedance spectroscopy (EIS) and potentiodynamic polarisation (PDP) indicate an increased impedance in the presence of both MBC and HBM and mixed-type inhibitors, respectively. Both inhibitors can mitigate corrosion in the range of 298-328 K. Values of free energy changes obtained from the Langmuir model suggest that the inhibitors suppress the corrosion process principally by chemisorption. The computational investigations were conducted to identify the factors connected with the anti-corrosive properties of the examined inhibitors.

Conceptual density functional theory based electronic structure principles

In this review article, we intend to highlight the basic electronic structure principles and various reactivity descriptors as defined within the premise of conceptual density functional theory (CDFT). Over the past several decades, CDFT has proven its worth in providing valuable insights into various static as well as time-dependent physicochemical problems. Herein, having briefly outlined the basics of CDFT, we describe various situations where CDFT based reactivity theory could be employed in order to gain insights into the underlying mechanism of several chemical processes.

Structural characterization and QSAR modeling of 1,2,4-triazole derivatives as α-glucosidase inhibitors

In order to identify potential new drugs classified as alpha-Glucosidase inhibitors (AGIs), we used quantum chemical descriptors and QSAR modeling to predict the biological activity of triazole derivatives as AGIs.

Exploring the antioxidant activity of thiaflavan compounds: a quantum chemical study

Conceptual DFT tools (HAT, SPLET, SET-PT, aromaticity index,…) have been used to explore the antioxidant activity of thiaflavan compounds, and predict which derivate should be the best one.

DFT studies on the diastereoselectivity and regioselectivity of multicomponent domino Knoevenagel/Diels–Alder reaction

In this study, mechanism and stereochemistry of multicomponent domino Knoevenagel/Diels–Alder reaction were investigated theoretically. Structures of reagents, transition states, intermediates, and products were optimized at M062X/6-31[Formula: see text]G(d,p) level of theory. Although the mechanism of this reaction involved several processes, including bond rotation, isomerization, asymmetric cycloaddition, acid-base, and nucleophile–electrophile competitions, critical processes were studied to provide a clearer picture of the mechanism of this valuable reaction in terms of stereochemistry considerations. According to the results, the ring closure step of reaction performed via a polar Dield-Alder mechanism, having enthalpy at approximately 40[Formula: see text]kcal/mol. The diastereoselectivity of the reaction was controlled by the interaction of dienophile with the less hindered face of diene through a more stable endo transition state ([Formula: see text] and 27.31 in methanol and gas phase, respectively). HSAB criteria could explain the regioselectivity of this reaction by considering the least softness difference ([Formula: see text]) for interacting C-atoms based on Hirshfeld populations. The result was the creation of cis-spiro cyclohexanone under kinetic and thermodynamic controls as a predominant diastereoselective and regioselective product.

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.

Substituent effects on the stability, physicochemical properties and chemical reactivity of nitroimidazole derivatives with potential antiparasitic effect: a computational study

Neglected tropical diseases caused by parasitic pathogens have caused an increase in research interest in drug discovery.

SEI-forming electrolyte additives for lithium-ion batteries: development and benchmarking of computational approaches

SEI-forming additives play an important role in lithium-ion batteries, and the key to improving battery functionality is to determine if, how, and when these additives are reduced. Here, we tested a number of computational approaches and methods to determine the best way to predict and describe the properties of the additives. A wide selection of factors were evaluated, including the influences of the solvent and lithium cation as well as the DFT functional and basis set used. An optimized computational methodology was employed to assess the usefulness of different descriptors.

Quantitative Structure-Property Relationship (QSPR) Models for a Local Quantum Descriptor: Investigation of the 4- and 3-Substituted-Cinnamic Acid Esterification

In this work, the theoretical description of the 4- and 3-substituted-cinnamic acid esterification with different electron donating and electron withdrawing groups was performed at the B3LYP and M06-2X levels, as a two-step process: the O-protonation and the nucleophile attack by ethanol. In parallel, an experimental work devoted to the synthesis and characterization of the substituted-cinnamate esters has also been performed. In order to quantify the substituents effects, quantitative structure–property relationship (QSPR) models based on the atomic charges, Fukui functions and the Frontier Effective-for-Reaction Molecular Orbitals (FERMO) energies were investigated. In fact, the Fukui functions, ƒ⁺C and ƒ⁻O, indicated poor correlations for each individual step, and in contrast with the general literature, the O-protonation step is affected both by the FERMO energies and the O-charges of the carbonyl group. Since the process was shown to not be totally described by either charge- or frontier-orbitals, it is proposed to be frontier-charge-miscere controlled. Moreover, the observed trend for the experimental reaction yields suggests that the electron withdrawing groups favor the reaction and the same was observed for Step 2, which can thus be pointed out as the determining step.

Chemical Reactivity of Isoproturon, Diuron, Linuron, and Chlorotoluron Herbicides in Aqueous Phase: A Theoretical Quantum Study Employing Global and Local Reactivity Descriptors

We have calculated global and local DFT reactivity descriptors for isoproturon, diuron, linuron, and chlorotoluron herbicides at the MP2/6-311++G(2d,2p) level of theory. The results suggest that, in aqueous conditions, chlorotoluron, linuron, and diuron herbicides may be degraded by elimination of urea moiety through electrophilic attacks. On the other hand, electrophilic, nucleophilic, and free radical attacks on isoproturon may cause the elimination of isopropyl fragment.

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.

Condensed-to-atoms electronic Fukui functions within the framework of spin-polarized density-functional theory

A simple formalism devised to calculate the condensed-to-atoms Fukui function [R. R. Contreras, P. Fuentealba, M. Galván, and P. Pérez, Chem. Phys. Lett. 304, 405 (1999)] has been further analyzed within a spin-polarized density-functional theory framework. The model is based on a frozen-core approximation to these local reactivity indices [M. Galván, A. Vela, and J. L. Gázquez, J. Phys. Chem. 92, 6470 (1988)], giving us an extended local reactivity description of systems based on the frontier spin-up and spin-down molecular orbitals. Degenerate molecular spin orbitals have been explicitly included in our model equations. Computational results for the nitric oxide (NO) and some simple carbene systems are presented in order to test the model. These quantities have been discussed in the context of changes both in charge density and spin density within the context of electron charge transfer or spin-polarization processes.