The role of cryptic ancestral symmetry in histone folding mechanisms across Eukarya and Archaea.

Histones compact and store DNA in both Eukarya and Archaea, forming heterodimers in Eukarya and homodimers in Archaea. Despite this, the folding mechanism of histones across species remains unclear. Our study addresses this gap by investigating 11 types of histone and histone-like proteins across humans, Drosophila, and Archaea through multiscale molecular dynamics (MD) simulations, complemented by NMR and circular dichroism experiments. We confirm and elaborate on the widely applied "folding upon binding" mechanism of histone dimeric proteins and report a new alternative conformation, namely, the inverted non-native dimer, which may be a thermodynamically metastable configuration. Protein sequence analysis indicated that the inverted conformation arises from the hidden ancestral head-tail sequence symmetry underlying all histone proteins, which is congruent with the previously proposed histone evolution hypotheses. Finally, to explore the potential formations of homodimers in Eukarya, we utilized MD-based AWSEM and AI-based AlphaFold-Multimer models to predict their structures and conducted extensive all-atom MD simulations to examine their respective structural stabilities. Our results suggest that eukaryotic histones may also form stable homodimers, whereas their disordered tails bring significant structural asymmetry and tip the balance towards the formation of commonly observed heterotypic dimers.

Myosin and [Formula: see text]-actinin regulation of stress fiber contractility under tensile stress.

Stress fibers are actomyosin bundles that regulate cellular mechanosensation and force transduction. Interacting with the extracellular matrix through focal adhesion complexes, stress fibers are highly dynamic structures regulated by myosin motors and crosslinking proteins. Under external mechanical stimuli such as tensile forces, the stress fiber remodels its architecture to adapt to external cues, displaying properties of viscoelastic materials. How the structural remodeling of stress fibers is related to the generation of contractile force is not well understood. In this work, we simulate mechanochemical dynamics and force generation of stress fibers using the molecular simulation platform MEDYAN. We model stress fiber as two connecting bipolar bundles attached at the ends to focal adhesion complexes. The simulated stress fibers generate contractile force that is regulated by myosin motors and [Formula: see text]-actinin crosslinkers. We find that stress fibers enhance contractility by reducing the distance between actin filaments to increase crosslinker binding, and this structural remodeling ability depends on the crosslinker turnover rate. Under tensile pulling force, the stress fiber shows an instantaneous increase of the contractile forces followed by a slow relaxation into a new steady state. While the new steady state contractility after pulling depends only on the overlap between actin bundles, the short-term contractility enhancement is sensitive to the tensile pulling distance. We further show that this mechanical response is also sensitive to the crosslinker turnover rate. Our results provide new insights into the stress fiber mechanics that have significant implications for understanding cellular adaptation to mechanical signaling.

Enabled primarily controls filopodial morphology, not actin organization, in the TSM1 growth cone in Drosophila.

Ena/VASP proteins are processive actin polymerases that are required throughout animal phylogeny for many morphogenetic processes, including axon growth and guidance. Here we use in vivo live imaging of morphology and actin distribution to determine the role of Ena in promoting the growth of the TSM1 axon of the Drosophila wing. Altering Ena activity causes stalling and misrouting of TSM1. Our data show that Ena has a substantial impact on filopodial morphology in this growth cone but exerts only modest effects on actin distribution. This is in contrast to the main regulator of Ena, Abl tyrosine kinase, which was shown previously to have profound effects on actin and only mild effects on TSM1 growth cone morphology. We interpret these data as suggesting that the primary role of Ena in this axon may be to link actin to the morphogenetic processes of the plasma membrane, rather than to regulate actin organization itself. These data also suggest that a key role of Ena, acting downstream of Abl, may be to maintain consistent organization and reliable evolution of growth cone structure, even as Abl activity varies in response to guidance cues in the environment.

A tug of war between filament treadmilling and myosin induced contractility generates actin rings.

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.

Nucleation causes an actin network to fragment into multiple high-density domains.

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 µm3). We find that the dynamics of Arp2/3 increase the abundance of short filaments and increases network treadmilling rate. By analyzing the density fields of F-actin, we find that at low Arp2/3 concentrations, F-actin is organized into a single connected and contractile domain, while at elevated Arp2/3 levels (10 nM and above), such high-density actin domains fragment into smaller domains spanning a wide range of volumes. These fragmented domains are extremely dynamic, continuously merging and splitting, owing to the high treadmilling rate of the underlying actin network. Treating the domain dynamics as a drift-diffusion process, we find that the fragmented state is stochastically favored, and the network state slowly drifts toward the fragmented state with considerable diffusion (variability) in the number of domains. We suggest that tuning the Arp2/3 concentration enables cells to transition from a globally coherent cytoskeleton, whose response involves the entire cytoplasmic network, to a fragmented cytoskeleton, where domains can respond independently to locally varying signals.

Computational simulations reveal that Abl activity controls cohesiveness of actin networks in growth cones.

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.

On Stretching, Bending, Shearing, and Twisting of Actin Filaments I: Variational Models.

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 (arXiv, DOI 10.48550/arXiv.2203.01284), approaches are presented for more efficient modeling that incorporates stretching, shearing, and twisting of actin filaments. Our modeling of a semiflexible 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 multiresolution 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.

Simulated actin reorganization mediated by motor proteins.

Cortical actin networks are highly dynamic and play critical roles in shaping the mechanical properties of cells. The actin cytoskeleton undergoes significant reorganization in many different contexts, including during directed cell migration and over the course of the cell cycle, when cortical actin can transition between different configurations such as open patched meshworks, homogeneous distributions, and aligned bundles. Several types of myosin motor proteins, characterized by different kinetic parameters, have been involved in this reorganization of actin filaments. Given the limitations in studying the interactions of actin with myosin in vivo, we propose stochastic agent-based models and develop a set of data analysis measures to assess how myosin motor proteins mediate various actin organizations. In particular, we identify individual motor parameters, such as motor binding rate and step size, that generate actin networks with different levels of contractility and different patterns of myosin motor localization, which have previously been observed experimentally. In simulations where two motor populations with distinct kinetic parameters interact with the same actin network, we find that motors may act in a complementary way, by tuning the actin network organization, or in an antagonistic way, where one motor emerges as dominant. This modeling and data analysis framework also uncovers parameter regimes where spatial segregation between motor populations is achieved. By allowing for changes in kinetic rates during the actin-myosin dynamic simulations, our work suggests that certain actin-myosin organizations may require additional regulation beyond mediation by motor proteins in order to reconfigure the cytoskeleton network on experimentally-observed timescales.

Understanding cytoskeletal avalanches using mechanical stability analysis.

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 a connection between the cytoquake phenomenon and the network's mechanical energy and can help guide future investigations of the cytoskeleton's structural susceptibility.

Membrane-MEDYAN: Simulating Deformable Vesicles Containing Complex Cytoskeletal Networks.

The plasma membrane defines the shape of the cell and plays an indispensable role in bridging intra- and extracellular 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 the 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 the 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 those 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.

A molecular dynamics study of protein denaturation induced by sulfonate-based surfactants.

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.

Binding Dynamics of Disordered Linker Histone H1 with a Nucleosomal Particle.

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 and conformation 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.

Tensile force-induced cytoskeletal remodeling: Mechanics before chemistry.

Understanding cellular remodeling in response to mechanical stimuli is a critical step in elucidating mechanical 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 data 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. Moreover, we performed computational simulations of active cytoskeletal networks under an external tensile 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 by initializing anisotropic filament orientations, then the chemical evolution of the network follows the anisotropic structures to further develop the bundle-like geometry. Our findings present an alternative two-step explanation for the formation of actin bundles due to mechanical stimulation and provide new insights into the mechanism of mechanotransduction.

AWSEM-Suite: a protein structure prediction server based on template-guided, coevolutionary-enhanced optimized folding landscapes.

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.

Gibbs free energy change of a discrete chemical reaction event.

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.

Minimal Cylinder Analysis Reveals the Mechanical Properties of Oncogenic Nucleosomes.

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 that uses strain fluctuations to determine the Young's modulus of nucleosomes from all-atom molecular dynamics simulations. Recently, we validated this computational tool against in vitro single-molecule nanoindentation of histone variant nucleosomes. In this report, we further extend minimal cylinder analysis 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. Here, we report that the heterotypic nucleosome has an intermediate elasticity (8.5 ± 0.5 MPa) compared to CENP-A (6.2 ± 0.4 MPa) and H3 (9.8 ± 0.7 MPa) and that the dynamics of both canonical and CENP-A nucleosomes are preserved and partitioned across the nucleosome pseudodyad. Furthermore, we investigate the mechanism by which the elasticity of these heterotypic nucleosomes augments cryptic binding surfaces. From these analyses, we predict that the heterotypic nucleosome is permissive to the binding of one copy of the kinetochore protein CENP-C while still retaining a closed DNA end configuration required for linker histone H1 to bind. We discuss that the ectopic deposition of CENP-A in cancer by H3.3 chaperones HIRA and DAXX may fortuitously result in hybrid nucleosome formation. Using these results, we propose biological outcomes that might arise when such heterotypic nucleosomes occupy large regions of the genome.

Intrinsic elasticity of nucleosomes is encoded by histone variants and calibrated by their binding partners.

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.

Turnover versus treadmilling in actin network assembly and remodeling.

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.

Remarkable structural transformations of actin bundles are driven by their initial polarity, motor activity, crosslinking, and filament treadmilling.

Bundled actin structures play a key role in maintaining cellular shape, in aiding force transmission to and from extracellular substrates, and in affecting cellular motility. Recent studies have also brought to light new details on stress generation, force transmission and contractility of actin bundles. In this work, we are primarily interested in the question of what determines the stability of actin bundles and what network geometries do unstable bundles eventually transition to. To address this problem, we used the MEDYAN mechano-chemical force field, modeling several micron-long actin bundles in 3D, while accounting for a comprehensive set of chemical, mechanical and transport processes. We developed a hierarchical clustering algorithm for classification of the different long time scale morphologies in our study. Our main finding is that initially unipolar bundles are significantly more stable compared with an apolar initial configuration. Filaments within the latter bundles slide easily with respect to each other due to myosin activity, producing a loose network that can be subsequently severely distorted. At high myosin concentrations, a morphological transition to aster-like geometries was observed. We also investigated how actin treadmilling rates influence bundle dynamics, and found that enhanced treadmilling leads to network fragmentation and disintegration, while this process is opposed by myosin and crosslinking activities. Interestingly, treadmilling bundles with an initial apolar geometry eventually evolve to a whole gamut of network morphologies based on relative positions of filament ends, such as sarcomere-like organization. We found that apolar bundles show a remarkable sensitivity to environmental conditions, which may be important in enabling rapid cytoskeletal structural reorganization and adaptation in response to intracellular and extracellular cues.

Quantifying dissipation in actomyosin networks.

Quantifying entropy production in various active matter phases will open new avenues for probing self-organization principles in these far-from-equilibrium systems. It has been hypothesized that the dissipation of free energy by active matter systems may be optimized, leading to system trajectories with histories of large dissipation and an accompanying emergence of ordered dynamical states. This interesting idea has not been widely tested. In particular, it is not clear whether emergent states of actomyosin networks, which represent a salient example of biological active matter, self-organize following the principle of dissipation optimization. In order to start addressing this question using detailed computational modelling, we rely on the MEDYAN simulation platform, which allows simulating active matter networks from fundamental molecular principles. We have extended the capabilities of MEDYAN to allow quantification of the rates of dissipation resulting from chemical reactions and relaxation of mechanical stresses during simulation trajectories. This is done by computing precise changes in Gibbs free energy accompanying chemical reactions using a novel formula and through detailed calculations of instantaneous values of the system's mechanical energy. We validate our approach with a mean-field model that estimates the rates of dissipation from filament treadmilling. Applying this methodology to the self-organization of small disordered actomyosin networks, we find that compact and highly cross-linked networks tend to allow more efficient transduction of chemical free energy into mechanical energy. In these simple systems, we observe that spontaneous network reorganizations tend to result in a decrease in the total dissipation rate to a low steady-state value. Future studies might carefully test whether the dissipation-driven adaptation hypothesis applies in this instance, as well as in more complex cytoskeletal geometries.