Molecular evolution of affinity and flexibility in the immune system

Accuracy and Reproducibility of Lipari-Szabo Order Parameters From Molecular Dynamics

The Lipari-Szabo generalized order parameter probes the picosecond to nanosecond time scale motions of a protein and is useful for rationalizing a multitude of biological processes such as protein recognition and ligand binding. Although these fast motions are an important and intrinsic property of proteins, it remains unclear what simulation conditions are most suitable to reproduce methyl symmetry axis side chain order parameter data (Saxis2) from molecular dynamics simulations. In this study, we show that, while Saxis2 tends to converge within tens of nanoseconds, it is essential to run 10 to 20 replicas starting from configurations close to the experimental structure to obtain the best agreement with experimental Saxis2 values. Additionally, in a comparison of force fields, AMBER ff14SB outperforms CHARMM36m in accurately capturing these fast time scale motions, and we suggest that the origin of this performance gap is likely attributed to differences in side chain torsional parametrization and not due to differences in the global protein conformations sampled by the force fields. This study provides insight into obtaining accurate and reproducible Saxis2 values from molecular simulations and underscores the necessity of using replica simulations to compute equilibrium properties.

A single-domain antibody library based on a stability-engineered human VH3 scaffold

Using conventional immunoglobulin G (IgG) molecules as therapeutic agents presents several well-known disadvantages owing to their large size and structural complexity, negatively impacting development and production efficiency. Single-domain antibodies (sdAbs) are the smallest functional antibody format (~ 15 kDa) and represent a viable alternative to IgG in many applications. However, unlike natural single-domain antibodies, such as camelid VHH, the variable domains of conventional antibodies show poor physicochemical properties when expressed as sdAbs. This report identified stable sdAb variants of human VH3-23 from a framework region 2-randomized human VH library by phage display selection under thermal challenge. Synthetic complementarity determining region diversity was introduced to one of the selected variants with high thermal stability, expression level, and monomeric content to construct a human VH sdAb library. The library was validated by panning against a panel of antigens, and target-specific binders were identified and characterized for their affinity and biophysical properties. The results of this study suggest that a synthetic sdAb library based on a stability-engineered human VH scaffold could be a facile source of high-quality sdAb for many practical applications.

Structural basis of broad SARS-CoV-2 cross-neutralization by affinity-matured public antibodies

Descendants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant now account for almost all SARS-CoV-2 infections. The Omicron variant and its sublineages have spike glycoproteins that are highly diverged from the pandemic founder and first-generation vaccine strain, resulting in significant evasion from monoclonal antibody therapeutics and vaccines. Understanding how commonly elicited antibodies can broaden to cross-neutralize escape variants is crucial. We isolate IGHV3-53, using "public" monoclonal antibodies (mAbs) from an individual 7 months post infection with the ancestral virus and identify antibodies that exhibit potent and broad cross-neutralization, extending to the BA.1, BA.2, and BA.4/BA.5 sublineages of Omicron. Deep mutational scanning reveals these mAbs' high resistance to viral escape. Structural analysis via cryoelectron microscopy of a representative broadly neutralizing antibody, CAB-A17, in complex with the Omicron BA.1 spike highlights the structural underpinnings of this broad neutralization. By reintroducing somatic hypermutations into a germline-reverted CAB-A17, we delineate the role of affinity maturation in the development of cross-neutralization by a public class of antibodies.

"Redirecting an anti-IL-1β antibody to bind a new, unrelated and computationally predicted epitope on hIL-17A"

Antibody engineering technology is at the forefront of therapeutic antibody development. The primary goal for engineering a therapeutic antibody is the generation of an antibody with a desired specificity, affinity, function, and developability profile. Mature antibodies are considered antigen specific, which may preclude their use as a starting point for antibody engineering. Here, we explore the plasticity of mature antibodies by engineering novel specificity and function to a pre-selected antibody template. Using a small, focused library, we engineered AAL160, an anti-IL-1β antibody, to bind the unrelated antigen IL-17A, with the introduction of seven mutations. The final redesigned antibody, 11.003, retains favorable biophysical properties, binds IL-17A with sub-nanomolar affinity, inhibits IL-17A binding to its cognate receptor and is functional in a cell-based assay. The epitope of the engineered antibody can be computationally predicted based on the sequence of the template antibody, as is confirmed by the crystal structure of the 11.003/IL-17A complex. The structures of the 11.003/IL-17A and the AAL160/IL-1β complexes highlight the contribution of germline residues to the paratopes of both the template and re-designed antibody. This case study suggests that the inherent plasticity of antibodies allows for re-engineering of mature antibodies to new targets, while maintaining desirable developability profiles.

Arginine cluster introduction on framework region in anti-lysozyme antibody improved association rate constant by changing conformational diversity of CDR loops

Antibodies are used for many therapeutic and biotechnological purposes. Because the affinity of an antibody to the antigen is critical for clinical efficacy of pharmaceuticals, many affinity maturation strategies have been developed. Although we previously reported an affinity maturation strategy in which the association rate of the antibody toward its antigen is improved by introducing a cluster of arginine residues into the framework region of the antibody, the detailed molecular mechanism responsible for this improvement has been unknown. In this study, we introduced five arginine residues into an anti-hen egg white lysozyme antibody (HyHEL10) Fab fragment to create the R5-mutant and comprehensively characterized the interaction between antibody and antigen using thermodynamic analysis, X-ray crystallography, and molecular dynamics (MD) simulations. Our results indicate that introduction of charged residues strongly enhanced the association rate, as previously reported, and the antibody-antigen complex structure was almost the same for the R5-mutant and wild-type Fabs. The MD simulations indicate that the mutation increased conformational diversity in complementarity-determining region loops and thereby enhanced the association rate. These observations provide the molecular basis of affinity maturation by R5 mutation.

Defining and Studying B Cell Receptor and TCR Interactions

BCRs (Abs) and TCRs (or adaptive immune receptors [AIRs]) are the means by which the adaptive immune system recognizes foreign and self-antigens, playing an integral part in host defense, as well as the emergence of autoimmunity. Importantly, the interaction between AIRs and their cognate Ags defies a simple key-in-lock paradigm and is instead a complex many-to-many mapping between an individual's massively diverse AIR repertoire, and a similarly diverse antigenic space. Understanding how adaptive immunity balances specificity with epitopic coverage is a key challenge for the field, and terms such as broad specificity, cross-reactivity, and polyreactivity remain ill-defined and are used inconsistently. In this Immunology Notes and Resources article, a group of experimental, structural, and computational immunologists define commonly used terms associated with AIR binding, describe methodologies to study these binding modes, as well as highlight the implications of these different binding modes for therapeutic design.

Surface patches induce nonspecific binding and phase separation of antibodies

Nonspecific interactions are a key challenge in the successful development of therapeutic antibodies. The tendency for nonspecific binding of antibodies is often difficult to reduce by rational design, and instead, it is necessary to rely on comprehensive screening campaigns. To address this issue, we performed a systematic analysis of the impact of surface patch properties on antibody nonspecificity using a designer antibody library as a model system and single-stranded DNA as a nonspecificity ligand. Using an in-solution microfluidic approach, we find that the antibodies tested bind to single-stranded DNA with affinities as high as KD = 1 µM. We show that DNA binding is driven primarily by a hydrophobic patch in the complementarity-determining regions. By quantifying the surface patches across the library, the nonspecific binding affinity is shown to correlate with a trade-off between the hydrophobic and total charged patch areas. Moreover, we show that a change in formulation conditions at low ionic strengths leads to DNA-induced antibody phase separation as a manifestation of nonspecific binding at low micromolar antibody concentrations. We highlight that phase separation is driven by a cooperative electrostatic network assembly mechanism of antibodies with DNA, which correlates with a balance between positive and negative charged patches. Importantly, our study demonstrates that both nonspecific binding and phase separation are controlled by the size of the surface patches. Taken together, these findings highlight the importance of surface patches and their role in conferring antibody nonspecificity and its macroscopic manifestation in phase separation.

New insights into the effect of mutations on affibody-Fc interaction, a molecular dynamics simulation approach

Staphylococcal protein A (SpA) domain B (the basis of affibody) has been widely used in affinity chromatography and found therapeutic applications against inflammatory diseases through targeting the Fc part of immunoglobulin G (IgG). We have performed extensive molecular dynamics simulation of 41 SpA mutants and compared their dynamics and conformations to wild type. The simulations revealed the molecular details of structural and dynamics changes that occurred due to introducing point mutations and helped to explain the SPR results. It was observed in some variants a point mutation caused extensive structural changes far from the mutation site, while an effect of some other mutations was limited to the site of the mutated residue. Also, the pattern of hydrogen bond networks and hydrophobic core arrangements were investigated. We figured out mutations that occurred at positions 128, 136, 150 and 153, affected two hydrophobic cores at the interface as well as mutations introduced at positions 129 and 154 interrupted two hydrogen bond networks of the interface, SPR data showed all of these mutations reduced binding affinity significantly. Overall, by scanning the SpA-Fc interface through the large numbers of introduced mutations, the new insights have been gained which would help to design high- affinity ligands of IgG.

Common framework mutations impact antibody interfacial dynamics and flexibility

Introduction

With the flood of engineered antibodies, there is a heightened need to elucidate the structural features of antibodies that contribute to specificity, stability, and breadth. While antibody flexibility and interface angle have begun to be explored, design rules have yet to emerge, as their impact on the metrics above remains unclear. Furthermore, the purpose of framework mutations in mature antibodies is highly convoluted.

Methods

To this end, a case study utilizing molecular dynamics simulations was undertaken to determine the impact framework mutations have on the VH-VL interface. We further sought to elucidate the governing mechanisms by which changes in the VH-VL interface angle impact structural elements of mature antibodies by looking at root mean squared deviations, root mean squared fluctuations, and solvent accessible surface area.

Results and discussion

Overall, our results suggest framework mutations can significantly shift the distribution of VH-VL interface angles, which leads to local changes in antibody flexibility through local changes in the solvent accessible surface area. The data presented herein highlights the need to reject the dogma of static antibody crystal structures and exemplifies the dynamic nature of these proteins in solution. Findings from this work further demonstrate the importance of framework mutations on antibody structure and lay the foundation for establishing design principles to create antibodies with increased specificity, stability, and breadth.

Molecular Dynamics Simulation for Investigating Antigen-Antibody Interaction

Molecular dynamics (MD) simulation is a computational method which elucidates the protein dynamics. Following analyses characterize the dynamics and structural change as well as interaction energy. To characterize the protein structure effectively, the internal angular coordinates are often useful. Directional analysis provides the averages and variances of those coordinates in a mathematically rigorous way. Here, we describe not only a standard MD simulation procedure for the antigen-antibody system but also an umbrella sampling method following a multistep targeted MD simulation (US/mTMD), which is useful for evaluating the free energy profile along the antigen-antibody dissociation coordinate.

Allosteric regulation in STAT3 interdomains is mediated by a rigid core: SH2 domain regulation by CCD in D170A variant

Signal Transducer and Activator of Transcription 3 (STAT3) plays a crucial role in cancer development and thus is a viable target for cancer treatment. STAT3 functions as a dimer mediated by phosphorylation of the SRC-homology 2 (SH2) domain, a key target for therapeutic drugs. While great efforts have been employed towards the development of compounds that directly target the SH2 domain, no compound has yet been approved by the FDA due to a lack of specificity and pharmacologic efficacy. Studies have shown that allosteric regulation of SH2 via the coiled-coil domain (CCD) is an alternative drug design strategy. Several CCD effectors have been shown to modulate SH2 binding and affinity, and at the time of writing at least one drug candidate has entered phase I clinical trials. However, the mechanism for SH2 regulation via CCD is poorly understood. Here, we investigate structural and dynamic features of STAT3 and compare the wild type to the reduced function variant D170A in order to delineate mechanistic differences and propose allosteric pathways. Molecular dynamics simulations were employed to explore conformational space of STAT3 and the variant, followed by structural, conformation, and dynamic analysis. The trajectories explored show distinctive conformational changes in the SH2 domain for the D170A variant, indicating long range allosteric effects. Multiple analyses provide evidence for long range communication pathways between the two STAT3 domains, which seem to be mediated by a rigid core which connects the CCD and SH2 domains via the linker domain (LD) and transmits conformational changes through a network of short-range interactions. The proposed allosteric mechanism provides new insight into the understanding of intramolecular signaling in STAT3 and potential pharmaceutical control of STAT3 specificity and activity.

Non-specificity as the sticky problem in therapeutic antibody development

Antibodies are highly potent therapeutic scaffolds with more than a hundred different products approved on the market. Successful development of antibody-based drugs requires a trade-off between high target specificity and target binding affinity. In order to better understand this problem, we here review non-specific interactions and explore their fundamental physicochemical origins. We discuss the role of surface patches - clusters of surface-exposed amino acid residues with similar physicochemical properties - as inducers of non-specific interactions. These patches collectively drive interactions including dipole-dipole, π-stacking and hydrophobic interactions to complementary moieties. We elucidate links between these supramolecular assembly processes and macroscopic development issues, such as decreased physical stability and poor in vivo half-life. Finally, we highlight challenges and opportunities for optimizing protein binding specificity and minimizing non-specificity for future generations of therapeutics.

Changes in antibody binding and functionality after humanizing a murine scFv anti-IFN-α2: From in silico studies to experimental analysis

The structural and dynamic changes introduced during antibody humanization continue to be a topic open to new contributions. For this reason, the study of structural and functional changes of a murine scFv (mu.scFv) anti-rhIFN-α2b after humanization was carried out. As it was shown by long molecular dynamics simulations and circular dichroism analysis, changes in primary sequence affected the tertiary structure of the humanized scFv (hz.scFv): the position of the variable domain of light chain (VL) respective to the variable domain of heavy chain (VH) in each scFv molecule was different. This change mainly impacted on conformation and dynamics of the complementarity-determining region 3 of VH (CDR-H3) which led to changes in the specificity and affinity of humanized scFv (hz.scFv). These observations agree with experimental results that showed a decrease in the antigen-binding strength of hz.scFv, and different capacities of these molecules to neutralize the in vitro rhIFN-α2b biological activity. Besides, experimental studies to characterize antigen-antibody binding showed that mu.scFv and hz.scFv bind to the same antigen area and recognize a conformational epitope, which is evidence of docking results. Finally, the differences between these molecules to neutralize the in vitro rhIFN-α2b biological activity were described as a consequence of the blockade of certain functionally relevant amino acids of the cytokine, after scFv binding. All these observations confirmed that humanization affected the affinity and specificity of hz.scFv and pointed out that two specific changes in the frameworks would be responsible.

From affinity selection to kinetic selection in Germinal Centre modelling

Affinity maturation is an evolutionary process by which the affinity of antibodies (Abs) against specific antigens (Ags) increases through rounds of B-cell proliferation, somatic hypermutation, and positive selection in germinal centres (GC). The positive selection of B cells depends on affinity, but the underlying mechanisms of affinity discrimination and affinity-based selection are not well understood. It has been suggested that selection in GC depends on both rapid binding of B-cell receptors (BcRs) to Ags which is kinetically favourable and tight binding of BcRs to Ags, which is thermodynamically favourable; however, it has not been shown whether a selection bias for kinetic properties is present in the GC. To investigate the GC selection bias towards rapid and tight binding, we developed an agent-based model of GC and compared the evolution of founder B cells with initially identical low affinities but with different association/dissociation rates for Ag presented by follicular dendritic cells in three Ag collection mechanisms. We compared an Ag collection mechanism based on association/dissociation rates of B-cell interaction with presented Ag, which includes a probabilistic rupture of bonds between the B-cell and Ag (Scenario-1) with a reference scenario based on an affinity-based Ag collection mechanism (Scenario-0). Simulations showed that the mechanism of Ag collection affects the GC dynamics and the GC outputs concerning fast/slow (un)binding of B cells to FDC-presented Ags. In particular, clones with lower dissociation rates outcompete clones with higher association rates in Scenario-1, while remaining B cells from clones with higher association rates reach higher affinities. Accordingly, plasma cell and memory B cell populations were biased towards B-cell clones with lower dissociation rates. Without such probabilistic ruptures during the Ag extraction process (Scenario-2), the selective advantage for clones with very low dissociation rates diminished, and the affinity maturation level of all clones decreased to the reference level.

Adaptive immunity is one of the vital defence mechanisms of the human body to fight virtually unlimited types of pathogens by producing antigen-specific high-affinity antibodies that bind to pathogens and neutralise them or mark them for further elimination. Affinity is a quantity used to measure and report the strength of interaction between antibodies and antigens that depends both on how fast antibodies bind to antigens (association rate) and how long the bond lasts (dissociation rate). The affinity of produced antibodies for a specific antigen increases in germinal centres through a process called affinity maturation, during which B cells with higher affinities have a competitive advantage and get positively selected to differentiate to antibody-producing plasma cells. Our research shows that the mechanism by which B cells capture Ag affects GC dynamics and GC output with respect to B-cell receptor kinetics. Notably, in a mechanism where rupture of CC-FDC bonds is possible during Ag extraction, B-cell clones with low dissociation rates outcompete clones with high association rates over time. Understanding how B cells get selected in germinal centres could help to develop an optimised and effective immune response against a disease through vaccination for a fast-operating and long-lasting immune response.

Conformational plasticity of the HIV-1 gp41 immunodominant region is recognized by multiple non-neutralizing antibodies

The early humoral immune response to acute HIV-1 infection is largely non-neutralizing. The principal target of these antibodies is the primary immunodominant region (PID) on the gp41 fusion protein. The PID is a highly conserved 15-residue region displayed on the surface of HIV-1 virions. In this study, we analyzed the humoral determinants of HIV-1 gp41 PID binding using biophysical, structural, and computational methods. In complex with a patient-derived near-germline antibody fragment, the PID motif adopts an elongated random coil, whereas the PID bound to affinity-matured Fab adopts a strand-turn-helix conformation. Molecular dynamics simulations showed that the PID is structurally plastic suggesting that the PID can form an ensemble of structural states recognized by various non-neutralizing antibodies, facilitating HIV-1 immunodominance observed in acute and chronic HIV-1 infections. An improved understanding of how the HIV-1 gp41 PID misdirects the early humoral response should guide the development of an effective HIV-1 vaccine.

The 15-amino-acid primary immunodominant (PID) region on HIV-1 gp41 adopts an ensemble of conformational states. This conformational plasticity is suggested to misdirect the early humoral immune response.

Selective Recognition of Carbohydrate Antigens by Germline Antibodies Isolated from AID Knockout Mice

Germline antibodies, the initial set of antibodies produced by the immune system, are critical for host defense, and information about their binding properties can be useful for designing vaccines, understanding the origins of autoantibodies, and developing monoclonal antibodies. Numerous studies have found that germline antibodies are polyreactive with malleable, flexible binding pockets. While insightful, it remains unclear how broadly this model applies, as there are many families of antibodies that have not yet been studied. In addition, the methods used to obtain germline antibodies typically rely on assumptions and do not work well for many antibodies. Herein, we present a distinct approach for isolating germline antibodies that involves immunizing activation-induced cytidine deaminase (AID) knockout mice. This strategy amplifies antigen-specific B cells, but somatic hypermutation does not occur because AID is absent. Using synthetic haptens, glycoproteins, and whole cells, we obtained germline antibodies to an assortment of clinically important tumor-associated carbohydrate antigens, including Lewis Y, the Tn antigen, sialyl Lewis C, and Lewis X (CD15/SSEA-1). Through glycan microarray profiling and cell binding, we demonstrate that all but one of these germline antibodies had high selectivity for their glycan targets. Using molecular dynamics simulations, we provide insights into the structural basis of glycan recognition. The results have important implications for designing carbohydrate-based vaccines, developing anti-glycan monoclonal antibodies, and understanding antibody evolution within the immune system.

Crystallography-based exploration of non-covalent interactions for the design and synthesis of coumarin for stronger protein binding

Protein molecules are a good target for the inhibition or promotion of biological processes. Different methods like QSAR and molecular docking have been developed to accurately design small binder molecules for target proteins. An alternative model has been developed wherein a statistical method is used to find the propensity of different non-covalent interactions between small molecules and amino acid residues of the protein. The results give hints as to the choice of substituents required at the SM to strongly bind to a protein. In this case, 75 different types of proteins bound with coumarin derivatives have been investigated and the non-covalent interactions observed between the basic coumarin moiety and amino acids have been analyzed. Density functional theory (DFT) calculations were used to identify the electronic features of coumarin to understand the feasibility of the observed non-covalent interactions and to find appropriate groups that can modulate these interactions. The binding affinity towards a protein (β-lactoglobulin (BLG)) and the stability of the protein complex have been investigated through docking and molecular dynamics of 100 ns, respectively. The modeled compounds were synthesized and investigated with regards to their interactions with the model carrier protein. The thermodynamics of the interactions were also investigated and the binding is governed by the Le Chatelier principle.

The role of cholesterol recognition (CARC/CRAC) mirror codes in the allosterism of the human organic cation transporter 2 (OCT2, SLC22A2)

The human organic cation transporter 2 (OCT2) is a multispecific transporter with cholesterol-dependent allosteric features. The present work elucidates the role of evolutionarily conserved cholesterol recognition/interaction amino acid consensus sequences (CRAC and CARC) in the allosteric binding to 1-methyl-4-phenylpyridinium (MPP⁺) in human embryonic kidney 293 cells stably or transiently expressing OCT2. Molecular blind simulations docked two mirroring cholesterol molecules in the 5th putative transmembrane domain, where a CARC and a CRAC sequence lie. The impact of the conserved amino acids that may constitute the CARC/CRAC mirror code was studied by alanine-scanning mutagenesis. At a saturating extracellular concentration of substrate, at which the impact of cholesterol depletion is maximal, five mutants transported MPP⁺ at a significantly lower rate than the wild-type OCT2 (WT), resembling the behavior of the WT upon cholesterol depletion. MPP⁺ influx rate as a function of the extracellular concentration of substrate was measured for the mutants R234A, R235A, L252A and R263A. R234A kinetic behavior was similar to that of the WT, whereas R235A, L252A and R263A activity shifted from allosteric to one-binding site kinetics, very much like the WT upon cholesterol depletion. The impact of cholesterol on protein thermal stability was assessed for WT, R234A and R263A. While the thermal stability of WT and R234A was improved by the supplementation with cholesterol, R263A was not sensitive to the presence of cholesterol. To conclude, the disruption of the CARC/CRAC mirror code in the 5th putative transmembrane domain is sufficient to abolish the allosteric interaction between OCT2 and MPP⁺.

Delicate balance among thermal stability, binding affinity, and conformational space explored by single-domain VHH antibodies

The high binding affinities and specificities of antibodies have led to their use as drugs and biosensors. Single-domain V_HH antibodies exhibit high specificity and affinity but have higher stability and solubility than conventional antibodies as they are single-domain proteins. In this work, based on physicochemical measurements and molecular dynamics (MD) simulations, we have gained insight that will facilitate rational design of single-chain V_HH antibodies. We first assessed two homologous V_HH antibodies by differential scanning calorimetry (DSC); one had a high (64.8 °C) and the other a low (58.6 °C) melting temperature. We then generated a series of the variants of the low stability antibody and analyzed their thermal stabilities by DSC and characterized their structures through MD simulations. We found that a single mutation that resulted in 8.2 °C improvement in melting temperature resulted in binding affinity an order of magnitude lower than the parent antibody, likely due to a shift of conformational space explored by the single-chain V_HH antibody. These results suggest that the delicate balance among conformational stability, binding capability, and conformational space explored by antibodies must be considered in design of fully functional single-chain V_HH antibodies.

Polyclonal antibody responses to HIV Env immunogens resolved using cryoEM

Engineered ectodomain trimer immunogens based on BG505 envelope glycoprotein are widely utilized as components of HIV vaccine development platforms. In this study, we used rhesus macaques to evaluate the immunogenicity of several stabilized BG505 SOSIP constructs both as free trimers and presented on a nanoparticle. We applied a cryoEM-based method for high-resolution mapping of polyclonal antibody responses elicited in immunized animals (cryoEMPEM). Mutational analysis coupled with neutralization assays were used to probe the neutralization potential at each epitope. We demonstrate that cryoEMPEM data can be used for rapid, high-resolution analysis of polyclonal antibody responses without the need for monoclonal antibody isolation. This approach allowed to resolve structurally distinct classes of antibodies that bind overlapping sites. In addition to comprehensive mapping of commonly targeted neutralizing and non-neutralizing epitopes in BG505 SOSIP immunogens, our analysis revealed that epitopes comprising engineered stabilizing mutations and of partially occupied glycosylation sites can be immunogenic.

Binding Insight of Anti-HIV Phytocompounds with Prime Targets of HIV: A Molecular Dynamics Simulation Analysis

BACKGROUND

Despite intense efforts, AIDS is difficult to tackle by current anti-retroviral therapy (ART) due to its side effects; therefore, there is an urgent need to discover potential, multitarget and low-cost anti-HIV compounds.

OBJECTIVE

We have shown that few phytocompounds can potentially inhibit the prime targets of HIV namely GP120 envelope protein, reverse transcriptase, protease, integrase and ribonulcease. In this study, top ranked prioritized compounds were subjected to Molecular Dynamics (MD) simulation in order to study the conformational dynamics and integrity of crucial interaction in the receptor sites.

METHODS

The system was built for selected protein-ligand complex using TIP3P water model and OPLS_2005 force field. Trajectories were recorded up to 20 ns simulation time in Desmond module of Schrödinger software.

RESULTS

As a result of a comprehensive analysis of molecular properties and dynamics of the complexes, it has been concluded that Chebulic acid, Curcumin and Mulberroside C could be developed as envelope glycoprotein GP120 inhibitor, reverse transcriptase inhibitor and protease inhibitor respectively. However, the fluctuation of Chebulic acid with respect to integrase and ribonuclease protein was higher during the simulation.

CONCLUSION

These findings can aid in the designing of the structural properties for more effective anti-HIV compounds against the given targets.

Ensembles in solution as a new paradigm for antibody structure prediction and design

The rise of antibodies as a promising and rapidly growing class of biotherapeutic proteins has motivated numerous studies to characterize and understand antibody structures. In the past decades, the number of antibody crystal structures increased substantially, which revolutionized the atomistic understanding of antibody functions. Even though numerous static structures are known, various biophysical properties of antibodies (i.e., specificity, hydrophobicity and stability) are governed by their dynamic character. Additionally, the importance of high-quality structures in structure–function relationship studies has substantially increased. These structure–function relationship studies have also created a demand for precise homology models of antibody structures, which allow rational antibody design and engineering when no crystal structure is available. Here, we discuss various aspects and challenges in antibody design and extend the paradigm of describing antibodies with only a single static structure to characterizing them as dynamic ensembles in solution.

Parameters and determinants of responses to selection in antibody libraries

The sequences of antibodies from a given repertoire are highly diverse at few sites located on the surface of a genome-encoded larger scaffold. The scaffold is often considered to play a lesser role than highly diverse, non-genome-encoded sites in controlling binding affinity and specificity. To gauge the impact of the scaffold, we carried out quantitative phage display experiments where we compare the response to selection for binding to four different targets of three different antibody libraries based on distinct scaffolds but harboring the same diversity at randomized sites. We first show that the response to selection of an antibody library may be captured by two measurable parameters. Second, we provide evidence that one of these parameters is determined by the degree of affinity maturation of the scaffold, affinity maturation being the process by which antibodies accumulate somatic mutations to evolve towards higher affinities during the natural immune response. In all cases, we find that libraries of antibodies built around maturated scaffolds have a lower response to selection to other arbitrary targets than libraries built around germline-based scaffolds. We thus propose that germline-encoded scaffolds have a higher selective potential than maturated ones as a consequence of a selection for this potential over the long-term evolution of germline antibody genes. Our results are a first step towards quantifying the evolutionary potential of biomolecules.

Antibodies in the immune system consist of a genetically encoded scaffold that exposes a few highly diverse, non-genetically encoded sites. This focused diversity is sufficient to produce antibodies that bind to any target molecule. To understand the role of the scaffold, which acquires hypermutations during the immune response, over the selective response, we analyze quantitative in vitro experiments where large antibody populations based on different scaffolds are selected against different targets. We show that selective responses are described statistically by two parameters, one of which depends on prior evolution of the scaffold as part of a previous response. Our work provides methods to assay whether naïve antibody scaffolds are endowed with a distinctively high selective potential.

Comparison of Huggins Coefficients and Osmotic Second Virial Coefficients of Buffered Solutions of Monoclonal Antibodies

The Huggins coefficient kH is a well-known metric for quantifying the increase in solution viscosity arising from intermolecular interactions in relatively dilute macromolecular solutions, and there has been much interest in this solution property in connection with developing improved antibody therapeutics. While numerous kH measurements have been reported for select monoclonal antibodies (mAbs) solutions, there has been limited study of kH in terms of the fundamental molecular interactions that determine this property. In this paper, we compare measurements of the osmotic second virial coefficient B22, a common metric of intermolecular and interparticle interaction strength, to measurements of kH for model antibody solutions. This comparison is motivated by the seminal work of Russel for hard sphere particles having a short-range "sticky" interparticle interaction, and we also compare our data with known results for uncharged flexible polymers having variable excluded volume interactions because proteins are polypeptide chains. Our observations indicate that neither the adhesive hard sphere model, a common colloidal model of globular proteins, nor the familiar uncharged flexible polymer model, an excellent model of intrinsically disordered proteins, describes the dependence of kH of these antibodies on B22. Clearly, an improved understanding of protein and ion solvation by water as well as dipole-dipole and charge-dipole effects is required to understand the significance of kH from the standpoint of fundamental protein-protein interactions. Despite shortcomings in our theoretical understanding of kH for antibody solutions, this quantity provides a useful practical measure of the strength of interprotein interactions at elevated protein concentrations that is of direct significance for the development of antibody formulations that minimize the solution viscosity.

Probing of the plasticity of the active site in pinene synthase elucidates its potential evolutionary mechanism

Terpenes form a class of highly diverse natural products. The diversity of terpenes is created by terpene synthases. During the reaction, carbocation intermediates form and their rearrangement could lead to the formation of various products. Terpene synthases determine the final product profile by controlling the conformation of the intermediate or stabilizing the carbocation. Pinene synthase is a monoterpene synthase catalyzing the cyclization of geranyl pyrophosphate (GPP) to form pinene. Our study suggests that by mutating the aromatic residue F482 to Ala, Val, Ile and Thr, the enzyme can be converted to sabinene synthase, with more than 90% of its total terpene products becoming sabinene, which indicates the aromaticity of this residue is essential for stabilizing the pinyl carbocation. We also identified a mutation S491A that could cause an about 29% increase in the overall activity of the enzyme without altering its produce selectivity. Molecular dynamic simulation indicates this mutation could decrease the flexibility of the enzyme when it forms a complex with the pinyl carbocation. Our study suggests the active pocket of pinene synthase has a certain level of plasticity, making it relatively easy to change the product selectivity or overall activity. This property could have an important implication in the evolution of terpene synthases and thereby terpene diversity, as by changing a few residues an enzyme could synthesize a completely different terpene product in a significant amount, which allows selection to take place.

Surprisingly Fast Interface and Elbow Angle Dynamics of Antigen-Binding Fragments

Fab consist of a heavy and light chain and can be subdivided into a variable (V H and V L ) and a constant region (C H 1 and C L ). The variable region contains the complementarity-determining region (CDR), which is formed by six hypervariable loops, shaping the antigen binding site, the paratope. Apart from the CDR loops, both the elbow angle and the relative interdomain orientations of the V H -V L and the C H 1-C L domains influence the shape of the paratope. Thus, characterization of the interface and elbow angle dynamics is essential to antigen specificity. We studied nine antigen-binding fragments (Fab) to investigate the influence of affinity maturation, antibody humanization, and different light-chain types on the interface and elbow angle dynamics. While the CDR loops reveal conformational transitions in the micro-to-millisecond timescale, both the interface and elbow angle dynamics occur on the low nanosecond timescale. Upon affinity maturation, we observe a substantial rigidification of the V H and V L interdomain and elbow-angle flexibility, reflected in a narrower and more distinct distribution. Antibody humanization describes the process of grafting non-human CDR loops onto a representative human framework. As the antibody framework changes upon humanization, we investigated if both the interface and the elbow angle distributions are changed or shifted. The results clearly showed a substantial shift in the relative V H -V L distributions upon antibody humanization, indicating that different frameworks favor distinct interface orientations. Additionally, the interface and elbow angle dynamics of five antibody fragments with different light-chain types are included, because of their strong differences in elbow angles. For these five examples, we clearly see a high variability and flexibility in both interface and elbow angle dynamics, highlighting the fact that Fab interface orientations and elbow angles interconvert between each other in the low nanosecond timescale. Understanding how the relative interdomain orientations and the elbow angle influence antigen specificity, affinity, and stability has broad implications in the field of antibody modeling and engineering.

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.

A proton-coupled folate transporter mutation causing hereditary folate malabsorption locks the protein in an inward-open conformation

The proton-coupled folate transporter (PCFT, SLC46A1) is required for folate intestinal absorption and transport across the choroid plexus. Recent work has identified a F392V mutation causing hereditary folate malabsorption. However, the residue properties responsible for this loss of function remains unknown. Using site-directed mutagenesis, we observed complete loss of function with charged (Lys, Asp, and Glu) and polar (Thr, Ser, and Gln) Phe-392 substitutions and minimal function with some neutral substitutions; however, F392M retained full function. Using the substituted-cysteine accessibility method (with N-biotinyl aminoethyl methanethiosulfonate labeling), Phe-392 mutations causing loss of function, although preserving membrane expression and trafficking, also resulted in loss of accessibility of the substituted cysteine in P314C-PCFT located within the aqueous translocation pathway. F392V function and accessibility of the P314C cysteine were restored by insertion of a G305L (suppressor) mutation. A S196L mutation localized in proximity to Gly-305 by homology modeling was inactive. However, when inserted into the inactive F392V scaffold, function was restored (mutually compensatory mutations), as was accessibility of the P314C cysteine residue. Reduced function, documented with F392H PCFT, was due to a 15-fold decrease in methotrexate influx V _max, accompanied by a decreased influx K_t (4.5-fold) and K_i (3-fold). The data indicate that Phe-392 is required for rapid oscillation of the carrier among its conformational states and suggest that this is achieved by dampening affinity of the protein for its folate substrates. F392V and other inactivating Phe-392 PCFT mutations lock the protein in its inward-open conformation. Reach (length) and hydrophobicity of Phe-392 appear to be features required for full activity.

Local and Global Rigidification Upon Antibody Affinity Maturation

During the affinity maturation process the immune system produces antibodies with higher specificity and activity through various rounds of somatic hypermutations in response to an antigen. Elucidating the affinity maturation process is fundamental in understanding immunity and in the development of biotherapeutics. Therefore, we analyzed 10 pairs of antibody fragments differing in their specificity and in distinct stages of affinity maturation using metadynamics in combination with molecular dynamics (MD) simulations. We investigated differences in flexibility of the CDR-H3 loop and global changes in plasticity upon affinity maturation. Among all antibody pairs we observed a substantial rigidification in flexibility and plasticity reflected in a substantial decrease of conformational diversity. To visualize and characterize these findings we used Markov-states models to reconstruct the kinetics of CDR-H3 loop dynamics and for the first time provide a method to define and localize surface plasticity upon affinity maturation.

Structural basis for Zika envelope domain III recognition by a germline version of a recurrent neutralizing antibody

Recent epidemics demonstrate the global threat of Zika virus (ZIKV), a flavivirus transmitted by mosquitoes. Although infection is usually asymptomatic or mild, newborns of infected mothers can display severe symptoms, including neurodevelopmental abnormalities and microcephaly. Given the large-scale spread, symptom severity, and lack of treatment or prophylaxis, a safe and effective ZIKV vaccine is urgently needed. However, vaccine design is complicated by concern that elicited antibodies (Abs) may cross-react with other flaviviruses that share a similar envelope protein, such as dengue virus, West Nile virus, and yellow fever virus. This cross-reactivity may worsen symptoms of a subsequent infection through Ab-dependent enhancement. To better understand the neutralizing Ab response and risk of Ab-dependent enhancement, further information on germline Ab binding to ZIKV and the maturation process that gives rise to potently neutralizing Abs is needed. Here we use binding and structural studies to compare mature and inferred-germline Ab binding to envelope protein domain III of ZIKV and other flaviviruses. We show that affinity maturation of the light-chain variable domain is important for strong binding of the recurrent VH3-23/VK1-5 neutralizing Abs to ZIKV envelope protein domain III, and identify interacting residues that contribute to weak, cross-reactive binding to West Nile virus. These findings provide insight into the affinity maturation process and potential cross-reactivity of VH3-23/VK1-5 neutralizing Abs, informing precautions for protein-based vaccines designed to elicit germline versions of neutralizing Abs.

[Semi-rational Design of Target-binding Small Proteins for Cancer Treatment]

Small proteins that have a high affinity for cancer cell surface markers can be promising cheap alternatives to antibodies (antibody mimetics). Various types of antibody mimetics have thus been extensively developed. We recently found that a target-binding peptide binds to its target molecule more strongly when it is structurally constrained. To apply this finding to the development of chemically synthesizable small antibody mimetics, we have established an efficient method of creating such proteins, named fluctuation-regulated affinity proteins (FLAPs). To identify desirable scaffolds, first, 13 human proteins (46-104 aa) were selected from the Protein Data Bank. Then, thirteen graft acceptor (GA) sites that efficiently immobilize the grafted peptide structure were identified from six small protein scaffolds using molecular dynamics simulation. To assess the designed antibody mimetics in vitro, human epidermal growth factor receptor 2 (HER2)-binding peptides were selected from the anti-HER2 antibody drugs trastuzumab and pertuzumab by calculating the binding energy, and these were then grafted into the GA sites of scaffolds to create 65 FLAP candidates. The FLAP candidates were expressed in bacteria as fusion proteins with Renilla luciferase (Rluc), and their relative binding affinity to HER2 was easily determined by measuring the Rluc bioluminescence intensity without protein purification. Finally, four out of the 65 showed specific binding to HER2 with a dissociation constant (K_D) of 24-65 nM, and these were used for the detection of HER2-expressing cancer cells. Our design strategy will promote the development of antibody mimetics for the effective treatment of cancers and other diseases.

Design Strategy to Create Antibody Mimetics Harbouring Immobilised Complementarity Determining Region Peptides for Practical Use

Monoclonal antibodies (mAbs) are attractive therapeutics for treating a wide range of human disorders, and bind to the antigen through their complementarity-determining regions (CDRs). Small stable proteins containing structurally retained CDRs are promising alternatives to mAbs. In this report, we present a method to create such proteins, named fluctuation-regulated affinity proteins (FLAPs). Thirteen graft acceptor (GA) sites that efficiently immobilise the grafted peptide structure were initially selected from six small protein scaffolds by computational identification. Five CDR peptides extracted by binding energy calculations from mAbs against breast cancer marker human epithelial growth factor receptor type 2 (HER2) were then grafted to the selected scaffolds. The combination of five CDR peptides and 13 GA sites in six scaffolds revealed that three of the 65 combinations showed specific binding to HER2 with dissociation constants (K_D) of 270–350 nM in biolayer interferometry and 24–65 nM in ELISA. The FLAPs specifically detected HER2-overexpressing cancer cells. Thus, the present strategy is a promising and practical method for developing small antibody mimetics.

Can we use structural knowledge to design a protective vaccine against HIV-1?

Human immunodeficiency virus/acquired immunodeficiency syndrome (HIV/AIDS) remains one of the most important current threats to global public health. Since its identification in the early 1980s, at least 75 million people have been infected with HIV-1, the virus that causes AIDS. Although antiretroviral drugs are effective at prolonging life after infection in the developed world, they are associated with significant side effects and are not in widespread use in the developing world. The best way to control the AIDS epidemic would be a vaccine that protects against infection by HIV-1. Most vaccines work by inducing antibodies in serum or mucosa that block infection or prevent invasion of the bloodstream. Here, I describe background related to my laboratory's attempts to develop an immunogen that would elicit protective antibodies against HIV-1.

Affinity Maturation Enhances Antibody Specificity but Compromises Conformational Stability

Monoclonal antibodies (mAbs) have recently emerged as one of the most promising classes of biotherapeutics. A potential advantage of B cell-derived mAbs as therapeutic agents is that they have been subjected to natural filtering mechanisms, which may enrich for B cell receptors (BCRs) with favorable biophysical properties. Here, we evaluated 400 human mAbs for polyreactivity, hydrophobicity, and thermal stability using high-throughput screening assays. Overall, mAbs derived from memory B cells and long-lived plasma cells (LLPCs) display reduced levels of polyreactivity, hydrophobicity, and thermal stability compared with naive B cell-derived mAbs. Somatic hypermutation (SHM) is inversely associated with all three biophysical properties, as well as BCR expression levels. Finally, the developability profiles of the human B cell-derived mAbs are comparable with those observed for clinical mAbs, suggesting their high therapeutic potential. The results provide insight into the biophysical consequences of affinity maturation and have implications for therapeutic antibody engineering and development.

Toward rational antibody design: recent advancements in molecular dynamics simulations

Because antibodies have become an important therapeutic tool, rational antibody design is a challenging issue involving various science and technology fields. From the computational aspect, many types of design-assist methods have been developed, but their accuracy is not fully satisfactory. Because of recent advancements in computational power, molecular dynamics (MD) simulation has become a helpful tool to trace the motion of proteins and to characterize their properties. Thus, MD simulation has been applied to various systems involving antigen-antibody complexes and has been shown to provide accurate insight into antigen-antibody interactions and dynamics at an atomic resolution. Therefore, it is highly possible that MD simulation will play several roles complementing the conventional antibody design. In this review, we address several important features of MD simulation in the context of rational antibody design.

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.

The Noisy and Marvelous Molecular World of Biology

At the molecular level biology is intrinsically noisy. The forces that regulate the myriad of molecular reactions in the cell are tiny, on the order of piconewtons (10−12 Newtons), yet they proceed in concerted action making life possible. Understanding how this is possible is one of the most fundamental questions biophysicists would like to understand. Single molecule experiments offer an opportunity to delve into the fundamental laws that make biological complexity surface in a physical world governed by the second law of thermodynamics. Techniques such as force spectroscopy, fluorescence, microfluidics, molecular sequencing, and computational studies project a view of the biomolecular world ruled by the conspiracy between the disorganizing forces due to thermal motion and the cosmic evolutionary drive. Here we will digress on some of the evidences in support of this view and the role of physical information in biology.

Computational design of structured loops for new protein functions

The ability to engineer the precise geometries, fine-tuned energetics and subtle dynamics that are characteristic of functional proteins is a major unsolved challenge in the field of computational protein design. In natural proteins, functional sites exhibiting these properties often feature structured loops. However, unlike the elements of secondary structures that comprise idealized protein folds, structured loops have been difficult to design computationally. Addressing this shortcoming in a general way is a necessary first step towards the routine design of protein function. In this perspective, we will describe the progress that has been made on this problem and discuss how recent advances in the field of loop structure prediction can be harnessed and applied to the inverse problem of computational loop design.

Characterizing the Diversity of the CDR-H3 Loop Conformational Ensembles in Relationship to Antibody Binding Properties

We present an approach to assess antibody CDR-H3 loops according to their dynamic properties using molecular dynamics simulations. We selected six antibodies in three pairs differing substantially in their individual promiscuity respectively specificity. For two pairs of antibodies crystal structures are available in different states of maturation and used as starting structures for the analyses. For a third pair we chose two antibody CDR sequences obtained from a synthetic library and predicted the respective structures. For all three pairs of antibodies we performed metadynamics simulations to overcome the limitations in conformational sampling imposed by high energy barriers. Additionally, we used classic molecular dynamics simulations to describe nano- to microsecond flexibility and to estimate up to millisecond kinetics of captured conformational transitions. The methodology represents the antibodies as conformational ensembles and allows comprehensive analysis of structural diversity, thermodynamics of conformations and kinetics of structural transitions. Referring to the concept of conformational selection we investigated the link between promiscuity and flexibility of the antibodies' binding interfaces. The obtained detailed characterization of the binding interface clearly indicates a link between structural flexibility and binding promiscuity for this set of antibodies.

Repertoire Analysis of Antibody CDR-H3 Loops Suggests Affinity Maturation Does Not Typically Result in Rigidification

Antibodies can rapidly evolve in specific response to antigens. Affinity maturation drives this evolution through cycles of mutation and selection leading to enhanced antibody specificity and affinity. Elucidating the biophysical mechanisms that underlie affinity maturation is fundamental to understanding B-cell immunity. An emergent hypothesis is that affinity maturation reduces the conformational flexibility of the antibody's antigen-binding paratope to minimize entropic losses incurred upon binding. In recent years, computational and experimental approaches have tested this hypothesis on a small number of antibodies, often observing a decrease in the flexibility of the complementarity determining region (CDR) loops that typically comprise the paratope and in particular the CDR-H3 loop, which contributes a plurality of antigen contacts. However, there were a few exceptions and previous studies were limited to a small handful of cases. Here, we determined the structural flexibility of the CDR-H3 loop for thousands of recent homology models of the human peripheral blood cell antibody repertoire using rigidity theory. We found no clear delineation in the flexibility of naïve and antigen-experienced antibodies. To account for possible sources of error, we additionally analyzed hundreds of human and mouse antibodies in the Protein Data Bank through both rigidity theory and B-factor analysis. By both metrics, we observed only a slight decrease in the CDR-H3 loop flexibility when comparing affinity matured antibodies to naïve antibodies, and the decrease was not as drastic as previously reported. Further analysis, incorporating molecular dynamics simulations, revealed a spectrum of changes in flexibility. Our results suggest that rigidification may be just one of many biophysical mechanisms for increasing affinity.

Role of framework mutations and antibody flexibility in the evolution of broadly neutralizing antibodies

Eliciting antibodies that are cross reactive with surface proteins of diverse strains of highly mutable pathogens (e.g., HIV, influenza) could be key for developing effective universal vaccines. Mutations in the framework regions of such broadly neutralizing antibodies (bnAbs) have been reported to play a role in determining their properties. We used molecular dynamics simulations and models of affinity maturation to study specific bnAbs against HIV. Our results suggest that there are different classes of evolutionary lineages for the bnAbs. If germline B cells that initiate affinity maturation have high affinity for the conserved residues of the targeted epitope, framework mutations increase antibody rigidity as affinity maturation progresses to evolve bnAbs. If the germline B cells exhibit weak/moderate affinity for conserved residues, an initial increase in flexibility via framework mutations may be required for the evolution of bnAbs. Subsequent mutations that increase rigidity result in highly potent bnAbs. Implications of our results for immunogen design are discussed.

Antibody Affinity Maturation by Computational Design

The immune systems protect our bodies from foreign molecules or antigens, where antibodies play important roles. Antibodies evolve over time upon antigen encounter by somatically mutating their genome sequences. The end result is a series of antibodies that display higher affinities and specificities to specific antigens. This process is called affinity maturation. Recent improvements in computer hardware and modeling algorithms now enable the rational design of protein structures and functions, and several works on computer-aided antibody design have been published. In this chapter, we briefly describe computational methods for antibody affinity maturation, focusing on methods for sampling antibody conformations and for scoring designed antibody variants. We also discuss lessons learned from the successful computer-aided design of antibodies.

Dynamics behind affinity maturation of an anti-HCMV antibody family influencing antigen binding

The investigation of antibody affinity maturation and its effects on antigen binding is important with respect to understanding the regulation of the immune response. To shed light on this crucial process, we analyzed two Igs neutralizing the human cytomegalovirus: the primary germline antibody M2J1 and its related mature antibody 8F9. Both antibodies target the AD-2S1 epitope of the gB envelope protein and are considered to establish similar interactions with the cognate antigen. We used molecular dynamics simulations to understand the effect of mutations on the antibody-antigen interactions. The results provide a qualitative explanation for the increased 8F9 peptide affinity compared with that of M2J1. The emerging atomistic-detailed description of these complexes reveals the molecular effects of the somatic hypermutations occurring during affinity maturation.

Genetic and structural analyses of affinity maturation in the humoral response to HIV-1

Most broadly neutralizing antibodies (BNAbs) elicited in response to HIV-1 infection are extraordinarily mutated. One goal of HIV-1 vaccine development is to induce antibodies that are similar to the most potent and broad BNAbs isolated from infected subjects. The most effective BNAbs have very high mutation frequencies, indicative of the long periods of continual activation necessary to acquire the BNAb phenotype through affinity maturation. Understanding the mutational patterns that define the maturation pathways in BNAb development is critical to vaccine design efforts to recapitulate through vaccination the successful routes to neutralization breadth and potency that have occurred in natural infection. Studying the mutational changes that occur during affinity maturation, however, requires accurate partitioning of sequence data into B-cell clones and identification of the starting point of a B-cell clonal lineage, the initial V(D)J rearrangement. Here, we describe the statistical framework we have used to perform these tasks. Through the recent advancement of these and similar computational methods, many HIV-1 ancestral antibodies have been inferred, synthesized and their structures determined. This has allowed, for the first time, the investigation of the structural mechanisms underlying the affinity maturation process in HIV-1 antibody development. Here, we review what has been learned from this atomic-level structural characterization of affinity maturation in HIV-1 antibodies and the implications for vaccine design.

Computer Simulations of Intrinsically Disordered Proteins

The investigation of intrinsically disordered proteins (IDPs) is a new frontier in structural and molecular biology that requires a new paradigm to connect structural disorder to function. Molecular dynamics simulations and statistical thermodynamics potentially offer ideal tools for atomic-level characterizations and thermodynamic descriptions of this fascinating class of proteins that will complement experimental studies. However, IDPs display sensitivity to inaccuracies in the underlying molecular mechanics force fields. Thus, achieving an accurate structural characterization of IDPs via simulations is a challenge. It is also daunting to perform a configuration-space integration over heterogeneous structural ensembles sampled by IDPs to extract, in particular, protein configurational entropy. In this review, we summarize recent efforts devoted to the development of force fields and the critical evaluations of their performance when applied to IDPs. We also survey recent advances in computational methods for protein configurational entropy that aim to provide a thermodynamic link between structural disorder and protein activity.

Structure-Based Reverse Vaccinology Failed in the Case of HIV Because it Disregarded Accepted Immunological Theory

Two types of reverse vaccinology (RV) should be distinguished: genome-based RV for bacterial vaccines and structure-based RV for viral vaccines. Structure-based RV consists in trying to generate a vaccine by first determining the crystallographic structure of a complex between a viral epitope and a neutralizing monoclonal antibody (nMab) and then reconstructing the epitope by reverse molecular engineering outside the context of the native viral protein. It is based on the unwarranted assumption that the epitope designed to fit the nMab will have acquired the immunogenic capacity to elicit a polyclonal antibody response with the same protective capacity as the nMab. After more than a decade of intensive research using this type of RV, this approach has failed to deliver an effective, preventive HIV-1 vaccine. The structure and dynamics of different types of HIV-1 epitopes and of paratopes are described. The rational design of an anti-HIV-1 vaccine is shown to be a misnomer since investigators who claim that they design a vaccine are actually only improving the antigenic binding capacity of one epitope with respect to only one paratope and not the immunogenic capacity of an epitope to elicit neutralizing antibodies. Because of the degeneracy of the immune system and the polyspecificity of antibodies, each epitope studied by the structure-based RV procedure is only one of the many epitopes that the particular nMab is able to recognize and there is no reason to assume that this nMab must have been elicited by this one epitope of known structure. Recent evidence is presented that the trimeric Env spikes of the virus possess such an enormous plasticity and intrinsic structural flexibility that it is it extremely difficult to determine which Env regions are the best candidate vaccine immunogens most likely to elicit protective antibodies.