Molecular dynamics of mesophilic-like mutants of a cold-adapted enzyme: insights into distal effects induced by the mutations

Networks and clusters of intramolecular interactions, as well as their "communication" across the three-dimensional architecture have a prominent role in determining protein stability and function. Special attention has been dedicated to their role in thermal adaptation. In the present contribution, seven previously experimentally characterized mutants of a cold-adapted α-amylase, featuring mesophilic-like behavior, have been investigated by multiple molecular dynamics simulations, essential dynamics and analyses of correlated motions and electrostatic interactions. Our data elucidate the molecular mechanisms underlying the ability of single and multiple mutations to globally modulate dynamic properties of the cold-adapted α-amylase, including both local and complex unpredictable distal effects. Our investigation also shows, in agreement with the experimental data, that the conversion of the cold-adapted enzyme in a warm-adapted variant cannot be completely achieved by the introduction of few mutations, also providing the rationale behind these effects. Moreover, pivotal residues, which are likely to mediate the effects induced by the mutations, have been identified from our analyses, as well as a group of suitable candidates for protein engineering. In fact, a subset of residues here identified (as an isoleucine, or networks of mesophilic-like salt bridges in the proximity of the catalytic site) should be considered, in experimental studies, to get a more efficient modification of the features of the cold-adapted enzyme.     

Role of lysine versus arginine in enzyme cold-adaptation: modifying lysine to homo-arginine stabilizes the cold-adapted alpha-amylase from Pseudoalteramonas haloplanktis

The cold-adapted alpha-amylase from Pseudoalteromonas haloplanktis (AHA) is a multidomain enzyme capable of reversible unfolding. Cold-adapted proteins, including AHA, have been predicted to be structurally flexible and conformationally unstable as a consequence of a high lysine-to-arginine ratio. In order to examine the role of low arginine content in structural flexibility of AHA, the amino groups of lysine were guanidinated to form homo-arginine (hR), and the structure-function-stability properties of the modified enzyme were analyzed by transverse urea gradient-gel electrophoresis. The extent of modification was monitored by MALDI-TOF-MS, and correlated to changes in activity and stability. Modifying lysine to hR produced a conformationally more stable and less active alpha-amylase. The k(cat) of the modified enzyme decreased with a concomitant increase in deltaH# and decrease in K(m). To interpret the structural basis of the kinetic and thermodynamic properties, the hR residues were modeled in the AHA X-ray structure and compared to the X-ray structure of a thermostable homolog. The experimental properties of the modified AHA were consistent with K106hR forming an intra-Domain B salt bridge to stabilize the active site and decrease the cooperativity of unfolding. Homo-Arg modification also appeared to alter Ca2+ and Cl- binding in the active site. Our results indicate that replacing lysine with hR generates mesophilic-like characteristics in AHA, and provides support for the importance of lysine residues in promoting enzyme cold adaptation. These data were consistent with computational analyses that show that AHA possesses a compositional bias that favors decreased conformational stability and increased flexibility.     

Temperature adaptation of proteins: engineering mesophilic-like activity and stability in a cold-adapted alpha-amylase

Two multiple mutants of a psychrophilic alpha-amylase were produced, bearing five mutations (each introducing additional weak interactions found in pig pancreatic alpha-amylase) with or without an extra disulfide bond specific to warm-blooded animals. Both multiple mutants display large modifications of stability and activity arising from synergic effects in comparison with single mutations. Newly introduced weak interactions and the disulfide bond confer mesophilic-like stability parameters, as shown by increases in the melting point t(m), in the calorimetric enthalpy DeltaH(cal) and in protection against heat inactivation, as well as by decreases in cooperativity and reversibility of unfolding. In addition, both kinetic and thermodynamic activation parameters of the catalyzed reaction are shifted close to the values of the porcine enzyme. This study confirms the central role of weak interactions in regulating the balance between stability and activity of an enzyme in order to adapt to the environmental temperature.     

Adaptation of a thermophilic enzyme, 3-isopropylmalate dehydrogenase, to low temperatures

Random mutagenesis coupled with screening of the active enzyme at a low temperature was applied to isolate cold-adapted mutants of a thermophilic enzyme. Four mutant enzymes with enhanced specific activities (up to 4.1-fold at 40 degrees C) at a moderate temperature were isolated from randomly mutated Thermus thermophilus 3-isopropylmalate dehydrogenase. Kinetic analysis revealed two types of cold-adapted mutants, i.e. k(cat)-improved and K(m)-improved types. The k(cat)-improved mutants showed less temperature-dependent catalytic properties, resulting in improvement of k(cat) (up to 7.5-fold at 40 degrees C) at lower temperatures with increased K(m) values mainly for NAD. The K(m)-improved enzyme showed higher affinities toward the substrate and the coenzyme without significant change in k(cat) at the temperatures investigated (30-70 degrees C). In k(cat)-improved mutants, replacement of a residue was found near the binding pocket for the adenine portion of NAD. Two of the mutants retained thermal stability indistinguishable from the wild-type enzyme. Extreme thermal stability of the thermophilic enzyme is not necessarily decreased to improve the catalytic function at lower temperatures. The present strategy provides a powerful tool for obtaining active mutant enzymes at lower temperatures. The results also indicate that it is possible to obtain cold-adapted mutant enzymes with high thermal stability.     

Optimisation of the surface electrostatics as a strategy for cold adaptation of uracil-DNA N-glycosylase (UNG) from Atlantic cod (Gadus morhua)

Cold-adapted enzymes are characterised by an increased catalytic efficiency and reduced temperature stability compared to their mesophilic counterparts. Lately, it has been suggested that an optimisation of the electrostatic surface potential is a strategy for cold adaptation for some enzymes. A visualisation of the electrostatic surface potential of cold-adapted uracil-DNA N-glycosylase (cUNG) from Atlantic cod indicates a more positively charged surface near the active site compared to human UNG (hUNG). In order to investigate the importance of the altered surface potential for the cold-adapted features of cod UNG, six mutants have been characterised and compared to cUNG and hUNG. The cUNG quadruple mutant (V171E, K185V, H250Q and H275Y) and four corresponding single mutants all comprise substitutions of residues present in the human enzyme. A human UNG mutant, E171V, comprises the equivalent residue found in cod UNG. In addition, crystal structures of the single mutants V171E and E171V have been determined. Results from the study show that a more negative electrostatic surface potential reduces the activity and increases the stability of cod UNG, and suggest an optimisation of the surface potential as a strategy for cold-adaptation of this enzyme. Val171 in cod UNG is especially important in this respect.     

2-Amino-N6-hydroxyadenine induces gene/point mutations and multiple-locus mutations, but not multilocus deletion mutations, in the ad-3 region of a two-component heterokaryon of Neurospora crassa.

The mutagenicity of 2-amino-N6-hydroxyadenine (AHA) has been studied in Neurospora crassa by treating a two-component heterokaryon (H-12) and recovering specific-locus mutations induced in the ad-3 region. This assay system permits the identification of ad-3A and/or ad-3B mutants resulting from gene/point mutations, multilocus deletion mutations, and multiple-locus mutations of various genotypes, involving one or both loci. Genetic characterization of the ad-3 mutants recovered from experiments with AHA in H-12 shows that 98.9% (270/273) of the ad-3 mutants are gene/point mutations (ad-3R), 1.1% (3/270) are unknowns, and none is a multilocus deletion mutation (ad-3IR). Among the gene/point mutations, 3.3% (9/273) are multiple-locus mutations (gene/point mutations with a closely-linked recessive lethal mutation [ad-3R + RLCL]). Another 25.3% (69/273) are multiple-locus mutations with a recessive lethal mutation located elsewhere in the genome [ad-3R + RL]. Heterokaryon tests for allelic complementation among the ad-3BR mutants showed that 90.8% (139/153) of the mutants were complementing, and 20.3% (31/153) were leaky. In addition, 32.5% (38/117) of the ad-3AR mutants were leaky. These data are consistent with the hypothesis that AHA produces specific-locus mutations in the ad-3 region of N. crassa by base-pair substitution. The data from the present experiments are compared with the data for 2-aminopurine (2AP)-induced ad-3 mutants in H-12 (de Serres and Brockman, 1991). Whereas, 2AP is a weak mutagen in H-12, AHA is extremely potent (Brockman et al., 1987). In contrast with 2AP, AHA induces ad-3 mutants exclusively by gene/point mutation in H-12. We conclude that whereas AHA induces ad-3 mutants predominantly by AT to GC base-pair transitions, 2AP induces ad-3 mutants by a wide variety of mechanisms including: (1) AT to GC and GC to AT base-pair transitions, (2) frameshift mutations, (3) other, as yet unidentified, intragenic alterations, (4) small multilocus deletion mutations, and (5) multiple-locus ad-3R mutations with closely linked recessive lethal mutations.     

Comparison of multiple molecular dynamics trajectories calculated for the drug-resistant HIV-1 integrase T66I/M154I catalytic domain

HIV-1 integrase (IN) is an essential enzyme for the viral replication and an interesting target for the design of new pharmaceuticals for multidrug therapy of AIDS. Single and multiple mutations of IN at residues T66, S153, or M154 confer degrees of resistance to several inhibitors that prevent the enzyme from performing its normal strand transfer activity. Four different conformations of IN were chosen from a prior molecular dynamics (MD) simulation on the modeled IN T66I/M154I catalytic core domain as starting points for additional MD studies. The aim of this article is to understand the dynamic features that may play roles in the catalytic activity of the double mutant enzyme in the absence of any inhibitor. Moreover, we want to verify the influence of using different starting points on the MD trajectories and associated dynamical properties. By comparison of the trajectories obtained from these MD simulations we have demonstrated that the starting point does not affect the conformational space explored by this protein and that the time of the simulation is long enough to achieve convergence for this system.     

In silico screening,molecular docking,and molecular dynamics studies of SNP-derived human P5CR mutants

Pyrroline-5-carboxylate reductase (P5CR) encoded by PYCR1 gene is a housekeeping enzyme that catalyzes the reduction of P5C to proline using NAD(P)H as the cofactor. In this study, we used in silico approaches to examine the role of nonsynonymous single-nucleotide polymorphisms in the PYCR1 gene and their putative functions in the pathogenesis of Cutis Laxa. Among the 348 identified SNPs, 15 were predicted to be potentially damaging by both SIFT and PolyPhen tools; of them two SNP‐derived mutations, R119G and G206W, have been previously reported to correlate with Cutis Laxa. These two mutations were therefore selected to be mapped to the wild‐type (WT) P5CR structure for further structural and functional analyses. The results of comparative computational analyses using I-Mutant and Autodock reveal reductions in both stability and cofactor binding affinity of these two mutants. Comparative molecular dynamics (MD) simulations were performed to evaluate the changes in dynamic properties of P5CR upon mutations. The results reveal that the two mutations enhance the rigidity of P5CR structure, especially that of cofactor binding site, which could result in decreased kinetics of cofactor entrance and egress. Comparison between the structural properties of the WT and mutants during MD simulations shows that the enhanced rigidity of mutants results most likely from the increased number of inter‐atomic interactions and the decreased number of dynamic hydrogen bonds. Our study provides novel insight into the deleterious effects of the R119G and G206W mutations on P5CR, and sheds light on the mechanisms by which these mutations mediate Cutis Laxa.     

Genotoxicity of 2-amino-6-N-hydroxyadenine (AHA) to mouse lymphoma and CHO cells.

2-Amino-6-N-hydroxyadenine (AHA) treated L5178Y/TK (+/-)-3.7.2C mouse lymphoma cells were evaluated for mutations at the tk, hgprt, and Na+/K+ ATPase loci, as well as for gross chromosome aberrations and induction of micronuclei. In addition, AHA was evaluated for its ability to induce HGPRT mutants in CHO cells. AHA was found to induce mutations at all evaluated loci and in both cell types. The TK mutants were primarily large colonies although a few small colonies were also induced, particularly at the higher concentrations. Preliminary cytogenetic analysis of AHA-treated mouse lymphoma cells indicated that some gross aberrations but not micronuclei were induced. The 20 small-colony TK mutants evaluated by banded karyotype indicate that only a small fraction (2 of 20) showed chromosome 11 abnormalities. From these studies, it appears that AHA may be one of a very few chemicals that is capable of inducing multi-locus point mutations, with only slight clastogenic activity. Particularly at the higher concentrations, some of the mutants may contain multi-locus point mutations that result in slow growth.     

Evaluation of Influence of Single Nucleotide Polymorphisms in Cytochrome P450 2B6 on Substrate Recognition Using Computational Docking and Molecular Dynamics Simulation

In this study, we investigated the influence of single nucleotide polymorphisms on the conformation of mutated cytochrome P450 (CYP) 2B6 proteins using molecular dynamics (MD) simulation. Some of these mutations influence drug metabolism activities, leading to individual variations in drug efficacy and pharmacokinetics. Using computational docking, we predicted the structure of the complex between the antimalarial agent artemether and CYP2B6 whose conformations were obtained by MD simulation. The simulation demonstrated that the entire structure of the protein changes even when a single residue is mutated. Moreover, the structural flexibility is affected by the mutations and it may influence the enzyme activity. The results suggest that some of the inactive mutants cannot recognize artemether due to structural changes caused by the mutation.     

The active site is the least stable structure in the unfolding pathway of a multidomain cold-adapted alpha-amylase

The cold-active alpha-amylase from the Antarctic bacterium Pseudoalteromonas haloplanktis (AHA) is the largest known multidomain enzyme that displays reversible thermal unfolding (around 30 degrees C) according to a two-state mechanism. Transverse urea gradient gel electrophoresis (TUG-GE) from 0 to 6.64 M was performed under various conditions of temperature (3 degrees C to 70 degrees C) and pH (7.5 to 10.4) in the absence or presence of Ca2+ and/or Tris (competitive inhibitor) to identify possible low-stability domains. Contrary to previous observations by strict thermal unfolding, two transitions were found at low temperature (12 degrees C). Within the duration of the TUG-GE, the structures undergoing the first transition showed slow interconversions between different conformations. By comparing the properties of the native enzyme and the N12R mutant, the active site was shown to be part of the least stable structure in the enzyme. The stability data supported a model of cooperative unfolding of structures forming the active site and independent unfolding of the other more stable protein domains. In light of these findings for AHA, it will be valuable to determine if active-site instability is a general feature of heat-labile enzymes from psychrophiles. Interestingly, the enzyme was also found to refold and rapidly regain activity after being heated at 70 degrees C for 1 h in 6.5 M urea. The study has identified fundamental new properties of AHA and extended our understanding of structure/stability relationships of cold-adapted enzymes.     

Human mitochondrial pyrophosphatase: cDNA cloning and analysis of the gene in patients with mtDNA depletion syndromes

Pyrophosphatases (PPases) catalyze the hydrolysis of inorganic pyrophosphate generated in several cellular enzymatic reactions. A novel human pyrophosphatase cDNA encoding a 334-amino-acid protein approximately 60% identical to the previously identified human cytosolic PPase was cloned and characterized. The novel enzyme, named PPase-2, was enzymatically active and catalyzed hydrolysis of pyrophosphate at a rate similar to that of the previously identified PPase-1. A functional mitochondrial import signal sequence was identified in the N-terminus of PPase-2, which targeted the enzyme to the mitochondrial matrix. The human pyrophosphatase 2 gene (PPase-2) was mapped to chromosome 4q25 and the 1.4-kb mRNA was ubiquitously expressed in human tissues, with highest levels in muscle, liver, and kidney. The yeast homologue of the mitochondrial PPase-2 is required for mitochondrial DNA maintenance and yeast cells lacking the enzyme exhibit mitochondrial DNA depletion. We sequenced the PPA2 gene in 13 patients with mitochondrial DNA depletion syndromes (MDS) of unknown cause to determine if mutations in the PPA2 gene of these patients were associated with this disease. No pathogenic mutations were identified in the PPA2 gene of these patients and we found no evidence that PPA2 gene mutations are a common cause of MDS in humans.     

Mutagenicity of 2-amino-N6-hydroxyadenine to TK6 human lymphoblast cells

TK6 human lymphoblast cells (tk +/-; hprt+) were treated with various concentrations of 2-amino-N6-hydroxyadenine (AHA) for 24 h. AHA was quite toxic to TK6 cells in the dose range 0-0.05 micrograms/ml, but additional toxicity was not observed between 0.05 and 0.10 micrograms/ml. AHA induced mutations at 2 distinct genetic loci: the autosomal thymidine kinase (tk) and the X-linked hypoxanthine-guanine phosphoribosyl transferase (hprt). Significant levels of both tk-NG mutants (normal growth rate of 16-18 h, colonies visible after 10-11 days incubation) and tk-SG mutants (slow growth rate of greater than 24 h, colonies visible after 18 days incubation) were induced. 15 hprt- mutants were isolated and analyzed by Southern blot. 8 of these had normal restriction fragment patterns after digestion with PstI, EcoRI, and HindIII, and were defined as 'point' mutations; the remaining 7 had partial deletions of the hprt gene. 32 tk- mutants were also isolated. 3 of 22 normal growth mutants and 6 of 10 slow growth mutants had lost the active tk allele. These data suggest that both point mutations and larger-scale alterations are induced by AHA.     

DROIDS 3.0—Detecting Genetic and Drug Class Variant Impact on Conserved Protein Binding Dynamics

The application of statistical methods to comparatively framed questions about the molecular dynamics (MD) of proteins can potentially enable investigations of biomolecular function beyond the current sequence and structural methods in bioinformatics. However, the chaotic behavior in single MD trajectories requires statistical inference that is derived from large ensembles of simulations representing the comparative functional states of a protein under investigation. Meaningful interpretation of such complex forms of big data poses serious challenges to users of MD. Here, we announce Detecting Relative Outlier Impacts from Molecular Dynamic Simulation (DROIDS) 3.0, a method and software package for comparative protein dynamics that includes maxDemon 1.0, a multimethod machine learning application that trains on large ensemble comparisons of concerted protein motions in opposing functional states generated by DROIDS and deploys learned classifications of these states onto newly generated MD simulations. Local canonical correlations in learning patterns generated from independent, yet identically prepared, MD validation runs are used to identify regions of functionally conserved protein dynamics. The subsequent impacts of genetic and/or drug class variants on conserved dynamics can also be analyzed by deploying the classifiers on variant MD simulations and quantifying how often these altered protein systems display opposing functional states. Here, we present several case studies of complex changes in functional protein dynamics caused by temperature, genetic mutation, and binding interactions with nucleic acids and small molecules. We demonstrate that our machine learning algorithm can properly identify regions of functionally conserved dynamics in ubiquitin and TATA-binding protein (TBP). We quantify the impact of genetic variation in TBP and drug class variation targeting the ATP-binding region of Hsp90 on conserved dynamics. We identify regions of conserved dynamics in Hsp90 that connect the ATP binding pocket to other functional regions. We also demonstrate that dynamic impacts of various Hsp90 inhibitors rank accordingly with how closely they mimic natural ATP binding.     

Cold adaptation of a psychrophilic chitinase: a mutagenesis study

The gene encoding chitinase ArChiB from the Antarctic Arthrobacter sp. TAD20 has been expressed in Escherichia coli and the recombinant enzyme purified to homogeneity. In an effort to engineer cold-adapted biocatalysts through rational redesign to operate at elevated temperatures, we performed several mutations aiming to increase the rigidity of the molecular edifice of the selected psychrophilic chitinase. The mutations were designed on the basis of a homology-based three-dimensional model of the enzyme, and included an attempt to introduce a salt bridge (mutant N198K) and replacements of selected Gly residues by either Pro (mutants G93P, G254P) or Gln (G406Q). Mutant N198K resulted in a more stable protein (DeltaTm = 0.6 degrees C). Mutant G93P exhibited a DeltaTm of 1.2 degrees C, while mutants G254P and G406Q exhibited decreased stability. We conclude that the effect of mutating Gly residues on enzyme stability is rather complex and can only be understood in the context of the structural environment. Kinetic and spectroscopic analysis of these enzyme variants revealed that the kinetic parameters kcat and Km have been significantly modified.     

A Bifunctional Enzyme in Pseudomonas aeruginosa: a New Pattern in the Organization of Enzymes Concerned with Phenylalanine and Tyrosine Biosynthesis

Two isozymes of chorismate mutase (CA mutase(1) and CA mutase(2)) and two isozymes of prephenate dehydratase (PPA dehydratase(1) and PPA dehydratase(2)) have been found in Pseudomonas aeruginosa. The activities CA mutase(2)-PPA dehydratase(2) catalyzing phenylalanine biosynthesis have been purified almost 40-fold and were found to be associated as a bifunctional enzyme or an enzyme complex. The enzymes specific for tyrosine biosynthesis did not appear to manifest such physical association. Thus, the organization of enzymes concerned with phenylalanine and tyrosine biosynthesis in P. aeruginosa is unique and is unlike most other organisms. Single site mutants have been isolated which have lost both CA mutase(2)-PPA dehydratase(2) activities resulting in a requirement for phenylalanine for growth. Single site revertants of these mutants regained both these activities simultaneously and were able to grow on minimal medium. A mutant, r(6), was also isolated which had normal CA mutase(2) but lacked PPA dehydratase(2) activity.     

Molecular dynamics simulation of a carboxy murine neuroglobin mutated on the proximal side: heme displacement and concomitant rearrangement in loop regions

Neuroglobin, a member of vertebrate globin family, is distributed primarily in the brain and retina. Considerable evidence has accumulated regarding its unique ligand-binding properties, neural-specific distribution, distinct expression regulation, and possible roles in processes such as neuron protection and enzymatic metabolism. Structurally, neuroglobin enjoys unique features, such as bis-histidyl coordination to heme iron in the absence of exogenous ligand, heme orientational heterogeneity, and a heme sliding mechanism accompanying ligand binding. In the present work, molecular dynamics (MD) simulations were employed to reveal functional and structural information in three carboxyl murine neuroglobin mutants with single point mutations F106Y, F106L and F106I, respectively. The MD simulation indicates a remarkable proximal effect on detectable displacement of heme and a larger tunnel in the protein matrix. In addition, the mutation at F106 confers on the CD region a very sensitive mobility in all three model structures. The dynamic features of neuroglobin demonstrate rearrangement of the inner space and highly active loop regions in solution. These imply that the conserved residue at the G5 site plays a key role in the physiological function of this unusual protein.     

A DNA ligase from the psychrophile Pseudoalteromonas haloplanktis gives insights into the adaptation of proteins to low temperatures.

The cloning, overexpression and characterization of a cold-adapted DNA ligase from the Antarctic sea water bacterium Pseudoalteromonas haloplanktis are described. Protein sequence analysis revealed that the cold-adapted Ph DNA ligase shows a significant level of sequence similarity to other NAD+-dependent DNA ligases and contains several previously described sequence motifs. Also, a decreased level of arginine and proline residues in Ph DNA ligase could be involved in the cold-adaptation strategy. Moreover, 3D modelling of the N-terminal domain of Ph DNA ligase clearly indicates that this domain is destabilized compared with its thermophilic homologue. The recombinant Ph DNA ligase was overexpressed in Escherichia coli and purified to homogeneity. Mass spectroscopy experiments indicated that the purified enzyme is mainly in an adenylated form with a molecular mass of 74 593 Da. Ph DNA ligase shows similar overall catalytic properties to other NAD+-dependent DNA ligases but is a cold-adapted enzyme as its catalytic efficiency (kcat/Km) at low and moderate temperatures is higher than that of its mesophilic counterpart E. coli DNA ligase. A kinetic comparison of three enzymes adapted to different temperatures (P. haloplanktis, E. coli and Thermus scotoductus DNA ligases) indicated that an increased kcat is the most important adaptive parameter for enzymatic activity at low temperatures, whereas a decreased Km for the nicked DNA substrate seems to allow T. scotoductus DNA ligase to work efficiently at high temperatures. Besides being useful for investigation of the adaptation of enzymes to extreme temperatures, P. haloplanktis DNA ligase, which is very efficient at low temperatures, offers a novel tool for biotechnology.     

Mutagenicity of 2-amino-N6-hydroxyadenine (AHA) at three loci in L5178Y/tk+/- mouse lymphoma cells: molecular and preliminary cytogenetic characterizations of AHA-induced tk-/- mutants.

2-Amino-N6-hydroxyadenine (AHA) is a remarkably efficient and specific inducer of point mutations in Neurospora, with few or no larger scale events being detected (de Serres et al., 1985). In the present studies, AHA is shown to be a potent point mutagen at the tk +/-, hprt+ and Na+/K+ ATPase loci in L5178Y/tk (+/-)-3.7.2C mouse lymphoma cells. Both large and small colony tk-/- mutants were analyzed at the molecular level and a preliminary assessment was made of small colony mutant karyotypes (230 bands/haploid metaphase cell; large colony mutants typically have normal karyotypes and were not analyzed). AHA induced greatly delayed (7-9 cell doublings) cytotoxicity, suggestive of a mutational mechanism (e.g., base-pair substitution) requiring DNA replication prior to its phenotypic expression. Approximately one-third of the tk -/- mutants formed small colonies, a phenotype which is typically associated with alterations to chromosome 11b, the site of the functional tkb allele in the parental cells. However, banded karyotypes have provided convincing evidence for alterations chromosome 11b in only 2 of the 7 small colony mutants analyzed. Southern blot analysis showed that 78% (21/27) of these small colony mutants have retained the Nco-1 6.3-kb band, which is diagnostic of the tkb allele. This makes AHA unique among the mutagens examined so far in inducing small colony mutants without inducing large losses of tkb DNA. Although a dose-dependent increase in the proportion of small colony mutants was noted, no significant dose-dependent differences were seen at the molecular level in the relatively few mutants analyzed. The majority of AHA-induced tk -/- mutants formed large colonies. Southern blot analysis showed that 86% (25/29) of these had retained the Nco-1 6.3-kb band which is diagnostic of the tkb allele. It is concluded that AHA induces primarily micromutations (less than 100 base pairs), probably through a base-pair substitution mechanism, at the tk, hprt and Na+/K+ ATPase loci in this system, with some larger scale damage (kilobases of DNA at the molecular level; chromosome 11b damage at the cytogenetic level) also occurring.