Role of Electrostatic Hotspots in the Selectivity of Complement Control Proteins Toward Human and Bovine Complement Inhibition.

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.

Systems Biology Modeling of the Complement System Under Immune Susceptible Pathogens.

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.

Factor H-Inspired Design of Peptide Biomarkers of the Complement C3d Protein.

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.

Spontaneous Embedding of DNA Mismatches Within the RNA:DNA Hybrid of CRISPR-Cas9.

CRISPR-Cas9 is the forefront technology for editing the genome. In this system, the Cas9 protein is programmed with guide RNAs to process DNA sequences that match the guide RNA forming an RNA:DNA hybrid structure. However, the binding of DNA sequences that do not fully match the guide RNA can limit the applicability of CRISPR-Cas9 for genome editing, resulting in the so-called off-target effects. Here, molecular dynamics is used to probe the effect of DNA base pair mismatches within the RNA:DNA hybrid in CRISPR-Cas9. Molecular simulations revealed that the presence of mismatched pairs in the DNA at distal sites with respect to the Protospacer Adjacent Motif (PAM) recognition sequence induces an extended opening of the RNA:DNA hybrid, leading to novel interactions established by the unwound nucleic acids and the protein counterpart. On the contrary, mismatched pairs upstream of the RNA:DNA hybrid are rapidly incorporated within the heteroduplex, with minor effect on the protein-nucleic acid interactions. As a result, mismatched pairs at PAM distal ends interfere with the activation of the catalytic HNH domain, while mismatches fully embedded in the RNA:DNA do not affect the HNH dynamics and enable its activation to cleave the DNA. These findings provide a mechanistic understanding to the intriguing experimental evidence that PAM distal mismatches hamper a proper function of HNH, explaining also why mismatches within the heteroduplex are much more tolerated. This constitutes a step forward in understanding off-target effects in CRISPR-Cas9, which encourages novel structure-based engineering efforts aimed at preventing the onset of off-target effects.

A computational model for the evaluation of complement system regulation under homeostasis, disease, and drug intervention.

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.

Quantitative Modeling of the Alternative Pathway of the Complement System.

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.