Lignin biosynthesis in Norway spruce: from a model system to the tree

Lignin is a hydrophobic polymer that is synthesised in the secondary cell walls of all vascular plants. It enables water conduction through the stem, supports the upright growth habit and protects against invading pathogens. In addition, lignin hinders the utilisation of the cellulosic cell walls of plants in pulp and paper industry and as forage. Lignin precursors are synthesised in the cytoplasm through the phenylpropanoid pathway, transported into the cell wall and oxidised by peroxidases or laccases to phenoxy radicals that couple to form the lignin polymer. This study was conducted to characterise the lignin biosynthetic pathway in Norway spruce (Picea abies (L.) Karst.). We focused on the less well-known polymerisation stage, to identify the enzymes and the regulatory mechanisms that are involved. Available data for lignin biosynthesis in gymnosperms is scarce and, for example, the latest improvements in precursor biosynthesis have only been verified in herbaceous plants. Therefore, we also wanted to study in detail the roles of individual gene family members during developmental and stress-induced lignification, using EST sequencing and real-time RT-PCR. We used, as a model, a Norway spruce tissue culture line that produces extracellular lignin into the culture medium, and showed that lignin polymerisation in the tissue culture depends on peroxidase activity. We identified in the culture medium a significant NADH oxidase activity that could generate H2O2 for peroxidases. Two basic culture medium peroxidases were shown to have high affinity to coniferyl alcohol. Conservation of the putative substrate-binding amino acids was observed when the spruce peroxidase sequences were compared with other peroxidases with high affinity to coniferyl alcohol. We also used different peroxidase fractions to produce synthetic in vitro lignins from coniferyl alcohol; however, the linkage pattern of the suspension culture lignin could not be reproduced in vitro with the purified peroxidases, nor with the full complement of culture medium proteins. This emphasised the importance of the precursor radical concentration in the reaction zone, which is controlled by the cells through the secretion of both the lignin precursors and the oxidative enzymes to the apoplast. In addition, we identified basic peroxidases that were reversibly bound to the lignin precipitate. They could be involved, for example, in the oxidation of polymeric lignin, which is required for polymer growth. The dibenzodioxocin substructure was used as a marker for polymer oxidation in the in vitro polymerisation studies, as it is a typical substructure in wood lignin and in the suspension culture lignin. Using immunolocalisation, we found the structure mainly in the S2+S3 layers of the secondary cell walls of Norway spruce tracheids. The structure was primarily formed during the late phases of lignification. Contrary to the earlier assumptions, it appears to be a terminal structure in the lignin macromolecule. Most lignin biosynthetic enzymes are encoded for by several genes, all of which may not participate in lignin biosynthesis. In order to identify the gene family members that are responsible for developmental lignification, ESTs were sequenced from the lignin-forming tissue culture and developing xylem of spruce. Expression of the identified lignin biosynthetic genes was studied using real-time RT-PCR. Candidate genes for developmental lignification were identified by a coordinated, high expression of certain genes within the gene families in all lignin-forming tissues. However, such coordinated expression was not found for peroxidase genes. We also studied stress-induced lignification either during compression wood formation by bending the stems or after Heterobasidion annosum infection. Based on gene expression profiles, stress-induced monolignol biosynthesis appeared similar to the developmental process, and only single PAL and C3H genes were specifically up-regulated by stress. On the contrary, the up-regulated peroxidase genes differed between developmental and stress-induced lignification, indicating specific responses.

Cholesterol Disturbances in Inherited Neurodegenerative Diseases : Studies on Npc1, Ppt1 and Ctsd deficient mice

The central nervous system (CNS) is the most cholesterol-rich organ in the body. Cholesterol is essential to CNS functions such as synaptogenesis and formation of myelin. Significant differences exist in cholesterol metabolism between the CNS and the peripheral organs. However, the regulation of cholesterol metabolism in the CNS is poorly understood compared to our knowledge of the regulation of cholesterol homeostasis in organs reached by cholesterol-carrying lipoprotein particles in the circulation. Defects in CNS cholesterol homeostasis have been linked to a variety of neurodegenerative diseases, including common diseases with complex pathogenetic mechanisms such as Alzheimer s disease. In spite of intense effort, the mechanisms which link disturbed cholesterol homeostasis to these diseases remain elusive. We used three inherited recessive neurodegenerative disorders as models in the studies included in this thesis: Niemann-Pick type C (NPC), infantile neuronal ceroid lipofuscinosis and cathepsin D deficiency. Of these three, NPC has previously been linked to disturbed intracellular cholesterol metabolism. Elucidating the mechanisms with which disturbances of cholesterol homeostasis link to neurodegeneration in recessive inherited disorders with known genetic lesions should shed light on how cholesterol is handled in the healthy CNS and help to understand how these and more complex diseases develop. In the first study we analyzed the synthesis of sterols and the assembly and secretion of lipoprotein particles in Npc1 deficient primary astrocytes. We found that both wild type and Npc1 deficient astrocytes retain significant amounts of desmosterol and other cholesterol precursor sterols as membrane constituents. No difference was observed in the synthesis of sterols and the secretion of newly synthesized sterols between Npc1 wild type, heterozygote or knockout astrocytes. We found that the incorporation of newly synthesized sterols into secreted lipoprotein particles was not inhibited by Npc1 mutation, and the lipoprotein particles were similar to those excreted by wild type astrocytes in shape and size. The bulk of cholesterol was found to be secreted independently of secreted NPC2. These observations demonstrate the ability of Npc1 deficient astrocytes to handle de novo sterols, and highlight the unique sterol composition in the developing brain. Infantile neuronal ceroid lipofuscinosis is caused by the deficiency of a functional Ppt1 enzyme in the cells. In the second study, global gene expression studies of approximately 14000 mouse genes showed significant changes in the expression of 135 genes in Ppt1 deficient neurons compared to wild type. Several genes encoding for enzymes of the mevalonate pathway of cholesterol biosynthesis showed increased expression. As predicted by the expression data, sterol biosynthesis was found to be upregulated in the knockout neurons. These data link Ppt1 deficiency to disturbed cholesterol metabolism in CNS neurons. In the third study we investigated the effect of cathepsin D deficiency on the structure of myelin and lipid homeostasis in the brain. Our proteomics data, immunohistochemistry and western blotting data showed altered levels of the myelin protein components myelin basic protein, proteolipid protein and 2 , 3 -cyclic nucleotide 3 phosphodiesterase in the brains of cathepsin D deficient mice. Electron microscopy revealed altered myelin structure in cathepsin D deficient brains. Additionally, plasmalogen-derived alkenyl chains and 20- and 24-carbon saturated and monounsaturated fatty acids typical for glycosphingolipids were found to be significantly reduced, but polyunsaturated species were significantly increased in the knockout brains, pointing to a decrease in white matter. The levels of ApoE and ABCA1 proteins linked to cholesterol efflux in the CNS were found to be altered in the brains of cathepsin D deficient mice, along with an accumulation of cholesteryl esters and a decrease in triglycerols. Together these data demonstrate altered myelin architecture in cathepsin D deficient mice and link cathepsin D deficiency to aberrant cholesterol metabolism and trafficking. Basic research into rare monogenic diseases sheds light on the underlying biological processes which are perturbed in these conditions and contributes to our understanding of the physiological function of healthy cells. Eventually, understanding gained from the study of disease models may contribute towards establishing treatment for these disorders and further our understanding of the pathogenesis of other, more complex and common diseases.

Antizyme Inhibitor 2 : A Novel Regulator of Cellular Polyamine Homeostasis

Polyamines are organic polycations that participate in various physiological functions, including cell proliferation, differentiation and apoptosis. Cellular polyamines originate from endogenous biosynthesis and exogenous sources. Their subcellular pool is under strict control, achieved by regulating their uptake and metabolism. Polyamine-induced proteins called antizymes (AZ) act as key regulators of intracellular polyamine concentration. They regulate both the transport of polyamines and the activity and degradation of ornithine decarboxylase (ODC), the rate-limiting enzyme in polyamine biosynthesis. AZs themselves are negatively regulated by antizyme inhibitor (AZIN). AZIN functions as a positive regulator of cellular polyamine homeostasis, which by binding to AZs reactivates ODC and induces the uptake of polyamines. In various pathological conditions, including cancer, polyamine levels are misregulated. Polyamine homeostasis has therefore become an attractive target for therapeutic interventions and it is thus crucial to characterize the molecular basis underlying the homeostatic regulation. A novel human AZIN-resembling protein was previously identified in our group. The purpose of this study was to elucidate the function and distribution of this protein, termed as an antizyme inhibitor 2 (AZIN2). According to my results, AZIN2 functions as a novel regulator of polyamine homeostasis. It shows no enzymatic activity, but instead it binds AZs and negates their activity, which subsequently leads to reactivation of ODC and inhibition of its degradation. Expression of AZIN2 is restricted to terminally differentiated cells, such as mast cells (MC) and neurosecretory cells. In these actively secreting cell types, AZIN2 localizes to subcellular vesicles or granules where its function is important for the vesicle-mediated secretion. In MCs, AZIN2 localizes to the serotonin-containing subset of MC granules, and its expression is coupled to MC activation. The functional role of polyamines as potential mediators of MC activity was also investigated, and it was observed that the secretion of serotonin is selectively dependent on activation of ODC. In neurosecretory cells, AZIN2-positive vesicles localize mainly to the trans-Golgi network (TGN). Depletion of AZIN2 or cellular polyamines causes selective fragmentation of the TGN and retards secretion of proteins. Since addition of exogenous polyamines reverses these effects, the data indicate that AZIN2 and its downstream effectors, polyamines, are functionally implicated in the regulation of secretory vesicle transport. My studies therefore reveal a novel function for polyamines as modulators of both constitutive and regulated secretion. Based on the results, I propose that the role of AZIN2 is to act as a local in situ activator of polyamine biosynthesis.

Interactions between transgenic trees and mycorrhizal and pathogenic fungi

The development of biotechnology techniques in plant breeding and the new commercial applications have raised public and scientific concerns about the safety of genetically modified (GM) crops and trees. To find out the feasibility of these new technologies in the breeding of commercially important Finnish hardwood species and to estimate the ecological risks of the produced transgenic plants, the experiments of this study have been conducted as a part of a larger project focusing on the risk assessment of GM-trees. Transgenic Betula pendula and Populus trees were produced via Agrobacterium mediated transformation. Stilbene synthase (STS) gene from pine (Pinus sylvestris) and chitinase gene from sugar beet (Beta vulgaris) were transferred to (hybrid) aspen and birch, respectively, to improve disease resistance against fungal pathogens. To modify lignin biosynthesis, a 4-coumarate:coenzyme A ligase (4CL) gene fragment in antisense orientation was introduced into two birch clones. In in vitro test, one transgenic aspen line expressing pine STS gene showed increased resistance to decay fungus Phellinus tremulae. In the field, chitinase transgenic birch lines were more susceptible to leaf spot (Pyrenopeziza betulicola) than the non-transgenic control clone while the resistance against birch rust (Melampsoridium betulinum) was improved. No changes in the content or composition of lignin were detected in the 4CL antisense birch lines. In order to evaluate the ecological effects of the produced GM trees on non-target organisms, an in vitro mycorrhiza experiment with Paxillus involutus and a decomposition experiment in the field were performed. The expression of a transgenic chitinase did not disturb the establishment of mycorrhizal symbiosis between birch and P. involutus in vitro. 4CL antisense transformed birch lines showed retarded root growth but were able to form normal ectomycorrhizal associations with the mycorrhizal fungus in vitro. 4CL lines also showed normal litter decomposition. Unexpected growth reductions resulting from the gene transformation were observed in chitinase transgenic and 4CL antisense birch lines. These results indicate that genetic engineering can provide a tool in increasing disease resistance in Finnish tree species. More extensive data with several ectomycorrhizal species is needed to evaluate the consequences of transgene expression on beneficial plant-fungus symbioses. The potential pleiotropic effects of the transgene should also be taken into account when considering the safety of transgenic trees.

Regulation of strawberry growth and development

Strawberries (Fragaria sp.) are adapted to diverse environmental conditions from the tropics to about 70ºN, so different responses to environmental conditions can be found. Most genotypes of garden strawberry (F. x ananassa Duch.) and woodland strawberry (F. vesca L.) are short-day (SD) plants that are induced to flowering by photoperiods under a critical limit, but also various photoperiod x temperature interactions can be found. In addition, continuously flowering everbearing (EB) genotypes are found. In addition to flowering, axillary bud differentiation in strawberry is regulated by photoperiod. In SD conditions, axillary buds differentiate to rosette-like structures called "branch crowns", whereas in long-day conditions (LD) they form runners, branches with 2 long internodes followed by a daughter plant (leaf rosette). The number of crown branches determines the yield of the plant, since inflorescences are formed from the apical meristems of the crown. Although axillary bud differentiation is an important developmental process in strawberries, its environmental and hormonal regulation has not been characterized in detail. Moreover, the genetic mechanisms underlying axillary bud differentiation and regulation of flowering time in these species are almost completely unresolved. These topics have been studied in this thesis in order to enhance strawberry research, cultivation and breeding. The results showed that 8-12 SD cycles suppressed runner initiation from the axillary buds of the garden strawberry cv. Korona with the concomitant induction of crown branching, and 3 weeks of SD was sufficient for the induction of flowering in the main crown. Furthermore, a second SD treatment given a few weeks after the first SD period can be used to induce flowering in the primary branch crowns and to induce the formation of secondary branches. Thus, artificial SD treatments effectively stimulate crown branching, providing one means for the increase of cropping (yield) potential in strawberry. It was also shown by growth regulation applications, quantitave hormone analysis and gene expression analysis that gibberellin (GA) is one of the key signals involved in the photoperiod control of shoot differentiation. The results indicate that photoperiod controls GA activity specifically in axillary buds, thereby determining bud fate. It was further shown that chemical control of GA biosynthesis by prohexadione-calcium can be utilized to prevent excessive runner formation and induce crown branching in strawberry fields. Moreover, ProCa increased berry yield up to 50%, showing that it is an easier and more applicable alternative to artificial SD treatments for controlling strawberry crown development and yield. Finally, flowering gene pathways in Fragaria were explored by searching for homologs of 118 Arabidopsis thaliana flowering-time genes. In total, 66 gene homologs were identified, and they distributed to all known flowering pathways, suggesting the presence of these pathways also in strawberry. Expression analysis of selected genes revealed that the mRNA of putative floral identity gene APETALA1 accumulated in the shoot apex of the EB genotype after the induction of flowering, whereas it was absent in vegetative SD genotype, indicating the usefulness of this gene product as the marker of floral initiation. The present data enables the further exploration of strawberry flowering pathways with genetic transformation, gene mapping and transcriptomics methods.

The effects of environmental factors on biomass and microcystin production by the freshwater cyanobacterial genera Microcystis and Anabaena

Several cyanobacterial genera produce the hepatotoxins, microcystins. Microcystins are produced only in cells that have microcystin synthetase gene (mcy) clusters, which encode enzyme complexes involved in microcystin biosynthesis. Microcystin-producing and nonmicrocystin-producing genotypes of single cyanobacterial genus may occur simultaneously in situ. Previously, the effects of environmental factors on the growth and microcystin production of cyanobacteria have mainly been studied by means of isolated cyanobacteria cultures in the laboratory. Studies in the field have been difficult, owing to the lack of methods to identify and quantify the different genotypes. In this study, genus-specific microcystin synthetase E (mcyE) gene primers were designed and a method to identify and quantify the mcyE copy numbers was developed and used in situ. Microcystis and Anabaena mcyE genes were observed in two Finnish lakes. Microcystis appeared to be the most abundant microcystin producer in Lake Tuusulanjärvi and in one basin of Lake Hiidenvesi. Because the most potent microcystin-producing genus of a lake can be identified, it will be possible in the future to design genus-targeted strategies for lake restoration. Effects of P and N concentrations on the biomass of microcystin-producing and nonmicrocystin-producing Microcystis strains and an Anabaena strain were studied in cultures. P and N concentrations and their combined effect increased cyanobacterial biomass of all Microcystis strains. The biomass of microcystin-producing Microcystis was higher than that of nonmicrocystin-producing strains at high nutrient concentrations. The P concentration increased Anabaena biomass, but the effect of N concentration was statistically insignificant for growth yield, probably due to the ability of the genus to fix molecular N2. P and N concentrations and combined nutrients caused an increase in cellular microcystin concentrations of the Microcystis strain cultivated in chemostat cultures. Cyanobacteria are able to hydrolyse nutrients from organic matter through extracellular enzyme activities. Leucine aminopeptidase (LAP) activity was observed in an axenic N2-fixing Anabaena strain grown in batch cultures. The P concentration caused a statistically significant increase in LAP activity, whereas the effect of N concentration was insignificant. The highest LAP activities were observed in the most eutrophic basins of Lake Hiidenvesi. LAP activity probably originated mostly from attached heterotrophic bacteria and less from cyanobacteria.

Lantibiotic Nisin and Its Detection Methods

The type A lantibiotic nisin produced by several Lactococcus lactis strains, and one Streptococcus uberis strainis a small antimicrobial peptide that inhibits the growth of a wide range of gram-positive bacteria, such as Bacillus, Clostridium, Listeria and Staphylococcus species. It is nontoxic to humans and used as a food preservative (E234) in more than 50 countries including the EU, the USA, and China. National legislations concerning maximum addition levels of nisin in different foods vary greatly. Therefore, there is a demand for non-laborious and sensitive methods to identify and quantify nisin reliably from different food matrices. The horizontal inhibition assay, based on the inhibitory effect of nisin to Micrococcus luteus is the base for most quantification methods developed so far. However, the sensitivity and accuracy of the agar diffusion method is affected by several parameters. Immunological tests have also been described. Taken into account the sensitivity of immunological methods to interfering substances within sample matrices, and possible cross-reactivities with lantibiotics structurally close to nisin, their usefulness for nisin detection from food samples remains limited. The proteins responsible for nisin biosynthesis, and producer self-immunity are encoded by genes arranged into two inducible operons, nisA/Z/QBTCIPRK and nisFEG, which also contain internal, constitutive promoters PnisI and PnisR. The transmembrane histidine kinase NisK and the response regulator NisR form a two-component signal transduction system, in which NisK autophosphorylates after exposure to extra cellular nisin, and subsequently transfers the phosphate to NisR. The phosphorylated NisR then relays the signal downstream by binding to two regulated promoters in the nisin gene cluster, i.e the nisA/Z/Qand the nisF promoters, thus activating transcription of the structural gene nisA/Z/Q and the downstream genes nisBTCIPRK from the nisA/Z/Q promoter, and the genes nisFEG from the nisF promoter. In this work two novel and highly sensitive nisin bioassays were developed. Both of these quantification methods were based on NisRK mediated, nisin induced Green Fluorescent Protein (GFP) fluorescence. The suitabilities of these assays for quantifica¬tion of nisin from food samples were evaluated in several food matrices. These bioassays had nisin sensitivities in the nanogram or picogram levels. In addition, shelf life of nisin in cooked sausages and retainment of the induction activity of nisin in intestinal chyme (intestinal content) was assessed.

Characterization of Actinobacteria Degrading and Tolerating Organic Pollutants

Species of the genera Rhodococcus, Gordonia and Mycobacterium are known as degraders of recalcitrant pollutants. These bacteria are good survivors in harsh environments. Due to such properties these organisms are able to occupy a wide range of environmental niches. The members of these taxa have been suggested as tools for biotechnical applications such as bioremediation and biosynthesis. At the same time several of the species are known as opportunistic human pathogens. Therefore, the detailed characterization of any isolate that has potential for biotechnological applications is very important. This thesis deals with several corynebacterial strains originating from different polluted environments: soil, water-damaged indoor walls, and drinking water distribution systems. A polyphasic taxonomic approach was applied for characterization of the isolates. We found that the strains degrading monoaromatic compounds belonged to Rhodococcus opacus, a species that has not been associated with any health problem. The taxonomic position of strain B293, used for many years in degradation research under different names, was clarified. We assigned it to the species Gordonia polyisoprenivorans. This species is classified under European Biohazard grouping 1, meaning that it is not considered a health hazard for humans. However, there are reports of catheter-associated bacteraemia caused by G. polyisoprenivorans. Our results suggested that the ability of the organism to grow on phthalate esters, used as softeners in medical plastics, may be associated with the colonization of catheters and other devices. In this thesis Mycobacterium lentiflavum, a new emerging opportunistic human pathogen, was isolated from biofilms growing in public drinking water distribution systems. Our report on isolation of M. lentiflavum from water supplies is the second report on this species from drinking water systems, which may thus constitute a reservoir of M. lentiflavum. Automated riboprinting was evaluated for its applicability in rapidly identifying environmental mycobacteria. The technique was found useful in the characterization of several species of rapidly and slowly growing environmental mycobacteria. The second aspect of this thesis refers to characterization of the degradation and tolerance power of several R. opacus, M. murale and G. polyisoprenivorans strains. R. opacus GM-14 utilizes a wide range of aromatic substrates, including benzene, 15 different halobenzenes, 18 phenols and 7 benzoates. This study revealed the high tolerance of R. opacus strains toward toxic hydrophobic compounds. R. opacus GM-14 grew in mineral medium to which benzene or monochlorobenzene was added in amounts of 13 or 3 g l-1, respectively. R. opacus GM-29 utilized toluene and benzene for growth. Strain GM-29 grew in mineral medium with 7 g l-1 of liquid toluene or benzene as the sole carbon source, corresponding to aqueous concentrations of 470 and 650 mg l-1, respectively. Most organic solvents, such as toluene and benzene, due to their high level of hydrophobicity, pass through the bacterial membrane, causing its disintegration. In this thesis the mechanisms of adaptation of rhodococci to toxic hydrophobic compounds were investigated. The rhodococcal strains increased the level of saturation of their cellular fatty acids in response to challenge with phenol, chlorophenol, benzene, chlorobenzene or toluene. The results indicated that increase in the saturation level of cellular fatty acids, particularly that in tuberculostearic acid, is part of the adaptation mechanism of strains GM-14 and GM-29 to the presence of toxic hydrophobic compounds.

Synthesis, Reactivity and Biological Activity of 17β-HSD1 Inhibitors Based on a Thieno[2,3-d]pryrimidin-4(3H)-one Core

Breast cancer is the most common cancer in women in Western countries. In the early stages of development most breast cancers are hormone-dependent, and estrogens, especially estradiol, have a pivotal role in their development and progression. One approach to the treatment of hormone-dependent breast cancers is to block the formation of the active estrogens by inhibiting the action of the steroid metabolising enzymes. 17beta-Hydroxysteroid dehydrogenase type 1 (17beta-HSD1) is a key enzyme in the biosynthesis of estradiol, the most potent female sex hormone. The 17beta-HSD1 enzyme catalyses the final step and converts estrone into the biologically active estradiol. Blocking 17beta-HSD1 activity with a specific enzyme inhibitor could provide a means to reduce circulating and tumour estradiol levels and thus promote tumour regression. In recent years 17beta-HSD1 has been recognised as an important drug target. Some inhibitors of 17beta-HSD1 have been reported, however, there are no inhibitors on the market nor have clinical trials been announced. The majority of known 17beta-HSD1 inhibitors are based on steroidal structures, while relatively little has been reported on non-steroidal inhibitors. As compared with 17beta-HSD1 inhibitors based on steroidal structures, non-steroidal compounds could have advantages of synthetic accessibility, drug-likeness, selectivity and non-estrogenicity. This study describes the synthesis of large group of novel 17beta-HSD1 inhibitors based on a non-steroidal thieno[2,3-d]pyrimidin-4(3H)-one core. An efficient synthesis route was developed for the lead compound and subsequently employed in the synthesis of thieno[2,3-d]pyrimidin-4(3H)-one based molecule library. The biological activities and binding of these inhibitors to 17beta-HSD1 and, finally, the quantitative structure activity relationship (QSAR) model are also reported. In this study, several potent and selective 17beta-HSD1 inhibitors without estrogenic activity were identified. This establishment of a novel class of inhibitors is a progressive achievement in 17beta-HSD1 inhibitor development. Furthermore, the 3D-QSAR model, constructed on the basis of this study, offers a powerful tool for future 17beta-HSD1 inhibitor development. As part of the fundamental science underpinning this research, the chemical reactivity of fused (di)cycloalkeno thieno[2,3-d]pyrimidin-4(3H)-ones with electrophilic reagents, i.e. Vilsmeier reagent and dimethylformamide dimethylacetal, was investigated. These findings resulted in a revision of the reaction mechanism of Vilsmeier haloformylation and further contributed to understanding the chemical reactivity of this compound class. This study revealed that the reactivity is dependent upon a stereoelectronic effect arising from different ring conformations.

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.