Motivation

Circular permutation is an important type of protein rearrangement. Natural circular permutations have implications for protein function, stability and evolution. Artificial circular permutations have also been used for protein studies. However, such relationships are difficult to detect for many sequence and structure comparison algorithms and require special consideration.

Results

We developed a new algorithm, called Combinatorial Extension for Circular Permutations (CE-CP), which allows the structural comparison of circularly permuted proteins. CE-CP was designed to be user friendly and is integrated into the RCSB Protein Data Bank. It was tested on two collections of circularly permuted proteins. Pairwise alignments can be visualized both in a desktop application or on the web using Jmol and exported to other programs in a variety of formats.

Availability and implementation

The CE-CP algorithm can be accessed through the RCSB website at http://www.rcsb.org/pdb/workbench/workbench.do. Source code is available under the LGPL 2.1 as part of BioJava 3 (http://biojava.org; http://github.com/biojava/biojava).

Contact

sbliven@ucsd.edu or info@rcsb.org.

Many proteins fold into highly regular and repetitive three dimensional structures. The analysis of structural patterns and repeated elements is fundamental to understand protein function and evolution. We present recent improvements to the CE-Symm tool for systematically detecting and analyzing the internal symmetry and structural repeats in proteins. In addition to the accurate detection of internal symmetry, the tool is now capable of i) reporting the type of symmetry, ii) identifying the smallest repeating unit, iii) describing the arrangement of repeats with transformation operations and symmetry axes, and iv) comparing the similarity of all the internal repeats at the residue level. CE-Symm 2.0 helps the user investigate proteins with a robust and intuitive sequence-to-structure analysis, with many applications in protein classification, functional annotation and evolutionary studies. We describe the algorithmic extensions of the method and demonstrate its applications to the study of interesting cases of protein evolution.

Few long-term multicenter investigations have evaluated the relationships between aggregate antimicrobial drug use in hospitals and bacterial resistance. We measured fluoroquinolone use from 1999 through 2003 in a network of US hospitals. The percentages of fluoroquinolone-resistant Pseudomonas aeruginosa and methicillin-resistant Staphylococcus aureus (MRSA) were obtained from yearly antibiograms at each hospital. Univariate linear regression showed significant associations between a hospital's volume of fluoroquinolone use and percent resistance in most individual study years (1999-2001 for P. aeruginosa, 1999-2002 for S. aureus). When the method of generalized estimating equations was used, a population-averaged longitudinal model incorporating total fluoroquinolone use and the previous year's resistance (to account for autocorrelation) did not show a significant effect of fluoroquinolone use on percent resistance for most drug-organism combinations, except for the relationship between levofloxacin use and percent MRSA. The ecologic relationship between fluoroquinolone use and resistance is complex and requires further study.

Symmetry is an important feature of protein tertiary and quaternary structures that has been associated with protein folding, function, evolution, and stability. Its emergence and ensuing prevalence has been attributed to gene duplications, fusion events, and subsequent evolutionary drift in sequence. This process maintains structural similarity and is further supported by this study. To further investigate the question of how internal symmetry evolved, how symmetry and function are related, and the overall frequency of internal symmetry, we developed an algorithm, CE-Symm, to detect pseudo-symmetry within the tertiary structure of protein chains. Using a large manually curated benchmark of 1007 protein domains, we show that CE-Symm performs significantly better than previous approaches. We use CE-Symm to build a census of symmetry among domain superfamilies in SCOP and note that 18% of all superfamilies are pseudo-symmetric. Our results indicate that more domains are pseudo-symmetric than previously estimated. We establish a number of recurring types of symmetry-function relationships and describe several characteristic cases in detail. With the use of the Enzyme Commission classification, symmetry was found to be enriched in some enzyme classes but depleted in others. CE-Symm thus provides a methodology for a more complete and detailed study of the role of symmetry in tertiary protein structure [availability: CE-Symm can be run from the Web at http://source.rcsb.org/jfatcatserver/symmetry.jsp. Source code and software binaries are also available under the GNU Lesser General Public License (version 2.1) at https://github.com/rcsb/symmetry. An interactive census of domains identified as symmetric by CE-Symm is available from http://source.rcsb.org/jfatcatserver/scopResults.jsp].

Background

The success of genome-scale models (GEMs) can be attributed to the high-quality, bottom-up reconstructions of metabolic, protein synthesis, and transcriptional regulatory networks on an organism-specific basis. Such reconstructions are biochemically, genetically, and genomically structured knowledge bases that can be converted into a mathematical format to enable a myriad of computational biological studies. In recent years, genome-scale reconstructions have been extended to include protein structural information, which has opened up new vistas in systems biology research and empowered applications in structural systems biology and systems pharmacology.

Results

Here, we present the generation, application, and dissemination of genome-scale models with protein structures (GEM-PRO) for Escherichia coli and Thermotoga maritima. We show the utility of integrating molecular scale analyses with systems biology approaches by discussing several comparative analyses on the temperature dependence of growth, the distribution of protein fold families, substrate specificity, and characteristic features of whole cell proteomes. Finally, to aid in the grand challenge of big data to knowledge, we provide several explicit tutorials of how protein-related information can be linked to genome-scale models in a public GitHub repository ( https://github.com/SBRG/GEMPro/tree/master/GEMPro_recon/).

Conclusions

Translating genome-scale, protein-related information to structured data in the format of a GEM provides a direct mapping of gene to gene-product to protein structure to biochemical reaction to network states to phenotypic function. Integration of molecular-level details of individual proteins, such as their physical, chemical, and structural properties, further expands the description of biochemical network-level properties, and can ultimately influence how to model and predict whole cell phenotypes as well as perform comparative systems biology approaches to study differences between organisms. GEM-PRO offers insight into the physical embodiment of an organism's genotype, and its use in this comparative framework enables exploration of adaptive strategies for these organisms, opening the door to many new lines of research. With these provided tools, tutorials, and background, the reader will be in a position to run GEM-PRO for their own purposes.

Unlabelled

BioJava is an open-source project for processing of biological data in the Java programming language. We have recently released a new version (3.0.5), which is a major update to the code base that greatly extends its functionality.

Results

BioJava now consists of several independent modules that provide state-of-the-art tools for protein structure comparison, pairwise and multiple sequence alignments, working with DNA and protein sequences, analysis of amino acid properties, detection of protein modifications and prediction of disordered regions in proteins as well as parsers for common file formats using a biologically meaningful data model.

Availability

BioJava is an open-source project distributed under the Lesser GPL (LGPL). BioJava can be downloaded from the BioJava website (http://www.biojava.org). BioJava requires Java 1.6 or higher. All inquiries should be directed to the BioJava mailing lists. Details are available at http://biojava.org/wiki/BioJava:MailingLists.

In the Florida Panhandle region, bottlenose dolphins (Tursiops truncatus) have been highly susceptible to large-scale unusual mortality events (UMEs) that may have been the result of exposure to blooms of the dinoflagellate Karenia brevis and its neurotoxin, brevetoxin (PbTx). Between 1999 and 2006, three bottlenose dolphin UMEs occurred in the Florida Panhandle region. The primary objective of this study was to determine if these mortality events were due to brevetoxicosis. Analysis of over 850 samples from 105 bottlenose dolphins and associated prey items were analyzed for algal toxins and have provided details on tissue distribution, pathways of trophic transfer, and spatial-temporal trends for each mortality event. In 1999/2000, 152 dolphins died following extensive K. brevis blooms and brevetoxin was detected in 52% of animals tested at concentrations up to 500 ng/g. In 2004, 105 bottlenose dolphins died in the absence of an identifiable K. brevis bloom; however, 100% of the tested animals were positive for brevetoxin at concentrations up to 29,126 ng/mL. Dolphin stomach contents frequently consisted of brevetoxin-contaminated menhaden. In addition, another potentially toxigenic algal species, Pseudo-nitzschia, was present and low levels of the neurotoxin domoic acid (DA) were detected in nearly all tested animals (89%). In 2005/2006, 90 bottlenose dolphins died that were initially coincident with high densities of K. brevis. Most (93%) of the tested animals were positive for brevetoxin at concentrations up to 2,724 ng/mL. No DA was detected in these animals despite the presence of an intense DA-producing Pseudo-nitzschia bloom. In contrast to the absence or very low levels of brevetoxins measured in live dolphins, and those stranding in the absence of a K. brevis bloom, these data, taken together with the absence of any other obvious pathology, provide strong evidence that brevetoxin was the causative agent involved in these bottlenose dolphin mortality events.

CD4+ T cells play a critical role in antitumor immunity via recognition of peptide antigens presented on MHC class II (MHC-II). Although some solid cancers can be induced to express MHC-II, the extent to which this enables direct recognition by tumor-specific CD4+ T cells is unclear. We isolated and characterized T cell antigen receptors (TCRs) from naturally primed CD4+ T cells specific for 2 oncoproteins, HPV-16 E6 and the activating KRASG12V mutation, from patients with head and neck squamous cell carcinoma and pancreatic ductal adenocarcinoma, respectively, and determined their ability to recognize autologous or human leukocyte antigen-matched antigen-expressing tumor cells. We found in both cases that the TCRs were capable of recognizing peptide-loaded target cells expressing the relevant MHC-II or B cell antigen-presenting cells (APCs) when the antigens were endogenously expressed and directed to the endosomal pathway but failed to recognize tumor cells expressing the source protein even after induction of surface MHC-II expression by IFN-γ or transduction with CIITA. These results suggest that priming and functional recognition of both a nuclear (E6) and a membrane-associated (KRAS) oncoprotein are predominantly confined to crosspresenting APCs rather than via direct recognition of tumor cells induced to express MHC-II.

Therapeutic benefit to immune checkpoint blockade (ICB) is currently limited to the subset of cancers thought to possess a sufficient tumor mutational burden (TMB) to allow for the spontaneous recognition of neoantigens (NeoAg) by autologous T cells. We explored whether the response to ICB of an aggressive low-TMB squamous cell tumor could be improved through combination immunotherapy using functionally defined NeoAg as targets for endogenous CD4+ and CD8+ T cells. We found that, whereas vaccination with CD4+ or CD8+ NeoAg alone did not offer prophylactic or therapeutic immunity, vaccines containing NeoAg recognized by both subsets overcame ICB resistance and led to the eradication of large established tumors that contained a subset of PD-L1+ tumor-initiating cancer stem cells (tCSC), provided the relevant epitopes were physically linked. Therapeutic CD4+/CD8+ T cell NeoAg vaccination produced a modified tumor microenvironment (TME) with increased numbers of NeoAg-specific CD8+ T cells existing in progenitor and intermediate exhausted states enabled by combination ICB-mediated intermolecular epitope spreading. We believe that the concepts explored herein should be exploited for the development of more potent personalized cancer vaccines that can expand the range of tumors treatable with ICB.

The central drift chamber for the BaBar detector at the SLAC B-factory is based on a hexagonal cell design with 7104 cells arranged in 40 layers and drift gas Helium:isobutane (80%:20%). Performance optimization and integration requirements led to an electronics design that mounts the amplifier-discriminator and digitizing circuitry directly on the endplate. High channel density is achieved using a 4-channel custom amplifier-discriminator IC and an 8-channel custom CMOS TDC/FADC IC on a single circuit board. Data read from the front ends are multiplexed on 4 fiber optic links, and prompt trigger data are sent out continuously on 24 links. Analysis of cosmic ray data demonstrates that the electronics design meets the performance goals for the BaBar drift chamber. The final electronics were installed on the drift chamber in July, 1998. Installation of BaBar on beamline is scheduled for March, 1999. © 1999 IEEE.