Modeling cellular processes in 3D

Recent advances in photonic imaging and fluorescent protein technology offer unprecedented views of molecular space-time dynamics in living cells. At the same time, advances in computing hardware and software enable modeling of ever more complex systems, from global climate to cell division. As modeling and experiment become more closely integrated we must address the issue of modeling cellular processes in 3D. Here, we highlight recent advances related to 3D modeling in cell biology. While some processes require full 3D analysis, we suggest that others are more naturally described in 2D or 1D. Keeping the dimensionality as low as possible reduces computational time and makes models more intuitively comprehensible; however, the ability to test full 3D models will build greater confidence in models generally and remains an important emerging area of cell biological modeling.     

miR-335 suppresses migration and invasion by targeting ROCK1 in osteosarcoma cells

Dynamic remodeling of the actin cytoskeleton is crucial for biological processes such as cell migration and cell spreading. S100A10 is a member of the S100 protein family and is involved in intracellular trafficking and cell migration. In this study, we examined the role of S100A10 in actin cytoskeletal organization and cell spreading. Depletion of S100A10 induced disruption of stress fiber formation and delay in cell spreading. Rac1 activation during spreading was suppressed by S100A10 knockdown, and exogenous expression of active Rac1 restored the ability of cells to spread in the absence of S100A10. Our results demonstrate the crucial role of S100A10 in actin dynamics promoting cell spreading via Rac1 activation.     

Vibrational normal modes and dynamical stability of DNA triplex poly(dA). 2poly(dT): S-type structure is more stable and in better agreement with observations in solution.

A normal-mode and statistical mechanical calculation was carried out to determine the vibrational normal modes, contribution of internal fluctuations to the free energy, and hydrogen bond disruption of DNA triplex poly(dA).2poly(dT). The calculation was performed on both the x-ray fiber diffraction model with a N-type sugar conformation, and a newly proposed model with a S-type sugar conformation. Our calculated normal modes for the S-type structure are in better agreement with observed IR spectra for samples in D2O solution. We also find that the contribution of internal fluctuations to free energy, premelting hydrogen bond disruption probability, and hydrogen bond melting temperatures for the Hoogsteen and Watson-Crick hydrogen bonds all show that the S-type structure is dynamically more stable than the N-type structure in a nominal solution environment. Therefore our calculation supports experimental findings that the triplex d(T)n.d(A)nd(T)n most likely adopts a S-type sugar conformation in solution or at high humidity. Our calculations, however, do not preclude the possibility of an N-type conformation at lower humidities.     

Structure-function studies of aromatase and its inhibitors: a progress report

The utilization of computer modeling, site-directed mutagenesis, inhibition kinetic analysis and reaction metabolite analysis allows us to better understand the structure–function relationship between aromatase and its inhibitors. Our results have helped in determining how steroidal and nonsteriodal aromatase inhibitors bind to the active site of the enzyme. This information has also aided in the understanding of the reaction mechanism of aromatase. Furthermore, our structure–function studies of aromatase have generated important information for predicting how environmental chemicals interact with the enzyme. During the last 2 years, a new aromatase computer model based on the X-ray structure of rabbit cytochrome P450 2C5 has been generated and used to evaluate the results obtained from new aromatase mutants produced in this laboratory. In addition, we have succeeded in the expression and purification of functionally active aromatase using an Escherichia coli expression method. The catalytic properties of this recombinant aromatase are similar to those properties exhibited by the human placental aromatase preparation and the mammalian cell-expressed enzyme. The E. coli expressed aromatase will be very useful for further structure–function studies of aromatase. Our laboratory has also evaluated the growth-inhibiting activity of aromatase inhibitors in estrogen receptor-positive breast cancer using three-dimensional cell cultures of aromatase-over expressing MCF-7 and T-47D cell lines (i.e. MCF-7aro and T-47Daro). Our results demonstrate that these three-dimensional cultures are valuable approaches to assess the growth-inhibiting activity of aromatase inhibitors. Finally, we have identified several phytochemicals to be potent inhibitors of aromatase. To demonstrate the impact of the phytochemicals on estrogen formation in vivo, we showed that the intake of anti-aromatase chemicals from red wine was capable of suppressing MCF-7aro-mediated tumor formation in nude mice and aromatase-induced hyperplasia in a transgenic mouse model in which aromatase is over-expressed in the mammary tissue.     

Myosin II and Arp2/3 cross-talk governs intracellular hydraulic pressure and lamellipodia formation

Human fibroblasts can switch between lamellipodia-dependent and -independent migration mechanisms on two-dimensional surfaces and in three-dimensional (3D) matrices. RhoA GTPase activity governs the switch from low-pressure lamellipodia to high-pressure lobopodia in response to the physical structure of the 3D matrix. Inhibiting actomyosin contractility in these cells reduces intracellular pressure and reverts lobopodia to lamellipodial protrusions via an unknown mechanism. To test the hypothesis that high pressure physically prevents lamellipodia formation, we manipulated pressure by activating RhoA or changing the osmolarity of the extracellular environment and imaged cell protrusions. We find RhoA activity inhibits Rac1-mediated lamellipodia formation through two distinct pathways. First, RhoA boosts intracellular pressure by increasing actomyosin contractility and water influx but acts upstream of Rac1 to inhibit lamellipodia formation. Increasing osmotic pressure revealed a second RhoA pathway, which acts through nonmuscle myosin II (NMII) to disrupt lamellipodia downstream from Rac1 and elevate pressure. Interestingly, Arp2/3 inhibition triggered a NMII-dependent increase in intracellular pressure, along with lamellipodia disruption. Together, these results suggest that actomyosin contractility and water influx are coordinated to increase intracellular pressure, and RhoA signaling can inhibit lamellipodia formation via two distinct pathways in high-pressure cells.     

A Spatial Model for Integrin Clustering as a Result of Feedback between Integrin Activation and Integrin Binding

Integrins are transmembrane adhesion receptors that bind extracellular matrix (ECM) proteins and signal bidirectionally to regulate cell adhesion and migration. In many cell types, integrins cluster at cell-ECM contacts to create the foundation for adhesion complexes that transfer force between the cell and the ECM. Even though the temporal and spatial regulation of these integrin clusters is essential for cell migration, how cells regulate their formation is currently unknown. It has been shown that integrin cluster formation is independent of actin stress fiber formation, but requires active (high-affinity) integrins, phosphoinositol-4,5-bisphosphate (PIP2), talin, and immobile ECM ligand. Based on these observations, we propose a minimal model for initial formation of integrin clusters, facilitated by localized activation and binding of integrins to ECM ligands as a result of biochemical feedback between integrin binding and integrin activation. By employing a diffusion-reaction framework for modeling these reactions, we show how spatial organization of bound integrins into clusters may be achieved by a local source of active integrins, namely protein complexes formed on the cytoplasmic tails of bound integrins. Further, we show how such a mechanism can turn small local increases in the concentration of active talin or active integrin into integrin clusters via positive feedback. Our results suggest that the formation of integrin clusters by the proposed mechanism depends on the relationships between production and diffusion of integrin-activating species, and that changes to the relative rates of these processes may affect the resulting properties of integrin clusters.     

Force Spectrum Microscopy Using Mitochondrial Fluctuations of Control and ATP-Depleted Cells

A single-cell assay of active and passive intracellular mechanical properties of mammalian cells could give significant insight into cellular processes. Force spectrum microscopy (FSM) is one such technique, which combines the spontaneous motion of probe particles and the mechanical properties of the cytoskeleton measured by active microrheology using optical tweezers to determine the force spectrum of the cytoskeleton. A simpler and noninvasive method to perform FSM would be very useful, enabling its widespread adoption. Here, we develop an alternative method of FSM using measurement of the fluctuating motion of mitochondria. Mitochondria of the C3H-10T1/2 cell line were labeled and tracked using confocal microscopy. Mitochondrial probes were selected based on morphological characteristics, and their mean-square displacement, creep compliance, and distributions of directional change were measured. We found that the creep compliance of mitochondria resembles that of particles in viscoelastic media. However, comparisons of creep compliance between controls and cells treated with pharmacological agents showed that perturbations to the actomysoin network had surprisingly small effects on mitochondrial fluctuations, whereas microtubule disruption and ATP depletion led to a significantly decreased creep compliance. We used properties of the distribution of directional change to identify a regime of thermally dominated fluctuations in ATP-depleted cells, allowing us to estimate the viscoelastic parameters for a range of timescales. We then determined the force spectrum by combining these viscoelastic properties with measurements of spontaneous fluctuations tracked in control cells. Comparisons with previous measurements made using FSM revealed an excellent match.     

Three-dimensional spherical gelatin bubble-based scaffold improves the myotube formation of H9c2 myoblasts

Microenvironmental factors including physical and chemical cues can regulate stem cells as well as terminally differentiated cells to modulate their biological function and differentiation. However, one of the physical cues, the substrate's dimensionality, has not been studied extensively. In this study, the flow-focusing method with a microfluidic device was used to generate gelatin bubbles to fabricate highly ordered three-dimensional (3D) scaffolds. Rat H9c2 myoblasts were seeded into the 3D gelatin bubble-based scaffolds and compared to those grown on 2D gelatin-coating substrates to demonstrate the influences of spatial cues on cell behaviors. Relative to cells on the 2D substrates, the H9c2 myoblasts were featured by a good survival and normal mitochondrial activity but slower cell proliferation within the 3D scaffolds. The cortical actin filaments of H9c2 cells were localized close to the cell membrane when cultured on the 2D substrates, while the F-actins distributed uniformly and occupied most of the cell cytoplasm within the 3D scaffolds. H9c2 myoblasts fused as multinuclear myotubes within the 3D scaffolds without any induction but cells cultured on the 2D substrates had a relatively lower fusion index even differentiation medium was provided. Although there was no difference in actin α 1 and myosin heavy chain 1, H9c2 cells had a higher myogenin messenger RNA level in the 3D scaffolds than those of on the 2D substrates. This study reveals that the dimensionality influences differentiation and fusion of myoblasts.     

Role for a Cindr-Arf6 axis in patterning emerging epithelia

Patterning of the Drosophila pupal eye is characterized by precise cell movements. In this paper, we demonstrate that these movements require an Arf regulatory cycle that connects surface receptors to actin-based movement. dArf6 activity-regulated by the Arf GTPase-activating proteins (ArfGAPs) dAsap and dArfGAP3 and the Arf GTP exchange factors Schizo and dPsd-promoted large cellular extensions; time-lapse microscopy indicated that these extensions presage cell rearrangements into correct epithelial niches. During this process, the Drosophila eye also requires interactions between surface Neph1/nephrin adhesion receptors Roughest and Hibris, which bind the adaptor protein Cindr (CD2AP). We provide evidence that Cindr forms a physical complex with dArfGAP3 and dAsap. Our data suggest this interaction sequesters ArfGAP function to liberate active dArf6 elsewhere in the cell. We propose that a Neph1/nephrin-Cindr/ArfGAP complex accumulates to limit local Arf6 activity and stabilize adherens junctions. Our model therefore links surface adhesion via an Arf6 regulatory cascade to dynamic modeling of the cytoskeleton, accounting for precise cell movements that organize the functional retinal field. Further, we demonstrate a similar relationship between the mammalian Cindr orthologue CD2AP and Arf6 activity in cell motility assays. We propose that this Cindr/CD2AP-mediated regulation of Arf6 is a widely used mechanism in emerging epithelia.     

D2 dopamine receptors interact directly with N-type calcium channels and regulate channel surface expression levels

N-type channels are located on dendrites and at pre-synaptic nerve terminals where they play a fundamental role in neurotransmitter release. They are potently regulated by the activation of a number of different types of pertussis toxin (PTX)-sensitive Gαi/o coupled receptors, which results in voltage-dependent inhibition of channel activity via Gβγ subunits. Using heterologous expression in HEK 293T cells, we show via whole cell patch clamp recordings that D2 receptors mediate both Gβγ (i.e. voltage-dependent) and voltage-independent inhibition of channel activity. Furthermore, using co-immunoprecipitation and pull down assays involving the intracellular regions of each protein, we show that D2 receptors and N-type channels form physical signaling complexes. Finally, we use confocal microscopy to demonstrate that D2 receptors regulate N-type channel trafficking to affect the number of calcium channels available at the plasma membrane. Taken together, these data provide evidence for multiple voltage-dependent and voltage-independent mechanisms by which D2 receptor subtypes influence N-type channel activity.     

Microrheology and ROCK signaling of human endothelial cells embedded in a 3D matrix

Cell function is profoundly affected by the geometry of the extracellular environment confining the cell. Whether and how cells plated on a two-dimensional matrix or embedded in a three-dimensional (3D) matrix mechanically sense the dimensionality of their environment is mostly unknown, partly because individual cells in an extended matrix are inaccessible to conventional cell-mechanics probes. Here we develop a functional assay based on multiple particle tracking microrheology coupled with ballistic injection of nanoparticles to measure the local intracellular micromechanical properties of individual cells embedded inside a matrix. With our novel assay, we probe the mechanical properties of the cytoplasm of individual human umbilical vein endothelial cells (HUVECs) embedded in a 3D peptide hydrogel in the presence or absence of vascular endothelial growth factor (VEGF). We found that VEGF treatment, which enhances endothelial migration, increases the compliance and reduces the elasticity of the cytoplasm of HUVECs in a matrix. This VEGF-induced softening response of the cytoplasm is abrogated by specific Rho-kinase (ROCK) inhibition. These results establish combined particle-tracking microrheology and ballistic injection as the first method able to probe the micromechanical properties and mechanical response to agonists and/or drug treatments of individual cells inside a matrix. These results suggest that ROCK plays an essential role in the regulation of the intracellular mechanical response to VEGF of endothelial cells in a 3D matrix.     

D2 dopamine receptors interact directly with N-type calcium channels and regulate channel surface expression levels

N-type channels are located on dendrites and at pre-synaptic nerve terminals where they play a fundamental role in neurotransmitter release. They are potently regulated by the activation of a number of different types of pertussis toxin (PTX)-sensitive G alpha(i/o) coupled receptors, which results in voltage-dependent inhibition of channel activity via G betagamma subunits. Using heterologous expression in HEK 293T cells, we show via whole cell patch clamp recordings that D2 receptors mediate both G betagamma (i.e., voltage-dependent) and voltage-independent inhibition of channel activity. Furthermore, using co-immunoprecipitation and pull down assays involving the intracellular regions of each protein, we show that D2 receptors and N-type channels form physical signaling complexes. Finally, we use confocal microscopy to demonstrate that D2 receptors regulate N-type channel trafficking to affect the number of calcium channels available at the plasma membrane. Taken together, these data provide evidence for multiple voltage-dependent and voltage-independent mechanisms by which D2 receptor subtypes influence N-type channel activity.     

The dynamics of camphor in the cytochrome P450 CYP101D2

The recent crystal structures of CYP101D2, a cytochrome P450 protein from the oligotrophic bacterium Novosphingobium aromaticivorans DSM12444 revealed that both the native (substrate‐free) and camphor‐soaked forms have open conformations. Furthermore, two other potential camphor‐binding sites were also identified from electron densities in the camphor‐soaked structure, one being located in the access channel and the other in a cavity on the surface near the F‐helix side of the F‐G loop termed the substrate recognition site. These latter sites may be key intermediate positions on the pathway for substrate access to or product egress from the active site. Here, we show via the use of unbiased atomistic molecular dynamics simulations that despite the open conformation of the native and camphor‐bound crystal structures, the underlying dynamics of CYP101D2 appear to be very similar to other CYP proteins. Simulations of the native structure demonstrated that the protein is capable of sampling many different conformational substates. At the same time, simulations with the camphor positioned at various locations within the access channel or recognition site show that movement towards the active site or towards bulk solvent can readily occur on a short timescale, thus confirming many previously reported in silico studies using steered molecular dynamics. The simulations also demonstrate how the fluctuations of an aromatic gate appear to control access to the active site. Finally, comparison of camphor‐bound simulations with the native simulations suggests that the fluctuations can be of similar level and thus are more representative of the conformational selection model rather than induced fit.     

Role of mitochondria and caspases in vitamin D-mediated apoptosis of MCF-7 breast cancer cells

Vitamin D(3) compounds are currently in clinical trials for human breast cancer and offer an alternative approach to anti-hormonal therapies for this disease. 1alpha,25-Dihydroxyvitamin D(3) (1alpha,25(OH)(2)D(3)), the active form of vitamin D(3), induces apoptosis in breast cancer cells and tumors, but the underlying mechanisms are poorly characterized. In these studies, we focused on the role of caspase activation and mitochondrial disruption in 1alpha,25(OH)(2)D(3)-mediated apoptosis in breast cancer cells (MCF-7) in vitro. The effect of 1alpha,25(OH)(2)D(3) on MCF-7 cells was compared with that of tumor necrosis factor alpha, which induces apoptosis via a caspase-dependent pathway. Our major findings are that 1alpha,25(OH)(2)D(3) induces apoptosis in MCF-7 cells by disruption of mitochondrial function, which is associated with Bax translocation to mitochondria, cytochrome c release, and production of reactive oxygen species. Moreover, we show that Bax translocation and mitochondrial disruption do not occur after 1alpha,25(OH)(2)D(3) treatment of a MCF-7 cell clone selected for resistance to 1alpha,25(OH)(2)D(3)-mediated apoptosis. These mitochondrial effects of 1alpha,25(OH)(2)D(3) do not require caspase activation, since they are not blocked by the cell-permeable caspase inhibitor z-Val-Ala-Asp-fluoromethylketone. Although caspase inhibition blocks 1alpha,25(OH)(2)D(3)-mediated events downstream of mitochondria such as poly(ADP-ribose) polymerase cleavage, external display of phosphatidylserine, and DNA fragmentation, MCF-7 cells still execute apoptosis in the presence of z-Val-Ala-Asp-fluoromethylketone, indicating that the commitment to 1alpha,25(OH)(2)D(3)-mediated cell death is caspase-independent.     

The mucin MUC4 and its membrane partner ErbB2 regulate biological properties of human CAPAN-2 pancreatic cancer cells via different signalling pathways

The mucin MUC4 and its membrane partner the ErbB2 oncogenic receptor are potential interacting partners in human pancreatic tumour development. However, the way they function is still largely unknown. In this work, we aimed to identify the cellular mechanisms and the intracellular signalling pathways under the control of both ErbB2 and MUC4 in a human pancreatic adenocarcinomatous cell line. Using co-immunoprecipitation and GST pull-down, we show that MUC4 and ErbB2 interact in the human pancreatic adenocarcinomatous cell line CAPAN-2 via the EGF domains of MUC4. Stable cell clones were generated in which either MUC4 or ErbB2 were knocked down (KD) by a shRNA approach. Biological properties of these cells were then studied in vitro and in vivo. Our results show that ErbB2-KD cells are more apoptotic and less proliferative (decreased cyclin D1 and increased p27kip1 expression) while migration and invasive properties were not altered. MUC4-KD clones were less proliferative with decreased cyclin D1 expression, G1 cell cycle arrest and altered ErbB2/ErbB3 expression. Their migration properties were reduced whereas invasive properties were increased. Importantly, inhibition of ErbB2 and MUC4 expression did not impair the same signalling pathways (inhibition of MUC4 expression affected the JNK pathway whereas that of ErbB2 altered the MAPK pathway). Finally, ErbB2-KD and MUC4-KD cells showed impaired tumour growth in vivo. Our results show that ErbB2 and MUC4, which interact physically, activate different intracellular signalling pathways to regulate biological properties of CAPAN-2 pancreatic cancer cells.     

Use of focused acoustics for cell disruption to provide ultra scale-down insights of microbial homogenization and its bioprocess impact--recovery of antibody fragments from rec E. coli

An ultra scale-down (USD) device that provides insight of how industrial homogenization impacts bioprocess performance is desirable in the biopharmaceutical industry, especially at the early stage of process development where only a small quantity of material is available. In this work, we assess the effectiveness of focused acoustics as the basis of an USD cell disruption method to mimic and study high-pressure, step-wise homogenization of rec Escherichia coli cells for the recovery of an intracellular protein, antibody fragment (Fab'). The release of both Fab' and of overall protein follows first-order reaction kinetics with respect to time of exposure to focused acoustics. The rate constant is directly proportional to applied electrical power input per unit volume. For nearly total protein or Fab' release (>99%), the key physical properties of the disruptate produced by focused acoustics, such as cell debris particle size distribution and apparent viscosity show good agreement with those for homogenates produced by high-pressure homogenization operated to give the same fractional release. The only key difference is observed for partial disruption of cells where focused acoustics yields a disruptate of lower viscosity than homogenization, evidently due to a greater extent of polynucleic acids degradation. Verification of this USD approach to cell disruption by high-pressure homogenization is achieved using USD centrifugation to demonstrate the same sedimentation characteristics of disruptates prepared using both the scaled-down focused acoustic and the pilot-scale homogenization methods for the same fraction of protein release.     

The distribution of vinculin to lipid rafts plays an important role in sensing stiffness of extracellular matrix

Extracellular matrix (ECM) stiffness regulates cell differentiation, survival, and migration. Our previous study has shown that the interaction of the focal adhesion protein vinculin with vinexin α plays a critical role in sensing ECM stiffness and regulating stiffness-dependent cell migration. However, the mechanism how vinculin–vinexin α interaction affects stiffness-dependent cell migration is unclear. Lipid rafts are membrane microdomains that are known to affect ECM-induced signals and cell behaviors. Here, we show that vinculin and vinexin α can localize to lipid rafts. Cell-ECM adhesion, intracellular tension, and a rigid ECM promote vinculin distribution to lipid rafts. The disruption of lipid rafts with Methyl-β-cyclodextrin impaired the ECM stiffness-mediated regulation of vinculin behavior and rapid cell migration on rigid ECM. These results indicate that lipid rafts play an important role in ECM-stiffness regulation of cell migration via vinculin.     

A quantitative analysis of the interplay of environment,neighborhood, and cell state in 3D spheroids

Cells react to their microenvironment by integrating external stimuli into phenotypic decisions via an intracellular signaling network. To analyze the interplay of environment, local neighborhood, and internal cell state effects on phenotypic variability, we developed an experimental approach that enables multiplexed mass cytometric imaging analysis of up to 240 pooled spheroid microtissues. We quantified the contributions of environment, neighborhood, and intracellular state to marker variability in single cells of the spheroids. A linear model explained on average more than half of the variability of 34 markers across four cell lines and six growth conditions. The contributions of cell‐intrinsic and environmental factors to marker variability are hierarchically interdependent, a finding that we propose has general implications for systems‐level studies of single‐cell phenotypic variability. By the overexpression of 51 signaling protein constructs in subsets of cells, we also identified proteins that have cell‐intrinsic and cell‐extrinsic effects. Our study deconvolves factors influencing cellular phenotype in a 3D tissue and provides a scalable experimental system, analytical principles, and rich multiplexed imaging datasets for future studies.     

Rab8 Binding to Immune Cell-Specific Adaptor LAX Facilitates Formation of trans-Golgi Network-Proximal CTLA-4 Vesicles for Surface Expression

Despite playing a central role in tolerance, little is known regarding the mechanism by which intracellular CTLA-4 is shuttled from the trans-Golgi network to the surfaces of T cells. In this context, Ras-related GTPase Rab8 plays an important role in the intracellular transport, while we have previously shown that CTLA-4 binds to the immune cell adaptor TRIM in T cells. In this study, we demonstrate that CTLA-4 forms a multimeric complex comprised of TRIM and related LAX that in turn binds to GTP bound Rab8 for post-Golgi transport to the cell surface. LAX bound via its N terminus to active GTP-Rab8, as well as the cytoplasmic tail of CTLA-4. TRIM required LAX for binding to Rab8 in a complex. Wild-type LAX or its N terminus (residues 1 to 77) increased CTLA-4 surface expression, whereas small interfering RNAs of Rab8 or LAX or disruption of LAX/Rab8 binding reduced numbers of CTLA-4-containing vesicles and its coreceptor surface expression. LAX also promoted the polarization of CTLA-4 and the reorientation of the microtubule-organizing center to the site of T-cell receptor engagement. Our results identify a novel CTLA-4/TRIM/LAX/Rab8 effector complex in the transport of CTLA-4 to the surfaces of T cells.