The role of leucine and isoleucine in tuning the hydropathy of class A GPCRs

Leucine and Isoleucine are two amino acids that differ only by the positioning of one methyl group. This small difference can have important consequences in α-helices, as the β-branching of Ile results in helix destabilization. We set out to investigate whether there are general trends for the occurrences of Leu and Ile residues in the structures and sequences of class A GPCRs (G protein-coupled receptors). GPCRs are integral membrane proteins in which α-helices span the plasma membrane seven times and which play a crucial role in signal transmission. We found that Leu side chains are generally more exposed at the protein surface than Ile side chains. We explored whether this difference might be attributed to different functions of the two amino acids and tested if Leu tunes the hydrophobicity of the transmembrane domain based on the Wimley-White whole-residue hydrophobicity scales. Leu content decreases the variation in hydropathy between receptors and correlates with the non-Leu receptor hydropathy. Both measures indicate that hydropathy is tuned by Leu. To test this idea further, we generated protein sequences with random amino acid compositions using a simple numerical model, in which hydropathy was tuned by adjusting the number of Leu residues. The model was able to replicate the observations made with class A GPCR sequences. We speculate that the hydropathy of transmembrane domains of class A GPCRs is tuned by Leu (and to some lesser degree by Lys and Val) to facilitate correct insertion into membranes and/or to stably anchor the receptors within membranes.     

Energetics, stability, and prediction of transmembrane helices

We show that the peptide backbone of an alpha-helix places a severe thermodynamic constraint on transmembrane (TM) stability. Neglect of this constraint by commonly used hydrophobicity scales underlies the notorious uncertainty of TM helix prediction by sliding-window hydropathy plots of membrane protein (MP) amino acid sequences. We find that an experiment-based whole-residue hydropathy scale (WW scale), which includes the backbone constraint, identifies TM helices of membrane proteins with an accuracy greater than 99 %. Furthermore, it correctly predicts the minimum hydrophobicity required for stable single-helix TM insertion observed in Escherichia coli. In order to improve membrane protein topology prediction further, we introduce the augmented WW (aWW) scale, which accounts for the energetics of salt-bridge formation. An important issue for genomic analysis is the ability of the hydropathy plot method to distinguish membrane from soluble proteins. We find that the method falsely predicts 17 to 43 % of a set of soluble proteins to be MPs, depending upon the hydropathy scale used.     

Membrane insertion and dissociation processes of a model transmembrane helix

An important subject for elucidating membrane protein (MP) folding is how transmembrane helices (TMHs) insert into and dissociate from membranes. We investigated helix dissociation kinetics and insertion topology by means of intervesicular transfer of the fluorophore-labeled completely hydrophobic model transmembrane helix NBD-(LALAAAA)(3)-NH(2) (NBD = 7-nitro-2-1,3-benzoxadiazol-4-yl). The peptide forms a topologically stable transmembrane helix, which is in a monomer-antiparallel dimer equilibrium [Yano, Y., Takemoto, T., Kobayashi, S., Yasui, H., Sakurai, H., Ohashi, W., Niwa, M., Futaki, S., Sugiura, Y., and Matsuzaki, K. (2002) Biochemistry 41, 3073-3080]. The helix transfer kinetics, representing the helix dissociation process, was monitored by fluorescence recovery of the quenched peptide in donor vesicles containing a quencher upon its transfer to acceptor vesicles without the quencher. The transfer kinetics and vesicle concentration dependence demonstrated that the transfer was mediated by monomer in the aqueous phase. Furthermore, the activation enthalpy was estimated to be +17.7 +/- 1.3 kcal mol(-1). Helix insertion topology, detected by chemical quenching of the NBD group in the outer leaflet by dithionite ions, was found to be controlled by transmembrane electric potential-helix macro dipole interaction. On the basis of these observations, a model for the helix insertion/dissociation processes was discussed.     

Membrane topology of epoxide hydrolase

The amino acid sequences of epoxide hydrolase from rat, rabbit and human have been subjected to hydropathy analysis and a novel model for the membrane topology of this enzyme is presented. The enzyme would appear to be retained in microsomal membranes by a single transmembrane segment located at the N-terminus and the majority (96%) of the protein is exposed at the cytosolic membrane surface. This model is significantly different from a scheme suggested by analysis of the rat enzyme alone which proposed six transmembrane domains (Porter et al. (1988) Arch. Biochem. Biophys. 248, 121-129). Experiments with rat microsomal membranes were conducted to distinguish between the two models and used proteolytic enzymes and non-permeant chemical probes. Epoxide hydrolase of intact and permeabilised membranes was resistant to digestion by a number of proteinases. However, this is likely to be related to a compact fold of the protein rather than membrane association since purified, delipidated enzyme preparations were also resistant to proteolysis. While the use of proteinases did not provide useful membrane topological information, experiments with the fluorescent probe, 3-azido-2,7-naphthalenedisulphonate strongly support the view that the majority of the protein is indeed exposed at the cytosolic surface of the membranes. The analysis illustrates the caution which must be employed in the formulation of topological models based on hydropathy plots alone and the value of considering homologous proteins.     

Hydrophobicity and prediction of the secondary structure of membrane proteins and peptides.

Reliability of the hydropathy method to predict the formation of membrane-spanning alpha-helices by integral membrane proteins and peptides whose structure is known from X-ray crystallography is analysed. It is shown that Kyte-Doolittle hydropathy plots do not predict accurately 22 transmembrane alpha-helices in the reaction centres (RC) of the photosynthetic bacteria Rhodopseudomonas viridis and Rhodobacter sphaeroides (R-26). The accuracy of prediction for these proteins was improved using an optimised Kyte-Doolittle hydrophobicity scale. However, this hydrophobicity scale did not improve the predictions for the alphabeta-peptides of the B800-850 (LH2) complexes of the photosynthetic bacteria Rhodopseudomonas acidophila and Rhodospirillum molischianum, which were excluded from the optimisation procedure. The best and worst predictions of membrane-spanning alpha-helices for the RC proteins and LH2 peptides, respectively, were obtained with a propensity scale (PRC) calculated from the amino acid sequences and X-ray data for the RC proteins. A propensity scale (PLH) obtained using the amino acid sequences and X-ray data for the alphabeta-peptides of the LH2 complexes did not give an acceptable prediction of the transmembrane segments in the LH2 peptides; moreover, it markedly contradicted the PRC scale. Amino acids have been concluded to have no significant preference to localisation in transmembrane segments. Therefore, the predictive ability of the hydropathy methodology appears to be limited: the number of transmembrane segments can be correctly calculated for the best case only, and the lengths and positions of membrane-spanning alpha-helices in a protein amino acid sequence can not be predicted exactly.     

The distal cytoplasmic tail of the influenza A M2 protein dynamically extends from the membrane

The influenza A M2 protein is a multifunctional membrane-associated homotetramer that orchestrates several essential events in the viral infection cycle. The monomeric subunits of the M2 homotetramer consist of an N-terminal ectodomain, a transmembrane domain, and a C-terminal cytoplasmic domain. The transmembrane domain forms a four-helix proton channel that promotes uncoating of virions upon host cell entry. The membrane-proximal region of the C-terminal domain forms a surface-associated amphipathic helix necessary for viral budding. The structure of the remaining ~34 residues of the distal cytoplasmic tail has yet to be fully characterized despite the functional significance of this region for influenza infectivity. Here, we extend structural and dynamic studies of the poorly characterized M2 cytoplasmic tail. We used SDSL-EPR to collect site-specific information on the mobility, solvent accessibility, and conformational properties of residues 61–70 of the full-length, cell-expressed M2 protein reconstituted into liposomes. Our analysis is consistent with the predominant population of the C-terminal tail dynamically extending away from the membranes surface into the aqueous medium. These findings provide insight into the hypothesis that the C-terminal domain serves as a sensor that regulates how M2 protein participates in critical events in the viral infection cycle.     

MD simulations of Mistic: conformational stability in detergent micelles and water

Mistic is an unusual membrane protein from Bacillus subtilis. It appears to fold and insert autonomously into a lipid bilayer and has been suggested as a tool that aids the targeting of eukaryotic membrane proteins to bacterial membranes. The NMR structure of Mistic in detergent (LDAO) micelles has revealed it to be a four alpha-helix bundle. From a structural perspective, Mistic does not resemble other membrane proteins. Its external surface is not very hydrophobic, and standard methods do not predict any of its helices to be in the transmembrane orientation. Molecular dynamics simulations (simulation times approximately 30 ns) in water and in detergent micelles have been used to explore the conformational stability of Mistic as a function of its environment. In water, the protein is stable, exhibiting no significant change in fold on a 30 ns time scale. In contrast, in three simulations in detergent micelles, the partial unfolding of Mistic occurred, whereby the H4 helix drifted away from the H1-H3 core. This was due to the penetration of detergent molecules between H4 and the remainder of the protein. This is unlike the behavior of several other membrane proteins, both alpha-helix bundles and beta-barrels, in comparable detergent micelle simulations. The unfolding of H4 from the H1-H3 core of Mistic could be partially reversed by a simulation in which the detergent molecules were removed, and the unfolded protein was simulated in water. These results suggest that Mistic may not be a stable integrated membrane protein but rather that it may undergo a conformational change upon interaction with a membrane or membrane-like environment.     

The determinants of hydrophobic mismatch response for transmembrane helices

Hydrophobic mismatch arises from a difference in the hydrophobic thickness of a lipid membrane and a transmembrane protein segment, and is thought to play an important role in the folding, stability and function of membrane proteins. We have investigated the possible adaptations that lipid bilayers and transmembrane α-helices undergo in response to mismatch, using fully-atomistic molecular dynamics simulations totaling 1.4 μs. We have created 25 different tryptophan-alanine-leucine transmembrane α-helical peptide systems, each composed of a hydrophobic alanine–leucine stretch, flanked by 1–4 tryptophan side chains, as well as the β-helical peptide dimer, gramicidin A. Membrane responses to mismatch include changes in local bilayer thickness and lipid order, varying systematically with peptide length. Adding more flanking tryptophan side chains led to an increase in bilayer thinning for negatively mismatched peptides, though it was also associated with a spreading of the bilayer interface. Peptide tilting, bending and stretching were systematic, with tilting dominating the responses, with values of up to ~ 45° for the most positively mismatched peptides. Peptide responses were modulated by the number of tryptophan side chains due to their anchoring roles and distributions around the helices. Potential of mean force calculations for local membrane thickness changes, helix tilting, bending and stretching revealed that membrane deformation is the least energetically costly of all mismatch responses, except for positively mismatched peptides where helix tilting also contributes substantially. This comparison of energetic driving forces of mismatch responses allows for deeper insight into protein stability and conformational changes in lipid membranes.     

Evidence that insertion of Tomato ringspot nepovirus NTB-VPg protein in endoplasmic reticulum membranes is directed by two domains: a C-terminal transmembrane helix and an N-terminal amphipathic helix

The NTB-VPg protein of Tomato ringspot nepovirus is an integral membrane protein found in association with endoplasmic reticulum (ER)-derived membranes active in virus replication. A transmembrane helix present in a hydrophobic region at the C terminus of the NTB domain was previously shown to traverse the membranes, resulting in the translocation of the VPg domain in the lumen. We have now conducted an in planta analysis of membrane-targeting domains within NTB-VPg using in-frame fusions to the green fluorescent protein (GFP). As expected, the entire NTB-VPg protein directed the GFP fluorescence to ER membranes. GFP fusion proteins containing the C-terminal 86 amino acids of NTB-VPg also associated with ER membranes, resulting in ER-specific glycosylation at a naturally occurring glycosylation site in the VPg domain. Deletion of the hydrophobic region prevented the membrane association. The N-terminal 80 amino acids of NTB were also sufficient to direct the GFP fluorescence to intracellular membranes. A putative amphipathic helix in this region was necessary and sufficient to promote membrane association of the fusion proteins. Using in vitro membrane association assays and glycosylation site mapping, we show that the N terminus of NTB can be translocated in the lumen at least in vitro. This translocation was dependent on the presence of the putative amphipathic helix, suggesting that oligomeric forms of this helix traverse the membrane. Taken together, our results suggest that at least two distinct elements play a key role in the insertion of NTB-VPg in the membranes: a C-terminal transmembrane helix and an N-terminal amphipathic helix. An updated model of the topology of the protein in the membrane is presented.     

Observations concerning topology and locations of helix ends of membrane proteins of known structure

Summary Hydropathy plots of amino acid sequences reveal the approximate locations of the transbilayer helices of membrane proteins of known structure and are thus used to predict the helices of proteins of unknown structure. Because the threedimensional structures of membrane proteins are difficult to obtain, it is important to be able to extract as much information as possible from hydropathy plots. We describe an augmented hydropathy plot analysis of the three membrane proteins of known structure, which should be useful for the systematic examination and comparison of membrane proteins of unknown structure. The sliding-window analysis utilizes the floating interfacial hydrophobicity scale [IFH(h)] of Jacobs and White (Jacobs, R.E., White, S.H., 1989.Biochemistry 28:3421–3437) and the reverse-turn (RT) frequencies of Levitt (Levitt, M., 1977,Biochemistry 17:4277–4285). The IFH(h) scale allows one to examine the consequences of different assumptions about the average hydrogen bond status (h=0 to 1) of polar side chains. Hydrophobicity plots of the three proteins show that (i) the intracellular helix-connecting links and chain ends can be distinguished from the extracellular ones and (ii) the main peaks of hydrophobicity are bounded by minor ones which bracket the helix ends. RT frequency plots show that (iii) the centers of helices are usually very close to wide-window minima of average RT frequency and (iv) helices are always bounded by narrowwindow maxima of average RT frequency. The analysis suggests that side-chain hydrogen bonding with membrane components during folding may play a key role in insertion.     

Interactions of hydrophobic peptides with lipid bilayers: Monte Carlo simulations with M2delta

We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.     

Interaction of transmembrane helices by a knobs-into-holes packing characteristic of soluble coiled coils

Membrane-embedded protein domains frequently exist as α-helical bundles, as exemplified by photosynthetic reaction centers, bacteriorhodopsin, and cytochrome C oxidase. The sidechain packing between their transmembrane helices was investigated by a nearest-neighbor analysis which identified sets of interfacial residues for each analyzed helix–helix interface. For the left-handed helix–helix pairs, the interfacial residues almost exclusively occupy positions a, d, e, or g within a heptad motif (abcdefg) which is repeated two to three times for each interacting helical surface. The connectivity between the interfacial residues of adjacent helices conforms to the knobs-into-holes type of sidechain packing known from soluble coiled coils. These results demonstrate on a quantitative basis that the geometry of sidechain packing is similar for left-handed helix–helix pairs embedded in membranes and coiled coils of soluble proteins. The transmembrane helix–helix interfaces studied are somewhat less compact and regular as compared to soluble coiled coils and tolerate all hydrophobic amino acid types to similar degrees. The results are discussed with respect to previous experimental findings which demonstrate that specific interactions between transmembrane helices are important for membrane protein folding and/or oligomerization. Proteins 31:150–159, 1998. © 1998 Wiley-Liss, Inc.     

NMR-Based Simulation Studies of Pf1 Coat Protein in Explicit Membranes

As time- and ensemble-averaged measures, NMR observables contain information about both protein structure and dynamics. This work represents a computational study to extract such information for membrane proteins from orientation-dependent NMR observables: solid-state NMR chemical shift anisotropy and dipolar coupling, and solution NMR residual dipolar coupling. We have performed NMR-restrained molecular dynamics simulations to refine the structure of the membrane-bound form of Pf1 coat protein in explicit lipid bilayers using the recently measured chemical shift anisotropy, dipolar coupling, and residual dipolar coupling data. From the simulations, we have characterized detailed protein-lipid interactions and explored the dynamics. All simulations are stable and the NMR restraints are well satisfied. The C-terminal transmembrane (TM) domain of Pf1 finds its optimal position in the membrane quickly (within 6 ns), illustrating efficient solvation of TM domains in explicit bilayer environments. Such rapid convergence also leads to well-converged interaction patterns between the TM helix and the membrane, which clearly show the interactions of interfacial membrane-anchoring residues with the lipids. For the N-terminal periplasmic helix of Pf1, we identify a stable, albeit dynamic, helix orientation parallel to the membrane surface that satisfies the amphiphatic nature of the helix in an explicit lipid bilayer. Such detailed information cannot be obtained solely from NMR observables. Therefore, the present simulations illustrate the usefulness of NMR-restrained MD refinement of membrane protein structure in explicit membranes.     

Cysteine-scanning mutagenesis and disulfide mapping studies of the conserved domain of the twin-arginine translocase TatB component

The cytoplasmic membrane protein TatB is an essential component of the Escherichia coli twin-arginine (Tat) protein translocation pathway. Together with the TatC component it forms a complex that functions as a membrane receptor for substrate proteins. Structural predictions suggest that TatB is anchored to the membrane via an N-terminal transmembrane alpha-helix that precedes an amphipathic alpha-helical section of the protein. From truncation analysis it is known that both these regions of the protein are essential for function. Here we construct 31 unique cysteine substitutions in the first 42 residues of TatB. Each of the substitutions results in a TatB protein that is competent to support Tat-dependent protein translocation. Oxidant-induced disulfide cross-linking shows that both the N-terminal and amphipathic helices form contacts with at least one other TatB protomer. For the transmembrane helix these contacts are localized to one face of the helix. Molecular modeling and molecular dynamics simulations provide insight into the possible structural basis of the transmembrane helix interactions. Using variants with double cysteine substitutions in the transmembrane helix, we were able to detect cross-links between up to five TatB molecules. Protein purification showed that species containing at least four cross-linked TatB molecules are found in correctly assembled TatBC complexes. Our results suggest that the transmembrane helices of TatB protomers are in the center rather than the periphery of the TatBC complex.     

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

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

Formation of an unfolding intermediate state of soluble chloride intracellular channel protein CLIC1 at acidic pH

CLIC proteins function as anion channels when their structures convert from a soluble form to an integral membrane form. While very little is known about the mechanism of the conversion process, channel formation and activity are highly pH-dependent. In this study, the structural properties and conformational stability of CLIC1 were determined as a function of pH in the absence of membranes to improve our understanding of how its conformation changes when the protein encounters the acidic environment at the surface of a membrane. Although the global conformation and size of CLIC1 are not significantly altered by pH in the range of 5.5-8.2, equilibrium unfolding studies reveal that the protein molecule becomes destabilized at low pH, resulting in the formation of a highly populated intermediate with a solvent-exposed hydrophobic surface. Unlike the intermediates formed by many soluble pore-forming proteins for their insertion into membranes, the CLIC1 intermediate is not a molten globule. Acid-induced destabilization and partial unfolding of CLIC1 involve helix alpha1 which is the major structural element of the transmembrane region. We propose that the acidic environment encountered by CLICs at the surface of membranes primes the transmembrane region in the N-domain, thereby lowering the energy barrier for the conversion of soluble CLICs to their membrane-inserted forms.     

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

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

A critical evaluation of the hydropathy profile of membrane proteins

New membrane-preference scales are introduced for categories of membrane proteins with different functions. A statistical analysis is carried out with several scales to verify the relative accuracy in the prediction of the transmembrane segments of polytopic membrane proteins. The correlation between some of the scales most used and those calculated here provides criteria for selecting the most appropriate methods for a given type of protein. The parameters used in the evaluation of the hydropathy profiles have been carefully ascertained in order to develop a reliable methodology for hydropathy analysis. Finally, an integrated hydropathy analysis using different methods has been applied to several sequences of related proteins. The above analysis indicates that (a) microsomal cytochrome P450 contains only one hydrophobic region at the N-terminus that is consistently predicted to transverse the membrane: (b) only four of the six or seven putative transmembrane helices of cytochrome oxidase subunit III are predicted and correspond to helices I, III, V and VI of the previous nomenclature; (c) the product of the mitochondrial ATPase-6 gene (or the chloroplast ATPase-IV gene) of F0-F1-ATPase shows that helix IV is not consistently predicted to traverse the membrane, suggesting a four-helix model for this family of proteins.     

NMR-Based Simulation Studies of Pf1 Coat Protein in Explicit Membranes

As time- and ensemble-averaged measures, NMR observables contain information about both protein structure and dynamics. This work represents a computational study to extract such information for membrane proteins from orientation-dependent NMR observables: solid-state NMR chemical shift anisotropy and dipolar coupling, and solution NMR residual dipolar coupling. We have performed NMR-restrained molecular dynamics simulations to refine the structure of the membrane-bound form of Pf1 coat protein in explicit lipid bilayers using the recently measured chemical shift anisotropy, dipolar coupling, and residual dipolar coupling data. From the simulations, we have characterized detailed protein-lipid interactions and explored the dynamics. All simulations are stable and the NMR restraints are well satisfied. The C-terminal transmembrane (TM) domain of Pf1 finds its optimal position in the membrane quickly (within 6 ns), illustrating efficient solvation of TM domains in explicit bilayer environments. Such rapid convergence also leads to well-converged interaction patterns between the TM helix and the membrane, which clearly show the interactions of interfacial membrane-anchoring residues with the lipids. For the N-terminal periplasmic helix of Pf1, we identify a stable, albeit dynamic, helix orientation parallel to the membrane surface that satisfies the amphiphatic nature of the helix in an explicit lipid bilayer. Such detailed information cannot be obtained solely from NMR observables. Therefore, the present simulations illustrate the usefulness of NMR-restrained MD refinement of membrane protein structure in explicit membranes.