Strong and weak hydrogen bonds in protein-ligand complexes of kinases: a comparative study

Strong and weak hydrogen bonds between protein and ligand are analyzed in a group of 233 X-ray crystal structures of the kinase family. These kinases are from both eukaryotic and prokaryotic organisms. The dataset comprises of 44 sub-families, out of which 35 are of human origin and the rest belong to other organisms. Interaction analysis was carried out in the active sites, defined here as a sphere of 10 A radius around the ligand. A majority of the interactions are observed between the main chain of the protein and the ligand atoms. As a donor, the ligand frequently interacts with amino acid residues like Leu, Glu and His. As an acceptor, the ligand interacts often with Gly, and Leu. Strong hydrogen bonds N-H...O, O-H...O, N-H...N and weak bonds C-H...O, C-H...N are common between the protein and ligand. The hydrogen bond donor capacity of Gly in N-H...O and C-H...O interactions is noteworthy. Similarly, the acceptor capacity of main chain Glu is ubiquitous in several kinase sub-families. Hydrogen bonds between protein and ligand form characteristic hydrogen bond patterns (supramolecular synthons). These synthon patterns are unique to each sub-family. The synthon locations are conserved across sub-families due to a higher percentage of conserved sequences in the active sites. The nature of active site water molecules was studied through a novel classification scheme, based on the extent of exposure of water molecules. Water which is least exposed usually participates in hydrogen bond formation with the ligand. These findings will help structural biologists, crystallographers and medicinal chemists to design better kinase inhibitors.     

In silico mutations and molecular dynamics studies on a winged bean chymotrypsin inhibitor protein

Winged bean chymotrypsin inhibitor (WCI) has an intruding residue Asn14 that plays a crucial role in stabilizing the reactive site loop conformation. This residue is found to be conserved in the Kunitz (STI) family of serine protease inhibitors. To understand the contribution of this scaffolding residue on the stability of the reactive site loop, it was mutated in silico to Gly, Ala, Ser, Thr, Leu and Val and molecular dynamics (MD) simulations were carried out on the mutants. The results of MD simulations reveal the conformational variability and range of motions possible for the reactive site loop of different mutants. The N-terminus side of the scissile bond, which is close to a beta-barrel, is conformationally less variable, while the C-terminus side, which is relatively far from any such secondary structural element, is more variable and needs stability through hydrogen-bonding interactions. The simulated structures of WCI and the mutants were docked in the peptide-binding groove of the cognate enzyme chymotrypsin and the ability to form standard hydrogen-bonding interactions at P3, P1 and P2' residues were compared. The results of the MD simulations coupled with docking studies indicate that hydrophobic residues like Leu and Val at the 14th position are disruptive for the integrity of the reactive site loop, whereas a residue like Thr, which can stabilize the C-terminus side of the scissile bond, can be predicted at this position. However, the size and charge of the Asn residue made it most suitable for the best maintenance of the integrity of the reactive site loop, explaining its conserved nature in the family.     

In silico 3D structure modeling and inhibitor binding studies of human male germ cell-associated kinase

Human male germ cell-associated kinase (hMAK) is an androgen-inducible gene in prostate epithelial cells, and it acts as a coactivator of androgen receptor signaling in prostate cancer. The 3D structure of the hMAK kinase was modeled based on the crystal structure of CDK2 kinase using comparative modeling methods, and the ATP-binding site was characterized. We have collected five inhibitors of hMAK from the literature and docked into the ATP-binding site of the kinase domain. Solvated interaction energies (SIE) of inhibitor binding are calculated from the molecular dynamics simulations trajectories of protein–inhibitor complexes. The contribution from each active site residue in hMAK toward inhibitor binding revealed the nature and extent of interactions between inhibitors and individual residues. The main chain atoms of Met79 invariably form hydrogen bonds with all five inhibitors. The amino acids Leu7, Val15, and Leu129 stabilize the inhibitors via CH–pi interactions. The Asp140 in the active site and Glu77 in hinge region show characteristic hydrogen bonding interactions with inhibitors. From SIE, the residue-wise interactions revealed the nature of non-bonding contacts and modifications required to increase the inhibitor activity. Our work provides 3D model structure of hMAK and molecular basis for the mechanisms of hMAK inhibition at atomic level that aid in designing new potent inhibitors.     

Structure prediction and R115866 binding study of human CYP26A1: homology modelling,fold recognition,molecular docking and MD simulations

Two three-dimensional (3D) models of human cytochrome P450 26A1 (CYP26A1) were constructed using the programs Modeller and Sybyl-GeneFold, respectively. After refinement by molecular mechanics and molecular dynamics (MD) simulations, the two models were validated by structure analysis-validation online server. Subsequently, a flexible docking study was performed on the model constructed by GeneFold with the potent and specific inhibitor R115866 to examine the enzyme–inhibitor interactions. From the docking results, we can see R115866 interacts with amino acid residues at the active site by multiple hydrophobic interactions including the side chains of His111, Trp112, Ser115, Val116, Leu125, Ser126, Leu221, Phe222, Glu296, Phe299, Gly300, Glu303, Thr304, Pro371 and the cofactor heme. Trp112 and Thr304 form hydrogen bonds with R115866 and play important roles in stabilising the complex. This constructed CYP26A1 model may provide an opportunity to understand the action mode of the enzyme and could be useful in designing novel retinoic acid metabolism blocking agents (RAMBAs).     

Insights from the molecular docking of withanolide derivatives to the target protein PknG from Mycobacterium tuberculosis

A crucial virulence factor for intracellular Mycobacterium tuberculosis survival is Protein kinase G (PknG), a eukaryotic-like serinethreonine protein kinase expressed by pathogenic mycobacteria that blocks the intracellular degradation of mycobacteria in lysosomes. Inhibition of PknG results in mycobacterial transfer to lysosomes. Withania somnifera, a reputed herb in ayurvedic medicine, comprises a large number of steroidal lactones known as withanolides which show various pharmacological activities. We describe the docking of 26 withanferin and 14 withanolides from Withania somnifera into the three dimensional structure of PknG of M. tuberculosis using GLIDE. The inhibitor binding positions and affinity were evaluated using scoring functions- Glidescore. The withanolide E, F and D and Withaferin - diacetate 2 phenoxy ethyl carbonate were identified as potential inhibitors of PknG. The available drug molecules and the ligand AX20017 showed hydrogen bond interaction with the aminoacid residues Glu233 and Val235. In addition to Val235 the other amino acids, Gly237, Gln238 and Ser239 are important for withanolide inhibitor recognition via hydrogen bonding mechanisms.     

Crystal Structures of the FAK Kinase in Complex with TAE226 and Related Bis-Anilino Pyrimidine Inhibitors Reveal a Helical DFG Conformation

Focal Adhesion Kinase (FAK) is a non-receptor tyrosine kinase required for cell migration, proliferation and survival. FAK overexpression has been documented in diverse human cancers and is associated with a poor clinical outcome. Recently, a novel bis-anilino pyrimidine inhibitor, TAE226, was reported to efficiently inhibit FAK signaling, arrest tumor growth and invasion and prolong the life of mice with glioma or ovarian tumor implants. Here we describe the crystal structures of the FAK kinase bound to TAE226 and three related bis-anilino pyrimidine compounds. TAE226 induces a conformation of the N-terminal portion of the kinase activation loop that is only observed in FAK, but is distinct from the conformation in both the active and inactive states of the kinase. This conformation appears to require a glycine immediately N-terminal to the “DFG motif”, which adopts a helical conformation stabilized by interactions with TAE226. The presence of a glycine residue in this position contributes to the specificity of TAE226 and related compounds for FAK. Our work highlights the fact that kinases can access conformational space that is not necessarily utilized for their native catalytic regulation, and that such conformations can explain and be exploited for inhibitor specificity.     

Refined X-ray crystal structures of the reactive site modified ovomucoid inhibitor third domains from silver pheasant (OMSVP3*) and from Japanese quail (OMJPQ3*).

Tetragonal and triclinic crystals of two ovomucoid inhibitor third domains from silver pheasant and Japanese quail, modified at their reactive site bonds Met18-Glu19 (OMSVP3*) and Lys18-Asp19 (OMJPQ3*), respectively, were obtained. Their molecular and crystal structures were solved using X-ray data to 2.5 A and 1.55 A by means of Patterson search methods using truncated models of the intact (virgin) inhibitors as search models. Both structures were crystallographically refined to R-values of 0.185 and 0.192, respectively, applying an energy restraint reciprocal space refinement procedure. Both modified inhibitors show large deviations from the intact derivatives only in the proteinase binding loops (Pro14 to Arg21) and in the amino-terminal segments (Leu1 to Val6). In the modified inhibitors the residues immediately adjacent to the cleavage site (in particular P2, P1, P1') are mobile and able to adapt to varying crystal environments. The charged end-groups, i.e. Met18 COO- and Glu19 NH3+ in OMSVP3*, and Lys18 COO- and Asp19 NH3+ in OMJPQ3*, do not form ion pairs with one another. The hydrogen bond connecting the side-chains of Thr17 and Glu19 (i.e. residues on either side of the scissile peptide bond) in OMSVP3 is broken in the modified form, and the hydrogen-bond interactions observed in the intact molecules between the Asn33 side-chain and the carbonyl groups of loop residues P2 and P1' are absent or weak in the modified inhibitors. The reactive site cleavage, however, has little effect on specific interactions within the protein scaffold such as the side-chain hydrogen bond between Asp27 and Tyr31 or the side-chain stacking of Tyr20 and Pro22. The conformational differences in the amino-terminal segment Leu1 to Val6 are explained by their ability to move freely, either to associate with segments of symmetry-related molecules under formation of a four-stranded beta-barrel (OMSVP3* and OMJPQ3) or to bind to surrounding molecules. Together with the results given in the accompanying paper, these findings probably explain why Khyd of small protein inhibitors of serine proteinases is generally found to be so small.     

NMR study of cyclic peptides with renin inhibitor activity

Several cyclic analogues of renin inhibitors, based on Glu-D-Phe-Lys motif have been investigated by NMR spectroscopy and molecular dynamics calculations (MD). The 15 membered macrocycle, resulting from Glu and Lys side-chain cyclization, exhibits conformational preference. The structural evidence from NMR shows the presence of hydrogen bond between Lys NH and Glu side-chain carbonyl, resulting in a 10 membered pseudo beta-turn-like structure. The structure of the cyclic moiety is similar in all the peptides, which takes at least two conformations around Calpha-Cbeta in Glu side chain. The restrained MD calculations further support such observations and show that the macrocycle is fairly rigid, with two conformations about the Glu Calpha-Cbeta bond. The linear peptide appendages, which are essential for activity in cyclic peptides, show an extended structure in the beta-region of Ramchandran plot. These calculations also demonstrate that for the most active peptide, two major conformers each exist about the Calpha-CO bond of the Lys, D-Trp and Leu residues. In this peptide, the cyclic moiety presents a negatively charged surface formed due to the carbonyl oxygens, which are thus available to form hydrogen bonds with the receptor. The linear fragment presents further binding sites with a surface which has the hydrophobic side chains of D-Trp, Leu and D-Met on one side and carbonyls on the other side.     

Insight into the binding mode of a novel LSD1 inhibitor by molecular docking and molecular dynamics simulations

Abstract

Lysine-specific demethylase (LSD1) is an important enzyme for histone lysine methylation. Downregulated LSD1 expression has been linked to cancer proliferation, migration and invasion, indicating that it is an important target for anti-cancer medication. In the present study, the binding modes of a recent reported new series of LSD1 inhibitor were analyzed by using molecular docking and molecular dynamics simulations. A binding mode of these inhibitors was proposed based on the results. According to this binding mode, Thr628 can form two important hydrogen bonds with these inhibitors. Moreover, if the inhibitors can form an additional hydrogen bond with hydroxyl group of Ser289, the potency of the inhibitor can be greatly improved, such as the best inhibitor (compound 12d) in this series. Hydrophobic interactions between the inhibitors and LSD1 are also key contributor here, such as the interaction between the hydrophobic groups (benzene rings) of the inhibitors and the hydrophobic residues of LSD1 (including Val288, Val317, Val811, Ala814, Leu659, Trp751 and Tyr761). Based on the results and analysis, it may provide some useful information for future novel LSD1 inhibitor design.     

Oral administration of FAK inhibitor TAE226 inhibits the progression of peritoneal dissemination of colorectal cancer

Peritoneal dissemination is one of the most terrible types of colorectal cancer progression. Focal adhesion kinase (FAK) plays a crucial role in the biological processes of cancer, such as cell attachment, migration, proliferation and survival, all of which are essential for the progression of peritoneal dissemination. Since we and other groups have reported that the inhibition of FAK activity exhibited a potent anticancer effect in several cancer models, we hypothesized that TAE226, a novel ATP-competitive tyrosine kinase inhibitor designed to target FAK, can prevent the occurrence and progression of peritoneal dissemination. In vitro, TAE226 greatly inhibited the proliferation and migration of HCT116 colon cancer cells, while their adhesion on the matrix surface was minimally inhibited when FAK activity and expression was suppressed by TAE226 and siRNA. In vivo, when HCT116 cells were intraperitoneally inoculated in mice, the cells could attach to the peritoneum and begin to grow within 24 h regardless of the pretreatment of cells with TAE226 or FAK-siRNA, suggesting that FAK is not essential, at least for the initial integrin-matrix contact. Interestingly, the treatment of mice before and after inoculation significantly suppressed cell attachment to the peritoneum. Furthermore, oral administration of TAE226 greatly reduced the size of disseminated tumors and prolonged survival in tumor-bearing mice. Taken together, a possible strategy for inhibiting peritoneal dissemination by targeting FAK with TAE226 appears to be applicable through anti-proliferative and anti-invasion/anti-migration mechanisms.     

Intermolecular insights into allosteric inhibition of histone lysine-specific demethylase 1

BackgroundHistone lysine-specific demethylase 1 (LSD1) has become a potential anticancer target for the novel drug discovery. Recent reports have shown that SP2509 and its derivatives strongly inhibit LSD1 as allosteric inhibitors. However, the binding mechanism of these allosteric inhibitors in the allosteric site of LSD1 is not known yet.MethodsThe stability and binding mechanism of allosteric inhibitors in the binding site of LSD1 were evaluated by molecular docking, ligand-based pharmacophore, molecular dynamics (MD) simulations, molecular mechanics generalized born surface area (MM/GBSA) analysis, quantum mechanics/molecular mechanics (QM/MM) calculation and Hirshfeld surface analysis.ResultsThe conformational geometry and the intermolecular interactions of allosteric inhibitors showed high binding affinity towards allosteric site of LSD1 with the neighboring amino acids (Gly358, Cys360, Leu362, Asp375 and Glu379). Meanwhile, MD simulations and MM/GBSA analysis were performed on selected allosteric inhibitors in complex with LSD1 protein, which confirmed the high stability and binding affinity of these inhibitors in the allosteric site of LSD1.ConclusionThe simulation results revealed the crucial factors accounting for allosteric inhibitors of LSD1, including different protein–ligand interactions, the positions and conformations of key residues, and the ligands flexibilities. Meanwhile, a halogen bond interaction between chlorine atom of ligand and key residues Trp531 and His532 was recurrent in our analysis confirming its importance.General significanceOverall, our research analyzed in depth the binding modes of allosteric inhibitors with LSD1 and could provide useful information for the design of novel allosteric inhibitors.     

Amino acid sequence of a protease inhibitor isolated from Sarcophaga bullata determined by mass spectrometry.

The amino acid sequence of a protease inhibitor isolated from the hemolymph of Sarcophaga bullata larvae was determined by tandem mass spectrometry. Homology considerations with respect to other protease inhibitors with known primary structures assisted in the choice of the procedure followed in the sequence determination and in the alignment of the various peptides obtained from specific chemical cleavage at cysteines and enzyme digests of the S. bullata protease inhibitor. The resulting sequence of 57 residues is as follows: Val Asp Lys Ser Ala Cys Leu Gln Pro Lys Glu Val Gly Pro Cys Arg Lys Ser Asp Phe Val Phe Phe Tyr Asn Ala Asp Thr Lys Ala Cys Glu Glu Phe Leu Tyr Gly Gly Cys Arg Gly Asn Asp Asn Arg Phe Asn Thr Lys Glu Glu Cys Glu Lys Leu Cys Leu.     

Multiple binding modes for palmitate to barley lipid transfer protein facilitated by the presence of proline 12

Molecular dynamics simulations have been used to characterise the binding of the fatty acid ligand palmitate in the barley lipid transfer protein 1 (LTP) internal cavity. Two different palmitate binding modes (1 and 2), with similar protein–ligand interaction energies, have been identified using a variety of simulation strategies. These strategies include applying experimental protein–ligand atom–atom distance restraints during the simulation, or protonating the palmitate ligand, or using the vacuum GROMOS 54B7 force‐field parameter set for the ligand during the initial stages of the simulations. In both the binding modes identified the palmitate carboxylate head group hydrogen bonds with main chain amide groups in helix A, residues 4 to 19, of the protein. In binding mode 1 the hydrogen bonds are to Lys 11, Cys 13, and Leu 14 and in binding mode 2 to Thr 15, Tyr 16, Val 17, Ser 24 and also to the OH of Thr 15. In both cases palmitate binding exploits irregularity of the intrahelical hydrogen‐bonding pattern in helix A of barley LTP due to the presence of Pro 12. Simulations of two variants of barley LTP, namely the single mutant Pro12Val and the double mutant Pro12Val Pro70Val, show that Pro 12 is required for persistent palmitate binding in the LTP cavity. Overall, the work identifies key MD simulation approaches for characterizing the details of protein–ligand interactions in complexes where NMR data provide insufficient restraints.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Computational modeling of human coreceptor CCR5 antagonist as a HIV-1 entry inhibitor: using an integrated homology modeling,docking, and membrane molecular dynamics simulation analysis approach

Chemokine receptor 5 (CCR5) is an integral membrane protein that is utilized during human immunodeficiency virus type-1 entry into host cells. CCR5 is a G-protein coupled receptor that contains seven transmembrane (TM) helices. However, the crystal structure of CCR5 has not been reported. A homology model of CCR5 was developed based on the recently reported CXCR4 structure as template. Automated docking of the most potent (14), medium potent (37), and least potent (25) CCR5 antagonists was performed using the CCR5 model. To characterize the mechanism responsible for the interactions between ligands (14, 25, and 37) and CCR5, membrane molecular dynamic (MD) simulations were performed. The position and orientation of ligands (14, 25, and 37) were found to be changed after MD simulations, which demonstrated the ability of this technique to identify binding modes. Furthermore, at the end of simulation, it was found that residues identified by docking were changed and some new residues were introduced in the proximity of ligands. Our results are in line with the majority of previous mutational reports. These results show that hydrophobicity is the determining factor of CCR5 antagonism. In addition, salt bridging and hydrogen bond contacts between ligands (14, 25, and 37) and CCR5 are also crucial for inhibitory activity. The residues newly identified by MD simulation are Ser160, Phe166, Ser180, His181, and Trp190, and so far no site-directed mutagenesis studies have been reported. To determine the contributions made by these residues, additional mutational studies are suggested. We propose a general binding mode for these derivatives based on the MD simulation results of higher (14), medium (37), and lower (25) potent inhibitors. Interestingly, we found some trend for these inhibitors such as, salt bridge interaction between basic nitrogen of ligand and acidic Glu283 seemed necessary for inhibitory activity. Also, two aromatic pockets (pocket I – TM1-3 and pocket II – TM3-6) were linked by the central polar region in TM7, and the simulated inhibitors show important interactions with the Trp86, Tyr89, Tyr108, Phe112, Ile198, Tyr251, Leu255, and Gln280 and Glu283 residues. These results shed light on the usage of MD simulation to identify more stable, optimal binding modes of the inhibitors.     

Molecular docking studies on quinazoline antifolate derivatives as human thymidylate synthase inhibitors

We have performed molecular docking on quinazoline antifolates complexed with human thymidylate synthase to gain insight into the structural preferences of these inhibitors. The study was conducted on a selected set of one hundred six compounds with variation in structure and activity. The structural analyses indicate that the coordinate bond interactions, the hydrogen bond interactions, the van der Waals interactions as well as the hydrophobic interactions between ligand and receptor are responsible simultaneously for the preference of inhibition and potency. In this study, fast flexible docking simulations were performed on quinazoline antifolates derivatives as human thymidylate synthase inhibitors. The results indicated that the quinazoline ring of the inhibitors forms hydrophobic contacts with Leu192, Leu221 and Tyr258 and stacking interaction is conserved in complex with the inhibitor and cofactor.     

Deciphering the mechanism behind the varied binding activities of COXIBs through Molecular Dynamic Simulations,MM-PBSA binding energy calculations and per-residue energy decomposition studies

COX-2 is a well-known drug target in inflammatory disorders. COX-1/COX-2 selectivity of NSAIDs is crucial in assessing the gastrointestinal side effects associated with COX-1 inhibition. Celecoxib, rofecoxib, and valdecoxib are well-known specific COX-2 inhibiting drugs. Recently, polmacoxib, a COX-2/CA-II dual inhibitor has been approved by the Korean FDA. These COXIBs have similar structure with diverse activity range. Present study focuses on unraveling the mechanism behind the 10-fold difference in the activities of these sulfonamide-containing COXIBs. In order to obtain insights into their binding with COX-2 at molecular level, molecular dynamics simulations studies, and MM-PBSA approaches were employed. Further, per-residue decomposition of these energies led to the identification of crucial amino acids and interactions contributing to the differential binding of COXIBs. The results clearly indicated that Leu338, Ser339, Arg499, Ile503, Phe504, Val509, and Ser516 (Leu352, Ser353, Arg513, Ile517, Phe518, Val523, and Ser530 in PGHS-1 numbering) were imperative in determining the activity of these COXIBs. The binding energies and energy contribution of various residues were similar in all the three simulations. The results suggest that hydrogen bond interaction between the hydroxyl group of Ser516 and five-membered ring of diarylheterocycles augments the affinity in COXIBs. The SAR of the inhibitors studied and the per-residue energy decomposition values suggested the importance of Ser516. Additionally, the positive binding energy obtained with Arg106 explains the binding of COXIBs in hydrophobic channel deep in the COX-2 active site. The findings of the present work would aid in the development of potent COX-2 inhibitors.     

A computational study of binding between 3-(4-fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide and tubulin

3-(4-Fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide (FFSA) is a potential tubulin polymerisation inhibitor. In this article, a theoretical study of the binding between FFSA and tubulin in colchicine site was carried out by molecular docking, molecular dynamics (MD) simulation and binding free energy calculations. The docking calculations preliminarily indicate that there are three possible binding modes 1, 2 and 3; MD simulations and binding free energy calculations identify that binding mode 2 is the most favourable, with the lowest binding free energy of ? 29.54 kcal/mol. Moreover, our valuable results for the binding are as follows: the inhibitor FFSA is suitably located at the colchicine site of tubulin, where it not only interacts with residues Leu248β, Lys254β, Leu255β, Lys352β, Met259β and Val181a by hydrophilic interaction, but also interacts with Val181α and Thr179α by hydrogen bond interaction. These two factors are both essential for FFSA strongly binding to tubulin. These theoretical results help understanding the action mechanism and designing new compounds with higher affinity to tubulin.     

Subsite specificity of memapsin 2 (beta-secretase): implications for inhibitor design

Memapsin 2 is the protease known as beta-secretase whose action on beta-amyloid precursor protein leads to the production of the beta-amyloid (Abeta) peptide. Since the accumulation of Abeta in the brain is a key event in the pathogenesis of Alzheimer's disease, memapsin 2 is an important target for the design of inhibitory drugs. Here we describe the residue preference for the subsites of memapsin 2. The relative k(cat)/K(M) values of residues in each of the eight subsites were determined by the relative initial cleavage rates of substrate mixtures as quantified by MALDI-TOF mass spectrometry. We found that each subsite can accommodate multiple residues. The S(1) subsite is the most stringent, preferring residues in the order of Leu > Phe > Met > Tyr. The preferences of other subsites are the following: S(2), Asp > Asn > Met; S(3), Ile > Val > Leu; S(4), Glu > Gln > Asp; S(1)', Met > Glu > Gln > Ala; S(2)', Val > Ile > Ala; S(3)', Leu > Trp > Ala; S(4)', Asp > Glu > Trp. In general, S subsites are more specific than the S' subsites. A peptide comprising the eight most favored residues (Glu-Ile-Asp-Leu-Met-Val-Leu-Asp) was found to be hydrolyzed with the highest k(cat)/K(M) value so far observed for memapsin 2. Residue preferences at four subsites were also studied by binding of memapsin 2 to a combinatorial inhibitor library. From 10 tight binding inhibitors, the consensus preferences were as follows: S(2), Asp and Glu; S(3), Leu and Ile; S(2)', Val; and S(3)', Glu and Gln. An inhibitor, OM00-3, Glu-Leu-Asp-LeuAla-Val-Glu-Phe (where the asterisk represents the hydroxyethylene tansition-state isostere), designed from the consensus residues, was found to be the most potent inhibitor of memapsin 2 so far reported (K(i) of 3.1 x 10(-10) M). A molecular model of OM00-3 binding to memapsin 2 revealed critical improvement of the interactions between inhibitor side chains with enzyme over a previous inhibitor, OM99-2 [Ghosh, A. K., et al. (2000) J. Am. Chem. Soc. 14, 3522-3523].