Novel Mass Spectrometric Analysis Methods for Anabolic Androgenic Steroids in Sports Drug Testing

The feasibility of different modern analytical techniques for the mass spectrometric detection of anabolic androgenic steroids (AAS) in human urine was examined in order to enhance the prevalent analytics and to find reasonable strategies for effective sports drug testing. A comparative study of the sensitivity and specificity between gas chromatography (GC) combined with low (LRMS) and high resolution mass spectrometry (HRMS) in screening of AAS was carried out with four metabolites of methandienone. Measurements were done in selected ion monitoring mode with HRMS using a mass resolution of 5000. With HRMS the detection limits were considerably lower than with LRMS, enabling detection of steroids at low 0.2-0.5 ng/ml levels. However, also with HRMS, the biological background hampered the detection of some steroids. The applicability of liquid-phase microextraction (LPME) was studied with metabolites of fluoxymesterone, 4-chlorodehydromethyltestosterone, stanozolol and danazol. Factors affecting the extraction process were studied and a novel LPME method with in-fiber silylation was developed and validated for GC/MS analysis of the danazol metabolite. The method allowed precise, selective and sensitive analysis of the metabolite and enabled simultaneous filtration, extraction, enrichment and derivatization of the analyte from urine without any other steps in sample preparation. Liquid chromatographic/tandem mass spectrometric (LC/MS/MS) methods utilizing electrospray ionization (ESI), atmospheric pressure chemical ionization (APCI) and atmospheric pressure photoionization (APPI) were developed and applied for detection of oxandrolone and metabolites of stanozolol and 4-chlorodehydromethyltestosterone in urine. All methods exhibited high sensitivity and specificity. ESI showed, however, the best applicability, and a LC/ESI-MS/MS method for routine screening of nine 17-alkyl-substituted AAS was thus developed enabling fast and precise measurement of all analytes with detection limits below 2 ng/ml. The potential of chemometrics to resolve complex GC/MS data was demonstrated with samples prepared for AAS screening. Acquired full scan spectral data (m/z 40-700) were processed by the OSCAR algorithm (Optimization by Stepwise Constraints of Alternating Regression). The deconvolution process was able to dig out from a GC/MS run more than the double number of components as compared with the number of visible chromatographic peaks. Severely overlapping components, as well as components hidden in the chromatographic background could be isolated successfully. All studied techniques proved to be useful analytical tools to improve detection of AAS in urine. Superiority of different procedures is, however, compound-dependent and different techniques complement each other.

From Polymorph Screening to Dissolution Testing - Solid Phase Analysis during Pharmaceutical Development and Manufacturing

Solid materials can exist in different physical structures without a change in chemical composition. This phenomenon, known as polymorphism, has several implications on pharmaceutical development and manufacturing. Various solid forms of a drug can possess different physical and chemical properties, which may affect processing characteristics and stability, as well as the performance of a drug in the human body. Therefore, knowledge and control of the solid forms is fundamental to maintain safety and high quality of pharmaceuticals. During manufacture, harsh conditions can give rise to unexpected solid phase transformations and therefore change the behavior of the drug. Traditionally, pharmaceutical production has relied on time-consuming off-line analysis of production batches and finished products. This has led to poor understanding of processes and drug products. Therefore, new powerful methods that enable real time monitoring of pharmaceuticals during manufacturing processes are greatly needed. The aim of this thesis was to apply spectroscopic techniques to solid phase analysis within different stages of drug development and manufacturing, and thus, provide a molecular level insight into the behavior of active pharmaceutical ingredients (APIs) during processing. Applications to polymorph screening and different unit operations were developed and studied. A new approach to dissolution testing, which involves simultaneous measurement of drug concentration in the dissolution medium and in-situ solid phase analysis of the dissolving sample, was introduced and studied. Solid phase analysis was successfully performed during different stages, enabling a molecular level insight into the occurring phenomena. Near-infrared (NIR) spectroscopy was utilized in screening of polymorphs and processing-induced transformations (PITs). Polymorph screening was also studied with NIR and Raman spectroscopy in tandem. Quantitative solid phase analysis during fluidized bed drying was performed with in-line NIR and Raman spectroscopy and partial least squares (PLS) regression, and different dehydration mechanisms were studied using in-situ spectroscopy and partial least squares discriminant analysis (PLS-DA). In-situ solid phase analysis with Raman spectroscopy during dissolution testing enabled analysis of dissolution as a whole, and provided a scientific explanation for changes in the dissolution rate. It was concluded that the methods applied and studied provide better process understanding and knowledge of the drug products, and therefore, a way to achieve better quality.

Liquid Chromatography-Mass Spectrometry in Studies of Drug Metabolism and Permeability

Poor pharmacokinetics is one of the reasons for the withdrawal of drug candidates from clinical trials. There is an urgent need for investigating in vitro ADME (absorption, distribution, metabolism and excretion) properties and recognising unsuitable drug candidates as early as possible in the drug development process. Current throughput of in vitro ADME profiling is insufficient because effective new synthesis techniques, such as drug design in silico and combinatorial synthesis, have vastly increased the number of drug candidates. Assay technologies for larger sets of compounds than are currently feasible are critically needed. The first part of this work focused on the evaluation of cocktail strategy in studies of drug permeability and metabolic stability. N-in-one liquid chromatography-tandem mass spectrometry (LC/MS/MS) methods were developed and validated for the multiple component analysis of samples in cocktail experiments. Together, cocktail dosing and LC/MS/MS were found to form an effective tool for increasing throughput. First, cocktail dosing, i.e. the use of a mixture of many test compounds, was applied in permeability experiments with Caco-2 cell culture, which is a widely used in vitro model for small intestinal absorption. A cocktail of 7-10 reference compounds was successfully evaluated for standardization and routine testing of the performance of Caco-2 cell cultures. Secondly, cocktail strategy was used in metabolic stability studies of drugs with UGT isoenzymes, which are one of the most important phase II drug metabolizing enzymes. The study confirmed that the determination of intrinsic clearance (Clint) as a cocktail of seven substrates is possible. The LC/MS/MS methods that were developed were fast and reliable for the quantitative analysis of a heterogenous set of drugs from Caco-2 permeability experiments and the set of glucuronides from in vitro stability experiments. The performance of a new ionization technique, atmospheric pressure photoionization (APPI), was evaluated through comparison with electrospray ionization (ESI), where both techniques were used for the analysis of Caco-2 samples. Like ESI, also APPI proved to be a reliable technique for the analysis of Caco-2 samples and even more flexible than ESI because of the wider dynamic linear range. The second part of the experimental study focused on metabolite profiling. Different mass spectrometric instruments and commercially available software tools were investigated for profiling metabolites in urine and hepatocyte samples. All the instruments tested (triple quadrupole, quadrupole time-of-flight, ion trap) exhibited some good and some bad features in searching for and identifying of expected and non-expected metabolites. Although, current profiling software is helpful, it is still insufficient. Thus a time-consuming largely manual approach is still required for metabolite profiling from complex biological matrices.

Expression of cytochrome P450 enzymes and metabolism of timolol in human ocular tissue

Glaucoma is a group of progressive optic neuropathies causing irreversible blindness if not diagnosed and treated in the early state of progression. Disease is often, but not always, associated with increased intraocular pressure (IOP), which is also the most important risk factor for glaucoma. Ophthlamic timolol preparations have been used for decades to lower increased intraocular pressure (IOP). Timolol is locally well tolerated but may cause e.g. cardiovascular and pulmonary adverse effects due to systemic absorption. It has been reported that approximately 80% of a topically administered eye drop is systemically absorbed. However, only limited information is available on timolol metabolism in the liver or especially in the human eye. The aim of this work was to investigate metabolism of timolol in human liver and human ocular tissues. The expression of drug metabolizing cytochrome P450 (CYP) enzymes in the human ciliary epithelial cells was studied. The metabolism of timolol and the interaction potential of timolol with other commercially available medicines were investigated in vitro using different liver preparations. The absorption of timolol to the aqueous humor from two commercially available products: 0.1% eye gel and 0.5% eye drops and the presence of timolol metabolites in the aqueous humor were investigated in a clinical trial. Timolol was confirmed to be metabolized mainly by CYP2D6 as previously suggested. Potent CYP2D6 inhibitors especially fluoxetine, paroxetine and quinidine inhibited the metabolism of timolol. The inhibition may be of clinical significance in patients using ophthalmic timolol products. CYP1A1 and CYP1B1 mRNAs were expressed in the human ciliary epithelial cells. CYP1B1 was also expressed at protein level and the expression was strongly induced by a known potent CYP1B1 inducer 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). The CYP1B1 induction is suggested to be mediated by aryl hydrocarbon receptor (AHR). Low levels of CYP2D6 mRNA splice variants were expressed in the human ciliary epithelial cells and very low levels of timolol metabolites were detected in the human aqueous humor. It seems that negligible amount of CYP2D6 protein is expressed in the human ocular tissues. Timolol 0.1% eye gel leads to aqueous humor concentration high enough to achieve therapeutic effect. Inter-individual variation in concentrations is low and intraocular as well as systemic safety can be increased when using this product with lower timolol concentration instead of timolol 0.5% eye drops.

Associations among Emdogain®, human oral carcinoma cells and oral proteolytic enzymes

The Enamel matrix derivative Emdogain® (EMD) is a commercially available tissue extract preparation of porcine enamel origin. Studies have shown EMD to be clinically useful in promoting periodontal regeneration. EMD has been widely used in periodontal therapy for over ten years, but the mechanism of its action and the exact composition are not completely clear. EMD is predominantly amelogenin (>90%). However, unlike amelogenin, EMD has a number of growth factor-like effects and it has been shown to enhance the proliferation, migration and other cellular functions of periodontal ligament fibroblasts and osteoblasts. In contrast, the effects of EMD on epithelial cell lines and in particular on oral malignant cells have not been adequately studied. In addition, EMD has effects on the production of cytokines by several oral cell lines and the product is in constant interaction with different oral enzymes. Regardless of the various unknown properties of EMD, it is said to be clinically safe in regenerative procedures, also in medically compromised patients. The aim of the study was to examine whether gingival crevicular fluid (GCF), which contains several different proteolysis enzymes, could degrade EMD and alter its biological functions. In addition, the objective was to study the effects of EMD on carcinogenesis-related factors, in particular MMPs, using in vitro and in vivo models. This study also aimed to contribute to the understanding of the composition of EMD. GCF was capable of degrading EMD, depending on the periodontal status, with markedly more degradation in all states of periodontal disease compared to healthy controls. EMD was observed to stimulate the migration of periodontal ligament fibroblasts (PLF), whereas EMD together with GCF could not stimulate this proliferation. In addition, recombinant amelogenin, the main component of EMD, decreased the migration of PLFs. A comparison of changes induced by EMD and TGF-β1 in the gene profiles of carcinoma cells showed TGF-β1 to regulate a greater number of genes than EMD. However, both of the study reagents enhanced the expression of MMP-10 and MMP-9. Furthermore, EMD was found to induce several factors closely related to carcinogenesis on gene, protein, cell and in vivo levels. EMD enhanced the production of MMP-2, MMP-9 and MMP-10 proteins by cultured carcinoma cells. In addition, EMD stimulated the migration and in vitro wound closure of carcinoma cells. EMD was also capable of promoting metastasis formation in mice. In conclusion, the diseased GCF, containing various proteases, causes degradation of EMD and decreased proliferation of PLFs. Thus, this in vitro study suggests that the regenerative effect of EMD may decrease due to proteases present in periodontal tissues during the inflammation and healing of the tissues in vivo. Furthermore, EMD was observed to enhance several carcinoma-related factors and in particular the production of MMPs by benign and malignant cell lines. These findings suggest that the clinical safety of EMD with regard to dysplastic mucosal lesions should be further investigated.

Defects in tricarboxylic acid cycle enzymes Fumarate hydratase and Succinate dehydrogenase in cancer

Hereditary leiomyomatosis and renal cell cancer (HLRCC) is a recently characterized cancer syndrome which predisposes to cutaneous and uterine leiomyomas as well as renal cell carcinoma (RCC). Uterine leiomyosarcoma (ULMS) has also been observed in certain Finnish HLRCC families. The predisposing gene for this syndrome, fumarate hydratase (FH), was identified in 2002. The well-known function of FH is in the tricarboxylic acid cycle (TCAC) in the energy metabolism of cells. As FH is a novel cancer gene, the role of FH mutations in tumours is in general unknown. Similarly, the mechanisms through which defective FH is associated with tumourigenesis are unclear. The loss of a wild type allele has been observed in virtually all HLRCC patients tumours and the FH enzyme activities are either totally lost or remarkably reduced in the tissues of mutation carrier patients. Therefore, FH is assumed to function as a tumour suppressor. Mutations in genes encoding subunits of other TCAC enzyme SDH have also been reported recently in tumours: mutations in SDHB, SDHC, and SDHD genes predispose to paraganglioma and pheochromocytoma. In the present study, mutations in the SDHB gene were observed to predispose to RCC. This was the first time that mutations in SDHB have been detected in extra-paraganglial tumours. Two different SDHB mutations were observed in two unrelated families. In the first family, the index patient was diagnosed with RCC at the age of 24 years. Additionally, his mother with a paraganglioma (PGL) of the heart and his maternal uncle with lung cancer were both carriers of the mutation. The RCC of the index patient and the PGL of his mother showed LOH. In the other family, an SDHB mutation was detected in two siblings who were both diagnosed with RCC at the ages of 24 and 26 years. Both of the siblings also suffered PGL. All these tumours showed LOH. Therefore, we concluded that mutations in SDHB predispose also for RCC in certain families. Several tumour types were analysed for FH mutations to define the role of FH mutations in these tumour types. In addition, patients with a putative cancer phenotype were analysed to identify new HLRCC families. Three FH variants were detected, of which two were novel. One of the variants was observed in a patient diagnosed with ULMS at the age of 41 years. However, LOH was not detected in the tumour tissue. The FH enzyme activity of the mutated protein was clearly reduced, being 43% of the activity of the normal protein. Together with the results from an earlier study we calculated that the prevalence of FH mutations in Finnish non-syndromic ULMS is around 2.4%. Therefore, FH mutations seem to have a minor role in the pathogenesis on non-syndromic ULMS. Two other germline variants were detected in a novel tumour type, ovarian mucinous cystadenoma. However, tumour tissues of the patients were not available for LOH studies and therefore LOH status remained unclear. Therefore, it is possible that FH mutations predispose also for ovarian tumours but further studies are needed to verify this result. A novel variant form of the FH gene (FHv) was identified and characterized in more detail. FHv contains an alternative first exon (1b), which appeared to function as 5 UTR sequence. The translation of FHv is initiated in vitro from exons two and three. The localization of FHv is both cytosolic and nuclear, in contrast to the localization of FH in mitochondria. FHv is expressed at low levels in all human tissues. Interestingly, the expression was induced after heat shock treatment and in chronic hypoxia. Therefore, FHv might have a role e.g. in the adaptation to unfavourable growth conditions. However, this remains to be elucidated.

Design and testing of stand-specific bucking instructions for use on modern cut-to-length harvesters

This study addresses three important issues in tree bucking optimization in the context of cut-to-length harvesting. (1) Would the fit between the log demand and log output distributions be better if the price and/or demand matrices controlling the bucking decisions on modern cut-to-length harvesters were adjusted to the unique conditions of each individual stand? (2) In what ways can we generate stand and product specific price and demand matrices? (3) What alternatives do we have to measure the fit between the log demand and log output distributions, and what would be an ideal goodness-of-fit measure? Three iterative search systems were developed for seeking stand-specific price and demand matrix sets: (1) A fuzzy logic control system for calibrating the price matrix of one log product for one stand at a time (the stand-level one-product approach); (2) a genetic algorithm system for adjusting the price matrices of one log product in parallel for several stands (the forest-level one-product approach); and (3) a genetic algorithm system for dividing the overall demand matrix of each of the several log products into stand-specific sub-demands simultaneously for several stands and products (the forest-level multi-product approach). The stem material used for testing the performance of the stand-specific price and demand matrices against that of the reference matrices was comprised of 9 155 Norway spruce (Picea abies (L.) Karst.) sawlog stems gathered by harvesters from 15 mature spruce-dominated stands in southern Finland. The reference price and demand matrices were either direct copies or slightly modified versions of those used by two Finnish sawmilling companies. Two types of stand-specific bucking matrices were compiled for each log product. One was from the harvester-collected stem profiles and the other was from the pre-harvest inventory data. Four goodness-of-fit measures were analyzed for their appropriateness in determining the similarity between the log demand and log output distributions: (1) the apportionment degree (index), (2) the chi-square statistic, (3) Laspeyres quantity index, and (4) the price-weighted apportionment degree. The study confirmed that any improvement in the fit between the log demand and log output distributions can only be realized at the expense of log volumes produced. Stand-level pre-control of price matrices was found to be advantageous, provided the control is done with perfect stem data. Forest-level pre-control of price matrices resulted in no improvement in the cumulative apportionment degree. Cutting stands under the control of stand-specific demand matrices yielded a better total fit between the demand and output matrices at the forest level than was obtained by cutting each stand with non-stand-specific reference matrices. The theoretical and experimental analyses suggest that none of the three alternative goodness-of-fit measures clearly outperforms the traditional apportionment degree measure. Keywords: harvesting, tree bucking optimization, simulation, fuzzy control, genetic algorithms, goodness-of-fit

Protein cross-linking with oxidative enzymes and transglutaminase : Effects in meat protein systems

The effects of tyrosinase, laccase and transglutaminase (TG) were studied in different meat protein systems. The study was focused on the effects of the enzymes on the gel formation properties of myofibrils, and on the textural and water-holding properties of the heated meat systems. The cross-linking efficiency of a novel Trichoderma reesei tyrosinase was compared to that of the commercial Agaricus bisporus tyrosinase. Trichoderma tyrosinase was found to be superior compared to the Agaricus enzyme in its protein cross-linking efficiency and in the incorporation of a small molecule into a complex proteinaceous substrate. Tyrosinase, laccase and TG all polymerised myofibrillar proteins, but laccase was also found to cause protein fragmentation. A positive connection between covalent cross-link and gel formation was observed with tyrosinase and TG. Laccase was able to increase the gel formation only slightly. With an excessive laccase dosage the gel formation declined due to protein fragmentation. Tyrosinase, laccase and TG had different effects on the texture and water-holding of the heated chicken breast meat homogenates. Tyrosinase improved the firmness of the homogenate gels free of phosphate and with a low amount of meat. TG improved the firmness of all studied homogenates. Laccase weakened the gel firmness of the low-meat, low-salt and low-salt/phosphate homogenates and maintained the firmness on the control level in the homogenate free of phosphate. Tyrosinase was the only enzyme capable of reducing the weight loss in the homogenates containing a low amount of meat and a low amount of NaCl. TG was the only enzyme that could positively affect the firmness of the homogenate gel containing both low NaCl and phosphate amounts. In pilot scale the test products were made of coarsely ground chicken breast fillet with a moderate amount of salt. Increasining the amount of meat, salt and TG contents favoured the development of firmness of the test products. The evaporation loss decreased slightly along with increasing TG and NaCl amounts in the experimental conditions used, indicating a positive interaction between these two factors. In this work it was shown that tyrosinase, laccase and TG affected the same myofibrillar proteins, i.e. myosin and troponin T. However, these enzymes had distinguishable effects on the gel formation of a myofibril system as well as on the textural and water-holding properties of the finely ground meat homogenates, reflecting distinctions at least in the reaction mechanisms and target amino acid availability in the protein substrates for these enzymes.

Characterization of lignin-modifying enzymes of the selective white-rot fungus Physisporinus rivulosus

White-rot fungi are wood degrading organisms that are able to decompose all wood polymers; lignin, cellulose and hemicellulose. Especially the selective white-rot fungi that decompose preferentially wood lignin are promising for biopulping applications. In biopulping the pretreatment of wood chips with white-rot fungi enhances the subsequent pulping step and substantially reduces the refining energy consumption in mechanical pulping. Because it is not possible to carry out biopulping in industrial scale as a closed process it has been necessary to search for new selective strains of white-rot fungi which naturally occur in Finland and cause selective white-rot of Finnish wood raw-material. In a screening of 300 fungal strains a rare polypore, Physisporinus rivulosus strain T241i isolated from a forest burn research site, was found to be a selective lignin degrader and promising for the use in biopulping. Since selective lignin degradation is apparently essential for biopulping, knowledge on lignin-modifying enzymes and the regulation of their production by a biopulping fungus is needed. White-rot fungal enzymes that participate in lignin degradation are laccase, lignin peroxidase (LiP), manganese peroxidase (MnP), versatile peroxidase (VP) and hydrogen peroxide forming enzymes. In this study, P. rivulosus was observed to produce MnP, laccase and oxalic acid during growth on wood chips. In liquid cultures manganese and veratryl alcohol increased the production of acidic MnP isoforms detected also in wood chip cultures. Laccase production by P. rivulosus was low unless the cultures were supplemented with sawdust and charred wood, the components of natural growth environment of the fungus. In white-rot fungi the lignin-modifying enzymes are typically present as multiple isoforms. In this study, two MnP encoding genes, mnpA and mnpB, were cloned and characterized from P. rivulosus T241i. Analysis of the N-terminal amino acid sequences of two purified MnPs and putative amino acid sequence of the two cloned mnp genes suggested that P. rivulosus possesses at least four mnp genes. The genes mnpA and mnpB markedly differ from each other by the gene length, sequence and intron-exon structure. In addition, their expression is differentially affected by the addition of manganese and veratryl alcohol. P. rivulosus produced laccase as at least two isoforms. The results of this study revealed that the production of MnP and laccase was differentially regulated in P. rivulosus, which ensures the efficient lignin degradation under a variety of environmental conditions.

Molecular Modelling of Estrogen-Producing Enzymes CYP450 Aromatase and 17beta-Hydroxysteroid Dehydrogenase Type 1

Breast cancer is the most common cancer in women in the western countries. Approximately two-thirds of breast cancer tumours are hormone dependent, requiring estrogens to grow. Estrogens are formed in the human body via a multistep route starting from cholesterol. The final steps in the biosynthesis include the CYP450 aromatase enzyme, converting the male hormones androgens (preferred substrate androstenedione ASD) into estrogens(estrone E1), and the 17beta-HSD1 enzyme, converting the biologically less active E1 into the active hormone 17beta-hydroxyestradiol E2. E2 is bound to the nuclear estrogen receptors causing a cascade of biochemical reactions leading to cell proliferation in normal tissue, and to tumour growth in cancer tissue. Aromatase and 17beta-HSD1 are expressed in or near the breast tumour, locally providing the tissue with estrogens. One approach in treating hormone dependent breast tumours is to block the local estrogen production by inhibiting these two enzymes. Aromatase inhibitors are already on the market in treating breast cancer, despite the lack of an experimentally solved structure. The structure of 17beta-HSD1, on the other hand, has been solved, but no commercial drugs have emerged from the drug discovery projects reported in the literature. Computer-assisted molecular modelling is an invaluable tool in modern drug design projects. Modelling techniques can be used to generate a model of the target protein and to design novel inhibitors for them even if the target protein structure is unknown. Molecular modelling has applications in predicting the activities of theoretical inhibitors and in finding possible active inhibitors from a compound database. Inhibitor binding at atomic level can also be studied with molecular modelling. To clarify the interactions between the aromatase enzyme and its substrate and inhibitors, we generated a homology model based on a mammalian CYP450 enzyme, rabbit progesterone 21-hydroxylase CYP2C5. The model was carefully validated using molecular dynamics simulations (MDS) with and without the natural substrate ASD. Binding orientation of the inhibitors was based on the hypothesis that the inhibitors coordinate to the heme iron, and were studied using MDS. The inhibitors were dietary phytoestrogens, which have been shown to reduce the risk for breast cancer. To further validate the model, the interactions of a commercial breast cancer drug were studied with MDS and ligand–protein docking. In the case of 17beta-HSD1, a 3D QSAR model was generated on the basis of MDS of an enzyme complex with active inhibitor and ligand–protein docking, employing a compound library synthesised in our laboratory. Furthermore, four pharmacophore hypotheses with and without a bound substrate or an inhibitor were developed and used in screening a commercial database of drug-like compounds. The homology model of aromatase showed stable behaviour in MDS and was capable of explaining most of the results from mutagenesis studies. We were able to identify the active site residues contributing to the inhibitor binding, and explain differences in coordination geometry corresponding to the inhibitory activity. Interactions between the inhibitors and aromatase were in agreement with the mutagenesis studies reported for aromatase. Simulations of 17beta-HSD1 with inhibitors revealed an inhibitor binding mode with hydrogen bond interactions previously not reported, and a hydrophobic pocket capable of accommodating a bulky side chain. Pharmacophore hypothesis generation, followed by virtual screening, was able to identify several compounds that can be used in lead compound generation. The visualisation of the interaction fields from the QSAR model and the pharmacophores provided us with novel ideas for inhibitor development in our drug discovery project.