Toward Accurate Simulation of Coupling between Protein Secondary Structure and Phase Separation

Benchmark of Coacervate Formation and Mechanism Exploration Using the Martini Force Field

Peptide-based coacervates are crucial for drug delivery due to their biocompatibility, versatility, high drug loading capacity, and cell penetration rates; however, their stability mechanism and phase behavior are not fully understood. Additionally, although Martini is one of the most famous force fields capable of describing coacervate formation with molecular details, a comprehensive benchmark of its accuracy has not been conducted. This research utilized the Martini 3.0 force field and machine learning algorithms to explore representative peptide-based coacervates, including those composed of polyaspartate (PAsp)/polyarginine (PArg), rmfp-1, sticker-and-spacer small molecules, and HBpep molecules. We identified key coacervate formation driving forces such as Coulomb, cation-π, and π-π interactions and established three criteria for determining coacervate formation in simulations. The results also indicate that while Martini 3.0 accurately captures coacervate formation trends, it tends to underestimate Coulomb interactions and overestimate π-π interactions. What is more, our study on drug encapsulation of HBpep and its derivative coacervates suggested that most loaded drugs were distributed on surfaces of HBpep clusters, awaiting experimental validation. This study employs simulation to enhance understanding of peptide-based coacervate phase behavior and stability mechanisms while also benchmarking Martini 3.0, thereby providing fundamental insights for future experimental and simulation investigations.

Tracking flaviviral protease conformational dynamics by tuning single-molecule nanopore tweezers

The flaviviral NS2B/NS3 protease is a conserved enzyme required for flavivirus replication. Its highly dynamic conformation poses major challenges but also offers opportunities for antiviral inhibition. Here, we established a nanopore tweezers-based platform to monitor NS2B/NS3 conformational dynamics in real time. Molecular simulations coupled with single-channel current recording measurements revealed that the protease could be captured in the middle of the ClyA nanopore lumen, stabilized mainly by dynamic electrostatic interactions. We designed a new Salmonella typhi ClyA nanopore with enhanced nanopore/protease interaction that can resolve the open and closed states at the single-molecule level for the first time. We demonstrated that the tailored ClyA could track the conformational transitions of the West Nile NS2B/NS3 protease and unravel the conformational energy landscape of various protease constructs through population and kinetic analysis. The new ClyA-protease platform paves a way to search for new allosteric inhibitors that target the NS2B and NS3 interface.

Microscopic Origins of Flow Activation Energy in Biomolecular Condensates

The material properties of biomolecular condensates govern their dynamics and functions by influencing the molecular diffusion rates and biochemical interactions. A recent report has identified a characteristic timescale of temperature-dependent viscosity in biomolecular condensates arising from an activated dissociation events collectively referred to as flow activation energy. The microscopic origin of this activation energy is a complex function of sequence, stoichiometry, and external conditions. In this study, we elucidate the microscopic origins of flow activation energy in single and multicomponent condensates formed from model peptide sequences with varying "sticker" and "spacer" motifs, incorporating RNA as a secondary component. We examined how condensate density, RNA stoichiometry, and peptide sequence patterning impact these properties through detailed sequence-resolved coarse-grained simulations. Our findings reveal that flow activation energy is closely tied to the lifetime of sticker-sticker interactions under specific conditions. However, the presence of multiple competing stickers may complicate this relationship leading to frustrated interactions in condensates and lowering of activation energy. The findings of this study should help to create predictive models of material properties of condensates, which in turn can facilitate a more profound understanding of functions and programmable design principles of biomolecular condensates.

An Intermediate Resolution Model of RNA Dynamics and Phase Separation with Explicit Mg 2 +

RNAs are major drivers of phase separation in the formation of biomolecular condensates. Recent studies suggest that RNAs can also undergo protein-free phase separation in the presence of divalent ions or crowding agents. Much remains to be understood regarding how the complex interplay of base stacking, base pairing, electrostatics, ion interactions, and structural propensities governs the phase behaviour of RNAs. Here we develop an intermediate resolution model for condensates of RNAs (iConRNA) that can capture key local and long-range structure features of dynamic RNAs and simulate their spontaneous phase transitions in the presence ofMg 2 + . Representing each nucleotide using 6 or 7 beads, iConRNA considers specific RNA base stacking and pairing interactions and includes explicitMg 2 + ions to studyMg 2 + -induced phase separation. Parametrized using theoretical and experimental data, the model can correctly reproduce the chain properties of A-form helical poly(rA) and coil poly(rU), and essential structures of several RNA hairpins. With an effectiveMg 2 + ion model, iConRNA simulations successfully recapitulate the experimentally observed lower critical solution temperature (LCST)-type phase separation of poly(rA) and the dissolution of poly(rU). Furthermore, the phase diagrams of CAG/CUG/CUU-repeat RNAs correctly reproduce the experimentally observed sequence- and length-dependence of phase separation propensity. These results suggest that iConRNA can be a viable tool for studying homotypic RNA and potentially heterotypic RNA-protein phase separations.

A coarse-grained model for disordered and multi-domain proteins

Many proteins contain more than one folded domain, and such modular multi-domain proteins help expand the functional repertoire of proteins. Because of their larger size and often substantial dynamics, it may be difficult to characterize the conformational ensembles of multi-domain proteins by simulations. Here, we present a coarse-grained model for multi-domain proteins that is both fast and provides an accurate description of the global conformational properties in solution. We show that the accuracy of a one-bead-per-residue coarse-grained model depends on how the interaction sites in the folded domains are represented. Specifically, we find excessive domain-domain interactions if the interaction sites are located at the position of the Cα atoms. We also show that if the interaction sites are located at the center of mass of the residue, we obtain good agreement between simulations and experiments across a wide range of proteins. We then optimize our previously described CALVADOS model using this center-of-mass representation, and validate the resulting model using independent data. Finally, we use our revised model to simulate phase separation of both disordered and multi-domain proteins, and to examine how the stability of folded domains may differ between the dilute and dense phases. Our results provide a starting point for understanding interactions between folded and disordered regions in proteins, and how these regions affect the propensity of proteins to self-associate and undergo phase separation.

Microscopic Origins of Flow Activation Energy in Biomolecular Condensates

Material properties of biomolecular condensates dictate their form and function, influencing the diffusion of regulatory molecules and the dynamics of biochemical reactions. The increasing quality and quantity of microrheology experiments on biomolecular condensates necessitate a deeper understanding of the molecular grammar that encodes their material properties. Recent reports have identified a characteristic timescale related to network relaxation dynamics in condensates, which governs their temperature-dependent viscoelastic properties. This timescale is intimately connected to an activated process involving the dissociation of sticker regions, with the energetic barrier referred to as flow activation energy. The microscopic origin of activation energy is a complex function of sequence patterns, component stoichiometry, and external conditions. This study elucidates the microscopic origins of flow activation energy in single and multicomponent condensates composed of model peptide sequences with varying sticker and spacer motifs, with RNA as a secondary component. We dissected the effects of condensate density, RNA stoichiometry, and peptide sequence patterning using extensive sequence-resolved coarse-grained simulations. We found that flow activation energy is closely linked to the lifetime of sticker-sticker pairs under certain conditions, though the presence of multiple competing stickers further complicates this relationship. The insights gained in this study should help establish predictive multiscale models for the material properties and serve as a valuable guide for the programmable design of condensates.

Unraveling mechanisms of protein encapsulation and release in coacervates via molecular dynamics and machine learning

Coacervates play a pivotal role in protein-based drug delivery research, yet their drug encapsulation and release mechanisms remain poorly understood. Here, we utilized the Martini model to investigate bovine serum albumin (BSA) protein encapsulation and release within polylysine/polyglutamate (PLys/PGlu) coacervates. Our findings emphasize the importance of ingredient addition sequence in coacervate formation and encapsulation rates, attributed to preference contact between oppositely charged proteins and poly(amino acid)s. Notably, coacervates composed of β-sheet poly(amino acid)s demonstrate greater BSA encapsulation efficiency due to their reduced entropy and flexibility. Furthermore, we examined the pH responsiveness of coacervates, shedding light on the dissolution process driven by Coulomb forces. By leveraging machine learning algorithms to analyze simulation results, our research advances the understanding of coacervate-based drug delivery systems, with the ultimate goal of optimizing therapeutic outcomes.

Multiscale simulations reveal the driving forces of p53C phase separation accelerated by oncogenic mutations

Liquid-Liquid phase separation (LLPS) of p53 to form liquid condensates has been implicated in cellular functions and dysfunctions. The p53 condensates may serve as amyloid fibril precursors to initiate p53 aggregation, which is associated with oncogenic gain-of-function and various human cancers. M237I and R249S mutations located in p53 core domain (p53C) have been detected respectively in glioblastomas and hepatocellular carcinoma. Interestingly, these p53C mutants can also undergo LLPS and liquid-to-solid phase transition, which are faster than wild type p53C. However, the underlying molecular basis governing the accelerated LLPS and liquid-to-solid transition of p53C remain poorly understood. Herein, we explore the M237I/R249S mutation-induced structural alterations and phase separation behavior of p53C by employing multiscale molecular dynamics simulations. All-atom simulations revealed conformational disruptions in the zinc-binding domain of the M237I mutant and in both loop3 and zinc-binding domain of the R249S mutant. The two mutations enhance hydrophobic exposure of those regions and attenuate intramolecular interactions, which may hasten the LLPS and aggregation of p53C. Martini 3 coarse-grained simulations demonstrated spontaneous phase separation of p53C and accelerated effects of M237I/R249S mutations on the phase separation of p53C. Importantly, we find that the regions with enhanced intermolecular interactions observed in coarse-grained simulations coincide with the disrupted regions with weakened intramolecular interactions observed in all-atom simulations, indicating that M237I/R249S mutation-induced local structural disruptions expedite the LLPS of p53C. This study unveils the molecular mechanisms underlying the two cancer-associated mutation-accelerated LLPS and aggregation of p53C, providing avenues for anticancer therapy by targeting the phase separation process.

Fundamental Aspects of Phase-Separated Biomolecular Condensates

Biomolecular condensates, formed through phase separation, are upending our understanding in much of molecular, cell, and developmental biology. There is an urgent need to elucidate the physicochemical foundations of the behaviors and properties of biomolecular condensates. Here we aim to fill this need by writing a comprehensive, critical, and accessible review on the fundamental aspects of phase-separated biomolecular condensates. We introduce the relevant theoretical background, present the theoretical basis for the computation and experimental measurement of condensate properties, and give mechanistic interpretations of condensate behaviors and properties in terms of interactions at the molecular and residue levels.

Vitamin C Drives Reentrant Actin Phase Transition: Biphasic Exocytosis Regulation Revealed by Single-Vesicle Electrochemistry

Unveiling molecular mechanisms that dominate protein phase dynamics has been a pressing need for deciphering the intricate intracellular modulation machinery. While ions and biomacromolecules have been widely recognized for modulating protein phase separations, effects of small molecules that essentially constitute the cytosolic chemical atmosphere on the protein phase behaviors are rarely understood. Herein, we report that vitamin C (VC), a key small molecule for maintaining a reductive intracellular atmosphere, drives reentrant phase transitions of myosin II/F-actin (actomyosin) cytoskeletons. The actomyosin bundle condensates dissemble in the low-VC regime and assemble in the high-VC regime in vitro or inside neuronal cells, through a concurrent myosin II protein aggregation-dissociation process with monotonic VC concentration increase. Based on this finding, we employ in situ single-cell and single-vesicle electrochemistry to demonstrate the quantitative modulation of catecholamine transmitter vesicle exocytosis by intracellular VC atmosphere, i.e., exocytotic release amount increases in the low-VC regime and decreases in the high-VC regime. Furthermore, we show how VC regulates cytomembrane-vesicle fusion pore dynamics through counteractive or synergistic effects of actomyosin phase transitions and the intracellular free calcium level on membrane tensions. Our work uncovers the small molecule-based reversive protein phase regulatory mechanism, paving a new way to chemical neuromodulation and therapeutic repertoire expansion.

Tuning single-molecule ClyA nanopore tweezers for real-time tracking of the conformational dynamics of West Nile viral NS2B/NS3 protease

The flaviviral NS2B/NS3 protease is a conserved enzyme required for flavivirus replication. Its highly dynamic conformation poses major challenges but also offers opportunities for antiviral inhibition. Here, we established a nanopore tweezers-based platform to monitor NS2B/NS3 conformational dynamics in real-time. Molecular simulations coupled with electrophysiology revealed that the protease could be captured in the middle of the ClyA nanopore lumen, stabilized mainly by dynamic electrostatic interactions. We designed a new Salmonella typhi ClyA nanopore with enhanced nanopore/protease interaction that can resolve the open and closed states at the single-molecule level for the first time. We demonstrated that the tailored ClyA could track the conformational transitions of the West Nile NS2B/NS3 protease and unravel the conformational energy landscape of various protease constructs through population and kinetic analysis. The new ClyA-protease platform paves a way to high-throughput screening strategies for discovering new allosteric inhibitors that target the NS2B and NS3 interface.

Backbone interactions and secondary structures in phase separation of disordered proteins

Intrinsically disordered proteins (IDPs) are one of the major drivers behind the formation and characteristics of biomolecular condensates. Due to their inherent flexibility, the backbones of IDPs are significantly exposed, rendering them highly influential and susceptible to biomolecular phase separation. In densely packed condensates, exposed backbones have a heightened capacity to interact with neighboring protein chains, which might lead to strong coupling between the secondary structures and phase separation and further modulate the subsequent transitions of the condensates, such as aging and fibrillization. In this mini-review, we provide an overview of backbone-mediated interactions and secondary structures within biomolecular condensates to underscore the importance of protein backbones in phase separation. We further focus on recent advances in experimental techniques and molecular dynamics simulation methods for probing and exploring the roles of backbone interactions and secondary structures in biomolecular phase separation involving IDPs.