Multiscale analysis of protein functions and stochastic modelling of gene transcriptional regulatory networks

Genomic and proteomic analyses have attracted a great deal of interests in biological research in recent years. Many methods have been applied to discover useful information contained in the enormous databases of genomic sequences and amino acid sequences. The results of these investigations inspire further research in biological fields in return. These biological sequences, which may be considered as multiscale sequences, have some specific features which need further efforts to characterise using more refined methods. This project aims to study some of these biological challenges with multiscale analysis methods and stochastic modelling approach. The first part of the thesis aims to cluster some unknown proteins, and classify their families as well as their structural classes. A development in proteomic analysis is concerned with the determination of protein functions. The first step in this development is to classify proteins and predict their families. This motives us to study some unknown proteins from specific families, and to cluster them into families and structural classes. We select a large number of proteins from the same families or superfamilies, and link them to simulate some unknown large proteins from these families. We use multifractal analysis and the wavelet method to capture the characteristics of these linked proteins. The simulation results show that the method is valid for the classification of large proteins. The second part of the thesis aims to explore the relationship of proteins based on a layered comparison with their components. Many methods are based on homology of proteins because the resemblance at the protein sequence level normally indicates the similarity of functions and structures. However, some proteins may have similar functions with low sequential identity. We consider protein sequences at detail level to investigate the problem of comparison of proteins. The comparison is based on the empirical mode decomposition (EMD), and protein sequences are detected with the intrinsic mode functions. A measure of similarity is introduced with a new cross-correlation formula. The similarity results show that the EMD is useful for detection of functional relationships of proteins. The third part of the thesis aims to investigate the transcriptional regulatory network of yeast cell cycle via stochastic differential equations. As the investigation of genome-wide gene expressions has become a focus in genomic analysis, researchers have tried to understand the mechanisms of the yeast genome for many years. How cells control gene expressions still needs further investigation. We use a stochastic differential equation to model the expression profile of a target gene. We modify the model with a Gaussian membership function. For each target gene, a transcriptional rate is obtained, and the estimated transcriptional rate is also calculated with the information from five possible transcriptional regulators. Some regulators of these target genes are verified with the related references. With these results, we construct a transcriptional regulatory network for the genes from the yeast Saccharomyces cerevisiae. The construction of transcriptional regulatory network is useful for detecting more mechanisms of the yeast cell cycle.

Expression of PSA-RP2, an alternatively spliced variant from the PSA gene, is increased in prostate cancer tissues but the protein is not secreted from prostate cancer cells

PSA-RP2 is a variant transcript expressed from the PSA gene that is conserved in gorillas, chimpanzees and humans suggesting a particular relevance for this transcript in these primates. We demonstrated by qRT-PCR that PSA-RP2 is upregulated in prostate cancer compared with benign prostatic hyperplasia tissues. The PSA-RP2 protein was not detected in seminal fluid and was cytoplasmically localised but not secreted from LNCaP or transfected PC3 prostate cells, despite secretion from transfected Cos-7 and HEK293 kidney cell lines. PSA-RP2-transfected PC3 cells showed slightly decreased proliferation and increased migration towards PC3-conditioned medium that could suggest a functional role in prostate cancer.

The adenine-rich tract in the 5ʹ-end of the hepatitis C virus ORF encodes a peptide regulating the binding of the C protein to RNA

Hepatitis C virus (HCV ) core (C) protein is thought to bind to viral RNA before it undergoes oligomerization leading to RNA encapsidation. Details of these events are so far unknown. The 5ʹ-terminal C protein coding sequence that includes an adenine (A)-rich tract is a part of an internal ribosome entry site(IRES). This nucleotide sequence but not the corresponding protein sequence is needed for proper initiation of translation of viral RNA by an IRES-dependent mechanism. In this study, we examined the importance of this sequence for the ability of the C protein to bind to viral RNA. Serially truncated C proteins with deletions from 10 up to 45 N-terminal amino acids were expressed in Escherichia coli, purified and tested for binding to viral RNA by a gel shift assay. The results showed that truncation of the C protein from its N-terminus by more than 10 amino acids abolished almost completely its expression in E. coli. The latter could be restored by adding a tag to the N-terminus of the protein. The tagged proteins truncated by 15 or more amino acids showed an anomalous migration in SDS-PAGE. Truncation by more than 20 amino acids resulted in a complete loss of ability of tagged C protein to bind to viral RNA. These results provide clues to the early events in the C protein - RNA interactions leading to C protein oligomerization, RNA encapsidation and virion assembly.

An investigation of the C77G and C772T variations within the human protein tyrosine phosphatase receptor type C gene for association with multiple sclerosis in an Australian population

Multiple sclerosis (MS) is a common cause of neurological disability in young adults. The disease generally manifests in early to middle adulthood and causes various neurological deficits. Autoreactive T lymphocytes and their associated antigens have long been presumed important features of MS pathogenesis. The Protein tyrosine phosphatase receptor type C gene (PTPRC) encodes the T-cell receptor CD45. Variations within PTPRC have been previously associated with diseases of autoimmune origin such as type 1 diabetes mellitus and Graves' disease. We set out to investigate two variants within the PTPRC gene, C77G and C772T in subjects with MS and matched healthy controls to determine whether significant differences exist in these markers in an Australian population. We employed high resolution melt analysis (HRM) and restriction length polymorphism (RFLP) techniques to determine genotypic and allelic frequencies. Our study found no significant difference between frequencies for PTPRC C77G by either genotype (Χ2 = 0.65, P = 0.72) or allele (Χ2 = 0.48, P = 0.49). Similarly, we did not find evidence to suggest an association between PTPRC C772T by genotype (Χ2 = 1.06, P = 0.59) or allele (Χ2 = 0.20, P = 0.66). Linkage disequilibrium (LD) analysis showed strong linkage disequilibrium between the two tested markers (D' = 0.9970, SD = 0.0385). This study reveals no evidence to suggest that these markers are associated with MS in the tested Australian Caucasian population. Although the PTPRC gene has a significant role in regulating CD4+ and CD8+ autoreactive T-cells, interferon-beta responsiveness, and potentially other important processes, our study does not support a role for the two tested variants of this gene in MS susceptibility in the Australian population.

G-protein β3 subunit gene (GNB3) variant in causation of essential hypertension

Essential hypertensives display enhanced signal transduction through pertussis toxin-sensitive G proteins. The T allele of a C825T variant in exon 10 of the G protein β3 subunit gene (GNB3) induces formation of a splice variant (Gβ3-s) with enhanced activity. The T allele of GNB3 was shown recently to be associated with hypertension in unselected German patients (frequency=0.31 versus 0.25 in control). To confirm and extend this finding in a different setting, we performed an association study in Australian white hypertensives. This involved an extensively examined cohort of 110 hypertensives, each of whom were the offspring of 2 hypertensive parents, and 189 normotensives whose parents were both normotensive beyond age 50 years. Genotyping was performed by polymerase chain reaction and digestion with BseDI, which either cut (C allele) or did not cut (T allele) the 268-bp polymerase chain reaction product. T allele frequency in the hypertensive group was 0.43 compared with 0.25 in the normotensive group (χ2=22; P=0.00002; odds ratio=2.3; 95% CI=1.7 to 3.3). The T allele tracked with higher pretreatment blood pressure: diastolic=105±7, 109±16, and 128±28 mm Hg (mean±SD) for CC, CT, and 7T, respectively (P=0.001 by 1-way ANOVA). Blood pressures were higher in female hypertensives with a T allele (P=0.006 for systolic and 0.0003 for diastolic by ANOVA) than they were in male hypertensives. In conclusion, the present study of a group with strong family history supports a role for a genetically determined, physiologically active splice variant of the G protein β3 subunit gene in the causation of essential hypertension.

Comparison of the coat protein, movement protein and RNA polymerase gene sequences of Australian, Chinese, and American isolates of barley yellow dwarf virus transmitted by Rhopalosiphum padi

Barley yellow dwarf luteovirus-GPV (BYDV-GPV) is a common problem in Chinese wheat crops but is unrecorded elsewhere. A defining characteristic of GPV is its capacity to be transmitted efficiently by both Schizaphis graminum and Rhopaloshiphum padi. This dual aphid species transmission contrasts with those of BYDV-RPV and BYDV-SGV, globally distributed viruses, which are efficiently transmitted only by Rhopaloshiphum padi and Schizaphis graminum respectively. The viral RNA sequences encoding the coat protein (22K) gene, the movement protein (17K) gene, the region surrounding the conserved GDD motif of the polymerase gene and the intergenic sequences between these genes were determined for GPV and an Australian isolate of BYDV-RPV (RPVa). In all three genes, the sequences of GPV and RPVa were more similar to those of an American isolate of BYDV-RPV (RPVu) than to any other luteovirus for which there is data available. RPVa and RPVu were very similar, especially their coat proteins which had 97% identity at the amino acid level. The coat protein of GPV had 76% and 78% amino acid identity with RPVa and RPVu respectively. The data suggest that RPVu and RPVa are correctly named as strains of the same serotype and that GPV is sufficiently different from either RPV strain to be considered a distinct BYDV type. The coat protein and movement protein genes of GPV are very dissimilar to SGV. The polymerase sequences of RPVu, RPVa and GPV show close affinities with those of the sobemo-like luteoviruses and little similarity with those of the carmo-like luteoviruses. The sequences of the coat proteins, movement proteins and the polymerase segments of BYDV serotypes, other than RPV and GPV, form a cluster that is separate from their counterpart sequences from dicot-infecting luteoviruses. The RPV and GPV isolates consistently fall within a dicot-infecting cluster. This suggests that RPV and GPV evolved from within this group of viruses. Since these other viruses all infect dicots it seems likely that their common ancestor infected a dicot and that RPV and GPV evolved from a virus that switched hosts from a dicot to a monocot.

Reduced excision repair cross-complementing 1 expression associates with enhanced papilloma formation and fibroblast transformation after genetic disruption of secreted protein acidic and rich in cysteine

SPARC (secreted protein acidic and rich in cysteine)/ osteonectin/BM-40 is a matricellular protein implicated in development, cell transformation and tumorigenesis. We have examined the role of SPARC in cell transformation induced chemically with 7,12-dimethylbenz[a]anthracene (DMBA) and 12- tetradecanoylphorbol-13-acetate (TPA) in embryonic fibroblasts and in the skin of mice. Embryonic fibroblasts from SPARCnull mice showed increases in cell proliferation, enhanced sensitivity to DMBA and a higher number of DMBA/TPA-induced transformation foci. The number of DMBA-DNA adducts was 9 times higher in SPARCnull fibroblasts and their stability was lower than wild-type fibroblasts, consistent with a reduction in excision repair cross-complementing 1 the nucleotide excision repair enzyme in these cells. The SPARCnull mice showed an increase in both the speed and number of papillomas forming after topical administration of DMBA/TPA to the skin. These papillomas showed reduced growth and reduced progression to a more malignant phenotype, indicating that the effect of SPARC on tumorigenesis depends upon the transformation stage and/or tissue context. These data reinforce a growing number of observations in which SPARC has shown opposite effects on different tumor types/stages.

UafB is a serine-rich repeat adhesin of Staphylococcus saprophyticus that mediates binding to fibronectin, fibrinogen and human uroepithelial cells

Staphylococcus saprophyticus is an important cause of urinary tract infection (UTI), particularly among young women, and is second only to uropathogenic Escherichia coli as the most frequent cause of UTI. The molecular mechanisms of urinary tract colonization by S. saprophyticus remain poorly understood. We have identified a novel 6.84 kb plasmid-located adhesin-encoding gene in S. saprophyticus strain MS1146 which we have termed uro-adherence factor B (uafB). UafB is a glycosylated serine-rich repeat protein that is expressed on the surface of S. saprophyticus MS1146. UafB also functions as a major cell surface hydrophobicity factor. To characterize the role of UafB we generated an isogenic uafB mutant in S. saprophyticus MS1146 by interruption with a group II intron. The uafB mutant had a significantly reduced ability to bind to fibronectin and fibrinogen. Furthermore, we show that a recombinant protein containing the putative binding domain of UafB binds specifically to fibronectin and fibrinogen. UafB was not involved in adhesion in a mouse model of UTI; however, we observed a striking UafB-mediated adhesion phenotype to human uroepithelial cells. We have also identified genes homologous to uafB in other staphylococci which, like uafB, appear to be located on transposable elements. Thus, our data indicate that UafB is a novel adhesin of S. saprophyticus that contributes to cell surface hydrophobicity, mediates adhesion to fibronectin and fibrinogen, and exhibits tropism for human uroepithelial cells.

Microparticle-mediated gene delivery for the enhanced expression of a 19-KDa fragment of merozoite surface protein 1 of Plasmodium falciparum

The 19 kDa carboxyl-terminal fragment of merozoite surface protein 1 (MSP119) is a major component of the invasion-inhibitory response in individual immunity to malaria. A novel ultrasonic atomization approach for the formulation of biodegradable poly(lactic-co-glycolic acid) (PLGA) microparticles of malaria DNA vaccines encoding MSP119 is presented here. After condensing the plasmid DNA (pDNA) molecules with a cationic polymer polyethylenimine (PEI), a 40 kHz ultrasonic atomization frequency was used to formulate PLGA microparticles at a flow rate of 18 mL h1. High levels of gene expression and moderate cytotoxicity in COS-7 cells were achieved with the condensed pDNA at a nitrogen to phosphate (N/P) ratio of 20, thus demonstrating enhanced cellular uptake and expression of the transgene. The ability of the microparticles to convey pDNA was examined by characterizing the formulated microparticles. The microparticles displayed Z-average hydrodynamic diameters of 1.50-2.10 lm and zeta potentials of 17.8-23.2 mV. The encapsulation efficiencies were between 78 and 83%, and 76 and 85% of the embedded malaria pDNA molecules were released under physiological conditions in vitro. These results indicate that PLGA-mediated microparticles can be employed as potential gene delivery systems to antigen-presenting cells in the prevention of malaria.

A novel ACVR1 mutation in the glycine/serine-rich domain found in the most benign case of a fibrodysplasia ossificans progressiva variant reported to date

Fibrodysplasia Ossificans Progressiva (FOP) is a rare, autosomal dominant condition, classically characterised by heterotopic ossification beginning in childhood and congenital great toe malformations; occurring in response to a c.617 G>A ACVR1 mutation in the functionally important glycine/serine-rich domain of exon 6. Here we describe a novel c.587 T>C mutation in the glycine/serine-rich domain of ACVR1, associated with delayed onset of heterotopic ossification and an exceptionally mild clinical course. Absence of great toe malformations, the presence of early ossification of the cervical spine facets joints, plus mild bilateral camptodactyly of the 5th fingers, together with a novel ACVR1 mutation, are consistent with the 'FOP-variant' syndrome. The c.587 T>C mutation replaces a conserved leucine with proline at residue 196. Modelling of the mutant protein reveals a steric clash with the kinase domain that will weaken interactions with FKBP12 and induce exposure of the glycine/serine-rich repeat. The mutant receptor is predicted to be hypersensitive to ligand stimulation rather than being constitutively active, consistent with the mild clinical phenotype. This case extends our understanding of the 'FOP-variant' syndrome.

G-protein beta3 subunit gene (GNB3) variant in causation of essential hypertension

-Essential hypertensives display enhanced signal transduction through pertussis toxin-sensitive G proteins. The T allele of a C825T variant in exon 10 of the G protein beta3 subunit gene (GNB3) induces formation of a splice variant (Gbeta3-s) with enhanced activity. The T allele of GNB3 was shown recently to be associated with hypertension in unselected German patients (frequency=0.31 versus 0.25 in control). To confirm and extend this finding in a different setting, we performed an association study in Australian white hypertensives. This involved an extensively examined cohort of 110 hypertensives, each of whom were the offspring of 2 hypertensive parents, and 189 normotensives whose parents were both normotensive beyond age 50 years. Genotyping was performed by polymerase chain reaction and digestion with BseDI, which either cut (C allele) or did not cut (T allele) the 268-bp polymerase chain reaction product. T allele frequency in the hypertensive group was 0.43 compared with 0.25 in the normotensive group (chi2=22; P=0.00002; odds ratio=2.3; 95% CI=1.7 to 3.3). The T allele tracked with higher pretreatment blood pressure: diastolic=105+/-7, 109+/-16, and 128+/-28 mm Hg (mean+/-SD) for CC, CT, and TT, respectively (P=0.001 by 1-way ANOVA). Blood pressures were higher in female hypertensives with a T allele (P=0.006 for systolic and 0.0003 for diastolic by ANOVA) than they were in male hypertensives. In conclusion, the present study of a group with strong family history supports a role for a genetically determined, physiologically active splice variant of the G protein beta3 subunit gene in the causation of essential hypertension.