Modeling and structural analysis of evolutionarily diverse S8 family serine proteases

Serine proteases are an abundant class of enzymes that are involved in a wide range of physiological processes and are classified into clans sharing structural homology. The active site of the subtilisin-like clan contains a catalytic triad in the order Asp, His, Ser (S8 family) or a catalytic tetrad in the order Glu, Asp and Ser (S53 family). The core structure and active site geometry of these proteases is of interest for many applications. The aim of this study was to investigate the structural properties of different S8 family serine proteases from a diverse range of taxa using molecular modeling techniques. In conjunction with 12 experimentally determined three-dimensional structures of S8 family members, our predicted structures from an archaeon, protozoan and a plant were used for analysis of the catalytic core. Amino acid sequences were obtained from the MEROPS database and submitted to the LOOPP server for threading based structure prediction. The predicted structures were refined and validated using PROCHECK, SCRWL and MODELYN. Investigation of secondary structures and electrostatic surface potential was performed using MOLMOL. Encompassing a wide range of taxa, our structural analysis provides an evolutionary perspective on S8 family serine proteases. Focusing on the common core containing the catalytic site of the enzyme, the analysis presented here is beneficial for future molecular modeling strategies and structure-based rational drug design.     

Prediction of secondary structures of 16S and 23S rRNA fromEscherichia coli

Small and large subunits ofEscherichia coli ribosome have three different rRNAs, the sequences of which are known. However, attempts by three groups to predict secondary structures of 16S and 23S rRNAs have certain common limitations namely, these structures are predicted assuming no interactions among various domains of the molecule and only 40% residues are involved in base pairing as against the experimental observation of 60 % residues in base paired state. Recent experimental studies have shown that there is a specific interaction between naked 16S and 23S rRNA molecules. This is significant because we have observed that the regions (oligonucleotides of length 9–10 residues), in 16S rRNA which are complementary to those in 23S rRNA do not have internal complementary sequences. Therefore, we have developed a simple graph theoretical approach to predict secondary structures of 16S and 23S rRNAs. Our method for model building not only uses complete sequence of 16S or 23S rRNA molecule along with other experimental observations but also takes into account the observation that specific recognition is possible through the complementary sequences between 16S and 23S rRNA molecules and, therefore, these parts of the molecules are not used for internal base pairing. The method used to predict secondary structures is discussed. A typical secondary structure of the complex between 16S and 23S rRNA molecules, obtained using our method, is presented and compared Briefly with earlier model Building studies.     

Crystal structures of S100A6 in the Ca(2+)-free and Ca(2+)-bound states: the calcium sensor mechanism of S100 proteins revealed at atomic resolution

S100A6 is a member of the S100 family of Ca(2+) binding proteins, which have come to play an important role in the diagnosis of cancer due to their overexpression in various tumor cells. We have determined the crystal structures of human S100A6 in the Ca(2+)-free and Ca(2+)-bound states to resolutions of 1.15 A and 1.44 A, respectively. Ca(2+) binding is responsible for a dramatic change in the global shape and charge distribution of the S100A6 dimer, leading to the exposure of two symmetrically positioned target binding sites. The results are consistent with S100A6, and most likely other S100 proteins, functioning as Ca(2+) sensors in a way analogous to the prototypical sensors calmodulin and troponin C. The structures have important implications for our understanding of target binding and cooperativity of Ca(2+) binding in the S100 family.     

Modeling of peptides connecting the ligand‐binding and transmembrane domains in the GluR2 glutamate receptor

Ligand‐gated Glutamate receptors (GluR) mediate synaptic signals in the nervous system. Ionotropic GluRs of AMPA type, the subject of this study, are tetrameric assemblies of monomer subunits, each of which is constructed in a modular fashion from functional subdomains. The extracellular ligand‐binding domain (LBD) changes its conformation upon binding of an agonist ligand followed by opening of a transmembrane (TM) ion channel. Peptides connecting the LBD and TM domains facilitate gating of the channel, and their structure and composition are important for the receptor functioning. In this study, we used replica exchange molecular dynamics (REMD) simulations to model S1M1 and S2M3 connecting peptides of the GluR2 receptor in two implicit solvents, water and interfacial water/lipid medium characterized by lower polarity. Propensity of these peptides to form helical structures was analyzed using helicity measure derived from the free energy of the simulated ensembles of structures. The S1M1 and S2M3 connecting peptides were not helical in our simulations in both dielectric environments in the absence of the rest of the protein. The structures of the LBD fragment with known high‐resolution α‐helical structure and of the TM3 helix were successfully predicted in the simulations, which in part validate our results. The S2M3 peptide, which is important in gating, formed a well‐defined coil structure and salt‐bridges with the S2 domain. The S1M1 peptide formed a loop structure via formation of internal salt‐bridges. Potential implications of these structures on function of the receptor are discussed. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Atomic structures at last: the ribosome in 2000

Last year, atomic structures of the 50S ribosomal subunit from Haloarcula marismortui and of the 30S ribosomal subunit from Thermus thermophilus were published. A year before that, a 7.8 A resolution electron density map of the 70S ribosome from T. thermophilus appeared. This information is revolutionizing our understanding of protein synthesis.     

SARS病毒S蛋白三维结构预测

蛋白质结构类型识别方法可以在没有序列同源性的蛋白质之间检测有没有结构相似性。利用蛋白质结构类型识别方法预测了SARS病毒S蛋白N端区域的结构。模建的SARS病毒S蛋白N端区域是一个全折叠的结构。     

Expression of S100A2 and S100P in human eccrine sweat glands and their application in differentiating secretory coil-like from duct-like structures in the 3D reconstituted eccrine sweat spheroids

Secretory coils and ducts are two components of eccrine sweat glands with different structures and functions. In our previous study, we combined keratins and α-SMA to distinguish between secretory coils and ducts. However, the key deficiency of the method was that none of the antibodies used was specific for ducts. In this study, we first examined the co-localization of K5/K7, α-SMA/K14, K7/S100P and α-SMA/S100A2 by double-immunofluorescence staining to confirm the localization of S100P and S100A2 in native human eccrine sweat glands, and second we identified secretory coil-like and duct-like structures in the 3D reconstituted eccrine sweat gland spheroids by double-immunofluorescence staining for K7/S100P and α-SMA/S100A2. In native human eccrine sweat glands, S100A2 immunoreactivity was confined to the outer layer and S100P to the inner layer of the duct. In 12-week Matrigel plugs containing eccrine sweat gland cells, double-immunofluorescence staining for K7/S100P and α-SMA/S100A2 could easily distinguish duct-like structures from secretory coil-like structures. We conclude that S100A2 and S100P can be used as specific duct markers in eccrine sweat glands, and combined use of S100P or S100A2 with keratins enables easy to distinction between secretory coils and ducts.     

Evaluation of the denaturing gradient gel electrophoresis-apparatus as a parameter influencing soil microbial community fingerprinting

We compared two denaturing gradient gel electrophoresis (DGGE) systems—DCode (Biorad, Hercules, CA, USA) and PhorU (Ingeny, Leiden, NL), performing community level 16S and 18S rRNA gene fragment-PCR-DGGE with total DNA extracted from upland pasture soil used for outdoor cattle husbandry. The methodological evaluation of the DGGE apparatus as parameter influencing DGGE fingerprinting, based on cluster analysis of soil bacterial and fungal community fingerprints, was made in terms of the resulting information about microbial community structures and their response to different degrees of cattle impact. Although the comparative DGGE analysis with different DGGE systems provided similar clustering of microbial community structures in correlation with the degree of cattle impact, our results suggest the DGGE system to be a factor influencing DGGE analysis. To our knowledge this is the first attempt to investigate the hypothetical impact of the DGGE system due to different technical characteristics, recommending the use of one and the same DGGE apparatus throughout an experiment, if the monitoring of microbial community structures requires multiple gel-to-gel analysis.     

An evaluation of ensign wasp classification (Hymenoptera: Evaniidae) based on molecular data and insights from ribosomal RNA secondary structure

Abstract.  Ensign wasps (Hymenoptera: Evaniidae) are colourful, frequently collected and easily distinguished from other parasitic Hymenoptera. Despite many fascinating biological attributes, this group of insects has been overlooked by ecologists and systematists. An imposing obstacle inhibiting research on these wasps is the current state of their chaotic and potentially flawed classification, which has more than 50% of all described species assigned to the genus Evania – a taxon long suspected of being polyphyletic. The generic classification has recently been redefined on the basis of morphological characters. We tested this reinterpreted classification by analysing sequence data from three genes [28S ribosomal RNA (rRNA), 16S rRNA and cytochrome oxidase I (COI)] under parsimony and Bayesian criteria. For the 28S and 16S rRNAs, we illustrate the predicted secondary structures and provide a series of summary statistics for them; information pertaining to these structures was incorporated into our phylogenetic analyses where appropriate. Phylogenetically, our results indicate that this new generic classification is relatively sound, but that more data are required to understand intergeneric relationships.     

Tagging ribosomal protein S7 allows rapid identification of mutants defective in assembly and function of 30 S subunits

Ribosomal protein S7 nucleates folding of the 16 S rRNA 3' major domain, which ultimately forms the head of the 30 S ribosomal subunit. Recent crystal structures indicate that S7 lies on the interface side of the 30 S subunit, near the tRNA binding sites of the ribosome. To map the functional surface of S7, we have tagged the protein with a Protein Kinase A recognition site and engineered alanine substitutions that target each exposed, conserved residue. We have also deleted conserved features of S7, using its structure to guide our design. By radiolabeling the tag sequence using Protein Kinase A, we are able to track the partitioning of each mutant protein into 30 S, 70 S, and polyribosome fractions in vivo. Overexpression of S7 confers a growth defect, and we observe a striking correlation between this phenotype and proficiency in 30 S subunit assembly among our collection of mutants. We find that the side chain of K35 is required for efficient assembly of S7 into 30 S subunits in vivo, whereas those of at least 17 other conserved exposed residues are not required. In addition, an S7 derivative lacking the N-terminal 17 residues causes ribosomes to accumulate on mRNA to abnormally high levels, indicating that our approach can yield interesting mutant ribosomes.     

Determination of base pairing in Escherichia coli and Bacillus stearothermophilus 5S RNAs by infrared spectroscopy.

The extent of base pairing in Escherichia coli and Bacillus stearothermophilus 5S RNAs was determined by infrared spectroscopy. From the infrared spectra taken at 20 degrees and 52 degrees C it is concluded that E. coli and B. stearothermophlius 5S RNAs possess a large number of base pairs (Table I). Comparison of our results with those previously published using other methods leads to the conclusion that the structures of prokaryotic 5S RNAs involve a large number of tertiary interactions, in which the base pairing is not necessarily solely of the Watson-Crick type.     

Improved model of a LexA repressor dimer bound to recA operator

A complete three dimensional model for the LexA repressor dimer bound to the recA operator site consistent with relevant biochemical and biophysical data for the repressor was proposed from our laboratory when no crystal structure of LexA was available. Subsequently, the crystal structures of four LexA mutants Delta(1-67) S119A, S119A, G85D and Delta(1-67) quadruple mutant in the absence of operator were reported. It is examined in this paper to what extent our previous model was correct and how, using the crystal structure of the operator-free LexA dimer we can predict an improved model of LexA dimer bound to recA operator. In our improved model, the C-domain dimerization observed repeatedly in the mutant operator-free crystals is retained but the relative orientation between the two domains within a LexA molecule changes. The crystal structure of wild type LexA with or without the recA operator cannot be solved as it autocleaves itself. We argue that the 'cleavable' cleavage site region found in the crystal structures is actually the more relevant form of the region in wild-type LexA since it agrees with the value of the pre-exponential Arrhenius factor for its autocleavage, absence of various types of trans-cleavages, difficulty in modifying the catalytic serine by diisopropyl flourophosphate and lack of cleavage at Arg 81 by trypsin; hence the concept of a 'conformational switch' inferred from the crystal structures is meaningless.     

The evolving genome of Salmonella enterica serovar Pullorum

Salmonella enterica serovar Pullorum is a fowl-adapted bacterial pathogen that causes dysentery (pullorum disease). Host adaptation and special pathogenesis make S. enterica serovar Pullorum an exceptionally good system for studies of bacterial evolution and speciation, especially regarding pathogen-host interactions and the acquisition of pathogenicity. We constructed a genome map of S. enterica serovar Pullorum RKS5078, using I-CeuI, XbaI, AvrII, and SpeI and Tn10 insertions. Pulsed-field gel electrophoresis was employed to separate the large DNA fragments generated by the endonucleases. The genome is 4,930 kb, which is similar to most salmonellas. However, the genome of S. enterica serovar Pullorum RKS5078 is organized very differently from the majority of salmonellas, with three major inversions and one translocation. This extraordinary genome structure was seen in most S. enterica serovar Pullorum strains examined, with different structures in a minority of S. enterica serovar Pullorum strains. We describe the coexistence of different genome structures among the same bacteria as genomic plasticity. Through comparisons with S. enterica serovar Typhimurium, we resolved seven putative insertions and eight deletions ranging in size from 12 to 157 kb. The genomic plasticity seen among S. enterica serovar Pullorum strains supported our hypothesis about its association with bacterial evolution: a large genomic insertion (157 kb in this case) disrupted the genomic balance, and rebalancing by independent recombination events in individual lineages resulted in diverse genome structures. As far as the structural plasticity exists, the S. enterica serovar Pullorum genome will continue evolving to reach a further streamlined and balanced structure.     

Automated motif extraction and classification in RNA tertiary structures

We used a novel graph-based approach to extract RNA tertiary motifs. We cataloged them all and clustered them using an innovative graph similarity measure. We applied our method to three widely studied structures: Haloarcula marismortui 50S (H.m 50S), Escherichia coli 50S (E. coli 50S), and Thermus thermophilus 16S (T.th 16S) RNAs. We identified 10 known motifs without any prior knowledge of their shapes or positions. We additionally identified four putative new motifs.     

一种预测单链RNA分子二级结构的计算机算法

本文给出了一个利用已知能量数据构成具有最小自由能的单链RNA分子二级结构的计算机算法,并给出了此算法的可行性证明和应用实例。     

Structural studies on site-directed mutants of domain 3 of Xenopus laevis oocyte 5 S ribosomal RNA

Base substitutions have been introduced into the highly conserved sequences of loops D and E within domain 3 of Xenopus laevis oocyte 5 S rRNA. The effects of these mutations on the solution structure of this 5 S rRNA have been studied by means of probing with nucleases, and with chemical reagents under native and semi-denaturing conditions. The data obtained with these mutants support the graphic model of Xenopus oocyte 5 S rRNA proposed by Westhof et al. In particular, our results rule out the existence of long-range base-pairing interactions between loop C and either loop D or loop E. The data also confirm that loops D and E in the wild-type 5 S RNA adopt unusual secondary structures and illustrate the importance of nucleotide sequence in the formation of intrinsic local loop conformations via non-canonical base-pairs and specific base-phosphate contacts. Consistent with this conclusion is our observation that the domain 3 fragment of Xenopus oocyte 5 S rRNA adopts the same conformation as the corresponding region in the full-length 5 S rRNA.     

Spina cortica and Tapetum spinosus, two new microstructures of flight feathers: description, function and distribution in modern birds

The importance of feathers for the avian group has made them one of the most studied epidermal structures both from the morphological and evolutionary point of view. Surprisingly, our observations by Scanning Electron Microscopy detected the presence of two structures widely distributed within different avian groups and not yet described. In this paper we describe these two new structures (Spina cortica and Tapetum spinosus) and map their distribution within modern birds. The S. cortica is a thorn-like microstructure that grows on the barb cortex and the T. spinosus is the assemblage of these thorns. The distribution of these new structures among birds and their morphological diversity could be of great interest to taxonomists and evolutionary biologists interested in the origin of bird flight.     

Models of voltage-dependent conformational changes in NaChBac channels

Models of the transmembrane region of the NaChBac channel were developed in two open/inactivated and several closed conformations. Homology models of NaChBac were developed using crystal structures of Kv1.2 and a Kv1.2/2.1 chimera as templates for open conformations, and MlotiK and KcsA channels as templates for closed conformations. Multiple molecular-dynamic simulations were performed to refine and evaluate these models. A striking difference between the S4 structures of the Kv1.2-like open models and MlotiK-like closed models is the secondary structure. In the open model, the first part of S4 forms an α-helix, and the last part forms a 3₁₀ helix, whereas in the closed model, the first part of S4 forms a 3₁₀ helix, and the last part forms an α-helix. A conformational change that involves this type of transition in secondary structure should be voltage-dependent. However, this transition alone is not sufficient to account for the large gating charge movement reported for NaChBac channels and for experimental results in other voltage-gated channels. To increase the magnitude of the motion of S4, we developed another model of an open/inactivated conformation, in which S4 is displaced farther outward, and a number of closed models in which S4 is displaced farther inward. A helical screw motion for the α-helical part of S4 and a simple axial translation for the 3₁₀ portion were used to develop models of these additional conformations. In our models, four positively charged residues of S4 moved outwardly during activation, across a transition barrier formed by highly conserved hydrophobic residues on S1, S2, and S3. The S4 movement was coupled to an opening of the activation gate formed by S6 through interactions with the segment linking S4 to S5. Consistencies of our models with experimental studies of NaChBac and K_v channels are discussed.     

RNA structure comparison,motif search and discovery using a reduced representation of RNA conformational space

Given the wealth of new RNA structures and the growing list of RNA functions in biology, it is of great interest to understand the repertoire of RNA folding motifs. The ability to identify new and known motifs within novel RNA structures, to compare tertiary structures with one another and to quantify the characteristics of a given RNA motif are major goals in the field of RNA research; however, there are few systematic ways to address these issues. Using a novel approach for visualizing and mathematically describing macromolecular structures, we have developed a means to quantitatively describe RNA molecules in order to rapidly analyze, compare and explore their features. This approach builds on the alternative eta,theta convention for describing RNA torsion angles and is executed using a new program called PRIMOS. Applying this methodology, we have successfully identified major regions of conformational change in the 50S and 30S ribosomal subunits, we have developed a means to search the database of RNA structures for the prevalence of known motifs and we have classified and identified new motifs. These applications illustrate the powerful capabilities of our new RNA structural convention, and they suggest future adaptations with important implications for bioinformatics and structural genomics.     

Structural organization of the 19S proteasome lid: insights from MS of intact complexes

The 26S proteasome contains a 19S regulatory particle that selects and unfolds ubiquitinated substrates for degradation in the 20S catalytic particle. To date there are no high-resolution structures of the 19S assembly, nor of the lid or base subcomplexes that constitute the 19S. Mass spectra of the intact lid complex from Saccharomyces cerevisiae show that eight of the nine subunits are present stoichiometrically and that a stable tetrameric subcomplex forms in solution. Application of tandem mass spectrometry to the intact lid complex reveals the subunit architecture, while the coupling of a cross-linking approach identifies further interaction partners. Taking together our results with previous analyses we are able to construct a comprehensive interaction map. In summary, our findings allow us to identify a scaffold for the assembly of the particle and to propose a regulatory mechanism that prevents exposure of the active site until assembly is complete. More generally, the results highlight the potential of mass spectrometry to add crucial insight into the structural organization of an endogenous, wild-type complex.