Structural and functional insights of the human peroxisomal ABC transporter ALDP
Abstract
Adrenoleukodystrophy protein (ALDP) is responsible for the transport of very-long-chain fatty acids (VLCFAs) and corresponding CoA-esters across the peroxisomal membrane. Dysfunction of ALDP leads to peroxisomal metabolic disorder exemplified by X-linked adrenoleukodystrophy (ALD). Hundreds of ALD-causing mutations have been identified on ALDP. However, the pathogenic mechanisms of these mutations are restricted to clinical description due to limited structural and biochemical characterization. Here we report the cryo-electron microscopy structure of human ALDP with nominal resolution at 3.4 Å. ALDP exhibits a cytosolic-facing conformation. Compared to other lipid ATP-binding cassette transporters, ALDP has two substrate binding cavities formed by the transmembrane domains. Such structural organization may be suitable for the coordination of VLCFAs. Based on the structure, we performed integrative analysis of the cellular trafficking, protein thermostability, ATP hydrolysis, and the transport activity of representative mutations. These results provide a framework for understanding the working mechanism of ALDP and pathogenic roles of disease-associated mutations.
Editor's evaluation
The adrenoleukodystrophy protein (ALDP) or ABCD1 is an ABC transporter that participates in the transport of free very long-chain fatty acids and their CoA esters across the peroxisomal membrane. By determining the cryo EM structure of human ABCD1 the study represents a valuable insight into its transport mechanism and the mechanistic basis for mutations causing the severe neurodegenerative disorder, X-linked adrenoleukodystrophy. The structure and functional studies of disease-causing mutations are solid and will appeal to the transporter and medical genetics communities.
https://doi.org/10.7554/eLife.75039.sa0Introduction
X-linked adrenoleukodystrophy (ALD) is an inherited disease characterized by progressive demyelination of the central nervous system (Berger et al., 2014). Loss of myelin slows down the transmission of nerve impulses and triggers neuroinflammation (Ferrer et al., 2010; Kettwig et al., 2021; Singh et al., 2009). Pathogenesis of ALD is tightly associated with mutations on the adrenoleukodystrophy protein (ALDP), which transports very-long-chain fatty acid-CoAs (VLCFA-CoAs) from cytosol into the peroxisome (Figure 1; Mosser et al., 1993). Over 900 disease-derived mutations have been identified on ALDP (https://adrenoleukodystrophy.info/mutations-biochemistry/mutations-biochemistry). ALDP, also known as ABCD1, belongs to the ATP-binding cassette sub-family D. The other three members are ALDP-related protein (ALDRP/ABCD2; Lombard-Platet et al., 1996), PMP70/ABCD3 (Kamijo et al., 1990), and a cobalamin transporter ABCD4 (Coelho et al., 2012). ABCD1–3 are distributed on peroxisomes, while ABCD4 is localized on lysosomes (Xu et al., 2019). ALDP and ALDRP transport VLCFA with different specificities (van Roermund et al., 2011). ALDRP may compensate for the function loss of ALDP (Holzinger et al., 1997; Liu et al., 1999; Lombard-Platet et al., 1996). PMP70 was reported to transport LCFA-CoA (Ranea-Robles et al., 2021; van Roermund et al., 2014).
Structure-based functional characterization of ALDP is eagerly required for understanding the VLCFA transport mechanism and the pathogenic roles of ALD-derived mutations. To address these questions, we first performed ATP hydrolysis assay using purified ALDP. In the presence of C22:0-CoA, ALDP exhibits robust ATP hydrolysis with the Vmax of ~193 ± 8.2 mol Pi min−1 mol−1 protein and Km value of ~0.17 µM C22:0-CoA (Figure 1B and C). The structure of human ALDP reveals an assembly of two identical subunits that exhibit a domain-swapped arrangement (Figure 1D and Figure 1—figure supplements 1–2 and Figure 1—source data 1). Each subunit contains six transmembrane helices (TMs; Figure 1D and Figure 1—figure supplement 2). A short helix located at the peroxisomal side (extracellular helix [EH]) is a featured structural element of ALDP based on sequence alignment (Figure 1—figure supplement 4). Two coiled-coil domain-like densities at the C-terminus were observed (C-terminal helix [CH]). The local resolution of these densities is estimated to be ~4.0 Å, which is not accurate to build an atomic model.
During submission of this study, four other groups independently reported ALDP structures in different states (Chen et al., 2022; Le et al., 2022; Wang et al., 2022; Li et al., 2021). Our structure is the ligand-free state as no substrates were supplemented during purification. We observed two unknown densities at the crevice of TMDs, which were also reported in the ligand-free structure by Chen et al., 2022. When superimposed to the oleoyl-CoA- or C22:0-CoA-bound structures, the unknown densities partially overlap with the substrates (Figure 1—figure supplement 3). Such observation indicates a probable extrusion of the unknown densities during substrate loading. Structure comparison of our structure with the C22:0-CoA-bound and ATP-bound structure reveals dramatic conformational changes during the transport cycle. For instance, the distance between TM4 and TM6 at the juxta-membrane region changes from 23.6 Å to 13.1 Å upon substrate binding and further gets closed to 5.7 Å in the presence of ATP. The nucleotide-binding domains (NBDs) undergo continuously spatial rearrangement. Consequently, the TM4-6 and TM1-2 bend to open the peroxisomal side exit for substrate (Figure 2). Conformational change of TMDs is accompanied with the movement of NBDs and the disappearance of EH or CH.
Structure determination of ALDP allows classification of ALD-associated mutations. 970 ALD-associated mutations in ALDP affecting 232 residues have been reported. Among these 232 residues, 88 harbor mutations to two or more types of amino acids, which we defined as hotspot residues. Based on the structural mapping, the 88 hotspot residues can be roughly classified into 3 groups (Figure 3B). The first group of 10 residues lines along the substrate-binding cavity. These residues likely play a role in substrate coordination and delivery. The second group of 35 residues is located in other region of the TMDs, which undergoes conformational changes. The third group of 43 residues is located on the NBDs that may have influence on ATP hydrolysis (Figure 3B and Figure 3—figure supplement 1).
To examine how ALD-associated mutations affect the function of ALDP, we selected those reported to have normal expressed levels but dysfunction for integrative analysis (Coll et al., 2005; Feigenbaum et al., 1996; Guimarães et al., 2002). These mutations include W339R, S342P, G343V, A396T, Q544R, and T693M (Imamura et al., 1997; Kemp et al., 2001; Lan et al., 2011; Liu et al., 2022; Pan et al., 2005; Takahashi et al., 2007; Watkins et al., 1995; Wichers et al., 1999). In addition, to evaluate the role of different structure elements, we designed four constructs based on structural changes. The S164/Y310 residues on TMDs, anchored by a hydrogen bond, both move quite large distance during the transport cycle but remain relatively static distance between themselves. The E291/R518 residues are mapped to the interface between TMD and NBD. The side chain of E291 stretches into the NBD of another protomer, while R518 binds to the α-phosphate of ATP. These two pairs of residues represent two kinds of dimeric interfaces. We mutated these residues into alanine to examine how disturbance of the dimeric assembly affects the activity. Another two constructs are truncation of EH (Δ364–374) or CH (Δ683–745).
We first measured the VLCFA-CoA transport. Different mutations are individually expressed in HEK293 cells. The cytosolic C22:0-CoA, C24:0-CoA, and C26:0-CoA were measured by quantitative mass spectrometry (Wang et al., 2019), and the transport ability was calculated (Figure 3C). Because the substrate transport is coupled with ATP hydrolysis, we also measured the ATPase activity of these mutants using purified proteins (Figure 3C). To exclude the possibility that the mutants may alter the cellular localization and thermostability compared to wild-type (WT), we investigated their subcellular localization and measured the Tm values of purified proteins. Compared with A95D, a disease-causing mutation known to not colocalize with peroxisome and reduce the C24:0 β-oxidation (Morita et al., 2013), all the selected mutants show puncta distribution and colocalize with catalase, suggesting their normal trafficking to the peroxisomes (Figure 3—figure supplement 2). The similar Tm values of purified mutant proteins rule out the instability-caused functional abrogation (Figure 3—figure supplement 3).
The transport activity and the ATP hydrolysis of the mutants are normalized to WT. All six disease-derived mutants exhibit decreased transport toward three species of VLCFA-CoA. Five mutants show abrogated ATP hydrolysis except for the mutation S342P. The ALD mutation T693M is mapped on the CH. Because the CH cannot be observed in the ATP-bound state, to evaluate the potential pathogenic role of the T693M mutant, we generated the homomeric model of ALDP by AlphaFold (hereafter as AF-model) (Evans et al., 2022; Jumper et al., 2021). The predicted model of ALDP displays similar conformation to the ATP-bound state. Albeit the part of CH in AF-model needs further validation, the prediction may provide potential clues for interpreting the T693M mutation at the current stage. In the AF-model, T693 is located near W664, which is surrounded by ALD mutants on the NBD. The distance between the side chains of T693 and W664 can reach in 3.4 Å (Figure 3—figure supplement 4). Replacement of Thr by Met residue with large side chain may affect the neighboring residues near W664 to disturb the ATP hydrolysis. For the four designed constructs, alanine substitution of S164/Y310 and E291/R518 has negligible influence on the ATP hydrolysis. By contrast, the lipid transport decreased because destabilization of the interaction between these residues may disturb the conformational changes during substrate translocation. As one featured structural element, deletion of the EH (Δ364–374) may generate structural constraints between TM5 and TM6 during the conformational change of ALDP (Chen et al., 2022; Le et al., 2022; Wang et al., 2022). Truncation of CH (Δ683–745) does not change the assembly of NBD but affects the conformational changes between the two protomers, thus the ATP hydrolysis level slightly decreased with abrogated lipid transport.
To gain insights into the transport of different lipids by ABC transporters, we compared ALDP with the phospholipid transporter ABCB4 (Nosol et al., 2021; Olsen et al., 2020). Typically, ALDP harbors two substrate-binding pockets that expand perpendicularly to the membrane, but ABCB4 binds one phosphatidylcholine (PC) (Figure 3—figure supplement 5). The dimension of one binding pocket is ~1695 Å3 for C22:0-CoA and ~1236 Å3 for PC calculated by ProteinPlus (Schöning-Stierand et al., 2020). The different formulas between VLCFA-CoA and PC may provide clues for the different cavities. Regardless of the flexibility, the length of C22:0-CoA (~55 Å) is larger than PC (~40 Å) (Figure 3—figure supplement 5B). The PC molecule has two acyl tails that can be coordinated by ABCB4. The VLCFAs have only one hydrophobic tail. More residues in a larger binding pocket of ALDP may help in stabilizing the flexible long tail of VLCFA-CoA. The EH and CH facilitate conformational changes to expand the binding pocket. Taken together, the integrative characterization of ALDP provides a framework to illustrate the molecular basis of VLCFA transport and the function of ALD-derived mutations.
Methods
Protein expression and purification
The coding sequence of human ALDP was subcloned into the pFastBac1 vector with a N-terminal Flag tag. The recombinant ALDP was expressed using the Bac-to-bac system (Invitrogen). Briefly, bacmids were generated in DH10Bac competent cells. The resulting baculoviruses were amplified in Sf9 insect cells (Invitrogen). ALDP was overexpressed in Sf9 insect cells grown in the serum-free medium (Gibco). 60 hr after P3 virus infection, transfected cell pellet was collected and homogenized in the buffer containing 25 mM Tris-HCl, pH 7.4, 150 mM NaCl supplemented with 1 mM phenylmethylsulfonyl fluoride, 1.3 μg ml−1 aprotinin, 0.7 μg ml−1 pepstatin, and 5 μg ml−1 leupeptin. After brief sonication, the suspension was supplemented with 2,2-didecylpropane-1,3-bis-β-D-maltopyranoside (LMNG, Anatrace) to a final concentration of 1% (w/v) and cholesteryl hemisuccinate Tris salt (CHS, Anatrace) to 0.1% (w/v). After incubation at 4°C for 2 hr, the mixture was centrifuged at 150,000 g for 30 min. The supernatant was mixed with anti-Flag M2 affinity gel (Sigma). The gel was rinsed with 25 mM Tris-HCl, pH 7.4, 150 mM NaCl, 0.005% LMNG, and 0.0005% CHS. The target protein was eluted with wash buffer plus 200 μg ml−1 flag peptide and further purified through gel-filtration chromatography in the buffer containing 25 mM Tris-HCl, pH 7.4, 150 mM NaCl, and 0.06% digitonin. Peak fractions were concentrated for cryo-EM analysis.
Preparation of the cryo-EM samples
Aliquots of 4 µl ALDP was dropped onto glow-discharged holey carbon grids (Quantifoil Au R1.2/1.3, 300 mesh). The grids were blotted for 3 s and flash-frozen in liquid ethane using Vitrobot Mark IV (FEI). The sample was imaged on an FEI 300 kV Titan Krios transmission electron microscope equipped with a Cs corrector and Gatan GIF Quantum energy filter (slit width 20 eV), recorded by a Gatan K2 Summit detector with a nominal magnification of ×640,000. A series of defocus values from –1.5 to –1.8 µm was used during data collection. Each image was dose-fractionated to 32 frames with a total electron dose of ~50 e− Å−2 and a total exposure time of 5.6 s. AutoEMation II (developed by J. Lei) (Lei and Frank, 2005) was used for fully automated data collection. All stacks were motion-corrected using MotionCor2 with a binning factor of 2 (Zheng et al., 2017), resulting in a pixel size of 1.0979 Å. The defocus values were estimated using Gctf (Zhang, 2016), and dose weighing was performed concurrently (Grant and Grigorieff, 2015).
Cryo-EM data processing and model building
In total, 1,259,972 particles were auto-picked from these 1922 movie stacks using Gautomatch (developed by Kai Zhang, http://www.mrc-lmb.cam.ac.uk/kzhang/Gautomatch) (Figure 1—figure supplement 1). After two-dimensional classification, 472,168 particles were selected to generate the initial model with a mask diameter of 200 Å in C1 symmetry. To avoid model bias, the initial model was low-pass filtered for the following three-dimensional (3D) classification. The 472,168 particles were subjected to 50 iterations of global angular search 3D classification. Each of the 50 iterations has one class and a step size of 7.5°. For each of the last five iterations (iterations 46–50) of the global search, the local angular search 3D classification was executed with a class number of 3, a step size of 3.75°, and a local search range of 15°. A total of 289,321 particles were selected from the local angular searching 3D classification after removing the redundant particles and particles from bad classes (Figure 1—figure supplement 1). The selected particles were subjected to 3D autorefinement, yielding a density map with average resolution at 3.6 Å. After multi-reference 3D classification, 135,662 particles were selected for autorefinement with C2 symmetry, resulting a map at 3.4 Å resolution. The structure of ABCD4 was used as a template for the model building of ALDP. The structure was refined in real space using PHENIX with secondary structure and geometry restraints (Adams et al., 2002). The atomic model was manually improved using COOT (Emsley and Cowtan, 2004).
Determination of cytosolic very-long-chain fatty acids
Sample preparation and liquid chromatography-mass spectrometry (LC-MS/MS) conditions as described (Wang et al., 2019). Chromatographic separation was achieved using Q Exactive HF LC-MS/MS system and Phenomenex Luna 5 µm C5 column (i.d. 100×2.0 mm). Mobile phase A contained HPLC-grade H2O-ACN 40/60 (v/v) with 10 mM ammonium acetate, and mobile B was isopropanol-ACN 90/10 (v/v). Both WT ALDP and mutant constructs were subcloned into a pCAG vector and expressed in HEK293F cells (Invitrogen). Transfected cells were harvested after 60 hr. 2.0×107 cells were collected and washed by PBS. 1 mL ACN and 10 μL internal standard mix solution were added to samples or quality control sample (expressed by empty vector). After sonication and centrifugation, 50 μL hydrochloric acid (12 mol/L) was added into the supernatant then VLCFA was hydrolyzed at 70°C for 1 hr. Each sample was extracted with 2 ml hexane, the dried residue was reconstituted in 100 μL methanol. 5 μL aliquots of the reconstitutes were loaded onto the LC–MS/MS system for analysis.
Immunofluorescence cytochemistry
Hela cells fixed in 4% Paraformaldehyde (20 min) were permeabilized with 0.1% Triton X-100 in PBS (30 min, room temperature), blocked with SuperBlock (30 min), before incubation with rabbit anti-catalase antibody (1:500; Abcam, ab16731). Immunoreactivity was visualized by AlexaFluor-488-conjugated goat anti-rabbit IgG (1:1000; 1 hr at room temperature).
ATPase activity assay
For determining ATP hydrolysis capacity, wtABCD1 and mutants were deleted N-terminal 54 residues to improve stability. The malachite green ATPase assay was chosen to exhibit the ABCD1 ATPase activity in the presence of C22:0-CoA (behenoyl coenzyme A, ammonium salt) as previously described (Baykov et al., 1988). Purified ABCD1 in LMNG/CHS micelles was pre-incubated on ice for 10 min in 20 mM Tris-HCl, pH 7.5, 50 mM KCl, 1 mM DTT, 2 mM MgCl2 with different concentrations of C22:0-CoA. The 100 µl reaction started by supplementing with 2 mM ATP and then carried out at 37°C for 30 min, after that terminated by adding 25 µl fresh Gold-mix solution (1 mM malachite green, 1.2% ammonium molybdate, and 0.15% Tween 20). The mixture was incubated at RT for 30 min before being detected absorbance in 96-well micro-plate at 630 nm. Statistical analysis was performed using GraphPad Prism 7.
The thermal stability assay
The thermal stability assay was performed by measuring intrinsic tryptophan fluorescence and light back-scattering for fragment screening by a Prometheus NT.48 device (NanoTemper Technologies GmbH, Munich, Germany) (Ahmad et al., 2021). wtABCD1 and mutations were purified in LMNG/CHS and concentrated to 0.5 mg/ml, after centrifuged at 14,000 × g for 15 min at 4°C, supernatant was loaded with standard capillaries (NanoTemper Technologies GmbH, Munich, Germany; Cat# PR-C002) into the Prometheus device and subjected to a linear thermal ramp (1°C/min, from 25 to 95°C) and collected fluorescence at 350 nm and 330 nm. Unfolding transition midpoints were determined automatically from the first derivative of the fluorescence ratio (F350/F330). Data analysis was performed using PR.Therm Control (NanoTemper Technologies GmbH) and GraphPad Prism 7.
Data availability
Cryo-EM data have been deposited in PDB under the accession code 7VR1. All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figure 1 and Figure 3—figure supplement 3.
-
RCSB Protein Data BankID 7VR1. Cryo-EM structure of the ATP-binding cassette sub-family D member 1 from Homo sapiens.
-
The ABCD1 Variant DatabaseID abcd1. The ABCD1 Variant Database.
References
-
PHENIX: building new software for automated crystallographic structure determinationActa Crystallographica. Section D, Biological Crystallography 58:1948–1954.https://doi.org/10.1107/s0907444902016657
-
Nano-differential scanning fluorimetry for screening in fragment-based lead discoveryJournal of Visualized Experiments 1:e469.https://doi.org/10.3791/62469
-
Mutations in ABCD4 cause a new inborn error of vitamin B12 metabolismNature Genetics 44:1152–1155.https://doi.org/10.1038/ng.2386
-
Coot: model-building tools for molecular graphicsActa Crystallographica. Section D, Biological Crystallography 60:2126–2132.https://doi.org/10.1107/S0907444904019158
-
Mutational and protein analysis of patients and heterozygous women with X-linked adrenoleukodystrophyAmerican Journal of Human Genetics 58:1135–1144.
-
Molecular characterization of 21 X-ALD Portuguese families: identification of eight novel mutations in the ABCD1 geneMolecular Genetics and Metabolism 76:62–67.https://doi.org/10.1016/s1096-7192(02)00023-9
-
Cdna cloning and mRNA expression of the human adrenoleukodystrophy related protein (ALDRP), a peroxisomal ABC transporterBiochemical and Biophysical Research Communications 239:261–264.https://doi.org/10.1006/bbrc.1997.7391
-
The 70-kDa peroxisomal membrane protein is a member of the MDR (P-glycoprotein) -related ATP-binding protein superfamilyThe Journal of Biological Chemistry 265:4534–4540.https://doi.org/10.1016/S0021-9258(19)39595-X
-
Targeted metabolomics revealed changes in phospholipids during the development of neuroinflammation in abcd1tm1kds mice and X-linked adrenoleukodystrophy patientsJournal of Inherited Metabolic Disease 44:1174–1185.https://doi.org/10.1002/jimd.12389
-
Molecular diagnosis of X-linked adrenoleukodystrophy: experience from a clinical genetic laboratory in mainland China with report of 13 novel mutationsClinica Chimica Acta; International Journal of Clinical Chemistry 412:970–974.https://doi.org/10.1016/j.cca.2011.01.036
-
Homo- and heterodimerization of peroxisomal ATP-binding cassette half-transportersThe Journal of Biological Chemistry 274:32738–32743.https://doi.org/10.1074/jbc.274.46.32738
-
Genetic analysis and prenatal diagnosis of 76 Chinese families with X-linked adrenoleukodystrophyMolecular Genetics & Genomic Medicine 10:e1844.https://doi.org/10.1002/mgg3.1844
-
Structure of the human lipid exporter Abcb4 in a lipid environmentNature Structural & Molecular Biology 27:62–70.https://doi.org/10.1038/s41594-019-0354-3
-
The peroxisomal transporter ABCD3 plays a major role in hepatic dicarboxylic fatty acid metabolism and lipid homeostasisJournal of Inherited Metabolic Disease 44:1419–1433.https://doi.org/10.1002/jimd.12440
-
ProteinsPlus: interactive analysis of protein-ligand binding interfacesNucleic Acids Research 48:W48–W53.https://doi.org/10.1093/nar/gkaa235
-
Differential substrate specificities of human ABCD1 and ABCD2 in peroxisomal fatty acid β-oxidationBiochimica et Biophysica Acta 1811:148–152.https://doi.org/10.1016/j.bbalip.2010.11.010
-
A role for the human peroxisomal half-transporter ABCD3 in the oxidation of dicarboxylic acidsBiochimica et Biophysica Acta 1841:563–568.https://doi.org/10.1016/j.bbalip.2013.12.001
-
Rapid liquid chromatography-tandem mass spectrometry to determine very-long-chain fatty acids in human and to establish reference intervals for the Chinese populationClinica Chimica Acta; International Journal of Clinical Chemistry 495:185–190.https://doi.org/10.1016/j.cca.2019.04.058
-
Altered expression of ALDP in X-linked adrenoleukodystrophyAmerican Journal of Human Genetics 57:292–301.
-
Cryo-Em structure of human lysosomal cobalamin exporter ABCD4Cell Research 29:1039–1041.https://doi.org/10.1038/s41422-019-0222-z
-
Gctf: real-time CTF determination and correctionJournal of Structural Biology 193:1–12.https://doi.org/10.1016/j.jsb.2015.11.003
Article and author information
Author details
Funding
National Natural Science Foundation of China (3217110084)
- Guanghui Yang
National Natural Science Foundation of China (32130081)
- Guanghui Yang
National Natural Science Foundation of China (32100759)
- Zhengxin Ying
Chinese Universities Scientific Fund (15050004)
- Guanghui Yang
Chinese Universities Scientific Fund (15051002)
- Guanghui Yang
Chinese Universities Scientific Fund (15050017)
- Guanghui Yang
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank the Tsinghua University Branch of China National Center for Protein Sciences (Beijing) for the cryo-EM facility and the computational facility support, Dr. Xiaomin Li, Dr. Fan Yang and Tao Liu for technical support in EM data acquisition. We thank Mr. Wei Gao and Ms. Yang Yue at Mo-In Biotechnology (Beijing) Co., Ltd for the data acquisition and analysis of the protein thermostability. This work was supported by National Natural Science Foundation of China (3217110084; 32130081; 32100759); Chinese Universities Scientific Fund (15050004; 15050017; 15051002); and Young Elite Scientists Sponsorship Program by China Association for Science and Technology.
Copyright
© 2022, Jia, Zhang 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
-
- 1,338
- views
-
- 282
- downloads
-
- 3
- 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
-
- Structural Biology and Molecular Biophysics
The functional effects of an RNA can arise from complex three-dimensional folds known as tertiary structures. However, predicting the tertiary structure of an RNA and whether an RNA adopts distinct tertiary conformations remains challenging. To address this, we developed BASH MaP, a single-molecule dimethyl sulfate (DMS) footprinting method and DAGGER, a computational pipeline, to identify alternative tertiary structures adopted by different molecules of RNA. BASH MaP utilizes potassium borohydride to reveal the chemical accessibility of the N7 position of guanosine, a key mediator of tertiary structures. We used BASH MaP to identify diverse conformational states and dynamics of RNA G-quadruplexes, an important RNA tertiary motif, in vitro and in cells. BASH MaP and DAGGER analysis of the fluorogenic aptamer Spinach reveals that it adopts alternative tertiary conformations which determine its fluorescence states. BASH MaP thus provides an approach for structural analysis of RNA by revealing previously undetectable tertiary structures.
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Riboswitches represent a class of non-coding RNA that possess the unique ability to specifically bind ligands and, in response, regulate gene expression. A recent report unveiled a type of riboswitch, known as the guanidine-IV riboswitch, which responds to guanidine levels to regulate downstream genetic transcription. However, the precise molecular mechanism through which the riboswitch senses its target ligand and undergoes conformational changes remain elusive. This gap in understanding has impeded the potential applications of this riboswitch. To bridge this knowledge gap, our study investigated the conformational dynamics of the guanidine-IV riboswitch RNA upon ligand binding. We employed single-molecule fluorescence resonance energy transfer (smFRET) to dissect the behaviors of the aptamer, terminator, and full-length riboswitch. Our findings indicated that the aptamer portion exhibited higher sensitivity to guanidine compared to the terminator and full-length constructs. Additionally, we utilized Position-specific Labelling of RNA (PLOR) combined with smFRET to observe, at the single-nucleotide and single-molecule level, the structural transitions experienced by the guanidine-IV riboswitch during transcription. Notably, we discovered that the influence of guanidine on the riboswitch RNA’s conformations was significantly reduced after the transcription of 88 nucleotides. Furthermore, we proposed a folding model for the guanidine-IV riboswitch in the absence and presence of guanidine, thereby providing insights into its ligand-response mechanism.