QATS: an ImageJ plugin for the quantification of toroidal nuclei in biological images

MOTIVATION

The toroidal nucleus is a novel chromosomal instability (CIN) biomarker which complements the micronucleus. Understanding the specific biological stresses leading to the formation of each CIN-associated phenotype requires the evaluation of large panels of biological images collected from different genetic backgrounds and environmental conditions. However, the quantification of toroidal nuclei is currently a manual process which is unviable on a large scale.

RESULTS

Here, we present QATS (QuAntification of Toroidal nuclei in biological imageS), a tool that automates the identification of toroidal nuclei, minimizing false positives while highly agreeing with the manual quantifications. Additionally, QATS identifies micronuclei for a convenient comparison of both CIN biomarkers. QATS is an open-source ImageJ plugin with a user-friendly interface that enables a wide scientific community to easily assess the frequency of CIN biomarkers for the determination of CIN levels in cellular models.

AVAILABILITY AND IMPLEMENTATION

QATS is an ImageJ plugin freely available at http://www.toroidalnucleus.org/qats. The user manual and the images used for the evaluation of QATS are included in the website. Supplementary data are available at Bioinformatics online.

Meta-analysis of dispensable essential genes and their interactions with bypass suppressors

Genes have been historically classified as essential or non-essential based on their requirement for viability. However, genomic mutations can sometimes bypass the requirement for an essential gene, challenging the binary classification of gene essentiality. Such dispensable essential genes represent a valuable model for understanding the incomplete penetrance of loss-of-function mutations often observed in natural populations. Here, we compiled data from multiple studies on essential gene dispensability in Saccharomyces cerevisiae to comprehensively characterize these genes. In analyses spanning different evolutionary timescales, dispensable essential genes exhibited distinct phylogenetic properties compared with other essential and non-essential genes. Integration of interactions with suppressor genes that can bypass the gene essentiality revealed the high functional modularity of the bypass suppression network. Furthermore, dispensable essential and bypass suppressor gene pairs reflected simultaneous changes in the mutational landscape of S. cerevisiae strains. Importantly, species in which dispensable essential genes were non-essential tended to carry bypass suppressor mutations in their genomes. Overall, our study offers a comprehensive view of dispensable essential genes and illustrates how their interactions with bypass suppressors reflect evolutionary outcomes.

Global analysis of suppressor mutations that rescue human genetic defects

BACKGROUND

Genetic suppression occurs when the deleterious effects of a primary "query" mutation, such as a disease-causing mutation, are rescued by a suppressor mutation elsewhere in the genome.

METHODS

To capture existing knowledge on suppression relationships between human genes, we examined 2,400 published papers for potential interactions identified through either genetic modification of cultured human cells or through association studies in patients.

RESULTS

The resulting network encompassed 476 unique suppression interactions covering a wide spectrum of diseases and biological functions. The interactions frequently linked genes that operate in the same biological process. Suppressors were strongly enriched for genes with a role in stress response or signaling, suggesting that deleterious mutations can often be buffered by modulating signaling cascades or immune responses. Suppressor mutations tended to be deleterious when they occurred in absence of the query mutation, in apparent contrast with their protective role in the presence of the query. We formulated and quantified mechanisms of genetic suppression that could explain 71% of interactions and provided mechanistic insight into disease pathology. Finally, we used these observations to predict suppressor genes in the human genome.

CONCLUSIONS

The global suppression network allowed us to define principles of genetic suppression that were conserved across diseases, model systems, and species. The emerging frequency of suppression interactions among human genes and range of underlying mechanisms, together with the prevalence of suppression in model organisms, suggest that compensatory mutations may exist for most genetic diseases.

K29-linked unanchored polyubiquitin chains disrupt ribosome biogenesis and direct ribosomal proteins to the Intranuclear Quality control compartment (INQ)

Ribosome assembly requires precise coordination between the production and assembly of ribosomal components. Mutations in ribosomal proteins that inhibit the assembly process or ribosome function are often associated with Ribosomopathies, some of which are linked to defects in proteostasis. In this study, we examine the interplay between several yeast proteostasis enzymes, including deubiquitylases (DUBs), Ubp2 and Ubp14, and E3 ligases, Ufd4 and Hul5, and we explore their roles in the regulation of the cellular levels of K29-linked unanchored polyubiquitin (polyUb) chains. Accumulating K29-linked unanchored polyUb chains associate with maturing ribosomes to disrupt their assembly, activate the Ribosome assembly stress response (RASTR), and lead to the sequestration of ribosomal proteins at the Intranuclear Quality control compartment (INQ). These findings reveal the physiological relevance of INQ and provide insights into mechanisms of cellular toxicity associated with Ribosomopathies.

A proteome-scale map of the SARS-CoV-2-human contactome

Understanding the mechanisms of coronavirus disease 2019 (COVID-19) disease severity to efficiently design therapies for emerging virus variants remains an urgent challenge of the ongoing pandemic. Infection and immune reactions are mediated by direct contacts between viral molecules and the host proteome, and the vast majority of these virus-host contacts (the 'contactome') have not been identified. Here, we present a systematic contactome map of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with the human host encompassing more than 200 binary virus-host and intraviral protein-protein interactions. We find that host proteins genetically associated with comorbidities of severe illness and long COVID are enriched in SARS-CoV-2 targeted network communities. Evaluating contactome-derived hypotheses, we demonstrate that viral NSP14 activates nuclear factor κB (NF-κB)-dependent transcription, even in the presence of cytokine signaling. Moreover, for several tested host proteins, genetic knock-down substantially reduces viral replication. Additionally, we show for USP25 that this effect is phenocopied by the small-molecule inhibitor AZ1. Our results connect viral proteins to human genetic architecture for COVID-19 severity and offer potential therapeutic targets.

Detection of Nuclear Biomarkers for Chromosomal Instability

Chromosomal instability (CIN) is a hallmark of cancer, which is characterized by the gain or loss of chromosomes as well as the rearrangement of the genetic material during cell division. Detection of mitotic errors such as misaligned chromosomes or chromosomal bridges (also known as lagging chromosomes) is challenging as it requires the analysis and manual discrimination of chromosomal aberrations in mitotic cells by molecular techniques. In interphase cells, more frequent in the cell population than mitotic cells, two distinct nuclear phenotypes are associated with CIN: the micronucleus and the toroidal nucleus. Several methods are available for the detection of micronuclei, but none for toroidal nuclei. Here, we provide a method to quantify the presence of both nuclear biomarkers for the evaluation of CIN status in non-mitotic cells particularly suited for genotoxicity screens.

Systematic analysis of bypass suppression of essential genes

Essential genes tend to be highly conserved across eukaryotes, but, in some cases, their critical roles can be bypassed through genetic rewiring. From a systematic analysis of 728 different essential yeast genes, we discovered that 124 (17%) were dispensable essential genes. Through whole-genome sequencing and detailed genetic analysis, we investigated the genetic interactions and genome alterations underlying bypass suppression. Dispensable essential genes often had paralogs, were enriched for genes encoding membrane-associated proteins, and were depleted for members of protein complexes. Functionally related genes frequently drove the bypass suppression interactions. These gene properties were predictive of essential gene dispensability and of specific suppressors among hundreds of genes on aneuploid chromosomes. Our findings identify yeast's core essential gene set and reveal that the properties of dispensable essential genes are conserved from yeast to human cells, correlating with human genes that display cell line-specific essentiality in the Cancer Dependency Map (DepMap) project.

Natural variants suppress mutations in hundreds of essential genes

The consequence of a mutation can be influenced by the context in which it operates. For example, loss of gene function may be tolerated in one genetic background, and lethal in another. The extent to which mutant phenotypes are malleable, the architecture of modifiers and the identities of causal genes remain largely unknown. Here, we measure the fitness effects of ~ 1,100 temperature‐sensitive alleles of yeast essential genes in the context of variation from ten different natural genetic backgrounds and map the modifiers for 19 combinations. Altogether, fitness defects for 149 of the 580 tested genes (26%) could be suppressed by genetic variation in at least one yeast strain. Suppression was generally driven by gain‐of‐function of a single, strong modifier gene, and involved both genes encoding complex or pathway partners suppressing specific temperature‐sensitive alleles, as well as general modifiers altering the effect of many alleles. The emerging frequency of suppression and range of possible mechanisms suggest that a substantial fraction of monogenic diseases could be managed by modulating other gene products.

Environmental robustness of the global yeast genetic interaction network

Phenotypes associated with genetic variants can be altered by interactions with other genetic variants (GxG), with the environment (GxE), or both (GxGxE). Yeast genetic interactions have been mapped on a global scale, but the environmental influence on the plasticity of genetic networks has not been examined systematically. To assess environmental rewiring of genetic networks, we examined 14 diverse conditions and scored 30,000 functionally representative yeast gene pairs for dynamic, differential interactions. Different conditions revealed novel differential interactions, which often uncovered functional connections between distantly related gene pairs. However, the majority of observed genetic interactions remained unchanged in different conditions, suggesting that the global yeast genetic interaction network is robust to environmental perturbation and captures the fundamental functional architecture of a eukaryotic cell.

Analysing the yeast complexome-the Complex Portal rising to the challenge

The EMBL-EBI Complex Portal is a knowledgebase of macromolecular complexes providing persistent stable identifiers. Entries are linked to literature evidence and provide details of complex membership, function, structure and complex-specific Gene Ontology annotations. Data are freely available and downloadable in HUPO-PSI community standards and missing entries can be requested for curation. In collaboration with Saccharomyces Genome Database and UniProt, the yeast complexome, a compendium of all known heteromeric assemblies from the model organism Saccharomyces cerevisiae, was curated. This expansion of knowledge and scope has led to a 50% increase in curated complexes compared to the previously published dataset, CYC2008. The yeast complexome is used as a reference resource for the analysis of complexes from large-scale experiments. Our analysis showed that genes coding for proteins in complexes tend to have more genetic interactions, are co-expressed with more genes, are more multifunctional, localize more often in the nucleus, and are more often involved in nucleic acid-related metabolic processes and processes where large machineries are the predominant functional drivers. A comparison to genetic interactions showed that about 40% of expanded co-complex pairs also have genetic interactions, suggesting strong functional links between complex members.

Pulvinar and pubic cartilage measurements to refine universal ultrasound screening for developmental dysplasia of the hip: Data from 1896 infant hips

PURPOSE

Screening modalities for Developmental Dysplasia of the Hip (DDH) and indications for treatment of mild forms remain controversial. Ultrasound (US) measurement of the pubofemoral distance (PFD > 6 mm, composed of the pubic cartilage and the pulvinar) can avoid late diagnoses of DDH. A thick pubic cartilage may nevertheless lead to false positives. The purpose of this study was to establish standard measurements of pubic cartilage and pulvinar, through universal US screening, to lower false positive results and thus any overtreatment.

METHODS

This is a single-center observational prospective study conducted from December 2016 to January 2018, on infants who underwent universal US screening for DDH. The only inclusion criterion was an adjusted age between 4 and 12 weeks when US was realized. PFD measurement was made using the Couture and Tréguier method. In addition, thicknesses of pubic cartilage and pulvinar were measured on the same US section, in millimeters.

RESULTS

Nine hundred and forty-eight patients, representing 1896 hips, were included. The average value of pubic cartilage thickness was 1.25 mm ± 0.58 mm, with an upper threshold of 2.39 mm (+1.96σ). The average value of pulvinar thickness was 2.67 mm ± 0.78 mm, with an upper threshold of 4.20 mm (+1.96σ). We found high inter-observer reproducibility in pubic cartilage measurements.

CONCLUSION

Systematic measurements of pubic cartilage and pulvinar may refine therapeutic decision by identifying false positives. Patients with increased PFD due to a thick pubic cartilage >2,39 mm, without an associated pulvinar enlargement (<4,20 mm), could be therefore only monitored and not overtreated.

Bioactivity Profile Similarities to Expand the Repertoire of COVID-19 Drugs

Until a vaccine becomes available, the current repertoire of drugs is our only therapeutic asset to fight the SARS-CoV-2 outbreak. Indeed, emergency clinical trials have been launched to assess the effectiveness of many marketed drugs, tackling the decrease of viral load through several mechanisms. Here, we present an online resource, based on small-molecule bioactivity signatures and natural language processing, to expand the portfolio of compounds with potential to treat COVID-19. By comparing the set of drugs reported to be potentially active against SARS-CoV-2 to a universe of 1 million bioactive molecules, we identify compounds that display analogous chemical and functional features to the current COVID-19 candidates. Searches can be filtered by level of evidence and mechanism of action, and results can be restricted to drug molecules or include the much broader space of bioactive compounds. Moreover, we allow users to contribute COVID-19 drug candidates, which are automatically incorporated to the pipeline once per day. The computational platform, as well as the source code, is available at https://sbnb.irbbarcelona.org/covid19.

A reference map of the human binary protein interactome

Global insights into cellular organization and genome function require comprehensive understanding of the interactome networks that mediate genotype-phenotype relationships^1,2. Here, we present a human “all-by-all” reference interactome map of human binary protein interactions, or “HuRI”. With ~53,000 high-quality protein-protein interactions (PPIs), HuRI has approximately four times more such interactions than high-quality curated interactions from small-scale studies. Integrating HuRI with genome³, transcriptome⁴, and proteome⁵ data enables the study of cellular function within most physiological or pathological cellular contexts. We demonstrate the utility of HuRI in identifying specific subcellular roles of PPIs. Inferred tissue-specific networks reveal general principles for the formation of cellular context-specific functions and elucidate potential molecular mechanisms underlying tissue-specific phenotypes of Mendelian diseases. HuRI represents a systematic proteome-wide reference linking genomic variation to phenotypic outcomes.

Systematic genetics and single-cell imaging reveal widespread morphological pleiotropy and cell-to-cell variability

Our ability to understand the genotype-to-phenotype relationship is hindered by the lack of detailed understanding of phenotypes at a single-cell level. To systematically assess cell-to-cell phenotypic variability, we combined automated yeast genetics, high-content screening and neural network-based image analysis of single cells, focussing on genes that influence the architecture of four subcellular compartments of the endocytic pathway as a model system. Our unbiased assessment of the morphology of these compartments-endocytic patch, actin patch, late endosome and vacuole-identified 17 distinct mutant phenotypes associated with ~1,600 genes (~30% of all yeast genes). Approximately half of these mutants exhibited multiple phenotypes, highlighting the extent of morphological pleiotropy. Quantitative analysis also revealed that incomplete penetrance was prevalent, with the majority of mutants exhibiting substantial variability in phenotype at the single-cell level. Our single-cell analysis enabled exploration of factors that contribute to incomplete penetrance and cellular heterogeneity, including replicative age, organelle inheritance and response to stress.

A PanorOmic view of personal cancer genomes

The massive molecular profiling of thousands of cancer patients has led to the identification of many tumor type specific driver genes. However, only a few (or none) of them are present in each individual tumor and, to enable precision oncology, we need to interpret the alterations found in a single patient. Cancer PanorOmics (http://panoromics.irbbarcelona.org) is a web-based resource to contextualize genomic variations detected in a personal cancer genome within the body of clinical and scientific evidence available for 26 tumor types, offering complementary cohort- and patient-centric views. Additionally, it explores the cellular environment of mutations by mapping them on the human interactome and providing quasi-atomic structural details, whenever available. This ‘PanorOmic’ molecular view of individual tumors, together with the appropriate genetic counselling and medical advice, should contribute to the identification of actionable alterations ultimately guiding the clinical decision-making process.

Features of the Chaperone Cellular Network Revealed through Systematic Interaction Mapping

A comprehensive view of molecular chaperone function in the cell was obtained through a systematic global integrative network approach based on physical (protein-protein) and genetic (gene-gene or epistatic) interaction mapping. This allowed us to decipher interactions involving all core chaperones (67) and cochaperones (15) of Saccharomyces cerevisiae. Our analysis revealed the presence of a large chaperone functional supercomplex, which we named the naturally joined (NAJ) chaperone complex, encompassing Hsp40, Hsp70, Hsp90, AAA+, CCT, and small Hsps. We further found that many chaperones interact with proteins that form foci or condensates under stress conditions. Using an in vitro reconstitution approach, we demonstrate condensate formation for the highly conserved AAA+ ATPases Rvb1 and Rvb2, which are part of the R2TP complex that interacts with Hsp90. This expanded view of the chaperone network in the cell clearly demonstrates the distinction between chaperones having broad versus narrow substrate specificities in protein homeostasis.

Systematic analysis of complex genetic interactions

To systematically explore complex genetic interactions, we constructed ~200,000 yeast triple mutants and scored negative trigenic interactions. We selected double-mutant query genes across a broad spectrum of biological processes, spanning a range of quantitative features of the global digenic interaction network and tested for a genetic interaction with a third mutation. Trigenic interactions often occurred among functionally related genes, and essential genes were hubs on the trigenic network. Despite their functional enrichment, trigenic interactions tended to link genes in distant bioprocesses and displayed a weaker magnitude than digenic interactions. We estimate that the global trigenic interaction network is ~100 times as large as the global digenic network, highlighting the potential for complex genetic interactions to affect the biology of inheritance, including the genotype-to-phenotype relationship.

A framework for exhaustively mapping functional missense variants

Although we now routinely sequence human genomes, we can confidently identify only a fraction of the sequence variants that have a functional impact. Here, we developed a deep mutational scanning framework that produces exhaustive maps for human missense variants by combining random codon mutagenesis and multiplexed functional variation assays with computational imputation and refinement. We applied this framework to four proteins corresponding to six human genes: UBE2I (encoding SUMO E2 conjugase), SUMO1 (small ubiquitin‐like modifier), TPK1 (thiamin pyrophosphokinase), and CALM1/2/3 (three genes encoding the protein calmodulin). The resulting maps recapitulate known protein features and confidently identify pathogenic variation. Assays potentially amenable to deep mutational scanning are already available for 57% of human disease genes, suggesting that DMS could ultimately map functional variation for all human disease genes.

Rate of adverse events of gastroduodenal snare polypectomy for non-flat polyp is low: A prospective and multicenter study

AIM

To evaluate the rate of adverse events (AEs) during consecutive gastric and duodenal polypectomies in several Spanish centers.

METHODS

Polypectomies of protruded gastric or duodenal polyps ≥ 5 mm using hot snare were prospectively included. Prophylactic measures of hemorrhage were allowed in predefined cases. AEs were defined and graded according to the lexicon recommended by the American Society for Gastrointestinal Endoscopy. Patients were followed for 48 h, one week and 1 mo after the procedure.

RESULTS

308 patients were included and a single polypectomy was performed in 205. Only 36 (11.7%) were on prior anticoagulant therapy. Mean polyp size was 15 ± 8.9 mm (5-60) and in 294 cases (95.4%) were located in the stomach. Hemorrhage prophylaxis was performed in 219 (71.1%) patients. Nine patients presented AEs (2.9%), and 6 of them were bleeding (n = 6, 1.9%) (in 5 out of 6 AE, different types of endoscopic treatment were performed). Other 24 hemorrhagic episodes could be managed without any change in the outcome of the endoscopy and, consequently, were considered incidents. We did not find any independent risk factor of bleeding.

CONCLUSION

Gastroduodenal polypectomy using prophylactic measures has a rate of AEs small enough to consider this procedure a safe and effective method for polyp resection independently of the polyp size and location.

Mechanisms of suppression: The wiring of genetic resilience

Recent analysis of genome sequences has identified individuals that are healthy despite carrying severe disease-associated mutations. A possible explanation is that these individuals carry a second genomic perturbation that can compensate for the detrimental effects of the disease allele, a phenomenon referred to as suppression. In model organisms, suppression interactions are generally divided into two classes: genomic suppressors which are secondary mutations in the genome that bypass a mutant phenotype, and dosage suppression interactions in which overexpression of a suppressor gene rescues a mutant phenotype. Here, we describe the general properties of genomic and dosage suppression, with an emphasis on the budding yeast. We propose that suppression interactions between genetic variants are likely relevant for determining the penetrance of human traits. Consequently, an understanding of suppression mechanisms may guide the discovery of protective variants in healthy individuals that carry disease alleles, which could direct the rational design of new therapeutics.

Exploring genetic suppression interactions on a global scale

Genetic suppression occurs when the phenotypic defects caused by a mutation in a particular gene are rescued by a mutation in a second gene. To explore the principles of genetic suppression, we examined both literature-curated and unbiased experimental data, involving systematic genetic mapping and whole-genome sequencing, to generate a large-scale suppression network among yeast genes. Most suppression pairs identified novel relationships among functionally related genes, providing new insights into the functional wiring diagram of the cell. In addition to suppressor mutations, we identified frequent secondary mutations,in a subset of genes, that likely cause a delay in the onset of stationary phase, which appears to promote their enrichment within a propagating population. These findings allow us to formulate and quantify general mechanisms of genetic suppression.

A global genetic interaction network maps a wiring diagram of cellular function

We generated a global genetic interaction network for Saccharomyces cerevisiae, constructing more than 23 million double mutants, identifying about 550,000 negative and about 350,000 positive genetic interactions. This comprehensive network maps genetic interactions for essential gene pairs, highlighting essential genes as densely connected hubs. Genetic interaction profiles enabled assembly of a hierarchical model of cell function, including modules corresponding to protein complexes and pathways, biological processes, and cellular compartments. Negative interactions connected functionally related genes, mapped core bioprocesses, and identified pleiotropic genes, whereas positive interactions often mapped general regulatory connections among gene pairs, rather than shared functionality. The global network illustrates how coherent sets of genetic interactions connect protein complex and pathway modules to map a functional wiring diagram of the cell.

A gene-centered C. elegans protein-DNA interaction network provides a framework for functional predictions

Transcription factors (TFs) play a central role in controlling spatiotemporal gene expression and the response to environmental cues. A comprehensive understanding of gene regulation requires integrating physical protein–DNA interactions (PDIs) with TF regulatory activity, expression patterns, and phenotypic data. Although great progress has been made in mapping PDIs using chromatin immunoprecipitation, these studies have only characterized ~10% of TFs in any metazoan species. The nematode C. elegans has been widely used to study gene regulation due to its compact genome with short regulatory sequences. Here, we delineated the largest gene‐centered metazoan PDI network to date by examining interactions between 90% of C. elegans TFs and 15% of gene promoters. We used this network as a backbone to predict TF binding sites for 77 TFs, two‐thirds of which are novel, as well as integrate gene expression, protein–protein interaction, and phenotypic data to predict regulatory and biological functions for multiple genes and TFs.

Comparison of two Dutch follow-up care models for spinal cord-injured patients and their impact on health problems, re-admissions and quality of care

Objective: To evaluate whether transmural care for people with spinal cord injury living in the community has more impact on health outcomes than traditional follow-up care within the Netherlands.

Design: Quasi-experiment with 12 months of follow-up.

Setting: Eight Dutch rehabilitation centres.

Subjects: Thirty-one patients who received transmural care in two `experimental' rehabilitation centres were compared with a matched sample of 31 patients having received `usual follow-up care' in six other rehabilitation centres.

Intervention: The core component of the transmural care consists of a transmural nurse, who `liaises' between former patients living in the community, primary care professionals and the rehabilitation team. The transmural care model provides activities to support patients and their family/partners and activities to promote continuity of care.

Main measures: The prevalence of pressure sores and urinary tract infections; the number and duration of re-admissions to hospital and rehabilitation centre due to pressure sores, bladder and bowel problems; and the experienced quality of follow-up care.

Results: The transmural care, as implemented, did not influence the health outcomes. The prevalence of pressure sores, urinary tract infections and the number of re-admissions (due to pressure sores, bladder and bowel problems) was respectively 13, 13 and 4 in the intervention group versus 14, 15 and 6 in the usual follow-up care group. Since the transmural care had been incompletely implemented and there were methodological and practical limitations, we formulated no final conclusions regarding its effectiveness.

Conclusion: Implementing the transmural care model strictly according to protocol may improve its effectiveness.

Genetic Interactions Implicating Postreplicative Repair in Okazaki Fragment Processing

Ubiquitination of the replication clamp proliferating cell nuclear antigen (PCNA) at the conserved residue lysine (K)164 triggers postreplicative repair (PRR) to fill single-stranded gaps that result from stalled DNA polymerases. However, it has remained elusive as to whether cells engage PRR in response to replication defects that do not directly impair DNA synthesis. To experimentally address this question, we performed synthetic genetic array (SGA) analysis with a ubiquitination-deficient K164 to arginine (K164R) mutant of PCNA against a library of S. cerevisiae temperature-sensitive alleles. The SGA signature of the K164R allele showed a striking correlation with profiles of mutants deficient in various aspects of lagging strand replication, including rad27Δ and elg1Δ. Rad27 is the primary flap endonuclease that processes 5' flaps generated during lagging strand replication, whereas Elg1 has been implicated in unloading PCNA from chromatin. We observed chronic ubiquitination of PCNA at K164 in both rad27Δ and elg1Δ mutants. Notably, only rad27Δ cells exhibited a decline in cell viability upon elimination of PRR pathways, whereas elg1Δ mutants were not affected. We further provide evidence that K164 ubiquitination suppresses replication stress resulting from defective flap processing during Okazaki fragment maturation. Accordingly, ablation of PCNA ubiquitination increased S phase checkpoint activation, indicated by hyperphosphorylation of the Rad53 kinase. Furthermore, we demonstrate that alternative flap processing by overexpression of catalytically active exonuclease 1 eliminates PCNA ubiquitination. This suggests a model in which unprocessed flaps may directly participate in PRR signaling. Our findings demonstrate that PCNA ubiquitination at K164 in response to replication stress is not limited to DNA synthesis defects but extends to DNA processing during lagging strand replication.

Transcription Factor Activity Mapping of a Tissue-Specific in vivo Gene Regulatory Network

A wealth of physical interaction data between transcription factors (TFs) and DNA has been generated, but these interactions often do not have apparent regulatory consequences. Thus, equating physical interaction data with gene regulatory networks (GRNs) is problematic. Here, we comprehensively assay TF activity, rather than binding, to construct a network of gene regulatory interactions in the C. elegans intestine. By manually observing the in vivo tissue-specific knockdown of 921 TFs on a panel of 19 fluorescent transcriptional reporters, we identified a GRN of 411 interactions between 19 promoters and 177 TFs. This GRN shows only modest overlap with physical interactions, indicating that many regulatory interactions are indirect. We applied nested effects modeling to uncover information flow between TFs in the intestine that converges on a small set of physical TF-promoter interactions. We found numerous cell nonautonomous regulatory interactions, illustrating tissue-to-tissue communication. Altogether, our study illuminates the complexity of gene regulation in the context of a living animal.

pyDockSAXS: protein-protein complex structure by SAXS and computational docking

Structural characterization of protein–protein interactions at molecular level is essential to understand biological processes and identify new therapeutic opportunities. However, atomic resolution structural techniques cannot keep pace with current advances in interactomics. Low-resolution structural techniques, such as small-angle X-ray scattering (SAXS), can be applied at larger scale, but they miss atomic details. For efficient application to protein–protein complexes, low-resolution information can be combined with theoretical methods that provide energetic description and atomic details of the interactions. Here we present the pyDockSAXS web server (http://life.bsc.es/pid/pydocksaxs) that provides an automatic pipeline for modeling the structure of a protein–protein complex from SAXS data. The method uses FTDOCK to generate rigid-body docking models that are subsequently evaluated by a combination of pyDock energy-based scoring function and their capacity to describe SAXS data. The only required input files are structural models for the interacting partners and a SAXS curve. The server automatically provides a series of structural models for the complex, sorted by the pyDockSAXS scoring function. The user can also upload a previously computed set of docking poses, which opens the possibility to filter the docking solutions by potential interface residues or symmetry restraints. The server is freely available to all users without restriction.

Almagate interference in breath test results for the diagnosis of Helicobacter pylori infection

BACKGROUND

Infection by Helicobacter pylori is common and affects both genders at any age. The 13C-urea breath test is a widely used test for the diagnosis of this infection. However, multiple drugs used for the treatment of Helicobacter pylori infection symptoms have interactions with this breath test that generate false negative results. This observational study was to assess the potential interaction between almagate and the breath test.

METHODS

Thirty subjects on almagate therapy who underwent a breath test were included. If the result was negative, almagate was withdrawn for a month and the breath test was then repeated.

RESULTS

In general, 51.9 % of assessed subjects had a negative result after the first test, and 100 % of these also had a negative result after the second test.

CONCLUSIONS

It was concluded that the use of almagate does not interfere in breath test results. These results provide a drug therapy option for the treatment of symptoms associated with Helicobacter pylori infection during the diagnostic process.

Expanding the frontiers of protein-protein modeling: from docking and scoring to binding affinity predictions and other challenges

In addition to protein-protein docking, this CAPRI edition included new challenges, like protein-water and protein-sugar interactions, or the prediction of binding affinities and ΔΔG changes upon mutation. Regarding the standard protein-protein docking cases, our approach, mostly based on the pyDock scheme, submitted correct models as predictors and as scorers for 67% and 57% of the evaluated targets, respectively. In this edition, available information on known interface residues hardly made any difference for our predictions. In one of the targets, the inclusion of available experimental small-angle X-ray scattering (SAXS) data using our pyDockSAXS approach slightly improved the predictions. In addition to the standard protein-protein docking assessment, new challenges were proposed. One of the new problems was predicting the position of the interface water molecules, for which we submitted models with 20% and 43% of the water-mediated native contacts predicted as predictors and scorers, respectively. Another new problem was the prediction of protein-carbohydrate binding, where our submitted model was very close to being acceptable. A set of targets were related to the prediction of binding affinities, in which our pyDock scheme was able to discriminate between natural and designed complexes with area under the curve = 83%. It was also proposed to estimate the effect of point mutations on binding affinity. Our approach, based on machine learning methods, showed high rates of correctly classified mutations for all cases. The overall results were highly rewarding, and show that the field is ready to move forward and face new interesting challenges in interactomics.

Extensive rewiring and complex evolutionary dynamics in a C. elegans multiparameter transcription factor network

Gene duplication results in two identical paralogs that diverge through mutation, leading to loss or gain of interactions with other biomolecules. Here, we comprehensively characterize such network rewiring for C. elegans transcription factors (TFs) within and across four newly delineated molecular networks. Remarkably, we find that even highly similar TFs often have different interaction degrees and partners. In addition, we find that most TF families have a member that is highly connected in multiple networks. Further, different TF families have opposing correlations between network connectivity and phylogenetic age, suggesting that they are subject to different evolutionary pressures. Finally, TFs that have similar partners in one network generally do not in another, indicating a lack of pressure to retain cross-network similarity. Our multiparameter analyses provide unique insights into the evolutionary dynamics that shaped TF networks.

pyDockWEB: a web server for rigid-body protein-protein docking using electrostatics and desolvation scoring

pyDockWEB is a web server for the rigid-body docking prediction of protein-protein complex structures using a new version of the pyDock scoring algorithm. We use here a new custom parallel FTDock implementation, with adjusted grid size for optimal FFT calculations, and a new version of pyDock, which dramatically speeds up calculations while keeping the same predictive accuracy. Given the 3D coordinates of two interacting proteins, pyDockWEB returns the best docking orientations as scored mainly by electrostatics and desolvation energy.

AVAILABILITY AND IMPLEMENTATION

The server does not require registration by the user and is freely accessible for academics at http://life.bsc.es/servlet/pydock.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Efficient Relaxation of Protein-Protein Interfaces by Discrete Molecular Dynamics Simulations

Protein-protein interactions are responsible for the transfer of information inside the cell and represent one of the most interesting research fields in structural biology. Unfortunately, after decades of intense research, experimental approaches still have difficulties in providing 3D structures for the hundreds of thousands of interactions formed between the different proteins in a living organism. The use of theoretical approaches like docking aims to complement experimental efforts to represent the structure of the protein interactome. However, we cannot ignore that current methods have limitations due to problems of sampling of the protein-protein conformational space and the lack of accuracy of available force fields. Cases that are especially difficult for prediction are those in which complex formation implies a non-negligible change in the conformation of the interacting proteins, i.e., those cases where protein flexibility plays a key role in protein-protein docking. In this work, we present a new approach to treat flexibility in docking by global structural relaxation based on ultrafast discrete molecular dynamics. On a standard benchmark of protein complexes, the method provides a general improvement over the results obtained by rigid docking. The method is especially efficient in cases with large conformational changes upon binding, in which structure relaxation with discrete molecular dynamics leads to a predictive success rate double that obtained with state-of-the-art rigid-body docking.

Cell-Dock: high-performance protein-protein docking

SUMMARY

The application of docking to large-scale experiments or the explicit treatment of protein flexibility are part of the new challenges in structural bioinformatics that will require large computer resources and more efficient algorithms. Highly optimized fast Fourier transform (FFT) approaches are broadly used in docking programs but their optimal code implementation leaves hardware acceleration as the only option to significantly reduce the computational cost of these tools. In this work we present Cell-Dock, an FFT-based docking algorithm adapted to the Cell BE processor. We show that Cell-Dock runs faster than FTDock with maximum speedups of above 200×, while achieving results of similar quality.

AVAILABILITY AND IMPLEMENTATION

The source code is released under GNU General Public License version 2 and can be downloaded from http://mmb.pcb.ub.es/~cpons/Cell-Dock.

CONTACT

djimenez@ac.upc.edu or juanf@bsc.es

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Effects of nucleotides adenosine monophosphate and adenosine triphosphate in combination with L-arginine on male rabbit corpus cavernosum tissue

Purines and more specifically adenosine monophosphate (AMP) and adenosine triphosphate (ATP) have a strong relaxant effect on smooth muscle cells of the dog, rabbit and human corpus cavernosum, to approximately the same degree as nitric oxide (NO). However, purines are considered as modulators of erectile function rather than key mediators. This suggests that the use of purines combined with NO donors could be effective to treat some specific erectile disorders. The relaxation induced by the combination of l-arginine (Arg), a natural substrate for NO synthase, was assessed with a purine-nucleotide (AMP, ATP) on a rabbit corpus cavernosum model, to determine if these substances could potentiate each other's effect. When a pre-contraction was induced by phenylephrine, AMP alone induced a 43% CC relaxation rate and ATP alone a 26% rate. The relaxation rate induced by Arg was lower in comparison (8% at 5.10(-4) m vs. 25% at AMP 5.10(-4) m and 15% at ATP 5.10(-4) m). NO synthase inhibitor n-nitro-l-arginine did not modify the relaxing effect provoked by AMP suggesting that the mechanism of action of this nucleotide does not involve the NO pathway. The combination of Arg at 5.10(-4) m with either AMP or ATP at different doses ranging from 5.10(-4) to 10(-3) m significantly enhanced the relaxing response reaching rates of 62 and 80% respectively, leading to a synergistic effect. The present data indicate that a 'NO donor' combined with an 'adenosine donor' could be an effective therapeutic approach.

Mummified fetus in the thoracic cavity of a domestic short-haired cat

This short communication describes the diagnosis, treatment, and clinical course of a domestic short-haired cat with diaphragmatic hernia in which the herniated structure in the thoracic cavity contained a mummified fetus. The cat was pregnant when rescued from the street and, days later, gave birth without abnormalities. Some months later, during an ovariohysterectomy, an abnormal localization of the uterus was observed, and at that time the case was referred to our center. A thoracic radiograph showed an abnormal thoracic mass cranial to the heart. The main suspicion was the presence of a thoracic hernia with the uterus herniated and containing a mummified fetus. A thoracotomy was performed to confirm the nature of the mass and reduce the diaphragmatic hernia. Although this clinical case is quite rare, a mummified fetus can be observed in thoracic hernias.

Community-wide assessment of protein-interface modeling suggests improvements to design methodology

The CAPRI (Critical Assessment of Predicted Interactions) and CASP (Critical Assessment of protein Structure Prediction) experiments have demonstrated the power of community-wide tests of methodology in assessing the current state of the art and spurring progress in the very challenging areas of protein docking and structure prediction. We sought to bring the power of community-wide experiments to bear on a very challenging protein design problem that provides a complementary but equally fundamental test of current understanding of protein-binding thermodynamics. We have generated a number of designed protein-protein interfaces with very favorable computed binding energies but which do not appear to be formed in experiments, suggesting that there may be important physical chemistry missing in the energy calculations. A total of 28 research groups took up the challenge of determining what is missing: we provided structures of 87 designed complexes and 120 naturally occurring complexes and asked participants to identify energetic contributions and/or structural features that distinguish between the two sets. The community found that electrostatics and solvation terms partially distinguish the designs from the natural complexes, largely due to the nonpolar character of the designed interactions. Beyond this polarity difference, the community found that the designed binding surfaces were, on average, structurally less embedded in the designed monomers, suggesting that backbone conformational rigidity at the designed surface is important for realization of the designed function. These results can be used to improve computational design strategies, but there is still much to be learned; for example, one designed complex, which does form in experiments, was classified by all metrics as a nonbinder.

H-bond network optimization in protein-protein complexes: are all-atom force field scores enough?

Structural prediction of protein-protein complexes given the structures of the two interacting compounds in their unbound state is a key problem in biophysics. In addition to the problem of sampling of near-native orientations, one of the modeling main difficulties is to discriminate true from false positives. Here, we present a hierarchical protocol for docking refinement able to discriminate near native poses from a group of docking candidates. The main idea is to combine an efficient sampling of the full system hydrogen bond network and side chains, together with an all-atom force field and a surface generalized born implicit solvent. We tested our method on a set of twenty two complexes containing a near-native solution within the top 100 docking poses, obtaining a near native solution as the top pose in 70% of the cases. We show that all atom force fields optimized H-bond networks do improve significantly state of the art scoring functions.

Prediction of protein-binding areas by small-world residue networks and application to docking

BACKGROUND

Protein-protein interactions are involved in most cellular processes, and their detailed physico-chemical and structural characterization is needed in order to understand their function at the molecular level. In-silico docking tools can complement experimental techniques, providing three-dimensional structural models of such interactions at atomic resolution. In several recent studies, protein structures have been modeled as networks (or graphs), where the nodes represent residues and the connecting edges their interactions. From such networks, it is possible to calculate different topology-based values for each of the nodes, and to identify protein regions with high centrality scores, which are known to positively correlate with key functional residues, hot spots, and protein-protein interfaces.

RESULTS

Here we show that this correlation can be efficiently used for the scoring of rigid-body docking poses. When integrated into the pyDock energy-based docking method, the new combined scoring function significantly improved the results of the individual components as shown on a standard docking benchmark. This improvement was particularly remarkable for specific protein complexes, depending on the shape, size, type, or flexibility of the proteins involved.

CONCLUSIONS

The network-based representation of protein structures can be used to identify protein-protein binding regions and to efficiently score docking poses, complementing energy-based approaches.