SARS冠状病毒M蛋白的生物信息学研究

针对GenBank上发布的来自不同国家地区的39条SARSCoV推测M蛋白,采用生物信息学软件分析其核酸和氨基酸序列,获得其分子生物学特征,确定突变位点,预测功能结构区、Motif及抗原决定簇,比较基因突变对这些功能结构的影响.结果表明:在39个病毒株M蛋白的666 bp中,共有18个病毒株在7个位点上发生了25次变异.在M蛋白序列上预测获得3个跨膜螺旋序列和一个可能的信号肽序列.氨基酸序列的变异主要发生在其跨膜和胞外区域,胞内区域相对较少.预测发现12个Motif和7个抗原决定簇.提示突变对M蛋白的结构功能区的影响不大,也未造成M蛋白的Motif的数量和构成发生改变.对抗原决定簇的影响也主要体现在序列成分构成的改变上,在设计疫苗时,应考虑由其导致的抗原特性改变.     

Interplay between metal binding and cis/trans isomerization in legume lectins: structural and thermodynamic study of P. angolensis lectin

The interplay between metal binding, carbohydrate binding activity, stability and structure of the lectin from Pterocarpus angolensis was investigated. Removal of the metals leads to a more flexible form of the protein with significantly less conformational stability. Crystal structures of this metal-free form show significant structural rearrangements, although some structural features that allow the binding of sugars are retained. We propose that substitution of an asparagine residue at the start of the C-terminal beta-strand of the legume lectin monomer hinders the trans-isomerization of the cis-peptide bond upon demetallization and constitutes an intramolecular switch governing the isomer state of the non-proline bond and ultimately the lectin phenotype.     

The periplasmic bacterial molecular chaperone SurA adapts its structure to bind peptides in different conformations to assert a sequence preference for aromatic residues

The periplasmic molecular chaperone protein SurA facilitates correct folding and maturation of outer membrane proteins in Gram-negative bacteria. It preferentially binds peptides that have a high fraction of aromatic amino acids. Phage display selections, isothermal titration calorimetry and crystallographic structure determination have been used to elucidate the basis of the binding specificity. The peptide recognition is imparted by the first peptidyl-prolyl isomerase (PPIase) domain of SurA. Crystal structures of complexes between peptides of sequence WEYIPNV and NFTLKFWDIFRK with the first PPIase domain of the Escherichia coli SurA protein at 1.3 A resolution, and of a complex between the dodecapeptide and a SurA fragment lacking the second PPIase domain at 3.4 A resolution, have been solved. SurA binds as a monomer to the heptapeptide in an extended conformation. It binds as a dimer to the dodecapeptide in an alpha-helical conformation, predicated on a substantial structural rearrangement of the SurA protein. In both cases, side-chains of aromatic residues of the peptides contribute a large fraction of the binding interactions. SurA therefore asserts a recognition preference for aromatic amino acids in a variety of sequence configurations by adopting alternative tertiary and quaternary structures to bind peptides in different conformations.     

Design and NMR conformational study of a beta-sheet peptide based on Betanova and WW domains

A good approach to test our current knowledge on formation of protein beta-sheets is de novo protein design. To obtain a three-stranded beta-sheet mini-protein, we have built a series of chimeric peptides by taking as a template a previously designed beta-sheet peptide, Betanova-LLM, and incorporating N- and/or C-terminal extensions taken from WW domains, the smallest natural beta-sheet domain that is stable in absence of disulfide bridges. Some Betanova-LLM strand residues were also substituted by those of a prototype WW domain. The designed peptides were cloned and expressed in Escherichia coli. The ability of the purified peptides to adopt beta-sheet structures was examined by circular dichroism (CD). Then, the peptide showing the highest beta-sheet population according to the CD spectra, named 3SBWW-2, was further investigated by 1H and 13C NMR. Based on NOE and chemical shift data, peptide 3SBWW-2 adopts a well defined three-stranded antiparallel beta-sheet structure with a disordered C-terminal tail. To discern between the contributions to beta-sheet stability of strand residues and the C-terminal extension, the structural behavior of a control peptide with the same strand residues as 3SBWW-2 but lacking the C-terminal extension, named Betanova-LYYL, was also investigated. beta-Sheet stability in these two peptides, in the parent Betanova-LLM and in WW-P, a prototype WW domain, decreased in the order WW-P > 3SBWW-2 > Betanova-LYYL > Betanova-LLM. Conclusions about the contributions to beta-sheet stability were drawn by comparing structural properties of these four peptides.     

Chemical cleavage of fusion proteins for high-level production of transmembrane peptides and protein domains containing conserved methionines

Due to their high hydrophobicity, it is a challenge to obtain high yields of transmembrane peptides for structural and functional characterization. In the present work, a robust method is developed for the expression, purification and reconstitution of transmembrane peptides, especially for those containing conserved methionines. By using a truncated and mutated glutathione-S-transferase construct as the carrier protein and hydroxylamine (which specifically cleaves the peptide bond between Asn and Gly) as the cleavage reagent, 10 mg of the first transmembrane helix of CorA, a Mg²⁺ transporter from Mycobacterium tuberculosis, can be conveniently obtained with high purity from 1 L of M9 minimal media under optimized conditions. The biophysical properties of the peptide were studied by circular dichroism and nuclear magnetic resonance spectroscopy, and the results show that this CorA peptide is well folded in detergent micelles and the secondary structure is very similar to that in recent crystal structures. In addition, this CorA construct is oligomeric in perfluoro-octanoic acid micelles. The compatibility with the transmembrane peptides containing conserved methionines, the high yield and the simple process make the present method competitive with other commonly used methods to produce such peptides for structural and functional studies.     

Structural investigation of cold activity and regulation of aspartate carbamoyltransferase from the extreme psychrophilic bacterium Moritella profunda

Aspartate carbamoyltransferase (EC 2.1.3.2) is extensively studied as a model for cooperativity and allosteric regulation. The structure of the Escherichia coli enzyme has been thoroughly analyzed by X-ray crystallography, and recently the crystal structures of two hyperthermophilic ATCases of the same structural class have been characterized. We here report the detailed functional and structural investigation of the ATCase from the psychrophilic deep sea bacterium Moritella profunda. Our analysis indicates that the enzyme conforms to the E. coli model in that two allosteric states exist that are influenced by similar homotropic interactions. The heterotropic properties differ in that CTP and UTP inhibit the holoenzyme, but ATP seems to exhibit a dual regulatory pattern, activating the enzyme at low concentrations and inhibiting it in the mM range. The crystal structure of the unliganded M. profunda ATCase shows resemblance to a more extreme T state reported previously for an E. coli ATCase mutant. A detailed molecular analysis reveals potential features of adaptation to cold activity and cold regulation. Moreover, M. profunda ATCase presents similarities with certain mutants of E. coli ATCase altered in their kinetic properties or temperature relationships. Finally, structural and functional comparison of ATCases across the full physiological temperature range agrees with an important, but fundamentally different role for electrostatics in protein adaptation at both extremes, i.e. an increased stability through the formation of ion pairs and ion pair networks at high physiological temperatures, and an increased flexibility through enhanced protein solvation at low temperatures.     

Modularity of the hydrophobic core and evolution of functional diversity in fold A glycosyltransferases

Hydrophobic cores are fundamental structural properties of proteins typically associated with protein folding and stability; however, how the hydrophobic core shapes protein evolution and function is poorly understood. Here, we investigated the role of conserved hydrophobic cores in fold-A glycosyltransferases (GT-As), a large superfamily of enzymes that catalyze formation of glycosidic linkages between diverse donor and acceptor substrates through distinct catalytic mechanisms (inverting versus retaining). Using hidden Markov models and protein structural alignments, we identify similarities in the phosphate-binding cassette (PBC) of GT-As and unrelated nucleotide-binding proteins, such as UDP-sugar pyrophosphorylases. We demonstrate that GT-As have diverged from other nucleotide-binding proteins through structural elaboration of the PBC and its unique hydrophobic tethering to the F-helix, which harbors the catalytic base (xED-Asp). While the hydrophobic tethering is conserved across diverse GT-A fold enzymes, some families, such as B3GNT2, display variations in tethering interactions and core packing. We evaluated the structural and functional impact of these core variations through experimental mutational analysis and molecular dynamics simulations and find that some of the core mutations (T336I in B3GNT2) increase catalytic efficiency by modulating the conformational occupancy of the catalytic base between “D-in” and acceptor-accessible “D-out” conformation. Taken together, our studies support a model of evolution in which the GT-A core evolved progressively through elaboration upon an ancient PBC found in diverse nucleotide-binding proteins, and malleability of this core provided the structural framework for evolving new catalytic and substrate-binding functions in extant GT-A fold enzymes.     

Structure prediction and evolution of a halo-acid dehalogenase of Burkholderia mallei

Environmental pollutants containing halogenated organic compounds e.g. haloacid, can cause a plethora of health problems. The structural and functional analyses of the gene responsible of their degradation are an important aspect for environmental studies and are important to human well-being. It has been shown that some haloacids are toxic and mutagenic. Microorganisms capable of degrading these haloacids can be found in the natural environment. One of these, a soil-borne Burkholderia mallei posses the ability to grow on monobromoacetate (MBA). This bacterium produces a haloacid dehalogenase that allows the cell to grow on MBA, a highly toxic and mutagenic environmental pollutant. For the structural and functional analysis, a 346 amino acid encoding protein sequence of haloacid dehalogenase is retrieve from NCBI data base. Primary and secondary structure analysis suggested that the high percentage of helices in the structure makes the protein more flexible for folding, which might increase protein interactions. The consensus protein sub-cellular localization predictions suggest that dehalogenase protein is a periplasmic protein 3D2GO server, suggesting that it is mainly employed in metabolic process followed by hydrolase activity and catalytic activity. The tertiary structure of protein was predicted by homology modeling. The result suggests that the protein is an unstable protein which is also an important characteristic of active enzyme enabling them to bind various cofactors and substrate for proper functioning. Validation of 3D structure was done using Ramachandran plot ProsA-web and RMSD score. This predicted information will help in better understanding of mechanism underlying haloacid dehalogenase encoding protein and its evolutionary relationship.     

Coupling efficiency of rhodopsin and transducin in bicelles

G protein coupled receptors (GPCRs) can be activated by various extracellular stimuli, including hormones, peptides, odorants, neurotransmitters, nucleotides, or light. After activation, receptors interact with heterotrimeric G proteins and catalyze GDP release from the Gα subunit, the rate limiting step in G protein activation, to form a high affinity nucleotide-free GPCR-G protein complex. In vivo, subsequent GTP binding reduces affinity of the Gα protein for the activated receptor. In this study, we investigated the biochemical and structural characteristics of the prototypical GPCR, rhodopsin, and its signaling partner, transducin (G(t)), in bicelles to better understand the effects of membrane composition on high affinity complex formation, stability, and receptor mediated nucleotide release. Our results demonstrate that the high-affinity complex (rhodopsin-G(t)(empty)) forms more readily and has dramatically increased stability when rhodopsin is integrated into bicelles of a defined composition. We increased the half-life of functional complex to 1 week in the presence of negatively charged phospholipids. These data suggest that a membrane-like structure is an important contributor to the formation and stability of functional receptor-G protein complexes and can extend the range of studies that investigate properties of these complexes.     

Domain atrophy creates rare cases of functional partial protein domains

BackgroundProtein domains display a range of structural diversity, with numerous additions and deletions of secondary structural elements between related domains. We have observed a small number of cases of surprising large-scale deletions of core elements of structural domains. We propose a new concept called domain atrophy, where protein domains lose a significant number of core structural elements.ResultsHere, we implement a new pipeline to systematically identify new cases of domain atrophy across all known protein sequences. The output of this pipeline was carefully checked by hand, which filtered out partial domain instances that were unlikely to represent true domain atrophy due to misannotations or un-annotated sequence fragments. We identify 75 cases of domain atrophy, of which eight cases are found in a three-dimensional protein structure and 67 cases have been inferred based on mapping to a known homologous structure. Domains with structural variations include ancient folds such as the TIM-barrel and Rossmann folds. Most of these domains are observed to show structural loss that does not affect their functional sites.ConclusionOur analysis has significantly increased the known cases of domain atrophy. We discuss specific instances of domain atrophy and see that there has often been a compensatory mechanism that helps to maintain the stability of the partial domain. Our study indicates that although domain atrophy is an extremely rare phenomenon, protein domains under certain circumstances can tolerate extreme mutations giving rise to partial, but functional, domains.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0655-8) contains supplementary material, which is available to authorized users.     

Evidence of colorectal cancer risk associated variant Lys25Ser in the proximity of human bone morphogenetic protein 2

Colorectal cancer (CRC) is the third most prevalent cancer and fourth leading cause of cancer-related deaths globally. It has been shown that the nsSNP variants play an important role in diseases, however it remained unclear how these variants are associated with the disease. Recently, several CRC risk associated SNPs have been discovered, however rs961253 (Lys25Arg at 20p12.3) located in the proximity of bone morphogenetic protein 2 (Bmp2) and fermitin family homolog 1 Fermt1 genes have been reported to be highly associated with the CRC risk. Here we provide evidence for the first time in silico biological functional and structural implications of non-synonymous (nsSNPs) CRC disease-associated variant Lys25Arg via molecular dynamic (MD) simulation. Protein structural analysis was performed with a particular variant allele (A/C, Lys25Arg) and compared with the predicted native protein structure. Our results showed that this nsSNP will cause changes in the protein structure and as a result is associated with the disease. In addition to the native and mutant 3D structures of CRC associated risk allele protein domain (CRAPD), they were also analyzed using solvent accessibility models for further protein stability confirmation. Taken together, this study confirmed that this variant has functional effect and structural impact on the CRAPD and may play an important role in CRC disease progression; hence it could be a reasonable approach for studying the effect of other deleterious variants in future studies.     

Structural Insights into Genetic Variants of Na+/Glucose Cotransporter SGLT1 Causing Glucose–Galactose Malabsorption: vSGLT as a Model Structure

Current advances in structural biology provide valuable insights into structure-function relationship of membrane transporters by solving crystal structures of bacterial homologs of human transporters. Therefore, scientists consider bacterial transporters as useful structural models for designing of drugs targeted in human diseases. The functional homology between Vibrio parahaemolyticus Na(+)/galactose transporter (vSGLT) and Na(+)/glucose cotransporter SGLT1 has been well established a decade ago. Now the crystal structure of vSGLT is considered quite valuable in explaining not only the cotransport mechanisms, but it also acts as a representative protein in understanding the protein stability and amino acid interactions within the core structure. We investigated the molecular mechanisms of genetic variations in SGLT1 that cause glucose-galactose malabsorption (GGM) defects using the crystal structure of vSGLT as a model sugar transporter. Our in silico mutagenesis and modeling analysis suggest that the GGM genetic variations lead to conformational changes either by structure destabilization or by formation of unnecessary interaction within the core structure of SGLT1 thereby explaining the genetic defects in Na(+) dependent sugar translocation across the cell membrane.     

Aromatic cross-strand ladders control the structure and stability of beta-rich peptide self-assembly mimics

Though β-rich self-assemblies comprise a major structural class of polypeptides, a detailed understanding of the determinants of their structure and stability is lacking. In particular, the roles of repetitive stretches of side chains running the long axis of these β-sheets, termed “cross-strand ladders,” remain poorly characterized due to the inherently insoluble and heterogeneous nature of self-assemblies. To overcome these experimental challenges, we have established a complementary experimental system termed “peptide self-assembly mimics” (PSAMs). The PSAMs capture a defined number of self-assembly-like peptide repeats within a soluble β-rich protein, making structural and energetic studies possible. In this work, we investigated the role of cross-strand ladders containing aromatic residues, which are prominent in self-assembling peptides. A combination of solution data and high-resolution crystal structures revealed that a single cross-strand ladder consisting solely of Tyr significantly stabilized, rigidified, and flattened the PSAM β-sheet. These characteristics would stabilize each β-sheet layer of a self-assembly and direct sheet conformations compatible with lamination. Our results therefore provide a rationale for the abundance of aromatic amino acids in fibril-forming peptides and establish important roles of cross-strand Tyr ladders in the structure and stability of β-rich peptide self-assemblies.     

Secondary structure and dynamics of micelle bound beta- and gamma-synuclein

We have used solution state NMR spectroscopy to characterize the secondary structure and backbone dynamics of the proteins beta- and gamma-synuclein in their detergent micelle-bound conformations. Comparison of the results with those previously obtained for the Parkinson's disease-linked protein alpha-synuclein shows that structural differences between the three homologous synuclein family members are directly related to variations in their primary amino acid sequences. An 11-residue deletion in the lipid-binding domain of beta-synuclein leads to the destabilization of an entire segment of the micelle-bound helical structure containing the deletion site. The acidic C-terminal tail region of gamma-synuclein, which displays extensive sequence divergence, is more highly disordered than the corresponding regions in the other two family members. The observed structural differences are likely to mediate functional variations between the three proteins, with differences between alpha- and beta-synuclein expected to revolve around their lipid interactions, while differences in gamma-synuclein function are expected to result from different protein-protein interactions mediated by its unique C-terminal tail.     

Thermal transitions in serum amyloid A in solution and on the lipid: implications for structure and stability of acute-phase HDL

Serum amyloid A (SAA) is an acute-phase protein that circulates mainly on plasma HDL. SAA interactions with its functional ligands and its pathogenic deposition in reactive amyloidosis depend, in part, on the structural disorder of this protein and its propensity to oligomerize. In vivo, SAA can displace a substantial fraction of the major HDL protein, apoA-I, and thereby influence the structural remodeling and functions of acute-phase HDL in ways that are incompletely understood. We use murine SAA1.1 to report the first structural stability study of human plasma HDL that has been enriched with SAA. Calorimetric and spectroscopic analyses of these and other SAA-lipid systems reveal two surprising findings. First, progressive displacement of the exchangeable fraction of apoA-I by SAA has little effect on the structural stability of HDL and its fusion and release of core lipids. Consequently, the major determinant for HDL stability is the nonexchangeable apoA-I. A structural model explaining this observation is proposed, which is consistent with functional studies in acute-phase HDL. Second, we report an α-helix folding/unfolding transition in SAA in the presence of lipid at near-physiological temperatures. This new transition may have potentially important implications for normal functions of SAA and its pathogenic misfolding.     

The binding ability of insulin-related peptides as the clue to the search for their functional role in phylogenesis. Introduction to the problem

The common plan of structure of the main peptides of the vertebrate insulin family—insulin itself, IGF-I, IGF-II, and relaxin—has distinct structural formations. Each of the peptides performs its characteristic function. However, overlapping of insulin and IGF-I actions and its stability in the vertebrate phylogenesis have formed the concept of their regulation of growth and metabolism as a function fixed in phylogenesis for a certain type of structures. At the same time, study of insulin-related peptides in invertebrates has revealed the wider spectrum, than in vertebrates, of biological effects; this indicated that the similarity of the total structure design is not sufficient for judging about their functional role. Functional possibilities of a regulatory peptide depend fundamentally on its capability for binding to the receptor realizing its biological action. However, the binding ability has a wider significance than merely transmission of biological signals. Thus, IGF-II when interacting with receptors realizing its biological effects also binds to the IGF-2 receptor limiting its action and, besides, to the binding proteins (BP) modulating its action. The entire cycle of interactions occurs in the body at different affinity levels. Meanwhile, insulin interacts neither with IGF-2 receptor nor with BP. In this case, specificity and sequence of interaction with each of receptors or with protein are due not to the general design of the peptide structure, but rather to structure of individual submolecular determinants—binding domains. The leading role in disclosure of composition and structure of these domains is played by the “mutant-ligand” approach evaluating affinity of modified analogs. To analyze role of structural elements of the binding domains, the author proposes the system of estimation of affinity of the studied analogs. The present work, alongside with consideration of methodical aspects of the forthcoming analysis, is an introduction to the problem of organization of the binding domains connected directly with functional role of peptides of the insulin type. The proposed analysis is due to necessity of specification of this organization both in one molecule and in different molecules with a similar plan of structure on the basis of not always unanimous literature data and of clarification of principles of structure of these domains.     

Quantification of tertiary structural conservation despite primary sequence drift in the globin fold.

The globin family of protein structures was the first for which it was recognized that tertiary structure can be highly conserved even when primary sequences have diverged to a virtually undetectable level of similarity. This principle of structural inertia in molecular evolution is now evident for many other protein families. We have performed a systematic comparison of the sequences and structures of 6 representative hemoglobin subunits as diverse in origin as plants, clams, and humans. Our analysis is based on a 97-residue helical core in common to all 6 structures. Amino acid sequence identities range from 12.4% to 42.3% in pairwise comparisons, and, despite these variations, the maximal RMS deviation in alpha-carbon positions is 3.02 A. Overall, sequence similarity and structural deviation are significantly anticorrelated, with a correlation coefficient of -0.71, but for a set of structures having under 20% pairwise identity, this anticorrelation falls to -0.38, which emphasizes the weak connection between a specific sequence and the tertiary fold. There is substantial variability in structure outside the helical core, and functional characteristics of these globins also differ appreciably. Nevertheless, despite variations in detail that the sequence dissimilarities and functional differences imply, the core structures of these globins remain remarkably preserved.     

Structural effects of Parkinson's disease linked DJ-1 mutations

Mutations in the protein DJ-1 are associated with familial forms of Parkinson's disease, indicating that DJ-1 may be involved in pathways related to the etiology of this disorder. Here we have used solution state NMR and circular dichroism spectroscopies to evaluate the extent of structural perturbations associated with five different Parkinson's disease linked DJ-1mutations: L166P, E64D, M26I, A104T, and D149A. Comparison of the data with those obtained for the wild-type protein shows that the L166P mutation leads to severe and global destabilization and unfolding of the protein structure, while the structure of the E64D mutation, as expected, is nearly unperturbed. Interestingly, the remaining three mutants all show different degrees of structural perturbation, which are accompanied by a reduction in the thermodynamic stability of the protein. The observed structural and thermodynamic differences are likely to underlie any functional variations between these mutants and the wild type, which in turn are likely responsible for the pathogenicity of these mutations.     

NMR and MD studies on the interaction between ligand peptides and alpha-bungarotoxin

The interaction between alpha-bungarotoxin and linear synthetic peptides, mimotope of the nicotinic acetylcholine receptor binding site, has been characterised extensively by several methods and a wealth of functional, kinetic and structural data are available. Hence, this system represents a suitable model to explore in detail the dynamics of a peptide-protein interaction. Here, the solution structure of a new complex of the protein toxin with a tridecapeptide ligand exhibiting high affinity has been determined by NMR. As observed for three other previously reported mimotope-alpha-bungarotoxin complexes, also in this case correlations between biological activity and kinetic data are not fully consistent with a static discussion of structural data. Molecular dynamics simulations of the four mimotope-toxin complexes indicate that a relevant contribution to the complex stability is given by the extent of the residual flexibility that the protein maintains upon peptide binding. This feature, limiting the entropy loss caused by protein folding and binding, ought to be generally considered in a rational design of specific protein ligands.     

Split Inteins介导的多肽链两端标记

多肽链多位点特异性标记有助于了解蛋白质的结构与功能,特别是在蛋白质的动态构象研究方面. 但是,现有的多肽链多位点特异性标记方法各有局限性,并且种类有限,所以有必要开发新的多肽链多位点特异性标记方法以满足研究需求. 本文以Diub (ubiquitin dimer) 蛋白为研究对象,借助S1和S11型2种不同的断裂蛋白质内含子 (split inteins) 的蛋白质反式剪接,将含有不同荧光基团的两种小肽成功剪接至靶蛋白的两端,最终达到对靶蛋白的末端标记目的.