Dimer opening of the nucleotide binding domains of ABC transporters after ATP hydrolysis

Atomic-level characterization of the conformational transition pathways in SARS-CoV-1 and SARS-CoV-2 spike proteins

Severe acute respiratory syndrome (SARS) coronaviruses 1 and 2 (SARS-CoV-1 and SARS-CoV-2) derive transmissibility from spike protein activation in the receptor binding domain (RBD) and binding to the host cell angiotensin converting enzyme 2 (ACE2). However, the mechanistic details that describe the large-scale conformational changes associated with spike protein activation or deactivation are still somewhat unknown. Here, we have employed an extensive set of nonequilibrium all-atom molecular dynamics (MD) simulations, utilizing a novel protocol, for the SARS-CoV-1 (CoV-1) and SARS-CoV-2 (CoV-2) prefusion spike proteins in order to characterize the conformational pathways associated with the active-to-inactive transition. Our results indicate that both CoV-1 and CoV-2 spike proteins undergo conformational transitions along pathways unique to each protein. We have identified a number of key residues that form various inter-domain saltbridges, suggesting a multi-stage conformational change along the pathways. We have also constructed the free energy profiles along the transition pathways for both CoV-1 and CoV-2 spike proteins. The CoV-2 spike protein must overcome larger free energy barriers to undergo conformational changes towards protein activation or deactivation, when compared to CoV-1.

ATP hydrolysis and nucleotide exit enhance maltose translocation in the MalFGK2E importer

ATP binding cassette (ABC) transporters employ ATP hydrolysis to harness substrate translocation across membranes. The Escherichia coli MalFGK₂E maltose importer is an example of a type I ABC importer and a model system for this class of ABC transporters. The MalFGK₂E importer is responsible for the intake of malto-oligossacharides in E.coli. Despite being extensively studied, little is known about the effect of ATP hydrolysis and nucleotide exit on substrate transport. In this work, we studied this phenomenon using extensive molecular dynamics simulations (MD) along with potential of mean force calculations of maltose transport across the pore, in the pre-hydrolysis, post-hydrolysis and nucleotide-free states. We concluded that ATP hydrolysis and nucleotide exit trigger conformational changes that result in the decrease of energetic barriers to maltose translocation towards the cytoplasm, with a concomitant increase of the energy barrier in the periplasmic side of the pore, contributing for the irreversibility of the process. We also identified key residues that aid in positioning and orientation of maltose, as well as a novel binding pocket for maltose in MalG. Additionally, ATP hydrolysis leads to conformations similar to the nucleotide-free state. This study shows the contribution of ATP hydrolysis and nucleotide exit in the transport cycle, shedding light on ABC type I importer mechanisms.

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.

Conformational changes in the nucleotide-binding domains of P-glycoprotein induced by ATP hydrolysis

P-glycoprotein (Pgp) is a member of the ABC transporter superfamily with high physiological importance. Pgp nucleotide-binding domains (NBDs) drive the transport cycle through ATP binding and hydrolysis. We use molecular dynamics simulations to investigate the ATP hydrolysis-induced conformational changes in NBDs. Five systems, including all possible ATP/ADP combinations in the NBDs and the APO system, are simulated. ATP/ADP exchange induces conformational changes mostly within the conserved signature motif of the NBDs, resulting in relative orientational changes in the NBDs. Nucleotide removal leads to additional orientational changes in the NBDs, allowing their dissociation. Furthermore, we capture putative hydrolysis-competent configurations in which the conserved glutamate in the Walker-B motif acts as a catalytic base capturing a water molecule likely initiating ATP hydrolysis.

Use of Molecular Dynamics Simulations in Structure-Based Drug Discovery

BACKGROUND

Traditional drug discovery is a lengthy process which involves a huge amount of resources. Modern-day drug discovers various multidisciplinary approaches amongst which, computational ligand and structure-based drug designing methods contribute significantly. Structure-based drug designing techniques require the knowledge of structural information of drug target and drug-target complexes. Proper understanding of drug-target binding requires the flexibility of both ligand and receptor to be incorporated. Molecular docking refers to the static picture of the drug-target complex(es). Molecular dynamics, on the other hand, introduces flexibility to understand the drug binding process.

OBJECTIVE

The aim of the present study is to provide a systematic review on the usage of molecular dynamics simulations to aid the process of structure-based drug design.

METHOD

This review discussed findings from various research articles and review papers on the use of molecular dynamics in drug discovery. All efforts highlight the practical grounds for which molecular dynamics simulations are used in drug designing program. In summary, various aspects of the use of molecular dynamics simulations that underline the basis of studying drug-target complexes were thoroughly explained.

RESULTS

This review is the result of reviewing more than a hundred papers. It summarizes various problems that use molecular dynamics simulations.

CONCLUSION

The findings of this review highlight how molecular dynamics simulations have been successfully implemented to study the structure-function details of specific drug-target complexes. It also identifies the key areas such as stability of drug-target complexes, ligand binding kinetics and identification of allosteric sites which have been elucidated using molecular dynamics simulations.

Conversion of chemical to mechanical energy by the nucleotide binding domains of ABCB1

P-glycoprotein (ABCB1) is an important component of barrier tissues that extrudes a wide range of chemically unrelated compounds. ABCB1 consists of two transmembrane domains forming the substrate binding and translocation domain, and of two cytoplasmic nucleotide binding domains (NBDs) that provide the energy by binding and hydrolyzing ATP. We analyzed the mechanistic and energetic properties of the NBD dimer via molecular dynamics simulations. We find that MgATP stabilizes the NBD dimer through strong attractive forces by serving as an interaction hub. The irreversible ATP hydrolysis step converts the chemical energy stored in the phosphate bonds of ATP into potential energy. Following ATP hydrolysis, interactions between the NBDs and the ATP hydrolysis products MgADP + Pi remain strong, mainly because Mg2+ forms stabilizing interactions with ADP and Pi. Despite these stabilizing interactions MgADP + Pi are unable to hold the dimer together, which becomes separated by avid interactions of MgADP + Pi with water. ATP binding to the open NBDs and ATP hydrolysis in the closed NBD dimer represent two steps of energy input, each leading to the formation of a high energy state. Relaxation from these high energy states occurs through conformational changes that push ABCB1 through the transport cycle.

An integrated transport mechanism of the maltose ABC importer

ATP-binding cassette (ABC) transporters use the energy of ATP hydrolysis to transport a large diversity of molecules actively across biological membranes. A combination of biochemical, biophysical, and structural studies has established the maltose transporter MalFGK2 as one of the best characterized proteins of the ABC family. MalF and MalG are the transmembrane domains, and two MalKs form a homodimer of nucleotide-binding domains. A periplasmic maltose-binding protein (MalE) delivers maltose and other maltodextrins to the transporter, and triggers its ATPase activity. Substrate import occurs in a unidirectional manner by ATP-driven conformational changes in MalK2 that allow alternating access of the substrate-binding site in MalF to each side of the membrane. In this review, we present an integrated molecular mechanism of the transport process considering all currently available information. Furthermore, we summarize remaining inconsistencies and outline possible future routes to decipher the full mechanistic details of transport by MalEFGK2 complex and that of related importer systems.

F508del disturbs the dynamics of the nucleotide binding domains of CFTR before and after ATP hydrolysis

The cystic fibrosis transmembrane conductance regulator (CFTR) channel is an ion channel responsible for chloride transport in epithelia and it belongs to the class of ABC transporters. The deletion of phenylalanine 508 (F508del) in CFTR is the most common mutation responsible for cystic fibrosis. Little is known about the effect of the mutation in the isolated nucleotide binding domains (NBDs), on dimer dynamics, ATP hydrolysis and even on nucleotide binding. Using molecular dynamics simulations of the human CFTR NBD dimer, we showed that F508del increases, in the prehydrolysis state, the inter-motif distance in both ATP binding sites (ABP) when ATP is bound. Additionally, a decrease in the number of catalytically competent conformations was observed in the presence of F508del. We used the subtraction technique to study the first 300 ps after ATP hydrolysis in the catalytic competent site and found that the F508del dimer evidences lower conformational changes than the wild type. Using longer simulation times, the magnitude of the conformational changes in both forms increases. Nonetheless, the F508del dimer shows lower C-α RMS values in comparison to the wild-type, on the F508del loop, on the residues surrounding the catalytic site and the portion of NBD2 adjacent to ABP1. These results provide evidence that F508del interferes with the NBD dynamics before and after ATP hydrolysis. These findings shed a new light on the effect of F508del on NBD dynamics and reveal a novel mechanism for the influence of F508del on CFTR.

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.

Dissecting the Forces that Dominate Dimerization of the Nucleotide Binding Domains of ABCB1

P-glycoprotein, also known as multidrug resistance protein 1 or ABCB1, can export a wide range of chemically unrelated compounds, including chemotherapeutic drugs. ABCB1 consists of two transmembrane domains that form the substrate binding and translocation domain, and of two cytoplasmic nucleotide binding domains (NBDs) that energize substrate transport by ATP binding and hydrolysis. ATP binding triggers dimerization of the NBDs, which switches the transporter from an inward facing to an outward facing transmembrane domain conformation. We performed MD simulations to study the dynamic behavior of the NBD dimer in the presence or absence of nucleotides. In the apo configuration, the NBDs were overall attractive to each other as shown in the potential of mean force profile, but the energy well was shallow and broad. In contrast, a sharp and deep energy minimum (∼-42 kJ/mol) was found in the presence of ATP, leading to a well-defined conformation. Motif interaction network analyses revealed that ATP stabilizes the NBD dimer by serving as the central hub for interdomain connections. Simulations showed that forces promoting dimerization are multilayered, dominated by electrostatic interactions between the nucleotide and conserved amino acids of the signature sequence and the Walker A motif. In addition, direct and water-bridged hydrogen bonds between NBDs provided conformation-defining interactions. Importantly, we characterized a largely unrecognized but essential contribution from hydrophobic interactions between the adenine moiety of the nucleotides and a hydrophobic surface of the X-loop to the stabilization of the nucleotide-bound NBD dimer. These hydrophobic interactions lead to a sharp energy minimum, thereby conformationally restricting the nucleotide-bound state.

Allosteric effects of ATP binding on the nucleotide-binding domain of a heterodimeric ATP-binding cassette transporter

ATP-binding cassette (ABC) exporters mediate vital transport of a variety of molecules across the lipid bilayer in all organisms. To explore the allosteric effect of ATP binding at the asymmetric ATPase sites, molecular dynamics simulations were performed on the nucleotide-binding domains (NBDs) of a heterodimeric exporter TM287/288 in 4 different ATP-bound states. The results showed that ATP bound at the degenerate site can maintain a semi-open conformation of NBD1-NBD2, which may be defective in ATP hydrolysis. By contrast, when bound at the consensus site, ATP can induce an intra-domain rotation of the α-helical subdomain towards the RecA-like subdomain of NBD2 at the degenerate site. The rotation of the α-helical subdomain rearranged the hydrogen bond networks at the NBD1-NBD2 interface, induced a significant conformational change in the D-loop at the degenerate site and inter- and intra-domain communications at both sites, and eventually elicited dimerization of NBD1-NBD2. These findings indicate that the asymmetric ATPase sites of the heterodimeric exporter are structurally and functionally distinct.

Cell envelope stress in mycobacteria is regulated by the novel signal transduction ATPase IniR in response to trehalose

The cell envelope of mycobacteria is a highly unique and complex structure that is functionally equivalent to that of Gram-negative bacteria to protect the bacterial cell. Defects in the integrity or assembly of this cell envelope must be sensed to allow the induction of stress response systems. The promoter that is specifically and most strongly induced upon exposure to ethambutol and isoniazid, first line drugs that affect cell envelope biogenesis, is the iniBAC promoter. In this study, we set out to identify the regulator of the iniBAC operon in Mycobacterium marinum using an unbiased transposon mutagenesis screen in a constitutively iniBAC-expressing mutant background. We obtained multiple mutants in the mce1 locus as well as mutants in an uncharacterized putative transcriptional regulator (MMAR_0612). This latter gene was shown to function as the iniBAC regulator, as overexpression resulted in constitutive iniBAC induction, whereas a knockout mutant was unable to respond to the presence of ethambutol and isoniazid. Experiments with the M. tuberculosis homologue (Rv0339c) showed identical results. RNAseq experiments showed that this regulatory gene was exclusively involved in the regulation of the iniBAC operon. We therefore propose to name this dedicated regulator iniBAC Regulator (IniR). IniR belongs to the family of signal transduction ATPases with numerous domains, including a putative sugar-binding domain. Upon testing different sugars, we identified trehalose as an activator and metabolic cue for iniBAC activation, which could also explain the effect of the mce1 mutations. In conclusion, cell envelope stress in mycobacteria is regulated by IniR in a cascade that includes trehalose.

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.

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.”

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.

Coupling between ATP hydrolysis and protein conformational change in maltose transporter

As the intracellular part of maltose transporter, MalK dimer utilizes the energy of ATP hydrolysis to drive protein conformational change, which then facilitates substrate transport. Free energy evaluation of the complete conformational change before and after ATP hydrolysis is helpful to elucidate the mechanism of chemical-to-mechanical energy conversion in MalK dimer, but is lacking in previous studies. In this work, we used molecular dynamics simulations to investigate the structural transition of MalK dimer among closed, semi-open and open states. We observed spontaneous structural transition from closed to open state in the ADP-bound system and partial closure of MalK dimer from the semi-open state in the ATP-bound system. Subsequently, we calculated the reaction pathways connecting the closed and open states for the ATP- and ADP-bound systems and evaluated the free energy profiles along the paths. Our results suggested that the closed state is stable in the presence of ATP but is markedly destabilized when ATP is hydrolyzed to ADP, which thus explains the coupling between ATP hydrolysis and protein conformational change of MalK dimer in thermodynamics. Proteins 2017; 85:207-220. © 2016 Wiley Periodicals, Inc.

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.

Toward Determining ATPase Mechanism in ABC Transporters: Development of the Reaction Path-Force Matching QM/MM Method

Adenosine triphosphate (ATP)-binding cassette (ABC) transporters are ubiquitous ATP-dependent membrane proteins involved in translocations of a wide variety of substrates across cellular membranes. To understand the chemomechanical coupling mechanism as well as functional asymmetry in these systems, a quantitative description of how ABC transporters hydrolyze ATP is needed. Complementary to experimental approaches, computer simulations based on combined quantum mechanical and molecular mechanical (QM/MM) potentials have provided new insights into the catalytic mechanism in ABC transporters. Quantitatively reliable determination of the free energy requirement for enzymatic ATP hydrolysis, however, requires substantial statistical sampling on QM/MM potential. A case study shows that brute force sampling of ab initio QM/MM (AI/MM) potential energy surfaces is computationally impractical for enzyme simulations of ABC transporters. On the other hand, existing semiempirical QM/MM (SE/MM) methods, although affordable for free energy sampling, are unreliable for studying ATP hydrolysis. To close this gap, a multiscale QM/MM approach named reaction path-force matching (RP-FM) has been developed. In RP-FM, specific reaction parameters for a selected SE method are optimized against AI reference data along reaction paths by employing the force matching technique. The feasibility of the method is demonstrated for a proton transfer reaction in the gas phase and in solution. The RP-FM method may offer a general tool for simulating complex enzyme systems such as ABC transporters.

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.

Conformational Changes of the Antibacterial Peptide ATP Binding Cassette Transporter McjD Revealed by Molecular Dynamics Simulations

The ATP binding cassette (ABC) transporters form one of the largest protein superfamilies. They use the energy of ATP hydrolysis to transport chemically diverse ligands across membranes. An alternating access mechanism in which the transporter switches between inward- and outward-facing conformations has been proposed to describe the translocation process. One of the main open questions in this process is the degree of opening of the transporter at different stages of the transport cycle, as crystal structures and biochemical data have suggested a wide range of distances between nucleotide binding domains. Recently, the crystal structure of McjD, an antibacterial peptide ABC transporter from Escherichia coli, revealed a new occluded intermediate state of the transport cycle. The transmembrane domain is closed on both sides of the membrane, forming a cavity that can accommodate its ligand, MccJ25, a lasso peptide of 21 amino acids. In this work, we investigate the degree of opening of the transmembrane cavity required for ligand translocation. By means of steered molecular dynamics simulations, the ligand was pulled from the internal cavity to the extracellular side. This resulted in an outward-facing state. Comparison with existing outward-facing crystal structures shows a smaller degree of opening in the simulations, suggesting that the large conformational changes in some crystal structures may not be necessary even for a large substrate like MccJ25.

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

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

The Nucleotide-Free State of the Multidrug Resistance ABC Transporter LmrA: Sulfhydryl Cross-Linking Supports a Constant Contact, Head-to-Tail Configuration of the Nucleotide-Binding Domains

ABC transporters are integral membrane pumps that are responsible for the import or export of a diverse range of molecules across cell membranes. ABC transporters have been implicated in many phenomena of medical importance, including cystic fibrosis and multidrug resistance in humans. The molecular architecture of ABC transporters comprises two transmembrane domains and two ATP-binding cassettes, or nucleotide-binding domains (NBDs), which are highly conserved and contain motifs that are crucial to ATP binding and hydrolysis. Despite the improved clarity of recent structural, biophysical, and biochemical data, the seemingly simple process of ATP binding and hydrolysis remains controversial, with a major unresolved issue being whether the NBD protomers separate during the catalytic cycle. Here chemical cross-linking data is presented for the bacterial ABC multidrug resistance (MDR) transporter LmrA. These indicate that in the absence of nucleotide or substrate, the NBDs come into contact to a significant extent, even at 4°C, where ATPase activity is abrogated. The data are clearly not in accord with an inward-closed conformation akin to that observed in a crystal structure of V. cholerae MsbA. Rather, they suggest a head-to-tail configuration ‘sandwich’ dimer similar to that observed in crystal structures of nucleotide-bound ABC NBDs. We argue the data are more readily reconciled with the notion that the NBDs are in proximity while undergoing intra-domain motions, than with an NBD ‘Switch’ mechanism in which the NBD monomers separate in between ATP hydrolysis cycles.

Identification of Possible Binding Sites for Morphine and Nicardipine on the Multidrug Transporter P-Glycoprotein Using Umbrella Sampling Techniques

The multidrug transporter P-glycoprotein (P-gp) is central to the development of multidrug resistance in cancer. While residues essential for transport and binding have been identified, the location, composition, and specificity of potential drug binding sites are uncertain. Here molecular dynamics simulations are used to calculate the free energy profile for the binding of morphine and nicardipine to P-gp. We show that morphine and nicardipine primarily interact with key residues implicated in binding and transport from mutational studies, binding at different but overlapping sites within the transmembrane pore. Their permeation pathways were distinct but involved overlapping sets of residues. The results indicate that the binding location and permeation pathways of morphine and nicardipine are not well separated and cannot be considered as unique. This has important implications for our understanding of substrate uptake and transport by P-gp. Our results are independent of the choice of starting structure and consistent with a range of experimental studies.

ATP-induced conformational changes of nucleotide-binding domains in an ABC transporter. Importance of the water-mediated entropic force

ATP binding cassette (ABC) proteins belong to a superfamily of active transporters. Recent experimental and computational studies have shown that binding of ATP to the nucleotide binding domains (NBDs) of ABC proteins drives the dimerization of NBDs, which, in turn, causes large conformational changes within the transmembrane domains (TMDs). To elucidate the active substrate transport mechanism of ABC proteins, it is first necessary to understand how the NBD dimerization is driven by ATP binding. In this study, we selected MalKs (NBDs of a maltose transporter) as a representative NBD and calculated the free-energy change upon dimerization using molecular mechanics calculations combined with a statistical thermodynamic theory of liquids, as well as a method to calculate the translational, rotational, and vibrational entropy change. This combined method is applied to a large number of snapshot structures obtained from molecular dynamics simulations containing explicit water molecules. The results suggest that the NBD dimerization proceeds with a large gain of water entropy when ATP molecules bind to the NBDs. The energetic gain arising from direct NBD-NBD interactions is canceled by the dehydration penalty and the configurational-entropy loss. ATP hydrolysis induces a loss of the shape complementarity between the NBDs, which leads to the dissociation of the dimer, due to a decrease in the water-entropy gain and an increase in the configurational-entropy loss. This interpretation of the NBD dimerization mechanism in concert with ATP, especially focused on the water-mediated entropy force, is potentially applicable to a wide variety of the ABC transporters.

In silico screening for inhibitors of p-glycoprotein that target the nucleotide binding domains

Multidrug resistances and the failure of chemotherapies are often caused by the expression or overexpression of ATP-binding cassette transporter proteins such as the multidrug resistance protein, P-glycoprotein (P-gp). P-gp is expressed in the plasma membrane of many cell types and protects cells from accumulation of toxins. P-gp uses ATP hydrolysis to catalyze the transport of a broad range of mostly hydrophobic compounds across the plasma membrane and out of the cell. During cancer chemotherapy, the administration of therapeutics often selects for cells which overexpress P-gp, thereby creating populations of cancer cells resistant to a variety of chemically unrelated chemotherapeutics. The present study describes extremely high-throughput, massively parallel in silico ligand docking studies aimed at identifying reversible inhibitors of ATP hydrolysis that target the nucleotide-binding domains of P-gp. We used a structural model of human P-gp that we obtained from molecular dynamics experiments as the protein target for ligand docking. We employed a novel approach of subtractive docking experiments that identified ligands that bound predominantly to the nucleotide-binding domains but not the drug-binding domains of P-gp. Four compounds were found that inhibit ATP hydrolysis by P-gp. Using electron spin resonance spectroscopy, we showed that at least three of these compounds affected nucleotide binding to the transporter. These studies represent a successful proof of principle demonstrating the potential of targeted approaches for identifying specific inhibitors of P-gp.

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.

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.

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.

Simulation studies of the mechanism of membrane transporters

Membrane transporters facilitate active transport of their specific substrates, often against their electrochemical gradients across the membrane, through coupling the process to various sources of cellular energy, for example, ATP binding and hydrolysis in primary transporters, and pre-established electrochemical gradient of molecular species other than the substrate in the case of secondary transporters. In order to provide efficient energy-coupling mechanisms, membrane transporters have evolved into molecular machines in which stepwise binding, translocation, and transformation of various molecular species are closely coupled to protein conformational changes that take the transporter from one functional state to another during the transport cycle. Furthermore, in order to prevent the formation of leaky states and to be able to pump the substrate against its electrochemical gradient, all membrane transporters use the widely-accepted "alternating access mechanism," which ensures that the substrate is only accessible from one side of the membrane at a given time, but relies on complex and usually global protein conformational changes that differ for each family of membrane transporters. Describing the protein conformational changes of different natures and magnitudes is therefore at the heart of mechanistic studies of membrane transporters. Here, using a number of membrane transporters from diverse families, we present common protocols used in setting up and performing molecular dynamics simulations of membrane transporters and in analyzing the results, in order to characterize relevant motions of the system. The emphasis will be on highlighting how optimal design of molecular dynamics simulations combined with mechanistically oriented analysis can shed light onto key functionally relevant protein conformational changes in this family of membrane proteins.

The power stroke driven by ATP binding in CFTR as studied by molecular dynamics simulations

Cystic fibrosis transmembrane conductance regulator (CFTR) is a chloride channel belonging to the ATP binding cassette (ABC) protein superfamily. Currently, it remains unclear how ATP binding causes the opening of the channel gate at the molecular level. To clarify this mechanism, we first constructed an atomic model of the inward-facing CFTR using the X-ray structures of other ABC proteins. Molecular dynamics (MD) simulations were then performed to explore the structure and dynamics of the inward-facing CFTR in a membrane environment. In the MgATP-bound state, two nucleotide-binding domains (NBDs) formed a head-to-tail type of dimer, in which the ATP molecules were sandwiched between the Walker A and signature motifs. Alternatively, one of the final MD structures in the apo state was similar to that of a "closed-apo" conformation found in the X-ray analysis of ATP-free MsbA. Principal component analysis for the MD trajectory indicated that NBD dimerization causes significant structural and dynamical changes in the transmembrane domains (TMDs), which is likely indicative of the formation of a chloride ion access path. This study suggests that the free energy gain from ATP binding acts as a driving force not only for NBD dimerization but also for NBD-TMD concerted motions.

Molecular dynamics simulations of membrane proteins

Membrane proteins control the traffic across cell membranes and thereby play an essential role in cell function from transport of various solutes to immune response via molecular recognition. Because it is very difficult to determine the structures of membrane proteins experimentally, computational methods have been increasingly used to study their structure and function. Here we focus on two classes of membrane proteins-ion channels and transporters-which are responsible for the generation of action potentials in nerves, muscles, and other excitable cells. We describe how computational methods have been used to construct models for these proteins and to study the transport mechanism. The main computational tool is the molecular dynamics (MD) simulation, which can be used for everything from refinement of protein structures to free energy calculations of transport processes. We illustrate with specific examples from gramicidin and potassium channels and aspartate transporters how the function of these membrane proteins can be investigated using MD simulations.

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.

Molecular-dynamics simulations of the ATP/apo state of a multidrug ATP-binding cassette transporter provide a structural and mechanistic basis for the asymmetric occluded state

ATP-binding cassette transporters use the energy of ATP hydrolysis to transport substrates across cellular membranes. They have two transmembrane domains and two cytosolic nucleotide-binding domains. Biochemical studies have characterized an occluded state of the transporter in which nucleotide is tenaciously bound in one active site, whereas the opposite active site is empty or binds nucleotide loosely. Here, we report molecular-dynamics simulations of the bacterial multidrug ATP-binding cassette transporter Sav1866. In two simulations of the ATP/apo state, the empty site opened substantially by way of rotation of the nucleotide-binding domain (NBD) core subdomain, whereas the ATP-bound site remained occluded and intact. We correlate our findings with elastic network and molecular-dynamics simulation analyses of the Sav1866 NBD monomer, and with existing experimental data, to argue that the observed transition is physiological, and that the final structure observed in the ATP/apo simulations corresponds to the tight/loose state of the NBD dimer characterized experimentally.

Conformational coupling of the nucleotide-binding and the transmembrane domains in ABC transporters

Basic architecture of ABC transporters includes two transmembrane domains (TMDs) and two nucleotide-binding domains (NBDs). Although the transport process takes place in the TMDs, which provide the substrate translocation pathway across the cell membrane and control its accessibility between the two sides of the membrane, the energy required for the process is provided by conformational changes induced in the NBDs by binding and hydrolysis of ATP. Nucleotide-dependent conformational changes in the NBDs, therefore, need to be coupled to structural changes in the TMDs. Using molecular dynamics simulations, we have investigated the structural elements involved in the conformational coupling between the NBDs and the TMDs in the Escherichia coli maltose transporter, an ABC importer for which an intact structure is available both in inward-facing and outward-facing conformations. The prevailing model of coupling is primarily based on a single structural motif, known as the coupling helices, as the main structural element for the NBD-TMD coupling. Surprisingly, we find that in the absence of the NBDs the coupling helices can be conformationally decoupled from the rest of the TMDs, despite their covalent connection. That is, the structural integrity of the coupling helices and their tight coupling to the core of the TMDs rely on the contacts provided by the NBDs. Based on the conformational and dynamical analysis of the simulation trajectories, we propose that the core coupling elements in the maltose transporter involve contributions from several structural motifs located at the NBD-TMD interface, namely, the EAA loops from the TMDs, and the Q-loop and the ENI motifs from the NBDs. These three structural motifs in small ABC importers show a high degree of correlation in motion and mediate the necessary conformational coupling between the core of TMDs and the helical subdomains of NBDs. A comprehensive analysis of the structurally known ABC transporters shows a high degree of conservation of the identified 3-motif coupling elements only in the subfamily of small ABC importers, suggesting a distinct mode of NBD-TMD coupling from the other two major ABC transporter folds, namely large ABC importers and ABC exporters.

Structural consequences of ATP hydrolysis on the ABC transporter NBD dimer: molecular dynamics studies of HlyB

ABC transporters are a large and important family of membrane proteins involved in substrate transport across the membrane. The transported substrates are quite diverse, ranging from monatomic ions to large biomolecules. Consequently, some ABC transporters are involved in biomedically relevant situations, from genetic diseases to multidrug resistance. The most conserved domains in ABC transporters are the nucleotide binding domains (NBDs), which form a dimer responsible for the binding and hydrolysis of ATP, concomitantly with substrate translocation. To elucidate how ATP hydrolysis structurally affects the NBD dimer, and consequently the transporter, we performed a molecular dynamics study on the NBD dimer of the HlyB ABC exporter. We have observed a change in the contact surface between the monomers after hydrolysis, even though we have not seen dimer opening in any of the five 100 ns simulations. We have also identified specific regions that respond to ATP hydrolysis, in particular the X-loop motif of ABC exporters, which has been shown to be in contact with the coupling helices of the transmembrane domains (TMDs). We propose that this motif is an important part of the NBD-TMD communication in ABC exporters. Through nonequilibrium analysis, we have also identified gradual conformational changes within a short time scale after ATP hydrolysis.

Inter-domain communication mechanisms in an ABC importer: a molecular dynamics study of the MalFGK2E complex

ATP-Binding Cassette transporters are ubiquitous membrane proteins that convert the energy from ATP-binding and hydrolysis into conformational changes of the transmembrane region to allow the translocation of substrates against their concentration gradient. Despite the large amount of structural and biochemical data available for this family, it is still not clear how the energy obtained from ATP hydrolysis in the ATPase domains is “transmitted” to the transmembrane domains. In this work, we focus our attention on the consequences of hydrolysis and inorganic phosphate exit in the maltose uptake system (MalFGK₂E) from Escherichia coli. The prime goal is to identify and map the structural changes occurring during an ATP-hydrolytic cycle. For that, we use extensive molecular dynamics simulations to study three potential intermediate states (with 10 replicates each): an ATP-bound, an ADP plus inorganic phosphate-bound and an ADP-bound state. Our results show that the residues presenting major rearrangements are located in the A-loop, in the helical sub-domain, and in the “EAA motif” (especially in the “coupling helices” region). Additionally, in one of the simulations with ADP we were able to observe the opening of the NBD dimer accompanied by the dissociation of ADP from the ABC signature motif, but not from its corresponding P-loop motif. This work, together with several other MD studies, suggests a common communication mechanism both for importers and exporters, in which ATP-hydrolysis induces conformational changes in the helical sub-domain region, in turn transferred to the transmembrane domains via the “coupling helices”.

ABC transporters are membrane proteins that couple ATP binding and hydrolysis with the active transport of substrates across membranes. These transporters form one of the largest families of membrane proteins and they can be found in all phyla of life. Moreover, some members of this family are involved in several genetic diseases (such as cystic fibrosis) and in multidrug resistance in bacteria, fungi and mammals. In this work, we use molecular dynamics simulations to study conformational changes due to ATP hydrolysis in an ABC transporter responsible for maltose uptake in E. coli. These conformational changes arising from one side of the protein (NBDs – Nucleotide Binding domains) where ATP binds, are propagated across the protein to more distant regions. Additionally, we can observe an NBD dimer interface dissociation event upon inorganic phosphate exit. These simulations together with other theoretical studies suggest that there is a general inter-domain communication mechanism common to importers and exporters.

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

ABC transporters play important roles in all types of organisms by participating in physiological and pathological processes. In order to modulate the function of ABC transporters, detailed knowledge regarding their structure and dynamics is necessary. Available structures of ABC proteins indicate three major conformations, a nucleotide-bound "bottom-closed" state with the two nucleotide binding domains (NBDs) tightly closed, and two nucleotide-free conformations, the "bottom-closed" and the "bottom-open", which differ in the extent of separation of the NBDs. However, it remains a question how the widely open conformation should be interpreted, and whether hydrolysis at one of the sites can drive conformational transitions while the NBDs remain in contact. To extend our knowledge, we have investigated the dynamic properties of the Sav1866 transporter using molecular dynamics (MD) simulations. We demonstrate that the replacement of one ATP by ADP alters the correlated motion patterns of the NBDs and the transmembrane domains (TMD). The results suggest that the hydrolysis of a single nucleotide could lead to extracellular closure, driving the transport cycle. Essential dynamics analysis of simulations suggests that single nucleotide hydrolysis can drive the system toward a "bottom-closed" apo conformation similar to that observed in the structure of the MsbA transporter. We also found significant structural instability of the "bottom-open" form of the transporters in simulations. Our results suggest that ATP hydrolysis at one of the sites promotes transport related conformational changes leading to the "bottom-closed" apo conformation, which could thus be physiologically more relevant for describing the structure of the apo state.

Mutant cycles at CFTR's non-canonical ATP-binding site support little interface separation during gating

Cystic fibrosis transmembrane conductance regulator (CFTR) is a chloride channel belonging to the adenosine triphosphate (ATP)-binding cassette (ABC) superfamily. ABC proteins share a common molecular mechanism that couples ATP binding and hydrolysis at two nucleotide-binding domains (NBDs) to diverse functions. This involves formation of NBD dimers, with ATP bound at two composite interfacial sites. In CFTR, intramolecular NBD dimerization is coupled to channel opening. Channel closing is triggered by hydrolysis of the ATP molecule bound at composite site 2. Site 1, which is non-canonical, binds nucleotide tightly but is not hydrolytic. Recently, based on kinetic arguments, it was suggested that this site remains closed for several gating cycles. To investigate movements at site 1 by an independent technique, we studied changes in thermodynamic coupling between pairs of residues on opposite sides of this site. The chosen targets are likely to interact based on both phylogenetic analysis and closeness on structural models. First, we mutated T460 in NBD1 and L1353 in NBD2 (the corresponding site-2 residues become energetically coupled as channels open). Mutation T460S accelerated closure in hydrolytic conditions and in the nonhydrolytic K1250R background; mutation L1353M did not affect these rates. Analysis of the double mutant showed additive effects of mutations, suggesting that energetic coupling between the two residues remains unchanged during the gating cycle. We next investigated pairs 460–1348 and 460–1375. Although both mutations H1348A and H1375A produced dramatic changes in hydrolytic and nonhydrolytic channel closing rates, in the corresponding double mutants these changes proved mostly additive with those caused by mutation T460S, suggesting little change in energetic coupling between either positions 460–1348 or positions 460–1375 during gating. These results provide independent support for a gating model in which ATP-bound composite site 1 remains closed throughout the gating cycle.

Conformational changes induced by ATP-hydrolysis in an ABC transporter: a molecular dynamics study of the Sav1866 exporter

ATP-Binding Cassette (ABC) transporters are ubiquitous membrane proteins that use energy from ATP binding or/and hydrolysis to actively transport allocrites across membranes. In this study, we identify ATP-hydrolysis induced conformational changes in a complete ABC exporter (Sav1866) from Staphylococcus aureaus, using molecular dynamics (MD) simulations. By performing MD simulations on the ATP and ADP+IP bound states, we identify the conformational consequences of hydrolysis, showing that the major rearrangements are not restricted to the NBDs, but extend to the transmembrane domains (TMDs) external regions. For the first time, to our knowledge, we see, within the context of a complete transporter, NBD dimer opening in the ADP+IP state in contrast with all ATP-bound states. This opening results from the dissociation of the ABC signature motif from the nucleotide. In addition, in both states, we observe the opening of a gate entrance in the intracellular loop region leading to the exposure of the TMDs internal cavity to the cytoplasm. To see if this opening was large enough to allow allocrite transport, the adiabatic energy profile for doxorubicin passage was determined. For both states, this profile, although an approximation, is overall downhill from the cytoplasmatic to the extracellular side, and the local energy barriers along the TMDs are relatively small, evidencing the exporter nature of Sav1866. The major difference between states is an energy barrier located in the cytoplasmic gate region, which becomes reduced upon hydrolysis, suggesting that allocrite passage is facilitated, and evidencing a possible molecular mechanism for the active transport in these proteins.

Capturing Functional Motions of Membrane Channels and Transporters with Molecular Dynamics Simulation

Conformational changes of proteins are involved in all aspects of protein function in biology. Almost all classes of proteins respond to changes in their environment, ligand binding, and interaction with other proteins and regulatory agents through undergoing conformational changes of various degrees and magnitudes. Membrane channels and transporters are the major classes of proteins that are responsible for mediating efficient and selective transport of materials across the cellular membrane. Similar to other proteins, they take advantage of conformational changes to make transitions between various functional states. In channels, large-scale conformational changes are mostly involved in the process of "gating", i.e., opening and closing of the pore of the channel protein in response to various signals. In transporters, conformational changes constitute various steps of the conduction process, and, thus, are more closely integrated in the transport process. Owing to significant progress in developing highly efficient parallel algorithms in molecular dynamics simulations and increased computational resources, and combined with the availability of high-resolution, atomic structures of membrane proteins, we are in an unprecedented position to use computer simulation and modeling methodologies to investigate the mechanism of function of membrane channels and transporters. While the entire transport cycle is still out of reach of current methodologies, many steps involved in the function of transport proteins have been characterized with molecular dynamics simulations. Here, we present several examples of such studies from our laboratory, in which functionally relevant conformational changes of membrane channels and transporters have been characterized using extended simulations.

Dynamics and structural changes induced by ATP binding in SAV1866, a bacterial ABC exporter

Multidrug transporters of the ATP-binding cassette family export a wide variety of compounds across membranes in both prokaryotes and eukaryotes, using ATP hydrolysis as energy source. Several of these membrane proteins are of clinical importance. Although biochemical and structural studies have provided insights into the mechanism underlying substrate transport, many key questions subsist regarding the molecular and structural nature of this mechanism. In particular, the detailed conformational changes occurring during the catalytic cycle are still elusive. We explored the conformational changes occurring upon ATP/Mg(2+) binding using molecular dynamics simulations starting from the nucleotide-bound structure of SAV1866 embedded in an explicit lipid bilayer. The removal of nucleotide revealed a major rearrangement in the outer membrane leaflet portion of the transmembrane domain (TMD) resulting in the closure of the central cavity at the extracellular side. This closure is similar to that observed in the crystal nucleotide-free structures. The interface of the nucleotide-binding domain dimer (NDB) is significantly more hydrated in the nucleotide-free trajectory though it is not disrupted. This finding suggests that the TMD closure could occur as a first step preceding the dissociation of the dimer. The transmission pathway of the signal triggered by the removal of ATP/Mg(2+) mainly involves the conserved Q-loop and X-loop as well as TM6.

Dynamics of alpha-helical subdomain rotation in the intact maltose ATP-binding cassette transporter

ATP-binding cassette (ABC) transporters are powered by a nucleotide-binding domain dimer that opens and closes during cycles of ATP hydrolysis. These domains consist of a RecA-like subdomain and an α-helical subdomain that is specific to the family. Many studies on isolated domains suggest that the helical subdomain rotates toward the RecA-like subdomain in response to ATP binding, moving the family signature motif into a favorable position to interact with the nucleotide across the dimer interface. Moreover, the transmembrane domains are docked into a cleft at the interface between these subdomains, suggesting a putative role of the rotation in interdomain communication. Electron paramagnetic resonance spectroscopy was used to study the dynamics of this rotation in the intact Escherichia coli maltose transporter MalFGK(2). This importer requires a periplasmic maltose-binding protein (MBP) that activates ATP hydrolysis by promoting the closure of the cassette dimer (MalK(2)). Whereas this rotation occurred during the transport cycle, it required not only trinucleotide, but also MBP, suggesting it is part of a global conformational change in the transporter. Interaction of AMP-PNP-Mg(2+) and a MBP that is locked in a closed conformation induced a transition from open MalK(2) to semiopen MalK(2) without significant subdomain rotation. Inward rotation of the helical subdomain and complete closure of MalK(2) therefore appear to be coupled to the reorientation of transmembrane helices and the opening of MBP, events that promote transfer of maltose into the transporter. After ATP hydrolysis, the helical subdomain rotates out as MalK(2) opens, resetting the transporter in an inward-facing conformation.