Table of contents

Amyloid fibrils characterized as highly intractable thread-like species are associated with many neurodegenerative diseases. Although neither the mechanism of amyloid formation nor the origin of amyloid toxicity is currently completely understood, the detailed three-dimensional atomic structures of the yeast protein Sup35 and Aβ amyloid protein determined by recent experiments provide the first and important step towards the comprehension of the pathogenesis and aggregation mechanisms of amyloid diseases. By analyzing these two amyloid peptides which have available crystal structures and other amyloid sequences with proposed structures using computational simulations, we delineate three common features in amyloid organizations and amyloid structures. These could contribute to an improved understanding of the molecular mechanism of amyloid formation, the nature of the aggregation driving forces that stabilize these structures and the development of potential therapeutic agents against amyloid diseases.

Computer simulations of β-hairpin folding are relatively difficult, especially those based on the explicit water model. This greatly limits the complete analysis and understanding of their folding mechanisms. In this paper, we use the generalized Born/solvent accessible implicit solvent model to simulate the folding processes of a nine-residue β-hairpin. We find that the β-hairpin can fold into its native structure very easily, even using the traditional molecular dynamics method. This allows us to extract 21 complete folding events and investigate the folding process sufficiently. Our results show that there exist four most stable states on the free energy landscape of the short peptide, one native state and three intermediates. We find that two of the non-native stable states have almost the same potential energy as the native state but with lower entropy. This suggests that the native state can be stabilized entropically. Furthermore, we find that the folding processes of this peptide have common features: to fold into its native state, the peptide undergoes a continuous collapsing–extending–recollapsing process to adjust the positions of the side chains in order to form the native middle inter-strand hydrogen bonds. The formations of these bonds are the key step of the folding process. Once these bonds are formed, the peptide can fold into the native state quickly.

A population of cells exhibits wide phenotypic variation even if it is genetically homogeneous. In particular, individual cells differ from one another in the amount of protein they express under a given regulatory system under fixed conditions. Here we study how protein distributions in a population of the yeast S. cerevisiae are shaped by a balance of processes: protein production—an intracellular process—and protein dilution due to cell division—a population process. We measure protein distributions by employing reporter green fluorescence protein (gfp) under the regulation of the yeast GAL system under conditions where it is metabolically essential. Cell populations are grown in chemostats, thus allowing control of the environment and stable measurements of distribution dynamics over many generations. Despite the essential functional role of the GAL system in a pure galactose medium, steady-state distributions are found to be universally broad, with exponential tails and a large standard-deviation-to-mean ratio. Under several different perturbations the dynamics of the distribution is observed to be asymmetric, with a much longer time to build a wide expression distribution from below compared with a fast relaxation of the distribution toward steady state from above. These results show that the main features of the protein distributions are largely determined by population effects and are less sensitive to the intracellular biochemical noise.

Gene arrangement into operons varies between bacterial species. Genes in a given system can be on one operon in some organisms and on several operons in other organisms. Existing theories explain why genes that work together should be on the same operon, since this allows for advantageous lateral gene transfer and accurate stoichiometry. But what causes the frequent separation into multiple operons of co-regulated genes that act together in a pathway? Here we suggest that separation is due to benefits made possible by differential regulation of each operon. We present a simple mathematical model for the optimal distribution of genes into operons based on a balance of the cost of operons and the benefit of regulation that provides 'just-when-needed' temporal order. The analysis predicts that genes are arranged such that genes on the same operon do not skip functional steps in the pathway. This prediction is supported by genomic data from 137 bacterial genomes. Our work suggests that gene arrangement is not only the result of random historical drift, genome re-arrangement and gene transfer, but has elements that are solutions of an evolutionary optimization problem. Thus gene functional order may be inferred by analyzing the operon structure across different genomes.

A key mediator of eukaryotic chemotaxis is the asymmetric accumulation of phosphatidylinositol-3,4,5-triphosphate (PIP3) on the cell membrane. Recent work has focused on understanding how a shallow external gradient of chemoattractant leads to an amplified internal gradient of PIP3. In this paper we dissect what fraction of this amplification is derived biochemically by the signal transduction network and how much arises entirely from the effects of cell morphology. Here we identify and formalize the role of morphology in signal detection and demonstrate its effects through simulation and experiments. Our key result is that an asymmetric distribution of membrane accounts for approximately one-half of the measured amplification from ligand concentration to PIP3 production. We also show that the underlying biochemical network behaves as a linear amplifier in the micropipette assay.

We consider a stochastic model of transcription factor (TF)-regulated gene expression. The model describes two genes, gene A and gene B, which synthesize the TFs and the target gene proteins, respectively. We show through analytic calculations that the TF fluctuations have a significant effect on the distribution of the target gene protein levels when the mean TF level falls in the highest sensitive region of the dose–response curve. We further study the effect of reducing the copy number of gene A from two to one. The enhanced TF fluctuations yield results different from those in the deterministic case. The probability that the target gene protein level exceeds a threshold value is calculated with the knowledge of the probability density functions associated with the TF and target gene protein levels. Numerical simulation results for a more detailed stochastic model are shown to be in agreement with those obtained through analytic calculations. The relevance of these results in the context of the genetic disorder haploinsufficiency is pointed out. Some experimental observations on the haploinsufficiency of the tumour suppressor gene, Nkx 3.1, are explained with the help of the stochastic model of TF-regulated gene expression.

Directed secretion of cytotoxins or cytokines by T cells during immune response depends on migration of the centrosome in the T cell to the interface with the target cell. The mechanism of the centrosome translocation has been elusive. The presented computational analysis demonstrates that the centrosome should be positioned at the interface if the T cell attempts simultaneously (a) to minimize its surface area, (b) to maximize the interface area, (c) to maintain the cell volume and (d) to straighten the microtubules. Live three-dimensional microscopy and measurements show that the optimal position of the centrosome is achieved in large part (by about 40%) via rolling of the entire T cell body on the target surface; this movement appears to entrain the centrosome. The theoretical and experimental results draw attention to the previously unrecognized role of the whole-cell structure and whole-cell movements in the T cell polarization.

Behavioral researchers have attached magnets to birds during orientation experiments, assuming that such magnets will disrupt their ability to obtain magnetic information. Here, we investigate the effect of an attached magnet on the ability to derive directional information from a radical-pair based compass mechanism. We outline in some detail the geometrical symmetries that would allow a bird to identify magnetic directions in a radical-pair based compass. We show that the artificial field through an attached magnet will quickly disrupt the birds' ability to distinguish pole-ward from equator-ward headings, but that much stronger fields are necessary to disrupt their ability to detect the magnetic axis. Together with estimates of the functional limits of a radical-pair based compass, our calculations suggest that artificial fields of comparable size to the geomagnetic field are not generally sufficient to render a radical-pair based compass non-functional.