Rational design of serine protease inhibitors

Proteases regulate a spectrum of diverse physiological processes, and dysregulation of proteolytic activity drives a plethora of pathological conditions. Understanding protease function is essential to appreciating many aspects of normal physiology and progression of disease. Consequently, development of potent and specific inhibitors of proteolytic enzymes is vital to provide tools for the dissection of protease function in biological systems and for the treatment of diseases linked to aberrant proteolytic activity. The studies in this thesis describe the rational design of potent inhibitors of three proteases that are implicated in disease development. Additionally, key features of the interaction of proteases and their cognate inhibitors or substrates are analysed and a series of rational inhibitor design principles are expounded and tested. Rational design of protease inhibitors relies on a comprehensive understanding of protease structure and biochemistry. Analysis of known protease cleavage sites in proteins and peptides is a commonly used source of such information. However, model peptide substrate and protein sequences have widely differing levels of backbone constraint and hence can adopt highly divergent structures when binding to a protease’s active site. This may result in identical sequences in peptides and proteins having different conformations and diverse spatial distribution of amino acid functionalities. Regardless of this, protein and peptide cleavage sites are often regarded as being equivalent. One of the key findings in the following studies is a definitive demonstration of the lack of equivalence between these two classes of substrate and invalidation of the common practice of using the sequences of model peptide substrates to predict cleavage of proteins in vivo. Another important feature for protease substrate recognition is subsite cooperativity. This type of cooperativity is commonly referred to as protease or substrate binding subsite cooperativity and is distinct from allosteric cooperativity, where binding of a molecule distant from the protease active site affects the binding affinity of a substrate. Subsite cooperativity may be intramolecular where neighbouring residues in substrates are interacting, affecting the scissile bond’s susceptibility to protease cleavage. Subsite cooperativity can also be intermolecular where a particular residue’s contribution to binding affinity changes depending on the identity of neighbouring amino acids. Although numerous studies have identified subsite cooperativity effects, these findings are frequently ignored in investigations probing subsite selectivity by screening against diverse combinatorial libraries of peptides (positional scanning synthetic combinatorial library; PS-SCL). This strategy for determining cleavage specificity relies on the averaged rates of hydrolysis for an uncharacterised ensemble of peptide sequences, as opposed to the defined rate of hydrolysis of a known specific substrate. Further, since PS-SCL screens probe the preference of the various protease subsites independently, this method is inherently unable to detect subsite cooperativity. However, mean hydrolysis rates from PS-SCL screens are often interpreted as being comparable to those produced by single peptide cleavages. Before this study no large systematic evaluation had been made to determine the level of correlation between protease selectivity as predicted by screening against a library of combinatorial peptides and cleavage of individual peptides. This subject is specifically explored in the studies described here. In order to establish whether PS-SCL screens could accurately determine the substrate preferences of proteases, a systematic comparison of data from PS-SCLs with libraries containing individually synthesised peptides (sparse matrix library; SML) was carried out. These SML libraries were designed to include all possible sequence combinations of the residues that were suggested to be preferred by a protease using the PS-SCL method. SML screening against the three serine proteases kallikrein 4 (KLK4), kallikrein 14 (KLK14) and plasmin revealed highly preferred peptide substrates that could not have been deduced by PS-SCL screening alone. Comparing protease subsite preference profiles from screens of the two types of peptide libraries showed that the most preferred substrates were not detected by PS SCL screening as a consequence of intermolecular cooperativity being negated by the very nature of PS SCL screening. Sequences that are highly favoured as result of intermolecular cooperativity achieve optimal protease subsite occupancy, and thereby interact with very specific determinants of the protease. Identifying these substrate sequences is important since they may be used to produce potent and selective inhibitors of protolytic enzymes. This study found that highly favoured substrate sequences that relied on intermolecular cooperativity allowed for the production of potent inhibitors of KLK4, KLK14 and plasmin. Peptide aldehydes based on preferred plasmin sequences produced high affinity transition state analogue inhibitors for this protease. The most potent of these maintained specificity over plasma kallikrein (known to have a very similar substrate preference to plasmin). Furthermore, the efficiency of this inhibitor in blocking fibrinolysis in vitro was comparable to aprotinin, which previously saw clinical use to reduce perioperative bleeding. One substrate sequence particularly favoured by KLK4 was substituted into the 14 amino acid, circular sunflower trypsin inhibitor (SFTI). This resulted in a highly potent and selective inhibitor (SFTI-FCQR) which attenuated protease activated receptor signalling by KLK4 in vitro. Moreover, SFTI-FCQR and paclitaxel synergistically reduced growth of ovarian cancer cells in vitro, making this inhibitor a lead compound for further therapeutic development. Similar incorporation of a preferred KLK14 amino acid sequence into the SFTI scaffold produced a potent inhibitor for this protease. However, the conformationally constrained SFTI backbone enforced a different intramolecular cooperativity, which masked a KLK14 specific determinant. As a consequence, the level of selectivity achievable was lower than that found for the KLK4 inhibitor. Standard mechanism inhibitors such as SFTI rely on a stable acyl-enzyme intermediate for high affinity binding. This is achieved by a conformationally constrained canonical binding loop that allows for reformation of the scissile peptide bond after cleavage. Amino acid substitutions within the inhibitor to target a particular protease may compromise structural determinants that support the rigidity of the binding loop and thereby prevent the engineered inhibitor reaching its full potential. An in silico analysis was carried out to examine the potential for further improvements to the potency and selectivity of the SFTI-based KLK4 and KLK14 inhibitors. Molecular dynamics simulations suggested that the substitutions within SFTI required to target KLK4 and KLK14 had compromised the intramolecular hydrogen bond network of the inhibitor and caused a concomitant loss of binding loop stability. Furthermore in silico amino acid substitution revealed a consistent correlation between a higher frequency of formation and the number of internal hydrogen bonds of SFTI-variants and lower inhibition constants. These predictions allowed for the production of second generation inhibitors with enhanced binding affinity toward both targets and highlight the importance of considering intramolecular cooperativity effects when engineering proteins or circular peptides to target proteases. The findings from this study show that although PS-SCLs are a useful tool for high throughput screening of approximate protease preference, later refinement by SML screening is needed to reveal optimal subsite occupancy due to cooperativity in substrate recognition. This investigation has also demonstrated the importance of maintaining structural determinants of backbone constraint and conformation when engineering standard mechanism inhibitors for new targets. Combined these results show that backbone conformation and amino acid cooperativity have more prominent roles than previously appreciated in determining substrate/inhibitor specificity and binding affinity. The three key inhibitors designed during this investigation are now being developed as lead compounds for cancer chemotherapy, control of fibrinolysis and cosmeceutical applications. These compounds form the basis of a portfolio of intellectual property which will be further developed in the coming years.

Defining the role of phytoene synthase in carotenoid accumulation of high provitamin A bananas

Vitamin A deficiency (VAD) is a serious problem in developing countries, affecting approximately 127 million children of preschool age and 7.2 million pregnant women each year. However, this deficiency is readily treated and prevented through adequate nutrition. This can potentially be achieved through genetically engineered biofortification of staple food crops to enhance provitamin A (pVA) carotenoid content. Bananas are the fourth most important food crop with an annual production of 100 million tonnes and are widely consumed in areas affected by VAD. However, the fruit pVA content of most widely consumed banana cultivars is low (~ 0.2 to 0.5 ìg/g dry weight). This includes cultivars such as the East African highland banana (EAHB), the staple crop in countries such as Uganda, where annual banana consumption is approximately 250 kg per person. This fact, in addition to the agronomic properties of staple banana cultivars such as vegetative reproduction and continuous cropping, make bananas an ideal target for pVA enhancement through genetic engineering. Interestingly, there are banana varieties known with high fruit pVA content (up to 27.8 ìg/g dry weight), although they are not widely consumed due to factors such as cultural preference and availability. The genes involved in carotenoid accumulation during banana fruit ripening have not been well studied and an understanding of the molecular basis for the differential capacity of bananas to accumulate carotenoids may impact on the effective production of genetically engineered high pVA bananas. The production of phytoene by the enzyme phytoene synthase (PSY) has been shown to be an important rate limiting determinant of pVA accumulation in crop systems such as maize and rice. Manipulation of this gene in rice has been used successfully to produce Golden Rice, which exhibits higher seed endosperm pVA levels than wild type plants. Therefore, it was hypothesised that differences between high and low pVA accumulating bananas could be due either to differences in PSY enzyme activity or factors regulating the expression of the psy gene. Therefore, the aim of this thesis was to investigate the role of PSY in accumulation of pVA in banana fruit of representative high (Asupina) and low (Cavendish) pVA banana cultivars by comparing the nucleic acid and encoded amino acid sequences of the banana psy genes, in vivo enzyme activity of PSY in rice callus and expression of PSY through analysis of promoter activity and mRNA levels. Initially, partial sequences of the psy coding region from five banana cultivars were obtained using reverse transcriptase (RT)-PCR with degenerate primers designed to conserved amino acids in the coding region of available psy sequences from other plants. Based on phylogenetic analysis and comparison to maize psy sequences, it was found that in banana, psy occurs as a gene family of at least three members (psy1, psy2a and psy2b). Subsequent analysis of the complete coding regions of these genes from Asupina and Cavendish suggested that they were all capable of producing functional proteins due to high conservation in the catalytic domain. However, inability to obtain the complete mRNA sequences of Cavendish psy2a, and isolation of two non-functional Cavendish psy2a coding region variants, suggested that psy2a expression may be impaired in Cavendish. Sequence analysis indicated that these Cavendish psy2a coding region variants may have resulted from alternate splicing. Evidence of alternate splicing was also observed in one Asupina psy1 coding region variant, which was predicted to produce a functional PSY1 isoform. The complete mRNA sequence of the psy2b coding regions could not be isolated from either cultivar. Interestingly, psy1 was cloned predominantly from leaf while psy2 was obtained preferentially from fruit, suggesting some level of tissue-specific expression. The Asupina and Cavendish psy1 and psy2a coding regions were subsequently expressed in rice callus and the activity of the enzymes compared in vivo through visual observation and quantitative measurement of carotenoid accumulation. The maize B73 psy1 coding region was included as a positive control. After several weeks on selection, regenerating calli showed a range of colours from white to dark orange representing various levels of carotenoid accumulation. These results confirmed that the banana psy coding regions were all capable of producing functional enzymes. No statistically significant differences in levels of activity were observed between banana PSYs, suggesting that differences in PSY activity were not responsible for differences in the fruit pVA content of Asupina and Cavendish. The psy1 and psy2a promoter sequences were isolated from Asupina and Cavendish gDNA using a PCR-based genome walking strategy. Interestingly, three Cavendish psy2a promoter clones of different sizes, representing possible allelic variants, were identified while only single promoter sequences were obtained for the other Asupina and Cavendish psy genes. Bioinformatic analysis of these sequences identified motifs that were previously characterised in the Arabidopsis psy promoter. Notably, an ATCTA motif associated with basal expression in Arabidopsis was identified in all promoters with the exception of two of the Cavendish psy2a promoter clones (Cpsy2apr2 and Cpsy2apr3). G1 and G2 motifs, linked to light-regulated responses in Arabidopsis, appeared to be differentially distributed between psy1 and psy2a promoters. In the untranscribed regulatory regions, the G1 motifs were found only in psy1 promoters, while the G2 motifs were found only in psy2a. Interestingly, both ATCTA and G2 motifs were identified in the 5’ UTRs of Asupina and Cavendish psy1. Consistent with other monocot promoters, introns were present in the Asupina and Cavendish psy1 5’ UTRs, while none were observed in the psy2a 5’ UTRs. Promoters were cloned into expression constructs, driving the â-glucuronidase (GUS) reporter gene. Transient expression of the Asupina and Cavendish psy1 and psy2a promoters in both Cavendish embryogenic cells and Cavendish fruit demonstrated that all promoters were active, except Cpsy2apr2 and Cpsy2apr3. The functional Cavendish psy2a promoter (Cpsy2apr1) appeared to have activity similar to the Asupina psy2a promoter. The activities of the Asupina and Cavendish psy1 promoters were similar to each other, and comparable to those of the functional psy2a promoters. Semi-quantitative PCR analysis of Asupina and Cavendish psy1 and psy2a transcripts showed that psy2a levels were high in green fruit and decreased during ripening, reinforcing the hypothesis that fruit pVA levels were largely dependent on levels of psy2a expression. Additionally, semi-quantitative PCR using intron-spanning primers indicated that high levels of unprocessed psy2a and psy2b mRNA were present in the ripe fruit of Cavendish but not in Asupina. This raised the possibility that differences in intron processing may influence pVA accumulation in Asupina and Cavendish. In this study the role of PSY in banana pVA accumulation was analysed at a number of different levels. Both mRNA accumulation and promoter activity of psy genes studied were very similar between Asupina and Cavendish. However, in several experiments there was evidence of cryptic or alternate splicing that differed in Cavendish compared to Asupina, although these differences were not conclusively linked to the differences in fruit pVA accumulation between Asupina and Cavendish. Therefore, other carotenoid biosynthetic genes or regulatory mechanisms may be involved in determining pVA levels in these cultivars. This study has contributed to an increased understanding of the role of PSY in the production of pVA carotenoids in banana fruit, corroborating the importance of this enzyme in regulating carotenoid production. Ultimately, this work may serve to inform future research into pVA accumulation in important crop varieties such as the EAHB and the discovery of avenues to improve such crops through genetic modification.

Proteasome inhibitors as anti-cancer agents

The ubiquitin (Ub)-proteasome pathway is the major nonlysosomal pathway of proteolysis in human cells and accounts for the degradation of most short-lived, misfolded or damaged proteins. This pathway is important in the regulation of a number of key biological regulatory mechanisms. Proteins are usually targeted for proteasome-mediated degradation by polyubiquitinylation, the covalent addition of multiple units of the 76 amino acid protein Ub, which are ligated to 1-amino groups of lysine residues in the substrate. Polyubiquitinylated proteins are degraded by the 26S proteasome, a large, ATP-dependent multicatalytic protease complex, which also regenerates monomeric Ub. The targets of this pathway include key regulators of cell proliferation and cell death. An alternative form of the proteasome, termed the immunoproteasome, also has important functions in the generation of peptides for presentation by MHC class I molecules. In recent years there has been a great deal of interest in the possibility that proteasome inhibitors, through elevation of the levels of proteasome targets, might prove useful as a novel class of anti-cancer drugs. Here we review the progress made to date in this area and highlight the potential advantages and weaknesses of this approach.

Viruses of banana in East Africa

Bananas are one of the world's most important food crops, providing sustenance and income for millions of people in developing countries and supporting large export industries. Viruses are considered major constraints to banana production, germplasm multiplication and exchange, and to genetic improvement of banana through traditional breeding. In Africa, the two most important virus diseases are bunchy top, caused by Banana bunchy top virus (BBTV), and banana streak disease, caused by Banana streak virus (BSV). BBTV is a serious production constraint in a number of countries within/bordering East Africa, such as Burundi, Democratic Republic of Congo, Malawi, Mozambique, Rwanda and Zambia, but is not present in Kenya, Tanzania and Uganda. Additionally, epidemics of banana streak disease are occurring in Kenya and Uganda. The rapidly growing tissue culture (TC) industry within East Africa, aiming to provide planting material to banana farmers, has stimulated discussion about the need for virus indexing to certify planting material as virus-free. Diagnostic methods for BBTV and BSV have been reported and, for BBTV, PCR-based assays are reliable and relatively straightforward. However for BSV, high levels of serological and genetic variability and the presence of endogenous virus sequences within the banana genome complicate diagnosis. Uganda has been shown to contain the greatest diversity in BSV isolates found anywhere in the world. A broad-spectrum diagnostic test for BSV detection, which can discriminate between endogenous and episomal BSV sequences, is a priority. This PhD project aimed to establish diagnostic methods for banana viruses, with a particular focus on the development of novel methods for BSV detection, and to use these diagnostic methods for the detection and characterisation of banana viruses in East Africa. A novel rolling-circle amplification (RCA) method was developed for the detection of BSV. Using samples of Banana streak MY virus (BSMYV) and Banana streak OL virus (BSOLV) from Australia, this method was shown to distinguish between endogenous and episomal BSV sequences in banana plants. The RCA assay was used to screen a collection of 56 banana samples from south-west Uganda for BSV. RCA detected at least five distinct BSV isolates in these samples, including BSOLV and Banana streak GF virus (BSGFV) as well as three BSV isolates (Banana streak Uganda-I, -L and -M virus) for which only partial sequences had been previously reported. These latter three BSV had only been detected using immuno-capture (IC)-PCR and thus were possible endogenous sequences. In addition to its ability to detect BSV, the RCA protocol was also demonstrated to detect other viruses within the family Caulimoviridae, including Sugar cane bacilliform virus, and Cauliflower mosaic virus. Using the novel RCA method, three distinct BSV isolates from both Kenya and Uganda were identified and characterised. The complete genome of these isolates was sequenced and annotated. All six isolates were shown to have a characteristic badnavirus genome organisation with three open reading frames (ORFs) and the large polyprotein encoded by ORF 3 was shown to contain conserved amino acid motifs for movement, aspartic protease, reverse transcriptase and ribonuclease H activities. As well, several sequences important for expression and replication of the virus genome were identified including the conserved tRNAmet primer binding site present in the intergenic region of all badnaviruses. Based on the International Committee on Taxonomy of Viruses (ICTV) guidelines for species demarcation in the genus Badnavirus, these six isolates were proposed as distinct species, and named Banana streak UA virus (BSUAV), Banana streak UI virus (BSUIV), Banana streak UL virus (BSULV), Banana streak UM virus (BSUMV), Banana streak CA virus (BSCAV) and Banana streak IM virus (BSIMV). Using PCR with species-specific primers designed to each isolate, a genotypically diverse collection of 12 virus-free banana cultivars were tested for the presence of endogenous sequences. For five of the BSV no amplification was observed in any cultivar tested, while for BSIMV, four positive samples were identified in cultivars with a B-genome component. During field visits to Kenya, Tanzania and Uganda, 143 samples were collected and assayed for BSV. PCR using nine sets of species-specific primers, and RCA, were compared for BSV detection. For five BSV species with no known endogenous counterpart (namely BSCAV, BSUAV, BSUIV, BSULV and BSUMV), PCR was used to detect 30 infections from the 143 samples. Using RCA, 96.4% of these samples were considered positive, with one additional sample detected using RCA which was not positive using PCR. For these five BSV, PCR and RCA were both useful for identifying infected samples, irrespective of the host cultivar genotype (Musa A- or B-genome components). For four additional BSV with known endogenous counterparts in the M. balbisiana genome (BSOLV, BSGFV, BSMYV and BSIMV), PCR was shown to detect 75 infections from the 143 samples. In 30 samples from cultivars with an A-only genome component there was 96.3% agreement between PCR positive samples and detection using RCA, again demonstrating either PCR or RCA are suitable methods for detection. However, in 45 samples from cultivars with some B-genome component, the level of agreement between PCR positive samples and RCA positive samples was 70.5%. This suggests that, in cultivars with some B-genome component, many infections were detected using PCR which were the result of amplification of endogenous sequences. In these latter cases, RCA or another method which discriminates between endogenous and episomal sequences, such as immuno-capture PCR, is needed to diagnose episomal BSV infection. Field visits were made to Malawi and Rwanda to collect local isolates of BBTV for validation of a PCR-based diagnostic assay. The presence of BBTV in samples of bananas with bunchy top disease was confirmed in 28 out of 39 samples from Malawi and all nine samples collected in Rwanda, using PCR and RCA. For three isolates, one from Malawi and two from Rwanda, the complete nucleotide sequences were determined and shown to have a similar genome organisation to previously published BBTV isolates. The two isolates from Rwanda had at least 98.1% nucleotide sequence identity between each of the six DNA components, while the similarity between isolates from Rwanda and Malawi was between 96.2% and 99.4% depending on the DNA component. At the amino acid level, similarities in the putative proteins encoded by DNA-R, -S, -M, - C and -N were found to range between 98.8% to 100%. In a phylogenetic analysis, the three East African isolates clustered together within the South Pacific subgroup of BBTV isolates. Nucleotide sequence comparison to isolates of BBTV from outside Africa identified India as the possible origin of East African isolates of BBTV.

Mastering the canonical loop of serine protease inhibitors : enhancing potency by optimising the internal hydrogen bond network

Background Canonical serine protease inhibitors commonly bind to their targets through a rigid loop stabilised by an internal hydrogen bond network and disulfide bond(s). The smallest of these is sunflower trypsin inhibitor (SFTI-1), a potent and broad-range protease inhibitor. Recently, we re-engineered the contact β-sheet of SFTI-1 to produce a selective inhibitor of kallikrein-related peptidase 4 (KLK4), a protease associated with prostate cancer progression. However, modifications in the binding loop to achieve specificity may compromise structural rigidity and prevent re-engineered inhibitors from reaching optimal binding affinity. Methodology/Principal Findings In this study, the effect of amino acid substitutions on the internal hydrogen bonding network of SFTI were investigated using an in silico screen of inhibitor variants in complex with KLK4 or trypsin. Substitutions favouring internal hydrogen bond formation directly correlated with increased potency of inhibition in vitro. This produced a second generation inhibitor (SFTI-FCQR Asn14) which displayed both a 125-fold increased capacity to inhibit KLK4 (Ki = 0.0386±0.0060 nM) and enhanced selectivity over off-target serine proteases. Further, SFTI-FCQR Asn14 was stable in cell culture and bioavailable in mice when administered by intraperitoneal perfusion. Conclusion/Significance These findings highlight the importance of conserving structural rigidity of the binding loop in addition to optimising protease/inhibitor contacts when re-engineering canonical serine protease inhibitors.

Identification of vitronectin as a novel insulin-like growth factor-II binding protein

We have previously reported the presence of a 70 kDa insulin-like growth factor (IGF)-II-specific binding protein in chicken serum using Western ligand blotting approaches. In order to ascertain the identity of this 70 kDa IGF-II binding species, the protein has been purified from chicken serum using a combination of ion-exchange and gel-permeation chromatography. Interestingly, amino acid sequencing of the purified protein revealed that it has the same N-terminal sequence as chicken vitronectin (VN). The protein has the ability to specifically bind IGF-II and not IGF-I as determined by ligand blotting, cross-linking and competitive binding assay approaches. In addition, the protein binds 125I-des(l-6)-IGF-II, suggesting that the interaction with IGF-II is different to those with other characterized IGF-binding proteins. Importantly, we have ascertained that both human and bovine VN also specifically bind IGF-II. These results are particularly relevant in the light of the recent report that the urokinase-type plasminogen activator receptor, a protein that also binds VN, has been shown to associate with the cation-independent mannose-6-phosphate/GF-II receptor and suggest a possible role for IGF-II in cell adhesion and invasion.

The roles of the preproghrelin-derived peptides - ghrelin, desacyl ghrelin and obestatin - in prostate cancer

Prostate cancer is the second most common cause of cancer related deaths in Western men. Despite the significant improvements in current treatment techniques, there is no cure for advanced metastatic, castrate-resistant disease. Early detection and prevention of progression to a castrate-resistant state may provide new strategies to improve survival. A number of growth factors have been shown to act in an autocrine/paracrine manner to modulate prostate cancer tumour growth. Our laboratory has previously shown that ghrelin and its receptors (the functional GHS-R1a and the non-functional GHS-R1b) are expressed in prostate cancer specimens and cell lines. We have shown that ghrelin increases cell proliferation in the PC3 and LNCaP prostate cancer cell lines through activation of ERK1/2, suggesting that ghrelin could regulate prostate cancer cell growth and play a role in the progression of the disease. Ghrelin is a 28 amino-acid peptide hormone, identified to be the natural ligand of the growth hormone secretagogue receptor (GHS-R1a). It is well characterised as a growth hormone releasing and as an orexigenic peptide that stimulates appetite and feeding and regulates energy expenditure and bodyweight. In addition to its orexigenic properties, ghrelin has been shown to play a regulatory role in a number of systems, including the reproductive, immune and cardiovascular systems and may play a role in a number of pathological conditions such as chronic heart failure, anorexia, cachexia, obesity, diabetes and cancer. In cancer, ghrelin and its receptor are expressed in a range of tumours and cancer cell lines and ghrelin has been demonstrated to modulate cell proliferation, apoptosis, migration and invasion in some cell types. The ghrelin gene (GHRL) encodes preproghrelin peptide, which is processed to produce three currently known functional peptides - ghrelin, desacyl ghrelin and obestatin. Prohormone convertases (PCs) have been shown to cleave the preproghrelin peptide into two primary products - the 28 amino acid peptide, ghrelin, and the remaining 117 amino acid C-terminal peptide, C-ghrelin. C-ghrelin can then be further processed to produce the 23 amino acid peptide, obestatin. Ghrelin circulates in two different forms - an octanoylated form (known as ghrelin) and a non-octanoylated form, desacyl ghrelin. The unique post-translational addition of octanoic acid to the serine 3 residue of the propeptide chain to form acylated ghrelin is catalysed by ghrelin O-acyltransferase (GOAT). This modification is necessary for binding of ghrelin to its only known functional receptor, the GHS-R1a. As desacyl ghrelin cannot bind and activate the GHS-R1a, it was initially thought to be an inactive peptide, despite the fact that it circulates at much higher levels than ghrelin. Further research has demonstrated that desacyl ghrelin is biologically active and shares some of the actions of ghrelin, as well as having some opposing and distinct roles. Interestingly, both ghrelin and desacyl ghrelin have been shown to modulate apoptosis, cell differentiation and proliferation in some cell types, and to stimulate cell proliferation through activation of ERK1/2 and PI3K/Akt pathways. The third known peptide product of the ghrelin preprohormone, obestatin, was initially thought to oppose the actions of ghrelin in appetite regulation and food intake and to mediate its effects through the G protein-coupled receptor 39 (GPR39). Subsequent research failed to reproduce the initial findings, however, and the possible anorexigenic effects of obestatin, as well as the identity of its receptor, remain unclear. Obestatin plays some important physiological roles, including roles in improving memory, the inhibition of thirst and anxiety, increased secretion of pancreatic juice, and regulation of cell proliferation, survival, apoptosis and differentiation. Preliminary studies have also shown that obestatin stimulates cell proliferation in some cell types through activation of ERK1/2, Akt and PKC pathways. Overall, however, at the commencement of this PhD project, relatively little was known regarding the functions and mechanisms of action of the preproghrelin-derived functional peptides in modulating prostate cancer cell proliferation. The roles of obestatin, and desacyl ghrelin as potential growth factors had not previously been investigated, and the potential expression and regulation of the preproghrelin processing enzymes, GOAT and prohormone convertases was unknown in prostate cancer cell lines. Therefore, the overall objectives of this study were to: 1. investigate the effects of obestatin on cell proliferation and signaling in prostate cancer cell lines 2. compare the effects of desacyl ghrelin and ghrelin on cell proliferation and signaling in prostate cancer cell lines 3. investigate whether prostate cancer cell lines possess the necessary enzymatic machinery to produce ghrelin and desacyl ghrelin and if these peptides can regulate GOAT expression Our laboratory has previously shown that ghrelin stimulates cell proliferation in the PC3 and LNCaP prostate cancer cell line through activation of the ERK1/2 pathway. In this study it has been demonstrated that treatments with either ghrelin, desacyl ghrelin or obestatin over 72 hours significantly increased cell proliferation in the PC3 prostate cancer cell line but had no significant effect in the RWPE-1 transformed normal prostate cell line. Ghrelin (1000nM) stimulated cell proliferation in the PC3 prostate cancer cell line by 31.66 6.68% (p<0.01) with the WST-1 method, and 13.55 5.68% (p<0.05) with the CyQUANT assay. Desacyl ghrelin (1000nM) increased cell proliferation in PC3 cells by 21.73 2.62% (p<0.01) (WST-1), and 15.46 7.05% (p<0.05) (CyQUANT) above untreated control. Obestatin (1000nM) induced a 28.37 7.47% (p<0.01) (WST-1) and 12.14 7.47% (p<0.05) (CyQUANT) significant increase in cell proliferation in the PC3 prostate cancer cell line. Ghrelin and desacyl ghrelin treatments stimulated Akt and ERK phosphorylation across a range of concentrations (p<0.01). Obestatin treatment significantly stimulated Akt, ERK and PKC phosphorylation (p<0.05). Through the use of specific inhibitors, the MAPK inhibitor U0126 and the Akt1/2 kinase inhibitor, it was demonstrated that ghrelin- and obestatin-induced cell proliferation in the PC3 prostate cancer cell line is mediated through activation of ERK1/2 and Akt pathways. Although desacyl ghrelin significantly stimulated Akt and ERK phosphorylation, U0126 failed to prevent desacyl ghrelin-induced cell proliferation suggesting ghrelin and desacyl ghrelin might act through different mechanisms to increase cell proliferation. Ghrelin and desacyl ghrelin have shown a proliferative effect in osteoblasts, pancreatic -cells and cardiomyocytes through activation of ERK1/2 and PI3K/Akt pathways. Here it has been shown that ghrelin and its non-acylated form exert the same function and stimulate cell proliferation in the PC3 prostate cancer cell line through activation of the Akt pathway. Ghrelin-induced proliferation was also mediated through activation of the ERK1/2 pathway, however, desacyl ghrelin seems to stimulate cell proliferation in an ERK1/2-independent manner. As desacyl ghrelin does not bind and activate GHSR1a, the only known functional ghrelin receptor, the finding that both ghrelin and desacyl ghrelin stimulate cell proliferation in the PC3 cell line suggests that these peptides could be acting through the yet unidentified alternative ghrelin receptor in this cell type. Obestatin treatment also stimulated PKC phosphorylation, however, a direct role for this pathway in stimulating cell proliferation could not be proven using available PKC pathway inhibitors, as they caused significant cell death over the extended timeframe of the cell proliferation assays. Obestatin has been shown to stimulate cell proliferation through activation of PKC isoforms in human retinal epithelial cells and in the human gastric cancer cell line KATO-III. We have demonstrated that all of the prostate-derived cell lines examined (PC3, LNCaP, DU145, 22Rv1, RWPE-1 and RWPE-2) expressed GOAT and at least one of the prohormone convertases, which are known to cleave the proghrelin peptide, PC1/3, PC2 and furin, at the mRNA level. These cells, therefore, are likely to possess the necessary machinery to cleave the preproghrelin protein and to produce the mature ghrelin and desacyl ghrelin peptides. In addition to prohormone convertases, the presence of octanoic acid is essential for acylated ghrelin production. In this study octanoic acid supplementation significantly increased cell proliferation in the PC3 prostate cancer cell line by over 20% compared to untreated controls (p<0.01), but surprisingly, not in the DU145, LNCaP or 22Rv1 prostate cancer cell lines or in the RWPE-1 and RWPE-2 prostate-derived cell lines. In addition, we demonstrated that exogenous ghrelin induced a statistically significant two-fold decrease in GOAT mRNA expression in the PC3 cell line (p<0.05), suggesting that ghrelin could pontentially downregulate its own acylation and, therefore, regulate the balance between ghrelin and desacyl ghrelin. This was not observed, however, in the DU145 and LNCaP prostate cancer cell lines. The GOAT-ghrelin system represents a direct link between ingested nutrients and regulation of ghrelin production and the ghrelin/desacyl ghrelin ratio. Regulation of ghrelin acylation is a potentially attractive and desirable tool for the development of better therapies for a number of pathological conditions where ghrelin has been shown to play a key role. The finding that desacyl ghrelin stimulates cell proliferation in the PC3 prostate cancer cell line, and responds to ghrelin in the same way, suggests that this cell line expresses an alternative ghrelin receptor. Although all the cell lines examined expressed both GHS-R1a and GHS-R1b mRNA, it remains uncertain whether these cell lines express the unidentified alternative ghrelin receptor. It is possible that the varied responses seen could be due to the expression of different ghrelin receptors in different cell lines. In addition to GOAT, prohormone convertases and octanoic acid availability may regulate the production of different peptides from the ghrelin preprohormone. The studies presented in this thesis provide significant new information regarding the roles and mechanisms of action of the preproghrelin-derived peptides, ghrelin, desacyl ghrelin and obestatin, in modulating prostate cancer cell line proliferation. A number of key questions remain to be resolved, however, including the identification of the alternative ghrelin/desacyl ghrelin receptor, the identification of the obestatin receptor, a clarification of the signaling mechanisms which mediate cell proliferation in response to obestatin treatment and a better understanding of the regulation at both the gene and post-translational levels of functional peptide generation. Further studies investigating the role of the ghrelin axis using in vivo prostate cancer models may be warranted. Until these issues are determined, the potential for the ghrelin axis, to be recognised as a novel useful target for therapy for cancer or other pathologies will be uncertain.

A butterfly eye's view of birds

The striking color patterns of butterflies and birds have long interested biologists. But how these animals see color is less well understood. Opsins are the protein components of the visual pigments of the eye. Color vision has evolved in butterflies through opsin gene duplications, through positive selection at individual opsin loci, and by the use of filtering pigments. By contrast, birds have retained the same opsin complement present in early-jawed vertebrates, and their visual system has diversified primarily through tuning of the short-wavelength-sensitive photoreceptors, rather than by opsin duplication or the use of filtering elements. Butterflies and birds have evolved photoreceptors that might use some of the same amino acid sites for generating similar spectral phenotypes across approximately 540 million years of evolution, when rhabdomeric and ciliary-type opsins radiated during the early Cambrian period. Considering the similarities between the two taxa, it is surprising that the eyes of birds are not more diverse. Additional taxonomic sampling of birds may help clarify this mystery.

Gene duplication is an evolutionary mechanism for expanding spectral diversity in the long-wavelength photopigments of butterflies

Butterfly long-wavelength (L) photopigments are interesting for comparative studies of adaptive evolution because of the tremendous phenotypic variation that exists in their wavelength of peak absorbance (lambda(max) value). Here we present a comprehensive survey of L photopigment variation by measuring lambda(max) in 12 nymphalid and 1 riodinid species using epi-microspectrophotometry. Together with previous data, we find that L photopigment lambda(max) varies from 510-565 nm in 22 nymphalids, with an even broader 505- to 600-nm range in riodinids. We then surveyed the L opsin genes for which lambda(max) values are available as well as from related taxa and found 2 instances of L opsin gene duplication within nymphalids, in Hermeuptychia hermes and Amathusia phidippus, and 1 instance within riodinids, in the metalmark butterfly Apodemia mormo. Using maximum parsimony and maximum likelihood ancestral state reconstructions to map the evolution of spectral shifts within the L photopigments of nymphalids, we estimate the ancestral pigment had a lambda(max) = 540 nm +/- 10 nm standard error and that blueshifts in wavelength have occurred at least 4 times within the family. We used ancestral state reconstructions to investigate the importance of several amino acid substitutions (Ile17Met, Ala64Ser, Asn70Ser, and Ser137Ala) previously shown to have evolved under positive selection that are correlated with blue spectral shifts. These reconstructions suggest that the Ala64Ser substitution has indeed occurred along the newly identified blueshifted L photopigment lineages. Substitutions at the other 3 sites may also be involved in the functional diversification of L photopigments. Our data strongly suggest that there are limits to the evolution of L photopigment spectral shifts among species with only one L opsin gene and that opsin gene duplication broadens the potential range of lambda(max) values.

Adaptive evolution of color vision as seen through the eyes of butterflies

Butterflies and primates are interesting for comparative color vision studies, because both have evolved middle- (M) and long-wavelength- (L) sensitive photopigments with overlapping absorbance spectrum maxima (lambda(max) values). Although positive selection is important for the maintenance of spectral variation within the primate pigments, it remains an open question whether it contributes similarly to the diversification of butterfly pigments. To examine this issue, we performed epimicrospectrophotometry on the eyes of five Limenitis butterfly species and found a 31-nm range of variation in the lambda(max) values of the L-sensitive photopigments (514-545 nm). We cloned partial Limenitis L opsin gene sequences and found a significant excess of replacement substitutions relative to polymorphisms among species. Mapping of these L photopigment lambda(max) values onto a phylogeny revealed two instances within Lepidoptera of convergently evolved L photopigment lineages whose lambda(max) values were blue-shifted. A codon-based maximum-likelihood analysis indicated that, associated with the two blue spectral shifts, four amino acid sites (Ile17Met, Ala64Ser, Asn70Ser, and Ser137Ala) have evolved substitutions in parallel and exhibit significant d(N)/d(S) >1. Homology modeling of the full-length Limenitis arthemis astyanax L opsin placed all four substitutions within the chromophore-binding pocket. Strikingly, the Ser137Ala substitution is in the same position as a site that in primates is responsible for a 5- to 7-nm blue spectral shift. Our data show that some of the same amino acid sites are under positive selection in the photopigments of both butterflies and primates, spanning an evolutionary distance >500 million years.

Cloning of a novel insulin-regulated ghrelin transcript in prostate cancer

Ghrelin is a multifunctional hormone, with roles in stimulating appetite and regulating energy balance, insulin secretion and glucose homeostasis. The ghrelin gene locus (GHRL) is highly complex and gives rise to a range of novel transcripts derived from alternative first exons and internally spliced exons. The wild-type transcript encodes a 117 amino acid preprohormone that is processed to yield the 28 amino acid peptide ghrelin. Here, we identified insulin-responsive transcription corresponding to cryptic exons in intron 2 of the human ghrelin gene. A transcript, termed in2c-ghrelin (intron 2-cryptic), was cloned from the testis and the LNCaP prostate cancer cell line. This transcript may encode an 83 AA preproghrelin isoform that codes for the ghrelin, but not obestatin. It is expressed in a limited number of normal tissues and in tumours of the prostate, testis, breast and ovary. Finally, we confirmed that in2c-ghrelin transcript expression, as well as the recently described in1-ghrelin transcript, is significantly upregulated by insulin in cultured prostate cancer cells. Metabolic syndrome and hyperinsulinaemia has been associated with prostate cancer risk and progression. This may be particularly significant after androgen deprivation therapy for prostate cancer, which induces hyperinsulinaemia, and this could contribute to castrate resistant prostate cancer growth. We have previously demonstrated that ghrelin stimulates prostate cancer cell line proliferation in vitro. This study is the first description of insulin regulation of a ghrelin transcript in cancer, and should provide further impetus for studies into the expression, regulation and function of ghrelin gene products.

An ensemble method for predicting subnuclear localizations from primary protein structures

Background Predicting protein subnuclear localization is a challenging problem. Some previous works based on non-sequence information including Gene Ontology annotations and kernel fusion have respective limitations. The aim of this work is twofold: one is to propose a novel individual feature extraction method; another is to develop an ensemble method to improve prediction performance using comprehensive information represented in the form of high dimensional feature vector obtained by 11 feature extraction methods. Methodology/Principal Findings A novel two-stage multiclass support vector machine is proposed to predict protein subnuclear localizations. It only considers those feature extraction methods based on amino acid classifications and physicochemical properties. In order to speed up our system, an automatic search method for the kernel parameter is used. The prediction performance of our method is evaluated on four datasets: Lei dataset, multi-localization dataset, SNL9 dataset and a new independent dataset. The overall accuracy of prediction for 6 localizations on Lei dataset is 75.2% and that for 9 localizations on SNL9 dataset is 72.1% in the leave-one-out cross validation, 71.7% for the multi-localization dataset and 69.8% for the new independent dataset, respectively. Comparisons with those existing methods show that our method performs better for both single-localization and multi-localization proteins and achieves more balanced sensitivities and specificities on large-size and small-size subcellular localizations. The overall accuracy improvements are 4.0% and 4.7% for single-localization proteins and 6.5% for multi-localization proteins. The reliability and stability of our classification model are further confirmed by permutation analysis. Conclusions It can be concluded that our method is effective and valuable for predicting protein subnuclear localizations. A web server has been designed to implement the proposed method. It is freely available at http://bioinformatics.awowshop.com/snlpr​ed_page.php.

A description of the Mei2-like protein family; structure, phylogenetic distribution and biological context

The Schizosaccharomyces pombe Mei2 gene encodes an RNA recognition motif (RRM) protein that stimulates meiosis upon binding a specific non-coding RNA and subsequent accumulation in a “mei2-dot” in the nucleus. We present here the first systematic characterization of the family of proteins with characteristic Mei2-like amino acid sequences. Mei2-like proteins are an ancient eukaryotic protein family with three identifiable RRMs. The C-terminal RRM (RRM3) is unique to Mei2-like proteins and is the most highly conserved of the three RRMs. RRM3 also contains conserved sequence elements at its C-terminus not found in other RRM domains. Single copy Mei2-like genes are present in some fungi, in alveolates such as Paramecium and in the early branching eukaryote Entamoeba histolytica, while plants contain small families of Mei2-like genes. While the C-terminal RRM is highly conserved between plants and fungi, indicating conservation of molecular mechanisms, plant Mei2-like genes have changed biological context to regulate various aspects of developmental pattern formation.

Metabolic and ecological study of environmental pentose utilizing bacteria (E-PUB)

Lignocellulosic materials, such as sugar cane bagasse, a waste product of the sugarcane processing industry, agricultural residues and herbaceous crops, may serve as an abundant and comparatively cheap feedstock for largescale industrial fermentation, resulting in the production of marketable end-products. However, the complex structure of lignocellulosic materials, the presence of various hexose and pentose sugars in the hemicellulose component, and the presence of various compounds that inhibit the organisms selected for the fermentation process, all constitute barriers that add to the production costs and make full scale industrial production economically less feasible. The work presented in this thesis was conducted in order to screen microorganisms for ability to utilize pentose sugars derived from the sugar mill industrial waste. A large number of individual bacterial strains were investigated from hemi-cellulose rich material collected at the Proserpine and Maryborough sugar mills, notably soil samples from the mill sites. The research conducted to isolation of six pentose-capable Gram-positive organisms from the actinomycetes group by using pentose as a sole carbon source in the cultivation process. The isolates were identified as Corynebacterium glutamicum, Actinomyces odontolyticus, Nocardia elegans, and Propionibacterium freudenreichii all of which were isolated from the hemicellulose-enriched soil. Pentose degrading microbes are very rare in the environment, so this was a significant discovery. Previous research indicated that microbes could degrade pentose after genetic modification but the microbes discovered in this research were able to naturally utilize pentose. Six isolates, identified as four different genera, were investigated for their ability to utilize single sugars as substrates (glucose, xylose, arabinose or ribose), and also dual sugars as substrates (a hexose plus a pentose). The results demonstrated that C. glutamicum, A. odontolyticus, N. elegans, and P. freudenreichii were pentose-capable (able to grow using xylose or other pentose sugar), and also showed diauxie growth characteristics during the dual-sugar (glucose, in combination with xylose, arabinose or ribose) carbon source tests. In addition, it was shown that the isolates displayed very small differences in growth rates when grown on dual sugars as compared to single sugars, whether pentose or hexose in nature. The anabolic characteristics of C. glutamicum, A. odontolyticus, N. elegans and P. freudenreichii were subsequently investigated by qualitative analysis of their end-products, using high performance liquid chromatography (HPLC). All of the organisms produced arginine and cysteine after utilization of the pentose substrates alone. In addition, P. freudenreichii produced alanine and glycine. The end-product profile arising from culture with dual carbon sources was also tested. Interestingly, this time the product was different. All of them produced the amino acid glycine, when grown on a combination substrate-mix of glucose with xylose, and also glucose with arabinose. Only N. elegans was able to break down ribose, either singly or in combination with glucose, and the end-product of metabolism of the glucose plus ribose substrate combination was glutamic acid. The ecological analysis of microbial abundance in sugar mill waste was performed using denaturing gradient gel electrophoresis (DGGE) and also the metagenomic microarray PhyloChip method. Eleven solid samples and seven liquid samples were investigated. A very complex bacterial ecosystem was demonstrated in the seven liquid samples after testing with the PhyloChip method. It was also shown that bagasse leachate was the most different, compared to all of the other samples, by virtue of its richness in variety of taxa and the complexity of its bacterial community. The bacterial community in solid samples from Proserpine, Mackay and Maryborough sugar mills showed huge diversity. The information found from 16S rDNA sequencing results was that the bacterial genera Brevibacillus, Rhodospirillaceae, Bacillus, Vibrio and Pseudomonas were present in greatest abundance. In addition, Corynebacterium was also found in the soil samples. The metagenomic studies of the sugar mill samples demonstrate two important outcomes: firstly that the bagasse leachate, as potentially the most pentose-rich sample tested, had the most complex and diverse bacterial community; and secondly that the pentose-capable isolates that were initially discovered at the beginning of this study, were not amongst the most abundant taxonomic groups discovered in the sugar mill samples, and in fact were, as suspected, very rare. As a bioprospecting exercise, therefore, the study has discovered organisms that are naturally present, but in very small numbers, in the appropriate natural environment. This has implications for the industrial application of E-PUB, in that a seeding process using a starter culture will be necessary for industrial purposes, rather than simply assuming that natural fermentation might occur.