Relation between sequence and structure in membrane proteins

Emulating proton-induced conformational changes in the vesicular monoamine transporter VMAT2 by mutagenesis

Neurotransporters located in synaptic vesicles are essential for communication between nerve cells in a process mediated by neurotransmitters. Vesicular monoamine transporter (VMAT), a member of the largest superfamily of transporters, mediates transport of monoamines to synaptic vesicles and storage organelles in a process that involves exchange of two H⁺ per substrate. VMAT transport is inhibited by the competitive inhibitor reserpine, a second-line agent to treat hypertension, and by the noncompetitive inhibitor tetrabenazine, presently in use for symptomatic treatment of hyperkinetic disorders. During the transport cycle, VMAT is expected to occupy at least three different conformations: cytoplasm-facing, occluded, and lumen-facing. The lumen- to cytoplasm-facing transition, facilitated by protonation of at least one of the essential membrane-embedded carboxyls, generates a binding site for reserpine. Here we have identified residues in the cytoplasmic gate and show that mutations that disrupt the interactions in this gate also shift the equilibrium toward the cytoplasm-facing conformation, emulating the effect of protonation. These experiments provide significant insight into the role of proton translocation in the conformational dynamics of a mammalian H⁺-coupled antiporter, and also identify key aspects of the mode of action and binding of two potent inhibitors of VMAT2: reserpine binds the cytoplasm-facing conformation, and tetrabenazine binds the lumen-facing conformation.

Warfarin resistance associated with genetic polymorphism of VKORC1: linking clinical response to molecular mechanism using computational modeling

The variable response to warfarin treatment often has a genetic basis. A protein homology model of human vitamin K epoxide reductase, subunit 1 (VKORC1), was generated to elucidate the mechanism of warfarin resistance observed in a patient with the Val66Met mutation. The VKORC1 homology model comprises four transmembrane (TM) helical domains and a half helical lid domain. Cys132 and Cys135, located in the N-terminal end of TM-4, are linked through a disulfide bond. Two distinct binding sites for warfarin were identified. Site-1, which binds vitamin K epoxide (KO) in a catalytically favorable orientation, shows higher affinity for S-warfarin compared with R-warfarin. Site-2, positioned in the domain occupied by the hydrophobic tail of KO, binds both warfarin enantiomers with similar affinity. Displacement of Arg37 occurs in the Val66Met mutant, blocking access of warfarin (but not KO) to Site-1, consistent with clinical observation of warfarin resistance.

On the role of a specific insert in acetate permeases (ActP) for tellurite uptake in bacteria: Functional and structural studies

The oxyanion tellurite (TeO₃^2-) is extremely toxic to bacterial cells. In Rhodobacter capsulatus, tellurite enters the cytosol by means of the high uptake-rate acetate permease RcActP2, encoded by one of the three actP genes present in this species (actP1, actP2 and actP3). Conversely, in Escherichia coli a low rate influx of the oxyanion is measured, which depends mainly on the phosphate transporter EcPitA, even though E. coli contains its own EcActP acetate permease. Here we report that when the actP2 gene from R. capsulatus is expressed in wild-type E. coli HB101 and in E. coli JW3460 ΔpitA mutant, the cellular intake of tellurite increases up to four times, suggesting intrinsic structural differences between EcActP and RcActP2. Indeed, a sequence analysis indicated the presence in RcActP2 of an insert of 15-16 residues, located between trans-membrane (TM) helices 6 and 7, which is absent in both EcActP and RcActP1. Based on this observation, the molecular models of homodimeric RcActP1 and RcActP2 were calculated and analyzed. In the RcActP2 model, the insert induces a perturbation in the conformation of the loop between TM helices 6 and 7, located at the RcActP2 dimerization interface. This perturbation opens a cavity on the periplasmic side that is closed, instead, in the RcActP1 model. This cavity also features an increase of the positive electric potential on the protein surface, an effect ascribed to specific residues Lys261, Lys281 and Arg560. We propose that this positively charged patch in RcActP2 is involved in recognition and translocation of the TeO₃^2- anion, attributing to RcActP2 a greater ability as compared to RcActP1 to transport this inorganic poison inside the cells.

The bacterial dicarboxylate transporter VcINDY uses a two-domain elevator-type mechanism

Secondary transporters use alternating access mechanisms to couple uphill substrate movement to downhill ion flux. Most known transporters utilize a “rocking bundle” motion, where the protein moves around an immobile substrate binding site. However, the glutamate transporter homolog, Glt_Ph, translocates its substrate binding site vertically across the membrane, an “elevator” mechanism. Here, we used the “repeat swap” approach to computationally predict the outward-facing state of the Na⁺/succinate transporter VcINDY, from Vibrio cholerae. Our model predicts a substantial “elevator”-like movement of vcINDY’s substrate binding site, with a vertical translation of ~15 Å and a rotation of ~43°; multiple disulfide crosslinks which completely inhibit transport provide experimental confirmation and demonstrate that such movement is essential. In contrast, crosslinks across the VcINDY dimer interface preserve transport, revealing an absence of large scale coupling between protomers.

TMFoldWeb: a web server for predicting transmembrane protein fold class

Background

Here we present TMFoldWeb, the web server implementation of TMFoldRec, a transmembrane protein fold recognition algorithm. TMFoldRec uses statistical potentials and utilizes topology filtering and a gapless threading algorithm. It ranks template structures and selects the most likely candidates and estimates the reliability of the obtained lowest energy model. The statistical potential was developed in a maximum likelihood framework on a representative set of the PDBTM database. According to the benchmark test the performance of TMFoldRec is about 77 % in correctly predicting fold class for a given transmembrane protein sequence.

Results

An intuitive web interface has been developed for the recently published TMFoldRec algorithm. The query sequence goes through a pipeline of topology prediction and a systematic sequence to structure alignment (threading). Resulting templates are ordered by energy and reliability values and are colored according to their significance level. Besides the graphical interface, a programmatic access is available as well, via a direct interface for developers or for submitting genome-wide data sets.

Conclusions

The TMFoldWeb web server is unique and currently the only web server that is able to predict the fold class of transmembrane proteins while assigning reliability scores for the prediction. This method is prepared for genome-wide analysis with its easy-to-use interface, informative result page and programmatic access.

Considering the info-communication evolution in the last few years, the developed web server, as well as the molecule viewer, is responsive and fully compatible with the prevalent tablets and mobile devices.

Reviewers

This article was reviewed by Dr. Michael Gromiha, Dr. Sandor Pongor and Dr. Frank Eisenhaber with Wing-Cheong Wong.

TMalphaDB and TMbetaDB: web servers to study the structural role of sequence motifs in α-helix and β-barrel domains of membrane proteins

Background

Membrane proteins represent over 25 % of human protein genes and account for more than 60 % of drug targets due to their accessibility from the extracellular environment. The increasing number of available crystal structures of these proteins in the Protein Data Bank permits an initial estimation of their structural properties.

Description

We have developed two web servers—TMalphaDB for α-helix bundles and TMbetaDB for β-barrels—to analyse the growing repertoire of available crystal structures of membrane proteins. TMalphaDB and TMbetaDB permit to search for these specific sequence motifs in a non-redundant structure database of transmembrane segments and quantify structural parameters such as ϕ and ψ backbone dihedral angles, χ₁ side chain torsion angle, unit bend and unit twist.

Conclusions

The structural information offered by TMalphaDB and TMbetaDB permits to quantify structural distortions induced by specific sequence motifs, and to elucidate their role in the 3D structure. This specific structural information has direct implications in homology modeling of the growing sequences of membrane proteins lacking experimental structure. TMalphaDB and TMbetaDB are freely available at http://lmc.uab.cat/TMalphaDB and http://lmc.uab.cat/TMbetaDB.

GPCRtm: An amino acid substitution matrix for the transmembrane region of class A G Protein-Coupled Receptors

Background

Protein sequence alignments and database search methods use standard scoring matrices calculated from amino acid substitution frequencies in general sets of proteins. These general-purpose matrices are not optimal to align accurately sequences with marked compositional biases, such as hydrophobic transmembrane regions found in membrane proteins. In this work, an amino acid substitution matrix (GPCRtm) is calculated for the membrane spanning segments of the G protein-coupled receptor (GPCR) rhodopsin family; one of the largest transmembrane protein family in humans with great importance in health and disease.

Results

The GPCRtm matrix reveals the amino acid compositional bias distinctive of the GPCR rhodopsin family and differs from other standard substitution matrices. These membrane receptors, as expected, are characterized by a high content of hydrophobic residues with regard to globular proteins. On the other hand, the presence of polar and charged residues is higher than in average membrane proteins, displaying high frequencies of replacement within themselves.

Conclusions

Analysis of amino acid frequencies and values obtained from the GPCRtm matrix reveals patterns of residue replacements different from other standard substitution matrices. GPCRs prioritize the reactivity properties of the amino acids over their bulkiness in the transmembrane regions. A distinctive role is that charged and polar residues seem to evolve at different rates than other amino acids. This observation is related to the role of the transmembrane bundle in the binding of ligands, that in many cases involve electrostatic and hydrogen bond interactions. This new matrix can be useful in database search and for the construction of more accurate sequence alignments of GPCRs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0639-4) contains supplementary material, which is available to authorized users.

Structural insights into thyroid hormone transport mechanisms of the L-type amino acid transporter 2

Thyroid hormones (THs) are transported across cell membranes by different transmembrane transporter proteins. In previous studies, we showed marked 3,3'-diiodothyronine (3,3'-T2) but moderate T3 uptake by the L-type amino acid transporter 2 (Lat2). We have now studied the structure-function relationships of this transporter and TH-like molecules. Our Lat2 homology model is based on 2 crystal structures of the homologous 12-transmembrane helix transporters arginine/agmatine antiporter and amino acid/polyamine/organocation transporter. Model-driven mutagenesis of residues lining an extracellular recognition site and a TH-traversing channel identified 9 sensitive residues. Using Xenopus laevis oocytes as expression system, we found that side chain shortening (N51S, N133S, N248S, and Y130A) expanded the channel and increased 3,3'-T2 transport. Side chain enlargements (T140F, Y130R, and I137M) decreased 3,3'-T2 uptake, indicating channel obstructions. The opposite results with mutations maintaining (F242W) or impairing (F242V) uptake suggest that F242 may have a gating function. Competitive inhibition studies of 14 TH-like compounds revealed that recognition by Lat2 requires amino and carboxylic acid groups. The size of the adjacent hydrophobic group is restricted. Bulky substituents in positions 3 and 5 of the tyrosine ring are allowed. The phenolic ring may be enlarged, provided that the whole molecule is flexible enough to fit into the distinctly shaped TH-traversing channel of Lat2. Taken together, the next Lat2 features were identified 1) TH recognition site; 2) TH-traversing channel in the center of Lat2; and 3) switch site that potentially facilitates intracellular substrate release. Together with identified substrate features, these data help to elucidate the molecular mechanisms and role of Lat2 in T2 transport.

De novo membrane protein structure prediction

Recent advances in identifying residue-residue contacts from large multiple sequence alignments have enabled impressive gains to be made in the field of protein structure prediction. In this chapter, we discuss these advances and provide a step-by-step guide to applying the latest tools to the de novo modelling of alpha-helical transmembrane proteins. As a practical example, we demonstrate the process of building an accurate 3D model of a G protein-coupled receptor, correctly orientated in the membrane, using only its primary protein sequence.

Ligand binding determinants for angiotensin II type 1 receptor from computer simulations

The ligand binding determinants for the angiotensin II type 1 receptor (AT1R), a G protein-coupled receptor (GPCR), have been characterized by means of computer simulations. As a first step, a pharmacophore model of various known AT1R ligands exhibiting a wide range of binding affinities was generated. Second, a structural model of AT1R was built making use of the growing set of crystal structures of GPCRs, which was further used for the docking of the AT1R ligands based on the devised pharmacophore model. Next, ligand-receptor-lipid bilayer systems were studied by means of molecular dynamics (MD) simulations. Overall, the present study has permitted, combining the pharmacophore model with binding free energy calculations obtained from the MD simulations, to propose the molecular mechanisms by which sartans interact with AT1R.

Structure-function analysis of the human ferroportin iron exporter (SLC40A1): effect of hemochromatosis type 4 disease mutations and identification of critical residues

Ferroportin (SLC40A1) is the only known iron exporter in mammals and is considered a key coordinator of the iron balance between intracellular and systemic iron homeostasis. However, the structural organization of ferroportin in the lipid bilayer remains controversial and very little is known about the mechanism underlying iron egress. In the present study, we have developed an approach based on comparative modeling, which has led to the construction of a model of the three-dimensional (3D) structure of ferroportin by homology to the crystal structure of a Major Facilitator Superfamily member (EmrD). This model predicts atomic details for the organization of ferroportin transmembrane helices and is in agreement with our current understanding of the ferroportin function and its interaction with hepcidin. Using in vitro experiments, we demonstrate that this model can be used to identify novel critical amino acids. In particular, we show that the tryptophan residue 42 (p.Trp42), which is localized within the extracellular end of the ferroportin pore, is likely involved in both the iron export function and in the mechanism of inhibition by hepcidin. Thus, our 3D model provides a new perspective for understanding the molecular basis of ferroportin functions and dysfunctions.