Antibody evolution constrains conformational heterogeneity by tailoring protein dynamics

Evolutionary trajectory of receptor binding specificity and promiscuity of the spike protein of SARS-CoV-2

SARS-CoV-2 infects cells by attachment to its receptor-the angiotensin converting enzyme 2 (ACE2). Regardless of the wealth of structural data, little is known about the physicochemical mechanism of interactions of the viral spike (S) protein with ACE2 and how this mechanism has evolved during the pandemic. Here, we applied experimental and computational approaches to characterize the molecular interaction of S proteins from SARS-CoV-2 variants of concern (VOC). Data on kinetics, activation-, and equilibrium thermodynamics of binding of the receptor binding domain (RBD) from VOC with ACE2 as well as data from computational protein electrostatics revealed a profound remodeling of the physicochemical characteristics of the interaction during the evolution. Thus, as compared to RBDs from Wuhan strain and other VOC, Omicron RBD presented as a unique protein in terms of conformational dynamics and types of non-covalent forces driving the complex formation with ACE2. Viral evolution resulted in a restriction of the RBD structural dynamics, and a shift to a major role of polar forces for ACE2 binding. Further, we investigated how the reshaping of the physicochemical characteristics of interaction affects the binding specificity of S proteins. Data from various binding assays revealed that SARS-CoV-2 Wuhan and Omicron RBDs manifest capacity for promiscuous recognition of unrelated human proteins, but they harbor distinct reactivity patterns. These findings might contribute for mechanistic understanding of the viral tropism and capacity to evade immune responses during evolution.

A Trait-Based Approach to Predicting Viral Host-Range Evolvability

Predicting the evolution of virus host range has proven to be extremely difficult, in part because of the sheer diversity of viruses, each with unique biology and ecological interactions. We have not solved this problem, but to make the problem more tractable, we narrowed our focus to three traits intrinsic to all viruses that may play a role in host-range evolvability: mutation rate, recombination rate, and phenotypic heterogeneity. Although each trait should increase evolvability, they cannot do so unbounded because fitness trade-offs limit the ability of all three traits to maximize evolvability. By examining these constraints, we can begin to identify groups of viruses with suites of traits that make them especially concerning, as well as ecological and environmental conditions that might push evolution toward accelerating host-range expansion.

Higher Affinity Antibodies Bind With Lower Hydration and Flexibility in Large Scale Simulations

We have carried out a long-timescale simulation study on crystal structures of nine antibody-antigen pairs, in antigen-bound and antibody-only forms, using molecular dynamics with enhanced sampling and an explicit water model to explore interface conformation and hydration. By combining atomic level simulation and replica exchange to enable full protein flexibility, we find significant numbers of bridging water molecules at the antibody-antigen interface. Additionally, a higher proportion of interactions excluding bulk waters and a lower degree of antigen bound CDR conformational sampling are correlated with higher antibody affinity. The CDR sampling supports enthalpically driven antibody binding, as opposed to entropically driven, in that the difference between antigen bound and unbound conformations do not correlate with affinity. We thus propose that interactions with waters and CDR sampling are aspects of the interface that may moderate antibody-antigen binding, and that explicit hydration and CDR flexibility should be considered to improve antibody affinity prediction and computational design workflows.

Allostery in the dynamic coactivator domain KIX occurs through minor conformational micro-states

The coactivator KIX of CBP uses two binding surfaces to recognize multiple activators and exhibits allostery in ternary complex formation. Activator•coactivator interactions are central to transcriptional regulation, yet the microscopic origins of allostery in dynamic proteins like KIX are largely unknown. Here, we investigate the molecular recognition and allosteric manifestations involved in two KIX ternary systems c-Myb•KIX•MLL and pKID•KIX•MLL. Exploring the hypothesis that binary complex formation prepays an entropic cost for positive cooperativity, we utilize molecular dynamics simulations, side chain methyl order parameters, and differential scanning fluorimetry (DSF) to explore conformational entropy changes in KIX. The protein's configurational micro-states from structural clustering highlight the utility of protein plasticity in molecular recognition and allostery. We find that apo KIX occupies a wide distribution of lowly-populated configurational states. Each binding partner has its own suite of KIX states that it selects, building a model of molecular recognition fingerprints. Allostery is maximized with MLL pre-binding, which corresponds to the observation of a significant reduction in KIX micro-states observed when MLL binds. With all binding partners, the changes in KIX conformational entropy arise predominantly from changes in the most flexible loop. Likewise, we find that a small molecule and mutations allosterically inhibit/enhance activator binding by tuning loop dynamics, suggesting that loop-targeting chemical probes could be developed to alter KIX•activator interactions. Experimentally capturing KIX stabilization is challenging, particularly because of the disordered nature of particular activators. However, DSF melting curves allow for inference of relative entropic changes that occur across complexes, which we compare to our computed entropy changes using simulation methyl order parameters.

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

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

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

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

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

Redesigning an antibody H3 loop by virtual screening of a small library of human germline-derived sequences

The design of superior biologic therapeutics, including antibodies and engineered proteins, involves optimizing their specific ability to bind to disease-related molecular targets. Previously, we developed and applied the Assisted Design of Antibody and Protein Therapeutics (ADAPT) platform for virtual affinity maturation of antibodies (Vivcharuk et al. in PLoS One 12(7):e0181490, 10.1371/journal.pone.0181490, 2017). However, ADAPT is limited to point mutations of hot-spot residues in existing CDR loops. In this study, we explore the possibility of wholesale replacement of the entire H3 loop with no restriction to maintain the parental loop length. This complements other currently published studies that sample replacements for the CDR loops L1, L2, L3, H1 and H2. Given the immense sequence space theoretically available to H3, we focused on the virtual grafting of over 5000 human germline-derived H3 sequences from the IGMT/LIGM database increasing the diversity of the sequence space when compared to using crystalized H3 loop sequences. H3 loop conformations are generated and scored to identify optimized H3 sequences. Experimental testing of high-ranking H3 sequences grafted into the framework of the bH1 antibody against human VEGF-A led to the discovery of multiple hits, some of which had similar or better affinities relative to the parental antibody. In over 75% of the tested designs, the re-designed H3 loop contributed favorably to overall binding affinity. The hits also demonstrated good developability attributes such as high thermal stability and no aggregation. Crystal structures of select re-designed H3 variants were solved and indicated that although some deviations from predicted structures were seen in the more solvent accessible regions of the H3 loop, they did not significantly affect predicted affinity scores.

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

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

Local and Global Rigidification Upon Antibody Affinity Maturation

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

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.

Involvement of Local, Rapid Conformational Dynamics in Binding of Flexible Recognition Motifs

Flexible protein sequences populate ensembles of rapidly interconverting states differentiated by small-scale fluctuations; however, elucidating whether and how the ensembles determine function experimentally is challenged by the combined high spatial and temporal resolution needed to capture the states. We used carbon-deuterium (C-D) bond vibrations incorporated as infrared probes to characterize with residue-specific detail the heterogeneity of states adopted by proline-rich (PR) sequences and assess their involvement in recognition of Src homology 3 domains. The C-D absorption envelopes provided evidence for two or three sub-populations at all proline residues. The changes in the subpopulations induced by binding generally reflected recognition by conformational selection but depended on the residue and the state of the ligand to illuminate distinct mechanisms among the PR ligands. Notably, the spectral data indicate that greater adaptability among the states is associated with reduced recognition specificity and that perturbation to the ensemble populations contributes to differences in binding entropy. Broadly, the study quantifies rapidly interconverting ensembles with residue-specific detail and implicates them in function.

Optimizing antibody affinity and stability by the automated design of the variable light-heavy chain interfaces

Antibodies developed for research and clinical applications may exhibit suboptimal stability, expressibility, or affinity. Existing optimization strategies focus on surface mutations, whereas natural affinity maturation also introduces mutations in the antibody core, simultaneously improving stability and affinity. To systematically map the mutational tolerance of an antibody variable fragment (Fv), we performed yeast display and applied deep mutational scanning to an anti-lysozyme antibody and found that many of the affinity-enhancing mutations clustered at the variable light-heavy chain interface, within the antibody core. Rosetta design combined enhancing mutations, yielding a variant with tenfold higher affinity and substantially improved stability. To make this approach broadly accessible, we developed AbLIFT, an automated web server that designs multipoint core mutations to improve contacts between specific Fv light and heavy chains (http://AbLIFT.weizmann.ac.il). We applied AbLIFT to two unrelated antibodies targeting the human antigens VEGF and QSOX1. Strikingly, the designs improved stability, affinity, and expression yields. The results provide proof-of-principle for bypassing laborious cycles of antibody engineering through automated computational affinity and stability design.

Structure and Dynamics of Stacking Interactions in an Antibody Binding Site

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

Dynamic Aspects of the Immunoglobulin Structure

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

Deciphering evolution of immune recognition in antibodies

Background

Antibody, the primary effector molecule of the immune system, evolves after initial encounter with the antigen from a precursor form to a mature one to effectively deal with the antigen. Antibodies of a lineage diverge through antigen-directed isolated pathways of maturation to exhibit distinct recognition potential. In the context of evolution in immune recognition, diversity of antigen cannot be ignored. While there are reports on antibody lineage, structural perspective with respect to diverse recognition potential in a lineage has never been studied. Hence, it is crucial to evaluate how maturation leads to topological tailoring within a lineage enabling them to interact with significantly distinct antigens.

Results

A data-driven approach was undertaken for the study. Global experimental mouse and human antibody-antigen complex structures from PDB were compiled into a coherent database of germline-linked antibodies bound with distinct antigens. Structural analysis of all lineages showed variations in CDRs of both H and L chains. Observations of conformational adaptation made from analysis of static structures were further evaluated by characterizing dynamics of interaction in two lineages, mouse V_H1–84 and human V_H5–51. Sequence and structure analysis of the lineages explained that somatic mutations altered the geometries of individual antibodies with common structural constraints in some CDRs. Additionally, conformational landscape obtained from molecular dynamics simulations revealed that incoming pathogen led to further conformational divergence in the paratope (as observed across datasets) even while maintaining similar overall backbone topology. MM-GB/SA analysis showed binding energies to be in physiological range. Results of the study are coherent with experimental observations.

Conclusions

The findings of this study highlight basic structural principles shaping the molecular evolution of a lineage for significantly diverse antigens. Antibodies of a lineage follow different developmental pathways while preserving the imprint of the germline. From the study, it can be generalized that structural diversification of the paratope is an outcome of natural selection of a conformation from an available ensemble, which is further optimized for antigen interaction. The study establishes that starting from a common lineage, antibodies can mature to recognize a wide range of antigens. This hypothesis can be further tested and validated experimentally.

Electronic supplementary material

The online version of this article (10.1186/s12900-018-0096-1) contains supplementary material, which is available to authorized users.

Allosteric Effects between the Antibody Constant and Variable Regions: A Study of IgA Fc Mutations on Antigen Binding

Therapeutic antibodies have shifted the paradigm of disease treatments from small molecules to biologics, especially in cancer therapy. Despite the increasing number of antibody candidates, much remains unknown about the antibody and how its various regions interact. Recent findings showed that the antibody constant region can govern localization effects that are useful in reducing side effects due to systemic circulation by the commonly used IgG isotypes. Given their localized mucosal effects, IgA antibodies are increasingly promising therapeutic biologics. While the antibody Fc effector cell activity has been a focus point, recent research showed that the Fc could also influence antigen binding, challenging the conventional idea of region-specific antibody functions. To investigate this, we analysed the IgA antibody constant region and its distal effects on the antigen binding regions using recombinant Pertuzumab IgA1 and IgA2 variants. We found that mutations in the C-region reduced Her2 binding experimentally, and computational structural analysis showed that allosteric communications were highly dependent on the antibody hinge, providing strong evidence that we should consider antibodies as whole proteins rather than a sum of functional regions.

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

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

The Unusual Genetics and Biochemistry of Bovine Immunoglobulins

Antibodies are the key circulating molecules that have evolved to fight infection by the adaptive immune system of vertebrates. Typical antibodies of most species contain six complementarity-determining regions (CDRs), where the third CDR of the heavy chain (CDR H3) has the greatest diversity and often makes the most significant contact with antigen. Generally, the process of V(D)J recombination produces a vast repertoire of antibodies; multiple V, D, and J gene segments recombine with additional junctional diversity at the V-D and D-J joints, and additional combinatorial possibilities occur through heavy- and light-chain pairing. Despite these processes, the overall structure of the resulting antibody is largely conserved, and binding to antigen occurs predominantly through the CDR loops of the immunoglobulin V domains. Bovines have deviated from this general paradigm by having few VH regions and thus little germline combinatorial diversity, but their antibodies contain long CDR H3 regions, with substantial diversity generated through somatic hypermutation. A subset of the repertoire comprises antibodies with ultralong CDR H3s, which can reach over 70 amino acids in length. Structurally, these unusual antibodies form a β-ribbon "stalk" and disulfide-bonded "knob" that protrude far from the antibody surface. These long CDR H3s allow cows to mount a particularly robust immune response when immunized with viral antigens, particularly to broadly neutralizing epitopes on a stabilized HIV gp140 trimer, which has been a challenge for other species. The unusual genetics and structural biology of cows provide for a unique paradigm for creation of immune diversity and could enable generation of antibodies against especially challenging targets and epitopes.

Antibody Affinity Maturation by Computational Design

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

Molecular rigidity and enthalpy-entropy compensation in DNA melting

Enthalpy-entropy compensation (EEC) is observed in diverse molecular binding processes of importance to living systems and manufacturing applications, but this widely occurring phenomenon is not sufficiently understood from a molecular physics standpoint. To gain insight into this fundamental problem, we focus on the melting of double-stranded DNA (dsDNA) since measurements exhibiting EEC are extensive for nucleic acid complexes and existing coarse-grained models of DNA allow us to explore the influence of changes in molecular parameters on the energetic parameters by using molecular dynamics simulations. Previous experimental and computational studies have indicated a correlation between EEC and changes in molecular rigidity in certain binding-unbinding processes, and, correspondingly, we estimate measures of DNA molecular rigidity under a wide range of conditions, along with resultant changes in the enthalpy and entropy of binding. In particular, we consider variations in dsDNA rigidity that arise from changes of intrinsic molecular rigidity such as varying the associative interaction strength between the DNA bases, the length of the DNA chains, and the bending stiffness of the individual DNA chains. We also consider extrinsic changes of molecular rigidity arising from the addition of polymer additives and geometrical confinement of DNA between parallel plates. All our computations confirm EEC and indicate that this phenomenon is indeed highly correlated with changes in molecular rigidity. However, two distinct patterns relating to how DNA rigidity influences the entropy of association emerge from our analysis. Increasing the intrinsic DNA rigidity increases the entropy of binding, but increases in molecular rigidity from external constraints decreases the entropy of binding. EEC arises in numerous synthetic and biological binding processes and we suggest that changes in molecular rigidity might provide a common origin of this ubiquitous phenomenon in the mutual binding and unbinding of complex molecules.

Therapeutic Antibodies to Ganglioside GD2 Evolved from Highly Selective Germline Antibodies

Antibodies play a crucial role in host defense and are indispensable research tools, diagnostics, and therapeutics. Antibody generation involves binding of genomically encoded germline antibodies followed by somatic hypermutation and in vivo selection to obtain antibodies with high affinity and selectivity. Understanding this process is critical for developing monoclonal antibodies, designing effective vaccines, and understanding autoantibody formation. Prior studies have found that antibodies to haptens, peptides, and proteins evolve from polyspecific germline antibodies. The immunological evolution of antibodies to mammalian glycans has not been studied. Using glycan microarrays, protein microarrays, cell binding studies, and molecular modeling, we demonstrate that therapeutic antibodies to the tumor-associated ganglioside GD2 evolved from highly specific germline precursors. The results have important implications for developing vaccines and monoclonal antibodies that target carbohydrate antigens. In addition, they demonstrate an alternative pathway for antibody evolution within the immune system that is distinct from the polyspecific germline pathway.

Applications of Normal Mode Analysis Methods in Computational Protein Design

Recent advances in coarse-grained normal mode analysis methods make possible the large-scale prediction of the effect of mutations on protein stability and dynamics as well as the generation of biologically relevant conformational ensembles. Given the interplay between flexibility and enzymatic activity, the combined analysis of stability and dynamics using the Elastic Network Contact Model (ENCoM) method has ample applications in protein engineering in industrial and medical applications such as in computational antibody design. Here, we present a detailed tutorial on how to perform such calculations using ENCoM.

Somatic Hypermutation-Induced Changes in the Structure and Dynamics of HIV-1 Broadly Neutralizing Antibodies

Antibody somatic hypermutation (SHM) and affinity maturation enhance antigen recognition by modifying antibody paratope structure to improve its complementarity with the target epitope. SHM-induced changes in paratope dynamics may also contribute to antibody maturation, but direct evidence of this is limited. Here, we examine two classes of HIV-1 broadly neutralizing antibodies (bNAbs) for SHM-induced changes in structure and dynamics, and delineate the effects of these changes on interactions with the HIV-1 envelope glycoprotein (Env). In combination with new and existing structures of unmutated and affinity matured antibody Fab fragments, we used hydrogen/deuterium exchange with mass spectrometry to directly measure Fab structural dynamics. Changes in antibody structure and dynamics were positioned to improve complementarity with Env, with changes in dynamics primarily observed at the paratope peripheries. We conclude that SHM optimizes paratope complementarity to conserved HIV-1 epitopes and restricts the mobility of paratope-peripheral residues to minimize clashes with variable features on HIV-1 Env.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Resolution of Site-Specific Conformational Heterogeneity in Proline-Rich Molecular Recognition by Src Homology 3 Domains

Conformational heterogeneity and dynamics are increasingly evoked in models of protein molecular recognition but are challenging to experimentally characterize. Here we combine the inherent temporal resolution of infrared (IR) spectroscopy with the spatial resolution afforded by selective incorporation of carbon-deuterium (C-D) bonds, which provide frequency-resolved absorptions within a protein IR spectrum, to characterize the molecular recognition of the Src homology 3 (SH3) domain of the yeast protein Sho1 with its cognate proline-rich (PR) sequence of Pbs2. The IR absorptions of C-D bonds introduced at residues along a peptide of the Pbs2 PR sequence report on the changes in the local environments upon binding to the SH3 domain. Interestingly, upon forming the complex the IR spectra of the peptides labeled with C-D bonds at either of the two conserved prolines of the PXXP consensus recognition sequence show more absorptions than there are C-D bonds, providing evidence for the population of multiple states. In contrast, the NMR spectra of the peptides labeled with (13)C at the same residues show only single resonances, indicating rapid interconversion on the NMR time scale. Thus, the data suggest that the SH3 domain recognizes its cognate peptide with a component of induced fit molecular recognition involving the adoption of multiples states, which have previously gone undetected due to interconversion between the populated states that is too fast to resolve using conventional methods.

How to assess the binding strength of antibodies elicited by vaccination against HIV and other viruses

Vaccines that protect against viral infections generally induce neutralizing antibodies. When vaccines are evaluated, the need arises to assess the affinity maturation of the antibody responses. Binding titers of polyclonal sera depend not only on the affinities of the constituent antibodies but also on their individual concentrations, which are difficult to ascertain. Therefore an assay based on chaotrope disruption of antibody-antigen complexes was designed for measuring binding strength. This assay works well with many viral antigens but gives differential results depending on the conformational dependence of epitopes on complex antigens such as the envelope glycoprotein of HIV-1. Kinetic binding assays might offer alternatives, since they can measure average off-rate constants for polyclonal antibodies in a serum. Here, potentials and fallacies of these techniques are discussed.

Using THz Spectroscopy, Evolutionary Network Analysis Methods, and MD Simulation to Map the Evolution of Allosteric Communication Pathways in c-Type Lysozymes

It is now widely accepted that protein function is intimately tied with the navigation of energy landscapes. In this framework, a protein sequence is not described by a distinct structure but rather by an ensemble of conformations. And it is through this ensemble that evolution is able to modify a protein's function by altering its landscape. Hence, the evolution of protein functions involves selective pressures that adjust the sampling of the conformational states. In this work, we focus on elucidating the evolutionary pathway that shaped the function of individual proteins that make-up the mammalian c-type lysozyme subfamily. Using both experimental and computational methods, we map out specific intermolecular interactions that direct the sampling of conformational states and accordingly, also underlie shifts in the landscape that are directly connected with the formation of novel protein functions. By contrasting three representative proteins in the family we identify molecular mechanisms that are associated with the selectivity of enhanced antimicrobial properties and consequently, divergent protein function. Namely, we link the extent of localized fluctuations involving the loop separating helices A and B with shifts in the equilibrium of the ensemble of conformational states that mediate interdomain coupling and concurrently moderate substrate binding affinity. This work reveals unique insights into the molecular level mechanisms that promote the progression of interactions that connect the immune response to infection with the nutritional properties of lactation, while also providing a deeper understanding about how evolving energy landscapes may define present-day protein function.

Methyl Relaxation Measurements Reveal Patterns of Fast Dynamics in a Viral RNA-Directed RNA Polymerase

Molecular dynamics (MD) simulations combined with biochemical studies have suggested the presence of long-range networks of functionally relevant conformational flexibility on the nanosecond time scale in single-subunit RNA polymerases in many RNA viruses. However, experimental verification of these dynamics at a sufficient level of detail has been lacking. Here we describe the fast, picosecond to nanosecond dynamics of an archetypal viral RNA-directed RNA polymerase (RdRp), the 75 kDa P2 protein from cystovirus ϕ12, using analyses of (1)H-(1)H dipole-dipole cross-correlated relaxation at the methyl positions of Ile (δ1), Leu, Val, and Met residues. Our results, which represent the most detailed experimental characterization of fast dynamics in a viral RdRp until date, reveal a highly connected dynamic network as predicted by MD simulations of related systems. Our results suggest that the entry portals for template RNA and substrate NTPs are relatively disordered, while conserved motifs involved in metal binding, nucleotide selection, and catalysis display greater rigidity. Perturbations at the active site through metal binding or functional mutation affect dynamics not only in the immediate vicinity but also at remote regions. Comparison with the limited experimental and extensive functional and in silico results available for homologous systems suggests conservation of the overall pattern of dynamics in viral RdRps.

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

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

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

Structural analysis of the unmutated ancestor of the HIV-1 envelope V2 region antibody CH58 isolated from an RV144 vaccine efficacy trial vaccinee

Human monoclonal antibody CH58 isolated from an RV144 vaccinee binds at Lys169 of the HIV-1 Env gp120 V2 region, a site of vaccine-induced immune pressure. CH58 neutralizes HIV-1 CRF_01 AE strain 92TH023 and mediates ADCC against CD4 + T cell targets infected with CRF_01 AE tier 2 virus. CH58 and other antibodies that bind to a gp120 V2 epitope have a second light chain complementarity determining region (LCDR2) bearing a glutamic acid, aspartic acid (ED) motif involved in forming salt bridges with polar, basic side amino acid side chains in V2. In an effort to learn how V2 responses develop, we determined the crystal structures of the CH58-UA antibody unliganded and bound to V2 peptide. The structures showed an LCDR2 structurally pre-conformed from germline to interact with V2 residue Lys169. LCDR3 was subject to conformational selection through the affinity maturation process. Kinetic analyses demonstrate that only a few contacts were responsible for a 2000-fold increase in K_D through maturation, and this effect was predominantly due to an improvement in off-rate. This study shows that preconformation and preconfiguration can work in concert to produce antibodies with desired immunogenic properties.

•

With only 2-3% mutation from germline, the HIV-1 antibody CH58 developed neutralizing and ADCC capabilities.

•

The LCDR2 Glu–Asp motif of the RV144 antibody CH58 is pre-conformed from germline to interact with the gp120 V2 loop.

•

Affinity and neutralization gains resulted from tuning local interactions rather than gross sequence or structure changes.

•

Structural analyses show the second light chain complementarity determining region Glu–Asp motif of the CH58 antibody isolated from an RV144 vaccinee is optimally pre-conformed from germline to interact with the gp120 V2 loop. The increased binding affinity and neutralization capacity of the mature antibody compared to its germline precursor were achieved with only 2–3% mutation from germline, and the fact that these gains appeared to be a result of the tuning of local interactions rather than gross sequential or conformational changes provides hope that a rational immunogen design for HIV-1 treatment may become a reality.

Adaptive mutations alter antibody structure and dynamics during affinity maturation

While adaptive mutations can bestow new functions on proteins via the introduction or optimization of reactive centers, or other structural changes, a role for the optimization of protein dynamics also seems likely but has been more difficult to evaluate. Antibody (Ab) affinity maturation is an example of adaptive evolution wherein the adaptive mutations may be identified and Abs may be raised to specific targets that facilitate the characterization of protein dynamics. Here, we report the characterization of three affinity matured Abs that evolved from a common germline precursor to bind the chromophoric antigen (Ag), 8-methoxypyrene-1,3,6-trisulfonate (MPTS). In addition to characterizing the sequence, molecular recognition, and structure of each Ab, we characterized the dynamics of each complex by determining their mechanical response to an applied force via three-pulse photon echo peak shift (3PEPS) spectroscopy and deconvoluting the response into elastic, anelastic, and plastic components. We find that for one Ab, affinity maturation was accomplished via the introduction of a single functional group that mediates a direct contact with MPTS and results in a complex with little anelasticity or plasticity. In the other two cases, more mutations were introduced but none directly contact MPTS, and while their effects on structure are subtle, their effects on anelasticity and plasticity are significant, with the level of plasticity correlated with specificity, suggesting that the optimization of protein dynamics may have contributed to affinity maturation. A similar optimization of structure and dynamics may contribute to the evolution of other proteins.

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

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

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.

Design, expression and characterization of a single chain anti-CD20 antibody; a germline humanized antibody derived from Rituximab

CD20 is a B cell lineage specific surface antigen involved in various B cell malignancies. So far, several murine and chimeric antibodies have been produced against this antigen among which Rituximab is a commercially approved antibody widely used in treatment of cancers associated with CD20 overexpression. The current study reports the production and characterization of a humanized single chain version of Rituximab through CDR grafting method. For either heavy or light chain variable domains, a human antibody with the highest sequence homology to Rituximab was selected from human germline sequences and used as framework donors. Vernier zone residues in framework regions were replaced with those of Rituximab to retain the antigen binding affinity of parental antibody. The reactivity of humanized single chain antibody with CD20 was examined by ELISA and dot blot assays. The ability of antibody to suppress the growth of CD20 overexpressing Raji cells was tested by MTT assay. Analysis of reactivity with CD20 antigen revealed that the humanized single chain antibody reacted to the target antigen with high affinity. Proliferation inhibition assay showed that humanized scFv could suppress the proliferation of Raji cells efficiently in a dose-dependent manner. This successful production of a humanized scFv with the ability to inhibit growth of CD20-expressing cancer cell may provide a promising alternative strategy for CD20 targeted therapy.

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.

Carbon-deuterium bonds as non-perturbative infrared probes of protein dynamics, electrostatics, heterogeneity, and folding

Vibrational spectroscopy is uniquely able to characterize protein dynamics and microenvironmental heterogeneity because it possesses an inherently high temporal resolution and employs probes of ultimately high structural resolution-the bonds themselves. The use of carbon-deuterium (C-D) bonds as vibrational labels circumvents the spectral congestion that otherwise precludes the use of vibrational spectroscopy to proteins and makes the observation of single vibrations within a protein possible while being wholly non-perturbative. Thus, C-D probes can be used to site-specifically characterize conformational heterogeneity and thermodynamic stability. C-D probes are also uniquely useful in characterizing the electrostatic microenvironment experienced by a specific residue side chain or backbone due to its effect on the C-D absorption frequency. In this chapter we describe the experimental procedures required to use C-D bonds and FT IR spectroscopy to characterize protein dynamics, structural and electrostatic heterogeneity, ligand binding, and folding.

Redistribution of flexibility in stabilizing antibody fragment mutants follows Le Châtelier's principle

Le Châtelier's principle is the cornerstone of our understanding of chemical equilibria. When a system at equilibrium undergoes a change in concentration or thermodynamic state (i.e., temperature, pressure, etc.), La Châtelier's principle states that an equilibrium shift will occur to offset the perturbation and a new equilibrium is established. We demonstrate that the effects of stabilizing mutations on the rigidity ⇔ flexibility equilibrium within the native state ensemble manifest themselves through enthalpy-entropy compensation as the protein structure adjusts to restore the global balance between the two. Specifically, we characterize the effects of mutation to single chain fragments of the anti-lymphotoxin-β receptor antibody using a computational Distance Constraint Model. Statistically significant changes in the distribution of both rigidity and flexibility within the molecular structure is typically observed, where the local perturbations often lead to distal shifts in flexibility and rigidity profiles. Nevertheless, the net gain or loss in flexibility of individual mutants can be skewed. Despite all mutants being exclusively stabilizing in this dataset, increased flexibility is slightly more common than increased rigidity. Mechanistically the redistribution of flexibility is largely controlled by changes in the H-bond network. For example, a stabilizing mutation can induce an increase in rigidity locally due to the formation of new H-bonds, and simultaneously break H-bonds elsewhere leading to increased flexibility distant from the mutation site via Le Châtelier. Increased flexibility within the VH β4/β5 loop is a noteworthy illustration of this long-range effect.

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.

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.

Directed polymerase evolution

Polymerases evolved in nature to synthesize DNA and RNA, and they underlie the storage and flow of genetic information in all cells. The availability of these enzymes for use at the bench has driven a revolution in biotechnology and medicinal research; however, polymerases did not evolve to function efficiently under the conditions required for some applications and their high substrate fidelity precludes their use for most applications that involve modified substrates. To circumvent these limitations, researchers have turned to directed evolution to tailor the properties and/or substrate repertoire of polymerases for different applications, and several systems have been developed for this purpose. These systems draw on different methods of creating a pool of randomly mutated polymerases and are differentiated by the process used to isolate the most fit members. A variety of polymerases have been evolved, providing new or improved functionality, as well as interesting new insight into the factors governing activity.

Experimental characterization of electrostatic and conformational heterogeneity in an SH3 domain

Electrostatic and conformational heterogeneity make central contributions to protein function, but their experimental characterization requires a combination of spatial and temporal resolution that is challenging to achieve. Src homology 3 (SH3) domains mediate protein-protein interactions, and NMR studies have demonstrated that most possess conformational heterogeneity, which could be critical for their function. Here, we use the IR absorptions of carbon-deuterium (C-D) bonds site-selectively incorporated throughout the N-terminal SH3 domain from the murine adapter protein Crk-II to characterize its different microenvironments with high spatial and temporal resolution. The C-D absorptions are only differentiated in the folded state of the protein where they show evidence of significant environmental heterogeneity. However, the spectra of the folded state are independent of temperature, and upon thermal denaturation the protein undergoes a single, global unfolding transition. While some evidence of conformational heterogeneity is found within the peptide backbone, the majority of the environmental heterogeneity appears to result from electrostatics.