A Rho family GTPase controls actin dynamics and tip growth via two counteracting downstream pathways in pollen tubes

Tip growth in neuronal cells, plant cells, and fungal hyphae is known to require tip-localized Rho GTPase, calcium, and filamentous actin (F-actin), but how they interact with each other is unclear. The pollen tube is an exciting model to study spatiotemporal regulation of tip growth and F-actin dynamics. An Arabidopsis thaliana Rho family GTPase, ROP1, controls pollen tube growth by regulating apical F-actin dynamics. This paper shows that ROP1 activates two counteracting pathways involving the direct targets of tip-localized ROP1: RIC3 and RIC4. RIC4 promotes F-actin assembly, whereas RIC3 activates Ca(2+) signaling that leads to F-actin disassembly. Overproduction or depletion of either RIC4 or RIC3 causes tip growth defects that are rescued by overproduction or depletion of RIC3 or RIC4, respectively. Thus, ROP1 controls actin dynamics and tip growth through a check and balance between the two pathways. The dual and antagonistic roles of this GTPase may provide a unifying mechanism by which Rho modulates various processes dependent on actin dynamics in eukaryotic cells.     

Relative contribution of cell contact pattern, specific PKC isoforms and gap junctional communication in tight junction assembly in the mouse early embryo

In mouse early development, cell contact patterns regulate the spatial organization and segregation of inner cell mass (ICM) and trophectoderm epithelium (TE) during blastocyst morphogenesis. Progressive membrane assembly of tight junctional (TJ) proteins in the differentiating TE during cleavage is upregulated by cell contact asymmetry (outside position) and suppressed within the ICM by cell contact symmetry (inside position). This is reversible, and immunosurgical isolation of the ICM induces upregulation of TJ assembly in a sequence that broadly mimics that occurring during blastocyst formation. The mechanism relating cell contact pattern and TJ assembly was investigated in the ICM model with respect to PKC-mediated signaling and gap junctional communication. Our results indicate that complete cell contact asymmetry is required for TJ biogenesis and acts upstream of PKC-mediated signaling. Specific inhibition of two PKC isoforms, PKCdelta and zeta, revealed that both PKC activities are required for membrane assembly of ZO-2 TJ protein, while only PKCzeta activity is involved in regulating ZO-1alpha+ membrane assembly, suggesting different mechanisms for individual TJ proteins. Gap junctional communication had no apparent influence on either TJ formation or PKC signaling but was itself affected by changes of cell contact patterns. Our data suggest that the dynamics of cell contact patterns coordinate the spatial organization of TJ formation via specific PKC signaling pathways during blastocyst biogenesis.     

Positive feedback interactions between microtubule and actin dynamics during cell motility

The migration of tissue cells requires interplay between the microtubule and actin cytoskeletal systems. Recent reports suggest that interactions of microtubules with actin dynamics creates a polarization of microtubule assembly behavior in cells, such that microtubule growth occurs at the leading edge and microtubule shortening occurs at the cell body and rear. Microtubule growth and shortening may activate Rac1 and RhoA signaling, respectively, to control actin dynamics. Thus, an actin-dependent gradient in microtubule dynamic-instability parameters in cells may feed back through the activation of specific signalling pathways to perpetuate the polarized actin-assembly dynamics required for cell motility.     

Master equation approach to the assembly of viral capsids

The distribution of inequivalent geometries occurring during self-assembly of the major capsid protein in thermodynamic equilibrium is determined based on a master equation approach. These results are implemented to characterize the assembly of SV40 virus and to obtain information on the putative pathways controlling the progressive build-up of the SV40 capsid. The experimental testability of the predictions is assessed and an analysis of the geometries of the assembly intermediates on the dominant pathways is used to identify targets for anti-viral drug design.     

SNARE assembly and disassembly exhibit a pronounced hysteresis.

SNARE proteins are essential for intracellular membrane fusion of eukaryotes. Their assembly into stable four-helix bundles bridges membranes and may provide the energy for initiating membrane fusion. In vitro, assembly of soluble SNARE fragments is accompanied by major structural rearrangements that can be described as a folding reaction. The pathways and the thermodynamics of SNARE protein interactions, however, are not known. Here we report that assembly and dissociation of two distantly related SNARE complexes exhibit a marked hysteresis. The assembled and disassembled native states are separated by a kinetic barrier and cannot equilibrate on biologically relevant timescales. We suggest that the hysteresis is a hallmark of all SNARE complexes and that complex assembly and disassembly follow different pathways that may be independently controlled.     

Stochastic dynamics of the nerve growth cone and its microtubules during neurite outgrowth

The controlled extension of neurites is essential not only for nervous system development, but also for effective nerve regeneration after injury. This process is critically dependent on microtubule assembly since axons fail to elongate in the presence of drugs which disrupt normal assembly dynamics. For this reason, neurite outgrowth is potentially controllable by manipulation of the assembly state of the intracellular array of microtubules. Therefore, understanding how microtubule assembly dynamics and neurite outgrowth are coupled, in the absence of drugs, can lend valuable insight into the control and guidance of the outgrowth process. In the present study we characterized the stochastic dynamics of neurite outgrowth and its corresponding microtubule array, which advances concomitantly with the advance of the nerve growth cone, the highly motile structure at the terminus of the growing neurite, using reported fluorescent microscopic image sequences (Tanaka and Kirschner, 1991, J. Cell Biol. 115:345-363). Although previously modeled as an uncorrelated random walk, the stochastic advance of the growth cone was found to be anticorrelated over a time scale of approximately 4 min, meaning that growth cone advances tended to be followed by growth cone retractions approximately 4 min later. The observed anticorrelation most likely reflects the periodic stops and starts of neurite outgrowth that have been reported anecdotally. A strikingly similar pattern of anticorrelation was also identified in the advance of the growth cone's microtubule array. Cross-correlation analysis showed that growth cone dynamics tended to precede microtubule dynamics on a time scale of approximately 0-2 min, while microtubules tended to precede growth cone dynamics on a approximately 0-20-s time scale, indicating a close temporal coupling between microtubule and growth cone dynamics. Finally, the scaling of the mean-squared displacements with time for both the growth cone and microtubules suggested a fractional Brownian motion model which accounts for the observed anticorrelation of growth cone and microtubule advance. (c) 1996 John Wiley & Sons, Inc.     

Distinct ECM mechanosensing pathways regulate microtubule dynamics to control endothelial cell branching morphogenesis

During angiogenesis, cytoskeletal dynamics that mediate endothelial cell branching morphogenesis during vascular guidance are thought to be regulated by physical attributes of the extracellular matrix (ECM) in a process termed mechanosensing. Here, we tested the involvement of microtubules in linking mechanosensing to endothelial cell branching morphogenesis. We used a recently developed microtubule plus end-tracking program to show that specific parameters of microtubule assembly dynamics, growth speed and growth persistence, are globally and regionally modified by, and contribute to, ECM mechanosensing. We demonstrated that engagement of compliant two-dimensional or three-dimensional ECMs induces local differences in microtubule growth speed that require myosin II contractility. Finally, we found that microtubule growth persistence is modulated by myosin II-mediated compliance mechanosensing when cells are cultured on two-dimensional ECMs, whereas three-dimensional ECM engagement makes microtubule growth persistence insensitive to changes in ECM compliance. Thus, compliance and dimensionality ECM mechanosensing pathways independently regulate specific and distinct microtubule dynamics parameters in endothelial cells to guide branching morphogenesis in physically complex ECMs.     

Energetics of quasiequivalence: computational analysis of protein-protein interactions in icosahedral viruses.

Quaternary structure polymorphism found in quasiequivalent virus capsids provides a static framework for studying the dynamics of protein interactions. The same protein subunits are found in different structural environments within these particles, and in some cases, the molecular switching required for the polymorphic quaternary interactions is obvious from high-resolution crystallographic studies. Employing atomic resolution structures, molecular mechanics, and continuum electrostatic methods, we have computed association energies for unique subunit interfaces of three icosahedral viruses, black beetle virus, southern bean virus, and human rhinovirus 14. To quantify the chemical determinants of quasiequivalence, the energetic contributions of individual residues forming quasiequivalent interfaces were calculated and compared. The potential significance of the differences in stabilities at quasiequivalent interfaces was then explored with the combinatorial assembly approach. The analysis shows that the unique association energies computed for each virus serve as a sensitive basis set that may determine distinct intermediates and pathways of virus capsid assembly. The pathways for the quasiequivalent viruses displayed isoenergetic oligomers at specific points, suggesting that these may determine the quaternary structure polymorphism required for the assembly of a quasiequivalent particle.     

An Assembly Funnel Makes Biomolecular Complex Assembly Efficient

Like protein folding and crystallization, the self-assembly of complexes is a fundamental form of biomolecular organization. While the number of methods for creating synthetic complexes is growing rapidly, most require empirical tuning of assembly conditions and/or produce low yields. We use coarse-grained simulations of the assembly kinetics of complexes to identify generic limitations on yields that arise because of the many simultaneous interactions allowed between the components and intermediates of a complex. Efficient assembly occurs when nucleation is fast and growth pathways are few, i.e. when there is an assembly “funnel”. For typical complexes, an assembly funnel occurs in a narrow window of conditions whose location is highly complex specific. However, by redesigning the components this window can be drastically broadened, so that complexes can form quickly across many conditions. The generality of this approach suggests assembly funnel design as a foundational strategy for robust biomolecular complex synthesis.     

NERDSS: A Nonequilibrium Simulator for Multibody Self-Assembly at the Cellular Scale

Currently, a significant barrier to building predictive models of cellular self-assembly processes is that molecular models cannot capture minutes-long dynamics that couple distinct components with active processes, whereas reaction-diffusion models cannot capture structures of molecular assembly. Here, we introduce the nonequilibrium reaction-diffusion self-assembly simulator (NERDSS), which addresses this spatiotemporal resolution gap. NERDSS integrates efficient reaction-diffusion algorithms into generalized software that operates on user-defined molecules through diffusion, binding and orientation, unbinding, chemical transformations, and spatial localization. By connecting the fast processes of binding with the slow timescales of large-scale assembly, NERDSS integrates molecular resolution with reversible formation of ordered, multisubunit complexes. NERDSS encodes models using rule-based formatting languages to facilitate model portability, usability, and reproducibility. Applying NERDSS to steps in clathrin-mediated endocytosis, we design multicomponent systems that can form lattices in solution or on the membrane, and we predict how stochastic but localized dephosphorylation of membrane lipids can drive lattice disassembly. The NERDSS simulations reveal the spatial constraints on lattice growth and the role of membrane localization and cooperativity in nucleating assembly. By modeling viral lattice assembly and recapitulating oscillations in protein expression levels for a circadian clock model, we illustrate the adaptability of NERDSS. NERDSS simulates user-defined assembly models that were previously inaccessible to existing software tools, with broad applications to predicting self-assembly in vivo and designing high-yield assemblies in vitro.     

Kinetically competing huntingtin aggregation pathways control amyloid polymorphism and properties

In polyglutamine (polyQ) containing fragments of the Huntington's disease protein huntingtin (htt), the N-terminal 17 amino acid htt(NT) segment serves as the core of α-helical oligomers whose reversible assembly locally concentrates the polyQ segments, thereby facilitating polyQ amyloid nucleation. A variety of aggregation inhibitors have been described that achieve their effects by neutralizing this concentrating function of the htt(NT) segment. In this paper we characterize the nature and limits of this inhibition for three means of suppressing htt(NT)-mediated aggregation. We show that the previously described action of htt(NT) peptide-based inhibitors is solely due to their ability to suppress the htt(NT)-mediated aggregation pathway. That is, under htt(NT) inhibition, nucleation of polyQ amyloid formation by a previously described alternative nucleation mechanism proceeds unabated and transiently dominates the aggregation process. Removal of the bulk of the htt(NT) segment by proteolysis or mutagenesis also blocks the htt(NT)-mediated pathway, allowing the alternative nucleation pathway to dominate. In contrast, the previously described immunoglobulin-based inhibitor, the antihtt(NT) V(L) 12.3 protein, effectively blocks both amyloid pathways, leading to stable accumulation of nonamyloid oligomers. These data show that the htt(NT)-dependent and -independent pathways of amyloid nucleation in polyQ-containing htt fragments are in direct kinetic competition. The results illustrate how amyloid polymorphism depends on assembly mechanism and kinetics and have implications for how the intracellular environment can influence aggregation pathways.     

Simulation study of the contribution of oligomer/oligomer binding to capsid assembly kinetics

The process by which hundreds of identical capsid proteins self-assemble into icosahedral structures is complex and poorly understood. Establishing constraints on the assembly pathways is crucial to building reliable theoretical models. For example, it is currently an open question to what degree overall assembly kinetics are dominated by one or a few most efficient pathways versus the enormous number theoretically possible. The importance of this question, however, is often overlooked due to the difficulties of addressing it in either theoretical or experimental practice. We apply a computer model based on a discrete-event simulation method to evaluate the contributions of nondominant pathways to overall assembly kinetics. This is accomplished by comparing two possible assembly models: one allowing growth to proceed only by the accretion of individual assembly subunits and the other allowing the binding of sterically compatible assembly intermediates any sizes. Simulations show that the two models perform almost identically under low binding rate conditions, where growth is strongly nucleation-limited, but sharply diverge under conditions of higher association rates or coat protein concentrations. The results suggest the importance of identifying the actual binding pattern if one is to build reliable models of capsid assembly or other complex self-assembly processes.     

The function of actin-binding proteins in pollen tube growth

Pollen tube growth is a key step in sexual reproduction of higher plants. The pollen tube is a typical example of tip-growing cells and shows a polarized cytoplasm. To develop and maintain polarized growth, pollen tubes need a carefully regulated actin cytoskeleton. It is well known that actin-binding proteins are responsible for the direct control of dynamic actin filaments and serve as a link between signal transduction pathways and dynamic actin changes in determining cellular architecture. Several of these classes have been identified in pollen tubes and their detailed characterisation is progressing rapidly. Here, we aim to survey what is known about the major actin-binding proteins that affect actin assembly and dynamics, and their higher-order organisation in pollen tube growth.     

Exploring the parameter space of complex self-assembly through virus capsid models

We use discrete event stochastic simulations to characterize the parameter space of a model of icosahedral viral capsid assembly as functions of monomer-monomer binding rates. The simulations reveal a parameter space characterized by three major assembly mechanisms, a standard nucleation-limited monomer-accretion pathway and two distinct hierarchical assembly pathways, as well as unproductive regions characterized by kinetically trapped species. Much of the productive parameter space also consists of border regions between these domains where hybrid pathways are likely to operate. A simpler octamer system studied for comparison reveals three analogous pathways, but is characterized by much lesser sensitivity to parameter variations in contrast to the sharp changes visible in the icosahedral model. The model suggests that modest changes in assembly conditions, consistent with expected differences between in vitro and in vivo assembly environments, could produce substantial shifts in assembly pathways. These results suggest that we must be cautious in drawing conclusions about in vivo capsid self-assembly dynamics from theoretical or in vitro models, as the nature of the basic assembly mechanisms accessible to a system can substantially differ between simple and complex model systems, between theoretical models and simulation results, and between in vitro and in vivo assembly conditions.     

Visualization and Analysis of Microtubule Dynamics Using Dual Color-Coded Display of Plus-End Labels

Investigating spatial and temporal control of microtubule dynamics in live cells is critical to understanding cell morphogenesis in development and disease. Tracking fluorescently labeled plus-end-tracking proteins over time has become a widely used method to study microtubule assembly. Here, we report a complementary approach that uses only two images of these labels to visualize and analyze microtubule dynamics at any given time. Using a simple color-coding scheme, labeled plus-ends from two sequential images are pseudocolored with different colors and then merged to display color-coded ends. Based on object recognition algorithms, these colored ends can be identified and segregated into dynamic groups corresponding to four events, including growth, rescue, catastrophe, and pause. Further analysis yields not only their spatial distribution throughout the cell but also provides measurements such as growth rate and direction for each labeled end. We have validated the method by comparing our results with ground-truth data derived from manual analysis as well as with data obtained using the tracking method. In addition, we have confirmed color-coded representation of different dynamic events by analyzing their history and fate. Finally, we have demonstrated the use of the method to investigate microtubule assembly in cells and provided guidance in selecting optimal image acquisition conditions. Thus, this simple computer vision method offers a unique and quantitative approach to study spatial regulation of microtubule dynamics in cells.     

Species assembly in model ecosystems, I: Analysis of the population model and the invasion dynamics

Recently we have introduced a simplified model of ecosystem assembly (Capitán et al., 2009) for which we are able to map out all assembly pathways generated by external invasions in an exact manner. In this paper we provide a deeper analysis of the model, obtaining analytical results and introducing some approximations which allow us to reconstruct the results of our previous work. In particular, we show that the population dynamics equations of a very general class of trophic-level structured food-web have an unique interior equilibrium point which is globally stable. We show analytically that communities found as end states of the assembly process are pyramidal and we find that the equilibrium abundance of any species at any trophic level is approximately inversely proportional to the number of species in that level. We also find that the per capita growth rate of a top predator invading a resident community is key to understand the appearance of complex end states reported in our previous work. The sign of these rates allows us to separate regions in the space of parameters where the end state is either a single community or a complex set containing more than one community. We have also built up analytical approximations to the time evolution of species abundances that allow us to determine, with high accuracy, the sequence of extinctions that an invasion may cause. Finally we apply this analysis to obtain the communities in the end states. To test the accuracy of the transition probability matrix generated by this analytical procedure for the end states, we have compared averages over those sets with those obtained from the graph derived by numerical integration of the Lotka-Volterra equations. The agreement is excellent.