The factors that explain the competition between intramolecular NO linkage photoisomerization and NO photorelease in five ruthenium nitrosyl complexes were investigated. By applying DFT-based methods, it was possible to characterize the ground states and lowest triplet potential energy surfaces of these species, and to establish that both photoisomerization and photorelease processes can occur in the lowest triplet state of each species. This work highlights the crucial role of the sideways-bonded isomer, a metastable state also known as the MS2 isomer, in the photochemical loss of NO, while the results obtained also indicate that the population of the triplet state of this isomer is compulsory for both processes and show how photoisomerization and photorelease interfere.
Herein we report a study of the switchable [3]rotaxane reported by Huang et al. (Appl Phys Lett 85(22):5391–5393, 1) that can be mounted to a surface to form a nanomechanical, linear, molecular motor. We demonstrate the application of semiempirical electronic structure theory to predict the average and instantaneous force generated by redox-induced ring shuttling. Detailed analysis of the geometric and electronic structure of the system reveals technical considerations essential to success of the approach. The force is found to be in the 100–200 pN range, consistent with published experimental estimates.
The ternary complexes ML???PyZX2???NH3 (ML?=?CuCl, CuCN, AgCN, and AuCN; Z?=?P, As, and Sb; X?=?H and F) have been investigated with quantum chemical calculations. The results showed that the existence of coordination interaction has a prominent enhancing effect on the strength of pnicogen bonding. Even in ML???PySbH2???NH3, ML???PyAsF2???NH3, and ML???PySbF2???NH3, the pnicogen bond varies from a purely closed-shell interaction to a partially covalent interaction. The coordination interaction results in the enlargement of the σ-hole on the pnicogen atom and thus the enhancement of pnicogen bonding. In addition, the contribution of orbital interaction is also important.
Theoretical calculations for the first tri-iron-based extended metal atom chain (EMAC) molecule are reported. The studied triple-high-spin (S?=?6) complex exhibits ferromagnetic ordering (according to Ising and spin-projection approximations), which renders it unique among all previously prepared and theoretically calculated EMAC compounds. This ordering originates from the prevailing ferromagnetic nearest-neighbor interactions, while the magnetic superexchange between terminal Fe2+ sites is weaker and antiferromagnetic. Calculations indicate that this linear chain system based on a tri-iron core shows potential for the development of spin-frustrated behavior, which could be achieved through rational modification of the equatorial and axial ligands.
A theoretical 1H NMR spectroscopy and thermodynamic analysis of the host–guest inclusion process involving the norfloxacin (NFX) into β-cyclodextrin (β-CD) was carried out. DFT structure and stabilization energies were obtained in both gas and aqueous phases. We could establish that the complex formation is enthalpy driven, and the hydrogen bonds established between NFX and β-CD play a major role in the complex stabilization. Besides, a theoretical 1H NMR analysis has shown to be a supplementary proceeding to predict appropriately the inclusion mode of norfloxacin molecule into the β-CD. In this work, a theoretical study of the NFX@β-CD complex is reported for the first time, seeking a deep understanding of topology and thermodynamics of the inclusion complex formation.
Density functional theory (DFT) calculations are performed to study the hydrogen-bonding in the DMSO-water and DMF-water complexes. Quantitative molecular electrostatic potential (MESP) and atoms-in-molecules (AIM) analysis are applied to quantify the relative complexation of DMSO and DMF with water molecules. The interaction energy of DMSO with water molecules was higher than in DMF-water complexes. The existence of cooperativity effect helps in the strong complex formation. A linear dependence was observed between the hydrogen bond energies EHB, and the total electron densities in the BCP’s of microsolvated complexes which supports the existence of cooperativity effect for the complexation process. Due to the stronger DMSO/DMF and water interaction, the water molecules in the formed complexes have a different structure than the isolated water clusters. NCI analysis shows that the steric area is more pronounced in DMF-water complex than the DMSO-water complex which accounts for the low stability of DMF-water complexes compared to the DMSO-water complex.
Graphical abstract NCI analysis shows that the steric area is more pronounced in DMF-water complex than the DMSO-water complex which accounts for the low stability of DMF-water complexes compared to the DMSO-water complex.
For the first time, the structures, stabilities and electronic properties of alkaline-earth metal doped B44 fullerenes were investigated by means of density functional theory calculations. Our results reveal that M@B44 (M = Ca, Sr, Ba) possess endohedral configurations as their lowest energy structures, whereas the exohedral form is favored when metal is Be or Mg. The large binding energies and sizable HOMO–LUMO gap energies of Ca@B44, Sr@B44 and Ba@B44 suggest the considerable possibility to achieve these novel endohedral borofullerenes experimentally. Born-Oppenheimer molecular dynamics (BO-MD) simulations at various temperatures further confirmed the extreme dynamic stabilities of these endohedral complexes. Their bonding patterns were also analyzed in detail. Finally, we simulated their infrared absorption spectra and 11B nuclear magnetic resonance spectra to help future structural characterization.
In this article, we explore the capacity of formed Schiff base complexes to trap metal atoms or ions, using their aromatic ends. The intrinsic geometry of each complex defines the process of substitution. Two cases were studied; one involving a trans Schiff base complex and the other considering how a salen ligand, with nickel systems traps chromium. We also assessed the nature of the new bonds and the frontier molecular orbitals.
Quantum chemical calculations were performed to investigate the stability of the ternary complexes BeH2···XMH3···NH3 (X?=?F, Cl, and Br; M?=?C, Si, and Ge) and the corresponding binary complexes at the atomic level. Our results reveal that the stability of the XMH3···BeH2 complexes is mainly due to both a strong beryllium bond and a weak tetrel–hydride interaction, while the XMH3···NH3 complexes are stabilized by a tetrel bond. The beryllium bond with a halogen atom as the electron donor has many features in common with a beryllium bond with an O or N atom as the electron donor, although they do exhibit some different characteristics. The stability of the XMH3···NH3 complex is dominated by the electrostatic interaction, while the orbital interaction also makes an important contribution. Interestingly, as the identities of the X and M atoms are varied, the strength of the tetrel bond fluctuates in an irregular manner, which can explained by changes in electrostatic potentials and orbital interactions. In the ternary systems, both the beryllium bond and the tetrel bond are enhanced, which is mainly ascribed to increased electrostatic potentials on the corresponding atoms and charge transfer. In particular, when compared to the strengths of the tetrel and beryllium bonds in the binary systems, in the ternary systems the tetrel bond is enhanced to a greater degree than the beryllium bond.
The structure and stability of various ternary complexes in which an extended aromatic system such as coronene interacts with ions/atoms/molecules on opposite faces of the π-electron cloud were investigated using ab initio calculations. By characterizing the nature of the intermolecular interactions using an energy decomposition analysis, it was shown that there is an interplay between various types of interactions and that there are co-operativity effects, particularly when different types of interactions coexist in the same system.
Graphical abstract Weak OH-π, π-π and van der Waals-π ternary systems are stabilized through dispersion interactions. Cation-π ternary systems are stabilized by through-space electrostatic interactions.
Quantum chemical computations were used for prediction of the structure and color of alizarin complex with alkali metal hydroxides in methanolic solutions. The color prediction relying on the single Gaussian-like band once again proved the usefulness of the PBE0 density functional due to the observed smallest color difference between computed and experimentally derived values. It was found that the alkali metal hydroxide molecules can bind to the two oxygen atoms of both hydroxyl groups of alizarin or to one of these atoms and the oxygen atom from the keto group in a complex with three methanol molecules. This means that two electronic transitions need to be taken into account when considering the spectra of the studied complexes. The resulting bond lengths and angles are correlated with the properties of the alkali metal atoms. The molar mass, the atomic radius, and the Pauling electronegativity of studied metals are quite accurate predictors of the geometric properties of hydroxide complexes with alizarin in methanol solution.
Graphical abstract The spectra of the neutral and monoanionic form of alizarin together with color changes resulting from addition of different metal hydroxides and represented in CIE color space
An experimentally determined structure for human CYP2J2—a member of the cytochrome P450 family with significant and diverse roles across a number of tissues—does not yet exist. Our understanding of how CYP2J2 accommodates its cognate substrates and how it might be inhibited by other ligands thus relies on our ability to computationally predict such interactions using modelling techniques. In this study we present a computational investigation of the binding of arachidonic acid (AA) to CYP2J2 using homology modelling, induced fit docking (IFD) and molecular dynamics (MD) simulations. Our study reveals a catalytically competent binding mode for AA that is distinct from a recently published study that followed a different computational pipeline. Our proposed binding mode for AA is supported by crystal structures of complexes of related enzymes to inhibitors, and evolutionary conservation of a residue whose role appears essential for placing AA in the right site for catalysis.
Graphical Abstract Arachidonic acid docked in the active site of CYP2J2 assumes a catalytically competent binding mode stabilised by hydrogen bonds to Arg117
The conversion of 2-phenylbenzimidazole using o-phenylenediamine and benzaldehyde can be improved significantly under β-cyclodextrin (β-CD). The density functional theory (DFT) method was applied to study the whole process. According to energy parameters (binding energy, deformation energy) and structural deformation, entry models and the reaction process can be pinpointed, with o-phenylenediamine embedding β-CD from a wide rim, and then benzaldehyde passing into the inclusion from the narrow rim. Subsequently, natural bonding orbital (NBO), Mulliken charge, frontier orbital, FuKui function and nuclear magnetic resonance (NMR) methods were employed to reveal the mechanism of electron transfer. The results illustrate that β-CD plays a catalytic role in synthesis reaction mechanism on the secondary side, improving the reactivity and selectivity of the process.
Graphical Abstract Density functional theory study of the effects of β-cyclodextrin in synthesis of 2-phenylbenzimidazole via benzaldehyde and o-phenylenediamine
A mechanistic investigation using Becke3LYP density functional theory (DFT) was carried out on the palladium-catalyzed amidition of bromobenzene and tBu-isocyanide. The whole catalytic cycle consists of five steps: oxidative addition, migratory insertion, anion exchange, reductive elimination, and hydrogen migration. The rate-determining step is oxidative addition, with a small Gibbs free energy of 14.6 kcal mol?1. In the migratory insertion step, tBu-isocyanide provides an important source of carboxy and amino groups to establish the amide group. For anion exchange, path 1a is suggested as the most favorable pathway with the help of the base, and water provides a source of oxygen which is perfectly in line with experimental observations. Finally, in the hydrogen migration step, we illustrate that the six-membered ring path is energetically favored due to the assisting influence of water. In addition, our calculations indicate that using dimethyl sulfoxide as a solvent does not change the rate-determining step.
Density functional theory calculations were carried out to investigate the formation mechanism of the thymine-thymine (6–4) dimer ((6–4)TT), which is one of the main DNA lesions induced by ultraviolet radiation and is closely related to skin cancers. The DNA backbone was found to have nonnegligible effects on the triplet reaction pathway, particularly the reaction steps involving substantial base rotations. The mechanism for the isomerization from (6–4)TT to its Dewar valence isomer (DewarTT) was also explored, confirming the necessity of absorbing a second photon. In addition, the solvation effects were examined and showed considerable influence on the potential energy surface.
Determination of electrophilic and nucleophilic sites of a molecule is the primary task to find the active sites of the lead molecule. In the present study, the active sites of busulfan have been predicted by molecular electrostatic potential surface and Fukui function analysis with the help of dispersion corrected density functional theory. Similarly, the identification of active binding sites of the proteins against lead compound plays a vital role in the field of drug discovery. Rigid and flexible molecular docking approaches are used for this purpose. For rigid docking, Hex 8.0.0 software employing fast Fourier transform (FFT) algorithm has been used. The partial flexible blind docking simulations have been performed with AutoDock 4.2 software; where a Lamarckian genetic algorithm is employed. The results showed that the most electrophilic atoms of busulfan bind with the targets. It is clear from the docking studies that busulfan has inhibition capability toward the targets 12CA and 1BZM.
