Sugar binding and protein conformational changes in lactose permease

Use of Molecular Dynamics Simulations in Structure-Based Drug Discovery

BACKGROUND

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

OBJECTIVE

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

METHOD

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

RESULTS

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

CONCLUSION

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

Structural and Functional Adaptability of Sucrose and Lactose Permeases from Escherichia coli to the Membrane Lipid Composition

The lipid environment in which membrane proteins are embedded can influence their structure and function. Lipid-protein interactions and lipid-induced conformational changes necessary for protein function remain intractable in vivo using high-resolution techniques. Using Escherichia coli strains in which the normal phospholipid composition can be altered or foreign lipids can be introduced, we established the importance of membrane lipid composition for the proper folding, assembly, and function of E. coli lactose (LacY) and sucrose (CscB) permeases. However, the molecular mechanism underlying the lipid dependence for active transport remains unknown. Herein, we demonstrate that the structure and function of CscB and LacY can be modulated by the composition of the lipid environment. Using a combination of assays (transport activity of the substrate, protein topology, folding, and assembly into the membrane), we found that alterations in the membrane lipid composition lead to lipid-dependent structural changes in CscB and LacY. These changes affect the orientation of residues involved in LacY proton translocation and impact the rates of protonation and deprotonation of E325 by affecting the arrangement of transmembrane domains in the vicinity of the R302-E325 charge pair. Furthermore, the structural changes caused by changes in membrane lipid composition can be altered by a single-point mutation, highlighting the adaptability of these transporters to their environment. Altogether, our results demonstrate that direct interactions between a protein and its lipid environment uniquely contribute to membrane protein organization and function. Because members of the major facilitator superfamily present with well-conserved functional architecture, we anticipate that our findings can be extrapolated to other membrane protein transporters.

Transport cycle of Escherichia coli lactose permease in a nonhomogeneous random walk model

We present Monte Carlo simulations for the transport cycle of Escherichia coli lactose permease (LacY), using as a starting point the model proposed by Kaback et al. [Nat. Rev. Mol. Cell Biol. 2, 610 (2001)NRMCBP1471-007210.1038/35085077], which is based on functional properties of mutants and x-ray structures. Kaback's model suggests the existence of six states for the whole cycle of lactose-H^{+} symport. However, the free-energy differences between these states have not yet been reported in the literature. Here, we analyzed the biochemical structure of each state and determined a range of possible values for each one of the five free-energy variations. Then, using the Metropolis algorithm in a nonhomogeneous random walk model, we tested all the possible combinations with these values to find the free-energy curve that best reproduces the dynamics of LacY. The agreement between our model predictions and the experimental data suggests that our free-energy curve is appropriate for describing the lactose-H^{+} symport. We found not only this curve, but also the time of occupancy of the permease at each conformation. In addition, we paved the way in this work to solve an open question related to this transport mechanism, which is the importance of protonation for lactose binding.

Uptake dynamics in the Lactose permease (LacY) membrane protein transporter

The sugar transporter Lactose permease (LacY) of Escherichia coli has become a prototype to understand the underlying molecular details of membrane transport. Crystal structures have trapped the protein in sugar-bound states facing the periplasm, but with narrow openings unable to accommodate sugar. Therefore, the molecular details of sugar uptake remain elusive. In this work, we have used extended simulations and metadynamics sampling to explore a putative sugar-uptake pathway and associated free energy landscape. We found an entrance at helix-pair 2 and 11, which involved lipid head groups and residues Gln 241 and Gln 359. Furthermore, the protein displayed high flexibility on the periplasmic side of Phe 27, which is located at the narrowest section of the pathway. Interactions to Phe 27 enabled passage into the binding site, which was associated with a 24 ± 4 kJ/mol binding free energy in excellent agreement with an independent binding free energy calculation and experimental data. Two free energy minima corresponding to the two possible binding poses of the lactose analog β-D-galactopyranosyl-1-thio-β-D-galactopyranoside (TDG) were aligned with the crystal structure-binding pocket. This work outlines the chemical environment of a putative periplasmic sugar pathway and paves way for understanding substrate affinity and specificity in LacY.

Exploration of conformational changes in lactose permease upon sugar binding and proton transfer through coarse-grained simulations

Escherichia coli lactose permease (LacY) actively transports lactose and other galactosides across cell membranes through lactose/H⁺ symport process. Lactose/H⁺ symport is a highly complex process that involves sugar translocation, H⁺ transfer, and large-scale protein conformational changes. The complete picture of lactose/H⁺ symport is largely unclear due to the complexity and multiscale nature of the process. In this work, we develop the force field for sugar molecules compatible with PACE, a hybrid and coarse-grained force field that couples the united-atom protein models with the coarse-grained MARTINI water/lipid. After validation, we implement the new force field to investigate the binding of a β-d-galactopyranosyl-1-thio- β-d-galactopyranoside (TDG) molecule to a wild-type LacY. Results show that the local interactions between TDG and LacY at the binding pocket are consistent with the X-ray experiment. Transitions from inward-facing to outward-facing conformations upon TDG binding and protonation of Glu269 have been achieved from ∼5.5 µs simulations. Both the opening of the periplasmic side and closure of the cytoplasmic side of LacY are consistent with double electron-electron resonance and thiol cross-linking experiments. Our analysis suggests that the conformational changes of LacY are a cumulative consequence of interdomain H-bonds breaking at the periplasmic side, interdomain salt-bridge formation at the cytoplasmic side, and the TDG orientational changes during the transition.

Effect of External Electric Field on Substrate Transport of a Secondary Active Transporter

Substrate transport across a membrane accomplished by a secondary active transporter (SAT) is essential to the normal physiological function of living cells. In the present research, a series of all-atom molecular dynamics (MD) simulations under different electric field (EF) strengths was performed to investigate the effect of an external EF on the substrate transport of an SAT. The results show that EF both affects the interaction between substrate and related protein's residues by changing their conformations and tunes the timeline of the transport event, which collectively reduces the height of energy barrier for substrate transport and results in the appearance of two intermediate conformations under the existence of an external EF. Our work spotlights the crucial influence of external EFs on the substrate transport of SATs and could provide a more penetrating understanding of the substrate transport mechanism of SATs.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.

In silico pharmacology: Drug membrane partitioning and crossing

Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; β-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.

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.

Probing the importance of lipid diversity in cell membranes via molecular simulation

Lipid membranes in prokaryotes and eukaryotes have a wide array of lipids that are necessary for proper membrane structure and function. In this paper, an introduction to lipid diversity in biology and a mini-review on how molecular simulations have been used to model biological membranes (primarily limited to one to three lipid types in most simulation-based models) is provided, which motivates the use of all-atom molecular dynamics (MD) simulations to study the effect of lipid diversity on properties of realistic membrane models of prokaryotes and eukaryotes. As an example, cytoplasmic membrane models of Escherichia coli were developed at different stages of the colony growth cycle (early-log, mid-log, stationary and overnight). The main difference between lipid compositions at each stage was the concentration of a cyclopropane-containing moiety on the sn-2 lipid acyl chain (cyC17:0). Triplicate MD simulations for each stage were run for 300 ns to study the influence of lipid diversity on the surface area per lipid, area compressibility modulus, deuterium order parameters, and electron density profiles. The overnight stage (also known as the death stage) had the highest average surface area per lipid, highest rigidity, and lowest bilayer thickness compare to other stages of E. coli cytoplasmic membrane. Although bilayer thickness did depend on the growth stage, the changes between these were small suggesting that the hydrophobic core of transmembrane proteins fit well with the membrane in all growth stages. Although it is still common practise in MD simulations of membrane proteins to use simple one- or two-component membranes, it can be important to use diverse lipid model membranes when membrane protein structure and function are influenced by changes in lipid membrane composition.

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.

•

The transport mechanism of LacY is hypothesized to involve an intermediate “occluded” state.

•

Such a state is observed in computer simulations of the conformational transitions.

•

Simulation data are validated with experimental double electron–electron spin resonance measurements.

•

The structural gating elements of LacY are identified.

•

Occluded LacY is similar to known occluded structures of homologous proteins.

Protonation drives the conformational switch in the multidrug transporter LmrP

Multidrug antiporters of the major facilitator superfamily couple proton translocation to the extrusion of cytotoxic molecules. The conformational changes that underlie the transport cycle and the structural basis of coupling of these transporters have not been elucidated. Here we used extensive double electron-electron resonance measurements to uncover the conformational equilibrium of LmrP, a multidrug transporter from Lactococcus lactis, and to investigate how protons and ligands shift this equilibrium to enable transport. We find that the transporter switches between outward-open and outward-closed conformations, depending on the protonation states of specific acidic residues forming a transmembrane protonation relay. Our data can be framed in a model of transport wherein substrate binding initiates the transport cycle by opening the extracellular side. Subsequent protonation of membrane-embedded acidic residues induces substrate release to the extracellular side and triggers a cascade of conformational changes that concludes in proton release to the intracellular side.

Structural basis for dynamic mechanism of proton-coupled symport by the peptide transporter POT

Proton-dependent oligopeptide transporters (POTs) are major facilitator superfamily (MFS) proteins that mediate the uptake of peptides and peptide-like molecules, using the inwardly directed H(+) gradient across the membrane. The human POT family transporter peptide transporter 1 is present in the brush border membrane of the small intestine and is involved in the uptake of nutrient peptides and drug molecules such as β-lactam antibiotics. Although previous studies have provided insight into the overall structure of the POT family transporters, the question of how transport is coupled to both peptide and H(+) binding remains unanswered. Here we report the high-resolution crystal structures of a bacterial POT family transporter, including its complex with a dipeptide analog, alafosfalin. These structures revealed the key mechanistic and functional roles for a conserved glutamate residue (Glu310) in the peptide binding site. Integrated structural, biochemical, and computational analyses suggested a mechanism for H(+)-coupled peptide symport in which protonated Glu310 first binds the carboxyl group of the peptide substrate. The deprotonation of Glu310 in the inward open state triggers the release of the bound peptide toward the intracellular space and salt bridge formation between Glu310 and Arg43 to induce the state transition to the occluded conformation.

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.

Proton-coupled dynamics in lactose permease

Lactose permease of Escherichia coli (LacY) catalyzes symport of a galactopyranoside and an H⁺ via an alternating access mechanism. The transition from an inward- to an outward-facing conformation of LacY involves sugar-release followed by deprotonation. Because the transition depends intimately upon the dynamics of LacY in a bilayer environment, molecular dynamics (MD) simulations may be the only means of following the accompanying structural changes in atomic detail. Here, we describe MD simulations of wild-type apo LacY in phosphatidylethanolamine (POPE) lipids that features two protonation states of the critical Glu325. While the protonated system displays configurational stability, deprotonation of Glu325 causes significant structural rearrangements that bring into proximity side chains important for H⁺ translocation and sugar binding and closes the internal cavity. Moreover, protonated LacY in phosphatidylcholine (DMPC) lipids shows that the observed dynamics are lipid-dependent. Together, the simulations describe early dynamics of the inward-to-outward transition of LacY that agree well with experimental data.

Simulations of substrate transport in the multidrug transporter EmrD

EmrD is a multidrug resistance (MDR) transporter from Escherichia coli, which is involved in the efflux of amphipathic compounds from the cytoplasm, and the first MDR member of the major facilitator superfamily to be crystallized. Molecular dynamics simulation of EmrD in a phospholipid bilayer was used to characterize the conformational dynamics of the protein. Motions that support a previously proposed lateral diffusion pathway for substrate from the cytoplasmic membrane leaflet into the EmrD central cavity were observed. In addition, the translocation pathway of meta-chloro carbonylcyanide phenylhydrazone (CCCP) was probed using both standard and steered molecular dynamics simulation. In particular, interactions of a few specific residues with CCCP have been identified. Finally, a large motion of two residues, Val 45 and Leu 233, was observed with the passage of CCCP into the periplasmic space, placing a lower bound on the extent of opening required at this end of the protein for substrate transport. Overall, our simulations probe details of the transport pathway, motions of EmrD at an atomic level of detail, and offer new insights into the functioning of MDR transporters.

Exploring substrate diffusion in channels using biased molecular dynamics simulations

Substrate transport and diffusion through membrane-bound channels are processes that can span a range of time scales, with only the fastest ones being amenable to most atomic-scale equilibrium molecular dynamics (MD) simulations. However, the application of forces within a simulation can greatly accelerate diffusion processes, revealing important structural and energetic features of the channel. Here, we demonstrate the use of two methods for applying biases to a substrate in a simulation, using the ammonia/ammonium transporter AmtB as an example. The first method, steered MD, applies a constant force or velocity constraint to the substrate, permitting the exploration of potential substrate pathways and the barriers encountered, although typically far outside of equilibrium. On the other hand, the second method, adaptive biasing forces, is quasi-equilibrium, permitting the derivation of a potential of mean force, which characterizes the free energy of the substrate during transport.

Lactose permease and the alternating access mechanism

Crystal structures of the lactose permease of Escherichia coli (LacY) reveal 12, mostly irregular transmembrane α-helices surrounding a large cavity open to the cytoplasm and a tightly sealed periplasmic side (inward-facing conformation) with the sugar-binding site at the apex of the cavity and inaccessible from the periplasm. However, LacY is highly dynamic, and binding of a galactopyranoside causes closing of the inward-facing cavity with opening of a complementary outward-facing cavity. Therefore, the coupled, electrogenic translocation of a sugar and a proton across the cytoplasmic membrane via LacY very likely involves a global conformational change that allows alternating access of sugar- and H(+)-binding sites to either side of the membrane. Here the various biochemical and biophysical approaches that provide strong support for the alternating access mechanism are reviewed. Evidence is also presented indicating that opening of the periplasmic cavity is probably the limiting step for binding and perhaps transport.

Nucleoside/nucleobase transporters: drug targets of the future?

BACKGROUND

Nucleoside/nucleobase transporters have been investigated since the 1960s. In particular, equilibrative nucleoside transporters were thought to be valuable drug targets, since they are involved in various kinds of viral and parasitic diseases as well as cancers.

DISCUSSION

In the postgenomic era multiple transporters, including different subtypes, have been cloned and characterized on the molecular level. In this article we summarize recent advances regarding structure, function and localization of nucleoside/nucleobase transporters as well as the pharmacological profile of selected drugs.

CONCLUSION

Knowledge of the different kinetic properties and structural features of nucleoside transporters can either be used for the rational design of therapeutics directly targeting the transporter itself or for the delivery of drugs using the transporter as a port of entry into the target cell. Equilibrative nucleoside transporters are of considerable pharmacological interest as drug targets for the development of drugs tailored to each patient's need for the treatment of cardiac disease, cancer and viral infections.

The alternating-access mechanism of MFS transporters arises from inverted-topology repeats

Lactose permease (LacY) is the prototype of the major facilitator superfamily (MFS) of secondary transporters. Available structures of LacY reveal a state in which the substrate is exposed to the cytoplasm but is occluded from the periplasm. However, the alternating-access transport mechanism requires the existence of a periplasm-facing state. We recently showed that inverted-topology structural repeats provide the foundation for the mechanisms of two transporter families with folds distinct from the MFS. Here, we generated a structural model of LacY by swapping the conformations of inverted-topology repeats identified in its two domains. The model exhibits all required properties of an outward-facing conformation, i.e., closure of the binding site to the cytoplasm and exposure to the periplasm. Furthermore, the model agrees with double electron-electron resonance distance changes, accessibility to cysteine-modifying reagents, cysteine cross-linking data, and a recent structure of a distantly related transporter. Analysis of the intradomain differences between the two states suggests a role for conserved sequence motifs in occluding the central pathway through kinking of the pore-lining helices. In addition, predicted re-pairing of critical salt-bridging residues in the binding sites agrees remarkably well with previous proposals, allowing a description of the proton/sugar transport mechanism. More fundamentally, our model demonstrates that inverted-topology repeats provide the foundation for the alternating-access mechanisms of MFS transporters.

Molecular dynamics simulation studies of GLUT4: substrate-free and substrate-induced dynamics and ATP-mediated glucose transport inhibition

Background

Glucose transporter 4 (GLUT4) is an insulin facilitated glucose transporter that plays an important role in maintaining blood glucose homeostasis. GLUT4 is sequestered into intracellular vesicles in unstimulated cells and translocated to the plasma membrane by various stimuli. Understanding the structural details of GLUT4 will provide insights into the mechanism of glucose transport and its regulation. To date, a crystal structure for GLUT4 is not available. However, earlier work from our laboratory proposed a well validated homology model for GLUT4 based on the experimental data available on GLUT1 and the crystal structure data obtained from the glycerol 3-phosphate transporter.

Methodology/Principal Findings

In the present study, the dynamic behavior of GLUT4 in a membrane environment was analyzed using three forms of GLUT4 (apo, substrate and ATP-substrate bound states). Apo form simulation analysis revealed an extracellular open conformation of GLUT4 in the membrane favoring easy exofacial binding of substrate. Simulation studies with the substrate bound form proposed a stable state of GLUT4 with glucose, which can be a substrate-occluded state of the transporter. Principal component analysis suggested a clockwise movement for the domains in the apo form, whereas ATP substrate-bound form induced an anti-clockwise rotation. Simulation studies suggested distinct conformational changes for the GLUT4 domains in the ATP substrate-bound form and favor a constricted behavior for the transport channel. Various inter-domain hydrogen bonds and switching of a salt-bridge network from E345-R350-E409 to E345-R169-E409 contributed to this ATP-mediated channel constriction favoring substrate occlusion and prevention of its release into cytoplasm. These data are consistent with the biochemical studies, suggesting an inhibitory role for ATP in GLUT-mediated glucose transport.

Conclusions/Significance

In the absence of a crystal structure for any glucose transporter, this study provides mechanistic details of the conformational changes in GLUT4 induced by substrate and its regulator.

The structural basis of secondary active transport mechanisms

Secondary active transporters couple the free energy of the electrochemical potential of one solute to the transmembrane movement of another. As a basic mechanistic explanation for their transport function the model of alternating access was put forward more than 40 years ago, and has been supported by numerous kinetic, biochemical and biophysical studies. According to this model, the transporter exposes its substrate binding site(s) to one side of the membrane or the other during transport catalysis, requiring a substantial conformational change of the carrier protein. In the light of recent structural data for a number of secondary transport proteins, we analyze the model of alternating access in more detail, and correlate it with specific structural and chemical properties of the transporters, such as their assignment to different functional states in the catalytic cycle of the respective transporter, the definition of substrate binding sites, the type of movement of the central part of the carrier harboring the substrate binding site, as well as the impact of symmetry on fold-specific conformational changes. Besides mediating the transmembrane movement of solutes, the mechanism of secondary carriers inherently involves a mechanistic coupling of substrate flux to the electrochemical potential of co-substrate ions or solutes. Mainly because of limitations in resolution of available transporter structures, this important aspect of secondary transport cannot yet be substantiated by structural data to the same extent as the conformational change aspect. We summarize the concepts of coupling in secondary transport and discuss them in the context of the available evidence for ion binding to specific sites and the impact of the ions on the conformational state of the carrier protein, which together lead to mechanistic models for coupling.

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.

Identification of specific lipid-binding sites in integral membrane proteins

Protein-lipid interactions are increasingly recognized as central to the structure and function of membrane proteins. However, with the exception of simplified models, specific protein-lipid interactions are particularly difficult to highlight experimentally. Here, we used molecular dynamics simulations to identify a specific protein-lipid interaction in lactose permease, a prototypical model for transmembrane proteins. The interactions can be correlated with the functional dependence of the protein to specific lipid species. The technique is simple and widely applicable to other membrane proteins, and a variety of lipid matrices can be used.

Relating protein conformational changes to packing efficiency and disorder

Changes in protein conformation play key roles in facilitating various biochemical processes, ranging from signaling and phosphorylation to transport and catalysis. While various factors that drive these motions such as environmental changes and binding of small molecules are well understood, specific causative effects on the structural features of the protein due to these conformational changes have not been studied on a large scale. Here, we study protein conformational changes in relation to two key structural metrics: packing efficiency and disorder. Packing has been shown to be crucial for protein stability and function by many protein design and engineering studies. We study changes in packing efficiency during conformational changes, thus extending the analysis from a static context to a dynamic perspective and report some interesting observations. First, we study various proteins that adopt alternate conformations and find that tendencies to show motion and change in packing efficiency are correlated: residues that change their packing efficiency show larger motions. Second, our results suggest that residues that show higher changes in packing during motion are located on the changing interfaces which are formed during these conformational changes. These changing interfaces are slightly different from shear or static interfaces that have been analyzed in previous studies. Third, analysis of packing efficiency changes in the context of secondary structure shows that, as expected, residues buried in helices show the least change in packing efficiency, whereas those embedded in bends are most likely to change packing. Finally, by relating protein disorder to motions, we show that marginally disordered residues which are ordered enough to be crystallized but have sequence patterns indicative of disorder show higher dislocation and a higher change in packing than ordered ones and are located mostly on the changing interfaces. Overall, our results demonstrate that between the two conformations, the cores of the proteins remain mostly intact, whereas the interfaces display the most elasticity, both in terms of disorder and change in packing efficiency. By doing a variety of tests, we also show that our observations are robust to the solvation state of the proteins.

Molecular dynamics simulations of membrane channels and transporters

Membrane transport constitutes one of the most fundamental processes in all living cells with proteins as major players. Proteins as channels provide highly selective diffusive pathways gated by environmental factors, and as transporters furnish directed, energetically uphill transport consuming energy. X-ray crystallography of channels and transporters furnishes a rapidly growing number of atomic resolution structures, permitting molecular dynamics (MD) simulations to reveal the physical mechanisms underlying channel and transporter function. Ever increasing computational power today permits simulations stretching up to 1 micros, that is, to physiologically relevant time scales. Membrane protein simulations presently focus on ion channels, on aquaporins, on protein-conducting channels, as well as on various transporters. In this review we summarize recent developments in this rapidly evolving field.

Identification of two distinct inactive conformations of the beta2-adrenergic receptor reconciles structural and biochemical observations

Fully understanding the mechanisms of signaling proteins such as G protein-coupled receptors (GPCRs) will require the characterization of their conformational states and the pathways connecting those states. The recent crystal structures of the beta(2)- and beta(1)-adrenergic receptors in a nominally inactive state constituted a major advance toward this goal, but also raised new questions. Although earlier biochemical observations had suggested that these receptors possessed a set of contacts between helices 3 and 6, known as the ionic lock, which was believed to form a molecular switch for receptor activation, the crystal structures lacked these contacts. The unexpectedly broken ionic lock has raised questions about the true conformation(s) of the inactive state and the role of the ionic lock in receptor activation and signaling. To address these questions, we performed microsecond-timescale molecular dynamics simulations of the beta(2)-adrenergic receptor (beta(2)AR) in multiple wild-type and mutant forms. In wild-type simulations, the ionic lock formed reproducibly, bringing the intracellular ends of helices 3 and 6 together to adopt a conformation similar to that found in inactive rhodopsin. Our results suggest that inactive beta(2)AR exists in equilibrium between conformations with the lock formed and the lock broken, whether or not the cocrystallized ligand is present. These findings, along with the formation of several secondary structural elements in the beta(2)AR loops during our simulations, may provide a more comprehensive picture of the inactive state of the beta-adrenergic receptors, reconciling the crystal structures with biochemical studies.

Using VMD: an introductory tutorial

VMD (Visual Molecular Dynamics) is a molecular visualization and analysis program designed for biological systems such as proteins, nucleic acids, lipid bilayer assemblies, etc. This unit will serve as an introductory VMD tutorial. We will present several step-by-step examples of some of VMD's most popular features, including visualizing molecules in three dimensions with different drawing and coloring methods, rendering publication-quality figures, animating and analyzing the trajectory of a molecular dynamics simulation, scripting in the text-based Tcl/Tk interface, and analyzing both sequence and structure data for proteins.

Gating at both ends and breathing in the middle: conformational dynamics of TolC

Drug extrusion via efflux through a tripartite complex (an inner membrane pump, an outer membrane protein, and a periplasmic protein) is a widely used mechanism in Gram-negative bacteria. The outer membrane protein (TolC in Escherichia coli; OprM in Pseudomonas aeruginosa) forms a tunnel-like pore through the periplasmic space and the outer membrane. Molecular dynamics simulations of TolC have been performed, and are compared to simulations of Y362F/R367S mutant, and to simulations of its homolog OprM. The results reveal a complex pattern of conformation dynamics in the TolC protein. Two putative gate regions, located at either end of the protein, can be distinguished. These regions are the extracellular loops and the mouth of the periplasmic domain, respectively. The periplasmic gate has been implicated in the conformational changes leading from the closed x-ray structure to a proposed open state of TolC. Between the two gates, a peristaltic motion of the periplasmic domain is observed, which may facilitate transport of the solutes from one end of the tunnel to the other. The motions observed in the atomistic simulations are also seen in coarse-grained simulations in which the protein tertiary structure is represented by an elastic network model.

Coarse-grained simulation: a high-throughput computational approach to membrane proteins

An understanding of the interactions of membrane proteins with a lipid bilayer environment is central to relating their structure to their function and stability. A high-throughput approach to prediction of membrane protein interactions with a lipid bilayer based on coarse-grained Molecular Dynamics simulations is described. This method has been used to develop a database of CG simulations (coarse-grained simulations) of membrane proteins (http://sbcb.bioch.ox.ac.uk/cgdb). Comparison of CG simulations and AT simulations (atomistic simulations) of lactose permease reveals good agreement between the two methods in terms of predicted lipid headgroup contacts. Both CG and AT simulations predict considerable local bilayer deformation by the voltage sensor domain of the potassium channel KvAP.

Amino acids that confer transport of raffinose and maltose sugars in the raffinose permease (RafB) of Escherichia coli as implicated by spontaneous mutations at Val-35, Ser-138, Ser-139, Gly-389 and Ile-391

In order to identify amino acid residues in the Escherichia coli raffinose-H(+) permease (RafB) that play a role in sugar selection and transport, we first incubated E. coli HS4006 containing plasmid pRU600 (expresses inducible raffinose permease and alpha-galactosidase) on maltose MacConkey indicator plates overnight. Initially, all colonies were white, indicating no fermentation of maltose. Upon further incubation, 100 mutants appeared red. pRU600 DNA was prepared from 55 mutants. Five mutants transferred the phenotype for fermentation of maltose (red). Plasmid DNA from five maltose-positive phenotype transformants was prepared and sequenced, revealing three distinct types of mutations. Two mutants exhibited Val-35-->Ala (MT1); one mutant had Ile-391-->Ser (MT2); and two mutants had Ser-138-->Asp, Ser-139-->Leu and Gly-389-->Ala (MT3). Transport studies of [(3)H]-maltose showed that cells harboring MT1, MT2 and MT3 had greater uptake (P <or= 0.05) than cells harboring wild-type RafB. However, [(14)C]-raffinose uptake was reduced in all mutant cells (P <or= 0.05) with MT1, MT2 and MT3 mutants compared to cells harboring wild-type RafB. Kinetic analysis showed enhanced apparent K (m) values for maltose and reduced V (max)/ K (m) ratios for raffinose compared to wild-type values. The apparent K (i) value of maltose for RafB indicates a competitive relationship between maltose and raffinose. Maltose "uphill" accumulation was greater for mutants (P <or= 0.05) than for cells with wild-type RafB. Thus, we implicate residues in RafB that are responsible for raffinose transport and suggest that the substituted residues in RafB dictate structures that enhance transport of maltose.

Substrate binding and formation of an occluded state in the leucine transporter

Translocation through the extracellular vestibule and binding of leucine in the leucine transporter (LeuT) have been studied with molecular dynamics simulations. More than 0.1 mus of all-atom molecular dynamics simulations have been performed on different combinations of LeuT, bound substrate, and bound structural Na(+) ions to describe molecular events involved in substrate binding and in the formation of the occluded state and to investigate the dynamics of this state. Three structural features are found to be directly involved in the initial steps of leucine transport: a Na(+) ion directly coordinated to leucine (Na-1), two aromatic residues closing the binding site toward the extracellular vestibule (Tyr-108 and Phe-253), and a salt bridge in the extracellular vestibule (Arg-30 and Asp-404). These features account for observed differences between simulations of LeuT with and without bound substrate and for a possible pathway for leucine binding and thereby formation of the occluded LeuT binding site.

Sugar binding induces an outward facing conformation of LacY

According to x-ray structure, the lactose permease (LacY) is a monomer organized into N- and C-terminal six-helix bundles that form a deep internal cavity open on the cytoplasmic side with a single sugar-binding site at the apex. The periplasmic side of the molecule is closed. During sugar/H(+) symport, a cavity facing the periplasmic side is thought to open with closure of the inward-facing cytoplasmic cavity so that the sugar-binding site is alternately accessible to either face of the membrane. Double electron-electron resonance (DEER) is used here to measure interhelical distance changes induced by sugar binding to LacY. Nitroxide-labeled paired-Cys replacements were constructed at the ends of transmembrane helices on the cytoplasmic or periplasmic sides of wild-type LacY and in the conformationally restricted mutant Cys-154-->Gly. Distances were then determined in the presence of galactosidic or nongalactosidic sugars. Strikingly, specific binding causes conformational rearrangement on both sides of the molecule. On the cytoplasmic side, each of six nitroxide-labeled pairs exhibits decreased interspin distances ranging from 4 to 21 A. Conversely, on the periplasmic side, each of three spin-labeled pairs shows increased distances ranging from 4 to 14 A. Thus, the inward-facing cytoplasmic cavity closes, and a cleft opens on the tightly packed periplasmic side. In the Cys-154-->Gly mutant, sugar-induced closing is observed on the cytoplasmic face, but little or no change occurs on periplasmic side. The DEER measurements in conjunction with molecular modeling based on the x-ray structure provide strong support for the alternative access model and reveal a structure for the outward-facing conformer of LacY.

Transport cycle intermediate in small multidrug resistance protein is revealed by substrate fluorescence

Efflux pumps of the small multidrug resistance family bind cationic, lipophilic antibiotics and transport them across the membrane in exchange for protons. The transport cycle must involve various conformational states of the protein needed for substrate binding, translocation, and release. A fluorescent substrate will therefore experience a significant change of environment while being transported, which influences its fluorescence properties. Thus the substrate itself can report intermediate states that form during the transport cycle. We show the existence of such a substrate-transporter complex for the EmrE homolog Mycobacterium tuberculosis TBsmr and its substrate ethidium bromide. The pH gradient needed for antiport has been generated by co-reconstituting TBsmr with bacteriorhodopsin. Sample illumination generates a DeltapH, which results in enhanced ethidium fluorescence intensity, which is abolished when DeltapH or DeltaPsi is collapsed or when the essential residue Glu-13 in TBsmr is exchanged with Ala. This observation shows the formation of a pH-dependent, transient substrate-protein complex between binding and release of ethidium. We have further characterized this state by determining the K(d), by inhibiting ethidium transport through titration with nonfluorescent substrate and by fluorescence anisotropy measurements. Our findings support a model with a single occluded intermediate state in which the substrate is highly immobile.

Secondary and tertiary structure elasticity of titin Z1Z2 and a titin chain model

The giant protein titin, which is responsible for passive elasticity in muscle fibers, is built from approximately 300 regular immunoglobulin-like (Ig) domains and FN-III repeats. While the soft elasticity derived from its entropic regions, as well as the stiff mechanical resistance derived from the unfolding of the secondary structure elements of Ig- and FN-III domains have been studied extensively, less is known about the mechanical elasticity stemming from the orientation of neighboring domains relative to each other. Here we address the dynamics and energetics of interdomain arrangement of two adjacent Ig-domains of titin, Z1, and Z2, using molecular dynamics (MD) simulations. The simulations reveal conformational flexibility, due to the domain-domain geometry, that lends an intermediate force elasticity to titin. We employ adaptive biasing force MD simulations to calculate the energy required to bend the Z1Z2 tandem open to identify energetically feasible interdomain arrangements of the Z1 and Z2 domains. The finding is cast into a stochastic model for Z1Z2 interdomain elasticity that is generalized to a multiple domain chain replicating many Z1Z2-like units and representing a long titin segment. The elastic properties of this chain suggest that titin derives so-called tertiary structure elasticity from bending and twisting of its domains. Finally, we employ steered molecular dynamics simulations to stretch individual Z1 and Z2 domains and characterize the so-called secondary structure elasticity of the two domains. Our study suggests that titin's overall elastic response at weak force stems from a soft entropic spring behavior (not described here), from tertiary structure elasticity with an elastic spring constant of approximately 0.001-1 pN/A and, at strong forces, from secondary structure elasticity.

Docking and homology modeling explain inhibition of the human vesicular glutamate transporters

As membrane transporter proteins, VGLUT1-3 mediate the uptake of glutamate into synaptic vesicles at presynaptic nerve terminals of excitatory neural cells. This function is crucial for exocytosis and the role of glutamate as the major excitatory neurotransmitter in the central nervous system. The three transporters, sharing 76% amino acid sequence identity in humans, are highly homologous but differ in regional expression in the brain. Although little is known regarding their three-dimensional structures, hydropathy analysis on these proteins predicts 12 transmembrane segments connected by loops, a topology similar to other members in the major facilitator superfamily, where VGLUT1-3 have been phylogenetically classified. In this work, we present a three-dimensional model for the human VGLUT1 protein based on its distant bacterial homolog in the same superfamily, the glycerol-3-phosphate transporter from Escherichia coli. This structural model, stable during molecular dynamics simulations in phospholipid bilayers solvated by water, reveals amino acid residues that face its pore and are likely to affect substrate translocation. Docking of VGLUT1 substrates to this pore localizes two different binding sites, to which inhibitors also bind with an overall trend in binding affinity that is in agreement with previously published experimental data.

Conformational Change in an MFS Protein: MD Simulations of LacY

Molecular dynamics simulations of lactose permease (LacY) in a phospholipid bilayer reveal the conformational dynamics of the protein. In inhibitor-bound simulations (i.e., those closest to the X-ray structure) the protein was stable, showing little conformational change over a 50 ns timescale. Movement of the bound inhibitor, TDG, to an alternative binding mode was observed, so that it interacted predominantly with the N-terminal domain and with residue E269 from the C-terminal domain. In multiple ligand-free simulations, a degree of domain closure occurred. This switched LacY to a state with a central cavity closed at both the intracellular and periplasmic ends. This may resemble a possible intermediate in the transport mechanism. Domain closure occurs by a combination of rigid-body movements of domains and of intradomain motions of helices, especially TM4, TM5, TM10, and TM11. A degree of intrahelix flexibility appears to be important in the conformational change.

OmpA: gating and dynamics via molecular dynamics simulations

Outer membrane proteins (OMPs) of Gram-negative bacteria have a variety of functions including passive transport, active transport, catalysis, pathogenesis and signal transduction. Whilst the structures of approximately 25 OMPs are currently known, there is relatively little known about their dynamics in different environments. The outer membrane protein, OmpA from Escherichia coli has been studied extensively in different environments both experimentally and computationally, and thus provides an ideal test case for the study of the dynamics and environmental interactions of outer membrane proteins. We review molecular dynamics simulations of OmpA and its homologues in a variety of different environments and discuss possible mechanisms of pore gating. The transmembrane domain of E. coli OmpA shows subtle differences in dynamics and interactions between a detergent micelle and a lipid bilayer environment. Simulations of the crystallographic unit cell reveal a micelle-like network of detergent molecules interacting with the protein monomers. Simulation and modelling studies emphasise the role of an electrostatic-switch mechanism in the pore-gating mechanism. Simulation studies have been extended to comparative models of OmpA homologues from Pseudomonas aeruginosa (OprF) and Pasteurella multocida (PmOmpA), the latter model including the periplasmic C-terminal domain.

Sugar transport across lactose permease probed by steered molecular dynamics

Escherichia coli lactose permease (LacY) transports sugar across the inner membrane of the bacterium using the proton motive force to accumulate sugar in the cytosol. We have probed lactose conduction across LacY using steered molecular dynamics, permitting us to follow molecular and energetic details of lactose interaction with the lumen of LacY during its permeation. Lactose induces a widening of the narrowest parts of the channel during permeation, the widening being largest within the periplasmic half-channel. During permeation, the water-filled lumen of LacY only partially hydrates lactose, forcing it to interact with channel lining residues. Lactose forms a multitude of direct sugar-channel hydrogen bonds, predominantly with residues of the flexible N-domain, which is known to contribute a major part of LacY's affinity for lactose. In the periplasmic half-channel lactose predominantly interacts with hydrophobic channel lining residues, whereas in the cytoplasmic half-channel key protein-substrate interactions are mediated by ionic residues. A major energy barrier against transport is found within a tight segment of the periplasmic half-channel where sugar hydration is minimal and protein-sugar interaction maximal. Upon unbinding from the binding pocket, lactose undergoes a rotation to permeate either half-channel with its long axis aligned parallel to the channel axis. The results hint at the possibility of a transport mechanism, in which lactose permeates LacY through a narrow periplasmic half-channel and a wide cytoplasmic half-channel, the opening of which is controlled by changes in protonation states of key protein side groups.

Lactose permease H+-lactose symporter: mechanical switch or Brownian ratchet?

Lactose permease structure is deemed consistent with a mechanical switch device for H(+)-coupled symport. Because the crystallography-assigned docking position of thiodigalactoside (TDG) does not make close contact with several amino acids essential for symport; the switch model requires allosteric interactions between the proton and sugar binding sites. The docking program, Autodock 3 reveals other lactose-docking sites. An alternative cotransport mechanism is proposed where His-322 imidazolium, positioned in the central pore equidistant (5-7 A) between six charged amino acids, Arg-302 and Lys-319 opposing Glu-269, Glu-325, Asp-237, and Asp-240, transfers a proton transiently to an H-bonded lactose hydroxyl group. Protonated lactose and its dissociation product H(3)O+ are repelled by reprotonated His-322 and drift in the electrostatic field toward the cytosol. This Brownian ratchet model, unlike the conventional carrier model, accounts for diminished symport by H322N mutant; how H322 mutants become uniporters; why exchanging Lys-319 with Asp-240 paradoxically inactivates symport; how some multiple mutants become revertant transporters; the raised export rate and affinity toward lactose of uncoupled mutants; the altered specificity toward lactose, melibiose, and galactose of some mutants, and the proton dissociation rate of H322 being 100-fold faster than the symport turnover rate.

Sugar binding in lactose permease: anomeric state of a disaccharide influences binding structure

Lactose permease in Escherichia coli (LacY) transports both anomeric states of disaccharides but has greater affinity for alpha-sugars. Molecular dynamics (MD) simulations are used to probe the protein-sugar interactions, binding structures, and global protein motions in response to sugar binding by investigating LacY (the experimental mutant and wild-type) embedded in a fully hydrated lipid bilayer. A total of 12 MD simulations of 20-25 ns each with beta(alpha)-d-galactopyranosyl-(1,1)-beta-d-galactopyranoside (betabeta-(Galp)(2)) and alphabeta-(Galp)(2) result in binding conformational families that depend on the anomeric state of the sugar. Both sugars strongly interact with Glu126 and alphabeta-(Galp)(2) has a greater affinity to this residue. Binding conformations are also seen that involve protein residues not observed in the crystal structure, as well as those involved in the proton translocation (Phe118, Asn119, Asn240, His322, Glu325, and Tyr350). Common to nearly all protein-sugar structures, water acts as a hydrogen bond bridge between the disaccharide and protein. The average binding energy is more attractive for alphabeta-(Galp)(2) than betabeta-(Galp)(2), i.e. -10.7(+/-0.7) and -3.1(+/-1.0) kcal/mol, respectively. Of the 12 helices in LacY, helix-IV is the least stable with betabeta-(Galp)(2) binding resulting in larger distortion than alphabeta-(Galp)(2).