<p>Peptide bond formation and translational termination on the ribosome have been simulated by molecular mechanics, free energy perturbation, empirical valence bond (MD/FEP/EVB) and automated docking methods. Recent X-ray crystallographic data is used here to calculate the entire free energy surface for the system complete with substrates, ribosomal groups, solvent molecules and ions. A reaction mechanism for peptide bond formation emerges that is found to be catalyzed by the ribosome, in agreement with kinetic data and activation entropy measurements. The results show a water mediated network of hydrogen bonds, capable of reducing the reorganization energy during peptidyl transfer. The predicted hydrogen bonds and the structure of the active site were later confirmed by new X-ray structures with proper transition states analogs. </p><p>Elongation termination on the ribosome is triggered by binding of a release factor (RF) protein followed by rapid release of the nascent peptide. The structure of the RF, bound to the ribosomal peptidyl transfer center (PTC), has not been resolved in atomic detail. Nor is the mechanism known, by which the hydrolysis proceeds. Using automated docking of a hepta-peptide RF fragment, containing the highly conserved GGQ motif, we identified a conformation capable of catalyzing peptide hydrolysis. The MD/FEP/EVB calculations also reproduce the slow spontaneous release when RF is absent, and rationalize available mutational data. The network of hydrogen bonds, the active site structure, and the reaction mechanism are found to be very similar for both peptidyl transfer and termination. </p><p>New structural data, placing a ribosomal protein (L27) in the PTC, motivated additional MD/FEP/EVB simulations to determine the effect of this protein on peptidyl transfer. The simulations predict that the protein N terminus interacts with the A-site substrate in a way that promotes binding. The catalytic effect of L27 in the ribosome, however, is shown to be marginal and it therefore seems valid to view the PTC as a ribozyme. Simulations with the model substrate puromycin (Pmn) predicts that protonation of the N terminus can reduce the rate of peptidyl transfer. This could explain the different pH-rate profiles measured for Pmn, compared to other substrates.</p>

Peptide bond formation and translational termination on the ribosome have been simulated by molecular mechanics, free energy perturbation, empirical valence bond (MD/FEP/EVB) and automated docking methods. Recent X-ray crystallographic data is used here to calculate the entire free energy surface for the system complete with substrates, ribosomal groups, solvent molecules and ions. A reaction mechanism for peptide bond formation emerges that is found to be catalyzed by the ribosome, in agreement with kinetic data and activation entropy measurements. The results show a water mediated network of hydrogen bonds, capable of reducing the reorganization energy during peptidyl transfer. The predicted hydrogen bonds and the structure of the active site were later confirmed by new X-ray structures with proper transition states analogs. Elongation termination on the ribosome is triggered by binding of a release factor (RF) protein followed by rapid release of the nascent peptide. The structure of the RF, bound to the ribosomal peptidyl transfer center (PTC), has not been resolved in atomic detail. Nor is the mechanism known, by which the hydrolysis proceeds. Using automated docking of a hepta-peptide RF fragment, containing the highly conserved GGQ motif, we identified a conformation capable of catalyzing peptide hydrolysis. The MD/FEP/EVB calculations also reproduce the slow spontaneous release when RF is absent, and rationalize available mutational data. The network of hydrogen bonds, the active site structure, and the reaction mechanism are found to be very similar for both peptidyl transfer and termination. New structural data, placing a ribosomal protein (L27) in the PTC, motivated additional MD/FEP/EVB simulations to determine the effect of this protein on peptidyl transfer. The simulations predict that the protein N terminus interacts with the A-site substrate in a way that promotes binding. The catalytic effect of L27 in the ribosome, however, is shown to be marginal and it therefore seems valid to view the PTC as a ribozyme. Simulations with the model substrate puromycin (Pmn) predicts that protonation of the N terminus can reduce the rate of peptidyl transfer. This could explain the different pH-rate profiles measured for Pmn, compared to other substrates.

Phosphoryl transfers are essential chemical reactions in key life processes, including energy production, signal transduction and protein synthesis. They are known for having extremely low reaction rates in aqueous solution, reaching the scale of millions of years. In order to make life possible, enzymes that catalyse phosphoryl transfer, phosphoryl transferases, have evolved to be tremendously proficient catalysts, increasing reaction rates to the millisecond timescale. Due to the nature of the electronic structure of phosphorus atoms, understanding how hydrolysis of phosphate esters occurs is a complex task. Experimental studies on the hydrolysis of phosphate monoesters with acidic leaving groups suggest a concerted mechanism with a loose, metaphosphate-like transition state. Theoretical studies have suggested two possible concerted pathways, either with loose or tight transition state geometries, plus the possibility of a stepwise mechanism with the formation of a phosphorane intermediate. Different pathways were shown to be energetically preferable depending on the acidity of the leaving group. Here we performed computational studies to revisit how this mechanistic shift occurs along a series of aryl phosphate monoesters, suggesting possible factors leading to such change. The fact that distinct pathways can occur in solution could mean that the same is possible for an enzyme active site. We performed simulations on the catalytic activity of β-phosphoglucomutase, suggesting that it is possible for two mechanisms to occur at the same time for the phosphoryl transfer. Curiously, several phosphoryl transferases were shown to be able to catalyse not only phosphate ester hydrolysis, but also the cleavage of other compounds. We modeled the catalytic mechanism of two highly promiscuous members of the alkaline phosphatase superfamily. Our model reproduces key experimental observables and shows that these enzymes are electrostatically flexible, employing the same set of residues to enhance the rates of different reactions, with different electrostatic contributions per residue.

Recent years have seen tremendous developments in methods for computational modeling of (bio-) molecular systems. Ever larger reactive systems are being studied with high accuracy approaches, and high-level QM/MM calculations are being routinely performed. However, applying high-accuracy methods to large biological systems is computationally expensive and becomes problematic when conformational sampling is needed. To address this challenge, classical force field based approaches such as free energy perturbation (FEP) and empirical valence bond calculations (EVB) have been employed in this work. Specifically: Force-field independent metal parameters have been developed for a range of alkaline earth and transition metal ions, which successfully reproduce experimental solvation free energies, metal-oxygen distances, and coordination numbers. These are valuable for the computational study of biological systems. Experimental studies have shown that the epoxide hydrolase from Solanum tuberosum (StEH1) is not only an enantioselective enzyme, but for smaller substrates, displays enantioconvergent behavior. For StEH1, two detailed studies, involving combined experimental and computational efforts have been performed: We first used trans-stilbene oxide to establish the basic reaction mechanism of this enzyme. Importantly, a highly conserved and earlier ignored histidine was identified to be important for catalysis. Following from this, EVB and experiment have been used to investigate the enantioconvergence of the StEH1-catalyzed hydrolysis of styrene oxide. This combined approach involved wildtype StEH1 and an engineered enzyme variant, and established a molecular understanding of enantioconvergent behavior of StEH1. A novel framework was developed for the Computer-Aided Directed Evolution of Enzymes (CADEE), in order to be able to quickly prepare, simulate, and analyze hundreds of enzyme variants. CADEE’s easy applicability is demonstrated in the form of an educational example. In conclusion, classical approaches are a computationally economical means to achieve extensive conformational sampling. Using the EVB approach has enabled me to obtain a molecular understanding of complex enzymatic systems. I have also increased the reach of the EVB approach, through the implementation of CADEE, which enables efficient and highly parallel in silico testing of hundreds-to-thousands of individual enzyme variants.

Enantioselective reactions are one of the ways to produce pure chiral compounds. Understanding the basis of this selectivity makes it possible to guide enzyme design towards more efficient catalysts. One approach to study enzymes involved in chiral chemistry is through the use of computational models that are able to simulate the chemical reaction taking place. The potato epoxide hydrolase is one enzyme that is known to be both highly enantioselective, while still being robust upon mutation of residues to change substrate scope. The enzyme was used to investigate the epoxide hydrolysis mechanism for a number of different substrates, using the EVB approach to the reaction both in solution and in several enzyme variants. In addition to this, work has been performed on new ways of performing simulations of divalent transition metals, as well as development of new simulation software.

In the search for a sustainable and environmentally friendly energy source, artificial photosynthesis has been proposed  as a promising solution. Using water as a substrate, solar energy can be utilised to store energy in the chemical form of hydrogen fuel. In part of this global scientific effort, this thesis work focuses on enabling molecular dynamics simulations of a particular set of ruthenium centred  water oxidation catalysts. The new catalysts show great success because of a binuclear  reaction  pathway in aqueous solution which makes it very interesting to model and investigate. Utilising quantum  mechanical tools, a set of molecular mechanics force field parameters for Ru-involved bonds, angles, torsions, and partial charges was successfully obtained and examined. The work allows future large scale simulation of water  oxidation and oxygen evolution in order to gain understanding and improve artificial photosynthesis.

The reaction mechanisms for the hemoglobin degrading enzymes in the Plasmodium falciparum malaria parasite, plasmepsin II (Plm II) and histo-aspartic protease (HAP), have been analyzed by molecular simulations. The reaction free energy profiles, calculated by the empirical valence bond (EVB) method in combination with molecular dynamics (MD) and free energy perturbation (FEP) simulations are in good agreement with experimental data. Additional computational methods, such as homology modelling and automated substrate docking, were necessary to generate a 3D model and a reactive substrate conformation before the reaction mechanism in HAP could be investigated. HAP is found to be an aspartic protease with a peptide cleaving mechanism similar to plasmepsin II. The major difference between these enzymes is that the negatively charged tetrahedral intermediate is stabilized by the charged histidine in HAP while in Plm II it is a neutral aspartic acid. Also the reaction mechanism for two other aspartic proteases, cathepsin D and HIV-1 protease, was simulated. These enzymes are relevant both for the inhibitor selectivity and for obtaining a general picture of catalysis in aspartic proteases. Another project involves inhibitor design towards plasmepsins. In particular, Plm II directed inhibitors based on the dihydroxyethylene scaffold have been characterized computationally. Molecular dynamics (MD) simulations were used to propagate the investigated system through time and to generate ensembles used for the calculation of free energies. The ligand binding affinities were calculated with the linear interaction energy (LIE) method. The most potent inhibitor had a Ki value of 6 nM and showed 78 % parasite inhibition when tested on red blood cells infected by malaria parasite P. falciparum. Citrate synthase is part of the citric acid cycle and is present in organisms that live in cold sea water as well as hot springs. The temperature adaptation of citrate synthase to cold and heat was investigated in terms of the difference in transition state stabilization between the psychrophilic, mesophilic and hyperthermophilic homologues. The EVB, FEP and MD methods were used to generate reaction free energy profiles. The investigated energetics points toward the electrostatic stabilization during the reaction as the major difference between the different citrate synthase homologues. The electrostatic stabilization of the transition state is most effective in the following order of the citrate synthase homologues: hyperthermophile, mesophile, psycrophile. This could be a general rule for temperature adaptation of enzyme catalysis.

Ground state destabilization is a hypothesis to explain enzyme catalysis. The most popular interpretation of it is the entropic effect, which states that enzymes accelerate biochemical reactions by bringing the reactants to a favorable position and orientation and the entropy cost of this is compensated by enthalpy of binding. Once the enzyme-substrate complex is formed, the reaction could proceed with negligible entropy cost. Deamination of cytidine catalyzed by E.coli cytidine deaminase appears to agree with this hypothesis. In this reaction, the chemical transformation occurs with a negligible entropy cost and the initial binding occurs with a large entropy penalty that is comparable to the entropic cost of the uncatalyzed reaction. Our calculations revealed that this reaction occurs with different mechanisms in the cytidine deaminase and water. The uncatalyzed reaction involves a concerted mechanism and the entropy cost of this reaction appears to be dominated by the reacting fragments and first solvation shell. The catalyzed reaction occurs via a stepwise mechanism in which a hydroxide ion acts as the nucleophile. In the active site, the entropy cost of hydroxide ion formation is eliminated due to pre-organization of the active site. Hence, the entropic effect in this reaction is due to a pre-organized active site rather than ground state destabilization. In the second part of this thesis, we investigated peptide bond formation and peptidyl-tRNA hydrolysis at the peptidyl transferase center of the ribosome. Peptidyl-tRNA hydrolysis occurs by nucleophilic attack of a water molecule on the ester carbon of peptidyl-tRNA. Our calculations showed that this reaction proceeds via a base catalyzed mechanism where the A76 O2’ is the general base and activates the nucleophilic water. Peptide bond formation occurs by nucleophilic attack of the α-amino group of aminoacyl-tRNA on the ester carbon of peptidyl-tRNA. For this reaction we investigated two mechanisms: i) the previously proposed proton shuttle mechanism which involves a zwitterionic tetrahedral intermediate, and ii) a general base mechanism that proceeds via a negatively charged tetrahedral intermediate. Although both mechanisms resulted in reasonable activation energies, only the proton shuttle mechanism found to be consistent with the pH dependence of peptide bond formation.

Enolate formation by proton abstraction from an sp3-hybridized carbon atom situated next to a carbonyl or carboxylate group is an abundant process in nature. Since the corresponding nonenzymatic process in water is slow and unfavorable due to high intrinsic free energy barriers and high substrate pKa s, enzymes catalyzing such reaction steps must overcome both kinetic and thermodynamic obstacles. Computer simulations were used to study enolate formation catalyzed by glyoxalase I (GlxI) and 3-oxo-Δ5-steroid isomerase (KSI). The results, which reproduce experimental kinetic data, indicate that for both enzymes the free energy barrier reduction originates mainly from the balancing of substrate and catalytic base pKas. This was found to be accomplished primarily by electrostatic interactions. The results also suggest that the remaining barrier reduction can be explained by the lower reorganization energy in the preorganized enzyme compared to the solution reaction. Moreover, it seems that quantum effects, arising from zero-point vibrations and proton tunnelling, do not contribute significantly to the barrier reduction in GlxI. For KSI, the formation of a low-barrier hydrogen bond between the enzyme and the enolate, which is suggested to stabilize the enolate, was investigated and found unlikely. The low pKa of the catalytic base in the nonpolar active site of KSI may possibly be explained by the presence of a water molecule not detected by experiments. The hemoglobin-degrading aspartic proteases plasmepsinI and plasmepsin II from Plasmodium falciparum have emerged as putative drug targets against malaria. A series of C2- symmetric compounds with a 1,2-dihydroxyethylene scaffold were investigated for plasmepsin affinity, using computer simulations and enzyme inhibition assays. The calculations correctly predicted the stereochemical preferences of the scaffold and the effect of chemical modifications. Calculated absolute binding free energies reproduced experimental data well. As these inhibitors have down to subnanomolar inhibition constants of the plasmepsins and no measurable affinity to human cathepsin D, they constitute promising lead compounds for further drug development.