A subset of skin tumor modifier loci determines survival time of tumor-bearing mice

Studies of mouse models of human cancer have established the existence of multiple tumor modifiers that influence parameters of cancer susceptibility such as tumor multiplicity, tumor size, or the probability of malignant progression. We have carried out an analysis of skin tumor susceptibility in interspecific Mus musculus/Mus spretus hybrid mice and have identified another seven loci showing either significant (six loci) or suggestive (one locus) linkage to tumor susceptibility or resistance. A specific search was carried out for skin tumor modifier loci associated with time of survival after development of a malignant tumor. A combination of resistance alleles at three markers [D6Mit15 (Skts12), D7Mit12 (Skts2), and D17Mit7 (Skts10)], all of which are close to or the same as loci associated with carcinoma incidence and/or papilloma multiplicity, is significantly associated with increased survival of mice with carcinomas, whereas the reverse combination of susceptibility alleles is significantly linked to early mortality caused by rapid carcinoma growth (χ2 = 25.22; P = 5.1 × 10−8). These data indicate that host genetic factors may be used to predict carcinoma growth rate and/or survival of individual backcross mice exposed to the same carcinogenic stimulus and suggest that mouse models may provide an approach to the identification of genetic modifiers of cancer survival in humans.

A new Eocene archaeocete (Mammalia, Cetacea) from India and the time of origin of whales

Himalayacetus subathuensis is a new pakicetid archaeocete from the Subathu Formation of northern India. The type dentary has a small mandibular canal indicating a lack of auditory specializations seen in more advanced cetaceans, and it has Pakicetus-like molar teeth suggesting that it fed on fish. Himalayacetus is significant because it is the oldest archaeocete known and because it was found in marine strata associated with a marine fauna. Himalayacetus extends the fossil record of whales about 3.5 million years back in geological time, to the middle part of the early Eocene [≈53.5 million years ago (Ma)]. Oxygen in the tooth-enamel phosphate has an isotopic composition intermediate between values reported for freshwater and marine archaeocetes, indicating that Himalayacetus probably spent some time in both environments. When the temporal range of Archaeoceti is calibrated radiometrically, comparison of likelihoods constrains the time of origin of Archaeoceti and hence Cetacea to about 54–55 Ma (beginning of the Eocene), whereas their divergence from extant Artiodactyla may have been as early as 64–65 Ma (beginning of the Cenozoic).

Planetary exploration in the time of astrobiology: Protecting against biological contamination

These are intriguing times in the exploration of other solar-system bodies. Continuing discoveries about life on Earth and the return of data suggesting the presence of liquid water environments on or under the surfaces of other planets and moons have combined to suggest the significant possibility that extraterrestrial life may exist in this solar system. Similarly, not since the Viking missions of the mid-1970s has there been as great an appreciation for the potential for Earth life to contaminate other worlds. Current plans for the exploration of the solar system include constraints intended to prevent biological contamination from being spread by solar-system exploration missions.

Probing the oligomeric structure of an enzyme by electrospray ionization time-of-flight mass spectrometry.

Electrospray ionization time-of-flight (ESI-TOF) mass spectrometry was used to study the quaternary structure of 4-oxalocrotonate tautomerase (EC 5.3.2; 4OT), and four analogues prepared by total chemical synthesis. Wild-type 4OT is a hexamer of 62 amino acid subunits and contains no cysteine residues. The analogues were: (desPro1)4OT, a truncated construct in which Pro1 was deleted; (Cpc1)4OT in which Pro1 was replaced with cyclopentane carboxylate; a derivative [Met(O)45]4OT in which Met45 was oxidized to the sulfoxide; and an analogue (Nle45)4OT in which Met45 was replaced with norleucine. ESI of (Nle45)4OT, (Cpc1)4OT, and 4OT from solution conditions under which the native enzyme was fully active (5 mM ammonium bicarbonate buffer, pH 7.5) gave the intact hexamer as the major species detected by TOF mass spectrometry. In contrast, analysis of [Met(O)45]4OT and (desPro1)4OT under similar conditions yielded predominantly monomer ions. The ESI-TOF measurements were consistent with structural data obtained from circular dichroism spectroscopy. In the context of kinetic data collected for 4OT and these analogues, ESI-TOF mass spectrometry also provided important evidence for the structural and mechanistic significance of the catalytically important Pro1 residue in 4OT.

DNA sequencing by delayed extraction-matrix-assisted laser desorption/ionization time of flight mass spectrometry.

Matrix-assisted laser desorption/ionization (MALDI) time of flight mass spectrometry was used to detect and order DNA fragments generated by Sanger dideoxy cycle sequencing. This was accomplished by improving the sensitivity and resolution of the MALDI method using a delayed ion extraction technique (DE-MALDI). The cycle sequencing chemistry was optimized to produce as much as 100 fmol of each specific dideoxy terminated fragment, generated from extension of a 13-base primer annealed on 40- and 50-base templates. Analysis of the resultant sequencing mixture by DE-MALDI identified the appropriate termination products. The technique provides a new non-gel-based method to sequence DNA which may ultimately have considerable speed advantages over traditional methodologies.

Brownian dynamics simulations of protein folding: Access to milliseconds time scale and beyond

Protein folding occurs on a time scale ranging from milliseconds to minutes for a majority of proteins. Computer simulation of protein folding, from a random configuration to the native structure, is nontrivial owing to the large disparity between the simulation and folding time scales. As an effort to overcome this limitation, simple models with idealized protein subdomains, e.g., the diffusion–collision model of Karplus and Weaver, have gained some popularity. We present here new results for the folding of a four-helix bundle within the framework of the diffusion–collision model. Even with such simplifying assumptions, a direct application of standard Brownian dynamics methods would consume 10,000 processor-years on current supercomputers. We circumvent this difficulty by invoking a special Brownian dynamics simulation. The method features the calculation of the mean passage time of an event from the flux overpopulation method and the sampling of events that lead to productive collisions even if their probability is extremely small (because of large free-energy barriers that separate them from the higher probability events). Using these developments, we demonstrate that a coarse-grained model of the four-helix bundle can be simulated in several days on current supercomputers. Furthermore, such simulations yield folding times that are in the range of time scales observed in experiments.

Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

Coalescent estimates of HIV-1 generation time in vivo

The generation time of HIV Type 1 (HIV-1) in vivo has previously been estimated using a mathematical model of viral dynamics and was found to be on the order of one to two days per generation. Here, we describe a new method based on coalescence theory that allows the estimate of generation times to be derived by using nucleotide sequence data and a reconstructed genealogy of sequences obtained over time. The method is applied to sequences obtained from a long-term nonprogressing individual at five sampling occasions. The estimate of viral generation time using the coalescent method is 1.2 days per generation and is close to that obtained by mathematical modeling (1.8 days per generation), thus strengthening confidence in estimates of a short viral generation time. Apart from the estimation of relevant parameters relating to viral dynamics, coalescent modeling also allows us to simulate the evolutionary behavior of samples of sequences obtained over time.

Stopped-flow NMR spectroscopy: real-time unfolding studies of 6-19F-tryptophan-labeled Escherichia coli dihydrofolate reductase.

Escherichia coli dihydrofolate reductase (DHFR; EC 1.5.1.3) contains five tryptophan residues that have been replaced with 6-19F-tryptophan. The 19F NMR assignments are known in the native, unliganded form and the unfolded form. We have used these assignments with stopped-flow 19F NMR spectroscopy to investigate the behavior of specific regions of the protein in real time during urea-induced unfolding. The NMR data show that within 1.5 sec most of the intensities of the native 19F resonances of the protein are lost but only a fraction (approximately 20%) of the intensities of the unfolded resonances appears. We postulate that the early disappearance of the native resonances indicates that most of the protein rapidly forms an intermediate in which the side chains have considerable mobility. Stopped-flow far-UV circular dichroism measurements indicate that this intermediate retains native-like secondary structure. Eighty percent of the intensities of the NMR resonances assigned to the individual tryptophans in the unfolded state appear with similar rate constants (k approximately 0.14 sec-1), consistent with the major phase of unfolding observed by stopped-flow circular dichroism (representing 80% of total amplitude). These data imply that after formation of the intermediate, which appears to represent an expanded structural form, all regions of the protein unfold at the same rate. Stopped-flow measurements of the fluorescence and circular dichroism changes associated with the urea-induced unfolding show a fast phase (half-time of about 1 sec) representing 20% of the total amplitude in addition to the slow phase mentioned above. The NMR data show that approximately 20% of the total intensity for each of the unfolded tryptophan resonances is present at 1.5 sec, indicating that these two phases may represent the complete unfolding of the two different populations of the native protein.

Modeling and optimization of populations subject to time-dependent mutation.

It has become clear that many organisms possess the ability to regulate their mutation rate in response to environmental conditions. So the question of finding an optimal mutation rate must be replaced by that of finding an optimal mutation schedule. We show that this task cannot be accomplished with standard population-dynamic models. We then develop a "hybrid" model for populations experiencing time-dependent mutation that treats population growth as deterministic but the time of first appearance of new variants as stochastic. We show that the hybrid model agrees well with a Monte Carlo simulation. From this model, we derive a deterministic approximation, a "threshold" model, that is similar to standard population dynamic models but differs in the initial rate of generation of new mutants. We use these techniques to model antibody affinity maturation by somatic hypermutation. We had previously shown that the optimal mutation schedule for the deterministic threshold model is phasic, with periods of mutation between intervals of mutation-free growth. To establish the validity of this schedule, we now show that the phasic schedule that optimizes the deterministic threshold model significantly improves upon the best constant-rate schedule for the hybrid and Monte Carlo models.