Search

Bridgehead alkenes are polycyclic molecules bearing at least one C=C bond that includes a bridgehead carbon atom. For small bicyclic systems, these bonds are highly strained due to geometric constraints placed on the sp2 hybridized carbon atoms. These small, strained molecules have been termed “anti-Bredt” alkenes. β-halo carbanions have served as convenient precursors to bridgehead alkenes in experimental studies. We observed that upon attempted computational geometric optimizations (ωB97X-D/aug-cc-pVDZ) of the precursors, spontaneous elimination of the halide occurs along with formation of the anti-Bredt alkene in many cases. Such computational eliminations were shown to faithfully mimic experimentally obtained results. Computational elimination was not observed for [1.1.1] or [2.1.1] frameworks, in agreement with predictions that these bridgehead alkenes are too strained to be formed. However, computational elimination from the [2.2.1] framework was observed to form 1-norbornene, a compound suggested in experimental work to be a reactive intermediate. Similarly, [3.1.1] frameworks and higher led to eliminations upon computational geometric optimization, in agreement with experimental findings. Natural bond order (NBO) calculations of the starting geometries proved to be excellent predictors as to whether elimination would take place. Those precursor compounds exhibiting delocalization energies in the order of 10 kcal/mol between the lone-pair electrons of the carbon atom and σ*C-Br were generally found to undergo elimination. Thus, computational optimization of β-halo substituted bicyclic precursor anions can be used to predict whether strained anti-Bredt alkenes are likely to be formed, thereby saving valuable time and costs in the experimental lab.

In this research work, we examined the decomposition mechanisms of N-substituted diacetamides. We focused on the substituent effect on the nitrogen lone-pair electron delocalization, with electron-withdrawing and electron donor groups. DFT functionals used the following: B1LYP, B3PW91, CAMB3LYP, LC-BLYP, and X3LYP. Dispersion corrections (d3bj) with Becke–Johnson damping were applied when necessary to improve non-covalent interactions in the transition state. Pople basis sets with higher angular moments and def2-TZVP basis sets were also applied and were crucial for obtaining consistent thermodynamic parameters. The proposed mechanism involves a six-membered transition state with the extraction of an α hydrogen. Several conformers of N-diacetamides were used to account for the decrease in entropy in the transition state in the rate-determining state. All calculations, including natural bond orbital (NBO) analyses, were performed using the Gaussian16 computational package and its GaussView 6.0 visualizer, along with VMD and GNUPLOT software. The isosurfaces and IBSIs were calculated using MultiWFN and IGMPlot, respectively.

The [ZnL(ONO2)2] 1 and [ZnL(NCS)2] 2 complexes were synthesized using self-assembly of the s-triazine tridentate ligand (L) with Zn(NO3)2·6H2O and Zn(ClO4)2·6H2O/NH4SCN, respectively. The Zn(II) is further coordinated by two nitrate and two isothiocyanate groups as monodentate ligands in 1 and 2, respectively. Both complexes have distorted square pyramidal coordination environments where the extent of distortion is found to be greater in 2 (τ5 = 0.41) than in 1 (τ5 = 0.28). Hirshfeld calculations explored the significant C···O, C···C, N···H, and O···H contacts in the molecular packing of both complexes. The energy framework analysis gave the total interaction energies of −317.8 and −353.5 kJ/mol for a single molecule in a 3.8 Å cluster of 1 and 2, respectively. The total energy diagrams exhibited a strong resemblance to the dispersion energy frameworks in both complexes. NBO charge analysis predicted the charges of the Zn(II) in complexes 1 and 2 to be 1.217 and 1.145 e, respectively. The electronic configuration of Zn1 is predicted to be [core] 4S0.32 3d9.98 4p0.45 4d0.02 5p0.01 for 1 and [core] 4S0.34 3d9.97 4p0.53 4d0.02 for 2. The increased occupancy of the valence orbitals is attributed to the donor→acceptor interactions from the ligand groups to Zn(II). The Zn(II) complexes were examined for their cytotoxic and antimicrobial activities. Both 1 and 2 have good cytotoxic efficiency towards HCT-116 and A-549 cancerous cell lines. We found that 1 is more active (IC50 = 29.53 ± 1.24 and 35.55 ± 1.69 µg/mL) than 2 (IC50 = 41.25 ± 2.91 and 55.05 ± 2.87 µg/mL) against both cell lines. Also, the selectivity indices for the Zn(II) complexes are higher than one, indicating their suitability for use as anticancer agents. In addition, both complexes have broad-spectrum antimicrobial activity (IC50 = 78–625 μg/mL) where the best result is found for 2 against P. vulgaris (IC50 = 78 μg/mL). Its antibacterial activity is found to be good compared to gentamycin (5 μg/mL) as a positive control against this microbe.

The Osme bond is defined as pairing a Group 8 metal atom as an electron acceptor in a noncovalent interaction with a nucleophile. DFT calculations with the ωB97XD functional consider MO4 (M = Ru, Os) as the Lewis acid, paired with a series of π electron donors C2H2, C2H4, C6H6, C4H5N, C4H4O, and C4H4S. The calculations establish interaction energies in the range between 9.5 and 26.4 kJ/mol. Os engages in stronger interactions than does Ru, and those involving more extensive π-systems within the aromatic rings form stronger bonds than do the smaller ethylene and acetylene. Extensive analysis questions the existence of a true Osme bond, as the bonding chiefly involves interactions with the three O atoms of MO4 that lie closest to the π-system, via π(C-C)→σ*(M-O) transfers. These interactions are supplemented by back donation from M-O bonds to the π*(CC) antibonding orbitals of the π-systems. Dispersion makes a large contribution to these interactions, higher than electrostatics and much greater than induction.

The present study elaborates on the synthesis, crystal structure, and computational studies of two new ionic liquids. In the crystal structure, [C5H12N][C21H14ClN2O2S] (4a), the anions form chains along the a-axis direction through C—H···π(ring) interactions. These are connected into layers that run approximately parallel to the ac plane by a variety of hydrogen bonds. In the compound structure, [C5H12N][C18H15N2O2S] (4b), the two ions are primarily associated by an N—H···N hydrogen bond. In the crystal structure, layers parallel to the bc plane are formed by pairs of C—H···O and N—H···S hydrogen bonds and by C—H···π(ring) interactions. A theoretical study reveals that 4a has lower energy than 4b and is more stable. The NBO and DOS studies further confine the liquids’ structural reactivity and electronic properties. The quantum theory of atoms in a molecule (QTAIM) analysis reveals the important non-covalent interactions among the fragments and charge transfer. The global reactivity descriptors indicate their molecular reactivity relationship with the presence of functional groups. The remarkable polarizability (αo) and hyperpolarizability (βo) values indicate their optical and nonlinear optical (NLO) properties. Furthermore, the analysis performed by CrystalExplorer shows the intermolecular interactions and reactive sites between cations and anions in ionic liquids. The 2D fingerprint plots and Hirshfeld surfaces indicate the major interactions of crystals with neighboring elements in crystal packing. For both compounds, the H···H interactions are significantly higher than the other element interactions.

In the present work, we describe the synthesis of new tetrahydroisoquinoline derivatives and the crystal structures of two of them. Density functional theory (DFT) investigations at the B3LYP/6-31+G(d,p) level provided their structural reactivity and nonlinear optical properties. The low HOMO-LUMO gaps (EH-L) suggest a soft nature and higher reactivity, while calculated global reactivity descriptors provide assessments of their reactivity and electronic stability. The calculated natural bonding molecular orbital (NBO) charges show excellent charge separation (charge transfer) and identify the donor and acceptor parts of the molecules. Density of states (DOS) analyses show the newly generated energy states and reduced band gaps, which impart higher conductive properties. For surface reactivity, 3D MESP surfaces are plotted and show electron-rich sites near the nitrogen atoms of the tetrahydroisoquinoline rings. Nonlinear optical (NLO) properties of the crystals are predicted from calculated polarizability (αo) and hyperpolarizability (βo) values. For IVb, the αo and βo values are 415.53 and 1003.44 au. The remarkable value (1003.44 au) of the hyperpolarizability (βo) shows IVb has excellent NLO properties. Structural activity relationship analysis suggests that nitric oxide synthases are better targets for both compounds, and they were further subjected to molecular docking simulations to understand the binding efficiency. In addition, ADMET analyses were carried out to understand the potential activity of the molecules as drug candidates.

Type I and II halogen bonds are well-recognized motifs that commonly occur within crystals. Quantum calculations are applied to examine whether such geometries might occur in their closely related chalcogen bond cousins. Homodimers are constructed of the R1R2C=Y and R1R2Y monomers, wherein Y represents a chalcogen atom, S, Se, or Te; R1 and R2 refer to either H or F. A Type II (T2) geometry wherein the lone pair of one Y is closely aligned with a σ-hole of its partner represents a stable arrangement for all except YH2, although not all such structures are true minima. The symmetric T1 geometry in which each Y atom serves as both electron donor and acceptor in the chalcogen bond is slightly higher in energy for R1R2C=Y, but the reverse is true for R1R2Y. Due to their deeper σ-holes, the latter molecules engage in stronger chalcogen bonds than do the former, with the exception of H2Y, whose dimers are barely bound. The interaction energies rise as the Y atom grows larger: S < Se < Te.

Amphotericin B (AmB) is an antibiotic with a wide spectrum of action and low multidrug resistance, although it exhibits self-aggregation, low specificity, and solubility in aqueous media. An alternative for its oral administration is its encapsulation in polymers modified with bioconjugates. The aim of the present computational research is to determine the affinity between AmB and six bioconjugates to define which one could be more suitable. The CAM-B3LYP-D3/6-31+G(d,p) method was used for all computational calculations. The dimerization enthalpy of the most stable and abundant systems at pH = 7 allows obtaining this affinity order: AmB_1,2-distearoyl-sn-glycerol-3-phosphorylethanolamine (DSPE) > AmB_γ-cyclodextrin > AmB_DSPEc > AmB_retinol > AmB_cholesterol > AmB_dodecanol, where DSPEc is a DSPE analog. Quantum theory of atoms in molecules, the non-covalent interactions index, and natural bond orbital analysis revealed the highest abundance of noncovalent interactions for AmB-DSPE (51), about twice the number of interactions of the other dimers. Depending on the interactions’ strength and abundance of the AmB-DSPE dimer, these are classified as strong: O-H---O (2), N-H---O (3) and weak: C-H---O (25), H---H (18), C-H---C (3). Although the C-H---O hydrogen bond is weak, the number of interactions involved in all dimers cannot be underestimated. Thus, non-covalent interactions drive the stabilization of copolymers, and from our analysis, the most promising candidates for encapsulating are DSPE and γ-cyclodextrin.

The structural, photophysical, and vibrational properties of a new hydroxyphenylamino Meldrum’s acid derivative, 3-((2-hydroxyphenylamino)methylene)-1,5-dioxaspiro[5.5]undecane-2,4-dione (HMD), were studied. The comparison of experimental and theoretical vibrational spectra can help understand basic vibration patterns and provides a better interpretation of IR spectra. The UV–Vis spectrum of HMD was computed using density functional theory (DFT)/B3LYP/6-311 G(d,p) basis set in the gas state, and the maximum wavelength was in accord with the experimental data. The molecular electrostatic potential (MEP) and Hirshfeld surface analysis confirmed O(1)–H(1A)···O(2) intermolecular hydrogen bonds in the HMD molecule. The natural bond orbital (NBO) analysis provided delocalizing interactions between π→π* orbitals and n→σ*/π* charge transfer transitions. Finally, the thermal gravimetric (TG)/differential scanning calorimeter (DSC) and the non-linear optical (NLO) properties of HMD were also reported.

In this work, the hydrogen bonds and halogen bonds in the complexes between Y2CTe (Y = H, F, CH3) and XF (X = F, Cl, Br, I) have been studied by quantum chemical calculations. We found three interesting abnormalities regarding the interactions. Firstly, the strength of halogen bonds increases in the order of IF < BrF < ClF < F2. Secondly, the halogen bonds formed by F2 are very strong, with an interaction energy in the range between −199.8 and −233.1 kJ/mol. Thirdly, all the halogen bonds are stronger than the hydrogen bonds in the systems we examined. All these results are against the general understanding of halogen bonds. These apparent abnormal properties are reconciled with the high polarizability of the Te atom and the strong inducing effect of F on the Te atom of Y2CTe. These findings provide a new perspective on halogen bonds. Additionally, we also proposed bonding distance-based methods to compare the strength of halogen/hydrogen bonds formed between different donor atoms and the same acceptor atom.

Neutral (n) and zwitterionic (z) forms of cysteine monomers are combined in this work to extensively explore the potential energy surfaces for the formation of cysteine dimers in aqueous environments represented by a continuum. A simulated annealing search followed by optimization and characterization of the candidate structures afforded a total of 746 structurally different dimers held together via 80 different types of intermolecular contacts in 2894 individual non-covalent interactions as concluded from Natural Bond Orbitals (NBO), Quantum Theory of Atoms in Molecules (QTAIM) and Non-Covalent Interactions (NCI) analyses. This large pool of interaction possibilities includes the traditional primary hydrogen bonds and salt bridges which actually dictate the structures of the dimers, as well as the less common secondary hydrogen bonds, exotic X⋯Y (X = C, N, O, S) contacts, and H⋯H dihydrogen bonds. These interactions are not homogeneous but have rather complex distributions of strengths, interfragment distances and overall stabilities. Judging by their Gibbs bonding energies, most of the structures located here are suitable for experimental detection at room conditions.

Bonding within the AsF3 crystal is analyzed via quantum chemical methods so as to identify and quantify the pnicogen bonds that are present. The structure of a finite crystal segment containing nine molecules is compared with that of a fully optimized cluster of the same size. The geometries are qualitatively different, with a much larger binding energy within the optimized nonamer. Although the total interaction energy of a central unit with the remaining peripheral molecules is comparable for the two structures, the binding of the peripherals with one another is far larger in the optimized cluster. This distinction of much stronger total binding within the optimized cluster is not limited to the nonamer but repeats itself for smaller aggregates as well. The average binding energy of the cluster rises quickly with size, asymptotically approaching a value nearly triple that of the dimer.

The abundance of applications of alginates in aqueous surroundings created by their interactions with water is a fascinating area of research. In this paper, computational analysis was used to evaluate the conformation, hydrogen bond network, and stabilities for putative intermolecular interactions between alginate dimers and water molecules. Two structural forms of alginate (alginic acid, alg, and sodium alginate, SA) were evaluated for their interactions with water molecules. The density functional theory (DFT-D3) method at the B3LYP functional and the basis set 6-31++G** was chosen for calculating the data. Hydrogen bonds were formed in the Alg-(H2O)n complexes, while the SA-(H2O)n complexes showed an increase in Van der Walls interactions and hydrogen bonds. Moreover, in the SA-(H2O)n complexes, metal-nonmetal bonds existed between the sodium atom in SA and the oxygen atom in water (Na…O). All computational data in this study demonstrated that alginate dimers and water molecules had moderate to high levels of interaction, giving more stability to their complex structure.

A theoretical investigation of mono-methyl derivatives of closo-borate anions of the general form [BnHnCH3]2– (n = 6, 10, 12) and monocarboranes [HCBnHnCH3]− (n = 5, 9, 11) was carried out. An analysis of the main bonding descriptors of exo-polyhedral B–C bonds was performed using the QTAIM (quantum theory of “Atoms in Molecules”), ELF (electron localisation function), NBOs (natural bond orbitals) analyses and several other approaches for the estimation of B–C bond orders (viz. Laplacian bond order (LBO), fuzzy bond order (FBO) and Mayer and Wiberg formalisms). Based on the data obtained on electron density descriptors, it can be concluded that orbital interaction increases with increasing boron cluster size. The present investigation provides a better understanding of exo-polyhedral B–C bond phenomena in boron cluster systems. The data obtained can be used to estimate B–C bond strength, which can be useful for studies devoted to the synthesis and properties of boron cluster systems.

Malondialdehyde (MDA) engages in a triel bond (TrB) with TrX3 (Tr = B and Al; X = H, F, Cl, and Br) in three modes, in which the hydroxyl O, carbonyl O, and central carbon atoms of MDA act as the electron donors, respectively. A H···X secondary interaction coexists with the TrB in the former two types of complexes. The carbonyl O forms a stronger TrB than the hydroxyl O, and both of them are better electron donors than the central carbon atom. The TrB formed by the hydroxyl O enhances the intramolecular H-bond in MDA and thus promotes proton transfer in MDA-BX3 (X = Cl and Br) and MDA-AlX3 (X = halogen), while a weakening H-bond and the inhibition of proton transfer are caused by the TrB formed by the carbonyl O. The TrB formed by the central carbon atom imposes little influence on the H-bond. The BH2 substitution on the central C-H bond can also realise the proton transfer in the triel-bonded complexes between the hydroxyl O and TrH3 (Tr = B and Al).

New Ni (II) and Cu (II) complexes with pyridoxal-semicarbazone were synthesized and their structures were solved by X-ray crystallography. This analysis showed the bis-ligand octahedral structure of [Ni(PLSC-H)2]·H2O and the dimer octahedral structure of [Cu(PLSC)(SO4)(H2O)]2·2H2O. Hirshfeld surface analysis was employed to determine the most important intermolecular interactions in the crystallographic structures. The structures of both complexes were further examined using density functional theory and natural bond orbital analysis. The photocatalytic decomposition of methylene blue in the presence of both compounds was investigated. Both compounds were active toward E. coli and S. aureus, with a minimum inhibition concentration similar to that of chloramphenicol. The obtained complexes led to the formation of free radical species, as was demonstrated in an experiment with dichlorofluorescein-diacetate. It is postulated that this is the mechanistic pathway of the antibacterial and photocatalytic activities. Cyclic voltammograms of the compounds showed the peaks of the reduction of metal ions. A molecular docking study showed that the Ni(II) complex exhibited promising activity towards Janus kinase (JAK), as a potential therapy for inflammatory diseases, cancers, and immunologic disorders.

The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car–Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion—the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.

The reaction of PtCl2 with s-triazine-type ligand (HTriaz) (1:1) in acetone under heating afforded a new [Pt(Triaz)Cl] complex. Single-crystal X-ray diffraction analysis showed that the ligand (HTriaz) is an NNO tridentate chelate via two N-atoms from the s-triazine and hydrazone moieties and one oxygen from the deprotonated phenolic OH. The coordination environment of the Pt(II) is completed by one Cl−1 ion trans to the Pt-N(hydrazone). Hirshfeld surface analysis showed that the most dominant interactions are the H···H, H···C and O···H intermolecular contacts. These interactions contributed by 60.9, 11.2 and 8.3% from the whole fingerprint area, respectively. Other minor contributions from the Cl···H, C···N, N···H and C···C contacts were also detected. Among these interactions, the most significant contacts are the O···H, H···C and H···H interactions. The amounts of the electron transfer from the ligand groups to Pt(II) metal center were predicted using NBO calculations. Additionally, the electronic spectra were assigned based on the TD-DFT calculations.

High level DLPNO–CCSD(T) electronic structure calculations with extended basis sets over B3LYP–D3 optimized geometries indicate that the three methyl groups in caffeine overcome steric hindrance to adopt uncommon conformations, each one placing a C–H bond on the same plane of the aromatic system, leading to the C–H bonds eclipsing one carbonyl group, one heavily delocalized C–N bond constituent of the fused double ring aromatic system, and one C–H bond from the imidazole ring. Deletion of indiscriminate and selective non-Lewis orbitals unequivocally show that hyperconjugation in the form of a bidirectional –CH3 ⇆ aromatic system charge transfer is responsible for these puzzling conformations. The structural preferences in caffeine are exclusively determined by orbital interactions, ruling out electrostatics, induction, bond critical points, and density redistribution because the steric effect, the allylic effect, the Quantum Theory of Atoms in Molecules (QTAIM), and the non-covalent interactions (NCI), all predict wrong energetic orderings. Tiny rotational barriers, not exceeding 1.3 kcal/mol suggest that at room conditions, each methyl group either acts as a free rotor or adopts fluxional behavior, thus preventing accurate determination of their conformations. In this context, our results supersede current experimental ambiguity in the assignation of methyl conformation in caffeine and, more generally, in methylated xanthines and their derivatives.

We unravel the potentialities of resonance Raman spectroscopy to detect ibuprofen in diluted aqueous solutions. In particular, we exploit a fully polarizable quantum mechanics/molecular mechanics (QM/MM) methodology based on fluctuating charges coupled to molecular dynamics (MD) in order to take into account the dynamical aspects of the solvation phenomenon. Our findings, which are discussed in light of a natural bond orbital (NBO) analysis, reveal that a selective enhancement of the Raman signal due to the normal mode associated with the C–C stretching in the ring, νC=C, can be achieved by properly tuning the incident wavelength, thus facilitating the recognition of ibuprofen in water samples.

The current work involves the synthesis of 2,2′-(6-(piperidin-1-yl)-1,3,5-triazine-2,4-diyl)bis(hydrazin-2-yl-1-ylidene))bis(methanylylidene))diphenol 4, characterization, and the DFT studies of the reported compound. The crystal unit cell parameters of 4 are a = 8.1139(2) Å, b = 11.2637(2) Å, c = 45.7836(8) Å. The unit cell volume is 4184.28(15) Å3 and Z = 4. It crystallized in the orthorhombic crystal system and Pbca space group. The O…H, N…H, C…H, H…H and C…C intermolecular contacts which affect the crystal stability were quantitatively analyzed using Hirshfeld calculations. Their percentages were calculated to be 9.8, 15.8, 23.7, 46.4, and 1.6% from the whole contacts occurred in the crystal, respectively. Conformational analysis was performed using DFT calculations for 17 suggested conformers and the most stable conformer was found to be the one which is stabilized by two intramolecular O-H…N hydrogen bonding interactions. This conclusion was further revealed by natural bond orbital calculations.

We report the results of a computational study of the hydrolysis reaction mechanism of N-acetyl-l-aspartyl-l-glutamate (NAAG) catalyzed by glutamate carboxypeptidase II. Analysis of both mechanistic and electronic structure aspects of this multistep reaction is in the focus of this work. In these simulations, model systems are constructed using the relevant crystal structure of the mutated inactive enzyme. After selection of reaction coordinates, the Gibbs energy profiles of elementary steps of the reaction are computed using molecular dynamics simulations with ab initio type QM/MM potentials (QM/MM MD). Energies and forces in the large QM subsystem are estimated in the DFT(PBE0-D3/6-31G**) approximation. The established mechanism includes four elementary steps with the activation energy barriers not exceeding 7 kcal/mol. The models explain the role of point mutations in the enzyme observed in the experimental kinetic studies; namely, the Tyr552Ile substitution disturbs the “oxyanion hole”, and the Glu424Gln replacement increases the distance of the nucleophilic attack. Both issues diminish the substrate activation in the enzyme active site. To quantify the substrate activation, we apply the QTAIM-based approaches and the NBO analysis of dynamic features of the corresponding enzyme-substrate complexes. Analysis of the 2D Laplacian of electron density maps allows one to define structures with the electron density deconcentration on the substrate carbon atom, i.e., at the electrophilic site of reactants. The similar electronic structure element in the NBO approach is a lone vacancy on the carbonyl carbon atom in the reactive species. The electronic structure patterns revealed in the NBO and QTAIM-based analyses consistently clarify the reactivity issues in this system.

Carbonate MCO3 (M = Zn, Cd) can act as both Lewis acid and base to engage in a spodium bond with nitrogen-containing bases (HCN, NHCH2, and NH3) and a chalcogen bond with SeHX (X = F, Cl, OH, OCH3, NH2, and NHCH3), respectively. There is also a weak hydrogen bond in the chalcogen-bonded dyads. Both chalcogen and hydrogen bonds become stronger in the order of F > Cl > OH > OCH3 > NH2 > NHCH3. The chalcogen-bonded dyads are stabilized by a combination of electrostatic and charge transfer interactions. The interaction energy of chalcogen-bonded dyad is less than −10 kcal/mol at most cases. Furthermore, the chalcogen bond can be strengthened through coexistence with a spodium bond in N-base-MCO3-SeHX. The enhancement of chalcogen bond is primarily attributed to the charge transfer interaction. Additionally, the spodium bond is also enhanced by the chalcogen bond although the corresponding enhancing effect is small.

X-ray structural determinations and computational studies were used to investigate halogen interactions in two halogenated oxindoles. Comparative analyses of the interaction energy and the interaction properties were carried out for Br···Br, C-H···Br, C-H···O and N-H···O interactions. Employing Møller–Plesset second-order perturbation theory (MP2) and density functional theory (DFT), the basis set superposition error (BSSE) corrected interaction energy (Eint(BSSE)) was determined using a supramolecular approach. The Eint(BSSE) results were compared with interaction energies obtained by Quantum Theory of Atoms in Molecules (QTAIM)-based methods. Reduced Density Gradient (RDG), QTAIM and Natural bond orbital (NBO) calculations provided insight into possible pathways for the intermolecular interactions examined. Comparative analysis employing the electron density at the bond critical points (BCP) and molecular electrostatic potential (MEP) showed that the interaction energies and the relative orientations of the monomers in the dimers may in part be understood in light of charge redistribution in these two compounds.

A new compounds named 3-4′-bipyrazoles 2 and 3 were synthesized in high chemical yield from a reaction of pyran-2,4-diketone 1 with aryl hydrazines under thermal conditions in MeOH. Compound 2 was unambiguously confirmed by single-crystal X-ray analysis. It crystalizes in a triclinic crystal system and space group P-1. Its crystal structure was found to be in good agreement with the spectral characterizations. With the aid of Hirshfeld calculations, the H…H (54.8–55.3%) and H…C (28.3–29.2%) intermolecular contacts are the most dominant, while the O…H (5.8–6.5%), N…H (3.8–4.6%) and C…C (3.0–4.9%) are less dominant. The compound has a polar nature with a net dipole moment of 6.388 Debye. The BD(2)C31-C32→BD*(2)N4-C34 (27.10 kcal/mol), LP(1)N5→BD*(2)C31-C32 (36.90 kcal/mol), BD(1)C32-C34→BD*(1)C18-C31 (6.78 kcal/mol) and LP(1)N4→BD*(1)N5-C31 (7.25 kcal/mol) are the strongest π→π*, n→π*, σ-σ* and n→σ* intramolecular charge transfer processes, respectively.

The photochemical behavior of the photosensitive first-line anticancer drug vemurafenib (VFB) is of great interest due to the impact of such behavior on its pharmacological activity. In this work, we computationally elucidated the mechanism of the photoinduced release of VFB from the 4,5-dimethoxy-2-nitrobenzene (DMNB) photoprotecting group by employing various density functional theory (DFT)/time-dependent DFT (TD-DFT) approaches. The computational investigations included a comparative assessment of the influence of the position of the photoprotecting group as a substituent on the thermodynamics and kinetics of the photouncaging reactions of two VFB-DMNB prodrugs, namely pyrrole (NP) and sulfonamide (NS). With the aid of the DFT calculations concerning the activation energy barrier (∆G‡), the obtained results suggest that the step of the photoinduced intramolecular proton transfer of the DMNB moiety is not detrimental concerning the overall reaction profile of the photouncaging reaction of both prodrugs. However, the obtained results suggested that the position of the substitution position of the DMNB photoprotecting group within the prodrug structure has a substantial impact on the photouncaging reaction. In particular, the DMNB-Ns-VFB prodrug exhibited a notable increase in ∆G‡ for the key step of ring opining within the DMNB moiety indicative of potentially hindered kinetics of the photouncaging process compared with DMNB-Np-VFB. Such an increase in ∆G‡ may be attributed to the electronic influence of the NP fragment of the prodrug. The results reported herein elaborate on the mechanism of the photoinduced release of an important anticancer drug from photoprotecting groups with the aim of enhancing our understanding of the photochemical behavior of such photosensitive pharmaceutical materials at the molecular level.

Bifurcated halogen bonds are constructed with FBr and FI as Lewis acids, paired with NH3 and NCH bases. The first type considered places two bases together with a single acid, while the reverse case of two acids sharing a single base constitutes the second type. These bifurcated systems are compared with the analogous H-bonds wherein FH serves as the acid. In most cases, a bifurcated system is energetically inferior to a single linear bond. There is a larger energetic cost to forcing the single σ-hole of an acid to interact with a pair of bases, than the other way around where two acids engage with the lone pair of a single base. In comparison to FBr and FI, the H-bonding FH acid is better able to participate in a bifurcated sharing with two bases. This behavior is traced to the properties of the monomers, in particular the specific shape of the molecular electrostatic potential, the anisotropy of the orbitals of the acid and base that interact directly with one another, and the angular extent of the total electron density of the two molecules.

The synthesis, characterization, structural and computational studies of Re(I) tricarbonyl bromo complexes bearing alkyl-substituted 1,4-diazabutadiene ligands, [Re(CO)3(1,4-DAB)Br], where 1,4-DAB = N,N-bis(2,4-dimethylbenzene)-1,4-diazabutadiene, 2,4-Me2DAB (1); N,N-bis(2,4-dimethylbenzene)-2,3-dimethyl-1,4-diazabutadiene, 2,4-Me2DABMe (2); N,N-bis(2,4,6-trimethylbenzene)-1,4-diazabutadiene, 2,4,6-Me3DAB (3); and N,N-bis(2,6-diisopropylbenzene)-1,4-diazabutadiene, 2,6-ipr2DAB (4) are reported. The complexes were characterized by different spectroscopic methods such as FT-IR, 1H-NMR, 13C-NMR, and elemental analyses and their solid-state structures were confirmed by X-ray diffraction. In each complex, the Re(I) centre shows a distorted octahedral shape with a facial geometry of carbonyl groups. The gas phase geometry of the complexes was identified by density functional theory. Interesting intermolecular n…π* interactions of complexes 1 and 3 were investigated by non-covalent interaction index (NCI), and natural bond orbital (NBO) analyses. The intramolecular n…σ*, σ…π*, π…σ* interactions were also studied in complexes 3 and 4.

The Zn(II) complex of salen-like scaffold [Zn(sal)](H2O) was synthesized and characterized by elemental analysis, IR, UV–Vis, and 1H-NMR spectroscopic techniques. The structure of complex was confirmed by single crystal X-ray diffraction studies. In the complex, Zn (II) was placed in the inner N2O2 compartment of the salen scaffold in square planar geometry and crystallized in the monoclinic space group P21/n. DFT and TDDFT calculations were performed to reproduce the experimentally observed structural and spectroscopic (IR and UV–vis) findings. The bonding of the Zn(II) framework in the [Zn(sal)](H2O) complex was explored in depth. The theoretical approaches employed were perturbation theory within the context of the natural bond orbital (NBO) framework, and quantum theory of atoms in molecule (QTAIM) and electron localization function (ELF) analysis. The study begins by delineating the difference between the NBO and QTAIM approaches. This paper thus exhibits the supportive nature of NBO theory and QTAIM in discussion of the bonding in the [Zn(sal)](H2O) complex, when both the methodologies are used in combination.

An organic supramolecular salt hydrate (imidazolium:N-phthalolylglycinate:H2O; IM+-NPG−-HYD) has been examined for its charge-transfer (CT) characteristics. Accordingly, IM+−NPG−−HYD has been characterized thoroughly using various spectroscopic techniques. Combined experimental and quantum chemical studies, along with wave function analysis, were performed to study the non-covalent interactions and their role in CT in the supramolecular salt hydrate. Notably, IM+−NPG−−HYD crystalizes in two configurations (A and B), both of which are held together via non-covalent interactions to result in a three-dimensional CT supramolecular assembly. The through-space CT occurs from NPG− (donor) to IM+ (acceptor), and this was mediated via non-covalent forces. We demonstrated the role of π−π stacking interactions (mixed-stacking donor-acceptor interactions) in the presence of charge-assisted hydrogen bonds in the regulation of CT properties in the self-assembly of the IM+−NPG−−HYD salt hydrate.

Gibbs free energies, based on DFT (Density Functional Theory) calculations, prove that enaminone (2-(anilinemethylidene)cyclohexane-1,3-dione) and ketamine (2-[(phenylimino)-methyl]cyclohexane-1,3-dione) are the most and least stable tautomeric forms of the studied systems, respectively. 1H and 13C NMR spectra prove that 2-(anilinemethylidene)cyclohexane-1,3-diones are the only tautomeric species present in dimethylsulfoxide solution (a very weak signal can be seen only for the p-methoxy derivatives). The zwitterionic character of these enaminones is strengthened by naphthoannulation and by the insertion of the electron-withdrawing substituent into the benzene ring (the latter weakens the intramolecular hydrogen bond in the compound). Substituent and naphtoannulation have no effect on the stability of the studied tautomers. Slight twisting of the benzene ring, with respect to the CArNC plane (seen in the crystalline state), was proven to also take place in vacuum and in solution.

Various types of σ-hole bond complexes were formed with FX, HFY, H2FZ, and H3FT (X = Cl, Br, I; Y = S, Se, Te; Z = P, As, Sb; T = Si, Ge, Sn) as Lewis acid. In order to examine their interactions with a protein, N-methylacetamide (NMA), a model of the peptide linkage was used as the base. These noncovalent bonds were compared by computational means with H-bonds formed by NMA with XH molecules (X = F, Cl, Br, I). In all cases, the A−F bond, which lies opposite the base and is responsible for the σ-hole on the A atom (A refers to the bridging atom), elongates and its stretching frequency undergoes a shift to the red with a band intensification, much as what occurs for the X−H bond in a H-bond (HB). Unlike the NMR shielding decrease seen in the bridging proton of a H-bond, the shielding of the bridging A atom is increased. The spectroscopic changes within NMA are similar for H-bonds and the other noncovalent bonds. The C=O bond of the amide is lengthened and its stretching frequency red-shifted and intensified. The amide II band shifts to higher frequency and undergoes a small band weakening. The NMR shielding of the O atom directly involved in the bond rises, whereas the C and N atoms both undergo a shielding decrease. The frequency shifts of the amide I and II bands of the base as well as the shielding changes of the three pertinent NMA atoms correlate well with the strength of the noncovalent bond.

Complexes were formed pairing FX, FHY, FH2Z, and FH3T (X = Cl, Br, I; Y = S, Se, Te; Z = P, As, Sb; T = Si, Ge, Sn) with NH3 in order to form an A⋯N noncovalent bond, where A refers to the central atom. Geometries, energetics, atomic charges, and spectroscopic characteristics of these complexes were evaluated via DFT calculations. In all cases, the A−F bond, which is located opposite the base and is responsible for the σ-hole on the A atom, elongates and its stretching frequency undergoes a shift to the red. This shift varies from 42 to 175 cm−1 and is largest for the halogen bonds, followed by chalcogen, tetrel, and then pnicogen. The shift also decreases as the central A atom is enlarged. The NMR chemical shielding of the A atom is increased while that of the F and electron donor N atom are lowered. Unlike the IR frequency shifts, it is the third-row A atoms that undergo the largest change in NMR shielding. The change in shielding of A is highly variable, ranging from negligible for FSnH3 all the way up to 1675 ppm for FBr, while those of the F atom lie in the 55−422 ppm range. Although smaller in magnitude, the changes in the N shielding are still easily detectable, between 7 and 27 ppm.

Using the second-order Møller−Plesset perturbation theory (MP2), together with Dunning’s all-electron correlation consistent basis set aug-cc-pVTZ, we show that the covalently bound oxygen atom present in a series of 21 prototypical monomer molecules examined does conceive a positive (or a negative) σ-hole. A σ-hole, in general, is an electron density-deficient region on a bound atom M along the outer extension of the R−M covalent bond, where R is the reminder part of the molecule, and M is the main group atom covalently bonded to R. We have also examined some exemplar 1:1 binary complexes that are formed between five randomly chosen monomers of the above series and the nitrogen- and oxygen-containing Lewis bases in N2, PN, NH3, and OH2. We show that the O-centered positive σ-hole in the selected monomers has the ability to form the chalcogen bonding interaction, and this is when the σ-hole on O is placed in the close proximity of the negative site in the partner molecule. Although the interaction energy and the various other 12 characteristics revealed from this study indicate the presence of any weakly bound interaction between the monomers in the six complexes, our result is strongly inconsistent with the general view that oxygen does not form a chalcogen-bonded interaction.

A coinage-metal bond has been predicted and characterized in the complexes of [1.1.1]propellane (P) and M2/MCl/MCH3 (M = Cu, Ag, and Au). The interaction energy varies between −16 and −47 kcal/mol, indicating that the bridgehead carbon atom of P has a good affinity for the coinage atom. The coinage-metal bond becomes stronger in the Ag < Cu < Au sequence. Relative to M2, both MCl and MCH3 engage in a stronger coinage-metal bond, both -Cl and -CH3 groups showing an electron-withdrawing property. The formation of coinage-metal bonding is mainly attributed to the donation orbital interactions from the occupied C-C orbital into the empty metal orbitals and a back-donation from the occupied d orbital of metal into the empty C-C anti-bonding orbital. In most complexes, the coinage-metal bond is dominated by electrostatic interaction, with moderate contribution of polarization. When P binds simultaneously with two coinage donors, negative cooperativity is found. Moreover, this cooperativity is prominent for the stronger coinage-metal bond.

When an N-base approaches the tetrel atom of TF4 (T = Si, Ge, Sn) the latter molecule deforms from a tetrahedral structure in the monomer to a trigonal bipyramid. The base can situate itself at either an axial or equatorial position, leading to two different equilibrium geometries. The interaction energies are considerably larger for the equatorial structures, up around 50 kcal/mol, which also have a shorter R(T··N) separation. On the other hand, the energy needed to deform the tetrahedral monomer into the equatorial structure is much higher than the equivalent deformation energy in the axial dimer. When these two opposite trends are combined, it is the axial geometry which is somewhat more stable than the equatorial, yielding binding energies in the 8⁻34 kcal/mol range. There is a clear trend of increasing interaction energy as the tetrel atom grows larger: Si < Ge < Sn, a pattern which is accentuated for the binding energies.

Intramolecular hydrogen bond (IMHB) interactions have attracted considerable attention due to their central role in molecular structure, chemical reactivity, and interactions of biologically active molecules. Precise correlations of the strength of IMHB’s with experimental parameters are a key goal in order to model compounds for drug discovery. In this work, we carry out an experimental (NMR) and theoretical (DFT) study of the IMHB in a series of structurally similar o-carbonyl hydroquinones. Geometrical parameters, as well as Natural Bond Orbital (NBO) and Quantum Theory of Atoms in Molecules (QTAIM) parameters for IMHB were compared with experimental NMR data. Three DFT functionals were employed to calculated theoretical parameters: B3LYP, M06-2X, and ωB97XD. O…H distance is the most suitable geometrical parameter to distinguish among similar IMHBs. Second order stabilization energies ΔEij(2) from NBO analysis and hydrogen bond energy (EHB) obtained from QTAIM analysis also properly distinguishes the order in strength of the studied IMHB. ΔEij(2) from NBO give values for the IMHB below 30 kcal/mol, while EHB from QTAIM analysis give values above 30 kcal/mol. In all cases, the calculated parameters using ωB97XD give the best correlations with experimental 1H-NMR chemical shifts for the IMHB, with R2 values around 0.89. Although the results show that these parameters correctly reflect the strength of the IMHB, when the weakest one is removed from the analysis, arguing experimental considerations, correlations improve significantly to values around 0.95 for R2.

The σ-hole tetrel bonds formed by a tetravalent molecule are compared with those involving a π-hole above the tetrel atom in a trivalent bonding situation. The former are modeled by TH4, TH3F, and TH2F2 (T = Si, Ge, Sn) and the latter by TH2=CH2, THF=CH2, and TF2=CH2, all paired with NH3 as Lewis base. The latter π-bonded complexes are considerably more strongly bound, despite the near equivalence of the σ and π-hole intensities. The larger binding energies of the π-dimers are attributed to greater electrostatic attraction and orbital interaction. Each progressive replacement of H by F increases the strength of the tetrel bond, whether σ or π. The magnitudes of the maxima of the molecular electrostatic potential in the two types of systems are not good indicators of either the interaction energy or even the full Coulombic energy. The geometry of the Lewis acid is significantly distorted by the formation of the dimer, more so in the case of the σ-bonded complexes, and this deformation intensifies the σ and π holes.

The synthesis and structural aspects of a new dinuclear silver (I) complex with malonamide type ligand (L) is reported. Each Ag ion in the [Ag2L2(NO3)2]·H2O complex is coordinated to two ligands, L, each acting as a bridged ligand via its two pyridine arms; Ag(I) acts as a connector between them. Two types of Ag-ligands close contacts were detected: Ag–N1, Ag–N4 from the two L units, and Ag–O5, Ag—O6 from the two nitrate anions, wherein both the nitrate ions are inside the cage formed by the [Ag2L2] unit. The coordination geometry around each Ag(I) is a distorted tetrahedron. The [Ag2L2(NO3)2] complex units are connected by weak intermolecular C—H…O interactions. The different intermolecular interactions were quantified using Hirshfeld surface analysis. Using two DFT methods (B3LYP and WB97XD), the nature and strength of the Ag–N and Ag–O interactions were described using atoms in molecules (AIM) and natural bond orbital (NBO) analyses. Topological parameters indicated that the strength of the two Ag–N bonds was similar, while that of the two Ag–O interactions were significantly different. Moreover, the Ag–N interactions have a predominant covalent character, while the Ag–O interactions are mainly ionic. The NBO analysis indicated that the most important anti-bonding Ag-orbital in these interactions has an s-orbital character.

Heterodimers are constructed containing imidazolium and its halogen-substituted derivatives as Lewis acid. N in its sp3, sp2 and sp hybridizations is taken as the electron-donating base. The halogen bond is strengthened in the Cl < Br < I order, with the H-bond generally similar in magnitude to the Br-bond. Methyl substitution on the N electron donor enhances the binding energy. Very little perturbation arises if the imidazolium is attached to a phenyl ring. The energetics are not sensitive to the hybridization of the N atom. More regular patterns appear in the individual phenomena. Charge transfer diminishes uniformly on going from amine to imine to nitrile, a pattern that is echoed by the elongation of the C-Z (Z=H, Cl, Br, I) bond in the Lewis acid. These trends are also evident in the Atoms in Molecules topography of the electron density. Molecular electrostatic potentials are not entirely consistent with energetics. Although I of the Lewis acid engages in a stronger bond than does H, it is the potential of the latter which is much more positive. The minimum on the potential of the base is most negative for the nitrile even though acetonitrile does not form the strongest bonds. Placing the systems in dichloromethane solvent reduces the binding energies but leaves intact most of the trends observed in vacuo; the same can be said of ∆G in solution.

A theoretical exploration about hydrogen bonding in a series of synthetic regioisomeric antitumor tricyclic hydroquinones is presented. The stabilization energy for the intramolecular hydrogen bond (IHB) formation in four structurally different situations were evaluated: (a) IHB between the proton of a phenolic hydroxyl group and an ortho-carbonyl group (forming a six-membered ring); (b) between the oxygen atom of a phenolic hydroxyl group and the proton of an hydroxyalkyl group (seven membered ring); (c) between the proton of a phenolic hydroxyl group with the oxygen atom of the hydroxyl group of a hydroxyalkyl moiety (seven-membered ring); and (d) between the proton of a phenolic hydroxyl group and an oxygen atom directly bonded to the aromatic ring in ortho position (five-membered ring). A conformational analysis for the rotation around the hydroxyalkyl substituent is also performed. It is observed that there is a correspondence between the conformational energies and the IHB. The strongest intramolecular hydrogen bonds are those involving a phenolic proton and a carbonyl oxygen atom, forming a six-membered ring, and the weakest are those involving a phenolic proton with the oxygen atom of the chromenone, forming five-membered rings. Additionally, the synthesis and structural assignment of two pairs of regioisomeric hydroquinones, by 2D-NMR experiments, are reported. These results can be useful in the design of biologically-active molecules.

Intermolecular interactions play a vital role in crystal structures. Therefore, we conducted a topological study, using Hirshfeld surfaces and atom in molecules (AIM) analysis, to decompose and analyze, respectively, the different intermolecular interactions in six hydrazone-diacetyl platinum(II) complexes. Using AIM and natural bond orbital (NBO) analyses, we determined the type, nature, and strength of the interactions. All the studied complexes contain C-H⋯O interactions, and the presence of bond critical points along the intermolecular paths underlines their significance. The electron densities (ρ(r)) at the bond critical points (0.0031–0.0156 e/a03) fall within the typical range for H-bonding interactions. Also, the positive values of the Laplacian of the electron density (∇2ρ(r)) revealed the depletion of electronic charge on the interatomic path, another characteristic feature of closed-shell interactions. The ratios of the absolute potential energy density to the kinetic energy density (|V(r)|/G(r)) and ρ(r) are highest for the O2⋯H15-N3 interaction in [Pt(COMe)2(2-pyCMe=NNH2)] (1); hence, this interaction has the highest covalent character of all the O⋯H intermolecular interactions. Interestingly, in [Pt(COMe)2(H2NN=CMe-CMe=NNH2)] (3), there are significant N-H⋯Pt interactions. Using the NBO method, the second-order interaction energies, E(2), of these interactions range from 3.894 to 4.061 kJ/mol. Furthermore, the hybrid Pt orbitals involved in these interactions are comprised of dxy, dxz, and s atomic orbitals.

The abundance of applications of alginates in aqueous surroundings created by their interactions with water is a fascinating area of research. In this paper, computational analysis was used to evaluate the conformation, hydrogen bond network, and stabilities for putative intermolecular interactions between alginate dimers and water molecules. Two structural forms of alginate (alginic acid, alg, and sodium alginate, SA) were evaluated for their interactions with water molecules. The density functional theory (DFT-D3) method at the B3LYP functional and the basis set 6-31++G** was chosen for calculating the data. Hydrogen bonds were formed in the Alg-(H 2 O) n complexes, while the SA-(H 2 O) n complexes showed an increase in Van der Walls interactions and hydrogen bonds. Moreover, in the SA-(H 2 O) n complexes, metal-nonmetal bonds existed between the sodium atom in SA and the oxygen atom in water (Na…O). All computational data in this study demonstrated that alginate dimers and water molecules had moderate to high levels of interaction, giving more stability to their complex structure.