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.

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

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.

Human islet amyloid polypeptide (hIAPP), a 37-residue protein cosecreted with insulin by β-cells in the pancreas, is known to form amyloid fibrils in type II diabetes patients. This fibril formation is also associated with β-cell death. However, rat IAPP (rIAPP) does not aggregate into fibrils, nor is it associated with β-cell toxicity. Determining solution properties of hIAPP experimentally is difficult because it aggregates quickly. Our study uses molecular dynamics simulation to explore and compare in-solution characteristics of hIAPP and rIAPP, as well as two single-point hIAPP mutants, hIAPP I26P and hIAPP S20G, which exhibit observed differences from hIAPP in aggregation propensities. We find that all four polypeptide monomers sample structured states in solution. More importantly, differences in the helicity over residues 7-16 may play an important role in early aggregation, although this region is outside of commonly assumed amyloidogenic region 20-29. The long-range contacts, though unexpected of IDPs, cause variation in sampled conformations among four polypeptides within same amino acid sequence. Our results also yield evidence that previously determined structures bound to micelles are also transiently sampled in the solution state. In particular, similarities found in region 8-17 together with the helical differences that we observe in region 7-16 lead us to suggest that the region 7-16 is potentially responsible for amyloidogenic behavior of amylin peptides. Our results also provide support for the proposed mechanism of fibril formation based on experimentally observed transient helices in amyloidogenic peptides.