Search

In our cells, DNA is folded and packed with the help of many proteins into chromatin whose basic unit is a nucleosome—DNA wrapped around octamer of histone proteins. The chain of nucleosomes is further folded and arranged into many layers and has a dynamic organization. How does the complex chromatin organization emerge from interactions among DNA, histones, and non-histone proteins have been a question of great interest. Here we review recent literature that investigated how nucleosome positioning and nucleosome-mediated interactions drive chromatin organization. Unlike our earlier understanding, chromatin is organized into 3D domains of various sizes having irregularly organized nucleosomes. These domains emerge due to heterogeneous nucleosome positioning and diverse inter-nucleosome interactions that vary in space and time. [ABSTRACT FROM AUTHOR]

Accounting for enzyme-mediated active sliding, disassembly, and sequence-dependent positioning of nucleosomes, we simulate nucleosome occupancy over cell-cycle-scale times using a stochastic kinetic model. We show that ATP-dependent active nucleosome sliding and nucleosome removal processes are essential to obtain in vivo-like nucleosome positioning. While active sliding leads to dense nucleosome filling, sliding events alone cannot ensure sequence-dependent nucleosome positioning: Active nucleosome removal is the crucial remodeling event that drives positioning. We also show that remodeling activity changes nucleosome dynamics from glassy to liquid-like, and that remodeling dramatically influences exposure dynamics of promoter regions. [ABSTRACT FROM AUTHOR]

During cell devision, maintaining the epigenetic information encoded in histone modification patterns is crucial for survival and identity of cells. The faithful inheritance of the histone marks from the parental to the daughter strands is a puzzle, given that each strand gets only half of the parental nucleosomes. Mapping DNA replication and reconstruction of modifications to equivalent problems in communication of information, we ask how well enzymes can recover the parental modifications, if they were ideal computing machines. Studying a parameter regime where realistic enzymes can function, our analysis predicts that enzymes may implement a critical threshold filling algorithm which fills unmodified regions of length at most k. This algorithm, motivated from communication theory, is derived from the maximum à posteriori probability (MAP) decoding which identifies the most probable modification sequence based on available observations. Simulations using our method produce modification patterns similar to what has been observed in recent experiments. We also show that our results can be naturally extended to explain inheritance of spatially distinct antagonistic modifications. Author summary: Chromatin is essentially the DNA that is folded and packaged with the help of proteins. While the nucleotide sequence in the DNA codes genetic information, the packaging of the DNA into chromatin encodes extra layer of information. This epigenetic code regulates reading of the genetic code and provides identity to cells—whether the cell is a skin cell or a brain cell, for example. In this work we examine a long standing puzzle, that is, how the epigenetic code in the form of histone modification patterns may get inherited, when a cell divides. Using theoretical arguments, we present an algorithm that enzymes could be executing so that the epigenetic code can be inherited with minimal error, soon after DNA replication. [ABSTRACT FROM AUTHOR]

Understanding kinetic control of biological processes is as important as identifying components that constitute pathways. Insulin signaling is central for almost all metazoans, and its perturbations are associated with various developmental disorders, metabolic diseases, and aging. While temporal phosphorylation changes and kinetic constants have provided some insights, constant or variable parameters that establish and maintain signal topology are poorly understood. Here, we report kinetic parameters that encode insulin concentration and nutrient-dependent flow of information using iterative experimental and mathematical simulation-based approaches. Our results illustrate how dynamics of distinct phosphorylation events collectively contribute to selective kinetic gating of signals and maximum connectivity of the signaling cascade under normo-insulinemic but not hyper-insulinemic states. In addition to identifying parameters that provide predictive value for maintaining the balance between metabolic and growth-factor arms, we posit a kinetic basis for the emergence of insulin resistance. Given that pulsatile insulin secretion during a fasted state precedes a fed response, our findings reveal rewiring of insulin signaling akin to memory and anticipation, which was hitherto unknown. Striking disparate temporal behavior of key phosphorylation events that destroy the topology under hyper-insulinemic states underscores the importance of unraveling regulatory components that act as bandwidth filters. In conclusion, besides providing fundamental insights, our study will help in identifying therapeutic strategies that conserve coupling between metabolic and growth-factor arms, which is lost in diseases and conditions of hyper-insulinemia. [ABSTRACT FROM AUTHOR]

We present a physics-based polymer model that can investigate 3D organization of chromatin accounting for DNA elasticity, DNA-bending due to nucleosomes, and 1D organization of nucleosomes along DNA. We find that the packing density of chromatin oscillates between densities corresponding to highly folded and extended configurations as we change the nucleosome organization (length of linker DNA). We compute the looping probability of chromatin and show that the presence of nucleosomes increases the looping probability of the chain compared to that of a bare DNA. We also show that looping probability has a large variability depending on the nature of nucleosome organization and density of linker histones. [ABSTRACT FROM AUTHOR]

Characterising complex kinetics of non-equilibrium self-assembly of bio-filaments is of general interest. Dynamic instability in microtubules, consisting of successive catastrophes and rescues, is observed to occur as a result of the non-equilibrium conversion of GTP-tubulin to GDP-tubulin. We study this phenomenon using a model for microtubule kinetics with GTP/GDP state-dependent polymerisation, depolymerisation and hydrolysis of subunits. Our results reveal a sharp switch-like transition in the mean velocity of the filaments, from a growth phase to a shrinkage phase, with an associated co-existence of the two phases. This transition is reminiscent of the discontinuous phase transition across the liquid-gas boundary. We probe the extent of discontinuity in the transition quantitatively using characteristic signatures such as bimodality in velocity distribution, variance and Binder cumulant, and also hysteresis behaviour of the system. We further investigate ageing behaviour in catastrophes of the filament, and find that the multi-step nature of catastrophes is intensified in the vicinity of the switching transition. This assumes importance in the context of Microtubule Associated Proteins which have the potential of altering kinetic parameter values. [ABSTRACT FROM AUTHOR]

Gal3p is an allosteric monomeric protein that activates the GAL genetic switch of Saccharomyces cerevisiae in response to galactose. Expression of constitutive mutant of Gal3p or overexpression of wild-type Gal3p activates the GAL switch in the absence of galactose. These data suggest that Gal3p exists as an ensemble of active and inactive conformations. Structural data have indicated that Gal3p exists in open (inactive) and closed (active) conformations. However, a mutant of Gal3p that predominantly exists in inactive conformation and is yet capable of responding to galactose has not been isolated. To understand the mechanism of allosteric transition, we have isolated a triple mutant of Gal3p with V273I, T404A, and N450D substitutions, which, upon overexpression, fails to activate the GAL switch on its own but activates the switch in response to galactose. Overexpression of Gal3p mutants with single or double mutations in any of the three combinations failed to exhibit the behavior of the triple mutant. Molecular dynamics analysis of the wild-type and the triple mutant along with two previously reported constitutive mutants suggests that the wild-type Gal3p may also exist in super-open conformation. Furthermore, our results suggest that the dynamics of residue F237 situated in the hydrophobic pocket located in the hinge region drives the transition between different conformations. Based on this study, we suggest that conformational selection mechanism is the driving force in the allosteric transition of Gal3p, which may have implications in other signaling pathways involving monomeric proteins. [ABSTRACT FROM AUTHOR]

Self-assembly of proteins into ordered, fibrillar structures is a commonly observed theme in biology. It has been observed that diverse set of proteins (e.g., alpha-synuclein, insulin, TATA-box binding protein, Sup35, p53), independent of their sequence, native structure, or function could self-assemble into highly ordered structures known as amyloids. What are the crucial features underlying amyloidogenesis that make it so generic? Using coarse-grained simulations of peptide self-assembly, we argue that variation in two physical parameters—bending stiffness of the polypeptide and strength of intermolecular interactions—can give rise to many of the structural features typically associated with amyloid self-assembly. We show that the interplay between these two factors gives rise to a rich phase diagram displaying high diversity in aggregated states. For certain parameters, we find a bimodal distribution for the order parameter implying the coexistence of ordered and disordered aggregates. Our findings may explain the experimentally observed variability including the “off-pathway” aggregated structures. Further, we demonstrate that sequence-dependence and protein-specific signatures could be mapped to our coarse-grained framework to study self-assembly behavior of realistic systems such as the STVIIE peptide and Aβ42. The work also provides certain guiding principles that could be used to design novel peptides with desired self-assembly properties, by tuning a few physical parameters. [ABSTRACT FROM AUTHOR]

Phase separation of biomolecules could be mediated by both specific and non-specific interactions. How the interplay between non-specific and specific interactions along with polymer entropy influences phase separation is an open question. We address this question by simulating self-associating molecules as polymer chains with a short core stretch that forms the specifically interacting functional interface and longer non-core regions that participate in non-specific/promiscuous interactions. Our results show that the interplay of specific (strength, ϵsp) and non-specific interactions (strength, ϵns) could result in phase separation of polymers and its transition to solid-like aggregates (mature state). In the absence of ϵns, the polymer chains do not dwell long enough in the vicinity of each other to undergo phase separation and transition into a mature state. On the other hand, in the limit of strong ϵns, the assemblies cannot transition into the mature state and form a non-specific assembly, suggesting an optimal range of interactions favoring mature multimers. In the scenario where only a fraction (Nfrac) of the non-core regions participate in attractive interactions, we find that slight modifications to either ϵns or Nfrac can result in dramatically altered self-assembled states. Using a combination of heterogeneous and homogeneous mix of polymers, we establish how this interplay between interaction energies dictates the propensity of biomolecules to find the correct binding partner at dilute concentrations in crowded environments. Author summary: Biological function relies on the ability of biomolecules to bind to specific interaction partners. In the crowded cellular milieu, the process of biomolecules binding to specific interaction partners to carry out a function is non-trivial. A mere diffusion-limited meeting of interaction partners in space does not ensure biomolecular function. Rather, even when in contact, these molecules have to find the correct orientations to carry out their function. Further, for dynamic biomolecules such as polymers, the functional configuration is often one of several possible configurations. In this scenario, assuming the functional configuration after binding involves overcoming a significant entropic barrier. Therefore, the interacting biomolecules have to dwell in contact long enough before they reorganize to find the functional orientation! While functional contacts offer an enthalpic gain, they are often only a small fraction of all protein-protein interactions. Therefore, in this paper, we study the role played by non-functional, promiscuous interactions in shaping the thermodynamics and kinetics of the formation of specific interactions. Our results suggest that there exists an optimal range of non-specific interaction strengths, which promotes the process of biomolecular complexes finding the functional configurations. [ABSTRACT FROM AUTHOR]

Microtubules are nano-machines that grow and shrink stochastically, making use of the coupling between chemical kinetics and mechanics of its constituent protofilaments (PFs). We investigate the stability and shrinkage of microtubules taking into account inter-protofilament interactions and bending interactions of intrinsically curved PFs. Computing the free energy as a function of PF tip position, we show that the competition between curvature energy, inter-PF interaction energy and entropy leads to a rich landscape with a series of minima that repeat over a length-scale determined by the intrinsic curvature. Computing Langevin dynamics of the tip through the landscape and accounting for depolymerization, we calculate the average unzippering and shrinkage velocities of GDP protofilaments and compare them with the experimentally known results. Our analysis predicts that the strength of the inter-PF interaction () has to be comparable to the strength of the curvature energy () such that , and questions the prevalent notion that unzippering results from the domination of bending energy of curved GDP PFs. Our work demonstrates how the shape of the free energy landscape is crucial in explaining the mechanism of MT shrinkage where the unzippered PFs will fluctuate in a set of partially peeled off states and subunit dissociation will reduce the length. [ABSTRACT FROM AUTHOR]

How cytoskeletal filaments collectively undergo growth and shrinkage is an intriguing question. Collective properties of multiple bio-filaments (actin or microtubules) undergoing hydrolysis have not been studied extensively earlier within simple theoretical frameworks. In this paper, we study the collective dynamical properties of multiple filaments under force, and demonstrate the distinct properties of a multi-filament system in comparison to a single filament. Comparing stochastic simulation results with recent experimental data, we show that multi-filament collective catastrophes are slower than catastrophes of single filaments. Our study also shows further distinctions as follows: (i) force-dependence of the cap-size distribution of multiple filaments are quantitatively different from that of single filaments, (ii) the diffusion constant associated with the system length fluctuations is distinct for multiple filaments, and (iii) switching dynamics of multiple filaments between capped and uncapped states and the fluctuations therein are also distinct. We build a unified picture by establishing interconnections among all these collective phenomena. Additionally, we show that the collapse times during catastrophes can be sharp indicators of collective stall forces exceeding the additive contributions of single filaments. [ABSTRACT FROM AUTHOR]

Physical properties of the extracellular matrix (ECM) are known to regulate cellular processes ranging from spreading to differentiation, with alterations in cell phenotype closely associated with changes in physical properties of cells themselves. When plated on substrates of varying stiffness, fibroblasts have been shown to exhibit stiffness matching property, wherein cell cortical stiffness increases in proportion to substrate stiffness up to 5 kPa, and subsequently saturates. Similar mechanoadaptation responses have also been observed in other cell types. Trypsin de-adhesion represents a simple experimental framework for probing the contractile mechanics of adherent cells, with de-adhesion timescales shown to scale inversely with cortical stiffness values. In this study, we combine experiments and computation in deciphering the influence of substrate properties in regulating de-adhesion dynamics of adherent cells. We first show that NIH 3T3 fibroblasts cultured on collagen-coated polyacrylamide hydrogels de-adhere faster on stiffer substrates. Using a simple computational model, we qualitatively show how substrate stiffness and cell-substrate bond breakage rate collectively influence de-adhesion timescales, and also obtain analytical expressions of de-adhesion timescales in certain regimes of the parameter space. Finally, by comparing stiffness-dependent experimental and computational de-adhesion responses, we show that faster de-adhesion on stiffer substrates arises due to force-dependent breakage of cell-matrix adhesions. In addition to illustrating the utility of employing trypsin de-adhesion as a biophysical tool for probing mechanoadaptation, our computational results highlight the collective interplay of substrate properties and bond breakage rate in setting de-adhesion timescales. [ABSTRACT FROM AUTHOR]

Why most of the in vivo experiments do not find the 30-nm chromatin fiber, well studied in vitro, is a puzzle. Two basic physical inputs that are crucial for understanding the structure of the 30-nm fiber are the stiffness of the linker DNA and the relative orientations of the DNA entering/exiting nucleosomes. Based on these inputs we simulate chromatin structure and show that the presence of non-histone proteins, which bind and locally bend linker DNA, destroys any regular higher order structures (e.g., zig-zag). Accounting for the bending geometry of proteins like nhp6 and HMG-B, our theory predicts phase-diagram for the chromatin structure as a function of DNA-bending non-histone protein density and mean linker DNA length. For a wide range of linker lengths, we show that as we vary one parameter, that is, the fraction of bent linker region due to non-histone proteins, the steady-state structure will show a transition from zig-zag to an irregular structure—a structure that is reminiscent of what is observed in experiments recently. Our theory can explain the recent in vivo observation of irregular chromatin having co-existence of finite fraction of the next-neighbor (i + 2) and neighbor (i + 1) nucleosome interactions. [ABSTRACT FROM AUTHOR]

[Display omitted] • Amyloid polymorphism is associated with distinct disease phenotypes in proteopathies. • Conformational heterogeneity during α-Syn aggregation forms the basis of polymorphism. • α-Syn forms polymorphs under the same aggregation conditions. • Polymorphs based on aggregation intermediates show distinct prion-like behavior. • α-Syn fibrils display a unique polymorphism-dependent structure-pathogenic relationship. α-Synuclein (α-Syn) amyloids in synucleinopathies are suggested to be structurally and functionally diverse, reminiscent of prion-like strains. The mechanism of how the aggregation of the same precursor protein results in the formation of fibril polymorphs remains elusive. Here, we demonstrate the structure–function relationship of two polymorphs, pre-matured fibrils (PMFs) and helix-matured fibrils (HMFs), based on α-Syn aggregation intermediates. These polymorphs display the structural differences as demonstrated by solid-state NMR and mass spectrometry studies and also possess different cellular activities such as seeding, internalization, and cell-to-cell transfer of aggregates. HMFs, with a compact core structure, exhibit low seeding potency but readily internalize and transfer from one cell to another. The less structured PMFs lack transcellular transfer ability but induce abundant α-Syn pathology and trigger the formation of aggresomes in cells. Overall, the study highlights that the conformational heterogeneity in the aggregation pathway may lead to fibril polymorphs with distinct prion-like behavior. [ABSTRACT FROM AUTHOR]

We investigate collective dynamics of branched actin networks growing against a rigid movable wall constrained by a resistive force. Computing the force velocity relations, we show that the stall force of such networks depends not only on the average number of filaments touching the wall, but also on the amount of fluctuation of the leading edge of the network. These differences arise due to differences in the network architecture, namely, distance between two adjacent branching points and the initial distance of the starting filament from the wall, with their relative magnitudes influencing the nature of the force velocity curves (convex versus concave). We also show that the introduction of branching results in nonmonotonic diffusion constant, a quantity that measures the growth in length fluctuation of the leading edge of the network, as a function of externally applied force. Together our results demonstrate how the collective dynamics of a branched network differs from that of a parallel filament network. [ABSTRACT FROM AUTHOR]

Protein aggregation and amyloid formation are known to play a role both in diseases and in biological functions. Transcription factor p53 plays a major role in tumor suppression by maintaining genomic stability. Recent studies have suggested that amyloid formation of p53 could lead to its loss of physiological function as a tumor suppressor. Here, we investigated the intrinsic amyloidogenic nature of wild-type p53 using sequence analysis. We used bioinformatics and aggregation prediction algorithms to establish the evolutionarily conserved nature of aggregation-prone sequences in wild-type p53. Further, we analyzed the amyloid forming capacity of conserved and aggregation-prone p53-derived peptides PILTIITL and YFTLQI in vitro using various biophysical techniques, including all atom molecular dynamics simulation. Finally, we probed the seeding ability of the PILTIITL peptide on p53 aggregation in vitro and in cells. Our data demonstrate the intrinsic amyloid forming ability of a sequence stretch of the p53 DNA binding domain (DBD) and its aggregation templating behavior on full-length and p53 core domain. Therefore, p53 aggregation, instigated through an amyloidogenic segment in its DBD, could be a putative driving force for p53 aggregation in vivo. [ABSTRACT FROM AUTHOR]

Collective dynamics and force generation by cytoskeletal filaments are crucial in many cellular processes. Investigating growth dynamics of a bundle of N independent cytoskeletal filaments pushing against a wall, we show that chemical switching (ATP/GTP hydrolysis) leads to a collective phenomenon that is currently unknown. Obtaining force-velocity relations for different models that capture chemical switching, we show, analytically and numerically, that the collective stall force of N filaments is greater than N times the stall force of a single filament. Employing an exactly solvable toy model, we analytically prove the above result for N = 2. We, further, numerically show the existence of this collective phenomenon, for , in realistic models (with random and sequential hydrolysis) that simulate actin and microtubule bundle growth. We make quantitative predictions for the excess forces, and argue that this collective effect is related to the non-equilibrium nature of chemical switching. [ABSTRACT FROM AUTHOR]

The storage of protein/peptide hormones within subcellular compartments and subsequent release are crucial for their native function, and hence these processes are intricately regulated in mammalian systems. Several peptide hormones were recently suggested to be stored as amyloids within endocrine secretory granules. This leads to an apparent paradox where storage requires formation of aggregates, and their function requires a supply of non-aggregated peptides on demand. The precise mechanism behind amyloid formation by these hormones and their subsequent release remain an open question. To address this, we examined aggregation and fibril reversibility of a cyclic peptide hormone somatostatin (SST)-14 using various techniques. After proving that SST gets stored as amyloid in vivo, we investigated the role of native structure in modulating its conformational dynamics and self-association by disrupting the disulfide bridge (Cys³-Cys14) in SST. Using two-dimensional NMR, we resolved the initial structure of somatostatin-14 leading to aggregation and further probed its conformational dynamics in silico. The perturbation in native structure (S-S cleavage) led to a significant increase in conformational flexibility and resulted in rapid amyloid formation. The fibrils formed by disulfide-reduced noncyclic SST possess greater resistance to denaturing conditions with decreased monomer releasing potency. MD simulations reveal marked differences in the intermolecular interactions in SST and noncyclic SST providing plausible explanation for differential aggregation and fibril reversibility observed experimentally in these structural variants. Our findings thus emphasize that subtle changes in the native structure of peptide hormone(s) could alter its conformational dynamics and amyloid formation, which might have significant implications on their reversible storage and secretion. [ABSTRACT FROM AUTHOR]

Glycosaminoglycans (GAGs) have been reported to play a significant role in amyloid formation of a wide range of proteins/peptides either associated with diseases or native biological functions. The exact mechanism by which GAGS influence amyloid formation is not clearly understood. Here, we studied two dosely related peptides, glucagon-like peptide (GLP and glucagon-like peptide 2 (GLPZ), for their amyloid formation in the presence and absence of the representative GAG heparin using various biophysical and computational approaches. We show that the aggregation and amyloid formation by these peptides follow distinct mechanisms, GLPI follows nucleation-dependent aggregation, whereas GLP2 forms amyloids without any significant lag time. Investigating the role of heparin, we also found that heparin interacts with GLP1, accelerates Its aggregation, and gets incorporated within its amyloid fibrils. In contrast, heparin neither affects the aggregation kinetics of GLP2 nor gets embedded within its fibrils. Furthermore, we found that heparin preferentially influences the stability of the GLPI fibrils over GLP2 fibrils. To understand the specific nature of the interaction of heparin with GLP1 and GLP2, we performed all-atom MD simulations. Our in silico results show that the basic-nonbasic-basic (B[c;3]motif of GLP1 (K28-G29- R30) facilitates the interaction between heparin and peptide monomers. However, the absence of such a motif in GLP2 could be the reason for a significantly lower strength of interaction between GLP2 and heparin. Our study not only helps to understand the role of hepann in inducing protein aggregation but also provides insight into the nature of heparin-protein interaction. [ABSTRACT FROM AUTHOR]

Microtubules consist of 13 protofilaments arranged in the form of a cylinder. The protofilaments are composed of longitudinally attached tubulin dimers that can exist in either a less curved state [GTP-bound tubulin (T)] or a more curved state [GDP-bound tubulin (D)]. Hydrolysis of T into D leaves the straight and laterally attached protofilaments of the microtubule in a mechanically stressed state, thus leading to their unzipping. The elastic energy in the unzipping protofilaments can be harnessed by a force transducer such as the Dam 1 -kinetochore ring complex in order to exert pulling force on chromosomes during cell division. In the present paper we develop a simple continuum model to obtain this pulling force as a function of the mechanical properties of protofilaments and the size of the Dami-kinetochore ring. We also extend this model to investigate the role played by the T subunits found at the plus end of the microtubule (the T cap) on the mechanical stability of microtubules. [ABSTRACT FROM AUTHOR]

Peptide/protein hormones could be stored as non-toxic amyloid-like structures in pituitary secretory granules. ACTH and βendorphin are two of the important peptide hormones that get co-stored in the pituitary secretory granules. Here, we study molecular interactions between ACTH and β-endorphin and their colocalization in the form of amyloid aggregates. Although ACTH is known to be a part of ACTH-β-endorphin aggregate, ACTH alone cannot aggregate into amyloid under various plausible conditions. Using all atom molecular dynamics simulation we investigate the early molecular interaction events in the ACTH-β-endorphin system, β-endorphin-only system and ACTH-only system. We find that β-endorphin and ACTH formed an interacting unit, whereas negligible interactions were observed between ACTH molecules in ACTH-only system. Our data suggest that ACTH is not only involved in interaction with β-endorphin but also enhances the stability of mixed oligomers of the entire system. [ABSTRACT FROM AUTHOR]

Molecular motors and cytoskeletal filaments work collectively most of the time under opposing forces. This opposing force may be due to cargo carried by motors or resistance coming from the cell membrane pressing against the cytoskeletal filaments. Some recent studies have shown that the collective maximum force (stall force) generated by multiple cytoskeletal filaments or molecular motors may not always be just a simple sum of the stall forces of the individual filaments or motors. To understand this excess or deficit in the collective force, we study a broad class of models of both cytoskeletal filaments and molecular motors. We argue that the stall force generated by a group of filaments or motors is additive, that is, the stall force of N number of filaments (motors) is N times the stall force of one filament (motor), when the system is reversible at stall. Conversely, we show that this additive property typically does not hold true when the system is irreversible at stall. We thus present a novel and unified understanding of the existing models exhibiting such non-addivity, and generalise our arguments by developing new models that demonstrate this phenomena. We also propose a quantity similar to thermodynamic efficiency to easily predict this deviation from stall-force additivity for filament and motor collectives. [ABSTRACT FROM AUTHOR]