Relation between sequence and structure in membrane proteins

RNA-protein interaction prediction using network-guided deep learning

Accurate computational determination of RNA-protein interactions remains challenging, particularly when encountering unknown RNAs and proteins. The limited number of RNAs and their flexibility constrained the effectiveness of the deep-learning models for RNA-protein interaction prediction. Here, we introduce ZHMolGraph, which integrates graph neural network and unsupervised large language models to predict RNA-protein interaction. We validate ZHMolGraph predictions on two benchmark datasets and outperform the current best methods. For the dataset of entirely unknown RNAs and proteins, ZHMolGraph shows an improvement in achieving high AUROC of 79.8% and AUPRC of 82.0%. This represents a substantial improvement of 7.1%-28.7% in AUROC and 4.6%-30.0% in AUPRC over other methods. We utilize ZHMolGraph to enhance the challenging SARS-CoV-2 RPI and unbound RNA-protein complex predictions. Such enhancements make ZHMolGraph a reliable option for genome-wide RNA-protein prediction. ZHMolGraph holds broad potential for modeling and designing RNA-protein complexes.

Improving the Annotations of JCVI-Syn3a Proteins

The JCVI-Syn3 organism is a minimal organism derived from Mycoplasma mycoides capri, which is capable of self-replication. While the ancestor has 863 genes, the synthetic progeny has only 473, with 434 of these coding for proteins. Despite initial efforts to understand all functions of the organism, a significant number of these protein-coding genes still have unknown functions, and subsequent studies have been only partially successful in elucidating their roles. In this study, we employ our innovative method PROST to identify homologs and better understand these previously unidentified genes. PROST employs protein language embeddings and enables the identification of remote homologs with as low as 16% sequence identity. PROST successfully finds functionally annotated homologs for 93% of the minimal genome with a high level of accuracy, both confirming previously identified functions, as well as proposing new functions for others. The results of our study can be accessed at https://bit.ly/prost-syn3a .

Living in trinity of extremes: Genomic and proteomic signatures of halophilic, thermophilic, and pH adaptation

Since nucleic acids and proteins of unicellular prokaryotes are directly exposed to extreme environmental conditions, it is possible to explore the genomic-proteomic compositional determinants of molecular mechanisms of adaptation developed by them in response to harsh environmental conditions. Using a wealth of currently available complete genomes/proteomes we were able to explore signatures of adaptation to three environmental factors, pH, salinity, and temperature, observing major trends in compositions of their nucleic acids and proteins. We derived predictors of thermostability, halophilic, and pH adaptations and complemented them by the principal components analysis. We observed a clear difference between thermophilic and salinity/pH adaptations, whereas latter invoke seemingly overlapping mechanisms. The genome-proteome compositional trade-off reveals an intricate balance between the work of base paring and base stacking in stabilization of coding DNA and r/tRNAs, and, at the same time, universal requirements for the stability and foldability of proteins regardless of the nucleotide biases. Nevertheless, we still found hidden fingerprints of ancient evolutionary connections between the nucleotide and amino acid compositions indicating their emergence, mutual evolution, and adjustment. The evolutionary perspective on the adaptation mechanisms is further studied here by means of the comparative analysis of genomic/proteomic traits of archaeal and bacterial species. The overall picture of genomic/proteomic signals of adaptation obtained here provides a foundation for future engineering and design of functional biomolecules resistant to harsh environments.

Unique Substrate Recognition and Sodium-Substrate Binding Stoichiometry in a Bacterial Serotonin Transporter, TuriSERT

All resolved high-resolution structures of the transporters in the neurotransmitter sodium symporter (NSS) family revealed that the NSS members share common structural and mechanistic features for substrate and ion binding and transport. However, a recently reported bacterial orthologue of the human serotonin transporter (hSERT), TuriSERT, possesses a structural characteristic specific for amino acid substrate binding but does transport a biogenic amine. The unique structural feature of TuriSERT requires a novel configuration for coordinating its substrate and ions. In the present study, we characterized TuriSERT expressed in Escherichia coli cells with a fluorescent substrate by biochemical, structural, and pharmacological approaches. Substrate transport by TuriSERT requires Na+ but not Cl-. Replacement of Asp262 by asparagine renders TuriSERT Cl--dependent. Substitutions of the corresponding Na1 residues did not alter Na+ dependence on substrate transport, whereas the mutation of a Na2 site residue led to a loss of transport activity, suggesting that Na+ binds only to the Na2 site in TuriSERT. In addition, substitutions of several residues essential for recognizing 5-hydroxytryptamine (5-HT) in hSERT had little effect on 5-HT displacement potency in transport assay for TuriSERT. In contrast, mutations of the residues that are proposed to coordinate with 5-HT in our docking model dramatically reduced 5-HT displacement. Furthermore, our results indicated that all tested antidepressants showed a weak inhibitory effect on TuriSERT. The present study demonstrated the existence of a unique substrate binding site and 1:1 stoichiometry of sodium-substrate binding in TuriSERT, a novel structural finding for the NSS transporters.

Exploring Flexibility and Folding Patterns Throughout Time in Voltage Sensors

The voltage-sensing domain (VSD) is a module capable of responding to changes in the membrane potential through conformational changes and facilitating electromechanical coupling to open a pore gate, activate proton permeation pathways, or promote enzymatic activity in some membrane-anchored phosphatases. To carry out these functions, this module acts cooperatively through conformational changes. The VSD is formed by four transmembrane segments (S1-S4) but the S4 segment is critical since it carries positively charged residues, mainly Arg or Lys, which require an aqueous environment for its proper function. The discovery of this module in voltage-gated ion channels (VGICs), proton channels (Hv1), and voltage sensor-containing phosphatases (VSPs) has expanded our understanding of the principle of modularity in the voltage-sensing mechanism of these proteins. Here, by sequence comparison and the evaluation of the relationship between sequence composition, intrinsic flexibility, and structural analysis in 14 selected representatives of these three major protein groups, we report five interesting differences in the folding patterns of the VSD both in prokaryotes and eukaryotes. Our main findings indicate that this module is highly conserved throughout the evolutionary scale, however: (1) segments S1 to S3 in eukaryotes are significantly more hydrophobic than those present in prokaryotes; (2) the S4 segment has retained its hydrophilic character; (3) in eukaryotes the extramembranous linkers are significantly larger and more flexible in comparison with those present in prokaryotes; (4) the sensors present in the kHv1 proton channel and the ciVSP phosphatase, both of eukaryotic origin, exhibit relationships of flexibility and folding patterns very close to the typical ones found in prokaryotic voltage sensors; and (5) archaeal channels KvAP and MVP have flexibility profiles which are clearly contrasting in the S3-S4 region, which could explain their divergent activation mechanisms. Finally, to elucidate the obscure origins of this module, we show further evidence for a possible connection between voltage sensors and TolQ proteins.

Molecular determinants of acrylamide neurotoxicity through covalent docking

Acrylamide (ACR) is formed during food processing by Maillard reaction between sugars and proteins at high temperatures. It is also used in many industries, from water waste treatment to manufacture of paper, fabrics, dyes and cosmetics. Unfortunately, cumulative exposure to acrylamide, either from diet or at the workplace, may result in neurotoxicity. Such adverse effects arise from covalent adducts formed between acrylamide and cysteine residues of several neuronal proteins via a Michael addition reaction. The molecular determinants of acrylamide reactivity and its impact on protein function are not completely understood. Here we have compiled a list of acrylamide protein targets reported so far in the literature in connection with neurotoxicity and performed a systematic covalent docking study. Our results indicate that acrylamide binding to cysteine is favored in the presence of nearby positively charged amino acids, such as lysines and arginines. For proteins with more than one reactive Cys, docking scores were able to discriminate between the primary ACR modification site and secondary sites modified only at high ACR concentrations. Therefore, docking scores emerge as a potential filter to predict Cys reactivity against acrylamide. Inspection of the ACR-protein complex structures provides insights into the putative functional consequences of ACR modification, especially for non-enzyme proteins. Based on our study, covalent docking is a promising computational tool to predict other potential protein targets mediating acrylamide neurotoxicity.

Computational Modeling of Antibody and T-Cell Receptor (CDR3 Loops)

Antibodies and T-cell receptors have been a subject of much interest due to their central role in the immune system and their potential applications in several biotechnological and medical applications from cancer therapy to vaccine development. A unique feature of these two lymphocyte receptors is their ability to bind a huge variety of different (pathogen) targets. This ability stems from six short loops in the binding domain that have hypervariable sequence due to genetic recombination mechanism. Particularly one of these loops, the third complementarity determining region (CDR3), has the highest sequence variability and a dominant role in binding the target. However, it has also been proven the most difficult to be modeled structurally, which is vitally important for downstream tasks such as binding prediction. This difficulty stems from its variability in sequence that both reduces the possibility of finding homologues and introduces unique structural features in the loop. We present here a general protocol for modeling such loops in antibodies and T-cell receptors. We also discuss the difficulties in loop modeling and the advantages and limitations of different modeling methods.

Identification of homologous GluN subunits variants accelerates GRIN variants stratification

The clinical spectrum of GRIN-related neurodevelopmental disorders (GRD) results from gene- and variant-dependent primary alterations of the NMDA receptor, disturbing glutamatergic neurotransmission. Despite GRIN gene variants' functional annotations being dually critical for stratification and precision medicine design, genetically diagnosed pathogenic GRIN variants currently outnumber their relative functional annotations. Based on high-resolution crystal 3D models and topological domains conservation between GluN1, GluN2A, and GluN2B subunits of the NMDAR, we have generated GluN1-GluN2A-GluN2B subunits structural superimposition model to find equivalent positions between GluN subunits. We have developed a GRIN structural algorithm that predicts functional changes in the equivalent structural positions in other GluN subunits. GRIN structural algorithm was computationally evaluated to the full GRIN missense variants repertoire, consisting of 4,525 variants. The analysis of this structure-based model revealed an absolute predictive power for GluN1, GluN2A, and GluN2B subunits, both in terms of pathogenicity-association (benign vs. pathogenic variants) and functional impact (loss-of-function, benign, gain-of-function). Further, we validated this computational algorithm experimentally, using an in silico library of GluN2B-equivalent GluN2A artificial variants, designed from pathogenic GluN2B variants. Thus, the implementation of the GRIN structural algorithm allows to computationally predict the pathogenicity and functional annotations of GRIN variants, resulting in the duplication of pathogenic GRIN variants assignment, reduction by 30% of GRIN variants with uncertain significance, and increase by 70% of functionally annotated GRIN variants. Finally, GRIN structural algorithm has been implemented into GRIN variants Database (http://lmc.uab.es/grindb), providing a computational tool that accelerates GRIN missense variants stratification, contributing to clinical therapeutic decisions for this neurodevelopmental disorder.

Molecular docking and biochemical validation of (-)-syringaresinol-4-O-β-D-apiofuranosyl-(1→2)-β-D-glucopyranoside binding to an allosteric site in monoamine transporters

Albizia julibrissin Durazz is one of the most common herbs used for depression and anxiety treatment, but its mechanism of action as an antidepressant or anxiolytic drug have not been fully understood. We previously isolated and identified one lignan glycoside compound from Albizia Julibrissin Durazz, (-)-syringaresinol-4-O-β-D-apiofuranosyl-(1→2)-β-D-glucopyranoside (SAG), that inhibited all three monoamine transporters with a mechanism of action different from that of the conventional antidepressants. In this study, we generated homology models for human dopamine transporter and human norepinephrine transporter, based on the X-ray structure of Drosophila dopamine transporter, and conducted the molecular docking of SAG to all three human monoamine transporters. Our computational results indicated that SAG binds to an allosteric site (S2) that has been demonstrated to be formed by an aromatic pocket positioned in the scaffold domain in the extracellular vestibule connected to the central site (S1) in these monoamine transporters. In addition, we demonstrated that SAG stabilizes a conformation of serotonin transporter with both the extracellular and cytoplasmic pathways closed. Furthermore, we performed mutagenesis of the residues in both the allosteric and orthosteric sites to biochemically validate SAG binding in all three monoamine transporters. Our results are consistent with the molecular docking calculation and support the association of SAG with the allosteric site. We expect that this herbal molecule could become a lead compound for the development of new therapeutic agents with a novel mechanism of action.

Discovery and characterization of a novel sub-group of UbiA-type terpene cyclases with a distinct motif I

A novel sub-group of UbiA-type TCs with a distinct motif I were identified from fungi to generate structurally diverse sesquiterpenes.

Nickel import and export in the human pathogen Helicobacter pylori, perspectives from molecular modelling

The uptake of essential metal ions and the ability to extrude them when their excess causes toxicity are crucial processes for all living beings. Nickel is a virulence factor for several human pathogens and in particular for the human gastric pathogen Helicobacter pylori because of its crucial role in the catalytic activity of two Ni-dependent enzymes, urease and hydrogenase. H. pylori requires efficient uptake mechanisms to import Ni(II) because of its scarcity in the human body, but the molecular details of Ni(II) homeostasis are not fully known. Here we offer a structural framework for the machinery of Ni(II) import/export in H. pylori, obtained through comparative modelling and macromolecular docking. The model structures reported in this perspective are initial steps towards the understanding of these processes at the molecular level and in the direction to exploit them to eradicate infections caused by this family of pathogens. The differences between the structural models obtained by using both the recently released neural network-based approach implemented in AlphaFold2 and a more classical user-driven modelling procedure are also discussed.

Photopharmacology of Ion Channels through the Light of the Computational Microscope

The optical control and investigation of neuronal activity can be achieved and carried out with photoswitchable ligands. Such compounds are designed in a modular fashion, combining a known ligand of the target protein and a photochromic group, as well as an additional electrophilic group for tethered ligands. Such a design strategy can be optimized by including structural data. In addition to experimental structures, computational methods (such as homology modeling, molecular docking, molecular dynamics and enhanced sampling techniques) can provide structural insights to guide photoswitch design and to understand the observed light-regulated effects. This review discusses the application of such structure-based computational methods to photoswitchable ligands targeting voltage- and ligand-gated ion channels. Structural mapping may help identify residues near the ligand binding pocket amenable for mutagenesis and covalent attachment. Modeling of the target protein in a complex with the photoswitchable ligand can shed light on the different activities of the two photoswitch isomers and the effect of site-directed mutations on photoswitch binding, as well as ion channel subtype selectivity. The examples presented here show how the integration of computational modeling with experimental data can greatly facilitate photoswitchable ligand design and optimization. Recent advances in structural biology, both experimental and computational, are expected to further strengthen this rational photopharmacology approach.

Structure-Aware Mycobacterium tuberculosis Functional Annotation Uncloaks Resistance, Metabolic, and Virulence Genes

Accurate and timely functional genome annotation is essential for translating basic pathogen research into clinically impactful advances. Here, through literature curation and structure-function inference, we systematically update the functional genome annotation of Mycobacterium tuberculosis virulent type strain H37Rv. First, we systematically curated annotations for 589 genes from 662 publications, including 282 gene products absent from leading databases. Second, we modeled 1,711 underannotated proteins and developed a semiautomated pipeline that captured shared function between 400 protein models and structural matches of known function on Protein Data Bank, including drug efflux proteins, metabolic enzymes, and virulence factors. In aggregate, these structure- and literature-derived annotations update 940/1,725 underannotated H37Rv genes and generate hundreds of functional hypotheses. Retrospectively applying the annotation to a recent whole-genome transposon mutant screen provided missing function for 48% (13/27) of underannotated genes altering antibiotic efficacy and 33% (23/69) required for persistence during mouse tuberculosis (TB) infection. Prospective application of the protein models enabled us to functionally interpret novel laboratory generated pyrazinamide (PZA)-resistant mutants of unknown function, which implicated the emerging coenzyme A depletion model of PZA action in the mutants’ PZA resistance. Our findings demonstrate the functional insight gained by integrating structural modeling and systematic literature curation, even for widely studied microorganisms. Functional annotations and protein structure models are available at https://tuberculosis.sdsu.edu/H37Rv in human- and machine-readable formats.

IMPORTANCEMycobacterium tuberculosis, the primary causative agent of tuberculosis, kills more humans than any other infectious bacterium. Yet 40% of its genome is functionally uncharacterized, leaving much about the genetic basis of its resistance to antibiotics, capacity to withstand host immunity, and basic metabolism yet undiscovered. Irregular literature curation for functional annotation contributes to this gap. We systematically curated functions from literature and structural similarity for over half of poorly characterized genes, expanding the functionally annotated Mycobacterium tuberculosis proteome. Applying this updated annotation to recent in vivo functional screens added functional information to dozens of clinically pertinent proteins described as having unknown function. Integrating the annotations with a prospective functional screen identified new mutants resistant to a first-line TB drug, supporting an emerging hypothesis for its mode of action. These improvements in functional interpretation of clinically informative studies underscore the translational value of this functional knowledge. Structure-derived annotations identify hundreds of high-confidence candidates for mechanisms of antibiotic resistance, virulence factors, and basic metabolism and other functions key in clinical and basic tuberculosis research. More broadly, they provide a systematic framework for improving prokaryotic reference annotations.

A Bioinformatics Approach for the Prediction of Immunogenic Properties and Structure of the SARS-COV-2 B.1.617.1 Variant Spike Protein

Background

B.1.617.1, a variant of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causing respiratory illness is responsible for the second wave of COVID-19 and associated with a high incidence of infectivity and mortality. To mitigate the B.1.617.1 variant of SARS-CoV-2, deciphering the protein structure and immunological responses by employing bioinformatics tools for data mining and analysis is pivotal.

Objectives

Here, an in silico approach was employed for deciphering the structure and immune function of the subunit of spike (S) protein of SARS-CoV-2 B.1.617.1 variant.

Methods

The partial amino acid sequence of SARS-CoV-2 B.1.617.1 variant S protein was analyzed, and its putative secondary and tertiary structure was predicted. Immunogenic analyses including B- and T-cell epitopes, interferon-gamma (IFN-γ) response, chemokine, and protective antigens for SARS-CoV 2 S proteins were predicted using appropriate tools.

Results

B.1.617.1 variant S protein sequence was found to be highly stable and amphipathic. ABCpred and CTLpred analyses led to the identification of two potential antigenic B cell and T cell epitopes with starting amino acid positions at 60 and 82 (for B cell epitopes) and 54 and 98 (for T cell epitopes) having prediction scores > 0.8. Further, RAMPAGE tool was used for determining the allowed and disallowed regions of the three-dimensional predicted structure of SARS-CoV-2 B.1.617.1 variant S protein.

Conclusion

Together, the in silico analysis revealed the predicted structure of partial S protein, immunogenic properties, and possible regions for S protein of SARS-CoV-2 and provides a valuable prelude for engineering the targeted vaccine or drug against B.1.617.1 variant of SARS-CoV-2.

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.

5-(Indol-2-yl)pyrazolo[3,4-b]pyridines as a New Family of TASK-3 Channel Blockers: A Pharmacophore-Based Regioselective Synthesis

TASK channels belong to the two-pore-domain potassium (K_2P) channels subfamily. These channels modulate cellular excitability, input resistance, and response to synaptic stimulation. TASK-channel inhibition led to membrane depolarization. TASK-3 is expressed in different cancer cell types and neurons. Thus, the discovery of novel TASK-3 inhibitors makes these bioactive compounds very appealing to explore new cancer and neurological therapies. TASK-3 channel blockers are very limited to date, and only a few heterofused compounds have been reported in the literature. In this article, we combined a pharmacophore hypothesis with molecular docking to address for the first time the rational design, synthesis, and evaluation of 5-(indol-2-yl)pyrazolo[3,4-b]pyridines as a novel family of human TASK-3 channel blockers. Representative compounds of the synthesized library were assessed against TASK-3 using Fluorometric imaging plate reader-Membrane Potential assay (FMP). Inhibitory properties were validated using two-electrode voltage-clamp (TEVC) methods. We identified one active hit compound (MM-3b) with our systematic pipeline, exhibiting an IC₅₀ ≈ 30 μM. Molecular docking models suggest that compound MM-3b binds to TASK-3 at the bottom of the selectivity filter in the central cavity, similar to other described TASK-3 blockers such as A1899 and PK-THPP. Our in silico and experimental studies provide a new tool to predict and design novel TASK-3 channel blockers.

Bitter Taste and Olfactory Receptors: Beyond Chemical Sensing in the Tongue and the Nose

Abstract

The Up-and-Coming-Scientist section of the current issue of the Journal of Membrane Biology features the invited essay by Dr. Mercedes Alfonso-Prieto, Assistant Professor at the Forschungszentrum Jülich (FZJ), Germany, and the Heinrich-Heine University Düsseldorf, Vogt Institute for Brain Research. Dr. Alfonso-Prieto completed her doctoral degree in chemistry at the Barcelona Science Park, Spain, in 2009, pursued post-doctoral research in computational molecular sciences at Temple University, USA, and then, as a Marie Curie post-doctoral fellow at the University of Barcelona, worked on computations of enzyme reactions and modeling of photoswitchable ligands targeting neuronal receptors. In 2016, she joined the Institute for Advanced Science and the Institute for Computational Biomedicine at the FZJ, where she pursues research on modeling and simulation of chemical senses. The invited essay by Dr. Alfonso-Prieto discusses state-of-the-art modeling of molecular receptors involved in chemical sensing – the senses of taste and smell. These receptors, and computational methods to study them, are the focus of Dr. Alfonso-Prieto’s research. Recently, Dr. Alfonso-Prieto and colleagues have presented a new methodology to predict ligand binding poses for GPCRs, and extensive computations that deciphered the ligand selectivity determinants of bitter taste receptors. These developments inform our current understanding of how taste occurs at the molecular level.

Graphic Abstract

Structural Phylogeny of Different Allergens May Reveal Common Epitopic Footprint

BACKGROUND

The incidence of allergy is increasing at an alarming rate for the last few decades.

OBJECTIVE

Our present study is focused on finding out the structurally homologous motifs present in different proteinaceous allergens Methods: Significant number of protein sequences and their corresponding structures of various pollen, fungal, bacterial, and food allergens were retrieved, and the sequence and structural identity were analyzed.

RESULTS

Intra- and inter-sequence and their structural analysis of the proteinaceous allergens, resulted in no significant relationships among them. A few, but not negligible number of high structural similarities were observed within different groups of allergens from fungus, angiosperms, and animals (Aves and Mammalia).

CONCLUSION

Our in silico study on thirty-six different allergens showed a significant level of structural similarities among themselves, regardless of their sequences.

Structural and molecular basis of cross-seeding barriers in amyloids

Neurodegenerative disorders are frequently associated with β-sheet-rich amyloid deposits. Amyloid-forming proteins can aggregate under different structural conformations known as strains, which can exhibit a prion-like behavior and distinct pathophenotypes. Precise molecular determinants defining strain specificity and cross-strain interactions (cross-seeding) are currently unknown. The HET-s prion protein from the fungus Podospora anserina represents a model system to study the fundamental properties of prion amyloids. Here, we report the amyloid prion structure of HELLF, a distant homolog of the model prion HET-s. We find that these two amyloids, sharing only 17% sequence identity, have nearly identical β-solenoid folds but lack cross-seeding ability in vivo, indicating that prion specificity can differ in extremely similar amyloid folds. We engineer the HELLF sequence to explore the limits of the sequence-to-fold conservation and to pinpoint determinants of cross-seeding and prion specificity. We find that amyloid fold conservation occurs even at an exceedingly low level of identity to HET-s (5%). Next, we derive a HELLF-based sequence, termed HEC, able to breach the cross-seeding barrier in vivo between HELLF and HET-s, unveiling determinants controlling cross-seeding at residue level. These findings show that virtually identical amyloid backbone structures might not be sufficient for cross-seeding and that critical side-chain positions could determine the seeding specificity of an amyloid fold. Our work redefines the conceptual boundaries of prion strain and sheds light on key molecular features concerning an important class of pathogenic agents.

TMSNP: a web server to predict pathogenesis of missense mutations in the transmembrane region of membrane proteins

The massive amount of data generated from genome sequencing brings tons of newly identified mutations, whose pathogenic/non-pathogenic effects need to be evaluated. This has given rise to several mutation predictor tools that, in general, do not consider the specificities of the various protein groups. We aimed to develop a predictor tool dedicated to membrane proteins, under the premise that their specific structural features and environment would give different responses to mutations compared to globular proteins. For this purpose, we created TMSNP, a database that currently contains information from 2624 pathogenic and 196 705 non-pathogenic reported mutations located in the transmembrane region of membrane proteins. By computing various conservation parameters on these mutations in combination with annotations, we trained a machine-learning model able to classify mutations as pathogenic or not. TMSNP (freely available at http://lmc.uab.es/tmsnp/) improves considerably the prediction power of commonly used mutation predictors trained with globular proteins.

Characterization of two sugar transporters responsible for efficient xylose uptake in an oleaginous yeast Candida tropicalis SY005

Microbial conversion of lignocellulosic feedstock to the target bioproduct requires efficient assimilation of its constituent sugars, a large part of which comprises of glucose and xylose. This study aims to identify and characterize sugar transporters capable of xylose uptake in an oleaginous strain of the industrially relevant yeast Candida tropicalis. In silico database mining resulted in two sugar transporter proteins- CtStp1 and CtStp2, containing conserved amino acid residues and motifs that have been previously reported to be involved in xylose transport in other organisms. Several softwares predicted the likelihood of 10-12 transmembrane (TM) helices to be present in both the Stps, while molecular modelling showed 12 TM helices that were organized into a typical structure found in the major facilitator superfamily of transporters. Docking with different sugars also predicted favorable interactions. Heterologous expression in a Saccharomyces cerevisiae strain harboring functional xylose metabolic genes validated the broad substrate specificity of the two Stps. Each transporter supported prominent growth of recombinant S. cerevisiae strains on six sugars including xylose at various concentrations. Expression of CtSTP1 and CtSTP2 along with the xylose metabolic genes in yeast transformants grown in presence of xylose was confirmed by transcript detection. Growth curve and sugar consumption profiles revealed uptake of both glucose and xylose simultaneously by the recombinant yeast strains, though CtStp1 showed relatively less effect of glucose repression in mixed sugars and was a better transporter of xylose than CtStp2. Such glucose-xylose utilizing efficient transporters can be effective tools for developing co-fermenting yeasts through genetic engineering in future, with noteworthy applications in renewable biomass utilization.

In Silico Investigation of Bitter Hop-Derived Compounds and Their Cognate Bitter Taste Receptors

The typical bitter taste of beer is caused by adding hops (Humulus lupulus L.) during the wort boiling process. The bitter taste of hop-derived compounds was found to be mediated by three bitter taste receptors: TAS2R1, TAS2R14, and TAS2R40. In this work, structural bioinformatics analyses were used to characterize the binding modes of trans-isocohumulone, trans-isohumulone, trans-isoadhumulone, cis-isocohumulone, cis-isohumulone, cis-isoadhumulone, cohumulone, humulone, adhumulone, and 8-prenylnaringenin into the orthosteric binding site of their cognate receptors. A conserved asparagine in transmembrane 3 was found to be essential for the recognition of hop-derived compounds, whereas the surrounding residues in the binding site of the three receptors encode the ligand specificity. Hop-derived compounds are renowned bioactive molecules and are considered as potential hit molecules for drug discovery to treat metabolic diseases. A chemoinformatics analysis revealed that hop-derived compounds cluster in a different region of the chemical space compared to known bitter food-derived compounds, pinpointing hop-derived compounds as a very peculiar class of bitter compounds.

The Molecular Basis of JAZ-MYC Coupling, a Protein-Protein Interface Essential for Plant Response to Stressors

The jasmonic acid (JA) signaling pathway is one of the primary mechanisms that allow plants to respond to a variety of biotic and abiotic stressors. Within this pathway, the JAZ repressor proteins and the basic helix-loop-helix (bHLH) transcription factor MYC3 play a critical role. JA is a volatile organic compound with an essential role in plant immunity. The increase in the concentration of JA leads to the decoupling of the JAZ repressor proteins and the bHLH transcription factor MYC3 causing the induction of genes of interest. The primary goal of this study was to identify the molecular basis of JAZ-MYC coupling. For this purpose, we modeled and validated 12 JAZ-MYC3 3D in silico structures and developed a molecular dynamics/machine learning pipeline to obtain two outcomes. First, we calculated the average free binding energy of JAZ-MYC3 complexes, which was predicted to be -10.94 +/-2.67 kJ/mol. Second, we predicted which ones should be the interface residues that make the predominant contribution to the free energy of binding (molecular hotspots). The predicted protein hotspots matched a conserved linear motif SL••FL•••R, which may have a crucial role during MYC3 recognition of JAZ proteins. As a proof of concept, we tested, both in silico and in vitro, the importance of this motif on PEAPOD (PPD) proteins, which also belong to the TIFY protein family, like the JAZ proteins, but cannot bind to MYC3. By mutating these proteins to match the SL••FL•••R motif, we could force PPDs to bind the MYC3 transcription factor. Taken together, modeling protein-protein interactions and using machine learning will help to find essential motifs and molecular mechanisms in the JA pathway.

Ligand Pose Predictions for Human G Protein-Coupled Receptors: Insights from the Amber-Based Hybrid Molecular Mechanics/Coarse-Grained Approach

Human G protein-coupled receptors (hGPCRs) are the most frequent targets of Food and Drug Administration (FDA)-approved drugs. Structural bioinformatics, along with molecular simulation, can support structure-based drug design targeting hGPCRs. In this context, several years ago, we developed a hybrid molecular mechanics (MM)/coarse-grained (CG) approach to predict ligand poses in low-resolution hGPCR models. The approach was based on the GROMOS96 43A1 and PRODRG united-atom force fields for the MM part. Here, we present a new MM/CG implementation using, instead, the Amber 14SB and GAFF all-atom potentials for proteins and ligands, respectively. The new implementation outperforms the previous one, as shown by a variety of applications on models of hGPCR/ligand complexes at different resolutions, and it is also more user-friendly. Thus, it emerges as a useful tool to predict poses in low-resolution models and provides insights into ligand binding similarly to all-atom molecular dynamics, albeit at a lower computational cost.

EncoMPASS: an online database for analyzing structure and symmetry in membrane proteins

The EncoMPASS online database (http://encompass.ninds.nih.gov) collects, organizes, and presents information about membrane proteins of known structure, emphasizing their structural similarities as well as their quaternary and internal symmetries. Unlike, e.g. SCOP, the EncoMPASS database does not aim for a strict classification of membrane proteins, but instead is organized as a protein chain-centric network of sequence and structural homologues. The online server for the EncoMPASS database provides tools for comparing the structural features of its entries, making it a useful resource for homology modeling and active site identification studies. The database can also be used for inferring functionality, which for membrane proteins often involves symmetry-related mechanisms. To this end, the online database also provides a comprehensive description of both the quaternary and internal symmetries in known membrane protein structures, with a particular focus on their orientation relative to the membrane.

Dynamical model of the CLC-2 ion channel reveals conformational changes associated with selectivity-filter gating

This work reports a dynamical Markov state model of CLC-2 "fast" (pore) gating, based on 600 microseconds of molecular dynamics (MD) simulation. In the starting conformation of our CLC-2 model, both outer and inner channel gates are closed. The first conformational change in our dataset involves rotation of the inner-gate backbone along residues S168-G169-I170. This change is strikingly similar to that observed in the cryo-EM structure of the bovine CLC-K channel, though the volume of the intracellular (inner) region of the ion conduction pathway is further expanded in our model. From this state (inner gate open and outer gate closed), two additional states are observed, each involving a unique rotameric flip of the outer-gate residue GLUex. Both additional states involve conformational changes that orient GLUex away from the extracellular (outer) region of the ion conduction pathway. In the first additional state, the rotameric flip of GLUex results in an open, or near-open, channel pore. The equilibrium population of this state is low (∼1%), consistent with the low open probability of CLC-2 observed experimentally in the absence of a membrane potential stimulus (0 mV). In the second additional state, GLUex rotates to occlude the channel pore. This state, which has a low equilibrium population (∼1%), is only accessible when GLUex is protonated. Together, these pathways model the opening of both an inner and outer gate within the CLC-2 selectivity filter, as a function of GLUex protonation. Collectively, our findings are consistent with published experimental analyses of CLC-2 gating and provide a high-resolution structural model to guide future investigations.

Intestinal serotonin and fluoxetine exposure modulate bacterial colonization in the gut

The gut microbiota regulates levels of serotonin (5-hydroxytryptamine (5-HT)) in the intestinal epithelium and lumen^1-5. However, whether 5-HT plays a functional role in bacteria from the gut microbiota remains unknown. We demonstrate that elevating levels of intestinal lumenal 5-HT by oral supplementation or genetic deficiency in the host 5-HT transporter (SERT) increases the relative abundance of spore-forming members of the gut microbiota, which were previously reported to promote host 5-HT biosynthesis. Within this microbial community, we identify Turicibacter sanguinis as a gut bacterium that expresses a neurotransmitter sodium symporter-related protein with sequence and structural homology to mammalian SERT. T. sanguinis imports 5-HT through a mechanism that is inhibited by the selective 5-HT reuptake inhibitor fluoxetine. 5-HT reduces the expression of sporulation factors and membrane transporters in T. sanguinis, which is reversed by fluoxetine exposure. Treating T. sanguinis with 5-HT or fluoxetine modulates its competitive colonization in the gastrointestinal tract of antibiotic-treated mice. In addition, fluoxetine reduces the membership of T. sanguinis in the gut microbiota of conventionally colonized mice. Host association with T. sanguinis alters intestinal expression of multiple gene pathways, including those important for lipid and steroid metabolism, with corresponding reductions in host systemic triglyceride levels and inguinal adipocyte size. Together, these findings support the notion that select bacteria indigenous to the gut microbiota signal bidirectionally with the host serotonergic system to promote their fitness in the intestine.

A large intragenic deletion in the CLCN1 gene causes Hereditary Myotonia in pigs

Mutations in the CLCN1 gene are the primary cause of non-dystrophic Hereditary Myotonia in several animal species. However, there are no reports of Hereditary Myotonia in pigs to date. Therefore, the objective of the present study was to characterize the clinical and molecular findings of Hereditary Myotonia in an inbred pedigree. The clinical, electromyographic, histopathological, and molecular findings were evaluated. Clinically affected pigs presented non-dystrophic recessive Hereditary Myotonia. Nucleotide sequence analysis of the entire coding region of the CLCN1 gene revealed the absence of the exons 15 and 16 in myotonic animals. Analysis of the genomic region flanking the deletion unveiled a large intragenic deletion of 4,165 nucleotides. Interestingly, non-related, non-myotonic pigs expressed transcriptional levels of an alternate transcript (i.e., X2) that was identical to the deleted X1 transcript of myotonic pigs. All myotonic pigs and their progenitors were homozygous recessive and heterozygous, respectively, for the 4,165-nucleotide deletion. This is the first study reporting Hereditary Myotonia in pigs and characterizing its clinical and molecular findings. Moreover, to the best of our knowledge, Hereditary Myotonia has never been associated with a genomic deletion in the CLCN1 gene in any other species.

Inter-residue interactions in alpha-helical transmembrane proteins

MOTIVATION

The number of available membrane protein structures has markedly increased in the last years and, in parallel, the reliability of the methods to detect transmembrane (TM) segments. In the present report, we characterized inter-residue interactions in α-helical membrane proteins using a dataset of 3462 TM helices from 430 proteins. This is by far the largest analysis published to date.

RESULTS

Our analysis of residue-residue interactions in TM segments of membrane proteins shows that almost all interactions involve aliphatic residues and Phe. There is lack of polar-polar, polar-charged and charged-charged interactions except for those between Thr or Ser sidechains and the backbone carbonyl of aliphatic and Phe residues. The results are discussed in the context of the preferences of amino acids to be in the protein core or exposed to the lipid bilayer and to occupy specific positions along the TM segment. Comparison to datasets of β-barrel membrane proteins and of α-helical globular proteins unveils the specific patterns of interactions and residue composition characteristic of α-helical membrane proteins that are the clue to understanding their structure.

AVAILABILITY AND IMPLEMENTATION

Results data and datasets used are available at http://lmc.uab.cat/TMalphaDB/interactions.php.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Dual binding mode of "bitter sugars" to their human bitter taste receptor target

The 25 human bitter taste receptors (hTAS2Rs) are responsible for detecting bitter molecules present in food, and they also play several physiological and pathological roles in extraoral compartments. Therefore, understanding their ligand specificity is important both for food research and for pharmacological applications. Here we provide a molecular insight into the exquisite molecular recognition of bitter β-glycopyranosides by one of the members of this receptor subclass, hTAS2R16. Most of its agonists have in common the presence of a β-glycopyranose unit along with an extremely structurally diverse aglycon moiety. This poses the question of how hTAS2R16 can recognize such a large number of "bitter sugars". By means of hybrid molecular mechanics/coarse grained molecular dynamics simulations, here we show that the three hTAS2R16 agonists salicin, arbutin and phenyl-β-D-glucopyranoside interact with the receptor through a previously unrecognized dual binding mode. Such mechanism may offer a seamless way to fit different aglycons inside the binding cavity, while maintaining the sugar bound, similar to the strategy used by several carbohydrate-binding lectins. Our prediction is validated a posteriori by comparison with mutagenesis data and also rationalizes a wealth of structure-activity relationship data. Therefore, our findings not only provide a deeper molecular characterization of the binding determinants for the three ligands studied here, but also give insights applicable to other hTAS2R16 agonists. Together with our results for other hTAS2Rs, this study paves the way to improve our overall understanding of the structural determinants of ligand specificity in bitter taste receptors.

Inhibition of the Alanine-Serine-Cysteine-1 Transporter by BMS-466442

Experiment and modeling were combined to understand inhibition of the alanine-serine-cysteine-1 (asc1) transporter. The structure-activity relationship (SAR) was explored with synthesis of analogues of BMS-466442. Direct target interaction and binding site location between TM helices 6 and 10 were confirmed via site directed mutagenesis. Computational modeling suggested the inhibitor binds via competitive occupation of the orthosteric site while also blocking the movement of TM helices that are required for transport.

The counterion-retinylidene Schiff base interaction of an invertebrate rhodopsin rearranges upon light activation

Animals sense light using photosensitive proteins-rhodopsins-containing a chromophore-retinal-that intrinsically absorbs in the ultraviolet. Visible light-sensitivity depends primarily on protonation of the retinylidene Schiff base (SB), which requires a negatively-charged amino acid residue-counterion-for stabilization. Little is known about how the most common counterion among varied rhodopsins, Glu181, functions. Here, we demonstrate that in a spider visual rhodopsin, orthologue of mammal melanopsins relevant to circadian rhythms, the Glu181 counterion functions likely by forming a hydrogen-bonding network, where Ser186 is a key mediator of the Glu181-SB interaction. We also suggest that upon light activation, the Glu181-SB interaction rearranges while Ser186 changes its contribution. This is in contrast to how the counterion of vertebrate visual rhodopsins, Glu113, functions, which forms a salt bridge with the SB. Our results shed light on the molecular mechanisms of visible light-sensitivity relevant to invertebrate vision and vertebrate non-visual photoreception.

Understanding Ligand Binding to G-Protein Coupled Receptors Using Multiscale Simulations

Human G-protein coupled receptors (GPCRs) convey a wide variety of extracellular signals inside the cell and they are one of the main targets for pharmaceutical intervention. Rational drug design requires structural information on these receptors; however, the number of experimental structures is scarce. This gap can be filled by computational models, based on homology modeling and docking techniques. Nonetheless, the low sequence identity across GPCRs and the chemical diversity of their ligands may limit the quality of these models and hence refinement using molecular dynamics simulations is recommended. This is the case for olfactory and bitter taste receptors, which constitute the first and third largest GPCR groups and show sequence identities with the available GPCR templates below 20%. We have developed a molecular dynamics approach, based on the combination of molecular mechanics and coarse grained (MM/CG), tailored to study ligand binding in GPCRs. This approach has been applied so far to bitter taste receptor complexes, showing significant predictive power. The protein/ligand interactions observed in the simulations were consistent with extensive mutagenesis and functional data. Moreover, the simulations predicted several binding residues not previously tested, which were subsequently verified by carrying out additional experiments. Comparison of the simulations of two bitter taste receptors with different ligand selectivity also provided some insights into the binding determinants of bitter taste receptors. Although the MM/CG approach has been applied so far to a limited number of GPCR/ligand complexes, the excellent agreement of the computational models with the mutagenesis and functional data supports the applicability of this method to other GPCRs for which experimental structures are missing. This is particularly important for the challenging case of GPCRs with low sequence identity with available templates, for which molecular docking shows limited predictive power.

Multiscale simulations on human Frizzled and Taste2 GPCRs

Recently, molecular dynamics simulations, from all atom and coarse grained to hybrid methods bridging the two scales, have provided exciting functional insights into class F (Frizzled and Taste2) GPCRs (about 40 members in humans). Findings include: (i) The activation of one member of the Frizzled receptors (FZD4) involves a bending of transmembrane helix TM7 far larger than that in class A GPCRs. (ii) The affinity of an anticancer drug targeting another member (Smoothened receptor) decreases in a specific drug-resistant variant, because the mutation ultimately disrupts the binding cavity and affects TM6. (iii) A novel two-state recognition mechanism explains the very large agonist diversity for at least one member of the Taste2 GPCRs, hTAS2R46.

Repository of Enriched Structures of Proteins Involved in the Red Blood Cell Environment (RESPIRE)

The Red Blood Cell (RBC) is a metabolically-driven cell vital for processes such a gas transport and homeostasis. RBC possesses at its surface exposing antigens proteins that are critical in blood transfusion. Due to their importance, numerous studies address the cell function as a whole but more and more details of RBC structure and protein content are now studied using massive state-of-the art characterisation techniques. Yet, the resulting information is frequently scattered in many scientific articles, in many databases and specialized web servers. To provide a more compendious view of erythrocytes and of their protein content, we developed a dedicated database called RESPIRE that aims at gathering a comprehensive and coherent ensemble of information and data about proteins in RBC. This cell-driven database lists proteins found in erythrocytes. For a given protein entry, initial data are processed from external portals and enriched by using state-of-the-art bioinformatics methods. As structural information is extremely useful to understand protein function and predict the impact of mutations, a strong effort has been put on the prediction of protein structures with a special treatment for membrane proteins. Browsing the database is available through text search for reference gene names or protein identifiers, through pre-defined queries or via hyperlinks. The RESPIRE database provides valuable information and unique annotations that should be useful to a wide audience of biologists, clinicians and structural biologists.

Database URL:http://www.dsimb.inserm.fr/respire

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.

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.

Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

Sulfate Transporters in Dissimilatory Sulfate Reducing Microorganisms: A Comparative Genomics Analysis

The first step in the sulfate reduction pathway is the transport of sulfate across the cell membrane. This uptake has a major effect on sulfate reduction rates. Much of the information available on sulfate transport was obtained by studies on assimilatory sulfate reduction, where sulfate transporters were identified among several types of protein families. Despite our growing knowledge on the physiology of dissimilatory sulfate-reducing microorganisms (SRM) there are no studies identifying the proteins involved in sulfate uptake in members of this ecologically important group of anaerobes. We surveyed the complete genomes of 44 sulfate-reducing bacteria and archaea across six phyla and identified putative sulfate transporter encoding genes from four out of the five surveyed protein families based on homology. We did not find evidence that ABC-type transporters (SulT) are involved in the uptake of sulfate in SRM. We speculate that members of the CysP sulfate transporters could play a key role in the uptake of sulfate in thermophilic SRM. Putative CysZ-type sulfate transporters were present in all genomes examined suggesting that this overlooked group of sulfate transporters might play a role in sulfate transport in dissimilatory sulfate reducers alongside SulP. Our in silico analysis highlights several targets for further molecular studies in order to understand this key step in the metabolism of SRMs.

Predicting ligand binding poses for low-resolution membrane protein models: Perspectives from multiscale simulations

Membrane receptors constitute major targets for pharmaceutical intervention. Drug design efforts rely on the identification of ligand binding poses. However, the limited experimental structural information available may make this extremely challenging, especially when only low-resolution homology models are accessible. In these cases, the predictions may be improved by molecular dynamics simulation approaches. Here we review recent developments of multiscale, hybrid molecular mechanics/coarse-grained (MM/CG) methods applied to membrane proteins. In particular, we focus on our in-house MM/CG approach. It is especially tailored for G-protein coupled receptors, the largest membrane receptor family in humans. We show that our MM/CG approach is able to capture the atomistic details of the receptor/ligand binding interactions, while keeping the computational cost low by representing the protein frame and the membrane environment in a highly simplified manner. We close this review by discussing ongoing improvements and challenges of the current implementation of our MM/CG code.

Comparative analysis of membrane protein structure databases

BACKGROUND

Membrane proteins play important roles in cell survival and cell communication, as they function as transporters, receptors, anchors and enzymes. They are also potential targets for drugs that block receptors or inhibit enzymes related to diseases. Although the number of known structures of membrane proteins is still small relative to the size of the proteome as a whole, many new membrane protein structures have been determined recently.

SCOPE OF THE ARTICLE

We compared and analyzed the widely used membrane protein databases, mpstruc, Orientations of Proteins in Membranes (OPM), and PDBTM, as well as the extended dataset of mpstruc based on sequence similarity, the PDB structures whose classification field indicates that they are "membrane proteins" and the proteins with Structural Classification of Proteins (SCOP) class-f domains. We evaluated the relationships between these databases or datasets based on the overlap in their contents and the degree of consistency in the structural, topological, and functional classifications and in the transmembrane domain assignment.

MAJOR CONCLUSIONS

The membrane databases differ from each other in their coverage, and in the criteria that they use for annotation and classification. To ensure the efficient use of these databases, it is important to understand their differences and similarities. The establishment of more detailed and consistent annotations for the sequence, structure, membrane association, and function of membrane proteins is still required.

GENERAL SIGNIFICANCE

Considering the recent growth of experimentally determined structures, a broad survey and cumulative analysis of the sum of knowledge as presented in the membrane protein structure databases can be helpful to elucidate structures and functions of membrane proteins. We also aim to provide a framework for future research and classification of membrane proteins.

In silico model of the human ClC-Kb chloride channel: pore mapping, biostructural pathology and drug screening

The human ClC-Kb channel plays a key role in exporting chloride ions from the cytosol and is known to be involved in Bartter syndrome type 3 when its permeation capacity is decreased. The ClC-Kb channel has been recently proposed as a potential therapeutic target to treat hypertension. In order to gain new insights into the sequence-structure-function relationships of this channel, to investigate possible impacts of amino-acid substitutions, and to design novel inhibitors, we first built a structural model of the human ClC-Kb channel using comparative modeling strategies. We combined in silico and in vitro techniques to analyze amino acids involved in the chloride ion pathway as well as to rationalize the possible role of several clinically observed mutations leading to the Bartter syndrome type 3. Virtual screening and drug repositioning computations were then carried out. We identified six novel molecules, including 2 approved drugs, diflusinal and loperamide, with Kd values in the low micromolar range, that block the human ClC-Kb channel and that could be used as starting point to design novel chemical probes for this potential therapeutic target.

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.

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.

Accessing gap-junction channel structure-function relationships through molecular modeling and simulations

Background

Gap junction channels (GJCs) are massive protein channels connecting the cytoplasm of adjacent cells. These channels allow intercellular transfer of molecules up to ~1 kDa, including water, ions and other metabolites. Unveiling structure-function relationships coded into the molecular architecture of these channels is necessary to gain insight on their vast biological function including electrical synapse, inflammation, development and tissular homeostasis. From early works, computational methods have been critical to analyze and interpret experimental observations. Upon the availability of crystallographic structures, molecular modeling and simulations have become a valuable tool to assess structure-function relationships in GJCs. Modeling different connexin isoforms, simulating the transport process, and exploring molecular variants, have provided new hypotheses and out-of-the-box approaches to the study of these important channels.

Methods

Here, we review foundational structural studies and recent developments on GJCs using molecular modeling and simulation techniques, highlighting the methods and the cross-talk with experimental evidence.

Results and discussion

By comparing results obtained by molecular modeling and simulations techniques with structural and functional information obtained from both recent literature and structural databases, we provide a critical assesment of structure-function relationships that can be obtained from the junction between theoretical and experimental evidence.