Chemically induced vesiculation as a platform for studying TMEM16F activity.

Calcium-activated phospholipid scramblase mediates the energy-independent bidirectional translocation of lipids across the bilayer, leading to transient or, in the case of apoptotic scrambling, sustained collapse of membrane asymmetry. Cells lacking TMEM16F-dependent lipid scrambling activity are deficient in generation of extracellular vesicles (EVs) that shed from the plasma membrane in a Ca2+-dependent manner, namely microvesicles. We have adapted chemical induction of giant plasma membrane vesicles (GPMVs), which require both TMEM16F-dependent phospholipid scrambling and calcium influx, as a kinetic assay to investigate the mechanism of TMEM16F activity. Using the GPMV assay, we identify and characterize both inactivating and activating mutants that elucidate the mechanism for TMEM16F activation and facilitate further investigation of TMEM16F-mediated lipid translocation and its role in extracellular vesiculation.

Atomistic insight into lipid translocation by a TMEM16 scramblase.

The transmembrane protein 16 (TMEM16) family of membrane proteins includes both lipid scramblases and ion channels involved in olfaction, nociception, and blood coagulation. The crystal structure of the fungal Nectria haematococca TMEM16 (nhTMEM16) scramblase suggested a putative mechanism of lipid transport, whereby polar and charged lipid headgroups move through the low-dielectric environment of the membrane by traversing a hydrophilic groove on the membrane-spanning surface of the protein. Here, we use computational methods to explore the membrane-protein interactions involved in lipid scrambling. Fast, continuum membrane-bending calculations reveal a global pattern of charged and hydrophobic surface residues that bends the membrane in a large-amplitude sinusoidal wave, resulting in bilayer thinning across the hydrophilic groove. Atomic simulations uncover two lipid headgroup-interaction sites flanking the groove. The cytoplasmic site nucleates headgroup-dipole stacking interactions that form a chain of lipid molecules that penetrate into the groove. In two instances, a cytoplasmic lipid interdigitates into this chain, crosses the bilayer, and enters the extracellular leaflet, and the reverse process happens twice as well. Continuum membrane-bending analysis carried out on homology models of mammalian homologs shows that these family members also bend the membrane-even those that lack scramblase activity. Sequence alignments show that the lipid-interaction sites are conserved in many family members but less so in those with reduced scrambling ability. Our analysis provides insight into how large-scale membrane bending and protein chemistry facilitate lipid permeation in the TMEM16 family, and we hypothesize that membrane interactions also affect ion permeation.

New Continuum Approaches for Determining Protein-Induced Membrane Deformations.

The influence of the membrane on transmembrane proteins is central to a number of biological phenomena, notably the gating of stretch activated ion channels. Conversely, membrane proteins can influence the bilayer, leading to the stabilization of particular membrane shapes, topological changes that occur during vesicle fission and fusion, and shape-dependent protein aggregation. Continuum elastic models of the membrane have been widely used to study protein-membrane interactions. These mathematical approaches produce physically interpretable membrane shapes, energy estimates for the cost of deformation, and a snapshot of the equilibrium configuration. Moreover, elastic models are much less computationally demanding than fully atomistic and coarse-grained simulation methodologies; however, it has been argued that continuum models cannot reproduce the distortions observed in fully atomistic molecular dynamics simulations. We suggest that this failure can be overcome by using chemically and geometrically accurate representations of the protein. Here, we present a fast and reliable hybrid continuum-atomistic model that couples the protein to the membrane. We show that the model is in excellent agreement with fully atomistic simulations of the ion channel gramicidin embedded in a POPC membrane. Our continuum calculations not only reproduce the membrane distortions produced by the channel but also accurately determine the channel's orientation. Finally, we use our method to investigate the role of membrane bending around the charged voltage sensors of the transient receptor potential cation channel TRPV1. We find that membrane deformation significantly stabilizes the energy of insertion of TRPV1 by exposing charged residues on the S4 segment to solution.

Continuum descriptions of membranes and their interaction with proteins: Towards chemically accurate models.

Biological membranes deform in response to resident proteins leading to a coupling between membrane shape and protein localization. Additionally, the membrane influences the function of membrane proteins. Here we review contributions to this field from continuum elastic membrane models focusing on the class of models that couple the protein to the membrane. While it has been argued that continuum models cannot reproduce the distortions observed in fully-atomistic molecular dynamics simulations, we suggest that this failure can be overcome by using chemically accurate representations of the protein. We outline our recent advances along these lines with our hybrid continuum-atomistic model, and we show the model is in excellent agreement with fully-atomistic simulations of the nhTMEM16 lipid scramblase. We believe that the speed and accuracy of continuum-atomistic methodologies will make it possible to simulate large scale, slow biological processes, such as membrane morphological changes, that are currently beyond the scope of other computational approaches. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

De novo design of pH-responsive self-assembling helical protein filaments.

Biological evolution has led to precise and dynamic nanostructures that reconfigure in response to pH and other environmental conditions. However, designing micrometre-scale protein nanostructures that are environmentally responsive remains a challenge. Here we describe the de novo design of pH-responsive protein filaments built from subunits containing six or nine buried histidine residues that assemble into micrometre-scale, well-ordered fibres at neutral pH. The cryogenic electron microscopy structure of an optimized design is nearly identical to the computational design model for both the subunit internal geometry and the subunit packing into the fibre. Electron, fluorescent and atomic force microscopy characterization reveal a sharp and reversible transition from assembled to disassembled fibres over 0.3 pH units, and rapid fibre disassembly in less than 1 s following a drop in pH. The midpoint of the transition can be tuned by modulating buried histidine-containing hydrogen bond networks. Computational protein design thus provides a route to creating unbound nanomaterials that rapidly respond to small pH changes.

Precisely patterned nanofibres made from extendable protein multiplexes.

Molecular systems with coincident cyclic and superhelical symmetry axes have considerable advantages for materials design as they can be readily lengthened or shortened by changing the length of the constituent monomers. Among proteins, alpha-helical coiled coils have such symmetric, extendable architectures, but are limited by the relatively fixed geometry and flexibility of the helical protomers. Here we describe a systematic approach to generating modular and rigid repeat protein oligomers with coincident C2 to C8 and superhelical symmetry axes that can be readily extended by repeat propagation. From these building blocks, we demonstrate that a wide range of unbounded fibres can be systematically designed by introducing hydrophilic surface patches that force staggering of the monomers; the geometry of such fibres can be precisely tuned by varying the number of repeat units in the monomer and the placement of the hydrophilic patches.

Local structural flexibility drives oligomorphism in computationally designed protein assemblies.

Many naturally occurring protein assemblies have dynamic structures that allow them to perform specialized functions. For example, clathrin coats adopt a wide variety of architectures to adapt to vesicular cargos of various sizes. Although computational methods for designing novel self-assembling proteins have advanced substantially over the past decade, most existing methods focus on designing static structures with high accuracy. Here we characterize the structures of three distinct computationally designed protein assemblies that each form multiple unanticipated architectures, and identify flexibility in specific regions of the subunits of each assembly as the source of structural diversity. Cryo-EM single-particle reconstructions and native mass spectrometry showed that only two distinct architectures were observed in two of the three cases, while we obtained six cryo-EM reconstructions that likely represent a subset of the architectures present in solution in the third case. Structural modeling and molecular dynamics simulations indicated that the surprising observation of a defined range of architectures, instead of non-specific aggregation, can be explained by constrained flexibility within the building blocks. Our results suggest that deliberate use of structural flexibility as a design principle will allow exploration of previously inaccessible structural and functional space in designed protein assemblies.

Natural and Designed Proteins Inspired by Extremotolerant Organisms Can Form Condensates and Attenuate Apoptosis in Human Cells.

Many organisms can survive extreme conditions and successfully recover to normal life. This extremotolerant behavior has been attributed in part to repetitive, amphipathic, and intrinsically disordered proteins that are upregulated in the protected state. Here, we assemble a library of approximately 300 naturally occurring and designed extremotolerance-associated proteins to assess their ability to protect human cells from chemically induced apoptosis. We show that several proteins from tardigrades, nematodes, and the Chinese giant salamander are apoptosis-protective. Notably, we identify a region of the human ApoE protein with similarity to extremotolerance-associated proteins that also protects against apoptosis. This region mirrors the phase separation behavior seen with such proteins, like the tardigrade protein CAHS2. Moreover, we identify a synthetic protein, DHR81, that shares this combination of elevated phase separation propensity and apoptosis protection. Finally, we demonstrate that driving protective proteins into the condensate state increases apoptosis protection, and highlights the ability of DHR81 condensates to sequester caspase-7. Taken together, this work draws a link between extremotolerance-associated proteins, condensate formation, and designing human cellular protection.

A multiscale model of mechanotransduction by the ankyrin chains of the NOMPC channel.

Our senses of touch and hearing are dependent on the conversion of external mechanical forces into electrical impulses by the opening of mechanosensitive channels in sensory cells. This remarkable feat involves the conversion of a macroscopic mechanical displacement into a subnanoscopic conformational change within the ion channel. The mechanosensitive channel NOMPC, responsible for hearing and touch in flies, is a homotetramer composed of four pore-forming transmembrane domains and four helical chains of 29 ankyrin repeats that extend 150 Å into the cytoplasm. Previous work has shown that the ankyrin chains behave as biological springs under extension and that tethering them to microtubules could be involved in the transmission of external forces to the NOMPC gate. Here we combine normal mode analysis (NMA), full-atom molecular dynamics simulations, and continuum mechanics to characterize the material properties of the chains under extreme compression and extension. NMA reveals that the lowest-frequency modes of motion correspond to fourfold symmetric compression/extension along the channel, and the lowest-frequency symmetric mode for the isolated channel domain involves rotations of the TRP domain, a putative gating element. Finite element modeling reveals that the ankyrin chains behave as a soft spring with a linear, effective spring constantof 22 pN/nm for deflections ≤15 Å. Force-balance analysis shows that the entire channel undergoes rigid body rotation during compression, and more importantly, each chain exerts a positive twisting moment on its respective linker helices and TRP domain. This torque is a model-independent consequence of the bundle geometry and would cause a clockwise rotation of the TRP domain when viewed from the cytoplasm. Force transmission to the channel for compressions >15 Å depends on the nature of helix-helix contact. Our work reveals that compression of the ankyrin chains imparts a rotational torque on the TRP domain, which potentially results in channel opening.

The Sixth Transmembrane Segment Is a Major Gating Component of the TMEM16A Calcium-Activated Chloride Channel.

Calcium-activated chloride channels (CaCCs) formed by TMEM16A or TMEM16B are broadly expressed in the nervous system, smooth muscles, exocrine glands, and other tissues. With two calcium-binding sites and a pore within each monomer, the dimeric CaCC exhibits voltage-dependent calcium sensitivity. Channel activity also depends on the identity of permeant anions. To understand how CaCC regulates neuronal signaling and how CaCC is, in turn, modulated by neuronal activity, we examined the molecular basis of CaCC gating. Here, we report that voltage modulation of TMEM16A-CaCC involves voltage-dependent occupancy of calcium- and anion-binding site(s) within the membrane electric field as well as a voltage-dependent conformational change intrinsic to the channel protein. These gating modalities all critically depend on the sixth transmembrane segment.

Transmembrane helix hydrophobicity is an energetic barrier during the retrotranslocation of integral membrane ERAD substrates.

Integral membrane proteins fold inefficiently and are susceptible to turnover via the endoplasmic reticulum-associated degradation (ERAD) pathway. During ERAD, misfolded proteins are recognized by molecular chaperones, polyubiquitinated, and retrotranslocated to the cytoplasm for proteasomal degradation. Although many aspects of this pathway are defined, how transmembrane helices (TMHs) are removed from the membrane and into the cytoplasm before degradation is poorly understood. In this study, we asked whether the hydrophobic character of a TMH acts as an energetic barrier to retrotranslocation. To this end, we designed a dual-pass model ERAD substrate, Chimera A*, which contains the cytoplasmic misfolded domain from a characterized ERAD substrate, Sterile 6* (Ste6p*). We found that the degradation requirements for Chimera A* and Ste6p* are similar, but Chimera A* was retrotranslocated more efficiently than Ste6p* in an in vitro assay in which retrotranslocation can be quantified. We then constructed a series of Chimera A* variants containing synthetic TMHs with a range of ΔG values for membrane insertion. TMH hydrophobicity correlated inversely with retrotranslocation efficiency, and in all cases, retrotranslocation remained Cdc48p dependent. These findings provide insight into the energetic restrictions on the retrotranslocation reaction, as well as a new computational approach to predict retrotranslocation efficiency.