Hinges, swivels and switches: the role of prolines in signalling via transmembrane alpha-helices

Two compound heterozygous variants in the CLN8 gene are responsible for neuronal cereidolipofuscinoses disorder in a child: a case report

Background

Neuronal Ceroid Lipofuscinosis (NCL) disorders, recognized as the primary cause of childhood dementia globally, constitute a spectrum of genetic abnormalities. CLN8, a subtype within NCL, is characterized by cognitive decline, motor impairment, and visual deterioration. This study focuses on an atypical case with congenital onset and a remarkably slow disease progression.

Methods

Whole-genome sequencing at 30× coverage was employed as part of a national genomics program to investigate the genetic underpinnings of rare diseases. This genomic approach aimed to challenge established classifications (vLINCL and EPMR) and explore the presence of a continuous phenotypic spectrum associated with CLN8.

Results

The whole-genome sequencing revealed two novel likely pathogenic mutations in the CLN8 gene on chromosome 8p23.3. These mutations were not previously associated with CLN8-related NCL. Contrary to established classifications (vLINCL and EPMR), our findings suggest a continuous phenotypic spectrum associated with CLN8. Pathological subcellular markers further validated the genomic insights.

Discussion

The identification of two previously undescribed likely pathogenic CLN8 gene mutations challenges traditional classifications and highlights a more nuanced phenotypic spectrum associated with CLN8. Our findings underscore the significance of genetic modifiers and interactions with unrelated genes in shaping variable phenotypic outcomes. The inclusion of pathological subcellular markers further strengthens the validity of our genomic insights. This research enhances our understanding of CLN8 disorders, emphasizing the need for comprehensive genomic analyses to elucidate the complexity of phenotypic presentations and guide tailored therapeutic strategies. The identification of new likely pathogenic mutations underscores the dynamic nature of CLN8-related NCL and the importance of individualized approaches to patient management.

Protein-Specific Signal Peptides for Mammalian Vector Engineering

Expression of recombinant proteins in mammalian cell factories relies on synthetic assemblies of genetic parts to optimally control flux through the product biosynthetic pathway. In comparison to other genetic part-types, there is a relative paucity of characterized signal peptide components, particularly for mammalian cell contexts. In this study, we describe a toolkit of signal peptide elements, created using bioinformatics-led and synthetic design approaches, that can be utilized to enhance production of biopharmaceutical proteins in Chinese hamster ovary cell factories. We demonstrate, for the first time in a mammalian cell context, that machine learning can be used to predict how discrete signal peptide elements will perform when utilized to drive endoplasmic reticulum (ER) translocation of specific single chain protein products. For more complex molecular formats, such as multichain monoclonal antibodies, we describe how a combination of in silico and targeted design rule-based in vitro testing can be employed to rapidly identify product-specific signal peptide solutions from minimal screening spaces. The utility of this technology is validated by deriving vector designs that increase product titers ≥1.8×, compared to standard industry systems, for a range of products, including a difficult-to-express monoclonal antibody. The availability of a vastly expanded toolbox of characterized signal peptide parts, combined with streamlined in silico/in vitro testing processes, will permit efficient expression vector re-design to maximize titers of both simple and complex protein products.

Using evolutionary data to make sense of macromolecules with a "face-lifted" ConSurf

The ConSurf web-sever for the analysis of proteins, RNA, and DNA provides a quick and accurate estimate of the per-site evolutionary rate among homologues. The analysis reveals functionally important regions, such as catalytic and ligand-binding sites, which often evolve slowly. Since the last report in 2016, ConSurf has been improved in multiple ways. It now has a user-friendly interface that makes it easier to perform the analysis and to visualize the results. Evolutionary rates are calculated based on a set of homologous sequences, collected using hidden Markov model-based search tools, recently embedded in the pipeline. Using these, and following the removal of redundancy, ConSurf assembles a representative set of effective homologues for protein and nucleic acid queries to enable informative analysis of the evolutionary patterns. The analysis is particularly insightful when the evolutionary rates are mapped on the macromolecule structure. In this respect, the availability of AlphaFold model structures of essentially all UniProt proteins makes ConSurf particularly relevant to the research community. The UniProt ID of a query protein with an available AlphaFold model can now be used to start a calculation. Another important improvement is the Python re-implementation of the entire computational pipeline, making it easier to maintain. This Python pipeline is now available for download as a standalone version. We demonstrate some of ConSurf's key capabilities by the analysis of caveolin-1, the main protein of membrane invaginations called caveolae.

The polar amino acid in the TatA transmembrane helix is not strictly necessary for protein function

The twin-arginine translocation (Tat) pathway utilizes the proton-motive force (pmf) to transport folded proteins across cytoplasmic membranes in bacteria and archaea, as well as across the thylakoid membrane in plants and the inner membrane in mitochondria. In most species, the minimal components required for Tat activity consist of three subunits, TatA, TatB, and TatC. Previous studies have shown that a polar amino acid is present at the N-terminus of the TatA transmembrane helix (TMH) across many different species. In order to systematically assess the functional importance of this polar amino acid in the TatA TMH in Escherichia coli, we examined a complete set of 19-amino-acid substitutions. Unexpectedly, although being preferred overall, our experiments suggest that the polar amino acid is not necessary for a functional TatA. Hydrophilicity and helix-stabilizing properties of this polar amino acid were found to be highly correlated with the Tat activity. Specifically, change in charge status of the amino acid side chain due to pH resulted in a shift in hydrophilicity, which was demonstrated to impact the Tat transport activity. Furthermore, we identified a four-residue motif at the N-terminus of the TatA TMH by sequence alignment. Using a biochemical approach, we found that the N-terminal motif was functionally significant, with evidence indicating a potential role in the preference for utilizing different pmf components. Taken together, these findings yield new insights into the functionality of TatA and its potential role in the Tat transport mechanism.

Investigation of novel de novo KCNC2 variants causing severe developmental and early-onset epileptic encephalopathy

Purpose The voltage-gated potassium channel K_v3.2, encoded by KCNC2, facilitates fast-spiking GABAergic interneurons to fire action potentials at high frequencies. It is pivotal to maintaining excitation/inhibition balance in mammalian brains. This study identified two novel de novo KCNC2 variants, p.Pro470Ser (P470S) and p.Phe382Leu (F382L), in patients with early onset developmental and epileptic encephalopathy (DEE). Methods To examine the molecular basis of DEE, we studied the functional characteristics of variant channels using patch-clamp techniques and computational modeling. Results Whole-cell patch clamp recordings from infected HEK293 cells revealed that channel activation and deactivation kinetics strongly decreased in both K_v3.2 P470S and F382L variant channels. This decrease also occurred in K_v3.2 p.Val471Leu (V471L) channels, known to be associated with DEE. In addition, K_v3.2 F382L and V471L variants exhibited a significant increase in channel conductance and a ∼20 mV negative shift in the threshold for voltage-dependent activation. Simulations of model GABAergic interneurons revealed that all variants decreased neuronal firing frequency. Thus, the variants' net loss-of-function effects disinhibited neural networks. Conclusion Our findings provide compelling evidence supporting the role of KCNC2 as a disease-causing gene in human neurodevelopmental delay and epilepsy.

Assembly of alpha-helical transmembrane pores through an intermediate state

Pore-forming alpha-helical proteins are well known for their dynamic assembly mechanism and it has been challenging to delineate the pore-forming structures in membranes. Previously, attempts have been made to elucidate their assembly mechanism and there is a large gap due to complex pathways by which these membrane-active pores impart their effect. Here we demonstrate a multi-step structural assembly pathway of alpha-helical peptide pores formed by a 37 amino acid synthetic peptide, pPorU, based on the natural porin from Corynebacterium urealyticum using single-channel electrical recordings. More specifically, we report detectable intermediate states during the membrane insertion and pore formation of pPorU. The fully assembled pore exhibited unusually large stable conductance, voltage-dependent gating, and functional blockage by cyclic sugars generally applicable to a range of transmembrane pores. Furthermore, we used rationally designed mutants to understand the role of specific amino acids in the assembly of these peptide pores. Mutant peptides that differ from wild-type peptides produced noisy and unstable intermediate states and low conductance pores, demonstrating sequence specificity in the pore-formation process supported by molecular dynamics simulations. We suggest that our study contributes to understanding the mechanism of action of naturally occurring alpha-helical pore-forming proteins and should be of broad interest to build peptide-based nanopore sensors.

Structural Characterization of Cytochrome c'β-Met from an Ammonia-Oxidizing Bacterium

The ammonia-oxidizing bacterium Nitrosomonas europaea expresses two cytochromes in the P460 superfamily that are predicted to be structurally similar. In one, cytochrome (cyt) P460, the substrate hydroxylamine (NH₂OH) is converted to nitric oxide (NO) and nitrous oxide (N₂O) requiring a unique heme-lysyl cross-link in the catalytic cofactor. In the second, cyt c'β-Met, the cross-link is absent, and the cytochrome instead binds H₂O₂ forming a ferryl species similar to compound II of peroxidases. Here, we report the 1.80 Å crystal structure of cyt c'β-Met─a well-expressed protein in N. europaea with a lysine to a methionine replacement at the cross-linking position. The structure of cyt c'β-Met is characterized by a large β-sheet typical of P460 members; however, several localized structural differences render cyt c'β-Met distinct. This includes a large lasso-like loop at the "top" of the cytochrome that is not observed in other structurally characterized members. Active site variation is also observed, especially in comparison to its closest homologue cyt c'β from the methane-oxidizing Methylococcus capsulatus Bath, which also lacks the cross-link. The phenylalanine "cap" which is presumed to control small ligand access to the distal heme iron is replaced with an arginine, reminiscent of the strictly conserved distal arginine in peroxidases and to the NH₂OH-oxidizing cytochromes P460. A critical proton-transferring glutamate residue required for NH₂OH oxidation is nevertheless missing in the active site. This in part explains the inability of cyt c'β-Met to oxidize NH₂OH. Our structure also rationalizes the absence of a methionyl cross-link, although the side chain's spatial position in the structure does not eliminate the possibility that it could form under certain conditions.

Deconstructing the transmembrane core of class A G protein-coupled receptors

Class A G protein-coupled receptors have evolved to recognize ligands ranging from small-molecule odorants to proteins. Although they are among the most diverse membrane receptors in eukaryotic organisms, they possess a highly conserved core within their seven-transmembrane helix framework. The conservation of the transmembrane core has led to the idea of a common mechanism by which ligand binding is coupled to the outward rotation of helix H6, the hallmark of an active receptor. Nevertheless, there is still no consensus on the mechanism of coupling or on the roles of specific residues within the core. Recent insights from crystallography and NMR spectroscopy provide a way to decompose the core into its essential structural and functional elements that shed new light on this important region.

Molecular and morphological approach to study the innexin gap junctions in Rhynchosciara americana

Gap junctions mediate communication between adjacent cells and are fundamental to the development and homeostasis in multicellular organisms. In invertebrates, gap junctions are formed by transmembrane proteins called innexins. Gap junctions allow the passage of small molecules through an intercellular channel, between a cell and another adjacent cell. The dipteran Rhynchosciara americana has contributed to studying the biology of invertebrates and the study of the interaction and regulation of genes during biological development. Therefore, this paper aimed to study the R. americana innexin-2 by molecular characterization, analysis of the expression profile and cellular localization. The molecular characterization results confirm that the message is from a gap junction protein and analysis of the expression and cellular localization profile shows that innexin-2 can participate in many physiological processes during the development of R. americana.

Structural determinants of the integrin transmembrane domain required for bidirectional signal transmission across the cell membrane

Studying the tight activity regulation of platelet-specific integrin α_IIbβ₃ is foundational and paramount to our understanding of integrin structure and activation. α_IIbβ₃ is essential for the aggregation and adhesion function of platelets in hemostasis and thrombosis. Structural and mutagenesis studies have previously revealed the critical role of α_IIbβ₃ transmembrane (TM) association in maintaining the inactive state. Gain-of-function TM mutations were identified and shown to destabilize the TM association leading to integrin activation. Studies using isolated TM peptides have suggested an altered membrane embedding of the β₃ TM α-helix coupled with α_IIbβ₃ activation. However, controversies remain as to whether and how the TM α-helices change their topologies in the context of full-length integrin in native cell membrane. In this study, we utilized proline scanning mutagenesis and cysteine scanning accessibility assays to analyze the structure and function correlation of the α_IIbβ₃ TM domain. Our identification of loss-of-function proline mutations in the TM domain suggests the requirement of a continuous TM α-helical structure in transmitting activation signals bidirectionally across the cell membrane, characterized by the inside-out activation for ligand binding and the outside-in signaling for cell spreading. Similar results were found for α_Lβ₂ and α₅β₁ TM domains, suggesting a generalizable mechanism. We also detected a topology change of β₃ TM α-helix within the cell membrane, but only under conditions of cell adhesion and the absence of α_IIb association. Our data demonstrate the importance of studying the structure and function of the integrin TM domain in the native cell membrane.

Motions around conserved helical weak spots facilitate GPCR activation

G protein‐coupled receptors (GPCRs) participate in most physiological processes and are important drug targets in many therapeutic areas. Recently, many GPCR X‐ray structures became available, facilitating detailed studies of their sequence‐structure‐mobility‐function relations. We show that the functional role of many conserved GPCR sequence motifs is to create weak spots in the transmembrane helices that provide the structural plasticity necessary for ligand binding and signaling. Different receptor families use different conserved sequence motifs to obtain similar helix irregularities that allow for the same motions upon GPCR activation. These conserved motions come together to facilitate the timely release of the conserved sodium ion to the cytosol. Most GPCR crystal structures could be determined only after stabilization of the transmembrane helices by mutations that remove weak spots. These mutations often lead to diminished binding of agonists, but not antagonists, which logically agrees with the fact that large helix rearrangements occur only upon agonist binding. Upon activation, six of the seven TM helices in GPCRs undergo helix motions and/or deformations facilitated by weak spots in these helices. The location of these weak spots is much more conserved than the sequence motifs that cause them. Knowledge about these weak spots helps understand the activation process of GPCRs and thus helps design medicines.

Adaptable and Multifunctional Ion-Conducting Aquaporins

Aquaporins function as water and neutral solute channels, signaling hubs, disease virulence factors, and metabolon components. We consider plant aquaporins that transport ions compared to some animal counterparts. These are candidates for important, as yet unidentified, cation and anion channels in plasma, tonoplast, and symbiotic membranes. For those individual isoforms that transport ions, water, and gases, the permeability spans 12 orders of magnitude. This requires tight regulation of selectivity via protein interactions and posttranslational modifications. A phosphorylation-dependent switch between ion and water permeation in AtPIP2;1 might be explained by coupling between the gates of the four monomer water channels and the central pore of the tetramer. We consider the potential for coupling between ion and water fluxes that could form the basis of an electroosmotic transducer. A grand challenge in understanding the roles of ion transporting aquaporins is their multifunctional modes that are dependent on location, stress, time, and development.

The role of prolines and glycine in the transmembrane domain of LAT

Linker for activation in T cells (LAT) is a critical regulator of T-cell development and function. It organises signalling events at the plasma membrane. However, the mechanism, which controls LAT localisation at the plasma membrane, is not fully understood. Here, we studied the impact of helix-breaking amino acids, two prolines and one glycine, in the transmembrane segment on localisation and function of LAT. Using in silico analysis, confocal and super-resolution imaging and flow cytometry, we demonstrate that central proline residue destabilises transmembrane helix by inducing a kink. The helical structure and dynamics are further regulated by glycine and another proline residue in the luminal part of LAT transmembrane domain. Replacement of these residues with aliphatic amino acids reduces LAT dependence on palmitoylation for sorting to the plasma membrane. However, surface expression of these mutants is not sufficient to recover function of nonpalmitoylated LAT in stimulated T cells. These data indicate that geometry and dynamics of LAT transmembrane segment regulate its localisation and function in immune cells.

GPR91, a critical signaling mechanism in modulating pathophysiologic processes in chronic illnesses

Succinate receptor GPR91 is one of G protein-coupled receptors (GPCRs), and is expressed in a variety of cell types and tissues. Succinate is its natural ligand, and its activation represents that an intrinsic metabolic intermediate exerts a regulatory role on many critical life processes involving pathophysiologic mechanisms, such as innate immunity, inflammation, tissue repair, and oncogenesis. With the illustration of 3-dimensional crystal structure of the receptor and discovery of its antagonists, it is possible to dissect the succinate-GPR91-G protein signaling pathways in different cell types under pathophysiological conditions. Deep understanding of the GPR91-ligand binding mode with various agonists and antagonists would aid in elucidating the molecular basis of a spectrum of chronic illnesses, such as hypertension, diabetes, and their renal and retina complications, metabolic-associated fatty liver diseases, such as nonalcoholic steatohepatitis and its fibrotic progression, inflammatory bowel diseases (Crohn's disease and ulcerative colitis), age-related macular degeneration, rheumatoid arthritis, and progressive behaviors of malignancies. With better delineation of critical regulatory role of the succinate-GPR91 axis in these illnesses, therapeutic intervention may be developed by specifically targeting this signaling pathway with small molecular antagonists or other strategies.

Gap Junction Coding Innexin in Lymnaea stagnalis: Sequence Analysis and Characterization in Tissues and the Central Nervous System

Connections between neurons called synapses are the key components underlying all nervous system functions of animals and humans. However, important genetic information on the formation and plasticity of one type, the electrical (gap junction-mediated) synapse, is understudied in many invertebrates. In the present study, we set forth to identify and characterize the gap junction-encoding gene innexin in the central nervous system (CNS) of the mollusk pond snail Lymnaea stagnalis. With PCR, 3′ and 5′ RACE, and BLAST searches, we identified eight innexin genes in the L. stagnalis genome, named Lst Inx1–Lst Inx8. Phylogenetic analysis revealed that the L. stagnalis innexin genes originated from a single copy in the common ancestor of molluskan species by multiple gene duplication events and have been maintained in L. stagnalis since they were generated. The paralogous innexin genes demonstrate distinct expression patterns among tissues. In addition, one paralog, Lst Inx1, exhibits heterogeneity in cells and ganglia, suggesting the occurrence of functional diversification after gene duplication. These results introduce possibilities to study an intriguing potential relationship between innexin paralog expression and cell-specific functional outputs such as heterogenic ability to form channels and exhibit synapse plasticity. The L. stagnalis CNS contains large neurons and functionally defined networks for behaviors; with the introduction of L. stagnalis in the gap junction gene field, we are providing novel opportunities to combine genetic research with direct investigations of functional outcomes at the cellular, synaptic, and behavioral levels.

Revealing the Positive Binding Cooperativity Mechanism between the Orthosteric and the Allosteric Antagonists of CCR2 by Metadynamics and Gaussian Accelerated Molecular Dynamics Simulations

CC chemokine receptor 2 (CCR2) and its endogenous CC chemokine ligands are associated with numerous inflammatory, neurodegenerative diseases, and cancer. CCR2 is becoming an attractive target in the treatment of autoimmune disease and neurodegenerative diseases. The orthosteric antagonist BMS-681 and allosteric antagonist CCR2-RA-[R] of CCR2 show positive binding cooperativity. We performed well-tempered metadynamics simulations and Gaussian accelerated MD simulations to reveal the influence of the orthosteric antagonist on the unbinding of allosteric antagonist of CCR2. We revealed different unbinding pathways of CCR2-RA-[R] in binary complex CCR2-VT5 and ternary complex CCR2-73R-VT5. The different unbinding pathways of CCR2-RA-[R] are due to the conformational dynamics of TM6. We obtained the significant conformational differences of the intracellular side of TM6 upon CCR2 binding to different ligands by GaMD simulation. The conformational dynamics of TM6 are consistent with the unbinding pathway analysis. GaMD simulations indicate that BMS-681 binding restricts the bend of intracellular side of TM6 by stabilizing the extracellular sides of TM6 and TM7. The charged residues Arg206^5.43 of TM5 and Glu291^7.39 of TM7 play key roles in stabling TM7 and TM6. TM6 and TM7 are crucial components in the orthosteric and allosteric binding sites. Our results illustrate the conformational details about the effect of the orthosteric antagonist on the allosteric antagonist of CCR2. The conformational dynamics of CCR2 upon binding to different ligands can provide a rational basis for development of allosteric ligands of CCR2.

GPR6 Structural Insights: Homology Model Construction and Docking Studies

GPR6 is an orphan G protein-coupled receptor that has been associated with the cannabinoid family because of its recognition of a sub-set of cannabinoid ligands. The high abundance of GPR6 in the central nervous system, along with high constitutive activity and a link to several neurodegenerative diseases make GPR6 a promising biological target. In fact, diverse research groups have demonstrated that GPR6 represents a possible target for the treatment of neurodegenerative disorders such as Parkinson's disease, Alzheimer's disease, and Huntington's disease. Several patents have claimed the use of a wide range of pyrazine derivatives as GPR6 inverse agonists for the treatment of Parkinson's disease symptoms and other dyskinesia syndromes. However, the full pharmacological importance of GPR6 has not yet been fully explored due to the lack of high potency, readily available ligands targeting GPR6. The long-term goal of the present study is to develop such ligands. In this paper, we describe our initial steps towards this goal. A human GPR6 homology model was constructed using a suite of computational techniques. This model permitted the identification of unique GPR6 structural features and the exploration of the GPR6 binding crevice. A subset of patented pyrazine analogs were docked in the resultant GPR6 inactive state model to validate the model, rationalize the structure-activity relationships from the reported patents and identify the key residues in the binding crevice for ligand recognition. We will take this structural knowledge into the next phase of GPR6 project, in which scaffold hopping will be used to design new GPR6 ligands.

The Structure and Function of Acylglycerophosphate Acyltransferase 4/ Lysophosphatidic Acid Acyltransferase Delta (AGPAT4/LPAATδ)

Lipid-modifying enzymes serve crucial roles in cellular processes such as signal transduction (producing lipid-derived second messengers), intracellular membrane transport (facilitating membrane remodeling needed for membrane fusion/fission), and protein clustering (organizing lipid domains as anchoring platforms). The lipid products crucial in these processes can derive from different metabolic pathways, thus it is essential to know the localization, substrate specificity, deriving products (and their function) of all lipid-modifying enzymes. Here we discuss an emerging family of these enzymes, the lysophosphatidic acid acyltransferases (LPAATs), also known as acylglycerophosphate acyltransferases (AGPATs), that produce phosphatidic acid (PA) having as substrates lysophosphatidic acid (LPA) and acyl-CoA. Eleven LPAAT/AGPAT enzymes have been identified in mice and humans based on sequence homologies, and their localization, specific substrates and functions explored. We focus on one member of the family, LPAATδ, a protein expressed mainly in brain and in muscle (though to a lesser extent in other tissues); while at the cellular level it is localized at the trans-Golgi network membranes and at the mitochondrial outer membranes. LPAATδ is a physiologically essential enzyme since mice knocked-out for Lpaatδ show severe dysfunctions including cognitive impairment, impaired force contractility and altered white adipose tissue. The LPAATδ physiological roles are related to the formation of its product PA. PA is a multifunctional lipid involved in cell signaling as well as in membrane remodeling. In particular, the LPAATδ-catalyzed conversion of LPA (inverted-cone-shaped lipid) to PA (cone-shaped lipid) is considered a mechanism of deformation of the bilayer that favors membrane fission. Indeed, LPAATδ is an essential component of the fission-inducing machinery driven by the protein BARS. In this process, a protein-tripartite complex (BARS/14-3-3γ/phosphoinositide kinase PI4KIIIβ) is recruited at the trans-Golgi network, at the sites where membrane fission is to occur; there, LPAATδ directly interacts with BARS and is activated by BARS. The resulting formation of PA is essential for membrane fission occurring at those spots. Also in mitochondria PA formation has been related to fusion/fission events. Since PA is formed by various enzymatic pathways in different cell compartments, the BARS-LPAATδ interaction indicates the relevance of lipid-modifying enzymes acting exactly where their products are needed (i.e., PA at the Golgi membranes).

Topology of TROL protein in thylakoid membranes of Arabidopsis thaliana

Thylakoid rhodanase-like protein (TROL) is a nuclear-encoded protein of thylakoid membranes required for tethering of ferredoxin:nicotinamide adenine dinucleotide phosphate (NADPH) oxydoreductase (FNR). It has been proposed that the dynamic interaction of TROL with flavoenzyme FNR, influenced by environmental light conditions, regulates the fate of photosynthetic electrons, directing them either to NADPH synthesis or to other acceptors, including reactive oxygen species detoxification pathways. Inside the chloroplasts, TROL has a dual localization: an inner membrane precursor form and a thylakoid membrane mature form, which has been confirmed by several large-scale chloroplast proteomics studies, as well as protein import experiments. Unlike the localization, the topology of TROL in the membranes, which is a prerequisite for further studies of its properties and function, has not been experimentally confirmed yet. Thermolysin was proven to be a valuable protease to probe the surface of chloroplasts and membranes in general. By treating the total chloroplast membranes using increasing protease concentration, sequential degradation of TROL was observed, indicating protected polypeptides of TROL and possible domain orientation. To further substantiate the obtained results, TROL-overexpressing Arabidopsis line (OX) and line in which the central rhodanase-like domain (RHO) has been partially deleted (ΔRHO), were used as well. While OX line showed the same degradation pattern of TROL as the wild-type, surprisingly, TROL from ΔRHO membranes was not detectable even at the lowest protease concentration applied, indicating the importance of this domain to the integrity of TROL. In conclusion, TROL is a polytopic protein with a stroma-exposed C-terminal FNR-binding region, and the thylakoid lumen-located RHO domain.

Proline residues in scavenger receptor-BI's C-terminal region support efficient cholesterol transport

High-density lipoproteins (HDLs) facilitate reverse cholesterol transport, a process in which HDL removes cholesterol from circulation and carries it to the liver for biliary excretion. Reverse cholesterol transport is also facilitated by HDL's high-affinity receptor, scavenger receptor-BI (SR-BI), by mechanisms that are not fully understood. To improve our understanding of SR-BI function, we previously solved the NMR (nuclear magnetic resonance) structure of a peptide encompassing amino acids 405–475 of SR-BI. This segment of SR-BI, that includes the functionally critical C-terminal transmembrane domain and part of the extracellular domain, also contains four conserved proline (Pro) residues. We hypothesized that these proline residues support SR-BI in a conformation that allows for efficient cholesterol transport. To test this, we generated individual Pro-to-alanine mutations in full-length SR-BI and transiently expressed the mutant receptors in COS-7 cells to measure the effects on SR-BI-mediated cholesterol transport functions. Our findings reveal that HDL cell association and uptake of HDL-cholesteryl esters are impaired by mutation of Pro-412, Pro-438, or the transmembrane proline kink residue (Pro-459). In addition, SR-BI-mediated cholesterol efflux and membrane cholesterol distribution are impaired by mutation of Pro-412 or Pro-438, indicating that these residues are essential for a fully functional SR-BI receptor. Furthermore, we demonstrate that Pro-408 is necessary for proper SR-BI expression, but mutation of Pro-408 does not cause SR-BI to become misfolded or rapidly degraded by the proteasome or the lysosome. We conclude that key proline residues play an important role in SR-BI function by allowing for the efficient transport of cholesterol between cells and HDL.

Reverse Engineering of a Thermosensing Regulator Switch

To address the mechanism of thermosensing and its implications for molecular engineering, we previously deconstructed the functional components of the bacterial thermosensor DesK, a histidine kinase with a five-span transmembrane domain that detects temperature changes. The system was first simplified by building a sensor that consists of a single chimerical transmembrane segment that retained full sensing capacity. Genetic and biophysical analysis of this minimal sensor enabled the identification of three modular components named determinants of thermodetection (DOTs). Here we combine and tune the DOTs to determine their contribution to activity. A transmembrane zipper represents the master DOT that drives a reversible and activating dimerization through the formation of hydrogen bonds. Our findings provide the mechanism and insights to construct a synthetic transmembrane helix based on a poly-valine scaffold that harbors the DOTs and regulates the activity. The construct constitutes a modular switch that may be exploited in biotechnology and genetic circuitry.

Helical Structure of Recombinant Melittin

Melittin is an extensively studied, 26-residue toxic peptide from honey bee venom. Because of its versatility in adopting a variety of secondary (helix or coil) and quaternary (monomer or tetramer) structures in various environments, melittin has been the focus of numerous investigations as a model peptide in protein folding studies as well as in studies involving binding to proteins, lipids, and polysaccharides. A significant body of evidence supports the view that melittin binds to these macromolecules in a predominantly helical conformation, but detailed structural knowledge of this conformation is lacking. In this report, we present nuclear magnetic resonance (NMR)-based structural insights into the helix formation of recombinant melittin in the presence of trifluoroethanol (TFE): a known secondary structure inducer in peptides. These studies were performed at neutral pH, with micromolar amounts of the peptide. Using nuclear Overhauser effect (NOE)-derived distance restraints from three-dimensional NMR spectra, we determined the atomic resolution solution NMR structure of recombinant melittin bearing a TFE-stabilized helix. To circumvent the complications with structure determination of small peptides with high conformational flexibility, we developed a workflow for enhancing proton NOEs by increasing the viscosity of the medium. In the TFE-containing medium, recombinant monomeric melittin forms a long, continuous helical structure, which consists of the N- and C-terminal α-helices and the noncanonical 3₁₀-helix in the middle. The noncanonical 3₁₀-helix is missing in the previously solved X-ray structure of tetrameric melittin and the NMR structure of melittin in methanol. Melittin's structure in TFE-containing medium provides insights into melittin's conformational transitions, which are relevant to the peptide's interactions with its biological targets.

Main-chain mutagenesis reveals intrahelical coupling in an ion channel voltage-sensor

Membrane proteins are universal signal decoders. The helical transmembrane segments of these proteins play central roles in sensory transduction, yet the mechanistic contributions of secondary structure remain unresolved. To investigate the role of main-chain hydrogen bonding on transmembrane function, we encoded amide-to-ester substitutions at sites throughout the S4 voltage-sensing segment of Shaker potassium channels, a region that undergoes rapid, voltage-driven movement during channel gating. Functional measurements of ester-harboring channels highlight a transitional region between α-helical and 3₁₀ segments where hydrogen bond removal is particularly disruptive to voltage-gating. Simulations of an active voltage sensor reveal that this region features a dynamic hydrogen bonding pattern and that its helical structure is reliant upon amide support. Overall, the data highlight the specialized role of main-chain chemistry in the mechanism of voltage-sensing; other catalytic transmembrane segments may enlist similar strategies in signal transduction mechanisms.

The helical transmembrane segments of membrane proteins play central roles in sensory transduction but the mechanistic basis for their function remains unresolved. Here the authors identify regions in the S4 voltage-sensing segment of Shaker potassium channels where local helical structure is reliant upon backbone amide support.

Influenza hemagglutinin membrane anchor

Viruses with membranes fuse them with cellular membranes, to transfer their genomes into cells at the beginning of infection. For Influenza virus, the membrane glycoprotein involved in fusion is the hemagglutinin (HA), the 3D structure of which is known from X-ray crystallographic studies. The soluble ectodomain fragments used in these studies lacked the "membrane anchor" portion of the molecule. Since this region has a role in membrane fusion, we have determined its structure by analyzing the intact, full-length molecule in a detergent micelle, using cryo-EM. We have also compared the structures of full-length HA-detergent micelles with full-length HA-Fab complex detergent micelles, to describe an infectivity-neutralizing monoclonal Fab that binds near the ectodomain membrane anchor junction. We determine a high-resolution HA structure which compares favorably in detail with the structure of the ectodomain seen by X-ray crystallography; we detect, clearly, all five carbohydrate side chains of HA; and we find that the ectodomain is joined to the membrane anchor by flexible, eight-residue-long, linkers. The linkers extend into the detergent micelle to join a central triple-helical structure that is a major component of the membrane anchor.

A statistical model for helices with applications

Motivated by a cutting edge problem related to the shape of α -helices in proteins, we formulate a parametric statistical model, which incorporates the cylindrical nature of the helix. Our focus is to detect a "kink," which is a drastic change in the axial direction of the helix. We propose a statistical model for the straight α -helix and derive the maximum likelihood estimation procedure. The cylinder is an accepted geometric model for α -helices, but our statistical formulation, for the first time, quantifies the uncertainty in atom positions around the cylinder. We propose a change point technique "Kink-Detector" to detect a kink location along the helix. Unlike classical change point problems, the change in direction of a helix depends on a simultaneous shift of multiple data points rather than a single data point, and is less straightforward. Our biological building block is crowdsourced data on straight and kinked helices; which has set a gold standard. We use this data to identify salient features to construct Kink-detector, test its performance and gain some insights. We find the performance of Kink-detector comparable to its computational competitor called "Kink-Finder." We highlight that identification of kinks by visual assessment can have limitations and Kink-detector may help in such cases. Further, an analysis of crowdsourced curved α -helices finds that Kink-detector is also effective in detecting moderate changes in axial directions.

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.

Structural basis for the interaction of the beta-secretase with copper

The β-secretase (BACE1) features a unique sulfur rich motif (M₄₆₂xxxC₄₆₆xxxM₄₇₀xxxC₄₇₄xxxC₄₇₈) in its transmembrane helix (BACE1-TM) which is characteristic for proteins involved in copper ion storage and transport. While this motif has been shown to promote BACE1-TM trimerization and binding of copper ions in vitro, the structural basis for the interaction of copper ions with the BACE1-TM is still not well understood. Using molecular dynamics (MD) simulations, we show that membrane embedded BACE1-TMs adopt a flexible trimeric structure that binds and conducts copper ions through variable coordination. In coarse-grained (CG) MD simulations, the spontaneous assembly of BACE1-TMs trimers results in a right-handed helix packing arrangement. In subsequent atomistic MD simulations the sulfur rich motif defines characteristic copper ion coordination sites along a constricted partially solvated axial pore. Sliding and tilting of BACE1-TMs along smooth A₄₅₉xxxA_463/464xxA₄₆₇ surfaces, facilitated by a central P472 induced kink, enables copper ions to alternate between different coordination sites, including the prominent C466 and M470. We shed light into the structural arrangement of BACE1-TM trimers and propose a mechanism for copper ion conduction that might also apply to other proteins involved in metal ion transport.

APLP1 is endoproteolytically cleaved by γ-secretase without previous ectodomain shedding

Regulated intramembrane proteolysis of the amyloid precursor protein (APP) and its homologs, the APP like proteins APLP1 and APLP2, is typically a two-step process, which is initiated by ectodomain-shedding of the substrates by α- or β-secretases. Growing evidence, however, indicates that the cleavage process for APLP1 is different than for APP. Here, we describe that full-length APLP1, but not APP or APLP2, is uniquely cleaved by γ-secretase without previous ectodomain shedding. The new fragment, termed sAPLP1γ, was exclusively associated with APLP1, not APP, APLP2. We provide an exact molecular analysis showing that sAPLP1γ was uniquely generated by γ-secretase from full-length APLP1. Mass spectrometry analysis showed that the sAPLP1γ fragment and the longest Aβ-like peptide share the C-terminus. This novel mechanism of γ-secretase action is consistent with an ϵ-cut based upon the nature of the reaction in APP. We further demonstrate that the APLP1 transmembrane sequence is the critical determinant for γ-shedding and release of full-length APLP1. Moreover, the APLP1 TMS is sufficient to convert larger type-I membrane proteins like APP into direct γ-secretase substrates. Taken together, the direct cleavage of APLP1 is a novel feature of the γ-secretase prompting a re-thinking of γ-secretase activity modulation as a therapeutic strategy for Alzheimer disease.

Transmembrane Signal Transduction in Two-Component Systems: Piston, Scissoring, or Helical Rotation?

Allosteric and transmembrane (TM) signaling are among the major questions of structural biology. Here, we review and discuss signal transduction in four-helical TM bundles, focusing on histidine kinases and chemoreceptors found in two-component systems. Previously, piston, scissors, and helical rotation have been proposed as the mechanisms of TM signaling. We discuss theoretically possible conformational changes and examine the available experimental data, including the recent crystallographic structures of nitrate/nitrite sensor histidine kinase NarQ and phototaxis system NpSRII:NpHtrII. We show that TM helices can flex at multiple points and argue that the various conformational changes are not mutually exclusive, and often are observed concomitantly, throughout the TM domain or in its part. The piston and scissoring motions are the most prominent motions in the structures, but more research is needed for definitive conclusions.

Conformational Change in the Ligand-Binding Pocket via a KISS1R Mutation (P147L) Leads to Isolated Gonadotropin-Releasing Hormone Deficiency

Context:

Kisspeptin receptor (KISS1R) is expressed in hypothalamic gonadotropin-releasing hormone neurons and responsible for pubertal onset and reproductive functions. KISS1R mutations remain a rare cause of congenital hypogonadotropic hypogonadism (CHH).

Objective:

The aim of this study was to identify the genetic cause of CHH in a patient and to functionally characterize a KISS1R mutation.

Design:

The patient was a 47-year-old Japanese man whose parents were first cousins. He lacked secondary sexual characteristics owing to normosmic CHH. Exon segments for the KISS1R gene in this patient were screened for mutations. Functional analyses were performed using HEK293 cells expressing KISS1R mutants. Molecular dynamics simulations were performed to compare the ligand-KISS1R mutant complex with those of wild-type KISS1R variants.

Results:

A homozygous mutation (c.440C>T, p.P147L) in KISS1R was identified. The P147L mutation did not affect either receptor expression level or subcellular localization in the recombinant expression system. Intracellular calcium measurements and cellular dielectric spectroscopy demonstrated that the P147L mutation impaired receptor function to an extent more severe than that of a previously reported L148S mutation. A receptor-ligand binding assay showed the P147L mutation causes a substantial loss of ligand-binding affinity. Molecular dynamics simulations revealed the P147L mutation decreases the contact surface area of the ligand-receptor complex in an expanded ligand-binding pocket.

Conclusion:

We identified a loss-of-function mutation in KISS1R associated with CHH. Our results demonstrated that the P147L mutation causes a severe phenotype and functional impairment resulting from the loss of ligand-binding affinity due to an expanded ligand-binding pocket.

Allosteric Fine-Tuning of the Binding Pocket Dynamics in the ITK SH2 Domain by a Distal Molecular Switch: An Atomistic Perspective

Although the regulation of function of proteins by allosteric interactions has been identified in many subcellular processes, molecular switches are also known to induce long-range conformational changes in proteins. A less well understood molecular switch involving cis-trans isomerization of a peptidyl-prolyl bond could induce a conformational change directly to the backbone that is propagated to other parts of the protein. However, these switches are elusive and hard to identify because they are intrinsic to biomolecules that are inherently dynamic. Here, we explore the conformational dynamics and free energy landscape of the SH2 domain of interleukin-2-inducible T-cell or tyrosine kinase (ITK) to fully understand the conformational coupling between the distal cis-trans molecular switch and its binding pocket of the phosphotyrosine motif. We use multiple microsecond-long all-atom molecular dynamics simulations in explicit water for over a total of 60 μs. We show that cis-trans isomerization of the Asn286-Pro287 peptidyl-prolyl bond is directly coupled to the dynamics of the binding pocket of the phosphotyrosine motif, in agreement with previous NMR experiments. Unlike the cis state that is localized and less dynamic in a single free energy basin, the trans state samples two distinct conformations of the binding pocket-one that recognizes the phosphotyrosine motif and the other that is somewhat similar to that of the cis state. The results provide an atomic-level description of a less well understood allosteric regulation by a peptidyl-prolyl cis-trans molecular switch that could aid in the understanding of normal and aberrant subcellular processes and the identification of these elusive molecular switches in other proteins.

Hydrogen/Deuterium Exchange Mass Spectrometry of Human Green Opsin Reveals a Conserved Pro-Pro Motif in Extracellular Loop 2 of Monostable Visual G Protein-Coupled Receptors

Opsins comprise the protein component of light sensitive G protein-coupled receptors (GPCRs) in the retina of the eye that are responsible for the transduction of light into a biochemical signal. Here, we used hydrogen/deuterium (H/D) exchange coupled with mass spectrometry to map conformational changes in green cone opsin upon light activation. We then compared these findings with those reported for rhodopsin. The extent of H/D exchange in green cone opsin was greater than in rhodopsin in the dark and bleached states, suggesting a higher structural heterogeneity for green cone opsin. Further analysis revealed that green cone opsin exists as a dimer in both dark (inactive) and bleached (active) states, and that the predicted glycosylation sites at N32 and N34 are indeed glycosylated. Comparison of deuterium uptake between inactive and active states of green cone opsin also disclosed a reduced solvent accessibility of the extracellular N-terminal region and an increased accessibility of the chromophore binding site. Increased H/D exchange at the extracellular side of transmembrane helix four (TM4) combined with an analysis of sequence alignments revealed a conserved Pro-Pro motif in extracellular loop 2 (EL2) of monostable visual GPCRs. These data present new insights into the locus of chromophore release at the extracellular side of TM4 and TM5 and provide a foundation for future functional evaluation.

Impact of geometry changes in the channel pore by the gating movements on the channel's conductance

Kv 1.2 are voltage-dependent potassium channels of great biological importance. Despite the existence of many reports considering structure - function relations of the Kv 1.2 channel's quantitative domains, some details of the voltage gating remain ambiguous, or even unknown. One of the examples is the range of the S4-S6 domains motions involved in channel activation and gating. Another important question is to what extent the channel geometry influences the observable channel conductance at different voltages, and what mechanism stands behind. Does the narrowing of the pore reduce the conductance by ohmic resistance growth? The answer is surprisingly negative. But it can be explained in an alternative way by considering the fluctuations. To address these problems, we formulate geometric models that mimic the generic features of voltage sensor movement and trigger the movement of the other domains involved in gating. We carry out a complete simulation of S4-S6 domains translations and tilts. The obtained pore profiles allow to estimate the (ohmic) conductance dependency on the voltage. From a family of analysed models, we choose the one most concurring with the experimental data. The results allow to suggest the most probable scenario of S4-S6 domains movement during channel activation by membrane depolarization.

Effects of isosteric substitutions on the conformational preference and cis–trans isomerization of proline-containing peptides

Isosteric substitutions of the peptide CO group by CS and CSe groups increased thetranspopulation and rotational barrier to the prolylcis–transisomerization of proline-containing peptides.

The Transmembrane Domain of a Bicomponent ABC Transporter Exhibits Channel-Forming Activity

Pseudomonas aeruginosa is an opportunistic pathogen that expresses two unique forms of lipopolysaccharides (LPSs) on its bacterial surface, the A- and B-bands. The A-band polysaccharides (A-band PSs) are thought to be exported into the periplasm via a bicomponent ATP-binding cassette (ABC) transporter located within the inner membrane. This ABC protein complex consists of the transmembrane (TMD) Wzm and nucleotide-binding (NBD) Wzt domain proteins. Here, we were able to probe ∼1.36 nS-average conductance openings of the Wzm-based protein complex when reconstituted into a lipid membrane buffered by a 200 mM KCl solution, demonstrating the large-conductance, channel-forming ability of the TMDs. In agreement with this finding, transmission electron microscopy (TEM) imaging revealed the ring-shaped structure of the transmembrane Wzm protein complex. As hypothesized, using liposomes, we demonstrated that Wzm interacts with Wzt. Further, the Wzt polypeptide indeed hydrolyzed ATP but exhibited a ∼75% reduction in the ATPase activity when its Walker A domain was deleted. The distribution and average unitary conductance of the TMD Wzm protein complex were altered by the presence of the NBD Wzt protein, confirming the regulatory role of the latter polypeptide. To our knowledge, the large-conductance, channel-like activity of the Wzm protein complex, although often hypothesized, has not previously been demonstrated. These results constitute a platform for future structural, biophysical, and functional explorations of this bicomponent ABC transporter.

Trans-membrane Signaling in Photosynthetic State Transitions: REDOX- AND STRUCTURE-DEPENDENT INTERACTION IN VITRO BETWEEN STT7 KINASE AND THE CYTOCHROME b6f COMPLEX

Trans-membrane signaling involving a serine/threonine kinase (Stt7 in Chlamydomonas reinhardtii) directs light energy distribution between the two photosystems of oxygenic photosynthesis. Oxidation of plastoquinol mediated by the cytochrome b₆f complex on the electrochemically positive side of the thylakoid membrane activates the kinase domain of Stt7 on the trans (negative) side, leading to phosphorylation and redistribution ("state transition") of the light-harvesting chlorophyll proteins between the two photosystems. The molecular description of the Stt7 kinase and its interaction with the cytochrome b₆f complex are unknown or unclear. In this study, Stt7 kinase has been cloned, expressed, and purified in a heterologous host. Stt7 kinase is shown to be active in vitro in the presence of reductant and purified as a tetramer, as determined by analytical ultracentrifugation, electron microscopy, and electrospray ionization mass spectrometry, with a molecular weight of 332 kDa, consisting of an 83.41-kDa monomer. Far-UV circular dichroism spectra show Stt7 to be mostly α-helical and document a physical interaction with the b₆f complex through increased thermal stability of Stt7 secondary structure. The activity of wild-type Stt7 and its Cys-Ser mutant at positions 68 and 73 in the presence of a reductant suggest that the enzyme does not require a disulfide bridge for its activity as suggested elsewhere. Kinase activation in vivo could result from direct interaction between Stt7 and the b₆f complex or long-range reduction of Stt7 by superoxide, known to be generated in the b₆f complex by quinol oxidation.

Helix perturbations in membrane proteins assist in inter-helical interactions and optimal helix positioning in the bilayer

Transmembrane (TM) helices in integral membrane proteins are primarily α-helical in structure. Here we analyze 1134 TM helices in 90 high resolution membrane proteins and find that apart from the widely prevalent α-helices, TM regions also contain stretches of 3₁₀ (3 to 8 residues) and π-helices (5 to 19 residues) with distinct sequence signatures. The various helix perturbations in TM regions comprise of helices with kinked geometry, as well as those with an interspersed 3₁₀/π-helical fragment and show high occurrence in a few membrane proteins. Proline is frequently present at sites of these perturbations, but it is neither a necessary nor a sufficient requirement. Helix perturbations are also conserved within a family of membrane proteins despite low sequence identity in the perturbed region. Furthermore, a perturbation influences the geometry of the TM helix, mediates inter-helical interactions within and across protein chains and avoids hydrophobic mismatch of the helix termini with the bilayer. An analysis of π-helices in the TM regions of the heme copper oxidase superfamily shows that interspersed π-helices can vary in length from 6 to 19 amino acids or be entirely absent, depending upon the protein function. The results presented here would be helpful for prediction of 3₁₀ and π-helices in TM regions and can assist the computational design of membrane proteins.

Free backbone carbonyls mediate rhodopsin activation

Conserved prolines in the transmembrane helices of G-protein-coupled receptors (GPCRs) are often considered to function as hinges that divide the helix into two segments capable of independent motion. Depending on their potential to hydrogen-bond, the free C=O groups associated with these prolines can facilitate conformational flexibility, conformational switching or stabilization of the receptor structure. To address the role of conserved prolines in family A GPCRs through solid-state NMR spectroscopy, we focus on bovine rhodopsin, a GPCR in the visual receptor subfamily. The free backbone C=O groups on helices H5 and H7 stabilize the inactive rhodopsin structure through hydrogen-bonds to residues on adjacent helices. In response to light-induced isomerization of the retinal chromophore, hydrogen-bonding interactions involving these C=O groups are released, thus facilitating repacking of H5 and H7 onto the transmembrane core of the receptor. These results provide insights into the multiple structural and functional roles of prolines in membrane proteins.

In silico Exploration of the Conformational Universe of GPCRs

The structural plasticity of G protein coupled receptors (GPCRs) leads to a conformational universe going from inactive to active receptor states with several intermediate states. Many of them have not been captured yet and their role for GPCR activation is not well understood. The study of this conformational space and the transition dynamics between different receptor populations is a major challenge in molecular biophysics. The rational design of effector molecules that target such receptor populations allows fine-tuning receptor signalling with higher specificity to produce drugs with safer therapeutic profiles. In this minireview, we outline highly conserved receptor regions which are considered determinant for the establishment of distinct receptor states. We then discuss in-silico approaches such as dimensionality reduction methods and Markov State Models to explore the GPCR conformational universe and exploit the obtained conformations through structure-based drug design.

Conformational changes of an ion-channel during gating and emerging electrophysiologic properties: Application of a computational approach to cardiac Kv7.1

Ion channels are the "building blocks" of the excitation process in excitable tissues. Despite advances in determining their molecular structure, understanding the relationship between channel protein structure and electrical excitation remains a challenge. The Kv7.1 potassium channel is an important determinant of the cardiac action potential and its adaptation to rate changes. It is subject to beta adrenergic regulation, and many mutations in the channel protein are associated with the arrhythmic long QT syndrome. In this theoretical study, we use a novel computational approach to simulate the conformational changes that Kv7.1 undergoes during activation gating and compute the resulting electrophysiologic function in terms of single-channel and macroscopic currents. We generated all possible conformations of the S4-S5 linker that couples the S3-S4 complex (voltage sensor domain, VSD) to the pore, and all associated conformations of VSD and the pore (S6). Analysis of these conformations revealed that VSD-to-pore mechanical coupling during activation gating involves outward translation of the voltage sensor, accompanied by a translation away from the pore and clockwise twist. These motions cause pore opening by moving the S4-S5 linker upward and away from the pore, providing space for the S6 tails to move away from each other. Single channel records, computed from the simulated motion trajectories during gating, have stochastic properties similar to experimentally recorded traces. Macroscopic current through an ensemble of channels displays two key properties of Kv7.1: an initial delay of activation and fast inactivation. The simulations suggest a molecular mechanism for fast inactivation; a large twist of the VSD following its outward translation results in movement of the base of the S4-S5 linker toward the pore, eliminating open pore conformations to cause inactivation.

Characterizing the Murine Leukemia Virus Envelope Glycoprotein Membrane-Spanning Domain for Its Roles in Interface Alignment and Fusogenicity

The membrane-proximal region of murine leukemia virus envelope (Env) is a critical modulator of its functionality. We have previously shown that the insertion of one amino acid (+1 leucine) within the membrane-spanning domain (MSD) abolished protein functionality in infectivity assays. However, functionality could be restored to this +1 leucine mutant by either inserting two additional amino acids (+3 leucine) or by deleting the cytoplasmic tail domain (CTD) in the +1 leucine background. We inferred that the ectodomain and CTD have protein interfaces that have to be in alignment for Env to be functional. Here, we made single residue deletions to the Env mutant with the +1 leucine insertion to restore the interface alignment (gain of functionality) and therefore define the boundaries of the two interfaces. We identified the glycine-proline pairs near the N terminus (positions 147 and 148) and the C terminus (positions 159 and 160) of the MSD as being the boundaries of the two interfaces. Deletions between these pairs restored function, but deletions outside of them did not. In addition, the vast majority of the single residue deletions regained function if the CTD was deleted. The exceptions were four hydroxyl-containing amino acid residues (T139, T140, S143, and T144) that reside in the ectodomain interface and the proline at position 148, which were all indispensable for functionality. We hypothesize that the hydroxyl-containing residues at positions T139 and S143 could be a driving force for stabilizing the ectodomain interface through formation of a hydrogen-bonding network.

IMPORTANCE

The membrane-proximal external region (MPER) and membrane-spanning domains (MSDs) of viral glycoproteins have been shown to be critical for regulating glycoprotein fusogenicity. However, the roles of these two domains are poorly understood. We report here that point deletions and insertions within the MPER or MSD result in functionally inactive proteins. However, when the C-terminal tail domain (CTD) is deleted, the majority of the proteins remain functional. The only residues that were found to be critical for function regardless of the CTD were four hydroxyl-containing amino acids located at the C terminus of the MPER (T139 and T140) and at the N terminus of the MSD (S143 and T144) and a proline near the beginning of the MSD (P148). We demonstrate that hydrogen-bonding at positions T139 and S143 is critical for protein function. Our findings provide novel insights into the role of the MPER in regulating fusogenic activity of viral glycoproteins.

A novel innexin2 forming membrane hemichannel exhibits immune responses and cell apoptosis in Scylla paramamosain

Innexins are a class of transmembrane proteins that are important for embryonic development, morphogenesis and electrical synapse formation. In the present study, a novel innexin2 gene from Scylla paramamosain was named Sp-inx2 and characterized. The complete cDNA and genomic DNA sequences of Sp-inx2 were revealed. Sp-inx2 mRNA transcripts were distributed in various tissues of S. paramamosain and were most abundant in the hemocytes. The Sp-inx2 was significantly upregulated in hemocyte, gill and hepatopancreas tissues with the challenge of either Vibrio alginolyticus, Vibrio parahaemolyticus or lipopolysaccharides (LPSs) when analyzed at 3 and 6 h using quantitative real-time PCR, suggesting that it could activate an immune response against the challenge of LPSs or Vibrio species. Using the chemical inhibitors carbenoxolone and probenecid, the absorption of the fluorescent dye Lucifer yellow decreased in the primary cultured hemocytes of crabs, thus confirming that hemichannels composed of Sp-inx2 existed in the crab hemocytes. With LPS stimulation, the level of mRNA transcripts and protein expression of Sp-inx2 in the same cultured hemocytes gradually increased from 6 to 48 h, while the activity of hemichannels was down-regulated at 6 and 12 h, demonstrating that LPSs could modulate the absorption activity of hemichannels in addition to its upregulation of Sp-inx2 gene expression. Furthermore, the dye uptake rate in HeLa cells in which Sp-inx2 was ectopically expressed increased dramatically but the increase was significantly down-regulated with the addition of 50 μg mL(-1) LPS, suggesting that the LPS stimulation could effectively reduce the activity of hemichannels. Interestingly, with the ectopic expression of Sp-inx2 in HeLa and EPC cells, apoptosis spontaneously occurred in both cultured cell lines when detected using TUNEL assay. In summary, a new Sp-inx2 gene was first characterized in a marine animal S. paramamosain and it had a function associated with immune response and cell apoptosis.

Mutational mapping of the transmembrane binding site of the G-protein coupled receptor TGR5 and binding mode prediction of TGR5 agonists

TGR5 (Gpbar-1, M-Bar) is a class A G-protein coupled bile acid-sensing receptor predominately expressed in brain, liver and gastrointestinal tract, and a promising drug target for the treatment of metabolic disorders. Due to the lack of a crystal structure of TGR5, the development of TGR5 agonists has been guided by ligand-based approaches so far. Three binding mode models of bile acid derivatives have been presented recently. However, they differ from one another in terms of overall orientation or with respect to the location and interactions of the cholane scaffold, or cannot explain all results from mutagenesis experiments. Here, we present an extended binding mode model based on an iterative and integrated computational and biological approach. An alignment of 68 TGR5 agonists based on this binding mode leads to a significant and good structure-based 3D QSAR model, which constitutes the most comprehensive structure-based 3D-QSAR study of TGR5 agonists undertaken so far and suggests that the binding mode model is a close representation of the "true" binding mode. The binding mode model is further substantiated in that effects predicted for eight mutations in the binding site agree with experimental analyses on the impact of these TGR5 variants on receptor activity. In the binding mode, the hydrophobic cholane scaffold of taurolithocholate orients towards the interior of the orthosteric binding site such that rings A and B are in contact with TM5 and TM6, the taurine side chain orients towards the extracellular opening of the binding site and forms a salt bridge with R79(EL1), and the 3-hydroxyl group forms hydrogen bonds with E169(5.44) and Y240(6.51). The binding mode thus differs in important aspects from the ones recently presented. These results are highly relevant for the development of novel, more potent agonists of TGR5 and should be a valuable starting point for the development of TGR5 antagonists, which could show antiproliferative effects in tumor cells.

Proline substitutions in a Mip-like peptidyl-prolyl cis-trans isomerase severely affect its structure, stability, shape and activity

FKBP22, an Escherichia coli-specific peptidyl-prolyl cis-trans isomerase, shows substantial homology with the Mip-like virulence factors. Mip-like proteins are homodimeric and possess a V-shaped conformation. Their N-terminal domains form dimers, whereas their C-terminal domains bind protein/peptide substrates and distinct inhibitors such as rapamycin and FK506. Interestingly, the two domains of the Mip-like proteins are separated by a lengthy, protease-susceptible α-helix. To delineate the structural requirement of this domain-connecting region in Mip-like proteins, we have investigated a recombinant FKBP22 (rFKBP22) and its three point mutants I65P, V72P and A82P using different probes. Each mutant harbors a Pro substitution mutation at a distinct location in the hinge region. We report that the three mutants are not only different from each other but also different from rFKBP22 in structure and activity. Unlike rFKBP22, the three mutants were unfolded by a non-two state mechanism in the presence of urea. In addition, the stabilities of the mutants, particularly I65P and V72P, differed considerably from that of rFKBP22. Conversely, the rapamycin binding affinity of no mutant was different from that of rFKBP22. Of the mutants, I65P showed the highest levels of structural/functional loss and dissociated partly in solution. Our computational study indicated a severe collapse of the V-shape in I65P due to the anomalous movement of its C-terminal domains. The α-helical nature of the domain-connecting region is, therefore, critical for the Mip-like proteins.

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

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

Marginally hydrophobic transmembrane α-helices shaping membrane protein folding

Cells have developed an incredible machinery to facilitate the insertion of membrane proteins into the membrane. While we have a fairly good understanding of the mechanism and determinants of membrane integration, more data is needed to understand the insertion of membrane proteins with more complex insertion and folding pathways. This review will focus on marginally hydrophobic transmembrane helices and their influence on membrane protein folding. These weakly hydrophobic transmembrane segments are by themselves not recognized by the translocon and therefore rely on local sequence context for membrane integration. How can such segments reside within the membrane? We will discuss this in the light of features found in the protein itself as well as the environment it resides in. Several characteristics in proteins have been described to influence the insertion of marginally hydrophobic helices. Additionally, the influence of biological membranes is significant. To begin with, the actual cost for having polar groups within the membrane may not be as high as expected; the presence of proteins in the membrane as well as characteristics of some amino acids may enable a transmembrane helix to harbor a charged residue. The lipid environment has also been shown to directly influence the topology as well as membrane boundaries of transmembrane helices-implying a dynamic relationship between membrane proteins and their environment.

Dynamic regulation of lipid-protein interactions

We review the importance of helix motions for the function of several important categories of membrane proteins and for the properties of several model molecular systems. For voltage-gated potassium or sodium channels, sliding, tilting and/or rotational movements of the S4 helix accompanied by a swapping of cognate side-chain ion-pair interactions regulate the channel gating. In the seven-helix G protein-coupled receptors, exemplified by the rhodopsins, collective helix motions serve to activate the functional signaling. Peptides which initially associate with lipid-bilayer membrane surfaces may undergo dynamic transitions from surface-bound to tilted-transmembrane orientations, sometimes accompanied by changes in the molecularity, formation of a pore or, more generally, the activation of biological function. For single-span membrane proteins, such as the tyrosine kinases, an interplay between juxtamembrane and transmembrane domains is likely to be crucial for the regulation of dimer assembly that in turn is associated with the functional responses to external signals. Additionally, we note that experiments with designed single-span transmembrane helices offer fundamental insights into the molecular features that govern protein-lipid interactions. This article is part of a Special Issue entitled: Lipid-protein interactions.

Conservation analysis of residues in the S4-S5 linker and the terminal part of the S5-P-S6 pore modulus in Kv and HCN channels: flexible determinants for the electromechanical coupling

Protein mobility is important to achieve protein function. Intrinsic flexibility associated with motion underlies this important issue and the analysis of side chain flexibility gives insights to understand it. In this work, the S5-P-S6 pore modulus (PM) of members of Kv and HCN channels was examined by a combination of sequence alignment, residue composition analysis, and intrinsic side chain flexibility. The PM sequences were organized as a database that was used to reveal and correlate the functional diversity of each analyzed family. Specifically, we focused our attention on the crucial role of the S4-S5 linker and its well-described interaction with the S6 T during the electromechanical coupling. Our analysis suggests the presence of a Gly-hinge in the middle of the S4-S5 linkers. This apparent Gly-hinge links a flexible N-terminal segment with a rigid C-terminal one, although in Kv7 channels, the latter segment is even more flexible. Instead, HCN channels exhibit a putative Thr-hinge and is rich in aromatic residues, in consequence, their linker is more rigid. Concerning S6, we confirm the presence of the two flexible kinks previously described and we provide the complete segmental flexibility profiles for the different families. Our results are discussed in terms of the relation between residue composition, conservation, and local conformational flexibility. This provides important insights to understand and differentiate the characteristic gating properties of these channels as well as their implications in cell physiology.