Solution conformation of endothelin determined by means of 1H-NMR spectroscopy and distance geometry calculations

The structure of endothelin-1 (ET-1), an endothelial cell-derived peptide with vasoconstricting activity, was determined in an aqueous solution by means of a combination of NMR and distance geometry calculations. The resulting structure is characterized by an alpha-helical conformation in the sequence region, Lys9-Cys15. Furthermore, an extended structure and a turn structure exist in the Cys1-Ser4 and Ser5-Asp8 regions respectively, and no preferred conformation was found for the C-terminal part of the peptide which was not uniquely constrained by the NMR data. These structural elements, the alpha-helical structure in the sequence portion, Cys-X-X-X-Cys, and the extended structure in Cys-X-Cys, are homologous to those found commonly in several neurotoxic peptides.     

Efficient parameterized algorithms for biopolymer structure-sequence alignment

Computational alignment of a biopolymer sequence (e.g., an RNA or a protein) to a structure is an effective approach to predict and search for the structure of new sequences. To identify the structure of remote homologs, the structure-sequence alignment has to consider not only sequence similarity, but also spatially conserved conformations caused by residue interactions and, consequently, is computationally intractable. It is difficult to cope with the inefficiency without compromising alignment accuracy, especially for structure search in genomes or large databases. This paper introduces a novel method and a parameterized algorithm for structure-sequence alignment. Both the structure and the sequence are represented as graphs, where, in general, the graph for a biopolymer structure has a naturally small tree width. The algorithm constructs an optimal alignment by finding in the sequence graph the maximum valued subgraph isomorphic to the structure graph. It has the computational time complexity O(k³N²) for the structure of N residues and its tree decomposition of width t. Parameter k, small in nature, is determined by a statistical cutoff for the correspondence between the structure and the sequence. This paper demonstrates a successful application of the algorithm to RNA structure search used for noncoding RNA identification. An application to protein threading is also discussed     

Solution structure of At3g04780.1-des15, an Arabidopsis thaliana ortholog of the C-terminal domain of human thioredoxin-like protein

The structure of At3g04780.1-des15, an Arabidopsis thaliana ortholog of the C-terminal domain of human thioredoxin-like protein, was determined by NMR spectroscopy. The structure is dominated by a beta-barrel sandwich. A two-stranded anti-parallel beta-sheet, which seals off one end of the beta-barrel, is flanked by two flexible loops rich in acidic amino acids. Although this fold often provides a ligand binding site, the structure did not reveal an appreciable cavity inside the beta-barrel. The three-dimensional structure of At3g04780.1-des15 provides an entry point for understanding its functional role and those of its mammalian homologs.     

Neuro-fuzzy structural classification of proteins for improved protein secondary structure prediction

Fourier transform infrared (FTIR) spectroscopy is a very flexible technique for characterization of protein secondary structure. Measurements can be carried out rapidly in a number of different environments based on only small quantities of proteins. For this technique to become more widely used for protein secondary structure characterization, however, further developments in methods to accurately quantify protein secondary structure are necessary. Here we propose a structural classification of proteins (SCOP) class specialized neural networks architecture combining an adaptive neuro-fuzzy inference system (ANFIS) with SCOP class specialized backpropagation neural networks for improved protein secondary structure prediction. Our study shows that proteins can be accurately classified into two main classes "all alpha proteins" and "all beta proteins" merely based on the amide I band maximum position of their FTIR spectra. ANFIS is employed to perform the classification task to demonstrate the potential of this architecture with moderately complex problems. Based on studies using a reference set of 17 proteins and an evaluation set of 4 proteins, improved predictions were achieved compared to a conventional neural network approach, where structure specialized neural networks are trained based on protein spectra of both "all alpha" and "all beta" proteins. The standard errors of prediction (SEPs) in % structure were improved by 4.05% for helix structure, by 5.91% for sheet structure, by 2.68% for turn structure, and by 2.15% for bend structure. For other structure, an increase of SEP by 2.43% was observed. Those results were confirmed by a "leave-one-out" run with the combined set of 21 FTIR spectra of proteins.     

Unusual mRNA pseudoknot structure is recognized by a protein translational repressor

Translation of ribosomal proteins in the alpha operon of E. coli is repressed by one of the encoded proteins, S4; it specifically recognizes an RNA fragment containing the translational initiation site for the first gene in the operon. RNA structure mapping experiments have suggested a pseudoknot structure for the S4 binding site: the loop of a hairpin is base paired to sequences downstream of the hairpin. Here, we systematically test this proposed structure by measuring S4 binding to an extensive set of site-directed mutations that create compensatory base pair changes in potential helices. The pseudoknot folding is confirmed, and two additional, unexpected interactions within the pseudoknot are also detected. The overall structure is an unusual "double pseudoknot" linking a hairpin upstream of the ribosome binding site with sequences 2-10 codons downstream of the initiation codon. Stabilization of this structure by S4 could account for translational repression.     

Structure-activity relationships of novel anti-malarial agents. Part 2: cinnamic acid derivatives

We have described compound 1 as a lead structure for a novel class of anti-malarial agents. Replacement of the 3-phenylpropionyl moiety of the lead structure 1 by a 4-propoxycinnamic acid residue resulted in a significant improvement in antimalarial activity. Compound 3q represents an important step in the development of lead structure 1 into an anti-malarial drug candidate.     

Pressure-induced changes in the solution structure of the GB1 domain of protein G

The solution structure of the GB1 domain of protein G at a pressure of 2 kbar is presented. The structure was calculated as a change from an energy-minimised low-pressure structure using (1)H chemical shifts. Two separate changes can be characterised: a compression/distortion, which is linear with pressure; and a stabilisation of an alternative folded state. On application of pressure, linear chemical shift changes reveal that the backbone structure changes by about 0.2 A root mean square, and is compressed by about 1% overall. The alpha-helix compresses, particularly at the C-terminal end, and moves toward the beta-sheet, while the beta-sheet is twisted, with the corners closest to the alpha-helix curling up towards it. The largest changes in structure are along the second beta-strand, which becomes more twisted. This strand is where the protein binds to IgG. Curved chemical shift changes with pressure indicate that high pressure also populates an alternative structure with a distortion towards the C-terminal end of the helix, which is likely to be caused by insertion of a water molecule.     

Structure stability analysis of industrial system and construction of eco-industrial chain: A case study of Lianyungang Xuwei New Area,China

Structure stability analysis is a vital prerequisite for the construction of eco-industrial chain. This study proposed four steps for applying an integrated structure stability index, which was derived to characterize both the diversity and equilibrium of eco-industrial chain, to structure stability analysis, and searching a final eco-industrial chain with the highest structure stability as follows: (1) analyzing links and link points of existing and planned industries through material and energy flow analysis; (2) identifying supplementary industries, which are integrated with pillar industries to enhance structure stability, through forward diffusion effect, backward diffusion effect and sideward diffusion effect analysis of pillar industries; (3) testing industrial structure stability by introducing an integrated structure stability index; (4) coupling of all definitized links and link points.These four steps were applied to a case study of Lianyungang Xuwei New Area, China. The results showed that: (1) there were 9 link points and 17 links of the planned five pillar industries in this area; (2) cement, building materials, chemical fertilizer, fine chemicals and marine chemicals with newly added 23 link points and 44 links should be considered as the supplementary industries; (3) the improvement in structure stability was confirmed by the integrated structure stability index increasing from 0.116 to 0.158 after adding supplementary industries; and (4) the final “steel–petrochemical-equipment manufacturing-logistics-IGCC polygeneration” eco-industrial chain was constructed by coupling 32 link points and 61 links. Therefore, by taking pillar industries as the kernel and by introducing an integrated structure stability index for the verification of stability improvement in industrial structure, these proposed four steps were proved to be an effective process and feasible method for enhancing the stability of eco-industrial system structure.     

Structure of the Bacillus subtilis quorum-sensing peptide pheromone ComX

The ComX pheromone is an extracellular signaling molecule that stimulates natural competence in response to crowding in the gram-positive bacterium Bacillus subtilis. The pheromone is formed by isoprenylation of an inactive precursor peptide, but its precise structure is not known. Here we report the structure of the ComX pheromone, showing that addition of a geranyl group to a tryptophan residue results in the formation of an unusual ring structure.     

Low Energy Conformations for S100 Binding Peptide from the Negative Regulatory Domain of p53

We have computed the low energy conformations of the negative regulatory domain of p53, residues 374–388 using Empirical Conformational Energies of Peptides Program including solvation and computed the statistical weights of distinct conformational states. We find that there are two high probability conformations, one an α-helix from Lys 374–Lys 381, followed by another helical structure involving Lys 382–Glu 388 (statistical weight of 0.48) and an all-α-helix for the entire sequence (statistical weight of 0.23). Both structure are superimposable on the NMR structure of this sequence bound to the S100 protein. The global minimum structure (statistical weight of 0.014) is a beta structure from Gly 374–Arg 379 followed by an α-helix from His 380–Glu 388. Based on these results, we propose a possible strategy for enhancement of p53 anti-tumor activity in cancer cells. Since the structure of this sequence bound to the sirtuin protein, Sir2, is a β-sheet, we further propose that the global minimum may be an intermediate on the α-β structure transition.     

Boltzmann probability of RNA structural neighbors and riboswitch detection

MOTIVATION: We describe algorithms implemented in a new software package, RNAbor, to investigate structures in a neighborhood of an input secondary structure S of an RNA sequence s. The input structure could be the minimum free energy structure, the secondary structure obtained by analysis of the X-ray structure or by comparative sequence analysis, or an arbitrary intermediate structure. RESULTS: A secondary structure T of s is called a delta-neighbor of S if T and S differ by exactly delta base pairs. RNAbor computes the number (N(delta)), the Boltzmann partition function (Z(delta)) and the minimum free energy (MFE(delta)) and corresponding structure over the collection of all delta-neighbors of S. This computation is done simultaneously for all delta < or = m, in run time O (mn3) and memory O(mn2), where n is the sequence length. We apply RNAbor for the detection of possible RNA conformational switches, and compare RNAbor with the switch detection method paRNAss. We also provide examples of how RNAbor can at times improve the accuracy of secondary structure prediction. AVAILABILITY: http://bioinformatics.bc.edu/clotelab/RNAbor/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Rethinking the importance of the structure of ecological networks under an environment‐dependent framework

A major quest in network and community ecology has been centered on understanding the importance of structural patterns in species interaction networks—the synthesis of who interacts with whom in a given location and time. In the past decades, much effort has been devoted to infer the importance of a particular structure by its capacity to tolerate an external perturbation on its structure or dynamics. Here, we demonstrate that such a perspective leads to inconsistent conclusions. That is, the importance of a network structure changes as a function of the external perturbations acting on a community at any given point in time. Thus, we discuss a research agenda to investigate the relative importance of the structure of ecological networks under an environment‐dependent framework. We hypothesize that only by studying systematically the link between network structure and community dynamics under an environment‐dependent framework, we can uncover the limits at which communities can tolerate environmental changes.     

Allosteric inhibition via R-state destabilization in ATP sulfurylase from Penicillium chrysogenum

The structure of the cooperative hexameric enzyme ATP sulfurylase from Penicillium chrysogenum bound to its allosteric inhibitor, 3'-phosphoadenosine-5'-phosphosulfate (PAPS), was determined to 2.6 A resolution. This structure represents the low substrate-affinity T-state conformation of the enzyme. Comparison with the high substrate-affinity R-state structure reveals that a large rotational rearrangement of domains occurs as a result of the R-to-T transition. The rearrangement is accompanied by the 17 A movement of a 10-residue loop out of the active site region, resulting in an open, product release-like structure of the catalytic domain. Binding of PAPS is proposed to induce the allosteric transition by destabilizing an R-state-specific salt linkage between Asp 111 in an N-terminal domain of one subunit and Arg 515 in the allosteric domain of a trans-triad subunit. Disrupting this salt linkage by site-directed mutagenesis induces cooperative inhibition behavior in the absence of an allosteric effector, confirming the role of these two residues.     

The quaternary structure of carbonmonoxy hemoglobin ypsilanti

We present a geometric analysis of the allosteric interface in the new Y state quaternary structure observed in liganded mutant hemoglobin Ypsilanti (β99 Asp → Tyr) by Smith, F.R., Lattman, E.E., Carter, C.W., Jr. (Proteins 10:81–91, 1991). The classical T to R quaternary structure change being a rotation of αβ dimers about an axis which is approximately parallel to the dimer axis of pseudosym-metry, the new quaternary structure is obtained by applying to R an additional rotation about an axis orthogonal to the first. This suggests that Y is a modified R state rather than an intermediate on the T to R pathway. Computer docking experiments designed to simulate the quaternary structure change support this suggestion. © 1993 Wiley-Liss, Inc.     

Rapid formation of secondary structure framework in protein folding studied by stopped-flow circular dichroism

Kinetic refolding reactions of ferricytochrome c and beta-lactoglobulin have been studied by stopped-flow circular dichroism by monitoring rapid ellipticity changes of peptide backbone and side-chain chromophores. In both proteins, a transient intermediate accumulates within the dead time of stopped-flow mixing (18 ms), and the intermediate has an appreciable amount of secondary structure but possesses an unfolded tertiary structure. It is suggested that the rapid formation of a secondary structure framework in protein folding is a common property observed in a variety of globular proteins.     

Potassium ions modulate a G-quadruplex-ribozyme's activity

Hepatitis delta virus ribozyme folds into a tightly packed tertiary structure. However, unlike other ribozymes, it does not appear to be able to follow alternative folding pathways. Molecular engineering of the hepatitis delta virus ribozyme led to the development of a ribozyme possessing an endoribonuclease activity that is under the control of a G-quadruplex structure (i.e., a G-quartzyme). This latter species represents an entirely new class of ribozyme. Mutants of this ribozyme were then generated in order to shed light on the modulation of the cleavage activity caused by the presence of the G-quadruplex structure. Kinetic characterization of the G-quartzyme was performed under various single turnover conditions. It was found to be active only in the presence of potassium cations that act as counter ions in the positioning of the four coplanar guanines that form the building block of the G-quadruplex structure. The G-quartzyme behaves as an allosteric ribozyme, with the potassium cations acting as positive effectors with a Hill coefficient of 2.9 +/- 0.2. The conformation transition caused by the presence of the potassium ions is supported by enzymatic and chemical probing of both the inactive (off) and active (on) structures. This study shows that it is possible to interfere with the tight structure of the hepatitis delta virus ribozyme by adding an unusual, stable structure. To our knowledge, the G-quartzyme is the sole ribozyme that exhibits a monovalent cation-dependent activity.     

An information measure of the quality of protein secondary structure prediction.

We describe an information-theory-based measure of the quality of secondary structure prediction (RELINFO). RELINFO has a simple yet intuitive interpretation: it represents the factor by which secondary structure choice at a residue has been restricted by a prediction scheme. As an alternative interpretation of secondary structure prediction, RELINFO complements currently used methods by providing an information-based view as to why a prediction succeeds and fails. To demonstrate this score's capabilities, we applied RELINFO to an analysis of a large set of secondary structure predictions obtained from the first five rounds of the Critical Assessment of Structure Prediction (CASP) experiment. RELINFO is compared with two other common measures: percent correct (Q3) and secondary structure overlap (SOV). While the correlation between Q3 and RELINFO is approximately 0.85, RELINFO avoids certain disadvantages of Q3, including overestimating the quality of a prediction. The correlation between SOV and RELINFO is approximately 0.75. The valuable SOV measure unfortunately suffers from a saturation problem, and perhaps has unfairly given the general impression that secondary structure prediction has reached its limit since SOV hasn't improved much over the recent rounds of CASP. Although not a replacement for SOV, RELINFO has greater dispersion. Over the five rounds of CASP assessed here, RELINFO shows that predictions targets have been more difficult in successive CASP experiments, yet the predictions quality has continued to improve measurably over each round. In terms of information, the secondary structure prediction quality has almost doubled from CASP1 to CASP5. Therefore, as a different perspective of accuracy, RELINFO can help to improve prediction of protein secondary structure by providing a measure of difficulty as well as final quality of a prediction.     

Assessment of protein model structure accuracy estimation in CASP13: Challenges in the era of deep learning

Scoring model structure is an essential component of protein structure prediction that can affect the prediction accuracy tremendously. Users of protein structure prediction results also need to score models to select the best models for their application studies. In Critical Assessment of techniques for protein Structure Prediction (CASP), model accuracy estimation methods have been tested in a blind fashion by providing models submitted by the tertiary structure prediction servers for scoring. In CASP13, model accuracy estimation results were evaluated in terms of both global and local structure accuracy. Global structure accuracy estimation was evaluated by the quality of the models selected by the global structure scores and by the absolute estimates of the global scores. Residue-wise, local structure accuracy estimations were evaluated by three different measures. A new measure introduced in CASP13 evaluates the ability to predict inaccurately modeled regions that may be improved by refinement. An intensive comparative analysis on CASP13 and the previous CASPs revealed that the tertiary structure models generated by the CASP13 servers show very distinct features. Higher consensus toward models of higher global accuracy appeared even for free modeling targets, and many models of high global accuracy were not well optimized at the atomic level. This is related to the new technology in CASP13, deep learning for tertiary contact prediction. The tertiary model structures generated by deep learning pose a new challenge for EMA (estimation of model accuracy) method developers. Model accuracy estimation itself is also an area where deep learning can potentially have an impact, although current EMA methods have not fully explored that direction.     

Descent of a split RNA

tmRNA is found in the usual one-piece form in most cyanobacteria, but a small clade has an unusual two-piece form, resulting from processing of a circularly permuted precursor RNA. Here, the secondary structure of the cyanobacterial one-piece tmRNA is established by phylogenetic sequence comparison; it deviates from most tmRNAs in having an extra (fifth) pseudoknot and a branched structure at the end of the tag reading frame. Patterns of sequence conservation between the two cyanobacterial tmRNA types suggest particular events in the evolution of the two-piece form. A simple gene permutation model of tandem duplication followed by loss of the outermost segments appears inadequate. An additional rearrangement is proposed, in which a structure impaired by deletion was regenerated at the expense of a neighboring structure.     

Unusual Armadillo Fold in the Human General Vesicular Transport Factor p115

The golgin family gives identity and structure to the Golgi apparatus and is part of a complex protein network at the Golgi membrane. The golgin p115 is targeted by the GTPase Rab1a, contains a large globular head region and a long region of coiled-coil which forms an extended rod-like structure. p115 serves as vesicle tethering factor and plays an important role at different steps of vesicular transport. Here we present the 2.2 Å-resolution X-ray structure of the globular head region of p115. The structure exhibits an armadillo fold that is decorated by elongated loops and carries a C-terminal non-canonical repeat. This terminal repeat folds into the armadillo superhelical groove and allows homodimeric association with important implications for p115 mediated multiple protein interactions and tethering.