What determines sub-diffusive behavior in crowded protein solutions?

The aqueous environment inside cells is densely packed. A typical cell has a macromolecular concentration in the range 90-450 g/L, with 5%-40% of its volume being occupied by macromolecules, resulting in what is known as macromolecular crowding. The space available for the free diffusion of metabolites and other macromolecules is thus greatly reduced, leading to so-called excluded volume effects. The slow diffusion of macromolecules under crowded conditions has been explained using transient complex formation. However, sub-diffusion noted in earlier works is not well characterized, particularly the role played by transient complex formation and excluded volume effects. We have used Brownian dynamics simulations to characterize the diffusion of chymotrypsin inhibitor 2 in protein solutions of bovine serum albumin and lysozyme at concentrations ranging from 50 to 300 g/L. The predicted changes in diffusion coefficient as a function of crowder concentration are consistent with NMR experiments. The sub-diffusive behavior observed in the sub-microsecond timescale can be explained in terms of a so-called cage effect, arising from rattling motion in a local molecular cage as a consequence of excluded volume effects. By selectively manipulating the nature of interactions between protein molecules, we determined that excluded volume effects induce sub-diffusive dynamics at sub-microsecond timescales. These findings may help to explain the diffusion-mediated effects of protein crowding on cellular processes.

Two-species totally asymmetric simple exclusion process model: From a simple description to intermittency and traveling traffic jams.

We extend the paradigmatic and versatile totally asymmetric simple exclusion process (TASEP) for stochastic 1D transport to allow for two different particle species, each having specific entry and exit rates. We offer a complete mean-field analysis, including a phase diagram, by mapping this model onto an effective one-species TASEP. Stochastic simulations confirm the results, but indicate deviations when the particle species have very different exit rates. We illustrate that this is due to a phenomenon of intermittency, and formulate a refined "intermittent" mean-field theory for this regime. We discuss how nonstationary effects may further enrich the phenomenology.

The molecular aetiology of tRNA synthetase depletion: induction of a GCN4 amino acid starvation response despite homeostatic maintenance of charged tRNA levels.

During protein synthesis, charged tRNAs deliver amino acids to translating ribosomes, and are then re-charged by tRNA synthetases (aaRS). In humans, mutant aaRS cause a diversity of neurological disorders, but their molecular aetiologies are incompletely characterised. To understand system responses to aaRS depletion, the yeast glutamine aaRS gene (GLN4) was transcriptionally regulated using doxycycline by tet-off control. Depletion of Gln4p inhibited growth, and induced a GCN4 amino acid starvation response, indicative of uncharged tRNA accumulation and Gcn2 kinase activation. Using a global model of translation that included aaRS recharging, Gln4p depletion was simulated, confirming slowed translation. Modelling also revealed that Gln4p depletion causes negative feedback that matches translational demand for Gln-tRNAGln to aaRS recharging capacity. This maintains normal charged tRNAGln levels despite Gln4p depletion, confirmed experimentally using tRNA Northern blotting. Model analysis resolves the paradox that Gln4p depletion triggers a GCN4 response, despite maintenance of tRNAGln charging levels, revealing that normally, the aaRS population can sequester free, uncharged tRNAs during aminoacylation. Gln4p depletion reduces this sequestration capacity, allowing uncharged tRNAGln to interact with Gcn2 kinase. The study sheds new light on mutant aaRS disease aetiologies, and explains how aaRS sequestration of uncharged tRNAs can prevent GCN4 activation under non-starvation conditions.

Definition of the Minimal Contents for the Molecular Simulation of the Yeast Cytoplasm.

The cytoplasm is a densely packed environment filled with macromolecules with hindered diffusion. Molecular simulation of the diffusion of biomolecules under such macromolecular crowding conditions requires the definition of a simulation cell with a cytoplasmic-like composition. This has been previously done for prokaryote cells (E. coli) but not for eukaryote cells such as yeast as a model organism. Yeast proteomics datasets vary widely in terms of cell growth conditions, the technique used to determine protein composition, the reported relative abundance of proteins, and the units in which abundances are reported. We determined that the gene ontology profiles of the most abundant proteins across these datasets are similar, but their abundances vary greatly. To overcome this problem, we chose five mass spectrometry proteomics datasets that fulfilled the following criteria: high internal consistency, consistency with published experimental data, and freedom from GFP-tagging artifacts. Using these datasets, the contents of a simulation cell containing a single 80S ribosome were defined, such that the macromolecular density and the mass ratio of ribosomal-to-cytoplasmic proteins were consistent with experiment and chosen datasets. Finally, multiple tRNAs were added, consistent with their experimentally-determined number in the yeast cell. The resulting composition can be readily used in molecular simulations representative of yeast cytoplasmic macromolecular crowding conditions to characterize a variety of phenomena, such as protein diffusion, protein-protein interactions and biological processes such as protein translation.

Destabilization of Eukaryote mRNAs by 5' Proximal Stop Codons Can Occur Independently of the Nonsense-Mediated mRNA Decay Pathway.

In eukaryotes, the binding of poly(A) binding protein (PAB) to the poly(A) tail is central to maintaining mRNA stability. PABP interacts with the translation termination apparatus, and with eIF4G to maintain 3'-5' mRNA interactions as part of an mRNA closed loop. It is however unclear how ribosome recycling on a closed loop mRNA is influenced by the proximity of the stop codon to the poly(A) tail, and how post-termination ribosome recycling affects mRNA stability. We show that in a yeast disabled for nonsense mediated mRNA decay (NMD), a PGK1 mRNA with an early stop codon at codon 22 of the reading frame is still highly unstable, and that this instability cannot be significantly countered even when 50% stop codon readthrough is triggered. In an NMD-deficient mutant yeast, stable reporter alleles with more 3' proximal stop codons could not be rendered unstable through Rli1-depletion, inferring defective Rli1 ribosome recycling is insufficient in itself to trigger mRNA instability. Mathematical modelling of a translation system including the effect of ribosome recycling and poly(A) tail shortening supports the hypothesis that impaired ribosome recycling from 5' proximal stop codons may compromise initiation processes and thus destabilize the mRNA. A model is proposed wherein ribosomes undergo a maturation process during early elongation steps, and acquire competency to re-initiate on the same mRNA as translation elongation progresses beyond the very 5' proximal regions of the mRNA.

Novel mRNA-specific effects of ribosome drop-off on translation rate and polysome profile.

The well established phenomenon of ribosome drop-off plays crucial roles in translational accuracy and nutrient starvation responses during protein translation. When cells are under stress conditions, such as amino acid starvation or aminoacyl-tRNA depletion due to a high level of recombinant protein expression, ribosome drop-off can substantially affect the efficiency of protein expression. Here we introduce a mathematical model that describes the effects of ribosome drop-off on the ribosome density along the mRNA and on the concomitant protein synthesis rate. Our results show that ribosome premature termination may lead to non-intuitive ribosome density profiles, such as a ribosome density which increases from the 5' to the 3' end. Importantly, the model predicts that the effects of ribosome drop-off on the translation rate are mRNA-specific, and we quantify their resilience to drop-off, showing that the mRNAs which present ribosome queues are much less affected by ribosome drop-off than those which do not. Moreover, among those mRNAs that do not present ribosome queues, resilience to drop-off correlates positively with the elongation rate, so that sequences using fast codons are expected to be less affected by ribosome drop-off. This result is consistent with a genome-wide analysis of S. cerevisiae, which reveals that under favourable growth conditions mRNAs coding for proteins involved in the translation machinery, known to be highly codon biased and using preferentially fast codons, are highly resilient to ribosome drop-off. Moreover, in physiological conditions, the translation rate of mRNAs coding for regulatory, stress-related proteins, is less resilient to ribosome drop-off. This model therefore allows analysis of variations in the translational efficiency of individual mRNAs by accounting for the full range of known ribosome behaviours, as well as explaining mRNA-specific variations in ribosome density emerging from ribosome profiling studies.

Iron-Catalyzed Synthesis of Sulfur-Containing Heterocycles.

An iron-catalyzed synthesis of sulfur- and sulfone-containing heterocycles is reported. The method is based on the cyclization of readily available substrates and proceeded with high efficiency and diastereoselectivity. A variety of sulfur-containing heterocycles bearing moieties suitable for subsequent functionalization are prepared. Illustrative examples of such postcyclization modifications are also presented.

Identification of the mRNA targets of tRNA-specific regulation using genome-wide simulation of translation.

tRNA gene copy number is a primary determinant of tRNA abundance and therefore the rate at which each tRNA delivers amino acids to the ribosome during translation. Low-abundance tRNAs decode rare codons slowly, but it is unclear which genes might be subject to tRNA-mediated regulation of expression. Here, those mRNA targets were identified via global simulation of translation. In-silico mRNA translation rates were compared for each mRNA in both wild-type and a [Formula: see text] sup70-65 mutant, which exhibits a pseudohyphal growth phenotype and a 75% slower CAG codon translation rate. Of 4900 CAG-containing mRNAs, 300 showed significantly reduced in silico translation rates in a simulated tRNA mutant. Quantitative immunoassay confirmed that the reduced translation rates of sensitive mRNAs were [Formula: see text] concentration-dependent. Translation simulations showed that reduced [Formula: see text] concentrations triggered ribosome queues, which dissipated at reduced translation initiation rates. To validate this prediction experimentally, constitutive gcn2 kinase mutants were used to reduce in vivo translation initiation rates. This repaired the relative translational rate defect of target mRNAs in the sup70-65 background, and ameliorated sup70-65 pseudohyphal growth phenotypes. We thus validate global simulation of translation as a new tool to identify mRNA targets of tRNA-specific gene regulation.

The fluorescent two-hybrid assay to screen for protein-protein interaction inhibitors in live cells: targeting the interaction of p53 with Mdm2 and Mdm4.

Protein-protein interactions (PPIs) are attractive but challenging targets for drug discovery. To overcome numerous limitations of the currently available cell-based PPI assays, we have recently established a fully reversible microscopy-assisted fluorescent two-hybrid (F2H) assay. The F2H assay offers a fast and straightforward readout: an interaction-dependent co-localization of two distinguishable fluorescent signals at a defined spot in the nucleus of mammalian cells. We developed two reversible F2H assays for the interactions between the tumor suppressor p53 and its negative regulators, Mdm2 and Mdm4. We then performed a pilot F2H screen with a subset of compounds, including small molecules (such as Nutlin-3) and stapled peptides. We identified five cell-penetrating compounds as potent p53-Mdm2 inhibitors. However, none exhibited intracellular activity on p53-Mdm4. Live cell data generated by the F2H assays enable the characterization of stapled peptides based on their ability to penetrate cells and disrupt p53-Mdm2 interaction as well as p53-Mdm4 interaction. Here, we show that the F2H assays enable side-by-side analysis of substances' dual Mdm2-Mdm4 activity. In addition, they are suitable for testing various types of compounds (e.g., small molecules and peptidic inhibitors) and concurrently provide initial data on cellular toxicity. Furthermore, F2H assays readily allow real-time visualization of PPI dynamics in living cells.

Controlling translation elongation efficiency: tRNA regulation of ribosome flux on the mRNA.

Gene expression can be regulated by a wide variety of mechanisms. One example concerns the growing body of evidence that the protein-production rate can be regulated at the level of translation elongation by controlling ribosome flux across the mRNA. Variations in the abundance of tRNA molecules cause different rates of translation of their counterpart codons. This, in turn, produces a variable landscape of translational rate across each and every mRNA, with the dynamic formation and deformation of ribosomal queues being regulated by both tRNA availability and the rates of translation initiation and termination. In the present article, a range of examples of tRNA control of gene expression are reviewed, and the use of mathematical modelling to develop a predictive understanding of the consequences of that regulation is discussed and explained. These findings encourage a view that predicting the protein-synthesis rate of each mRNA requires a holistic understanding of how each stage of translation, including elongation, contributes to the overall protein-production rate.

Synthesis of substituted indenones and indanones by a Suzuki-Miyaura coupling/acid-promoted cyclisation sequence.

A one-pot Suzuki-Miyaura cross-coupling/acid-catalyzed cyclisation leading to indenones and indanones in modest to good yields is reported.

Cobalt-catalyzed diastereoselective synthesis of C-furanosides. Total synthesis of (-)-isoaltholactone.

An array of C-aryl and C-vinyl furanosides were prepared in good yields and diastereoselectivities from C-halogeno furanosides either with aryl Grignard or with vinyl Grignard using the convenient Co(acac)3/TMEDA catalytic system. This method is illustrated by the total synthesis of the (-)-isoaltholactone.

Tandem Suzuki-Miyaura coupling/acid-catalyzed cyclization between vinyl ether boronates and vinyl halides: a concise approach to polysubstituted furans.

Polysubstituted 2-(ω-hydroxyalkyl)furans were prepared by tandem Suzuki-Miyaura coupling/acid-catalyzed cyclization starting from appropriately substituted 3-haloallylic alcohols and dihydrofuran-, dihydropyran- or glycal-derived pinacol boronates.

Ribosome traffic on mRNAs maps to gene ontology: genome-wide quantification of translation initiation rates and polysome size regulation.

To understand the complex relationship governing transcript abundance and the level of the encoded protein, we integrate genome-wide experimental data of ribosomal density on mRNAs with a novel stochastic model describing ribosome traffic dynamics during translation elongation. This analysis reveals that codon arrangement, rather than simply codon bias, has a key role in determining translational efficiency. It also reveals that translation output is governed both by initiation efficiency and elongation dynamics. By integrating genome-wide experimental data sets with simulation of ribosome traffic on all Saccharomyces cerevisiae ORFs, mRNA-specific translation initiation rates are for the first time estimated across the entire transcriptome. Our analysis identifies different classes of mRNAs characterised by their initiation rates, their ribosome traffic dynamics, and by their response to ribosome availability. Strikingly, this classification based on translational dynamics maps onto key gene ontological classifications, revealing evolutionary optimisation of translation responses to be strongly influenced by gene function.

A yeast tRNA mutant that causes pseudohyphal growth exhibits reduced rates of CAG codon translation.

In Saccharomyces cerevisiae, the SUP70 gene encodes the CAG-decoding tRNA(Gln)(CUG). A mutant allele, sup70-65, induces pseudohyphal growth on rich medium, an inappropriate nitrogen starvation response. This mutant tRNA is also a UAG nonsense suppressor via first base wobble. To investigate the basis of the pseudohyphal phenotype, 10 novel sup70 UAG suppressor alleles were identified, defining positions in the tRNA(Gln)(CUG) anticodon stem that restrict first base wobble. However, none conferred pseudohyphal growth, showing altered CUG anticodon presentation cannot itself induce pseudohyphal growth. Northern blot analysis revealed the sup70-65 tRNA(Gln)(CUG) is unstable, inefficiently charged, and 80% reduced in its effective concentration. A stochastic model simulation of translation predicted compromised expression of CAG-rich ORFs in the tRNA(Gln)(CUG)-depleted sup70-65 mutant. This prediction was validated by demonstrating that luciferase expression in the mutant was 60% reduced by introducing multiple tandem CAG (but not CAA) codons into this ORF. In addition, the sup70-65 pseudohyphal phenotype was partly complemented by overexpressing CAA-decoding tRNA(Gln)(UUG), an inefficient wobble-decoder of CAG. We thus show that introducing codons decoded by a rare tRNA near the 5' end of an ORF can reduce eukaryote translational expression, and that the mutant tRNA(CUG)(Gln) constitutive pseudohyphal differentiation phenotype correlates strongly with reduced CAG decoding efficiency.

Regulation of release factor expression using a translational negative feedback loop: a systems analysis.

The essential eukaryote release factor eRF1, encoded by the yeast SUP45 gene, recognizes stop codons during ribosomal translation. SUP45 nonsense alleles are, however, viable due to the establishment of feedback-regulated readthrough of the premature termination codon; reductions in full-length eRF1 promote tRNA-mediated stop codon readthrough, which, in turn, drives partial production of full-length eRF1. A deterministic mathematical model of this eRF1 feedback loop was developed using a staged increase in model complexity. Model predictions matched the experimental observation that strains carrying the mutant SUQ5 tRNA (a weak UAA suppressor) in combination with any of the tested sup45(UAA) nonsense alleles exhibit threefold more stop codon readthrough than that of an SUQ5 yeast strain. The model also successfully predicted that eRF1 feedback control in an SUQ5 sup45(UAA) mutant would resist, but not completely prevent, imposed changes in eRF1 expression. In these experiments, the introduction of a plasmid-borne SUQ5 copy into a sup45(UAA) SUQ5 mutant directed additional readthrough and full-length eRF1 expression, despite feedback. Secondly, induction of additional sup45(UAA) mRNA expression in a sup45(UAA) SUQ5 strain also directed increased full-length eRF1 expression. The autogenous sup45 control mechanism therefore acts not to precisely control eRF1 expression, but rather as a damping mechanism that only partially resists changes in release factor expression level. The validated model predicts that the degree of feedback damping (i.e., control precision) is proportional to eRF1 affinity for the premature stop codon. The validated model represents an important tool to analyze this and other translational negative feedback loops.

Diastereoselective metal-catalyzed synthesis of C-aryl and C-vinyl glycosides.

Cobalt, the catalyst of choice: The diastereoselective cobalt-catalyzed cross-coupling of 1-bromo glycosides and aryl or vinyl Grignard reagents is described. A convenient and inexpensive catalyst, [Co(acac)(3)]/tmeda (acac = acetylacetonate, tmeda = N,N'-tetramethylethylenediamine), gives full α selectivity in the mannose and galactose series, and an α selectivity in the glucose series with α/ß ratios of 1.3:1-3:1.

Two-step one-pot synthesis of benzoannulated spiroacetals by Suzuki-Miyaura coupling/acid-catalyzed spiroacetalization.

Substituted benzoannulated spiroacetals were prepared from (2-haloaryl)alkyl alcohols and dihydropyranyl or dihydrofuranyl pinacol boronates using a Suzuki-Miyaura coupling followed by an acid-catalyzed spirocyclization. Application of the reaction to a glycal boronate provides an approach to annulated spiroacetals in enantiopure form.

Self-directed student research through analysis of microarray datasets: a computer-based functional genomics practical class for masters-level students.

This report describes a linked series of Masters-level computer practical workshops. They comprise an advanced functional genomics investigation, based upon analysis of a microarray dataset probing yeast DNA damage responses. The workshops require the students to analyse highly complex transcriptomics datasets, and were designed to stimulate active learning through experience of current research methods in bioinformatics and functional genomics. They seek to closely mimic a realistic research environment, and require the students first to propose research hypotheses, then test those hypotheses using specific sections of the microarray dataset. The complexity of the microarray data provides students with the freedom to propose their own unique hypotheses, tested using appropriate sections of the microarray data. This research latitude was highly regarded by students and is a strength of this practical. In addition, the focus on DNA damage by radiation and mutagenic chemicals allows them to place their results in a human medical context, and successfully sparks broad interest in the subject material. In evaluation, 79% of students scored the practical workshops on a five-point scale as 4 or 5 (totally effective) for student learning. More broadly, the general use of microarray data as a "student research playground" is also discussed.

Analysing GCN4 translational control in yeast by stochastic chemical kinetics modelling and simulation.

BACKGROUND: The yeast Saccharomyces cerevisiae responds to amino acid starvation by inducing the transcription factor Gcn4. This is mainly mediated via a translational control mechanism dependent upon the translation initiation eIF2·GTP·Met-tRNAiMet ternary complex, and the four short upstream open reading frames (uORFs) in its 5' mRNA leader. These uORFs act to attenuate GCN4 mRNA translation under normal conditions. During amino acid starvation, levels of ternary complex are reduced. This overcomes the GCN4 translation attenuation effect via a scanning/reinitiation control mechanism dependent upon uORF spacing. RESULTS: Using published experimental data, we have developed and validated a probabilistic formulation of GCN4 translation using the Chemical Master Equation (Model 1). Model 1 explains GCN4 translation's nonlinear dependency upon uORF placements, and predicts that an as yet unidentified factor, which was proposed to regulate GCN4 translation under some conditions, only has pronounced effects upon GCN4 translation when intercistronic distances are unnaturally short. A simpler Model 2 that does not include this unidentified factor could well represent the regulation of a natural GCN4 mRNA. Using parameter values optimised for this algebraic Model 2, we performed stochastic simulations by Gillespie algorithm to investigate the distribution of ribosomes in different sections of GCN4 mRNA under distinct conditions. Our simulations demonstrated that ribosomal loading in the 5'-untranslated region is mainly determined by the ratio between the rates of 5'-initiation and ribosome scanning, but was not significantly affected by rate of ternary complex binding. Importantly, the translation rate for codons starved of cognate tRNAs is predicted to be the most significant contributor to the changes in ribosomal loading in the coding region under repressing and derepressing conditions. CONCLUSIONS: Our integrated probabilistic Models 1 and 2 explained GCN4 translation and helped to elucidate the role of a yet unidentified factor. The ensuing stochastic simulations evaluated different factors that may impact on the translation of GCN4 mRNA, and integrated translation status with ribosomal density.