Probing antibody internal dynamics with fluorescence anisotropy and molecular dynamics simulations

Viscosity Prediction of High-Concentration Antibody Solutions with Atomistic Simulations

The computational prediction of the viscosity of dense protein solutions is highly desirable, for example, in the early development phase of high-concentration biopharmaceutical formulations where the material needed for experimental determination is typically limited. Here, we use large-scale atomistic molecular dynamics (MD) simulations with explicit solvation to de novo predict the dynamic viscosities of solutions of a monoclonal IgG1 antibody (mAb) from the pressure fluctuations using a Green-Kubo approach. The viscosities at simulated mAb concentrations of 200 and 250 mg/mL are compared to the experimental values, which we measured with rotational rheometry. The computational viscosity of 24 mPa·s at the mAb concentration of 250 mg/mL matches the experimental value of 23 mPa·s obtained at a concentration of 213 mg/mL, indicating slightly different effective concentrations (or activities) in the MD simulations and in the experiments. This difference is assigned to a slight underestimation of the effective mAb-mAb interactions in the simulations, leading to a too loose dynamic mAb network that governs the viscosity. Taken together, this study demonstrates the feasibility of all-atom MD simulations for predicting the properties of dense mAb solutions and provides detailed microscopic insights into the underlying molecular interactions. At the same time, it also shows that there is room for further improvements and highlights challenges, such as the massive sampling required for computing collective properties of dense biomolecular solutions in the high-viscosity regime with reasonable statistical precision.

How neutron scattering techniques benefit investigating structures and dynamics of monoclonal antibody

Over the past several decades, great progresses have been made for the pharmaceutical industry of monoclonal antibody (mAb). More and more mAb products were approved for human therapeutics. This review describes the state of art of utilizing neutron scattering to investigate mAbs, in the aspects of structures, dynamics, physicochemical stability, functionality, etc. Firstly, brief histories of mAbs and neutron scattering, as well as some basic knowledges and principles of neutron scattering were introduced. Then specific examples were demonstrated. For the structure and structural evolution investigation of in dilute and concentrated mAbs solution, in situ small angle neutron scattering (SANS) was frequently utilized. Neutron reflectometry (NR) is powerful to probe the absorption behaviors of mAbs on various surfaces and interfaces. While for dynamic investigation, quasi-elastic scattering techniques such as neutron spin echo (NSE) demonstrate the capabilities. With this review, how to utilize and take advantages of neutron scattering on investigating structures and dynamics of mAbs were demonstrated and discussed.

Non-covalent Fc-Fab interactions significantly alter internal dynamics of an IgG1 antibody

The fragment-antigen-binding arms (Fab1 and Fab2) in a canonical immunoglobulin G (IgG) molecule have identical sequences and hence are always expected to exhibit symmetric conformations and dynamics. Using long all atom molecular simulations of a human IgG1 crystal structure 1HZH, we demonstrate that the translational and rotational dynamics of Fab1 and Fab2 also strongly depend on their interactions with each other and with the fragment-crystallizable (Fc) region. We show that the Fab2 arm in the 1HZH structure is non-covalently bound to the Fc region via long-lived hydrogen bonds, involving its light chain and both heavy chains of the Fc region. These highly stable interactions stabilize non-trivial conformer states with constrained fluctuations. We observe subtle modifications in Fab1 dynamics in response to Fab2-Fc interactions that points to novel allosteric interactions between the Fab arms. These results yield novel insights into the inter- and intra-fragment motions of immunoglobulins which could help us better understand the relation between their structure and function.

Modeling and Molecular Dynamic Simulation of F(ab')2 Fragment of Nimotuzumab for Lung Cancer Diagnostics

Lung cancer is one of the leading causes of cancer-related deaths in the world among both men and women. Several studies in the literature report that overexpression and mutation of the epidermal growth factor receptor (EGFR) are implicated in the pathogenesis of some lung cancers. Nimotuzumab is a humanized monoclonal antibody (mAb) that inhibits EGF binding because it binds to the extracellular domain of the EGFR. Nimotuzumab requires bivalent binding for stable attachment to cellular surface, which leads to nimotuzumab selectively binding to cells that express mAbs of moderate to high EGFR levels, and this could explain its low toxicity. This property has an advantage for development of nimotuzumab as a therapeutic and diagnostic agent. Monoclonal antibodies are large in size (150 kDa), thus penetrating slowly and residing in the blood for extended periods of time (from days to weeks); their use in imaging studies can result in low signal-to-background ratios and poor image quality. A reduction in the size of the immunoglobulin molecule has also been proposed as a means for increasing tumor penetration by mAbs. Nevertheless, it is known that the penetration of mAb into tumor cell is slow, due to its high molecular weight. Therefore, mAb is not very attractive to be used for imaging diagnostic purpose because of its kinetics and potential to elicit antibody response. The objective of this research was to study the homology modeling of a simpler functional molecule based on nimotuzumab, which consists of 2 antigen-binding fragments (Fab), namely, F(ab′)₂, using MODELER. The crystal structure of Fab of nimotuzumab from protein data bank was used as a template to construct the model of F(ab′)₂. Molecular dynamic simulation was performed to evaluate the stability of F(ab′)₂ and conformational changes of F(ab′)₂ in simulation. The result showed the dynamic behavior of antigen-binding site region of F(ab′)₂ throughout simulation. This result is expected to be useful in the further development of F(ab′)₂ fragment nimotuzumab as a lung cancer diagnostic.

A universal in silico V(D)J recombination strategy for developing humanized monoclonal antibodies

BACKGROUND

Humanization of mouse monoclonal antibodies (mAbs) is crucial for reducing their immunogenicity in humans. However, humanized mAbs often lose their binding affinities. Therefore, an in silico humanization method that can prevent the loss of the binding affinity of mAbs is needed.

METHODS

We developed an in silico V(D)J recombination platform in which we used V(D)J human germline gene sequences to design five humanized candidates of anti-tumor necrosis factor (TNF)-α mAbs (C1-C5) by using different human germline templates. The candidates were subjected to molecular dynamics simulation. In addition, the structural similarities of their complementarity-determining regions (CDRs) to those of original mouse mAbs were estimated to derive the weighted interatomic root mean squared deviation (wRMSDi) value. Subsequently, the correlation of the derived wRMSDi value with the half maximal effective concentration (EC50) and the binding affinity (KD) of the humanized anti-TNF-α candidates was examined. To confirm whether our in silico estimation method can be used for other humanized mAbs, we tested our method using the anti-epidermal growth factor receptor (EGFR) a4.6.1, anti-glypican-3 (GPC3) YP9.1 and anti-α4β1 integrin HP1/2L mAbs.

RESULTS

The R2 value for the correlation between the wRMSDi and log(EC50) of the recombinant Remicade and those of the humanized anti-TNF-α candidates was 0.901, and the R2 value for the correlation between wRMSDi and log(KD) was 0.9921. The results indicated that our in silico V(D)J recombination platform could predict the binding affinity of humanized candidates and successfully identify the high-affinity humanized anti-TNF-α antibody (Ab) C1 with a binding affinity similar to that of the parental chimeric mAb (5.13 × 10-10). For the anti-EGFR a4.6.1, anti-GPC3 YP9.1, and anti-α4β1 integrin HP1/2L mAbs, the wRMSDi and log(EC50) exhibited strong correlations (R2 = 0.9908, 0.9999, and 0.8907, respectively).

CONCLUSIONS

Our in silico V(D)J recombination platform can facilitate the development of humanized mAbs with low immunogenicity and high binding affinities. This platform can directly transform numerous mAbs with therapeutic potential to humanized or even human therapeutic Abs for clinical use.

IgG1 conformational behavior: elucidation of the N-glycosylation role via molecular dynamics

Currently, monoclonal antibodies (mAbs) are the most used biopharmaceuticals for human therapy. One of the key aspects in their development is the control of effector functions mediated by the interaction between fragment crystallizable (Fc) and Fcγ receptors, which is a secondary mechanism of the action of biotherapeutics. N-glycosylation at the Fc portion can regulate these mechanisms, and much experimental evidence suggests that modifications of glycosidic chains can affect antibody binding to FcγRIIIa, consequently impacting the immune response. In this work, we try to elucidate via in silico procedures the structural role exhibited by glycans, particularly fucose, in mAb conformational freedom that can potentially affect the receptor recognition. By using adalimumab, a marketed IgG1, as a general template, after rebuilding its three-dimensional (3D) structure through homology modeling approaches, we carried out molecular dynamics simulations of three differently glycosylated species: aglycosylated, afucosylated, and fucosylated antibody. Trajectory analysis showed different dynamical behaviors and pointed out that sugars can influence the overall 3D structure of the antibody. As a result, we propose a putative structural mechanism by which the presence of fucose introduces conformational constraints in the whole antibody and not only in the Fc domain, preventing a conformation suitable for the interaction with the receptor. As secondary evidence, we observed a high flexibility of the antibodies that is translated into an asymmetric behavior of Fab portions shown by all the simulated biopolymers, making the dynamical asymmetry a new, to our knowledge, molecular aspect that may be further investigated. In conclusion, these findings can help understand the contribution of sugars on the structural architecture of mAbs, paving the way to novel strategies of pharmaceutical development.

Communication pathways bridge local and global conformations in an IgG4 antibody

The affinity of an antibody for its antigen is primarily determined by the specific sequence and structural arrangement of the complementarity-determining regions (CDRs). Recent evidence, however, points toward a nontrivial relation between the CDR and distal sites: variations in the binding strengths have been observed upon mutating residues separated from the paratope by several nanometers, thus suggesting the existence of a communication network within antibodies, whose extension and relevance might be deeper than insofar expected. In this work, we test this hypothesis by means of molecular dynamics (MD) simulations of the IgG4 monoclonal antibody pembrolizumab, an approved drug that targets the programmed cell death protein 1 (PD-1). The molecule is simulated in both the apo and holo states, totalling 4 μs of MD trajectory. The analysis of these simulations shows that the bound antibody explores a restricted range of conformations with respect to the apo one, and that the global conformation of the molecule correlates with that of the CDR. These results support the hypothesis that pembrolizumab featues a multi-scale hierarchy of intertwined global and local conformational changes. The analysis pipeline developed in this work is general, and it can help shed further light on the mechanistic aspects of antibody function.

Stress-dependent flexibility of a full-length human monoclonal antibody: Insights from molecular dynamics to support biopharmaceutical development

After several decades of advancements in drug discovery, product development of biopharmaceuticals remains a time- and resource-consuming endeavor. One of the main reasons is associated to the lack of fundamental understanding of conformational dynamics of such biologic entities, and how they respond to various stresses encountered during manufacturing. In this work, we have studied the conformational dynamics of human IgG₁κ b12 monoclonal antibody (mAb) using molecular dynamics simulations. The hundreds of nanoseconds long trajectories reveal that b12 mAb is highly flexible. Its variable domains show greater conformational fluctuations than the constant domains. Additionally, it collapses towards a more globular shape in response to thermal stress, leading to decrease in the total solvent exposed surface area and radius of gyration. This behavior is more pronounced for the deglycosylated b12 mAb, and it appears to correlate with increase in inter-domain contacts between specific regions of the antibody. Conformational fluctuations also cause temporary formation and disruption of hydrophobic and charged patches on the antibody surface, which is particularly important for the prediction of CMC properties during development phases of antibody-based biotherapeutics. The insights gained through these simulations may help the development of biologic drugs, especially with regards to manufacturing processes where antibodies may undergo significant thermal stress.

Comparison of the pH- and thermally-induced fluctuations of a therapeutic antibody Fab fragment by molecular dynamics simulation

Successful development of protein therapeutics depends critically on achieving stability under a range of conditions. A deeper understanding of the drivers of instability across different stress conditions, will enable the engineering of more robust protein scaffolds. We compared the impacts of low pH and high temperature stresses on the structure of a humanized antibody fragment (Fab) A33, using atomistic molecular dynamics simulations, using a recent 2.5 Å crystal structure. This revealed that low-pH induced the loss of native contacts in the domain C_L. By contrast, thermal stress led to 5-7% loss of native contacts in all four domains, and simultaneous loss of >30% of native contacts in the V_L-V_H and C_L-C_H interfaces. This revealed divergent destabilising pathways under the two different stresses. The underlying cause of instability was probed using FoldX and Rosetta mutation analysis, and packing density calculations. These agreed that mutations in the C_L domain, and C_L-C_H1 interface have the greatest potential for stabilisation of Fab A33. Several key salt bridge losses underpinned the conformational change in C_L at low pH, whereas at high temperature, salt bridges became more dynamic, thus contributing to an overall destabilization. Lastly, the unfolding events at the two stress conditions exposed different predicted aggregation-prone regions (APR) to solvent, which would potentially lead to different aggregation mechanisms. Overall, our results identified the early stages of unfolding and stability-limiting regions of Fab A33, and the V_H and C_L domains as interesting future targets for engineering stability to both pH- and thermal-stresses simultaneously.

Antigen binding allosterically promotes Fc receptor recognition

A key question in immunology is whether antigen recognition and Fc receptor (FcR) binding are allosterically linked. This question is also relevant for therapeutic antibody design. Antibody Fab and Fc domains are connected by flexible unstructured hinge region. Fc chains have conserved glycosylation sites at Asn297, with each conjugated to a core heptasaccharide and forming biantennary Fc glycan. The glycans modulate the Fc conformations and functions. It is well known that the antibody Fab and Fc domains and glycan affect antibody activity, but whether these elements act independently or synergistically is still uncertain. We simulated four antibody complexes: free antibody, antigen-bound antibody, FcR-bound antibody, and an antigen-antibody-FcR complex. We found that, in the antibody's "T/Y" conformation, the glycans, and the Fc domain all respond to antigen binding, with the antibody population shifting to two dominant clusters, both with the Fc-receptor binding site open. The simulations reveal that the Fc-glycan-receptor complexes also segregate into two conformational clusters, one corresponding to the antigen-free antibody-FcR baseline binding, and the other with an antigen-enhanced antibody-FcR interaction. Our study confirmed allosteric communications in antibody-antigen recognition and following FcR activation. Even though we observed allosteric communications through the IgG domains, the most important mechanism that we observed is the communication via population shift, stimulated by antigen binding and propagating to influence FcR recognition.

Thinking outside the Laboratory: Analyses of Antibody Structure and Dynamics within Different Solvent Environments in Molecular Dynamics (MD) Simulations

Monoclonal antibodies (mAbs) have revolutionized the biomedical field, directly influencing therapeutics and diagnostics in the biopharmaceutical industry, while continuing advances in computational efficiency have enabled molecular dynamics (MD) simulations to provide atomistic insight into the structure and function of mAbs. Despite the success of MD tools, further optimizations are still required to enhance the computational efficiency of complex mAb simulations. This issue can be tackled by changing the way the solvent system is modelled to reduce the number of atoms to be tracked but must be done without compromising the accuracy of the simulations. In this work, the structure of the IgG_2a antibody was analyzed in three solvent systems: explicit water and ions, implicit water and ions, and implicit water and explicit ions. Root-mean-square distance (RMSD), root-mean-square fluctuations (RMSF), and interchain angles were used to quantify structural changes. The explicit system provides the most atomistic detail but is ~6 times slower in its exploration of configurational space and required ~4 times more computational time on our supercomputer than the implicit simulations. Overall, the behavior of the implicit and explicit simulations is quantifiably similar, with the inclusion of explicit ions in the implicit simulation stabilizing the antibody to reproduce well the statistical fluctuations of the fully explicit system. Therefore, this approach holds promise to maximize the use of computational resources to explore antibody behavior.

The molecular mechanism of two coreceptor binding site antibodies X5 and 17b neutralizing HIV-1: Insights from molecular dynamics simulation

The coreceptor binding site of gp120 plays an important role in HIV entry into host cell. X5 and 17b are typical coreceptor binding site antibodies with the ability to broadly neutralize HIV. Thus, here, to study the neutralizing mechanism of two antibodies and identify the source of two antibodies with different neutralizing ability, we performed molecular dynamics simulations for the complexes of X5 and 17b with gp120 and CD4. The simulation results indicate X5 and 17b mainly affects CD4 and coreceptor binding sites. Specifically, for CD4 binding site (CD4bs), the binding of antibodies has different effects on CD4bs with and without CD4. However, for coreceptor binding sites, the binding of the antibodies has consistent influence on the region adjacent to loop V3 despite of the simulated systems with or without CD4. The binding of the antibodies enhances the interactions of gp120 region adjacent to loop V3 with other region of gp120, which are unfavorable for conformational rearrangements of the region adjacent to loop V3 and further binding the coreceptor. Additionally, the interactions of loop V3 and bridging sheet with X5 lead to the close motion of loop V3 in X5 bound form, which further influences the rearrangements in gp120.

Fast antibody fragment motion: flexible linkers act as entropic spring

A flexible linker region between three fragments allows antibodies to adjust their binding sites to an antigen or receptor. Using Neutron Spin Echo Spectroscopy we observed fragment motion on a timescale of 7 ns with motional amplitudes of about 1 nm relative to each other. The mechanistic complexity of the linker region can be described by a spring model with Brownian motion of the fragments in a harmonic potential. Displacements, timescale, friction and force constant of the underlying dynamics are accessed. The force constant exhibits a similar strength to an entropic spring, with friction of the fragment matching the unbound state. The observed fast motions are fluctuations in pre-existing equilibrium configurations. The Brownian motion of domains in a harmonic potential is the appropriate model to examine functional hinge motions dependent on the structural topology and highlights the role of internal forces and friction to function.

Simulation and Theory of Antibody Binding to Crowded Antigen-Covered Surfaces

In this paper we introduce a fully flexible coarse-grained model of immunoglobulin G (IgG) antibodies parametrized directly on cryo-EM data and simulate the binding dynamics of many IgGs to antigens adsorbed on a surface at increasing densities. Moreover, we work out a theoretical model that allows to explain all the features observed in the simulations. Our combined computational and theoretical framework is in excellent agreement with surface-plasmon resonance data and allows us to establish a number of important results. (i) Internal flexibility is key to maximize bivalent binding, flexible IgGs being able to explore the surface with their second arm in search for an available hapten. This is made clear by the strongly reduced ability to bind with both arms displayed by artificial IgGs designed to rigidly keep a prescribed shape. (ii) The large size of IgGs is instrumental to keep neighboring molecules at a certain distance (surface repulsion), which essentially makes antigens within reach of the second Fab always unoccupied on average. (iii) One needs to account independently for the thermodynamic and geometric factors that regulate the binding equilibrium. The key geometrical parameters, besides excluded-volume repulsion, describe the screening of free haptens by neighboring bound antibodies. We prove that the thermodynamic parameters govern the low-antigen-concentration regime, while the surface screening and repulsion only affect the binding at high hapten densities. Importantly, we prove that screening effects are concealed in relative measures, such as the fraction of bivalently bound antibodies. Overall, our model provides a valuable, accurate theoretical paradigm beyond existing frameworks to interpret experimental profiles of antibodies binding to multi-valent surfaces of different sorts in many contexts.

Antibodies are the main working horses of the human immune system. Remarkably, no matter the size or the shape of the pathological intruders, these extremely flexible three-lobe molecules are able to form a complex, thus eliciting an immune response. What makes antibodies so effective? To answer this and other questions, we have developed a simplified computational scheme to simulate the dynamics of many antibodies interacting with each other and with antigens. Coarse-grained models are a great opportunity, as they give access to a true multi-scale approach to biologically relevant problems. In this work, our innovative method allowed us to simulate the binding process of many antibodies to surface-adsorbed antigens. This led us to elucidate and quantify many important physical aspects of their biological function in agreement with experiments, such as the role of their flexibility and crowding effects at the hapten-covered surface, which were shown to finely regulate the avidity.

Molecular basis of high viscosity in concentrated antibody solutions: Strategies for high concentration drug product development

Effective translation of breakthrough discoveries into innovative products in the clinic requires proactive mitigation or elimination of several drug development challenges. These challenges can vary depending upon the type of drug molecule. In the case of therapeutic antibody candidates, a commonly encountered challenge is high viscosity of the concentrated antibody solutions. Concentration-dependent viscosity behaviors of mAbs and other biologic entities may depend on pairwise and higher-order intermolecular interactions, non-native aggregation, and concentration-dependent fluctuations of various antibody regions. This article reviews our current understanding of molecular origins of viscosity behaviors of antibody solutions. We discuss general strategies and guidelines to select low viscosity candidates or optimize lead candidates for lower viscosity at early drug discovery stages. Moreover, strategies for formulation optimization and excipient design are also presented for candidates already in advanced product development stages. Potential future directions for research in this field are also explored.

Correlations between changes in conformational dynamics and physical stability in a mutant IgG1 mAb engineered for extended serum half-life

This study compares the local conformational dynamics and physical stability of an IgG1 mAb (mAb-A) with its corresponding YTE (M255Y/S257T/T259E) mutant (mAb-E), which was engineered for extended half-life in vivo. Structural dynamics was measured using hydrogen/deuterium (H/D) exchange mass spectrometry while protein stability was measured with differential scanning calorimetry (DSC) and size exclusion chromatography (SEC). The YTE mutation induced differences in H/D exchange kinetics at both pH 6.0 and 7.4. Segments covering the YTE mutation sites and the FcRn binding epitopes showed either subtle or no observable differences in local flexibility. Surprisingly, several adjacent segments in the CH2 and distant segments in the VH, CH1, and VL domains had significantly increased flexibility in the YTE mutant. Most notable among the observed differences is increased flexibility of the 244-254 segment of the CH2 domain, where increased flexibility has been shown previously to correlate with decreased conformational stability and increased aggregation propensity in other IgG1 mAbs (e.g., presence of destabilizing additives as well as upon de-glycosylation or methionine oxidation). DSC analysis showed decreases in both thermal onset (Tonset) and unfolding (Tm1) temperatures of 7°C and 6.7°C, respectively, for the CH2 domain of the YTE mutant. In addition, mAb-E aggregated faster than mAb-A under accelerated stability conditions as measured by SEC analysis. Hence, the relatively lower physical stability of the YTE mutant correlates with increased local flexibility of the 244-254 segment, providing a site-directed mutant example that this segment of the CH2 domain is an aggregation hot spot in IgG1 mAbs.

Fluorescence anisotropy (polarization): from drug screening to precision medicine

INTRODUCTION

Fluorescence anisotropy (FA) is one of the major established methods accepted by industry and regulatory agencies for understanding the mechanisms of drug action and selecting drug candidates utilizing a high-throughput format.

AREAS COVERED

This review covers the basics of FA and complementary methods, such as fluorescence lifetime anisotropy and their roles in the drug discovery process. The authors highlight the factors affecting FA readouts, fluorophore selection and instrumentation. Furthermore, the authors describe the recent development of a successful, commercially valuable FA assay for long QT syndrome drug toxicity to illustrate the role that FA can play in the early stages of drug discovery.

EXPERT OPINION

Despite the success in drug discovery, the FA-based technique experiences competitive pressure from other homogeneous assays. That being said, FA is an established yet rapidly developing technique, recognized by academic institutions, the pharmaceutical industry and regulatory agencies across the globe. The technical problems encountered in working with small molecules in homogeneous assays are largely solved, and new challenges come from more complex biological molecules and nanoparticles. With that, FA will remain one of the major work-horse techniques leading to precision (personalized) medicine.

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.

In silico selection of therapeutic antibodies for development: viscosity, clearance, and chemical stability

For mAbs to be viable therapeutics, they must be formulated to have low viscosity, be chemically stable, and have normal in vivo clearance rates. We explored these properties by observing correlations of up to 60 different antibodies of the IgG1 isotype. Unexpectedly, we observe significant correlations with simple physical properties obtainable from antibody sequences and by molecular dynamics simulations of individual antibody molecules. mAbs viscosities increase strongly with hydrophobicity and charge dipole distribution and decrease with net charge. Fast clearance correlates with high hydrophobicities of certain complementarity determining regions and with high positive or high negative net charge. Chemical degradation from tryptophan oxidation correlates with the average solvent exposure time of tryptophan residues. Aspartic acid isomerization rates can be predicted from solvent exposure and flexibility as determined by molecular dynamics simulations. These studies should aid in more rapid screening and selection of mAb candidates during early discovery.

Effects of galactosylation in immunoglobulin G from all-atom molecular dynamics simulations

Explicit water atomistic molecular dynamics simulations (200 ns, ∼330,000 atoms) were performed to study the effects of galactosylation in the Fc domain of immunoglobulin G1. Two glycoforms were simulated to observe changes in protein-carbohydrate interactions and carbohydrate structure. A high degree of flexibility was observed in the small hinge region of the protein, while large domains remained stable. The hinge region flexibility allowed both translation and rotation of the domains relative to each other, resulting in a large number of possible conformations available. The distributions of rotational orientations between the Fab1 and Fab2 domains showed that while these domains are able to orient themselves rather freely pointing in space they rotated in unison to remain rotationally oriented at specific angles. Additionally, removing specific terminal galactose residues increased the mobility of the carbohydrate, resulting in different protein-carbohydrate interactions. Glycosylation has been suggested as a route to improve the aggregation resistance of monoclonal antibodies for therapeutic treatments to aid the immune system. The results presented here may provide insight into the search for IgG molecules with increased aggregation resistance to be used as monoclonal antibodies.

Molecular dynamics simulation of the crystallizable fragment of IgG1-insights for the design of Fcabs

An interesting format in the development of therapeutic monoclonal antibodies uses the crystallizable fragment of IgG1 as starting scaffold. Engineering of its structural loops allows generation of an antigen binding site. However, this might impair the molecule’s conformational stability, which can be overcome by introducing stabilizing point mutations in the CH3 domains. These point mutations often affect the stability and unfolding behavior of both the CH2 and CH3 domains. In order to understand this cross-talk, molecular dynamics simulations of the domains of the Fc fragment of human IgG1 are reported. The structure of human IgG1-Fc obtained from X-ray crystallography is used as a starting point for simulations of the wild-type protein at two different pH values. The stabilizing effect of a single point mutation in the CH3 domain as well as the impact of the hinge region and the glycan tree structure connected to the CH2 domains is investigated. Regions of high local flexibility were identified as potential sites for engineering antigen binding sites. Obtained data are discussed with respect to the available X-ray structure of IgG1-Fc, directed evolution approaches that screen for stability and use of the scaffold IgG1-Fc in the design of antigen binding Fc proteins.