Hydrophobicity Plots Research Articles

English

Scientific efforts in the eradication of neglected tropical diseases, such as those caused by the parasitic helminthes, can be improved if a database of key virulence factors directly implicated in pathogenesis is available. As a first step towards creating SchistoVir, a database of virulence protein factors in schistosomes, in this study, we curated, annotated and aligned sequences of twenty virulence factors identified from the literature, using several bioinformatics tools including UniProtKB, SchistoDB, VirulentPred, InterProScan, ProtScale, MotifScan, TDRtarget, SignalP, MODBASE, PDB and MUSCLE. Among the protein entries, the most frequently occurring amino acid residues were lysine, serine, leucine, glutamine, glycine and cysteine in order of magnitude. Although sequence repeat regions (SRRs) of significant value were identified manually in fifty percent of the proteins (while dipeptide repeats (DiPs) and single amino acid repeats (SAARs) were not), nevertheless, seventy-two percent of the protein entries were classified as virulent by the prediction model, VirulentPred. Most of the entries (eighty percent) did not have target compounds based on the database of available chemical compounds at TDRtargets. Fourteen of the twenty entries (seventy percent) had more than 30 consecutively negative amino acid residues based on the ProtScale’s Kyte and Doolittle hydrophobicity plot. Hence, they would be hydrophobic enough to be transmembrane in location or secretory in nature. Only 7 (tyrosinase, serine protease1, Tspan-1, VAL4, cathepsin b and L and calreticulin) had cleavage sites and signal peptides, while none had a significant signal anchor probability. The annotations and characterization provided by this work and the development of a SchistoVir database will aid in further research of schistosome pathogenesis and control.   Key words: Protein database, bioinformatics tools, virulence proteins/factors, annotation, schistosomes.

MALDI mass sequencing and biochemical characterization of Setaria cervi protein tyrosine phosphatase

A 30-kDa acid phosphatase with protein tyrosine phosphatase activity was identified in Setaria cervi (ScPTP). The enzyme was purified to homogeneity using three-step column chromatography. Matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) analysis of purified ScPTP yielded a total of eight peptides matching most closely to phosphoprotein phosphatase of Ricinus communis (RcPP). A hydrophilicity plot of RcPP revealed the presence of these peptides in the hydrophilic region, suggesting their antigenic nature. The substrate specificity of ScPTP with ortho-phospho-L-tyrosine and inhibition with sodium orthovanadate and ammonium molybdate affirmed it as a protein tyrosine phosphatase. ScPTP was also found to be tartrate resistant. The Km and Vmax were 6.60mM and 83.3μM/ml/min, respectively, with pNPP and 8.0mM and 111μM/ml/min, respectively, with ortho-phospho-L-tyrosine as the substrate. The Ki value with sodium orthovanadate was calculated to be 16.10mM. Active site modification with DEPC, EDAC and pHMB suggested the presence of histidine, cysteine and aspartate at its active site. Thus, on the basis of MALDI-TOF and biochemical studies, it was confirmed that purified acid phosphatase is a PTP.

Investigating the interaction between the violaxanthin cycle enzyme zeaxanthin epoxidase and the thylakoid membrane

In the present study the interaction between the violaxanthin cycle enzyme zeaxanthin epoxidase (ZEP) and the thylakoid membrane was investigated. Isolated, active thylakoid membranes of spinach (Spinacia oleracea L.) were subjected to different salt and detergent treatments that are generally used to isolate peripheral and integral membrane proteins. These salt and detergent treatments included the use of the salts NaBr, Na2CO3 and Tris and the detergents octylglucoside (OG) and dodecylmaltoside (DM). After the treatments the activity of the ZEP was determined in washed thylakoid membranes. To obtain additional information about the mode of ZEP binding to the membrane a hydrophobicity plot based on the amino acid sequence of the protein was constructed. The plot was then compared to a diagram obtained for the photosystem II antenna Lhcb1 protein whose integration into the thylakoid membrane is known. The results of the salt and detergent treatments of the thylakoid membrane suggest that the ZEP is a peripheral, rather weakly bound membrane protein. Results from the hydrophobicity plots indicate the existence of specialized protein domains which may realize the partial integration and binding of the ZEP to the thylakoid membrane.

Comparative modeling of DszC, an enzyme in biodesulfurization, and performing in silico point mutation for increasing tendency to oil.

Desulfurization protein named DszC from Rhodococcus erythropolis is the key enzyme for biodesulforization of dibenzothiophene (DBT) in 4S pathway, which is a pathway with four enzymes. DszC enzyme biodesulfurizes DBT and its derivatives in oil components and biphasic systems. It functions well at the oil- water interface. In this study point mutation performed in DszC enzyme regarding to increase protein hydrophobicity and stability for application in immobilized form. 3D model of DszC predicted using Phyre2, SAM-T08 and M4t servers. I-Mutant 2 server used to determine potential spots for point mutation, and Molegro Virtual Docker (MVD) used for performing point mutation on 3D model. Hydrophobicity plots generated by Bioedit version 7.0.8.0 in Kyte-Doolittle scale indicated that protein hydrophobicity is increased after mutation. Also protein stability increased 26.11 units in scale of DDC2.

Bioinformatics analysis of the structure and linear B–cell epitopes of aquaporin–3 from Schistosoma japonicum

Objective To analyze the structure of aquaporins-3(AQP-3) from Schistosoma japonicum(SJAQP-3) using bioinformatical methods, and to provid of references for vaccine targets research. Methods Protparam, BepiPred, TMHMM Server, MLRC, Geno3d, DNA star software packages were used to predict the physical and chemical properties, hydrophilicity plot, flexibility regions, antigenic index, surface probability plot, secondary structure, and tertiary structure of amino acid sequence of SJAQP-3. Results SJAQP-3 had six transmembrane regions and two half-spanning helices that form a central channel. The half-spanning helices fold into the centre of the channel. Either of the half-spanning helix had a conserved motif of NPA common to all aquaporins. Predicted linear B-Cell epitopes were most likely at the N-terminal amino acid residues of 5aa-7aa, 59aa-62aa, 225aa-230aa, 282aa -288aa, 294aa -298aa and 305aa -307aa area. 59aa- 62aa, 225aa-230aa located outside the membrane, the others located inside the cell. Conclusions SJAQP-3 is a integral membrane protein in Schistosoma japonicum tegument. There are six potential epitopes in SJAQP-3. It might be a potential molecular target for the development of vaccines.

Cloning, sequencing, and in silico characterization of Omp 28 of Salmonella Typhi (strain MTCC 733) to develop r-DNA vaccine for typhoid fever.

Background:Typhoid is one of the most important diseases of human beings caused by Salmonella Typhi. There are many vaccine reported against Salmonella Typhi, but search for new candidate vaccine antigens is still going on because presently available vaccines have several limitations such as short-term immunity, high cost, and allergic reaction. Several approaches such as subunit vaccines, Vi polysaccharide, mutant vaccines, and r-DNA vaccines have been tested. r- DNA vaccines have shown some promising potential (targeted Omp). Omp 28 had shown very promising results and suggests that it should be used in further studies of animal protection against the disease.Objective:Cloning, Sequencing and In silico analysis of Omp 28 gene to develop r-DNA vaccine of S. Typhi.Materials and Methods:Omp 28 is made up of three identical subunits of 9.6 kDa showing PCR amplicon of 330 bp which has been cloned in the pJET vector. Recombinant clones has been sequenced, and data submitted to NCBI. Secondary structure was deduced by the Chou Fasman and Garnier method. The sequence of Omp 28 was studied for antigenic indexing, epitope mapping, and MHC mapping using various bioinformatics tool.Results and Conclusion:The sequence of Omp 28 has been assigned accession no GQ 907044.1 by NCBI. Secondary structure has shown it has more alpha region. Hydrophobic plot and surface probability plot shows most amino acids are surface exposed which is a requirement to develop a r-DNA vaccine. Antigenic sites are located within surface exposed regions and eight antigenic determinants are present in Omp 28. On Prosite analysis of Protein shown two motifs i.e. anaphylatoxin domain signature motif at position 219-252 and other one was iron sulphur binding region signature motif at position 36-44. On epitope analysis total six major B cell epitopes were observed which can provoke humoral immunity. On T cell epitope mapping several major epitopes has been found in case of MHC class I and MHC class II. It indicates that Omp 28 can provoke cell mediated as well as humoral immunity and can be proven a promising candidates of Salmonella Typhi.

Docking studies of doxycycline on pathogenic Leptospira Species with common pharmacopore

Every disease can be cured if its roots are traced; means there is cure for every disease. Leptospirosis is a hidden and emerging public health problem as it causes fever, acute renal failure and jaundice. The major outer membrane lipoprotein “LipL32” is expressed during infection which is highly conserved among pathogenic leptospira. It is absent in non-pathogenic avirulent leptospira. Insilico comparison of amino acid sequences of LipL32 showed 97.8% average sequence identity among the pathogenic species. In this study some 3D model were built for LipL32 protein using molecular modelling. Protein hydrophobicity plots were utilized to locate the antigenic sites. Docking studies of retrieved structures and modeled LipL32 proteins with doxycycline revealed that doxycycline can inhibit LipL32 proteins of Leptospira santarosai, Leptospira weilii and Leptospira kirschner as they have bound docking with doxycycline. They share a common pharmacopore within a genus and a common drug may be utilized to combat all the related species.

Phylogenetic analysis of the full genome of rabbit hepatitis E virus (rbHEV) and molecular biologic study on the possibility of cross species transmission of rbHEV

Rabbit HEV strains were isolated in China from farmed rabbits, which may be a reservoir of a novel genotype of HEV. However, there is little information regarding host range and zoonotic potential of rabbit HEV. A previous study from this group identified 25 specific nucleotide positions as possibly being involved in specifying the host range of HEV and/or in determining the severity of hepatitis E disease, whereas this previous analysis did not extend to the new rabbit HEVs. Consequently the complete genome of the rabbit HEV strain ch-bj-n1 was amplified as overlapping fragments using a combination of reverse-transcription-nested polymerase chain reaction and rapid amplification of cDNA ends. The resulting sequence together with two other rabbit HEVs were compared with 91 other HEV isolates representative of genotypes 1–4. Phylogenetic analysis showed that the novel rabbit HEV strain should be classified into a new HEV genotype which has a closer evolutionary relationship with genotype 3 than genotypes 1, 2 and 4, with 20 of the 25 specific nucleotide positions identified as being present in rabbit HEV genome. Alpha-helix, hydrophilicity, and surface probability plots together with antigenic index analysis of the ORF2 protein were all undertaken to investigate the antigenicity of rabbit HEV. The alpha-helix plot and antigenicity index of the rabbit ORF2 protein was nearly identical to that of genotype 1 HEV isolates from human. Furthermore, 20 possible host range determinants in rabbit HEV were identical to those isolates from zoonotic groups, which may suggest a potential possibility of cross species transmission of rabbit HEV.

New User-Friendly Approach to Obtain an Eisenberg Plot and Its Use as a Practical Tool in Protein Sequence Analysis

The Eisenberg plot or hydrophobic moment plot methodology is one of the most frequently used methods of bioinformatics. Bioinformatics is more and more recognized as a helpful tool in Life Sciences in general, and recent developments in approaches recognizing lipid binding regions in proteins are promising in this respect. In this study a bioinformatics approach specialized in identifying lipid binding helical regions in proteins was used to obtain an Eisenberg plot. The validity of the Heliquest generated hydrophobic moment plot was checked and exemplified. This study indicates that the Eisenberg plot methodology can be transferred to another hydrophobicity scale and renders a user-friendly approach which can be utilized in routine checks in protein–lipid interaction and in protein and peptide lipid binding characterization studies. A combined approach seems to be advantageous and results in a powerful tool in the search of helical lipid-binding regions in proteins and peptides. The strength and limitations of the Eisenberg plot approach itself are discussed as well. The presented approach not only leads to a better understanding of the nature of the protein–lipid interactions but also provides a user-friendly tool for the search of lipid-binding regions in proteins and peptides.

Identification of Key Amino Acid Residues That Determine the Ability of High Risk HPV16-E7 to Dysregulate Major Histocompatibility Complex Class I Expression

High risk human Papillomavirus (HPV) types are the major causative agents of cervical cancer. Reduced expression of major histocompatibility complex class I (MHC I) on HPV-infected cells might be responsible for insufficient T cell response and contribute to HPV-associated malignancy. The viral gene product required for subversion of MHC I synthesis is the E7 oncoprotein. Although it has been suggested that high and low risk HPVs diverge in their ability to dysregulate MHC I expression, it is not known what sequence determinants of HPV-E7 are responsible for this important functional difference. To investigate this, we analyzed the capability to affect MHC I of a set of chimeric E7 variants containing sequence elements from either high risk HPV16 or low risk HPV11. HPV16-E7, but not HPV11-E7, causes significant diminution of mRNA synthesis and surface presentation of MHC I, which depend on histone deacetylase activity. Our experiments demonstrate that the C-terminal region within the zinc finger domain of HPV-E7 is responsible for the contrasting effects of HPV11- and HPV16-E7 on MHC I. By using different loss- and gain-of-function mutants of HPV11- and HPV16-E7, we identify for the first time a residue variation at position 88 that is highly critical for HPV16-E7-mediated suppression of MHC I. Furthermore, our studies suggest that residues at position 78, 80, and 88 build a minimal functional unit within HPV16-E7 required for binding and histone deacetylase recruitment to the MHC I promoter. Taken together, our data provide new insights into how high risk HPV16-E7 dysregulates MHC I for immune evasion.

The human polycystin-2 protein represents an integral membrane protein with six membrane-spanning domains and intracellular N- and C-termini

PKD2 is one of the two genes mutated in ADPKD (autosomal-dominant polycystic kidney disease). The protein product of PKD2, polycystin-2, functions as a non-selective cation channel in the endoplasmic reticulum and possibly at the plasma membrane. Hydrophobicity plots and its assignment to the TRP (transient receptor potential) family of cation channels suggest that polycystin-2 contains six transmembrane domains and that both the N- and C-termini extend into the cytoplasm. However, no experimental evidence for this model has so far been provided. To determine the orientation of the different loops of polycystin-2, we truncated polycystin-2 within the predicted loops 1-5 and tagged the constructs at the C-terminus with an HA (haemagglutinin) epitope. After transient expression and selective membrane permeabilization, immunofluorescence staining for the HA epitope revealed that loops 1, 3 and 5 extend into the lumen of the endoplasmic reticulum or the extracellular space, whereas loops 2 and 4 extend into the cytoplasm. This approach also confirmed the cytoplasmic orientation of the N- and C-termini of polycystin-2. In accordance with the immunofluorescence data, protease protection assays from microsomal preparations yielded protected fragments when polycystin-2 was truncated in loops 1, 3 and 5, whereas no protected fragments could be detected when polycystin-2 was truncated in loops 2 and 4. The results of the present study therefore provide the first experimental evidence for the topological orientation of polycystin-2.

Suppression subtractive hybridization analysis of low-protein diet- and vitamin D-induced gene expression from rat kidney inner medullary base

Protein restriction and hypercalcemia result in a urinary concentrating defect in rats and humans. Previous tubular perfusion studies show that there is an increased active urea transport activity in the initial inner medullary (IM) collecting duct in low-protein diet (LPD) and vitamin D (Vit D) animal models. To investigate the possible mechanisms that cause the urinary concentrating defect and to clone the new active urea transporter, we employed a modified two-tester suppression subtractive hybridization (ttSSH) approach and examined gene expression induced by LPD and Vit D in kidney IM base. Approximately 600 clones from the subtracted library were randomly selected; 150 clones were further confirmed to be the true positive genes by slot blot hybridization with subtracted probes from LPD and Vit D and sent for DNA sequencing. We identified 10 channel/transporter genes that were upregulated in IM base in LPD and Vit D animal models; 8 were confirmed by real-time PCR. These genes include aquaporin 2 (AQP2), two-pore calcium channel protein 2, brain-specific organic cation transporter, Na(+)- and H(+)-coupled glutamine transporter, and solute carrier family 25. Nine genes are totally new, and twelve are uncharacterized hypothetical proteins. Among them, four genes were shown to be new transmembrane proteins as judged by Kyte-Doolittle hydrophobic plot analysis. ttSSH provides a useful method to identify new genes from two conditioned populations.

The Molecular Chaperone Hsp90 Modulates Intermediate Steps of Amyloid Assembly of the Parkinson-related Protein α-Synuclein

Alpha-synuclein is an intrinsically unstructured protein that binds to membranes, forms fibrils, and is involved in neurodegeneration. We used a reconstituted in vitro system to show that the molecular chaperone Hsp90 influenced alpha-synuclein vesicle binding and amyloid fibril formation, two processes that are tightly coupled to alpha-synuclein folding. Binding of Hsp90 to monomeric alpha-synuclein occurred in the low micromolar range, involving regions of alpha-synuclein that are critical for vesicle binding and amyloidogenesis. As a consequence, both processes were affected. In the absence of ATP, the accumulation of non-amyloid alpha-synuclein oligomers prevailed over fibril formation, whereas ATP favored fibril growth. This suggests that Hsp90 modulates the assembly of alpha-synuclein in an ATP-dependent manner. We propose that Hsp90 affects these folding processes by restricting conformational fluctuations of alpha-synuclein.

A theoretical analysis of secondary structural characteristics of anticancer peptides

Here, cluster analysis showed that a database of 158 anticancer peptides formed 21 clusters based on net positive charge, hydrophobicity and amphiphilicity. In general, these clusters showed similar median toxicities (P = 0.176) against eukaryotic cell lines and no single combination of these properties was found optimal for efficacy. The database contained 14 peptides, which showed selectivity for tumour cell lines only (ACP(CT)), 123 peptides with general toxicity to eukaryotic cells (ACP(GT)) and 21 inactive peptides (ACP(I)). Hydrophobic arc size analysis showed that there was no significant difference across the datasets although peptides with wide hydrophobic arcs (>270 degrees) appeared to be associated with decreased toxicity. Extended hydrophobic moment plot analysis predicted that over 50% of ACP(CT) and ACP(GT) peptides would be surface active, which led to the suggestion that amphiphilicity is a key driver of the membrane interactions for these peptides but probably plays a role in their efficacy rather than their selectivity. This analysis also predicted that only 14% of ACP(CT) peptides compared to 45% of ACP(GT) peptides were candidates for tilted peptide formation, which led to the suggestion that the absence of this structure may support cancer cell selectivity. However, these analyses predicted that ACP(I) peptides, which possess no anticancer activity, would also form surface active and tilted alpha-helices, clearly showing that other factors are involved in determining the efficacy and selectivity of ACPs.

GLUT8, the enigmatic intracellular hexose transporter

GLUT8 is a class III sugar transporter predominantly expressed in testis and brain. In contrast to the class I and class II transporters, hydrophobicity plots predict a short extracellular loop between transmembrane domain (TM)1 and TM2 and a long extracellular loop between TM9 and TM10 that contains the only N-glycosylation site. In vitro translated GLUT8 migrates as a 35-kDa protein that is glycosylated in the presence of microsomal membranes. In heterologous expression systems, glucose transport activity (Km of 2 mM) was inhibited by fructose and galactose. The transporter carries an NH2-terminal endosomal/lysosomal targeting motif ([DE]XXXL[LI]). Accordingly, constitutive GLUT8 has been found to be associated with endosomes and lysosomes but also with membranes of the endoplasmic reticulum. A similar distribution was detected after overexpression of wild-type or tagged GLUT8 in different cell systems. In these cells, none of the conventional signals tested induced a translocation of GLUT8 to the plasma membrane. Therefore, GLUT8 appears to catalyze transport of sugars or sugar derivatives through intracellular membranes. Slc2a8 knockout mice were viable, developed normally, and showed mild alterations in brain (increased proliferation of neuronal cells in dentate gyrus of the hippocampus, hyperactivity), heart (impaired transmission of electrical wave through the atrium), and sperm cells (reduced number of motile sperm cells associated with reduced mitochondrial membrane potential and ATP levels in sperm). The links between molecular function, cellular localization and phenotype of the knockout mouse is unclear and remains to be determined.

Stepwise prediction and statistical screening: B-cell epitopes on neuraminidase of human avian H5N1 virus.

The B-cell epitopes of virus are associated with the antiviral drug and the vaccine screening. As the nucleotide sequences of neuraminidase (NA) of stain GD-01-06 were sequenced, we predicted the α-helix and β-fold structure and the indexes of the flexible regions of secondary structure of NA with methods of the Hydrophilicity plot by Kyte-Doolittle, the Surface probability plot by Emini and the Antigenic index by Jameson-Wolf, and then screened statistically the parameters to predict B-cell epitopes by the Hierarchical cluster and the Bivariate correlation and the quartiles with SPSS 13.0. The impact of variation of amino acids in NA on its epitopes was analyzed. The predictive results were evaluated by Wu’s Antigenic Index and SWISS-MODEL. We found that the most possible epitopes on NA were located within or nearby its N-terminal Nos. 120–137, 81–84, 408–415, 273–282, 429–432, 356–368, 46–55, 146–155, 341–350 and 198–209, which were the dominant regions of NA epitopes. Peptide 120–137 including the glycoprotein domain (NGT126–128) was first chosen as the B-cell epitopes on NA. NA in H5N1 strain isolated after 2003 lacked in No. 53 amino acid (I), resulting in an increase in the surface flexible region of NA in GD-01-06 and an enlargement to their epitope regions (VEP46–48 → VEPISNTNFL46–55). Conclusively, prediction of the B-cell epitopes on the NA based on multiple parameters is useful for researches on the molecular immunology and drug screening and immuno-prophylaxis. A deletion of No. 53 amino acid (I) in NA in strain GD-01-06 might increase its antigenicity.

Expression and Characterization of theStreptomyces coelicolorSerine/Threonine Protein Kinase PkaD

We identified and characterized the gene encoding a new eukaryotic-type protein kinase from Streptomyces coelicolor A3(2) M145. PkaD, consisting of 598 amino acid residues, contained the catalytic domain of eukaryotic protein kinases in the N-terminal region. A hydrophobicity plot indicated the presence of a putative transmembrane spanning sequence downstream of the catalytic domain, suggesting that PkaD is a transmembrane protein kinase. The recombinant PkaD was found to be phosphorylated at the threonine and tyrosine residues. In S. coelicolor A3(2), pkaD was transcribed as a monocistronic mRNA, and it was expressed constitutively throughout the life cycle. Disruption of chromosomal pkaD resulted in a significant loss of actinorhodin production. This result implies the involvement of pkaD in the regulation of secondary metabolism.

Molecular analysis of physiological responses to changes in nitrogen in a marine macroalga, Porphyra yezoensis (Rhodophyta)

The rhodophyte seaweed Porphyra yezoensis, known more commonly world-wide as "nori", is an important commercial crop in Japan. Cultivation of nori in Japan is often affected by outbreaks of "iroochi", a discoloration of the thalli due to a decrease in inorganic nutrients in the culture area that in turn decreases the amount of photosynthetic pigments in the thalli. Treating thalli with inorganic nitrogen can reverse iroochi. In this paper, we report on the characterization of three P. yezoensis genes, a nitrate transporter (PyNRT2) and two urea transporters (PyUT1 and PyUT2), which may be involved in reversing iroochi. The predicted length of the PyNRT2 protein was 479 amino acids (AA), while PyUT1 and PyUT2 were 740 and 680 AA, respectively. PyNRT2 was more similar to NRT2 from a chromophyte than to NRTs from Chlamydomonas and higher plants. The two P. yezoensis UTs had 56% AA identity to each other, and showed the closest relationship to higher plant and yeast DUR3 proteins which formed a subfamily of the sodium-solute symporter protein family. Hydrophobicity plots of the AA sequences showed that the PyNRT2, PyUT1, and PyUT2 included 12, 15, and 16 transmembrane domains, respectively. Southern blot analysis indicated that the P. yezoensis genome has a single NRT2-encoding gene and at least four UT-encoding genes. Expression analysis of PyNRT2 and PyUT genes showed that the messenger RNA level of the PyNRT2 gene reached a maximum after 48 h in the nitrate starvation condition and was then restored to the constitutive level, while expression of the PyUT genes was induced in proportion to treatment times in the nitrate starvation condition. These results suggest that the PyNRT2 and PyUT are responsible for the high-affinity nitrate/urea transport systems that operate under low external nitrate concentrations.

Analysis and prediction of protective continuous B-cell epitopes on pathogen proteins

BackgroundThe application of peptide based diagnostics and therapeutics mimicking part of protein antigen is experiencing renewed interest. So far selection and design rationale for such peptides is usually driven by T-cell epitope prediction, available experimental and modelled 3D structure, B-cell epitope predictions such as hydrophilicity plots or experience. If no structure is available the rational selection of peptides for the production of functionally altering or neutralizing antibodies is practically impossible. Specifically if many alternative antigens are available the reduction of required synthesized peptides until one successful candidate is found is of central technical interest. We have investigated the integration of B-cell epitope prediction with the variability of antigen and the conservation of patterns for post-translational modification (PTM) prediction to improve over state of the art in the field. In particular the application of machine-learning methods shows promising results.ResultsWe find that protein regions leading to the production of functionally altering antibodies are often characterized by a distinct increase in the cumulative sum of three presented parameters. Furthermore the concept to maximize antigenicity, minimize variability and minimize the likelihood of post-translational modification for the identification of relevant sites leads to biologically interesting observations. Primarily, for about 50% of antigen the approach works well with individual area under the ROC curve (AROC) values of at least 0.65. On the other hand a significant portion reveals equivalently low AROC values of < = 0.35 indicating an overall non-Gaussian distribution. While about a third of 57 antigens are seemingly intangible by our approach our results suggest the existence of at least two distinct classes of bioinformatically detectable epitopes which should be predicted separately. As a side effect of our study we present a hand curated dataset for the validation of protectivity classification. Based on this dataset machine-learning methods further improve predictive power to a class separation in an equilibrated dataset of up to 83%.ConclusionWe present a computational method to automatically select and rank peptides for the stimulation of potentially protective or otherwise functionally altering antibodies. It can be shown that integration of variability, post-translational modification pattern conservation and B-cell antigenicity improve rational selection over random guessing. Probably more important, we find that for about 50% of antigen the approach works substantially better than for the overall dataset of 57 proteins. Essentially as a side effect our method optimizes for presumably best applicable peptides as they tend to be likely unmodified and as invariable as possible which is answering needs in diagnosis and treatment of pathogen infection. In addition we show the potential for further improvement by the application of machine-learning methods, in particular Random Forests.

Relative Hydrophobicity/Hydrophilicity of Fructose, Glucose, Sucrose, and Trehalose as Probed by 1-Propanol: A Differential Approach in Solution Thermodynamics

Earlier, we developed the 1-propanol probing methodology that could separately evaluate the relative hydrophobicity and hydrophilicity scale of a given sample. We applied it here to fructose, glucose, sucrose, and trehalose and evaluated the same for these saccharides. We then construct a two-dimensional map with the hydrophobicity and the hydrophilicity axes and plot the above saccharides together with other nonelectrolytes subjected earlier to the same analysis. We point out that these saccharides together with other so-called "osmolytes" that accumulate in vivo under H2O stress occupy a small specific region near H2O in this map.