On the nature of Parr functions to predict the most reactive sites along organic polar reactions

Revisiting nucleophilicity: an index for chemical reactivity from a CDFT approach

CONTEXT

Understanding and predicting the nucleophilic reactivity are paramount in elucidating organic chemical reactions and designing new synthetic pathways. In this study, we propose a nucleophilicity index within the framework of Conceptual Density Functional Theory (CDFT). Through rigorous theoretical formulations, we introduce an original quantum reactivity descriptor that captures the nucleophilic propensity of molecules based on their electronic structure and chemical environment. Subsequently, this proposed index is applied to a series of nucleophiles (pyrrolidines derivatives), spanning a diverse range of chemical functionalities. Our computational assessments reveal insightful correlations between the predicted nucleophilicity index and experimental observations of nucleophilic behavior. Thereby, they offer a promising avenue for advancing the understanding of organic reactivity and guiding synthetic efforts.

METHODS

Experimentally, Mayr's experimental parameters accounting for nucleophilicity were selected for the pyrrolidines. This study used DFT calculations at the B3LYP/Aug-cc-pVTZ level of theory using the Gaussian 16 program. Geometry optimization was thus performed, and the methodology employed for the computation of quantum reactivity descriptor is presented. Solvent effect was also taken into account using IEFPCM, and empirical dispersion correction (GD3) was employed.

Inhibition of corrosion of an aluminum alloy by rosemary and eucalyptus extracted oils in 1 M hydrochloric acid medium: an experimental and theoretical study

In this study, we explored aluminum corrosion inhibition field of study in a 1 M HCl solution, harnessing the power of essential oils extracted from rosemary and eucalyptus plants. Our exploration gives a comprehensive analysis of the pivotal factors that shape the corrosion inhibition process. Our scientific journey was marked by a deliberate and systematic approach, encompassing the utilization of gravimetric analysis (weight loss), electrochemical potentiodynamic polarization, and the sophisticated electrochemical impedance spectrometry (EIS) techniques. Our findings unveiled promising and nuanced outcomes, particularly in the area of the electrochemical technique. This method demonstrated remarkable inhibition efficiencies, ranging from 42% to an impressive 92% for rosemary essential oil and from 37 to 84% for eucalyptus essential oil. These results unveiled a dynamic relationship between essential oil concentration and inhibition efficiency, a revelation that further deepens our understanding of the corrosion inhibition process. The inhibition efficiency increased with higher concentrations of essential oil but decreased with elevated temperatures. Furthermore, our analysis traversed into the realms of potentiodynamic and thermodynamic insights. These analytical techniques unearthed the complex mechanisms at play, explaining the pathway followed by the studied inhibitors. They exhibited their prowess by forming protective films on the metal surface, acting as vigilant protectors against the relentless forces of corrosion. Complementing our experimental findings, our study of computational chemistry through density functional theory (DFT) unveiled remarkable insights. It elucidated the spontaneous adsorption process of inhibitor molecules onto the aluminum surface in the presence of H2O solvent. This computational harmony with our experimental results strengthened our confidence in the robustness of our findings. One of the key findings of this study was the superior inhibitory power of camphor in rosemary EO and β-myrcene in eucalyptus essential oil EO, respectively, attributed to the distinctive characteristics of the active sites found in each compound. The inhibitory effectiveness followed the order β-myrcene > camphor > borneol > α-pinene > bornyl acetate > p-cymene > 1,8-cineole. These compounds, notable for their distinct active sites, emerged as exceptional agents in the pursuit of effective corrosion inhibition.

A molecular electron density theory study on the Chichibabin reaction: The origin of regioselectivity

A Molecular Electron Density Theory (MEDT) study was performed on the nucleophilic amination of pyridine and benzene to elucidate the observed regioselectivity in the Chichibabin reaction. For this purpose, three possible reaction paths were considered between NaNH₂ and the 2-, 3- and 4-positions of pyridine. The reaction of NaNH₂ and benzene was also investigated. The results indicated that in each reaction, the first step involves the nucleophilic attack of the NH₂‾ ion toward the aromatic system to produce an anionic intermediate. The second step, proceeds via hydride elimination of the anionic intermediate to produce the corresponding aromatic amine. This step is more favourable both kinetically and thermodynamically for the reaction at the 2-position of pyridine relative to the 4-position. In contrast to the Parr functions analysis which indicates that the 4-position is activated more electrophilically in pyridine, the PES analysis reveals that the 2-position attack with the NH₂‾ ion is more favourable both thermodynamically and kinetically. The results for the reaction of sodium amide and benzene indicated that this reaction is not feasible due to the high activation barriers. In consistent with the experimental results, the ELF analysis indicated that in the first step of the reaction between pyridine and sodium amide, the formation of the C2-N bond takes place via a pseudoradical coupling between the C2 carbon atom of pyridine and the N atom of NaNH₂. These variations occur in the second-half of the first step of the reaction. NCI analysis on the reaction revealed that the more attractive forces between the fragments, make the transition states for the 2-position attack more stable relative to the others. Thus, the NCI analysis can satisfactorily describe the observed regioselectivity in the Chichibabin reaction.

Novel Quinazolinone–Isoxazoline Hybrids: Synthesis, Spectroscopic Characterization, and DFT Mechanistic Study

Quinazolinone and isoxazoline systems have attracted much attention due to their interesting pharmacological properties. The association of these two pharmacophores in a single hybrid structure can boost the biological activity or bring a new one. Inspired by this new paradigm, in the present work we report the synthesis and spectroscopic characterization of new quinazolinone–isoxazoline hybrids. The target compounds were obtained via 1,3-dipolar cycloaddition reactions of arylnitriloxides and N-allylquinazolinone. The synthesized compounds were characterized using spectroscopic techniques such as IR, 1D NMR (1H and 13C), 2D NMR (COSY and HSQC), and high-resolution mass spectrometry (HRMS). The spectral data show that this reaction leads only to the 3,5-disubstituted isoxazoline regioisomer, and that the observed regiochemistry is not affected by the nature of the substituents in the phenyl ring of the dipole. In addition, a theoretical study was performed using density functional theory (DFT) to support the experimental results in regard to the regiochemistry of the studied reactions. The computational mechanistic study was in good agreement with the experimental data.

Mesoporous Silica Modified with 2-Phenylimidazo[1,2-a] pyridine-3-carbaldehyde as an Effective Adsorbent for Cu(II) from Aqueous Solutions: A Combined Experimental and Theoretical Study

In this study, we will present an efficient and selective adsorbent for the removal of Cu(II) ions from aqueous solutions. The silica-based adsorbent is functionalized by 2-phenylimidazo[1,2-a] pyridine-3-carbaldehyde (SiN-imd-py) and the characterization was carried out by applying various techniques including FT-IR, SEM, TGA and elemental analysis. The SiN-imd-py adsorbent shows a good selectivity and high adsorption capacity towards Cu(II) and reached 100 mg/g at pH = 6 and T = 25 °C. This adsorption capacity is important compared to other similar adsorbents which are currently published. The adsorption mechanism, thermodynamics, reusability and the effect of different experimental conditions, such as contact time, pH and temperature, on the adsorption process, were also investigated. In addition, a theoretical study was carried out to understand the adsorption mechanism and the active sites of the adsorbent, as well as the stability of the complex formed and the nature of the bonds.

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⁻¹, 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 Ca_v3.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.

Understanding the Mechanism of Diels-Alder Reactions with Anionic Dienophiles: A Systematic Comparison of [ECX]- (E = P, As; X = O, S, Se) Anions

While Diels–Alder (DA) reactions involving neutral or cationic dienophiles are well-known, the characteristics of the analogous reactions with anionic dienophiles are practically unexplored. Herein we present the first comparative computational investigations on the characteristics of DA cycloadditions with anionic dienophiles on the basis of the reactions of [ECX]⁻ anions (E = P, As; X = O, S, Se) with 2H-pyran-2-one. All of these reactions were found to be both kinetically and thermodynamically feasible, enabling synthetic access toward 2-phosphaphenolate and arsaphenolate derivatives in the future. This study also reveals that the [ECO]⁻ anions show clear regioselectivity, while for [ECS]⁻ and [ECSe]⁻ anions, the two possible reaction channels have very similar energetics. Additionally, the activation barriers for the [ECO]⁻ anions are lower than those of the heavier analogues. The observed differences can be traced back to the starkly differing nucleophilic character of the pnictogen center in the anions, leading to a barrier-lowering effect in the case of the [ECO]⁻ anions. Furthermore, analysis of the geometries and electron distributions of the corresponding transition states revealed structure–property relationships, and thus a direct comparison of the cycloaddition reactivity of these anions was achieved. Along one of the two pathways, a good correlation was found between the activation barriers and suitable nucleophilicity descriptors (nucleophilic Parr function and global nucleophilicity). Additionally, the tendency of the reaction energies can be explained by the changing aromaticity of the products.

In contrast to the phosphaethynolate [PCO]⁻ anion, the cycloaddition reactivity of the heavier congeners ([ECX]⁻, where E = P, As and X = O, S, Se) is unexplored. In this computational study, the Diels−Alder reaction between the known [ECX]⁻ anions and 2-pyrone was employed to compare the reactivity patterns. The first activation barrier of these reactions correlates with the nucleophilicity of the anions, indicating a barrier-lowering effect. The feasibility of the studied reactions, leading to P and As heterocycles, was also explored.

Computational peptidology approach to the study of the chemical reactivity and bioactivity properties of Aspergillipeptide D, a cyclopentapeptide of marine origin

Aspergillipeptide D is a cyclic pentapeptide isolated from the marine gorgonian Melitodes squamata-derived fungus Aspergillus sp. SCSIO 41501 that it has been shown to present moderate activity against herpes virus simplex type 1 (HSV-1). Thus, this paper presents the results of a computational study of this cyclopentapeptide's chemical reactivity and bioactivity properties using a CDFT-based computational peptidology (CDFT-CP) methodology, which is derived from combining chemical reactivity descriptors derived from Conceptual Density Functional Theory (CDFT) and some Cheminformatics tools which may be used. This results in an improvement of the virtual screening procedure by a similarity search allowing the identification and validation of the known ability of the peptide to act as a possible useful drug. This was followed by an examination of the drug's bioactivity and pharmacokinetics indices in relation to the ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) characteristics. The findings provide further evidence of the MN12SX density functional's superiority in proving the Janak and Ionization Energy theorems using the proposed KID approach. This has proven to be beneficial in accurately predicting CDFT reactivity characteristics, which aid in the understanding of chemical reactivity. The Computational Pharmacokinetics study revealed the potential ability of Aspergillipeptide D as a therapeutic drug through the interaction with different target receptors. The ADMET indices confirm this assertion through the absence of toxicity and good absorption and distribution properties.

A DFT study of NHC-catalyzed reactions between 2-bromo-2-enals and acylhydrazones: mechanisms, and chemo- and stereoselectivities

The reaction mechanisms and origins of the chemo- and stereoselectivities of NHC-catalyzed [4 + 2] annulation of 2-bromo-2-enals and acylhydrazones.

Mechanism and regio- and stereoselectivity in an NHC-catalyzed Mannich/lactamization domino reaction

Density functional theory (DFT) calculations (M06-2X) have been employed to disclose the mechanisms and regio- and stereo-selectivities of the N-heterocyclic carbene (NHC)-catalyzed reaction of 2-benzothiazolimines and α-chloroaldehydes. The preferred mechanism is initiated by the nucleophilic attack of NHC on α-chloroaldehyde (first step), followed by 1,2-proton transfer which was assisted by the Brønsted acid DABCO·H+ to generate the Breslow intermediate (second step). The cleavage of the C-Cl bond (third step) and deprotonation (fourth step) form the enolate intermediate. This further reacts with 2-benzothiazolimine which leads to the formation of a new C-C bond (fifth step). Subsequent cyclization takes place via the formation of a new C-N bond (sixth step). Catalyst regeneration completes the whole catalytic cycle and affords the final product (seventh step). The DFT results indicate that the fifth step determines the stereochemistry of the reaction and leads to benzothiazolopyrimidinone with the SS configuration, which agrees well with experimental observations. Intramolecular cyclization is found to be the regioselectivity-determining step, for which the [4+2] annulation pathway is more preferred than that via [2+2] annulation, which again agrees well with experimental observations. Based on the mechanism proposed, the origins of regio- and stereoselectivities have also been investigated by performing distortion/interaction, natural bond orbital (NBO) and non-covalent interaction (NCI) analyses. The mechanistic insights gained in this work should be helpful in the rational design of potential catalysts for analogous reactions.

Structure-directed formation of the dative/covalent bonds in complexes with C70piperidine

The combined experimental-computational study has been performed to investigate the complexes formed between C70 carbon allotrope and piperidine. The results of FT-IR, H-NMR, and C-NMR measurements, together with the calculations based on the DFT approach and molecular dynamics simulations, prove the existence of dative/covalent bonding in C70piperidine complexes. The dative bond forms not only at the region of five- and six-membered rings, observed previously with C60, but also at the region formed of six-membered rings. The structure, i.e., nonplanarity, explains the observed dative bond formation. New findings on the character of interaction of secondary amines with C70 bring new aspects for the rational design of modified fullerenes and their applications in electrocatalysis, spintronics, and energy storage.

Unifying Conceptual Density Functional and Valence Bond Theory: The Hardness-Softness Conundrum Associated with Protonation Reactions and Uncovering Complementary Reactivity Modes

In this study, we address the long-standing issue-arising prominently from conceptual density functional theory (CDFT)-of the relative importance of electrostatic, i.e., "hard-hard", versus spin-pairing, i.e., "soft-soft", interactions in determining regiochemical preferences. We do so from a valence bond (VB) perspective and demonstrate that VB theory readily enables a clear-cut resolution of both of these contributions to the bond formation/breaking process. Our calculations indicate that appropriate local reactivity descriptors can be used to gauge the magnitude of both interactions individually, e.g., Fukui functions or HOMO/LUMO orbitals for the spin-pairing/(frontier) orbital interactions and molecular electrostatic potentials (and/or partial charges) for the electrostatic interactions. In contrast to previous reports, we find that protonation reactions cannot generally be classified as either charge- or frontier orbital-controlled; instead, our results indicate that these two bonding contributions generally interplay in more subtle patterns, only giving the impression of a clear-cut dichotomy. Finally, we demonstrate that important covalent, i.e., spin pairing, reactivity modes can be missed when only a single spin-pairing/orbital interaction descriptor is considered. This study constitutes an important step in the unification of CDFT and VB theory.

Virtual Screening of Marine Natural Compounds by Means of Chemoinformatics and CDFT-Based Computational Peptidology

This work presents the results of a computational study of the chemical reactivity and bioactivity properties of the members of the theopapuamides A-D family of marine peptides by making use of our proposed methodology named Computational Peptidology (CP) that has been successfully considered in previous studies of this kind of molecular system. CP allows for the determination of the global and local descriptors that come from Conceptual Density Functional Theory (CDFT) that can give an idea about the chemical reactivity properties of the marine natural products under study, which are expected to be related to their bioactivity. At the same time, the validity of the procedure based on the adoption of the KID (Koopmans In DFT) technique, as well as the MN12SX/Def2TZVP/H2O model chemistry is successfully verified. Together with several chemoinformatic tools that can be used to improve the process of virtual screening, some additional properties of these marine peptides are identified related to their ability to behave as useful drugs. With the further objective of analyzing their bioactivity, some useful parameters for future QSAR studies, their predicted biological targets, and the ADMET (Absorption, Distribution, Metabolism, Excretion and Toxicity) parameters related to the theopapuamides A-D pharmacokinetics are also reported.

Understanding the Mechanism and Selectivities of the Reaction of Meta-Chloroperbenzoic Acid and Dibromocarbene with β-Himachalene: A DFT Study

This study was performed to understand the site selectivity in the reaction between β-himachalene and meta-chloroperbenzoic acid (m-CPBA) in the first step followed by the addition of dibromocarbene (CBr2) to the main monoepoxidation product Pα formed in the first reaction. Calculations were performed using the Becke three-parameter hybrid exchange functional and the Lee–Yang–Parr correlation functional (B3LYP) with the 6-311 + G (d, p) basis set. Transition states were located by QST2, and their highlighting was validated by the existence of only one imaginary frequency in the Hessian matrix. The action of m-CPBA on β-himachalene was analyzed on the two double bonds of β-himachalene whose theoretical calculations show that the attack affects the most substituted double bond on α side containing hydrogen of ring junction. The obtained Pα product thereafter treated with dibromocarbene leads via an exothermic reaction to the six-membered ring double bond position of α-monoepoxide. The major products Pαα are kinetically and thermodynamically favored with a high stereoselectivity in perfect correlation with the experimental observations.

Conceptual DFT-Based Computational Peptidology of Marine Natural Compounds: Discodermins A-H

A methodology based on the concepts that arise from Density Functional Theory named Conceptual Density Functional Theory (CDFT) was chosen for the calculation of some global and local reactivity descriptors of the Discodermins A-H family of marine peptides through the consideration of the KID (Koopmans in DFT) technique that was successfully used in previous studies of this kind of molecular systems. The determination of active sites of the studied molecules for different kinds of reactivities was achieved by resorting to some CDFT-based descriptors like the Fukui functions as well as the Parr functions derived from Molecular Electron Density Theory (MEDT). A few properties identified with their ability to behave as a drug and the bioactivity of the peptides considered in this examination were acquired by depending on a homology model by studying the correlation with the known bioactivity of related molecules in their interaction with various biological receptors. With the further object of analyzing their bioactivity, some parameters of usefulness for future QSAR studies, their predicted biological targets, and the ADME (Absorption, Distribution, Metabolism, and Excretion) parameters related to the Discodermins A-H pharmacokinetics are also reported.

Efficient and Environmentally Friendly Adsorbent Based on β-Ketoenol-Pyrazole-Thiophene for Heavy-Metal Ion Removal from Aquatic Medium: A Combined Experimental and Theoretical Study

A new sustainable and environmentally friendly adsorbent based on a β-ketoenol-pyrazole-thiophene receptor grafted onto a silica surface was developed and applied to the removal of heavy-metal ions (Pb(II), Cu(II), Zn(II), and Cd(II)) from aquatic medium. The new material SiNPz-Th was well characterized and confirms the success of covalent binding of the receptor on the silica surface. The effect of environmental parameters on adsorption including pH, contact time, temperature, and the initial concentration were investigated. The maximum adsorption capacities of SiNPz-Th for Pb(II), Cu(II), Zn(II), and Cd(II) ions were 102.20, 76.42, 68.95, and 32.68 mg/g, respectively, at 30 min and pH = 6. The adsorption isotherms, kinetics, and thermodynamic process were investigated and showed efficiency and selectivity toward Pb(II) and good regeneration performance. Density functional theory, noncovalent-interaction, and quantum theory of atoms in molecules calculations were used to study and to gain a deeper understanding of both the adsorption mechanism and selectivity of metal ions onto the adsorbent. Accordingly, metal ions such as Pb(II), Cu(II), and Zn(II) were bidentate coordinated with the adsorbent by nitrogen and oxygen atoms of the Schiff base C=N and hydroxyl group -OH, respectively, to form stable complexes. Whereas Cd(II) was coordinated in a monodentate fashion with oxygen atom of the hydroxyl group. Furthermore, the affinity of SiNPz-Th toward the metal ions was decreased in the order of Pb(II) > Cu(II) > Zn(II) > Cd(II), in good agreement with the experimental results. All these results highlight that SiNPz-Th has good potential to be an advanced adsorbent for the removal of lead ions from real water.

Predicting the chemical reactivity of organic materials using a machine-learning approach

Stability and compatibility between chemical components are essential parameters that need to be considered in the selection of functional materials in configuring a system. In configuring devices such as batteries or solar cells, not only the functionality of individual constituting materials such as electrodes or electrolyte but also an appropriate combination of materials which do not undergo unwanted side reactions is critical in ensuring their reliable performance in long-term operation. While the universal theory that can predict the general chemical reactivity between materials is long awaited and has been the subject of studies with a rich history, traditional ways proposed to date have been mostly based on simple electronic properties of materials such as electronegativity, ionization energy, electron affinity and hardness/softness, and could be applied to only a small group of materials. Moreover, prediction has often been far from accurate and has failed to offer general implications; thus it was practically inadequate as a selection criterion from a large material database, i.e. data-driven material discovery. Herein, we propose a new model for predicting the general reactivity and chemical compatibility among a large number of organic materials, realized by a machine-learning approach. As a showcase, we demonstrate that our new implemented model successfully reproduces previous experimental results reported on side-reactions occurring in lithium–oxygen electrochemical cells. Furthermore, the mapping of chemical stability among more than 90 available electrolyte solvents and the representative redox mediators is realized by this approach, presenting an important guideline in the development of stable electrolyte/redox mediator couples for lithium–oxygen batteries.

Stability and compatibility between chemical components are essential parameters that need to be considered in the selection of functional materials in configuring a system.

1,2,4-Oxadiazole-5-ones as analogues of tamoxifen: synthesis and biological evaluation

A series of 2,3,4-triaryl-substituted 1,2,4-oxadiazole-5-ones have been prepared as fixed-ring analogues of tamoxifen (TAM), a drug inhibitor of Estradiol Receptor (ER) used in breast cancer therapy, by an efficient synthetic protocol based on a 1,3-dipolar cycloaddition of nitrones to isocyanates. Some of the newly synthesized compounds (14d-f, 14h and 14k) show a significant cytotoxic effect in a human breast cancer cell line (MCF-7) possessing IC50 values between 15.63 and 31.82 μM. In addition, compounds 14d-f, 14h and 14k are able to increase the p53 expression levels, activating also the apoptotic pathway. Molecular modeling studies of novel compounds performed on the crystal structure of ER reveal the presence of strong hydrophobic interactions with the aromatic rings of the ligands similar to TAM. These data suggest that 1,2,4-oxadiazole-5-ones can be considered analogues of TAM, and that their anticancer activity might be partially due to ER inhibition.

Calculation of the Global and Local Conceptual DFT Indices for the Prediction of the Chemical Reactivity Properties of Papuamides A-F Marine Drugs

A well-behaved model chemistry previously validated for the study of the chemical reactivity of peptides was considered for the calculation of the molecular properties and structures of the Papuamide family of marine peptides. A methodology based on Conceptual Density Functional Theory (CDFT) was chosen for the determination of the reactivity descriptors. The molecular active sites were associated with the active regions of the molecules related to the nucleophilic and electrophilic Parr functions. Finally, the drug-likenesses and the bioactivity scores for the Papuamide peptides were predicted through a homology methodology relating them with the calculated reactivity descriptors, while other properties such as the pKas were determined following a methodology developed by our group.

Chemical reactivity and bioactivity properties of the Phallotoxin family of fungal peptides based on Conceptual Peptidology and DFT study

A methodology based on the concepts that arise from Density Functional Theory (CDFT) was chosen for the calculation of the global and local reactivity descriptors of the Phallotoxin family of fungal peptides. The determination of the active sites for the molecules has been achieved by resorting some descriptors within Molecular Electron Density Theory (MEDT) like the Dual Descriptor and the Parr functions. Phallosacin has been found as the most reactive of the peptides on the basis of the calculated Global Reactivity Descriptors. The pKas of the seven studied peptides were established using a proposed relationship between this property and the calculated Global Hardness. The bioactivity properties of the peptides considered in this study were obtained by resorting to a homology model by comparison with the bioactivity of related molecules in their interaction with different receptors.

A molecular electron density theory study on the [3+2] cycloaddition reaction of 5,5-dimethyl-1-pyrroline N-oxide with 2-cyclopentenone

In the present work, the [3 + 2] cycloaddition reaction of 5,5-dimethyl-1-pyrroline N-oxide (Nit-5) and 2-cyclopentenone (CPN-6), experimentally reported by Tamura et al., was theoretically studied using the newly introduced molecular electron density theory (MEDT). Based on the experimental findings, this reaction takes place in an O3-C4 regio- and an exo-stereospecific fashion to give corresponding [3 + 2] exo cycloadduct as the sole product. The results of the potential energy surface analysis indicated that the experimentally reported product is more favorable both thermodynamically and kinetically relative to other possible adducts. In complete agreement with the experimental outcomes, the conceptual density functional theory reactivity indices explained the reactivity and regioselectivity of the reaction. Calculation of global electron density transfer of the energetically most preferred transition state indicated that the electron density fluxes from Nit-5 as a nucleophilic species toward CPN-6 as an electrophilic species. Analysis of the molecular electrostatic potential map of the most favorable transition state showed that approach of Nit-5 and CPN-6 locates the oppositely charged regions over each other leading to attractive forces between two reagents rationalizing the exo stereoselectivity predominance. The molecular mechanism of the reactions was specified using electron localization function analysis over some relevant points along the intrinsic reaction coordinate profile of the most favorable transition state and the results indicated that this zwitterionic-type [3 + 2] cycloaddition reaction proceeds through a two-stage one-step mechanism. In fact, while the O3-C4 single bond is initialy formed between two fragments through donation of some electron density from the O3 oxygen lone electron-pairs of Nit-5 toward the C4 carbon atom of CPN-6, the delayed C1-C5 single bond begins to form via C1- to -C5 coupling of pseudodiracal centers created on theses atoms over the course of reaction.

Computational Peptidology Assisted by Conceptual Density Functional Theory for the Study of Five New Antifungal Tripeptides

A well-behaved model chemistry previously validated for the study of the chemical reactivity of peptides was considered for the calculation of the molecular properties and structures of a group of five new antifungal tripeptides, namely (2R)-2-[(2S)-2-[(2S)-2-amino-3-phenylpropanamido]propanamido]-5-[(diaminomethylidene)amino]pentanoic acid, (2S)-2-[(2S)-2-[(2S)-2-amino-3-phenyl propanamido]propanamido]-3-(4-hydroxyphenyl)propanoic acid, (2S)-2-[(2S)-2-[(2S)-2-amino-3-phenylpropanamido]-3-methylbutanamido]-3-(4-hydroxyphenyl)propanoic acid, (2R)-2-[(2S)-2-[(2S)-2-amino-3-phenylpropanamido]-3-(1H-indol-3-yl)propanamido]-3-sulfanylpropanoic acid, and (2S)-2-[(2S)-2-[(2S)-2-amino-3-phenylpropanamido]-3-(1H-indol-3-yl)propanamido]-3-(4-hydroxyphenyl)propanoic acid, according to their amino acid sequences. A methodology based on conceptual density functional theory was chosen for the determination of the reactivity descriptors. The molecular active sites were associated with the active regions of the molecules that were associated with the nucleophilic, electrophilic, and radical Fukui functions. Additionally, the pK a values for the different peptides are predicted with great accuracy, which constitutes a useful knowledge for the process of drug design. Finally, the bioactivity scores for the new antifungal peptides are predicted through a homology methodology relating them with the calculated reactivity descriptors.

CDFT-Based Reactivity Descriptors as a Useful MEDT Chemoinformatics Tool for the Study of the Virotoxin Family of Fungal Peptides

Virotoxins are monocyclic peptides formed by at least five different compounds: alaviroidin, viroisin, deoxoviroisin, viroidin and deoxovirodin. These are toxic peptides singularly found in Amanita virosa mushrooms. Here we perform computational studies on the structural and electronic conformations of these peptides using the MN12SX/Def2TZVP/H2O chemistry model to investigate their chemical reactivity. CDFT-based descriptors (for Conceptual Density Functional Theory) (e.g., Parr functions and Nucleophilicity) are also considered. At the same time, other properties (e.g., pKas) will be determined and used to study virotoxins solubility and to inform decisions about repurposing these agents in medicinal chemistry.

Quantum Reality in the Selective Reduction of a Benzofuran System

Two 2,3-disubstituted benzofurans (1 and 2), analogs of gamma-aminobutyric acid (GABA), were synthesized to obtain their 2,3-dihydro derivatives from the Pd/C-driven catalytic reduction of the double bond in the furanoid ring. The synthesis produced surprising by-products. Therefore, theoretical calculations of global and local reactivity were performed based on Pearson's hard and soft acids and bases (HSAB) principle to understand the regioselectivity that occurred in the reduction of the olefinic carbons of the compounds. Local electrophilicity (ω_k) was the most useful parameter for explaining the selectivity of the polar reactions. This local parameter was defined with the condensed Fukui function and redefined with the electrophilic (P_k⁺) Parr function. The similar patterns of both resulting sets of values helped to demonstrate the electrophilic behavior (soft acid) of the olefinic carbons in these compounds. The theoretical calculations, nuclear magnetic resonance, and resonance hybrids showed the moieties in each compound that are most susceptible to reduction.