Explicit influence of water microsolvation on charge transfer and dynamics in ground and excited electronic states of molecular systems
- Modern computational molecular quantum chemical studies, such as the present one, typically employ a wide range of theoretical techniques. The latter are often rather complicated and one should not generally expect that an experimental scientist in the area of physical chemistry, a potential reader of this work, should be familiar with all these techniques. To simplify the reading of the Thesis and to make it self-sufficient, it is supplied with an overview of the employed theoretical methodologies (Chapter 1). The overview explains basic quantum-chemical terminology referred to throughout the Thesis, introduces theoretical foundations of the methods and outlines their properties and limitations. In Part 1.1 of Chapter 1, methods for the solution of the molecular Schrödinger equation are introduced, while in the subsequent Parts 1.2 and 1.3 methods for the solution of the electronic Schrödinger equation are presented to find the ground and excited states, respectively. Part 1.4 is dedicated to basis-set effects which are omnipresent in electronic-structure calculations. It contains a number of unusual insights and concepts proposed by the author and, thus, may be insightful also to experts in quantum chemistry.
In Chapter 2, the phenomenon of acetone-water proton exchange catalyzed by tubular as well as amorphous aggregates of calixhydroquinone (CHQ) macromolecules, which has been observed previously in NMR experiments (Ref. D1D), is investigated by means of correlated quantum-chemical methods. The first part of the study (Section 2.3-2.7) considers concerted proton transfer, assisted by several initially neutral OH-groups in the hydrogen-bonded networks of CHQ aggregates. The second part of the study (Section 2.8-2.13) is dedicated to a second mechanism of proton exchange: step-wise proton transfer via formation of ionic intermediates resulting from CHQ pre-dissociation. CHQ application-specific as well as general conclusions, relevant to the main topic of the Thesis (i.e. influence of specific microsolvation on the considered proton transfer processes), are presented in Section 2.14.
The phenomenon of dual fluorescence observed in clusters of methyl 4-N,N-dimethylaminobenzoate ester (DMABME) and two water molecules in the gas phase, is studied in Chapter 3. Experimentally, the dual fluorescence was detected in experiments combining optical and ground-state ion-depletion infrared spectroscopies in ultracold molecular beams (Ref. D2D). In Section 3.3, calculated ground-state infrared spectra are presented that allow to identify the structures of those isomers, which are present in the gas-phase, as well as the structure of the isomer responsible for dual fluorescence. To further understand the reaction mechanism of dual fluorescence, excited-state potential energy surfaces of the identified isomers were computed along the relevant twisted intermolecular charge-transfer formation coordinate and the mechanism of energy dissipation in these complexes was investigated (Section 3.4-3.5) (Ref. D3D). A brief summary of the main results of this chapter and conclusions are given in Section 3.6. Finally, in Section 3.7 a complementary benchmark study of the quality of ground-state potential energy surfaces of prototypical hydrogen-bonded systems (ammonia-water and formic acid-water dimers) obtained at the level of BSSE-corrected MP2 combined with moderate basis sets, has been conducted. The quality of potential energy surfaces was tested with respect to basis-set size, level of electron correlation and anharmonicity effects and the applied methodology to identify the IR spectrum of hydrated DMABME complexes (Section 3.3) has been found to be sufficient to uniquely assign the IR spectra.
Development of decoration and preferential-etching methods for delineation of crystal defects in semiconductor materials
- Silicon wafers such as Silicon on Insulator (SOI) and strained silicon on Insulator (sSOI) are the essential and basic materials of advanced microelectronic devices. However, they often show various kinds of crystal defects which impair the function of these devices. The most efficient method to date, for detecting such defects and for determining their density, is to delineate them by etching the wafers with a suitable etching solution and characterise them via light optical microscopy. Etch pits are formed at defect sites which are etched at a faster rate than at the perfect lattice. The standard etching solution used for SOI and sSOI is a dilute version of Secco. As Secco contains carcinogenic and environmentally hazardous chromium (VI), the use of which is or will be restricted by law in many countries, suitable chromium (VI)-free etching solutions like Organic Peracid Etches (OPE), modified Chemical Polishing Etches (CP) like CP4 mod and mixtures with organic oxidizing agents like chloranil (CA) have been developed for the successful delineation of various types of crystal defects.
However there are still nanometer-sized defects which are hard to detect or escape detection by this method. Copper decoration is a well known method to magnify these defects. It consists in applying a copper nitrate solution to the back of the SOI or sSOI wafer. On annealing, copper diffuses through the substrate and the BOX (buried oxide) to the SOI/sSOI film and on quenching to room temperature, copper precipitates as copper silicide, SiCu3, foremost at crystal defects where the lattice strain is greater than at perfect lattice sites. These silicides increase the volume in these parts of the crystal lattice and defect magnification occurs. A considerable disadvantage of this method is its tendency for artefact formation, when the copper concentration used is too high, with the copper precipitating at the film surface. The consequence is a higher density of etch pits whereby true defect etch pits cannot be differentiated from those caused by artefacts.
The aim of this thesis is to show that the processes of decorating and etching can be combined successfully to delineate all crystal defects in SOI and sSOI. An ideal result would have been to find a copper decoration procedure that decorates all existing crystal defects at a copper concentration that avoids artefact formation.
Di-μ-bromido-bis-[(diethyl ether-κO)(2,4,6-tri-methyl-phen-yl)magnesium]: the mesityl Grignard reagent
- The crystal structure of the title compound, [Mg2Br2(C9H11)2(C4H10O)2], features a centrosymmetric two-centre magnesium complex with half a mol-ecule in the asymmetric unit. The Mg atom is in a considerably distorted Br2CO coordination. Bond lengths and angles are comparable with already published values. The crystal packing is stabilized by C-H⋯π inter-actions linking the complexes into sheets parallel to (0-11).
EGI user forum 2011 : book of abstracts
WeNMR: the tale of virtual research community in NMR and structural biology
Cell-free expression and molecular modeling of the γ-secretase complex and G-protein-coupled receptors
- Alzheimer’s disease (AD), which was first reported more than a century ago by Alhzeimer, is one of the commonest forms of dementia which affects >30 million people globally (>8 million in Europe). The origin and pathogenesis of AD is poorly understood and there is no cure available for the disease. AD is characterized by the accumulation of senile plaques composed of amyloid beta peptides (Ab 37-43) which is formed by the gamma secretase (GS) complex by cleaving amyloid precursor protein. Therefore GS can be an attractive drug target. Since GS processes several other substrates like Notch, CD44 and Cadherins, nonspecific inhibition of GS has many side effects. Due to the lack of crystal structure of GS, which is attributed to the extreme difficulties in purifying it, molecular modeling can be useful to understand its architecture. So far only low resolution cryoEM structures of the complex has been solved which only provides a rough structure of the complex at low 12-15 A resolution Furthermore the activity of GS in vitro can be achieved by means of cell-free (CF) expression.
GS comprises catalytic subunits namely presenilins and supporting elements containing Pen-2, Aph-1 and Nicastrin. The origin of AD is hidden in the regulated intramembrnae proteolysis (RIP) which is involved in various physiological processes and also in leukemia. So far growth factors, cytokines, receptors, viral proteins, cell adhesion proteins, signal peptides and GS has been shown to undergo RIP. During RIP, the target proteins undergo extracellular shredding and intramembrane proteolysis.
This thesis is based on molecular modeling, molecular dynamics (MD) simulations, cell-free (CF) expression, mass spectrometry, NMR, crystallization, activity assay etc of the components of GS complex and G-protein coupled receptors (GPCRs).
First I validated the NMR structure of PS1 CTF in detergent micelles and lipid bilayers using coarse-grained MD simulations using MARTINI forcefield implemented in Gromacs. CTF was simulated in DPC micelles, DPPC and DLPC lipid bilayer. Starting from random configuration of detergent and lipids, micelle and lipid bilyer were formed respectively in presence of CTF and it was oriented properly to the micelle and bilyer during the simulation. Around DPC molecules formed micelle around CTF in agreement of the experimental results in which 80-85 DPC molecules are required to form micelles. The structure obtained in DPC was similar to that of NMR structure but differed in bilayer simulations showed the possibility of substrate docking in the conserved PAL motif. Simulations of CTF in implicit membrane (IMM1) in CHAMM yielded similar structure to that from coarse grained MD.
I performed cell-free expression optimization, crystallization and NMR spectroscopy of Pen-2 in various detergent micelles. Additionally Pen-2 was modeled by a combination of rosetta membrane ab-initio method, HHPred distant homology modeling and incorporating NMR constraints. The models were validated by all atom and coarse grained MD simulations both in detergent micelles and POPC/DPPC lipid bilayers using MARTINI forcefield.
GS operon consisting of all four subunits was co-expressed in CF and purified. The presence of of GS subunits after pull-down with Aph-1 was determined by western blotting (Pen-2) and mass spectrometry (Presenilin-1 and Aph-1). I also studied interactions of especially PS1 CTF, APP and NTF by docking and MD.
I also made models and interfaces of Pen-2 with PS1 NTF and checked their stability by MD simulations and compared with experimental results. The goal is to model the interfaces between GS subunits using molecular modeling approaches based on available experimental data like cross-linking, mutations and NMR structure of C-terminal fragment of PS1 and transmembrane part of APP. The obtained interfaces of GS subunits may explain its catalysis mechanism which can be exploited for novel lead design. Due to lack of crystal/NMR structure of the GS subunits except the PS1 CTF, it is not possible to predict the effect of mutations in terms of APP cleavage. So I also developed a sequence based approach based on machine learning using support vector machine to predict the effect of PS1 CTF L383 mutations in terms of Aβ40/Aβ42 ratio with 88% accuracy. Mutational data derived from the Molgen database of Presenilin 1 mutations was using for training.
GPCRs (also called 7TM receptors) form a large superfamily of membrane proteins, which can be activated by small molecules, lipids, hormones, peptides, light, pain, taste and smell etc. Although 50% of the drugs in market target GPCRs , only few are targeted therapeutically. Such wide range of targets is due to involvement of GPCRs in signaling pathways related to many diseases i.e. dementia (like Alzheimer's disease), metabolic (like diabetes) including endocrinological disorders, immunological including viral infections, cardiovascular, inflammatory, senses disorders, pain and cancer.
Cannabinoid and adrenergic receptors belong to the class A (similar to rhodopsin) GPCRs. Docking of agonists and antagonists to CB1 and CB2 cannabinoid receptors revealed the importance of a centrally located rotamer toggle switch, and its possible role in the mechanism of agonist/antagonist recognition. The switch is composed of two residues, F3.36 and W6.48, located on opposite transmembrane helices TM3 and TM6 in the central part of the membranous domain of cannabinoid receptors. The CB1 and CB2 receptor models were constructed based on the adenosine A2A receptor template. The two best scored conformations of each receptor were used for the docking procedure. In all poses (ligand-receptor conformations) characterized by the lowest ligand-receptor intermolecular energy and free energy of binding the ligand type matched the state of the rotamer toggle switch: antagonists maintained an inactive state of the switch, whereas agonists changed it. In case of agonists of β2AR, the (R,R) and (S,S) stereoisomers of fenoterol, the molecular dynamics simulations provided evidence of different binding modes while preserving the same average position of ligands in the binding site. The (S,S) isomer was much more labile in the binding site and only one stable hydrogen bond was created. Such dynamical binding modes may also be valid for ligands of cannabinoid receptors because of the hydrophobic nature of their ligand-receptor interactions. However, only very long molecular dynamics simulations could verify the validity of such binding modes and how they affect the process of activation.
Human N-formyl peptide receptors (FPRs) are G protein-coupled receptors (GPCRs) involved in many physiological processes, including host defense against bacterial infection and resolving inflammation. The three human FPRs (FPR1, FPR2 and FPR3) share significant sequence homology and perform their action via coupling to Gi protein. Activation of FPRs induces a variety of responses, which are dependent on the agonist, cell type, receptor subtype, and also species involved. FPRs are expressed mainly by phagocytic leukocytes. Together, these receptors bind a large number of structurally diverse groups of agonistic ligands, including N-formyl and nonformyl peptides of different composition, that chemoattract and activate phagocytes. For example, N-formyl-Met-Leu-Phe (fMLF), an FPR1 agonist, activates human phagocyte inflammatory responses, such as intracellular calcium mobilization, production of cytokines, generation of reactive oxygen species, and chemotaxis. This ligand can efficiently activate the major bactericidal neutrophil functions and it was one of the first characterized bacterial chemotactic peptides. Whereas fMLF is by far the most frequently used chemotactic peptide in studies of neutrophil functions, atomistic descriptions for fMLF-FPR1 binding mode are still scarce mainly because of the absence of a crystal structure of this receptor. Elucidating the binding modes may contribute to designing novel and more efficient non-peptide FPR1 drug candidates. Molecular modeling of FPR1, on the other hand, can provide an efficient way to reveal details of ligand binding and activation of the receptor. However, recent modelings of FPRs were confined only to bovine rhodopsin as a template.
To locate specific ligand-receptor interactions based on a more appropriate template than rhodopsin we generated the homology models of FPR1 using the crystal structure of the chemokine receptor CXCR4, which shares over 30% sequence identity with FPR1 and is located in the same γ branch of phylogenetic tree of GPCRs (rhodopsin is located in α branch). Docking and model refinement procedures were pursued afterward. Finally, 40 ns full-atom MD simulations were conducted for the Apo form as well as for complexes of fMLF (agonist) and tBocMLF (antagonist) with FPR1 in the membrane. Based on locations of the N- and C-termini of the ligand the FPR1 extracellular pocket can be divided into two zones, namely, the anchor and activation regions. The formylated M1 residue of fMLF bound to the activation region led to a series of conformational changes of conserved residues. Internal water molecules participating in extended hydrogen bond networks were found to play a crucial role in transmitting the agonist-receptor interactions. A mechanism of initial steps of the activation concurrent with ligand binding is proposed.
I accurately predicted the structure and ligand binding pose of dopamine receptor 3 (RMSD to the crystal structure: 2.13 Å) and chemokine receptor 4 (CXCR4, RMSD to the crystal structure 3.21 Å) in GPCR-Dock 2010 competition. The homology model of the dopamine receptor 3 was 8 th best overall in the competition.
Insights into a Lipid Regulator by Solid-State MAS NMR: Kinetic and Structure-Functional Studies on Diacylglycerol Kinase
Sandra Johanna Ullrich
- In this thesis the integral membrane protein diacylglycerol kinase (DAGK) from E.coli is investigated with solid-state NMR. The aim is to gain an insight into the enzyme’s mechanism through integration of kinetic, structural and dynamic data. The biological function of DAGK is the transfer of the γ-phosphate group from Mg*ATP to diacylglycerol (DAG) building phosphatidic acid (PA) as port of the membrane-derived oligosaccharide cycle[31,34]. Surprisingly, DAGK does not share structural or sequential similarities with other kinases. Typical sequence motives found in other kinases, which catalyze phosphoryl transfer reactions, are not found. In its physiological form DAGK is a homo-trimer with nine transmembrane helices, three catalytic centers and a size of 39.6 kDa.
First, the set-up of a real-time 31P MAS NMR experiment is shown. This experiment allows measuring in real-time the simultaneous ATP hydrolysis in the aqueous phase and lipid substrate phos-phorylation in the membrane phase with atomic resolution under magic angle spinning. After fast transfer of the sample into the NMR spectrometer the enzymatic reaction is started with a temperature jump. This approach of real-time MAS NMR in a dual-phase system was demonstrated for the lipid substrate analogs dioleoyl- (DOG) and dibutyrylglycerol (DBG), with a C8 and C4 aliphatic chain, respectively. The combination of 31P direct and cross polarization functions as a dynamic filter. In the 31P direct polarized experiment nuclei in both phases are detected, while in the 31P cross polar-ized experiment, only nuclei in the membrane phase are detected. Rates for substrate turnover, i.e. degradation of γP-, βP, αP-ATP and build-up of βP-, αP-ADP, free phosphate as side reaction, and PA are obtained, which reveal a Michaelis-Menten behavior with regard to Mg*ATP and DBG. Here Mg*ATP and DBG follow a random-equilibrium model, where every substrate can bind indepen-dently from the other substrate. Analyses of the peak integrals from educts and products of the enzymatic reaction, revealed the stoichiometry of the reaction: 1.5 ATP molecules are used to phos-phorylate one DBG molecule. The excess of ATP is attributed to the basal ATPase activity. Further-more, experiments with ATPγS, usually regarded as a non-hydrolysable ATP-analog, where carried out. Surprisingly, DAGK hydrolyzes ATPγS and also transfers the thio-phosphate group to the lipid acceptor DBG, which points to a certain degree of plasticity in the active center. A phosphorylated enzyme intermediate was not detected. These results suggest the building of a ternary complex of Mg*ATP, DBG and DAGK performing a direct-phosphoryl transfer reaction, without passing through a phosphorylated enzyme intermediate. Experiments with the transition state analog ortho-vanadate (Vi) showed a decoupling of the ATP hydrolysis activity from lipid substrate phosphorylation. This indicates a specific transfer site for the γ-phosphate group from ATP to DAG, which can be blocked by Vi.
A general disadvantage of NMR spectroscopy compared to other spectroscopic methods is its inherent low sensitivity. One possible starting point for the improvement of signal-to-noise per unit time is the reduction of the spin-lattice relaxation time of protons. Usually 95 % of the experi-mental time is required for the relaxation of the 1H to equilibrium. The addition of paramagnetic species can be used to reduce the 1H T1. In a comprehensive study four different paramagnetic agents were tested: Cu2+-EDTA, Cu2+-EDTA-tag, Gd3+-TTAHA and Gd3+-DOTA. The titration of these paramagnetic complexes showed the principle feasibility of this approach, but differences between the tested species exist. The most promising complex is Gd3+-DOTA which, at a concentration of 2 mM, causes a 10-time improvement of signal-to-noise ratio per unit time. This allowed measuring 2D 13C-13C correlation spectra of proteoliposomes in one tenth of the usual required experimental time (i.e. 10 hours vs. 4 days) with good signal-to-noise.
For the investigation of structural or dynamic changes in the protein upon substrate interaction with MAS NMR, the spectral properties CP efficiency and resolution of the DAGK in liposomes needed to be improved. The most critical step during sample preparation is the reconstitution of the membrane protein from detergent micelles into a membrane of synthetic lipids under detergent removal. For this procedure the important criteria are enzymatic activity, measured in a coupled ATPase assay, and homogeneity of the proteoliposomes, which was tested e.g. on a discontinuous sucrose step gradient. Therefore an extensive study was carried out, in which different detergents, lipids and lipid mixtures, techniques for detergent removal and different protein-to-lipid ratios were tested. A direct correlation between high ATPase activity and good resolution was not found. Moreover, active DAGK in a mixture of DMPC and cholesterol, which emulates the membrane features of a membrane containing DAG, showed the best CP efficiency and resolution.
The assignment of the protein backbone and amino acid side chains the first mandatory step towards the investigation of structural and dynamical features influencing and defining the enzymatic mechanism by MAS NMR. As the assignment procedure is very time consuming for a total protein, a special labeling scheme for DAGK was developed, which allows assigning most of the protein areas presumably involved in enzyme catalysis. The assignment of DAGK with solution NMR was not transferable to the MAS NMR spectra. Most important for the assignment process were the unique pairs, two consecutive amino acids which only appear once in the amino acid sequence. These unique pairs served as anchor points. Five different multinuclear MAS NMR experiments (DARR, NCO, NCA, NCACX, NCOCX) were required for the sequential assignment. It was possible to assign 35 % of the total amino acid sequence with one sample and 8 experiments acquired at 850 MHz. The secondary structure analysis showed subtle differences to the DAGK assignment with solution NMR, which can be attributed to the different environment in lipid bilayers and detergent micelles.
Data about structural and dynamical changes under substrate interaction can reveal details about the enzymatic mechanism. Therefore changes in chemical shift in 2D heteronuclear correlation experiments in the apo-state and under substrate saturated conditions with the substrates Mg*AMP-PNP, a non-hydrolysable ATP-analog, DOG, a mixture of Mg*AMP-PNP and DOG as well as inhibited by Vi were recorded. The most significant peak changes were observed at the interface membrane-cytoplasm as well as the the N-terminal amphipathic helix. The residues revealing chemical shift perturbations correlate with conserved residues or such residues, for which importance for catalysis and/or folding could be shown in mutation studies. Especially noticeable were the changes at the amino acids Asn 72, Lys 64, His 87, Tyr 86 and Asp 95.
Beside changes of the chemical shift, changes of line width or signal doubling were observable. These changes can point to a correlation with dynamic reorientations in the μs-ms time regime, which are most relevant for enzymatic processes. The protein backbone dynamics in the apo-state as well as saturated with the substrates or inhibited with Vi were investigated with a 15N-CODEX experiment, which is based on the reorientation of the CSA tensor upon dynamical changes. Specific effects of the different substrates or analogs on the protein backbone dynamic were revealed complementing the structural data and the chemical shift perturbation experiments.
Biophysical Studies of Lipid Membranes by Solid State NMR and Molecular Dynamics Simulations
- Biological membranes separate the cell interior from the outside and have diverse functions from signal transduction, apoptosis to transportations of ions and small molecules in and out of the cell. Most of these functions are fulfilled by proteins incorporated in the membrane. However, lipids as the main component of membrane not only serve as structural element for bilayer formation but they are also directly involved e.g. signalling processes and bilayer properties are important to mediate protein interactions. To fully understand the role of lipids, it is necessary to develop a molecular understanding of how certain membrane components modify bulk bilayer structure and dynamics. Membranes are known to have many different motions in different conditions and time scales. Temperature, pH, water content and many other conditions change membrane dynamics in a high degree. In addition to this, time scales of motions in membranes vary from ns to ms range corresponding to fast motion and slow motion, respectively. Therefore, membranes are needed to be studied systematically by varying the conditions and using methods to investigate motions in various time scales separately. The aim of this study was therefore perform a combined solid-state NMR / molecular dynamics study on model membranes. Different substrates, such as potential drugs, polarizing agents and signaling lipids were incorporated into bilayers and their location within the membrane and their effect onto the membrane was probed. NSAIDs (non-steroidal anti-inflammatory drugs), pirinixic acid derivatives, ceramides and polarizing agents were the substrates for membranes in this study. There were several experimental methods that were applied in order to investigate effects of these substrates on membrane dynamics. Different kind of phospholipids including POPC, DMPC and DPPC were used. In addition to experimental work, with the information gathered from solid state NMR experiments molecular dynamics simulations were performed to obtain more information about the membranes at the molecular level. As a result, combination of solid-state NMR with molecular dynamics simulations provides very systematic way of investigating membrane dynamics in a large range of time scales.
Pirinixic acid derivatives were special interest of this study because of their activity on peroxisome proliferator-activated receptor (PPAR) as an agonist as well as on enzymes of microsomal prostaglandin E2 synthase-1 (PGE2s) -1 and 5-lipoxygenase (5-LO) as dual inhibitor. Two potent pirinixic acid derivatives, 2-(4-chloro-6-(quinolin-6-ylamino)pyrimidin-2-ylthio)octanoic acid (compound 2) and 2-(4-chloro-6-(quinolin-6-ylamino)pyrimidin-2-ylthio)octanoate (compound 3), have been worked and their insertion depts were investigated by combining of solid state NMR and molecular dynamics simulations. Both experimental and theoretical results pointed out that compound 3 was inserted the phospholipid bilayer more deeply than 2. NSAIDs – lipid mixtures have been also studied here. It is known that consumption of NSAIDs as in mixture with lipids results much fewer side effects than consumption of the drugs alone. Thus, it is crucial to understand interactions of NSAIDs with lipids and investigate the possible complex formation of drugs with lipids. In this study, interactions of three widely used NSAIDs, ibuprofen, diclofenac and piroxicam, with DPPC were investigated by solid-state NMR. 1H and 31P NMR results depicted that ibuprofen and diclofenac had interactions with lipids, which is an indication of drug-lipid complex formation whereas piroxicam didn’t show any interactions with lipids suggesting that no complex formation occurred in the case of piroxicam. Ceramides are known to play key roles in many cell processes and many studies showed that the functions of ceramides are related with the ceramide effects on biological membranes. Therefore, in this study, influences of ceramides on biophysics of lipid bilayers were investigated by using various solid state NMR techniques and molecular dynamics simulations. Results from molecular dynamics simulations clearly showed that ceramide and lipids have strong interactions. More evidences about ceramide-lipid interactions were provided from 1H and 14N NMR results. In addition, it was indicated by both simulation and experimental methods that ceramide increased the rigidity of DMPC by increasing chain order parameters. BTbk is a biradical, which is used as polarizing agent for dynamic nuclear polarization (DNP) experiments and found to be more efficient than other widely used polarizing agents such as TOTAPOL. Since it is a hydrophobic compound, which prefers to stay inside lipid bilayer it is important to investigate the location and orientation of bTbk along the bilayer in order to understand its enhancement profile in DNP measurements. In this study, both NMR relaxation time measurements and molecular dynamics simulations revealed that bTbk tends to stay more close to hydrophobic chain of lipids than the interfacial part of lipids at bilayer surface.
In the first part of this work, a brief introduction on lipid membranes as well as a theoretical summary on both methods of solid-state NMR and molecular dynamics simulations is given. Then, in the second part methodology is introduced for both solid-state NMR spectrometer and theoretical calculations. Afterwards, results of different membrane systems are discussed in the following parts for both solid state NMR and MD. Finally, in the last part, a summary and the conclusion of the overall results together with some future plans are explained.
Function of plant photosystem II subunits in photoprotection
- Plants absorb sunlight via photosynthetic pigments and convert light energy intochemical energy in the process of photosynthesis. These pigments are mainly bound to antenna protein complexes that funnel the excitation energy to the photosynthetic reaction centres. The peripheral antenna of plant photosystem II (PSII) consists of the major light-harvesting complex of PSII (LHC-II) and the minor LHCs CP29, CP26 and CP24. Light intensity can change frequently and plants need to adapt to high-light conditions in order to avoid photodamage. When more photons are absorbed than can be utilised by the photosynthetic machinery, excessive excitation energy is dissipated as heat by short-term adaptation processes collectively known as non-photochemical quenching (NPQ). A decrease in PSII antenna chlorophyll (Chl) fluorescence yield and a reduction in the average Chl fluorescence lifetime are associated with NPQ. The main component of NPQ is the so-called energy-dependent quenching (qE), and it is triggered by the rapid drop in thylakoid lumenal pH resulting from the plant’s photosynthetic activity. This process is thought to take place at the PSII antenna complexes, which therefore not only capture and transfer light energy but are also involved in balancing the energy flow. The decrease in lumenal pH acivates the enzyme violaxanthin de-epoxidase (VDE), which converts the xanthophyll violaxanthin (Vio) into zeaxanthin (Zea) in the xanthophyll cycle. In addition, the PSII subunit PsbS was discovered to be essential for qE by screening qE-deficient Arabidopsis thaliana mutants. This membrane protein is considered a member of the LHC superfamily, which also includes LHC-II and the minor LHCs. Previous studies on PsbS isolated either from native source or refolded in vitro have produced inconsistent results on its pigment binding capacity. Interestingly, a pH-dependent change in the quaternary structure of PsbS under high light conditions has been reported. This observed dimer-tomonomer transition very likely follows the protonation of lumenal glutamates upon the drop in pH and is accompanied by a change in PSII supercomplex localisation. PsbS dimers are preferentially found in association with the PSII core, whereas PsbS monomers co-localise with LHC-II.Despite the identification of !pH, Zea and PsbS as key players in qE, both the nature of the quencher(s) as well as the underlying molecular mechanism leading to excess energy dissipation still remain unknown. Several models have been put forward to explain the reversible switch in the antenna from an energy-transmitting to a quenched state. Proposals include a simple pigment exchange of Vio for Zea, and aggregation or an internal conformational change of LHC-II. Charge transfer (CT)quenching in the minor LHCs or quenching by carotenoid dark state (Car S1)-Chl interactions have also been suggested. However, none of these qE models has so far been capable of accommodating all the physiological observations and available experimental data. Most importantly, the function of PsbS remains an enigma. A recent qE model suggested that monomerisation of PsbS enables the protein to transiently bind a carotenoid and form a quenching unit with a Chl of a PSII LHC. In view of the various proposed qE mechanisms, this thesis aimed at understanding the interplay of the different qE components and the contribution of the PSII subunits LHC-II, the minor LHCs and PsbS to qE. The initial approach was to investigate the properties of the PSII subunits in the most simple in vitro model system, namely in detergent solution. For this purpose, LHC-II was isolated either from native source or refolded from recombinantly produced protein. Investigation of the minor LHCs and PsbS required heterologous expression and refolding. In addition, experiments were performed on aggregated LHC-II. Aggregates of LHC-II have been used as a popular model system for qE because they exhibit highly quenched Chl fluorescence. At the final stage of this doctoral work, a more sophisticated model system to approximate the thylakoid membrane was developed by reconstitution of the PSII subunits LHC-II and PsbS into liposomes. This system not only allowed for investigation of these membrane proteins in their native environment, but also for mimicking the xanthophyll cycle by distribution of Zea within the membrane as well as !pH by outside buffer exchange. The role of Zea in qE was first investigated with detergent solubilised antenna proteins. The requirement of this xanthophyll for qE is well-known, but the specific contribution to the molecular quenching mechansim is unclear. Previous work had shown that replacement of Vio for Zea in LHC-II was not sufficient to induce Chl fluorescence quenching in Zea-LHC-II, as suggested by the so-called molecular gearshift mechanism. However, by means of selective two-photon excitation spectroscopy, an increase in electronic interactions between Car S1 and Chls was observed for LHC-II upon lowering the pH of the detergent buffer. Electronic Car S1-Chl coupling became even stronger when Zea-LHC-II was probed. The extent of Car S1-Chl coupling correlated directly with the extent of Chl fluorescence quenching, in a similar way as observed previously in live plants under high-light conditions. However, very similar results were obtained with LHC-II aggregates. This implied that the increase in electronic interactions and fluorescence quenching was independent of Zea and low pH. Further experiments on aggregates of LHC-II Chl mutants indicated that the targeted pigments were also not essential for the observed effects. It is proposed that the same molecular mechanism causes an increase in electronic Car S1-Chl interactions and Chl fluorescence quenching in Zea-LHC-II at low pH as well as in aggregated LHC-II. Most likely, surface exposed pigments form random quenching centres in both cases. On the other hand, it was possible that Zea could act as a direct quencher of excess excitation energy in the minor LHCs. However, enrichment of refolded CP29, CP26 and CP24 with Zea did not lead to a change in the Chl excited state lifetime. Formation of a carotenoid radical cation, previously implied in CT quenching, was also not observed, although artificial generation of such a radical cation was principally possible as shown for CP29. During the course of this work, a study reporting the formation of Zea radical cations in minor LHCs was published. Therefore, Zea-enriched minor LHCs were again investigated on the experimental apparatus used in the reported study. Indeed, the presence of at least one carotenoid radical cation for each minor complex was detected. It is suggested that either the preparation method of incubating the refolded minor LHCs with Zea in contrast to refolding the complexes with only Zea and lutein causes the observed differences or that the observed spectral radical cation signatures are due to experimental artifacts. While the experiments with LHC-II and the minor LHCs gave useful insights into the putative qE mechanism, the quencher site and the mode of action of Zea could still not be unambiguously identified. Most importantly, these studies could not explain the function of the qE keyplayer PsbS. Therefore, the focus of the work was shifted to PsbS protein production, purification and characterisation. In view of inconsistent reports on the pigment binding capacity of this PSII subunit, refolding trials with and without photosynthetic pigments were conducted. The formation of a specific pigmentprotein complex typical for other LHCs was not observed and neither was the earlier reported “activation” of Zea for qE by binding to this protein. Nevertheless, PsbS refolded without pigments displayed secondary structure content in agreement with previous studies, indicating pigment-independent folding. Reconstitution of pigmentfree, refolded PsbS into liposomes confirmed that the protein is stable in the absence of pigments. Zea distributed in PsbS-containing liposomes also showed no spectral alteration that would indicate its “activation”. With the ability to reconstitute PsbS, it was then possible to proceed to modelling qE in a proteoliposome system. For this purpose, PsbS was co-reconstituted with LHC-II, which has been reported to interact with PsbS. One-photon excitation (OPE) and two-photon excitation (TPE) spectroscopy measurements were performed on LHC-II- and LHC-II/PsbS-containing liposomes. This enabled both quantification of Chl fluorescence quenching as well as determination of the extent of electronic Car S1-Chl interactions. The effect of Zea was investigated by incorporating it in the proteoliposome membrane. It was shown that Zea alone was not able to induce significant Chl fluorescence quenching when only LHC-II was present. However, when LHC-II and PsbS were co-reconstituted, pronounced Chl fluorescence quenching and an increase in electronic Car S1-Chl interactions were observed and both effects were enhanced when Zea was present. Western blot analysis indicated the presence of a LHC-II/PsbS-heterodimer in these proteoliposomes. In addition to the OPE and TPE measurements, the average Chl fluorescence lifetime was determined in detergent-free buffer at neutral pH and directly after buffer exchange to low pH. No significant changes in the average lifetime were observed for LHC-II proteoliposomes when either Zea was present or after exchange for low pH buffer. This indicated that Zea alone cannot act as a direct quencher, which concurs with the OPE measurements. Moreover, the complex was also properly reconstituted as no aggregation or significant Chl fluorescence quenching were observed. The average lifetime was not significantly affected in LHC-II/PsbS-proteoliposomes, independent of Zea or pH. However, a shortlived component in the presence of a long-lived component was not resolvable with the time resolution of the fluorescence lifetime apparatus.
Implications for qE model systems and the in vivo quenching mechanism are discussed based on the experiments in detergent solution, on LHC-II aggregates and with the proteoliposome model system.
"Wir brauchen den Vergleich mit anderen Forschungseinrichtungen nicht zu scheuen"