Double-mutant Cycle Research Articles

Docking of μ-Conotoxin GIIIA in the Sodium Channel Outer Vestibule

μ-Conotoxin GIIIA (μ-CTX) is a high-affinity ligand for the outer vestibule of selected isoforms of the voltage-gated Na+ channel. The detailed bases for the toxin’s high affinity binding and isoform selectivity are unclear. The outer vestibule is lined by four pore-forming (P) loops, each with an acidic residue near the mouth of the vestibule. μ-CTX has seven positively charged residues that may interact with these acidic P-loop residues. Using pair-wise alanine replacement of charged toxin and channel residues, in conjunction with double mutant cycle analysis, we determined coupling energies for specific interactions between each P-loop acidic residue and selected toxin residues to systematically establish quantitative restraints on the toxin orientation in the outer vestibule. Xenopus oocytes were injected with the mutant or native Na+ channel mRNA, and currents measured by two-electrode voltage clamp. Mutant cycle analysis revealed novel, strong, toxin-channel interactions between K9/E403, K11/D1241, K11/D1532, and R19/D1532. Experimentally determined coupling energies for interacting residue pairs provided restraints for molecular dynamics simulations of μ-CTX docking. Our simulations suggest a refined orientation of the toxin in the pore, with toxin basic side-chains playing key roles in high-affinity binding. This modeling also provides a set of testable predictions toxin-channel interactions, hitherto undescribed, that may contribute to high-affinity binding and channel isoform selectivity.

Influence of hPin1 WW N‐terminal domain boundaries on function, protein stability, and folding

An N-terminally truncated and cooperatively folded version (residues 6-39) of the human Pin1 WW domain (hPin1 WW hereafter) has served as an excellent model system for understanding triple-stranded beta-sheet folding energetics. Here we report that the negatively charged N-terminal sequence (Met1-Ala-Asp-Glu-Glu5) previously deleted, and which is not conserved in highly homologous WW domain family members from yeast or certain fungi, significantly increases the stability of hPin1 WW (approximately 4 kJ mol(-1) at 65 degrees C), in the context of the 1-39 sequence based on equilibrium measurements. N-terminal truncations and mutations in conjunction with a double mutant cycle analysis and a recently published high-resolution X-ray structure of the hPin1 cis/trans-isomerase suggest that the increase in stability is due to an energetically favorable ionic interaction between the negatively charged side chains in the N terminus of full-length hPin1 WW and the positively charged epsilon-ammonium group of residue Lys13 in beta-strand 1. Our data therefore suggest that the ionic interaction between Lys13 and the charged N terminus is the optimal solution for enhanced stability without compromising function, as ascertained by ligand binding studies. Kinetic laser temperature-jump relaxation studies reveal that this stabilizing interaction has not formed to a significant extent in the folding transition state at near physiological temperature, suggesting a differential contribution of the negatively charged N-terminal sequence to protein stability and folding rate. As neither the N-terminal sequence nor Lys13 are highly conserved among WW domains, our data further suggest that caution must be exercised when selecting domain boundaries for WW domains for structural, functional, or thermodynamic studies.

Role of Aromatic Side Chains in the Folding and Thermodynamic Stability of Integral Membrane Proteins

Aromatic residues are frequently found in helical and beta-barrel integral membrane proteins enriched at the membrane-water interface. Although the importance of these residues in membrane protein folding has been rationalized by thermodynamic partition measurements using peptide model systems, their contribution to the stability of bona fide membrane proteins has never been demonstrated. Here, we have investigated the contribution of interfacial aromatic residues to the thermodynamic stability of the beta-barrel outer membrane protein OmpA from Escherichia coli in lipid bilayers by performing extensive mutagenesis and equilibrium folding experiments. Isolated interfacial tryptophanes contribute -2.0 kcal/mol, isolated interfacial tyrosines contribute -2.6 kcal/mol, and isolated interfacial phenylalanines contribute -1.0 kcal/mol to the stability of this protein. These values agree well with the prediction from the Wimley-White interfacial hydrophobicity scale, except for tyrosine residues, which contribute more than has been expected from the peptide models. Double mutant cycle analysis reveals that interactions between aromatic side chains become significant when their centroids are separated by less than 6 A but are nearly insignificant above 7 A. Aromatic-aromatic side chain interactions are on the order of -1.0 to -1.4 kcal/mol and do not appear to depend on the type of aromatic residue. These results suggest that the clustering of aromatic side chains at membrane interfaces provides an additional heretofore not yet recognized driving force for the folding and stability of integral membrane proteins.

On the Dynamic Nature of the Transition State for Protein–Protein Association as Determined by Double-mutant Cycle Analysis and Simulation

The process of protein–protein association starts with their random collision, which may develop into an encounter complex followed by a transition state and final complex formation. Here we aim to experimentally characterize the nature of the transition state of protein–protein association for three different protein–protein interactions; Barnase-Barstar, TEM1-BLIP and IFNα2-IFNAR2, and use the data to model the transition state structures. To model the transition state, we determined inter-protein distance-constraints of the activated complex by using double mutant cycles (DMC) assuming that interacting residues are spatially close. Significant ΔΔ G ‡ int values were obtained only between residues on Barnase and Barstar. However, introducing specific mutations that optimize the charge complementarity between BLIP and TEM1 resulted in the introduction of significant ΔΔ G ‡ int values also between residues of these two proteins. While electrostatic interactions make major contributions towards stabilizing the transition state, we show two examples where steric hindrance exerts an effect on the transition state as well. To model the transition-state structures from the experimental ΔΔ G ‡ int values, we introduced a method for structure perturbation, searching for those inter-protein orientations that best support the experimental ΔΔ G ‡ int values. Two types of transition states were found, specific as observed for Barnase-Barstar and the electrostatically optimized TEM1-BLIP mutants, and diffusive as shown for wild-type TEM1-BLIP and IFNα2-IFNAR2. The specific transition states are characterized by defined inter-protein orientations, which cannot be modeled for the diffusive transition states. Mutations introduced through rational design can change the transition state from diffusive to specific. Together, these data provide a structural view of the mechanism allowing rates of association to differ by five orders of magnitude between different protein complexes.

A Thermodynamic Ligand Binding Study of the Third PDZ Domain (PDZ3) from the Mammalian Neuronal Protein PSD-95

The thermodynamic parameters associated with the binding of several series of linear peptides to the third PDZ domain (PDZ3) of the postsynaptic density 95 protein (PSD-95) have been measured using isothermal titration calorimetry (ITC). Two strategies were pursued in developing these binding ligands: (1) systematic N-terminal truncation of sequences derived from the C-terminal regions of identified PDZ3-binding proteins (CRIPT, neuroligin-1, and citron) and (2) selective mutation of specific positions within a consensus hexapeptide (KKETEV) known to bind PDZ3. Each synthetically prepared peptide was used to titrate PDZ3, which yielded the changes in Gibbs free energy (DeltaG), enthalpy (DeltaH), and entropy (TDeltaS) for the binding event. Selected peptides were subjected to additional analysis, which entailed (1) measuring the change in heat capacity (DeltaCp) upon association, to assess the character of the binding interface, and (2) constructing thermodynamic double mutant cycles, to determine the presence of cooperative effects. From the first series, the CRIPT protein proved to be the better source for higher affinity sequences. From the second series, enhanced binding was associated with peptides that closely adhered to the established motif for class I PDZ domain C-termini, X-(T/S)-X-(V/I/L), and more specifically to a narrower motif of X-T-X-V. Further, in both series a length of six residues was necessary and sufficient to capture maximal affinity. In addition, there were significant influences upon binding by modifying the abutting "X" positions. The cumulative results provide greater detail into the specific nature of ligand binding to PDZ3 and will assist in the development of selective molecular probes for the study of this and structurally homologous PDZ domains.

Chemical Double‐Mutant Cycles: Dissecting Non‐Covalent Interactions

AbstractChemInform is a weekly Abstracting Service, delivering concise information at a glance that was extracted from about 200 leading journals. To access a ChemInform Abstract, please click on HTML or PDF.

Reconciling the Solution and X-ray Structures of the Villin Headpiece Helical Subdomain: Molecular Dynamics Simulations and Double Mutant Cycles Reveal a Stabilizing Cation−π Interaction

The 36-residue helical subdomain of the villin headpiece, HP36, is one of the smallest cooperatively folded proteins, folding on the microsecond time scale. The domain is an extraordinarily popular model system for both experimental and computational studies of protein folding. The structure of HP36 has been determined using X-ray crystallography and NMR spectroscopy, with the resulting structures exhibiting differences in helix packing, van der Waals contacts, and hydrogen bonding. It is important to determine the solution structure of HP36 with as much accuracy as possible since this structure is widely used as a reference for simulations and experiments. We complement the existing data by using all-atom molecular dynamics simulations with explicit solvent to evaluate which of the experimental models is the better representation of HP36 in solution. After simulation for 50 ns initiated with the NMR structure, we observed that the protein spontaneously adopts structures with a backbone conformation, core packing, and C-capping motif on the third helix that are more consistent with the crystal structure. We also examined hydrogen bonding and side chain packing interactions between D44 and R55 and between F47 and R55, respectively, which were observed in the crystal structure but not in the NMR-based solution structure. Simulations showed large fluctuations in the distance between D44 and R55, while the distance between F47 and R55 remained stable, suggesting the formation of a cation-pi interaction between those residues. Experimental double mutant cycles confirmed that the F47-R55 pair has a larger energetic coupling than the D44-R55 interaction. Overall, these combined experimental and computational studies show that the X-ray crystal structure is the better reference structure for HP36 in solution at neutral pH. Our analysis also shows how detailed molecular dynamics simulations combined with experimental validation can help bridge the gap between NMR and crystallographic methods.

Double-Mutant Cycle Scanning of the Interaction of a Peptide Ligand and Its G Protein-Coupled Receptor

The interaction between the yeast G protein-coupled receptor (GPCR), Ste2p, and its alpha-factor tridecapeptide ligand was subjected to double-mutant cycle scanning analysis by which the pairwise interaction energy of each ligand residue with two receptor residues, N205 and Y266, was determined. The mutations N205A and Y266A were previously shown to result in deficient signaling but cause only a 2.5-fold and 6-fold decrease, respectively, in the affinity for alpha-factor. The analysis shows that residues at the amine terminus of alpha-factor interact strongly with N205 and Y266 whereas residues in the center and at the carboxyl terminus of the peptide interact only weakly if at all with these receptor residues. Multiple-mutant thermodynamic cycle analysis was used to assess whether the energies of selected pairwise interactions between residues of the alpha-factor peptide changed upon binding to Ste2p. Strong positive cooperativity between residues 1 through 4 of alpha-factor was observed during receptor binding. In contrast, no thermodynamic evidence was found for an interaction between a residue near the carboxyl terminus of alpha-factor (position 11) and one at the N-terminus (position 3). The study shows that multiple-mutant cycle analyses of the binding of an alanine-scanned peptide to wild-type and mutant GPCRs can provide detailed information on contributions of inter- and intramolecular interactions to the binding energy and potentially prove useful in developing 3D models of ligand docked to its receptor.

Thermodynamic stability of a cold‐adapted protein, type III antifreeze protein, and energetic contribution of salt bridges

A thermodynamic analysis of a cold-adapted protein, type III anti-freeze protein (AFP), was carried out. The results indicate that the folding equilibrium of type III AFP is a reversible, unimolecular, two-state process with no populated intermediates. Compared to most mesophilic proteins whose folding is two-state, the psychrophilic type III AFP has a much lower thermodynamic stability at 25 degrees C, approximately 3 kcal/mol, and presents a remarkably downshifted stability-temperature curve, reaching a maximum of 5 kcal/mol around 0 degrees C. Type III AFPs contain few and non-optimally distributed surface charges relative to their mesophilic homologs, the C-terminal domains of sialic acid synthases. We used thermodynamic double mutant cycles to evaluate the energetic role of every surface salt bridge in type III AFP. Two isolated salt bridges provided no contribution to stability, while the Asp36-Arg39 salt bridge, involved in a salt bridge network with the C-terminal carboxylate, had a substantial contribution (approximately 1 kcal/mol). However, this contribution was more than counteracted by the destabilizing effect of the Asp36 carboxylate itself, whose removal led to a net 30% increase in stability at 25 degrees C. This study suggests that type III AFPs may have evolved for a minimally acceptable stability at the restricted, low temperature range (around 0 degrees C) at which AFPs must function. In addition, it indicates that salt bridge networks are used in nature also for the stability of psychrophilic proteins, and has led to a type III AFP variant of increased stability that could be used for biotechnological purposes.

Folding Cooperativity in RNA and DNA Is Dependent on Position in the Helix

Secondary structural motifs play essential roles in the folding and function of RNA and DNA molecules. Previous work from our lab compared the folding of small DNA and RNA hairpin loops containing a sheared GA pair [Moody, E. M., Feerar, J. C., and Bevilacqua, P. C. (2004) Biochemistry 43, 7992-7998]. We found that the small DNA hairpins fold in a highly cooperative manner with indirect coupling, while their RNA counterparts fold in a much less cooperative fashion and display direct coupling. Herein, we extend this study to the double-stranded helix. We carried out double mutant cycles on base pairs having identical nearest-neighbor contexts but located in either external or internal helical registers. In the external register, both RNA and DNA exhibit extensive folding cooperativity between the penultimate and terminal base pair, which is independent of mismatch identity. In contrast, DNA exhibits virtually no folding cooperativity in the center of the helix, while RNA maintains substantial coupling, which is dependent on mismatch identity. Two models account for these non-nearest-neighbor effects: one involves the unfavorable entropy of helix initiation common to DNA and RNA, and the other involves steric and electrostatic strain peculiar to RNA. These data show that RNA can display cooperativity less than, greater than, or equal to that of DNA depending on context and position.

Does the A·T or G·C Base-Pair Possess Enhanced Stability? Quantifying the Effects of CH···O Interactions and Secondary Interactions on Base-Pair Stability Using a Phenomenological Analysis and ab Initio Calculations

An empirically based relationship between overall complex stability (-DeltaG degrees ) and various possible component interactions is developed to probe the question of whether the A.T/U and G.C base-pairs exhibit enhanced stability relative to similarly hydrogen-bonded complexes. This phenomenological approach suggests ca. 2-2.5 kcal mol-1 in additional stability for A.T owing to a group interaction containing a CH...O contact. Pairing geometry and the role of the CH...O interaction in the A.T base-pair were also probed using MP2/6-31+G(d,p) calculations and a double mutant cycle. The ab initio studies indicated that Hoogsteen geometry is preferred over Watson-Crick geometry in A.T by ca. 1 kcal mol-1. Factors that might contribute to the preference for Hoogsteen geometry are a shorter CH...O contact, a favorable alignment of dipoles, and greater distances between secondary repulsive sites. The CH...O interaction was also investigated in model complexes of adenine with ketene and isocyanic acid. The ab initio calculations support the result of the phenomenological approach that the A.T base-pair does have enhanced stability relative to hydrogen-bonded complexes with just N-H...N and N-H...O hydrogen bonds.

Substituent effects on aromatic stacking interactions

Synthetic supramolecular zipper complexes have been used to quantify substituent effects on the free energies of aromatic stacking interactions. The conformational properties of the complexes have been characterised using NMR spectroscopy in CDCl(3), and by comparison with the solid state structures of model compounds. The structural similarity of the complexes makes it possible to apply the double mutant cycle method to evaluate the magnitudes of 24 different aromatic stacking interactions. The major trends in the interaction energy can be rationalised using a simple model based on electrostatic interactions between the pi-faces of the two aromatic rings. However, electrostatic interactions between the substituents of one ring and the pi-face of the other make an additional contribution, due to the slight offset in the stacking geometry. This property makes aromatic stacking interactions particularly sensitive to changes in orientation as well as the nature and location of substituents.

Role of Pairwise Interactions between M1 and M2 Domains of the Nicotinic Receptor in Channel Gating

The adult form of the nicotinic acetylcholine receptor (AChR) consists of five subunits ( α 2 βɛδ), each having four transmembrane domains (M1–M4). The atomic model of the nicotinic acetylcholine receptor shows that the pore-lining M2 domains make no extensive contacts with the rest of the transmembrane domains. However, there are several sites where close appositions between segments occur. It has been suggested that the pair αM1-F15′ and αM2-L11′ is one of the potential interactions between segments. To determine experimentally if these residues are interacting and to explore if this interhelical interaction is essential for channel gating, we combined mutagenesis with single-channel kinetic analysis. Mutations in αM1-F15′ lead to profound changes in the opening rate and slighter changes in the closing rate. Channel gating is impaired as the volume of the residue increases. Rate-equilibrium linear free-energy relationship analysis reveals an ∼70% open-state-like environment for αM1-F15′ at the transition state of the gating reaction, suggesting that it moves early during the gating process. Replacing the residue at αM1-15′ by that at αM2-11′ and vice versa profoundly alters gating, but the combination of the two mutations restores gating to near normal, indicating that αM1-F15′ and αM2-L11′ are interchangeable. Double-mutant cycle analysis shows that these residues are energetically coupled. Thus, the interaction between M1 and M2 plays a key role in channel gating.

Chemical double-mutant cycles: dissecting non-covalent interactions

Thermodynamic double-mutant cycles and triple-mutant boxes are widely employed for the experimental quantification of non-covalent interactions and cooperative effects in proteins. This review describes the application of these powerful methodologies to the study of non-covalent interactions in synthetic systems.

Electrostatic couplings in OmpA ion-channel gating suggest a mechanism for pore opening

The molecular forces that drive structural transitions between the open and closed states of channels and transporters are not well understood. The gate of the OmpA channel is formed by the central Glu52-Arg138 salt bridge, which can open to form alternate ion pairs with Lys82 and Glu128. To gain deeper insight into the channel-opening mechanism, we measured interaction energies between the relevant side chains by double-mutant cycle analysis and correlated these with the channel activities of corresponding point mutants. The closed central salt bridge has a strong interaction energy of -5.6 kcal mol(-1), which can be broken by forming the open-state salt bridge Glu52-Lys82 (DeltaDeltaG(Inter) = -3.5 kcal mol(-1)) and a weak interaction between Arg138 and Glu128 (DeltaDeltaG(Inter) = -0.6 kcal mol(-1)). A covalent disulfide bond in place of the central salt bridge completely blocks the channel. Growth assays indicate that this gating mechanism could physiologically contribute to the osmoprotection of Escherichia coli cells from environmental stress.

Stability of 100 Homo and Heterotypic Coiled−Coil a−a‘ Pairs for Ten Amino Acids (A, L, I, V, N, K, S, T, E, and R)

We present the thermal stability monitored by circular dichroism (CD) spectroscopy at 222 nm of 100 heterodimers that contain all possible coiled-coil a-a' pairs for 10 amino acids (I, V, L, N, A, K S, T, E, and R). This includes the stability of 36 heterodimers for 6 amino acids (I, V, L, N, A, and K) previously described and 64 new heterodimers including the 4 amino acids (S, T, E, and R). We have calculated a double mutant alanine thermodynamic cycle to determine a-a' pair coupling energies to evaluate which a-a' pairs encourage specific dimerization partners. The four new homotypic a-a' pairs (T-T, S-S, R-R, E-E) are repulsive relative to A-A and have destabilizing coupling energies. Among the 90 heterotypic a-a' pairs, the stabilizing coupling energies contain lysine or arginine paired with either an aliphatic or a polar amino acid. The range in coupling energies for each amino acid reveals its potential to regulate dimerization specificity. The a-a' pairs containing isoleucine and asparagine have the greatest range in coupling energies and thus contribute dramatically to dimerization specificity, which is to encourage homodimerization. In contrast, the a-a' pairs containing charged amino acids (K, R, and E) show the least range in coupling energies and promiscuously encourage heterodimerization.

The Structural Context of Disease-causing Mutations in Gap Junctions

Gap junctions form intercellular channels that mediate metabolic and electrical signaling between neighboring cells in a tissue. Lack of an atomic resolution structure of the gap junction has made it difficult to identify interactions that stabilize its transmembrane domain. Using a recently computed model of this domain, which specifies the locations of each amino acid, we postulated the existence of several interactions and tested them experimentally. We introduced mutations within the transmembrane domain of the gap junction-forming protein connexin that were previously implicated in genetic diseases and that apparently destabilized the gap junction, as evidenced here by the absence of the protein from the sites of cell-cell apposition. The model structure helped identify positions on adjacent helices where second-site mutations restored membrane localization, revealing possible interactions between residue pairs. We thus identified two putative salt bridges and one pair involved in packing interactions in which one disease-causing mutation suppressed the effects of another. These results seem to reveal some of the physical forces that underlie the structural stability of the gap junction transmembrane domain and suggest that abrogation of such interactions bring about some of the effects of disease-causing mutations.

Pairwise Coupling Analysis of Helical Junction Hydrogen Bonding Interactions in Luteoviral RNA Pseudoknots

A 28-nucleotide mRNA pseudoknot that overlaps the P1 and P2 genes of sugarcane yellow leaf virus (ScYLV) stimulates -1 ribosomal frameshifting. The in vitro frameshifting efficiency is decreased >or=8-fold upon substitution of the 3'-most loop 2 nucleotide (C27) with adenosine, which accepts a hydrogen bond from the 2'-OH group of C14 in stem S1. The solution structures of the wild-type (WT) and C27A ScYLV RNA pseudoknots show that while the RNAs adopt virtually identical overall structures, there are significant structural differences at the helical junctions of the two RNAs. Specifically, C8(+) in loop L1 in the C8(+).(G12.C28) L1-S2 major groove base triple is displaced by approximately 2.3 A relative to the accepting stem 2 base pair (G12.C28) in the C27A RNA. Here, we use a double mutant cycle approach to analyze the pairwise coupling of the C8(+).(G12.C28)...C27.(C14-G7) and ...A27.(C14-G7) hydrogen bonds in the WT and C27A ScYLV RNAs, respectively, and compare these findings with previous results from the beet western yellows virus (BWYV) RNA. We find that the pairwise coupling free energy (delta(AB)(i)) is favorable for the WT RNA (-0.7 +/- 0.1 kcal/mol), thus revealing that formation of these two hydrogen bonds is positively cooperative. In contrast, delta(AB)(i) is 0.9 +/- 0.4 kcal/mol for the poorly functional C27A ScYLV RNA, indicative of nonadditive hydrogen bond formation. These results reveal that cooperative hydrogen bond formation across the helical stem junction in H-type pseudoknots correlates with enhanced frameshift stimulation by luteoviral mRNA pseudoknots.

The Energetic Network of Hotspot Residues between Cdc25B Phosphatase and its Protein Substrate

We have investigated the functional network of hotspot residues at the remote docking site of two cell cycle regulators, namely Cdc25B phosphatase and its native protein substrate Cdk2-pTpY/CycA. Specifically, we have studied the roles of energetically important residues (Arg488, Arg492, Tyr497 on Cdc25B and Asp206 and Asp210 on Cdk2-pTpY/CycA) by generating a diverse set of substitutions and performing double and triple mutant cycle analyses. This transient protein–protein interaction is particularly well-suited for this mutagenic approach because various control experiments ensure that the effect of each mutation is limited to the interaction of interest. We find binary coupling energies for ion pairs and hydrogen bonds ranging from 0.7 kcal/mol to 3.9 kcal/mol and ternary coupling energies of 1.9 kcal/mol and 2.8 kcal/mol. Overall our biochemical analyses are in good agreement with the docked structure of the complex and suggest the following roles for the individual hotspot residues on Cdc25B. The most important contributor, Arg492, forms a specific and tight bidentate interaction with Asp206 and a weaker interaction with Asp210 that cannot be replaced by a Lys. Although Tyr497 does not directly participate in this ionic network, it is important for buttressing Arg492 using both its hydrophobic (aromatic ring) and hydrophilic characteristics (hydrogen bonding). Arg488 participates less specifically in the electrostatic network with Asp206 and Asp210 of the protein substrate as it can partially be replaced by Lys. Our data provide insight how a cluster of residues in a docking site remote from the site of the chemical reaction can bring about efficient and specific substrate recognition.

Improving the Stability of an Antibody Variable Fragment by a Combination of Knowledge-based Approaches: Validation and Mechanisms

Numerous approaches have been described to obtain variable fragments of antibodies (Fv or scFv) that are sufficiently stable for their applications. Here, we combined several knowledge-based methods to increase the stability of pre-existing scFvs by design. Firstly, the consensus sequence approach was used in a non-stringent way to predict a large basic set of potentially stabilizing mutations. These mutations were then prioritized by other methods of design, mainly the formation of additional hydrogen bonds, an increase in the hydrophilicity of solvent exposed residues, and previously described mutations in other antibodies. We validated this combined method with antibody mAbD1.3, directed against lysozyme. Fourteen potentially stabilizing mutations were designed and introduced into scFvD1.3 by site-directed mutagenesis, either individually or in combinations. We characterized the effects of the mutations on the thermodynamic stability of scFvD1.3 by experiments of unfolding with urea, monitored by spectrofluorometry, and tested the additivity of their effects by double-mutant cycles. We also quantified the individual contributions of the resistance to denaturation ([urea] 1/2) and cooperativity of unfolding ( m) to the variations of stability and the energy of coupling between mutations by a novel approach. Most mutations (75%) were stabilizing and none was destabilizing. The progressive recombination of the mutations into the same molecule of scFvD1.3 showed that their effects were mostly additive or synergistic, provided a large overall increase in protein stability (9.1 kcal/mol), and resulted in a highly stable scFvD1.3 derivative. The mechanisms of the mutations and of their combinations involved variations in the resistance to denaturation, cooperativity of unfolding, and likely residual structures of the denatured state, which was constrained by two disulfide bonds. This combined method should be applicable to any recombinant antibody fragment, through a single step of mutagenesis.