CLICK--topology-independent comparison of biomolecular 3D structures

SupeRNAlign: a new tool for flexible superposition of homologous RNA structures and inference of accurate structure-based sequence alignments

RNA has been found to play an ever-increasing role in a variety of biological processes. The function of most non-coding RNA molecules depends on their structure. Comparing and classifying macromolecular 3D structures is of crucial importance for structure-based function inference and it is used in the characterization of functional motifs and in structure prediction by comparative modeling. However, compared to the numerous methods for protein structure superposition, there are few tools dedicated to the superimposing of RNA 3D structures. Here, we present SupeRNAlign (v1.3.1), a new method for flexible superposition of RNA 3D structures, and SupeRNAlign-Coffee—a workflow that combines SupeRNAlign with T-Coffee for inferring structure-based sequence alignments. The methods have been benchmarked with eight other methods for RNA structural superposition and alignment. The benchmark included 151 structures from 32 RNA families (with a total of 1734 pairwise superpositions). The accuracy of superpositions was assessed by comparing structure-based sequence alignments to the reference alignments from the Rfam database. SupeRNAlign and SupeRNAlign-Coffee achieved significantly higher scores than most of the benchmarked methods: SupeRNAlign generated the most accurate sequence alignments among the structure superposition methods, and SupeRNAlign-Coffee performed best among the sequence alignment methods.

In silico methods for design of biological therapeutics

It has been twenty years since the first rationally designed small molecule drug was introduced into the market. Since then, we have progressed from designing small molecules to designing biotherapeutics. This class of therapeutics includes designed proteins, peptides and nucleic acids that could more effectively combat drug resistance and even act in cases where the disease is caused because of a molecular deficiency. Computational methods are crucial in this design exercise and this review discusses the various elements of designing biotherapeutic proteins and peptides. Many of the techniques discussed here, such as the deterministic and stochastic design methods, are generally used in protein design. We have devoted special attention to the design of antibodies and vaccines. In addition to the methods for designing these molecules, we have included a comprehensive list of all biotherapeutics approved for clinical use. Also included is an overview of methods that predict the binding affinity, cell penetration ability, half-life, solubility, immunogenicity and toxicity of the designed therapeutics. Biotherapeutics are only going to grow in clinical importance and are set to herald a new generation of disease management and cure.

Insights into a key sulfite scavenger enzyme sulfite oxidase (SOX) gene in plants

Sulfite oxidase (SOX) is a crucial molybdenum cofactor-containing enzyme in plants that re-oxidizes the sulfite back to sulfate in sulfite assimilation pathway. However, studies of this crucial enzyme are quite limited hence this work was attempted to understand the SOXs in four plant species namely, Arabidopsis thaliana, Solanum lycopersicum, Populus trichocarpa and Brachypodium distachyon. Herein studied SOX enzyme was characterized with both oxidoreductase molybdopterin binding and Mo-co oxidoreductase dimerization domains. The alignment and motif analyses revealed the highly conserved primary structure of SOXs. The phylogeny constructed with additional species demonstrated a clear divergence of monocots, dicots and lower plants. In addition, to further understand the phylogenetic relationship and make a functional inference, a structure-based phylogeny was constructed using normalized RMSD values in five superposed models from four modelled plant SOXs herein and one previously characterized chicken SOX structure. The plant and animal SOXs showed a clear divergence and also implicated their functional divergences. Based on tree topology, monocot B. distachyon appeared to be diverged from other dicots, pointing out a possible monocot-dicot split. The expression patterns of sulfite scavengers including SOX were differentially modulated under cold, heat, salt and high light stresses. Particularly, they tend to be up-regulated under high light and heat while being down-regulated under cold and salt stresses. The presence of cis-regulatory motifs associated with different stresses in upstream regions of SOX genes was thus justified. The protein-protein interaction network of AtSOX and network enrichment with gene ontology (GO) terms showed that most predicted proteins, including sulfite reductase, ATP sulfurylases and APS reductases were among prime enzymes involved in sulfite pathway. Finally, SOX-sulfite docked structures indicated that arginine residues particularly Arg374 is crucial for SOX-sulfite binding and additional two other residues such as Arg51 and Arg103 may be important for SOX-sulfite bindings in plants.

Redesigning pH optimum of Geobacillus sp. TF16 endoxylanase through in silico designed DNA swapping strategy

Thermoalkaliphilic xylanases are highly desired and of great importance due to their vast potential in paper pulp and bleaching processes. Here, we report rapid, cost-effective, and result-oriented combinatorial potential of in silico DNA swapping strategy to engineer the pH optimum of industrially crucial enzymes, particularly engineering of Geobacillus sp. TF16 endoxylanase for alkaline environments. The 3D structures of Geobacillus sp. TF16 and donor Bacillus halodurans C-125 endoxylanases were firstly predicted, analyzed, and compared for their similarities before any in silico design of mutants. Reasonably, to improve its alkaline pH tolerance, the corresponding regions in Geobacillus sp.TF16 endoxylanase were further engineered by swapping with negatively-charged amino acid-rich regions from B. halodurans C-125 endoxylanase. Through only two of four in silico-designed mutants, the optimum pH of GeoTF16 endoxylanase was improved from 8.5 to 10.0. Moreover, as compared to GeoTF16 parental enzyme, both GeoInt3 and GeoInt4 mutants revealed (i) enhanced biobleaching performance, (ii) improved adaptability to alkaline conditions, and (iii) better activity for broader pH range. Unlike GeoTF16 losing activity at pH 11.0 completely, GeoInt4 retained 60% and 40% of its activity at pH 11.0 and 12.0, respectively. Thus, GeoInt4 stands out as a more competent biocatalyst that is suitable for alkaline environments of diverse industrial applications. The current study represents an efficient protein engineering strategy to adapt industrial catalysts to diverse processing conditions. Further comprehensive and fine-tuned research efforts may result in biotechnologically more promising outcomes.

Discovery of Rab1 binding sites using an ensemble of clustering methods

Targeting non-native-ligand binding sites for potential investigative and therapeutic applications is an attractive strategy in proteins that share common native ligands, as in Rab1 protein. Rab1 is a subfamily member of Rab proteins, which are members of Ras GTPase superfamily. All Ras GTPase superfamily members bind to native ligands GTP and GDP, that switch on and off the proteins, respectively. Rab1 is physiologically essential for autophagy and transport between endoplasmic reticulum and Golgi apparatus. Pathologically, Rab1 is implicated in human cancers, a neurodegenerative disease, cardiomyopathy, and bacteria-caused infectious diseases. We have performed structural analyses on Rab1 protein using a unique ensemble of clustering methods, including multi-step principal component analysis, non-negative matrix factorization, and independent component analysis, to better identify representative Rab1 proteins than the application of a single clustering method alone does. We then used the identified representative Rab1 structures, resolved in multiple ligand states, to map their known and novel binding sites. We report here at least a novel binding site on Rab1, involving Rab1-specific residues that could be further explored for the rational design and development of investigative probes and/or therapeutic small molecules against the Rab1 protein. Proteins 2017; 85:859-871. © 2016 Wiley Periodicals, Inc.

Genome-wide exploration of metal tolerance protein (MTP) genes in common wheat (Triticum aestivum): insights into metal homeostasis and biofortification

Metal transport process in plants is a determinant of quality and quantity of the harvest. Although it is among the most important of staple crops, knowledge about genes that encode for membrane-bound metal transporters is scarce in wheat. Metal tolerance proteins (MTPs) are involved in trace metal homeostasis at the sub-cellular level, usually by providing metal efflux out of the cytosol. Here, by using various bioinformatics approaches, genes that encode for MTPs in the hexaploid wheat genome (Triticum aestivum, abbreviated as Ta) were identified and characterized. Based on the comparison with known rice MTPs, the wheat genome contained 20 MTP sequences; named as TaMTP1-8A, B and D. All TaMTPs contained a cation diffusion facilitator (CDF) family domain and most members harbored a zinc transporter dimerization domain. Based on motif, phylogeny and alignment analysis, A, B and D genomes of TaMTP3-7 sequences demonstrated higher homology compared to TaMTP1, 2 and 8. With reference to their rice orthologs, TaMTP1s and TaMTP8s belonged to Zn-CDFs, TaMTP2s to Fe/Zn-CDFs and TaMTP3-7s to Mn-CDFs. Upstream regions of TaMTP genes included diverse cis-regulatory motifs, indicating regulation by developmental stage, tissue type and stresses. A scan of the coding sequences of 20 TaMTPs against published miRNAs predicted a total of 14 potential miRNAs, mainly targeting the members of most diverged groups. Expression analysis showed that several TaMTPs were temporally and spatially regulated during the developmental time-course. In grains, MTPs were preferentially expressed in the aleurone layer, which is known as a reservoir for high concentrations of iron and zinc. The work identified and characterized metal tolerance proteins in common wheat and revealed a potential involvement of MTPs in providing a sink for trace element storage in wheat grains.

RNA-TVcurve: a Web server for RNA secondary structure comparison based on a multi-scale similarity of its triple vector curve representation

Background

RNAs have been found to carry diverse functionalities in nature. Inferring the similarity between two given RNAs is a fundamental step to understand and interpret their functional relationship. The majority of functional RNAs show conserved secondary structures, rather than sequence conservation. Those algorithms relying on sequence-based features usually have limitations in their prediction performance. Hence, integrating RNA structure features is very critical for RNA analysis. Existing algorithms mainly fall into two categories: alignment-based and alignment-free. The alignment-free algorithms of RNA comparison usually have lower time complexity than alignment-based algorithms.

Results

An alignment-free RNA comparison algorithm was proposed, in which novel numerical representations RNA-TVcurve (triple vector curve representation) of RNA sequence and corresponding secondary structure features are provided. Then a multi-scale similarity score of two given RNAs was designed based on wavelet decomposition of their numerical representation. In support of RNA mutation and phylogenetic analysis, a web server (RNA-TVcurve) was designed based on this alignment-free RNA comparison algorithm. It provides three functional modules: 1) visualization of numerical representation of RNA secondary structure; 2) detection of single-point mutation based on secondary structure; and 3) comparison of pairwise and multiple RNA secondary structures. The inputs of the web server require RNA primary sequences, while corresponding secondary structures are optional. For the primary sequences alone, the web server can compute the secondary structures using free energy minimization algorithm in terms of RNAfold tool from Vienna RNA package.

Conclusion

RNA-TVcurve is the first integrated web server, based on an alignment-free method, to deliver a suite of RNA analysis functions, including visualization, mutation analysis and multiple RNAs structure comparison. The comparison results with two popular RNA comparison tools, RNApdist and RNAdistance, showcased that RNA-TVcurve can efficiently capture subtle relationships among RNAs for mutation detection and non-coding RNA classification. All the relevant results were shown in an intuitive graphical manner, and can be freely downloaded from this server. RNA-TVcurve, along with test examples and detailed documents, are available at: http://ml.jlu.edu.cn/tvcurve/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1481-7) contains supplementary material, which is available to authorized users.

Nup358 binds to AGO proteins through its SUMO-interacting motifs and promotes the association of target mRNA with miRISC

MicroRNA (miRNA)-guided mRNA repression, mediated by the miRNA-induced silencing complex (miRISC), is an important component of post-transcriptional gene silencing. However, how miRISC identifies the target mRNA in vivo is not well understood. Here, we show that the nucleoporin Nup358 plays an important role in this process. Nup358 localizes to the nuclear pore complex and to the cytoplasmic annulate lamellae (AL), and these structures dynamically associate with two mRNP granules: processing bodies (P bodies) and stress granules (SGs). Nup358 depletion disrupts P bodies and concomitantly impairs the miRNA pathway. Furthermore, Nup358 interacts with AGO and GW182 proteins and promotes the association of target mRNA with miRISC A well-characterized SUMO-interacting motif (SIM) in Nup358 is sufficient for Nup358 to directly bind to AGO proteins. Moreover, AGO and PIWI proteins interact with SIMs derived from other SUMO-binding proteins. Our study indicates that Nup358-AGO interaction is important for miRNA-mediated gene silencing and identifies SIM as a new interacting motif for the AGO family of proteins. The findings also support a model wherein the coupling of miRISC with the target mRNA could occur at AL, specialized domains within the ER, and at the nuclear envelope.

Structure-based classification of FAD binding sites: A comparative study of structural alignment tools

A total of six different structural alignment tools (TM-Align, TriangleMatch, CLICK, ProBis, SiteEngine and GA-SI) were assessed for their ability to perform two particular tasks: (i) discriminating FAD (flavin adenine dinucleotide) from non-FAD binding sites, and (ii) performing an all-to-all comparison on a set of 883 FAD binding sites for the purpose of classifying them. For the first task, the consistency of each alignment method was evaluated, showing that every method is able to distinguish FAD and non-FAD binding sites with a high Matthews correlation coefficient. Additionally, GA-SI was found to provide alignments different from those of the other approaches. The results obtained for the second task revealed more significant differences among alignment methods, as reflected in the poor correlation of their results and highlighted clearly by the independent evaluation of the structural superimpositions generated by each method. The classification itself was performed using the combined results of all methods, using the best result found for each comparison of binding sites. A number of different clustering methods (Single-linkage, UPGMA, Complete-linkage, SPICKER and k-Means clustering) were also used. The groups of similar binding sites (proteins) or clusters generated by the best performing method were further analyzed in terms of local sequence identity, local structural similarity and conservation of analogous contacts with the FAD ligands. Each of the clusters was characterized by a unique set of structural features or patterns, demonstrating that the groups generated truly reflect the structural diversity of FAD binding sites. Proteins 2016; 84:1728-1747. © 2016 Wiley Periodicals, Inc.

Topology independent comparison of RNA 3D structures using the CLICK algorithm

RNA molecules are attractive therapeutic targets because non-coding RNA molecules have increasingly been found to play key regulatory roles in the cell. Comparing and classifying RNA 3D structures yields unique insights into RNA evolution and function. With the rapid increase in the number of atomic-resolution RNA structures, it is crucial to have effective tools to classify RNA structures and to investigate them for structural similarities at different resolutions. We previously developed the algorithm CLICK to superimpose a pair of protein 3D structures by clique matching and 3D least squares fitting. In this study, we extend and optimize the CLICK algorithm to superimpose pairs of RNA 3D structures and RNA-protein complexes, independent of the associated topologies. Benchmarking Rclick on four different datasets showed that it is either comparable to or better than other structural alignment methods in terms of the extent of structural overlaps. Rclick also recognizes conformational changes between RNA structures and produces complementary alignments to maximize the extent of detectable similarity. Applying Rclick to study Ribonuclease III protein correctly aligned the RNA binding sites of RNAse III with its substrate. Rclick can be further extended to identify ligand-binding pockets in RNA. A web server is developed at http://mspc.bii.a-star.edu.sg/minhn/rclick.html.

De Novo modeling in cryo-EM density maps with Pathwalking

As electron cryo-microscopy (cryo-EM) can now frequently achieve near atomic resolution, accurate interpretation of these density maps in terms of atomistic detail has become paramount in deciphering macromolecular structure and function. However, there are few software tools for modeling protein structure from cryo-EM density maps in this resolution range. Here, we present an extension of our original Pathwalking protocol, which can automatically trace a protein backbone directly from a near-atomic resolution (3-6Å) density map. The original Pathwalking approach utilized a Traveling Salesman Problem solver for backbone tracing, but manual adjustment was still required during modeling. In the new version, human intervention is minimized and we provide a more robust approach for backbone modeling. This includes iterative secondary structure identification, termini detection and the ability to model multiple subunits without prior segmentation. Overall, the new Pathwalking procedure provides a more complete and robust tool for annotating protein structure function in near-atomic resolution density maps.

ILP-2 modeling and virtual screening of an FDA-approved library:a possible anticancer therapy

BACKGROUND/AIM

The members of the inhibitors of apoptosis protein (IAP) family inhibit diverse components of the caspase signaling pathway, notably caspase 3, 7, and 9. ILP-2 (BIRC-8) is the most recently identified member of the IAPs, mainly interacting with caspase 9. This interaction would eventually lead to death resistance in the case of cancerous cells. Therefore, structural modeling of ILP-2 and finding applicable inhibitors of its interaction with caspase 9 are a compelling challenge.

MATERIALS AND METHODS

Three main protein modeling approaches along with various model refinement measures were harnessed to achieve a reliable 3D model, using state-of-the-art software. Thereafter, the selected model was employed to perform virtual screening of an FDA approved library.

RESULTS

A model built by a combinatorial approach (homology and ab initio approaches) was chosen as the best model. Model refinement processes successfully bolstered the model quality. Virtual screening of the compound library introduced several high affinity inhibitor candidates that interact with functional residues of ILP2.

CONCLUSION

Given the 3D structure of the ILP2 molecule, we found promising inhibitory molecules. In addition to high affinity towards the ILP2 molecule, these molecules interact with residues that play pivotal rules in ILP2-caspase interaction. These molecules would inhibit ILP2-caspase interaction and consequently would lead to reactivated cell apoptosis through the caspases pathway.

Geomfinder: a multi-feature identifier of similar three-dimensional protein patterns: a ligand-independent approach

Background

Since the structure of proteins is more conserved than the sequence, the identification of conserved three-dimensional (3D) patterns among a set of proteins, can be important for protein function prediction, protein clustering, drug discovery and the establishment of evolutionary relationships. Thus, several computational applications to identify, describe and compare 3D patterns (or motifs) have been developed. Often, these tools consider a 3D pattern as that described by the residues surrounding co-crystallized/docked ligands available from X-ray crystal structures or homology models. Nevertheless, many of the protein structures stored in public databases do not provide information about the location and characteristics of ligand binding sites and/or other important 3D patterns such as allosteric sites, enzyme-cofactor interaction motifs, etc. This makes necessary the development of new ligand-independent methods to search and compare 3D patterns in all available protein structures.

Results

Here we introduce Geomfinder, an intuitive, flexible, alignment-free and ligand-independent web server for detailed estimation of similarities between all pairs of 3D patterns detected in any two given protein structures. We used around 1100 protein structures to form pairs of proteins which were assessed with Geomfinder. In these analyses each protein was considered in only one pair (e.g. in a subset of 100 different proteins, 50 pairs of proteins can be defined). Thus: (a) Geomfinder detected identical pairs of 3D patterns in a series of monoamine oxidase-B structures, which corresponded to the effectively similar ligand binding sites at these proteins; (b) we identified structural similarities among pairs of protein structures which are targets of compounds such as acarbose, benzamidine, adenosine triphosphate and pyridoxal phosphate; these similar 3D patterns are not detected using sequence-based methods; (c) the detailed evaluation of three specific cases showed the versatility of Geomfinder, which was able to discriminate between similar and different 3D patterns related to binding sites of common substrates in a range of diverse proteins.

Conclusions

Geomfinder allows detecting similar 3D patterns between any two pair of protein structures, regardless of the divergency among their amino acids sequences. Although the software is not intended for simultaneous multiple comparisons in a large number of proteins, it can be particularly useful in cases such as the structure-based design of multitarget drugs, where a detailed analysis of 3D patterns similarities between a few selected protein targets is essential.

Electronic supplementary material

The online version of this article (doi:10.1186/s13321-016-0131-9) contains supplementary material, which is available to authorized users.

Identification and Comparative Analysis of H2O2-Scavenging Enzymes (Ascorbate Peroxidase and Glutathione Peroxidase) in Selected Plants Employing Bioinformatics Approaches

Among major reactive oxygen species (ROS), hydrogen peroxide (H2O2) exhibits dual roles in plant metabolism. Low levels of H2O2 modulate many biological/physiological processes in plants; whereas, its high level can cause damage to cell structures, having severe consequences. Thus, steady-state level of cellular H2O2 must be tightly regulated. Glutathione peroxidases (GPX) and ascorbate peroxidase (APX) are two major ROS-scavenging enzymes which catalyze the reduction of H2O2 in order to prevent potential H2O2-derived cellular damage. Employing bioinformatics approaches, this study presents a comparative evaluation of both GPX and APX in 18 different plant species, and provides valuable insights into the nature and complex regulation of these enzymes. Herein, (a) potential GPX and APX genes/proteins from 18 different plant species were identified, (b) their exon/intron organization were analyzed, (c) detailed information about their physicochemical properties were provided, (d) conserved motif signatures of GPX and APX were identified, (e) their phylogenetic trees and 3D models were constructed, (f) protein-protein interaction networks were generated, and finally (g) GPX and APX gene expression profiles were analyzed. Study outcomes enlightened GPX and APX as major H2O2-scavenging enzymes at their structural and functional levels, which could be used in future studies in the current direction.

Withanolide A offers neuroprotection, ameliorates stress resistance and prolongs the life expectancy of Caenorhabditis elegans

Withanolide A (steroidal lactone) forms the major constituent of the most popular herbal drug in Ayurvedic medicine, Ashwagandha. It has been used since ancient times as an alternative medicine for the treatment of a variety of age related disorders. Here we provide multiple lines of evidence indicating that Withanolide A improves healthspan, delays age-associated physiological changes and also extends the lifespan of Caenorhabditis elegans. We also report several neuroprotective benefits of this natural product, including its anti-amyloidogenic effects, alleviation of α-synuclein aggregation and neuroprotection through modulation of neural mediators like acetylcholine. We observed that Withanolide A mediates lifespan extension and promotes stress resistance via insulin/insulin-like growth factor signaling pathway. Such findings could be helpful to develop a therapeutic medicine from this natural product for the prevention or reversal of age-related ailments and to improve the survival of patients suffering from Alzheimer's or Parkinson's disease.

Molecular docking of Glycine max and Medicago truncatula ureases with urea; bioinformatics approaches

Urease (EC 3.5.1.5) is a nickel-dependent metalloenzyme catalyzing the hydrolysis of urea into ammonia and carbon dioxide. It is present in many bacteria, fungi, yeasts and plants. Most species, with few exceptions, use nickel metalloenzyme urease to hydrolyze urea, which is one of the commonly used nitrogen fertilizer in plant growth thus its enzymatic hydrolysis possesses vital importance in agricultural practices. Considering the essentiality and importance of urea and urease activity in most plants, this study aimed to comparatively investigate the ureases of two important legume species such as Glycine max (soybean) and Medicago truncatula (barrel medic) from Fabaceae family. With additional plant species, primary and secondary structures of 37 plant ureases were comparatively analyzed using various bioinformatics tools. A structure based phylogeny was constructed using predicted 3D models of G. max and M. truncatula, whose crystallographic structures are not available, along with three additional solved urease structures from Canavalia ensiformis (PDB: 4GY7), Bacillus pasteurii (PDB: 4UBP) and Klebsiella aerogenes (PDB: 1FWJ). In addition, urease structures of these species were docked with urea to analyze the binding affinities, interacting amino acids and atom distances in urease-urea complexes. Furthermore, mutable amino acids which could potentially affect the protein active site, stability and flexibility as well as overall protein stability were analyzed in urease structures of G. max and M. truncatula. Plant ureases demonstrated similar physico-chemical properties with 833-878 amino acid residues and 89.39-90.91 kDa molecular weight with mainly acidic (5.15-6.10 pI) nature. Four protein domain structures such as urease gamma, urease beta, urease alpha and amidohydro 1 characterized the plant ureases. Secondary structure of plant ureases also demonstrated conserved protein architecture, with predominantly α-helix and random coil structures. In structure-based phylogeny, plant ureases from G. max, M. truncatula and C. ensiformis were clearly diverged from bacterial ureases of B. pasteurii and K. aerogenes. Glu, Thr, His and Gly were commonly found as interacting residues in most urease-urea docking complexes while Glu was available in all docked structures. Besides, Ala and Arg residues, which are reported in active-site architecture of plant and bacterial ureases were present in G. max urea-urease complex but not present in others. Moreover, Arg435 and Arg437 in M. truncatula and G. max, respectively were identified as highly mutable hotspot residues residing in amidohydro 1 domain of enzyme. In addition, a number of stabilizing residues were predicted upon mutation of these hotspot residues however Cys and Thr made strong implications since they were also found in codon-aligned sequences as substitutions of hotspot residues. Comparative analyses of primary sequence and secondary structure in 37 different plants demonstrated quite conserved natures of ureases in plant kingdom. Structure-based phylogeny indicated the presence of a possible prokaryote-eukaryote split and implicated the subjection of bacterial ureases to heavy selection in prokaryotic evolution compared to plants. Urea-urease docking complexes suggested that different species could share common interacting residues as well as may have some other uncommon residues at species-dependent way. In silico mutation analyses identified mutable amino acids, which were predicted to reside in catalytic site of enzyme therefore mutagenesis at these sites seemed to have adverse effects on enzyme efficiency or function. This study findings will become valuable preliminary resource for future studies to further understand the primary, secondary and tertiary structures of urease sequences in plants as well as it will provide insights about various binding features of urea-urease complexes.

In silico analysis of Mn transporters (NRAMP1) in various plant species

Manganese (Mn) is an essential micronutrient in plant life cycle. It may be involved in photosynthesis, carbohydrate and lipid biosynthesis, and oxidative stress protection. Mn deficiency inhibits the plant growth and development, and causes the various plant symptoms such as interveinal chlorosis and tissue necrosis. Despite its importance in plant life cycle, we still have limited knowledge about Mn transporters in many plant species. Therefore, this study aimed to identify and characterize high affinity Arabidopsis Mn root transporter NRAMP1 orthologs in 17 different plant species. Various in silico methods and digital gene expression data were used in identification and characterization of NRAMP1 homologs; physico-chemical properties of sequences were calculated, putative transmembrane domains (TMDs) and conserved motif signatures were determined, phylogenetic tree was constructed, 3D models and interactome map were generated, and gene expression data was analyzed. 49 NRAMP1 homologs were identified from proteome datasets of 17 plant species using AtNRAMP1 as query. Identified sequences were characterized with a NRAMP domain structure, 10-12 putative TMDs with cytosolic N- and C-terminuses, and 10-14 exons encoding a protein of 500-588 amino acids and 53.8-64.3 kDa molecular weight with basic characteristics. Consensus transport residues, GQSSTITGTYAGQY(/F)V(/I)MQGFLD(/E/N) between TMD-8 and 9 were identified in all sequences but putative N-linked glycosylation sites were not highly conserved. In phylogeny, NRAMP1 sequences demonstrated divergence in lower and higher plants as well as in monocots and dicots. Despite divergence of lower plant Physcomitrella patens in phylogeny, it showed similarity in superposed 3D models. Phylogenetic distribution of AtNRAMP1 and 6 homologs inferred a functional relationship to NRAMP6 sequences in Mn transport, while distribution of OsNRAMP1 and 5 homologs implicated an involvement of NRAMP1 sequences in Mn transport or a cross-talk between in Fe-Mn homeostasis. Interactome analysis further confirmed this cross-talk between Mn and Fe pathways. Gene expression profile of AtNRAMP1 under Fe-, K-, P- and S-deficiencies, and cold, drought, heat and salt stresses revealed various proteins involving in transcription regulation, cofactor biosynthesis, diverse developmental roles, carbohydrate metabolism, oxidation-reduction reactions, cellular signaling and protein degradation pathways. Mn deficiency or toxicity could cause serious adverse effects in plants as well as in humans. To reduce these adversities mainly rely on understanding the molecular mechanisms underlying Mn uptake from the soil. However, we still have limited knowledge regarding the structural and functional roles of Mn transporters in many plant species. Therefore, identification and characterization of Mn root uptake transporter, NRAMP1 orthologs in various plant species will provide valuable theoretical knowledge to better understand Mn transporters as well as it may become an insight for future studies aiming to develop genetically engineered and biofortified plants.

Solution Structural Studies of GTP:Adenosylcobinamide-Phosphateguanylyl Transferase (CobY) from Methanocaldococcus jannaschii

GTP:adenosylcobinamide-phosphate (AdoCbi-P) guanylyl transferase (CobY) is an enzyme that transfers the GMP moiety of GTP to AdoCbi yielding AdoCbi-GDP in the late steps of the assembly of Ado-cobamides in archaea. The failure of repeated attempts to crystallize ligand-free (apo) CobY prompted us to explore its 3D structure by solution NMR spectroscopy. As reported here, the solution structure has a mixed α/β fold consisting of seven β-strands and five α-helices, which is very similar to a Rossmann fold. Titration of apo-CobY with GTP resulted in large changes in amide proton chemical shifts that indicated major structural perturbations upon complex formation. However, the CobY:GTP complex as followed by 1H-15N HSQC spectra was found to be unstable over time: GTP hydrolyzed and the protein converted slowly to a species with an NMR spectrum similar to that of apo-CobY. The variant CobYG153D, whose GTP complex was studied by X-ray crystallography, yielded NMR spectra similar to those of wild-type CobY in both its apo- state and in complex with GTP. The CobYG153D:GTP complex was also found to be unstable over time.

Fast and accurate non-sequential protein structure alignment using a new asymmetric linear sum assignment heuristic

MOTIVATION

The three dimensional tertiary structure of a protein at near atomic level resolution provides insight alluding to its function and evolution. As protein structure decides its functionality, similarity in structure usually implies similarity in function. As such, structure alignment techniques are often useful in the classifications of protein function. Given the rapidly growing rate of new, experimentally determined structures being made available from repositories such as the Protein Data Bank, fast and accurate computational structure comparison tools are required. This paper presents SPalignNS, a non-sequential protein structure alignment tool using a novel asymmetrical greedy search technique.

RESULTS

The performance of SPalignNS was evaluated against existing sequential and non-sequential structure alignment methods by performing trials with commonly used datasets. These benchmark datasets used to gauge alignment accuracy include (i) 9538 pairwise alignments implied by the HOMSTRAD database of homologous proteins; (ii) a subset of 64 difficult alignments from set (i) that have low structure similarity; (iii) 199 pairwise alignments of proteins with similar structure but different topology; and (iv) a subset of 20 pairwise alignments from the RIPC set. SPalignNS is shown to achieve greater alignment accuracy (lower or comparable root-mean squared distance with increased structure overlap coverage) for all datasets, and the highest agreement with reference alignments from the challenging dataset (iv) above, when compared with both sequentially constrained alignments and other non-sequential alignments.

AVAILABILITY AND IMPLEMENTATION

SPalignNS was implemented in C++. The source code, binary executable, and a web server version is freely available at: http://sparks-lab.org

CONTACT

yaoqi.zhou@griffith.edu.au.

MultiSETTER: web server for multiple RNA structure comparison

Background

Understanding the architecture and function of RNA molecules requires methods for comparing and analyzing their tertiary and quaternary structures. While structural superposition of short RNAs is achievable in a reasonable time, large structures represent much bigger challenge. Therefore, we have developed a fast and accurate algorithm for RNA pairwise structure superposition called SETTER and implemented it in the SETTER web server. However, though biological relationships can be inferred by a pairwise structure alignment, key features preserved by evolution can be identified only from a multiple structure alignment. Thus, we extended the SETTER algorithm to the alignment of multiple RNA structures and developed the MultiSETTER algorithm.

Results

In this paper, we present the updated version of the SETTER web server that implements a user friendly interface to the MultiSETTER algorithm. The server accepts RNA structures either as the list of PDB IDs or as user-defined PDB files. After the superposition is computed, structures are visualized in 3D and several reports and statistics are generated.

Conclusion

To the best of our knowledge, the MultiSETTER web server is the first publicly available tool for a multiple RNA structure alignment. The MultiSETTER server offers the visual inspection of an alignment in 3D space which may reveal structural and functional relationships not captured by other multiple alignment methods based either on a sequence or on secondary structure motifs.

Approaching rational epitope vaccine design for hepatitis C virus with meta-server and multivalent scaffolding

Development of a prophylactic vaccine against hepatitis C virus (HCV) has been hampered by the extraordinary viral diversity and the poor host immune response. Scaffolding, by grafting an epitope onto a heterologous protein scaffold, offers a possible solution to epitope vaccine design. In this study, we designed and characterized epitope vaccine antigens for the antigenic sites of HCV envelope glycoproteins E1 (residues 314–324) and E2 (residues 412–423), for which neutralizing antibody-bound structures are available. We first combined six structural alignment algorithms in a “scaffolding meta-server” to search for diverse scaffolds that can structurally accommodate the HCV epitopes. For each antigenic site, ten scaffolds were selected for computational design, and the resulting epitope scaffolds were analyzed using structure-scoring functions and molecular dynamics simulation. We experimentally confirmed that three E1 and five E2 epitope scaffolds bound to their respective neutralizing antibodies, but with different kinetics. We then investigated a “multivalent scaffolding” approach by displaying 24 copies of an epitope scaffold on a self-assembling nanoparticle, which markedly increased the avidity of antibody binding. Our study thus demonstrates the utility of a multi-scale scaffolding strategy in epitope vaccine design and provides promising HCV immunogens for further assessment in vivo.

Computational tools for epitope vaccine design and evaluation

Rational approaches will be required to develop universal vaccines for viral pathogens such as human immunodeficiency virus, hepatitis C virus, and influenza, for which empirical approaches have failed. The main objective of a rational vaccine strategy is to design novel immunogens that are capable of inducing long-term protective immunity. In practice, this requires structure-based engineering of the target neutralizing epitopes and a quantitative readout of vaccine-induced immune responses. Therefore, computational tools that can facilitate these two areas have played increasingly important roles in rational vaccine design in recent years. Here we review the computational techniques developed for protein structure prediction and antibody repertoire analysis, and demonstrate how they can be applied to the design and evaluation of epitope vaccines.

A glycine insertion in the estrogen-related receptor (ERR) is associated with enhanced expression of three cytochrome P450 genes in transgenic Drosophila melanogaster

Insecticide-resistant Drosophila melanogaster strains represent a resource for the discovery of the underlying molecular mechanisms of cytochrome P450 constitutive over-expression, even if some of these P450s are not directly involved in the resistance phenotype. For example, in select 4,4'-dichlorodiphenyltrichloroethane (DDT) resistant strains the glucocorticoid receptor-like (GR-like) potential transcription factor binding motifs (TFBMs) have previously been shown to be associated with constitutively differentially-expressed cytochrome P450s, Cyp12d1, Cyp6g2 and Cyp9c1. However, insects are not known to have glucocorticoids. The only ortholog to the mammalian glucocorticoid receptor (GR) in D. melanogaster is an estrogen-related receptor (ERR) gene, which has two predicted alternative splice isoforms (ERRa and ERRb). Sequencing of ERRa and ERRb in select DDT susceptible and resistant D. melanogaster strains has revealed a glycine (G) codon insertion which was only observed in the ligand binding domain of ERR from the resistant strains tested (ERR-G). Transgenic flies, expressing the ERRa-G allele, constitutively over-expressed Cyp12d1, Cyp6g2 and Cyp9c1. Only Cyp12d1 and Cyp6g2 were over-expressed in the ERRb-G transgenic flies. Phylogenetic studies show that the G-insertion appeared to be located in a less conserved domain in ERR and this insertion is found in multiple species across the Sophophora subgenera.

A distinct sodium channel voltage-sensor locus determines insect selectivity of the spider toxin Dc1a

β-Diguetoxin-Dc1a (Dc1a) is a toxin from the desert bush spider Diguetia canities that incapacitates insects at concentrations that are non-toxic to mammals. Dc1a promotes opening of German cockroach voltage-gated sodium (Na_v) channels (BgNa_v1), whereas human Na_v channels are insensitive. Here, by transplanting commonly targeted S3b-S4 paddle motifs within BgNa_v1 voltage sensors into K_v2.1, we find that Dc1a interacts with the domain II voltage sensor. In contrast, Dc1a has little effect on sodium currents mediated by PaNa_v1 channels from the American cockroach even though their domain II paddle motifs are identical. When exploring regions responsible for PaNa_v1 resistance to Dc1a, we identified two residues within the BgNa_v1 domain II S1–S2 loop that when mutated to their PaNa_v1 counterparts drastically reduce toxin susceptibility. Overall, our results reveal a distinct region within insect Na_v channels that helps determine Dc1a sensitivity, aconcept that will be valuable for the design of insect-selective insecticides.

Protruding knob-like proteins violate local symmetries in an icosahedral marine virus

Marine viruses play crucial roles in shaping the dynamics of oceanic microbial communities and in the carbon cycle on Earth. Here we report a 4.7-Å structure of a cyanobacterial virus, Syn5, by electron cryo-microscopy and modelling. A Cα backbone trace of the major capsid protein (gp39) reveals a classic phage protein fold. In addition, two knob-like proteins protruding from the capsid surface are also observed. Using bioinformatics and structure analysis tools, these proteins are identified to correspond to gp55 and gp58 (each with two copies per asymmetric unit). The non 1:1 stoichiometric distribution of gp55/58 to gp39 breaks all expected local symmetries and leads to non-quasi-equivalence of the capsid subunits, suggesting a role in capsid stabilization. Such a structural arrangement has not yet been observed in any known virus structures.

Multiple functional roles of the accessory I-domain of bacteriophage P22 coat protein revealed by NMR structure and CryoEM modeling

Some capsid proteins built on the ubiquitous HK97-fold have accessory domains imparting specific functions. Bacteriophage P22 coat protein has a unique insertion domain (I-domain). Two prior I-domain models from subnanometer cryoelectron microscopy (cryoEM) reconstructions differed substantially. Therefore, the I-domain's nuclear magnetic resonance structure was determined and also used to improve cryoEM models of coat protein. The I-domain has an antiparallel six-stranded β-barrel fold, not previously observed in HK97-fold accessory domains. The D-loop, which is dynamic in the isolated I-domain and intact monomeric coat protein, forms stabilizing salt bridges between adjacent capsomers in procapsids. The S-loop is important for capsid size determination, likely through intrasubunit interactions. Ten of 18 coat protein temperature-sensitive-folding substitutions are in the I-domain, indicating its importance in folding and stability. Several are found on a positively charged face of the β-barrel that anchors the I-domain to a negatively charged surface of the coat protein HK97-core.

A high-affinity protein binder that blocks the IL-6/STAT3 signaling pathway effectively suppresses non-small cell lung cancer

Interleukin-6 (IL-6) is a multifunctional cytokine that regulates immune responses for host defense and tumorigenic process. Upregulation of IL-6 is known to constitutively phosphorylate signal transducer and activator of transcription 3 (STAT3), leading to activation of multiple oncogene pathways and inflammatory cascade. Here, we present the development of a high-affinity protein binder, termed repebody, which effectively suppresses non-small cell lung cancer in vivo by blocking the IL-6/STAT3 signaling. We selected a repebody that prevents human IL-6 (hIL-6) from binding to its receptor by a competitive immunoassay, and modulated its binding affinity for hIL-6 up to a picomolar range by a modular approach that mimics the combinatorial assembly of diverse modules to form antigen-specific receptors in nature. The resulting repebody was highly specific for hIL-6, effectively inhibiting the STAT3 phosphorylation in a dose- and binding affinity-response manner in vitro. The repebody was shown to have a remarkable suppression effect on the growth of tumors and STAT3 phosphorylation in xenograft mice with non-small cell lung cancer by blocking the hIL-6/STAT3 signaling. Structural analysis of the repebody and IL-6 complex revealed that the repebody binds the site 2a of hIL-6, overlapping a number of epitope residues at site 2a with gp130, and consequently causes a steric hindrance to the formation of IL-6/IL-6Rα complex. Our results suggest that high-affinity repebody targeting the IL-6/STAT3 pathway can be developed as therapeutics for non-small cell lung cancer.

Influence of hydrophobic and electrostatic residues on SARS-coronavirus S2 protein stability: insights into mechanisms of general viral fusion and inhibitor design

Severe acute respiratory syndrome (SARS) is an acute respiratory disease caused by the SARS-coronavirus (SARS-CoV). SARS-CoV entry is facilitated by the spike protein (S), which consists of an N-terminal domain (S1) responsible for cellular attachment and a C-terminal domain (S2) that mediates viral and host cell membrane fusion. The SARS-CoV S2 is a potential drug target, as peptidomimetics against S2 act as potent fusion inhibitors. In this study, site-directed mutagenesis and thermal stability experiments on electrostatic, hydrophobic, and polar residues to dissect their roles in stabilizing the S2 postfusion conformation was performed. It was shown that unlike the pH-independent retroviral fusion proteins, SARS-CoV S2 is stable over a wide pH range, supporting its ability to fuse at both the plasma membrane and endosome. A comprehensive SARS-CoV S2 analysis showed that specific hydrophobic positions at the C-terminal end of the HR2, rather than electrostatics are critical for fusion protein stabilization. Disruption of the conserved C-terminal hydrophobic residues destabilized the fusion core and reduced the melting temperature by 30°C. The importance of the C-terminal hydrophobic residues led us to identify a 42-residue substructure on the central core that is structurally conserved in all existing CoV S2 fusion proteins (root mean squared deviation = 0.4 Å). This is the first study to identify such a conserved substructure and likely represents a common foundation to facilitate viral fusion. We have discussed the role of key residues in the design of fusion inhibitors and the potential of the substructure as a general target for the development of novel therapeutics against CoV infections.

A Computational Approach towards the Understanding of Plasmodium falciparum Multidrug Resistance Protein 1

The emergence of drug resistance in Plasmodium falciparum tremendously affected the chemotherapy worldwide while the intense distribution of chloroquine-resistant strains in most of the endemic areas added more complications in the treatment of malaria. The situation has even worsened by the lack of molecular mechanism to understand the resistance conferred by Plasmodia species. Recent studies have suggested the association of antimalarial resistance with P. falciparum multidrug resistance protein 1 (PfMDR1), an ATP-binding cassette (ABC) transporter and a homologue of human P-glycoprotein 1 (P-gp1). The present study deals about the development of PfMDR1 computational model and the model of substrate transport across PfMDR1 with insights derived from conformations relative to inward- and outward-facing topologies that switch on/off the transportation system. Comparison of ATP docked positions and its structural motif binding properties were found to be similar among other ATPases, and thereby contributes to NBD domains dimerization, a unique structural agreement noticed in Mus musculus Pgp and Escherichia coli MDR transporter homolog (MsbA). The interaction of leading antimalarials and phytochemicals within the active pocket of both wild-type and mutant-type PfMDR1 demonstrated the mode of binding and provided insights of less binding affinity thereby contributing to parasite's resistance mechanism.

Validated near-atomic resolution structure of bacteriophage epsilon15 derived from cryo-EM and modeling

High-resolution structures of viruses have made important contributions to modern structural biology. Bacteriophages, the most diverse and abundant organisms on earth, replicate and infect all bacteria and archaea, making them excellent potential alternatives to antibiotics and therapies for multidrug-resistant bacteria. Here, we improved upon our previous electron cryomicroscopy structure of Salmonella bacteriophage epsilon15, achieving a resolution sufficient to determine the tertiary structures of both gp7 and gp10 protein subunits that form the T = 7 icosahedral lattice. This study utilizes recently established best practice for near-atomic to high-resolution (3-5 Å) electron cryomicroscopy data evaluation. The resolution and reliability of the density map were cross-validated by multiple reconstructions from truly independent data sets, whereas the models of the individual protein subunits were validated adopting the best practices from X-ray crystallography. Some sidechain densities are clearly resolved and show the subunit-subunit interactions within and across the capsomeres that are required to stabilize the virus. The presence of the canonical phage and jellyroll viral protein folds, gp7 and gp10, respectively, in the same virus suggests that epsilon15 may have emerged more recently relative to other bacteriophages.

Depth: a web server to compute depth, cavity sizes, detect potential small-molecule ligand-binding cavities and predict the pKa of ionizable residues in proteins

Residue depth accurately measures burial and parameterizes local protein environment. Depth is the distance of any atom/residue to the closest bulk water. We consider the non-bulk waters to occupy cavities, whose volumes are determined using a Voronoi procedure. Our estimation of cavity sizes is statistically superior to estimates made by CASTp and VOIDOO, and on par with McVol over a data set of 40 cavities. Our calculated cavity volumes correlated best with the experimentally determined destabilization of 34 mutants from five proteins. Some of the cavities identified are capable of binding small molecule ligands. In this study, we have enhanced our depth-based predictions of binding sites by including evolutionary information. We have demonstrated that on a database (LigASite) of ∼200 proteins, we perform on par with ConCavity and better than MetaPocket 2.0. Our predictions, while less sensitive, are more specific and precise. Finally, we use depth (and other features) to predict pK_as of GLU, ASP, LYS and HIS residues. Our results produce an average error of just <1 pH unit over 60 predictions. Our simple empirical method is statistically on par with two and superior to three other methods while inferior to only one. The DEPTH server (http://mspc.bii.a-star.edu.sg/depth/) is an ideal tool for rapid yet accurate structural analyses of protein structures.

Structural basis for the recruitment of the human CCR4-NOT deadenylase complex by tristetraprolin

Tristetraprolin (TTP) is an RNA-binding protein that controls the inflammatory response by limiting the expression of several proinflammatory cytokines. TTP post-transcriptionally represses gene expression by interacting with AU-rich elements (AREs) in 3' untranslated regions of target mRNAs and subsequently engenders their deadenylation and decay. TTP accomplishes these tasks, at least in part, by recruiting the multisubunit CCR4-NOT deadenylase complex to the mRNA. Here we identify an evolutionarily conserved C-terminal motif in human TTP that directly binds a central domain of CNOT1, a core subunit of the CCR4-NOT complex. A high-resolution crystal structure of the TTP-CNOT1 complex was determined, providing the first structural insight, to our knowledge, into an ARE-binding protein bound to the CCR4-NOT complex. Mutations at the CNOT1-TTP interface impair TTP-mediated deadenylation, demonstrating the significance of this interaction in TTP-mediated gene silencing.

Structural phylogeny by profile extraction and multiple superimposition using electrostatic congruence as a discriminator

Phylogenetic analysis of proteins using multiple sequence alignment (MSA) assumes an underlying evolutionary relationship in these proteins which occasionally remains undetected due to considerable sequence divergence. Structural alignment programs have been developed to unravel such fuzzy relationships. However, none of these structure based methods have used electrostatic properties to discriminate between spatially equivalent residues. We present a methodology for MSA of a set of related proteins with known structures using electrostatic properties as an additional discriminator (STEEP). STEEP first extracts a profile, then generates a multiple structural superimposition providing a consolidated spatial framework for comparing residues and finally emits the MSA. Residues that are aligned differently by including or excluding electrostatic properties can be targeted by directed evolution experiments to transform the enzymatic properties of one protein into another. We have compared STEEP results to those obtained from a MSA program (ClustalW) and a structural alignment method (MUSTANG) for chymotrypsin serine proteases. Subsequently, we used PhyML to generate phylogenetic trees for the serine and metallo-β-lactamase superfamilies from the STEEP generated MSA, and corroborated the accepted relationships in these superfamilies. We have observed that STEEP acts as a functional classifier when electrostatic congruence is used as a discriminator, and thus identifies potential targets for directed evolution experiments. In summary, STEEP is unique among phylogenetic methods for its ability to use electrostatic congruence to specify mutations that might be the source of the functional divergence in a protein family. Based on our results, we also hypothesize that the active site and its close vicinity contains enough information to infer the correct phylogeny for related proteins.

Gorgon and pathwalking: macromolecular modeling tools for subnanometer resolution density maps

The complex interplay of proteins and other molecules, often in the form of large transitory assemblies, are critical to cellular function. Today, X-ray crystallography and electron cryo-microscopy (cryo-EM) are routinely used to image these macromolecular complexes, though often at limited resolutions. Despite the rapidly growing number of macromolecular structures, few tools exist for modeling and annotating structures in the range of 3-10 Å resolution. To address this need, we have developed a number of utilities specifically targeting subnanometer resolution density maps. As part of the 2010 Cryo-EM Modeling Challenge, we demonstrated two of our latest de novo modeling tools, Pathwalking and Gorgon, as well as a tool for secondary structure identification (SSEHunter) and a new rigid-body/flexible fitting tool in Gorgon. In total, we submitted 30 structural models from ten different subnanometer resolution data sets in four of the six challenge categories. Each of our utlities produced accurate structural models and annotations across the various density maps. In the end, the utilities that we present here offer users a robust toolkit for analyzing and modeling protein structure in macromolecular assemblies at non-atomic resolutions.

Constructing and validating initial Cα models from subnanometer resolution density maps with pathwalking

A significant number of macromolecular structures solved by electron cryo-microscopy and X-ray crystallography obtain resolutions of 3.5-6Å, at which direct atomistic interpretation is difficult. To address this, we developed pathwalking, a semi-automated protocol to enumerate reasonable Cα models from near-atomic resolution density maps without a structural template or sequence-structure correspondence. Pathwalking uses an approach derived from the Traveling Salesman Problem to rapidly generate an ensemble of initial models for individual proteins, which can later be optimized to produce full atomic models. Pathwalking can also be used to validate and identify potential structural ambiguities in models generated from near-atomic resolution density maps. In this work, examples from the EMDB and PDB are used to assess the broad applicability and accuracy of our method. With the growing number of near-atomic resolution density maps from cryo-EM and X-ray crystallography, pathwalking can become an important tool in modeling protein structures.

NASSAM: a server to search for and annotate tertiary interactions and motifs in three-dimensional structures of complex RNA molecules

Similarities in the 3D patterns of RNA base interactions or arrangements can provide insights into their functions and roles in stabilization of the RNA 3D structure. Nucleic Acids Search for Substructures and Motifs (NASSAM) is a graph theoretical program that can search for 3D patterns of base arrangements by representing the bases as pseudo-atoms. The geometric relationship of the pseudo-atoms to each other as a pattern can be represented as a labeled graph where the pseudo-atoms are the graph’s nodes while the edges are the inter-pseudo-atomic distances. The input files for NASSAM are PDB formatted 3D coordinates. This web server can be used to identify matches of base arrangement patterns in a query structure to annotated patterns that have been reported in the literature or that have possible functional and structural stabilization implications. The NASSAM program is freely accessible without any login requirement at http://mfrlab.org/grafss/nassam/.