Conformational substates in enzyme mechanism: the 120 K structure of alpha-lytic protease at 1.5 A resolution.

Insight into the dynamic properties of alpha-lytic protease (alpha LP) has been obtained through the use of low-temperature X-ray crystallography and multiple-conformation refinement. Previous studies of alpha LP have shown that the residues around the active site are able to move significantly to accommodate substrates of different sizes. Here we show a link between the ability to accommodate ligands and the dynamics of the binding pocket. Although the structure of alpha LP at 120 K has B-factors with a uniformly low value of 4.8 A2 for the main chain, four regions stand out as having significantly higher B-factors. Because thermal motion should be suppressed at cryogenic temperatures, the high B-factors are interpreted as the result of trapped conformational substates. The active site residues that are perturbed during accommodation of different substrates are precisely those showing conformational substates, implying that substrate binding selects a subset of conformations from the ensemble of accessible states. To better characterize the precise nature of these substates, a protein model consisting of 16 structures has been refined and evaluated. The model reveals a number of features that could not be well-described by conventional B-factors: for example, 40% of the main-chain residue conformations are distributed asymmetrically or in discrete clusters. Furthermore, these data demonstrate an unexpected correlation between motions on either side of the binding pocket that we suggest is a consequence of "dynamic close packing." These results provide strong evidence for the role of protein dynamics in substrate binding and are consistent with the results of dynamic studies of ligand binding in myoglobin and ribonuclease A.     

The folding landscape of an alpha-lytic protease variant reveals the role of a conserved beta-hairpin in the development of kinetic stability

Most secreted bacterial proteases, including alpha-lytic protease (alphaLP), are synthesized with covalently attached pro regions necessary for their folding. The alphaLP folding landscape revealed that its pro region, a potent folding catalyst, is required to circumvent an extremely large folding free energy of activation that appears to be a consequence of its unique unfolding transition. Remarkably, the alphaLP native state is thermodynamically unstable; a large unfolding free energy barrier is solely responsible for the persistence of its native state. Although alphaLP folding is well characterized, the structural origins of its remarkable folding mechanism remain unclear. A conserved beta-hairpin in the C-terminal domain was identified as a structural element whose formation and positioning may contribute to the large folding free energy barrier. In this article, we characterize the folding of an alphaLP variant with a more favorable beta-hairpin turn conformation (alphaLP(beta-turn)). Indeed, alphaLP(beta-turn) pro region-catalyzed folding is faster than that for alphaLP. However, instead of accelerating spontaneous folding, alphaLP(beta-turn) actually unfolds more slowly than alphaLP. Our data support a model where the beta-hairpin is formed early, but its packing with a loop in the N-terminal domain happens late in the folding reaction. This tight packing at the domain interface enhances the kinetic stability of alphaLP(beta-turn), to nearly the same degree as the change between alphaLP and a faster folding homolog. However, alphaLP(beta-turn) has impaired proteolytic activity that negates the beneficial folding properties of this variant. This study demonstrates the evolutionary limitations imposed by the simultaneous optimization of folding and functional properties.     

Molecular mechanisms of resistance: free energy calculations of mutation effects on inhibitor binding to HIV-1 protease

The changes in the inhibitor binding constants due to the mutation of isoleucine to valine at position 84 of HIV-1 protease are calculated using molecular dynamics simulations. The calculations are done for three potent inhibitors--KNI-272, L-735,524 (indinavir or MK-639), and Ro 31-8959 (saquinavir). The calculations agree with the experimental data both in terms of an overall trend and in the magnitude of the resulting free energy change. HIV-1 protease is a homodimer, so each mutation causes two changes in the enzyme. The decrease in the binding free energy from each mutated side chain differs among the three inhibitors and correlates well with the size of the cavities induced in the protein interior near the mutated residue. The cavities are created as a result of a mutation to a smaller side chain, but the cavities are less than would be predicted from the wild-type structures, indicating that there is significant relaxation to partially fill the cavities.     

Linkers for improved cleavage of fusion proteins with an engineered alpha-lytic protease.

Addition of an N-terminal fusion partner can greatly aid the expression and purification of a recombinant protein in Escherichia coli. We investigated two genetically engineered proteases designed to remove the fusion partner after the protein of interest has been expressed. Recombinant human insulin-like growth factor-II (hIGF-II) has been produced from E. coli-derived fusion proteins using a novel enzymatic cleavage system that uses a mutant of alpha-lytic protease. Initially, two potential fusion protein linkers were designed, Pro-Ala-Pro-His (PAPH) and Pro-Ala-Pro-Met (PAPM), and were tested as substrates in the form of synthetic dodecapeptides. Using mass spectrometry and reverse-phase HPLC, the position of cleavage was confirmed and the kinetics of synthetic peptide cleavage were examined. Use of the linkers in hIGF-II fusion proteins produced in E. coli was then evaluated. The fusion proteins constructed consist of the first 11 amino acids of porcine growth hormone linked N-terminally to hIGF-II by six amino acids that include the dipeptide Val-Asn followed by a variable tetrapeptide protease cleavage motif. Mass spectrometry and N-terminal sequencing confirmed that proteolytic cleavage of the fusion proteins had occurred at the predicted sites. Using the fusion proteins as substrates, the cleavage of the rationally designed motifs by the alpha-lytic protease mutant was compared. The fusion protein containing the motif PAPM had a k(cat)/K(M) ratio indicating a 1.6-fold preference over the PAPH fusion protein for cleavage by this enzyme. Furthermore, when hIGF-II fusion proteins containing the designed cleavable linkers were processed with the engineered alpha-lytic protease, they gave greatly improved yields of native hIGF-II compared to an analogous fusion protein cleaved by H64A subtilisin. Comparison of the peptide and protein cleavage studies shows that the efficient proteolysis of the cleavage motifs is an inherent property of the designed sequences and is not determined by secondary or tertiary structure in the fusion proteins.     

Free energy calculations show that acidic P1 variants undergo large pKa shifts upon binding to trypsin

Serine proteinases and their protein inhibitors belong to one of the most comprehensively studied models of protein-protein interactions. It is well established that the narrow trypsin specificity is caused by the presence of a negatively charged aspartate at the specificity pocket. X-ray crystallography as well as association measurements revealed, surprisingly, that BPTI with glutamatic acid as the primary binding (P1) residue was able to bind to trypsin. Previous free energy calculations showed that there was a substantially unfavorable binding free energy associated with accommodation of ionized P1 Glu at the S1-site of trypsin. In this study, the binding of P1 Glu to trypsin has been systematically investigated in terms of the protonation states of P1 Glu and Asp189, the orientation of Gln192, as well as the possible presence of counterions using the linear interaction energy (LIE) approach and the free energy perturbation (FEP) method. Twenty-four conceivable binding arrangements were evaluated and quantitative agreement with experiments is obtained when the P1 Glu binds in its protonated from. The results suggest that P1 Glu is one of the variants of BPTI that inhibit trypsin strongest at low pH, contrary to the specificity profile of trypsin, suggesting a new regulation mechanism of trypsin-like enzymes.     

Disabling the folding catalyst is the last critical step in alpha-lytic protease folding

Alpha-Lytic protease (alphaLP) is an extracellular bacterial pro-protease marked by extraordinary conformational rigidity and a highly cooperative barrier to unfolding. Although these properties successfully limit its proteolytic destruction, thereby extending the functional lifetime of the protease, they come at the expense of foldability (t(1/2) = 1800 yr) and thermodynamic stability (native alphaLP is less stable than the unfolded species). Efficient folding has required the coevolution of a large N-terminal pro region (Pro) that rapidly catalyzes alphaLP folding (t(1/2) = 23 sec) and shifts the thermodynamic equilibrium in favor of folded protease through tight native-state binding. Release of active alphaLP from this stabilizing, but strongly inhibitory, complex requires the proteolytic destruction of Pro. alphaLP is capable of initiating Pro degradation via cleavage of a flexible loop within the Pro C-terminal domain. This single cleavage event abolishes Pro catalysis while maintaining strong native-state binding. Thus, the loop acts as an Achilles' heel by which the Pro foldase machinery can be safely dismantled, preventing Pro-catalyzed unfolding, without compromising alphaLP native-state stability. Once the loop is cleaved, Pro is rapidly degraded, releasing active alphaLP.     

Insights into the functional role of protonation states in the HIV-1 protease-BEA369 complex: molecular dynamics simulations and free energy calculations

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method combined with molecular dynamics (MD) simulations were used to investigate the functional role of protonation in human immunodeficiency virus type 1 (HIV-1) protease complexed with the inhibitor BEA369. Our results demonstrate that protonation of two aspartic acids (Asp25/Asp25′) has a strong influence on the dynamics behavior of the complex, the binding free energy of BEA369, and inhibitor–residue interactions. Relative binding free energies calculated using the MM-PBSA method show that protonation of Asp25 results in the strongest binding of BEA369 to HIV-1 protease. Inhibitor–residue interactions computed by the theory of free energy decomposition also indicate that protonation of Asp25 has the most favorable effect on binding of BEA369. In addition, hydrogen-bond analysis based on the trajectories of the MD simulations shows that protonation of Asp25 strongly influences the water-mediated link of a conserved water molecule, Wat301. We expect that the results of this study will contribute significantly to binding calculations for BEA369, and to the design of high affinity inhibitors.     

Enzyme specificity under dynamic control II: Principal component analysis of α-lytic protease using global and local solvent boundary conditions

The contributions of conformational dynamics to substrate specificity have been examined by the application of principal component analysis to molecular dynamics trajectories of alpha-lytic protease. The wild-type alpha-lytic protease is highly specific for substrates with small hydrophobic side chains at the specificity pocket, while the Met190-->Ala binding pocket mutant has a much broader specificity, actively hydrolyzing substrates ranging from Ala to Phe. Based on a combination of multiconformation analysis of cryo-X-ray crystallographic data, solution nuclear magnetic resonance (NMR), and normal mode calculations, we had hypothesized that the large alteration in specificity of the mutant enzyme is mainly attributable to changes in the dynamic movement of the two walls of the specificity pocket. To test this hypothesis, we performed a principal component analysis using 1-nanosecond molecular dynamics simulations using either a global or local solvent boundary condition. The results of this analysis strongly support our hypothesis and verify the results previously obtained by in vacuo normal mode analysis. We found that the walls of the wild-type substrate binding pocket move in tandem with one another, causing the pocket size to remain fixed so that only small substrates are recognized. In contrast, the M190A mutant shows uncoupled movement of the binding pocket walls, allowing the pocket to sample both smaller and larger sizes, which appears to be the cause of the observed broad specificity. The results suggest that the protein dynamics of alpha-lytic protease may play a significant role in defining the patterns of substrate specificity. As shown here, concerted local movements within proteins can be efficiently analyzed through a combination of principal component analysis and molecular dynamics trajectories using a local solvent boundary condition to reduce computational time and matrix size.     

Free energy barrier estimation of unfolding the alpha-helical surfactant-associated polypeptide C

Molecular dynamics simulations were conducted to estimate the free energy barrier of unfolding surfactant-associated polypeptide C (SP-C) from an alpha-helical conformation. Experimental studies indicate that while the helical fold of SP-C is thermodynamically stable in phospholipid micelles, it is metastable in a mixed organic solvent of CHCl3/CH3OH/0.1 M HCl at 32:64:5 (v/v/v), in which it undergoes an irreversible transformation to an insoluble aggregate that contains beta-sheet. On the basis of experimental observations, the free energy barrier was estimated to be approximately 100 kJ/mole by applying Eyring's transition state theory to the experimental rate of unfolding [Protein Sci 1998;7:2533-2540]. These studies prompted us to carry out simulations to investigate the unwinding process of two helical turns encompassing residues 25-32 in water and in methanol. The results give an upper bound estimation for the free energy barrier of unfolding of SP-C of approximately 20 kJ/mole. The results suggest a need to reconsider the applicability of a single-mode activated process theory to protein unfolding.     

Reaction path and free energy calculations of the transition between alternate conformations of HIV-1 protease

Two different structures of ligand-free HIV protease have been determined by X-ray crystallography. These structures differ in the position of two 12 residue, β-hairpin regions (or “flaps”) which cap the active site. The movements of the flaps must be involved in the binding of substrates since, in either conformation, the flaps block the binding site. One of these structures is similar to structures of the ligand-bound enzyme; however, the importance of both structures to enzyme function is unclear. This transformation takes place on a time scale too long for conventional molecular dynamics simulations, so the process was studied by first identifying a reaction path between the two structures and then calculating the free energy along this path using umbrella sampling. For the ligand-free enzyme, it is found that the two structures are nearly equally stable, with the ligand-bound-type structure being less stable, consistent with X-ray crystallography data. The more stable open structure does not have a lower potential energy, but is stabilized by entropy. The transition occurs through a collapse and reformation of the β-sheet structure of the conformationally flexible, glycine-rich flap ends. Additionally, some problems in studying conformational changes in proteins through the use of a single reaction path are addressed. Proteins 32:7–16, 1998. © 1998 Wiley-Liss, Inc.     

Probing structural requirements of fMLP receptor: on the size of the hydrophobic pocket corresponding to residue 2 of the tripeptide.

The conformationally constrained f-L-Met-Ac(n)c-L-Phe-OMe (n = 4,9-12) tripeptides, analogues of the chemoattractant f-L-Met-L-Leu-L-Phe-OH, were synthesized in solution by classical methods and fully characterized. These compounds and the published f-L-Met-Xxx-L-Phe-OMe (Xxx = Aib and Ac(n)c where n = 3, 5-8) analogues were compared to determine the combined effect of backbone preferred conformation and side-chain bulkiness at position 2 on the relation of 3D-structure to biological activity. A conformational study of all the analogues was performed in solution by FT-IR absorption and 1H-NMR techniques. In parallel, each peptide was tested for its ability to induce chemotaxis, superoxide anion production and lysozyme secretion from human neutrophils. The biological and conformational data are discussed in relation to the proposed model of the chemotactic receptor on neutrophils, in particular of the hydrophobic pocket accommodating residue 2 of the tripeptide.     

pKa calculations for class A beta-lactamases: influence of substrate binding

Beta-Lactamases are responsible for bacterial resistance to beta-lactams and are thus of major clinical importance. However, the identity of the general base involved in their mechanism of action is still unclear. Two candidate residues, Glu166 and Lys73, have been proposed to fulfill this role. Previous studies support the proposal that Glu166 acts during the deacylation, but there is no consensus on the possible role of this residue in the acylation step. Recent experimental data and theoretical considerations indicate that Lys73 is protonated in the free beta-lactamases, showing that this residue is unlikely to act as a proton abstractor. On the other hand, it has been proposed that the pKa of Lys73 would be dramatically reduced upon substrate binding and would thus be able to act as a base. To check this hypothesis, we performed continuum electrostatic calculations for five wild-type and three beta-lactamase mutants to estimate the pKa of Lys73 in the presence of substrates, both in the Henri-Michaelis complex and in the tetrahedral intermediate. In all cases, the pKa of Lys73 was computed to be above 10, showing that it is unlikely to act as a proton abstractor, even when a beta-lactam substrate is bound in the enzyme active site. The pKa of Lys234 is also raised in the tetrahedral intermediate, thus confirming a probable role of this residue in the stabilization of the tetrahedral intermediate. The influence of the beta-lactam carboxylate on the pKa values of the active-site lysines is also discussed.     

A poke in the eye: inhibiting HIV-1 protease through its flap-recognition pocket

A novel mechanism of inhibiting HIV-1 protease (HIVp) is presented. Using computational solvent mapping to identify complementary interactions and the Multiple Protein Structure method to incorporate protein flexibility, we generated a receptor-based pharmacophore model of the flexible flap region of the semiopen, apo state of HIVp. Complementary interactions were consistently observed at the base of the flap, only within a cleft with a specific structural role. In the closed, bound state of HIVp, each flap tip docks against the opposite monomer, occupying this cleft. This flap-recognition site is filled by the protein and cannot be identified using traditional approaches based on bound, closed structures. Virtual screening and dynamics simulations show how small molecules can be identified to complement this cleft. Subsequent experimental testing confirms inhibitory activity of this new class of inhibitor. This may be the first new inhibitor class for HIVp since dimerization inhibitors were introduced 17 years ago.     

Thrombin inhibitors with novel P1 binding pocket functionality: free energy of binding analysis

The high incidence of thrombembolic diseases justifies the development of new antithrombotics. The search for a direct inhibitor has resulted in the synthesis of a considerable number of low molecular weight molecules that inhibit human α-thrombin potently. However, efforts to develop an orally active drug remain in progress as the most active inhibitors with a highly basic P1 moiety exhibit an unsatisfactory bioavailability profile. In our previous work we solved several X-ray structures of human α-thrombin in complexes with (1) novel bicyclic arginine mimetics attached to the glycylproline amide and pyridinone acetamide scaffold and (2) inhibitors with a novel aza scaffold and with charged or neutral P1 moieties. In the present contribution, we correlate the structures of the complex between these inhibitors and the protein with the calculated free energy of binding. The energy of solvation was calculated using the Poisson–Boltzmann approach. In particular, the requirements for successful recognition of an inhibitor at the protein’s active site pocket S1 are discussed. Figure We report here on free energy of binding analysis of thrombin inhibitors with novel aza scaffold and novel bicyclic arginine mimetics in S1 pocket of thrombin     

Recombinant HIV1 protease secreted by Saccharomyces cerevisiae correctly processes myristylated gag polyprotein

The protease of the human immunodeficiency virus type I (HIV1) was expressed both intracellularly and extracellularly in Saccharomyces cerevisiae. Intracellular expression of the protease was achieved by fusing a 179 amino acid precursor form of the protease to human superoxide dismutase (hSOD). Self-processing of the viral enzyme from the hybrid precursor was demonstrated to occur within the yeast host. Secretion of the protease was achieved by fusing the leader sequence of yeast alpha-factor to the precursor form of the protease or to the 99 amino acid mature form of the protease. Authentic and active forms of the retroviral enzyme were detected in yeast supernatants of cells expressing the precursor or the mature form of the protease. A D25E active site variant of the retroviral enzyme exhibited diminished autocatalytic activity when expressed intracellularly or secreted from yeast. The wild-type protease was active in an in vitro assay on the natural substrate, myristylated gag precursor, Pr53gag. Correct processing of Pr53gag at the Tyr 138-Pro 139 junction was confirmed by amino terminal sequence analysis of the resulting capsid protein (CA, p24). The secreted protease was purified to homogeneity from yeast media using preparative isoelectric focusing and reverse-phase HPLC. Amino terminal sequence analysis showed a sequence beginning at amino acid 1 of the mature enzyme (Pro) and another sequence beginning at amino acid 6 (Trp). This shorter sequence may represent a natural autolytic product of the protease.     

Development and binding characteristics of phosphonate inhibitors of SplA protease from Staphylococcus aureus

Staphylococcus aureus is responsible for a variety of human infections, including life-threatening, systemic conditions. Secreted proteome, including a range of proteases, constitutes the major virulence factor of the bacterium. However, the functions of individual enzymes, in particular SplA protease, remain poorly characterized. Here, we report development of specific inhibitors of SplA protease. The design, synthesis, and activity of a series of α-aminoalkylphosphonate diaryl esters and their peptidyl derivatives are described. Potent inhibitors of SplA are reported, which may facilitate future investigation of physiological function of the protease. The binding modes of the high-affinity compounds Cbz-Phe^P-(OC₆H₄−4-SO₂CH₃)₂ and Suc-Val-Pro-Phe^P-(OC₆H₅)₂ are revealed by high-resolution crystal structures of complexes with the protease. Surprisingly, the binding mode of both compounds deviates from previously characterized canonical interaction of α-aminoalkylphosphonate peptidyl derivatives and family S1 serine proteases.     

Alpha1-microglobulin chromophores are located to three lysine residues semiburied in the lipocalin pocket and associated with a novel lipophilic compound

Alpha1-microglobulin (alpha1m) is an electrophoretically heterogeneous plasma protein. It belongs to the lipocalin superfamily, a group of proteins with a three-dimensional (3D) structure that forms an internal hydrophobic ligand-binding pocket. Alpha1m carries a covalently linked unidentified chromophore that gives the protein a characteristic brown color and extremely heterogeneous optical properties. Twenty-one different colored tryptic peptides corresponding to residues 88-94, 118-121, and 122-134 of human alpha1m were purified. In these peptides, the side chains of Lys92, Lys118, and Lys130 carried size heterogeneous, covalently attached, unidentified chromophores with molecular masses between 122 and 282 atomic mass units (amu). In addition, a previously unknown uncolored lipophilic 282 amu compound was found strongly, but noncovalently associated with the colored peptides. Uncolored tryptic peptides containing the same Lys residues were also purified. These peptides did not carry any additional mass (i.e., chromophore) suggesting that only a fraction of the Lys92, Lys118, and Lys130 are modified. The results can explain the size, charge, and optical heterogeneity of alpha1m. A 3D model of alpha1m, based on the structure of rat epididymal retinoic acid-binding protein (ERABP), suggests that Lys92, Lys118, and Lys130 are semiburied near the entrance of the lipocalin pocket. This was supported by the fluorescence spectra of alpha1m under native and denatured conditions, which indicated that the chromophores are buried, or semiburied, in the interior of the protein. In human plasma, approximately 50% of alpha1m is complex bound to IgA. Only the free alpha1m carried colored groups, whereas alpha1m linked to IgA was uncolored.     

Investigation on the binding mechanism of loratinib with the c-ros oncogene 1 (ROS1) receptor tyrosine kinase via molecular dynamics simulation and binding free energy calculations

The c-ros oncogene 1 (ROS1) has proven to be an important cancer target for the treatment of various human cancers. The anaplastic lymphoma kinase inhibitor crizotinib has been granted approval for the treatment of patients with ROS1 positive metastatic non-small-cell lung cancer by the Food and Drug Administration on 2016. However, serious resistance due to the secondary mutation of glycine 2032 to arginine (G2032R) was developed in clinical studies. Loratinib (PF-06463922), a macrocyclic analog of crizotinib, showed significantly improved inhibitory activity against wild–type (WT) ROS1 and ROS1^G2032R mutant. To provide insights into the inhibition mechanism, molecular dynamics simulations and free energy calculations were carried out for the complexes of loratinib with WT and G2032R mutated ROS1. The apo-ROS1^WT and apo-ROS1^G2032R systems showed similar RMSF distributions, while ROS1^G2032R-loratinib showed significantly higher than that of WT ROS1-loratinib, which revealed that the binding of loratinib to ROS1^G2032R significantly interfered the ?uctuation of protein. Calculations of binding free energies indicate that G2032R mutation significantly reduces the binding affinity of loratinib for ROS1, which arose mostly from the increase of conformation entropy and the decrease of solvation energy. Furthermore, detailed per-residue binding free energies highlighted the increased and decreased contributions of some residues in the G2032R mutated systems. The present study revealed the detailed inhibitory mechanism of loratinib as potent WT and G2032R mutated ROS1 inhibitor, which was expected to provide a basis for rational drug design.     

Anatomy of a conformational transition of beta-strand 6 in soybean beta-amylase caused by substrate (or inhibitor) binding to the catalytical site.

A computational study of the five soybean beta-amylase X-ray structure reported so far revealed a peculiar conformational transition after substrate (or inhibitor) binding, which affects a segment of the beta-strand 6 (residues 341-343) in the (beta/alpha)8 molecular scaffold. Backbone distortions that involve considerable changes in the phi and psi angles were observed, as well as two sharp rotamer transitions for the Thr342 and Cys343 side chains. These changes caused the outermost CA-layer (at the C-terminal side of the barrel), which is involved in the catalysis, to shrink. Our observations strongly suggest that the 341FTC343 residue conformations in the free enzyme are not optimal for protein stability. Furthermore, as a result of conformational transitions in the ligand-binding process, there is a negative enthalpy change for these residues (-27 and -34 kcal/mol, after substrate or inhibitor binding, respectively). These findings support the proposed "stability-function" hypothesis for proteins that recognize a ligand (Shoichet BK, Baase WA, Kuroki R, Matthews BW. 1995. A relationship between protein stability and protein function. Proc Natl Acad Sci USA 92:452-456). They are also in good agreement with other experimental results in the literature that describe the role of the 341-343 segment in beta-amylase activity. Site-directed mutagenesis focused on these residues could be useful for undertaking functional studies of beta-amylase.     

Study of the insulin dimerization: binding free energy calculations and per-residue free energy decomposition

A calculation of the binding free energy for the dimerization of insulin has been performed using the molecular mechanics-generalized Born surface area approach. The calculated absolute binding free energy is -11.9 kcal/mol, in approximate agreement with the experimental value of -7.2 kcal/mol. The results show that the dimerization is mainly due to nonpolar interactions. The role of the hydrogen bonds between the 2 monomers appears to give the direction of the interactions. A per-atom decomposition of the binding free energy has been performed to identify the residues contributing most to the self association free energy. Residues B24-B26 are found to make the largest favorable contributions to the dimerization. Other residues situated at the interface between the 2 monomers were found to make favorable but smaller contributions to the dimerization: Tyr B16, Val B12, and Pro B28, and to an even lesser extent, Gly B23. The energy decomposition on a per-residue basis is in agreement with experimental alanine scanning data. The results obtained from a single trajectory (i.e., the dimer trajectory is also used for the monomer analysis) and 2 trajectories (i.e., separate trajectories are used for the monomer and dimer) are similar.