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.
Ecdysone receptor (EcR) is a significant target in the identification of new environmentally friendly pesticides. There are two types of ecdysone agonists: steroidal ecdysone agonists and dibenzoylhydrazines (DBHs). In this study, various modeling methods (homology modeling, molecular docking, MD simulation, binding free energy calculation, and per-residue binding free energy decomposition) were utilized to study the different binding mechanisms of two types of ecdysone agonists. Our theoretical results indicated that the relative binding potencies of DBHs can be ranked sufficiently accurately using the MOE docking method. However, MM/PBSA calculations more accurately predicted the binding affinities between steroidal ecdysone agonists and EcR-LBD. To identify the key residues involved in ecdysone agonist binding, the binding free energy (ΔGBind) was decomposed into the energy contributions of individual residues. The results revealed that nine residues—Ile339, Thr343, Met380, Met381, Tyr403, Tyr408, Asp419, Gln503, and Asn504—determined the binding affinities of the DBHs. Glu309, Met342, Arg383, Arg387, and Leu396 were important influences on the binding affinities of the steroidal ecdysone agonists.
Nowadays, breast cancer is one of the most widespread malignancies in women, and the second leading cause of cancer death among women. The progesterone receptor (PR) is one of the treatment targets in breast cancer, and can be blocked with selective progesterone receptor modulators (SPRMs). Since administration of chemical drugs can cause serious side effects, and patients, especially those undergoing long-term treatment, can suffer harmful consequences, there is an urgent need to discover novel potent drugs. Large-scale structural diversity is a feature of natural compounds. Accordingly, in the present study, we selected a library of 20,000 natural compounds from the ZINC database, and screened them against the PR for binding affinity and efficacy. In addition, we evaluated the pharmacodynamics and ADMET properties of the compounds and performed molecular docking. Moreover, molecular dynamics (MD) simulation was carried out in order to examine the stability of the protein. In addition, principal component analysis (PCA) was performed to study the motions of the protein. Finally, the MMPBSA method was applied in order to estimate the binding free energy. Our docking results reveal that compounds ZINC00936598, ZINC00869973 and ZINC01020370 have the highest binding energy into the PR binding site, comparable with that of Levonorgestrel (positive control). Moreover, RMSD, RMSF, Rg and H-bond analysis demonstrate that the lead compounds preserve stability in complex with PR during simulation. Our PCA analysis results were in accordance with MD results and the binding free energies support the docking results. This study paves the way for discovery of novel drugs from natural sources and with optimal efficacy, targeting the PR.
The main purpose of this study was to address the performance of virtual screening methods based on ligands and the protein structure of acetylcholinesterase (AChE) in order to retrieve novel human AChE (hAChE) inhibitors. In addition, a protocol was developed to identify novel hit compounds and propose new promising AChE inhibitors from the ZINC database with 10 million commercially available compounds. In this sense, 3D similarity searches using rapid overlay of chemical structures and similarity analysis through comparison of electrostatic overlay of docked hits were used to retrieve AChE inhibitors from collected databases. Molecular dynamics simulation of 100 ns was carried out to study the best docked compounds from similarity searches. Some key residues were identified as crucial for the dual binding mode of inhibitor with the interaction site. All results indicated the relevant use of EON and docking strategy for identifying novel hit compounds as promising potential anticholinesterase candidates, and seven new structures were selected as potential hAChE inhibitors.
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.
Bone morphogenetic proteins (BMPs) are a family of more than 30 ligands and several receptors, such as activin like kinases (ALKs) and bone morphogenetic protein receptor (BMPR). Physiological significance of these proteins lies in their prominent role during homeostasis, apoptosis, tissue remodeling, embryonic patterning, and normal development. Fibrodysplasia ossificans progressive (FOP) is one among several other diseases caused by impaired BMP signaling. FOP is caused by the pathogenicity of activating mutation of ALK2. In order to treat FOP, a search for good inhibitors of ALK2 based on dorsomorphin and LDN substitution, which in essence is a ligand based search of inhibitors, is in progress. Contributing to this area of research we identified several lead molecules based on protein structure using virtual screening. After virtual screening of a huge library of small molecules and ab initio calculation of selected molecules for drug efficacy, we did molecular dynamic simulation of lead molecules and protein complexes. We identified five potential drug molecules that show very stable binding on the same binding site as LDN-213844. We also ranked these lead molecules based on MM/PBSA binding energy. This study provides a basis to think beyond the pyrimidine nucleus of dorsomorphin/LDN and design new chemical derivatives for effective treatment of FOP.
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
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.
DSS1 is a small acidic intrinsically disordered protein (IDP) that can fold upon binding with PCID2 TREX-2. The resulting complex plays a key role in mRNA export. However, the binding mechanism between DSS1 and PCID2 is unsolved. Here, three independent 500-ns molecular dynamics (MD) simulations were performed to study the DSS1–PCID2 binding mechanism by comparing apo-PCID2 and bound PCID2. The results show that the conformational variation of bound PCID2 is smaller than that of apo-PCID2, especially in the binding domain of two helices (helix IV and VIII). The probability of coil formation between helix III and helix IV of bound PCID2 increases, and a short anti-parallel β-sheet forms upon DSS1 binding. The decomposition of binding free energy into protein and residue pairs suggests that electrostatic and hydrophobic interactions play key roles in the recognition between DSS1 and PCID2. There is a hydrophobic core of seven residues in DSS1 favorable to the binding of PCID2. These analytical methods can be used to reveal the recognition mechanisms of other IDPs and their partners.
The activation of human epidermal growth factor receptor (hEGFR) involves a large conformational change in its soluble extracellular domains (sECD, residues 1–620), from a tethered to an extended conformation upon binding of ligands, such as EGF. It has been reported that this dynamic process is pH-dependent, that is, hEGFR can be activated by EGF at high pH to form an extended dimer but remains as an inactive monomer at low pH. In this paper, we perform all-atom molecular dynamics (MD) simulations starting from the tethered conformation of sECD:EGF complex, at pH 5.0 and 8.5, respectively. Simulation results indicate that sECD:EGF shows different dynamic properties between the two pHs, and the complex may have a higher tendency of activation at pH 8.5. Twenty residues, including 13 histidines, in sECD:EGF have different protonation states between the two pHs (calculated by the H++ server). The charge distribution at pH 8.5 is more favorable for forming an extended conformation toward the active state of sECD than that at pH 5.0. Our study may shed light on the mechanism of pH dependence of hEGFR activation.
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.
Odorant binding proteins (OBPs) are important in insect olfactory recognition. These proteins bind specifically to insect semiochemicals and induce their seeking, mating, and alarm behaviors. Molecular docking and molecular dynamics simulations were performed to provide computational insight into the interaction mode between AgamOBP7 and novel (E)-β-farnesene (EBF) analogues with an aromatic ring. The ligand-binding cavity in OBP7 was found to be mostly hydrophobic due to the presence of several nonpolar residues. The interactions between the EBF analogues and the hydrophobic residues in the binding cavity increased in strength as the distance between them decreased. The EBF analogues with an N-methyl formamide or ester linkage had higher docking scores than those with an amide linkage. Moreover, delocalized π–π and electrostatic interactions were found to contribute significantly to the binding between the ligand benzene ring and nearby protein residues. To design new compounds with higher activity, four EBF analogues D1–D4 with a benzene ring were synthesized and evaluated based on their docking scores and binding affinities. D2, which had an N-methyl formamide group linkage, exhibited stronger binding than D1, which had an amide linkage. D4 exhibited particularly strong binding due to multiple hydrophobic interactions with the protein. This study provides crucial foundations for designing novel EBF analogues based on the OBP structure.
As one of the most important antiviral drugs against 2009 influenza A (H1N1), will zanamivir be effective for the possible drug resistant mutants? To answer this question, we combined multiple molecular dynamics simulations and molecular mechanics generalized Born surface area (MM-GBSA) calculations to study the efficiency of zanamivir over the most frequent drug-resistant strains of neuraminidase including R293K, R152K, E119A/D and H275Y mutants. The calculated results indicate that the modeled mutants of the 2009-H1N1 strains except H275Y will be significantly resistant to zanamivir. The resistance to zanamivir is mainly caused by the loss of polar interactions. The identified potential resistance sites in this study will be useful for the development of new effective anti-influenza drugs and to avoid the occurrence of the state without effective drugs to new mutant influenza strains.
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.
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.
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.
