Identification of novel allosteric modulator binding sites in NMDA receptors: A molecular modeling study

The dysfunction of N-methyl-d-Aspartate receptors (NMDARs), a subtype of glutamate receptors, is correlated with schizophrenia, stroke, and many other neuropathological disorders. However, not all NMDAR subtypes equally contribute towards these disorders. Since NMDARs composed of different GluN2 subunits (GluN2A-D) confer varied physiological properties and have different distributions in the brain, pharmacological agents that target NMDARs with specific GluN2 subunits have significant potential for therapeutic applications. In our previous research, we have identified a family of novel allosteric modulators that differentially potentiate and/or inhibit NMDARs of differing GluN2 subunit composition. To further elucidate their molecular mechanisms, in the present study, we have identified four potential binding sites for novel allosteric modulators by performing molecular modeling, docking, and in silico mutations. The molecular determinants of the modulator binding sites (MBS), analysis of particular MBS electrostatics, and the specific loss or gain of binding after mutations have revealed modulators that have strong potential affinities for specific MBS on given subunits and the role of key amino acids in either promoting or obstructing modulator binding. These findings will help design higher affinity GluN2 subunit-selective pharmaceuticals, which are currently unavailable to treat psychiatric and neurological disorders.     

The role of LasR active site amino acids in the interaction with the Acyl Homoserine Lactones (AHLs) analogues: A computational study

The present work combines molecular docking calculations, 3D-QSAR, molecular dynamics simulations and free binding energy calculations (MM/PBSA and MM/GBSA) in a set of 28 structural analogues of acyl homoserine lactones with Quorum Sensing antagonist activity. The aim of this work is to understand how ligand binds and is affected by the molecular microenvironment in the active site of the LasR receptor for pseudomonas aeruginosa. We also study the stability of the interaction to find key structural characteristics that explain the antagonist activities of this set of ligands. This information is relevant for the rational modification or design of molecules and their identification as powerful LasR modulators.The analysis of molecular docking simulations shows that the 28 analogues have a similar binding mode compared to the native ligand. The carbonyl groups belonging to the lactone ring and the amide group of the acyl chain are oriented towards the amino acids forming hydrogen bond like interactions. The difference in antagonist activity is due to location and orientation of the LasR side chains within the hydrophobic pocket in its binding site. Additionally, we carried out molecular dynamics simulations to understand the conformational changes in the ligand-receptor interaction and the stability of each complex. Results show a direct relationship among the interaction energies of the ligands and the activities as an antagonist of the LasR receptor.     

Investigating detailed interactions between novel PAR1 antagonist F16357 and the receptor using docking and molecular dynamic simulations

Currently, Vorapaxar is the only recently FDA-approved antiplatelet drug targeting Protease-activated receptor 1 (PAR1). However, a novel antagonist, F16357, has been shown to prevent painful bladder syndrome, also known as interstitial cystitis (IC). Unfortunately, there is no high resolution structure of the F16357-receptor complex, hindering its optimization as a therapeutic agent. In this study, we used docking and molecular dynamic (MD) simulations to investigate the detailed interactions between F16357 and PAR1 at a molecular level. The recently solved crystal structure of human PAR1 complexed with Vorapaxar was used in our docking of F16357 into the binding pocket of the receptor. To enhance binding pose selection, F16357 was docked first without constraints and then with a positional constraint to invert its orientation to become similar to that of Vorapaxar. The three systems, with crystal Vorapaxar, F16357 and an inverted F16357, were subjected to 3.0 μs MD simulations. The MM-GBSA binding energy analysis showed that F16357 binds more strongly in a pose obtained from an unrestrained docking than in the inverted pose from a restrained docking; and Vorapaxar binds more strongly than F17357. This ordering is consistent with the experimental pIC₅₀ values. Our structural data showed subtle changes in the binding pose between Vorapaxar and F16357. Transmembrane helices 1, 2, 5, and 7 were most significantly affected; most notably a large kink at F279^5.47 in TM helix 5 of the Vorapaxar complex was completely absent in the F16357 complex. The results of this study facilitate the future development of other therapeutic PAR1 antagonists.     

Modeling the protonation states of β-secretase binding pocket by molecular dynamics simulations and docking studies

β-secretase (BACE1) is an aspartyl protease that processes the β-amyloid peptide in the human brain in patients with Alzheimer’s disease. There are two catalytic aspartates (ASP32 and ASP228) in the active domain of BACE1. Although it is believed that the net charge of the Asp dyad is −1, the exact protonation state still remains a matter of debate. We carried out molecular dynamic (MD) simulations for the four protonation states of BACE1 proteins. We applied Glide docking studies to 21 BACE1 inhibitors against the MD extracted conformations. The dynamic results infer that the protein/ligand complex remains stable during the entire simulation course for HD32D228 model. The results show that the hydrogen bonds between the inhibitor and the Asp dyad are maintained in the 10,000th ps snapshot of HD32D228 model. Our results also reveal the significant loop residues in maintaining the active binding conformation in the HD32D228 model. Molecular docking results show that the HD32D228 model provided the best enrichment factor score, suggesting that this model was able to recognize the most active compounds. Our observations provide an evidence for the preference of the anionic state (HD32D228) in BACE1 binding site and are in accord with reported computational data. The protonation state study would provide significant information to assign the correct protonation state for structure-based drug design and docking studies targeting the BACE1 proteins as a tactic to develop potential AD inhibitors.     

Structural dynamics of Casein Kinase I (CKI) from malarial parasite Plasmodium falciparum (Isolate 3D7): Insights from theoretical modelling and molecular simulations

The protein kinases (PKs), belonging to serine/threonine kinase (STKs), are important drug targets for a wide spectrum of diseases in human. Among protein kinases, the Casein Kinases (CKs) are vastly expanded in various organisms, where, the malarial parasite Plasmodium falciparum possesses a single member i.e., PfCKI, which can phosphorylate various proteins in parasite extracts in vitro condition. But, the structure-function relationship of PfCKI and dynamics of ATP binding is yet to be understood. Henceforth, an attempt was made to study the dynamics, stability, and ATP binding mechanisms of PfCKI through computational modelling, docking, molecular dynamics (MD) simulations, and MM/PBSA binding free energy estimation. Bi-lobed catalytic domain of PfCKI shares a high degree of secondary structure topology with CKI domains of rice, human, and mouse indicating co-evolution of these kinases. Molecular docking study revealed that ATP binds to the active site where the glycine-rich ATP-binding motif (G16-X-G18-X-X-G21) along with few conserved residues plays a crucial role maintaining stability of the complex. Structural superposition of PfCKI with close structural homologs depicted that the location and length of important loops are different, indicating the dynamic properties of these loops among CKIs, which is consistent with principal component analysis (PCA). PCA displayed that the overall global motion of ATP-bound form is comparatively higher than that of apo form. The present study provides insights into the structural features of PfCKI, which could contribute towards further understanding of related protein structures, dynamics of catalysis and phosphorylation mechanism in these important STKs from malarial parasite in near future.     

Binding mode analyses and pharmacophore model development for sulfonamide chalcone derivatives, a new class of alpha-glucosidase inhibitors

Sulfonamide chalcone derivatives are a new class of non-saccharide compounds that effectively inhibit glucosidases which are the major targets in the treatment of Type 2 diabetes and HIV infection. Our aim is to explore their binding mode of interaction at the active site by comparing with the sugar derivatives and to develop a pharmacophore model which would represent the critical features responsible for alpha-glucosidase inhibitory activity. The homology modeled structure of Saccharomyces cerevisiae alpha-glucosidase was built and used for molecular docking of non-sugar/sugar derivatives. The validated docking results projected the crucial role of NH group in the binding of sugar/non-sugar derivatives to the active site. Ligplot analyses revealed that Tyr71, and Phe177 form hydrophobic interactions with sugar/non-sugar derivatives by holding the terminal glycosidic ring mimics. Molecular dynamic (MD) simulation studies were performed for protein alone and with chalcone derivative to prove its binding mechanism as shown by docking/Ligplot results. It would also help to substantiate the homology modeled structure stability. With the knowledge of the crucial interactions between ligand and protein from docking and MD simulation studies, features for pharmacophore model development were chosen. The CATALYST/HipHop was used to generate a five featured pharmacophore model with a training set of five non-sugar derivatives. As validation, all the crucial features of the model were perfectly mapped onto the 3D structures of the sugar derivatives as well as the newly tested non-sugar derivatives. Thus, it can be useful in virtual screening for finding new non-sugar derivatives as alpha-glucosidase inhibitors.     

Binding mode analyses and pharmacophore model development for sulfonamide chalcone derivatives,a new class of α-glucosidase inhibitors

Sulfonamide chalcone derivatives are a new class of non-saccharide compounds that effectively inhibit glucosidases which are the major targets in the treatment of Type 2 diabetes and HIV infection. Our aim is to explore their binding mode of interaction at the active site by comparing with the sugar derivatives and to develop a pharmacophore model which would represent the critical features responsible for α-glucosidase inhibitory activity. The homology modeled structure of Saccharomyces cerevisiae α-glucosidase was built and used for molecular docking of non-sugar/sugar derivatives. The validated docking results projected the crucial role of NH group in the binding of sugar/non-sugar derivatives to the active site. Ligplot analyses revealed that Tyr71, and Phe177 form hydrophobic interactions with sugar/non-sugar derivatives by holding the terminal glycosidic ring mimics. Molecular dynamic (MD) simulation studies were performed for protein alone and with chalcone derivative to prove its binding mechanism as shown by docking/Ligplot results. It would also help to substantiate the homology modeled structure stability. With the knowledge of the crucial interactions between ligand and protein from docking and MD simulation studies, features for pharmacophore model development were chosen. The CATALYST/HipHop was used to generate a five featured pharmacophore model with a training set of five non-sugar derivatives. As validation, all the crucial features of the model were perfectly mapped onto the 3D structures of the sugar derivatives as well as the newly tested non-sugar derivatives. Thus, it can be useful in virtual screening for finding new non-sugar derivatives as α-glucosidase inhibitors.     

Conformational and energy evaluations of novel peptides binding to dengue virus envelope protein

Effective novel peptide inhibitors which targeted the domain III of the dengue envelope (E) protein by blocking dengue virus (DENV) entry into target cells, were identified. The binding affinities of these peptides towards E-protein were evaluated by using a combination of docking and explicit solvent molecular dynamics (MD) simulation methods. The interactions of these complexes were further investigated by using the Molecular Mechanics-Poisson Boltzmann Surface Area (MMPBSA) and Molecular Mechanics Generalized Born Surface Area (MMGBSA) methods. Free energy calculations of the peptides interacting with the E-protein demonstrated that van der Waals (vdW) and electrostatic interactions were the main driving forces stabilizing the complexes. Interestingly, calculated binding free energies showed good agreement with the experimental dissociation constant (K_d) values. Our results also demonstrated that specific residues might play a crucial role in the effective binding interactions. Thus, this study has demonstrated that a combination of docking and molecular dynamics simulations can accelerate the identification process of peptides as potential inhibitors of dengue virus entry into host cells.     

Prediction of the binding mode of biarylpropylsulfonamide allosteric AMPA receptor modulators based on docking, GRID molecular interaction fields and 3D-QSAR analysis

A novel approach of combining flexible molecular docking, GRID molecular interaction fields, analysis of ligand-protein hydrogen bond interactions, conformational energy penalties and 3D-QSAR analysis was used to propose a binding mode in the dimer interface of the iGluR2 receptor for the biarylpropylsulfonamide class of positive allosteric AMPA modulators. Possible binding poses were generated by flexible molecular docking. GRID molecular interaction fields of the binding site, ligand-protein hydrogen bonding interactions and conformational energy penalties were used to select the most likely binding mode. The selected binding poses were subjected to a 3D-QSAR analysis using previously published activity data. The resulting model (2 LVs, R2=0.89, q2=0.61) predicted the activities of the compounds in the test set with a standard deviation on error of prediction of 0.17. The proposed binding mode was validated by interpretation of the PLS-coefficient regions from the 3D-QSAR analysis in terms of interactions between the receptor and the modulators.     

Theoretical studies on binding modes of copper-based nucleases with DNA

In the present work, molecular simulations were performed for the purpose of predicting the binding modes of four types of copper nucleases (a total 33 compounds) with DNA. Our docking results accurately predicted the groove binding and electrostatic interaction for some copper nucleases with B-DNA. The intercalation modes were also reproduced by “gap DNA”. The obtained results demonstrated that the ligand size, length, functional groups and chelate ring size bound to the copper center could influence the binding affinities of copper nucleases. The binding affinities obtained from the docking calculations herein also replicated results found using MM-PBSA approach. The predicted DNA binding modes of copper nucleases with DNA will ultimately help us to better understand the interaction of copper compounds with DNA.     

In silico identification and computational analysis of the nucleotide binding site in the C-terminal domain of Hsp90

Hsp90 contains two distinct Nucleotide Binding Sites (NBS), in its N-terminal domain (NTD) and C-terminal domain (CTD), respectively. The NTD site belongs to the GHKL super-family of ATPases and has been the subject of extensive characterization. However, a structure of the nucleotide-bound form of CTD is still unavailable. In this study molecular modeling was employed to incorporate experimental data using partial constructs of the CTD, from work published by many research groups, onto existing structural models of its apo- form. Our attempts to locate potential nucleotide ligand-binding sites or cavities yielded one major candidate—a structurally unconventional site—exhibiting the requisite shape and volume for accommodation of tri-phosphate nucleotides. Its structure was refined by molecular dynamics (MD)-based techniques. We reproducibly docked the Mg²⁺ complexed form of ATP, GTP, CTP, TTP and UTP to this putative NBS. These docking simulations and calculated ligand-binding scores are in general agreement with published data about experimentally measured binding to the CTD. The overall pattern of interactions between residues lining the site and docked nucleotides is conserved and broadly similar to that of other nucleotide-binding sites. Our docking simulations suggest that nucleotide binding stabilizes the only structurally labile region, thereby providing a rationale for the increased resistance to thermal denaturation and proteolysis. The docked nucleotides do not intrude onto the surface of residues involved in dimerization or chaperoning. Our molecular modeling permitted recognition of larger structural changes in the nucleotide-bound CTD dimer, including stabilization of helix-2 in both chains and intra- and inter- chain interactions between three residues (I613, Q617, R620).     

Computational determination of binding structures and free energies of glucose 6-phosphate dehydrogenase with novel steroid inhibitors

Glucose 6-phosphate dehydrogenase (G6PD), the first and the rate-limiting enzyme in the pentose phosphate pathway (PPP), catalyzes the oxidation of G6P to 6-phosphogluconolactone and the reduction of NADP⁺ to NADPH. Its key role in cancer promotes the development of a potent and selective inhibitor that might increase cancer cell death when combined with radiotherapy. In the present study, we investigated the detailed binding modes and binding free energies for G6PD interacting with a promising series of recently developed inhibitors, i.e., the steroid derivatives, by performing molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations. The docking indicates that the inhibitors occupy the binding sites of both G6P and NADP⁺. The calculated binding free energies on the basis of the MD-simulated enzyme–inhibitor complexes are in good agreement with the experimental activity data for all of the examined inhibitors. The valuable insights into the detailed enzyme–inhibitor binding including the important intermolecular interactions, e.g., the hydrogen bond interaction and the hydrophobic interaction, have been provided. The computational results provide new insights into future rational design of more potent inhibitors of G6PD as a treatment for cancer.     

Comparison of homology models and X-ray structures of the nuclear receptor CAR: assessing the structural basis of constitutive activity

The constitutive androstane receptor (CAR) possesses an intrinsic basal activity whose structural basis has been analysed during the last decade. Recently, we published a homology model of the CAR ligand binding domain (LBD) based on the X-ray structures of the closely related pregnane X (PXR) and vitamin D (VDR) receptor. A detailed analysis of the homology model and molecular dynamics (MD) simulations afforded us to propose a potential mechanism underlying the constitutive activity of CAR. Almost simultaneously, X-ray structures of human and mouse CAR LBD were released. In the present study, a detailed analysis and comparison of homology model and X-ray structures is carried out in order to evaluate the quality and reliability of our homology modelling procedure. The hypothesis of the constitutive activity which we proposed on the basis of our modelling results was tested for consistency with the crystal structures. In addition, the features stated to be essential for the basal activity based on the X-ray data were investigated by means of molecular dynamics simulations. Our results show that the homology modelling procedure was able to predict the CAR LBD structure with high accuracy. Structural features that have been revealed as critical for constitutive activity in the model are also observed in the X-ray structures. Furthermore, the MD simulations of the CAR X-ray structures and a detailed analysis of other NRs clarify the role of distinct structural features that have been assigned an important role for the constitutive activity.     

Designing of multi-targeted molecules using combination of molecular screening and in silico drug cardiotoxicity prediction approaches

We have previously investigated and reported a set of phenol- and indole-based derivatives at the binding pockets of carbonic anhydrase isoenzymes using in silico and in vitro analyses. In this study, we extended our analysis to explore multi-targeted molecules from this set of compounds. Thus, 26 ligands are screened at the binding sites of 229 proteins from 5 main enzyme family classes using molecular docking algorithms. Derived docking scores are compared with reported results of ligands at carbonic anhydrase I and II isoenzymes. Results showed potency of multi-targeted drugs of a few compounds from investigated ligand set. These promising ligands are then tested in silico for their cardiotoxicity risks. Results of this work can be used to improve the desired effects of these compounds by molecular engineering studies. In addition these results may lead to further investigation of studied molecules by medicinal chemists to explore different therapeutic aims.     

Exploring the molecular basis of the enantioselective binding of penicillin G acylase towards a series of 2-aryloxyalkanoic acids: a docking and molecular dynamics study

In the present paper, molecular modeling studies were undertaken in order to shed light on the molecular basis of the observed enantioselectivity of penicillin G acylase (PGA), a well known enzyme for its industrial applications, towards 16 racemic 2-aryloxyalkanoic acids, which have been reported to affect several biological systems. With this intention docking calculations and MD simulations were performed. Docking results indicated that the (S)-enantiomers establish several electrostatic interactions with SerB1, SerB386 and ArgB263 of PGA. Conversely, the absence of specific polar interactions between the (R)-enantiomers and ArgB263 seems to be the main reason for the different binding affinities observed between the two enantiomers. Results of molecular dynamics simulations demonstrated that polar interactions are responsible for both the ligand affinity and PGA enantiospecificity. Modeling calculations provided possible explanations for the observed enantioselectivity of the enzyme that rationalize available experimental data and could be the basis for future protein engineering efforts.     

Identification of non-resistant ROS-1 inhibitors using structure based pharmacophore analysis

Proto-oncogene receptor tyrosine kinase ROS-1 plays a key role in regulating a variety of cancers mainly non-small cell lung cancer (NSCLC). The marketed ROS-1 inhibitors such as Crizotinib suffer from the tribulations of growing resistance due to mutations primarily Gly2032Arg in the ROS-1 protein. To curb the problem of resistance, researchers have developed inhibitors such as Lorlatinib against the mutant protein. The present study was designed to identify inhibitors against wild type (WT) as well as mutant ROS-1 protein that will offer a broader spectrum of activity. Exploring crystal structure of ROS-1 complexed with Lorlatinib, receptor-ligand pharmacophore model was developed using Discovery Studio (DS) software. The developed pharmacophore model consisted of one hydrogen bond acceptor (HBA), one hydrogen bond donor (HBD) and two hydrophobic features (HY), subsequently utilized for virtual screening of commercially available databases and the retrieved hits were further subjected to fitness score and Lipinski’s filter. Thereafter, the retrieved hits were docked in WT and mutated (Gly2032Arg) proteins of ROS-1. Total five molecules were retrieved with good docking scores and good binding interactions within the active site of WT and mutated ROS-1. The binding energies of the ligand-receptor complexes were predicted via calculation of MM-GBSA score. To predict the stability of the ligand receptor complexes with mutant and wild type proteins, molecular dynamic simulation was performed. Thus, these identified hits show good binding affinities with WT and mutant ROS-1 proteins that may be further evaluated for their in-vitro/in-vivo activity.     

Computational docking,molecular dynamics simulation and subsite structure analysis of a maltogenic amylase from Bacillus lehensis G1 provide insights into substrate and product specificity

Maltogenic amylase (MAG1) from Bacillus lehensis G1 displayed the highest hydrolysis activity on β-cyclodextrin (β-CD) to produce maltose as a main product and exhibited high transglycosylation activity on malto-oligosaccharides with polymerization degree of three and above. These substrate and product specificities of MAG1 were elucidated from structural point of view in this study. A three-dimensional structure of MAG1 was constructed using homology modeling. Docking of β-CD and malto-oligosaccharides was then performed in the MAG1 active site. An aromatic platform in the active site was identified which is responsible in substrate recognition especially in determining the enzyme’s preference toward β-CD. Molecular dynamics (MD) simulation showed MAG1 structure is most stable when docked with β-CD and least stable when docked with maltose. The docking analysis and MD simulation showed that the main subsites for substrate stabilization in the active site are −2, −1, +1 and +2. A bulky residue, Trp359 at the +2 subsite was identified to cause steric interference to the bound linear malto-oligosaccharides thus prevented it to occupy subsite +3, which can only be reached by a highly bent glucose molecule such as β-CD. The resulted modes of binding from docking simulation show a good correlation with the experimentally determined hydrolysis pattern. The subsite structure generated from this study led to a possible mode of action that revealed how maltose was mainly produced during hydrolysis. Furthermore, maltose only occupies subsite +1 and +2, therefore could not be hydrolyzed or transglycosylated by the enzyme. This important knowledge has paved the way for a novel structure-based molecular design for modulation of its catalytic activities.