Multiple analysis of three common genetic alterations associated with thrombophilia

We have previously reported the use of a novel mini-sequencing protocol for detection of the factor V Leiden variant, the first nucleotide change (FNC) technology. This technology is based on a single nucleotide extension of a primer, which is hybridized immediately adjacent to the site of mutation. The extended nucleotide that carries a reporter molecule (fluorescein) has the power to discriminate the genotype at the site of mutation. More recently, the prothrombin 20210 and thermolabile methylene tetrahydrofolate reductase (MTHFR) 677 variants have been identified as possible risk factors associated with thrombophilia. This study describes the use of the FNC technology in a combined assay to detect factor V, prothrombin and MTHFR variants in a population of Australian blood donors, and describes the objective numerical methodology used to determine genotype cut-off values for each genetic variation. Using FNC to test 500 normal blood donors, the incidence of Factor V Leiden was 3.6% (all heterozygous), that of prothrombin 20210 was 2.8% (all heterozygous) and that of MTHFR was 10% (homozygous). The combined FNC technology offers a simple, rapid, automatable DNA-based test for the detection of these three important mutations that are associated with familial thrombophilia. (C) 2000 Lippincott Williams and Wilkins.

Exponentiated gradient algorithms for conditional random fields and max-margin Markov Networks

Log-linear and maximum-margin models are two commonly-used methods in supervised machine learning, and are frequently used in structured prediction problems. Efficient learning of parameters in these models is therefore an important problem, and becomes a key factor when learning from very large data sets. This paper describes exponentiated gradient (EG) algorithms for training such models, where EG updates are applied to the convex dual of either the log-linear or max-margin objective function; the dual in both the log-linear and max-margin cases corresponds to minimizing a convex function with simplex constraints. We study both batch and online variants of the algorithm, and provide rates of convergence for both cases. In the max-margin case, O(1/ε) EG updates are required to reach a given accuracy ε in the dual; in contrast, for log-linear models only O(log(1/ε)) updates are required. For both the max-margin and log-linear cases, our bounds suggest that the online EG algorithm requires a factor of n less computation to reach a desired accuracy than the batch EG algorithm, where n is the number of training examples. Our experiments confirm that the online algorithms are much faster than the batch algorithms in practice. We describe how the EG updates factor in a convenient way for structured prediction problems, allowing the algorithms to be efficiently applied to problems such as sequence learning or natural language parsing. We perform extensive evaluation of the algorithms, comparing them to L-BFGS and stochastic gradient descent for log-linear models, and to SVM-Struct for max-margin models. The algorithms are applied to a multi-class problem as well as to a more complex large-scale parsing task. In all these settings, the EG algorithms presented here outperform the other methods.

Linkage to chromosome 2q32.2-q33.3 in familial serrated neoplasia (Jass syndrome)

Abstract Causative genetic variants have to date been identified for only a small proportion of familial colorectal cancer (CRC). While conditions such as Familial Adenomatous Polyposis and Lynch syndrome have well defined genetic causes, the search for variants underlying the remainder of familial CRC is plagued by genetic heterogeneity. The recent identification of families with a heritable predisposition to malignancies arising through the serrated pathway (familial serrated neoplasia or Jass syndrome) provides an opportunity to study a subset of familial CRC in which heterogeneity may be greatly reduced. A genome-wide linkage screen was performed on a large family displaying a dominantly-inherited predisposition to serrated neoplasia genotyped using the Affymetrix GeneChip Human Mapping 10 K SNP Array. Parametric and nonparametric analyses were performed and resulting regions of interest, as well as previously reported CRC susceptibility loci at 3q22, 7q31 and 9q22, were followed up by finemapping in 10 serrated neoplasia families. Genome-wide linkage analysis revealed regions of interest at 2p25.2-p25.1, 2q24.3-q37.1 and 8p21.2-q12.1. Finemapping linkage and haplotype analyses identified 2q32.2-q33.3 as the region most likely to harbour linkage, with heterogeneity logarithm of the odds (HLOD) 2.09 and nonparametric linkage (NPL) score 2.36 (P = 0.004). Five primary candidate genes (CFLAR, CASP10, CASP8, FZD7 and BMPR2) were sequenced and no segregating variants identified. There was no evidence of linkage to previously reported loci on chromosomes 3, 7 and 9.

Mutation analysis of the CDKN2A promoter in Australian melanoma families

Approximately 50% of all melanoma families worldwide show linkage to 9p21-22, but only about half of these have been shown to contain germ line CDKN2A mutations. It has been hypothesized that a proportion of these families carry mutations in the noncoding regions of CDKN2A. Several Canadian families have been reported to carry a mutation in the 5' UTR, at position -34 relative to the start site, which gives rise to a novel AUG translation initiation codon that markedly decreases translation from the wild-type AUG (Liu et al., 1999). Haplotype sharing in these Canadian families suggested that this mutation is of British origin. We sequenced 1,327 base pairs (bp) of CDKN2A, making up 1,116 bp of the 5' UTR and promoter, all of exon 1, and 61 bp of intron 1, in at least one melanoma case from 110 Australian families with three or more affected members known not to carry mutations within the p16 coding region. In addition, 431 bp upstream of the start codon was sequenced in an additional 253 affected probands from two-case melanoma families for which the CDKN2A mutation status was unknown. Several known polymorphisms at positions -33, -191, -493, and -735 were detected, in addition to four novel variants at positions 120, -252, -347, and -981 relative to the start codon. One of the probands from a two-case family was found to have the previously reported Q50R mutation. No family member was found to carry the mutation at position -34 or any other disease-associated mutation. For further investigation of noncoding CDKN2A mutations that may affect transcription, allele-specific expression analysis was carried out in 31 of the families with at least three affected members who showed either complete or "indeterminate" 9p haplotype sharing without CDKN2A exonic mutations. Reverse transcription polymerase chain reaction and automated sequencing showed expression of both CDKN2A alleles in all family members tested. The lack of CDKN2A promoter mutations and the absence of transcriptional silencing in the germ line of this cohort of families suggest that mutations in the promoter and 5' UTR play a very limited role in melanoma predisposition.

Quality of life and genetics in men with prostate cancer

As family history has been established as a risk factor for prostate cancer, attempts have been made to isolate predisposing genetic variants that are related to hereditary prostate cancer. With many genetic variants still to be identified and investigated, it is not yet possible to fully understand the impact of genetic variants on prostate cancer development. The high survival rates among men with prostate cancer have meant that other issues, such as quality of life (QoL), have also become important. Through their effect on a person’s health, a range of inherited genetic variants may potentially influence QoL in men with prostate cancer, even prior to treatment. Until now, limited research has been conducted on the relationship between genetics and QoL. Thus, this study contributes to an emerging field by aiming to identify certain genetic variants related to the QoL found in men with prostate cancer. It is hoped that this study may lead to future research that will identify men who have an increased risk of a poor QoL following prostate cancer treatment, which will aid in developing treatments that are individually tailored to support them. Previous studies have established that genetic variants of Vascular Endothelial Growth Factor (VEGF) and Insulin-like Growth Factor 1 (IGF-1) may play a role in prostate cancer development. VEGF and IGF-1 have also been reported to be associated with QoL in people with ovarian cancer and colorectal cancer, respectively. This study completed a series of secondary analyses using two major data-sets (from 850 men newly diagnosed with prostate cancer, and approximately 550 men from the general Queensland population), in which genetic variants of VEGF and IGF-1 were investigated for associations with prostate cancer susceptibility and QoL. The first aim of this research was to investigate genetic variants in the VEGF and IGF-I gene for an association with the risk of prostate cancer. It was found that one IGF-1 genetic variant (rs35765) had a statistically significant association with prostate cancer (p = 0.04), and one VEGF genetic variant (rs2146323) had a statistically significant association with advanced prostate cancer (p = 0.02). The estimates suggest that carriers of the CA and AA genotype for rs35765 may have a reduced risk of developing prostate cancer (Odds Ratio (OR) = 0.72, 95% Confidence Interval (CI) = 0.55, 0.95, OR = 0.60, 95% CI = 0.26, 1.39, respectively). Meanwhile, carriers of the CA and AA genotype for rs2146323 may be at increased risk of advanced prostate cancer, which was determined by a Gleason score of above 7 (OR = 1.72, 95% CI = 1.12, 2.63, OR = 1.90, 95% CI = 1.08, 3.34, respectively). Utilising the widely used short-form health survey, the SF-36v2, the second aim of this study was to investigate the relationship between prostate cancer and QoL prior to treatment. Assessing QoL at this time-point was important as little research has been conducted to evaluate if prostate cancer affects QoL regardless of treatment. The analyses found that mean SF-36v2 scale scores related to physical health were higher by at least 0.3 Standard Deviations (SD) among men with prostate cancer than the general population comparison group. This difference was considered clinically significant (defined by group differences in mean SF-36v2 scores by at least 0.3 SD). These differences were also statistically significant (p<0.05). Mean QoL scale scores related to mental health were similar between men with prostate cancer and those from the general population comparison group. The third aim of this study was to investigate genetic variants in the VEGF and IGF-1 gene for an association with QoL in prostate cancer patients prior to their treatment. It was essential to evaluate these relationships prior to treatment, before the involvement of these genes was potentially interrupted by treatment. The analyses found that some genetic variants had a small clinically significant association (0.3 SD) to some QoL domains experienced by these men. However, most relationships were not statistically significant (p>0.05). Most of the associations found identified that a small sub-group of men with prostate cancer (approximately 2%) reported, on average, a slightly better QoL than the majority of the prostate cancer patients. The fourth aim of this research was to investigate whether associations between genetic variants in VEGF and IGF-1 and QoL were specific to men with prostate cancer, or were also applicable to the general male population. It was found that twenty out of one-hundred relationships between the genetic variants of VEGF and IGF-1 and QoL health-measures and scales examined differed between these groups. In the majority of the relationships involving VEGF SNPs that differed, a clinically significant difference (0.3 or more SD) between mean scores among the genotype groups in prostate cancer patients was found, while mean scores among men from the general-population comparison group were similar. For example, prostate cancer participants who carried at least one T allele (CT or TT genotype) for rs3024994 had a clinically significant higher (0.3 SD) mean QoL score in terms of the role-physical scale, than participants who carried the CC genotype. This was not seen among men from the general population sample, as the mean score was similar between genotype groups. The opposite was seen in regards to the IGF-1 SNPs examined. Overall, these relationships were not considered to directly impact on the clinical options for men with prostate cancer. As this study utilised secondary data from two separate studies, there are a number of important limitations that should be acknowledged including issues of multiple comparisons, power, and missing or unavailable data. It is recommended that this study be replicated as a better-designed study that takes greater consideration of the many factors involved in prostate cancer and QoL. Investigation into other genetic variants of VEGF or IGF-1 is also warranted, as is consideration of other genes and their relationship with QoL. Through identifying certain genetic variants that have a modest association to prostate cancer, this project adds to the knowledge surrounding VEGF and IGF-1 and their role in prostate cancer susceptibility. Importantly, this project has also introduced the potential role genetics plays in QoL, through investigating the relationships between genetic variants of VEGF and IGF-1 and QoL.

Gait energy volumes and frontal gait recognition using depth images

Gait energy images (GEIs) and its variants form the basis of many recent appearance-based gait recognition systems. The GEI combines good recognition performance with a simple implementation, though it suffers problems inherent to appearance-based approaches, such as being highly view dependent. In this paper, we extend the concept of the GEI to 3D, to create what we call the gait energy volume, or GEV. A basic GEV implementation is tested on the CMU MoBo database, showing improvements over both the GEI baseline and a fused multi-view GEI approach. We also demonstrate the efficacy of this approach on partial volume reconstructions created from frontal depth images, which can be more practically acquired, for example, in biometric portals implemented with stereo cameras, or other depth acquisition systems. Experiments on frontal depth images are evaluated on an in-house developed database captured using the Microsoft Kinect, and demonstrate the validity of the proposed approach.

Structural Systematics of the Anhydrous 1:1 Proton-Transfer Compounds of 3,5-Dinitrosalicylic Acid with Aniline and Monosubstituted Anilines

The crystal structures of the proton-transfer compounds of 3,5-dinitrosalicylic acid (DNSA) with a series of aniline-type Lewis bases [aniline, 2-hydroxyaniline, 2-methoxyaniline, 3-methoxyaniline, 4-fluoroaniline, 4-chloroaniline and 2-aminoaniline] have been determined and their hydrogen-bonding systems analysed. All are anhydrous 1:1 salts: [(C6H8N)+(C7H3N2O7)-], (1), [(C6H8NO)+(C7H3N2O7)-], (2), [(C7H10NO)+(C7H3N2O7)-], (3), [(C7H10NO)+(C7H3N2O7)-], (4), [(C6H7FN)+(C7H3N2O7)-], (5), [(C6H7ClN)+(C7H3N2O7)-], (6), and [(C6H9N2)+(C7H3N2O7)-], (7) respectively. Crystals of 1 and 6 are triclinic, space group P-1 while the remainder are monoclinic with space group either P21/n (2, 4, 5 and 7) or P21 (3). Unit cell dimensions and contents are: for 1, a = 7.2027(17), b = 7.5699(17), c = 12.9615(16) Å, α = 84.464(14), β = 86.387(15), γ = 75.580(14)o, Z = 2; for 2, a = 7.407(3), b = 6.987(3), c = 27.653(11) Å, β = 94.906(7)o, Z = 4; for 3, a = 8.2816(18), b = 23.151(6), c = 3.9338(10), β = 95.255(19)o, Z = 2; for 4, a = 11.209(2), b = 8.7858(19), c = 15.171(3) Å, β = 93.717(4)o, Z = 4; for 5, a = 26.377(3), b = 10.1602(12), c = 5.1384(10) Å, β = 91.996(13)o, Z = 4; for 6, a = 11.217(3), b = 14.156(5), c = 4.860(3) Å, α = 99.10(4), β = 96.99(4), γ = 76.35(2)o, Z = 2; for 7, a = 12.830(4), b = 8.145(3), c = 14.302(4) Å, β = 102.631(6)o, Z = 4. In all compounds at least one primary linear intermolecular N+-H…O(carboxyl) hydrogen-bonding interaction is present which, together with secondary hydrogen bonding results in the formation of mostly two-dimensional network structures, exceptions being with compounds 4 and 5 (one-dimensional) and compound 6 (three-dimensional). In only two cases [compounds 1 and 4], are weak cation-anion or cation-cation π-π interactions found while weak aromatic C-H…O interactions are insignificant. The study shows that all compounds fit the previously formulated classification scheme for primary and secondary interactive modes for proton-transfer compounds of 3,5-dinitrosalicylic acid but there are some unusual variants.

A conditional retransmission enabled transport protocol for real-time networked control systems

Real-time networked control systems (NCSs) over data networks are being increasingly implemented on a massive scale in industrial applications. Along with this trend, wireless network technologies have been promoted for modern wireless NCSs (WNCSs). However, popular wireless network standards such as IEEE 802.11/15/16 are not designed for real-time communications. Key issues in real-time applications include limited transmission reliability and poor transmission delay performance. Considering the unique features of real-time control systems, this paper develops a conditional retransmission enabled transport protocol (CRETP) to improve the delay performance of the transmission control protocol (TCP) and also the reliability performance of the user datagram protocol (UDP) and its variants. Key features of the CRETP include a connectionless mechanism with acknowledgement (ACK), conditional retransmission and detection of ineffective data packets on the receiver side.

Classification of packet contents for malware detection

Many existing schemes for malware detection are signature-based. Although they can effectively detect known malwares, they cannot detect variants of known malwares or new ones. Most network servers do not expect executable code in their in-bound network traffic, such as on-line shopping malls, Picasa, Youtube, Blogger, etc. Therefore, such network applications can be protected from malware infection by monitoring their ports to see if incoming packets contain any executable contents. This paper proposes a content-classification scheme that identifies executable content in incoming packets. The proposed scheme analyzes the packet payload in two steps. It first analyzes the packet payload to see if it contains multimedia-type data (such as . If not, then it classifies the payload either as text-type (such as or executable. Although in our experiments the proposed scheme shows a low rate of false negatives and positives (4.69% and 2.53%, respectively), the presence of inaccuracies still requires further inspection to efficiently detect the occurrence of malware. In this paper, we also propose simple statistical and combinatorial analysis to deal with false positives and negatives.

Ensuring Correctness During Process Configuration via Partner Synthesis

Variants of the same process can be encountered within one organization or across different organizations. For example, different municipalities, courts, and rental agencies all need to support highly similar processes. In fact, procurement and sales processes can be found in almost any organization. However, despite these similarities, there is also the need to allow for local variations in a controlled manner. Therefore, many academics and practitioners have advocated the use of configurable process models (sometimes referred to as reference models). A configurable process model describes a family of similar process models in a given domain. Such a model can be configured to obtain a specific process model that is subsequently used to handle individual cases, for instance, to process customer orders. Process configuration is notoriously difficult as there may be all kinds of interdependencies between configuration decisions. In fact, an incorrect configuration may lead to behavioral issues such as deadlocks and livelocks. To address this problem, we present a novel verification approach inspired by the “operating guidelines” used for partner synthesis. We view the configuration process as an external service, and compute a characterization of all such services which meet particular requirements via the notion of configuration guideline. As a result, we can characterize all feasible configurations (i. e., configurations without behavioral problems) at design time, instead of repeatedly checking each individual configuration while configuring a process model.

Evidence for CRHR1 in multiple sclerosis using supervised machine learning and meta-analysis in 12 566 individuals

The primary genetic risk factor in multiple sclerosis (MS) is the HLA-DRB1*1501 allele; however, much of the remaining genetic contribution to MS has yet to be elucidated. Several lines of evidence support a role for neuroendocrine system involvement in autoimmunity which may, in part, be genetically determined. Here, we comprehensively investigated variation within eight candidate hypothalamic-pituitary-adrenal (HPA) axis genes and susceptibility to MS. A total of 326 SNPs were investigated in a discovery dataset of 1343 MS cases and 1379 healthy controls of European ancestry using a multi-analytical strategy. Random Forests, a supervised machine-learning algorithm, identified eight intronic SNPs within the corticotrophin-releasing hormone receptor 1 or CRHR1 locus on 17q21.31 as important predictors of MS. On the basis of univariate analyses, six CRHR1 variants were associated with decreased risk for disease following a conservative correction for multiple tests. Independent replication was observed for CRHR1 in a large meta-analysis comprising 2624 MS cases and 7220 healthy controls of European ancestry. Results from a combined meta-analysis of all 3967 MS cases and 8599 controls provide strong evidence for the involvement of CRHR1 in MS. The strongest association was observed for rs242936 (OR = 0.82, 95% CI = 0.74-0.90, P = 9.7 × 10-5). Replicated CRHR1 variants appear to exist on a single associated haplotype. Further investigation of mechanisms involved in HPA axis regulation and response to stress in MS pathogenesis is warranted. © The Author 2010. Published by Oxford University Press. All rights reserved.

The ghrelin signalling system is involved in the consumption of sweets

The gastric-derived orexigenic peptide ghrelin affects brain circuits involved in energy balance as well as in reward. Indeed, ghrelin activates an important reward circuit involved in natural- as well as drug-induced reward, the cholinergic-dopaminergic reward link. It has been hypothesized that there is a common reward mechanism for alcohol and sweet substances in both animals and humans. Alcohol dependent individuals have higher craving for sweets than do healthy controls and the hedonic response to sweet taste may, at least in part, depend on genetic factors. Rat selectively bred for high sucrose intake have higher alcohol consumption than non-sucrose preferring rats and vice versa. In the present study a group of alcohol-consuming individuals selected from a population cohort was investigated for genetic variants of the ghrelin signalling system in relation to both their alcohol and sucrose consumption. Moreover, the effects of GHS-R1A antagonism on voluntary sucrose- intake and operant self-administration, as well as saccharin intake were investigated in preclinical studies using rodents. The effects of peripheral grelin administration on sucrose intake were also examined. Here we found associations with the ghrelin gene haplotypes and increased sucrose consumption, and a trend for the same association was seen in the high alcohol consumers. The preclinical data show that a GHS-R1A antagonist reduces the intake and self-administration of sucrose in rats as well as saccharin intake in mice. Further, ghrelin increases the intake of sucrose in rats. Collectively, our data provide a clear indication that the GHS-R1A antagonists reduces and ghrelin increases the intake of rewarding substances and hence, the central ghrelin signalling system provides a novel target for the development of drug strategies to treat addictive behaviours. © 2011 Landgren et al.

Evaluation of novel Streptococcus pyogenes vaccine candidates incorporating multiple conserved sequences from the C-repeat region of the M-protein

A major challenge for Streptococcus pyogenes vaccine development is the identification of epitopes that confer protection from infection by multiple S. pyogenes M-types. Here we have identified and characterised the distribution of common variant sequences from individual repeat units of the C-repeat region (CRR) of M-proteins representing 77 different M-types. Three polyvalent fusion vaccine candidates (SV1, SV2 and SV3) incorporating the most common variants were subsequently expressed and purified, and demonstrated to be alpha-helical by Circular Dichroism (CD), a secondary conformational characteristic of the CRR in the M-protein. Antibodies raised against each of these constructs recognise M-proteins that vary in their CRR, and bind to the surface of multiple S. pyogenes isolates. Antibodies raised against SV1, containing five variant sequences, also kill heterologous S. pyogenes isolates in in vitro bactericidal assays. Further structural characterisation of this construct demonstrated the conformation of SV1 was stable at different pHs, and thermal unfolding of SV1 a reversible process. Our findings demonstrate that linkage of multiple variant sequences into a single recombinant construct overcomes the need to embed the variant sequences in foreign helix promoting flanking sequences for conformational stability, and demonstrates the viability of the polyvalent candidates as global S. pyogenes vaccine candidates.

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.

Bayesian mixtures for modelling complex medical data : a case study in Parkinson’s disease

Mixture models are a flexible tool for unsupervised clustering that have found popularity in a vast array of research areas. In studies of medicine, the use of mixtures holds the potential to greatly enhance our understanding of patient responses through the identification of clinically meaningful clusters that, given the complexity of many data sources, may otherwise by intangible. Furthermore, when developed in the Bayesian framework, mixture models provide a natural means for capturing and propagating uncertainty in different aspects of a clustering solution, arguably resulting in richer analyses of the population under study. This thesis aims to investigate the use of Bayesian mixture models in analysing varied and detailed sources of patient information collected in the study of complex disease. The first aim of this thesis is to showcase the flexibility of mixture models in modelling markedly different types of data. In particular, we examine three common variants on the mixture model, namely, finite mixtures, Dirichlet Process mixtures and hidden Markov models. Beyond the development and application of these models to different sources of data, this thesis also focuses on modelling different aspects relating to uncertainty in clustering. Examples of clustering uncertainty considered are uncertainty in a patient’s true cluster membership and accounting for uncertainty in the true number of clusters present. Finally, this thesis aims to address and propose solutions to the task of comparing clustering solutions, whether this be comparing patients or observations assigned to different subgroups or comparing clustering solutions over multiple datasets. To address these aims, we consider a case study in Parkinson’s disease (PD), a complex and commonly diagnosed neurodegenerative disorder. In particular, two commonly collected sources of patient information are considered. The first source of data are on symptoms associated with PD, recorded using the Unified Parkinson’s Disease Rating Scale (UPDRS) and constitutes the first half of this thesis. The second half of this thesis is dedicated to the analysis of microelectrode recordings collected during Deep Brain Stimulation (DBS), a popular palliative treatment for advanced PD. Analysis of this second source of data centers on the problems of unsupervised detection and sorting of action potentials or "spikes" in recordings of multiple cell activity, providing valuable information on real time neural activity in the brain.