Detecting pore-lining regions in transmembrane protein sequences

Structural modelling and preventive strategy targeting of WSSV hub proteins to combat viral infection in shrimp Penaeus monodon

White spot syndrome virus (WSSV) presents a considerable peril to the aquaculture sector, leading to notable financial consequences on a global scale. Previous studies have identified hub proteins, including WSSV051 and WSSV517, as essential binding elements in the protein interaction network of WSSV. This work further investigates the functional structures and potential applications of WSSV hub complexes in managing WSSV infection. Using computational methodologies, we have successfully generated comprehensive three-dimensional (3D) representations of hub proteins along with their three mutual binding counterparts, elucidating crucial interaction locations. The results of our study indicate that the WSSV051 hub protein demonstrates higher binding energy than WSSV517. Moreover, a unique motif, denoted as "S-S-x(5)-S-x(2)-P," was discovered among the binding proteins. This pattern perhaps contributes to the detection of partners by the hub proteins of WSSV. An antiviral strategy targeting WSSV hub proteins was demonstrated through the oral administration of dual hub double-stranded RNAs to the black tiger shrimp, Penaeus monodon, followed by a challenge assay. The findings demonstrate a decrease in shrimp mortality and a cessation of WSSV multiplication. In conclusion, our research unveils the structural features and dynamic interactions of hub complexes, shedding light on their significance in the WSSV protein network. This highlights the potential of hub protein-based interventions to mitigate the impact of WSSV infection in aquaculture.

The Plasmodiophora brassicae Golgi-localized UPF0016 protein PbGDT1 mediates calcium but not manganese transport in yeast and Nicotiana benthamiana

Manganese and calcium homeostasis and signalling, in eukaryotic organisms, are regulated through membrane located pumps, channels and exchangers, including the Mn2+/Ca2+ uncharacterized protein family 0016 (UPF0016). Here we show that Plasmodiophora brassicae PbGDT1 is a member of the UPF0016 and an ortholog of Saccharomyces cerevisiae Gdt1p (GCR Dependent Translation Factor 1) protein involved in manganese homeostasis as well as the calcium mediated stress response in yeast. PbGDT1 complemented the ScGdt1p and ScPMR1 (Ca2+ ATPase) double null mutant under elevated calcium stress but not under elevated manganese conditions. In both yeast and Nicotiana benthamiana, PbGDT1 localizes to the Golgi apparatus, with additional ER association in N. benthamiana. Expression of PbGDT1 in N. benthamiana, suppresses BAX-triggered cell death, further highlighting the importance of calcium homeostasis in maintaining cell physiology and integrity in a stress environment.

Halocin H4 is activated through cleavage by halolysin HlyR4

Halocins are antimicrobial peptides secreted by haloarchaea capable of inhibiting the growth of other haloarchaea or bacteria. Halocin H4 (HalH4) is secreted by the model halophilic archaeon Haloferax mediterranei ATCC 33500. Despite attempts to express halH4 heterologously in Escherichia coli and subsequent careful renaturation procedures commonly employed for haloarchaeal proteins, no active halocin was obtained. However, it was discovered that the antihaloarchaeal activity of this halocin could be activated through cleavage by halolysin R4 (HlyR4), a serine protease also secreted by Hfx. mediterranei ATCC 33500. Replacement of the cysteine at the number 115 amino acid with glycine and deletion of the internal trans-membrane region (15 aa) markedly abolished HalH4's antihaloarchaeal activity. Compared to the N-terminus, the C-terminal amino acid sequence was found to be more crucial for HalH4 to exert its antihaloarchaeal activity. Mass spectrometry analysis revealed that the biologically active antihaloarchaeal peptide produced after hydrolytic cleavage by HlyR4 was the C-terminus of HalH4, suggesting a potential mechanism of action involving pore formation within competitor species' cell membranes. Taken together, this study offers novel insights into the interplay between halocins and secreted proteases, as well as their contribution to antagonistic interaction within haloarchaea.

IMPORTANCE

The antihaloarchaeal function of halocin H4 (HalH4) can be activated by extracellular proteases from haloarchaea, as demonstrated in this study. Notably, we report the first instance of halocin activation through proteolytic cleavage, highlighting its significance in the field. The C-terminus of HalH4 (CTH4) has been identified as the antihaloarchaeal peptide present in hydrolysates generated by HlyR4. The CTH4 exhibited inhibitory activity against a range of haloarchaeal species (Haloarchaeobius spp., Haloarcula spp., Haloferax spp., Halorubellus spp., and Halorubrum spp.), as well as selected bacterial species (Aliifodinibius spp. and Salicola spp.), indicating its broad-spectrum inhibitory potential across domains. The encoding gene of halocin HalH4, halH4, from the model halophilic archaeon Haloferax mediterranei ATCC 33500 can be expressed in Escherichia coli without codon optimization.

Verticillium dahliae VdPBP1 Transcription Factor Is Required for Hyphal Growth, Virulence, and Microsclerotia Formation

Verticillium dahliae, a fungal pathogen that affects more than 200 plant species, including tomatoes, requires specific proteins for its early steps in plant infection. One such crucial protein, VdPBP1, exhibits high expression in the presence of tomato roots. Its 313-amino acid C-terminal section restores adhesion in nonadhesive Saccharomyces cerevisiae strains. To uncover its role, we employed a combination of bioinformatics, genetics, and morphological analyses. Our findings underscore the importance of VdPBP1 in fungal growth and pathogenesis. Bioinformatic analysis revealed that the VdPBP1 gene consists of four exons and three introns, encoding a 952-codon reading frame. The protein features a 9aaTAD domain, LsmAD, and PAB1 DNA-binding sites, as well as potential nuclear localization and transmembrane helix signals. Notably, the deletion of a 1.1 kb fragment at the gene's third end impedes microsclerotia formation and reduces pathogenicity. Mutants exhibit reduced growth and slower aerial mycelial development compared to the wild type. The VdPBP1 deletion strain does not induce disease symptoms in tomato plants. Furthermore, VdPBP1 deletion correlates with downregulated microsclerotia formation-related genes, and promoter analysis reveals regulatory elements, including sites for Rfx1, Mig1, and Ste12 proteins. Understanding the regulation and target genes of VdPBP1 holds promise for managing Verticillium wilt disease and related fungal pathogens.

Brewpitopes: a pipeline to refine B-cell epitope predictions during public health emergencies

The application of B-cell epitope identification to develop therapeutic antibodies and vaccine candidates is well established. However, the validation of epitopes is time-consuming and resource-intensive. To alleviate this, in recent years, multiple computational predictors have been developed in the immunoinformatics community. Brewpitopes is a pipeline that curates bioinformatic B-cell epitope predictions obtained by integrating different state-of-the-art tools. We used additional computational predictors to account for subcellular location, glycosylation status, and surface accessibility of the predicted epitopes. The implementation of these sets of rational filters optimizes in vivo antibody recognition properties of the candidate epitopes. To validate Brewpitopes, we performed a proteome-wide analysis of SARS-CoV-2 with a particular focus on S protein and its variants of concern. In the S protein, we obtained a fivefold enrichment in terms of predicted neutralization versus the epitopes identified by individual tools. We analyzed epitope landscape changes caused by mutations in the S protein of new viral variants that were linked to observed immune escape evidence in specific strains. In addition, we identified a set of epitopes with neutralizing potential in four SARS-CoV-2 proteins (R1AB, R1A, AP3A, and ORF9C). These epitopes and antigenic proteins are conserved targets for viral neutralization studies. In summary, Brewpitopes is a powerful pipeline that refines B-cell epitope bioinformatic predictions during public health emergencies in a high-throughput capacity to facilitate the optimization of experimental validation of therapeutic antibodies and candidate vaccines.

Dimeric Transmembrane Structure of the SARS-CoV-2 E Protein

The SARS-CoV-2 E protein is a transmembrane (TM) protein with its N-terminus exposed on the external surface of the virus. At debate is its oligomeric state, let alone its function. Here, the TM structure of the E protein is characterized by oriented sample and magic angle spinning solid-state NMR in lipid bilayers and refined by molecular dynamics simulations. This protein was previously found to be a pentamer, with a hydrophobic pore that appears to function as an ion channel. We identify only a front-to-front, symmetric helix-helix interface, leading to a dimeric structure that does not support channel activity. The two helices have a tilt angle of only 6°, resulting in an extended interface dominated by Leu and Val sidechains. While residues Val14-Thr35 are almost all buried in the hydrophobic region of the membrane, Asn15 lines a water-filled pocket that potentially serves as a drug-binding site. The E and other viral proteins may adopt different oligomeric states to help perform multiple functions.

Large-scale ORF screening based on LC-MS to discover novel lncRNA-encoded peptides responding to ionizing radiation and microgravity

In the human genome, 98% of genes can be transcribed into non-coding RNAs (ncRNAs), among which lncRNAs and their encoded peptides play important roles in regulating various aspects of cellular processes and may serve as crucial factors in modulating the biological effects induced by ionizing radiation and microgravity. Unfortunately, there are few reports in space radiation biology on lncRNA-encoded peptides below 10kD due to limitations in detection techniques. To fill this gap, we integrated a variety of methods based on genomics and peptidomics, and discovered 22 lncRNA-encoded small peptides that are sensitive to space radiation and microgravity, which have never been reported before. We concurrently validated the transmembrane helix, subcellular localization, and biological function of these small peptides using bioinformatics and molecular biology techniques. More importantly, we found that these small peptides function independently of the lncRNAs that encode them. Our findings have uncovered a previously unknown human proteome encoded by 'non-coding' genes in response to space conditions and elucidated their involvement in biological processes, providing valuable strategies for individual protection mechanisms for astronauts who carry out deep space exploration missions in space radiation environments.

Lagovirus Non-structural Protein p23: A Putative Viroporin That Interacts With Heat Shock Proteins and Uses a Disulfide Bond for Dimerization

The exact function(s) of the lagovirus non-structural protein p23 is unknown as robust cell culture systems for the Rabbit haemorrhagic disease virus (RHDV) and other lagoviruses have not been established. Instead, a range of in vitro and in silico models have been used to study p23, revealing that p23 oligomerizes, accumulates in the cytoplasm, and possesses a conserved C-terminal region with two amphipathic helices. Furthermore, the positional homologs of p23 in other caliciviruses have been shown to possess viroporin activity. Here, we report on the mechanistic details of p23 oligomerization. Site-directed mutagenesis revealed the importance of an N-terminal cysteine for dimerization. Furthermore, we identified cellular interactors of p23 using stable isotope labeling with amino acids in cell culture (SILAC)-based proteomics; heat shock proteins Hsp70 and 110 interact with p23 in transfected cells, suggesting that they ‘chaperone’ p23 proteins before their integration into cellular membranes. We investigated changes to the global transcriptome and proteome that occurred in infected rabbit liver tissue and observed changes to the misfolded protein response, calcium signaling, and the regulation of the endoplasmic reticulum (ER) network. Finally, flow cytometry studies indicate slightly elevated calcium concentrations in the cytoplasm of p23-transfected cells. Taken together, accumulating evidence suggests that p23 is a viroporin that might form calcium-conducting channels in the ER membranes.

Diversity and conservation of plant small secreted proteins associated with arbuscular mycorrhizal symbiosis

Arbuscular mycorrhizal symbiosis (AMS) is widespread mutualistic association between plants and fungi, which plays an essential role in nutrient exchange, enhancement in plant stress resistance, development of host, and ecosystem sustainability. Previous studies have shown that plant small secreted proteins (SSPs) are involved in beneficial symbiotic interactions. However, the role of SSPs in the evolution of AMS has not been well studied yet. In this study, we performed computational analysis of SSPs in 60 plant species and identified three AMS-specific ortholog groups containing SSPs only from at least 30% of the AMS species in this study and three AMS-preferential ortholog groups containing SSPs from both AMS and non-AMS species, with AMS species containing significantly more SSPs than non-AMS species. We found that independent lineages of monocot and eudicot plants contained genes in the AMS-specific ortholog groups and had significant expansion in the AMS-preferential ortholog groups. Also, two AMS-preferential ortholog groups showed convergent changes, between monocot and eudicot species, in gene expression in response to arbuscular mycorrhizal fungus Rhizophagus irregularis. Furthermore, conserved cis-elements were identified in the promoter regions of the genes showing convergent gene expression. We found that the SSPs, and their closely related homologs, in each of three AMS-preferential ortholog groups, had some local variations in the protein structural alignment. We also identified genes co-expressed with the Populus trichocarpa SSP genes in the AMS-preferential ortholog groups. This first plant kingdom-wide analysis on SSP provides insights on plant-AMS convergent evolution with specific SSP gene expression and local diversification of protein structures.

New Insights into the Chloroplast Outer Membrane Proteome and Associated Targeting Pathways

Plastids are a dynamic class of organelle in plant cells that arose from an ancient cyanobacterial endosymbiont. Over the course of evolution, most genes encoding plastid proteins were transferred to the nuclear genome. In parallel, eukaryotic cells evolved a series of targeting pathways and complex proteinaceous machinery at the plastid surface to direct these proteins back to their target organelle. Chloroplasts are the most well-characterized plastids, responsible for photosynthesis and other important metabolic functions. The biogenesis and function of chloroplasts rely heavily on the fidelity of intracellular protein trafficking pathways. Therefore, understanding these pathways and their regulation is essential. Furthermore, the chloroplast outer membrane proteome remains relatively uncharted territory in our understanding of protein targeting. Many key players in the cytosol, receptors at the organelle surface, and insertases that facilitate insertion into the chloroplast outer membrane remain elusive for this group of proteins. In this review, we summarize recent advances in the understanding of well-characterized chloroplast outer membrane protein targeting pathways as well as provide new insights into novel targeting signals and pathways more recently identified using a bioinformatic approach. As a result of our analyses, we expand the known number of chloroplast outer membrane proteins from 117 to 138.

Identification of Key Amino Acids that Impact Organic Solute Transporter α/β (OSTα/β)

Organic solute transporter α/β (OSTα/β) is a bidirectional bile acid transporter localized on the basolateral membrane of hepatic, intestinal, and renal epithelial cells. OSTα/β plays a critical role in intestinal bile acid reabsorption and is upregulated in hepatic diseases characterized by elevated bile acids, whereas genetic variants in SLC51A/B have been associated with clinical cholestasis. OSTα/β also transports and is inhibited by commonly used medications. However, there is currently no high-resolution structure of OSTα/β, and structure-function data for OSTα, the proposed substrate-binding subunit, are lacking. The present study addressed this knowledge gap and identified amino acids in OSTα that are important for bile acid transport. This was accomplished using computational modeling and site-directed mutagenesis of the OSTα subunit to generate OSTα/β mutant cell lines. Out of the 10 OSTα/β mutants investigated, four (S228K, T229S, Q269E, Q269K) exhibited decreased [3H]-taurocholate (TCA) uptake (ratio of geometric means relative to OSTα/β wild type (WT) of 0.76, 0.75, 0.79, and 0.13, respectively). Three OSTα/β mutants (S228K, Q269K, E305A) had reduced [3H]-TCA efflux % (ratio of geometric means relative to OSTα/β WT of 0.86, 0.65, and 0.79, respectively). Additionally, several OSTα/β mutants demonstrated altered expression and cellular localization when compared with OSTα/β WT. In summary, we identified OSTα residues (Ser228, Thr229, Gln269, Glu305) in predicted transmembrane domains that affect expression of OSTα/β and may influence OSTα/β-mediated bile acid transport. These data advance our understanding of OSTα/β structure/function and can inform future studies designed to gain further insight into OSTα/β structure or to identify additional OSTα/β substrates and inhibitors. SIGNIFICANCE STATEMENT: OSTα/β is a clinically important transporter involved in enterohepatic bile acid recycling with currently no high-resolution protein structure and limited structure-function data. This study identified four OSTα amino acids (Ser228, Thr229, Gln269, Glu305) that affect expression of OSTα/β and may influence OSTα/β-mediated bile acid transport. These data can be utilized to inform future investigation of OSTα/β structure and refine molecular modeling approaches to facilitate the identification of substrates and/or inhibitors of OSTα/β.

Titratable transmembrane residues and a hydrophobic plug are essential for manganese import via the Bacillus anthracis ABC transporter MntBC-A

All extant life forms require trace transition metals (e.g., Fe^2/3+, Cu^1/2+, and Mn²⁺) to survive. However, as these are environmentally scarce, organisms have evolved sophisticated metal uptake machineries. In bacteria, high-affinity import of transition metals is predominantly mediated by ABC transporters. During bacterial infection, sequestration of metal by the host further limits the availability of these ions, and accordingly, bacterial ABC transporters (importers) of metals are key virulence determinants. However, the structure–function relationships of these metal transporters have not been fully elucidated. Here, we used metal-sensitivity assays, advanced structural modeling, and enzymatic assays to study the ABC transporter MntBC-A, a virulence determinant of the bacterial human pathogen Bacillus anthracis. We find that despite its broad metal-recognition profile, MntBC-A imports only manganese, whereas zinc can function as a high-affinity inhibitor of MntBC-A. Computational analysis shows that the transmembrane metal permeation pathway is lined with six titratable residues that can coordinate the positively charged metal, and mutagenesis studies show that they are essential for manganese transport. Modeling suggests that access to these titratable residues is blocked by a ladder of hydrophobic residues, and ATP-driven conformational changes open and close this hydrophobic seal to permit metal binding and release. The conservation of this arrangement of titratable and hydrophobic residues among ABC transporters of transition metals suggests a common mechanism. These findings advance our understanding of transmembrane metal recognition and permeation and may aid the design and development of novel antibacterial agents.

Calicivirus Non-structural Proteins: Potential Functions in Replication and Host Cell Manipulation

The Caliciviridae are a family of viruses with a single-stranded, non-segmented RNA genome of positive polarity. The ongoing discovery of caliciviruses has increased the number of genera in this family to 11 (Norovirus, Nebovirus, Sapovirus, Lagovirus, Vesivirus, Nacovirus, Bavovirus, Recovirus, Salovirus, Minovirus, and Valovirus). Caliciviruses infect a wide range of hosts that include fishes, amphibians, reptiles, birds, and marine and land mammals. All caliciviruses have a genome that encodes a major and a minor capsid protein, a genome-linked viral protein, and several non-structural proteins. Of these non-structural proteins, only the helicase, protease, and RNA-dependent RNA polymerase share clear sequence and structural similarities with proteins from other virus families. In addition, all caliciviruses express two or three non-structural proteins for which functions have not been clearly defined. The sequence diversity of these non-structural proteins and a multitude of processing strategies suggest that at least some have evolved independently, possibly to counteract innate and adaptive immune responses in a host-specific manner. Studying these proteins is often difficult as many caliciviruses cannot be grown in cell culture. Nevertheless, the study of recombinant proteins has revealed many of their properties, such as intracellular localization, capacity to oligomerize, and ability to interact with viral and/or cellular proteins; the release of non-structural proteins from transfected cells has also been investigated. Here, we will summarize these findings and discuss recent in silico studies that identified previously overlooked putative functional domains and structural features, including transmembrane domains that suggest the presence of viroporins.

Advances and perspectives in discovery and functional analysis of small secreted proteins in plants

Small secreted proteins (SSPs) are less than 250 amino acids in length and are actively transported out of cells through conventional protein secretion pathways or unconventional protein secretion pathways. In plants, SSPs have been found to play important roles in various processes, including plant growth and development, plant response to abiotic and biotic stresses, and beneficial plant-microbe interactions. Over the past 10 years, substantial progress has been made in the identification and functional characterization of SSPs in several plant species relevant to agriculture, bioenergy, and horticulture. Yet, there are potentially a lot of SSPs that have not been discovered in plant genomes, which is largely due to limitations of existing computational algorithms. Recent advances in genomics, transcriptomics, and proteomics research, as well as the development of new computational algorithms based on machine learning, provide unprecedented capabilities for genome-wide discovery of novel SSPs in plants. In this review, we summarize known SSPs and their functions in various plant species. Then we provide an update on the computational and experimental approaches that can be used to discover new SSPs. Finally, we discuss strategies for elucidating the biological functions of SSPs in plants.

Modeling of full-length Piezo1 suggests importance of the proximal N-terminus for dome structure

Piezo1 forms a mechanically activated calcium-permeable nonselective cation channel that is functionally important in many cell types. Structural data exist for C-terminal regions, but we lack information about N-terminal regions and how the entire channel interacts with the lipid bilayer. Here, we use computational approaches to predict the three-dimensional structure of the full-length Piezo1 and simulate it in an asymmetric membrane. A number of novel insights are suggested by the model: 1) Piezo1 creates a trilobed dome in the membrane that extends beyond the radius of the protein, 2) Piezo1 changes the lipid environment in its vicinity via preferential interactions with cholesterol and phosphatidylinositol 4,5-bisphosphate (PIP₂) molecules, and 3) cholesterol changes the depth of the dome and PIP₂ binding preference. In vitro alteration of cholesterol concentration inhibits Piezo1 activity in a manner complementing some of our computational findings. The data suggest the importance of N-terminal regions of Piezo1 for dome structure and membrane cholesterol and PIP₂ interactions.

Sequence Features of Mitochondrial Transporter Protein Families

Mitochondrial carriers facilitate the transfer of small molecules across the inner mitochondrial membrane (IMM) to support mitochondrial function and core cellular processes. In addition to the classical SLC25 (solute carrier family 25) mitochondrial carriers, the past decade has led to the discovery of additional protein families with numerous members that exhibit IMM localization and transporter-like properties. These include mitochondrial pyruvate carriers, sideroflexins, and mitochondrial cation/H+ exchangers. These transport proteins were linked to vital physiological functions and disease. Their structures and transport mechanisms are, however, still largely unknown and understudied. Protein sequence analysis per se can often pinpoint hotspots that are of functional or structural importance. In this review, we summarize current knowledge about the sequence features of mitochondrial transporters with a special focus on the newly included SLC54, SLC55 and SLC56 families of the SLC solute carrier superfamily. Taking a step further, we combine sequence conservation analysis with transmembrane segment and secondary structure prediction methods to extract residue positions and sequence motifs that likely play a role in substrate binding, binding site gating or structural stability. We hope that our review will help guide future experimental efforts by the scientific community to unravel the transport mechanisms and structures of these novel mitochondrial carriers.

The Puzzle of Metabolite Exchange and Identification of Putative Octotrico Peptide Repeat Expression Regulators in the Nascent Photosynthetic Organelles of Paulinella chromatophora

The endosymbiotic acquisition of mitochondria and plastids more than one billion years ago was central for the evolution of eukaryotic life. However, owing to their ancient origin, these organelles provide only limited insights into the initial stages of organellogenesis. The cercozoan amoeba Paulinella chromatophora contains photosynthetic organelles-termed chromatophores-that evolved from a cyanobacterium ∼100 million years ago, independently from plastids in plants and algae. Despite the more recent origin of the chromatophore, it shows tight integration into the host cell. It imports hundreds of nucleus-encoded proteins, and diverse metabolites are continuously exchanged across the two chromatophore envelope membranes. However, the limited set of chromatophore-encoded solute transporters appears insufficient for supporting metabolic connectivity or protein import. Furthermore, chromatophore-localized biosynthetic pathways as well as multiprotein complexes include proteins of dual genetic origin, suggesting that mechanisms evolved that coordinate gene expression levels between chromatophore and nucleus. These findings imply that similar to the situation in mitochondria and plastids, also in P. chromatophora nuclear factors evolved that control metabolite exchange and gene expression in the chromatophore. Here we show by mass spectrometric analyses of enriched insoluble protein fractions that, unexpectedly, nucleus-encoded transporters are not inserted into the chromatophore inner envelope membrane. Thus, despite the apparent maintenance of its barrier function, canonical metabolite transporters are missing in this membrane. Instead we identified several expanded groups of short chromatophore-targeted orphan proteins. Members of one of these groups are characterized by a single transmembrane helix, and others contain amphipathic helices. We hypothesize that these proteins are involved in modulating membrane permeability. Thus, the mechanism generating metabolic connectivity of the chromatophore fundamentally differs from the one for mitochondria and plastids, but likely rather resembles the poorly understood mechanism in various bacterial endosymbionts in plants and insects. Furthermore, our mass spectrometric analysis revealed an expanded family of chromatophore-targeted helical repeat proteins. These proteins show similar domain architectures as known organelle-targeted expression regulators of the octotrico peptide repeat type in algae and plants. Apparently these chromatophore-targeted proteins evolved convergently to plastid-targeted expression regulators and are likely involved in gene expression control in the chromatophore.

Insights on the Quest for the Structure-Function Relationship of the Mitochondrial Pyruvate Carrier

Simple Summary

The atomic structure of a biological macromolecule determines its function. Discovering how one or more amino acid chains fold and interact to form a protein complex is critical, from understanding the most fundamental cellular processes to developing new therapies. However, this is far from a straightforward task, especially when studying a membrane protein. The functional link between the oligomeric state and complex composition of the proteins involved in the active mitochondrial transport of cytosolic pyruvate is a decades-old question but remains urgent. We present a brief historical review beginning with the identification of the so-called mitochondrial pyruvate carrier (MPC) proteins, followed by a rigorous conceptual analysis of technical approaches in more recent biochemical studies that seek to isolate and reconstitute the functional MPC complex(es) in vitro. We correlate these studies with early kinetic observations and current experimental and computational knowledge to assess their main contributions, identify gaps, resolve ambiguities, and better define the research goal.

Abstract

The molecular identity of the mitochondrial pyruvate carrier (MPC) was presented in 2012, forty years after the active transport of cytosolic pyruvate into the mitochondrial matrix was first demonstrated. An impressive amount of in vivo and in vitro studies has since revealed an unexpected interplay between one, two, or even three protein subunits defining different functional MPC assemblies in a metabolic-specific context. These have clear implications in cell homeostasis and disease, and on the development of future therapies. Despite intensive efforts by different research groups using state-of-the-art computational tools and experimental techniques, MPCs’ structure-based mechanism remains elusive. Here, we review the current state of knowledge concerning MPCs’ molecular structures by examining both earlier and recent studies and presenting novel data to identify the regulatory, structural, and core transport activities to each of the known MPC subunits. We also discuss the potential application of cryogenic electron microscopy (cryo-EM) studies of MPC reconstituted into nanodiscs of synthetic copolymers for solving human MPC2.

Prediction of Ion Channels and their Types from Protein Sequences: Comprehensive Review and Comparative Assessment

BACKGROUND

Ion channels are a large and growing protein family. Many of them are associated with diseases, and consequently, they are targets for over 700 drugs. Discovery of new ion channels is facilitated with computational methods that predict ion channels and their types from protein sequences. However, these methods were never comprehensively compared and evaluated.

OBJECTIVE

We offer first-of-its-kind comprehensive survey of the sequence-based predictors of ion channels. We describe eight predictors that include five methods that predict ion channels, their types, and four classes of the voltage-gated channels. We also develop and use a new benchmark dataset to perform comparative empirical analysis of the three currently available predictors.

RESULTS

While several methods that rely on different designs were published, only a few of them are currently available and offer a broad scope of predictions. Support and availability after publication should be required when new methods are considered for publication. Empirical analysis shows strong performance for the prediction of ion channels and modest performance for the prediction of ion channel types and voltage-gated channel classes. We identify a substantial weakness of current methods that cannot accurately predict ion channels that are categorized into multiple classes/types.

CONCLUSION

Several predictors of ion channels are available to the end users. They offer practical levels of predictive quality. Methods that rely on a larger and more diverse set of predictive inputs (such as PSIONplus) are more accurate. New tools that address multi-label prediction of ion channels should be developed.

New Insights into the Structure-Function Relationship of the Endosomal-Type Na+, K+/H+ Antiporter NHX6 from Mulberry (Morus notabilis)

The endosomal-type Na⁺, K⁺/H⁺ antiporters (NHXs) play important roles in K⁺, vesicle pH homeostasis, and protein trafficking in plant. However, the structure governing ion transport mechanism and the key residues related to the structure–function of the endosomal-type NHXs remain unclear. Here, the structure-function relationship of the only endosomal-type NHX from mulberry, MnNHX6, was investigated by homology modeling, mutagenesis, and localization analyses in yeast. The ectopic expression of MnNHX6 in arabidopsis and Nhx1 mutant yeast can enhance their salt tolerance. MnNHX6’s three-dimensional structure, established by homology modeling, was supported by empirical, phylogenetic, and experimental data. Structure analysis showed that MnNHX6 contains unusual 13 transmembrane helices, but the structural core formed by TM5-TM12 assembly is conserved. Localization analysis showed that MnNHX6 has the same endosomal localization as yeast Nhx1/VPS44, and Arg402 is important for protein stability of MnNHX6. Mutagenesis analysis demonstrated MnNHX6 contains a conserved cation binding mechanism and a similar charge-compensated pattern as NHE1, but shares a different role in ion selectivity than the vacuolar-type NHXs. These results improve our understanding of the role played by the structure–function related key residues of the plant endosomal-type NHXs, and provide a basis for the ion transport mechanism study of endosomal-type NHXs.

Prediction and validation of the structural features of Ov58GPCR, an immunogenic determinant of Onchocerca volvulus

Onchocerciasis is a severely debilitating yet neglected tropical disease (NTD) that creates social stigma, generates and perpetuates poverty, and leads ultimately in some cases to irreversible unilateral or bilateral blindness if untreated. Consequently, the disease is a major impediment to socioeconomic development. Many control programs have been launched for the disease with moderate successes achieved. This mitigated hit is partially due to the lingering need for reliable, non-invasive and easily applicable tools for mapping endemic regions and post-elimination surveillance. In this work, bioinformatics analyses combined with immunological assays were applied in a bid to develop potential tools for diagnosis and assessing the success of drug treatment programs. We report that (i) the O. volvulus antigen, Ov58GPCR is a G-protein coupled receptor (GPCR) conserved in related nematodes, (ii) synthetic peptides predicted to be in the extracellular domain (ECD) of Ov58GPCR are indeed immunogenic epitopes in actively-infected individuals, (iii) synthetic peptide cocktails discriminate between actively-infected individuals, treated individuals and healthy African controls, (iv) polyclonal antibodies against one of the peptides or against the bacterially-expressed ECD reacted specifically with the native antigen of O. volvulus total and surface extracts, (v) Ov58GPCR is transcribed in both larvae and adult parasite stages, (vi) IgG and IgE responses to the recombinant ECD decline with ivermectin treatment. All these findings suggest that the extracellular domain and synthetic peptides of Ov58GPCR, as well as the specific immune response generated could be harnessed in the context of disease diagnosis and surveillance.

Aluminum-Activated Malate Transporters Can Facilitate GABA Transport

Plant aluminum-activated malate transporters (ALMTs) are currently classified as anion channels; they are also known to be regulated by diverse signals, leading to a range of physiological responses. Gamma-aminobutyric acid (GABA) regulation of anion flux through ALMT proteins requires a specific amino acid motif in ALMTs that shares similarity with a GABA binding site in mammalian GABA_A receptors. Here, we explore why TaALMT1 activation leads to a negative correlation between malate efflux and endogenous GABA concentrations ([GABA]_i) in both wheat (Triticum aestivum) root tips and in heterologous expression systems. We show that TaALMT1 activation reduces [GABA]_i because TaALMT1 facilitates GABA efflux but GABA does not complex Al³⁺ TaALMT1 also leads to GABA transport into cells, demonstrated by a yeast complementation assay and via ¹⁴C-GABA uptake into TaALMT1-expressing Xenopus laevis oocytes; this was found to be a general feature of all ALMTs we examined. Mutation of the GABA motif (TaALMT1^F213C) prevented both GABA influx and efflux, and resulted in no correlation between malate efflux and [GABA]_i We conclude that ALMTs are likely to act as both GABA and anion transporters in planta. GABA and malate appear to interact with ALMTs in a complex manner to regulate each other's transport, suggestive of a role for ALMTs in communicating metabolic status.

Topological analysis of DPY19L3, a human C-mannosyltransferase

C-mannosylation is a rare type of protein glycosylation, the functions and mechanisms of which remain unclear. Recently, we identified DPY19L3 as a C-mannosyltransferase of R-spondin1 in human cells. DPY19L3 is predicted to be a multipass transmembrane protein that localizes to the endoplasmic reticulum (ER); however, its structure is undetermined. In this study, we propose a topological structure of DPY19L3 by in silico analysis and experimental methods such as redox-sensitive luciferase assay and introduction of N-glycosylation sites, suggesting that DPY19L3 comprises 11 transmembrane regions and two re-entrant loops with the N- and C-terminal ends facing the cytoplasm and ER lumen, respectively. Furthermore, DPY19L3 has four predicted N-glycosylation sites, and we have demonstrated that DPY19L3 is N-glycosylated at Asn¹¹⁸ and Asn⁷⁰⁴ but not Asn³¹⁹ and Asn⁴³⁹ , supporting our topological model. By mass spectrometry, we measured the C-mannosyltransferase activity of N-glycosylation-defective mutants of DPY19L3 and isoform2, a splice variant, which lacks the C-terminal luminal region of DPY19L3. Isoform2 does not possess C-mannosyltransferase activity, indicating the importance of the C-terminal region; however, N-glycosylations of DPY19L3 do not have any roles for its enzymatic activity. These novel findings on DPY19L3 provide important insights into the mechanism of C-mannosylation.

Membrane topology and identification of key residues of EaDAcT, a plant MBOAT with unusual substrate specificity

Euonymus alatus diacylglycerol acetyltransferase (EaDAcT) catalyzes the transfer of an acetyl group from acetyl-CoA to the sn-3 position of diacylglycerol to form 3-acetyl-1,2-diacyl-sn-glycerol (acetyl-TAG). EaDAcT belongs to a small, plant-specific subfamily of the membrane bound O-acyltransferases (MBOAT) that acylate different lipid substrates. Sucrose gradient density centrifugation revealed that EaDAcT colocalizes to the same fractions as an endoplasmic reticulum (ER)-specific marker. By mapping the membrane topology of EaDAcT, we obtained an experimentally determined topology model for a plant MBOAT. The EaDAcT model contains four transmembrane domains (TMDs), with both the N- and C-termini orientated toward the lumen of the ER. In addition, there is a large cytoplasmic loop between the first and second TMDs, with the MBOAT signature region of the protein embedded in the third TMD close to the interface between the membrane and the cytoplasm. During topology mapping, we discovered two cysteine residues (C187 and C293) located on opposite sides of the membrane that are important for enzyme activity. In order to identify additional amino acid residues important for acetyltransferase activity, we isolated and characterized acetyltransferases from other acetyl-TAG-producing plants. Among them, the acetyltransferase from Euonymus fortunei possessed the highest activity in vivo and in vitro. Mutagenesis of conserved amino acids revealed that S253, H257, D258 and V263 are essential for EaDAcT activity. Alteration of residues unique to the acetyltransferases did not alter the unique acyl donor specificity of EaDAcT, suggesting that multiple amino acids are important for substrate recognition.

The elusive MAESTRO gene: Its human reproductive tissue-specific expression pattern

The encoded transcript of the Maestro—Male-specific Transcription in the developing Reproductive Organs (MRO) gene exhibits sexual dimorphic expression during murine gonadal development. The gene has no homology to any known gene and its expression pattern, protein function or structure are still unknown. Previously, studying gene expression in human ovarian cumulus cells, we found increased expression of MRO in lean-type Polycystic Ovarian Syndrome (PCOS) subjects, as compared to controls. In this study, we examined the MRO splice variants and protein expression pattern in various human tissues and cells. We found a differential expression pattern of the MRO 5’-UTR region in luteinized granulosa-cumulus cells and in testicular tissues as compared to non-gonadal tissues. Our study also shows a punctate nuclear expression pattern and disperse cytoplasmic expression pattern of the MRO protein in human granulosa-cumulus cells and in testicular germ cells, which was later validated by western blotting. The tentative and unique features of the protein hampered our efforts to gain more insight about this elusive protein. A better understanding of the tissue-specific MRO isoforms expression patterns and the unique structure of the protein may provide important insights into the function of this gene and possibly to the pathophysiology of PCOS.

Characterization of Ghrelin O-Acyltransferase (GOAT) in goldfish (Carassius auratus)

Ghrelin is the only known hormone posttranslationally modified with an acylation. This modification is crucial for most of ghrelin’s physiological effects and is catalyzed by the polytopic enzyme ghrelin O-acyltransferase (GOAT). The aim of this study was to characterize GOAT in a teleost model, goldfish (Carassius auratus). First, the full-length cDNA sequence was obtained by RT-PCR and rapid amplification of cDNA ends methods. Two highly homologous cDNAs of 1491 and 1413 bp, respectively, named goat-V1 and goat-V2 were identified. Deduced protein sequences (393 and 367 amino acids, respectively) are predicted to present 11 and 9 transmembrane regions, respectively, and both contain two conserved key residues proposed to be involved in catalysis: asparagine 273 and histidine 304. RT-qPCR revealed that both forms of goat mRNAs show a similar widespread tissue distribution, with the highest expression in the gastrointestinal tract and gonads and less but considerable expression in brain, pituitary, liver and adipose tissue. Immunostaining of intestinal sections showed the presence of GOAT immunoreactive cells in the intestinal mucosa, some of which colocalize with ghrelin. Using an in vitro approach, we observed that acylated ghrelin downregulates GOAT gene and protein levels in cultured intestine in a time-dependent manner. Finally, we found a rhythmic oscillation of goat mRNA expression in the hypothalamus, pituitary and intestinal bulb of goldfish fed at midday, but not at midnight. Together, these findings report novel data characterizing GOAT, and offer new information about the ghrelinergic system in fish.

Computational comparison of a calcium-dependent jellyfish protein (apoaequorin) and calmodulin-cholesterol in short-term memory maintenance

Memory reconsolidation and maintenance depend on calcium channels and on calcium/calmodulin-dependent kinases regulating protein turnover in the hippocampus. Ingestion of a jellyfish protein, apoaequorin, reportedly protects and/or improves verbal learning in adults and is currently widely advertised for use by the elderly. Apoaequorin is a member of the EF-hand calcium binding family of proteins that includes calmodulin. Calmodulin-1 (148 residues) differs from Apoaequorin (195 residues) in that it contains four rather than three Ca²⁺-binding sites and three rather than four cholesterol-binding (CRAC, CARC) domains. All three cholesterol-binding CARC domains in calmodulin have a high interaction affinity for cholesterol compared to only two high affinity CARC domains in apoaequorin. Both calmodulin and apoaequorin can form dimers with a potential of eight bound Ca²⁺ ions and six high affinity-bound cholesterol molecules in calmodulin with six bound Ca²⁺ ions and a mixed population of eight cholesterols bound to both CARC and CRAC domains in apoaqueorin. MEMSAT-SVM analysis indicates that both calmodulin and apoaqueorin have a pore-lining region. The Peptide-Cutter algorithm predicts that calmodulin-1 contains 11 trypsin-specific cleavage sites (compared to 21 in apoaqueorin), four of which are potentially blocked by cholesterol and three are within the Ca-binding domains and/or the pore-lining region. Three are clustered between the third and fourth Ca²⁺-binding sites. Only calmodulin pore-lining regions contain Ca²⁺ binding sites and as dimers may insert into the plasma membrane of neural cells and act as Ca²⁺ channels. In a dietary supplement, bound cholesterol may protect both apoaequorin and calmodulin from proteolysis in the gut as well as facilitate uptake across the blood-brain barrier. Our results suggest that a physiological calmodulin-cholesterol complex, not cholesterol-free jellyfish protein, may better serve as a dietary supplement to facilitate memory maintenance.

Studies on the Virome of the Entomopathogenic Fungus Beauveria bassiana Reveal Novel dsRNA Elements and Mild Hypervirulence

The entomopathogenic fungus Beauveria bassiana has a wide host range and is used as a biocontrol agent against arthropod pests. Mycoviruses have been described in phytopathogenic fungi while in entomopathogenic fungi their presence has been reported only rarely. Here we show that 21.3% of a collection of B. bassiana isolates sourced from worldwide locations, harbor dsRNA elements. Molecular characterization of these elements revealed the prevalence of mycoviruses belonging to the Partitiviridae and Totiviridae families, the smallest reported virus to date, belonging to the family Narnaviridae, and viruses unassigned to a family or genus. Of particular importance is the discovery of members of a newly proposed family Polymycoviridae in B. bassiana. Polymycoviruses, previously designated as tetramycoviruses, consist of four non-conventionally encapsidated capped dsRNAs. The presence of additional non-homologous genomic segments in B. bassiana polymycoviruses and other fungi illustrates the unprecedented dynamic nature of the viral genome. Finally, a comparison of virus-free and virus-infected isogenic lines derived from an exemplar B. bassiana isolate revealed a mild hypervirulent effect of mycoviruses on the growth of their host isolate and on its pathogenicity against the greater wax moth Galleria mellonella, highlighting for the first time the potential of mycoviruses as enhancers of biocontrol agents.

Molecular Properties of Globin Channels and Pores: Role of Cholesterol in Ligand Binding and Movement

Globins contain one or more cavities that control or affect such functions as ligand movement and ligand binding. Here we report that the extended globin family [cytoglobin (Cygb); neuroglobin (Ngb); myoglobin (Mb); hemoglobin (Hb) subunits Hba(α); and Hbb(β)] contain either a transmembrane (TM) helix or pore-lining region as well as internal cavities. Protein motif/domain analyses indicate that Ngb and Hbb each contain 5 cholesterol- binding (CRAC/CARC) domains and 1 caveolin binding motif, whereas the Cygb dimer has 6 cholesterol-binding domains but lacks caveolin-binding motifs. Mb and Hba each exhibit 2 cholesterol-binding domains and also lack caveolin-binding motifs. The Hb αβ-tetramer contains 14 cholesterol-binding domains. Computer algorithms indicate that Cygb and Ngb cavities display multiple partitions and C-terminal pore-lining regions, whereas Mb has three major cavities plus a C-terminal pore-lining region. The Hb tetramer exhibits a large internal cavity but the subunits differ in that they contain a C-terminal TM helix (Hba) and pore-lining region (Hbb). The cavities include 43 of 190 Cygb residues, 38 of 151 of Ngb residues, 55 of 154 Mb residues, and 137 of 688 residues in the Hb tetramer. Each cavity complex includes 6 to 8 residues of the TM helix or pore-lining region and CRAC/CARC domains exist within all cavities. Erythrocyte Hb αβ-tetramers are largely cytosolic but also bind to a membrane anion exchange protein, "band 3," which contains a large internal cavity and 12 TM helices (5 being pore-lining regions). The Hba TM helix may be the erythrocyte membrane "band 3" attachment site. "Band 3" contributes 4 caveolin binding motifs and 10 CRAC/CARC domains. Cholesterol binding may create lipid-disordered phases that alter globin cavities and facilitate ligand movement, permitting ion channel formation and conformational changes that orchestrate anion and ligand (O2, CO2, NO) movement within the large internal cavities and channels of the globins.

A putative amino acid transporter determines sensitivity to the two-peptide bacteriocin plantaricin JK

Lactobacillus plantarum produces a number of antimicrobial peptides (bacteriocins) that mostly target closely related bacteria. Although bacteriocins are important for the ecology of these bacteria, very little is known about how the peptides target sensitive cells. In this work, a putative membrane protein receptor of the two-peptide bacteriocin plantaricin JK was identified by comparing Illumina sequence reads from plantaricin JK-resistant mutants to a crude assembly of the sensitive wild-type Weissella viridescens genome using the polymorphism discovery tool VAAL. Ten resistant mutants harbored altogether seven independent mutations in a gene encoding an APC superfamily protein with 12 transmembrane helices. The APC superfamily transporter thus is likely to serve as a target for plantaricin JK on sensitive cells.

Evolution of the α-Subunit of Na/K-ATPase from Paramecium to Homo sapiens: Invariance of Transmembrane Helix Topology

Na/K-ATPase is a key plasma membrane enzyme involved in cell signaling, volume regulation, and maintenance of electrochemical gradients. The α-subunit, central to these functions, belongs to a large family of P-type ATPases. Differences in transmembrane (TM) helix topology, sequence homology, helix–helix contacts, cell signaling, and protein domains of Na/K-ATPase α-subunit were compared in fungi (Beauveria), unicellular organisms (Paramecia), primitive multicellular organisms (Hydra), and vertebrates (Xenopus, Homo sapiens), and correlated with evolution of physiological functions in the α-subunit. All α-subunits are of similar length, with groupings of four and six helices in the N- and C-terminal regions, respectively. Minimal homology was seen for protein domain patterns in Paramecium and Hydra, with high correlation between Hydra and vertebrates. Paramecium α-subunits display extensive disorder, with minimal helix contacts. Increases in helix contacts in Hydra approached vertebrates. Protein motifs known to be associated with membrane lipid rafts and cell signaling reveal significant positional shifts between Paramecium and Hydra vulgaris, indicating that regional membrane fluidity changes occur during evolution. Putative steroid binding sites overlapping TM-3 occurred in all species. Sites associated with G-protein-receptor stimulation occur both in vertebrates and amphibia but not in Hydra or Paramecia. The C-terminus moiety “KETYY,” necessary for the Na⁺ activation of pump phosphorylation, is not present in unicellular species indicating the absence of classical Na⁺/K⁺-pumps. The basic protein topology evolved earliest, followed by increases in protein domains and ordered helical arrays, correlated with appearance of α-subunit regions known to involve cell signaling, membrane recycling, and ion channel formation.

Plant abiotic stress-related RCI2/PMP3s: multigenes for multiple roles

RCI2 / PMP3 s participate in abiotic stress responses and impact the expression of other genes. Their multifunctionality is determined by differential expression and by distinct activities of their structurally different proteins. In plants, RCI2/PMP3 genes, which encode small membrane proteins of the PMP3 family, are closely associated with abiotic stress responses. Their involvement in mediating stress tolerance is supported by genetic evidence and overexpression studies. RCI2/PMP3s occur as multigenes in plant genomes and their encoded proteins belong to distinct and conserved structural groups. In addition, different isoforms appear to be targeted to the plasma membrane or to distinct endomembrane compartments in cells. Several studies have revealed that RCI2/PMP3 proteins participate in cell ion homeostasis, and in regulation of membrane stability and polarization. They also appear to potentiate plant transcriptional responses to abiotic stresses. However, their mechanisms of action remain unknown. This paper reviews the current knowledge of the multiple roles of plant RCI2/PMP3 genes resulting from their differential expression under normal and stress conditions. The structural diversity of RCI2/PMP3 proteins is analyzed and evidence supporting their functional specialization and possible activity mechanisms is examined. Finally, strategies are discussed for exploiting new and established technologies to overcome the difficulties posed by the multigene status of RCI2s and the integral membrane character of their proteins, enabling the probing of their individual functions and collective significance.

TMEM175 Is an Organelle K(+) Channel Regulating Lysosomal Function

Potassium is the most abundant ion to face both plasma and organelle membranes. Extensive research over the past seven decades has characterized how K(+) permeates the plasma membrane to control fundamental processes such as secretion, neuronal communication, and heartbeat. However, how K(+) permeates organelles such as lysosomes and endosomes is unknown. Here, we directly recorded organelle K(+) conductance and discovered a major K(+)-selective channel KEL on endosomes and lysosomes. KEL is formed by TMEM175, a protein with unknown function. Unlike any of the ∼80 plasma membrane K(+) channels, TMEM175 has two repeats of 6-transmembrane-spanning segments and has no GYG K(+) channel sequence signature-containing, pore-forming P loop. Lysosomes lacking TMEM175 exhibit no K(+) conductance, have a markedly depolarized ΔΨ and little sensitivity to changes in [K(+)], and have compromised luminal pH stability and abnormal fusion with autophagosomes during autophagy. Thus, TMEM175 comprises a K(+) channel that underlies the molecular mechanism of lysosomal K(+) permeability.

The Regulatory Domain of Squalene Monooxygenase Contains a Re-entrant Loop and Senses Cholesterol via a Conformational Change

Squalene monooxygenase (SM) is an important control point in cholesterol synthesis beyond 3-hydroxy-3-methylglutaryl-CoA reductase. Although it is known to associate with the endoplasmic reticulum, its topology has not been determined. We have elucidated the membrane topology of the sterol-responsive domain of SM comprising the first 100 amino acids fused to GFP (SM N100-GFP) by determining the accessibility of 16 introduced cysteines to the cysteine-reactive, membrane-impermeable reagent PEG-maleimide. We have identified a region integrally associated with the endoplasmic reticulum membrane that is likely to interact with cholesterol or respond to cholesterol-induced membrane effects. By comparing cysteine accessibility with and without cholesterol treatment, we further present evidence to suggest that cholesterol induces a conformational change in SM N100-GFP. This change is likely to lead to its targeted degradation by the ubiquitin-proteasome system because degradation is blunted by treatment with the chemical chaperone glycerol, which retains SM N100-GFP in its native conformation. Furthermore, degradation can be disrupted by insertion of two N-terminal myc tags, implicating the N terminus in this process. Together, this information provides new molecular insights into the regulation of this critical control point in cholesterol synthesis.

Transmembrane helices in "classical" nuclear reproductive steroid receptors: a perspective

Abstract Steroid receptors of the nuclear receptor superfamily are proposed to be either: 1) located in the cytosol and moved to the cell nucleus upon activation, 2) tethered to the inside of the plasma membrane, or 3) retained in the nucleus until free steroid hormone enters and activates specific receptors. Using computational methods to analyze peptide receptor topology, we find that the “classical” nuclear receptors for progesterone (PRB/PGR), androgen (ARB/AR) and estrogen (ER1/ESR1) contain two transmembrane helices (TMH) within their ligand-binding domains (LBD).The MEMSAT-SVM algorithm indicates that ARB and ER2 (but not PRB or ER1) contain a pore-lining (channel-forming) region which may merge with other pore-lining regions to form a membrane channel. ER2 lacks a TMH, but contains a single pore-lining region. The MemBrain algorithm predicts that PRB, ARB and ER1 each contain one TMH plus a half TMH separated by 51 amino acids.ER2 contains two half helices. The TM-2 helices of ARB, ER1 and ER2 each contain 9-13 amino acid motifs reported to translocate the receptor to the plasma membrane, as well as cysteine palmitoylation sites. PoreWalker analysis of X-ray crystallographic data identifies a pore or channel within the LBDs of ARB and ER1 and predicts that 70 and 72 residues are pore-lining residues, respectively. The data suggest that (except for ER2), cytosolic receptors become anchored to the plasma membrane following synthesis. Half-helices and pore-lining regions in turn form functional ion channels and/or facilitate passive steroid uptake into the cell. In perspective, steroid-dependent insertion of “classical” receptors containing pore-lining regions into the plasma membrane may regulate permeability to ions such as Ca²⁺, Na⁺ or K⁺, as well as facilitate steroid translocation into the nucleus.

Computational analysis of the extracellular domain of the Ca²⁺-sensing receptor: an alternate model for the Ca²⁺ sensing region

The extracellular Ca(2+) sensing receptor (CaSR) belongs to Class C G-protein-coupled receptors (GPCRs) which include receptors for amino acids, γ-aminobutyric acid and glutamate neurotransmitters. CaSR has been described as having an extended sequence containing a Ca(2+) binding pocket within an extracellular amino (N)-terminal domain, called a Venus Fly Trap (VFT) module. CaSR is thought to consist of three domains: 1) a Ca(2+-)sensory domain, 2) a region containing 7 transmembrane (TM) helices, and 3) a carboxy (C)-terminal tail. We find that SPOCTOPUS (a combination of hidden Markov models and artificial neural networks) predicts that Homo sapiens CaSR contains two additional TM helices ((190)D - G(210); (262)S-E(282)), with the second TM helix containing a pore-lining region ((265)K - I(280)). This predicts that the putative Ca(2+) sensory domain is within an extracellular loop, N-terminal to the highly conserved heptahelical bundle. This loop contains both the cysteine-rich domain ((537)V - C(598)) and a 14 residue "linker" sequence ((599)I - F(612)) thought to support signal transmission to the heptahelical bundle. Thus domain 1 may contain a 189 residue N-terminal extracellular region followed successively by TM-1, a short intracellular loop, TM-2 and a 329 residue extracellular loop; rather than the proposed 620 residue VFT module based on crystallography of the N-terminal region of mGluR1. Since the topologies of the two proteins differ, the published CaSR VFT model is questionable. CaSR also contains multiple caveolin-binding motifs and cholesterol-binding (CRAC/CARC) domains, facilitating localization to plasma membrane lipid rafts. Ion sensing may involve combination of pore-lining regions from CaSR dimers and CaSR-bound caveolins to form ion channels capable of monitoring ionized Ca(2+) levels.

De novo membrane protein structure prediction

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

Membrane topology of hedgehog acyltransferase

Hedgehog acyltransferase (Hhat) is a multipass transmembrane enzyme that mediates the covalent attachment of the 16-carbon fatty acid palmitate to the N-terminal cysteine of Sonic Hedgehog (Shh). Palmitoylation of Shh by Hhat is critical for short and long range signaling. Knowledge of the topological organization of Hhat transmembrane helices would enhance our understanding of Hhat-mediated Shh palmitoylation. Bioinformatics analysis of transmembrane domains within human Hhat using 10 different algorithms resulted in highly consistent predictions in the C-terminal, but not in the N-terminal, region of Hhat. To empirically determine the topology of Hhat, we designed and exploited Hhat constructs containing either terminal or 12 different internal epitope tags. We used selective permeabilization coupled with immunofluorescence as well as a protease protection assay to demonstrate that Hhat contains 10 transmembrane domains and 2 re-entrant loops. The invariant His and highly conserved Asp residues within the membrane-bound O-acyltransferase (MBOAT) homology domain are segregated on opposite sides of the endoplasmic reticulum membrane. The localization of His-379 on the lumenal membrane surface is consistent with a role for this invariant residue in catalysis. Analysis of the activity and stability of the Hhat constructs revealed that the C-terminal MBOAT domain is especially sensitive to manipulation. Moreover, there was remarkable similarity in the overall topological organization of Hhat and ghrelin O-acyltransferase, another MBOAT family member. Knowledge of the topological organization of Hhat could serve as an important tool for further design of selective Hhat inhibitors.

The pore-lining regions in cytochrome c oxidases: A computational analysis of caveolin, cholesterol and transmembrane helix contributions to proton movement

Cytochrome c oxidase (CcO) is the terminal enzyme in the electron transfer chain. CcO catalyzes a four electron reduction of O2 to water at a catalytic site formed by high-spin heme (a3) and copper atoms (CuB). While it is recognized that proton movement is coupled to oxygen reduction, the proton channel(s) have not been well defined. Using computational methods developed to study protein topology, membrane channels and 3D packing arrangements within transmembrane (TM) helix arrays, we find that subunit-1 (COX-1), subunit-2 (COX-2) and subunit-3 (COX-3) contribute 139, 46 and 25 residues, respectively, to channel formation between the mitochondrial matrix and intermembrane space. Nine of 12 TM helices in COX-1, both helices in COX-2 and 5 of the 6 TM helices in COX-3 are pore-lining regions (possible channel formers). Heme a3 and the CuB sites (as well as the CuA center of COX-2) are located within the channel that includes TM-6, TM-7, TM-10 and TM-11 of COX-1 and are associated with multiple cholesterol and caveolin-binding (CB) motifs. Sequence analysis identifies five CB motifs within COX-1, two within COX-2 and four within COX-3; each caveolin containing a pore-lining helix C-terminal to a TM helix-turn-helix. Channel formation involves interaction between multiple pore-lining regions within protein subunits and/or dimers. PoreWalker analysis lends support to the D-channel model of proton translocation. Under physiological conditions, caveolins may introduce channel formers juxtaposed to those in COX-1, COX-2 and COX-3, which together with cholesterol may form channel(s) essential for proton translocation through the inner mitochondrial membrane.

Molecular evolution of a novel family of putative calcium transporters

The UPF0016 family is a group of uncharacterized membrane proteins, well conserved through evolution and defined by the presence of one or two copies of an E-Φ-G-D-(KR)-(ST) consensus motif. Our previous results have shown that two members of this family, the human TMEM165 and the budding yeast Gdt1p, are functionally related and might form a new group of cation/Ca2+ exchangers. Most members of the family are made of two homologous clusters of three transmembrane spans, separated by a central loop and assembled with an opposite orientation in the membrane. However, some bacterial members of the family have only one cluster of transmembrane domains. Among these 'single-domain membrane proteins' some cyanobacterial members were found as pairs of adjacent genes within the genome, but each gene was slightly different. We performed a bioinformatic analysis to propose the molecular evolution of the UPF0016 family and the emergence of the antiparallel topology. Our hypotheses were confirmed experimentally using functional complementation in yeast. This suggests an important and conserved function for UPF0016 proteins in a fundamental cellular process. We also show that members of the UPF0016 family share striking similarities, but no primary sequence homology, with members of the cation/Ca2+ exchangers (CaCA) superfamily. Such similarities could be an example of convergent evolution, supporting the previous hypothesis that members of the UPF0016 family are cation/Ca2+ exchangers.

Identification of key residues and regions important for porcupine-mediated Wnt acylation

Wnts comprise a family of lipid-modified, secreted signaling proteins that control embryogenesis, as well as tissue homeostasis in adults. Post-translational attachment of palmitoleate (C16:1) to a conserved Ser in Wnt proteins is catalyzed by Porcupine (Porcn), a member of the membrane bound O-acyltransferase (MBOAT) family, and is required for Wnt secretion and signaling. Moreover, genetic alterations in the PORCN gene lead to focal dermal hypoplasia, an X-linked developmental disorder. Despite its physiological importance, the biochemical mechanism governing Wnt acylation by Porcn is poorly understood. Here, we use a cell-based fatty acylation assay that is a direct readout of Porcn acyltransferase activity to perform structure-function analysis of highly conserved residues in Porcn and Wnt3a. In total, 16-point mutations in Porcn and 13 mutations in Wnt3a were generated and analyzed. We identified key residues within Porcn required for enzymatic activity, stability, and Wnt3a binding and mapped these active site residues to predicted transmembrane domain 9. Analysis of focal dermal hypoplasia-associated mutations in Porcn revealed that loss of enzymatic activity arises from altered stability. A consensus sequence within Wnt3a was identified (CXCHGXSXXCXXKXC) that contains residues that mediate Porcn binding, fatty acid transfer, and Wnt signaling. We also showed that Ser or Thr, but not Cys, can serve as a fatty acylation site in Wnt, establishing Porcn as an O-acyltransferase. This analysis sheds light into the mechanism by which Porcn transfers fatty acids to Wnt proteins and provides insight into the mechanisms of fatty acid transfer by MBOAT family members.

Genome segregation and packaging machinery in Acanthamoeba polyphaga mimivirus is reminiscent of bacterial apparatus

Genome packaging is a critical step in the virion assembly process. The putative ATP-driven genome packaging motor of Acanthamoeba polyphaga mimivirus (APMV) and other nucleocytoplasmic large DNA viruses (NCLDVs) is a distant ortholog of prokaryotic chromosome segregation motors, such as FtsK and HerA, rather than other viral packaging motors, such as large terminase. Intriguingly, APMV also encodes other components, i.e., three putative serine recombinases and a putative type II topoisomerase, all of which are essential for chromosome segregation in prokaryotes. Based on our analyses of these components and taking the limited available literature into account, here we propose for the first time a model for genome segregation and packaging in APMV that can possibly be extended to NCLDV subfamilies, except perhaps Poxviridae and Ascoviridae. This model might represent a unique variation of the prokaryotic system acquired and contrived by the large DNA viruses of eukaryotes. It is also consistent with previous observations that unicellular eukaryotes, such as amoebae, are melting pots for the advent of chimeric organisms with novel mechanisms.

IMPORTANCE

Extremely large viruses with DNA genomes infect a wide range of eukaryotes, from human beings to amoebae and from crocodiles to algae. These large DNA viruses, unlike their much smaller cousins, have the capability of making most of the protein components required for their multiplication. Once they infect the cell, these viruses set up viral replication centers, known as viral factories, to carry out their multiplication with very little help from the host. Our sequence analyses show that there is remarkable similarity between prokaryotes (bacteria and archaea) and large DNA viruses, such as mimivirus, vaccinia virus, and pandoravirus, in the way that they process their newly synthesized genetic material to make sure that only one copy of the complete genome is generated and is meticulously placed inside the newly synthesized viral particle. These findings have important evolutionary implications about the origin and evolution of large viruses.

Membrane protein orientation and refinement using a knowledge-based statistical potential

BACKGROUND

Recent increases in the number of deposited membrane protein crystal structures necessitate the use of automated computational tools to position them within the lipid bilayer. Identifying the correct orientation allows us to study the complex relationship between sequence, structure and the lipid environment, which is otherwise challenging to investigate using experimental techniques due to the difficulty in crystallising membrane proteins embedded within intact membranes.

RESULTS

We have developed a knowledge-based membrane potential, calculated by the statistical analysis of transmembrane protein structures, coupled with a combination of genetic and direct search algorithms, and demonstrate its use in positioning proteins in membranes, refinement of membrane protein models and in decoy discrimination.

CONCLUSIONS

Our method is able to quickly and accurately orientate both alpha-helical and beta-barrel membrane proteins within the lipid bilayer, showing closer agreement with experimentally determined values than existing approaches. We also demonstrate both consistent and significant refinement of membrane protein models and the effective discrimination between native and decoy structures. Source code is available under an open source license from http://bioinf.cs.ucl.ac.uk/downloads/memembed/.

Architectural organization of the metabolic regulatory enzyme ghrelin O-acyltransferase

Ghrelin O-acyltransferase (GOAT) is a polytopic integral membrane protein required for activation of ghrelin, a secreted metabolism-regulating peptide hormone. Although GOAT is a potential therapeutic target for the treatment of obesity and diabetes and plays a key role in other physiologic processes, little is known about its structure or mechanism. GOAT is a member of the membrane-bound O-acyltransferase (MBOAT) family, a group of polytopic integral membrane proteins involved in lipid-biosynthetic and lipid-signaling reactions from prokaryotes to humans. Here we use phylogeny and a variety of bioinformatic tools to predict the topology of GOAT. Using selective permeabilization indirect immunofluorescence microscopy in combination with glycosylation shift immunoblotting, we demonstrate that GOAT contains 11 transmembrane helices and one reentrant loop. Development of the V5Glyc tag, a novel, small, and sensitive dual topology reporter, facilitated these experiments. The MBOAT family invariant residue His-338 is in the ER lumen, consistent with other family members, but conserved Asn-307 is cytosolic, making it unlikely that both are involved in catalysis. Photocross-linking of synthetic ghrelin analogs and inhibitors demonstrates binding to the C-terminal region of GOAT, consistent with a role of His-338 in the active site. This knowledge of GOAT architecture is important for a deeper understanding of the mechanism of GOAT and other MBOATs and could ultimately advance the discovery of selective inhibitors for these enzymes.

Proton-dependent gating and proton uptake by Wzx support O-antigen-subunit antiport across the bacterial inner membrane

Wzx flippases are crucial for bacterial cell surface polysaccharide assembly as they transport undecaprenyl pyrophosphate-linked sugar repeat units from the cytoplasmic to the periplasmic leaflets of the inner membrane (IM) for final assembly. Our recently reported three-dimensional (3D) model structure of Wzx from Pseudomonas aeruginosa PAO1 (Wzx_Pa) displayed a cationic internal vestibule and functionally essential acidic amino acids within transmembrane segment bundles. Herein, we examined the intrinsic transport function of Wzx_Pa following its purification and reconstitution in phospholipid liposomes. Wzx_Pa was capable of mediating anion flux, consistent with its cationic interior. This flux was electrogenic and modified by extraliposomal pH. Mutation of the above-mentioned acidic residues (E61, D269, and D359) reduced proton (H⁺)-modified anion flux, showing the role of these amino acid side chains in H⁺-dependent transport. Wzx also mediated acidification of the proteoliposome interior in the presence of an outward anion gradient. These results indicate H⁺-dependent gating and H⁺ uptake by Wzx_Pa and allow for the first H⁺-dependent antiport mechanism to be proposed for lipid-linked oligosaccharide translocation across the bacterial IM.

Many bacterial cell surface polysaccharides that are important for survival and virulence are synthesized at the periplasmic leaflet of the inner membrane (IM) using precursors produced in the cytoplasm. Wzx flippases are responsible for translocation of lipid-linked sugar repeat units across the IM and had been previously suggested to simply facilitate passive substrate diffusion. Through our characterization of purified Wzx in a reconstitution system described herein, we have observed protein-dependent intrinsic transport producing a change in the electrical potential of the system, with H⁺ identified as the coupling ion. These results provide the first evidence for coupled (i.e., secondary active) transport by these proteins and, in conjunction with structural data, allow for an antiport mechanism to be proposed for the directed transport of lipid-linked sugar substrates across the IM. These findings bring our understanding of lipid-linked polysaccharide transporter proteins more in line with the efflux pumps to which they are evolutionarily related.

Scalable web services for the PSIPRED Protein Analysis Workbench

Here, we present the new UCL Bioinformatics Group's PSIPRED Protein Analysis Workbench. The Workbench unites all of our previously available analysis methods into a single web-based framework. The new web portal provides a greatly streamlined user interface with a number of new features to allow users to better explore their results. We offer a number of additional services to enable computationally scalable execution of our prediction methods; these include SOAP and XML-RPC web server access and new HADOOP packages. All software and services are available via the UCL Bioinformatics Group website at http://bioinf.cs.ucl.ac.uk/.

A simple method for predicting transmembrane proteins based on wavelet transform

The increasing protein sequences from the genome project require theoretical methods to predict transmembrane helical segments (TMHs). So far, several prediction methods have been reported, but there are some deficiencies in prediction accuracy and adaptability in these methods. In this paper, a method based on discrete wavelet transform (DWT) has been developed to predict the number and location of TMHs in membrane proteins. PDB coded as 1KQG is chosen as an example to describe the prediction process by this method. 80 proteins with known 3D structure from Mptopo database are chosen at random as data sets (including 325 TMHs) and 80 sequences are divided into 13 groups according to their function and type. TMHs prediction is carried out for each group of membrane protein sequences and obtain satisfactory result. To verify the feasibility of this method, 80 membrane protein sequences are treated as test sets, 308 TMHs can be predicted and the prediction accuracy is 96.3%. Compared with the main prediction results of seven popular prediction methods, the obtained results indicate that the proposed method in this paper has higher prediction accuracy.