Probing the Effect of Point Mutations at Protein-Protein Interfaces with Free Energy Calculations

We have studied the effect of point mutations of the primary binding residue (P1) at the protein-protein interface in complexes of chymotrypsin and elastase with the third domain of the turkey ovomucoid inhibitor and in trypsin with the bovine pancreatic trypsin inhibitor, using molecular dynamics simulations combined with the linear interaction energy (LIE) approach. A total of 56 mutants have been constructed and docked into their host proteins. The free energy of binding could be reliably calculated for 52 of these mutants that could unambiguously be fitted into the binding sites. We find that the predicted binding free energies are in very good agreement with experimental data with mean unsigned errors between 0.50 and 1.03 kcal/mol. It is also evident that the standard LIE model used to study small drug-like ligand binding to proteins is not suitable for protein-protein interactions. Three different LIE models were therefore tested for each of the series of protein-protein complexes included, and the best models for each system turn out to be very similar. The difference in parameterization between small drug-like compounds and protein point mutations is attributed to the preorganization of the binding surface. Our results clearly demonstrate the potential of free energy calculations for probing the effect of point mutations at protein-protein interfaces and for exploring the principles of specificity of hot spots at the interface

Accurate computation of thermodynamic activation parameters in the chorismate mutase reaction from empirical valence bond simulations

Chorismate mutase enzymes have long served as model systems for benchmarking new methods and tools in computational chemistry. Despite the enzymes’ prominence in literature, the extent of the roles activation enthalpy and entropy play in catalyzing the conversion of chorismite to prephenate is still subject to debate. Knowledge of these parameters is a key piece in fully understanding the mechanism of chorismite mutases. Within this study, we utilize EVB/MD free energy perturbations at a range of temperatures, allowing us to extract activation enthalpies and entropies from an Arrhenius plot of reaction free energies of the reaction catalyzed by a monofunctional Bacillus subtilis chorismate mutase and the promiscuous enzyme isochorismate pyruvate lyase of Pseudomonas aeruginosa. In comparison to the uncatalyzed reaction, our results show that both enzyme-catalyzed reactions exhibit a substantial reduction in activation enthalpy, while the effect on activation entropy is relatively minor, demonstrating that enzyme-catalyzed chorismate mutase reactions are enthalpically-driven. Furthermore, we observe that the monofunctional chorismate mutase from B. subtilis more efficiently catalyzes this reaction than its promiscuous counterpart. This is supported by a structural analysis of the reaction pathway at the transition state, from which we identified key residues explaining the enthalpically-driven nature of the reactions, and also the difference in efficiencies between the two enzymes

Structure and Mechanism of a Cold-Adapted Bacterial Lipase

The structural origin of enzyme cold-adaptation has been the subject of considerable research efforts in recent years. Comparative studies of orthologous mesophilic-psychrophilic enzyme pairs found in nature are an obvious strategy for solving this problem, but they often suffer from relatively low sequence identity of the enzyme pairs. Small bacterial lipases adapted to distinctly different temperatures appear to provide an excellent model system for these types of studies, as they may show a very high degree of sequence conservation. Here, we report the first crystal structures of lipase A from the psychrophilic bacterium Bacillus pumilus, which confirm the high structural similarity to the mesophilic Bacillus subtilis enzyme, as indicated by their 81% sequence identity. We further employ extensive QM/MM calculations to delineate the catalytic reaction path and its energetics. The computational prediction of a rate-limiting deacylation step of the enzymatic ester hydrolysis reaction is verified by stopped-flow experiments, and steady-state kinetics confirms the psychrophilic nature of the B. pumilus enzyme. These results provide a useful benchmark for examining the structural basis of cold-adaptation and should now make it possible to disentangle the effects of the 34 mutations between the two enzymes on catalytic properties and thermal stability

Trypsin specificity as elucidated by LIE calculations, X-ray structures, and association constant measurements

The variation in inhibitor specificity for five different amine inhibitors bound to CST, BT, and the cold-adapted AST has been studied by use of association constant measurements, structural analysis of high-resolution crystal structures, and the LIE method. Experimental data show that AST binds the 1BZA and 2BEA inhibitors 0.8 and 0.5 kcal/mole more strongly than BT. However, structural interactions and orientations of the inhibitors within the S1 site have been found to be virtually identical in the three enzymes studied. For example, the four water molecules in the inhibitor-free structures of AST and BT are channeled into similar positions in the S1 site, and the nitrogen atom(s) of the inhibitors are found in two cationic binding sites denoted Position1 and Position2. The hydrophobic binding contributions for all five inhibitors, estimated by the LIE calculations, are also in the same order (−2.1 ± 0.2 kcal/mole) for all three enzymes. Our hypothesis is therefore that the observed variation in inhibitor binding arises from different electrostatic interactions originating from residues outside the S1 site. This is well illustrated by AST, in which Asp 150 and Glu 221B, despite some distance from the S1 binding site, lower the electrostatic potential of the S1 site and thus enhance substrate binding. Because the trends in the experimentally determined binding energies were reproduced by the LIE calculations after adding the contribution from long-range interactions, we find this method very suitable for rational studies of protein–substrate interactions

Resolving the energy paradox of chaperone/usher-mediated fibre assembly

Periplasmic chaperone/usher machineries are used for assembly of filamentous adhesion organelles of Gram-negative pathogens in a process that has been suggested to be driven by folding energy. Structures of mutant chaperone–subunit complexes revealed a final folding transition (condensation of the subunit hydrophobic core) on the release of organelle subunit from the chaperone–subunit pre-assembly complex and incorporation into the final fibre structure. However, in view of the large interface between chaperone and subunit in the pre-assembly complex and the reported stability of this complex, it is difficult to understand how final folding could release sufficient energy to drive assembly. In the present paper, we show the X-ray structure for a native chaperone–fibre complex that, together with thermodynamic data, shows that the final folding step is indeed an essential component of the assembly process. We show that completion of the hydrophobic core and incorporation into the fibre results in an exceptionally stable module, whereas the chaperone–subunit pre-assembly complex is greatly destabilized by the high-energy conformation of the bound subunit. This difference in stabilities creates a free energy potential that drives fibre formation