Flow cytometric analysis of human lysozyme production in recombinantSaccharomyces cerevisiae

Flow cytometric techniques were used to investigate cell size, protein content and cell cycle behavior of recombinantSaccharomyces cerevisiae strains producing human lysozyme (HLZ). Two different signal sequences, the native yeastMFα1 signal sequence and the rat α-amylase signal sequence, were used for secretion of HLZ. The strain containing the rat α-amylase signal sequence showed a higher level of internal lysozyme and lower specific growth rates. Flow cytometric analysis of the total protein content and cell size showed the strain harboring the native yeast signal sequence had a higher total protein content than the strain containing the rat α-amylase signal sequence. Cell cycle analysis indicated that the two lysozyme producing recombinant strains had an increased number of cells in the G₂+M phase of the yeast cell cycle compared with the host strain SEY2102.     

A stopped-flow fluorescence study of the native and modified lysozyme

The protein folding kinetics of hen egg white lysozyme (HEWL) was studied using experimental and bioinformatics tools. The structure of the transition state in the unfolding pathway of lysozyme was determined with stopped-flow kinetics using intact HEWL and its chemically modified derivative, in which six lysine residues have been modified. The overall consistency of φ-value (φ ≈ 1) indicates that lysine side chains interactions are subject to breaking in the structure of the transition state. Following experimental evidences, multiple sequence alignment of lysozyme family in vertebrates and exact structural examination of lysozyme, showed that the α-helix in the structure of lysozyme has critical role in the unfolding kinetics.     

Influence of assignment on the prediction of transmembrane helices in protein structures

α-Helical transmembrane proteins (TMP_α) are composed of a series of helices embedded in the lipid bilayer. Due to technical difficulties, few 3D structures are available. Therefore, the design of structural models of TMP_α is of major interest. We study the secondary structures of TMP_α by analyzing the influence of secondary structures assignment methods (SSAMs). For this purpose, a published and updated benchmark databank of TMP_α is used and several SSAMs (9) are evaluated. The analysis of the results points to significant differences in SSA depending on the methods used. Pairwise comparisons between SSAMs led to more than 10% of disagreement. Helical regions corresponding to transmembrane zones are often correctly characterized. The study of the sequence–structure relationship shows very limited differences with regard to the structural disagreement. Secondary structure prediction based on Bayes’ rule and using only a single sequence give correct prediction rates ranging from 78 to 81%. A structural alphabet approach gives a slightly better prediction, i.e., only 2% less than the best equivalent approach, whereas the prediction rate with a very different assignment bypasses 86%. This last result highlights the importance of the correct assignment choice to evaluate the prediction assessment.     

Genes for two homologous G-protein α subunits map to different human chromosomes

Summary Signal transduction across biological membranes is modulated by a family of related GTP-binding proteins termed G proteins. These G proteins have a heterotrimeric structure composed of α, β, and γ subunits. The α subunits of the G proteins bind GTP and appear to determine the biochemical specificity of the protein. We have recently cloned and characterized cDNA encoding two G-protein α subunits, α_i and α_h. The former is a substrate for ADP-ribosylation by pertussis toxin. The protein corresponding to α_h has not yet been identified. These cDNAs encode proteins, which demonstrate 90% sequence identity to one another and also show marked similarity to other G proteins. The present studies were designed to determine whether the genes for these related proteins are clustered on a single human chromosome. Genomic DNA isolated from a panel of mouse-human hybrid cell lines was analyzed by hybridization to cDNAs for α_i and α_h. Based on the distribution patterns of α_i and α_h in cell hybrids, the gene for α_i was assigned to human chromosome 7, and the gene for α_h assigned to chromosome 12. These data suggest that the G-protein gene family may be distributed over at least two human chromosomes.     

Sequential nearest-neighbor effects on computed <Superscript>13</Superscript>C<Superscript>α</Superscript> chemical shifts

To evaluate sequential nearest-neighbor effects on quantum-chemical calculations of ¹³C^α chemical shifts, we selected the structure of the nucleic acid binding (NAB) protein from the SARS coronavirus determined by NMR in solution (PDB id 2K87). NAB is a 116-residue α/β protein, which contains 9 prolines and has 50% of its residues located in loops and turns. Overall, the results presented here show that sizeable nearest-neighbor effects are seen only for residues preceding proline, where Pro introduces an overestimation, on average, of 1.73 ppm in the computed ¹³C^α chemical shifts. A new ensemble of 20 conformers representing the NMR structure of the NAB, which was calculated with an input containing backbone torsion angle constraints derived from the theoretical ¹³C^α chemical shifts as supplementary data to the NOE distance constraints, exhibits very similar topology and comparable agreement with the NOE constraints as the published NMR structure. However, the two structures differ in the patterns of differences between observed and computed ¹³C^α chemical shifts, Δ _ca,i , for the individual residues along the sequence. This indicates that the Δ _ca,i -values for the NAB protein are primarily a consequence of the limited sampling by the bundles of 20 conformers used, as in common practice, to represent the two NMR structures, rather than of local flaws in the structures.     

Solution NMR structures reveal a distinct architecture and provide first structures for protein domain family PF04536

The protein family (Pfam) PF04536 is a broadly conserved domain family of unknown function (DUF477), with more than 1,350 members in prokaryotic and eukaryotic proteins. High-quality NMR structures of the N-terminal domain comprising residues 41–180 of the 684-residue protein CG2496 from Corynebacterium glutamicum and the N-terminal domain comprising residues 35–182 of the 435-residue protein PG0361 from Porphyromonas gingivalis both exhibit an α/β fold comprised of a four-stranded β-sheet, three α-helices packed against one side of the sheet, and a fourth α-helix attached to the other side. In spite of low sequence similarity (18%) assessed by structure-based sequence alignment, the two structures are globally quite similar. However, moderate structural differences are observed for the relative orientation of two of the four helices. Comparison with known protein structures reveals that the α/β architecture of CG2496(41–180) and PG0361(35–182) has previously not been characterized. Moreover, calculation of surface charge potential and identification of surface clefts indicate that the two domains very likely have different functions.     

Orangutan α-satellite monomers are closely related to the human consensus sequence

Alpha-satellite is a family of tandemly repeated DNA found at the centromeric regions of all human and primate chromosomes. Human α-satellite subsets are largely chromosome-specific and have been further grouped into four suprachromosomal families (SFs), each characterized by a unique set of monomeric types. Although chimpanzee and gorilla α-satellites share sufficient sequence similarity to fit the established SFs, the assumption that the derived human α-satellite consensus and monomeric types represent the sequence of ancestral repeats remains unestablished. By using oligonucleotide primers specific for a conserved region of human α-satellite DNA, we have PCR amplified, cloned, and characterized α-satellite sequences from the orangutan genome. Nucleotide sequence analysis demonstrated that orangutan α-satellite is formed by a single monomeric type that is significantly closer in percentage of sequence identity (mean = 92%, range = 89–96%) to the overall consensus of human α-satellite than to the monomeric types corresponding to the four SFs. Use of cloned sequences as hybridization probes to orangutan genomic DNA digested with a panel of restriction enzymes showed that most orangutan α-satellite subsets are characterized by a monomeric construction. The subset homologous to clone PPY2-5 is organized in distinct higher-order repeat structures consisting of 18 adjacent monomers. By FISH two clones, PPY3-4 and PPY3-5, proved to be specific for the α-satellite on the orangutan homologs of human Chromosomes (Chrs) 10 and 8, respectively. Our data indicate that there was an ancestral monomeric type displaying high sequence similarity to the overall human consensus from which the different great ape and human subsets and SFs may have originated. Received: 24 November 1997 / Accepted: 29 January 1998     

How uniform is the peptide plane geometry? A high-accuracy NMR study of dipolar C<Superscript>α</Superscript>–C′/H<Superscript>N</Superscript>–N cross-correlated relaxation

Highly precise and accurate measurements of very small NMR cross-correlated relaxation rates, namely those between protein H_i^N–N_i and C_i−1^α–C_i−1′ dipoles, are demonstrated with an error of 0.03 s⁻¹ for GB3. Because the projection angles between the two dipole vectors are very close to the magic angle the rates range only from −0.2 to +0.2 s⁻¹. Small changes of the average vector orientations have a dramatic impact on the relative values. The rates suggest deviation from idealized peptide plane geometry caused by twists around the C′–N bonds and/or pyramidalization of the nitrogen atoms. A clear alternating pattern along the sequence is observed in β strands 1, 3 and 4 of GB3, where the side chains of almost all residues with large positive rates are solvent exposed. In the α helix all rates are relatively large and positive. Some of the currently most accurate structures of GB3 determined by both high resolution X-ray crystallography and NMR are in satisfactory agreement with the experimental rates in the helix and β strand 3, but not in the loops and the two central strands of the sheet for which no alternating pattern is predicted.     

Success of the test task performance by students with different spectral α-rhythm characteristics in the baseline electroencephalogram

The subjects were divided into two groups according to the α₁ and α₂ spectral powers in the occipital derivations of the EEG recorded in the initial state with their eyes closed. Group I included subjects whose α₁-rhythm spectral power (7–10 Hz) was more than 70% of the total α-rhythm band power. Group II included subjects whose α₂-rhythm spectral power (10–13 Hz) was more than 70% of the total α-rhythm band power. It was established that, in the tasks requiring prediction of the subsequent result (memorizing a certain sequence of signals and its subsequent reproduction on the monitor screen), group I subjects differed from group II subjects in fewer sequence errors and a greater number of accurate predictions. In group II subjects, a decrease in the ϑ-band spectral power in the EEG of the central and frontal cortical areas was observed compared to the baseline. Therefore, the EEG ϑ-rhythm power at the memorizing stage was lower in them than in group II subjects. The results suggest that the baseline characteristics of ϑ-activity can be regarded as prognostic criteria of similar types of activity.     

Methods of NMR structure refinement: molecular dynamics simulations improve the agreement with measured NMR data of a C-terminal peptide of GCN4-p1

The C-terminal trigger sequence is essential in the coiled-coil formation of GCN4-p1; its conformational properties are thus of importance for understanding this process at the atomic level. A solution NMR model structure of a peptide, GCN4p16–31, encompassing the GCN4-p1 trigger sequence was proposed a few years ago. Derived using a standard single-structure refinement protocol based on 172 nuclear Overhauser effect (NOE) distance restraints, 14 hydrogen-bond and 11 ϕ torsional-angle restraints, the resulting set of 20 NMR model structures exhibits regular α-helical structure. However, the set slightly violates some measured NOE bounds and does not reproduce all 15 measured ³J(H_N-H_Cα)-coupling constants, indicating that different conformers of GCN4p16–31 might be present in solution. With the aim to resolve structures compatible with all NOE upper distance bounds and ³J-coupling constants, we executed several structure refinement protocols employing unrestrained and restrained molecular dynamics (MD) simulations with two force fields. We find that only configurational ensembles obtained by applying simultaneously time-averaged NOE distance and ³J-coupling constant restraining with either force field reproduce all the experimental data. Additionally, analyses of the simulated ensembles show that the conformational variability of GCN4p16–31 in solution admitted by the available set of 187 measured NMR data is larger than represented by the set of the NMR model structures. The conformations of GCN4p16–31 in solution differ in the orientation not only of the side-chains but also of the backbone. The inconsistencies between the NMR model structures and the measured NMR data are due to the neglect of averaging effects and the inclusion of hydrogen-bond and torsional-angle restraints that have little basis in the primary, i.e. measured NMR data.     

Distinct interactions between the human adrenergic β2 receptor and Gαs—an in silico study

The aim of this study was to perform an in silico analysis of the interaction of the human β₂ adrenergic receptor with Gα_s. In a first step, a systematic surface-interaction-scan between the inactive or active human β₂ adrenergic receptor and Gα_s was performed in order to gain knowledge about energetically preferred areas on the potential energy surface. Subsequently, two energetically favored regions for the active human β₂ adrenergic receptor–Gα_s complex were identified. Two representative complex structures were put into a POPC (1-palmitoyl-2-oleoyl-phosphatidylcholine) bilayer and solvated in order to perform molecular dynamic simulations. The simulations revealed that both conformations, which have comparable potential energy, are stable. A mean number of about 14 hydrogen bonds was observed between the active receptor and Gα_s for both conformations. Based on these results, two energetically favored β₂–Gα_scomplexes can be proposed.     

Domain Evolution in the α-Amylase Family

The available amino acid sequences of the α-amylase family (glycosyl hydrolase family 13) were searched to identify their domain B, a distinct domain that protrudes from the regular catalytic (β/α)₈-barrel between the strand β3 and the helix α3. The isolated domain B sequences were inspected visually and also analyzed by Hydrophobic Cluster Analysis (HCA) to find common features. Sequence analyses and inspection of the few available three-dimensional structures suggest that the secondary structure of domain B varies with the enzyme specificity. Domain B in these different forms, however, may still have evolved from a common ancestor. The largest number of different specificities was found in the group with structural similarity to domain B from Bacillus cereus oligo-1,6-glucosidase that contains an α-helix succeeded by a three-stranded antiparallel β-sheet. These enzymes are α-glucosidase, cyclomaltodextrinase, dextran glucosidase, trehalose-6-phosphate hydrolase, neopullulanase, and a few α-amylases. Domain B of this type was observed also in some mammalian proteins involved in the transport of amino acids. These proteins show remarkable similarity with (β/α)₈-barrel elements throughout the entire sequence of enzymes from the oligo-1,6-glucosidase group. The transport proteins, in turn, resemble the animal 4F2 heavy-chain cell surface antigens, for which the sequences either lack domain B or contain only parts thereof. The similarities are compiled to indicate a possible route of domain evolution in the α-amylase family. Received: 4 December 1996 / Accepted: 13 March 1997     

The α-gliadin gene family. II. DNA and protein sequence variation, subfamily structure, and origins of pseudogenes

 The derived amino-acid sequences of all reported α-gliadin clones are compared and analyzed, and the patterns of sequence change within the α-gliadin family are examined. The most variable sequences are two polyglutamine domains. These two domains are characteristic features of the α-gliadin storage proteins and account for most of the variation in protein size of this otherwise highly conserved protein family. In addition, their encoding DNA sequences form microsatellites. Single-base substitutions in the α-gliadin genes show a preponderance of transitions, including the C to T substitution which contributes to the generation of stop codons, and consequently to the observation that approximately 50% of the α-gliadin genes are pseudogenes. In one unusual gene, a microsatellite has expanded to 321 bp as compared to the normal 36–72 bp, and may result from similar mechanisms that produce polyglutamine-associated genetic diseases in humans. A comparison of the 27 reported sequences show several α-gliadin gene subfamilies, at least some of which are genome specific. Received: 1 October 1996/Accepted: 20 December 1996     

Close Evolutionary Relatedness of α-Amylases from Archaea and Plants

The amino acid sequences of 22 α-amylases from family 13 of glycosyl hydrolases were analyzed with the aim of revealing the evolutionary relationships between the archaeal α-amylases and their eubacterial and eukaryotic counterparts. Two evolutionary distance trees were constructed: (i) the first one based on the alignment of extracted best-conserved sequence regions (58 residues) comprising β2, β3, β4, β5, β7, and β8 strand segments of the catalytic (α/β)₈-barrel and a short conserved stretch in domain B protruding out of the barrel in the β3 →α3 loop, and (ii) the second one based on the alignment of the substantial continuous part of the (α/β)₈-barrel involving the entire domain B (consensus length: 386 residues). With regard to archaeal α-amylases, both trees compared brought, in fact, the same results; i.e., all family 13 α-amylases from domain Archaea were clustered with barley pI isozymes, which represent all plant α-amylases. The enzymes from Bacillus licheniformis and Escherichia coli, representing liquefying and cytoplasmic α-amylases, respectively, seem to be the further closest relatives to archaeal α-amylases. This evolutionary relatedness clearly reflects the discussed similarities in the amino acid sequences of these α-amylases, especially in the best-conserved sequence regions. Since the results for α-amylases belonging to all three domains (Eucarya, Eubacteria, Archaea) offered by both evolutionary trees are very similar, it is proposed that the investigated conserved sequence regions may indeed constitute the ``sequence fingerprints' of a given α-amylase. Received: 3 June 1998 / Accepted: 20 August 1998     

Preferred structures of constrained peptides from achiral α,α-dialkyiated glycyl residues with acyclic side chains

Conformational energy computations of the monopeptides from three achiral α,α-dialkylated glycyl residues with acyclic side chains (namely α,α-dimethyl-; α,α-diethyl-; and α, α-di-n-propylglycines) are reported as a function of the relevant N-C^α-C′ bond angle. In parallel, experimental studies were performed in the solid state (infrared absorption and X-ray diffraction) and in solution (infrared absorption and proton magnetic resonance) on the corresponding protected homo-peptide series (the former series to the dodecamer, the other two series to the pentamers). The results obtained unequivocally indicate that the preference from a helical to a fully extended conformation increases as side-chain bulkiness increases. The longest homo-peptides from α,α-dimethylglycine form stable 3₁₀-helices. A picture of the mode of self-association of the helical structures has also been determined. The results of the theoretical analyses fit well with the experimentally observed conformational properties in the solid state and in chloroform solution.     

Genomic organization of the rat α2u-globulin gene cluster

The α_2u-globulins are a group of similar proteins, belonging to the lipocalin superfamily of proteins, that are synthesized in a subset of secretory tissues in rats. The many α_2u-globulin isoforms are encoded by a multigene family that exhibits extensive homology. Despite a high degree of sequence identity, individual family members show diverse expression patterns involving complex hormonal, tissue-specific, and developmental regulation. Analysis suggests that there are approximately 20 α_2u-globulin genes in the rat genome. We have used fluorescence in situ hybridization (FISH) to show that the α_2u-globulin genes are clustered at a single site on rat Chromosome (Chr) 5 (5q22-24). Southern blots of rat genomic DNA separated by pulsed field gel electrophoresis indicated that the α_2u-globulin genes are contained on two NruI fragments with a total size of 880 kbp. Analysis of three P1 clones containing α_2u-globulin genes indicated that the α_2u-globulin genes are tandemly arranged in a head-to-tail fashion. The organization of the α_2u-globulin genes in the rat as a tandem array of single genes differs from the homologous major urinary protein genes in the mouse, which are organized as tandem arrays of divergently oriented gene pairs. The structure of these gene clusters may have consequences for the proposed function, as a pheromone transporter, for the protein products encoded by these genes. Received: 12 August 1998 / Accepted: 13 January 1999     

Conformational Studies of the 313-320 and 313-332 Peptide Fragments Derived from the αIIb Subunit of Integrin Receptor with Molecular Dynamics Simulations

The peptide sequence YMESRADRKLAEVGRVYLFL, derived from 313-332 region of the α_IIb, has been identified as a potent inhibitor of platelet aggregation and fibrinogen binding to α_IIbβ₃. More detailed studies have revealed that the Y₃₁₃MESRADR₃₂₀ sequence is the shortest octapeptide with strong inhibitory activity. This work provides insight of the solution conformation of these peptides, by performing extensive molecular dynamics simulations of 100 ns. The 8mer peptide has no stable conformation in water while the 20mer peptide retains a relative conformational stability. Analysis of side chain orientation of the RAD fragment revealed the synplanar arrangement of guanidinium and β-carboxylic groups providing a framework for explaining the similar biological activity of the two peptides, despite their differences in sequence and conformation.     

Detecting Natural Selection at the Molecular Level: A Reexamination of Some “Classic” Examples of Adaptive Evolution

An important criterion used to detect adaptive evolution in DNA sequence data is ω_i > 1, where ω_i is the ratio of nonsynonymous to synonymous substitution rates in lineage i. However, the evaluation of multiple ω_i within a phylogenetic tree can easily inflate the statistical type I error rate. We developed two rigorous methods of analysis that avoid this and other potential pitfalls. We applied these methods to four published examples of adaptive evolution. One case was strongly supported by our reanalysis (abalone sperm lysin), and one was weakly supported (baboon α-globin), but two examples (primate lysozyme and Antarctic fish β-globin) did not show significant evidence of adaptive evolution. Our first method is a “bottom-up” hierarchical maximum likelihood approach, which (1) tests for significant heterogeneity in ω across the phylogeny, (2) locates its source using a sequence of planned comparisons, and (3) tests homogeneous groups of ω for ω > 1, using a modified level of significance that incorporates the pretesting. The second method is a “top-down” log-linear analysis based on estimates of nonsynonymous and synonymous substitutions in pairs of lineages. The log-linear test is applied to pairs of lineages joined at progressively deeper nodes. For each pair, the analysis simultaneously tests for adaptive evolution (ω > 1), a shift in natural selection (ω₁ ≠ω₂), and unequal evolution rate (the relative rate test). In both tests, we emphasized that the criterion ω₁ ≠ ω₂ is an important additional indicator of a phylogenetic shift in the balance between natural selection and genetic drift between two related lineages. [Reviewing Editor: Dr. John Huelsenbeck]     

Isolation and characterization of a new phycoerythrin from the cyanobacterium <Emphasis Type="Italic">Synechococcus</Emphasis> sp. ECS-18

A new phycoerythrin, SCH-phycoerythrin, was purified from Synechococcus sp. ECS-18 by DEAE-Sephacel anion exchange chromatography and Sephacryl S-300 gel filtration. The protein pigment had an absorbance maximum at 542 nm and a fluorescence maximum at 565 nm. The native molecular mass was approximately 219 kDa as determined by gel filtration, and sodium dodecyl sulfate polyacrylamide gel electrophoresis demonstrated the presence of two subunits, with molecular mass of 19 and 17.9 kDa. These observations are consistent with the (αβ)₆ subunit composition that is characteristic of phycoerythrins. The α- and β-subunits showed immunological identity by Ouchterlony double immunodiffusion with an anti-phycoerythrin antiserum. The DNA sequence of the SCH-phycoerythrin gene was determined by PCR amplification using primers based on the conserved N-terminal amino acid sequence of the α- and β-subunits of phycoerythrins.