PocketMatch: A new algorithm to compare binding sites in protein structures

Background: Recognizing similarities and deriving relationships among protein molecules is a fundamental
requirement in present-day biology. Similarities can be present at various levels which can be detected through comparison of protein sequences or their structural folds. In some cases similarities obscure at these levels could be present merely in the substructures at their binding sites. Inferring functional similarities between protein molecules by comparing their binding sites is still largely exploratory and not as yet a routine protocol. One of
the main reasons for this is the limitation in the choice of appropriate analytical tools that can compare binding sites with high sensitivity. To benefit from the enormous amount of structural data that is being rapidly accumulated, it is essential to have high throughput tools that enable large scale binding site comparison.

Results: Here we present a new algorithm PocketMatch for comparison of binding sites in a frame invariant
manner. Each binding site is represented by 90 lists of sorted distances capturing shape and chemical nature of the site. The sorted arrays are then aligned using an incremental alignment method and scored to obtain PMScores for pairs of sites. A comprehensive sensitivity analysis and an extensive validation of the algorithm have been carried out. Perturbation studies where the geometry of a given site was retained but the residue types were changed randomly, indicated that chance similarities were virtually non-existent. Our analysis also demonstrates that shape information alone is insufficient to discriminate between diverse binding sites, unless
combined with chemical nature of amino acids.

Conclusions: A new algorithm has been developed to compare binding sites in accurate, efficient and
high-throughput manner. Though the representation used is conceptually simplistic, we demonstrate that along
with the new alignment strategy used, it is sufficient to enable binding comparison with high sensitivity. Novel methodology has also been presented for validating the algorithm for accuracy and sensitivity with respect to geometry and chemical nature of the site. The method is also fast and takes about 1/250th second for one comparison on a single processor. A parallel version on BlueGene has also been implemented

Ecosystem vulnerability to alien and invasive species: a case study on marine habitats along the Italian coast

1. Spread of alien species (AS) is a serious threat to marine habitats and analysis of principal descriptors of their occurrence is pivotal to set reliable conservation strategies. 2. In order to assess the susceptibility of marine habitats to biological invasions, a dataset was gathered of the occurrence of 3899 species from 29 phyla, taken from 93 marine sites located along the Italian coast in the period 2000\u20132012.3. In total, 61 AS belonging to 11 phyla have been recorded. Invertebrates were the most represented (63%). Alien species were found in all the habitats examined (EUNIS, level 2), although they showed highest abundance in benthic habitats. Most of the AS were associated with a single EUNIS habitat, while some of them were present in more than one habitat. Trans-habitat occurrence suggests the potential invasiveness of AS. 4. According to statistical analysis, AS recorded could have been more numerous, since some of the marine habitats seemed to be still unsaturated. The model that best describes the spread of AS takes account of both native species richness (Rn) and EUNIS habitat type as explanatory variables. The number of observed AS was directly related to Rn and it was highest in rocky circalittoral and infralittoral habitats. 5. The results of this macro-ecological study focus on the importance of performing large-scale studies, since adopting ecosystem approaches to marine invasion management seems especially fruitful. 6. The results, moreover, highlight the importance of AS monitoring of different habitats, from those subjected to anthropogenic pressure, historically considered to be hubs of introduction of AS, to the most biologically rich and diverse marine habitats. Indeed, it is necessary to set monitoring strategies to detect the introduction, the distribution and persistence of AS over time. These recommendations are especially significant in the light of the strategic plans currently under formulation in Mediterranean countries with regard to AS monitoring

Genetic markers in the study of Anisakis typica (Diesing, 1860): larval identification and genetic relationships with other species of Anisakis Dujardin, 1845 (Nematoda: Anisakidae)

Genetic variation at 21 gene-enzyme systems was studied in a sample of an adult population of Anisakis typica (Diesing, 1860) recovered in the dolphin Sotalia fluviatilis from the Atlantic coast of Brazil. The characteristic alleles, detected in this population, made it possible to identify as A. typica, Anisakis larvae with a Type I morphol ogy (sensu Berland, 1961) from various fishes: Thunnus thynnus and Auxis thazard from Brazil waters, Trachurus picturatus and Scomber japonicus from Madeiran waters, Scomberomorus commerson, Euthynnus affinis, Sarda orientalis and Coryphaena hippurus from the Somali coast of the Indian Ocean, and Merluccius merluccius from the Eastern Mediterranean. Characteristic allozymes are given for the identification, at any life-stage and in both sexes, of A. typica and the other Anisakis species so far studied genetically. The distribution of A. typica in warmer temperate and tropical waters is confirmed; the definitive hosts so far identified for this species belong to delphinids, phocoenids and pontoporids. The present findings represent the first established records of intermediate/paratenic hosts of A. typica and extend its range to Somali waters of the Indian Ocean and to the Eastern Mediterranean Sea. A remarkable genetic homogeneity was observed in larval and adult samples of A. typica despite their different geographical origin; interpopulation genetic distances were low, ranging from DNei = 0.004 (Eastern Mediterranean versus Somali) to DNei = 0.010 (Brazilian versus Somali). Accordingly, indirect estimates of gene flow gave a rather high average value of Nm = 6.00. Genetic divergence of A. typica was, on average, DNei = 1.12 from the members of the A. simplex complex (A. simplex s.s, A. pegreffii, A. simplex C) and DNei = 1.41 from A. ziphidarum, which all share Type I larvae; higher values were found from both A. physeteris (DNei = 2.77) and A. brevispiculata (DNei = 2.52), which have Type II larvae (sensu Berland, 1961). Genetic relationships among these species are shown using multidimensional scaling ordination (MDS). The genus Anisakis appears to be phylogenetically heterogeneous and includes two distinct groups of species, which are morphologically and genetically differentiated.info:eu-repo/semantics/publishedVersio

Comparative and Functional Genomics of Rhodococcus opacus PD630 for Biofuels Development

The Actinomycetales bacteria Rhodococcus opacus PD630 and Rhodococcus jostii RHA1 bioconvert a diverse range of organic substrates through lipid biosynthesis into large quantities of energy-rich triacylglycerols (TAGs). To describe the genetic basis of the Rhodococcus oleaginous metabolism, we sequenced and performed comparative analysis of the 9.27 Mb R. opacus PD630 genome. Metabolic-reconstruction assigned 2017 enzymatic reactions to the 8632 R. opacus PD630 genes we identified. Of these, 261 genes were implicated in the R. opacus PD630 TAGs cycle by metabolic reconstruction and gene family analysis. Rhodococcus synthesizes uncommon straight-chain odd-carbon fatty acids in high abundance and stores them as TAGs. We have identified these to be pentadecanoic, heptadecanoic, and cis-heptadecenoic acids. To identify bioconversion pathways, we screened R. opacus PD630, R. jostii RHA1, Ralstonia eutropha H16, and C. glutamicum 13032 for growth on 190 compounds. The results of the catabolic screen, phylogenetic analysis of the TAGs cycle enzymes, and metabolic product characterizations were integrated into a working model of prokaryotic oleaginy.Cambridge-MIT InstituteMassachusetts Institute of Technology. (Seed Grant program)Shell Oil CompanyNational Institute of Allergy and Infectious Diseases (U.S.)United States. National Institutes of HealthNational Institutes of Health. Department of Health and Human Services (Contract No. HHSN272200900006C

Research for PECH Committee - Discard ban, landing obligation and MSY in the Western Mediterranean Sea - the Spanish case

The demersal fisheries in the Mediterranean Sea are heavily overfished but the landing obligation will not help to reach MSY because it will not decrease fishing mortality. The new proposal of the Commission introduces total allowable effort as a new way to regulate Western Mediterranean demersal fisheries by significantly reducing fishing time. However, this new management measure must be complemented with increased gear selectivity, implementation of closed areas and local co-management plans. Different approaches to reduce fishing mortality may have different socio-economic impact

High-throughput Screening and Sensitized Bacteria Identify an M. tuberculosis Dihydrofolate Reductase Inhibitor with Whole Cell Activity

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, is a bacterial pathogen that claims roughly 1.4 million lives every year. Current drug regimens are inefficient at clearing infection, requiring at least 6 months of chemotherapy, and resistance to existing agents is rising. There is an urgent need for new drugs that are more effective and faster acting. The folate pathway has been successfully targeted in other pathogens and diseases, but has not yielded a lead drug against tuberculosis. We developed a high-throughput screening assay against Mtb dihydrofolate reductase (DHFR), a critical enzyme in the folate pathway, and screened a library consisting of 32,000 synthetic and natural product-derived compounds. One potent inhibitor containing a quinazoline ring was identified. This compound was active against the wild-type laboratory strain H37Rv (MIC99 = 207 µM). In addition, an Mtb strain with artificially lowered DHFR levels showed increased sensitivity to this compound (MIC99 = 70.7 µM), supporting that the inhibition was target-specific. Our results demonstrate the potential to identify Mtb DHFR inhibitors with activity against whole cells, and indicate the power of using a recombinant strain of Mtb expressing lower levels of DHFR to facilitate the discovery of antimycobacterial agents. With these new tools, we highlight the folate pathway as a potential target for new drugs to combat the tuberculosis epidemic

The essential mycobacterial genes, fabG1 and fabG4, encode 3-oxoacyl-thioester reductases that are functional in yeast mitochondrial fatty acid synthase type 2

Mycobacterium tuberculosis represents a severe threat to human health worldwide. Therefore, it is important to expand our knowledge of vital mycobacterial processes, such as that effected by fatty acid synthase type 2 (FASII), as well as to uncover novel ones. Mycobacterial FASII undertakes mycolic acid biosynthesis, which relies on a set of essential enzymes, including 3-oxoacyl-AcpM reductase FabG1/Rv1483. However, the M. tuberculosis genome encodes four additional FabG homologs, designated FabG2–FabG5, whose functions have hitherto not been characterized in detail. Of the four candidates, FabG4/Rv0242c was recently shown to be essential for the survival of M. bovis BCG. The present work was initiated by assessing the suitability of yeast oar1Δ mutant cells lacking mitochondrial 3-oxoacyl-ACP reductase activity to act as a surrogate system for expressing FabG1/MabA directed to the mitochondria. Mutant yeast cells producing this targeted FabG1 variant were essentially wild type for all of the chronicled phenotype characteristics, including respiratory growth on glycerol medium, cytochrome assembly and lipoid acid production. This indicated that within the framework of de novo fatty acid biosynthesis in yeast mitochondria, FabG1 was able to act on shorter (C4) acyl substrates than was previously proposed (C8–20) during mycolic acid biosynthesis in M. tuberculosis. Thereafter, FabG2–FabG5 were expressed as mitochondrial proteins in the oar1Δ strain, and FabG4 was found to complement the mutant phenotype and contain high levels of 3-oxoacyl-thioester reductase activity. Hence, like FabG1, FabG4 is also an essential, physiologically functional 3-oxoacyl-thioester reductase, albeit the latter’s involvement in mycobacterial FASII remains to be explored

An Anaerobic-Type α-Ketoglutarate Ferredoxin Oxidoreductase Completes the Oxidative Tricarboxylic Acid Cycle of Mycobacterium tuberculosis

Aerobic organisms have a tricarboxylic acid (TCA) cycle that is functionally distinct from those found in anaerobic organisms. Previous reports indicate that the aerobic pathogen Mycobacterium tuberculosis lacks detectable α-ketoglutarate (KG) dehydrogenase activity and drives a variant TCA cycle in which succinyl-CoA is replaced by succinic semialdehyde. Here, we show that M. tuberculosis expresses a CoA-dependent KG dehydrogenase activity, albeit one that is typically found in anaerobic bacteria. Unlike most enzymes of this family, the M. tuberculosis KG: ferredoxin oxidoreductase (KOR) is extremely stable under aerobic conditions. This activity is absent in a mutant strain deleted for genes encoding a previously uncharacterized oxidoreductase, and this strain is impaired for aerobic growth in the absence of sufficient amounts of CO2. Interestingly, inhibition of the glyoxylate shunt or exclusion of exogenous fatty acids alleviates this growth defect, indicating the presence of an alternate pathway that operates in the absence of β-oxidation. Simultaneous disruption of KOR and the first enzyme of the succinic semialdehyde pathway (KG decarboxylase; KGD) results in strict dependence upon the glyoxylate shunt for growth, demonstrating that KG decarboxylase is also functional in M. tuberculosis intermediary metabolism. These observations demonstrate that unlike most organisms M. tuberculosis utilizes two distinct TCA pathways from KG, one that functions concurrently with β-oxidation (KOR-dependent), and one that functions in the absence of β-oxidation (KGD-dependent). As these pathways are regulated by metabolic cues, we predict that their differential utilization provides an advantage for growth in different environments within the host

Genetic Basis of Virulence Attenuation Revealed by Comparative Genomic Analysis of Mycobacterium tuberculosis Strain H37Ra versus H37Rv

Tuberculosis, caused by Mycobacterium tuberculosis, remains a leading infectious disease despite the availability of chemotherapy and BCG vaccine. The commonly used avirulent M. tuberculosis strain H37Ra was derived from virulent strain H37 in 1935 but the basis of virulence attenuation has remained obscure despite numerous studies. We determined the complete genomic sequence of H37Ra ATCC25177 and compared that with its virulent counterpart H37Rv and a clinical isolate CDC1551. The H37Ra genome is highly similar to that of H37Rv with respect to gene content and order but is 8,445 bp larger as a result of 53 insertions and 21 deletions in H37Ra relative to H37Rv. Variations in repetitive sequences such as IS6110 and PE/PPE/PE-PGRS family genes are responsible for most of the gross genetic changes. A total of 198 single nucleotide variations (SNVs) that are different between H37Ra and H37Rv were identified, yet 119 of them are identical between H37Ra and CDC1551 and 3 are due to H37Rv strain variation, leaving only 76 H37Ra-specific SNVs that affect only 32 genes. The biological impact of missense mutations in protein coding sequences was analyzed in silico while nucleotide variations in potential promoter regions of several important genes were verified by quantitative RT-PCR. Mutations affecting transcription factors and/or global metabolic regulations related to in vitro survival under aging stress, and mutations affecting cell envelope, primary metabolism, in vivo growth as well as variations in the PE/PPE/PE-PGRS family genes, may underlie the basis of virulence attenuation. These findings have implications not only for improved understanding of pathogenesis of M. tuberculosis but also for development of new vaccines and new therapeutic agents