The QUEST system for quantitative analysis of two-dimensional gels

The strategies and methods used by the QUEST system for two-dimensional gel analysis are described, and the performance of the system is evaluated. Radiolabeled proteins, resolved on two-dimensional gels and detected using calibrated exposures to film, are quantified in units of disintegrations per minute or as a fraction of the total protein radioactivity applied to the gel. Spot quantitation and resolution of overlapping spots is performed by two-dimensional gaussian fitting. Pattern matching is carried out for groups of gels called matchsets, and within each matchset every gel is matched to every other gel. During the matching process, spots are automatically added to each pattern at positions where unmatched spots were detected in other patterns. This results in enhanced accuracy for both spot detection and for matching. The spot fitting procedure is repeated after matching. Tests show that up to 97% of spots in each pattern can be matched and that fewer than 1% of the spots are matched inconsistently. Approximately 2000 proteins are detected from typical gels. Of these 1600 are high quality spots. Tests to measure the coefficient of variation of spot quantitation versus spot quality show that the average coefficient of variation for high quality spots is 21%. The intensities of the detected proteins range from 4 to 20,000 ppm of total protein synthesis. The QUEST analysis system has been used to build a quantitative database for the proteins of normal and transformed REF52 cells, as presented in the accompanying reports (Garrels, J., and Franza, B. R., Jr. (1989) J. Biol. Chem. 264, 5283-5298, 5299-5312).     

A combined proteomic and genetic analysis of the highly variable platelet proteome: From plasmatic proteins and SNPs

High biological variation in protein expression represents a major challenge in clinical proteomics. In a study based on 2D-DIGE, we found that the standardised abundance of only a few proteins varied by more than 50%. While some of the highest variable proteins in platelets of 52 healthy elderly were of plasmatic origin, such as albumin or haptoglobin, absence of several other high-abundant plasma proteins strongly suggests that plasma-derived proteins represent an integral part of the platelet proteome. Amongst the highly variable platelet-derived proteins, two spots were both identified as GSTO1 and assigned to either the wild-type or mutant isoform of SNP A140D. Remarkably, when the spots were considered within the respective genotype groups, their CV decreased to about the median variation. Albeit 2D-DIGE allowed correct genotyping, two individuals seemed to be GSTO1*A140 deficient. Probing 2D-Western blots with novel mAb, however, detected A140 protein as additional spot at pH 8.1, caused by the SNPs E155del and E208K. In contrast to previous studies, we show that GSTO1 protein is expressed in vivo, despite the deletion E155. Our data indicate that incorporation of exogenous proteins and genetic polymorphisms of endogenous proteins represent the main source of extreme biological variation in the platelet proteome.     

Interindividual Variation in the Proteome of Human Peripheral Blood Mononuclear Cells

Peripheral blood mononuclear cells (PBMCs) are main actors in inflammatory processes and linked to many diseases, including rheumatoid arthritis, atherosclerosis, asthma, HIV and cancer. Moreover, they seem an interesting ‘surrogate tissue’ that can be used in biomarker discovery. In order to get a good experimental design for quantitative expression studies, the knowledge of the interindividual variation is an essential part. Therefore, PBMCs were isolated from 24 healthy volunteers (15 males, 9 females, ages 63–86) with no clinical signs of inflammation. The extracted proteins were separated using the two dimensional difference in gel electrophoresis technology (2D-DIGE), and the gel images were processed with the DeCyder 2D software. Protein spots present in at least 22 out of 24 healthy volunteers were selected for further statistical analysis. Determination of the coefficient of variation (CV) of the normalized spot volume values of these proteins, reveals that the total variation of the PBMC proteome varies between 12,99% to 148,45%, with a mean value of 28%. A supplemental look at the causes of technical variation showed that the isolation of PBMCs from whole blood is the factor which influences the experimental variance the most. This isolation should be handled with extra care and an additional washing step would be beneficial. Knowing the extent of variation, we show that at least 10 independent samples per group are needed to obtain statistical powerful data. This study demonstrates the importance of considering variance of a human population for a good experimental design for future protein profiling or biomarker studies.     

Assessing estradiol in biobehavioral studies using saliva and blood spots: simple radioimmunoassay protocols, reliability, and comparative validity

We developed simple, reliable, and highly sensitive assay modifications of commercially available radioimmunoassay kits to measure estradiol in saliva and blood spot specimens. The saliva assay has average intra- and interassay coefficients of variation (CV) of 6.45 and 9.01%, with average analytical and serial dilution recoveries 100.65 and 89.25%. The blood spot assay has average intra- and interassay CVs of 7.57 and 8.22%, with analytical and serial dilution recoveries of 80.50 and 108.50%. The analytical sensitivity ranges of the saliva (0.25-7.50 pg/ml) and blood spot (2. 00-375 pg/ml) assays are sufficient to determine levels in the majority of pre- and postpubertal males and females. Blood spot assay results are correlated with serum estradiol levels for adult males, r (17) = 0.73, and females, r (18) = 0.96. In contrast, the serum-saliva correlation is only modest for adult females, r (14) = 0.60, and not significant for adult males. Substitution of blood spot assay results for serum values underestimates the known serum estradiol-behavior correlation by only 3.45%, whereas substitution of saliva assay results for serum values underestimates the association by 37.55%. The findings have important implications for the use and potential misuse of noninvasive measures of estradiol in studies of health and human development.     

Comparative proteomics analysis of human lung squamous carcinoma

Two-dimensional polyacrylamide gel electrophoresis (2-DE) profiles of human lung squamous carcinoma tissue and paired surrounding normal bronchial epithelial tissue were compared. Selected differential protein-spots were identified with peptide mass fingerprinting based on matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) and database searching. Well-resolved and reproducible 2-DE patterns of both the tumor and the normal tissues were acquired. The average deviations of spot position were 0.873+/-0.125mm in IEF direction and 1.025+/-0.213mm in SDS-PAGE direction, respectively. For the tumor tissues, a total of 1349+/-67 spots were detected and 1235+/-48 spots were matched with an average matching rate of 91.5%. For the corresponding normal tissues, a total of 1297+/-73 spots were detected and 1183+/-56 spots were matched with an average matching rate of 91.2%. A total of 1069+/-45 spots were matched between the tumor and the normal tissues. Forty differential proteins between tumor and normal tissues were characterized. Some proteins were the products of oncogenes and others were involved in the regulation of cell cycle and signal transduction. These data are valuable for mass identification of differentially expressed proteins involved in lung carcinogenesis, establishing human lung cancer proteome database and screening molecular marker to further study human lung squamous carcinoma.     

Variation in protein levels obtained from human blood cells and biofluids for platelet, peripheral blood mononuclear cell, plasma, urine and saliva proteomics

Blood cells and biofluid proteomics are emerging as a valuable tool to assess effects of interventions on health and disease. This study is aimed to assess the amount and variability of proteins from platelets, peripheral blood mononuclear cells (PBMC), plasma, urine and saliva from ten healthy volunteers for proteomics analysis, and whether protein yield is affected by prolonged fasting. Volunteers provided blood, saliva and morning urine samples once a week for 4 weeks after an overnight fast. Volunteers were fasted for a further 24 h after the fourth sampling before providing their final samples. Each 10 mL whole blood provided 400–1,500 μg protein from platelets, and 100–600 μg from PBMC. 30 μL plasma depleted of albumin and IgG provided 350–650 μg protein. A sample of morning urine provided 0.9–8.6 mg protein/dL, and a sample of saliva provided 70–950 μg protein/mL. None of these yields were influenced by the degree of fasting (overnight or 36 h). In conclusion, in contrast to the yields from plasma, platelets and PBMC, the protein yields of urine and saliva samples were highly variable within and between subjects. Certain disease conditions may cause higher or lower PBMC counts and thus protein yields, or increased urinary protein levels.     

Proteomic comparison of two-dimensional gel electrophoresis profiles from human lung squamous carcinoma and normal bronchial epithelial tissues

Differential proteome profiles of human lung squamous carcinoma tissue compared to paired tumor-adjacent normal bronchial epithelial tissue were established and analyzed by means of immobilized pH gradient-based two-dimensional polyacrylamide gel electrophoresis (2-D PAGE) and matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF-MS). The results showed that well-resolved, reproducible 2-DE patterns of human lung squamous carcinoma and adjacent normal bronchial epithelial tissues were obtained under the condition of 0.75-mg protein-load. The average deviation of spot position was 0.733+/-0.101 mm in IEF direction, and 0.925+/-0.207 mm in SDS-PAGE direction. For tumor tissue, a total of 1241+/-88 spots were detected, 987+/-65 spots were matched with an average matching rate of 79.5%. For control, a total of 1190+/-72 spots were detected, and 875+/-48 spots were matched with an average matching rate of 73.5%. A total of 864+/-34 spots were matched between tumors and controls. Forty-three differential proteins were characterized: some proteins were related to oncogenes, and others involved in the regulation of cell cycle and signal transduction. It is suggested that the differential proteomic approach is valuable for mass identification of differentially expressed proteins involved in lung carcinogenesis. These data will be used to establish human lung cancer proteome database to further study human lung squamous carcinoma.     

Matrix protein microarrays for spatially and compositionally controlled microspot thrombosis under laminar flow

Microarraying allows the spatial and compositional control of surfaces, typically for the purpose of binding reactions. Collagen and/or von Willebrand Factor (vWF) in 5% glycerol was contact printed onto glass slides to create defined microspots (176-microm diameter) of adsorbed protein without sample dehydration. The arrays were mounted on flow chambers allowing video microscopy during perfusion (wall shear rate of 100-500 s(-1)) of recalcified corn trypsin inhibitor-treated whole blood or platelet rich plasma and subsequent array scanning via anti-GPIbalpha and anti-fibrin(ogen) immunofluorescence. To mimic the subendothelial matrix, vWF was microarrayed over sonicated type I collagen microspots. For whole blood perfusion (500 s(-1), 10 min) over collagen, vWF, and collagen/vWF microspots, the amount of platelet deposition on the collagen/vWF spots was approximately 2 times greater in comparison to the collagen spots and approximately 18 times greater in comparison to the vWF spots. The amount of fibrin(ogen) deposition on the collagen/vWF spots was approximately 2 times greater in comparison to the collagen spots and approximately 4 times greater in comparison to the vWF spots. This protocol allowed for highly uniform (CV = 18%) and precisely located thrombus formation at a density of >or=400 spots/cm(2). Microarrays are ideal for the combinatorial assembly of adhesive and procoagulant proteins to study thrombosis as well as to study axial and lateral transport effects between discrete microspots of distinct composition.     

Quantification of cationic anti-malaria agent methylene blue in different human biological matrices using cation exchange chromatography coupled to tandem mass spectrometry

Selective and sensitive methods for the determination of the cationic dye and anti-malarial methylene blue in human liquid whole blood, dried whole blood (paper spot), and plasma depending on protein precipitation and cation exchange chromatography coupled to electrospray ionisation (ESI) tandem mass spectrometry (MS/MS) have been developed, validated according to FDA standards, and applied to samples of healthy individuals and malaria patients within clinical studies. Acidic protein precipitation with acetonitrile and trifluoroacetic acid was used for liquid whole blood and plasma. For the extraction of methylene blue from paper spots aqueous acetonitrile was used. Sample extracts were chromatographed on a mixed mode column (cation exchange/reversed phase, Uptisphere MM1) using an aqueous ammonium acetate/acetonitrile gradient. Methylene blue was quantified with MS/MS in the selected reaction monitoring mode using ESI and methylene violet 3RAX as internal standard. Depending on the sample volume (whole blood and plasma 250 microL, and 100 microL on paper spots) the method was linear at least within 75 and 10,000 ng/mL and the limit of quantification in all matrices was 75 ng/mL. Batch-to-batch accuracies of the whole blood, plasma, and paper spot methods varied between -4.5 and +6.6%, -3.7 and +7.5%, and -5.8 and +11.1%, respectively, with corresponding precision ranging from 3.8 to 11.8% CV. After a single oral dose (500 mg) methylene blue concentrations were detectable for 72 h in plasma. The methods were applied within clinical studies to samples from healthy individuals and malaria patients from Burkina Faso.     

Protein subunit interfaces: a statistical analysis of hot spots in Sm proteins

The distinguishing property of Sm protein associations is very high stability. In order to understand this property, we analyzed the interfaces and compared the properties of Sm protein interfaces with those of a test set, the Binding Interface Database (BID). The comparison revealed that the main differences between the interfaces of Sm proteins and those of the BID set are the content of charged residues, the coordination numbers of the residues, knowledge-based pair potentials, and the conservation scores of hot spots. In Sm proteins, the interfaces have more hydrophobic and fewer charged residues than the surfaces, which is also the case for the BID test set and other proteins. However, in the interfaces, the content of charged residues in Sm proteins (26%) is substantially larger than that in the BID set (22%). Hot spots are residues that make up a small fraction of the interfaces, but they contribute most of the binding energy. These residues are critical to protein–protein interactions. Analyses of knowledge-based pair potentials of hot spot and non-hot spot residues in Sm proteins show that they are significantly different; their mean values are 31.5 and 11.3, respectively. In the BID set, this difference is smaller; in this case, the mean values for hot spot and non-hot spot residues are 20.7 and 12.4, respectively. Hence, the pair potentials of hot spots differ significantly for the Sm and BID data sets. In the interfaces of Sm proteins, the amino acids are tightly packed, and the coordination numbers are larger in Sm proteins than in the BID set for both hot spots and non-hot spots. At the same time, the coordination numbers are higher for hot spots; the average coordination number of the hot spot residues in Sm proteins is 7.7, while it is 6.1 for the non-hot spot residues. The difference in the calculated average conservation score for hot spots and non-hot spots in Sm proteins is significantly larger than it is in the BID set. In Sm proteins, the average conservation score for the hot spots is 7.4. Hot spots are surrounded by residues that are moderately conserved (5.9). The average conservation score for the other interface residues is 5.6. The conservation scores in the BID set do not show a significant distinction between hot and non-hot spots: the mean values for hot and non-hot spot residues are 5.5 and 5.2, respectively. These data show that structurally conserved residues and hot spots are significantly correlated in Sm proteins.     

Improved image analysis workflow for 2-D gels enables large-scale 2-D gel-based proteomics studies--COPD biomarker discovery study

2-D gel electrophoresis has been used for more than three decades to study the protein complement of organisms, tissues, and cells. Three issues are holding back large-scale proteomics studies: low-throughput, high technical variation, and study designs lacking statistical power. We identified image analysis as the central factor connecting these three issues. By developing an improved image analysis workflow we shortened project timelines, decreased technical variation, and thus enabled large-scale proteomics studies that are statistically powered. Rather than detecting protein spots on each gel image and matching spots across gel images, the improved workflow is based on aligning images first, then creating a consensus spot pattern and finally propagating the consensus spot pattern to all gel images for quantitation. This results in a data table without gaps. As an example we show here a study aimed at discovering circulating biomarkers for chronic obstructive pulmonary disease (COPD). Eight candidate biomarkers were identified by comparing plasma from 24 smokers with COPD and 24 smokers without COPD. Among the candidates are proteins such as plasma retinal-binding protein (RETB) and fibrinogen that had previously been linked to the disease and are frequently monitored in COPD patients, as well as other proteins such as apolipoprotein E (ApoE), inter-alpha-trypsininhibitor heavy chain H4 (ITIH4), and glutathione peroxidase.     

Alternative and effective proteomic analysis in Arabidopsis

Various functional genomics platforms are required to define the phenotype associated with a mutant. Global protein analyses may be included in any study. We describe here a rapid method of protein sample preparation and analysis, suitable for all laboratories and using Arabidopsis plantlets as the starting material. This reliable and reproducible method for high yield protein extraction from small amounts of material can be used on even the most recalcitrant tissues. The proteins extracted are suitable for many types of protein analysis, including nondenaturing investigations. This method was validated by a rigorous 2-DE approach, coupled with unambiguous LC-MS/MS identifications featuring strong sequence coverage (average of 26% with eight different peptides/spot protein). The reproducibility of the method was demonstrated by multiple protein identifications from identical series of spots. An interactive map (http://www.isv.cnrsgif.fr/gel2d/), including 435 protein variants showed that (i) 38% of the proteins were yet unreported, (ii) reduced subfractionation, (iii) had frequent protein modifications (average of two spots/protein entry), and (iv) underwent no major proteolytic events other than leader peptide cleavage. Finally, a simple mobility shift method for the large subunit of RuBisCo (LS) in the first dimension made it possible to characterize previously masked protein spots.     

Evaluation of individual protein errors in silver-stained two-dimensional gels

The relationship between spot volume and variation for all protein spots observed on large format 2D gels when utilising silver stain technology and a model system based on mammalian NSO cell extracts is reported. By running multiple gels we have shown that the reproducibility of data generated in this way is dependent on individual protein spot volumes, which in turn are directly correlated with the coefficient of variation. The coefficients of variation across all observed protein spots were highest for low abundant proteins which are the primary contributors to process error, and lowest for more abundant proteins. Using the relationship between spot volume and coefficient of variation we show it is necessary to calculate variation for individual protein spot volumes. The inherent limitations of silver staining therefore mean that errors in individual protein spot volumes must be considered when assessing significant changes in protein spot volume and not global error.     

Proteomic analysis of blood cells in fish exposed to chemotherapeutics: evidence for long term effects

Proteomics technology are increasingly used in ecotoxicological studies to characterize and monitor biomarkers of exposure. The present study aims at identifying long term effects of malachite green (MG) exposure on the proteome of peripheral blood mononuclear cells (PBMC) from the Asian catfish, Pangasianodon hypophthalmus. A common (0.1 ppm) concentration for therapeutic treatment was applied twice with a 72 h interval. PBMC were collected directly at the end of the second bath of MG (T1) and after 1 month of decontamination (T2). Analytical 2D-DIGE gels were run and a total of 2551±364 spots were matched. Among them, MG induced significant changes in abundance of 116 spots with no recovery after one month of decontamination. Using LC-MS/MS and considering single identification per spot, we could identify 25 different proteins. Additionally, MG residues were measured in muscle and in blood indicating that leuco-MG has almost totally disappeared after one month of decontamination. This work highlights long term effects of MG treatment on the PBMC proteome from fish intended for human consumption.     

Changes in the blood plasma protein composition in an experiment with seven-day dry immersion

This study was designed to track changes in the blood plasma proteome caused by exposure of healthy humans to dry immersion. The blood plasma specimens of five healthy volunteers aged 23–29 years obtained seven days before immersion, on day 7 of exposure in the immersion bath, and seven days after completion of the experiment were analyzed. After the extraction of major and the concentration of minor proteins, depleted blood plasma fractions were subjected to two-dimensional electrophoresis with the subsequent identification of significantly different protein spots with the method of mass-spectrometric analysis of peptide fragments. A significant change in intensity was observed for 25 protein spots with a coefficient of variation of >60%. Cluster analysis confirmed substantial changes in the blood plasma proteome in the period of readaptation of the body to the conditions of normal vital activity upon completion of immersion. On day 7 of readaptation, the apolipoprotein A-1, A-IV, and E concentrations were significantly higher compared to the baseline data and the measurements on day 7 of immersion. However, the levels of α-, β-fibrinogen, the complement C4B factor, and serum amyloid P were found to be lower.     

Triton X-114 cloud point extraction to subfractionate blood plasma proteins for two-dimensional gel electrophoresis

A simple and reproducible procedure for enrichment of a plasma protein subfraction suitable for two-dimensional polyacrylamide gel electrophoresis (2DE) was developed, using a Triton X-114-based cloud point extraction (CPE). Appropriate conditions for such a CPE procedure were found by SDS-PAGE to be a plasma protein concentration of about 10 mg/ml in 3% (w/v) Triton X-114. 2DE of proteins obtained by CPE of 400 μl of human plasma revealed about 200 spots constituting a spot pattern very different from the pattern of total plasma. The CPE procedure only had a limited contribution to the technical variation. Identification of about 60 spots, representing only 22 proteins, revealed that several proteins in the obtained subfraction were present in more isoforms or modifications. Among these were apolipoproteins (A-1, D, E, L1, and M), haptoglobin-related protein, phosphatidylcholine-sterol acyltransferase, serum amyloid A, and serum paraoxonase/arylesterase 1, which are proteins of a hydrophobic nature, as in plasma they relate to lipoprotein particles. Thus, Triton X-114-based CPE is a simple plasma prefractionation tool, attractive for detailed 2DE studies of hydrophobic plasma proteins and their isoforms or modifications.     

Optimizing search conditions for the mass fingerprint-based identification of proteins

The two central problems in protein identification by searching a protein sequence collection with MS data are the optimal use of experimental information to allow for identification of low abundance proteins and the accurate assignment of the probability that a result is false. For comprehensive MS-based protein identification, it is necessary to choose an appropriate algorithm and optimal search conditions. We report a systematic study of the quality of PMF-based protein identifications under different sequence collection search conditions using the Probability algorithm, which assigns the statistical significance to each result. We employed 2244 PMFs from 2-DE-separated human blood plasma proteins, and performed identification under various search constraints: mass accuracy (0.01-0.3 Da), maximum number of missed cleavage sites (0-2), and size of the sequence collection searched (5.6 x 10(4)-1.8 x 10(5)). By counting the number of significant results (significance levels 0.05, 0.01, and 0.001) for each condition, we demonstrate the search condition impact on the successful outcome of proteome analysis experiments. A mass correction procedure utilizing mass deviations of albumin matching peptides was tested in an attempt to improve the statistical significance of identifications and iterative searching was employed for identification of multiple proteins from each PMF.     

Using standard positions and image fusion to create proteome maps from collections of two-dimensional gel electrophoresis images

Databases for two-dimensional protein gels pose new challenges in extracting meaningful information from large numbers of experiments. In order to create expression profiles, positions of corresponding protein spots across all gel images have to be established. In larger gel sets errors may accumulate rapidly during this spot matching process, effectively limiting the number of samples available for data mining. Here we present a novel approach for organizing spot data based on the concept of a standard position for a protein species. Standard positions are meaningful average positions that are determined using all occurrences of a protein species. They can be extended to spots that are not annotated via interpolation. The standard position of a spot can serve as a unifying index across all gels in a database, thus allowing creation and analysis of expression profiles that span the whole collection. The standard position gives a much more accurate estimation of a spot's position on a gel than can be obtained using theoretical isoelectric point and molecular weight. Positional indexing is a complement to a priori identifications (e.g. by mass spectrometry or Edman degradation). Moreover it can be used in advance to select spots that are worth identifying because they show relevant expression profiles. Furthermore, we show how to combine all spots that occur on any of the gels into one synthetic but nevertheless realistic-looking image. This composite image is produced such that all spots have their standard positions. It can serve as a proteome reference map for an organism. As an application, we have computed a reference map from 23 gel images of Bacillus subtilis, using an enhanced prerelease version of the gel analysis software Delta2D (DECODON, Greifswald, Germany).     

Proteomic analysis of near-isogenic sunflower varieties differing in seed oil traits

Near-isogenic sunflower lines containing 25% (inbred RHA280) and 48% (RHA801) oil by seed dry mass were comparatively analyzed in biological triplicate at 18 days after flowering using two-dimensional (both pI 3-10 and 4-7) Difference Gel Electrophoresis. Additionally, two inbred lines varying in oleic acid content, HA89 (18% oleic) and HA341 (89% oleic), were also analyzed in the same manner. Statistical analyses of these sunflower lines was performed beginning with fitting a mixed effects linear model to the log-transformed optical volume of each spot to account for gel variation, followed by testing the significance between varieties for mean transformed optical spot volumes. The p-values from the spot analysis procedures were then used to find the cutoff point for differential expression using a 10% false-discovery rate (FDR). Comparison of the oil content and oleic acid composition lines revealed 77 and 42 protein spots below the 10% FDR cutoff, respectively, and were therefore declared differentially expressed. Liquid chromatography-tandem mass spectrometry analysis of each of these protein spots resulted in assignments for 44 and 17 spots, respectively. Fructokinase, plastid phosphoglycerate kinase, and enolase proteins were determined to be up-regulated in the high oil line, while phosphofructokinase, cytosolic phosphoglucomutase, and cytsolic phosphoglycerate kinase were up-regulated in the low oil variety. Additionally, four activities involved in amino acid synthesis were up-regulated in the low oil variety in addition to 12S storage proteins and a protein similar to legumin storage protein. Interestingly, two 2-DE spots identified as 14-3-3 proteins were found to be up-regulated in high oleic acid variety. Alteration of glycolytic and amino acid biosynthetic enzymes, as well as storage protein levels, suggests seed oil content is tightly linked to carbohydrate metabolism and protein synthesis in a complex manner.     

Hierarchical grid transformation for image warping in the analysis of two-dimensional electrophoresis gels

Two-dimensional electrophoresis is a widely used method for separating a large number of proteins from complex protein mixtures and for revealing differential patterns of protein expressions. In the computer-assisted proteome research, the comparison of protein separation profiles involves several heuristic steps, ranging from protein spot detection to matching of unknown spots. An important prerequisite for efficient protein spot matching is the image warping step, where the geometric relationship between the gel profiles is modeled on the basis of a given set of known corresponding spots, so-called landmarks, and the locations of unknown spots are predicted using the optimized model. Traditionally, polynomial functions together with least squares optimization has been used, even though this approach is known to be incapable of modeling all the complex distortions inherent in electrophoretic data. To satisfy the need of more flexible gel distortion correction, a hierarchical grid transformation method with stochastic optimization is presented. The method provides an adaptive multiresolution model between the gels, and good correction performance in the practical cross-validation tests suggests that automatic warping of gel images could be based on this approach. We believe that the proposed model also has significance in the ultimate comparison of corresponding protein spots since the matching process should benefit from the closeness of the true spot pairs.