Effects of protein flexibility on the site of metabolism prediction for CYP2A6 substrates

Structure-based prediction for the site of metabolism (SOM) of a compound metabolized by human cytochrome P450s (CYPs) is highly beneficial in drug discovery and development. However, the flexibility of the CYPs’ active site remains a huge challenge for accurate SOM prediction. Compared with other CYPs, the active site of CYP2A6 is relatively small and rigid. To address the impact of the flexibility of CYP2A6 active site residues on the SOM prediction for substrates, in this work, molecular dynamics (MD) simulations and molecular docking were used to predict the SOM of 96 CYP2A6 substrates. Substrates with known SOM were docked into the snapshot structures from MD simulations and the crystal structures of CYP2A6. Compared to the crystal structures, the protein structures obtained from MD simulations showed more accurate prediction for SOM. Our results indicated that the flexibility of the active site of CYP2A6 significantly affects the SOM prediction results. Further analysis for the 40 substrates with definite K_m values showed that the prediction accuracy for the low K_m substrates is comparable to that by ligand-based methods.     

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.     

Computational analysis of Asp519 and Glu665 mutations of coagulation factor FVIIIa: Implications for enhanced binding affinity of A2-domain

The A2-domain of blood coagulation factor VIIIa is non-covalently bound to the A1 and A3 domains via weak intermolecular interactions. Functional instability due to rapid dissociation of A2-domain from the active FVIII in blood presents a major hurdle for the therapeutic applications of FVIIIa to treat Hemophilia-A. To identify the ideal hot-spot residues at the interface of A2 and A1/A3 domains that could enhance the structural stability of FVIIIa, we performed a comprehensive computational mutagenesis study of two A2-domain residues, Asp519 and Glu665, that interface the A1 and A3-domains. Each residue was mutated to 15 uncharged amino-acids and the mutant structures were refined by MD simulations. Based on the estimated relative binding affinities of mutant structures, we predict that the mutation of Asp519 to Leu, Gln, Thr, Val and the mutation of Glu665 to Val, Ile, Met, Asn and Trp enhance the A2-domain binding affinity by more than 20 kcal/mol, compared to the WT structure. We anticipate that these predictions will be valuable for enzymatic studies towards the rational design of FVIIIa synthetic constructs with improved A2-domain binding affinity.     

QM/MM investigation of the reaction rates of substrates of 2,3-dimethylmalate lyase: A catabolic protein isolated from Aspergillus niger

Aspergillus niger is an industrially important microorganism used in the production of citric acid. It is a common cause of food spoilage and represents a health issue for patients with compromised immune systems. Recent studies on Aspergillus niger have revealed details on the isocitrate lyase (ICL) superfamily and its role in catabolism, including (2R, 3S)-dimethylmalate lyase (DMML). Members of this and related lyase super families are of considerable interest as potential treatments for bacterial and fungal infections, including Tuberculosis. In our efforts to better understand this class of protein, we investigate the catalytic mechanism of DMML, studying five different substrates and two different active site metals configurations using molecular dynamics (MD) and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. We show that the predicted barriers to reaction for the substrates show good agreement with the experimental k_cat values. This results help to confirm the validity of the proposed mechanism and open up the possibility of developing novel mechanism based inhibitors specifically for this target.     

Molecular docking and molecular dynamics simulation studies of Trypanosoma cruzi triosephosphate isomerase inhibitors. Insights into the inhibition mechanism and selectivity

Trypanosoma cruzi (T. cruzi) triosephosphate isomerase (TcTIM) is a glycolytic enzyme essential for parasite survival and has been considered an interesting target for the development of new antichagasic compounds. The homodimeric enzyme is catalytically active only as a dimer. Interestingly, significant differences exist between the human and parasite TIMs interfaces with a sequence identity of 52%. Therefore, compounds able to specifically disrupt TcTIM but not Homo sapiens TIM (hTIM) dimer interface could become selective antichagasic drugs. In the present work, the binding modes of 1,2,4-thiadiazol, phenazine and 1,2,6-thiadiazine derivatives to TcTIM were investigated using molecular docking combined with molecular dynamics (MD) simulations. The results show that phenazine and 1,2,6-thiadiazine derivatives, 2 and 3, act as dimer-disrupting inhibitors of TcTIM having also allosteric effects in the conformation of the active site. On the other hand, the 1,2,4-thiadiazol derivative 1 binds into the active site causing a significant decrease in enzyme mobility in both monomers. The loss of conformational flexibility upon compound 1 binding suggests that this inhibitor could be preventing essential motions of the enzyme required for optimal activity. The lack of inhibitory activity of 1 against hTIM was also investigated and seems to be related with the high mobility of hTIM which would hinder the formation of a stable ligand–enzyme complex. This work has contributed to understand the mechanism of action of this kind of inhibitors and could result of great help for future rational novel drug design.     

Structural analysis of a novel N-carbamoyl-d-amino acid amidohydrolase from a Brazilian Bradyrhizobium japonicum strain: In silico insights by molecular modelling,docking and molecular dynamics

In this work we performed several in silico analyses to describe the relevant structural aspects of an enzyme N-Carbamoyl-d-amino acid amidohydrolase (d-NCAase) encoded on the genome of the Brazilian strain CPAC 15 (=SEMIA 5079) of Bradyrhizobium japonicum, a nonpathogenic species belonging to the order Rhizobiales. d-NCAase has wide applications particularly in the pharmaceutical industry, since it catalyzes the production of d-amino acids such as D-p-hydroxyphenylglycine (D-HPG), an intermediate in the synthesis of β-lactam antibiotics. We applied a homology modelling approach and 50 ns of molecular dynamics simulations to predict the structure and the intersubunit interactions of this novel d-NCAase. Also, in order to evaluate the substrate binding site, the model was subjected to 50 ns of molecular dynamics simulations in the presence of N-Carbamoyl-d-p-hydroxyphenylglycine (Cp-HPG) (a d-NCAase canonical substrate) and water-protein/water-substrate interactions analyses were performed. Overall, the structural analysis and the molecular dynamics simulations suggest that d-NCAase of B. japonicum CPAC-15 has a homodimeric structure in solution. Here, we also examined the substrate specificity of the catalytic site of our model and the interactions with water molecules into the active binding site were comprehensively discussed. Also, these simulations showed that the amino acids Lys123, His125, Pro127, Cys172, Asp174 and Arg176 are responsible for recognition of ligand in the active binding site through several chemical associations, such as hydrogen bonds and hydrophobic interactions. Our results show a favourable environment for a reaction of hydrolysis that transforms N-Carbamoyl-d-p-hydroxyphenylglycine (Cp-HPG) into the active compound D-p-hydroxyphenylglycine (D-HPG). This work envisage the use of d-NCAase from the Brazilian Bradyrhizobium japonicum strain CPAC-15 (=SEMIA 5079) for the industrial production of D-HPG, an important intermediate for semi-synthesis of β-lactam antibiotics such as penicillins, cephalosporins and amoxicillin.     

Effects of mutations on active site conformation and dynamics of RNA-dependent RNA polymerase from Coxsackievirus B3

Recent crystal structures of RNA-dependent RNA polymerase (3D^pol) from Coxsackievirus B3 (CVB3) revealed that a tyrosine mutation at Phe364 (F364Y) resulted in structures with open active site whereas a hydrophobic mutation at Phe364 (F364A) led to conformations with closed active site. Besides, the crystal structures showed that the F364W mutation had no preference between the open and closed active sites, similar to wild-type. In this paper, we present a molecular dynamics (MD) study on CVB3 3D^pol in order to address some important questions raised by experiments. First, MD simulations of F364Y and F364A were carried out to explore how these mutations at Phe364 influence active site dynamics and conformations. Second, MD simulations of wild-type and mutants were performed to discover the connection between active site dynamics and polymerase function. MD simulations reveal that the effect of mutations on active site dynamics is associated with the interaction between the structural motifs A and D in CVB3 3D^pol. Interestingly, we discover that the active site state is influenced by the formation of a hydrogen bond between backbone atoms of Ala231 (in motif A) and Ala358 (in motif D), which has never been revealed before.     

Molecular modeling and molecular dynamics simulation study of archaeal leucyl-tRNA synthetase in complex with different mischarged tRNA in editing conformation

Aminoacyl-tRNA synthetases (aaRSs) play important roles in maintaining the accuracy of protein synthesis. Some aaRSs accomplish this via editing mechanisms, among which leucyl-tRNA synthetase (LeuRS) edits non-cognate amino acid norvaline mainly by post-transfer editing. However, the molecular basis for this pathway for eukaryotic and archaeal LeuRS remain unclear. In this study, a complex of archaeal P. horikoshii LeuRS (PhLeuRS) with misacylated tRNA^Leu was modeled wherever tRNA’s acceptor stem was oriented directly into the editing site. To understand the distinctive features of organization we reconstructed a complex of PhLeuRS with tRNA and visualize post-transfer editing interactions mode by performing molecular dynamics (MD) simulation studies. To study molecular basis for substrate selectivity by PhLeuRS’s editing site we utilized MD simulation of the entire LeuRS complexes using a diverse charged form of tRNAs, namely norvalyl-tRNA^Leu and isoleucyl-tRNA^Leu. In general, the editing site organization of LeuRS from P.horikoshii has much in common with bacterial LeuRS. The MD simulation results revealed that the post-transfer editing substrate norvalyl-A76, binds more strongly than isoleucyl-A76. Moreover, the branched side chain of isoleucine prevents water molecules from being closer and hence the hydrolysis reaction slows significantly. To investigate a possible mechanism of the post-transfer editing reaction, by PhLeuRS we have determined that two water molecules (the attacking and assisting water molecules) are localized near the carbonyl group of the amino acid to be cleaved off. These water molecules approach the substrate from the opposite side to that observed for Thermus thermophilus LeuRS (TtLeuRS). Based on the results obtained, it was suggested that the post-transfer editing mechanism of PhLeuRS differs from that of prokaryotic TtLeuRS.     

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.     

Effects of different solvents on the conformations of apoptotic cytochrome c: Structural insights from molecular dynamics simulation

Cytochrome c (cyt-c) upon binding with cardiolipin acquires peroxidase activity and is strictly connected to the pathogenesis of many human diseases including neurodegenerative and cardiovascular diseases. Interaction of cyt-c with cardiolipin mimics partial unfolding/conformational changes of cyt-c in different solvent environments. Dynamic pictures of these conformational changes of cyt-c are crucial in understanding their physiological roles in mitochondrial functions. Therefore, atomistic molecular dynamics (MD) simulations have been carried out to investigate the effect of different solvents (water, urea/water, MeOH and DMSO) on the structure and conformations of apoptotic cyt-c (Fe³⁺). Our study demonstrates that the structural changes in the protein are solvent dependent. The structural differences are observed majorly on the β-sheets and α-helical conformations and the degree of their perturbation are specific to the solvent. Although a complete loss of β-sheets (0%) is observed in MeOH and DMSO, by contrast, well preserved β-sheets (3.84%) are observed in water and urea/water. A significant decrease in the α-helical contents is observed in MeOH (41.34%) and water (42.46%), a negligible alteration in DMSO (44.25%) and well preserved α-helical (45.19%) contents in urea/water. The distances between the residues critical for electron transfer are decreased considerably for DMSO. Further, the reduction in residue flexibility and the conformational space indicate that the collective motions of cyt-c are reduced when compared to other cosolvents. Essential dynamics analysis implies that the overall motions of cyt-c in water, MeOH and urea/water are involved in three to four eigenvectors and in first eigenvector in DMSO. Overall, we believe that MD simulations of cyt-c in different solvents can provide a detailed microscopic understanding of the physiological roles, electron transport and peroxidase function in the early events of apoptosis which are hard to probe experiments.     

Flap opening mechanism of HIV-1 protease

The active site of aspartic proteases, such as HIV-1 protease (PR), is covered by one or more flaps, which restrict access to the active site. For HIV-1 PR, X-ray diffraction studies suggested that in the free enzyme the two flaps are packed onto each other loosely in a semi-open conformation, while molecular dynamics (MD) studies observed that the flaps can also separate into open conformations. In this study, the mechanism of flap opening and the structure and dynamics of HIV-1 PR with semi-open and open flap conformations were investigated using molecular dynamics simulations. The flaps showed complex dynamic behavior as two distinct mechanisms of flap opening and various stable flap conformations (semi-open, open and curled) were observed during the simulations. A network of weakly polar interactions between the flaps were proposed to be responsible for stabilizing the semi-open flap conformation. It is hypothesized that such interactions could be responsible for making flap opening a highly sensitive gating mechanism which control access to the active site.     

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.     

Molecular basis for benzimidazole resistance from a novel β-tubulin binding site model

Benzimidazole-2-carbamate derivatives (BzCs) are the most commonly used antiparasitic drugs for the treatment of protozoan and helminthic infections. BzCs inhibit the microtubule polymerization mechanism through binding selectively to the β-tubulin subunit in which mutations have been identified that lead to drug resistance. Currently, the lack of crystallographic structures of β-tubulin in parasites has limited the study of the binding site and the analysis of the resistance to BzCs. Recently, our research group has proposed a model to explain the interaction between the BzCs and a binding site in the β-tubulin. Herein, we report the homology models of two susceptible (Haemonchus contortus and Giardia intestinalis) parasites and one unsusceptible (Entamoeba histolytica) generated using the structure of the mammal Ovis aries, considered as a low susceptible organism, as a template. Additionally, the mechanism by which the principal single point mutations Phe167Tyr, Glu198Ala and Phe200Tyr could lead to resistance to BzCs is analyzed. Molecular docking and molecular dynamics studies were carried out in order to evaluate the models and the ligand–protein complexes’ behaviors. This study represents a first attempt towards understanding, at the molecular level, the structural composition of β-tubulin in all organisms, also suggesting possible resistance mechanisms. Furthermore, these results support the importance of benzimidazole derivative optimization in order to design more potent and selective (less toxic) molecules for the treatment of parasitic diseases.     

Molecular dynamics simulations of mutated Mycobacterium tuberculosis l-alanine dehydrogenase to illuminate the role of key residues

l-Alanine dehydrogenase from Mycobacterium tuberculosis (l-MtAlaDH) catalyzes the NADH-dependent interconversion of l-alanine and pyruvate, and it is considered to be a potential target for the treatment of tuberculosis. The experiment has verified that amino acid replacement of the conserved active-site residues which have strong stability and no great changes in biological evolutionary process, such as His96 and Asp270, could lead to inactive mutants [Ågren et al., J. Mol. Biol. 377 (2008) 1161–1173]. However, the role of these conserved residues in catalytic reaction still remains unclear. Based on the crystal structures, a series of mutant structures were constructed to investigate the role of the conserved residues in enzymatic reaction by using molecular dynamics simulations. The results show that whatever the conserved residues were mutated, the protein can still convert its conformation from open state to closed state as long as NADH is present in active site. Asp270 maintains the stability of nicotinamide ring and ribose of NADH through hydrogen bond interactions, and His96 is helpful to convert the protein conformation by interactions with Gln271, whereas, they would lead to the structural rearrangement in active site and lose the catalytic activity when they were mutated. Additionally, we deduce that Met301 plays a major role in catalytic reaction due to fixing the nicotinamide ring of NADH to prevent its rotation, and we propose that Met301 would be mutated to the hydrophobic residue with large steric hindrance in side chain to test the activity of the protein in future experiment.     

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.     

Interactions of acetylcholine binding site residues contributing to nicotinic acetylcholine receptor gating: Role of residues Y93, Y190, K145 and D200

The nicotinic acetylcholine receptor exhibits multiple conformational states, resting (channel closed), active (channel open) and desensitized (channel closed). The resting state may be distinguished from the active and desensitized states by the orientation of loop C in the extracellular ligand binding domain (LBD). Homology modeling was used to generate structures of the Torpedo californica α₂βδγ nAChR that initially represent the resting state (loop C open) and the desensitized state (loop C closed). Molecular dynamics (MD) simulations were performed on the extracellular LBD on each nAChR conformational state, with and without the agonist anabaseine present in each binding site (the αγ and the αδ sites). Three MD simulations of 10 ns each were performed for each of the four conditions. Comparison of dynamics revealed that in the presence of agonist, loop C was drawn inward and attains a more stable conformation. Examination of side-chain interactions revealed that residue αY190 exhibited hydrogen-bonding interactions either with residue αY93 in the ligand binding site or with residue αK145 proximal to the binding site. αK145 also exhibited side chain (salt bridge) interactions with αD200 and main chain interactions with αY93. Residues αW149, αY198, γY116/δT119, γL118/δL121 and γL108/δL111 appear to play the role of stabilizing ligand in the binding site. In MD simulations for the desensitized state, the effect of ligand upon the interactions among αK145, αY190, and αY93 as well as ligand-hydrogen-bonding to αW149 were more pronounced at the αγ interface than at the αδ interface. Differences in affinity for the desensitized state were determined experimentally to be 10-fold. The changes in side chain interactions observed for the two conformations and induced by ligand support a model wherein hydrogen bond interactions between αD200 and αY93 are broken and rearrange to form a salt-bridge between αK145 and αD200 and hydrogen bond interactions between αY93 and αY190 and between αK145 and αY190.     

Molecular dynamics simulations of isoleucine-release pathway in GAF domain of N-CodY from Bacillus Subtilis

The GAF domain located in the N-terminal motifs of CodY (N-CodY) is responsible for increasing the affinity of CodY to its target sites on DNA by its interaction with the branched chain amino acids (BCAAs) involving isoleucine, leucine and valine. The study of the interaction of GAF domain with isoleucine gains much attention in recent years, but the mechanism of isoleucine release still remains unclear. In this paper, a conventional molecular dynamics (MD) and force probe molecular dynamics (FPMD) simulations have been performed with the aim to understand how the isoleucine ligand escapes from the GAF domain of N-CodY from Bacillus subtilis. The MD results reveal that the ligand release is a gradual process, which is accompanied by the movement of the loop between β3 and β4. During the periods of ligand escaping from the bottom to the top of binding pocket, isoleucine forms hydrogen bonds one after another with series of residues, such as ARG61, THR96, PHE98, VAL100, GLU101 and ASN102, under the mediation of hydrophobic contacts. The FPMD results show that the easiest way to pull ligand out of the cavity is along x direction (i.e. the direction is opposite to MET62).     

Nanosecond molecular dynamics simulations of Cdc25B and its complex with a 1,4-naphthoquinone inhibitor: implications for rational inhibitor design

Cdc25 phosphatases have been considered as attractive drug targets for anticancer therapies due to the correlation of their overexpression with a wide variety of cancers. To gain insight into designing new potent inhibitors, we investigate the dynamic properties of Cdc25B and its complex with a 1,4-naphtoquinone inhibitor NSC 95397 by means of molecular dynamics simulations in aqueous solution. It is shown from the calculated dynamic properties that the malleability of the residues 530–532 residing at the start of C-terminal region around the active site should be responsible for the catalytic action of Cdc25B. However, binding of the inhibitor in the active site leads to a substantial decrease in the motional amplitude of the flexible residues, due to the hydrophobic interactions with the side chain of Met531. The simulation results also indicate that at least four hydrogen bonds are involved in the enzyme-inhibitor complex. Among them, the hydrogen bond between the side chain carboxylate group of Glu478 and one of the hydroxyl groups of the inhibitor is found to be the most significant binding force stabilizing the inhibitor in the active site. This result supports the previous experimental implication that the possession of a single hydroxyl group is sufficient for the inhibitory activity of 1,4-naphthoquinone inhibitors.     

Binding and discerning interactions of PTP1B allosteric inhibitors: Novel insights from molecular dynamics simulations

The α7 helix is either disordered or missing in the three co-crystal structures of allosteric inhibitors with protein tyrosine phosphatase 1B (PTP1B). It was modeled in each complex using the open form of PTP1B structure and studied using molecular dynamics (MD) simulations for 25 ns. B-factor analysis of the residues sheds light on its disordered nature in the co-crystal structures. Further, the ability of inhibitors to act as allosteric inhibitor was studied and established using novel hydrogen bond criteria. The MD simulations were utilized to determine the relative importance of electrostatic and hydrophobic component in to the binding of inhibitors. It was revealed that the hydrophobic interactions predominantly drive the molecular recognition of these inhibitors. Per residue energy decomposition analysis attributed dissimilar affinities of three inhibitors to the several hydrogen bonds and non-bonded interactions. Among the secondary structure elements that surround the allosteric site, helices α6, α7 and loop α6–α7 were notorious in providing variable affinities to the inhibitors. A novel hydrophobic pocket lined by the α7 helix residues Val287, Asn289 and Trp291 was identified in the allosteric site. This study provides useful insights for the rational design of high affinity PTP1B allosteric inhibitors.     

Study of the interaction of Huperzia saururus Lycopodium alkaloids with the acetylcholinesterase enzyme

In the present study, we describe and compare the binding modes of three Lycopodium alkaloids (sauroine, 6-hydroxylycopodine and sauroxine; isolated from Huperzia saururus) and huperzine A with the enzyme acetylcholinesterase. Refinement and rescoring of the docking poses (obtained with different programs) with an all atom force field helped to improve the quality of the protein–ligand complexes. Molecular dynamics simulations were performed to investigate the complexes and the alkaloid's binding modes. The combination of the latter two methodologies indicated that binding in the active site is favored for the active compounds. On the other hand, similar binding energies in both the active and the peripheral sites were obtained for sauroine, thus explaining its experimentally determined lack of activity. MM-GBSA predicted the order of binding energies in agreement with the experimental IC₅₀ values.