Towards Accurate Computations of Cofactor-Containing Biosystems

Abstract

In the cell, cofactor proteins take on a wide range of tasks. Photosystem II, for example, generates oxygen during photosynthesis and channelrhodopsin from the alga Chlamydomonas reinhardtii makes the organism move in response to light. For these proteins to function properly, the presence of their cofactor molecules is indispensable. Theoretical attempts to studying cofactor- containing systems like photosystem II might struggle with the lack of reliable force field parameters. For a retinal cofactor, on the other hand, quantum mechanical description of its electrostatic interactions is necessary to correctly characterize its behaviour inside of the protein. Here, we used two different theoretical approaches to study cofactors of two proteins: classical mechanics and force field optimization to describe the iron- containing cofactors of photosystem II, and combined quantum and classical mechanics to analyse the environmental effects on proton transfer energetics in channelrhodopsin. We successfully derived new CHARMM force field parameters for the hæm and non-hæm iron complexes of photosystem II. By comparing to quantum mechanical data and by using test simulations, we showed that these new parameters provide a greatly improved description of intramolecular and non-bonded interactions. The parameters presented here will facilitate reliable all-atom simulations of proteins that contain hæm and non-hæm iron complexes. In channelrhodopsin, deprotonation of the Schiff base only occurs after the light-induced isomerization of the retinal. We used the Weighted Histogram Approach Method to generate proton transfer pathways to understand what prevents retinal deprotonation in the dark state of the protein and identified three important determinants that ensure the chromophore’s continued protonation. These determinants were the presence of a positively charged lysine (K132), water molecules in the active site and the protonation state of the Schiff base counterion glutamate-162. Ultimately, we demonstrated that a combined quantum and classical mechanical approach can be applied to proteins embedded in a hydrated lipid membrane, and we created a protocol that can be transferred to other proteins to study their proton transfer pathways and energetics while taking the protein and lipid environment into account.

Kofaktorproteine übernehmen innerhalb der Zelle eine Vielfalt von Aufgaben. So generiert Photosystem II Sauerstoff im Rahmen der Photosynthese und Kanalrhodopsin kontrolliert phototaxische Reaktionen der Alge Chlamydomonas reinhardtii. Um die korrekte Funktionsweise solcher Proteine zu garantieren, ist das Vorhandensein von bestimmten Kofaktormolekülen unabdingbar. Versuche, Kofaktorsysteme wie Photosystem II zu untersuchen, können daran scheitern, dass keine geeigneten Kraftfeldparameter zur Verfügung stehen. Für Retinal als Kofaktor wiederum ist eine quantenmechanische Beschreibung der elektrostatischen Interaktionen notwendig, um zu gewährleisten, dass sein Verhalten innerhalb des Proteins korrekt dargestellt wird. In dieser Doktorarbeit wurden zwei verschiedene theoretische Ansätze angewandt, um Kofaktoren zweier Proteine zu untersuchen: klassische Mechanik und Kraftfeldoptimierung zur Beschreibung der eisenhaltigen Kofaktoren von Photosystem II sowie kombinierte Quantenmechanik/Molekularmechanik zur Analyse der Auswirkungen, die die Proteinumgebung auf den Protonentransfer in Kanalrhodopsin hat. Wir konnten neue CHARMM -Kraftfeldparameter für die Häm- und Nichthämeisenkomplexe in Photosystem II herleiten. Indem wir quantenmechanische Daten zum Vergleich heranzogen und Testsimulation nutzten, konnten wir zeigen, dass die neuen Parameter eine stark verbesserte Beschreibung der intra- und intermolekularen Interaktion bieten und zuverlässige Ganz-Atom-Simulationen von Proteinen ermöglichen, die Häm- und Nichthämeisenkomplexe enthalten. In Kanalrhodopsin deprotoniert die Schiffsche Base erst, nachdem die lichtinduzierte Isomerisierung des Retinals eingetreten ist. Wir haben den gewichteten Histogramm-Ansatz (Weighted Histogram Approach Method) genutzt, um Protonentransferpfade herzuleiten, um zu verstehen, welche Faktoren die Deprotonierung im Grundzustand des Proteins verhindern. Wir konnten drei wichtige Bedingungsfaktoren identifiziere: das Vorhandensein eines positiv geladenen Lysins (K132), Wassermoleküle im aktiven Zentrum sowie der Protonierungszustand von Glutamat-162, einem Gegenion der Schiffschen Base. Letzten Endes konnten wir demonstrieren, dass ein kombiniert quantenmechanisch/molekularmechanischer Ansatz zur Beschreibung von Proteinen in einer hydratisierten Lipidmembran genutzt werden kann. Außerdem haben wir ein Protokoll entwickelt, welches zur Untersuchung von Protonentransferpfaden und -energien unter Berücksichtigung der Protein- und Lipidumgebung herangezogen werden kann.