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.     

Modeling of glycerol-3-phosphate transporter suggests a potential 'tilt' mechanism involved in its function

Many major facilitator superfamily (MFS) transporters have similar 12-transmembrane alpha-helical topologies with two six-helix halves connected by a long loop. In humans, these transporters participate in key physiological processes and are also, as in the case of members of the organic anion transporter (OAT) family, of pharmaceutical interest. Recently, crystal structures of two bacterial representatives of the MFS family--the glycerol-3-phosphate transporter (GlpT) and lac-permease (LacY)--have been solved and, because of assumptions regarding the high structural conservation of this family, there is hope that the results can be applied to mammalian transporters as well. Based on crystallography, it has been suggested that a major conformational "switching" mechanism accounts for ligand transport by MFS proteins. This conformational switch would then allow periodic changes in the overall transporter configuration, resulting in its cyclic opening to the periplasm or cytoplasm. Following this lead, we have modeled a possible "switch" mechanism in GlpT, using the concept of rotation of protein domains as in the DynDom program17 and membranephilic constraints predicted by the MAPAS program.(23) We found that the minima of energies of intersubunit interactions support two alternate positions consistent with their transport properties. Thus, for GlpT, a "tilt" of 9 degrees -10 degrees rotation had the most favorable energetics of electrostatic interaction between the two halves of the transporter; moreover, this confirmation was sufficient to suggest transport of the ligand across the membrane. We conducted steered molecular dynamics simulations of the GlpT-ligand system to explore how glycerol-3-phosphate would be handled by the "tilted" structure, and obtained results generally consistent with experimental mutagenesis data. While biochemical data remain most consistent with a single-site alternating access model, our results raise the possibility that, while the "rocker switch" may apply to certain MFS transporters, intermediate "tilted" states may exist under certain circumstances or as transitional structures. Although wet lab experimental confirmation is required, our results suggest that transport mechanisms in this transporter family should probably not be assumed to be conserved simply based on standard structural homology considerations. Furthermore, steered molecular dynamics elucidating energetic interactions of ligands with amino acid residues in an appropriately modeled transporter may have predictive value in understanding the impact of mutations and/or polymorphisms on transporter function.     

Lactose permease as a paradigm for membrane transport proteins (Review)

Our structural knowledge of the major facilitator superfamily (MFS) has dramatically increased in the past year with three structures of proteins from the MFS (oxalate/formate antiporter; lactose/proton symporter and the P_i/glycerol-3-phosphate antiporter). All three structures revealed 12 transmembrane helices forming two distinct domains and could imply that members of the MFS have preserved both secondary as well as tertiary structural elements during evolution. Lactose permease, a particularly well-studied member of the MFS, has been extensively explored by a number of molecular biological, biochemical and biophysical approaches. In this review, we take a closer look at the structure of LacY and incorporate a wealth of biochemical and biophysical data in order to propose a possible mechanism for lactose/proton symport. In addition, we make some brief comparisons between the structures of LacY and GlpT.     

Lactose permease as a paradigm for membrane transport proteins (Review)

Our structural knowledge of the major facilitator superfamily (MFS) has dramatically increased in the past year with three structures of proteins from the MFS (oxalate/formate antiporter; lactose/proton symporter and the P(i)/glycerol-3-phosphate antiporter). All three structures revealed 12 transmembrane helices forming two distinct domains and could imply that members of the MFS have preserved both secondary as well as tertiary structural elements during evolution. Lactose permease, a particularly well-studied member of the MFS, has been extensively explored by a number of molecular biological, biochemical and biophysical approaches. In this review, we take a closer look at the structure of LacY and incorporate a wealth of biochemical and biophysical data in order to propose a possible mechanism for lactose/proton symport. In addition, we make some brief comparisons between the structures of LacY and GlpT.     

3D model of the Escherichia coli multidrug transporter MdfA reveals an essential membrane-embedded positive charge

MdfA is an Escherichia coli multidrug transporter of the major facilitator superfamily (MFS) of secondary transporters. Although several aspects of multidrug recognition by MdfA have been characterized, better understanding the detailed mechanism of its function requires structural information. Previous studies have modeled the 3D structures of MFS proteins, based on the X-ray structure of LacY and GlpT. However, because of poor sequence homology, between LacY, GlpT, and MdfA additional constraints were required for a reliable homology modeling. Using an algorithm that predicts the angular orientation of each transmembrane helix (TM) (kPROT), we obtained a remarkably similar pattern for the 12 TMs of MdfA and those of GlpT and LacY, suggesting that they all have similar helix packing. Consequently, a 3D model was constructed for MdfA by structural alignment with LacY and GlpT, using the kPROT results as an additional constraint. Further refinement and a preliminary evaluation of the model were achieved by correlated mutation analysis and the available experimental data. Surprisingly, in addition to the previously characterized membrane-embedded glutamate at position 26, the model suggests that Asp34 and Arg112 are located within the membrane, on the same face of the cavity as Glu26. Importantly, Arg112 is evolutionarily conserved in secondary drug transporters, and here we show that a positive charge at this position is absolutely essential for multidrug transport by MdfA.     

Structural modeling of dual-affinity purified Pho84 phosphate transporter

The phosphate transporter Pho84 of Saccharomyces cerevisiae is predicted to contain 12 transmembrane (TM) regions, divided into two partially duplicated parts of 6 TM segments. The three-dimensional (3D) organization of the Pho84 protein has not yet been determined. However, the 3D crystal structure of the Escherichia coli MFS glycerol-3-phosphate/phosphate antiporter, GlpT, and lactose transporter, LacY, has recently been determined. On the basis of extensive prediction and fold recognition analyses (at the MetaServer), GlpT was proposed as the best structural template on which the arrangement of TM segments of the Pho84 transporter was fit, using the comparative structural modeling program MODELLER. To initiate an evaluation of the appropriateness of the Pho84 model, we have performed two direct tests by targeting spin labels to putative TM segments 8 and 12. Electron paramagnetic resonance spectroscopy was then applied on purified and spin labeled Pho84. The line shape from labels located at both positions is consistent with the structural environment predicted by the template-generated model, thus supporting the model.     

Proline-induced hinges in transmembrane helices: possible roles in ion channel gating

A number of ion channels contain transmembrane (TM) alpha-helices that contain proline-induced molecular hinges. These TM helices include the channel-forming peptide alamethicin (Alm), the S6 helix from voltage-gated potassium (Kv) channels, and the D5 helix from voltage-gated chloride (CLC) channels. For both Alm and KvS6, experimental data implicate hinge-bending motions of the helix in an aspect of channel gating. We have compared the hinge-bending motions of these TM helices in bilayer-like environments by multi-nanosecond MD simulations in an attempt to describe motions of these helices that may underlie possible modes of channel gating. Alm is an alpha-helical channel-forming peptide, which contains a central kink associated with a Gly-x-x-Pro motif in its sequence. Simulations of Alm in a TM orientation for 10 ns in an octane slab indicate that the Gly-x-x-Pro motif acts as a molecular hinge. The S6 helix from Shaker Kv channels contains a Pro-Val-Pro motif. Modeling studies and recent experimental data suggest that the KvS6 helix may be kinked in the vicinity of this motif. Simulations (10 ns) of an isolated KvS6 helix in an octane slab and in a POPC bilayer reveal hinge-bending motions. A pattern-matching approach was used to search for possible hinge-bending motifs in the TM helices of other ion channel proteins. This uncovered a conserved Gly-x-Pro motif in TM helix D5 of CLC channels. MD simulations of a model of hCLC1-D5 spanning an octane slab suggest that this channel also contains a TM helix that undergoes hinge-bending motion. In conclusion, our simulations suggest a model in which hinge-bending motions of TM helices may play a functional role in the gating mechanisms of several different families of ion channels.     

Salt-bridge dynamics control substrate-induced conformational change in the membrane transporter GlpT

Active transport of substrates across cytoplasmic membranes is of great physiological, medical and pharmaceutical importance. The glycerol-3-phosphate (G3P) transporter (GlpT) of the E. coli inner membrane is a secondary active antiporter from the ubiquitous major facilitator superfamily that couples the import of G3P to the efflux of inorganic phosphate (P_i) down its concentration gradient. Integrating information from a novel combination of structural, molecular dynamics simulations and biochemical studies, we identify the residues involved directly in binding of substrate to the inward-facing conformation of GlpT, thus defining the structural basis for the substrate-specificity of this transporter. The substrate binding mechanism involves protonation of a histidine residue at the binding site. Furthermore, our data suggest that the formation and breaking of inter- and intradomain salt bridges control the conformational change of the transporter that accompanies substrate translocation across the membrane. The mechanism we propose may be a paradigm for organophosphate:phosphate antiporters.     

Structural Basis of Substrate Selectivity in the Glycerol-3-Phosphate: Phosphate Antiporter GlpT

Major facilitators represent the largest superfamily of secondary active transporter proteins and catalyze the transport of an enormous variety of small solute molecules across biological membranes. However, individual superfamily members, although they may be architecturally similar, exhibit strict specificity toward the substrates they transport. The structural basis of this specificity is poorly understood. A member of the major facilitator superfamily is the glycerol-3-phosphate (G3P) transporter (GlpT) from the Escherichia coli inner membrane. GlpT is an antiporter that transports G3P into the cell in exchange for inorganic phosphate (P_i). By combining large-scale molecular-dynamics simulations, mutagenesis, substrate-binding affinity, and transport activity assays on GlpT, we were able to identify key amino acid residues that confer substrate specificity upon this protein. Our studies suggest that only a few amino acid residues that line the transporter lumen act as specificity determinants. Whereas R45, K80, H165, and, to a lesser extent Y38, Y42, and Y76 contribute to recognition of both free P_i and the phosphate moiety of G3P, the residues N162, Y266, and Y393 function in recognition of only the glycerol moiety of G3P. It is the latter interactions that give the transporter a higher affinity to G3P over P_i.     

Eukaryotic major facilitator superfamily transporter modeling based on the prokaryotic GlpT crystal structure (Review)

The major facilitator superfamily (MFS) of transporters represents the largest family of secondary active transporters and has a diverse range of substrates. With structural information for four MFS transporters, we can see a strong structural commonality suggesting, as predicted, a common architecture for MFS transporters. The rate for crystal structure determination of MFS transporters is slow, making modeling of both prokaryotic and eukaryotic transporters more enticing. In this review, models of eukaryotic transporters Glut1, G6PT, OCT1, OCT2 and Pho84, based on the crystal structures of the prokaryotic GlpT, based on the crystal structure of LacY are discussed. The techniques used to generate the different models are compared. In addition, the validity of these models and the strategy of using prokaryotic crystal structures to model eukaryotic proteins are discussed. For comparison, E. coli GlpT was modeled based on the E. coli LacY structure and compared to the crystal structure of GlpT demonstrating that experimental evidence is essential for accurate modeling of membrane proteins.     

A structural model for the osmosensor, transporter, and osmoregulator ProP of Escherichia coli

Transporter ProP of Escherichia coli, a member of the major facilitator superfamily (MFS), acts as an osmosensor and an osmoregulator in cells and after purification and reconstitution in proteoliposomes. H(+)-osmoprotectant symport via ProP is activated when medium osmolality is elevated with membrane impermeant osmolytes. The three-dimensional structure of ProP was modeled with the crystal structure of MFS member GlpT as a template. This GlpT structure represents the inward (or cytoplasm)-facing conformation predicted by the alternating access model for transport. LacZ-PhoA fusion analysis and site-directed fluorescence labeling substantiated the membrane topology and orientation predicted by this model and most hydropathy analyses. The model predicts the presence of a proton pathway within the N-terminal six-helix bundle of ProP (as opposed to the corresponding pathway found within the C-terminal helix bundle of its paralogue, LacY). Replacement of residues within the N-terminal helix bundle impaired the osmotic activation of ProP, providing the first indication that residues outside the C-terminal domain are involved in osmosensing. Some residues that were accessible from the periplasmic side, as predicted by the structural model, were more susceptible to covalent labeling in permeabilized membrane fractions than in intact bacteria. These residues may be accessible from the cytoplasmic side in structures not represented by our current model, or their limited exposure in vivo may reflect constraints on transporter structure that are related to its osmosensory mechanism.     

Functionally important amino acids in rice sucrose transporter OsSUT1

Six conserved, charged amino acids within membrane spans in rice sucrose transporter OsSUT1 were identified using a three-dimensional structural model based on the crystal structures of three major facilitator superfamily (MFS) proteins: LacY, GlpT, and EmrD. These positions in OsSUT1 were selected for mutagenesis and biochemical assays. Among the six mutants, D177N completely lost transport function, D331N retained only a small fraction of sucrose uptake activity (2.3% of that of the wild type), and R335H and E336Q also displayed a substantial decrease in transport activity. D329N functioned as well as wild-type OsSUT1. R188K did not transport sucrose but showed a H(+) leak that was inhibited by sucrose, indicating that R188K had uncoupled sucrose and H(+) translocation. This demonstrates that charged amino acids within membrane spans are important for the transport mechanism of OsSUT1 as they are in lactose permease.     

Identification of proline residues in the core cytoplasmic and transmembrane regions of multidrug resistance protein 1 (MRP1/ABCC1) important for transport function, substrate specificity, and nucleotide interactions

Multidrug resistance protein 1 (MRP1/ABCC1) is an ATP-binding cassette transporter that confers resistance to drugs and mediates the transport of organic anions. MRP1 has a core structure of two membrane spanning domains (MSDs) each followed by a nucleotide binding domain. This core structure is preceded by a third MSD with five transmembrane (TM) helices, whereas MSD2 and MSD3 each contain six TM helices. We investigated the consequences of Ala substitution of 18 Pro residues in both the non-membrane and TM regions of MSD2 and MSD3 on MRP1 expression and organic anion transport function. All MRP1-Pro mutants except P1113A were expressed in human embryonic kidney cells at levels comparable with wild-type MRP1. In addition, five mutants containing substitutions of Pro residues in or proximal to the TM helices of MSD2 (TM6-Pro(343), TM8-Pro(448), TM10-Pro(557), and TM11-Pro(595)) and MSD3 (TM14-Pro(1088)) exhibited significantly reduced transport of five organic anion substrates. In contrast, mutation of Pro(1150) in the cytoplasmic loop (CL7) linking TM15 to TM16 caused a substantial increase in 17beta-estradiol-17-beta-(D-glucuronide) and methotrexate transport, whereas transport of other organic anions was reduced or unchanged. Significant substrate-specific changes in the ATP dependence of transport and binding by the P1150A mutant were also observed. Our findings demonstrate the importance of TM6, TM8, TM10, TM11, and TM14 in MRP1 transport function and suggest that CL7 may play a differential role in coupling the activity of the nucleotide binding domains to the translocation of different substrates across the membrane.     

Structure-Function Studies of the SLC17 Transporter Sialin Identify Crucial Residues and Substrate-induced Conformational Changes

Salla disease and infantile sialic acid storage disorder are human diseases caused by loss of function of sialin, a lysosomal transporter that mediates H⁺-coupled symport of acidic sugars N-acetylneuraminic acid and glucuronic acid out of lysosomes. Along with the closely related vesicular glutamate transporters, sialin belongs to the SLC17 transporter family. Despite their critical role in health and disease, these proteins remain poorly understood both structurally and mechanistically. Here, we use substituted cysteine accessibility screening and radiotracer flux assays to evaluate experimentally a computationally generated three-dimensional structure model of sialin. According to this model, sialin consists of 12 transmembrane helices (TMs) with an overall architecture similar to that of the distantly related glycerol 3-phosphate transporter GlpT. We show that TM4 in sialin lines a large aqueous cavity that forms a part of the substrate permeation pathway and demonstrate substrate-induced alterations in accessibility of substituted cysteine residues in TM4. In addition, we demonstrate that one mutant, F179C, has a dramatically different effect on the apparent affinity and transport rate for N-acetylneuraminic acid and glucuronic acid, suggesting that it may be directly involved in substrate recognition and/or translocation. These findings offer a basis for further defining the transport mechanism of sialin and other SLC17 family members.     

Mechanics of channel gating of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is a key molecule involved in the propagation of signals in the central nervous system and peripheral synapses. Although numerous computational and experimental studies have been performed on this receptor, the structural dynamics of the receptor underlying the gating mechanism is still unclear. To address the mechanical fundamentals of nAChR gating, both conventional molecular dynamics (CMD) and steered rotation molecular dynamics (SRMD) simulations have been conducted on the cryo-electron microscopy (cryo-EM) structure of nAChR embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer and water molecules. A 30-ns CMD simulation revealed a collective motion amongst C-loops, M1, and M2 helices. The inward movement of C-loops accompanying the shrinking of acetylcholine (ACh) binding pockets induced an inward and upward motion of the outer β-sheet composed of β9 and β10 strands, which in turn causes M1 and M2 to undergo anticlockwise motions around the pore axis. Rotational motion of the entire receptor around the pore axis and twisting motions among extracellular (EC), transmembrane (TM), and intracellular MA domains were also detected by the CMD simulation. Moreover, M2 helices undergo a local twisting motion synthesized by their bending vibration and rotation. The hinge of either twisting motion or bending vibration is located at the middle of M2, possibly the gate of the receptor. A complementary twisting-to-open motion throughout the receptor was detected by a normal mode analysis (NMA). To mimic the pulsive action of ACh binding, nonequilibrium MD simulations were performed by using the SRMD method developed in one of our laboratories. The result confirmed all the motions derived from the CMD simulation and NMA. In addition, the SRMD simulation indicated that the channel may undergo an open-close (O ↔ C) motion. The present MD simulations explore the structural dynamics of the receptor under its gating process and provide a new insight into the gating mechanism of nAChR at the atomic level.     

Single replacement constructs of all hydroxyl, basic, and acidic amino acids identify new function and structure-sensitive regions of the mitochondrial phosphate transport protein

The phosphate transport protein (PTP) catalyzes the proton cotransport of phosphate into the mitochondrial matrix. It functions as a homodimer, and thus residues of the phosphate and proton pores are somewhat scattered throughout the primary sequence. With 71 new single mutation per subunit PTPs, all its hydroxyl, basic, and acidic residues have now been replaced to identify these essential residues. We assayed the initial rate of pH gradient-dependent unidirectional phosphate transport activity and the liposome incorporation efficiency (LIE) of these mutants. Single mutations of Thr79, Tyr83, Lys90, Tyr94, and Lys98 inactivate transport. The spacings between these residues imply that they are located along the same face of transmembrane (TM) helix B, requiring an extension of its current model C-terminal domain by 10 residues. This extension superimposes very well onto the shorter bovine PTP helix B, leaving a 15-residue hydrophobic extension of the yeast helix B N-terminus. This is similar to the helix D and F regions of the yeast PTP. Only one transport-inhibiting mutation is located within loops: Ser158Thr in the matrix loop between helices C and D. All other transport-inhibiting mutations are located within the TM helices. Mutations that yield LIEs of <6% are all, except for four, within helices. The four exceptions are Tyr12Ala near the PTP N-terminus and Arg159Ala, Glu163Gln, and Glu164Gln in the loop between helices C and D. The PTP C-terminal segment beyond Thr214 at the N-terminus of helix E has 11 mutations with LIEs >20% and none with LIE <6%. Mutations with LIEs >20% are located near the ends of all the TM helices except TM helix D. Only a few mutations alter PTP structure (LIE) and also affect PTP transport activity. A novel observation is that Ser4Ala blocks the formation of PTP bacterial inclusion bodies.     

Mitochondrial phosphate transport protein. Reversions of inhibitory conservative mutations identify four helices and a nonhelix protein segment with transmembrane interactions and Asp39, Glu137, and Ser158 as nonessential for transport

The mitochondrial phosphate transport protein (PTP) has six (A--F) transmembrane (TM) helices per subunit of functional homodimer with all mutations referring to the subunit of the homodimer. In earlier studies, conservative replacements of several residues located either at the matrix end (Asp39/helix A, Glu137/helix C, Asp236/helix E) or at the membrane center (His32/helix A, Glu136/helix C) of TM helices yielded inactive single mutation PTPs. Some of these residues were suggested to act as phosphate ligands or as part of the proton cotransport path. We now show that the mutation Ser158Thr, not part of a TM helix but located near the center of the matrix loop (Ile141--Ser171) between TM helices C and D, inactivates PTP and is thus also functionally relevant. On the other side of the membrane, the single mutation Glu192Asp at the intermembrane space end of TM helix D yields a PTP with 33% wild-type activity. We constructed double mutants by adding this mutation to the six transport-inactivating mutations. Transport was detected only in those with Asp39Asn, Glu137Gln, or Ser158Thr. We conclude that TM helix D can interact with TM helices A and C and matrix loop Ile141--Ser171 and that Asp39, Glu137, and Ser158 are not essential for phosphate transport. Since our results are consistent with residues present in all 12 functionally identified members of the mitochondrial transport protein (MTP) family, they lead to a general rule that specifies MTP residue types at 7 separate locations. The conformations of all the double mutation PTPs (except that with the matrix loop Ser158Thr) are significantly different from those of the single mutation PTPs, as indicated by their very low liposome incorporation efficiency and their requirement for less detergent (Triton X-100) to stay in solution. These dramatic conformational differences also suggest an interaction between TM helices D and E. The results are discussed in terms of TM helix movements and changes in the PTP monomer/dimer ratio.     

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of α-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 × 5 μs) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2× rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of ∼16°. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.     

Transmembrane helices can induce domain formation in crowded model membranes

We studied compositionally heterogeneous multi-component model membranes comprised of saturated lipids, unsaturated lipids, cholesterol, and α-helical TM protein models using coarse-grained molecular dynamics simulations. Reducing the mismatch between the length of the saturated and unsaturated lipid tails reduced the driving force for segregation into liquid-ordered (l(o)) and liquid-disordered (l(d)) lipid domains. Cholesterol depletion had a similar effect, and binary lipid mixtures without cholesterol did not undergo large-scale phase separation under the simulation conditions. The phase-separating ternary dipalmitoyl-phosphatidylcholine (DPPC)/dilinoleoyl-PC (DLiPC)/cholesterol bilayer was found to segregate into l(o) and l(d) domains also in the presence of a high concentration of ΤΜ helices. The l(d) domain was highly crowded with TM helices (protein-to-lipid ratio ~1:5), slowing down lateral diffusion by a factor of 5-10 as compared to the dilute case, with anomalous (sub)-diffusion on the μs time scale. The membrane with the less strongly unsaturated palmitoyl-linoleoyl-PC instead of DLiPC, which in the absence of TM α-helices less strongly deviated from ideal mixing, could be brought closer to a miscibility critical point by introducing a high concentration of TM helices. Finally, the 7-TM protein bacteriorhodopsin was found to partition into the l(d) domains irrespective of hydrophobic matching. These results show that it is possible to directly study the lateral reorganization of lipids and proteins in compositionally heterogeneous and crowded model biomembranes with coarse-grained molecular dynamics simulations, a step toward simulations of realistic, compositionally complex cellular membranes. This article is part of a Special Issue entitled: Protein Folding in Membranes.