Polypeptide amino acid composition and isoelectric point. II. Comparison between experiment and theory

Purification and biofabrication of 5-aminolevulinic acid for photodynamic therapy against pathogens and cancer cells

5-Aminolevulinic acid (5-ALA) is a non-proteinogenic amino acid which has involved in heme metabolism of organisms, and has been widely applied in agriculture, and medical fields nowadays. 5-ALA is used in the elimination of pathogens or cancer cells by photodynamic therapy (PDT) owing to the photosensitizer reaction which releases the reactive oxygen species (ROS). Currently, biofabrication of 5-ALA is regarded as the most efficient and eco-friendly approach, but the complicated ingredient of medium causes the nuisance process of purification, resulting in low recovery and high producing cost. In this study, hydrogen chloride, sodium acetate, and ammonia were examined to maximize the recovery of 5-ALA from ion-exchange chromatography (IEC), thus a 92% recovery in 1 M ammonia at pH 9.5 was obtained. Afterward, the activated carbon was used for decolorization to further remove the pigments from the eluent. Four organic solvents, i.e., diethyl ether, methanol, ethanol, and acetone were compared to extract and form 5-ALA precipitation. The purified 5-ALA was verified to eliminate 74% of A549 human lung cancer and 83% of A375 melanoma skin cancer cell. Moreover, Proteus hauseri, Aeromonas hydrophila, Bacillus cereus, and Staphylococcus aureus were killed via anti-microbial PDT with 1% 5-ALA and reached 100% killing rate at optimal condition. With the addition of 0.05% 5-ALA during the culture, the growth of microalgae Chlorella sorokiniana was improved to against a common aquatic pathogen, A. hydrophila. The broad application of 5-ALA was demonstrated in this study for the first time.

Proteome-pI 2.0: proteome isoelectric point database update

Proteome-pI 2.0 is an update of an online database containing predicted isoelectric points and pKa dissociation constants of proteins and peptides. The isoelectric point-the pH at which a particular molecule carries no net electrical charge-is an important parameter for many analytical biochemistry and proteomics techniques. Additionally, it can be obtained directly from the pKa values of individual charged residues of the protein. The Proteome-pI 2.0 database includes data for over 61 million protein sequences from 20 115 proteomes (three to four times more than the previous release). The isoelectric point for proteins is predicted by 21 methods, whereas pKa values are inferred by one method. To facilitate bottom-up proteomics analysis, individual proteomes were digested in silico with the five most commonly used proteases (trypsin, chymotrypsin, trypsin + LysC, LysN, ArgC), and the peptides' isoelectric point and molecular weights were calculated. The database enables the retrieval of virtual 2D-PAGE plots and customized fractions of a proteome based on the isoelectric point and molecular weight. In addition, isoelectric points for proteins in NCBI non-redundant (nr), UniProt, SwissProt, and Protein Data Bank are available in both CSV and FASTA formats. The database can be accessed at http://isoelectricpointdb2.org.

IPC 2.0: prediction of isoelectric point and pKa dissociation constants

The isoelectric point is the pH at which a particular molecule is electrically neutral due to the equilibrium of positive and negative charges. In proteins and peptides, this depends on the dissociation constant (pK_a) of charged groups of seven amino acids and NH⁺ and COO⁻ groups at polypeptide termini. Information regarding isoelectric point and pK_a is extensively used in two-dimensional gel electrophoresis (2D-PAGE), capillary isoelectric focusing (cIEF), crystallisation, and mass spectrometry. Therefore, there is a strong need for the in silico prediction of isoelectric point and pK_a values. In this paper, I present Isoelectric Point Calculator 2.0 (IPC 2.0), a web server for the prediction of isoelectric points and pK_a values using a mixture of deep learning and support vector regression models. The prediction accuracy (RMSD) of IPC 2.0 for proteins and peptides outperforms previous algorithms: 0.848 versus 0.868 and 0.222 versus 0.405, respectively. Moreover, the IPC 2.0 prediction of pK_a using sequence information alone was better than the prediction from structure-based methods (0.576 versus 0.826) and a few folds faster. The IPC 2.0 webserver is freely available at www.ipc2-isoelectric-point.org

Role of pKA in Charge Regulation and Conformation of Various Peptide Sequences

Peptides containing amino acids with ionisable side chains represent a typical example of weak ampholytes, that is, molecules with multiple titratable acid and base groups, which generally exhibit charge regulating properties upon changes in pH. Charged groups on an ampholyte interact electrostatically with each other, and their interaction is coupled to conformation of the (macro)molecule, resulting in a complex feedback loop. Their charge-regulating properties are primarily determined by the pKA of individual ionisable side-chains, modulated by electrostatic interactions between the charged groups. The latter is determined by the amino acid sequence in the peptide chain. In our previous work we introduced a simple coarse-grained model of a flexible peptide. We validated it against experiments, demonstrating its ability to quantitatively predict charge on various peptides in a broad range of pH. In the current work, we investigated two types of peptide sequences: diblock and alternating, each of them consisting of an equal number of amino acids with acid and base side-chains. We showed that changing the sequence while keeping the same overall composition has a profound effect on the conformation, whereas it practically does not affect total charge on the peptide. Nevertheless, the sequence significantly affects the charge state of individual groups, showing that the zero net effect on the total charge is a consequence of unexpected cancellation of effects. Furthermore, we investigated how the difference between the pKA of acid and base side chains affects the charge and conformation of the peptide, showing that it is possible to tune the charge-regulating properties by following simple guiding principles based on the pKA and on the amino acid sequence. Our current results provide a theoretical basis for understanding of the complex coupling between the ionisation and conformation in flexible polyampholytes, including synthetic polymers, biomimetic materials and biological molecules, such as intrinsically disordered proteins, whose function can be regulated by changes in the pH.

Accurate Estimation of pKb Values for Amino Groups from Surface Electrostatic Potential (VS,min) Calculations: The Isoelectric Points of Amino Acids as a Case Study

Theoretical calculation of equilibrium dissociation constants is a very computationally demanding and time-consuming process since it requires an extremely accurate computation of the solvation free energy changes for each of the species involved. By correlating the minimum surface electrostatic potential (V_S,min) on the nitrogen atom of several aliphatic amino groups-calculated at the density functional theory (DFT) ωB97X-D/cc-pVDZ level of theory-we obtained regression models for each kind of substitution pattern from which we interpolate their corresponding pK_b values with remarkable accuracy: primary R² = 0.9519; secondary R² = 0.9112; and tertiary R² = 0.8172 (N = 20 for each family). These models were validated with tests sets (N = 5) with mean absolute error (MAE) values of 0.1213 (primary), 0.4407 (secondary), and 0.3057 (tertiary). Combining this ansatz with another previously reported by our group to estimate pK_a values [Caballero-García, G.; et al. Molecules 2019, 24(1), 79] we are able to reproduce the isoelectric points of 13 amino acids with no titrable side chains with MAE = 0.4636 pI units.

Identification of duck liver-expressed antimicrobial peptide 2 and characterization of its bactericidal activity

Objective

This study was conducted to identify duck liver-expressed antimicrobial peptide 2 (LEAP-2) and demonstrate its antimicrobial activity against various pathogens.

Methods

Tissue samples were collected from 6 to 8-week-old Pekin ducks (Anas platyrhynchos domesticus), total RNA was extracted, and cDNA was synthesized. To confirm the duck LEAP-2 transcript expression levels, quantitative real-time polymerase chain reaction was conducted. Two kinds of peptides (a linear peptide and a disulfide-type peptide) were synthesized to compare the antimicrobial activity. Then, antimicrobial activity assay and fluorescence microscopic analysis were conducted to demonstrate duck LEAP-2 bactericidal activity.

Results

The duck LEAP-2 peptide sequence showed high identity with those of other avian species (>85%), as well as more than 55% of identity with mammalian sequences. LEAP-2 mRNA was highly expressed in the liver with duodenum next, and then followed by lung, spleen, bursa and jejunum and was the lowest in the muscle. Both of LEAP-2 peptides efficiently killed bacteria, although the disulfide-type LEAP-2 showed more powerful bactericidal activity. Also, gram-positive bacteria was more susceptible to duck LEAP-2 than gram-negative bacteria. Using microscopy, we confirmed that LEAP-2 peptides could kill bacteria by disrupting the bacterial cell envelope.

Conclusion

Duck LEAP-2 showed its antimicrobial activity against both gram-positive and gram-negative bacteria. Disulfide bonds were important for the powerful killing effect by disrupting the bacterial cell envelope. Therefore, duck LEAP-2 can be used for effective antibiotics alternatives.

Protein surface functionalisation as a general strategy for facilitating biomimetic mineralisation of ZIF-8

The surface charge and chemistry of a protein determines its ability to facilitate biomimetic mineralisation.

The durability of enzymes in harsh conditions can be enhanced by encapsulation within metal–organic frameworks (MOFs) via a process called biomimetic mineralisation. Herein we show that the surface charge and chemistry of a protein determines its ability to seed MOF growth. We demonstrate that chemical modification of amino acids on the protein surface is an effective method for systematically controlling biomimetic mineralisation by zeolitic imidazolate framework-8 (ZIF-8). Reaction of surface lysine residues with succinic (or acetic) anhydride facilitates biomimetic mineralisation by increasing the surface negative charge, whereas reaction of surface carboxylate moieties with ethylenediamine affords a more positively charged protein and hinders the process. Moreover, computational studies confirm that the surface electrostatic potential of a protein is a good indicator of its ability to induce biomimetic mineralisation. This study highlights the important role played by protein surface chemistry in encapsulation and outlines a general method for facilitating the biomimetic mineralisation of proteins.

Design of a software for calculating isoelectric point of a polypeptide according to their net charge using the graphical programming language LabVIEW

A software to calculate the net charge and to predict the isoelectric point (pI) of a polypeptide is developed in this work using the graphical programming language LabVIEW. Through this instrument the net charges of the ionizable residues of the polypeptide chains of the proteins are calculated at different pH values, tabulated, pI is predicted and an Excel (-xls) type file is generated. In this work, the experimental values of the pIs (pI) of different proteins are compared with the values of the pIs (pI) calculated graphically, achieving a correlation coefficient (R) of 0.934746 which represents a good reliability for a p < 0.01. In this way the generated program can constitute an instrument applicable in the laboratory, facilitating the calculation to graduate students and junior researchers. © 2017 by The International Union of Biochemistry and Molecular Biology, 46(1):39-46, 2018.

A combined experimental and molecular simulation study of factors influencing interaction of quinoa proteins-carrageenan

The interaction between quinoa proteins isolate (QP isolate) and the negatively charged polysaccharide ι-Carragennan (Carr) as a function of pH was studied. Experimental measurements as turbidity, hydrophobic surface, ζ-potential, and hydrodynamic size were carried out. Associative interaction between QP and Carr was found in the pH range between 1 and 2.9. When both molecules are negatively charged (pH>5,5), a pure Coulombic repulsion regime is observed and the self-association of QP due to the Carr exclusion is proposed. In the intermediate pH range, the experimental data suggests that the charge regulation mechanism can overcome the electrostatic repulsion that may take place (and an attraction between QP and Carr can still be observed). Computational simulations by means of free energy derivatives using the Monte Carlo method were carried out to better understand the interaction mechanism between QP and Carr. QP was modeled as a single protein using one of the major proteins, Chenopodin (Ch), and Carr was modeled as a negatively charged polyelectrolyte (NCP) chain, both in the cell model framework. Simulation results showed attractive interactions in agreement with the experimental data.

Protein charge distribution in proteomes and its impact on translation

As proteins are synthesized, the nascent polypeptide must pass through a negatively charged exit tunnel. During this stage, positively charged stretches can interact with the ribosome walls and slow the translation. Therefore, charged polypeptides may be important factors that affect protein expression. To determine the frequency and distribution of positively and negatively charged stretches in different proteomes, the net charge was calculated for every 30 consecutive amino acid residues, which corresponds to the length of the ribosome exit tunnel. The following annotated and reviewed proteins in the UniProt database (Swiss-Prot) were analyzed: 551,705 proteins from different organisms and a total of 180 million protein segments. We observed that there were more negative than positive stretches and that super-charged positive sequences (i.e., net charges ≥ 14) were underrepresented in the proteomes. Overall, the proteins were more positively charged at their N-termini and C-termini, and this feature was present in most organisms and subcellular localizations. To investigate whether the N-terminal charges affect the elongation rates, previously published ribosomal profiling data obtained from S. cerevisiae, without translation-interfering drugs, were analyzed. We observed a nonlinear effect of the charge on the ribosome occupancy in which values ≥ +5 and ≤ -6 showed increased and reduced ribosome densities, respectively. These groups also showed different distributions across 80S monosomes and polysomes. Basic polypeptides are more common within short proteins that are translated by monosomes, whereas negative stretches are more abundant in polysome-translated proteins. These findings suggest that the nascent peptide charge impacts translation and can be one of the factors that regulate translation efficiency and protein expression.

Which factors shape the sequence of amino acids that will form a protein? The biochemical features of amino acids, such as their charge and hydrophobicity, are important drivers of protein tridimensional folding, which creates interaction sites for binding other molecules and directs proteins to specific cellular compartments. These features all impact the activity of the proteins after they are produced. Another less obvious factor that influences the protein’s primary structure may be how efficiently a given amino acid sequence is produced by the ribosome. It is known that a repetitive stretch of positively charged amino acids may interact with the negative charges in the ribosome exit tunnel, slowing, or even halting, translation. By analyzing the charge of protein stretches in different organisms, we observed that proteins tend to present positively charged stretches at their extremities, and high charge values can slow (for positive charges) or speed (for negative charges) translation. An interesting consequence of this trend is that proteins that are translated in high quantities by several ribosomes at the same RNA (polysomes) tend to have more negatively charged stretches than proteins that are translated by a single ribosome per RNA (monosomes).

Proteome-pI: proteome isoelectric point database

Proteome-pI is an online database containing information about predicted isoelectric points for 5029 proteomes calculated using 18 methods. The isoelectric point, the pH at which a particular molecule carries no net electrical charge, is an important parameter for many analytical biochemistry and proteomics techniques, especially for 2D gel electrophoresis (2D-PAGE), capillary isoelectric focusing, liquid chromatography–mass spectrometry and X-ray protein crystallography. The database, available at http://isoelectricpointdb.org allows the retrieval of virtual 2D-PAGE plots and the development of customised fractions of proteome based on isoelectric point and molecular weight. Moreover, Proteome-pI facilitates statistical comparisons of the various prediction methods as well as biological investigation of protein isoelectric point space in all kingdoms of life. For instance, using Proteome-pI data, it is clear that Eukaryotes, which evolved tight control of homeostasis, encode proteins with pI values near the cell pH. In contrast, Archaea living frequently in extreme environments can possess proteins with a wide range of isoelectric points. The database includes various statistics and tools for interactive browsing, searching and sorting. Apart from data for individual proteomes, datasets corresponding to major protein databases such as UniProtKB/TrEMBL and the NCBI non-redundant (nr) database have also been precalculated and made available in CSV format.

IPC - Isoelectric Point Calculator

BACKGROUND

Accurate estimation of the isoelectric point (pI) based on the amino acid sequence is useful for many analytical biochemistry and proteomics techniques such as 2-D polyacrylamide gel electrophoresis, or capillary isoelectric focusing used in combination with high-throughput mass spectrometry. Additionally, pI estimation can be helpful during protein crystallization trials.

RESULTS

Here, I present the Isoelectric Point Calculator (IPC), a web service and a standalone program for the accurate estimation of protein and peptide pI using different sets of dissociation constant (pKa) values, including two new computationally optimized pKa sets. According to the presented benchmarks, the newly developed IPC pKa sets outperform previous algorithms by at least 14.9 % for proteins and 0.9 % for peptides (on average, 22.1 % and 59.6 %, respectively), which corresponds to an average error of the pI estimation equal to 0.87 and 0.25 pH units for proteins and peptides, respectively. Moreover, the prediction of pI using the IPC pKa's leads to fewer outliers, i.e., predictions affected by errors greater than a given threshold.

CONCLUSIONS

The IPC service is freely available at http://isoelectric.ovh.org Peptide and protein datasets used in the study and the precalculated pI for the PDB and some of the most frequently used proteomes are available for large-scale analysis and future development.

REVIEWERS

This article was reviewed by Frank Eisenhaber and Zoltán Gáspári.

Photoreponsive Hybrid Nanoparticles with Inherent FRET Activity

The photoactivated inherent fluorescence resonance energy transfer (FRET) properties of a hard-and-soft hybrid nanosystem comprising poly(1'-(2-methacryloxyethyl)-3',3'-dimethyl-6-nitrospiro-(2H-1-benzopyran-2,2'-indoline))-co-poly[2-(dimethylamino)ethyl methacrylate] (PSPMA-co-PDMAEMA) random copolymer brushes on silica nanoparticles are described. This unique FRET process is switched on by the simultaneous generation of isomer X and merocyanine (MC), which are bipolar in nature and comprise donor-acceptor dyads, from a single spiropyran (SP) chromophore upon UV irradiation. These X-MC species exhibit sufficient lifetimes to allow the read-out of the FRET process. The phenomenon is gradually switched off because of the thermal relaxation of the bipolar chromophores. This inherent property of the nanoemitters is employed in the development of biosensors of high specificity by monitoring variations in the FRET efficiency and lifetime of the hybrids in the presence of biological substances. More specifically, bovine serum albumin (BSA) augments the formation of MC species and retards the MC photobleaching process, leading to the enhancement of the FRET efficiency and lifetime, respectively. On the other hand, amino acid l-histidine further retards the MC thermal relaxation and prolongs the FRET process. We envisage that this platform opens new perspectives in the development of novel, optical nanosensors for applications in various fields including healthcare products and environmental monitoring.

Accurate estimation of isoelectric point of protein and peptide based on amino acid sequences

Motivation: In any macromolecular polyprotic system—for example protein, DNA or RNA—the isoelectric point—commonly referred to as the pI—can be defined as the point of singularity in a titration curve, corresponding to the solution pH value at which the net overall surface charge—and thus the electrophoretic mobility—of the ampholyte sums to zero. Different modern analytical biochemistry and proteomics methods depend on the isoelectric point as a principal feature for protein and peptide characterization. Protein separation by isoelectric point is a critical part of 2-D gel electrophoresis, a key precursor of proteomics, where discrete spots can be digested in-gel, and proteins subsequently identified by analytical mass spectrometry. Peptide fractionation according to their pI is also widely used in current proteomics sample preparation procedures previous to the LC-MS/MS analysis. Therefore accurate theoretical prediction of pI would expedite such analysis. While such pI calculation is widely used, it remains largely untested, motivating our efforts to benchmark pI prediction methods.

Results: Using data from the database PIP-DB and one publically available dataset as our reference gold standard, we have undertaken the benchmarking of pI calculation methods. We find that methods vary in their accuracy and are highly sensitive to the choice of basis set. The machine-learning algorithms, especially the SVM-based algorithm, showed a superior performance when studying peptide mixtures. In general, learning-based pI prediction methods (such as Cofactor, SVM and Branca) require a large training dataset and their resulting performance will strongly depend of the quality of that data. In contrast with Iterative methods, machine-learning algorithms have the advantage of being able to add new features to improve the accuracy of prediction.

Contact: yperez@ebi.ac.uk

Availability and Implementation: The software and data are freely available at https://github.com/ypriverol/pIR.

Supplementary information:Supplementary data are available at Bioinformatics online.

On-line capillary isoelectric focusing hyphenated to native electrospray ionization mass spectrometry for the characterization of interferon-γ and variants

The on-line hyphenation of Capillary IsoElectric Focusing (CIEF) with ElectroSpray Ionization Mass Spectrometry (ESI/MS) has been carried out in a non-denaturing detection mode at the CIEF-MS interface.

Use of SciDBMaker as Tool for the Design of Specialized Biological Databases

The exponential growth of molecular biology research in recent decades has brought concomitant growth in the number and size of genomic and proteomic databases used to interpret experimental findings. Particularly, growth of protein sequence records created the need for smaller and manually annotated databases. Since scientists are continually developing new specific databases to enhance their understanding of biological processes, the authors created SciDBMaker to provide a tool for easy building of new specialized protein knowledge bases. This chapter also suggests best practices for specialized biological databases design, and provides examples for the implementation of these practices.

Calibrative approaches to protein solubility modeling of a mutant series using physicochemical descriptors

A set of physicochemical properties describing a protein of known structure is employed for a calibrative approach to protein solubility. Common hydrodynamic and electrophoretic properties routinely measured in the bio-analytical laboratory such as zeta potential, dipole moment, the second osmotic virial coefficient are first estimated in silico as a function a pH and solution ionic strength starting with the protein crystal structure. The utility of these descriptors in understanding the solubility of a series of ribonuclease Sa mutants is investigated. A simple two parameter model was trained using solubility data of the wild type protein measured at a restricted number of solution pHs. Solubility estimates of the mutants demonstrate that zeta potential and dipole moment may be used to rationalize solubility trends over a wide pH range. Additionally a calibrative model based on the protein's second osmotic virial coefficient, B₂₂ was developed. A modified DVLO type potential along with a simplified representation of the protein allowed for efficient computation of the second viral coefficient. The standard error of prediction for both models was on the order of 0.3 log S units. These results are very encouraging and demonstrate that these models may be trained with a small number of samples and employed extrapolatively for estimating mutant solubilities.

The isoelectric region of proteins: a systematic analysis

Background

Binding of proteins in ion exchange chromatography is dominated by electrostatic interactions and can be tuned by adjusting pH and ionic strength of the solvent. Therefore, the isoelectric region (IER), the pH region of almost zero charge near the pI, has been used to predict the binding properties of proteins.

Principal findings

Usually the IER is small and binding and elution is carried out at pH values near to the pI. However, some proteins with an extended IER have been shown to bind and elute far away from its pI. To analyze factors that mediate the size of the IER and to identify proteins with an extended IER, two protein families consisting of more than 7000 proteins were systematically investigated. Most proteins were found to have a small IER and thus are expected to bind or elute near to their pI, while only a small fraction of less than 2% had a large IER.

Conclusions

Only four factors, the number of histidines, the pI, the number of titratable amino acids and the ratio of acidic to basic residues, are sufficient to reliably classify proteins by their IER based on their sequence only, and thus to predict their binding and elution behaviour in ion exchange chromatography.

Genomic and bioinformatics analysis of human adenovirus type 37: new insights into corneal tropism

Background

Human adenovirus type 37 (HAdV-37) is a major etiologic agent of epidemic keratoconjunctivitis, a common and severe eye infection associated with long-term visual morbidity due to persistent corneal inflammation. While HAdV-37 has been known for over 20 years as an important cause, the complete genome sequence of this serotype has yet to be reported. A detailed bioinformatics analysis of the genome sequence of HAdV-37 is extremely important to understanding its unique pathogenicity in the eye.

Results

We sequenced and annotated the complete genome of HAdV-37, and performed genomic and bioinformatics comparisons with other HAdVs to identify differences that might underlie the unique corneal tropism of HAdV-37. Global pairwise genome alignment with HAdV-9, a human species D adenovirus not associated with corneal infection, revealed areas of non-conserved sequence principally in genes for the virus fiber (site of host cell binding), penton (host cell internalization signal), hexon (principal viral capsid structural protein), and E3 (site of several genes that mediate evasion of the host immune system). Phylogenetic analysis revealed close similarities between predicted proteins from HAdV-37 of species D and HAdVs from species B and E. However, virtual 2D gel analyses of predicted viral proteins uncovered unexpected differences in pI and/or size of specific proteins thought to be highly similar by phylogenetics.

Conclusion

This genomic and bioinformatics analysis of the HAdV-37 genome provides a valuable tool for understanding the corneal tropism of this clinically important virus. Although disparities between HAdV-37 and other HAdV within species D in genes encoding structural and host receptor-binding proteins were to some extent expected, differences in the E3 region suggest as yet unknown roles for this area of the genome. The whole genome comparisons and virtual 2D gel analyses reported herein suggest potent areas for future studies.

SciDBMaker: new software for computer-aided design of specialized biological databases

Background

The exponential growth of research in molecular biology has brought concomitant proliferation of databases for stocking its findings. A variety of protein sequence databases exist. While all of these strive for completeness, the range of user interests is often beyond their scope. Large databases covering a broad range of domains tend to offer less detailed information than smaller, more specialized resources, often creating a need to combine data from many sources in order to obtain a complete picture. Scientific researchers are continually developing new specific databases to enhance their understanding of biological processes.

Description

In this article, we present the implementation of a new tool for protein data analysis. With its easy-to-use user interface, this software provides the opportunity to build more specialized protein databases from a universal protein sequence database such as Swiss-Prot. A family of proteins known as bacteriocins is analyzed as 'proof of concept'.

Conclusion

SciDBMaker is stand-alone software that allows the extraction of protein data from the Swiss-Prot database, sequence analysis comprising physicochemical profile calculations, homologous sequences search, multiple sequence alignments and the building of new and more specialized databases. It compiles information with relative ease, updates and compares various data relevant to a given protein family and could solve the problem of dispersed biological search results.

Improved application of the oscillating method for the isoelectric point determination of protein: Potential connection with protein data banks

The oscillating method (OM) for the theoretical determination of the pI values, one by one, of proteins and other macromolecules has been previously published [Sillero and Maldonado, Comput. Biol. Med 36 (2006) 157-166]. An improved application of the method, here named as improved oscillating method (IOM), allows the pI determination of group of proteins. This characteristic may be useful to explore the pI value and electric charge of family of enzymes. As an example the pI values of 1630 enzymes collected in a Swiss-Prot data bank (www.expasy.org), as belonging to the enzymes ligases (EC 6. 2. 1. *) is presented. The method also permits the determination of the pI value of any group of proteins stored in data banks provided that they can be supplied to the program in a FASTA format. Free access to IOM can be reached at http://www.bq.uam.es/otros/pical3.zip.

Synthesis, characterization, and DNA adsorption studies of ampholytic model conetworks based on cross-linked star copolymers

Five model conetworks based on cross-linked star ampholytic copolymers were synthesized by group transfer polymerization. The ampholytic copolymers were based on two hydrophilic monomers: the positively ionizable 2-(dimethylamino)ethyl methacrylate (DMAEMA) and the negatively ionizable methacrylic acid (MAA). Ethylene glycol dimethacrylate was used as the cross-linker. These five ampholytic model conetworks were isomers based on equimolar DMAEMA-MAA copolymer stars of different architectures: heteroarm (two), star block (two), and statistical. The two networks based on the homopolymer stars were also synthesized. The MAA units were introduced via the polymerization of tetrahydropyranyl methacrylate and the acid hydrolysis of the latter after network formation. All the precursors to the (co)networks were characterized in terms of their molecular weights using gel permeation chromatography (GPC). The mass of the extractables from the (co)networks was measured and characterized in terms of molecular weight and composition using GPC and proton nuclear magnetic resonance (1H NMR) spectroscopy, respectively. The degrees of swelling (DS) of all the ampholytic conetworks were measured as a function of pH and were found to present a minimum at a pH value which was taken as the isoelectric point, pI. The DS and the pI values did not present a dependence on conetwork architecture. Finally, DNA adsorption studies onto the ampholyte conetworks indicated that DNA binding was governed by electrostatics.

Patterned Biofunctional Poly(acrylic acid) Brushes on Silicon Surfaces

Protein patterning was carried out using a simple procedure based on photolithography wherein the protein was not subjected to UV irradiation and high temperatures or contacted with denaturing solvents or strongly acidic or basic solutions. Self-assembled monolayers of poly(ethylene glycol) (PEG) on silicon surfaces were exposed to oxygen plasma through a patterned photoresist. The etched regions were back-filled with an initiator for surface-initiated atom transfer radical polymerization (ATRP). ATRP of sodium acrylate was readily achieved at room temperature in an aqueous medium. Protonation of the polymer resulted in patterned poly(acrylic acid) (PAA) brushes. A variety of biomolecules containing amino groups could be covalently tethered to the dense carboxyl groups of the brush, under relatively mild conditions. The PEG regions surrounding the PAA brush greatly reduced nonspecific adsorption. Avidin was covalently attached to PAA brushes, and biotin-tagged proteins could be immobilized through avidin-biotin interaction. Such an immobilization method, which is based on specific interactions, is expected to better retain protein functionality than direct covalent binding. Using biotin-tagged bovine serum albumin (BSA) as a model, a simple strategy was developed for immobilization of small biological molecules using BSA as linkages, while BSA can simultaneously block nonspecific interactions.

Protein pI shifts due to posttranslational modifications in the separation and characterization of proteins

Proteins from breast cancer cell lines are characterized using a 2-D liquid separation technique in which protein pI is used as the first-dimension separation parameter. To effect this protein pI separation, chromatofocusing(CF) is employed whereby a pH gradient is generated on-column using a weak anion exchange medium with the intact proteins fractionated and collected every 0.2 pH unit. It is demonstrated that the pI for expressed intact proteins as generated by CF is an important parameter for identification and characterization of the actual protein modifications occurring in the cancer cell. For most proteins, the experimentally determined pI is very close to that predicted by the databases. In other cases, however, where the pI is observed to be shifted from the expected value, it is shown that this shift is often correlated to protein modifications. The modifications that cause such shifts include truncations and deletions often observed in cancer cells or phosphorylations that can shift the pI by several pH units. It is also shown that the effects of phosphorylation on the observed shift can vary depending upon the protein and the amount of phosphorylation. Moreover, large changes in the pI are often observed for proteins with a pI above 7.0 upon phosphorylation, whereas little change is observed for proteins with a pI of approximately 5.0. The expressed protein's pI value thus becomes an important parameter together with the intact MW value, peptide map, and MS/MS results for identification of the presence and type of posttranslational modifications occurring in the cancer cell.

Simple proteomic checks for detecting noncoding RNA

Proper validation can accelerate sequence-based discovery of proteins and protein-coding genes. Databases currently contain a backlog of experimentally unverified gene models and tentative assignments of observed transcripts to coding or noncoding RNA. We present and apply a general principle, founded on base composition and the genetic code and validated here by bulk 2-D gels, that can improve the reliability of such classifications and of the algorithms or pipelines that lead to them.

InSilicoSpectro: An Open-Source Proteomics Library

We present a new proteomics open-source project, InSilicoSpectro, aimed at implementing recurrent computations that are necessary for proteomics data analysis. Illustrative examples are mass list file format conversions, protein sequence digestion, theoretical peptide and fragment mass computations, graphical display, matching with experimental data, isoelectric point estimation, and peptide retention time prediction. The project library is written in Perl, a widely used scripting language in bioinformatics, and it offers a unique framework of integrated objects to implement complex proteomics data analyses. For instance, only a few lines of code are required to digest a protein with fixed and variable modifications, label peptides with 18O, compute the fragmentation spectra and display their match with experimental spectra. We believe that InSilicoSpectro will be of great help to bioinformaticians, without detailed knowledge of proteomics specifics, and to mass spectrometrists with computer programming interest as well.

Expression, purification, and characterization of Clostridium botulinum type B light chain

A full-length synthetic gene encoding the light chain of botulinum neurotoxin serotype B, approximately 50 kDa (BoNT/B LC), has been cloned into a bacterial expression vector pET24a+. BoNT/B LC was expressed in Escherichia coli BL21.DE3.pLysS and isolated from the soluble fraction. The resultant protein was purified to homogeneity by cation chromatography and was determined to be >98% pure as assessed by SDS-polyacrylamide gel stained with SilverXpress and analyzed by densitometry. Mass spectroscopic analysis indicated the protein to be 50.8 kDa, which equaled the theoretically expected mass. N-terminal sequencing of the purified protein showed the sequence corresponded to the known reported sequence. The recombinant BoNT/B light chain was found to be highly stable, catalytically active, and has been used to prepare antisera that neutralizes against BoNT/B challenge. Characterization of the protein including pH, temperature, and the stability of the protein in the presence or absence of zinc is described within. The influence of pH differences, buffer, and added zinc on secondary and tertiary structure of BoNT/B light chain was analyzed by circular dichroism and tryptophan fluorescence measurements. Optimal conditions for obtaining maximum metalloprotease activity and stabilizing the protein for long term storage were determined. We further analyzed the thermal denaturation of BoNT/B LC as a function of temperature to probe the pH and added zinc effects on light chain stability. The synthetic BoNT/B LC has been found to be highly active on its substrate (vesicle associated membrane protein-2) and, therefore, can serve as a useful reagent for BoNT/B research.

“Plasmo2D”: An Ancillary Proteomic Tool to Aid Identification of Proteins from Plasmodium falciparum

Bioinformatics tools to aid gene and protein sequence analysis have become an integral part of biology in the post-genomic era. Release of the Plasmodium falciparum genome sequence has allowed biologists to define the gene and the predicted protein content as well as their sequences in the parasite. Using pI and molecular weight as characteristics unique to each protein, we have developed a bioinformatics tool to aid identification of proteins from Plasmodium falciparum. The tool makes use of a Virtual 2-DE generated by plotting all of the proteins from the Plasmodium database on a pI versus molecular weight scale. Proteins are identified by comparing the position of migration of desired protein spots from an experimental 2-DE and that on a virtual 2-DE. The procedure has been automated in the form of user-friendly software called "Plasmo2D". The tool can be downloaded from http://144.16.89.25/Plasmo2D.zip.

Lysozyme dimerization: Brownian dynamics simulation

The lysozyme dimerization reaction has been studied within the framework of encounter-complex (EC) formation theory using the MacroDox software package. Two types of energetically favorite ECs were determined. In the first of them, active-center amino acids of lysozyme take part in the complex formation or the second molecule blocks accessibility to active center sterically. Epitope amino-acid residues are involved in the complex of type II. The existence of both types of complexes does not contradict experimental data. Dimer-formation rate constants for different kinds of EC were calculated. Increasing the pH from 2.0 to 10.0 decreases the total positive lysozyme charge and eliminates the unfavorable repulsive electrostatic interaction. The rate constant of EC formation is inversely proportional to the protein total charge. The association rate constant was also enhanced by an increase of ionic strength that screened repulsive electrostatic interaction between positively charged proteins.

Isoelectric point determination of proteins and other macromolecules: oscillating method

A program written in Visual Basic has been developed to calculate the isoelectric point of proteins and other macromolecules bearing acid-basic residues. The pI value can be theoretically calculated with the precision required. The computer automatically supplies a representation of the charge of the protein versus pH values. The corresponding values can also be obtained, on command, in the form of table.

Cationic double-hydrophilic model networks: synthesis, characterization, modeling and protein adsorption studies

Group transfer polymerization (GTP) was used for the preparation of eight networks based on two hydrophilic monomers, 2-(dimethylamino)ethyl methacrylate (DMAEMA) and poly(ethylene glycol) methacrylate (PEGMA). Ethylene glycol dimethacrylate (EGDMA) served as the cross-linker, whereas 1,4-bis(methoxytrimethylsiloxymethylene)cyclohexane (MTSMC) was used as a bifunctional initiator. Seven of the networks had linear segments of accurate molecular weight between the cross-links, i.e., they were model networks, whereas the eighth was an equimolar randomly cross-linked network. Five of the seven model networks were based on ABA triblock copolymers with PEGMA midblocks and DMAEMA endblocks, in which the DMAEMA/PEGMA ratio was varied. The remaining two model networks were equimolar isomers, the one based on BAB triblocks (with a DMAEMA midblock) and the other based on the statistical copolymer. The degrees of swelling of all of the networks were measured as a function of pH and were found to increase below pH 7. The degrees of swelling at low pH values increased with the percentage of the DMAEMA monomer, which is ionized under these conditions. These swelling results were confirmed qualitatively by theoretical calculations. Finally, the pH-dependence of the adsorption of the proteins pepsin, bovine serum albumin, and lysozyme onto one of the model networks was studied.

JVirGel: Calculation of virtual two-dimensional protein gels

We developed JVirGel, a collection of tools for the simulation and analysis of proteomics data. The software creates and visualizes virtual two-dimensional (2D) protein gels based on the migration behaviour of proteins in dependence of their theoretical molecular weights in combination with their calculated isoelectric points. The utilization of all proteins of an organism of interest deduced from genes of the corresponding genome project in combination with the elimination of obvious membrane proteins permits the creation of an optimized calculated proteome map. The electrophoretic separation behaviour of single proteins is accessible interactively in a Java(TM) applet (small application in a web browser) by selecting a pI/MW range and an electrophoretic timescale of interest. The calculated pattern of protein spots helps to identify unknown proteins and to localize known proteins during experimental proteomics approaches. Differences between the experimentally observed and the calculated migration behaviour of certain proteins provide first indications for potential protein modification events. When possible, the protein spots are directly linked via a mouse click to the public databases SWISS-PROT and PRODORIC. Additionally, we provide tools for the serial calculation and visualization of specific protein properties like pH dependent charge curves and hydrophobicity profiles. These values are helpful for the rational establishment of protein purification procedures. The proteomics tools are available on the World Wide Web at http://prodoric.tu-bs.de/proteomics.php.

Identification of extrinsic blood coagulation pathway inhibitors from the tick Ornithodoros savignyi (Acari: Argasidae)

The salt BaSO(4) selectively adsorbs two proteins from crude Ornithodoros savignyi salivary gland extract. They co-purify during reversed-phase HPLC, but can be separated by hydrophobic-interaction chromatography. Their molecular masses are 9333 and 9173Da. The 9.3kDa protein was designated BSAP1 and the 9.1kDa protein BSAP2. Their amino acid compositions show significant differences, in particular the presence of seven and eight cysteine residues in BSAP1 and BSAP2, respectively. The proteins do not contain gamma-carboxyglutamic acid, hydroxyproline, or hydroxylysine. The proteins do not inhibit the intrinsic coagulation cascade, but inhibit the extrinsic pathway. The observed inhibition is not due to inhibition of factor VII. Both proteins bind to membranes. BSAP1 binds neutral and negatively charged membranes more strongly than BSAP2. Its affinity for negative membranes is, however, much lower than for neutral membranes. In contrast, BSAP2 binds both membranes equally strongly. The binding of the proteins to the membranes was significantly lowered upon pre-incubation with Ca(2+).