Refine
Year of publication
Document Type
- Doctoral Thesis (95)
Has Fulltext
- yes (95)
Is part of the Bibliography
- no (95)
Keywords
- Heterologe Genexpression (3)
- Endothelin (2)
- G-Protein gekoppelte Rezeptoren (2)
- Gentherapie (2)
- HIV (2)
- Membrane Proteins (2)
- NMR-Spektroskopie (2)
- gene therapy (2)
- 5-Lipoxygenase (1)
- APOBEC3G (1)
Institute
- Biochemie und Chemie (63)
- Biochemie, Chemie und Pharmazie (28)
- Pharmazie (5)
- Biowissenschaften (1)
- Georg-Speyer-Haus (1)
The focus of this research was to understand the molecular mechanism that lies behind the insertion of tail-anchored membrane proteins into the ER membrane of yeast cells. State-of-art instruments such as LILBID, and Cryo-EM, combined with the introduction of direct electron detectors, were used to analyze the proteins that capture tail-anchored proteins near the ER membrane and help their releases from a chaperone, an ATPase named Get3. Get3 escorts TA proteins to the ER membrane, where both Get3 and the TA proteins interact sequentially to Get3 membrane bound receptors Get1 and Get2. Get1 and Get2 are homologs of mammalian WRB and CAML.
The native host was used to separately produce Get1, Get2, and the Get2/Get1 single chain constructs. The studies showed that when Get1 is expressed alone, Get1 does not seems to be located in the ER membrane but rather in microbodies like shape organelles (or peroxisome). Interestingly, Get1 seems to be located in the ER membrane when it is linked to Get2 as single chain construct.
The localization study of Get2/Get1 fused to GFP shows from the fluorescence intensity that Get2/Get1.GFP has a tube-like morphology or membrane-enclosed sacs (cisterna), implying that Get2/Get1 is actually targeted to the ER membrane and is likely functional. In other words, Get1 and Get2 stabilize each other in the ER membrane.
The expression of Get2/Get1 was found to be already optimum when expressed as single chain construct because the fluorescence counts did not improve when additives such as DMSO or histidine were added. However, when Get1 and Get2 are expressed separately, additives improve their protein production yield. In 1 liter culture, Get1 yield is increased by about 3 mg and Get2 by 1.8 mg. This can be explained by the space that Get1 and Get2 should occupy within the ER membrane as they must coexist with other membrane components to maintain the homeostasis of the cell. Hence, if there were no gain for single chain construct expression, it meant that Get2/Get1 was already well expressed on its own in ER membrane and has reached its optimum expression without the help of additives. The Get2/Get1 overexpression is more stable, tolerated and less toxic for the cells to express it at a high level.
DDM has proved to be the best detergent from the detergents tested to solubilize Get1, Get2, and Get2/Get1.
Thereafter, Get1, Get2 (data not shown), and Get2/Get1 were successfully purified in DDM micelles.
Furthermore, for the first time using LILBID, the actual study has shown that Get1 and Get2 are predominantly a heterotetramer (2xGet1 and 2xGet2) but higher oligomerization may exist as well.
Get3 binds to Get1 in a biphasic way with a specific strong binding of an affinity of 57 nM and the second of 740 nM nonspecific indicative of heterogeneity within the interaction between Get1 and Get3. This heterogeneity is caused by the presence of different conformation of either protein. However, in order to characterize a high-resolution structure model of a specific target one needs highly homogenous and identical molecules of the target protein or complex in solution. The homogeneity increases the chances of growing crystals during crystallography as the good homogeneity will likely generate a perfect packing of unit cells stack (also known as crystal lattice) in the three-dimensional spaces. The same truth goes for the single particles analysis Cryo-EM, especially for smaller complexes where having less or no conformation alterations of specific targets will enable the researcher to classify the particles in 2D and 3D, therefore improving the signal-to-noise-ratio that will ultimately lead to high-resolution structure determination.
Get1, Get2/Get1 and chimeric variants (tGet2/Get1, T4l.Get2/Get1, T4l.Get2.apocyte.Get1) were crystallized but none of the crystals could diffract due to heterogeneity.
This heterogeneity was not only occurring upon the binding of Get3 to its membrane receptors, but seems to be already present within the receptors themselves through possibly different conformation.
In this Ph.D. thesis, the heterogeneity of purified Get2 and Get1 as complex or individually in detergent is then, so far, the limiting factor for obtaining a high-resolution structure model of Get1 and Get2. As mentioned above, the heterogeneity observed was not due to the quality of the sample preparation but rather to the effect of different conformations that could have been native, or just because of the micelle used, as it was proven by the 3-D heterogeneity classification by Cryo-EM.
In general, crosslinking is one way to keep the integrity of protein complexes, however it appeared not to improve the sample quality when it was analyzed in micelles. Often the integrity of some membrane proteins is affected when they are solubilized and purified in detergents.
Finally, in this study, the structural map of Get2 and Get1 complex linked with chimeric protein T4 lysozyme and apocytochrome C b562RIL gene was obtained at 10 Å. However, this single chain construct has a density map corresponding to heterodimer species (one Get1 and Get2). Therefore, based on those data the tertiary structure of Get2/Get1 in micelle is poorly defined. It could be that the membrane extraction in DDM and the purification destabilizes the structure of the complex.
Die Tumorprotein-Familie des Proteins p53 besteht aus drei Familienmitgliedern p53, p63 und p73 mit diversen Funktionen als Transkriptionsfaktoren. p53 war das erste Mitglied dieser Familie, das im Jahre 1979 entdeckt wurde und wurde zunächst als krebsverursachendes Protein eingeordnet, weil es in vielen Tumorgeweben in erhöhter Menge vorgefunden wurde. Es wurde allerdings festgestellt, dass der Großteil dieser gefundenen p53-Proteine funktionsunfähig durch Mutationen in ihrer Aminosäuresequenz waren. Unmutiertes p53 hingegen führt zu einem Stopp von Zellteilung oder sogar Zelltod, sofern die Zellen genetischem Stress durch Strahlung oder mutagene Chemikalien ausgesetzt sind. Heute wird p53 als eines der wichtigsten Tumor-Unterdrückungsproteine betrachtet. Die beiden anderen Familienmitglieder p63 und p73 existieren in einer Vielzahl von Isoformen. Neben carboxyterminaler alternativer mRNA-Prozessierung (α, β, γ, usw. Isoformen) führen zwei unabhängige Promotoren auch zu zwei unterschiedlichen Aminotermini. Hier wird zwischen ΔN- und TA-Isoformen unterschieden. Im Falle von p63 treten zwei dominante Isoformen auf, ΔNp63α und TAp63α. Während ΔNp63α eine Rolle in der Differenzierung von Haut spielt, wurde TAp63α bisher ausschließlich in Eizellen gefunden. Dort hat es die Funktion eines Sensors, der die genetische Integrität der weiblichen Keimbahn sicherstellt. Es liegt in Eizellen in hoher Konzentration vor, allerdings in einer komplett inaktiven Form. Werden Schäden im der Erbgut der Eizelle festgestellt, so wird das Protein aktiviert und kann so den Prozess des Zelltods der Eizelle einleiten. Mutationen oder das Fehlen des p63-Genes führen zu Missbildungen während der Entwicklung und zu unvollständig ausgebildeter Haut. Im Falle von p73 gibt es ebenfalls mehrere Isoformen, wobei die Funktionen und Relevanzen der einzelnen Isoformen bisher nicht komplett geklärt werden konnten. Eine p73-negative Maus hat einen diffusen Phänotyp, der sich durch niedrige Intelligenz, fast sterile Männchen und chronische bronchiale Infektion auszeichnet. Generell sind alle Mitglieder der p53-Familie tetramere Proteine und sind nur in diesem Zustand auch aktiv. Die einzige Ausnahme stellt, wie oben beschrieben, TAp63α dar, das in einem inaktiven dimeren Zustand vorliegt und nur durch Modifikation durch zwei unabhängige Kinasen aktiviert werden kann. Dabei geht es in den tetrameren Zustand über und ist daraufhin aktiv.
Alle drei Proteine haben (anhand ihrer längsten Isoform beschrieben) eine konservierte Domänenstruktur. Am Aminoterminus befindet sich zunächst die transaktivierende-Domäne (TAD), die für Interaktionen mit transkriptionellen Koaktivatioren relevant ist. Danach folgt die stark konservierte Desoxyribonukleinsäure (DNA) bindende Domäne (DBD). Sie stellt sicher, dass der Transkriptionsfaktor sequenzspezifisch an der richtigen Stelle auf die DNA bindet. Weitergehend folgt die Tetramerisierungsdomäne (TD), welche den oligomeren Zustand des Proteins herstellt. Im Falle von p53 endet das Protein an dieser Stelle, bei p63 und p73 folgen noch das Sterile-Alpha-Motiv (SAM) und die Transkription-inhibierende Domäne (TID). Die SAM Domäne wird generell als Interaktionsdomäne beschrieben, es konnte allerdings bis dato kein Interaktionspartner gefunden werden. Die TID hat einen negativen Einfluss auf die transkriptionelle Aktivität der Proteine. Im Falle von TAp63α interagiert sie zusätzlich mit der TAD um den Dimeren Zustand zu stabilisieren.
Histon Acetylasen
Die Acetylierung von Histonen ist neben deren Methylierung die wichtigste Modifikation. Sie ist essenziell für die Transkription innerhalb aller eukaryontischen Lebewesen, da sie durch die Modifikation von Histonen die DNA für die DNA-Polymerase II zugänglich macht. Es gibt insgesamt fünf verschiedene, nicht näher miteinander verwandte Familien von Histonacetylasen. Diese Studie beschäftigt sich ausschließlich mit der KAT3 Familie, bestehend aus den Proteinen p300 und CBP. Beide sind hochgradig konserviert, in gefalteten Bereichen der Proteine erreicht die Sequenzidentität fast 100%. Beide Proteine scheinen sehr ähnliche Aufgaben zu erfüllen, die jedoch nicht komplett identisch sind. Die Fehlfunktion von einem Allel von CBP führt zum Krankheitsbild des Rubinstein-Taybi-Syndrom (RTS), während ein Mangel an p300 sich in Mäusen auf das Gedächtnis auswirkt. Der komplette Verlust beider Allele eines der Proteine ist immer tödlich, genauso wie auch Verlust jeweils eines Allels bei beiden Proteinen. Insgesamt vier unabhängige Domänen in p300/CBP sind in der Lange die transaktivierende Domänen der p53-Familie zu binden. Bei zwei der Domänen handelt es sich um Zinkfinger-Proteine (Taz1 und Taz2), die anderen beiden sind kleine, ausschließlich α-helikale Domänen (Kix und IBiD).
Diese Studie beschäftigt sich mit der Lösung von Strukturen von der transaktivierenden Domäne von p63 und p73 mit der p300-Domäne Taz2. Außerdem wurden die Auswirkungen von direkten Acetylierungen von TAp63α charakterisiert und der Effekt von einem potenten p300/CBP Inhibitor auf Oozyten unter genotoxischem Stress analysiert. Zusätzlich wurde die Phosphorylierungskinetiken von Tap63α wärend der Aktivierung durch Kinasen untersucht.
...
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)[6] as port of the membrane-derived oligosaccharide cycle[31,34]. Surprisingly, DAGK does not share structural or sequential similarities with other kinases[12]. Typical sequence motives found in other kinases, which catalyze phosphoryl transfer reactions, are not found[13]. 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[56]. 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[209]. 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[233]. 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[55], 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[132] was not transferable to the MAS NMR spectra. Most important for the assignment process were the unique pairs[335], 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[132], 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[8]. 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[350]. 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.
Integral membrane proteins (IMPs) account for 20-40% of all open reading frames in fully sequenced genomes and they are target of approximately 60% of all modern drugs. So far, cellular expression systems are often very insufficient for the high-level production of IMPs. Toxic effects, instability or formation of inclusion bodies are frequently observed effects that prevent the synthesis of sufficient amounts of functional protein. I have successfully established an individual cell-free (CF) expression system to overcome these IMP synthesis difficulties. The CF system was established in two different expression modes. If no hydrophobic compartment is provided, the IMPs precipitate in the reaction mixture. Interestingly, these insoluble proteins are found to differ from inclusion bodies as they readily solubilize in mild detergents and the bacterial small multi drug transporter EmrE, expressed in the insoluble mode was shown to reconstitute into liposomes in an active form. Alternatively, IMPs can be synthesized in a soluble way by supplementing the CF system with detergents. A comprehensive overview of 24 commonly used detergents was provided by analyzing their impact on the CF system as well as their ability to keep three structurally very different proteins in solution. The class of long chain polyoxyethylene-alkyl-ethers turned out to be most suitable for soluble expression of a-helical EmrE, the bacterial b-barrel type nucleoside transporter Tsx and the porcine vasopressin receptor type 2, resulting in several mg of protein per mL of reaction mixture. So far IMPs have almost completely been excluded from solution nuclear magnetic resonance (NMR) analyses. I could demonstrate that CF expression enables efficient isotopic labeling of IMPs for NMR analysis and further facilitates selective labeling strategies with combinations of 13C and 15N enriched amino acids that have not been feasible before. Four different G-protein coupled receptors (GPCRs) were successfully CF expressed in preparative scale and for the human endothelin B receptor (ETB), ligand binding ability was observed. A series of truncated ETB derivatives containing nested terminal deletions have been CF produced and functionally characterized. The core area essential for Endothelin-1 binding as well as a central region responsible for ETB oligomer formation was confined to a 39 amino acid fragment including the proposed transmembrane segment 1. The binding constant (KD) of ETB was determined to 6 nM for circular ET-1 by SPR and 29 nM for linear ET-1 by TIRFS. This data indicate a large potential of the established individual CF expression system for functional IMP synthesis.
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.