Molecular dynamics simulations of the bacterial ABC transporter SAV1866 in the closed form

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.

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.

Dawn of a New Era for Membrane Protein Design

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

Dimethyl phthalate damages Staphylococcus aureus by changing the cell structure, inducing oxidative stress and inhibiting energy metabolism

Dimethyl phthalate (DMP), used as a plasticizer in industrial products, exists widely in air, water and soil. Staphylococcus aureus is a typical model organism representing Gram-positive bacteria. The molecular mechanisms of DMP toxicology in S. aureus were researched by proteomic and transcriptomic analyses. The results showed that the cell wall, membrane and cell surface characteristics were damaged and the growth was inhibited in S. aureus by DMP. Oxidative stress was induced by DMP in S. aureus. The activities of succinic dehydrogenase (SDH) and ATPase were changed by DMP, which could impact energy metabolism. Based on proteomic and transcriptomic analyses, the oxidative phosphorylation pathway was enhanced and the glycolysis/gluconeogenesis and pentose phosphate pathways were inhibited in S. aureus exposed to DMP. The results of real-time reverse transcription quantitative PCR (RT-qPCR) further confirmed the results of the proteomic and transcriptomic analyses. Lactic acid, pyruvic acid and glucose were reduced by DMP in S. aureus, which suggested that DMP could inhibit energy metabolism. The results indicated that DMP damaged the cell wall and membrane, induced oxidative stress, and inhibited energy metabolism and activation in S. aureus.

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.

Efflux mechanism and pathway of verapamil pumping by human P-glycoprotein

Multidrug resistance (MDR) caused by overexpressed permeability-glycoprotein (P-gp) in cancer cells is the main barrier for the cure of cancers. P-gp can pump many chemotherapeutic drugs, which is a viable target to overcome P-gp-mediated MDR by efficient inhibitors of P-gp. However, limited understanding of the efflux mechanism by human P-gp hinders the development of efficient inhibitors. Herein, the transport of a P-gp inhibitor, verapamil, by human P-gp has been investigated using targeted molecular dynamics simulations and energetics analysis based on our previous research on the transport of a drug (doxorubicin). The energetics analysis identifies that the driving forces for the transport of verapamil are electrostatic repulsions contributed by the positively charged residues in the initial stage and then hydrophobic interactions contributed by the important residues in the later stage. This scenario is generally consistent with that in the transport of doxorubicin. However, the positively charged residues and the important residues for the transport of verapamil are incompletely consistent with the relative residues for the transport of doxorubicin. Moreover, the binding free energy contributions of the positively charged residues for the transport of verapamil are generally higher than them for the transport of doxorubicin, while the important residues constitute significantly different binding free energy compositions in the transports of the two substrates. Consequently, the pathway for the transport of verapamil is identified, which shares only two residues (F336 and M986) with the pathway of doxorubicin. This may imply the weak competitiveness of verapamil with doxorubicin in the substrate efflux. Taken together, this work provided new insights into the efflux mechanisms by human P-gp and would be beneficial in the design of potent P-gp inhibitors.

ATP-Dependent Signaling in Simulations of a Revised Model of Cystic Fibrosis Transmembrane Conductance Regulator (CFTR)

Cystic fibrosis transmembrane conductance regulator (CFTR) is a member of the ATP-binding cassette (ABC) transporter superfamily that has uniquely evolved to function as a chloride channel. It binds and hydrolyzes ATP at its nucleotide binding domains to form a pore providing a diffusive pathway within its transmembrane domains. CFTR is the only known protein from the ABC superfamily with channel activity, and its dysfunction causes the disease cystic fibrosis. While much is known about the functional aspects of CFTR, significant gaps remain, such as the structure-function relationship underlying signaling of ATP binding. In the present work, we refined an existing homology model using an intermediate-resolution (9 Å) published cryo-electron microscopy map. The newly derived models have been simulated in equilibrium molecular dynamics simulations for a total of 2.5 μs in multiple ATP-occupancy states. Putative conformational movements connecting ATP binding with pore formation are elucidated and quantified. Additionally, new interdomain interactions between E543, K968, and K1292 have been identified and confirmed experimentally; these interactions may be relevant for signaling ATP binding and hydrolysis to the transmembrane domains and induction of pore opening.

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.

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

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

Response of Pseudomonas fluorescens to dimethyl phthalate

Dimethyl phthalate (DMP) is a ubiquitous pollutant that is very harmful to organisms due to its mutagenicity, teratogenicity and carcinogenicity. Pseudomonas fluorescens (P. fluorescens) is one of the most important bacteria in the environment. In this study, the response of P. fluorescens to DMP was investigated. It was found that DMP greatly inhibited the growth and glucose utilization of P. fluorescens when the concentration of DMP was ranged from 20 to 40 mg/l. The surface hydrophobicity and membrane permeability of P. fluorescens were also increased by DMP. DMP could lead to the deformations of cell membrane and the mis-opening of membrane channels. RNA-Seq and RT-qPCR results showed that the expression of some genes in P. fluorescens were altered, including the genes involved in energy metabolism, ATP-binding cassette (ABC) transporting and two-component systems. Additionally, the productions of lactic acid and pyruvic acid were reduced and the activity of hexokinase was inhibited in P. fluorescens by DMP. Clearly, the results suggested that DMP contamination could alter the biological function of P. fluorescens in the environment.

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.

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.

Maltose-binding protein effectively stabilizes the partially closed conformation of the ATP-binding cassette transporter MalFGK2

Maltose transporter MalFGK₂is a type-I importer in the ATP-binding cassette (ABC) transporter superfamily.

Conformational Dynamics and Protein-Substrate Interaction of ABC Transporter BtuCD at the Occluded State Revealed by Molecular Dynamics Simulations

ATP-binding cassette (ABC) transporters are ubiquitous in all three kingdoms of life and are implicated in many clinically relevant physiological processes. They couple the energy released by ATP hydrolysis to facilitate substrate translocation across cell membranes. The crystal structures of type II ABC importers have revealed their unique transmembrane domain architecture consisting of 10 transmembrane helices and their structurally conserved nucleotide-binding domains among all ABC transporters. However, molecular details of the interactions between the importers and their substrate remain largely elusive. Taking vitamin B₁₂ importer BtuCD as an exemplar of type II importers, we investigated the dynamics of its occluded state and the detailed protein-substrate interactions using molecular dynamics simulation. Our trajectories show that the importer accommodates the substrate through a nonspecific binding mode as the substrate undergoes evident vertical and tilt motions inside the translocation cavity. Extensive hydrogen bond and hydrophobic interactions were observed between the substrate and the importer; however, most of these interactions are weak, with <38% occurrence. The presence of substrate leads to enlargement of the translocation cavity, especially at its cytoplasmic end, which may activate cytoplasmic regions and probably facilitate the transportation. The perturbations caused by periplasmic binding protein and nucleotides were also investigated. The study provides deeper insight into the translocation mechanism of BtuCD.

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.

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.

Release of Entropic Spring Reveals Conformational Coupling Mechanism in the ABC Transporter BtuCD-F

Substrate translocation by ATP-binding cassette (ABC) transporters involves coupling of ATP binding and hydrolysis in the nucleotide-binding domains (NBDs) to conformational changes in the transmembrane domains. We used molecular dynamics simulations to investigate the atomic-level mechanism of conformational coupling in the ABC transporter BtuCD-F, which imports vitamin B12 across the inner membrane of Escherichia coli. Our simulations show how an engineered disulfide bond across the NBD dimer interface reduces conformational fluctuations and hence configurational entropy. As a result, the disulfide bond is under substantial mechanical stress. Releasing this entropic spring, as is the case in the wild-type transporter, combined with analyzing the pairwise forces between individual residues, unravels the coupling mechanism. The identified pathways along which force is propagated from the NBDs via the coupling helix to the transmembrane domains are composed of highly conserved residues, underlining their functional relevance. This study not only reveals the details of conformational coupling in BtuCD-F, it also provides a promising approach to other long-range conformational couplings, e.g., in ABC exporters or other ATP-driven molecular machines.

Characterizing the binding interactions between P-glycoprotein and eight known cardiovascular transport substrates

The multidrug efflux pump P-glycoprotein (Pgp) is upregulated in cardiomyocytes following chronic ischemia from infarction and hypoxia caused by sleep apnea. This report summarizes the molecular dynamic studies performed on eight cardiovascular drugs to determine their corresponding binding sites on mouse Pgp. Selected Pgp transport ligands include: Amiodarone, Bepridil, Diltiazem, Dipyridamole, Nicardipine, Nifedipine, Propranolol, and Quinidine. Extensive molecular dynamic equilibration simulations were performed to determine drug docking interactions. Distinct binding sites were not observed, but rather a binding belt was seen with multiple residues playing a role in each studied drug's stable docking. Three key drug–protein interactions were identified: hydrogen bonding, hydrophobic packing, and the formation of a “cage” of aromatic residues around the drug. After drug stabilization, water molecules were observed to leak into the binding belt and condense around the drug. Water influx into the binding domain of Pgp may play a role in catalytic transition and drug expulsion. The cytoplasmic recruitment theory was also tested, and the drugs were observed to interact with conserved loops of residues with a strong affinity. A free energy change of astronomical value is required to recruit the drug from the cytoplasm to the binding belt within the transmembrane domain of Pgp.

Molecular insight into conformational transmission of human P-glycoprotein

P-glycoprotein (P-gp), a kind of ATP-binding cassette transporter, can export candidates through a channel at the two transmembrane domains (TMDs) across the cell membranes using the energy released from ATP hydrolysis at the two nucleotide-binding domains (NBDs). Considerable evidence has indicated that human P-gp undergoes large-scale conformational changes to export a wide variety of anti-cancer drugs out of the cancer cells. However, molecular mechanism of the conformational transmission of human P-gp from the NBDs to the TMDs is still unclear. Herein, targeted molecular dynamics simulations were performed to explore the atomic detail of the conformational transmission of human P-gp. It is confirmed that the conformational transition from the inward- to outward-facing is initiated by the movement of the NBDs. It is found that the two NBDs move both on the two directions (x and y). The movement on the x direction leads to the closure of the NBDs, while the movement on the y direction adjusts the conformations of the NBDs to form the correct ATP binding pockets. Six key segments (KSs) protruding from the TMDs to interact with the NBDs are identified. The relative movement of the KSs along the y axis driven by the NBDs can be transmitted through α-helices to the rest of the TMDs, rendering the TMDs to open towards periplasm in the outward-facing conformation. Twenty eight key residue pairs are identified to participate in the interaction network that contributes to the conformational transmission from the NBDs to the TMDs of human P-gp. In addition, 9 key residues in each NBD are also identified. The studies have thus provided clear insight into the conformational transmission from the NBDs to the TMDs in human P-gp.

Mechanisms of CFTR functional variants that impair regulated bicarbonate permeation and increase risk for pancreatitis but not for cystic fibrosis

CFTR is a dynamically regulated anion channel. Intracellular WNK1-SPAK activation causes CFTR to change permeability and conductance characteristics from a chloride-preferring to bicarbonate-preferring channel through unknown mechanisms. Two severe CFTR mutations (CFTRsev) cause complete loss of CFTR function and result in cystic fibrosis (CF), a severe genetic disorder affecting sweat glands, nasal sinuses, lungs, pancreas, liver, intestines, and male reproductive system. We hypothesize that those CFTR mutations that disrupt the WNK1-SPAK activation mechanisms cause a selective, bicarbonate defect in channel function (CFTRBD) affecting organs that utilize CFTR for bicarbonate secretion (e.g. the pancreas, nasal sinus, vas deferens) but do not cause typical CF. To understand the structural and functional requirements of the CFTR bicarbonate-preferring channel, we (a) screened 984 well-phenotyped pancreatitis cases for candidate CFTRBD mutations from among 81 previously described CFTR variants; (b) conducted electrophysiology studies on clones of variants found in pancreatitis but not CF; (c) computationally constructed a new, complete structural model of CFTR for molecular dynamics simulation of wild-type and mutant variants; and (d) tested the newly defined CFTRBD variants for disease in non-pancreas organs utilizing CFTR for bicarbonate secretion. Nine variants (CFTR R74Q, R75Q, R117H, R170H, L967S, L997F, D1152H, S1235R, and D1270N) not associated with typical CF were associated with pancreatitis (OR 1.5, p = 0.002). Clones expressed in HEK 293T cells had normal chloride but not bicarbonate permeability and conductance with WNK1-SPAK activation. Molecular dynamics simulations suggest physical restriction of the CFTR channel and altered dynamic channel regulation. Comparing pancreatitis patients and controls, CFTRBD increased risk for rhinosinusitis (OR 2.3, p<0.005) and male infertility (OR 395, p<<0.0001). WNK1-SPAK pathway-activated increases in CFTR bicarbonate permeability are altered by CFTRBD variants through multiple mechanisms. CFTRBD variants are associated with clinically significant disorders of the pancreas, sinuses, and male reproductive system.

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.

Modeling the conformational changes underlying channel opening in CFTR

Mutations in the gene encoding the cystic fibrosis transmembrane conductance regulator protein (CFTR) cause cystic fibrosis (CF), the most common life-shortening genetic disease among Caucasians. Although general features of the structure of CFTR have been predicted from homology models, the conformational changes that result in channel opening and closing have yet to be resolved. We created new closed- and open-state homology models of CFTR, and performed targeted molecular dynamics simulations of the conformational transitions in a channel opening event. The simulations predict a conformational wave that starts at the nucleotide binding domains and ends with the formation of an open conduction pathway. Changes in side-chain interactions are observed in all major domains of the protein, and experimental confirmation was obtained for a novel intra-protein salt bridge that breaks near the end of the transition. The models and simulation add to our understanding of the mechanism of ATP-dependent gating in this disease-relevant ion channel.

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.

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.