Effects of DNA replication on mRNA noise.

There are several sources of fluctuations in gene expression. Here we study the effects of time-dependent DNA replication, itself a tightly controlled process, on noise in mRNA levels. Stochastic simulations of constitutive and regulated gene expression are used to analyze the time-averaged mean and variation in each case. The simulations demonstrate that to capture mRNA distributions correctly, chromosome replication must be realistically modeled. Slow relaxation of mRNA from the low copy number steady state before gene replication to the high steady state after replication is set by the transcript's half-life and contributes significantly to the shape of the mRNA distribution. Consequently both the intrinsic kinetics and the gene location play an important role in accounting for the mRNA average and variance. Exact analytic expressions for moments of the mRNA distributions that depend on the DNA copy number, gene location, cell doubling time, and the rates of transcription and degradation are derived for the case of constitutive expression and subsequently extended to provide approximate corrections for regulated expression and RNA polymerase variability. Comparisons of the simulated models and analytical expressions to experimentally measured mRNA distributions show that they better capture the physics of the system than previous theories.

Heterogeneity in protein expression induces metabolic variability in a modeled Escherichia coli population.

Stochastic gene expression can lead to phenotypic differences among cells even in isogenic populations growing under macroscopically identical conditions. Here, we apply flux balance analysis in investigating the effects of single-cell proteomics data on the metabolic behavior of an in silico Escherichia coli population. We use the latest metabolic reconstruction integrated with transcriptional regulatory data to model realistic cells growing in a glucose minimal medium under aerobic conditions. The modeled population exhibits a broad distribution of growth rates, and principal component analysis was used to identify well-defined subpopulations that differ in terms of their pathway use. The cells differentiate into slow-growing acetate-secreting cells and fast-growing CO2-secreting cells, and a large population growing at intermediate rates shift from glycolysis to Entner-Doudoroff pathway use. Constraints imposed by integrating regulatory data have a large impact on NADH oxidizing pathway use within the cell. Finally, we find that stochasticity in the expression of only a few genes may be sufficient to capture most of the metabolic variability of the entire population.

Leucyl-tRNA synthetase editing domain functions as a molecular rheostat to control codon ambiguity in Mycoplasma pathogens.

Mycoplasma leucyl-tRNA synthetases (LeuRSs) have been identified in which the connective polypeptide 1 (CP1) amino acid editing domain that clears mischarged tRNAs are missing (Mycoplasma mobile) or highly degenerate (Mycoplasma synoviae). Thus, these enzymes rely on a clearance pathway called pretransfer editing, which hydrolyzes misactivated aminoacyl-adenylate intermediate via a nebulous mechanism that has been controversial for decades. Even as the sole fidelity pathway for clearing amino acid selection errors in the pathogenic M. mobile, pretransfer editing is not robust enough to completely block mischarging of tRNA(Leu), resulting in codon ambiguity and statistical proteins. A high-resolution X-ray crystal structure shows that M. mobile LeuRS structurally overlaps with other LeuRS cores. However, when CP1 domains from different aminoacyl-tRNA synthetases and origins were fused to this common LeuRS core, surprisingly, pretransfer editing was enhanced. It is hypothesized that the CP1 domain evolved as a molecular rheostat to balance multiple functions. These include distal control of specificity and enzyme activity in the ancient canonical core, as well as providing a separate hydrolytic active site for clearing mischarged tRNA.

Naturally occurring aminoacyl-tRNA synthetases editing-domain mutations that cause mistranslation in Mycoplasma parasites.

Mycoplasma parasites escape host immune responses via mechanisms that depend on remarkable phenotypic plasticity. Identification of these mechanisms is of great current interest. The aminoacyl-tRNA synthetases (AARSs) attach amino acids to their cognate tRNAs, but occasionally make errors that substitute closely similar amino acids. AARS editing pathways clear errors to avoid mistranslation during protein synthesis. We show here that AARSs in Mycoplasma parasites have point mutations and deletions in their respective editing domains. The deleterious effect on editing was confirmed with a specific example studied in vitro. In vivo mistranslation was determined by mass spectrometric analysis of proteins produced in the parasite. These mistranslations are uniform cases where the predicted closely similar amino acid replaced the correct one. Thus, natural AARS editing-domain mutations in Mycoplasma parasites cause mistranslation. We raise the possibility that these mutations evolved as a mechanism for antigen diversity to escape host defense systems.

Cytochrome c(2) Exit Strategy: Dissociation Studies and Evolutionary Implications.

Small, water-soluble, type c cytochromes form a transient network connecting major bioenergetic membrane protein complexes in both photosynthesis and respiration. In the photosynthesis cycle of Rhodobacter sphaeroides, cytochrome c2 (cyt c2) docks to the reaction center (RC), undergoes electron transfer, and exits for the cytochrome bc1 complex. Translations of cyt c2 about the RC-cyt c2 docking interface and surrounding membrane reveal possible exit pathways. A pathway at a minimal elevation allowed by the architecture of the RC is analyzed using both an all-atom steered molecular dynamics simulation of the RC-cyt c2 complex and a bioinformatic analysis of the structures and sequences of cyt c. The structure-based phylogenetic analysis allows for the identification of structural elements that have evolved to satisfy the requirements of having multiple functional partners. The patterns of evolutionary variation obtained from the phylogenetic analysis of both docking partners of cyt c2 reveal conservation of key residues involved in the interaction interfaces that would be candidates for further experimental studies. Additionally, using the molecular mechanics Poisson-Boltzmann surface area method we calculate that the binding free energy of reduced cyt c2 to the RC is nearly 6 kcal/mol more favorable than with oxidized cyt c2. The redox-dependent variations lead to changes in structural flexibility, behavior of the interfacial water molecules, and eventually changes in the binding free energy of the complex.

A network of conserved interactions regulates the allosteric signal in a glutamine amidotransferase.

We have combined equilibrium and steered molecular dynamics (SMD) simulations with principal component and correlation analyses to probe the mechanism of allosteric regulation in imidazole glycerol phosphate (IGP) synthase. An evolutionary analysis of IGP synthase revealed a conserved network of interactions leading from the effector binding site to the glutaminase active site, forming conserved communication pathways between the remote active sites. SMD simulations of the undocking of the ribonucleotide effector N1-[(5'-phosphoribulosyl)-formino]-5'-aminoimidazole carboxamide ribonucleotide (PRFAR) resulted in a large scale hinge-opening motion at the interface. Principal component analysis and a correlation analysis of the equilibration protein motion indicate that the dynamics involved in the allosteric transition are mediated by coupled motion between sites that are more than 25 A apart. Furthermore, conserved residues at the substrate-binding site, within the barrel, and at the interface were found to exhibit highly correlated motion during the allosteric transition. The coupled motion between PRFAR unbinding and the directed opening of the interface is interpreted in combination with kinetic assays for the wild-type and mutant systems to develop a model of allosteric regulation in IGP synthase that is monitored and investigated with atomic resolution.

The evolutionary history of Cys-tRNACys formation.

The recent discovery of an alternate pathway for indirectly charging tRNA(Cys) has stimulated a re-examination of the evolutionary history of Cys-tRNA(Cys) formation. In the first step of the pathway, O-phosphoseryl-tRNA synthetase charges tRNA(Cys) with O-phosphoserine (Sep), a precursor of the cognate amino acid. In the following step, Sep-tRNA:Cys-tRNA synthase (SepCysS) converts Sep to Cys in a tRNA-dependent reaction. The existence of such a pathway raises several evolutionary questions, including whether the indirect pathway is a recent evolutionary invention, as might be implied from its localization to the Euryarchaea, or, as evidence presented here indicates, whether this pathway is more ancient, perhaps already in existence at the time of the last universal common ancestral state. A comparative phylogenetic approach is used, combining evolutionary information from protein sequences and structures, that takes both the signature of horizontal gene transfer and the recurrence of the full canonical phylogenetic pattern into account, to document the complete evolutionary history of cysteine coding and understand the nature of this process in the last universal common ancestral state. Resulting from the historical study of tRNA(Cys) aminoacylation and the integrative perspective of sequence, structure, and function are 3D models of O-phosphoseryl-tRNA synthetase and SepCysS, which provide experimentally testable predictions regarding the identity and function of key active-site residues in these proteins. The model of SepCysS is used to suggest a sulfhydrylation reaction mechanism, which is predicted to occur at the interface of a SepCysS dimer.

Reaction coupling through interdomain contacts in imidazole glycerol phosphate synthase.

Imidazole glycerol phosphate (IGP) synthase, a triad glutamine amidotransferase, catalyzes the fifth step in the histidine biosynthetic pathway, where ammonia from glutamine is incorporated into N1-[(5'-phosphoribulosyl)-formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (PRFAR) to yield IGP and 5'-(5-aminoimidazole-4-carboxamide) ribonucleotide (AICAR). The triad family of glutamine amidotransferases is formed by the coupling of two disparate subdomains, an acceptor domain and a glutamine hydrolysis domain. Each of the enzymes in this family share a common glutaminase domain for which the glutaminase activity is tightly regulated by an acceptor substrate domain. In IGP synthase the glutaminase and PRFAR binding sites are separated by 30 A. Using kinetic analyses of site-specific mutants and molecular dynamic simulations, we have determined that an interdomain salt bridge in IGP synthase between D359 and K196 (approximately 16 A from the PRFAR binding site) plays a key role in mediating communication between the two active sites. This interdomain contact modulates the glutaminase loop containing the histidine and glutamic acid of the catalytic triad to control glutamine hydrolysis. We propose this to be a general principle of catalytic coupling that may be applied to the entire triad glutamine amidotransferase family.

Structural elements in IGP synthase exclude water to optimize ammonia transfer.

In the complex pathway of histidine biosynthesis, a key branch point linking amino acid and purine biosynthesis is catalyzed by the bifunctional enzyme imidazole glycerol phosphate (IGP) synthase. The first domain of IGP synthase, a triad glutamine amidotransferase, hydrolyzes glutamine to form glutamate and ammonia. Its activity is tightly regulated by the binding of the substrate PRFAR to its partner synthase domain. Recent crystal structures and molecular dynamics simulations strongly suggest that the synthase domain, a (beta/alpha)(8) barrel protein, mediates the insertion of ammonia and ring formation in IGP by channeling ammonia from one remote active site to the other. Here, we combine both mutagenesis experiments and computational investigations to gain insight into the transfer of ammonia and the mechanism of conduction. We discover an alternate route for the entrance of ammonia into the (beta/alpha)(8) barrel and argue that water acts as both agonist and antagonist to the enzymatic function. Our results indicate that the architecture of the two subdomains, most notably the strict conservation of key residues at the interface and within the (beta/alpha)(8) barrel, has been optimized to allow the efficient passage of ammonia, and not water, between the two remote active sites.

Optimizing physical energy functions for protein folding.

We optimize a physical energy function for proteins with the use of the available structural database and perform three benchmark tests of the performance: (1) recognition of native structures in the background of predefined decoy sets of Levitt, (2) de novo structure prediction using fragment assembly sampling, and (3) molecular dynamics simulations. The energy parameter optimization is based on the energy landscape theory and uses a Monte Carlo search to find a set of parameters that seeks the largest ratio deltaE(s)/DeltaE for all proteins in a training set simultaneously. Here, deltaE(s) is the stability gap between the native and the average in the denatured states and DeltaE is the energy fluctuation among these states. Some of the energy parameters optimized are found to show significant correlation with experimentally observed quantities: (1) In the recognition test, the optimized function assigns the lowest energy to either the native or a near-native structure among many decoy structures for all the proteins studied. (2) Structure prediction with the fragment assembly sampling gives structure models with root mean square deviation less than 6 A in one of the top five cluster centers for five of six proteins studied. (3) Structure prediction using molecular dynamics simulation gives poorer performance, implying the importance of having a more precise description of local structures. The physical energy function solely inferred from a structural database neither utilizes sequence information from the family of the target nor the outcome of the secondary structure prediction but can produce the correct native fold for many small proteins.

Is the olfactory receptor a metalloprotein?

The sense of smell is arguably our most primal faculty and also the least understood. Even our own olfactorily impaired species is capable of detecting approximately 10,000 distinct scents [Buck, L. & Axel, R. (1991) Cell 65, 175-187]. To achieve that amazing diversity, mammals have approximately 1,000 olfactory genes, which accounts for approximately 3% of their entire genome [Mombaerts, P. (1999) Science 286, 707-711]. The olfactory receptors (ORs) are believed to be seven-helix transmembrane proteins, with an odorant-binding site on the periplasmic domain and a G protein-binding site on the cytoplasmic domain. Odorants first bind to an OR, which then undergoes some structural change that triggers the G protein activation and the following cascade of events leading to nerve cell activity. The structural details of ORs, however, remain to be determined. In this paper, we will describe a hypothesis in which metal ions play an important role for odorant recognition. We analyze the predicted structure and consensus sequence of the ORs and propose a metal-binding site in the loop between fourth and fifth helix (4-5 loop). We have prepared synthetically a pentapeptide that contains this putative binding site and find that it not only has high affinity for binding Cu(II) and Zn(II) ions, but that it also undergoes a dramatic transition to an alpha-helical structure upon metal ion binding. Based on these observations, we propose a "shuttlecock" mechanism for the possible structural change in ORs upon odorant binding. This mechanism involves membrane penetration of the 4-5 loop after residue charge neutralization by metal ion binding.