Tentonin 3 is a pore-forming subunit of a slow inactivation mechanosensitive channel

Mechanically activating (MA) channels transduce numerous physiological functions. Tentonin 3/TMEM150C (TTN3) confers MA currents with slow inactivation kinetics in somato- and barosensory neurons. However, questions were raised about its role as a Piezo1 regulator and its potential as a channel pore. Here, we demonstrate that purified TTN3 proteins incorporated into the lipid bilayer displayed spontaneous and pressure-sensitive channel currents. These MA currents were conserved across vertebrates and differ from Piezo1 in activation threshold and pharmacological response. Deep neural network structure prediction programs coupled with mutagenetic analysis predicted a rectangular-shaped, tetrameric structure with six transmembrane helices and a pore at the inter-subunit center. The putative pore aligned with two helices of each subunit and had constriction sites whose mutations changed the MA currents. These findings suggest that TTN3 is a pore-forming subunit of a distinct slow inactivation MA channel, potentially possessing a tetrameric structure.

Synthesis and anthelmintic activity of novel thiosemicarbazide and 1,2,4-triazole derivatives: In vitro, in vivo, and in silico study

INTRODUCTION

Intestinal parasitic infections are neglected diseases and, due to the increasing resistance of parasites to available drugs, they pose an increasing therapeutic challenge. Therefore, there is a great need for finding new compounds with antiparasitic activity.

OBJECTIVES

In this work, new thiosemicarbazide and 1,2,4-triazole derivatives were synthesized and tested for their anthelmintic activity.

METHODS

The synthesis was carried out by classical methods of organic chemistry. Anthelmintic activity tests were carried out in vitro (Rhabditis sp., Haemonchus contortus, Strongylidae sp.) in vivo (Heligmosomoides polygyrus/bakeri), and in silico analysis was performed.

RESULTS

Quinoline-6-carboxylic acid derivative compounds were designed and synthesized. The highest activity in the screening tests in the Rhabditis model was demonstrated by compound II-1 with a methoxyphenyl substituent LC50 = 0.3 mg/mL. In the next stage of the research, compound II-1 was analyzed in the H. contortus model. The results showed that compound II-1 was active and had ovicidal (percentage of dead eggs > 45 %) and larvicidal (percentage of dead larvae > 75 %) properties. Studies in the Strongylidae sp. model confirmed the ovicidal activity of compound II-1 (percentage of dead eggs ≥ 55 %). In vivo studies conducted in the H. polygyrus/bakeri nematode model showed that the number of nematodes decreased by an average of 30 % under the influence of compound II-1. In silico studies have shown two possible modes of action of compound II-1, i.e. inhibition of tubulin polymerization and SDH. The test compound did not show any systemic toxic effects. Its influence on drug metabolism related to the activity of cytochrome CYP450 enzymes was also investigated.

CONCLUSION

The results obtained in the in vitro, in vivo, and in silico studies indicate that the test compound can be described as a HIT, which in the future may be used in the treatment of parasitic diseases in humans and animals.

Challenges of Protein-Protein Docking of the Membrane Proteins

Despite the recent advances in the determination of high-resolution membrane protein (MP) structures, the structural and functional characterization of MPs remains extremely challenging, mainly due to the hydrophobic nature, low abundance, poor expression, purification, and crystallization difficulties associated with MPs. Whereby the major challenges/hurdles for MP structure determination are associated with the expression, purification, and crystallization procedures. Although there have been significant advances in the experimental determination of MP structures, only a limited number of MP structures (approximately less than 1% of all) are available in the Protein Data Bank (PDB). Therefore, the structures of a large number of MPs still remain unresolved, which leads to the availability of widely unplumbed structural and functional information related to MPs. As a result, recent developments in the drug discovery realm and the significant biological contemplation have led to the development of several novel, low-cost, and time-efficient computational methods that overcome the limitations of experimental approaches, supplement experiments, and provide alternatives for the characterization of MPs. Whereby the fine tuning and optimizations of these computational approaches remains an ongoing endeavor.Computational methods offer a potential way for the elucidation of structural features and the augmentation of currently available MP information. However, the use of computational modeling can be extremely challenging for MPs mainly due to insufficient knowledge of (or gaps in) atomic structures of MPs. Despite the availability of numerous in silico methods for 3D structure determination the applicability of these methods to MPs remains relatively low since all methods are not well-suited or adequate for MPs. However, sophisticated methods for MP structure predictions are constantly being developed and updated to integrate the modifications required for MPs. Currently, different computational methods for (1) MP structure prediction, (2) stability analysis of MPs through molecular dynamics simulations, (3) modeling of MP complexes through docking, (4) prediction of interactions between MPs, and (5) MP interactions with its soluble partner are extensively used. Towards this end, MP docking is widely used. It is notable that the MP docking methods yet few in number might show greater potential in terms of filling the knowledge gap. In this chapter, MP docking methods and associated challenges have been reviewed to improve the applicability, accuracy, and the ability to model macromolecular complexes.

Disrupted Choline Clearance and Sustained Acetylcholine Release In Vivo by a Common Choline Transporter Coding Variant Associated with Poor Attentional Control in Humans

Transport of choline via the neuronal high-affinity choline transporter (CHT; SLC5A7) is essential for cholinergic terminals to synthesize and release acetylcholine (ACh). In humans, we previously demonstrated an association between a common CHT coding substitution (rs1013940; Ile89Val) and reduced attentional control as well as attenuated frontal cortex activation. Here, we used a CRISPR/Cas9 approach to generate mice expressing the I89V substitution and assessed, in vivo, CHT-mediated choline transport, and ACh release. Relative to wild-type (WT) mice, CHT-mediated clearance of choline in male and female mice expressing one or two Val89 alleles was reduced by over 80% in cortex and over 50% in striatum. Choline clearance in CHT Val89 mice was further reduced by neuronal inactivation. Deficits in ACh release, 5 and 10 min after repeated depolarization at a low, behaviorally relevant frequency, support an attenuated reloading capacity of cholinergic neurons in mutant mice. The density of CHTs in total synaptosomal lysates and neuronal plasma-membrane-enriched fractions was not impacted by the Val89 variant, indicating a selective impact on CHT function. When challenged with a visual disruptor to reveal attentional control mechanisms, Val89 mice failed to adopt a more conservative response bias. Structural modeling revealed that Val89 may attenuate choline transport by altering conformational changes of CHT that support normal transport rates. Our findings support the view that diminished sustained cholinergic signaling capacity underlies perturbed attentional performance in individuals expressing CHT Val89. The CHT Val89 mouse serves as a valuable model to study heritable risk for cognitive disorders arising from cholinergic dysfunction.SIGNIFICANCE STATEMENT Acetylcholine (ACh) signaling depends on the functional capacity of the neuronal choline transporter (CHT). Previous research demonstrated that humans expressing the common CHT coding variant Val89 exhibit attentional vulnerabilities and attenuated fronto-cortical activation during attention. Here, we find that mice engineered to express the Val89 variant exhibit reduced CHT-mediated choline clearance and a diminished capacity to sustain ACh release. Additionally, Val89 mice lack cognitive flexibility in response to an attentional challenge. These findings provide a mechanistic and cognitive framework for interpreting the attentional phenotype associated with the human Val89 variant and establish a model that permits a more invasive interrogation of CNS effects as well as the development of therapeutic strategies for those, including Val89 carriers, with presynaptic cholinergic perturbations.

Structural Insights into the Substrate Transport Mechanisms in GTR Transporters through Ensemble Docking

Glucosinolate transporters (GTRs) are part of the nitrate/peptide transporter (NPF) family, members of which also transport specialized secondary metabolites as substrates. Glucosinolates are defense compounds derived from amino acids. We selected 4-methylthiobutyl (4MTB) and indol-3-ylmethyl (I3M) glucosinolates to study how GTR1 from Arabidopsis thaliana transports these substrates in computational simulation approaches. The designed pipeline reported here includes massive docking of 4MTB and I3M in an ensemble of GTR1 conformations (in both inward and outward conformations) extracted from molecular dynamics simulations, followed by clustered and substrate–protein interactions profiling. The identified key residues were mutated, and their role in substrate transport was tested. We were able to identify key residues that integrate a major binding site of these substrates, which is critical for transport activity. In silico approaches employed here represent a breakthrough in the plant transportomics field, as the identification of key residues usually takes a long time if performed from a purely wet-lab experimental perspective. The inclusion of structural bioinformatics in the analyses of plant transporters significantly speeds up the knowledge-gaining process and optimizes valuable time and resources.

Chloride-dependent conformational changes in the GlyT1 glycine transporter

The human GlyT1 glycine transporter requires chloride for its function. However, the mechanism by which Cl^- exerts its influence is unknown. To examine the role that Cl^- plays in the transport cycle, we measured the effect of Cl^- on both glycine binding and conformational changes. The ability of glycine to displace the high-affinity radioligand [³H]CHIBA-3007 required Na⁺ and was potentiated over 1,000-fold by Cl^- We generated GlyT1b mutants containing reactive cysteine residues in either the extracellular or cytoplasmic permeation pathways and measured changes in the reactivity of those cysteine residues as indicators of conformational changes in response to ions and substrate. Na⁺ increased accessibility in the extracellular pathway and decreased it in the cytoplasmic pathway, consistent with stabilizing an outward-open conformation as observed in other members of this transporter family. In the presence of Na⁺, both glycine and Cl^- independently shifted the conformation of GlyT1b toward an outward-closed conformation. Together, Na⁺, glycine, and Cl^- stabilized an inward-open conformation of GlyT1b. We then examined whether Cl^- acts by interacting with a conserved glutamine to allow formation of an ion pair that stabilizes the closed state of the extracellular pathway. Molecular dynamics simulations of a GlyT1 homolog indicated that this ion pair is formed more frequently as that pathway closes. Mutation of the glutamine blocked the effect of Cl^-, and substituting it with glutamate or lysine resulted in outward- or inward-facing transporter conformations, respectively. These results provide an unexpected insight into the role of Cl^- in this family of transporters.

Structural Determinants of the Dopamine Transporter Regulation Mediated by G Proteins

Dopamine clearance in the brain is controlled by the dopamine transporter (DAT), a protein residing in the plasma membrane, which drives reuptake of extracellular dopamine into presynaptic neurons. Studies have revealed that the βγ subunits of heterotrimeric G proteins modulate DAT function through a physical association with the C-terminal region of the transporter. Regulation of neurotransmitter transporters by Gβγ subunits is unprecedented in the literature; therefore, it is interesting to investigate the structural details of this particular protein-protein interaction. Here, we refined the crystal structure of the Drosophila melanogaster DAT (dDAT), modeling de novo the N- and C-terminal domains; subsequently, we used the full-length dDAT structure to generate a comparative model of human DAT (hDAT). Both proteins were assembled with Gβ1γ2 subunits employing protein-protein docking, and subsequent molecular dynamics simulations were run to identify the specific interactions governing the formation of the hDAT:Gβγ and dDAT:Gβγ complexes. A [L/F]R[Q/E]R sequence motif containing the residues R588 in hDAT and R587 in dDAT was found as key to bind the Gβγ subunits through electrostatic interactions with a cluster of negatively charged residues located at the top face of the Gβ subunit. Alterations of DAT function have been associated with multiple devastating neuropathological conditions; therefore, this work represents a step toward better understanding DAT regulation by signaling proteins, allowing us to predict therapeutic target regions.

Studies of structural determinants of substrate binding in the Creatine Transporter (CreaT, SLC6A8) using molecular models

Creatine is a crucial metabolite that plays a fundamental role in ATP homeostasis in tissues with high-energy demands. The creatine transporter (CreaT, SLC6A8) belongs to the solute carrier 6 (SLC6) transporters family, and more particularly to the GABA transporters (GATs) subfamily. Understanding the molecular determinants of specificity within the SLC6 transporters in general, and the GATs in particular is very challenging due to the high similarity of these proteins. In the study presented here, our efforts focused on finding key structural features involved in binding selectivity for CreaT using structure-based computational methods. Due to the lack of three-dimensional structures of SLC6A8, our approach was based on the realization of two reliable homology models of CreaT using the structures of two templates, i.e. the human serotonin transporter (hSERT) and the prokaryotic leucine transporter (LeuT). Our models reveal that an optimal complementarity between the shape of the binding site and the size of the ligands is necessary for transport. These findings provide a framework for a deeper understanding of substrate selectivity of the SLC6 family and other LeuT fold transporters.

Identification of a crucial amino acid implicated in the hydroxylation/desaturation ratio of CpFAH12 bifunctional hydroxylase

Claviceps purpurea bifunctional Δ12-hydroxylase/desaturase, CpFAH12, and monofunctional desaturase CpFAD2, share 86% of sequence identity. To identify the underlying determinants of the hydroxylation/desaturation specificity, chimeras of these two enzymes were tested for their fatty acid production in an engineered Yarrowia lipolytica strain. It reveals that transmembrane helices are not involved in the hydroxylation/desaturation specificity whereas all cytosolic domains have an impact on it. Especially, replacing the CpFAH12 cytosolic part near the second histidine-box by the corresponding CpFAD2 part annihilates all hydroxylation activity. Further mutagenesis experiments within this domain identified isoleucine 198 as the crucial element for the hydroxylation activity of CpFAH12. Monofunctional variants performing the only desaturation were obtained when this position was exchanged by the threonine of CpFAD2. Saturation mutagenesis at this position showed modulation in the hydroxylation/desaturation specificity in the different variants. The WT enzyme was demonstrated as the most efficient for ricinoleic acid production and some variants showed a better desaturation activity. A model based on the recently discovered membrane desaturase structures indicate that these changes in specificity are more likely due to modifications in the di-iron center geometry rather than changes in the substrate binding mode.

A structural model of the human serotonin transporter in an outward-occluded state

The human serotonin transporter hSERT facilitates the reuptake of its endogenous substrate serotonin from the synaptic cleft into presynaptic neurons after signaling. Reuptake regulates the availability of this neurotransmitter and therefore hSERT plays an important role in balancing human mood conditions. In 2016, the first 3D structures of this membrane transporter were reported in an inhibitor-bound, outward-open conformation. These structures revealed valuable information about interactions of hSERT with antidepressant drugs. Nevertheless, the question remains how serotonin facilitates the specific conformational changes that open and close pathways from the synapse and to the cytoplasm as required for transport. Here, we present a serotonin-bound homology model of hSERT in an outward-occluded state, a key intermediate in the physiological cycle, in which the interactions with the substrate are likely to be optimal. Our approach uses two template structures and includes careful refinement and comprehensive computational validation. According to microsecond-long molecular dynamics simulations, this model exhibits interactions between the gating residues in the extracellular pathway, and these interactions differ from those in an outward-open conformation of hSERT bound to serotonin. Moreover, we predict several features of this state by monitoring the intracellular gating residues, the extent of hydration, and, most importantly, protein-ligand interactions in the central binding site. The results illustrate common and distinct characteristics of these two transporter states and provide a starting point for future investigations of the transport mechanism in hSERT.

A Structural Model of the Inactivation Gate of Voltage-Activated Potassium Channels

After opening, the Shaker voltage-gated potassium (K_V) channel rapidly inactivates when one of its four N-termini enters and occludes the channel pore. Although it is known that the tip of the N-terminus reaches deep into the central cavity, the conformation adopted by this domain during inactivation and the nature of its interactions with the rest of the channel remain unclear. Here, we use molecular dynamics simulations coupled with electrophysiology experiments to reveal the atomic-scale mechanisms of inactivation. We find that the first six amino acids of the N-terminus spontaneously enter the central cavity in an extended conformation, establishing hydrophobic contacts with residues lining the pore. A second portion of the N-terminus, consisting of a long 24 amino acid α-helix, forms numerous polar contacts with residues in the intracellular entryway of the T1 domain. Double mutant cycle analysis revealed a strong relationship between predicted interatomic distances and empirically observed thermodynamic coupling, establishing a plausible model of the transition of K_V channels to the inactivated state.

Current Advances in Studying Clinically Relevant Transporters of the Solute Carrier (SLC) Family by Connecting Computational Modeling and Data Science

Organic anion and cation transporting proteins (OATs, OATPs, and OCTs), as well as the Multidrug and Toxin Extrusion (MATE) transporters of the Solute Carrier (SLC) family are playing a pivotal role in the discovery and development of new drugs due to their involvement in drug disposition, drug-drug interactions, adverse drug effects and related toxicity. Computational methods to understand and predict clinically relevant transporter interactions can provide useful guidance at early stages in drug discovery and design, especially if they include contemporary data science approaches. In this review, we summarize the current state-of-the-art of computational approaches for exploring ligand interactions and selectivity for these drug (uptake) transporters. The computational methods discussed here by highlighting interesting examples from the current literature are ranging from semiautomatic data mining and integration, to ligand-based methods (such as quantitative structure-activity relationships, and combinatorial pharmacophore modeling), and finally structure-based methods (such as comparative modeling, molecular docking, and molecular dynamics simulations). We are focusing on promising computational techniques such as fold-recognition methods, proteochemometric modeling or techniques for enhanced sampling of protein conformations used in the context of these ADMET-relevant SLC transporters with a special focus on methods useful for studying ligand selectivity.

A Membrane Burial Potential with H-Bonds and Applications to Curved Membranes and Fast Simulations

We use the statistics of a large and curated training set of transmembrane helical proteins to develop a knowledge-based potential that accounts for the dependence on both the depth of burial of the protein in the membrane and the degree of side-chain exposure. Additionally, the statistical potential includes depth-dependent energies for unsatisfied backbone hydrogen bond donors and acceptors, which are found to be relatively small, ∼2 RT. Our potential accurately places known proteins within the bilayer. The potential is applied to the mechanosensing MscL channel in membranes of varying thickness and curvature, as well as to the prediction of protein structure. The potential is incorporated into our new Upside molecular dynamics algorithm. Notably, we account for the exchange of protein-lipid interactions for protein-protein interactions as helices contact each other, thereby avoiding overestimating the energetics of helix association within the membrane. Simulations of most multimeric complexes find that isolated monomers and the oligomers retain the same orientation in the membrane, suggesting that the assembly of prepositioned monomers presents a viable mechanism of oligomerization.

Structural models of the NaPi-II sodium-phosphate cotransporters

Progress towards understanding the molecular mechanisms of phosphate homeostasis through sodium-dependent transmembrane uptake has long been stymied by the absence of structural information about the NaPi-II sodium-phosphate transporters. For many other coupled transporters, even those unrelated to NaPi-II, internal repeated elements have been revealed as a key feature that is inherent to their function. Here, we review recent structure prediction studies for NaPi-II transporters. Attempts to identify structural templates for NaPi-II transporters have leveraged the structural repeat perspective to uncover an otherwise obscured relationship with the dicarboxylate-sodium symporters (DASS). This revelation allowed the prediction of three-dimensional structural models of human NaPi-IIa and flounder NaPi-IIb, whose folds were evaluated by comparison with available biochemical data outlining the transmembrane topology and solvent accessibility of various regions of the protein. Using these structural models, binding sites for sodium and phosphate were proposed. The predicted sites were tested and refined based on detailed electrophysiological and biochemical studies and were validated by comparison with subsequently reported structures of transporters belonging to the AbgT family. Comparison with the DASS transporter VcINDY suggested a conformational mechanism involving a large, two-domain structural change, known as an elevator-like mechanism. These structural models provide a foundation for further studies into substrate binding, conformational change, kinetics, and energetics of sodium-phosphate transport. We discuss future opportunities, as well as the challenges that remain.

Structural relationship between the putative hair cell mechanotransduction channel TMC1 and TMEM16 proteins

The hair cell mechanotransduction (MET) channel complex is essential for hearing, yet it's molecular identity and structure remain elusive. The transmembrane channel-like 1 (TMC1) protein localizes to the site of the MET channel, interacts with the tip-link responsible for mechanical gating, and genetic alterations in TMC1 alter MET channel properties and cause deafness, supporting the hypothesis that TMC1 forms the MET channel. We generated a model of TMC1 based on X-ray and cryo-EM structures of TMEM16 proteins, revealing the presence of a large cavity near the protein-lipid interface that also harbors the Beethoven mutation, suggesting that it could function as a permeation pathway. We also find that hair cells are permeable to 3 kDa dextrans, and that dextran permeation requires TMC1/2 proteins and functional MET channels, supporting the presence of a large permeation pathway and the hypothesis that TMC1 is a pore forming subunit of the MET channel complex.

The ins and outs of vesicular monoamine transporters

Yaffe et al. review structure-guided studies that have provided insight into the mechanism of proton-monoamine antiport by VMATs.

The H⁺-coupled vesicular monoamine transporter (VMAT) is a transporter essential for life. VMAT mediates packaging of the monoamines serotonin, dopamine, norepinephrine, and histamine from the neuronal cytoplasm into presynaptic vesicles, which is a key step in the regulated release of neurotransmitters. However, a detailed understanding of the mechanism of VMAT function has been limited by the lack of availability of high-resolution structural data. In recent years, a series of studies guided by homology models has revealed significant insights into VMAT function, identifying residues that contribute to the binding site and to specific steps in the transport cycle. Moreover, to characterize the conformational transitions that occur upon binding of the substrate and coupling ion, we have taken advantage of the unique and powerful pharmacology of VMAT as well as of mutants that affect the conformational equilibrium of the protein and shift it toward defined conformations. This has allowed us to identify an important role for the proton gradient in driving a shift from lumen-facing to cytoplasm-facing conformations.

Iodide Binding in Sodium-Coupled Cotransporters

Several apical iodide translocation pathways have been proposed for iodide efflux out of thyroid follicular cells, including a pathway mediated by the sodium-coupled monocarboxylate transporter 1 (SMCT1), which remains controversial. Herein, we evaluate structural and functional similarities between SMCT1 and the well-studied sodium-iodide symporter (NIS) that mediates the first step of iodide entry into the thyroid. Free-energy calculations using a force field with electronic polarizability verify the presence of a conserved iodide-binding pocket between the TM2, TM3, and TM7 segments in hNIS, where iodide is coordinated by Phe67, Gln72, Cys91, and Gln94. We demonstrate the mutation of residue Gly93 of hNIS to a larger amino acid expels the side chain of a critical tryptophan residue (Trp255) into the interior of the binding pocket, partially occluding the iodide binding site and reducing iodide affinity, which is consistent with previous reports associating mutation of this residue with iodide uptake deficiency and hypothyroidism. Furthermore, we find that the position of Trp255 in this hNIS mutant mirrors that of Trp253 in wild-type hSMCT1, where a threonine (Thr91) occupies the position homologous to that occupied by glycine in wild-type hNIS (Gly93). Correspondingly, mutation of Thr91 to glycine in hSMCT1 makes the pocket structure more like that of wild-type hNIS, increasing its iodide affinity. These results suggest that wild-type hSMCT1 in the inward-facing conformation may bind iodide only very weakly, which may have implications for its ability to transport iodide.

Cation Interactions and Membrane Potential Induce Conformational Changes in NaPi-IIb

Voltage-dependence of Na(+)-coupled phosphate cotransporters of the SLC34 family arises from displacement of charges intrinsic to the protein and the binding/release of one Na(+) ion in response to changes in the transmembrane electric field. Candidate coordination residues for the cation at the Na1 site were previously predicted by structural modeling using the x-ray structure of dicarboxylate transporter VcINDY as template and confirmed by functional studies. Mutations at Na1 resulted in altered steady-state and presteady-state characteristics that should be mirrored in the conformational changes induced by membrane potential changes. To test this hypothesis by functional analysis, double mutants of the flounder SLC34A2 protein were constructed that contain one of the Na1-site perturbing mutations together with a substituted cysteine for fluorophore labeling, as expressed in Xenopus oocytes. The locations of the mutations were mapped onto a homology model of the flounder protein. The effects of the mutagenesis were characterized by steady-state, presteady-state, and fluorometric assays. Changes in fluorescence intensity (ΔF) in response to membrane potential steps were resolved at three previously identified positions. These fluorescence data corroborated the altered presteady-state kinetics upon perturbation of Na1, and furthermore indicated concomitant changes in the microenvironment of the respective fluorophores, as evidenced by changes in the voltage dependence and time course of ΔF. Moreover, iodide quenching experiments indicated that the aqueous nature of the fluorophore microenvironment depended on the membrane potential. These findings provide compelling evidence that membrane potential and cation interactions induce significant large-scale structural rearrangements of the protein.

Exploring the binding properties and structural stability of an opsin in the chytrid Spizellomyces punctatus using comparative and molecular modeling

Background

Opsin proteins are seven transmembrane receptor proteins which detect light. Opsins can be classified into two types and share little sequence identity: type 1, typically found in bacteria, and type 2, primarily characterized in metazoa. The type 2 opsins (Rhodopsins) are a subfamily of G-protein coupled receptors (GPCRs), a large and diverse class of seven transmembrane proteins and are generally restricted to metazoan lineages. Fungi use light receptors including opsins to sense the environment and transduce signals for developmental or metabolic changes. Opsins characterized in the Dikarya (Ascomycetes and Basidiomycetes) are of the type 1 bacteriorhodopsin family but the early diverging fungal lineages have not been as well surveyed. We identified by sequence similarity a rhodopsin-like GPCR in genomes of early diverging chytrids and examined the structural characteristics of this protein to assess its likelihood to be homologous to animal rhodopsins and bind similar chromophores.

Methods

We used template-based structure modeling, automated ligand docking, and molecular modeling to assess the structural and binding properties of an identified opsin-like protein found in Spizellomyces punctatus, a unicellular, flagellated species belonging to Chytridiomycota, one of the earliest diverging fungal lineages. We tested if the sequence and inferred structure were consistent with a solved crystal structure of a type 2 rhodopsin from the squid Todarodes pacificus.

Results

Our results indicate that the Spizellomyces opsin has structural characteristics consistent with functional animal type 2 rhodopsins and is capable of maintaining a stable structure when associated with the retinaldehyde chromophore, specifically the 9-cis-retinal isomer. Together, these results support further the homology of Spizellomyces opsins to animal type 2 rhodopsins.

Discussion

This represents the first test of structure/function relationship of a type 2 rhodopsin identified in early branching fungal lineages, and provides a foundation for future work exploring pathways and components of photoreception in early fungi.

An Extra Amino Acid Residue in Transmembrane Domain 10 of the γ-Aminobutyric Acid (GABA) Transporter GAT-1 Is Required for Efficient Ion-coupled Transport

The GABA transporter GAT-1 mediates electrogenic transport of its substrate together with sodium and chloride. It is a member of the neurotransmitter:sodium:symporters, which are crucial for synaptic transmission. Compared with all other neurotransmitter:sodium:symporters, GAT-1 and other members of the GABA transporter subfamily all contain an extra amino acid residue at or near a conserved glycine in transmembrane segment 10. Therefore, we studied the functional impact of deletion and replacement mutants of Gly-457 and its two adjacent residues in GAT-1. The glycine replacement mutants were devoid of transport activity, but remarkably the deletion mutant was active, as were mutants obtained by deleting positions on either side of Gly-457. However, the inward rectification of GABA-induced transport currents by all three deletion mutants was diminished, and the charge-to-flux ratio was increased by more than 2.5-fold, both of which indicate substantial uncoupled transport. These observations suggest that the deletions render the transporters less tightly packed. Consistent with this interpretation, the inactive G457A mutant was partially rescued by removing the adjacent serine residue. Moreover, the activity of several gating mutants was also partially rescued upon deletion of Gly-457. Structural modeling showed that the stretch surrounding Gly-457 is likely to form a π-helix. Our data indicate that the "extra" residue in transmembrane domain 10 of the GABA transporter GAT-1 provides extra bulk, probably in the form of a π-helix, which is required for stringent gating and tight coupling of ion and substrate fluxes in the GABA transporter family.

Transmembrane topology of the arsenite permease Acr3 from Saccharomyces cerevisiae

Acr3 is a plasma membrane transporter, a member of the bile/arsenite/riboflavin transporter (BART) superfamily, which confers high-level resistance to arsenicals in the yeast Saccharomyces cerevisiae. We have previously shown that the yeast Acr3 acts as a low affinity As(III)/H⁺ and Sb(III)/H⁺ antiporter. We have also identified several amino acid residues that are localized in putative transmembrane helices (TM) and appeared to be critical for the Acr3 activity. In the present study, the topology of Acr3 was investigated by insertion of glycosylation and factor Xa protease cleavage sites at predicted hydrophilic regions. The analysis of the glycosylation pattern and factor Xa cleavage products of resulting Acr3 fusion constructs provide evidence supporting a topological model of Acr3 with 10 TM segments and cytoplasmically oriented N- and C-terminal domains. Next, we investigated the role of the hydrophilic loop connecting TM8 and TM9, the large size of which is unique to members of the yeast Acr3 family of metalloid transporters. We found that a 28 amino acid deletion in this region does not affect Acr3 folding, trafficking substrate binding, or transport activity. Finally, we constructed a homology-based structural model of Acr3 using the crystal structure of the Yersinia frederiksenii homologue of the human bile acid sodium symporter ASBT.

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.

SLC Transporters: Structure, Function, and Drug Discovery

The human Solute Carrier (SLC) transporters are important targets for drug development. Structure-based drug discovery for SLC transporters requires the description of their structure, dynamics, and mechanism of interaction with small molecule ligands and ions. The recent determination of atomic structures of human SLC transporters and their homologs, combined with improved computational power and prediction methods have led to an increased applicability of structure-based drug design methods for human SLC members. In this review, we provide an overview of the SLC transporters' structures and transport mechanisms. We then describe computational techniques, such as homology modeling and virtual screening that are emerging as key tools to discover chemical probes for human SLC members. We illustrate the utility of these methods by presenting case studies in which rational integration of computation and experiment was used to characterize SLC members that transport key nutrients and metabolites, including the amino acid transporters LAT-1 and ASCT2, the SLC13 family of citric acid cycle intermediate transporters, and the glucose transporter GLUT1. We conclude with a brief discussion about future directions in structure-based drug discovery for the human SLC superfamily, one of the most structurally and functionally diverse protein families in human.

Mechanism of Paroxetine (Paxil) Inhibition of the Serotonin Transporter

The serotonin transporter (SERT) is an integral membrane protein that exploits preexisting sodium-, chloride-, and potassium ion gradients to catalyze the thermodynamically unfavorable movement of synaptic serotonin into the presynaptic neuron. SERT has garnered significant clinical attention partly because it is the target of multiple psychoactive agents, including the antidepressant paroxetine (Paxil), the most potent selective serotonin reuptake inhibitor known. However, the binding site and orientation of paroxetine in SERT remain controversial. To provide molecular insight, we constructed SERT homology models based on the Drosophila melanogaster dopamine transporter and docked paroxetine to these models. We tested the predicted binding configurations with a combination of radioligand binding and flux assays on wild-type and mutant SERTs. Our data suggest that the orientation of paroxetine, specifically its fluorophenyl ring, in SERT’s substrate binding site directly depends on this pocket’s charge distribution, and thereby provide an avenue toward understanding and enhancing high-affinity antidepressant activity.

Structure-Activity Relationship Studies of N- and C-Terminally Modified Secretin Analogs for the Human Secretin Receptor

The pleiotropic role of human secretin (hSCT) validates its potential use as a therapeutic agent. Nevertheless, the structure of secretin in complex with its receptor is necessary to develop a suitable therapeutic agent. Therefore, in an effort to design a three-dimensional virtual homology model and identify a peptide agonist and/or antagonist for the human secretin receptor (hSR), the significance of the primary sequence of secretin peptides in allosteric binding and activation was elucidated using virtual docking, FRET competitive binding and assessment of the cAMP response. Secretin analogs containing various N- or C-terminal modifications were prepared based on previous findings of the role of these domains in receptor binding and activation. These analogs exhibited very low or no binding affinity in a virtual model, and were found to neither exhibit in vitro binding nor agonistic or antagonistic properties. A parallel analysis of the analogs in the virtual model and in vitro studies revealed instability of these peptide analogs to bind and activate the receptor.

Identification of the first sodium binding site of the phosphate cotransporter NaPi-IIa (SLC34A1)

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.

Structural insights into the mechanism of activation of the TRPV1 channel by a membrane-bound tarantula toxin

Venom toxins are invaluable tools for exploring the structure and mechanisms of ion channels. Here, we solve the structure of double-knot toxin (DkTx), a tarantula toxin that activates the heat-activated TRPV1 channel. We also provide improved structures of TRPV1 with and without the toxin bound, and investigate the interactions of DkTx with the channel and membranes. We find that DkTx binds to the outer edge of the external pore of TRPV1 in a counterclockwise configuration, using a limited protein-protein interface and inserting hydrophobic residues into the bilayer. We also show that DkTx partitions naturally into membranes, with the two lobes exhibiting opposing energetics for membrane partitioning and channel activation. Finally, we find that the toxin disrupts a cluster of hydrophobic residues behind the selectivity filter that are critical for channel activation. Collectively, our findings reveal a novel mode of toxin-channel recognition that has important implications for the mechanism of thermosensation.

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.

Insight into the Meligethes aeneus voltage-sensitive sodium channel structure and an attempt to select the best pyrethroid ligands

BACKGROUND

The voltage-sensitive sodium channel (VSSC) is a target for the pharmacological action of pyrethroids which are used in controlling pests, including those of agricultural importance. Among these is the pollen beetle (Meligethes aeneus F.) - the most serious pest of Brassica napus. Owing to the heavy use of pyrethroids, a widespread build-up of resistance has occurred. The main cause of pyrethroid insensitivity in M. aeneus is considered to be an increased oxidative metabolism; however, the additional mechanism of resistance associated with mutations in the VSSC might contribute to this phenomenon.

RESULTS

We generated a VSSC 3D model to study the docking affinities of pyrethroids to their target site within the channel. Our goal was to identify the pyrethroids for which docking affinity scores were high and not affected by potential mutations in the VSSC. We found that the docking scores of cypermethrin are hardly influenced by the appearance of point mutations. Additionally, tau-fluvalinate, deltamethrin and bifenthrin are VSSC ligands with high affinity scores.

CONCLUSIONS

Our docking models suggest that point mutations in the VSSC binding pocket might affect the stability of ligand interactions and change the pattern of ligand docking locations, which might have a potential effect on VSSC gating properties.

Repeat-swap homology modeling of secondary active transporters: updated protocol and prediction of elevator-type mechanisms

Secondary active transporters are critical for neurotransmitter clearance and recycling during synaptic transmission and uptake of nutrients. These proteins mediate the movement of solutes against their concentration gradients, by using the energy released in the movement of ions down pre-existing concentration gradients. To achieve this, transporters conform to the so-called alternating-access hypothesis, whereby the protein adopts at least two conformations in which the substrate binding sites are exposed to one or other side of the membrane, but not both simultaneously. Structures of a bacterial homolog of neuronal glutamate transporters, Glt_Ph, in several different conformational states have revealed that the protein structure is asymmetric in the outward- and inward-open states, and that the conformational change connecting them involves a elevator-like movement of a substrate binding domain across the membrane. The structural asymmetry is created by inverted-topology repeats, i.e., structural repeats with similar overall folds whose transmembrane topologies are related to each other by two-fold pseudo-symmetry around an axis parallel to the membrane plane. Inverted repeats have been found in around three-quarters of secondary transporter folds. Moreover, the (a)symmetry of these systems has been successfully used as a bioinformatic tool, called “repeat-swap modeling” to predict structural models of a transporter in one conformation using the known structure of the transporter in the complementary conformation as a template. Here, we describe an updated repeat-swap homology modeling protocol, and calibrate the accuracy of the method using Glt_Ph, for which both inward- and outward-facing conformations are known. We then apply this repeat-swap homology modeling procedure to a concentrative nucleoside transporter, VcCNT, which has a three-dimensional arrangement related to that of Glt_Ph. The repeat-swapped model of VcCNT predicts that nucleoside transport also occurs via an elevator-like mechanism.

TMFoldRec: a statistical potential-based transmembrane protein fold recognition tool

Background

Transmembrane proteins (TMPs) are the key components of signal transduction, cell-cell adhesion and energy and material transport into and out from the cells. For the deep understanding of these processes, structure determination of transmembrane proteins is indispensable. However, due to technical difficulties, only a few transmembrane protein structures have been determined experimentally. Large-scale genomic sequencing provides increasing amounts of sequence information on the proteins and whole proteomes of living organisms resulting in the challenge of bioinformatics; how the structural information should be gained from a sequence.

Results

Here, we present a novel method, TMFoldRec, for fold prediction of membrane segments in transmembrane proteins. TMFoldRec based on statistical potentials was tested on a benchmark set containing 124 TMP chains from the PDBTM database. Using a 10-fold jackknife method, the native folds were correctly identified in 77 % of the cases. This accuracy overcomes the state-of-the-art methods. In addition, a key feature of TMFoldRec algorithm is the ability to estimate the reliability of the prediction and to decide with an accuracy of 70 %, whether the obtained, lowest energy structure is the native one.

Conclusion

These results imply that the membrane embedded parts of TMPs dictate the TM structures rather than the soluble parts. Moreover, predictions with reliability scores make in this way our algorithm applicable for proteome-wide analyses.

Availability

The program is available upon request for academic use.

Electronic supplementary material

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

Sequence-structure relationship study in all-α transmembrane proteins using an unsupervised learning approach

Transmembrane proteins (TMPs) are major drug targets, but the knowledge of their precise topology structure remains highly limited compared with globular proteins. In spite of the difficulties in obtaining their structures, an important effort has been made these last years to increase their number from an experimental and computational point of view. In view of this emerging challenge, the development of computational methods to extract knowledge from these data is crucial for the better understanding of their functions and in improving the quality of structural models. Here, we revisit an efficient unsupervised learning procedure, called Hybrid Protein Model (HPM), which is applied to the analysis of transmembrane proteins belonging to the all-α structural class. HPM method is an original classification procedure that efficiently combines sequence and structure learning. The procedure was initially applied to the analysis of globular proteins. In the present case, HPM classifies a set of overlapping protein fragments, extracted from a non-redundant databank of TMP 3D structure. After fine-tuning of the learning parameters, the optimal classification results in 65 clusters. They represent at best similar relationships between sequence and local structure properties of TMPs. Interestingly, HPM distinguishes among the resulting clusters two helical regions with distinct hydrophobic patterns. This underlines the complexity of the topology of these proteins. The HPM classification enlightens unusual relationship between amino acids in TMP fragments, which can be useful to elaborate new amino acids substitution matrices. Finally, two challenging applications are described: the first one aims at annotating protein functions (channel or not), the second one intends to assess the quality of the structures (X-ray or models) via a new scoring function deduced from the HPM classification.

Assessing the local structural quality of transmembrane protein models using statistical potentials (QMEANBrane)

Motivation: Membrane proteins are an important class of biological macromolecules involved in many cellular key processes including signalling and transport. They account for one third of genes in the human genome and >50% of current drug targets. Despite their importance, experimental structural data are sparse, resulting in high expectations for computational modelling tools to help fill this gap. However, as many empirical methods have been trained on experimental structural data, which is biased towards soluble globular proteins, their accuracy for transmembrane proteins is often limited.

Results: We developed a local model quality estimation method for membrane proteins (‘QMEANBrane’) by combining statistical potentials trained on membrane protein structures with a per-residue weighting scheme. The increasing number of available experimental membrane protein structures allowed us to train membrane-specific statistical potentials that approach statistical saturation. We show that reliable local quality estimation of membrane protein models is possible, thereby extending local quality estimation to these biologically relevant molecules.

Availability and implementation: Source code and datasets are available on request.

Contact:torsten.schwede@unibas.ch

Supplementary Information:Supplementary data are available at Bioinformatics online.

Computational modeling of membrane proteins

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1-2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug-specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans-membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α-helical MPs as well as β-barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge-based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade.

Improving homology modeling of G-protein coupled receptors through multiple-template derived conserved inter-residue interactions

Evidenced by the three-rounds of G-protein coupled receptors (GPCR) Dock competitions, improving homology modeling methods of helical transmembrane proteins including the GPCRs, based on templates of low sequence identity, remains an eminent challenge. Current approaches addressing this challenge adopt the philosophy of "modeling first, refinement next". In the present work, we developed an alternative modeling approach through the novel application of available multiple templates. First, conserved inter-residue interactions are derived from each additional template through conservation analysis of each template-target pairwise alignment. Then, these interactions are converted into distance restraints and incorporated in the homology modeling process. This approach was applied to modeling of the human β2 adrenergic receptor using the bovin rhodopsin and the human protease-activated receptor 1 as templates and improved model quality was demonstrated compared to the homology model generated by standard single-template and multiple-template methods. This method of "refined restraints first, modeling next", provides a fast and complementary way to the current modeling approaches. It allows rational identification and implementation of additional conserved distance restraints extracted from multiple templates and/or experimental data, and has the potential to be applicable to modeling of all helical transmembrane proteins.

Functionally important carboxyls in a bacterial homologue of the vesicular monoamine transporter (VMAT)

Transporters essential for neurotransmission in mammalian organisms and bacterial multidrug transporters involved in antibiotic resistance are evolutionarily related. To understand in more detail the evolutionary aspects of the transformation of a bacterial multidrug transporter to a mammalian neurotransporter and to learn about mechanisms in a milieu amenable for structural and biochemical studies, we identified, cloned, and partially characterized bacterial homologues of the rat vesicular monoamine transporter (rVMAT2). We performed preliminary biochemical characterization of one of them, Brevibacillus brevis monoamine transporter (BbMAT), from the bacterium B. brevis. BbMAT shares substrates with rVMAT2 and transports them in exchange with >1H(+), like the mammalian transporter. Here we present a homology model of BbMAT that has the standard major facilitator superfamily fold; that is, with two domains of six transmembrane helices each, related by 2-fold pseudosymmetry whose axis runs normal to the membrane and between the two halves. The model predicts that four carboxyl residues, a histidine, and an arginine are located in the transmembrane segments. We show here that two of the carboxyls are conserved, equivalent to the corresponding ones in rVMAT2, and are essential for H(+)-coupled transport. We conclude that BbMAT provides an excellent experimental paradigm for the study of its mammalian counterparts and bacterial multidrug transporters.

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.

In silico modeling of human α2C-adrenoreceptor interaction with filamin-2

Vascular smooth muscle α2C-adrenoceptors (α2C-ARs) mediate vasoconstriction of small blood vessels, especially arterioles. Studies of endogenous receptors in human arteriolar smooth muscle cells (referred to as microVSM) and transiently transfected receptors in heterologous HEK293 cells show that the α2C-ARs are perinuclear receptors that translocate to the cell surface under cellular stress and elicit a biological response. Recent studies in microVSM unraveled a crucial role of Rap1A-Rho-ROCK-F-actin pathways in receptor translocation, and identified protein-protein interaction of α2C-ARs with the actin binding protein filamin-2 as an essential step in the process. To better understand the molecular nature and specificity of this interaction, in this study, we constructed comparative models of human α2C-AR and human filamin-2 proteins. Finally, we performed in silico protein-protein docking to provide a structural platform for the investigation of human α2C-AR and filamin-2 interactions. We found that electrostatic interactions seem to play a key role in this complex formation which manifests in interactions between the C-terminal arginines of α2C-ARs (particularly R454 and R456) and negatively charged residues from filamin-2 region between residues 1979 and 2206. Phylogenetic and sequence analysis showed that these interactions have evolved in warm-blooded animals.

Structural and functional analysis of a platelet-activating lysophosphatidylcholine of Trypanosoma cruzi

Background

Trypanosoma cruzi is the causative agent of the life-threatening Chagas disease, in which increased platelet aggregation related to myocarditis is observed. Platelet-activating factor (PAF) is a potent intercellular lipid mediator and second messenger that exerts its activity through a PAF-specific receptor (PAFR). Previous data from our group suggested that T. cruzi synthesizes a phospholipid with PAF-like activity. The structure of T. cruzi PAF-like molecule, however, remains elusive.

Methodology/Principal findings

Here, we have purified and structurally characterized the putative T. cruzi PAF-like molecule by electrospray ionization-tandem mass spectrometry (ESI-MS/MS). Our ESI-MS/MS data demonstrated that the T. cruzi PAF-like molecule is actually a lysophosphatidylcholine (LPC), namely sn-1 C18:1(delta 9)-LPC. Similar to PAF, the platelet-aggregating activity of C18:1-LPC was abrogated by the PAFR antagonist, WEB 2086. Other major LPC species, i.e., C16:0-, C18:0-, and C18:2-LPC, were also characterized in all T. cruzi stages. These LPC species, however, failed to induce platelet aggregation. Quantification of T. cruzi LPC species by ESI-MS revealed that intracellular amastigote and trypomastigote forms have much higher levels of C18:1-LPC than epimastigote and metacyclic trypomastigote forms. C18:1-LPC was also found to be secreted by the parasite in extracellular vesicles (EV) and an EV-free fraction. A three-dimensional model of PAFR was constructed and a molecular docking study was performed to predict the interactions between the PAFR model and PAF, and each LPC species. Molecular docking data suggested that, contrary to other LPC species analyzed, C18:1-LPC is predicted to interact with the PAFR model in a fashion similar to PAF.

Conclusions/Significance

Taken together, our data indicate that T. cruzi synthesizes a bioactive C18:1-LPC, which aggregates platelets via PAFR. We propose that C18:1-LPC might be an important lipid mediator in the progression of Chagas disease and its biosynthesis could eventually be exploited as a potential target for new therapeutic interventions.

Chagas disease, caused by the parasite Trypanosoma cruzi, was exclusively confined to Latin America but it has recently spread to other regions of the world. Chagas disease affects 8–10 million people and kills thousands of them every year. Lysophosphatidylcholine (LPC) is a major bioactive phospholipid of human plasma low-density lipoproteins (LDL). Platelet-activating factor (PAF) is a phospholipid similar to LPC and a potent intercellular mediator. Both PAF and LPC have been reported to act on mammalian cells through PAF receptor (PAFR). Previous data from our group suggested that T. cruzi produces a phospholipid with PAF activity. Here, we describe the structural and functional analysis of different species of LPC from T. cruzi, including a LPC with a fatty acid chain of 18 carbon atoms and one double bond (C18:1-LPC). We also show that C18:1-LPC is able to induce rabbit platelet aggregation, which is abrogated by a PAFR antagonist. In addition, a three-dimensional model of human PAFR was constructed. Contrary to other T. cruzi LPC molecules, C18:1-LPC is predicted to interact with the PAFR model in a fashion similar to PAF. Further studies are needed to validate the biosynthesis of T. cruzi C18:1-LPC as a potential drug target in Chagas disease.

ProQM-resample: improved model quality assessment for membrane proteins by limited conformational sampling

Summary: Model Quality Assessment Programs (MQAPs) are used to predict the quality of modeled protein structures. These usually use two approaches: methods using consensus of many alternative models and methods requiring only a single model to do its prediction. The consensus methods are useful to improve overall accuracy; however, they frequently fail to pick out the best possible model and cannot be used to generate and score new structures. Single-model methods, on the other hand, do not have these inherent shortcomings and can be used to both sample new structures and improve existing consensus methods. Here, we present ProQM-resample, a membrane protein-specific single-model MQAP, that couples side-chain resampling with MQAP rescoring by ProQM to improve model selection. The side-chain resampling is able to improve side-chain packing for 96% of all models, and improve model selection by 24% as measured by the sum of the Z-score for the first-ranked model (from 25.0 to 31.1), even better than the state-of-the-art consensus method Pcons. The improved model selection can be attributed to the improved side-chain quality, which enables the MQAP to rescue good backbone models with poor side-chain packing.

Availability and implementation:http://proqm.wallnerlab.org/download/.

Contact:bjornw@ifm.liu.se

Supplementary information:Supplementary data are available at Bioinformatics online.

Energetic and molecular water permeation mechanisms of the human red blood cell urea transporter B

Urea transporter B (UT-B) is a passive membrane channel that facilitates highly efficient permeation of urea. In red blood cells (RBC), while the major function of UT-B is to transport urea, it is assumed that this protein is able to conduct water. Here, we have revisited this last issue by studying RBCs and ghosts from human variants with defects of aquaporin 1 (AQP1) or UT-B. We found that UT-B's osmotic water unit permeability (pfunit) is similar to that of AQP1. The determination of diffusional permeability coefficient (Pd) allowed the calculation of the Pf/Pd ratio, which is consistent with a single-file water transport. Molecular dynamic simulations of water conduction through human UT-B confirmed the experimental finding. From these results, we propose an atomistic description of water-protein interactions involved in this permeation. Inside the UT-B pore, five water molecules were found to form a single-file and move rapidly along a channel by hydrogen bond exchange involving two critical threonines. We further show that the energy barrier for water located in the central region coincides with a water dipole reorientation, which can be related to the proton exclusion observed experimentally. In conclusion, our results indicate that UT-B should be considered as a new member of the water channel family.

Improved model quality assessment using ProQ2

Background

Employing methods to assess the quality of modeled protein structures is now standard practice in bioinformatics. In a broad sense, the techniques can be divided into methods relying on consensus prediction on the one hand, and single-model methods on the other. Consensus methods frequently perform very well when there is a clear consensus, but this is not always the case. In particular, they frequently fail in selecting the best possible model in the hard cases (lacking consensus) or in the easy cases where models are very similar. In contrast, single-model methods do not suffer from these drawbacks and could potentially be applied on any protein of interest to assess quality or as a scoring function for sampling-based refinement.

Results

Here, we present a new single-model method, ProQ2, based on ideas from its predecessor, ProQ. ProQ2 is a model quality assessment algorithm that uses support vector machines to predict local as well as global quality of protein models. Improved performance is obtained by combining previously used features with updated structural and predicted features. The most important contribution can be attributed to the use of profile weighting of the residue specific features and the use features averaged over the whole model even though the prediction is still local.

Conclusions

ProQ2 is significantly better than its predecessors at detecting high quality models, improving the sum of Z-scores for the selected first-ranked models by 20% and 32% compared to the second-best single-model method in CASP8 and CASP9, respectively. The absolute quality assessment of the models at both local and global level is also improved. The Pearson’s correlation between the correct and local predicted score is improved from 0.59 to 0.70 on CASP8 and from 0.62 to 0.68 on CASP9; for global score to the correct GDT_TS from 0.75 to 0.80 and from 0.77 to 0.80 again compared to the second-best single methods in CASP8 and CASP9, respectively. ProQ2 is available at http://proq2.wallnerlab.org.