ATP hydrolysis at one of the two sites in ABC transporters initiates transport related conformational transitions

ABC transporter inhibition by beauvericin partially overcomes drug resistance in Leishmania tropica

Leishmaniasis is a neglected tropical disease infecting the world's poorest populations. Miltefosine (ML) remains the primary oral drug against the cutaneous form of leishmaniasis. The ATP-binding cassette (ABC) transporters are key players in the xenobiotic efflux, and their inhibition could enhance the therapeutic index. In this study, the ability of beauvericin (BEA) to overcome ABC transporter-mediated resistance of Leishmania tropica to ML was assessed. In addition, the transcription profile of genes involved in resistance acquisition to ML was inspected. Finally, we explored the efflux mechanism of the drug and inhibitor. The efficacy of ML against all developmental stages of L. tropica in the presence or absence of BEA was evaluated using an absolute quantification assay. The expression of resistance genes was evaluated, comparing susceptible and resistant strains. Finally, the mechanisms governing the interaction between the ABC transporter and its ligands were elucidated using molecular docking and dynamic simulation. Relative quantification showed that the expression of the ABCG sub-family is mostly modulated by ML. In this study, we used BEA to impede resistance of Leishmania tropica. The IC50 values, following BEA treatment, were significantly reduced from 30.83, 48.17, and 16.83 µM using ML to 8.14, 11.1, and 7.18 µM when using a combinatorial treatment (ML + BEA) against promastigotes, axenic amastigotes, and intracellular amastigotes, respectively. We also demonstrated a favorable BEA-binding enthalpy to L. tropica ABC transporter compared to ML. Our study revealed that BEA partially reverses the resistance development of L. tropica to ML by blocking the alternate ATP hydrolysis cycle.

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.

Structure and Function of Hepatobiliary ATP Binding Cassette Transporters

The liver is beyond any doubt the most important metabolic organ of the human body. This function requires an intensive crosstalk within liver cellular structures, but also with other organs. Membrane transport proteins are therefore of upmost importance as they represent the sensors and mediators that shuttle signals from outside to the inside of liver cells and/or vice versa. In this review, we summarize the known literature of liver transport proteins with a clear emphasis on functional and structural information on ATP binding cassette (ABC) transporters, which are expressed in the human liver. These primary active membrane transporters form one of the largest families of membrane proteins. In the liver, they play an essential role in for example bile formation or xenobiotic export. Our review provides a state of the art and comprehensive summary of the current knowledge of hepatobiliary ABC transporters. Clearly, our knowledge has improved with a breath-taking speed over the last few years and will expand further. Thus, this review will provide the status quo and will lay the foundation for new and exciting avenues in liver membrane transporter research.

The transport pathway in the ABCG2 protein and its regulation revealed by molecular dynamics simulations

Atomic-level structural insight on the human ABCG2 membrane protein, a pharmacologically important transporter, has been recently revealed by several key papers. In spite of the wealth of structural data, the pathway of transmembrane movement for the large variety of structurally different ABCG2 substrates and the physiological lipid regulation of the transporter has not been elucidated. The complex molecular dynamics simulations presented here may provide a breakthrough in understanding the steps of the substrate transport process and its regulation by cholesterol. Our analysis revealed drug binding cavities other than the central binding site and delineated a putative dynamic transport pathway for substrates with variable structures. We found that membrane cholesterol accelerated drug transport by promoting the closure of cytoplasmic protein regions. Since ABCG2 is present in all major biological barriers and drug-metabolizing organs, influences the pharmacokinetics of numerous clinically applied drugs, and plays a key role in uric acid extrusion, this information may significantly promote a reliable prediction of clinically important substrate characteristics and drug-drug interactions.

ATP-binding Cassette Exporters: Structure and Mechanism with a Focus on P-glycoprotein and MRP1

BACKGROUND

Proteins that belong to the ATP-binding cassette superfamily include transporters that mediate the efflux of substrates from cells. Among these exporters, P-glycoprotein and MRP1 are involved in cancer multidrug resistance, protection from endo and xenobiotics, determination of drug pharmacokinetics, and the pathophysiology of a variety of disorders.

OBJECTIVE

To review the information available on ATP-binding cassette exporters, with a focus on Pglycoprotein, MRP1 and related proteins. We describe tissue localization and function of these transporters in health and disease, and discuss the mechanisms of substrate transport. We also correlate recent structural information with the function of the exporters, and discuss details of their molecular mechanism with a focus on the nucleotide-binding domains.

METHODS

Evaluation of selected publications on the structure and function of ATP-binding cassette proteins.

CONCLUSIONS

Conformational changes on the nucleotide-binding domains side of the exporters switch the accessibility of the substrate-binding pocket between the inside and outside, which is coupled to substrate efflux. However, there is no agreement on the magnitude and nature of the changes at the nucleotide- binding domains side that drive the alternate-accessibility. Comparison of the structures of Pglycoprotein and MRP1 helps explain differences in substrate selectivity and the bases for polyspecificity. P-glycoprotein substrates are hydrophobic and/or weak bases, and polyspecificity is explained by a flexible hydrophobic multi-binding site that has a few acidic patches. MRP1 substrates are mostly organic acids, and its polyspecificity is due to a single bipartite binding site that is flexible and displays positive charge.

Quantitative comparison of ABC membrane protein type I exporter structures in a standardized way

An increasing number of ABC membrane protein structures are determined by cryo-electron microscopy and X-ray crystallography, consequently identifying differences between their conformations has become an arising issue. Therefore, we propose to define standardized measures for ABC Type I exporter structure characterization. We set conformational vectors, conftors, which describe the relative orientation of domains and can highlight structural differences. In addition, continuum electrostatics calculations were performed to characterize the energetics of membrane insertion illuminating functionally crucial regions. In summary, the proposed metrics contribute to deeper understanding of ABC membrane proteins' structural features, structure validation, and analysis of movements observed in a molecular dynamics trajectory. Moreover, the concept of standardized metrics can be applied not only to ABC membrane protein structures (http://conftors.hegelab.org).

Comparison of mechanistic transport cycle models of ABC exporters

ABC (ATP binding cassette) transporters, ubiquitous in all kingdoms of life, carry out essential substrate transport reactions across cell membranes. Their transmembrane domains bind and translocate substrates and are connected to a pair of nucleotide binding domains, which bind and hydrolyze ATP to energize import or export of substrates. Over four decades of investigations into ABC transporters have revealed numerous details from atomic-level structural insights to their functional and physiological roles. Despite all these advances, a comprehensive understanding of the mechanistic principles of ABC transporter function remains elusive. The human multidrug resistance transporter ABCB1, also referred to as P-glycoprotein (P-gp), is one of the most intensively studied ABC exporters. Using ABCB1 as the reference point, we aim to compare the dominating mechanistic models of substrate transport and ATP hydrolysis for ABC exporters and to highlight the experimental and computational evidence in their support. In particular, we point out in silico studies that enhance and complement available biochemical data. “This article is part of a Special Issue entitled: Beyond the Structure Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.”

Structural basis for antibacterial peptide self-immunity by the bacterial ABC transporter McjD

Certain pathogenic bacteria produce and release toxic peptides to ensure either nutrient availability or evasion from the immune system. These peptides are also toxic to the producing bacteria that utilize dedicated ABC transporters to provide self‐immunity. The ABC transporter McjD exports the antibacterial peptide MccJ25 in Escherichia coli. Our previously determined McjD structure provided some mechanistic insights into antibacterial peptide efflux. In this study, we have determined its structure in a novel conformation, apo inward‐occluded and a new nucleotide‐bound state, high‐energy outward‐occluded intermediate state, with a defined ligand binding cavity. Predictive cysteine cross‐linking in E. coli membranes and PELDOR measurements along the transport cycle indicate that McjD does not undergo major conformational changes as previously proposed for multi‐drug ABC exporters. Combined with transport assays and molecular dynamics simulations, we propose a novel mechanism for toxic peptide ABC exporters that only requires the transient opening of the cavity for release of the peptide. We propose that shielding of the cavity ensures that the transporter is available to export the newly synthesized peptides, preventing toxic‐level build‐up.

Luminescence resonance energy transfer spectroscopy of ATP-binding cassette proteins

The ATP-binding cassette (ABC) superfamily includes regulatory and transport proteins. Most human ABC exporters pump substrates out of cells using energy from ATP hydrolysis. Although major advances have been made toward understanding the molecular mechanism of ABC exporters, there are still many issues unresolved. During the last few years, luminescence resonance energy transfer has been used to detect conformational changes in real time, with atomic resolution, in isolated ABC nucleotide binding domains (NBDs) and full-length ABC exporters. NBDs are particularly interesting because they provide the power stroke for substrate transport. Luminescence resonance energy transfer (LRET) is a spectroscopic technique that can provide dynamic information with atomic-resolution of protein conformational changes under physiological conditions. Using LRET, it has been shown that NBD dimerization, a critical step in ABC proteins catalytic cycle, requires binding of ATP to two nucleotide binding sites. However, hydrolysis at just one of the sites can drive dissociation of the NBD dimer. It was also found that the NBDs of the bacterial ABC exporter MsbA reconstituted in a lipid bilayer membrane and studied at 37°C never separate as much as suggested by crystal structures. This observation stresses the importance of performing structural/functional studies of ABC exporters under physiologic conditions. This article is part of a Special Issue entitled: Beyond the Structure-Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.

Unidirectional Transport Mechanism in an ATP Dependent Exporter

ATP-binding cassette (ABC) transporters use the energy of ATP binding and hydrolysis to move a large variety of compounds across biological membranes. P-glycoprotein, involved in multidrug resistance, is the most investigated eukaryotic family member. Although a large number of biochemical and structural approaches have provided important information, the conformational dynamics underlying the coupling between ATP binding/hydrolysis and allocrite transport remains elusive. To tackle this issue, we performed molecular dynamic simulations for different nucleotide occupancy states of Sav1866, a prokaryotic P-glycoprotein homologue. The simulations reveal an outward-closed conformation of the transmembrane domain that is stabilized by the binding of two ATP molecules. The hydrolysis of a single ATP leads the X-loop, a key motif of the ATP binding cassette, to interfere with the transmembrane domain and favor its outward-open conformation. Our findings provide a structural basis for the unidirectionality of transport in ABC exporters and suggest a ratio of one ATP hydrolyzed per transport cycle.

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.

ATP Hydrolysis Induced Conformational Changes in the Vitamin B12 Transporter BtuCD Revealed by MD Simulations

ATP binding cassette (ABC) transporters utilize the energy of ATP hydrolysis to uni-directionally transport substrates across cell membrane. ATP hydrolysis occurs at the nucleotide-binding domain (NBD) dimer interface of ABC transporters, whereas substrate translocation takes place at the translocation pathway between the transmembrane domains (TMDs), which is more than 30 angstroms away from the NBD dimer interface. This raises the question of how the hydrolysis energy released at NBDs is "transmitted" to trigger the conformational changes at TMDs. Using molecular dynamics (MD) simulations, we studied the post-hydrolysis state of the vitamin B12 importer BtuCD. Totally 3-μs MD trajectories demonstrate a predominantly asymmetric arrangement of the NBD dimer interface, with the ADP-bound site disrupted and the ATP-bound site preserved in most of the trajectories. TMDs response to ATP hydrolysis by separation of the L-loops and opening of the cytoplasmic gate II, indicating that hydrolysis of one ATP could facilitate substrate translocation by opening the cytoplasmic end of translocation pathway. It was also found that motions of the L-loops and the cytoplasmic gate II are coupled with each other through a contiguous interaction network involving a conserved Asn83 on the extended stretch preceding TM3 helix plus the cytoplasmic end of TM2/6/7 helix bundle. These findings entail a TMD-NBD communication mechanism for type II ABC importers.

Jump into a New Fold-A Homology Based Model for the ABCG2/BCRP Multidrug Transporter

ABCG2/BCRP is a membrane protein, involved in xenobiotic and endobiotic transport in key pharmacological barriers and drug metabolizing organs, in the protection of stem cells, and in multidrug resistance of cancer. Pharmacogenetic studies implicated the role of ABCG2 in response to widely used medicines and anticancer agents, as well as in gout. Its Q141K variant exhibits decreased functional expression thus increased drug accumulation and decreased urate secretion. Still, there has been no reliable molecular model available for this protein, as the published structures of other ABC transporters could not be properly fitted to the ABCG2 topology and experimental data. The recently published high resolution structure of a close homologue, the ABCG5-ABCG8 heterodimer, revealed a new ABC transporter fold, unique for ABCG proteins. Here we present a structural model of the ABCG2 homodimer based on this fold and detail the experimental results supporting this model. In order to describe the effect of mutations on structure and dynamics, and characterize substrate recognition and cholesterol regulation we performed molecular dynamics simulations using full length ABCG2 protein embedded in a membrane bilayer and in silico docking simulations. Our results show that in the Q141K variant the introduced positive charge diminishes the interaction between the nucleotide binding and transmembrane domains and the R482G variation alters the orientation of transmembrane helices. Moreover, the R482 position, which plays a role the substrate specificity of the transporter, is located in one of the substrate binding pockets identified by the in silico docking calculations. In summary, the ABCG2 model and in silico simulations presented here may have significant impact on understanding drug distribution and toxicity, as well as drug development against cancer chemotherapy resistance or gout.

Homology modelling of human P-glycoprotein

P-glycoprotein (P-gp) is an ATP-binding cassette transporter that exports a huge range of compounds out of cells and is thus one of the key proteins in conferring multi-drug resistance in cancer. Understanding how it achieves such a broad specificity and the series of conformational changes that allow export to occur form major, on-going, research objectives around the world. Much of our knowledge to date has been derived from mutagenesis and assay data. However, in recent years, there has also been great progress in structural biology and although the structure of human P-gp has not yet been solved, there are now a handful of related structures on which homology models can be built to aid in the interpretation of the vast amount of experimental data that currently exists. Many models for P-gp have been built with this aim, but the situation is complicated by the apparent flexibility of the system and by the fact that although many potential templates exist, there is large variation in the conformational state in which they have been crystallized. In this review, we summarize how homology modelling has been used in the past, how models are typically selected and finally illustrate how MD simulations can be used as a means to give more confidence about models that have been generated via this approach.

In silico pharmacology: Drug membrane partitioning and crossing

Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; β-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.

Access Path to the Ligand Binding Pocket May Play a Role in Xenobiotics Selection by AhR

Understanding of multidrug binding at the atomic level would facilitate drug design and strategies to modulate drug metabolism, including drug transport, oxidation, and conjugation. Therefore we explored the mechanism of promiscuous binding of small molecules by studying the ligand binding domain, the PAS-B domain of the aryl hydrocarbon receptor (AhR). Because of the low sequence identities of PAS domains to be used for homology modeling, structural features of the widely employed HIF-2α and a more recent suitable template, CLOCK were compared. These structures were used to build AhR PAS-B homology models. We performed molecular dynamics simulations to characterize dynamic properties of the PAS-B domain and the generated conformational ensembles were employed in in silico docking. In order to understand structural and ligand binding features we compared the stability and dynamics of the promiscuous AhR PAS-B to other PAS domains exhibiting specific interactions or no ligand binding function. Our exhaustive in silico binding studies, in which we dock a wide spectrum of ligand molecules to the conformational ensembles, suggest that ligand specificity and selection may be determined not only by the PAS-B domain itself, but also by other parts of AhR and its protein interacting partners. We propose that ligand binding pocket and access channels leading to the pocket play equally important roles in discrimination of endogenous molecules and xenobiotics.

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

ATP-binding cassette (ABC) transporters are constructed from some common structural units: the highly conserved nucleotide-binding domains (NBDs), which work as a nucleotide-dependent engine for driving substrate transport, the diverse transmembrane domains (TMDs), which create the translocation pathway, and the coupling helices (CHs), which are located at the NBD-TMD interface. Although the CHs are believed to be essential for NBD-TMD communication, their roles remain unclear. In this study, we performed enzyme assays and molecular dynamics (MD) simulations of the ABC transporter MsbA and two MsbA mutants in which the amino acid residues of one of the CHs were mutated to alanines: (i) wild type (Wt), (ii) CH1 mutant (Mt1), and (iii) CH2 mutant (Mt2). The experiments show that the CH2 mutation decreases the ATPase activity (kcat) compared with that of the Wt (a decrease of 32%), and a nearly equal degree of decrease in the ATP binding affinity (Km) was observed for both Mt1 and Mt2. The MD simulations successfully accounted for several structural and dynamical origins for these experimental observations. In addition, on the basis of collective motion and morphing analyses, we propose that the reverse-rotational motions and noddinglike motions between the NBDs and TMDs are indispensable for the conformational transition between the inward- and outward-facing conformations. In particular, CH2 is significantly important for the occurrence of the noddinglike motion. These findings provide important insights into the structure-function relationship of ABC transporters.

Structure and mechanism of ATP-dependent phospholipid transporters

BACKGROUND

ATP-binding cassette (ABC) transporters and P4-ATPases are two large and seemingly unrelated families of primary active pumps involved in moving phospholipids from one leaflet of a biological membrane to the other.

SCOPE OF REVIEW

This review aims to identify common mechanistic features in the way phospholipid flipping is carried out by two evolutionarily unrelated families of transporters.

MAJOR CONCLUSIONS

Both protein families hydrolyze ATP, although they employ different mechanisms to use it, and have a comparable size with twelve transmembrane segments in the functional unit. Further, despite differences in overall architecture, both appear to operate by an alternating access mechanism and during transport they might allow access of phospholipids to the internal part of the transmembrane domain. The latter feature is obvious for ABC transporters, but phospholipids and other hydrophobic molecules have also been found embedded in P-type ATPase crystal structures. Taken together, in two diverse groups of pumps, nature appears to have evolved quite similar ways of flipping phospholipids.

GENERAL SIGNIFICANCE

Our understanding of the structural basis for phospholipid flipping is still limited but it seems plausible that a general mechanism for phospholipid flipping exists in nature. This article is part of a Special Issue entitled Structural biochemistry and biophysics of membrane proteins.

Asymmetric perturbations of signalling oligomers

This review focuses on rapid and reversible noncovalent interactions for symmetric oligomers of signalling proteins. Symmetry mismatch, transient symmetry breaking and asymmetric perturbations via chemical (ligand binding) and physical (electric or mechanic) effects can initiate the signalling events. Advanced biophysical methods can reveal not only structural symmetries of stable membrane-bound signalling proteins but also asymmetric functional transition states. Relevant techniques amenable to distinguish between symmetric and asymmetric architectures are discussed including those with the capability of capturing low-populated transient conformational states. Typical examples of signalling proteins are overviewed for symmetry breaking in dimers (GPCRs, growth factor receptors, transcription factors); trimers (acid-sensing ion channels); tetramers (voltage-gated cation channels, ionotropic glutamate receptor, CNG and CHN channels); pentameric ligand-gated and mechanosensitive channels; higher order oligomers (gap junction channel, chaperonins, proteasome, virus capsid); as well as primary and secondary transporters. In conclusion, asymmetric perturbations seem to play important functional roles in a broad range of communicating networks.

Tertiary structure prediction and identification of druggable pocket in the cancer biomarker - Osteopontin-c

Background

Osteopontin (Eta, secreted sialoprotein 1, opn) is secreted from different cell types including cancer cells. Three splice variant forms namely osteopontin-a, osteopontin-b and osteopontin-c have been identified. The main astonishing feature is that osteopontin-c is found to be elevated in almost all types of cancer cells. This was the vital point to consider it for sequence analysis and structure predictions which provide ample chances for prognostic, therapeutic and preventive cancer research.

Methods

Osteopontin-c gene sequence was determined from Breast Cancer sample and was translated to protein sequence. It was then analyzed using various software and web tools for binding pockets, docking and druggability analysis. Due to the lack of homological templates, tertiary structure was predicted using ab-initio method server – I-TASSER and was evaluated after refinement using web tools. Refined structure was compared with known bone sialoprotein electron microscopic structure and docked with CD44 for binding analysis and binding pockets were identified for drug designing.

Results

Signal sequence of about sixteen amino acid residues was identified using signal sequence prediction servers. Due to the absence of known structures of similar proteins, three dimensional structure of osteopontin-c was predicted using I-TASSER server. The predicted structure was refined with the help of SUMMA server and was validated using SAVES server. Molecular dynamic analysis was carried out using GROMACS software. The final model was built and was used for docking with CD44. Druggable pockets were identified using pocket energies.

Conclusions

The tertiary structure of osteopontin-c was predicted successfully using the ab-initio method and the predictions showed that osteopontin-c is of fibrous nature comparable to firbronectin. Docking studies showed the significant similarities of QSAET motif in the interaction of CD44 and osteopontins between the normal and splice variant forms of osteopontins and binding pockets analyses revealed several pockets which paved the way to the identification of a druggable pocket.

Towards understanding promiscuity in multidrug efflux pumps

Drug export from cells is a major factor in the acquisition of cellular resistance to antimicrobial and cancer chemotherapy, and poses a significant threat to future clinical management of disease. Many of the proteins that catalyse drug efflux do so with remarkably low substrate specificity, a phenomenon known as multidrug transport. For these reasons we need a greater understanding of drug recognition and transport in multidrug pumps to inform research that attempts to circumvent their action. Structural and computational studies have been heralded as being great strides towards a full elucidation of multidrug recognition and transport. In this review we summarise these advances and ask how close we are to a molecular understanding of this remarkable phenomenon.

Nonequilibrium gating of CFTR on an equilibrium theme

Malfunction of cystic fibrosis transmembrane conductance regulator (CFTR), a member of the ABC protein superfamily that functions as an ATP-gated chloride channel, causes the lethal genetic disease, cystic fibrosis. This review focuses on the most recent findings on the gating mechanism of CFTR. Potential clinical relevance and implications to ABC transporter function are also discussed.

Mechanistic picture for conformational transition of a membrane transporter at atomic resolution

During their transport cycle, ATP-binding cassette (ABC) transporters undergo large-scale conformational changes between inward- and outward-facing states. Using an approach based on designing system-specific reaction coordinates and using nonequilibrium work relations, we have performed extensive all-atom molecular dynamics simulations in the presence of explicit membrane/solvent to sample a large number of mechanistically distinct pathways for the conformational transition of MsbA, a bacterial ABC exporter whose structure has been solved in multiple functional states. The computational approach developed here is based on (i) extensive exploration of system-specific biasing protocols (e.g., using collective variables designed based on available low-resolution crystal structures) and (ii) using nonequilibrium work relations for comparing the relevance of the transition pathways. The most relevant transition pathway identified using this approach involves several distinct stages reflecting the complex nature of the structural changes associated with the function of the protein. The opening of the cytoplasmic gate during the outward- to inward-facing transition of apo MsbA is found to be disfavored when the periplasmic gate is open and facilitated by a twisting motion of the nucleotide-binding domains that involves a dramatic change in their relative orientation. These results highlight the cooperativity between the transmembrane and the nucleotide-binding domains in the conformational transition of ABC exporters. The approach introduced here provides a framework to study large-scale conformational changes of other membrane transporters whose computational investigation at an atomic resolution may not be currently feasible using conventional methods.

Hydrolysis at one of the two nucleotide-binding sites drives the dissociation of ATP-binding cassette nucleotide-binding domain dimers

The functional unit of ATP-binding cassette (ABC) transporters consists of two transmembrane domains and two nucleotide-binding domains (NBDs). ATP binding elicits association of the two NBDs, forming a dimer in a head-to-tail arrangement, with two nucleotides "sandwiched" at the dimer interface. Each of the two nucleotide-binding sites is formed by residues from the two NBDs. We recently found that the prototypical NBD MJ0796 from Methanocaldococcus jannaschii dimerizes in response to ATP binding and dissociates completely following ATP hydrolysis. However, it is still unknown whether dissociation of NBD dimers follows ATP hydrolysis at one or both nucleotide-binding sites. Here, we used luminescence resonance energy transfer to study heterodimers formed by one active (donor-labeled) and one catalytically defective (acceptor-labeled) NBD. Rapid mixing experiments in a stop-flow chamber showed that NBD heterodimers with one functional and one inactive site dissociated at a rate indistinguishable from that of dimers with two hydrolysis-competent sites. Comparison of the rates of NBD dimer dissociation and ATP hydrolysis indicated that dissociation followed hydrolysis of one ATP. We conclude that ATP hydrolysis at one nucleotide-binding site drives NBD dimer dissociation.

Substrate versus inhibitor dynamics of P-glycoprotein

By far the most studied multidrug resistance protein is P-glycoprotein. Despite recent structural data, key questions about its function remain. P-glycoprotein (P-gp) is flexible and undergoes large conformational changes as part of its function and in this respect, details not only of the export cycle, but also the recognition stage are currently lacking. Given the flexibility, molecular dynamics (MD) simulations provide an ideal tool to examine this aspect in detail. We have performed MD simulations to examine the behaviour of P-gp. In agreement with previous reports, we found that P-gp undergoes large conformational changes which tended to result in the nucleotide-binding domains coming closer together. In all simulations, the approach of the NBDs was asymmetrical in agreement with previous observations for other ABC transporter proteins. To validate the simulations, we make extensive comparison to previous cross-linking data. Our results show very good agreement with the available data. We then went on to compare the influence of inhibitor compounds bound with simulations of a substrate (daunorubicin) bound. Our results suggest that inhibitors may work by keeping the NBDs apart, thus preventing ATP-hydrolysis. On the other hand, repeat simulations of daunorubicin (substrate) in one particular binding pose suggest that the approach of the NBDs is not impaired and that the structure would be still be competent to perform ATP hydrolysis, thus providing a model for inhibition or substrate transport. Finally we compare the latter to earlier QSAR data to provide a model for the first part of substrate transport within P-gp.

On the origin of large flexibility of P-glycoprotein in the inward-facing state

P-glycoprotein (Pgp) is one of the most biomedically relevant transporters in the ATP binding cassette (ABC) superfamily due to its involvement in developing multidrug resistance in cancer cells. Employing molecular dynamics simulations and double electron-electron resonance spectroscopy, we have investigated the structural dynamics of membrane-bound Pgp in the inward-facing state and found that Pgp adopts an unexpectedly wide range of conformations, highlighted by the degree of separation between the two nucleotide-binding domains (NBDs). The distance between the two NBDs in the equilibrium simulations covers a range of at least 20 Å, including, both, more open and more closed NBD configurations than the crystal structure. The double electron-electron resonance measurements on spin-labeled Pgp mutants also show wide distributions covering both longer and shorter distances than those observed in the crystal structure. Based on structural and sequence analyses, we propose that the transmembrane domains of Pgp might be more flexible than other structurally known ABC exporters. The structural flexibility of Pgp demonstrated here is not only in close agreement with, but also helps rationalize, the reported high NBD fluctuations in several ABC exporters and possibly represents a fundamental difference in the transport mechanism between ABC exporters and ABC importers. In addition, during the simulations we have captured partial entrance of a lipid molecule from the bilayer into the lumen of Pgp, reaching the putative drug binding site. The location of the protruding lipid suggests a putative pathway for direct drug recruitment from the membrane.

Visualizing functional motions of membrane transporters with molecular dynamics simulations

Computational modeling and molecular simulation techniques have become an integral part of modern molecular research. Various areas of molecular sciences continue to benefit from, indeed rely on, the unparalleled spatial and temporal resolutions offered by these technologies, to provide a more complete picture of the molecular problems at hand. Because of the continuous development of more efficient algorithms harvesting ever-expanding computational resources, and the emergence of more advanced and novel theories and methodologies, the scope of computational studies has expanded significantly over the past decade, now including much larger molecular systems and far more complex molecular phenomena. Among the various computer modeling techniques, the application of molecular dynamics (MD) simulation and related techniques has particularly drawn attention in biomolecular research, because of the ability of the method to describe the dynamical nature of the molecular systems and thereby to provide a more realistic representation, which is often needed for understanding fundamental molecular properties. The method has proven to be remarkably successful in capturing molecular events and structural transitions highly relevant to the function and/or physicochemical properties of biomolecular systems. Herein, after a brief introduction to the method of MD, we use a number of membrane transport proteins studied in our laboratory as examples to showcase the scope and applicability of the method and its power in characterizing molecular motions of various magnitudes and time scales that are involved in the function of this important class of membrane proteins.

Discrete molecular dynamics: an efficient and versatile simulation method for fine protein characterization

Until now it has been impractical to observe protein folding in silico for proteins larger than 50 residues. Limitations of both force field accuracy and computational efficiency make the folding problem very challenging. Here we employ discrete molecular dynamics (DMD) simulations with an all-atom force field to fold fast-folding proteins. We extend the DMD force field by introducing long-range electrostatic interactions to model salt-bridges and a sequence-dependent semiempirical potential accounting for natural tendencies of certain amino acid sequences to form specific secondary structures. We enhance the computational performance by parallelizing the DMD algorithm. Using a small number of commodity computers, we achieve sampling quality and folding accuracy comparable to the explicit-solvent simulations performed on high-end hardware. We demonstrate that DMD can be used to observe equilibrium folding of villin headpiece and WW domain, study two-state folding kinetics, and sample near-native states in ab initio folding of proteins of ∼100 residues.

Perspectives on the structure-function of ABC transporters: the Switch and Constant Contact models

ABC transporters constitute one of the largest protein families across the kingdoms of archaea, eubacteria and eukarya. They couple ATP hydrolysis to vectorial translocation of diverse substrates across membranes. The ABC transporter architecture comprises two transmembrane domains and two cytosolic ATP-binding cassettes. During 2002-2012, nine prokaryotic ABC transporter structures and two eukaryotic structures have been solved to medium resolution. Despite a wealth of biochemical, biophysical, and structural data, fundamental questions remain regarding the coupling of ATP hydrolysis to unidirectional substrate translocation, and the mechanistic suite of steps involved. The mechanics of the ATP cassette dimer is defined most popularly by the 'Switch Model', which proposes that hydrolysis in each protomer is sequential, and that as the sites are freed of nucleotide, the protomers lose contact across a large solvent-filled gap of 20-30 Å; as captured in several X-ray solved structures. Our 'Constant Contact' model for the operational mechanics of ATP binding and hydrolysis in the ATP-binding cassettes is derived from the 'alternating sites' model, proposed in 1995, and which requires an intrinsic asymmetry in the ATP sites, but does not require the partner protomers to lose contact. Thus one of the most debated issues regarding the function of ABC transporters is whether the cooperative mechanics of ATP hydrolysis requires the ATP cassettes to separate or remain in constant contact and this dilemma is discussed at length in this review.

Asymmetry in the homodimeric ABC transporter MsbA recognized by a DARPin

ABC transporters harness the energy from ATP binding and hydrolysis to translocate substrates across the membrane. Binding of two ATP molecules at the nucleotide binding domains (NBDs) leads to the formation of an outward-facing state. The conformational changes required to reset the transporter to the inward-facing state are initiated by sequential hydrolysis of the bound nucleotides. In a homodimeric ABC exporter such as MsbA responsible for lipid A transport in Escherichia coli, sequential ATP hydrolysis implies the existence of an asymmetric conformation. Here we report the in vitro selection of a designed ankyrin repeat protein (DARPin) specifically binding to detergent-solubilized MsbA. Only one DARPin binds to the homodimeric transporter in the absence as well as in the presence of nucleotides, suggesting that it recognizes asymmetries in MsbA. DARPin binding increases the rate of ATP hydrolysis by a factor of two independent of the substrate-induced ATPase stimulation. Electron paramagnetic resonance (EPR) measurements are found to be in good agreement with the available crystal structures and reveal that DARPin binding does not affect the large nucleotide-driven conformational changes of MsbA. The binding epitope was mapped by cross-linking and EPR to the membrane-spanning part of the transmembrane domain (TMD). Using cross-linked DARPin-MsbA complexes, 8-azido-ATP was found to preferentially photolabel one chain of the homodimer, suggesting that the asymmetries captured by DARPin binding at the TMDs are propagated to the NBDs. This work demonstrates that in vitro selected binders are useful tools to study the mechanism of membrane proteins.