Analysis of the structural and functional roles of coupling helices in the ATP-binding cassette transporter MsbA through enzyme assays and molecular dynamics simulations

Probing the allosteric NBD-TMD crosstalk in the ABC transporter MsbA by solid-state NMR

The ABC transporter MsbA plays a critical role in Gram-negative bacteria in the regulation of the outer membrane by translocating core-LPS across the inner membrane. Additionally, a broad substrate specificity for lipophilic drugs has been shown. The allosteric interplay between substrate binding in the transmembrane domains and ATP binding and turnover in the nucleotide-binding domains must be mediated via the NBD/TMD interface. Previous studies suggested the involvement of two intracellular loops called coupling helix 1 and 2 (CH1, CH2). Here, we demonstrate by solid-state NMR spectroscopy that substantial chemical shift changes within both CH1 and CH2 occur upon substrate binding, in the ATP hydrolysis transition state, and upon inhibitor binding. CH2 is domain-swapped within the MsbA structure, and it is noteworthy that substrate binding induces a larger response in CH2 compared to CH1. Our data demonstrate that CH1 and CH2 undergo structural changes as part of the TMD-NBD cross-talk.

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.

MsbA: an ABC transporter paradigm

ATP-binding cassette (ABC) transporters play an important role in various cellular processes. They display a similar architecture and share a mechanism which couples ATP hydrolysis to substrate transport. However, in the light of current data and recent experimental progress, this protein superfamily appears as multifaceted as their broad substrate range. Among the prokaryotic ABC transporters, MsbA can serve as a paradigm for research in this field. It is located in the inner membrane of Gram-negative bacteria and functions as a floppase for the lipopolysaccharide (LPS) precursor core-LPS, which is involved in the biogenesis of the bacterial outer membrane. While MsbA shows high similarity to eukaryotic ABC transporters, its expression in Gram-negative bacteria makes it conveniently accessible for many experimental approaches from spectroscopy to 3D structure determination. As an essential protein for bacterial membrane integrity, MsbA has also become an attractive target for the development of novel antibiotics. Furthermore, it serves as a model for multidrug efflux pumps. Here we provide an overview of recent findings and their relevance to the field, highlight the potential of methods such as solid-state NMR and EPR spectroscopy and provide a perspective for future work.

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.

Development of an ATP force field for coarse-grained simulation of ATPases and its application to the maltose transporter

The biological functions of ATPases, such as myosin, kinesin, and ABC transporter, are due to large conformational motions driven by energy obtained from ATP. Elucidation of the mechanisms underlying these ATP-driven movements is one of the greatest challenges in computational chemistry. It has been shown that the MARTINI coarse-grained method is a promising tool for the investigation of large conformational motions in various proteins. However, this method has not yet been applied to ATPases because of the lack of a force field for the ATP molecule. Here, we developed force field parameters for the ATP molecule and conducted simulations using these parameters for the subunits (MalK₂ ) and the full-length structure (MalFGK₂ -E) of a maltose transporter. It was found for both targets that the dimerization of the nucleotide binding domains (NBDs) is induced upon ATP binding. Moreover, for the full-length transporter, the conformational transition from the pre-translocation state to the outward-facing state was observed and was accompanied by an initial transport motion of the substrate. It is expected that coarse-grained simulations utilizing the parameters for the ATP molecule developed here will serve as a powerful tool for investigating other ATPases as well. © 2019 Wiley Periodicals, Inc.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Relationship between conformation shift and disease related variation sites in ATP-binding cassette transporter proteins

Transport of small molecules across the cell membrane is a crucial biological mechanism for the maintenance of the cell activity. ABC transporter family is a huge group in the transporter membrane proteins and actively transports the substrates using the energy derived from ATP hydrolysis. In humans, there are 48 distinct genes for ABC transporters. A variation of a single amino acid in the amino acid sequence of ABC transporter has been known to be linked with certain disease. The mechanism of the onset of the disease by the variation is, however, still unclear. Recent progress in the method to measure the structures of huge membrane proteins has enabled determination of the 3D structures of ABC transporters and the accumulation of coordinate data of ABC transporter has enabled us to obtain clues for the onset of the disease caused by a single variation of amino acid residue. We compared the structures of ABC transporter in apo and ATP-binding forms and found a possible conformation shift around pivot-like residues in the transmembrane domains. When this conformation change in ABC transporter and the location of pathogenic variation were compared, we found a reasonable match between the two, explaining the onset of the disease by the variation. They likely cause impairment of the pivot-like movement, weakening of ATP binding and weakening of membrane surface interactions. These findings will give a new interpretation of the variations on ABC transporter genes and pave a way to analyse the effect of variation on protein structure and function.

Exploring movement and energy in human P-glycoprotein conformational rearrangement

Human P-glycoprotein (P-gp), a kind of ATP-Binding Cassette transporter, can export a diverse variety of anti-cancer drugs out of the tumor cell. Its overexpression is one of the main reasons for the multidrug resistance (MDR) of tumor cells. It has been confirmed that during the substrate transport process, P-gp experiences a large-scale structural rearrangement from the inward- to outward-facing states. However, the mechanism of how the nucleotide-binding domains (NBDs) control the transmembrane domains (TMDs) to open towards the periplasm in the outward-facing state has not yet been fully characterized. Herein, targeted molecular dynamics simulations were performed to explore the conformational rearrangement of human P-gp. The results show that the allosteric process proceeds in a coupled way, and first the transition is driven by the NBDs, and then transmitted to the cytoplasmic parts of TMDs, finally to the periplasmic parts. The trajectories show that besides the translational motions, the NBDs undergo a rotation movement, which mainly occurs in xy plane and ensures the formation of the correct ATP-binding pockets. The analyses on the interaction energies between the six structure segments (cICLs) from the TMDs and NBDs reveal that their subtle energy differences play an important role in causing the periplasmic parts of the transmembrane helices to separate from each other in the established directions and in appropriate amplitudes. This conclusion can explain the two experimental phenomena about human P-gp in some extent. These studies have provided a detailed exploration into human P-gp rearrangement process and given an energy insight into the TMD reorientation during P-gp transition.

Molecular dynamics simulation study on the structural instability of the most common cystic fibrosis-associated mutant ΔF508-CFTR

Cystic fibrosis transmembrane conductance regulator (CFTR) is an anion channel that belongs to the ATP binding cassette protein superfamily. Deletion of phenylalanine at position 508 (ΔF508) is the most common CF-associated mutation and is present in nearly 90% of CF patients. Currently, atomistic level studies are insufficient for understanding the mechanism by which the deletion of a single amino acid causes greatly reduced folding as well as trafficking and gating defects. To clarify this mechanism, we first constructed an atomic model of the inward-facing ΔF508-CFTR and performed allatom molecular dynamics (MD) simulations of the protein in a membrane environment. All of the computational methodologies used are based on those developed in our previous study for wild-type CFTR. Two important findings were obtained. First, consistent with several previous computational results, the deletion of F508 causes a disruption of a hydrophobic cluster located at the interface between the nucleotide binding domain 1 (NBD1) and intracellular loop 4 (ICL4). This exerts unfavorable influences on the correlated motion between ICLs and transmembrane domains (TMDs), likely resulting in gating defects. Second, the F508 deletion affected the NBD1-NBD2 interface via allosteric communication originating from the correlated motion between NBDs and ICLs. As a result, several unusual inter-residue interactions are caused at the NBD1-NBD2 interface. In other words, their correct dimerization is impaired. This study provided insight into the atomic-level details of structural and dynamics changes caused by the ΔF508 mutation and thus provides good insight for drug design.

Analysis of an ATP-induced conformational transition of ABC transporter MsbA using a coarse-grained model

Upon the binding of ATP molecules to nucleotide binding domains (NBDs), ATP-binding cassette (ABC) exporters undergo a conformational transition from an inward-facing (IF) to an outward-facing (OF) state. This molecular event is a typical example of chemo-mechanical coupling. However, the underlying mechanism remains unclear. In this study, we analyzed the IF→OF transition of a representative ABC exporter, MsbA, by solving the equation of motion under an elastic network model (ENM). ATP was represented as a single node in ENM or replaced by external forces. When two ATP nodes were added to the ENM of the IF state protein, the two NBDs dimerized; subsequently, the two transmembrane domains opened toward the extracellular side, resulting in the formation of the OF structure. Such a conformational transition was also reproduced by applying external forces, which caused the rotational motion of the NBDs instead of the addition of ATP nodes. The process of the conformational transition was analyzed in detail using cross-correlation maps for node-node interactions. More importantly, it was revealed that the ATP binding energy is converted into distortion energy of several transmembrane helices. These results are useful for understanding the chemo-mechanical coupling in ABC transporters.

Probing native metal ion association sites through quenching of fluorophores in the nucleotide-binding domains of the ABC transporter MsbA

ATP-binding cassette (ABC) transporters are ubiquitously present in prokaryotic and eukaryotic cells. Binding of ATP to the nucleotide-binding domains (NBDs) elicits major conformational changes of the transporters resulting in the transport of the substrate across the membrane. The availability of a crystal structure of the NBDs enabled us to elucidate the local structure and small-scale dynamics in the NBDs. Here, we labeled the ABC transporter MsbA, a homodimeric flippase from Escherichia coli, with a fluorescent probe, Alexa532, within the NBDs. ATP application elicited collisional quenching, whereas no quenching was observed after the addition of ATP analogs or ATP hydrolysis inhibitors. The Alexa532-conjugated MsbA variants exhibited transition metal ion Förster resonance energy transfer (tmFRET) after the addition of Ni²⁺, and ATP decreased this Ni²⁺-mediated FRET of the NBDs. Structure modeling developed from crystallographic data and examination of tmFRET measurements of MsbA variants in the absence of ATP revealed the presence of metal ion-associated pockets (MiAPs) in the NBDs. Three histidines were predicted to participate in chelating Ni²⁺ in the two possible MiAPs. Performing histidine-substitution experiments with the NBDs showed that the dissociation constant for Ni²⁺ of MiAP2 was smaller than that of MiAP1. The structural allocation of the MiAPs was further supported by showing that the addition of Cu²⁺ resulted in higher quenching than Ni²⁺ Taken together, the present study showed that the NBDs contain two native binding sites for metal ions and ATP addition affects the Ni²⁺-binding activity of the MiAPs.

Structural Dynamics of the Heterodimeric ABC Transporter TM287/288 Induced by ATP and Substrate Binding

TM287/288 is a heterodimeric ATP-binding cassette (ABC) transporter, which harnesses the energy of ATP binding and hydrolysis at the nucleotide-binding domains (NBDs) to transport a wide variety of molecules through the transmembrane domains (TMDs) by alternating inward- and outward-facing conformations. Here, we conducted multiple 100 ns molecular dynamics simulations of TM287/288 in different ATP- and substrate-bound states to elucidate the effects of ATP and substrate binding. As a result, the binding of two ATP molecules to the NBDs induced the formation of the consensus ATP-binding pocket (ABP2) or the NBD dimerization, whereas these processes did not occur in the presence of a single ATP molecule or when the protein was in its apo state. Moreover, binding of the substrate to the TMDs enhanced the formation of ABP2 through allosteric TMD-NBD communication. Furthermore, in the apo state, α-helical subdomains of the NBDs approached each other, acquiring a conformation with core half-pockets exposed to the solvent, appropriate for ATP binding. We propose a "core-exposed" model for this novel conformation found in the apo state of ABC transporters. These findings provide important insights into the structural dynamics of ABC transporters.

Investigating the dynamic nature of the ABC transporters: ABCB1 and MsbA as examples for the potential synergies of MD theory and EPR applications

ABC transporters are primary active transporters found in all kingdoms of life. Human multidrug resistance transporter ABCB1, or P-glycoprotein, has an extremely broad substrate spectrum and confers resistance against chemotherapy drug treatment in cancer cells. The bacterial ABC transporter MsbA is a lipid A flippase and a homolog to the human ABCB1 transporter, with which it partially shares its substrate spectrum. Crystal structures of MsbA and ABCB1 have been solved in multiple conformations, providing a glimpse into the possible conformational changes the transporter could be going through during the transport cycle. Crystal structures are inherently static, while a dynamic picture of the transporter in motion is needed for a complete understanding of transporter function. Molecular dynamics (MD) simulations and electron paramagnetic resonance (EPR) spectroscopy can provide structural information on ABC transporters, but the strength of these two methods lies in the potential to characterise the dynamic regime of these transporters. Information from the two methods is quite complementary. MD simulations provide an all atom dynamic picture of the time evolution of the molecular system, though with a narrow time window. EPR spectroscopy can probe structural, environmental and dynamic properties of the transporter in several time regimes, but only through the attachment sites of an exogenous spin label. In this review the synergistic effects that can be achieved by combining the two methods are highlighted, and a brief methodological background is also presented.

Analysis of the Free Energy Landscapes for the Opening-Closing Dynamics of the Maltose Transporter ATPase MalK2 Using Enhanced-Sampling Molecular Dynamics Simulation

Protein dynamics are considered significant for many physiological processes, such as metabolism, biomolecular recognition, and the regulation of several vital cellular processes. Due to their flexibility, proteins may stay in different substates with or without the existence of the cognate substrates. To describe these phenomena, two models have been proposed: the "induced fit" and the "conformational selection" mechanisms. In this study, we used MalK2, the subunits that mainly include the nucleotide-binding domains (NBDs) of the maltose transporter from Escherichia coli, as a target to understand the NBD dimerization mechanism. Accelerated and conventional molecular dynamics have been performed. The results revealed that Mg-ATP binding to MalK2 led to a significant change in the free energy profile and thus stabilized the closed conformation. On the contrary, when Mg-ATP was removed, the open conformation would be favored. The fact that ligand binding induces a drastic free energy change leads to a significant inference: MalK2 dimerization would occur through the induced-fit mechanism rather than the conformational selection mechanism. This study sheds new light on the NBD dimerization mechanism and would be of wide applicability to other ABC transporters.