Estimation of the hydrophobic effect in an antigen-antibody protein-protein interface

Engineering antigen-presenting cells for immunotherapy of autoimmunity

Autoimmune diseases are burdensome conditions that affect a significant fraction of the global population. The hallmark of autoimmune disease is a host's immune system being licensed to attack its tissues based on specific antigens. There are no cures for autoimmune diseases. The current clinical standard for treating autoimmune diseases is the administration of immunosuppressants, which weaken the immune system and reduce auto-inflammatory responses. However, people living with autoimmune diseases are subject to toxicity, fail to mount a sufficient immune response to protect against pathogens, and are more likely to develop infections. Therefore, there is a concerted effort to develop more effective means of targeting immunomodulatory therapies to antigen-presenting cells, which are involved in modulating the immune responses to specific antigens. In this review, we highlight approaches that are currently in development to target antigen-presenting cells and improve therapeutic outcomes in autoimmune diseases.

Computational Mutagenesis of Antibody Fragments: Disentangling Side Chains from ΔΔG Predictions

The development of highly potent antibodies and antibody fragments as binding agents holds significant implications in fields such as biosensing and biotherapeutics. Their binding strength is intricately linked to the arrangement and composition of residues at the binding interface. Computational techniques offer a robust means to predict the three-dimensional structure of these complexes and to assess the affinity changes resulting from mutations. Given the interdependence of structure and affinity prediction, our objective here is to disentangle their roles. We aim to evaluate independently six side-chain reconstruction methods and ten binding affinity estimation techniques. This evaluation was pivotal in predicting affinity alterations due to single mutations, a key step in computational affinity maturation protocols. Our analysis focuses on a data set comprising 27 distinct antibody/hen egg white lysozyme complexes, each with crystal structures and experimentally determined binding affinities. Using six different side-chain reconstruction methods, we transformed each structure into its corresponding mutant via in silico single-point mutations. Subsequently, these structures undergo minimization and molecular dynamics simulation. We therefore estimate ΔΔG values based on the original crystal structure, its energy-minimized form, and the ensuing molecular dynamics trajectories. Our research underscores the critical importance of selecting reliable side-chain reconstruction methods and conducting thorough molecular dynamics simulations to accurately predict the impact of mutations. In summary, our study demonstrates that the integration of conformational sampling and scoring is a potent approach to precisely characterizing mutation processes in single-point mutagenesis protocols and crucial for computational antibody design.

Contribution of Signaling Partner Association to Strigolactone Receptor Selectivity

The parasitic plant witchweed, Striga hermonthica, results in agricultural losses of billions of dollars per year. It perceives its host via plant hormones called strigolactones, which act as germination stimulants for witchweed. Strigolactone signaling involves substrate binding to the strigolactone receptor, followed by substrate hydrolysis and a conformational change from an inactive, or open state, to an active, or closed state. In the active state, the receptor associates with a signaling partner, MAX2. Recently, it was shown that this MAX2 association process acts as a strong contributor to the uniquely high signaling activity observed in ShHTL7; however, it is unknown why ShHTL7 has enhanced MAX2 association affinity. Using an umbrella sampling molecular dynamics approach, we characterized the association processes of AtD14, ShHTL7, a mutant of ShHTL7, and ShHTL6 with MAX2 homologue OsD3. From these results, we show that ShHTL7 has an enhanced standard binding free energy of OsD3 compared to those of the other receptors. Additionally, our results suggest that the overall topology of the T2/T3 helix region is likely an important modulator of MAX2 binding. Thus, differences in MAX2 association, modulated by differences in the T2/T3 helix region, are a contributor to differences in signaling activity between different strigolactone receptors.

Current advances in biopharmaceutical informatics: guidelines, impact and challenges in the computational developability assessment of antibody therapeutics

Therapeutic monoclonal antibodies and their derivatives are key components of clinical pipelines in the global biopharmaceutical industry. The availability of large datasets of antibody sequences, structures, and biophysical properties is increasingly enabling the development of predictive models and computational tools for the "developability assessment" of antibody drug candidates. Here, we provide an overview of the antibody informatics tools applicable to the prediction of developability issues such as stability, aggregation, immunogenicity, and chemical degradation. We further evaluate the opportunities and challenges of using biopharmaceutical informatics for drug discovery and optimization. Finally, we discuss the potential of developability guidelines based on in silico metrics that can be used for the assessment of antibody stability and manufacturability.

Structural dissimilarity from self drives neoepitope escape from immune tolerance

T-cell recognition of peptides incorporating nonsynonymous mutations, or neoepitopes, is a cornerstone of tumor immunity and forms the basis of new immunotherapy approaches including personalized cancer vaccines. Yet as they are derived from self-peptides, the means through which immunogenic neoepitopes overcome immune self-tolerance are often unclear. Here we show that a point mutation in a non-major histocompatibility complex anchor position induces structural and dynamic changes in an immunologically active ovarian cancer neoepitope. The changes pre-organize the peptide into a conformation optimal for recognition by a neoepitope-specific T-cell receptor, allowing the receptor to bind the neoepitope with high affinity and deliver potent T-cell signals. Our results emphasize the importance of structural and physical changes relative to self in neoepitope immunogenicity. Considered broadly, these findings can help explain some of the difficulties in identifying immunogenic neoepitopes from sequence alone and provide guidance for developing novel, neoepitope-based personalized therapies.

Conformational Plasticity in Broadly Neutralizing HIV-1 Antibodies Triggers Polyreactivity

Human high-affinity antibodies to pathogens often recognize unrelated ligands. The molecular origin and the role of this polyreactivity are largely unknown. Here, we report that HIV-1 broadly neutralizing antibodies (bNAbs) are frequently polyreactive, cross-reacting with non-HIV-1 molecules, including self-antigens. Mutating bNAb genes to increase HIV-1 binding and neutralization also results in de novo polyreactivity. Unliganded paratopes of polyreactive bNAbs with improved HIV-1 neutralization exhibit a conformational flexibility, which contributes to enhanced affinity of bNAbs to various HIV-1 envelope glycoproteins and non-HIV antigens. Binding adaptation of polyreactive bNAbs to the divergent ligands mainly involves hydrophophic interactions. Plasticity of bNAbs' paratopes may, therefore, facilitate accommodating divergent viral variants, but it simultaneously triggers promiscuous binding to non-HIV-1 antigens. Thus, a certain level of polyreactivity can be a mark of adaptable antibodies displaying optimal pathogens' recognition.

Assessment of Solvated Interaction Energy Function for Ranking Antibody-Antigen Binding Affinities

Affinity modulation of antibodies and antibody fragments of therapeutic value is often required in order to improve their clinical efficacies. Virtual affinity maturation has the potential to quickly focus on the critical hotspot residues without the combinatorial explosion problem of conventional display and library approaches. However, this requires a binding affinity scoring function that is capable of ranking single-point mutations of a starting antibody. We focus here on assessing the solvated interaction energy (SIE) function that was originally developed for and is widely applied to scoring of protein-ligand binding affinities. To this end, we assembled a structure-function data set called Single-Point Mutant Antibody Binding (SiPMAB) comprising several antibody-antigen systems suitable for this assessment, i.e., based on high-resolution crystal structures for the parent antibodies and coupled with high-quality binding affinity measurements for sets of single-point antibody mutants in each system. Using this data set, we tested the SIE function with several mutation protocols based on the popular methods SCWRL, Rosetta, and FoldX. We found that the SIE function coupled with a protocol limited to sampling only the mutated side chain can reasonably predict relative binding affinities with a Spearman rank-order correlation coefficient of about 0.6, outperforming more aggressive sampling protocols. Importantly, this performance is maintained for each of the seven system-specific component subsets as well as for other relevant subsets including non-alanine and charge-altering mutations. The transferability and enrichment in affinity-improving mutants can be further enhanced using consensus ranking over multiple methods, including the SIE, Talaris, and FOLDEF energy functions. The knowledge gained from this study can lead to successful prospective applications of virtual affinity maturation.

Current Experimental Methods for Characterizing Protein-Protein Interactions

Protein molecules often interact with other partner protein molecules in order to execute their vital functions in living organisms. Characterization of protein-protein interactions thus plays a central role in understanding the molecular mechanism of relevant protein molecules, elucidating the cellular processes and pathways relevant to health or disease for drug discovery, and charting large-scale interaction networks in systems biology research. A whole spectrum of methods, based on biophysical, biochemical, or genetic principles, have been developed to detect the time, space, and functional relevance of protein-protein interactions at various degrees of affinity and specificity. This article presents an overview of these experimental methods, outlining the principles, strengths and limitations, and recent developments of each type of method.

Current Experimental Methods for Characterizing Protein-Protein Interactions

Protein molecules often interact with other partner protein molecules in order to execute their vital functions in living organisms. Characterization of protein-protein interactions thus plays a central role in understanding the molecular mechanism of relevant protein molecules, elucidating the cellular processes and pathways relevant to health or disease for drug discovery, and charting large-scale interaction networks in systems biology research. A whole spectrum of methods, based on biophysical, biochemical, or genetic principles, have been developed to detect the time, space, and functional relevance of protein-protein interactions at various degrees of affinity and specificity. This article presents an overview of these experimental methods, outlining the principles, strengths and limitations, and recent developments of each type of method.

Molecular basis of in vitro affinity maturation and functional evolution of a neutralizing anti-human GM-CSF antibody

X-ray structure analysis of 4 antibody Fab fragments, each in complex with human granulocyte macrophage colony stimulating factor (GM-CSF), was performed to investigate the changes at the protein-protein binding interface during the course of in vitro affinity maturation by phage display selection. The parental antibody MOR03929 was compared to its derivatives MOR04252 (CDR-H2 optimized), MOR04302 (CDR-L3 optimized) and MOR04357 (CDR-H2 and CDR-L3 optimized). All antibodies bind to a conformational epitope that can be divided into 3 sub-epitopes. Specifically, MOR04357 binds to a region close to the GM-CSF N-terminus (residues 11-24), a short second sub-epitope (residues 83-89) and a third at the C-terminus (residues 112-123). Modifications introduced during affinity maturation in CDR-H2 and CDR-L3 led to the establishment of additional hydrogen bonds and van der Waals contacts, respectively, providing a rationale for the observed improvement in binding affinity and neutralization potency. Once GM-CSF is complexed to the antibodies, modeling predicts a sterical clash with GM-CSF binding to GM-CSF receptor α and β chain. This predicted mutually exclusive binding was confirmed by a GM-CSF receptor α chain ligand binding inhibition assay. Finally, high throughput sequencing of clones obtained after affinity maturation phage display pannings revealed highly selected consensus sequences for CDR-H2 as well for CDR-L3, which are in accordance with the sequence of the highest affinity antibody MOR04357. The resolved crystal structures highlight the criticality of these strongly selected residues for high affinity interaction with GM-CSF.

Biosensor Regeneration: A Review of Common Techniques and Outcomes

Biosensors are ideally portable, low-cost tools for the rapid detection of pathogens, proteins, and other analytes. The global biosensor market is currently worth over 10 billion dollars annually and is a burgeoning field of interdisciplinary research that is hailed as a potential revolution in consumer, healthcare, and industrial testing. A key barrier to the widespread adoption of biosensors, however, is their cost. Although many systems have been validated in the laboratory setting and biosensors for a range of analytes are proven at the concept level, many have yet to make a strong commercial case for their acceptance. Though it is true with the development of cheaper electrodes, circuits, and components that there is a downward pressure on costs, there is also an emerging trend toward the development of multianalyte biosensors that is pushing in the other direction. One way to reduce the cost that is suitable for certain systems is to enable their reuse, thus reducing the cost per test. Regenerating biosensors is a technique that can often be used in conjunction with existing systems in order to reduce costs and accelerate the commercialization process. This article discusses the merits and drawbacks of regeneration schemes that have been proven in various biosensor systems and indicates parameters for successful regeneration based on a systematic review of the literature. It also outlines some of the difficulties encountered when considering the role of regeneration at the point of use. A brief meta-analysis has been included in this review to develop a working definition for biosensor regeneration, and using this analysis only ∼60% of the reported studies analyzed were deemed a success. This highlights the variation within the field and the need to normalize regeneration as a standard process across the field by establishing a consensus term.

Simulation of protein association: Kinetic pathways towards crystal contacts

We conducted molecular dynamics simulations combined with distance-based umbrella sampling and forward flux sampling to investigate the early stages of protein crystallization. Formation of contacts with long-range interactions and/or an exposed position on the protein surface was kinetically preferred over more stable hydrophobic contacts with a shorter attractive range, while the thermodynamic stability of the protein crystal was provided by hydrophobic interactions. Contacts with a large interaction area showed complex dissociation pathways that were not detected by distance-based umbrella sampling. Instead, forward flux sampling simulations of contact dissociation identified long-range attractive interactions.

Inferring the microscopic surface energy of protein-protein interfaces from mutation data

Mutations at protein-protein recognition sites alter binding strength by altering the chemical nature of the interacting surfaces. We present a simple surface energy model, parameterized with empirical ΔΔG values, yielding mean energies of -48 cal mol(-1) Å(-2) for interactions between hydrophobic surfaces, -51 to -80 cal mol(-1) Å(-2) for surfaces of complementary charge, and 66-83 cal mol(-1) Å(-2) for electrostatically repelling surfaces, relative to the aqueous phase. This places the mean energy of hydrophobic surface burial at -24 cal mol(-1) Å(-2) . Despite neglecting configurational entropy and intramolecular changes, the model correlates with empirical binding free energies of a functionally diverse set of rigid-body interactions (r = 0.66). When used to rerank docking poses, it can place near-native solutions in the top 10 for 37% of the complexes evaluated, and 82% in the top 100. The method shows that hydrophobic burial is the driving force for protein association, accounting for 50-95% of the cohesive energy. The model is available open-source from http://life.bsc.es/pid/web/surface_energy/ and via the CCharpPPI web server http://life.bsc.es/pid/ccharppi/.

Antibodies: Computer-Aided Prediction of Structure and Design of Function

With the advent of high-throughput sequencing, and the increased availability of experimental structures of antibodies and antibody-antigen complexes, comes the improvement of computational approaches to predict the structure and design the function of antibodies and antibody-antigen complexes. While antibodies pose formidable challenges for protein structure prediction and design due to their large size and highly flexible loops in the complementarity-determining regions, they also offer exciting opportunities: the central importance of antibodies for human health results in a wealth of structural and sequence information that—as a knowledge base—can drive the modeling algorithms by limiting the conformational and sequence search space to likely regions of success. Further, efficient experimental platforms exist to test predicted antibody structure or designed antibody function, thereby leading to an iterative feedback loop between computation and experiment. We briefly review the history of computer-aided prediction of structure and design of function in the antibody field before we focus on recent methodological developments and the most exciting application examples.

Aggregation-resistant domain antibodies engineered with charged mutations near the edges of the complementarity-determining regions

Antibodies commonly contain hydrophobic residues within their complementarity-determining regions (CDRs) that mediate binding to target antigens. Unfortunately, hydrophobic CDRs can also promote antibody aggregation, which is especially concerning for therapeutic antibodies due to the immunogenicity of antibody aggregates. Here we investigate how the sequences of CDRs within single-domain (V(H)) antibodies specific for the Alzheimer's amyloid β peptide can be engineered to resist aggregation without reducing binding affinity. We find that domain antibodies containing clusters of hydrophobic residues within their third CDR (CDR3) are prone to aggregate within days at 25°C and minutes above 70°C. However, inserting two or more negatively charged residues at each edge of CDR3 potently suppresses antibody aggregation without altering binding affinity. We also find that inserting charged mutations at one edge of CDR3 (N- or C-terminal) prevents aggregation, but only if such mutations are located at the edge closest to most hydrophobic portion of CDR3. In contrast, charged mutations outside of CDR3 fail to suppress aggregation. Our findings demonstrate that the sequence of CDR loops can be engineered in a systematic manner to improve antibody solubility without altering binding affinity or specificity.

Subtype- and antigenic site-specific differences in biophysical influences on evolution of influenza virus hemagglutinin

Background

Influenza virus undergoes rapid evolution by both antigenic shift and antigenic drift. Antibodies, particularly those binding near the receptor-binding site of hemagglutinin (HA) or the neuraminidase (NA) active site, are thought to be the primary defense against influenza infection, and mutations in antibody binding sites can reduce or eliminate antibody binding. The binding of antibodies to their cognate antigens is governed by such biophysical properties of the interacting surfaces as shape, non-polar and polar surface area, and charge.

Methods

To understand forces shaping evolution of influenza virus, we have examined HA sequences of human influenza A and B viruses, assigning each amino acid values reflecting total accessible surface area, non-polar and polar surface area, and net charge due to the side chain. Changes in each of these values between neighboring sequences were calculated for each residue and mapped onto the crystal structures.

Results

Areas of HA showing the highest frequency of pairwise changes agreed well with previously identified antigenic sites in H3 and H1 HAs, and allowed us to propose more detailed antigenic maps and novel antigenic sites for H1 and influenza B HA. Changes in biophysical properties differed between HAs of different subtypes, and between different antigenic sites of the same HA. For H1, statistically significant differences in several biophysical quantities compared to residues lying outside antigenic sites were seen for some antigenic sites but not others. Influenza B antigenic sites all show statistically significant differences in biophysical quantities for all antigenic sites, whereas no statistically significant differences in biophysical quantities were seen for any antigenic site is seen for H3. In many cases, residues previously shown to be under positive selection at the genetic level also undergo rapid change in biophysical properties.

Conclusions

The biophysical consequences of amino acid changes introduced by antigenic drift vary from subtype to subtype, and between different antigenic sites. This suggests that the significance of antibody binding in selecting new variants may also be variable for different antigenic sites and influenza subtypes.

Generation of recombinant antibodies and means for increasing their affinity

Highly specific interaction with foreign molecules is a unique feature of antibodies. Since 1975, when Keller and Milstein proposed the method of hybridoma technology and prepared mouse monoclonal antibodies, many antibodies specific to various antigens have been obtained. Recent development of methods for preparation of recombinant DNA libraries and in silico bioinformatics approaches for protein structure analysis makes possible antibody preparation using gene engineering approaches. The development of gene engineering methods allowed creating recombinant antibodies and improving characteristics of existing antibodies; this significantly extends the applicability of antibodies. By modifying biochemical and immunochemical properties of antibodies by changing their amino acid sequences it is possible to create antibodies with properties optimal for certain tasks. For example, application of recombinant technologies resulted in antibody preparation of high affinity significantly exceeding the initial affinity of natural antibodies. In this review we summarize information about the structure, modes of preparation, and application of recombinant antibodies and their fragments and also consider the main approaches used to increase antibody affinity.

High-resolution, three-dimensional modeling of human leukocyte antigen class I structure and surface electrostatic potential reveals the molecular basis for alloantibody binding epitopes

The potential of human leukocyte antigens (HLA) to stimulate humoral alloimmunity depends on the orientation, accessibility and physiochemical properties of polymorphic amino acids. We have generated high-resolution structural and physiochemical models of all common HLA class I alleles and analyzed the impact of amino acid polymorphisms on surface electrostatic potential. Atomic resolution three-dimensional structural models of HLA class I molecules were generated using the MODELLER computer algorithm. The molecular surface electrostatic potential was calculated using the DelPhi program. To confirm that electrostatic surface topography reflects known HLA B cell epitopes, we examined Bw4 and Bw6 and ascertained the impact of amino acid polymorphisms on their tertiary and physiochemical composition. The HLA protein structures generated performed well when subjected to stereochemical and energy-based testing for structural integrity. The electrostatic pattern and conformation of Bw4 and Bw6 epitopes are maintained among HLA molecules even when expressed in a different structural context. Importantly, variation in epitope amino acid composition does not always translate into a different electrostatic motif, providing an explanation for serologic cross-reactivity. Mutations of critical amino acids that abrogate antibody binding also induce distinct changes in epitope electrostatic properties. In conclusion, high-resolution structural modeling provides a physiochemical explanation for serologic patterns of antibody binding and provides novel insights into HLA immunogenicity.

Peptides and proteins as a continuing exciting source of inspiration for peptidomimetics

Despite their enormous diversity in biological function and structure, peptides and proteins are endowed with properties that have induced and stimulated the development of peptidomimetics. Clearly, peptides can be considered as the "stem" of a phylogenetic molecular development tree from which branches of oligomeric peptidomimetics such as peptoids, peptidosulfonamides, urea peptidomimetics, as well as β-peptides have sprouted. It is still a challenge to efficiently synthesize these oligomeric species, and study their structural and biological properties. Combining peptides and peptidomimetics led to the emergence of peptide-peptidomimetic hybrids in which one or more (proteinogenic) amino acid residues have been replaced with these mimetic residues. In scan-like approaches, the influence of these replacements on biological activity can then be studied, to evaluate to what extent a peptide can be transformed into a peptidomimetic structure while maintaining, or even improving, its biological properties. A central issue, especially with the smaller peptides, is the lack of secondary structure. Important approaches to control secondary structure include the introduction of α,α-disubstituted amino acids, or (di)peptidomimetic structures such as the Freidinger lactam. Apart from intra-amino acid constraints, inter-amino acid constraints for formation of a diversity of cyclic peptides have shaped a thick branch. Apart from the classical disulfide bridges, the repertoire has been extended to include sulfide and triazole bridges as well as the single-, double- and even triple-bond replacements, accessible by the extremely versatile ring-closing alkene/alkyne metathesis approaches. The latter approach is now the method of choice for the secondary structure that presents the greatest challenge for structural stabilization: the α-helix.

Assessing energetic contributions to binding from a disordered region in a protein-protein interaction

Many functional proteins are at least partially disordered prior to binding. Although the structural transitions upon binding of disordered protein regions can influence the affinity and specificity of protein complexes, their precise energetic contributions to binding are unknown. Here, we use a model protein-protein interaction system in which a locally disordered region has been modified by directed evolution to quantitatively assess the thermodynamic and structural contributions to binding of disorder-to-order transitions. Through X-ray structure determination of the protein binding partners before and after complex formation and isothermal titration calorimetry of the interactions, we observe a correlation between protein ordering and binding affinity for complexes along this affinity maturation pathway. Additionally, we show that discrepancies between observed and calculated heat capacities based on buried surface area changes in the protein complexes can be explained largely by heat capacity changes that would result solely from folding the locally disordered region. Previously developed algorithms for predicting binding energies of protein-protein interactions, however, are unable to correctly model the energetic contributions of the structural transitions in our model system. While this highlights the shortcomings of current computational methods in modeling conformational flexibility, it suggests that the experimental methods used here could provide training sets of molecular interactions for improving these algorithms and further rationalizing molecular recognition in protein-protein interactions.

Systematic mutation and thermodynamic analysis of central tyrosine pairs in polyspecific NKG2D receptor interactions

The homodimeric, activating natural killer cell receptor NKG2D interacts with multiple monomeric ligands polyspecifically, yet without central conformational flexibility. Crystal structures of multiple NKG2D-ligand interactions have identified the NKG2D tyrosine pair Tyr 152 and Tyr 199 as forming multiple specific but diverse interactions with MICA and related proteins. Here we systematically altered each tyrosine to tryptophan, phenylalanine, isoleucine, leucine, valine, serine, and alanine to measure the effect of mutation on affinity and thermodynamics for binding a range of similar ligands: MICA, the higher-affinity ligand MICB, and MICdesign, a high-affinity version of MICA that shares all NKG2D contact residues with MICA. Affinity and residue size were related: tryptophan could often substitute for tyrosine without loss of affinity; loss of the tyrosine hydroxyl through mutation to phenylalanine was tolerated more at position 152 than 199; and the smallest residues coincide with lowest affinities in general. NKG2D mutant van't Hoff binding thermodynamics generally show that substitution of other residues for tyrosine causes a moderate positive or flat van't Hoff slope consistent with moderate loss of binding enthalpy. One set of NKG2D mutations caused MICA to adopt a positive van't Hoff slope corresponding to absorption of heat, and another set caused MICB to adopt a negative slope of greater heat release than wild-type. MICdesign shared one example of the first set with MICA and one of the second set with MICB. When the NKG2D mutation affinities were arranged according to change in nonpolar surface area and compared to results from specific antibody-antigen and protein-peptide interactions, it was found that hydrophobic surface loss in NKG2D reduced binding affinity less than reported in the other contexts. The hydrophobic effect at the center of the NKG2D binding appears more similar to that at the periphery of an antibody-antigen binding site than at its center. Therefore the polyspecific NKG2D binding site is more tolerant of structural alteration in general than either an antibody-antigen or protein-peptide binding site, and this tolerance may adapt NKG2D to a broad range of protein surfaces with micromolar affinity.

Stacking and energetic contribution of aromatic islands at the binding interface of antibody proteins

Background

The enrichment and importance of some aromatic residues, such as Tyr and Trp, have been widely noticed at the binding interfaces of antibodies from many experimental and statistical results, some of which were even identified as “hot spots” contributing significantly greater to the binding affinity than other amino acids. However, how these aromatic residues influence the immune binding still deserves further investigation. A large-scale examination was done regarding the local spatial environment around the interfacial Tyr or Trp residues. Energetic contribution of these Tyr and Trp residues to the binding affinity was then studied regarding 82 representative antibody interfaces covering 509 immune complexes from the PDB database and IMGT/3Dstructure-DB.

Results

The connectivity analysis of interfacial residues showed that Tyr and Trp tended to cluster into the spatial Aromatic Islands (AI) rather than being distributed randomly at the antibody interfaces. Out of 82 antibody-antigen complexes, 72% (59) interfaces were found to contain AI with more than 3 aromatic residues. The statistical test against an empirical distribution indicated that the existence of AI was significant in about 60% representative antibody interfaces. Secondly, the loss of solvent accessible surface area (SASA) for side chains of aromatic residues between actually crowded state and independent state was nicely correlated with the AI size increasing in a linearly positive way which indicated that the aromatic side chains in AI tended to take a compact and ordered stacking conformation at the interfaces. Interestingly, the SASA loss of AI was also correlated roughly with the averaged gap of binding free energy between the theoretical and experimental data for immune complexes.

Conclusions

The results of our study revealed the wide existence and statistical significance of “Aromatic Island” (AI) composed of the spatially clustered Tyr and Trp residues at the antibody interfaces. The regular arrangement and stacking of aromatic side chains in AI could probably produce extra cooperative effects to the binding affinity which was firstly observed through the large-scale data analysis. The finding in this work not only provides insights into the functional role of aromatic residues in the antibody-antigen interaction, but also may facilitate the antibody engineering and potential clinical applications.

Naturally occurring variability in the envelope glycoprotein of HIV-1 and development of cell entry inhibitors

Naturally occurring genetic variability across HIV-1 subtypes causes amino acid polymorphisms in encoded HIV-1 proteins including the envelope glycoproteins associated with viral entry. The effects of amino acid polymorphisms on the mechanism of HIV-1 entry into cells, a process initiated by the binding of the viral envelope glycoprotein gp120 to the cellular CD4 receptor, are largely unknown. In this study, we demonstrate that amino acid polymorphisms affect the structural stability and domain cooperativity of gp120 and that those differences are reflected in the binding mechanism of the viral envelope glycoprotein to the cell surface receptor and coreceptor. Moreover, subtype differences also affect the binding behavior of experimental HIV cell entry inhibitors. While gp120-A has a slightly lower denaturation temperature than gp120-B, the most notable stability difference is that for gp120-B the van't Hoff to calorimetric enthalpy ratio (DeltaH(vH)/DeltaH) is 0.95 whereas for gp120-A is 0.6, indicative of more cooperative domain/domain interactions in gp120-B, as this protein more closely approaches a two-state transition. Isothermal titration calorimetry demonstrates that CD4 and 17b (a surrogate antibody for the chemokine coreceptor) exhibit 7- and 3-fold weaker binding affinities for gp120-A. The binding of these proteins as well as that of the experimental entry inhibitor NBD-556 induces smaller conformational changes in gp120-A as evidenced by significantly smaller binding enthalpies and binding entropies. Together, these results describe the effects of gp120 polymorphisms on binding to host cell receptors and emphasize that guidelines for developing future entry inhibitors must recognize and deal with genomic differences between HIV strains.

Contribution of asparagine residues to the stabilization of a proteinaceous antigen-antibody complex, HyHEL-10-hen egg white lysozyme

Many germ line antibodies have asparagine residues at specific sites to achieve specific antigen recognition. To study the role of asparagine residues in the stabilization of antigen-antibody complexes, we examined the interaction between hen egg white lysozyme (HEL) and the corresponding HyHEL-10 variable domain fragment (Fv). We introduced Ala and Asp substitutions into the Fv side chains of L-Asn-31, L-Asn-32, and L-Asn-92, which interact directly with residues in HEL via hydrogen bonding in the wild-type Fv-HEL complex, and we investigated the interactions between these mutant antibodies and HEL. Isothermal titration calorimetric analysis showed that all the mutations decreased the negative enthalpy change and decreased the association constants of the interaction. Structural analyses showed that the effects of the mutations on the structure of the complex could be compensated for by conformational changes and/or by gains in other interactions. Consequently, the contribution of two hydrogen bonds was minor, and their abolition by mutation resulted in only a slight decrease in the affinity of the antibody for its antigen. By comparison, the other two hydrogen bonds buried at the interfacial area had large enthalpic advantage, despite entropic loss that was perhaps due to stiffening of the interface by the bonds, and were crucial to the strength of the interaction. Deletion of these strong hydrogen bonds could not be compensated for by other structural changes. Our results suggest that asparagine can provide the two functional groups for strong hydrogen bond formation, and their contribution to the antigen-antibody interaction can be attributed to their limited flexibility and accessibility at the complex interface.

Binding site on the transferrin receptor for the parvovirus capsid and effects of altered affinity on cell uptake and infection

Canine parvovirus (CPV) and its relative feline panleukopenia virus (FPV) bind the transferrin receptor type 1 (TfR) to infect their host cells but show differences in the interactions with the feline and canine TfRs that determine viral host range and tissue tropism. We changed apical and protease-like domain residues by introducing point mutations and adding or removing glycosylation signals, and we then examined the interactions of those mutant TfRs with the capsids. Most substitutions had little effect on virus binding and uptake. However, mutations of several sites in the apical domain of the receptor either prevented binding to the capsids or reduced the affinity of receptor binding to various degrees. Glycans within the virus binding face of the apical domain also controlled capsid binding. CPV, but not the related feline parvovirus, could use receptors containing a canine TfR-specific glycosylation to mediate efficient infection, while addition of other N-linked glycosylation sites into the virus binding face of the feline apical domain reduced or eliminated both binding and infection. Replacement of critical feline TfR residue 221 with every amino acid had effects on binding and infection which were significantly associated with the biochemical properties of the residue replaced. Receptors with reduced affinities mostly showed proportional changes in their ability to mediate infection. Testing feline TfR variants for their binding and uptake patterns in cells showed that low-affinity versions bound fewer capsids and also differed in attachment to the cell surface and filopodia, but transport to the perinuclear endosome was similar.

Predicting HLA class I alloantigen immunogenicity from the number and physiochemical properties of amino acid polymorphisms

BACKGROUND

Knowledge of the human leukocyte antigen (HLA) amino acid (AA) sequence combined with crystallographic structural data may enable prediction of the relative immunogenicity of individual donor/recipient HLA mismatches.

METHODS

Multiple sera from 32 highly sensitized patients awaiting kidney transplantation were screened using Luminex/single-antigen beads to determine the HLA-specific antibody levels against mismatched HLA class I specificities. A computer program was developed to allow intralocus and interlocus comparison of mismatched HLA-A and -B specificities with corresponding recipient HLA class I type, and to determine the number, position, and physiochemical disparity (hydrophobicity and electrostatic charge) of polymorphic AA.

RESULTS

HLA-specific antibody was detected against 1666 (85%) of the 1964 mismatched HLA specificities evaluated, with a close correlation between increasing number of AA polymorphisms and the presence and magnitude of the alloantibody response (P<0.0001). Hydrophobicity and electrostatic charge disparity scores were independent predictors of alloantibody production (adjusted P=0.0009 and P=0.0005, respectively). Mismatched specificities with physiochemical scores within the first decile of the scale led to weak alloantibody responses (median fluorescence intensity 2330), whereas those with scores above the sixth decile led to strong alloantibody production (median fluorescence intensity >10,000).

CONCLUSION

Differences in AA number, hydrophobicity, and electrostatic charge between HLA class I specificities enable prediction of donor HLA class I types with low immunogenicity for a given recipient.

Association energetics of cross-reactive and specific antibodies

HyHEL-8, HyHEL-10, and HyHEL-26 (HH8, HH10, and HH26, respectively) are murine monoclonal IgG(1) antibodies which share over 90% variable-region amino acid sequence identity and recognize identical structurally characterized epitopes on hen egg white lysozyme (HEL). Previous immunochemical and surface plasmon resonance-based studies have shown that these antibodies differ widely in their tolerance of mutations in the epitope. While HH8 is the most cross-reactive, HH26 is rigidified by a more extensive network of intramolecular salt links and is highly specific, with both association and dissociation rates strongly affected by epitope mutations. HH10 is of intermediate specificity, and epitope mutations produce changes primarily in the dissociation rate. Calorimetric characterization of the association energetics of these three antibodies with the native antigen HEL and with Japanese quail egg white lysozyme (JQL), a naturally occurring avian variant, shows that the energetics of interaction correlate with cross-reactivity and specificity. These results suggest that the greater cross-reactivity of HH8 may be mediated by a combination of conformational flexibility and less specific intermolecular interactions. Thermodynamic calculations suggest that upon association HH8 incurs the largest configurational entropic penalty and also the smallest loss of enthalpic driving force with variant antigen. Much smaller structural perturbations are expected in the formation of the less flexible HH26 complex, and the large loss of enthalpic driving force observed with variant antigen reflects its specificity. The observed thermodynamic parameters correlate well with the observed functional behavior of the antibodies and illustrate fundamental differences in thermodynamic characteristics between cross-reactive and specific molecular recognition.

Different thermodynamic binding mechanisms and peptide fine specificities associated with a panel of structurally similar high-affinity T cell receptors

To understand the mechanisms that govern T cell receptor (TCR)-peptide MHC (pMHC) binding and the role that different regions of the TCR play in affinity and antigen specificity, we have studied the TCR from T cell clone 2C. High-affinity mutants of the 2C TCR that bind QL9-L(d) as a strong agonist were generated previously by site-directed mutagenesis of complementarity determining regions (CDRs) 1beta, 2alpha, 3alpha, or 3beta. We performed isothermal titration calorimetry to assess whether they use similar thermodynamic mechanisms to achieve high affinity for QL9-L(d). Four of the five TCRs examined bound to QL9-L(d) in an enthalpically driven, entropically unfavorable manner. In contrast, the high-affinity CDR1beta mutant resembled the wild-type 2C TCR interaction, with favorable entropy. To assess fine specificity, we measured the binding and kinetics of these mutants for both QL9-L(d) and a single amino acid peptide variant of QL9, called QL9-Y5-L(d). While 2C and most of the mutants had equal or higher affinity for the Y5 variant than for QL9, mutant CDR1beta exhibited 8-fold lower affinity for Y5 compared to QL9. To examine possible structural correlates of the thermodynamic and fine specificity signatures of the TCRs, the structure of unliganded QL9-L(d) was solved and compared to structures of the 2C TCR/QL9-L(d) complex and three high-affinity TCR/QL9-L(d) complexes. Our findings show that the QL9-L(d) complex does not undergo major conformational changes upon binding. Thus, subtle changes in individual CDRs account for the diverse thermodynamic and kinetic binding mechanisms and for the different peptide fine specificities.

Binding thermodynamics of the N-terminal peptide of the CCR5 coreceptor to HIV-1 envelope glycoprotein gp120

The initial events of HIV-1 cell infection involve the sequential binding of the viral envelope glycoprotein gp120 to the cellular CD4 receptor and the coreceptor, usually CCR5 or CXCR4. Binding to the coreceptor triggers the chain of events that culminates with the entry of the virus into the cell. In this process, the interaction of gp120 with the tyrosine-sulfated N-terminus of CCR5 is critical; however, this interaction has never been characterized at a quantitative or thermodynamic level. Here, we present the first thermodynamic analysis of the interaction of gp120 with the N-terminal peptide of the CCR5 coreceptor. Microcalorimetric titrations demonstrate that measurable binding of S22 peptide, a 22-amino acid tyrosine-sulfated peptide corresponding to the CCR5 N-terminus, requires prior binding of CD4 to gp120. The S22 peptide binds to the gp120-CD4 complex with a binding affinity of 4.5 x 10(5) M(-1) (K(d) = 2.2 microM) in an enthalpically and entropically favorable process. An identical peptide lacking the sulfated tyrosine residues is unable to bind the gp120-CD4 complex. These results indicate that the sulfated tyrosines contribute close to -3.5 kcal/mol to the Gibbs energy of binding. Furthermore, the S22 peptide is a competitive inhibitor of the 17b HIV-1 neutralizing antibody, which is known to bind to the CCR5 coreceptor site in gp120. Together, these results point toward compounds containing sulfated aromatic groups as potential inhibitors of viral entry. In analogy to existing inhibitors that bind to the CCR5 coreceptor directly, these compounds will accomplish the same result by binding to the coreceptor site in gp120.

Structural insight into the kinetics and DeltaCp of interactions between TEM-1 beta-lactamase and beta-lactamase inhibitory protein (BLIP)

In a previous study, we examined thermodynamic parameters for 20 alanine mutants in beta-lactamase inhibitory protein (BLIP) for binding to TEM-1 beta-lactamase. Here we have determined the structures of two thermodynamically distinctive complexes of BLIP mutants with TEM-1 beta-lactamase. The complex BLIP Y51A-TEM-1 is a tight binding complex with the most negative binding heat capacity change (DeltaG = approximately -13 kcal mol(-1) and DeltaCp = approximately -0.8 kcal mol(-1) K(-1)) among all of the mutants, whereas BLIP W150A-TEM-1 is a weak complex with one of the least negative binding heat capacity changes (DeltaG = approximately -8.5 kcal mol(-1) and DeltaCp = approximately -0.27 kcal mol(-1) K(-1)). We previously determined that BLIP Tyr51 is a canonical and Trp150 an anti-canonical TEM-1-contact residue, where canonical refers to the alanine substitution resulting in a matched change in the hydrophobicity of binding free energy. Structure determination indicates a rearrangement of the interactions between Asp49 of the W150A BLIP mutant and the catalytic pocket of TEM-1. The Asp49 of W150A moves more than 4 angstroms to form two new hydrogen bonds while losing four original hydrogen bonds. This explains the anti-canonical nature of the Trp150 to alanine substitution, and also reveals a strong long distance coupling between Trp150 and Asp49 of BLIP, because these two residues are more than 25 angstroms apart. Kinetic measurements indicate that the mutations influence the dissociation rate but not the association rate. Further analysis of the structures indicates that an increased number of interface-trapped water molecules correlate with poor interface packing in a mutant. It appears that the increase of interface-trapped water molecules is inversely correlated with negative binding heat capacity changes.

Structural basis of antibody-antigen interactions

Antibody molecules can be regarded as products of a protein engineering system for the generation of a virtually unlimited repertoire of complementary molecular surfaces. This extreme structural heterogeneity is required for recognition of the nearly infinite array of antigenic determinants. This chapter discusses the structures of antibodies and their specific recognition of antigens, the binding energetics of these interactions, the cross-reactivity and specificity of antibody-antigen interactions, the role of conformational flexibility in antigen recognition, and the structural basis of the antibody affinity maturation process.

E9-Im9 colicin DNase-immunity protein biomolecular association in water: a multiple-copy and accelerated molecular dynamics simulation study

Protein-protein transient and dynamic interactions underlie all biological processes. The molecular dynamics (MD) of the E9 colicin DNase protein, its Im9 inhibitor protein, and their E9-Im9 recognition complex are investigated by combining multiple-copy (MC) MD and accelerated MD (aMD) explicit-solvent simulation approaches, after validation with crystalline-phase and solution experiments. Im9 shows higher flexibility than its E9 counterpart. Im9 displays a significant reduction of backbone flexibility and a remarkable increase in motional correlation upon E9 association. Im9 loops 23-31 and 54-64 open with respect to the E9-Im9 X-ray structure and show high conformational diversity. Upon association a large fraction (approximately 20 nm2) of E9 and Im9 protein surfaces become inaccessible to water. Numerous salt bridges transiently occurring throughout our six 50 ns long MC-MD simulations are not present in the X-ray model. Among these Im9 Glu31-E9 Arg96 and Im9 Glu41-Lys89 involve interface interactions. Through the use of 10 ns of Im9 aMD simulation, we reconcile the largest thermodynamic impact measured for Asp51Ala mutation with Im9 structure and dynamics. Lys57 acts as an essential molecular switch to shift Im9 surface loop towards an ideal configuration for E9 inhibition. This is achieved by switching Asp60-Lys57 and Asp62-Lys57 hydrogen bonds to Asp51-Lys57 salt bridge. E9-Im9 recognition involves shifts of conformational distributions, reorganization of intramolecular hydrogen bond patterns, and formation of new inter- and intramolecular interactions. The description of key transient biological interactions can be significantly enriched by the dynamic and atomic-level information provided by computer simulations.

Protein–protein interaction and quaternary structure

Protein–protein recognition plays an essential role in structure and function. Specific non-covalent interactions stabilize the structure of macromolecular assemblies, exemplified in this review by oligomeric proteins and the capsids of icosahedral viruses. They also allow proteins to form complexes that have a very wide range of stability and lifetimes and are involved in all cellular processes. We present some of the structure-based computational methods that have been developed to characterize the quaternary structure of oligomeric proteins and other molecular assemblies and analyze the properties of the interfaces between the subunits. We compare the size, the chemical and amino acid compositions and the atomic packing of the subunit interfaces of protein–protein complexes, oligomeric proteins, viral capsids and protein–nucleic acid complexes. These biologically significant interfaces are generally close-packed, whereas the non-specific interfaces between molecules in protein crystals are loosely packed, an observation that gives a structural basis to specific recognition. A distinction is made within each interface between a core that contains buried atoms and a solvent accessible rim. The core and the rim differ in their amino acid composition and their conservation in evolution, and the distinction helps correlating the structural data with the results of site-directed mutagenesis and in vitro studies of self-assembly.

Structural mobility of the monomeric C-terminal domain of the HIV-1 capsid protein

The capsid protein of HIV-1 (p24) (CA) forms the mature capsid of the human immunodeficiency virus. Capsid assembly involves hexamerization of the N-terminal domain and dimerization of the C-terminal domain of CA (CAC), and both domains constitute potential targets for anti-HIV therapy. CAC homodimerization occurs mainly through its second helix, and it is abolished when its sole tryptophan is mutated to alanine. This mutant, CACW40A, resembles a transient monomeric intermediate formed during dimerization. Its tertiary structure is similar to that of the subunits in the dimeric, non-mutated CAC, but the segment corresponding to the second helix samples different conformations. The present study comprises a comprehensive examination of the CACW40A internal dynamics. The results obtained, with movements sampling a wide time regime (from pico- to milliseconds), demonstrate the high flexibility of the whole monomeric protein. The conformational exchange phenomena on the micro-to-millisecond time scale suggest a role for internal motions in the monomer-monomer interactions and, thus, flexibility of the polypeptide chain is likely to contribute to the ability of the protein to adopt different conformational states, depending on the biological environment.

Critical contribution of aromatic rings to specific recognition of polyether rings. The case of ciguatoxin CTX3C-ABC and its specific antibody 1C49

To address how proteins recognize polyether toxin compounds, we focused on the interaction between the ABC ring compound of ciguatoxin 3C and its specific antibody, 1C49. Surface plasmon resonance analyses indicated that Escherichia coli-expressed variable domain fragments (Fv) of 1C49 had the high affinity constants and slow dissociation constants typical of antigen-antibody interactions. Linear van't Hoff analyses suggested that the interaction is enthalpy-driven. We resolved the crystal structure of 1C49 Fv bound to ABC ring compound of ciguatoxin 3C at a resolution of 1.7A. The binding pocket of the antibody had many aromatic rings and bound the antigen by shape complementarity typical of hapten-antibody interactions. Three hydrogen bonds and many van der Waals interactions were present. We mutated several residues of the antibody to Ala, and we used surface plasmon resonance to analyze the interactions between the mutated antibodies and the antigen. This analysis identified Tyr-91 and Trp-96 in the light chain as hot spots for the interaction, and other residues made incremental contributions by conferring enthalpic advantages and reducing the dissociation rate constant. Systematic mutation of Tyr-91 indicated that CH-pi and pi-pi interactions between the aromatic ring at this site and the antigen made substantial contributions to the association, and van der Waals interactions inhibited dissociation, suggesting that aromaticity and bulkiness are critical for the specific recognition of polyether compounds by proteins.

The PD-1/PD-L1 complex resembles the antigen-binding Fv domains of antibodies and T cell receptors

Signaling through the programmed death 1 (PD-1) inhibitory receptor upon binding its ligand, PD-L1, suppresses immune responses against autoantigens and tumors and plays an important role in the maintenance of peripheral immune tolerance. Release from PD-1 inhibitory signaling revives "exhausted" virus-specific T cells in chronic viral infections. Here we present the crystal structure of murine PD-1 in complex with human PD-L1. PD-1 and PD-L1 interact through the conserved front and side of their Ig variable (IgV) domains, as do the IgV domains of antibodies and T cell receptors. This places the loops at the ends of the IgV domains on the same side of the PD-1/PD-L1 complex, forming a surface that is similar to the antigen-binding surface of antibodies and T cell receptors. Mapping conserved residues allowed the identification of residues that are important in forming the PD-1/PD-L1 interface. Based on the structure, we show that some reported loss-of-binding mutations involve the PD-1/PD-L1 interaction but that others compromise protein folding. The PD-1/PD-L1 interaction described here may be blocked by antibodies or by designed small-molecule drugs to lower inhibitory signaling that results in a stronger immune response. The immune receptor-like loops offer a new surface for further study and potentially the design of molecules that would affect PD-1/PD-L1 complex formation and thereby modulate the immune response.

Recommendations on measurement and analysis of results obtained on biological substances using isothermal titration calorimetry (IUPAC Technical Report)

Isothermal titration calorimetry (ITC) is widely used to determine the thermodynamics of biological interactions including protein-protein, small molecule-protein, protein-DNA, small molecule-DNA, and antigen-antibody interactions. An ITC measurement consists of monitoring the transfer of heat between an analyte solution in a sample vessel and a reference solution in a reference vessel upon injection of a small aliquot of titrant solution into the sample vessel at a fixed ITC operating temperature. A binding isotherm is generated from the heat-transferred-per-injection data and values for the binding constants, the apparent binding enthalpies, and the apparent ratio of the amount of titrant to analyte for the binding reaction are then determined from fits of a binding model, whether it is a single site, identical multi-site, or an interacting multi-site binding model, to the binding isotherm. Prior to the fitting procedure, corrections should be made for contributions from extraneous heat of mixing determined separately from injections of the titrant into just the dialysate/buffer solution. Ultra-high binding constants, which cannot be directly determined from an ITC measurement, can be determined by a displacement ITC method where injections of the tight-binding titrant into a solution of a weaker-binding titrant-analyte complex displaces the weaker-binding titrant from the complex. The Michaelis and catalytic constants can be determined for an enzyme reaction from injections of a substrate or enzyme titrant into an enzyme or substrate analyte solution. Several binding reactions are suggested to check the operating performance of the ITC. The reporting of ITC results must be specific with regard to the composition of the titrant and the analyte solutions, the temperature, and the model used in the analysis.

Amino acid side chain interactions in the presence of salts

The effects of salt on the intermolecular interactions between polar/charged amino acids are investigated through molecular dynamics simulations. The mean forces and associated potentials are calculated for NaCl salt in the 0-2 M concentration range at 298 K. It is found that the addition of salt may stabilize or destabilize the interactions, depending on the nature of the interacting molecules. The degree of (de)stabilization is quantified, and the origin of the salt-dependent modulation is discussed based upon an analysis of solvent density profiles. To gain insight into the molecular origin of the salt modulation, spatial distribution functions (sdf's) are calculated, revealing a high degree of solvent structuredness in all cases. The peaks in the sdf's are consistent with long-range hydrogen-bonding networks connecting the solute hydrophilic groups, and that contribute to their intermolecular solvent-induced forces. The restructuring of water around the solutes as they dissociate from close contact is analyzed. This analysis offers clues on how the solvent structure modulates the effective intermolecular interactions in complex solutes. This modulation results from a critical balance between bulk electrostatic forces and those exerted by (i) the water molecules in the structured region between the monomers, which is disrupted by ions that transiently enter the hydration shells, and (ii) the ions in the hydration shells in direct interactions with the solutes. The implications of these findings in protein/ligand (noncovalent) association/dissociation mechanisms are briefly discussed.

Thermodynamic studies on the interaction of antibodies with beta-amyloid peptide

Antibodies against beta-amyloid peptides (Abetas) are considered an important therapeutic opportunity in Alzheimer's disease. Despite the vast interest in Abeta no thermodynamic data on the interaction of antibodies with Abeta are available as yet. In the present study we use isothermal titration calorimetry (ITC) and surface plasmon resonance to provide a quantitative thermodynamic analysis of the interaction between soluble monomeric Abeta(1-40) and mouse monoclonal antibodies (mAb). Using four different antibodies directed against the N-terminal, middle, and C-terminal Abeta epitopes, we measured the thermodynamic parameters for the binding to Abeta. Each antibody species was found to have two independent and equal binding sites for Abeta with binding constants in the range of 10(7) to 10(8) M(-1). The binding reaction was essentially enthalpy driven with a reaction enthalpy of DeltaH(0)(Abeta) approximately -19 to -8 kcal/mol, indicating the formation of tight complexes. The loss in conformational freedom was supported by negative values for the reaction entropy DeltaS(0)(Abeta). We also measured the heat capacity change of the 1mAb:2Abeta reaction. DeltaC(0)(p, abeta) was large and negative but could not be explained exclusively by the hydrophobic effect. The free energy of binding was found to be linearly correlated with the size of the epitope.

Flexibility and Conformational Entropy in Protein-Protein Binding

To better understand the interplay between protein-protein binding and protein dynamics, we analyzed molecular dynamics simulations of 17 protein-protein complexes and their unbound components. Complex formation does not restrict the conformational freedom of the partner proteins as a whole, but, rather, it leads to a redistribution of dynamics. We calculate the change in conformational entropy for seven complexes with quasiharmonic analysis. We see significant loss, but also increased or unchanged conformational entropy. Where comparison is possible, the results are consistent with experimental data. However, stringent error estimates based on multiple independent simulations reveal large uncertainties that are usually overlooked. We observe substantial gains of pseudo entropy in individual partner proteins, and we observe that all complexes retain residual stabilizing intermolecular motions. Consequently, protein flexibility has an important influence on the thermodynamics of binding and may disfavor as well as favor association. These results support a recently proposed unified model for flexible protein-protein association.

Structural basis of affinity maturation and intramolecular cooperativity in a protein-protein interaction

Although protein-protein interactions are involved in nearly all cellular processes, general rules for describing affinity and selectivity in protein-protein complexes are lacking, primarily because correlations between changes in protein structure and binding energetics have not been well determined. Here, we establish the structural basis of affinity maturation for a protein-protein interaction system that we had previously characterized energetically. This model system exhibits a 1500-fold affinity increase. Also, its affinity maturation is restricted by negative intramolecular cooperativity. With three complex and six unliganded variant X-ray crystal structures, we provide molecular snapshots of protein interface remodeling events that span the breadth of the affinity maturation process and present a comprehensive structural view of affinity maturation. Correlating crystallographically observed structural changes with measured energetic changes reveals molecular bases for affinity maturation, intramolecular cooperativity, and context-dependent binding.

Membrane chromatographic immunoassay method for rapid quantitative analysis of specific serum antibodies

This paper discusses a membrane chromatographic immunoassay method for rapid detection and quantitative analysis of specific serum antibodies. A type of polyvinylidine fluoride (PVDF) microfiltration membrane was used in the method for its ability to reversibly and specifically bind IgG antibodies from antiserum samples by hydrophobic interaction. Using this form of selective antibody binding and enrichment an affinity membrane with antigen binding ability was obtained in-situ. This was done by passing a pulse of diluted antiserum sample through a stack of microporous PVDF membranes. The affinity membrane thus formed was challenged with a pulse of antigen solution and the amount of antigen bound was accurately determined using chromatographic methods. The antigen binding correlated well with the antibody loading on the membrane. This method is direct, rapid and accurate, does not involve any chemical reaction, and uses very few reagents. Moreover, the same membrane could be repeatedly used for sequential immunoassays on account of the reversible nature of the antibody binding. Proof of concept of this method is provided using human hemoglobin as model antigen and rabbit antiserum against human hemoglobin as the antibody source.

Design and validation of a synthetic VH repertoire with tailored diversity for protein recognition

Previous studies have indicated differences in the specificity-determining residues (SDRs) of antibodies that recognize haptens, peptides, or proteins. Here, we designed a V(H) repertoire based on the human scaffold 3-23/J(H)4 and diversification of high and medium-usage SDRs of anti-protein and anti-peptide antibodies. The repertoire was synthesized by overlapping polymerase chain reaction (PCR) and combined with the V(L) chain of the anti-hen egg-white lysozyme (HEL) antibody D1.3. The resulting chimeric single-chain Fv fragments (scFvs) phage-displayed library was panned in HEL-coated immunotubes. After two rounds of selection under non-stringent conditions, that is, trypsinization after 2 h of incubation at room temperature, 63 of 167 clones analyzed (38%) were found to express scFvs specific to HEL. Twenty clones were characterized by DNA sequencing resulting in 10 unique scFvs. Interestingly, the panel of unique scFvs was highly diverse, with V(H) sequences differing in 16 of the 17 positions variegated in the repertoire. Thus, diverse chemico-physical and structural solutions were selected from the library, even when the V(H) repertoire was constrained by the V(L) chain of D1.3 to yield binders against a definite region of HEL surface. The more often selected scFvs, namely H6-1 and B7-1, which differed in eight SDRs, showed levels of expression in E. coli TG1 strain, 6 and 10 times higher than the parental D1.3 Fv fragment, respectively. Dissociation constants (K(Ds)) measured in the BIAcore were 11 and 6.6 nM for H6-1 and B7-1, respectively. These values compared well to the K(D) of 4.7 nM measured for D1.3, indicating that the V(H) repertoire here designed is a valuable source of diverse, well-expressed and high affinity V(H) domains.

A high-throughput fluorescence polarization assay specific to the CD4 binding site of HIV-1 glycoproteins based on a fluorescein-labelled CD4 mimic

The three-dimensional structure of CD4M33, a mimic of the host-cell receptor-antigen CD4 and a powerful inhibitor of CD4-gp120 (viral envelope glycoprotein 120) interaction and HIV-1 entry into cells [Martin, Stricher, Misse, Sironi, Pugniere, Barthe, Prado-Gotor, Freulon, Magne, Roumestand et al. (2003) Nat. Biotechnol. 21, 71-76], was solved by 1H-NMR and its structure was modelled in its complex with gp120. In this complex, CD4M33 binds in a CD4-like mode and inserts its unnatural and prominent Bip23 (biphenylalanine-23) side-chain into the gp120 interior 'Phe43 cavity', thus filling its volume. CD4M33 was specifically labelled with fluorescein and shown by fluorescence anisotropy to bind to different gp120 glycoproteins with dissociation constants in the nanomolar range. Fluorescent CD4M33 was also used in a miniaturized 384-well-plate assay to study direct binding to a large panel of gp120 glycoproteins and in a competition assay to study binding of CD4 or other ligands targeting the CD4 binding site of gp120. Furthermore, by using the fluorescently labelled CD4M33 and the [Phe23]M33 mutant, which possesses a natural Phe23 residue and thus cannot penetrate the gp120 Phe43 cavity, we show that a recently discovered small-molecule-entry inhibitor, BMS-378806, does not target the CD4 binding site nor the Phe43 cavity of gp120. The fluorescently labelled CD4M33 mimic, its mutants and their derivatives represent useful tools with which to discover new molecules which target the CD4 binding site and/or the Phe43 cavity of gp120 glycoproteins in a high-throughput fluorescence-polarization assay and to characterize their mechanism of action.

Magnitude of the hydrophobic effect at central versus peripheral sites in protein-protein interfaces

Hydrophobic interactions are essential for stabilizing protein-protein complexes, whose interfaces generally consist of a central cluster of hot spot residues surrounded by less important peripheral residues. According to the O-ring hypothesis, a condition for high affinity binding is solvent exclusion from interacting residues. This hypothesis predicts that the hydrophobicity at the center is significantly greater than at the periphery, which we estimated at 21 cal mol(-1) A(-2). To measure the hydrophobicity at the center, structures of an antigen-antibody complex where a buried phenylalanine was replaced by smaller hydrophobic residues were determined. By correlating structural changes with binding free energies, we estimate the hydrophobicity at this central site to be 46 cal mol(-1) A(-2), twice that at the periphery. This context dependence of the hydrophobic effect explains the clustering of hot spots at interface centers and has implications for hot spot prediction and the design of small molecule inhibitors.