A unified analysis of evolutionary and population constraint in protein domains highlights structural features and pathogenic sites

Protein evolution is constrained by structure and function, creating patterns in residue conservation that are routinely exploited to predict structure and other features. Similar constraints should affect variation across individuals, but it is only with the growth of human population sequencing that this has been tested at scale. Now, human population constraint has established applications in pathogenicity prediction, but it has not yet been explored for structural inference. Here, we map 2.4 million population variants to 5885 protein families and quantify residue-level constraint with a new Missense Enrichment Score (MES). Analysis of 61,214 structures from the PDB spanning 3661 families shows that missense depleted sites are enriched in buried residues or those involved in small-molecule or protein binding. MES is complementary to evolutionary conservation and a combined analysis allows a new classification of residues according to a conservation plane. This approach finds functional residues that are evolutionarily diverse, which can be related to specificity, as well as family-wide conserved sites that are critical for folding or function. We also find a possible contrast between lethal and non-lethal pathogenic sites, and a surprising clinical variant hot spot at a subset of missense enriched positions

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

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

Effects of common mutations in the SARS-CoV-2 Spike RBD and its ligand the human ACE2 receptor on binding affinity and kinetics

The interaction between the SARS-CoV-2 virus Spike protein receptor binding domain (RBD) and the ACE2 cell surface protein is required for viral infection of cells. Mutations in the RBD are present in SARS-CoV-2 variants of concern that have emerged independently worldwide. For example, the B.1.1.7 lineage has a mutation (N501Y) in its Spike RBD that enhances binding to ACE2. There are also ACE2 alleles in humans with mutations in the RBD binding site. Here we perform a detailed affinity and kinetics analysis of the effect of five common RBD mutations (K417N, K417T, N501Y, E484K, and S477N) and two common ACE2 mutations (S19P and K26R) on the RBD/ACE2 interaction. We analysed the effects of individual RBD mutations and combinations found in new SARS-CoV-2 Alpha (B.1.1.7), Beta (B.1.351), and Gamma (P1) variants. Most of these mutations increased the affinity of the RBD/ACE2 interaction. The exceptions were mutations K417N/T, which decreased the affinity. Taken together with other studies, our results suggest that the N501Y and S477N mutations enhance transmission primarily by enhancing binding, the K417N/T mutations facilitate immune escape, and the E484K mutation enhances binding and immune escape

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

Development of a 2,4-diaminothiazole series for the treatment of human African trypanosomiasis highlights the importance of static-cidal screening of analogues

While treatment options for human African trypanosomiasis (HAT) have improved significantly, there is still a need for new drugs with eradication now a realistic possibility. Here, we report the development of 2,4-diaminothiazoles that demonstrate significant potency against Trypanosoma brucei, the causative agent of HAT. Using phenotypic screening to guide structure-activity relationships, potent drug-like inhibitors were developed. Proof of concept was established in an animal model of the hemolymphatic stage of HAT. To treat the meningoencephalitic stage of infection, compounds were optimized for pharmacokinetic properties, including blood-brain barrier penetration. However, in vivo efficacy was not achieved, in part due to compounds evolving from a cytocidal to a cytostatic mechanism of action. Subsequent studies identified a nonessential kinase involved in the inositol biosynthesis pathway as the molecular target of these cytostatic compounds. These studies highlight the need for cytocidal drugs for the treatment of HAT and the importance of static-cidal screening of analogues

Mutations involving the SRY-related gene SOX8 are associated with a spectrum of human reproductive anomalies.

© The Author(s) 2018. Published by Oxford University Press. All rights reserved. SOX8 is an HMG-box transcription factor closely related to SRY and SOX9. Deletion of the gene encoding Sox8 in mice causes reproductive dysfunction but the role of SOX8 in humans is unknown. Here, we show that SOX8 is expressed in the somatic cells of the early developing gonad in the human and influences human sex determination. We identified two individuals with 46, XY disorders/differences in sex development (DSD) and chromosomal rearrangements encompassing the SOX8 locus and a third individual with 46, XY DSD and a missense mutation in the HMG-box of SOX8. In vitro functional assays indicate that this mutation alters the biological activity of the protein. As an emerging body of evidence suggests that DSDs and infertility can have common etiologies, we also analysed SOX8 in a cohort of infertile men (n=274) and two independent cohorts of women with primary ovarian insufficiency (POI; n=153 and n=104). SOX8 mutations were found at increased frequency in oligozoospermic men (3.5%; P < 0.05) and POI (5.06%; P=4.5×10 -5 ) as compared with fertile/normospermic control populations (0.74%). The mutant proteins identified altered SOX8 biological activity as compared with the wild-type protein. These data demonstrate that SOX8 plays an important role in human reproduction and SOX8 mutations contribute to a spectrum of phenotypes including 46, XY DSD, male infertility and 46, XX POI.Link_to_subscribed_fulltex