Search

The complement system is assembled from a network of proteins that function to bring about the first line of defense of the body against invading pathogens. However, complement deficiencies or invasive pathogens can hijack complement to subsequently increase susceptibility of the body to infections. Moreover, 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 mathematical 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 and (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 <<1.0% of the pathogen surface. Conversely, the MAC production in nasopharynx occupies about 1.5-10% of the total N. meningitidis surface, thus making nasal MAC levels at least about eight orders of magnitude higher. Altogether, we predict complement-imbalance, favoring overactivation, is associated with nasopharynx homeostasis. Conversely, orientating toward complement-balance may cause disruption to the nasopharynx homeostasis. Thus, for sporadic meningococcal disease, our model predicts rising nasal levels of complement regulators as early infection biomarkers., Competing Interests: Conflict of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Narkhede, Gautam, Hsu, Rodriguez, Zewde, Harrison, Arantes, Gaieb, Gorham, Kieslich, Morikis, Sahu and Palermo.)

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.

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., Competing Interests: The authors declare no competing financial interest., (Copyright © 2020 American Chemical Society.)

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.

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 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., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

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., (Copyright © 2018 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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.

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 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., Competing Interests: The authors declare the following competing financial interest(s): D.M., P.T., R.D.G., R.R.M., and R.E.S.H. are coinventors in a patent and/or patent application(s) for compstatin family peptides.

The complement system is an intricate defense network that rapidly removes invading pathogens. Although many complement regulators are present to protect host cells under homeostasis, the impairment of Factor H (FH) regulatory mechanism has been associated with several autoimmune and inflammatory diseases. To understand the dynamics involved in the pivotal balance between activation and regulation, we have developed a comprehensive computational model of the alternative and classical pathways of the complement system. The model is composed of 290 ordinary differential equations with 142 kinetic parameters that describe the state of complement system under homeostasis and disorder through FH impairment. We have evaluated the state of the system by generating concentration-time profiles for the biomarkers C3, C3a-desArg, C5, C5a-desArg, Factor B (FB), Ba, Bb, and fC5b-9 that are influenced by complement dysregulation. We show that FH-mediated disorder induces substantial levels of complement activation compared to homeostasis, by generating reduced levels of C3 and FB, and to a lesser extent C5, and elevated levels of C3a-desArg, Ba, Bb, C5a-desArg, and fC5b-9. These trends are consistent with clinically observed biomarkers associated with complement-mediated diseases. Furthermore, we introduced therapy states by modeling known inhibitors of the complement system, a compstatin variant (C3 inhibitor) and eculizumab (C5 inhibitor). Compstatin demonstrates strong restorative effects for early-stage biomarkers, such as C3a-desArg, FB, Ba, and Bb, and milder restorative effects for late-stage biomarkers, such as C5a-desArg and fC5b-9, whereas eculizumab has strong restorative effects on late-stage biomarkers, and negligible effects on early-stage biomarkers. These results highlight the need for patient-tailored therapies that target early complement activation at the C3 level, or late-stage propagation of the terminal cascade at the C5 level, depending on the specific FH-mediated disease and the manifestations of a patient's genetic profile in complement regulatory function., Competing Interests: The authors have received a research grant from Achillion Pharmaceuticals, a company that specializes on complement-based drug discovery, for a disease-specific computational model. DM has received lecture/travel honorarium from Achillion Pharmaceuticals. DM is a co-inventor in a patent on compstatin peptides (Title: Compstatin Analogs; US Patent 9512180), and in provisional and patent applications on compstatin peptides. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

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., (Copyright © 2017 Elsevier Inc. All rights reserved.)

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., (Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Structure and dynamics are essential elements of protein function. Protein structure is constantly fluctuating and undergoing conformational changes, which are captured by molecular dynamics (MD) simulations. We introduce a computational framework that provides a compact representation of the dynamic conformational space of biomolecular simulations. This method presents a systematic approach designed to reduce the large MD simulation spatiotemporal datasets into a manageable set in order to guide our understanding of how protein mechanics emerge from side chain organization and dynamic reorganization. We focus on the detection of side chain interactions that undergo rearrangements mediating global domain motions and vice versa. Side chain rearrangements are extracted from side chain interactions that undergo well-defined abrupt and persistent changes in distance time series using Gaussian mixture models, whereas global domain motions are detected using dynamic cross-correlation. Both side chain rearrangements and global domain motions represent the dynamic components of the protein MD simulation, and are both mapped into a network where they are connected based on their degree of coupling. This method allows for the study of allosteric communication in proteins by mapping out the protein dynamics into an intramolecular network to reduce the large simulation data into a manageable set of communities composed of coupled side chain rearrangements and global domain motions. This computational framework is suitable for the study of tightly packed proteins, such as G protein-coupled receptors, and we present an application on a seven microseconds MD trajectory of CC chemokine receptor 7 (CCR7) bound to its ligand CCL21.

Purpose: To redesign a complement-inhibiting peptide with the potential to become a therapeutic for dry and wet age-related macular degeneration (AMD)., Methods: We present a new potent peptide (Peptide 2) of the compstatin family. The peptide is developed by rational design, based on a mechanistic binding hypothesis, and structural and physicochemical properties derived from molecular dynamics (MD) simulation. The inhibitory activity, efficacy, and solubility of Peptide 2 are evaluated using a hemolytic assay, a human RPE cell-based assay, and ultraviolet (UV) absorption properties, respectively, and compared to the respective properties of its parent peptide (Peptide 1)., Results: The sequence of Peptide 2 contains an arginine-serine N-terminal extension (a characteristic of parent Peptide 1) and a novel 8-polyethylene glycol (PEG) block C-terminal extension. Peptide 2 has significantly improved aqueous solubility compared to Peptide 1 and comparable complement inhibitory activity. In addition, Peptide 2 is more efficacious in inhibiting complement activation in a cell-based model that mimics the pathobiology of dry AMD., Conclusions: We have designed a new peptide analog of compstatin that combines N-terminal polar amino acid extensions and C-terminal PEGylation extensions. This peptide demonstrates significantly improved aqueous solubility and complement inhibitory efficacy, compared to the parent peptide. The new peptide overcomes the aggregation limitation for clinical translation of previous compstatin analogs and is a candidate to become a therapeutic for the treatment of AMD.

Biased ligand binding to G protein-coupled receptors enables functional selectivity of intracellular effectors to mediate cellular function. Despite the significant advances made in characterizing the conformational states (transmembrane helical arrangements) capable of discriminating between G protein and arrestin binding, the role of the ligand in stabilizing such conformations remains unclear. To address this issue, we simulate microsecond dynamics of CC chemokine receptor 7 (CCR7) bound to its native biased ligands, CCL19 and CCL21, and detect a series of molecular switches that are mediated by various ligand-induced allosteric events. These molecular switches involve three tyrosine residues (Y112(3.32), Y255(6.51), and Y288(7.39)), three phenylalanine residues (F116(3.36), F208(5.47), and F248(6.44)), and a polar interaction between Q252(6.48) and R294(7.45) in the transmembrane domain of CCR7. Conformational changes within these switches, particularly hydrogen bond formation between Y112(3.32) and Y255(6.51), lead to global helical movements in the receptor's transmembrane helices and contribute to the transitioning of the receptor to distinct states. Ligand-induced helical movements in the receptor highlight the ability of biased ligands to stabilize the receptor in different states through a dynamic network of allosteric events.

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 complement system is an integral part of innate immunity that detects and eliminates invading pathogens through a cascade of reactions. The destructive effects of the complement activation on host cells are inhibited through versatile regulators that are present in plasma and bound to membranes. Impairment in the capacity of these regulators to function in the proper manner results in autoimmune diseases. To better understand the delicate balance between complement activation and regulation, we have developed a comprehensive quantitative model of the alternative pathway. Our model incorporates a system of ordinary differential equations that describes the dynamics of the four steps of the alternative pathway under physiological conditions: (i) initiation (fluid phase), (ii) amplification (surfaces), (iii) termination (pathogen), and (iv) regulation (host cell and fluid phase). We have examined complement activation and regulation on different surfaces, using the cellular dimensions of a characteristic bacterium (E. coli) and host cell (human erythrocyte). In addition, we have incorporated neutrophil-secreted properdin into the model highlighting the cross talk of neutrophils with the alternative pathway in coordinating innate immunity. Our study yields a series of time-dependent response data for all alternative pathway proteins, fragments, and complexes. We demonstrate the robustness of alternative pathway on the surface of pathogens in which complement components were able to saturate the entire region in about 54 minutes, while occupying less than one percent on host cells at the same time period. Our model reveals that tight regulation of complement starts in fluid phase in which propagation of the alternative pathway was inhibited through the dismantlement of fluid phase convertases. Our model also depicts the intricate role that properdin released from neutrophils plays in initiating and propagating the alternative pathway during bacterial infection.

The exocyst complex regulates the last steps of exocytosis, which is essential to organisms across kingdoms. In humans, its dysfunction is correlated with several significant diseases, such as diabetes and cancer progression. Investigation of the dynamic regulation of the evolutionarily conserved exocyst-related processes using mutants in genetically tractable organisms such as Arabidopsis thaliana is limited by the lethality or the severity of phenotypes. We discovered that the small molecule Endosidin2 (ES2) binds to the EXO70 (exocyst component of 70 kDa) subunit of the exocyst complex, resulting in inhibition of exocytosis and endosomal recycling in both plant and human cells and enhancement of plant vacuolar trafficking. An EXO70 protein with a C-terminal truncation results in dominant ES2 resistance, uncovering possible distinct regulatory roles for the N terminus of the protein. This study not only provides a valuable tool in studying exocytosis regulation but also offers a potentially new target for drugs aimed at addressing human disease.

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.

Background: To extend the potency of vaccines against infectious diseases, vaccines should be able to exploit multiple arms of the immune system. One component of the immune system that is under-used in vaccine design is the subset of B cells known to be capable of responding to repetitive antigenic epitopes and differentiate into plasma cells even in the absence of T cell help (T-independent, TI)., Results: To target vaccine responses from T-independent B cells, we reengineered a bacterial Flagellin (FliC) by replacing its exposed D3 domain with a viral envelope protein from Dengue virus (DENV2). The resulting hybrid FliC protein (hFliC) was able to form stable filaments decorated with conformationally intact DENV2 envelope domains. These filaments were not only capable of inducing a T cell-dependent (TD) humoral antibody response, but also significant IgM and IgG3 antibody response in a helper T cell repertoire-restricted transgenic mouse model., Conclusions: Our results provide proof-of-principle demonstration that a reengineered hybrid FliC could be used as a platform for polymeric subunit vaccines, enhancing T cell-dependent and possibly inducing T-independent antibody responses from B-1 B cells as well.

As a part of innate immunity, the complement system relies on activation of the alternative pathway (AP). While feed-forward amplification generates an immune response towards foreign surfaces, the process requires regulation to prevent an immune response on the surface of host cells. Factor H (FH) is a complement protein secreted by native cells to negatively regulate the AP. In terms of structure, FH is composed of 20 complement-control protein (CCP) modules that are structurally homologous but vary in composition and function. Mutations in these CCPs have been linked to states of autoimmunity. In particular, several mutations in CCP 19-20 are correlated to atypical hemolytic uremic syndrome (aHUS). From crystallographic structures there are three putative binding sites of CCP 19-20 on C3d. Since there has been some controversy over the primary mode of binding from experimental studies, we approach characterization of binding using computational methods. Specifically, we compare each binding mode in terms of electrostatic character, structural stability, dissociative and associative properties, and predicted free energy of binding. After a detailed investigation, we found two of the three binding sites to be similarly stable while varying in the number of contacts to C3d and in the energetic barrier to complex dissociation. These sites are likely physiologically relevant and may facilitate multivalent binding of FH CCP 19-20 to C3b and either C3d or host glycosaminoglycans. We propose thermodynamically stable binding with modules 19 and 20, the latter driven by electrostatics, acting synergistically to increase the apparent affinity of FH for host surfaces., (© 2015 The Protein Society.)

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., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

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.

Adenoviruses bind to a variety of human cells to cause infection. Both the B2 adenovirus 11 and B1 adenovirus 21 use protein knobs to bind to complement regulator CD46(SCR1-2) in order to gain entry into host cells. In each complex, the two proteins are highly negatively charged but bind to each other at an interface with oppositely charged surface patches. We computationally generated single-alanine mutants of charged residues in the complexes CD46(SCR1-2)-Ad11k and CD46(SCR1-2)-Ad21k. We used electrostatic clustering and Poisson-Boltzmann free energy calculations to propose a hypothesis on the role of electrostatics in association. Our results delineate specific interfacial electrostatic interactions that are critical for association in both CD46(SCR1-2)-Ad11k and CD46(SCR1-2)-Ad21k. These results will serve as a predictive tool in the selection of mutants with desired binding affinity in experimental mutagenesis studies. This study will also serve as a foundation for the design of inhibitors to treat adenovirus infections.

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., (Copyright © 2014 The Authors. Published by Elsevier B.V. All rights reserved.)

The complement system surveillance in the host is effective in controlling viral propagation. Consequently, to subvert this effector mechanism, viruses have developed a series of adaptations. One among these is encoding mimics of host regulators of complement activation (RCA) which help viruses to avoid being labeled as 'foreign' and protect them from complement-mediated neutralization and complement-enhanced antiviral adaptive immunity. In this review, we provide an overview on the structure, function and evolution of viral RCA proteins., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

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.

Staphylococcus aureus expresses numerous virulence factors that aid in immune evasion. The four-domain staphylococcal immunoglobulin binding (Sbi) protein interacts with complement component 3 (C3) and its thioester domain (C3d)-containing fragments. Recent structural data suggested two possible modes of binding of Sbi domain IV (Sbi-IV) to C3d, but the physiological binding mode remains unclear. We used a computational approach to provide insight into the C3d-Sbi-IV interaction. Molecular dynamics (MD) simulations showed that the first binding mode (PDB code 2WY8) is more robust than the second (PDB code 2WY7), with more persistent polar and nonpolar interactions, as well as conserved interfacial solvent accessible surface area. Brownian dynamics and steered MD simulations revealed that the first binding mode has faster association kinetics and maintains more stable intermolecular interactions compared to the second binding mode. In light of available experimental and structural data, our data confirm that the first binding mode represents Sbi-IV interaction with C3d (and C3) in a physiological context. Although the second binding mode is inherently less stable, we suggest a possible physiological role. Both binding sites may serve as a template for structure-based design of novel complement therapeutics., (Copyright © 2014 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

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.

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., (© 2012 John Wiley & Sons A/S.)

Background: The V3 loop of the glycoprotein gp120 of HIV-1 plays an important role in viral entry into cells by utilizing as coreceptor CCR5 or CXCR4, and is implicated in the phenotypic tropisms of HIV viruses. It has been hypothesized that the interaction between the V3 loop and CCR5 or CXCR4 is mediated by electrostatics. We have performed hierarchical clustering analysis of the spatial distributions of electrostatic potentials and charges of V3 loop structures containing consensus sequences of HIV-1 subtypes., Results: Although the majority of consensus sequences have a net charge of +3, the spatial distribution of their electrostatic potentials and charges may be a discriminating factor for binding and infectivity. This is demonstrated by the formation of several small subclusters, within major clusters, which indicates common origin but distinct spatial details of electrostatic properties. Some of this information may be present, in a coarse manner, in clustering of sequences, but the spatial details are largely lost. We show the effect of ionic strength on clustering of electrostatic potentials, information that is not present in clustering of charges or sequences. We also make correlations between clustering of electrostatic potentials and net charge, coreceptor selectivity, global prevalence, and geographic distribution. Finally, we interpret coreceptor selectivity based on the N6X7T8|S8X9 sequence glycosylation motif, the specific positive charge location according to the 11/24/25 rule, and the overall charge and electrostatic potential distribution., Conclusions: We propose that in addition to the sequence and the net charge of the V3 loop of each subtype, the spatial distributions of electrostatic potentials and charges may also be important factors for receptor recognition and binding and subsequent viral entry into cells. This implies that the overall electrostatic potential is responsible for long-range recognition of the V3 loop with coreceptors CCR5/CXCR4, whereas the charge distribution contributes to the specific short-range interactions responsible for the formation of the bound complex. We also propose a scheme for coreceptor selectivity based on the sequence glycosylation motif, the 11/24/25 rule, and net charge.

The interaction between complement fragment C3d and complement receptor 2 (CR2) is a key aspect of complement immune system activation, and is a component in a link between innate and adaptive immunities. The complement immune system is an ancient mechanism for defense, and can be found in species that have been on Earth for the last 600 million years. However, the link between the complement system and adaptive immunity, which is formed through the association of the B-cell co-receptor complex, including the C3d-CR2 interaction, is a much more recent adaptation. Human C3d and CR2 have net charges of -1 and +7 respectively, and are believed to have evolved favoring the role of electrostatics in their functions. To investigate the role of electrostatics in the function and evolution of human C3d and CR2, we have applied electrostatic similarity methods to identify regions of evolutionarily conserved electrostatic potential based on 24 homologues of complement C3d and 4 homologues of CR2. We also examine the effects of structural perturbation, as introduced through molecular dynamics and mutations, on spatial distributions of electrostatic potential to identify perturbation resistant regions, generated by so-called electrostatic "hot-spots". Distributions of electrostatic similarity based on families of perturbed structures illustrate the presence of electrostatic "hot-spots" at the two functional sites of C3d, while the surface of CR2 lacks electrostatic "hot-spots" despite its excessively positive nature. We propose that the electrostatic "hot-spots" of C3d have evolved to optimize its dual-functionality (covalently attaching to pathogen surfaces and interaction with CR2), which are both necessary for the formation B-cell co-receptor complexes. Comparison of the perturbation resistance of the electrostatic character of the homologues of C3d suggests that there was an emergence of a new role of electrostatics, and a transition in the function of C3d, after the divergence of jawless fish.

The V3 loop of the glycoprotein 120 (gp120) is a contact point for cell entry of HIV-1 leading to infection. Despite sequence variability and lack of specific structure, the highly flexible V3 loop possesses a well-defined role in recognizing and selecting cell-bound coreceptors CCR5 and CXCR4 through a mechanism of charge complementarity. We have performed two independent molecular dynamics (MD) simulations to gain insights into the dynamic character of two V3 loops with slightly different sequences, but significantly different starting crystallographic structures. We have identified highly populated trajectory-specific salt bridges between oppositely charged stem residues Arg9 and Glu25 or Asp29. The two trajectories share nearly identical correlated motions within the simulations, despite their different overall structures. High occupancy salt bridges play a key role in the major cross-correlated motions in both trajectories, and may be responsible for transient structural stability in preparation for coreceptor binding. In addition, the two V3 loops visit conformations with similarities in spatial distributions of electrostatic potentials, despite their inherent flexibility, which may play a role in coreceptor recognition. It is plausible that cooperativity between overall electrostatic potential, charged residue interactions, and correlated motions could be associated with a coreceptor selection and binding.

For systems involving highly and oppositely charged proteins, electrostatic forces dominate association and contribute to biomolecular complex stability. Using experimental or theoretical alanine-scanning mutagenesis, it is possible to elucidate the contribution of individual ionizable amino acids to protein association. We evaluated our electrostatic free energy calculations by comparing calculated and experimental data for alanine mutants of five protein complexes. We calculated Poisson-Boltzmann electrostatic free energies based on a thermodynamic cycle, which incorporates association in a reference (Coulombic) and solvated (solution) state, as well as solvation effects. We observe that Coulombic and solvation free energy values correlate with experimental data in highly and oppositely charged systems, but not in systems comprised of similarly charged proteins. We also observe that correlation between solution and experimental free energies is dependent on dielectric coefficient selection for the protein interior. Free energy correlations improve as protein dielectric coefficient increases, suggesting that the protein interior experiences moderate dielectric screening, despite being shielded from solvent. We propose that higher dielectric coefficients may be necessary to more accurately predict protein-protein association. Additionally, our data suggest that Coulombic potential calculations alone may be sufficient to predict relative binding of protein mutants., (2011 Wiley Periodicals, Inc.)