Die Bestimmung von freien Energieunterschieden ist essentiell für die Untersuchung von zahlreichen Prozessen, wie Wirkmechanismen von Medikamenten, den Verlauf von enzymatischen Reaktionen oder die Löslichkeit von Chemikalien. Mittels Molekulardynamiksimulationen sind freie Energierechnungen in der Lage freie Energieunterschiede mit hoher Genauigkeit zu bestimmen. Diese Genauigkeit ist
jedoch mit einem gewaltigen Rechenaufwand verbunden, sodass eine Bestimmung oft Tage oder Wochen dauert. Aus diesem Grunde werden erhebliche Anstrengungen zur Optimierung dieser Techniken unternommen.
Der erste Teil dieser Dissertation beschreibt die Anwendung der Bennett's Acceptance Ratio Methode (BAR) auf Probleme, wo konventionelle freie Energierechnungen nicht anwendbar sind. Dies illustriert die Vielseitigkeit dieser Methode. Desweiteren demonstrieren wir eine Erweiterung von BAR zur Behandlung von gebiasten Simulationen, die nicht der klassischen Boltzmann-Verteilung gehorchen. Wir bezeichnen diese Methode ergo als Non-Boltzmann Bennett (NBB). Im Rahmen von einigen praktischen Anwendungen wird anschliessend gezeigt, wie eine kreative Wahl des gebiasten Zustands die Effizienz von freien Energierechnungen erhöhen kann.
Im zweiten Teil werden BAR und NBB zur Bestimmung von Hydrationsenergien verwendet. Speziell in der Proteinfaltung oder beim Binden von Liganden spielen die Energiekosten für die (De-)Solvatation eine erhebliche Rolle. Unglücklicherweise können die Hydrationsenergien von Aminosäuren nicht experimentell bestimmt werden, weshalb oft auf Schätzungen mittels Seitenkettenanaloga zurückgegriffen wird. Die Annahme, dass Seitenketten repräsentativ für volle Aminosäuren sind, ist jedoch bisher nicht wissenschaftlich getestet worden. Daher wurden die Hydrationsenergien sowohl von Aminosäuren, als auch von Seitenkettenanaloga bestimmt. Es zeigt sich dabei eine erhebliche Diskrepanz, was sich auf zwei Effekte zurückführen läßt: Solventexklusion und Selbstsolvatation. Während beim Erstem der Zugang zum Lösungsmittel sterisch versperrt wird, ensteht zweiteres durch Wechselwirkungen des Aminosäurerückgrats mit polaren Gruppen der Seitenkette. Da viele Techniken in der computergestützten Chemie Selbstlösung nicht berücksichtigen, hat dies schwere Auswirkungen auf die Genauigkeit. Wir illustrieren dies anhand von impliziten Solventmethoden und diskutieren den Einfluss unserer Ergebnisse auf Proteinstudien.The determination of free energy differences is fundamental to the study of several processes such as the binding of drugs to proteins, the paths of enzymatic reactions or the solubility of chemical compounds. By employing molecular dynamics simulations, free energy calculations are capable to compute such free energy differences with high accuracy. However, this accuracy comes at excessive computational costs, often requiring days or weeks to obtain exact results. Thus, considerable effort still has to be invested in the optimization of such techniques.
The first half of this dissertation focuses on the application of Bennett's Acceptance Ratio method (BAR) to problems where standard methods to compute free energy differences are not feasible. This highlights the unique versatility of BAR. Furthermore, we demonstrate how to extend BAR in order to make use of non-Boltzmann probability distributions in biased simulations. We refer to this method as Non-Boltzmann Bennett (NBB). The NBB method is illustrated by several examples that demonstrate how a creative choice of the biased state can also improve the efficiency of free energy simulations.
The second half is concerned with the application of BAR and NBB to the study of hydration free energies. Especially in protein folding or ligand binding (de)-solvation penalties can contribute considerably to the free energy difference. Unfortunately, hydration free energies of amino acids cannot be measured experimentally. Thus, approximations based on side chain analog data are used instead. However, the assumption that side chain analogs are representative for full amino acids has never been thoroughly tested. We, therefore, computed both relative and absolute solvation free energies of amino acids and side chain analogs, showing that the results can deviate considerably due to two effects: Solvent exclusion and self-solvation. While the former accounts for the reduction of solute--solvent interactions due to sterical occlusions, the latter arises from interactions between the backbone and the polar functional groups of the side chains. Since several techniques in computational chemistry do not account for self-solvation, this finding has severe consequences. We illustrate this for several implicit solvent models and briefly discuss the implication of our results for the field of protein science