Crystal structure of a hypothetical protein,TM841 of Thermotoga maritima,reveals its function as a fatty acid-binding protein

We determined the three-dimensional (3D) crystal structure of protein TM841, a protein product from a hypothetical open-reading frame in the genome of the hyperthermophile bacterium Thermotoga maritima, to 2.0 A resolution. The protein belongs to a large protein family, DegV or COG1307 of unknown function. The 35 kDa protein consists of two separate domains, with low-level structural resemblance to domains from other proteins with known 3D structures. These structural homologies, however, provided no clues for the function of TM841. But the electron density maps revealed clear density for a bound fatty-acid molecule in a pocket between the two protein domains. The structure indicates that TM841 has the molecular function of fatty-acid binding and may play a role in the cellular functions of fatty acid transport or metabolism.     

Efficient RMSD measures for the comparison of two molecular ensembles. Root-mean-square deviation

Quantitative measures are presented for comparing the conformations of two molecular ensembles. The measures are based on Kabsch's formula for the root-mean-square deviation (RMSD) and the covariance matrix of atomic positions of isotropically distributed ensembles (IDE). By using a Taylor series expansion, it is shown that the RMSD can be expressed solely in terms of the IDE matrices. A fast approximate method is introduced for the pairwise RMSD determination whose computational cost scales linearly with the number of structures. A similarity measure for two structural ensembles that is based on the trace metric of the differences of powers of the IDE matrices is presented. The measures are illustrated for conformational ensembles generated by a molecular dynamics computer simulation of a partially folded A-state analog of ubiquitin.     

Determination of molecular alignment tensors without backbone resonance assignment: Aid to rapid analysis of protein-protein interactions

Based on high-resolution structures of the free molecules accurate determination of structures of protein complexes by NMR spectroscopy is possible using residual dipolar couplings. In order, however, to be able to apply these methods, protein backbone resonances have to be assigned first. This NMR assignment process is particularly difficult and time consuming for protein sizes above 20 kDa. Here we show that, when NMR resonances belonging to a specific amino acid type are selected either by amino acid specific labeling, by their characteristic C/C chemical shifts or by dedicated NMR experiments, molecular alignment tensors of proteins up to 80 kDa can be determined without prior backbone resonance assignment. This offers the opportunity to greatly accelerate determination of three-dimensional structures of protein-protein and protein-ligand complexes, and validation of multimeric states of proteins. Moreover, exhaustive back-calculation can be performed using only ¹D_NH couplings. Therefore, it avoids expensive ¹³C-labeling and it gives access to orientational information for large proteins that strongly aggregate at concentrations above 50 M, i.e., experimental conditions where 3D triple resonance experiments are not sensitive enough to allow backbone resonance assignment.     

On the evaluation and optimization of protein X-ray structures for pKa calculations

The calculation of the physical properties of a protein from its X-ray structure is of importance in virtually every aspect of modern biology. Although computational algorithms have been developed for calculating everything from the dynamics of a protein to its binding specificity, only limited information is available on the ability of these methods to give accurate results when used with a particular X-ray structure. We examine the ability of a pKa calculation algorithm to predict the proton-donating residue in the catalytic mechanism of hen egg white lysozyme. We examine the correlation between the ability of the pKa calculation method to obtain the correct result and the overall characteristics of 41 X-ray structures such as crystallization conditions, resolution, and the output of structure validation software. We furthermore examine the ability of energy minimizations (EM), molecular dynamics (MD) simulations, and structure-perturbation methods to optimize the X-ray structures such that these give correct results with the pKa calculation algorithm. We propose a set of criteria for identifying the proton donor in a catalytic mechanism, and demonstrate that the application of these criteria give highly accurate prediction results when using unmodified X-ray structures. More specifically, we are able to successfully identify the proton donor in 85% of the X-ray structures when excluding structures with crystal contacts near the active site. Neither the use of the overall characteristics of the X-ray structures nor the optimization of the structure by EM, MD, or other methods improves the results of the pKa calculation algorithm. We discuss these results and their implications for the design of structure-based energy calculation algorithms in general.     

Solution structure of Vibrio cholerae protein VC0424: a variation of the ferredoxin-like fold

The structure of Vibrio cholerae protein VC0424 was determined by NMR spectroscopy. VC0424 belongs to a conserved family of bacterial proteins of unknown function (COG 3076). The structure has an alpha-beta sandwich architecture consisting of two layers: a four-stranded antiparallel beta-sheet and three side-by-side alpha-helices. The secondary structure elements have the order alphabetaalphabetabetaalphabeta along the sequence. This fold is the same as the ferredoxin-like fold, except with an additional long N-terminal helix, making it a variation on this common motif. A cluster of conserved surface residues on the beta-sheet side of the protein forms a pocket that may be important for the biological function of this conserved family of proteins.     

Relationship between local structural entropy and protein thermostability

We developed a technique to compute structural entropy directly from protein sequences. We explored the possibility of using structural entropy to identify residues involved in thermal stabilization of various protein families. Examples include methanococcal adenylate kinase, Ribonuclease HI and holocytochrome c(551). Our results show that the positions of the largest structural entropy differences between wild type and mutant usually coincide with the residues relevant to thermostability. We also observed a good linear relationship between the average structural entropy and the melting temperatures for adenylate kinase and its chimeric constructs. To validate this linear relationship, we compiled a large dataset comprised of 1153 sequences and found that most protein families still display similar linear relationships. Our results suggest that the multitude of interactions involved in thermal stabilization may be generalized into the tendency of proteins to maintain local structural conservation. The linear relationship between structural entropy and protein thermostability should be useful in the study of protein thermal stabilization.     

Crystal structure of the YGR205w protein from Saccharomyces cerevisiae: close structural resemblance to E. coli pantothenate kinase

The protein product of the YGR205w gene of Saccharomyces cerevisiae was targeted as part of our yeast structural genomics project. YGR205w codes for a small (290 amino acids) protein with unknown structure and function. The only recognizable sequence feature is the presence of a Walker A motif (P loop) indicating a possible nucleotide binding/converting function. We determined the three-dimensional crystal structure of Se-methionine substituted protein using multiple anomalous diffraction. The structure revealed a well known mononucleotide fold and strong resemblance to the structure of small metabolite phosphorylating enzymes such as pantothenate and phosphoribulo kinase. Biochemical experiments show that YGR205w binds specifically ATP and, less tightly, ADP. The structure also revealed the presence of two bound sulphate ions, occupying opposite niches in a canyon that corresponds to the active site of the protein. One sulphate is bound to the P-loop in a position that corresponds to the position of beta-phosphate in mononucleotide protein ATP complex, suggesting the protein is indeed a kinase. The nature of the phosphate accepting substrate remains to be determined.     

Solution structure of the highly acidic protein HI1450 from Haemophilus influenzae, a putative double-stranded DNA mimic

The solution structure of the acidic protein HI1450 from Haemophilus influenzae has been determined by NMR spectroscopy. HI1450 has homologues in ten other bacterial species including Escherichia coli, Vibrio cholerae, and Yersinia pestis but there are no functional assignments for any members of the family. Thirty-one of the amino acids in this 107-residue protein are aspartates or glutamates, contributing to an unusually low pI of 3.72. The secondary structure elements are arranged in an alpha-alpha-beta-beta-beta-beta order with the two alpha helices packed against the same side of an anti-parallel four-stranded beta meander. Two large loops, one between beta1 and beta2 and the other between beta2 and beta3 bend almost perpendicularly across the beta-strands in opposite directions on the non-helical side of the beta-sheet to form a conserved hydrophobic cavity. The HI1450 structure has some similarities to the structure of the double-stranded DNA (dsDNA) mimic uracil DNA glycosylase inhibitor (Ugi) including the distribution of surface charges and the position of the hydrophobic cavity. Based on these similarities, as well as having a comparable molecular surface to dsDNA, we propose that HI1450 may function as a dsDNA mimic in order to inhibit or regulate an as yet unidentified dsDNA binding protein.