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.
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.
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.
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.
In this study, the doped defects in nitromethane crystals were investigated using first-principles calculations for the first time. We introduce dopant atoms in the interstitial sites of the nitromethane lattice, aiming to study the effects of element-doping on the structural properties, electronic properties, and sensitivity characteristics. The obtained results show that doped defects obviously affect the neighboring nitromethane molecules. The modification of electronic properties shows that the band gaps are significantly influenced by doped defects. Partial density of states and population analysis further reveal the mechanism for sensitivity control of nitromethane. It is shown that the new electronic states were introduced in the forbidden bands and the doped defects resulted in charge redistributions in the systems.
In this work, we address the effects of molecular doping on the electronic properties of fluorinated and chlorinated silicon nanowires (SiNWs), in comparison with those corresponding to hydrogen-passivated SiNWs. Adsorption of n-type dopant molecules on hydrogenated and halogenated SiNWs and their chemisorption energies, formation energies, and electronic band gap are studied by using density functional theory calculations. The results show that there are considerable charge transfers and strong covalent interactions between the dopant molecules and the SiNWs. Moreover, the results show that the energy band gap of SiNWs changes due to chemical surface doping and it can be further tuned by surface passivation. We conclude that a molecular based ex-situ doping, where molecules are adsorbed on the surface of the SiNW, can be an alternative path to conventional doping.
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 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.
The local and condensed Fukui functions as well as the principle of hard and soft acids and bases were used to study the addition of free radicals to the exocyclic and endocyclic double bonds of seven monocyclic monoterpenes of formula C10H16. The results obtained showed that, in general, the most reactive double bond was the one with the most substituents on the double-bonded carbon atoms, and that the reaction of a double bond with an electrophile is a soft–soft interaction. The effects of substituents on the double-bonded carbon atoms and the stabilization of the monoterpenes were interpreted by invoking hyperconjugated structures, which led us to propose a simple rule: the larger the value of the Fukui function for the double bond, the greater the hyperconjugative stabilization and the susceptibility of the double bond to electrophilic attack. In general, our results are in good accordance with relevant experimental and theoretical results published in the literature.
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.
Present molecular dynamics simulations indicate that the methanol component in a methanol/water mixture is more likely to be trapped in a cyclic peptide nanotube (CPNT), while water molecules tend to be present at the channel mouths as transient guests. Channel water resides mainly between methanol and the CPNT wall, resulting in a distinct decrease in the H-bond number per channel methanol. Six designed CPNTs with different channel diameters and outer surface characteristics all possess distinct selectivity to methanol over water. Of these, the amphipathic 8?×?(AQ)4-CPNT exhibits the best performance. Results in this study provide basic information for the application of a CPNT to enrich methanol from a methanol/water mixture.
The number of hydrogen bonds and detailed information on the interlayer spacing of graphene oxide (GO) confined water molecules were calculated through experiments and molecular dynamics simulations. Experiments play a crucial role in the modeling strategy and verification of the simulation results. The binding of GO and water molecules is essentially controlled by hydrogen bond networks involving functional groups and water molecules confined in the GO layers. With the increase in the water content, the clusters of water molecules are more evident. The water molecules bounding to GO layers are transformed to a free state, making the removal of water molecules from the system difficult at low water contents. The diffuse behaviors of the water molecules are more evident at high water contents. With an increase in the water content, the functional groups are surrounded by fewer water molecules, and the distance between the functional groups and water molecules increases. As a result, the water molecules adsorbed into the GO interlamination will enlarge the interlayer spacing. The interlayer spacing is also affected by the number of GO layers. These results were confirmed by the calculations of number of hydrogen bonds, water state, mean square displacement, radial distribution function, and interlayer spacing of hydrated GO.
Graphical Abstract This work research the interaction between GO functional groups and confined water molecules. The state of water molecules and interlayer spacing of graphene oxide were proved to be related to the number of hydrogen bonds.
The low bending rigidity of graphene facilitates the formation of folds into the structure. This curvature change affects the reactivity and electron transport of the sheet. One novel extension of this is the intercalation of small molecules into these folds. We construct a model incorporating a single-walled carbon nanotube into a sheet of folded graphene. Variational calculus techniques are employed to determine the minimum energy structure and the resulting curves are shown to agree well with molecular dynamics study.
Graphical Abstract Using calculus of variations, the elastic bending energy and van der Waals energy are minimised giving rise to Euler-Lagrange equation for which analytical solutions are derived to determine the optimal curved sturctures of graphene wrapped around carbon nanotubes . Overall agreement between the analytical solutions (with different values of bending rigidities) and results from molecular dynamics simulations (grey) is shown here for (6,6), (8,8) and (10,10) armchair nanotubes, respectively.
High-level ab initio calculations on the complexes between noble gas atoms (He, Ne, Ar, Kr, and Xe) and dihalogen molecules (F2, Cl2, Br2, and I2) reveal trends, both in interaction energies and the energy difference between the linear and T-shaped structures, that can be explained well in terms of dispersion interactions enhanced by polar flattening of the halogens. The partial discrepancies with experimental findings are discussed.
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.
A series of three star-shaped compounds containing both donor (carbazole) and acceptor (2,4,6-triphenyl-1,3,5-triazine) moieties linked through various linking bridges was studied theoretically at the linear response TD-DFT level of theory to describe their absorption and fluorescence spectra. The concept of a localized charge-transfer excited state has been applied successfully to explain the observed strong solvatochromic effect in the emission spectra of the studied molecules, which can be utilized for the fabrication of color tunable solution-processable OLEDs. The concept is in particularly applicable to donor–acceptor species with a C3 symmetry point group where the static dipole moment changes dramatically upon electronic excitation. An important peculiarity of the studied molecules is that they are characterized by non-zero values of the HOMO and LUMO orbitals in the same common part of molecular space that provides a large electric dipole transition moment for both light absorption and emission.
To determine the kinase inhibitory potential of natural products that could be utilized in lung cancer therapy in the near future, a pharmacophore-based activity profiling protocol using parallel pharmacophore-based virtual screening of ZINC—a natural product database—was employed. The work presented here is based on the previously explored fact that pharmacophore-based parallel screening is a reliable in silico protocol to predict the possible biological activities of any compound, or any compound library, by screening it with a number of pharmacophore models. The present study involves ligand-based pharmacophore modeling of various kinases, including EGFR (T790 M), cMET, ErbB2, FGFR and ALK, which are well established targets of normal as well resistant lung cancer. The generated pharmacophore models were then utilized for parallel and cross screening. The profiled molecules for each target were then validated using molecular docking and molecular dynamic simulations. The results show that kinase inhibitory activity profiling of some natural product molecules was successfully achieved.
A dispersion correction is introduced and tested for MNDO. The shift in electron density caused by the interaction between oscillating dipoles in the London picture of dispersion is mimicked by adding a small r?7-dependent attractive nucleus–electron potential to the core Hamiltonian. This potential results in a shift in electron density similar to that used by Feynman to explain dispersion. The resulting parameterized self-consistent and inherently multicenter treatment (MNDO-F) gives good results for CHNO compounds that do not exhibit hydrogen bonds, which MNDO cannot reproduce. This “Feynman” dispersion correction is also applicable to Hartree–Fock and density functional theory.
The lateral pressure exerted on cell membrane is of great importance to signal transduction. Here, we perform molecular dynamics simulation to explore how lateral pressure affects the biophysical properties of lipid bilayer as well as nitrogen distribution in the membrane. Our results show that both physical properties of cell membrane and nitrogen distribution are highly sensitive to the lateral pressure. With the increasing lateral pressure, area per lipid drops and thickness of membrane increases obviously, while nitrogen molecules are more congested in the center of lipid bilayer than those under lower lateral pressure. These results suggest that the mechanism of nitrogen narcosis may be related to the lateral pressure.
Graphical Abstract (a) Schematic illustration of pressure exerted on the membrane. Carbon (cyan), oxygen (red), nitrogen (blue) and phosphorous (tan) in head groups are shown by using big spheres, while the lipid tails are shown as dynamicbonds. Nitrogen molecules are depicted as the yellow balls. Water molecules are omitted for clarity. (b) Probability density of the z direction of nitrogen molecules, where the z axis is normal to the lipid bilayer. The center of the lipid bilayer is fixed at z?=?0 nm
