The Dundee Resource for Sequence Analysis and Structure Prediction

Classification of likely functional class for ligand binding sites identified from fragment screening

Fragment screening is used to identify binding sites and leads in drug discovery, but it is often unclear which binding sites are functionally important. Here, data from 37 experiments, and 1309 protein structures binding to 1601 ligands were analysed. A method to group ligands by binding sites is introduced and sites clustered according to profiles of relative solvent accessibility. This identified 293 unique ligand binding sites, grouped into four clusters (C1-4). C1 includes larger, buried, conserved, and population missense-depleted sites, enriched in known functional sites. C4 comprises smaller, accessible, divergent, missense-enriched sites, depleted in functional sites. A site in C1 is 28 times more likely to be functional than one in C4. Seventeen sites, which to the best of our knowledge are novel, in 13 proteins are identified as likely to be functionally important with examples from human tenascin and 5-aminolevulinate synthase highlighted. A multi-layer perceptron, and K-nearest neighbours model are presented to predict cluster labels for ligand binding sites with an accuracy of 96% and 100%, respectively, so allowing functional classification of sites for proteins not in this set. Our findings will be of interest to those studying protein-ligand interactions and developing new drugs or function modulators.

A genomic platform for surveillance and antigen discovery in Plasmodium spp. using long-read amplicon sequencing

Many vaccine candidate proteins in the malaria parasite Plasmodium falciparum are under strong immunological pressure and confer antigenic diversity. We present a sequencing and data analysis platform for the genomic surveillance of the insertion or deletion (indel)-rich antigens merozoite surface protein 1 (MSP1), MSP2, glutamate-rich protein (GLURP), and CSP from P. falciparum using long-read circular consensus sequencing (CCS) in multiclonal malaria isolates. Our platform uses 40 PCR primers per gene to asymmetrically barcode and identify multiclonal infections in pools of up to 384 samples. With msp2, we validated the method using 235 mock infections combining 10 synthetic variants at different concentrations and infection complexities. We applied this strategy to P. falciparum isolates from a longitudinal cohort in Tanzania. Finally, we constructed an analysis pipeline that streamlines the processing and interpretation of epidemiological and antigenic diversity data from demultiplexed FASTQ files. This platform can be easily adapted to other polymorphic antigens of interest in Plasmodium or any other human pathogen.

The adenomatous polyposis coli protein 30 years on

Mutations in the gene encoding the Adenomatous polyposis coli protein (APC) were discovered as driver mutations in colorectal cancers almost 30 years ago. Since then, the importance of APC in normal tissue homeostasis has been confirmed in a plethora of other (model) organisms spanning a large evolutionary space. APC is a multifunctional protein, with roles as a key scaffold protein in complexes involved in diverse signalling pathways, most prominently the Wnt signalling pathway. APC is also a cytoskeletal regulator with direct and indirect links to and impacts on all three major cytoskeletal networks. Correspondingly, a wide range of APC binding partners have been identified. Mutations in APC are extremely strongly associated with colorectal cancers, particularly those that result in the production of truncated proteins and the loss of significant regions from the remaining protein. Understanding the complement of its role in health and disease requires knowing the relationship between and regulation of its diverse functions and interactions. This in turn requires understanding its structural and biochemical features. Here we set out to provide a brief overview of the roles and function of APC and then explore its conservation and structure using the extensive sequence data, which is now available, and spans a broad range of taxonomy. This revealed conservation of APC across taxonomy and new relationships between different APC protein families.

Missense variants in human ACE2 strongly affect binding to SARS-CoV-2 Spike providing a mechanism for ACE2 mediated genetic risk in Covid-19: A case study in affinity predictions of interface variants

SARS-CoV-2 Spike (Spike) binds to human angiotensin-converting enzyme 2 (ACE2) and the strength of this interaction could influence parameters relating to virulence. To explore whether population variants in ACE2 influence Spike binding and hence infection, we selected 10 ACE2 variants based on affinity predictions and prevalence in gnomAD and measured their affinities and kinetics for Spike receptor binding domain through surface plasmon resonance (SPR) at 37°C. We discovered variants that reduce and enhance binding, including three ACE2 variants that strongly inhibited (p.Glu37Lys, ΔΔG = –1.33 ± 0.15 kcal mol^-1 and p.Gly352Val, predicted ΔΔG = –1.17 kcal mol^-1) or abolished (p.Asp355Asn) binding. We also identified two variants with distinct population distributions that enhanced affinity for Spike. ACE2 p.Ser19Pro (ΔΔG = 0.59 ± 0.08 kcal mol^-1) is predominant in the gnomAD African cohort (AF = 0.003) whilst p.Lys26Arg (ΔΔG = 0.26 ± 0.09 kcal mol^-1) is predominant in the Ashkenazi Jewish (AF = 0.01) and European non-Finnish (AF = 0.006) cohorts. We compared ACE2 variant affinities to published SARS-CoV-2 pseudotype infectivity data and confirmed that ACE2 variants with reduced affinity for Spike can protect cells from infection. The effect of variants with enhanced Spike affinity remains unclear, but we propose a mechanism whereby these alleles could cause greater viral spreading across tissues and cell types, as is consistent with emerging understanding regarding the interplay between receptor affinity and cell-surface abundance. Finally, we compared mCSM-PPI2 ΔΔG predictions against our SPR data to assess the utility of predictions in this system. We found that predictions of decreased binding were well-correlated with experiment and could be improved by calibration, but disappointingly, predictions of highly enhanced binding were unreliable. Recalibrated predictions for all possible ACE2 missense variants at the Spike interface were calculated and used to estimate the overall burden of ACE2 variants on Covid-19.

One of the first things the SARS-CoV-2 virus does to invade human cells is bind to a cell surface receptor called angiotensin-converting enzyme 2 (ACE2). The virus attaches to this receptor through its Spike protein and knowledge from other viruses tells us that the strength of this interaction influences how infectious and or virulent it is. We hypothesised that the Spike-ACE2 affinity might vary in people who have different amino acids in the part of ACE2 where Spike binds and consequently might be protected–or more at risk–from the virus. To test this idea, we measured the affinity of several ACE2 mutants, representing different versions found in humans, for the Spike protein and we found that some strengthened the interactions alongside others that weakened it. Most of these variants are rare, but two are present in over 1 in 1,000 individuals in certain populations and so might be important for the epidemiology of COVID-19. We then used computational methods to predict the affinity of even more ACE2 mutants than we could test in the lab and again found many that might alter this interaction. These data may help identify people who are at higher or lower risk from COVID-19.

Ankyrin repeats in context with human population variation

Ankyrin protein repeats bind to a wide range of substrates and are one of the most common protein motifs in nature. Here, we collate a high-quality alignment of 7,407 ankyrin repeats and examine for the first time, the distribution of human population variants from large-scale sequencing of healthy individuals across this family. Population variants are not randomly distributed across the genome but are constrained by gene essentiality and function. Accordingly, we interpret the population variants in context with evolutionary constraint and structural features including secondary structure, accessibility and protein-protein interactions across 383 three-dimensional structures of ankyrin repeats. We find five positions that are highly conserved across homologues and also depleted in missense variants within the human population. These positions are significantly enriched in intra-domain contacts and so likely to be key for repeat packing. In contrast, a group of evolutionarily divergent positions are found to be depleted in missense variants in human and significantly enriched in protein-protein interactions. Our analysis also suggests the domain has three, not two surfaces, each with different patterns of enrichment in protein-substrate interactions and missense variants. Our findings will be of interest to those studying or engineering ankyrin-repeat containing proteins as well as those interpreting the significance of disease variants.

The prefusion structure of herpes simplex virus glycoprotein B

Cell entry of enveloped viruses requires specialized viral proteins that mediate fusion with the host membrane by substantial structural rearrangements from a metastable pre- to a stable postfusion conformation. This metastability renders the herpes simplex virus 1 (HSV-1) fusion glycoprotein B (gB) highly unstable such that it readily converts into the postfusion form, thereby precluding structural elucidation of the pharmacologically relevant prefusion conformation. By identification of conserved sequence signatures and molecular dynamics simulations, we devised a mutation that stabilized this form. Functionally locking gB allowed the structural determination of its membrane-embedded prefusion conformation at sub-nanometer resolution and enabled the unambiguous fit of all ectodomains. The resulting pseudo-atomic model reveals a notable conservation of conformational domain rearrangements during fusion between HSV-1 gB and the vesicular stomatitis virus glycoprotein G, despite their very distant phylogeny. In combination with our comparative sequence-structure analysis, these findings suggest common fusogenic domain rearrangements in all class III viral fusion proteins.

Nucleophosmin in Its Interaction with Ligands

Nucleophosmin (NPM1) is a mainly nucleolar protein that shuttles between nucleoli, nucleoplasm and cytoplasm to fulfill its many functions. It is a chaperone of both nucleic acids and proteins and plays a role in cell cycle control, centrosome duplication, ribosome maturation and export, as well as the cellular response to a variety of stress stimuli. NPM1 is a hub protein in nucleoli where it contributes to nucleolar organization through heterotypic and homotypic interactions. Furthermore, several alterations, including overexpression, chromosomal translocations and mutations are present in solid and hematological cancers. Recently, novel germline mutations that cause dyskeratosis congenita have also been described. This review focuses on NPM1 interactions and inhibition. Indeed, the list of NPM1 binding partners is ever-growing and, in recent years, many studies contributed to clarifying the structural basis for NPM1 recognition of both nucleic acids and several proteins. Intriguingly, a number of natural and synthetic ligands that interfere with NPM1 interactions have also been reported. The possible role of NPM1 inhibitors in the treatment of multiple cancers and other pathologies is emerging as a new therapeutic strategy.

The majority of the matrix protein TapA is dispensable for Bacillus subtilis colony biofilm architecture

Biofilm formation is a co-operative behaviour, where microbial cells become embedded in an extracellular matrix. This biomolecular matrix helps manifest the beneficial or detrimental outcome mediated by the collective of cells. Bacillus subtilis is an important bacterium for understanding the principles of biofilm formation. The protein components of the B. subtilis matrix include the secreted proteins BslA, which forms a hydrophobic coat over the biofilm, and TasA, which forms protease-resistant fibres needed for structuring. TapA is a secreted protein also needed for biofilm formation and helps in vivo TasA-fibre formation but is dispensable for in vitro TasA-fibre assembly. We show that TapA is subjected to proteolytic cleavage in the colony biofilm and that only the first 57 amino acids of the 253-amino acid protein are required for colony biofilm architecture. Through the construction of a strain which lacks all eight extracellular proteases, we show that proteolytic cleavage by these enzymes is not a prerequisite for TapA function. It remains unknown why TapA is synthesised at 253 amino acids when the first 57 are sufficient for colony biofilm structuring; the findings do not exclude the core conserved region of TapA having a second role beyond structuring the B. subtilis colony biofilm.

The Dundee Resource for Sequence Analysis and Structure Prediction

The Dundee Resource for Sequence Analysis and Structure Prediction (DRSASP; http://www.compbio.dundee.ac.uk/drsasp.html) is a collection of web services provided by the Barton Group at the University of Dundee. DRSASP's flagship services are the JPred4 webserver for secondary structure and solvent accessibility prediction and the JABAWS 2.2 webserver for multiple sequence alignment, disorder prediction, amino acid conservation calculations, and specificity-determining site prediction. DRSASP resources are available through conventional web interfaces and APIs but are also integrated into the Jalview sequence analysis workbench, which enables the composition of multitool interactive workflows. Other existing Barton Group tools are being brought under the banner of DRSASP, including NoD (Nucleolar localization sequence detector) and 14-3-3-Pred. New resources are being developed that enable the analysis of population genetic data in evolutionary and 3D structural contexts. Existing resources are actively developed to exploit new technologies and maintain parity with evolving web standards. DRSASP provides substantial computational resources for public use, and since 2016 DRSASP services have completed over 1.5 million jobs.