STING Millennium: A web-based suite of programs for comprehensive and simultaneous analysis of protein structure and sequence

STING Millennium Suite (SMS) is a new web-based suite of programs and databases providing visualization and a complex analysis of molecular sequence and structure for the data deposited at the Protein Data Bank (PDB). SMS operates with a collection of both publicly available data (PDB, HSSP, Prosite) and its own data (contacts, interface contacts, surface accessibility). Biologists find SMS useful because it provides a variety of algorithms and validated data, wrapped-up in a user friendly web interface. Using SMS it is now possible to analyze sequence to structure relationships, the quality of the structure, nature and volume of atomic contacts of intra and inter chain type, relative conservation of amino acids at the specific sequence position based on multiple sequence alignment, indications of folding essential residue (FER) based on the relationship of the residue conservation to the intra-chain contacts and Calpha-Calpha and Cbeta-Cbeta distance geometry. Specific emphasis in SMS is given to interface forming residues (IFR)-amino acids that define the interactive portion of the protein surfaces. SMS may simultaneously display and analyze previously superimposed structures. PDB updates trigger SMS updates in a synchronized fashion. SMS is freely accessible for public data at http://www.cbi.cnptia.embrapa.br, http://mirrors.rcsb.org/SMS and http://trantor.bioc.columbia.edu/SMS.     

STING Contacts: a web-based application for identification and analysis of amino acid contacts within protein structure and across protein interfaces

Amino acid contacts in terms of atomic interactions are essential factors to be considered in the analysis of the structure of a protein and its complexes. Consequently, molecular biologists do require specific tools for the identification and visualization of all such contacts. Graphical contacts (GC) and interface forming residue graphical contacts (IFRgc) presented here, calculate atomic contacts among amino acids based on a table of predefined pairs of the atom types and their distances, and then display them using number of different forms. The inventory of currently listed contact types by GC and IFRgc include hydrogen bonds (in nine different flavors), hydrophobic interactions, charge-charge interactions, aromatic stacking and disulfide bonds. Such extensive catalog of the interactions, representing the forces that govern protein folding, stability and binding, is the key feature of these two applications. GC and IFRgc are part of STING Millennium Suite. AVAILABILITY: http://sms.cbi.cnptia.embrapa.br/SMS, http://trantor.bioc.columbia.edu/SMS, http://mirrors.rcsb.org//SMS, http://www.es.embnet.org/SMS and http://www.ar.embnet.org/SMS (Options: Graphical Contacts and IFR Graphical Contacts).     

Defining 3D residue environment in protein structures using SCORPION and FORMIGA

SUMMARY: Two web-based applications to analyze amino acids three-dimensional (3D) local environment within protein structures-SCORPION and FORMIGA-are presented. SCORPION and FORMIGA produce a graphical presentation for simple statistical data showing the frequency of residue occurrence within a given sphere (defined here as the 3D contacts). The center of that sphere is placed at the Calpha and at the last heavy atom in the side chain of the selected amino acid. Further depth of detail is given in terms of a secondary structure to which the profiled amino acid belongs. Results obtained with those two applications are relevant for estimating the importance of the amino acid 3D local environment for protein folding and stability. Effectively, SCORPION and FORMIGA construct knowledge-based force fields. The difference between SCORPION and FORMIGA is in that the latter operates on protein interfaces, while the former only functions for a single protein chain. Both applications are implemented as stand-alone components of STING Millennium Suite. AVAILABILITY: http://sms.cbi.cnptia.embrapa.br/SMS, http://trantor.bioc.columbia.edu/SMS, http://mirrors.rcsb.org/SMS, http://www.es.embnet.org/SMS and http://www.ar.embnet.org/SMS. [options: Scorpion, Formiga]     

VoroMQA: Assessment of protein structure quality using interatomic contact areas

In the absence of experimentally determined protein structure many biological questions can be addressed using computational structural models. However, the utility of protein structural models depends on their quality. Therefore, the estimation of the quality of predicted structures is an important problem. One of the approaches to this problem is the use of knowledge‐based statistical potentials. Such methods typically rely on the statistics of distances and angles of residue‐residue or atom‐atom interactions collected from experimentally determined structures. Here, we present VoroMQA (Voronoi tessellation‐based Model Quality Assessment), a new method for the estimation of protein structure quality. Our method combines the idea of statistical potentials with the use of interatomic contact areas instead of distances. Contact areas, derived using Voronoi tessellation of protein structure, are used to describe and seamlessly integrate both explicit interactions between protein atoms and implicit interactions of protein atoms with solvent. VoroMQA produces scores at atomic, residue, and global levels, all in the fixed range from 0 to 1. The method was tested on the CASP data and compared to several other single‐model quality assessment methods. VoroMQA showed strong performance in the recognition of the native structure and in the structural model selection tests, thus demonstrating the efficacy of interatomic contact areas in estimating protein structure quality. The software implementation of VoroMQA is freely available as a standalone application and as a web server at http://bioinformatics.lt/software/voromqa . Proteins 2017; 85:1131–1145. © 2017 Wiley Periodicals, Inc.     

Parameter optimized surfaces (POPS): analysis of key interactions and conformational changes in the ribosome

We present a new method for the calculation of solvent accessible surface areas at the atomic and residue levels, which we call parameter optimized surfaces (POPS-A and POPS-R ). Atomic and residue areas (the latter simulated with a single sphere centered at the C^αs atom for amino acids and at the P atom for nucleotides) have been optimized versus accurate all-atoms methods. We concentrated on an analytical formula for the approximation of solvent accessibilities. The formula is simple, easily derivable and fast to compute, therefore it is practical for use in molecular dynamics simulations as an approximation to the first solvation shell. The residue based approach POPS-R has been derived as a useful tool for the analysis of large macromolecular assemblies like the ribosome, and is especially suited for use in refinement of low resolution structures. The structures of the 70S, 50S and 30S ribosomes have been analyzed in detail and most of the interactions within the subunits and at their interfaces were clearly identified. Some interesting differences between 30S alone and within the 70S have been highlighted. Owing to the presence of the P-tRNA in the 70S ribosome, localized conformational rearrangements occur within the subunits, exposing Arg and Lys residues to negatively charged binding sites of P-tRNA. POPS-R also allows for estimates of the loss of free energy of solvation upon complex formation, particularly useful in designing new protein–RNA complexes and in suggesting more focused experimental work.     

Molecular mechanism study of several inhibitors binding to BRD9 bromodomain based on molecular simulations

Bromodomain-containing protein 9 (BRD9) has been employed as a potential target for anticancer drugs in recent years. In this work, molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and per residue energy decomposition approaches were performed to elucidate the different binding modes between four pyridinone-like scaffold inhibitors and BRD9 bromodomain. Analysis results indicate that non-polar contribution mainly deriving from van der Waals energy is a critical impact on binding affinity of inhibitors against BRD9. Some key residues Phe44, Phe47, Val49, and Ile53 (at ZA loop) enhance the binding energy of inhibitors in BRD9 by means of providing hydrophobic interactions. Moreover, it is observed that BRD9 is anchored by the formation of a stable hydrogen bond between the carbonyl of the inhibitors and the residue Asn100 (at BC loop), and a strong π–π stacking interaction formed between the residue Tyr106 (at BC loop) and the inhibitors. The existence of dimethoxyphenyl structure and the aromatic ring merged to pyridinone scaffold are useful to enhance the BRD9 binding affinity. These findings should guide the rational design of more prospective inhibitors targeting BRD9.

Communicated by Ramaswamy H. Sarma     

Functional Reconstitution of Plasma Membrane H+-ATPase from Mung Bean (Vigna radiata L.) Hypocotyls in Liposomes Prepared with Various Molecular Species of Phospholipids

The activity of solubilized plasma membrane ATPase is affectedby the nature of exogenously added molecular species of phospholipids.To examine the role of the polar head group and of the molecularspecies of phospholipids in H⁺-pumping, the ATPase solubilizedfrom plasma membranes of mung bean (Vigna radiata L.) hypocotylswas reconstituted in liposomes prepared with a variety of phospholipids. The extent of activation of solubilized plasma membrane ATPasedue to the addition of 1-palmitoyl 2-oleoyl-phospholipids (PO-phospholipids)and asolectin decreased in the following order: POPS POPC asolectin POPG > POPE > POPA (see List of Abbreviations). H⁺-pumpinginto proteoliposomes reconstituted with asolectin and plasmamembrane ATPase was demonstrated by quinacrine fluorescencequenching in the presence of ATP-MgSO₄. H⁺-pumping was inhibitedby VO₄ and gramicidin D. When plasma membrane ATPase was reconstitutedin liposomes prepared with various PO-phospholipids, the abilityof PO-phospholipids to support H⁺-pumping into the proteoliposomesdecreased in the following order: POPG POPS > asolectin POPC. POPE and POPA failed to support any H⁺-pumping. A remarkablyhigh rate of H⁺-pumping was observed in proteoliposomes preparedwith 1-saturated 2-unsaturated fatty acids, such as POPC, butH⁺-pumping could hardly be detected in proteoliposomes preparedwith 1-, 2-unsaturated or 1-, 2-saturated fatty acids, suchas PSPC or DLPC. ATPase activity in proteoliposomes was dependenton the species of PO-phospholipids used for reconstitution anddecreased in the following order: POPS > POPG > POPC asolectin > POPA > POPE. DLPC (see List of Abbreviations)which includes a 1-, 2-unsaturated fatty acid supported onlymarkedly depressed activity. Both H⁺-pumping and the hydrolysis of ATP by the plasma membraneATPase are strongly affected by the polar head group and compositionof the fatty acyl chain of phospholipids used to prepare liposomesfor reconstitution of the ATPase. (Received May 31, 1991; Accepted September 18, 1991)     

Rapid assessment of correlated amino acids from pair-to-pair (P2P) substitution matrices

Identification of correlated amino acids in proteins has been a topic of broad interest in view of its functional implications and importance in protein design. A new set of pair-to-pair (P2P) substitution matrices for amino acids was recently introduced as a useful tool for inferring information on such correlated sites. We present a website developed for automated application of these matrices for analysis of query sequences. The site offers options for graphical analysis of correlations, as well as visualization of correlated amino acids on representative, structurally characterized, members of the examined family of sequences. Availability: http://www.ccbb.pitt.edu/p2p.     

sGAL: a computational method for finding surface exposed sites in proteins suitable for Cys-mediated cross-linking

sGAL is a computer program designed to find pairs of sites suitable for introducing chemical cross-links into proteins. sGAL takes a protein structure file in PDB format as input, truncates each residue sequentially to its gamma side chain atom to mimic mutation to Cys, and calculates the exposed surface area of the gamma atom. The user then inputs the minimum and maximum lengths of the cross-linker. sGAL provides as output pairs of residues that would have exposed gamma atom separations that fall within this range. Furthermore, if a line joining the pair of gamma atoms contacts more than a given number of buried atoms, that pair is discarded. In this way, sites for which the protein would sterically interfere with cross-linking are avoided. AVAILABILITY: http://www.chem.utoronto.ca/staff/GAW/links.html; (Surface Racer is also required see: http://monte.biochem.wisc.edu/~tsodikov/surface.html).     

CODEHOP (COnsensus-DEgenerate Hybrid Oligonucleotide Primer) PCR primer design

We have developed a new primer design strategy for PCR amplification of distantly related gene sequences based on consensus-degenerate hybrid oligonucleotide primers (CODEHOPs). An interactive program has been written to design CODEHOP PCR primers from conserved blocks of amino acids within multiply-aligned protein sequences. Each CODEHOP consists of a pool of related primers containing all possible nucleotide sequences encoding 3-4 highly conserved amino acids within a 3' degenerate core. A longer 5' non-degenerate clamp region contains the most probable nucleotide predicted for each flanking codon. CODEHOPs are used in PCR amplification to isolate distantly related sequences encoding the conserved amino acid sequence. The primer design software and the CODEHOP PCR strategy have been utilized for the identification and characterization of new gene orthologs and paralogs in different plant, animal and bacterial species. In addition, this approach has been successful in identifying new pathogen species. The CODEHOP designer (http://blocks.fhcrc.org/codehop.html) is linked to BlockMaker and the Multiple Alignment Processor within the Blocks Database World Wide Web (http://blocks.fhcrc.org).     

Disulfide by Design: a computational method for the rational design of disulfide bonds in proteins

SUMMARY: Disulfide by Design is a program for the design of novel disulfide bonds in proteins. Protein structure files in PDB format are analyzed to identify residue pairs that are likely to form a disulfide bond if the respective amino acids are mutated to cysteines. The output displays residue pairs having the appropriate geometric characteristics for disulfide formation and provides automated generation of modified PDB files including modeled disulfides. Validation demonstrates a high level of accuracy for the algorithm. AVAILABILITY: http://www.ehscenter.org/dbd/ Supplementary information: http://www.ehscenter.org/dbd/     

ConSSeq: a web-based application for analysis of amino acid conservation based on HSSP database and within context of structure

SUMMARY: A web-based application to analyze protein amino acids conservation-Consensus Sequence (ConSSeq) is presented. ConSSeq graphically represents information about amino acid conservation based on sequence alignments reported in homology-derived structures of proteins. Beyond the relative entropy for each position in the alignment, ConSSeq also presents the consensus sequence and information about the amino acids, which are predominant at each position of the alignment. ConSSeq is part of the STING Millennium Suite and is implemented as a Java Applet. AVAILABILITY: http://sms.cbi.cnptia.embrapa.br/SMS/STINGm/consseq/, http://trantor.bioc.columbia.edu/SMS/STINGm/consseq/, http://mirrors.rcsb.org//SMS/STINGm/consseq/, http://www.es.embnet.org/SMS/STINGm/consseq/ and http://www.ar.embnet.org/SMS/STINGm/consseq/     

Automated analysis of interatomic contacts in proteins.

MOTIVATION: New software has been designed to assist the molecular biologist in understanding the structural consequences of modifying a ligand and/or protein. RESULTS: Tools are described for the analysis of ligand-protein contacts (LPC software) and contacts of structural units (CSU software) such as helices, sheets, strands and residues. Our approach is based on a detailed analysis of interatomic contacts and interface complementarity. For any ligand or structural unit, these software automatically: (i) calculate the solvent-accessible surface of every atom; (ii) determine the contacting residues and type of interaction they undergo (hydrophobic-hydrophobic, aromatic-aromatic, etc.); (iii) indicate all putative hydrogen bonds. LPC software further predicts changes in binding strength following chemical modification of the ligand. AVAILABILITY: Both LPC and CSU can be accessed through the PDB and are integrated in the 3DB Atlas page of all PDB files. For any given file, the tools can also be accessed at http://www.pdb.bnl. gov/pdb-bin/lpc?PDB_ID= and http://www.pdb.bnl. gov/pdb-bin/csu?PDB_ID= with the four-letter PDB code added at the end in each case. Finally, LPC and CSU can be accessed at: http://sgedg.weizmann.ac.il/lpc and http://sgedg.weizmann.ac.il/csu.     

A New Supervised Over-Sampling Algorithm with Application to Protein-Nucleotide Binding Residue Prediction

Protein-nucleotide interactions are ubiquitous in a wide variety of biological processes. Accurately identifying interaction residues solely from protein sequences is useful for both protein function annotation and drug design, especially in the post-genomic era, as large volumes of protein data have not been functionally annotated. Protein-nucleotide binding residue prediction is a typical imbalanced learning problem, where binding residues are extremely fewer in number than non-binding residues. Alleviating the severity of class imbalance has been demonstrated to be a promising means of improving the prediction performance of a machine-learning-based predictor for class imbalance problems. However, little attention has been paid to the negative impact of class imbalance on protein-nucleotide binding residue prediction. In this study, we propose a new supervised over-sampling algorithm that synthesizes additional minority class samples to address class imbalance. The experimental results from protein-nucleotide interaction datasets demonstrate that the proposed supervised over-sampling algorithm can relieve the severity of class imbalance and help to improve prediction performance. Based on the proposed over-sampling algorithm, a predictor, called TargetSOS, is implemented for protein-nucleotide binding residue prediction. Cross-validation tests and independent validation tests demonstrate the effectiveness of TargetSOS. The web-server and datasets used in this study are freely available at http://www.csbio.sjtu.edu.cn/bioinf/TargetSOS/.     

Fatty acids from oleaginous yeasts and yeast-like fungi and their potential applications

Purpose: Oleaginous yeasts, fatty acids biosynthesis and regulation in the oleaginous yeasts and the fatty acids from the oleaginous yeasts and their applications are reviewed in this article.

Results: Oleaginous yeasts such as Rhodosporidium toruloides, Yarrowia lipolytica, Rhodotorula mucilaginosa, and Aureobasidium melanogenum, which can accumulate over 50% lipid of their cell dry weight, have many advantages over other oleaginous microorganisms. The fatty acids from the oleaginous yeasts have many potential applications. Many oleaginous yeasts have now been genetically modified to over-produce fatty acids and their derivatives. The most important features of the oleaginous yeasts are that they have special enzymatic systems for enhanced biosynthesis and regulation of fatty acids in their lipid particles. Recently, some oleaginous yeasts such as R. toruloides have been found to have a unique fatty acids synthetase and other oleaginous yeasts such as A. melanogenum have a unique highly reducing polyketide synthase (HR-PKS) involved in the biosynthesis of hydroxyl fatty acids.

Conclusions: It is necessary to further enhance lipid biosynthesis using metabolic engineering and explore new applications of fatty acids in biotechnology.