Metals in proteins from sulphate-reducing bacteria: adenylate kinase and ATP sulfurylase. Proteins containing cobalt, zinc, iron (II) ions

A Thesis submitted at the Faculty Science and Technology of the New University of Lisbon for a degree in Doctor of Philosophy in Biochemistry with specialization in Physical Biochemistry

Conservation of complex knotting and slipknotting patterns in proteins.

While analyzing all available protein structures for the presence of knots and slipknots, we detected a strict conservation of complex knotting patterns within and between several protein families despite their large sequence divergence. Because protein folding pathways leading to knotted native protein structures are slower and less efficient than those leading to unknotted proteins with similar size and sequence, the strict conservation of the knotting patterns indicates an important physiological role of knots and slipknots in these proteins. Although little is known about the functional role of knots, recent studies have demonstrated a protein-stabilizing ability of knots and slipknots. Some of the conserved knotting patterns occur in proteins forming transmembrane channels where the slipknot loop seems to strap together the transmembrane helices forming the channel.

Identifying knots in proteins.

Polypeptide chains form open knots in many proteins. How these knotted proteins fold and finding the evolutionary advantage provided by these knots are among some of the key questions currently being studied in the protein folding field. The detection and identification of protein knots are substantial challenges. Different methods and many variations of them have been employed, but they can give different results for the same protein. In the present article, we review the various knot identification algorithms and compare their relative strengths when applied to the study of knots in proteins. We show that the statistical approach based on the uniform closure method is advantageous in comparison with other methods used to characterize protein knots.

Emergence of novel domains in proteins

Proteins are composed of a combination of discrete, well-defined, sequence domains, associated with specific functions that have arisen at different times during evolutionary history. The emergence of novel domains is related to protein functional diversification and adaptation. But currently little is known about how novel domains arise and how they subsequently evolve. To gain insights into the impact of recently emerged domains in protein evolution we have identified all human young protein domains that have emerged in approximately the past 550 million years. We have classified them into vertebrate-specific and mammalian-specific groups, and compared them to older domains. We have found 426 different annotated young domains, totalling 995 domain occurrences, which represent about 12.3% of all human domains. We have observed that 61.3% of them arose in newly formed genes, while the remaining 38.7% are found combined with older domains, and have very likely emerged in the context of a previously existing protein. Young domains are preferentially located at the N-terminus of the protein, indicating that, at least in vertebrates, novel functional sequences often emerge there. Furthermore, young domains show significantly higher non-synonymous to synonymous substitution rates than older domains using human and mouse orthologous sequence comparisons. This is also true when we compare young and old domains located in the same protein, suggesting that recently arisen domains tend to evolve in a less constrained manner than older domains. We conclude that proteins tend to gain domains over time, becoming progressively longer. We show that many proteins are made of domains of different age, and that the fastest evolving parts correspond to the domains that have been acquired more recently.

Knot localization in proteins.

The backbones of proteins form linear chains. In the case of some proteins, these chains can be characterized as forming linear open knots. The knot type of the full chain reveals only limited information about the entanglement of the chain since, for example, subchains of an unknotted protein can form knots and subchains of a knotted protein can form different types of knots than the entire protein. To understand fully the entanglement within the backbone of a given protein, a complete analysis of the knotting within all of the subchains of that protein is necessary. In the present article, we review efforts to characterize the full knotting complexity within individual proteins and present a matrix that conveys information about various aspects of protein knotting. For a given protein, this matrix identifies the precise localization of knotted regions and shows the knot types formed by all subchains. The pattern in the matrix can be considered as a knotting fingerprint of that protein. We observe that knotting fingerprints of distantly related knotted proteins are strongly conserved during evolution and discuss how some characteristic motifs in the knotting fingerprints are related to the structure of the knotted regions and their possible biological role.

Statistical Methods for Conservation and Alignment Quality in Proteins

Construction of multiple sequence alignments is a fundamental task in Bioinformatics. Multiple sequence alignments are used as a prerequisite in many Bioinformatics methods, and subsequently the quality of such methods can be critically dependent on the quality of the alignment. However, automatic construction of a multiple sequence alignment for a set of remotely related sequences does not always provide biologically relevant alignments.Therefore, there is a need for an objective approach for evaluating the quality of automatically aligned sequences. The profile hidden Markov model is a powerful approach in comparative genomics. In the profile hidden Markov model, the symbol probabilities are estimated at each conserved alignment position. This can increase the dimension of parameter space and cause an overfitting problem. These two research problems are both related to conservation. We have developed statistical measures for quantifying the conservation of multiple sequence alignments. Two types of methods are considered, those identifying conserved residues in an alignment position, and those calculating positional conservation scores. The positional conservation score was exploited in a statistical prediction model for assessing the quality of multiple sequence alignments. The residue conservation score was used as part of the emission probability estimation method proposed for profile hidden Markov models. The results of the predicted alignment quality score highly correlated with the correct alignment quality scores, indicating that our method is reliable for assessing the quality of any multiple sequence alignment. The comparison of the emission probability estimation method with the maximum likelihood method showed that the number of estimated parameters in the model was dramatically decreased, while the same level of accuracy was maintained. To conclude, we have shown that conservation can be successfully used in the statistical model for alignment quality assessment and in the estimation of emission probabilities in the profile hidden Markov models.

Alterations in proteins of bone marrow extracellular matrix in undernourished mice

The objective of the present study was to determine the effect of protein malnutrition on the glycoprotein content of bone marrow extracellular matrix (ECM). Two-month-old male Swiss mice were submitted to protein malnutrition with a low-protein diet containing 4% casein as compared to 20% casein in the control diet. When the experimental group had attained a 20% loss of their original body weight, we extracted the ECM proteins from bone marrow with PBS buffer, and analyzed ECM samples by SDS-PAGE (7.5%) and ECL Western blotting. Quantitative differences were observed between control and experimental groups. Bone marrow ECM from undernourished mice had greater amounts of extractable fibronectin (1.6-fold increase) and laminin (4.8-fold increase) when compared to the control group. These results suggest an association between fluctuations in the composition of the hematopoietic microenvironment and altered hematopoiesis observed in undernourished mice.

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favourability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB) versus Thr(HB)-Asn(nHB). Significantly (p <= 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi(1) rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi(1) conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi(1) and chi(2) conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favourability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http:// www.rubic.rdg.ac.uk/betapairprefsparallel/). (c) 2005 Elsevier Ltd. All rights reserved.

Analysis of void volumes in proteins and application to stability of the p53 tumour suppressor protein

We have developed a new method for the analysis of voids in proteins (defined as empty cavities not accessible to solvent). This method combines analysis of individual discrete voids with analysis of packing quality. While these are different aspects of the same effect, they have traditionally been analysed using different approaches. The method has been applied to the calculation of total void volume and maximum void size in a non-redundant set of protein domains and has been used to examine correlations between thermal stability and void size. The tumour-suppressor protein p53 has then been compared with the non-redundant data set to determine whether its low thermal stability results from poor packing. We found that p53 has average packing, but the detrimental effects of some previously unexplained mutations to p53 observed in cancer can be explained by the creation of unusually large voids. (C) 2004 Elsevier Ltd. All rights reserved.

Prediction and functional analysis of native disorder in proteins from the Three Kingdoms of Life

An automatic method for recognizing natively disordered regions from amino acid sequence is described and benchmarked against predictors that were assessed at the latest critical assessment of techniques for protein structure prediction (CASP) experiment. The method attains a Wilcoxon score of 90.0, which represents a statistically significant improvement on the methods evaluated on the same targets at CASP. The classifier, DISOPRED2, was used to estimate the frequency of native disorder in several representative genomes from the three kingdoms of life. Putative, long (>30 residue) disordered segments are found to occur in 2.0% of archaean, 4.2% of eubacterial and 33.0% of eukaryotic proteins. The function of proteins with long predicted regions of disorder was investigated using the gene ontology annotations supplied with the Saccharomyces genome database. The analysis of the yeast proteome suggests that proteins containing disorder are often located in the cell nucleus and are involved in the regulation of transcription and cell signalling. The results also indicate that native disorder is associated with the molecular functions of kinase activity and nucleic acid binding.

OX133, a monoclonal antibody recognizing protein-bound N-ethylmaleimide for the identification of reduced disulfide bonds in proteins

In vivo, enzymatic reduction of some protein disulfide bonds, allosteric disulfide bonds, provides an important level of structural and functional regulation. The free cysteine residues generated can be labeled by maleimide reagents, including biotin derivatives, allowing the reduced protein to be detected or purified. During the screening of monoclonal antibodies for those specific for the reduced forms of proteins, we isolated OX133, a unique antibody that recognizes polypeptide resident, N-ethylmaleimide (NEM)-modified cysteine residues in a sequence-independent manner. OX133 offers an alternative to biotin-maleimide reagents for labeling reduced/alkylated antigens and capturing reduced/alkylated proteins with the advantage that NEM-modified proteins are more easily detected in mass spectrometry, and may be more easily recovered than is the case following capture with biotin based reagents.

Direct Observation of Millisecond to Second Motions in Proteins by Dipolar CODEX NMR Spectroscopy

We present a site-resolved study of stow (ms to s) motions in a protein in the solid (microcrystalline) state performed with the use of a modified version of the centerband-only detection of exchange (CODEX) NMR experiment. CODEX was originally based on measuring changes in molecular orientation by means of the chemical shift anisotropy (CSA) tensor, and in our modification, angular reorientations of internuclear vectors are observed. The experiment was applied to the study of stow (15)N-(1)H motions of the SH3 domain of chicken a-spectrin. The protein was perdeuterated with partial back-exchange of protons at labile sites. This allowed indirect (proton) detection of (15)N nuclei and thus a significant enhancement of sensitivity. The diluted proton system also made negligible proton-driven spin diffusion between (15)N nuclei, which interferes with the molecular exchange (motion) and hampers the acquisition of dynamic parameters. The experiment has shown that approximately half of the peaks in the 2D (15)N-(1)H correlation spectrum exhibit exchange in a different extent. The correlation time of the slow motion for most peaks is 1 to 3 s. This is the first NMR study of the internal dynamics of proteins in the solid state on the millisecond to second time scale with site-specific spectral resolution that provides both time-scale and geometry information about molecular motions.