Search

In most eukaryotic cells, actin filaments assemble into a shell-like actin cortex under the plasma membrane, controlling cellular morphology, mechanics, and signaling. The actin cortex is highly polymorphic, adopting diverse forms such as the ring-like structures found in podosomes, axonal rings, and immune synapses. The biophysical principles that underlie the formation of actin rings and cortices remain unknown. Using a molecular simulation platform called MEDYAN, we discovered that varying the filament treadmilling rate and myosin concentration induces a finite size phase transition in actomyosin network structures. We found that actomyosin networks condense into clusters at low treadmilling rates or high myosin concentrations but form ring-like or cortex-like structures at high treadmilling rates and low myosin concentrations. This mechanism is supported by our corroborating experiments on live T cells, which exhibit ring-like actin networks upon activation by stimulatory antibody. Upon disruption of filament treadmilling or enhancement of myosin activity, the pre-existing actin rings are disrupted into actin clusters or collapse towards the network center respectively. Our analyses suggest that the ring-like actin structure is a preferred state of low mechanical energy, which is, importantly, only reachable at sufficiently high treadmilling rates.

Histones are the dominant proteins to compact and store DNA in both Eukarya and Archaea. For a long time, histones are observed to exist in the unit of dimers but diverge into different formats such as heterodimers in Eukarya or homodimers in Archaea. Here, by studying 11 types of histone proteins, both monomers and their dimeric complexes, using multiscale molecular dynamics (MD) simulations combined with NMR and circular dichroism experiments, we confirm the widely applied “folding upon binding” mechanism of histone structures. A histone dimer appears to form the longestα2 helices followed by other shorter helices and inter-molecular tertiary structures. We report an alternative conformation, namely, the inverted non-native dimer, which has a minimum free energy state. Protein sequence analysis indicates that the inverted conformation can be attributed to a hidden head-tail sequence symmetry underlying all histone proteins. This finding strongly support previously proposed histone evolution hypotheses. Finally, we separately used the MD-based AWSEM and AI-based AlphaFold-Multimer model to predict eukaryotic histone homodimer structures and performed extensive allatom MD simulations to examine their structural stabilities. Our results suggest that eukaryotic histones can also form stable homodimers, whereas their disordered tails— the structurally asymmetrical region—may tip the balance towards the formation of heterotypic dimers.

Extensive studies of growing axons have revealed many individual components and protein interactions that guide neuronal morphogenesis. Despite this, however, we lack any clear picture of the emergent mechanism by which this nanometer-scale biochemistry generates the multimicron-scale morphology and cell biology of axon growth and guidance in vivo. To address this, we studied the downstream effects of the Abl signaling pathway using a computer simulation software (MEDYAN) that accounts for mechanochemical dynamics of active polymers. Previous studies implicate two Abl effectors, Arp2/3 and Enabled, in Abl-dependent axon guidance decisions. We now find that Abl alters actin architecture primarily by activating Arp2/3, while Enabled plays a more limited role. Our simulations show that simulations mimicking modest levels of Abl activity bear striking similarity to actin profiles obtained experimentally from live imaging of actin in wild-type axons in vivo. Using a graph theoretical filament-filament contact analysis, moreover, we find that networks mimicking hyperactivity of Abl (enhanced Arp2/3) are fragmented into smaller domains of actin that interact weakly with each other, consistent with the pattern of actin fragmentation observed upon Abl overexpression in vivo. Two perturbative simulations further confirm that high-Arp2/3 actin networks are mechanically disconnected and fail to mount a cohesive response to perturbation. Taken together, these data provide a molecular-level picture of how the large-scale organization of the axonal cytoskeleton arises from the biophysics of actin networks.

Mechanochemical simulations of actomyosin networks are traditionally based on one-dimensional models of actin filaments having zero width. Here, and in the follow up paper, approaches are presented for more efficient modelling which incorporates stretching, bending, shearing and twisting of actin filaments. Our modelling of a semi-flexible filament with a small but finite width is based on the Cosserat theory of elastic rods, which allows for six degrees of freedom at every point on the filament's backbone. In the variational models presented in this paper, a small and discrete set of parameters is used to describe a smooth filament shape having all degrees of freedom allowed in the Cosserat theory. Two main approaches are introduced: one where polynomial spline functions describe the filament's configuration, and one in which geodesic curves in the space of the configurational degrees of freedom are used. We find that in the latter representation the strain energy function can be calculated without resorting to a small-angle expansion, so it can describe arbitrarily large filament deformations without systematic error. These approaches are validated by a dynamical model of a Cosserat filament, which can be further extended by using multi-resolution methods to allow more detailed monomer-based resolution in certain parts of the actin filament, as introduced in the follow up paper. The presented framework is illustrated by showing how torsional compliance in a finite-width filament can induce broken chiral symmetry in the structure of a cross-linked bundle., 49 pages, 9 figures

Actin networks rely on nucleation mechanisms to generate new filaments because spontaneous nucleation is kinetically disfavored. Branching nucleation of actin filaments by actin-related protein (Arp2/3), in particular, is critical for actin self-organization. In this study, we use the simulation platform for active matter MEDYAN to generate 2000 s long stochastic trajectories of actin networks, under varying Arp2/3 concentrations, in reaction volumes of biologically meaningful size (20 μm

Extensive studies of growing axons have revealed many individual components and protein interactions that guide neuronal morphogenesis. Despite this, however, we lack any clear picture of the emergent mechanism by which this nanometer-scale biochemistry generates the multi-micron scale morphology and cell biology of axon growth and guidance in vivo. To address this, we studied the downstream effects of the Abl signaling pathway using a computer simulation software (MEDYAN) that accounts for mechanochemical dynamics of active polymers. Previous studies implicate two Abl effectors, Arp2/3 and Enabled, in Abl-dependent axon guidance decisions. We now find that Abl alters actin architecture primarily by activating Arp2/3, while Enabled plays a more limited role. Our simulations show that simulations mimicking modest levels of Abl activity bear striking similarity to actin profiles obtained experimentally from live-imaging of actin in wild type axons in vivo. Using a graph-theoretical filament-filament contact analysis, moreover, we find that networks mimicking hyperactivity of Abl (enhanced Arp2/3) are fragmented into smaller domains of actin that interact weakly with each other, consistent with the pattern of actin fragmentation observed upon Abl overexpression in vivo. Two perturbative simulations further confirm that high Arp2/3 actin networks are mechanically disconnected and fail to mount a cohesive response to perturbation. Taken together, these data provide a molecular-level picture of how the large-scale organization of the axonal cytoskeleton arises from the biophysics of actin networks.Highlight summaryHow do single-molecule dynamics produce multi-micron scale changes in actin organization in an extending axon? Comparison of computational simulations to in vivo data suggests that Abl kinase and Arp2/3 expand actomyosin networks by fragmenting into multiple domains, thus toggling the axon between states of local vs global internal connectivity.

Eukaryotic cells are mechanically supported by a polymer network called the cytoskeleton, which consumes chemical energy to dynamically remodel its structure. Recent experiments in vivo have revealed that this remodeling occasionally happens through anomalously large displacements, reminiscent of earthquakes or avalanches. These cytoskeletal avalanches might indicate that the cytoskeleton's structural response to a changing cellular environment is highly sensitive, and they are therefore of significant biological interest. However, the physics underlying "cytoquakes" is poorly understood. Here, we use agent-based simulations of cytoskeletal self-organization to study fluctuations in the network's mechanical energy. We robustly observe non-Gaussian statistics and asymmetrically large rates of energy release compared to accumulation in a minimal cytoskeletal model. The large events of energy release are found to correlate with large, collective displacements of the cytoskeletal filaments. We also find that the changes in the localization of tension and the projections of the network motion onto the vibrational normal modes are asymmetrically distributed for energy release and accumulation. These results imply an avalanche-like process of slow energy storage punctuated by fast, large events of energy release involving a collective network rearrangement. We further show that mechanical instability precedes cytoquake occurrence through a machine learning model that dynamically forecasts cytoquakes using the vibrational spectrum as input. Our results provide the first connection between the cytoquake phenomenon and the network's mechanical energy and can help guide future investigations of the cytoskeleton's structural susceptibility., 35 pages, 18 figures

The plasma membrane defines the shape of the cell and plays an indispensable role in bridging intra- and extra-cellular environments. Mechanochemical interactions between plasma membrane and cytoskeleton are vital for cell biomechanics and mechanosensing. A computational model that comprehensively captures the complex, cell-scale cytoskeleton-membrane dynamics is still lacking. In this work, we introduce a triangulated membrane model that accounts for membrane’s elastic properties, as well as for membrane-filament steric interactions. The corresponding force-field was incorporated into the active biological matter simulation platform, MEDYAN (“Mechanochemical Dynamics of Active Networks”). Simulations using the new model shed light on how actin filament bundling affects generation of tubular membrane protrusions. In particular, we used membrane-MEDYAN simulations to investigate protrusion initiation and dynamics while varying geometries of filament bundles, membrane rigidities and local G-Actin concentrations. We found that bundles’ protrusion propensities sensitively depend on the synergy between bundle thickness and inclination angle at which the bundle approaches the membrane. The new model paves the way for simulations of biological systems involving intricate membrane-cytoskeleton interactions, such as occurring at the leading edge and the cortex, eventually helping to uncover the fundamental principles underlying the active matter organization in the vicinity of the membrane.

Microsecond timescale explicit-solvent atomistic simulations were carried out to investigate how anionic surfactants modulate protein structure and dynamics. We found that lysozyme undergoes near-complete denaturation at the high concentration (> 0.1 M) of sodium pentadecyl sulfonate (SPDS), while only partial denaturation occurs at the concentration slightly below 0.1 M. In large part, protein denaturation is structurally manifested by disappearance of helical segments and loss of tertiary interactions. The computational prediction of the extent of burial of cysteine residues was experimentally validated by measuring the accessibility of the respective sulfhydryl groups. Overall, our work indicates an interesting synergy between electrostatic and hydrophobic contributions to lysozyme’s denaturation process by anionic surfactants. In fact, first disulfide bridges and hydrogen bonds from protein surface to SPDS head groups loosen the protein globule followed by fuller denaturation via insertion of the surfactant’s hydrophobic tails into the protein core.

Actin networks are highly dynamic cytoskeletal structures that continuously undergo structural remodeling. One prominent way to probe these processes is via Fluorescence Recovery After Photobleaching (FRAP), which can be used to estimate the rate of turnover for filamentous actin monomers. It is thought that head-to-tail treadmilling and de novo filament nucleation constitute two primary mechanisms underlying turnover kinetics. More generally, these self-assembly activities are responsible for many important cellular functions such as force generation, cellular shape dynamics, and cellular motility. In what relative proportions filament treadmilling and de novo filament nucleation contribute to actin network turnover is still not fully understood. We used an advanced stochastic reaction-diffusion model in three dimensions, MEDYAN, to study turnover dynamics of actin networks containing Arp2/3, formin and capping protein at experimentally meaningful length- and time-scales. Our results reveal that, most commonly, treadmilling of older filaments is the main contributor to actin network turnover. On the other hand, although turnover and treadmilling are often used interchangeably, we show clear instances where this assumption would not be justified, for example, finding that rapid turnover is accompanied by slow treadmilling in highly dendritic Arp2/3 networks.

In most eukaryotic cells, actin filaments assemble into a shell-like actin cortex under the plasma membrane, controlling cellular morphology, mechanics, and signaling. The actin cortex is highly polymorphic, adopting diverse forms such as the ring-like structures found in podosomes, axonal rings, and immune synapses. The biophysical principles that underlie the formation of actin rings and cortices remain unknown. Using a molecular simulation platform, called MEDYAN, we discovered that varying the filament treadmilling rate and myosin concentration induces a finite size phase transition in actomyosin network structures. We found that actomyosin networks condense into clusters at low treadmilling rates or high myosin concentration but form ring-like or cortex-like structures at high treadmilling rates and low myosin concentration. This mechanism is supported by our corroborating experiments on live T cells, which exhibit ring-like actin networks upon activation by stimulatory antibody. Upon disruption of filament treadmilling or enhancement of myosin activity, the pre-existing actin rings are disrupted into actin clusters or collapse towards the network center respectively. Our analyses suggest that the ring-like actin structure is a preferred state of low mechanical energy, which is, importantly, only reachable at sufficiently high treadmilling rates.

Microsecond timescale explicit-solvent atomistic simulations were carried out to investigate how anionic surfactants modulate protein structure and dynamics. We found that lysozyme undergoes near-complete denaturation at the high concentration ( 0.1 M) of sodium pentadecyl sulfonate (SPDS), while only partial denaturation occurs at the concentration slightly below 0.1 M. In large part, protein denaturation is structurally manifested by disappearance of helical segments and loss of tertiary interactions. The computational prediction of the extent of burial of cysteine residues was experimentally validated by measuring the accessibility of the respective sulfhydryl groups. Overall, our work indicates an interesting synergy between electrostatic and hydrophobic contributions to lysozyme's denaturation process by anionic surfactants. In fact, first disulfide bridges and hydrogen bonds from protein surface to SPDS head groups loosen the protein globule followed by fuller denaturation via insertion of the surfactant's hydrophobic tails into the protein core.

β- and γ-cytoplasmic actins are ubiquitously expressed in every cell type and are nearly identical at the amino acid level but play vastly different roles in vivo. Their essential roles in embryogenesis and mesenchymal cell migration critically depend on the nucleotide sequences of their genes, rather than their amino acid sequences; however, it is unclear which gene elements underlie this effect. Here we address the specific role of the coding sequence in β- and γ-cytoplasmic actins’ intracellular functions, using stable polyclonal populations of immortalized mouse embryonic fibroblasts with exogenously expressed actin isoforms and their ‘codon-switched’ variants. When targeted to the cell periphery using β-actin 3′UTR; β-actin and γ-actin have differential effects on cell migration. These effects directly depend on the coding sequence. Single-molecule measurements of actin isoform translation, combined with fluorescence recovery after photobleaching, demonstrate a pronounced difference in β- and γ-actins’ translation elongation rates in cells, leading to changes in their dynamics at focal adhesions, impairments in actin bundle formation, and reduced cell anchoring to the substrate during migration. Our results demonstrate that coding sequence-mediated differences in actin translation play a key role in cell migration., eLife digest Most mammalian cells make both β- and γ-actin, two proteins which shape the cell’s internal skeleton and its ability to migrate. The molecules share over 99% of their sequence, yet they play distinct roles. In fact, deleting the β-actin gene in mice causes death in the womb, while the animals can survive with comparatively milder issues without their γ-actin gene. How two similar proteins can have such different biological roles is a long-standing mystery. A closer look could hold some clues: β- and γ-actin may contain the same blocks (or amino acids), but the genetic sequences that encode these proteins differ by about 13%. This is because different units of genetic information – known as synonymous codons – can encode the same amino acid. These ‘silent substitutions’ have no effect on the sequence of the proteins, yet a cell reads synonymous codons (and therefore produces proteins) at different speeds. To find out the impact of silent substitutions, Vedula et al. swapped the codons for the two proteins, forcing mouse cells to produce β-actin using γ-actin codons, and vice versa. Cells with non-manipulated γ-actin and those with β-actin made using γ-actin codons could move much faster than cells with β-actin. This suggested that silent substitutions were indeed affecting the role of the protein. Vedula et al. found that cells read γ-codons – and therefore made γ-actin – much more slowly than β-codons: this also affected how quickly the protein could be dispatched where it was needed in the cell. Slower production meant that bundles of γ-actin were shorter, which allowed cells to move faster by providing a weaker anchoring system. Overall, this work provides new links between silent substitutions and protein behavior, a relatively new research area which is likely to shed light on other protein families.

Simulating soft matter systems such as the cytoskeleton can enable deep understanding of experimentally observed phenomena. One challenge of modeling such systems is realistic description of the steric repulsion between nearby polymers. Previous models of the polymeric excluded volume interaction have the deficit of being non-analytic, being computationally expensive, or allowing polymers to erroneously cross each other. A recent solution to these issues, implemented in the MEDYAN simulation platform, uses analytical expressions obtained from integrating an interaction kernel along the lengths of two polymer segments to describe their repulsion. Here, we extend this model by re-deriving it for lower-dimensional geometrical configurations, deriving similar expressions using a steeper interaction kernel, comparing it to other commonly used potentials, and showing how to parameterize these models. We also generalize this new integrated style of potential by introducing a segmental Lennard-Jones potential, which enables modelling both attractive and repulsive interactions in semi-flexible polymer networks. These results can be further generalized to facilitate the development of effective interaction potentials for other finite elements in simulations of soft-matter systems., Comment: 18 pages, 6 figures

The accurate and reliable prediction of the 3D structures of proteins and their assemblies remains difficult even though the number of solved structures soars and prediction techniques improve. In this study, a free and open access web server, AWSEM-Suite, whose goal is to predict monomeric protein tertiary structures from sequence is described. The model underlying the server’s predictions is a coarse-grained protein force field which has its roots in neural network ideas that has been optimized using energy landscape theory. Employing physically motivated potentials and knowledge-based local structure biasing terms, the addition of homologous template and co-evolutionary restraints to AWSEM-Suite greatly improves the predictive power of pure AWSEM structure prediction. From the independent evaluation metrics released in the CASP13 experiment, AWSEM-Suite proves to be a reasonably accurate algorithm for free modeling, standing at the eighth position in the free modeling category of CASP13. The AWSEM-Suite server also features a front end with a user-friendly interface. The AWSEM-Suite server is a powerful tool for predicting monomeric protein tertiary structures that is most useful when a suitable structure template is not available. The AWSEM-Suite server is freely available at: https://awsem.rice.edu.

In modeling the interior of cells by simulating a reaction–diffusion master equation over a grid of compartments, one employs the assumption that the copy numbers of various chemical species are small, discrete quantities. We show that, in this case, textbook expressions for the change in Gibbs free energy accompanying a chemical reaction or diffusion between adjacent compartments are inaccurate. We derive exact expressions for these free energy changes for the case of discrete copy numbers and show how these expressions reduce to traditional expressions under a series of successive approximations leveraging the relative sizes of the stoichiometric coefficients and the copy numbers of the solutes and solvent. Numerical results are presented to corroborate the claim that if the copy numbers are treated as discrete quantities, then only these more accurate expressions lead to correct behavior. Thus, the newly derived expressions are critical for correctly computing entropy production in mesoscopic simulations based on the reaction–diffusion master equation formalism.

Understanding cellular remodeling in response to mechanical stimuli is a critical step in elucidating mechano-activation of biochemical signaling pathways. Experimental evidence indicates that external stress-induced subcellular adaptation is accomplished through dynamic cytoskeletal reorganization. To study the interactions between subcellular structures involved in transducing mechanical signals, we combined experimental and computational simulations to evaluate real-time mechanical adaptation of the actin cytoskeletal network. Actin cytoskeleton was imaged at the same time as an external tensile force was applied to live vascular smooth muscle cells using a fibronectin-functionalized atomic force microscope probe. In addition, we performed computational simulations of active cytoskeletal networks under a tensile external force. The experimental data and simulation results suggest that mechanical structural adaptation occurs before chemical adaptation during filament bundle formation: actin filaments first align in the direction of the external force, initializing anisotropic filament orientations, then the chemical evolution of the network follows the anisotropic structures to further develop the bundle-like geometry. This finding presents an alternative, novel explanation for the stress fiber formation and provides new insight into the mechanism of mechanotransduction.Author SummaryRemodeling the cytoskeletal network in response to external force is key to mechanosensing and locomotion. Despite much focus on cytoskeletal remodeling in recent years, a comprehensive understanding of actin remodeling in real-time in cells under mechanical stimuli is still lacking. We integrated stress-induced 3D actin imaging and 3D computational simulations of actin cytoskeleton to study how the actin cytoskeleton form bundles and how these bundles evolve over time upon external tensile stress. We found a rapid actin alignment and a slower bundle evolution leading to denser bundles. Based on these results, we propose a “mechanics before chemistry” model of actin cytoskeleton remodeling under external force.

A long time ago, Kuhn predicted that long polymers should approach a limit where their global motion is controlled by solvent friction alone, with ruggedness of their energy landscapes having no consequences for their dynamics. In contrast, internal friction effects are important for polymers of modest length. Internal friction in proteins, in particular, affects how fast they fold or find their binding targets and, as such, has attracted much recent attention. Here we explore the molecular origins of internal friction in unfolded proteins using atomistic simulations, coarse-grained models and analytic theory. We show that the characteristic internal friction timescale is directly proportional to the timescale of hindered dihedral rotations within the polypeptide chain, with a proportionality coefficient b that is independent of the chain length. Such chain length independence of b provides experimentally testable evidence that internal friction arises from concerted, crankshaft-like dihedral rearrangements. In accord with phenomenological models of internal friction, we find the global reconfiguration timescale of a polypeptide to be the sum of solvent friction and internal friction timescales. At the same time, the time evolution of inter-monomer distances within polypeptides deviates both from the predictions of those models and from a simple, one-dimensional diffusion model.

Significance Nucleosomes are the base units that organize eukaryotic genomes. Besides the canonical histone, histone variants create unique local chromatin domains that fine-tune transcription, replication, DNA damage repair, and faithful chromosome segregation. We developed computational and single-molecule nanoindentation tools to determine mechanical properties of histone variant nucleosomes. We found that the CENP-A nucleosome variant is more elastic than the canonical H3 nucleosome but becomes stiffer when bound to its partner CENP-C. In addition, CENP-C induces cross-array clustering, creating a chromatin state that is less accessible. These data suggest that innate material properties of nucleosomes can influence the ultimate chromatin state, and thereby influence biological outcomes., Histone variants fine-tune transcription, replication, DNA damage repair, and faithful chromosome segregation. Whether and how nucleosome variants encode unique mechanical properties to their cognate chromatin structures remains elusive. Here, using in silico and in vitro nanoindentation methods, extending to in vivo dissections, we report that histone variant nucleosomes are intrinsically more elastic than their canonical counterparts. Furthermore, binding proteins, which discriminate between histone variant nucleosomes, suppress this innate elasticity and also compact chromatin. Interestingly, when we overexpress the binding proteins in vivo, we also observe increased compaction of chromatin enriched for histone variant nucleosomes, correlating with diminished access. Taken together, these data suggest a plausible link between innate mechanical properties possessed by histone variant nucleosomes, the adaptability of chromatin states in vivo, and the epigenetic plasticity of the underlying locus.

Histone variants regulate replication, transcription, DNA damage repair, and chromosome segregation. Though widely accepted as a paradigm, it has not been rigorously demonstrated that histone variants encode unique mechanical properties. Here, we present a new theoretical approach called Minimal Cylinder Analysis (MCA) to determine the Young's modulus of nucleosomes from all-atom Molecular Dynamics (MD) simulations. Recently, we validated this computational analysis against in vitro single-molecule nanoindentation of histone variant nucleosomes. In this report, we further extend MCA to study the biophysical properties of hybrid nucleosomes that are known to exist in human cancer cells and contain H3 histone variants CENP-A and H3.3. We investigate the mechanism by which the elasticity of these heterotypic nucleosomes augments cryptic binding surfaces. Further, we derive biological predictions that might arise when such heterotypic nucleosomes take over large parts of the genome., 37 pages, 3 figures

In eukaryotes, DNA is packaged within nucleosomes. The DNA of each nucleosome is typically centered around an octameric histone protein core: one central tetramer plus two separate dimers. Studying the assembly mechanisms of histones is essential for understanding the dynamics of entire nucleosomes and higher-order DNA packaging. Here we investigate canonical histone assembly and that of the centromere-specific histone variant CENP-A using molecular dynamics simulations. We quantitatively characterize their thermodynamical and dynamical features, showing that two H3/H4 dimers form a structurally floppy, weakly-bound complex, the latter exhibiting large instability around the central interface manifested via a swiveling motion of two halves. This finding is consistent with the recently observed DNA handedness flipping of the tetrasome. In contrast, the variant CENP-A encodes distinctive stability to its tetramer with a rigid but twisted interface compared to the crystal structure, implying diverse structural possibilities of the histone variant. Interestingly, the observed tetramer dynamics alter significantly and appear to reach a new balance when H2A/H2B dimers are present. Furthermore, we found that the preferred structure for the (CENP-A/H4)2 tetramer is incongruent with the octameric structure, explaining many of the unusual dynamical behaviors of the CENP-A nucleosome. In all, these data reveal key mechanistic insights and structural details for the assembly of canonical and variant histone tetramers and octamers, providing theoretical quantifications and physical interpretations for longstanding and recent experimental observations. Based on these findings, we propose different chaperone-assisted binding and nucleosome assembly mechanisms for the canonical and CENP-A histone oligomers.

Linker histone H1 is an essential regulatory protein for many critical biological processes, such as eukaryotic chromatin packaging and gene expression. Mis-regulation of H1s is commonly observed in tumor cells, where the balance between different H1 subtypes has been shown to alter the cancer phenotype. Consisting of a rigid globular domain and two highly charged terminal domains, H1 can bind to multiple sites on a nucleosomal particle to alter chromatin hierarchical condensation levels. In particular, the disordered H1 amino- and carboxyl-terminal domains (NTD/CTD) are believed to enhance this binding affinity, but their detailed dynamics and functions remain unclear. In this work, we used a coarse-grained computational model AWSEM-DNA to simulate the H1.0b-nucleosome complex, namely chromatosome. Our results demonstrate that H1 disordered domains restrict the dynamics of both globular H1 and linker DNA arms, resulting in a more compact and rigid chromatosome particle. Furthermore, we identified regions of H1 disordered domains that are tightly tethered to DNA near the entry-exit site. Overall, our study elucidates at near atomic resolution the way the disordered linker histone H1 modulates nucleosome’s structural preferences and conformational dynamics.

Histone proteins are essential for the organization, expression, and inheritance of genetic material for eukaryotic cells. A centromere-specific H3 histone variant, centromere protein A (CENP-A), shares about 50% amino acid sequence identity with H3. CENP-A is required for packaging the centromere and for the proper separation of chromosomes during mitosis. Despite their distinct biological functions, previously reported crystal structures of the CENP-A/H4 and H3/H4 dimers reveal a high degree of similarity. In this work, we characterize the structural dynamics of CENP-A/H4 and H3/H4 dimers based on a dual-resolution approach, using both microsecond-scale explicit-solvent all-atom and coarse-grained (CG) molecular dynamics (MD) simulations. Our data show that the H4 histone is significantly more rigid compared with the H3 histone and its variant CENP-A, hence, serving as a reinforcing structural element within the histone core. We report that the CENP-A/H4 dimer is significantly more dynamic than its canonical counterpart H3/H4, and our results provide a physical explanation for this flexibility. Further, we observe that the centromere-specific chaperone Holliday Junction Recognition Protein (HJURP) stabilizes the CENP-A/H4 dimer by forming a specific electrostatic interaction network. Finally, replacing CENP-A S68 with E68 disrupts the binding interface between CENP-A and HJURP in all-atom MD simulation, and consistently, in vivo experiments demonstrate that replacing CENP-A S68 with E68 disrupts CENP-A’s localization to the centromere. Based on all our results, we propose that, during the CENP-A/H4 deposition process, the chaperone HJURP protects various substructures of the dimer, serving both as a folding and binding chaperone.

The nucleosome is the fundamental unit for packaging the genome. A detailed molecular picture for its conformational dynamics is crucial for understanding transcription and gene regulation. We investigate the disassembly of single nucleosomes using a predictive coarse-grained protein DNA model with transferable force fields. This model quantitatively describes the thermodynamic stability of both the histone core complex and the nucleosome and predicts rates of transient nucleosome opening that match experimental measurements. Quantitative characterization of the free-energy landscapes reveals the mechanism of nucleosome unfolding in which DNA unwinding and histone protein disassembly are coupled. The interfaces between H2A-H2B dimers and the (H3-H4)2 tetramer are first lost when the nucleosome opens releasing a large fraction but not all of its bound DNA. For the short strands studied in single molecule experiments, the DNA unwinds asymmetrically from the histone proteins, with only one of its two ends preferentially exposed. The detailed molecular mechanism revealed in this work provides a structural basis for interpreting experimental studies of nucleosome unfolding.

Although protein folding reactions are usually studied under static external conditions, it is likely that proteins fold in a locally fluctuating cellular environment in vivo. To mimic such behavior in in vitro experiments, the local temperature of the solvent can be modulated either harmonically or using correlated noise. In this study, coarse-grained molecular simulations are used to investigate these possibilities, and it is found that both periodic and correlated random fluctuations of the environment can indeed accelerate folding kinetics if the characteristic frequencies of the applied fluctuations are commensurate with the internal timescale of the folding reaction; this is consistent with the phenomenon of stochastic resonance observed in many other condensed-matter processes. To test this theoretical prediction, the folding dynamics of phosphoglycerate kinase under harmonic temperature fluctuations are experimentally probed using Förster resonance energy transfer fluorescence measurements. To analyze these experiments, a combination of theoretical approaches is developed, including stochastic simulations of folding kinetics and an analytical mean-field kinetic theory. The experimental observations are consistent with the theoretical predictions of stochastic resonance in phosphoglycerate kinase folding. When combined with an alternative experiment on the protein VlsE using a power spectrum analysis, elaborated in Dave et al., ChemPhysChem 2016, 10.1002/cphc.201501041, the overall data overwhelmingly point to the experimental confirmation of stochastic resonance in protein folding dynamics.

Histone variants fine-tune transcription, replication, DNA damage repair, and faithful chromosome segregation. Whether and how nucleosome variants encode unique mechanical properties to their cognate chromatin structures remains elusive. Here, using novelin silicoandin vitronanoindentation methods, extending toin vivodissections, we report that histone variant nucleosomes are intrinsically more elastic than their canonical counterparts. Furthermore, binding proteins which discriminate between histone variant nucleosomes suppress this innate elasticity and also compact chromatin. Interestingly, when we overexpress the binding proteinsin vivo, we also observe increased compaction of chromatin enriched for histone variant nucleosomes, correlating with diminished access. Together, these data suggest a plausible link between innate mechanical properties possessed by histone variant nucleosomes, the adaptability of chromatin statesin vivo, and the epigenetic plasticity of the underlying locus.SignificanceNucleosomes are the base unit which organize eukaryotic genomes. Besides the canonical histone, histone variants create unique local chromatin domains that fine-tune transcription, replication, DNA damage repair, and faithful chromosome segregation. We developed computational and single-molecule nanoindentation tools to determine mechanical properties of histone variant nucleosomes. We found that the CENP-A nucleosome variant is more elastic than the canonical H3 nucleosome but becomes stiffer when bound to its partner CENP-C. In addition, CENP-C induces cross-array clustering, creating a chromatin state that less accessible. These data suggest that innate material properties of nucleosomes can influence the ultimate chromatin state, thereby influence biological outcomes.

The associative memory, water-mediated, structure and energy model (AWSEM) has been successfully used to study protein folding, binding, and aggregation problems. In this work, we introduce AWSEM-IDP, a new AWSEM branch for simulating intrinsically disordered proteins (IDPs), where the weights of the potentials determining secondary structure formation have been finely tuned, and a novel potential is introduced that helps to precisely control both the average extent of protein chain collapse and the chain's fluctuations in size. AWSEM-IDP can efficiently sample large conformational spaces, while retaining sufficient molecular accuracy to realistically model proteins. We applied this new model to two IDPs, demonstrating that AWSEM-IDP can reasonably well reproduce higher-resolution reference data, thus providing the foundation for a transferable IDP force field. Finally, we used thermodynamic perturbation theory to show that, in general, the conformational ensembles of IDPs are highly sensitive to fine-tuning of force field parameters.

Current understanding of how contractility emerges in disordered actomyosin networks of non-muscle cells is still largely based on the intuition derived from earlier works on muscle contractility. This view, however, largely overlooks the free energy gain following passive cross-linker binding, which, even in the absence of active fluctuations, provides a thermodynamic drive towards highly overlapping filamentous states. In this work, we shed light on this phenomenon, showing that passive cross-linkers, when considered in the context of two anti-parallel filaments, generate noticeable contractile forces. However, as binding free energy of cross-linkers is increased, a sharp onset of kinetic arrest follows, greatly diminishing effectiveness of this contractility mechanism, allowing the network to contract only with weakly resisting tensions at its boundary. We have carried out stochastic simulations elucidating this mechanism, followed by a mean-field treatment that predicts how contractile forces asymptotically scale at small and large binding energies, respectively. Furthermore, when considering an active contractile filament pair, based on non-muscle myosin II, we found that the non-processive nature of these motors leads to highly inefficient force generation, due to recoil slippage of the overlap during periods when the motor is dissociated. However, we discovered that passive cross-linkers can serve as a structural ratchet during these unbound motor time spans, resulting in vast force amplification. Our results shed light on the non-equilibrium effects of transiently binding proteins in biological active matter, as observed in the non-muscle actin cytoskeleton, showing that highly efficient contractile force dipoles result from synergy of passive cross-linker and active motor dynamics, via a ratcheting mechanism on a funneled energy landscape., 13 pages, 6 figures