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.