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.