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

Signal Propagation in the ATPase Domain of Mycobacterium tuberculosis DNA Gyrase from Dynamical-Nonequilibrium Molecular Dynamics Simulations

DNA gyrases catalyze negative supercoiling of DNA, are essential for bacterial DNA replication, transcription, and recombination, and are important antibacterial targets in multiple pathogens, including Mycobacterium tuberculosis, which in 2021 caused >1.5 million deaths worldwide. DNA gyrase is a tetrameric (A2B2) protein formed from two subunit types: gyrase A (GyrA) carries the breakage-reunion active site, whereas gyrase B (GyrB) catalyzes ATP hydrolysis required for energy transduction and DNA translocation. The GyrB ATPase domains dimerize in the presence of ATP to trap the translocated DNA (T-DNA) segment as a first step in strand passage, for which hydrolysis of one of the two ATPs and release of the resulting inorganic phosphate is rate-limiting. Here, dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations of the dimeric 43 kDa N-terminal fragment of M. tuberculosis GyrB show how events at the ATPase site (dissociation/hydrolysis of bound nucleotides) are propagated through communication pathways to other functionally important regions of the GyrB ATPase domain. Specifically, our simulations identify two distinct pathways that respectively connect the GyrB ATPase site to the corynebacteria-specific C-loop, thought to interact with GyrA prior to DNA capture, and to the C-terminus of the GyrB transduction domain, which in turn contacts the C-terminal GyrB topoisomerase-primase (TOPRIM) domain responsible for interactions with GyrA and the centrally bound G-segment DNA. The connection between the ATPase site and the C-loop of dimeric GyrB is consistent with the unusual properties of M. tuberculosis DNA gyrase relative to those from other bacterial species.

Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

Solid-State NMR Reveals Asymmetric ATP Hydrolysis in the Multidrug ABC Transporter BmrA

The detailed mechanism of ATP hydrolysis in ATP-binding cassette (ABC) transporters is still not fully understood. Here, we employed ³¹P solid-state NMR to probe the conformational changes and dynamics during the catalytic cycle by locking the multidrug ABC transporter BmrA in prehydrolytic, transition, and posthydrolytic states, using a combination of mutants and ATP analogues. The ³¹P spectra reveal that ATP binds strongly in the prehydrolytic state to both ATP-binding sites as inferred from the analysis of the nonhydrolytic E504A mutant. In the transition state of wild-type BmrA, the symmetry of the dimer is broken and only a single site is tightly bound to ADP:Mg²⁺:vanadate, while the second site is more ‘open’ allowing exchange with the nucleotides in the solvent. In the posthydrolytic state, weak binding, as characterized by chemical exchange with free ADP and by asymmetric ³¹P–³¹P two-dimensional (2D) correlation spectra, is observed for both sites. Revisiting the ¹³C spectra in light of these findings confirms the conformational nonequivalence of the two nucleotide-binding sites in the transition state. Our results show that following ATP binding, the symmetry of the ATP-binding sites of BmrA is lost in the ATP-hydrolysis step, but is then recovered in the posthydrolytic ADP-bound state.

The fatty acid site is coupled to functional motifs in the SARS-CoV-2 spike protein and modulates spike allosteric behaviour

The SARS-CoV-2 spike protein is the first contact point between the SARS-CoV-2 virus and host cells and mediates membrane fusion. Recently, a fatty acid binding site was identified in the spike (Toelzer et al. Science 2020). The presence of linoleic acid at this site modulates binding of the spike to the human ACE2 receptor, stabilizing a locked conformation of the protein. Here, dynamical-nonequilibrium molecular dynamics simulations reveal that this fatty acid site is coupled to functionally relevant regions of the spike, some of them far from the fatty acid binding pocket. Removal of a ligand from the fatty acid binding site significantly affects the dynamics of distant, functionally important regions of the spike, including the receptor-binding motif, furin cleavage site and fusion-peptide-adjacent regions. Simulations of the D614G mutant show differences in behaviour between these clinical variants of the spike: the D614G mutant shows a significantly different conformational response for some structural motifs relevant for binding and fusion. The simulations identify structural networks through which changes at the fatty acid binding site are transmitted within the protein. These communication networks significantly involve positions that are prone to mutation, indicating that observed genetic variation in the spike may alter its response to linoleate binding and associated allosteric communication.

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.

Allosteric communication in class A β-lactamases occurs via cooperative coupling of loop dynamics

Understanding allostery in enzymes and tools to identify it offer promising alternative strategies to inhibitor development. Through a combination of equilibrium and nonequilibrium molecular dynamics simulations, we identify allosteric effects and communication pathways in two prototypical class A β-lactamases, TEM-1 and KPC-2, which are important determinants of antibiotic resistance. The nonequilibrium simulations reveal pathways of communication operating over distances of 30 Å or more. Propagation of the signal occurs through cooperative coupling of loop dynamics. Notably, 50% or more of clinically relevant amino acid substitutions map onto the identified signal transduction pathways. This suggests that clinically important variation may affect, or be driven by, differences in allosteric behavior, providing a mechanism by which amino acid substitutions may affect the relationship between spectrum of activity, catalytic turnover, and potential allosteric behavior in this clinically important enzyme family. Simulations of the type presented here will help in identifying and analyzing such differences.

Dynamical nonequilibrium molecular dynamics reveals the structural basis for allostery and signal propagation in biomolecular systems

ABSTRACT

A dynamical approach to nonequilibrium molecular dynamics (D-NEMD), proposed in the 1970s by Ciccotti et al., is undergoing a renaissance and is having increasing impact in the study of biological macromolecules. This D-NEMD approach, combining MD simulations in stationary (in particular, equilibrium) and nonequilibrium conditions, allows for the determination of the time-dependent structural response of a system using the Kubo-Onsager relation. Besides providing a detailed picture of the system's dynamic structural response to an external perturbation, this approach also has the advantage that the statistical significance of the response can be assessed. The D-NEMD approach has been used recently to identify a general mechanism of inter-domain signal propagation in nicotinic acetylcholine receptors, and allosteric effects in β -lactamase enzymes, for example. It complements equilibrium MD and is a very promising approach to identifying and analysing allosteric effects. Here, we review the D-NEMD approach and its application to biomolecular systems, including transporters, receptors, and enzymes.

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.

Characterization of the nucleotide-binding domain NsrF from the BceAB-type ABC-transporter NsrFP from the human pathogen Streptococcus agalactiae

Treatment of bacterial infections is a great challenge of our era due to the various resistance mechanisms against antibiotics. Antimicrobial peptides are considered to be potential novel compound as antibiotic treatment. However, some bacteria, especially many human pathogens, are inherently resistant to these compounds, due to the expression of BceAB-type ABC transporters. This rather new transporter family is not very well studied. Here, we report the first full characterization of the nucleotide binding domain of a BceAB type transporter from Streptococcus agalactiae, namely SaNsrF of the transporter SaNsrFP, which confers resistance against nisin and gallidermin. We determined the NTP hydrolysis kinetics and used molecular modeling and simulations in combination with small angle X-ray scattering to obtain structural models of the SaNsrF monomer and dimer. The fact that the SaNsrF_H202A variant displayed no ATPase activity was rationalized in terms of changes of the structural dynamics of the dimeric interface. Kinetic data show a clear preference for ATP as a substrate, and the prediction of binding modes allowed us to explain this selectivity over other NTPs.

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.

Identification of the Initial Steps in Signal Transduction in the α4β2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations

Nicotinic acetylcholine receptors (nAChRs) modulate synaptic transmission in the nervous system. These receptors have emerged as therapeutic targets in drug discovery for treating several conditions, including Alzheimer's disease, pain, and nicotine addiction. In this in silico study, we use a combination of equilibrium and nonequilibrium molecular dynamics simulations to map dynamic and structural changes induced by nicotine in the human α4β2 nAChR. They reveal a striking pattern of communication between the extracellular binding pockets and the transmembrane domains (TMDs) and show the sequence of conformational changes associated with the initial steps in this process. We propose a general mechanism for signal transduction for Cys-loop receptors: the mechanistic steps for communication proceed firstly through loop C in the principal subunit, and are subsequently transmitted, gradually and cumulatively, to loop F of the complementary subunit, and then to the TMDs through the M2-M3 linker.

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.

Unexplored Nucleotide Binding Modes for the ABC Exporter MsbA

The ATP-binding cassette (ABC) transporter MsbA is an ATP-driven lipid-A flippase. It belongs to the ABC protein superfamily whose members are characterized by conserved motifs in their nucleotide binding domains (NBDs), which are responsible for ATP hydrolysis. Recently, it was found that MsbA could catalyze a reverse adenylate kinase (rAK)-like reaction in addition to ATP hydrolysis. Both reactions are connected and mediated by the same conserved NBD domains. Here, the structural foundations underlying the nucleotide binding to MsbA were therefore explored using a concerted approach based on conventional- and DNP-enhanced solid-state NMR, pulsed-EPR, and MD simulations. MsbA reconstituted into lipid bilayers was trapped in various catalytic states corresponding to intermediates of the coupled ATPase-rAK mechanism. The analysis of nucleotide-binding dependent chemical shift changes, and the detection of through-space contacts between bound nucleotides and MsbA within these states provides evidence for an additional nucleotide-binding site in close proximity to the Q-loop and the His-Switch. By replacing Mg²⁺ with Mn²⁺ and employing pulsed EPR spectroscopy, evidence is provided that this newly found nucleotide binding site does not interfere with the coordination of the required metal ion. Molecular dynamic (MD) simulations of nucleotide and metal binding required for the coupled ATPase-rAK mechanism have been used to corroborate these experimental findings and provide additional insight into nucleotide location, orientation, and possible binding modes.

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.

How Intrinsic Dynamics Mediates the Allosteric Mechanism in the ABC Transporter Nucleotide Binding Domain Dimer

A protein's architecture facilitates specific motions-intrinsic dynamic modes-that are employed to effect function. Here we used molecular dynamics (MD) simulations to investigate the dynamics of the MJ0796 ABC transporter nucleotide-binding domain (NBD). ABC transporter NBDs form a rotationally symmetric dimer whereby two equivalent active sites are formed at their interface; in complex with a dimer of transmembrane domains they hydrolyze ATP to energize translocation of substrates across cellular membranes. Our data suggest the ABC NBD's ensemble of functional states can be understood predominately in terms of conformational changes between its major subdomains, occurring along two orthogonal dynamic modes. The data show that ligands and oligomeric interactions modulate the equilibrium conformation of the NBD with respect to these motions, suggesting that allostery is achieved by affecting the energetic profile along these two modes. The observed dynamics and allostery integrate consonantly and logically within a mechanistic framework for the ABC NBD dimer, which is supported by a large body of experimental and theoretical data, providing a higher resolution view of the enzyme's dynamic cycle. Our study shows how valuable mechanistic inferences can be derived from accessible short-time scale MD simulations of an enzyme's substructures.

Large-scale classification of P-glycoprotein inhibitors using SMILES-based descriptors

P-glycoprotein (Pgp) inhibition has been considered as an effective strategy towards combating multidrug-resistant cancers. Owing to the substrate promiscuity of Pgp, the classification of its interacting ligands is not an easy task and is an ongoing issue of debate. Chemical structures can be represented by the simplified molecular input line entry system (SMILES) in the form of linear string of symbols. In this study, the SMILES notations of 2254 Pgp inhibitors including 1341 active, and 913 inactive compounds were used for the construction of a SMILE-based classification model using CORrelation And Logic (CORAL) software. The model provided an acceptable predictive performance as observed from statistical parameters consisting of accuracy, sensitivity and specificity that afforded values greater than 70% and MCC value greater than 0.6 for training, calibration and validation sets. In addition, the CORAL method highlighted chemical features that may contribute to increased and decreased Pgp inhibitory activities. This study highlights the potential of CORAL software for rapid screening of prospective compounds from a large chemical space and provides information that could aid in the design and development of potential Pgp inhibitors.

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.

Conformational Motions and Functionally Key Residues for Vitamin B12 Transporter BtuCD-BtuF Revealed by Elastic Network Model with a Function-Related Internal Coordinate

BtuCD-BtuF from Escherichia coli is a binding protein-dependent adenosine triphosphate (ATP)-binding cassette (ABC) transporter system that uses the energy of ATP hydrolysis to transmit vitamin B12 across cellular membranes. Experimental studies have showed that during the transport cycle, the transporter undergoes conformational transitions between the "inward-facing" and "outward-facing" states, which results in the open-closed motions of the cytoplasmic gate of the transport channel. The opening-closing of the channel gate play critical roles for the function of the transporter, which enables the substrate vitamin B12 to be translocated into the cell. In the present work, the extent of opening of the cytoplasmic gate was chosen as a function-related internal coordinate. Then the mean-square fluctuation of the internal coordinate, as well as the cross-correlation between the displacement of the internal coordinate and the movement of each residue in the protein, were calculated based on the normal mode analysis of the elastic network model to analyze the function-related motions encoded in the structure of the system. In addition, the key residues important for the functional motions of the transporter were predicted by using a perturbation method. In order to facilitate the calculations, the internal coordinate was introduced as one of the axes of the coordinate space and the conventional Cartesian coordinate space was transformed into the internal/Cartesian space with linear approximation. All the calculations were carried out in this internal/Cartesian space. Our method can successfully identify the functional motions and key residues for the transporter BtuCD-BtuF, which are well consistent with the experimental observations.

Do adsorbed drugs onto P-glycoprotein influence its efflux capability?

The membrane biophysical aspects by which multidrug resistance (MDR) relate to the ABC transporter function still remain largely unknown. Notwithstanding the central role that efflux pumps like P-glycoprotein have in MDR onset, experimental studies classified additionally the lipid micro-environment where P-gp is inserted as a determinant for the increased efflux capability demonstrated in MDR cell lines. Recently, a nonlinear model for drug-membrane interactions showed that, upon drug adsorption, long-range mechanical alterations are predicted to affect the P-gp ATPase function at external drug concentrations of ∼10-100 μM. However, our results also show that drug adsorption may also occur at P-gp nucleotide-binding domains where conformational changes drive the efflux cycle. Thus, we assessed the effect of drug adsorption to both protein-water and lipid-water interfaces by means of molecular dynamics simulations. The results show that free energies of adsorption are lower for modulators in both lipid/water and protein/water interfaces. Important differences in drug-protein interactions, protein dynamics and membrane biophysical characteristics were observed between the different classes. Therefore, we hypothesize that drug adsorption to the protein and lipid-water interface accounts for a complex network of events that affect the ability of transporters to efflux drugs.

Cystic Fibrosis Transmembrane Conductance Regulator (CFTR): CLOSED AND OPEN STATE CHANNEL MODELS

The cystic fibrosis transmembrane conductance regulator (CFTR) is a member of the ATP-binding cassette (ABC) transporter superfamily. CFTR controls the flow of anions through the apical membrane of epithelia. Dysfunctional CFTR causes the common lethal genetic disease cystic fibrosis. Transitions between open and closed states of CFTR are regulated by ATP binding and hydrolysis on the cytosolic nucleotide binding domains, which are coupled with the transmembrane (TM) domains forming the pathway for anion permeation. Lack of structural data hampers a global understanding of CFTR and thus the development of "rational" approaches directly targeting defective CFTR. In this work, we explored possible conformational states of the CFTR gating cycle by means of homology modeling. As templates, we used structures of homologous ABC transporters, namely TM(287-288), ABC-B10, McjD, and Sav1866. In the light of published experimental results, structural analysis of the transmembrane cavity suggests that the TM(287-288)-based CFTR model could correspond to a commonly occupied closed state, whereas the McjD-based model could represent an open state. The models capture the important role played by Phe-337 as a filter/gating residue and provide structural information on the conformational transition from closed to open channel.

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.

Strategies for improving potassium use efficiency in plants

Potassium is a macronutrient that is crucial for healthy plant growth. Potassium availability, however, is often limited in agricultural fields and thus crop yields and quality are reduced. Therefore, improving the efficiency of potassium uptake and transport, as well as its utilization, in plants is important for agricultural sustainability. This review summarizes the current knowledge on the molecular mechanisms involved in potassium uptake and transport in plants, and the molecular response of plants to different levels of potassium availability. Based on this information, four strategies for improving potassium use efficiency in plants are proposed; 1) increased root volume, 2) increasing efficiency of potassium uptake from the soil and translocation in planta, 3) increasing mobility of potassium in soil, and 4) molecular breeding new varieties with greater potassium efficiency through marker assisted selection which will require identification and utilization of potassium associated quantitative trait loci.

Allosteric coupling between the intracellular coupling helix 4 and regulatory sites of the first nucleotide-binding domain of CFTR

Cystic fibrosis is caused by mutations in CFTR (cystic fibrosis transmembrane conductance regulator), leading to folding and processing defects and to chloride channel gating misfunction. CFTR is regulated by ATP binding to its cytoplasmic nucleotide-binding domains, NBD1 and NBD2, and by phosphorylation of the NBD1 regulatory insert (RI) and the regulatory extension (RE)/R region. These regulatory effects are transmitted to the rest of the channel via NBD interactions with intracellular domain coupling helices (CL), particularly CL4. Using a sensitive method for detecting inter-residue correlations between chemical shift changes in NMR spectra, an allosteric network was revealed within NBD1, with a construct lacking RI. The CL4-binding site couples to the RI-deletion site and the C-terminal residues of NBD1 that precede the R region in full-length CFTR. Titration of CL4 peptide into NBD1 perturbs the conformational ensemble in these sites with similar titration patterns observed in F508del, the major CF-causing mutant, and in suppressor mutants F494N, V510D and Q637R NBD1, as well as in a CL4-NBD1 fusion construct. Reciprocally, the C-terminal mutation, Q637R, perturbs dynamics in these three sites. This allosteric network suggests a mechanism synthesizing diverse regulatory NBD1 interactions and provides biophysical evidence for the allosteric coupling required for CFTR function.

Putative dioxygen-binding sites and recognition of tigecycline and minocycline in the tetracycline-degrading monooxygenase TetX

Expression of the aromatic hydroxylase TetX under aerobic conditions confers bacterial resistance against tetracycline antibiotics. Hydroxylation inactivates and degrades tetracyclines, preventing inhibition of the prokaryotic ribosome. X-ray crystal structure analyses of TetX in complex with the second-generation and third-generation tetracyclines minocycline and tigecycline at 2.18 and 2.30 Å resolution, respectively, explain why both clinically potent antibiotics are suitable substrates. Both tetracyclines bind in a large tunnel-shaped active site in close contact to the cofactor FAD, pre-oriented for regioselective hydroxylation to 11a-hydroxytetracyclines. The characteristic bulky 9-tert-butylglycylamido substituent of tigecycline is solvent-exposed and does not interfere with TetX binding. In the TetX-minocycline complex a second binding site for a minocycline dimer is observed close to the active-site entrance. The pocket is formed by the crystal packing arrangement on the surface of two neighbouring TetX monomers. Crystal structure analysis at 2.73 Å resolution of xenon-pressurized TetX identified two adjacent Xe-binding sites. These putative dioxygen-binding cavities are located in the substrate-binding domain next to the active site. Molecular-dynamics simulations were performed in order to characterize dioxygen-diffusion pathways to FADH2 at the active site.

Substrate versus inhibitor dynamics of P-glycoprotein

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

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.

A conserved interdomain communication pathway of pseudosymmetrically distributed residues affects substrate specificity of the fungal multidrug transporter Cdr1p

Understanding the communication pathways between remote sites in proteins is of key importance for understanding their function and mechanism of action. These remain largely unexplored among the pleiotropic drug resistance (PDR) representatives of the ubiquitous superfamily of ATP-binding cassette (ABC) transporters. To identify functionally coupled residues important for the polyspecific transport by the fungal ABC multidrug transporter Cdr1p a new selection strategy, towards increased resistance to a preferred substrate of the homologous Snq2p, was applied to a library of randomly generated mutants. The single amino acid substitutions, located pseudosymmetrically in each domain of the internally duplicated protein: the H-loop of the N-terminal nucleotide binding domain (NBD1) (C363R) and in the C-terminal NBD2 region preceding Walker A (V885G). The central regions of the first transmembrane helices 1 and 7 of both transmembrane domains were also affected by the G521S/D and A1208V substitutions respectively. Although the mutants were expressed at a similar level and located correctly to the plasma membrane, they selectively affected transport of multiple drugs, including azole antifungals. The synergistic effects of combined mutations on drug resistance, drug dependent ATPase activity and transport support the view inferred from the statistical coupling analysis (SCA) of aminoacid coevolution and mutational analysis of other ABC transporter families that these residues are an important part of the conserved, allosterically coupled interdomain communication network. Our results shed new light on the communication between the pseudosymmetrically arranged domains in a fungal PDR ABC transporter and reveal its profound influence on substrate specificity.

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.

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.