Site-specific analysis of mutational effects in proteins

Epitopes, paratopes, and other topes 30 years on: Understanding what we are talking about

The question of which protein antigens, such as HLA class I or class II molecules, will bind, and how well, to a given antibody is often assumed to depend exclusively on the details of protein surface structure. These structures are usually based on static models resulting from X-ray crystallography. While these notions are useful, the ultimate causal factors determining how well a given antigen binds a given antibody are based in thermodynamics and can include atomic mobility and the time-varying conformations of proteins. In this article, fundamental biophysical principles of antibody-antigen interaction are discussed, concepts critical for a deeper understanding of the pertinent molecular phenomena are highlighted, and common misunderstandings are identified and debunked.

Role of sequence and position of the cleavage sites in prothrombin activation

In the penultimate step of the coagulation cascade, the multidomain vitamin-K-dependent zymogen prothrombin is converted to thrombin by the prothrombinase complex composed of factor Xa, cofactor Va, and phospholipids. Activation of prothrombin requires cleavage at two residues, R271 and R320, along two possible pathways generating either the intermediate prethrombin-2 (following initial cleavage at R271) or meizothrombin (following initial cleavage at R320). The former pathway is preferred in the absence of and the latter in the presence of cofactor Va. Several mechanisms have been proposed to explain this preference, but the role of the sequence and position of the sites of cleavage has not been thoroughly investigated. In this study, we engineered constructs where the sequences ²⁶¹DEDSDRAIEGRTATSEYQT²⁷⁹ and ³¹⁰RELLESYIDGRIVEGSDAE³²⁸ were swapped between the R271 and R320 sites. We found that in the absence of cofactor Va, the wild-type sequence at the R271 site is cleaved preferentially regardless of its position at the R271 or R320 site, whereas in the presence of cofactor Va, the R320 site is cleaved preferentially regardless of its sequence. Additional single-molecule FRET measurements revealed that the environment of R271 changes significantly upon cleavage at R320 due to the conformational transition from the closed form of prothrombin to the open form of meizothrombin. Detailed kinetics of cleavage at the R271 site were monitored by a newly developed assay based on loss of FRET. These findings show how sequence and position of the cleavage sites at R271 and R320 dictate the preferred pathway of prothrombin activation.

Conformational Dynamics and Cooperativity Drive the Specificity of a Protein-Ligand Interaction

Molecular recognition is critical for the fidelity of signal transduction in biology. Conversely, the disruption of protein-protein interactions can lead to disease. Thus, comprehension of the molecular determinants of specificity is essential for understanding normal biological signaling processes and for the development of precise therapeutics. Although high-resolution structures have provided atomic details of molecular interactions, much less is known about the influence of cooperativity and conformational dynamics. Here, we used the Tiam2 PSD-95/Dlg/ZO-1 (PDZ) domain and a quadruple mutant (QM), engineered by swapping the identity of four residues important for specificity in the Tiam1 PDZ into the Tiam2 PDZ domain, as a model system to investigate the role of cooperativity and dynamics in PDZ ligand specificity. Surprisingly, equilibrium binding experiments found that the ligand specificity of the Tiam2 QM was switched to that of the Tiam1 PDZ. NMR-based studies indicated that Tiam2 QM PDZ, but not other mutants, had extensive microsecond to millisecond motions distributed throughout the entire domain suggesting structural cooperativity between the mutated residues. Thermodynamic analyses revealed energetic cooperativity between residues in distinct specificity subpockets that was dependent upon the identity of the ligand, indicating a context-dependent binding mechanism. Finally, isothermal titration calorimetry experiments showed distinct entropic signatures along the mutational trajectory from the Tiam2 wild-type to the QM PDZ domain. Collectively, our studies provide unique insights into how structure, conformational dynamics, and thermodynamics combine to modulate ligand-binding specificity and have implications for the evolution, regulation, and design of protein-ligand interactions.

Investigating Substitutions in Antibody-Antigen Complexes Using Molecular Dynamics: A Case Study with Broad-spectrum, Influenza A Antibodies

In studying the binding of host antibodies to the surface antigens of pathogens, the structural and functional characterization of antibody–antigen complexes by X-ray crystallography and binding assay is important. However, the characterization requires experiments that are typically time consuming and expensive: thus, many antibody–antigen complexes are under-characterized. For vaccine development and disease surveillance, it is often vital to assess the impact of amino acid substitutions on antibody binding. For example, are there antibody substitutions capable of improving binding without a loss of breadth, or antigen substitutions that lead to antigenic escape? The questions cannot be answered reliably from sequence variation alone, exhaustive substitution assays are usually impractical, and alanine scans provide at best an incomplete identification of the critical residue–residue interactions. Here, we show that, given an initial structure of an antibody bound to an antigen, molecular dynamics simulations using the energy method molecular mechanics with Generalized Born surface area (MM/GBSA) can model the impact of single amino acid substitutions on antibody–antigen binding energy. We apply the technique to three broad-spectrum antibodies to influenza A hemagglutinin and examine both previously characterized and novel variant strains observed in the human population that may give rise to antigenic escape. We find that in some cases the impact of a substitution is local, while in others it causes a reorientation of the antibody with wide-ranging impact on residue–residue interactions: this explains, in part, why the change in chemical properties of a residue can be, on its own, a poor predictor of overall change in binding energy. Our estimates are in good agreement with experimental results—indeed, they approximate the degree of agreement between different experimental techniques. Simulations were performed on commodity computer hardware; hence, this approach has the potential to be widely adopted by those undertaking infectious disease research. Novel aspects of this research include the application of MM/GBSA to investigate binding between broadly binding antibodies and a viral glycoprotein; the development of an approach for visualizing substrate–ligand interactions; and the use of experimental assay data to rescale our predictions, allowing us to make inferences about absolute, as well as relative, changes in binding energy.

A compare-and-contrast NMR dynamics study of two related RRMs: U1A and SNF

The U1A/U2B″/SNF family of small nuclear ribonucleoproteins uses a phylogenetically conserved RNA recognition motif (RRM1) to bind RNA stemloops in U1 and/or U2 small nuclear RNA (snRNA). RRMs are characterized by their α/β sandwich topology, and these RRMs use their β-sheet as the RNA binding surface. Unique to this RRM family is the tyrosine-glutamine-phenylalanine (YQF) triad of solvent-exposed residues that are displayed on the β-sheet surface; the aromatic residues form a platform for RNA nucleobases to stack. U1A, U2B″, and SNF have very different patterns of RNA binding affinity and specificity, however, so here we ask how YQF in Drosophila SNF RRM1 contributes to RNA binding, as well as to domain stability and dynamics. Thermodynamic double-mutant cycles using tyrosine and phenylalanine substitutions probe the communication between those two residues in the free and bound states of the RRM. NMR experiments follow corresponding changes in the glutamine side-chain amide in both U1A and SNF, providing a physical picture of the RRM1 β-sheet surface. NMR relaxation and dispersion experiments compare fast (picosecond to nanosecond) and intermediate (microsecond-to-millisecond) dynamics of U1A and SNF RRM1. We conclude that there is a network of amino acid interactions involving Tyr-Gln-Phe in both SNF and U1A RRM1, but whereas mutations of the Tyr-Gln-Phe triad result in small local responses in U1A, they produce extensive microsecond-to-millisecond global motions throughout SNF that alter the conformational states of the RRM.

Protein thermodynamics and the cognitive ecology of biomedicine

Assessments of scientific contributions critically influence decisions about grant funding and academic promotion. Unfortunately, there is a tendency for more junior and less assertive individuals to receive less credit than deserved. Acknowledgement of the complexity of relationships among researchers and the different modes of contributing to scientific progress could improve this situation. The thermodynamics of ligand binding is arguably among the most quantitative and empirically validated theoretical frameworks that permit precise apportionment of "credit" to multiple interacting entities that collectively account for a biologically relevant outcome, in this case, receptor-ligand complex formation. The process for assigning credit for research advances to individual researchers might benefit from emulating this thermodynamic thought process by recognizing that contributions of equal quantitative significance can be of different types and can originate through indirect effects. If the hypothesis that some categories of research contribution are frequently under-valued is correct, then calling attention to this state of affairs and providing an alternative way to conceptualize the task of credit attribution has the potential to begin altering the status quo. A beginning step to improving our credit attribution process would be the empirical investigation of accounts of contributions to particular scientific advances from all research team members.

Why Ser and not Thr brokers catalysis in the trypsin fold

Although Thr is equally represented as Ser in the human genome and as a nucleophile is as good as Ser, it is never found in the active site of the large family of trypsin-like proteases that utilize the Asp/His/Ser triad. The molecular basis of the preference of Ser over Thr in the trypsin fold was investigated with X-ray structures of the thrombin mutant S195T free and bound to an irreversible active site inhibitor. In the free form, the methyl group of T195 is oriented toward the incoming substrate in a conformation seemingly incompatible with productive binding. In the bound form, the side chain of T195 is reoriented for efficient substrate acylation without causing steric clash within the active site. Rapid kinetics prove that this change is due to selection of an active conformation from a preexisting ensemble of reactive and unreactive rotamers whose relative distribution determines the level of activity of the protease. Consistent with these observations, the S195T substitution is associated with a weak yet finite activity that allows identification of an unanticipated important role for S195 as the end point of allosteric transduction in the trypsin fold. The S195T mutation abrogates the Na(+)-dependent enhancement of catalytic activity in thrombin, activated protein C, and factor Xa and significantly weakens the physiologically important allosteric effects of thrombomodulin on thrombin and of cofactor Va on factor Xa. The evolutionary selection of Ser over Thr in trypsin-like proteases was therefore driven by the need for high catalytic activity and efficient allosteric regulation.

Molecular mechanisms of prolactin and its receptor

Prolactin and the prolactin receptors are members of a family of hormone/receptor pairs which include GH, erythropoietin, and other ligand/receptor pairs. The mechanisms of these ligand/receptor pairs have broad similarities, including general structures, ligand/receptor stoichiometries, and activation of several common signaling pathways. But significant variations in the structural and mechanistic details are present among these hormones and their type 1 receptors. The prolactin receptor is particularly interesting because it can be activated by three sequence-diverse human hormones: prolactin, GH, and placental lactogen. This system offers a unique opportunity to compare the detailed molecular mechanisms of these related hormone/receptor pairs. This review critically evaluates selected literature that informs these mechanisms, compares the mechanisms of the three lactogenic hormones, compares the mechanism with those of other class 1 ligand/receptor pairs, and identifies information that will be required to resolve mechanistic ambiguities. The literature describes distinct mechanistic differences between the three lactogenic hormones and their interaction with the prolactin receptor and describes more significant differences between the mechanisms by which other related ligands interact with and activate their receptors.

Cooperative tertiary interaction network guides RNA folding

Noncoding RNAs form unique 3D structures, which perform many regulatory functions. To understand how RNAs fold uniquely despite a small number of tertiary interaction motifs, we mutated the major tertiary interactions in a group I ribozyme by single-base substitutions. The resulting perturbations to the folding energy landscape were measured using SAXS, ribozyme activity, hydroxyl radical footprinting, and native PAGE. Double- and triple-mutant cycles show that most tertiary interactions have a small effect on the stability of the native state. Instead, the formation of core and peripheral structural motifs is cooperatively linked in near-native folding intermediates, and this cooperativity depends on the native helix orientation. The emergence of a cooperative interaction network at an early stage of folding suppresses nonnative structures and guides the search for the native state. We suggest that cooperativity in noncoding RNAs arose from natural selection of architectures conducive to forming a unique, stable fold.

Binding techniques to study the allosteric energy cycle

Thermodynamic principles of cooperativity and allostery have long been used as a starting point to begin understanding the interplay between ligand binding events. Understanding the nature of allosteric effects requires an experimental technique that can be used to quantify ligand binding energies and simultaneously give experimental insights into the conformational dynamics at play upon ligand binding. CD spectroscopy provides macroscopic information about the relative secondary and tertiary structures present in a protein. Here, we use this spectroscopic technique with thermal shift assays wherein ligand binding constants can be quantified based on their stabilizing effect against thermally induced protein denaturation. Binding constants for two ligands are used to determine a pairwise coupling free energy which defines the shared energy that favors or opposes binding of the second ligand binding event in an allosteric system. In CD-based thermal shift assays, temperature is the driving force for protein unfolding and can also influence protein conformational dynamics present in the unbound protein or ligand-bound proteins. Dihydrofolate reductase (DHFR) and glutamate dehydrogenase (GDH) are proposed as example test systems. NADP and methotrexate bind DHFR with positive cooperativity. Mammalian GDH exhibits negative cooperativity with respect to binding of NAD and NADPH coenzyme molecules, activation by ADP, and inhibition by GTP.

Functional epitopes for site 1 of human prolactin

Human prolactin (hPRL) binds two human prolactin receptor molecules, creating active heterotrimeric complexes. Receptors bind dissimilar hormone surfaces termed site 1 and site 2 in an obligate ordered process. We sought to map the functional epitopes in site 1 of hPRL. Extensive alanine mutagenesis (102 of the 199 residues) showed approximately 40% of these mutant hPRLs changed the ΔG for site 1 receptor binding. Six of these residues are within 3.5 Å of the receptor and form the site 1 functional epitopes. We identified a set of noncovalent interactions between these six residues and the receptor. We identified a second group of site 1 residues that are between 3.5 and 5 Å from the receptor where alanine mutations reduced the affinity. This second group has noncovalent interactions with other hormone residues and stabilized the topology of the functional epitopes by linking these to the body of the protein. Finally, we identified a third group of residues that are outside site 1 (>5 Å) and extend to site 2 and whose mutation to alanine significantly weakened receptor binding at site 1 of prolactin. These three groups of residues form a contiguous structural motif between sites 1 and 2 of human prolactin and may constitute structural features that functionally couple sites 1 and 2. This work identifies the residues that form the functional epitopes for site 1 of human prolactin and also identifies a set of residues that support the concept that sites 1 and 2 are functionally coupled by an allosteric mechanism.

Attributing functions to genes and gene products

A major focus of modern biochemical, biophysical and cell biological research is the attribution of function to elements of structure: gene products, genes and higher-order cellular structures. Misunderstandings and controversies can arise in connection with such assignments, in part because of the logical complexity inherent in the relating of structure to function and the failure to distinguish clearly among the different senses in which function can be imputed to elements of structure. I explore distinct ways in which functions are connected to structures and factors that contribute to the context-dependence of such associations so that the multiple senses of function can be made explicit.

Engineering protein allostery: 1.05 A resolution structure and enzymatic properties of a Na+-activated trypsin

Some trypsin-like proteases are endowed with Na(+)-dependent allosteric enhancement of catalytic activity, but this important mechanism has been difficult to engineer in other members of the family. Replacement of 19 amino acids in Streptomyces griseus trypsin targeting the active site and the Na(+)-binding site were found necessary to generate efficient Na(+) activation. Remarkably, this property was linked to the acquisition of a new substrate selectivity profile similar to that of factor Xa, a Na(+)-activated protease involved in blood coagulation. The X-ray crystal structure of the mutant trypsin solved to 1.05 A resolution defines the engineered Na(+) site and active site loops in unprecedented detail. The results demonstrate that trypsin can be engineered into an efficient allosteric protease, and that Na(+) activation is interwoven with substrate selectivity in the trypsin scaffold.

Effect of somatic mutation on DNA binding properties of anti-DNA autoantibodies

Autoantibodies that bind DNA are a hallmark of systemic lupus erythematosus. A subset of autoantibody*DNA complexes localize to kidney tissue and lead to damage and even death. 11F8, 9F11, and 15B10 are clonally related anti-DNA autoantibodies isolated from an autoimmune mouse. 11F8 binds ssDNA in a sequence-specific manner and causes tissue damage, while 9F11 and 15B10 bind ssDNA non-specifically and are benign. Among these antibodies, DNA binding properties are mediated by five amino acid differences in primary sequence. Thermodynamic and kinetic parameters associated with recognition of structurally different DNA sequences were determined for each antibody to provide insight toward recognition strategies, and to explore a link between binding properties and disease pathogenesis. A model of 11F8 bound to its high affinity consensus sequence provides a foundation for understanding the differences in thermodynamic and kinetic parameters between the three mAbs. Our data suggest that 11F8 utilizes the proposed ssDNA recognition motif including (Y32)V(L), a hydrogen bonding residue at (91)V(L), and an aromatic residue at the tip of the third heavy chain complementarity determining region. Interestingly, a somatic mutation to arginine at (31)V(H) in 11F8 may afford additional binding site contacts including (R31)V(H), (R96)V(H), and (R98)V(H) that could determine specificity.

A study of collective atomic fluctuations and cooperativity in the U1A-RNA complex based on molecular dynamics simulations

Cooperative interactions play an important role in recognition and binding in macromolecular systems. In this study, we find that cross-correlated atomic fluctuations can be used to identify cooperative networks in a protein-RNA system. The dynamics of the RRM-containing protein U1A-stem loop 2 RNA complex have been calculated theoretically from a 10 ns molecular dynamics (MD) simulation. The simulation was analyzed by calculating the covariance matrix of all atomic fluctuations. These matrix elements are then presented in the form of a two-dimensional grid, which displays fluctuations on a per residue basis. The results indicate the presence of strong, selective cross-correlated fluctuations throughout the RRM in U1A-RNA. The atomic fluctuations correspond well with previous biophysical studies in which a multiplicity of cooperative networks have been reported and indicate that the various networks identified in separate individual experiments are fluctuationally correlated into a hyper-network encompassing most of the RRM. The calculated results also correspond well with independent results from a statistical covariance analysis of 330 aligned RRM sequences. This method has significant implications as a predictive tool regarding cooperativity in the protein-nucleic acid recognition process.

Hirudin Binding Reveals Key Determinants of Thrombin Allostery

Thrombin exists in two allosteric forms, slow (S) and fast (F), that recognize natural substrates and inhibitors with significantly different affinities. Because under physiologic conditions the two forms are almost equally populated, investigation of thrombin function must address the contribution from the S and F forms and the molecular origin of their differential recognition of ligands. Using a panel of 79 Ala mutants, we have mapped for the first time the epitopes of thrombin recognizing a macromolecular ligand, hirudin, in the S and F forms. Hirudin binding is a relevant model for the interaction of thrombin with fibrinogen and PAR1 and is likewise influenced by the allosteric S-->F transition. The epitopes are nearly identical and encompass two hot spots, one in exosite I and the other in the Na+ site at the opposite end of the protein. The higher affinity of the F form is due to the preferential interaction of hirudin with Lys-36, Leu-65, Thr-74, and Arg-75 in exosite I; Gly-193 in the oxyanion hole; and Asp-221 and Asp-222 in the Na+ site. Remarkably, no correlation is found between the energetic and structural involvements of thrombin residues in hirudin recognition, which invites caution in the analysis of protein-protein interactions in general.

Molecular dissection of Na+ binding to thrombin

Na(+) binding near the primary specificity pocket of thrombin promotes the procoagulant, prothrombotic, and signaling functions of the enzyme. The effect is mediated allosterically by a communication between the Na(+) site and regions involved in substrate recognition. Using a panel of 78 Ala mutants of thrombin, we have mapped the allosteric core of residues that are energetically linked to Na(+) binding. These residues are Asp-189, Glu-217, Asp-222, and Tyr-225, all in close proximity to the bound Na(+). Among these residues, Asp-189 shares with Asp-221 the important function of transducing Na(+) binding into enhanced catalytic activity. None of the residues of exosite I, exosite II, or the 60-loop plays a significant role in Na(+) binding and allosteric transduction. X-ray crystal structures of the Na(+)-free (slow) and Na(+)-bound (fast) forms of thrombin, free or bound to the active site inhibitor H-d-Phe-Pro-Arg-chloromethyl-ketone, document the conformational changes induced by Na(+) binding. The slow --> fast transition results in formation of the Arg-187:Asp-222 ion pair, optimal orientation of Asp-189 and Ser-195 for substrate binding, and a significant shift of the side chain of Glu-192 linked to a rearrangement of the network of water molecules that connect the bound Na(+) to Ser-195 in the active site. The changes in the water network and the allosteric core explain the thermodynamic signatures linked to Na(+) binding and the mechanism of thrombin activation by Na(+). The role of the water network uncovered in this study establishes a new paradigm for the allosteric regulation of thrombin and other Na(+)-activated enzymes involved in blood coagulation and the immune response.

Dissecting substrate recognition by thrombin using the inactive mutant S195A

The catalytically inactive mutant S195A was used to study the interaction of thrombin with substrates under equilibrium conditions. By monitoring changes in intrinsic fluorescence, we measured dissociation constants for a variety of synthetic substrates, PAR peptides and the inhibitor PPACK. The S195A mutant retains the Na(+)-binding properties of the wild type, and substrate binding to the mutant is enhanced by the presence of Na(+). Temperature dependence studies allowed calculation of the thermodynamic parameters of substrate binding at the active site and showed a negligible deltaC(p). Titration of synthetic substrates carrying substitutions at the P1-P3 positions revealed energetics consistent with the specificity hierarchy identified in hydrolysis by the wild type. Titration with PAR peptides, which interact with both the active site and exosite I of thrombin, also showed consistency with the results obtained with the wild type at steady state. These findings demonstrate that inactive mutants of enzymes make it possible to dissect the equilibrium components linked to substrate binding and complement information on the kinetic properties of the wild type.

Characterization of the interactions between Asp141 and Phe236 in the Mn2+-l-malate binding of pigeon liver malic enzyme

The cytosolic malic enzyme from pigeon liver is very sensitive to the metal-catalysed oxidation systems. Our previous studies using the Cu2+-ascorbate as the oxidation system showed that the enzyme was oxidized and cleaved at several positions, including Asp141. The recently resolved crystal structure of pigeon liver malic enzyme revealed that Asp141 was near to the metal-binding site, but was not a direct metal ligand. However, Asp141 is located next to Phe236, which directly follows the metal ligands Glu234 and Asp235. Mutation at Asp141 caused a drastic effect on the metal-binding affinity of the enzyme. Since Asp141 and Phe236 are highly conserved in most species of malic enzyme, we used a double-mutant cycle to study the possible interactions between these two residues. Four single mutants [D141A (Asp141-->Ala), D141N, F236A and F236L] and four double mutants (D141A/F236A, D141N/F236A, D141A/F236L and D141N/F236L), plus the wild-type enzyme were successfully cloned, expressed and purified to homogeneity. The secondary, tertiary and quaternary structures of these mutants, as assessed by CD, fluorescence and analytical ultracentrifuge techniques, were similar to that of the wild-type enzyme. Initial velocity experiments were performed to derive the various kinetic parameters, which were used to analyse further the free energy change and the coupling energy (DeltaDeltaG(int)) between any two residues. The dissociation constants for Mn2+ ( K (d,Mn)) of the D141A and F236A mutants were increased by approx. 6- and 65-fold respectively, compared with that of the wild-type enzyme. However, the K (d,Mn) for the double mutant D141A/F236A was only increased by 150-fold. A coupling energy of -2.12 kcal/mol was obtained for Asp141 and Phe236. We suggest that Asp141 is involved in the second sphere of the metal-binding network of the enzyme.

The thrombin epitope recognizing thrombomodulin is a highly cooperative hot spot in exosite I

The functional epitope of thrombin recognizing thrombomodulin was mapped using Ala-scanning mutagenesis of 54 residues located around the active site, the Na(+) binding loop, the 186-loop, the autolysis loop, exosite I, and exosite II. The epitope for thrombomodulin binding is shaped as a hot spot in exosite I, centered around the buried ion quartet formed by Arg(67), Lys(70), Glu(77), and Glu(80), and capped by the hydrophobic residues Tyr(76) and Ile(82). The hot spot is a much smaller subset of the structural epitope for thrombomodulin binding recently documented by x-ray crystallography. Interestingly, the contribution of each residue of the epitope to the binding free energy shows no correlation with the change in its accessible surface area upon formation of the thrombin-thrombomodulin complex. Furthermore, residues of the epitope are strongly coupled in the recognition of thrombomodulin, as seen for the interaction of human growth hormone and insulin with their receptors. Finally, the Ala substitution of two negatively charged residues in exosite II, Asp(100) and Asp(178), is found unexpectedly to significantly increase thrombomodulin binding.

Characterization of the active sites in decay-accelerating factor

Decay-accelerating factor (DAF) is a complement regulator that dissociates autologous C3 convertases, which assemble on self cell surfaces. Its activity resides in the last three of its four complement control protein repeats (CCP2-4). Previous modeling on the nuclear magnetic resonance structure of CCP15-16 in the serum C3 convertase regulator factor H proposed a positively charged surface area on CCP2 extending into CCP3, and hydrophobic moieties between CCPs 2 and 3 as being primary convertase-interactive sites. To map the residues providing for the activity of DAF, we analyzed the functions of 31 primarily alanine substitution mutants based in part on this model. Replacing R69, R96, R100, and K127 in the positively charged CCP2-3 groove or hydrophobic F148 and L171 in CCP3 markedly impaired the function of DAF in both activation pathways. Significantly, mutations of K126 and F169 and of R206 and R212 in downstream CCP4 selectively reduced alternative pathway activity without affecting classical pathway activity. Rhesus macaque DAF has all the above human critical residues except for F169, which is an L, and its CCPs exhibited full activity against the human classical pathway C3 convertase. The recombinants whose function was preferentially impaired against the alternative pathway C3bBb compared with the classical pathway C4b2a were tested in classical pathway C5 convertase (C4b2a3b) assays. The effects on C4b2a and C4b2a3b were comparable, indicating that DAF functions similarly on the two enzymes. When CCP2-3 of DAF were oriented according to the crystal structure of CCP1-2 of membrane cofactor protein, the essential residues formed a contiguous region, suggesting a similar spatial relationship.

Dimensions of antigen recognition and levels of immunological specificity

Although recognition and specificity are among the most fundamental concepts in immunology, there is a common tendency to equate these notions with the fit, especially in terms of molecular shape, between interacting molecules. Even in the case of monovalent recognition, there are factors that contribute to the energetics of the interaction that are not readily accounted for by detailed structural analysis of the interacting (epitopic and paratopic) molecular surfaces. Consequently, recognition involves more than just the three spatial dimensions and time. Factors such as solute-solvent interactions, molecular crowding, and confinement, not directly related to the details of the intermolecular interface, can play crucial roles in determining both intrinsic affinity and differential intrinsic affinity. Furthermore, stating that a given structural subunit (e.g., amino acid) is recognized in a given noncovalent interaction does not clarify whether the structural subunit in question participates in the interaction through van der Waals contact, contribution to intrinsic affinity, or differential contribution to relative intrinsic affinities for two or more different ligands. Additional factors become relevant in considering the specificity exhibited in multivalent interactions, cell activation, and activation of the whole immune system. Therefore, specificity as defined for a monovalent binding event can diverge from specificity as it is defined for higher-order interactions. A corollary of this conclusion is that the composition of epitopes and paratopes, defined in terms of the structural elements for which substitutions have an effect on the specificity-defining measurement, can differ in different contexts despite complete conservation of the structures that physically make direct contact. An analysis of specificity at the organismal level suggests that the immune system does not recognize or respond to substances that correspond precisely to either nonself substances or to dangerous substances. An alternative notion for the molecular origins of immunological discrimination does not require that there be any single reason for immune responsiveness. This concept of what the immune system recognizes and responds to derives from the recognition that the ultimate function of the immune system is to contribute to survival and reproductive success through any available means.

Rational design of a potent anticoagulant thrombin

Thrombin acts as a procoagulant when it cleaves fibrinogen and promotes the formation of a fibrin clot and functions as an anticoagulant when it activates protein C with the assistance of the cofactor thrombomodulin. The dual function of thrombin in the blood poses the challenge to turn the enzyme into a potent anticoagulant by selectively abrogating fibrinogen cleavage. Using functional and structural data, we have rationally designed a thrombin mutant, W215A/E217A, that cleaves fibrinogen with a value of k(cat)/K(m) about 20,000-fold slower than wild-type but activates protein C in the presence of thrombomodulin with a specificity comparable with wild-type. This mutant demonstrates for the first time that the relative specificity of thrombin toward fibrinogen and protein C can be completely reversed.

Non-ionic detergent affects the conformation of a functionally active mutant of Bcl-X(L)

We found that a mutant, Bcl-X(L)(F131V), which was previously reported to have impaired binding capacity, can bind to Bax almost as strongly as wild-type Bcl-X(L). In the absence of detergent, the Bcl-X(L)(F131V) mutant adopts the same conformation as wild-type Bcl-X(L), as determined by circular dichroism spectroscopy and size-exclusion chromatography. However, non-ionic detergent induces a conformational change in the Bcl-X(L)(F131V) mutant and causes it to lose Bax-binding capacity. Wild-type Bcl-X(L), on the other hand, is more resistant to detergent-induced effects and retains its ability to bind Bax in the presence of detergent. Since it has been shown that the Bcl-X(L)(F131V) mutant has nearly the same anti-apoptotic activity as wild-type Bcl-X(L), it would be likely that the Bcl-X(L)(F131V) mutant can adopt the wild-type conformation, rather than the detergent-induced conformational state and can bind to Bax in vivo. Therefore, our data demonstrated that non-ionic detergent can have unpredicted effects on protein conformation, differential effects on wild-type and mutant Bcl-X(L) proteins in this case and may cause complications in the interpretation of in vitro binding studies.

Binding sites in Escherichia coli dihydrofolate reductase communicate by modulating the conformational ensemble

To explore how distal mutations affect binding sites and how binding sites in proteins communicate, an ensemble-based model of the native state was used to define the energetic connectivities between the different structural elements of Escherichia coli dihydrofolate reductase. Analysis of this model protein has allowed us to identify two important aspects of intramolecular communication. First, within a protein, pair-wise couplings exist that define the magnitude and extent to which mutational effects propagate from the point of origin. These pair-wise couplings can be identified from a quantity we define as the residue-specific connectivity. Second, in addition to the pair-wise energetic coupling between residues, there exists functional connectivity, which identifies energetic coupling between entire functional elements (i.e., binding sites) and the rest of the protein. Analysis of the energetic couplings provides access to the thermodynamic domain structure in dihydrofolate reductase as well as the susceptibility of the different regions of the protein to both small-scale (e.g., point mutations) and large-scale perturbations (e. g., binding ligand). The results point toward a view of allosterism and signal transduction wherein perturbations do not necessarily propagate through structure via a series of conformational distortions that extend from one active site to another. Instead, the observed behavior is a manifestation of the distribution of states in the ensemble and how the distribution is affected by the perturbation.

Distal recognition site for classical pathway convertase located in the C345C/netrin module of complement component C5

Previous studies focused on indels in the complement C345 protein family identified a number of potential protein-protein interaction sites in components C3 and C5. Here, one of these sites in C5, near the alpha-chain C terminus, was examined by alanine-scanning mutagenesis at 16 of the 18 non-alanine residues in the sequence KEALQIKYNFSF RYIYPLD. Alanine substitutions affected activities in the highly variable manner characteristic of binding sites. Substitutions at the lysine or either phenylalanine residue in the central KYNFSF sequence had the greatest effects, yielding mutants with <20% of the normal activity. These three mutants were also resistant to the classical pathway (CP) C5 convertase, with sensitivities roughly proportional to their hemolytic activities, but had normal susceptibilities to the cobra venom factor (CVF)-dependent convertase. Synthetic peptide MGKEALQIKYNFS-NH2 was found similarly to inhibit CP but not CVF convertase activation, and the effects of alanine substitutions in this peptide largely reflected those of the equivalent mutations in C5. These results indicate that residues KYNFSF form a novel, distal binding site for the CP, but not CVF convertase. This site lies approximately 880 residues downstream of the convertase cleavage site within a module that has been independently named C345C and NTR; this module is found in diverse proteins including netrins and tissue inhibitors of metalloproteinases.

The autolysis loop of activated protein C interacts with factor Va and differentiates between the Arg506 and Arg306 cleavage sites

The anticoagulant human plasma serine protease, activated protein C (APC), inactivates blood coagulation factors Va (FVa) and VIIIa. The so-called autolysis loop of APC (residues 301-316, equivalent to chymotrypsin [CHT] residues 142-153) has been hypothesized to bind FVa. In this study, site-directed mutagenesis was used to probe the role of the charged residues in this loop in interactions between APC and FVa. Residues Arg306 (147 CHT), Glu307, Lys308, Glu309, Lys311, Arg312, and Arg314 were each individually, or in selected combinations, mutated to Ala. The purified recombinant protein C mutants were characterized using activated partial thromboplastin time (APTT) clotting assays and FVa inactivation assays. Mutants 306A, 308A, 311A, 312A, and 314A had mildly reduced anticoagulant activity. Based on FVa inactivation assays and APTT assays using purified Gln506-FVa and plasma containing Gln506-FV, it appeared that these mutants were primarily impaired for cleavage of FVa at Arg506. Studies of the quadruple APC mutant (306A, 311A, 312A, and 314A) suggested that the autolysis loop provides for up to 15-fold discrimination of the Arg506 cleavage site relative to the Arg306 cleavage site. This study shows that the loop on APC of residues 306 to 314 defines an FVa binding site and accounts for much of the difference in cleavage rates at the 2 major cleavage sites in FVa.

The autolysis loop of activated protein C interacts with factor Va and differentiates between the Arg506 and Arg306 cleavage sites

The anticoagulant human plasma serine protease, activated protein C (APC), inactivates blood coagulation factors Va (FVa) and VIIIa. The so-called autolysis loop of APC (residues 301-316, equivalent to chymotrypsin [CHT] residues 142-153) has been hypothesized to bind FVa. In this study, site-directed mutagenesis was used to probe the role of the charged residues in this loop in interactions between APC and FVa. Residues Arg306 (147 CHT), Glu307, Lys308, Glu309, Lys311, Arg312, and Arg314 were each individually, or in selected combinations, mutated to Ala. The purified recombinant protein C mutants were characterized using activated partial thromboplastin time (APTT) clotting assays and FVa inactivation assays. Mutants 306A, 308A, 311A, 312A, and 314A had mildly reduced anticoagulant activity. Based on FVa inactivation assays and APTT assays using purified Gln506-FVa and plasma containing Gln506-FV, it appeared that these mutants were primarily impaired for cleavage of FVa at Arg506. Studies of the quadruple APC mutant (306A, 311A, 312A, and 314A) suggested that the autolysis loop provides for up to 15-fold discrimination of the Arg506 cleavage site relative to the Arg306 cleavage site. This study shows that the loop on APC of residues 306 to 314 defines an FVa binding site and accounts for much of the difference in cleavage rates at the 2 major cleavage sites in FVa.

The role of individual lysine residues in the basic patch on turnip cytochrome f for electrostatic interactions with plastocyanin in vitro

The role of electrostatic interactions in determining the rate of electron transfer between cytochrome f and plastocyanin has been examined in vitro with mutants of turnip cytochrome f and mutants of pea and spinach plastocyanins. Mutation of lysine residues Lys58, Lys65 and Lys187 of cytochrome f to neutral or acidic residues resulted in decreased binding constants and decreased rates of electron transfer to wild-type pea plastocyanin. Interaction of the cytochrome f mutant K187E with the pea plastocyanin mutant D51K gave a further decrease in electron transfer rate, indicating that a complementary charge pair at these positions could not compensate for the decreased overall charge on the proteins. Similar results were obtained with the interaction of the cytochrome f mutant K187E with single, double and triple mutants of residues in the acidic patches of spinach plastocyanin. These results suggest that the lysine residues of the basic patch on cytochrome f are predominantly involved in long-range electrostatic interactions with plastocyanin. However, analysis of the data using thermodynamic cycles provided evidence for the interaction of Lys187 of cytochrome f with Asp51, Asp42 and Glu43 of plastocyanin in the complex, in agreement with a structural model of a cytochrome f-plastocyanin complex determined by NMR.

Identification and modification of the uridine-binding site of the UDP-GalNAc (GlcNAc) pyrophosphorylase

UDP-GalNAc pyrophosphorylase (UDP-GalNAcPP; AGX1) catalyzes the synthesis of UDP-GalNAc from UTP and GalNAc-1-P. The 475-amino acid protein (57 kDa protein) also synthesizes UDP-GlcNAc at about 25% the rate of UDP-GalNAc. The cDNA for this enzyme, termed AGX1, was cloned in Escherichia coli, and expressed as an active enzyme that cross-reacted with antiserum against the original pig liver UDP-HexNAcPP. In the present study, we incubated recombinant AGX1 with N(3)-UDP-[(32)P]GlcNAc and N(3)-UDP-[(32)P]GalNAc probes to label the nucleotide-binding site. Proteolytic digestions of the labeled enzyme and analysis of the resulting peptides indicated that both photoprobes cross-linked to one 24-amino acid peptide located between residues Val(216) and Glu(240). Four amino acids in this peptide were found to be highly conserved among closely related enzymes, and each of these was individually modified to alanine. Mutation of Gly(222) to Ala in the peptide almost completely eliminated UDP-GlcNAc and UDP-GalNAc synthesis, while mutation of Gly(224) to Ala, almost completely eliminated UDP-GalNAc synthesis, but UDP-GlcNAc was only diminished by 50%. Both of these mutations also resulted in almost complete loss of the ability of the mutated proteins to cross-link N(3)-UDP-[(32)P]GlcNAc or N(3)-UDP-[(32)P]GalNAc. On the other hand, mutations of either Pro(220) or Tyr(227) to Ala did not greatly affect enzymatic activity, although there was some reduction in the ability of these proteins to cross-link the photoaffinity probes. We also mutated Gly(111) to Ala since this amino acid was reported to be necessary for catalysis (Mio, T., Yabe, T., Arisawa, M., and Yamada-Okabe, H. (1998) J. Biol. Chem. 273, 14392-14397). The Gly(111) to Ala mutant lost all enzymatic activity, but interestingly enough, this mutant protein still cross-linked the radioactive N(3)-UDP-GlcNAc although not nearly as well as the wild type. On the other hand, mutation of Arg(115) to Ala had no affect on enzymatic activity although it also reduced the amount of cross-linking of N(3)-UDP-[(32)P]GlcNAc. These studies help to define essential amino acids at or near the nucleotide-binding site and the catalytic site, as well as peptides involved in binding and catalysis.

Investigation of phosphotyrosine recognition by the SH2 domain of the Src kinase

The binding of tyrosine phosphorylated targets by SH2 domains is required for propagation of many cellular signals in higher eukaryotes; however, the determinants of phosphotyrosine (pTyr) recognition by SH2 domains are not well understood. In order to identify the attributes of pTyr required for high affinity interaction with SH2 domains, the binding of the SH2 domain of the Src kinase (Src SH2 domain) to a dephosphorylated peptide, a phosphoserine-containing peptide, and the amino acid pTyr was studied using titration calorimetry and compared with the binding of a high affinity tyrosyl phosphopeptide. The dephosphorylated peptide and the phosphoserine containing peptide both bind extremely weakly to the Src SH2 domain (DeltaGo (dephosphorylated)=-3.6 kcal/mol, DeltaGo (phosphoserine) >-3.7 kcal/mol); however, the DeltaGo value of pTyr binding is more favorable (-4.7 kcal/mol, or 50 % of the entire binding free energy of a high affinity tyrosyl phosphopeptide). These results indicate that both the phosphate and the tyrosine ring of the pTyr are critical determinants of high affinity binding. Alanine mutagenesis was also used to evaluate the energetic contribution to binding of ten residues located in the pTyr-binding site. Mutation of the strictly conserved Arg betaB5 resulted in a large increase in DeltaGo (DeltaDeltaGo=3.2 kcal/mol) while elimination of the other examined residues each resulted in a significantly smaller (DeltaDeltaGo<1.4 kcal/mol) reduction in affinity, indicating that Arg betaB5 is the single most important determinant of pTyr recognition. However, mutation of Cys betaC3, a residue unique to the Src SH2 domain, surprisingly increased affinity by eightfold (DeltaDeltaGo=-1.1 kcal/mol). Using a double mutant cycle analysis, it was revealed that residues of the pTyr-binding pocket are not coupled to the peptide residues C-terminal to the pTyr. In addition, comparison of each residue's DeltaDeltaGo value upon mutation with that residue's sequence conservation among SH2 domains revealed only a modest correlation between a residue's energetic contribution to pTyr recognition and its conservation throughout evolution. The results of this investigation highlight the importance of a single critical interaction, the buried ionic bond between the phosphate of the pTyr and Arg betaB5 of the SH2 domain, driving the binding of SH2 domains to tyrosine phosphorylated targets.

Incorporation of noncoded amino acids into the N-terminal domain 1-47 of hirudin yields a highly potent and selective thrombin inhibitor

Hirudin is an anticoagulant polypeptide isolated from a medicinal leech that inhibits thrombin with extraordinary potency (Kd = 0.2-1.0 pM) and selectivity. Hirudin is composed of a compact N-terminal region (residues 1-47, cross-linked by three disulfide bridges) that binds to the active site of thrombin, and a flexible C-terminal tail (residues 48-64) that interacts with the exosite I of the enzyme. To minimize the sequence of hirudin able to bind thrombin and also to improve its therapeutic profile, several N-terminal fragments have been prepared as potential anticoagulants. However, the practical use of these fragments has been impaired by their relatively poor affinity for the enzyme, as given by the increased value of the dissociation constant (Kd) of the corresponding thrombin complexes (Kd = 30-400 nM). The aim of the present study is to obtain a derivative of the N-terminal domain 1-47 of hirudin displaying enhanced inhibitory potency for thrombin compared to the natural product. In this view, we have synthesized an analogue of fragment 1-47 of hirudin HM2 in which Val1 has been replaced by tert-butylglycine, Ser2 by Arg, and Tyr3 by beta-naphthylalanine, to give the BugArgNal analogue. The results of chemical and conformational characterization indicate that the synthetic peptide is able to fold efficiently with the correct disulfide topology (Cys6-Cys14, Cys16-Cys28, Cys22-Cys37), while retaining the conformational properties of the natural fragment. Thrombin inhibition data indicate that the effects of amino acid replacements are perfectly additive if compared to the singly substituted analogues (De Filippis V, Quarzago D, Vindigni A, Di Cera E, Fontana A, 1998, Biochemistry 37:13507-13515), yielding a molecule that inhibits the fast or slow form of thrombin by 2,670- and 6,818-fold more effectively than the natural fragment, and that binds exclusively at the active site of the enzyme with an affinity (Kd,fast = 15.4 pM, Kd,slow = 220 pM) comparable to that of full-length hirudin (Kd,fast = 0.2 pM, Kd,slow = 5.5 pM). Moreover, BugArgNal displays absolute selectivity for thrombin over the other physiologically important serine proteases trypsin, plasmin, factor Xa, and tissue plasminogen activator, up to the highest concentration of inhibitor tested (10 microM).

Interaction between the peripheral site residues of human butyrylcholinesterase, D70 and Y332, in binding and hydrolysis of substrates

Human butyrylcholinesterase displays substrate activation with positively charged butyrylthiocholine (BTC) as the substrate. Peripheral anionic site (PAS) residues D70 and Y332 appear to be involved in the initial binding of charged substrates and in activation control. To determine the contribution of PAS residues to binding and hydrolysis of quaternary substrates and activation control, the single mutants D70G/Y and Y332F/A/D and the double mutants Y332A/D70G and Y332D/D70Y were studied. Steady-state hydrolysis of the charged substrates, BTC and succinyldithiocholine, and the neutral ester o-nitrophenyl butyrate was measured. In addition, inhibition of wild-type and mutant enzymes by tetramethylammonium was investigated, at low concentrations of BTC. Single and double mutants of D70 and Y332 showed little or no substrate activation, suggesting that both residues were important for activation control. The effects of double mutations on D70 and Y332 were complex. Double-mutant cycle analysis provided evidence for interaction between these residues. The category of interaction (either synergistic, additive, partially additive or antagonistic) was found to depend on the nature of the substrate and on measured binding or kinetic parameters. This complexity reflects both the cross-talk between residues involved in the sequential formation of productive Michaelian complexes and the effect of peripheral site residues on catalysis. It is concluded that double mutations on the PAS induce a conformational change in the active site gorge of butyrylcholinesterase that can alter both substrate binding and enzyme acylation.

RNA recognition by the human U1A protein is mediated by a network of local cooperative interactions that create the optimal binding surface

One of the most common structural motifs in RNA-binding proteins is the RNA-binding domain (RBD). These domains share a common alpha/beta sandwich tertiary fold, and are highly conserved, though they bind diverse RNA targets with a wide range of binding affinities. The N-terminal RNA-binding domain (RBD1) of the human U1A protein binds specifically to stem/loop II of the U1 snRNA with sub-nanomolar affinity. Solvent-exposed aromatic residues on the beta-sheet surface are highly conserved among RBD domains; in RBD1, these are Tyr13 and Phe56, with a unique Gln at position 54. Effects of substitutions at these positions were examined using energetic pairwise coupling to describe the communication between these residues in both the free and RNA-bound states of the protein. 15N NMR experiments were used to determine effects of the beta-sheet substitutions on the structural and dynamic properties of this domain. The combination of thermodynamic pairwise coupling and 15N-backbone dynamics provides direct evidence for local cooperative interactions among Y13, Q54, and F56, and a non-conserved loop that directly affect RNA-binding. The results describe how conserved and non-conserved regions of an RBD can communicate with each other to mediate recognition of the RNA.

Role of P225 and the C136-C201 disulfide bond in tissue plasminogen activator

The protease domain of tissue plasminogen activator (tPA), a key fibrinolytic enzyme, was expressed in Escherichia coli with a yield of 1 mg per liter of media. The recombinant protein was titrated with the Erythrina caraffa trypsin inhibitor (ETI) and characterized in its interaction with plasminogen and the natural inhibitor plasminogen activator inhibitor-1 (PAI-1). Analysis of the catalytic properties of tPA using a library of chromogenic substrates carrying substitutions at P1, P2, and P3 reveals a strong preference for Arg over Lys at P1, unmatched by other serine proteases like thrombin or trypsin. In contrast to these proteases and plasmin, tPA shows little or no preference for Pro over Gly at P2. A specific inhibition of tPA by Cu2+ was discovered. The divalent cation presumably binds to H188 near D189 in the primary specificity pocket and inhibits substrate binding in a competitive manner with a Kd = 19 microM. In an attempt to engineer Na+ binding and enhanced catalytic activity in tPA, P225 was replaced with Tyr, the residue present in Na+-dependent allosteric serine proteases. The P225Y mutation did not result in cation binding, but caused a significant loss of specificity (up to 100-fold) toward chromogenic substrates and plasminogen and considerably reduced the inhibition by PAI-1 and ETI. Interestingly, the P225Y substitution enhanced the ability of Cu2+ to inhibit the enzyme. Elimination of the C136-C201 disulfide bond, that is absent in all Na+-dependent allosteric serine proteases, significantly enhanced the yield (5 mg per liter of media) of expression in E. coli, but caused no changes in the properties of the enzyme whether residue 225 was Pro or Tyr. These findings point out an unanticipated crucial role for residue 225 in controlling the catalytic activity of tPA, and suggest that engineering of a Na+-dependent allosteric enhancement of catalytic activity in this enzyme, must involve substantial changes in the region homologous to the Na+ binding site of allosteric serine proteases.

RNA binding mediates the local cooperativity between the beta-sheet and the C-terminal tail of the human U1A RBD1 protein

Pairwise coupling theory is applied here to determine the energetic interactions between two elements of the N-terminal RNA binding domain (RBD) of the human U1A protein. The novel application of the theory to this system incorporates both measurements of protein stability and RNA binding to define thermodynamic cycles. In this first example of the application, two regions of the protein are selected for study: tyrosine 13, one of the conserved aromatic residues on the surface of the beta-sheet, and the C-terminal tail of the RBD. The six initial pairwise coupling free energies derived from this system describe the communication between these positions, both in the free and RNA-bound states of the protein. The results show that in the absence of RNA, these two elements of the protein act independently. However, when RNA is bound, there is indirect coupling between Tyr13 and the tail, mediated through the RNA. Subsequent thermodynamic cycles involving additional perturbations to the C-terminal tail further define the communication between the C terminus and the beta-sheet. This work demonstrates the general applicability of the pairwise coupling theory to protein:nucleic acid interactions, and illustrates the necessity of such analyses to describe the network of energetic interactions that comprise RNA recognition by this RBD.