Changing Cellular Location of CheZ Predicted by Molecular Simulations

In the chemotaxis pathway of the bacterium Escherichia coli, signals are carried from a cluster of receptors to the flagellar motors by the diffusion of the protein CheY-phosphate (CheYp) through the cytoplasm. A second protein, CheZ, which promotes dephosphorylation of CheYp, partially colocalizes with receptors in the plasma membrane. CheZ is normally dimeric in solution but has been suggested to associate into highly active oligomers in the presence of CheYp. A model is presented here and supported by Brownian dynamics simulations, which accounts for these and other experimental data: A minority component of the receptor cluster (dimers of CheA(short)) nucleates CheZ oligomerization and CheZ molecules move from the cytoplasm to a bound state at the receptor cluster depending on the current level of cellular stimulation. The corresponding simulations suggest that dynamic CheZ localization will sharpen cellular responses to chemoeffectors, increase the range of detectable ligand concentrations, and make adaptation more precise and robust. The localization and activation of CheZ constitute a negative feedback loop that provides a second tier of adaptation to the system. Subtle adjustments of this kind are likely to be found in many other signaling pathways

Modelling the Establishment of PAR Protein Polarity in the One-Cell C. elegans Embryo

At the one-cell stage, the C. elegans embryo becomes polarized along the anterior-posterior axis. The PAR proteins form complementary anterior and posterior domains in a dynamic process driven by cytoskeletal rearrangement. Initially, the PAR proteins are uniformly distributed throughout the embryo. Following a cue from fertilization, cortical actomyosin contracts towards the anterior pole. PAR-3/PAR-6/PKC-3 (the anterior PAR proteins) become restricted to the anterior cortex. PAR-1 and PAR-2 (the posterior PAR proteins) become enriched in the posterior cortical region. We present a mathematical model of this polarity establishment process, in which we take a novel approach to combine reaction-diffusion dynamics of the PAR proteins coupled to a simple model of actomyosin contraction. We show that known interactions between the PAR proteins are sufficient to explain many aspects of the observed cortical PAR dynamics in both wild-type and mutant embryos. However, cytoplasmic PAR protein polarity, which is vital for generating daughter cells with distinct molecular components, cannot be properly explained within such a framework. We therefore consider additional mechanisms that can reproduce the proper cytoplasmic polarity. In particular we predict that cytoskeletal asymmetry in the cytoplasm, in addition to the cortical actomyosin asymmetry, is a critical determinant of PAR protein localization.Comment: 28 pages, 8 figure

Overview of mathematical approaches used to model bacterial chemotaxis I: the single cell

Mathematical modeling of bacterial chemotaxis systems has been influential and insightful in helping to understand experimental observations. We provide here a comprehensive overview of the range of mathematical approaches used for modeling, within a single bacterium, chemotactic processes caused by changes to external gradients in its environment. Specific areas of the bacterial system which have been studied and modeled are discussed in detail, including the modeling of adaptation in response to attractant gradients, the intracellular phosphorylation cascade, membrane receptor clustering, and spatial modeling of intracellular protein signal transduction. The importance of producing robust models that address adaptation, gain, and sensitivity are also discussed. This review highlights that while mathematical modeling has aided in understanding bacterial chemotaxis on the individual cell scale and guiding experimental design, no single model succeeds in robustly describing all of the basic elements of the cell. We conclude by discussing the importance of this and the future of modeling in this area

The proteins encoded by the pogo-like Lemi1 element bind the TIRs and subterminal repeated motifs of the Arabidopsis Emigrant MITE: consequences for the transposition mechanism of MITEs

MITEs (miniature inverted-repeated transposable elements) are a particular class of defective DNA transposons usually present within genomes as high copy number populations of highly homogeneous elements. Although an active MITE, the mPing element, has recently been characterized in rice, the transposition mechanism of MITEs remains unknown. It has been proposed that transposases of related transposons could mobilize MITEs in trans. Moreover, it has also been proposed that the presence of conserved terminal inverted-repeated (TIR) sequences could be the only requirement of MITEs for mobilization, allowing divergent or unrelated elements to be mobilized by a particular transposase. We present here evidence for a recent mobility of the Arabidopsis Emigrant MITE and we report on the capacity of the proteins encoded by the related Lemi1 transposon, a pogo-related element, to specifically bind Emigrant elements. This suggests that Lemi1 could mobilize Emigrant elements and makes the Lemi1/Emigrant couple an ideal system to study the transposition mechanism of MITEs. Our results show that Lemi1 proteins bind Emigrant TIRs but also bind cooperatively to subterminal repeated motifs. The requirement of internal sequences for the formation of proper DNA/protein structure could affect the capacity of divergent MITEs to be mobilized by distantly related transposases

Identification of multiple binding sites for the THAP domain of the Galileo transposase in the long terminal inverted-repeats

Galileo is a DNA transposon responsible for the generation of several chromosomal inversions in Drosophila. In contrast to other members of the P-element superfamily, it has unusually long terminal inverted-repeats (TIRs) that resemble those of Foldback elements. To investigate the function of the long TIRs we derived consensus and ancestral sequences for the Galileo transposase in three species of Drosophilids. Following gene synthesis, we expressed and purified their constituent THAP domains and tested their binding activity towards the respective Galileo TIRs. DNase I footprinting located the most proximal DNA binding site about 70 bp from the transposon end. Using this sequence we identified further binding sites in the tandem repeats that are found within the long TIRs. This suggests that the synaptic complex between Galileo ends may be a complicated structure containing higher-order multimers of the transposase. We also attempted to reconstitute Galileo transposition in Drosophila embryos but no events were detected. Thus, although the limited numbers of Galileo copies in each genome were sufficient to provide functional consensus sequences for the THAP domains, they do not specify a fully active transposase. Since the THAP recognition sequence is short, and will occur many times in a large genome, it seems likely that the multiple binding sites within the long, internally repetitive, TIRs of Galileo and other Foldback-like elements may provide the transposase with its binding specificity

Transposition of the human Hsmar1 transposon: rate-limiting steps and the importance of the flanking TA dinucleotide in second strand cleavage

Hsmar1 is a member of the mariner family of DNA transposons. Although widespread in nature, their molecular mechanism remains obscure. Many other cut-and-paste elements use a hairpin intermediate to cleave the two strands of DNA at each transposon end. However, this intermediate is absent in mariner, suggesting that these elements use a fundamentally different mechanism for second-strand cleavage. We have taken advantage of the faithful and efficient in vitro reaction provided by Hsmar1 to characterize the products and intermediates of transposition. We report different factors that particularly affect the reaction, which are the reaction pH and the transposase concentration. Kinetic analysis revealed that first-strand nicking and integration are rapid. The rate of the reaction is limited in part by the divalent metal ion-dependent assembly of a complex between transposase and the transposon end(s) prior to the first catalytic step. Second-strand cleavage is the rate-limiting catalytic step of the reaction. We discuss our data in light of a model for the two metal ion catalytic mechanism and propose that mariner excision involves a significant conformational change between first- and second-strand cleavage at each transposon end. Furthermore, this conformational change requires specific contacts between transposase and the flanking TA dinucleotide

Meredys, a multi-compartment reaction-diffusion simulator using multistate realistic molecular complexes

<p>Abstract</p> <p>Background</p> <p>Most cellular signal transduction mechanisms depend on a few molecular partners whose roles depend on their position and movement in relation to the input signal. This movement can follow various rules and take place in different compartments. Additionally, the molecules can form transient complexes. Complexation and signal transduction depend on the specific states partners and complexes adopt. Several spatial simulator have been developed to date, but none are able to model reaction-diffusion of realistic multi-state transient complexes.</p> <p>Results</p> <p><it>Meredys </it>allows for the simulation of multi-component, multi-feature state molecular species in two and three dimensions. Several compartments can be defined with different diffusion and boundary properties. The software employs a Brownian dynamics engine to simulate reaction-diffusion systems at the reactive particle level, based on compartment properties, complex structure, and hydro-dynamic radii. Zeroth-, first-, and second order reactions are supported. The molecular complexes have realistic geometries. Reactive species can contain user-defined feature states which can modify reaction rates and outcome. Models are defined in a versatile NeuroML input file. The simulation volume can be split in subvolumes to speed up run-time.</p> <p>Conclusions</p> <p><it>Meredys </it>provides a powerful and versatile way to run accurate simulations of molecular and sub-cellular systems, that complement existing multi-agent simulation systems. <it>Meredys </it>is a Free Software and the source code is available at <url>http://meredys.sourceforge.net/</url>.</p

Crosstalk between transposase subunits during cleavage of the mariner transposon

Mariner transposition is a complex reaction that involves three recombination sites and six strand breaking and joining reactions. This requires precise spatial and temporal coordination between the different components to ensure a productive outcome and minimize genomic instability. We have investigated how the cleavage events are orchestrated within the mariner transpososome. We find that cleavage of the non-transferred strand is completed at both transposon ends before the transferred strand is cleaved at either end. By introducing transposon-end mutations that interfere with cleavage, but leave transpososome assembly unaffected, we demonstrate that a structural transition preceding transferred strand cleavage is coordinated between the two halves of the transpososome. Since mariner lacks the DNA hairpin intermediate, this transition probably reflects a reorganization of the transpososome to allow the access of different monomers onto the second pair of strands, or the relocation of the DNA within the same active site between two successive hydrolysis events. Communication between transposase subunits also provides a failsafe mechanism that restricts the generation of potentially deleterious double-strand breaks at isolated sites. Finally, we identify transposase mutants that reveal that the conserved WVPHEL motif provides a structural determinant of the coordination mechanism

Solution conformations of early intermediates in Mos1 transposition

DNA transposases facilitate genome rearrangements by moving DNA transposons around and between genomes by a cut-and-paste mechanism. DNA transposition proceeds in an ordered series of nucleoprotein complexes that coordinate pairing and cleavage of the transposon ends and integration of the cleaved ends at a new genomic site. Transposition is initiated by transposase recognition of the inverted repeat sequences marking each transposon end. Using a combination of solution scattering and biochemical techniques, we have determined the solution conformations and stoichiometries of DNA-free Mos1 transposase and of the transposase bound to a single transposon end. We show that Mos1 transposase is an elongated homodimer in the absence of DNA and that the N-terminal 55 residues, containing the first helix-turn-helix motif, are required for dimerization. This arrangement is remarkably different from the compact, crossed architecture of the dimer in the Mos1 paired-end complex (PEC). The transposase remains elongated when bound to a single-transposon end in a pre-cleavage complex, and the DNA is bound predominantly to one transposase monomer. We propose that a conformational change in the single-end complex, involving rotation of one half of the transposase along with binding of a second transposon end, could facilitate PEC assembly