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.

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.

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.

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.

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.

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.

TMPfold: A Web Tool for Predicting Stability of Transmembrane α-Helix Association.

Estimating energies of transmembrane (TM) α-helix association is essential for understanding folding of membrane proteins and formation of their functional assemblies. A new physics-based method was developed and implemented in the TMPfold web server for the calculation of the free energy of TM helix association (ΔGasc) in TM α-bundles of known structure. The method was verified using the experimental ΔGasc values for 36 TM complexes, including dimers of 10 glycophorin A mutants. The calculated free energy changes (ΔΔGasc) caused by mutations in TM helices correlated with experimental changes in the stability of 42 mutants of bacteriorhodopsin and 25 mutants of rhomboid protease. TMPfold was applied for evaluation of ΔGasc in 554 PDB structures of 85 seven-helical TM proteins and identification of stable two-helical folding intermediates. The proposed tentative paths of cotranslational helix assembly of several polytopic proteins were consistent with experimental studies of their folding. TMPfold is accessible at (https://opm.phar.umich.edu/tmpfold_server).

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/ .

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).

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.

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.

TMDOCK: An Energy-Based Method for Modeling α-Helical Dimers in Membranes.

TMDOCK is a novel computational method for the modeling of parallel homodimers formed by transmembrane (TM) α-helices. Three-dimensional (3D) models of dimers are generated by threading a target amino acid sequence through several structural templates, followed by local energy minimization. This is the first method that identifies helix dimerization modes and ranks them based on the calculated free energy of α-helix association. Free energy components include van der Waals, hydrogen bonding, and dipole interactions; side-chain conformational entropy; and solvation energy in the anisotropic lipid environment. TMDOCK reproduced 26 experimental dimeric structures formed by TM α-helices of 21 single-pass membrane proteins (including 4 mutants) with Cα atom rmsd from 1.0 to 3.3Å. Assessment of dimerization heterogeneity of these TM domains demonstrated that 7 of them have a unique dimer structure, 12 have at least 2 alternative conformations, and 2 have a large number of different association modes. All unique experimental structures of proteins from the first group and eight structures from the second group were reproduced in computations as top-ranked models. A fast version of the method is available through the web server (http://membranome.org/tm_server.php).

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.

Life at the border: adaptation of proteins to anisotropic membrane environment.

This review discusses main features of transmembrane (TM) proteins which distinguish them from water-soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co-translational and post-translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large-scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen-bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region-between double bonds and carbonyl groups of lipids. These "midpolar" regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side-chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein-lipid binding.

Solvation models and computational prediction of orientations of peptides and proteins in membranes.

Membrane-associated peptides and proteins function in the highly heterogeneous environment of the lipid bilayer whose physico-chemical properties change non-monotonically along the bilayer normal. To simulate insertion of peptides and proteins into membranes and correctly reproduce the energetics of this process, an adequate solvation model and physically realistic representation of the lipid bilayer should be employed. We present a brief overview of the existing solvation models and their application for prediction of binding affinities and orientations of proteins in membranes. Particular emphasis is placed on the recently proposed PPM method, the corresponding web server, and the OPM database that were designed for positioning in membranes of integral and peripheral proteins with known three-dimensional structures.

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.

Antimicrobial action of the cyclic peptide bactenecin on Burkholderia pseudomallei correlates with efficient membrane permeabilization.

Burkholderia pseudomallei is a category B agent that causes Melioidosis, an acute and chronic disease with septicemia. The current treatment regimen is a heavy dose of antibiotics such as ceftazidime (CAZ); however, the risk of a relapse is possible. Peptide antibiotics are an alternative to classical antibiotics as they exhibit rapid action and are less likely to result in the development of resistance. The aim of this study was to determine the bactericidal activity against B. pseudomallei and examine the membrane disrupting abilities of the potent antimicrobial peptides: bactenecin, RTA3, BMAP-18 and CA-MA. All peptides exhibited >97% bactericidal activity at 20 µM, with bactenecin having slightly higher activity. Long term time-kill assays revealed a complete inhibition of cell growth at 50 µM bactenecin and CA-MA. All peptides inhibited biofilm formation comparable to CAZ, but exhibited faster kinetics (within 1 h). Bactenecin exhibited stronger binding to LPS and induced perturbation of the inner membrane of live cells. Interaction of bactenecin with model membranes resulted in changes in membrane fluidity and permeability, leading to leakage of dye across the membrane at levels two-fold greater than that of other peptides. Modeling of peptide binding on the membrane showed stable and deep insertion of bactenecin into the membrane (up to 9 Å). We propose that bactenecin is able to form dimers or large ß-sheet structures in a concentration dependent manner and subsequently rapidly permeabilize the membrane, leading to cytosolic leakage and cell death in a shorter period of time compared to CAZ. Bactenecin might be considered as a potent antimicrobial agent for use against B. pseudomallei.

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.

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.

Cytotoxic potency of small macrocyclic knot proteins: structure-activity and mechanistic studies of native and chemically modified cyclotides.

The cyclotides are a family of circular and knotted proteins of natural origin with extreme enzymatic and thermal stability. They have a wide range of biological activities that make them promising tools for pharmaceutical and crop-protection applications. The cyclotides are divided into two subfamilies depending on the presence (Möbius) or absence (bracelet) of a cis-Pro peptide bond. In the current work we report a series of experiments to give further insight into the structure-activity relationship of cyclotides in general, and the differences between subfamilies and the role of their hydrophobic surface in particular. Selective chemical modifications of Glu, Arg, Lys and Trp residues was tested for cytotoxic activity: derivatives in which the Trp residue was modified showed low effect, demonstrating the existence of a connection between hydrophobicity and activity. However, over the full set of cyclotides examined, there was no strong correlation between the cytotoxic activity and their hydrophobicity. Instead, it seems more like that the distribution of charged and hydrophobic residues determines the ultimate degree of potency. Furthermore, we found that while the Glu residue is very important in maintaining the activity of the bracelet cyclotide cycloviolacin O2, it is much less important in the Möbius cyclotides. Despite these differences between cyclotide subfamilies, a systematic test of mixtures of cyclotides revealed that they act in an additive way.