Search

We examined keratin aggregate formation and the possible mechanisms involved. With this aim, we observed the effect that different ratios between mutant and wild-type keratins expressed in cultured keratinocytes may have on aggregate formation in vitro, as well as how keratin aggregate formation affects the mechanical properties of cells at the cell cortex. To this end we prepared clones with expression rates as close as possible to 25%, 50% and 100% of the EGFP-K14 proteins (either WT or R125P and V270M mutants). Our results showed that only in the case of the 25% EGFP-K14 R125P mutant significant differences could be seen. Namely, we observed in this case the largest accumulation of keratin aggregates and a significant reduction in cell stiffness. To gain insight into the possible mechanisms behind this observation, we extended our previous mathematical model of keratin dynamics by implementing a more complex reaction network that considers the coexistence of wild-type and mutant keratins in the cell. The new model, consisting of a set of coupled, non-linear, ordinary differential equations, allowed us to draw conclusions regarding the relative amounts of intermediate filaments and aggregates in cells, and suggested that aggregate formation by asymmetric binding between wild-type and mutant keratins could explain the data obtained on cells grown in culture.

Keratin intermediate filaments are the principal structural element of epithelial cells. Their importance in providing bulk cellular stiffness is well recognized, but their role in the mechanics of cell cortex is less understood. In this study, we therefore compared the cortical stiffness of three keratinocyte lines: primary wild type cells (NHEK2), immortalized wild type cells (NEB1) and immortalized mutant cells (KEB7). The cortical stiffness was measured by lateral indentation of cells with AOD-steered optical tweezers without employing any moving mechanical elements. The method was validated on fixed cells and Cytochalasin-D treated cells to ensure that the observed variations in stiffness within a single cell line were not a consequence of low measurement precision. The measurements of the cortical stiffness showed that primary wild type cells were significantly stiffer than immortalized wild type cells, which was also detected in previous studies of bulk elasticity. In addition, a small difference between the mutant and the wild type cells was detected, showing that mutation of keratin impacts also the cell cortex. Thus, our results indicate that the role of keratins in cortical stiffness is not negligible and call for further investigation of the mechanical interactions between keratins and elements of the cell cortex.

We analyze the effect that the geometrical place of anastomosis in the circulatory tree has on blood flow. We introduce an idealized model that consists of a symmetric network for the arterial and venous vascular trees. We consider that the network contains a viscoelastic fluid with the rheological characteristics of blood, and analyze the network hydrodynamic response to a time-dependent periodic pressure gradient. This response is a measurement of the resistance to flow: the larger the response, the smaller the resistance to flow. We find that for networks whose vessels have the same radius and length, the outer the level of the branching tree in which anastomosis occurs, the larger the network response. Moreover, when anastomosis is incorporated in the form of bypasses that bridge vessels at different bifurcation levels, the further apart are the levels bridged by the bypass, the larger the response is. Furthermore, we apply the model to the available information for the dog circulatory system and find that the effect that anastomosis causes at different bifurcation levels is strongly determined by the structure of the underlying network without anastomosis. We rationalize our results by introducing two idealized models and approximated analytical expressions that allow us to argue that, to a large extent, the response of the network with anastomosis is determined locally. We have also considered the influence of the myogenic effect. This one has a large quantitative impact on the network response. However, the qualitative behavior of the network response with anastomosis is the same with or without consideration of the myogenic effect. That is, it depends on the structure that the underlying vessel network has in a small neighborhood around the place where anastomosis occurs. This implies that whenever there is an underlying tree-like network in an in vivo vasculature, our model is able to interpret the anastomotic effect.

Tumor growth is becoming a central problem in biophysics both from its social and medical interest and, more fundamentally, because it is a remarkable example of an emergent complex system. Focusing on the description of the spatial and dynamical features of tumor growth, in this paper we review recent tumor modeling approaches using a technique borrowed from materials science: the phase-field models. These models allow us, with a large degree of generality, to identify the paramount mechanisms causing the uncontrolled growth of tumor cells as well as to propose new guidelines for experimentation both in simulation and in the laboratory. We finish by discussing open directions of research in phase-field modeling of tumor growth to catalyze the interest of physicists and mathematicians in this emergent field.

Systematic Monte Carlo simulations of simple lattice models show that the final stage of protein folding is an ordered process where native contacts get locked (i.e., the residues come into contact and remain in contact for the duration of the folding process) in a well-defined order. The detailed study of the folding dynamics of protein-like sequences designed as to exhibit different contact energy distributions, as well as different degrees of sequence optimization (i.e., participation of non-native interactions in the folding process), reveals significant differences in the corresponding locking scenarios—the collection of native contacts and their average locking times, which are largely ascribable to the dynamics of non-native contacts. Furthermore, strong evidence for a positive role played by non-native contacts at an early folding stage was also found. Interestingly, for topologically simple target structures, a positive interplay between native and non-native contacts is observed also toward the end of the folding process, suggesting that non-native contacts may indeed affect the overall folding process. For target models exhibiting clear two-state kinetics, the relation between the nucleation mechanism of folding and the locking scenario is investigated. Our results suggest that the stabilization of the folding transition state can be achieved through the establishment of a very small network of native contacts that are the first to lock during the folding process.

In the phase separation occurring at the miscibility gap (at the spinodal region) of an alloy a discrete symmetry is spontaneously broken and a domain wall network is formed. Field theory simulations are often used to study the dynamics of topological defects networks appearing in different physical contexts. In this work, we focus on the dynamics of the two immiscible liquids appearing on the phase diagram of the Bi–Zn system, one of the basic systems of lead free solders. We use phase field simulations to quantitatively simulate the dynamics of the two liquids separation in the Bi–Zn system, at different temperatures and for different concentrations. We obtain the miscibility gap curve and the domain morphologies of the system as a function of time, temperature and component concentrations using simulations.

Systematic Monte Carlo simulations of simple lattice models show that the final stage of protein folding is an ordered process where native contacts get locked (i.e., the residues come into contact and remain in contact for the duration of the folding process) in a well-defined order. The detailed study of the folding dynamics of protein-like sequences designed as to exhibit different contact energy distributions, as well as different degrees of sequence optimization (i.e., participation of non-native interactions in the folding process), reveals significant differences in the corresponding locking scenarios--the collection of native contacts and their average locking times, which are largely ascribable to the dynamics of non-native contacts. Furthermore, strong evidence for a positive role played by non-native contacts at an early folding stage was also found. Interestingly, for topologically simple target structures, a positive interplay between native and non-native contacts is observed also toward the end of the folding process, suggesting that non-native contacts may indeed affect the overall folding process. For target models exhibiting clear two-state kinetics, the relation between the nucleation mechanism of folding and the locking scenario is investigated. Our results suggest that the stabilization of the folding transition state can be achieved through the establishment of a very small network of native contacts that are the first to lock during the folding process.

In this work we model the extent of hypoxia in the diabetic retina as a function of the area affected by vessel disruption. We find two regimes that differ on the ratio between the area of disrupted vasculature and the area of tissue in hypoxia. In the first regime the hypoxia is localized in the vicinity of the vascular disruption, while in the second regime there is a generalized hypoxia in the affected tissue. The transition between these two regimes occurs when the tissue area affected by individual sites of vessel damage is on the order of the square of the characteristic irrigation length in the tissue (the maximum distance that an irrigated point in the tissue is from an existing vessel). We observe that very high levels of hypoxia are correlated with the rupture of larger vessels in the retina, and with smaller radii of individual sites of vesseldamage. Based on this property of vascular networks, we propose a novel mechanism for the transition between the nonproliferative and the roliferative stages in diabetic retinopathy.

We carry out systematic Monte Carlo simulations of Gō lattice proteins to investigate and compare the folding processes of two model proteins whose native structures differ from each other due to the presence of a trefoil knot located near the terminus of one of the protein chains. We show that the folding time of the knotted fold is larger than that of the unknotted protein and that this difference in folding time is particularly striking in the temperature region below the optimal folding temperature. Both proteins display similar folding transition temperatures, which is indicative of similar thermal stabilities. By using the folding probability reaction coordinate as an estimator of folding progression we have found out that the formation of the knot is mainly a late folding event in our shallow knot system.

Glioblastoma multiforme (GBM) is the most aggressive and invasive malignant brain cancer. GBM is characterized by a dramatic metabolic imbalance leading to increased secretion of the pro-angiogenic factor VEGF and subsequent abnormal tumor vascularization. In 2009, FDA approved the intravenous administration of bevacizumab, an anti-VEGF monoclonal antibody, as a therapeutic agent for patients with GBM. However, the number of systemic side effects and reduced accessibility of bevacizumab to the central nervous system and consequently to the GBM tumor mass limited its effectiveness in improving patient survival. In this study, we combined experimental and computational modelling to quantitatively characterize the dynamics of VEGF secretion and turnover in GBM and in normal brain cells and simultaneous monitoring of vessel growth. We showed that sequestration of VEGF inside GBM cells, can be used as a novel target for improved bevacizumab-based therapy. We have engineered the VEGF nanotrapper, a cargo system that allows cellular uptake of bevacizumab and inhibits VEGF secretion required for angiogenesis activation and development. Here, we show the therapeutic efficacy of this nanocargo in reducing vascularization and tumor cell mass of GBM in vitro and in vivo cancer models.