Benchmarking Molecular Dynamics Force Fields for All-Atom Simulations of Biological Condensates

TIDGN: A Transfer Learning Framework for Predicting Interactions of Intrinsically Disordered Proteins with High Conformational Dynamics

Interactions between intrinsically disordered proteins (IDPs) are crucial for biological processes, such as intracellular liquid-liquid phase separation (LLPS). Experiments (e.g., NMR) and simulations used to study IDP interactions encounter a variety of difficulties, highlighting the necessity to develop relevant machine learning methods. However, reliable machine learning methods face the challenge resulting from the scarcity of available training data. In this work, we propose a transfer learning-based invariant geometric dynamic graph model, named TIDGN, for predicting IDP interactions. The model consists of a pretraining task module and a downstream task module. The pretraining task module learns the dynamic structural encoding of IDP monomers, which is then used by the downstream task module for interaction site prediction. The IDP monomer structure data set and the IDP interaction event data set are constructed using all-atom molecular dynamics (MD) simulations. The transfer learning strategy effectively enhances the model's performance. Both homotypic interactions and heterotypic interactions between two IDPs are considered in this work. Interestingly, TIDGN performs well for the heterotypic interaction prediction. Additionally, the feature ablation analysis emphasizes the importance of invariant geometric graph features. Taken together, our work demonstrates that the integration of transfer learning and the invariant geometric graph network offers a promising approach for addressing data scarcity challenges of IDP interaction prediction.

Toward Predictive Coarse-Grained Simulations of Biomolecular Condensates

Phase separation is a fundamental process that enables cellular organization by forming biomolecular condensates. These assemblies regulate diverse functions by creating distinct environments, influencing reaction kinetics, and facilitating processes such as genome organization, signal transduction, and RNA metabolism. Recent studies highlight the complexity of condensate properties, shaped by intrinsic molecular features and external factors such as temperature and pH. Molecular simulations serve as an effective approach to establishing a comprehensive framework for analyzing these influences, offering high-resolution insights into condensate stability, dynamics, and material properties. This review evaluates recent advancements in biomolecular condensate simulations, with a particular focus on coarse-grained 1-bead-per-amino-acid (1BPA) protein models, and emphasizes OpenABC, a tool designed to simplify and streamline condensate simulations. OpenABC supports the implementation of various coarse-grained force fields, enabling their performance evaluation. Our benchmarking identifies inconsistencies in phase behavior predictions across force fields, even though these models accurately capture single-chain statistics. This finding underscores the need for enhanced force field accuracy, achievable through enriched training data sets, many-body potentials, and advanced optimization techniques. Such refinements could significantly improve the predictive capacity of coarse-grained models, bridging molecular details with emergent condensate behaviors.

Insight into the thermo-responsive phase behavior of the P1 domain of α-synuclein using atomistic simulations

Biomolecular condensate formation driven by intrinsically disordered proteins (IDPs) is regulated by interactions between various domains of the proteins. Such condensates are implicated in various neurodegenerative diseases. The presynaptic intrinsically disordered protein, α-Syn is involved in the pathogenesis of Parkinson's disease. The central non-amyloid β-component (NAC) domain in the protein is considered to be a major driver of pathogenic aggregation, although recent studies have suggested that the P1 domain from the flanking N-terminal region can act as a 'master controller' for α-Syn function and aggregation. To gain molecular insight into the phase behavior of the P1 domain itself, we investigate how assemblies of P1 (residues 36-42) chains phase separate with varying temperatures using all-atom molecular dynamics simulations. The simulations reveal that P1 is able to phase separate above a lower critical solution temperature. Formation of a condensed phase is driven by exclusion of water molecules by the hydrophobic chains. P1 chain density in the condensate is determined by weak multi-chain interactions between the residues. Moreover, tyrosine (Y39) is involved in the formation of strongest contacts between residue pairs in the dense phase. These results provide a detailed picture of condensate formation by a key segment of the α-Syn molecule.

Application of modern artificial intelligence techniques in the development of organic molecular force fields

The molecular force field (FF) determines the accuracy of molecular dynamics (MD) and is one of the major bottlenecks that limits the application of MD in molecular design. Recently, artificial intelligence (AI) techniques, such as machine-learning potentials (MLPs), have been rapidly reshaping the landscape of MD. Meanwhile, organic molecular systems feature unique characteristics, and require more careful treatment in both model construction, optimization, and validation. While an accurate and generic organic molecular force field is still missing, significant progress has been made with the facilitation of AI, warranting a promising future. In this review, we provide an overview of the various types of AI techniques used in molecular FF development and discuss both the advantages and weaknesses of these methodologies. We show how AI methods provide unprecedented capabilities in many tasks such as potential fitting, atom typification, and automatic optimization. Meanwhile, it is also worth noting that more efforts are needed to improve the transferability of the model, develop a more comprehensive database, and establish more standardized validation procedures. With these discussions, we hope to inspire more efforts to solve the existing problems, eventually leading to the birth of next-generation generic organic FFs.

Emerging regulatory mechanisms and functions of biomolecular condensates: implications for therapeutic targets

Cells orchestrate their processes through complex interactions, precisely organizing biomolecules in space and time. Recent discoveries have highlighted the crucial role of biomolecular condensates-membrane-less assemblies formed through the condensation of proteins, nucleic acids, and other molecules-in driving efficient and dynamic cellular processes. These condensates are integral to various physiological functions, such as gene expression and intracellular signal transduction, enabling rapid and finely tuned cellular responses. Their ability to regulate cellular signaling pathways is particularly significant, as it requires a careful balance between flexibility and precision. Disruption of this balance can lead to pathological conditions, including neurodegenerative diseases, cancer, and viral infections. Consequently, biomolecular condensates have emerged as promising therapeutic targets, with the potential to offer novel approaches to disease treatment. In this review, we present the recent insights into the regulatory mechanisms by which biomolecular condensates influence intracellular signaling pathways, their roles in health and disease, and potential strategies for modulating condensate dynamics as a therapeutic approach. Understanding these emerging principles may provide valuable directions for developing effective treatments targeting the aberrant behavior of biomolecular condensates in various diseases.

Recent Progress in Modeling and Simulation of Biomolecular Crowding and Condensation Inside Cells

Macromolecular crowding in the cellular cytoplasm can potentially impact diffusion rates of proteins, their intrinsic structural stability, binding of proteins to their corresponding partners as well as biomolecular organization and phase separation. While such intracellular crowding can have a large impact on biomolecular structure and function, the molecular mechanisms and driving forces that determine the effect of crowding on dynamics and conformations of macromolecules are so far not well understood. At a molecular level, computational methods can provide a unique lens to investigate the effect of macromolecular crowding on biomolecular behavior, providing us with a resolution that is challenging to reach with experimental techniques alone. In this review, we focus on the various physics-based and data-driven computational methods developed in the past few years to investigate macromolecular crowding and intracellular protein condensation. We review recent progress in modeling and simulation of biomolecular systems of varying sizes, ranging from single protein molecules to the entire cellular cytoplasm. We further discuss the effects of macromolecular crowding on different phenomena, such as diffusion, protein-ligand binding, and mechanical and viscoelastic properties, such as surface tension of condensates. Finally, we discuss some of the outstanding challenges that we anticipate the community addressing in the next few years in order to investigate biological phenomena in model cellular environments by reproducing in vivo conditions as accurately as possible.

Nonspecific Yet Selective Interactions Contribute to Small Molecule Condensate Binding

Biomolecular condensates are essential in various cellular processes, and their misregulation has been demonstrated to underlie disease. Small molecules that modulate condensate stability and material properties offer promising therapeutic approaches, but mechanistic insights into their interactions with condensates remain largely lacking. We employ a multiscale approach to enable long-time, equilibrated all-atom simulations of various condensate-ligand systems. Systematic characterization of the ligand binding poses reveals that condensates can form diverse and heterogeneous chemical environments with one or multiple chains to bind small molecules. Unlike traditional protein-ligand interactions, these chemical environments are dominated by nonspecific hydrophobic interactions. Nevertheless, the chemical environments feature unique amino acid compositions and physicochemical properties that favor certain small molecules over others, resulting in varied ligand partitioning coefficients within condensates. Notably, different condensates share similar sets of chemical environments but at different populations. This population shift drives ligand selectivity toward specific condensates. Our approach can enhance the interpretation of experimental screening data and may assist in the rational design of small molecules targeting specific condensates.

Non-specific yet selective interactions contribute to small molecule condensate partitioning behavior

Biomolecular condensates are essential in various cellular processes, and their misregulation has been demonstrated to be underly disease. Small molecules that modulate condensate stability and material properties offer promising therapeutic approaches, but mechanistic insights into their interactions with condensates remain largely lacking. We employ a multiscale approach to enable long-time, equilibrated all-atom simulations of various condensate-ligand systems. Systematic characterization of the ligand binding poses reveals that condensates can form diverse and heterogeneous chemical environments with one or multiple chains to bind small molecules. Unlike traditional protein-ligand interactions, these chemical environments are dominated by non-specific hydrophobic interactions. Nevertheless, the chemical environments feature unique amino acid compositions and physicochemical properties that favor certain small molecules over others, resulting in varied ligand partitioning coefficients within condensates. Notably, different condensates share similar sets of chemical environments but at different populations. This population shift drives ligand selectivity towards specific condensates. Our approach can enhance the interpretation of experimental screening data and may assist in the rational design of small molecules targeting specific condensates.

Computational Chemistry in Structure-Based Solute Carrier Transporter Drug Design: Recent Advances and Future Perspectives

Solute carrier transporters (SLCs) are a class of important transmembrane proteins that are involved in the transportation of diverse solute ions and small molecules into cells. There are approximately 450 SLCs within the human body, and more than a quarter of them are emerging as attractive therapeutic targets for multiple complex diseases, e.g., depression, cancer, and diabetes. However, only 44 unique transporters (∼9.8% of the SLC superfamily) with 3D structures and specific binding sites have been reported. To design innovative and effective drugs targeting diverse SLCs, there are a number of obstacles that need to be overcome. However, computational chemistry, including physics-based molecular modeling and machine learning- and deep learning-based artificial intelligence (AI), provides an alternative and complementary way to the classical drug discovery approach. Here, we present a comprehensive overview on recent advances and existing challenges of the computational techniques in structure-based drug design of SLCs from three main aspects: (i) characterizing multiple conformations of the proteins during the functional process of transportation, (ii) identifying druggability sites especially the cryptic allosteric ones on the transporters for substrates and drugs binding, and (iii) discovering diverse small molecules or synthetic protein binders targeting the binding sites. This work is expected to provide guidelines for a deep understanding of the structure and function of the SLC superfamily to facilitate rational design of novel modulators of the transporters with the aid of state-of-the-art computational chemistry technologies including artificial intelligence.

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.

Martini on the Rocks: Can a Coarse-Grained Force Field Model Crystals?

Computational chemistry is an important tool in numerous scientific disciplines, including drug discovery and structural biology. Coarse-grained models offer simple representations of molecular systems that enable simulations of large-scale systems. Because there has been an increase in the adoption of such models for simulations of biomolecular systems, critical evaluation is warranted. Here, the stability of the amyloid peptide and organic crystals is evaluated using the Martini 3 coarse-grained force field. The crystals change shape drastically during the simulations. Radial distribution functions show that the distance between backbone beads in β-sheets increases by ∼1 Å, breaking the crystals. The melting points of organic compounds are much too low in the Martini force field. This suggests that Martini 3 lacks the specific interactions needed to accurately simulate peptides or organic crystals without imposing artificial restraints. The problems may be exacerbated by the use of the 12-6 potential, suggesting that a softer potential could improve this model for crystal simulations.

AMBERff at Scale: Multimillion-Atom Simulations with AMBER Force Fields in NAMD

All-atom molecular dynamics (MD) simulations are an essential structural biology technique with increasing application to multimillion-atom systems, including viruses and cellular machinery. Classical MD simulations rely on parameter sets, such as the AMBER family of force fields (AMBERff), to accurately describe molecular motion. Here, we present an implementation of AMBERff for use in NAMD that overcomes previous limitations to enable high-performance, massively parallel simulations encompassing up to two billion atoms. Single-point potential energy comparisons and case studies on model systems demonstrate that the implementation produces results that are as accurate as running AMBERff in its native engine.

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

Intrinsically disordered proteins (IDPs) frequently mediate phase separation that underlies the formation of a biomolecular condensate. Together with theory and experiment, efficient coarse-grained (CG) simulations have been instrumental in understanding the sequence-specific phase separation of IDPs. However, the widely used Cα-only models are limited in capturing the peptide nature of IDPs, particularly backbone-mediated interactions and effects of secondary structures, in phase separation. Here, we describe a hybrid resolution (HyRes) protein model toward a more accurate description of the backbone and transient secondary structures in phase separation. With an atomistic backbone and coarse-grained side chains, HyRes can semiquantitatively capture the residue helical propensity and overall chain dimension of monomeric IDPs. Using GY-23 as a model system, we show that HyRes is efficient enough for the direct simulation of spontaneous phase separation and, at the same time, appears accurate enough to resolve the effects of single His to Lys mutations. HyRes simulations also successfully predict increased β-structure formation in the condensate, consistent with available experimental CD data. We further utilize HyRes to study the phase separation of TPD-43, where several disease-related mutants in the conserved region (CR) have been shown to affect residual helicities and modulate the phase separation propensity as measured by the saturation concentration. The simulations successfully recapitulate the effect of these mutants on the helicity and phase separation propensity of TDP-43 CR. Analyses reveal that the balance between backbone and side chain-mediated interactions, but not helicity itself, actually determines phase separation propensity. These results support that HyRes represents an effective protein model for molecular simulation of IDP phase separation and will help to elucidate the coupling between transient secondary structures and phase separation.

RNA regulates cohesiveness and porosity of a biological condensate

Biological condensates have emerged as key elements of a biological cell function, concentrating disparate biomolecules to accomplish specific biological tasks. RNA was identified as a key ingredient of such condensates, however, its effect on the physical properties of the condensate was found to depend on the condensate's composition while its effect on the microstructure has remained elusive. Here, we characterize the physical properties and the microstructure of a protein-RNA condensate by means of large-scale coarse-grained (CG) molecular dynamics simulations. By developing a custom CG model of RNA compatible with a popular CG model of proteins, we systematically investigate the structural, thermodynamic, and kinetic properties of condensate droplets containing thousands of individual protein and RNA molecules over a range of temperatures. While we find RNA to increase the condensate's cohesiveness, its effect on the condensate's fluidity is more nuanced with longer molecules compacting the condensate and making it less fluid. We show that a biological condensate has a sponge-like morphology of interconnected channels of size that increases with temperature and decreases in the presence of RNA. Our results suggest that longer RNA form a dynamic scaffold within a condensate, regulating not only its fluidity but also permeability to intruder molecules.

Developing an Implicit Solvation Machine Learning Model for Molecular Simulations of Ionic Media

Molecular dynamics (MD) simulations of biophysical systems require accurate modeling of their native environment, i.e., aqueous ionic solution, as it critically impacts the structure and function of biomolecules. On the other hand, the models should be computationally efficient to enable simulations of large spatiotemporal scales. Here, we present the deep implicit solvation model for sodium chloride solutions that satisfies both requirements. Owing to the use of the neural network potential, the model can capture the many-body potential of mean force, while the implicit water treatment renders the model inexpensive. We demonstrate our approach first for pure ionic solutions with concentrations ranging from physiological to 2 M. We then extend the model to capture the effective ion interactions in the vicinity and far away from a DNA molecule. In both cases, the structural properties are in good agreement with all-atom MD, showcasing a general methodology for the efficient and accurate modeling of ionic media.

"Boundary residues" between the folded RNA recognition motif and disordered RGG domains are critical for FUS-RNA binding

Fused in sarcoma (FUS) is an abundant RNA-binding protein, which drives phase separation of cellular condensates and plays multiple roles in RNA regulation. The RNA-binding ability of FUS protein is crucial to its cellular function. Here, our molecular simulation study on the FUS-RNA complex provides atomic resolution insights into the observations from biochemical studies and also illuminates our understanding of molecular driving forces that mediate the structure, stability, and interaction of the RNA recognition motif (RRM) and RGG domains of FUS with a stem-loop junction RNA. We observe clear cooperativity and division of labor among the ordered (RRM) and disordered domains (RGG1 and RGG2) of FUS that leads to an organized and tighter RNA binding. Irrespective of the length of RGG2, the RGG2-RNA interaction is confined to the stem-loop junction and the proximal stem regions. On the other hand, the RGG1 interactions are primarily with the longer RNA stem. We find that the C terminus of RRM, which make up the "boundary residues" that connect the folded RRM with the long disordered RGG2 stretch of the protein, plays a critical role in FUS-RNA binding. Our study provides high-resolution molecular insights into the FUS-RNA interactions and forms the basis for understanding the molecular origins of full-length FUS interaction with RNA.

In-Cell Dynamics: The Next Focus of All-Atom Simulations

The cell is a crowded space where large biomolecules and metabolites are in continuous motion. Great strides have been made in in vitro studies of protein dynamics, folding, and protein-protein interactions, and much new data are emerging of how they differ in the cell. In this Perspective, we highlight the current progress in atomistic modeling of in-cell environments, both bacteria and mammals, with emphasis on classical all-atom molecular dynamics simulations. These simulations have been recently used to capture and characterize functional and non-functional protein-protein interactions, protein folding dynamics of small proteins with varied topologies, and dynamics of metabolites. We further discuss the challenges and efforts for updating modern force fields critical to the progress of cellular environment simulations. We also briefly summarize developments in relevant state-of-the-art experimental techniques. As computational and experimental methodologies continue to progress and produce more directly comparable data, we are poised to capture the complex atomistic picture of the cell.

Multiscale simulations reveal TDP-43 molecular-level interactions driving condensation

The RNA-binding protein TDP-43 is associated with mRNA processing and transport from the nucleus to the cytoplasm. TDP-43 localizes in the nucleus as well as accumulating in cytoplasmic condensates such as stress granules. Aggregation and formation of amyloid-like fibrils of cytoplasmic TDP-43 are hallmarks of numerous neurodegenerative diseases, most strikingly present in >90% of amyotrophic lateral sclerosis (ALS) patients. If excessive accumulation of cytoplasmic TDP-43 causes, or is caused by, neurodegeneration is presently not known. In this work, we use molecular dynamics simulations at multiple resolutions to explore TDP-43 self- and cross-interaction dynamics. A full-length molecular model of TDP-43, all 414 amino acids, was constructed from select structures of the protein functional domains (N-terminal domain, and two RNA recognition motifs, RRM1 and RRM2) and modeling of disordered connecting loops and the low complexity glycine-rich C-terminus domain. All-atom CHARMM36m simulations of single TDP-43 proteins served as guides to construct a coarse-grained Martini 3 model of TDP-43. The Martini model and a coarser implicit solvent C⍺ model, optimized for disordered proteins, were subsequently used to probe TDP-43 interactions; self-interactions from single-chain full-length TDP-43 simulations, cross-interactions from simulations with two proteins and simulations with assemblies of dozens to hundreds of proteins. Our findings illustrate the utility of different modeling scales for accessing TDP-43 molecular-level interactions and suggest that TDP-43 has numerous interaction preferences or patterns, exhibiting an overall strong, but dynamic, association and driving the formation of biomolecular condensates.

Viscosity Prediction of High-Concentration Antibody Solutions with Atomistic Simulations

The computational prediction of the viscosity of dense protein solutions is highly desirable, for example, in the early development phase of high-concentration biopharmaceutical formulations where the material needed for experimental determination is typically limited. Here, we use large-scale atomistic molecular dynamics (MD) simulations with explicit solvation to de novo predict the dynamic viscosities of solutions of a monoclonal IgG1 antibody (mAb) from the pressure fluctuations using a Green-Kubo approach. The viscosities at simulated mAb concentrations of 200 and 250 mg/mL are compared to the experimental values, which we measured with rotational rheometry. The computational viscosity of 24 mPa·s at the mAb concentration of 250 mg/mL matches the experimental value of 23 mPa·s obtained at a concentration of 213 mg/mL, indicating slightly different effective concentrations (or activities) in the MD simulations and in the experiments. This difference is assigned to a slight underestimation of the effective mAb-mAb interactions in the simulations, leading to a too loose dynamic mAb network that governs the viscosity. Taken together, this study demonstrates the feasibility of all-atom MD simulations for predicting the properties of dense mAb solutions and provides detailed microscopic insights into the underlying molecular interactions. At the same time, it also shows that there is room for further improvements and highlights challenges, such as the massive sampling required for computing collective properties of dense biomolecular solutions in the high-viscosity regime with reasonable statistical precision.

Thermodynamic forces from protein and water govern condensate formation of an intrinsically disordered protein domain

Liquid-liquid phase separation (LLPS) can drive a multitude of cellular processes by compartmentalizing biological cells via the formation of dense liquid biomolecular condensates, which can function as membraneless organelles. Despite its importance, the molecular-level understanding of the underlying thermodynamics of this process remains incomplete. In this study, we use atomistic molecular dynamics simulations of the low complexity domain (LCD) of human fused in sarcoma (FUS) protein to investigate the contributions of water and protein molecules to the free energy changes that govern LLPS. Both protein and water components are found to have comparably sizeable thermodynamic contributions to the formation of FUS condensates. Moreover, we quantify the counteracting effects of water molecules that are released into the bulk upon condensate formation and the waters retained within the protein droplets. Among the various factors considered, solvation entropy and protein interaction enthalpy are identified as the most important contributions, while solvation enthalpy and protein entropy changes are smaller. These results provide detailed molecular insights on the intricate thermodynamic interplay between protein- and solvation-related forces underlying the formation of biomolecular condensates.

Engineering Biological Nanopore Approaches toward Protein Sequencing

Biotechnological innovations have vastly improved the capacity to perform large-scale protein studies, while the methods we have for identifying and quantifying individual proteins are still inadequate to perform protein sequencing at the single-molecule level. Nanopore-inspired systems devoted to understanding how single molecules behave have been extensively developed for applications in genome sequencing. These nanopore systems are emerging as prominent tools for protein identification, detection, and analysis, suggesting realistic prospects for novel protein sequencing. This review summarizes recent advances in biological nanopore sensors toward protein sequencing, from the identification of individual amino acids to the controlled translocation of peptides and proteins, with attention focused on device and algorithm development and the delineation of molecular mechanisms with the aid of simulations. Specifically, the review aims to offer recommendations for the advancement of nanopore-based protein sequencing from an engineering perspective, highlighting the need for collaborative efforts across multiple disciplines. These efforts should include chemical conjugation, protein engineering, molecular simulation, machine-learning-assisted identification, and electronic device fabrication to enable practical implementation in real-world scenarios.

Accurate Simulation of Coupling between Protein Secondary Structure and Liquid-Liquid Phase Separation

Intrinsically disordered proteins (IDPs) frequently mediate liquid-liquid phase separation (LLPS) that underlies the formation of membraneless organelles. Together with theory and experiment, efficient coarse-grained (CG) simulations have been instrumental in understanding sequence-specific phase separation of IDPs. However, the widely-used Cα-only models are severely limited in capturing the peptide nature of IDPs, including backbone-mediated interactions and effects of secondary structures, in LLPS. Here, we describe a hybrid resolution (HyRes) protein model for accurate description of the backbone and transient secondary structures in LLPS. With an atomistic backbone and coarse-grained side chains, HyRes accurately predicts the residue helical propensity and chain dimension of monomeric IDPs. Using GY-23 as a model system, we show that HyRes is efficient enough for direct simulation of spontaneous phase separation, and at the same time accurate enough to resolve the effects of single mutations. HyRes simulations also successfully predict increased beta-sheet formation in the condensate, consistent with available experimental data. We further utilize HyRes to study the phase separation of TPD-43, where several disease-related mutants in the conserved region (CR) have been shown to affect residual helicities and modulate LLPS propensity. The simulations successfully recapitulate the effect of these mutants on the helicity and LLPS propensity of TDP-43 CR. Analyses reveal that the balance between backbone and sidechain-mediated interactions, but not helicity itself, actually determines LLPS propensity. We believe that the HyRes model represents an important advance in the molecular simulation of LLPS and will help elucidate the coupling between IDP transient secondary structures and phase separation.