共查询到20条相似文献,搜索用时 46 毫秒
1.
Joseph Wong David Baddeley Eric?A. Bushong Zeyun Yu Mark?H. Ellisman Masahiko Hoshijima Christian Soeller 《Biophysical journal》2013,104(11):L22-L24
We conducted super-resolution light microscopy (LM) imaging of the distribution of ryanodine receptors (RyRs) and caveolin-3 (CAV3) in mouse ventricular myocytes. Quantitative analysis of data at the surface sarcolemma showed that 4.8% of RyR labeling colocalized with CAV3 whereas 3.5% of CAV3 was in areas with RyR labeling. These values increased to 9.2 and 9.0%, respectively, in the interior of myocytes where CAV3 was widely expressed in the t-system but reduced in regions associated with junctional couplings. Electron microscopic (EM) tomography independently showed only few couplings with caveolae and little evidence for caveolar shapes on the t-system. Unexpectedly, both super-resolution LM and three-dimensional EM data (including serial block-face scanning EM) revealed significant increases in local t-system diameters in many regions associated with junctions. We suggest that this regional specialization helps reduce ionic accumulation and depletion in t-system lumen during excitation-contraction coupling to ensure effective local Ca2+ release. Our data demonstrate that super-resolution LM and volume EM techniques complementarily enhance information on subcellular structure at the nanoscale.The contraction of cardiac ventricular myocytes depends on the rapid cell-wide transient increase in intracellular [Ca2+] upon depolarization of the cell-membrane potential. The cardiac ryanodine receptor (RyR) (1), which is the intracellular Ca2+ release channel in the sarcoplasmic reticulum (SR), plays a central role in shaping Ca2+ transients. RyRs form clusters of various sizes (2,3) with the majority located within junctions between the SR and the surface membrane and its cytoplasmic extension, the transverse tubular (t-) system. It has been suggested that some RyR clusters are associated with caveolae, a specialized signaling microdomain of the surface membrane. Previous studies were complicated by the limited resolution of optical imaging methods of ∼250 nm, much larger than the nanometer scale of RyRs and caveolae. Accordingly, these studies report varying colocalization between RyRs and caveolin-3 (CAV3), a caveolar marker also expressed in the t-system (4,5).In this work, we investigated the relative distribution of CAV3 and RyRs in mouse ventricular myocytes both in the cytosol and near the cell surface with super-resolution fluorescence microscopy that achieves a resolution approaching 30 nm. Our data revealed unexpected local t-system swellings near junctional couplings, which was supported by two different three-dimensional electron microscopy (EM) modalities with <10-nm resolution: EM tomography and serial block-face scanning EM (SBFSEM).Super-resolution images of CAV3 and RyR labeling at the surface sarcolemma of mouse myocytes showed little overlap, suggesting that few RyRs were in couplings with caveolae (Fig. 1
A, for detailed methods, see the Supporting Material). Only ∼4.8% of RyR labeling was associated with CAV3 positive areas and ∼3.5% of CAV3 associated with RyR positive areas (n = 6 cells from three animals, Fig. 1
B, see also Table S1 in the Supporting Material), broadly consistent with previous data in rats (6). To support this finding, EM tomography was applied to mouse ventricular tissue that included a part of the surface sarcolemma, to our knowledge for the first time. Segmentation of peripheral couplings (containing RyR foot structures) and surface caveolae (∼60 nm in diameter and often interconnected) confirmed that the great majority of peripheral couplings were in regions devoid of caveolae (Fig. 1
C). A few junctional couplings containing feet were between caveolae and subsarcolemmal SR (Fig. 1
D, see also Fig. S1 and Movie S1 in the Supporting Material). We conducted a similar analysis in the cytosol where CAV3 expression occurs in the t-system (5) and RyRs are abundant in dyadic junctions between the t-system and SR terminal cisterns.Open in a separate windowFigure 1Colocalization of CAV3 and RyRs at the surface sarcolemma. (A) Super-resolution micrograph of the distribution of CAV3 (green) and RyRs (red) at the surface of a mouse cardiac myocyte. (B) Analysis of the association of CAV3 with RyRs. The fraction of RyR labeling within CAV3 positive areas was ∼4.8% (front data) whereas ∼3.5% of CAV3 was found in RyR-positive membrane areas. (C) Segmented EM tomogram containing a patch of surface sarcolemma (light blue) and associated caveolae (green) as well as peripheral couplings (red). (D) Detailed view of a region with abundant caveolae. (Arrows) Couplings with caveolae.As shown in Fig. 2
A, the spatial distribution of CAV3 and RyR clusters in super-resolution micrographs taken several microns below the surface sarcolemma is consistent with this view. The association of the two labels is slightly increased (as compared to the surface), according to distance analysis with 9% of CAV3 and 9.2% of RyR labeling associating with each other (Fig. 2
B, n = 6 cells from three animals). The similarity of manually traced t-system in EM tomograms (Fig. 2
C) and super-resolved CAV3 labeling suggested that CAV3 is widely distributed in the t-system except for regions where dyadic membrane junctions occur as CAV3 labeling was much weaker in regions with strong RyR labeling. It was notable that the t-system diameter appeared to increase at regions of strong RyR labeling (Fig. 2
D), broadly consistent with the behavior seen in tomograms (Fig. 2
C). This was confirmed by a quantitative analysis of t-tubule diameters in dyadic versus extradyadic regions on the basis of CAV3 and RyR labeling, with full-width at quarter-maximum mean diameters increasing from ∼150 nm distal to dyads, to ∼190 nm (using CAV3 signal only) or ∼280 nm (using CAV3 and RyR signal) near dyads (Fig. 2, G and H, see also Methods in the Supporting Material). The combined RyR and CAV3 signals seemed to be a better representation of the entire t-system lumen near junctions (see Fig. S2).Open in a separate windowFigure 2Distribution of CAV3 and RyRs in the cell interior. (A) Super-resolution micrograph of CAV3 (green) and RyR (red) distribution at t-system. (Arrow) Direction of longitudinal cell axis. (B) Distance analysis of the CAV3 and RyR association (N = 6 cells per group). (C) Segmented EM tomogram of a similar region with three-dimensional mesh models of t-system membrane (green) and dyadic couplings (red). (D) This image illustrates the tracing (white path) of t-tubules. The label distribution was extracted and linearized along the path (E) to calculate a mask that shows the full width at quarter-maximum diameter along tubules, CAV3 (green) and RyR (red) (F). (G) Histograms of local diameters extracted from traced t-tubules. (H) Mean diameters in junctional (dyad) and nonjunctional (ex-dyad) regions. See main text and the Supporting Material for details. **p < 0.01.Taken together, super-resolution imaging and EM tomography strongly support the presence of local t-system dilations in regions where the t-system opposes SR at dyads and such t-system bulges are connected by narrower tubule segments. Further support was provided by SBFSEM, another volume EM technique to study larger cell volumes (albeit at the expense of a slightly lower resolution). SBFSEM clearly showed local t-system dilations were regularly involved in the architecture of most (but not all) dyads (Fig. 3, see also Fig. S3 and Movie S2), as also observed in full three-dimensional super-resolution images (see Fig. S3
C).Open in a separate windowFigure 3Segmented SBFSEM data showing t-system dilations near dyadic junctions. (A) The overview shows t-system membranes (green) and jSR (red) in a mouse myocyte. (B, enlarged inset from panel A) Thin connecting tubules (arrows) and regular swellings in junctional regions at z-lines.Our data identify local dilations of the t-system associated with dyads in mouse cardiac myocytes. Frequent tubule distensions had been observed especially at the intersections of transverse and axial tubules (7), and constrictions were seen in rabbit myocytes although their relationship to dyads was unknown (8). The increased local t-system lumen near junctions may help reduce the predicted ionic accumulation/depletion during excitation-contraction coupling (9). Alternatively, it might simply be secondary to increasing local membrane area and allow the formation of large area junctions that harbor many RyRs. In connection with this point, it would be interesting to investigate the t-system near junctions in species that have larger average tubule diameters (e.g., human and rabbit (10)), or if this architecture changes in mouse heart failure models where t-tubule diameters are often increased.Most peripheral couplings were in regions void of surface caveolae, although a small number of RyR clusters were in junctional couplings between subsarcolemmal SR and caveolae as shown both by the low colocalization between CAV3 and RyRs as well as direct evidence from EM tomography. Similarly, a relatively small fraction of CAV3 colocalized with RyR clusters in the t-system although CAV3 was expressed widely in the t-system. A structural role of CAV3 in the t-system is still unclear—t-tubules in tomogram data did not reveal distinct caveolae shapes on the t-system membrane (see Fig. S4), although this might change in pathology (11). In any case, the t-system exhibits high curvature orthogonal to the tubule axis, which may be supported by CAV3 oligomerization. In addition, the presence of CAV3 in the t-system may be important for regulating other signaling systems (e.g., adrenergic signaling).Finally, our data demonstrate that complementary data from optical super-resolution and three-dimensional EM images assists data interpretation and reliability. We suggest that truly correlative optical and EM imaging approaches should provide further information and improve our knowledge of the basis of cardiac excitation-contraction coupling. 相似文献
2.
Sara Kaliman Christina Jayachandran Florian Rehfeldt Ana-Sun?ana Smith 《Biophysical journal》2014,106(7):L25-L28
It is well established that MDCK II cells grow in circular colonies that densify until contact inhibition takes place. Here, we show that this behavior is only typical for colonies developing on hard substrates and report a new growth phase of MDCK II cells on soft gels. At the onset, the new phase is characterized by small, three-dimensional droplets of cells attached to the substrate. When the contact area between the agglomerate and the substrate becomes sufficiently large, a very dense monolayer nucleates in the center of the colony. This monolayer, surrounded by a belt of three-dimensionally packed cells, has a well-defined structure, independent of time and cluster size, as well as a density that is twice the steady-state density found on hard substrates. To release stress in such dense packing, extrusions of viable cells take place several days after seeding. The extruded cells create second-generation clusters, as evidenced by an archipelago of aggregates found in a vicinity of mother colonies, which points to a mechanically regulated migratory behavior.Studying the growth of cell colonies is an important step in the understanding of processes involving coordinated cell behavior such as tissue development, wound healing, and cancer progression. Apart from extremely challenging in vivo studies, artificial tissue models are proven to be very useful in determining the main physical factors that affect the cooperativity of cells, simply because the conditions of growth can be very well controlled. One of the most established cell types in this field of research is the Madin-Darby canine kidney epithelial cell (MDCK), originating from the kidney distal tube (1). A great advantage of this polarized epithelial cell line is that it retained the ability for contact inhibition (2), which makes it a perfect model system for studies of epithelial morphogenesis.Organization of MDCK cells in colonies have been studied in a number of circumstances. For example, it was shown that in three-dimensional soft Matrigel, MDCK cells form a spherical enclosure of a lumen that is enfolded by one layer of polarized cells with an apical membrane exposed to the lumen side (3). These structures can be altered by introducing the hepatocyte growth factor, which induces the formation of linear tubes (4). However, the best-studied regime of growth is performed on two-dimensional surfaces where MDCK II cells form sheets and exhibit contact inhibition. Consequently, the obtained monolayers are well characterized in context of development (5), mechanical properties (6), and obstructed cell migration (7–9).Surprisingly, in the context of mechanics, several studies of monolayer formation showed that different rigidities of polydimethylsiloxane gels (5) and polyacrylamide (PA) gels (9) do not influence the nature of monolayer formation nor the attainable steady-state density. This is supposedly due to long-range forces between cells transmitted by the underlying elastic substrate (9). These results were found to agree well with earlier works on bovine aortic endothelial cells (10) and vascular smooth muscle cells (11), both reporting a lack of sensitivity of monolayers to substrate elasticity. Yet, these results are in stark contrast with single-cell experiments (12–15) that show a clear response of cell morphology, focal adhesions, and cytoskeleton organization to substrate elasticity. Furthermore, sensitivity to the presence of growth factors that are dependent on the elasticity of the substrate in two (16) and three dimensions (4) makes this result even more astonishing. Therefore, we readdress the issue of sensitivity of tissues to the elasticity of the underlying substrate and show that sufficiently soft gels induce a clearly different tissue organization.We plated MDCK II cells on soft PA gels (Young’s modulus E = 0.6 ± 0.2 kPa), harder PA gels (E = 5, 11, 20, 34 kPa), and glass, all coated with Collagen-I. Gels were prepared following the procedure described in Rehfeldt et al. (17); rigidity and homogeneity of the gels was confirmed by bulk and microrheology (see the Supporting Material for comparison). Seeding of MDCK II cells involved a highly concentrated solution dropped in the middle of a hydrated gel or glass sample. For single-cell experiments, cells were dispersed over the entire dish. Samples were periodically fixed up to Day 12, stained for nuclei and actin, and imaged with an epifluorescence microscope. Details are described in the Supporting Material.On hard substrates and glass it was found previously that the area of small clusters expands exponentially until the movement of the edge cannot keep up with the proliferation in the bulk (5). Consequently, the bulk density increases toward the steady state, whereas the density of the edge remains low. At the same time, the colony size grows subexponentially (5). This is what we denote “the classical regime of growth”. Our experiments support these observations for substrates with E ≥ 5 kPa. Specifically, on glass, colonies start as small clusters of very low density of 700 ± 200 cells/mm2 (Fig. 1, A and B), typically surrounded by a strong actin cable (Fig. 1, B and C). Interestingly, the spreading area of single cells (Fig. 1
A) on glass was found to be significantly larger, i.e., (2.0 ± 0.9) × 10−3 mm2. After Day 4 (corresponding cluster area of 600 ± 100 mm2), the density in the center of the colony reached the steady state with 6,800 ± 500 cells/mm2, whereas the mean density of the edge profile grew to 4,000 ± 500 cells/mm2. This density was retained until Day 12 (cluster area 1800 ± 100 mm2), which is in agreement with previous work (9).Open in a separate windowFigure 1Early phase of cluster growth on hard substrates. (A) Well-spread single cells, and small clusters with a visible actin cable 6 h after seeding. (B) Within one day, clusters densify and merge, making small colonies. (C) Edge of clusters from panel B.In colonies grown on 0.6 kPa gels, however, we encounter a very different growth scenario. The average spreading area of single cells is (0.34 ± 0.3) × 10−3 mm2, which is six times smaller than on glass substrates (Fig. 2
A). Clusters of only few cells show that cells have a preference for cell-cell contacts (a well-established flat contact zone can be seen at the cell-cell interface in Fig. 2
A) rather than for cell-substrate contacts (contact zone is diffusive and the shape of the cells appears curved). The same conclusion emerges from the fact that dropletlike agglomerates, resting on the substrate, form spontaneously (Fig. 2
A), and that attempts to seed one single cluster of 90,000 cells fail, resulting in a number of three-dimensional colonies (Fig. 2
A). When the contact area with the substrate exceeds 4.7 × 10−3 mm2, a monolayer appears in the center of such colonies (Fig. 2
B). The colonies can merge, and if individual colonies are small, the collapse into a single domain is associated with the formation of transient irregular structures (Fig. 2
B). Ultimately, large elliptical colonies (average major/minor axis of e = 1.8 ± 0.6) with a smooth edge are formed (Fig. 2
C), unlike on hard substrates where circular clusters (e = 1.06 ± 0.06) with a ragged edge comprise the characteristic phenotype.Open in a separate windowFigure 2Early phase of cluster growth on soft substrates. (A) Twelve hours after seeding, single cells remain mostly round and small. They are found as individual, or within small, three-dimensional structures (top). The latter nucleate a monolayer in their center (bottom), if the contact area with the substrate exceeds ∼5 × 10−3 mm2. (B) Irregularly-shaped clusters appear due to merging of smaller droplets. A stable monolayer surrounded by a three-dimensional belt of densely packed cells is clearly visible, even in larger structures. (C) All colonies are recorded on Day 4.Irrespective of cluster size, in the new regime of growth, the internal structure is built of two compartments (Fig. 2
B):
- 1.The first is the edge (0.019 ± 0.05-mm wide), a three-dimensional structure of densely packed cells. This belt is a signature of the new regime because on hard substrates the edge is strictly two-dimensional (Fig. 1 C).
- 2.The other is the centrally placed monolayer with a spatially constant density that is very weakly dependent on cluster size and age (Fig. 3). The mean monolayer density is 13,000 ± 2,000 cells/mm2, which is an average over 130 clusters that are up to 12 days old and have a size in the range of 10−3 to 10 mm2, each shown by a data point in Fig. 3. This density is twice the steady-state density of the bulk tissue in the classical regime of growth.Open in a separate windowFigure 3Monolayer densities in colonies grown on 0.6 kPa substrates, as a function of the cluster size and age. Each cluster is represented by a single data point signifying its mean monolayer density. (Black lines) Bulk and (red dashed lines) edge of steady-state densities from monolayers grown on glass substrates. Error bars are omitted for clarity, but are discussed in the Supporting Material.
3.
FtsZ, a bacterial homolog of eukaryotic tubulin, assembles into the Z ring required for cytokinesis. In Escherichia coli, FtsZ interacts directly with FtsA and ZipA, which tether the Z ring to the membrane. We used three-dimensional structured illumination microscopy to compare the localization patterns of FtsZ, FtsA, and ZipA at high resolution in Escherichia coli cells. We found that FtsZ localizes in patches within a ring structure, similar to the pattern observed in other species, and discovered that FtsA and ZipA mostly colocalize in similar patches. Finally, we observed similar punctate and short polymeric structures of FtsZ distributed throughout the cell after Z rings were disassembled, either as a consequence of normal cytokinesis or upon induction of an endogenous cell division inhibitor.The assembly of the bacterial tubulin FtsZ has been well studied in vitro, but the fine structure of the cytokinetic Z ring it forms in vivo is not well defined. Super-resolution microscopy methods including photoactivated localization microscopy (PALM) and three-dimensional-structured illumination microscopy (3D-SIM) have recently provided a more detailed view of Z-ring structures. Two-dimensional PALM showed that Z rings in Escherichia coli are likely composed of loosely-bundled dynamic protofilaments (1,2). Three-dimensional PALM studies of Caulobacter crescentus initially showed that Z rings were comprised of loosely bundled protofilaments forming a continuous but dynamic ring (1–3). However, a more recent high-throughput study showed that the Z rings of this bacterium are patchy or discontinuous (4), similar to Z rings of Bacillus subtilis and Staphylococcus aureus using 3D-SIM (5). Strauss et al. (5) also demonstrated that the patches in B. subtilis Z rings are highly dynamic.Assembly of the Z ring is modulated by several proteins that interact directly with FtsZ and enhance assembly or disassembly (6). For example, FtsA and ZipA promote ring assembly in E. coli by tethering it to the cytoplasmic membrane (7,8). SulA is an inhibitor of FtsZ assembly, induced only after DNA damage, which sequesters monomers of FtsZ to prevent its assembly into a Z ring (9). Our initial goals were to visualize Z rings in E. coli using 3D-SIM, and then examine whether any FtsZ polymeric structures remain after SulA induction. We also asked whether FtsA and ZipA localized in patchy patterns similar to those of FtsZ.We used a DeltaVision OMX V4 Blaze microscope (Applied Precision, GE Healthcare, Issaquah, WA) to view the high-resolution localization patterns of FtsZ in E. coli cells producing FtsZ-GFP (Fig. 1). Three-dimensional views were reconstructed using softWoRx software (Applied Precision). To rule out GFP artifacts, we also visualized native FtsZ from a wild-type strain (WM1074) by immunofluorescence (IF).Open in a separate windowFigure 1Localization of FtsZ in E. coli. (A) Cell with a Z ring labeled with FtsZ-GFP. (B) Rotated view of Z ring in panel A. (C) Cell with a Z ring labeled with DyLight 550 (Thermo Fisher Scientific, Waltham, MA). (D) Rotated view of Z ring in panel C. (B1 and D1) Three-dimensional surface intensity plots of Z rings in panels B and D, respectively. (E) A dividing cell producing FtsZ-GFP. The cell outline is shown in the schematic. (Asterisk) Focus of FtsZ localization; (open dashed ovals) filamentous structures of FtsZ. Three-dimensional surface intensity plots were created using the software ImageJ (19). Scale bars, 1 μm.Both FtsZ-GFP (Fig. 1, A, B, and B1) and IF staining for FtsZ (Fig. 1, C, D, and D1) consistently localized to patches around the ring circumference, similar to the B. subtilis and C. crescentus FtsZ patterns (4,5). Analysis of fluorescence intensities (see Fig. S1, A and B, in the Supporting Material) revealed that the majority of Z rings contain one or more gaps in which intensity decreases to background levels (82% for FtsZ-GFP and 69% for IF). Most rings had 3–5 areas of lower intensity, but only a small percentage of these areas had fluorescence below background intensity (34% for FtsZ-GFP and 21% for IF), indicating that the majority of areas with lower intensity contain at least some FtsZ.To elucidate how FtsZ transitions from a disassembled ring to a new ring, we imaged a few dividing daughter cells before they were able to form new Z rings (Fig. 1
E). Previous conventional microscopy had revealed dynamic FtsZ helical structures (10), but the resolution had been insufficient to see further details. Here, FtsZ visualized in dividing cells by 3D-SIM localized throughout as a mixture of patches and randomly-oriented short filaments (asterisk and dashed oval in Fig. 1, respectively). These structures may represent oligomeric precursors of Z ring assembly.To visualize FtsZ after Z-ring disassembly another way, we overproduced SulA, a protein that blocks FtsZ assembly. We examined E. coli cells producing FtsZ-GFP after induction of sulA expression from a pBAD33-sulA plasmid (pWM1736) with 0.2% arabinose. After 30 min of sulA induction, Z rings remained intact in most cells (Fig. 2
A and data not shown). The proportion of cellular FtsZ-GFP in the ring before and after induction of sulA was consistent with previous data (data not shown) (1,11).Open in a separate windowFigure 2Localization of FtsZ after overproduction of SulA. (A) Cell producing FtsZ-GFP after 0.2% arabinose induction of SulA for 30 min. (B) After 45 min. (B1) Magnified cell shown in panel B. (C) Cell producing native FtsZ labeled with AlexaFluor 488 (Life Technologies, Carlsbad, CA) 30 min after induction; (D) 45 min after induction. (D1) Magnified cell shown in panel D. Scale bars, 1 μm. (Asterisk) Focus of FtsZ localization; (open dashed ovals) filamentous structures of FtsZ.Notably, after 45 min of sulA induction, Z rings were gone (Fig. 2, B and B1), replaced by numerous patches and randomly-oriented short filaments (asterisk and dashed ovals in Fig. 2), similar to those observed in a dividing cell. FtsZ normally rapidly recycles from free monomers to ring-bound polymers (11), but a critical concentration of SulA reduces the pool of available FtsZ monomers, resulting in breakdown of the Z ring (9). The observed FtsZ-GFP patches and filaments are likely FtsZ polymers that disassemble before they can organize into a ring.We confirmed this result by overproducing SulA in wild-type cells and detecting FtsZ localization by IF (Fig. 2, C, D, and D1). The overall fluorescence patterns in cells producing FtsZ-GFP versus cells producing only native FtsZ were similar (Fig. 2, B1 and D1), although we observed fewer filaments with IF, perhaps because FtsZ-GFP confers slight resistance to SulA, or because the increased amount of FtsZ in FtsZ-GFP producing cells might titrate the SulA more effectively.Additionally, we wanted to observe the localization patterns of the membrane tethers FtsA and ZipA. Inasmuch as both proteins bind to the same C-terminal conserved tail of FtsZ (12–14), they would be expected to colocalize with the circumferential FtsZ patches in the Z ring. We visualized FtsA using protein fusions to mCherry and GFP (data not shown) as well as IF using a wild-type strain (WM1074) (Fig. 3
A). We found that the patchy ring pattern of FtsA localization was similar to the FtsZ pattern. ZipA also displayed a similar patchy localization in WM1074 by IF (Fig. 3
B).Open in a separate windowFigure 3Localization of FtsA (A) and ZipA (B) by IF using AlexaFluor 488. (C) FtsA-GFP ring. (D) Same cell shown in panel C with ZipA labeled with DyLight 550. (C1 and D1) Three-dimensional surface intensity plots of FtsA ring from panel C or ZipA ring from panel D, respectively. (E) Merged image of FtsA (green) and ZipA (red) from the ring shown in panels C and D. (F) Intensity plot of FtsA (green) and ZipA (red) of ring shown in panel E. The plot represents intensity across a line drawn counterclockwise from the top of the ring around the circumference, then into its lumen. Red/green intensity plot and three-dimensional surface intensity plots were created using the software ImageJ (19). Scale bar, 1 μm.To determine whether FtsA and ZipA colocalized to these patches, we used a strain producing FtsA-GFP (WM4679) for IF staining of ZipA using a red secondary antibody. FtsA-GFP (Fig. 3
C) and ZipA (Fig. 3
D) had similar patterns of fluorescence, although the three-dimensional intensity profiles (Fig. 3, C1 and D1) reveal slight differences in intensity that are also visible in a merged image (Fig. 3
E). Quantitation of fluorescence intensities around the circumference of the rings revealed that FtsA and ZipA colocalized almost completely in approximately half of the rings analyzed (Fig. 3
F, and see Fig. S2
A), whereas in the other rings there were significant differences in localization in one or more areas (see Fig. S2
B). FtsA and ZipA bind to the same C-terminal peptide of FtsZ and may compete for binding. Cooperative self-assembly of FtsA or ZipA might result in large-scale differential localization visible by 3D-SIM.In conclusion, our 3D-SIM analysis shows that the patchy localization of FtsZ is conserved in E. coli and suggests that it may be widespread among bacteria. After disassembly of the Z ring either in dividing cells or by excess levels of the cell division inhibitor SulA, FtsZ persisted as patches and short filamentous structures. This is consistent with a highly dynamic population of FtsZ monomers and oligomers outside the ring, originally observed as mobile helices in E. coli by conventional fluorescence microscopy (10) and by photoactivation single-molecule tracking (15). FtsA and ZipA, which bind to the same segment of FtsZ and tether it to the cytoplasmic membrane, usually display a similar localization pattern to FtsZ and each other, although in addition to the differences we detect by 3D-SIM, there are also likely differences that are beyond its ∼100-nm resolution limit in the X,Y plane.As proposed previously (16), gaps between FtsZ patches may be needed to accommodate a switch from a sparse Z ring to a more condensed ring, which would provide force to drive ring constriction (17). If this model is correct, the gaps should close upon ring constriction, although this may be beyond the resolution of 3D-SIM in constricted rings. Another role for patches could be to force molecular crowding of low-abundance septum synthesis proteins such as FtsI, which depend on FtsZ/FtsA/ZipA for their recruitment, into a few mobile supercomplexes.How are FtsZ polymers organized within the Z-ring patches? Recent polarized fluorescence data suggest that FtsZ polymers are oriented both axially and circumferentially within the Z ring in E. coli (18). The seemingly random orientation of the non-ring FtsZ polymeric structures we observe here supports the idea that there is no strong constraint requiring FtsZ oligomers to follow a circumferential path around the cell cylinder. The patches of FtsZ in the unperturbed E. coli Z ring likely represent randomly oriented clusters of FtsZ filaments that are associated with ZipA, FtsA, and essential septum synthesis proteins. New super-resolution microscopy methods should continue to shed light on the in vivo organization of these protein assemblies. 相似文献
4.
Andrew?K. Lewis Zachary?M. James Jesse?E. McCaffrey Anthony?R. Braun Christine?B. Karim David?D. Thomas Jonathan?N. Sachs 《Biophysical journal》2014,106(6):L21-L24
It has long been presumed that activation of the apoptosis-initiating Death Receptor 5, as well as other structurally homologous members of the TNF-receptor superfamily, relies on ligand-stabilized trimerization of noninteracting receptor monomers. We and others have proposed an alternate model in which the TNF-receptor dimer—sitting at the vertices of a large supramolecular receptor network of ligand-bound receptor trimers—undergoes a closed-to-open transition, propagated through a scissorslike conformational change in a tightly bundled transmembrane (TM) domain dimer. Here we have combined electron paramagnetic resonance spectroscopy and potential-of-mean force calculations on the isolated TM domain of the long isoform of DR5. The experiments and calculations both independently validate that the opening transition is intrinsic to the physical character of the TM domain dimer, with a significant energy barrier separating the open and closed states.Death receptor 5 (DR5) is a member of the tumor necrosis factor receptor (TNFR) superfamily that mediates apoptosis when bound by its cognate ligand, TNF-related apoptosis-inducing ligand (1). Upregulated in cancer cells, DR5 is among the most actively pursued anticancer targets (2). TNF-related apoptosis-inducing ligand binds to preassembled DR5 trimers at their extracellular domains, causing the formation of oligomeric ligand-receptor networks that are held together by receptor dimers (3). In the long-isoform of DR5, this dimer is crosslinked via ligand-induced disulfide bond formation between two transmembrane (TM) domain α-helices at Cys-209, and is further stabilized by a GxxxG motif one helix-turn downstream (3).Our recent study of the structurally homologous TNFR1 showed that receptor activation involves a conformational change that propagates from the extracellular domain to the cytosolic domain through a separation (or opening) of the TM domains of the dimer (4). We have therefore hypothesized that the activation of DR5, and indeed all structurally homologous TNF-receptors, involves a scissorslike opening of the TM domain dimer (Fig. 1).Open in a separate windowFigure 1Activation model of the DR5-L TM dimer. The sequence and positions of the disulfide bond and TOAC spin label (top), along with our previously published model (bottom, left) are shown. We propose an activation model (bottom, right) in which the transmembrane dimer pivots at its disulfide bond to reach an active open conformation.Using electron paramagnetic resonance (EPR) spectroscopy, a technique that has been used previously to study TM helix architecture and dynamics (5,6), and potential-of-mean force (PMF) calculations (7,8), this study addresses the question of whether the isolated disulfide-linked DR5-L TM domain dimer occupies distinct open and closed states (Fig. 1), and how its dynamic behavior contributes to the free-energy landscape of the opening transition of the full-length receptor.The DR5-L TM domain was synthesized with TOAC, an amino acid with a nitroxide spin label rigidly fixed to the α-carbon (9), incorporated at position 32 (Fig. 1), with some minor modification to facilitate EPR measurements. Previous work confirmed that this peptide forms disulfide-linked dimers (e.g., via comparison to 2-ME treated sample) and a negligible population of higher-order oligomers (further supported by model fitting of the EPR data below). For peptide work, residues were renumbered such that Thr-204 corresponds to Thr-1, and so on. The cytosolic Cys-29 (which we previously showed does not participate in a disulfide bond in cells) was replaced with serine to prevent the formation of antiparallel disulfide-linked dimers, and Trp-34 was replaced with tyrosine to prevent intrinsic fluorescence in fluorescence studies (not published). Continuous-wave (CW) dipolar EPR (sensitive only to spin-spin distances <25 Å) was used to measure TOAC-TOAC distances within the TM dimers and revealed an ordered Gaussian distribution centered at 16 Å (full width half-maximum (FWHM) = 4 Å), corresponding to a closed state (Fig. 2
A). Double electron-electron resonance (DEER) (sensitive to spin-spin distances from 15 to 60 Å) also detected a short distance consistent with the dipolar EPR data, along with a longer, disordered component (32.9 Å, FWHM = 28 Å) (Fig. 2
B). Together, these measurements indicate the presence of a compact, ordered closed state and a broader, disordered open state. EPR on oriented membranes also indicated two structural states. Global fitting revealed two populations of spin-label tilt angles (orientation of the nitroxide principal axis relative to the membrane normal): a narrow conformation (24°, FWHM = 20°), and a disordered conformation (50°, FWHM = 48°) (Fig. 2
C). This bimodal orientational distribution (Fig. 2
C) is remarkably consistent with the bimodal distance distribution (Fig. 2
B).Open in a separate windowFigure 2EPR spectra (left) of 32-TOAC-DR5 in lipid, and resulting structural distributions (right). (A) CW dipolar EPR spectra (left) of dimer (1 mM diamide) and monomer (1 mM 2-mercaptoethanol). Best-fit spin-spin distance distribution was a single Gaussian centered at 16 ± 2 Å (right). (B) The DEER waveform (left) of 32-TOAC-DR5 dimer was best fit (right) to a two-Gaussian distribution. The short distance was constrained to agree with the CW data, because DEER has poor sensitivity for distances <20 Å. The long-distance distribution is centered at 32.9 Å and is much broader. (C) CW EPR spectra (left) of 32-TOAC-DR5, with the membrane-normal oriented parallel (red) and perpendicular (blue) to the field. Simultaneous (global) fitting of these spectra reveals narrow and broad components (right). (In panels B and C, the overall distribution is plotted as black, while the closed and open components are plotted as green and magenta, respectively.)We subsequently conducted a PMF calculation (10) using the DR5-L TM dimer starting configuration developed by our group previously (3), embedded in a DMPC bilayer, with the Leu-32/Leu-32 Cα distance as the reaction coordinate. Three calculations were run from independent starting configurations, each using 50 windows spaced in 0.5° increments, and run for 20 ns at each window (totaling 3 μs). Each of the calculations yielded a similar result, and the averaged free energy curve (Fig. 3
A) agrees remarkably well with our EPR measurements: a narrow distribution at the closed conformation (∼16 Å, Fig. 3
B) separated by an ∼3 kcal/mol energy barrier from a broad distribution of accessible open conformations at ∼27 Å, (Fig. 3
C). Each of the three individual PMF plots can be found in Fig. S1 in the Supporting Material.Open in a separate windowFigure 3(A) PMF calculation of the DR5 TM domain dimer along the Leu-32/Leu-32 distance reaction coordinate. The PMF calculation reveals a narrow closed state and a broader open state separated by a free energy barrier. Representative snapshots of the (B) closed state and (C) open state.In the closed state, the helices are tightly packed at the GxxxG interfacial motif and all the way down the juxtaposed helix faces at residues Ala-18, Leu-22, Ala-25, and Val-26. The tight packing is aided by kinking and twisting of the two helices around their common axis, increasing the interacting surface area. In the open conformations, the Ala-18, Leu-22, Ala-25, and Val-26 pairs are dissociated and, interestingly, the GxxxG motif at Gly-10 and Gly-14 remains tightly packed. The open state energy well is only slightly less favorable than the closed state (by ∼2 kcal/mol), and its free energy profile is relatively broad and flat. The increased crossing angle in the open state is facilitated by straightening of the helix kink and is not accommodated by a change in bilayer thickness (see Fig. S3, A and B).The observed change in helix-helix distance (11 Å between the two minima in the PMF) is extremely close to that observed previously in live-cell FRET studies of a constitutively active form of TNFR1 (∼8 Å change between states using large fluorescence probes at the cytosolic domains) (4). The change observed in the EPR data (17 Å) may be an overestimate because the measurement is made between TOAC spin labels that likely protrude from the two helices, depending on rotational orientation. These results collectively show that activation of these receptors requires a small, but clearly significant conformational opening of the TM domains. One important note is that our EPR experiments recapitulate the equilibrium distribution of the two states despite there being no driving force to traverse the barrier between them (∼3 kcal/mol in the closed-to-open transition and ∼1 kcal/mol in the open-to-closed transition, Fig. 3). We do not interpret the results to mean that the dimer necessarily traverses these barriers at 4°C. Rather, there likely exist multiple reaction paths for dimerization of the abstracted TM domains. Finally, in the context of the full-length receptor, how the ligand induces a conformational change capable of overcoming the closed-to-open barrier remains an important question.Whether the observed structural transition in the TM domain dimer of the long-isoform of DR5 is a ubiquitous conformational switch that acts over the entire TNFR superfamily remains unknown. Vilar et al. (11) first proposed a similar scissors-model for activation of p75 neurotrophin receptor, which has a cysteine at the center of its TM helix. The short isoform of DR5 lacks a TM domain cysteine, but does form noncovalent dimers in cells, with likely TM domain dimer contacts (3). Among the other closely related and structurally homologous members of the TNFR superfamily, TNFR1 contains a cysteine at the center of the TM domain, but lacks any discernible small residue motifs (e.g., GxxxG). TNFR2 lacks a TM cysteine on the extracellular side, but does have a GxxxG motif positioned similarly to that of DR5. On the other hand, Death Receptor 4, whose functional distinction from DR5 has remained somewhat elusive, lacks both a cysteine and any recognizable small-residue hydrophobic motif.In summary, we have extended recent findings that point to the TM domain of DR5 as an essential structural component in the conformational change associated with activation. Our findings that the DR5-L TM domain occupies distinct open and closed states, separated by a substantial energy barrier, points the way to further studies across the TNF-receptor superfamily. 相似文献
5.
6.
Magnus Andersson Jakob?P. Ulmschneider Martin?B. Ulmschneider Stephen?H. White 《Biophysical journal》2013,104(6):L12-L14
The distribution of peptide conformations in the membrane interface is central to partitioning energetics. Molecular-dynamics simulations enable characterization of in-membrane structural dynamics. Here, we describe melittin partitioning into dioleoylphosphatidylcholine lipids using CHARMM and OPLS force fields. Although the OPLS simulation failed to reproduce experimental results, the CHARMM simulation reported was consistent with experiments. The CHARMM simulation showed melittin to be represented by a narrow distribution of folding states in the membrane interface.Unstructured peptides fold into the membrane interface because partitioned hydrogen-bonded peptide bonds are energetically favorable compared to free peptide bonds (1–3). This folding process is central to the mechanisms of antimicrobial and cell-penetrating peptides, as well as to lipid interactions and stabilities of larger membrane proteins (4). The energetics of peptide partitioning into membrane interfaces can be described by a thermodynamic cycle (Fig. 1). State A is a theoretical state representing the fully unfolded peptide in water, B is the unfolded peptide in the membrane interface, C is the peptide in water, and D is the folded peptide in the membrane. The population of peptides in solution (State C) is best described as an ensemble of folded and unfolded conformations, whereas the population of peptides in State D generally is assumed to have a single, well-defined helicity, as shown in Fig. 1
A (5). Given that, in principle, folding in solution and in the membrane interface should follow the same basic rules, peptides in state D could reasonably be assumed to also be an ensemble. A fundamental question (5) is therefore whether peptides in state D can be correctly described as having a single helicity. Because differentiating an ensemble of conformations and a single conformation may be an impossible experimental task (5), molecular-dynamics (MD) simulations provide a unique high-resolution view of the phenomenon.Open in a separate windowFigure 1Thermodynamic cycles for peptide partitioning into a membrane interface. States A and B correspond to the fully unfolded peptide in solution and membrane interface, respectively. The folded peptide in solution is best described as an ensemble of unfolded and folded conformations (State C). State D is generally assumed to be one of peptides with a narrow range of conformations, but the state could actually be an ensemble of states as in the case of State C.Melittin is a 26-residue, amphipathic peptide that partitions strongly into membrane interfaces and therefore has become a model system for describing folding energetics (3,6–8). Here, we describe the structural dynamics of melittin in a dioleoylphosphatidylcholine (DOPC) bilayer by means of two extensive MD simulations using two different force fields.We extended a 12-ns equilibrated melittin-DOPC system (9) by 17 μs using the Anton specialized hardware (10) with the CHARMM22/36 protein/lipid force field and CMAP correction (11,12) (see Fig. S1 and Fig. S2 in the Supporting Material). To explore force-field effects, a similar system was simulated for 2 μs using the OPLS force field (13) (see Methods in the Supporting Material). In agreement with x-ray diffraction measurements on melittin in DOPC multilayers (14), melittin partitioned spontaneously into the lipid headgroups at a position below the phosphate groups at similar depth as glycerol/carbonyl groups (Fig. 2).Open in a separate windowFigure 2Melittin partitioned into the polar headgroup region of the lipid bilayer. (A) Snapshot of the simulation cell showing two melittin molecules (MLT1 and MLT2, in yellow) at the lipid-water interface. (B) Density cross-section of the simulation cell extracted from the 17-μs simulation. The peptides are typically located below the lipid phosphate (PO4) groups, in a similar depth as the glycerol/carbonyl (G/C) groups.To describe the secondary structure for each residue, we defined helicity by backbone dihedral angles (φ, ψ) within 30° from the ideal α-helical values (–57°, –47°). The per-residue helicity in the CHARMM simulation displays excellent agreement with amide exchange rates from NMR measurements that show a proline residue to separate two helical segments, which are unfolded below Ala5 and above Arg22 (15) (Fig. 3
A). In contrast, the OPLS simulation failed to reproduce the per-residue helicity except for a short central segment (see Fig. S3).Open in a separate windowFigure 3Helicity and conformational distribution of melittin as determined via MD simulation. (A) Helicity per residue for MLT1 and MLT2. (B) Corresponding evolution of the helicity. (C) Conformational distributions over the entire 17-μs simulation.Circular dichroism experiments typically report an average helicity of ∼70% for melittin at membrane interfaces (3,6,16,17), but other methods yield average helicities as high as 85% (15,18). Our CHARMM simulations are generally consistent with the experimental results, especially amide-exchange measurements (15); melittin helicity averaged to 78% for MLT1, whereas MLT2 transitioned from 75% to 89% helicity at t ≈ 8 μs, with an overall average helicity of 82% (Fig. 3
B). However, in the OPLS simulation, melittin steadily unfolds over the first 1.3 μs, after which the peptide remains only partly folded, with an average helicity of 33% (see Fig. S3). Similar force-field-related differences in peptide helicity were recently reported, albeit at shorter timescales (19). Although suitable NMR data are not presently available, we have computed NMR quadrupolar splittings for future reference (see Fig. S4).To answer the question asked in this article—whether the conformational space of folded melittin in the membrane interface can be described by a narrow distribution—the helicity distributions for the equilibrated trajectories are shown in Fig. 3
C. Whereas MLT1 in the CHARMM simulation produces a single, narrow distribution of the helicity, MLT2 has a bimodal distribution as a consequence of the folding event at t ≈ 8 μs (Fig. 3
C). We note that CHARMM force fields have a propensity for helix-formation and this transition might therefore be an artifact. We performed a cluster analysis to describe the structure of the peptide in the membrane interface. The four most populated conformations in the CHARMM simulation are shown in Fig. 4. The dominant conformation for both peptides was a helix kinked at G12 and unfolded at the last 5–6 residues of the C-terminus. The folding transition of MLT2 into a complete helix is visible by the 48% occupancy of a fully folded helix.Open in a separate windowFigure 4Conformational clusters of the two melittin peptides (MLT1 and MLT2) from the 17-μs CHARMM simulation in DOPC. Clustering is based on Cα-RMSD with a cutoff criterion of 2 Å.We conclude that the general assumption when calculating folding energetics holds: Folded melittin partitioned into membrane interfaces can be described by a narrow distribution of conformations. Furthermore, extended (several microsecond) simulations are needed to differentiate force-field effects. Although the CHARMM and OPLS simulations would seem to agree for the first few hundred nanoseconds, the structural conclusions differ drastically with longer trajectories, with CHARMM parameters being more consistent with experiments. However, as implied by the difference in substate distributions between MLT1 and MLT2, 17 μs might not be sufficient to observe the fully equilibrated partitioning process. The abrupt change in MLT2 might indicate that the helicity will increase to greater than experimentally observed in a sufficiently long simulation. On the other hand, it could be nothing more than a transient fluctuation. Increased sampling will provide further indicators of convergence of the helix partitioning process. 相似文献
7.
The T-cell actin cytoskeleton mediates adaptive immune system responses to peptide antigens by physically directing the motion and clustering of T-cell receptors (TCRs) on the cell surface. When TCR movement is impeded by externally applied physical barriers, the actin network exhibits transient enrichment near the trapped receptors. The coordinated nature of the actin density fluctuations suggests that they are composed of filamentous actin, but it has not been possible to eliminate de novo polymerization at TCR-associated actin polymerizing factors as an alternative cause. Here, we use a dual-probe cytoskeleton labeling strategy to distinguish between stable and polymerizing pools of actin. Our results suggest that TCR-associated actin consists of a relatively high proportion of the stable cytoskeletal fraction and extends away from the cell membrane into the cell. This implies that actin enrichment at mechanically trapped TCRs results from three-dimensional bunching of the existing filamentous actin network.The T-cell actin cytoskeleton is critical for proper antigen recognition by the mammalian adaptive immune system. During T-cell receptor (TCR) triggering by antigen peptides presented on major histocompatibility proteins (pMHCs) on the surfaces of antigen-presenting cells (APCs), the T-cell actin cytoskeleton adopts a pattern of centrosymmetric retrograde flow (1–3). This simultaneously promotes further TCR triggering (4) and rearranges various T-cell membrane proteins and their APC counterparts into an organized cell-cell interface termed the immunological synapse (IS) (5–7). During this process, TCRs form microclusters that move to the center of the IS in an actin-dependent manner (8,9). When engineered physical barriers interrupt the centripetal motion of TCR clusters, actin flow slows near the pinned microclusters, and the cytoskeletal network transiently accumulates and dissipates at the sites (10,11). The amplitude and duration of the induced cytoskeletal fluctuations are much greater than would be expected for a random distribution of independent objects, indicating that the actin in the local environment is coordinated. Whether this coordination arises from a rearrangement in the existing F-actin network or represents de novo polymerization of the cytoskeleton, as predicted by the association of TCRs with actin polymerizing factors (12), remains unclear. Here, we use a dual-probe cytoskeleton labeling approach that has previously been applied to distinguish between stable and dynamic populations of actin by exploiting the different relative affinities of monomeric actin and actin-binding proteins toward each population (13). This strategy reveals that TCR-associated actin is composed primarily of the stable cytoskeletal fraction and that local enrichment results from three-dimensional bunching of the existing filamentous actin network.Primary T cells from mice transgenic for the AND TCR were triggered using synthetic APCs consisting of supported lipid bilayers functionalized with pMHC and the integrin ligand intercellular adhesion molecule 1. Nanopatterned metal grids on the bilayer substrate acted as diffusion barriers that prevented lateral transport of TCR-pMHC complexes (14,15). Transient enrichment of actin at TCR clusters trapped at these barriers was visualized using fluorescent fusions of actin itself (mKate2-β-actin) and the F-actin binding domain of utrophin (EGFP-UtrCH). Such a dual-probe strategy theoretically allows for discrimination between different pools of actin: dynamic populations characterized by high polymerization and/or short filament fragments tend to be relatively better labeled by direct actin fusions whereas stable populations composed of longer filaments can support higher labeling by fluorescent fusions of F-actin binding proteins. This visualization method has been validated in Xenopus oocytes, where it distinguishes actin populations during wound healing (13). It has not been explicitly applied to T cells; however, simultaneous labeling of the Jurkat cell cytoskeleton using EGFP-actin and Alexa 568-phalloidin reveals distinct populations of actin consistent with the results expected from Xenopus (13,16).Our results show that the T-cell periphery is relatively enriched in mKate2-β-actin (Fig. 1
C, box 1), while EGFP-UtrCH dominates toward the center of the IS (Fig. 1
C, box 2). We infer from this probe distribution that the cytoskeleton at the T-cell periphery is composed of short fragments and is a site of active polymerization, whereas at the center of the IS, actin filaments are longer and predominantly stable. This is consistent with previous models of the T-cell actin network (3,16). An effective way to highlight each of these cytoskeletal regions is to consider the relative ratios of the two probes at each location. In this case, a high UtrCH/actin ratio corresponds to stable actin, and a high actin/UtrCH ratio corresponds to dynamic actin (Fig. 1
D). When T cells are treated with cytochalasin D, an inhibitor of actin polymerization, the overall UtrCH/actin ratio of the cell decreases as would be expected from a general decrease in polymerized actin (see Movie S7 and Movie S8 in the Supporting Material). However, it should be noted that photobleaching can also shift the UtrCH/actin ratio over time. We limit quantitative analysis of the ratio to its spatial gradients at a single time point, but such analysis is possible in systems that permit rigorous calibration for probe expression and photobleaching.Open in a separate windowFigure 1Ratiometric imaging of the cytoskeleton in live T cells distinguishes between dynamic and stable actin populations. (A) mKate2-β-actin, (B) EGFP-UtrCH, and (C) merged images of a triggered T cell show different actin pools. The cutouts in panel C correspond to (1) a region high in dynamic actin featuring short, polymerizing filaments and/or actin monomers and (2) a region with a stable actin population featuring longer filaments to which UtrCH can bind. (D) The UtrCH/actin ratio image highlights pools of relatively high UtrCH (red) or actin (blue). (Scale bars: 5 μm.)Actin enrichment at trapped TCR clusters incorporates both mKate2-β-actin (Fig. 2, A and C) and EGFP-UtrCH (Fig. 2, B and C). The relative UtrCH/actin ratio at these sites (Fig. 2
D, box 2) is quite high relative to nearby background areas (Fig. 2
D, box 1), indicating that the actin is derived primarily from the stable actin population.Open in a separate windowFigure 2Receptor-induced cytoskeletal enrichment at sites of pinned TCRs corresponds to a primarily stable actin fraction. (A) mKate2-β-actin, (B) EGFP-UtrCH, and (C) merged images of a triggered T cell interacting with a nanopatterned supported lipid bilayer show actin enrichment corresponding to putative sites of pinned TCRs. (D) The UtrCH/actin ratio is high at sites displaying actin enrichment, indicating a primarily stable actin fraction in (1) these regions compared to (2) nearby background areas. (Scale bars: 5 μm.)The three-dimensional distribution of TCR-associated actin was analyzed in dual-labeled live T cells using a spinning disk confocal microscope. The recordings show actin extending away from the cell membrane in the vicinity of trapped TCRs, while the rest of the actin cytoskeleton remains relatively flat (Fig. 3 and see Fig. S1 in the Supporting Material). These protrusions of actin away from the membrane surface are predominantly composed of stable, filamentous actin, as indicated by their relatively high UtrCH/actin ratio (Fig. 3
B).Open in a separate windowFigure 3Three-dimensional ratiometric imaging shows that actin enrichment extends away from the cell membrane. Single planes from (A) merged mKate2-β-actin and EGFP-UtrCH and (B) UtrCH/actin ratio three-dimensional stacks show actin enrichment at the cell membrane. Cutouts represent Z projections passing through sites of (1) enrichment and (2) nearby background regions. The color distribution in panel B is analogous to that in Figs. 1D and and22D, and is omitted for clarity. (Scale bar: 5 μm in the x axis only. Scale box: 1 μm.)Our interpretation of these results is that the filamentous actin network is relatively dense at sites of pinned TCRs. This is the simplest explanation out of several possibilities, one of which is formin-mediated mKate2-β-actin-deficient actin nucleation (17). Filament bunching at pinned TCRs can arise from consistent biophysical properties without assuming heterogeneity between the biochemistry of these receptors and other actin-associated proteins such as those at the cell edge, where locally high probe ratios are absent.Although TCRs are intentionally trapped as part of this experimental strategy, it is likely APCs can naturally impede TCR ligand mobilities under certain circumstances, and this has been shown to impact T-cell signaling (18,19). Actin architecture near cell surface proteins has been extensively studied in focal adhesions of fibroblasts (20), but the lack of stress fibers in T cells makes it unlikely that the two structures are similar. Thus, receptor-induced cytoskeletal enrichment at TCR clusters adds to the catalog of actin behaviors in situ, which is conveniently probed by techniques such as ratiometric dual-probe imaging in live cells. These techniques can be coupled to various spatial analysis algorithms to further extend their utility. 相似文献
8.
Fabian G?ttfert Christian?A. Wurm Veronika Mueller Sebastian Berning Volker?C. Cordes Alf Honigmann Stefan?W. Hell 《Biophysical journal》2013,105(1):L01-L03
We report on a fiber laser-based stimulated emission-depletion microscope providing down to ∼20 nm resolution in raw data images as well as 15–19 nm diameter probing areas in fluorescence correlation spectroscopy. Stimulated emission depletion pulses of nanosecond duration and 775 nm wavelength are used to silence two fluorophores simultaneously, ensuring offset-free colocalization analysis. The versatility of this superresolution method is exemplified by revealing the octameric arrangement of Xenopus nuclear pore complexes and by quantifying the diffusion of labeled lipid molecules in artificial and living cell membranes.Since its first demonstration in (live) cell imaging (1), stimulated emission depletion (STED) fluorescence microscopy has been realized in many variants. Particularly, the key phenomenon employed in this method, namely switching fluorophores transiently off by stimulated emission, has been accomplished with laser pulses varying from picoseconds to nanoseconds in duration, and from kHz to MHz in repetition rate. Because continuous-wave beams are suitable as well (2), STED microscopy has been implemented with rather different laser systems, ranging from model-locked femtosecond to continuous-wave laser diodes (3,4). Although it underscores the versatility of STED to modulate the fluorescence capability of a fluorophore, this wide range of options may confuse adopters when balancing simplicity, applicability, and resolution gain. The situation is exacerbated when implementing pairs of excitation and STED beams for dual-color colocalization studies (5,6).Here we report on a simple arrangement providing dual-color STED nanoscopy (Fig. 1) and molecular diffusion quantification down to ∼20 nm in (living) cells. The presented dual-channel STED microscope utilizes a single fiber laser providing a 20-MHz train of 775 nm wavelength pulses of 1.2-ns duration. This compact laser source enables STED on fluorophores emitting in the orange to red range. Specifically, we applied this laser on the orange dyes Atto590 and Atto594 (excitation: 595 nm; detection: 620 ± 20 nm), and the red dyes KK114 and Abberior Star635P (excitation: 640 nm; detection: 670 ± 20 nm). Although the spectra of the dyes are partially overlapping, the individual color channels can be separated without data processing (see Fig. S1 and Fig. S2 in the Supporting Material). Both channels are recorded simultaneously within 50 ns, using temporally interleaved pulsed excitation in combination with time-gated detection (5,7,8).Open in a separate windowFigure 1Fluorescence nanoscopy of protein complexes with a compact near-infrared nanosecond-pulsed STED microscope. (A) STED reveals immunolabeled subunits in amphibian NPC; raw data smoothed with a Gaussian filter extending over 14 nm in FWHM. The diameter of the octameric gp210 ring is established as ∼160 nm. Scale bar, 500 nm. (B) Individual NPC image showing eight antibody-labeled gp210 homodimers as 20–40 nm sized units and a 80 nm-sized localization of the subunits in the central channel.Because in STED microscopy, the STED doughnuts firmly determine the position of the fluorescently active molecules, the use of a single doughnut for both fluorophores guarantees that the two color channels are almost perfectly coaligned. The use of the doughnut even counteracts misalignments of the confocal excitation and detection channels (Fig. 2, and see Fig. S3), making STED microscopy particularly powerful for colocalization measurements.Open in a separate windowFigure 2Determination of the colocalization accuracy. Xenopus A6 cells, labeled with an antiserum against multiple NUP subunits in the central NPC channel and two secondary antibodies decorated with the fluorophores Abberior STAR635P and Atto594 were imaged by STED microscopy. (A) Upon overlaying both channels, a high degree of colocalization is directly visible. Scale bar, 200 nm. (B) Quantification of the colocalization by cross correlation of much larger images (see Fig. S3). The correlation is maximal for zero displacement of the images, proving colocalization. (C) Confocal image of monocolored fluorescent beads taken with improperly coaligned excitation beams (left). Improper coalignment spoils the colocalization accuracy in confocal imaging; the two channels should be perfectly coaligned, but they show a false offset as indicated by the color difference. The offset is quantified by the cross correlation of the two channels (right). (D) The STED image of the same beads (left) not only shows 10-fold improved resolution over the confocal image in panel C, but also improved colocalization, again quantified by cross correlation (right). Thus, by predetermining the position of emission, the STED doughnut counteracts errors induced by imperfect coalignment of the two confocal color channels (for details, see Fig. S3). Scale bars = 100 nm.The cross section for stimulated emission is lower at 775 nm as compared to that found at somewhat shorter wavelengths (5), yet STED pulse energies of ∼7 nJ in the focus are sufficient to yield a resolution of ∼30 nm and ∼20 nm in the orange and red channels, respectively (see Fig. S4). In addition, due to the lower peak intensity, the 1.2 ns pulses are likely to induce less nonlinear absorption and hence less photostress as compared to their more commonly used <0.2 ns counterparts (8,9). On the other hand, the pulses are only 2–4 times shorter than the typical lifetime of the excited state, which lessens their STED efficiency. This slight reduction is neutralized here by detecting photons emitted ∼1 ns after excitation (5,7,8).The potential of this straightforward implementation of STED microscopy is evident when imaging immunolabeled nuclear pore complexes (NPCs) of cultured Xenopus cells. Contrary to the confocal recording, STED microscopy reveals subunits of this protein complex, specifically the typical eightfold symmetry of its peripheral transmembrane protein gp210, along with a set of proteins in the central pore channel (Fig. 1, and see Fig. S5 and Fig. S6). Unlike in stochastic superresolution imaging of gp210 (10), the color channels are inherently coaligned and simultaneously recorded simply by executing a single scan. Apart from a weak smoothing and background subtraction applied to enhance image contrast, the images are raw.Because fluorescence off-switching by STED is an instant process, STED microscopy can be employed to study fast spatial translocations, such as the diffusion of molecules on the nanoscale (3). To benchmark the performance of our setup, we analyzed the diffusion of a fluorescent glycerophospholipid analog (11) by fluorescence correlation spectroscopy (FCS) in membranes of living mammalian PtK2-cells (Fig. 3). STED allowed us to reduce the diameter of the probed area from the 250 nm-sized diffraction limit down to 19 nm (FWHM), representing σ = 8 nm in standard deviation of a Gaussian fit. The attained subdiffraction area is 2.5 times smaller as compared to what has been reported in living cells to date (4). In model membranes, the smallest diameter was 15 nm (σ = 6.4 nm).Open in a separate windowFigure 3Nanoscale molecular diffusion analyzed by STED FCS. (A) For moderate and larger STED beam power PSTED, the resolution scales inversely with its square-root, attaining 15 nm in FWHM of the distribution of fluorescence emission in space, describing the measurement area. Note the relatively small threshold power PS = 1.4 mW, which implies that a large resolution gain is already attained for PSTED < 100 mW. (Inset) The resolution was determined by measuring the transit time of a fluorescent phospholipid-analog (DSPE-PEG-KK114) in a lipid model membrane through the detection area by FCS. (B) In living mammalian Ptk2-cells, the transit time of the lipid analog scales linearly with the detection area, revealing a diffusion constant Dlat = 0.33 μm2/s, and showing that this lipid analog diffuses largely freely in the plasma membrane down to <20 nm scales.In both measurements, the molecular transit time depends linearly on the probed area, indicating that the labeled lipid molecules diffuse essentially freely down to spatial scales of 20 nm. Accordingly, the anomaly exponent α was close to 1 with values of α > 0.85, showing only minor deviations from free diffusion (see Fig. S7). Because the diameter is inversely proportional to the square-root of the STED beam power, the resolution can be adapted to a particular application need (Fig. 3, A and B).In summary, our arrangement provides up-to-date STED microscopy resolution in offset-free colocalization recordings. The ready-to-use near-infrared laser pulses keep undesired single and multiphoton absorption low and leave the visible spectrum amenable for further studies. 相似文献
9.
George?W. Ashdown Andrew Cope Paul?W. Wiseman Dylan?M. Owen 《Biophysical journal》2014,107(9):L21-L23
We combine total internal reflection fluorescence structured illumination microscopy with spatiotemporal image correlation spectroscopy to quantify the flow velocities and directionality of filamentous-actin at the T cell immunological synapse. These techniques demonstrate it is possible to image retrograde flow of filamentous-actin at superresolution and provide flow quantification in the form of velocity histograms and flow vector maps. The flow was found to be retrograde and radially directed throughout the periphery of T-cells during synapse formation.Many biological processes are now being visualized with the use of superresolution fluorescence microscopy techniques. However, localization-based techniques primarily rely on fixed or slow moving samples to permit the collection of structural information. The 10-fold gains in resolution afforded by these superresolution techniques are usually possible through sacrificing the factors that originally made microscopy such a powerful tool: the ability to image live cells. In the case of stimulated emission depletion imaging, the scanning approach associated with this technique may fail to detect faster molecular events when imaging whole cellular regions.Structured illumination microscopy (SIM) is an alternative to these methods (1). It increases the resolution of conventional fluorescence microscopy twofold; it has the advantage of using a wide-field system, providing fast acquisition speeds of whole cells with relatively low laser powers; and it is compatible with standard fluorophores. By using a physical grating to produce interference patterns from a laser, periodic illumination is created. This patterned illumination causes information from higher spatial frequencies to be downmodulated (i.e., shifted) into the optical transfer function (support region) of the lens, resulting in higher-resolution spatial information being captured than is ordinarily obtainable.To quantify the directional motion of intracellular molecules, spatiotemporal image correlation spectroscopy (STICS (2)) was applied. Using spatial image correlation in time, STICS measures the similarity of pixels with those surrounding in lagging frames via a correlation function. The correlation function provides information on both flow velocities and directionality, while discounting static structures through the immobile object filter, achieved by subtracting a moving average of pixel intensities.The formation of an immunological synapse between T cells and antigen-presenting cells is a process requiring many dynamic (3) and subdiffraction-limited clustering events (4–6) to take place. The polymerization of actin is important for the spreading of cells over their target antigen-presenting cells (7), as well as cell mobility and migration (8). Retrograde flow of densely meshed cortical actin is observed at the basal membrane of synapse-forming T cells, where it may have a role in the corralling and clustering of signaling molecules at the plasma membrane (9), as well as at the leading edge of migrating cells (10). Filamentous actin is an extremely dynamic (7), densely packed, and thin (7-nm) structure (11,12).Here, we perform STICS on SIM data acquired on a total internal reflection fluorescence (TIRF) microscope system, which generated an evanescent field of 75-nm depth for excitation. To our knowledge, this is the first demonstration of an image correlation approach to quantify molecular dynamics on subresolution length scales using wide-field microscopy. To demonstrate the technique, we analyze two-dimensional actin flows in CD4+ T cells during immunological synapse formation, performed after cross-linking of antigen T cell receptors on a coverslip coated with specific antibodies.Fig. 1
a shows a schematic of the TIRF SIM setup. Excitation light (488 nm) passes through a polarizing module and then a phase-grating block, producing diffracted beams. These are then passed through a diffraction filter module to isolate the −1 and +1 order laser beams. These first-order laser beams are angled through the objective to produce total internal reflection conditions at the glass-water interface. The two evanescent waves interfere at the sample, producing structured illumination. The setup then produces lateral and rotational shifts through three orientations, producing nine raw images containing higher spatial frequencies than can normally be acquired by an objective using standard light microscopy. Fig. 1
b demonstrates the increased resolution obtained from TIRF SIM. Shown are the collected Fourier frequencies compared to those of a conventional microscope (dotted red line). Resolution was also measured using sparse 100-nm diameter fluorescent beads. Fig. 1
c shows a magnified image of these beads from which a line profile was obtained (yellow arrow). The full width at half-maximum of this profile (Fig. 1
d) gives a lateral resolution for the system of 120 nm.Open in a separate windowFigure 1(a) Schematic of the TIRF SIM setup. (b) Demonstration of the doubling of spatial resolution of collected frequencies through a Fourier transform (superimposed red circle demonstrating regular spatial frequency limits). (c) SIM reconstructed image of 100-nm bead (scale bar 0.5 μm). (d) (Plotted line) Bead showing full width at half-maximum of 120 nm.We then applied STICS analysis to quantify actin flow in T cell synapses acquired using TIRF SIM (Fig. 2). Fig. 2
a shows a schematic of the STICS analysis. From the raw data, immobile objects are first filtered by subtracting a moving average of the pixel values. Vector maps were obtained from correlation analysis of the time-series as previously published in Hebert et al. (2) and Brown et al. (13). Fig. 2
b shows a reconstructed TIRF SIM image of a mature T cell immunological synapse, representative of a time-point derived from the time series acquired at 1.28 fps (see Movie S1 in the Supporting Material). From this reconstructed image, two representative regions have been selected. In these regions, pseudo-colored actin flow vectors are overlaid onto the fluorescence intensity image. These range in magnitude from 0.01 μm/min (blue) to 5.61 μm/min (red). It can be observed that all flow vectors are directed radially toward the synapse center. A histogram of this flow is shown in Fig. 2
c. The histogram shows a peak retrograde flow velocity of 1.91 ± 1.27 μm/min. These data are representative of n = 7 T-cell synapses imaged by TIRF SIM.Open in a separate windowFigure 2(a) STICS analysis, performed by isolating mobile from immobile structures through a moving average filter (i) and binning a subset of pixels into blocks of superpixels (ii); the STICS software correlates spatial fluorescence fluctuations through time (iii). The code then outputs vector maps showing directionality and flow velocities. (b) TIRF SIM image of actin flow in a T cell 5 min after contact with a stimulatory coverslip. (Zoomed regions) Retrograde actin flow at the synapse periphery. (c) Histograms showing flow speed statistics of vectors from T-cell synapses (n = 7). 相似文献
10.
Emmanuel Adu-Gyamfi Michelle?A. Digman Enrico Gratton Robert?V. Stahelin 《Biophysical journal》2012,103(9):L41-L43
The Ebola virus causes severe hemorrhagic fever and has a mortality rate that can be as high as 90%, yet no vaccines or approved therapeutics, to our knowledge, are available. To replicate and egress the infected host cell the Ebola virus uses VP40, its major matrix protein to assemble at the inner leaflet of the plasma membrane. The assembly and budding of VP40 from the plasma membrane of host cells seem still poorly understood. We investigated the assembly and egress of VP40 at the plasma membrane of human cells using single-particle tracking. Our results demonstrate that actin coordinates the movement and assembly of VP40, a critical step in viral egress. These findings underscore the ability of single-molecule techniques to investigate the interplay of VP40 and host proteins in viral replication.The actin cortex below the plasma membrane of mammalian cells is essential for maintenance of cell shape and for cell movement. This cortex has also been found to play an essential role in the replication process of a number of viruses including West Nile virus (1), respiratory syncytial virus (2), influenza (3), and vaccinia virus (4). Additionally, actin has been found to play a central role in the assembly and budding of HIV-1 (5) whereas Marburg virus has been shown to use actin-enriched filopodia to exit the host cell (6). Actin has also been found to be packaged into Ebola-virus-like particles (VLPs) (7). Ebola virus, which causes severe hemorrhagic fever, harbors a single-stranded negative-sense RNA genome encoding seven proteins. Of these seven proteins, VP40 is the most abundantly expressed and has been found to play a central role in the budding of the virus from the plasma membrane (8). Whereas actin has been found in Ebola VLPs (7), the role of actin in Ebola VP40 assembly is still seemingly unknown. Here, we have used Raster image correlation spectroscopy (RICS) (9) and three-dimensional single-particle tracking (see Fig. S1 in the Supporting Material) (10) to investigate the dynamics of Ebola VP40 and actin. We report that preassembled VLPs (pVLPs) of Ebola VP40 require actin for directed movement and assembly.Ebola VP40 has been demonstrated to colocalize with actin and actin is found in VP40 VLPs (7), suggesting an important role for actin in the replication cycle of the virus. To confirm the colocalization between VP40 and actin in HEK293 and CHO-K1 cells, we used confocal microscopy to examine the distribution of EGFP-VP40 and mCherry-actin. EGFP-VP40 and mCherry-actin displayed colocalization at the plasma membrane of HEK293 and CHO-K1 cells (see Fig. S2
A), which was markedly reduced in response to treatment with LAT-A (see Fig. S2
B and Fig. S3
A), an actin polymerization inhibitor. VP40 plasma membrane localization was not disrupted by LAT-A treatment (not unexpected, as VP40 is a lipid-binding protein (11) where high affinity for the PM drives its cellular localization (E. Adu-Gyamfi and R. V. Stahelin, unpublished)). To test whether this VP40-actin interaction is important to viral egress, we detected EGFP-VP40 with an anti-EGFP antibody used to measure VLPs formed from cells expressing EGFP-VP40. This was also performed to assess the effect of pharmacological treatment on EGFP-VP40-expressing cells with LAT-A or with the microtubule polymerization inhibitor nocodazole (see Fig. S3
B). LAT-A treatment led to a significant reduction in VLP formation whereas nocodazole did not display detectable effects.To test whether the VP40 and actin are engaged in synchronized movement, we performed time-lapse imaging in both the green and red channels. We observed that the pVLPs move with actin fibers extending from the plasma membrane (see Movie S1 in the Supporting Material). The movement was rapid, and caused smaller particles to merge into larger filamentous forms. To further demonstrate that the motion of actin and VP40 spatially overlapped, we used RICS to obtain correlation maps of EGFP-VP40 and mCherry-actin (Fig. 1). The spatial cross-correlation map indicated significant overlap of VP40 and actin movement (Fig. 2, A–C) at the plasma membrane (Fig. 1 and see Fig. S6), but not in the cytosol (see Fig. S5 and Fig. S7). In contrast, EGFP-VP40 and mCherry-α-tubulin (see Fig. S8, Fig. S9, and Fig. S10) displayed no significant spatial cross-correlation at the plasma membrane (Fig. S11) or other regions of the cell (see Fig. S12), supporting the VLP egress data where inhibition of microtubule polymerization did not influence viral egress.Open in a separate windowFigure 1EGFP-VP40 and mCherry-actin RICS analysis at the membrane. (A) HEK293 cells expressing EGFP-VP40 and mCherry-actin were imaged for 100 frames at 256 × 256 pixels. (White scale bar = 2 μm.) (B) Average intensity image of EGFP-VP40 across the 100 collected frames. (Pink box) Used to select a region of interest to yield the (C) average EGFP-VP40 intensity image. (D) Average intensity image of mCherry-actin taken for 100 frames at 256 × 256 pixels was used to select the same region of interest as in panel B (pink box) to yield the (E) average intensity image of the mCherry-actin signal in this region. (F) The two-dimensional spatial cross-correlation analysis of panels C and E demonstrates significant cross-correlation of VP40 and actin signals.Open in a separate windowFigure 2Three-dimensional RICS correlation maps of VP40 and actin cross-correlate at the plasma membrane. (A) EGFP-VP40 and (B) mCherry-actin (Fig. 1 and see Fig. S6 in the Supporting Material) RICS autocorrelation functions. (C) Appreciable cross-correlation is observed for EGFP-VP40 and mCherry-actin at the plasma membrane.To test whether the motion of the pVLPs is directed by actin, we applied the three-dimensional orbital tracking method first introduced by Levi et. al. (10). Tracking of isolated particles (Fig. 3
A) in five different cells allowed determination of the pVLPs trajectories (Fig. 3
D), which suggested that the VP40 particles undergo a directed motion. To verify this, we plotted the mean-square displacement (MSD) curves for the pVLPs (Fig. 3
C), which confirmed the trajectory was characteristic of directed motion. Analysis of the intensity profile of the dynamic VP40 particles suggested that the intensity of the particle changes with respect to time. Bleaching is expected if the molecule is exposed to the laser beam for an extended period of time; however, an increase in intensities was observed along the trajectory of the green channel due to addition of VP40 molecules. This suggests that the movement of the particles along actin fibers promote multimerization and maturation of the pVLPs. When actin polymerization was inhibited in four different cells with LAT-B, the rapid movement (see Fig. S13) and the directed trajectories of the pVLPs were lost (Fig. 3, E and F). This was reflected in a change from directed motion to movement indicative of random then constrained diffusion (Fig. 3, E and F).Open in a separate windowFigure 3Actin directs the movement of VP40 particles. HEK293 cells transfected with EGFP-VP40 were imaged with an electronic zoom of 2000 mV, corresponding to 72 nm/pixel in both X and Y. (A) An isolated and representative VP40 particle (highlighted by white box, inset) was tracked as described in the Supporting Material. (B) Intensity profile of the pVLP in A demonstrates increases in EGFP-VP40 intensity along the trajectory. (C) MSD of the pVLP, which follows a ballistic motion with a velocity of 0.067 ± 0.01 μm2 s−1. (D) The three-dimensional trajectory of the particle shown in panels A–C. (E) MSD curve of VP40 particles yields random then constrained diffusion after LAT-B treatment with a mean velocity of 0.017 ± 0.006 μm2 s−1 (p < 0.001). (F) Three-dimensional trajectory of the same particle shown in panel E displays a random then constrained diffusion.Taken together, our findings demonstrate that the movement of the pVLPs is driven by actin. Analysis of the pVLPs trajectories also suggests that the motion of pVLPs on actin enables further addition of VP40 molecules. These findings raise important questions regarding contemporary understanding of Ebola assembly and egress. VP40 lacks a consensus actin-binding motif, suggesting an adaptor protein such as an actin motor protein may function in this process. For instance, Myo10 has been found to be essential to Marburg virus release (6); however, Marburg VP40-Myo10 direct interactions were not observed, suggesting other cellular adaptor proteins may function in this process. Given the pathogenic nature of the Ebola virus and the necessity of VP40 to the assembly and egress of the virus (8), the VP40-actin coordination represents, to us, a novel target for therapeutic development. 相似文献
11.
Charlotte?K. Colenso Yang Cao Richard?B. Sessions Jules?C. Hancox Christopher?E. Dempsey 《Biophysical journal》2014,107(10):L25-L28
Relaxation of a hERG K+ channel model during molecular-dynamics simulation in a hydrated POPC bilayer was accompanied by transitions of an arginine gating charge across a charge transfer center in two voltage sensor domains. Inspection of the passage of arginine side chains across the charge transfer center suggests that the unique hydration properties of the arginine guanidine cation facilitates charge transfer during voltage sensor responses to changes in membrane potential, and underlies the preference of Arg over Lys as a mobile charge carrier in voltage-sensitive ion channels.The response of voltage-sensitive ion channels to changes in membrane potential is mediated by voltage sensor domains (VSD) containing a transmembrane helical segment (S4) with a repeating motif of positively charged and hydrophobic amino acids (Fig. 1) (1,2). Changes in membrane potential drive the S4 helix through the membrane plane with the charged side chains (largely arginine) on S4 swapping Glu/Asp carboxylate partners that lie on less mobile elements of the VSD (2). Movement of S4 is coupled to the ion-conducting pore to transmit changes in membrane potential to channel gating (3).Open in a separate windowFigure 1Structures of the VSD of membrane domains before MD in a POPC bilayer. The S2 (pink) and S4 (blue) helices of the VSD of the hERG model (A) and Kv1.2/2.1 chimera structure (B) are highlighted. (C) Sequence alignment of S2 and S4 among homologous voltage-sensitive K+ channels.The VSD charge-pairing motif of K+ and Na+ channels is best represented in VSD states at zero membrane potential (S4 helix up) for which crystal structures exist for Kv1.2 (4), Kv1.2/2.1 chimera (5), and Nav channels (6,7). In these states, positively charged residues on the intra- and extracellular sections of the S4 helix are separated by a hydrophobic charge-transfer center (CTC) (1) or plug (8) containing a highly conserved Phe residue (Fig. 1). This plug restricts water incursion across the VSD, focusing the electric field across a narrow region near the bilayer center. In voltage-driven transitions between S4 down- and up-states, positively charged S4 side chains move across the CTC.The ether-à-go-go (eag) and eag-related family of voltage-sensitive K+ channels likely share similar charge pairing interactions with VSDs in other channels (9,10). However, eag VSDs contain an extra negative charge on S2 (underlined in Fig. 1
C) so that in hERG, Asp residues (D460 and D466) lie approximately one helical turn above and below the conserved charge-transfer center Phe (F463) (Fig. 1). This eag-specific motif might be expected to facilitate transfer of Arg side chains through the CTC and to stabilize the voltage sensor (VS) in the up state. We recently described an open state (VS-up) hERG model built on the crystal structure template of the Kv1.2/2.1 chimera and molecular-dynamics (MD) simulation of this model in a hydrated POPC bilayer (11). We have inspected an extended version of this simulation and identified transitions of a gating charge into the CTC despite the absence of a membrane potential change. These transitions are absent in equivalent MD simulations of the chimera structure in a POPC bilayer.Fig. 1 shows a single VS from starting structures of the hERG model and the chimera structure in a hydrated POPC bilayer, after restrained MD to anneal the protein-lipid interface (see Methods in the Supporting Material). Because the hERG model is constructed on the chimera structure according to the alignment in Fig. 1 the pattern of pairing between S4 charges and acidic VS side chains is equivalent in the hERG model and chimera structure.The arrangement of charge-paired side chains remains constant during MD in all subunits of the chimera (e.g., Fig. 2
E and see Fig. S2 in the Supporting Material). However, in two subunits of the hERG model the R534 side chain moves toward the extracellular side of the bilayer, sliding into the CTC to form a charge interaction with the extra Asp residue (D460 in hERG) that lies just above F463 (Fig. 2, A–C). This transition is facilitated by changes in side-chain rotamers of R534 and F463 as the planar Arg guanidine group rotates past the F463 ring, and the availability of D460 as a counterion for the R534 guanidine (Fig. 2). Movement of an Arg guanidine past the Phe side chain of the CTC is similar to that described in steered MD of an isolated VS domain (12).Open in a separate windowFigure 2Movement of the R534 side chain across the CTC in chain a of the hERG model simulation (A). Similar transitions are observed in chains a and b (panels B and C), but not chains c (D) or d (not shown), where the R534 side chain remains close to D466. In all subunits of the Kv1.2/2.1 chimera simulation, charge pairing of the starting structure (Fig. 1B) was maintained throughout (e.g., panel E and see Fig. S2 in the Supporting Material). (Black and blue lines) Distances from the Arg CZ or Lys ε atom to the two O atoms, respectively, of Asp or Glu.Mason et al. (13) have shown, using neutron scattering, that the low charge density guanidine cation (Gdm+) corresponding to the Arg side chain is poorly hydrated above and below the molecular plane. This property may underlie the universal preference for Arg (over Lys) in voltage sensor charge transfer. Although the poorly-hydrated surfaces of Gdm+ interact favorably with nonpolar (especially planar) surfaces (14,15), Gdm+ retains in-plane hydrogen bonding (13). In the transition of R534 across the CTC, in-plane solvation of the guanidine side chain is provided initially by D466, D501, and water molecules below the CTC, and during and after the transition by D501 and D460 side chains and waters above the CTC (Fig. 3, A and B). Complete transfer of the R534 side chain across the CTC was not observed, but would be expected to involve movement of the guanidine group away from H-bonding distance with D501.Open in a separate windowFigure 3In-plane solvation of R534 guanidine in the charge transfer center during the hERG model MD (A). (Dotted lines) H-bond distances of <2.5 Å. The right-hand group consists of top-down (B) and end-on (C) views of the distribution of oxygen atoms around the side chain of hERG R534 at 20-ns intervals during MD (subunit a). (D) End-on view of equivalent atom distributions around the K302 side chain during the Kv1.2/2.1 chimera MD (subunit c). (Red spheres, water O; pink, Asp OD1 and OD2; purple:, Glu OE1 and OE2.)The atom distribution around the R534 side chain during MD (Fig. 3, B and C) conforms to the experimental Gdm+ hydration structure (13), with H-bonding to waters and side-chain Asp O atoms exclusively in the guanidine plane. The passage of Gdm+ through the CTC is facilitated by the hydrophobic nature of Gdm+ above and below the molecular plane (13), which allows interaction with the nonpolar groups (especially F463) in the CTC (Fig. 3
A and see Fig. S3). This contrasts with the solvation properties of the Lys amino group (e.g., K302 of the Kv1.2/2.1 chimera (Fig. 1), which has a spherical distribution of H-bonding and charge-neutralizing oxygen atoms (Fig. 3
D and see Fig. S4).To further test these interpretations, we ran additional MD simulations of the isolated hERG VS domain model and an R534K mutant in a hydrated POPC bilayer. Again, the R534 side chain entered the CTC in the wild-type model simulation whereas the K534 side chain did not (see Fig. S5). Inspection of the atom distributions in Fig. 3
D (and see Fig. S4) indicates that the pocket below the conserved Phe of the CTC is particularly favorable for a Lys side chain, with waters and acidic side chains that satisfy the spherical solvation requirements of the terminal amino group, and nonpolar side chains that interact with the aliphatic part of the side chain.The occurrence of transitions of the R534 side chain through the CTC in the hERG model, in the absence of a change in membrane potential, indicates a relaxation from a less-stable starting structure. However, the path of the R534 side chain provides useful molecular-level insight into the nature of charge transfer in voltage sensors. How do these observations accord with broader evidence of charge transfer in voltage-sensitive channels in general, and hERG in particular? Studies with fluorinated analogs of aromatic side chains equivalent to F463 of hERG or F233 of the chimera indicate the absence of a significant role for cation-π interactions involving the CTC aromatic group in K+ and Nav channels, although a planar side chain is preferred in some cases (1,16). In hERG, F463 can be replaced by M, L, or V with small effects on channel gating (17), indicating that the hERG CTC requires only a bulky nonpolar side chain to seal the hydrophobic center of the VS and allow passage of the Arg side chain through the CTC. Both absence of requirement for cation-π interactions, and accommodation of nonplanar hydrophobic side chains in a functional hERG CTC, are broadly consistent with the interpretation that it is the poorly-hydrated nature of the Arg guanidine group above and below the molecular plane (together with its tenacious proton affinity (18)) that governs its role in carrying gating charge in voltage sensors.While the simulations suggest that R534 may interact with D460 in the open channel state, the possibility that the extra carboxylate side chain above the CTC might facilitate gating charge transfer is seemingly inconsistent with the slow activation of hERG, although hERG D460C does activate even more slowly than the WT channel (9). However, S4 movement in hERG occurs in advance of channel opening (19), and slow gating is partly mediated by interactions involving hERG cytoplasmic domains (20); thus, slow S4 movement may not be an inherent property of the hERG voltage sensor. Recent studies show that when hERG gating is studied at very low [Ca2+] (50 μM) and low [H+] (pH 8.0), the channel is strongly sensitized in the direction of the open state; this effect is reduced in hERG D460C (and hERG D509C) (10). These observations support a role for the extra hERG Asp residues in binding Ca2+ (and H+) (10), allowing the channel to be allosterically responsive to changes in pH and [Ca2+]. A true comparison of a hERG model with experimental channel gating might involve studies on a channel lacking cytoplasmic domains that modulate gating, and using conditions (high pH and low [Ca2+]) that leave the eag-specific Asp residues unoccupied. This could reveal the inherent current-voltage relationships and kinetics of the hERG voltage sensor. 相似文献
12.
Jeremy?S. Treger Michael?F. Priest Raymond Iezzi Francisco Bezanilla 《Biophysical journal》2014,107(6):L09-L12
Clinical methods used to assess the electrical activity of excitable cells are often limited by their poor spatial resolution or their invasiveness. One promising solution to this problem is to optically measure membrane potential using a voltage-sensitive dye, but thus far, none of these dyes have been available for human use. Here we report that indocyanine green (ICG), an infrared fluorescent dye with FDA approval as an intravenously administered contrast agent, is voltage-sensitive. The fluorescence of ICG can follow action potentials in artificial neurons and cultured rat neurons and cardiomyocytes. ICG also visualized electrical activity induced in living explants of rat brain. In humans, ICG labels excitable cells and is routinely visualized transdermally with high spatial resolution. As an infrared voltage-sensitive dye with a low toxicity profile that can be readily imaged in deep tissues, ICG may have significant utility for clinical and basic research applications previously intractable for potentiometric dyes.Voltage-sensitive dyes provide a way to observe cellular electrical activity without the physical limitations imposed by electrodes. Although these dyes can monitor membrane potential with a resolution of a few microns from large populations of cells (1), there are three obstacles that prevent the use of these dyes in many research settings, including clinical research:
- 1.Most voltage-sensitive dyes use visible wavelengths of light that prevent imaging of tissues beneath the skin.
- 2.Many of these dyes produce significant toxicity or off-target effects (2).
- 3.Before this report, to our knowledge, no voltage-sensitive dyes have ever been available for administration in humans, which has limited their value in biomedically focused research.
13.
Jing Xu Stephen?J. King Maryse Lapierre-Landry Brian Nemec 《Biophysical journal》2013,105(10):L23-L25
Although the disease-relevant microtubule-associated protein tau is known to severely inhibit kinesin-based transport in vitro, the potential mechanisms for reversing this detrimental effect to maintain healthy transport in cells remain unknown. Here we report the unambiguous upregulation of multiple-kinesin travel distance despite the presence of tau, via decreased single-kinesin velocity. Interestingly, the presence of tau also modestly reduced cargo velocity in multiple-kinesin transport, and our stochastic simulations indicate that the tau-mediated reduction in single-kinesin travel underlies this observation. Taken together, our observations highlight a nontrivial interplay between velocity and travel distance for kinesin transport, and suggest that single-kinesin velocity is a promising experimental handle for tuning the effect of tau on multiple-kinesin travel distance.Conventional kinesin is a major microtubule-based molecular motor that enables long-range transport in living cells. Although traditionally investigated in the context of single-motor experiments, two or more kinesin motors are often linked together to transport the same cargo in vivo (1–4). Understanding the control and regulation of the group function of multiple kinesins has important implications for reversing failure modes of transport in a variety of human diseases, particularly neurodegenerative diseases. Tau is a disease-relevant protein enriched in neurons (5,6). The decoration of microtubules with tau is known to strongly inhibit kinesin transport in vitro (7–9), but how kinesin-based transport is maintained in the presence of high levels of tau, particularly in healthy neurons, remains an important open question. To date, no mechanism has been directly demonstrated to reverse the inhibitory effect of tau on kinesin-based transport. Here we present a simple in vitro study that demonstrates the significant upregulation of multiple-kinesin travel distance with decreasing ATP concentration, despite the presence of tau.This investigation was motivated by our recent finding that single-kinesin velocity is a key controller for multiple-kinesin travel distance along bare microtubules (10). The active stepping of each kinesin motor is stimulated by ATP (11), and each kinesin motor remains strongly bound to the microtubule between successive steps (10,11). As demonstrated for bare microtubules (10), with decreasing ATP concentrations, each microtubule-bound kinesin experiences a decreased stepping rate per unit time and spends an increased fraction of time in the strongly bound state; additional unbound kinesins on the same cargo have more time to bind to the microtubule before cargo travel terminates. Thus, reductions in single-kinesin velocity increase the probability that at least one kinesin motor will remain bound to the microtubule per unit time, thereby increasing the travel distance of each cargo (10). Because this effect only pertains to the stepping rate of each individual kinesin and does not address the potential presence of roadblocks such as tau on the microtubules, we hypothesized in this study that single-kinesin velocity may be exploited to relieve the impact of tau on multiple-kinesin travel distance.We focused our in vitro investigation on human tau 23 (htau23, or 3RS tau), an isoform of tau that exhibits the strongest inhibitory effect on kinesin-based transport (7–9). Importantly, htau23 does not alter the stepping rate of individual kinesins (7,9), supporting our hypothesis and enabling us to decouple single-kinesin velocity from the potential effects of tau. We carried out multiple-kinesin motility experiments using polystyrene beads as in vitro cargos (8,10), ATP concentration as an in vitro handle to controllably tune single-kinesin velocity (10,11), and three input kinesin concentrations to test the generality of potential findings for multiple-kinesin transport. Combined with previous two-kinesin studies (10,12), our measurements of travel distance (Fig. 1
A) indicate that the lowest kinesin concentration employed (0.8 nM) corresponds to an average of ∼2–3 kinesins per cargo. Note that in the absence of tau, the observed decrease in bead velocity at the higher kinesin concentrations (Fig. 1
A) is consistent with a recent in vitro finding (13). At 1 mM ATP, htau23 reduced kinesin-based travel distance by a factor of two or more (Fig. 1, A and B). This observation is in good agreement with previous reports (7,8).Open in a separate windowFigure 1Distributions of multiple-kinesin travel distances measured at three experimental conditions, to verify the effect of tau (A and B) and to investigate the impact of single-kinesin velocity on the tau effect (B and C). Shaded bars at 8.7 μm indicate counts of travel exceeding the field of view. The mean travel distance (d; ± standard error of mean, SEM), sample size (n), and corresponding mean velocity (v; ± SEM) are indicated. MT denotes microtubule. Mean travel distance increased substantially at 20 μM ATP (C), despite the presence of htau23. This effect persisted across all three kinesin concentrations tested (left to right).Consistent with our hypothesis, reducing the available ATP concentration to 20 μM increased the multiple-kinesin travel distance by >1.4-fold for all three input kinesin concentrations (Fig. 1, B and C), despite the presence of htau23. The corresponding reduction in single-kinesin velocity with decreasing ATP concentration (10,11) is reflected in the ∼3.4-fold reduction in the measured bead velocities (Fig. 1, B and C). Therefore, the strong negative relationship between single-kinesin velocity and multiple-kinesin travel distance occurs not only for bare microtubules (10), but also for tau-decorated microtubules.What causes the observed increase in travel distance at the lower ATP concentration (Fig. 1, B and C)? In addition to the mechanism discussed above for the case of bare microtubules (10), an intriguing mechanism was suggested by recent studies of tau-microtubule interactions in which htau23 was observed to dynamically diffuse along microtubule lattices (14,15): reducing the stepping rate of a microtubule-bound kinesin may effectively increase the probability that a tau roadblock can diffuse away before the kinesin takes its next step.Perhaps surprisingly, although htau23 does not impact single-kinesin velocity (7,9), we observed a modest reduction in the average velocity of multiple-kinesin transport in experiments using tau-decorated microtubules (Fig. 1, A and B). This decreased velocity reflects a substantially larger variance in the instantaneous velocity for bead trajectories in the presence of htau23 (see Fig. S1 in the Supporting Material), as quantified by parsing each bead trajectory into a series of constant-velocity segments using a previously developed automatic software incorporating Bayesian statistics (16).To test the possibility that single-kinesin travel distance impacts multiple-kinesin velocity, we performed stochastic simulations (see the Supporting Material) that assumed N identical kinesin motors available for transport and included kinesin’s detachment kinetics (17). Previously, this model successfully captured multiple-dynein travel distances in vivo using single-dynein characteristics measured in vitro (18). In this study, we introduced one (and only one) free parameter to reflect the probability of each bound kinesin encountering tau at each step. When encountering tau, each kinesin has a 54% probability of detaching from the microtubule (interpolated from Fig. 2A of Dixit et al. (7)); the undetached kinesin is assumed to remain engaged in transport and completes its step along the microtubule despite the presence of tau.Remarkably, our simple simulation suggested that the tau-mediated reduction in single-kinesin travel is sufficient to reduce multiple-kinesin velocity (Fig. 2
A). The majority of the velocity decrease is predicted to occur at the transition from single-kinesin to two-kinesin transport (Fig. 2). Further decreases in cargo velocity with increasing motor number are predicted to be modest and largely independent of tau (Fig. 2
B). The results of our simulation remain qualitatively the same when evaluated at two bounds (40 and 65%) encompassing the interpolated 54% probability of kinesin detaching at tau (see Fig. S2).Open in a separate windowFigure 2Stochastic simulations predict a tau-dependent reduction in multiple-kinesin velocity, assuming that the only effect of tau protein is to prematurely detach kinesin from the microtubule (or, to reduce single-kinesin travel distance). (A) Average velocity of cargos carried by the indicated number of kinesins was evaluated at 1 mM ATP, and for four probabilities that a kinesin may encounter tau at each step. Mean velocity was evaluated using 600 simulated trajectories for all simulation conditions. Error bars indicate SEM. (B) Change in cargo velocity with each additional kinesin (ΔVel/kinesin) as a function of tau-encounter probability. These values were calculated from cargo velocities shown in panel A. Error bars indicate SEM.We note that our simple simulations do not consider the possibility that kinesin may pause in front of a tau roadblock, as previously reported in Dixit et al. (7). We omitted this consideration because the interaction strength between kinesin and the microtubule in such a paused state is unknown. In a multiple-motor geometry, could a paused kinesin be dragged along by the other motors bound to the same cargo? Could a tau roadblock be forcefully swept off the microtubule surface by the collective motion of the cargo-motor complex? Significant experimental innovations are necessary to specifically address these questions in future multiple-motor assays and to guide modeling efforts. Nonetheless, our simple simulation demonstrates that reducing single-kinesin travel distance is sufficient to decrease multiple-kinesin travel distance.Taken together, our observations highlight a nontrivial interplay between velocity and travel distance for kinesin-based transport in the presence of tau. We uncover a previously unexplored dual inhibition of tau on kinesin-transport: in addition to limiting cargo travel distance, the tau-mediated reduction in single-kinesin travel distance also leads to a modest reduction in multiple-kinesin velocity. We provide what we believe to be the first demonstration of the unambiguous upregulation of multiple-kinesin travel distance despite the presence of tau, via reducing single-kinesin velocity, suggesting a mechanism that could be harnessed for future therapeutic interventions in diseases that result from aberrant kinesin-based transport. 相似文献
14.
Andreas Anderluh Enrico Klotzsch Jonas Ries Alexander?W.A.F. Reismann Stefan Weber Martin F?lser Florian Koban Michael Freissmuth Harald?H. Sitte Gerhard?J. Schütz 《Biophysical journal》2014,106(9):L33-L35
Transmembrane proteins are synthesized and folded in the endoplasmic reticulum (ER), an interconnected network of flattened sacs or tubes. Up to now, this organelle has eluded a detailed analysis of the dynamics of its constituents, mainly due to the complex three-dimensional morphology within the cellular cytosol, which precluded high-resolution, single-molecule microscopy approaches. Recent evidences, however, pointed out that there are multiple interaction sites between ER and the plasma membrane, rendering total internal reflection microscopy of plasma membrane proximal ER regions feasible. Here we used single-molecule fluorescence microscopy to study the diffusion of the human serotonin transporter at the ER and the plasma membrane. We exploited the single-molecule trajectories to map out the structure of the ER close to the plasma membrane at subdiffractive resolution. Furthermore, our study provides a comparative picture of the diffusional behavior in both environments. Under unperturbed conditions, the majority of proteins showed similar mobility in the two compartments; at the ER, however, we found an additional 15% fraction of molecules moving with 25-fold faster mobility. Upon degradation of the actin skeleton, the diffusional behavior in the plasma membrane was strongly influenced, whereas it remained unchanged in the ER.Live-cell microscopy and three-dimensional electron tomography has boosted our understanding of endoplasmic reticulum (ER) dynamics and morphology. Proteins have been identified which regulate the formation of cisternae versus tubelike membranes, and the contacts between ER and the various cellular organelles have been studied in detail (1). Little information, however, is available when it comes to protein dynamics and organization within the ER membrane. Its complex three-dimensional topology hampers standard diffraction-limited fluorescence microscopy approaches: in fluorescence recovery after photobleaching, for example, the obtained diffusion coefficients can be several-folds off, if the ER morphology is not correctly taken into account (2). A method is therefore needed which allows for resolving molecular movements on length scales below the typical dimensions of the ER structures.In principle, single-molecule tracking would provide the required spatial resolution due to the high precision in localizing the moving point emitters: localization errors of <40 nm can be easily achieved (3). This technique has given rise to multiple studies, in which the paths of the diffusing objects were used to make conclusions on the properties of the environment; particularly, the plasma membrane has become a favorite target for such investigations, yielding precise determinations of the diffusion coefficients of a variety of membrane proteins or lipids (4).Here, we report what is, to our knowledge, the first application of single-molecule tracking for a comparative study of the diffusion dynamics of a membrane protein at the ER versus the plasma membrane. As the protein of interest, we chose the human serotonin transporter (SERT): it is a polytopic membrane protein containing 12 transmembrane domains, with both C- and N-termini residing in the cytoplasm. Stable SERT oligomers of various degrees were observed to coexist in the plasma membrane (5). Functionally, SERT (6) is a pivotal element in shaping serotonergic neurotransmission: SERT-mediated high-affinity uptake of released serotonin clears the synaptic cleft and supports refilling of vesicular stores (7). Wild-type SERT (SERT-wt) is efficiently targeted to the presynaptic plasma membrane, whereas the truncation of its C-terminus (SERT-ΔC30) retains the mutant protein in the ER (8). The N-terminal mGFP- and eYFP-fusion constructs of the two versions of SERT thus allowed us to specifically address SERT located at the ER (eYFP-SERT-ΔC30) or at the plasma membrane (mGFP-SERT-wt (7)).Our experiments were performed at 37°C on proteins heterologously expressed in CHO cells. Total internal reflection (TIR) illumination afforded a reduction in background fluorescence and allowed for selective imaging of single mGFP-SERT-wt molecules at the cells’ plasma membrane or single eYFP-SERT-ΔC30 molecules at plasma membrane-proximal ER (Fig. 1 and see the Supporting Material). TIR was particularly crucial for single-molecule imaging of the ER-retained mutant, where out-of-focus background would surpass the weak single-molecule signals in epi-illumination.Open in a separate windowFigure 1Schematics of the plasma membrane (PM) and a part of the ER containing mGFP-SERT-wt or the ER-retained eYFP-SERT-ΔC30 mutant, respectively. Both can be excited by total internal reflection fluorescence (TIRF) excitation. Experiments were carried out either on cells expressing mGFP-SERT-wt or eYFP-SERT-ΔC30.For both mutants, the majority of molecules were mobile: in fluorescence-recovery-after-photobleaching experiments we observed a mobile fraction of 82 ± 8% for mGFP-SERT-wt and 91 ± 4% for eYFP-SERT-ΔC30. For single-molecule tracking, the high surface density of signals was reduced by completely photobleaching a rectangular part of the cell in epi-illumination; after a brief recovery period, a few single-molecule signals had entered the bleached area and could be monitored and tracked at high signal/noise using TIR excitation. Samples were illuminated stroboscopically for till = 2 ms, and movies of 500 frames were recorded with a delay of tdel = 6 ms; the short delay times ensured that even rapidly diffusing molecules hardly reached the borders of the ER tubes between two consecutive frames. This illumination protocol was run for 20 times per cell, yielding ∼2500 trajectories per cell.The single-molecule localizations were first used to map those areas that are accessible to the diffusing proteins. eYFP-SERT-ΔC30 showed distinct hotspots, representing plasma membrane-proximal ER, excitable by the evanescent field (Fig. 2
A). These hotspots hardly moved within the timescale of an experiment (tens of minutes, see Fig. S1 in the Supporting Material); indeed, remarkable ER stability was previously observed using superresolution microscopy (9). In contrast, a rather homogeneous distribution was observed for mGFP-SERT-wt in the plasma membrane (Fig. 2
B).Open in a separate windowFigure 2Superresolution and tracking data at the ER and the plasma membrane. Superresolution images are shown for the ER-retained SERT mutant eYFP-SERT-ΔC30 (A) and for mGFP-SERT-wt in the plasma membrane (B). (C and D) Diffusion coefficients of eYFP-SERT-ΔC30 (C) and mGFP-SERT-wt (D) are shown as normalized histograms before (blue) and after (red) Cytochalasin D treatment. Data were fitted by Gaussian mobility distributions (see Table S1 in the Supporting Material for the fit results).Next, we compared the mobility of the observed proteins. Single-molecule localizations were linked to trajectories as described in Gao and Kilfoil (10), and the apparent diffusion coefficient, D, of each molecule was estimated from the first two points of the mean-square displacement membrane. The distribution of log10
D showed a pronounced single peak (Fig. 2
D). It could be well fitted by a linear combination of two Gaussian functions, with the major fraction (85%) characterized by Dwt = 0.30 μm2/s; a broad shoulder to the left indicates the presence of proteins that are immobilized during the observation period. In contrast, the mobility of the ER-retained mutant showed a substantially different distribution, containing two clearly visible peaks (Fig. 2 C). We fitted the data with a three-component Gaussian model: the main fraction (82%) behaved similar to SERT at the plasma membrane, with DΔC30 = 0.32 μm2/s. In addition, a large fraction (15%) with high mobility of DΔC30 = 7.8 μm2/s and a minor fraction (3%) with low mobility was observed. The proteins responded as expected to degradation of the actin membrane skeleton (red bars in Fig. 2, C and D): at the plasma membrane, the mobility of mGFP-SERT-wt increased 4.6-fold (mean values), whereas at the ER membrane there was only a minor change for eYFP-SERT-ΔC30 mobility (1.06-fold increase; note that the ER is not connected to actin filaments (11)).The observation of a high mobility subfraction at the ER membrane is surprising. In general, the presence of obstacles—irrespective of whether randomly distributed or clustered, mobile or immobile—reduces the diffusivity of mobile tracers in a membrane (12). It is generally assumed that the high protein density in cell membranes is responsible for the rather low fluidity when compared to synthetic membranes (compare, e.g., Saxton and Jacobson (13) with Weiss et al. (14)). Interestingly, the observed diffusion constant of 7.8 μm2/s is of similar order as the mobility determined for various proteins in synthetic lipid membranes (14). It is thus tempting to hypothesize the presence of extended protein-depleted regions of higher fluidity within the ER membrane; such membrane domains were indeed observed already at the plasma membrane (15). We were also concerned, however, that protein degradation fragments could have contributed to our data: the three-dimensional mobility of an 85-kDa protein is ∼10 μm2/s (16), similar to the high mobility diffusion constant of eYFP-SERT-ΔC30.We tested the two explanations by analyzing the spatial distribution of fast (DΔC30 > 1 μm2/s) versus slow trajectories (DΔC30 < 1 μm2/s) of eYFP-SERT-ΔC30 (Fig. 3). Both types of trajectories clustered in the same regions, and no segregation into ER subdomains was observable at the resolved length scales. This finding—on the one hand—disfavors freely diffusing protein fragments as the origin of the high mobility fraction. On the other hand, it calls for further experiments to identify the origin of the fast and the slow mobility subfraction. Interestingly, when analyzing all eYFP-SERT-ΔC30 trajectories we found that 80% of the molecules showed diffusion confined to domains of 230-nm radius (see Fig. S2). This size is clearly smaller than the lateral extensions of the visible ER regions observed in Fig. 3. The finding indicates domain formation at the ER membrane; domains are averaged out in Fig. 3 due to the long recording times. Note that free diffusion was observed for mGFP-SERT-wt at the plasma membrane (5).Open in a separate windowFigure 3Ripley’s K function analysis of the different mobility fractions in the ER. For the cell presented in Fig. 2, the first position of every slow (D < 1 μm2/s; red) and fast (D > 1 μm2/s; blue) trajectory was plotted in panel A. Contour lines indicate regions of ER attachment to the plasma membrane. In panel B, the point-correlation function L(r)−r is plotted for the slow (red) and fast (blue) fraction. Furthermore, the correlation between fast versus slow is plotted (green). All three curves show a peak at ∼450 nm, which agrees with the extensions of the ER attachment zones.In conclusion, we have shown that single-molecule tracking is feasible for constituents of the ER membrane. We found a surprising diffusion behavior of SERT resulting in the following:
- 1.A slow fraction showing mobility reminiscent of protein diffusion in the plasma membrane, likely reflecting SERT diffusing in protein-crowded regions of the ER membrane; and
- 2.A fast fraction showing 25-fold faster diffusion kinetics.
15.
16.
We use all-atom molecular dynamics simulations on a massive scale to compute the standard binding free energy of the 13-residue antimicrobial peptide indolicidin to a lipid bilayer. The analysis of statistical convergence reveals systematic sampling errors that correlate with reorganization of the bilayer on the microsecond timescale and persist throughout a total of 1.4 ms of sampling. Consistent with experimental observations, indolicidin induces membrane thinning, although the simulations significantly overestimate the lipophilicity of the peptide.Antimicrobial peptides are a component of the innate immune system of eukaryotes (1). As such, they must interact with pathogenic membranes, either during translocation or by disrupting their structural integrity (2). Here we examine the binding of the 13-residue cationic antimicrobial peptide indolicidin (3) (ILPWKWPWWPWRR-NH2) to a lipid membrane as a first step towards elucidating its mechanism of action.Molecular solutes interact with lipid membranes in many cellular processes (4). Computational approaches such as molecular dynamics simulations have been widely used to characterize these interactions (5). However, molecular dynamics simulations can require unfeasibly long times to reach equilibrium (6). Therefore, it is common to compute equilibrium properties of solute insertion into lipid bilayers using umbrella sampling (7) simulations in which the solute is restrained along the bilayer normal using harmonic restraining potentials, or umbrellas, centered at values distributed between bulk water and the bilayer center.It is often assumed that equilibrium properties rapidly attain convergence in umbrella sampling simulations; accordingly, convergence measures are rarely published (8). However, we have recently shown that umbrella sampling simulations require up to 100 ns per umbrella (3 μs in total) to eliminate systematic sampling errors in the standard free energy of binding, Δ, of an arginine side-chain analog from bulk water to a lipid bilayer (8). The fact that umbrella sampling has been used to investigate the bilayer insertion of substantially larger solutes (9) motivates a systematic evaluation of statistical sampling convergence of Δ for indolicidin in a lipid bilayer.To estimate Δ of indolicidin to a lipid bilayer, we conducted 60 sets of umbrella-sampling simulations while systematically varying the initial conformation. In each umbrella sampling simulation, each umbrella was simulated for 1.5 μs, yielding a total simulation time of 1.4 ms and 60 independent free energy or potential of mean force (PMF) profiles from bulk water to the center of a POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine) lipid bilayer.The PMF profiles indicate that indolicidin strongly binds to the bilayer, partitioning inside the lipid headgroups (Fig. 1, A and E). Importantly, the mean estimate of Δ decays exponentially with equilibration time teq, indicating that systematic sampling errors in individual simulations continued to decrease throughout the 1.5-μs interval as rare events led to more favorable states (Fig. 1
B). The low frequency of transitions to more favorable states exacerbates the requirement for massive sampling using multiple independent simulations.Open in a separate windowFigure 1PMF for indolicidin partitioning into a POPC bilayer. (A) Average PMF from 60 independent umbrella-sampling simulations based on 1 < t ≤ 1.5 μs/umbrella. (B) Average Δ from 50-ns time intervals per umbrella (teq < t ≤ teq+50 ns) as a function of equilibration time, teq. (Solid line) Single exponential fit to the mean over 0.5 < teq ≤ 1.5 μs. (C) Mean values of Δ from the 10-μs/umbrella simulations (crosses) together with the mean values of Δ (triangles) and exponential fit from panel B. PMF and Δ profiles obtained from each of the 60 independent simulations are shown in Fig. S1 in the Supporting Material. (D–F) Representative conformations after 1.5 μs of simulation at = (D) 3 nm, (E) 1.2 nm, and (F) 0.0 nm. To see this figure in color, go online.Computational limitations precluded extending all 60 sets of umbrella-sampling simulations to even longer times. Instead, we identified the two simulations at each umbrella that appeared to be most representative of equilibrium and extended each to 10 μs per umbrella (see Methods in the Supporting Material). The resulting estimates of Δ continued to decrease until teq = 4 μs (68 μs in total), after which they stabilized at the asymptotic limit of the exponential fit of the shorter simulation data (Δ = −26 ± 5 kcal/mol; Fig. 1
C).As indolicidin approaches the bilayer, it is drawn closer (Fig. 2
A) as salt bridges form between the peptide and the phospholipid headgroups (Fig. 2
B), inducing their protrusion (Figs. 1
D and and22
C). At large separation distances, this state is attained only when the peptide becomes highly extended (Fig. 2, D and E). As indolicidin is inserted more deeply, the surface of the lipid bilayer invaginates (Figs. 1
E and and22
C), maintaining peptide-lipid salt bridges (Fig. 2
B) and leading to the formation of a pore when the solute is near the bilayer center (Figs. 1
F and and22
C, and see Fig. S2, Fig. S3, Fig. S4, and Fig. S5 in the Supporting Material). These Boltzmann-weighted ensemble averages may not be mechanistically representative of nonequilibrium binding events (8,10).Open in a separate windowFigure 2Slow equilibration of bilayer and peptide. (A–D) Color quantifies conformational reorganization for teq < t ≤ teq + 100 ns as a function of teq and ||. (A) Deviation of insertion depth, z, from , Δz ≡ z − ; (B) number of peptide-lipid salt bridges, NSB; (C) volume change of the bilayer’s proximal leaflet in the radial vicinity of the solute, Vε; and (D) peptide end-to-end distance (EED). There is no sampling for t > 0.5 μs at || ≥ 4.5 nm. (E) Representative time-series of a trajectory at = 3.9 nm. (F) Representative conformation at 10 μs for || = 1.2 nm. To see this figure in color, go online.The reorganization of the peptide, the bilayer, and the ionic interactions between them became more pronounced with increasing simulation time at peptide insertion depths shallower than the global free energy minimum (|| >1.4 nm; Fig. 1
A and Fig. 2, B–D). These conformational transitions are likely the source of the systematic drift of Δ. Reorganization of the bilayer also controls the rate of equilibration during membrane insertion of an arginine side chain (8,9) and a cyclic arginine nonamer (11), suggesting that slow reorganization of lipids around cationic solutes presents a general impediment to simulation convergence.Consistent with the perturbation of membrane thickness observed by in situ atomic force microscopy (12), our results suggest that indolicidin insertion induces local thinning of the bilayer (Fig. 1, E and F, and Fig. 2, C and F). The different conformational ensembles sampled by the peptide in water and in the lipid bilayer (Fig. 2
D) are consistent with the observations that indolicidin is disordered in solution (13) and adopts stable conformations in the presence of detergent (14). Although the peptide’s conformation continued to change when it was deeply inserted (Fig. 2
D), the amount of water in the bilayer’s hydrophobic core converged relatively rapidly (see Fig. S2). Indolicidin can induce the formation of hydrated, porelike defects (see Fig. S2, Fig. S3, Fig. S4, and Fig. S5) but does not act as a chloride carrier (see Fig. S6 and Fig. S7). Future studies of the mechanism of indolicidin action will examine the effect of multiple peptide binding.The PMF profile presented in this Letter is strikingly different from that computed by Yeh et al. (15) using different force field parameters for indolicidin partitioning into a DMPC (1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine) bilayer, from which the binding free energy was estimated to be 0 kcal/mol (15). However, that study comprised only 25 ns per umbrella and likely suffers from systematic sampling errors induced by initial conditions (see Fig. S8).Our estimate of the binding affinity is much larger than the values obtained for indolicidin and large unilamellar POPC vesicles using isothermal titration calorimetry, −7.4 kcal/mol (16), and equilibrium dialysis, −8.8 kcal/mol (13). Such a discrepancy suggests that the relative accuracy of binding free energies for amino-acid side-chain analogs (8,9,17) does not necessarily extend to polypeptides. Although more work is needed to elucidate the source of this discrepancy, this study underlines the importance of attaining convergence before evaluating force-field accuracy.Importantly, this work also highlights the extensive sampling required to remove systematic errors induced by initial conditions in atomistic simulations of peptides in membranes. Slow equilibration of the system is due to rare transitions across hidden free energy barriers involving reorganization of the membrane. Two simple recommendations are 1), evaluating the time-dependence of ensemble averages, and 2), conducting multiple simulations with different initial conditions. We have recently shown that by using enhanced sampling techniques it is possible to identify the locations of hidden free energy barriers without a priori knowledge (9). Future research will examine strategies for speeding up the crossing of these barriers, such as optimized order parameters including bilayer reorganization and enhanced sampling techniques including a random walk along the order parameter (9). 相似文献
17.
Amy?P. Guilfoyle Chandrika?N. Deshpande Josep Font Sadurni Miriam-Rose Ash Samuel Tourle Gerhard Schenk Megan?J. Maher Mika Jormakka 《Biophysical journal》2014,107(12):L45-L48
The release of GDP from GTPases signals the initiation of a GTPase cycle, where the association of GTP triggers conformational changes promoting binding of downstream effector molecules. Studies have implicated the nucleotide-binding G5 loop to be involved in the GDP release mechanism. For example, biophysical studies on both the eukaryotic Gα proteins and the GTPase domain (NFeoB) of prokaryotic FeoB proteins have revealed conformational changes in the G5 loop that accompany nucleotide binding and release. However, it is unclear whether this conformational change in the G5 loop is a prerequisite for GDP release, or, alternatively, the movement is a consequence of release. To gain additional insight into the sequence of events leading to GDP release, we have created a chimeric protein comprised of Escherichia coli NFeoB and the G5 loop from the human Giα1 protein. The protein chimera retains GTPase activity at a similar level to wild-type NFeoB, and structural analyses of the nucleotide-free and GDP-bound proteins show that the G5 loop adopts conformations analogous to that of the human nucleotide-bound Giα1 protein in both states. Interestingly, isothermal titration calorimetry and stopped-flow kinetic analyses reveal uncoupled nucleotide affinity and release rates, supporting a model where G5 loop movement promotes nucleotide release.The hydrolysis of guanosine triphosphate (GTP) by GTPases, such as the oncoprotein p21 Ras and heterotrimeric Gα proteins, is a critical regulatory activity for cell growth and proliferation (1). Aberrant GTPases are consequently often implicated in tumorigenesis, developmental disorders, and metabolic diseases (2). Critical for the initiation of a GTPase cycle is the release of guanosine diphosphate (GDP), which allows GTP to bind and switch the protein from an inactive to an active conformation. The GTP is subsequently hydrolyzed to GDP and inorganic phosphate, returning the GTPase to an inactive conformation (3).Given that the release of GDP is the fundamental step in the initiation of a GTPase cycle, the detailed mechanism by which it is released has been under intense scrutiny. Studies using double electron-electron resonance, deuterium-exchange, Rosetta energy analysis, and electron paramagnetic resonance, have shown that the mechanism involves conformational changes in the nucleotide-coordinating G5 loop, one of five nucleotide recognition motifs (4, 5, 6, 7, 8, 9, 10, 11). Structural studies of eukaryotic Gα proteins and the intracellular TEES-type GTPase domain of the prokaryotic iron transporter FeoB (NFeoB) have also illustrated distinct conformations of the G5 loop, depending on the nucleotide-bound state (9, 12).Recently, we reported mutational studies of the G5 loop of Escherichia coli NFeoB, which illustrated a correlation between the sequence composition of the loop and the intrinsic GDP release rate (13). However, despite these observations, it is unclear whether the observed conformational changes in the G5 loop are a prerequisite for GDP release, or if the movement is a consequence of GDP release. To address this fundamental question, in this study we have used a combination of protein engineering and biophysical methods.Initially, to assess the relevance of conformational flexibility in the G5 loop, we aimed to create a protein chimera combining sequence and structural characteristics of both fast and slow GDP-releasing GTPases. We thus engineered a protein chimera using E. coli NFeoB as the scaffold (a protein with fast intrinsic GDP release) and substituted the G5 loop with that of a slow GDP-releasing protein (the human Giα1 protein; Gene ID 2770; Fig. 1
A (5)). GTP hydrolysis assays comparing wild-type (wt) NFeoB (wtNFeoB) and the protein chimera (ChiNFeoB) validated the integrity of the GTPase activities of both proteins (kcat = 0.46 and 0.36 min−1, respectively). To further assess the ChiNFeoB protein, we determined its crystal structure at 2.2 Å resolution (see Table S1 in the Supporting Material). The ChiNFeoB structure contains two molecules in the asymmetric unit, with molecule A bound to GDP. They are essentially identical to the nucleotide-bound wtNFeoB structure (root-mean-square deviation of 1.2 Å over 226 Cα atoms; Fig. 2).Open in a separate windowFigure 1Chimera model and structural comparison. (A) Illustration highlighting the chimera sequence change. (Orange) Sequence of the extended G5 loop from Giα1, which replaced the NFeoB sequence (gray). (B–F) Structural comparison of the G5 loop between (B) WT apo (PDB:3HYR) and nucleotide-bound (PDB:3HYT) NFeoB structures. (C) NFeoB nucleotide-bound and Giα1 (PDB:2ZJZ). (D) Nucleotide-bound NFeoB and chimera (Chi_GDP). (E) Nucleotide-bound chimera and Giα1. (F) Nucleotide-free (Chi_apo) and bound chimera protein. (G) Overview of the nucleotide binding site and structural overlay of chimera and Giα1 structures. To see this figure in color, go online.Open in a separate windowFigure 2Superimposition of nucleotide-bound NFeoB and chimera protein, with thermodynamic parameters. To see this figure in color, go online.However, the ChiNFeoB structure, when compared to the wtNFeoB structure, revealed an alteration in the conformation of the G5 loop, showing an extra turn on the N-terminal end of the α6 helix. This is structurally distinct from the wtFeoB protein, but with a conformation similar to that of the Giα1 protein (PDB:2ZJZ; Fig. 1, B–F). As in the crystal structures of wtNFeoB and Giα1, ChiNFeoB residues implicated in coordination of the nucleotide base maintain their positions in the G5 loop relative to GDP. In particular, residues Ala∗150 and Thr∗151 (NFeoB numbering, the asterisk indicates Giα1 chimera residue) are involved in electrostatic interactions with the nucleotide base moiety, analogous to the structures of both wtNFeoB and Giα1 (Fig. 1
G). Serendipitously, the second molecule in the asymmetric unit of ChiNFeoB (molecule B) was present in the nucleotide-free state. The two molecules (GDP-bound and nucleotide-free) are nearly identical (the superposition of molecules A and B yields a root-mean-square deviation of 0.36 Å over 229 Cα atoms), with the G5 loop adopting a nearly indistinguishable conformation compared to that of the GDP-bound molecule A (Fig. 1
F).Importantly, this conformation is independent of the crystallographic packing, inasmuch as the loop is not involved in any crystal contacts. In contrast, the structures of nucleotide-bound and nucleotide-free wtNFeoB illustrated a large conformational change in the G5 loop (Fig. 1
B). Hence, the substitution in the chimera extends the secondary structure of the α6 helix, and as hypothesized, the engineered ChiNFeoB protein has a G5 loop structure that is more conformationally stable than that of wtNFeoB.We subsequently measured the affinity of the ChiNFeoB protein for GDP using isothermal titration calorimetry (ITC). Nonlinear regression was used to attain the thermodynamic parameters (including the GDP binding affinity, Ka; the corresponding dissociation constant (Kd) was calculated from the equation Kd = 1/Ka). Interestingly, these measurements revealed the ChiNFeoB protein to have an almost 10-fold reduced affinity for GDP (82 vs. 9 μM measured for the WT protein; Fig. 2). In contrast, in a recent alanine scanning mutagenesis study of the G5 loop we observed a fivefold increase in affinity for GDP in a Ser150Ala mutant (2 μM) (14). This mutant protein has a coordination environment for the GDP base analogous to that of the ChiNFeoB protein (Fig. 1
A), indicating that it is not the presence of an alanine at position 150 that causes the reduced GDP affinity observed for the chimera protein. Instead, the analysis by ITC and comparison with previous mutagenesis studies indicates that the GDP binding site is less accessible in the ChiNFeoB protein, likely due to the introduction of conformational rigidity that accompanies the extension of secondary structural elements within the loop (Fig. 1
D).To further evaluate the functional characteristics of the chimera protein, we used stopped-flow fluorescence assays to determine the rate of nucleotide dissociation (koff) and association (kon) for the ChiNFeoB protein. The association rate for the GTP analog mant-GMPPNP was determined from the slope of a linear plot of protein concentration versus the observed association constant (kobs). The kon for the chimera was determined to be 3.20 μM−1 min−1 (Supporting Material), the dissociation rate (koff) of GDP for the chimera was determined to be 16.6 s−1 (vs. 144 s−1 for wtNFeoB; Designation mGMPPNP mGDP Protein kona (μM−1 min−1) koffb (min−1) Kdc (μM) kond (μM−1 min−1) koffe (s−1) NFeoB 8.1 ± 0.1 78.6 ± 1.6 9.7 15.9 144.7 ± 2.0 Chimera 3.2 ± 0.1 208.2 ± 1.3 65.1 0.2 16.61 ± 0.50