Structural interactions of a voltage sensor toxin with lipid membranes

Protein toxins from tarantula venom alter the activity of diverse ion channel proteins, including voltage, stretch, and ligand-activated cation channels. Although tarantula toxins have been shown to partition into membranes, and the membrane is thought to play an important role in their activity, the structural interactions between these toxins and lipid membranes are poorly understood. Here, we use solid-state NMR and neutron diffraction to investigate the interactions between a voltage sensor toxin (VSTx1) and lipid membranes, with the goal of localizing the toxin in the membrane and determining its influence on membrane structure. Our results demonstrate that VSTx1 localizes to the headgroup region of lipid membranes and produces a thinning of the bilayer. The toxin orients such that many basic residues are in the aqueous phase, all three Trp residues adopt interfacial positions, and several hydrophobic residues are within the membrane interior. One remarkable feature of this preferred orientation is that the surface of the toxin that mediates binding to voltage sensors is ideally positioned within the lipid bilayer to favor complex formation between the toxin and the voltage sensor.Protein toxins from venomous organisms have been invaluable tools for studying the ion channel proteins they target. For example, in the case of voltage-activated potassium (Kv) channels, pore-blocking scorpion toxins were used to identify the pore-forming region of the channel (1, 2), and gating modifier tarantula toxins that bind to S1–S4 voltage-sensing domains have helped to identify structural motifs that move at the protein–lipid interface (3–5). In many instances, these toxin–channel interactions are highly specific, allowing them to be used in target validation and drug development (6–8).Tarantula toxins are a particularly interesting class of protein toxins that have been found to target all three families of voltage-activated cation channels (3, 9–12), stretch-activated cation channels (13–15), as well as ligand-gated ion channels as diverse as acid-sensing ion channels (ASIC) (16–21) and transient receptor potential (TRP) channels (22, 23). The tarantula toxins targeting these ion channels belong to the inhibitor cystine knot (ICK) family of venom toxins that are stabilized by three disulfide bonds at the core of the molecule (16, 17, 24–31). Although conventional tarantula toxins vary in length from 30 to 40 aa and contain one ICK motif, the recently discovered double-knot toxin (DkTx) that specifically targets TRPV1 channels contains two separable lobes, each containing its own ICK motif (22, 23).One unifying feature of all tarantula toxins studied thus far is that they act on ion channels by modifying the gating properties of the channel. The best studied of these are the tarantula toxins targeting voltage-activated cation channels, where the toxins bind to the S3b–S4 voltage sensor paddle motif (5, 32–36), a helix-turn-helix motif within S1–S4 voltage-sensing domains that moves in response to changes in membrane voltage (37–41). Toxins binding to S3b–S4 motifs can influence voltage sensor activation, opening and closing of the pore, or the process of inactivation (4, 5, 36, 42–46). The tarantula toxin PcTx1 can promote opening of ASIC channels at neutral pH (16, 18), and DkTx opens TRPV1 in the absence of other stimuli (22, 23), suggesting that these toxin stabilize open states of their target channels.For many of these tarantula toxins, the lipid membrane plays a key role in the mechanism of inhibition. Strong membrane partitioning has been demonstrated for a range of toxins targeting S1–S4 domains in voltage-activated channels (27, 44, 47–50), and for GsMTx4 (14, 50), a tarantula toxin that inhibits opening of stretch-activated cation channels in astrocytes, as well as the cloned stretch-activated Piezo1 channel (13, 15). In experiments on stretch-activated channels, both the d- and l-enantiomers of GsMTx4 are active (14, 50), implying that the toxin may not bind directly to the channel. In addition, both forms of the toxin alter the conductance and lifetimes of gramicidin channels (14), suggesting that the toxin inhibits stretch-activated channels by perturbing the interface between the membrane and the channel. In the case of Kv channels, the S1–S4 domains are embedded in the lipid bilayer and interact intimately with lipids (48, 51, 52) and modification in the lipid composition can dramatically alter gating of the channel (48, 53–56). In one study on the gating of the Kv2.1/Kv1.2 paddle chimera (53), the tarantula toxin VSTx1 was proposed to inhibit Kv channels by modifying the forces acting between the channel and the membrane. Although these studies implicate a key role for the membrane in the activity of Kv and stretch-activated channels, and for the action of tarantula toxins, the influence of the toxin on membrane structure and dynamics have not been directly examined. The goal of the present study was to localize a tarantula toxin in membranes using structural approaches and to investigate the influence of the toxin on the structure of the lipid bilayer.     

Kinetics of nucleotide-dependent structural transitions in the kinesin-1 hydrolysis cycle

To dissect the kinetics of structural transitions underlying the stepping cycle of kinesin-1 at physiological ATP, we used interferometric scattering microscopy to track the position of gold nanoparticles attached to individual motor domains in processively stepping dimers. Labeled heads resided stably at positions 16.4 nm apart, corresponding to a microtubule-bound state, and at a previously unseen intermediate position, corresponding to a tethered state. The chemical transitions underlying these structural transitions were identified by varying nucleotide conditions and carrying out parallel stopped-flow kinetics assays. At saturating ATP, kinesin-1 spends half of each stepping cycle with one head bound, specifying a structural state for each of two rate-limiting transitions. Analysis of stepping kinetics in varying nucleotides shows that ATP binding is required to properly enter the one-head–bound state, and hydrolysis is necessary to exit it at a physiological rate. These transitions differ from the standard model in which ATP binding drives full docking of the flexible neck linker domain of the motor. Thus, this work defines a consensus sequence of mechanochemical transitions that can be used to understand functional diversity across the kinesin superfamily.Kinesin-1 is a motor protein that steps processively toward microtubule plus-ends, tracking single protofilaments and hydrolyzing one ATP molecule per step (1–6). Step sizes corresponding to the tubulin dimer spacing of 8.2 nm are observed when the molecule is labeled by its C-terminal tail (7–10) and to a two-dimer spacing of 16.4 nm when a single motor domain is labeled (4, 11, 12), consistent with the motor walking in a hand-over-hand fashion. Kinesin has served as an important model system for advancing single-molecule techniques (7–10) and is clinically relevant for its role in neurodegenerative diseases (13), making dissection of its step a popular ongoing target of study.Despite decades of work, many essential components of the mechanochemical cycle remain disputed, including (i) how much time kinesin-1 spends in a one-head–bound (1HB) state when stepping at physiological ATP concentrations, (ii) whether the motor waits for ATP in a 1HB or two-heads–bound (2HB) state, and (iii) whether ATP hydrolysis occurs before or after tethered head attachment (4, 11, 14–20). These questions are important because they are fundamental to the mechanism by which kinesins harness nucleotide-dependent structural changes to generate mechanical force in a manner optimized for their specific cellular tasks. Addressing these questions requires characterizing a transient 1HB state in the stepping cycle in which the unattached head is located between successive binding sites on the microtubule. This 1HB intermediate is associated with the force-generating powerstroke of the motor and underlies the detachment pathway that limits motor processivity. Optical trapping (7, 19, 21, 22) and single-molecule tracking studies (4, 8–11) have failed to detect this 1HB state during stepping. Single-molecule fluorescence approaches have detected a 1HB intermediate at limiting ATP concentrations (11, 12, 14, 15), but apart from one study that used autocorrelation analysis to detect a 3-ms intermediate (17), the 1HB state has been undetectable at physiological ATP concentrations.Single-molecule microscopy is a powerful tool for studying the kinetics of structural changes in macromolecules (23). Tracking steps and potential substeps for kinesin-1 at saturating ATP has until now been hampered by the high stepping rates of the motor (up to 100 s⁻¹), which necessitates high frame rates, and the small step size (8.2 nm), which necessitates high spatial precision (7). Here, we apply interferometric scattering microscopy (iSCAT), a recently established single-molecule tool with high spatiotemporal resolution (24–27) to directly visualize the structural changes underlying kinesin stepping. By labeling one motor domain in a dimeric motor, we detect a 1HB intermediate state in which the tethered head resides over the bound head for half the duration of the stepping cycle at saturating ATP. We further show that at physiological stepping rates, ATP binding is required to enter this 1HB state and that ATP hydrolysis is required to exit it. This work leads to a significant revision of the sequence and kinetics of mechanochemical transitions that make up the kinesin-1 stepping cycle and provides a framework for understanding functional diversity across the kinesin superfamily.     

α-Synuclein assembles into higher-order multimers upon membrane binding to promote SNARE complex formation

Physiologically, α-synuclein chaperones soluble NSF attachment protein receptor (SNARE) complex assembly and may also perform other functions; pathologically, in contrast, α-synuclein misfolds into neurotoxic aggregates that mediate neurodegeneration and propagate between neurons. In neurons, α-synuclein exists in an equilibrium between cytosolic and membrane-bound states. Cytosolic α-synuclein appears to be natively unfolded, whereas membrane-bound α-synuclein adopts an α-helical conformation. Although the majority of studies showed that cytosolic α-synuclein is monomeric, it is unknown whether membrane-bound α-synuclein is also monomeric, and whether chaperoning of SNARE complex assembly by α-synuclein involves its cytosolic or membrane-bound state. Here, we show using chemical cross-linking and fluorescence resonance energy transfer (FRET) that α-synuclein multimerizes into large homomeric complexes upon membrane binding. The FRET experiments indicated that the multimers of membrane-bound α-synuclein exhibit defined intermolecular contacts, suggesting an ordered array. Moreover, we demonstrate that α-synuclein promotes SNARE complex assembly at the presynaptic plasma membrane in its multimeric membrane-bound state, but not in its monomeric cytosolic state. Our data delineate a folding pathway for α-synuclein that ranges from a monomeric, natively unfolded form in cytosol to a physiologically functional, multimeric form upon membrane binding, and show that only the latter but not the former acts as a SNARE complex chaperone at the presynaptic terminal, and may protect against neurodegeneration.α-Synuclein is an abundant presynaptic protein that physiologically acts to promote soluble NSF attachment protein receptor (SNARE) complex assembly in vitro and in vivo (1–3). Point mutations in α-synuclein (A30P, E46K, H50Q, G51D, and A53T) as well as α-synuclein gene duplications and triplications produce early-onset Parkinson''s disease (PD) (4–10). Moreover, α-synuclein is a major component of intracellular protein aggregates called Lewy bodies, which are pathological hallmarks of neurodegenerative disorders such as PD, Lewy body dementia, and multiple system atrophy (11–14). Strikingly, neurotoxic α-synuclein aggregates propagate between neurons during neurodegeneration, suggesting that such α-synuclein aggregates are not only intrinsically neurotoxic but also nucleate additional fibrillization (15–18).α-Synuclein is highly concentrated in presynaptic terminals where α-synuclein exists in an equilibrium between a soluble and a membrane-bound state, and is associated with synaptic vesicles (19–22). The labile association of α-synuclein with membranes (23, 24) suggests that binding of α-synuclein to synaptic vesicles, and its dissociation from these vesicles, may regulate its physiological function. Membrane-bound α-synuclein assumes an α-helical conformation (25–32), whereas cytosolic α-synuclein is natively unfolded and monomeric (refs. 25, 26, 31, and 32; however, see refs. 33 and 34 and Discussion for a divergent view). Membrane binding by α-synuclein is likely physiologically important because in in vitro experiments, α-synuclein remodels membranes (35, 36), influences lipid packing (37, 38), and induces vesicle clustering (39). Moreover, membranes were found to be important for the neuropathological effects of α-synuclein (40–44).However, the relation of membrane binding to the in vivo function of α-synuclein remains unexplored, and it is unknown whether α-synuclein binds to membranes as a monomer or oligomer. Thus, in the present study we have investigated the nature of the membrane-bound state of α-synuclein and its relation to its physiological function in SNARE complex assembly. We found that soluble monomeric α-synuclein assembles into higher-order multimers upon membrane binding and that membrane binding of α-synuclein is required for its physiological activity in promoting SNARE complex assembly at the synapse.     

Dopamine release from transplanted neural stem cells in Parkinsonian rat striatum in vivo

Embryonic stem cell-based therapies exhibit great potential for the treatment of Parkinson’s disease (PD) because they can significantly rescue PD-like behaviors. However, whether the transplanted cells themselves release dopamine in vivo remains elusive. We and others have recently induced human embryonic stem cells into primitive neural stem cells (pNSCs) that are self-renewable for massive/transplantable production and can efficiently differentiate into dopamine-like neurons (pNSC–DAn) in culture. Here, we showed that after the striatal transplantation of pNSC–DAn, (i) pNSC–DAn retained tyrosine hydroxylase expression and reduced PD-like asymmetric rotation; (ii) depolarization-evoked dopamine release and reuptake were significantly rescued in the striatum both in vitro (brain slices) and in vivo, as determined jointly by microdialysis-based HPLC and electrochemical carbon fiber electrodes; and (iii) the rescued dopamine was released directly from the grafted pNSC–DAn (and not from injured original cells). Thus, pNSC–DAn grafts release and reuptake dopamine in the striatum in vivo and alleviate PD symptoms in rats, providing proof-of-concept for human clinical translation.Parkinson’s disease (PD) is a chronic progressive neurodegenerative disorder characterized by the specific loss of dopaminergic neurons in the substantia nigra pars compacta and their projecting axons, resulting in loss of dopamine (DA) release in the striatum (1). During the last two decades, cell-replacement therapy has proven, at least experimentally, to be a potential treatment for PD patients (2–7) and in animal models (8–15). The basic principle of cell therapy is to restore the DA release by transplanting new DA-like cells. Until recently, obtaining enough transplantable cells was a major bottleneck in the practicability of cell therapy for PD. One possible source is embryonic stem cells (ESCs), which can develop infinitely into self-renewable pluripotent cells with the potential to generate any type of cell, including DA neurons (DAns) (16, 17).Recently, several groups including us have introduced rapid and efficient ways to generate primitive neural stem cells (pNSCs) from human ESCs using small-molecule inhibitors under chemically defined conditions (12, 18, 19). These cells are nonpolarized neuroepithelia and retain plasticity upon treatment with neuronal developmental morphogens. Importantly, pNSCs differentiate into DAns (pNSC–DAn) with high efficiency (∼65%) after patterning by sonic hedgehog (SHH) and fibroblast growth factor 8 (FGF8) in vitro, providing an immediate and renewable source of DAns for PD treatment. Importantly, the striatal transplantation of human ESC-derived DA-like neurons, including pNSC–DAn, are able to relieve the motor defects in a PD rat model (11–13, 15, 19–23). Before attempting clinical translation of pNSC–DAn, however, there are two fundamental open questions. (i) Can pNSC–DAn functionally restore the striatal DA levels in vivo? (ii) What cells release the restored DA, pNSC–DAn themselves or resident neurons/cells repaired by the transplants?Regarding question 1, a recent study using nafion-coated carbon fiber electrodes (CFEs) reported that the amperometric current is rescued in vivo by ESC (pNSC–DAn-like) therapy (19). Both norepinephrine (NE) and serotonin are present in the striatum (24, 25). However, CFE amperometry/chronoamperometry alone cannot distinguish DA from other monoamines in vivo, such as NE and serotonin (Fig. S1) (see also refs. 26–28). Considering that the compounds released from grafted ESC-derived cells are unknown, the work of Kirkeby et al. was unable to determine whether DA or other monoamines are responsible for the restored amperometric signal. Thus, the key question of whether pNSC–DAn can rescue DA release needs to be reexamined for the identity of the restored amperometric signal in vivo.Regarding question 2, many studies have proposed that DA is probably released from the grafted cells (8, 12, 13, 20), whereas others have proposed that the grafted stem cells might restore striatal DA levels by rescuing injured original cells (29, 30). Thus, whether the grafted cells are actually capable of synthesizing and releasing DA in vivo must be investigated to determine the future cellular targets (residual cells versus pNSC–DAn) of treatment.To address these two mechanistic questions, advanced in vivo methods of DA identification and DA recording at high spatiotemporal resolution are required. Currently, microdialysis-based HPLC (HPLC) (31–33) and CFE amperometric recordings (34, 35) have been used independently by different laboratories to assess evoked DA release from the striatum in vivo. The major advantage of microdialysis-based HPLC is to identify the substances secreted in the cell-grafted striatum (33), but its spatiotemporal resolution is too low to distinguish the DA release site (residual cells or pNSC–DAn). In contrast, the major advantage of CFE-based amperometry is its very high temporal (ms) and spatial (μm) resolution, making it possible to distinguish the DA release site (residual cells or pNSC–DAn) in cultured cells, brain slices, and in vivo (34–39), but it is unable to distinguish between low-level endogenous oxidizable substances (DA versus serotonin and NE) in vivo.In the present study, we developed a challenging experimental paradigm of combining the two in vivo methods, microdialysis-based HPLC and CFE amperometry, to identify the evoked substance as DA and its release site as pNSC–DAn in the striatum of PD rats.     

From the Cover: The point of no return in vetoing self-initiated movements

In humans, spontaneous movements are often preceded by early brain signals. One such signal is the readiness potential (RP) that gradually arises within the last second preceding a movement. An important question is whether people are able to cancel movements after the elicitation of such RPs, and if so until which point in time. Here, subjects played a game where they tried to press a button to earn points in a challenge with a brain–computer interface (BCI) that had been trained to detect their RPs in real time and to emit stop signals. Our data suggest that subjects can still veto a movement even after the onset of the RP. Cancellation of movements was possible if stop signals occurred earlier than 200 ms before movement onset, thus constituting a point of no return.It has been repeatedly shown that spontaneous movements are preceded by early brain signals (1–8). As early as a second before a simple voluntary movement, a so-called readiness potential (RP) is observed over motor-related brain regions (1–3, 5). The RP was found to precede the self-reported time of the “‘decision’ to act” (ref. 3, p. 623). Similar preparatory signals have been observed using invasive electrophysiology (8, 9) and functional MRI (7, 10), and have been demonstrated also for choices between multiple-response options (6, 7, 10), for abstract decisions (10), for perceptual choices (11), and for value-based decisions (12). To date, the exact nature and causal role of such early signals in decision making is debated (12–20).One important question is whether a person can still exert a veto by inhibiting the movement after onset of the RP (13, 18, 21, 22). One possibility is that the onset of the RP triggers a causal chain of events that unfolds in time and cannot be cancelled. The onset of the RP in this case would be akin to tipping the first stone in a row of dominoes. If there is no chance of intervening, the dominoes will gradually fall one-by-one until the last one is reached. This has been coined a ballistic stage of processing (23, 24). A different possibility is that participants can still terminate the process, akin to taking out a domino at some later stage in the chain and thus preventing the process from completing. Here, we directly tested this in a real-time experiment that required subjects to terminate their decision to move once a RP had been detected by a brain–computer interface (BCI) (25–31).     

Dynamic localization of Mps1 kinase to kinetochores is essential for accurate spindle microtubule attachment

The spindle assembly checkpoint (SAC) is a conserved signaling pathway that monitors faithful chromosome segregation during mitosis. As a core component of SAC, the evolutionarily conserved kinase monopolar spindle 1 (Mps1) has been implicated in regulating chromosome alignment, but the underlying molecular mechanism remains unclear. Our molecular delineation of Mps1 activity in SAC led to discovery of a previously unidentified structural determinant underlying Mps1 function at the kinetochores. Here, we show that Mps1 contains an internal region for kinetochore localization (IRK) adjacent to the tetratricopeptide repeat domain. Importantly, the IRK region determines the kinetochore localization of inactive Mps1, and an accumulation of inactive Mps1 perturbs accurate chromosome alignment and mitotic progression. Mechanistically, the IRK region binds to the nuclear division cycle 80 complex (Ndc80C), and accumulation of inactive Mps1 at the kinetochores prevents a dynamic interaction between Ndc80C and spindle microtubules (MTs), resulting in an aberrant kinetochore attachment. Thus, our results present a previously undefined mechanism by which Mps1 functions in chromosome alignment by orchestrating Ndc80C–MT interactions and highlight the importance of the precise spatiotemporal regulation of Mps1 kinase activity and kinetochore localization in accurate mitotic progression.Faithful distribution of the duplicated genome into two daughter cells during mitosis depends on proper kinetochore–microtubule (MT) attachments. Defects in kinetochore–MT attachments result in chromosome missegregation, causing aneuploidy, a hallmark of cancer (1, 2). To ensure accurate chromosome segregation, cells use the spindle assembly checkpoint (SAC) to monitor kinetochore biorientation and to control the metaphase-to-anaphase transition. Cells enter anaphase only after the SAC is satisfied, requiring that all kinetochores be attached to MTs and be properly bioriented (3, 4). The core components of SAC signaling include mitotic arrest deficient-like 1 (Mad1), Mad2, Mad3/BubR1 (budding uninhibited by benzimidazole-related 1), Bub1, Bub3, monopolar spindle 1 (Mps1), and aurora B. The full SAC function requires the correct centromere/kinetochore localization of all SAC proteins (5).Among the SAC components, Mps1 was identified originally in budding yeast as a gene required for duplication of the spindle pole body (6). Subsequently, Mps1 orthologs were found in various species, from fungi to mammals. The stringent requirement of Mps1 for SAC activity is conserved in evolution (6–13). Human Mps1 kinase (also known as “TTK”) is expressed in a cell-cycle–dependent manner and has highest expression levels and activity during mitosis. Its localization is also dynamic (8, 14). Although the molecular mechanism remains unclear, Mps1 is required to recruit Mad1 and Mad2 to unattached kinetochores, supporting its essential role in SAC activity (15–18). It also is clear that aurora B kinase activity and the outer-layer kinetochore protein nuclear division cycle 80 (Ndc80)/Hec1 are required for Mps1 localization to kinetochores, as evidenced by recent work, including ours (17, 19–24). How Mps1 activates the SAC is now becoming clear. Mps1 recruits Bub1/Bub3 and BubR1/Bub3 to kinetochores through phosphorylation of KNL1, the kinetochore receptor protein of Bub1 and BubR1 (25–30).Despite much progress in understanding Mps1 functions, it remains unclear how Mps1 is involved in regulating chromosome alignment. In budding yeast mitosis, Mps1 regulates mitotic chromosome alignment by promoting kinetochore biorientation independently of Ipl1 (aurora B in humans) (31), but in budding yeast meiosis Mps1 must collaborate with Ipl1 to mediate meiotic kinetochore biorientation (32). In humans, Mps1 regulates chromosomal alignment by modulating aurora B kinase activity (33), but recent chemical biology studies show that Mps1 kinase activity is important for proper chromosome alignment and segregation, independently of aurora B (22, 34–36). Therefore whether Mps1 regulates chromosome alignment through modulation of aurora B kinase activity is still under debate (37).In this study, we reexamined the function of human Mps1 in chromosome alignment. We found that chromosomal alignment is largely achieved in Mps1 knockdown cells, provided that cells are arrested in metaphase in the presence of MG132, a proteasome inhibitor. However, disrupting Mps1 activity via small molecule inhibitors perturbs chromosomal alignment, even in the presence of MG132. This chromosome misalignment is caused by the abnormal accumulation of inactive Mps1 in the kinetochore and the subsequent failure of correct kinetochore–MT attachments. Further, we demonstrate that inactive Mps1 does not depend on the previously reported tetratricopeptide repeat (TPR) domain for localizing to kinetochores, and we identify a previously unidentified region adjacent to the C terminus of the TPR domain that is responsible for localizing inactive Mps1 to kinetochores. Thus, our work highlights that Mps1 kinase activity is necessary in regulating chromosome alignment and that it must be tightly regulated in space and time to ensure proper localization of Mps1 at kinetochores.     

Modulation of folding energy landscape by charge–charge interactions: Linking experiments with computational modeling

The kinetics of folding–unfolding of a structurally diverse set of four proteins optimized for thermodynamic stability by rational redesign of surface charge–charge interactions is characterized experimentally. The folding rates are faster for designed variants compared with their wild-type proteins, whereas the unfolding rates are largely unaffected. A simple structure-based computational model, which incorporates the Debye–Hückel formalism for the electrostatics, was used and found to qualitatively recapitulate the experimental results. Analysis of the energy landscapes of the designed versus wild-type proteins indicates the differences in refolding rates may be correlated with the degree of frustration of their respective energy landscapes. Our simulations indicate that naturally occurring wild-type proteins have frustrated folding landscapes due to the surface electrostatics. Optimization of the surface electrostatics seems to remove some of that frustration, leading to enhanced formation of native-like contacts in the transition-state ensembles (TSE) and providing a less frustrated energy landscape between the unfolded and TS ensembles. Macroscopically, this results in faster folding rates. Furthermore, analyses of pairwise distances and radii of gyration suggest that the less frustrated energy landscapes for optimized variants are a result of more compact unfolded and TS ensembles. These findings from our modeling demonstrates that this simple model may be used to: (i) gain a detailed understanding of charge–charge interactions and their effects on modulating the energy landscape of protein folding and (ii) qualitatively predict the kinetic behavior of protein surface electrostatic interactions.The energy landscape theory provides a conceptual framework to describe the ensemble nature of the protein folding process (1–3). However, a more detailed understanding of contributions from specific types of interactions remains an active area of research (4, 5). Particularly, the question of how interactions between charged residues modulate the funneled energy landscape is not well explored. These interactions are long-range and thus can alter the conformational ensemble at every step of the folding process. The interactions between charged residues are also nonspecific and either attractive or repulsive and therefore their potential effects on the folding energy landscape can be highly complex (6, 7). Traditionally, the modulation of electrostatic interactions in proteins was done by changing the pH or to a lesser degree changing the ionic strength of the solution (8, 9). Such approaches are complicated by the difficulties of predicting the titration properties of individual amino acid residues in the context of ensembles of protein conformations that are sampled during the folding reaction (10). A more attractive approach is to modulate electrostatic interactions via substitutions that perturb the thermodynamic and kinetic properties of proteins using simple and computationally tractable model systems. Previously, we have shown that the stability of a diverse set of globular proteins can be modulated by rationally redesigning surface charge–charge interactions (11–26). These redesigned proteins are ideally suited to probe the role of electrostatic interactions in modulating the folding energy landscape. The redesigned variants have higher thermodynamic stability than their wild-type proteins. However, because the redesigned proteins contain very few substitutions (less than 5% of total), and because all of the substitutions are on the protein surface, they do not disrupt the native contacts that are important for defining funneled energy landscape (13, 27). Finally, the properties of these proteins can be compared at the same pH, largely eliminating the need to compute titration profiles. In this work, we used four of these redesigned proteins to experimentally probe their folding kinetics and compared them to their corresponding wild-type proteins. The experimental thermodynamic and kinetic data were further rationalized by molecular dynamics simulations using a structure-based model that incorporates the Debye–Hückel formalism to describe interactions between charges. We found that this model qualitatively predicts experimental thermodynamics and kinetics for all four studied proteins and provides insights of how charge–charge interactions modulate the protein folding energy landscape.     

Therapeutic mechanisms of high-frequency stimulation in Parkinson’s disease and neural restoration via loop-based reinforcement

High-frequency deep brain stimulation (HFS) is clinically recognized to treat parkinsonian movement disorders, but its mechanisms remain elusive. Current hypotheses suggest that the therapeutic merit of HFS stems from increasing the regularity of the firing patterns in the basal ganglia (BG). Although this is consistent with experiments in humans and animal models of Parkinsonism, it is unclear how the pattern regularization would originate from HFS. To address this question, we built a computational model of the cortico-BG-thalamo-cortical loop in normal and parkinsonian conditions. We simulated the effects of subthalamic deep brain stimulation both proximally to the stimulation site and distally through orthodromic and antidromic mechanisms for several stimulation frequencies (20–180 Hz) and, correspondingly, we studied the evolution of the firing patterns in the loop. The model closely reproduced experimental evidence for each structure in the loop and showed that neither the proximal effects nor the distal effects individually account for the observed pattern changes, whereas the combined impact of these effects increases with the stimulation frequency and becomes significant for HFS. Perturbations evoked proximally and distally propagate along the loop, rendezvous in the striatum, and, for HFS, positively overlap (reinforcement), thus causing larger poststimulus activation and more regular patterns in striatum. Reinforcement is maximal for the clinically relevant 130-Hz stimulation and restores a more normal activity in the nuclei downstream. These results suggest that reinforcement may be pivotal to achieve pattern regularization and restore the neural activity in the nuclei downstream and may stem from frequency-selective resonant properties of the loop.High-frequency (i.e., above 100 Hz) deep brain stimulation (HFS) of the basal ganglia (BG) and thalamus is clinically recognized to treat movement disorders in Parkinson’s disease (PD) (1–4), but its therapeutic mechanisms remain unclear (5, 6).Early hypotheses about HFS were derived from the rate-based model of the BG function (7, 8) and postulated the disruption of the output of the BG-thalamic system via either the inactivation of neurons in the stimulated site (target) (9–15), which would provide an effect similar to a surgical lesion, or the abnormal excitation of axons projecting out of the target (16–19), which would disrupt the neuronal activity in the structures downstream, including any pathophysiological activity (20).More recently, an ever-growing number of experiments in PD humans and animal models of Parkinsonism has indicated that HFS affects the firing patterns of the neurons rather than the mean firing rate both in the target and the structures downstream (18, 19, 21–31) and it replaces repetitive low-frequency (i.e., ≤50 Hz) bursting patterns with regularized (i.e., more tonic) patterns at higher frequencies (25, 26). It has been proposed that increased pattern regularity of neurons in the target may be therapeutic (5, 32–37), but it is still unknown how this regularity comes about with HFS.It has been suggested that an increased pattern regularity can deplete the information content of the target output and this lack of information would act as an “information lesion” (33) and prevent the pathological activity from being transmitted within the BG-thalamic system (22, 33, 36). As a result, an information lesion in the target [typically, one among the subthalamic nucleus (STN), internal globus pallidus (GPi), or thalamus] would have effects similar to those of a destructive lesion in the same site, which has been reported to alleviate the movement disorders (38).Instead, studies (32, 34, 35, 37) have suggested that an increased pattern regularity of the BG output partly compensates the PD-evoked impairment of the information-processing capabilities of the thalamo-cortical system, and this restores a more faithful thalamic relay of the sensorimotor information (35, 39).Although intriguing, these hypotheses remain elusive on (i) the neuronal mechanisms that would elicit pattern regularization (e.g., why regularization would be relevant only for HFS) and (ii) the effects that increased regularity would have on the cortico-BG-thalamo-cortical loop.It has been hypothesized that pattern regularization occurs because axons projecting out of the target follow the pattern of the stimulus pulses (40, 41) and, given the segregated organization of the BG-thalamic connections (42), it has been assumed that pattern regularization percolates straightforward from the target to the structures immediately downstream (34, 36). However, this representation of the pattern regularization as a “local” effect can hardly be reconciled with the fact that HFS of any structure of the cortico-BG-thalamo-cortical loop is therapeutic for at least some movement disorders (1–4, 43–47), nor does it explain why stimulation at frequencies above 160–180 Hz is not necessarily therapeutic despite the fact that the regularity of the axonal patterns may increase (48, 49). Moreover, coherence in the 8–30-Hz band among neurons across different structures may decrease under HFS but not for lower frequencies (26, 50–52), which suggests the emergence of diffused changes in neuronal activity that would be hardly accounted for with purely local effects.There is emerging evidence, instead, that HFS affects multiple structures simultaneously. First, it has been shown that deep brain stimulation (DBS) may antidromically activate afferent axons and fibers of passage (53–59), thus reaching structures not immediately downstream. Second, studies (57, 58) observed in 6-hydroxydopamine (6-OHDA)-intoxicated rats that the antidromic effects increase with the stimulation frequency and peak around 110–130 Hz. Third, it has been shown in 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-intoxicated nonhuman primates (NHPs) that STN DBS may evoke similar poststimulus responses in different BG structures, both downstream from and upstream to the STN (5, 27, 28, 30, 60). Finally, it has been reported that the cortico-BG-thalamo-cortical system consists of multiple sets of reentrant, interconnected, and partially overlapping neuronal loops (5, 42, 61, 62), which means that the structures upstream to the target (e.g., the striatum) may play an important role in the therapeutic mechanisms of HFS.Altogether, these results suggest that (A) pattern regularization is a global effect that exploits the closed-loop nature of the cortico-BG-thalamo-cortical system and selectively emerges only for specific HFS values, and that (B) the therapeutic merit of pattern regularization has to deal with the restoration of a more normal functionality of the entire cortico-BG-thalamo-cortical loop rather than with variations in the information content of one specific structure.We explored hypotheses (A) and (B) and assessed the system-wide effects of DBS by constructing a computational model of the cortico-BG-thalamo-cortical loop in both normal and parkinsonian conditions and by simulating the effects of STN DBS both at low (20–80 Hz) and high (100–180 Hz) frequencies. The model includes populations of single-compartment neurons and interneurons from motor cortex, striatum, GPi, and thalamus according to a network topology derived from the NHP anatomy, and it simulates both the orthodromic and antidromic effects of DBS. As a result, this model reproduced both average activity and discharge patterns of single units in NHP and rats under normal and parkinsonian conditions, with and without DBS, for all modeled structures.We show through numerical simulation that hypothesis (A) is significantly contributed by reinforcement mechanisms in the striatum. These mechanisms are selectively elicited by HFS, facilitate the percolation of regularized discharge patterns from the striatum to the GPi, and have a primary role in (B), because the percolated striato-pallidal input combines with the local effects of STN DBS to restore the thalamic relay function (63).     

Microscopic identification of the order parameter governing liquid–liquid transition in a molecular liquid

A liquid–liquid transition (LLT) in a single-component substance is an unconventional phase transition from one liquid to another. LLT has recently attracted considerable attention because of its fundamental importance in our understanding of the liquid state. To access the order parameter governing LLT from a microscopic viewpoint, here we follow the structural evolution during the LLT of an organic molecular liquid, triphenyl phosphite (TPP), by time-resolved small- and wide-angle X-ray scattering measurements. We find that locally favored clusters, whose characteristic size is a few nanometers, are spontaneously formed and their number density monotonically increases during LLT. This strongly suggests that the order parameter of LLT is the number density of locally favored structures and of nonconserved nature. We also show that the locally favored structures are distinct from the crystal structure and these two types of orderings compete with each other. Thus, our study not only experimentally identifies the structural order parameter governing LLT, but also may settle a long-standing debate on the nature of the transition in TPP, i.e., whether the transition is LLT or merely microcrystal formation.Liquid-liquid transition (LLT) is an intriguing phenomenon in which a liquid transforms into another one via a first-order transition. This means that there can be more than two liquid states for a single-component substance. Despite its counterintuitive nature, there have recently been many pieces of experimental and numerical evidence for the existence of LLT, for various liquids such as water (1–5), aqueous solutions (6–8), triphenyl phosphite (9–12), l-butanol (13), phosphorus (14), silicon (15, 16), germanium (17), and Y₂O₃–Al₂O₃ (18, 19). This suggests that the LLT may be rather universally observed for various types of liquids. However, none of the LLTs reported so far is free from criticisms (20, 21), mainly because these LLTs take place under experimentally difficult conditions [e.g., at high temperature and pressure (14, 15, 17–19)] or in a supercooled state below the melting point (1–3, 5–7, 9, 10), where the transition is inevitably contaminated by microcrystal formation. The latter is not limited to experiments but arises in numerical simulations, often causing many controversies [LLT (22–25) vs. crystallization (26–28)]. For ST2 water, however, this issue has recently been settled by an extensive simulation study by Palmer et al. (4).One of the hottest and long-standing debates is on the nature of the transition found in a molecular liquid, triphenyl phosphite (TPP), by Kivelson and his coworkers (29). The transition is very easy to access experimentally, because it takes place at ambient pressure and at a temperature range between 230 and 210 K and the transformation speed is slow enough to follow the kinetics. Since the finding of this transition (29, 30), many researchers thus have been interested in this intriguing phenomenon and there have been hot discussions on the nature of the transition (20, 21). Some people interpreted this as a liquid-associated phenomenon (9, 10, 31, 32), but others interpret it differently. All of the controversies come from the fact that this transition accompanies microcrystal formation and thus the final state, which is called “glacial phase,” often contains microcrystallites. This led many researchers to explain the transition by non-LLT scenarios, which include a defect-ordered phase scenario predicted by a frustration limited domain theory (29, 30, 33, 34), a microcrystallization scenario (35–38), and a liquid-crystal or plastic-crystal phase scenario (39). Each scenario captures a certain feature of the glacial phase, but fails in explaining all of the experimental results in a consistent manner. Similar situations are often seen in other candidates of LLTs, such as l-butanol [LLT (13) vs. microcrystallization (40–43)], confined water [LLT (5) vs. other phenomena (44–46)], and aqueous solutions [LLT (6, 7) vs. microcrystallization (8, 28, 47, 48)]. For TPP, however, some pieces of experimental evidence supportive of the LLT scenario rather than the microcrystallization scenario have recently been reported (11, 12).We propose a two-order-parameter (TOP) model of a liquid to explain LLT (20, 49). The main point of this model is that it is necessary to consider the spatiotemporal hierarchical nature of a liquid to understand LLT. More specifically, we argue that in addition to density order parameter ρ describing a gas–liquid transition, we need an additional scalar order parameter S, which is the number density of locally favored structures (LFS). In this model, LLT is a consequence of the cooperative ordering of the scalar nonconserved order parameter S, i.e., the cooperative formation of LFS. In other words, LLT is regarded as a gas–liquid-like transition of LFS: one liquid is a gas state of LFS (low-S state), and the other is its liquid state (high-S state). Recently, it was proposed by Anisimov and coworkers (50, 51) that the thermodynamic ordering field conjugate to the order parameter is the conversion equilibrium constant, which further characterizes the nature of LLT. We explained our experimental observation of LLT in TPP in terms of this model (9, 10). We also studied the phase transition dynamics and the physical and chemical properties of the second liquid state (liquid II), which were also explained by the model (20, 21).However, we have not had any direct experimental evidence for the formation of such LFS up to now; thus, an open question is, what is the relevant order parameter governing LLT, although the link of the order parameter to the enthalpy (9, 10), the refractive index (or, density) (9, 10, 29, 30), and the polarity associated with local molecular ordering (12) has been suggested for LLT in TPP. There have been structural studies on LLT by X-ray and neutron scattering measurements, focusing on local liquid structures at an inter- and intramolecular scale (36, 38, 52–54) and mesoscopic structures (34, 55). However, there has been no experimental evidence for the presence of locally favored structures, which characterize the liquid state uniquely, or the order parameter has still not been identified from a microscopic viewpoint.Here we study the structural change of TPP during LLT by time-resolved small- and wide-angle X-ray scattering measurements, which cover a length scale from a single molecule size ( ～  1 nm) to more than tens of nanometers. We show, to our knowledge, the first direct evidence for the presence of LFS and the temporal increase upon the liquid I-to-liquid II transformation. Furthermore, we also find an indication of the formation of microcrystallites during LLT. However, we reveal that LFS and microcrystallites have different sizes and growth kinetics, indicating that although they sometimes appear simultaneously during the process of LLT, LLT itself is driven by the formation of LFS and not by that of microcrystallites. We also discover that LFS are destroyed upon crystallization, clearly indicating not only that these two types of orderings are competing with each other but also that LFS is a structure unique to the liquid state. Our findings provide a comprehensive view on the long-standing controversy on the origin of the glacial phase, which was discovered by Kivelson and his coworkers (29, 30), and show that the fraction of LFS may be the relevant order parameter of LLT. This suggests that a liquid can have a spatiotemporal hierarchical structure at a low temperature, contrary to the common picture of a high-temperature liquid where the structure is random and homogeneous beyond the molecular size.     

Structural insights into the interaction of IL-33 with its receptors

Interleukin (IL)-33 is an important member of the IL-1 family that has pleiotropic activities in innate and adaptive immune responses in host defense and disease. It signals through its ligand-binding primary receptor ST2 and IL-1 receptor accessory protein (IL-1RAcP), both of which are members of the IL-1 receptor family. To clarify the interaction of IL-33 with its receptors, we determined the crystal structure of IL-33 in complex with the ectodomain of ST2 at a resolution of 3.27 Å. Coupled with structure-based mutagenesis and binding assay, the structural results define the molecular mechanism by which ST2 specifically recognizes IL-33. Structural comparison with other ligand–receptor complexes in the IL-1 family indicates that surface-charge complementarity is critical in determining ligand-binding specificity of IL-1 primary receptors. Combined crystallography and small-angle X-ray–scattering studies reveal that ST2 possesses hinge flexibility between the D3 domain and D1D2 module, whereas IL-1RAcP exhibits a rigid conformation in the unbound state in solution. The molecular flexibility of ST2 provides structural insights into domain-level conformational change of IL-1 primary receptors upon ligand binding, and the rigidity of IL-1RAcP explains its inability to bind ligands directly. The solution architecture of IL-33–ST2–IL-1RAcP complex from small-angle X-ray–scattering analysis resembles IL-1β–IL-1RII–IL-1RAcP and IL-1β–IL-1RI–IL-1RAcP crystal structures. The collective results confer IL-33 structure–function relationships, supporting and extending a general model for ligand–receptor assembly and activation in the IL-1 family.Interleukin (IL)-33 has important roles in initiating a type 2 immune response during infectious, inflammatory, and allergic diseases (1–5). It was initially identified as a nuclear factor in endothelial cells and named NF-HEV (nuclear factor from high endothelial venules) (6, 7). In 2005, it was rediscovered as a new member of the IL-1 family and an extracellular ligand for the orphan IL-1 receptor family member ST2 (8). As an extracellular cytokine, IL-33 is involved in the polarization of Th2 cells and activation of mast cells, basophils, eosinophils, and natural killer cells (1–3). Recent studies also discovered that the type 2 innate lymphoid cells (ILC2s) are major target cells of IL-33 (9, 10). ILC2s express a high level of ST2 and secrete large amounts of Th2 cytokines, most notably IL-5 and IL-13, when stimulated with IL-33 (11–13). Activation of ILC2s is essential in the initiation of the type 2 immune response against helminth infection and during allergic diseases such as asthma (9, 10).IL-33 does not have a signal peptide and is synthesized with an N-terminal propeptide upstream of the IL-1–like cytokine domain. It is preferentially and constitutively expressed in the nuclei of structural and lining cells, particularly in epithelial and endothelial cells (14, 15). Tissue damage caused by pathogen invasion or allergen exposure may lead to the release of IL-33 into extracellular environment from necrotic cells, which functions as an endogenous danger signal or alarmin (14, 16). Full-length human IL-33 consists of 270 residues and is biologically active (17, 18). It is also a substrate of serine proteases released by inflammatory cells recruited to the site of injury (18, 19). The proteases elastase, cathespin G, and proteinase 3 cleave full-length IL-33 to release N-terminal–truncated mature forms containing the IL-1–like cytokine domain: IL-33_95–270, l-33_99–270, and IL-33_109–270 (18). These mature IL-33 forms process a 10-fold greater potency to activate ST2 than full-length IL-33 (18). Caspase-1 was also suggested to cleave IL-33 to generate an active IL-33_112–270 that is the commercially available mature IL-33 form (8). However, it was later demonstrated that this cleavage site does not exist and cleavage by caspases at other sites actually inactivates IL-33 (17, 20, 21).The signaling of IL-33 depends on its binding to the primary receptor ST2 and subsequent recruitment of accessory receptor IL-1RAcP (8, 22, 23). The ligand-binding–induced receptor heterodimerization results in the juxtaposition of the intracellular toll/interleukin-1 receptor (TIR) domains of both receptors, which is necessary and sufficient to activate NF-κB and MAPK pathways in the target cells (24). Previously, we determined the complex structure of IL-1β with its decoy receptor IL-1RII and accessory receptor IL-1RAcP (25). Based on this structure and other previous studies, we proposed a general structural model for the assembly and activation of IL-1 family of cytokines with their receptors (25). In this model, ligand recognition relies on interaction of IL-1 cytokine with its primary receptor: IL-1α and IL-1β with IL-1RI; IL-33 with ST2; IL-18 with IL-18Rα; and IL-36α, IL-36β, and IL-36γ with IL-1Rrp2 (26–28). The binding forms a composite surface to recruit accessory receptor IL-1RAcP shared by IL-1α, IL-1β, IL-33, IL-36α, IL-36β, and IL-36γ and IL-18Rβ by IL-18 (26, 27). This general structural model is further supported by the subsequent structural determination of IL-1β with IL-1RI and IL-1RAcP (29). However, there are still many key missing parts in the general structural model of ligand–receptor interaction in the IL-1 family. For example, the structural basis for specific recognition of IL-33 by ST2 and IL-18 by IL-18Rα, and promiscuous recognition of IL-36α, IL-36β, and IL-36γ by IL-1Rrp2 remains elusive. The proposed general model also needs further confirmation from structural studies of other signaling complexes in the IL-1 family. To address these issues, we studied the interaction of IL-33 with its receptors by a combination of X-ray crystallography and small-angle X-ray–scattering (SAXS) methods.