Studying Conformational Properties of Transmembrane Domain of KCNE3 in a Lipid Bilayer Membrane Using Molecular Dynamics Simulations

KCNE3 is a single-pass integral membrane protein that regulates numerous voltage-gated potassium channel functions such as KCNQ1. Previous solution NMR studies suggested a moderate degree of curved α-helical structure in the transmembrane domain (TMD) of KCNE3 in lyso-myristoylphosphatidylcholine (LMPC) micelles and isotropic bicelles with the residues T71, S74 and G78 situated along the concave face of the curved helix. During the interaction of KCNE3 and KCNQ1, KCNE3 pushes its transmembrane domain against KCNQ1 to lock the voltage sensor in its depolarized conformation. A cryo-EM study of KCNE3 complexed with KCNQ1 in nanodiscs suggested a deviation of the KCNE3 structure from its independent structure in isotropic bicelles. Despite the biological significance of KCNE3 TMD, the conformational properties of KCNE3 are poorly understood. Here, all atom molecular dynamics (MD) simulations were utilized to investigate the conformational dynamics of the transmembrane domain of KCNE3 in a lipid bilayer containing a mixture of POPC and POPG lipids (3:1). Further, the effect of the interaction impairing mutations (V72A, I76A and F68A) on the conformational properties of the KCNE3 TMD in lipid bilayers was investigated. Our MD simulation results suggest that the KCNE3 TMD adopts a nearly linear α helical structural conformation in POPC-POPG lipid bilayers. Additionally, the results showed no significant change in the nearly linear α-helical conformation of KCNE3 TMD in the presence of interaction impairing mutations within the sampled time frame. The KCNE3 TMD is more stable with lower flexibility in comparison to the N-terminal and C-terminal of KCNE3 in lipid bilayers. The overall conformational flexibility of KCNE3 also varies in the presence of the interaction-impairing mutations. The MD simulation data further suggest that the membrane bilayer width is similar for wild-type KCNE3 and KCNE3 containing mutations. The Z-distance measurement data revealed that the TMD residue site A69 is close to the lipid bilayer center, and residue sites S57 and S82 are close to the surfaces of the lipid bilayer membrane for wild-type KCNE3 and KCNE3 containing interaction-impairing mutations. These results agree with earlier KCNE3 biophysical studies. The results of these MD simulations will provide complementary data to the experimental outcomes of KCNE3 to help understand its conformational dynamic properties in a more native lipid bilayer environment.

How T118M peripheral myelin protein 22 predisposes humans to Charcot-Marie-Tooth disease

Data from gnomAD indicate that a missense mutation encoding the T118M variation in human peripheral myelin protein 22 (PMP22) is found in roughly one of every 75 genomes of western European lineage (1:120 in the overall human population). It is unusual among PMP22 variants that cause Charcot-Marie-Tooth (CMT) disease in that it is not 100% penetrant. Here, we conducted cellular and biophysical studies to determine why T118M PMP22 predisposes humans to CMT, but with only incomplete penetrance. We found that T118M PMP22 is prone to mistraffic but differs even from the WT protein in that increased expression levels do not result in a reduction in trafficking efficiency. Moreover, the T118M mutant exhibits a reduced tendency to form large intracellular aggregates relative to other disease mutants and even WT PMP22. NMR spectroscopy revealed that the structure and dynamics of T118M PMP22 resembled those of WT. These results show that the main consequence of T118M PMP22 in WT/T118M heterozygous individuals is a reduction in surface-trafficked PMP22, unaccompanied by formation of toxic intracellular aggregates. This explains the incomplete disease penetrance and the mild neuropathy observed for WT/T118M CMT cases. We also analyzed BioVU, a biobank linked to deidentified electronic medical records, and found a statistically robust association of the T118M mutation with the occurrence of long and/or repeated episodes of carpal tunnel syndrome. Collectively, our results illuminate the cellular effects of the T118M PMP22 variation leading to CMT disease and indicate a second disorder for which it is a risk factor.

Purification and membrane interactions of human KCNQ1100-370 potassium ion channel

KCNQ1 (Kv7.1 or KvLQT1) is a voltage-gated potassium ion channel that is involved in the ventricular repolarization following an action potential in the heart. It forms a complex with KCNE1 in the heart and is the pore forming subunit of slow delayed rectifier potassium current (I_ks). Mutations in KCNQ1, leading to a dysfunctional channel or loss of activity have been implicated in a cardiac disorder, long QT syndrome. In this study, we report the overexpression, purification, biochemical characterization of human KCNQ1_100-370, and lipid bilayer dynamics upon interaction with KCNQ1_100-370. The recombinant human KCNQ1 was expressed in Escherichia coli and purified into n-dodecylphosphocholine (DPC) micelles. The purified KCNQ1_100-370 was biochemically characterized by SDS-PAGE electrophoresis, western blot and nano-LC-MS/MS to confirm the identity of the protein. Circular dichroism (CD) spectroscopy was utilized to confirm the secondary structure of purified protein in vesicles. Furthermore, ³¹P and ²H solid-state NMR spectroscopy in DPPC/POPC/POPG vesicles (MLVs) indicated a direct interaction between KCNQ_100-370 and the phospholipid head groups. Finally, a visual inspection of KCNQ1_100-370 incorporated into MLVs was confirmed by transmission electron microscopy (TEM). The findings of this study provide avenues for future structural studies of the human KCNQ1 ion channel to have an in depth understanding of its structure-function relationship.

Comparing the structural dynamics of the human KCNE3 in reconstituted micelle and lipid bilayered vesicle environments

KCNE3 is a single transmembrane protein of the KCNE family that modulates the function and trafficking of several voltage-gated potassium channels, including KCNQ1. Structural studies of KCNE3 have been previously conducted in a wide range of model membrane mimics. However, it is important to assess the impact of the membrane mimics used on the observed conformation and dynamics. In this study, we have optimized a method for the reconstitution of the KCNE3 into POPC/POPG lipid bilayer vesicles for electron paramagnetic resonance (EPR) spectroscopy. Our CD spectroscopic data suggested that the degree of regular secondary structure for KCNE3 protein reconstituted into lipid bilayered vesicle is significantly higher than in DPC detergent micelles. Electron paramagnetic resonance (EPR) spectroscopy in combination with site-directed spin labeling (SDSL) was used to probe the structural dynamics of S49C, M59C, L67C, V85C, and S101C mutations of KCNE3 in both DPC micelles and in POPC/POPG lipid bilayered vesicles. Our CW-EPR power saturation data suggested that the site S74C is buried inside the lipid bilayered membrane while the site V85C is located outside the membrane, in contrast to DPC micelle results. These results suggest that the KCNE3 micelle structures need to be refined using data obtained in the lipid bilayered vesicles in order to ascertain the native structure of KCNE3. This work will provide guidelines for detailed structural studies of KCNE3 in a more native membrane environment and comparing the lipid bilayer results to the isotropic bicelle structure and to the KCNQ1-bound cryo-EM structure.

Investigating Structural Dynamics of KCNE3 in Different Membrane Environments Using Molecular Dynamics Simulations

KCNE3 is a potassium channel accessory transmembrane protein that regulates the function of various voltage-gated potassium channels such as KCNQ1. KCNE3 plays an important role in the recycling of potassium ion by binding with KCNQ1. KCNE3 can be found in the small intestine, colon, and in the human heart. Despite its biological significance, there is little information on the structural dynamics of KCNE3 in native-like membrane environments. Molecular dynamics (MD) simulations are a widely used as a tool to study the conformational dynamics and interactions of proteins with lipid membranes. In this study, we have utilized all-atom molecular dynamics simulations to characterize the molecular motions and the interactions of KCNE3 in a bilayer composed of: a mixture of POPC and POPG lipids (3:1), POPC alone, and DMPC alone. Our MD simulation results suggested that the transmembrane domain (TMD) of KCNE3 is less flexible and more stable when compared to the N- and C-termini of KCNE3 in all three membrane environments. The conformational flexibility of N- and C-termini varies across these three lipid environments. The MD simulation results further suggested that the TMD of KCNE3 spans the membrane width, having residue A69 close to the center of the lipid bilayers and residues S57 and S82 close to the lipid bilayer membrane surfaces. These results are consistent with previous biophysical studies of KCNE3. The outcomes of these MD simulations will help design biophysical experiments and complement the experimental data obtained on KCNE3 to obtain a more detailed understanding of its structural dynamics in the native membrane environment.

Predicting the functional impact of KCNQ1 variants with artificial neural networks

Recent advances in experimental and computational protein structure determination have provided access to high-quality structures for most human proteins and mutants thereof. However, linking changes in structure in protein mutants to functional impact remains an active area of method development. If successful, such methods can ultimately assist physicians in taking appropriate treatment decisions. This work presents three artificial neural network (ANN)-based predictive models that classify four key functional parameters of KCNQ1 variants as normal or dysfunctional using PSSM-based evolutionary and/or biophysical descriptors. Recent advances in predicting protein structure and variant properties with artificial intelligence (AI) rely heavily on the availability of evolutionary features and thus fail to directly assess the biophysical underpinnings of a change in structure and/or function. The central goal of this work was to develop an ANN model based on structure and physiochemical properties of KCNQ1 potassium channels that performs comparably or better than algorithms using only on PSSM-based evolutionary features. These biophysical features highlight the structure-function relationships that govern protein stability, function, and regulation. The input sensitivity algorithm incorporates the roles of hydrophobicity, polarizability, and functional densities on key functional parameters of the KCNQ1 channel. Inclusion of the biophysical features outperforms exclusive use of PSSM-based evolutionary features in predicting activation voltage dependence and deactivation time. As AI is increasingly applied to problems in biology, biophysical understanding will be critical with respect to 'explainable AI', i.e., understanding the relation of sequence, structure, and function of proteins. Our model is available at www.kcnq1predict.org.

High-Content Imaging Platform to Discover Chemical Modulators of Plasma Membrane Rafts

Plasma membrane organization profoundly impacts cellular functionality. A well-known mechanism underlying this organization is through nanoscopic clustering of distinct lipids and proteins in membrane rafts. Despite their physiological importance, rafts remain a difficult-to-study aspect of membrane organization, in part because of the paucity of chemical tools to experimentally modulate their properties. Methods to selectively target rafts for therapeutic purposes are also currently lacking. To tackle these problems, we developed a high-throughput screen and an accompanying image analysis pipeline to identify small molecules that enhance or inhibit raft formation. Cell-derived giant plasma membrane vesicles were used as the experimental platform. A proof-of-principle screen using a bioactive lipid library demonstrates that this method is robust and capable of validating established raft modulators including C6- and C8-ceramide, miltefosine, and epigallocatechin gallate as well as identifying new ones. The platform we describe here represents a powerful tool to discover new chemical approaches to manipulate rafts and their components.

A Model for the Signal Initiation Complex Between Arrestin-3 and the Src Family Kinase Fgr

Arrestins regulate a wide range of signaling events, most notably when bound to active G protein-coupled receptors (GPCRs). Among the known effectors recruited by GPCR-bound arrestins are Src family kinases, which regulate cellular growth and proliferation. Here, we focus on arrestin-3 interactions with Fgr kinase, a member of the Src family. Previous reports demonstrated that Fgr exhibits high constitutive activity, but can be further activated by both arrestin-dependent and arrestin-independent pathways. We report that arrestin-3 modulates Fgr activity with a hallmark bell-shaped concentration-dependence, consistent with a role as a signaling scaffold. We further demonstrate using NMR spectroscopy that a polyproline motif within arrestin-3 interacts directly with the SH3 domain of Fgr. To provide a framework for this interaction, we determined the crystal structure of the Fgr SH3 domain at 1.9 Å resolution and developed a model for the GPCR-arrestin-3-Fgr complex that is supported by mutagenesis. This model suggests that Fgr interacts with arrestin-3 at multiple sites and is consistent with the locations of disease-associated Fgr mutations. Collectively, these studies provide a structural framework for arrestin-dependent activation of Fgr.

Compendium of causative genes and their encoded proteins for common monogenic disorders

A compendium is presented of inherited monogenic disorders that have a prevalence of >1:20,000 in the human population, along with their causative genes and encoded proteins. "Simple" monogenic diseases are those for which the clinical features are caused by mutations impacting a single gene, usually in a manner that alters the sequence of the encoded protein. Of course, for a given "monogenic disorder", there is sometimes more than one potential disease gene, mutations in any one of which is sufficient to cause phenotypes of that disorder. Disease-causing mutations for monogenic disorders are usually passed on from generation to generation in a Mendelian fashion, and originate from spontaneous (de novo) germline founder mutations. In the past monogenic disorders have often been written off as targets for drug discovery because they sometimes are assumed to be rare disorders, for which the meager projected financial payoff of drug discovery and development has discouraged investment. However, not all monogenic diseases are rare. Here, we report that that currently available data identifies 72 disorders with a prevalence of at least 1 in 20,000 humans. For each, we tabulate the gene(s) for which mutations cause the spectrum of phenotypes associated with that disorder. We also identify the gene and protein that most commonly causes each disease. 34 of these disorders are caused exclusively by mutations in only a single gene and encoded protein.

The transmembrane amyloid precursor C99 protein exhibits non-specific interaction with tau

Historically, the two most prominent proteins in Alzheimer's disease (AD) research have been the amyloid precursor protein (APP) and the microtubule assembly protein tau. In the classical model for the etiology of AD, amyloid-β (Aβ)-an APP derivative and hyperphosphorylated tau form aggregates in the brain that underlie the pathogenesis of the disease. However, the connection between Aβ and tau pathologies remains unclear. Several studies have provided evidence that the presence of Aβ can induce or enhance neurofibrillary tangle formation by tau. Others have reported a direct interaction between tau and short fragments of the APP transmembrane domain, C99. Structural studies of C99 show that these in vitro tau-binding fragments of C99 are buried in the lipid bilayer and are likely unavailable to bind tau in vivo. Given the importance of APP and tau in AD, we sought to characterize the potential interaction of the Aβ precursor, full length C99, and tau in vitro using NMR spectroscopy. We found that C99 and soluble tau interact only weakly and, most likely, non-specifically.

Recombinant SARS-CoV-2 envelope protein traffics to the trans-Golgi network following amphipol-mediated delivery into human cells

The severe acute respiratory syndrome coronavirus 2 envelope protein (S2-E) is a conserved membrane protein that is important for coronavirus (CoV) assembly and budding. Here, we describe the recombinant expression and purification of S2-E in amphipol-class amphipathic polymer solutions, which solubilize and stabilize membrane proteins, but do not disrupt membranes. We found that amphipol delivery of S2-E to preformed planar bilayers results in spontaneous membrane integration and formation of viroporin cation channels. Amphipol delivery of the S2-E protein to human cells results in plasma membrane integration, followed by retrograde trafficking to the trans-Golgi network and accumulation in swollen perinuclear lysosomal-associated membrane protein 1-positive vesicles, likely lysosomes. CoV envelope proteins have previously been proposed to manipulate the luminal pH of the trans-Golgi network, which serves as an accumulation station for progeny CoV particles prior to cellular egress via lysosomes. Delivery of S2-E to cells will enable chemical biological approaches for future studies of severe acute respiratory syndrome coronavirus 2 pathogenesis and possibly even development of "Trojan horse" antiviral therapies. Finally, this work also establishes a paradigm for amphipol-mediated delivery of membrane proteins to cells.

Glycosylation limits forward trafficking of the tetraspan membrane protein PMP22

Peripheral myelin protein 22 (PMP22) folds and trafficks inefficiently, with only 20% of newly expressed protein trafficking to the cell surface. This behavior is exacerbated in many of the mutants associated with Charcot–Marie–Tooth disease, motivating further study. Here we characterized the role of N-glycosylation in limiting PMP22 trafficking. We first eliminated N-glycosylation using an N41Q mutation, which resulted in an almost 3-fold increase in trafficking efficiency of wildtype (WT) PMP22 and a 10-fold increase for the severely unstable L16P disease mutant in HEK293 cells, with similar results in Schwann cells. Total cellular levels were also much higher for the WT/N41Q mutant, although not for the L16P/N41Q form. Depletion of oligosaccharyltransferase OST-A and OST-B subunits revealed that WT PMP22 is N-glycosylated posttranslationally by OST-B, whereas L16P is cotranslationally glycosylated by OST-A. Quantitative proteomic screens revealed similarities and differences in the interactome for WT, glycosylation-deficient, and unstable mutant forms of PMP22 and also suggested that L16P is sequestered at earlier stages of endoplasmic reticulum quality control. CRISPR knockout studies revealed a role for retention in endoplasmic reticulum sorting receptor 1 (RER1) in limiting the trafficking of all three forms, for UDP-glucose glycoprotein glucosyltransferase 1 (UGGT1) in limiting the trafficking of WT and L16P but not N41Q, and calnexin (CNX) in limiting the trafficking of WT and N41Q but not L16P. This work shows that N-glycosylation is a limiting factor to forward trafficking PMP22 and sheds light on the proteins involved in its quality control.

The C99 domain of the amyloid precursor protein resides in the disordered membrane phase

Processing of the amyloid precursor protein (APP) via the amyloidogenic pathway is associated with the etiology of Alzheimer's disease. The cleavage of APP by β-secretase to generate the transmembrane 99-residue C-terminal fragment (C99) and subsequent processing of C99 by γ-secretase to yield amyloid-β (Aβ) peptides are essential steps in this pathway. Biochemical evidence suggests that amyloidogenic processing of C99 occurs in cholesterol- and sphingolipid-enriched liquid-ordered phase membrane rafts. However, direct evidence that C99 preferentially associates with these rafts has remained elusive. Here, we tested this by quantifying the affinity of C99-GFP for raft domains in cell-derived giant plasma membrane vesicles (GPMVs). We found that C99 was essentially excluded from ordered domains in vesicles from HeLa cells, undifferentiated SH-SY5Y cells, or SH-SY5Y-derived neurons; instead, ∼90% of C99 partitioned into disordered domains. The strong association of C99 with disordered domains occurred independently of its cholesterol-binding activity or homodimerization, or of the presence of the familial Alzheimer disease Arctic mutation (APP E693G). Finally, through biochemical studies we confirmed previous results, which showed that C99 is processed in the plasma membrane by α-secretase, in addition to the well-known γ-secretase. These findings suggest that C99 itself lacks an intrinsic affinity for raft domains, implying that either i) amyloidogenic processing of the protein occurs in disordered regions of the membrane, ii) processing involves a marginal subpopulation of C99 found in rafts, or iii) as-yet-unidentified protein-protein interactions with C99 in living cells drive this protein into membrane rafts to promote its cleavage therein.

Structural determinants of cholesterol recognition in helical integral membrane proteins

Cholesterol is an integral component of mammalian membranes. It has been shown to modulate membrane fluidity and dynamics and alter integral membrane protein function. However, understanding the molecular mechanisms of how cholesterol impacts protein function is complicated by limited and conflicting structural data. Because of the nature of the crystallization and cryo-EM structure determination, it is difficult to distinguish between specific and biologically relevant interactions and a nonspecific association. The only widely recognized search algorithm for cholesterol-integral-membrane-protein interaction sites is sequence based, i.e., searching for the so-called "Cholesterol Recognition/interaction Amino acid Consensus" motif. Although these motifs are present in numerous integral membrane proteins, there is inconclusive evidence to support their necessity or sufficiency for cholesterol binding. Here, we leverage the increasing number of experimental cholesterol-integral-membrane-protein structures to systematically analyze putative interaction sites based on their spatial arrangement and evolutionary conservation. This analysis creates three-dimensional representations of general cholesterol interaction sites that form clusters across multiple integral membrane protein classes. We also classify cholesterol-integral-membrane-protein interaction sites as either likely-specific or nonspecific. Information gleaned from our characterization will eventually enable a structure-based approach to predict and design cholesterol-integral-membrane-protein interaction sites.

Disease-linked supertrafficking of a potassium channel

Gain-of-function (GOF) mutations in the voltage-gated potassium channel subfamily Q member 1 (KCNQ1) can induce cardiac arrhythmia. In this study, it was tested whether any of the known human GOF disease mutations in KCNQ1 act by increasing the amount of KCNQ1 that reaches the cell surface-"supertrafficking." Seven of the 15 GOF mutants tested were seen to surface traffic more efficiently than the WT channel. Among these, we found that the levels of R231C KCNQ1 in the plasma membrane were fivefold higher than the WT channel. This was shown to arise from the combined effects of enhanced efficiency of translocon-mediated membrane integration of the S4 voltage-sensor helix and from enhanced post-translational folding/trafficking related to the energetic linkage of C231 with the V129 and F166 side chains. Whole-cell electrophysiology recordings confirmed that R231C KCNQ1 in complex with the voltage-gated potassium channel-regulatory subfamily E member 1 not only exhibited constitutive conductance but also revealed that the single-channel activity of this mutant is only 20% that of WT. The GOF phenotype associated with R231C therefore reflects the effects of supertrafficking and constitutive channel activation, which together offset reduced channel activity. These investigations show that membrane protein supertrafficking can contribute to human disease.

Delivery of recombinant SARS-CoV-2 envelope protein into human cells

SARS-CoV-2 envelope protein (S2-E) is a conserved membrane protein that is essential to coronavirus assembly and budding. Here, we describe the recombinant expression and purification of S2-E into amphipol-class amphipathic polymer solutions. The physical properties of amphipols underpin their ability to solubilize and stabilize membrane proteins without disrupting membranes. Amphipol delivery of S2-E to pre-formed planar bilayers results in spontaneous membrane integration and formation of viroporin ion channels. Amphipol delivery of the S2-E protein to human cells results in membrane integration followed by retrograde trafficking to a location adjacent to the endoplasmic reticulum-to-Golgi intermediate compartment (ERGIC) and the Golgi, which are the sites of coronavirus replication. Delivery of S2-E to cells enables both chemical biological approaches for future studies of SARS-CoV-2 pathogenesis and development of "Trojan Horse" anti-viral therapies. This work also establishes a paradigm for amphipol-mediated delivery of membrane proteins to cells.

Collision-Induced Unfolding Differentiates Functional Variants of the KCNQ1 Voltage Sensor Domain

The KCNQ1 voltage-gated potassium channel regulates the repolarization of cardiac cells, and a plurality of point mutations in its voltage-sensing domain (VSD) are associated with toxic gain or loss of pore function, resulting in disease. As is the case with many disease-associated membrane proteins, there are hundreds of human variants of interest identified for KCNQ1; however, a significant portion of these variants have not been characterized in relation to their functional and disease associations. Additionally, as the VSD consists of four transmembrane helices, studies into dynamic structural differences among KCNQ1 VSD variants are hindered by the current limitations and deficits in the high-resolution structure determination of membrane proteins. Here, we use native ion mobility-mass spectrometry and collision-induced unfolding (CIU) to address the need for a high throughput-compatible method for the structural characterization of membrane protein variants of unknown significance using the KCNQ1 VSD as a model system. We perform CIU on wild-type and three mutant KCNQ1 VSD forms associated with the toxic gain or loss of function and show through both automated feature detection and comprehensive difference analysis of the CIU data sets that the variants are clearly grouped by function and disease association. We also construct a classification scheme based on the CIU data sets, which is able to differentiate the variant functional groups and classify a recently characterized variant to its correct grouping. Further, we probe the stability of the KCNQ1 VSD variants when liberated from C12E8 micelles at pH 8.0 and find preliminary evidence that the R231C mutation associated with the gain of the pore function is destabilized relative to the wild-type and loss of function variants.

Allosteric mechanism for KCNE1 modulation of KCNQ1 potassium channel activation

The function of the voltage-gated KCNQ1 potassium channel is regulated by co-assembly with KCNE auxiliary subunits. KCNQ1-KCNE1 channels generate the slow delayed rectifier current, I_Ks, which contributes to the repolarization phase of the cardiac action potential. A three amino acid motif (F57-T58-L59, FTL) in KCNE1 is essential for slow activation of KCNQ1-KCNE1 channels. However, how this motif interacts with KCNQ1 to control its function is unknown. Combining computational modeling with electrophysiological studies, we developed structural models of the KCNQ1-KCNE1 complex that suggest how KCNE1 controls KCNQ1 activation. The FTL motif binds at a cleft between the voltage-sensing and pore domains and appears to affect the channel gate by an allosteric mechanism. Comparison with the KCNQ1-KCNE3 channel structure suggests a common transmembrane-binding mode for different KCNEs and illuminates how specific differences in the interaction of their triplet motifs determine the profound differences in KCNQ1 functional modulation by KCNE1 versus KCNE3.

Direct relationship between increased expression and mistrafficking of the Charcot-Marie-Tooth-associated protein PMP22

Charcot-Marie-Tooth disease (CMT) is a neuropathy of the peripheral nervous system that afflicts ∼1:2500 people. The most common form of this disease (CMT1A, 1:4000) is associated with duplication of chromosome fragment 17p11.2-12, which results in a third WT PMP22 allele. In rodent models overexpressing the PMP22 (peripheral myelin protein 22) protein and in dermal fibroblasts from patients with CMT1A, PMP22 aggregates have been observed. This suggests that overexpression of PMP22 under CMT1A conditions overwhelms the endoplasmic reticulum quality control system, leading to formation of cytotoxic aggregates. In this work, we used a single-cell flow-cytometry trafficking assay to quantitatively examine the relationship between PMP22 expression and trafficking efficiency in individual cells. We observed that as expression of WT or disease variants of PMP22 is increased, the amount of intracellular PMP22 increases to a greater extent than the amount of surface-trafficked protein. This was true for both transiently transfected cells and PMP22 stable expressing cells. Our results support the notion that overexpression of PMP22 in CMT1A leads to a disproportionate increase in misfolding and mistrafficking of PMP22, which is likely a contributor to disease pathology and progression.

Bicelles Rich in both Sphingolipids and Cholesterol and Their Use in Studies of Membrane Proteins

How the distinctive lipid composition of mammalian plasma membranes impacts membrane protein structure is largely unexplored, partly because of the dearth of isotropic model membrane systems that contain abundant sphingolipids and cholesterol. This gap is addressed by showing that sphingomyelin and cholesterol-rich (SCOR) lipid mixtures with phosphatidylcholine can be cosolubilized by n-dodecyl-β-melibioside to form bicelles. Small-angle X-ray and neutron scattering, as well as cryo-electron microscopy, demonstrate that these assemblies are stable over a wide range of conditions and exhibit the bilayered-disc morphology of ideal bicelles even at low lipid-to-detergent mole ratios. SCOR bicelles are shown to be compatible with a wide array of experimental techniques, as applied to the transmembrane human amyloid precursor C99 protein in this medium. These studies reveal an equilibrium between low-order oligomer structures that differ significantly from previous experimental structures of C99, providing an example of how ordered membranes alter membrane protein structure.

Structures Illuminate Cardiac Ion Channel Functions in Health and in Long QT Syndrome

The cardiac action potential is critical to the production of a synchronized heartbeat. This electrical impulse is governed by the intricate activity of cardiac ion channels, among them the cardiac voltage-gated potassium (Kv) channels KCNQ1 and hERG as well as the voltage-gated sodium (Nav) channel encoded by SCN5A. Each channel performs a highly distinct function, despite sharing a common topology and structural components. These three channels are also the primary proteins mutated in congenital long QT syndrome (LQTS), a genetic condition that predisposes to cardiac arrhythmia and sudden cardiac death due to impaired repolarization of the action potential and has a particular proclivity for reentrant ventricular arrhythmias. Recent cryo-electron microscopy structures of human KCNQ1 and hERG, along with the rat homolog of SCN5A and other mammalian sodium channels, provide atomic-level insight into the structure and function of these proteins that advance our understanding of their distinct functions in the cardiac action potential, as well as the molecular basis of LQTS. In this review, the gating, regulation, LQTS mechanisms, and pharmacological properties of KCNQ1, hERG, and SCN5A are discussed in light of these recent structural findings.

Structure and physiological function of the human KCNQ1 channel voltage sensor intermediate state

Voltage-gated ion channels feature voltage sensor domains (VSDs) that exist in three distinct conformations during activation: resting, intermediate, and activated. Experimental determination of the structure of a potassium channel VSD in the intermediate state has previously proven elusive. Here, we report and validate the experimental three-dimensional structure of the human KCNQ1 voltage-gated potassium channel VSD in the intermediate state. We also used mutagenesis and electrophysiology in Xenopus laevisoocytes to functionally map the determinants of S4 helix motion during voltage-dependent transition from the intermediate to the activated state. Finally, the physiological relevance of the intermediate state KCNQ1 conductance is demonstrated using voltage-clamp fluorometry. This work illuminates the structure of the VSD intermediate state and demonstrates that intermediate state conductivity contributes to the unusual versatility of KCNQ1, which can function either as the slow delayed rectifier current (IKs) of the cardiac action potential or as a constitutively active epithelial leak current.

High-Throughput Functional Evaluation of KCNQ1 Decrypts Variants of Unknown Significance

BACKGROUND

The explosive growth in known human gene variation presents enormous challenges to current approaches for variant classification that have implications for diagnosis and treatment of many genetic diseases. For disorders caused by mutations in cardiac ion channels as in congenital arrhythmia syndromes, in vitro electrophysiological evidence has high value in discriminating pathogenic from benign variants, but these data are often lacking because assays are cost, time, and labor intensive.

METHODS

We implemented a strategy for performing high-throughput functional evaluations of ion channel variants that repurposed an automated electrophysiological recording platform developed previously for drug discovery.

RESULTS

We demonstrated the success of this approach by evaluating 78 variants in KCNQ1, a major gene involved in genetic disorders of cardiac arrhythmia susceptibility. We benchmarked our results with traditional electrophysiological approaches and observed a high level of concordance. This strategy also enabled studies of dominant-negative behavior of variants exhibiting severe loss-of-function. Overall, our results provided functional data useful for reclassifying >65% of the studied KCNQ1 variants.

CONCLUSIONS

Our results illustrate an efficient and high-throughput paradigm linking genotype to function for a human cardiac ion channel that will enable data-driven classification of large numbers of variants and create new opportunities for precision medicine.

Upgraded molecular models of the human KCNQ1 potassium channel

The voltage-gated potassium channel KCNQ1 (K_V7.1) assembles with the KCNE1 accessory protein to generate the slow delayed rectifier current, I_KS, which is critical for membrane repolarization as part of the cardiac action potential. Loss-of-function (LOF) mutations in KCNQ1 are the most common cause of congenital long QT syndrome (LQTS), type 1 LQTS, an inherited genetic predisposition to cardiac arrhythmia and sudden cardiac death. A detailed structural understanding of KCNQ1 is needed to elucidate the molecular basis for KCNQ1 LOF in disease and to enable structure-guided design of new anti-arrhythmic drugs. In this work, advanced structural models of human KCNQ1 in the resting/closed and activated/open states were developed by Rosetta homology modeling guided by newly available experimentally-based templates: X. leavis KCNQ1 and various resting voltage sensor structures. Using molecular dynamics (MD) simulations, the capacity of the models to describe experimentally established channel properties including state-dependent voltage sensor gating charge interactions and pore conformations, PIP2 binding sites, and voltage sensor–pore domain interactions were validated. Rosetta energy calculations were applied to assess the utility of each model in interpreting mutation-evoked KCNQ1 dysfunction by predicting the change in protein thermodynamic stability for 50 experimentally characterized KCNQ1 variants with mutations located in the voltage-sensing domain. Energetic destabilization was successfully predicted for folding-defective KCNQ1 LOF mutants whereas wild type-like mutants exhibited no significant energetic frustrations, which supports growing evidence that mutation-induced protein destabilization is an especially common cause of KCNQ1 dysfunction. The new KCNQ1 Rosetta models provide helpful tools in the study of the structural basis for KCNQ1 function and can be used to generate hypotheses to explain KCNQ1 dysfunction.

Life During Wartime: A Personal Recollection of the Circa 1990 Prestegard Lab and Its Contributions to Membrane Biophysics

A subjective account is presented of challenges and excitement of being a postdoctoral trainee in the lab of James H. Prestegard at Yale University in New Haven, Connecticut from 1989 to 1991. This includes accounts of the early development of bicelles and of oriented sample NMR results that contributed to our modern understanding of the properties of the water-lipid interface of disordered phase biological membranes.

A unified structural model of the mammalian translocator protein (TSPO)

The translocator protein (TSPO), previously known as the peripheral benzodiazepine receptor (PBR), is a membrane protein located on the outer mitochondrial membrane. Experimentally-derived structures of mouse TSPO (mTSPO) and its homologs from bacterial species have been determined by NMR spectroscopy and X-ray crystallography, respectively. These structures and ligand interactions within the TSPO binding pocket display distinct differences. Here, we leverage experimental and computational studies to derive a unified structural model of mTSPO in the presence and absence of the TSPO ligand, PK11195, and study the effects of DPC detergent micelles on the TSPO structure and ligand binding. From this work, we conclude that that the lipid-mimetic system used to solubilize mTSPO for NMR studies thermodynamically destabilizes the protein, introduces structural perturbations, and alters the characteristics of ligand binding. Furthermore, we used Rosetta to construct a unified mTSPO model that reconciles deviating features of the mammalian and bacterial TSPO. These deviating features are likely a consequence of the detergent system used for structure determination of mTSPO by NMR. The unified mTSPO model agrees with available experimental NMR data, appears to be physically realistic (i.e. thermodynamically not frustrated as judged by the Rosetta energy function), and simultaneously shares the structural features observed in sequence-conserved regions of the bacterial proteins. Finally, we identified the binding site for an imaging ligand VUIIS8310 that is currently positioned for clinical translation using NMR spectroscopy and propose a computational model of the VUIIS8310-mTSPO complex.

Peripheral myelin protein 22 modulates store-operated calcium channel activity, providing insights into Charcot-Marie-Tooth disease etiology

Charcot-Marie-Tooth (CMT) disease is a peripheral neuropathy associated with gene duplication and point mutations in the peripheral myelin protein 22 (PMP22) gene. However, the role of PMP22 in Schwann cell physiology and the mechanisms by which PMP22 mutations cause CMT are not well-understood. On the basis of homology between PMP22 and proteins associated with modulation of ion channels, we hypothesized that PMP22 alters ion channel activity. Using whole-cell electrophysiology, we show here that heterologous PMP22 expression increases the amplitude of currents similar to those ascribed to store-operated calcium (SOC) channels, particularly those involving transient receptor canonical channel 1 (TrpC1). These channels help replenish Ca²⁺ in the endoplasmic reticulum (ER) following stimulus-induced depletion. Currents with similar properties were recorded in WT but not pmp22 ^-/- mouse Schwann cells. Heterologous expression of the CMT-associated PMP22_L16P variant, which fails to reach the plasma membrane and localizes to the ER, led to larger currents than WT PMP22. Similarly, Schwann cells isolated from Trembler J (TrJ; PMP22_L16P) mice had larger currents than WT littermates. Calcium imaging in live nerves and cultured Schwann cells revealed elevated intracellular Ca²⁺ in TrJ mice compared with WT. Moreover, we found that PMP22 co-immunoprecipitated with stromal interaction molecule 1 (STIM1), the Ca²⁺ sensor SOC channel subunit in the ER. These results suggest that in the ER, PMP22 interacts with STIM1 and increases Ca²⁺ influx through SOC channels. Excess or mutant PMP22 in the ER may elevate intracellular Ca²⁺ levels, which could contribute to CMT pathology.

Folding and Misfolding of Human Membrane Proteins in Health and Disease: From Single Molecules to Cellular Proteostasis

Advances over the past 25 years have revealed much about how the structural properties of membranes and associated proteins are linked to the thermodynamics and kinetics of membrane protein (MP) folding. At the same time biochemical progress has outlined how cellular proteostasis networks mediate MP folding and manage misfolding in the cell. When combined with results from genomic sequencing, these studies have established paradigms for how MP folding and misfolding are linked to the molecular etiologies of a variety of diseases. This emerging framework has paved the way for the development of a new class of small molecule "pharmacological chaperones" that bind to and stabilize misfolded MP variants, some of which are now in clinical use. In this review, we comprehensively outline current perspectives on the folding and misfolding of integral MPs as well as the mechanisms of cellular MP quality control. Based on these perspectives, we highlight new opportunities for innovations that bridge our molecular understanding of the energetics of MP folding with the nuanced complexity of biological systems. Given the many linkages between MP misfolding and human disease, we also examine some of the exciting opportunities to leverage these advances to address emerging challenges in the development of therapeutics and precision medicine.

The vexing complexity of the amyloidogenic pathway

The role of the amyloidogenic pathway in the etiology of Alzheimer's disease (AD), particularly the common sporadic late onset forms of the disease, is controversial. To some degree, this is a consequence of the failure of drug and therapeutic antibody trials based either on targeting the proteases in this pathway or its amyloid end products. Here, we explore the formidable complexity of the biochemistry and cell biology associated with this pathway. For example, we review evidence that the immediate precursor of amyloid-β, the C99 domain of the amyloid precursor protein (APP), may itself be toxic. We also review important new results that appear to finally establish a direct genetic link between mutations in APP and the sporadic forms of AD. Based on the complexity of amyloidogenesis, it seems possible that a major contributor to the failure of related drug trials is that we have an incomplete understanding of this pathway and how it is linked to Alzheimer's pathogenesis. If so, this highlights a need for further characterization of this pathway, not its abandonment.

NMR resonance assignments and secondary structure of a mutant form of the human KCNE1 channel accessory protein that exhibits KCNE3-like function

KCNQ1 (Q1) is a voltage-gated potassium channel that is modulated by members of the KCNE family, the best-characterized being KCNE1 (E1) and KCNE3 (E3). The Q1/E1 complex generates a channel with delayed activation and increased conductance. This complex is expressed in cardiomyocytes where it provides the I_Ks current that is critical for the repolarization phase of the cardiac action potential. The Q1/E3 complex, on the other hand, is expressed in epithelial cells of the colon and stomach, where it serves as a constitutively active leak channel to help maintain water and ion homeostasis. Studies show the single transmembrane segments (TMS) present in both E1 and E3 are essential to their distinct functions. More specifically, residues FTL located near the middle of the E1 TMS are essential for the delayed activation of Q1, while the corresponding TVG sites in E3 are critical for constitutive activation of the channel. Swapping these three residues leads to the switching of the functional properties for both Q1/E1_FTL→TVG and Q1/E3_TVG→FTL complexes. This work details the backbone assignments and chemical shifts for the E1_FTL→TVG mutant, as determined using a suite of 3D NMR experiments along with specific and inverse amino acid isotopic labeling. The completed assignments can be used, in conjunction with other NMR experiments, to generate a 3D structure of E1_FTL→TVG. The results of TALOS-N analysis of the chemical shifts are reported here. The E1_FTL→TVG structure will be compared to the already available E1 and E3 structures to determine the roles that their TMS triplet motifs play in each protein to dictate their distinct channel-modulatory functions.

Protein structure aids predicting functional perturbation of missense variants in SCN5A and KCNQ1

Rare variants in the cardiac potassium channel K_V7.1 (KCNQ1) and sodium channel Na_V1.5 (SCN5A) are implicated in genetic disorders of heart rhythm, including congenital long QT and Brugada syndromes (LQTS, BrS), but also occur in reference populations. We previously reported two sets of Na_V1.5 (n = 356) and K_V7.1 (n = 144) variants with in vitro characterized channel currents gathered from the literature. Here we investigated the ability to predict commonly reported Na_V1.5 and K_V7.1 variant functional perturbations by leveraging diverse features including variant classifiers PROVEAN, PolyPhen-2, and SIFT; evolutionary rate and BLAST position specific scoring matrices (PSSM); and structure-based features including “functional densities” which is a measure of the density of pathogenic variants near the residue of interest. Structure-based functional densities were the most significant features for predicting Na_V1.5 peak current (adj. R² = 0.27) and K_V7.1 + KCNE1 half-maximal voltage of activation (adj. R² = 0.29). Additionally, use of structure-based functional density values improves loss-of-function classification of SCN5A variants with an ROC-AUC of 0.78 compared with other predictive classifiers (AUC = 0.69; two-sided DeLong test p = .01). These results suggest structural data can inform predictions of the effect of uncharacterized SCN5A and KCNQ1 variants to provide a deeper understanding of their burden on carriers.

Probing the Dynamics and Structural Topology of the Reconstituted Human KCNQ1 Voltage Sensor Domain (Q1-VSD) in Lipid Bilayers Using Electron Paramagnetic Resonance Spectroscopy

KCNQ1 (Kv7.1 or KvLQT1) is a potassium ion channel protein found in the heart, ear, and other tissues. In complex with the KCNE1 accessory protein, it plays a role during the repolarization phase of the cardiac action potential. Mutations in the channel have been associated with several diseases, including congenital deafness and long QT syndrome. Nuclear magnetic resonance (NMR) structural studies in detergent micelles and a cryo-electron microscopy structure of KCNQ1 from Xenopus laevis have shown that the voltage sensor domain (Q1-VSD) of the channel has four transmembrane helices, S1-S4, being overall structurally similar with other VSDs. In this study, we describe a reliable method for the reconstitution of Q1-VSD into (POPC/POPG) lipid bilayer vesicles. Site-directed spin labeling electron paramagnetic resonance spectroscopy was used to probe the structural dynamics and topology of several residues of Q1-VSD in POPC/POPG lipid bilayer vesicles. Several mutants were probed to determine their location and corresponding immersion depth (in angstroms) with respect to the membrane. The dynamics of the bilayer vesicles upon incorporation of Q1-VSD were studied using ³¹P solid-state NMR spectroscopy by varying the protein:lipid molar ratios confirming the interaction of the protein with the bilayer vesicles. Circular dichroism spectroscopic data showed that the α-helical content of Q1-VSD is higher for the protein reconstituted in vesicles than in previous studies using DPC detergent micelles. This study provides insight into the structural topology and dynamics of Q1-VSD reconstituted in a lipid bilayer environment, forming the basis for more advanced structural and functional studies.

LCP1 preferentially binds clasped αMβ2 integrin and attenuates leukocyte adhesion under flow

Integrins are α/β heterodimers that interconvert between inactive and active states. In the active state the α/β cytoplasmic domains recruit integrin-activating proteins and separate the transmembrane and cytoplasmic (TMcyto) domains (unclasped TMcyto). Conversely, in the inactive state the α/β TMcyto domains bind integrin-inactivating proteins, resulting in the association of the TMcyto domains (clasped TMcyto). Here, we report the isolation of integrin cytoplasmic tail interactors using either lipid bicelle-incorporated integrin TMcyto domains (α5, αM, αIIb, β1, β2 and β3 integrin TMcyto) or a clasped, lipid bicelle-incorporated αMβ2 TMcyto. Among the proteins found to preferentially bind clasped rather than the isolated αM and β2 subunits was L-plastin (LCP1, also known as plastin-2), which binds to and maintains the inactive state of αMβ2 integrin in vivo and thereby regulates leukocyte adhesion to integrin ligands under flow. Our findings offer a global view on cytoplasmic proteins interacting with different integrins and provide evidence for the existence of conformation-specific integrin interactors.

Bexarotene Binds to the Amyloid Precursor Protein Transmembrane Domain, Alters Its α-Helical Conformation, and Inhibits γ-Secretase Nonselectively in Liposomes

Bexarotene is a pleiotropic molecule that has been proposed as an amyloid-β (Aβ)-lowering drug for the treatment of Alzheimer's disease (AD). It acts by upregulation of an apolipoprotein E (apoE)-mediated Aβ clearance mechanism. However, whether bexarotene induces removal of Aβ plaques in mouse models of AD has been controversial. Here, we show by NMR and CD spectroscopy that bexarotene directly interacts with and stabilizes the transmembrane domain α-helix of the amyloid precursor protein (APP) in a region where cholesterol binds. This effect is not mediated by changes in membrane lipid packing, as bexarotene does not share with cholesterol the property of inducing phospholipid condensation. Bexarotene inhibited the intramembrane cleavage by γ-secretase of the APP C-terminal fragment C99 to release Aβ in cell-free assays of the reconstituted enzyme in liposomes, but not in cells, and only at very high micromolar concentrations. Surprisingly, in vitro, bexarotene also inhibited the cleavage of Notch1, another major γ-secretase substrate, demonstrating that its inhibition of γ-secretase is not substrate specific and not mediated by acting via the cholesterol binding site of C99. Our data suggest that bexarotene is a pleiotropic molecule that interfere with Aβ metabolism through multiple mechanisms.

Predicting the Functional Impact of KCNQ1 Variants of Unknown Significance

BACKGROUND

An emerging standard-of-care for long-QT syndrome uses clinical genetic testing to identify genetic variants of the KCNQ1 potassium channel. However, interpreting results from genetic testing is confounded by the presence of variants of unknown significance for which there is inadequate evidence of pathogenicity.

METHODS AND RESULTS

In this study, we curated from the literature a high-quality set of 107 functionally characterized KCNQ1 variants. Based on this data set, we completed a detailed quantitative analysis on the sequence conservation patterns of subdomains of KCNQ1 and the distribution of pathogenic variants therein. We found that conserved subdomains generally are critical for channel function and are enriched with dysfunctional variants. Using this experimentally validated data set, we trained a neural network, designated Q1VarPred, specifically for predicting the functional impact of KCNQ1 variants of unknown significance. The estimated predictive performance of Q1VarPred in terms of Matthew's correlation coefficient and area under the receiver operating characteristic curve were 0.581 and 0.884, respectively, superior to the performance of 8 previous methods tested in parallel. Q1VarPred is publicly available as a web server at http://meilerlab.org/q1varpred.

CONCLUSIONS

Although a plethora of tools are available for making pathogenicity predictions over a genome-wide scale, previous tools fail to perform in a robust manner when applied to KCNQ1. The contrasting and favorable results for Q1VarPred suggest a promising approach, where a machine-learning algorithm is tailored to a specific protein target and trained with a functionally validated data set to calibrate informatics tools.