Search for: All records

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. 

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. 

Abstract Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., “Gaussian approximation”). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica‐exchange GaMD (rex‐GaMD) and replica‐exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new “selective GaMD” algorithms including the Ligand GaMD (LiGaMD) and Peptide GaMD (Pep‐GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small‐molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein–protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design. This article is categorized under:Structure and Mechanism > Computational Biochemistry and BiophysicsMolecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo MethodsMolecular and Statistical Mechanics > Free Energy Methods