Search

The complement system is assembled from a network of proteins that function to bring about the body’s first line of defense against invading pathogens. However, complement deficiencies or invasive pathogens can hijack complement to subsequently increase the body’s susceptibility to infections. What’s more, invasive pathogens are increasingly becoming resistant to the currently available therapies. Hence, it is important to gain insights into the highly dynamic interaction between complement and invading microbes in the frontlines of immunity. Here, we developed a computational model of the complement system composed of 670 ordinary differential equations with 328 kinetic parameters, which describes all three complement pathways (alternative, classical, and lectin) and includes description of mannose-binding lectin, collectins, ficolins, factor H-related proteins, immunoglobulin M, and pentraxins. Additionally, we incorporate two pathogens: (type 1) complement susceptible pathogen, (type 2) Neisseria meningitidis located in either nasopharynx or bloodstream. In both cases, we generate time profiles of the pathogen surface occupied by complement components and the membrane attack complex (MAC). Our model shows both pathogen types in bloodstream are saturated by complement proteins whereas MACs occupy <

Central nervous system (CNS) injury and infection can result in profound tissue remodeling in the brain, the mechanism and purpose of which is poorly understood. Infection with the protozoan parasite Toxoplasma gondii causes chronic infection and inflammation in the brain parenchyma. Control of parasite replication requires the continuous presence of IFNγ-producing T cells to keep T. gondii in its slowly replicating cyst form. During infection, a network of extracellular matrix fibers, revealed using multiphoton microscopy, forms in the brain. The origin and composition of these structures are unknown but the fibers have been observed to act as a substrate for migrating T cells. In this study, we show a critical regulator of extracellular matrix (ECM) remodeling, Secreted Protein, Acidic, Rich in Cysteine (SPARC), is upregulated in the brain during the early phases of infection in the frontal cortex. In the absence of SPARC, a reduced and disordered fibrous network, increased parasite burden, and reduced antigen-specific T cell entry into the brain points to a role for SPARC in T cell recruitment to and migration within the brain. We also report SPARC can directly bind to CCR7 ligands CCL19 and CCL21 but not CXCL10, and enhance migration toward a chemokine gradient. Measurement of T cell behavior points to tissue remodeling being important for access of immune cells to the brain and facilitating cellular locomotion. Together, these data identify SPARC as an important regulatory component of immune cell trafficking and access to the inflamed CNS.

Poxviruses are dangerous pathogens, which can cause fatal infection in unvaccinated individuals. The causative agent of smallpox in humans, variola virus, is closely related to the bovine vaccinia virus, yet the molecular basis of their selectivity is currently incompletely understood. Here, we examine the role of the electrostatics in the selectivity of the smallpox protein SPICE and vaccinia protein VCP toward the human and bovine complement protein C3b, a key component of the complement immune response. Electrostatic calculations, in-silico alanine-scan and electrostatic hotspot analysis, as introduced by Kieslich and Morikis (PLoS Comput. Biol. 2012), are used to assess the electrostatic complementarity and to identify sites resistant to local perturbation where the electrostatic potential is likely to be evolutionary conserved. The calculations suggest that the bovine C3b is electrostatically prone to selectively bind its VCP ligand. On the other hand, the human isoform of C3b exhibits a lower electrostatic complementarity toward its SPICE ligand. Yet, the human C3b displays a highly preserved electrostatic core, which suggests that this isoform could be less selective in binding different ligands like SPICE and the human Factor H. This is supported by experimental cofactor activity assays revealing that the human C3b is prone to bind both SPICE and Factor H, which exhibit diverse electrostatic properties. Additional investigations considering mutants of SPICE and VCP that revert their selectivity reveal an “electrostatic switch” into the central modules of the ligands, supporting the critical role of the electrostatics in the selectivity. Taken together, these evidences provide insights into the selectivity mechanism of the complement regulator proteins encoded by the variola and vaccinia viruses to circumvent the complement immunity and exert their pathogenic action. These fundamental aspects are valuable for the development of novel vaccines and therapeutic strategies.

© 2020 Bernardino de la Serna, Mellado, Dustin, Garcia-Parajo and Morikis., The immune system comprises a collection of specialized cells, tissues, and organs that protect the organisms against pathogens and can survey cancer cells. Immune responses are precisely coordinated events that take place in complex, specialized tissue microenvironments. For an integrated view of innate and adaptive immune responses at the molecular level, we ideally need a better understanding of how immune cells communicate and fulfill their tasks in vivo, following events spatially and temporally. Conventional biochemical and genetic methods consider the cell as an individual entity and ligand/receptor pairs as isolated systems. Often, the data obtained refers to the average behavior of a pool of cells and/or receptors removed from their real-life context. The use of new technologies, particularly real-time imaging approaches, is showing us that biological responses are very dynamic and extremely dependent on the context in which they take place and are therefore much more diverse than we initially thought. The combination of these new approaches is radically transforming and enriching immunology, as demonstrated by the increasing number of publications in which physical and/or engineering tools are applied to study the immune response. Whilst scientists are often questioned for the discipline their research is best framed in, we rather think that one scientific discipline cannot be reduced to the terms of another. However, defining and naming cross-disciplinary fields sets our minds on common ground and helps establish a fluent communication to eventually produce groundbreaking, beautiful pieces of science. For instance, ImmunoPhysics was probably first coined by Prof. Morikis a couple of decades ago (https://www.biophysics.org/profiles/dimitrios-morikis); nevertheless, ImmunoPhysics has not become widely regarded as a discipline, despite the continuously growing body of research that requires physical approaches to resolving immunological questions. Hence, with this special issue, we wanted to open a scientific platform compiling ImmunoPhysics and ImmunoEngineering research breakthroughs and future perspectives. Sadly, towards the end of this fascinating journey Prof. Morikis passed away; thus now, with this special issue, we would also like to pay tribute to his fundamental contributions to the field.

A single nucleotide polymorphism, tyrosine at position 402 to histidine (Y402H), within the gene encoding complement factor H (FH) predisposes individuals to acquiring age-related macular degeneration (AMD) after aging. This polymorphism occurs in short consensus repeat (SCR) 7 of FH and results in decreased binding affinity of SCR6-8 for heparin. As FH is responsible for regulating the complement system, decreased affinity for heparin results in decreased regulation on surfaces of self. To understand the involvement of the Y402H polymorphism in AMD, we leverage methods from bioinformatics and computational biophysics to quantify structural and dynamical differences between SCR7 isoforms that contribute to decreased pattern recognition in SCR7(H402). Our data from molecular and Brownian dynamics simulations suggest a revised mechanism for decreased heparin binding. In this model, transient contacts not observed in structures for SCR7 are predicted to occur in molecular dynamics simulations between coevolved residues Y402 and I412, stabilizing SCR7(Y402) in a conformation that promotes association with heparin. H402 in the risk isoform is less likely to form a contact with I412 and samples a larger conformational space than Y402. We observe energy minima for sidechains of Y402 and R404 from SCR7(Y402) that are predicted to associate with heparin at a rate constant faster than energy minima for sidechains of H402 and R404 from SCR7(H402). As both carbohydrate density and degree of sulfation decrease with age in Bruch’s membrane of the macula, the decreased heparin recognition of SCR7(H402) may contribute to the pathogenesis of AMD.

The complement system is our first line of defense against foreign pathogens, but when it is not properly regulated, complement is implicated in the pathology of several autoimmune and inflammatory disorders. Compstatin is a peptidic complement inhibitor that acts by blocking the cleavage of complement protein C3 to the proinflammatory fragment C3a and opsonin fragment C3b. In this study, we aim to identify druglike small-molecule complement inhibitors with physicochemical, geometric, and binding properties similar to those of compstatin. We employed two approaches using various high-throughput virtual screening methods, which incorporate molecular dynamics (MD) simulations, pharmacophore model design, energy calculations, and molecular docking and scoring. We have generated a library of 274 chemical compounds with computationally predicted binding affinities for the compstatin binding site of C3. We have tested subsets of these chemical compounds experimentally for complement inhibitory activity, using hemolytic assays, and for binding affinity, using microscale thermophoresis. As a result, although none of the compounds showed inhibitory activity, compound 29 was identified to exhibit weak competitive binding against a potent compstatin analogue, therefore validating our computational approaches. Additional docking and MD simulation studies suggest that compound 29 interacts with C3 residues, which have been shown to be important in binding of compstatin to the C3c fragment of C3. Compound 29 is amenable to physicochemical optimization to acquire inhibitory properties. Additionally, it is possible that some of the untested compounds will demonstrate binding and inhibition in future experimental studies.

C3b, the central component of the alternative pathway (AP) of the complement system, coexists as a mixture of conformations in solution. These conformational changes can affect interactions with other proteins and complement regulators. Here we combine a computational model for electrostatic interactions within C3b with molecular imaging to study the conformation of C3b. The computational analysis shows that the TED domain in C3b is tethered ionically to the macroglobulin (MG) ring. Monovalent counterion concentration affects the magnitude of electrostatic forces anchoring the TED domain to the rest of the C3b molecule in a thermodynamic model. This is confirmed by observing NaCl concentration dependent conformational changes using single molecule electron microscopy (EM). We show that the displacement of the TED domain is compatible with C3b binding to Factor B (FB), suggesting that the regulation of the C3bBb convertase could be affected by conditions that promote movement in the TED domain. Our molecular model also predicts mutations that could alter the positioning of the TED domain, including the common R102G polymorphism, a risk variant for developing age-related macular degeneration. The common C3b isoform, C3bS, and the risk isoform, C3bF, show distinct energetic barriers to displacement in the TED that are related to a network of electrostatic interactions at the interface of the TED and MG-ring domains of C3b. These computational predictions agree with experimental evidence that shows differences in conformation observed in C3b isoforms purified from homozygous donors. Altogether, we reveal an ionic, reversible attachment of the TED domain to the MG ring that may influence complement regulation in some mutations and polymorphisms of C3b. Work in this report has been funded by the Spanish Ministry of Economy and Competitiveness (SAF2011-26583 to SRC and SAF2014-52301-R to OL), the Fundación Renal Iñigo Alvarez de Toledo and the Seventh Framework Programme European Union Project EURenOmics (305608) to SRC. In addition, this work has been supported by a grant from the Autonomous Region of Madrid (S2010/BMD-2316) to SRC and OL.

[Image: see text] C3d is a hallmark protein of the complement system, whose presence is critical to measure the progression of several immune diseases. Here, we propose to directly target C3d through small peptides mimicking the binding of its natural ligand, the complement regulator Factor H (FH). Through iterative computational analysis and binding affinity experiments, we establish a rationale for the structure-based design of FH-inspired peptides, leading to low-micromolar affinity for C3d and stable binding over microsecond-length simulations. Our FH-inspired peptides call now for further optimization toward high-affinity binding and suggest that small peptides are promising as novel C3d biomarkers and therapeutic tools.

Background Engineered vaccine proteins incorporating both antigen and adjuvant components are constructed with the aim of combining functions to induce effective protective immunity. Bacterial flagellin is a strong candidate for an engineered vaccine scaffold as it is known to provide adjuvant activity through its TLR5 and inflammasome activation. Moreover, polymerized flagellin filaments can elicit a more robust immunoglobulin response than monomeric flagellin, and the multimeric antigen form can also promote T cell-independent antibody responses. Here, we aim to produce and test a covalently stabilized polymerized flagellar filament, providing additional immune efficacy through stabilization of its polymeric filament structure, as well as stabilization for long-term storage. Results Computational modeling of monomer packing in flagellin filaments helped identify amino acids with proximity to neighboring flagella protofilaments. Paired cysteine substitutions were made at amino acids predicted to form inter-monomer disulfide cross-links, and these substitutions were capable of forming flagella when transfected into a flagellin-negative strain of Salmonella enterica subspecies Typhimurium. Interestingly, each paired substitution stabilized different helical conformational polymorphisms; the stabilized filaments lost the ability to transition between conformations, reducing bacterial motility. More importantly, the paired substitutions enabled extensive disulfide cross links and intra-filament multimer formation, and in one of the three variants, permitted filament stability in high acidic and temperature conditions where wild-type filaments would normally rapidly depolymerize. In addition, with regard to potential adjuvant activity, all crosslinked flagella filaments were able to induce wild-type levels of epithelial NF-κB in a cell reporter system. Finally, bacterial virulence was unimpaired in epithelial adherence and invasion, and the cysteine substitutions also appeared to increase bacterial resistance to oxidizing and reducing conditions. Conclusions We identified amino acid pairs, with cysteine substitutions, were able to form intermolecular disulfide bonds that stabilized the resulting flagellar filaments in detergent, hydrochloric acid, and high temperatures while retaining its immunostimulatory function. Flagellar filaments with disulfide-stabilized protofilaments introduce new possibilities for the application of flagella as a vaccine adjuvant. Specifically, increased stability and heat tolerance permits long-term storage in a range of temperature environments, as well as delivery under a range of clinical conditions. Electronic supplementary material The online version of this article (10.1186/s12896-019-0545-3) contains supplementary material, which is available to authorized users.

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.

Electrostatic effects are ubiquitous in protein interactions and are found to be pervasive in the complement system as well. The interaction between complement fragment C3d and complement receptor 2 (CR2) has evolved to become a link between innate and adaptive immunity. Electrostatic interactions have been suggested to be the driving factor for the association of the C3d:CR2 complex. In this study, we investigate the effects of ionic strength and mutagenesis on the association of C3d:CR2 through Brownian dynamics simulations. We demonstrate that the formation of the C3d:CR2 complex is ionic strength-dependent, suggesting the presence of long-range electrostatic steering that accelerates the complex formation. Electrostatic steering occurs through the interaction of an acidic surface patch in C3d and the positively charged CR2 and is supported by the effects of mutations within the acidic patch of C3d that slow or diminish association. Our data are in agreement with previous experimental mutagenesis and binding studies and computational studies. Although the C3d acidic patch may be locally destabilizing because of unfavorable Coulombic interactions of like charges, it contributes to the acceleration of association. Therefore, acceleration of function through electrostatic steering takes precedence to stability. The site of interaction between C3d and CR2 has been the target for delivery of CR2-bound nanoparticle, antibody, and small molecule biomarkers, as well as potential therapeutics. A detailed knowledge of the physicochemical basis of C3d:CR2 association may be necessary to accelerate biomarker and drug discovery efforts.

The G protein-coupled receptor (GPCR) CXCR4 is a co-receptor for HIV and is involved in cancers and autoimmune diseases. We characterized five purine or quinazoline core polyamine pharmacophores used for targeting CXCR4 dysregulation in diseases. All were neutral antagonists for wild-type (WT) CXCR4 and two were biased antagonists with effects on β-arrestin-2 only at high concentrations. These compounds displayed various activities for a constitutively active mutant (CAM). We use the IT1t-CXCR4 crystal structure and molecular dynamics (MD) simulations to develop two hypotheses for the activation of the N119(3.35)A CAM. The N119(3.35)A mutation facilitates increased coupling of TM helices III and VI. IT1t deactivates the CAM by disrupting the coupling between TM helices III and VI, mediated primarily by residue F87(2.53). Mutants of F87(2.53) in N119(3.35)A CXCR4 precluded constitutive signaling and prevented inverse agonism. This work characterizes CXCR4 ligands and provides a mechanism for N119(3.35)A constitutive activation.

The complex between complement system proteins C5b and C6 is the cornerstone for the assembly of the membrane attack complex (MAC, also known as C5b6789n). MAC is the terminal product of three converging pathways of the complement system and functions as a pore forming complex on cell surfaces, as a response of the immune system in fighting pathogens. However, when proper regulation of the complement system is compromised, MAC also attacks host tissues and contributes to several complement-mediated autoimmune diseases. We performed a molecular dynamics and electrostatics study to elucidate the mechanism of interaction between C5b and C6 and the formation of the C5b6 complex. The C5b-C6 interface consists of three binding sites stabilized predominantly by van der Waals interactions, and several critical salt bridges and hydrogen bonds. We discuss differences between domains C5d and C3d that lead to mono-functionality of C5d in acting as the scaffold for MAC formation, as opposed to dual functionality of C3d in acting as an opsonin for phagocytosis and as a link between innate and adaptive immunity, based on a comparative sequence, structural, and physicochemical analysis. We also extended our analysis to pathway dynamics to demonstrate the significance of consumption-production rates of C5b, C6 and C5b6 that lead towards MAC formation. Finally, we propose that C5d is a target for drug discovery, aiming to the inhibition of the MAC formation in autoimmune diseases originating from MAC-mediated host cell lysis.

Despite the similar enzyme cascade in the Ubiquitin and Ubiquitin-like peptide(Ubl) conjugation, the involvement of single or heterodimer E1 activating enzyme has been a mystery. Here, by using a quantitative Förster Resonance Energy Transfer (FRET) technology, aided with Analysis of Electrostatic Similarities Of Proteins (AESOP) computational framework, we elucidate in detail the functional properties of each subunit of the E1 heterodimer activating-enzyme for NEDD8, UBA3 and APPBP1. In contrast to SUMO activation, which requires both subunits of its E1 heterodimer AOS1-Uba2 for its activation, NEDD8 activation requires only one of two E1 subunits, UBA3. The other subunit, APPBP1, only contributes by accelerating the activation reaction rate. This discovery implies that APPBP1 functions mainly as a scaffold protein to enhance molecular interactions and facilitate catalytic reaction. These findings for the first time reveal critical new mechanisms and a potential evolutionary pathway for Ubl activations. Furthermore, this quantitative FRET approach can be used for other general biochemical pathway analysis in a dynamic mode.

The spherical form (s-form) of a genetically-modified gold-binding M13 bacteriophage was investigated as a scaffold for gold synthesis. Repeated mixing of the phage with chloroform caused a 15-fold contraction from a nearly one micron long filament to an approximately 60 nm diameter spheroid. The geometry of the viral template and the helicity of its major coat protein were monitored throughout the transformation process using electron microscopy and circular dichroism spectroscopy, respectively. The transformed virus, which retained both its gold-binding and mineralization properties, was used to assemble gold colloid clusters and synthesize gold nanostructures. Spheroid-templated gold synthesis products differed in morphology from filament-templated ones. Spike-like structures protruded from the spherical template while isotropic particles developed on the filamentous template. Using inductively coupled plasma-mass spectroscopy (ICP-MS), gold ion adsorption was found to be comparatively high for the gold-binding M13 spheroid, and likely contributed to the dissimilar gold morphology. Template contraction was believed to modify the density, as well as the avidity of gold-binding peptides on the scaffold surface. The use of the s-form of the M13 bacteriophage significantly expands the templating capabilities of this viral platform and introduces the potential for further morphological control of a variety of inorganic material systems.

Complement activation plays a major role in many acute and chronic inflammatory conditions. C3d, a terminal product of complement activation, remains covalently attached to cells and is an excellent biomarker of complement-mediated inflammation. We employed a virtual high-throughput screening protocol to identify molecules with predicted binding to complement C3d and with intrinsic fluorescence properties to enable detection. Pharmacophore models were developed based on known C3d-ligand interactions and information from computational analysis of structural and molecular dynamics data. Iterative pharmacophore-based virtual screening was performed to identify druglike molecules with physicochemical similarity to the natural C3d ligand CR2. Hits from the pharmacophore screens were docked to C3d and ranked based on predicted binding free energies. Top-ranked molecules were selected for experimental validation of binding affinity to C3d, using microscale thermophoresis, and for their suitability to become molecular imaging agents, using fluorescence spectroscopy. This work serves as a foundation for identifying additional fluorescent molecules with high-affinity for C3d that will subsequently be explored as noninvasive in vivo diagnostics of complement-mediated inflammation, for spatiotemporal monitoring of disease progression, and for targeting therapeutics to sites of inflammation.

The C3d:CR2(SCR1-2) interaction plays an important role in bridging innate and adaptive immunity, leading to enhanced antibody production at sites of complement activation. Over the past decade, there has been much debate over the binding mode of this interaction. An initial cocrystal structure (PDB: 1GHQ) was published in 2001, in which the only interactions observed were between the SCR2 domain of CR2 and a side-face of C3d whereas a cocrystal structure (PDB: 3OED) published in 2011 showed both the SCR1 and SCR2 domains of CR2 interacting with an acidic patch on the concave surface of C3d. The initial 1GHQ structure is at odds with the majority of existing biochemical data and the publication of the 3OED structure renewed uncertainty regarding the physiological relevance of 1GHQ, suggesting that crystallization may have been influenced by the presence of zinc acetate in the crystallization process. In our study, we used a variety of computational approaches to gain insight into the binding mode between C3d and CR2 and demonstrate that the binding site at the acidic patch (3OED) is electrostatically more favorable, exhibits better structural and dissociative stability, specifically at the SCR1 domain, and has higher binding affinity than the 1GHQ binding mode. We also observe that nonphysiological zinc ions enhance the formation of the C3d:CR2 complex at the side face of C3d (1GHQ) through increases in electrostatic favorability, intermolecular interactions, dissociative character and overall energetic favorability. These results provide a theoretical basis for the association of C3d:CR2 at the acidic cavity of C3d and provide an explanation for binding of CR2 at the side face of C3d in the presence of nonphysiological zinc ions.

Author(s): Winicki, Nolan; Harrison, Reed; Perry, Jefferson; Morikis, Dimitrios | Abstract: Cyclin-dependent protein kinase regulatory subunit 1 (Cks1) is involved in cell cycle progressionthrough interactions with cyclin-dependent kinases (CDK) and ubiquitination of cyclin-dependentkinase inhibitors (CKI). Dysfunction of CDK dependent associations can affect the entrance of acell into mitosis, particularly the G1-S phase transition. Abnormal assistance from Cks1 with themultiprotein complex SCF (Skp2) in the ubiquitination of CDKN1B (p27Kip1) can disrupt the mitoticregulatory protein levels escalating to cancer development. In this study, we use computationalmethods to investigate interactions in three different complexes sharing Cks1 in order to create atargeted pharmaceutical solution that will assist in regulation. Our analysis is performed using thecrystal structure of Cks1 in three complexes to include involvement between one ubiquitin ligaseand two CDKs. The assessment is based on the intermolecular electrostatic interactions, such ashydrogen bonds and charge-charge interactions. We observe that charge and hydrogen bonding playsa significant role in the stability between Cks1 and adjacent proteins in each complex. Due to largebinding interfaces and varied distribution of charge across the main contact regions, we decided thata condensed pocket on Cks1 that interacts with phosphorylated p27Kip1 should be selected for furtherin-depth examination.

Electric fields often play a role in guiding the association of protein complexes. Such interactions can be further engineered to accelerate complex association, resulting in protein systems with increased productivity. This is especially true for enzymes where reaction rates are typically diffusion limited. To facilitate quantitative comparisons of electrostatics in protein families and to describe electrostatic contributions of individual amino acids, we previously developed a computational framework called AESOP. We now implement this computational tool in Python with increased usability and the capability of performing calculations in parallel. AESOP utilizes PDB2PQR and Adaptive Poisson-Boltzmann Solver to generate grid-based electrostatic potential files for protein structures provided by the end user. There are methods within AESOP for quantitatively comparing sets of grid-based electrostatic potentials in terms of similarity or generating ensembles of electrostatic potential files for a library of mutants to quantify the effects of perturbations in protein structure and protein-protein association.

SUMOylation, one of the most important protein post-translational modifications, plays critical roles in a variety of physiological and pathological processes. SENP (Sentrin/SUMO-specific protease), a family of SUMO-specific proteases, is responsible for the processing of pre-SUMO and removal of SUMO from conjugated substrates. SUMO4, the latest discovered member in the SUMO family, has been found as a type 1 diabetes susceptibility gene and its maturation is not understood so far. Despite the 14 amino acid differences between pre-SUMO4 and SUMO2, pre-SUMO4 is not processed by SENP2 but pre-SUMO2 does. A novel interdisciplinary approach involving computational modeling and a FRET-based protease assay was taken to engineer pre-SUMO4 as a substrate of SENP2. Given the difference in net charge between pre-SUMO4 and pre-SUMO2, the computational framework analysis of electrostatic similarities of proteins was applied to determine the contribution of each ionizable amino acid in a model of SENP2-(pre-SUMO4) binding, and to propose pre-SUMO4 mutations. The specificities of the SENP2 toward different pre-SUMO4 mutants were determined using a quantitative FRET assay by characterizing the catalytic efficiencies (kcat/KM). A single amino acid mutation made pre-SUMO4 amenable to SENP2 processing and a combination of two amino acid mutations made it highly accessible as SENP2 substrate. The combination of the two approaches provides a powerful protein engineering tool for future SUMOylation studies.

Despite its sequence variability and structural flexibility, the V3 loop of the HIV-1 envelope glycoprotein gp120 is capable of recognizing cell-bound coreceptors CCR5 and CXCR4 and infecting cells. Viral selection of CCR5 is associated with the early stages of infection, and transition to selection of CXCR4 indicates disease progression. We have developed a predictive statistical model for coreceptor selectivity that uses the discrete property of net charge and the binary coreceptor preference markers of the N(6)X(7)[T/S](8)X(9) glycosylation motif and 11/24/25 positive amino acid rule. The model is based on analysis of 2,054 V3 loop sequences from patient data and allows us to infer the most likely state of the disease from physicochemical characteristics of the sequences. The performance of the model is comparable to established sequence-based predictive methods, and may be used in combination with other methods as a supportive diagnostic for coreceptor selection. This model may be used for personalized medical decisions in administering coreceptor-specific therapies.

Ligand-binding to G protein-coupled receptors (GPCRs) acts as the local driving force that initiates signal transduction through the receptor and mediates its conformational transitions and interactions with various intracellular effectors. In a recent study, We have shown that the binding of ligands CCL19 and CCL21 to CCR7 induces biased triggering of side chain-based molecular switches, which coordinate concerted transmembrane helical domain motions and transitioning of the receptor to distinct conformational states (Gaieb, Z., D.D. Lo, and D. Morikis. 2016. Molecular Mechanism of Biased Ligand Conformational Changes in CC Chemokine Receptor 7. Journal of Chemical Information and Modeling. 56: 1808–1822, DOI: 10.1021/acs.jcim.6b00367). To complement our previous study, we compare the results of the free (apo) CCR7 microsecond molecular dynamics simulations to those of the ligand-bound CCR7, and show that the apo receptor is found in conformational heterogeneity that only exhibits random fluctuations and lacks the coordinated helical motions seen in ligand-bound receptors. We conclude that ligand binding is responsible for coordinating the stochastic conformational nature of CCR7 into specific conformational states, initiated and propagated by specific physicochemical events.

Targeting the complement component 3a receptor (C3aR) with selective agonists or antagonists is believed to be a viable therapeutic option for several diseases such as stroke, heart attack, reperfusion injuries, and rheumatoid arthritis. We designed a number of agonists, partial agonists, and antagonists of C3aR using our two-stage de novo protein design framework. Of the peptides tested using a degranulation assay in C3aR-transfected rat basophilic leukemia cells, two were prominent agonists (EC(50) values of 25.3 and 66.2 nM) and two others were partial agonists (IC(50) values of 15.4 and 26.1 nM). Further testing of these lead compounds in a calcium flux assay in U937 cells yielded similar results although with reduced potencies compared to transfected cells. The partial agonists also displayed full antagonist activity when tested in a C3aR inhibition assay. In addition, the electrostatic potential profile was shown to potentially discriminate between full agonists and partial agonists.

We report the computational and rational design of new generations of potential peptide-based inhibitors of the complement protein C3 from the compstatin family. The binding efficacy of the peptides is tested by extensive molecular dynamics-based structural and physicochemical analysis, using 32 atomic detail trajectories in explicit water for 22 peptides bound to human, rat or mouse target protein C3, with a total of 257 ns. The criteria for the new design are: (i) optimization for C3 affinity and for the balance between hydrophobicity and polarity to improve solubility compared to known compstatin analogs; and (ii) development of dual specificity, human-rat/mouse C3 inhibitors, which could be used in animal disease models. Three of the new analogs are analyzed in more detail as they possess strong and novel binding characteristics and are promising candidates for further optimization. This work paves the way for the development of an improved therapeutic for age-related macular degeneration, and other complement system-mediated diseases, compared to known compstatin variants.

The complement system is a component of innate immunity and is activated by a cascade of protein interactions whose function is vital to our ability to fight infection. When proper regulation fails, the complement system is unable to recognize "self" from "nonself" and, therefore, attacks own tissues leading to autoimmune diseases. The central protein of the complement system is C3, which is the convergence point of three independently activated but communicating pathways. Regulation of C3 occurs through modular proteins which consist of many repeats of complement control protein (CCP) modules. CCP modules have diverse sequences, similar structures, and diverse physicochemical compositions, with excess of charge being a predominant characteristic. The goal of our study is to understand the electrostatic mechanism that underlies the interaction between the C3d domain of C3 and the fourth module of the complement regulator Factor H (FH4). We have performed a computational alanine scan in which we have replaced every ionizable amino acid, one at a time, with an alanine to generate a family of mutants for the C3d-FH4 complex. We have used Poisson-Boltzmann electrostatic calculations in combination with clustering of spatial distributions of electrostatic potentials and free energy calculations to delineate the contribution of each replaced amino acid to the C3d-FH4 interaction. We have analyzed our data in view of a two-step model which separates association into long-range recognition and short-range binding and we have identified key amino acids that contribute to association. We discuss the complex role of C3d in binding FH4 and the bacterial proteins Efb/Ehp from Staphylococcus aureus.

The peptide compstatin and its derivatives inhibit the complement-component protein C3 in primate mammals and are potential therapeutic agents against the unregulated activation of complement in humans, but are inactive against C3 from lower mammals. Recent molecular dynamics (MD) simulations showed that the most potent compstatin analog comprised entirely of natural amino acids (W4A9) had a smaller affinity for rat C3, due to reproducible changes in the rat protein structure with respect to the human protein, which eliminated or weakened specific protein-ligand interactions seen in the human C3:W4A9 complex. Here, we study by MD simulations three W4A9 complexes with the mouse C3 protein, and two "transgenic" mouse derivatives, containing a small number (6-9) of human C3 substitutions. The mouse complex experiences the conformational changes and affinity reduction of the rat complex. In the "transgenic" complexes, the conformation remains closer to that of the human complex, the protein-ligand interactions are improved, and the affinity for compstatin becomes "human-like." The present work creates new avenues for a compstatin-sensitive animal model. A similar strategy, involving the comparison of a series of complexes by MD simulations, could be used to design "transgenic" sequences in other systems.

Virulence factors EfbC and Ehp from Staphylococcus aureus are potent inhibitors of complement activation. Both are excessively charged and bind to complement protein C3d at an acidic interface. We computationally generated single-alanine mutants of charged residues in the C3d–EfbC and C3d–Ehp complexes, and utilized electrostatic clustering and Poisson–Boltzmann free energy calculations to evaluate the role of electrostatics in association. Our results indicate that both interfacial electrostatic interactions and electrostatic potential distribution are crucial for C3d–EfbC and C3d–Ehp association. The results presented herein serve as a predictive tool in the selection of mutants with desired binding and immunological activity, and will narrow the search for viable candidates for further computational and experimental analyses. This study serves as the foundation for development of an inhibitor of the C3d–EfbC and C3d–Ehp interactions to combat bacterial infection and design of a complement inhibitor using EfbC and Ehp as a model.

Molecular similarity is an important tool in protein and drug design for analyzing the quantitative relationships between physicochemical properties of two molecules. We present a family of similarity measures which exploits the ability of wavelet transformation to analyze the spectral components of physicochemical properties and suggests a sensitive way for measuring similarities of biological molecules. In order to investigate how effective wavelet-based similarity measures were against conventional measures, we defined several patterns which involve scalar or topological changes in the distribution of electrostatic properties. The wavelet-based measures were more successful in discriminating these patterns in contrast to the current state-of-art similarity measures. We also present the validity of wavelet-based similarity measures through the hierarchical clustering of two protein datasets consisting of families of homologous domains and alanine scan mutants. This type of similarity analysis is useful for protein structure-function studies and protein design.

Factor H (FH) contributes to the regulation of the complement system by binding to polyanionic surfaces and the proteins C3b/C3c/C3d. This implicates charge and electrostatic interactions in recognition and binding of FH. Despite the large amount of experimental and pathology data the exact mechanism at molecular level is not yet known. We have implemented a computational framework for comparative analysis of the charge and electrostatic diversity of FH modules and C3b domains to identify electrostatic hotspots and predict potential binding sites. Our electrostatic potential clustering analysis shows that charge distributions and electrostatic potential distributions are more useful in understanding C3b–FH interactions than net charges alone. We present a model of non-specific electrostatic interactions of FH with polyanion-rich surfaces and specific interactions with C3b, using our computational data and existing experimental data. We discuss the electrostatic contributions to the formation of the C3b–FH complex and the competition between FH and Factor Bb (Bb) for binding to C3b. We also discuss the significance of mutations of charged amino acids in the pathobiology of FH-mediated disease, such as age-related macular degeneration, atypical hemolytic uremic syndrome, and dense deposit disease. Our data can be used to guide future experimental studies.

Compstatin family peptides are potent inhibitors of the complement system and promising drug candidates against diseases involving under-regulated complement activation. Compstatin is a 13-residue cyclized peptide that inhibits cleavage of complement protein C3, preventing downstream complement activation. We present three new compstatin variants, characterized by tryptophan replacement at positions 1 and/or 13. Peptide design was based on physicochemical reasoning and was inspired by earlier work which identified tryptophan substitutions at positions 1 and 13 in peptides with predicted C3c binding abilities (Bellows, M. L.; Fung, H. K.; Taylor, M. S.; Floudas, C. A.; Lopez de Victoria, A.; Morikis, D. (2010) Biophys J 98: 2337–2346). The new variants preserve distinct polar and nonpolar surfaces of compstatin, but have altered local interaction capabilities with C3. All three peptides exhibited potent C3 binding by surface plasmon resonance (SPR) and potent complement inhibition by ELISAs. We also present ELISA data and detailed SPR kinetic data of three peptides from previous computational design.

Electrostatically-driven association of proteins is important to many biological functions, and understanding which amino acid residues contribute to these interactions is crucial to protein design. Theoretical calculations that are used to elucidate the role of electrostatics in association are typically based on a single experimentally determined protein structure, while an underlying rigid-body assumption is adopted. The goal of this study was to investigate the role of conformational fluctuations on electrostatic interaction energies, as applied to the electrostatic analysis of barnase–barstar. For our calculations, we apply theoretical alanine-scan mutagenesis to introduce charge perturbations by replacing every ionizable residue with alanine, one at a time. Electrostatic clustering and free energy calculations based on the Poisson–Boltzmann method are used to evaluate the effects of each perturbation. Molecular dynamics simulations are performed for the barnase–barstar parent complex and seven experimental alanine mutations from the literature, in order to introduce relaxation before and after mutation. We discuss the effects of dynamics, in the form of pre- and post-mutation relaxation, on electrostatic clustering and free energies of association in light of experimental data. We also examine the utility of nine electrostatic similarity methods for clustering of barnase alanine-scan mutants. Our calculations suggest that the rigid-body assumption is reasonable for electrostatic calculations of barnase–barstar.

The development of compounds to regulate the activation of the complement system in non-primate species is of profound interest because it can provide models for human diseases. The peptide compstatin inhibits protein C3 in primate mammals and is a potential therapeutic agent against unregulated activation of complement in humans but is inactive against nonprimate species. Here, we elucidate this species specificity of compstatin by molecular dynamics simulations of complexes between the most potent natural compstatin analog and human or rat C3. The results are compared against an experimental conformation of the human complex, determined recently by X-ray diffraction at 2.4-A resolution. The human complex simulations provide information on the relative contributions to stability of specific C3 and compstatin residues. In the rat simulations, the protein undergoes reproducible conformational changes, which eliminate or weaken specific interactions and reduce the complex stability. The simulation insights can be used to design improved compstatin-based inhibitors for human C3 and active inhibitors against lower mammals.

Two de novo protein design frameworks are applied to the discovery of new compstatin variants. One is based on sequence selection and fold specificity, whereas the other approach is based on sequence selection and approximate binding affinity calculations. The proposed frameworks were applied to a complex of C3c with compstatin variant E1 and new variants with improved binding affinities are predicted and experimentally validated. The computational studies elucidated key positions in the sequence of compstatin that greatly affect the binding affinity. Positions 4 and 13 were found to favor Trp, whereas positions 1, 9, and 10 are dominated by Asn, and position 11 consists mainly of Gln. A structural analysis of the C3c-bound peptide analogs is presented.

Actin-depolymerizing-factor (ADF)/cofilins have emerged as key regulators of cytoskeletal dynamics in cell motility, morphogenesis, endocytosis, and cytokinesis. The activities of ADF/cofilins are regulated by membrane phospholipid PI(4,5)P(2) in vitro and in cells, but the mechanism of the ADF/cofilin-PI(4,5)P(2) interaction has remained controversial. Recent studies suggested that ADF/cofilins interact with PI(4,5)P(2) through a specific binding pocket, and that this interaction is dependent on pH. Here, we combined systematic mutagenesis with biochemical and spectroscopic methods to elucidate the phosphoinositide-binding mechanism of ADF/cofilins. Our analysis revealed that cofilin does not harbor a specific PI(4,5)P(2)-binding pocket, but instead interacts with PI(4,5)P(2) through a large, positively charged surface of the molecule. Cofilin interacts simultaneously with multiple PI(4,5)P(2) headgroups in a cooperative manner. Consequently, interactions of cofilin with membranes and actin exhibit sharp sensitivity to PI(4,5)P(2) density. Finally, we show that cofilin binding to PI(4,5)P(2) is not sensitive to changes in the pH at physiological salt concentration, although the PI(4,5)P(2)-clustering activity of cofilin is moderately inhibited at elevated pH. Collectively, our data demonstrate that ADF/cofilins bind PI(4,5)P(2) headgroups through a multivalent, cooperative mechanism, and suggest that the actin filament disassembly activity of ADF/cofilin can be accurately regulated by small changes in the PI(4,5)P(2) density at cellular membranes.