Regulation of Neurite Growth by Inorganic Pyrophosphatase 1 via JNK Dephosphorylation

Neural cell differentiation during development is controlled by multiple signaling pathways, in which protein phosphorylation and dephosphorylation play an important role. In this study, we examined the role of pyrophosphatase1 (PPA1) in neuronal differentiation using the loss and gain of function analysis. Neuronal differentiation induced by external factors was studied using a mouse neuroblastoma cell line (N1E115). The neuronal like differentiation in N1E115 cells was determined by morphological analysis based on neurite growth length. In order to analyze the loss of the PPA1 function in N1E115, si-RNA specifically targeting PPA1 was generated. To study the effect of PPA1 overexpression, an adenoviral gene vector containing the PPA1 gene was utilized to infect N1E115 cells. To address the need for pyrophosphatase activity in PPA1, D117A PPA1, which has inactive pyrophosphatase, was overexpressed in N1E115 cells. We used valproic acid (VPA) as a neuronal differentiator to examine the effect of PPA1 in actively differentiated N1E115 cells. Si-PPA1 treatment reduced the PPA1 protein level and led to enhanced neurite growth in N1E115 cells. In contrast, PPA1 overexpression suppressed neurite growth in N1E115 cells treated with VPA, whereas this effect was abolished in D117A PPA1. PPA1 knockdown enhanced the JNK phosphorylation level, and PPA1 overexpression suppressed it in N1E115 cells. It seems that recombinant PPA1 can dephosphorylate JNK while no alteration of JNK phosphorylation level was seen after treatment with recombinant PPA1 D117A. Enhanced neurite growth by PPA1 knockdown was also observed in rat cortical neurons. Thus, PPA1 may play a role in neuronal differentiation via JNK dephosphorylation.     

Model of Growth Cone Membrane Polarization via Microtubule Length Regulation

We present a mathematical model of membrane polarization in growth cones. We proceed by coupling an active transport model of cytosolic proteins along a two-dimensional microtubule (MT) network with a modified Dogterom-Leibler model of MT growth. In particular, we consider a Rac1-stathmin-MT pathway in which the growth and catastrophe rates of MTs are regulated by cytosolic stathmin, while the stathmin is regulated by Rac1 at the membrane. We use regular perturbation theory and numerical simulations to determine the steady-state stathmin concentration, the mean MT length distribution, and the resulting distribution of membrane-bound proteins. We thus show how a nonuniform Rac1 distribution on the membrane generates a polarized distribution of membrane proteins. The mean MT length distribution and hence the degree of membrane polarization are sensitive to the precise form of the Rac1 distribution and parameters such as the catastrophe-promoting constant and tubulin association rate. This is a consequence of the fact that the lateral diffusion of stathmin tends to weaken the effects of Rac1 on the distribution of mean MT lengths.     

Computational modeling of biomechanics and biorheology of heated red blood cells

Because of their compromised deformability, heat denatured erythrocytes have been used as labeled probes to visualize spleen tissue or to assess the ability of the spleen to retain stiff red blood cells (RBCs) for over three decades, e.g., see Looareesuwan et al. N. Engl. J. Med. (1987). Despite their good accessibility, it is still an open question how heated RBCs compare to certain diseased RBCs in terms of their biomechanical and biorheological responses, which may undermine their effective usage and even lead to misleading experimental observations. To help answering this question, we perform a systematic computational study of the hemorheological properties of heated RBCs with several physiologically relevant static and hemodynamic settings, including optical-tweezers test, relaxation of prestretched RBCs, RBC traversal through a capillary-like channel and a spleen-like slit, and a viscometric rheology test. We show that our in silico RBC models agree well with existing experiments. Moreover, under static tests, heated RBCs exhibit deformability deterioration comparable to certain disease-impaired RBCs such as those in malaria. For RBC traversal under confinement (through microchannel or slit), heated RBCs show prolonged transit time or retention depending on the level of confinement and heating procedure, suggesting that carefully heat-treated RBCs may be useful for studying splenic- or vaso-occlusion in vascular pathologies. For the rheology test, we expand the existing bulk viscosity data of heated RBCs to a wider range of shear rates (1–1000 s⁻¹) to represent most pathophysiological conditions in macro- or microcirculation. Although heated RBC suspension shows elevated viscosity comparable to certain diseased RBC suspensions under relatively high shear rates (100–1000 s⁻¹), they underestimate the elevated viscosity (e.g., in sickle cell anemia) at low shear rates (<10 s⁻¹). Our work provides mechanistic rationale for selective usage of heated RBC as a potentially useful model for studying the abnormal traversal dynamics and hemorheology in certain blood disorders.     

Self-Reproduction of Fatty Acid Vesicles: A Combined Experimental and Simulation Study

Dilution of a fatty acid micellar solution at basic pH toward neutrality results in spontaneous formation of vesicles with a broad size distribution. However, when vesicles of a defined size are present before dilution, the size distribution of the newly formed vesicles is strongly biased toward that of the seed vesicles. This so-called matrix effect is believed to be a key feature of early life. Here we reproduced this effect for oleate micelles and seed vesicles of either oleate or dioleoylphosphatidylcholine. Fluorescence measurements showed that the vesicle contents do not leak out during the replication process. We hypothesized that the matrix effect results from vesicle fission induced by an imbalance of material across both leaflets of the vesicle upon initial insertion of fatty acids into the outer leaflet of the seed vesicle. This was supported by experiments that showed a significant increase in vesicle size when the equilibration of oleate over both leaflets was enhanced by either slowing down the rate of fatty acid addition or increasing the rate of fatty acid transbilayer movement. Coarse-grained molecular-dynamics simulations showed excellent agreement with the experimental results and provided further mechanistic details of the replication process.     

Computational Tool to Study Perturbations in Muscle Regulation and Its Application to Heart Disease

Striated muscle contraction occurs when myosin thick filaments bind to thin filaments in the sarcomere and generate pulling forces. This process is regulated by calcium, and it can be perturbed by pathological conditions (e.g., myopathies), physiological adaptations (e.g., β-adrenergic stimulation), and pharmacological interventions. Therefore, it is important to have a methodology to robustly determine the impact of these perturbations and statistically evaluate their effects. Here, we present an approach to measure the equilibrium constants that govern muscle activation, estimate uncertainty in these parameters, and statistically test the effects of perturbations. We provide a MATLAB-based computational tool for these analyses, along with easy-to-follow tutorials that make this approach accessible. The hypothesis testing and error estimation approaches described here are broadly applicable, and the provided tools work with other types of data, including cellular measurements. To demonstrate the utility of the approach, we apply it to elucidate the biophysical mechanism of a mutation that causes familial hypertrophic cardiomyopathy. This approach is generally useful for studying muscle diseases and therapeutic interventions that target muscle contraction.     

The Regulative Role of Neurite Mechanical Tension in Network Development

A bewildering series of dynamical processes take part in the development of the nervous system. Neuron branching dynamics, the continuous formation and elimination of neural interconnections, are instrumental in constructing distinct neuronal networks, which are the functional building blocks of the nervous system. In this study, we investigate and validate the important regulative role of mechanical tension in determining the final morphology of neuronal networks. To single out the mechanical effect, we cultured relatively large invertebrate neurons on clean quartz surfaces. Applied to these surfaces were isolated anchoring sites consisting of carbon nanotube islands to which the cells and the neurites could mechanically attach. Inspection of branching dynamics and network wiring upon development revealed an innate selection mechanism in which one axon branch wins over another. The apparent mechanism entails the build-up of mechanical tension in developing axons. The tension is maintained by the attachment of the growth cone to the substrate or, alternatively, to the neurites of a target neuron. The induced tension promotes the stabilization of one set of axon branches while causing retraction or elimination of axon collaterals. We suggest that these findings represent a crucial, early step that precedes the formation of synapses and regulates neuronal interconnections. Mechanical tension serves as a signal for survival of the axonal branch and perhaps for the subsequent formation of synapses.     

A New Computational Method for Membrane Compressibility: Bilayer Mechanical Thickness Revisited

Because lipid bilayers can bend and stretch in ways similar to thin elastic sheets, physical models of bilayer deformation have utilized mechanical constants such as the moduli for bending rigidity $\left({\kappa }_C\right)$ and area compressibility (K_A). However, the use of these models to quantify the energetics of membrane deformation associated with protein-membrane interactions, and the membrane response to stress is often hampered by the shortage of experimental data suitable for the estimation of the mechanical constants of various lipid mixtures. Although computational tools such as molecular dynamics simulations can provide alternative means to estimate K_A values, current approaches suffer significant technical limitations. Here, we present a novel, to our knowledge, computational framework that allows for a direct estimation of K_A values for individual bilayer leaflets. The theory is based on the concept of elasticity and derives K_A from real-space analysis of local thickness fluctuations sampled in molecular dynamics simulations. We explore and validate the model on a large set of single and multicomponent bilayers of different lipid compositions and sizes, simulated at different temperatures. The calculated bilayer compressibility moduli agree with values estimated previously from experiments and those obtained from a standard computational method based on a series of constrained tension simulations. We further validate our framework in a comparison with an existing polymer brush model and confirm the polymer brush model’s predicted linear relationship with proportionality coefficient of 24, using elastic parameters calculated from the simulation trajectories. The robustness of the results that emerge from the method allows us to revisit the origins of the bilayer mechanical (compressible) thickness and in particular its dependence on acyl-chain unsaturation and the presence of cholesterol.     

Ensembles of Breathing Nucleosomes: A Computational Study

About three-fourths of the human DNA molecules are wrapped into nucleosomes, protein spools with DNA. Nucleosomes are highly dynamic, transiently exposing their DNA through spontaneous unspooling. Recent experiments allowed to observe the DNA of an ensemble of such breathing nucleosomes through x-ray diffraction with contrast matching between the solvent and the protein core. In this study, we calculate such an ensemble through a Monte Carlo simulation of a coarse-grained nucleosome model with sequence-dependent DNA mechanics. Our analysis gives detailed insights into the sequence dependence of nucleosome breathing observed in the experiment and allows us to determine the adsorption energy of the DNA bound to the protein core as a function of the ionic strength. Moreover, we predict the breathing behavior of other potentially interesting sequences and compare the findings to earlier related experiments.     

Lipid Raft-Mediated Regulation of G-Protein Coupled Receptor Signaling by Ligands which Influence Receptor Dimerization: A Computational Study

G-protein coupled receptors (GPCRs) are the largest family of cell surface receptors; they activate heterotrimeric G-proteins in response to ligand stimulation. Although many GPCRs have been shown to form homo- and/or heterodimers on the cell membrane, the purpose of this dimerization is not known. Recent research has shown that receptor dimerization may have a role in organization of receptors on the cell surface. In addition, microdomains on the cell membrane termed lipid rafts have been shown to play a role in GPCR localization. Using a combination of stochastic (Monte Carlo) and deterministic modeling, we propose a novel mechanism for lipid raft partitioning of GPCRs based on reversible dimerization of receptors and then demonstrate that such localization can affect GPCR signaling. Modeling results are consistent with a variety of experimental data indicating that lipid rafts have a role in amplification or attenuation of G-protein signaling. Thus our work suggests a new mechanism by which dimerization-inducing or inhibiting characteristics of ligands can influence GPCR signaling by controlling receptor organization on the cell membrane.     

The Effect of Perfusion on Soft Tissue Mechanical Properties: A Computational Model

Some simple finite element models were constructed to investigate the magnitude and character of changes in mechanical properties of very soft tissues due to the loss of perfusion. Changes in the apparent stress-strain curve were used to characterise the effect of simulated blood perfusion pressure on the engineering stress-strain curve. The results indicated that the blood to tissue volume ratio and the perfusion pressure have the strongest effect on the effective stress-strain response of a representative tissue cell. Tissue viscoelasticity increased the sensitivity of the system to perfusion pressure changes.     

Varicones and Growth Cones: Two Neurite Terminals in PC12 Cells

The rat adrenal pheochromocytoma PC12 cell line is one of the traditional models for the study of neurite outgrowth and growth cone behavior. To clarify to what extent PC12 neurite terminals can be compared to neuronal growth cones, we have analyzed their morphology and protein distribution in fixed PC12 cells by immunocytochemistry. Our results show that that PC12 cells display a special kind of neurite terminal that includes a varicosity in close association with a growth cone. This hybrid terminal, or “varicone”, is characterized by the expression of specific markers not typically present in neuronal growth cones. For example, we show that calpain-2 is a specific marker of varicones and can be detected even before the neurite develops. Our data also shows that a fraction of PC12 neurites end in regular growth cones, which we have compared to hippocampal neurites as a control. We also report the extraordinary incidence of varicones in the literature referred to as “growth cones”. In summary, we provide evidence of two different kinds of neurite terminals in PC12 cells, including a PC12-specific terminal, which implies that care must be taken when using them as a model for neuronal growth cones or neurite outgrowth.     

Functional Analysis of Metabolic Channeling and Regulation in Lignin Biosynthesis: A Computational Approach

Lignin is a polymer in secondary cell walls of plants that is known to have negative impacts on forage digestibility, pulping efficiency, and sugar release from cellulosic biomass. While targeted modifications of different lignin biosynthetic enzymes have permitted the generation of transgenic plants with desirable traits, such as improved digestibility or reduced recalcitrance to saccharification, some of the engineered plants exhibit monomer compositions that are clearly at odds with the expected outcomes when the biosynthetic pathway is perturbed. In Medicago, such discrepancies were partly reconciled by the recent finding that certain biosynthetic enzymes may be spatially organized into two independent channels for the synthesis of guaiacyl (G) and syringyl (S) lignin monomers. Nevertheless, the mechanistic details, as well as the biological function of these interactions, remain unclear. To decipher the working principles of this and similar control mechanisms, we propose and employ here a novel computational approach that permits an expedient and exhaustive assessment of hundreds of minimal designs that could arise in vivo. Interestingly, this comparative analysis not only helps distinguish two most parsimonious mechanisms of crosstalk between the two channels by formulating a targeted and readily testable hypothesis, but also suggests that the G lignin-specific channel is more important for proper functioning than the S lignin-specific channel. While the proposed strategy of analysis in this article is tightly focused on lignin synthesis, it is likely to be of similar utility in extracting unbiased information in a variety of situations, where the spatial organization of molecular components is critical for coordinating the flow of cellular information, and where initially various control designs seem equally valid.     

微RNA介导巨噬细胞极化的转录后调控

巨噬细胞极化是根据周围刺激环境做出表型调节的一个过程.一般极化为2个表型,分别为经典激活的M1巨噬细胞和替代激活的M2巨噬细胞.简而言之,M1巨噬细胞的特征是促炎和抗肿瘤;M2巨噬细胞是抗炎和促肿瘤.巨噬细胞极化被认为是人体生理和病理的关键调节器,其发挥作用的有效性依赖于关键因子的协调表达,而这些关键因子的表达在转录后...