Evolutionary conservation of RecA genes in relation to protein structure and function

Intracellular Ebola virus nucleocapsid assembly revealed by in situ cryo-electron tomography

Filoviruses, including the Ebola and Marburg viruses, cause hemorrhagic fevers with up to 90% lethality. The viral nucleocapsid is assembled by polymerization of the nucleoprotein (NP) along the viral genome, together with the viral proteins VP24 and VP35. We employed cryo-electron tomography of cells transfected with viral proteins and infected with model Ebola virus to illuminate assembly intermediates, as well as a 9 Å map of the complete intracellular assembly. This structure reveals a previously unresolved third and outer layer of NP complexed with VP35. The intrinsically disordered region, together with the C-terminal domain of this outer layer of NP, provides the constant width between intracellular nucleocapsid bundles and likely functions as a flexible tether to the viral matrix protein in the virion. A comparison of intracellular nucleocapsids with prior in-virion nucleocapsid structures reveals that the nucleocapsid further condenses vertically in the virion. The interfaces responsible for nucleocapsid assembly are highly conserved and offer targets for broadly effective antivirals.

Conditional protein splicing of the Mycobacterium tuberculosis RecA intein in its native host

The recA gene, encoding Recombinase A (RecA) is one of three Mycobacterium tuberculosis (Mtb) genes encoding an in-frame intervening protein sequence (intein) that must splice out of precursor host protein to produce functional protein. Ongoing debate about whether inteins function solely as selfish genetic elements or benefit their host cells requires understanding of interplay between inteins and their hosts. We measured environmental effects on native RecA intein splicing within Mtb using a combination of western blots and promoter reporter assays. RecA splicing was stimulated in bacteria exposed to DNA damaging agents or by treatment with copper in hypoxic, but not normoxic, conditions. Spliced RecA was processed by the Mtb proteasome, while free intein was degraded efficiently by other unknown mechanisms. Unspliced precursor protein was not observed within Mtb despite its accumulation during ectopic expression of Mtb recA within E. coli. Surprisingly, Mtb produced free N-extein in some conditions, and ectopic expression of Mtb N-extein activated LexA in E. coli. These results demonstrate that the bacterial environment greatly impacts RecA splicing in Mtb, underscoring the importance of studying intein splicing in native host environments and raising the exciting possibility of intein splicing as a novel regulatory mechanism in Mtb.

A comparative genomic and phenotypic study of Vibrio cholerae model strains using hybrid sequencing

Next-generation sequencing methods have become essential for studying bacterial biology and pathogenesis, often depending on high-quality, closed genomes. In this study, we utilized a hybrid sequencing approach to assemble the genome of C6706, a widely used Vibrio cholerae model strain. We present a manually curated annotation of the genome, enhancing user accessibility by linking each coding sequence to its counterpart in N16961, the first sequenced V. cholerae isolate and a commonly used reference genome. Comparative genomic analysis between V. cholerae C6706 and N16961 uncovered multiple genetic differences in genes associated with key biological functions. To determine whether these genetic variations result in phenotypic differences, we compared several phenotypes relevant to V. cholerae pathogenicity like genetic stability, acid sensitivity, biofilm formation and motility. Notably, V. cholerae N16961 exhibited greater motility and reduced biofilm formation compared to V. cholerae C6706. These phenotypic differences appear to be mediated by variations in quorum sensing and cyclic di-GMP signalling pathways between the strains. This study provides valuable insights into the regulation of biofilm formation and motility in V. cholerae.

Conditional protein splicing of the Mycobacterium tuberculosis RecA intein in its native host

The recA gene, encoding Recombinase A (RecA) is one of three Mycobacterium tuberculosis (Mtb) genes encoding an in-frame intervening protein sequence (intein) that must splice out of precursor host protein to produce functional protein. Ongoing debate about whether inteins function solely as selfish genetic elements or benefit their host cells requires understanding of interplay between inteins and their hosts. We measured environmental effects on native RecA intein splicing within Mtb using a combination of western blots and promoter reporter assays. RecA splicing was stimulated in bacteria exposed to DNA damaging agents or by treatment with copper in hypoxic, but not normoxic, conditions. Spliced RecA was processed by the Mtb proteasome, while free intein was degraded efficiently by other unknown mechanisms. Unspliced precursor protein was not observed within Mtb despite its accumulation during ectopic expression of Mtb recA within E. coli. Surprisingly, Mtb produced free N-extein in some conditions, and ectopic expression of Mtb N-extein activated LexA in E. coli. These results demonstrate that the bacterial environment greatly impacts RecA splicing in Mtb, underscoring the importance of studying intein splicing in native host environments and raising the exciting possibility of intein splicing as a novel regulatory mechanism in Mtb.

Homology recognition without double-stranded DNA-strand separation in D-loop formation by RecA

RecA protein and RecA/Rad51 orthologues are required for homologous recombination and DNA repair in all living creatures. RecA/Rad51 catalyzes formation of the D-loop, an obligatory recombination intermediate, through an ATP-dependent reaction consisting of two phases: homology recognition between double-stranded (ds)DNA and single-stranded (ss)DNA to form a hybrid-duplex core of 6-8 base pairs and subsequent hybrid-duplex/D-loop processing. How dsDNA recognizes homologous ssDNA is controversial. The aromatic residue at the tip of the β-hairpin loop (L2) was shown to stabilize dsDNA-strand separation. We tested a model in which dsDNA strands were separated by the aromatic residue before homology recognition and found that the aromatic residue was not essential to homology recognition, but was required for D-loop processing. Contrary to the model, we found that the double helix was not unwound even a single turn during search for sequence homology, but rather was unwound only after the homologous sequence was recognized. These results suggest that dsDNA recognizes its homologous ssDNA before strand separation. The search for homologous sequence with homologous ssDNA without dsDNA-strand separation does not generate stress within the dsDNA; this would be an advantage for dsDNA to express homology-dependent functions in vivo and also in vitro.

Evolutionary conservation of sequence motifs at sites of protein modification

Gene duplications are common in biology and are likely to be an important source of functional diversification and specialization. The yeast Saccharomyces cerevisiae underwent a whole-genome duplication event early in evolution, and a substantial number of duplicated genes have been retained. We identified more than 3500 instances where only one of two paralogous proteins undergoes posttranslational modification despite having retained the same amino acid residue in both. We also developed a web-based search algorithm (CoSMoS.c.) that scores conservation of amino acid sequences based on 1011 wild and domesticated yeast isolates and used it to compare differentially modified pairs of paralogous proteins. We found that the most common modifications-phosphorylation, ubiquitylation, and acylation but not N-glycosylation-occur in regions of high sequence conservation. Such conservation is evident even for ubiquitylation and succinylation, where there is no established 'consensus site' for modification. Differences in phosphorylation were not associated with predicted secondary structure or solvent accessibility but did mirror known differences in kinase-substrate interactions. Thus, differences in posttranslational modification likely result from differences in adjoining amino acids and their interactions with modifying enzymes. By integrating data from large-scale proteomics and genomics analysis, in a system with such substantial genetic diversity, we obtained a more comprehensive understanding of the functional basis for genetic redundancies that have persisted for 100 million years.

Identification of amino acid domains of Borrelia burgdorferi P66 that are surface exposed and important for localization, oligomerization, and porin function of the protein

P66, a bifunctional integral outer membrane protein, is necessary for Borrelia burgdorferi to establish initial infection and to disseminate in mice. The integrin binding function of P66 facilitates extravasation and dissemination, but the role of its porin function during murine infection has not been investigated. A limitation to studying P66 porin function during mammalian infection has been the lack of structural information for P66. In this study, we experimentally characterized specific domains of P66 with regard to structure and function. First, we aligned the amino acid sequences of P66 from Lyme disease-causing Borrelia and relapsing fever-causing Borrelia to identify conserved and unique domains between these disease-causing clades. Then, we examined whether specific domains of P66 are exposed on the surface of the bacteria by introducing c-Myc epitope tags into each domain of interest. The c-Myc epitope tag inserted C-terminally to E33 (highly conserved domain), to T187 (integrin binding region domain and a non-conserved domain), and to E334 (non-conserved domain) were all detected on the surface of Borrelia burgdorferi. The c-Myc epitope tag inserted C-terminally to E33 and D303 in conserved domains disrupted P66 oligomerization and porin function. In a murine model of infection, the E33 and D303 mutants exhibited decreased infectivity and dissemination. Taken together, these results suggest the importance of these conserved domains, and potentially P66 porin function, in vivo.

A Genome-Phenome Association study in native microbiomes identifies a mechanism for cytosine modification in DNA and RNA

Shotgun metagenomic sequencing is a powerful approach to study microbiomes in an unbiased manner and of increasing relevance for identifying novel enzymatic functions. However, the potential of metagenomics to relate from microbiome composition to function has thus far been underutilized. Here, we introduce the Metagenomics Genome-Phenome Association (MetaGPA) study framework, which allows linking genetic information in metagenomes with a dedicated functional phenotype. We applied MetaGPA to identify enzymes associated with cytosine modifications in environmental samples. From the 2365 genes that met our significance criteria, we confirm known pathways for cytosine modifications and proposed novel cytosine-modifying mechanisms. Specifically, we characterized and identified a novel nucleic acid-modifying enzyme, 5-hydroxymethylcytosine carbamoyltransferase, that catalyzes the formation of a previously unknown cytosine modification, 5-carbamoyloxymethylcytosine, in DNA and RNA. Our work introduces MetaGPA as a novel and versatile tool for advancing functional metagenomics.

Many industrial processes, such as starch processing and oil refinement, use chemicals that cause harm to the environment. These can often be switched to more sustainable biological processes that are powered by proteins called enzymes.

Enzymes are micro-factories that speed up biochemical reactions in most living things. Communities of microorganisms (also known as microbiomes) are an amazing but often untapped resource for discovering enzymes that can be harnessed for industrial purposes. To gain a better picture of the microbes present within a population, researchers often extract and sequence the genetic material of all microorganisms in an environmental sample, also known as the metagenome. While current methods for analyzing the metagenome are good at identifying new species, they often provide limited information about the microorganism’s functional role within the community. This makes it difficult to find new enzymes that may be useful for industry.

Here, Yang, Lin et al. have developed a new technique called Metagenomics Genome-Phenome Association, or MetaGPA for short. The method works in a similar way to genome-wide association studies (GWAS) which are used to identify genes involved in human disease. However, instead of disease associated genes in humans, MetaGPA finds microbial genes that are associated with a biological process useful for biotechnology.

Like GWAS, the new approach created by Yang, Lin et al. compares two groups: the first contains microorganisms that carry out a specific process, and the second contains all organisms in the microbiome. The metagenome of each group is extracted and a computational pipeline is then applied to identify genes, including those coding for enzymes, that are found more often in the group performing the desired task.

To test the technique, Yang, Lin et al. used MetGPA to find new enzymes involved in DNA modification. Microbiome samples were collected from coastal water and sewage, and the computational pipeline was applied to discover genes that are associated with this process. Further analysis revealed that one of the identified genes codes for an enzyme that introduces a previously unknown change to DNA.

MetaGPA could be applied to other processes and microbiomes, and, if successful, may help researchers to identify more diverse enzymes than is currently available. This could scale up the discovery of new enzymes that can be used to power industrial reactions.

Biochemical characterization of Recombinase A from Wolbachia endosymbiont of filarial nematode Brugia malayi (wBmRecA)

Lymphatic filariasis is a debilitating disease that affects over 890 million people in 49 countries. A lack of vaccines, non-availability of adulticidal drugs, the threat of emerging drug resistance against available chemotherapeutics and an incomplete understanding of the immunobiology of the disease have sustained the problem. Characterization of Wolbachia proteins, the bacterial endosymbiont which helps in the growth and development of filarial worms, regulates fecundity in female worms and mediates immunopathogenesis of Lymphatic Filariasis, is an important approach to gain insights into the immunopathogenesis of the disease. In this study, we carried out extensive biochemical characterization of Recombinase A from Wolbachia of the filarial nematode Brugia malayi (wBmRecA) using an Electrophoretic Mobility Shift Assay, an ATP binding and hydrolysis assay, DNA strand exchange reactions, DAPI displacement assay and confocal microscopy, and evaluated anti-filarial activity of RecA inhibitors. Confocal studies showed that wBmRecA was expressed and localised within B. malayi microfilariae (Mf) and uteri and lateral chord of adult females. Recombinant wBmRecA was biochemically active and showed intrinsic binding capacity towards both single-stranded DNA and double-stranded DNA that were enhanced by ATP, suggesting ATP-induced cooperativity. wBmRecA promoted ATP hydrolysis and DNA strand exchange reactions in a concentration-dependent manner, and its binding to DNA was sensitive to temperature, pH and salt concentration. Importantly, the anti-parasitic drug Suramin, and Phthalocyanine tetrasulfonate (PcTs)-based inhibitors Fe-PcTs and 3,4-Cu-PcTs, inhibited wBmRecA activity and affected the motility and viability of Mf. The addition of Doxycycline further enhanced microfilaricidal activity of wBmRecA, suggesting potential synergism. Taken together, the omnipresence of wBmRecA in B. malayi life stages and the potent microfilaricidal activity of RecA inhibitors suggest an important role of wBmRecA in filarial pathogenesis.

Targeting evolution to inhibit antibiotic resistance

Drug-resistant bacterial infections have led to a global health crisis. Although much effort is placed on the development of new antibiotics or variants that are less subject to existing resistance mechanisms, history shows that this strategy by itself is unlikely to solve the problem of drug resistance. Here, we discuss inhibiting evolution as a strategy that, in combination with antibiotics, may resolve the problem. Although mutagenesis is the main driver of drug resistance development, attacking the drivers of genetic diversification in pathogens has not been well explored. Bacteria possess active mechanisms that increase the rate of mutagenesis, especially at times of stress, such as during replication within eukaryotic host cells, or exposure to antibiotics. We highlight how the existence of these promutagenic proteins (evolvability factors) presents an opportunity that can be capitalized upon for the effective inhibition of drug resistance development. To help move this idea from concept to execution, we first describe a set of criteria that an 'optimal' evolvability factor would likely have to meet to be a viable therapeutic target. We then discuss the intricacies of some of the known mutagenic mechanisms and evaluate their potential as drug targets to inhibit evolution. In principle, and as suggested by recent studies, we argue that the inhibition of these and other evolvability factors should reduce resistance development. Finally, we discuss the challenges of transitioning anti-evolution drugs from the laboratory to the clinic.

Cancer-associated mutations in the ribosomal protein L5 gene dysregulate the HDM2/p53-mediated ribosome biogenesis checkpoint

Perturbations in ribosome biogenesis have been associated with cancer. Such aberrations activate p53 through the RPL5/RPL11/5S rRNA complex-mediated inhibition of HDM2. Studies using animal models have suggested that this signaling pathway might constitute an important anticancer barrier. To gain a deeper insight into this issue in humans, here we analyze somatic mutations in RPL5 and RPL11 coding regions, reported in The Cancer Genome Atlas and International Cancer Genome Consortium databases. Using a combined computational and statistical approach, complemented by a range of biochemical and functional analyses in human cancer cell models, we demonstrate the existence of several mechanisms by which RPL5 mutations may impair wild-type p53 upregulation and ribosome biogenesis. Unexpectedly, the same approach provides only modest evidence for a similar role of RPL11, suggesting that RPL5 represents a preferred target during human tumorigenesis in cancers with wild-type p53. Furthermore, we find that several functional cancer-associated RPL5 somatic mutations occur as rare germline variants in general population. Our results shed light on the so-far enigmatic role of cancer-associated mutations in genes encoding ribosomal proteins, with implications for our understanding of the tumor suppressive role of the RPL5/RPL11/5S rRNA complex in human malignancies.

A multilocus sequence typing scheme for Mycobacterium abscessus complex (MAB-multilocus sequence typing) using whole-genome sequencing data

Background

Mycobacterium abscessus is a rapid growing nontuberculous mycobacteria (NTM) and a clinically significant pathogen capable of causing variable infections in humans that are difficult to treat and may require months of therapy/surgical interventions. Like other NTMs, M. abscessus can be associated with outbreaks leading to complex investigations and treatment of affected cases. Typing schemes for bacterial pathogens provide numerous applications; including identifying chain of transmission and tracking genomic evolution, are lacking or limited for many NTMs including M. abscessus.

Methods

We extended the existing scheme from PubMLST using whole-genome data for M. abscessus by extracting data for 15 genetic regions within the M. abscessus genome. A total of 168 whole genomes and 11 gene sequences were used to build this scheme (MAB-multilocus sequence typing [MLST]).

Results

All seven genes from the PubMLST scheme, namely argH, cya, gnd, murC, pta, purH, and rpoB, were expanded by 10, 14, 13, 10, 13, 10, and 9 alleles, respectively. Another eight novel genes were added including hsp 65, erm(41), arr, rrs, rrl, gyrA, gyrB, and recA with 16, 16, 25, 7, 32, 35, 29, and 15 alleles, respectively, with 85 unique sequence types identified among all isolates.

Conclusion

MAB-MLST can provide identification of M. abscessus complex to the subspecies level based on three genes and can provide antimicrobial resistance susceptibility prediction based on results from seven genes. MAB-MLST generated a total of 85 STs, resulting in subtyping of 90 additional isolates that could not be genotyped using PubMLST and yielding results comparable to whole-genome sequencing (WGS). Implementation of a Galaxy-based data analysis tool, MAB-MLST, that simplifies the WGS data and yet maintains a high discriminatory index that can aid in deciphering an outbreak has vast applicability for routine diagnostics.

In vitro role of Rad54 in Rad51-ssDNA filament-dependent homology search and synaptic complexes formation

Homologous recombination (HR) uses a homologous template to accurately repair DNA double-strand breaks and stalled replication forks to maintain genome stability. During homology search, Rad51 nucleoprotein filaments probe and interact with dsDNA, forming the synaptic complex that is stabilized on a homologous sequence. Strand intertwining leads to the formation of a displacement-loop (D-loop). In yeast, Rad54 is essential for HR in vivo and required for D-loop formation in vitro, but its exact role remains to be fully elucidated. Using electron microscopy to visualize the DNA-protein complexes, here we find that Rad54 is crucial for Rad51-mediated synaptic complex formation and homology search. The Rad54−K341R ATPase-deficient mutant protein promotes formation of synaptic complexes but not D-loops and leads to the accumulation of stable heterologous associations, suggesting that the Rad54 ATPase is involved in preventing non-productive intermediates. We propose that Rad51/Rad54 form a functional unit operating in homology search, synaptic complex and D-loop formation.

Homologous recombination uses a template to accurately repair DNA double-strand breaks and stalled replication forks to maintain genome stability. Here authors use electron microscopy to investigate the role of Rad54 in homology search and synaptic complex formation.

Human NK cell receptor KIR2DS4 detects a conserved bacterial epitope presented by HLA-C

Natural killer (NK) cells have an important role in immune defense against viruses and cancer. Activation of human NK cell cytotoxicity toward infected or tumor cells is regulated by killer cell immunoglobulin-like receptors (KIRs) that bind to human leukocyte antigen class I (HLA-I). Combinations of KIR with HLA-I are genetically associated with susceptibility to disease. KIR2DS4, an activating member of the KIR family with poorly defined ligands, is a receptor of unknown function. Here, we show that KIR2DS4 has a strong preference for rare peptides carrying a Trp at position 8 (p8) of 9-mer peptides bound to HLA-C*05:01. The complex of a peptide bound to HLA-C*05:01 with a Trp at p8 was sufficient for activation of primary KIR2DS4+ NK cells, independent of activation by other receptors and of prior NK cell licensing. HLA-C*05:01+ cells that expressed the peptide epitope triggered KIR2DS4+ NK cell degranulation. We show an inverse correlation of the worldwide allele frequency of functional KIR2DS4 with that of HLA-C*05:01, indicative of functional interaction and balancing selection. We found a highly conserved peptide sequence motif for HLA-C*05:01-restricted activation of human KIR2DS4+ NK cells in bacterial recombinase A (RecA). KIR2DS4+ NK cells were stimulated by RecA epitopes from multiple human pathogens, including Helicobacter, Chlamydia, Brucella, and Campylobacter. We predict that over 1,000 bacterial species could activate NK cells through KIR2DS4, and propose that human NK cells also contribute to immune defense against bacteria through recognition of a conserved RecA epitope presented by HLA-C*05:01.

Ureaplasma diversum protein interaction networks: evidence of horizontal gene transfer and evolution of reduced genomes among Mollicutes

Ureaplasma diversum is a member of the Mollicutes class responsible for urogenital tract infection in cattle and small ruminants. Studies indicate that the process of horizontal gene transfer, the exchange of genetic material among different species, has a crucial role in mollicute evolution, affecting the group's characteristic genomic reduction process and simplification of metabolic pathways. Using bioinformatics tools and the STRING database of known and predicted protein interactions, we constructed the protein-protein interaction network of U. diversum and compared it with the networks of other members of the Mollicutes class. We also investigated horizontal gene transfer events in subnetworks of interest involved in purine and pyrimidine metabolism and urease function, chosen because of their intrinsic importance for host colonization and virulence. We identified horizontal gene transfer events among Mollicutes and from Ureaplasma to Staphylococcus aureus and Corynebacterium, bacterial groups that colonize the urogenital niche. The overall tendency of genome reduction and simplification in the Mollicutes is echoed in their protein interaction networks, which tend to be more generalized and less selective. Our data suggest that the process was permitted (or enabled) by an increase in host dependence and the available gene repertoire in the urogenital tract shared via horizontal gene transfer.

The Use of Biosensors to Explore the Potential of Probiotic Strains to Reduce the SOS Response and Mutagenesis in Bacteria

A model system based on the Escherichia coli MG1655 (pRecA-lux) Lux-biosensor was used to evaluate the ability of the fermentates of eight probiotic strains to reduce the SOS response stimulated by ciprofloxacin in bacteria and mutagenesis mediated by it. Preliminary attempts to estimate the chemical nature of active components of the fermentates were conducted.

Uncovering missing pieces: duplication and deletion history of arrestins in deuterostomes

BACKGROUND

The cytosolic arrestin proteins mediate desensitization of activated G protein-coupled receptors (GPCRs) via competition with G proteins for the active phosphorylated receptors. Arrestins in active, including receptor-bound, conformation are also transducers of signaling. Therefore, this protein family is an attractive therapeutic target. The signaling outcome is believed to be a result of structural and sequence-dependent interactions of arrestins with GPCRs and other protein partners. Here we elucidated the detailed evolution of arrestins in deuterostomes.

RESULTS

Identity and number of arrestin paralogs were determined searching deuterostome genomes and gene expression data. In contrast to standard gene prediction methods, our strategy first detects exons situated on different scaffolds and then solves the problem of assigning them to the correct gene. This increases both the completeness and the accuracy of the annotation in comparison to conventional database search strategies applied by the community. The employed strategy enabled us to map in detail the duplication- and deletion history of arrestin paralogs including tandem duplications, pseudogenizations and the formation of retrogenes. The two rounds of whole genome duplications in the vertebrate stem lineage gave rise to four arrestin paralogs. Surprisingly, visual arrestin ARR3 was lost in the mammalian clades Afrotheria and Xenarthra. Duplications in specific clades, on the other hand, must have given rise to new paralogs that show signatures of diversification in functional elements important for receptor binding and phosphate sensing.

CONCLUSION

The current study traces the functional evolution of deuterostome arrestins in unprecedented detail. Based on a precise re-annotation of the exon-intron structure at nucleotide resolution, we infer the gain and loss of paralogs and patterns of conservation, co-variation and selection.

RecA Inhibitors Potentiate Antibiotic Activity and Block Evolution of Antibiotic Resistance

Antibiotic resistance arises from the maintenance of resistance mutations or genes acquired from the acquisition of adaptive de novo mutations or the transfer of resistance genes. Antibiotic resistance is acquired in response to antibiotic therapy by activating SOS-mediated DNA repair and mutagenesis and horizontal gene transfer pathways. Initiation of the SOS pathway promotes activation of RecA, inactivation of LexA repressor, and induction of SOS genes. Here, we have identified and characterized phthalocyanine tetrasulfonic acid RecA inhibitors that block antibiotic-induced activation of the SOS response. These inhibitors potentiate the activity of bactericidal antibiotics, including members of the quinolone, β-lactam, and aminoglycoside families in both Gram-negative and Gram-positive bacteria. They reduce the ability of bacteria to acquire antibiotic resistance mutations and to transfer mobile genetic elements conferring resistance. This study highlights the advantage of including RecA inhibitors in bactericidal antibiotic therapies and provides a new strategy for prolonging antibiotic shelf life.

Predicting Protein-Protein Interaction Sites Using Sequence Descriptors and Site Propensity of Neighboring Amino Acids

Information about the interface sites of Protein–Protein Interactions (PPIs) is useful for many biological research works. However, despite the advancement of experimental techniques, the identification of PPI sites still remains as a challenging task. Using a statistical learning technique, we proposed a computational tool for predicting PPI interaction sites. As an alternative to similar approaches requiring structural information, the proposed method takes all of the input from protein sequences. In addition to typical sequence features, our method takes into consideration that interaction sites are not randomly distributed over the protein sequence. We characterized this positional preference using protein complexes with known structures, proposed a numerical index to estimate the propensity and then incorporated the index into a learning system. The resulting predictor, without using structural information, yields an area under the ROC curve (AUC) of 0.675, recall of 0.597, precision of 0.311 and accuracy of 0.583 on a ten-fold cross-validation experiment. This performance is comparable to the previous approach in which structural information was used. Upon introducing the B-factor data to our predictor, we demonstrated that the AUC can be further improved to 0.750. The tool is accessible at http://bsaltools.ym.edu.tw/predppis.

Investigating direct interaction between Escherichia coli topoisomerase I and RecA

Protein-protein interactions are of special importance in cellular processes, including replication, transcription, recombination, and repair. Escherichia coli topoisomerase I (EcTOP1) is primarily involved in the relaxation of negative DNA supercoiling. E. coli RecA, the key protein for homologous recombination and SOS DNA-damage response, has been shown to stimulate the relaxation activity of EcTOP1. The evidence for their direct protein-protein interaction has not been previously established. We report here the direct physical interaction between E. coli RecA and topoisomerase I. We demonstrated the RecA-topoisomerase I interaction via pull-down assays, and surface plasmon resonance measurements. Molecular docking supports the observation that the interaction involves the topoisomerase I N-terminal domains that form the active site. Our results from pull-down assays showed that ATP, although not required, enhances the RecA-EcTOP1 interaction. We propose that E. coli RecA physically interacts with topoisomerase I to modulate the chromosomal DNA supercoiling.

Anionic Phospholipids Stabilize RecA Filament Bundles in Escherichia coli

We characterize the interaction of RecA with membranes in vivo and in vitro and demonstrate that RecA binds tightly to the anionic phospholipids cardiolipin (CL) and phosphatidylglycerol (PG). Using computational models, we identify two regions of RecA that interact with PG and CL: (1) the N-terminal helix and (2) loop L2. Mutating these regions decreased the affinity of RecA to PG and CL in vitro. Using 3D super-resolution microscopy, we demonstrate that depleting Escherichia coli PG and CL altered the localization of RecA foci and hindered the formation of RecA filament bundles. Consequently, E. coli cells lacking aPLs fail to initiate a robust SOS response after DNA damage, indicating that the membrane acts as a scaffold for nucleating the formation of RecA filament bundles and plays an important role in the SOS response.

Highly thermostable RadA protein from the archaeon Pyrococcus woesei enhances specificity of simplex and multiplex PCR assays

The radA gene of the hyperthermophilic archaeon Pyrococcus woesei (Thermococcales) was cloned and overexpressed in Escherichia coli. The 1050-bp gene codes for a 349-amino-acid polypeptide with an M r of 38,397 which shows 100 % positional amino acid identity to Pyrococcus furiosus RadA and 27.1 % to the E. coli RecA protein. Recombinant RadA was overproduced in Escherichia coli as a His-tagged fusion protein and purified to electrophoretic homogeneity using a simple procedure consisting of ammonium sulfate precipitation and metal-affinity chromatography. In solution RadA exists as an undecamer (11-mer). The protein binds both to ssDNA and dsDNA. RadA has been found to be highly thermostable, it remains almost unaffected by a 4-h incubation at 94 °C. The addition of the RadA protein to either simplex or multiplex PCR assays, significantly improves the specificity of DNA amplification by eliminating non-specific products. Among applications tested the RadA protein proved to be useful in allelic discrimination assay of HADHA gene associated with long-chain 3-hydroxylacyl-CoA dehydrogenase deficiency that in infancy may lead to hypotonia, serious heart and liver problems and even sudden death.

An integrative approach to the study of filamentous oligomeric assemblies, with application to RecA

Oligomeric macromolecules in the cell self-organize into a wide variety of geometrical motifs such as helices, rings or linear filaments. The recombinase proteins involved in homologous recombination present many such assembly motifs. Here, we examine in particular the polymorphic characteristics of RecA, the most studied member of the recombinase family, using an integrative approach that relates local modes of monomer/monomer association to the global architecture of their screw-type organization. In our approach, local modes of association are sampled via docking or Monte Carlo simulations. This enables shedding new light on fiber morphologies that may be adopted by the RecA protein. Two distinct RecA helical morphologies, the so-called "extended" and "compressed" forms, are known to play a role in homologous recombination. We investigate the variability within each form in terms of helical parameters and steric accessibility. We also address possible helical discontinuities in RecA filaments due to multiple monomer-monomer association modes. By relating local interface organization to global filament morphology, the strategies developed here to study RecA self-assembly are particularly well suited to other DNA-binding proteins and to filamentous protein assemblies in general.

Loop L1 governs the DNA-binding specificity and order for RecA-catalyzed reactions in homologous recombination and DNA repair

In all organisms, RecA-family recombinases catalyze homologous joint formation in homologous genetic recombination, which is essential for genome stability and diversification. In homologous joint formation, ATP-bound RecA/Rad51-recombinases first bind single-stranded DNA at its primary site and then interact with double-stranded DNA at another site. The underlying reason and the regulatory mechanism for this conserved binding order remain unknown. A comparison of the loop L1 structures in a DNA-free RecA crystal that we originally determined and in the reported DNA-bound active RecA crystals suggested that the aspartate at position 161 in loop L1 in DNA-free RecA prevented double-stranded, but not single-stranded, DNA-binding to the primary site. This was confirmed by the effects of the Ala-replacement of Asp-161 (D161A), analyzed directly by gel-mobility shift assays and indirectly by DNA-dependent ATPase activity and SOS repressor cleavage. When RecA/Rad51-recombinases interact with double-stranded DNA before single-stranded DNA, homologous joint-formation is suppressed, likely by forming a dead-end product. We found that the D161A-replacement reduced this suppression, probably by allowing double-stranded DNA to bind preferentially and reversibly to the primary site. Thus, Asp-161 in the flexible loop L1 of wild-type RecA determines the preference for single-stranded DNA-binding to the primary site and regulates the DNA-binding order in RecA-catalyzed recombinase reactions.

Local packing density is the main structural determinant of the rate of protein sequence evolution at site level

Functional and biophysical constraints result in site-dependent patterns of protein sequence variability. It is commonly assumed that the key structural determinant of site-specific rates of evolution is the Relative Solvent Accessibility (RSA). However, a recent study found that amino acid substitution rates correlate better with two Local Packing Density (LPD) measures, the Weighted Contact Number (WCN) and the Contact Number (CN), than with RSA. This work aims at a more thorough assessment. To this end, in addition to substitution rates, we considered four other sequence variability scores, four measures of solvent accessibility (SA), and other CN measures. We compared all properties for each protein of a structurally and functionally diverse representative dataset of monomeric enzymes. We show that the best sequence variability measures take into account phylogenetic tree topology. More importantly, we show that both LPD measures (WCN and CN) correlate better than all of the SA measures, regardless of the sequence variability score used. Moreover, the independent contribution of the best LPD measure is approximately four times larger than that of the best SA measure. This study strongly supports the conclusion that a site's packing density rather than its solvent accessibility is the main structural determinant of its rate of evolution.

Discovery and characterization of RecA protein of thermophilic bacterium Thermus thermophilus MAT72 phage Tt72 that increases specificity of a PCR-based DNA amplification

The recA gene of newly discovered Thermus thermophilus MAT72 phage Tt72 (Myoviridae) was cloned and overexpressed in Escherichia coli. The 1020-bp gene codes for a 339-amino-acid polypeptide with an Mr of 38,155 which shows 38.7% positional identity to the E. coli RecA protein. When expressed in E. coli, the Tt72 recA gene did not confer the ability to complement the ultraviolet light (254nm) sensitivity of an E. coli recA mutant. Tt72 RecA protein has been purified with good yield to catalytic and electrophoretic homogeneity using a three-step chromatography procedure. Biochemical characterization indicated that the protein can pair and promote ATP-dependent strand exchange reaction resulting in formation of a heteroduplex DNA at 60°C under conditions otherwise optimal for E. coli RecA. When the Tt72 RecA protein was included in a standard PCR-based DNA amplification reaction, the specificity of the PCR assays was significantly improved by eliminating non-specific products.

Functional diversification of ROK-family transcriptional regulators of sugar catabolism in the Thermotogae phylum

Large and functionally heterogeneous families of transcription factors have complex evolutionary histories. What shapes specificities toward effectors and DNA sites in paralogous regulators is a fundamental question in biology. Bacteria from the deep-branching lineage Thermotogae possess multiple paralogs of the repressor, open reading frame, kinase (ROK) family regulators that are characterized by carbohydrate-sensing domains shared with sugar kinases. We applied an integrated genomic approach to study functions and specificities of regulators from this family. A comparative analysis of 11 Thermotogae genomes revealed novel mechanisms of transcriptional regulation of the sugar utilization networks, DNA-binding motifs and specific functions. Reconstructed regulons for seven groups of ROK regulators were validated by DNA-binding assays using purified recombinant proteins from the model bacterium Thermotoga maritima. All tested regulators demonstrated specific binding to their predicted cognate DNA sites, and this binding was inhibited by specific effectors, mono- or disaccharides from their respective sugar catabolic pathways. By comparing ligand-binding domains of regulators with structurally characterized kinases from the ROK family, we elucidated signature amino acid residues determining sugar-ligand regulator specificity. Observed correlations between signature residues and the sugar-ligand specificities provide the framework for structure functional classification of the entire ROK family.

Genome wide cloning of maize meiotic recombinase Dmc1 and its functional structure through molecular phylogeny

The development of meiotic division and associated genetic recombination paved the way for evolutionary changes. However, the secondary and tertiary structure and functional domains of many of the proteins involved in genetic recombination have not been studied in detail. We used the human Dmc1 gene product along with secondary and tertiary domain structures of Escherichia coli RecA protein to help determine the molecular structure and function of maize Dmc1, which is required for synaptonemal complex formation and cell cycle progression. The maize recombinase Dmc1 gene was cloned and characterized, using rice Dmc1 cDNA as an orthologue. The deduced amino acid sequence was used for elaborating its 3-D structure, and functional analysis was made with the CDD software, showing significant identity of the Dmc1 gene product in Zea mays with that of Homo sapiens. Based on these results, the domains and motives of WalkerA and WalkerB as ATP binding sites, a multimer site (BRC) interface, the putative ssDNA binding L1 and L2 loops, the putative dsDNA binding helix-hairpin-helix, a polymerization motif, the subunit rotation motif, and a small N-terminal domain were proposed for maize recombinase Dmc1.

Molecular evolution of the transmembrane domains of G protein-coupled receptors

G protein-coupled receptors (GPCRs) are a superfamily of integral membrane proteins vital for signaling and are important targets for pharmaceutical intervention in humans. Previously, we identified a group of ten amino acid positions (called key positions), within the seven transmembrane domain (7TM) interhelical region, which had high mutual information with each other and many other positions in the 7TM. Here, we estimated the evolutionary selection pressure at those key positions. We found that the key positions of receptors for small molecule natural ligands were under strong negative selection. Receptors naturally activated by lipids had weaker negative selection in general when compared to small molecule-activated receptors. Selection pressure varied widely in peptide-activated receptors. We used this observation to predict that a subgroup of orphan GPCRs not under strong selection may not possess a natural small-molecule ligand. In the subgroup of MRGX1-type GPCRs, we identified a key position, along with two non-key positions, under statistically significant positive selection.

Separation of recombination and SOS response in Escherichia coli RecA suggests LexA interaction sites

RecA plays a key role in homologous recombination, the induction of the DNA damage response through LexA cleavage and the activity of error-prone polymerase in Escherichia coli. RecA interacts with multiple partners to achieve this pleiotropic role, but the structural location and sequence determinants involved in these multiple interactions remain mostly unknown. Here, in a first application to prokaryotes, Evolutionary Trace (ET) analysis identifies clusters of evolutionarily important surface amino acids involved in RecA functions. Some of these clusters match the known ATP binding, DNA binding, and RecA-RecA homo-dimerization sites, but others are novel. Mutation analysis at these sites disrupted either recombination or LexA cleavage. This highlights distinct functional sites specific for recombination and DNA damage response induction. Finally, our analysis reveals a composite site for LexA binding and cleavage, which is formed only on the active RecA filament. These new sites can provide new drug targets to modulate one or more RecA functions, with the potential to address the problem of evolution of antibiotic resistance at its root.

In eubacteria, genome integrity is in large part orchestrated by RecA, which directly participates in recombination, induction of DNA damage response through LexA repressor cleavage and error-prone DNA synthesis. Yet, most of the interaction sites necessary for these vital processes are largely unknown. By comparing divergences among RecA sequences and computing putative functional regions, we discovered four functional sites of RecA. Targeted point-mutations were then tested for both recombination and DNA damage induction and reveal distinct RecA functions at each one of these sites. In particular, one new set of mutants is deficient in promoting LexA cleavage and yet maintains the ability to induce the DNA damage response. These results reveal specific amino acid determinants of the RecA–LexA interaction and suggest that LexA binds RecA_i and RecA_i+6 at a composite site on the RecA filament, which could explain the role of the active filament during LexA cleavage.

Relative von Neumann entropy for evaluating amino acid conservation

The Shannon entropy is a common way of measuring conservation of sites in multiple sequence alignments, and has also been extended with the relative Shannon entropy to account for background frequencies. The von Neumann entropy is another extension of the Shannon entropy, adapted from quantum mechanics in order to account for amino acid similarities. However, there is yet no relative von Neumann entropy defined for sequence analysis. We introduce a new definition of the von Neumann entropy for use in sequence analysis, which we found to perform better than the previous definition. We also introduce the relative von Neumann entropy and a way of parametrizing this in order to obtain the Shannon entropy, the relative Shannon entropy and the von Neumann entropy at special parameter values. We performed an exhaustive search of this parameter space and found better predictions of catalytic sites compared to any of the previously used entropies.

Modulating cellular recombination potential through alterations in RecA structure and regulation

The wild-type Escherichia coli RecA protein is a recombinase platform with unrealized recombination potential. We have explored the factors affecting recombination during conjugation with a quantitative assay. Regulatory proteins that affect RecA function have the capacity to increase or decrease recombination frequencies by factors up to sixfold. Autoinhibition by the RecA C-terminus can affect recombination frequency by factors up to fourfold. The greatest changes in recombination frequency measured here are brought about by point mutations in the recA gene. RecA variants can increase recombination frequencies by more than 50-fold. The RecA protein thus possesses an inherently broad functional range. The RecA protein of E. coli (EcRecA) is not optimized for recombination function. Instead, much of the recombination potential of EcRecA is structurally suppressed, probably reflecting cellular requirements. One point mutation in EcRecA with a particularly dramatic effect on recombination frequency, D112R, exhibits an enhanced capacity to load onto SSB-coated ssDNA, overcome the effects of regulatory proteins such as PsiB and RecX, and to pair homologous DNAs. Comparisons of key RecA protein mutants reveal two components to RecA recombination function - filament formation and the inherent DNA pairing activity of the formed filaments.

A comparative study of conservation and variation scores

BACKGROUND

Conservation and variation scores are used when evaluating sites in a multiple sequence alignment, in order to identify residues critical for structure or function. A variety of scores are available today but it is not clear how different scores relate to each other.

RESULTS

We applied 25 conservation and variation scores to alignments from the Catalytic Site Atlas (CSA). We calculated distances among scores based on correlation coefficients, and constructed a dendrogram of the scores by average linking cluster analysis. The cluster analysis showed that most scores fall into one of two groups--substitution matrix based group and frequency based group respectively. We also evaluated the scores' performance in predicting catalytic sites and found that frequency based scores generally perform best.

CONCLUSIONS

Conservation and variation scores can be classified into mainly two large groups. When using a score to predict catalytic sites, frequency based scores that also consider a background distribution are most successful.

UV hyper-resistance in Prochlorococcus MED4 results from a single base pair deletion just upstream of an operon encoding nudix hydrolase and photolyase

Exposure to solar radiation can cause mortality in natural communities of pico-phytoplankton, both at the surface and to a depth of at least 30 m. DNA damage is a significant cause of death, mainly due to cyclobutane pyrimidine dimer formation, which can be lethal if not repaired. While developing a UV mutagenesis protocol for the marine cyanobacterium Prochlorococcus, we isolated a UV-hyper-resistant variant of high light-adapted strain MED4. The hyper-resistant strain was constitutively upregulated for expression of the mutT-phrB operon, encoding nudix hydrolase and photolyase, both of which are involved in repair of DNA damage that can be caused by UV light. Photolyase (PhrB) breaks pyrimidine dimers typically caused by UV exposure, using energy from visible light in the process known as photoreactivation. Nudix hydrolase (MutT) hydrolyses 8-oxo-dGTP, an aberrant form of GTP that results from oxidizing conditions, including UV radiation, thus impeding mispairing and mutagenesis by preventing incorporation of the aberrant form into DNA. These processes are error-free, in contrast to error-prone SOS dark repair systems that are widespread in bacteria. The UV-hyper-resistant strain contained only a single mutation: a 1 bp deletion in the intergenic region directly upstream of the mutT-phrB operon. Two subsequent enrichments for MED4 UV-hyper-resistant strains from MED4 wild-type cultures gave rise to strains containing this same 1 bp deletion, affirming its connection to the hyper-resistant phenotype. These results have implications for Prochlorococcus DNA repair mechanisms, genome stability and possibly lysogeny.

Binding selectivity of RecA to a single stranded DNA, a computational approach

Homologous recombination (HR) is the major DNA double strand break repair pathway which maintains the genomic integrity. It is fundamental for the survivability and functionality of all organisms. One of the initial steps in HR is the formation of the nucleoprotein filament composed by a single stranded DNA chain surrounded by the recombinases protein. The filament orchestrates the search for an undamaged homologue, as a template for the repair process. Our theoretical study was aimed at elucidating the selectivity of the interaction between a monomer of the recombinases enzyme in the Escherichia coli, EcRecA, the bacterial homologue of human Rad51, with a series of oligonucleotides of nine bases length. The complex, equilibrated for 20 ns with Langevian dynamics, was inserted in a periodic box with a 8 Å buffer of water molecules explicitly described by the TIP3P model. The absolute binding free energies are calculated in an implicit solvent using the Poisson-Boltzmann (PB) and the generalized Born (GB) solvent accessible surface area, using the MM-PB(GB)SA model. The solute entropic contribution is also calculated by normal mode analysis. The results underline how a significant contribution of the binding free energy is due to the interaction with the Arg196, a critical amino acid for the activity of the enzyme. The study revealed how the binding affinity of EcRecA is significantly higher toward dT₉ rather than dA₉, as expected from the experimental results.

RecA-mediated cleavage of lambda cI repressor accepts repressor dimers: probable role of prolyl cis-trans isomerization and catalytic involvement of H163, K177, and K232 of RecA

The lambda cI repressor is found to be cleaved in the presence of activated RecA in its DNA-bound dimeric form at a rate similar to that in the absence of operator DNA in contrast to previous studies inferring repressor monomer as a preferred substrate. Though activated RecA does not possess any measurable isomerase activity against a standard peptide substrate, prolyl isomerase inhibitors cyclosporin A and rapamycin do inhibit RecA-mediated cleavage. Histidine and lysine to a smaller extent, are shown to cleave cI repressor in a non-enzymatic fashion whereas arginine and glutamate do not. When activated RecA filament is covalently modified by using an excess of diethyl pyrocarbonate or maleic anhydride, RecA-mediated cleavage of cI repressor is inhibited. Combining our chemical modification data with model building and earlier mutagenesis data, it is argued that H163, K177, and K232 in RecA are crucial residues involved in cI repressor cleavage by combining with the catalytic Ser149 and K192 in the repressor. It is suggested by model building that subunits n, n+4, and n+5 in the RecA filament contribute one loop each for holding the C-terminal domain of the repressor during cleavage within the RecA helical groove, explaining why its ADP-form is inactive and its ATP-form is active regarding repressor cleavage.

Prediction of catalytic residues using the variation of stereochemical properties

In this paper, we investigate a simple protein sequence conservation measure which takes amino acid similarity into account. Instead of grouping 20 amino acids into disjoint sets in previous methods, we consider ten overlapping classes. The method is based on the assumption that a column in a multiple sequence alignment is evolved from an identical column in the evolutionary history. Two ten-dimensional vectors are constructed for each position to denote frequencies of ten classes in a column and the corresponding hypothetical identical column. Then the cosine function of the angle between these two vectors is considered as a measure of divergence of stereochemical properties at this position. This divergence, combining with other conservation scores, is used as conservation measure of the column. Finally, we evaluate our methods by identifying catalytic sites, using rank analysis criterion and receiver operator characteristic analysis criterion.

The N-terminal domain of Escherichia coli RecA have multiple functions in promoting homologous recombination

Escherichia coli RecA mediates homologous recombination, a process essential to maintaining genome integrity. In the presence of ATP, RecA proteins bind a single-stranded DNA (ssDNA) to form a RecA-ssDNA presynaptic nucleoprotein filament that captures donor double-stranded DNA (dsDNA), searches for homology, and then catalyzes the strand exchange between ssDNA and dsDNA to produce a new heteroduplex DNA. Based upon a recently reported crystal structure of the RecA-ssDNA nucleoprotein filament, we carried out structural and functional studies of the N-terminal domain (NTD) of the RecA protein. The RecA NTD was thought to be required for monomer-monomer interaction. Here we report that it has two other distinct roles in promoting homologous recombination. It first facilitates the formation of a RecA-ssDNA presynaptic nucleoprotein filament by converting ATP to an ADP-Pi intermediate. Then, once the RecA-ssDNA presynaptic nucleoprotein filament is stably assembled in the presence of ATPγS, the NTD is required to capture donor dsDNA. Our results also suggest that the second function of NTD may be similar to that of Arg243 and Lys245, which were implicated earlier as binding sites of donor dsDNA. A two-step model is proposed to explain how a RecA-ssDNA presynaptic nucleoprotein filament interacts with donor dsDNA.

The RAD51 gene family in bread wheat is highly conserved across eukaryotes, with RAD51A upregulated during early meiosis

The RADiation sensitive protein 51 (RAD51) recombinase is a eukaryotic homologue of the bacterial Recombinase A (RecA). It is required for homologous recombination of DNA during meiosis where it plays a role in processes such as homology searching and strand invasion. RAD51 is well conserved in eukaryotes with as many as four paralogues identified in vertebrates and some higher plants. Here we report the isolation and preliminary characterisation of four RAD51 gene family members in hexaploid (bread) wheat (Triticum aestivum L.). RAD51A1, RAD51A2 and RAD51D were located on chromosome group 7, and RAD51C was on chromosome group 2. Q-PCR gene expression profiling revealed that RAD51A1 was upregulated during meiosis with lower expression levels seen in mitotic tissue, and bioinformatics analysis demonstrated the evolutionary linkages of this gene family to other eukaryotic RAD51 sequences. Western blot analysis of heterologously expressed RAD51 from bread wheat has shown that it is detectable using anti-human RAD51 antibodies and that molecular modelling of the same protein revealed structural conservation when compared with yeast, human, Arabidopsis and maize RAD51A orthologues. This report has widened the knowledge base of this important protein family in plants, and highlighted the high level of structural conservation among RAD51 proteins from various species.

Molecular Design and Functional Organization of the RecA Protein

The bacterial RecA protein participates in a remarkably diverse set of functions, all of which are involved in the maintenance of genomic integrity. RecA is a central component in both the catalysis of recombinational DNA repair and the regulation of the cellular SOS response. Despite the mechanistic differences of its functions, all require formation of an active RecA/ATP/DNA complex. RecA is a classic allosterically regulated enzyme, and ATP binding results in a dramatic increase in DNA binding affinity and a cooperative assembly of RecA subunits to form an ordered, helical nucleoprotein filament. The molecular events that underlie this ATP-induced structural transition are becoming increasingly clear. This review focuses on descriptions of our current understanding of the molecular design and allosteric regulation of RecA. We present a comprehensive list of all published recA mutants and use the results of various genetic and biochemical studies, together with available structural information, to develop ideas regarding the design of RecA functional domains and their catalytic organization.

The contrasting properties of conservation and correlated phylogeny in protein functional residue prediction

Background

Amino acids responsible for structure, core function or specificity may be inferred from multiple protein sequence alignments where a limited set of residue types are tolerated. The rise in available protein sequences continues to increase the power of techniques based on this principle.

Results

A new algorithm, SMERFS, for predicting protein functional sites from multiple sequences alignments was compared to 14 conservation measures and to the MINER algorithm. Validation was performed on an automatically generated dataset of 1457 families derived from the protein interactions database SNAPPI-DB, and a smaller manually curated set of 148 families. The best performing measure overall was Williamson property entropy, with ROC_0.1 scores of 0.0087 and 0.0114 for domain and small molecule contact prediction, respectively. The Lancet method performed worse than random on protein-protein interaction site prediction (ROC_0.1 score of 0.0008). The SMERFS algorithm gave similar accuracy to the phylogenetic tree-based MINER algorithm but was superior to Williamson in prediction of non-catalytic transient complex interfaces. SMERFS predicts sites that are significantly more solvent accessible compared to Williamson.

Conclusion

Williamson property entropy is the the best performing of 14 conservation measures examined. The difference in performance of SMERFS relative to Williamson in manually defined complexes was dependent on complex type. The best choice of analysis method is therefore dependent on the system of interest. Additional computation employed by Miner in calculation of phylogenetic trees did not produce improved results over SMERFS. SMERFS performance was improved by use of windows over alignment columns, illustrating the necessity of considering the local environment of positions when assessing their functional significance.