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

Cholesterol metabolism: A strategy for overcoming drug resistance in tumors

Despite significant advancements in targeted tumor therapies, the emergence of drug resistance remains a complex challenge. Cholesterol accumulation within tumor cells plays a crucial role in mediating drug resistance through various mechanisms, including altered membrane dynamics, enhanced drug efflux, and activation of survival signaling pathways. Targeting cholesterol metabolism presents an innovative strategy to enhance therapeutic sensitivity, particularly in breast cancer. Consequently, ongoing preclinical studies and clinical trials involving cholesterol-lowering agents indicate a promising direction for improving treatment outcomes in tumors. The combination of these agents with existing therapeutic regimens may lead to enhanced efficacy, highlighting the necessity for continued research in this vital area. This review examines the impact of cholesterol metabolism on drug resistance in tumors, particularly solid tumors, identifies therapeutic targets in this metabolic pathway (with a special focus on breast cancer), and discusses recent advances in cholesterol-lowering drugs in preclinical, as well as those that have entered clinical trials.

Dissecting Large-Scale Structural Transitions in Membrane Transporters Using Advanced Simulation Technologies

Membrane transporters are integral membrane proteins that act as gatekeepers of the cell, controlling fundamental processes such as recruitment of nutrients and expulsion of waste material. At a basic level, transporters operate using the "alternating access model," in which transported substances are accessible from only one side of the membrane at a time. This model usually involves large-scale structural changes in the transporter, which often cannot be captured using unbiased, conventional molecular simulation techniques. In this article, we provide an overview of some of the major simulation techniques that have been applied to characterize the structural dynamics and energetics involved in the transition of membrane transporters between their functional states. After briefly introducing each technique, we discuss some of their advantages and limitations and provide some recent examples of their application to membrane transporters.

ATP-binding cassette (ABC) transporters: structures and roles in bacterial pathogenesis

Adenosine triphosphate (ATP)-binding cassette (ABC) transporter systems are divided into importers and exporters that facilitate the movement of diverse substrate molecules across the lipid bilayer, against the concentration gradient. These transporters comprise two highly conserved nucleotide-binding domains (NBDs) and two transmembrane domains (TMDs). Unlike ABC exporters, prokaryotic ABC importers require an additional substrate-binding protein (SBP) as a recognition site for specific substrate translocation. The discovery of a large number of ABC systems in bacterial pathogens revealed that these transporters are crucial for the establishment of bacterial infections. The existing literature has highlighted the roles of ABC transporters in bacterial growth, pathogenesis, and virulence. These roles include importing essential nutrients required for a variety of cellular processes and exporting outer membrane-associated virulence factors and antimicrobial substances. This review outlines the general structures and classification of ABC systems to provide a comprehensive view of the activities and roles of ABC transporters associated with bacterial virulence and pathogenesis during infection.

Deciphering the Interdomain Coupling in a Gram-Negative Bacterial Membrane Insertase

YidC is a membrane protein that plays an important role in inserting newly generated proteins into lipid membranes. The Sec-dependent complex is responsible for inserting proteins into the lipid bilayer in bacteria. YidC facilitates the insertion and folding of membrane proteins, both in conjunction with the Sec complex and independently. Additionally, YidC acts as a chaperone during the folding of proteins. Multiple investigations have conclusively shown that Gram-positive bacterial YidC has Sec-independent insertion mechanisms. Through the use of microsecond-level all-atom molecular dynamics (MD) simulations, we have carried out an in-depth investigation of the YidC protein originating from Gram-negative bacteria. This research sheds light on the significance of multiple domains of the YidC structure at a detailed molecular level by utilizing equilibrium MD simulations. Specifically, multiple models of YidC embedded in the lipid bilayer were constructed to characterize the critical role of the C2 loop and the periplasmic domain (PD) present in Gram-negative YidC, which is absent in its Gram-positive counterpart. Based on our results, the C2 loop plays a role in the overall stabilization of the protein, most notably in the transmembrane (TM) region, and it also has an allosteric influence on the PD region. We have found critical inter- and intradomain interactions that contribute to the stability of the protein and its function. Finally, our study provides a hypothetical Sec-independent insertion mechanism for Gram-negative bacterial YidC.

Lipid-dependent conformational dynamics of bacterial ATP-binding cassette transporter Sav1866

Sav1866, a bacterial ATP-binding cassette (ABC) exporter, plays a crucial role in cellular processes by facilitating the efflux of a diverse range of substrates, including drugs, chemotherapeutic agents, peptides, and lipids. This efflux activity significantly impacts the effectiveness of various therapies against bacterial infections. In our recent investigation, we focused on understanding the conformational dynamics of Sav1866 within different lipid environments. Specifically, we explored its behavior in environments composed of DMPC and POPE lipids, which exhibit crucial distinctions not only in their headgroup polarity but also in the length and saturation of their hydrophobic tails. Our extensive set of equilibrium microsecond-level all-atom molecular dynamics (MD) simulations revealed significant distinctions in transporter behavior influenced by these lipid compositions. We observed a rapid transition to an occluded-inward-facing (IF-occ) conformation in POPE environments, contrasting with a channel-like behavior in DMPC environments, deviating from the expected alternating access mechanism (AAM). These findings underscore the significant impact of lipid compositions on ABC transporter function, offering new perspectives on membrane transport mechanisms.

Dynamic framework for large-scale modeling of membranes and peripheral proteins

In this chapter, we present a novel computational framework to study the dynamic behavior of extensive membrane systems, potentially in interaction with peripheral proteins, as an alternative to conventional simulation methods. The framework effectively describes the complex dynamics in protein-membrane systems in a mesoscopic particle-based setup. Furthermore, leveraging the hydrodynamic coupling between the membrane and its surrounding solvent, the coarse-grained model grounds its dynamics in macroscopic kinetic properties such as viscosity and diffusion coefficients, marrying the advantages of continuum- and particle-based approaches. We introduce the theoretical background and the parameter-space optimization method in a step-by-step fashion, present the hydrodynamic coupling method in detail, and demonstrate the application of the model at each stage through illuminating examples. We believe this modeling framework to hold great potential for simulating membrane and protein systems at biological spatiotemporal scales, and offer substantial flexibility for further development and parametrization.

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.

The Important Role of Transporter Structures in Drug Disposition, Efficacy, and Toxicity

The ATP-binding cassette (ABC) and solute carrier (SLC) transporters are critical determinants of drug disposition, clinical efficacy, and toxicity as they specifically mediate the influx and efflux of various substrates and drugs. ABC transporters can modulate the pharmacokinetics of many drugs via mediating the translocation of drugs across biologic membranes. SLC transporters are important drug targets involved in the uptake of a broad range of compounds across the membrane. However, high-resolution experimental structures have been reported for a very limited number of transporters, which limits the study of their physiologic functions. In this review, we collected structural information on ABC and SLC transporters and described the application of computational methods in structure prediction. Taking P-glycoprotein (ABCB1) and serotonin transporter (SLC6A4) as examples, we assessed the pivotal role of structure in transport mechanisms, details of ligand-receptor interactions, drug selectivity, the molecular mechanisms of drug-drug interactions, and differences caused by genetic polymorphisms. The data collected contributes toward safer and more effective pharmacological treatments. SIGNIFICANCE STATEMENT: The experimental structure of ATP-binding cassette and solute carrier transporters was collected, and the application of computational methods in structure prediction was described. P-glycoprotein and serotonin transporter were used as examples to reveal the pivotal role of structure in transport mechanisms, drug selectivity, the molecular mechanisms of drug-drug interactions, and differences caused by genetic polymorphisms.

cpSRP43 Is Both Highly Flexible and Stable: Structural Insights Using a Combined Experimental and Computational Approach

The novel multidomain protein, cpSRP43, is a unique subunit of the post-translational chloroplast signal recognition particle (cpSRP) targeting pathway in higher plants. The cpSRP pathway is responsible for targeting and insertion of light-harvesting chlorophyll a/b binding proteins (LHCPs) to the thylakoid membrane. Upon emergence into the stroma, LHCPs form a soluble transit complex with the cpSRP heterodimer, which is composed of cpSRP43 and cpSRP54. cpSRP43 is irreplaceable as a chaperone to LHCPs in their translocation to the thylakoid membrane and remarkable in its ability to dissolve aggregates of LHCPs without the need for external energy input. In previous studies, cpSRP43 has demonstrated significant flexibility and interdomain dynamics. In this study, we explore the structural stability and flexibility of cpSRP43 using a combination of computational and experimental techniques and find that this protein is concurrently highly stable and flexible. In addition to microsecond-level unbiased molecular dynamics (MD), biased MD simulations based on system-specific collective variables are used along with biophysical experimentation to explain the basis of the flexibility and stability of cpSRP43, showing that the free and cpSRP54-bound cpSRP43 has substantially different conformations and conformational dynamics.

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.

On the interplay between lipids and asymmetric dynamics of an NBS degenerate ABC transporter

Multidrug resistance-associated proteins are ABC C-family exporters. They are crucial in pharmacology as they transport various substrates across membranes. However, the role of the degenerate nucleotide-binding site (NBS) remains unclear likewise the interplay with the surrounding lipid environment. Here, we propose a dynamic and structural overview of MRP1 from ca. 110 μs molecular dynamics simulations. ATP binding to NBS1 is likely maintained along several transport cycles. Asymmetric NBD behaviour is ensured by lower signal transduction from NBD1 to the rest of the protein owing to the absence of ball-and-socket conformation between NBD1 and coupling helices. Even though surrounding lipids play an active role in the allosteric communication between the substrate-binding pocket and NBDs, our results suggest that lipid composition has a limited impact, mostly by affecting transport kinetics. We believe that our work can be extended to other degenerate NBS ABC proteins and provide hints for deciphering mechanistic differences among ABC transporters.

Differential dynamics and direct interaction of bound ligands with lipids in multidrug transporter ABCG2

ABCG2 is an ATP-binding cassette (ABC) transporter that extrudes a wide range of xenobiotics and drugs from the cell and contributes to multidrug resistance in cancer cells. Following our recent structural characterization of topotecan-bound ABCG2, here, we present cryo-EM structures of ABCG2 under turnover conditions in complex with a special modulator and slow substrate, tariquidar, in nanodiscs. The structures reveal that similar to topotecan, tariquidar induces two distinct ABCG2 conformations under turnover conditions (turnover-1 and turnover-2). μs-scale molecular dynamics simulations of drug-bound and apo ABCG2 in native-like lipid bilayers, in both topotecan- and tariquidar-bound states, characterize the ligand size as a major determinant of its binding stability. The simulations highlight direct lipid-drug interactions for the smaller topotecan, which exhibits a highly dynamic binding mode. In contrast, the larger tariquidar occupies most of the available volume in the binding pocket, thus leaving little space for lipids to enter the cavity and interact with it. Similarly, when simulating ABCG2 in the apo inward-open state, we also observe spontaneous penetration of phospholipids into the binding cavity. The captured phospholipid diffusion pathway into ABCG2 offers a putative general path to recruit any hydrophobic/amphiphilic substrates directly from the membrane. Our simulations also reveal that ABCG2 rejects cholesterol as a substrate, which is omnipresent in plasma membranes that contain ABCG2. At the same time, cholesterol is found to prohibit the penetration of phospholipids into ABCG2. These molecular findings have direct functional ramifications on ABCG2's function as a transporter.

Computing Individual Area per Head Group Reveals Lipid Bilayer Dynamics

Lipid bilayers express a range of phases from solid-like to gel-like to liquid-like as a function of temperature and lipid surface concentration. The area occupied per lipid head group serves as one useful indicator of the bilayer phase, in conjunction with the two-dimensional radial distribution function (i.e., structure factor) within the bilayer. Typically, the area per head group is determined by dividing the bilayer area equally among all head groups. Such an approach is less satisfactory for a multicomponent set of diverse lipids. In this work, area determination is performed on a lipid-by-lipid basis by attributing to a lipid the volume that surrounds each atom. Voronoi tessellation provides this division of the interfacial region on a per-atom basis. The method is applied to a multicomponent system of water, NaCl, and 19 phospholipid types that was devised recently [Langmuir2022, 38, 9481-9499] as a computational representation of the Gram-positive Staphylococcus aureus phospholipid bilayer. Results demonstrate that lipids and water molecules occupy similar extents of area within the interfacial region; ascribing area only to head groups implicitly incorporates assumptions about head group hydration. Results further show that lipid tails provide non-negligible contributions to area on the membrane side of the bilayer-water interface. Results for minimum and maximum area of individual lipids reveal that spontaneous fluctuations displace head groups more than 10 Å from the interfacial region during an NPT simulation at 310 K, leading to a zero contribution to total area at some times. Total area fluctuations and fluctuations per individual lipid relax with a correlation time of ∼10 ns. The method complements density profile as an approach to quantify the structure and dynamics of computational lipid bilayers.

An investigation of the YidC-mediated membrane insertion of Pf3 coat protein using molecular dynamics simulations

YidC is a membrane protein that facilitates the insertion of newly synthesized proteins into lipid membranes. Through YidC, proteins are inserted into the lipid bilayer via the SecYEG-dependent complex. Additionally, YidC functions as a chaperone in protein folding processes. Several studies have provided evidence of its independent insertion mechanism. However, the mechanistic details of the YidC SecY-independent protein insertion mechanism remain elusive at the molecular level. This study elucidates the insertion mechanism of YidC at an atomic level through a combination of equilibrium and non-equilibrium molecular dynamics (MD) simulations. Different docking models of YidC-Pf3 in the lipid bilayer were built in this study to better understand the insertion mechanism. To conduct a complete investigation of the conformational difference between the two docking models developed, we used classical molecular dynamics simulations supplemented with a non-equilibrium technique. Our findings indicate that the YidC transmembrane (TM) groove is essential for this high-affinity interaction and that the hydrophilic nature of the YidC groove plays an important role in protein transport across the cytoplasmic membrane bilayer to the periplasmic side. At different stages of the insertion process, conformational changes in YidC's TM domain and membrane core have a mechanistic effect on the Pf3 coat protein. Furthermore, during the insertion phase, the hydration and dehydration of the YidC's hydrophilic groove are critical. These results demonstrate that Pf3 coat protein interactions with the membrane and YidC vary in different conformational states during the insertion process. Finally, this extensive study directly confirms that YidC functions as an independent insertase.

Principles of Alternating Access in LeuT-fold Transporters: Commonalities and Divergences

Found in all domains of life, transporters belonging to the LeuT-fold class mediate the import and exchange of hydrophilic and charged compounds such as amino acids, metals, and sugar molecules. Nearly two decades of investigations on the eponymous bacterial transporter LeuT have yielded a library of high-resolution snapshots of its conformational cycle linked by solution-state experimental data obtained from multiple techniques. In parallel, its topology has been observed in symporters and antiporters characterized by a spectrum of substrate specificities and coupled to gradients of distinct ions. Here we review and compare mechanistic models of transport for LeuT, its well-studied homologs, as well as functionally distant members of the fold, emphasizing the commonalities and divergences in alternating access and the corresponding energy landscapes. Our integrated summary illustrates how fold conservation, a hallmark of the LeuT fold, coincides with divergent choreographies of alternating access that nevertheless capitalize on recurrent structural motifs. In addition, it highlights the knowledge gap that hinders the leveraging of the current body of research into detailed mechanisms of transport for this important class of membrane proteins.

Elucidating the molecular basis of spontaneous activation in an engineered mechanosensitive channel

Mechanosensitive channel of large conductance (MscL) detects and responds to changes in the pressure profile of cellular membranes and transduces the mechanical energy into electrical and/or chemical signals. MscL can be activated using ultrasonic or chemical activation methods to improve the absorption of medicines and bioactive compounds into cells. However, re-engineering chemical signals such as pH change can trigger channel activation in MscL. This study elucidates the activation mechanism of an engineered MscL at an atomic level through a combination of equilibrium and non-equilibrium (NE) molecular dynamics (MD) simulations. Comparing the wild-type (WT) and engineered MscL activation processes suggests that the two systems are likely associated with different active states and different transition pathways. These findings indicate that (1) periplasmic loops play a key role in the activation process of MscL, (2) the loss of various backbone-backbone hydrogen bonds and salt bridge interactions in the engineered MscL channel causes the spontaneous opening of the channel, and (3) the most significant interactions lost during the activation process are between the transmembrane helices 1 and 2 in engineered MscL channel. The orientation-based biasing approach for producing and optimizing an open MscL model used in this work is a promising way to characterize unknown protein functional states and investigate the activation processes in ion channels and transmembrane proteins in general. This work paves the way for a computational framework for engineering more efficient pH-sensing mechanosensitive channels.

Sampling alternative conformational states of transporters and receptors with AlphaFold2

Equilibrium fluctuations and triggered conformational changes often underlie the functional cycles of membrane proteins. For example, transporters mediate the passage of molecules across cell membranes by alternating between inward- and outward-facing states, while receptors undergo intracellular structural rearrangements that initiate signaling cascades. Although the conformational plasticity of these proteins has historically posed a challenge for traditional de novo protein structure prediction pipelines, the recent success of AlphaFold2 (AF2) in CASP14 culminated in the modeling of a transporter in multiple conformations to high accuracy. Given that AF2 was designed to predict static structures of proteins, it remains unclear if this result represents an underexplored capability to accurately predict multiple conformations and/or structural heterogeneity. Here, we present an approach to drive AF2 to sample alternative conformations of topologically diverse transporters and G-protein-coupled receptors that are absent from the AF2 training set. Whereas models of most proteins generated using the default AF2 pipeline are conformationally homogeneous and nearly identical to one another, reducing the depth of the input multiple sequence alignments by stochastic subsampling led to the generation of accurate models in multiple conformations. In our benchmark, these conformations spanned the range between two experimental structures of interest, with models at the extremes of these conformational distributions observed to be among the most accurate (average template modeling score of 0.94). These results suggest a straightforward approach to identifying native-like alternative states, while also highlighting the need for the next generation of deep learning algorithms to be designed to predict ensembles of biophysically relevant states.

MLCV: Bridging Machine-Learning-Based Dimensionality Reduction and Free-Energy Calculation

Importance-sampling algorithms leaning on the definition of a model reaction coordinate (RC) are widely employed to probe processes relevant to chemistry and biology alike, spanning time scales not amenable to common, brute-force molecular dynamics (MD) simulations. In practice, the model RC often consists of a handful of collective variables (CVs) chosen on the basis of chemical intuition. However, constructing manually a low-dimensional RC model to describe an intricate geometrical transformation for the purpose of free-energy calculations and analyses remains a daunting challenge due to the inherent complexity of the conformational transitions at play. To solve this issue, remarkable progress has been made in employing machine-learning techniques, such as autoencoders, to extract the low-dimensional RC model from a large set of CVs. Implementation of the differentiable, nonlinear machine-learned CVs in common MD engines to perform free-energy calculations is, however, particularly cumbersome. To address this issue, we present here a user-friendly tool (called MLCV) that facilitates the use of machine-learned CVs in importance-sampling simulations through the popular Colvars module. Our approach is critically probed with three case examples consisting of small peptides, showcasing that through hard-coded neural network in Colvars, deep-learning and enhanced-sampling can be effectively bridged with MD simulations. The MLCV code is versatile, applicable to all the CVs available in Colvars, and can be connected to any kind of dense neural networks. We believe that MLCV provides an effective, powerful, and user-friendly platform accessible to experts and nonexperts alike for machine-learning (ML)-guided CV discovery and enhanced-sampling simulations to unveil the molecular mechanisms underlying complex biochemical processes.

Toward Convergence in Free Energy Calculations for Protein Conformational Changes: A Case Study on the Thin Gate of Mhp1 Transporter

It has been challenging to obtain reliable free energies for protein conformational changes from all-atom molecular dynamics simulations, despite the availability of many enhanced sampling techniques. To alleviate the difficulties associated with the enormous complexity of the conformational space, here we propose a few practical strategies for such calculations, including (1) a stringent method to examine convergence by comparing independent simulations starting from different initial coordinates, (2) adoption of multistep schemes in which the complete conformational change consists of multiple transition steps, each sampled using a distinct reaction coordinate, and (3) application of boundary restraints to simplify the conformational space. We demonstrate these strategies on the conformational changes between the outward-facing and outward-occluded states of the Mhp1 membrane transporter, obtaining the equilibrium thermodynamics of the relevant metastable states, the kinetic rates between these states, and the reactive trajectories that reveal the atomic details of spontaneous transitions. Our approaches thus promise convergent and reliable calculations to examine intuition-based hypotheses and to eventually elucidate the underlying molecular mechanisms of reversible conformational changes in complex protein systems.

Mechanistic Picture for Chemomechanical Coupling in a Bacterial Proton-Coupled Oligopeptide Transporter from Streptococcus Thermophilus

Proton-coupled oligopeptide transporters (POTs) use the proton electrochemical gradient to transport peptides across the cell membrane. Despite the significant biological and biomedical relevance of these proteins, a detailed mechanistic picture for chemomechanical couplings involved in substrate/proton transport and protein structural changes is missing. Therefore, we performed microsecond-level molecular dynamics simulations of bacterial POT PepT_St, which shares ∼80% sequence identity with the human POT, PepT1, in the substrate-binding region. Three different conformational states of PepT_St were simulated, including (i) occluded, apo, (ii) inward-facing, apo, and (iii) inward-facing_occluded, Leu-Ala bound. We propose that the interaction of R33 with E299 and E300 acts as a conformational switch (i.e., to trigger the conformational change from an inward- to outward-facing state) in the substrate transport. Additionally, we propose that E299 and E400 disengage from interacting with the substrate either through protonation or through coordination with a cation for the substrate to get transported. This study provides clues to understand the chemomechanical couplings in POTs and paves the way to decipher the molecular-level underpinnings of the structure-function relationship in this important family of transporters.

Selective Nutrient Transport in Bacteria: Multicomponent Transporter Systems Reign Supreme

Multicomponent transporters are used by bacteria to transport a wide range of nutrients. These systems use a substrate-binding protein to bind the nutrient with high affinity and then deliver it to a membrane-bound transporter for uptake. Nutrient uptake pathways are linked to the colonisation potential and pathogenicity of bacteria in humans and may be candidates for antimicrobial targeting. Here we review current research into bacterial multicomponent transport systems, with an emphasis on the interaction at the membrane, as well as new perspectives on the role of lipids and higher oligomers in these complex systems.

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.

Identification and characterisation of putative drug binding sites in human ATP-binding cassette B5 (ABCB5) transporter

The human ATP-binding cassette B5 (ABCB5) transporter, a member of the ABC transporter superfamily, is linked to chemoresistance in tumour cells by drug effluxion. However, little is known about its structure and drug-binding sites. In this study, we generated an atomistic model of the full-length human ABCB5 transporter with the highest quality using the X-ray crystal structure of mouse ABCB1 (Pgp1), a close homologue of ABCB5 and a well-studied member of the ABC family. Molecular dynamics simulations were used to validate the atomistic model of ABCB5 and characterise its structural properties in model cell membranes. Molecular docking simulations of known ABCB5 substrates such as taxanes, anthracyclines, camptothecin and etoposide were then used to identify at least three putative binding sites for chemotherapeutic drugs transported by ABCB5. The location of these three binding sites is predicted to overlap with the corresponding binding sites in Pgp1. These findings will serve as the basis for future in vitro studies to validate the nature of the identified substrate-binding sites in the full-length ABCB5 transporter.

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.

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.

The Role of a Crystallographically Unresolved Cytoplasmic Loop in Stabilizing the Bacterial Membrane Insertase YidC2

YidC, a bacterial member of the YidC/Alb3/Oxa1 insertase family, mediates membrane protein assembly and insertion. Cytoplasmic loops are known to have functional significance in membrane proteins such as YidC. Employing microsecond-level molecular dynamics (MD) simulations, we show that the crystallographically unresolved C2 loop plays a crucial role in the structural dynamics of Bacillus halodurans YidC2. We have modeled the C2 loop and used all- atom MD simulations to investigate the structural dynamics of YidC2 in its apo form, both with and without the C2 loop. The C2 loop was found to stabilize the entire protein and particularly the C1 region. C2 was also found to stabilize the alpha-helical character of the C-terminal region. Interestingly, the highly polar or charged lipid head groups of the simulated membranes were found to interact with and stabilize the C2 loop. These findings demonstrate that the crystallographically unresolved loops of membrane proteins could be important for the stabilization of the protein despite the apparent lack of structure, which could be due to the absence of the relevant lipids to stabilize them in crystallographic conditions.

An Effective Electric Dipole Model for Voltage-induced Gating Mechanism of Lysenin

Lysenin is a pore-forming toxin, which self-inserts open channels into sphingomyelin containing membranes and is known to be voltage regulated. The mechanistic details of its voltage gating mechanism, however, remains elusive despite much recent efforts. Here, we have employed a novel combination of experimental and computational techniques to examine a model for voltage gating, that is based on the existence of an “effective electric dipole” inspired by recent reported structures of lysenin. We support this mechanism by the observations that (i) the charge-reversal and neutralization substitutions in lysenin result in changing its electrical gating properties by modifying the strength of the dipole, and (ii) an increase in the viscosity of the solvent increases the drag force and slows down the gating. In addition, our molecular dynamics (MD) simulations of membrane-embedded lysenin provide a mechanistic picture for lysenin conformational changes, which reveals, for the first time, the existence of a lipid-dependent bulge region in the pore-forming module of lysenin, which may explain the gating mechanism of lysenin at a molecular level.

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.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.