Four-body potentials reveal protein-specific correlations to stability changes caused by hydrophobic core mutations

Mutational experiments show how changes in the hydrophobic cores of proteins affect their stabilities. Here, we estimate these effects computationally, using four-body likelihood potentials obtained by simplicial neighborhood analysis of protein packing (SNAPP). In this procedure, the volume of a known protein structure is tiled with tetrahedra having the center of mass of one amino acid side-chain at each vertex. Log-likelihoods are computed for the 8855 possible tetrahedra with equivalent compositions from structural databases and amino acid frequencies. The sum of these four-body potentials for tetrahedra present in a given protein yields the SNAPP score. Mutations change this sum by changing the compositions of tetrahedra containing the mutated residue and their related potentials. Linear correlation coefficients between experimental mutational stability changes, Delta(DeltaG(unfold)), and those based on SNAPP scoring range from 0.70 to 0.94 for hydrophobic core mutations in five different proteins. Accurate predictions for the effects of hydrophobic core mutations can therefore be obtained by virtual mutagenesis, based on changes to the total SNAPP likelihood potential. Significantly, slopes of the relation between Delta(DeltaG(unfold)) and DeltaSNAPP for different proteins are statistically distinct, and we show that these protein-specific effects can be estimated using the average SNAPP score per residue, which is readily derived from the analysis itself. This result enhances the predictive value of statistical potentials and supports previous suggestions that "comparable" mutations in different proteins may lead to different Delta(DeltaG(unfold)) values because of differences in their flexibility and/or conformational entropy.     

Comparison of protein structural profiles by interactive computer graphics

This paper describes a novel computer graphics tool for predicting protein structures. The method is based on structural profiles; which are plots of hydrophobicity, parameters used for secondary structure prediction, or other residue-specific traits against sequence number. Similar structural profiles can indicate similar tertiary structures, in the absence of sequence homology. The profiles of reference proteins, with known structure, can be used for prediction. In the method presented here, structural profiles are compared by interactive computer graphics, using the program Multiplot. As a test, a structural profile comparison of several proteins known to have similar 3D structures is presented. Comparison of structural profiles detects similar folding of the two domains of rhodanese, which was not easily detected by sequence homology.     

Alanine-scanning mutagenesis of the beta-sheet region of phage T4 lysozyme suggests that tertiary context has a dominant effect on beta-sheet formation

In general, alpha-helical conformations in proteins depend in large part on the amino acid residues within the helix and their proximal interactions. For example, an alanine residue has a high propensity to adopt an alpha-helical conformation, whereas that of a glycine residue is low. The sequence preferences for beta-sheet formation are less obvious. To identify the factors that influence beta-sheet conformation, a series of scanning polyalanine mutations were made within the strands and associated turns of the beta-sheet region in T4 lysozyme. For each construct the stability of the folded protein was reduced substantially, consistent with removal of native packing interactions. However, the crystal structures showed that each of the mutants retained the beta-sheet conformation. These results suggest that the structure of the beta-sheet region of T4 lysozyme is maintained to a substantial extent by tertiary interactions with the surrounding parts of the protein. Such tertiary interactions may be important in determining the structures of beta-sheets in general.     

Molecular dynamics study of the stabilities of consensus designed ankyrin repeat proteins

Two designed ankyrin repeat (AR) proteins (E3_5 and E3_19) are high homologous (with about 87% sequence identity) and their crystal structures have a Calpha atom-positional root-mean-square difference of about 0.14 nm. However, it was found that E3_5 is considerably more stable than E3_19 in guanidinium hydrochloride and thermal denaturation experiments. With the goal of providing insights into the various factors contributing to the stabilities of the designed AR proteins and suggesting possible mutations to enhance their stabilities, homology modeling and molecular dynamics (MD) simulations with explicit solvent have been performed. Because the crystal structure of E3_19 was solved later than that of E3_5, a homology model of E3_19 based on the crystal structure of E3_5 was also used in the simulations. E3_5 shows a very stable trajectory in both crystal and solution simulations. In contrast, the C-terminal repeat of E3_19 unfolds in the simulations starting from either the modeled structure or the crystal structure, although it has a sequence identical to that of E3_5. A continuum electrostatic model was used to estimate the effect of single mutations on protein stability and to study the interaction between the internal ARs and the C-terminal capping AR. Mutations involving charged residues were found to have large effects on stability. Due to the difference in charge distribution in the internal ARs of E3_19 and E3_5, their interaction with the C-terminal capping AR is less favorable in E3_19. The simulation trajectories suggest that the stability of the designed AR proteins can be increased by optimizing the electrostatic interactions within and between the different repeats.     

Incorporating protein conformational flexibility into the calculation of pH-dependent protein properties.

A method for combining calculations of residue pKa's with changes in the position of polar hydrogens has been developed. The Boltzmann distributions of proton positions in hydroxyls and neutral titratable residues are found in the same Monte Carlo sampling procedure that determines the amino acid ionization states at each pH. Electrostatic, Lennard-Jones potentials, and torsion angle energies are considered at each proton position. Many acidic and basic residues are found to have significant electrostatic interactions with either a water- or hydroxyl-containing side chain. Protonation state changes are coupled to reorientation of the neighboring hydroxyl dipoles, resulting in smaller free energy differences between neutral and ionized residues than when the protein is held rigid. Multiconformation pH titration gives better agreement with the experimental pKa's for triclinic hen egg lysozyme than conventional rigid protein calculations. The hydroxyl motion significantly increases the protein dielectric response, making it sensitive to the composition of the local protein structure. More than one conformer per residue is often found at a given pH, providing information about the distribution of low-energy lysozyme structures.     

Structural features that govern enzymatic activity in carbonic anhydrase from a low-temperature adapted fish, Chionodraco hamatus

The carbonic anhydrase (CA) family of zinc metalloenzymes includes many known isozymes that have different subcellular distributions. The study described here focuses on identification of the structural features that define low-temperature adaptation in a Chionodraco hamatus protein, both for the reaction center, at an atomic level, and for the tertiary structure of the protein. To this aim, an x-ray absorption near-edge spectroscopy/Minuit x-ray absorption near-edge spectroscopy analysis of the reaction center was undertaken for both a structurally characterized human CAII and CA of C. hamatus. Higher structural levels were analyzed by sequence comparison and homology modeling. To establish whether the structural insights acquired in fish CAs are general, theoretical models were generated by homology modeling for three temperate-climate-adapted fish CAs. The measured structural differences between the two proteins are discussed in terms of the differences in the electrostatic potential between human CAII and CA of C. hamatus. We conclude that modulation of the interaction between the catalytic water molecule and the zinc ion could depend on the effect of the electrostatic potential distribution.     

农杆菌双组分信号转导系统中VirG蛋白结构与转化功能的分析

农杆菌介导的转化方法在植物、真菌基因转化上有着重要意义,其中VirG蛋白起着重要作用.利用多种序列、结构分析对比软件,对根癌农杆菌和发根农杆菌4种代表性VirG蛋白基本性质、一级结构、二级结构、三级结构和表面静电势进行了预测比较和分析,发现N端差异较大.来自发根农杆菌的VirG蛋白更是在N端缺失了一段氨基酸,这可能是发根农杆菌转化效率低于根癌农杆菌的一个重要原因.结合文献关于定点突变的研究结果,经过总结分析VirG蛋白各位点突变对功能的影响,发现只有诸如N54和I77等少数位点的突变才对功能有促进作用.这些研究结果为提高农杆菌,尤其是发根农杆菌的侵染能力和范围提供了理论依据.     

The Protein Mutant Database.

Currently the protein mutant database (PMD) contains over 81 000 mutants, including artificial as well as natural mutants of various proteins extracted from about 10 000 articles. We recently developed a powerful viewing and retrieving system (http://pmd.ddbj.nig.ac.jp), which is integrated with the sequence and tertiary structure databases. The system has the following features: (i) mutated sequences are displayed after being automatically generated from the information described in the entry together with the sequence data of wild-type proteins integrated. This is a convenient feature because it allows one to see the position of altered amino acids (shown in a different color) in the entire sequence of a wild-type protein; (ii) for those proteins whose 3D structures have been experimentally determined, a 3D structure is displayed to show mutation sites in a different color; (iii) a sequence homology search against PMD can be carried out with any query sequence; (iv) a summary of mutations of homologous sequences can be displayed, which shows all the mutations at a certain site of a protein, recorded throughout the PMD.     

Annotating nucleic acid-binding function based on protein structure

Many of the targets of structural genomics will be proteins with little or no structural similarity to those currently in the database. Therefore, novel function prediction methods that do not rely on sequence or fold similarity to other known proteins are needed. We present an automated approach to predict nucleic-acid-binding (NA-binding) proteins, specifically DNA-binding proteins. The method is based on characterizing the structural and sequence properties of large, positively charged electrostatic patches on DNA-binding protein surfaces, which typically coincide with the DNA-binding-sites. Using an ensemble of features extracted from these electrostatic patches, we predict DNA-binding proteins with high accuracy. We show that our method does not rely on sequence or structure homology and is capable of predicting proteins of novel-binding motifs and protein structures solved in an unbound state. Our method can also distinguish NA-binding proteins from other proteins that have similar, large positive electrostatic patches on their surfaces, but that do not bind nucleic acids.     

Structure of tobacco genes encoding pathogenesis-related proteins from the PR-1 group.

Infection of Samsun NN tobacco with tobacco mosaic virus (TMV) was found to induce the synthesis of mRNA encoding a basic protein with a 67% amino acid sequence homology to the known acidic pathogenesis-related (PR) proteins 1a, 1b and 1c. By Southern blot hybridization it was shown that the tobacco genome contains at least eight genes for acidic PR-1 proteins and a similar number of genes encoding the basic homologues. Clones corresponding to three of the genes for acidic PR-1 proteins were isolated from a genomic library of Samsun NN tobacco. The nucleotide sequence of these genes and their flanking sequences were determined. One clone was found to correspond to the PR-1a gene; the two other clones do not correspond to known TMV-induced PR-1 mRNA's and may represent silent genes. Compared to the PR-1a gene, these genes contain an insertion or deletion in the putative promoter region and mutations affecting the PR-1 reading frame.     

基于SwissPdbViewer和MATLAB对东北虎Sox蛋白三级结构建模及分析

SOX蛋白具有一个与DNA特异结合的高保守HMG-box结构域。为研究东北虎SOX蛋白三级结构的分子机理,利用MATLAB的Bioinformatics工具从GenBank中下载东北虎SOX蛋白序列信息,以三级结构已知的SOX2为模板,联合SwissPdbViewer与MATLAB,采用同源建模方法对SOX蛋白HMG-box进行建模、预测;利用MATLAB的Visualization Tool分析预测结果的三维结构。结果显示PtSox蛋白的HMG-box由3个α-螺旋和2个loop区构成;热稳定性分析表明PtSox蛋白loop区的热力学结构不稳定;表面静电分布显示出PtSox蛋白C-端的中间有一个可能与其它小分子或蛋白质的相互作用位点的N/C腔,上述空间结构可能与其活性与功能的调控有关。     

Sequence-structure specificity of a knowledge based energy function at the secondary structure level

MOTIVATION: This paper investigates the sequence-structure specificity of a representative knowledge based energy function by applying it to threading at the level of secondary structures of proteins. Assessing the strengths and weaknesses of an energy function at this fundamental level provides more detailed and insightful information than at the tertiary structure level and the results obtained can be useful in tertiary level threading. RESULTS: We threaded each of the 293 non-redundant proteins onto the secondary structures contained in its respective native protein (host template). We also used 68 pairs of proteins with similar folds and low sequence identity. For each pair, we threaded the sequence of one protein onto the secondary structures of the other protein. The discerning power of the total energy function and its one-body, pairwise, and mutation components is studied. We then applied our energy function to a recent study which demonstrated how a designed 11-amino acid sequence can replace distinct segments (one segment is an alpha-helix, the other is a beta-sheet) of a protein without changing its fold. We conducted random mutations of the designed sequence to determine the patterns for favorable mutations. We also studied the sequence-structure specificity at the boundaries of a secondary structure. Finally, we demonstrated how to speed up tertiary level threading by filtering out alignments found to be energetically unfavorable during the secondary structure threading. AVAILABILITY: The program is available on request from the authors. CONTACT: xud@ornl.gov     

Knowledge based modelling of homologous proteins, Part I: Three-dimensional frameworks derived from the simultaneous superposition of multiple structures

An approach is described for modelling the three-dimensional structure of a protein from the tertiary structures of several homologous proteins that have been determined by X-ray analysis. A method is developed for the simultaneous superposition of several protein molecules and for the calculation of an 'average structure' or 'framework'. Investigation of the convergence properties of this method, in the case of both weighted and unweighted least squares, demonstrates that both give a unique answer and the latter is robust for an homologous family of proteins. Multi-dimensional scaling is used to subgroup of the proteins with respect to structural homology. The framework calculated on the basis of the family of homologous proteins, or of an appropriate subgroup, is used to align fragments of the known protein structures of high sequence homology with the unknown. This alignment provides a basis for model building the tertiary structure. Different techniques for using the framework to model the mainchain of various globins and an immunoglobulin domain in the structurally conserved regions are investigated.     

Mth10b, a unique member of the Sac10b family, does not bind nucleic acid

The Sac10b protein family is regarded as a group of nucleic acid-binding proteins that are highly conserved and widely distributed within archaea. All reported members of this family are basic proteins that exist as homodimers in solution and bind to DNA and/or RNA without apparent sequence specificity in vitro. Here, we reported a unique member of the family, Mth10b from Methanobacterium thermoautotrophicum ΔH, whose amino acid sequence shares high homology with other Sac10b family proteins. However, unlike those proteins, Mth10b is an acidic protein; its potential isoelectric point is only 4.56, which is inconsistent with the characteristics of a nucleic acid-binding protein. In this study, Mth10b was expressed in Escherichia coli and purified using a three-column chromatography purification procedure. Biochemical characterization indicated that Mth10b should be similar to typical Sac10b family proteins with respect to its secondary and tertiary structure and in its preferred oligomeric forms. However, an electrophoretic mobility shift analysis (EMSA) showed that neither DNA nor RNA bound to Mth10b in vitro, indicating that either Mth10b likely has a physiological function that is distinct from those of other Sac10b family members or nucleic acid-binding ability may not be a fundamental factor to the actual function of the Sac10b family.     

Computer modeling studies of the structure of a repressor

Due to advances in molecular biology the DNA sequences of structural genes coding for proteins are often known before a protein is characterized or even isolated. The function of a protein whose amino acid sequence has been deduced from a DNA sequence may not even be known. This has created greater interest in the development of methods to predict the tertiary structures of proteins. The a priori prediction of a protein's structure from its amino acid sequence is not yet possible. However, since proteins with similar amino acid sequences are observed to have similar three-dimensional structures, it is possible to use an analogy with a protein of known structure to draw some conclusions about the structure and properties of an uncharacterized protein. The process of predicting the tertiary structure of a protein relies very much upon computer modeling and analysis of the structure. The prediction of the structure of the bacteriophage 434 cro repressor is used as an example illustrating current procedures.     

Search for fucose binding domains in recently sequenced hypothetical proteins using molecular modeling techniques and structural analysis

The crystal structure of a fucose-binding lectin from the bacteria Pseudomonas aeruginosa in complex with α-L-fucose has been recently determined. It is a tetramer; each monomer displays a nine-stranded, antiparallel, β-sandwiched arrangement and contains two calcium ions that mediate the binding of fucose in a recognition mode unique among protein-carbohydrate interactions. In search of this type of unique interactions in other newly discovered protein sequences, we have used molecular modeling techniques to predict and analyze the 3-D structures of some proteins, which exhibited reasonable degree of homology with the amino acid sequence of the bacterial protein. A BLAST search with the sequence of Pseudomonas aeruginosa as query in the non-redundant sequence database identified four proteins from different species, three organisms from bacteria and one from archaea. We have modeled the structures of these proteins as well as those of the complexes with carbohydrates and studied the nature of physicochemical forces involved in the complex formation both in presence and absence of calcium. The calcium-binding loops have been found to be highly conserved both in terms of primary and tertiary structures in these proteins, although a less acidic character is observed in Photorhabdus lectin due to the absence of two aspartic acid residues on the calcium-binding loop which also resulted in lower binding affinity. All these structures exhibited highly negative electrostatic environment in the vicinity of the calcium-binding loops which was essential for neutralizing the positive charges of two closely situated Ca⁺² ions. The comparison of the binding affinities of some monosaccharides other than fucose, e.g. mannose and fructose, showed higher binding energies confirming the fucose specificity of these proteins.     

Solution structure of switch Arc,a mutant with 3(10) helices replacing a wild-type beta-ribbon

Adjacent N11L and L12N mutations in the antiparallel beta-ribbon of Arc repressor result in dramatic changes in local structure in which each beta-strand is replaced by a right-handed helix. The full solution structure of this "switch" Arc mutant shows that irregular 3(10) helices compose the new secondary structure. This structural metamorphosis conserves the number of main-chain and side-chain to main-chain hydrogen bonds and the number of fully buried core residues. Apart from a slight widening of the interhelical angle between alpha-helices A and B and changes in side-chain conformation of a few core residues in Arc, no large-scale structural adjustments in the remainder of the protein are necessary to accommodate the ribbon-to-helix change. Nevertheless, some changes in hydrogen-exchange rates are observed, even in regions that have very similar structures in the two proteins. The surface of switch Arc is packed poorly compared to wild-type, leading to approximately 1000A(2) of additional solvent-accessible surface area, and the N termini of the 3(10) helices make unfavorable head-to-head electrostatic interactions. These structural features account for the positive m value and salt dependence of the ribbon-to-helix transition in Arc-N11L, a variant that can adopt either the mutant or wild-type structures. The tertiary fold is capped in different ways in switch and wild-type Arc, showing how stepwise evolutionary transformations can arise through small changes in amino acid sequence.     

Evolutionary conservativeness of electric field in the Cu,Zn superoxide dismutase active site. Evidence for co-ordinated mutation of charged amino acid residues.

Equipotential lines were calculated, using the Poisson-Boltzmann equation, for six Cu,Zn superoxide dismutases with different protein electric charge and various degrees of sequence homology, namely those from ox, pig, sheep, yeast, and the isoenzymes A and B from the amphibian Xenopus laevis. The three-dimensional structures of the porcine and ovine superoxide dismutases were obtained by molecular modelling reconstruction using the structure of the highly homologous bovine enzyme as a template. The three-dimensional structure of the evolutionary distant yeast Cu,Zn superoxide dismutase was recently resolved by us, while computer-modelled structures are available for X. laevis isoenzymes. The six proteins display large differences in the net protein charge and distribution of electrically charged surface residues but the trend of the equipotential lines in the proximity of the active sites was found to be constant in all cases. These results are in line with the very similar catlytic rate constants experimentally measured for the corresponding enzyme activities. This analysis shows that electrostatic guidance for the enzyme-substrate interaction in Cu,Zn superoxide dismutases is related to a spatial distribution of charges, arranged so as to maintain, in the area surrounding the active sites, an identical electrostatic potential distribution, which is conserved in the evolution of this protein family.     

Prediction of protein stability changes for single-site mutations using support vector machines

Accurate prediction of protein stability changes resulting from single amino acid mutations is important for understanding protein structures and designing new proteins. We use support vector machines to predict protein stability changes for single amino acid mutations leveraging both sequence and structural information. We evaluate our approach using cross-validation methods on a large dataset of single amino acid mutations. When only the sign of the stability changes is considered, the predictive method achieves 84% accuracy-a significant improvement over previously published results. Moreover, the experimental results show that the prediction accuracy obtained using sequence alone is close to the accuracy obtained using tertiary structure information. Because our method can accurately predict protein stability changes using primary sequence information only, it is applicable to many situations where the tertiary structure is unknown, overcoming a major limitation of previous methods which require tertiary information. The web server for predictions of protein stability changes upon mutations (MUpro), software, and datasets are available at http://www.igb.uci.edu/servers/servers.html.     

Redox properties of Arabidopsis cytochrome c6 are independent of the loop extension specific to higher plants

Cytochrome c6 (cytc6) from Arabidopsis differs from the cyanobacterial and algal homologues in several redox properties. It is possible that these differences might be due to the presence of a 12 amino acid residue loop extension common to higher plant cytc6 proteins. However, homology modelling suggests this is not the case. We report experiments to test if differences in biochemical properties could be due to this extension. Analysis of mutant forms of Arabidopsis cytc6 in which the entire extension was lacking, or a pair of cysteine residues in the extension had been exchanged for serine, revealed no significant effect of these changes on either the redox potential of the haem group or the reactivity towards Photosystem I (PSI). We conclude that the differences in properties are due to more subtle unidentified differences in structure, and that the sequence extension in the higher plant proteins has a function yet to be identified.