Computational design of peptides that target transmembrane helices

The role of EscD in supporting EscC polymerization in the type III secretion system of enteropathogenic Escherichia coli

The type III secretion system (T3SS) is a multi-protein complex that plays a central role in the virulence of many Gram-negative bacterial pathogens. In enteropathogenic Escherichia coli, a prevalent cause of diarrheal diseases, the needle complex base of the T3SS is formed by multi-rings: two concentric inner-membrane rings made by the two oligomerizing proteins (EscD and EscJ), and an outer ring made of a single oligomerizing protein (EscC). Although the oligomerization activity of these proteins is critical for their function and can, therefore, affect the virulence of the pathogen, the mechanisms underlying the oligomerization of these proteins have yet to be identified. In this study, we report that the proteins forming the inner-membrane T3SS rings, EscJ and EscD proteins, are crucial for the oligomerization of EscC. Moreover, we elucidate the oligomerization process of EscD and determine the contribution of individual regions of the protein to its self-oligomerization activity. We show that the oligomerization motif of EscD is located at its N-terminal portion and that its transmembrane domain can self-oligomerize, thus contributing to the self-oligomerization of the full-length EscD.

Toward high-resolution computational design of the structure and function of helical membrane proteins

The computational design of α-helical membrane proteins is still in its infancy but has already made great progress. De novo design allows stable, specific and active minimal oligomeric systems to be obtained. Computational reengineering can improve the stability and function of naturally occurring membrane proteins. Currently, the major hurdle for the field is the experimental characterization of the designs. The emergence of new structural methods for membrane proteins will accelerate progress.

Inducible epitope imprinting: 'generating' the required binding site in membrane receptors for targeted drug delivery

We report an inducible epitope imprinting strategy that as a template, a flexible peptide chain can have a disordered-to-ordered conformational change by suitable inducement through a molecular imprinting procedure, and the formed nanoparticles can, in turn, remold the original peptide into the expected conformation and specifically bind to the corresponding protein.

Ionic Strength, Surface Charge, and Packing Density Effects on the Properties of Peptide Self-Assembled Monolayers

The various environmental parameters of packing density, ionic strength, and solution charge were examined for their effects on the properties of the immobilized peptide mimotope CH19 (CGSGSGSQLGPYELWELSH) that binds with the therapeutic antibody Trastuzumab (Herceptin) on a gold substrate. The immobilization of CH19 onto gold was examined with a quartz crystal microbalance (QCM). The QCM data showed the presence of intermolecular interactions resulting in the increase of viscoelastic properties of the peptide self-assembled monolayer (SAM). The CH19 SAM was diluted with CS7 (CGSGSGS) to decrease the packing density as CH19/CS7. The packing density and ionic strength parameters were evaluated by atomic force microscopy (AFM), ellipsometry, and QCM. AFM and ellipsometry showed a distinct conformational difference between CH19 and CH19/CS7, indicating a relationship between packing density and conformational state of the immobilized peptide. The CH19 SAM thickness was 40 Å with a rough topology, while the CH19/CS7 SAM thickness was 20 Å with a smooth topology. The affinity studies showed that the affinity of CH19 and CH19/CS7 to Trastuzumab were both on the order of 10⁷ M^-1 in undiluted PBS buffer, while the dilution of the buffer by 1000× increased both SAMs affinities to Trastuzumab to the order of 10¹⁵ M^-2 and changed the binding behavior from noncooperative to cooperative binding. This indicated that ionic strength had a more pronounced effect on binding properties of the CH19 SAM than packing density. Electrochemical impedance spectroscopy (EIS) was conducted on the CH19/CS7 SAM, which showed an increase in impedance after each EIS measurement cycle. Cyclic voltammetry on the CH19/CS7 SAM decreased impedance to near initial values. The impact of the packing density, buffer ionic strength, and local charge perturbation of the peptide SAM properties was interpreted based on the titratable sites in CH19 that could participate in the proton transfer and water equilibrium.

Achievements and Challenges in Computational Protein Design

Computational protein design (CPD), a yet evolving field, includes computer-aided engineering for partial or full de novo designs of proteins of interest. Designs are defined by a requested structure, function, or working environment. This chapter describes the birth and maturation of the field by presenting 101 CPD examples in a chronological order emphasizing achievements and pending challenges. Integrating these aspects presents the plethora of CPD approaches with the hope of providing a "CPD 101". These reflect on the broader structural bioinformatics and computational biophysics field and include: (1) integration of knowledge-based and energy-based methods, (2) hierarchical designated approach towards local, regional, and global motifs and the integration of high- and low-resolution design schemes that fit each such region, (3) systematic differential approaches towards different protein regions, (4) identification of key hot-spot residues and the relative effect of remote regions, (5) assessment of shape-complementarity, electrostatics and solvation effects, (6) integration of thermal plasticity and functional dynamics, (7) negative design, (8) systematic integration of experimental approaches, (9) objective cross-assessment of methods, and (10) successful ranking of potential designs. Future challenges also include dissemination of CPD software to the general use of life-sciences researchers and the emphasis of success within an in vivo milieu. CPD increases our understanding of protein structure and function and the relationships between the two along with the application of such know-how for the benefit of mankind. Applied aspects range from biological drugs, via healthier and tastier food products to nanotechnology and environmentally friendly enzymes replacing toxic chemicals utilized in the industry.

The Framework of Computational Protein Design

Computational protein design (CPD) has established itself as a leading field in basic and applied science with a strong coupling between the two. Proteins are computationally designed from the level of amino acids to the level of a functional protein complex. Design targets range from increased thermo- (or other) stability to specific requested reactions such as protein-protein binding, enzymatic reactions, or nanotechnology applications. The design scheme may encompass small regions of the proteins or the entire protein. In either case, the design may aim at the side-chains or at the full backbone conformation. Herein, the main framework for the process is outlined highlighting key elements in the CPD iterative cycle. These include the very definition of CPD, the diverse goals of CPD, components of the CPD protocol, methods for searching sequence and structure space, scoring functions, and augmenting the CPD with other optimization tools. Taken together, this chapter aims to introduce the framework of CPD.

Enhanced cellular uptake and tumor penetration of nanoparticles by imprinting the “hidden” part of membrane receptors for targeted drug delivery

Biting deep into the target cell membrane: nanoparticles targeting the transmembrane domain by trifluoroethanol-assisted epitope imprinting.

Therapeutic design of peptide modulators of protein-protein interactions in membranes

Membrane proteins play the central roles in a variety of cellular processes, ranging from nutrient uptake and signalling, to cell-cell communication. Their biological functions are directly related to how they fold and assemble; defects often lead to disease. Protein-protein interactions (PPIs) within the membrane are therefore of great interest as therapeutic targets. Here we review the progress in the application of membrane-insertable peptides for the disruption or stabilization of membrane-based PPIs. We describe the design and preparation of transmembrane peptide mimics; and of several categories of peptidomimetics used for study, including d-enantiomers, non-natural amino acids, peptoids, and β-peptides. Further aspects of the review describe modifications to membrane-insertable peptides, including lipidation and cyclization via hydrocarbon stapling. These approaches provide a pathway toward the development of metabolically stable, non-toxic, and efficacious peptide modulators of membrane-based PPIs. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.

Conformational Changes in the Epidermal Growth Factor Receptor: Role of the Transmembrane Domain Investigated by Coarse-Grained MetaDynamics Free Energy Calculations

The epidermal growth factor receptor (EGFR) is a dimeric membrane protein that regulates key aspects of cellular function. Activation of the EGFR is linked to changes in the conformation of the transmembrane (TM) domain, brought about by changes in interactions of the TM helices of the membrane lipid bilayer. Using an advanced computational approach that combines Coarse-Grained molecular dynamics and well-tempered MetaDynamics (CG-MetaD), we characterize the large-scale motions of the TM helices, simulating multiple association and dissociation events between the helices in membrane, thus leading to a free energy landscape of the dimerization process. The lowest energy state of the TM domain is a right-handed dimer structure in which the TM helices interact through the N-terminal small-X₃-small sequence motif. In addition to this state, which is thought to correspond to the active form of the receptor, we have identified further low-energy states that allow us to integrate with a high level of detail a range of previous experimental observations. These conformations may lead to the active state via two possible activation pathways, which involve pivoting and rotational motions of the helices, respectively. Molecular dynamics also reveals correlation between the conformational changes of the TM domains and of the intracellular juxtamembrane domains, paving the way for a comprehensive understanding of EGFR signaling at the cell membrane.

Regulation of integrins in platelets

Blood platelets prevent bleeding after trauma by forming occlusive aggregates at sites of vascular injury. Platelet aggregation is mediated by the integrin heterodimer αIIbβ3 and occurs when platelet agonists generated at the injury site convert αIIbβ3 from its resting to its active conformation. Active αIIbβ3 is then able to bind macromolecular ligands such as fibrinogen that crosslink adjacent platelets into hemostatic aggregates. Platelets circulate in a plasma milieu containing high concentrations of the principal αIIbβ3 ligand fibrinogen. Thus, αIIbβ3 activity is tightly regulated to prevent the spontaneous formation of platelet aggregates. αIIbβ3 activity is regulated at least three levels. First, intramolecular interactions involving motifs located in the membrane-proximal stalk regions, transmembrane domains, and the membrane-proximal cytosolic tails of αIIb and β3 maintain αIIbβ3 in its inactive conformation. Transmembrane domain interactions appear particularly important because disrupting these interactions causes constitutive αIIbβ3 activation. Second, the agonist-stimulated binding of the cytosolic proteins talin and kindlin-3 to the β3 cytosolic tail rapidly causes αIIbβ3 activation by disrupting the intramolecular interactions constraining αIIbβ3 activity. Third, the strength of ligand binding to active αIIbβ3 seems to be allosterically regulated. Thus, αIIbβ3 exists in a minimum of three interconvertible states: an inactive (resting) state that does not interact with ligands and two active ligand binding states that differ in their affinity for fibrinogen and in the mechanical stability of fibrinogen complexes they form.

Fluorophores, environments, and quantification techniques in the analysis of transmembrane helix interaction using FRET

Förster resonance energy transfer (FRET) has been widely used as a spectroscopic tool in vitro to study the interactions between transmembrane (TM) helices in detergent and lipid environments. This technique has been instrumental to many studies that have greatly contributed to quantitative understanding of the physical principles that govern helix-helix interactions in the membrane. These studies have also improved our understanding of the biological role of oligomerization in membrane proteins. In this review, we focus on the combinations of fluorophores used, the membrane mimetic environments, and measurement techniques that have been applied to study model systems as well as biological oligomeric complexes in vitro. We highlight the different formalisms used to calculate FRET efficiency and the challenges associated with accurate quantification. The goal is to provide the reader with a comparative summary of the relevant literature for planning and designing FRET experiments aimed at measuring TM helix-helix associations.

Protein design: Past, present, and future

Building on the pioneering work of Ho and DeGrado (J Am Chem Soc 1987, 109, 6751-6758) in the late 1980s, protein design approaches have revealed many fundamental features of protein structure and stability. We are now in the era that the early work presaged - the design of new proteins with practical applications and uses. Here we briefly survey some past milestones in protein design, in addition to highlighting recent progress and future aspirations.

Functional characterization of a melittin analog containing a non-natural tryptophan analog

Tryptophan (Trp) is a naturally occurring amino acid, which exhibits fluorescence emission properties that are dependent on the polarity of the local environment around the Trp side chain. However, this sensitivity also complicates interpretation of fluorescence emission data. A non-natural analogue of tryptophan, β-(1-azulenyl)-L-alanine, exhibits fluorescence insensitive to local solvent polarity and does not impact the structure or characteristics of several peptides examined. In this study, we investigated the effect of replacing Trp with β-(1-azulenyl)-L-alanine in the well-known bee-venom peptide melittin. This peptide provides a model framework for investigating the impact of replacing Trp with β-(1-azulenyl)-L-alanine in a functional peptide system that undergoes significant shifts in Trp fluorescence emission upon binding to lipid bilayers. Microbiological methods including assessment of the antimicrobial activity by minimal inhibitory concentration (MIC) assays and bacterial membrane permeability assays indicated little difference between the Trp and the β-(1-azulenyl)-L-alanine-substituted versions of melittin. Circular dichroism spectroscopy showed both that peptides adopted the expected α-helical structures when bound to phospholipid bilayers and electrophysiological analysis indicated that both created membrane disruptions leading to significant conductance increases across model membranes. Both peptides exhibited a marked protection of the respective fluorophores when bound to bilayers indicating a similar membrane-bound topology. As expected, while fluorescence quenching and CD indicate the peptides are stably bound to lipid vesicles, the peptide containing β-(1-azulenyl)-L-alanine exhibited no fluorescence emission shift upon binding while the natural Trp exhibited >10 nm shift in emission spectrum barycenter. Taken together, the β-(1-azulenyl)-L-alanine can serve as a solvent insensitive alternative to Trp that does not have significant impacts on structure or function of membrane interacting peptides.

Mechanistic and quantitative insight into cell surface targeted molecular imaging agent design

Molecular imaging agent design involves simultaneously optimizing multiple probe properties. While several desired characteristics are straightforward, including high affinity and low non-specific background signal, in practice there are quantitative trade-offs between these properties. These include plasma clearance, where fast clearance lowers background signal but can reduce target uptake, and binding, where high affinity compounds sometimes suffer from lower stability or increased non-specific interactions. Further complicating probe development, many of the optimal parameters vary depending on both target tissue and imaging agent properties, making empirical approaches or previous experience difficult to translate. Here, we focus on low molecular weight compounds targeting extracellular receptors, which have some of the highest contrast values for imaging agents. We use a mechanistic approach to provide a quantitative framework for weighing trade-offs between molecules. Our results show that specific target uptake is well-described by quantitative simulations for a variety of targeting agents, whereas non-specific background signal is more difficult to predict. Two in vitro experimental methods for estimating background signal in vivo are compared – non-specific cellular uptake and plasma protein binding. Together, these data provide a quantitative method to guide probe design and focus animal work for more cost-effective and time-efficient development of molecular imaging agents.

PepComposer: computational design of peptides binding to a given protein surface

There is a wide interest in designing peptides able to bind to a specific region of a protein with the aim of interfering with a known interaction or as starting point for the design of inhibitors. Here we describe PepComposer, a new pipeline for the computational design of peptides binding to a given protein surface. PepComposer only requires the target protein structure and an approximate definition of the binding site as input. We first retrieve a set of peptide backbone scaffolds from monomeric proteins that harbor the same backbone arrangement as the binding site of the protein of interest. Next, we design optimal sequences for the identified peptide scaffolds. The method is fully automatic and available as a web server at http://biocomputing.it/pepcomposer/webserver.

Designed multi-stranded heme binding β-sheet peptides in membrane

Designed peptides demonstrating well-defined structures and functioning in membrane environment are of significant interest in developing novel proteins for membrane active biological processes including enzymes, electron transfer, ion channels and energy conversion. Heme proteins' ability to carry out multiple functions in nature has inspired the design of several helical heme binding peptides and proteins soluble in water and also recently in membrane. Naturally occurring β-sheet proteins are both water and membrane soluble, and are known to bind heme, however, designed heme binding β-sheet proteins are yet to be reported, plausibly because of the complex folding and difficulty in introducing heme binding sites in the β-sheet structures. Here, we describe the design, NMR structures and biochemical functional characterization of four stranded and six stranded membrane soluble β-sheet peptides that bind heme and di-heme, respectively. The designed peptides contain either DP-G or DP-DA residues for the nucleation of β-turns intended to stabilize multi-stranded β-sheet topologies and ligate heme with bis-His coordination between adjacent antiparallel β-strands. Furthermore, we have optimized a high affinity heme binding pocket, Kd ∼ nM range, in the adjacent β-strands by utilizing a series of four stranded β-sheet peptides employing β- and ω-amino acids. We find that there is a progressive increase in cofactor binding affinity in the designed peptides with the alkyl chain length of ω-amino acids. Notably, the six stranded β-sheet peptide binds two molecules of heme in a cooperative fashion. The designed peptides perform peroxidase activity with varying ability and efficiently carried out electron transfer with membrane associated protein cytochrome c. The current study demonstrates the designing of functional β-sheet proteins in a membrane environment and expands the repertoire of heme protein design.

Directly Activating the Integrin αIIbβ3 Initiates Outside-In Signaling by Causing αIIbβ3 Clustering

αIIbβ3 activation in platelets is followed by activation of the tyrosine kinase c-Src associated with the carboxyl terminus of the β3 cytosolic tail. Exogenous peptides designed to interact with the αIIb transmembrane (TM) domain activate single αIIbβ3 molecules in platelets by binding to the αIIb TM domain and causing separation of the αIIbβ3 TM domain heterodimer. Here we asked whether directly activating single αIIbβ3 molecules in platelets using the designed peptide anti-αIIb TM also initiates αIIbβ3-mediated outside-in signaling by causing activation of β3-associated c-Src. Anti-αIIb TM caused activation of β3-associated c-Src and the kinase Syk, but not the kinase FAK, under conditions that precluded extracellular ligand binding to αIIbβ3. c-Src and Syk are activated by trans-autophosphorylation, suggesting that activation of individual αIIbβ3 molecules can initiate αIIbβ3 clustering in the absence of ligand binding. Consistent with this possibility, incubating platelets with anti-αIIb TM resulted in the redistribution of αIIbβ3 from a homogenous ring located at the periphery of discoid platelets into nodular densities consistent with clustered αIIbβ3. Thus, these studies indicate that not only is resting αIIbβ3 poised to undergo a conformational change that exposes its ligand-binding site, but it is poised to rapidly assemble into intracellular signal-generating complexes as well.

Neutralization of pro-inflammatory monocytes by targeting TLR2 dimerization ameliorates colitis

Monocytes have emerged as critical driving force of acute inflammation. Here, we show that inhibition of Toll-like receptor 2(TLR2) dimerization by a TLR2 transmembrane peptide (TLR2-p) ameliorated DSS-induced colitis by interfering specifically with the activation of Ly6C(+) monocytes without affecting their recruitment to the colon. We report that TLR2-p directly interacts with TLR2 within the membrane, leading to inhibition of TLR2-TLR6/1 assembly induced by natural ligands. This was associated with decreased levels of extracellular signal-regulated kinases (ERK) signaling and reduced secretion of pro-inflammatory cytokines, such as interleukin (IL)-6, IL-23, IL-12, and IL-1β. Altogether, our study provides insights into the essential role of TLR2 dimerization in the activation of pathogenic pro-inflammatory Ly6C(hi) monocytes and suggests that inhibition of this aggregation by TLR2-p might have therapeutic potential in the treatment of acute gut inflammation.

Drugging Membrane Protein Interactions

The majority of therapeutics target membrane proteins, accessible on the surface of cells, to alter cellular signaling. Cells use membrane proteins to transduce signals into cells, transport ions and molecules, bind cells to a surface or substrate, and catalyze reactions. Newly devised technologies allow us to drug conventionally "undruggable" regions of membrane proteins, enabling modulation of protein-protein, protein-lipid, and protein-nucleic acid interactions. In this review, we survey the state of the art of high-throughput screening and rational design in drug discovery, and we evaluate the advances in biological understanding and technological capacity that will drive pharmacotherapy forward against unorthodox membrane protein targets.

An overview of peptide and peptoid foldamers in medicinal chemistry

INTRODUCTION

Foldamers are artificial self-organizing systems with various critical properties: i) a stable and designable secondary structure; ii) a larger molecular surface as compared with ordinary organic drug molecules; iii) appropriate control of the orientation of the side-chain functional groups; iv) resistance against proteolytic degradation, which leads to potentially increased oral bioavailability and a longer serum half-life relative to ordinary α-peptides; and v) the lower conformational freedom may result in increased receptor binding in comparison with the natural analogs.

AREAS COVERED

This article covers the general properties and types of foldamers. This includes highlighted examples of medicinal chemical applications, including antibacterial and cargo molecules, anti-Alzheimer compounds and protein-protein interaction modifiers.

EXPERT OPINION

Various new foldamers have been created with a range of structures and biological applications. Membrane-acting antibacterial foldamers have been introduced. A general property of these structures is their amphiphilic nature. The amphiphilicity can be stationary or induced by the membrane binding. Cell-penetrating foldamers have been described which serve as cargo molecules, and foldamers have been used as autophagy inducers. Anti-Alzheimer compounds too have been created and the greatest breakthrough was attained via the modification of protein-protein interactions. This can serve as the chemical and pharmaceutical basis for the relevance of foldamers in the future.

Biologically active LIL proteins built with minimal chemical diversity

We have constructed 26-amino acid transmembrane proteins that specifically transform cells but consist of only two different amino acids. Most proteins are long polymers of amino acids with 20 or more chemically distinct side-chains. The artificial transmembrane proteins reported here are the simplest known proteins with specific biological activity, consisting solely of an initiating methionine followed by specific sequences of leucines and isoleucines, two hydrophobic amino acids that differ only by the position of a methyl group. We designate these proteins containing leucine (L) and isoleucine (I) as LIL proteins. These proteins functionally interact with the transmembrane domain of the platelet-derived growth factor β-receptor and specifically activate the receptor to transform cells. Complete mutagenesis of these proteins identified individual amino acids required for activity, and a protein consisting solely of leucines, except for a single isoleucine at a particular position, transformed cells. These surprisingly simple proteins define the minimal chemical diversity sufficient to construct proteins with specific biological activity and change our view of what can constitute an active protein in a cellular context.

Evolutionary-guided de novo structure prediction of self-associated transmembrane helical proteins with near-atomic accuracy

How specific protein associations regulate the function of membrane receptors remains poorly understood. Conformational flexibility currently hinders the structure determination of several classes of membrane receptors and associated oligomers. Here we develop EFDOCK-TM, a general method to predict self-associated transmembrane protein helical (TMH) structures from sequence guided by co-evolutionary information. We show that accurate intermolecular contacts can be identified using a combination of protein sequence covariation and TMH binding surfaces predicted from sequence. When applied to diverse TMH oligomers, including receptors characterized in multiple conformational and functional states, the method reaches unprecedented near-atomic accuracy for most targets. Blind predictions of structurally uncharacterized receptor tyrosine kinase TMH oligomers provide a plausible hypothesis on the molecular mechanisms of disease-associated point mutations and binding surfaces for the rational design of selective inhibitors. The method sets the stage for uncovering novel determinants of molecular recognition and signalling in single-spanning eukaryotic membrane receptors.

Computational design and experimental characterization of peptides intended for pH-dependent membrane insertion and pore formation

There are many opportunities to use macromolecules, such as peptides and oligonucleotides, for intracellular applications. Despite this, general methods for delivering these molecules to the cytosol in a safe and efficient manner are not available. Efforts to develop a variety of intracellular drug delivery systems such as viral vectors, lipoplexes, nanoparticles, and amphiphilic peptides have been made, but various challenges such as delivery efficiency, toxicity, and controllability remain. A central challenge is the ability to selectively perturb, not destroy, the membrane to facilitate cargo introduction. Herein, we describe our efforts to design and characterize peptides that form pores inside membranes at acidic pH, so-called pH-switchable pore formation (PSPF) peptides, as a potential means for facilitating cargo translocation through membranes. Consistent with pore formation, these peptides exhibit low-pH-triggered selective release of ATP and miRNA, but not hemoglobin, from red blood cells. Consistent with these observations, biophysical studies (tryptophan fluorescence, circular dichroism, size-exclusion chromatography, analytical ultracentrifugation, and attenuated total reflectance Fourier transformed infrared spectroscopy) show that decreased pH destabilizes the PSPF peptides in aqueous systems while promoting their membrane insertion. Together, these results suggest that reduced pH drives insertion of PSPF peptides into membranes, leading to target-specific escape through a proposed pore formation mechanism.

Computational modeling of membrane proteins

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1-2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug-specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans-membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α-helical MPs as well as β-barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge-based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade.

Computational design of selective peptides to discriminate between similar PDZ domains in an oncogenic pathway

Reagents that target protein-protein interactions to rewire signaling are of great relevance in biological research. Computational protein design may offer a means of creating such reagents on demand, but methods for encoding targeting selectivity are sorely needed. This is especially challenging when targeting interactions with ubiquitous recognition modules--for example, PDZ domains, which bind C-terminal sequences of partner proteins. Here we consider the problem of designing selective PDZ inhibitor peptides in the context of an oncogenic signaling pathway, in which two PDZ domains (NHERF-2 PDZ2-N2P2 and MAGI-3 PDZ6-M3P6) compete for a receptor C-terminus to differentially modulate oncogenic activities. Because N2P2 has been shown to increase tumorigenicity and M3P6 to decreases it, we sought to design peptides that inhibit N2P2 without affecting M3P6. We developed a structure-based computational design framework that models peptide flexibility in binding yet is efficient enough to rapidly analyze tradeoffs between affinity and selectivity. Designed peptides showed low-micromolar inhibition constants for N2P2 and no detectable M3P6 binding. Peptides designed for reverse discrimination bound M3P6 tighter than N2P2, further testing our technology. Experimental and computational analysis of selectivity determinants revealed significant indirect energetic coupling in the binding site. Successful discrimination between N2P2 and M3P6, despite their overlapping binding preferences, is highly encouraging for computational approaches to selective PDZ targeting, especially because design relied on a homology model of M3P6. Still, we demonstrate specific deficiencies of structural modeling that must be addressed to enable truly robust design. The presented framework is general and can be applied in many scenarios to engineer selective targeting.

Dimerization of the EphA1 receptor tyrosine kinase transmembrane domain: Insights into the mechanism of receptor activation

EphA1 is a receptor tyrosine kinase (RTK) that plays a key role in developmental processes, including guidance of the migration of axons and cells in the nervous system. EphA1, in common with other RTKs, contains an N-terminal extracellular domain, a single transmembrane (TM) α-helix, and a C-terminal intracellular kinase domain. The TM helix forms a dimer, as seen in recent NMR studies. We have modeled the EphA1 TM dimer using a multiscale approach combining coarse-grain (CG) and atomistic molecular dynamics (MD) simulations. The one-dimensional potential of mean force (PMF) for this system, based on interhelix separation, has been calculated using CG MD simulations. This provides a view of the free energy landscape for helix–helix interactions of the TM dimer in a lipid bilayer. The resulting PMF profiles suggest two states, consistent with a rotation-coupled activation mechanism. The more stable state corresponds to a right-handed helix dimer interacting via an N-terminal glycine zipper motif, consistent with a recent NMR structure (2K1K). A second metastable state corresponds to a structure in which the glycine zipper motif is not involved. Analysis of unrestrained CG MD simulations based on representative models from the PMF calculations or on the NMR structure reveals possible pathways of interconversion between these two states, involving helix rotations about their long axes. This suggests that the interaction of TM helices in EphA1 dimers may be intrinsically dynamic. This provides a potential mechanism for signaling whereby extracellular events drive a shift in the repopulation of the underlying TM helix dimer energy landscape.

Talin-driven inside-out activation mechanism of platelet αIIbβ3 integrin probed by multimicrosecond, all-atom molecular dynamics simulations

Platelet aggregation is the consequence of the binding of extracellular bivalent ligands such as fibrinogen and von Willebrand factor to the high affinity, active state of integrin αIIbβ3. This state is achieved through a so-called "inside-out" mechanism characterized by the membrane-assisted formation of a complex between the F2 and F3 subdomains of intracellular protein talin and the integrin β3 tail. Here, we present the results of multi-microsecond, all-atom molecular dynamics simulations carried on the complete transmembrane (TM) and C-terminal (CT) domains of αIIbβ3 integrin in an explicit lipid-water environment, and in the presence or absence of the talin-1 F2 and F3 subdomains. These large-scale simulations provide unprecedented molecular-level insights into the talin-driven inside-out activation of αIIbβ3 integrin. Specifically, they suggest a preferred conformation of the complete αIIbβ3 TM/CT domains in a lipid-water environment, and testable hypotheses of key intermolecular interactions between αIIbβ3 integrin and the F2/F3 domains of talin-1. Notably, not only do these simulations give support to a stable left-handed reverse turn conformation of the αIIb juxtamembrane motif rather than a helical turn, but they raise the question as to whether TM helix separation is required for talin-driven integrin activation.

Photoactive Spatial Proximity Probes for Binding Pairs with Epigenetic Marks

A new strategy for encoding polypeptide libraries with photolabile tags is developed. The photoassisted assay, based on conditional release of encoding tags only from bound pairs, can differentiate between peptides which have minor differences in a form of post-translational modifications with epigenetic marks. The encoding strategy is fully compatible with automated peptide synthesis. The encoding pendants are compact and do not perturb potential binding interactions.

The HIV-1 envelope transmembrane domain binds TLR2 through a distinct dimerization motif and inhibits TLR2-mediated responses

HIV-1 uses a number of means to manipulate the immune system, to avoid recognition and to highjack signaling pathways. HIV-1 infected cells show limited Toll-Like Receptor (TLR) responsiveness via as yet unknown mechanisms. Using biochemical and biophysical approaches, we demonstrate that the trans-membrane domain (TMD) of the HIV-1 envelope (ENV) directly interacts with TLR2 TMD within the membrane milieu. This interaction attenuates TNFα, IL-6 and MCP-1 secretion in macrophages, induced by natural ligands of TLR2 both in in vitro and in vivo models. This was associated with decreased levels of ERK phosphorylation. Furthermore, mutagenesis demonstrated the importance of a conserved GxxxG motif in driving this interaction within the membrane milieu. The administration of the ENV TMD in vivo to lipotechoic acid (LTA)/Galactosamine-mediated septic mice resulted in a significant decrease in mortality and in tissue damage, due to the weakening of systemic macrophage activation. Our findings suggest that the TMD of ENV is involved in modulation of the innate immune response during HIV infection. Furthermore, due to the high functional homology of viral ENV proteins this function may be a general character of viral-induced immune modulation.

To understand viral pathology and the tools needed to eliminate infection, it is important to understand how viral immune evasion occurs. One such mode of inhibition is the decreased responsiveness of Toll-Like Receptors (TLRs). To date, the exact mechanism inducing this inhibition is not clear. In this study, we utilized a multidisciplinary approach and report on direct modulation of TLR2 activity by the envelope trans-membrane protein of HIV-1 through trans-membrane domain interactions. This interaction resulted in a decreased response in vitro of TLR2 to its natural ligand LTA. Through mutagenesis analysis we show that the GxxxG motif is the driving force of this interaction. Interestingly, the inhibitory effect was also highly effective in protecting mice from lethal effects in a sepsis-like model. Our findings implicate that ENV participates in innate immune impairment, which may occur during viral entry and at latent stages. Furthermore, due to the high functional homology between viral ENV proteins, this function may exhibit a general character of viral-induced immune modulation.

Transmembrane recognition of the semaphorin co-receptors neuropilin 1 and plexin A1: coarse-grained simulations

The cancer associated class 3 semaphorins require direct binding to neuropilins and association to plexins to trigger cell signaling. Here, we address the role of the transmembrane domains of neuropilin 1 and plexin A1 for the dimerization of the two receptors by characterizing the assembly in lipid bilayers using coarse-grained molecular dynamics simulations. From experimental evidence using a two-hybrid system showing the biochemical association of the two receptors transmembrane domains, we performed molecular simulations in DOPC and POPC demonstrating spontaneously assembly to form homodimers and heterodimers with a very high propensity for right-handed packing of the helices. Inversely, left-handed packing was observed with a very low propensity. This mode of packing was observed uniquely when the plexin A1 transmembrane domain was involved in association. Potential of mean force calculations were used to predict a hierarchy of self-association for the monomers: the two neuropilin 1 transmembrane domains strongly associated, neuropilin 1 and plexin A1 transmembrane domains associated less and the two plexin A1 transmembrane domains weakly but significantly associated. We demonstrated that homodimerization and heterodimerization are driven by GxxxG motifs, and that the sequence context modulates the packing mode of the plexin A1 transmembrane domains. This work presents major advances towards our understanding of membrane signaling platforms assembly through membrane domains and provides exquisite information for the design of antagonist drugs defining a novel class of therapeutic agents.

Designed di-heme binding helical transmembrane protein

De novo designing of functional membrane proteins is fundamental in terms of understanding the structure, folding, and stability of membrane proteins. In this work, we report the design and characterization of a transmembrane protein, termed HETPRO (HEme-binding Transmembrane PROtein), that binds two molecules of heme in a membrane and catalyzes oxidation/reduction reactions. The primary structure of HETPRO has been optimized in a guided fashion, from an antimicrobial peptide, for transmembrane orientation, defined 3D structure, and functions. HETPRO assembles into a tetrameric form, from an apo dimeric helical structure, in complex with cofactor in detergent micelles. The NMR structure of the apo HETPRO in micelles reveals an antiparallel helical dimer that inserts into the nonpolar core of detergent micelles. The well-defined structure of HETPRO and its ability to bind to heme moieties could be utilized to develop a functional membrane protein mimic for electron transport and photosystems.

Predicting the Impact of Missense Mutations on Protein-Protein Binding Affinity

The crucial prerequisite for proper biological function is the protein's ability to establish highly selective interactions with macromolecular partners. A missense mutation that alters the protein binding affinity may cause significant perturbations or complete abolishment of the function, potentially leading to diseases. The availability of computational methods to evaluate the impact of mutations on protein-protein binding is critical for a wide range of biomedical applications. Here, we report an efficient computational approach for predicting the effect of single and multiple missense mutations on protein-protein binding affinity. It is based on a well-tested simulation protocol for structure minimization, modified MM-PBSA and statistical scoring energy functions with parameters optimized on experimental sets of several thousands of mutations. Our simulation protocol yields very good agreement between predicted and experimental values with Pearson correlation coefficients of 0.69 and 0.63 and root-mean-square errors of 1.20 and 1.90 kcal mol^-1 for single and multiple mutations, respectively. Compared with other available methods, our approach achieves high speed and prediction accuracy and can be applied to large datasets generated by modern genomics initiatives. In addition, we report a crucial role of water model and the polar solvation energy in estimating the changes in binding affinity. Our analysis also reveals that prediction accuracy and effect of mutations on binding strongly depends on the type of mutation and its location in a protein complex.

De novo selection of oncogenes

All cellular proteins are derived from preexisting ones by natural selection. Because of the random nature of this process, many potentially useful protein structures never arose or were discarded during evolution. Here, we used a single round of genetic selection in mouse cells to isolate chemically simple, biologically active transmembrane proteins that do not contain any amino acid sequences from preexisting proteins. We screened a retroviral library expressing hundreds of thousands of proteins consisting of hydrophobic amino acids in random order to isolate four 29-aa proteins that induced focus formation in mouse and human fibroblasts and tumors in mice. These proteins share no amino acid sequences with known cellular or viral proteins, and the simplest of them contains only seven different amino acids. They transformed cells by forming a stable complex with the platelet-derived growth factor β receptor transmembrane domain and causing ligand-independent receptor activation. We term this approach de novo selection and suggest that it can be used to generate structures and activities not observed in nature, create prototypes for novel research reagents and therapeutics, and provide insight into cell biology, transmembrane protein-protein interactions, and possibly virus evolution and the origin of life.

Hepatitis C virus RNA replication and virus particle assembly require specific dimerization of the NS4A protein transmembrane domain

Hepatitis C virus (HCV) NS4A is a single-pass transmembrane (TM) protein essential for viral replication and particle assembly. The sequence of the NS4A TM domain is highly conserved, suggesting that it may be important for protein-protein interactions. To test this hypothesis, we measured the potential dimerization of the NS4A TM domain in a well-characterized two-hybrid TM protein interaction system. The NS4A TM domain exhibited a strong homotypic interaction that was comparable in affinity to glycophorin A, a well-studied human blood group antigen that forms TM homodimers. Several mutations predicted to cluster on a common surface of the NS4A TM helix caused significant reductions in dimerization, suggesting that these residues form an interface for NS4A dimerization. Mutations in the NS4A TM domain were further examined in the JFH-1 genotype 2a replicon system; importantly, all mutations that destabilized NS4A dimers also caused defects in RNA replication and/or virus assembly. Computational modeling of NS4A TM interactions suggests a right-handed dimeric interaction of helices with an interface that is consistent with the mutational effects. Furthermore, defects in NS4A oligomerization and virus particle assembly of two mutants were rescued by NS4A A15S, a TM mutation recently identified through forward genetics as a cell culture-adaptive mutation. Together, these data provide the first example of a functionally important TM dimer interface within an HCV nonstructural protein and reveal a fundamental role of the NS4A TM domain in coordinating HCV RNA replication and virus particle assembly.