Expansion of Intrinsically Disordered Proteins Increases the Range of Stability of Liquid-Liquid Phase Separation.

Proteins containing intrinsically disordered regions (IDRs) are ubiquitous within biomolecular condensates, which are liquid-like compartments within cells formed through liquid-liquid phase separation (LLPS). The sequence of amino acids of a protein encodes its phase behaviour, not only by establishing the patterning and chemical nature (e.g., hydrophobic, polar, charged) of the various binding sites that facilitate multivalent interactions, but also by dictating the protein conformational dynamics. Besides behaving as random coils, IDRs can exhibit a wide-range of structural behaviours, including conformational switching, where they transition between alternate conformational ensembles. Using Molecular Dynamics simulations of a minimal coarse-grained model for IDRs, we show that the role of protein conformation has a non-trivial effect in the liquid-liquid phase behaviour of IDRs. When an IDR transitions to a conformational ensemble enriched in disordered extended states, LLPS is enhanced. In contrast, IDRs that switch to ensembles that preferentially sample more compact and structured states show inhibited LLPS. This occurs because extended and disordered protein conformations facilitate LLPS-stabilising multivalent protein-protein interactions by reducing steric hindrance; thereby, such conformations maximize the molecular connectivity of the condensed liquid network. Extended protein configurations promote phase separation regardless of whether LLPS is driven by homotypic and/or heterotypic protein-protein interactions. This study sheds light on the link between the dynamic conformational plasticity of IDRs and their liquid-liquid phase behaviour

Valency and Binding Affinity Variations Can Regulate the Multilayered Organization of Protein Condensates with Many Components.

Biomolecular condensates, which assemble via the process of liquid-liquid phase separation (LLPS), are multicomponent compartments found ubiquitously inside cells. Experiments and simulations have shown that biomolecular condensates with many components can exhibit multilayered organizations. Using a minimal coarse-grained model for interacting multivalent proteins, we investigate the thermodynamic parameters governing the formation of multilayered condensates through changes in protein valency and binding affinity. We focus on multicomponent condensates formed by scaffold proteins (high-valency proteins that can phase separate on their own via homotypic interactions) and clients (proteins recruited to condensates via heterotypic scaffold-client interactions). We demonstrate that higher valency species are sequestered to the center of the multicomponent condensates, while lower valency proteins cluster towards the condensate interface. Such multilayered condensate architecture maximizes the density of LLPS-stabilizing molecular interactions, while simultaneously reducing the surface tension of the condensates. In addition, multilayered condensates exhibit rapid exchanges of low valency proteins in and out, while keeping higher valency proteins-the key biomolecules involved in condensate nucleation-mostly within. We also demonstrate how modulating the binding affinities among the different proteins in a multicomponent condensate can significantly transform its multilayered structure, and even trigger fission of a condensate into multiple droplets with different compositions.Engineering and Physical Sciences Research Council (EPSRC) scholarship to Ignacio Sanchez-Burgo

Size conservation emerges spontaneously in biomolecular condensates formed by scaffolds and surfactant clients

Funder: Junior Research Fellow at Kings CollegeFunder: Ernest Oppenheimer Memorial Trust; doi: http://dx.doi.org/10.13039/501100009978Abstract: Biomolecular condensates are liquid-like membraneless compartments that contribute to the spatiotemporal organization of proteins, RNA, and other biomolecules inside cells. Some membraneless compartments, such as nucleoli, are dispersed as different condensates that do not grow beyond a certain size, or do not present coalescence over time. In this work, using a minimal protein model, we show that phase separation of binary mixtures of scaffolds and low-valency clients that can act as surfactants—i.e., that significantly reduce the droplet surface tension—can yield either a single drop or multiple droplets that conserve their sizes on long timescales (herein ‘multidroplet size-conserved’ scenario’), depending on the scaffold to client ratio. Our simulations demonstrate that protein connectivity and condensate surface tension regulate the balance between these two scenarios. The multidroplet size-conserved scenario spontaneously arises at increasing surfactant-to-scaffold concentrations, when the interfacial penalty for creating small liquid droplets is sufficiently reduced by the surfactant proteins that are preferentially located at the interface. In contrast, low surfactant-to-scaffold concentrations enable continuous growth and fusion of droplets without restrictions. Overall, our work proposes one thermodynamic mechanism to help rationalize how size-conserved coexisting condensates can persist inside cells—shedding light on the roles of protein connectivity, binding affinity, and droplet composition in this process

RNA length has a non-trivial effect in the stability of biomolecular condensates formed by RNA-binding proteins.

Funder: Oppenheimer FellowshipFunder: Roger Ekins FellowshipFunder: Derek Brewer Emmanuel College scholarshipBiomolecular condensates formed via liquid-liquid phase separation (LLPS) play a crucial role in the spatiotemporal organization of the cell material. Nucleic acids can act as critical modulators in the stability of these protein condensates. To unveil the role of RNA length in regulating the stability of RNA binding protein (RBP) condensates, we present a multiscale computational strategy that exploits the advantages of a sequence-dependent coarse-grained representation of proteins and a minimal coarse-grained model wherein proteins are described as patchy colloids. We find that for a constant nucleotide/protein ratio, the protein fused in sarcoma (FUS), which can phase separate on its own-i.e., via homotypic interactions-only exhibits a mild dependency on the RNA strand length. In contrast, the 25-repeat proline-arginine peptide (PR25), which does not undergo LLPS on its own at physiological conditions but instead exhibits complex coacervation with RNA-i.e., via heterotypic interactions-shows a strong dependence on the length of the RNA strands. Our minimal patchy particle simulations suggest that the strikingly different effect of RNA length on homotypic LLPS versus RBP-RNA complex coacervation is general. Phase separation is RNA-length dependent whenever the relative contribution of heterotypic interactions sustaining LLPS is comparable or higher than those stemming from protein homotypic interactions. Taken together, our results contribute to illuminate the intricate physicochemical mechanisms that influence the stability of RBP condensates through RNA inclusion

Chromatin Unfolding by Epigenetic Modifications Explained by Dramatic Impairment of Internucleosome Interactions: A Multiscale Computational Study.

Histone tails and their epigenetic modifications play crucial roles in gene expression regulation by altering the architecture of chromatin. However, the structural mechanisms by which histone tails influence the interconversion between active and inactive chromatin remain unknown. Given the technical challenges in obtaining detailed experimental characterizations of the structure of chromatin, multiscale computations offer a promising alternative to model the effect of histone tails on chromatin folding. Here we combine multimicrosecond atomistic molecular dynamics simulations of dinucleosomes and histone tails in explicit solvent and ions, performed with three different state-of-the-art force fields and validated by experimental NMR measurements, with coarse-grained Monte Carlo simulations of 24-nucleosome arrays to describe the conformational landscape of histone tails, their roles in chromatin compaction, and the impact of lysine acetylation, a widespread epigenetic change, on both. We find that while the wild-type tails are highly flexible and disordered, the dramatic increase of secondary-structure order by lysine acetylation unfolds chromatin by decreasing tail availability for crucial fiber-compacting internucleosome interactions. This molecular level description of the effect of histone tails and their charge modifications on chromatin folding explains the sequence sensitivity and underscores the delicate connection between local and global structural and functional effects. Our approach also opens new avenues for multiscale processes of biomolecular complexes.This work was supported by the European Union Seventh Framework Programme (FP7/2007–2013) [275096 to R.C.-G. and M.O.]; the European Union’s Horizon 2020 research and innovation programme under a Marie Sklodowska-Curie grant [654812 to G.P. and M.V.]; Sara Borrell Fellowships [to G.P. and M.O.]; the Spanish MINECO [BIO2012–32868 to M.O.]; the Spanish National Institute of Bioinformatics (INB) [to M.O.]; the European Research Council (ERC) [Advanced Investigator Grant to M.O.]; the National Science Foundation [MCB0316771 to T. S.]; the National Institutes of Health [R01 GM55164 to T. S.]; Philip Morris USA [to T. S.]; and Philip Morris International [to T.S.]. M.O. is an ICREA-Academia fellow

Deoxyribonucleic Acid Encoded and Size-Defined π-Stacking of Perylene Diimides.

Funder: University of CambridgeNatural photosystems use protein scaffolds to control intermolecular interactions that enable exciton flow, charge generation, and long-range charge separation. In contrast, there is limited structural control in current organic electronic devices such as OLEDs and solar cells. We report here the DNA-encoded assembly of π-conjugated perylene diimides (PDIs) with deterministic control over the number of electronically coupled molecules. The PDIs are integrated within DNA chains using phosphoramidite coupling chemistry, allowing selection of the DNA sequence to either side, and specification of intermolecular DNA hybridization. In this way, we have developed a "toolbox" for construction of any stacking sequence of these semiconducting molecules. We have discovered that we need to use a full hierarchy of interactions: DNA guides the semiconductors into specified close proximity, hydrophobic-hydrophilic differentiation drives aggregation of the semiconductor moieties, and local geometry and electrostatic interactions define intermolecular positioning. As a result, the PDIs pack to give substantial intermolecular π wave function overlap, leading to an evolution of singlet excited states from localized excitons in the PDI monomer to excimers with wave functions delocalized over all five PDIs in the pentamer. This is accompanied by a change in the dominant triplet forming mechanism from localized spin-orbit charge transfer mediated intersystem crossing for the monomer toward a delocalized excimer process for the pentamer. Our modular DNA-based assembly reveals real opportunities for the rapid development of bespoke semiconductor architectures with molecule-by-molecule precision.ERC Horizon 2020 (grant agreement No 670405 and No 803326) EPSRC Tier-2 capital grant EP/P020259/1. Winton Advanced Research Programme for the Physics of Sustainability. Simons Foundation (Grant 601946). Swedish research council, Vetenskapsrådet 2018-0023

Integrating transposable elements in the 3D genome

Chromosome organisation is increasingly recognised as an essential component of genome regulation, cell fate and cell health. Within the realm of transposable elements (TEs) however, the spatial information of how genomes are folded is still only rarely integrated in experimental studies or accounted for in modelling. Whilst polymer physics is recognised as an important tool to understand the mechanisms of genome folding, in this commentary we discuss its potential applicability to aspects of TE biology. Based on recent works on the relationship between genome organisation and TE integration, we argue that existing polymer models may be extended to create a predictive framework for the study of TE integration patterns. We suggest that these models may offer orthogonal and generic insights into the integration profiles (or "topography") of TEs across organisms. In addition, we provide simple polymer physics arguments and preliminary molecular dynamics simulations of TEs inserting into heterogeneously flexible polymers. By considering this simple model, we show how polymer folding and local flexibility may generically affect TE integration patterns. The preliminary discussion reported in this commentary is aimed to lay the foundations for a large-scale analysis of TE integration dynamics and topography as a function of the three-dimensional host genome

Bimolecular reaction rates from ring polymer molecular dynamics.

We describe an efficient procedure for calculating the rates of bimolecular chemical reactions in the gas phase within the ring polymer molecular dynamics approximation. A key feature of the procedure is that it does not require that one calculate the absolute quantum mechanical partition function of the reactants or the transition state: The rate coefficient only depends on the ratio of these two partition functions which can be obtained from a thermodynamic integration along a suitable reaction coordinate. The procedure is illustrated with applications to the three-dimensional H + H(2), Cl + HCl, and F + H(2) reactions, for which well-converged quantum reactive scattering results are computed for comparison. The ring polymer rate coefficients agree with these exact results at high temperatures and are within a factor of 3 of the exact results at temperatures in the deep quantum tunneling regime, where the classical rate coefficients are too small by several orders of magnitude. This is probably already good enough to encourage future applications of the ring polymer theory to more complex chemical reactions, which it is capable of treating in their full dimensionality. However, there is clearly some scope for improving on the ring polymer approximation at low temperatures, and we end by suggesting a way in which this might be accomplished

Bimolecular reaction rates from ring polymer molecular dynamics: application to H + CH4 → H2 + CH3.

In a recent paper, we have developed an efficient implementation of the ring polymer molecular dynamics (RPMD) method for calculating bimolecular chemical reaction rates in the gas phase, and illustrated it with applications to some benchmark atom-diatom reactions. In this paper, we show that the same methodology can readily be used to treat more complex polyatomic reactions in their full dimensionality, such as the hydrogen abstraction reaction from methane, H + CH(4) → H(2) + CH(3). The present calculations were carried out using a modified and recalibrated version of the Jordan-Gilbert potential energy surface. The thermal rate coefficients obtained between 200 and 2000 K are presented and compared with previous results for the same potential energy surface. Throughout the temperature range that is available for comparison, the RPMD approximation gives better agreement with accurate quantum mechanical (multiconfigurational time-dependent Hartree) calculations than do either the centroid density version of quantum transition state theory (QTST) or the quantum instanton (QI) model. The RPMD rate coefficients are within a factor of 2 of the exact quantum mechanical rate coefficients at temperatures in the deep tunneling regime. These results indicate that our previous assessment of the accuracy of the RPMD approximation for atom-diatom reactions remains valid for more complex polyatomic reactions. They also suggest that the sensitivity of the QTST and QI rate coefficients to the choice of the transition state dividing surface becomes more of an issue as the dimensionality of the reaction increases

Bimolecular reaction rates from ring polymer molecular dynamics: application to H + CH4 → H2 + CH3.

In a recent paper, we have developed an efficient implementation of the ring polymer molecular dynamics (RPMD) method for calculating bimolecular chemical reaction rates in the gas phase, and illustrated it with applications to some benchmark atom-diatom reactions. In this paper, we show that the same methodology can readily be used to treat more complex polyatomic reactions in their full dimensionality, such as the hydrogen abstraction reaction from methane, H + CH(4) → H(2) + CH(3). The present calculations were carried out using a modified and recalibrated version of the Jordan-Gilbert potential energy surface. The thermal rate coefficients obtained between 200 and 2000 K are presented and compared with previous results for the same potential energy surface. Throughout the temperature range that is available for comparison, the RPMD approximation gives better agreement with accurate quantum mechanical (multiconfigurational time-dependent Hartree) calculations than do either the centroid density version of quantum transition state theory (QTST) or the quantum instanton (QI) model. The RPMD rate coefficients are within a factor of 2 of the exact quantum mechanical rate coefficients at temperatures in the deep tunneling regime. These results indicate that our previous assessment of the accuracy of the RPMD approximation for atom-diatom reactions remains valid for more complex polyatomic reactions. They also suggest that the sensitivity of the QTST and QI rate coefficients to the choice of the transition state dividing surface becomes more of an issue as the dimensionality of the reaction increases