RESUMO
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.
Assuntos
Proteínas , Humanos , Domínios Proteicos , Proteínas/metabolismo , Ligação Proteica , Sequência de BasesRESUMO
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.
Assuntos
Descoberta de Drogas , Proteínas , Humanos , Ligantes , Sítios de Ligação , Proteínas/metabolismo , Descoberta de Drogas/métodosRESUMO
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.
Assuntos
Citostáticos , Tripanossomicidas , Trypanosoma brucei brucei , Tripanossomíase Africana , Animais , Humanos , Tripanossomíase Africana/tratamento farmacológico , Tripanossomicidas/uso terapêutico , Tripanossomicidas/farmacocinética , Citostáticos/uso terapêutico , Barreira HematoencefálicaRESUMO
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.
Assuntos
Enzima de Conversão de Angiotensina 2/genética , COVID-19/genética , Mutação de Sentido Incorreto , Glicoproteína da Espícula de Coronavírus/metabolismo , Enzima de Conversão de Angiotensina 2/metabolismo , Predisposição Genética para Doença , Humanos , Ligação ProteicaRESUMO
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.
Assuntos
Enzima de Conversão de Angiotensina 2/metabolismo , Mutação , SARS-CoV-2/genética , Glicoproteína da Espícula de Coronavírus/genética , Glicoproteína da Espícula de Coronavírus/metabolismo , Enzima de Conversão de Angiotensina 2/química , Sítios de Ligação , COVID-19/virologia , Humanos , Cinética , Ligantes , Simulação de Dinâmica Molecular , Ligação Proteica , Domínios Proteicos , Domínios e Motivos de Interação entre Proteínas , SARS-CoV-2/classificação , SARS-CoV-2/metabolismo , Glicoproteína da Espícula de Coronavírus/químicaRESUMO
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.
Assuntos
Repetição de Anquirina , Variação Genética , Humanos , Interações Hidrofóbicas e Hidrofílicas , Modelos Moleculares , Mutação de Sentido Incorreto , Ligação Proteica , Proteínas/química , Proteínas/genéticaRESUMO
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.
Assuntos
Biologia Computacional/métodos , Proteínas/química , Análise de Sequência de Proteína/métodos , Estrutura Secundária de Proteína , Alinhamento de Sequência , Software , NavegadorRESUMO
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.
Assuntos
Deficiência Intelectual/genética , N-Acetilglucosaminiltransferases/química , N-Acetilglucosaminiltransferases/genética , Mutação Puntual , Humanos , Simulação de Dinâmica Molecular , N-Acetilglucosaminiltransferases/metabolismo , Conformação Proteica , Domínios Proteicos , Repetições de TetratricopeptídeosRESUMO
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.
Assuntos
Proteínas de Bactérias/química , Bacterioclorofilas/química , Complexos de Proteínas Captadores de Luz/química , Cristalização , Conformação MolecularRESUMO
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.
Assuntos
Proteínas/química , Tetrapirróis/química , Modelos Moleculares , Conformação MolecularRESUMO
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.
Assuntos
Citocromos/química , Heme/química , Proteobactérias/enzimologia , Teoria Quântica , Físico-Química , Citocromos/metabolismo , Heme/metabolismo , Modelos Moleculares , Oxirredução , Proteobactérias/metabolismoRESUMO
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.