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.

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.

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.

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.

As the COVID-19 pandemic has progressed, new variants of the virus SARS-CoV-2 have emerged that are more infectious than the original form. The variants known as Alpha, Beta and Gamma have mutations in a protein on the virus's surface that is vital for attaching to cells and infecting them. This protein, called Spike, carries out its role by binding to ACE2, a protein on the surface of human cells. Mutations on Spike are found on the region where it binds to ACE2. The interaction between these two proteins appears to be important to the behaviour of SARS-CoV-2, but the impact of individual mutations in Spike is unknown. In addition, some people have different variants of ACE2 with mutations in the region that interacts with Spike, but it is not known whether this affects these people's risk of contracting COVID-19. To answer these questions, Barton et al. measured the precise effect of mutations in Spike and ACE2 on the strength of the interaction between the two proteins. The experiments showed that three of the five common Spike mutations in the Alpha, Beta and Gamma SARS-CoV-2 variants strengthened binding to ACE2. The two mutations that weakened binding were only found together with other mutations that strengthened binding. This meant that the Spike proteins in all three of these SARS-CoV-2 variants bind to ACE2 more strongly than the original form. The experiments also showed that two common variants of ACE2 also increased the strength of binding to Spike. Interestingly, one of these ACE2 variants reversed the effect of a specific SARS-CoV-2 mutation, suggesting that carriers would be resistant to SARS-CoV-2 variants with this mutation. Identifying the precise effects of Spike mutations on ACE2 binding helps understand why new variants of SARS-CoV-2 spread more rapidly. This could help to identify concerning new variants before they spread widely and inform the response by health authorities. The finding that two common ACE2 variants bind more strongly to Spike suggests that people with these mutations could be more susceptible to SARS-CoV-2.

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 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.

Disease related single point mutations alter the global dynamics of a tetratricopeptide (TPR) α-solenoid domain.

Tetratricopeptide repeat (TPR) proteins belong to the class of α-solenoid proteins, in which repetitive units of α-helical hairpin motifs stack to form superhelical, often highly flexible structures. TPR domains occur in a wide variety of proteins, and perform key functional roles including protein folding, protein trafficking, cell cycle control and post-translational modification. Here, we look at the TPR domain of the enzyme O-linked GlcNAc-transferase (OGT), which catalyses O-GlcNAcylation of a broad range of substrate proteins. A number of single-point mutations in the TPR domain of human OGT have been associated with the disease Intellectual Disability (ID). By extended steered and equilibrium atomistic simulations, we show that the OGT-TPR domain acts as an elastic nanospring, and that each of the ID-related local mutations substantially affect the global dynamics of the TPR domain. Since the nanospring character of the OGT-TPR domain is key to its function in binding and releasing OGT substrates, these changes of its biomechanics likely lead to defective substrate interaction. We find that neutral mutations in the human population, selected by analysis of the gnomAD database, do not incur these changes. Our findings may not only help to explain the ID phenotype of the mutants, but also aid the design of TPR proteins with tailored biomechanical properties.

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

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.

Contribution of bacteriochlorophyll conformation to the distribution of site-energies in the FMO protein.

The structural data for the Fenna-Matthews-Olson (FMO) protein indicate that the bacteriochlorophylls (BChls) display a significant degree of conformational heterogeneity of their peripheral substituents and the protein-induced nonplanar skeletal deformations of the tetrapyrrole macrocycle. As electronic properties of chromophores are altered by such differences, a conformational effect may influence the site-energies of specific pigments and thus play a role in mediating the excitation energy transfer dynamics, but this has not yet been established. The difficulty of assessing this question is shown to be partly the result of the inability of the sequential truncation approach usually employed to account for interactions between the conformations of the macrocycle and its substituents and an alternative approach is suggested. By assigning the BChl atoms to meaningful atom groups and performing all possible permutations of partial optimizations in a full-factorial design, where each group is either frozen in the crystal geometry or optimized in vacuo, followed by excited state calculations on each resulting structure (PM6//ZIndo/S), the specific effects of the conformations of each BChl component as well as mutual interactions between the molecular fragments on the site-energy can be delineated. This factorial relaxation procedure gives different estimates of the macrocycle conformational perturbation than the approach of sequentially truncating the BChl periphery. The results were evaluated in the context of published site-energies for the FMO pigments from three species to identify how conformational effects contribute to their distribution and instances of cross-species conservation and functional divergence of the BChl nonplanarity conformational contribution are described.

Conformational control of cofactors in nature - the influence of protein-induced macrocycle distortion on the biological function of tetrapyrroles.

Tetrapyrrole-containing proteins are one of the most fundamental classes of enzymes in nature and it remains an open question to give a chemical rationale for the multitude of biological reactions that can be catalyzed by these pigment-protein complexes. There are many fundamental processes where the same (i.e., chemically identical) porphyrin cofactor is involved in chemically quite distinct reactions. For example, heme is the active cofactor for oxygen transport and storage (hemoglobin, myoglobin) and for the incorporation of molecular oxygen in organic substrates (cytochrome P450). It is involved in the terminal oxidation (cytochrome c oxidase) and the metabolism of H2O2 (catalases and peroxidases) and catalyzes various electron transfer reactions in cytochromes. Likewise, in photosynthesis the same chlorophyll cofactor may function as a reaction center pigment (charge separation) or as an accessory pigment (exciton transfer) in light harvesting complexes (e.g., chlorophyll a). Whilst differences in the apoprotein sequences alone cannot explain the often drastic differences in physicochemical properties encountered for the same cofactor in diverse protein complexes, a critical factor for all biological functions must be the close structural interplay between bound cofactors and the respective apoprotein in addition to factors such as hydrogen bonding or electronic effects. Here, we explore how nature can use the same chemical molecule as a cofactor for chemically distinct reactions using the concept of conformational flexibility of tetrapyrroles. The multifaceted roles of tetrapyrroles are discussed in the context of the current knowledge on distorted porphyrins. Contemporary analytical methods now allow a more quantitative look at cofactors in protein complexes and the development of the field is illustrated by case studies on hemeproteins and photosynthetic complexes. Specific tetrapyrrole conformations are now used to prepare bioengineered designer proteins with specific catalytic or photochemical properties.

Expanding the clinical spectrum of hereditary fibrosing poikiloderma with tendon contractures, myopathy and pulmonary fibrosis due to FAM111B mutations.

BACKGROUND: Hereditary Fibrosing Poikiloderma (HFP) with tendon contractures, myopathy and pulmonary fibrosis (POIKTMP [MIM 615704]) is a very recently described entity of syndromic inherited poikiloderma. Previously by using whole exome sequencing in five families, we identified the causative gene, FAM111B (NM_198947.3), the function of which is still unknown. Our objective in this study was to better define the specific features of POIKTMP through a larger series of patients. METHODS: Clinical and molecular data of two families and eight independent sporadic cases, including six new cases, were collected. RESULTS: Key features consist of: (i) early-onset poikiloderma, hypotrichosis and hypohidrosis; (ii) multiple contractures, in particular triceps surae muscle contractures; (iii) diffuse progressive muscular weakness; (iv) pulmonary fibrosis in adulthood and (v) other features including exocrine pancreatic insufficiency, liver impairment and growth retardation. Muscle magnetic resonance imaging was informative and showed muscle atrophy and fatty infiltration. Histological examination of skeletal muscle revealed extensive fibroadipose tissue infiltration. Microscopy of the skin showed a scleroderma-like aspect with fibrosis and alterations of the elastic network. FAM111B gene analysis identified five different missense variants (two recurrent mutations were found respectively in three and four independent families). All the mutations were predicted to localize in the trypsin-like cysteine/serine peptidase domain of the protein. We suggest gain-of-function or dominant-negative mutations resulting in FAM111B enzymatic activity changes. CONCLUSIONS: HFP with tendon contractures, myopathy and pulmonary fibrosis, is a multisystemic disorder due to autosomal dominant FAM111B mutations. Future functional studies will help in understanding the specific pathological process of this fibrosing disorder.

Severe dermatitis, multiple allergies, and metabolic wasting syndrome caused by a novel mutation in the N-terminal plakin domain of desmoplakin.

BACKGROUND: Severe dermatitis, multiple allergies, and metabolic wasting (SAM) syndrome is a recently recognized syndrome caused by mutations in the desmoglein 1 gene (DSG1). To date, only 3 families have been reported. OBJECTIVE: We studied a new case of SAM syndrome known to have no mutations in DSG1 to detail the clinical, histopathologic, immunofluorescent, and ultrastructural phenotype and to identify the underlying molecular mechanisms in this rare genodermatosis. METHODS: Histopathologic, electron microscopy, and immunofluorescent studies were performed. Whole-exome sequencing data were interrogated for mutations in desmosomal and other skin structural genes, followed by Sanger sequencing of candidate genes in the patient and his parents. RESULTS: No mutations were identified in DSG1; however, a novel de novo heterozygous missense c.1757A>C mutation in the desmoplakin gene (DSP) was identified in the patient, predicting the amino acid substitution p.His586Pro in the desmoplakin polypeptide. CONCLUSIONS: SAM syndrome can be caused by mutations in both DSG1 and DSP. Knowledge of this genetic heterogeneity is important for both analysis of patients and genetic counseling of families. This condition and these observations reinforce the importance of heritable skin barrier defects, in this case desmosomal proteins, in the pathogenesis of atopic disease.

Computational quantification of the physicochemical effects of heme distortion: redox control in the reaction center cytochrome subunit of Blastochloris viridis.

A facile, experimentally calibrated computational procedure is described that affords the relative ordering of heme cofactor reduction potentials with respect to intrinsic shifts brought about by apoprotein induced heme-macrocycle distortion. The method utilizes heme-Fe partial atomic charges and is useful with the computationally inexpensive B3LYP/3-21g method calculated for simplified heme models extracted from the Protein Data Bank incorporating only the effects of varying macrocycle conformations and thereby delineating their physicochemical effects. The procedure was successfully calibrated using the atomic coordinates and published midpoint potentials from the heme cofactors in wild-type and a series of heme-NO and -O(2) binding domain mutants and thus confirmed the sole conformational modulation of the redox potentials in these complexes. This technique was also applied to the reaction center tetraheme cytochrome subunit of Blastochloris viridis to build upon previous work elucidating the role that conformational control plays in photosynthetic systems, and it was found that this effect may account for up to 70% (54mv) of the observed differences in the reduction potentials of the four hemes. We validate the approach using larger basis sets up to and including the triple-ζ, doubly polarized and augmented 6-311+g** basis and discuss the specific conformational origins of the effect.

Conformational control of cofactors in nature--functional tetrapyrrole conformations in the photosynthetic reaction centers of purple bacteria.

Chemically identical tetrapyrrole cofactors such as hemes and chlorophylls participate in functionally diverse biological roles. An analysis of the available protein structural data for the bacteriochlorophylls in the photosynthetic reaction center gives statistically reliable evidence of the hypothesis that the protein induced cofactor conformation is a modulator of the bio-molecular function of each reaction center. The results serve as a general model to illustrate conformational control of tetrapyrrole cofactors in other proteins.