André H. Juffer

Proteome-wide Analysis of Lysine Acetylation Suggests its Broad Regulatory Scope in Saccharomyces cerevisiae

Post-translational modification of proteins by lysine acetylation plays important regulatory roles in living cells. The budding yeast Saccharomyces cerevisiae is a widely used unicellular eukaryotic model organism in biomedical research. S. cerevisiae contains several evolutionary conserved lysine acetyltransferases and deacetylases. However, only a few dozen acetylation sites in S. cerevisiae are known, presenting a major obstacle for further understanding the regulatory roles of acetylation in this organism. Here we use high resolution mass spectrometry to identify about 4000 lysine acetylation sites in S. cerevisiae. Acetylated proteins are implicated in the regulation of diverse cytoplasmic and nuclear processes including chromatin organization, mitochondrial metabolism, and protein synthesis. Bioinformatic analysis of yeast acetylation sites shows that acetylated lysines are significantly more conserved compared with nonacetylated lysines. A large fraction of the conserved acetylation sites are present on proteins involved in cellular metabolism, protein synthesis, and protein folding. Furthermore, quantification of the Rpd3-regulated acetylation sites identified several previously known, as well as new putative substrates of this deacetylase. Rpd3 deficiency increased acetylation of the SAGA (Spt-Ada-Gcn5-Acetyltransferase) complex subunit Sgf73 on K33. This acetylation site is located within a critical regulatory domain in Sgf73 that interacts with Ubp8 and is involved in the activation of the Ubp8-containing histone H2B deubiquitylase complex. Our data provides the first global survey of acetylation in budding yeast, and suggests a wide-ranging regulatory scope of this modification. The provided dataset may serve as an important resource for the functional analysis of lysine acetylation in eukaryotes.

Expanding the lysyl hydroxylase toolbox: new insights into the localization and activities of lysyl hydroxylase 3 (LH3).

Hydroxylysine and its glycosylated forms, galactosylhydroxylysine and glucosylgalactosylhydroxylysine, are post‐translational modifications unique to collagenous sequences. They are found in collagens and in many proteins having a collagenous domain in their structure. Since the last published reviews, significant new data have accumulated regarding these modifications. One of the lysyl hydroxylase isoforms, lysyl hydroxylase 3 (LH3), has been shown to possess three catalytic activities required sequentially to produce hydroxylysine and its glycosylated forms, that is, the lysyl hydroxylase (LH), galactosyltransferase (GT), and glucosyltransferase (GGT) activities. Studies on mouse models have revealed the importance of these different activities of LH3 in vivo. LH3 is the main molecule responsible for GGT activity in mouse embryos. A lack of this activity causes intracellular accumulation of type IV collagen, which disrupts the formation of basement membranes (BMs) during mouse embryogenesis and leads to embryonic lethality. The specific inactivation of the LH activity of LH3 causes minor alterations in the structure of the BM and collagen fibril organization, but does not affect the lifespan of mutated mice. Recent data from zebrafish demonstrate that growth cone migration depends critically on the LH3 glycosyltransferase domain. LH3 is located in the ER loosely associated with the membranes, but, unlike the other isoforms, LH3 is also found in the extracellular space in some tissues. LH3 is able to adjust the amount of hydroxylysine and hydroxylysine‐linked carbohydrates of extracellular proteins in their native conformation, suggesting that it may have a role in matrix remodeling. J. Cell. Physiol. 212: 323–329, 2007.

Molecular Dynamics Study of Peptide-Bilayer Adsorption

Two 6-ns simulations of the somatostatin analog sandostatin and a 1-palmityl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer are presented. In the first simulation, the peptide was placed in a region of bulk water density and allowed to spontaneously move toward and bind to the bilayer surface. An attractive force between the peptide and bilayer drove the binding process, which was opposed by a significant frictional force caused by the solvent (water). During the approach of the peptide toward the bilayer the area of the interacting surface between the species was inversely proportional to the distance between them, supporting the application of such a relationship in continuum calculations of peptide-bilayer binding free energies. In the second simulation, the N-terminus of the surface-bound peptide was deprotonated. Consistent with experiment, this strengthened interactions between the peptide and the bilayer. Details of both peptide-bilayer complexes, including the orientation, percent buried surface area, and orientation of the lipid headgroups are in good agreement with those obtained from experiment. The location of the different side chains in the bilayer is in direct correlation with an interfacial hydrophobicity scale developed using model peptides. The aromatic side chains of the Phe and Trp residues all lie flat with respect to the bilayer surface in both complexes. Changes in lipid and water ordering due to peptide binding suggest a possible domination of lipophobic over hydrophobic effects, as proposed by other workers. Where appropriate, peptide and lipid properties in the bound states are compared with separate simulations of sandostatin and the bilayer in water, respectively, so as to monitor the response of the system to the binding process.

Characterization of Collagenous Peptides Bound to Lysyl Hydroxylase Isoforms

Lysyl hydroxylase (LH, EC 1.14.11.4) is the enzyme catalyzing the formation of hydroxylysyl residues in collagens and other proteins with collagenous domains. Although lower species, such as Caenorhabditis elegans, have only one LH orthologue, LH activity in higher species, such as human, rat, and mouse, is present in three molecules, LH1, LH2, and LH3, encoded by three different genes. In addition, LH2 is present in two alternatively spliced forms (LH2a, LH2b). To understand the functions of the four molecular forms of LH in vertebrates, we analyzed differences in the binding and hydroxylation of various collagenous peptides by the LH isoforms. Nine-amino acid-long synthetic peptides on Pepspot were used for the binding analysis and an activity assay to measure hydroxylation. Our data with 727 collagenous peptides indicated that a positive charge on the peptide and specific amino acid residues in close proximity to the lysyl residues in the collagenous sequences are the key factors promoting peptide binding to the LH isoforms. The data suggest that the LH binding site is not a deep hydrophobic pocket but is open and hydrophilic where acidic amino acids play an important role in the binding. The data do not indicate strict sequence specificity for the LH isoforms, but the data indicated that there was a clear preference for some sequences to be bound and hydroxylated by a certain isoform.

Incorporating the effect of ionic strength in free energy calculations using explicit ions.

The incorporation of explicit ions to mimic the effect of ionic strength or to neutralize the overall charge on a system in free energy calculations using molecular dynamics simulations is investigated. The difference in the free energy of hydration between two triosephosphate isomerase inhibitors calculated at five different ion concentrations is used as an example. We show that the free energy difference can be highly sensitive to the presence of explicit ions even in cases where the mutation itself does not involve a change in the overall charge. The effect is most significant if the molecule carries a net charge close to the site mutated. Furthermore, it is shown that the introduction of a small number of ions can lead to very severe sampling problems suggesting that in practical calculations convergence can best be achieved by incorporating either no counterions or by simulating at high ionic strength to ensure sufficient sampling of the ion distribution.

Structural Enzymological Studies of 2-Enoyl Thioester Reductase of the Human Mitochondrial Fas II Pathway: New Insights Into its Substrate Recognition Properties.

Structural and kinetic properties of the human 2-enoyl thioester reductase [mitochondrial enoyl-coenzyme A reductase (MECR)/ETR1] of the mitochondrial fatty acid synthesis (FAS) II pathway have been determined. The crystal structure of this dimeric enzyme (at 2.4 A resolution) suggests that the binding site for the recognition helix of the acyl carrier protein is in a groove between the two adjacent monomers. This groove is connected via the pantetheine binding cleft to the active site. The modeled mode of NADPH binding, using molecular dynamics calculations, suggests that Tyr94 and Trp311 are critical for catalysis, which is supported by enzyme kinetic data. A deep, water-filled pocket, shaped by hydrophobic and polar residues and extending away from the catalytic site, was recognized. This pocket can accommodate a fatty acyl tail of up to 16 carbons. Mutagenesis of the residues near the end of this pocket confirms the importance of this region for the binding of substrate molecules with long fatty acyl tails. Furthermore, the kinetic analysis of the wild-type MECR/ETR1 shows a bimodal distribution of catalytic efficiencies, in agreement with the notion that two major products are generated by the mitochondrial FAS II pathway.

A previously unobserved conformation for the human Pex5p receptor suggests roles for intrinsic flexibility and rigid domain motions in ligand binding

BackgroundThe C-terminal tetratricopeptide (TPR) repeat domain of Pex5p recognises proteins carrying a peroxisomal targeting signal type 1 (PTS1) tripeptide in their C-terminus. Previously, structural data have been obtained from the TPR domain of Pex5p in both the liganded and unliganded states, indicating a conformational change taking place upon cargo protein binding. Such a conformational change would be expected to play a major role both during PTS1 protein recognition as well as in cargo release into the peroxisomal lumen. However, little information is available on the factors that may regulate such structural changes.ResultsWe have used a range of biophysical and computational methods to further analyse the conformational flexibility and ligand binding of Pex5p. A new crystal form for the human Pex5p C-terminal domain (Pex5p(C)) was obtained in the presence of Sr2+ ions, and the structure presents a novel conformation, distinct from all previous liganded and apo crystal structures for Pex5p(C). The difference relates to a near-rigid body movement of two halves of the molecule, and this movement is different from that required to reach a ring-like conformation upon PTS1 ligand binding. The bound Sr2+ ion changes the dynamic properties of Pex5p(C) affecting its conformation, possibly by making the Sr2+-binding loop – located near the hinge region for the observed domain motions – more rigid.ConclusionThe current data indicate that Pex5p(C) is able to sample a range of conformational states in the absence of bound PTS1 ligand. The domain movements between various apo conformations are distinct from those involved in ligand binding, although the differences between all observed conformations so far can be characterised by the movement of the two halves of Pex5p(C) as near-rigid bodies with respect to each other.

Calculation of affinities of peptides for proteins

Several methodologies were employed to calculate the Gibbs standard free energy of binding for a collection of protein–ligand complexes, where the ligand is a peptide and the protein is representative for various protein families. Almost 40 protein–ligand complexes were employed for a continuum approach, which considers the protein and the peptide at the atomic level, but includes solvent as a polarizable continuum. Five protein–ligand complexes were employed for an all‐atom approach that relies on a combination of the double decoupling method with thermodynamic integration and molecular dynamics. These affinities were also computed by means of the linear interaction energy method. Although it generally proved rather difficult to predict the absolute free energies correctly, for some protein families the experimental ranking order was correctly reproduced by the continuum and all‐atom approach. Considerable attention has also been given to correctly analyze the affinities of charged peptides, where it is required to judge the effect of one or more ions that are being decoupled in an all‐atom approach to preserve electroneutrality. The various methods are further judged upon their merits.

Atomic resolution view into the structure–function relationships of the human myelin peripheral membrane protein P2

The structure of the human myelin peripheral membrane protein P2 has been refined at 0.93 Å resolution. In combination with functional experiments in vitro, in vivo and in silico, the fine details of the structure–function relationships in P2 are emerging.

Theoretical investigations of prostatic acid phosphatase

The phosphotyrosyl protein phosphatase activity of prostatic acid phosphatase (PAP) has been well established. It has also been suggested that PAP partly regulates the activity of growth factor receptors by dephosphorylating the autophosphorylysable tyrosines in them. We studied the binding of the peptides from epidermal growth factor receptor (EGFR) and its homolog (ErbB‐2), corresponding to their autophosphorylation sites, to PAP using theoretical modeling and molecular dynamics (MD) simulation methods. Nine different peptides, each with a phosphotyrosine residue, were docked on human PAP. The binding energies of these peptide–PAP complexes were calculated theoretically and compared to experimentally obtained affinities. The peptide AceDNLpYYWDNH2 from ErbB‐2(1197–1203) showed the most favorable free energy of binding when estimated theoretically. The results demonstrate that the presence of another tyrosine residue proximate to C‐terminal of autophosphorylysable Tyr enhances the binding affinity considerably. The presence of a bulky group instead prevents the binding, as is observed in case of peptide AceNLYpYWDQNH2 which failed to bind, both in theoretical calculations and experiments. Thus we demonstarted that PAP could potentially bind to EGFR and Erbb‐2 and dephosphorylate them. Thus it could be involved in the regulation of the function of such receptors. In addition, complexes of a peptide from AngiotensinII and phosphotyrosine(pY) with human PAP were also modeled. The effects of different protonation states of the titratable active site residues on ligand (pY) binding have also been investigated. For a favorable binding His12 and Asp258 should be neutral, His257 should be positively charged and the phosphate group of the ligand should be in PO  43− state. Furthermore, the analysis of protein motion as observed during simulations suggests the loop–loop contact in the PAP dimer to be of importance in cooperativity. Proteins 2005.