Modeling and simulation of ion-coupled and ATP-driven membrane proteins

In silico characterization of residues essential for substrate binding of human cystine transporter, xCT

xCT is a sodium-independent amino acid antiporter that imports L-cystine and exports L-glutamate in a 1:1 ratio. It is a component of heterodimeric amino acid transporter system Xc- working at the cross-roads of maintaining neurological processes and regulating antioxidant defense. The transporter has 12 transmembrane domains with intracellular N- and C-termini, and like other transporter proteins can undergo various conformational changes while switching the ligand accessibilities from intracellular to extracellular site. In the present study, we generated two homology models of human xCT in two distinct conformations: inward-facing occluded state and outward-facing open state. Our results indicated the substrate translocation channel composed of transmembrane helices TMs 1, 3, 6, 8, and 10. We docked anionic L-cystine and L-glutamate within the cavities to assess the two distinct binding scenarios for xCT as antiporter. We also assessed the interactions between the ligands and transporter and observed that ligands bind to similar residues within the channel. Using MM-PBSA/MM-GBSA approach, we computed the binding energies of these ligands to different conformational states. Cystine and glutamate bind xCT with favorable binding energies, with more favorable binding observed in inward occluded state than in outward open state. We further computed the residue-wise decomposition of these binding energies and identified the residues as essential for substrate binding/permeation. Filtering the residues that form favorable energetic contributions to the ligand binding in both the states, our studies suggest T56, A60, R135, A138, V141, Y244, A247, F250, S330, L392, and R396 as critical residues for ligand binding as well as ligand transport for any conformational state adopted by xCT during its transport cycle. .Graphical Abstract.

Interpretation of spectroscopic data using molecular simulations for the secondary active transporter BetP

Mechanistic understanding of dynamic membrane proteins such as transporters, receptors, and channels requires accurate depictions of conformational ensembles, and the manner in which they interchange as a function of environmental factors including substrates, lipids, and inhibitors. Spectroscopic techniques such as electron spin resonance (ESR) pulsed electron-electron double resonance (PELDOR), also known as double electron-electron resonance (DEER), provide a complement to atomistic structures obtained from x-ray crystallography or cryo-EM, since spectroscopic data reflect an ensemble and can be measured in more native solvents, unperturbed by a crystal lattice. However, attempts to interpret DEER data are frequently stymied by discrepancies with the structural data, which may arise due to differences in conditions, the dynamics of the protein, or the flexibility of the attached paramagnetic spin labels. Recently, molecular simulation techniques such as EBMetaD have been developed that create a conformational ensemble matching an experimental distance distribution while applying the minimal possible bias. Moreover, it has been proposed that the work required during an EBMetaD simulation to match an experimentally determined distribution could be used as a metric with which to assign conformational states to a given measurement. Here, we demonstrate the application of this concept for a sodium-coupled transport protein, BetP. Because the probe, protein, and lipid bilayer are all represented in atomic detail, the different contributions to the work, such as the extent of protein backbone movements, can be separated. This work therefore illustrates how ranking simulations based on EBMetaD can help to bridge the gap between structural and biophysical data and thereby enhance our understanding of membrane protein conformational mechanisms.

Molecular Modeling of Drug-Transporter Interactions-An International Transporter Consortium Perspective

Membrane transporters play diverse roles in the pharmacokinetics and pharmacodynamics of small-molecule drugs. Understanding the mechanisms of drug-transporter interactions at the molecular level is, therefore, essential for the design of drugs with optimal therapeutic effects. This white paper examines recent progress, applications, and challenges of molecular modeling of membrane transporters, including modeling techniques that are centered on the structures of transporter ligands, and those focusing on the structures of the transporters. The goals of this article are to illustrate current best practices and future opportunities in using molecular modeling techniques to understand and predict transporter-mediated effects on drug disposition and efficacy.Membrane transporters from the solute carrier (SLC) and ATP-binding cassette (ABC) superfamilies regulate the cellular uptake, efflux, and homeostasis of many essential nutrients and significantly impact the pharmacokinetics of drugs; further, they may provide targets for novel therapeutics as well as facilitate prodrug approaches. Because of their often broad substrate selectivity they are also implicated in many undesirable and sometimes life-threatening drug-drug interactions (DDIs).^5,6.

Mechanism of Substrate Translocation in an Alternating Access Transporter

Transporters shuttle molecules across cell membranes by alternating among distinct conformational states. Fundamental questions remain about how transporters transition between states and how such structural rearrangements regulate substrate translocation. Here, we capture the translocation process by crystallography and unguided molecular dynamics simulations, providing an atomic-level description of alternating access transport. Simulations of a SWEET-family transporter initiated from an outward-open, glucose-bound structure reported here spontaneously adopt occluded and inward-open conformations. Strikingly, these conformations match crystal structures, including our inward-open structure. Mutagenesis experiments further validate simulation predictions. Our results reveal that state transitions are driven by favorable interactions formed upon closure of extracellular and intracellular "gates" and by an unfavorable transmembrane helix configuration when both gates are closed. This mechanism leads to tight allosteric coupling between gates, preventing them from opening simultaneously. Interestingly, the substrate appears to take a "free ride" across the membrane without causing major structural rearrangements in the transporter.

Sodium recognition by the Na+/Ca2+ exchanger in the outward-facing conformation

Na(+)/Ca(2+) exchangers (NCXs) are ubiquitous membrane transporters with a key role in Ca(2+) homeostasis and signaling. NCXs mediate the bidirectional translocation of either Na(+) or Ca(2+), and thus can catalyze uphill Ca(2+) transport driven by a Na(+) gradient, or vice versa. In a major breakthrough, a prokaryotic NCX homolog (NCX_Mj) was recently isolated and its crystal structure determined at atomic resolution. The structure revealed an intriguing architecture consisting of two inverted-topology repeats, each comprising five transmembrane helices. These repeats adopt asymmetric conformations, yielding an outward-facing occluded state. The crystal structure also revealed four putative ion-binding sites, but the occupancy and specificity thereof could not be conclusively established. Here, we use molecular-dynamics simulations and free-energy calculations to identify the ion configuration that best corresponds to the crystallographic data and that is also thermodynamically optimal. In this most probable configuration, three Na(+) ions occupy the so-called Sext, SCa, and Sint sites, whereas the Smid site is occupied by one water molecule and one H(+), which protonates an adjacent aspartate side chain (D240). Experimental measurements of Na(+)/Ca(2+) and Ca(2+)/Ca(2+) exchange by wild-type and mutagenized NCX_Mj confirm that transport of both Na(+) and Ca(2+) requires protonation of D240, and that this side chain does not coordinate either ion at Smid. These results imply that the ion exchange stoichiometry of NCX_Mj is 3:1 and that translocation of Na(+) across the membrane is electrogenic, whereas transport of Ca(2+) is not. Altogether, these findings provide the basis for further experimental and computational studies of the conformational mechanism of this exchanger.

Structural fold and binding sites of the human Na⁺-phosphate cotransporter NaPi-II

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.

Molecular modeling and ligand docking for solute carrier (SLC) transporters

Solute Carrier (SLC) transporters are membrane proteins that transport solutes, such as ions, metabolites, peptides, and drugs, across biological membranes, using diverse energy coupling mechanisms. In human, there are 386 SLC transporters, many of which contribute to the absorption, distribution, metabolism, and excretion of drugs and/or can be targeted directly by therapeutics. Recent atomic structures of SLC transporters determined by X-ray crystallography and NMR spectroscopy have significantly expanded the applicability of structure-based prediction of SLC transporter ligands, by enabling both comparative modeling of additional SLC transporters and virtual screening of small molecules libraries against experimental structures as well as comparative models. In this review, we begin by describing computational tools, including sequence analysis, comparative modeling, and virtual screening, that are used to predict the structures and functions of membrane proteins such as SLC transporters. We then illustrate the applications of these tools to predicting ligand specificities of select SLC transporters, followed by experimental validation using uptake kinetic measurements and other assays. We conclude by discussing future directions in the discovery of the SLC transporter ligands.

Investigation of the sodium-binding sites in the sodium-coupled betaine transporter BetP

Sodium-coupled substrate transport plays a central role in many biological processes. However, despite knowledge of the structures of several sodium-coupled transporters, the location of the sodium-binding site(s) often remains unclear. Several of these structures have the five transmembrane-helix inverted-topology repeat, LeuT-like (FIRL) fold, whose pseudosymmetry has been proposed to facilitate the alternating-access mechanism required for transport. Here, we provide biophysical, biochemical, and computational evidence for the location of the two cation-binding sites in the sodium-coupled betaine symporter BetP. A recent X-ray structure of BetP in a sodium-bound closed state revealed that one of these sites, equivalent to the Na2 site in related transporters, is located between transmembrane helices 1 and 8 of the FIRL-fold; here, we confirm the location of this site by other means. Based on the pseudosymmetry of this fold, we hypothesized that the second site is located between the equivalent helices 6 and 3. Molecular dynamics simulations of the closed-state structure suggest this second sodium site involves two threonine sidechains and a backbone carbonyl from helix 3, a phenylalanine from helix 6, and a water molecule. Mutating the residues proposed to form the two binding sites increased the apparent K(m) and K(d) for sodium, as measured by betaine uptake, tryptophan fluorescence, and (22)Na(+) binding, and also diminished the transient currents measured in proteoliposomes using solid supported membrane-based electrophysiology. Taken together, these results provide strong evidence for the identity of the residues forming the sodium-binding sites in BetP.

Structural genomics plucks high-hanging membrane proteins

Recent years have seen the establishment of structural genomics centers that explicitly target integral membrane proteins. Here, we review the advances in targeting these extremely high-hanging fruits of structural biology in high-throughput mode. We observe that the experimental determination of high-resolution structures of integral membrane proteins is increasingly successful both in terms of getting structures and of covering important protein families, for example, from Pfam. Structural genomics has begun to contribute significantly toward this progress. An important component of this contribution is the set up of robotic pipelines that generate a wealth of experimental data for membrane proteins. We argue that prediction methods for the identification of membrane regions and for the comparison of membrane proteins largely suffice to meet the challenges of target selection for structural genomics of membrane proteins. In contrast, we need better methods to prioritize the most promising members in a family of closely related proteins and to annotate protein function from sequence and structure in absence of homology.

Atomistic models of ion and solute transport by the sodium-dependent secondary active transporters

The recent determination of high-resolution crystal structures of several transporters offers unprecedented insights into the structural mechanisms behind secondary transport. These proteins utilize the facilitated diffusion of the ions down their electrochemical gradients to transport the substrate against its concentration gradient. The structural studies revealed striking similarities in the structural organization of ion and solute binding sites and a well-conserved inverted-repeat topology between proteins from several gene families. In this paper we will overview recent atomistic simulations applied to study the mechanisms of selective binding of ion and substrate in LeuT, Glt, vSGLT and hSERT as well as its consequences for the transporter conformational dynamics. This article is part of a Special Issue entitled: Membrane protein structure and function.