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.
A post-calculation correction is established for PM7 band gaps of transition-metal oxides. The correction is based on the charge on the metal cation of interest, as obtained from MOPAC PM7 calculations. Application of the correction reduces the average error in the PM7 band gap from ~3 eV to ~1 eV. The residual error after correction is shown to be uncorrelated to the Hartree–Fock method upon which PM7 is based.
Graphical Abstract Comparison between calculated band gaps and experimental band gaps for binary oxides. The orange crosses are for corrected PM7 band gaps. Blue squares are uncorrected values. The orange crosses fall closer to the diagonal dashed line, showing an overall improvement of the accuracy of calculated values
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.
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 the development of osteoarthritis, aggrecan degrades prior to cartilage destruction. Aggrecanase-1 (ADAMTS-4) is considered to be the major enzyme responsible for cleaving the Glu373–Ala374 bond in the interglobular domain of aggrecan in humans. Therefore, inhibitors of ADAMTS-4 have therapeutic potential in the treatment of osteoarthritis. In the present work, we developed a chemical feature based pharmacophore model of ADAMTS-4 inhibitors using the HipHop module within the Catalyst program package in order to elucidate the structure–activity relationship and to carry out in-silico screening. The Maybridge database was screened using Hypo1 as a 3D query, and the best-fit hits that followed Lipinski’s rule of five were subsequently screened to select the compounds. The hit compounds were then docked into the active site of ADAMTS-4, and interactions were visualized to determine the potential lead molecules. After subjecting all of the hits to various screening and filtering processes, 13 compounds were finally evaluated for their in vitro inhibitory activities. This study resulted in the identification of two lead compounds with potent inhibitory effects on ADAMTS-4 activity, with IC50 values of 0.042 μM and 0.028 μM, respectively. These results provide insight into the pharmacophoric requirements for the development of more potent ADAMTS-4 inhibitors.
Graphical Abstract The aggrecan-degrading metalloprotease ADAMTS-4 has been identified as a novel therapeutic target for osteoarthritis. In this work, we used HipHop-based pharmacophore modeling and virtual screening of the Maybridge database to identify novel ADAMTS-4 inhibitors. These novel lead compounds act as potent and specific inhibitors for the ADAMTS-4 enzyme and could have therapeutic potential in the treatment of OA
RET (rearranged during transfection) tyrosine kinase is a promising target for several human cancers. Abt-348, Birb-796, Motesanib and Sorafenib are DFG-out multi-kinase inhibitors that have been reported to inhibit RET activity with good IC50 values. Although the DFG-out conformation has attracted great interest in the design of type II inhibitors, the structural requirements for binding to the RET DFG-out conformation remains unclear. Herein, the DFG-out conformation of RET was determined by homology modelling, the four inhibitors were docked, and the binding modes investigated by molecular dynamics simulation. Binding free energies were calculated using the molecular mechanics/Poisson-Bolzmann surface area (MM/PBSA) method. The trends in predicted binding free affinities correlated well with experimental data and were used to explain the activity difference of the studied inhibitors. Per-residue energy decomposition analyses provided further information on specific interaction properties. Finally, we also conducted a detailed e-pharmacophore modelling of the different RET-inhibitor complexes, explaining the common and specific pharmacophore features of the different complexes. The results reported herein will be useful in future rational design of novel DFG-out RET inhibitors.
Graphical AbstractLeft Ribbon representation of DFG-out RET tyrosine kinase structure showing key residues of RET interacting with inhibitors. Right e-Pharmacophore hypothesis for RET-Abt-348 generated from the complex structure
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.
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.
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.
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 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.
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.
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.
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.
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.
Designed multi-target ligand (DML) is an emerging strategy for the development of new drugs and involves the engagement of multiple targets with the same moiety. In the context of NSAIDs it has been suggested that targeting the thromboxane prostanoid (TP) receptor along with cyclooxygenase-2 (COX-2) may help to overcome cardiovascular (CVS) complications associated with COXIBs. In the present work, azaisoflavones were studied for their COX-2 and TP receptor binding activities using structure based drug design (SBDD) techniques. Flavonoids were selected as a starting point based on their known COX-2 inhibitory and TP receptor antagonist activity. Iterative design and docking studies resulted in the evolution of a new class scaffold replacing the benzopyran-4-one ring of flavonoids with quinolin-4-one. The docking and binding parameters of these new compounds are found to be promising in comparison to those of selective COX-2 inhibitors, such as SC-558 and celecoxib. Owing to the lack of structural information, a model for the TP receptor was generated using a threading base alignment method with loop optimization performed using an ab initio method. The model generated was validated against known antagonists for TP receptor using docking/MMGBSA. Finally, the molecules that were designed for selective COX-2 inhibition were docked into the active site of the TP receptor. Iterative structural modifications and docking on these molecules generated a series which displays optimum docking scores and binding interaction for both targets. Molecular dynamics studies on a known TP receptor antagonist and a designed molecule show that both molecules remain in contact with protein throughout the simulation and interact in similar binding modes.
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.
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.
