Successful molecular dynamics simulation of the zinc-bound farnesyltransferase using the cationic dummy atom approach

Farnesyltransferase (FT) inhibitors can suppress tumor cell proliferation without substantially interfering with normal cell growth, thus holding promise for cancer treatment. A structure-based approach to the design of improved FT inhibitors relies on knowledge of the conformational flexibility of the zinc-containing active site of FT. Although several X-ray structures of FT have been reported, detailed information regarding the active site conformational flexibility of the enzyme is still not available. Molecular dynamics (MD) simulations of FT can offer the requisite information, but have not been applied due to a lack of effective methods for simulating the four-ligand coordination of zinc in proteins. Here, we report in detail the problems that occurred in the conventional MD simulations of the zinc-bound FT and a solution to these problems by employing a simple method that uses cationic dummy atoms to impose orientational requirement for zinc ligands. A successful 1.0 ns (1.0 fs time step) MD simulation of zinc-bound FT suggests that nine conserved residues (Asn127alpha, Gln162alpha, Asn165alpha, Gln195alpha, His248beta, Lys294beta, Leu295beta, Lys353beta, and Ser357beta) in the active site of mammalian FT are relatively mobile. Some of these residues might be involved in the ligand-induced active site conformational rearrangement upon binding and deserve attention in screening and design of improved FT inhibitors for cancer chemotherapy.     

Structure of malonic acid-based inhibitors bound to human neutrophil collagenase. A new binding mode explains apparently anomalous data.

Matrix metalloproteinases (MMPs) are a family of zinc endopeptidases, which have been implicated in various disease processes. Various classes of MMP inhibitors, including hydroxamic acids, phosphinic acids, and thiols, have been previously described. Most of these mimic peptides, and most likely bind analogous to the corresponding peptide substrates. Among the hydroxamic acids, malonic acid derivatives have been used as MMP inhibitors, although optimization of their inhibition potency was not successful. Here we report the design of malonic acid-based inhibitors using the X-ray structure of a collagenase/inhibitor complex, which revealed a nonsubstrate-like binding mode. The proposed beta-type turn-like conformation for the improved inhibitors was confirmed by X-ray crystallography. The observation of nonsubstrate-like binding confirms the original strategy for structure-based modeling of improved malonic acid inhibitors, and explains kinetic data that are inconsistent with substrate-like binding. Detailed interactions for the improved inhibitors seen in the crystal structure also suggest possibilities for further modifications in cycles of structure based drug design. Indeed, we have designed nonpeptidic inhibitors with approximately 500-fold improved inhibition based on these structures.     

X-ray structure of human stromelysin catalytic domain complexed with nonpeptide inhibitors: implications for inhibitor selectivity.

Effective inhibitors of matrix metalloproteinases (MMPs), a family of connective tissue-degrading enzymes, could be useful for the treatment of diseases such as cancer, multiple sclerosis, and arthritis. Many of the known MMP inhibitors are derived from peptide substrates, with high potency in vitro but little selectivity among MMPs and poor bioavailability. We have discovered nonpeptidic MMP inhibitors with improved properties, and report here the crystal structures of human stromelysin-1 catalytic domain (SCD) complexed with four of these inhibitors. The structures were determined and refined at resolutions ranging from 1.64 to 2.0 A. Each inhibitor binds in the active site of SCD such that a bulky diphenyl piperidine moiety penetrates a deep, predominantly hydrophobic S'1 pocket. The active site structure of the SCD is similar in all four inhibitor complexes, but differs substantially from the peptide hydroxamate complex, which has a smaller side chain bound in the S'1 pocket. The largest differences occur in the loop forming the "top" of this pocket. The occupation of these nonpeptidic inhibitors in the S'1 pocket provides a structural basis to explain their selectivity among MMPs. An analysis of the unique binding mode predicts structural modifications to design improved MMP inhibitors.     

A molecular dynamics study of the interaction of D-peptide amyloid inhibitors with their target sequence reveals a potential inhibitory pharmacophore conformation

The self-assembly of soluble proteins and peptides into β-sheet-rich oligomeric structures and insoluble fibrils is a hallmark of a large number of human diseases known as amyloid diseases. Drugs that are able to interfere with these processes may be able to prevent and/or cure these diseases. Experimental difficulties in the characterization of the intermediates involved in the amyloid formation process have seriously hampered the application of rational drug design approaches to the inhibition of amyloid formation and growth. Recently, short model peptide systems have proved useful in understanding the relationship between amino acid sequence and amyloid formation using both experimental and theoretical approaches. Moreover, short d-peptide sequences have been shown to specifically interfere with those short amyloid stretches in proteins, blocking oligomer formation or disassembling mature fibrils. With the aim of rationalizing which interactions drive the binding of inhibitors to nascent β-sheet oligomers, in this study, we have carried out extensive molecular dynamics simulations of the interaction of selected d-peptide sequences with oligomers of the target model sequence STVIIE. Structural analysis of the simulations helped to identify the molecular determinants of an inhibitory core whose conformational and physicochemical properties are actually shared by nonpeptidic small-molecule inhibitors of amyloidogenesis. Selection of one of these small molecules and experimental validation against our model system proved that it was indeed an effective inhibitor of fibril formation by the STVIIE sequence, supporting theoretical predictions. We propose that the inhibitory determinants derived from this work be used as structural templates in the development of pharmacophore models for the identification of novel nonpeptidic inhibitors of aggregation.     

Computational studies of the farnesyltransferase ternary complex part I: substrate binding

Farnesyltransferase (FTase) catalyzes the transfer of farnesyl from farnesyl diphosphate (FPP) to cysteine residues at or near the C-terminus of protein acceptors with a CaaX motif (a, aliphatic; X, Met). Farnesylation is a critical modification to many switch proteins involved in the extracellular signal transduction pathway, which facilitates their fixation on the cell membrane where the extracellular signal is processed. The malfunction caused by mutations in these proteins often causes uncontrolled cell reproduction and leads to tumor formation. FTase inhibitors have been extensively studied as potential anticancer agents in recent years, several of which have advanced to different phases of clinical trials. However, the mechanism of this biologically important enzyme has not been firmly established. Understanding how FTase recruits the FPP substrate is the first and foremost step toward further mechanistic investigations and the design of effective FTase inhibitors. Molecular dynamic simulations were carried out on the ternary structure of FTase complexed with the FPP substrate and an acetyl-capped tetrapeptide (acetyl-CVIM), which revealed that the FPP substrate maintains an inactive conformation and the binding of the diphosphate group can be largely attributed to residues R291beta, K164alpha, K294beta, and H248beta. The FPP substrate assumes an extended conformation in the binding site with restricted rotation of the backbone dihedral angles; however, it does not have a well-defined conformation when unbound in solution. This is evident from multinanosecond MD simulations of the FPP substrate in a vacuum and solution. Our conclusion is further supported by theoretical J coupling calculations. Our results on the FPP binding are in good agreement with previous experimental kinetic studies on FTase mutants. The hypothetical conformational activation of the FPP substrate is currently under investigation.     

A novel bisphosphonate inhibitor of squalene synthase combined with a statin or a nitrogenous bisphosphonate in vitro

Statins and nitrogenous bisphosphonates (NBP) inhibit 3-hydroxy-3-methylglutaryl-coenzyme-A reductase (HMGCR) and farnesyl diphosphate synthase (FDPS), respectively, leading to depletion of farnesyl diphosphate (FPP) and disruption of protein prenylation. Squalene synthase (SQS) utilizes FPP in the first committed step from the mevalonate pathway toward cholesterol biosynthesis. Herein, we have identified novel bisphosphonates as potent and specific inhibitors of SQS, including the tetrasodium salt of 9-biphenyl-4,8-dimethyl-nona-3,7-dienyl-1,1-bisphosphonic acid (compound 5). Compound 5 reduced cholesterol biosynthesis and lead to a substantial intracellular accumulation of FPP without reducing cell viability in HepG2 cells. At high concentrations, lovastatin and zoledronate impaired protein prenylation and decreased cell viability, which limits their potential use for cholesterol depletion. When combined with lovastatin, compound 5 prevented lovastatin-induced FPP depletion and impairment of protein farnesylation. Compound 5 in combination with the NBP zoledronate completely prevented zoledronate-induced impairment of both protein farnesylation and geranylgeranylation. Cotreatment of cells with compound 5 and either lovastatin or zoledronate was able to significantly prevent the reduction of cell viability caused by lovastatin or zoledronate alone. The combination of an SQS inhibitor with an HMGCR or FDPS inhibitor provides a rational approach for reducing cholesterol synthesis while preventing nonsterol isoprenoid depletion.     

HIV-1 protease inhibitors: enthalpic versus entropic optimization of the binding affinity

Existing experimental as well as computational screening methods select potential ligands or drug candidates on the basis of binding affinity. Since the binding affinity is a function of the enthalpy (DeltaH) and entropy (DeltaS) changes, it is apparent that improved binding can be achieved in different ways: by optimizing DeltaH, DeltaS, or a combination of both. However, the behavior of enthalpically or entropically optimized inhibitors is fundamentally different, including their response to mutations that may elicit drug resistance. In the design of HIV-1 protease inhibitors, high binding affinity has usually been achieved by preshaping lead compounds to the geometry of the binding site and by incorporating a high degree of hydrophobicity. The thermodynamic consequence of that approach is that the binding affinity of the resulting inhibitors becomes entropically favorable but enthalpically unfavorable. Specifically, the resulting high binding affinity is due to an increased solvation entropy (hydrophobic effect) combined with a reduced loss of conformational entropy of the inhibitor upon binding (structural rigidity). Here we report that tripeptide inhibitors derived from the transframe region of Gag-Pol (Glu-Asp-Leu and Glu-Asp-Phe) bind to the HIV-1 protease with a favorable enthalpy change. This behavior is qualitatively different from that of known inhibitors and points to new strategies for inhibitor design. Since the binding affinities of enthalpically favorable and enthalpically unfavorable inhibitors have opposite temperature dependence, it is possible to design fast screening protocols that simultaneously select inhibitors on the basis of affinity and enthalpy.     

Predicting inactive conformations of protein kinases using active structures: conformational selection of type-II inhibitors

Protein kinases have been found to possess two characteristic conformations in their activation-loops: the active DFG-in conformation and the inactive DFG-out conformation. Recently, it has been very interesting to develop type-II inhibitors which target the DFG-out conformation and are more specific than the type-I inhibitors binding to the active DFG-in conformation. However, solving crystal structures of kinases with the DFG-out conformation remains a challenge, and this seriously hampers the application of the structure-based approaches in development of novel type-II inhibitors. To overcome this limitation, here we present a computational approach for predicting the DFG-out inactive conformation using the DFG-in active structures, and develop related conformational selection protocols for the uses of the predicted DFG-out models in the binding pose prediction and virtual screening of type-II ligands. With the DFG-out models, we predicted the binding poses for known type-II inhibitors, and the results were found in good agreement with the X-ray crystal structures. We also tested the abilities of the DFG-out models to recognize their specific type-II inhibitors by screening a database of small molecules. The AUC (area under curve) results indicated that the predicted DFG-out models were selective toward their specific type-II inhibitors. Therefore, the computational approach and protocols presented in this study are very promising for the structure-based design and screening of novel type-II kinase inhibitors.     

A Combination of Receptor-Based Pharmacophore Modeling & QM Techniques for Identification of Human Chymase Inhibitors

Inhibition of chymase is likely to divulge therapeutic ways for the treatment of cardiovascular diseases, and fibrotic disorders. To find novel and potent chymase inhibitors and to provide a new idea for drug design, we used both ligand-based and structure-based methods to perform the virtual screening(VS) of commercially available databases. Different pharmacophore models generated from various crystal structures of enzyme may depict diverse inhibitor binding modes. Therefore, multiple pharmacophore-based approach is applied in this study. X-ray crystallographic data of chymase in complex with different inhibitors were used to generate four structure–based pharmacophore models. One ligand–based pharmacophore model was also developed from experimentally known inhibitors. After successful validation, all pharmacophore models were employed in database screening to retrieve hits with novel chemical scaffolds. Drug-like hit compounds were subjected to molecular docking using GOLD and AutoDock. Finally four structurally diverse compounds with high GOLD score and binding affinity for several crystal structures of chymase were selected as final hits. Identification of final hits by three different pharmacophore models necessitates the use of multiple pharmacophore-based approach in VS process. Quantum mechanical calculation is also conducted for analysis of electrostatic characteristics of compounds which illustrates their significant role in driving the inhibitor to adopt a suitable bioactive conformation oriented in the active site of enzyme. In general, this study is used as example to illustrate how multiple pharmacophore approach can be useful in identifying structurally diverse hits which may bind to all possible bioactive conformations available in the active site of enzyme. The strategy used in the current study could be appropriate to design drugs for other enzymes as well.     

Bioactive conformation of a potent stromelysin inhibitor determined by X-nucleus filtered and multidimensional NMR spectroscopy

The biologically active conformation of a novel, very potent, nonpeptidic stromelysin inhibitor was determined by X-nucleus filtered and multidimensional NMR spectroscopy. This bound conformer was subsequently docked into the stromelysin catalytic domain (SCD) using intermolecular distance constraints derived from NOE data. The complex showed the S1′ pocket of stromelysin to be the major site of enzyme-inhibitor interaction with other portions of the inhibitor spanning the S2′ and S1 binding sites. Theoretical predictions of SCD-inhibitor binding from molecular modeling studies were consistent with the NMR data. Comparison of modeled enzyme-inhibitor complexes for stromelysin and collagenase revealed an alternate binding mode for the inhibitor in collagenase, suggesting a similar binding interaction might also be possible for stromelysin. The NMR results, however, revealed a single SCD-inhibitor binding mode and provided a structural template for the design of more potent stromelysin inhibitors.     

Allosteric model of maraviroc binding to CC chemokine receptor 5 (CCR5)

Maraviroc is a nonpeptidic small molecule human immunodeficiency virus type 1 (HIV-1) entry inhibitor that has just entered the therapeutic arsenal for the treatment of patients. We recently demonstrated that maraviroc binding to the HIV-1 coreceptor, CC chemokine receptor 5 (CCR5), prevents it from binding the chemokine CCL3 and the viral envelope glycoprotein gp120 by an allosteric mechanism. However, incomplete knowledge of ligand-binding sites and the lack of CCR5 crystal structures have hampered an in-depth molecular understanding of how the inhibitor works. Here, we addressed these issues by combining site-directed mutagenesis (SDM) with homology modeling and docking. Six crystal structures of G-protein-coupled receptors were compared for their suitability for CCR5 modeling. All CCR5 models had equally good geometry, but that built from the recently reported dimeric structure of the other HIV-1 coreceptor CXCR4 bound to the peptide CVX15 (Protein Data Bank code 3OE0) best agreed with the SDM data and discriminated CCR5 from non-CCR5 binders in a virtual screening approach. SDM and automated docking predicted that maraviroc inserts deeply in CCR5 transmembrane cavity where it can occupy three different binding sites, whereas CCL3 and gp120 lie on distinct yet overlapped regions of the CCR5 extracellular loop 2. Data suggesting that the transmembrane cavity remains accessible for maraviroc in CCL3-bound and gp120-bound CCR5 help explain our previous observation that the inhibitor enhances dissociation of preformed ligand-CCR5 complexes. Finally, we identified residues in the predicted CCR5 dimer interface that are mandatory for gp120 binding, suggesting that receptor dimerization might represent a target for new CCR5 entry inhibitors.     

Azole-based inhibitors of AKT/PKB for the treatment of cancer

Through a combination of screening and structure-based rational design, we have discovered a series of N¹-(5-(heterocyclyl)-thiazol-2-yl)-3-(4-trifluoromethylphenyl)-1,2-propanediamines that were developed into potent ATP competitive inhibitors of AKT. Studies of linker strand-binding adenine isosteres identified SAR trends in potency and selectivity that were consistent with binding interactions observed in structures of the inhibitors bound to AKT1 and to the counter-screening target PKA. One compound was shown to have acceptable pharmacokinetic properties and to be a potent inhibitor of AKT signaling and of in vivo xenograft tumor growth in a preclinical model of glioblastoma.     

Crystal structure of avian aminoimidazole-4-carboxamide ribonucleotide transformylase in complex with a novel non-folate inhibitor identified by virtual ligand screening

Aminoimidazole-4-carboxamide ribonucleotide transformylase (AICAR Tfase), one of the two folate-dependent enzymes in the de novo purine biosynthesis pathway, is a promising target for anti-neoplastic chemotherapy. Although classic antifolates, such as methotrexate, have been developed as anticancer agents, their general toxicity and drug resistance are major issues associated with their clinical use and future development. Identification of inhibitors with novel scaffolds could be an attractive alternative. We present here the crystal structure of avian AICAR Tfase complexed with the first non-folate based inhibitor identified through virtual ligand screening of the National Cancer Institute Diversity Set. The inhibitor 326203-A (2-[5-hydroxy-3-methyl-1-(2-methyl-4-sulfophenyl)-1H-pyrazol-4-ylazo]-4-sulfo-benzoic acid) displayed competitive inhibition against the natural cofactor, 10-formyl-tetrahydrofolate, with a K(i) of 7.1 mum. The crystal structure of AICAR Tfase with 326203-A at 1.8 A resolution revealed a unique binding mode compared with antifolate inhibitors. The inhibitor also accessed an additional binding pocket that is not occupied by antifolates. The sulfonate group of 326203-A appears to form the dominant interaction of the inhibitor with the proposed oxyanion hole through interaction with a helix dipole and Lys(267). An aromatic interaction with Phe(316) also likely contributes to favorable binding. Based on these structural insights, several inhibitors with improved potency were subsequently identified in the National Cancer Institute Compound Library and the Available Chemical Directory by similarity search and molecular modeling methods. These results provide further support for our combined virtual ligand screening rational design approach for the discovery of novel, non-folate-based inhibitors of AICAR Tfase.     

Inhibitors of the Yersinia protein tyrosine phosphatase through high throughput and virtual screening approaches

The bacterial protein tyrosine phosphatase YopH is an essential virulence determinant in Yersinia pestis and a potential antibacterial drug target. Here we report our studies of screening for small molecule inhibitors of YopH using both high throughput and in silico approaches. The identified inhibitors represent a diversity of chemotypes and novel pTyr mimetics, providing a starting point for further development and fragment-based design of multi-site binding inhibitors. We demonstrate that the applications of high throughput and virtual screening, when guided by structural binding mode analysis, is an effective approach for identifying potent and selective inhibitors of YopH and other protein phosphatases for rational drug design.     

Discovery and design of benzimidazolone based inhibitors of p38 MAP kinase

A new class of benzimidazolone p38 MAP kinase inhibitors was discovered through high-throughput screening. X-ray crystallographic data of the lead molecule with p38 were used to design analogues with improved binding affinity and potency in a cell assay of LPS-induced TNF production. Herein, we report the SAR of this new class of p38 inhibitors.     

Structure modeling, ligand binding, and binding affinity calculation (LR-MM-PBSA) of human heparanase for inhibition and drug design

Human heparanase is an endo-beta-D-glycosidase that cleaves heparan sulphate (HS) chains in the extracellular matrix and basement membrane. It is known that the cleavage of HS by heparanase results in cell invasion and metastasis of cancer. Therefore, heparanase is considered an important target for cancer drug development. The three-dimensional structure of heparanase would be useful in the rational design of inhibitors targeted to the enzyme; however, the three-dimensional structure has not yet been determined. In our effort to design inhibitors, we developed a three-dimensional structure of heparanase using a homology-modeling approach. The homology-built structure is consistent to previous bioinformatics and site-mutation experimental results. The heparanase features a (alpha/beta)(8) TIM-barrel fold with two glutamate residues (Glu225 and Glu343) located in the active-site cleft. This feature supports the putative mechanism of proton donor and nucleophilic sites. Docking simulations yielded 41 complex structures, which indicate that the bound inhibitor could block ligand binding into the catalytic site. A free energy of binding model was established for 25 heparanase inhibitors with a training set of 25 heparanase inhibitors using the linear response MM-PBSA approach (LR-MM-PBSA). The correlation between calculated and experimental activity was 0.79 and the reliability of the model was validated with leave-one-out cross-validation method. Its predictive capability was further validated using a test set of 16 inhibitors similar to the training set of inhibitors. The correlation between the predicted and observed activities is significantly improved by the protein "induced-fit" that accounts for the flexibility of the receptor. These interaction and pharmacophore elements provide a unique insight to the rational design of new ligands targeted to the enzyme.     

Structure and orientation of peptide inhibitors bound to beta-amyloid fibrils

Polymerization of the soluble beta-amyloid peptide into highly ordered fibrils is hypothesized to be a causative event in the development of Alzheimer's disease. Understanding the interactions of Abeta with inhibitors on an atomic level is fundamental for the development of diagnostics and therapeutic approaches, and can provide, in addition, important indirect information of the amyloid fibril structure. We have shown recently that trRDCs can be measured in solution state NMR for peptide ligands binding weakly to amyloid fibrils. We present here the structures for two inhibitor peptides, LPFFD and DPFFL, and their structural models bound to fibrillar Abeta(14-23) and Abeta(1-40) based on transferred nuclear Overhauser effect (trNOE) and transferred residual dipolar coupling (trRDC) data. In a first step, the inhibitor peptide structure is calculated on the basis of trNOE data; the trRDC data are then validated on the basis of the trNOE-derived structure using the program PALES. The orientation of the peptide inhibitors with respect to Abeta fibrils is obtained from trRDC data, assuming that Abeta fibrils orient such that the fibril axis is aligned in parallel with the magnetic field. The trRDC-derived alignment tensor of the peptide ligand is then used as a restraint for molecular dynamics docking studies. We find that the structure with the lowest rmsd value is in agreement with a model in which the inhibitor peptide binds to the long side of an amyloid fibril. Especially, we detect interactions involving the hydrophobic core, residues K16 and E22/D23 of the Abeta sequence. Structural differences are observed for binding of the inhibitor peptide to Abeta14-23 and Abeta1-40 fibrils, respectively, indicating different fibril structure. We expect this approach to be useful in the rational design of amyloid ligands with improved binding characteristics.     

Synthesis of C-ring-modified blebbistatin derivatives and evaluation of their myosin II ATPase inhibitory potency

(S)-Blebbistatin is a micromolar myosin II ATPase inhibitor that is extensively used in research. In search of analogs with improved potency, we have synthesized for the first time C-ring modified analogs. We introduced hydroxymethyl or allyloxymethyl functionalities in search of additional favorable interactions and a more optimal filling of the binding pocket. Unfortunately, the resulting compounds did not significantly inhibit the ATPase activity of rabbit skeletal-muscle myosin II. This and earlier reports suggest that rational design of potent myosin II inhibitors based on the architecture of the blebbistatin binding pocket is an ineffective strategy.     

Crystal structures of thrombin complexed to a novel series of synthetic inhibitors containing a 5,5-trans-lactone template.

The binding modes of four active site-directed, acylating inhibitors of human alpha-thrombin have been determined using X-ray crystallography. These inhibitors (GR157368, GR166081, GR167088, and GR179849) are representatives of a series utilizing a novel 5, 5-trans-lactone template to specifically acylate Ser195 of thrombin, resulting in an acyl complex. In each case the crystal structure of the complex reveals a binding mode which is consistent with the formation of a covalent bond between the ring-opened lactone of the inhibitor and residue Ser195. Improvements in potency and selectivity of these inhibitors for thrombin are rationalized on the basis of the observed protein/inhibitor interactions identified in these complexes. Occupation of the thrombin S2 and S3 pockets is shown to be directly correlated with improved binding and a degree of selectivity. The binding mode of GR179849 to thrombin is compared with the thrombin/PPACK complex [Bode, W., Turk, D., and Karshikov, A. (1992) Protein Sci. 1, 426-471] as this represents the archetypal binding mode for a thrombin inhibitor. This series of crystal structures is the first to be reported of synthetic, nonpeptidic acylating inhibitors bound to thrombin and provides details of the molecular recognition features that resulted in nanomolar potency.     

Natural iminosugar (+)-lentiginosine inhibits ATPase and chaperone activity of hsp90

Heat shock protein 90 (Hsp90) is a significant target in the development of rational cancer therapy due to its role at the crossroads of multiple signaling pathways associated with cell proliferation and cell viability. The relevance of Hsp90 as a therapeutic target for numerous diseases states has prompted the identification and optimization of novel Hsp90 inhibitors as an emerging therapeutic strategy. We performed a screening aimed to identify novel Hsp90 inhibitors among several natural compounds and we focused on the iminosugar (+)-lentiginosine, a natural amyloglucosidases inhibitor, for its peculiar bioactivity profile. Characterization of Hsp90 inhibition was performed using a panel of chemical and biological approaches, including limited proteolysis, biochemical and cellular assays. Our result suggested that the middle domain of Hsp90, as opposed to its ATP-binding pocket, is a promising binding site for new classes of Hsp90 inhibitors with multi-target anti-cancer potential.