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.
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.
The aldol reaction in the presence of L-proline acting as an organocatalyst is a well-known example of asymmetric synthesis. Many theoretical and experimental studies have been carried out to probe the mechanism of this reaction. In this work, two levels of density functional theory in the gas phase and DMSO were used to elucidate the best pathways for this reaction, with the enamine and enol considered intermediates and L-proline considered either a reactant or a facilitator. The calculations indicated that both intermediates are formed simultaneously in the reaction medium. Interestingly, the formation of the enamine intermediate predominates in DMSO at room temperature, whereas the enol becomes the predominant intermediate upon the addition of water.
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.
A perfectly planar Al13+ cluster (CI) and a quasi-planar Al13+ cluster (CII) have been found for the first time. Both clusters have a triangular core surrounded by a set of ten Al atoms in the form of a ring. These cationic clusters have substantial aromatic character. The planar CI cluster has local antiaromatic patches within global aromatic sea. It is doubly aromatic having both σ and π aromatic character. The quasi-planar CII cluster is also aromatic but it has more σ-delocalization.
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.
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 present study reports the geometries, electronic structures, growth behavior, and stabilities of neutral and ionized copper-doped germanium clusters containing 1–20 Ge atoms within the framework of linear combination of atomic orbitals density functional theory (DFT) under the spin-polarized generalized gradient approximation. It was found that Cu-capped Gen (or Cu-substituted Gen+1) and Cu-encapsulated Gen clusters mostly occur in the ground state at a particular cluster size (n). In order to explain the relative stabilities of the ground-state clusters, parameters such as the average binding energy per atom (BE), the embedding energy (EE), and the fragmentation energy (FE) of the clusters were calculated, and the resulting values are discussed. To explain the chemical stabilities of the clusters, parameters such as the energy gap between the highest occupied and the lowest unoccupied molecular orbitals (the HOMO–LUMO gap), the ionization energy (IP), the electron affinity (EA), the chemical potential (μ), the chemical hardness (η), and the polarizability were calculated, and the resulting values are also discussed. Natural atomic orbital (NAO) and natural bond orbital (NBO) analyses were also used to determine the electron-counting rule that should be applied to the most stable Ge10Cu cluster. Finally, the relevance of the calculated results to the design of Ge-based superatoms is discussed.
Figure Contributions of the valance orbitals of the Ge and Cu atom(s) to the HOMO of the ground-state icosahedral Ge10Cu cluster obtained from NBO analysis. The numbers below the clusters represent the occupancies of the HOMO 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.
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.
The 1A1 ground and the first 1B2 excited states of the methylenecyclopropene (triafulvene) are described by localized wave functions, based on 20 structures valence bond structures. The results are compared to CASSCF(4,4) calculations for both the energetics and the dipole moment. Additional calculations with partial electronic delocalization are presented, and it is shown that the dipole moment modification does not correspond to a situation where the antiaromatic situation prevails (with 4n electrons in the cycle). Part of the analysis uses a “trust factor” that helps to decide if a wave function is appropriate to describe a given state. The trust factor compares the VB wave function to the CASSCF’s with their overlap. Finally, the valence bond density is used to produce density maps that illustrate the electron transfer upon excitation.
Graphical Abstract A projector-based method compares CASSCF wave functions to local wave functions, including Lewis structures as shown in the picture. A “trust factor” (τ) is obtained. Both the ground state and the first excited state of the methylenecyclopropene are discussed
To explore the adsorption mechanism of NO, NH3, N2 on a carbon surface, and the effect of basic and acidic functional groups, density functional theory was employed to investigate the interactions between these molecules and carbon surfaces. Molecular electrostatic potential, Mulliken population analyses, reduced density gradient, and Mayer bond order analyses were used to clarify the adsorption mechanism. The results indicate that van der Waals interactions are responsible for N2 physisorption, and N2 is the least likely to adsorb on a carbon surface. Modification of carbon materials to decorate basic or acidic functional groups could enhance the NH3 physisorption because of hydrogen bonding or electrostatic interactions, however, NO physisorption on a carbon surface is poor. Zig-zag sites are more reactive than armchair sites when these gas molecules absorb on the edge sites of carbon surface.
Mechanisms for the activation of water, ammonia, and other small molecules by the PCcarbeneP nickel pincer complex were studied computationally with the aid of density functional theory. The calculation results indicate that the strongly donating, nucleophilic carbene center can engage in a variety of heterolytic splitting of E?H (E=H, C, N, O) bonds, some of which are reversible. The cleavage of E?H bonds across the Ni=C bond represents a new mode of bond activation by ligand cooperativity in nickel pincer complex. On the basis of the calculations, we also demonstrate that reversible H2 activation across the Ir=C bond via the PCcarbeneP iridium pincer complex was observed in the experiments, while other E?H (E=C, N, O) bonds were not activated. Our calculations are in good agreement with experimental observations and could provide new insights into ligand cooperativity in nickel pincer complexes.
Molecularly imprinted polymers can be anticipated as synthetic imitation of natural antibodies, receptors and enzymes. In case of successful imprinting the selectivity and affinity of the imprint for substrate molecules are comparable with those of natural counterparts. The selection of the optimal functional monomer, monomer/template ratio as well as choosing of polymerization solvent is crucial determinants of the successful imprinting. In the present study the simulation approach to the development of molecular imprinting polymers for the extraction of new protein kinase ATP-competitive inhibitors is presented. By imprinting tri-O-acetyladenosine into polymer matrix the synthetic reproduction of adenosine triphosphate binding site to protein kinases can be fabricated and further used for adenosine triphosphate analogs screening in different sources. The optimized geometrical structure and energy of the pre-polymerization complexes of tri-O-acetyladenosine (template) with three different monomers—methacrylic acid, 3-vinyl benzoic acid and acrylamide in vacuum were calculated using hybrid quantum mechanical/molecular mechanical (QM/MM) approach. These calculations demonstrate that methacrylic acid forms the most stable complex with template, the next is 3-vinyl benzoic acid complex and the third—acrylamide one. The bond energies of the complexes are shown to increase monotonically as more monomers are linked to the template. The same conclusions are made from purely quantum self-consistent field calculations of pre-polymerization complex energy and structure. Hybrid calculation is shown to be effective and can substantially accelerate the development of the imprinting technology.
Vitamin C is one of the most abundant exogenous antioxidants in the cell, and it is of the utmost importance to elucidate its mechanism of action against radicals. In this study, the reactivity of vitamin C toward OH and \( {HO}_2/{O}_2^{-} \) radicals in aqueous medium was analyzed by ab initio molecular dynamics using CPMD code. The simulations led to results similar to those of static studies or experiments for the pair of \( {HO}_2/{O}_2^{-} \) radicals but bring new insights for the reactivity with hydroxyl radical: the reaction takes place before the formation of an adduct and consists of two steps: first an electron is transferred to hydroxyl radical and then the ascorbyl radical loses a proton.
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.
Using density functional theory (DFT) and molecular dynamics (MD), we studied the interaction of a titanium atom with a half of a C60 fullerene (i.e., C30), formed from the corannulene structure with a pentagonal base. We considered atmospheric pressure and 300 K. We found that the most stable adsorption of the titanium atom on C30 occurs in the concave surface of the molecule. Afterward, we investigated the interaction of the system C30-titanium with carbon monoxide and carbon dioxide molecules, respectively. We found that each of these molecules is chemisorbed, with no dissociation. The value of the adsorption energy for the carbon monoxide molecule varies from ?0.897 to ?1.673 eV, and for the carbon dioxide molecule, it is between ?1.065 and ?1.274 eV. These values depend on the initial orientation of these molecules with respect to TiC30.
The family of human epidermal growth factor receptors (HER) is involved in tumor cell growth. Homodimerization and heterodimerization of the HER family are important for activation of these receptors. The structures of homodimer conformation are well characterized, while the structures of heterodimer conformations, especially between HER1 and HER2, are not completely understood. In this study, two models of possible asymmetric HER1/HER2 kinase domains were built. Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) methods were applied to examine the possibility of these two-heterodimer interactions. From our results, it could be concluded that the HER2 kinase domain prefers to serve as the receiver rather than the activator. Key binding residues of this dimer complex at N lobe of HER2 is ALA683 and at C lobe of HER1 are GLU914, GLU917, and ASP930. This study will be useful in allowing us to predict and be able to control activity of this enzyme in disease in the future.
Graphical abstract A model of the asymmetric dimer of HER1–HER2 heterodimer with key intereaction residues in (a) HER1A and (b) HER2R by molecular dynamic simulation