Effects of calreticulin mutations on cell transformation and immunity.

Myeloproliferative neoplasms (MPNs) are cancers involving dysregulated production and function of myeloid lineage hematopoietic cells. Among MPNs, Essential thrombocythemia (ET), Polycythemia Vera (PV) and Myelofibrosis (MF), are driven by mutations that activate the JAK-STAT signalling pathway. Somatic mutations of calreticulin (CRT), an endoplasmic reticulum (ER)-localized lectin chaperone, are driver mutations in approximately 25% of ET and 35% of MF patients. The MPN-linked mutant CRT proteins have novel frameshifted carboxy-domain sequences and lack an ER retention motif, resulting in their secretion. Wild type CRT is a regulator of ER calcium homeostasis and plays a key role in the assembly of major histocompatibility complex (MHC) class I molecules, which are the ligands for antigen receptors of CD8+ T cells. Mutant CRT-linked oncogenesis results from the dysregulation of calcium signalling in cells and the formation of stable complexes of mutant CRT with myeloproliferative leukemia (MPL) protein, followed by downstream activation of the JAK-STAT signalling pathway. The intricate participation of CRT in ER protein folding, calcium homeostasis and immunity suggests the involvement of multiple mechanisms of mutant CRT-linked oncogenesis. In this review, we highlight recent findings related to the role of MPN-linked CRT mutations in the dysregulation of calcium homeostasis, MPL activation and immunity.

Structural Modeling of Cytokine-Receptor-JAK2 Signaling Complexes Using AlphaFold Multimer.

Homodimeric class 1 cytokine receptors include the erythropoietin (EPOR), thrombopoietin (TPOR), granulocyte colony-stimulating factor 3 (CSF3R), growth hormone (GHR), and prolactin receptors (PRLR). These cell-surface single-pass transmembrane (TM) glycoproteins regulate cell growth, proliferation, and differentiation and induce oncogenesis. An active TM signaling complex consists of a receptor homodimer, one or two ligands bound to the receptor extracellular domains, and two molecules of Janus Kinase 2 (JAK2) constitutively associated with the receptor intracellular domains. Although crystal structures of soluble extracellular domains with ligands have been obtained for all of the receptors except TPOR, little is known about the structure and dynamics of the complete TM complexes that activate the downstream JAK-STAT signaling pathway. Three-dimensional models of five human receptor complexes with cytokines and JAK2 were generated here by using AlphaFold Multimer. Given the large size of the complexes (from 3220 to 4074 residues), the modeling required a stepwise assembly from smaller parts, with selection and validation of the models through comparisons with published experimental data. The modeling of active and inactive complexes supports a general activation mechanism that involves ligand binding to a monomeric receptor followed by receptor dimerization and rotational movement of the receptor TM α-helices, causing proximity, dimerization, and activation of associated JAK2 subunits. The binding mode of two eltrombopag molecules to the TM α-helices of the active TPOR dimer was proposed. The models also help elucidate the molecular basis of oncogenic mutations that may involve a noncanonical activation route. Models equilibrated in explicit lipids of the plasma membrane are publicly available.

Comparative Molecular Dynamics Simulation Studies of Realistic Eukaryotic, Prokaryotic, and Archaeal Membranes.

We present a comparative all-atom molecular dynamics simulation study of 18 biomembrane systems with lipid compositions corresponding to eukaryotic, bacterial, and archaebacterial membranes together with three single-component lipid bilayers. A total of 105 lipid types used in this study include diverse sterols and glycerol-based lipids with acyl chains of various lengths, unsaturation degrees, and branched or cyclic moieties. Our comparative analysis provides deeper insight into the influences of sterols and lipid unsaturation on the structural and mechanical properties of these biomembranes, including water permeation into the membrane hydrocarbon core. For sterol-containing membranes, sterol fraction is correlated with the membrane thickness, the area compressibility modulus, and lipid order but anticorrelated with the area per lipid and sterol tilt angles. Similarly, for all 18 biomembranes, lipid order is correlated with the membrane thickness and area compressibility modulus. Sterols and lipid unsaturation produce opposite effects on membrane thickness, but only sterols influence water permeation into the membrane. All membrane systems are accessible for public use in CHARMM-GUI Archive. They can be used as templates to expedite future modeling of realistic cell membranes with transmembrane and peripheral membrane proteins to study their structure, dynamics, molecular interactions, and function in a nativelike membrane environment.

Thermodynamics-Based Molecular Modeling of α-Helices in Membranes and Micelles.

The Folding of Membrane-Associated Peptides (FMAP) method was developed for modeling α-helix formation by linear peptides in micelles and lipid bilayers. FMAP 2.0 identifies locations of α-helices in the amino acid sequence, generates their three-dimensional models in planar bilayers or spherical micelles, and estimates their thermodynamic stabilities and tilt angles, depending on temperature and pH. The method was tested for 723 peptides (926 data points) experimentally studied in different environments and for 170 single-pass transmembrane (TM) proteins with available crystal structures. FMAP 2.0 detected more than 95% of experimentally observed α-helices with an average error in helix end determination of around 2, 3, 4, and 5 residues per helix for peptides in water, micelles, bilayers, and TM proteins, respectively. Helical and nonhelical residue states were predicted with an accuracy from 0.86 to 0.96, and the Matthews correlation coefficient was from 0.64 to 0.88 depending on the environment. Experimental micelle- and membrane-binding energies and tilt angles of peptides were reproduced with a root-mean-square deviation of around 2 kcal/mol and 7°, respectively. The TM and non-TM states of hydrophobic and pH-triggered α-helical peptides in various lipid bilayers were reproduced in more than 95% of cases. The FMAP 2.0 web server (https://membranome.org/fmap) is publicly available to explore the structural polymorphism of antimicrobial, cell-penetrating, fusion, and other membrane-binding peptides, which is important for understanding the mechanisms of their biological activities.

Phospholipid Component Defines Pharmacokinetic and Pharmacodynamic Properties of Synthetic High-Density Lipoproteins.

Synthetic high-density lipoprotein (sHDL) nanoparticles composed of apolipoprotein A-I mimetic peptide and phospholipids have been shown to reduce atherosclerosis in animal models. Cholesterol is mobilized from atheroma macrophages by sHDL into the blood compartment and delivered to the liver for elimination. Historically, sHDL drug discovery efforts were focused on optimizing peptide sequences for interaction with cholesterol cellular transporters rather than understanding how both sHDL components, peptide and lipid, influence its pharmacokinetic and pharmacodynamic profiles. We designed two sets of sHDL having either identical phospholipid but variable peptide sequences with different plasma stability or identical peptide and phospholipids with variable fatty acid chain length and saturation. We found that sHDL prepared with proteolytically stable 22A-P peptide had 2-fold longer circulation half-time relative to the less stable 22A peptide. Yet, longer half-life did not translate into any improvement in cholesterol mobilization. In contrast, sHDL with variable phospholipid compositions showed significant differences in phospholipid PK, with distearoyl phosphatidylcholine-based sHDL demonstrating the longest half-life of 6.0 hours relative to 1.0 hour for palmitoyl-oleoyl phosphatidylcholine-based sHDL. This increase in half-life corresponded to an approx. 6.5-fold increase in the area under the curve for the mobilized cholesterol. Therefore, the phospholipid component in sHDL plays a major role in cholesterol mobilization in vivo and should not be overlooked in the design of future sHDL. SIGNIFICANCE STATEMENT: The phospholipid composition in sHDL plays a critical role in determining half-life and cholesterol mobilization in vivo.

Membranome 2.0: database for proteome-wide profiling of bitopic proteins and their dimers.

Motivation: Structural studies of TM domains of single-spanning (bitopic) membrane proteins are impeded by their instability, flexibility and heterogeneity. The new computational method TMDOCK allows reliable modeling of homodimers of transmembrane (TM) α-helices on a proteomic scale. Results: 3D models of 2129 parallel homodimers formed by TM α-helices of bitopic proteins from six evolutionarily distant organisms were modeled by TMDOCK, verified through experimental data available for nearly 600 proteins, and included in the Membranome database (v.2.0) along with related information to facilitate structural and evolutionary analysis of bitopic proteins. Availability and implementation: http://membranome.org. Contact: almz@umich.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Physics-Based Method for Modeling Passive Membrane Permeability and Translocation Pathways of Bioactive Molecules.

Assessment of permeability is a critical step in the drug development process for selection of drug candidates with favorable ADME properties. We have developed a novel physics-based method for fast computational modeling of passive permeation of diverse classes of molecules across lipid membranes. The method is based on heterogeneous solubility-diffusion theory and operates with all-atom 3D structures of solutes and the anisotropic solvent model of the lipid bilayer characterized by transbilayer profiles of dielectric and hydrogen bonding capacity parameters. The optimal translocation pathway of a solute is determined by moving an ensemble of representative conformations of the molecule through the dioleoyl-phosphatidylcholine (DOPC) bilayer and optimizing their rotational orientations in every point of the transmembrane trajectory. The method calculates (1) the membrane-bound state of the solute molecule; (2) free energy profile of the solute along the permeation pathway; and (3) the permeability coefficient obtained by integration over the transbilayer energy profile and assuming a constant size-dependent diffusivity along the membrane normal. The accuracy of the predictions was evaluated against experimental permeability coefficients measured in pure lipid membranes (for 78 compounds, R2 was 0.88 and rmse was 1.15 log units), PAMPA-DS (for 280 compounds, R2 was 0.75 and rmse was 1.59 log units), BBB (for 182 compounds, R2 was 0.69 and rmse was 0.87 log units), and Caco-2/MDCK assays (for 165 compounds, R2 was 0.52 and rmse was 0.89 log units).

PerMM: A Web Tool and Database for Analysis of Passive Membrane Permeability and Translocation Pathways of Bioactive Molecules.

The PerMM web server and database were developed for quantitative analysis and visualization of passive translocation of bioactive molecules across lipid membranes. The server is the first physics-based web tool that calculates membrane binding energies and permeability coefficients of diverse molecules through artificial and natural membranes (phospholipid bilayers, PAMPA-DS, blood-brain barrier, and Caco-2/MDCK cell membranes). It also visualizes the transmembrane translocation pathway as a sequence of translational and rotational positions of a permeant as it moves across the lipid bilayer, along with the corresponding changes in solvation energy. The server can be applied for prediction of permeability coefficients of compounds with diverse chemical scaffolds to facilitate selection and optimization of potential drug leads. The complementary PerMM database allows comparison of computationally and experimentally determined permeability coefficients for more than 500 compounds in different membrane systems. The website and database are freely accessible at https://permm.phar.umich.edu/ .

Membranome: a database for proteome-wide analysis of single-pass membrane proteins.

The Membranome database was developed to assist analysis and computational modeling of single-pass (bitopic) transmembrane (TM) proteins and their complexes by providing structural information about these proteins on a genomic scale. The database currently collects data on >6000 bitopic proteins from Homo sapiens, Arabidopsis thaliana, Dictyostelium discoideum, Saccharomyces cerevisiae, Escherichia coli and Methanocaldococcus jannaschii It presents the following data: (i) hierarchical classification of bitopic proteins into 15 functional classes, 689 structural superfamilies and 1404 families; (ii) 446 complexes of bitopic proteins with known three-dimensional (3D) structures classified into 129 families; (iii) computationally generated three-dimensional models of TM α-helices positioned in membranes; (iv) amino acid sequences, domain architecture, functional annotation and available experimental structures of bitopic proteins; (v) TM topology and intracellular localization, (vi) physical interactions between proteins from the database along with links to other resources. The database is freely accessible at http://membranome.org There is a variety of options for browsing, sorting, searching and retrieval of the content, including downloadable coordinate files of TM domains with calculated membrane boundaries.

Evolution and adaptation of single-pass transmembrane proteins.

A comparative analysis of 6039 single-pass (bitopic) membrane proteins from six evolutionarily distant organisms was performed based on data from the Membranome database. The observed repertoire of bitopic proteins is significantly enlarged in eukaryotic cells and especially in multicellular organisms due to the diversification of enzymes, emergence of proteins involved in vesicular trafficking, and expansion of receptors, structural, and adhesion proteins. The majority of bitopic proteins in multicellular organisms are located in the plasma membrane (PM) and involved in cell communication. Bitopic proteins from different membranes significantly diverge in terms of their biological functions, size, topology, domain architecture, physical properties of transmembrane (TM) helices and propensity to form homodimers. Most proteins from eukaryotic PM and endoplasmic reticulum (ER) have the N-out topology. The predicted lengths of TM helices and hydrophobic thicknesses, stabilities and hydrophobicities of TM α-helices are the highest for proteins from eukaryotic PM, intermediate for proteins from prokaryotic cells, ER and Golgi apparatus, and lowest for proteins from mitochondria, chloroplasts, and peroxisomes. Tyr and Phe residues accumulate at the cytoplasmic leaflet of PM and at the outer leaflet of membranes of bacteria, Golgi apparatus, and nucleus. The propensity for dimerization increases from unicellular to multicellular eukaryotes, from enzymes to receptors, and from intracellular membrane proteins to PM proteins. More than half of PM proteins form homodimers with a 2:1 ratio of right-handed to left-handed helix packing arrangements. The inverse ratio (1:2) was observed for dimers from the ER, Golgi and vesicles.

Effect of Site-Specific Intermolecular Lysine-Tryptophan Interactions on the Aggregation of Gramicidin-Based Peptides Leading to Pore Formation in Lipid Membranes.

In contrast to the parent pentadecapeptide gramicidin A (gA), some of its cationic analogs have been shown previously to form large-diameter pores in lipid membranes. These pores are permeable to fluorescent dyes, which allows one to monitor pore formation by using the fluorescence de-quenching assay. According to the previously proposed model, the gA analog with lysine substituted for alanine at position 3, [Lys3]gA, forms pores by a homopentameric assembly of gramicidin double-stranded ß-helical dimers. Here, we studied the newly synthesized analogs of [Lys3]gA with single, double and triple substitutions of isoleucines for tryptophans at positions 9, 11, 13, and 15. Replacement of any of the tryptophans of [Lys3]gA with isoleucine resulted in suppression of the pore-forming activity of the peptide, the effect being significantly dependent on the position of tryptophans. In particular, the peptide with a single substitution of tryptophan 13 showed much lower activity than the analogs with single substitutions at positions 9, 11, or 15. Of the peptides with double substitutions, the strongest suppression of the leakage was observed with tryptophans 13 and 15. In the case of triple substitutions, only the peptide retaining tryptophan 11 exhibited noticeable activity. It is concluded that tryptophans 11 and 13 contribute most to pore stabilization in the membrane, whereas tryptophan 9 is not so important for pore formation. Cation-π interactions between the lysine and tryptophan residues of the peptide are suggested to be crucial for the formation of the [Lys3]gA pore.

Use of Functional Polymorphisms To Elucidate the Peptide Binding Site of TAP Complexes.

TAP1/TAP2 complexes translocate peptides from the cytosol to the endoplasmic reticulum lumen to enable immune surveillance by CD8(+) T cells. Peptide transport is preceded by peptide binding to a cytosol-accessible surface of TAP1/TAP2 complexes, but the location of the TAP peptide-binding pocket remains unknown. Guided by the known contributions of polymorphic TAP variants to peptide selection, we combined homology modeling of TAP with experimental measurements to identify several TAP residues that interact with peptides. Models for peptide-TAP complexes were generated, which indicate bent conformation for peptides. The peptide binding site of TAP is located at the hydrophobic boundary of the cytosolic membrane leaflet, with striking parallels to the glutathione binding site of NaAtm1, a transporter that functions in bacterial heavy metal detoxification. These studies illustrate the conservation of the ligand recognition modes of bacterial and mammalians transporters involved in peptide-guided cellular surveillance.

Gramicidin A disassembles large conductive clusters of its lysine-substituted derivatives in lipid membranes.

N-terminally substituted lysine derivatives of gramicidin A (gA), [Lys1]gA and [Lys3]gA, but not glutamate- or aspartate-substituted peptides have been previously shown to cause the leakage of carboxyfluorescein from liposomes. Here, the leakage induction was also observed for [Arg1]gA and [Arg3]gA, while [His1]gA and [His3]gA were inactive at neutral pH. The Lys3-containing analogue with all tryptophans replaced by isoleucines did not induce liposome leakage, similar to gA. This suggests that the presence of both tryptophans and N-terminal cationic residues is critical for pore formation. Remarkably, the addition of gA blocked the leakage induced by [Lys3]gA. By examining with fluorescence correlation spectroscopy the peptide-induced leakage of fluorescent markers from liposomes, we estimated the diameter of pores responsible for the leakage to be about 1.6 nm. Transmission electron cryo-microscopy imaging of liposomes with [Lys3]gA showed that the liposomal membranes contained high electron density particles with a size of about 40 Å, suggesting the formation of peptide clusters. No such clusterization was observed in liposomes incorporating gA or a mixture of gA with [Lys3]gA. Three-dimensional reconstruction of the clusters was compatible with their pentameric arrangement. Based on experimental data and computational modeling, we suggest that the large pore formed by [Lys3]gA represents a barrel-stave oligomeric cluster formed by antiparallel double-stranded helical dimers (DH). In a tentative model, the pentamer of dimers may be stabilized by aromatic Trp-Trp and cation-π Trp-Lys interactions between the neighboring DHs. The inhibiting effect of gA on the [Lys3]gA-induced leakage can be attributed to breaking of cation-π interactions, which prevents peptide clusterization and pore formation.

Structural adaptations of proteins to different biological membranes.

To gain insight into adaptations of proteins to their membranes, intrinsic hydrophobic thicknesses, distributions of different chemical groups and profiles of hydrogen-bonding capacities (α and ß) and the dipolarity/polarizability parameter (π*) were calculated for lipid-facing surfaces of 460 integral α-helical, ß-barrel and peripheral proteins from eight types of biomembranes. For comparison, polarity profiles were also calculated for ten artificial lipid bilayers that have been previously studied by neutron and X-ray scattering. Estimated hydrophobic thicknesses are 30-31Å for proteins from endoplasmic reticulum, thylakoid, and various bacterial plasma membranes, but differ for proteins from outer bacterial, inner mitochondrial and eukaryotic plasma membranes (23.9, 28.6 and 33.5Å, respectively). Protein and lipid polarity parameters abruptly change in the lipid carbonyl zone that matches the calculated hydrophobic boundaries. Maxima of positively charged protein groups correspond to the location of lipid phosphates at 20-22Å distances from the membrane center. Locations of Tyr atoms coincide with hydrophobic boundaries, while distributions maxima of Trp rings are shifted by 3-4Å toward the membrane center. Distributions of Trp atoms indicate the presence of two 5-8Å-wide midpolar regions with intermediate π* values within the hydrocarbon core, whose size and symmetry depend on the lipid composition of membrane leaflets. Midpolar regions are especially asymmetric in outer bacterial membranes and cell membranes of mesophilic but not hyperthermophilic archaebacteria, indicating the larger width of the central nonpolar region in the later case. In artificial lipid bilayers, midpolar regions are observed up to the level of acyl chain double bonds.

Translation of structure-activity relationships from cyclic mixed efficacy opioid peptides to linear analogues.

Most opioid analgesics used in the treatment of pain are mu opioid receptor (MOR) agonists. While effective, there are significant drawbacks to opioid use, including the development of tolerance and dependence. However, the coadministration of a MOR agonist with a delta opioid receptor (DOR) antagonist slows the development of MOR-related side effects, while maintaining analgesia. We have previously reported a series of cyclic mixed efficacy MOR agonist/DOR antagonist ligands. Here we describe the transfer of key features from these cyclic analogs to linear sequences. Using the linear MOR/DOR agonist, Tyr-DThr-Gly-Phe-Leu-Ser-NH2 (DTLES), as a lead scaffold, we replaced Phe(4) with bulkier and/or constrained aromatic residues shown to confer DOR antagonism in our cyclic ligands. These replacements failed to confer DOR antagonism in the DTLES analogs, presumably because the more flexible linear ligands can adopt binding poses that will fit in the narrow binding pocket of the active conformations of both MOR and DOR. Nonetheless, the pharmacological profile observed in this series, high affinity and efficacy for MOR and DOR with selectivity relative to KOR, has also been shown to reduce the development of unwanted side effects. We further modified our lead MOR/DOR agonist with a C-terminal glucoserine to improve bioavailability. The resulting ligand displayed high efficacy and potency at both MOR and DOR and no efficacy at KOR.

OPM database and PPM web server: resources for positioning of proteins in membranes.

The Orientations of Proteins in Membranes (OPM) database is a curated web resource that provides spatial positions of membrane-bound peptides and proteins of known three-dimensional structure in the lipid bilayer, together with their structural classification, topology and intracellular localization. OPM currently contains more than 1200 transmembrane and peripheral proteins and peptides from approximately 350 organisms that represent approximately 3800 Protein Data Bank entries. Proteins are classified into classes, superfamilies and families and assigned to 21 distinct membrane types. Spatial positions of proteins with respect to the lipid bilayer are optimized by the PPM 2.0 method that accounts for the hydrophobic, hydrogen bonding and electrostatic interactions of the proteins with the anisotropic water-lipid environment described by the dielectric constant and hydrogen-bonding profiles. The OPM database is freely accessible at http://opm.phar.umich.edu. Data can be sorted, searched or retrieved using the hierarchical classification, source organism, localization in different types of membranes. The database offers downloadable coordinates of proteins and peptides with membrane boundaries. A gallery of protein images and several visualization tools are provided. The database is supplemented by the PPM server (http://opm.phar.umich.edu/server.php) which can be used for calculating spatial positions in membranes of newly determined proteins structures or theoretical models.

Structural modeling of cytokine-receptor-JAK2 signaling complexes using AlphaFold Multimer.

Homodimeric class 1 cytokine receptors include the erythropoietin (EPOR), thrombopoietin (TPOR), granulocyte colony-stimulating factor 3 (CSF3R), growth hormone (GHR), and prolactin receptors (PRLR). They are cell-surface single-pass transmembrane (TM) glycoproteins that regulate cell growth, proliferation, and differentiation and induce oncogenesis. An active TM signaling complex consists of a receptor homodimer, one or two ligands bound to the receptor extracellular domains and two molecules of Janus Kinase 2 (JAK2) constitutively associated with the receptor intracellular domains. Although crystal structures of soluble extracellular domains with ligands have been obtained for all the receptors except TPOR, little is known about the structure and dynamics of the complete TM complexes that activate the downstream JAK-STAT signaling pathway. Three-dimensional models of five human receptor complexes with cytokines and JAK2 were generated using AlphaFold Multimer. Given the large size of the complexes (from 3220 to 4074 residues), the modeling required a stepwise assembly from smaller parts with selection and validation of the models through comparisons with published experimental data. The modeling of active and inactive complexes supports a general activation mechanism that involves ligand binding to a monomeric receptor followed by receptor dimerization and rotational movement of the receptor TM α-helices causing proximity, dimerization, and activation of associated JAK2 subunits. The binding mode of two eltrombopag molecules to TM α-helices of the active TPOR dimer was proposed. The models also help elucidating the molecular basis of oncogenic mutations that may involve non-canonical activation route. Models equilibrated in explicit lipids of the plasma membrane are publicly available.

Spatial arrangement of proteins in planar and curved membranes by PPM 3.0.

Cellular protrusions, invaginations, and many intracellular organelles have strongly curved membrane regions. Transmembrane and peripheral membrane proteins that induce, sense, or stabilize such regions cannot be properly fitted into a single flat bilayer. To treat such proteins, we developed a new method and a web tool, PPM 3.0, for positioning proteins in curved or planar, single or multiple membranes. This method determines the energetically optimal spatial position, the hydrophobic thickness, and the radius of intrinsic curvature of a membrane-deforming protein structure by arranging it in a single or several sphere-shaped or planar membrane sections. In addition, it can define the lipid-embedded regions of a protein that simultaneously spans several membranes or determine the optimal position of a peptide in a spherical micelle. The PPM 3.0 web server operates with 17 types of biological membranes and 4 types of artificial bilayers. It is publicly available at https://opm.phar.umich.edu/ppm_server3. PPM 3.0 was applied to identify and characterize arrangements in membranes of 128 proteins with a significant intrinsic curvature, such as BAR domains, annexins, Piezo, and MscS mechanosensitive channels, cation-chloride cotransporters, as well as mitochondrial ATP synthases, calcium uniporters, and TOM complexes. These proteins form large complexes that are mainly localized in mitochondria, plasma membranes, and endosomes. Structures of bacterial drug efflux pumps, AcrAB-TolC, MexAB-OrpM, and MacAB-TolC, were positioned in both membranes of the bacterial cell envelop, while structures of multimeric gap-junction channels were arranged in two opposed cellular membranes.

Membranome 3.0: Database of single-pass membrane proteins with AlphaFold models.

The Membranome database provides comprehensive structural information on single-pass (i.e., bitopic) membrane proteins from six evolutionarily distant organisms, including protein-protein interactions, complexes, mutations, experimental structures, and models of transmembrane α-helical dimers. We present a new version of this database, Membranome 3.0, which was significantly updated by revising the set of 5,758 bitopic proteins and incorporating models generated by AlphaFold 2 in the database. The AlphaFold models were parsed into structural domains located at the different membrane sides, modified to exclude low-confidence unstructured terminal regions and signal sequences, validated through comparison with available experimental structures, and positioned with respect to membrane boundaries. Membranome 3.0 was re-developed to facilitate visualization and comparative analysis of multiple 3D structures of proteins that belong to a specified family, complex, biological pathway, or membrane type. New tools for advanced search and analysis of proteins, their interactions, complexes, and mutations were included. The database is freely accessible at https://membranome.org.

Membrane topology of the colicin E1 channel using genetically encoded fluorescence.

The membrane topology of the colicin E1 channel domain was studied by fluorescence resonance energy transfer (FRET). The FRET involved a genetically encoded fluorescent amino acid (coumarin) as the donor and a selectively labeled cysteine residue tethered with DABMI (4-(dimethylamino)phenylazophenyl-4'-maleimide) as the FRET acceptor. The fluorescent coumarin residue was incorporated into the protein via an orthogonal tRNA/aminoacyl-tRNA synthetase pair that allowed selective incorporation into any site within the colicin channel domain. Each variant harbored a stop (TAG) mutation for coumarin incorporation and a cysteine (TGT) mutation for DABMI attachment. Six interhelical distances within helices 1-6 were determined using FRET analysis for both the soluble and membrane-bound states. The FRET data showed large changes in the interhelical distances among helices 3-6 upon membrane association providing new insight into the membrane-bound structure of the channel domain. In general, the coumarin-DABMI FRET interhelical efficiencies decreased upon membrane binding, building upon the umbrella model for the colicin channel. A tentative model for the closed state of the channel domain was developed based on current and previously published FRET data. The model suggests circular arrangement of helices 1-7 in a clockwise direction from the extracellular side and membrane interfacial association of helices 1, 6, 7, and 10 around the central transmembrane hairpin formed by helices 8 and 9.