Modulation of FGF pathway signaling and vascular differentiation using designed oligomeric assemblies.

Many growth factors and cytokines signal by binding to the extracellular domains of their receptors and driving association and transphosphorylation of the receptor intracellular tyrosine kinase domains, initiating downstream signaling cascades. To enable systematic exploration of how receptor valency and geometry affect signaling outcomes, we designed cyclic homo-oligomers with up to 8 subunits using repeat protein building blocks that can be modularly extended. By incorporating a de novo-designed fibroblast growth factor receptor (FGFR)-binding module into these scaffolds, we generated a series of synthetic signaling ligands that exhibit potent valency- and geometry-dependent Ca2+ release and mitogen-activated protein kinase (MAPK) pathway activation. The high specificity of the designed agonists reveals distinct roles for two FGFR splice variants in driving arterial endothelium and perivascular cell fates during early vascular development. Our designed modular assemblies should be broadly useful for unraveling the complexities of signaling in key developmental transitions and for developing future therapeutic applications.

Blueprinting extendable nanomaterials with standardized protein blocks.

A wooden house frame consists of many different lumber pieces, but because of the regularity of these building blocks, the structure can be designed using straightforward geometrical principles. The design of multicomponent protein assemblies, in comparison, has been much more complex, largely owing to the irregular shapes of protein structures1. Here we describe extendable linear, curved and angled protein building blocks, as well as inter-block interactions, that conform to specified geometric standards; assemblies designed using these blocks inherit their extendability and regular interaction surfaces, enabling them to be expanded or contracted by varying the number of modules, and reinforced with secondary struts. Using X-ray crystallography and electron microscopy, we validate nanomaterial designs ranging from simple polygonal and circular oligomers that can be concentrically nested, up to large polyhedral nanocages and unbounded straight 'train track' assemblies with reconfigurable sizes and geometries that can be readily blueprinted. Because of the complexity of protein structures and sequence-structure relationships, it has not previously been possible to build up large protein assemblies by deliberate placement of protein backbones onto a blank three-dimensional canvas; the simplicity and geometric regularity of our design platform now enables construction of protein nanomaterials according to 'back of an envelope' architectural blueprints.

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.

Blueprinting expandable nanomaterials with standardized protein building blocks.

A wooden house frame consists of many different lumber pieces, but because of the regularity of these building blocks, the structure can be designed using straightforward geometrical principles. The design of multicomponent protein assemblies in comparison has been much more complex, largely due to the irregular shapes of protein structures 1 . Here we describe extendable linear, curved, and angled protein building blocks, as well as inter-block interactions that conform to specified geometric standards; assemblies designed using these blocks inherit their extendability and regular interaction surfaces, enabling them to be expanded or contracted by varying the number of modules, and reinforced with secondary struts. Using X-ray crystallography and electron microscopy, we validate nanomaterial designs ranging from simple polygonal and circular oligomers that can be concentrically nested, up to large polyhedral nanocages and unbounded straight "train track" assemblies with reconfigurable sizes and geometries that can be readily blueprinted. Because of the complexity of protein structures and sequence-structure relationships, it has not been previously possible to build up large protein assemblies by deliberate placement of protein backbones onto a blank 3D canvas; the simplicity and geometric regularity of our design platform now enables construction of protein nanomaterials according to "back of an envelope" architectural blueprints.

Modulation of FGF pathway signaling and vascular differentiation using designed oligomeric assemblies.

Growth factors and cytokines signal by binding to the extracellular domains of their receptors and drive association and transphosphorylation of the receptor intracellular tyrosine kinase domains, initiating downstream signaling cascades. To enable systematic exploration of how receptor valency and geometry affects signaling outcomes, we designed cyclic homo-oligomers with up to 8 subunits using repeat protein building blocks that can be modularly extended. By incorporating a de novo designed fibroblast growth-factor receptor (FGFR) binding module into these scaffolds, we generated a series of synthetic signaling ligands that exhibit potent valency- and geometry-dependent Ca2+ release and MAPK pathway activation. The high specificity of the designed agonists reveal distinct roles for two FGFR splice variants in driving endothelial and mesenchymal cell fates during early vascular development. The ability to incorporate receptor binding domains and repeat extensions in a modular fashion makes our designed scaffolds broadly useful for probing and manipulating cellular signaling pathways.

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.

Membrane Protein Properties Revealed through Data-Rich Electrostatics Calculations.

The electrostatic properties of membrane proteins often reveal many of their key biophysical characteristics, such as ion channel selectivity and the stability of charged membrane-spanning segments. The Poisson-Boltzmann (PB) equation is the gold standard for calculating protein electrostatics, and the software APBSmem enables the solution of the PB equation in the presence of a membrane. Here, we describe significant advances to APBSmem, including full automation of system setup, per-residue energy decomposition, incorporation of PDB2PQR, calculation of membrane-induced pKa shifts, calculation of non-polar energies, and command-line scripting for large-scale calculations. We highlight these new features with calculations carried out on a number of membrane proteins, including the recently solved structure of the ion channel TRPV1 and a large survey of 1,614 membrane proteins of known structure. This survey provides a comprehensive list of residues with large electrostatic penalties for being embedded in the membrane, potentially revealing interesting functional information.

Interrelation between external oscillatory muscle coupling amplitude and in vivo intramedullary pressure related bone adaptation.

Interstitial bone fluid flow (IBFF) is suggested as a communication medium that bridges external physical signals and internal cellular activities in the bone, which thus regulates bone remodeling. Intramedullary pressure (ImP) is one main regulatory factor of IBFF and bone adaptation related mechanotransduction. Our group has recently observed that dynamic hydraulic stimulation (DHS), as an external oscillatory muscle coupling, was able to induce local ImP with minimal bone strain as well as to mitigate disuse bone loss. The current study aimed to evaluate the dose dependent relationship between DHS's amplitude, i.e., 15 and 30mmHg, and in vivo ImP induction, as well as this correlation on bone's phenotypic change. Simultaneous measurements of ImP and DHS cuff pressures were obtained from rats under DHS with various magnitudes and a constant frequency of 2Hz. ImP inductions and cuff pressures upon DHS loading showed a positively proportional response over the amplitude sweep. The relationship between ImP and DHS cuff pressure was evaluated and shown to be proportional, in which ImP was raised with increases of DHS cuff pressure amplitudes (R(2)=0.98). A 4-week in vivo experiment using a rat hindlimb suspension model demonstrated that the mitigation effect of DHS on disuse trabecular bone was highly dose dependent and related to DHS's amplitude, where a higher ImP led to a higher bone volume. This study suggested that sufficient physiological DHS is needed to generate ImP. Oscillatory DHS, potentially induces local fluid flow, has shown dose dependence in attenuation of disuse osteopenia.

Dynamic hydraulic fluid stimulation regulated intramedullary pressure.

Physical signals within the bone, i.e. generated from mechanical loading, have the potential to initiate skeletal adaptation. Strong evidence has pointed to bone fluid flow (BFF) as a media between an external load and the bone cells, in which altered velocity and pressure can ultimately initiate the mechanotransduction and the remodeling process within the bone. Load-induced BFF can be altered by factors such as intramedullary pressure (ImP) and/or bone matrix strain, mediating bone adaptation. Previous studies have shown that BFF induced by ImP alone, with minimum bone strain, can initiate bone remodeling. However, identifying induced ImP dynamics and bone strain factor in vivo using a non-invasive method still remains challenging. To apply ImP as a means for alteration of BFF, it was hypothesized that non-invasive dynamic hydraulic stimulation (DHS) can induce local ImP with minimal bone strain to potentially elicit osteogenic adaptive responses via bone-muscle coupling. The goal of this study was to evaluate the immediate effects on local and distant ImP and strain in response to a range of loading frequencies using DHS. Simultaneous femoral and tibial ImP and bone strain values were measured in three 15-month-old female Sprague Dawley rats during DHS loading on the tibia with frequencies of 1Hz to 10Hz. DHS showed noticeable effects on ImP induction in the stimulated tibia in a nonlinear fashion in response to DHS over the range of loading frequencies, where they peaked at 2Hz. DHS at various loading frequencies generated minimal bone strain in the tibiae. Maximal bone strain measured at all loading frequencies was less than 8µÎµ. No detectable induction of ImP or bone strain was observed in the femur. This study suggested that oscillatory DHS may regulate the local fluid dynamics with minimal mechanical strain in the bone, which serves critically in bone adaptation. These results clearly implied DHS's potential as an effective, non-invasive intervention for osteopenia and osteoporosis treatments.