Accurate Computation of Thermodynamic Activation Parameters in the Chorismate Mutase Reaction from Empirical Valence Bond Simulations.

Chorismate mutase (CM) enzymes have long served as model systems for benchmarking new methods and tools in computational chemistry. Despite the enzymes' prominence in the literature, the extent of the roles that activation enthalpy and entropy play in catalyzing the conversion of chorismate to prephenate is still subject to debate. Knowledge of these parameters is a key piece in fully understanding the mechanism of chorismate mutases. Within this study, we utilize EVB/MD free energy perturbation calculations at a range of temperatures, allowing us to extract activation enthalpies and entropies from an Arrhenius plot of activation free energies of the reaction catalyzed by a monofunctional Bacillus subtilis CM 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 CM reactions are enthalpically driven. Furthermore, we observe that the monofunctional CM 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.

Structure and mechanism of a phage-encoded SAM lyase revises catalytic function of enzyme family.

The first S-adenosyl methionine (SAM) degrading enzyme (SAMase) was discovered in bacteriophage T3, as a counter-defense against the bacterial restriction-modification system, and annotated as a SAM hydrolase forming 5'-methyl-thioadenosine (MTA) and L-homoserine. From environmental phages, we recently discovered three SAMases with barely detectable sequence similarity to T3 SAMase and without homology to proteins of known structure. Here, we present the very first phage SAMase structures, in complex with a substrate analogue and the product MTA. The structure shows a trimer of alpha-beta sandwiches similar to the GlnB-like superfamily, with active sites formed at the trimer interfaces. Quantum-mechanical calculations, thin-layer chromatography, and nuclear magnetic resonance spectroscopy demonstrate that this family of enzymes are not hydrolases but lyases forming MTA and L-homoserine lactone in a unimolecular reaction mechanism. Sequence analysis and in vitro and in vivo mutagenesis support that T3 SAMase belongs to the same structural family and utilizes the same reaction mechanism.

Bacteria can be infected by viruses known as bacteriophages. These viruses inject their genetic material into bacterial cells and use the bacteria's own machinery to build the proteins they need to survive and infect other cells. To protect themselves, bacteria produce a molecule called S-adenosyl methionine, or SAM for short, which deposits marks on the bacteria's DNA. These marks help the bacteria distinguish their own genetic material from the genetic material of foreign invaders: any DNA not bearing the mark from SAM will be immediately broken down by the bacterial cell. This system helps to block many types of bacteriophage infections, but not all. Some bacteriophages carry genes that code for enzymes called SAMases, which can break down SAM, switching off the bacteria's defenses. The most well-known SAMase was first discovered in the 1960s in a bacteriophage called T3. Chemical studies of this SAMase suggested that it works as a 'hydrolase', meaning that it uses water to break SAM apart. New SAMases have since been discovered in bacteriophages from environmental water samples, which, despite being able to degrade SAM, are genetically dissimilar to one another and the SAMase in T3. This brings into question whether these enzymes all use the same mechanism to break SAM down. To gain a better understanding of how these SAMases work, Guo, Söderholm, Kanchugal, Isaksen et al. solved the crystal structure of one of the newly discovered enzymes called Svi3-3. This revealed three copies of the Svi3-3 enzyme join together to form a unit that SAM binds to at the border between two of the enzymes. Computer simulations of this structure suggested that Svi3-3 holds SAM in a position where it cannot interact with water, and that once in the grip of the SAMase, SAM instead reacts with itself and splits into two. Experiments confirmed these predictions for Svi3-3 and the other tested SAMases. Furthermore, the SAMase from bacteriophage T3 was also found to degrade SAM using the same mechanism. This shows that this group of SAMases are not hydrolases as originally thought, but in fact 'lyases': enzymes that break molecules apart without using water. These findings form a starting point for further investigations into how SAM lyases help bacteriophages evade detection. SAM has various different functions in other living organisms, and these lyases could be used to modulate the levels of SAM in future studies investigating its role.

Towards Rational Computational Engineering of Psychrophilic Enzymes.

Cold-adapted enzymes from psychrophilic species achieve their high catalytic efficiency at low temperature by a different partitioning of the activation free energy into its enthalpic and entropic components, compared to orthologous mesophilic enzymes. Their lower activation enthalpy, partly compensated by an increased entropic penalty, has been suggested to originate from changes in flexibility of the protein surface. Multiple sequence alignments of psychrophilic and mesophilic enzymes also show characteristic motifs located in surface loops of the protein. Here, we use computer simulations to examine the effects of a number of designed surface mutations of psychrophilic and mesophilic elastases on the temperature dependence of the catalyzed peptide cleavage reaction. For each of 14 mutant enzyme variants we report calculations of their thermodynamic activation parameters. The results show that substitution of psychrophilic loop residues into the mesophilic enzyme consistently changes both the activation parameters and loop flexibilities towards the former, and vice versa for opposite substitutions.

QresFEP: An Automated Protocol for Free Energy Calculations of Protein Mutations in Q.

Predicting the effect of single-point mutations on protein stability or protein-ligand binding is a major challenge in computational biology. Free energy calculations constitute the most rigorous approach to this problem, though the estimation of converged values for amino acid mutations remains challenging. To overcome this limitation, we developed tailored protocols to calculate free energy shifts associated with single-point mutations. We herein describe the QresFEP protocol, which includes an extension of our recent protocols to cover all amino acids mutations, based on the latest versions of the OPLS-AA force field. QresFEP is implemented in an application programming interface framework and the graphic interface QGui, for the molecular dynamics software Q. The complete protocol is benchmarked in several model systems, optimizing a number of sampling parameters and the implementation of Zwanzig's exponential formula and Bennet's acceptance ratio methods. QresFEP shows an excellent performance on estimating the hydration free energies of amino acid side-chain mimics, including their charged analogues. We also examined its performance on a protein-ligand binding problem of pharmaceutical relevance, the antagonism of neuropeptide Y1 G protein-coupled receptor. Here, the calculations show very good agreement with the experimental effect of 16 mutations on the binding of antagonists BIBP3226, in line with our recent applications in this field. Finally, the characterization of 43 mutations of T4-lysozyme reveals the capacity of our protocol to assess variations of the thermal stability of proteins, achieving a similar performance to alternative free energy perturbation (FEP) approaches. In summary, QresFEP is a robust, versatile, and user-friendly computational FEP protocol to examine biochemical effects of single-point mutations with high accuracy.

Catalytic Adaptation of Psychrophilic Elastase.

The class I pancreatic elastase from Atlantic salmon is considered to be a cold-adapted enzyme in view of the cold habitat, the reduced thermostability of the enzyme, and the fact that it is faster than its mesophilic porcine counterpart at room temperature. However, no experimental characterization of its catalytic properties at lower temperatures has actually been reported. Here we use extensive computer simulations of its catalytic reaction, at different temperatures and with different peptide substrates, to compare its characteristics with those of porcine pancreatic elastase, with which it shares 67% sequence identity. We find that both enzymes have a preference for smaller aliphatic residues at the P1 position, while the reaction rate with phenylalanine at P1 is predicted to be substantially lower. With the former class of substrates, the calculated reaction rates for salmon enzyme are consistently higher than those of the porcine ortholog at all temperatures examined, and the difference is most pronounced at the lowest temperature. As observed for other cold-adapted enzymes, this is caused by redistribution of the activation free energy in terms of enthalpy and entropy and can be linked to differences in the mobility of surface-exposed loops in the two enzymes. Such mobility changes are found to be reflected by characteristic sequence conservation patterns in psychrophilic and mesophilic species. Hence, calculations of mutations in a single surface loop show that the temperature dependence of the catalytic reaction is altered in a predictable way.

A comparative study of cold- and warm-adapted Endonucleases A using sequence analyses and molecular dynamics simulations.

The psychrophilic and mesophilic endonucleases A (EndA) from Aliivibrio salmonicida (VsEndA) and Vibrio cholera (VcEndA) have been studied experimentally in terms of the biophysical properties related to thermal adaptation. The analyses of their static X-ray structures was no sufficient to rationalize the determinants of their adaptive traits at the molecular level. Thus, we used Molecular Dynamics (MD) simulations to compare the two proteins and unveil their structural and dynamical differences. Our simulations did not show a substantial increase in flexibility in the cold-adapted variant on the nanosecond time scale. The only exception is a more rigid C-terminal region in VcEndA, which is ascribable to a cluster of electrostatic interactions and hydrogen bonds, as also supported by MD simulations of the VsEndA mutant variant where the cluster of interactions was introduced. Moreover, we identified three additional amino acidic substitutions through multiple sequence alignment and the analyses of MD-based protein structure networks. In particular, T120V occurs in the proximity of the catalytic residue H80 and alters the interaction with the residue Y43, which belongs to the second coordination sphere of the Mg2+ ion. This makes T120V an amenable candidate for future experimental mutagenesis.

Entropy and Enzyme Catalysis.

The role played by entropy for the enormous rate enhancement achieved by enzymes has been debated for many decades. There are, for example, several confirmed cases where the activation free energy is reduced by around 10 kcal/mol due to entropic effects, corresponding to a rate enhancement of ∼107 compared to the uncatalyzed reaction. However, despite substantial efforts from both the experimental and theoretical side, no real consensus has been reached regarding the origin of such large entropic contributions to enzyme catalysis. Another remarkable instance of entropic effects is found in enzymes that are adapted by evolution to work at low temperatures, near the freezing point of water. These cold-adapted enzymes invariably show a more negative entropy and a lower enthalpy of activation than their mesophilic orthologs, which counteracts the exponential damping of reaction rates at lower temperature. The structural origin of this universal phenomenon has, however, remained elusive. The basic problem with connecting macroscopic thermodynamic quantities, such as activation entropy and enthalpy derived from Arrhenius plots, to the 3D protein structure is that the underlying detailed (microscopic) energetics is essentially inaccessible to experiment. Moreover, attempts to calculate entropy contributions by computer simulations have mostly focused only on substrate entropies, which do not provide the full picture. We have recently devised a new approach for accessing thermodynamic activation parameters of both enzyme and solution reactions from computer simulations, which turns out to be very successful. This method is analogous to the experimental Arrhenius plots and directly evaluates the temperature dependence of calculated reaction free energy profiles. Hence, by extensive molecular dynamics simulations and calculations of up to thousands of independent free energy profiles, we are able to extract activation parameters with sufficient precision for making direct comparisons to experiment. We show here that the agreement with the measured quantities, for both enzyme catalyzed and spontaneous solution reactions, is quite remarkable. Importantly, we can now address some of the most spectacular entropy effects in enzymes and clarify their detailed microscopic origin. Herein, we discuss as examples the conversion of cytidine to uridine catalyzed by cytidine deaminase and reactions taking place on the ribosome, namely, peptide bond formation and GTP hydrolysis by elongation factor Tu. It turns out that the large entropy contributions to catalysis in these cases can now be rationalized by our computational approach. Finally, we address the problem of cold adaptation of enzyme reaction rates and prove by computational experiments that the universal activation enthalpy-entropy phenomenon originates from mechanical properties of the outer protein surface.

Thermodynamics of the Purine Nucleoside Phosphorylase Reaction Revealed by Computer Simulations.

Enzymes are able to catalyze chemical reactions by reducing the activation free energy, yielding significant increases in the reaction rates. This can thermodynamically be accomplished by either reducing the activation enthalpy or increasing the activation entropy. The effect of remote mutations on the thermodynamic activation parameters of human purine nucleoside phosphorylase is examined using extensive molecular dynamics and free energy simulations. More than 2700 independent reaction free energy profiles for six different temperatures have been calculated to obtain high-precision computational Arrhenius plots. On the basis of these, the activation enthalpies and entropies were computed from linear regression of the plots with ΔG⧧ as a function of 1/T, and the obtained thermodynamic activation parameters are in very good agreement with those from experiments. The Arrhenius plots immediately show that the 6-oxopurines (INO and GUO) have identical slopes, whereas the 6-aminopurine (ADO) has a significantly different slope, indicating that the substrate specificity is related to the difference in thermodynamic activation parameters. Furthermore, the calculations show that the human PNP specificity for 6-oxopurines over 6-aminopurines originates from significant differences in electrostatic preorganization. The effect of the remote double mutation, K22E and H104R (E:R), has also been examined, as it alters human PNP toward the bovine PNP. These residues are situated on the protein surface, 28-35 Å from the active site, and the mutation alters the enthalpy-entropy balance with little effect on the catalytic rates. It is thus quite remarkable that the empirical valence bond method can reproduce the enthalpies and entropies induced by these long-range mutations.

Enzyme surface rigidity tunes the temperature dependence of catalytic rates.

The structural origin of enzyme adaptation to low temperature, allowing efficient catalysis of chemical reactions even near the freezing point of water, remains a fundamental puzzle in biocatalysis. A remarkable universal fingerprint shared by all cold-active enzymes is a reduction of the activation enthalpy accompanied by a more negative entropy, which alleviates the exponential decrease in chemical reaction rates caused by lowering of the temperature. Herein, we explore the role of protein surface mobility in determining this enthalpy-entropy balance. The effects of modifying surface rigidity in cold- and warm-active trypsins are demonstrated here by calculation of high-precision Arrhenius plots and thermodynamic activation parameters for the peptide hydrolysis reaction, using extensive computer simulations. The protein surface flexibility is systematically varied by applying positional restraints, causing the remarkable effect of turning the cold-active trypsin into a variant with mesophilic characteristics without changing the amino acid sequence. Furthermore, we show that just restraining a key surface loop causes the same effect as a point mutation in that loop between the cold- and warm-active trypsin. Importantly, changes in the activation enthalpy-entropy balance of up to 10 kcal/mol are almost perfectly balanced at room temperature, whereas they yield significantly higher rates at low temperatures for the cold-adapted enzyme.

Computer Simulations Reveal Substrate Specificity of Glycosidic Bond Cleavage in Native and Mutant Human Purine Nucleoside Phosphorylase.

Purine nucleoside phosphorylase (PNP) catalyzes the reversible phosphorolysis of purine ribonucleosides and 2'-deoxyribonucleosides, yielding the purine base and (2'-deoxy)ribose 1-phosphate as products. While this enzyme has been extensively studied, several questions with respect to the catalytic mechanism have remained largely unanswered. The role of the phosphate and key amino acid residues in the catalytic reaction as well as the purine ring protonation state is elucidated using density functional theory calculations and extensive empirical valence bond (EVB) simulations. Free energy surfaces for adenosine, inosine, and guanosine are fitted to ab initio data and yield quantitative agreement with experimental data when the surfaces are used to model the corresponding enzymatic reactions. The cognate substrates 6-aminopurines (inosine and guanosine) interact with PNP through extensive hydrogen bonding, but the substrate specificity is found to be a direct result of the electrostatic preorganization energy along the reaction coordinate. Asn243 has previously been identified as a key residue providing substrate specificity. Mutation of Asn243 to Asp has dramatic effects on the substrate specificity, making 6-amino- and 6-oxopurines equally good as substrates. The principal effect of this particular mutation is the change in the electrostatic preorganization energy between the native enzyme and the Asn243Asp mutant, clearly favoring adenosine over inosine and guanosine. Thus, the EVB simulations show that this particular mutation affects the electrostatic preorganization of the active site, which in turn can explain the substrate specificity.

Qgui: A high-throughput interface for automated setup and analysis of free energy calculations and empirical valence bond simulations in biological systems.

Structural information and activity data has increased rapidly for many protein targets during the last decades. In this paper, we present a high-throughput interface (Qgui) for automated free energy and empirical valence bond (EVB) calculations that use molecular dynamics (MD) simulations for conformational sampling. Applications to ligand binding using both the linear interaction energy (LIE) method and the free energy perturbation (FEP) technique are given using the estrogen receptor (ERα) as a model system. Examples of free energy profiles obtained using the EVB method for the rate-limiting step of the enzymatic reaction catalyzed by trypsin are also shown. In addition, we present calculation of high-precision Arrhenius plots to obtain the thermodynamic activation enthalpy and entropy with Qgui from running a large number of EVB simulations.

Protein surface softness is the origin of enzyme cold-adaptation of trypsin.

Life has effectively colonized most of our planet and extremophilic organisms require specialized enzymes to survive under harsh conditions. Cold-loving organisms (psychrophiles) express heat-labile enzymes that possess a high specific activity and catalytic efficiency at low temperatures. A remarkable universal characteristic of cold-active enzymes is that they show a reduction both in activation enthalpy and entropy, compared to mesophilic orthologs, which makes their reaction rates less sensitive to falling temperature. Despite significant efforts since the early 1970s, the important question of the origin of this effect still largely remains unanswered. Here we use cold- and warm-active trypsins as model systems to investigate the temperature dependence of the reaction rates with extensive molecular dynamics free energy simulations. The calculations quantitatively reproduce the catalytic rates of the two enzymes and further yield high-precision Arrhenius plots, which show the characteristic trends in activation enthalpy and entropy. Detailed structural analysis indicates that the relationship between these parameters and the 3D structure is reflected by significantly different internal protein energy changes during the reaction. The origin of this effect is not localized to the active site, but is found in the outer regions of the protein, where the cold-active enzyme has a higher degree of softness. Several structural mechanisms for softening the protein surface are identified, together with key mutations responsible for this effect. Our simulations further show that single point-mutations can significantly affect the thermodynamic activation parameters, indicating how these can be optimized by evolution.

Targeting the S1 and S3 subsite of trypsin with unnatural cationic amino acids generates antimicrobial peptides with potential for oral administration.

This study investigates how the S1 and S3 site of trypsin can be challenged with cationic amino acid analogues to yield active antimicrobial peptides with stability toward tryptic degradation. It is shown that unnatural analogues can be incorporated to generate stable peptides with maintained bioactivity to allow for a potential oral uptake. Selected peptides were studied using isothermal calorimetry and computational methods. Both stable and unstable peptides were found to bind stoichiometrically to trypsin with dissociation constants ranging 2-60 µM, suggesting several different binding modes. The stability of selected peptides was analyzed in whole organ extracts and the incorporation of homoarginine and 2-amino-(3-guanidino)propanoic acid resulted in a 14- and 50-fold increase in duodenal stability. In addition, a 40- and 70-fold increase in stomach stability is also reported. Overall, these results illustrate how the incorporation of cationic side chains can be employed to generate bioactive peptides with significant systemic stability.

In vitro characterization of human peptide transporter hPEPT1 interactions and passive permeation studies of short cationic antimicrobial peptides.

The present study assesses the permeation of cationic antimicrobial di- and tripeptides derived from lactoferricin via interaction with the human intestinal peptide transporter hPEPT1 and via passive routes. While some tested peptides displayed moderate affinity (0.6 and 2.7 mM) for interaction with hPEPT1, none served as substrate for hPEPT1 expressed by Xenopus laevis oocytes. It is shown that structural strategies employed to generate sufficient biological activity and metabolic stability such as introduction of large hydrophobic unnatural amino acids and different C-terminal modifications counteracted hPEPT1 mediated uptake. Most of the included peptides were nevertheless shown to permeate at rates suggesting moderate to excellent human oral absorption in the applied phospholipid vesicle-based passive permeation assay. Although the main factor governing passive permeation appears to be the hydrophobicity, peptide structure was also important and the overall permeation behavior was difficult to predict. Comparisons with a theoretical prediction model were also performed.

Unnatural amino acid side chains as S1, S1', and S2' probes yield cationic antimicrobial peptides with stability toward chymotryptic degradation.

This work describes how the systematic incorporation of a range of unnatural amino acid derivatives in the P1, P1', and P2' positions allows for the generation of short lactoferricin based cationic antimicrobial peptides with a stability toward chymotryptic degradation. The necessary pharmacophore sets the peptides up for degradation by chymotrypsin, and a heavily truncated native tripeptide was rapidly digested despite its short sequence. Degradation studies indicated that increased half-lives could be obtained by altering the binding to each subsite surrounding the active site without sacrificing the antimicrobial activity. Important structural and mechanistic features were revealed in a fashion not feasible through the use of native peptide substrates. The results, which are generally applicable to a range of relevant peptides, further show that not only the S1 pocket, but also to the notably less studied S1' site can be used to control the proteolytic stability by incorporating different analogues of tryptophan and arginine.