Simulated mutagenesis of the hypervariable loops of a llama VHH domain for the recovery of canonical conformations

General Trends of the Camelidae Antibody VHHs Domain Dynamics

Conformational flexibility plays an essential role in antibodies' functional and structural stability. They facilitate and determine the strength of antigen-antibody interactions. Camelidae express an interesting subtype of single-chain antibody, named Heavy Chain only Antibody. They have only one N-terminal Variable domain (V_HH) per chain, composed of Frameworks (FRs) and Complementarity Determining regions (CDRs) like their VH and VL counterparts in IgG. Even when expressed independently, V_HH domains display excellent solubility and (thermo)stability, which helps them to retain their impressive interaction capabilities. Sequence and structural features of V_HH domains contributing to these abilities have already been studied compared to classical antibodies. To have the broadest view and understand the changes in dynamics of these macromolecules, large-scale molecular dynamics simulations for a large number of non-redundant V_HH structures have been performed for the first time. This analysis reveals the most prevalent movements in these domains. It reveals the four main classes of V_HHs dynamics. Diverse local changes were observed in CDRs with various intensities. Similarly, different types of constraints were observed in CDRs, while FRs close to CDRs were sometimes primarily impacted. This study sheds light on the changes in flexibility in different regions of V_HH that may impact their in silico design.

Computational affinity maturation of camelid single-domain intrabodies against the nonamyloid component of alpha-synuclein

Improving the affinity of protein-protein interactions is a challenging problem that is particularly important in the development of antibodies for diagnostic and clinical use. Here, we used structure-based computational methods to optimize the binding affinity of V_HNAC1, a single-domain intracellular antibody (intrabody) from the camelid family that was selected for its specific binding to the nonamyloid component (NAC) of human α-synuclein (α-syn), a natively disordered protein, implicated in the pathogenesis of Parkinson's disease (PD) and related neurological disorders. Specifically, we performed ab initio modeling that revealed several possible modes of V_HNAC1 binding to the NAC region of α-syn as well as mutations that potentially enhance the affinity between these interacting proteins. While our initial design strategy did not lead to improved affinity, it ultimately guided us towards a model that aligned more closely with experimental observations, revealing a key residue on the paratope and the participation of H4 loop residues in binding, as well as confirming the importance of electrostatic interactions. The binding activity of the best intrabody mutant, which involved just a single amino acid mutation compared to parental V_HNAC1, was significantly enhanced primarily through a large increase in association rate. Our results indicate that structure-based computational design can be used to successfully improve the affinity of antibodies against natively disordered and weakly immunogenic antigens such as α-syn, even in cases such as ours where crystal structures are unavailable.

Tilting the balance between canonical and noncanonical conformations for the H1 hypervariable loop of a llama VHH through point mutations

Nanobodies are single-domain antibodies found in camelids. These are the smallest naturally occurring binding domains and derive functionality via three hypervariable loops (H1-H3) that form the binding surface. They are excellent candidates for antibody engineering because of their favorable characteristics like small size, high solubility, and stability. To rationally engineer antibodies with affinity for a specific target, the hypervariable loops can be tailored to obtain the desired binding surface. As a first step toward such a goal, we consider the design of loops with a desired conformation. In this study, we focus on the H1 loop of the anti-hCG llama nanobody that exhibits a noncanonical conformation. We aim to "tilt" the stability of the H1 loop structure from a noncanonical conformation to a (humanized) type 1 canonical conformation by studying the effect of selected mutations to the amino acid sequence of the H1, H2, and proximal residues. We use all-atomistic, explicit-solvent, biased molecular dynamic simulations to simulate the wild-type and mutant loops in a prefolded framework. We thus find mutants with increasing propensity to form a stable type 1 canonical conformation of the H1 loop. Free energy landscapes reveal the existence of conformational isomers of the canonical conformation that may play a role in binding different antigenic surfaces. We also elucidate the approximate mechanism and kinetics of transitions between such conformational isomers by using a Markovian model. We find that a particular three-point mutant has the strongest thermodynamic propensity to form the H1 type 1 canonical structure but also to exhibit transitions between conformational isomers, while a different, more rigid three-point mutant has the strongest propensity to be kinetically trapped in such a canonical structure.

Characterizing the structural behavior of selected Aβ-42 monomers with different solubilities

The conformational behavior of the wild-type amyloid β-42 (Aβ-42) monomer and two of its mutants was explored via all-atom replica exchange molecular dynamics simulations in explicit solvent, to identify structural features that may promote or deter early-stage oligomerization. The markers used for this purpose indicate that while the three peptides are relatively flexible they have distinct preferential structures and degree of rigidity. In particular, we found that one mutant that remains in the monomeric state in experiments displays a characteristic N-terminal structure that significantly enhances its rigidity. This finding is consistent with various studies that have detected a reduction in oligomerization frequency and Aβ-related toxicity upon sequence-specific antibody or ligand binding to the N-terminal tail of wild-type monomers, likely leading to the stabilization of this region. In general, our results highlight a potential role of the N-terminal segment on Aβ oligomerization and give insights into specific interactions that may be responsible for promoting the pronounced structural changes observed upon introducing point mutations on the wild-type Aβ-42 peptide.

Theoretical aspects of immunity

The immune system recognizes a myriad of invading pathogens and their toxic products. It does so with a finite repertoire of antibodies and T cell receptors. We here describe theories that quantify the dynamics of the immune system. We describe how the immune system recognizes antigens by searching the large space of receptor molecules. We consider in some detail the theories that quantify the immune response to influenza and dengue fever. We review theoretical descriptions of the complementary evolution of pathogens that occurs in response to immune system pressure. Methods including bioinformatics, molecular simulation, random energy models, and quantum field theory contribute to a theoretical understanding of aspects of immunity.