Comparison of different schemes to treat long-range electrostatic interactions in molecular dynamics simulations of a protein crystal

Eight molecular dynamics simulations of a ubiquitin crystal unit cell were performed to investigate the effect of different schemes to treat the long-range electrostatic interactions as well as the need to include counter ions. A crystal system was chosen as the test system, because the higher charge density compared with a protein in solution makes it more sensitive to the way of treating the electrostatic interactions. Three different schemes of treating the long-range interactions were compared: straight cutoff, reaction-field approximation, and a lattice-sum method (P3M). For each of these schemes, two simulations were performed, one with and one without the counter ions. Two additional simulations with a reaction-field force and different initial placements of the counter ions were performed to examine the effect of the initial positions of the ions. The inclusion of long-range electrostatic interactions using either a reaction-field or a lattice-sum method proved to be necessary for the simulation of crystals. These two schemes did not differ much in their ability to reproduce the crystallographic structure. The inclusion of counter ions, on the other hand, seems not necessary for obtaining a stable simulation. The initial positions of the ions have a visible but small effect on the simulation.     

The use of immobilized ubiquitin for biosensor analysis of the mitochondrial subinteractome

Protein ubiquitination is an important mechanism responsible not only for specific labeling of proteins for their subsequent degradation; it also determines localization of proteins in the cell and regulation of protein-protein interactions. In the context of protein-protein interactions binding of (mono/poly)ubiquitinated molecules to proteins containing specific ubiquitin binding domains plays the decisive role. Formation of the ubiquitin interactome has been demonstrated for cytosol. Involvement of mitochondria and associated extramitochondrial proteins into such interactions still requires detailed investigation. In this study using an optical biosensor we have demonstrated binding of proteins of mouse brain mitochondrial lysates to immobilized monomeric ubiquitin. Model purified proteins, which are known to be associated with the outer mitochondrial compartment (glyceraldehyde-3-phosphate dehydorgenase, creatine phosphokinase), interacted with immobilized ubiquitin as well as with each other. This suggests that (poly)ubiquitinated chains may be involved in protein-protein interactions between ubiquitinated and non-ubiquitinated proteins and thus may contribute to formation of (mitochondrial) ubiquitin subinteractome.     

Data handling for interactive metabolomics: tools for studying the dynamics of metabolome-macromolecule interactions

All published metabolomics studies investigate changes in either absolute or relative quantities of metabolites. However, blood plasma, one of the most commonly studied biofluids for metabolomics applications, is a complex, heterogeneous mixture of lipoproteins, proteins, small organic molecules and ions which together undergo a variety of possible molecular interactions including metal complexation, chemical exchange processes, micellular compartmentation of metabolites, enzyme-mediated biotransformations and small-molecule-macromolecule binding. In particular, many low molecular weight (MW) compounds (including drugs) can exist both ‘free’ in solution and bound to proteins or within organised aggregates of macromolecules. To study the effects of e.g. disease on these interactions we suggest that new approaches are needed. We have developed a technique termed ‘interactive metabolomics’ or i-metabolomics. i-metabolomics can be defined as: “The study of interactions between low MW biochemicals and macromolecules in heterogeneous biosamples such as blood plasma, without pre-selection of the components of interest”. Standard 1D NMR experiments commonly used in metabolomics allow metabolite concentration differences between samples to be investigated because the intensity of each peak depends on the concentration of the compound in question. On the other hand, the instrument can be set-up to measure molecular interactions by monitoring the diffusion coefficients of molecules. According to the Stokes–Einstein equation, the diffusion coefficient of a molecule is inversely proportional to its effective size, as represented by the hydrodynamic radius. Therefore, when low MW compounds are non-covalently bound to proteins, the observed diffusion coefficient for the compound will be intermediate between those of its free and bound forms. By measuring diffusion by NMR, the degree of protein binding can be estimated for either low MW endogenous biochemicals or xenobiotics. This type of experiment is referred to as either Diffusion-Ordered Spectroscopy (DOSY) or Diffusion-Edited Spectroscopy, depending on the type of post-acquisition data processing applied to the spectra. Results presented in this paper demonstrate approaches for the non-selective modelling of metabolite-macromolecule interactions (i-metabolomics), whilst additionally highlighting some of the all too frequently ignored issues associated with interpretation of data derived from profiling of blood plasma.     

Data handling for interactive metabolomics: tools for studying the dynamics of metabolome-macromolecule interactions

All published metabolomics studies investigate changes in either absolute or relative quantities of metabolites. However, blood plasma, one of the most commonly studied biofluids for metabolomics applications, is a complex, heterogeneous mixture of lipoproteins, proteins, small organic molecules and ions which together undergo a variety of possible molecular interactions including metal complexation, chemical exchange processes, micellular compartmentation of metabolites, enzyme-mediated biotransformations and small-molecule-macromolecule binding. In particular, many low molecular weight (MW) compounds (including drugs) can exist both ‘free’ in solution and bound to proteins or within organised aggregates of macromolecules. To study the effects of e.g. disease on these interactions we suggest that new approaches are needed. We have developed a technique termed ‘interactive metabolomics’ or i-metabolomics. i-metabolomics can be defined as: “The study of interactions between low MW biochemicals and macromolecules in heterogeneous biosamples such as blood plasma, without pre-selection of the components of interest”. Standard 1D NMR experiments commonly used in metabolomics allow metabolite concentration differences between samples to be investigated because the intensity of each peak depends on the concentration of the compound in question. On the other hand, the instrument can be set-up to measure molecular interactions by monitoring the diffusion coefficients of molecules. According to the Stokes–Einstein equation, the diffusion coefficient of a molecule is inversely proportional to its effective size, as represented by the hydrodynamic radius. Therefore, when low MW compounds are non-covalently bound to proteins, the observed diffusion coefficient for the compound will be intermediate between those of its free and bound forms. By measuring diffusion by NMR, the degree of protein binding can be estimated for either low MW endogenous biochemicals or xenobiotics. This type of experiment is referred to as either Diffusion-Ordered Spectroscopy (DOSY) or Diffusion-Edited Spectroscopy, depending on the type of post-acquisition data processing applied to the spectra. Results presented in this paper demonstrate approaches for the non-selective modelling of metabolite-macromolecule interactions (i-metabolomics), whilst additionally highlighting some of the all too frequently ignored issues associated with interpretation of data derived from profiling of blood plasma.

     

High rates of calcium-free diffusion in the cytosol of living cells

Calcium (Ca²⁺) is a universal second messenger that participates in the regulation of innumerous physiological processes. The way in which local elevations of the cytosolic Ca²⁺ concentration spread in space and time is key for the versatility of the signals. Ca²⁺ diffusion in the cytosol is hindered by its interaction with proteins that act as buffers. Depending on the concentrations and the kinetics of the interactions, there is a large range of values at which Ca²⁺ diffusion can proceed. Having reliable estimates of this range, particularly of its highest end, which corresponds to the ions free diffusion, is key to understand how the signals propagate. In this work, we present the first experimental results with which the Ca²⁺-free diffusion coefficient is directly quantified in the cytosol of living cells. By means of fluorescence correlation spectroscopy experiments performed in Xenopus laevis oocytes and in cells of Saccharomyces cerevisiae, we show that the ions can freely diffuse in the cytosol at a higher rate than previously thought.     

On the dielectrically observable consequences of the diffusional motions of lipids and proteins in membranes

A system consisting of any array of cylindrical, polytopic membrane proteins (or protein complexes) possessed of a permanent dipole moment and immersed in a closed, spherical phospholipid bilayer sheet is considered. It is assumed that rotation of the protein (complex) in a plane normal to the membrane, if occurring, is restricted by viscous drag alone. Lateral diffusion is assumed either to be free and random or to be partially constrained by barriers of an unspecified nature. The dielectric relaxation times calculated for membrane protein rotation in a suspension of vesicles of the above type are much longer than those observed with globular proteins in aqueous solution, and fall in the mid-to-high audio frequency range. If the long range lateral diffusion of (charged) membrane protein complexes is essentially unrestricted, as in the "fluid mosaic" membrane model, dielectric relaxation times for lateral motions will lie, except in the case of the very smallest vesicles, in the sub-audio (ELF) range. If, in contrast, the lateral diffusion of membrane protein complexes is partially restricted by "barriers" or "long-range" interactions (of unspecified nature), significant dielectric dispersions may be expected in both audio- and radio-frequency ranges, the critical (characteristic) frequencies depending upon the average distance moved before a barrier is encountered. Similar analyses are given for rotational and translational motions of phospholipids. At very low frequencies, a dispersion due to vesicle orientation might in principle also be observed; the dielectrically observable extent of this rotation will depend, inter alia, upon the charge mobility and disposition of the membrane protein complexes, as well as, of course, on the viscosity of the aqueous phase. The role of electroosmotic interactions between double layer ions (and water dipoles) and proteins raised above the membrane surface is considered. In some cases, it seems likely that such interactions serve to raise the dielectric increment, relative to that which might otherwise have been expected, of dispersions due to protein motions in membranes. Depending upon the tortuosity of the ion-relaxation pathways, such a relaxation mechanism might lead to almost any characteristic frequency, and, even in the absence of protein/lipid motions, would cause dielectric spectra to be much broader than one might expect from a simple, macroscopic treatment.     

An enquiry into metabolite domains

It is currently assumed that two or more pools of the same metabolite can coexist in the cytosolic compartment of mammalian cells. These pools are thought to be generated by the differential subcellular location of enzymes and transporters, much in the way calcium microdomains arise by the combined workings of channels, buffers, and pumps. With the aim of estimating the amplitude and spatial dimensions of these metabolite pools, we developed an analytical tool based on Brownian diffusion and the turnover numbers of the proteins involved. The outcome of the analysis is that ATP, glucose, pyruvate, lactate, and glutamate cannot be concentrated at their sources to an extent that would affect their downstream targets. For these metabolites, and others produced by slow enzymes or transporters and present at micromolar concentrations or higher, the cytosol behaves as a well-mixed, homogenous compartment. In contrast, the analysis showed microdomains known to be generated by calcium channels and revealed that calcium and pH nanodomains are to be found in the vicinity of slow enzymes and transporters in the steady state. The analysis can be readily applied to any other molecule, provided knowledge is available about rate of production, average concentration, and diffusion coefficient. Our main conclusion is that the notion of cytosolic compartmentation of metabolites needs reevaluation, as it seems to be in conflict with the underlying physical chemistry.     

Comparison of force fields to study the zinc-finger containing protein NPL4, a target for disulfiram in cancer therapy

Molecular dynamics (MD) simulations are a powerful approach to studying the structure and dynamics of proteins related to health and disease. Advances in the MD field allow modeling proteins with high accuracy. However, modeling metal ions and their interactions with proteins is still challenging. NPL4 is a zinc-binding protein and works as a cofactor for p97 to regulate protein homeostasis. NPL4 is of biomedical importance and has been proposed as the target of disulfiram, a drug recently repurposed for cancer treatment. Experimental studies proposed that the disulfiram metabolites, bis-(diethyldithiocarbamate)‑copper and cupric ions, induce NPL4 misfolding and aggregation. However, the molecular details of their interactions with NPL4 and consequent structural effects are still elusive. Here, biomolecular simulations can help to shed light on the related structural details. To apply MD simulations to NPL4 and its interaction with copper the first important step is identifying a suitable force field to describe the protein in its zinc-bound states. We examined different sets of non-bonded parameters because we want to study the misfolding mechanism and cannot rule out that the zinc may detach from the protein during the process and copper replaces it. We investigated the force-field ability to model the coordination geometry of the metal ions by comparing the results from MD simulations with optimized geometries from quantum mechanics (QM) calculations using model systems of NPL4. Furthermore, we investigated the performance of a force field including bonded parameters to treat copper ions in NPL4 that we obtained based on QM calculations.     

Extensive in vivo metabolite-protein interactions revealed by large-scale systematic analyses

Natural small compounds comprise most cellular molecules and bind proteins as substrates, products, cofactors, and ligands. However, a large-scale investigation of in?vivo protein-small metabolite interactions has not been performed. We developed a mass spectrometry assay for the large-scale identification of in?vivo protein-hydrophobic small metabolite interactions in yeast and analyzed compounds that bind ergosterol biosynthetic proteins and protein kinases. Many of these proteins bind small metabolites; a few interactions were previously known, but the vast majority are new. Importantly, many key regulatory proteins such as protein kinases bind metabolites. Ergosterol was found to bind many proteins and may function as a general regulator. It is required for the activity of Ypk1, a mammalian AKT/SGK kinase homolog. Our study defines potential key regulatory steps in lipid biosynthetic pathways and suggests that small metabolites may play a more general role as regulators of protein activity and function than previously appreciated.     

Alcohol-induced conformational changes of ubiquitin

Ubiquitin has been found to be soluble in ethylene glycol and alcohols as the perchlorate or hydrochloride salt. When the effect of alcohol on the structure of ubiquitin is examined, two reversible conformational transitions are observed. Upon lowering the dielectric constant of aqueous alcohol solutions of ubiquitin from 80 to 45, the native structure of ubiquitin is converted to a form consistent with 50% helical structure. This conformational change results in a change in exposure to solvent of the single methionine and the single tyrosine residues of ubiquitin. In agreement with crystallographic results, these residues are buried in the native conformation but become fully exposed to solvent upon undergoing this transition. Further lowering of the dielectric constant to 20 results in the accumulation of a conformation with almost complete helical structure. Thus, hydrophobic interactions cause facile conformational changes in the ubiquitin structure. These results are discussed in terms of a preferential solvation model. It is shown that the results obtained with different alcohols can be normalized by the use of a dielectric constant scale. This normalization corrects for the different molar volumes of different alcohols, allows comparison of results obtained with different alcohols, and should be useful in studying this phenomenon with different proteins.     

Dual Mechanism of Ion Permeation through VDAC Revealed with Inorganic Phosphate Ions and Phosphate Metabolites

In the exchange of metabolites and ions between the mitochondrion and the cytosol, the voltage-dependent anion channel (VDAC) is a key element, as it forms the major transport pathway for these compounds through the mitochondrial outer membrane. Numerous experimental studies have promoted the idea that VDAC acts as a regulator of essential mitochondrial functions. In this study, using a combination of molecular dynamics simulations, free-energy calculations, and electrophysiological measurements, we investigated the transport of ions through VDAC, with a focus on phosphate ions and metabolites. We showed that selectivity of VDAC towards small anions including monovalent phosphates arises from short-lived interactions with positively charged residues scattered throughout the pore. In dramatic contrast, permeation of divalent phosphate ions and phosphate metabolites (AMP and ATP) involves binding sites along a specific translocation pathway. This permeation mechanism offers an explanation for the decrease in VDAC conductance measured in the presence of ATP or AMP at physiological salt concentration. The binding sites occur at similar locations for the divalent phosphate ions, AMP and ATP, and contain identical basic residues. ATP features a marked affinity for a central region of the pore lined by two lysines and one arginine of the N-terminal helix. This cluster of residues together with a few other basic amino acids forms a “charged brush” which facilitates the passage of the anionic metabolites through the pore. All of this reveals that VDAC controls the transport of the inorganic phosphates and phosphate metabolites studied here through two different mechanisms.     

The steroid ligands of estrogen binding proteins in Xenopus laevis liver cytosol are primarily metabolites of 17β-estradiol

The interactions in vitro between [³H]estradiol and liver proteins from Xenopus laevis have been examined to determine if the binding reaction meets criteria of steroid-receptors which may function in the induction of vitellogenesis. Estrogenic hormones associated with proteins in serum and liver cytosol from Xenopus laevis. However, the interactions between soluble liver proteins and estrogens apparently do not result from serum contamination of liver as specific binding was distinguishable by ligand affinity and by differential mobility on polyacrylamide gels. Steroid ligands bound by liver proteins during incubation in vitro were examined by solubility and by thin-layer chromatography. Only a small percentage (13%) of the bound radioactive ligand was recovered as the original tritium-labeled steroid, 17β-estradiol. The major ligand was recovered as a water-soluble metabolite of estradiol which was identified tentatively as an estradiol-glucoside. To investigate whether the protein-bound estradiol metabolite(s) merely masks a small amount of authentic estradiol-receptor complexes or if the metabolite could be an intermediate in estrogen function, isolated liver nuclei were incubated with liver cytosol containing ³H-labeled steroid-protein complexes or with serum protein-bound [³H]estradiol. Nuclei preferentially accumulated ³H-labelea steroids from liver cytosol protein-steroid complexes relative to [³H]estradiol from serum proteins. However, analysis of the steroids recovered in the nuclei after incubation with liver cytosol revealed that both 17β-[³H]estradiol and the ³H-labeled water-soluble metabolite were retained in vitro by nuclei.     

An atomistic view of the YiiP structural changes upon zinc(II) binding

BackgroundYiiP is a bacterial zinc-for-proton antiporter belonging to the cation diffusion facilitator family. The zinc(II) ions are transported across the cell membrane, from the cytosol to the extracellular space.MethodsWe performed atomistic molecular dynamics simulations of the YiiP dimer with zinc(II) ions in solution to elucidate how the metal ions interact with the protein while moving from the cytosol to the transport site.ResultsWe observed that of the two cavities of the dimer, only one was accessible from the cytosol during transport. Zinc(II) binding to D49 of the transport site triggered a rearrangement of the transmembrane domain that closed the accessible cavity. Finally, we analyzed the free-energy profiles of metal transit in the channel and observed the existence of a high barrier preventing release from the transport site.ConclusionsThe observed dynamics is consistent with the dimer-dimer interface forming a stable scaffold against which the rest of the trans-membrane rearranges.General significanceZinc(II) transporters are present in all kingdoms of life. The present study highlights structural features that might be of general relevance.     

Serpula lacrymans reacts with a general,unspecialized chemical response during interaction with mycoparasitic Trichoderma spp. and bacteria

Serpula lacrymans causes structural damage in timber via cellulose and hemicellulose decomposition. Serpula lacrymans shares its habitat with bacteria and fungi, interactions that often result in stress and changed behaviour of the partners. A typical response to environmental stress is the production of bioactive metabolites and pigments. The diversity and function of these metabolites in inter- and intra-kingdom interactions, is largely elusive. Using dual culture approaches and secondary metabolite fingerprinting with high performance thin layer chromatography (HPTLC) it was demonstrated that stress induces an unspecific vs. Similar metabolite patterns were identified when S. lacrymans was co-cultured with different bacteria and mycoparasitic fungi (Trichoderma spp.). By analysing the metabolites produced along the gradient from the contact zone to more remote parts of the mycelium, we show a systemic reaction of S. lacrymans and that direct contact with other microbes is not a necessity to induce secondary metabolites.     

Reaction within the coulomb-cage; science in retrospect

The Coulomb-cage is defined as the space where the electrostatic interaction between two bodies is more intensive than the thermal energy (k_BT). For small molecule, the Coulomb-cage is a small sphere, extending only few water molecules towards the bulk and its radius is sensitive to the ionic strength of the solution. For charged proteins or membranal structures, the Coulomb-cage can engulf large fraction of the surface and provides a preferred pathway for ion propagation along the surface. Similarly, electrostatic potential at the inner space of a channel can form preferential trajectories passage for ions. The dynamics of ions inside the Coulomb-cage of ions was formulated by the studies of proton-anion recombination of excited photoacids. In the present article, we recount the study of intra- Coulomb-cage reaction taking place on the surface of macro-molecular bodies like micelles, membranes, proteins and intra-protein cavities. The study progressed stepwise, tracing the dynamics of a proton ejected from a photo-acid molecule located at defined sites (on membrane, inter-membrane space, active site of enzyme, inside Large Pore Channels etc.). Accumulation of experimental observations encouraged us to study of the reaction mechanism by molecular dynamics simulations of ions within the Coulomb-cage of proteins surface or inside large pores. The intra-Coulomb-cage proton transfer events follows closely the fine structure of the electrostatic field inside the cage and reflects the shape of nearby dielectric boundaries, the temporal ordering of the solvent molecules and the structural fluctuations of the charged side chains. The article sums some 40 years of research, which in retrospect clarifies the intra-Coulomb-cage reaction mechanism.     

Solution structure of a CUE-ubiquitin complex reveals a conserved mode of ubiquitin binding

Monoubiquitination serves as a regulatory signal in a variety of cellular processes. Monoubiquitin signals are transmitted by binding to a small but rapidly expanding class of ubiquitin binding motifs. Several of these motifs, including the CUE domain, also promote intramolecular monoubiquitination. The solution structure of a CUE domain of the yeast Cue2 protein in complex with ubiquitin reveals intermolecular interactions involving conserved hydrophobic surfaces, including the Leu8-Ile44-Val70 patch on ubiquitin. The contact surface extends beyond this patch and encompasses Lys48, a site of polyubiquitin chain formation. This suggests an occlusion mechanism for inhibiting polyubiquitin chain formation during monoubiquitin signaling. The CUE domain shares a similar overall architecture with the UBA domain, which also contains a conserved hydrophobic patch. Comparative modeling suggests that the UBA domain interacts analogously with ubiquitin. The structure of the CUE-ubiquitin complex may thus serve as a paradigm for ubiquitin recognition and signaling by ubiquitin binding proteins.     

Role of the malate-aspartate shuttle on the metabolic response to myocardial ischemia

The malate-aspartate (M-A) shuttle provides an important mechanism to regulate glycolysis and lactate metabolism in the heart by transferring reducing equivalents from cytosol into mitochondria. However, experimental characterization of the M-A shuttle has been incomplete because of limitations in quantifying cytosolic and mitochondrial metabolites. In this study, we developed a multi-compartment model of cardiac metabolism with detailed presentation of the M-A shuttle to quantitatively predict non-observable fluxes and metabolite concentrations under normal and ischemic conditions in vivo. Model simulations predicted that the M-A shuttle is functionally localized to a subdomain that spans the mitochondrial and cytosolic spaces. With the onset of ischemia, the M-A shuttle flux rapidly decreased to a new steady state in proportion to the reduction in blood flow. Simulation results suggest that the reduced M-A shuttle flux during ischemia was not due to changes in shuttle-associated enzymes and transporters. However, there was a redistribution of shuttle-associated metabolites in both cytosol and mitochondria. Therefore, the dramatic acceleration in glycolysis and the switch to lactate production that occur immediately after the onset of ischemia is mediated by reduced M-A shuttle flux through metabolite redistribution of shuttle associated species across the mitochondrial membrane.     

IONIC RADIUS AND IONIC EFFICIENCY

1. It has been shown in preceding papers that when we separate solutions of electrolytes from pure water by collodion membranes the ions with the same sign of charge as the membrane increase while the ions with the opposite sign of charge diminish the rate of diffusion of water from the side of pure water to the side of solution; and that the accelerating and depressing effects of these ions on the rate of diffusion of water increase with their valency. 2. It is shown in this paper that aside from the valency a second quantity of the ion plays a rôle in this effect, namely the radius, which in a monatomic ion means the distance between the central positive nucleus and the outermost ring or shell of electrons of the ion. In monatomic, monovalent anions the radius increases in the order Cl< Br< I (being largest in I), while for monatomic, monovalent cations it increases in the order Li< Na< K< Rb (being largest in Rb). 3. It is shown that the accelerating as well as the depressing effect of the anions mentioned increases directly with the order of magnitude of their radius and that the efficiency is greatest in the case of I which has the largest radius; while the accelerating as well as the depressing effect of cations increases inversely with the order of magnitude of their radius, Li with the smallest radius having the greatest efficiency. 4. This is intelligible on the assumption that the action of the ions is electrostatic in character, in the case of cations due to the electrostatic effect of the excess charge of their positive nucleus, and in the case of anions due to the excess charge of their captured electron. The electrostatic effect of the positive nucleus of a cation on the membrane (or any other body) will be the greater the smaller the ionic radius of the cation; and the electrostatic effect of an excess electron will be the greater the further its distance from its own positive nucleus. 5. It is suggested that this rule may possibly include polyatomic, monovalent ions (e.g. NO₃, CNS, etc.) when we replace these polyatomic ions by monatomic models in which the radius is calculated in such a way as to give the model the same electrostatic effect which the polyatomic ion possesses. 6. These conclusions are in harmony with the fact that the efficiency of ions increases also with their valency. 7. It is suggested that these rules concerning the influence of the ionic radius can possibly be demonstrated in other phenomena, depending on the electrostatic effect of ions.     

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.