Lucy Forrest | National Institutes of Health (original) (raw)
Papers by Lucy Forrest
Biophysical Journal, 2014
Biochemistry, 2014
Uptake of neurotransmitters by sodium-coupled monoamine transporters of the NSS family is require... more Uptake of neurotransmitters by sodium-coupled monoamine transporters of the NSS family is required for termination of synaptic transmission. Transport is tightly regulated by protein−protein interactions involving the small cytoplasmic segments at the amino-and carboxy-terminal ends of the transporter. Although structures of homologues provide information about the transmembrane regions of these transporters, the structural arrangement of the terminal domains remains largely unknown. Here, we combined molecular modeling, biochemical, and biophysical approaches in an iterative manner to investigate the structure of the 82-residue N-terminal and 30-residue C-terminal domains of human serotonin transporter (SERT). Several secondary structures were predicted in these domains, and structural models were built using the Rosetta fragment-based methodology. One-dimensional 1 H nuclear magnetic resonance and circular dichroism spectroscopy supported the presence of helical elements in the isolated SERT N-terminal domain. Moreover, introducing helix-breaking residues within those elements altered the fluorescence resonance energy transfer signal between terminal cyan fluorescent protein and yellow fluorescent protein tags attached to full-length SERT, consistent with the notion that the fold of the terminal domains is relatively well-defined. Full-length models of SERT that are consistent with these and published experimental data were generated. The resultant models predict confined loci for the terminal domains and predict that they move apart during the transport-related conformational cycle, as predicted by structures of homologues and by the "rocking bundle" hypothesis, which is consistent with spectroscopic measurements. The models also suggest the nature of binding to regulatory interaction partners. This study provides a structural context for functional and regulatory mechanisms involving SERT terminal domains.
Biophysical Journal, 2015
Transporters of the SLC34 family (NaPi-IIa,b,c) catalyze uptake of inorganic phosphate (Pi) in re... more Transporters of the SLC34 family (NaPi-IIa,b,c) catalyze uptake of inorganic phosphate (Pi) in renal and intestinal epithelia. The transport cycle requires three Na(+) ions and one divalent Pi to bind before a conformational change enables translocation, intracellular release of the substrates, and reorientation of the empty carrier. The electrogenic interaction of the first Na(+) ion with NaPi-IIa/b at a postulated Na1 site is accompanied by charge displacement, and Na1 occupancy subsequently facilitates binding of a second Na(+) ion at Na2. The voltage dependence of cotransport and presteady-state charge displacements (in the absence of a complete transport cycle) are directly related to the molecular architecture of the Na1 site. The fact that Li(+) ions substitute for Na(+) at Na1, but not at the other sites (Na2 and Na3), provides an additional tool for investigating Na1 site-specific events. We recently proposed a three-dimensional model of human SLC34a1 (NaPi-IIa) including the binding sites Na2, Na3, and Pi based on the crystal structure of the dicarboxylate transporter VcINDY. Here, we propose nine residues in transmembrane helices (TM2, TM3, and TM5) that potentially contribute to Na1. To verify their roles experimentally, we made single alanine substitutions in the human NaPi-IIa isoform and investigated the kinetic properties of the mutants by voltage clamp and (32)P uptake. Substitutions at five positions in TM2 and one in TM5 resulted in relatively small changes in the substrate apparent affinities, yet at several of these positions, we observed significant hyperpolarizing shifts in the voltage dependence. Importantly, the ability of Li(+) ions to substitute for Na(+) ions was increased compared with the wild-type. Based on these findings, we adjusted the regions containing Na1 and Na3, resulting in a refined NaPi-IIa model in which five positions (T200, Q206, D209, N227, and S447) contribute directly to cation coordination at Na1.
Biophysical Journal, 2014
Phosphate plays essential biological roles and its plasma level in humans requires tight control ... more Phosphate plays essential biological roles and its plasma level in humans requires tight control to avoid bone loss (insufficiency) or vascular calcification (excess). Intestinal absorption and renal reabsorption of phosphate are mediated by members of the SLC34 family of sodium-coupled transporters (NaPi-IIa,b,c) whose membrane expression is regulated by various hormones, circulating proteins, and phosphate itself. Consequently, NaPi-II proteins are also potentially important pharmaceutical targets for controlling phosphate levels. Their crucial role in Pi homeostasis is underscored by pathologies resulting from naturally occurring SLC34 mutations and SLC34 knockout animals. SLC34 isoforms have been extensively studied with respect to transport mechanism and structure-function relationships; however, the three-dimensional structure is unknown. All SLC34 transporters share a duplicated motif comprising a glutamine followed by a stretch of threonine or serine residues, suggesting the presence of structural repeats as found in other transporter families. Nevertheless, standard bioinformatic approaches fail to clearly identify a suitable template for molecular modeling. Here, we used hydrophobicity profiles and hidden Markov models to define a structural repeat common to all SLC34 isoforms. Similar approaches identify a relationship with the core regions in a crystal structure of Vibrio cholerae Na(+)-dicarboxylate transporter VcINDY, from which we generated a homology model of human NaPi-IIa. The aforementioned SLC34 motifs in each repeat localize to the center of the model, and were predicted to form Na(+) and Pi coordination sites. Functional relevance of key amino acids was confirmed by biochemical and electrophysiological analysis of expressed, mutated transporters. Moreover, the validity of the predicted architecture is corroborated by extensive published structure-function studies. The model provides key information for elucidating the transport mechanism and predicts candidate substrate binding sites.
Biophysical Journal, 2004
This article presents the results of a detailed analysis of helix-helix interactions in membrane ... more This article presents the results of a detailed analysis of helix-helix interactions in membrane and soluble proteins. A data set of interacting pairs of helices in membrane proteins of known structure was constructed and a structure alignment algorithm was used to identify pairs of helices in soluble proteins that superimpose well with pairs of helices in the membrane-protein data set. Most helix pairs in membrane proteins are found to have a significant number of structural homologs in soluble proteins, although in some cases, primarily involving irregular helices, no close homologs exist. An analysis of geometric relationships between interacting helices in the two sets of proteins identifies some differences in the distributions of helix length, interfacial area, packing angle, and distance between the polypeptide backbones. However, a subset of soluble-protein helix pairs that are close structural homologs to membrane-protein helix pairs exhibits distributions that mirror those observed in membrane proteins. The larger average interface size and smaller distance of closest approach seen for helices in membrane proteins appears due in part to a relative enrichment of alanines and glycines, particularly as components of the AxxxA and GxxxG motifs. It is argued that membrane helices are not on average more tightly packed than helices in soluble proteins; they are simply able to approach each other more closely. This enables them to interact over longer distances, which may in turn facilitate their remaining in contact over much of the width of the lipid bilayer. The close structural similarity seen between some pairs of helices in membrane and soluble proteins suggests that packing patterns observed in soluble proteins may be useful in the modeling of membrane proteins. Moreover, there do not appear to be fundamental differences between the magnitude of the forces that drive helix packing in membrane and soluble proteins, suggesting that strategies to make membrane proteins more soluble by mutating surface residues are likely to encounter success, at least in some cases.
Biophysical Journal, 2008
The recent structural elucidation of about one dozen channels (in which we include transporters) ... more The recent structural elucidation of about one dozen channels (in which we include transporters) has provided further evidence that these membrane proteins typically undergo large movements during their function. However, it is still not well understood how these proteins achieve the necessary trade-off between stability and mobility. To identify specific structural properties of channels, we compared the helix-packing and hydrogen-bonding patterns of channels with those of membrane coils; the latter is a class of membrane proteins whose structures are expected to be more rigid. We describe in detail how in channels, helix pairs are usually arranged in packing motifs with large crossing angles (jtj % 40°), where the (small) side chains point away from the packing core and the backbones of the two helices are in close contact. We found that this contributes to a significant enrichment of Ca-HÁÁÁO bonds and to a packing geometry where right-handed parallel (t ¼ À40°6 10°) and antiparallel (t ¼ 1140°6 25°) arrangements are equally preferred. By sharp contrast, the interdigitation and hydrogen bonding of side chains in helix pairs of membrane coils results in narrowly distributed left-handed antiparallel arrangements with crossing angles t ¼ À160°6 10°(jtj % 20°). In addition, we show that these different helix-packing modes of the two types of membrane proteins correspond to specific hydrogen-bonding patterns. In particular, in channels, three times as many of the hydrogenbonded helix pairs are found in parallel right-handed motifs than are non-hydrogen-bonded helix pairs. Finally, we discuss how the presence of weak hydrogen bonds, water-containing cavities, and right-handed crossing angles may facilitate the required conformational flexibility between helix pairs of channels while maintaining sufficient structural stability.
Biophysical Journal, 2000
The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-... more The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-selective channels. The native form of the protein is 97 residues long, although peptides representing the transmembrane section display ion channel activity, which (like the native channel) is blocked by the antiviral drug amantadine. As a small ion channel, M2 may provide useful insights into more complex channel systems. Models of tetrameric bundles of helices containing either 18 or 22 residues have been simulated while embedded in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylcholine bilayer. Several different starting models have been used. These suggest that the simulation results, at least on a nanosecond time scale, are sensitive to the exact starting structure. Electrostatics calculations carried out on a ring of four ionizable aspartate residues at the N-terminal mouth of the channel suggest that at any one time, only one will be in a charged state. Helix bundle models were mostly stable over the duration of the simulation, and their helices remained tilted relative to the bilayer normal. The M2 helix bundles form closed channels that undergo breathing motions, alternating between a tetramer and a dimer-of-dimers structure. Under these conditions either the channel forms a pocket of trapped waters or it contains a column of waters broken predominantly at the C-terminal mouth of the pore. These waters exhibit restricted motion in the pore and are effectively "frozen" in a way similar to those seen in previous simulations of a proton channel formed by a four-helix bundle of a synthetic leucine-serine peptide (Randa et al.
Biophysical journal, Jan 15, 2006
In this study, we investigate the extent to which techniques for homology modeling that were deve... more In this study, we investigate the extent to which techniques for homology modeling that were developed for water-soluble proteins are appropriate for membrane proteins as well. To this end we present an assessment of current strategies for homology modeling of membrane proteins and introduce a benchmark data set of homologous membrane protein structures, called HOMEP. First, we use HOMEP to reveal the relationship between sequence identity and structural similarity in membrane proteins. This analysis indicates that homology modeling is at least as applicable to membrane proteins as it is to water-soluble proteins and that acceptable models (with C alpha-RMSD values to the native of 2 A or less in the transmembrane regions) may be obtained for template sequence identities of 30% or higher if an accurate alignment of the sequences is used. Second, we show that secondary-structure prediction algorithms that were developed for water-soluble proteins perform approximately as well for mem...
Journal of Molecular Biology, 2003
The light-driven photocycle of rhodopsin begins the photoreceptor cascade that underlies visual r... more The light-driven photocycle of rhodopsin begins the photoreceptor cascade that underlies visual response. In a sequence of events, the retinal covalently attached to the rhodopsin protein undergoes a conformational change that communicates local changes to a global conformational change throughout the whole protein. In turn, the large-scale protein change then activates G-proteins and signal amplification throughout the cell. The nature
The influenza M2 protein is a simple membrane protein, containing a single transmembrane helix. I... more The influenza M2 protein is a simple membrane protein, containing a single transmembrane helix. It is representative of a very large family of single-transmembrane helix proteins. The functional protein is a tetramer, with the four transmembrane helices forming a proton-permeable channel across the bilayer. Two independently derived models of the M2 channel domain are compared, in order to assess the success of applying molecular modelling approaches to simple membrane proteins.
Nature Communications, 2014
The Na(+)-coupled betaine symporter BetP shares a highly conserved fold with other sequence unrel... more The Na(+)-coupled betaine symporter BetP shares a highly conserved fold with other sequence unrelated secondary transporters, for example, with neurotransmitter symporters. Recently, we obtained atomic structures of BetP in distinct conformational states, which elucidated parts of its alternating-access mechanism. Here, we report a structure of BetP in a new outward-open state in complex with an anomalous scattering substrate, adding a fundamental piece to an unprecedented set of structural snapshots for a secondary transporter. In combination with molecular dynamics simulations these structural data highlight important features of the sequential formation of the substrate and sodium-binding sites, in which coordinating water molecules play a crucial role. We observe a strictly interdependent binding of betaine and sodium ions during the coupling process. All three sites undergo progressive reshaping and dehydration during the alternating-access cycle, with the most optimal coordination of all substrates found in the closed state.
Proteins: Structure, Function, and Bioinformatics, 2004
In the current report, we provide a quantitative analysis of the convergence of the sampling of c... more In the current report, we provide a quantitative analysis of the convergence of the sampling of conformational space accomplished in molecular dynamics simulations of membrane proteins of duration in the order of 10 nanoseconds. A set of proteins of diverse size and topology is considered, ranging from helical pores such as gramicidin and small -barrels such as OmpT, to larger and more complex structures such as rhodopsin and FepA. Principal component analysis of the C ␣ -atom trajectories was employed to assess the convergence of the conformational sampling in both the transmembrane domains and the whole proteins, while the time-dependence of the average structure was analyzed to obtain single-domain information. The membrane-embedded regions, particularly those of small or structurally simple proteins, were found to achieve reasonable convergence. By contrast, extramembranous domains lacking secondary structure are often markedly under-sampled, exhibiting a continuous structural drift. This drift results in a significant imprecision in the calculated B-factors, which detracts from any quantitative comparison to experimental data. In view of such limitations, we suggest that similar analyses may be valuable in simulation studies of membrane protein dynamics, in order to attach a level of confidence to any biologically relevant observations. Proteins 2004;57: 783-791.
Proteins: Structure, Function, and Bioinformatics, 2007
Proteins: Structure, Function, and Genetics, 2000
Multiple nanosecond duration molecular dynamics simulations on the pore-lining M2 helix of the ni... more Multiple nanosecond duration molecular dynamics simulations on the pore-lining M2 helix of the nicotinic acetylcholine receptor reveal how its structure and dynamics change as a function of environment. In water, the M2 helix partially unfolds to form a molecular hinge in the vicinity of a central Leu residue that has been implicated in the mechanism of ion channel gating. In a phospholipid bilayer, either as a single transmembrane helix, or as part of a pentameric helix bundle, the M2 helix shows less flexibility, but still exhibits a kink in the vicinity of the central Leu. The single M2 helix tilts relative to the bilayer normal by 12 degrees, in agreement with recent solid state NMR data (Opella et al., Nat Struct Biol 6:374-379, 1999). The pentameric helix bundle, a model for the pore domain of the nicotinic receptor and for channels formed by M2 peptides in a bilayer, is remarkably stable over a 2-ns MD simulation in a bilayer, provided one adjusts the pK(A)s of ionizable residues to their calculated values (when taking their environment into account) before starting the simulation. The resultant transbilayer pore shows fluctuations at either mouth which transiently close the channel. Proteins 2000;39:47-55.
Proteins: Structure, Function, and Bioinformatics, 2005
The addition of hydrogen atoms to models or experimental structures of proteins that contain only... more The addition of hydrogen atoms to models or experimental structures of proteins that contain only non-hydrogen atoms is a common step in crystallographic structure refinement, in theoretical studies of proteins, and in protein structure prediction. Accurate prediction of the hydrogen positions is essential, since they constitute around half of the atoms in proteins and hence contribute significantly to their energetics. Many computational tools exist for predicting hydrogen positions, although to date no quantitative comparison has been made of their accuracy or efficiency. Here we take advantage of the recent increase in ultra-high-resolution X-ray crystal structures (< 0.9 A resolution), as well as of a number of relatively high-resolution neutron diffraction structures (< 1.8 A resolution), to compare the quality of the predictions generated by a large set of commonly used methods. These include CHARMM, CNS, GROMACS, MCCE, MolProbity, WHAT IF, and X-PLOR. The hydrogen atoms that lack a rotational degree of freedom are mostly, but not always, accurately predicted. For hydrogens with a rotational degree of freedom, all the methods give much less accurate predictions. The predictions for the hydroxyl hydrogens are analyzed in detail, particularly those buried within the protein, and some explanation is provided for the errors observed. The results provide a means to make informed decisions regarding the choice and implementation of methodologies for placing hydrogens on structures of proteins. They also point to shortcomings in current force fields and suggest the need for improved descriptions of hydrogen bonding energetics.
Proceedings of the National Academy of Sciences, 2007
The recent determination of the crystal structure of the leucine transporter from Aquifex aeolicu... more The recent determination of the crystal structure of the leucine transporter from Aquifex aeolicus (aaLeuT) has provided significant insights into the function of neurotransmitter:sodium symporters. Transport by aaLeuT is Cl ؊ independent, whereas many neurotransmitter:sodium symporters from higher organisms depend on Cl ؊ ions. However, the only Cl ؊ ion identified in the aaLeuT structure interacts with nonconserved residues in extracellular loops, and thus the relevance of this binding site is unclear. Here, we use calculations of pK As and homology modeling to predict the location of a functionally important Cl ؊ binding site in serotonin transporter and other Cl ؊ -dependent transporters. We validate our model through the site-directed mutagenesis of residues predicted to coordinate the Cl ؊ ion and through the observation of sequence conservation patterns in other Cl ؊ -dependent transporters. The proposed site is located midway across the membrane and is formed by residues from transmembrane helices 2, 6, and 7. It is close to the Na1 sodium binding site, thus providing an explanation for the coupling of Cl ؊ and Na ؉ ions during transport. Other implications of the model are also discussed.
Proceedings of the National Academy of Sciences, 2008
Crystal structures of LeuT, a bacterial homologue of mammalian neurotransmitter transporters, sho... more Crystal structures of LeuT, a bacterial homologue of mammalian neurotransmitter transporters, show a molecule of bound substrate that is essentially exposed to the extracellular space but occluded from the cytoplasm. Thus, there must exist an alternate conformation for LeuT in which the substrate is accessible to the cytoplasm and a corresponding mechanism that switches accessibility from one side of the membrane to the other. Here, we identify the cytoplasmic accessibility pathway of the alternate conformation in a mammalian serotonin transporter (SERT) (a member of the same transporter family as LeuT). We also propose a model for the cytoplasmic-facing state that exploits the internal pseudosymmetry observed in the crystal structure. LeuT contains two structurally similar repeats (TMs1-5 and TMs 6 -10) that are inverted with respect to the plane of the membrane. The conformational differences between them result in the formation of the extracellular pathway. Our model for the cytoplasm-facing state exchanges the conformations of the two repeats and thus exposes the substrate and ion-binding sites to the cytoplasm. The conformational change that connects the two states primarily involves the tilting of a 4-helix bundle composed of transmembrane helices 1, 2, 6, and 7. Switching the tilt angle of this bundle is essentially equivalent to switching the conformation of the two repeats. Extensive mutagenesis of SERT and accessibility measurements, using cysteine reagents, are accommodated by our model. These observations may be of relevance to other transporter families, many of which contain internal inverted repeats.
Physiology, 2009
Crystal structures of the bacterial amino acid transporter LeuT have provided the basis for under... more Crystal structures of the bacterial amino acid transporter LeuT have provided the basis for understanding the conformational changes associated with substrate translocation by a multitude of transport proteins with the same fold. Biochemical and modeling studies led to a "rocking bundle" mechanism for LeuT that was validated by subsequent transporter structures. These advances suggest how coupled solute transport might be defined by the internal symmetry of proteins containing inverted structural repeats.
![Research paper thumbnail of Corrigendum to “Structural Asymmetry in a Trimeric Na+/Betaine Symporter, BetP, from Corynebacterium glutamicum” [J. Mol. Biol. 407/3 (2011) 368–381]](https://attachments.academia-assets.com/42864830/thumbnails/1.jpg)
](Journal of Molecular Biology, 2011
Journal of Molecular Biology, 2011
The Na + -coupled symporter BetP catalyzes the uptake of the compatible solute betaine in the soi... more The Na + -coupled symporter BetP catalyzes the uptake of the compatible solute betaine in the soil bacterium Corynebacterium glutamicum. BetP also senses hyperosmotic stress and regulates its own activity in response to stress level. We determined a three-dimensional (3D) map (at 8 Å in-plane resolution) of a constitutively active mutant of BetP in a C. glutamicum membrane environment by electron cryomicroscopy of two-dimensional crystals. The map shows that the constitutively active mutant, which lacks the C-terminal domain involved in osmosensing, is trimeric like wild-type BetP. Recently, we reported the X-ray crystal structure of BetP at 3.35 Å, in which all three protomers displayed a substrate-occluded state. Rigid-body fitting of this trimeric structure to the 3D map identified the periplasmic and cytoplasmic sides of the membrane. Fitting of an X-ray monomer to the individual protomer maps allowed assignment of transmembrane helices and of the substrate pathway, and revealed differences in trimer architecture from the X-ray structure in the tilt angle of each protomer with respect to the membrane. The three protomer maps showed pronounced differences around the substrate pathway, suggesting three different conformations within the same trimer. Two of those protomer maps closely match those of the atomic structures of the outward-facing and inward-facing states of the hydantoin transporter Mhp1, suggesting that the BetP protomer conformations reflect key states of the transport cycle. Thus, the asymmetry in the two-dimensional maps may reflect cooperativity of conformational changes within the BetP trimer, which potentially increases the rate of glycine betaine uptake.
Biophysical Journal, 2014
Biochemistry, 2014
Uptake of neurotransmitters by sodium-coupled monoamine transporters of the NSS family is require... more Uptake of neurotransmitters by sodium-coupled monoamine transporters of the NSS family is required for termination of synaptic transmission. Transport is tightly regulated by protein−protein interactions involving the small cytoplasmic segments at the amino-and carboxy-terminal ends of the transporter. Although structures of homologues provide information about the transmembrane regions of these transporters, the structural arrangement of the terminal domains remains largely unknown. Here, we combined molecular modeling, biochemical, and biophysical approaches in an iterative manner to investigate the structure of the 82-residue N-terminal and 30-residue C-terminal domains of human serotonin transporter (SERT). Several secondary structures were predicted in these domains, and structural models were built using the Rosetta fragment-based methodology. One-dimensional 1 H nuclear magnetic resonance and circular dichroism spectroscopy supported the presence of helical elements in the isolated SERT N-terminal domain. Moreover, introducing helix-breaking residues within those elements altered the fluorescence resonance energy transfer signal between terminal cyan fluorescent protein and yellow fluorescent protein tags attached to full-length SERT, consistent with the notion that the fold of the terminal domains is relatively well-defined. Full-length models of SERT that are consistent with these and published experimental data were generated. The resultant models predict confined loci for the terminal domains and predict that they move apart during the transport-related conformational cycle, as predicted by structures of homologues and by the "rocking bundle" hypothesis, which is consistent with spectroscopic measurements. The models also suggest the nature of binding to regulatory interaction partners. This study provides a structural context for functional and regulatory mechanisms involving SERT terminal domains.
Biophysical Journal, 2015
Transporters of the SLC34 family (NaPi-IIa,b,c) catalyze uptake of inorganic phosphate (Pi) in re... more Transporters of the SLC34 family (NaPi-IIa,b,c) catalyze uptake of inorganic phosphate (Pi) in renal and intestinal epithelia. The transport cycle requires three Na(+) ions and one divalent Pi to bind before a conformational change enables translocation, intracellular release of the substrates, and reorientation of the empty carrier. The electrogenic interaction of the first Na(+) ion with NaPi-IIa/b at a postulated Na1 site is accompanied by charge displacement, and Na1 occupancy subsequently facilitates binding of a second Na(+) ion at Na2. The voltage dependence of cotransport and presteady-state charge displacements (in the absence of a complete transport cycle) are directly related to the molecular architecture of the Na1 site. The fact that Li(+) ions substitute for Na(+) at Na1, but not at the other sites (Na2 and Na3), provides an additional tool for investigating Na1 site-specific events. We recently proposed a three-dimensional model of human SLC34a1 (NaPi-IIa) including the binding sites Na2, Na3, and Pi based on the crystal structure of the dicarboxylate transporter VcINDY. Here, we propose nine residues in transmembrane helices (TM2, TM3, and TM5) that potentially contribute to Na1. To verify their roles experimentally, we made single alanine substitutions in the human NaPi-IIa isoform and investigated the kinetic properties of the mutants by voltage clamp and (32)P uptake. Substitutions at five positions in TM2 and one in TM5 resulted in relatively small changes in the substrate apparent affinities, yet at several of these positions, we observed significant hyperpolarizing shifts in the voltage dependence. Importantly, the ability of Li(+) ions to substitute for Na(+) ions was increased compared with the wild-type. Based on these findings, we adjusted the regions containing Na1 and Na3, resulting in a refined NaPi-IIa model in which five positions (T200, Q206, D209, N227, and S447) contribute directly to cation coordination at Na1.
Biophysical Journal, 2014
Phosphate plays essential biological roles and its plasma level in humans requires tight control ... more Phosphate plays essential biological roles and its plasma level in humans requires tight control to avoid bone loss (insufficiency) or vascular calcification (excess). Intestinal absorption and renal reabsorption of phosphate are mediated by members of the SLC34 family of sodium-coupled transporters (NaPi-IIa,b,c) whose membrane expression is regulated by various hormones, circulating proteins, and phosphate itself. Consequently, NaPi-II proteins are also potentially important pharmaceutical targets for controlling phosphate levels. Their crucial role in Pi homeostasis is underscored by pathologies resulting from naturally occurring SLC34 mutations and SLC34 knockout animals. SLC34 isoforms have been extensively studied with respect to transport mechanism and structure-function relationships; however, the three-dimensional structure is unknown. All SLC34 transporters share a duplicated motif comprising a glutamine followed by a stretch of threonine or serine residues, suggesting the presence of structural repeats as found in other transporter families. Nevertheless, standard bioinformatic approaches fail to clearly identify a suitable template for molecular modeling. Here, we used hydrophobicity profiles and hidden Markov models to define a structural repeat common to all SLC34 isoforms. Similar approaches identify a relationship with the core regions in a crystal structure of Vibrio cholerae Na(+)-dicarboxylate transporter VcINDY, from which we generated a homology model of human NaPi-IIa. The aforementioned SLC34 motifs in each repeat localize to the center of the model, and were predicted to form Na(+) and Pi coordination sites. Functional relevance of key amino acids was confirmed by biochemical and electrophysiological analysis of expressed, mutated transporters. Moreover, the validity of the predicted architecture is corroborated by extensive published structure-function studies. The model provides key information for elucidating the transport mechanism and predicts candidate substrate binding sites.
Biophysical Journal, 2004
This article presents the results of a detailed analysis of helix-helix interactions in membrane ... more This article presents the results of a detailed analysis of helix-helix interactions in membrane and soluble proteins. A data set of interacting pairs of helices in membrane proteins of known structure was constructed and a structure alignment algorithm was used to identify pairs of helices in soluble proteins that superimpose well with pairs of helices in the membrane-protein data set. Most helix pairs in membrane proteins are found to have a significant number of structural homologs in soluble proteins, although in some cases, primarily involving irregular helices, no close homologs exist. An analysis of geometric relationships between interacting helices in the two sets of proteins identifies some differences in the distributions of helix length, interfacial area, packing angle, and distance between the polypeptide backbones. However, a subset of soluble-protein helix pairs that are close structural homologs to membrane-protein helix pairs exhibits distributions that mirror those observed in membrane proteins. The larger average interface size and smaller distance of closest approach seen for helices in membrane proteins appears due in part to a relative enrichment of alanines and glycines, particularly as components of the AxxxA and GxxxG motifs. It is argued that membrane helices are not on average more tightly packed than helices in soluble proteins; they are simply able to approach each other more closely. This enables them to interact over longer distances, which may in turn facilitate their remaining in contact over much of the width of the lipid bilayer. The close structural similarity seen between some pairs of helices in membrane and soluble proteins suggests that packing patterns observed in soluble proteins may be useful in the modeling of membrane proteins. Moreover, there do not appear to be fundamental differences between the magnitude of the forces that drive helix packing in membrane and soluble proteins, suggesting that strategies to make membrane proteins more soluble by mutating surface residues are likely to encounter success, at least in some cases.
Biophysical Journal, 2008
The recent structural elucidation of about one dozen channels (in which we include transporters) ... more The recent structural elucidation of about one dozen channels (in which we include transporters) has provided further evidence that these membrane proteins typically undergo large movements during their function. However, it is still not well understood how these proteins achieve the necessary trade-off between stability and mobility. To identify specific structural properties of channels, we compared the helix-packing and hydrogen-bonding patterns of channels with those of membrane coils; the latter is a class of membrane proteins whose structures are expected to be more rigid. We describe in detail how in channels, helix pairs are usually arranged in packing motifs with large crossing angles (jtj % 40°), where the (small) side chains point away from the packing core and the backbones of the two helices are in close contact. We found that this contributes to a significant enrichment of Ca-HÁÁÁO bonds and to a packing geometry where right-handed parallel (t ¼ À40°6 10°) and antiparallel (t ¼ 1140°6 25°) arrangements are equally preferred. By sharp contrast, the interdigitation and hydrogen bonding of side chains in helix pairs of membrane coils results in narrowly distributed left-handed antiparallel arrangements with crossing angles t ¼ À160°6 10°(jtj % 20°). In addition, we show that these different helix-packing modes of the two types of membrane proteins correspond to specific hydrogen-bonding patterns. In particular, in channels, three times as many of the hydrogenbonded helix pairs are found in parallel right-handed motifs than are non-hydrogen-bonded helix pairs. Finally, we discuss how the presence of weak hydrogen bonds, water-containing cavities, and right-handed crossing angles may facilitate the required conformational flexibility between helix pairs of channels while maintaining sufficient structural stability.
Biophysical Journal, 2000
The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-... more The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-selective channels. The native form of the protein is 97 residues long, although peptides representing the transmembrane section display ion channel activity, which (like the native channel) is blocked by the antiviral drug amantadine. As a small ion channel, M2 may provide useful insights into more complex channel systems. Models of tetrameric bundles of helices containing either 18 or 22 residues have been simulated while embedded in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylcholine bilayer. Several different starting models have been used. These suggest that the simulation results, at least on a nanosecond time scale, are sensitive to the exact starting structure. Electrostatics calculations carried out on a ring of four ionizable aspartate residues at the N-terminal mouth of the channel suggest that at any one time, only one will be in a charged state. Helix bundle models were mostly stable over the duration of the simulation, and their helices remained tilted relative to the bilayer normal. The M2 helix bundles form closed channels that undergo breathing motions, alternating between a tetramer and a dimer-of-dimers structure. Under these conditions either the channel forms a pocket of trapped waters or it contains a column of waters broken predominantly at the C-terminal mouth of the pore. These waters exhibit restricted motion in the pore and are effectively "frozen" in a way similar to those seen in previous simulations of a proton channel formed by a four-helix bundle of a synthetic leucine-serine peptide (Randa et al.
Biophysical journal, Jan 15, 2006
In this study, we investigate the extent to which techniques for homology modeling that were deve... more In this study, we investigate the extent to which techniques for homology modeling that were developed for water-soluble proteins are appropriate for membrane proteins as well. To this end we present an assessment of current strategies for homology modeling of membrane proteins and introduce a benchmark data set of homologous membrane protein structures, called HOMEP. First, we use HOMEP to reveal the relationship between sequence identity and structural similarity in membrane proteins. This analysis indicates that homology modeling is at least as applicable to membrane proteins as it is to water-soluble proteins and that acceptable models (with C alpha-RMSD values to the native of 2 A or less in the transmembrane regions) may be obtained for template sequence identities of 30% or higher if an accurate alignment of the sequences is used. Second, we show that secondary-structure prediction algorithms that were developed for water-soluble proteins perform approximately as well for mem...
Journal of Molecular Biology, 2003
The light-driven photocycle of rhodopsin begins the photoreceptor cascade that underlies visual r... more The light-driven photocycle of rhodopsin begins the photoreceptor cascade that underlies visual response. In a sequence of events, the retinal covalently attached to the rhodopsin protein undergoes a conformational change that communicates local changes to a global conformational change throughout the whole protein. In turn, the large-scale protein change then activates G-proteins and signal amplification throughout the cell. The nature
The influenza M2 protein is a simple membrane protein, containing a single transmembrane helix. I... more The influenza M2 protein is a simple membrane protein, containing a single transmembrane helix. It is representative of a very large family of single-transmembrane helix proteins. The functional protein is a tetramer, with the four transmembrane helices forming a proton-permeable channel across the bilayer. Two independently derived models of the M2 channel domain are compared, in order to assess the success of applying molecular modelling approaches to simple membrane proteins.
Nature Communications, 2014
The Na(+)-coupled betaine symporter BetP shares a highly conserved fold with other sequence unrel... more The Na(+)-coupled betaine symporter BetP shares a highly conserved fold with other sequence unrelated secondary transporters, for example, with neurotransmitter symporters. Recently, we obtained atomic structures of BetP in distinct conformational states, which elucidated parts of its alternating-access mechanism. Here, we report a structure of BetP in a new outward-open state in complex with an anomalous scattering substrate, adding a fundamental piece to an unprecedented set of structural snapshots for a secondary transporter. In combination with molecular dynamics simulations these structural data highlight important features of the sequential formation of the substrate and sodium-binding sites, in which coordinating water molecules play a crucial role. We observe a strictly interdependent binding of betaine and sodium ions during the coupling process. All three sites undergo progressive reshaping and dehydration during the alternating-access cycle, with the most optimal coordination of all substrates found in the closed state.
Proteins: Structure, Function, and Bioinformatics, 2004
In the current report, we provide a quantitative analysis of the convergence of the sampling of c... more In the current report, we provide a quantitative analysis of the convergence of the sampling of conformational space accomplished in molecular dynamics simulations of membrane proteins of duration in the order of 10 nanoseconds. A set of proteins of diverse size and topology is considered, ranging from helical pores such as gramicidin and small -barrels such as OmpT, to larger and more complex structures such as rhodopsin and FepA. Principal component analysis of the C ␣ -atom trajectories was employed to assess the convergence of the conformational sampling in both the transmembrane domains and the whole proteins, while the time-dependence of the average structure was analyzed to obtain single-domain information. The membrane-embedded regions, particularly those of small or structurally simple proteins, were found to achieve reasonable convergence. By contrast, extramembranous domains lacking secondary structure are often markedly under-sampled, exhibiting a continuous structural drift. This drift results in a significant imprecision in the calculated B-factors, which detracts from any quantitative comparison to experimental data. In view of such limitations, we suggest that similar analyses may be valuable in simulation studies of membrane protein dynamics, in order to attach a level of confidence to any biologically relevant observations. Proteins 2004;57: 783-791.
Proteins: Structure, Function, and Bioinformatics, 2007
Proteins: Structure, Function, and Genetics, 2000
Multiple nanosecond duration molecular dynamics simulations on the pore-lining M2 helix of the ni... more Multiple nanosecond duration molecular dynamics simulations on the pore-lining M2 helix of the nicotinic acetylcholine receptor reveal how its structure and dynamics change as a function of environment. In water, the M2 helix partially unfolds to form a molecular hinge in the vicinity of a central Leu residue that has been implicated in the mechanism of ion channel gating. In a phospholipid bilayer, either as a single transmembrane helix, or as part of a pentameric helix bundle, the M2 helix shows less flexibility, but still exhibits a kink in the vicinity of the central Leu. The single M2 helix tilts relative to the bilayer normal by 12 degrees, in agreement with recent solid state NMR data (Opella et al., Nat Struct Biol 6:374-379, 1999). The pentameric helix bundle, a model for the pore domain of the nicotinic receptor and for channels formed by M2 peptides in a bilayer, is remarkably stable over a 2-ns MD simulation in a bilayer, provided one adjusts the pK(A)s of ionizable residues to their calculated values (when taking their environment into account) before starting the simulation. The resultant transbilayer pore shows fluctuations at either mouth which transiently close the channel. Proteins 2000;39:47-55.
Proteins: Structure, Function, and Bioinformatics, 2005
The addition of hydrogen atoms to models or experimental structures of proteins that contain only... more The addition of hydrogen atoms to models or experimental structures of proteins that contain only non-hydrogen atoms is a common step in crystallographic structure refinement, in theoretical studies of proteins, and in protein structure prediction. Accurate prediction of the hydrogen positions is essential, since they constitute around half of the atoms in proteins and hence contribute significantly to their energetics. Many computational tools exist for predicting hydrogen positions, although to date no quantitative comparison has been made of their accuracy or efficiency. Here we take advantage of the recent increase in ultra-high-resolution X-ray crystal structures (< 0.9 A resolution), as well as of a number of relatively high-resolution neutron diffraction structures (< 1.8 A resolution), to compare the quality of the predictions generated by a large set of commonly used methods. These include CHARMM, CNS, GROMACS, MCCE, MolProbity, WHAT IF, and X-PLOR. The hydrogen atoms that lack a rotational degree of freedom are mostly, but not always, accurately predicted. For hydrogens with a rotational degree of freedom, all the methods give much less accurate predictions. The predictions for the hydroxyl hydrogens are analyzed in detail, particularly those buried within the protein, and some explanation is provided for the errors observed. The results provide a means to make informed decisions regarding the choice and implementation of methodologies for placing hydrogens on structures of proteins. They also point to shortcomings in current force fields and suggest the need for improved descriptions of hydrogen bonding energetics.
Proceedings of the National Academy of Sciences, 2007
The recent determination of the crystal structure of the leucine transporter from Aquifex aeolicu... more The recent determination of the crystal structure of the leucine transporter from Aquifex aeolicus (aaLeuT) has provided significant insights into the function of neurotransmitter:sodium symporters. Transport by aaLeuT is Cl ؊ independent, whereas many neurotransmitter:sodium symporters from higher organisms depend on Cl ؊ ions. However, the only Cl ؊ ion identified in the aaLeuT structure interacts with nonconserved residues in extracellular loops, and thus the relevance of this binding site is unclear. Here, we use calculations of pK As and homology modeling to predict the location of a functionally important Cl ؊ binding site in serotonin transporter and other Cl ؊ -dependent transporters. We validate our model through the site-directed mutagenesis of residues predicted to coordinate the Cl ؊ ion and through the observation of sequence conservation patterns in other Cl ؊ -dependent transporters. The proposed site is located midway across the membrane and is formed by residues from transmembrane helices 2, 6, and 7. It is close to the Na1 sodium binding site, thus providing an explanation for the coupling of Cl ؊ and Na ؉ ions during transport. Other implications of the model are also discussed.
Proceedings of the National Academy of Sciences, 2008
Crystal structures of LeuT, a bacterial homologue of mammalian neurotransmitter transporters, sho... more Crystal structures of LeuT, a bacterial homologue of mammalian neurotransmitter transporters, show a molecule of bound substrate that is essentially exposed to the extracellular space but occluded from the cytoplasm. Thus, there must exist an alternate conformation for LeuT in which the substrate is accessible to the cytoplasm and a corresponding mechanism that switches accessibility from one side of the membrane to the other. Here, we identify the cytoplasmic accessibility pathway of the alternate conformation in a mammalian serotonin transporter (SERT) (a member of the same transporter family as LeuT). We also propose a model for the cytoplasmic-facing state that exploits the internal pseudosymmetry observed in the crystal structure. LeuT contains two structurally similar repeats (TMs1-5 and TMs 6 -10) that are inverted with respect to the plane of the membrane. The conformational differences between them result in the formation of the extracellular pathway. Our model for the cytoplasm-facing state exchanges the conformations of the two repeats and thus exposes the substrate and ion-binding sites to the cytoplasm. The conformational change that connects the two states primarily involves the tilting of a 4-helix bundle composed of transmembrane helices 1, 2, 6, and 7. Switching the tilt angle of this bundle is essentially equivalent to switching the conformation of the two repeats. Extensive mutagenesis of SERT and accessibility measurements, using cysteine reagents, are accommodated by our model. These observations may be of relevance to other transporter families, many of which contain internal inverted repeats.
Physiology, 2009
Crystal structures of the bacterial amino acid transporter LeuT have provided the basis for under... more Crystal structures of the bacterial amino acid transporter LeuT have provided the basis for understanding the conformational changes associated with substrate translocation by a multitude of transport proteins with the same fold. Biochemical and modeling studies led to a "rocking bundle" mechanism for LeuT that was validated by subsequent transporter structures. These advances suggest how coupled solute transport might be defined by the internal symmetry of proteins containing inverted structural repeats.
Journal of Molecular Biology, 2011
Journal of Molecular Biology, 2011
The Na + -coupled symporter BetP catalyzes the uptake of the compatible solute betaine in the soi... more The Na + -coupled symporter BetP catalyzes the uptake of the compatible solute betaine in the soil bacterium Corynebacterium glutamicum. BetP also senses hyperosmotic stress and regulates its own activity in response to stress level. We determined a three-dimensional (3D) map (at 8 Å in-plane resolution) of a constitutively active mutant of BetP in a C. glutamicum membrane environment by electron cryomicroscopy of two-dimensional crystals. The map shows that the constitutively active mutant, which lacks the C-terminal domain involved in osmosensing, is trimeric like wild-type BetP. Recently, we reported the X-ray crystal structure of BetP at 3.35 Å, in which all three protomers displayed a substrate-occluded state. Rigid-body fitting of this trimeric structure to the 3D map identified the periplasmic and cytoplasmic sides of the membrane. Fitting of an X-ray monomer to the individual protomer maps allowed assignment of transmembrane helices and of the substrate pathway, and revealed differences in trimer architecture from the X-ray structure in the tilt angle of each protomer with respect to the membrane. The three protomer maps showed pronounced differences around the substrate pathway, suggesting three different conformations within the same trimer. Two of those protomer maps closely match those of the atomic structures of the outward-facing and inward-facing states of the hydantoin transporter Mhp1, suggesting that the BetP protomer conformations reflect key states of the transport cycle. Thus, the asymmetry in the two-dimensional maps may reflect cooperativity of conformational changes within the BetP trimer, which potentially increases the rate of glycine betaine uptake.