Chromatin dynamics controls epigenetic domain formation.

In multicellular organisms, nucleosomes carry epigenetic information that defines distinct patterns of gene expression, which are inherited over multiple generations. The enhanced capacity for information storage arises by nucleosome modifications, which are triggered by enzymes. Modified nucleosomes can transfer the mark to others that are in proximity by a positive-feedback (modification begets modification) mechanism. We created a generic polymer model, referred to as 3DSpreader, in which each bead, representing a nucleosome, stochastically switches between unmodified (U) and modified (M) states depending on the states of the neighbors. Modification begins at a specific nucleation site (NS) that is permanently in the M state, and could spread to other loci that is dictated by chromatin dynamics. Transfer of marks among the non-nucleation loci occurs stochastically as chromatin evolves in time. If the spreading rate is slower than the chromatin relaxation rate, which is biologically pertinent, then finite-sized domains form, driven by contacts between nucleosomes through a three-dimensional looping mechanism. Surprisingly, simulations based on the 3DSpreader model result in finite bounded domains that arise without the need for any boundary elements. Maintenance of spatially and temporally stable domains requires the presence of the NS, whose removal eliminates finite-sized modified domains. The theoretical predictions are in excellent agreement with experimental data for H3K9me3 spreading in mouse embryonic stem cells.

Critical structural fluctuations of proteins upon thermal unfolding challenge the Lindemann criterion.

Internal subnanosecond timescale motions are key for the function of proteins, and are coupled to the surrounding solvent environment. These fast fluctuations guide protein conformational changes, yet their role for protein stability, and for unfolding, remains elusive. Here, in analogy with the Lindemann criterion for the melting of solids, we demonstrate a common scaling of structural fluctuations of lysozyme protein embedded in different environments as the thermal unfolding transition is approached. By combining elastic incoherent neutron scattering and advanced molecular simulations, we show that, although different solvents modify the protein melting temperature, a unique dynamical regime is attained in proximity of thermal unfolding in all solvents that we tested. This solvation shell-independent dynamical regime arises from an equivalent sampling of the energy landscape at the respective melting temperatures. Thus, we propose that a threshold for the conformational entropy provided by structural fluctuations of proteins exists, beyond which thermal unfolding is triggered.

Capillary assembly of microscale ellipsoidal, cuboidal, and spherical particles at interfaces.

Micron-sized anisotropic particles with homogeneous surface properties at a fluid interface can deform the interface due to their shape. The particles thereby create excess interfacial area and interact in order to minimize this area, which lowers the total interfacial energy. We present a systematic investigation of the interface deformations around single ellipsoidal particles and cuboidal particles with rounded edges in the near field for various contact angles and particle aspect ratios. The correlation of these deformations with capillary bond energies-the interaction energies of two particles at contact-quantifies the relation between the interactions and the near-field deformations. We characterize the interactions using effective power laws and investigate how anisotropic particles self-assemble by capillary forces. Interface deformations and particle interactions for cuboidal particles are weaker compared with those for ellipsoidal particles with the same aspect ratios. For both particle shapes, the bound state in side-by-side orientation is most stable, while the interaction in tip-to-side orientation is repulsive. Furthermore, we find capillary attraction between spherical and ellipsoidal particles. Our calculations therefore suggest cluster formation of spherical and ellipsoidal particles, which elucidates the role of spherical particles as stoppers for the growth of worm-like chains of ellipsoidal particles. The interaction between spherical and ellipsoidal particles might also explain the suppression of the "coffee-ring effect" that has been observed for evaporating droplets with mixtures of spherical and ellipsoidal particles. In general, our calculations of the near-field interactions complement previous calculations in the far field and help to predict colloidal assembly and rheological properties of particle-laden interfaces.

Specific Interactions and Environment Flexibility Tune Protein Stability under Extreme Crowding.

Macromolecular crowding influences protein mobility and stability in vivo. A precise description of the crowding effect on protein thermal stability requires the estimate of the combined effects of excluded volume, specific protein-environment interactions, as well as the thermal response of the crowders. Here, we explore an ideal model system, the lysozyme protein in powder state, to dissect the factors controlling the melting of the protein under extreme crowding. By deploying state-of-the art molecular simulations, supported by calorimetric experiments, we assess the role of the environment flexibility and of intermolecular electrostatic interactions. In particular, we show that the temperature-dependent flexibility of the macromolecular crowders, along with specific interactions, significantly alleviates the stabilizing contributions of the static volume effect.

Configurational Disorder of Water Hydrogen-Bond Network at the Protein Dynamical Transition.

We introduce a novel strategy to quantify the disorder of extended water-water hydrogen-bond (HB) networks sampled in particle-based computer simulations. The method relies on the conformational clustering of the HB connectivity states. We successfully applied it to unveil the fine relationship among the protein dynamical transition in hydrated powder, which marks the activation of protein flexibility at Td ≈ 240 K, and the sudden increase in the configurational disorder of the water HB network enveloping the proteins. Our finding links, in the spirit of the Adam-Gibbs relationship, the diffusivity of protein atoms, as quantified by the hydrogen mean-square displacements, and the thermodynamic solvent configurational entropy.

Thermal activation of 'allosteric-like' large-scale motions in a eukaryotic Lactate Dehydrogenase.

Conformational changes occurring during the enzymatic turnover are essential for the regulation of protein functionality. Individuating the protein regions involved in these changes and the associated mechanical modes is still a challenge at both experimental and theoretical levels. We present here a detailed investigation of the thermal activation of the functional modes and conformational changes in a eukaryotic Lactate Dehydrogenase enzyme (LDH). Neutron Spin Echo spectroscopy and Molecular Dynamics simulations were used to uncover the characteristic length- and timescales of the LDH nanoscale motions in the apo state. The modes involving the catalytic loop and the mobile region around the binding site are activated at room temperature, and match the allosteric reorganisation of bacterial LDHs. In a temperature window of about 15 degrees, these modes render the protein flexible enough and capable of reorganising the active site toward reactive configurations. On the other hand an excess of thermal excitation leads to the distortion of the protein matrix with a possible anti-catalytic effect. Thus, the temperature activates eukaryotic LDHs via the same conformational changes observed in the allosteric bacterial LDHs. Our investigation provides an extended molecular picture of eukaryotic LDH's conformational landscape that enriches the static view based on crystallographic studies alone.

Stability and Function at High Temperature. What Makes a Thermophilic GTPase Different from Its Mesophilic Homologue.

Comparing homologous enzymes adapted to different thermal environments aids to shed light on their delicate stability/function trade-off. Protein mechanical rigidity was postulated to secure stability and high-temperature functionality of thermophilic proteins. In this work, we challenge the corresponding-state principle for a pair of homologous GTPase domains by performing extensive molecular dynamics simulations, applying conformational and kinetic clustering, as well as exploiting an enhanced sampling technique (REST2). While it was formerly shown that enhanced protein flexibility and high temperature stability can coexist in the apo hyperthermophilic variant, here we focus on the holo states of both homologues by mimicking the enzymatic turnover. We clearly show that the presence of the ligands affects the conformational landscape visited by the proteins, and that the corresponding state principle applies for some functional modes. Namely, in the hyperthermophilic species, the flexibility of the effector region ensuring long-range communication and of the P-loop modulating ligand binding are recovered only at high temperature.