Coarse‐grained simulations of proton‐dependent conformational changes in lactose permease

During lactose/H⁺ symport, the Escherichia coli lactose permease (LacY) undergoes a series of global conformational transitions between inward‐facing (open to cytoplasmic side) and outward‐facing (open to periplasmic side) states. However, the exact local interactions and molecular mechanisms dictating those large‐scale structural changes are not well understood. All‐atom molecular dynamics simulations have been performed to investigate the molecular interactions involved in conformational transitions of LacY, but the simulations can only explore early or partial global structural changes because of the computational limits (< 100 ns). In this work, we implement a hybrid force field that couples the united‐atom protein models with the coarse‐grained MARTINI water/lipid, to investigate the proton‐dependent dynamics and conformational changes of LacY. The effects of the protonation states on two key glutamate residues (Glu325 and Glu269) have been studied. Our results on the salt‐bridge dynamics agreed with all‐atom simulations at early short time period, validating our simulations. From our microsecond simulations, we were able to observe the complete transition from inward‐facing to outward‐facing conformations of LacY. Our results showed that all helices have participated during the global conformational transitions and helical movements of LacY. The inter‐helical distances measured in our simulations were consistent with the double electron‐electron resonance experiments at both cytoplasmic and periplasmic sides. Our simulations indicated that the deprotonation of Glu325 induced the opening of the periplasmics side and partial closure of the cytoplasmic side of LacY, while protonation of the Glu269 caused a stable cross‐domain salt‐bridge (Glu130‐Arg344) and completely closed the cytoplasmic side. Proteins 2016; 84:1067–1074. © 2016 Wiley Periodicals, Inc.     

Helix dynamics in a membrane transport protein: comparative simulations of the glycerol-3-phosphate transporter and its constituent helices

The glycerol-3-phosphate transporter (GlpT) is a member of the major facilitator superfamily (MFS). GlpT is an organic phosphate/inorganic phosphate antiporter. It shares a similar fold with other MFS transporters (e.g. LacY and EmrD) consisting of 12 transmembrane (TM) helices which form two domains (each of six TM helices) surrounding a central ligand-binding cavity. The TM helices (especially the cavity-lining helices) contain a large number of proline and glycine residues, which may aid in the conformational changes believed to underline the transport mechanism. Molecular dynamics simulations in a phospholipid bilayer have been used to compare the conformational properties of the isolated TM helices with those in the intact GlpT protein. Analysis of these simulations focuses on the role of proline-induced flexibility in the TM helices. Our results are consistent with the proposed rocker switch mechanism for transport by GlpT. In particular, the simulations highlight the cavity-lining helices (H4, H5, H10 and H11) as being significantly flexible, suggesting that the transport mechanism may involve intra-helix motions in addition to pseudo-rigid body motions of the N- and C-terminal domains relative to one another.     

Asymmetric stability among the transmembrane helices of lactose permease

Combining structure determinations from nuclear magnetic resonance (NMR) data and molecular dynamics simulations (MD) under the same environmental conditions revealed a startling asymmetry in the intrinsic conformational stability of secondary structure in the transmembrane domain of lactose permease (LacY). Eleven fragments, corresponding to transmembrane segments (TMs) of LacY, were synthesized, and their secondary structure in solution was determined by NMR. Eight of the TMs contained significant regions of helical structure. MD simulations, both in DMSO and in a DMPC bilayer, showed sites of local stability of helical structure in these TMs, punctuated by regions of conformational instability, in substantial agreement with the NMR data. Mapping the stable regions onto the crystal structure of LacY reveals a marked asymmetry, contrasting with the pseudosymmetry in the static structure: the secondary structure in the C-terminal half is more stable than in the N-terminal half. The relative stability of secondary structure is likely exploited in the transport mechanism of LacY. Residues supporting proton conduction are in more stable regions of secondary structure, while residues key to substrate binding are found in considerably unstable regions of secondary structure.     

Probing the Periplasmic-Open State of Lactose Permease in Response to Sugar Binding and Proton Translocation

Based on the crystal structure of lactose permease (LacY) open to the cytoplasm, a hybrid molecular simulation approach with self-guided Langevin dynamics is used to describe conformational changes that lead to a periplasmic-open state. This hybrid approach consists of implicit (IM) and explicit (EX) membrane simulations and requires self-guided Langevin dynamics to enhance protein motions during the IM simulations. The pore radius of the lumen increases by 3.5 Å on the periplasmic side and decreases by 2.5 Å on the cytoplasmic side (relative to the crystal structure), suggesting a lumen that is fully open to the periplasm to allow for extracellular sugar transport and closed to the cytoplasm. Based on our simulations, the mechanism that triggers this conformational change to the periplasmic-open state is the protonation of Glu269 and binding of the disaccharide. Then, helix packing is destabilized by breaking of several side chains involved in hydrogen bonding (Asn245, Ser41, Glu374, Lys42, and Gln242). For the periplasmic-open conformations obtained from our simulations, helix-helix distances agree well with experimental measurements using double electron-electron resonance, fluorescence resonance energy transfer, and varying sized cross-linkers. The periplasmic-open conformations are also in compliance with various substrate accessibility/reactivity measurements that indicate an opening of the protein lumen on the periplasmic side on sugar binding. The comparison with these measurements suggests a possible incomplete closure of the cytoplasmic half in our simulations. However, the closure is sufficient to prevent the disaccharide from transporting to the cytoplasm, which is in accordance with the well-established alternating access model. Ser53, Gln60, and Phe354 are determined to be important in sugar transport during the periplasmic-open stage of the sugar transport cycle and the sugar is found to undergo an orientational change in order to escape the protein lumen.     

Phosphatidylethanolamine and monoglucosyldiacylglycerol are interchangeable in supporting topogenesis and function of the polytopic membrane protein lactose permease

To determine the specific role lipids play in membrane protein topogenesis in vivo, the orientation with respect to the membrane bilayer of Escherichia coli lactose permease (LacY) transmembrane (TM) domains and their flanking extramembrane domains was compared after assembly in native membranes and membranes with genetically modified lipid content using the substituted cysteine accessibility method for determining TM domain mapping. LacY assembled in the absence of the major membrane lipid phosphatidylethanolamine (PE) does not carry out uphill transport of substrate and displays an inverted orientation for the N-terminal six-TM domain helical bundle (Bogdanov, M., Heacock, P. N., and Dowhan, W. (2002) EMBO J. 21, 2107-2116). Strikingly, the replacement of PE in vivo by the foreign lipid monoglucosyldiacylglycerol (MGlcDAG), synthesized by the Acholeplasma laidlawii MGlcDAG synthase, restored uphill transport and supported the wild type TM topology of the N-terminal helical bundle of LacY. An interchangeable role in defining membrane protein TM domain orientation and supporting function is played by the two most abundant lipids, PE and MGlcDAG, in gram-negative and gram-positive bacteria, respectively. Therefore, these structurally diverse lipids endow the membrane with similar properties necessary for the proper organization of protein domains in LacY that are highly sensitive to lipids as topological determinants.     

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

A novel mechanism for membrane modulation of transmembrane protein structure, and consequently function, is suggested in which mismatch between the hydrophobic surface of the protein and the hydrophobic interior of the lipid bilayer induces a flexing or bending of a transmembrane segment of the protein. Studies on model hydrophobic transmembrane peptides predict that helices tilt to submerge the hydrophobic surface within the lipid bilayer to satisfy the hydrophobic effect if the helix length exceeds the bilayer width. The hydrophobic surface of transmembrane helix 1 (TM1) of lactose permease, LacY, is accessible to the bilayer, and too long to be accommodated in the hydrophobic portion of a typical lipid bilayer if oriented perpendicular to the membrane surface. Hence, nuclear magnetic resonance (NMR) data and molecular dynamics simulations show that TM1 from LacY may flex as well as tilt to satisfy the hydrophobic mismatch with the bilayer. In an analogous study of the hydrophobic mismatch of TM7 of bovine rhodopsin, similar flexing of the transmembrane segment near the conserved NPxxY sequence is observed. As a control, NMR data on TM5 of lacY, which is much shorter than TM1, show that TM5 is likely to tilt, but not flex, consistent with the close match between the extent of hydrophobic surface of the peptide and the hydrophobic thickness of the bilayer. These data suggest mechanisms by which the lipid bilayer in which the protein is embedded modulates conformation, and thus function, of integral membrane proteins through interactions with the hydrophobic transmembrane helices.     

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

A novel mechanism for membrane modulation of transmembrane protein structure, and consequently function, is suggested in which mismatch between the hydrophobic surface of the protein and the hydrophobic interior of the lipid bilayer induces a flexing or bending of a transmembrane segment of the protein. Studies on model hydrophobic transmembrane peptides predict that helices tilt to submerge the hydrophobic surface within the lipid bilayer to satisfy the hydrophobic effect if the helix length exceeds the bilayer width. The hydrophobic surface of transmembrane helix 1 (TM1) of lactose permease, LacY, is accessible to the bilayer, and too long to be accommodated in the hydrophobic portion of a typical lipid bilayer if oriented perpendicular to the membrane surface. Hence, nuclear magnetic resonance (NMR) data and molecular dynamics simulations show that TM1 from LacY may flex as well as tilt to satisfy the hydrophobic mismatch with the bilayer. In an analogous study of the hydrophobic mismatch of TM7 of bovine rhodopsin, similar flexing of the transmembrane segment near the conserved NPxxY sequence is observed. As a control, NMR data on TM5 of lacY, which is much shorter than TM1, show that TM5 is likely to tilt, but not flex, consistent with the close match between the extent of hydrophobic surface of the peptide and the hydrophobic thickness of the bilayer. These data suggest mechanisms by which the lipid bilayer in which the protein is embedded modulates conformation, and thus function, of integral membrane proteins through interactions with the hydrophobic transmembrane helices.     

Molecular genetic and biochemical approaches for defining lipid-dependent membrane protein folding

We provide an overview of lipid-dependent polytopic membrane protein folding and topogenesis. Lipid dependence of this process was determined by employing Escherichia coli cells in which specific lipids can be eliminated, substituted, tightly titrated or controlled temporally during membrane protein synthesis and assembly. The secondary transport protein lactose permease (LacY) was used to establish general principles underlying the molecular basis of lipid-dependent effects on protein domain folding, protein transmembrane domain (TM) orientation, and function. These principles were then extended to several other secondary transport proteins of E. coli. The methods used to follow proper conformational organization of protein domains and the topological organization of protein TMs in whole cells and membranes are described. The proper folding of an extramembrane domain of LacY that is crucial for energy dependent uphill transport function depends on specific lipids acting as non-protein molecular chaperones. Correct TM topogenesis is dependent on charge interactions between the cytoplasmic surface of membrane proteins and a proper balance of the membrane surface net charge defined by the lipid head groups. Short-range interactions between the nascent protein chain and the translocon are necessary but not sufficient for establishment of final topology. After release from the translocon short-range interactions between lipid head groups and the nascent protein chain, partitioning of protein hydrophobic domains into the membrane bilayer, and long-range interactions within the protein thermodynamically drive final membrane protein organization. Given the diversity of membrane lipid compositions throughout nature, it is tempting to speculate that during the course of evolution the physical and chemical properties of proteins and lipids have co-evolved in the context of the lipid environment of membrane systems in which both are mutually dependent on each other for functional organization of proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.     

Large-scale conformational dynamics of the HIV-1 integrase core domain and its catalytic loop mutants

HIV-1 integrase is one of the three essential enzymes required for viral replication and has great potential as a novel target for anti-HIV drugs. Although tremendous efforts have been devoted to understanding this protein, the conformation of the catalytic core domain around the active site, particularly the catalytic loop overhanging the active site, is still not well characterized by experimental methods due to its high degree of flexibility. Recent studies have suggested that this conformational dynamics is directly correlated with enzymatic activity, but the details of this dynamics is not known. In this study, we conducted a series of extended-time molecular dynamics simulations and locally enhanced sampling simulations of the wild-type and three loop hinge mutants to investigate the conformational dynamics of the core domain. A combined total of >480 ns of simulation data was collected which allowed us to study the conformational changes that were not possible to observe in the previously reported short-time molecular dynamics simulations. Among the main findings are a major conformational change (>20 A) in the catalytic loop, which revealed a gatinglike dynamics, and a transient intraloop structure, which provided a rationale for the mutational effects of several residues on the loop including Q(148), P(145), and Y(143). Further, clustering analyses have identified seven major conformational states of the wild-type catalytic loop. Their implications for catalytic function and ligand interaction are discussed. The findings reported here provide a detailed view of the active site conformational dynamics and should be useful for structure-based inhibitor design for integrase.     

Molecular dynamics simulations of the ligand-binding domain of an N-methyl-D-aspartate receptor

The mechanism of partial agonism at N-methyl-D-aspartate receptors is an unresolved issue, especially with respect to the role of protein dynamics. We have performed multiple molecular dynamics simulations (7 x 20 ns) to examine the behavior of the ligand-binding core of the NR1 subunit with a series of ligands. Our results show that water plays an important role in stabilizing different conformations of the core and how a closed cleft conformation of the protein might be stabilized in the absence of ligands. In the case of ligand-bound simulations with both full and partial agonists, we observed that ligands within the binding cleft may undergo distinct conformational changes, without grossly influencing the degree of cleft closure within the ligand-binding domain. In agreement with recently published crystallographic data, we also observe similar changes in backbone torsions corresponding to the hinge region between the two lobes for the partial agonist, D-cycloserine. This observation rationalizes the classification of D-cycloserine as a partial agonist and should provide a basis with which to predict partial agonism in this class of receptor by analyzing the behavior of these torsions with other potential ligands.     

Lipid-protein interactions as determinants of membrane protein structure and function

To determine how the lipid environment affects membrane protein structure and function, strains of Escherichia coli were developed in which normal phospholipid composition can be altered or foreign lipids can be introduced. The properties of LacY (lactose permease) were investigated as a function of lipid environment. Assembly of LacY in membranes lacking PE (phosphatidylethanolamine) results in misorientation of the N-terminal six-TM (transmembrane domain) helical bundle with loss of energy-dependent uphill transport and retention of energy-independent downhill transport. Post-assembly introduction of PE results in nearly native orientation of TMs and restoration of uphill transport. Foreign lipids with no net charge can substitute for PE in supporting native LacY topology, but restoration of uphill transport is dependent on native topology and the proper folding of a solvent-exposed domain. Increasing the positive charge density of the cytoplasmically exposed surface of LacY counters TM misorientation in the absence of neutral lipids, demonstrating that charge interactions between these domains and the surface of the membrane bilayer are determinants of TM orientation. Therefore membrane protein organization or reorganization is determined either during initial assembly or post-insertionally through direct interactions between the protein and the lipid environment, which affects the topogenic potency of opposing charged residues as topological signals independent of the translocon.     

Transmembrane signaling of chemotaxis receptor tar: insights from molecular dynamics simulation studies

Transmembrane signaling of chemotaxis receptors has long been studied, but how the conformational change induced by ligand binding is transmitted across the bilayer membrane is still elusive at the molecular level. To tackle this problem, we carried out a total of 600-ns comparative molecular dynamics simulations (including model-building simulations) of the chemotaxis aspartate receptor Tar (a part of the periplasmic domain/transmembrane domain/HAMP domain) in explicit lipid bilayers. These simulations reveal valuable insights into the mechanistic picture of Tar transmembrane signaling. The piston-like movement of a transmembrane helix induced by ligand binding on the periplasmic side is transformed into a combination of both longitudinal and transversal movements of the helix on the cytoplasmic side as a result of different protein-lipid interactions in the ligand-off and ligand-on states of the receptor. This conformational change alters the dynamics and conformation of the HAMP domain, which is presumably a mechanism to deliver the signal from the transmembrane domain to the cytoplasmic domain. The current results are consistent with the previously suggested dynamic bundle model in which the HAMP dynamics change is a key to the signaling. The simulations provide further insights into the conformational changes relevant to the HAMP dynamics changes in atomic detail.     

Transmembrane signal transduction of the alpha(IIb)beta(3) integrin

Integrins are composed of noncovalently bound dimers of an alpha- and a beta-subunit. They play an important role in cell-matrix adhesion and signal transduction through the cell membrane. Signal transduction can be initiated by the binding of intracellular proteins to the integrin. Binding leads to a major conformational change. The change is passed on to the extracellular domain through the membrane. The affinity of the extracellular domain to certain ligands increases; thus at least two states exist, a low-affinity and a high-affinity state. The conformations and conformational changes of the transmembrane (TM) domain are the focus of our interest. We show by a global search of helix-helix interactions that the TM section of the family of integrins are capable of adopting a structure similar to the structure of the homodimeric TM protein Glycophorin A. For the alpha(IIb)beta(3) integrin, this structural motif represents the high-affinity state. A second conformation of the TM domain of alpha(IIb)beta(3) is identified as the low-affinity state by known mutational and nuclear magnetic resonance (NMR) studies. A transition between these two states was determined by molecular dynamics (MD) calculations. On the basis of these calculations, we propose a three-state mechanism.     

Structure and Function in Homodimeric Enzymes: Simulations of Cooperative and Independent Functional Motions

Large-scale conformational change is a common feature in the catalytic cycles of enzymes. Many enzymes function as homodimers with active sites that contain elements from both chains. Symmetric and anti-symmetric cooperative motions in homodimers can potentially lead to correlated active site opening and/or closure, likely to be important for ligand binding and release. Here, we examine such motions in two different domain-swapped homodimeric enzymes: the DcpS scavenger decapping enzyme and citrate synthase. We use and compare two types of all-atom simulations: conventional molecular dynamics simulations to identify physically meaningful conformational ensembles, and rapid geometric simulations of flexible motion, biased along normal mode directions, to identify relevant motions encoded in the protein structure. The results indicate that the opening/closure motions are intrinsic features of both unliganded enzymes. In DcpS, conformational change is dominated by an anti-symmetric cooperative motion, causing one active site to close as the other opens; however a symmetric motion is also significant. In CS, we identify that both symmetric (suggested by crystallography) and asymmetric motions are features of the protein structure, and as a result the behaviour in solution is largely non-cooperative. The agreement between two modelling approaches using very different levels of theory indicates that the behaviours are indeed intrinsic to the protein structures. Geometric simulations correctly identify and explore large amplitudes of motion, while molecular dynamics simulations indicate the ranges of motion that are energetically feasible. Together, the simulation approaches are able to reveal unexpected functionally relevant motions, and highlight differences between enzymes.     

Essential dynamics sampling study of adenylate kinase: comparison to citrate synthase and implication for the hinge and shear mechanisms of domain motions

Essential dynamics sampling simulations of the domain conformations of unliganded Escherichia coli adenylate kinase have been performed to determine whether the ligand-induced closed-domain conformation is accessible to the open unliganded enzyme. Adenylate kinase is a three- domain protein with a central CORE domain and twoflanking domains, the LID and the NMPbind domains. The sampling simulations were applied to the CORE and NMPbind domain pair and the CORE and LID domain pair separately. One aim is to compare the results to those of a similar study on the enzyme citrate synthase to determine whether a similar domain-locking mechanism operates in adenylate kinase. Although for adenylate kinase the simulations suggest that the closed-domain conformation of the unliganded enzyme is at a slightly higher free energy than the open for both domain pairs, the results are radically different to those found for citrate synthase. In adenylate kinase the targeted domain conformations could always be achieved, whereas this was not the case in citrate synthase due to an apparent free-energy barrier between the open and closed conformations. Adenylate kinase has been classified as a protein that undergoes closure through a hinge mechanism, whereas citrate synthase has been assigned to the shear mechanism. This was quantified here in terms of the change in the number of interdomain contacting atoms upon closure which showed a considerable increase in adenylate kinase. For citrate synthase this number remained largely the same, suggesting that the domain faces slide over each other during closure. This suggests that shear and hinge mechanisms of domain closure may relate to the existence or absence of an appreciable barrier to closure for the unliganded protein, as the latter can hinge comparatively freely, whereas the former must follow a more constrained path. In general though it appears a bias toward keeping the unliganded enzyme in the open-domain conformation may be a common feature of domain enzymes.     

The mechanism of nucleotide‐binding domain dimerization in the intact maltose transporter as studied by all‐atom molecular dynamics simulations

The Escherichia coli maltose transporter MalFGK₂‐E belongs to the protein superfamily of ATP‐binding cassette (ABC) transporters. This protein is composed of heterodimeric transmembrane domains (TMDs) MalF and MalG, and the homodimeric nucleotide‐binding domains (NBDs) MalK₂. In addition to the TMDs and NBDs, the periplasmic maltose binding protein MalE captures maltose and shuttle it to the transporter. In this study, we performed all‐atom molecular dynamics (MD) simulations on the maltose transporter and found that both the binding of MalE to the periplasmic side of the TMDs and binding of ATP to the MalK₂ are necessary to facilitate the conformational change from the inward‐facing state to the occluded state, in which MalK₂ is completely dimerized. MalE binding suppressed the fluctuation of the TMDs and MalF periplasmic region (MalF‐P2), and thus prevented the incorrect arrangement of the MalF C‐terminal (TM8) helix. Without MalE binding, the MalF TM8 helix showed a tendency to intrude into the substrate translocation pathway, hindering the closure of the MalK₂. This observation is consistent with previous mutagenesis experimental results on MalF and provides a new point of view regarding the understanding of the substrate translocation mechanism of the maltose transporter.     

Effects of clinically relevant MPL mutations in the transmembrane domain revealed at the atomic level through computational modeling

Background

Mutations in the thrombopoietin receptor (MPL) may activate relevant pathways and lead to chronic myeloproliferative neoplasms (MPNs). The mechanisms of MPL activation remain elusive because of a lack of experimental structures. Modern computational biology techniques were utilized to explore the mechanisms of MPL protein activation due to various mutations.

Results

Transmembrane (TM) domain predictions, homology modeling, ab initio protein structure prediction, and molecular dynamics (MD) simulations were used to build structural dynamic models of wild-type and four clinically observed mutants of MPL. The simulation results suggest that S505 and W515 are important in keeping the TM domain in its correct position within the membrane. Mutations at either of these two positions cause movement of the TM domain, altering the conformation of the nearby intracellular domain in unexpected ways, and may cause the unwanted constitutive activation of MPL''s kinase partner, JAK2.

Conclusions

Our findings represent the first full-scale molecular dynamics simulations of the wild-type and clinically observed mutants of the MPL protein, a critical element of the MPL-JAK2-STAT signaling pathway. In contrast to usual explanations for the activation mechanism that are based on the relative translational movement between rigid domains of MPL, our results suggest that mutations within the TM region could result in conformational changes including tilt and rotation (azimuthal) angles along the membrane axis. Such changes may significantly alter the conformation of the adjacent and intrinsically flexible intracellular domain. Hence, caution should be exercised when interpreting experimental evidence based on rigid models of cytokine receptors or similar systems.     

Coupling between the voltage-sensing and pore domains in a voltage-gated potassium channel

Voltage-dependent potassium (Kv), sodium (Nav), and calcium channels open and close in response to changes in transmembrane (TM) potential, thus regulating cell excitability by controlling ion flow across the membrane. An outstanding question concerning voltage gating is how voltage-induced conformational changes of the channel voltage-sensing domains (VSDs) are coupled through the S4-S5 interfacial linking helices to the opening and closing of the pore domain (PD). To investigate the coupling between the VSDs and the PD, we generated a closed Kv channel configuration from Aeropyrum pernix (KvAP) using atomistic simulations with experiment-based restraints on the VSDs. Full closure of the channel required, in addition to the experimentally determined TM displacement, that the VSDs be displaced both inwardly and laterally around the PD. This twisting motion generates a tight hydrophobic interface between the S4-S5 linkers and the C-terminal ends of the pore domain S6 helices in agreement with available experimental evidence.     

Interactions of the pleckstrin homology domain with phosphatidylinositol phosphate and membranes: characterization via molecular dynamics simulations

The mechanism of interaction of pleckstrin homology (PH) domains with phosphatidylinositol 4,5-bisphosphate (PIP 2)-containing lipid bilayers remains uncertain. While crystallographic studies have emphasized PH-inositol 1,4,5-trisphosphate (IP 3) interactions, biophysical studies indicate a degree of less specific protein-bilayer interactions. We have used molecular dynamics simulations to characterize the interactions of the PH domain from phospholipase C-delta1 with IP 3 and with PIP 2, the latter in lipid bilayers and in detergent micelles. Simulations of the PH domain in water reveal a reduction in protein flexibility when IP 3 is bound. Simulations of the PH domain bound to PIP 2 in lipid bilayers indicate a tightening of ligand-protein interactions relative to the PH-IP 3 complex, alongside formation of H-bonds between PH side chains and lipid (PC) headgroups, and a degree of penetration of hydrophobic side chains into the core of the bilayer. Comparison with simulations of the PH-bound domain to a PC bilayer in the absence of PIP 2 suggests that the presence of PIP 2 increases the extent of PH-membrane interactions. Thus, comparative molecular dynamics simulations reveal how a PI-binding domain undergoes changes in conformational dynamics on binding to a PIP 2-containing membrane and how interactions additional to those with the PI headgroup are formed.     

Flexible Gates Generate Occluded Intermediates in the Transport Cycle of LacY

The major facilitator superfamily (MFS) transporter lactose permease (LacY) alternates between cytoplasmic and periplasmic open conformations to co-transport a sugar molecule together with a proton across the plasma membrane. Indirect experimental evidence suggested the existence of an occluded transition intermediate of LacY, which would prevent leaking of the proton gradient. As no experimental structure is known, the conformational transition is not fully understood in atomic detail. We simulated transition events from a cytoplasmic open conformation to a periplasmic open conformation with the dynamic importance sampling molecular dynamics method and observed occluded intermediates. Analysis of water permeation pathways and the electrostatic free-energy landscape of a solvated proton indicated that the occluded state contains a solvated central cavity inaccessible from either side of the membrane. We propose a pair of geometric order parameters that capture the state of the pathway through the MFS transporters as shown by a survey of available crystal structures and models. We present a model for the occluded state of apo-LacY, which is similar to the occluded crystal structures of the MFS transporters EmrD, PepT_So, NarU, PiPT and XylE. Our simulations are consistent with experimental double electron spin–spin distance measurements that have been interpreted to show occluded conformations. During the simulations, a salt bridge that has been postulated to be involved in driving the conformational transition formed. Our results argue against a simple rigid-body domain motion as implied by a strict “rocker-switch mechanism” and instead hint at an intricate coupling between two flexible gates.