De novo peptide design with C3a receptor agonist and antagonist activities: theoretical predictions and experimental validation

In response to complement anaphylatoxin peptides C3a and C5a, human vascular endothelial cells migrate and mediate the activation of B-cells and polarization of T-cells

The vascular endothelium has been discovered in the past several years to be important in shaping the cellular immune response. During the immune response the vascular endothelium is constantly perturbed by biologically potent molecules, including the complement activation peptides, C3a and C5a. Despite the importance of C3a and C5a in inflammation and immunity, their role in modulating lymphocyte function via activation of vascular endothelial cells is unknown. Accordingly, we investigated the regulated expression of the C3a and C5a receptors (complement anaphylatoxin C3a receptor [C3aR] and complement anaphylatoxin C5a receptor 1 [C5aR1]) on human umbilical vascular endothelial cells (HUVECs) and examined how C3a or C5a activation of HUVECs affects the activation and polarization of lymphatic cells. Our findings demonstrated that C3a and C5a increase C3aR and C5aR1 expression by HUVECs as well as directing their cellular transmigration and spreading through transwell filters. Moreover, C3a- or C5a-stimulated endothelial cells: (1) caused activation of B-lymphoblasts with significant increase in Fas Ligand (CD95L) (FasL), CD69, and IL-R1 expression, and (2) skewed T-lymphoblast cells toward a Th1 subtype, (CD4⁺ /CCR5⁺ ) that correlated with significant increase of IFN-γ. Collectively, these data indicate that C3a and C5a signaling is important in the activation and polarization of lymphocytes as they traffic through the vascular endothelium during the immune response.

Rigorous Computational and Experimental Investigations on MDM2/MDMX-Targeted Linear and Macrocyclic Peptides

There is interest in peptide drug design, especially for targeting intracellular protein–protein interactions. Therefore, the experimental validation of a computational platform for enabling peptide drug design is of interest. Here, we describe our peptide drug design platform (CMDInventus) and demonstrate its use in modeling and predicting the structural and binding aspects of diverse peptides that interact with oncology targets MDM2/MDMX in comparison to both retrospective (pre-prediction) and prospective (post-prediction) data. In the retrospective study, CMDInventus modules (CMDpeptide, CMDboltzmann, CMDescore and CMDyscore) were used to accurately reproduce structural and binding data across multiple MDM2/MDMX data sets. In the prospective study, CMDescore, CMDyscore and CMDboltzmann were used to accurately predict binding affinities for an Ala-scan of the stapled α-helical peptide ATSP-7041. Remarkably, CMDboltzmann was used to accurately predict the results of a novel D-amino acid scan of ATSP-7041. Our investigations rigorously validate CMDInventus and support its utility for enabling peptide drug design.

The Complement Anaphylatoxins C5a and C3a Suppress IFN-β Production in Response to Listeria monocytogenes by Inhibition of the Cyclic Dinucleotide-Activated Cytosolic Surveillance Pathway

Listeria monocytogenes is an intracellular Gram-positive bacterium that induces expression of type I IFNs (IFN-α/IFN-β) during infection. These cytokines are detrimental to the host during infection by priming leukocytes to undergo L. monocytogenes-mediated apoptosis. Our previous studies showed that C5aR1^-/- and C3aR^-/- mice are highly susceptible to L. monocytogenes infection as a result of increased IFN-β-mediated apoptosis of major leukocyte cell populations, including CD4⁺ and CD8⁺ T cells. However, the mechanisms by which C3a and C5a modulate IFN-β expression during L. monocytogenes infection were not examined in these initial investigations. Accordingly, we report in this article that C5a and C3a suppress IFN-β production in response to L. monocytogenes via cyclic di-AMP (c-di-AMP), a secondary messenger molecule of L. monocytogenes, in J774A.1 macrophage-like cells and in bone marrow-derived dendritic cells (BMDCs). Moreover, C5a and C3a suppress IFN-β production by acting through their respective receptors, because no inhibition was seen in C5aR1^-/- or C3aR^-/- BMDCs, respectively. C5a and C3a suppress IFN-β production in a manner that is dependent on Bruton's tyrosine kinase, p38 MAPK, and TANK-binding kinase 1 (TBK1), as demonstrated by the individual use of Bruton's tyrosine kinase, p38 MAPK, and TBK1 inhibitors. Pretreatment of cells with C5a and C3a reduced the expression of the IFN-β signaling molecules DDX41, STING, phosphorylated TBK1, and phosphorylated p38 MAPK in wild-type BMDCs following treatment with c-di-AMP. Collectively, these data demonstrate that C3a and C5a, via direct signaling through their specific receptors, suppress IFN-β expression by modulation of a distinct innate cytosolic surveillance pathway involving DDX41, STING, and other downstream molecular targets of L. monocytogenes-generated c-di-AMP.

Discovery of functionally selective C5aR2 ligands: novel modulators of C5a signalling

The complement cascade is comprised of a highly sophisticated network of innate immune proteins that are activated in response to invading pathogens or tissue injury. The complement activation peptide, C5a, binds two seven transmembrane receptors, namely the C5a receptor 1 (C5aR1) and C5a receptor 2 (C5aR2, or C5L2). C5aR2 is a non-G-protein-signalling receptor whose biological role remains controversial. Some of this controversy arises owing to the lack of selective ligands for C5aR2. In this study, a library of 61 peptides based on the C-terminus of C5a was assayed for the ability to selectively modulate C5aR2 function. Two ligands (P32 and P59) were identified as functionally selective C5aR2 ligands, exhibiting selective recruitment of β-arrestin 2 via C5aR2, partial inhibition of C5a-induced ERK1/2 activation and lipopolysaccharide-stimulated interleukin-6 release from human monocyte-derived macrophages. Importantly, neither ligand could induce ERK1/2 activation or inhibit C5a-induced ERK1/2 activation via C5aR1 directly. Finally, P32 inhibited C5a-mediated neutrophil mobilisation in wild-type, but not C5aR2(-/-) mice. These functionally selective ligands for C5aR2 are novel tools that can selectively modulate C5a activity in vitro and in vivo, and thus will be valuable tools to interrogate C5aR2 function.

De novo design of a four-fold symmetric TIM-barrel protein with atomic-level accuracy

Despite efforts for over 25 years, de novo protein design has not succeeded in achieving the TIM-barrel fold. Here we describe the computational design of four-fold symmetrical (β/α)8 barrels guided by geometrical and chemical principles. Experimental characterization of 33 designs revealed the importance of side chain-backbone hydrogen bonds for defining the strand register between repeat units. The X-ray crystal structure of a designed thermostable 184-residue protein is nearly identical to that of the designed TIM-barrel model. PSI-BLAST searches do not identify sequence similarities to known TIM-barrel proteins, and sensitive profile-profile searches indicate that the design sequence is distant from other naturally occurring TIM-barrel superfamilies, suggesting that Nature has sampled only a subset of the sequence space available to the TIM-barrel fold. The ability to design TIM barrels de novo opens new possibilities for custom-made enzymes.

Rational, computer-enabled peptide drug design: principles, methods, applications and future directions

Peptides provide promising templates for developing drugs to occupy a middle space between small molecules and antibodies and for targeting 'undruggable' intracellular protein-protein interactions. Importantly, rational or in cerebro design, especially when coupled with validated in silico tools, can be used to efficiently explore chemical space and identify islands of 'drug-like' peptides to satisfy diverse drug discovery program objectives. Here, we consider the underlying principles of and recent advances in rational, computer-enabled peptide drug design. In particular, we consider the impact of basic physicochemical properties, potency and ADME/Tox opportunities and challenges, and recently developed computational tools for enabling rational peptide drug design. Key principles and practices are spotlighted by recent case studies. We close with a hypothetical future case study.

Crystal structure of designed PX domain from cytokine-independent survival kinase and implications on evolution-based protein engineering

The Phox homology domain (PX domain) is a phosphoinositide-binding structural domain that is critical in mediating protein and cell membrane association and has been found in more than 100 eukaryotic proteins. The abundance of PX domains in nature offers an opportunity to redesign the protein using EvoDesign, a computational approach to design new sequences based on structure profiles of multiple evolutionarily related proteins. In this study, we report the X-ray crystallographic structure of a designed PX domain from the cytokine-independent survival kinase (CISK), which has been implicated as functioning in parallel with PKB/Akt in cell survival and insulin responses. Detailed data analysis of the designed CISK-PX protein demonstrates positive impacts of knowledge-based secondary structure and solvation predictions and structure-based sequence profiles on the efficiency of the evolutionary-based protein design method. The structure of the designed CISK-PX domain is close to the wild-type (1.54 Å in Cα RMSD), which was accurately predicted by I-TASSER based fragment assembly simulations (1.32 Å in Cα RMSD). This study represents the first successfully designed conditional peripheral membrane protein fold and has important implications in the examination and experimental validation of the evolution-based protein design approaches.

Is the complement activation product C3a a proinflammatory molecule? Re-evaluating the evidence and the myth

The complement activation product C3a is often described as a proinflammatory mediator, alongside its downstream cousin, C5a. However, emerging studies show that C3a has several anti-inflammatory facets in vivo. For example, in the acute inflammatory response, C3a acts in direct opposition to C5a, through preventing the accumulation of neutrophils in inflamed tissues by independently regulating their mobilization. This acute, protective, and opposing activity of C3a to C5a is also illustrated in models of septicemia. In this article, we reinvestigate the discovery and original classification of C3a as a proinflammatory mediator and highlight the emerging studies demonstrating anti-inflammatory effects for C3a in the immune response. It is our hope that this review illuminates these apparently contradictory roles for C3a and challenges the general dogma surrounding C3a, which, historically, has ubiquitously been described as a proinflammatory mediator. In light of this, we urge investigators to use "inflammatory modulator" as the descriptor for C3a.

New compstatin peptides containing N-terminal extensions and non-natural amino acids exhibit potent complement inhibition and improved solubility characteristics

Compstatin peptides are complement inhibitors that bind and inhibit cleavage of complement C3. Peptide binding is enhanced by hydrophobic interactions; however, poor solubility promotes aggregation in aqueous environments. We have designed new compstatin peptides derived from the W4A9 sequence (Ac-ICVWQDWGAHRCT-NH2, cyclized between C2 and C12), based on structural, computational, and experimental studies. Furthermore, we developed and utilized a computational framework for the design of peptides containing non-natural amino acids. These new compstatin peptides contain polar N-terminal extensions and non-natural amino acid substitutions at positions 4 and 9. Peptides with α-modified non-natural alanine analogs at position 9, as well as peptides containing only N-terminal polar extensions, exhibited similar activity compared to W4A9, as quantified via ELISA, hemolytic, and cell-based assays, and showed improved solubility, as measured by UV absorbance and reverse-phase HPLC experiments. Because of their potency and solubility, these peptides are promising candidates for therapeutic development in numerous complement-mediated diseases.

Forcefield_NCAA: ab initio charge parameters to aid in the discovery and design of therapeutic proteins and peptides with unnatural amino acids and their application to complement inhibitors of the compstatin family

We describe the development and testing of ab initio derived, AMBER ff03 compatible charge parameters for a large library of 147 noncanonical amino acids including β- and N-methylated amino acids for use in applications such as protein structure prediction and de novo protein design. The charge parameter derivation was performed using the RESP fitting approach. Studies were performed assessing the suitability of the derived charge parameters in discriminating the activity/inactivity between 63 analogs of the complement inhibitor Compstatin on the basis of previously published experimental IC50 data and a screening procedure involving short simulations and binding free energy calculations. We found that both the approximate binding affinity (K*) and the binding free energy calculated through MM-GBSA are capable of discriminating between active and inactive Compstatin analogs, with MM-GBSA performing significantly better. Key interactions between the most potent Compstatin analog that contains a noncanonical amino acid are presented and compared to the most potent analog containing only natural amino acids and native Compstatin. We make the derived parameters and an associated web interface that is capable of performing modifications on proteins using Forcefield_NCAA and outputting AMBER-ready topology and parameter files freely available for academic use at http://selene.princeton.edu/FFNCAA . The forcefield allows one to incorporate these customized amino acids into design applications with control over size, van der Waals, and electrostatic interactions.

Derivation of ligands for the complement C3a receptor from the C-terminus of C5a

The complement cascade is a highly sophisticated network of proteins that are well regulated and directed in response to invading pathogens or tissue injury. Complement C3a and C5a are key mediators produced by this cascade, and their dysregulation has been linked to a plethora of inflammatory and autoimmune diseases. Consequently, this has stimulated interest in the development of ligands for the receptors for these complement peptides, C3a receptor, and C5a1 (C5aR/CD88). In this study we used computational methods to design novel C5a1 receptor ligands. However, functional screening in human monocyte-derived macrophages using the xCELLigence label-free platform demonstrated altered specificity of our ligands. No agonist/antagonist activity was observed at C5a1, but we instead saw that the ligands were able to partially agonize the closely related complement receptor C3a receptor. This was verified in the presence of C3a receptor antagonist SB 290157 and in a stable cell line expressing either C5a1 or C3a receptor alone. C3a agonism has been suggested to be a potential treatment of acute neutrophil-driven traumatic pathologies, and may have great potential as a therapeutic avenue in this arena.

De novo design and experimental characterization of ultrashort self-associating peptides

Self-association is a common phenomenon in biology and one that can have positive and negative impacts, from the construction of the architectural cytoskeleton of cells to the formation of fibrils in amyloid diseases. Understanding the nature and mechanisms of self-association is important for modulating these systems and in creating biologically-inspired materials. Here, we present a two-stage de novo peptide design framework that can generate novel self-associating peptide systems. The first stage uses a simulated multimeric template structure as input into the optimization-based Sequence Selection to generate low potential energy sequences. The second stage is a computational validation procedure that calculates Fold Specificity and/or Approximate Association Affinity (K*association) based on metrics that we have devised for multimeric systems. This framework was applied to the design of self-associating tripeptides using the known self-associating tripeptide, Ac-IVD, as a structural template. Six computationally predicted tripeptides (Ac-LVE, Ac-YYD, Ac-LLE, Ac-YLD, Ac-MYD, Ac-VIE) were chosen for experimental validation in order to illustrate the self-association outcomes predicted by the three metrics. Self-association and electron microscopy studies revealed that Ac-LLE formed bead-like microstructures, Ac-LVE and Ac-YYD formed fibrillar aggregates, Ac-VIE and Ac-MYD formed hydrogels, and Ac-YLD crystallized under ambient conditions. An X-ray crystallographic study was carried out on a single crystal of Ac-YLD, which revealed that each molecule adopts a β-strand conformation that stack together to form parallel β-sheets. As an additional validation of the approach, the hydrogel-forming sequences of Ac-MYD and Ac-VIE were shuffled. The shuffled sequences were computationally predicted to have lower K*association values and were experimentally verified to not form hydrogels. This illustrates the robustness of the framework in predicting self-associating tripeptides. We expect that this enhanced multimeric de novo peptide design framework will find future application in creating novel self-associating peptides based on unnatural amino acids, and inhibitor peptides of detrimental self-aggregating biological proteins.

Elucidating a key anti-HIV-1 and cancer-associated axis: the structure of CCL5 (Rantes) in complex with CCR5

CCL5 (RANTES) is an inflammatory chemokine which binds to chemokine receptor CCR5 and induces signaling. The CCL5:CCR5 associated chemotactic signaling is of critical biological importance and is a potential HIV-1 therapeutic axis. Several studies provided growing evidence for the expression of CCL5 and CCR5 in non-hematological malignancies. Therefore, the delineation of the CCL5:CCR5 complex structure can pave the way for novel CCR5-targeted drugs. We employed a computational protocol which is primarily based on free energy calculations and molecular dynamics simulations, and report, what is to our knowledge, the first computationally derived CCL5:CCR5 complex structure which is in excellent agreement with experimental findings and clarifies the functional role of CCL5 and CCR5 residues which are associated with binding and signaling. A wealth of polar and non-polar interactions contributes to the tight CCL5:CCR5 binding. The structure of an HIV-1 gp120 V3 loop in complex with CCR5 has recently been derived through a similar computational protocol. A comparison between the CCL5 : CCR5 and the HIV-1 gp120 V3 loop : CCR5 complex structures depicts that both the chemokine and the virus primarily interact with the same CCR5 residues. The present work provides insights into the blocking mechanism of HIV-1 by CCL5.

Molecular recognition of CCR5 by an HIV-1 gp120 V3 loop

The binding of protein HIV-1 gp120 to coreceptors CCR5 or CXCR4 is a key step of the HIV-1 entry to the host cell, and is predominantly mediated through the V3 loop fragment of HIV-1 gp120. In the present work, we delineate the molecular recognition of chemokine receptor CCR5 by a dual tropic HIV-1 gp120 V3 loop, using a comprehensive set of computational tools predominantly based on molecular dynamics simulations and free energy calculations. We report, what is to our knowledge, the first complete HIV-1 gp120 V3 loop : CCR5 complex structure, which includes the whole V3 loop and the N-terminus of CCR5, and exhibits exceptional agreement with previous experimental findings. The computationally derived structure sheds light into the functional role of HIV-1 gp120 V3 loop and CCR5 residues associated with the HIV-1 coreceptor activity, and provides insights into the HIV-1 coreceptor selectivity and the blocking mechanism of HIV-1 gp120 by maraviroc. By comparing the binding of the specific dual tropic HIV-1 gp120 V3 loop with CCR5 and CXCR4, we observe that the HIV-1 gp120 V3 loop residues 13-21, which include the tip, share nearly identical structural and energetic properties in complex with both coreceptors. This result paves the way for the design of dual CCR5/CXCR4 targeted peptides as novel potential anti-AIDS therapeutics.

Molecular recognition of CXCR4 by a dual tropic HIV-1 gp120 V3 loop

HIV-1 cell entry is initiated by the interaction of the viral envelope glycoprotein gp120 with CD4, and chemokine coreceptors CXCR4 and CCR5. The molecular recognition of CXCR4 or CCR5 by the HIV-1 gp120 is mediated through the V3 loop, a fragment of gp120. The binding of the V3 loop to CXCR4 or CCR5 determines the cell tropism of HIV-1 and constitutes a key step before HIV-1 cell entry. Thus, elucidating the molecular recognition of CXCR4 by the V3 loop is important for understanding HIV-1 viral infectivity and tropism, and for the design of HIV-1 inhibitors. We employed a comprehensive set of computational tools, predominantly based on free energy calculations and molecular-dynamics simulations, to investigate the molecular recognition of CXCR4 by a dual tropic V3 loop. We report what is, to our knowledge, the first HIV-1 gp120 V3 loop:CXCR4 complex structure. The computationally derived structure reveals an abundance of polar and nonpolar intermolecular interactions contributing to the HIV-1 gp120:CXCR4 binding. Our results are in remarkable agreement with previous experimental findings. Therefore, this work sheds light on the functional role of HIV-1 gp120 V3 loop and CXCR4 residues associated with HIV-1 coreceptor activity.

Elucidating a key component of cancer metastasis: CXCL12 (SDF-1α) binding to CXCR4

The chemotactic signaling induced by the binding of chemokine CXCL12 (SDF-1α) to chemokine receptor CXCR4 is of significant biological importance and is a potential therapeutic axis against HIV-1. However, as CXCR4 is overexpressed in certain cancer cells, the CXCL12:CXCR4 signaling is involved in tumor metastasis, progression, angiogenesis, and survival. Motivated by the pivotal role of the CXCL12:CXCR4 axis in cancer, we employed a comprehensive set of computational tools, predominantly based on free energy calculations and molecular dynamics simulations, to obtain insights into the molecular recognition of CXCR4 by CXCL12. We report, what is to our knowledge, the first computationally derived CXCL12:CXCR4 complex structure which is in remarkable agreement with experimental findings and sheds light into the functional role of CXCL12 and CXCR4 residues which are associated with binding and signaling. Our results reveal that the CXCL12 N-terminal domain is firmly bound within the CXCR4 transmembrane domain, and the central 24–50 residue domain of CXCL12 interacts with the upper N-terminal domain of CXCR4. The stability of the CXCL12:CXCR4 complex structure is attributed to an abundance of nonpolar and polar intermolecular interactions, including salt bridges formed between positively charged CXCL12 residues and negatively charged CXCR4 residues. The success of the computational protocol can mainly be attributed to the nearly exhaustive docking conformational search, as well as the heterogeneous dielectric implicit water-membrane-water model used to simulate and select the optimum conformations. We also recently utilized this protocol to elucidate the binding of an HIV-1 gp120 V3 loop in complex with CXCR4, and a comparison between the molecular recognition of CXCR4 by CXCL12 and the HIV-1 gp120 V3 loop shows that both CXCL12 and the HIV-1 gp120 V3 loop share the same CXCR4 binding pocket, as they mostly interact with the same CXCR4 residues.

De novo peptide design and experimental validation of histone methyltransferase inhibitors

Histones are small proteins critical to the efficient packaging of DNA in the nucleus. DNA–protein complexes, known as nucleosomes, are formed when the DNA winds itself around the surface of the histones. The methylation of histone residues by enhancer of zeste homolog 2 (EZH2) maintains gene repression over successive cell generations. Overexpression of EZH2 can silence important tumor suppressor genes leading to increased invasiveness of many types of cancers. This makes the inhibition of EZH2 an important target in the development of cancer therapeutics. We employed a three-stage computational de novo peptide design method to design inhibitory peptides of EZH2. The method consists of a sequence selection stage and two validation stages for fold specificity and approximate binding affinity. The sequence selection stage consists of an integer linear optimization model that was solved to produce a rank-ordered list of amino acid sequences with increased stability in the bound peptide-EZH2 structure. These sequences were validated through the calculation of the fold specificity and approximate binding affinity of the designed peptides. Here we report the discovery of novel EZH2 inhibitory peptides using the de novo peptide design method. The computationally discovered peptides were experimentally validated in vitro using dose titrations and mechanism of action enzymatic assays. The peptide with the highest in vitro response, SQ037, was validated in nucleo using quantitative mass spectrometry-based proteomics. This peptide had an IC₅₀ of 13.5 M, demonstrated greater potency as an inhibitor when compared to the native and K27A mutant control peptides, and demonstrated competitive inhibition versus the peptide substrate. Additionally, this peptide demonstrated high specificity to the EZH2 target in comparison to other histone methyltransferases. The validated peptides are the first computationally designed peptides that directly inhibit EZH2. These inhibitors should prove useful for further chromatin biology investigations.

Forcefield_PTM: Ab Initio Charge and AMBER Forcefield Parameters for Frequently Occurring Post-Translational Modifications

In this work, we introduce Forcefield_PTM, a set of AMBER forcefield parameters consistent with ff03 for 32 common post-translational modifications. Partial charges were calculated through ab initio calculations and a two-stage RESP-fitting procedure in an ether-like implicit solvent environment. The charges were found to be generally consistent with others previously reported for phosphorylated amino acids, and trimethyllysine, using different parameterization methods. Pairs of modified and their corresponding unmodified structures were curated from the PDB for both single and multiple modifications. Background structural similarity was assessed in the context of secondary and tertiary structures from the global dataset. Next, the charges derived for Forcefield_PTM were tested on a macroscopic scale using unrestrained all-atom Langevin molecular dynamics simulations in AMBER for 34 (17 pairs of modified/unmodified) systems in implicit solvent. Assessment was performed in the context of secondary structure preservation, stability in energies, and correlations between the modified and unmodified structure trajectories on the aggregate. As an illustration of their utility, the parameters were used to compare the structural stability of the phosphorylated and dephosphorylated forms of OdhI. Microscopic comparisons between quantum and AMBER single point energies along key χ torsions on several PTMs were performed and corrections to improve their agreement in terms of mean squared errors and squared correlation coefficients were parameterized. This forcefield for post-translational modifications in condensed-phase simulations can be applied to a number of biologically relevant and timely applications including protein structure prediction, protein and peptide design, docking, and to study the effect of PTMs on folding and dynamics. We make the derived parameters and an associated interactive webtool capable of performing post-translational modifications on proteins using Forcefield_PTM available at http://selene.princeton.edu/FFPTM.

Protein folding and de novo protein design for biotechnological applications

In the postgenomic era, the medical/biological fields are advancing faster than ever. However, before the power of full-genome sequencing can be fully realized, the connection between amino acid sequence and protein structure, known as the protein folding problem, needs to be elucidated. The protein folding problem remains elusive, with significant difficulties still arising when modeling amino acid sequences lacking an identifiable template. Understanding protein folding will allow for unforeseen advances in protein design; often referred to as the inverse protein folding problem. Despite challenges in protein folding, de novo protein design has recently demonstrated significant success via computational techniques. We review advances and challenges in protein structure prediction and de novo protein design, and highlight their interplay in successful biotechnological applications.

An evolution-based approach to De Novo protein design and case study on Mycobacterium tuberculosis

Computational protein design is a reverse procedure of protein folding and structure prediction, where constructing structures from evolutionarily related proteins has been demonstrated to be the most reliable method for protein 3-dimensional structure prediction. Following this spirit, we developed a novel method to design new protein sequences based on evolutionarily related protein families. For a given target structure, a set of proteins having similar fold are identified from the PDB library by structural alignments. A structural profile is then constructed from the protein templates and used to guide the conformational search of amino acid sequence space, where physicochemical packing is accommodated by single-sequence based solvation, torsion angle, and secondary structure predictions. The method was tested on a computational folding experiment based on a large set of 87 protein structures covering different fold classes, which showed that the evolution-based design significantly enhances the foldability and biological functionality of the designed sequences compared to the traditional physics-based force field methods. Without using homologous proteins, the designed sequences can be folded with an average root-mean-square-deviation of 2.1 Å to the target. As a case study, the method is extended to redesign all 243 structurally resolved proteins in the pathogenic bacteria Mycobacterium tuberculosis, which is the second leading cause of death from infectious disease. On a smaller scale, five sequences were randomly selected from the design pool and subjected to experimental validation. The results showed that all the designed proteins are soluble with distinct secondary structure and three have well ordered tertiary structure, as demonstrated by circular dichroism and NMR spectroscopy. Together, these results demonstrate a new avenue in computational protein design that uses knowledge of evolutionary conservation from protein structural families to engineer new protein molecules of improved fold stability and biological functionality.

The goal of computational protein design is to create new protein sequences of desirable structure and biological function. Most protein design methods are developed to search for sequences with the lowest free-energy based on physics-based force fields following Anfinsen's thermodynamic hypothesis. A major obstacle of such approaches is the inaccuracy of the force-field design, which cannot accurately describe atomic interactions or correctly recognize protein folds. We propose a novel method which uses evolutionary information, in the form of sequence profiles from structure families, to guide the sequence design. Since sequence profiles are generally more accurate than physics-based potentials in protein fold recognition, a unique advantage lies on that it targets the design procedure to a family of protein sequence profiles to enhance the robustness of designed sequences. The method was tested on 87 proteins and the designed sequences can be folded by I-TASSER to models with an average RMSD 2.1 Å. As a case study of large-scale application, the method is extended to redesign all structurally resolved proteins in the human pathogenic bacteria, Mycobacterium tuberculosis. Five sequences varying in fold and sizes were characterized by circular dichroism and NMR spectroscopy experiments and three were shown to have ordered tertiary structure.

Novel compstatin family peptides inhibit complement activation by drusen-like deposits in human retinal pigmented epithelial cell cultures

We have used a novel human retinal pigmented epithelial (RPE) cell-based model that mimics drusen biogenesis and the pathobiology of age-related macular degeneration to evaluate the efficacy of newly designed peptide inhibitors of the complement system. The peptides belong to the compstatin family and, compared to existing compstatin analogs, have been optimized to promote binding to their target, complement protein C3, and to enhance solubility by improving their polarity/hydrophobicity ratios. Based on analysis of molecular dynamics simulation data of peptide-C3 complexes, novel binding features were designed by introducing intermolecular salt bridge-forming arginines at the N-terminus and at position -1 of N-terminal dipeptide extensions. Our study demonstrates that the RPE cell assay has discriminatory capability for measuring the efficacy and potency of inhibitory peptides in a macular disease environment.

Protein WISDOM: a workbench for in silico de novo design of biomolecules

The aim of de novo protein design is to find the amino acid sequences that will fold into a desired 3-dimensional structure with improvements in specific properties, such as binding affinity, agonist or antagonist behavior, or stability, relative to the native sequence. Protein design lies at the center of current advances drug design and discovery. Not only does protein design provide predictions for potentially useful drug targets, but it also enhances our understanding of the protein folding process and protein-protein interactions. Experimental methods such as directed evolution have shown success in protein design. However, such methods are restricted by the limited sequence space that can be searched tractably. In contrast, computational design strategies allow for the screening of a much larger set of sequences covering a wide variety of properties and functionality. We have developed a range of computational de novo protein design methods capable of tackling several important areas of protein design. These include the design of monomeric proteins for increased stability and complexes for increased binding affinity. To disseminate these methods for broader use we present Protein WISDOM (http://www.proteinwisdom.org), a tool that provides automated methods for a variety of protein design problems. Structural templates are submitted to initialize the design process. The first stage of design is an optimization sequence selection stage that aims at improving stability through minimization of potential energy in the sequence space. Selected sequences are then run through a fold specificity stage and a binding affinity stage. A rank-ordered list of the sequences for each step of the process, along with relevant designed structures, provides the user with a comprehensive quantitative assessment of the design. Here we provide the details of each design method, as well as several notable experimental successes attained through the use of the methods.

Systematic optimization model and algorithm for binding sequence selection in computational enzyme design

A systematic optimization model for binding sequence selection in computational enzyme design was developed based on the transition state theory of enzyme catalysis and graph-theoretical modeling. The saddle point on the free energy surface of the reaction system was represented by catalytic geometrical constraints, and the binding energy between the active site and transition state was minimized to reduce the activation energy barrier. The resulting hyperscale combinatorial optimization problem was tackled using a novel heuristic global optimization algorithm, which was inspired and tested by the protein core sequence selection problem. The sequence recapitulation tests on native active sites for two enzyme catalyzed hydrolytic reactions were applied to evaluate the predictive power of the design methodology. The results of the calculation show that most of the native binding sites can be successfully identified if the catalytic geometrical constraints and the structural motifs of the substrate are taken into account. Reliably predicting active site sequences may have significant implications for the creation of novel enzymes that are capable of catalyzing targeted chemical reactions.

Anaphylatoxins in organ transplantation

C3a and C5a (also called anaphylatoxins) are inflammatory peptides generated during complement activation. They do not only play important roles in innate immunity through the initiation and regulation of inflammatory responses, but also significantly influence adaptive immune responses. Organ transplantation triggers an initial inflammatory response and subsequent to the specific immune response (also called the alloimmune response), both of which contribute to graft rejection. Emerging evidence suggests that anaphylatoxins, particularly C5a, are significantly involved in both inflammatory and alloimmune responses following organ transplantation, thus influencing graft outcome. This review will provide the information on our current understanding of the roles for anaphylatoxins in ischemia-reperfusion injury, graft rejection, and transplant tolerance, and the therapeutic potential of targeting anaphylatoxin receptors in organ transplantation.

The receptor for complement component C3a mediates protection from intestinal ischemia-reperfusion injuries by inhibiting neutrophil mobilization

C3a is a key complement activation fragment, yet its neutrophil-expressed receptor (C3aR) still has no clearly defined role. In this study, we used a neutrophil-dependent mouse model of intestinal ischemia-reperfusion (IR) injury to explore the role of C3aR in acute tissue injuries. C3aR deficiency worsened intestinal injury, which corresponded with increased numbers of tissue-infiltrating neutrophils. Circulating neutrophils were significantly increased in C3aR(-/-) mice after intestinal ischemia, and C3aR(-/-) mice also mobilized more circulating neutrophils after granulocyte colony-stimulating factor infusion compared with WT mice, indicating a specific role for C3aR in constraining neutrophil mobilization in response to intestinal injury. In support of this role, C3aR(-/-) mice reconstituted with WT bone marrow reversed IR pathology back to WT levels. Complement C5a receptor (C5aR) antagonism in C3aR(-/-) mice also rectified the worsened pathology after intestinal IR injury but had no effect on circulating neutrophils, highlighting the opposing roles of C3a and C5a in disease pathogenesis. Finally, we found that using a potent C3a agonist to activate C3aR in vivo reduced neutrophil mobilization and ameliorated intestinal IR pathology in WT, but not C3aR(-/-), mice. This study identifies a role for C3aR in regulating neutrophil mobilization after acute intestinal injury and highlights C3aR agonism as a potential treatment option for acute, neutrophil-driven pathologies.