All-atom simulations with explicit solvent can provide a detailed and accurate description of dynamics and mechanisms in biomolecular systems, including intrinsically disordered proteins (IDPs) and their condensates. However, interactions involving charged residues and ions remain a persistent source of systematic error. Here we introduce an osmometry-guided optimization strategy that directly targets residue–residue, residue–ion and ion–ion interactions. Osmotic pressure provides key experimental information on molecular interactions and can be calculated directly and rapidly from simulations, enabling iterative force-field optimization. The resulting parameters improve agreement with single-molecule FRET data for IDPs, NMR relaxation data for an IDP–folded-domain complex, and chain dynamics and dimensions in biomolecular condensates of charged IDPs. For such condensates, simulations with our osmometry-optimized force field provide the missing link for predicting condensate dynamics across length and time scales. The strategy is broadly extensible to other interaction classes, including those governing protein–nucleic-acid assemblies.
An improved force field is developed, derived from its parent, Amber ff24EXP-GA, and its evaluation against Amber ff14SB and other contemporary force fields, such as CHARMM36m, in capturing the empirically determined conformational properties of unfolded systems: short peptides that serve as model systems for IDPs, and...
Condensates formed by oppositely charged intrinsically disordered proteins provide model systems for understanding how transient electrostatic interactions govern structure and dynamics in biomolecular assemblies. Here we investigate a nearly charge-neutral condensate composed of 50 Prothymosin alpha (ProTalpha) and 40...
Ionic liquids are salts that exist in the liquid state at room temperature and exhibit high viscosity because of their strong electrostatic interactions. It was difficult to reproduce their viscosity by molecular dynamics simulations with conventional nonpolarizable force fields; however, recent development of force fi...
This work investigates how protein folding landscapes are altered inside condensates, using the protein α-helix as a model folded domain and develops a chemically specific, residue-resolution model for quantification of α-helical folding and applies it to characterize diverse helices within condensates of varying physi...
Nathaniel Hess, Jerelle A. Joseph· Journal of the American Chem...· 1 citation
The results show that while the commonly used 10 Å buffer may be sufficient to maintain the stability of the double-stranded nucleic acid, larger cells are required to capture the conformational dynamics of single-stranded structures.
Nainsy Baghel, Pranchal Shrivastava, R. Mehra· bioRxiv· 0 citations
Aromatic organic solutes in water exhibit a delicate balance between hydrophobic solvation and directional O-H$\cdots \pi$ hydrogen bonds, yet widely used force fields and state-of-the-art density functional approaches struggle to provide a consistent picture of these pivotal interactions. We introduce a data-efficient...
N. Stolte, H. Forbert, Yu. Lysogorskiy et al.· 0 citations
A new machine-learning framework aims to improve the success rate of computational protein design while moving away from results that reproduce sequences found in nature.
