Random mutagenesis of Thermus aquaticus DNA polymerase I: concordance of immutable sites in vivo with the crystal structure.

Expression of Thermus aquaticus (Taq) DNA polymerase I (pol I) in Escherichia, coli complements the growth defect caused by a temperature-sensitive mutation in the host pol I. We replaced the nucleotide sequence encoding amino acids 659-671 of the O-helix of Taq DNA pol I, corresponding to the substrate binding site, with an oligonucleotide containing random nucleotides. Functional Taq pol I mutants were selected based on colony formation at the nonpermissive temperature. By using a library with 9% random substitutions at each of 39 positions, we identified 61 active Taq pol I mutants, each of which contained from one to four amino acid substitutions. Some amino acids, such as alanine-661 and threonine-664, were tolerant of several or even many diverse replacements. In contrast, no replacements or only conservative replacements were identified at arginine-659, lysine-663, and tyrosine-671. By using a library with totally random nucleotides at five different codons (arginine-659, arginine-660, lysine-663, phenylalanine-667, and glycine-668), we confirmed that arginine-659 and lysine-663 were immutable, and observed that only tyrosine substituted for phenylalanine-667. The two immutable residues and the two residues that tolerate only highly conservative replacements lie on the side of O-helix facing the incoming deoxynucleoside triphosphate, as determined by x-ray analysis. Thus, we offer a new approach to assess concordance of the active conformation of an enzyme, as interpreted from the crystal structure, with the active conformation inferred from in vivo function.

Heterodimer formation and activity in the human enzyme galactose-1-phosphate uridylyltransferase.

One of the fundamental questions concerning expression and function of dimeric enzymes involves the impact of naturally occurring mutations on subunit assembly and heterodimer activity. This question is of particular interest for the human enzyme galactose-l-phosphate uridylyl-transferase (GALT), impairment of which results in the inherited metabolic disorder galactosemia, because many if not most patients studied to date are compound heterozygotes rather than true molecular homozygotes. Furthermore, the broad range of phenotypic severity observed in these patients raises the possibility that allelic combination, not just allelic constitution, may play some role in determining outcome. In the work described herein, we have selected two distinct naturally occurring null mutations of GALT, Q188R and R333W, and asked the questions (i) what are the impacts of these mutations on subunit assembly, and (ii) if heterodimers do form, are they active? To answer these questions, we have established a yeast system for the coexpression of epitope-tagged alleles of human GALT and investigated both the extent of specific GALT subunit interactions and the activity of defined heterodimer pools. We have found that both homodimers and heterodimers do form involving each of the mutant subunits tested and that both heterodimer pools retain substantial enzymatic activity. These results are significant not only in terms of their implications for furthering our understanding of galactosemia and GALT holoenzyme structure-function relationships but also because the system described may serve as a model for similar studies of other complexes composed of multiple subunits.

Structure and function in rhodopsin: correct folding and misfolding in point mutants at and in proximity to the site of the retinitis pigmentosa mutation Leu-125-->Arg in the transmembrane helix C.

L125R is a mutation in the transmembrane helix C of rhodopsin that is associated with autosomal dominant retinitis pigmentosa. To probe the orientation of the helix and its packing in the transmembrane domain, we have prepared and studied the mutations E122R, I123R, A124R, S127R, L125F, and L125A at, and in proximity to, the above mutation site. Like L125R, the opsin expressed in COS-1 cells from E122R did not bind 11-cis-retinal, whereas those from I123R and S127R formed the rhodopsin chromophore partially. A124R opsin formed the rhodopsin chromophore (lambda max 495 nm) in the dark, but the metarhodopsin II formed on illumination decayed about 6.5 times faster than that of the wild type and was defective in transducin activation. The mutant opsins from L125F and L125A bound 11-cis-retinal only partially, and in both cases, the mixtures of the proteins produced were separated into retinal-binding and non-retinal-binding (misfolded) fractions. The purified mutant rhodopsin from L125F showed lambda max at 500 nm, whereas that from L125A showed lambda max at 503 nm. The mutant rhodopsin L125F showed abnormal bleaching behavior and both mutants on illumination showed destabilized metarhodopsin II species and reduced transducin activation. Because previous results have indicated that misfolding in rhodopsin is due to the formation of a disulfide bond other than the normal disulfide bond between Cys-110 and Cys-187 in the intradiscal domain, we conclude from the misfolding in mutants L125F and L125A that the folding in vivo in the transmembrane domain is coupled to that in the intradiscal domain.

Structural changes of tumor necrosis factor alpha associated with membrane insertion and channel formation.

Low pH enhances tumor necrosis factor alpha (TNF)-induced cytolysis of cancer cells and TNF-membrane interactions that include binding, insertion, and ion-channel formation. We have also found that TNF increases Na+ influx in cells. Here, we examined the structural features of the TNF-membrane interaction pathway that lead to channel formation. Fluorometric studies link TNF's acid-enhanced membrane interactions to rapid but reversible acquisition of hydrophobic surface properties. Intramembranous photolabeling shows that (i) protonation of TNF promotes membrane insertion, (ii) the physical state of the target bilayer affects the kinetics and efficiency of TNF insertion, and (iii) binding and insertion of TNF are two distinct events. Acidification relaxes the trimeric structure of soluble TNF so that the cryptic carboxyl termini, centrally located at the base of the trimer cone, become susceptible to carboxypeptidase Y. After membrane insertion, TNF exhibits a trimeric configuration in which the carboxyl termini are no longer exposed; however, the proximal salt-bridged Lys-11 residues as well as regional surface amino acids (Glu-23, Arg-32, and Arg-44) are notably more accessible to proteases. The sequenced cleavage products bear the membrane-restricted photoreactive probe, proof that surface-cleaved TNF has an intramembranous disposition. In summary, the trimer's structural plasticity is a major determinant of its channel-forming ability. Channel formation occurs when cracked or partially splayed trimers bind and penetrate the bilayer. Reannealing leads to a slightly relaxed trimeric structure. The directionality of bilayer penetration conforms with x-ray data showing that receptor binding to the monomer interfaces of TNF poises the tip of the trimeric cone directly above the target cell membrane.

Three-dimensional structure of recombinant human osteogenic protein 1: structural paradigm for the transforming growth factor beta superfamily.

We report the three-dimensional structure of osteogenic protein 1 (OP-1, also known as bone morphogenetic protein 7) to 2.8-A resolution. OP-1 is a member of the transforming growth factor beta (TGF-beta) superfamily of proteins and is able to induce new bone formation in vivo. Members of this superfamily share sequence similarity in their C-terminal regions and are implicated in embryonic development and adult tissue repair. Our crystal structure makes possible the structural comparison between two members of the TGF-beta superfamily. We find that although there is limited sequence identity between OP-1 and TGF-beta 2, they share a common polypeptide fold. These results establish a basis for proposing the OP-1/TGF-beta 2 fold as the primary structural motif for the TGF-beta superfamily as a whole. Detailed comparison of the OP-1 and TGF-beta 2 structures has revealed striking differences that provide insights into how these growth factors interact with their receptors.

Bathorhodopsin structure in the room-temperature rhodopsin photosequence: picosecond time-resolved coherent anti-Stokes Raman scattering.

Structural changes in the retinal chromophore during the formation of the bathorhodopsin intermediate (bathoRT) in the room-temperature rhodopsin (RhRT) photosequence (i.e., vision) are examined using picosecond time-resolved coherent anti-Stokes Raman scattering. Specifically, the retinal structure assignable to bathoRT following 8-ps excitation of RhRT is measured via vibrational Raman spectroscopy at a 200-ps time delay where the only intermediate present is bathoRT. Significant differences are observed between the C=C stretching frequencies of the retinal chromophore at low temperature where bathorhodopsin is stabilized and at room temperature where bathorhodopsin is a transient species in the RhRT photosequence. These vibrational data are discussed in terms of the formation of bathoRT, an important step in the energy storage/transduction mechanism of RhRT.

Projection structure of frog rhodopsin in two crystal forms.

Rhodopsin is the G protein-coupled receptor that upon light activation triggers the visual transduction cascade. Rod cell outer segment disc membranes were isolated from dark-adapted frog retinas and were extracted with Tween detergents to obtain two-dimensional rhodopsin crystals for electron crystallography. When Tween 80 was used, tubular structures with a p2 lattice (a = 32 A, b = 83 A, gamma = 91 degrees) were formed. The use of a Tween 80/Tween 20 mixture favored the formation of larger p22(1)2(1) lattices (a = 40 A, b = 146 A, gamma = 90 degrees). Micrographs from frozen hydrated frog rhodopsin crystals were processed, and projection structures to 7-A resolution for the p22(1)2(1) form and to 6-A resolution for the p2 form were calculated. The maps of frog rhodopsin in both crystal forms are very similar to the 9-A map obtained previously for bovine rhodopsin and show that the arrangement of the helices is the same. In a tentative topographic model, helices 4, 6, and 7 are nearly perpendicular to the plane of the membrane. In the higher-resolution projection maps of frog rhodopsin, helix 5 looks more tilted than it appeared previously. The quality of the two frog rhodopsin crystals suggests that they would be suitable to obtain a three-dimensional structure in which all helices would be resolved.

Analysis of homeodomain function by structure-based design of a transcription factor.

The homeodomain is a 60-amino acid module which mediates critical protein-DNA and protein-protein interactions for a large family of regulatory proteins. We have used structure-based design to analyze the ability of the Oct-1 homeodomain to nucleate an enhancer complex. The Oct-1 protein regulates herpes simplex virus (HSV) gene expression by participating in the formation of a multiprotein complex (C1 complex) which regulates alpha (immediate early) genes. We recently described the design of ZFHD1, a chimeric transcription factor containing zinc fingers 1 and 2 of Zif268, a four-residue linker, and the Oct-1 homeodomain. In the presence of alpha-transinduction factor and C1 factor, ZFHD1 efficiently nucleates formation of the C1 complex in vitro and specifically activates gene expression in vivo. The sequence specificity of ZFHD1 recruits C1 complex formation to an enhancer element which is not efficiently recognized by Oct-1. ZFHD1 function depends on the recognition of the Oct-1 homeodomain surface. These results prove that the Oct-1 homeodomain mediates all the protein-protein interactions that are required to efficiently recruit alpha-transinduction factor and C1 factor into a C1 complex. The structure-based design of transcription factors should provide valuable tools for dissecting the interactions of DNA-bound domains in other regulatory circuits.

X-ray structure of clotting factor IXa: active site and module structure related to Xase activity and hemophilia B.

Hereditary deficiency of factor IXa (fIXa), a key enzyme in blood coagulation, causes hemophilia B, a severe X chromosome-linked bleeding disorder afflicting 1 in 30,000 males; clinical studies have identified nearly 500 deleterious variants. The x-ray structure of porcine fIXa described here shows the atomic origins of the disease, while the spatial distribution of mutation sites suggests a structural model for factor X activation by phospholipid-bound fIXa and cofactor VIIIa. The 3.0-A-resolution diffraction data clearly show the structures of the serine proteinase module and the two preceding epidermal growth factor (EGF)-like modules; the N-terminal Gla module is partially disordered. The catalytic module, with covalent inhibitor D-Phe-1I-Pro-2I-Arg-3I chloromethyl ketone, most closely resembles fXa but differs significantly at several positions. Particularly noteworthy is the strained conformation of Glu-388, a residue strictly conserved in known fIXa sequences but conserved as Gly among other trypsin-like serine proteinases. Flexibility apparent in electron density together with modeling studies suggests that this may cause incomplete active site formation, even after zymogen, and hence the low catalytic activity of fIXa. The principal axes of the oblong EGF-like domains define an angle of 110 degrees, stabilized by a strictly conserved and fIX-specific interdomain salt bridge. The disorder of the Gla module, whose hydrophobic helix is apparent in electron density, can be attributed to the absence of calcium in the crystals; we have modeled the Gla module in its calcium form by using prothrombin fragment 1. The arched module arrangement agrees with fluorescence energy transfer experiments. Most hemophilic mutation sites of surface fIX residues occur on the concave surface of the bent molecule and suggest a plausible model for the membrane-bound ternary fIXa-FVIIIa-fX complex structure: fIXa and an equivalently arranged fX arch across an underlying fVIIIa subdomain from opposite sides; the stabilizing fVIIIa interactions force the catalytic modules together, completing fIXa active site formation and catalytic enhancement.

The three-dimensional structure of NAD(P)H:quinone reductase, a flavoprotein involved in cancer chemoprotection and chemotherapy: mechanism of the two-electron reduction.

Quinone reductase [NAD(P)H:(quinone acceptor) oxidoreductase, EC 1.6.99.2], also called DT diaphorase, is a homodimeric FAD-containing enzyme that catalyzes obligatory NAD(P)H-dependent two-electron reductions of quinones and protects cells against the toxic and neoplastic effects of free radicals and reactive oxygen species arising from one-electron reductions. These two-electron reductions participate in the reductive bioactivation of cancer chemotherapeutic agents such as mitomycin C in tumor cells. Thus, surprisingly, the same enzymatic reaction that protects normal cells activates cytotoxic drugs used in cancer chemotherapy. The 2.1-A crystal structure of rat liver quinone reductase reveals that the folding of a portion of each monomer is similar to that of flavodoxin, a bacterial FMN-containing protein. Two additional portions of the polypeptide chains are involved in dimerization and in formation of the two identical catalytic sites to which both monomers contribute. The crystallographic structures of two FAD-containing enzyme complexes (one containing NADP+, the other containing duroquinone) suggest that direct hydride transfers from NAD(P)H to FAD and from FADH2 to the quinone [which occupies the site vacated by NAD(P)H] provide a simple rationale for the obligatory two-electron reductions involving a ping-pong mechanism.

Structure of human T-cell receptors specific for an immunodominant myelin basic protein peptide: positioning of T-cell receptors on HLA-DR2/peptide complexes.

T-cell receptors (TCRs) recognize peptide bound within the relatively conserved structural framework of major histocompatibility complex (MHC) class I or class II molecules but can discriminate between closely related MHC molecules. The structural basis for the specificity of ternary complex formation by the TCR and MHC/peptide complexes was examined for myelin basic protein (MBP)-specific T-cell clones restricted by different DR2 subtypes. Conserved features of this system allowed a model for positioning of the TCR on DR2/peptide complexes to be developed: (i) The DR2 subtypes that presented the immunodominant MBP peptide differed only at a few polymorphic positions of the DR beta chain. (ii) TCR recognition of a polymorphic residue on the helical portion of the DR beta chain (position DR beta 67) was important in determining the MHC restriction. (iii) The TCR variable region (V) alpha 3.1 gene segment was used by all of the T-cell clones. TCR V beta usage was more diverse but correlated with the MHC restriction--i.e., with the polymorphic DR beta chains. (iv) Two clones with conserved TCR alpha chains but different TCR beta chains had a different MHC restriction but a similar peptide specificity. The difference in MHC restriction between these T-cell clones appeared due to recognition of a cluster of polymorphic DR beta-chain residues (DR beta 67-71). MBP-(85-99)-specific TCRs therefore appeared to be positioned on the DR2/peptide complex such that the TCR beta chain contacted the polymorphic DR beta-chain helix while the conserved TCR alpha chain contacted the nonpolymorphic DR alpha chain.

One sequence, two folds: a metastable structure of CD2.

When expressed as part of a glutathione S-transferase fusion protein the NH2-terminal domain of the lymphocyte cell adhesion molecule CD2 is shown to adopt two different folds. The immunoglobulin superfamily structure of the major (85%) monomeric component has previously been determined by both x-ray crystallography and NMR spectroscopy. We now describe the structure of a second, dimeric, form present in about 15% of recombinant CD2 molecules. After denaturation and refolding in the absence of the fusion partner, dimeric CD2 is converted to monomer, illustrating that the dimeric form represents a metastable folded state. The crystal structure of this dimeric form, refined to 2.0-A resolution, reveals two domains with overall similarity to the IgSF fold found in the monomer. However, in the dimer each domain is formed by the intercalation of two polypeptide chains. Hence each domain represents a distinct folding unit that can assemble in two different ways. In the dimer the two domains fold around a hydrophilic interface believed to mimic the cell adhesion interaction at the cell surface, and the formation of dimer can be regulated by mutating single residues at this interface. This unusual misfolded form of the protein, which appears to result from inter- rather than intramolecular interactions being favored by an intermediate structure formed during the folding process, illustrates that evolution of protein oligomers is possible from the sequence for a single protein domain.

On the possibility of coherently stimulated recombination and cosmological structure generation: recombination instability.

Possible instabilities during cosmological recombination may produce an epoch of nonlinear density growth and fractal-like structural patterns out to the horizon scale at that epoch (approximately 200 Mpc today). With this motivation, we examine the consequences of the change in effective radiative recombination reaction rate coefficients produced by intense stimulated emission. The proton-electron recombination is considered as a natural laser, leading to the formation of spatially nonuniform distributions of neutral matter earlier than the recombination epoch.

On the possibility of coherently stimulated recombination and cosmological structure generation: cosmological consequences.

Given a specific physical mechanism for instabilities during cosmological recombination discussed in an earlier paper, we examine the nonlinear growth of density structures to form fractal-like structural patterns out to the horizon scale at that epoch (approximately 200 Mpc today). A model for such fractal patterns is presented. Such effects could explain observed large-scale structure patterns and the formation of objects at high z, while keeping microwave background anisotropies at the observed minimal levels. We also discuss possible microwave background implications of such a transition and note a potentially observable spectral signature at lambda approximately 0.18 mm as well as a weak line near the peak in the microwave background.

Hairpins are formed by the single DNA strands of the fragile X triplet repeats: structure and biological implications.

Inordinate expansion and hypermethylation of the fragile X DNA triplet repeat, (GGC)n.(GCC)n, are correlated with the ability of the individual G- and C-rich single strands to form hairpin structures. Two-dimensional NMR and gel electrophoresis studies show that both the G- and C-rich single strands form hairpins under physiological conditions. This propensity of hairpin formation is more pronounced for the C-rich strand than for the G-rich strand. This observation suggests that the C-rich strand is more likely to form hairpin or "slippage" structure and show asymmetric strand expansion during replication. NMR data also show that the hairpins formed by the C-rich strands fold in such a way that the cytosine at the CpG step of the stem is C.C paired. The presence of a C.C mismatch at the CpG site generates local flexibility, thereby providing analogs of the transition to the methyltransferase. In other words, the hairpins of the C-rich strand act as better substrates for the human methyltransferase than the Watson-Crick duplex or the G-rich strand. Therefore, hairpin formation could account for the specific methylation of the CpG island in the fragile X repeat that occurs during inactivation of the FMR1 gene during the onset of the disease.