Linking genome and proteome by mass spectrometry: Large-scale identification of yeast proteins from two dimensional gels

The function of many of the uncharacterized open reading frames discovered by genomic sequencing can be determined at the level of expressed gene products, the proteome. However, identifying the cognate gene from minute amounts of protein has been one of the major problems in molecular biology. Using yeast as an example, we demonstrate here that mass spectrometric protein identification is a general solution to this problem given a completely sequenced genome. As a first screen, our strategy uses automated laser desorption ionization mass spectrometry of the peptide mixtures produced by in-gel tryptic digestion of a protein. Up to 90% of proteins are identified by searching sequence data bases by lists of peptide masses obtained with high accuracy. The remaining proteins are identified by partially sequencing several peptides of the unseparated mixture by nanoelectrospray tandem mass spectrometry followed by data base searching with multiple peptide sequence tags. In blind trials, the method led to unambiguous identification in all cases. In the largest individual protein identification project to date, a total of 150 gel spots—many of them at subpicomole amounts—were successfully analyzed, greatly enlarging a yeast two-dimensional gel data base. More than 32 proteins were novel and matched to previously uncharacterized open reading frames in the yeast genome. This study establishes that mass spectrometry provides the required throughput, the certainty of identification, and the general applicability to serve as the method of choice to connect genome and proteome.

The high mobility group protein Abf2p influences the level of yeast mitochondrial DNA recombination intermediates in vivo

Abf2p is a high mobility group (HMG) protein found in yeast mitochondria that is required for the maintenance of wild-type (ρ+) mtDNA in cells grown on fermentable carbon sources, and for efficient recombination of mtDNA markers in crosses. Here, we show by two-dimensional gel electrophoresis that Abf2p promotes or stabilizes Holliday recombination junction intermediates in ρ+ mtDNA in vivo but does not influence the high levels of recombination intermediates readily detected in the mtDNA of petite mutants (ρ−). mtDNA recombination junctions are not observed in ρ+ mtDNA of wild-type cells but are elevated to detectable levels in cells with a null allele of the MGT1 gene (Δmgt1), which codes for a mitochondrial cruciform-cutting endonuclease. The level of recombination intermediates in ρ+ mtDNA of Δmgt1 cells is decreased about 10-fold if those cells contain a null allele of the ABF2 gene. Overproduction of Abf2p by ≥ 10-fold in wild-type ρ+ cells, which leads to mtDNA instability, results in a dramatic increase in mtDNA recombination intermediates. Specific mutations in the two Abf2p HMG boxes required for DNA binding diminishes these responses. We conclude that Abf2p functions in the recombination of ρ+ mtDNA.

Identification of a new mammalian centrin gene, more closely related to Saccharomyces cerevisiae CDC31 gene

Among the numerous centrin isoforms identified by two-dimensional gel electrophoresis in human cells, an acidic and slow-migrating isoform is particularly enriched in a centrosome fraction. We report here that this isoform specifically reacts with antibodies raised against Saccharomyces cerevisiae Cdc31p and is present, as other centrin isoforms, in the distal lumen of centrioles. It is encoded by a new centrin gene, which we propose to name HsCEN3 (Homo sapiens centrin gene 3). This gene is more closely related to the yeast CDC31 gene, and shares less identity with algae centrin than HsCEN1 and HsCEN2. A murine CDC31-related gene was also found that shows 98% identity and 100% similarity with HsCEN3, demonstrating a higher interspecies conservation than the murine centrin gene MmCEN1 (Mus musculus centrin gene 1) with either HsCEN1, or HsCEN2. Finally, immunological data suggest that a CDC31-related gene could exist in amphibians and echinoderms as well. All together, our data suggest the existence of two divergent protein subfamilies in the current centrin family, which might be involved in distinct centrosome-associated functions. The possible implication of this new mammalian centrin gene in centrosome duplication is discussed.

Short DNA Fragments without Sequence Similarity Are Initiation Sites for Replication in the Chromosome of the Yeast Yarrowia lipolytica

We have previously shown that both a centromere (CEN) and a replication origin are necessary for plasmid maintenance in the yeast Yarrowia lipolytica (Vernis et al., 1997). Because of this requirement, only a small number of centromere-proximal replication origins have been isolated from Yarrowia. We used a CEN-based plasmid to obtain noncentromeric origins, and several new fragments, some unique and some repetitive sequences, were isolated. Some of them were analyzed by two-dimensional gel electrophoresis and correspond to actual sites of initiation (ORI) on the chromosome. We observed that a 125-bp fragment is sufficient for a functional ORI on plasmid, and that chromosomal origins moved to ectopic sites on the chromosome continue to act as initiation sites. These Yarrowia origins share an 8-bp motif, which is not essential for origin function on plasmids. The Yarrowia origins do not display any obvious common structural features, like bent DNA or DNA unwinding elements, generally present at or near eukaryotic replication origins. Y. lipolytica origins thus share features of those in the unicellular Saccharomyces cerevisiae and in multicellular eukaryotes: they are discrete and short genetic elements without sequence similarity.

Proteomic analysis of the bacterial cell cycle

A global approach was used to analyze protein synthesis and stability during the cell cycle of the bacterium Caulobacter crescentus. Approximately one-fourth (979) of the estimated C. crescentus gene products were detected by two-dimensional gel electrophoresis, 144 of which showed differential cell cycle expression patterns. Eighty-one of these proteins were identified by mass spectrometry and were assigned to a wide variety of functional groups. Pattern analysis revealed that coexpression groups were functionally clustered. A total of 48 proteins were rapidly degraded in the course of one cell cycle. More than half of these unstable proteins were also found to be synthesized in a cell cycle-dependent manner, establishing a strong correlation between rapid protein turnover and the periodicity of the bacterial cell cycle. This is, to our knowledge, the first evidence for a global role of proteolysis in bacterial cell cycle control.

Molecular Cloning and Expression Analysis of the Mitochondrial Pyruvate Dehydrogenase from Maize1

Four cDNAs, one encoding an α-subunit and three encoding β-subunits of the mitochondrial pyruvate dehydrogenase, were isolated from maize (Zea mays L.) libraries. The deduced amino acid sequences of both α- and β-subunits are approximately 80% identical with Arabidopsis and pea (Pisum sativum L.) homologs. The mature N terminus was determined for the β-subunit by microsequencing the protein purified from etiolated maize shoot mitochondria and was resolved by two-dimensional gel electrophoresis. This single isoelectric species comprised multiple isoforms. Both α- and β-subunits are encoded by multigene families in maize, as determined by Southern-blot analyses. RNA transcripts for both α- and β-subunits were more abundant in roots than in young leaves or etiolated shoots. Pyruvate dehydrogenase activity was also higher in roots (5-fold) compared with etiolated shoots and leaves. Both subunits were present at similar levels in all tissues examined, indicating coordinated gene regulation. The protein levels were highest in heterotrophic organs and in pollen, which contained about 2-fold more protein than any other organ examined. The relative abundance of these proteins in nonphotosynthetic tissues may reflect a high cellular content of mitochondria, a high level of respiratory activity, or an extra plastidial requirement for acetate.

Evidence that the 50-kDa substrate of brefeldin A-dependent ADP-ribosylation binds GTP and is modulated by the G-protein beta gamma subunit complex.

Brefeldin A, a fungal metabolite that inhibits membrane transport, induces the mono(ADP-ribosyl)ation of two cytosolic proteins of 38 and 50 kDa as judged by SDS/PAGE. The 38-kDa substrate has been previously identified as glyceraldehyde-3-phosphate dehydrogenase (GAPDH). We report that the 50-kDa BFA-induced ADP-ribosylated substrate (BARS-50) has native forms of 170 and 130 kDa, as determined by gel filtration of rat brain cytosol, indicating that BARS-50 might exist as a multimeric complex. BARS-50 can bind GTP, as indicated by blot-overlay studies with [alpha-32P]GTP and by photoaffinity labeling with guanosine 5'-[gamma-32P] [beta,gamma-(4-azidoanilido)]triphosphate. Moreover, ADP-ribosylation of BARS-50 was completely inhibited by the beta gamma subunit complex of G proteins, while the ADP-ribosylation of GAPDH was unmodified, indicating that this effect was due to an interaction of the beta gamma complex with BARS-50, rather than with the ADP-ribosylating enzyme. Two-dimensional gel electrophoresis and immunoblot analysis shows that BARS-50 is a group of closely related proteins that appear to be different from all the known GTP-binding proteins.

In vivo estradiol-dependent dephosphorylation of the repressor MDBP-2-H1 correlates with the loss of in vitro preferential binding to methylated DNA.

We have previously shown that estradiol treatment of roosters resulted in a rapid loss of binding activity of the repressor MDBP-2-H1 (a member of the histone H1 family) to methylated DNA that was not due to a decrease in MDBP-2-H1 concentration. Here we demonstrate that MDBP-2-H1 from rooster liver nuclear extracts is a phosphoprotein. Phosphoamino acid analysis reveals that the phosphorylation occurs exclusively on serine residues. Two-dimensional gel electrophoresis and tryptic phosphopeptide analysis show that MDBP-2-H1 is phosphorylated at several sites. Treatment of roosters with estradiol triggers a dephosphorylation of at least two sites in the protein. Phosphatase treatment of purified rooster MDBP-2-H1 combined with gel mobility shift assay indicates that phosphorylation of MDBP-2-H1 is essential for the binding to methylated DNA and that the dephosphorylation can occur on the protein bound to methylated DNA causing its release from DNA. Thus, these results suggest that in vivo modification of the phosphorylation status of MDBP-2-H1 caused by estradiol treatment may be a key step for the down regulation of its binding to methylated DNA.

Rapid mass spectrometric peptide sequencing and mass matching for characterization of human melanoma proteins isolated by two-dimensional PAGE.

We report a general mass spectrometric approach for the rapid identification and characterization of proteins isolated by preparative two-dimensional polyacrylamide gel electrophoresis. This method possesses the inherent power to detect and structurally characterize covalent modifications. Absolute sensitivities of matrix-assisted laser desorption ionization and high-energy collision-induced dissociation tandem mass spectrometry are exploited to determine the mass and sequence of subpicomole sample quantities of tryptic peptides. These data permit mass matching and sequence homology searching of computerized peptide mass and protein sequence data bases for known proteins and design of oligonucleotide probes for cloning unknown proteins. We have identified 11 proteins in lysates of human A375 melanoma cells, including: alpha-enolase, cytokeratin, stathmin, protein disulfide isomerase, tropomyosin, Cu/Zn superoxide dismutase, nucleoside diphosphate kinase A, galaptin, and triosephosphate isomerase. We have characterized several posttranslational modifications and chemical modifications that may result from electrophoresis or subsequent sample processing steps. Detection of comigrating and covalently modified proteins illustrates the necessity of peptide sequencing and the advantages of tandem mass spectrometry to reliably and unambiguously establish the identity of each protein. This technology paves the way for studies of cell-type dependent gene expression and studies of large suites of cellular proteins with unprecedented speed and rigor to provide information complementary to the ongoing Human Genome Project.

Detection of competing DNA structures by thermal gradient gel electrophoresis: from self-association to triple helix formation by (G,A)-containing oligonucleotides

Sequence-specific recognition of DNA can be achieved by triple helix-forming oligonucleotides that bind to the major groove of double-helical DNA. These oligonucleotides have been used as sequence-specific DNA ligands for various purposes, including sequence-specific gene regulation in the so-called ‘antigene strategy’. In particular, (G,A)-containing oligonucleotides can form stable triple helices under physiological conditions. However, triplex formation may be in competition with self-association of these oligonucleotides. For biological applications it would be interesting to identify the conditions under which one structure is favoured as compared to the other(s). Here we have directly studied competition between formation of a parallel (G,A) homoduplex and that of a triple helix by a 13 nt (G,A)-containing oligonucleotide. Temperature gradient gel electrophoresis allows simultaneous detection of competition between the two structures, because of their different temperature dependencies and gel electrophoretic mobilities, and characterisation of this competition.

Membrane-bound state of the colicin E1 channel domain as an extended two-dimensional helical array

Atomic level structures have been determined for the soluble forms of several colicins and toxins, but the structural changes that occur after membrane binding have not been well characterized. Changes occurring in the transition from the soluble to membrane-bound state of the C-terminal 190-residue channel polypeptide of colicin E1 (P190) bound to anionic membranes are described. In the membrane-bound state, the α-helical content increases from 60–64% to 80–90%, with a concomitant increase in the average length of the helical segments from 12 to 16 or 17 residues, close to the length required to span the membrane bilayer in the open channel state. The average distance between helical segments is increased and interhelix interactions are weakened, as shown by a major loss of tertiary structure interactions, decreased efficiency of fluorescence resonance energy transfer from an energy donor on helix V of P190 to an acceptor on helix IX, and decreased resonance energy transfer at higher temperatures, not observed in soluble P190, implying freedom of motion of helical segments. Weaker interactions are also shown by a calorimetric thermal transition of low cooperativity, and the extended nature of the helical array is shown by a 3- to 4-fold increase in the average area subtended per molecule to 4,200 Å2 on the membrane surface. The latter, with analysis of the heat capacity changes, implies the absence of a developed hydrophobic core in the membrane-bound P190. The membrane interfacial layer thus serves to promote formation of a highly helical extended two-dimensional flexible net. The properties of the membrane-bound state of the colicin channel domain (i.e., hydrophobic anchor, lengthened and loosely coupled α-helices, and close association with the membrane interfacial layer) are plausible structural features for the state that is a prerequisite for voltage gating, formation of transmembrane helices, and channel opening.

Two-dimensional freezing in the liquid–vapor interface of a dilute Pb:Ga alloy

We report the results of x-ray reflectivity and grazing incidence x-ray diffraction studies of the liquid–vapor interface of a dilute alloy of Pb in Ga over the temperature range of 23–76°C. Our data show that the liquid–vapor interface of this alloy is stratified for several atomic diameters into the bulk liquid and that a monolayer of Pb forms the outermost stratum of the interface. Over the temperature range of 23–56°C, the monolayer of Pb is in an ordered hexagonal phase. At about 58°C, this monolayer undergoes a first-order transition to a hexatic phase, which remains stable to 76°C. An analogy between the observed transition and the first-order melting transition in a one-component classical plasma is suggested.

Two-dimensional IR correlation spectroscopy: Sequential events in the unfolding process of the λ Cro-V55C repressor protein

A question often posed in protein folding/unfolding studies is whether the process is fully cooperative or whether it contains sequential elements. To address this question, one needs tools capable of resolving different events. It seems that, at least in certain cases, two-dimensional (2D) IR correlation spectroscopy can provide answers to this question. To illustrate this point, we have turned to the Cro-V55C dimer of the λ Cro repressor, a protein known to undergo thermal unfolding in two discrete steps through a stable equilibrium intermediate. The secondary structure of this intermediate is compatible with that of a partially unfolded protein and involves a reorganization of the N terminus, whereas the antiparallel β-ribbon formed by the C-terminal part of each subunit remains largely intact. To establish whether the unfolding process involves sequential events, we have performed a 2D correlation analysis of IR spectra recorded over the temperature range of 20–95°C. The 2D IR correlation analysis indeed provides evidence for a sequential formation of the stable intermediate, which is created in three (closely related) steps. A first step entails the unfolding of the short N-terminal β-strand, followed by the unfolding of the α-helices in a second step, and the third step comprises the reorganization of the remaining β-sheet and of some unordered segments in the protein. The complete unfolding of the stable intermediate at higher temperatures also undergoes sequential events that ultimately end with the breaking of the H bonds between the two β-strands at the dimer interface.

Systematic derivation of partition functions for ligand binding to two-dimensional lattices

The Ising problem consists in finding the analytical solution of the partition function of a lattice once the interaction geometry among its elements is specified. No general analytical solution is available for this problem, except for the one-dimensional case. Using site-specific thermodynamics, it is shown that the partition function for ligand binding to a two-dimensional lattice can be obtained from those of one-dimensional lattices with known solution. The complexity of the lattice is reduced recursively by application of a contact transformation that involves a relatively small number of steps. The transformation implemented in a computer code solves the partition function of the lattice by operating on the connectivity matrix of the graph associated with it. This provides a powerful new approach to the Ising problem, and enables a systematic analysis of two-dimensional lattices that model many biologically relevant phenomena. Application of this approach to finite two-dimensional lattices with positive cooperativity indicates that the binding capacity per site diverges as Na (N = number of sites in the lattice) and experiences a phase-transition-like discontinuity in the thermodynamic limit N → ∞. The zeroes of the partition function tend to distribute on a slightly distorted unit circle in complex plane and approach the positive real axis already for a 5×5 square lattice. When the lattice has negative cooperativity, its properties mimic those of a system composed of two classes of independent sites with the apparent population of low-affinity binding sites increasing with the size of the lattice, thereby accounting for a phenomenon encountered in many ligand-receptor interactions.

The two-dimensional IR nonlinear spectroscopy of a cyclic penta-peptide in relation to its three-dimensional structure

A form of two-dimensional (2D) vibrational spectroscopy, which uses two ultrafast IR laser pulses, is used to examine the structure of a cyclic penta-peptide in solution. Spectrally resolved cross peaks occur in the off-diagonal region of the 2D IR spectrum of the amide I region, analogous to those in 2D NMR spectroscopy. These cross peaks measure the coupling between the different amide groups in the structure. Their intensities and polarizations relate directly to the three-dimensional structure of the peptide. With the help of a model coupling Hamiltonian, supplemented by density functional calculations, the spectra of this penta-peptide can be regenerated from the known solution phase structure. This 2D-IR measurement, with an intrinsic time resolution of less than 1 ps, could be used in all time regimes of interest in biology.