Search

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.

Condensation/dissolution has become a widely acknowledged biological macromolecular assembly phenomenon in subcellular compartmentalization. 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 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 (LC) domain of FUsed in Sarcoma (FUS) protein using a MARTINI 3 force field by systematically modifying (increasing) the protein-water interaction strength and varying 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. Inter-chain 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 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 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.

Studying the DNA hybridization equilibrium via brute force molecular dynamics (MD) or commonly used advanced sampling approches is notoriously difficult at atomistic lengthscale. However, besides providing a more realistic modeling of this microscopic phenomenon, atomistic resolution is a necessity for some fundamental research questions, such as the ones related to DNA’s chirality. Here, we describe an order parameter-based advanced sampling technique to calculate the free energy surface of hybridization and estimate melting temperature of DNA oligomers at atomistic resolution, using a native topology-based order parameter. We show that the melting temperatures estimated from our atomistic simulations follow an order consistent with the predictions from melting experiments and those from the nearest neighbor model, for a range of DNA sequences of different GC content. Moreover, free energy surfaces and melting temperatures are calculated to be identical for D- and L-enantiomers of Drew-Dickerson dodecamer.Abstract FigureGraphical TOC Entry

Liquid-liquid critical point of water Despite a broad range of experimental observations that indirectly point out the possible existence of a liquid-liquid phase transition in deeply supercooled water, no unambiguous experiment has shown this yet. The challenges of performing experiments under such conditions are associated with the inevitable rapid crystallization of the metastable liquid state, and therefore computational simulations become a crucial alternative. Using two accurate classical models of water and long (more than 40 microseconds) isobaricisothermal molecular dynamics simulations, Debenedetti et al. provide strong computational evidence of the presence of a second metastable liquid-liquid critical point in the deeply supercooled region that has critical behavior consistent with the three-dimensional Ising universality class. Science this issue p. 289

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.

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.

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.

Significance TDP-43 is an essential RNA-binding protein that assembles into protein inclusions in >95% of cases of amyotrophic lateral sclerosis (ALS). A partially helical region in the predominantly disordered C-terminal domain harbors several mutations associated with ALS and is important for TDP-43 function and liquid–liquid phase separation. We directly demonstrate that this helical region undergoes large structural changes upon helix–helix dimerization and that point mutations can enhance helix–helix assembly. Furthermore, we demonstrate that these point variants can be used to control the material properties of phase-separated TDP-43 constructs in cells and can enhance TDP-43 RNA-splicing function. Therefore, engineered forms of the TDP-43 helical domain could be used to control in-cell phase separation, dynamic assembly, 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.

Many nanomaterials are promising biosensor elements due to capabilities for transduction of biomolecular interactions into an electrical or optical signal. Certain nanomaterials have intrinsic charges and are thus susceptible to electrostatic forces that may enhance or attenuate their response. Here, using a combination of experimental and computational approaches, we found and characterized a critical role of the solvent salt conditions in determining the extent of optical changes due to anionic analyte interaction with an intrinsically responsive nanomaterial, single-walled carbon nanotubes. Using a well-characterized model anionic analyte, we found that monovalent salts enabled greater optical changes of a polyanion-bound carbon nanotube in a highly dose-dependent manner but not with a neutrally charged analyte. Molecular dynamics simulations were used to derive a quantitative understanding of this mechanism from a free-energy perspective. We also show that salt can be used to enhance the sensitivity o...

Neuronal inclusions of aggregated RNA‐binding protein fused in sarcoma (FUS) are hallmarks of ALS and frontotemporal dementia subtypes. Intriguingly, FUS9s 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.

Liquid-liquid phase separation (LLPS) is involved in the formation of membraneless organelles (MLOs) associated with RNA processing. Present in several MLOs, TDP-43 undergoes LLPS and is linked to the pathogenesis of amyotrophic lateral sclerosis (ALS). While some disease variants of 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 in a dimeric 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 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 mini-gene assay. Therefore, TDP-43 helical region serves as a short but uniquely tunable module that shows promise as for controlling assembly and function in cellular and synthetic biology applications of LLPS.