Tracking multiple conformations occurring on angstrom-and-millisecond scales in single amino-acid-transporter molecules
Abstract
Most membrane protein molecules undergo conformational changes as they transition from one functional state to another one. An understanding of the mechanism underlying these changes requires the ability to resolve individual conformational states, whose changes often occur on millisecond and angstrom scales. Tracking such changes and acquiring a sufficiently large amount of data remain challenging. Here, we use the amino-acid transporter AdiC as an example to demonstrate the application of a high-resolution fluorescence-polarization-microscopy method in tracking multistate conformational changes of a membrane protein. We have successfully resolved four conformations of AdiC by monitoring the emission-polarization changes of a fluorophore label and quantified their probabilities in the presence of a series of concentrations of its substrate arginine. The acquired data are sufficient for determining all equilibrium constants that fully establish the energetic relations among the four states. The KD values determined for arginine in four individual conformations are statistically comparable to the previously reported overall KD determined using isothermal titration calorimetry. This demonstrated strong resolving power of the present polarization-microscopy method will enable an acquisition of the quantitative information required for understanding the expected complex conformational mechanism underlying the transporter’s function, as well as those of other membrane proteins.
Editor's evaluation
This study uses a single-molecule polarization microscopy approach to identify the different conformation states that the arginine/agmatine antiporter AdiC transitions through during transport. Four states are identified and proposed to correspond to the key steps in the transport cycle, including inward-open, inward occluded, outward occluded and outward open, setting the stage for measurements of equilibrium constants and kinetics associated with transport. This is a cutting-edge and challenging approach that offers the potential for obtaining direct information of protein conformational equilibria that will be of interest to anyone studying membrane transport mechanisms.
https://doi.org/10.7554/eLife.82175.sa0Introduction
Biological membranes enclose individual cells and thereby separate them from their environments. Proteins embedded in the membrane play important roles such as the receptor-mediated signal transduction, and the ion-channel- or transporter-mediated movement of inorganic ions or organic molecules across the membrane. To accomplish these tasks, the protein molecules undergo many necessary conformational changes. An understanding of the mechanism underlying these changes requires the ability to resolve individual conformational states. However, three-dimensional (3D) protein-conformational changes are not only rapid but also occur usually on an angstrom scale. For these reasons, at the single-molecule level, the tasks to reliably resolve the multistate conformational changes of typical proteins, which occur in four dimensions (4D) on angstrom-and-millisecond scales, remain extremely challenging, not to mention the need for acquiring a sufficiently large amount of data for examining complex conformational mechanisms. For example, protein conformational changes have been probed using a fluorescence resonance energy transfer (FRET)-based method, which allows the deduction of distance between a fluorescence donor and an acceptor, generally on a nanometer scale (Stryer, 1978). This deduction requires knowing the actual fluorescence transfer efficiency, which is often difficult to assess because it depends on the local environment and, in most cases, on the relative orientation of the fluorophore pair. When a protein exists in multiple states, resolution and unambiguous identification of conformational states are even more challenging.
Importantly, while a protein molecule undergoes conformational changes, some secondary structures, e.g., α-helices, which are spatially constrained by other secondary structures, inevitably adopt unique spatial orientations in each conformational state. Thus, this feature offers the opportunity to track these conformational states, without the need of determining their detailed 3D features, by monitoring such an α-helix’s spatial orientation defined in terms of the inclination and rotation angles (θ and φ; Figure 1A) with a method of adequate resolution.
One effective way to track the orientation change of a protein is to monitor the emission polarization change of a bifunctional rhodamine attached to an α-helix (Figure 1B) using a polarization microscope (Sase et al., 1997; Warshaw et al., 1998; Ha et al., 1998; Adachi et al., 2000; Sosa et al., 2001; Forkey et al., 2003; Beausang et al., 2008; Rosenberg et al., 2005; Forkey et al., 2005; Fourkas, 2001; Ohmachi et al., 2012; Lippert et al., 2017; Lewis and Lu, 2019c). The polarization of individual emitted photons, unlike their travel direction, is not meaningfully affected by the diffraction caused by a so-called polarization-preserving objective. The documented resolutions of such polarization-based detection of rotation motion within a protein had been ≥25°, estimated on the basis of 2.5 times of the standard deviation (σ) of angle measurements (Rosenberg et al., 2005; Forkey et al., 2005; Ohmachi et al., 2012; Lippert et al., 2017). Thus, this technique has been limited to the investigation of proteins that undergo large-angle changes, e.g., the lever-arm of a myosin that rotates as much as 80°. Recently, our group assembled a polarization microscope with four polarized-emission-recording channels and tracked the orientation change of an isolated, soluble domain of the MthK K+ channel via a bifunctional fluorophore label attached to an α-helix within the protein (Lewis and Lu, 2019c; Lewis and Lu, 2019a; Lewis and Lu, 2019b). By finding optimal hardware, devising necessary numerical corrections for certain system parameters, and developing essential analyses, we have achieved an effective σ as low as 2°, translating to 5° resolution for detecting changes in both θ and φ. For reference, the estimated median radius of proteins is ~20 Å (Brocchieri and Karlin, 2005; Erickson, 2009), and a rotation of 5° or 10° of a site 20 Å away from the origin would lead to an 1.7 or 3.5 Å change in the chord distance. Thus, the capability to resolve this small angle change allows one to track protein-conformational changes that occur on an angstrom scale in proteins of typical sizes.
While we have succeeded in tracking the conformational changes in the artificially excised, soluble gating-domain protein of MthK, our main goal here is to further develop this method for tracking the conformational changes in a membrane protein. In principle, we can use a protein artificially engineered for this purpose, or a biological protein minimally or heavily modified to enable the required measurements. Simply put, what is required here is merely a membrane protein that adopts multiple conformational states, which can be reported by a fluorophore label, so that we can determine whether our method can (i) resolve its conformational changes by tracking orientation of the dipole of an attached fluorophore at the expected resolution set by the signal-to-noise ratio (SNR), (ii) collect a required large amount of data, (iii) analyze the data acquired from individual protein molecules that are attached to the support glass with some flexibility, and (iv) demonstrate how to extract all the equilibrium constants that fully determine the energetic relations among the resolved states, relations that constrain an equilibrium model.
We chose the bacterial transporter AdiC protein for its expected adaptation of multiple states that we can use as a convenient preparation to test the resolving capability of our method in a membrane protein (Gong et al., 2003; Iyer et al., 2003). An understanding of the detailed biological mechanism of AdiC is not the primary goal here.
For a general context, AdiC, a member of the amino-acid and polyamine organocation (APC) transporter superfamily (Jack et al., 2000; Casagrande et al., 2008; Bosshart and Fotiadis, 2019), is a critical component of a proton-extruding system in pathogenic enterobacteria, e.g., Escherichia, Salmonella, and Shigella, which helps the bacteria to survive the insult inflicted by a host’s highly acidic gastric defense barrier with a pH value of as low as 2 (Gong et al., 2003; Foster, 2004; Fang et al., 2007; Iyer et al., 2003; Krammer and Prévost, 2019). AdiC facilitates the movement of arginine (Arg+) into and agmatine (Agm2+) out of bacteria, along their gradients. Inside bacteria, Arg+ is rapidly decarboxylated to Agm2+ by the enzyme AdiA, consuming a proton (Gong et al., 2003; Foster, 2004; Fang et al., 2007; Iyer et al., 2003; Tsai and Miller, 2013). An exchange between extracellular Arg+ of a single positive charge and intracellular Agm2+ of two charges effectively extrudes H+. The direction of the net exchange of the two substrates is dictated by their natural energy gradient. However, as an intrinsic property, AdiC can facilitate the movement of a given substrate in either direction, where the transport of one type of substrate does not markedly depend on which type of substrate is on the opposite side.
Results
Sample preparations
To monitor the conformational changes of AdiC, we chose to attach a fluorophore label to the surface-exposing helix 6a in one of its two subunits (Figure 1; Gao et al., 2009; Fang et al., 2009; Gao et al., 2010; Kowalczyk et al., 2011; Ilgü et al., 2016, Ilgü et al., 2021). Hereafter, unless specified otherwise, AdiC simply refers to one of its two functionally independent subunits (Figure 1D and E). Helix 6A is spatially constrained by other secondary structures, and moves along with them, adopting differing spatial orientations between two different known structural states of AdiC (Figure 1E). A bifunctional rhodamine molecule was attached via two mutant cysteine residues, spaced seven residues apart, to helix 6a in the region extracellular to the substrate-binding site to avoid affecting the binding affinity (see ‘Discussion’). Such an attachment aligned the fluorophore dipole along the axis of the helix (Corrie et al., 1998; Figure 1B). Under the same labeling condition, there was little detectable fluorescent labeling in the absence of the mutant cysteine residues (Figure 1—figure supplement 1A and B). To minimize the background labeling, we removed two native cysteine residues. Removal of native cysteine residues in AdiC has been shown to have very limited impacts on its function (Tsai et al., 2012).
For microscopic examination, individual AdiC protein molecules were inserted into nanodiscs (Ritchie et al., 2009; Denisov et al., 2019). For attaching the protein molecules to streptavidin adhered to the polylysine-coated surface of a piece of coverslip glass, the protein was made to contain a biotin-moiety covalently linked to the N-terminus and the streptavidin-binding tags linked to the N- and C-termini in each of its two functionally independent and structurally symmetric subunits, totaling six sites available for the binding of streptavidin molecules (Figure 1C and Figure 1—figure supplement 1C; ‘Materials and methods’). Assessed with isothermal titration calorimetry (ITC), the protein resulting from the cDNA construct genetically engineered for the present purpose exhibited a KD of 104 μM for Arg+ (Figure 1—figure supplement 1D), which is within the previously reported range of 32–204 μM for AdiC (Fang et al., 2007; Casagrande et al., 2008; Gao et al., 2010; Tsai et al., 2012; Wang et al., 2014). To minimize the potential constraint of the tag-mediated attachment on the conformational movement of the protein, we included multiple flexible spacer sequences between (i) the Avi and Strep II tags, (ii) these tags and the N-terminus of AdiC, (iii) the thrombin cleavage and Strep II sequences, and (iv) these two sequences and the C-terminus (Figure 1—figure supplement 1C). Consequently, the inclination angle θ varied considerably among molecules. Thus, we needed to develop the analytical method discussed below to spatially align individual molecules.
The density of protein molecules on a cover slip was sufficiently low such that individual fluorescent particles could be readily resolved spatially on microscopic images. The probability of individual AdiC molecules being attached with a single fluorophore was optimized with an empirical protein-to-fluorophore ratio during the labeling procedure, and such molecules were identified during the offline analysis on the basis of a single-step bleaching of fluorescence.
Reconstitution of the protein into lipid-containing nanodiscs, tags, attachment, mutations, and the fluorophore label are all necessary components of the present method, potentially contributing to system errors. Currently, there are no other techniques that can be used to assess, under the same conditions, the ultimate impact of these factors on AdiC, which highlights the need to develop the present method. In terms of energetic impact, a practically relevant question would be whether the values estimated using the present method with all those potentially impactful factors are comparable to those values estimated for the wild-type protein by means of an already established method, such as ITC (see ‘Discussion’). As such, we would only know in the end whether the estimates by our method were valid, underscoring the importance of a judicious deliberation in designing the preparation and experiments on the basis of available information, which will help to maximize the chance to obtain a successful outcome.
Fluorescence intensity recordings
The fluorophores attached to individual AdiC protein molecules were excited by an evenescent field generated by a circularly polarized laser beam under a total internal reflection (TIR) condition. The light emitted from individual attached fluorophores was captured via the objective of a TIR fluorescence (TIRF) microscope. To assess the polarization of the captured fluorescence light, we first split it into two equal portions and then further split one portion into 0° and 90° polarized components ( and ) and the other portion into 45° and 135° components ( and ) (Figures 2 and 3; Lewis and Lu, 2019c). These polarized components of fluorescence were recorded with an electron-multiplying charge-coupled device (EMCCD) camera. Effectively, the intensity counts recorded in the four polarization channels encode the full 3D information regarding the orientation of the fluorophore dipole.
We integrated the individual intensity images (Figure 3A) and plotted the resulting values against time. As an example, a set of intensity traces for a molecule, each examined in the absence or the presence of 0.75 mM Arg+, is shown in Figure 3B and C. From these traces, we calculated the trace of total emitted intensity (Itot) using Equation 12 (equations with a number greater than 6 are given in ‘Materials and methods’; the bleaching step is shown in Figure 3—figure supplement 1).
From Equations 1 and 2 derived for ideal conditions, one can see that the orientation of the tracked fluorophore would be specifically reflected by the relative intensities, or underlying photon counts, of its four polarized components (Lewis and Lu, 2019c; for the solution of three channels, see Fourkas, 2001).
Given that either angle would be a function of an intensity ratio, its resolution should be primarily limited by the SNR of intensities. Thus, one can visually notice intramolecular motions that occur on the angstrom scale from the relative variations in the four intensities (Figure 3A, Figure 3—video 1).
Detection of intensity changes
As a protein molecule transitions from one conformation to another conformation, the orientation of the attached fluorophore changes with respect to the fixed polarization angles of the two polarized beam splitters. Consequently, the intensities of the four polarized fluorescence components undergo characteristic changes. For example, when the fluorophore’s mean dipole vector moved such that its φ angle increased from 0 to 90°, the intensity recorded in the 0° channel would decrease, whereas that in 90° would increase accordingly; the intensities of 45° and 135° channels would also characteristically change in directions opposite to each other. Compared to many other types of fluorescence-based methods, such as FRET, the expected concurrent, characteristic changes in all four channels here, which serve as a critical constraint in analysis, markedly increase the confidence of the detection of transitions, given that such types of changes are not expected for statistically random changes, namely, the so-called noise.
By examining the characteristic changes in the four fluorescence components with a so-called changepoint algorithm (Chen and Gupta, 2001), we detected the time point where a change in the fluorophore’s orientation occurred, which was brought about by the underlying protein conformational change. To do so, we assessed the change in the number of recorded photons () per unit time with 95% confidence by evaluating the likelihood of two alternative possibilities that a change did or did not occur within a given time interval. This assessment was based on concurrent changes in the recorded among all four polarized components, redundancy that markedly increased the confidence that identified transitions were genuine.
During individual consecutive 10 ms recording intervals, the mean for all four polarization components together was 92 (‘Materials and methods’), with an effective SNR of 7. As shown in Figure 3B and C, individual detected transitions are demarcated by the vertical lines in the black traces superimposed on the data traces. From the polarization properties of 92 photons on average, we could reliably detect the individual transitions in the fluorophore orientation among different conformational states at the intended time resolution. Thus, this method offers an exquisitely sensitive detection of changes in the fluorophore orientation.
Angle calculation and state identification
Resolution of individual conformational states in terms of θ and φ angles requires a higher SNR and thus a much greater number of photons than what is required for detecting the fluorophore’s orientation changes. One way to solve this problem would be to increase the number of emitted photons by raising the intensity of the excitation laser, but a strong excitation intensity would undesirably shorten the lifetime of the fluorophore. An alternative way to increase SNR is to use the total number of photons recorded from each of the four channels within the duration of an event when an examined protein molecule adopts a specific conformation, dubbed dwell time. Practically, given that angles were calculated from ratios of intensities, a calculation using the total or the average number of photons recorded during an event would yield the same result.
Our ability to identify individual state-transition points was a prerequisite for us to perform such a summation, or averaging, of the number of recorded photons for a given event because individual dwell times were demarcated by these points. As shown in Figure 3B and C, while the vertical lines in the black traces, superimposed on the observed intensities traces, indicate the individual orientation-transition points, the average intensities over individual dwell times are shown as the horizontal lines in the black traces.
From these black traces of , , , and , we calculated θL and φL traces that are color coded (Figure 3D and E); these two angles are defined in the standard laboratory (L) frame of reference for microscopy studies where the z-axis is defined as being parallel to the optical axis of the objective and the x-y plane parallel to the sample coverslip. All angle calculations were done using expanded versions of Equations 1 and 2 (Equations 11 and 13), which contain four necessary, predetermined system parameters. These previously established parameters numerically correct for incomplete photon collection (α) (Axelrod, 1979), depolarization caused by imperfect extinction ratios of the polarizers (f) (Lewis and Lu, 2019c) and by fast wobble motion of the fluorophore dipole (δ) (Forkey et al., 2005), and slightly different intensity-recording efficiency of the four polarization channels (g). In the angle calculation, we used the information from all photons (1740 photons on average) recorded from all four polarization channels during individual dwell times to determine the corresponding (mean) angles over these individual durations with adequately low σ. The relations between σ for the two angles and SNR among individual analyzed particles are shown in Figure 4, in which the highest value of <4° is translated to a resolution of better than 10°.
Subsequently, conformational state populations could be identified from both angles together to increase resolvability and confidence, without any preconceived kinetic model. For the ease in technique development, in a previous study of the isolated gating ring of the MthK channel, the orientation of individual protein molecules was strictly constrained such that their central axis was well aligned with the optical axis. Under such a stringent condition and for that specific protein, its conformations could be resolved from the differences in the θ angle alone. As such, the conformational states could be directly sorted in terms of θ and φ, which is mathematically a 2D operation. Because one angle was adequate there, the sorting effectively became a 1D operation. Here, we need to rely on the combined differences on a 3D sphere. The asymmetric characteristics of θ and φ, illustrated below, would introduce some biases in state identification if these two angles were used in a 2D sorting operation. The angles θ and φ are related to x,y,z coordinates in the corresponding Cartesian system as
and
where contains the information of all x, y, and z, whereas contains that of only x and y. Furthermore, the two angles are related to x and y differently. To ensure equal weighting of the contributions from all three dimensions in the sorting operation, we mapped the information of the two angles on a unit sphere defined in a Cartesian coordinate system and performed the sorting in terms of x, y, and z coordinates on the basis of the shortest distance principle.
This 3D-sorting task would be extremely challenging for us to do manually on subjective visual guidance, if possible. Thus, to effectively sort individual events into distinct populations, we developed a k-means-clustering algorithm to achieve the solution for the shortest-distance case, which was optimized with two coupled algorithms (simulated annealing plus Nelder–Mead downhill simplex) (Press et al., 2007). The program examined from the case of a one-state distribution to the case of k-state distributions, one at a time in a series of separate operations. When two consecutive events were determined to belong to the same state distribution, the transition between them identified by the changepoint algorithm would most likely be part of the expected false-positive outcomes. These events would then be merged to form a single event for the state distribution. This operation reduces the false-positive identification of transition points.
For each successful operation with a specifically set number of states, the resolution criterion was ensured, that is, peaks of all populations were separated by at least 2.5σ. By this common criterion, the highest number of resolvable states was 4 over the examined Arg+ concentration range (Figures 5 and 6); in these two figures, θ and φ are the angle coordinates of a local framework transformed from θL and φL of the standard laboratory framework, as described below. The conformational states are numerated according to the sequence in the solution to achieve the shortest pathway among the four states, denoted as C1, C2, C3, and C4 (‘Materials and methods’). To illustrate the sorting result, x, y, and z coordinates of individual events in the Cartesian system or their angle coordinates in a spherical system are presented in Figure 7, in which they are color-coded according to states. Thus, the minimal model must have four states for either the apo or the substrate-bound forms, totaling eight states.
-
Figure 5—source data 1
- https://cdn.elifesciences.org/articles/82175/elife-82175-fig5-data1-v1.xlsx
-
Figure 5—source data 2
- https://cdn.elifesciences.org/articles/82175/elife-82175-fig5-data2-v1.xlsx
-
Figure 5—source data 3
- https://cdn.elifesciences.org/articles/82175/elife-82175-fig5-data3-v1.xlsx
-
Figure 5—source data 4
- https://cdn.elifesciences.org/articles/82175/elife-82175-fig5-data4-v1.xlsx
Spatial alignment of individual molecules
For greater likelihood to accurately estimate angle values, we need to determine their mean values from individual events of numerous molecules that had already been analyzed individually for necessary precision of angle measurements and thus spatial resolution. This operation requires all molecules be aligned in the same orientation. Unfortunately, individual molecules on the coverslip, and consequently the tracked helix, were randomly oriented in the x-y plane relative to the x-axis, that is, their φL varying randomly among different molecules. This problem previously prevented us from building the distribution of φL for a given state among individual molecules (Lewis and Lu, 2019c). Furthermore, each dimer molecule of AdiC is anchored to the cover slip coated with streptavidin through two available biotin moieties and four streptavidin-binding tags, each of which was covalently linked to an N- or C-terminus of the polypeptides of two subunits. As mentioned above, these terminal binding regions were of some flexibility. Consequently, the twofold symmetry axis of the individual dimer molecules was not aligned with the optical (z) axis of the microscope framework. As such, the orientation of one molecule varied considerably from that of another (e.g., Figure 7A versus Figure 7B or Figure 7E versus Figure 7F). Resolving this issue would also pave the way for studying individual randomly oriented molecules in the future.
To align all the molecules during analysis, we mathematically rotate them from the laboratory frame of reference into a local coordinate system defined on the basis of the spatial features of the tracked helix in the protein (Figure 7—figure supplement 1A). This system is defined such that the tracked helix in C1 is always aligned with the local x-axis, that is, mean φ1=0°, and the helix in the plane defined by C1 and C4 is always in the local x-y plane, that is, mean θ1 or θ4=90°. These two features fully define the x,y,z-axes, and thus θ and φ in the local framework of AdiC. Evaluated in this local framework (Figures 3F, G, 7C, D, G and H), individual events of the same state observed with all molecules under a given condition could be used to build a single distribution (Figure 5). From the distribution for each state, we determined the mean θi and φi. To specify the six mean direct-angle Ωi,j values among the tracked dipole orientations in the four states, which were calculated using Equation 33, we used four unit-vectors to specify the orientations in the local framework (Figure 8A). The mean θi, φi, and Ωi,j values are all summarized in Figure 8B and C.
Following the rotation operation, we could also build a distribution of mean angle values of individual molecules to illustrate the molecule-to-molecule variation in θi and φi (Figure 9). If the contribution of each molecule were weighted by the length of the trace, the mean value of each distribution in Figure 9 should statistically be the same as the corresponding one in Figure 5, built with the values of individual events from all molecules. Although the ensemble mean angle values should be more accurate statistically, the apparent angle resolution would be poorer. Indeed, compared with σ for individual particles (Figure 4), the ensemble angle distributions built with data determined from molecule-by-molecule analyses have much larger σ (8–11°, Figures 5 and 8B), except for those of φ1 and θ1 or θ4 which are normalized such that their mean values equal 0° and 90°. Thus, the initial molecule-by-molecule data analysis, as shown in Figure 3, was essential for resolving individual conformational states in terms of θ and φ.
Determination of equilibrium constants among the Apo and Arg+-bound conformational states
As an essential evaluation of the basic energetic properties of the protein’s individual states, we examined individual AdiC molecules in the presence of a large series of Arg+ concentrations. Both sides of an AdiC molecule facing the same solution allowed us to determine the probabilities of individual conformational states under conditions where the system as a whole was in equilibrium, through which we could straightforwardly obtain equilibrium constants (Ki,j) among the states, including dissociation constants (KDi) to be compared in ‘Discussion’ with previously reported KD values.
Shown in Figure 5 are plots of the probability density distributions of θi and φi of each of the four states in a series of Arg+ concentrations, built with the data from numerous individually analyzed particles. These angle distributions contain the information regarding the probabilities of occupying each of the four states in the presence of the corresponding Arg+ concentrations (Figure 10A). All four states appeared in the absence and the presence of Arg+. With increasing concentration, the probability of C3 increased whereas that of C2 decreased. In contrast, the probability of C1 or C4 exhibited relatively small changes. Nonetheless, these observations indicate that all states bind Arg+ because any state that could not bind Arg+ would practically vanish in a saturating concentration of Arg+. The observation of all four conformations in the absence or the presence of Arg+ is consistent with the Monod–Wyman–Changeux model for ligand-dependent conformational mechanism of allosteric proteins in that all conformations of a protein occur in the absence of ligand, and the binding of a ligand does not create a new conformation but merely energetically stabilizes a spontaneous conformation that is ready to capture the ligand (Monod et al., 1965).
Here, a model with a minimum number of eight states is required to account for the observed conformational behavior of AdiC: a set of four without ligand bound (C1, C2, C3, and C4) and another set of four with the Arg+ ligand (L) bound (C1.L, C2.L, C3.L, and C4.L) (Figure 10B). The probability pi of occupying state i is expressed as
where C1 is used as a reference for the other states (Ci) to define equilibrium constants:
where K1,1 is defined to equal one. Each equilibrium constant reflects the free energy difference between a pair of states that may or may not be connected kinetically. In principle, the equilibrium constants for the apo or Arg+-bound states are constrained by the state probabilities under the condition of zero or saturating Arg+, and by the so-called midpoint positions of the curves (Figure 10A). In practice, we determined these parameters by fitting Equation 5 to the four plots in Figure 10A simultaneously as a global fit, all summarized in Table 1. Together, these constants would fully define the energetic relations among the eight states. If needed, the remaining six equilibrium constants (, , , , , and ) could be calculated from the ones given in Table 1. Note that even though the concentration-dependent plot of the probability of C1 or C4 is relatively flat, their values are determinable. This is because the model can be fully specified by a certain combination of merely seven equilibrium constants, e.g., three of the six values, three of the six values, and only one of the four values. As such, in the present so-called overdetermined case, the values of the relatively flat traces are fully constrained by those more curved traces in a global fit to all four plots simultaneously.
Discussion
In this study, on the basis of θ or φ alone, the smallest angle changes resolved, namely, those between C2 and C3, are ≥10° (Figure 8B), which requires a σ ≤ 4°. Thus, it is justified to use ultimately the intensity data with SNR of ≥ 5, which corresponds to σ for the θ and φ distributions of < 4°, translated to a minimum resolution of better than 10° (Figure 4). Note that our actual use of the combined information of θ and φ should lead to higher effective resolution.
The present examination of the properties of conformational changes was performed under equilibrium conditions, which allowed us to determine straightforwardly the KD values of four conformations from the dependence of their probabilities on the ligand concentration. The resulting values of KD1 through KD4 ranges from 25 to 69 μM (Table 1), statistically comparable with the previously reported overall KD determined by ITC: a relatively narrow range of 32–93 μM for one lab without lab-to-lab variations (Fang et al., 2007; Tsai et al., 2012), or a wider range of 32–204 μM when data from other labs are included (Casagrande et al., 2008; Gao et al., 2010; Wang et al., 2014). Thus, the KD values estimated by the present method are valid. Maneuvers such as introducing cysteine mutations and attaching the fluorophore to the AdiC molecule, which are part of what contributes to system errors of the present method, have no markedly consequential energetic impacts on its affinity for Arg+. Such a finding is not particularly surprising because the chosen labeling part is on the surface of the protein such that it is not at, but external to, the ligand-binding site. Thus, as an advantage of the present method, one can and should choose to attach the label to the surface of a moving part of the protein, a part that is not of functional activity.
Generally, a transporter is expected to adopt four main types of structure-function states in terms of the accessibility of its external and internal sides to substrates, dubbed the externally open (Eo), externally occluded (Ex), internally open (Io), and internally occluded (Ix) states (Post et al., 1972; Gao et al., 2010; Krammer and Prévost, 2019). By this definition, when a transport molecule adopts the Eo and Io states, it is accessible only to extracellular and intracellular ligands, respectively. The molecule in the Ex or Ix state is inaccessible to ligands from either side. Thus far, the crystal structures of AdiC in the Eo and Ex state have been solved (Gao et al., 2009; Fang et al., 2009; Gao et al., 2010; Kowalczyk et al., 2011; Ilgü et al., 2016, Ilgü et al., 2021). Additionally, the structures of BasC and ApcT (Shaffer et al., 2009; Errasti-Murugarren et al., 2019), which share the same fold with AdiC, were solved in Io and Ix states, respectively.
Intriguingly, as shown in the results, the present study of fluorescence polarization resolved four macroscopic conformational-state populations, which was performed without any preconceived number of states. The probabilities of C2 and C3 clearly vary with the concentration of Arg+, and these two states should thus be open states (Figure 10A). In contrast, the probabilities of C1 and C4 vary little, and thus are not open states directly accessible to ligands, which are consistent with occluded states in adopting the prevailing nomenclature. Under the present conditions, the comparable probabilities of C2 and C3 in a practically saturating concentration of Arg+ are consistent with the following experimental findings obtained from previous flux assays: both sides of AdiC have comparable Km, and the apparent maximal net flux rates of radioactive Arg+, which were separately measured for the two opposite directions, are also comparable (Tsai et al., 2012; confer Krammer et al., 2016).
In the previous study of the isolated gating ring of the MthK channel, we could easily relate the states identified in the polarization study to the structural states on the basis of the matching relations of θ angles alone, which is effectively an 1D operation. For ease of communication, the former types of state will be referred to as the conformational states. However, to relate two sets of structural and conformational states in a 3D operation would require a quantitative approach. Thus, as an additional exercise, we illustrate below such an approach that we have developed.
Thus far, only the structures of the Eo and Ex states of AdiC have been solved, so for this exercise we used the structures of the transporters BasC and ApcT in the Io and Ix states as proxies for those of AdiC because they all have the same structural fold (Shaffer et al., 2009; Gao et al., 2010; Errasti-Murugarren et al., 2019; Ilgü et al., 2021). This exercise can be readily repeated when the Io and Ix structures of AdiC become available. The helix 6 in the two AdiC structures and that of its counterparts in BasC and ApcT adopts different orientations. To illustrate these differences, the four structures are aligned spatially, and their relevant portions are shown in Figure 11A and B. The orientation of the helix 6 or its counterpart is represented by a vector that is color-coded for a specific structural state (open arrow heads, Figure 11C). There are 24 possible ways of relating the four structural states to the four conformational states. To determine the best matching relation, we fit the unit-vectors (open heads, Figure 11C) and consequently their θ and φ (open circles, Figure 11D), which specify the helix’s orientations in the four structural states, to those of the conformational states in the local framework (closed heads or circles), while the internal spatial relations among four states in either set of compared data were fixed. To our surprise, in the statistically best matching case (Figure 11C), the four mean orientations of the fluorophore attached to the helix 6 apparently well match the orientations of this helix in the structures of AdiC, BasC, and ApcT captured in four different states (Figure 8B and C).
For further illustration, we calculated the combined least-square () values between the two sets of open and closed arrowheads in Figure 11C for all 24 combinations, and plotted the values in Figure 11E. In the statistically best matching case based on the largest , C1 would correspond to Ex, C2 to Eo, C3 to Io, and C4 to Ix, in which the values θi and φi and Ω are comparable between the corresponding conformational and structural states (Figure 8B and C). In terms of open versus occluded states, the assignments of C2 and C3 as open states, and C1 and C4 as occluded states are consistent with those based on the Arg+ concentration dependence, discussed above. In terms of sidedness, judging from the structures of AdiC, that is, Eo and Ex, the spatial orientations of helix 6A in them (θ and φ) are compatible with those of C2 and C1.
In summary, using a state-of-the-art fluorescence polarization microscopy system assembled in house, we have, for the first time, successfully tracked conformational changes in a single integral membrane protein molecule among four states that occur on angstrom and millisecond scales. The resolving power of this technique allowed us to experimentally establish a fully determined quantitative model of eight states required to account for the energetics of the observed conformational changes of an AdiC transport molecule with or without a ligand bound. This capability of resolving and tracking the conformational states of individual molecules forms the foundation for performing the kinetic study to acquire the necessary dynamic information for understanding the transporting mechanism, and for creating an experiment-based 4D model to quantitatively account for the complex spatiotemporal behaviors of a transporter molecule, which can be applied to the investigation of other membrane proteins.
Materials and methods
Materials
Detergent n-dodecyl-β-D-maltopyranoside (DDM) was purchased from Anatrace, 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) from Avanti Polar Lipid Inc, bifunctional rhodamine bis-((N-iodoacetyl)-piperazinyl)-sulfonerhodamine from Invitrogen (B10621), strep-tactin resin from IBA, and cover slip (#1.5) and microscope slide glass from Fisher Scientific or VWR. Unless specified otherwise, all other reagents were purchased from Sigma, Thermo Fisher Scientific, or EMD Millipore.
Cloning, protein expression, and purification
Request a detailed protocolA double-stranded DNA fragment, synthesized by Integrated DNA Technologies (IDT), contains, from the N-terminus to the C-terminus, Avi tag for recognition by biotin ligase, a Strep II tag, a linker (GGGSGGGS), the gene of AdiC of Escherichia coli, a linker (GGGS), a thrombin protease recognition site, a C-termini Strep II tag, and a stop codon, which was cloned into the pET28b vector. The removal of two native cysteines (C238A and C281A) and introduction of the double G188C and S195C cysteine mutations in helix 6a for attaching bifunctional rhodamine were carried out using the QuickChange technique and verified by DNA sequencing.
The AdiC protein was expressed in E. coli BL21(DE3) cells transformed with AdiC-gene-containing plasmids. The transformed cells were grown in Luria Broth at 37°C to an A600 of ~1.0. Protein expression was induced with 0.5 mM iso-propyl β-D-thiogalactopyranoside (IPTG) at 22°C overnight. The cells were harvested and resuspended in a solution containing 100 mM NaCl, 50 mM tris-(hydroxymethyl)-aminomethane (Tris) titrated to pH 8.0, 1 mM phenylmethylsulfonyl fluoride (PMSF), 4 mM tris-(2-carboxyethyl)-phosphine (TCEP), 1 µg/ml leupeptin, and 1 µg/ml pepstatin A. To extract membrane proteins, 40 mM DDM was added to the cell suspension; the conical tube containing this suspension was placed on a rotating device in a 4°C cold room for 2–3 hr. The cell suspension was then sonicated, the cell lysate was centrifuged at 12,000 × g for 20 min, and the resulting supernatant was loaded onto a gravity flow column packed with strep-tactin super-flow resin (IBA). The column was washed with a wash buffer (WB) containing 100 mM NaCl, 50 mM Tris titrated to pH 8.0, 2 mM TCEP, and 2 mM DDM. AdiC protein was eluted from the column using WB added with 10 mM desthiobiotin. The AdiC-containing fractions were pooled and then concentrated with an Amicon Ultra concentrator (50K MWCO) before a further purification with a size-exclusion FPLC column (superdex 200 10/30, GE) equilibrated in WB.
Labeling and reconstitution of AdiC into nanodiscs
Request a detailed protocolPurified AdiC was mixed with POPC and the MSP2N2 membrane scaffold protein purified as described (Ritchie et al., 2009; Denisov et al., 2019) in an 1:1500:10 molar ratio of the AdiC dimer:POPC:MSP2N2. The mixture was incubated at 4°C for at least 2 hr. Nanodiscs were assembled during a dialysis of the mixture against a buffer containing 100 mM NaCl, 20 mM Tris titrated to pH 8.0, and 0.5 mM TCEP at 4°C overnight. After the dialysis, the nanodiscs were labeled with biotin using the BirA Enzyme (Avidity) following the protocol provided by the manufacturer, and then mixed with bifunctional rhodamine in a 2:1 molar ratio of the AdiC dimer:bifunctional rhodamine, a ratio to maximize the chance of only one AdiC monomer in each dimer to be labeled; this mixture was incubated at room temperature for more than 4 hr. The remaining free dye was removed first with Biobeads (SM-2, Bio-Rad) and then by size-exclusion chromatography (Superose 6 10/30 or Superderx 200 10/300, GE) equilibrated with a solution containing 100 mM NaCl and 50 mM Tris titrated to pH 8.0. As expected, the peak of the protein containing the two mutant cysteine residues, detected at 280 nm, co-migrated with that of the fluorophore dye detected at 550 nm. In contrast, no notable absorbance peak at 550 nm co-migrated with the peak of the protein without these mutant cysteine residues. Thus, there was no detectable background labeling (Figure 1—figure supplement 1A and B). The final product of nanodiscs harboring labeled AdiC was aliquoted, flash frozen in liquid nitrogen, and stored in a –80°C freezer.
Sample preparation for data collection with polarization TIRF microscope
Request a detailed protocolFor adhesion of streptavidin, one side of a cover slip is exposed to 0.01% poly-L-lysine solution for 1 hr, before being rinsed with distilled water and air dried. Prior to an experiment, a cover slip was attached, via thin transparent Scotch adhesive tapes placed on its left and right edges, to a microscopy slide, with the poly-L-lysine-coated side of the cover slip facing the bottom side of the slide. Solutions were to be placed in the space between the two pieces of glass created by the adhesive tapes that acted as a spacer. The poly-L-lysine-coated side of the cover slip was exposed to 5 mg/ml streptavidin (Promega) in this space for 15 min, and the remaining free streptavidin was washed away with a solution containing 50 mM HEPES titrated to pH 7.5 and 100 mM NaCl.
An AdiC-containing nanodisc sample was diluted to 30–100 pM, estimated from an evaluation of the absorbance of the sample at 550 nm wavelength against the extinction coefficient of bifunctional rhodamine. The diluted sample was flowed into the space between the assembled cover slip and slide. After allowing a biotin-moiety covalently linked to the N-terminus and the streptavidin-binding tags linked to the N- and C-termini in each of AdiC subunit (i.e., totaling six available attachment points per AdiC dimer) to bind to streptavidin on the cover slip, the space was thoroughly washed to remove unattached AdiC with a solution (pH 5) containing 100 mM NaCl, 100 mM dithiothreitol (DTT, Fisher, BP172), and 50 mM acetic acid titrated to pH 5, without or with arginine at a specific concentration. DDT was used to scavenge oxygen to minimize its adverse impact on the fluorophore’s emitting intensity and lifetime.
Fluorescence polarization microscope and intensity recording
Request a detailed protocolAs previously described, the fluorescence polarization microscope was built from a Nikon TIRF microscope (model Ti-E) (Lewis and Lu, 2019c). To produce an evanescent field at the surface of the sample coverslip, a 140 mW linearly polarized laser beam (532 nm) generated from a 500 mW laser (Crystalaser CL532-500-S) was directed to pass through a ¼ λ-plate, which transformed the linear polarization to circular polarization. After passing through a polarization-preserving, high numerical aperture 100× objective (Nikon Achromatic, NA = 1.49), the beam approached the coverslip with an incident angle of 68°, the so-called critical angle that leads to TIR required for the formation of an evanescent field (Axelrod et al., 1984). The emission of polarized fluorescent light from individual fluorophores excited by the evanescent field was directed to a 50:50 non-polarizing beam splitter (Thorlabs CM1-BS013) after passing through the objective (Figure 2), and then to a 540/593 nm bandpass filter (Semrock FF01-593/40-25) that prevents the propagation of excitation light. One resulting beam was further split by a glass (N-SF1) polarizing beam splitter (Thorlabs CM1 PBS251) along 0° and 90° and the other by a wire-grid polarizing beam splitter (Thorlabs WP25M-Vis) along 45° and 135°. These four emission intensity components, labeled as , , , and , were individually directed onto four designated sectors in the CCD grid of an EMCCD camera (Andor iXon Ultra 897), where the four intensities from a given fluorophore appeared in the corresponding positions of the four sections.
Here, fluorescence intensities from individual bifunctional rhodamine molecules, each attached to helix 6A in AdiC, were collected with the microscope and captured every 10 ms with an EMCCD camera at the room temperature 22°C. Following extraction of temporal information from the intensities with the changepoint analysis described below, we applied a Gaussian filter (with a corner frequency of 7.5 Hz) to all four intensity traces to reduce high-frequency noise, where the rise time was 22 ms. From these filtered intensities, we calculated angles as described below.
The experiments were performed on five separate occasions. Data collected among these separate collections are statistically comparable and were pooled together, resulting in sufficiently narrow distributions as illustrated in Figure 5. The width of the distributions reflects both technical and biological variations. Outlier data were excluded on the following basis. First, while fluorescence intensity is expected to vary among different polarization directions, the total intensity should not exhibit large variations unrelated to protein-conformational changes, such as more than one step bleaching of fluorescence. Second, for a given recording, at least 15 events are required to obtain a 95% confidence level for state identification, so any short traces with less than 15 events were excluded on the assumption that the short and long traces belong to the same distribution. All traces used contain four states and have an average length of ~4 s. There were no traces with an expected resolvable SNR that did not exhibit transitioning events. Third, for event detection and state identification, an SNR greater than 5 is necessary for the required minimum angle resolution. Thus, any set of intensity traces with this ratio less than 5 were excluded. The sample sizes were estimated on the basis of previous studies (Lewis and Lu, 2019c; Lewis and Lu, 2019a; Lewis and Lu, 2019b) to yield sufficiently small standard errors of mean to obtain accurate estimate of the mean. Practically, the error bars are comparable to the sizes of the symbols of data as illustrated in Figures 6 and 10A. The 95% confidence intervals are provided for all determined equilibrium constants in Table 1.
System parameters
Request a detailed protocolTheoretical aspects of the three-channel polarized emission system have been previously described (Fourkas, 2001), and extended to and practically implemented with four channels (Ohmachi et al., 2012; Lippert et al., 2017; Lewis and Lu, 2019c). Briefly, the intensity of a given polarized component with angle ψ, Iψ, collected by the microscope’s objective from the emission of a fluorophore, which is excited by an evenescent field generated by a circularly polarized laser beam under a TIR condition, is dependent on the fluorophore’s orientation. As such, Iψ is defined by the spherical coordinates θ and φ in a manner that
where the polarization angle ψ = 0°, 45°, 90°, or 135°. The factor fψ corrects for systematic reduction in the maximal achievable anisotropy of the light per channel ψ. The coefficients X1, X2, and X3 correct for the incomplete collection of photons by a microscope objective with collection half-angle α:
Ideally, when presented with a beam of non-polarized light, the system splits that beam into four of equal intensity. Small deviations from this theoretical equality are corrected by normalizing each intensity of a given channel (ψ°) to that of the 90° channel chosen as the reference here:
Besides those three types of system parameters, the coefficient X4 corrects for the fast diffusive motion (‘wobble’) of the probe relative to the attached protein, which is measured in terms of the half-angle δ of the wobble cone:
The parameter δ was experimentally estimated to be 22.5° in a separated macroscopic anisotropy study of the proteins as previously described (Lewis and Lu, 2019c).
Analytic solutions have been found for θ, φ and Itot from Equation 7 (Lewis and Lu, 2019c):
The angle θ was defined from 0° to 90° and φ from 0° to 180° (originally calculated as –90° to 90°) with degenerate solutions of θ or 180°- θ and φ or 180° + φ. The solutions for individual events of a given molecule, which were expected to be limited within an appropriately defined quarter sphere because of the expected small angle changes of about 10°–40°, were chosen on the basis of minimizing the variance of the distribution of the resulting angles for the events. When θ is near 90°, the noise of intensities might prevent a solution of θ, in which case we simply set the θ values to 90°. We only analyzed data from the particles with no more than a few percent of such data points.
Calibrations of the camera
Request a detailed protocolPhotons hitting the CCD chip of an EMCCD camera result in the release of individual electrons in accordance with the photoelectric effect. These electrons pass through multiple layers, at each of which they have a probability to cause the release of additional electrons, effectively amplifying the original signal with a gain (G) (Lidke et al., 2005; Heintzmann, 2016). The relationship of the N photons released over a given time interval Δt and the recorded intensity I is expressed as
As an inherent feature of EMCCD cameras, the offset is the intensity recorded by the camera when the shutter is fully closed. These parameters are estimated by analyzing the relationship between photon intensities recorded at multiple laser intensities versus the corresponding SNR. For photon count N, the standard deviation σ due to shot noise is given by the square root of N. In addition, EMCCD cameras add an additional multiplicative noise that effectively scales the σ due to shot noise by a factor of √2 that relates the SNR to N via the expression
Substituting Equation 15 into Equation 14 yields
On a plot of I·Δt versus SNR2, the slope is twice the gain, whereas the y-intercept is the offset. For our system, the value of the gain was found to be 146 and that of the offset was 220.
Furthermore, when the signal from an interested fluorophore is split from the background signal, the combined intensity signal recorded is related to N by the following relation
rearranged to
where is the number of photons underlying the background signal, which needs to be subtracted. The quantum efficiency (QE) of the Andor Ixon EMCCD camera is estimated as ~0.95 in the visible light range, meaning that 95% of photons contacting the CCD chip are detected. Therefore, the actual photon count is calculated as Ndet = N / QE, and the effective SNR is also scaled accordingly.
Detection of event transitions
Request a detailed protocolThe photon flux intensity recorded in a given channel changes when the orientation of the fluorophore is changed. To detect these transitions within the measured polarized intensities, we adopted a version of the changepoint algorithm (Chen and Gupta, 2001; Beausang et al., 2008). Such a process was based on calculating a log likelihood ratio over a period of time to determine the maximal ratio that identified the point where the perceived change of photon-release rate occurred, that is, the time at which the fluorophore transitions from one orientation to another. This method has previously been applied to analyzing the photon-arrival time captured on a continuous basis with photon-multiplier-based multi-channel recordings (Beausang et al., 2008), which was adapted for analyzing photons collected over a fixed time interval with an EMCCD camera (Lewis and Lu, 2019c).
When a camera is used as a detector, photons emerging from each channel are effectively binned over each frame. A series of consecutive k frames with a constant exposure time (Δt) is expressed as
For a given frame i, the intensity is defined by the rate of photon release:
where ni is the number of photons release within frame i. The cumulative distribution, mj, is built by adding the number of photons for the successive time frames:
If during an interval T the rate changes from I1 to I2 at the time point τ = iΔt, and the number of emitted photons prior to this change is m (Equation 21) and N – m after the change, then the likelihood ratio of a transition occurring at frame i in the log form is described by
where h equals 4, the number of emission channels in the system. We set the threshold of significance for LLR at the level that limits the false-positive events to 5% on the basis of simulation studies for the corresponding levels of SNR, and the resulting false-negative events were about 1%.
The program was started by identifying one transition over the entire trace. If a change-point X was identified, it would then search for additional transitions between the start of the trace and point X and between X and the end. This iterative search with successively shortened stretches continued until no more transitions were identified.
State identification
Request a detailed protocolFollowing the detection of intensity transition time points using the changepoint method and subsequent calculation of angles, the states of individual events, each demarcated by two consecutive transition timepoints, were identified on the basis of x,y,z values, calculated from the corresponding θ and φ values, along with of a unity value, in accordance with the transformation relations between the Cartesian and spherical coordinate systems:
Given the of a unity value, which carries no information regarding spatial orientation, for all cases, x,y,z in all cases would always be on a unit sphere and fully encode the orientation information specified by θ and φ.
The identification of the states is done by using a ‘nearest-neighbor’ method in a Cartesian coordinate system, where a given event was assigned to the closest state distribution. Closeness was determined by the minimum distance di,k between the x,y,z coordinates of the ith event and the distribution means , and , where k is the number of potential states, calculated as:
Minimization was performed by using a k-means clustering algorithm optimized with two coupled algorithms (simulated annealing plus Nelder–Mead downhill simplex)(Press et al., 2007), as we described previously in more detail (Lewis and Lu, 2019c). The resulting four conformational state distributions are denoted as C1 – C4.
States must be indexed so that they are always numerated in the same sequential order among individual molecules, a prerequisite for correctly relating the states identified here to those identified crystallographically. This indexing is based on the shortest, if not straightest, path distance between the first and last states, C1 and CM (where M = 4 in our case) as they detour through the remaining states (Figure 7—figure supplement 1B). This path length is calculated as the sum of the cumulative distances between two adjacent states, denoted as dtot:
where the number 1 indicates a chosen starting state and the distance di,j between states indicated by positions i and j in a Cartesian coordinate system is defined as:
If the states were located along a perfect line, the total path distance would equal the distance between states 1 and M, that is, d1,M, where d1,M = dtot. However, if they were not on a line, dtot would be greater than d1,M by Δd such that
Upon calculating for each of the 24 possible sequences relating the four states, the shortest path can be found on the basis of the smallest value of , which has two solutions with equal path length, that is, 1-2-3-4 versus 4-3-2-1 (Figure 7—figure supplement 1C vs. Figure 7—figure supplement 1D). The sequence of C1-C2-C3-C4 shown in Figure 7—figure supplement 1C was consistently chosen.
Transformation of the laboratory framework to a local framework
Request a detailed protocolAs explained in the text, individual molecules did not have the same orientation. Thus, a direct comparison among them requires all molecules be rotated from the laboratory frame of reference into a common local frame of reference (Figure 7—figure supplement 1A). In the local frame, the x,y-plane is defined by the vectors representing C1 (V1) and C4 (V4) and the x-axis is defined by V1. The z-axis is perpendicular to the x,y-plane. The x, y, and z-axes, represented by unit vectors, are defined by:
where
θ and φ in the local frame of reference are then calculated as:
where
Calculation of the direct angle change Ω between two states
Request a detailed protocolThe direct angle change Ωi,j between two states i and j, as represented by the vectors Vi and Vj defined above, can be calculated from the relation:
Data availability
Source-data files for all relevant figures are provided.
References
-
Total internal reflection fluorescenceAnnual Review of Biophysics and Bioengineering 13:247–268.https://doi.org/10.1146/annurev.bb.13.060184.001335
-
Secondary active transportersSub-Cellular Biochemistry 92:275–299.https://doi.org/10.1007/978-3-030-18768-2_9
-
Protein length in eukaryotic and prokaryotic proteomesNucleic Acids Research 33:3390–3400.https://doi.org/10.1093/nar/gki615
-
Projection structure of a member of the amino acid/polyamine/organocation transporter superfamilyThe Journal of Biological Chemistry 283:33240–33248.https://doi.org/10.1074/jbc.M806917200
-
On change point detection and estimationCommunications in Statistics - Simulation and Computation 30:665–697.https://doi.org/10.1081/SAC-100105085
-
A homobifunctional rhodamine for labeling proteins with defined orientations of a fluorophoreBioconjugate Chemistry 9:160–167.https://doi.org/10.1021/bc970174e
-
Nanodiscs as a new tool to examine lipid-protein interactionsMethods in Molecular Biology 2003:645–671.https://doi.org/10.1007/978-1-4939-9512-7_25
-
A bacterial arginine-agmatine exchange transporter involved in extreme acid resistanceThe Journal of Biological Chemistry 282:176–182.https://doi.org/10.1074/jbc.M610075200
-
Escherichia coli acid resistance: tales of an amateur acidophileNature Reviews. Microbiology 2:898–907.https://doi.org/10.1038/nrmicro1021
-
Hindered rotational diffusion and rotational jumps of single moleculesPhysical Review Letters 80:2093–2096.https://doi.org/10.1103/PhysRevLett.80.2093
-
Arginine-Agmatine antiporter in extreme acid resistance in Escherichia coliJournal of Bacteriology 185:6556–6561.https://doi.org/10.1128/JB.185.22.6556-6561.2003
-
Function and regulation of acid resistance antiportersThe Journal of Membrane Biology 252:465–481.https://doi.org/10.1007/s00232-019-00073-6
-
Energetics of ångström-scale conformational changes in an RCK domain of the mthk K+ channelNature Structural & Molecular Biology 26:808–815.https://doi.org/10.1038/s41594-019-0275-1
-
Integrating spatiotemporal features of a ligand-regulated, multi-state allosteric proteinNature Structural & Molecular Biology 26:816–822.https://doi.org/10.1038/s41594-019-0276-0
-
Resolution of ångström-scale protein conformational changes by analyzing fluorescence anisotropyNature Structural & Molecular Biology 26:802–807.https://doi.org/10.1038/s41594-019-0274-2
-
The role of photon statistics in fluorescence anisotropy imagingIEEE Transactions on Image Processing 14:1237–1245.https://doi.org/10.1109/tip.2005.852458
-
On the nature of allosteric transitions: a plausible modelJournal of Molecular Biology 12:88–118.https://doi.org/10.1016/s0022-2836(65)80285-6
-
Activation by adenosine triphosphate in the phosphorylation kinetics of sodium and potassium ion transport adenosine triphosphataseThe Journal of Biological Chemistry 247:6530–6540.
-
Chapter 11-reconstitution of membrane proteins in phospholipid bilayer nanodiscsMethods in Enzymology 464:211–231.https://doi.org/10.1016/S0076-6879(09)64011-8
-
Rotational motions of macro-molecules by single-molecule fluorescence microscopyAccounts of Chemical Research 38:583–593.https://doi.org/10.1021/ar040137k
-
Fluorescence energy transfer as a spectroscopic rulerAnnual Review of Biochemistry 47:819–846.https://doi.org/10.1146/annurev.bi.47.070178.004131
Article and author information
Author details
Funding
National Institute of Diabetes and Digestive and Kidney Diseases (DK125521)
- Zhe Lu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This study was supported by the grant DK125521 from the National Institute of Diabetes and Digestive and Kidney Diseases.
Copyright
© 2023, Zhou et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 485
- views
-
- 59
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Neuroscience
The assembly and maintenance of neural circuits is crucial for proper brain function. Although the assembly of brain circuits has been extensively studied, much less is understood about the mechanisms controlling their maintenance as animals mature. In the olfactory system, the axons of olfactory sensory neurons (OSNs) expressing the same odor receptor converge into discrete synaptic structures of the olfactory bulb (OB) called glomeruli, forming a stereotypic odor map. The OB projection neurons, called mitral and tufted cells (M/Ts), have a single dendrite that branches into a single glomerulus, where they make synapses with OSNs. We used a genetic method to progressively eliminate the vast majority of M/T cells in early postnatal mice, and observed that the assembly of the OB bulb circuits proceeded normally. However, as the animals became adults the apical dendrite of remaining M/Ts grew multiple branches that innervated several glomeruli, and OSNs expressing single odor receptors projected their axons into multiple glomeruli, disrupting the olfactory sensory map. Moreover, ablating the M/Ts in adult animals also resulted in similar structural changes in the projections of remaining M/Ts and axons from OSNs. Interestingly, the ability of these mice to detect odors was relatively preserved despite only having 1–5% of projection neurons transmitting odorant information to the brain, and having highly disrupted circuits in the OB. These results indicate that a reduced number of projection neurons does not affect the normal assembly of the olfactory circuit, but induces structural instability of the olfactory circuitry of adult animals.
-
- Neuroscience
Specialized chemosensory signals elicit innate social behaviors in individuals of several vertebrate species, a process that is mediated via the accessory olfactory system (AOS). The AOS comprising the peripheral sensory vomeronasal organ has evolved elaborate molecular and cellular mechanisms to detect chemo signals. To gain insight into the cell types, developmental gene expression patterns, and functional differences amongst neurons, we performed single-cell transcriptomics of the mouse vomeronasal sensory epithelium. Our analysis reveals diverse cell types with gene expression patterns specific to each, which we made available as a searchable web resource accessed from https://www.scvnoexplorer.com. Pseudo-time developmental analysis indicates that neurons originating from common progenitors diverge in their gene expression during maturation with transient and persistent transcription factor expression at critical branch points. Comparative analysis across two of the major neuronal subtypes that express divergent GPCR families and the G-protein subunits Gnai2 or Gnao1, reveals significantly higher expression of endoplasmic reticulum (ER) associated genes within Gnao1 neurons. In addition, differences in ER content and prevalence of cubic membrane ER ultrastructure revealed by electron microscopy, indicate fundamental differences in ER function.