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.
Unknown force-field parameters for metal organic beryllium complexes used in emitting and electron transporting layers of OLED structures are determined. These parameters can be used for the predictive atomistic simulations of the structure and properties of amorphous organic layers containing beryllium complexes. The parameters are found for the AMBER force field using a relaxed scan procedure and quantum-mechanical DFT calculations of potential energy curves for specific internal (angular) coordinates in a series of three Be complexes (Bebq2; Be(4-mpp)2; Bepp2). The obtained parameters are verified in calculations of some molecular and crystal structures available from either quantum-mechanical DFT calculations or experimental data.
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.
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.
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 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.
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.
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.
The absorption and emission spectra of dichlorvos and the dichlorvos-MAA complex in methanol, water, and chloroform in the molecularly imprinted recognition were investigated systematically. The M06-2X results revealed that: 1) the hydroxyl groups in polar solvents such as methanol and water may markedly influence the weak interactions, and then alter the adsorption and emission spectra; 2) the electronic excitation in absorption spectra of dichlorvos is dominated by the configuration HOMO?→?LUMO, but in the most stable dichlorvos-MAA it becomes the ππ* excitation of HOMO?→?LUMO?+?1; 3) Mulliken charges reveal that dichlorvos almost dissociates to Cl- and a cation in its S1 excitation state; 4) the phosphorescence spectra of dichlorvos-MAA are relatively weak.
Graphical Abstract The absorption and emission spectra of dichlorvos and the dichlorvos-MAA complex in the molecularly imprinted recognition of dichlorvos were investigated systematically in methanol, water, and chloroform as solvents.
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
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.
Phosphatidylinositols and their phosphorylated derivatives, phosphoinositides, play a central role in regulating diverse cellular functions. These phospholipids have been shown to interact with the hydrophobic TH domain of the tumor necrosis factor (TNF)-α-induced protein 8 (TIPE) family of proteins. However, the precise mechanism of interaction of these lipids is unclear. Here we report the binding mode and interactions of these phospholipids in the TH domain, as elucidated using molecular docking and simulations. Results indicate that phosphoinositides bind to the TH domain in a similar way by inserting their lipid tails in the hydrophobic cavity. The exposed head group is stabilized by interactions with critical positively charged residues on the surface of these proteins. Further MD simulations confirmed the dynamic stability of these lipids in the TH domain. This computational analysis thus provides insight into the binding mode of phospholipids in the TH domain of the TIPE family of proteins.
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.
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.
Bond critical points (BCPs) in the quantum theory of atoms in molecules (QTAIM) are shown to be a consequence of the molecular topology, symmetry, and the Poincaré-Hopf relationship, which defines the numbers of critical points of different types in a scalar field. BCPs can be induced by a polarizing field or by addition of a single non-bonded atom to a molecule. BCPs and their associated bond paths are therefore suggested not to be a suitable means of identifying chemical bonds, or even attractive intermolecular interactions.
Polychlorinated dibenzothiophene (PCDT) and polychlorinated thianthrene (PCTA) are sulfur analogues of dioxins, such as polychlorinated dibenzo-p-dioxins and polychlorinated dibenzofurans (PCDD/F). In this work, we present a detailed mechanistic and kinetic analysis of PCDT and PCTA formation from the combustion of 2,4,5-trichlorothiophenol. It is shown that the formation of these persistent organic pollutants is more favourable, both kinetically and thermodynamically, than their analogous dioxin counterparts. This is rationalised in terms of the different influences of the S–H and O–H moieties in the 2,4,5-trichlorothiophenol and 2,4,5-trichlorophenol precursors. Kinetic parameters also indicate that the yield of PCDT should exceed that of PCDD. Finally, we demonstrate here that the degree and pattern of chlorination on the 2,4,5-trichlorothiophenol precursor leads to subtle thermodynamic and kinetic changes to the PCDT/PCTA formation mechanisms.
Graphical abstract Formation mechanisms of persistant organic pollutants, PCDT and PCTA, from 2,4,5-trichlorothiophenol combustion, has been investigated using density functional theory
The structures and energies for the Huisgen 1,3-dipolar cycloaddition reactions of methyl and ethyl azides with some cyclooctynes and dibenzocyclooctynes were computed at the B3LYP/6-311++G(d,p) level. The activation strain model (ASM) and quantitative molecular orbital (MO) theory were used to investigate the reactivity and regiochemistry in these reactions. The energy decomposition analysis (EDA) was used to identify the intrinsic electronic factor that lead to the preferential formation of 1,7-regiochemistry products. The reactivity order agrees with formation of more synchronous transition states and lower distortion energies. For the reaction of N3Met with azacyclooctyne, the 1,7-regiochemistry preference is attributed to a lower FMO gap and a higher contribution of the polarization term of the interaction energy than for the 1,8-transition state. For the reaction with aza-dibenzocyclooctyne, the 1,7-preference is due to a lower strain energy and a more pronounced contribution of the exchange term of the interaction energy.
Graphical Abstract In the reactions between methyl and ethyl azides with azacyclooctynes the regiochemistry is governed by the intrinsic electronic factors
We have provided a straightforward route to screen a series of metal-free doped bulk h-BN as potential visible-light photocatalysts via the first-principle computations. Various nonmetal dopants are considered including Si, P, C, S, Cl, O, and F atoms according to increasing electronegativity. The results show that the introduction of nonmetals leads to small lattice distortions but significant modifications of band structures, electron transition paths and chemical bonding interactions. Generally, all doped h-BN except Si doping have an active response to the visible-light, and dopants with higher electronegativity can significantly narrow the band gaps, which could induce easier optical transition under visible-light excitation. Based on the electronic structures and absorption spectra, three different mechanisms of enhanced visible-light response for the doping effect are proposed. It is expected that F, Cl, and S-doped h-BN could be used as potential efficient visible-light driven photocatalysts. This study could aid in the design of novel efficient h-BN photocatalysts.
