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.

On the environmental genetics of depressive symptoms : Focus on two serotonergic genes

Depression is a complex psychiatric disorder influenced by several genes, environmental factors, and their interplay. Serotonin receptor 2A (HTR2A) and tryptophan hydroxylase 1 (TPH1) genes have been implicated in vulnerability to depression and other psychiatric disorders, but the results have been inconsistent. The present study examined whether these two genes moderated the influence of different depressogenic environmental factors on subthreshold depressive symptoms (assessed on a modified version of Beck s Depression Inventory, BDI) and depression-related temperament, i.e., harm avoidance (assessed on the Temperament and Character Inventory, TCI). The environmental factors included measures of childhood and adolescence exposure, i.e., maternal nurturance and parental socioeconomic status, and adulthood social circumstances, i.e., perceived social support and urban/rural residence. The participants were two randomly selected subsamples (n = 1246, n = 341) from the longitudinal population-based Cardiovascular Risk in Young Finns study (n = 3596). Childhood environmental factors were assessed when the participants were 3 to 18 years of age, and three years after the baseline. Adulthood environmental factors and outcome measures were assessed 17 and 21 years later when the participants were 21 to 39 years of age. The T102C polymorphism of the HTR2A gene moderated the association between childhood maternal nurturance and adulthood depressive symptoms, such that exposure to high maternal nurturance predicted low depressive symptoms among individuals carrying the T/T or T/C genotypes, but not among those carrying the C/C genotype. Likewise, high parental SES predicted low adulthood harm avoidance in individuals carrying the T/T or T/C genotype, but not in C/C-genotype carriers. Individuals carrying the T/T or T/C genotype were also sensitive to urban/rural residence, such that they had lower depressive symptoms in urban than in rural areas, whereas those carrying the C/C genotype were not sensitive to urban/rural residence difference. HTR2A did not moderate the influence of social support. TheA779C/A218C haplotype of the TPH1 gene was not involved in the association between childhood environment and adulthood outcomes. However, individuals carrying A alleles of the TPH1 haplotype were more vulnerable to the lack of adulthood social support in terms of high depressive symptoms than their counterparts carrying no A alleles. Furthermore, individuals living in remote rural areas and carrying the A/A haplotype had higher depressive symptoms than those carrying other genotypes of the TPH1. The findings suggest that the HTR2A and TPH1 genes may be involved in the development of depression by influencing individual s sensitivity to depressogenic environmental influences.

Circulating Glucocorticoid Bioactivity and Serum Glucocorticoid-Responsive Biomarkers during Steroid Therapy in Children and Adolescents with Inflammatory Bowel Disease

Objective: Glucocorticoid therapy is used worldwide to treat various inflammatory and immune conditions, including inflammatory bowel disease (IBD). In IBD, 80% of the patients obtain a positive response to the therapy; however the development of glucocorticoid-related side-effects is common. Our aim was therefore to study the possibility of optimizing glucocorticoid therapy in children and adolescents with IBD by measuring circulating glucocorticoid bioactivity (GBA) and serum glucocorticoid-responsive biomarkers in patients receiving steroid treatment for active disease. Methods: A total of sixty-nine paediatric IBD patients from the Paediatric Outpatient Clinics of the University Hospitals of Helsinki and Tampere participated in the studies. Control patients included 101 non-IBD patients and 41 disease controls in remission. In patients with active disease, blood samples were withdrawn before the glucocorticoid therapy was started, at 2-4 weeks after the initiation of the steroid and at 1-month intervals thereafter. Clinical response to glucocorticoid treatment and the development of steroid adverse events was carefully registered. GBA was analyzed with a COS-1 cell bioassay. The measured glucocorticoid therapy-responsive biomarkers included adipocyte-derived adiponectin and leptin, bone turnover-related collagen markers amino-terminal type I procollagen propeptide (PINP) and carboxyterminal telopeptide of type I collagen (ICTP) as well as insulin-like growth factor 1 (IGF-1) and sex hormone-binding globulin (SHBG), and inflammatory marker high-sensitivity C-reactive protein (hs-CRP). Results: The most promising marker for glucocorticoid sensitivity was serum adiponectin that associated with steroid therapy–related adverse events. Serum leptin indicated a similar trend. In contrast, circulating GBA rose in all subjects receiving glucocorticoid treatment but did not associate with the clinical response to steroids or with glucocorticoid therapy-related side-effects. Of notice, young patients (<10 years) showed similar GBA levels than older patients, despite receiving higher weight-adjusted doses of glucocorticoid. Markers of bone formation were lower in children with active IBD than in the control patients, probably reflecting the suppressive effect of the active inflammation. The onset of the glucocorticoid therapy further suppressed bone turnover. Inflammatory marker hs-CRP decreased readily after the initiation of the steroid, however the decrease did not associate with the clinical response to glucocorticoids. Conclusions: This is the first study to show that adipocyte-derived adiponectin associates with steroid therapy-induced side-effects. Further studies are needed, but it is possible that the adiponectin measurement could aid the recognition of glucocorticoid-sensitive patients in the future. GBA and the other markers reflecting glucocorticoid activity in different tissues changed during the treatment, however their change did not correlate with the therapeutic response to steroids or with the development of glucocorticoid-related side effects and therefore cannot guide the therapy in these patients. Studies such as as the present one that combine clinical data with newly developed biomolecular technology are needed to step-by-step build a general picture of the glucocorticoid actions in different tissues.

Exhaled nitric oxide : Variability and association with bronchial hyperresponsiveness and atopy

Airway inflammation is a key feature of bronchial asthma. In asthma management, according to international guidelines, the gold standard is anti-inflammatory treatment. Currently, only conventional procedures (i.e., symptoms, use of rescue medication, PEF-variability, and lung function tests) were used to both diagnose and evaluate the results of treatment with anti-inflammatory drugs. New methods for evaluation of degree of airway inflammation are required. Nitric oxide (NO) is a gas which is produced in the airways of healthy subjects and especially produced in asthmatic airways. Measurement of NO from the airways is possible, and NO can be measured from exhaled air. Fractional exhaled NO (FENO) is increased in asthma, and the highest concentrations are measured in asthmatic patients not treated with inhaled corticosteroids (ICS). Steroid-treated patients with asthma had levels of FENO similar to those of healthy controls. Atopic asthmatics had higher levels of FENO than did nonatopic asthmatics, indicating that level of atopy affected FENO level. Associations between FENO and bronchial hyperresponsiveness (BHR) occur in asthma. The present study demonstrated that measurement of FENO had good reproducibility, and the FENO variability was reasonable both short- and long-term in both healthy subjects and patients with respiratory symptoms or asthma. We demonstrated the upper normal limit for healthy subjects, which was 12 ppb calculated from two different healthy study populations. We showed that patients with respiratory symptoms who did not fulfil the diagnostic criteria of asthma had FENO values significantly higher than in healthy subjects, but significantly lower than in asthma patients. These findings suggest that BHR to histamine is a sensitive indicator of the effect of ICS and a valuable tool for adjustment of corticosteroid treatment in mild asthma. The findings further suggest that intermittent treatment periods of a few weeks’ duration are insufficient to provide long-term control of BHR in patients with mild persistent asthma. Moreover, during the treatment with ICS changes in BHR and changes in FENO were associated. FENO level was associated with BHR measured by a direct (histamine challenge) or indirect method (exercise challenge) in steroid-naïve symptomatic, non-smoking asthmatics. Although these associations could be found only in atopics, FENO level in nonatopic asthma was also increased. It can thus be concluded that assessment of airway inflammation by measuring FENO can be useful for clinical purposes. The methodology of FENO measurements is now validated. Especially in those patients with respiratory symptoms who did not fulfil the diagnostic criteria of asthma, FENO measurement can aid in treatment decisions. Serial measurement of FENO during treatment with ICS can be a complementary or an alternative method for evaluation in patients with asthma.

Silva Fennica Vol. 21, 1987.

Julkaistu Silva Fennica Vol. 21(4) -numeron liitteenä.

Effects of gender, and SLCO1B1 and CYP7A1 polymorphisms on plasma bile acids in humans and the pharmacokinetics of therapeutic ursodeoxycholic acid

Bile acids are important steroid-derived molecules essential for fat absorption in the small intestine. They are produced in the liver and secreted into the bile. Bile acids are transported by bile flow to the small intestine, where they aid the digestion of lipids. Most bile acids are reabsorbed in the small intestine and return to the liver through the portal vein. The whole recycling process is referred to as the enterohepatic circulation, during which only a small amount of bile acids are removed from the body via faeces. The enterohepatic circulation of bile acids involves the delicate coordination of a number of bile acid transporters expressed in the liver and the small intestine. Organic anion transporting polypeptide 1B1 (OATP1B1), encoded by the solute carrier organic anion transporter family, member 1B1 (SLCO1B1) gene, mediates the sodium independent hepatocellular uptake of bile acids. Two common SNPs in the SLCO1B1 gene are well known to affect the transport activity of OATP1B1. Moreover, bile acid synthesis is an important elimination route for cholesterol. Cholesterol 7α-hydroxylase (CYP7A1) is the rate-limiting enzyme of bile acid production. The aim of this thesis was to investigate the effects of SLCO1B1 polymorphism on the fasting plasma levels of individual endogenous bile acids and a bile acid synthesis marker, and the pharmacokinetics of exogenously administered ursodeoxycholic acid (UDCA). Furthermore, the effects of CYP7A1 genetic polymorphism and gender on the fasting plasma concentrations of individual endogenous bile acids and the bile acid synthesis marker were evaluated. Firstly, a high performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) method for the determination of bile acids was developed (Study I). A retrospective study examined the effects of SLCO1B1 genetic polymorphism on the fasting plasma concentrations of individual bile acids and a bile acid synthesis marker in 65 healthy subjects (Study II). In another retrospective study with 143 healthy individuals, the effects of CYP7A1 genetic polymorphism and gender as well as SLCO1B1 polymorphism on the fasting plasma levels of individual bile acids and the bile acid synthesis marker were investigated (Study III). The effects of SLCO1B1 polymorphism on the pharmacokinetics of exogenously administered UDCA were evaluated in a prospective genotype panel study including 27 healthy volunteers (Study IV). A robust, sensitive and simple HPLC-MS/MS method was developed for the simultaneous determination of 16 individual bile acids in human plasma. The method validation parameters for all the analytes met the requirements of the FDA (Food and Drug Administration) bioanalytical guidelines. This HPLC-MS/MS method was applied in Studies II-IV. In Study II, the fasting plasma concentrations of several bile acids and the bile acid synthesis marker seemed to be affected by SLCO1B1 genetic polymorphism, but these findings were not replicated in Study III with a larger sample size. Moreover, SLCO1B1 polymorphism had no effect on the pharmacokinetic parameters of exogenously administered UDCA. Furthermore, no consistent association was observed between CYP7A1 genetic polymorphism and the fasting plasma concentrations of individual bile acids or the bile acid synthesis marker. In contrast, gender had a major effect on the fasting plasma concentrations of several bile acids and also total bile acids. In conclusion, gender, but not SLCO1B1 or CYP7A1 polymorphisms, has a major effect on the fasting plasma concentrations of individual bile acids. Moreover, the common genetic polymorphism of CYP7A1 is unlikely to influence the activity of CYP7A1 under normal physiological conditions. OATP1B1 does not play an important role in the in vivo disposition of exogenously administered UDCA.