Atomistic Mechanism of Large-Scale Conformational Transition in a Heterodimeric ABC Exporter

Coupling the role of lipids to the conformational dynamics of the ABC transporter P-glycoprotein

The ATP-binding cassette transporter P-glycoprotein (P-gp) is a multidrug efflux pump that is overexpressed in a variety of cancers and associated with the drug-resistance phenomenon. P-gp structures were previously determined in detergent and in nanodiscs, in which different transmembrane helix conformations were found, "straight" and "kinked," respectively, indicating a possible role of the lipid environment on the P-gp structural ensemble. Here, we investigate the dynamic conformational ensembles and protein-lipid interactions of two human P-gp inward-open conformers, straight and kinked, employing all-atom molecular dynamics (MD) simulations in asymmetric multicomponent lipid bilayers that mimic the highly specialized hepatocyte membrane in which P-gp is expressed. The two conformers are found to differ in terms of the accessibility of the substrate cavity. The MD simulations show how cholesterol and different lipid species wedge, snorkel, and partially enter into the cavity of the straight P-gp conformer solved in detergent. However, access to the cavity of the kinked P-gp conformer solved in nanodiscs is restricted. Furthermore, the volume and dynamic fluctuations of the substrate cavity largely differ between the two P-gp conformers and are modulated by the presence (or absence) of cholesterol in the membrane and/or of ATP. From the mechanistic perspective, the findings indicate that the straight conformer likely precedes the kinked conformer in the functional working cycle of P-gp, with the latter conformation representing a post substrate-bound state. The inaccessibility of the main transmembrane cavity in the kinked conformer might be crucial in preventing substrate disengagement and transport withdrawal. Remarkably, in our unbiased MD simulations, one transmembrane helix (TM10) of the straight conformer underwent a spontaneous transition to a kinked conformation, underlining the relevance of both conformations in a native phospholipid environment and revealing structural descriptors defining the transition between the two P-gp conformers.

Membrane-assisted tariquidar access and binding mechanisms of human ATP-binding cassette transporter P-glycoprotein

The human multidrug transporter P-glycoprotein (P-gp) is physiologically essential and of key relevance to biomedicine. Recent structural studies have shed light on the mode of inhibition of the third-generation inhibitors for human P-gp, but the molecular mechanism by which these inhibitors enter the transmembrane sites remains poorly understood. In this study, we utilized all-atom molecular dynamics (MD) simulations to characterize human P-gp dynamics under a potent inhibitor, tariquidar, bound condition, as well as the atomic-level binding pathways in an explicit membrane/water environment. Extensive unbiased simulations show that human P-gp remains relatively stable in tariquidar-free and bound states, while exhibiting a high dynamic binding mode at either the drug-binding pocket or the regulatory site. Free energy estimations by partial nudged elastic band (PNEB) simulations and Molecular Mechanics Generalized Born Surface Area (MM/GBSA) method identify two energetically favorable binding pathways originating from the cytoplasmic gate with an extended tariquidar conformation. Interestingly, free tariquidar in the lipid membrane predominantly adopts extended conformations similar to those observed at the regulatory site. These results suggest that membrane lipids may preconfigure tariquidar into an active ligand conformation for efficient binding to the regulatory site. However, due to its conformational plasticity, tariquidar ultimately moves toward the drug-binding pocket in both pathways, explaining how it acts as a substrate at low concentrations. Our molecular findings propose a membrane-assisted mechanism for the access and binding of the third-generation inhibitors to the binding sites of human P-gp, and offer deeper insights into the molecule design of more potent inhibitors against P-gp-mediated drug resistance.

Asymmetric conformations and lipid interactions shape the ATP-coupled cycle of a heterodimeric ABC transporter

Here we used cryo-electron microscopy (cryo-EM), double electron-electron resonance spectroscopy (DEER), and molecular dynamics (MD) simulations, to capture and characterize ATP- and substrate-bound inward-facing (IF) and occluded (OC) conformational states of the heterodimeric ATP binding cassette (ABC) multidrug exporter BmrCD in lipid nanodiscs. Supported by DEER analysis, the structures reveal that ATP-powered isomerization entails changes in the relative symmetry of the BmrC and BmrD subunits that propagates from the transmembrane domain to the nucleotide binding domain. The structures uncover asymmetric substrate and Mg2+ binding which we hypothesize are required for triggering ATP hydrolysis preferentially in one of the nucleotide-binding sites. MD simulations demonstrate that multiple lipid molecules differentially bind the IF versus the OC conformation thus establishing that lipid interactions modulate BmrCD energy landscape. Our findings are framed in a model that highlights the role of asymmetric conformations in the ATP-coupled transport with general implications to the mechanism of ABC transporters.

Computational and structural insights into the pre- and post-hydrolysis states of bovine multidrug resistance-associated protein 1

ATP-binding cassette C-family drug membrane transporters play an important role in local pharmacokinetics, that is, drug concentration in cellular compartments. From the structural point of view, only the bovine ortholog of the multidrug resistance-associated protein 1 (bMRP1) has been resolved. We here used μs-scaled molecular dynamics simulations to investigate the structure and dynamics of the bovine multidrug resistance-associated protein 1 in pre- and post-hydrolysis functional states. The present work aims to examine the slight but likely relevant structural differences between pre- and post-hydrolysis states of outward-facing conformations as well as the interactions between the multidrug resistance-associated protein 1 and the surrounding lipid bilayer. Global conformational dynamics show unfavourable extracellular opening associated with nucleotide-binding domain dimerization indicating that the post-hydrolysis state adopts a close-cleft conformation rather than an outward-open conformation. Our present simulations also highlight persistent interactions with annular cholesterol molecules and the expected active role of lipid bilayer in the allosteric communication between distant domains of the transporter.

Time-Resolved Mn2+ -NO and NO-NO Distance Measurements Reveal That Catalytic Asymmetry Regulates Alternating Access in an ABC Transporter

ATP-binding cassette (ABC) transporters shuttle diverse substrates across biological membranes. Transport is often achieved through a transition between an inward-facing (IF) and an outward-facing (OF) conformation of the transmembrane domains (TMDs). Asymmetric nucleotide-binding sites (NBSs) are present among several ABC subfamilies and their functional role remains elusive. Here we addressed this question using concomitant NO-NO, Mn2+ -NO, and Mn2+ -Mn2+ pulsed electron-electron double-resonance spectroscopy of TmrAB in a time-resolved manner. This type-IV ABC transporter undergoes a reversible transition in the presence of ATP with a significantly faster forward transition. The impaired degenerate NBS stably binds Mn2+ -ATP, and Mn2+ is preferentially released at the active consensus NBS. ATP hydrolysis at the consensus NBS considerably accelerates the reverse transition. Both NBSs fully open during each conformational cycle and the degenerate NBS may regulate the kinetics of this process.

The Alternating Access Mechanism in Mammalian Multidrug Resistance Transporters and Their Bacterial Homologs

Multidrug resistance (MDR) proteins belonging to the ATP-Binding Cassette (ABC) transporter group play a crucial role in the export of cytotoxic drugs across cell membranes. These proteins are particularly fascinating due to their ability to confer drug resistance, which subsequently leads to the failure of therapeutic interventions and hinders successful treatments. One key mechanism by which multidrug resistance (MDR) proteins carry out their transport function is through alternating access. This mechanism involves intricate conformational changes that enable the binding and transport of substrates across cellular membranes. In this extensive review, we provide an overview of ABC transporters, including their classifications and structural similarities. We focus specifically on well-known mammalian multidrug resistance proteins such as MRP1 and Pgp (MDR1), as well as bacterial counterparts such as Sav1866 and lipid flippase MsbA. By exploring the structural and functional features of these MDR proteins, we shed light on the roles of their nucleotide-binding domains (NBDs) and transmembrane domains (TMDs) in the transport process. Notably, while the structures of NBDs in prokaryotic ABC proteins, such as Sav1866, MsbA, and mammalian Pgp, are identical, MRP1 exhibits distinct characteristics in its NBDs. Our review also emphasizes the importance of two ATP molecules for the formation of an interface between the two binding sites of NBD domains across all these transporters. ATP hydrolysis occurs following substrate transport and is vital for recycling the transporters in subsequent cycles of substrate transportation. Specifically, among the studied transporters, only NBD2 in MRP1 possesses the ability to hydrolyze ATP, while both NBDs of Pgp, Sav1866, and MsbA are capable of carrying out this reaction. Furthermore, we highlight recent advancements in the study of MDR proteins and the alternating access mechanism. We discuss the experimental and computational approaches utilized to investigate the structure and dynamics of MDR proteins, providing valuable insights into their conformational changes and substrate transport. This review not only contributes to an enhanced understanding of multidrug resistance proteins but also holds immense potential for guiding future research and facilitating the development of effective strategies to overcome multidrug resistance, thus improving therapeutic interventions.

ATP binding and ATP hydrolysis in full-length MsbA monitored via time-resolved Fourier transform infrared spectroscopy

The essential Escherichia coli ATPase MsbA is a lipid flippase that serves as a prototype for multi drug resistant ABC transporters. Its physiological function is the transport of lipopolisaccharides to build up the outer membranes of gram negative bacteria. Although several structural and biochemical studies of MsbA have been conducted previously, a detailed picture of the dynamic processes that link ATP hydrolysis to allocrit transport remains elusive. We report here for the first time time-resolved Fourier transform infrared (FTIR) spectroscopic measurements of the ATP binding and ATP hydrolysis reaction of full-length MsbA and determined reaction rates at 288 K of k ₁ = 0.49 ± 0.28 s^-1 and k ₂ = 0.014 ± 0.003 s^-1, respectively. We further verified these rates with photocaged NPEcgAppNHp where only nucleotide binding was observable and the negative mutant MsbA-H537A that showed slow hydrolysis (k ₂ < 2 × 10^-4 s^-1). Besides single turnover kinetics, FTIR measurements also deliver IR signatures of all educts, products and the protein. ADP remains protein-bound after ATP hydrolysis. In addition, the spectral changes observed for the two variants MsbA-S378A and MsbA-S482A correlated with the loss of hydrogen bonding to the γ-phosphate of ATP. This study paves the way for FTIR-spectroscopic investigations of allocrite transport in full-length MsbA.

The effect of spin-lattice relaxation on DEER background decay

The common approach to background removal in double electron-electron resonance (DEER) measurements on frozen solutions with a three-dimensional homogeneous distribution of doubly labeled biomolecules is to fit the background to an exponential decay function. Excluded volume effects or distribution in a dimension lower than three, such as proteins in a membrane, can lead to a stretched exponential decay. In this work, we show that in cases of spin labels with short spin-lattice relaxation time, up to an order of magnitude longer than the DEER trace length, relevant for metal-based spin labels, spin flips that take place during the DEER evolution time affect the background decay shape. This was demonstrated using a series of temperature-dependent DEER measurements on frozen solutions of a nitroxide radical, a Gd(III) complex, Cu(II) ions, and a bis-Gd(III) model complex. As expected, the background decay was exponential for the nitroxide, whereas deviations were noted for Gd(III) and Cu(II). Based on the theoretical approach of Keller et al. (Phys. Chem. Chem. Phys. 21 (2019) 8228-8245), which addresses the effect of spin-lattice relaxation-induced spin flips during the evolution time, we show that the background decay can be fitted to an exponent including a linear and quadratic term in t, which is the position of the pump pulse. Analysis of the data in terms of the probability of spontaneous spin flips induced by spin-lattice relaxation showed that this approach worked well for the high temperature range studied for Gd(III) and Cu(II). At the low temperature range, the spin flips that occured during the DEER evolution time for Gd(III) exceeded the measured spin-lattice relaxation rate and include contributions from spin flips due to another mechanisms, most likely nuclear spin diffusion.

The ABC transporter MsbA adopts the wide inward-open conformation in E. coli cells

Membrane proteins are currently investigated after detergent extraction from native cellular membranes and reconstitution into artificial liposomes or nanodiscs, thereby removing them from their physiological environment. However, to truly understand the biophysical properties of membrane proteins in a physiological environment, they must be investigated within living cells. Here, we used a spin-labeled nanobody to interrogate the conformational cycle of the ABC transporter MsbA by double electron-electron resonance. Unexpectedly, the wide inward-open conformation of MsbA, commonly considered a nonphysiological state, was found to be prominently populated in Escherichia coli cells. Molecular dynamics simulations revealed that extensive lateral portal opening is essential to provide access of its large natural substrate core lipid A to the binding cavity. Our work paves the way to investigate the conformational landscape of membrane proteins in cells.

Dawn of a New Era for Membrane Protein Design

A major advancement has recently occurred in the ability to predict protein secondary structure from sequence using artificial neural networks. This new accessibility to high-quality predicted structures provides a big opportunity for the protein design community. It is particularly welcome for membrane protein design, where the scarcity of solved structures has been a major limitation of the field for decades. Here, we review the work done to date on the membrane protein design and set out established and emerging tools that can be used to most effectively exploit this new access to structures.

The Complex Roles of Adenosine Triphosphate in Bioenergetics

ATP is generally defined as the "energy currency" of the cell. Its phosphoanhydride P-O bonds are often considered to be "high energy" linkages that release free energy when broken, and its hydrolysis is described as "strongly exergonic". However, breaking bonds cannot release energy and ATP hydrolysis in motor and active transport proteins is not "strongly exergonic". So, the relevance of ATP resides elsewhere. As important as the nucleotide are the proteins that undergo functionally relevant conformational changes upon both ATP binding and release of ADP and inorganic phosphate. ATP phosphorylates proteins for signaling, active transport, and substrates in condensation reactions. The ensuing dephosphorylation has different consequences in each case. In signaling and active transport the phosphate group is hydrolyzed whereas in condensation reactions the phosphoryl fragment acts as a dehydrating agent. As it will be discussed in this article, ATP does much more than simply contribute free energy to biological processes.

Substrate-bound and substrate-free outward-facing structures of a multidrug ABC exporter

Multidrug ABC transporters translocate drugs across membranes by a mechanism for which the molecular features of drug release are so far unknown. Here, we resolved three ATP-Mg²⁺-bound outward-facing conformations of the Bacillus subtilis (homodimeric) BmrA by x-ray crystallography and single-particle cryo-electron microscopy (EM) in detergent solution, one of them with rhodamine 6G (R6G), a substrate exported by BmrA when overexpressed in B. subtilis. Two R6G molecules bind to the drug-binding cavity at the level of the outer leaflet, between transmembrane (TM) helices 1-2 of one monomer and TM5'-6' of the other. They induce a rearrangement of TM1-2, highlighting a local flexibility that we confirmed by hydrogen/deuterium exchange and molecular dynamics simulations. In the absence of R6G, simulations show a fast postrelease occlusion of the cavity driven by hydrophobicity, while when present, R6G can move within the cavity, maintaining it open.

Substrate-bound and substrate-free outward-facing structures of a multidrug ABC exporter

Multidrug ABC transporters translocate drugs across membranes by a mechanism for which the molecular features of drug release are so far unknown. Here, we resolved three ATP-Mg²⁺-bound outward-facing conformations of the Bacillus subtilis (homodimeric) BmrA by x-ray crystallography and single-particle cryo-electron microscopy (EM) in detergent solution, one of them with rhodamine 6G (R6G), a substrate exported by BmrA when overexpressed in B. subtilis. Two R6G molecules bind to the drug-binding cavity at the level of the outer leaflet, between transmembrane (TM) helices 1-2 of one monomer and TM5'-6' of the other. They induce a rearrangement of TM1-2, highlighting a local flexibility that we confirmed by hydrogen/deuterium exchange and molecular dynamics simulations. In the absence of R6G, simulations show a fast postrelease occlusion of the cavity driven by hydrophobicity, while when present, R6G can move within the cavity, maintaining it open.

Protein functional dynamics from the rigorous global analysis of DEER data: Conditions, components, and conformations

The potential of spin labeling to reveal the dynamic dimension of macromolecules has been recognized since the dawn of the methodology in the 1960s. However, it was the development of pulsed electron paramagnetic resonance spectroscopy to detect dipolar coupling between spin labels and the availability of turnkey instrumentation in the 21st century that realized the full promise of spin labeling. Double electron-electron resonance (DEER) spectroscopy has seen widespread applications to channels, transporters, and receptors. In these studies, distance distributions between pairs of spin labels obtained under different biochemical conditions report the conformational states of macromolecules, illuminating the key movements underlying biological function. These experimental studies have spurred the development of methods for the rigorous analysis of DEER spectroscopic data along with methods for integrating these distributions into structural models. In this tutorial, we describe a model-based approach to obtaining a minimum set of components of the distance distribution that correspond to functionally relevant protein conformations with a set of fractional amplitudes that define the equilibrium between these conformations. Importantly, we review and elaborate on the error analysis reflecting the uncertainty in the various parameters, a critical step in rigorous structural interpretation of the spectroscopic data.

Structural dynamics of ABC transporters: molecular simulation studies

The biological activities of living organisms involve various inputs and outputs. The ATP-driven substances (biomolecules) responsible for these kinds of activities through membrane (i.e. uptake and efflux of substrates) include ATP-binding cassette (ABC) transporters, some of which play important roles in multidrug resistance. The basic architecture of ABC transporters comprises transmembrane domains (TMDs) and nucleotide-binding domains (NBDs). The functional dynamics (substrate transport) of ABC transporters are realized by concerted motions, such as NBD dimerization, mechanical transmission via coupling helices (CHs), and the translocation of substrates through TMDs, which are induced by the binding and/or hydrolysis of ATP molecules and substrates. In this mini-review, we briefly discuss recent progresses in the structural dynamics as revealed by molecular simulation studies at all-atom (AA), coarse-grained (CG), and quantum mechanics/molecular mechanics (QM/MM) levels.

DEER-PREdict: Software for efficient calculation of spin-labeling EPR and NMR data from conformational ensembles

Owing to their plasticity, intrinsically disordered and multidomain proteins require descriptions based on multiple conformations, thus calling for techniques and analysis tools that are capable of dealing with conformational ensembles rather than a single protein structure. Here, we introduce DEER-PREdict, a software program to predict Double Electron-Electron Resonance distance distributions as well as Paramagnetic Relaxation Enhancement rates from ensembles of protein conformations. DEER-PREdict uses an established rotamer library approach to describe the paramagnetic probes which are bound covalently to the protein.DEER-PREdict has been designed to operate efficiently on large conformational ensembles, such as those generated by molecular dynamics simulation, to facilitate the validation or refinement of molecular models as well as the interpretation of experimental data. The performance and accuracy of the software is demonstrated with experimentally characterized protein systems: HIV-1 protease, T4 Lysozyme and Acyl-CoA-binding protein. DEER-PREdict is open source (GPLv3) and available at github.com/KULL-Centre/DEERpredict and as a Python PyPI package pypi.org/project/DEERPREdict.

Functional Oligomeric Forms of Uncoupling Protein 2: Strong Evidence for Asymmetry in Protein and Lipid Bilayer Systems

Stoichiometry of uncoupling proteins (UCPs) and their coexistence as functional monomeric and associated forms in lipid membranes remain intriguing open questions. In this study, tertiary and quaternary structures of UCP2 were analyzed experimentally and through molecular dynamics (MD) simulations. UCP2 was overexpressed in the inner membrane of Escherichia coli, then purified and reconstituted in lipid vesicles. Structure and proton transport function of UCP2 were characterized by circular dichroism (CD) spectroscopy and fluorescence methods. Findings suggest a tetrameric functional form for UCP2. MD simulations conclude that tetrameric UCP2 is a dimer of dimers, is more stable than its monomeric and dimeric forms, is asymmetrical and induces asymmetry in the membrane's lipid structure, and a biphasic on-off switch between the dimeric units is its possible mode of transport. MD simulations also show that the water density inside the UCP2 monomer is asymmetric, with the cytoplasmic side having a higher water density and a wider radius. In contrast, the structurally comparable adenosine 5'-diphosphate (ADP)/adenosine 5'-triphosphate (ATP) carrier (AAC1) did not form tetramers, implying that tetramerization cannot be generalized to all mitochondrial carriers.

From in vitro towards in situ: structure-based investigation of ABC exporters by electron paramagnetic resonance spectroscopy

ATP-binding cassette (ABC) exporters have been studied now for more than four decades, and recent structural investigation has produced a large number of protein database entries. Yet, important questions about how ABC exporters function at the molecular level remain debated, such as which are the molecular recognition hotspots and the allosteric couplings dynamically regulating the communication between the catalytic cycle and the export of substrates. This conundrum mainly arises from technical limitations confining all research to in vitro analysis of ABC transporters in detergent solutions or embedded in membrane-mimicking environments. Therefore, a largely unanswered question is how ABC exporters operate in situ, namely in the native membrane context of a metabolically active cell. This review focuses on novel mechanistic insights into type I ABC exporters gained through a unique combination of structure determination, biochemical characterization, generation of conformation-specific nanobodies/sybodies and double electron-electron resonance.

The role of the degenerate nucleotide binding site in type I ABC exporters

ATP‐binding cassette (ABC) transporters are fascinating molecular machines that are capable of transporting a large variety of chemically diverse compounds. The energy required for translocation is derived from binding and hydrolysis of ATP. All ABC transporters share a basic architecture and are composed of two transmembrane domains and two nucleotide binding domains (NBDs). The latter harbor all conserved sequence motifs that hallmark the ABC transporter superfamily. The NBDs form the nucleotide binding sites (NBSs) in their interface. Transporters with two active NBSs are called canonical transporters, while ABC exporters from eukaryotic organisms, including humans, frequently have a degenerate NBS1 containing noncanonical residues that strongly impair ATP hydrolysis. Here, we summarize current knowledge on degenerate ABC transporters. By integrating structural information with biophysical and biochemical evidence of asymmetric function, we develop a model for the transport cycle of degenerate ABC transporters. We will elaborate on the unclear functional advantages of a degenerate NBS.

Structures of ABC transporters: handle with care

In the past two decades, the ATP‐binding cassette (ABC) transporters' field has undergone a structural revolution. The importance of structural biology to the development of the field of ABC transporters cannot be overstated, as the ensemble of structures not only revealed the architecture of ABC transporters but also shaped our mechanistic view of these remarkable molecular machines. Nevertheless, we advocate that the mechanistic interpretation of the structures is not trivial and should be carried out with prudence. Herein, we bring several examples of structures of ABC transporters that merit re‐interpretation via careful comparison to experimental data. We propose that it is of the upmost importance to place new structures within the context of the available experimental data.

Thermodynamic Basis for Conformational Coupling in an ATP-Binding Cassette Exporter

ATP-binding cassette (ABC) transporters constitute one of the largest protein superfamilies, and they mediate the transport of diverse substrates across the membrane. The molecular mechanism for transducing the energy from ATP binding and hydrolysis into the conformational changes remains elusive. Here, we determined the thermodynamics underlying the ATP-induced global conformational switching for the ABC exporter TmrAB using temperature-resolved pulsed electron-electron double resonance (PELDOR or DEER) spectroscopy. We show that a strong entropy-enthalpy compensation mechanism enables the closure of the nucleotide-binding domains (NBDs) over a wide temperature range. This is mechanically coupled with an outward opening of the transmembrane domains (TMDs) accompanied by an entropy gain. The conserved catalytic glutamate plays a key role in the overall energetics. Our results reveal the thermodynamic basis for the chemomechanical energy coupling in an ABC exporter and present a new strategy to explore the energetics of similar membrane protein complexes.

Capturing Substrate Translocation in an ABC Exporter at the Atomic Level

ATP-binding cassette (ABC) transporters chemomechanically couple ATP binding and hydrolysis to large-scale conformational changes, ultimately leading to substrate translocation across biological membranes. Despite recent progress in the structure determination of substrate-bound ABC exporters, the inherently dynamic mechanism of substrate transport remains unclear at the atomic level. In this work, we capture substrate translocation in the heterodimeric ABC exporter TM287/288 from the hyperthermophilic bacterium Thermotoga maritima using all-atom molecular dynamics (MD) simulations. Unguided multimicrosecond simulations at 375 K show how the drugs daunorubicin and verapamil, which were initially docked into the ABC transporter, get translocated through the exporter by following its large-scale alternating access conformational transitions between an inward-facing (IF) and an outward-facing (OF) conformation. Triggered by the affinity difference due to differential solvation of the binding cavity in the IF and OF conformations, the substrates unbind from the OF transporter and partition into the lipid bilayer. While daunorubicin is stably inserted into the outer leaflet of the bilayer, verapamil dynamically flip flops between the bilayer leaflets, possibly rendering its net transport futile.

Computational Dissection of Membrane Transport at a Microscopic Level

Membrane transporters are key gatekeeper proteins at cellular membranes that closely control the traffic of materials. Their function relies on structural rearrangements of varying degrees that facilitate substrate translocation across the membrane. Characterizing these functionally important molecular events at a microscopic level is key to our understanding of membrane transport, yet challenging to achieve experimentally. Recent advances in simulation technology and computing power have rendered molecular dynamics (MD) simulation a powerful biophysical tool to investigate a wide range of dynamical events spanning multiple spatial and temporal scales. Here, we review recent studies of diverse membrane transporters using computational methods, with an emphasis on highlighting the technical challenges, key lessons learned, and new opportunities to illuminate transporter structure and function.

Spin-labeled nanobodies as protein conformational reporters for electron paramagnetic resonance in cellular membranes

Nanobodies are emerging tools in a variety of fields such as structural biology, cell imaging, and drug discovery. Here we pioneer the use of their spin-labeled variants as reporters of conformational dynamics of membrane proteins using DEER spectroscopy. At the example of the bacterial ABC transporter TM287/288, we show that two gadolinium-labeled nanobodies allow us to quantify, via analysis of the modulation depth of DEER traces, the fraction of transporters adopting the outward-facing state under different experimental conditions. Additionally, we quantitatively follow the interconversion from the outward- to the inward-facing state in the conformational ensemble under ATP turnover conditions. We finally show that the specificity of the nanobodies for the target protein allows the direct attainment of structural information on the wild-type TM287/288 expressed in cellular membranes without the need to purify or label the investigated membrane protein.

Chemo-Mechanical Coupling in the Transport Cycle of a Heme ABC Transporter

The heme importer from pathogenic bacteria is a member of the ATP-binding cassette (ABC) transporter family, which uses the energy of ATP-binding and hydrolysis for extensive conformational changes. Previous studies have indicated that conformational changes after heme translocation are triggered by ATP-binding to nucleotide binding domains (NBDs) and then, in turn, induce conformational transitions of the transmembrane domains (TMDs). In this study, we applied a template-based iterative all-atom molecular dynamics (MD) simulation to predict the ATP-bound outward-facing conformation of the Burkholderia cenocepacia heme importer BhuUV-T. The resulting model showed a stable conformation of the TMD with the cytoplasmic gate in the closed state and the periplasmic gate in the open state. Furthermore, targeted MD simulation predicted the intermediate structure of an occluded form (Occ) with bound ATP, in which both ends of the heme translocation channel are closed. The MD simulation of the predicted Occ revealed that Ser147 on the ABC signature motifs (LSGG[Q/E]) of NBDs occasionally flips and loses the active conformation required for ATP-hydrolysis. The flipping motion was found to be coupled to the inter-NBD distance. Our results highlight the functional significance of the signature motif of ABC transporters in regulation of ATPase and chemo-mechanical coupling mechanism.

The extracellular gate shapes the energy profile of an ABC exporter

ABC exporters harness the energy of ATP to pump substrates across membranes. Extracellular gate opening and closure are key steps of the transport cycle, but the underlying mechanism is poorly understood. Here, we generated a synthetic single domain antibody (sybody) that recognizes the heterodimeric ABC exporter TM287/288 exclusively in the presence of ATP, which was essential to solve a 3.2 Å crystal structure of the outward-facing transporter. The sybody binds to an extracellular wing and strongly inhibits ATPase activity by shifting the transporter's conformational equilibrium towards the outward-facing state, as shown by double electron-electron resonance (DEER). Mutations that facilitate extracellular gate opening result in a comparable equilibrium shift and strongly reduce ATPase activity and drug transport. Using the sybody as conformational probe, we demonstrate that efficient extracellular gate closure is required to dissociate the NBD dimer after ATP hydrolysis to reset the transporter back to its inward-facing state.

Lipid-Dependent Alternating Access Mechanism of a Bacterial Multidrug ABC Exporter

By undergoing conformational changes, active membrane transporters alternate between an inward-facing (IF) and an outward-facing (OF) state to transport their substrates across cellular membrane. The conformational landscape of membrane transporters, however, could be influenced by their environment, and the dependence of the alternating access mechanism on the lipid composition has not been understood at the molecular level. We have performed an extensive set of microsecond-level all-atom molecular dynamics (MD) simulations on bacterial ATP binding cassette (ABC) exporter Sav1866 in six different phosphocholine (PC) and phosphoethanolamine (PE) lipid membrane environments. This study mainly focuses on the energetically downhill OF-to-IF conformational transition of Sav1866 upon the ATP hydrolysis. We observe that the transporter undergoes large-scale conformational changes in the PE environment, particularly in the POPE lipids, resulting in an IF-occluded conformation, a transition that does not occur when the transporter is embedded in any of the PC lipid bilayers. We propose that the PE lipids facilitate the closing of the protein on the periplasmic side due to their highly polar headgroups that mediate the interaction of the two transmembrane (TM) bundles by a network of lipid-lipid and lipid-protein hydrogen bonds. POPE lipids in particular facilitate the closure of periplasmic gate by promoting a hinge formation in TM helices and an interbundle salt bridge formation. This study explains how the alternating access mechanism and the flippase activity in ABC exporters could be lipid-dependent.

EPR Techniques to Probe Insertion and Conformation of Spin-Labeled Proteins in Lipid Bilayers

Electron paramagnetic resonance (EPR) spectroscopy of spin-labeled membrane proteins is a valuable biophysical technique to study structural details and conformational transitions of proteins close to their physiological environment, for example, in liposomes, membrane bilayers, and nanodiscs. Unlike in nuclear magnetic resonance (NMR) spectroscopy, having only one or few specific side chains labeled at a time with paramagnetic probes makes the size of the object under investigation irrelevant in terms of technique sensitivity. As a drawback, extensive site-directed mutagenesis is required in order to analyze the properties of the protein under investigation. EPR can provide detailed information on side chain dynamics of large membrane proteins or protein complexes embedded in membranes with an exquisite sensitivity for flexible regions and on water accessibility profiles across the membrane bilayer. Moreover, distances between the two spin-labeled side chains in membrane proteins can be detected with high precision at cryogenic temperatures. The application of EPR to membrane proteins still presents some challenges in terms of sample preparation, sensitivity and data interpretation, thus it is difficult to give ready-to-go methodological recipes. However, new technological developments (arbitrary waveform generators) and new spin labels spectroscopically orthogonal to nitroxides increased the range of applicability from in vitro toward in-cell EPR experiments. This chapter is an updated version of the one published in the first edition of the book and describes the state of the art in the application of nitroxide-based site-directed spin labeling EPR to membrane proteins, addressing new tools such as arbitrary waveform generators and spectroscopically orthogonal labels, such as Gd(III)-based labels. We will present challenges in sample preparation and data analysis for functional and structural membrane protein studies using site-directed spin labeling techniques and give experimental details on EPR techniques providing information on side chain dynamics and water accessibility using nitroxide probes. An updated optimal Q-band DEER setup for nitroxide probes will be described, and its extension to gadolinium-containing samples will be addressed.

Molecular Mechanism of ATP Hydrolysis in an ABC Transporter

Hydrolysis of nucleoside triphosphate (NTP) plays a key role for the function of many biomolecular systems. However, the chemistry of the catalytic reaction in terms of an atomic-level understanding of the structural, dynamic, and free energy changes associated with it often remains unknown. Here, we report the molecular mechanism of adenosine triphosphate (ATP) hydrolysis in the ATP-binding cassette (ABC) transporter BtuCD-F. Free energy profiles obtained from hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations show that the hydrolysis reaction proceeds in a stepwise manner. First, nucleophilic attack of an activated lytic water molecule at the ATP γ-phosphate yields ADP + HPO₄ ^2- as intermediate product. A conserved glutamate that is located very close to the γ-phosphate transiently accepts a proton and thus acts as catalytic base. In the second step, the proton is transferred back from the catalytic base to the γ-phosphate, yielding ADP + H₂PO₄ ^-. These two chemical reaction steps are followed by rearrangements of the hydrogen bond network and the coordination of the Mg²⁺ ion. The rate constant estimated from the computed free energy barriers is in very good agreement with experiments. The overall free energy change of the reaction is close to zero, suggesting that phosphate bond cleavage itself does not provide a power stroke for conformational changes. Instead, ATP binding is essential for tight dimerization of the nucleotide-binding domains and the transition of the transmembrane domains from inward- to outward-facing, whereas ATP hydrolysis resets the conformational cycle. The mechanism is likely relevant for all ABC transporters and might have implications also for other NTPases, as many residues involved in nucleotide binding and hydrolysis are strictly conserved.

Combining theoretical and experimental data to decipher CFTR 3D structures and functions

Cryo-electron microscopy (cryo-EM) has recently provided invaluable experimental data about the full-length cystic fibrosis transmembrane conductance regulator (CFTR) 3D structure. However, this experimental information deals with inactive states of the channel, either in an apo, quiescent conformation, in which nucleotide-binding domains (NBDs) are widely separated or in an ATP-bound, yet closed conformation. Here, we show that 3D structure models of the open and closed forms of the channel, now further supported by metadynamics simulations and by comparison with the cryo-EM data, could be used to gain some insights into critical features of the conformational transition toward active CFTR forms. These critical elements lie within membrane-spanning domains but also within NBD1 and the N-terminal extension, in which conformational plasticity is predicted to occur to help the interaction with filamin, one of the CFTR cellular partners.