Modelling of factor Xa-inhibitor complexes: a computational flexible docking approach

An in-silico layer-by-layer adsorption study of the interaction between Rebaudioside A and the T1R2 human sweet taste receptor: modelling and biosensing perspectives

The human sweet taste receptor (T1R2) monomer-a member of the G-protein coupled receptor family that detects a wide variety of chemically and structurally diverse sweet tasting molecules, is known to pose a significant threat to human health. Protein that lack crystal structure is a challenge in structure-based protein design. This study focused on the interaction of the T1R2 monomer with rebaudioside A (Reb-A), a steviol glycoside with potential use as a natural sweetener using in-silico and biosensing methods. Herein, homology modelling, docking studies, and molecular dynamics simulations were applied to elucidate the interaction between Reb-A and the T1R2 monomer. In addition, the electrochemical sensing of the immobilised T1R2-Reb-A complex with zinc oxide nanoparticles (ZnONPs) and graphene oxide (GO) were assessed by testing the performance of multiwalled carbon nanotube (MWCNT) as an adsorbent experimentally. Results indicate a strong interaction between Reb-A and the T1R2 receptor, revealing the stabilizing interaction of the amino acids with the Reb-A by hydrogen bonds with the hydroxyl groups of the glucose moieties, along with a significant amount of hydrophobic interactions. Moreover, the presence of the MWCNT as an anchor confirms the adsorption strength of the T1R2-Reb-A complex onto the GO nanocomposite and supported with electrochemical measurements. Overall, this study could serve as a cornerstone in the development of electrochemical immunosensor for the detection of Reb-A, with applications in the food industry.

Design, synthesis, structure, in vitro cytotoxic activity evaluation and docking studies on target enzyme GSK-3β of new indirubin-3'-oxime derivatives

The addition of chalcone and amine components into indirubin-3′-oxime resulted in 15 new derivatives with high yields. Structures of new derivatives were also elucidated through 1D, 2D-NMR and HR-MS(ESI) spectra and X-ray crystallography. All designed compounds were screened for cytotoxic activity against four human cancer cell lines (HepG2, LU-1, SW480 and HL-60) and one human normal kidney cell line (HEK-293). Compound 6f exhibited the most marked cytotoxicity meanwhile cytotoxicity of compounds 6e, 6h and 6l was more profound toward cancer cell lines than toward normal cell. These new derivatives were further analyzed via molecular docking studies on GSK-3β enzyme. Docking analysis shows that most of the derivatives exhibited potential inhibition activity against GSK-3β with characteristic interacting residues in the binding site. The fast pulling of ligand scheme was then employed to refine the binding affinity and mechanism between ligands and GSK-3β enzyme. The computational results are expected to contribute to predicting enzyme target of the trial inhibitors and their possible interaction, from which the design of new cytotoxic agents could be created in the future.

An unexpected dynamic binding mode between coagulation factor X and Rivaroxaban reveals importance of flexibility in drug binding

Rivaroxaban (RIV) is a direct oral anticoagulant (DOAC) targeting activated coagulation factor X (FXa). An earlier study reported the F174A mutant of FXa resistant to a RIV-like inhibitor, Apixaban. In current study, the detailed molecular mechanism of the resistance has been explored by molecular dynamics simulations on the impaired interactions between RIV and FXa in the damaged S4 pocket of F174A mutant. Besides, an unexpected relative stable binding mode of S1'S1 was revealed, which required dynamic motions of Gln192 and Gln61 to allow the morpholinone moiety of RIV to shift into the S1' pocket and form strong interactions. These dynamic motions of RIV and critical residues might be important in drug design for direct inhibitors of coagulation factors.

Quantitative Correlation of Conformational Binding Enthalpy with Substrate Specificity of Serine Proteases

Members of the same protease family show different substrate specificity, even if they share identical folds, depending on the physiological processes they are part of. Here, we investigate the key factors for subpocket and global specificity of factor Xa, elastase, and granzyme B which despite all being serine proteases and sharing the chymotrypsin-fold show distinct substrate specificity profiles. We determined subpocket interaction potentials with GRID for static X-ray structures and an in silico generated ensemble of conformations. Subpocket interaction potentials determined for static X-ray structures turned out to be insufficient to explain serine protease specificity for all subpockets. Therefore, we generated conformational ensembles using molecular dynamics simulations. We identified representative binding site conformations using distance-based hierarchical agglomerative clustering and determined subpocket interaction potentials for each representative conformation of the binding site. Considering the differences in subpocket interaction potentials for these representative conformations as well as their abundance allowed us to quantitatively explain subpocket specificity for the nonprime side for all three example proteases on a molecular level. The methods to identify key regions determining subpocket specificity introduced in this study are directly applicable to other serine proteases, and the results provide starting points for new strategies in rational drug design.

Selection of peptide inhibitors for double-stranded RNA-dependent protein kinase PKR

Protein kinase inhibitors have been developed and applied as antitumor drugs. The majority of these inhibitors are derived from ATP analogs with limited specificity towards the kinase target. Here we present our proof-of-principle study on peptide inhibitors for kinases. Two peptides were selected by phage display against double-stranded RNA-dependent protein kinase (PKR). In vitro assay revealed that these peptides exhibit an inhibitory effect on PKR-catalyzed phosphorylation of the alpha subunit of eukaryotic initiation factor 2 (eIF2α). The peptides also interrupt PKR activity in cells infected by viruses, as PKR activation is one of the hallmarks of host response to viral infection. Kinetic study revealed that one of the peptides, named P1, is a competitive inhibitor for PKR, while the other, named P2, exhibits a more complicated pattern of inhibition on PKR activity. Fragment-based docking of the PKR-peptide complex suggests that P1 occupies the substrate pocket of PKR and thus inhibits the binding between PKR and eIF2α, whereas P2 sits near the substrate pocket. The computational model of PKR-peptide complex agrees with their kinetic behavior. We surmise that peptide inhibitors for kinases have higher specificity than ATP analogs, and that they provide promising leads for the optimization of kinase inhibitors.

Virtual screening: an in silico tool for interlacing the chemical universe with the proteome

In silico screening both in the forward (traditional virtual screening) and reverse sense (inverse virtual screening (IVS)) are helpful techniques for interlacing the chemical universe of small molecules with the proteome. The former, which is using a protein structure and a large chemical database, is well-known by the scientific community. We have chosen here to provide an overview on the latter, focusing on validation and target prioritization strategies. By comparing it to complementary or alternative wet-lab approaches, we put IVS in the broader context of chemical genomics, target discovery and drug design. By giving examples from the literature and an own example on how to validate the approach, we provide guidance on the issues related to IVS.

Withanone as an inhibitor of survivin: a potential drug candidate for cancer therapy

Survivin, the smallest inhibitor of apoptosis protein, which has been reported to be highly expressed in almost all known cancers, plays a dual role in survival as well as the proliferation of cancer cells. It inhibits apoptosis by inhibiting caspases as well as facilitating mitosis by becoming a part of chromosomal passenger complex through its BIR5 domain. Docking studies carried out with herbal ligand withanone derived from roots of Withania somnifera have shown strong binding affinity of -19.1088 kJ/mol with BIR5 domain of survivin and in turn interferes with inhibitory action against caspases and may lead to apoptosis. Binding of withanone at BIR5 domain of survivin may also interfere with chromosomal passenger complex and lead to halt the mitotic process within the cancer cell. Docking studies support various experimental outcomes and suggest withanone as a potential candidate molecule in cancer therapy.

Ashwagandha derived withanone targets TPX2-Aurora A complex: computational and experimental evidence to its anticancer activity

Cancer is largely marked by genetic instability. Specific inhibition of individual proteins or signalling pathways that regulate genetic stability during cell division thus hold a great potential for cancer therapy. The Aurora A kinase is a Ser/Thr kinase that plays a critical role during mitosis and cytokinesis and is found upregulated in several cancer types. It is functionally regulated by its interactions with TPX2, a candidate oncogene. Aurora A inhibitors have been proposed as anticancer drugs that work by blocking its ATP binding site. This site is common to other kinases and hence these inhibitors lack specificity for Aurora A inhibition in particular, thus advocating the need of some alternative inhibition route. Previously, we identified TPX2 as a cellular target for withanone that selectively kill cancer cells. By computational approach, we found here that withanone binds to TPX2-Aurora A complex. In experiment, withanone treatment to cancer cells indeed resulted in dissociation of TPX2-Aurora A complex and disruption of mitotic spindle apparatus proposing this as a mechanism of the anticancer activity of withanone. From docking analysis, non-formation/disruption of the active TPX2-Aurora A association complex could be discerned. Our MD simulation results suggesting the thermodynamic and structural stability of TPX2-Aurora A in complex with withanone further substantiates the binding. We report a computational rationale of the ability of naturally occurring withanone to alter the kinase signalling pathway in an ATP-independent manner and experimental evidence in which withanone cause inactivation of the TPX2-Aurora A complex. The study demonstrated that TPX2-Aurora A complex is a target of withanone, a potential natural anticancer drug.

A comparative reverse docking strategy to identify potential antineoplastic targets of tea functional components and binding mode

The main functional components of green tea, such as epigallocatechin gallate (EGCG), epigallocatechin (EGC), epicatechin gallate (ECG) and epicatechin (EC), are found to have a broad antineoplastic activity. The discovery of their targets plays an important role in revealing the antineoplastic mechanism. Therefore, to identify potential target proteins for tea polyphenols, we have taken a comparative virtual screening approach using two reverse docking systems, one based on Autodock software and the other on Tarfisdock. Two separate in silico workflows were implemented to derive a set of target proteins related to human diseases and ranked by the binding energy score. Several conventional clinically important proteins with anti-tumor effects are screened out from the PDTD protein database as the potential receptors by both procedures. To further analyze the validity of docking results, we study the binding mode of EGCG and the potential target protein Leukotriene A4 hydrolase in detail. We indicate that interactions mediated by electrostatic and hydrogen bond play a key role in ligand binding. EGCG binds to the enzyme with certain orientation and conformation that is suitable for nucleophilic attacks by several electrical residues inside the enzyme's activity cavity. This study provides useful information for studying the antitumor mechanism of tea's functional components. The comparative reverse docking strategy presented generates a tractable set of antineoplastic proteins for future experimental validation as drug targets against tumors.

Hsp90/Cdc37 chaperone/co-chaperone complex, a novel junction anticancer target elucidated by the mode of action of herbal drug Withaferin A

BACKGROUND

HSPs (Heat shock proteins) are highly conserved ubiquitous proteins among species which are involved in maintaining appropriate folding and conformation of other proteins and are thus referred to as molecular chaperones. Hsp90 (Heat-shock protein 90 kDa) is one of a group of molecular chaperones responsible for managing protein folding and quality control in cell environment. However it is also involved in the maturation and stabilization of a wide range of oncogenic client proteins which are crucial for oncogenesis and malignant progression. Hsp90 requires a series of co-chaperones to assemble into a super-chaperone complex for its function. These co-chaperones bind and leave the complex at various stages to regulate the chaperoning process. Arresting the chaperone cycle at these stages by targeting different co-chaperone/Hsp90 interactions seems to be quite a viable alternative and is likely to achieve similar consequences as that of Hsp90 direct inhibition with added favors of high specificity and reduced side effect profile. The study conducted here is an attempt to explore the potential of Withania somnifera's major constituent WA (Withaferin A) in attenuating the Hsp90/Cdc37 chaperone/co-chaperone interactions for enhanced tumor arresting activity and to elucidate the underlying mode of action using computational approaches.

RESULTS

Formation of active Hsp90/Cdc37 complex is one of the essential steps for facilitation of chaperone client interaction, non-assembly of which can lead to prevention of the chaperone-client association resulting in apoptosis of tumor cells. From our flexible docking analysis of WA into active Hsp90/Cdc37 complex in which key interfacing residues of the complex were kept flexible, disruption of the active association complex can be discerned. While docking of WA into segregated Hsp90 leaves the interface residues untouched. Thus the molecular docking analysis of WA into Hsp90 and active Hsp90/Cdc37 complex conducted in this study provides significant evidence in support of the proposed mechanism of chaperone assembly suppression by inhibition or disruption of active Hsp90/Cdc37 complex formation being accounted by non-assembly of the catalytically active Hsp90/Cdc37 complex. Results from the molecular dynamics simulations in water show that the trajectories of the protein complexed with ligand WA are stable over a considerably long time period of 4 ns, with the energies of the complex being lowered in comparison to the un-docked association complex, suggesting the thermodynamic stability of WA complexed Hsp90/Cdc37.

CONCLUSIONS

The molecular chaperone Hsp90 has been a promising target for cancer therapy. Cancer is a disease marked by genetic instability. Thus specific inhibition of individual proteins or signalling pathways holds a great potential for subversion of this genetic plasticity of cancers. This study is a step forward in this direction. Our computational analysis provided a rationalization to the ability of naturally occurring WA to alter the chaperone signalling pathway. The large value of binding energy involved in binding of WA to the active Hsp90/Cdc37 complex consolidates the thermodynamic stability of the binding. Our docking results obtained substantiate the hypothesis that WA has the potential to inhibit the association of chaperone (Hsp90) to its co-chaperone (Cdc37) by disrupting the stability of attachment of Hsp90 to Cdc37. Conclusively our results strongly suggest that withaferin A is a potent anticancer agent as ascertained by its potent Hsp90-client modulating capability.

Inhibition of the NEMO/IKKβ association complex formation, a novel mechanism associated with the NF-κB activation suppression by Withania somnifera's key metabolite withaferin A

Background

Nuclear Factor kappa B (NF-κB) is a transcription factor involved in the regulation of cell signaling responses and is a key regulator of cellular processes involved in the immune response, differentiation, cell proliferation, and apoptosis. The constitutive activation of NF-κB contributes to multiple cellular outcomes and pathophysiological conditions such as rheumatoid arthritis, asthma, inflammatory bowel disease, AIDS and cancer. Thus there lies a huge therapeutic potential beneath inhibition of NF-κB signalling pathway for reducing these chronic ailments. Withania somnifera, a reputed herb in ayurvedic medicine, comprises a large number of steroidal lactones known as withanolides which show plethora of pharmacological activities like anti- inflammatory, antitumor, antibacterial, antioxidant, anticonvulsive, and immunosuppressive. Though a few studies have been reported depicting the effect of WA (withaferin A) on suppression of NF-κB activation, the mechanism behind this is still eluding the researchers. The study conducted here is an attempt to explore NF-κB signalling pathway modulating capability of Withania somnifera’s major constituent WA and to elucidate its possible mode of action using molecular docking and molecular dynamics simulations studies.

Results

Formation of active IKK (IκB kinase) complex comprising NEMO (NF-κB Essential Modulator) and IKKβ subunits is one of the essential steps for NF-κB signalling pathway, non-assembly of which can lead to prevention of the above mentioned vulnerable disorders. As observed from our semi-flexible docking analysis, WA forms strong intermolecular interactions with the NEMO chains thus building steric as well as thermodynamic barriers to the incoming IKKβ subunits, which in turn pave way to naive complex formation capability of NEMO with IKKβ. Docking of WA into active NEMO/IKKβ complex using flexible docking in which key residues of the complex were kept flexible also suggest the disruption of the active complex. Thus the molecular docking analysis of WA into NEMO and active NEMO/IKKβ complex conducted in this study provides significant evidence in support of the proposed mechanism of NF-κB activation suppression by inhibition or disruption of active NEMO/IKKβ complex formation being accounted by non-assembly of the catalytically active NEMO/IKKβ complex. Results from the molecular dynamics simulations in water show that the trajectories of the native protein and the protein complexed with WA are stable over a considerably long time period of 2.6 ns.

Conclusions

NF-κB is one of the most attractive topics in current biological, biochemical, and pharmacological research, and in the recent years the number of studies focusing on its inhibition/regulation has increased manifolds. Small ligands (both natural and synthetic) are gaining particular attention in this context. Our computational analysis provided a rationalization of the ability of naturally occurring withaferin A to alter the NF-κB signalling pathway along with its proposed mode of inhibition of the pathway. The absence of active IKK multisubunit complex would prevent degradation of IκB proteins, as the IκB proteins would not get phosphorylated by IKK. This would ultimately lead to non-release of NF-κB and its further translocation to the nucleus thus arresting its nefarious acts. Conclusively our results strongly suggest that withaferin A is a potent anticancer agent as ascertained by its potent NF-κB modulating capability. Moreover the present MD simulations made clear the dynamic structural stability of NEMO/IKKβ in complex with the drug WA, together with the inhibitory mechanism.

Stabilizing of a globular protein by a highly complex water network: a molecular dynamics simulation study on factor Xa

The role of water molecules is increasingly attracting attention in structural biology, and many studies have demonstrated their crucial contribution to the stability and function of proteins. Here, we present molecular dynamics studies on factor Xa (fXa) to investigate the effect of water molecules in this serine protease. fXa is a key enzyme in the blood coagulation cascade, and thus, an important target for antithrombotic drugs. A reasonable representation of the structure is crucial for an investigation at the molecular level and, thus, a prerequisite for structure-based drug design. Simulations of well-resolved fXa X-ray structures with different sets of water molecules show the importance of a well-determined water set for the simulation. We discuss implications of different water sets on the structure and dynamics of fXa.

Computational modeling of the potential interactions of the proteasome beta5 subunit and catechol-O-methyltransferase-resistant EGCG analogs

(-)-Epigallocatechin gallate [(-)-EGCG] has been implicated in cancer chemoprevention and has been shown as an inhibitor of tumor proteasomal chymotrypsin-like activity in vitro and in vivo. However, EGCG is subjected to rapid biotransforming modifications such as methylation by catechol-Omicron-methyltransferase (COMT) that limits its action. We recently reported that structure 7, an EGCG analog which should be resistant to COMT-mediated methylation and inactivation in cells, was able to inhibit the activity of purified 20S proteasome and cellular 26S proteasome. However, the involved molecular mechanism is unknown. Herein, we applied computational solution to understand the possible interaction between EGCG analogs including structure 7 and the proteasome beta5 subunit which is responsible for the chymotrypsin-like activity. We report that the ester carbonyls at C2 and C3 carbon atoms may be the active sites for nucleophilic attack in structure 7 and 5. Equally spaced carbon atoms in COMT-resistant structure 7 give more stable conformation and lower docked free energy than other EGCG analogs. The absence of a second gallate group in structure 16 and 21 significantly decreases the ability to inhibit the proteasome.

Prediction of factor Xa inhibitors by machine learning methods

Factor Xa (FXa) inhibitors have been explored as anticoagulants for treatment and prevention of thrombotic diseases. Molecular docking, pharmacophore, quantitative structure-activity relationships, and support vector machines (SVM) have been used for computer prediction of FXa inhibitors. These methods achieve promising prediction accuracies of 69-80% for FXa inhibitors and 85-99% for non-inhibitors. Prediction performance, particularly for inhibitors, may be further improved by exploring methods applicable to more diverse range of compounds and by using more appropriate set of molecular descriptors. We tested the capability of several machine learning methods (C4.5 decision tree, k-nearest neighbor, probabilistic neural network, and support vector machine) by using a much more diverse set of 1098 compounds (360 inhibitors and 738 non-inhibitors) than those in other studies. A feature selection method was used for selecting molecular descriptors appropriate for distinguishing FXa inhibitors and non-inhibitors. The prediction accuracies of these methods are 89.1-97.5% for FXa inhibitors and 92.3-98.1% for non-inhibitors. In particular, compared to other studies, support vector machine gives a substantially improved accuracy of 94.6% for FXa non-inhibitors and maintains a comparable accuracy of 98.1% for inhibitors, based-on a more rigorous test with more diverse range of compounds. Our study suggests that machine learning methods such as SVM are useful for facilitating the prediction of FXa inhibitors.

Structure-based activity prediction for an enzyme of unknown function

With many genomes sequenced, a pressing challenge in biology is predicting the function of the proteins that the genes encode. When proteins are unrelated to others of known activity, bioinformatics inference for function becomes problematic. It would thus be useful to interrogate protein structures for function directly. Here, we predict the function of an enzyme of unknown activity, Tm0936 from Thermotoga maritima, by docking high-energy intermediate forms of thousands of candidate metabolites. The docking hit list was dominated by adenine analogues, which appeared to undergo C6-deamination. Four of these, including 5-methylthioadenosine and S-adenosylhomocysteine (SAH), were tested as substrates, and three had substantial catalytic rate constants (10(5) M(-1 )s(-1)). The X-ray crystal structure of the complex between Tm0936 and the product resulting from the deamination of SAH, S-inosylhomocysteine, was determined, and it corresponded closely to the predicted structure. The deaminated products can be further metabolized by T. maritima in a previously uncharacterized SAH degradation pathway. Structure-based docking with high-energy forms of potential substrates may be a useful tool to annotate enzymes for function.

E230Q mutation of the catalytic subunit of cAMP-dependent protein kinase affects local structure and the binding of peptide inhibitor

The active site of the mammalian cAMP-dependent protein kinase catalytic subunit (C-subunit) has a cluster of nonconserved acidic residues-Glu127, Glu170, Glu203, Glu230, and Asp241-that are crucial for substrate recognition and binding. Studies have shown that the Glu230 to Gln mutant (E230Q) of the enzyme has physical properties similar to the wild-type enzyme and has decreased affinity for a short peptide substrate, Kemptide. However, recent experiments intended to crystallize ternary complex of the E230Q mutant with MgATP and protein kinase inhibitor (PKI) could only obtain crystals of the apo-enzyme of E230Q mutant. To deduce the possible mechanism that prevented ternary complex formation, we used the relaxed-complex method (Lin, J.-H., et al. J Am Chem Soc 2002, 24, 5632-5633) to study PKI binding to the E230Q mutant C-subunit. In the E230Q mutant, we observed local structural changes of the peptide binding site that correlated closely to the reduced PKI affinity. The structural changes occurred in the F-to-G helix loop and appeared to hinder PKI binding. Reduced electrostatic potential repulsion among Asp241 from the helix loop section and the other acidic residues in the peptide binding site appear to be responsible for the structural change.

Exploring structurally conserved solvent sites in protein families

Protein-bound water molecules are important components of protein structure, and therefore, protein function and energetics. Although structural conservation of solvent has been studied in a few protein families, a lack of suitable computational tools has hindered more comprehensive analyses. Herein we present a semiautomated computational approach for identifying solvent sites that are conserved among proteins sharing a common three-dimensional structure. This method is tested on six protein families: (1) monodomain cytochrome c, (2) fatty-acid binding protein, (3) lactate/malate dehydrogenase, (4) parvalbumin, (5) phospholipase A2, and (6) serine protease. For each family, the method successfully identified previously known conserved solvent sites. Moreover, the method discovered 22 novel conserved solvent sites, some of which have higher degrees of conservation than the previously known sites. All six families studied had solvent sites with more than 90% conservation and these sites were invariably located in regions of the protein with very high sequence conservation. These results suggest that highly conserved solvent sites, by virtue of their proximity to conserved residues, should be considered as one of the defining three-dimensional structural characteristics of protein families and folds.

Decoys for docking

Molecular docking is widely used to predict novel lead compounds for drug discovery. Success depends on the quality of the docking scoring function, among other factors. An imperfect scoring function can mislead by predicting incorrect ligand geometries or by selecting nonbinding molecules over true ligands. These false-positive hits may be considered "decoys". Although these decoys are frustrating, they potentially provide important tests for a docking algorithm; the more subtle the decoy, the more rigorous the test. Indeed, decoy databases have been used to improve protein structure prediction algorithms and protein-protein docking algorithms. Here, we describe 20 geometric decoys in five enzymes and 166 "hit list" decoys-i.e., molecules predicted to bind by our docking program that were tested and found not to do so-for beta-lactamase and two cavity sites in lysozyme. Especially in the cavity sites, which are very simple, these decoys highlight particular weaknesses in our scoring function. We also consider the performance of five other widely used docking scoring functions against our geometric and hit list decoys. Intriguingly, whereas many of these other scoring functions performed better on the geometric decoys, they typically performed worse on the hit list decoys, often highly ranking molecules that seemed to poorly complement the model sites. Several of these "hits"from the other scoring functions were tested experimentally and found, in fact, to be decoys. Collectively, these decoys provide a tool for the development and improvement of molecular docking scoring functions. Such improvements may, in turn, be rapidly tested experimentally against these and related experimental systems, which are well-behaved in assays and for structure determination.

Dietary flavonoids as proteasome inhibitors and apoptosis inducers in human leukemia cells

It has been shown that proteasome activity is required for cancer cell survival and consumption of fruits and vegetables is associated with decreased cancer risk. Previously, we reported that grape extract could inhibit proteasome activity and induce apoptosis in tumor cells. In this study, we examined the flavonoids apigenin, quercetin, kaempferol and myricetin for their proteasome-inhibitory and apoptosis-inducing abilities in human tumor cells. We report that apigenin and quercetin are much more potent than kaempferol and myricetin at: (i) inhibiting chymotrypsin-like activity of purified 20S proteasome and of 26S proteasome in intact leukemia Jurkat T cells; (ii) accumulating putative ubiquitinated forms of two proteasome target proteins, Bax and Inhibitor of nuclear factor kappabeta-alpha in Jurkat T cells and (iii) inducing activation of caspase-3 and cleavage of poly(ADP-ribose) polymerase in Jurkat T cells. The proteasome-inhibitory abilities of these compounds correlated with their apoptosis-inducing potencies. Results from computational modeling of the potential interactions of these flavonoids to the chymotrypsin site (beta5 subunit) of the proteasome were consistent with the obtained proteasome-inhibitory activities. We found that the C(4) carbon may be a site of nucleophilic attack by the OH group of N-terminal threonine of proteasomal beta5 subunit and that the C(3) hydroxyl may alter the ability of these flavonoids to inhibit the proteasome. Finally, apigenin neither effectively inhibited the proteasome activity nor induced apoptosis in non-transformed human natural killer cells. Our results suggested that the proteasome may be a target of these dietary flavonoids in human tumor cells and that inhibition of the proteasome by flavonoids may be one of the mechanisms responsible for their cancer-preventive effects.

Automated docking of ligands to an artificial active site: augmenting crystallographic analysis with computer modeling

The W191G cavity of cytochrome c peroxidase is useful as a model system for introducing small molecule oxidation in an artificially created cavity. A set of small, cyclic, organic cations was previously shown to bind in the buried, solvent-filled pocket created by the W191G mutation. We docked these ligands and a set of non-binders in the W191G cavity using AutoDock 3.0. For the ligands, we compared docking predictions with experimentally determined binding energies and X-ray crystal structure complexes. For the ligands, predicted binding energies differed from measured values by +/- 0.8 kcal/mol. For most ligands, the docking simulation clearly predicted a single binding mode that matched the crystallographic binding mode within 1.0 A RMSD. For 2 ligands, where the docking procedure yielded an ambiguous result, solutions matching the crystallographic result could be obtained by including an additional crystallographically observed water molecule in the protein model. For the remaining 2 ligands, docking indicated multiple binding modes, consistent with the original electron density, suggesting disordered binding of these ligands. Visual inspection of the atomic affinity grid maps used in docking calculations revealed two patches of high affinity for hydrogen bond donating groups. Multiple solutions are predicted as these two sites compete for polar hydrogens in the ligand during the docking simulation. Ligands could be distinguished, to some extent, from non-binders using a combination of two trends: predicted binding energy and level of clustering. In summary, AutoDock 3.0 appears to be useful in predicting key structural and energetic features of ligand binding in the W191G cavity.

Docking studies and model development of tea polyphenol proteasome inhibitors: applications to rational drug design

Previously, we demonstrated that natural and synthetic ester bond-containing green tea polyphenols were potent and specific non-peptide proteasome inhibitors. However, the molecular mechanism of inhibition is currently unknown. Here, we report that inhibition of the chymotrypsin activity of the 20S proteasome by (-)-epigallocatechin-3-gallate (EGCG) is time-dependent and irreversible, implicating acylation of the beta5-subunit's catalytic N-terminal threonine (Thr 1). This knowledge is used, along with in silico docking experiments, to aid in the understanding of binding and inhibition. On the basis of these docking experiments, we propose that (-)-EGCG binds the chymotrypsin site in an orientation and conformation that is suitable for a nucleophilic attack by Thr 1. Consistently, the distance from the electrophilic carbonyl carbon of (-)-EGCG to the hydroxyl group of Thr 1 was measured as 3.18 A. Furthermore, the A ring of (-)-EGCG acts as a tyrosine mimic, binding to the hydrophobic S1 pocket of the beta5-subunit. In the process, the (-)-EGCG scissile bond may become strained, which could lower the activation energy for attack by the hydroxyl group of Thr 1. This model is validated by comparison of predicted and actual activities of several EGCG analogs, either naturally occurring, previously synthesized, or rationally synthesized.

Inhibition of the proteasome activity, a novel mechanism associated with the tumor cell apoptosis-inducing ability of genistein

Epidemiological studies have suggested that increased soy consumption is associated with reduced cancer occurrence. Genistein, a soy isoflavone, has been reported to inhibit the growth of human tumor cells although the involved molecular mechanisms are not clearly defined. Here we report that genistein inhibits the proteasomal chymotrypsin-like activity in vitro and in vivo. Computational docking studies suggest that the interaction of genistein with the proteasomal beta 5 subunit is responsible for inhibition of the chymotrypsin-like activity. Inhibition of the proteasome by genistein in prostate cancer LNCaP and breast cancer MCF-7 cells is associated with accumulation of ubiquitinated proteins and three known proteasome target proteins, the cyclin-dependent kinase inhibitor p27(Kip1), inhibitor of nuclear factor-kappa B (I kappa B-alpha), and the pro-apoptotic protein Bax. Genistein-mediated proteasome inhibition was accompanied by induction of apoptosis in these solid tumor cells. Finally, genistein induced proteasome inhibition and apoptosis selectively in simian virus 40-transformed human fibroblasts, but not in their parental normal counterpart. Our results suggest that the proteasome is a potential target of genistein in human tumor cells and that inhibition of the proteasome activity by genistein might contribute to its cancer-preventive properties.

Binding affinity prediction of novel estrogen receptor ligands using receptor-based 3-D QSAR methods

We have recently reported the development of a 3-D QSAR model for estrogen receptor ligands showing a significant correlation between calculated molecular interaction fields and experimentally measured binding affinity. The ligand alignment obtained from docking simulations was taken as basis for a comparative field analysis applying the GRID/GOLPE program. Using the interaction field derived with a water probe and applying the smart region definition (SRD) variable selection procedure, a significant and robust model was obtained (q(2)(LOO)=0.921, SDEP=0.345). To further analyze the robustness and the predictivity of the established model several recently developed estrogen receptor ligands were selected as external test set. An excellent agreement between predicted and experimental binding data was obtained indicated by an external SDEP of 0.531. Two other traditionally used prediction techniques were applied in order to check the performance of the receptor-based 3-D QSAR procedure. The interaction energies calculated on the basis of receptor-ligand complexes were correlated with experimentally observed affinities. Also ligand-based 3-D QSAR models were generated using program FlexS. The interaction energy-based model, as well as the ligand-based 3-D QSAR models yielded models with lower predictivity. The comparison with the interaction energy-based model and with the ligand-based 3-D QSAR models, respectively, indicates that the combination of receptor-based and 3-D QSAR methods is able to improve the quality of prediction.

Models of the extracellular domain of the nicotinic receptors and of agonist- and Ca2+-binding sites

We constructed a three-dimensional model of the amino-terminal extracellular domain of three major types of nicotinic acetylcholine receptor, (alpha7)5, (alpha4)2(beta2)3, and (alpha1)2beta1gammadelta, on the basis of the recent x-ray structure determination of the molluscan acetylcholine-binding protein. Comparative analysis of the three models reveals that the agonist-binding pocket is much more conserved than the overall structure. Differences exist, however, in the side chains of several residues. In particular, a phenylalanine residue, present in beta2 but not in alpha7, is proposed to contribute to the high affinity for agonists in receptors containing the beta2 subunit. The semiautomatic docking of agonists in the ligand-binding pocket of (alpha7)5 led to positions consistent with labeling and mutagenesis experiments. Accordingly, the quaternary ammonium head group of nicotine makes a pi-cation interaction with W148 (alpha7 numbering), whereas the pyridine ring is close to both the cysteine pair 189-190 and the complementary component of the binding site. The intrinsic affinities inferred from docking give a rank order epibatidine > nicotine > acetylcholine, in agreement with experimental values. Finally, our models offer a structural basis for potentiation by external Ca2+.

Structure-based 3D QSAR and design of novel acetylcholinesterase inhibitors

The paper describes the construction, validation and application of a structure-based 3D QSAR model of novel acetylcholinesterase (AChE) inhibitors. Initial use was made of four X-ray structures of AChE complexed with small, non-specific inhibitors to create a model of the binding of recently developed aminopyridazine derivatives. Combined automated and manual docking methods were applied to dock the co-crystallized inhibitors into the binding pocket. Validation of the modelling process was achieved by comparing the predicted enzyme-bound conformation with the known conformation in the X-ray structure. The successful prediction of the binding conformation of the known inhibitors gave confidence that we could use our model to evaluate the binding conformation of the aminopyridazine compounds. The alignment of 42 aminopyridazine compounds derived by the docking procedure was taken as the basis for a 3D QSAR analysis applying the GRID/GOLPE method. A model of high quality was obtained using the GRID water probe, as confirmed by the cross-validation method (q2LOO = 0.937, q2L50%O = 0.910). The validated model, together with the information obtained from the calculated AChE-inhibitor complexes, were considered for the design of novel compounds. Seven designed inhibitors which were synthesized and tested were shown to be highly active. After performing our modelling study the X-ray structure of AChE complexed with donepezil, an inhibitor structurally related to the developed aminopyirdazines, has been made available. The good agreement found between the predicted binding conformation of the aminopyridazines and the one observed for donepezil in the crystal structure further supports our developed model.

Molecular mechanism for agonist-promoted alpha(2A)-adrenoceptor activation by norepinephrine and epinephrine

We present a mechanism for agonist-promoted alpha(2A)-adrenergic receptor (alpha(2A)-AR) activation based on structural, pharmacological, and theoretical evidence of the interactions between phenethylamine ligands and alpha(2A)-AR. In this study, we have: 1) isolated enantiomerically pure phenethylamines that differ both in their chirality about the beta-carbon, and in the presence/absence of one or more hydroxyl groups: the beta-OH and the catecholic meta- and para-OH groups; 2) used [(3)H]UK-14,304 [5-bromo-N-(4,5-dihydro-1H-imidazol-2-yl)-6-quinoxalinamine; agonist] and [(3)H]RX821002 [2-(2-methoxy-1,4-benzodioxan-2-yl)-2-imidazoline; antagonist] competition binding assays to determine binding affinities of these ligands to the high- and low-affinity forms of alpha(2A)-AR; 3) tested the ability of the ligands to promote receptor activation by measuring agonist-induced stimulation of [(35)S]GTPgammaS binding in isolated cell membranes; and 4) used automated docking methods and our alpha(2A)-AR model to predict the binding modes of the ligands inside the alpha(2A)-AR binding site. The ligand molecules are sequentially missing different functional groups, and we have correlated the structural features of the ligands and ligand-receptor interactions with experimental ligand binding and receptor activation data. Based on the analysis, we show that structural rearrangements in transmembrane helix (TM) 5 could take place upon binding and subsequent activation of alpha(2A)-AR by phenethylamine agonists. We suggest that the following residues are important in phenethylamine interactions with alpha(2A)-AR: Asp113 (D(3.32)), Val114 (V(3.33)), and Thr118 (T(3.37)) in TM3; Ser200 (S(5.42)), Cys201 (C(5.43)), and Ser204 (S(5.46)) in TM5; Phe391 (F(6.52)) and Tyr394 (Y(6.55)) in TM6; and Phe411 (F(7.38)) and Phe412 (F(7.39)) in TM7.

Factor Xa: simulation studies with an eye to inhibitor design

Factor Xa is a serine protease which activates thrombin and plays a key regulatory role in the blood-coagulation cascade. Factor Xa is at the crossroads of the extrinsic and intrinsic pathways of coagulation and, hence, has become an important target for the design of anti-thrombotics (inhibitors). It is not known to be involved in other processes than hemostasis and its binding site is different to that of other serine proteases, thus facilitating selective inhibition. The design of high-affinity selective inhibitors of factor Xa requires knowledge of the structural and dynamical characteristics of its active site. The three-dimensional structure of factor Xa was resolved by X-ray crystallography and refined at 2.2 A resolution by Padmanabhan and collaborators. In this article we present results from molecular dynamics simulations of the catalytic domain of factor Xa in aqueous solution. The simulations were performed to characterise the mobility and flexibility of the residues delimiting the unoccupied binding site of the enzyme, and to determine hydrogen bonding propensities (with protein and with solvent atoms) of those residues in the active site that could interact with a substrate or a potential inhibitor. The simulation data is aimed at facilitating the design of high-affinity selective inhibitors of factor Xa.

Receptor-based 3D QSAR analysis of estrogen receptor ligands--merging the accuracy of receptor-based alignments with the computational efficiency of ligand-based methods

One of the major challenges in computational approaches to drug design is the accurate prediction of binding affinity of biomolecules. In the present study several prediction methods for a published set of estrogen receptor ligands are investigated and compared. The binding modes of 30 ligands were determined using the docking program AutoDock and were compared with available X-ray structures of estrogen receptor-ligand complexes. On the basis of the docking results an interaction energy-based model, which uses the information of the whole ligand-receptor complex, was generated. Several parameters were modified in order to analyze their influence onto the correlation between binding affinities and calculated ligand-receptor interaction energies. The highest correlation coefficient (r2 = 0.617, q2Loo = 0.570) was obtained considering protein flexibility during the interaction energy evaluation. The second prediction method uses a combination of receptor-based and 3D quantitative structure-activity relationships (3D QSAR) methods. The ligand alignment obtained from the docking simulations was taken as basis for a comparative field analysis applying the GRID/GOLPE program. Using the interaction field derived with a water probe and applying the smart region definition (SRD) variable selection, a significant and robust model was obtained (r2 = 0.991, q2LOO = 0.921). The predictive ability of the established model was further evaluated by using a test set of six additional compounds. The comparison with the generated interaction energy-based model and with a traditional CoMFA model obtained using a ligand-based alignment (r2 = 0.951, q2L00 = 0.796) indicates that the combination of receptor-based and 3D QSAR methods is able to improve the quality of the underlying model.

DX-9065a, a specific factor Xa inhibitor, as a universal anticoagulant for blood collection tubes

DX-9065a is a direct and selective factor Xa inhibitor. We evaluated the usefulness of DX-9065a in anticoagulating specimens for routine laboratory tests. Results using blood anticoagulated with DX9065a corresponded well with results with blood treated with ethylendiamine tetraacetic acid (EDTA) in the complete blood count (CBC), including white blood cell (WBC) differential count and morphology of blood cells. Clinical chemistry results from DX-9065a-treated samples were similar to results obtained with serum specimens except for leucine aminopeptidase (LAP) and cholinesterase (Ch-E). Thus, DX-9065a may be a useful universal anticoagulant for laboratory medicine.