Significant discrepancies between van't Hoff and calorimetric enthalpies.

In this paper we show that the usual assumption in studies of the temperature variation of equilibrium constants for equilibria of the form A+B <-->AB that a plot of ln K vs. 1/T (K = equilibrium constant, T = temperature in degrees kelvin) is a straight line with slope equal to -delta HvH/R (delta HvH = van't Hoff or apparent enthalpy, R = gas constant) is not valid in many cases. In all the cases considered here, delta HvH is temperature dependent and is significantly different from the true or calorimetrically measured enthalpy, and the respective values for delta Cp are also significantly different.

Alanine is helix-stabilizing in both template-nucleated and standard peptide helices

Alanine-based peptides of defined sequence and length show measurable helix contents, allowing them to be used as a model system both for analyzing the mechanism of helix formation and for investigating the contributions of side-chain interactions to protein stability. Extensive characterization of many peptide sequences with varying amino acid contents indicates that the favorable helicity of alanine-based peptides can be attributed to the large helix-stabilizing propensity of alanine. Based on their analysis of alanine-rich sequences N-terminally linked to a synthetic helix-inducing template, Kemp and coworkers [Kemp, D. S., Boyd, J. G. & Muendel, C. C. (1991) Nature (London) 352, 451–454; Kemp, D. S., Oslick, S. L. & Allen, T. J. (1996) J. Am. Chem. Soc. 118, 4249–4255] argue that alanine is helix-indifferent, however, and that the favorable helix contents of alanine-based peptides must have some other explanation. Here, we show that the helix contents of template-nucleated sequences are influenced strongly by properties of the template–helix junction. A model in which the helix propensities of residues at the template–peptide junction are treated separately brings the results from alanine-based peptides and template-nucleated helices into agreement. The resulting model provides a physically plausible resolution of the discrepancies between the two systems and allows the helix contents of both template-nucleated and standard peptide helices to be predicted by using a single set of helix propensities. Helix formation in both standard peptides and template–peptide conjugates can be attributed to the large intrinsic helix-forming tendency of alanine.

Lack of cell cycle regulation of telomerase activity in human cells

Conflicting reports have appeared concerning the cell cycle regulation of telomerase activity and its possible repression during quiescence and cell differentiation. We have reexamined these issues in an attempt to uncover the basis for the discrepancies. Variations in extracted telomerase activity during the cell cycle are not observed in cells sorted on the basis of DNA content. Variations are observed in cells synchronized using some biochemical cell cycle inhibitors, but only with those agents where cellular toxicity is evident. A progressive decline in telomerase activity is observed in cells whose growth rate is reduced from seven to eight population doublings per week to one to two doublings per week. Telomerase is largely absent in cells that truly exit the cell cycle and do not divide over the 7-day period. Although it is not necessary for all cell types to regulate telomerase in the same way, we conclude that in the immortal cultured cell lines examined, extracted telomerase activity does not change significantly during progression through the stages of the cell cycle. Telomerase activity generally correlates with growth rate and is repressed in cells that exit the cell cycle and become quiescent.

Clustering of meiotic double-strand breaks on yeast chromosome III

In the yeast Saccharomyces cerevisiae, meiotic recombination is initiated by transient DNA double-strand breaks (DSBs) that are repaired by interaction of the broken chromosome with its homologue. To identify a large number of DSB sites and gain insight into the control of DSB formation at both the local and the whole chromosomal levels, we have determined at high resolution the distribution of meiotic DSBs along the 340 kb of chromosome III. We have found 76 DSB regions, mostly located in intergenic promoter-containing intervals. The frequency of DSBs varies at least 50-fold from one region to another. The global distribution of DSB regions along chromosome III is nonrandom, defining large (39–105 kb) chromosomal domains, both hot and cold. The distribution of these localized DSBs indicates that they are likely to initiate most crossovers along chromosome III, but some discrepancies remain to be explained.

Computation, prediction, and experimental tests of fitness for bacteriophage T7 mutants with permuted genomes

We created a simulation based on experimental data from bacteriophage T7 that computes the developmental cycle of the wild-type phage and also of mutants that have an altered genome order. We used the simulation to compute the fitness of more than 105 mutants. We tested these computations by constructing and experimentally characterizing T7 mutants in which we repositioned gene 1, coding for T7 RNA polymerase. Computed protein synthesis rates for ectopic gene 1 strains were in moderate agreement with observed rates. Computed phage-doubling rates were close to observations for two of four strains, but significantly overestimated those of the other two. Computations indicate that the genome organization of wild-type T7 is nearly optimal for growth: only 2.8% of random genome permutations were computed to grow faster, the highest 31% faster, than wild type. Specific discrepancies between computations and observations suggest that a better understanding of the translation efficiency of individual mRNAs and the functions of qualitatively “nonessential” genes will be needed to improve the T7 simulation. In silico representations of biological systems can serve to assess and advance our understanding of the underlying biology. Iteration between computation, prediction, and observation should increase the rate at which biological hypotheses are formulated and tested.

Contribution of S opponent cells to color appearance

We measured the regions in isoluminant color space over which observers perceive red, yellow, green, and blue and examined the extent to which the colors vary in perceived amount within these regions. We compared color scaling of various isoluminant stimuli by using large spots, which activate all cone types, to that with tiny spots in the central foveola, where S cones, and thus S opponent (So) cell activity, are largely or entirely absent. The addition of So input to that from the L and M opponent cells changes the chromatic appearance of all colors, affecting each primary color in different chromatic regions in the directions and by the amount predicted by our color model. Shifts from white to the various chromatic stimuli we used produced sinusoidal variations in cone activation as a function of color angle for each cone type and in the responses of lateral geniculate cells. However, psychophysical color-scaling functions for 2° spots were nonsinusoidal, being much more peaked. The color-scaling functions are well fit by sine waves raised to exponents between 1 and 3. The same is true for the color responses of a large subpopulation of striate cortex cells. The narrow color tuning, the discrepancies between the spectral loci of the peaks of the color-scaling curves and those of lateral geniculate cells, and the changes in color appearance produced by eliminating So input provide evidence for a cortical processing stage at which the color axes are rotated by a combination of the outputs of So cells with those of L and M opponent cells in the manner that we postulated earlier. There seems to be an expansive response nonlinearity at this stage.

NMR reveals hydrogen bonds between oxygen and distal histidines in oxyhemoglobin

Compared with free heme, the proteins hemoglobin (Hb) and myoglobin (Mb) exhibit greatly enhanced affinity for oxygen relative to carbon monoxide. This physiologically vital property has been attributed to either steric hindrance of CO or stabilization of O2 binding by a hydrogen bond with the distal histidine. We report here the first direct evidence of such a hydrogen bond in both α- and β-chains of oxyhemoglobin, as revealed by heteronuclear NMR spectra of chain-selectively labeled samples. Using these spectra, we have assigned the imidazole ring 1H and 15N chemical shifts of the proximal and distal histidines in both carbonmonoxy- and oxy-Hb. Because of their proximity to the heme, these chemical shifts are extremely sensitive to the heme pocket conformation. Comparison of the measured chemical shifts with values predicted from x-ray structures suggests differences between the solution and crystal structures of oxy-Hb. The chemical shift discrepancies could be accounted for by very small displacements of the proximal and distal histidines. This suggests that NMR could be used to obtain very high-resolution heme pocket structures of Hb, Mb, and other heme proteins.

Cosmic velocity fields and their interpretation

We review the current status of our knowledge of cosmic velocity fields, on both small and large scales. A new statistic is described that characterizes the incoherent, thermal component of the velocity field on scales less than 2h−1 Mpc (h is H0/100 km·s−1·Mpc−1, where H0 is the Hubble constant and 1 Mpc = 3.09 × 1022 m) and smaller. The derived velocity is found to be quite stable across different catalogs and is of remarkably low amplitude, consistent with an effective Ω ∼ 0.15 on this scale. We advocate the use of this statistic as a standard diagnostic of the small-scale kinetic energy of the galaxy distribution. The analysis of large-scale flows probes the velocity field on scales of 10–60 h−1 Mpc and should be adequately described by linear perturbation theory. Recent work has focused on the comparison of gravity or density fields derived from whole-sky redshift surveys of galaxies [e.g., the Infrared Astronomical Satellite (IRAS)] with velocity fields derived from a variety of sources. All the algorithms that directly compare the gravity and velocity fields suggest low values of the density parameter, while the POTENT analysis, using the same data but comparing the derived IRAS galaxy density field with the Mark-III derived matter density field, leads to much higher estimates of the inferred density. Since the IRAS and Mark-III fields are not fully consistent with each other, the present discrepancies might result from the very different weighting applied to the data in the competing methods.

Gene recognition via spliced sequence alignment.

Gene recognition is one of the most important problems in computational molecular biology. Previous attempts to solve this problem were based on statistics, and applications of combinatorial methods for gene recognition were almost unexplored. Recent advances in large-scale cDNA sequencing open a way toward a new approach to gene recognition that uses previously sequenced genes as a clue for recognition of newly sequenced genes. This paper describes a spliced alignment algorithm and software tool that explores all possible exon assemblies in polynomial time and finds the multiexon structure with the best fit to a related protein. Unlike other existing methods, the algorithm successfully recognizes genes even in the case of short exons or exons with unusual codon usage; we also report correct assemblies for genes with more than 10 exons. On a test sample of human genes with known mammalian relatives, the average correlation between the predicted and actual proteins was 99%. The algorithm correctly reconstructed 87% of genes and the rare discrepancies between the predicted and real exon-intron structures were caused either by short (less than 5 amino acids) initial/terminal exons or by alternative splicing. Moreover, the algorithm predicts human genes reasonably well when the homologous protein is nonvertebrate or even prokaryotic. The surprisingly good performance of the method was confirmed by extensive simulations: in particular, with target proteins at 160 accepted point mutations (PAM) (25% similarity), the correlation between the predicted and actual genes was still as high as 95%.

The effects of sequence context on DNA curvature.

Recent experiments have exposed significant discrepancies between experimental data and predictive models for DNA structure. These results strongly suggest that DNA structural parameters incorporated in the models are not always sufficient to account for the influence of sequence context and of specific ion effects. In an attempt to evaluate these two effects, we have investigated repetitive DNA sequences with the sequence motif GAGAG.CTCTC located in different helical phasing arrangements with respect to poly(A) tracts and GGGCCC.GGGCCC sequence motifs. Methods used are ligase-mediated cyclization and gel mobility experiments along with DNase I cutting and chemical probe studies. The results provide new evidence for curvature in poly(A) tracts. They also show that the sequence context in which bending and flexible sequence elements are found is an important aspect of sequence-dependent DNA conformation. Although dinucleotide models generally have good predictive power, this work demonstrates that in some instances sequence elements larger than the dinucleotide must be taken into account, and hence it provides a starting point for the appropriate modification and refinement of existing structural models for DNA.

Measurement of the binding of tyrosyl phosphopeptides to SH2 domains: a reappraisal.

Src homology 2 (SH2) domain-mediated interactions with phosphotyrosine residues are critical in many intracellular signal transduction pathways. Attempts to understand the determinants of specificity and selectivity of these interactions have prompted many binding studies that have used several techniques. Some discrepancies, in both the absolute and relative values of the dissociation constants for particular interactions, are apparent. To establish the correct dissociation constants and to understand the origin of these differences, we have analyzed three previously determined interactions using the techniques of surface plasmon resonance and isothermal titration calorimetry. We find that the binding of SH2 domains to phosphopeptides is weaker than generally presumed. A phosphopeptide based on the hamster polyoma middle tumor antigen interacts with the SH2 domain from Src with an equilibrium dissociation constant (Kd) of 600 nM; a phosphopeptide based on one binding site from the platelet-derived growth factor receptor binds to the N-terminal SH2 domain of the 1-phosphatidylinositol 3-kinase p85 subunit with a Kd of 300 nM; and a phosphopeptide based on the C terminus of Lck binds to the SH2 domain of Lck with a Kd of 4 microM. In addition, we demonstrate that avidity effects that result from the dimerization of glutathione S-transferase fusion proteins with SH2 domains could be responsible for overestimates of affinities for these interactions previously studied by surface plasmon resonance.