A Pareto-Optimal Refinement Method for Protein Design Scaffolds

Computational design of protein function involves a search for amino acids with the lowest energy subject to a set of constraints specifying function. In many cases a set of natural protein backbone structures, or “scaffolds”, are searched to find regions where functional sites (an enzyme active site, ligand binding pocket, protein – protein interaction region, etc.) can be placed, and the identities of the surrounding amino acids are optimized to satisfy functional constraints. Input native protein structures almost invariably have regions that score very poorly with the design force field, and any design based on these unmodified structures may result in mutations away from the native sequence solely as a result of the energetic strain. Because the input structure is already a stable protein, it is desirable to keep the total number of mutations to a minimum and to avoid mutations resulting from poorly-scoring input structures. Here we describe a protocol using cycles of minimization with combined backbone/sidechain restraints that is Pareto-optimal with respect to RMSD to the native structure and energetic strain reduction. The protocol should be broadly useful in the preparation of scaffold libraries for functional site design.     

A Membrane Topology Model for Human Interferon Inducible Transmembrane Protein 1

InterFeron Inducible TransMembrane proteins 1–3 (IFITM1, IFITM2 and IFITM3) are a family of proteins capable of inhibiting the cellular entry of numerous human and animal viruses. IFITM1-3 are unique amongst the currently described viral restriction factors in their apparent ability to block viral entry. This restrictive property is dependant on the localisation of the proteins to plasma and endosomal membranes, which constitute the main portals of viral entry into cells. The topology of the IFITM proteins within cell membranes is an unresolved aspect of their biology. Here we present data from immunofluorescence microscopy, protease cleavage, biotin-labelling and immuno-electron microscopy assays, showing that human IFITM1 has a membrane topology in which the N-terminal domain resides in the cytoplasm, and the C-terminal domain is extracellular. Furthermore, we provide evidence that this topology is conserved for all of the human interferon-induced IFITM proteins. This model is consistent with that recently proposed for murine IFITM3, but differs from that proposed for murine IFITM1.     

A Phase-Field Model for Articular Cartilage Regeneration in Degradable Scaffolds

Degradable scaffolds represent a promising solution for tissue engineering of damaged or degenerated articular cartilage which due to its avascular nature, is characterized by a low self-repair capacity. To estimate the articular cartilage regeneration process employing degradable scaffolds, we propose a mathematical model as the extension of Olson and Haider’s work (Int. J. Pure Appl. Math. 53:333–353, 2009). The simulated tissue engineering procedure consists in (i) the explant of a cylindrical sample, (ii) the removal of the inner core region, and (iii) the filling of the inner region with hydrogels, degradable scaffolds enriched with nutrients, such as oxygen and glucose. The phase-field model simulates the cartilage regeneration process at the scaffold-cartilage interface. It embeds reaction-diffusion equations, which are used to model the nutrient and regenerated extracellular matrix. The equations are solved using an unconditionally stable hybrid numerical scheme. Cartilage repair processes with full-thickness defects, which are controlled by properties of hydrogel materials and cartilage explant culture based on biological interest are observed. The implemented mathematical model shows the capability to simulate cartilage repairing processes, which can be virtually controlled evaluating hydrogel and cartilage material properties including nutrient supply and defected magnitude. In particular, the adopted methodology is able to explain the regeneration time of cartilage within hydrogel environments. With the numerical scheme, the numerical simulations are demonstrated for the potential improvement of hydrogel structures.     

A Simple Detection Method for Low-Affinity Membrane Protein Interactions by Baculoviral Display

Background

Membrane protein interactions play an important role in cell-to-cell recognition in various biological activities such as in the immune or neural system. Nevertheless, there has remained the major obstacle of expression of the membrane proteins in their active form. Recently, we and other investigators found that functional membrane proteins express on baculovirus particles (budded virus, BV). In this study, we applied this BV display system to detect interaction between membrane proteins important for cell-to-cell interaction in immune system.

Methodology/Principal Findings

We infected Sf9 cells with recombinant baculovirus encoding the T cell membrane protein CD2 or its ligand CD58 and recovered the BV. We detected specific interaction between CD2-displaying BV and CD58-displaying BV by an enzyme-linked immunosorbent assay (ELISA). Using this system, we also detected specific interaction between two other membrane receptor-ligand pairs, CD40-CD40 ligand (CD40L), and glucocorticoid-induced TNFR family-related protein (GITR)-GITR ligand (GITRL). Furthermore, we observed specific binding of BV displaying CD58, CD40L, or GITRL to cells naturally expressing their respective receptors by flowcytometric analysis using anti-baculoviral gp64 antibody. Finally we isolated CD2 cDNA from a cDNA expression library by magnetic separation using CD58-displayng BV and anti-gp64 antibody.

Conclusions

We found the BV display system worked effectively in the detection of the interaction of membrane proteins. Since various membrane proteins and their oligomeric complexes can be displayed on BV in the native form, this BV display system should prove highly useful in the search for natural ligands or to develop screening systems for therapeutic antibodies and/or compounds.     

Membrane Morphologies Induced by Arc-Shaped Scaffolds Are Determined by Arc Angle and Coverage

The intricate shapes of biological membranes such as tubules and membrane stacks are induced by proteins. In this article, we systematically investigate the membrane shapes induced by arc-shaped scaffolds such as proteins and protein complexes with coarse-grained modeling and simulations. We find that arc-shaped scaffolds induce membrane tubules at membrane coverages larger than a threshold of ～40%, irrespective of their arc angle. The membrane morphologies at intermediate coverages below this tubulation threshold, in contrast, strongly depend on the arc angle. Scaffolds with arc angles of about 60°, akin to N-BAR domains, do not change the membrane shape at coverages below the tubulation threshold, whereas scaffolds with arc angles larger than about 120° induce double-membrane stacks at intermediate coverages. The scaffolds stabilize the curved membrane edges that connect the membrane stacks, as suggested for complexes of reticulon proteins. Our results provide general insights on the determinants of membrane shaping by arc-shaped scaffolds.     

Sensing Membrane Stresses by Protein Insertions

Protein domains shallowly inserting into the membrane matrix are ubiquitous in peripheral membrane proteins involved in various processes of intracellular membrane shaping and remodeling. It has been suggested that these domains sense membrane curvature through their preferable binding to strongly curved membranes, the binding mechanism being mediated by lipid packing defects. Here we make an alternative statement that shallow protein insertions are universal sensors of the intra-membrane stresses existing in the region of the insertion embedding rather than sensors of the curvature per se. We substantiate this proposal computationally by considering different independent ways of the membrane stress generation among which some include changes of the membrane curvature whereas others do not alter the membrane shape. Our computations show that the membrane-binding coefficient of shallow protein insertions is determined by the resultant stress independently of the way this stress has been produced. By contrast, consideration of the correlation between the insertion binding and the membrane curvature demonstrates that the binding coefficient either increases or decreases with curvature depending on the factors leading to the curvature generation. To validate our computational model, we treat quantitatively the experimental results on membrane binding by ALPS1 and ALPS2 motifs of ArfGAP1.     

Lipid Segregation and Membrane Budding Induced by the Peripheral Membrane Binding Protein Annexin A2

The formation of dynamic membrane microdomains is an important phenomenon in many signal transduction and membrane trafficking events. It is driven by intrinsic properties of membrane lipids and integral as well as membrane-associated proteins. Here we analyzed the ability of one peripherally associated membrane protein, annexin A2 (AnxA2), to induce the formation of phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂)-rich domains in giant unilamellar vesicles (GUVs) of complex lipid composition. AnxA2 is a cytosolic protein that can bind PI(4,5)P₂ and other acidic phospholipids in a Ca²⁺-dependent manner and that has been implicated in cellular membrane dynamics in endocytosis and exocytosis. We show that AnxA2 binding to GUVs induces lipid phase separation and the recruitment of PI(4,5)P₂, cholesterol and glycosphingolipids into larger clusters. This property is observed for the full-length monomeric protein, a mutant derivative comprising the C-terminal protein core domain and for AnxA2 residing in a heterotetrameric complex with its intracellular binding partner S100A10. All AnxA2 derivatives inducing PI(4,5)P₂ clustering are also capable of forming interconnections between PI(4,5)P₂-rich microdomains of adjacent GUVs. Furthermore, they can induce membrane indentations rich in PI(4,5)P₂ and inward budding of these membrane domains into the lumen of GUVs. This inward vesiculation is specific for AnxA2 and not shared with other PI(4,5)P₂-binding proteins such as the pleckstrin homology (PH) domain of phospholipase Cδ1. Together our results indicate that annexins such as AnxA2 can efficiently induce membrane deformations after lipid segregation, a mechanism possibly underlying annexin functions in membrane trafficking.     

Mechanism for Active Membrane Fusion Triggering by Morbillivirus Attachment Protein

The paramyxovirus entry machinery consists of two glycoproteins that tightly cooperate to achieve membrane fusion for cell entry: the tetrameric attachment protein (HN, H, or G, depending on the paramyxovirus genus) and the trimeric fusion protein (F). Here, we explore whether receptor-induced conformational changes within morbillivirus H proteins promote membrane fusion by a mechanism requiring the active destabilization of prefusion F or by the dissociation of prefusion F from intracellularly preformed glycoprotein complexes. To properly probe F conformations, we identified anti-F monoclonal antibodies (MAbs) that recognize conformation-dependent epitopes. Through heat treatment as a surrogate for H-mediated F triggering, we demonstrate with these MAbs that the morbillivirus F trimer contains a sufficiently high inherent activation energy barrier to maintain the metastable prefusion state even in the absence of H. This notion was further validated by exploring the conformational states of destabilized F mutants and stabilized soluble F variants combined with the use of a membrane fusion inhibitor (3g). Taken together, our findings reveal that the morbillivirus H protein must lower the activation energy barrier of metastable prefusion F for fusion triggering.     

Membrane-Mediated Interaction between Strongly Anisotropic Protein Scaffolds

Specialized proteins serve as scaffolds sculpting strongly curved membranes of intracellular organelles. Effective membrane shaping requires segregation of these proteins into domains and is, therefore, critically dependent on the protein-protein interaction. Interactions mediated by membrane elastic deformations have been extensively analyzed within approximations of large inter-protein distances, small extents of the protein-mediated membrane bending and small deviations of the protein shapes from isotropic spherical segments. At the same time, important classes of the realistic membrane-shaping proteins have strongly elongated shapes with large and highly anisotropic curvature. Here we investigated, computationally, the membrane mediated interaction between proteins or protein oligomers representing membrane scaffolds with strongly anisotropic curvature, and addressed, quantitatively, a specific case of the scaffold geometrical parameters characterizing BAR domains, which are crucial for membrane shaping in endocytosis. In addition to the previously analyzed contributions to the interaction, we considered a repulsive force stemming from the entropy of the scaffold orientation. We computed this interaction to be of the same order of magnitude as the well-known attractive force related to the entropy of membrane undulations. We demonstrated the scaffold shape anisotropy to cause a mutual aligning of the scaffolds and to generate a strong attractive interaction bringing the scaffolds close to each other to equilibrium distances much smaller than the scaffold size. We computed the energy of interaction between scaffolds of a realistic geometry to constitute tens of k_BT, which guarantees a robust segregation of the scaffolds into domains.     

A Self-Scaffolding Model for G Protein Signaling

Activation of heterotrimeric G proteins is generally believed to induce dissociation of Gα and Gβγ subunits, which are then free to bind to and change the catalytic activity of a variety of intracellular enzymes. We have previously found that in cells, Gαq subunits remain complexed with its major effector, phospholipase Cβ1, through the activation cycle. To determine whether this behavior may be operative in other systems, we carried out Förster resonance energy transfer studies and found that eYFP-Gαi and eCFP-Gβγ remain associated after stimulation in HEK293 cells. We also found that the level of Forster resonance energy transfer between Alexa546-phospholipase Cβ2 and eGFP-Gβγ is significant and unchanged upon activation in HEK293 cells, thus showing that these proteins can localize into stable signaling complexes. To understand the basis for this stabilization, we carried out in vitro studies using a series of single-Cys mutants labeled with fluorescence tags and monitored their interaction with Gβγ subunits and changes in their fluorescence properties and accessibility upon activation and Gβγ binding. Our studies suggest a significant change in the orientation between G protein subunits upon activation that allows the G proteins to remain complexed while activating effectors.     

A Multi-label Classifier for Prediction Membrane Protein Functional Types in Animal

Membrane protein is an important composition of cell membrane. Given a membrane protein sequence, how can we identify its type(s) is very important because the type keeps a close correlation with its functions. According to previous studies, membrane protein can be divided into the following eight types: single-pass type I, single-pass type II, single-pass type III, single-pass type IV, multipass, lipid-anchor, GPI-anchor, peripheral membrane protein. With the avalanche of newly found protein sequences in the post-genomic age, it is urgent to develop an automatic and effective computational method to rapid and reliable prediction of the types of membrane proteins. At present, most of the existing methods were based on the assumption that one membrane protein only belongs to one type. Actually, a membrane protein may simultaneously exist at two or more different functional types. In this study, a new method by hybridizing the pseudo amino acid composition with multi-label algorithm called LIFT (multi-label learning with label-specific features) was proposed to predict the functional types both singleplex and multiplex animal membrane proteins. Experimental result on a stringent benchmark dataset of membrane proteins by jackknife test show that the absolute-true obtained was 0.6342, indicating that our approach is quite promising. It may become a useful high-through tool, or at least play a complementary role to the existing predictors in identifying functional types of membrane proteins.     

A Role for Weak Electrostatic Interactions in Peripheral Membrane Protein Binding

Bacillus thuringiensis phosphatidylinositol-specific phospholipase C (BtPI-PLC) is a secreted virulence factor that binds specifically to phosphatidylcholine (PC) bilayers containing negatively charged phospholipids. BtPI-PLC carries a negative net charge and its interfacial binding site has no obvious cluster of basic residues. Continuum electrostatic calculations show that, as expected, nonspecific electrostatic interactions between BtPI-PLC and membranes vary as a function of the fraction of anionic lipids present in the bilayers. Yet they are strikingly weak, with a calculated ΔG_el below 1 kcal/mol, largely due to a single lysine (K44). When K44 is mutated to alanine, the equilibrium dissociation constant for small unilamellar vesicles increases more than 50 times (∼2.4 kcal/mol), suggesting that interactions between K44 and lipids are not merely electrostatic. Comparisons of molecular-dynamics simulations performed using different lipid compositions reveal that the bilayer composition does not affect either hydrogen bonds or hydrophobic contacts between the protein interfacial binding site and bilayers. However, the occupancies of cation-π interactions between PC choline headgroups and protein tyrosines vary as a function of PC content. The overall contribution of basic residues to binding affinity is also context dependent and cannot be approximated by a rule-of-thumb value because these residues can contribute to both nonspecific electrostatic and short-range protein-lipid interactions. Additionally, statistics on the distribution of basic amino acids in a data set of membrane-binding domains reveal that weak electrostatics, as observed for BtPI-PLC, might be a less unusual mechanism for peripheral membrane binding than is generally thought.     

分析超速离心技术研究膜蛋白TmrAB

去垢剂在膜蛋白的提取纯化过程中起到必要的作用,对膜蛋白的聚合状态、结晶条件以及理化性质等方面都有较大影响.分析超速离心技术(analytical ultracentrifuge,AUC)通过测定溶液中膜蛋白-去垢剂复合物在离心场中的沉降运动轨迹,可以分析获得其沉降系数、摩尔质量、流体力学半径、结合常数等水力学和热力学性质,进而判断膜蛋白-去垢剂复合物的均一性及聚合状态.本文以嗜热菌来源的ATP结合转运蛋白(ABC transporter)TmrAB作为研究对象,利用分析超速离心技术结合分子排阻层析和冷冻电镜负染技术,研究其均一性、聚合状态以及去垢剂与膜蛋白的摩尔比.结果显示,在8倍临界胶束浓度(critical micelle concentration,CMC)的DDM条件下,TmrAB性质均一,并以异二聚体的单体形式存在,DDM与Tmr AB的摩尔比为116∶1.本研究表明,分析超速离心技术是一种测定膜蛋白分子质量、研究膜蛋白聚合状态的可靠手段.     

Assembling a Correctly Folded and Functional Heptahelical Membrane Protein by Protein Trans-splicing

Protein trans-splicing using split inteins is well established as a useful tool for protein engineering. Here we show, for the first time, that this method can be applied to a membrane protein under native conditions. We provide compelling evidence that the heptahelical proteorhodopsin can be assembled from two separate fragments consisting of helical bundles A and B and C, D, E, F, and G via a splicing site located in the BC loop. The procedure presented here is on the basis of dual expression and ligation in vivo. Global fold, stability, and photodynamics were analyzed in detergent by CD, stationary, as well as time-resolved optical spectroscopy. The fold within lipid bilayers has been probed by high field and dynamic nuclear polarization-enhanced solid-state NMR utilizing a ¹³C-labeled retinal cofactor and extensively ¹³C-¹⁵N-labeled protein. Our data show unambiguously that the ligation product is identical to its non-ligated counterpart. Furthermore, our data highlight the effects of BC loop modifications onto the photocycle kinetics of proteorhodopsin. Our data demonstrate that a correctly folded and functionally intact protein can be produced in this artificial way. Our findings are of high relevance for a general understanding of the assembly of membrane proteins for elucidating intramolecular interactions, and they offer the possibility of developing novel labeling schemes for spectroscopic applications.     

Vaccinia Viral Protein A27 Is Anchored to the Viral Membrane via a Cooperative Interaction with Viral Membrane Protein A17

The vaccinia viral protein A27 in mature viruses specifically interacts with heparan sulfate for cell surface attachment. In addition, A27 associates with the viral membrane protein A17 to anchor to the viral membrane; however, the specific interaction between A27 and A17 remains largely unclear. To uncover the active binding sites and the underlying binding mechanism, we expressed and purified the N-terminal (18–50 residues) and C-terminal (162–203 residues) fragments of A17, which are denoted A17-N and A17-C. Through surface plasmon resonance, the binding affinity of A27/A17-N (K_A = 3.40 × 10⁸ m⁻¹) was determined to be approximately 3 orders of magnitude stronger than that of A27/A17-C (K_A = 3.40 × 10⁵ m⁻¹), indicating that A27 prefers to interact with A17-N rather than A17-C. Despite the disordered nature of A17-N, the A27-A17 interaction is mediated by a specific and cooperative binding mechanism that includes two active binding sites, namely ³²SFMPK³⁶ (denoted as F1 binding) and ²⁰LDKDLFTEEQ²⁹ (F2). Further analysis showed that F1 has stronger binding affinity and is more resistant to acidic conditions than is F2. Furthermore, A27 mutant proteins that retained partial activity to interact with the F1 and F2 sites of the A17 protein were packaged into mature virus particles at a reduced level, demonstrating that the F1/F2 interaction plays a critical role in vivo. Using these results in combination with site-directed mutagenesis data, we established a computer model to explain the specific A27-A17 binding mechanism.