Molecular description of flexibility in an antibody combining site

Interplay of the mechanical and structural properties of DNA nanostructures determines their electrostatic interactions with lipid membranes

Nucleic acids and lipids function in close proximity in biological processes, as well as in nanoengineered constructs for therapeutic applications. As both molecules carry a rich charge profile, and frequently coexist in complex ionic solutions, the electrostatics surely play a pivotal role in interactions between them. Here we discuss how each component of a DNA/ion/lipid system determines its electrostatic attachment. We examine membrane binding of a library of DNA molecules varying from nanoengineered DNA origami through plasmids to short DNA domains, demonstrating the interplay between the molecular structure of the nucleic acid and the phase of lipid bilayers. Furthermore, the magnitude of DNA/lipid interactions is tuned by varying the concentration of magnesium ions in the physiologically relevant range. Notably, we observe that the structural and mechanical properties of DNA are critical in determining its attachment to lipid bilayers and demonstrate that binding is correlated positively with the size, and negatively with the flexibility of the nucleic acid. The findings are utilized in a proof-of-concept comparison of membrane interactions of two DNA origami designs - potential nanotherapeutic platforms - showing how the results can have a direct impact on the choice of DNA geometry for biotechnological applications.

Significance of Receptor Mobility in Multivalent Binding on Lipid Membranes

Numerous key biological processes rely on the concept of multivalency, where ligands achieve stable binding only upon engaging multiple receptors. These processes, like viral entry or immune synapse formation, occur on the diffusive cellular membrane. One crucial, yet underexplored aspect of multivalent binding is the mobility of coupled receptors. Here, we discuss the consequences of mobility in multivalent processes from four perspectives: (I) The facilitation of receptor recruitment by the multivalent ligand due to their diffusivity prior to binding. (II) The effects of receptor preassembly, which allows their local accumulation. (III) The consequences of changes in mobility upon the formation of receptor/ligand complex. (IV) The changes in the diffusivity of lipid environment surrounding engaged receptors. We demonstrate how understanding mobility is essential for fully unravelling the principles of multivalent membrane processes, leading to further development in studies on receptor interactions, and guide the design of new generations of multivalent ligands.

Molecular dynamics analysis of the binding of human interleukin-6 with interleukin-6 α-receptor

Human interleukin-6 (hIL-6) is a multifunctional cytokine that regulates immune and inflammatory responses in addition to metabolic and regenerative processes and cancer. hIL-6 binding to the IL-6 receptor (IL-6Rα) induces homodimerization and recruitment of the glycoprotein (gp130) to form a hexameric signaling complex. Anti-IL-6 and IL-6R antibodies are clinically approved inhibitors of IL-6 signaling pathway for treating rheumatoid arthritis and Castleman's disease, respectively. There is a potential to develop novel small molecule IL-6 antagonists derived from understanding the structural basis for IL-6/IL-6Rα interactions. Here, we combine homology modeling with extensive molecular dynamics (MD) simulations to examine the association of hIL-6 with IL-6Rα. A comparison with MD of apo hIL-6 reveals that the binding of hIL-6 to IL-6Rα induces structural and dynamic rearrangements in the AB loop region of hIL-6, disrupting intraprotein contacts and increasing the flexibility of residues 48 to 58 of the AB loop. In contrast, due to the involvement of residues 59 to 78 in forming contacts with the receptor, these residues of the AB loop are observed to rigidify in the presence of the receptor. The binary complex is primarily stabilized by two pairs of salt bridges, Arg181 (hIL-6)- Glu182 (IL-6Rα) and Arg184 (hIL-6)- Glu183 (IL-6Rα) as well as hydrophobic and aromatic stacking interactions mediated essentially by Phe residues in both proteins. An interplay of electrostatic, hydrophobic, hydrogen bonding, and aromatic stacking interactions facilitates the formation of the hIL-6/IL-6Rα complex.

Structure and Dynamics of Stacking Interactions in an Antibody Binding Site

For years, antibodies (Abs) have been used as a paradigm for understanding how protein structure contributes to molecular recognition. However, with the ability to evolve Abs that recognize specific chromophores, they also have great potential as models for how protein dynamics contribute to molecular recognition. We previously raised murine Abs to different chromophores and, with the use of three-pulse photon echo peak shift spectroscopy, demonstrated that the immune system is capable of producing Abs with widely varying flexibility. We now report the characterization of the complexes formed between two Abs, 5D11 and 10A6, and the chromophoric ligand that they were evolved to recognize, 8-methoxypyrene-1,3,6-trisulfonic acid (MPTS). The sequences of the Ab genes indicate that they evolved from a common precursor. We also used a variety of spectroscopic methods to probe the photophysics and dynamics of the Ab-MPTS complexes and found that they are similar to each other but distinct from previously characterized anti-MPTS Abs. Structural studies revealed that this difference likely results from a unique mode of binding in which MPTS is sandwiched between the side chain of Phe^H98, which interacts with the chromophore via T-stacking, and the side chain of Trp^L91, which interacts with the chromophore via parallel stacking. The T-stacking interaction appears to mediate relaxation on the picosecond time scale, while the parallel stacking appears to mediate relaxation on an ultrafast, femtosecond time scale, which dominates the response. The anti-MPTS Abs thus not only demonstrate the simultaneous use of the two limiting modes of stacking for molecular recognition, but also provide a unique opportunity to characterize how dynamics might contribute to molecular recognition. Both types of stacking are common in proteins and protein complexes where they may similarly contribute to dynamics and molecular recognition.

Dynamic Aspects of the Immunoglobulin Structure

Immunoglobulin (Ig) molecules are composed of Fab and Fc portions tethered by a hinge region that enables them to rotate and flex, relative to each other. Variable (V) and constant (C) domains of the Fab are connected by a flexible elbow region that is responsible for the movements of the V and C heterodimers. Significant movements of Fc domains have also been documented. The Ig portion's rotational freedom greatly enhances its ability to react with antigens and cell receptors, often simultaneously. The antigen-combining site also displays a dynamic structure. The ability of its various parts to change position greatly facilitates their complexation with various antigenic compounds.

Effective binding to protein antigens by antibodies from antibody libraries designed with enhanced protein recognition propensities

Antibodies provide immune protection by recognizing antigens of diverse chemical properties, but elucidating the amino acid sequence-function relationships underlying the specificity and affinity of antibody-antigen interactions remains challenging. We designed and constructed phage-displayed synthetic antibody libraries with enriched protein antigen-recognition propensities calculated with machine learning predictors, which indicated that the designed single-chain variable fragment variants were encoded with enhanced distributions of complementarity-determining region (CDR) hot spot residues with high protein antigen recognition propensities in comparison with those in the human antibody germline sequences. Antibodies derived directly from the synthetic antibody libraries, without affinity maturation cycles comparable to those in in vivo immune systems, bound to the corresponding protein antigen through diverse conformational or linear epitopes with specificity and affinity comparable to those of the affinity-matured antibodies from in vivo immune systems. The results indicated that more densely populated CDR hot spot residues were sustainable by the antibody structural frameworks and could be accompanied by enhanced functionalities in recognizing protein antigens. Our study results suggest that synthetic antibody libraries, which are not limited by the sequences found in antibodies in nature, could be designed with the guidance of the computational machine learning algorithms that are programmed to predict interaction propensities to molecules of diverse chemical properties, leading to antibodies with optimal characteristics pertinent to their medical applications.

The Potential Role of Solvation in Antibody Recognition of the Lewis Y Antigen

Solvents play an important role in protein folding, protein-protein associations, stability, and specificity of recognition as in the case of antibody-antigen interactions through hydrogen bonds. One of the underappreciated features of protein-associated waters is that it weakens inter- and intra-molecular interactions by modulating electrostatic interactions and influencing conformational changes. Such observations demonstrate the direct relationship between macroscopic solvent effects on protein-protein interactions and atom-scale solvent-protein interactions. Although crystallographic solvents do explain some aspects of solvent-mediated interactions, molecular simulation allows the study of the dynamic role of solvents. Thus, analysis of conformations from molecular simulations are employed to understand the role of solvent on the inherent polyspecificity of a Lewis Y reactive germline gene relative to its expanded hybridomas and a humanized anti-Lewis Y antibody. Our analysis reveals that solvent mediates critical contacts through charged residues to facilitate cross-reactivity to carbohydrate antigens, but also increases the flexibility of some anti-Lewis Y antibodies concomitant with mutations (amino acid substitutions) to the germline antibody. Such flexibility might better allow for recognition and binding of internal structures of extended carbohydrate structures on tumor cells.

Defining the complementarities between antibodies and haptens to refine our understanding and aid the prediction of a successful binding interaction

BACKGROUND

Low molecular weight haptens (<1000 Da) cannot be recognized by the immune system unless conjugated to larger carrier molecules. Antibodies to these exceptionally small antigens can still be generated with exquisite sensitivity. A detailed understanding at the molecular level of this incredible ability of antibodies to recognize haptens, is still limited compared to other antigen classes.

METHODS

Different hapten targets with a broad range of structural flexibility and polarity were conjugated to carrier proteins, and utilized in sheep immunization. Three antibody libraries were constructed and used as potential pools to isolate specific antibodies to each target. The isolated antibodies were analysed in term of CDR length, canonical structure, and binding site shape and electrostatic potential.

RESULTS

The simple, chemically naïve structure of squalane (SQA) was recognized with micromolar sensitivity. An increase in structural rigidity of the hydrophobic and cyclic coprostane (COP) did not improve this binding sensitivity beyond the micromolar range, whilst the polar etioporphyrin (POR) was detected with nanomolar sensitivity. Homoserine lactone (HSL) molecules, which combine molecular flexibility and polarity, generated super-sensitive (picomolar) interactions. To better understand this range of antibody-hapten interactions, analyses were extended to examine the binding loop canonical structures and CDR lengths of a series of anti-hapten clones. Analyses of the pre and post- selection (panning of the phage displayed libraries) sequences revealed more conserved sites (123) within the post-selection sequences, when compared to their pre-selection counterparts (28). The strong selection pressure, generated by panning against these haptens resulted in the isolation of antibodies with significant sequence conservation in the FW regions, and suitable binding site cavities, representing only a relatively small subset of the available full repertoire sequence and structural diversity. As part of this process, the important influence of CDR H2 on antigen binding was observed through its direct interaction with individual antigens and indirect impact on the orientation and the pocket shape, when combined with CDRs H3 and L3. The binding pockets also displayed electrostatic surfaces that were complementary to the hydrophobic nature of COP, SQA, and POR, and the negatively charged HSL.

CONCLUSIONS

The best binding antibodies have shown improved capacity to recognize these haptens by establishing complementary binding pockets in terms of size, shape, and electrostatic potential.

Rigidity Emerges during Antibody Evolution in Three Distinct Antibody Systems: Evidence from QSFR Analysis of Fab Fragments

The effects of somatic mutations that transform polyspecific germline (GL) antibodies to affinity mature (AM) antibodies with monospecificity are compared among three GL-AM Fab pairs. In particular, changes in conformational flexibility are assessed using a Distance Constraint Model (DCM). We have previously established that the DCM can be robustly applied across a series of antibody fragments (VL to Fab), and subsequently, the DCM was combined with molecular dynamics (MD) simulations to similarly characterize five thermostabilizing scFv mutants. The DCM is an ensemble based statistical mechanical approach that accounts for enthalpy/entropy compensation due to network rigidity, which has been quite successful in elucidating conformational flexibility and Quantitative Stability/Flexibility Relationships (QSFR) in proteins. Applied to three disparate antibody systems changes in QSFR quantities indicate that the VH domain is typically rigidified, whereas the VL domain and CDR L2 loop become more flexible during affinity maturation. The increase in CDR H3 loop rigidity is consistent with other studies in the literature. The redistribution of conformational flexibility is largely controlled by nonspecific changes in the H-bond network, although certain Arg to Asp salt bridges create highly localized rigidity increases. Taken together, these results reveal an intricate flexibility/rigidity response that accompanies affinity maturation.

Antibodies are protective proteins used by the immune system to recognize and neutralize foreign objects through interactions with a specific part of the target, called an antigen. Antibody structures are Y-shaped, contain multiple protein chains, and include two antigen-binding sites. The binding sites are located at the end of the Fab fragments, which are the upward facing arms of the Y-structure. The Fab fragments maintain binding affinity by themselves, and are thus often used as surrogates to student antibody-antigen interactions. High affinity antibodies are produced during the course of an immune response by successive mutations to germline gene-encoded antibodies. Germline antibodies are more likely to be polyspecific, whereas the affinity maturation process yields monoclonal antibodies that bind specifically to the target antigen. In this work, we use a computational Distance Constraint Model to characterize how mechanical properties change as three disparate germline antibodies are converted to affinity mature. Our results reveal a rich set of mechanical responses throughout the Fab structure. Nevertheless, increased rigidity in the VH domain is consistently observed, which is consistent with the transition from polyspecificity to monospecificity. That is, flexible antibody structures are able to recognize multiple antigens, while increased affinity and specificity is achieved—in part—by structural rigidification.

Key mutations stabilize antigen-binding conformation during affinity maturation of a broadly neutralizing influenza antibody lineage

Affinity maturation, the process in which somatic hypermutation and positive selection generate antibodies with increasing affinity for an antigen, is pivotal in acquired humoral immunity. We have studied the mechanism of affinity gain in a human B-cell lineage in which two main maturation pathways, diverging from a common ancestor, lead to three mature antibodies that neutralize a broad range of H1 influenza viruses. Previous work showed that increased affinity in the mature antibodies derives primarily from stabilization of the CDR H3 loop in the antigen-binding conformation. We have now used molecular dynamics simulations and existing crystal structures to identify potentially key maturation mutations, and we have characterized their effects on the CDR H3 loop and on antigen binding using further simulations and experimental affinity measurements, respectively. In the two maturation pathways, different contacts between light and heavy chains stabilize the CDR H3 loop. As few as two single-site mutations in each pathway can confer substantial loop stability, but none of them confers experimentally detectable stability on its own. Our results support models of the germinal center reaction in which two or more mutations can occur without concomitant selection and show how divergent pathways have yielded functionally equivalent antibodies. Proteins 2014; 83:771–780. © 2014 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Antigen exposure leads to rigidification of germline antibody combining site

Immune complexes involving diverse antigens and corresponding antibodies were analyzed for mapping conformational transitions of an antibody before antigen binding, upon antigen binding and after antigen release. Molecular dynamics simulations of the two comprehensive datasets consisting of the antigen-free and antigen-bound structures of the germline antibodies 36-65 and BBE6.12H3 provided mechanistic model of antigen encounter by primary antibodies. While native germline antibodies exhibit substantial mobility in the antigen-combining sites, their antigen-bound states exhibit relatively rigid conformations, even in the absence of the antigen suggesting preservation of the structural state after antigen release. It is proposed that acquired rigidity by a germline antibody upon antigen binding may be the first step in affinity maturation in favor of that antigen.

Autoreactivity and exceptional CDR plasticity (but not unusual polyspecificity) hinder elicitation of the anti-HIV antibody 4E10

The broadly-neutralizing anti-HIV antibody 4E10 recognizes an epitope in the membrane-proximal external region of the HIV envelope protein gp41. Previous attempts to elicit 4E10 by vaccination with envelope-derived or reverse-engineered immunogens have failed. It was presumed that the ontogeny of 4E10-equivalent responses was blocked by inherent autoreactivity and exceptional polyreactivity. We generated 4E10 heavy-chain knock-in mice, which displayed significant B cell dysregulation, consistent with recognition of autoantigen/s by 4E10 and the presumption that tolerance mechanisms may hinder the elicitation of 4E10 or 4E10-equivalent responses. Previously proposed candidate 4E10 autoantigens include the mitochondrial lipid cardiolipin and a nuclear splicing factor, 3B3. However, using carefully-controlled assays, 4E10 bound only weakly to cardiolipin-containing liposomes, but also bound negatively-charged, non-cardiolipin-containing liposomes comparably poorly. 4E10/liposome binding was predominantly mediated by electrostatic interactions rather than presumed hydrophobic interactions. The crystal structure of 4E10 free of bound ligands showed a dramatic restructuring of the combining site, occluding the HIV epitope binding site and revealing profound flexibility, but creating an electropositive pocket consistent with non-specific binding of phospholipid headgroups. These results strongly suggested that antigens other than cardiolipin mediate 4E10 autoreactivity. Using a synthetic peptide library spanning the human proteome, we determined that 4E10 displays limited and focused, but unexceptional, polyspecificity. We also identified a novel autoepitope shared by three ER-resident inositol trisphosphate receptors, validated through binding studies and immunohistochemistry. Tissue staining with 4E10 demonstrated reactivity consistent with the type 1 inositol trisphosphate receptor as the most likely candidate autoantigen, but is inconsistent with splicing factor 3B3. These results demonstrate that 4E10 recognition of liposomes competes with MPER recognition and that HIV antigen and autoepitope recognition may be distinct enough to permit eliciting 4E10-like antibodies, evading autoimmunity through directed engineering. However, 4E10 combining site flexibility, exceptional for a highly-matured antibody, may preclude eliciting 4E10 by conventional immunization strategies.

Human germline antibody gene segments encode polyspecific antibodies

Structural flexibility in germline gene-encoded antibodies allows promiscuous binding to diverse antigens. The binding affinity and specificity for a particular epitope typically increase as antibody genes acquire somatic mutations in antigen-stimulated B cells. In this work, we investigated whether germline gene-encoded antibodies are optimal for polyspecificity by determining the basis for recognition of diverse antigens by antibodies encoded by three VH gene segments. Panels of somatically mutated antibodies encoded by a common VH gene, but each binding to a different antigen, were computationally redesigned to predict antibodies that could engage multiple antigens at once. The Rosetta multi-state design process predicted antibody sequences for the entire heavy chain variable region, including framework, CDR1, and CDR2 mutations. The predicted sequences matched the germline gene sequences to a remarkable degree, revealing by computational design the residues that are predicted to enable polyspecificity, i.e., binding of many unrelated antigens with a common sequence. The process thereby reverses antibody maturation in silico. In contrast, when designing antibodies to bind a single antigen, a sequence similar to that of the mature antibody sequence was returned, mimicking natural antibody maturation in silico. We demonstrated that the Rosetta computational design algorithm captures important aspects of antibody/antigen recognition. While the hypervariable region CDR3 often mediates much of the specificity of mature antibodies, we identified key positions in the VH gene encoding CDR1, CDR2, and the immunoglobulin framework that are critical contributors for polyspecificity in germline antibodies. Computational design of antibodies capable of binding multiple antigens may allow the rational design of antibodies that retain polyspecificity for diverse epitope binding.

Computer-aided antibody design

Recent clinical trials using antibodies with low toxicity and high efficiency have raised expectations for the development of next-generation protein therapeutics. However, the process of obtaining therapeutic antibodies remains time consuming and empirical. This review summarizes recent progresses in the field of computer-aided antibody development mainly focusing on antibody modeling, which is divided essentially into two parts: (i) modeling the antigen-binding site, also called the complementarity determining regions (CDRs), and (ii) predicting the relative orientations of the variable heavy (V(H)) and light (V(L)) chains. Among the six CDR loops, the greatest challenge is predicting the conformation of CDR-H3, which is the most important in antigen recognition. Further computational methods could be used in drug development based on crystal structures or homology models, including antibody-antigen dockings and energy calculations with approximate potential functions. These methods should guide experimental studies to improve the affinities and physicochemical properties of antibodies. Finally, several successful examples of in silico structure-based antibody designs are reviewed. We also briefly review structure-based antigen or immunogen design, with application to rational vaccine development.

Cytokinergic IgE Action in Mast Cell Activation

Some 10 years ago it emerged that at sufficiently high concentrations certain monoclonal mouse IgEs exert previously unsuspected effects on mast cells. Thus they can both promote survival and induce activation of mast cells without the requirement for antigens. This was a wake up call that appears to have been missed (or dismissed) by the majority of immunologists. The structural attributes responsible for the potency of the so-called “highly cytokinergic” or HC IgEs have not yet been determined, but the events that ensue when such IgEs bind to the high-affinity receptor, FcεRI, on mast cells have been thoroughly studied, and are strikingly similar to those engendered by antigens when they form cross-linked complexes with the receptors. We review the evidence for the cytokinergic activity of IgE, and the structural features and known properties of immunoglobulins, and of IgE in particular, most likely to be implicated in the phenomenon. We suggest that IgEs with cytokinergic activity may be generated by local germinal center reactions in the target organs of allergy. We consider also the important implications that the existence of cytokinergic IgE may have for a fuller understanding of adaptive immunity and of the action of IgE in asthma and other diseases.

The Fundamental Role of Flexibility on the Strength of Molecular Binding

Non-covalent molecular association underlies a diverse set of biologically and technologically relevant phenomena, including the action of drugs on their biomolecular targets and self- and supra-molecular assembly processes. Computer models employed to model binding frequently use interaction potentials with atomistic detail while neglecting the thermal molecular motions of the binding species. However, errors introduced by this simplification and, more broadly, the thermodynamic consequences of molecular flexibility on binding, are little understood. Here, we isolate the fundamental relationship of molecular flexibility to binding thermodynamics via simulations of simplified molecules with a wide range of flexibilities but the same interaction potential. Disregarding molecular motion is found to generate large errors in binding entropy, enthalpy and free energy, even for molecules that are nearly rigid. Indeed, small decreases in rigidity markedly reduce affinity for highly rigid molecules. Remarkably, precisely the opposite occurs for more flexible molecules, for which increasing flexibility leads to stronger binding affinity. We also find that differences in flexibility suffice to generate binding specificity: for example, a planar surface selectively binds rigid over flexible molecules. Intriguingly, varying molecular flexibility while keeping interaction potentials constant leads to near-linear enthalpy-entropy compensation over a wide range of flexibilities, with the unexpected twist that increasing flexibility produces opposite changes in entropy and enthalpy for molecules in the flexible versus the rigid regime. Molecular flexibility is thus a crucial determinant of binding affinity and specificity and variations in flexibility can lead to strong yet non-intuitive consequences.

Filtering artifacts from lifetime distributions when maximizing entropy using a bootstrapped model

The maximum entropy method (MEM) has been used in many studies to reliably recover effective lifetimes from kinetics, whether measured experimentally or simulated computationally. Here, recent claims made by Mulligan et al. regarding MEM analyses of kinetics (Anal. Biochem. 421 (2012) 181-190) are shown to be unfounded. Their assertion that their software allows "analysis of datasets too noisy to process by existing iterative search algorithms" is refuted with a MEM analysis of their triexponential test case with increased noise. In addition, it is shown that lifetime distributions recovered from noisy kinetics data with the MEM can be improved by using a simple filter when bootstrapping the prior model. When deriving the bootstrapped model from the lifetime distribution obtained using a uniform model, only the slower processes are represented as Gaussians in the bootstrapped model. Using this new approach, results are clearly superior to those of Mulligan et al. despite the presence of increased noise. In a second example, ambiguity in the interpretation of Poisson kinetics in the presence of scattered excitation light is resolved by filtering the prior model.

Antigen-antibody interactions and structural flexibility of a femtomolar-affinity antibody

The femtomolar-affinity mutant antibody (4M5.3) generated by directed evolution is interesting because of the potential of antibody engineering. In this study, the mutant and its wild type (4-4-20) were compared in terms of antigen-antibody interactions and structural flexibility to elucidate the effects of directed evolution. For this purpose, multiple steered molecular dynamics (SMD) simulations were performed. The pulling forces of SMD simulations elucidated the regions that form strong attractive interactions in the binding pocket. Structural analysis in these regions showed two important mutations for improving attractive interactions. First, mutation of Tyr102(H) to Ser (sequence numbering of Protein Data Bank entry 1FLR ) played a role in resolving the steric hindrance on the pathway of the antigen in the binding pocket. Second, mutation of Asp31(H) to His played a role in resolving electrostatic repulsion. Potentials of mean force (PMFs) of both the wild type and the mutant showed landscapes that do not include obvious intermediate states and go directly to the bound state. These landscapes were regarded as funnel-like binding free energy landscapes. Furthermore, the structural flexibility based on the fluctuations of the positions of atoms was analyzed. It was shown that the fluctuations in the positions of the antigen and residues in contact with antigen tend to be smaller in the mutant than in the wild type. This result suggested that structural flexibility decreases as affinity is improved by directed evolution. This suggestion is similar to the relationship between affinity and flexibility for in vivo affinity maturation, which was suggested by Romesberg and co-workers [Jimenez, R., et al. (2003) Proc. Natl. Acad. Sci. U.S.A.100, 92-97]. Consequently, the relationship was found to be applicable up to femotomolar affinity levels.