An agent-based model for mRNA export through the nuclear pore complex.

mRNA export from the nucleus is an essential step in the expression of every protein- coding gene in eukaryotes, but many aspects of this process remain poorly understood. The density of export receptors that must bind an mRNA to ensure export, as well as how receptor distribution affects transport dynamics, is not known. It is also unclear whether the rate-limiting step for transport occurs at the nuclear basket, in the central channel, or on the cytoplasmic face of the nuclear pore complex. Using previously published biophysical and biochemical parameters of mRNA export, we implemented a three-dimensional, coarse-grained, agent-based model of mRNA export in the nanosecond regime to gain insight into these issues. On running the model, we observed that mRNA export is sensitive to the number and distribution of transport receptors coating the mRNA and that there is a rate-limiting step in the nuclear basket that is potentially associated with the mRNA reconfiguring itself to thread into the central channel. Of note, our results also suggest that using a single location-monitoring mRNA label may be insufficient to correctly capture the time regime of mRNA threading through the pore and subsequent transport. This has implications for future experimental design to study mRNA transport dynamics

Role of KASH domain lengths in the regulation of LINC complexes.

The linker of the nucleoskeleton and cytoskeleton (LINC) complex is formed by the conserved interactions between Sad-1 and UNC-84 (SUN) and Klarsicht, ANC-1, SYNE homology (KASH) domain proteins, providing a physical coupling between the nucleoskeleton and cytoskeleton that mediates the transfer of physical forces across the nuclear envelope. The LINC complex can perform distinct cellular functions by pairing various KASH domain proteins with the same SUN domain protein. For example, in Caenorhabditis elegans, SUN protein UNC-84 binds to two KASH proteins UNC-83 and ANC-1 to mediate nuclear migration and anchorage, respectively. In addition to distinct cytoplasmic domains, the luminal KASH domain also varies among KASH domain proteins of distinct functions. In this study, we combined in vivo C. elegans genetics and in silico molecular dynamics simulations to understand the relation between the length and amino acid composition of the luminal KASH domain, and the function of the SUN-KASH complex. We show that longer KASH domains can withstand and transfer higher forces and interact with the membrane through a conserved membrane proximal EEDY domain that is unique to longer KASH domains. In agreement with our models, our in vivo results show that swapping the KASH domains of ANC-1 and UNC-83, or shortening the KASH domain of ANC-1, both result in a nuclear anchorage defect in C. elegans

Agent-Based Modeling in Molecular Systems Biology.

Molecular systems orchestrating the biology of the cell typically involve a complex web of interactions among various components and span a vast range of spatial and temporal scales. Computational methods have advanced our understanding of the behavior of molecular systems by enabling us to test assumptions and hypotheses, explore the effect of different parameters on the outcome, and eventually guide experiments. While several different mathematical and computational methods are developed to study molecular systems at different spatiotemporal scales, there is still a need for methods that bridge the gap between spatially-detailed and computationally-efficient approaches. In this review, we summarize the capabilities of agent-based modeling (ABM) as an emerging molecular systems biology technique that provides researchers with a new tool in exploring the dynamics of molecular systems/pathways in health and disease

Quality control of mRNAs at the entry of the nuclear pore: Cooperation in a complex molecular system.

Despite extensive research on how mRNAs are quality controlled prior to export into the cytoplasm, the exact underlying mechanisms are still under debate. Specifically, it is unclear how quality control proteins at the entry of the nuclear pore complex (NPC) distinguish normal and aberrant mRNAs. While some of the involved components are suggested to act as switches and recruit different factors to normal versus aberrant mRNAs, some experimental and computational evidence suggests that the combined effect of the regulated stochastic interactions between the involved components could potentially achieve an efficient quality control of mRNAs. In this review, we present a state-of-the-art portrait of the mRNA quality control research and discuss the current hypotheses proposed for dynamics of the cooperation between the involved components and how it leads to their shared goal: mRNA quality control prior to export into the cytoplasm

Higher nucleoporin-Importinβ affinity at the nuclear basket increases nucleocytoplasmic import.

Several in vitro studies have shown the presence of an affinity gradient in nuclear pore complex proteins for the import receptor Importinβ, at least partially contributing to nucleocytoplasmic transport, while others have historically argued against the presence of such a gradient. Nonetheless, the existence of an affinity gradient has remained an uncharacterized contributing factor. To shed light on the affinity gradient theory and better characterize how the existence of such an affinity gradient between the nuclear pore and the import receptor may influence the nucleocytoplasmic traffic, we have developed a general-purpose agent based modeling (ABM) framework that features a new method for relating rate constants to molecular binding and unbinding probabilities, and used our ABM approach to quantify the effects of a wide range of forward and reverse nucleoporin-Importinβ affinity gradients. Our results indicate that transport through the nuclear pore complex is maximized with an effective macroscopic affinity gradient of 2000 µM, 200 µM and 10 µM in the cytoplasmic, central channel and nuclear basket respectively. The transport rate at this gradient is approximately 10% higher than the transport rate for a comparable pore lacking any affinity gradient, which has a peak transport rate when all nucleoporins have an affinity of 200 µM for Importinβ. Furthermore, this optimal ratio of affinity gradients is representative of the ratio of affinities reported for the yeast nuclear pore complex--suggesting that the affinity gradient seen in vitro is highly optimized

Agent-Based Modeling in Molecular Systems Biology.

Molecular systems orchestrating the biology of the cell typically involve a complex web of interactions among various components and span a vast range of spatial and temporal scales. Computational methods have advanced our understanding of the behavior of molecular systems by enabling us to test assumptions and hypotheses, explore the effect of different parameters on the outcome, and eventually guide experiments. While several different mathematical and computational methods are developed to study molecular systems at different spatiotemporal scales, there is still a need for methods that bridge the gap between spatially-detailed and computationally-efficient approaches. In this review, we summarize the capabilities of agent-based modeling (ABM) as an emerging molecular systems biology technique that provides researchers with a new tool in exploring the dynamics of molecular systems/pathways in health and disease

On the nuclear pore complex and its emerging role in cellular mechanotransduction.

The nuclear pore complex (NPC) is a large protein assembly that perforates the nuclear envelope and provides a sole gateway for traffic between the cytoplasm and the nucleus. The NPC controls the nucleocytoplasmic transport by selectively allowing cargoes such as proteins and mRNA to pass through its central channel, thereby playing a vital role in protecting the nuclear component and regulating gene expression and protein synthesis. The selective transport through the NPC originates from its exquisite molecular structure featuring a large scaffold and the intrinsically disordered central channel domain, but the exact mechanism underlying the selective transport remains elusive and is the subject of various, often conflicting, hypotheses. Moreover, recent studies have suggested a new role for the NPC as a mechanosensor, where the NPC changes its channel diameter depending on the nuclear envelope tension, altering the molecular transportability through this nanopore. In this mini-review, we summarize the current understandings of the selective nature of the NPC and discuss its emerging role in cellular mechanotransduction

Localized lipid packing of transmembrane domains impedes integrin clustering.

Integrin clustering plays a pivotal role in a host of cell functions. Hetero-dimeric integrin adhesion receptors regulate cell migration, survival, and differentiation by communicating signals bidirectionally across the plasma membrane. Thus far, crystallographic structures of integrin components are solved only separately, and for some integrin types. Also, the sequence of interactions that leads to signal transduction remains ambiguous. Particularly, it remains controversial whether the homo-dimerization of integrin transmembrane domains occurs following the integrin activation (i.e. when integrin ectodomain is stretched out) or if it regulates integrin clustering. This study employs molecular dynamics modeling approaches to address these questions in molecular details and sheds light on the crucial effect of the plasma membrane. Conducting a normal mode analysis of the intact αllbβ3 integrin, it is demonstrated that the ectodomain and transmembrane-cytoplasmic domains are connected via a membrane-proximal hinge region, thus merely transmembrane-cytoplasmic domains are modeled. By measuring the free energy change and force required to form integrin homo-oligomers, this study suggests that the β-subunit homo-oligomerization potentially regulates integrin clustering, as opposed to α-subunit, which appears to be a poor regulator for the clustering process. If α-subunits are to regulate the clustering they should overcome a high-energy barrier formed by a stable lipid pack around them. Finally, an outside-in activation-clustering scenario is speculated, explaining how further loading the already-active integrin affects its homo-oligomerization so that focal adhesions grow in size

A strain-based finite element model for calcification progression in aortic valves.

Calcific aortic valve disease (CAVD) is a serious disease affecting the aging population. A complex interaction between biochemicals, cells, and mechanical cues affects CAVD initiation and progression. In this study, motivated by the progression of calcification in regions of high strain, we developed a finite element method (FEM) based spatial calcification progression model. Several cardiac cycles of transient structural FEM simulations were simulated. After each simulation cycle, calcium deposition was placed in regions of high circumferential strain. Our results show the radial expansion of calcification as spokes starting from the attachment region, agreeing very well with the reported clinical data

On the significance of microtubule flexural behavior in cytoskeletal mechanics.

Quantitative description of cell mechanics has challenged biological scientists for the past two decades. Various structural models have been attempted to analyze the structure of the cytoskeleton. One important aspect that has been largely ignored in all these modeling approaches is related to the flexural and buckling behavior of microtubular filaments. The objective of this paper is to explore the influence of this flexural and buckling behavior in cytoskeletal mechanics.In vitro the microtubules are observed to buckle in the first mode, reminiscent of a free, simply-supported beam. In vivo images of microtubules, however, indicate that the buckling mostly occurs in higher modes. This buckling mode switch takes place mostly because of the lateral support of microtubules via their connections to actin and intermediate filaments. These lateral loads are exerted throughout the microtubule length and yield a considerable bending behavior that, unless properly accounted for, would produce erroneous results in the modeling and analysis of the cytoskeletal mechanics.One of the promising attempts towards mechanical modeling of the cytoskeleton is the tensegrity model, which simplifies the complex network of cytoskeletal filaments into a combination merely of tension-bearing actin filaments and compression-bearing microtubules. Interestingly, this discrete model can qualitatively explain many experimental observations in cell mechanics. However, evidence suggests that the simplicity of this model may undermine the accuracy of its predictions, given the model's underlying assumption that "every single member bears solely either tensile or compressive behavior," i.e. neglecting the flexural behavior of the microtubule filaments. We invoke an anisotropic continuum model for microtubules and compare the bending energy stored in a single microtubule with its axial strain energy at the verge of buckling. Our results suggest that the bending energy can exceed the axial energy of microtubules by 40 folds. A modification to tensegrity model is, therefore, proved necessary in order to take into account the flexural response of microtubules. The concept of "bendo-tensegrity" is proposed as a modification to contemporary cytoskeletal tensegrity models