Krishna Neupane

Direct observation of multiple misfolding pathways in a single prion protein molecule

Protein misfolding is a ubiquitous phenomenon associated with a wide range of diseases. Single-molecule approaches offer a powerful tool for deciphering the mechanisms of misfolding by measuring the conformational fluctuations of a protein with high sensitivity. We applied single-molecule force spectroscopy to observe directly the misfolding of the prion protein PrP, a protein notable for having an infectious misfolded state that is able to propagate by recruiting natively folded PrP. By measuring folding trajectories of single PrP molecules held under tension in a high-resolution optical trap, we found that the native folding pathway involves only two states, without evidence for partially folded intermediates that have been proposed to mediate misfolding. Instead, frequent but fleeting transitions were observed into off-pathway intermediates. Three different misfolding pathways were detected, all starting from the unfolded state. Remarkably, the misfolding rate was even higher than the rate for native folding. A mutant PrP with higher aggregation propensity showed increased occupancy of some of the misfolded states, suggesting these states may act as intermediates during aggregation. These measurements of individual misfolding trajectories demonstrate the power of single-molecule approaches for characterizing misfolding directly by mapping out nonnative folding pathways.

Energy landscape analysis of native folding of the prion protein yields the diffusion constant, transition path time, and rates.

Protein folding is described conceptually in terms of diffusion over a configurational free-energy landscape, typically reduced to a one-dimensional profile along a reaction coordinate. In principle, kinetic properties can be predicted directly from the landscape profile using Kramers theory for diffusive barrier crossing, including the folding rates and the transition time for crossing the barrier. Landscape theory has been widely applied to interpret the time scales for protein conformational dynamics, but protein folding rates and transition times have not been calculated directly from experimentally measured free-energy profiles. We characterized the energy landscape for native folding of the prion protein using force spectroscopy, measuring the change in extension of a single protein molecule at high resolution as it unfolded/refolded under tension. Key parameters describing the landscape profile were first recovered from the distributions of unfolding and refolding forces, allowing the diffusion constant for barrier crossing and the transition path time across the barrier to be calculated. The full landscape profile was then reconstructed from force-extension curves, revealing a double-well potential with an extended, partially unfolded transition state. The barrier height and position were consistent with the previous results. Finally, Kramers theory was used to predict the folding rates from the landscape profile, recovering the values observed experimentally both under tension and at zero force in ensemble experiments. These results demonstrate how advances in single-molecule theory and experiment are harnessing the power of landscape formalisms to describe quantitatively the mechanics of folding.

Single-molecule force spectroscopy of the add adenine riboswitch relates folding to regulatory mechanism.

Riboswitches regulate gene expression via ligand binding to an aptamer domain which induces conformational changes in a regulatory expression platform. By unfolding and refolding single add adenine riboswitch molecules in an optical trap, an integrated picture of the folding was developed and related to the regulatory mechanism. Force-extension curves (FECs) and constant-force folding trajectories measured on the aptamer alone revealed multiple partially-folded states, including several misfolded states not on the native folding pathway. All states were correlated to key structural components and interactions within hierarchical folding pathways. FECs of the full-length riboswitch revealed that the thermodynamically stable conformation switches upon ligand binding from a structure repressing translation to one permitting it. Along with rapid equilibration of the two structures in the absence of adenine, these results support a thermodynamically-controlled regulatory mechanism, in contrast with the kinetic control of the closely-related pbuE adenine riboswitch. Comparison of the folding of these riboswitches revealed many similarities arising from shared structural features but also essential differences related to their different regulatory mechanisms.

Direct observation of transition paths during the folding of proteins and nucleic acids.

How biomolecules fold In order to fold, biomolecules must search a conformational energy landscape to find low-energy states. There are peaks in the landscape where the molecules must occupy unstable conformations for a short time. Neupane et al. used optical tweezers to observe these transition paths directly for single nucleic acid and protein molecules (see the Perspective by Wolynes). They measured a distribution of times taken to cross the transition path and found that the shape of the distribution agrees well with theory that assumes one-dimensional diffusion over the landscape. Science, this issue p. 239; see also p. 150 Optical tweezers reveal how a protein and a DNA hairpin cross the barrier between the folded and unfolded states. [Also see Perspective by Wolynes] Transition paths, the fleeting trajectories through the transition states that dominate the dynamics of biomolecular folding reactions, encapsulate the critical information about how structure forms. Owing to their brief duration, however, they have not previously been observed directly. We measured transition paths for both nucleic acid and protein folding, using optical tweezers to observe the microscopic diffusive motion of single molecules traversing energy barriers. The average transit times and the shapes of the transit-time distributions agreed well with theoretical expectations for motion over the one-dimensional energy landscapes reconstructed for the same molecules, validating the physical theory of folding reactions. These measurements provide a first look at the critical microscopic events that occur during folding, opening exciting new avenues for investigating folding phenomena.

Diverse metastable structures formed by small oligomers of α-synuclein probed by force spectroscopy.

Oligomeric aggregates are widely suspected as toxic agents in diseases caused by protein aggregation, yet they remain poorly characterized, partly because they are challenging to isolate from a heterogeneous mixture of species. We developed an assay for characterizing structure, stability, and kinetics of individual oligomers at high resolution and sensitivity using single-molecule force spectroscopy, and applied it to observe the formation of transient structured aggregates within single oligomers of α-synuclein, an intrinsically-disordered protein linked to Parkinson’s disease. Measurements of the molecular extension as the proteins unfolded under tension in optical tweezers revealed that even small oligomers could form numerous metastable structures, with a surprisingly broad range of sizes. Comparing the structures formed in monomers, dimers and tetramers, we found that the average mechanical stability increased with oligomer size. Most structures formed within a minute, with size-dependent rates. These results provide a new window onto the complex α-synuclein aggregation landscape, characterizing the microscopic structural heterogeneity and kinetics of different pathways.

Transition-Path Probability as a Test of Reaction-Coordinate Quality Reveals DNA Hairpin Folding Is a One-Dimensional Diffusive Process.

Chemical reactions are typically described in terms of progress along a reaction coordinate. However, the quality of reaction coordinates for describing reaction dynamics is seldom tested experimentally. We applied a framework for gauging reaction-coordinate quality based on transition-path analysis to experimental data for the first time, looking at folding trajectories of single DNA hairpin molecules measured under tension applied by optical tweezers. The conditional probability for being on a reactive transition path was compared with the probability expected for ideal diffusion over a 1D energy landscape based on the committor function. Analyzing measurements and simulations of hairpin folding where end-to-end extension is the reaction coordinate, after accounting for instrumental effects on the analysis, we found good agreement between transition-path and committor analyses for model two-state hairpins, demonstrating that folding is well-described by 1D diffusion. This work establishes transition-path analysis as a powerful new tool for testing experimental reaction-coordinate quality.

Quantifying Instrumental Artifacts in Folding Kinetics Measured by Single-Molecule Force Spectroscopy

Force spectroscopy is commonly used to measure the kinetics of processes occurring in single biological molecules. These measurements involve attaching the molecule of interest to micron-sized or larger force probes via compliant linkers. Recent theoretical work has described how the properties of the probes and linkers can alter the observed kinetics from the intrinsic behavior of the molecule in isolation. We applied this theory to estimate the errors in measurements of folding made using optical tweezers. Errors in the folding rates arising from instrument artifacts were only ∼20% for constant-force measurements of DNA hairpins with typical choices of linker length and probe size. Measurements of transition paths using a constant trap position at high trap stiffness were also found to be in the low-artifact limit. These results indicate that typical optical trap measurements of kinetics reflect the dynamics of the molecule fairly well, and suggest practical limitations on experimental design to ensure reliable kinetic measurements.

Pharmacological chaperone reshapes the energy landscape for folding and aggregation of the prion protein

The development of small-molecule pharmacological chaperones as therapeutics for protein misfolding diseases has proven challenging, partly because their mechanism of action remains unclear. Here we study Fe-TMPyP, a tetrapyrrole that binds to the prion protein PrP and inhibits misfolding, examining its effects on PrP folding at the single-molecule level with force spectroscopy. Single PrP molecules are unfolded with and without Fe-TMPyP present using optical tweezers. Ligand binding to the native structure increases the unfolding force significantly and alters the transition state for unfolding, making it more brittle and raising the barrier height. Fe-TMPyP also binds the unfolded state, delaying native refolding. Furthermore, Fe-TMPyP binding blocks the formation of a stable misfolded dimer by interfering with intermolecular interactions, acting in a similar manner to some molecular chaperones. The ligand thus promotes native folding by stabilizing the native state while also suppressing interactions driving aggregation.

Direct measurement of sequence-dependent transition path times and conformational diffusion in DNA duplex formation

Significance The transition path time in folding reactions—the time required to cross through the high-energy transition states—offers a powerful tool for studying folding mechanisms because it is uniquely sensitive to the microscopic motions along the transition paths that dominate the folding dynamics. We used transition path times to probe helix formation in nucleic acids and elucidate differences between the formation of canonical G:C and A:T base pairs. Stem-loop hairpins of the same size but increased G:C content had faster transition path times, revealing that the conformational diffusion coefficient determining the timescale for molecular dynamics was higher for formation of G:C base pairs than A:T base pairs. Such sequence dependence, not previously observed, reflects differences in the internal friction for reconfiguring the bases. The conformational diffusion coefficient, D, sets the timescale for microscopic structural changes during folding transitions in biomolecules like nucleic acids and proteins. D encodes significant information about the folding dynamics such as the roughness of the energy landscape governing the folding and the level of internal friction in the molecule, but it is challenging to measure. The most sensitive measure of D is the time required to cross the energy barrier that dominates folding kinetics, known as the transition path time. To investigate the sequence dependence of D in DNA duplex formation, we measured individual transition paths from equilibrium folding trajectories of single DNA hairpins held under tension in high-resolution optical tweezers. Studying hairpins with the same helix length but with G:C base-pair content varying from 0 to 100%, we determined both the average time to cross the transition paths, τtp, and the distribution of individual transit times, PTP(t). We then estimated D from both τtp and PTP(t) from theories assuming one-dimensional diffusive motion over a harmonic barrier. τtp decreased roughly linearly with the G:C content of the hairpin helix, being 50% longer for hairpins with only A:T base pairs than for those with only G:C base pairs. Conversely, D increased linearly with helix G:C content, roughly doubling as the G:C content increased from 0 to 100%. These results reveal that G:C base pairs form faster than A:T base pairs because of faster conformational diffusion, possibly reflecting lower torsional barriers, and demonstrate the power of transition path measurements for elucidating the microscopic determinants of folding.

Testing Kinetic Identities Involving Transition-Path Properties Using Single-Molecule Folding Trajectories

Recent advances in single-molecule assays have allowed individual transition paths during the folding of single molecules to be observed directly. We used the transition paths of DNA hairpins having different sequences, measured with high-resolution optical tweezers, to test theoretical relations between the properties of the transition paths and the folding kinetics. We showed that folding and unfolding rates were related to the average transition-path times, as expected from theory, for all hairpins studied. We also found that the probability distribution of transition-path occupancies agreed with the profile of the average velocity along the transition paths for each of the hairpins, as expected theoretically. Finally, we used the latter result to show that the committor probability recovered from the velocity profile matches the committor measured empirically. These results validate the proposed kinetic identities.