Molecular evolution of affinity and flexibility in the immune system

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

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

Structural basis for germline antibody recognition of HIV-1 immunogens

Efforts to elicit broadly neutralizing antibodies (bNAbs) against HIV-1 require understanding germline bNAb recognition of HIV-1 envelope glycoprotein (Env). The VRC01-class bNAb family derived from the VH1-2*02 germline allele arose in multiple HIV-1–infected donors, yet targets the CD4-binding site on Env with common interactions. Modified forms of the 426c Env that activate germline-reverted B cell receptors are candidate immunogens for eliciting VRC01-class bNAbs. We present structures of germline-reverted VRC01-class bNAbs alone and complexed with 426c-based gp120 immunogens. Germline bNAb–426c gp120 complexes showed preservation of VRC01-class signature residues and gp120 contacts, but detectably different binding modes compared to mature bNAb-gp120 complexes. Unlike typical antibody-antigen interactions, VRC01–class germline antibodies exhibited preformed antigen-binding conformations for recognizing immunogens. Affinity maturation introduced substitutions increasing induced-fit recognition and electropositivity, potentially to accommodate negatively-charged complex-type N-glycans on gp120. These results provide general principles relevant to the unusual evolution of VRC01–class bNAbs and guidelines for structure-based immunogen design.

DOI:http://dx.doi.org/10.7554/eLife.13783.001

When human immunodeficiency virus-1 (HIV-1) infects humans it can cause a serious disease that damages the immune system. Currently there is no cure for this disease and there are no vaccines available to halt the spread of the virus. Researchers are hoping to be able to develop a single vaccine that can protect individuals against every form (or strain) of HIV-1, but this has proved difficult because many different versions of the virus exist.

An effective vaccine triggers long-lasting immunity to a particular virus or microbe by activating the production of proteins called antibodies that identify and help to destroy the threat. Research has shown that most individuals infected with HIV-1 produce antibodies that can only recognize a few HIV strains. However, there are rare individuals who produce “broadly neutralizing antibodies”; that is, antibodies that can recognize and help to kill 90% or more of HIV-1 strains. Understanding how broadly neutralizing antibodies are produced in infected individuals may aid the development of a vaccine that can protect others from the many circulating strains of HIV.

When an individual encounters a virus, immature antibodies are modified to generate mature antibodies that bind more effectively to specific virus proteins. Here, Scharf et al. investigated how a class of broadly neutralizing antibodies called VRC01-class antibodies, which bind to an HIV protein called gp120, are produced. The experiments used a technique called X-ray crystallography to reveal the three-dimensional structures of immature versions of these antibodies when they are bound to gp120.

Scharf et al. discovered that, unlike most antibodies, the overall final structure of VRC01 antibodies is formed before the antibody matures. Instead of making large changes to the structure of these antibodies, the maturation process makes VRC01-class antibodies become more positively charged, which allows them to bind to gp120 proteins on a wider variety of HIV viruses. These findings suggest that it may be possible to use modified gp120 proteins in vaccines to trigger the production of broadly neutralizing antibodies against HIV.

DOI:http://dx.doi.org/10.7554/eLife.13783.002

Ancestral Protein Reconstruction Yields Insights into Adaptive Evolution of Binding Specificity in Solute-Binding Proteins

The promiscuous functions of proteins are an important reservoir of functional novelty in protein evolution, but the molecular basis for binding promiscuity remains elusive. We used ancestral protein reconstruction to experimentally characterize evolutionary intermediates in the functional expansion of the polar amino acid-binding protein family, which has evolved to bind a variety of amino acids with high affinity and specificity. High-resolution crystal structures of an ancestral arginine-binding protein in complex with l-arginine and l-glutamine show that the promiscuous binding of l-glutamine is enabled by multi-scale conformational plasticity, water-mediated interactions, and selection of an alternative conformational substate productive for l-glutamine binding. Evolution of specialized glutamine-binding proteins from this ancestral protein was achieved by displacement of water molecules from the protein-ligand interface, reducing the entropic penalty associated with the promiscuous interaction. These results provide a structural and thermodynamic basis for the co-option of a promiscuous interaction in the evolution of binding specificity.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Vaccine-elicited antibody that neutralizes H5N1 influenza and variants binds the receptor site and polymorphic sites

Antigenic drift of circulating seasonal influenza viruses necessitates an international vaccine effort to reduce the impact on human health. A critical feature of the seasonal vaccine is that it stimulates an already primed immune system to diversify memory B cells to recognize closely related, but antigenically distinct, influenza glycoproteins (hemagglutinins). Influenza pandemics arise when hemagglutinins to which no preexisting adaptive immunity exists acquire the capacity to infect humans. Hemagglutinin 5 is one subtype to which little preexisting immunity exists and is only a few acquired mutations away from the ability to transmit efficiently between ferrets, and possibly humans. Here, we describe the structure and molecular mechanism of neutralization by H5.3, a vaccine-elicited antibody that neutralizes hemagglutinin 5 viruses and variants with expanded host range. H5.3 binds in the receptor-binding site, forming contacts that recapitulate many of the sialic acid interactions, as well as multiple peripheral interactions, yet is not sensitive to mutations that alter sialic acid binding. H5.3 is highly specific for a subset of H5 strains, and this specificity arises from interactions to the periphery of the receptor-binding site. H5.3 is also extremely potent, despite retaining germ line-like conformational flexibility.

Specificity, polyspecificity, and heterospecificity of antibody-antigen recognition

The concept of antibody specificity is analyzed and shown to reside in the ability of an antibody to discriminate between two antigens. Initially, antibody specificity was attributed to sequence differences in complementarity determining regions (CDRs), but as increasing numbers of crystallographic antibody-antigen complexes were elucidated, specificity was analyzed in terms of six antigen-binding regions (ABRs) that only roughly correspond to CDRs. It was found that each ABR differs significantly in its amino acid composition and tends to bind different types of amino acids at the surface of proteins. In spite of these differences, the combined preference of the six ABRs does not allow epitopes to be distinguished from the rest of the protein surface. These findings explain the poor success of past and newly proposed methods for predicting protein epitopes. Antibody polyspecificity refers to the ability of one antibody to bind a large variety of epitopes in different antigens, and this property explains how the immune system develops an antibody repertoire that is able to recognize every antigen the system is likely to encounter. Antibody heterospecificity arises when an antibody reacts better with another antigen than with the one used to raise the antibody. As a result, an antibody may sometimes appear to have been elicited by an antigen with which it is unable to react. The implications of antibody polyspecificity and heterospecificity in vaccine development are pointed out.

Identification of Tyr residues that enhance folate substrate binding and constrain oscillation of the proton-coupled folate transporter (PCFT-SLC46A1)

The proton-coupled folate transporter (PCFT) mediates intestinal folate absorption and transport of folates across the choroid plexus. This study focuses on the role of Tyr residues in PCFT function. The substituted Cys-accessibility method identified four Tyr residues (Y291, Y362, Y315, and Y414) that are accessible to the extracellular compartment; three of these (Y291, Y362, and Y315) are located within or near the folate binding pocket. When the Tyr residues were replaced with Cys or Ala, these mutants showed similar (up to 6-fold) increases in influx Vmax and Kt/Ki for [(3)H]methotrexate and [(3)H]pemetrexed. When the Tyr residues were replaced with Phe, these changes were moderated or absent. When Y315A PCFT was used as representative of the mutants and [(3)H]pemetrexed as the transport substrate, this substitution did not increase the efflux rate constant. Furthermore, neither influx nor efflux mediated by Y315A PCFT was transstimulated by the presence of substrate in the opposite compartment; however, substantial bidirectional transstimulation of transport was mediated by wild-type PCFT. This resulted in a threefold greater efflux rate constant for cells that express wild-type PCFT than for cells that express Y315 PCFT under exchange conditions. These data suggest that these Tyr residues, possibly through their rigid side chains, secure the carrier in a high-affinity state for its folate substrates. However, this may be achieved at the expense of constraining the carrier's mobility, thereby decreasing the rate at which the protein oscillates between its conformational states. The Vmax generated by these Tyr mutants may be so rapid that further augmentation during transstimulation may not be possible.

Back to the Future: Nuclear Magnetic Resonance and Bioinformatics Studies on Intrinsically Disordered Proteins

From the 1970s to the present, regions of missing electron density in protein structures determined by X-ray diffraction and the characterization of the functions of these regions have suggested that not all protein regions depend on prior 3D structure to carry out function. Motivated by these observations, in early 1996 we began to use bioinformatics approaches to study these intrinsically disordered proteins (IDPs) and IDP regions. At just about the same time, several laboratory groups began to study a collection of IDPs and IDP regions using nuclear magnetic resonance. The temporal overlap of the bioinformatics and NMR studies played a significant role in the development of our understanding of IDPs. Here the goal is to recount some of this history and to project from this experience possible directions for future work.

Antigen exposure leads to rigidification of germline antibody combining site

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

Origination of the Protein Fold Repertoire from Oily Pluripotent Peptides

While the repertoire of protein folds that exists today underlies most of life’s capabilities, our mechanistic picture of protein fold origination is incomplete. This paper discusses a hypothetical mechanism for the emergence of the protein fold repertoire from highly dynamic and collapsed peptides, exemplified by peptides with high oil content or hydrophobicity. These peptides are called pluripotent to emphasize their capacity to evolve into numerous folds transiently available to them. As evidence, the paper will discuss previous simulation work on the superior fold evolvability of oily peptides, trace (“fossil”) evidence within proteomes seen today, and a general relationship between protein dynamism and evolvability. Aside from implications on the origination of protein folds, the hypothesis implies that the vanishing utility of a random peptide in protein origination may be relatively exaggerated, as some random peptides with a certain composition (e.g., oily) may fare better than others. In later sections, the hypothesis is discussed in the context of existing discussions regarding the spontaneous origination of biomolecules.

Dynamic New World: Refining Our View of Protein Structure, Function and Evolution

Proteins are crucial to the functioning of all lifeforms. Traditional understanding posits that a single protein occupies a single structure ("fold"), which performs a single function. This view is radically challenged with the recognition that high structural dynamism-the capacity to be extra "floppy"-is more prevalent in functional proteins than previously assumed. As reviewed here, this dynamic take on proteins affects our understanding of protein "structure", function, and evolution, and even gives us a glimpse into protein origination. Specifically, this review will discuss historical developments concerning protein structure, and important new relationships between dynamism and aspects of protein sequence, structure, binding modes, binding promiscuity, evolvability, and origination. Along the way, suggestions will be provided for how key parts of textbook definitions-that so far have excluded membership to intrinsically disordered proteins (IDPs)-could be modified to accommodate our more dynamic understanding of proteins.

Backbone flexibility of CDR3 and immune recognition of antigens

Conformational entropy is an important component of protein-protein interactions; however, there is no reliable method for computing this parameter. We have developed a statistical measure of residual backbone entropy in folded proteins by using the ϕ-ψ distributions of the 20 amino acids in common secondary structures. The backbone entropy patterns of amino acids within helix, sheet or coil form clusters that recapitulate the branching and hydrogen bonding properties of the side chains in the secondary structure type. The same types of residues in coil and sheet have identical backbone entropies, while helix residues have much smaller conformational entropies. We estimated the backbone entropy change for immunoglobulin complementarity-determining regions (CDRs) from the crystal structures of 34 low-affinity T-cell receptors and 40 high-affinity Fabs as a result of the formation of protein complexes. Surprisingly, we discovered that the computed backbone entropy loss of only the CDR3, but not all CDRs, correlated significantly with the kinetic and affinity constants of the 74 selected complexes. Consequently, we propose a simple algorithm to introduce proline mutations that restrict the conformational flexibility of CDRs and enhance the kinetics and affinity of immunoglobulin interactions. Combining the proline mutations with rationally designed mutants from a previous study led to 2400-fold increase in the affinity of the A6 T-cell receptor for Tax-HLAA2. However, this mutational scheme failed to induce significant binding changes in the already-high-affinity C225-Fab/huEGFR interface. Our results will serve as a roadmap to formulate more effective target functions to design immune complexes with improved biological functions.

Large-scale analysis of somatic hypermutations in antibodies reveals which structural regions, positions and amino acids are modified to improve affinity

The principles of affinity maturation of antibodies (Abs), which underlies B cell-mediated immunity, are still under debate. It is unclear whether the antigen (Ag) binding site is a preferred target for mutations, and what the role of activation-induced deaminase (AID) hotspots is in this process. Here we report a structural analysis of 3495 residues that have been replaced through somatic hypermutations (SHMs) in 196 Abs. We show that there is no correlation between the propensity of an amino acid to be in AID hotspot and the probability that it is replaced during the SHM process. Although AID hotspots may be necessary to enable SHMs, they are not a major driving force in determining which residues are mutated. We identified Ab positions that are highly mutated and significantly affect binding. The effect of mutation on binding energy is a major factor in determining which structural regions of the Ab are mutated. There is a clear preference for mutations at the Ag-binding site. However, positions outside this region that also affect binding are often preferred targets for SHMs. As for amino acid preferences, a general trend during SHM is to make Ab-Ag interfaces more similar to protein-protein interfaces in general. In different regions of the Ab, there are different sets of preferences for amino acid substitution. This mapping improves our understanding of Ab affinity maturation and may assist in Ab engineering.

Bond elasticity controls molecular recognition specificity in antibody-antigen binding

Force-spectroscopy experiments make it possible to characterize single ligand-receptor pairs. Here we measure the spectrum of bond strengths and flexibilities in antibody-antigen interactions using optical tweezers. We characterize the mechanical evolution of polyclonal antibodies generated under infection and the ability of a monoclonal antibody to cross-react against different antigens. Our results suggest that bond flexibility plays a major role in remodeling antibody-antigen bonds in order to improve recognition during the maturation of the humoral immune system.

Energetic and dynamic aspects of the affinity maturation process: characterizing improved variants from the bevacizumab antibody with molecular simulations

Antibody affinity maturation is one of the fundamental processes of immune defense against invading pathogens. From the biological point of view, the clonal selection hypothesis represents the most accepted mechanism to explain how mutations increasing the affinity for target antigens are introduced and selected in antibody molecules. However, understanding at the molecular level how protein modifications, such as point mutation, can modify and modulate the affinity of an antibody for its antigen is still a major open issue in molecular biology. In this paper, we address various aspects of this problem by analyzing and comparing atomistic simulations of 17 variants of the bevacizumab antibody, all directed against the common target protein VEGF-A. In particular, we examine MD-based descriptors of the internal energetics and dynamics of mutated antibodies and their possible correlations with experimentally determined affinities for the antigens. Our results show that affinity improvement is correlated with a variation of the internal stabilization energy of the antibody molecule when bound to the antigen, compensated by the variation in the interaction energy between the antigen and the antibody, paralleled by an overall modulation of internal coordination within the antibody molecular structure. A possible model of the mechanism of rigidification and of the main residues involved is proposed. Overall, our results can help in understanding the molecular determinants of antigen recognition and have implications in the rational design of new antibodies with optimized affinities.

Structure and dynamics of the gp120 V3 loop that confers noncompetitive resistance in R5 HIV-1(JR-FL) to maraviroc

Maraviroc, an (HIV-1) entry inhibitor, binds to CCR5 and efficiently prevents R5 human immunodeficiency virus type 1 (HIV-1) from using CCR5 as a coreceptor for entry into CD4(+) cells. However, HIV-1 can elude maraviroc by using the drug-bound form of CCR5 as a coreceptor. This property is known as noncompetitive resistance. HIV-1(V3-M5) derived from HIV-1(JR-FLan) is a noncompetitive-resistant virus that contains five mutations (I304V/F312W/T314A/E317D/I318V) in the gp120 V3 loop alone. To obtain genetic and structural insights into maraviroc resistance in HIV-1, we performed here mutagenesis and computer-assisted structural study. A series of site-directed mutagenesis experiments demonstrated that combinations of V3 mutations are required for HIV-1(JR-FLan) to replicate in the presence of 1 µM maraviroc, and that a T199K mutation in the C2 region increases viral fitness in combination with V3 mutations. Molecular dynamic (MD) simulations of the gp120 outer domain V3 loop with or without the five mutations showed that the V3 mutations induced (i) changes in V3 configuration on the gp120 outer domain, (ii) reduction of an anti-parallel β-sheet in the V3 stem region, (iii) reduction in fluctuations of the V3 tip and stem regions, and (iv) a shift of the fluctuation site at the V3 base region. These results suggest that the HIV-1 gp120 V3 mutations that confer maraviroc resistance alter structure and dynamics of the V3 loop on the gp120 outer domain, and enable interactions between gp120 and the drug-bound form of CCR5.

Ensemble-based characterization of unbound and bound states on protein energy landscape

Physicochemical description of numerous cell processes is fundamentally based on the energy landscapes of protein molecules involved. Although the whole energy landscape is difficult to reconstruct, increased attention to particular targets has provided enough structures for mapping functionally important subspaces associated with the unbound and bound protein structures. The subspace mapping produces a discrete representation of the landscape, further called energy spectrum. We compiled and characterized ensembles of bound and unbound conformations of six small proteins and explored their spectra in implicit solvent. First, the analysis of the unbound-to-bound changes points to conformational selection as the binding mechanism for four proteins. Second, results show that bound and unbound spectra often significantly overlap. Moreover, the larger the overlap the smaller the root mean square deviation (RMSD) between the bound and unbound conformational ensembles. Third, the center of the unbound spectrum has a higher energy than the center of the corresponding bound spectrum of the dimeric and multimeric states for most of the proteins. This suggests that the unbound states often have larger entropy than the bound states. Fourth, the exhaustively long minimization, making small intrarotamer adjustments (all-atom RMSD ≤ 0.7 Å), dramatically reduces the distance between the centers of the bound and unbound spectra as well as the spectra extent. It condenses unbound and bound energy levels into a thin layer at the bottom of the energy landscape with the energy spacing that varies between 0.8-4.6 and 3.5-10.5 kcal/mol for the unbound and bound states correspondingly. Finally, the analysis of protein energy fluctuations showed that protein vibrations itself can excite the interstate transitions, including the unbound-to-bound ones.

Fragment-based drug discovery using a multidomain, parallel MD-MM/PBSA screening protocol

We have developed a rigorous computational screening protocol to identify novel fragment-like inhibitors of N(5)-CAIR mutase (PurE), a key enzyme involved in de novo purine synthesis that represents a novel target for the design of antibacterial agents. This computational screening protocol utilizes molecular docking, graphics processing unit (GPU)-accelerated molecular dynamics, and Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) free energy estimations to investigate the binding modes and energies of fragments in the active sites of PurE. PurE is a functional octamer comprised of identical subunits. The octameric structure, with its eight active sites, provided a distinct advantage in these studies because, for a given simulation length, we were able to place eight separate fragment compounds in the active sites to increase the throughput of the MM/PBSA analysis. To validate this protocol, we have screened an in-house fragment library consisting of 352 compounds. The theoretical results were then compared with the results of two experimental fragment screens, Nuclear Magnetic Resonance (NMR) and Surface Plasmon Resonance (SPR) binding analyses. In these validation studies, the protocol was able to effectively identify the competitive binders that had been independently identified by experimental testing, suggesting the potential utility of this method for the identification of novel fragments for future development as PurE inhibitors.

Investigating allostery in molecular recognition: insights from a computational study of multiple antibody-antigen complexes

Antibody-antigen recognition plays a key role in the immune response against pathogens. Here, we have investigated various aspects of this problem by analyzing a large and diverse set of antibodies and their respective complexes with protein antigens through atomistic simulations. Common features of antibody response to the presence of antigens are elucidated by the analysis of the proteins' internal dynamics and coordination in different ligand states, combined with the analysis of the interaction networks implicated in the stabilization of functional structures. The use of a common structural reference reveals preferential changes in the dynamic coordination and intramolecular interaction networks induced by antigen binding and shared by all antibodies. Such changes propagate from the binding region through the whole immunoglobulin domains. Overall, complexed antibodies show more diffuse networks of nonbonded interactions and a general higher internal dynamic coordination, which preferentially involve the immunoglobulin (Ig) domains of the heavy chain. The combined results provide atomistic insights into the correlations between the modulation of conformational dynamics, structural stability, and allosteric signal transduction. In particular, the results suggest that specific networks of residues, shared among all the analyzed proteins, define the molecular pathways by which antibody structures respond to antigen binding. Our studies may have implications in practical use, such as the rational design of antibodies with specifically modulated antigen-binding affinities.

A single E105K mutation far from the active site of influenza B virus neuraminidase contributes to reduced susceptibility to multiple neuraminidase-inhibitor drugs

Drugs inhibiting the enzymatic activity of influenza virus neuraminidase (NA) are the cornerstone of therapy for influenza virus infection. The emergence of drug-resistant variants may limit the benefits of antiviral therapy. Here we report the recovery of an influenza B virus with reduced susceptibilities to NA inhibitors from a human patient with no history of antiviral drug treatment. The virus, designated B/Kochi/61/2011, was isolated by inoculating Madin-Darby canine kidney (MDCK) cells with respiratory specimens from the patient. NA inhibition assays demonstrated that the B/Kochi/61/2011 isolate showed a remarkable reduction in susceptibility to peramivir. The isolate also exhibited low to moderately reduced sensitivity to oseltamivir, laninamivir, and zanamivir. A sequence analysis of viruses propagated in MDCK cells revealed that the isolate contained a mutation (E105K) not previously associated with reduced susceptibility to NA inhibitors. However, pyrosequencing analysis showed that the NA E105K mutation was below a detectable level in the original clinical specimens, suggesting that the mutant virus may be preferably selected during propagation in MDCK cells. Analysis of the three-dimensional model of E105 and K105 NAs with peramivir suggested that the E105K mutation at the monomer-monomer interface of the NA tetramer may destabilize the tetrameric form of NA, leading to decreased susceptibility to NA inhibitors. These results have implications for understanding the mechanism of resistance against NA-inhibitor drugs.

The role of flexibility and conformational selection in the binding promiscuity of PDZ domains

In molecular recognition, it is often the case that ligand binding is coupled to conformational change in one or both of the binding partners. Two hypotheses describe the limiting cases involved; the first is the induced fit and the second is the conformational selection model. The conformational selection model requires that the protein adopts conformations that are similar to the ligand-bound conformation in the absence of ligand, whilst the induced-fit model predicts that the ligand-bound conformation of the protein is only accessible when the ligand is actually bound. The flexibility of the apo protein clearly plays a major role in these interpretations. For many proteins involved in signaling pathways there is the added complication that they are often promiscuous in that they are capable of binding to different ligand partners. The relationship between protein flexibility and promiscuity is an area of active research and is perhaps best exemplified by the PDZ domain family of proteins. In this study we use molecular dynamics simulations to examine the relationship between flexibility and promiscuity in five PDZ domains: the human Dvl2 (Dishevelled-2) PDZ domain, the human Erbin PDZ domain, the PDZ1 domain of InaD (inactivation no after-potential D protein) from fruit fly, the PDZ7 domain of GRIP1 (glutamate receptor interacting protein 1) from rat and the PDZ2 domain of PTP-BL (protein tyrosine phosphatase) from mouse. We show that despite their high structural similarity, the PDZ binding sites have significantly different dynamics. Importantly, the degree of binding pocket flexibility was found to be closely related to the various characteristics of peptide binding specificity and promiscuity of the five PDZ domains. Our findings suggest that the intrinsic motions of the apo structures play a key role in distinguishing functional properties of different PDZ domains and allow us to make predictions that can be experimentally tested.

Proteins that are capable of binding to many different ligands are said to have broad specificity. This is sometimes also referred to as promiscuity. Whether a protein is promiscuous or not can sometimes be readily explained by the structure of the protein and the ligand in terms of electrostatic and steric effects. Sometimes however, this simple interpretation can struggle to explain the experimentally observed data. A prominent case in point is the PDZ domains. These small protein domains bind to unstructured regions of other proteins and are involved in many signaling pathways. Some PDZ domains appear to be more promiscuous than others, but this has been difficult to explain purely on the basis of the composition of residues in the binding groove. In this work we examine the dynamics and conformational flexibility of five key PDZ domains and demonstrate that despite similar folds, these proteins can exhibit quite different dynamics. Furthermore the difference in the dynamic behavior appears to correlate with the observed promiscuity. Our findings suggest that knowledge of the dynamic behavior of the PDZs can be used to rationalize the extent of expected promiscuity. Such knowledge will be critical for drug design against PDZ domains.

Analysis of the relationships between evolvability, thermodynamics, and the functions of intrinsically disordered proteins/regions

The evolvability of proteins is not only restricted by functional and structural importance, but also by other factors such as gene duplication, protein stability, and an organism's robustness. Recently, intrinsically disordered proteins (IDPs)/regions (IDRs) have been suggested to play a role in facilitating protein evolution. However, the mechanisms by which this occurs remain largely unknown. To address this, we have systematically analyzed the relationship between the evolvability, stability, and function of IDPs/IDRs. Evolutionary analysis shows that more recently emerged IDRs have higher evolutionary rates with more functional constraints relaxed (or experiencing more positive selection), and that this may have caused accelerated evolution in the flanking regions and in the whole protein. A systematic analysis of observed stability changes due to single amino acid mutations in IDRs and ordered regions shows that while most mutations induce a destabilizing effect in proteins, mutations in IDRs cause smaller stability changes than in ordered regions. The weaker impact of mutations in IDRs on protein stability may have advantages for protein evolvability in the gain of new functions. Interestingly, however, an analysis of functional motifs in the PROSITE and ELM databases showed that motifs in IDRs are more conserved, characterized by smaller entropy and lower evolutionary rate, than in ordered regions. This apparently opposing evolutionary effect may be partly due to the flexible nature of motifs in IDRs, which require some key amino acid residues to engage in tighter interactions with other molecules. Our study suggests that the unique conformational and thermodynamic characteristics of IDPs/IDRs play an important role in the evolvability of proteins to gain new functions.

Computer-aided antibody design

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

Cytokinergic IgE Action in Mast Cell Activation

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

Molecular dynamics simulation in virus research

Virus replication in the host proceeds by chains of interactions between viral and host proteins. The interactions are deeply influenced by host immune molecules and anti-viral compounds, as well as by mutations in viral proteins. To understand how these interactions proceed mechanically and how they are influenced by mutations, one needs to know the structures and dynamics of the proteins. Molecular dynamics (MD) simulation is a powerful computational method for delineating motions of proteins at an atomic-scale via theoretical and empirical principles in physical chemistry. Recent advances in the hardware and software for biomolecular simulation have rapidly improved the precision and performance of this technique. Consequently, MD simulation is quickly extending the range of applications in biology, helping to reveal unique features of protein structures that would be hard to obtain by experimental methods alone. In this review, we summarize the recent advances in MD simulations in the study of virus–host interactions and evolution, and present future perspectives on this technique.

Structural dynamics of HIV-1 envelope Gp120 outer domain with V3 loop

BACKGROUND

The net charge of the hypervariable V3 loop on the HIV-1 envelope gp120 outer domain plays a key role in modulating viral phenotype. However, the molecular mechanisms underlying the modulation remain poorly understood.

METHODOLOGY/PRINCIPAL FINDINGS

By combining computational and experimental approaches, we examined how V3 net charge could influence the phenotype of the gp120 interaction surface. Molecular dynamics simulations of the identical gp120 outer domain, carrying a V3 loop with net charge of +3 or +7, showed that the V3 change alone could induce global changes in fluctuation and conformation of the loops involved in binding to CD4, coreceptor and antibodies. A neutralization study using the V3 recombinant HIV-1 infectious clones showed that the virus carrying the gp120 with +3 V3, but not with +7 V3, was resistant to neutralization by anti-CD4 binding site monoclonal antibodies. An information entropy study shows that otherwise variable surface of the gp120 outer domain, such as V3 and a region around the CD4 binding loop, are less heterogeneous in the gp120 subpopulation with +3 V3.

CONCLUSIONS/SIGNIFICANCE

These results suggest that the HIV-1 gp120 V3 loop acts as an electrostatic modulator that influences the global structure and diversity of the interaction surface of the gp120 outer domain. Our findings will provide a novel structural basis to understand how HIV-1 adjusts relative replication fitness by V3 mutations.

The Fundamental Role of Flexibility on the Strength of Molecular Binding

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

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

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

Structural elucidation of the mechanistic basis of degeneracy in the primary humoral response

The mechanistic basis for efficient combating of the infinite range of foreign Ags by the limited repertoire of naive Abs expressed on primary B cell surfaces during their first encounter was addressed through elegantly designed crystallographic analyses. Resolution of the discrepancy arising from the limited number of possible germline Ab receptors on primary B cells for recognizing the unlimited pool of possible Ags has been attempted by invoking the degenerate recognition potential of the germline Abs. Structural analyses of germline mAb BBE6.12H3 in an Ag-free state, as well as bound to four different peptide Ags, established the correlation of its degenerate specificity with conformational versatility of the paratope. Six distinct paratope topologies observed for a single germline mAb provided a quantitative description of the primary Ag recognition repertoire at the tertiary structural level. Each of the four different peptide Ags was bound specifically to a distinct conformation of the paratope, which was also different from that of the Ag-free states of the same germline mAb. A minimal conserved motif in the pristine Ag-combining site essential for multispecificity and Ag binding-mediated change in the elbow angle of Fab was also discernible. It is proposed that the generation of a primary Ab repertoire involves large, yet finite, germline Ab clones, each capable of adopting discrete conformations, which in turn exhibit diverse binding modes.

A universal combinatorial design of antibody framework to graft distinct CDR sequences: a bioinformatics approach

Antibody (Ab) humanization is crucial to generate clinically relevant biologics from hybridoma-derived monoclonal antibodies (mAbs). In this study, we integrated antibody structural information from the Protein Data Bank with known back-to-mouse mutational data to build a universal consensus of framework positions (10 heavy and 7 light) critical for the preservation of the functional conformation of the Complimentarity Determining Region of antibodies. On the basis of FR consensus, we describe here a universal combinatorial library suitable for humanizing exogenous antibodies by CDR-grafting. The six CDRs of the murine anti-human EGFR Fab M225 were grafted onto a distinct (low FR sequence similarity to M225) human FR sequence that incorporates at the 17 FR consensus positions the permutations of the naturally observed amino acid diversities. Ten clones were selected from the combinatorial library expressing phage-displayed humanized M225 Fabs. Surprisingly, 2 of the 10 clones were found to bind EGFR with stronger affinity than M225. Cell-based assays demonstrated that the 10 selected clones retained epitope specificity by blocking EGFR phosphorylation and thus hindering cellular proliferation. Our results suggest that there is a universal and structurally rigid near-CDR set of FR positions that cooperatively support the binding conformation of CDRs.

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.

Highly specific off-target binding identified and eliminated during the humanization of an antibody against FGF receptor 4

Off-target binding can significantly affect the pharmacokinetics (PK), tissue distribution, efficacy and toxicity of a therapeutic antibody. Herein we describe the development of a humanized anti- fibroblast growth factor receptor 4 (FGFR4) antibody as a potential therapeutic for hepatocellular carcinoma (HCC). A chimeric anti FGFR4 monoclonal antibody (chLD1) was previously shown to block ligand binding and to inhibit FGFR4 mediated signaling as well as tumor growth in vivo. A humanized version of chLD1, hLD1.vB, had similar binding affinity and in vitro blocking activity, but it exhibited rapid clearance, poor target tissue biodistribution and limited efficacy when compared to chLD1 in a HUH7 human HCC xenograft mouse model. These problems were traced to instability of the molecule in rodent serum. Size exclusion high performance liquid chromatography, immunoprecipitation and mass spectral sequencing identified a specific interaction between hLD1.vB and mouse complement component 3 (C3). A PK study in C3 knock-out mice further confirmed this specific interaction. Subsequently, an affinity-matured variant derived from hLD1.vB (hLD1.v22), specifically selected for its lack of binding to mouse C3 was demonstrated to have a PK profile and in vivo efficacy similar to that of chLD1 in mice. Although reports of non-specific off-target binding have been observed for other antibodies, this represents the first report identifying a specific off-target interaction that affected disposition and biological activity. Screens developed to identify general non-specific interactions are likely to miss the rare and highly specific cross-reactivity identified in this study, thus highlighting the importance of animal models as a proxy for avoiding unexpected clinical outcomes.

Rational design of T cell receptors with enhanced sensitivity for antigen

Enhancing the affinity of therapeutic T cell receptors (TCR) without altering their specificity is a significant challenge for adoptive immunotherapy. Current efforts have primarily relied on empirical approaches. Here, we used structural analyses to identify a glycine-serine variation in the TCR that modulates antigen sensitivity. A G at position 107 within the CDR3β stalk is encoded within a single mouse and human TCR, TRBV13-2 and TRBV12-5 respectively. Most TCR bear a S107. The S hydroxymethyl side chain intercalates into the core of the CDR3β loop, stabilizing it. G107 TRBV possess a gap in their CDR3β where this S hydroxymethyl moiety would fit. We predicted based on modeling and molecular dynamics simulations that a G107S substitution would increase CDR3β stability and thereby augment receptor sensitivity. Experimentally, a G107S replacement led to an ∼10–1000 fold enhanced antigen sensitivity in 3 of 4 TRBV13-2⁺ TCR tested. Analysis of fine specificity indicated a preserved binding orientation. These results support the feasibility of developing high affinity antigen specific TCR for therapeutic purposes through the identification and manipulation of critical framework residues. They further indicate that amino acid variations within TRBV not directly involved in ligand contact can program TCR sensitivity, and suggest a role for CDR3 stability in this programming.

Enzyme catalysis from improved packing in their transition-state structures

The binding of ligands to proteins can be enhanced through improved packing within the proteins that may, or may not, occur with conformational change. Enzymes can similarly improve their catalytic magic through better packing in the transition state (TS) for reaction. In principle, the improved packing demands no more than the minute shortening of non-covalent interactions throughout much of the structure of the protein (positively cooperative binding). Improved protein packing can account for the remarkably high biotin/streptavidin affinity, and perhaps also for a major part of the catalytic function of hypoxanthine-guanine phosphoribosyltransferase and purine nucleoside phosphorylase (PNP). As successive NAD(+) molecules bind to the glyceraldehyde phosphate dehydrogenase tetramer, they do so with positively cooperative binding (using the term as applied in crystallization and protein folding) that decreases at each step. This binding is negatively cooperative in the usage stemming from Monod and co-workers.

Effects of somatic mutations on CDR loop flexibility during affinity maturation

Prior studies suggest that antibody affinity maturation is achieved, in part, via prearranging the CDRs for binding. The implication is that the entropy cost of binding is reduced and that this rigidification occurs as a consequence of somatic mutations during maturation. However, how these mutations modulate CDR flexibility is unclear. Here, molecular dynamics simulations captured CDR flexibility differences between four mature antibodies (7G12, AZ28, 28B4, and 48G7) and their germline predecessors. Analysis of their trajectories: (1) rationalized how mutations during affinity maturation restrict CDR motility, (2) captured the equilibrium between bound and unbound conformations for the H3 loop of unliganded 7G12, and (3) predicted a set of new mutations that, according to our simulations, should diminish binding by increasing flexibility.

Protein flexibility and conformational entropy in ligand design targeting the carbohydrate recognition domain of galectin-3

Rational drug design is predicated on knowledge of the three-dimensional structure of the protein-ligand complex and the thermodynamics of ligand binding. Despite the fundamental importance of both enthalpy and entropy in driving ligand binding, the role of conformational entropy is rarely addressed in drug design. In this work, we have probed the conformational entropy and its relative contribution to the free energy of ligand binding to the carbohydrate recognition domain of galectin-3. Using a combination of NMR spectroscopy, isothermal titration calorimetry, and X-ray crystallography, we characterized the binding of three ligands with dissociation constants ranging over 2 orders of magnitude. (15)N and (2)H spin relaxation measurements showed that the protein backbone and side chains respond to ligand binding by increased conformational fluctuations, on average, that differ among the three ligand-bound states. Variability in the response to ligand binding is prominent in the hydrophobic core, where a distal cluster of methyl groups becomes more rigid, whereas methyl groups closer to the binding site become more flexible. The results reveal an intricate interplay between structure and conformational fluctuations in the different complexes that fine-tunes the affinity. The estimated change in conformational entropy is comparable in magnitude to the binding enthalpy, demonstrating that it contributes favorably and significantly to ligand binding. We speculate that the relatively weak inherent protein-carbohydrate interactions and limited hydrophobic effect associated with oligosaccharide binding might have exerted evolutionary pressure on carbohydrate-binding proteins to increase the affinity by means of conformational entropy.

Using NMR to study fast dynamics in proteins: methods and applications

Proteins exist not as singular structures with precise coordinates, but rather as fluctuating bodies that move rapidly through an enormous number of conformational substates. These dynamics have important implications for understanding protein function and for structure-based drug design. NMR spectroscopy is particularly well suited to characterize the dynamics of proteins and other molecules in solution at atomic resolution. Here, NMR relaxation methods for characterizing thermal motions on the picosecond-nanosecond (ps-ns) timescale are reviewed. Motion on this timescale can be conveniently captured by the Lipari-Szabo order parameter, S², a bond-specific measure of restriction of motion. Approaches for determining order parameters are discussed, as are recent examples from the literature that link ps-ns dynamics with conformational entropy, allostery, and protein function in general.

The intrinsic protein flexibility of endogenous protease inhibitor TIMP-1 controls its binding interface and affects its function

Protein flexibility is thought to play key roles in numerous biological processes, including antibody affinity maturation, signal transduction, and enzyme catalysis, yet only limited information is available regarding the molecular details linking protein dynamics with function. A single point mutation at the distal site of the endogenous tissue inhibitor of metalloproteinase 1 (TIMP-1) enables this clinical target protein to tightly bind and inhibit membrane type 1 matrix metalloproteinase (MT1-MMP) by increasing only the association constant. The high-resolution X-ray structure of this complex determined at 2 A could not explain the mechanism of enhanced binding and pointed to a role for protein conformational dynamics. Molecular dynamics (MD) simulations reveal that the high-affinity TIMP-1 mutants exhibit significantly reduced binding interface flexibility and more stable hydrogen bond networks. This was accompanied by a redistribution of the ensemble of substrates to favorable binding conformations that fit the enzyme catalytic site. Apparently, the decrease in backbone flexibility led to a lower entropy cost upon formation of the complex. This work quantifies the effect of a single point mutation on the protein conformational dynamics and function of TIMP-1. Here we argue that controlling the intrinsic protein dynamics of MMP endogenous inhibitors may be utilized for rationalizing the design of selective novel protein inhibitors for this class of enzymes.

Class-specific restrictions define primase interactions with DNA template and replicative helicase

Bacterial primase is stimulated by replicative helicase to produce RNA primers that are essential for DNA replication. To identify mechanisms regulating primase activity, we characterized primase initiation specificity and interactions with the replicative helicase for gram-positive Firmicutes (Staphylococcus, Bacillus and Geobacillus) and gram-negative Proteobacteria (Escherichia, Yersinia and Pseudomonas). Contributions of the primase zinc-binding domain, RNA polymerase domain and helicase-binding domain on de novo primer synthesis were determined using mutated, truncated, chimeric and wild-type primases. Key residues in the β4 strand of the primase zinc-binding domain defined class-associated trinucleotide recognition and substitution of these amino acids transferred specificity across classes. A change in template recognition provided functional evidence for interaction in trans between the zinc-binding domain and RNA polymerase domain of two separate primases. Helicase binding to the primase C-terminal helicase-binding domain modulated RNA primer length in a species-specific manner and productive interactions paralleled genetic relatedness. Results demonstrated that primase template specificity is conserved within a bacterial class, whereas the primase-helicase interaction has co-evolved within each species.

Molecular description of flexibility in an antibody combining site

Mature antibodies (Abs) that are exquisitely specific for virtually any foreign molecule may be produced by affinity maturation of naïve (or germline) Abs. However, the finite number of germline Abs available suggests that, in contrast to mature Abs, germline Abs must be broadly polyspecific so that they are able to recognize a wide range of ligands. Thus, affinity maturation must play a role in mediating Ab specificity. One biophysical property that distinguishes polyspecificity from specificity is protein flexibility; a flexible combining site is able to adopt different conformations that recognize different foreign molecules (or antigens), while a rigid combining site is locked into a conformation that is specific for a given antigen. Recent studies (Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 8821-8826) have examined, at the atomic level, the structural properties that mediate changes in flexibility at four stages of affinity maturation in the 4-4-20 Ab. These studies employed molecular dynamics simulations to reveal a network of residue interactions that mediate the flexibility changes accompanying maturation. The flexibility of the Ab combining sites in these molecular systems was originally measured using three-pulse photon echo spectroscopy (3PEPS). The present investigation extends this work by providing a concrete link between structural properties of the Ab molecules and features of the spectroscopic measurements used to characterize their flexibility. Results obtained from the simulations are in good qualitative agreement with the experimental measurements and indicate that the spectroscopic signal is sensitive to protein dynamics distributed throughout the entire combining site. Thus, the simulations provide a molecular-level interpretation of the changes induced by affinity maturation of the Ab. The results suggest that 3PEPS spectroscopy in combination with molecular dynamics simulations can provide a detailed description of protein dynamics and, in this case, how it is evolved for biological function.

Computing free energies of protein conformations from explicit solvent simulations

We report a fully general technique addressing a long standing challenge of calculating conformational free energy differences between various states of a polymer chain from simulations using explicit solvent force fields. The main feature of our method is a special mapping variable, a path coordinate, which continuously connects two conformations. The path variable has been designed to preserve locality in the phase space near the path endpoints. We avoid the problem of sampling the unfolded states by creating an artificial confinement "tube" in the phase space that prevents the molecule from unfolding without affecting the calculation of the desired free energy difference. We applied our technique to compute the free energy difference between two native-like conformations of the small protein Trp-cage using the CHARMM force field with explicit solvent. We verified this result by comparing it with an independent, significantly more expensive calculation. Overall, the present study suggests that the new method of computing free energy differences between polymer chain conformations is accurate and highly computationally efficient.

Free energy simulations reveal a double mutant avian H5N1 virus hemagglutinin with altered receptor binding specificity

Historically, influenza pandemics have been triggered when an avian influenza virus or a human/avian reassorted virus acquires the ability to replicate efficiently and become transmissible in the human population. Most critically, the major surface glycoprotein hemagglutinin (HA) must adapt to the usage of human-like (alpha-2,6-linked) sialylated glycan receptors. Therefore, identification of mutations that can switch the currently circulating H5N1 HA receptor binding specificity from avian to human might provide leads to the emergence of pandemic H5N1 viruses. To define such mutations in the H5 subtype, here we provide a computational framework that combines molecular modeling with extensive free energy simulations. Our results show that the simulated binding affinities are in good agreement with currently available experimental data. Moreover, we predict that one double mutation (V135S and A138S) in HA significantly enhances alpha-2,6-linked receptor recognition by the H5 subtype. Our simulations indicate that this double mutation in H5N1 HA increases the binding affinity to alpha-2,6-linked sialic acid receptors by 2.6 +/- 0.7 kcal/mol per HA monomer that primarily arises from the electrostatic interactions. Further analyses reveal that introduction of this double mutation results in a conformational change in the receptor binding pocket of H5N1 HA. As a result, a major rearrangement occurs in the hydrogen-bonding network of HA with the human receptor, making the human receptor binding pattern of double mutant H5N1 HA surprisingly similar to that observed in human H1N1 HA. These large scale molecular simulations on single and double mutants thus provide new insights into our understanding toward human adaptation of the avian H5N1 virus.