Search

Colloidal systems possess unique features to investigate the governing principles behind liquid-to-solid transitions. The phase diagram and crystallization landscape of colloidal particles can be finely tuned by the range, number and angular distribution of attractive interactions between the constituent particles. In this work, we present a computational study of colloidal patchy particles with high-symmetry bonding—six patches displaying octahedral symmetry—that can crystallize into distinct competing ordered phases: a cubic simple (CS) lattice, a body-centered cubic (BCC) phase, and two face-centered cubic (FCC) solids (orientationally ordered and disordered). We investigate the underlying mechanisms by which these competing crystals emerge from a disordered fluid at different pressures. Strikingly, we identify instances where the structure of the crystalline embryo corresponds to the stable solid, while in others it corresponds to a metastable crystal whose nucleation is enabled by its lower interfacial free energy with the liquid. Moreover, we find the exceptional phenomenon that, due to a subtle balance between volumetric enthalpy and interfacial free energy, the CS phase nucleates via crystalline cubic nuclei rather than through spherical clusters as the majority of crystal solids in nature. Finally, by examining growth beyond the nucleation stage, we uncover a series of alternating one-phase and two-phase crystallization mechanisms, depending on whether or not the same phase that nucleates keeps growing. Taken together, we show that an octahedral distribution of attractive sites in colloidal particles results in an extremely rich crystallization landscape where subtle differences in pressure crucially determine the crystallizing polymorph.

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. We find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

Various physics- and data-driven sequence-dependent protein coarse-grained models have been developed to study biomolecular phase separation and elucidate the dominant physicochemical driving forces. Here we present Mpipi, a multiscale coarse-grained model that describes almost quantitatively the change in protein critical temperatures as a function of amino acid sequence. The model is parameterized from both atomistic simulations and bioinformatics data and accounts for the dominant role of π–π and hybrid cation–π/π–π interactions and the much stronger attractive contacts established by arginines than lysines. We provide a comprehensive set of benchmarks for Mpipi and seven other residue-level coarse-grained models against experimental radii of gyration and quantitative in vitro phase diagrams, demonstrating that Mpipi predictions agree well with experiments on both fronts. Moreover, Mpipi can account for protein–RNA interactions, correctly predicts the multiphase behavior of a charge-matched poly-arginine/poly-lysine/RNA system, and recapitulates experimental liquid–liquid phase separation trends for sequence mutations on FUS, DDX4 and LAF-1 proteins. Combining bioinformatics data and atomistic simulations, this study develops a sequence-dependent coarse-grained model for biomolecular phase separation. This model achieves a quantitative agreement with experimental observations. Extensive benchmarks exemplify its performance.

Despite the importance of ice nucleation, this process has been barely explored at negative pressures. Here, we study homogeneous ice nucleation in stretched water by means of molecular dynamics seeding simulations using the TIP4P/Ice model. We observe that the critical nucleus size, interfacial free energy, free energy barrier, and nucleation rate barely change between isobars from −2600 to 500 bars when they are represented as a function of supercooling. This allows us to identify universal empirical expressions for homogeneous ice nucleation in the pressure range from −2600 to 500 bars. We show that this universal behavior arises from the pressure dependence of the interfacial free energy, which we compute by means of the mold integration technique, finding a shallow minimum around −2000 bars. Likewise, we show that the change in the interfacial free energy with pressure is proportional to the excess entropy and the slope of the melting line, exhibiting in the latter a reentrant behavior also at the same negative pressure. Finally, we estimate the excess internal energy and the excess entropy of the ice Ih–water interface.

Biomolecular condensates are important contributors to the internal organization of the cell material. While initially described as liquid-like droplets, the term biomolecular condensates is now used to describe a diversity of condensed phase assemblies with material properties extending from low to high viscous liquids, gels, and even glasses. Because the material properties of condensates are determined by the intrinsic behaviour of their molecules, characterising such properties is integral to rationalising the molecular mechanisms that dictate their functions and roles in health and disease. Here, we apply and compare three distinct computational methods to measure the viscoelasticity of biomolecular condensates in molecular simulations. These methods are the shear stress relaxation modulus integration (SSRMI), the oscillatory shear (OS) technique, and the bead tracking (BT) method. We find that, although all of these methods provide consistent results for the viscosity of the condensates, the SSRMI and OS techniques outperform the BT method in terms of computational efficiency and statistical uncertainty. We, thus, apply the SSRMI and OS techniques for a set of 12 different protein/RNA systems using a sequence-dependent high-resolution coarse-grained model. Our results reveal a strong correlation between condensate viscosity and density, as well as with protein/RNA length and the number of stickersvs.spacers in the amino-acid protein sequence. Moreover, we couple the SSRMI and the OS technique to nonequilibrium molecular dynamics simulations that mimic the progressive liquid-to-gel transition of protein condensates due to the accumulation of inter-proteinβ-sheets. We compare the behaviour of three different protein condensates—i.e., those formed by either hnRNPA1, FUS, or TDP-43 proteins—whose liquid-to-gel transitions are associated with the onset of amyotrophic lateral sclerosis and frontotemporal dementia. We find that both SSRMI and OS techniques successfully predict the transition from functional liquid-like behaviour to kinetically arrested states once the network of inter-proteinβ-sheets has percolated through the condensates. Overall, our work provides a comparison of different modelling rheological techniques to assess the viscosity of biomolecular condensates, a critical magnitude that provides information on the behaviour of biomolecules inside condensates.

Funder: European Research Council (ERC), The protein high mobility group A1 (HMGA1) is an important regulator of chromatin organization and function. However, the mechanisms by which it exerts its biological function are not fully understood. Here, we report that the HMGA isoform, HMGA1a, nucleates into foci that display liquid-like properties in the nucleus, and that the protein readily undergoes phase separation to form liquid condensates in vitro. By bringing together machine-leaning modelling, cellular and biophysical experiments and multiscale simulations, we demonstrate that phase separation of HMGA1a is promoted by protein-DNA interactions, and has the potential to be modulated by post-transcriptional effects such as phosphorylation. We further show that the intrinsically disordered C-terminal tail of HMGA1a significantly contributes to its phase separation through electrostatic interactions via AT hooks 2 and 3. Our work sheds light on HMGA1 phase separation as an emergent biophysical factor in regulating chromatin structure.

Costimulation blockade (CoB)–based immunosuppression offers the promise of improved transplantation outcomes with reduced drug toxicity. However, it is hampered by early acute rejections, mediated at least in part by differentiated, CoB-resistant T cells, such as CD57+PD1− CD4 T cells. In this study, we characterize these cells pretransplant, determine their fate posttransplant, and examine their proliferative capacity in vitro in humans. Our studies show that CD57+PD1− CD4 T cells are correlated with increasing age and CMV infection pretransplant, and persist for up to 1 y posttransplant. These cells are replication incompetent alone but proliferated in the presence of unsorted PBMCs in a contact-independent manner. When stimulated, cells sorted by CD57/PD1 status upregulate markers of activation with proliferation. Up to 85% of CD57+PD1− cells change expression of CD57/PD1 with stimulation, typically, upregulating PD1 and downregulating CD57. PD1 upregulation is accentuated in the presence of rapamycin but prevented by tacrolimus. These data support a general theory of CoB-resistant cells as Ag-experienced, costimulation-independent cells and suggest a mechanism for the synergy of belatacept and rapamycin, with increased expression of the activation marker PD1 potentiating exhaustion of CoB-resistant cells.

Funder: Ernest Oppenheimer Memorial Trust, Biomolecular condensates, some of which are liquid-like during health, can age over time becoming gel-like pathological systems. One potential source of loss of liquid-like properties during ageing of RNA-binding protein condensates is the progressive formation of inter-protein β-sheets. To bridge microscopic understanding between accumulation of inter-protein β-sheets over time and the modulation of FUS and hnRNPA1 condensate viscoelasticity, we develop a multiscale simulation approach. Our method integrates atomistic simulations with sequence-dependent coarse-grained modelling of condensates that exhibit accumulation of inter-protein β-sheets over time. We reveal that inter-protein β-sheets notably increase condensate viscosity but does not transform the phase diagrams. Strikingly, the network of molecular connections within condensates is drastically altered, culminating in gelation when the network of strong β-sheets fully percolates. However, high concentrations of RNA decelerate the emergence of inter-protein β-sheets. Our study uncovers molecular and kinetic factors explaining how the accumulation of inter-protein β-sheets can trigger liquid-to-solid transitions in condensates, and suggests a potential mechanism to slow such transitions down.

Freezing of water is the most common liquid-to-crystal phase transition on Earth; however, despite its critical implications on climate change and cryopreservation among other disciplines, its characterization through experimental and computational techniques remains elusive. In this work, we make use of computer simulations to measure the nucleation rate ( J) of water at normal pressure under different supercooling conditions, ranging from 215 to 240 K. We employ two different water models: mW, a coarse-grained potential for water, and TIP4P/ICE, an atomistic nonpolarizable water model that provides one of the most accurate representations of the different ice phases. To evaluate J, we apply the Lattice Mold technique, a computational method based on the use of molds to induce the nucleus formation from the metastable liquid under conditions at which observing spontaneous nucleation would be unfeasible. With this method, we obtain estimates of the nucleation rate for ice Ih and Ic and a stacking mixture of ice Ih/Ic, reaching consensus with most of the previously reported rates, although differing with some others. Furthermore, we confirm that the predicted nucleation rates obtained by the TIP4P/ICE model are in better agreement with experimental data than those obtained through the mW potential. Taken together, our study provides a reliable methodology to measure nucleation rates in a simple and computationally efficient manner that contributes to benchmarking the freezing behavior of two popular water models.

Hard-sphere crystallization has been widely investigated over the last six decades by means of colloidal suspensions and numerical methods. However, some aspects of its nucleation behaviour are still under debate. Here, we provide a detailed computational characterisation of the polymorphic nucleation competition between the face-centered cubic (fcc) and the hexagonal-close packed (hcp) hard-sphere crystal phases. By means of several state-of-the-art simulation techniques, we evaluate the melting pressure, chemical potential difference, interfacial free energy and nucleation rate of these two polymorphs, as well as of a random stacking mixture of both crystals. Our results highlight that, despite the fact that both polymorphs have very similar stability, the interfacial free energy of the hcp phase could be marginally higher than that of the fcc solid, which in consequence, mildly decreases its propensity to nucleate from the liquid compared to the fcc phase. Moreover, we analyse the abundance of each polymorph in grown crystals from different types of inserted nuclei: fcc, hcp and stacking disordered fcc/hcp seeds, as well as from those spontaneously emerged from brute force simulations. We find that post-critical crystals fundamentally grow maintaining the polymorphic structure of the critical nucleus, at least until moderately large sizes, since the only crystallographic orientation that allows stacking close-packed disorder is the fcc (111) plane, or equivalently the hcp (0001) one. Taken together, our results contribute with one more piece to the intricate puzzle of colloidal hard-sphere crystallization.

10 pags., 7 figs., 2 tabs., The seeding method is an approximate approach to investigate nucleation that combines molecular dynamics simulations with classical nucleation theory. Recently, this technique has been successfully implemented in a broad range of nucleation studies. However, its accuracy is subject to the arbitrary choice of the order parameter threshold used to distinguish liquid-like from solid-like molecules. We revisit here the crystallization of NaCl from a supersaturated brine solution and show that consistency between seeding and rigorous methods, like Forward Flux Sampling (from previous work) or spontaneous crystallization (from this work), is achieved by following a mislabelling criterion to select such threshold (i.e. equaling the fraction of the mislabelled particles in the bulk parent and nucleating phases). This work supports the use of seeding to obtain fast and reasonably accurate nucleation rate estimates and the mislabelling criterion as one giving the relevant cluster size for classical nucleation theory in crystallization studies., This project has been funded by grants FIS2016-78117-P and PID2019-105898GB-C21 of MEC. E. G. N. thanks Agencia Estatal de Investigacion and Fondo Europeo de Desarrollo Regional (FEDER), Grant No FIS2017-89361-C3-2-P. C. P. L. thanks Ministerio de Educacion y Formacio´n Profesional for a predoctoral Formacion Profesorado Universitario Grant No. FPU18/03326 and also Ayuntamiento de Madrid for a Residencia de Estudiantes grant. The authors acknowledge the computer resources and technical assistance provided by RES. PMdH acknowledges financial support from the FPI grant no. BES6712017-080074. J. R. E. acknowledges funding from the Oppenheimer Research fellowship and the Roger Ekins Research Fellowship of Emmanuel College. M. M. C. thanks financial support from PID2019-105898GA-C22 of the MICINN and CAM and UPM through the CavitieS project No. APOYOJOVENES-01HQ1S-129-B5E4MM from ‘‘Accion financiada por la Comunidad de Madrid en el marco del Convenio Plurianual con la Universidad Politecnica de Madrid en la linea de actuacion estimulo a la investigacion de jovenes doctores’’. The authors gratefully acknowledge Universidad Politecnica de Madrid (www.upm.es) for providing computing resources on Magerit Supercomputer.

A biosolid was produced from the composted digestate of the organic fraction of urban solid wastes (OFUSW). During a first hydrolytic-acidogenic stage, the milled OFUSW mass is transformed into volatile fatty acids and other soluble or suspended substances in the leachate while the residual solids constitute the digestate. A biosolid was produced from the composted digestate of the organic fraction of urban solid wastes (OFUSW). During a first hydrolytic-acidogenic stage, the milled OFUSW mass is transformed into volatile fatty acids and other soluble or suspended substances in the leachate while the residual solids constitute the digestate. To obtain a composting mass (CM) with a C/N between 25 and 35 and a matrix to allow aeration, mixtures of digestate (fD), OFUSW (fF) and sawdust (fS) were prepared using a simplex centroid design. Results were adjusted to a multiple regression model with volatile solids degradation efficiency (ηVS) as a response variable. It was found that the operation zone with the highest ηVS, were CM with the following compositions: fD: 0.425-0.625, fF: 0.275-0.450; and fS: 0.1-0.15. Three assays were carried out in a bench scale reactor, the final biosolids had a pH = 7.5, a C/N = 15 and a germination index (GI) = 84%. A linear correlation between ηvs and CO2 production was proposed for process control. A rapid CM stabilization was reached, up to ηVS = 35% in 12 days and the biosolids with high GI, can be used to improve soils.

Crystallographic data comes from a space-time average over all the unit cells within the crystal, so dynamic phenomena do not contribute significantly to the diffraction data. Many efforts have been made to reconstitute the movement of the macromolecules and explore the microstates that the confined proteins can adopt in the crystalline network. We explored different strategies to simulate a heart fatty acid binding protein (H-FABP) crystal by means of Molecular Dynamics (MD) simulations. We evaluate the effect of introducing restraints according to experimental isotropic B-factors and we analyzed the H-FABP motions in the crystal using Principal Component Analysis (PCA), isotropic and anisotropic B-factors. We compared the behavior of the protein simulated in the crystal confinement versus in solution, and we observed the effect of that confinement in the mobility of the protein residues. Restraining one-third of Cα atoms based on experimental B-factors produce lower B-factors than simulations without restraints, showing that the position restraint of the atoms with the lowest experimental B-factor is a good strategy to maintain the geometry of the crystal with an obvious decrease in the degrees of motion of the protein. PCA shows that, as position restraint reduces the conformational space explored by the system, the motion of the crystal is better recovered, for an essential subspace of the same size, in the simulations without restraints. Restraining only one Cα seems to be a good balance between giving flexibility to the system and preserving its structure. Fil: Espinosa Silva, Yanis Ricardo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Física de Líquidos y Sistemas Biológicos. Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Instituto de Física de Líquidos y Sistemas Biológicos; Argentina. Universidad Industrial Santander; Colombia Fil: Alvarez, Hugo Ariel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Física de Líquidos y Sistemas Biológicos. Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Instituto de Física de Líquidos y Sistemas Biológicos; Argentina. Universidad Nacional Arturo Jauretche; Argentina Fil: Howard, Eduardo Ignacio. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Física de Líquidos y Sistemas Biológicos. Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Instituto de Física de Líquidos y Sistemas Biológicos; Argentina. Universidad Tecnológica Nacional; Argentina Fil: Carlevaro, Carlos Manuel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Física de Líquidos y Sistemas Biológicos. Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Instituto de Física de Líquidos y Sistemas Biológicos; Argentina. Universidad Tecnológica Nacional; Argentina

The synthesis of triarylmethanes via Pd-catalyzed Suzuki-Miyaura reactions between diarylmethyl 2,3,4,5,6-pentafluorobenzoates and aryl boronic acids is described. The system operates at mild conditions and has a broad substrate scope, including the coupling of diphenylmethanol derivatives that do not contain extended aromatic substituents. This is significant as these substrates, which result in the types of triarylmethane products that are prevalent in pharmaceuticals, have not previously been compatible with systems for diarylmethyl ester coupling. Further, the reaction can be performed stereospecifically to generate stereo-inverted products. On the basis of DFT calculations, it is proposed that the oxidative addition of the diarylmethyl 2,3,4,5,6-pentafluorobenzoate substrate occurs via an S(N)2 pathway, which results in the inverted products. Mechanistic studies indicate that oxidative addition of the diarylmethyl 2,3,4,5,6-pentafluorobenzoate substrates to (IPr)Pd(0) results in the selective cleavage of the O–C(benzyl) bond in part because of a stabilizing η(3)-interaction between the benzyl ligand and Pd. This is in contrast to previously described Pd-catalyzed Suzuki-Miyaura reactions involving phenyl esters, which involve selective cleavage of the C(acyl)–O bond, because there is no stabilizing η(3)-interaction. It is anticipated that this fundamental knowledge will aid the development of new catalytic systems, which use esters as electrophiles in cross-coupling reactions.

Biomolecular condensates formed by the process of liquid–liquid phase separation (LLPS) play diverse roles inside cells, from spatiotemporal compartmentalisation to speeding up chemical reactions. Upon maturation, the liquid-like properties of condensates, which underpin their functions, are gradually lost, eventually giving rise to solid-like states with potential pathological implications. Enhancement of inter-protein interactions is one of the main mechanisms suggested to trigger the formation of solid-like condensates. To gain a molecular-level understanding of how the accumulation of stronger interactions among proteins inside condensates affect the kinetic and thermodynamic properties of biomolecular condensates, and their shapes over time, we develop a tailored coarse-grained model of proteins that transition from establishing weak to stronger inter-protein interactions inside condensates. Our simulations reveal that the fast accumulation of strongly binding proteins during the nucleation and growth stages of condensate formation results in aspherical solid-like condensates. In contrast, when strong inter-protein interactions appear only after the equilibrium condensate has been formed, or when they accumulate slowly over time, with respect to the time needed for droplets to fuse and grow, spherical solid-like droplets emerge. By conducting atomistic potential-of-mean-force simulations of NUP-98 peptides—prone to forming inter-protein β -sheets—we observe that formation of inter-peptide β -sheets increases the strength of the interactions consistently with the loss of liquid-like condensate properties we observe at the coarse-grained level. Overall, our work aids in elucidating fundamental molecular, kinetic, and thermodynamic mechanisms linking the rate of change in protein interaction strength to condensate shape and maturation during ageing., Adiran Garaizar is funded by the EPRSC Doctoral Programme Training number EP/N509620/.

The insertion of carbon dioxide into metal element σ-bonds is an important elementary step in many catalytic reactions for carbon dioxide valorization. Here, the insertion of carbon dioxide into a family of group 10 alkyl complexes of the type (RPBP)M(CH3) (RPBP = B(NCH2PR2)2C6H4−; R = Cy or tBu; M = Ni or Pd) to generate κ1-acetate complexes of the form (RPBP)M{OC(O)CH3} is investigated. This involved the preparation and characterization of a number of new complexes supported by the unusual RPBP ligand, which features a central boryl donor that exerts a strong trans-influence, and the identification of a new decomposition pathway that results in C–B bond formation. In contrast to other group 10 methyl complexes supported by pincer ligands, carbon dioxide insertion into (RPBP)M(CH3) is facile and occurs at room temperature because of the high trans-influence of the boryl donor. Given the mild conditions for carbon dioxide insertion, we perform a rare kinetic study on carbon dioxide insertion into a late-transition metal alkyl species using (tBuPBP)Pd(CH3). These studies demonstrate that the Dimroth–Reichardt parameter for a solvent correlates with the rate of carbon dioxide insertion and that Lewis acids do not promote insertion. DFT calculations indicate that insertion into (tBuPBP)M(CH3) (M = Ni or Pd) proceeds via an SE2 mechanism and we compare the reaction pathway for carbon dioxide insertion into group 10 methyl complexes with insertion into group 10 hydrides. Overall, this work provides fundamental insight that will be valuable for the development of improved and new catalysts for carbon dioxide utilization.

Funder: Oppenheimer Fellowship, Funder: Roger Ekins Fellowship, Funder: Derek Brewer Emmanuel College scholarship, Biomolecular 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.

Liquid–liquid phase separation underlies the formation of biological condensates. Physically, such systems are microemulsions which have a general propensity to fuse and coalesce; however, many condensates persist as independent droplets inside cells. This stability is crucial for their functioning, but the physicochemical mechanisms that control the emulsion stability of condensates remain poorly understood. Here, by combining single-condensate zeta potential measurements, optical microscopy, tweezer experiments, and multiscale molecular modelling, we investigate how the forces that sustain condensates impact their stability against fusion. By comparing PR25:PolyU and FUS condensates, we show that a higher condensate surface charge correlates with a lower fusion propensity, and that this behavior can be inferred from their zeta potentials. We reveal that overall stabilization against fusion stems from a combination of repulsive forces between condensates and the effects that surface electrostatics have on lowering surface tension, thus shedding light on the molecular determinants of condensate coalescence.

The insertion of carbon dioxide into metal element σ-bonds is an important elementary step in many catalytic reactions for carbon dioxide valorization. Here, the insertion of carbon dioxide into a family of group 10 alkyl complexes of the type (

The protein high mobility group A1 (HMGA1) is an important regulator of chromatin organization and function. However, the mechanisms by which it exerts its biological function are not fully understood. Here, we report that the HMGA isoform, HMGA1a, nucleates into foci that display liquid-like properties in the nucleus, and that the protein readily undergoes phase separation to form liquid condensates in vitro. By bringing together machine-leaning modelling, cellular and biophysical experiments and multiscale simulations, we demonstrate that phase separation of HMGA1a is critically promoted by protein–DNA interactions, and has the potential to be modulated by post-transcriptional effects such as phosphorylation. We further show that the intrinsically disordered C-terminal tail of HMGA1a significantly contributes to its phase separation through cation–π and electrostatic interactions. Our work sheds light on HMGA1 phase separation as an emergent biophysical factor in regulating chromatin structure.

Phase-separated biomolecular condensates that contain multiple coexisting phases are widespread in vitro and in cells. Multiphase condensates emerge readily within multi-component mixtures of biomolecules (e.g. proteins and nucleic acids) when the different components present sufficient physicochemical diversity (e.g. in inter-molecular forces, structure, and chemical composition) to sustain separate coexisting phases. Because such diversity is highly coupled to the solution conditions (e.g. temperature, pH, salt, composition), it can manifest itself immediately from the nucleation and growth stages of condensate formation, develop spontaneously due to external stimuli, or progressively as the condensates age. Here, we investigate thermodynamic factors that can explain the intrinsic transformation of single-component condensates into multiphase architectures during the nonequilibrium process of ageing. We develop a multiscale model that integrates atomistic simulations of proteins, sequence-dependent coarse-grained simulations of condensates, and a minimal model of dynamically ageing condensates with non-conservative inter-molecular forces. Our nonequilibrium simulations of condensate ageing predict that single-component condensates that are initially homogeneous and liquid-like can transform into gel-core/liquid-shell or liquid-core/gel-shell multiphase condensates as they age, due to gradual and irreversible enhancement of inter-protein interactions. The type of multiphase architecture is determined by the ageing mechanism, the molecular organization of the gel and liquid phases, and the chemical make up of the protein. Notably, we predict that inter-protein disorder-to-order transitions within the prion-like domains of intracellular proteins could lead to the required non-conservative enhancement of inter-molecular interactions. Our study, therefore, predicts a potential mechanismSignificance StatementBiomolecular condensates are highly diverse systems spanning not only homogeneous liquid droplets, but also gels, glasses, and even multiphase architectures that contain various coexisting liquid-like and/or gel-like inner phases. Multiphase architectures form when the different biomolecular components in a multi-component condensate establish sufficiently imbalanced inter-molecular forces to sustain different coexisting phases. While such a requirement seems, at first glance, impossible to fulfil for a condensate formed exclusively of chemically-identical proteins (i.e., single-component), our simulations predict conditions under which this may be possible. During condensate ageing, a sufficiently large imbalance in inter-molecular interactions can emerge intrinsically from the accumulation of protein structural transitions—driving even single-component condensates into nonequilibrium liquid-core/gel-shell or gel-core/liquid-shell multiphase architectures.

Multivalent proteins and nucleic acids can self-assemble into biomolecular condensates that contribute to compartmentalize the cell interior. Computer simulations offer a unique view to elucidate the mechanisms and key intermolecular interactions behind the dynamic formation and dissolution of these condensates. In this work, we present a novel approach to include explicit water and salt in sequence-dependent coarse-grained (CG) models for proteins and RNA, enabling the study of biomolecular condensate formation in a salt-dependent manner. Our framework combines a reparameterized version of the HPS protein force field with the monoatomic mW water model and the mW-ion potential for NaCl. We show how our CG model qualitatively captures the experimental radius of the gyration trend of a subset of intrinsically disordered proteins and reproduces the experimental protein concentration and water percentage of the human fused in sarcoma (FUS) low-complexity-domain droplets at physiological salt concentration. Moreover, we perform seeding simulations as a function of salt concentration for two antagonist systems: the engineered peptide PR25 and poly-uridine/poly-arginine mixtures, finding good agreement with their reported in vitro phase behavior with salt concentration in both cases. Taken together, our work represents a step forward towards extending sequence-dependent CG models to include water and salt, and to consider their key role in biomolecular condensate self-assembly.

One of the key mechanisms employed by cells to control their spatiotemporal organization is the formation and dissolution of phase-separated condensates. The balance between condensate assembly and disassembly can be critically regulated by the presence of RNA. In this work, we use a chemically-accurate sequence-dependent coarse-grained model for proteins and RNA to unravel the impact of RNA in modulating the transport properties and stability of biomolecular condensates. We explore the phase behavior of several RNA-binding proteins such as FUS, hnRNPA1, and TDP-43 proteins along with that of their corresponding prion-like domains and RNA recognition motifs from absence to moderately high RNA concentration. By characterizing the phase diagram, key molecular interactions, surface tension, and transport properties of the condensates, we report a dual RNA-induced behavior: on the one hand, RNA enhances phase separation at low concentration as long as the RNA radius of gyration is comparable to that of the proteins, whereas at high concentration, it inhibits the ability of proteins to self-assemble independently of its length. On the other hand, along with the stability modulation, the viscosity of the condensates can be considerably reduced at high RNA concentration as long as the length of the RNA chains is shorter than that of the proteins. Conversely, long RNA strands increase viscosity even at high concentration, but barely modify protein self-diffusion which mainly depends on RNA concentration and on the effect RNA has on droplet density. On the whole, our work rationalizes the different routes by which RNA can regulate phase separation and condensate dynamics, as well as the subsequent aberrant rigidification implicated in the emergence of various neuropathologies and age-related diseases.

In this work, a low-cost and rapid electrochemical resistive DNA biosensor based on the current relaxation method is described. A DNA probe, complementary to the specific human papillomavirus type 16 (HPV-16) sequence, was immobilized onto a screen-printed gold electrode. DNA hybridization was detected by applying a potential step of 30 mV to the system, composed of an external capacitor and the modified electrode DNA/gold, for 750 µs and then relaxed back to the OCP, at which point the voltage and current discharging curves are registered for 25 ms. From the discharging curves, the potential and current relaxation were evaluated, and by using Ohm’s law, the charge transfer resistance through the DNA-modified electrode was calculated. The presence of a complementary sequence was detected by the change in resistance when the ssDNA is transformed in dsDNA due to the hybridization event. The target DNA concentration was detected in the range of 5 to 20 nM. The results showed a good fit to the regression equation ΔRtotal(Ω)=2.99&nbsp, ×&nbsp, [DNA]+81.55, and a detection limit of 2.39 nM was obtained. As the sensing approach uses a direct current, the electronic architecture of the biosensor is simple and allows for the separation of faradic and nonfaradaic contributions. The simple electrochemical resistive biosensor reported here is a good candidate for the point-of-care diagnosis of HPV at a low cost and in a short detection time.

Intracellular liquid-liquid phase separation enables the formation of biomolecular condensates, such as ribonucleoprotein granules, which play a crucial role in the spatiotemporal organization of biomolecules (e.g., proteins and RNAs). Here, we introduce a patchy-particle polymer model to investigate liquid-liquid phase separation of protein-RNA mixtures. We demonstrate that at low to moderate concentrations, RNA enhances the stability of RNA-binding protein condensates because it increases the molecular connectivity of the condensed-liquid phase. Importantly, we find that RNA can also accelerate the nucleation stage of phase separation. Additionally, we assess how the capacity of RNA to increase the stability of condensates is modulated by the relative protein-protein/protein-RNA binding strengths. We find that phase separation and multiphase organization of multicomponent condensates is favored when the RNA binds with higher affinity to the lower-valency proteins in the mixture than to the cognate higher-valency proteins. Collectively, our results shed light on the roles of RNA in ribonucleoprotein granule formation and the internal structuring of stress granules.

Liquid–liquid phase separation of proteins underpins the formation of membraneless compartments in living cells. Elucidating the molecular driving forces underlying protein phase transitions is therefore a key objective for understanding biological function and malfunction. Here we show that cellular proteins, which form condensates at low salt concentrations, including FUS, TDP-43, Brd4, Sox2, and Annexin A11, can reenter a phase-separated regime at high salt concentrations. By bringing together experiments and simulations, we demonstrate that this reentrant phase transition in the high-salt regime is driven by hydrophobic and non-ionic interactions, and is mechanistically distinct from the low-salt regime, where condensates are additionally stabilized by electrostatic forces. Our work thus sheds light on the cooperation of hydrophobic and non-ionic interactions as general driving forces in the condensation process, with important implications for aberrant function, druggability, and material properties of biomolecular condensates., Elucidating the molecular driving forces underlying liquid–liquid phase separation is a key objective for understanding biological function and malfunction. Here the authors show that a wide range of cellular proteins, including FUS, TDP-43, Brd4, Sox2, and Annexin A11, which form condensates at low salt concentrations, can reenter a phase-separated regime at high salt concentrations.

Hard-sphere model systems are well-suited in both experiment and simulations to investigate fundamental aspects of the crystallization of fluids. In experiments on colloidal models of hard-sphere fluids, the fluid is unavoidably in contact with the walls of the sample cell, where heterogeneous crystallization may take place. In this work we use simulations to investigate the competition between homogeneous and heterogeneous crystallization. We report simulations of wall-induced nucleation for different confining walls. Combining the results of these simulations with earlier studies of homogeneous crystallization allows us to assess the competition between homogeneous and heterogeneous nucleation as a function of the wall type, fluid density and the system size. On flat walls, heterogeneous nucleation will typically overwhelm homogeneous nucleation. Even for surfaces randomly coated with spheres with a diameter that was some three times larger than that of the fluid spheres - as has been used in some experiments - heterogeneous nucleation is likely to be dominant for volume fractions smaller than ∼0.535. Only for a disordered coating that has the same structure as the liquid did we find that nucleation was likely to occur in the bulk. Hence, such coatings might be used to suppress heterogeneous nucleation in experiments. Finally, we report the apparent homogeneous nucleation rate taking into account the formation of crystallites both in the bulk and at the walls. We find that the apparent overall nucleation rates coincide with those reported in "homogeneous nucleation" experiments. This suggests that heterogeneous nucleation at the walls could partly explain the large discrepancies found between experimental measurements and simulation estimates of the homogeneous nucleation rate.

Colloidal particles have been extensively used to comprehend the main principles governing liquid-crystal nucleation. Multiple mechanisms and frameworks have been proposed, through either experiments or computational approaches, to rationalise the ubiquitous formation of colloidal crystals. In this work, we elucidate the nucleation scenario behind the crystallization of oppositely charged colloids. By performing molecular dynamics simulations of colloidal electrolytes in combination with the Seeding technique, we evaluate the fundamental factors, such as the nucleation rate, free energy barrier, surface tension and kinetic pre-factor, that determine the liquid-to-solid transition of several crystalline polymorphs. Our results show that at a high packing fraction, there is a cross-over between the nucleation of the CsCl structure and that of a substitutionally disordered fcc phase, despite the CuAu crystal being the most stable phase. We demonstrate that the crucial factor in determining which phase nucleates the fastest is the free energy cost of the cluster formation rather than their kinetic ability to grow from the liquid. While at a low packing fraction, the stable phase, CsCl, is the one that nucleates and subsequently grows, we show how at moderate and high packing fractions, a disordered fcc phase subsequently grows regardless of the nature of the nucleating phase, termed parasitic crystallization. Taken together, our results provide a panoramic perspective of the complex nucleation scenario of oppositely charged colloids at moderate temperature and rationalise the different thermodynamic and kinetic aspects behind it.

Proteins containing intrinsically disordered regions (IDRs) are ubiquitous within biomolecular condensates, which are liquid-like compartments within cells formed through liquid&ndash, 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&ndash, 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&ndash, 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&ndash, protein interactions. This study sheds light on the link between the dynamic conformational plasticity of IDRs and their liquid&ndash, liquid phase behaviour.

Many cellular proteins have the ability to demix spontaneously from solution to form liquid condensates. These phase-separated structures form membraneless compartments in living cells and have wide-ranging roles in health and disease. Elucidating the molecular driving forces underlying liquid–liquid phase separation (LLPS) of proteins has thus become a key objective for understanding biological function and malfunction. Here we show that proteins implicated in cellular phase separation, such as FUS, TDP-43, and Annexin A11, which form condensates at low salt concentrations via homotypic multivalent interactions, also have the ability to undergo LLPS at high salt concentrations by reentering into a phase-separated regime. Through a combination of experiments and simulations, we demonstrate that phase separation in the high-salt regime is mainly driven by hydrophobic and non-ionic interactions. As such, it is mechanistically distinct from the low-salt regime, where condensates are stabilized by a broad mix of electrostatic, hydrophobic, and non-ionic forces. Our work thus expands the molecular grammar of interactions governing LLPS of cellular proteins and provides a new view on hydrophobicity and non-ionic interactions as non-specific driving forces for the condensation process, with important implications for the aberrant function, druggability, and material properties of biomolecular condensates. One Sentence Summary Proteins implicated in cellular phase separation can undergo a salt-mediated reentrant liquid–liquid phase transition.

To predict the radiative forcing of clouds it is necessary to know the rate with which ice homogeneously nucleates in supercooled water. Such rate is often measured in drops to avoid the presence of impurities. At large supercooling small (nanoscopic) drops must be used to prevent simultaneous nucleation events. The pressure inside such drops is larger than the atmospheric one by virtue of the Laplace equation. In this work, we take into account such pressure raise in order to predict the nucleation rate in droplets using the TIP4P/Ice water model. We start from a recent estimate of the maximum drop size that can be used at each supercooling avoiding simultaneous nucleation events [Espinosa et al. J. Chem. Phys., 2016]. We then evaluate the pressure inside the drops with the Laplace equation. Finally, we obtain the rate as a function of the supercooling by interpolating our previous results for 1 and 2000 bar [Espinosa et al. Phys. Rev. Lett. 2016] using the Classical Nucleation Theory expression for the rate. This requires, in turn, interpolating the ice-water interfacial free energy and chemical potential difference. The TIP4P/Ice rate curve thus obtained is in good agreement with most droplet-based experiments. In particular, we find a good agreement with measurements performed using nanoscopic drops, that are currently under debate. The successful comparison between model and experiments suggests that TIP4P/Ice is a reliable model to study the water-to-ice transition and that Classical Nucleation Theory is a good framework to understand it., paper