In vivo assembly of aspartate transcarbamoylase from fragmented and circularly permuted catalytic polypeptide chains

Previous studies on Escherichia coli aspartate transcarbamoylase (ATCase) demonstrated that active, stable enzyme was formed in vivo from complementing polypeptides of the catalytic (c) chain encoded by gene fragments derived from the pyrBI operon. However, the enzyme lacked the allosteric properties characteristic of wild-type ATCase. In order to determine whether the loss of homotropic and heterotropic properties was attributable to the location of the interruption in the polypeptide chain rather than to the lack of continuity, we constructed a series of fragmented genes so that the breaks in the polypeptide chains would be dispersed in different domains and diverse regions of the structure. Also, analogous molecules containing circularly permuted c chains with altered termini were constructed for comparison with the ATCase molecules containing fragmented c chains. Studies were performed on four sets of ATCase molecules containing cleaved c chains at positions between residues 98 and 99, 121 and 122, 180 and 181, and 221 and 222; the corresponding circularly permuted chains had N termini at positions 99, 122, 181, and 222. All of the ATCase molecules containing fragmented or circularly permuted c chains exhibited the homotropic and heterotropic properties characteristic of the wild-type enzyme. Hill coefficients (n(H:)) and changes in them upon the addition of ATP and CTP were similar to those observed with wild-type ATCase. In addition, the conformational changes revealed by the decrease in sedimentation coefficient upon the addition of a bisubstrate analog were virtually identical to that for the wild-type enzyme. Differential scanning calorimetry showed that neither the breakage of the polypeptide chains nor the newly formed covalent bond between the termini in the wild-type enzyme had a significant impact on the thermal stability of the assembled dodecamers. The studies demonstrate that continuity of the polypeptide chain within structural domains is not essential for the assembly, activity, and allosteric properties of ATCase.     

240s loop interactions stabilize the T state of Escherichia coli aspartate transcarbamoylase

Here the functional and structural importance of interactions involving the 240s loop of the catalytic chain for the stabilization of the T state of aspartate transcarbamoylase were tested by replacement of Lys-244 with Asn and Ala. For the K244A and K244N mutant enzymes, the aspartate concentration required to achieve half-maximal specific activity was reduced to 8.4 and 4.0 mm, respectively, as compared with 12.4 mM for the wild-type enzyme. Both mutant enzymes exhibited dramatic reductions in homotropic cooperativity and the ability of the heterotropic effectors to modulate activity. Small angle x-ray scattering studies showed that the unligated structure of the mutant enzymes, and the structure of the mutant enzymes ligated with N-phosphonacetyl-L-aspartate, were similar to that observed for the unligated and N-phosphonacetyl-L-aspartateligated wild-type enzyme. A saturating concentration of carbamoyl phosphate alone has little influence on the small angle x-ray scattering of the wild-type enzyme. However, carbamoyl phosphate was able to shift the structure of the two mutant enzymes dramatically toward R, establishing that the mutations had destabilized the T state of the enzyme. The x-ray crystal structure of K244N enzyme showed that numerous local T state stabilizing interactions involving 240s loop residues were lost. Furthermore, the structure established that the mutation induced additional alterations at the subunit interfaces, the active site, the relative position of the domains of the catalytic chains, and the allosteric domain of the regulatory chains. Most of these changes reflect motions toward the R state structure. However, the K244N mutation alone only changes local conformations of the enzyme to an R-like structure, without triggering the quaternary structural transition. These results suggest that loss of cooperativity and reduction in heterotropic effects is due to the dramatic destabilization of the T state of the enzyme by this mutation in the 240s loop of the catalytic chain.     

Mechanism of insulin chain combination. Asymmetric roles of A-chain alpha-helices in disulfide pairing

The A and B chains of insulin combine to form native disulfide bridges without detectable isomers. The fidelity of chain combination thus recapitulates the folding of proinsulin, a precursor protein in which the two chains are tethered by a disordered connecting peptide. We have recently shown that chain combination is blocked by seemingly conservative substitutions in the C-terminal alpha-helix of the A chain. Such analogs, once formed, nevertheless retain high biological activity. By contrast, we demonstrate here that chain combination is robust to non-conservative substitutions in the N-terminal alpha-helix. Introduction of multiple glycine substitutions into the N-terminal segment of the A chain (residues A1-A5) yields analogs that are less stable than native insulin and essentially without biological activity. (1)H NMR studies of a representative analog lacking invariant side chains Ile(A2) and Val(A3) (A chain sequence GGGEQCCTSICSLYQLENYCN; substitutions are italicized and cysteines are underlined) demonstrate local unfolding of the A1-A5 segment in an otherwise native-like structure. That this and related partial folds retain efficient disulfide pairing suggests that the native N-terminal alpha-helix does not participate in the transition state of the reaction. Implications for the hierarchical folding mechanisms of proinsulin and insulin-like growth factors are discussed.     

Importance of domain closure for the catalysis and regulation of Escherichia coli aspartate transcarbamoylase

Two hybrid versions of Escherichia coli aspartate transcarbamoylase were studied to determine the influence of domain closure on the homotropic and heterotropic properties of the enzyme. Each hybrid holoenzyme had one wild-type and one inactive catalytic subunit. In the first case the inactive catalytic subunit had Arg-54 replaced by alanine. The holoenzyme with this mutation in all six catalytic chains exhibits a 17,000-fold reduction in activity, no loss in substrate affinity, and an R state structurally identical to that of the wild-type enzyme. In the second case, the inactive catalytic subunit had Arg-105 replaced by alanine. The holoenzyme with this mutation in all six catalytic chains exhibits a 1,100-fold reduction in activity, substantial loss in substrate affinity, and loss of the ability to be converted to the R state. Thus, the R54A substitution results in a holoenzyme that can undergo closure of the catalytic chain domains to form the high activity, high affinity active site and to undergo the allosteric transition, whereas the R105A substitution results in a holoenzyme that can neither undergo domain closure nor the allosteric transition. The hybrid holoenzyme with one wild-type and one R54A catalytic subunit exhibited the same maximal velocity per active site as the wild-type holoenzyme, reduced cooperativity, and normal heterotropic interactions. The hybrid with one wild-type and one R105A catalytic subunit exhibited significantly reduced maximal velocity per active site as compared with the wild-type holoenzyme, reduced cooperativity, and substantially reduced heterotropic interactions. Small angle x-ray scattered was used to verify that the R105A-containing hybrid could attain an R state structure. These results indicate the global nature of the conformational changes associated with the allosteric transition in the enzyme. If one catalytic subunit cannot undergo domain closure to create the active sites, then the entire molecule cannot attain the high activity, high activity R state.     

Stochasticity in Protein Levels Drives Colinearity of Gene Order in Metabolic Operons of Escherichia coli

In bacterial genomes, gene order is not random. This is most evident when looking at operons, these often encoding enzymes involved in the same metabolic pathway or proteins from the same complex. Is gene order within operons nonrandom, however, and if so why? We examine this issue using metabolic operons as a case study. Using the metabolic network of Escherichia coli, we define the temporal order of reactions. We find a pronounced trend for genes to appear in operons in the same order as they are needed in metabolism (colinearity). This is paradoxical as, at steady state, enzymes abundance should be independent of order within the operon. We consider three extensions of the steady-state model that could potentially account for colinearity: (1) increased productivity associated with higher expression levels of the most 5′ genes, (2) a faster metabolic processing immediately after up-regulation, and (3) metabolic stalling owing to stochastic protein loss. We establish the validity of these hypotheses by employing deterministic and stochastic models of enzyme kinetics. The stochastic stalling hypothesis correctly and uniquely predicts that colinearity is more pronounced both for lowly expressed operons and for genes that are not physically adjacent. The alternative models fail to find any support. These results support the view that stochasticity is a pervasive problem to a cell and that gene order evolution can be driven by the selective consequences of fluctuations in protein levels.     

Variation in 16S-23S rRNA intergenic spacer regions among Bacillus subtilis 168 isolates

The genome of the Bacillus subtilis 168-type strain contains 10 ribosomal RNA (rRNA) operons. In the intergenic spacer region (ISR) between the 16S and 23S rRNA genes, five rRNA operons, rrnI-H-G and rrnJ-W, lack a trinucleotide signature region. Precise determination of molecular weight (MW), using electrospray mass spectrometry (MS), of the polymerase chain reaction (PCR) products from a segment of the ISR from the 168-type strain and B. subtilis 168-like strain 23071 demonstrated 114 and 111 basepair (bp) PCR products (due to the presence or absence of the insert in the operons) as predicted from sequence. However, PCR of the ISR segment for five other B. subtilis 168 isolates generated only a 114 bp PCR product, suggesting the presence of the trinucleotide signature region in all rRNA operons for these strains. Additional genetic variability between the seven B. subtilis 168 isolates was demonstrated by restriction fragment length polymorphism (RFLP) of the rRNA operons, with three distinct patterns found upon Southern blot analysis. The 168-type strain and three others (23066, 23067, and 23071) exhibited the same Southern pattern. Thus, operon deletion is not responsible for the absence of a 111 bp product on MS analysis for strains 23066 and 23067. Restriction analysis confirmed the presence of the trinucleotide signature region in the ISR of all rRNA operons for five B. subtilis 168 isolates; sequencing of rrnW/H from a representative strain also upheld this finding. These results help provide a better understanding of variations in sequence, operon number and chromosomal organization, both within a genome and among isolates of B. subtilis subgroup 168. It is also hypothesized that the presence of the trinucleotide insert in certain rRNA operons may play a role in rRNA maturation and protein synthesis.     

Biodegradation of halogenated organic compounds.

In this review we discuss the degradation of chlorinated hydrocarbons by microorganisms, emphasizing the physiological, biochemical, and genetic basis of the biodegradation of aliphatic, aromatic, and polycyclic compounds. Many environmentally important xenobiotics are halogenated, especially chlorinated. These compounds are manufactured and used as pesticides, plasticizers, paint and printing-ink components, adhesives, flame retardants, hydraulic and heat transfer fluids, refrigerants, solvents, additives for cutting oils, and textile auxiliaries. The hazardous chemicals enter the environment through production, commercial application, and waste. As a result of bioaccumulation in the food chain and groundwater contamination, they pose public health problems because many of them are toxic, mutagenic, or carcinogenic. Although synthetic chemicals are usually recalcitrant to biodegradation, microorganisms have evolved an extensive range of enzymes, pathways, and control mechanisms that are responsible for catabolism of a wide variety of such compounds. Thus, such biological degradation can be exploited to alleviate environmental pollution problems. The pathways by which a given compound is degraded are determined by the physical, chemical, and microbiological aspects of a particular environment. By understanding the genetic basis of catabolism of xenobiotics, it is possible to improve the efficacy of naturally occurring microorganisms or construct new microorganisms capable of degrading pollutants in soil and aquatic environments more efficiently. Recently a number of genes whose enzyme products have a broader substrate specificity for the degradation of aromatic compounds have been cloned and attempts have been made to construct gene cassettes or synthetic operons comprising these degradative genes. Such gene cassettes or operons can be transferred into suitable microbial hosts for extending and custom designing the pathways for rapid degradation of recalcitrant compounds. Recent developments in designing recombinant microorganisms and hybrid metabolic pathways are discussed.     

Coupling oscillations and switches in genetic networks

Switches (bistability) and oscillations (limit cycle) are omnipresent in biological networks. Synthetic genetic networks producing bistability and oscillations have been designed and constructed experimentally. However, in real biological systems, regulatory circuits are usually interconnected and the dynamics of those complex networks is often richer than the dynamics of simple modules. Here we couple the genetic Toggle switch and the Repressilator, two prototypic systems exhibiting bistability and oscillations, respectively. We study two types of coupling. In the first type, the bistable switch is under the control of the oscillator. Numerical simulation of this system allows us to determine the conditions under which a periodic switch between the two stable steady states of the Toggle switch occurs. In addition we show how birhythmicity characterized by the coexistence of two stable small-amplitude limit cycles, can easily be obtained in the system. In the second type of coupling, the oscillator is placed under the control of the Toggleswitch. Numerical simulation of this system shows that this construction could for example be exploited to generate a permanent transition from a stable steady state to self-sustained oscillations (and vice versa) after a transient external perturbation. Those results thus describe qualitative dynamical behaviors that can be generated through the coupling of two simple network modules. These results differ from the dynamical properties resulting from interlocked feedback loops systems in which a given variable is involved at the same time in both positive and negative feedbacks. Finally the models described here may be of interest in synthetic biology, as they give hints on how the coupling should be designed to get the required properties.     

Mapping of genes involved in macromolecular synthesis on the chromosome of Streptomyces coelicolor A3(2).

The genes for the beta, beta', and seven sigma factor subunits of RNA polymerase, for elongation factors EF-Tu1 and EF-Tu3, and for six rRNA operons were mapped on the combined genetic and physical map of the Streptomyces coelicolor chromosome. Like the previously mapped tRNA genes, the RNA polymerase and rRNA genes map to scattered positions. The lack of rRNA operons in the immediate vicinity of the origin of replication (oriC) and the absence of tRNA genes in any of the rRNA operons are novel features of the Streptomyces chromosome.     

DNA shuffling: Modifying the hand that nature dealt

Summary DNA shuffling is a technique being utilized for in vitro recombination of a single gene or pools of homologous genes. The genes are fragmented into randomly sized pieces, and polymerase chain reaction (PCR) reassembly of full-length genes from the fragments, via self-priming, yields recombination due to PCR template switching. After these PCR products are screened and the interesting products sequenced, improved clones are reshuffled to recombine useful mutations in additive or synergistic ways, in effect mimicking the process of natural sexual recombination. Proteins can be ‘bred’ with the appropriate individual properties and then their ‘progeny’ screened for the desired combination of traits. DNA shuffling is a powerful tool enabling rapid and directed evolution of new genes, operons and whole viral genomes.     

In vivo formation of allosteric aspartate transcarbamoylase containing circularly permuted catalytic polypeptide chains: implications for protein folding and assembly.

Because the N- and C-terminal amino acids of the catalytic (c) polypeptide chains of Escherichia coli aspartate transcarbamoylase (ATCase) are in close proximity to each other, it has been possible to form in vivo five different active ATCase variants in which the terminal regions of the wild-type c chains are linked in a continuous polypeptide chain and new termini are introduced elsewhere in either of the two structural domains of the c chain. These circularly permuted (cp) chains were produced by constructing tandem pyrB genes, which encode the c chain of ATCase, followed by application of PCR. Chains expressed in this way assemble efficiently in vivo to form active, stable ATCase variants. Three such variants have been purified and shown to have the kinetic and physical properties characteristic of wild-type ATCase composed of two catalytic (C) trimers and three regulatory (R) dimers. The values of Vmax for cpATCase122, cpATCase222, and cpATCase281 ranged from 16-21 mumol carbamoylaspartate per microgram per h, compared with 15 for wild-type ATCase, and the values for K0.5 for the variants were 4-17 mM aspartate, whereas wild-type ATCase exhibited a value of 6 mM. Hill coefficients for the three variants varied from 1.8 to 2.1, compared with 1.4 for the wild-type enzyme. As observed with wild-type ATCase, ATP activated the variants containing the circularly permuted chains, as shown by the lowering of K0.5 for aspartate and a decrease in the Hill coefficient (nH). In contrast, CTP caused both an increase in K0.5 and nH for the variants, just as observed with wild-type ATCase. Thus, the enzyme containing the permuted chains with widely diverse N- and C-termini exhibited the homotropic and heterotropic effects characteristic of wild-type ATCase. The decrease in the sedimentation coefficient of the variants caused by the binding of the bisubstrate ligand N-(phosphonacetyl)-L-aspartate (PALA) was also virtually identical to that obtained with wild-type ATCase, thereby indicating that these altered ATCase molecules undergo the analogous ligand-promoted allosteric transition from the taut (T) state to the relaxed (R) conformation. These ATCase molecules with new N- and C-termini widely dispersed throughout the c chains are valuable models for studying in vivo and in vitro folding of polypeptide chains.     

Physical mapping of stable RNA genes in Bacillus subtilis using polymerase chain reaction amplification from a yeast artificial chromosome library.

A new approach for mapping genes which utilizes yeast artificial chromosome clones carrying parts of the Bacillus subtilis genome and the polymerase chain reaction technique is described. This approach was used to physically map stable RNA genes of B. subtilis. Results from over 400 polymerase chain reactions carried out with the yeast artificial chromosome clone library, using primers specific for the genes of interest and designed from published sequences, were collected. The locations of 10 known rRNA gene regions (rrnO, rrnA, rrnE, rrnD, rrnB, rrnJ-rrnW, and rrnI-rrnH-rrnG) have been determined by this method, and these results correlate with those observed by standard genetic mapping. All rRNA operons, except rrnB, are found between 0 and 90 degrees, while rrnB has been placed in the area of 270 degrees on the chromosome map. Also localized were the tRNA gene clusters associated with the following ribosomal operons: rrnB (21 tRNAs), rrnJ (9 tRNAs), rrnD (16 tRNAs), and rrnO and rrnA (2 internal tRNAs). A previously unmapped four-tRNA gene cluster, trnY, a tRNA gene region that is not associated with a ribosomal operon, was found near the origin of replication. The P-RNA gene, important for processing of tRNAs, was found between map locations 197 and 204 degrees.     

Ligand binding and allostery can emerge simultaneously

A heterotropic allosteric effect involves an effector molecule that is distinct from the substrate or ligand of the protein. How heterotropic allostery originates is an unanswered question. We have previously created several heterotropic allosteric enzymes by recombining the genes for TEM1 beta-lactamase (BLA) and maltose binding protein (MBP) to create BLAs that are positively or negatively regulated by maltose. We show here that one of these engineered enzymes has approximately 10(6) M(-1) affinity for Zn(2+), a property that neither of the parental proteins possesses. Furthermore, Zn(2+) is a negative effector that noncompetitively switches off beta-lactam hydrolysis activity. Mutagenesis experiments indicate that the Zn(2+)-binding site does not involve a histidine or a cysteine, which is atypical of natural Zn(2+)-binding sites. These studies also implicate helices 1 and 12 of the BLA domain in allosteric signal propagation. These results support a model for the evolution of heterotropic allostery in which effector affinity and allosteric signaling emerge simultaneously.     

Different functional modulation by heterotropic ligands (2,3-diphosphoglycerate and chlorides) of the two haemoglobins from fallow-deer (Dama dama).

Two haemoglobin components have been identified and purified from fallow-deer (Dama dama) erythrocytes. They are present in similar amounts and the two tetrameric molecules share the same alpha chain, while two different beta chains are detected in the two components. The beta chains differ by 14 residues, even though they both have 145 amino-acid residues, which account for a molecular mass of 16,023 and 16,064 Da, respectively, while alpha chain has 141 residues, yielding a molecular mass of 15,142 Da. Compared with human Hb, the N-terminal region of both beta chains shows deletion of Val beta 1 and the replacement of His beta 2 by a methionyl residue, a modification which is common to most ruminant haemoglobins. Although both isolated components show a low intrinsic affinity for oxygen, meaningful differences between the two haemoglobins have been found with respect to the effect of heterotropic effectors, such as 2,3-diphosphoglycerate and chloride ions. In view of the high sequence homology between the two components, the different effect of heterotropic ligands has been tentatively correlated to possible localized structural variations between beta chains of the two haemoglobin components.     

On the evolution of mutation in changing environments: recombination and phenotypic switching

Phenotypic switching has been observed in laboratory studies of yeast and bacteria, in which the rate of such switching appears to adjust to match the frequency of environmental changes. Among possible mechanisms of switching are epigenetic influences on gene expression and variation in levels of methylation; thus environmental and/or genetic factors may contribute to the rate of switching. Most previous analyses of the evolution of phenotypic switching have compared exponential growth rates of noninteracting populations, and recombination has been ignored. Our genetic model of the evolution of switching rates is framed in terms of a mutation-modifying gene, environments that cause periodic changes in fitness, and recombination between the mutation modifier and the gene under selection. Exact results are obtained for all recombination rates and symmetric fitnesses that strongly generalize earlier results obtained under complete linkage and strong constraints on the relation between fitness and period of switching. Our analytical and numerical results suggest a general principle that recombination reduces the stable rate of switching in symmetric and asymmetric fitness regimes and when the period of switching is random. As the recombination rate increases, it becomes less likely that there is a stable nonzero rate of switching.