Negative second virial coefficients as predictors of protein crystal growth: Evidence from sedimentation equilibrium studies that refutes the designation of those light scattering parameters as osmotic virial coefficients

The effects of ammonium sulphate concentration on the osmotic second virial coefficient (B-AA/M-A) for equine serum albumin (pH 5.6, 20 degrees C) have been examined by sedimentation equilibrium. After an initial steep decrease with increasing ammonium sulphate concentration, B-AA/M-A assumes an essentially concentration-independent magnitude of 8-9 ml/g. Such behaviour conforms with the statistical-mechanical prediction that a sufficient increase in ionic strength should effectively eliminate the contributions of charge interactions to B-AA/M-A but have no effect on the covolume contribution (8.4 ml/g for serum albumin). A similar situation is shown to apply to published sedimentation equilibrium data for lysozyme (pH 4.5). Although termed osmotic second virial coefficients and designated as such (B-22), the negative values obtained in published light scattering studies of both systems have been described incorrectly because of the concomitant inclusion of the protein-salt contribution to thermodynamic nonideality of the protein. Those negative values are still valid predictors of conditions conducive to crystal growth inasmuch as they do reflect situations in which there is net attraction between protein molecules. However, the source of attraction responsible for the negative virial coefficient stems from the protein-salt rather than the protein-protein contribution, which is necessarily positive. (c) 2005 Elsevier B.V. All rights reserved.

Peptide surfactants (Pepfactants) for switchable foams and emulsions

We have developed a series of sustainable peptide surfactants (Pepfactants) capable of stabilizing foams and emulsions in a stimuli-responsive manner, based on the reversible formation of a mechanically strong interfacial film. Under conditions where the interfacially adsorbed peptide forms a mechanically strong film state, foam or emulsion stabilization occurs as a direct result of the film strength. Under conditions where the interfacially adsorbed peptide forms a mobile detergent state, foam or emulsion stabilization is either reduced, or does not occur. Preformed foams or emulsions stabilized by Pepfactants undergo rapid phase coalescence when the film state is converted to the detergent state. Switching between film and detergent states is readily and reversibly achieved by a change in the bulk solution composition, such as a change in pH, or the addition or sequestering of metal ions. Copyright © 2007 Curtin University of Technology and John Wiley & Sons, Ltd.

Mutations of beta-arrestin 2 that limit self-association also interfere with interactions with the beta2-adrenoceptor and the ERK1/2 MAPKs:implications for beta2-adrenoceptor signalling via the ERK1/2 MAPKs

FRET (fluorescence resonance energy transfer) and co-immunoprecipitation studies confirmed the capacity of beta-arrestin 2 to self-associate. Amino acids potentially involved in direct protein-protein interaction were identified via combinations of spot-immobilized peptide arrays and mapping of surface exposure. Among potential key amino acids, Lys(285), Arg(286) and Lys(295) are part of a continuous surface epitope located in the polar core between the N- and C-terminal domains. Introduction of K285A/R286A mutations into beta-arrestin 2-eCFP (where eCFP is enhanced cyan fluorescent protein) and beta-arrestin 2-eYFP (where eYFP is enhanced yellow fluorescent protein) constructs substantially reduced FRET, whereas introduction of a K295A mutation had a more limited effect. Neither of these mutants was able to promote beta2-adrenoceptor-mediated phosphorylation of the ERK1/2 (extracellular-signal-regulated kinase 1/2) MAPKs (mitogen-activated protein kinases). Both beta-arrestin 2 mutants displayed limited capacity to co-immunoprecipitate ERK1/2 and further spot-immobilized peptide arrays indicated each of Lys(285), Arg(286) and particularly Lys(295) to be important for this interaction. Direct interactions between beta-arrestin 2 and the beta2-adrenoceptor were also compromised by both K285A/R286A and K295A mutations of beta-arrestin 2. These were not non-specific effects linked to improper folding of beta-arrestin 2 as limited proteolysis was unable to distinguish the K285A/R286A or K295A mutants from wild-type beta-arrestin 2, and the interaction of beta-arrestin 2 with JNK3 (c-Jun N-terminal kinase 3) was unaffected by the K285A/R286A or L295A mutations. These results suggest that amino acids important for self-association of beta-arrestin 2 also play an important role in the interaction with both the beta2-adrenoceptor and the ERK1/2 MAPKs. Regulation of beta-arrestin 2 self-association may therefore control beta-arrestin 2-mediated beta2-adrenoceptor-ERK1/2 MAPK signalling.

Gold nanoparticles decorated with oligo(ethylene glycol) thiols:protein resistance and colloidal stability

The interactions between proteins and gold colloids functionalized with protein-resistant oligo(ethylene glycol) (OEG) thiol, HS(CH(2))(11) (OCH(2)CH(2))(6)OMe (EG(6)OMe), in aqueous solution have been studied by small-angle X-ray scattering (SAXS) and UV-vis spectroscopy. The mean size, 2R, and the size distribution of the decorated gold colloids have been characterized by SAXS. The monolayer-protected gold colloids have no correlations due to the low volume fraction in solution and are stable in a wide range of temperatures (5-70 degrees C, pH (1.3-12.4), and ionic strength (0-1.0 M). In contrast, protein (bovine serum albumin) solutions with concentrations in the range of 60-200 mg/mL (4.6-14.5 vol show a pronounced correlation peak in SAXS, which results from the repulsive electrostatic interaction between charged proteins. These protein interactions show significant dependence on ionic strength, as would be expected for an electrostatic interaction (Zhang et al. J. Phys. Chem. B 2007, 111, 251). For a mixture of proteins and gold colloids, the protein-protein interaction changes little upon mixing with OEG-decorated gold colloids. In contrast, the colloid-colloid interaction is found to be strongly dependent on the protein concentration and the size of the colloid itself. Adding protein to a colloidal solution results in an attractive depletion interaction between functionalized gold colloids, and above a critical protein concentration, c*, the colloids form aggregates and flocculate. Adding salt to such mixtures enhances the depletion effect and decreases the critical protein concentration. The aggregation is a reversible process (i.e., diluting the solution leads to dissolution of aggregates). The results also indicate that the charge of the OEG self-assembled monolayer at a curved interface has a rather limited effect on the colloidal stabilization and the repulsive interaction with proteins.

Protein interactions studied by SAXS:effect of ionic strength and protein concentration for BSA in aqueous solutions

We have studied a series of samples of bovine serum albumin (BSA) solutions with protein concentration, c, ranging from 2 to 500 mg/mL and ionic strength, I, from 0 to 2 M by small-angle X-ray scattering (SAXS). The scattering intensity distribution was compared to simulations using an oblate ellipsoid form factor with radii of 17 x 42 x 42 A, combined with either a screened Coulomb, repulsive structure factor, S-SC(q), or an attractive square-well structure factor, S-SW(q). At pH = 7, BSA is negatively charged. At low ionic strength, I <0.3 M, the total interaction exhibits a decrease of the repulsive interaction when compared to the salt-free solution, as the net surface charge is screened, and the data can be fitted by assuming an ellipsoid form factor and screened Coulomb interaction. At moderate ionic strength (0.3-0.5 M), the interaction is rather weak, and a hard-sphere structure factor has been used to simulate the data with a higher volume fraction. Upon further increase of the ionic strength (I >= 1.0 M), the overall interaction potential was dominated by an additional attractive potential, and the data could be successfully fitted by an ellipsoid form factor and a square-well potential model. The fit parameters, well depth and well width, indicate that the attractive potential caused by a high salt concentration is weak and long-ranged. Although the long-range, attractive potential dominated the protein interaction, no gelation or precipitation was observed in any of the samples. This is explained by the increase of a short-range, repulsive interaction between protein molecules by forming a hydration layer with increasing salt concentration. The competition between long-range, attractive and short-range, repulsive interactions accounted for the stability of concentrated BSA solution at high ionic strength.

Peptide binding to the HLA-DRB1 supertype:a proteochemometrics analysis

A proteochemometrics approach was applied to a set of 2666 peptides binding to 12 HLA-DRB1 proteins. Sequences of both peptide and protein were described using three z-descriptors. Cross terms accounting for adjacent positions and for every second position in the peptides were included in the models, as well as cross terms for peptide/protein interactions. Models were derived based on combinations of different blocks of variables. These models had moderate goodness of fit, as expressed by r2, which ranged from 0.685 to 0.732; and good cross-validated predictive ability, as expressed by q2, which varied from 0.678 to 0.719. The external predictive ability was tested using a set of 356 HLA-DRB1 binders, which showed an r2(pred) in the range 0.364-0.530. Peptide and protein positions involved in the interactions were analyzed in terms of hydrophobicity, steric bulk and polarity.

Histidine hydrogen bonding in MHC at pH 5 and pH 7 modeled by molecular docking and molecular dynamics simulations

Hydrogen bonds play important roles in maintaining the structure of proteins and in the formation of most biomolecular protein-ligand complexes. All amino acids can act as hydrogen bond donors and acceptors. Among amino acids, Histidine is unique, as it can exist in neutral or positively charged forms within the physiological pH range of 5.0 to 7.0. Histidine can thus interact with other aromatic residues as well as forming hydrogen bonds with polar and charged residues. The ability of His to exchange a proton lies at the heart of many important functional biomolecular interactions, including immunological ones. By using molecular docking and molecular dynamics simulation, we examine the influence of His protonation/deprotonation on peptide binding affinity to MHC class II proteins from locus HLA-DP. Peptide-MHC interaction underlies the adaptive cellular immune response, upon which the next generation of commercially-important vaccines will depend. Consistent with experiment, we find that peptides containing protonated His residues bind better to HLA-DP proteins than those with unprotonated His. Enhanced binding at pH 5.0 is due, in part, to additional hydrogen bonds formed between peptide His+ and DP proteins. In acidic endosomes, protein His79β is predominantly protonated. As a result, the peptide binding cleft narrows in the vicinity of His79β, which stabilizes the peptide - HLA-DP protein complex. © 2014 Bentham Science Publishers.

Efeitos in vivo e in vitro de agonistas parciais do receptor NOP: implicações para o tratamento da ansiedade, depressão e mania

Introduction: This study aimed to investigate the effects of the two peptide NOP partial agonists (UFP-113 and [F/G]N/OFQ(1-13)NH2) and the non peptide NOP partial agonist (AT-090) in the mouse emotional behavior as well as in the intracellular transduction pathways following the receptor binding. Methods: Male Swiss or CD-1 mice were used in this study together with NOP(+/+) and NOP(-/-) mice. The elevated plus maze (EPM) was used to evaluate the effects of compounds on anxiety-like behaviors. Diazepam and the NOP agonists, N/OFQ and Ro 65-6570, were used as positive controls in the EPM. NOP(+/+) and NOP(-/-) mice were used to evaluate the selectivity of those compounds that induced anxiolytic-like behaviors. The forced swim test (FST) was used to evaluate the effects of compounds on depressive-like behaviors. Nortriptyline and the NOP antagonists, UFP-101 and SB-612111, were used as positive controls in the FST. The effects of N/OFQ, UFP-101, SB-612111, UFP-113, [F/G]N/OFQ(1-13)NH2, and AT-090 were assessed in the methylphenidate-induced hyperlocomotion (MIH) test; in this assay valproate was used as positive control. The G protein and β-arrestin 2 transduction pathways of NOP receptor agonists (N/OFQ and Ro 65-6570), antagonist (UFP-101), and partial agonists (UFP-113, [F/G]N/OFQ(1-13)NH2, and AT-090) were also evaluated using an innovative assay that measures a bioluminescence resonance energy transfer process. For this, cell lines permanently co-expressing the NOP receptor coupled to luciferase (energy donor), and green fluorescent protein (energy acceptor) coupled to one of the effector proteins (G protein or β-arrestin 2) were used. Results: Diazepam (1 mg/kg), N/OFQ (1 nmol), Ro 65-6570 (0.1 mg/kg), and AT-090 (0.01 mg/kg) induced anxiolytic-like effect in mice in the EPM. The effects of Ro 65-6570 and AT-090 were selective to NOP receptor. UFP-113 (0.01-1 nmol) and [F/G]N/OFQ(1-13)NH2 (0.1-3 nmol) were inactive in the EPM. In the FST, nortriptyline (30 mg/kg), UFP-101 (10 nmol), SB-612111 (10 mg/kg), UFP-113 (0.01 and 0.1 nmol), and [F/G]N/OFQ(1-13)NH2 (0.3 and 1 nmol) induced antidepressant-like effects, while AT-090 (0.001-0.1 mg/kg) was inactive in this assay. The effects of UFP-113 and [F/G]N/OFQ(1-13)NH2 were selective to NOP receptor. Valproate (400 mg/kg) counteracted methylphenidate (MPH, 10 mg/kg)-induced hyperlocomotion in mice in the open field. N/OFQ (1 nmol), UFP-113 (0.01-0.1 nmol), and [F/G]N/OFQ(1-13)NH2 (1 nmol) were also able to reduce the MPH-induced hyperlocomotion, without changing the locomotor activity per se. The effect of UFP-113 was selective to NOP receptor. The UFP-101 (10 nmol), SB-612111 (10 mg/kg), and AT-090 (0.001-0.03 mg/kg) did not change the hyperlocomotor effect of methylphenidate. In vitro, N/OFQ and Ro 65-6570 behaved as NOP full agonists for G-protein and β-arrestin 2 pathways. AT-090 behaved as NOP receptor partial agonist for both transduction pathways, while UFP-113 and [F/G]N/OFQ(1-13)NH2 behaved as partial agonists and antagonists of NOP receptor for NOP/G protein and NOP/β-arrestin 2, respectively. UFP-101 behaved as NOP receptor antagonist for both transduction pathways. Conclusion: NOP ligands producing same effects on NOP/G protein interaction (partial agonism), but with opposite effects on β-arrestin 2 recruitment (partial agonism vs antagonism), can promote different in vivo effects on anxiety and mood as it was observed in the behavioral tests. This work corroborates the potential of NOP receptor as an innovative pharmacological target for the treatment of emotional disorders.

Efeitos in vivo e in vitro de agonistas parciais do receptor NOP: implicações para o tratamento da ansiedade, depressão e mania

Introduction: This study aimed to investigate the effects of the two peptide NOP partial agonists (UFP-113 and [F/G]N/OFQ(1-13)NH2) and the non peptide NOP partial agonist (AT-090) in the mouse emotional behavior as well as in the intracellular transduction pathways following the receptor binding. Methods: Male Swiss or CD-1 mice were used in this study together with NOP(+/+) and NOP(-/-) mice. The elevated plus maze (EPM) was used to evaluate the effects of compounds on anxiety-like behaviors. Diazepam and the NOP agonists, N/OFQ and Ro 65-6570, were used as positive controls in the EPM. NOP(+/+) and NOP(-/-) mice were used to evaluate the selectivity of those compounds that induced anxiolytic-like behaviors. The forced swim test (FST) was used to evaluate the effects of compounds on depressive-like behaviors. Nortriptyline and the NOP antagonists, UFP-101 and SB-612111, were used as positive controls in the FST. The effects of N/OFQ, UFP-101, SB-612111, UFP-113, [F/G]N/OFQ(1-13)NH2, and AT-090 were assessed in the methylphenidate-induced hyperlocomotion (MIH) test; in this assay valproate was used as positive control. The G protein and β-arrestin 2 transduction pathways of NOP receptor agonists (N/OFQ and Ro 65-6570), antagonist (UFP-101), and partial agonists (UFP-113, [F/G]N/OFQ(1-13)NH2, and AT-090) were also evaluated using an innovative assay that measures a bioluminescence resonance energy transfer process. For this, cell lines permanently co-expressing the NOP receptor coupled to luciferase (energy donor), and green fluorescent protein (energy acceptor) coupled to one of the effector proteins (G protein or β-arrestin 2) were used. Results: Diazepam (1 mg/kg), N/OFQ (1 nmol), Ro 65-6570 (0.1 mg/kg), and AT-090 (0.01 mg/kg) induced anxiolytic-like effect in mice in the EPM. The effects of Ro 65-6570 and AT-090 were selective to NOP receptor. UFP-113 (0.01-1 nmol) and [F/G]N/OFQ(1-13)NH2 (0.1-3 nmol) were inactive in the EPM. In the FST, nortriptyline (30 mg/kg), UFP-101 (10 nmol), SB-612111 (10 mg/kg), UFP-113 (0.01 and 0.1 nmol), and [F/G]N/OFQ(1-13)NH2 (0.3 and 1 nmol) induced antidepressant-like effects, while AT-090 (0.001-0.1 mg/kg) was inactive in this assay. The effects of UFP-113 and [F/G]N/OFQ(1-13)NH2 were selective to NOP receptor. Valproate (400 mg/kg) counteracted methylphenidate (MPH, 10 mg/kg)-induced hyperlocomotion in mice in the open field. N/OFQ (1 nmol), UFP-113 (0.01-0.1 nmol), and [F/G]N/OFQ(1-13)NH2 (1 nmol) were also able to reduce the MPH-induced hyperlocomotion, without changing the locomotor activity per se. The effect of UFP-113 was selective to NOP receptor. The UFP-101 (10 nmol), SB-612111 (10 mg/kg), and AT-090 (0.001-0.03 mg/kg) did not change the hyperlocomotor effect of methylphenidate. In vitro, N/OFQ and Ro 65-6570 behaved as NOP full agonists for G-protein and β-arrestin 2 pathways. AT-090 behaved as NOP receptor partial agonist for both transduction pathways, while UFP-113 and [F/G]N/OFQ(1-13)NH2 behaved as partial agonists and antagonists of NOP receptor for NOP/G protein and NOP/β-arrestin 2, respectively. UFP-101 behaved as NOP receptor antagonist for both transduction pathways. Conclusion: NOP ligands producing same effects on NOP/G protein interaction (partial agonism), but with opposite effects on β-arrestin 2 recruitment (partial agonism vs antagonism), can promote different in vivo effects on anxiety and mood as it was observed in the behavioral tests. This work corroborates the potential of NOP receptor as an innovative pharmacological target for the treatment of emotional disorders.

An allosteric role for receptor activity-modifying proteins in defining GPCR pharmacology

G protein-coupled receptors are allosteric proteins that control transmission of external signals to regulate cellular response. Although agonist binding promotes canonical G protein signalling transmitted through conformational changes, G protein-coupled receptors also interact with other proteins. These include other G protein-coupled receptors, other receptors and channels, regulatory proteins and receptor-modifying proteins, notably receptor activity-modifying proteins (RAMPs). RAMPs have at least 11 G protein-coupled receptor partners, including many class B G protein-coupled receptors. Prototypic is the calcitonin receptor, with altered ligand specificity when co-expressed with RAMPs. To gain molecular insight into the consequences of this protein–protein interaction, we combined molecular modelling with mutagenesis of the calcitonin receptor extracellular domain, assessed in ligand binding and functional assays. Although some calcitonin receptor residues are universally important for peptide interactions (calcitonin, amylin and calcitonin gene-related peptide) in calcitonin receptor alone or with receptor activity-modifying protein, others have RAMP-dependent effects, whereby mutations decreased amylin/calcitonin gene-related peptide potency substantially only when RAMP was present. Remarkably, the key residues were completely conserved between calcitonin receptor and AMY receptors, and between subtypes of AMY receptor that have different ligand preferences. Mutations at the interface between calcitonin receptor and RAMP affected ligand pharmacology in a RAMP-dependent manner, suggesting that RAMP may allosterically influence the calcitonin receptor conformation. Supporting this, molecular dynamics simulations suggested that the calcitonin receptor extracellular N-terminal domain is more flexible in the presence of receptor activity-modifying protein 1. Thus, RAMPs may act in an allosteric manner to generate a spectrum of unique calcitonin receptor conformational states, explaining the pharmacological preferences of calcitonin receptor-RAMP complexes. This provides novel insight into our understanding of G protein-coupled receptor-protein interaction that is likely broadly applicable for this receptor class.

Comparison of module detection algorithms in protein networks and investigation of the biological meaning of predicted modules

Background
It is generally acknowledged that a functional understanding of a biological system can only be obtained by an understanding of the collective of molecular interactions in form of biological networks. Protein networks are one particular network type of special importance, because proteins form the functional base units of every biological cell. On a mesoscopic level of protein networks, modules are of significant importance because these building blocks may be the next elementary functional level above individual proteins allowing to gain insight into fundamental organizational principles of biological cells.
Results
In this paper, we provide a comparative analysis of five popular and four novel module detection algorithms. We study these module prediction methods for simulated benchmark networks as well as 10 biological protein interaction networks (PINs). A particular focus of our analysis is placed on the biological meaning of the predicted modules by utilizing the Gene Ontology (GO) database as gold standard for the definition of biological processes. Furthermore, we investigate the robustness of the results by perturbing the PINs simulating in this way our incomplete knowledge of protein networks.
Conclusions
Overall, our study reveals that there is a large heterogeneity among the different module prediction algorithms if one zooms-in the biological level of biological processes in the form of GO terms and all methods are severely affected by a slight perturbation of the networks. However, we also find pathways that are enriched in multiple modules, which could provide important information about the hierarchical organization of the system

Studies of Protein-Protein and Protein-Water Interactions by Small Angle X-Ray Scattering, Terahertz Spectroscopy, ASMOS, And Computer Simulation

The protein folding problem has been one of the most challenging subjects in biological physics due to its complexity. Energy landscape theory based on statistical mechanics provides a thermodynamic interpretation of the protein folding process. We have been working to answer fundamental questions about protein-protein and protein-water interactions, which are very important for describing the energy landscape surface of proteins correctly. At first, we present a new method for computing protein-protein interaction potentials of solvated proteins directly from SAXS data. An ensemble of proteins was modeled by Metropolis Monte Carlo and Molecular Dynamics simulations, and the global X-ray scattering of the whole model ensemble was computed at each snapshot of the simulation. The interaction potential model was optimized and iterated by a Levenberg-Marquardt algorithm. Secondly, we report that terahertz spectroscopy directly probes hydration dynamics around proteins and determines the size of the dynamical hydration shell. We also present the sequence and pH-dependence of the hydration shell and the effect of the hydrophobicity. On the other hand, kinetic terahertz absorption (KITA) spectroscopy is introduced to study the refolding kinetics of ubiquitin and its mutants. KITA results are compared to small angle X-ray scattering, tryptophan fluorescence, and circular dichroism results. We propose that KITA monitors the rearrangement of hydrogen bonding during secondary structure formation. Finally, we present development of the automated single molecule operating system (ASMOS) for a high throughput single molecule detector, which levitates a single protein molecule in a 10 µm diameter droplet by the laser guidance. I also have performed supporting calculations and simulations with my own program codes.

The bacterial helicase-primase interaction: a common structural/functional module

The lack of a high-resolution structure for the bacterial helicase-primase complex and the fragmented structural information for the individual proteins have been hindering our detailed understanding of this crucial binary protein interaction. Two new structures for the helicase-interacting domain of the bacterial primases from Escherichia coli and Bacillus stearothermophilus have recently been solved and both revealed a unique and surprising structural similarity to the amino-terminal domain of the helicase itself. In this minireview, the current data are discussed and important new structural and functional aspects of the helicase-primase interaction are highlighted. An attractive structural model with direct biological significance for the function of this complex and also for the development of new antibacterial compounds is examined.

Fuzzy methods for analysis of microarrays and networks

Bioinformatics involves analyses of biological data such as DNA sequences, microarrays and protein-protein interaction (PPI) networks. Its two main objectives are the identification of genes or proteins and the prediction of their functions. Biological data often contain uncertain and imprecise information. Fuzzy theory provides useful tools to deal with this type of information, hence has played an important role in analyses of biological data. In this thesis, we aim to develop some new fuzzy techniques and apply them on DNA microarrays and PPI networks. We will focus on three problems: (1) clustering of microarrays; (2) identification of disease-associated genes in microarrays; and (3) identification of protein complexes in PPI networks. The first part of the thesis aims to detect, by the fuzzy C-means (FCM) method, clustering structures in DNA microarrays corrupted by noise. Because of the presence of noise, some clustering structures found in random data may not have any biological significance. In this part, we propose to combine the FCM with the empirical mode decomposition (EMD) for clustering microarray data. The purpose of EMD is to reduce, preferably to remove, the effect of noise, resulting in what is known as denoised data. We call this method the fuzzy C-means method with empirical mode decomposition (FCM-EMD). We applied this method on yeast and serum microarrays, and the silhouette values are used for assessment of the quality of clustering. The results indicate that the clustering structures of denoised data are more reasonable, implying that genes have tighter association with their clusters. Furthermore we found that the estimation of the fuzzy parameter m, which is a difficult step, can be avoided to some extent by analysing denoised microarray data. The second part aims to identify disease-associated genes from DNA microarray data which are generated under different conditions, e.g., patients and normal people. We developed a type-2 fuzzy membership (FM) function for identification of diseaseassociated genes. This approach is applied to diabetes and lung cancer data, and a comparison with the original FM test was carried out. Among the ten best-ranked genes of diabetes identified by the type-2 FM test, seven genes have been confirmed as diabetes-associated genes according to gene description information in Gene Bank and the published literature. An additional gene is further identified. Among the ten best-ranked genes identified in lung cancer data, seven are confirmed that they are associated with lung cancer or its treatment. The type-2 FM-d values are significantly different, which makes the identifications more convincing than the original FM test. The third part of the thesis aims to identify protein complexes in large interaction networks. Identification of protein complexes is crucial to understand the principles of cellular organisation and to predict protein functions. In this part, we proposed a novel method which combines the fuzzy clustering method and interaction probability to identify the overlapping and non-overlapping community structures in PPI networks, then to detect protein complexes in these sub-networks. Our method is based on both the fuzzy relation model and the graph model. We applied the method on several PPI networks and compared with a popular protein complex identification method, the clique percolation method. For the same data, we detected more protein complexes. We also applied our method on two social networks. The results showed our method works well for detecting sub-networks and give a reasonable understanding of these communities.

A highly divergent South African geminivirus species illuminates the ancient evolutionary history of this family

Background. We have characterised a new highly divergent geminivirus species, Eragrostis curvula streak virus (ECSV), found infecting a hardy perennial South African wild grass. ECSV represents a new genus-level geminivirus lineage, and has a mixture of features normally associated with other specific geminivirus genera. Results. Whereas the ECSV genome is predicted to express a replication associated protein (Rep) from an unspliced complementary strand transcript that is most similar to those of begomoviruses, curtoviruses and topocuviruses, its Rep also contains what is apparently a canonical retinoblastoma related protein interaction motif such as that found in mastreviruses. Similarly, while ECSV has the same unusual TAAGATTCC virion strand replication origin nonanucleotide found in another recently described divergent geminivirus, Beet curly top Iran virus (BCTIV), the rest of the transcription and replication origin is structurally more similar to those found in begomoviruses and curtoviruses than it is to those found in BCTIV and mastreviruses. ECSV also has what might be a homologue of the begomovirus transcription activator protein gene found in begomoviruses, a mastrevirus-like coat protein gene and two intergenic regions. Conclusion. Although it superficially resembles a chimaera of geminiviruses from different genera, the ECSV genome is not obviously recombinant, implying that the features it shares with other geminiviruses are those that were probably present within the last common ancestor of these viruses. In addition to inferring how the ancestral geminivirus genome may have looked, we use the discovery of ECSV to refine various hypotheses regarding the recombinant origins of the major geminivirus lineages. © 2009 Varsani et al; licensee BioMed Central Ltd.