Refine
Document Type
- Doctoral Thesis (6) (remove)
Language
- English (6) (remove)
Has Fulltext
- yes (6)
Is part of the Bibliography
- no (6) (remove)
Keywords
- Proteine (6) (remove)
Institute
In the present work, the problem of protein folding is addressed from the point of view of equilibrium thermodynamics. The conformation of a globular protein in solution at common temperatures is quite complicated without any geometrical symmetry, but it is an ordered state in the sense of its biological activity. This complicated conformation of a single protein molecule is destroyed upon increasing the temperature or by the addition of appropriate chemical agents, as is revealed by the loss of its activity and change of the physical properties, and so on. Once the complicated native structures having biological activity are lost, it would be natural to suppose that the native structure could hardly be restored. Nevertheless, pioneers, such as Anson and Mirsky, recognized as early as in 1925 that this was not always the case. If one defines the folded and unfolded states of a protein as two distinct phases of a system, then under the variation of temperature the system is transformed from one phase state into another and vice versa. The process of protein folding is accompanied by the release or absorption of a certain amount of energy, corresponding to the first-oder-type phase transitions in the bulk. Knowing the partition function of the system one can evaluate its energy and heat capacity under different temperatures. This task was performed in this work. The results of the developed statistical mechanics model were compared with the results of molecular dynamic simulations of alanine poylpeptides. In particular, the dependencies on temperature of the total energy of the system and heat capacity were compared for alanine polypeptides consisting of 21, 30, 40, 50 and 100 amino acids. The good correspondence of the results of the theoretical model with the results of molecular dynamics simulations allowed to validate the assumptions made about the system and to establish the accuracy range of the theory. In order to perform the comparison of the results of theoretical model and the molecular dynamics simulations it is necessary to perform the efficient analysis of the results of molecular dynamics simulations. This task was also addressed in the present work. In particular, different ways to obtain dependence of the heat capacity on temperature from molecular dynamics simulations are discussed and the most efficient one is proposed. The present thesis reports the result of molecular dynamic simulations for not only alanine polypeptides by also for valine and leucine polypeptides. In valine and leucine polypeptides, it is also possible to observe the helix↔random coil transitions with the increase of temperature. The current thesis presents a work that starts with the investigation of the fundamental degrees of freedom in polypeptides that are responsible for the conformational transitions. Then this knowledge is applied for the statistical mechanics description of helix↔coil transitions in polypeptides. Finally, the theoretical formalism is generalized for the case of proteins in water environment and the comparison of the results of the statistical mechanics model with the experimental measurements of the heat capacity on temperature dependencies for two globular proteins is performed. The presented formalism is based on fundamental physical properties of the system and provides the possibility to describe the folding↔unfolding transitions quantitatively. The combination of these two facts is the major novelty of the presented approach in comparison to the existing ones. The “transparent” physical nature of the formalism provides a possibility to further apply it to a large variety of systems and processes. For instance, it can be used for investigation of the influence of the mutations in the proteins on their stability. This task is of primary importance for design of novel proteins and drug delivering molecules in medicine. It can provide further insights into the problem of protein aggregation and formation of amyloids. The problem of protein aggregation is closely associated with various illnesses such as Alzheimer and mad cow disease. With certain modifications, the presented theoretical method can be applied to the description of the protein crystallization process, which is important for the determination of the structure of proteins with X-Rays. There many other possible applications of the ideas described in the thesis. For instance, the similar formalism can be developed for the description of melting and unzipping of DNA, growth of nanotubes, formation of fullerenes, etc.
Transport of proteins into or across cellular membranes is mediated by the conserved and ubiquitous Sec-machinery. The Sec-homologue in the inner membrane of Escherichia coli is SecYEG. Sec-mediated insertion of numerous membrane proteins is aided by YidC, another protein integral to the inner membrane of Escherichia coli. YidC fulfils in addition the integration of a variety of membrane proteins Sec-independently. It belongs to a conserved but structurally uncharacterised family of proteins important for membrane protein biogenesis and comprises homologues in mitochondria and chloroplasts. By modification of a former crystallisation protocol two-dimensional crystals of SecYEG were grown in presence of the signal sequence peptide of LamB. Recording of structural data by electron cryo-microscopy and calculation of a difference structure comparing a former SecYEG projection structure with the one of SecYEG crystallised in presence of the substrate revealed several new and vacant densities. These hint to signal peptide binding close to the translocation pore and to significant rearrangements in proximity to the lateral exit site for transmembrane domains in SecYEG. The difference structure suggests that dimeric SecYEG is an asymmetric molecule consisting of one active and one inactive SecYEG monomer. Detergent removal from a mixture of purified YidC and lipids produced two-dimensional crystals that were highly dependent on the ionic strength and lipid composition for their growth. Electron cryo-microscopy on the frozen-hydrated crystals and image processing visualised structural details at about 10 Å resolution. Averaging two alternative projection structures in p2 and p121_a symmetry, respectively, yielded essentially the same features. Four YidC monomers form one unit cell (dimensions 82 x 71 Å, included angle 85 ° and 90 °, respectively) and seem to be arranged as two sets of dimers integrated in an anti-parallel fashion into the membrane. An area of low density in the centre of each YidC monomer resembles possibly a constriction of the membrane, which could have particular relevance for the integration of substrate proteins into the lipid bilayer.
Results were presented from Brownian dynamics simulations for cyt c molecules approximated as spherical particles with diameter 2R ' 3.3 nm interacting with a charged planar membrane surface. Using the well-known Ermak-McCammon algorithm of ref. [36, 37] for solving the Langevin equations (see Chapter 2), a new computer program in C++ was developed. An overview of the way it is implemented is given in Chapter 3. The program in its current state is able to compute the trajectories (translation and rotation) of hundreds of spherical particles in systems with typical dimensions of 103 − 1003 nm3 . As explained in the introductory Chapter 1 the motivation for studying the dynamics of cyt c molecules in such systems came from the progress in the research of photosynthetic bacteria, e.g. While the internal processes of energy transduction (light harvesting, channelling to RC, charge separation) are quite well understood, the dynamics of soluble cyt c as an electron transporter in this context is not yet clear. In many textbooks one can find illustrations where a single cyt c is responsible for the electron transport between two integral membrane proteins (the reaction centre RC and the bc1 complex). But as pointed out in publications like refs. [49], [59], [60], [61] or [62] biological cells are crowded with different molecules. Consequently, one can assume that the electron transport between two integral membrane proteins is not simply taken on by one single cyt c molecule. Instead it is likely that many of these particles are located in a cyt c pool above the membrane and that they perform the electron transport in turns. Thus, it is desirable to have a simulation package that is able to compute the trajectories of many proteins. Note that the detailed processes of electron transfer and binding to membrane proteins are not modelled here. The details of these processes are quite complicated so that we refrained from including them in the coarse-grained simulations. Here, the actual binding is simply defined by a particle distance zb from the membrane which marks the beginning of the attractive potential. ...
Ubiquitin and the ubiquitin-like protein ATG8 are covalently attached to their respective targets via a coordinated cascade involving E1 activating, E2 conjugating and E3 ligating enzymes. Whereas ubiquitin is conferred to proteins as mono- and/or polymer(s) to alter their stability, localization and/or activity, the ubiquitin-like modifier (UBL) ATG8 is conjugated to the phospholipid phosphatidylethanolamine (PE). The best understood function of ATG8 is during autophagy where ATG8-PE conjugates are incorporated into both layers of incipient autophagosomes and serve as multipurpose docking sites for autophagosomal cargo receptors as well as regulatory factors (termed adaptors) that drive formation and maturation of autophagosomes. Mammalian cells harbor six ATG8 family members that can be subclassified into the LC3- and GABARAP-family and that can all be lipidated. However, it is currently unclear to what extent these proteins are functionally redundant or fulfil unique roles.
Cullin-RING ligase complexes (CRLs) are modular E3 ubiquitin ligases that comprise a RING-finger protein that associates with the ubiquitin-charged E2 enzyme, a substrate recruiting module as well as a cullin scaffold as a linker between RING protein and substrate adaptor. Whereas SCF (SKP1-CUL1-F-box protein) complexes, the most studied CRLs, harbor cullin-1 (CUL1) as scaffold and F-box proteins as substrate binding modules, CUL3-containing CRL complexes employ cullin-3 (CUL3), RING-box protein 1 (RBX1) and BTB proteins as substrate adaptors. Here, the BTB domain serves as binding interface for CUL3 and is usually complemented by an additional protein-protein interaction domain such as MATH or Kelch that mediates binding to the substrate of the E3 ligase complex.
Besides ubiquitylation, guanine nucleotide binding is another common way to regulate protein activity and signaling in cells. Here, small Rho GTPases cycle between active and inactive states by binding of the guanine nucleotides GTP or GDP with the help of regulatory proteins. Whereas GTPase-activating proteins (GAP) render RAC1 inactive by facilitating GTP hydrolysis, guanine exchange factors (GEF) such as T-lymphoma invasion and metastasis-inducing protein 1 (TIAM1) activate RAC1 by stimulating the exchange of GDP to GTP. Local control of RAC1 activity is essential to allow a specific cellular response to stimuli such as growth factors or migratory impulses.
This study reports an unexpected link between the GABARAP subfamily of mammalian ATG8 proteins, the ubiquitin proteasome system and RAC1 through the ubiquitylation of the RAC1 GEF TIAM1. The Kelch repeat and BTB domain-containing proteins 6 (KBTBD6) and 7 (KBTBD7) were established as heterodimeric substrate adaptors for CUL3. Interestingly, a thorough proteomic analysis revealed a number of putative substrates but, out of 11 substrate candidates tested, only the RAC1 GEF TIAM1 appeared to be influenced by depletion of CUL3KBTBD6/KBTBD7. Binding studies showed that KBTBD7 binds TIAM1 via the Kelch repeats and that this binding was markedly enhanced when CUL3 activation was abolished upon treatment with the neddylation inhibitor MLN4924. Also, total TIAM1 abundance was increased upon CUL3KBTBD6/KBTBD7 depletion and accumulation of TIAM1 upon proteasome inhibition suggested that TIAM1 is degraded via the proteasome. In vivo ubiquitylation assays and denaturing immunoprecipitations as well as mass spectrometrical analysis confirmed that CUL3KBTBD6/KBTBD7 ubiquitylates TIAM1 at two distinct lysines (K1404 and K1420) close to its C-terminus.
Previously, KBTBD6 and KBTBD7 were found as interactors of several members of the human ATG8 family of proteins in a proteomic study analyzing the human autophagy network. This association was confirmed in the present work. Furthermore, peptide array technology and mutational analysis revealed that KBTBD6 and KBTBD7 employ a classical ATG8-family interacting motif (AIM; also referred to as LC3-interacting region or LIR) as binding interface. The AIMs of KBTBD6 (W-V-R-V) and KBTBD7 (W-V-Q-V) fulfil the consensus AIM sequence motif (F/W/Y1-X2-X3-I/L/V4) and are preceded by several acidic residues and serines. A series of structural and cell biological experiments revealed a binding preference for the GABARAP subfamily of human ATG8 proteins and most importantly, a requirement of the GABARAP-KBTBD6 and -KBTBD7 interaction for TIAM1 ubiquitylation. The finding that TIAM1 binding to KBTBD6 and KBTBD7 AIM mutants was diminished raised the possibility that GABARAP binding mediates the recruitment of CUL3KBTBD6/KBTBD7 to membranes where TIAM1 is localized. Interestingly, colocalization of KBTBD6, GABARAPL1 and TIAM1 in punctuate structures could be observed. Since only a very small fraction of GABARAPL1 colocalized with LC3B, and colocalization between KBTBD6 and LC3B was not observed, these vesicular structures are most likely distinct from autophagosomes. Furthermore, TIAM1 ubiquitylation was reduced when GABARAP, but not LC3B, was depleted or when lipidation of GABARAP was prevented.
Stabilization of TIAM1 upon KBTBD6 and/or KBTBD7 depletion led to elevated TIAM1-dependent RAC1 activity, altered actin morphology with increased cortical actin and loss of vinculin foci. Re-introduction of wild-type KBTBD6 or KBTBD7 but not AIM mutants reverted all these phenotypes. Moreover, depletion of KBTBD6 or KBTBD7 in human breast cancer cells massively increased their invasiveness, whereas TIAM1 knockdown had the opposite outcome. All physiological effects of KBTBD6 and KBTBD7 depletion were inhibited by additional depletion of TIAM1 or RAC1 confirming that the phenotypes observed are indeed mediated by the CUL3KBTBD6/KBTBD7-TIAM1-RAC1 signaling pathway. Intriguingly, KBTBD6 and KBTBD7 were not subject to autophagosomal degradation, thereby establishing a new function for GABARAP proteins beyond autophagosomal degradation in providing a signaling platform for recruitment of the E3 ligase CUL3KBTBD6/KBTBD7 in close proximity to its substrate TIAM1, enabling localized ubiquitylation.
Local restricted control of RAC1 activity by ubiquitylation has been described for TIAM1-RAC1 signaling previously. Examples are HECT, UBA and WWE domain-containing protein 1 (HUWE1)-mediated TIAM1 ubiquitylation that occurs predominantly at cell-cell-junctions in response to hepatocyte growth factor stimulation in MDCKII cells or inhibition of RAC1 activity by the RAC1 GAP protein BCR (breakpoint cluster region) at the leading edge of astrocytes through binding to the TIAM1-Par (polarity) complex. SCFBTRC mediates ubiquitylation of TIAM1 in response to mitogens or DNA damage, though it has not been explored whether this regulation is spatially restricted. Thus, this study adds a novel layer of complexity to the spatial regulation of RAC1 signaling by implicating membrane-bound human ATG8 proteins in this process.
Also, this study is the first report specifically implicating the GABARAP proteins in cellular signaling events. It will be interesting to explore whether the concept of localized signaling mediated by GABARAPs applies to other substrates of CUL3KBTBD6/KBTBD7 and membranerelated signaling processes in which GABARAP proteins are involved. Controlling RAC1 activity at GABARAP-decorated membranes might also be important for trafficking events or autophagy since it was described that RAC1 has an inhibitory function on autophagy. Therefore, spatial restricted ubiquitylation of TIAM1 resulting in specific deactivation of RAC1 could promote the autophagic process when locally needed. Although the catalytic mTOR inhibitor Torin1 and the lysosomal H+ ATPase inhibitor BafilomycinA1 promoted TIAM1 ubiquitylation by increasing the pool of membrane-conjugated GABARAP, but other signals that stimulate GABARAP-KBTBD6/KBTBD7 association and subsequent TIAM1 ubiquitylation are to be identified. Besides, determining the KBTBD6/KBTBD7 binding site in TIAM1 or uncovering a deubiquitylating enzyme (DUB) that locally counteracts the ubiquitylation of TIAM1 will enable a better comprehension of the complete localized signaling cascade.
The heat stress response is characterized by the presence of heat stress transcription factors (Hsfs) which mediate transcription of heat stress genes. In tomato (Lycopersicon peruvianum) cell cultures the simultaneous expression of four Hsfs, which are either constitutively (HsfA1 and HsfA3) or heat-stress inducible (HsfA2 and HsfB1) expressed, results in a complex network with dynamically changing cellular levels, intracellular localization and functional interactions. In order to examine the relevance of their multiplicity as well as to get more insights into the complexity of the plant heat stress response, the individual tomato Hsfs were investigated with respect to their protein interactions in vitro and in vivo. To this aim, I used pull-down assays as well as yeast assays to study the following aspects: 1. Oligomeric state of Hsfs: the results show that all class A Hsfs (HsfA1, HsfA2 and HsfA3) are trimeric proteins and interact with each other via the oligomerization (HR-A/B) domain. The similarity of their HRA/B regions allows formation of homo- and heterooligomeric complexes between all class A Hsfs. This special property was investigated by mutational studies with HsfA2 indicating that the linker and the HR-B regions are the minimal part required for Hsf/Hsf interactions. The conserved hydrophobic amino acid residues of the HR-B region are most important whereas the amino acid residues of the linker may provide higher flexibility to the HR-B region. Another investigated factor was HsfB1. HsfB1 is a member of class B Hsfs, which are characterized by an oligomerization domain without the 21 amino acid residues linker inserted between the HR-A and HR-B regions. It has a low activator potential and exists exclusively as dimer. HsfB1 can not physically interact with class A Hsfs. However, HsfB1 and HsfA1, binding to adjacent HSE sites, are assumed to cause strong synergistic effects in gene activation. 2. Potential HsfB1 interacting proteins: we searched for HsfB1 interacting proteins by using recombinant His-tagged proteins with HsfB1 as baits in pull-down assays. Histones H2A, H2B and H4 were identified by means of Peptide Mass Finger Printing and N-terminal sequencing analyses. The three histones represent the major proteins in tomato whole cell extracts retrieved by HsfB1. 3. HsfA2/small heat stress proteins (sHsps) interaction: pull-down and yeast two-hybrid assays were used to study the specific interaction of HsfA2 with tomato class II sHsp. This interaction occurs via the oligomerization domain of HsfA2. Other members of the plant Hsp20 family, including class I sHsp, do not interact with HsfA2. Heterooligomers of HsfA2 with class II sHsp may represent precursor forms of the plant higher molecular weight cytoplasmic complexes of heat stress granules, which form during heat stress. The findings presented in this thesis are a contribution to support the concept of a Hsfs network via protein-protein interactions. These data, together with information obtained from other studies, are used to propose a tentative model of the complex Hsfs network controlling the plant heat stress response.
Specific functions of biological systems often require conformational transitions of macromolecules. Thus, being able to describe and predict conformational changes of biological macromolecules is not only important for understanding their impact on biological function, but will also have implications for the modelling of (macro)molecular complex formation and in structure-based drug design approaches. The “conformational selection model” provides the foundation for computational investigations of conformational fluctuations of the unbound protein state. These fluctuations may reveal conformational states adopted by the bound proteins. The aim of this work is to incorporate directional information in a geometry-based approach, in order to sample biologically relevant conformational space extensively. Interestingly, coarse-grained normal mode (CGNM) approaches, e.g., the elastic network model (ENM) and rigid cluster normal mode analysis (RCNMA), have emerged recently and provide directions of intrinsic motions in terms of harmonic modes (also called normal modes). In my previous work and in other studies it has been shown that conformational changes upon ligand binding occur along a few low-energy modes of unbound proteins and can be efficiently calculated by CGNM approaches. In order to explore the validity and the applicability of CGNM approaches, a large-scale comparison of essential dynamics (ED) modes from molecular dynamics (MD) simulations and normal modes from CGNM was performed over a dataset of 335 proteins. Despite high coarse-graining, low frequency normal modes from CGNM correlate very well with ED modes in terms of directions of motions (average maximal overlap is 0.65) and relative amplitudes of motions (average maximal overlap is 0.73). In order to exploit the potential of CGNM approaches, I have developed a three-step approach for efficient exploration of intrinsic motions of proteins. The first two steps are based on recent developments in rigidity and elastic network theory. Initially, static properties of the protein are determined by decomposing the protein into rigid clusters using the graph-theoretical approach FIRST at an all-atom representation of the protein. In a second step, dynamic properties of the molecule are revealed by the rotations-translations of blocks approach (RTB) using an elastic network model representation of the coarse-grained protein. In the final step, the recently introduced idea of constrained geometric simulations of diffusive motions in proteins is extended for efficient sampling of conformational space. Here, the low-energy (frequency) normal modes provided by the RCNMA approach are used to guide the backbone motions. The NMSim approach was validated on hen egg white lysozyme by comparing it to previously mentioned simulation methods in terms of residue fluctuations, conformational space explorations, essential dynamics, sampling of side-chain rotamers, and structural quality. Residue fluctuations in NMSim generated ensemble is found to be in good agreement with MD fluctuations with a correlation coefficient of around 0.79. A comparison of different geometry-based simulation approaches shows that FRODA is restricted in sampling the backbone conformational space. CONCOORD is restricted in sampling the side-chain conformational space. NMSim sufficiently samples both the backbone and the side-chain conformations taking experimental structures and conformations from the state of the art MD simulation as reference. The NMSim approach is also applied to a dataset of proteins where conformational changes have been observed experimentally, either in domain or functionally important loop regions. The NMSim simulations starting from the unbound structures are able to reach conformations similar to ligand bound conformations (RMSD < 2.4 Å) in 4 out of 5 cases of domain moving proteins. In these four cases, good correlation coefficients (R > 0.7) between the RMS fluctuations derived from NMSim generated structures and two experimental structures are observed. Furthermore, intrinsic fluctuations in NMSim simulation correlate with the region of loop conformational changes observed upon ligand binding in 2 out of 3 cases. The NMSim generated pathway of conformational change from the unbound structure to the ligand bound structure of adenylate kinase is validated by a comparison to experimental structures reflecting different states of the pathway as proposed by previous studies. Interestingly, the generated pathway confirms that the LID domain closure precedes the closing of the NMPbind domain, even if no target conformation is provided in NMSim. Hence, the results in this study show that, incorporating directional information in the geometry-based approach NMSim improves the sampling of biologically relevant conformational space and provides a computationally efficient alternative to state of the art MD simulations.