Optimizing the Martini 3 Force Field Reveals the Effects of the Intricate Balance between Protein-Water Interaction Strength and Salt Concentration on Biomolecular Condensate Formation.

Condensation/dissolution has become a widely acknowledged biological macromolecular assembly phenomenon in subcellular compartmentalization. The MARTINI force field offers a coarse-grained protein model with a resolution that preserves molecular details with an explicit (CG) solvent. Despite its relatively higher resolution, it can still achieve condensate formation in a reasonable computing time with explicit solvent and ionic species. Therefore, it is highly desirable to tune this force field to be able to reproduce the experimentally observed properties of the condensate formation. In this work, we studied the condensate formation of the low-sequence complexity domain of fused in sarcoma protein using a MARTINI 3 force field by systematically modifying (increasing) the protein-water interaction strength and varying the salt concentration. We found that the condensate formation is sensitive both to the protein-water interaction strength and the presence of salt. While the unmodified MARTINI force field yields a complete collapse of proteins into one dense phase (i.e., no dilute phase), we reported a range of modified protein-water interaction strength that is capable of capturing the experimentally found transfer free energy between dense and dilute phases. We also found that the condensates lose their spherical shape upon the addition of salt, especially when the protein-water interactions are weak. Interchain amino acid contact map analysis showed one explanation for this observation: the protein-protein contact fraction reduces as salt is added to systems (when the protein-water interactions are weak), consistent with electrostatic screening effects. This reduction might be responsible for the condensates becoming nonspherical upon the addition of salt by reducing the need for minimizing the interfacial area. However, as the protein-water interactions become stronger to the extent that makes the transfer free energy agree well with experimentally observed transfer free energy, we found an increase in the protein-protein contact fraction upon the addition of salt, consistent with the salting-out effects. Therefore, we concluded that there is an intricate balance between screening effects and salting-out effects upon the addition of salt and this balance is highly sensitive to the strength of protein-water interactions.

Comparison of On-the-Fly Probability Enhanced Sampling and Parallel Tempering Combined with Metadynamics for Atomistic Simulations of RNA Tetraloop Folding.

Atomistic simulations with reliable models are extremely useful in providing exquisitely detailed pictures of biomolecular phenomena that are not always accessible to experiments. One such biomolecular phenomenon is RNA folding, which often requires exhaustive simulations with combined advanced sampling techniques. In this work, we employed the multithermal-multiumbrella on-the-fly probability enhanced sampling (MM-OPES) technique and compared it against combined parallel tempering and metadynamics simulations. We found that MM-OPES simulations were successful in reproducing the free energy surfaces from combined parallel tempering and metadynamics simulations. Importantly, we also investigated a wide range of temperature sets (minimum and maximum temperatures) for MM-OPES simulations in order to identify some guidelines for deciding the temperature limits for an accurate and efficient exploration of the free energy landscapes. We found that most temperature sets yielded almost the same accuracy in reproducing the free energy surface at the ambient conditions as long as (i) the maximum temperature is reasonably high, (ii) the temperature at which we run the simulation is reasonably high (in our simulations, running temperature is defined as [minimum temperature + maximum temperature]/2), and (iii) the effective sample size at the temperature of interest is statistically reasonable. In terms of the computational cost, all MM-OPES simulations were nearly 4 times less costly than the combined parallel tempering and metadynamics simulations. We concluded that the demanding combined parallel tempering and metadynamics simulations can safely be replaced with approximately 4 times less costly MM-OPES simulations (with carefully selected temperature limits) to obtain the same information.

A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States.

An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. We isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.

Thermodynamics of DNA Hybridization from Atomistic Simulations.

Studying DNA hybridization equilibrium at atomistic length scales, either via molecular dynamics (MD) or through commonly used advanced sampling approaches, is notoriously difficult. In this work, we describe an order-parameter-based advanced sampling technique to calculate the free energy of hybridization, and estimate the melting temperature of DNA oligomers at atomistic resolution. The free energy landscapes are reported as a function of a native-topology-based order parameter for the Drew-Dickerson dodecamer and for a range of DNA decamer sequences of different GC content. Our estimated melting temperatures match the experimental numbers within ±15 °C. As a test of the numerical reliability of the procedures employed, it was verified that the predicted free energy surfaces and melting temperatures of the d- and l-enantiomers of the Drew-Dickerson dodecamer were indistinguishable within numerical accuracy.

Second critical point in two realistic models of water.

The hypothesis that water has a second critical point at deeply supercooled conditions was formulated to provide a thermodynamically consistent interpretation of numerous experimental observations. A large body of work has been devoted to verifying or falsifying this hypothesis, but no unambiguous experimental proof has yet been found. Here, we use histogram reweighting and large-system scattering calculations to investigate computationally two molecular models of water, TIP4P/2005 and TIP4P/Ice, widely regarded to be among the most accurate classical force fields for this substance. We show that both models have a metastable liquid-liquid critical point at deeply supercooled conditions and that this critical point is consistent with the three-dimensional Ising universality class.

The Handedness of DNA Assembly around Carbon Nanotubes Is Determined by the Chirality of DNA.

Single-stranded DNA chains enable postsynthesis sorting of single-walled carbon nanotubes (CNTs) according to their diameter and helicity by wrapping helically around CNT surfaces. Both DNA chains and CNTs in these CNT-DNA conjugates are intrinsically chiral. Using a single-stranded DNA chain in both of its chiral realizations, we systematically study cross-chiral interactions between DNA and CNTs by varying the helicity of CNTs within a relatively narrow range of diameters. We find that regardless of the helicity or handedness of the carbon nanotube, the chirality of DNA dictates the handedness of its predominant helical wrap around carbon nanotubes.

TDP-43 α-helical structure tunes liquid-liquid phase separation and function.

Liquid-liquid phase separation (LLPS) is involved in the formation of membraneless organelles (MLOs) associated with RNA processing. The RNA-binding protein TDP-43 is present in several MLOs, undergoes LLPS, and has been linked to the pathogenesis of amyotrophic lateral sclerosis (ALS). While some ALS-associated mutations in TDP-43 disrupt self-interaction and function, here we show that designed single mutations can enhance TDP-43 assembly and function via modulating helical structure. Using molecular simulation and NMR spectroscopy, we observe large structural changes upon dimerization of TDP-43. Two conserved glycine residues (G335 and G338) are potent inhibitors of helical extension and helix-helix interaction, which are removed in part by variants at these positions, including the ALS-associated G335D. Substitution to helix-enhancing alanine at either of these positions dramatically enhances phase separation in vitro and decreases fluidity of phase-separated TDP-43 reporter compartments in cells. Furthermore, G335A increases TDP-43 splicing function in a minigene assay. Therefore, the TDP-43 helical region serves as a short but uniquely tunable module where application of biophysical principles can precisely control assembly and function in cellular and synthetic biology applications of LLPS.

Computational investigation of retro-isomer equilibrium structures: Intrinsically disordered, foldable, and cyclic peptides.

We use all-atom modeling and advanced-sampling molecular dynamics simulations to investigate quantitatively the effect of peptide bond directionality on the equilibrium structures of four linear (two foldable, two disordered) and two cyclic peptides. We find that the retro forms of cyclic and foldable linear peptides adopt distinctively different conformations compared to their parents. While the retro form of a linear intrinsically disordered peptide with transient secondary structure fails to reproduce a secondary structure content similar to that of its parent, the retro form of a shorter disordered linear peptide shows only minor differences compared to its parent.

Molecular interactions underlying liquid-liquid phase separation of the FUS low-complexity domain.

The low-complexity domain of the RNA-binding protein FUS (FUS LC) mediates liquid-liquid phase separation (LLPS), but the interactions between the repetitive SYGQ-rich sequence of FUS LC that stabilize the liquid phase are not known in detail. By combining NMR and Raman spectroscopy, mutagenesis, and molecular simulation, we demonstrate that heterogeneous interactions involving all residue types underlie LLPS of human FUS LC. We find no evidence that FUS LC adopts conformations with traditional secondary structure elements in the condensed phase; rather, it maintains conformational heterogeneity. We show that hydrogen bonding, π/sp2, and hydrophobic interactions all contribute to stabilizing LLPS of FUS LC. In addition to contributions from tyrosine residues, we find that glutamine residues also participate in contacts leading to LLPS of FUS LC. These results support a model in which FUS LC forms dynamic, multivalent interactions via multiple residue types and remains disordered in the densely packed liquid phase.

Evolution of All-Atom Protein Force Fields to Improve Local and Global Properties.

Experimental studies on intrinsically disordered and unfolded proteins have shown that in isolation they typically have low populations of secondary structure and exhibit distance scalings suggesting that they are at near-theta-solvent conditions. Until recently, however, all-atom force fields failed to reproduce these fundamental properties of intrinsically disordered proteins (IDPs). Recent improvements by refining against ensemble-averaged experimental observables for polypeptides in aqueous solution have addressed deficiencies including secondary structure bias, global conformational properties, and thermodynamic parameters of biophysical reactions such as folding and collapse. To date, studies utilizing these improved all-atom force fields have mostly been limited to a small set of unfolded or disordered proteins. Here, we present data generated for a diverse library of unfolded or disordered proteins using three progressively improved generations of Amber03 force fields, and we explore how global and local properties are affected by each successive change in the force field. We find that the most recent force field refinements significantly improve the agreement of the global properties such as radii of gyration and end-to-end distances with experimental estimates. However, these global properties are largely independent of the local secondary structure propensity. This result stresses the need to validate force fields with reference to a combination of experimental data providing information about both local and global structure formation.

Effect of heterochiral inversions on the structure of a ß-hairpin peptide.

We study computationally a family of ß-hairpin peptides with systematically introduced chiral inversions, in explicit water, and we investigate the extent to which the backbone structure is able to fold in the presence of heterochiral perturbations. In contrast to the recently investigated case of a helical peptide, we do not find a monotonic change in secondary structure content as a function of the number of L- to D-inversions. The effects of L- to D-inversions are instead found to be highly position-specific. Additionally, in contrast to the helical peptide, some inversions increase the stability of the folded peptide: in such cases, we compute an increase in ß-sheet content in the aqueous solution equilibrium ensemble. However, the tertiary structures of the stable (folded) configurations for peptides for which inversions cause an increase in ß-sheet content show differences from one another, as well as from the native fold of the nonchirally perturbed ß-hairpin. Our results suggest that although some chiral perturbations can increase folding stability, chirally perturbed proteins may still underperform functionally, given the relationship between structure and function.

Computational Investigation of the Effect of Backbone Chiral Inversions on Polypeptide Structure.

Studying a set of helix-folding polyalanine peptides with systematically inserted chiral inversions in explicit water, we investigate quantitatively the effect of chiral perturbations on the structural ensembles of the peptides, thereby assessing the extent to which the backbone structure is able to fold in the presence of systematic heterochiral perturbations. Starting from the homochiral l-Ala20 peptide, we invert the backbone chiralities of Ala residues one by one along a specific perturbation pathway, until reaching the homochiral d-Ala20 peptide. Analysis of the helical contents of the simulated structural ensembles of the peptides shows that even a single inversion in the middle of the peptide completely breaks the helical structure in its vicinity and drastically reduces the helical content of the peptide. Further inversions in the middle of the peptide monotonically decrease the original helical content, that is, the right-handed helical content for l-Ala, and increase the helical content of the opposite chirality. Further analysis of the peptide ensembles using several size- and shape-related order parameters also indicate the drastic global changes in the peptide structure due to the local effects caused by the chiral inversions, such as formation of a reverse turn. However, the degree of the structural changes introduced by opposite chirality substitutions depends on the position of the inversion.

Computational modeling highlights the role of the disordered Formin Homology 1 domain in profilin-actin transfer.

Formins accelerate actin polymerization, assumed to occur through flexible Formin Homology 1 (FH1) domain-mediated transfer of profilin-actin to the barbed end. To study FH1 properties and address sequence effects, including varying length/distribution of profilin-binding proline-rich motifs, we performed all-atom simulations of a set of representative FH1 domains of formins: mouse mDia1 and mDia2, budding yeast Bni1 and Bnr1, and fission yeast Cdc12, For3, and Fus1. We find FH1 has flexible regions between high-propensity polyproline helix regions. A coarse-grained model retaining sequence specificity, assuming rigid polyproline segments, describes their size. Multiple bound profilins or profilin-actin complexes expand mDia1-FH1, which may be important in cells. Simulations of the barbed end bound to Bni1-FH1-FH2 dimer show that the leading FH1 can better transfer profilin or profilin-actin, with decreasing probability as the distance from FH2 increases.

Inferring properties of disordered chains from FRET transfer efficiencies.

Förster resonance energy transfer (FRET) is a powerful tool for elucidating both structural and dynamic properties of unfolded or disordered biomolecules, especially in single-molecule experiments. However, the key observables, namely, the mean transfer efficiency and fluorescence lifetimes of the donor and acceptor chromophores, are averaged over a broad distribution of donor-acceptor distances. The inferred average properties of the ensemble therefore depend on the form of the model distribution chosen to describe the distance, as has been widely recognized. In addition, while the distribution for one type of polymer model may be appropriate for a chain under a given set of physico-chemical conditions, it may not be suitable for the same chain in a different environment so that even an apparently consistent application of the same model over all conditions may distort the apparent changes in chain dimensions with variation of temperature or solution composition. Here, we present an alternative and straightforward approach to determining ensemble properties from FRET data, in which the polymer scaling exponent is allowed to vary with solution conditions. In its simplest form, it requires either the mean FRET efficiency or fluorescence lifetime information. In order to test the accuracy of the method, we have utilized both synthetic FRET data from implicit and explicit solvent simulations for 30 different protein sequences, and experimental single-molecule FRET data for an intrinsically disordered and a denatured protein. In all cases, we find that the inferred radii of gyration are within 10% of the true values, thus providing higher accuracy than simpler polymer models. In addition, the scaling exponents obtained by our procedure are in good agreement with those determined directly from the molecular ensemble. Our approach can in principle be generalized to treating other ensemble-averaged functions of intramolecular distances from experimental data.

Highly Disordered Amyloid-ß Monomer Probed by Single-Molecule FRET and MD Simulation.

Monomers of amyloid-ß (Aß) protein are known to be disordered, but there is considerable controversy over the existence of residual or transient conformations that can potentially promote oligomerization and fibril formation. We employed single-molecule Förster resonance energy transfer (FRET) spectroscopy with site-specific dye labeling using an unnatural amino acid and molecular dynamics simulations to investigate conformations and dynamics of Aß isoforms with 40 (Aß40) and 42 residues (Aß42). The FRET efficiency distributions of both proteins measured in phosphate-buffered saline at room temperature show a single peak with very similar FRET efficiencies, indicating there is apparently only one state. 2D FRET efficiency-donor lifetime analysis reveals, however, that there is a broad distribution of rapidly interconverting conformations. Using nanosecond fluorescence correlation spectroscopy, we measured the timescale of the fluctuations between these conformations to be ∼35 ns, similar to that of disordered proteins. These results suggest that both Aß40 and Aß42 populate an ensemble of rapidly reconfiguring unfolded states, with no long-lived conformational state distinguishable from that of the disordered ensemble. To gain molecular-level insights into these observations, we performed molecular dynamics simulations with a force field optimized to describe disordered proteins. We find, as in experiments, that both peptides populate configurations consistent with random polymer chains, with the vast majority of conformations lacking significant secondary structure, giving rise to very similar ensemble-averaged FRET efficiencies.

Mechanistic View of hnRNPA2 Low-Complexity Domain Structure, Interactions, and Phase Separation Altered by Mutation and Arginine Methylation.

hnRNPA2, a component of RNA-processing membraneless organelles, forms inclusions when mutated in a syndrome characterized by the degeneration of neurons (bearing features of amyotrophic lateral sclerosis [ALS] and frontotemporal dementia), muscle, and bone. Here we provide a unified structural view of hnRNPA2 self-assembly, aggregation, and interaction and the distinct effects of small chemical changes-disease mutations and arginine methylation-on these assemblies. The hnRNPA2 low-complexity (LC) domain is compact and intrinsically disordered as a monomer, retaining predominant disorder in a liquid-liquid phase-separated form. Disease mutations D290V and P298L induce aggregation by enhancing and extending, respectively, the aggregation-prone region. Co-aggregating in disease inclusions, hnRNPA2 LC directly interacts with and induces phase separation of TDP-43. Conversely, arginine methylation reduces hnRNPA2 phase separation, disrupting arginine-mediated contacts. These results highlight the mechanistic role of specific LC domain interactions and modifications conserved across many hnRNP family members but altered by aggregation-causing pathological mutations.

Interplay Between Membrane Composition and Structural Stability of Membrane-Bound hIAPP.

Amyloid aggregates are characteristic of many serious diseases such as Alzheimer's disease, Parkinson's, and type 2 diabetes and commonly involve intrinsically disordered proteins (IDPs), those that populate an ensemble of conformations rather than a single folded structure. Human islet amyloid polypeptide (hIAPP or amylin) is an amyloidogenic IDP implicated in pancreatic ß-cell death during the pathogenesis of type 2 diabetes. The target of amylin's toxic activity is thought to be the cell's lipid membrane, which may also act as a catalyst for aggregation. Since amylin is intrinsically disordered, differing environments can have a large impact on its equilibrium conformational ensemble. We apply atomistic molecular dynamics simulations on multiple systems containing a full-length amylin monomer and a lipid bilayer to study the changes induced by the membrane. We observe stabilized helical conformations structurally similar to those determined by NMR experiments conducted in similar environments. We also find that bilayers of different compositions result in greatly different equilibrium ensembles of amylin. Finally, we discuss how a mixed bilayer containing zwitterionic and anionic lipid headgroups can allow for greater preference toward conformations which are adsorbed below the membrane surface through rearrangement of lipids for more favorable protein-lipid interactions.

A Carbon Nanotube Reporter of miRNA Hybridization Events In Vivo.

MicroRNAs and other small oligonucleotides in biofluids are promising disease biomarkers, yet conventional assays require complex processing steps that are unsuitable for point-of-care testing or for implantable or wearable sensors. Single-walled carbon nanotubes are an ideal material for implantable sensors, owing to their emission in the near-infrared spectral region, photostability and exquisite sensitivity. Here, we report an engineered carbon-nanotube-based sensor capable of real-time optical quantification of hybridization events of microRNA and other oligonucleotides. The mechanism of the sensor arises from competitive effects between displacement of both oligonucleotide charge groups and water from the nanotube surface, which result in a solvatochromism-like response. The sensor, which allows for detection via single-molecule sensor elements and for multiplexing by using multiple nanotube chiralities, can monitor toehold-based strand-displacement events, which reverse the sensor response and regenerate the sensor complex. We also show that the sensor functions in whole urine and serum, and can non-invasively measure DNA and microRNA after implantation in live mice.

Phosphorylation of the FUS low-complexity domain disrupts phase separation, aggregation, and toxicity.

Neuronal inclusions of aggregated RNA-binding protein fused in sarcoma (FUS) are hallmarks of ALS and frontotemporal dementia subtypes. Intriguingly, FUS's nearly uncharged, aggregation-prone, yeast prion-like, low sequence-complexity domain (LC) is known to be targeted for phosphorylation. Here we map in vitro and in-cell phosphorylation sites across FUS LC We show that both phosphorylation and phosphomimetic variants reduce its aggregation-prone/prion-like character, disrupting FUS phase separation in the presence of RNA or salt and reducing FUS propensity to aggregate. Nuclear magnetic resonance spectroscopy demonstrates the intrinsically disordered structure of FUS LC is preserved after phosphorylation; however, transient domain collapse and self-interaction are reduced by phosphomimetics. Moreover, we show that phosphomimetic FUS reduces aggregation in human and yeast cell models, and can ameliorate FUS-associated cytotoxicity. Hence, post-translational modification may be a mechanism by which cells control physiological assembly and prevent pathological protein aggregation, suggesting a potential treatment pathway amenable to pharmacologic modulation.