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.
In this work, through a docking analysis of compounds from the ZINC chemical library on human β-tubulin using high performance computer cluster, we report new polycyclic aromatic compounds that bind with high energy on the colchicine binding site of β-tubulin, suggesting three new key amino acids. However, molecular dynamic analysis showed low stability in the interaction between ligand and receptor. Results were confirmed experimentally in in vitro and in vivo models that suggest that molecular dynamics simulation is the best option to find new potential β-tubulin inhibitors.
We employ computational methods to investigate the possibility of using electron-donating or electron-withdrawing substituents to reduce the free-energy barriers of the thermal isomerizations that limit the rotational frequencies achievable by synthetic overcrowded alkene-based molecular motors. Choosing as reference systems one of the fastest motors known to date and two variants thereof, we consider six new motors obtained by introducing electron-donating methoxy and dimethylamino or electron-withdrawing nitro and cyano substituents in conjugation with the central olefinic bond connecting the two (stator and rotator) motor halves. Performing density functional theory calculations, we then show that electron-donating (but not electron-withdrawing) groups at the stator are able to reduce the already small barriers of the reference motors by up to 18 kJ mol?1. This result outlines a possible strategy for improving the rotational frequencies of motors of this kind. Furthermore, exploring the origin of the catalytic effect, it is found that electron-donating groups exert a favorable steric influence on the thermal isomerizations, which is not manifested by electron-withdrawing groups. This finding suggests a new mechanism for controlling the critical steric interactions of these motors.
Graphical Abstract The introduction of electron-donating groups in one of the fastest rotary molecular motors known to date is found to reduce the free-energy barriers of the thermal steps that limit the rotational frequencies by up to 18 kJ mol?1.
In this work, we demonstrate that the inclusion of long-range interactions has a significant impact on the estimation of ligand–protein binding energies. Within the scope of the electrostatically embedded adaptation of the molecular fragmentation with conjugated caps (EE-AMFCC) scheme, we unveil the role played by long-range contributions in distinct levels of quantum mechanical calculations. As a prototypical system, we consider ibuprofen coupled to the human serum albumin. In particular, we show that some relevant ligand–residue interaction energies can only be accurately captured in density functional theory (DFT) approaches when the electrostatic background is properly represented by an explicit point charge distribution.
Graphical Abstract (left) The binding site FA3/FA4 of HSA containing the attached IBU. (right) Absolute value of difference between the biding energies calculated including the electrostatic embedding and the energies calculated without the electrostatic embedding using the HF, B3LYP, CAM-B3LYP, and MP2 methodologies
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.
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.
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.
Coarse-grained force field (CGFF) methods were applied to study the self-assembly of sodium dodecyl sulfate with fragrance additives. The CGFF parameters were parameterized and validated using experimental and all-atom simulation data. Direct molecular dynamics simulations were carried out to characterize the initial aggregation, partitioning of fragrances, and chemical potentials of the surfactant and fragrance molecules in aggregates of different sizes. The equilibrium critical micelle concentrations (CMCs) and micelle size distributions, which could not be obtained by direct simulation, were predicted using the calculated chemical potentials in combination with a thermodynamic model. The predicted partitioning of fragrances, CMCs, micelle sizes, and micelle structures agree well with previously reported experimental data.
Magnetic shielding constants for an isolated fullerene C60, cucurbituril CB[9], and the host-guest complex C60@CB[9] were calculated as a function of separation of the monomers. Our results in the gas phase and water indicate a significant variation of the magnetic properties for all atoms of the monomers in the complex and after liberation of fullerene C60 from the interior of the CB[9] cavity. The interaction between the two monomers results in a charge transfer that collaborates with a redistribution of electron density to deshield the monomers.
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.
Calcium complexes with bidentate carbonyl ligands are important in biological systems, medicine and industry, where the concentration of Ca2+ is controlled using chelating ligands. The exchange of two water molecules of [Ca(H2O)6]2+ for one bidentate monosubstituted and homo disubstituted dicarbonyl ligand was investigated using the B3LYP/6-311++G(d,p) method. The ligand substituents NH2, OCH3, OH, CH3, H, F, Cl, CN and NO2 are functional groups with distinct electron-donating and -withdrawing effects that bond directly to the sp2 C atom of the carbonyl group. The geometry, charge and energy characteristics of the complexes were analyzed to help understand the effects of substituents, spacer length and chelation. Coordination strength was quantified in terms of the enthalpy and free energy of the exchange reaction. The most negative enthalpies were calculated for the coordination of bidentate ligands containing three to five methylene group spacers between carbonyls. The chelate effect contribution was analyzed based on the thermochemistry. The electronic character of the substituent modulates the strength of binding to the metal cation, as ligands containing electron-donor substituents coordinate stronger than those with electron-acceptor substituents. This is reflected in the geometric (bond length and chelating angle), electronic (atomic charges) and energetic (components of the total interacting energy) characteristics of the complexes. Energy decomposition analysis (EDA)—an approach for partitioning of the energy into its chemical origins—shows that the electrostatic component of the coordination is predominant, and yields relevant contribution of the covalent term, especially for the electron-withdrawing substituted ligands. The chelate effect of the bidentate ligands was noticeable when compared with substitution by two monodentate ligands.
Graphical abstract The affinity of 18 bidentate carbonyl ligands toward the [Ca(H2O)4]2+ cation is evaluated in terms of energetic, geometric and electronic parameters of the isolated ligands and the substituted aqua complexes. The electronic effects—inductive and mesomeric—intrinsic to the molecular structure of each ligand are found to modulate the strength of the metal-ligand interaction. The effects of polysubstitution, chelation and the length of the alkyl spacers between the anchor points of the ligand are also analyzed.
The present paper reports the analysis of surface decoration on the structural, electronic, and optical properties of (n,0) ZnO nanotubes, performed by means of a density function theory based ab-initio approach. Fe functionalization induced buckling in ZnO nanotubes affects its electronic and optical properties. Increase in Fe functionalization leads to better stability of ZnO nanotube and shows enhanced metallic character. The possibility of its use in optoelectronics has been analyzed in terms of dielectric constant, absorption coefficient, and refractive index. In another observation, the high sensitivity of the HCN molecule for the Fe-incorporated ZnO nanotube suggests it as a potential gas sensor.
In this article, we explore, both theoretically and experimentally, the general reactivity of alkyl hydrogeno-phenylphosphinates with alcohols. We show that alcohol molecules act exclusively as nucleophilic species, and add to alkyl hydrogeno-phenylphosphinates, leading to pentacoordinated intermediates. These intermediates are shown to subsequently competitively undergo alcohol eliminations and/or Berry pseudorotations. This offers several possible routes for racemizations and/or alcohol exchange reactions. Transition standard Gibbs free energies predicted from DFT calculations for the overall alcohol exchange mechanism are shown to be compatible with those experimentally measured in case ethanol reacts with ethyl hydrogeno-phenylphosphinate (134.5～136.0 kJ mol?1 at 78 °C).
Catalytic fields illustrate topology of the optimal charge distribution of a molecular environment reducing the activation energy for any process involving barrier crossing, like chemical reaction, bond rotation etc. Until now, this technique has been successfully applied to predict catalytic effects resulting from intermolecular interactions with individual water molecules constituting the first hydration shell, aminoacid mutations in enzymes or Si→Al substitutions in zeolites. In this contribution, hydrogen to fluorine (H→F) substitution effects for two model reactions have been examined indicating qualitative applicability of the catalytic field concept in the case of systems involving intramolecular interactions.
In an effort to replace the widely used ruthenium metal complexes with low-cost, earth abundant iron complexes as photosensitizers for dye-sensitized solar cell (DSSC) applications, herein we report the computational design of heteroleptic iron complexes (FC1–3) coordinated with benzimidazole-phenylcarbene (C^N) ligands. DFT and TDDFT calculations predicted the stronger σ-donating and π-accepting nature of phenyl carbene ligands substituted with electron-withdrawing CF3, donating –N(CH3)2, and benzothiazine annulation than the imidazole carbene ligands (FC4); consequently, the metal-ligand bond distances and interactions that influence the ordering of charge transfer states with respect to metal centered states are altered in FC1–3 complexes. Detailed analysis based on energy decomposition analysis, spin density distribution analysis, and ab initio ligand field theory parameters were enabled to understand the nature of heteroleptic ligand interactions with the rest of the metal complex. The results from the study shed light on the judicious choice of ligands, as the same non-innocent ligand that is experimentally proven as favorable for Ru-dyes (TC1) is found to be detrimental for Fe-dyes (FC1). Among the complexes studied, the FC3 complex is a promising sensitizer for DSSC with 1,3MLCT energy level well separated from 3,5MC, thereby preventing the deactivation of MLCT. The outcome of the study is therefore an important step toward the development of photosensitizers based on iron metal.
Graphical abstract Potential photosensitzers based on earth-abundant, low cost iron metal have been designed for dye sensitized solar cell applications.
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.
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.
The methylxanthines have attracted interest due to the changes on their biological activities and physicochemical properties in terms of the number and position of the methyl groups present in the xanthine moiety. We report a theoretical study of the influence of the methyl substituent in the basicity and reactivity of xanthine and its methylated derivatives. Our results provide that when the xanthine increases the number of methyl substituents, the gas phase basicity increases (reactivity to proton increases), and the global hardness decreases. The result is in agreement with the maximum hardness principle (MHP) that states, “at equilibrium, chemical systems are as hard as possible” (Pearson, R.G., J. Chem. Educ., 1987, 64, 561–567, and Parr R.G., Chattaraj P.K., J. Am. Chem. Soc. 1991, 113, 1854–1855).
