Atomic Decomposition Scheme of Noncovalent Interactions Applied to Host–Guest Assemblies

Theoretical Study of Cyanidin-Resveratrol Copigmentation by the Functional Density Theory

Anthocyanins are colored water-soluble plant pigments. Upon consumption, anthocyanins are quickly absorbed and can penetrate the blood-brain barrier (BBB). Research based on population studies suggests that including anthocyanin-rich sources in the diet lowers the risk of neurodegenerative diseases. The copigmentation caused by copigments is considered an effective way to stabilize anthocyanins against adverse environmental conditions. This is attributed to the covalent and noncovalent interactions between colored forms of anthocyanins (flavylium ions and quinoidal bases) and colorless or pale-yellow organic molecules (copigments). The present work carried out a theoretical study of the copigmentation process between cyanidin and resveratrol (CINRES). We used three levels of density functional theory: M06-2x/6-31g+(d,p) (d3bj); ωB97X-D/6-31+(d,p); APFD/6-31+(d,p), implemented in the Gaussian16W package. In a vacuum, the CINRES was found at a copigmentation distance of 3.54 Å between cyanidin and resveratrol. In water, a binding free energy ∆G was calculated, rendering -3.31, -1.68, and -6.91 kcal/mol, at M06-2x/6-31g+(d,p) (d3bj), ωB97X-D/6-31+(d,p), and APFD/6-31+(d,p) levels of theory, respectively. A time-dependent density functional theory (TD-DFT) was used to calculate the UV spectra of the complexes and then compared to its parent molecules, resulting in a lower energy gap at forming complexes. Excited states' properties were analyzed with the ωB97X-D functional. Finally, Shannon aromaticity indices were calculated and isosurfaces of non-covalent interactions were evaluated.

Volatile lanthanide complexes with fluorinated heptadentate ligands

Understanding factors that influence the volatility of lanthanide complexes remains an important goal for applications such as gas-phase f-metal separations and the synthesis of lanthanide-containing thin films. Lanthanide complexes often exhibit volatility differences that depend on the ability of ligands to saturate the lanthanide coordination sphere and attenuate intermolecular bonding in the solid state. This can make it difficult to assess how electronic factors associated with differing ligand substituents influence volatility. Here we describe the synthesis, structures, and thermal properties of a series of volatile lanthanide complexes (Ln = Nd, Er, and Yb) containing N4O3 ligands decorated with different alkyl and fluoroalkyl substituents (CF3, CF2CF2CF3, Me, and tBu). These ligands completely enveloped the tested lanthanides to form monomeric complexes with 7-coordinate distorted capped octahedral coordination geometries, as determined using single-crystal X-ray diffraction. Thermogravimetric analysis and bulk sublimation studies show how metal encapsulation affords complexes with the same volatility regardless of metal size, even with lanthanide ions with significantly different radii such as Nd3+ and Yb3+. Most notably, the results show that increasing ligand fluorination, a strategy often used to increase the volatility of metal complexes, is not always beneficial and can significantly attenuate the volatility of lanthanide complexes depending on location with respect to other substituents in the ligand framework. A pair-wise model based on density functional theory shows that the net intermolecular interactions in the unit cell can still be stronger when fluorination is present. In other words, even if individual interactions between neighboring molecules are weaker, the total number of interactions in the solid arising from the nature of crystal packing is equally important to consider.

A novel alpha-amylase inhibitor-based spirooxindole-pyrrolidine-clubbed thiochromene-pyrzaole pharmacophores: Unveiling the [3+2] cycloaddition reaction by molecular electron density theory

A novel spirooxindole-pyrrolidine clubbed thiochromene and pyrazole motifs were synthesized by [3+2] cycloaddition (32CA) reactions in one step process starting from the ethylene-based thiochromene and pyrazole scaffolds with the secondary amino-acids and substituted isatins in high yield. The 32CA reaction of AY 10 with ethylene derivative 6 has also been studied with Molecular Electron Density Theory. The high nucleophilic character of AY 10, N = 4.39 eV, allows explaining that the most favorable TS-on is 13.9 kcal mol-1 below the separated reagent. This 32CA, which takes place through a non-concerted one-step mechanism, presents a total ortho regio- and endo stereoselectivity, which is controlled by the formation of two intramolecular H… O hydrogen bonds. The design of spirooxindole-pyrrolidines engrafted thiochromene and pyrazole was tested for alpha-amylase inhibition and show a high efficacy in nanoscale range of reactivity. The key interaction between the most active hybrids and the receptor was studied by molecular docking. The physiochemical properties of the designed spirooxindole-pyrrolidines were carried out by in silico ADMET prediction. The newly synthesized most potent hybrid could be considered as a lead compound for drug discovery development for type 2 diabetes mellitus (T2DM).

On the Importance of Noncovalent Interactions in the Stabilization of Nonconventional Compounds Using Bulky Groups

In this article, we studied the capability of bulky groups to contribute to the stabilization of a given compound in addition to the well-known steric effect related to substituents due to their composition (alkyl chains and aromatic groups, among others). For this purpose, the recently synthesized 1-bora-3-boratabenzene anion which contains large substituents was analyzed by means of the independent gradient model (IGM), natural population analysis (NPA) at the TPSS/def2-TZVP level, force field-based energy decomposition analysis (EDA-FF) applying the universal force field (UFF), and molecular dynamics calculations under the GFN2-xTB approach. The results indicate that the bulky groups should not only be considered for their steric effects but also for their ability to stabilize a system that could be very reactive.

IGMPlot: A program to identify, characterize, and quantify molecular interactions

We describe the development and features of a program called IGMPlot, which is based on the independent gradient model (IGM) and its local descriptor δ g $$ \delta g $$ . The IGM approach analyzes the gradient of the electron density (ED) in a molecular system to identify regions of space where chemical interactions take place. IGMPlot is intended for use by both experimental scientists and theoretical chemists. It is standalone software written in C++, with versions available for multiple platforms. Some key features are: probing and quantifying interactions between two given molecular fragments, determining bond strength (IBSI), estimating the atomic contributions to an intermolecular interaction and preparing data to build 2D and 3D representations of interaction regions. The software has been updated to include new features: critical point analysis of the ED, assessing ED asymmetry of a given bond (PDA) and a new descriptor called q g $$ qg $$ designed to enhance the IGM- δ g $$ \delta g $$ analysis. The program can be found at: http://igmplot.univ-reims.fr.

A New Motif in Halogen Bonding: Cooperative Intermolecular S-Br⋅⋅⋅O, O⋅⋅⋅F, and F⋅⋅⋅F Associations in the Crystal Packing of α,ω-Di(sulfonyl bromide) Perfluoroalkanes

We recently reported the first examples of S-Cl⋅⋅⋅O halogen bonding complemented by short F⋅⋅⋅F contacts between neighboring chains that resulted in stabilized crystals of ClSO₂ (CF₂ )₄ SO₂ Cl and ClSO₂ (CF₂ )₆ SO₂ Cl. More recently, other researchers studied our crystallographic data further using an Independent Gradient Model (IGM), and they suggested if one goes beyond IUPAC's proposed 'less than the sum of the van der Waals radii' criterion that even more noncovalent interactions between fluorine atoms on neighboring chains as well as Cl⋅⋅⋅Cl, Cl⋅⋅⋅S, O⋅⋅⋅F, and O⋅⋅⋅S attractive interactions can be found. With that said, we have prepared samples of the related BrSO₂ (CF₂ )_n SO₂ Br derivatives (where n=4, 6, 8, and others), which give rise to even stronger S-Br⋅⋅⋅O halogen bonding interactions complemented minimally by O⋅⋅⋅F and F⋅⋅⋅F intermolecular interactions as shown by X-ray crystallography and computational chemistry using IGM isosurface plots. Additional spectroscopic characterization (multinuclear NMR, FT-IR, and MS) of the disulfonyl bromide derivatives BrSO₂ (CF₂ )₄ SO₂ Br, BrSO₂ (CF₂ )₆ SO₂ Br, and BrSO₂ (CF₂ )₈ SO₂ Br has also been obtained as well as some preliminary spectroscopic evidence for BrSO₂ (CF₂ )₂ SO₂ Br and BrSO₂ CF₂ O(CF₂ )₂ OCF₂ SO₂ Br. The implication of these results toward the preparation of the corresponding disulfonyl iodides is discussed.

Extracting Quantitative Information at Quantum Mechanical Level from Noncovalent Interaction Index Analyses

The noncovalent interaction (NCI) index is nowadays a well-known strategy to detect NCIs in molecular systems. Even though it initially provided only qualitative descriptions, the technique has been recently extended to also extract quantitative information. To accomplish this task, integrals of powers of the electron distribution were considered, with the requirement that the overall electron density can be clearly decomposed as sum of distinct fragment contributions to enable the definition of the (noncovalent) integration region. So far, this was done by only exploiting approximate promolecular electron densities, which are given by the sum of spherically averaged atomic electron distributions and thus represent too crude approximations. Therefore, to obtain more quantum mechanically (QM) rigorous results from NCI index analyses, in this work, we propose to use electron densities obtained through the transfer of extremely localized molecular orbitals (ELMOs) or through the recently developed QM/ELMO embedding technique. Although still approximate, the electron distributions resulting from the abovementioned methods are fully QM and, above all, are again partitionable into subunit contributions, which makes them completely suitable for the NCI integral approach. Therefore, we benchmarked the integrals resulting from NCI index analyses (both those based on the promolecular densities and those based on ELMO electron distributions) against interaction energies computed at a high quantum chemical level (in particular, at the coupled cluster level). The performed test calculations have indicated that the NCI integrals based on ELMO electron densities outperform the promolecular ones. Furthermore, it was observed that the novel quantitative NCI-(QM/)ELMO approach can be also profitably exploited both to characterize and evaluate the strength of specific interactions between ligand subunits and protein residues in protein-ligand complexes and to follow the evolution of NCIs along trajectories of molecular dynamics simulations. Although further methodological improvements are still possible, the new quantitative ELMO-based technique could be already exploited in situations in which fast and reliable assessments of NCIs are crucial, such as in computational high-throughput screenings for drug discovery.

Revisiting the Sweet Taste Receptor T1R2-T1R3 through Molecular Dynamics Simulations Coupled with a Noncovalent Interactions Analysis

It is nowadays widely accepted that sweet taste perception is elicited by the activation of the heterodimeric complex T1R2-T1R3, customarily known as sweet taste receptor (STR). However, the interplay between STR and sweeteners has not yet been fully clarified. Here through a methodology coupling molecular dynamics and the independent gradient model (igm) approach we determine the main interacting signatures of the closed (active) conformation of the T1R2 Venus flytrap domain (VFD) toward aspartame. The igm methodology provides a rapid and reliable quantification of noncovalent interactions through a score (Δg^inter) based on the attenuation of the electronic density gradient when two molecular fragments approach each other. Herein, this approach is coupled to a 100 ns molecular dynamics simulation (MD-igm) to explore the ligand-cavity contacts on a per-residue basis as well as a series of key inter-residue interactions that stabilize the closed form of VFD. We also apply an atomic decomposition scheme of noncovalent interactions to quantify the contribution of the ligand segments to the noncovalent interplay. Finally, a series of structural modification on aspartame are conducted in order to obtain guidelines for the rational design of novel sweeteners. Given that innovative methodologies to reliably quantify the extent of ligand-protein coupling are strongly demanded, this approach combining a noncovalent analysis and MD simulations represents a valuable contribution, that can be easily applied to other relevant biomolecular systems.

A Theoretical Study of the Halogen Bond between Heteronuclear Halogen and Benzene

Halogen bonds play an important role in many fields, such as biological systems, drug design and crystal engineering. In this work, the structural characteristics of the halogen bond between heteronuclear halogen XD (ClF, BrCl, IBr, ICl, BrF and IF) and benzene were studied using density functional theory. The structures of the complexes between heteronuclear halogen and benzene have Cs symmetry. The interaction energies of the complexes between heteronuclear halogen XD (ClF, BrCl, IBr, ICl, BrF and IF) and benzene range from -27.80 to -37.18 kJ/mol, increasing with the increases in the polarity between the atoms of X and D, and are proportional to the angles of a between the Z axis and the covalent bond of heteronuclear halogen. The electron density (ρ) and corresponding Laplacian (∇²ρ) values indicate that the interaction of the heteronuclear halogen and benzene is a typical long-range weak interaction similar to a hydrogen bond. Independent gradient model analysis suggests that the van der Waals is the main interaction between the complexes of heteronuclear halogen and benzene. Symmetry-adapted perturbation theory analysis suggests that the electrostatic interaction is the dominant part in the complexes of C₆H₆⋯ClF, C₆H₆⋯ICl, C₆H₆⋯BrF and C₆H₆⋯IF, and the dispersion interaction is the main part in the complexes of C₆H₆⋯BrCl, C₆H₆⋯IBr.

A Molecular Electron Density Theory Study of the [3+2] Cycloaddition Reaction of an Azomethine Ylide with an Electrophilic Ethylene Linked to Triazole and Ferrocene Units

The [3+2] cycloaddition (32CA) reaction of an azomethine ylide (AY) with an electrophilic ethylene linked to triazole and ferrocene units has been studied within the Molecular Electron Density Theory (MEDT) at the ωB97X-D/6-311G(d,p) level. The topology of the electron localization function (ELF) of this AY allows classifying it as a pseudo(mono)radical species characterized by the presence of two monosynaptic basins, integrating a total of 0.76 e, at the C1 carbon. While the ferrocene ethylene has a strong electrophilic character, the AY is a supernucleophile, suggesting that the corresponding 32CA reaction has a high polar character and a low activation energy. The most favorable ortho/endo reaction path presents an activation enthalpy of 8.7 kcal·mol^-1, with the 32CA reaction being exergonic by -42.1 kcal·mol^-1. This reaction presents a total endo stereoselectivity and a total ortho regioselectivity. Analysis of the global electron density transfer (GEDT) at the most favorable TS-on (0.23 e) accounts for the high polar character of this 32CA reaction, classified as forward electron density flux (FEDF). The formation of two intermolecular hydrogen bonds between the two interacting frameworks at the most favorable TS-on accounts for the unexpected ortho regioselectivity experimentally observed.

Meta-Hybrid Density Functional Theory Prediction of the Reactivity, Stability, and IGM of Azepane, Oxepane, Thiepane, and Halogenated Cycloheptane

The application of plain cycloalkanes and heterocyclic derivatives in the synthesis of valuable natural products and pharmacologically active intermediates has increased tremendously in recent times with much attention being paid to the lower cycloalkane members. The structural and molecular properties of higher seven-membered and nonaromatic heterocyclic derivatives are less known despite their stable nature and vast application; thus, an insight into their structural and electronic properties is still needed. Appropriate quantum chemical calculations utilizing the ab initio (MP2) method, meta-hybrid (M06-2X) functional, and long-range-separated functionals (ωB97XD) have been utilized in this work to investigate the structural reactivity, stability, and behavior of substituents on cycloheptane (CHP) and its derivatives: azepane, oxepane, thiepane, fluorocycloheptane (FCHP), bromocycloheptane (BrCHP), and chlorocycloheptane (ClCHP). Molecular global reactivity descriptors such as Fukui function, frontier molecular orbitals (FMOs), and molecular electrostatic potential were computed and compared with lower members. The results of two population methods CHELPG and Atomic Dipole Corrected Hirshfeld Charges (ADCH) were equally compared to scrutinize the charge distribution in the molecules. The susceptibility of intramolecular interactions between the substituents and cycloalkane ring is revealed by natural bond orbital analysis and intramolecular weak interactions by the independent gradient model (IGM). Other properties such as atomic density of states, intrinsic bond strength index (IBSI), and dipole moments are considered. It is acclaimed that the strain effect is a major determinant effect in the energy balance of cyclic molecules; thus, the ring strain energies and validation of spectroscopic specificities with reference to the X-ray crystallographic data are also considered.

Structural study of 1- and 2-naphthol: new insights into the non-covalent H-H interaction in cis-1-naphthol

Previous microwave studies of naphthol monomers were supplemented by measuring spectra of all ¹³C mono-substituted isotopologues of the cis- and trans-conformers of 1-naphthol and 2-naphthol in their natural abundances. The resulting data were utilized to determine substitution structures and so-called semi-experimental effective structures. Results from electronic structure calculations show that the OH group of cis-1-naphthol points ≈6° out of plane, which is consistent with the inertial defect data of cis- and trans-1-naphthol. The non-planarity of cis-1-naphthol is a result of a close-contact H-atom-H-atom interaction. This type of H-H interaction has been the subject of much controversy in the past and we provide here an in-depth theoretical analysis of it. The naphthol system is particularly well-suited for such analysis as it provides internal standards with its four different isomers. The methods used include quantum theory of atoms in molecules, non-covalent interactions, independent gradient model, local vibrational mode, charge model 5, and natural bond orbital analyses. We demonstrate that the close-contact H-H interaction is neither a purely attractive nor repulsive interaction, but rather a mixture of the two.

Independent gradient model based on Hirshfeld partition: A new method for visual study of interactions in chemical systems

The powerful independent gradient model (IGM) method has been increasingly popular in visual analysis of intramolecular and intermolecular interactions in recent years. However, we frequently observed that there is an evident shortcoming of IGM map in graphically studying weak interactions, that is its isosurfaces are usually too bulgy; in these cases, not only the graphical effect is poor, but also the color on some areas on the isosurfaces is inappropriate and may lead to erroneous analysis conclusions. In addition, the IGM method was originally proposed based on promolecular density, which is quite crude and does not take actual electronic structure into account. In this article, we propose an improvement version of IGM, namely IGM based on Hirshfeld partition of molecular density (IGMH), which replaces the free-state atomic densities involved in the IGM method with the atomic densities derived by Hirshfeld partition of actual molecular electron density. This change makes IGM have more rigorous physical background. A large number of application examples in this article, including molecular and periodic systems, weak and chemical bond interactions, fully demonstrate the important value of IGMH in intuitively understanding interactions in chemical systems. Comparisons also showed that the IGMH usually has markedly better graphical effect than IGM and overcomes known problems in IGM. Currently IGMH analysis has been supported in our wavefunction analysis code Multiwfn (http://sobereva.com/multiwfn). We hope that IGMH will become a new useful method among chemists for exploring interactions in wide variety of chemical systems.

Reactivity of Coinage Metal Hydrides for the Production of H2 Molecules

The formation of molecular hydrogen as well as the possibility of using coinage metal hydrides as a prospective complex to produce hydrogen was presented in this work. Therefore, the reactions involving the interaction between two coinage metal hydrides, MH (M=Cu, Ag and Au, homo and heterodimers), were studied. The free energy profiles corresponding to aforementioned complexation were analysed by means of ab initio methods of quantum chemistry. The characteristics of these intermediates, final complexes and the electron density properties of the established interactions were discussed.

Evaluation of Electron Density Shifts in Noncovalent Interactions

In the present paper, we report the quantitative evaluation of the electron density shift (EDS) maps within different complexes. Values associated with the total EDS maps exhibited good correlation with different quantities such as interaction energies, Eint, intermolecular distances, bond critical points, and LMOEDA energy decomposition terms. Besides, EDS maps at different cutoffs were also evaluated and related with the interaction energies values. Finally, EDS maps and their corresponding values are found to correlate with Eint within systems with cooperative effects. To our knowledge, this is the first time that the EDS has been quanitatively evaluated.

Pyridazinone derivatives as potential anti-inflammatory agents: synthesis and biological evaluation as PDE4 inhibitors

Cyclic nucleotide phosphodiesterase type 4 (PDE4), which controls the intracellular level of cyclic adenosine monophosphate (cAMP), has aroused scientific attention as a suitable target for anti-inflammatory therapy of respiratory diseases. This work describes the development and characterization of pyridazinone derivatives bearing an indole moiety as potential PDE4 inhibitors and their evaluation as anti-inflammatory agents. Among these derivatives, 4-(5-methoxy-1H-indol-3-yl)-6-methylpyridazin-3(2H)-one possesses promising activity, and selectivity towards PDE4B isoenzymes and is able to regulate potent pro-inflammatory cytokine and chemokine production by human primary macrophages.

Exploration of adsorption mechanism of 2-phosphonobutane-1,2,4-tricarboxylic acid onto kaolinite and montmorillonite via batch experiment and theoretical studies

Two clay minerals, kaolinite (Kaol) and montmorillonite (Mt) with different crystal structures were chosen to investigate the comparative adsorption of 2-phosphonobutane-1,2,4-tricarboxylic acid (PBTC) through batch control experiments and theoretical studies. The systematical isotherm and kinetic studies agreed with Langmuir model and pseudo-second-order model, confirming a monolayer and chemisorption interaction process, respectively. The maximum removal capacities of Kaol and Mt for PBTC were 72.297 mg/g and 121.163 mg/g at pH=3.0 and T=298 K, respectively. Furthermore, the adsorption mechanisms were investigated by molecular dynamic (MD) simulations and density functional theory (DFT). The Interface force field (IFF) was firstly introduced into Materials Studio package to explore the microscopic mechanism of clay mineral interface. The dynamics behaviors verified that the oxygen (O) atom of carboxyl group has stronger affinity at the external surface of Mt, which consistent with the experimental data well. For DFT calculations, quantitative analysis around molecular van der Waals (vdW) surface was adopted to predict reactive sites for the electrophilic reaction. Independent Gradient Model (IGM) and Hirshfeld surface analyses in Multiwfn indicated that the high adsorption effect mainly attributes to hydrogen bond action. These findings improve our ability to explore the related properties occurring at the interface of different clay minerals.

A Step toward the Quantification of Noncovalent Interactions in Large Biological Systems: The Independent Gradient Model-Extremely Localized Molecular Orbital Approach

The independent gradient model (IGM) is a recent electron density-based computational method that enables to detect and quantify covalent and noncovalent interactions. When applied to large systems, the original version of the technique still relies on promolecular electron densities given by the sum of spherically averaged atomic electron distributions, which leads to approximate evaluations of the inter- and intramolecular interactions occurring in systems of biological interest. To overcome this drawback and perform IGM analyses based on quantum mechanically rigorous electron densities also for macromolecular systems, we coupled the IGM approach with the recently constructed libraries of extremely localized molecular orbitals (ELMOs) that allow fast and reliable reconstructions of polypeptide and protein electron densities. The validation tests performed on small polypeptides and peptide dimers have shown that the novel IGM-ELMO strategy provides results that are systematically closer to the fully quantum mechanical ones and outperforms the IGM method based on the crude promolecular approximation, but still keeping a quite low computational cost. The results of the test calculations carried out on proteins have also confirmed the trends observed for the IGM analyses conducted on small systems. This makes us envisage the future application of the novel IGM-ELMO approach to unravel complicated noncovalent interaction networks (e.g., in protein-protein contacts) or to rationally design new drugs through molecular docking calculations and virtual high-throughput screenings.

The role of zero-field splitting and π-stacking interaction of different nitrogen-donor ligands on the optical properties of luminescent rhenium tricarbonyl complexes

In this work, a systematic evaluation of the role of zero-field splitting (ZFS), and the geometric arrangement of different nitrogen-donor ligands, including π-stacking interactions, in five selected rhenium luminescent complexes was performed.

Exploring Epigenetic Marks by Analysis of Noncovalent Interactions

Epigenetic marks are modest chemical modifications on DNA and histone proteins that regulate the activation or silencing of genes through modulation of the intermolecular interactions between the DNA strands and the protein machinery. The process is complex and not always well understood. One of the systems studied in greater detail is the epigenetic mark on H3K9: lysine 9 of histone 3. The degree of methylation or acetylation of this histone is linked to silencing or activation of the corresponding gene, but it is not clear which effect each mark has in gene expression. We shed light on this particular methylation process by using density functional theory (DFT) calculations. First, we built a model consisting of a DNA double strand containing three base pairs and a sequence of three amino acids of the histone's tail. Then, we computed the modulation introduced into the intermolecular interactions by each epigenetic modification: from mono- to trimethylation and acetylation. The calculations show that whereas acetylation and trimethylation result in a reduction of the DNA-peptide interaction; non-, mono-, and dimethylation increase the intermolecular interactions. Such observations compare well with the findings reported in the literature, and highlight the correlation between the balance of intermolecular forces and biological properties, simultaneously advancing quantum-mechanical studies of large biochemical systems at molecular level through the use of DFT methods.

Understand the Specific Regio- and Enantioselectivity of Fluostatin Conjugation in the Post-Biosynthesis

Fluostatins, benzofluorene-containing aromatic polyketides in the atypical angucycline family, conjugate into dimeric and even trimeric compounds in the post-biosynthesis. The formation of the C-C bond involves a non-enzymatic stereospecific coupling reaction. In this work, the unusual regio- and enantioselectivities were rationalized by density functional theory calculations with the M06-2X (SMD, water)/6-311 + G(d,p)//6-31G(d) method. These DFT calculations reproduce the lowest energy C1-(R)-C10'-(S) coupling pathway observed in a nonenzymatic reaction. Bonding of the reactive carbon atoms (C1 and C10') of the two reactant molecules maximizes the HOMO-LUMO interactions and Fukui function involving the highest occupied molecular orbital (HOMO) of nucleophile p-QM and lowest unoccupied molecular orbital (LUMO) of electrophile FST2- anion. In particular, the significant π-π stacking interactions of the low-energy pre-reaction state are retained in the lowest energy pathway for C-C coupling. The distortion/interaction-activation strain analysis indicates that the transition state (TScp-I) of the lowest energy pathway involves the highest stabilizing interactions and small distortion among all possible C-C coupling reactions. One of the two chiral centers generated in this step is lost upon aromatization of the phenol ring in the final difluostatin products. Thus, the π-π stacking interactions between the fluostatin 6-5-6 aromatic ring system play a critical role in the stereoselectivity of the nonenzymatic fluostatin conjugation.

New Way for Probing Bond Strength

The covalent chemical bond is intimately linked to electron sharing between atoms. The recent independent gradient model (IGM) and its δg descriptor provide a way to quantify locally this electron density interpenetration from wavefunction calculations. Each bond has its own IGM-δg^pair signature. The present work establishes for the first time a strong link between this bond signature and the physically grounded bond force constant concept. Analyzing a large set of compounds and bonds, the intrinsic bond strength index (IBSI) emerges from the IGM formulation. Our study shows that the IBSI does not belong to the class of conventional bond orders (like Mulliken, Wiberg, Mayer, delocalization index, or electron localization function-ELF), but is rather a new complementary index, related to the bond strength. A fundamental outcome of this research is a novel index allowing to range all two-center chemical bonds by their intrinsic strength in molecular situation. We believe that the IBSI is a powerful and robust tool for interpretation accessible to a wide community of chemists (organic, inorganic chemistry, including transition-metal complexes and reaction mechanisms).

The usefulness of energy decomposition schemes to rationalize host–guest interactions

The findings reported here reveal the robustness and practical application of EDA-NOCV in rationalizing molecular recognition situations in host–guest systems.