Hairpin Folding Research Articles

Molecular Compressive Force Sensor for Mapping Forces at the Cell-Substrate Interface.

Mechanical forces are crucial for biological processes such as T cell antigen recognition. A suite of molecular tension probes to measure pulling forces have been reported over the past decade; however, there are no reports of molecular probes for measuring compressive forces, representing a gap in the current mechanobiology toolbox. To address this gap, we report a molecular compression reporter using pseudostable hairpins (M-CRUSH). The design principle was based on a pseudostable DNA structure that folds in response to an external compressive force. We created a library of DNA stem-loop hairpins with varying thermodynamic stability, and then used Förster Resonance Energy Transfer (FRET) to quantify hairpin folding stability as a function of temperature and crowding. We identified an optimal pseudostable DNA hairpin highly sensitive to molecular crowding that displayed a shift in melting temperature (Tm) of 7 °C in response to a PEG crowding agent. When immobilized on surfaces, this optimized DNA hairpin showed a 29 ± 6% increase in FRET index in response to 25% w/w PEG 8K. As a proof-of-concept demonstration, we employed M-CRUSH to map the compressive forces generated by primary naïve T cells. We noted dynamic compressive forces that were highly sensitive to antigen presentation and coreceptor engagement. Importantly, mechanical forces are generated by cytoskeletal protrusions caused by acto-myosin activity. This was confirmed by treating cells with cytoskeletal inhibitors, which resulted in a lower FRET response when compared to untreated cells. Furthermore, we showed that M-CRUSH signal is dependent on probe density with greater density probes showing enhanced signal. Finally, we demonstrated that M-CRUSH probes are modular and can be applied to different cell types by displaying a compressive signal observed under human platelets. M-CRUSH offers a powerful tool to complement tension sensors and map out compressive forces in living systems.

Quantifying the Properties of Nonproductive Attempts at Thermally Activated Energy-Barrier Crossing through Direct Observation

Thermally activated energy-barrier crossing is ubiquitous in physical, chemical, and biological processes. Most barrier-crossing attempts have insufficient energy to overcome the barrier; hence, productive transition paths that successfully cross the barrier are very rare compared to nonproductive fluctuations that enter the barrier region but return without crossing it. Recent experimental advances have yielded important insights into transition paths, but nonproductive attempts remain little studied experimentally or theoretically, even though they can reveal information about parts of the reaction energy landscape not visited during transition paths. Observing the diffusive dynamics of a bead hopping between bistable optical traps as a model system, we measured the duration, maximum position along the reaction coordinate, and occupancy statistics of unsuccessful crossing attempts. Experimental results agreed quantitatively with expectations of an analytical framework we derived from committor theory. Applying these analyses to a more complex example, DNA hairpin folding under tension, we found that some properties differed from those of transition paths, such as the asymmetric occupancies for folding and unfolding attempts, whereas others were similar, such as the diffusion coefficient reflecting landscape roughness. These results show how nonproductive crossing attempts can be detected and analyzed rigorously, enabling characterization of the full dynamics within the transition region. Published by the American Physical Society 2024

Molecular Latent Space Simulators for Distributed and Multimolecular Trajectories.

All atom molecular dynamics (MD) simulations offer a powerful tool for molecular modeling, but the short time steps required for numerical stability of the integrator place many interesting molecular events out of reach of unbiased simulations. The popular and powerful Markov state modeling (MSM) approach can extend these time scales by stitching together multiple short discontinuous trajectories into a single long-time kinetic model but necessitates a configurational coarse-graining of the phase space that entails a loss of spatial and temporal resolution and an exponential increase in complexity for multimolecular systems. Latent space simulators (LSS) present an alternative formalism that employs a dynamical, as opposed to configurational, coarse graining comprising three back-to-back learning problems to (i) identify the molecular system's slowest dynamical processes, (ii) propagate the microscopic system dynamics within this slow subspace, and (iii) generatively reconstruct the trajectory of the system within the molecular phase space. A trained LSS model can generate temporally and spatially continuous synthetic molecular trajectories at orders of magnitude lower cost than MD to improve sampling of rare transition events and metastable states to reduce statistical uncertainties in thermodynamic and kinetic observables. In this work, we extend the LSS formalism to short discontinuous training trajectories generated by distributed computing and to multimolecular systems without incurring exponential scaling in computational cost. First, we develop a distributed LSS model over thousands of short simulations of a 264-residue proteolysis-targeting chimera (PROTAC) complex to generate ultralong continuous trajectories that identify metastable states and collective variables to inform PROTAC therapeutic design and optimization. Second, we develop a multimolecular LSS architecture to generate physically realistic ultralong trajectories of DNA oligomers that can undergo both duplex hybridization and hairpin folding. These trajectories retain thermodynamic and kinetic characteristics of the training data while providing increased precision of folding populations and time scales across simulation temperature and ion concentration.

CAG-Repeat RNA Hairpin Folding and Recruitment to Nuclear Speckles with a Pivotal Role of ATP as a Cosolute.

A hallmark of Huntington's disease (HD) is a prolonged polyglutamine sequence in the huntingtin protein and, correspondingly, an expanded cytosine, adenine, and guanine (CAG) triplet repeat region in the mRNA. A majority of studies investigating disease pathology were concerned with toxic huntingtin protein, but the mRNA moved into focus due to its recruitment to RNA foci and emerging novel therapeutic approaches targeting the mRNA. A hallmark of CAG-RNA is that it forms a stable hairpin in vitro which seems to be crucial for specific protein interactions. Using in-cell folding experiments, we show that the CAG-RNA is largely destabilized in cells compared to dilute buffer solutions but remains folded in the cytoplasm and nucleus. Surprisingly, we found the same folding stability in the nucleoplasm and in nuclear speckles under physiological conditions suggesting that CAG-RNA does not undergo a conformational transition upon recruitment to the nuclear speckles. We found that the metabolite adenosine triphosphate (ATP) plays a crucial role in promoting unfolding, enabling its recruitment to nuclear speckles and preserving its mobility. Using in vitro experiments and molecular dynamics simulations, we found that the ATP effects can be attributed to a direct interaction of ATP with the nucleobases of the CAG-RNA rather than ATP acting as "a fuel" for helicase activity. ATP-driven changes in CAG-RNA homeostasis could be disease-relevant since mitochondrial function is affected in HD disease progression leading to a decline in cellular ATP levels.

Structural basis for intrinsic transcription termination.

Efficient and accurate termination is required for gene transcription in all living organisms1,2. Cellular RNA polymerases in both bacteria and eukaryotes can terminate their transcription through a factor-independent termination pathway3,4-called intrinsic termination transcription in bacteria-in which RNA polymerase recognizes terminator sequences, stops nucleotide addition and releases nascent RNA spontaneously. Here we report a set of single-particle cryo-electron microscopy structures of Escherichia coli transcription intrinsic termination complexes representing key intermediate states of the event. The structures show how RNA polymerase pauses at terminator sequences, how the terminator RNA hairpin folds inside RNA polymerase, and how RNA polymerase rewinds the transcription bubble to release RNA and then DNA. These macromolecular snapshots define a structural mechanism for bacterial intrinsic termination and a pathway for RNA release and DNA collapse that isrelevant for factor-independent termination by all RNA polymerases.

Impact of cholesterol and Lumacaftor on the folding of CFTR helical hairpins

Cystic fibrosis (CF) is caused by mutations in the gene that codes for the chloride channel cystic fibrosis transmembrane conductance regulator (CFTR). Recent advances in CF treatment have included use of small-molecule drugs known as modulators, such as Lumacaftor (VX-809), but their detailed mechanism of action and interplay with the surrounding lipid membranes, including cholesterol, remain largely unknown. To examine these phenomena and guide future modulator development, we prepared a set of wild type (WT) and mutant helical hairpin constructs consisting of CFTR transmembrane (TM) segments 3 and 4 and the intervening extracellular loop (termed TM3/4 hairpins) that represent minimal membrane protein tertiary folding units. These hairpin variants, including CF-phenotypic loop mutants E217G and Q220R, and membrane-buried mutant V232D, were reconstituted into large unilamellar phosphatidylcholine (POPC) vesicles, and into corresponding vesicles containing 70 mol% POPC +30 mol% cholesterol, and studied by single-molecule FRET and circular dichroism experiments. We found that the presence of 30 mol% cholesterol induced an increase in helicity of all TM3/4 hairpins, suggesting an increase in bilayer cross-section and hence an increase in the depth of membrane insertion compared to pure POPC vesicles. Importantly, when we added the corrector VX-809, regardless of the presence or absence of cholesterol, all mutants displayed folding and helicity largely indistinguishable from the WT hairpin. Fluorescence spectroscopy measurements suggest that the corrector alters lipid packing and water accessibility. We propose a model whereby VX-809 shields the protein from the lipid environment in a mutant-independent manner such that the WT scaffold prevails. Such ‘normalization’ to WT conformation is consistent with the action of VX-809 as a protein-folding chaperone.

Organized states arising from compression of single semiflexible polymer chains in nanochannels.

We use molecular dynamics simulation to probe the nonequilibrium physics of single nanochannel-confined semiflexible polymers in a homogeneous flow field. The flow field compresses the polymer against the end of the nanochannel, simulating an experiment of a nanochannel confined chain compressed against a slit barrier. The flow-based compression gives rise to a packing of the chain against the channel end that possesses a striking organization, consisting of interweaving of folds and circular coils. For stiff chains at low flow, we find that the organization is dominated by repeated hairpin folds. For stiff chains at higher flow, we observe that circular coils arise along with the folds, with folding and coiling domains becoming interwoven at the highest flow speeds. Chain organization is retained even when the chain persistence length is on order of the channel width. We show that the global polymer organization, consisting of a number of defined folds and coiled loops, arises from the minimization of the total chain free energy.

RNA hairpin folding dynamics in a cell

RNA folding plays an essential role for temperature dependent functions in cells. For example, RNA thermometers regulate translation of heat shock and virulence genes by folding and unfolding. Here, we employed the fourU RNA thermometer motif of Salmonella, a hairpin—structured RNA to study in-cell folding stability of an RNA hairpin and its in vitro stability which is perturbed by crowding and sticking agents. To discover the solvation environment that accounts for cellular crowding and sticking we conduct a two-dimensional scan of mixtures of artificial crowding and sticking agents, PEG10K (polyethylene glycol, 0-300 mg/ml) and M-PER buffer (Mammalian Protein Extraction Reagent, 0-80%).

Observing the base-by-base search for native structure along transition paths during the folding of single nucleic acid hairpins

Biomolecular folding involves searching among myriad possibilities for the native conformation, but the elementary steps expected from theory for this search have never been detected directly. We probed the dynamics of folding at high resolution using optical tweezers, measuring individual trajectories as nucleic acid hairpins passed through the high-energy transition states that dominate kinetics and define folding mechanisms. We observed brief but ubiquitous pauses in the transition states, with a dwell time distribution that matched microscopic theories of folding quantitatively. The sequence dependence suggested that pauses were dominated by microbarriers from nonnative conformations during the search by each nucleotide residue for the native base-pairing conformation. Furthermore, the pauses were position dependent, revealing subtle local variations in energy-landscape roughness and allowing the diffusion coefficient describing the microscopic dynamics within the barrier to be found without reconstructing the shape of the energy landscape. These results show how high-resolution measurements can elucidate key microscopic events during folding to test fundamental theories of folding.

Tumuc1: A New Accurate DNA Force Field Consistent with High-Level Quantum Chemistry.

An accurate molecular mechanics force field forms the basis of Molecular Dynamics simulations to obtain a realistic view of the structure and dynamics of biomolecules such as DNA. Although frequently updated to improve agreement with available experimental data, DNA force fields still rely in part on parameters introduced more than 20 years ago. We have developed an entirely new DNA force field, Tumuc1, derived from quantum mechanical calculations to obtain a consistent set of bonded parameters and partial atomic charges. The performance of the force field was extensively tested on a variety of DNA molecules. It excels in accuracy of B-DNA simulations but also performs very well on other types of DNA structures and structure formation processes such as hairpin folding, duplex formation, and dynamics of DNA-protein complexes. It can complement existing force fields in order to provide an increasingly accurate description of the structure and dynamics of DNA during simulation studies.

Folding in Place: Design of β-Strap Motifs to Stabilize the Folding of Hairpins with Long Loops.

Despite their pivotal role in defining antibody affinity and protein function, β-hairpins harboring long noncanonical loops remain synthetically challenging because of the large entropic penalty associated with their conformational folding. Little is known about the contribution and impact of stabilizing motifs on the folding of β-hairpins with loops of variable length and plasticity. Here, we report a design of minimalist β-straps (strap = strand + cap) that offset the entropic cost of long-loop folding. The judicious positioning of noncovalent interactions (hydrophobic cluster and salt-bridge) within the novel 8-mer β-strap design RW(V/H)W···WVWE stabilizes hairpins with up to 10-residue loops of varying degrees of plasticity (Tm up to 52 °C; 88 ± 1% folded at 18 °C). This "hyper" thermostable β-strap outperforms the previous gold-standard technology of β-strand-β-cap (16-mer) and provides a foundation for producing new classes of long hairpins as a viable and practical alternative to macrocyclic peptides.

Complex Conformational Dynamics of the Heart Failure-Associated Pre-miRNA-377 Hairpin Revealed by Single-Molecule Optical Tweezers.

Pre-miRNA-377 is a hairpin-shaped regulatory RNA associated with heart failure. Here, we use single-molecule optical tweezers to unzip pre-miRNA-377 and study its stability and dynamics. We show that magnesium ions have a strong stabilizing effect, and that sodium ions stabilize the hairpin more than potassium ions. The hairpin unfolds in a single step, regardless of buffer composition. Interestingly, hairpin folding occurs either in a single step (type 1) or through the formation of intermediates, in multiple steps (type 2) or gradually (type 3). Type 3 occurs only in the presence of both sodium and magnesium, while type 1 and 2 take place in all buffers, with type 1 being the most prevalent. By reducing the size of the native hairpin loop from fourteen to four nucleotides, we demonstrate that the folding heterogeneity originates from the large size of the hairpin loop. Further, while efficient pre-miRNA-377 binders are lacking, we demonstrate that the recently developed C2 ligand displays bimodal activity: it enhances the mechanical stability of the pre-miRNA-377 hairpin and perturbs its folding. The knowledge regarding pre-miRNA stability and dynamics that we provide is important in understanding its regulatory function and how it can be modulated to achieve a therapeutic effect, e.g., in heart failure treatment.

Conformational Dynamics of Poly(T) Single-Stranded DNA at the Single-Molecule Level.

The conformational dynamics of single-stranded DNA (ss-DNA) are implicated in the mechanisms of several key biological processes such as DNA replication and damage repair and have been modeled with those of semiflexible or flexible polymer. The high flexibility and customizability of ss-DNA also make it an excellent polymeric material for materials engineering. Polythymidine (poly(T)) is an excellent model ss-DNA as a flexible polymer since it does not form any secondary structure. However, only limited experimental results have been reported of poly(T) conformational dynamics with a very short length that is not relevant to the aforementioned processes and applications. Here, we provide the first experimental results of the conformational dynamics of poly(T) with lengths in the range of 130-170 nucleotides at the single-molecule level. Our experiments are based on single-molecule FRET and a DNA hairpin structure of which the folding kinetics are governed by the conformational dynamics of poly(T). We found that the folding kinetics deviate far from those of a flexible polymer model with a harmonic bending potential. To this end, we derived a simple model for the kinetics of DNA hairpin folding from the self-avoiding-walk (SAW). Our model describes the conformational dynamics of poly(T) very well and enables estimation of the conformational dimensionality. The estimated dimensionalities suggest that ss-DNA is completely flexible at 100 mM or a higher NaCl concentration, but not at 50 mM. These results will help understand the conformational dynamics of ss-DNA implicated in several key biological processes and maximize the utility of ss-DNA for materials engineering. Also, our system and model provide an excellent platform to investigate the conformational dynamics of ss-DNA.

Correlative Single-Molecule Force-FRET Measurements of DNA Hairpin Conformational Dynamics

The biological function of many macromolecules such as proteins and nucleic acids is tightly coupled to their conformation and their conformational dynamics. Therefore, studying how macromolecules fold and undergo conformational changes is crucial to understand the underlying mechanisms of biological processes. In particular DNA/RNA hairpins play essential regulatory roles in vivo, since they can regulate the interaction with nucleic acid processing proteins in processes such as transcription, recombination or replication. Therefore, understanding the principles of their structure formation is important for elucidating those biological roles. Single-molecule force spectroscopy represents an ideal tool for such studies because of its unique capability to isolate individual biomolecules and observe conformational transitions in real-time, while providing the control over the population distribution as opposed to approaches in solution. Here we report a study on FRET-labelled DNA hairpin folding transitions using high-resolution optical tweezers combined with fluorescence microscopy. We present equilibrium measurements at unprecedent trap distance stability, which allows to measure conformational state distributions over hours without altering the energy landscape. Measuring over longer time periods while keeping the trap distance within nanometer precision is crucial to capture the rare conformational states that might be overlooked in measurements over a timescale that is limited by instrumental drift. These measurements were enabled by integrating trap-distance and temperature stabilization feedback loops to achieve the ultimate stability in distance between the traps. Kinetic, energetic and structural properties of these states could be determined quantitatively from a correlative information source – force and fluorescence data, thereby providing new insights into folding pathways and energy landscape. This study shows how correlating force spectroscopy measurements with fluorescence microscopy could complement the investigation of the prevalence and energy landscape of certain conformational states with extra three-dimensional structural information of the macromolecule of interest.

Modeling Loop Composition and Ion Concentration Effects in RNA Hairpin Folding Stability

The ability to accurately predict RNA hairpin structure and stability for different loop sequences and salt conditions is important for understanding, modeling, and designing larger RNA folds. However, traditional RNA secondary structure models cannot treat loop-sequence and ionic effects on RNA hairpin folding. Here, we describe a general, three-dimensional (3D) conformation-based computational method for modeling salt concentration-dependent conformational distributions and the detailed 3D structures for a set of three RNA hairpins that contain a variable, 15-nucleotide loop sequence. For a given RNA sequence, the new, to our knowledge, method integrates a Vfold2D two-dimensional structure folding model with IsRNA coarse-grained molecular dynamics 3D folding simulations and Monte Carlo tightly bound ion estimations of ion-mediated electrostatic interactions. The model predicts free-energy landscapes for the different RNA hairpin-forming sequences with variable salt conditions. The theoretically predicted results agree with the experimental fluorescence measurements, validating the strategy. Furthermore, the theoretical model goes beyond the experimental results by enabling in-depth 3D structural analysis, revealing energetic mechanisms for the sequence- and salt-dependent folding stability. Although the computational framework presented here is developed for RNA hairpin systems, the general method may be applied to investigate other RNA systems, such as multiway junctions or pseudoknots in mixed metal ion solutions.

BamA is required for autotransporter secretion

BackgroundIn Gram-negative bacteria, type Va and Vc autotransporters are proteins that contain both a secreted virulence factor (the “passenger” domain) and a β-barrel that aids its export. While it is known that the folding and insertion of the β-barrel domain utilize the β-barrel assembly machinery (BAM) complex, how the passenger domain is secreted and folded across the membrane remains to be determined. The hairpin model states that passenger domain secretion occurs independently through the fully-formed and membrane-inserted β-barrel domain via a hairpin folding intermediate. In contrast, the BamA-assisted model states that the passenger domain is secreted through a hybrid of BamA, the essential subunit of the BAM complex, and the β-barrel domain of the autotransporter. MethodsTo ascertain the models' plausibility, we have used molecular dynamics to simulate passenger domain secretion for two autotransporters, EspP and YadA. ResultsWe observed that each protein's β-barrel is unable to accommodate the secreting passenger domain in a hairpin configuration without major structural distortions. Additionally, the force required for secretion through EspP's β-barrel is more than that through the BamA β-barrel. ConclusionsSecretion of autotransporters most likely occurs through an incompletely formed β-barrel domain of the autotransporter in conjunction with BamA. General significanceSecretion of virulence factors is a process used by practically all pathogenic Gram-negative bacteria. Understanding this process is a necessary step towards limiting their infectious capacity.

Memory Effects in Single-Molecule Force Spectroscopy Measurements of Biomolecular Folding

Folding is generally assumed to be a Markov process, without memory. When the molecular motion is coupled to that of a probe as in single-molecule force spectroscopy (SMFS) experiments, however, theory predicts that the coupling to a second Markov process should induce memory when monitoring a projection of the full multi-dimensional motion onto a reduced coordinate. We developed a method to evaluate the time constant of the induced memory from its effects on the autocorrelation function, which can be readily determined from experimental data. Applying this method to both simulated SMFS measurements and experimental trajectories of DNA hairpin folding measured by optical tweezers, we validated the prediction that the linker induces memory. For these measurements, the timescale of the induced memory was found to be similar to the time required for the force probe to respond to changes in the molecule, and in the regime where the experimentally observed dynamics were not significantly perturbed by probe-molecule coupling artifacts. Memory effects are thus a general feature of SMFS measurements induced by the mechanical connection between the molecule and force probe that should be considered when interpreting experimental data.

Tension propagation during DNA hairpin folding

Characterization of DNA/RNA folding dynamics is an interesting problem involving the complex force balance between the relaxation of the growing duplex and the stretched single-stranded segments. In a previous paper, we reported that the helicity plays an important role in determining the folding time for a chain of size N, where is the Flory exponent. Here we analyze this process in further detail by using molecular dynamics, with particular emphasis on tension propagation along a single strand on the unfolded segment. We directly observe that the single-strand segments are always stretched during the folding process by the tension induced by base pair formation, propagating ahead of the y -junction. Our molecular dynamics simulations verify the existence of a stem-flower structure in the unfolded segments, with a power-law dependence of the stem length Ns on the duplex length n. For the longest hairpin structures considered, we also find that the power-law regime terminates before the first half of the folding is completed. We demonstrate that the structure’s helicity is an integral aspect of the folding dynamics by comparing our findings with those obtained from a hypothetical non-helical model.

Spontaneous Tilting Transition in Liquid-Crystalline Polymer Brushes

The equilibrium phase behavior of main-chain liquid-crystalline polymer brushes in good solvent is investigated using self-consistent field theory. The calculation assumes semidilute brushes, where the grafting density is sufficient to cause overlap among the polymers but low enough to allow the solvent to be treated implicitly. The polymers are modeled as semiflexible wormlike chains with Maier–Saupe interactions. Previous calculations based on freely jointed chains have predicted that extended brushes collapse into folded nematic brushes, as the liquid-crystalline interaction strength is increased. For the more sophisticated model of wormlike chains, which penalizes hairpin folds, we find that extended brushes become unstable with respect to tilting prior to the onset of folding. This implies that the brushes spontaneously collapse by a symmetry-breaking transition, where the chains tilt rather than fold.

Super-Resolution Optical Lithography with DNA.

We developed an efficient, versatile, and accessible super-resolution microscopy method to construct a nanoparticle assembly at a spatial resolution below the optical diffraction limit. The method utilizes DNA and a photoactivated DNA cross-linker. Super-resolution optical techniques have been used only as a means to make measurements below the light diffraction limit. Furthermore, no optical technique is currently available to construct nanoparticle assemblies with a precisely designed shape and internal structure at a resolution of a few tens of nanometers (nm). Here we demonstrate that we can fulfill this deficiency by utilizing spontaneous structural dynamics of DNA hairpins combined with single-molecule fluorescence resonance energy transfer (smFRET) microscopy and a photoactivated DNA cross-linker. The stochastic fluorescence blinking due to the spontaneous folding and unfolding motions of DNA hairpins enables us to precisely localize a folded hairpin and solidify it only when it is within a predesigned target area whose size is below the diffraction limit. As the method is based on an optical microscope and an easily clickable DNA cross-linking reagent, it will provide an efficient means to create large nanoparticle assemblies with a shape and internal structure at an optical super-resolution, opening a wide window of opportunities toward investigating their photophysical and optoelectronic properties and developing novel devices.