Extracting kinetics from single-molecule force spectroscopy: nanopore unzipping of DNA hairpins

Single-molecule force experiments provide powerful new tools to explore biomolecular interactions. Here, we describe a systematic procedure for extracting kinetic information from force-spectroscopy experiments, and apply it to nanopore unzipping of individual DNA hairpins. Two types of measurements are considered: unzipping at constant voltage, and unzipping at constant voltage-ramp speeds. We perform a global maximum-likelihood analysis of the experimental data at low-to-intermediate ramp speeds. To validate the theoretical models, we compare their predictions with two independent sets of data, collected at high ramp speeds and at constant voltage, by using a quantitative relation between the two types of measurements. Microscopic approaches based on Kramers theory of diffusive barrier crossing allow us to estimate not only intrinsic rates and transition state locations, as in the widely used phenomenological approach based on Bell's formula, but also free energies of activation. The problem of extracting unique and accurate kinetic parameters of a molecular transition is discussed in light of the apparent success of the microscopic theories in reproducing the experimental data.     

Disease‐linked TDP‐43 hyperphosphorylation suppresses TDP‐43 condensation and aggregation

Post‐translational modifications (PTMs) have emerged as key modulators of protein phase separation and have been linked to protein aggregation in neurodegenerative disorders. The major aggregating protein in amyotrophic lateral sclerosis and frontotemporal dementia, the RNA‐binding protein TAR DNA‐binding protein (TDP‐43), is hyperphosphorylated in disease on several C‐terminal serine residues, a process generally believed to promote TDP‐43 aggregation. Here, we however find that Casein kinase 1δ‐mediated TDP‐43 hyperphosphorylation or C‐terminal phosphomimetic mutations reduce TDP‐43 phase separation and aggregation, and instead render TDP‐43 condensates more liquid‐like and dynamic. Multi‐scale molecular dynamics simulations reveal reduced homotypic interactions of TDP‐43 low‐complexity domains through enhanced solvation of phosphomimetic residues. Cellular experiments show that phosphomimetic substitutions do not affect nuclear import or RNA regulatory functions of TDP‐43, but suppress accumulation of TDP‐43 in membrane‐less organelles and promote its solubility in neurons. We speculate that TDP‐43 hyperphosphorylation may be a protective cellular response to counteract TDP‐43 aggregation.     

Expression of Smad4 in the FaDu cell line partially restores TGF-beta growth inhibition but is not sufficient to regulate fibronectin expression or suppress tumorigenicity

Mutations of the Smad4 gene, a member of a group of TGF-beta signal transduction components, occur in several types of cancer suggesting that its inactivation significantly affects TGF-beta responsiveness in these tumors. To further investigate the role of Smad4 with respect to TGF-beta signaling and carcinogenesis, we re-expressed the Smad4 gene in the Smad4-deficient cancer cell line FaDu by microcell-mediated chromosome transfer (MMCT) and retroviral infection to closely approximate physiological protein levels. The Smad4-expressing FaDu clones were then evaluated for TGF-beta responsiveness to assess the role of Smad4 in TGF-beta-induced growth inhibition and target gene regulation. We found that the re-expression of the Smad4 gene by either method partially restored TGF-beta responsiveness in FaDu cells with respect to both growth inhibition and expression of p21WAF1/CIP1 and p15INK4B. However, only the microcell hybrids showed growth retardation in organotypic raft culture and an enhanced ability to upregulate fibronectin. In contrast, the re-expression of Smad4 by either method failed to suppress tumorigenicity. These results suggest that in addition to a homozygous deletion of Smad4, FaDu cells contain additional defects within the TGF-beta signaling pathway, thereby limiting the extent of TGF-beta responsiveness upon Smad4 re-expression and perhaps accounting for the inability to induce p15INK4B to a high level. They also demonstrate the advantages of providing a physiological extracellular environment, when assessing TGFbeta responsiveness.     

Helix nucleation kinetics from molecular simulations in explicit solvent

We study the reversible folding/unfolding of short Ala and Gly-based peptides by molecular dynamics simulations of all-atom models in explicit water solvent. A kinetic analysis shows that the formation of a first alpha-helical turn occurs within 0.1-1 ns, in agreement with the analyses of laser temperature jump experiments. The unfolding times exhibit Arrhenius temperature dependence. For a rapidly nucleating all-Ala peptide, the helix nucleation time depends only weakly on temperature. For a peptide with enthalpically competing turn-like structures, helix nucleation exhibits an Arrhenius temperature dependence, corresponding to the unfolding of enthalpic traps in the coil ensemble. An analysis of structures in a "transition-state ensemble" shows that helix-to-coil transitions occur predominantly through breaking of hydrogen bonds at the helix ends, particularly at the C-terminus. The temperature dependence of the transition-state ensemble and the corresponding folding/unfolding pathways illustrate that folding mechanisms can change with temperature, possibly complicating the interpretation of high-temperature unfolding simulations. The timescale of helix formation is an essential factor in molecular models of protein folding. The rapid helix nucleation observed here suggests that transient helices form early in the folding event.     

Hydrophobicity maps of the N-peptide coiled coil of HIV-1 gp41

Blocking HIV-1 viral entry into the host cell offers a promising new strategy for interfering with the HIV-1 life cycle. A major target of inhibitor design is to prevent binding of fusogenic gp41 C-peptides to the trimeric coiled coil of fusion-active N-peptides. Here, we map the hydrophobic character of the binding surface of the IQN17 peptide, a soluble analogue of the N-peptide coiled coil. The local binding affinity for a hydrophobic probe is determined by three methods: a hydrophobic force field, and molecular dynamics in solution analyzed by test particle insertion and inhomogeneous information theory. The regions of highest calculated hydrophobicity overlap with the positions of the hydrophobic anchor residues of the native C-peptides, and of two known inhibitors. Additional binding sites not exploited by these inhibitors are identified, and modifications for enhancing their binding affinity are suggested.     

Cysteine cross-linking in native membranes establishes the transmembrane architecture of Ire1

The ER is a key organelle of membrane biogenesis and crucial for the folding of both membrane and secretory proteins. Sensors of the unfolded protein response (UPR) monitor the unfolded protein load in the ER and convey effector functions for maintaining ER homeostasis. Aberrant compositions of the ER membrane, referred to as lipid bilayer stress, are equally potent activators of the UPR. How the distinct signals from lipid bilayer stress and unfolded proteins are processed by the conserved UPR transducer Ire1 remains unknown. Here, we have generated a functional, cysteine-less variant of Ire1 and performed systematic cysteine cross-linking experiments in native membranes to establish its transmembrane architecture in signaling-active clusters. We show that the transmembrane helices of two neighboring Ire1 molecules adopt an X-shaped configuration independent of the primary cause for ER stress. This suggests that different forms of stress converge in a common, signaling-active transmembrane architecture of Ire1.     

Ion transport through membrane-spanning nanopores studied by molecular dynamics simulations and continuum electrostatics calculations

Narrow hydrophobic regions are a common feature of biological channels, with possible roles in ion-channel gating. We study the principles that govern ion transport through narrow hydrophobic membrane pores by molecular dynamics simulation of model membranes formed of hexagonally packed carbon nanotubes. We focus on the factors that determine the energetics of ion translocation through such nonpolar nanopores and compare the resulting free-energy barriers for pores with different diameters corresponding to the gating regions in closed and open forms of potassium channels. Our model system also allows us to compare the results from molecular dynamics simulations directly to continuum electrostatics calculations. Both simulations and continuum calculations show that subnanometer wide pores pose a huge free-energy barrier for ions, but a small increase in the pore diameter to approximately 1 nm nearly eliminates that barrier. We also find that in those wider channels the ion mobility is comparable to that in the bulk phase. By calculating local electrostatic potentials, we show that the long range Coulomb interactions of ions are strongly screened in the wide water-filled channels. Whereas continuum calculations capture the overall energetics reasonably well, the local water structure, which is not accounted for in this model, leads to interesting effects such as the preference of hydrated ions to move along the pore wall rather than through the center of the pore.