Variational umbrella seeding for calculating nucleation barriers.

In this work, we introduce variational umbrella seeding, a novel technique for computing nucleation barriers. This new method, a refinement of the original seeding approach, is far less sensitive to the choice of order parameter for measuring the size of a nucleus. Consequently, it surpasses seeding in accuracy and umbrella sampling in computational speed. We test the method extensively and demonstrate excellent accuracy for crystal nucleation of nearly hard spheres and two distinct models of water: mW and TIP4P/ICE. This method can easily be extended to calculate nucleation barriers for homogeneous melting, condensation, and cavitation.

Synergistic epistasis among cancer drivers can rescue early tumors from the accumulation of deleterious passengers.

Epistasis among driver mutations is pervasive and explains relevant features of cancer, such as differential therapy response and convergence towards well-characterized molecular subtypes. Furthermore, a growing body of evidence suggests that tumor development could be hampered by the accumulation of slightly deleterious passenger mutations. In this work, we combined empirical epistasis networks, computer simulations, and mathematical models to explore how synergistic interactions among driver mutations affect cancer progression under the burden of slightly deleterious passengers. We found that epistasis plays a crucial role in tumor development by promoting the transformation of precancerous clones into rapidly growing tumors through a process that is analogous to evolutionary rescue. The triggering of epistasis-driven rescue is strongly dependent on the intensity of epistasis and could be a key rate-limiting step in many tumors, contributing to their unpredictability. As a result, central genes in cancer epistasis networks appear as key intervention targets for cancer therapy.

Relationships between Patient-Reported Outcome Measures and Clinical Measures in Naïve Neovascular Age-Related Macular Degeneration Patients Treated with Intravitreal Ranibizumab.

Our objective was to evaluate changes in patient-reported outcome measures using the NEI-VFQ 25 questionnaire during a treat and extend regimen in naive neovascular Age-Related Macular Degeneration patients, and its correlation with anatomical and functional data. We conducted a prospective observational study. Patients underwent a treat and extend regimen with intravitreal ranibizumab for neovascular Age-Related Macular Degeneration. Initial response was evaluated at 4th month, and subsequently in every follow-up visit. If a clinical response was achieved, the injection interval was extended in two-week increments, up to a maximum of 12 weeks. Quality of life was assessed using the NEI-VFQ 25 questionnaire at baseline, 4th months, and 12th months. Patients were categorized as good or poor responders based on Best corrected visual acuity, central foveal thickness, intraretinal fluid, or subretinal fluid. Treatment with ranibizumab led to a significant improvement in quality of life, with a mean increase in NEI-VFQ 25 score of 4.27 points in the 12th month. No significant differences in improvement were observed between good and poor responders. Quality of life scores in neovascular Age-Related Macular Degeneration patients improved with intravitreal treatment regardless of the clinical response. The early response following the loading phase could indicate better quality of life after one year of treatment, with Best corrected visual acuity being the clinical parameter with the greatest influence on quality of life.

Forehead Reduction Surgery with Hairline Advancement.

Hairline reduction surgery, also known as aesthetic forehead reduction, is a surgical procedure that aims to reduce the upper facial third and improve facial harmony. This article describes the anatomy of hairline advancement surgery and the surgical technique used by the author. The study included 21 patients from 2019 to 2023, and the forehead reduction length was on average 22.48 ± 2.64 mm. The most common complaint was hypoesthesia of the scalp, present in 100% of patients, resolving in all cases by 2 months after surgery. Forehead reduction surgery is among the procedures that provide more satisfaction to patients due to the great changes that it can achieve.

Principles of assembly and regulation of condensates of Polycomb repressive complex 1 through phase separation.

Polycomb repressive complex 1 (PRC1) undergoes phase separation to form Polycomb condensates that are multi-component hubs for silencing Polycomb target genes. In this study, we demonstrate that formation and regulation of PRC1 condensates are consistent with the scaffold-client model, where the Chromobox 2 (CBX2) protein behaves as the scaffold while the other PRC1 proteins are clients. Such clients induce a re-entrant phase transition of CBX2 condensates. The composition of the multi-component PRC1 condensates (1) determines the dynamic properties of the scaffold protein; (2) selectively promotes the formation of CBX4-PRC1 condensates while dissolving condensates of CBX6-, CBX7-, and CBX8-PRC1; and (3) controls the enrichment of CBX4-, CBX7-, and CBX8-PRC1 in CBX2-PRC1 condensates and the exclusion of CBX6-PRC1 from CBX2-PRC1 condensates. Our findings uncover how multi-component PRC1 condensates are assembled via an intricate scaffold-client mechanism whereby the properties of the PRC1 condensates are sensitively regulated by its composition and stoichiometry.

Peptide-conjugated antimiRs improve myotonic dystrophy type 1 phenotypes by promoting endogenous MBNL1 expression.

Myotonic dystrophy type 1 (DM1) is a rare neuromuscular disease caused by a CTG repeat expansion in the DMPK gene that generates toxic RNA with a myriad of downstream alterations in RNA metabolism. A key consequence is the sequestration of alternative splicing regulatory proteins MBNL1/2 by expanded transcripts in the affected tissues. MBNL1/2 depletion interferes with a developmental alternative splicing switch that causes the expression of fetal isoforms in adults. Boosting the endogenous expression of MBNL proteins by inhibiting the natural translational repressors miR-23b and miR-218 has previously been shown to be a promising therapeutic approach. We designed antimiRs against both miRNAs with a phosphorodiamidate morpholino oligonucleotide (PMO) chemistry conjugated to cell-penetrating peptides (CPPs) to improve delivery to affected tissues. In DM1 cells, CPP-PMOs significantly increased MBNL1 levels. In some candidates, this was achieved using concentrations less than two orders of magnitude below the median toxic concentration, with up to 5.38-fold better therapeutic window than previous antagomiRs. In HSALR mice, intravenous injections of CPP-PMOs improve molecular, histopathological, and functional phenotypes, without signs of toxicity. Our findings place CPP-PMOs as promising antimiR candidates to overcome the treatment delivery challenge in DM1 therapy.

The liquid-to-solid transition of FUS is promoted by the condensate surface.

A wide range of macromolecules can undergo phase separation, forming biomolecular condensates in living cells. These membraneless organelles are typically highly dynamic, formed reversibly, and carry out essential functions in biological systems. Crucially, however, a further liquid-to-solid transition of the condensates can lead to irreversible pathological aggregation and cellular dysfunction associated with the onset and development of neurodegenerative diseases. Despite the importance of this liquid-to-solid transition of proteins, the mechanism by which it is initiated in normally functional condensates is unknown. Here we show, by measuring the changes in structure, dynamics, and mechanics in time and space, that single-component FUS condensates do not uniformly convert to a solid gel, but rather that liquid and gel phases coexist simultaneously within the same condensate, resulting in highly inhomogeneous structures. Furthermore, our results show that this transition originates at the interface between the condensate and the dilute continuous phase, and once initiated, the gelation process propagates toward the center of the condensate. To probe such spatially inhomogeneous rheology during condensate aging, we use a combination of established micropipette aspiration experiments together with two optical techniques, spatial dynamic mapping and reflective confocal dynamic speckle microscopy. These results reveal the importance of the spatiotemporal dimension of the liquid-to-solid transition and highlight the interface of biomolecular condensates as a critical element in driving pathological protein aggregation.

Mechano-dependent sorbitol accumulation supports biomolecular condensate.

Biomolecular condensates regulate a wide range of cellular functions from signaling to RNA metabolism 1, 2 , yet, the physiologic conditions regulating their formation remain largely unexplored. Biomolecular condensate assembly is tightly regulated by the intracellular environment. Changes in the chemical or physical conditions inside cells can stimulate or inhibit condensate formation 3-5 . However, whether and how the external environment of cells can also regulate biomolecular condensation remain poorly understood. Increasing our understanding of these mechanisms is paramount as failure to control condensate formation and dynamics can lead to many diseases 6, 7 . Here, we provide evidence that matrix stiffening promotes biomolecular condensation in vivo . We demonstrate that the extracellular matrix links mechanical cues with the control of glucose metabolism to sorbitol. In turn, sorbitol acts as a natural crowding agent to promote biomolecular condensation. Using in silico simulations and in vitro assays, we establish that variations in the physiological range of sorbitol, but not glucose, concentrations, are sufficient to regulate biomolecular condensates. Accordingly, pharmacologic and genetic manipulation of intracellular sorbitol concentration modulates biomolecular condensates in breast cancer - a mechano-dependent disease. We propose that sorbitol is a mechanosensitive metabolite enabling protein condensation to control mechano-regulated cellular functions. Altogether, we uncover molecular driving forces underlying protein phase transition and provide critical insights to understand the biological function and dysfunction of protein phase separation.

Location and Concentration of Aromatic-Rich Segments Dictates the Percolating Inter-Molecular Network and Viscoelastic Properties of Ageing Condensates.

Maturation of functional liquid-like biomolecular condensates into solid-like aggregates has been linked to the onset of several neurodegenerative disorders. Low-complexity aromatic-rich kinked segments (LARKS) contained in numerous RNA-binding proteins can promote aggregation by forming inter-protein ß-sheet fibrils that accumulate over time and ultimately drive the liquid-to-solid transition of the condensates. Here, atomistic molecular dynamics simulations are combined with sequence-dependent coarse-grained models of various resolutions to investigate the role of LARKS abundance and position within the amino acid sequence in the maturation of condensates. Remarkably, proteins with tail-located LARKS display much higher viscosity over time than those in which the LARKS are placed toward the center. Yet, at very long timescales, proteins with a single LARKS-independently of its location-can still relax and form high viscous liquid condensates. However, phase-separated condensates of proteins containing two or more LARKS become kinetically trapped due to the formation of percolated ß-sheet networks that display gel-like behavior. Furthermore, as a work case example, they demonstrate how shifting the location of the LARKS-containing low-complexity domain of FUS protein toward its center effectively precludes the accumulation of ß-sheet fibrils in FUS-RNA condensates, maintaining functional liquid-like behavior without ageing.

Time-Dependent Material Properties of Aging Biomolecular Condensates from Different Viscoelasticity Measurements in Molecular Dynamics Simulations.

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 behavior of their molecules, characterizing such properties is integral to rationalizing 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 Green-Kubo (GK) relation, 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 GK and OS techniques outperform the BT method in terms of computational efficiency and statistical uncertainty. We thus apply the GK and OS techniques for a set of 12 different protein/RNA systems using a sequence-dependent 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 stickers vs spacers in the amino acid protein sequence. Moreover, we couple the GK 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 interprotein ß-sheets. We compare the behavior 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 the GK and OS techniques successfully predict the transition from functional liquid-like behavior to kinetically arrested states once the network of interprotein ß-sheets has percolated through the condensates. Overall, our work provides a comparison of different modeling rheological techniques to assess the viscosity of biomolecular condensates, a critical magnitude that provides information on the behavior of biomolecules inside condensates.

A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water.

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.

On the possible locus of the liquid-liquid critical point in real water from studies of supercooled water using the TIP4P/Ice model.

One of the most accepted hypothesis to explain the anomalous behavior of water is the presence of a critical point between two liquids, the liquid-liquid critical point (LLCP), buried within the deep supercooled regime. Unfortunately, such hypothesis is hard to be experimentally confirmed due to fast freezing. Here, we show that the TIP4P/Ice water potential shifted by 400 bar can reproduce with unprecedented accuracy the experimental isothermal compressibility of water and its liquid equation of state for a wide pressure and temperature range. We find, both by extrapolation of response function maxima and by a Maxwell construction, that the location of the model LLCP is consistent with previous calculations. According to the pressure shift needed to recover the experimental behavior of supercooled water, we estimate the experimental LLCP to be located around 1250 bar and 195 K. We use the model to estimate the ice nucleation rate (J) in the vicinity of the hypothesized LLCP experimental location and obtain J = 1024 m-3 s-1. Thereby, experiments where the ratio between the cooling rate and the sample volume is equal or larger than the estimated nucleation rate could probe liquid-liquid equilibrium before freezing. Such conditions are not accessible in common experiments with microdroplets cooled at a few kelvin per second, but they could be, for instance, using nanodroplets of around 50 nm radius observed in a millisecond timescale.

BlockmiR AONs as Site-Specific Therapeutic MBNL Modulation in Myotonic Dystrophy 2D and 3D Muscle Cells and HSALR Mice.

The symptoms of Myotonic Dystrophy Type 1 (DM1) are multi-systemic and life-threatening. The neuromuscular disorder is rooted in a non-coding CTG microsatellite expansion in the DM1 protein kinase (DMPK) gene that, upon transcription, physically sequesters the Muscleblind-like (MBNL) family of splicing regulator proteins. The high-affinity binding occurring between the proteins and the repetitions disallow MBNL proteins from performing their post-transcriptional splicing regulation leading to downstream molecular effects directly related to disease symptoms such as myotonia and muscle weakness. In this study, we build on previously demonstrated evidence showing that the silencing of miRNA-23b and miRNA-218 can increase MBNL1 protein in DM1 cells and mice. Here, we use blockmiR antisense technology in DM1 muscle cells, 3D mouse-derived muscle tissue, and in vivo mice to block the binding sites of these microRNAs in order to increase MBNL translation into protein without binding to microRNAs. The blockmiRs show therapeutic effects with the rescue of mis-splicing, MBNL subcellular localization, and highly specific transcriptomic expression. The blockmiRs are well tolerated in 3D mouse skeletal tissue inducing no immune response. In vivo, a candidate blockmiR also increases Mbnl1/2 protein and rescues grip strength, splicing, and histological phenotypes.

Direct Calculation of the Interfacial Free Energy between NaCl Crystal and Its Aqueous Solution at the Solubility Limit.

Salty water is the most abundant electrolyte aqueous mixture on Earth, however, very little is known about the NaCl-saturated solution interfacial free energy (γ_{s}). Here, we provide the first direct estimation of γ_{s} for several NaCl crystallographic planes by means of the mold integration technique, a highly efficient computational method to evaluate interfacial free energies with anisotropic crystal resolution. Making use of the JC-SPC/E model, one of the most benchmarked force fields for NaCl water solutions, we measure γ_{s} of four different crystal planes, (100), (110), (111), and (112[over ¯]) with the saturated solution at normal conditions. We find high anisotropy between the different crystal orientations with values ranging from 100 to 150 mJ m^{-2}, and the average value of the distinct planes being γ[over ¯]_{s}=137(20) mJ m^{-2}. This value for the coexistence interfacial free energy is in reasonable agreement with previous extrapolations from nucleation studies. Our Letter represents a milestone in the computational calculation of interfacial free energies between ionic crystals and aqueous solutions.

Surfactants or scaffolds? RNAs of varying lengths control the thermodynamic stability of condensates differently.

Biomolecular condensates, thought to form via liquid-liquid phase separation of intracellular mixtures, are multicomponent systems that can include diverse types of proteins and RNAs. RNA is a critical modulator of RNA-protein condensate stability, as it induces an RNA concentration-dependent reentrant phase transition-increasing stability at low RNA concentrations and decreasing it at high concentrations. Beyond concentration, RNAs inside condensates can be heterogeneous in length, sequence, and structure. Here, we use multiscale simulations to understand how different RNA parameters interact with one another to modulate the properties of RNA-protein condensates. To do so, we perform residue/nucleotide resolution coarse-grained molecular dynamics simulations of multicomponent RNA-protein condensates containing RNAs of different lengths and concentrations, and either FUS or PR25 proteins. Our simulations reveal that RNA length regulates the reentrant phase behavior of RNA-protein condensates: increasing RNA length sensitively rises the maximum value that the critical temperature of the mixture reaches, and the maximum concentration of RNA that the condensate can incorporate before beginning to become unstable. Strikingly, RNAs of different lengths are organized heterogeneously inside condensates, which allows them to enhance condensate stability via two distinct mechanisms: shorter RNA chains accumulate at the condensate's surface acting as natural biomolecular surfactants, while longer RNA chains concentrate inside the core to saturate their bonds and enhance the density of molecular connections in the condensate. Using a patchy particle model, we additionally demonstrate that the combined impact of RNA length and concentration on condensate properties is dictated by the valency, binding affinity, and polymer length of the various biomolecules involved. Our results postulate that diversity on RNA parameters within condensates allows RNAs to increase condensate stability by fulfilling two different criteria: maximizing enthalpic gain and minimizing interfacial free energy; hence, RNA diversity should be considered when assessing the impact of RNA on biomolecular condensates regulation.

Protocol for comparing gene-level selection on coding mutations between two groups of samples with Coselens.

The study of genes that evolve under conditional selection can shed light on the genomic underpinnings of adaptation, revealing epistasis and phenotypic plasticity. This protocol describes how to use the Coselens package to compare gene-level selection between two groups of samples. After installing Coselens and preparing the datasets, a typical run on a laptop takes less than 10 min. Coselens is best suited to analyze somatic mutations and data from experimental evolution, for which independently evolved samples are available. For complete details on the use and execution of this protocol, please refer to Iranzo et al. (2022).1.

CRISPR-Cas9 editing of a TNPO3 mutation in a muscle cell model of limb-girdle muscular dystrophy type D2.

A single-nucleotide deletion in the stop codon of the nuclear import receptor transportin-3 (TNPO3), also involved in human immunodeficiency virus type 1 (HIV-1) infection, causes the ultrarare autosomal dominant disease limb-girdle muscular dystrophy D2 (LGMDD2) by extending the wild-type protein. Here, we generated a patient-derived in vitro model of LGMDD2 as an immortalized myoblast cell line carrying the TNP O 3 mutation. The cell model reproduced critical molecular alterations seen in patients, such as TNP O 3 overexpression, defects in terminal muscle markers, and autophagy overactivation. Correction of the TNP O 3 mutation via CRISPR-Cas9 editing caused a significant reversion of the pathological phenotypes in edited cells, including a complete absence of the mutant TNPO3 protein, as detected with a polyclonal antibody specific against the abnormal 15-aa peptide. Transcriptomic analyses found that 15% of the transcriptome was differentially expressed in model myotubes. CRISPR-Cas9-corrected cells showed that 44% of the alterations were rescued toward normal levels. MicroRNAs (miRNAs) analyses showed that around 50% of miRNAs with impaired expression because of the disease were recovered on the mutation edition. In summary, this work provides proof of concept of the potential of CRISPR-Cas9-mediated gene editing of TNP O 3 as a therapeutic approach and describes critical reagents in LGMDD2 research.

The Chromatin Regulator HMGA1a Undergoes Phase Separation in the Nucleus.

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.