Correlating protein function and stability through the analysis of single amino acid substitutions

Mutations resulting in the disruption of protein function are the underlying causes of many genetic diseases. Some mutations affect the number of expressed proteins while others alter the activity on a per-molecule basis. Single amino acid substitutions as caused by non-synonymous Single Nucleotide Polymorphisms (nsSNPs) often disrupt function by altering protein structure and/or stability, but can also wreak havoc by directly impacting functional binding sites. Given the experimental three-dimensional (3D) structure of a protein, we can try to differentiate between the "effect on structure/stability" and the "effect on binding". However, experimental 3D structures are available for only 1% of all known proteins; the magnitude of stability change caused by a given mutation is more widely available. Here, we analyze to which extent the functional effect of a mutation can be predicted from the effect on protein stability. We find that simple sequence-based methods succeed in predicting functional effects of nsSNPs. In fact, such methods consistently outperform approaches that predict functional change through the application of binary thresholds to stability change. We also observed that if stability is affected, functional change is easier to predict than when stability is not affected. Our results confirmed that stability change is somehow related to function change. However, we also show that the knowledge of stability changes in no way suffices to predict functional changes and that many function changing mutations have no effect on stability

Altering the stability of the Cdc8 overlap region modulates the ability of this tropomyosin to bind cooperatively to actin and regulate myosin.

Tropomyosin (Tm) is an evolutionarily conserved ?-helical coiled-coil protein, dimers of which form end-to-end polymers capable of associating with and stabilising actin-filaments and regulate myosin function. The fission yeast, Schizosaccharomyces pombe, possesses a single essential Tm, Cdc8, which can be acetylated on its amino terminal methionine to increase its affinity for actin and enhance its ability to regulate myosin function. We have designed and generated a number of novel Cdc8 mutant proteins with amino terminal substitutions to explore how stability of the Cdc8-polymer overlap region affects the regulatory function of this Tm. By correlating the stability of each protein, its propensity to form stable polymers, its ability to associate with actin and to regulate myosin, we have shown the stability of the amino terminal of the Cdc8 ?-helix is crucial for Tm function. In addition we have identified a novel Cdc8 mutant with increased amino-terminal stability, dimers of which are capable of forming Tm-polymers significantly longer than the wild-type protein. This protein had a reduced affinity for actin with respect to wild type, and was unable to regulate actomyosin interactions. The data presented here are consistent with acetylation providing a mechanism for modulating the formation and stability of Cdc8 polymers within the fission yeast cell. The data also provide evidence for a mechanism in which Tm dimers form end-to-end polymers on the actin-filament, consistent with a cooperative model for Tm binding to actin

Effect of bet missense mutations on bromodomain function, inhibitor binding and stability

Lysine acetylation is an important epigenetic mark regulating gene transcription and chromatin structure. Acetylated lysine residues are specifically recognized by bromodomains, small protein interaction modules that read these modification in a sequence and acetylation dependent way regulating the recruitment of transcriptional regulators and chromatin remodelling enzymes to acetylated sites in chromatin. Recent studies revealed that bromodomains are highly druggable protein interaction domains resulting in the development of a large number of bromodomain inhibitors. BET bromodomain inhibitors received a lot of attention in the oncology field resulting in the rapid translation of early BET bromodomain inhibitors into clinical studies. Here we investigated the effects of mutations present as polymorphism or found in cancer on BET bromodomain function and stability and the influence of these mutants on inhibitor binding. We found that most BET missense mutations localize to peripheral residues in the two terminal helices. Crystal structures showed that the three dimensional structure is not compromised by these mutations but mutations located in close proximity to the acetyl-lysine binding site modulate acetyl-lysine and inhibitor binding. Most mutations affect significantly protein stability and tertiary structure in solution, suggesting new interactions and an alternative network of protein-protein interconnection as a consequence of single amino acid substitution. To our knowledge this is the first report studying the effect of mutations on bromodomain function and inhibitor binding

Extreme genetic fragility of the HIV-1 capsid

Genetic robustness, or fragility, is defined as the ability, or lack thereof, of a biological entity to maintain function in the face of mutations. Viruses that replicate via RNA intermediates exhibit high mutation rates, and robustness should be particularly advantageous to them. The capsid (CA) domain of the HIV-1 Gag protein is under strong pressure to conserve functional roles in viral assembly, maturation, uncoating, and nuclear import. However, CA is also under strong immunological pressure to diversify. Therefore, it would be particularly advantageous for CA to evolve genetic robustness. To measure the genetic robustness of HIV-1 CA, we generated a library of single amino acid substitution mutants, encompassing almost half the residues in CA. Strikingly, we found HIV-1 CA to be the most genetically fragile protein that has been analyzed using such an approach, with 70% of mutations yielding replication-defective viruses. Although CA participates in several steps in HIV-1 replication, analysis of conditionally (temperature sensitive) and constitutively non-viable mutants revealed that the biological basis for its genetic fragility was primarily the need to coordinate the accurate and efficient assembly of mature virions. All mutations that exist in naturally occurring HIV-1 subtype B populations at a frequency >3%, and were also present in the mutant library, had fitness levels that were >40% of WT. However, a substantial fraction of mutations with high fitness did not occur in natural populations, suggesting another form of selection pressure limiting variation in vivo. Additionally, known protective CTL epitopes occurred preferentially in domains of the HIV-1 CA that were even more genetically fragile than HIV-1 CA as a whole. The extreme genetic fragility of HIV-1 CA may be one reason why cell-mediated immune responses to Gag correlate with better prognosis in HIV-1 infection, and suggests that CA is a good target for therapy and vaccination strategies

The loss and gain of functional amino acid residues is a common mechanism causing human inherited disease

Elucidating the precise molecular events altered by disease-causing genetic variants represents a major challenge in translational bioinformatics. To this end, many studies have investigated the structural and functional impact of amino acid substitutions. Most of these studies were however limited in scope to either individual molecular functions or were concerned with functional effects (e.g. deleterious vs. neutral) without specifically considering possible molecular alterations. The recent growth of structural, molecular and genetic data presents an opportunity for more comprehensive studies to consider the structural environment of a residue of interest, to hypothesize specific molecular effects of sequence variants and to statistically associate these effects with genetic disease. In this study, we analyzed data sets of disease-causing and putatively neutral human variants mapped to protein 3D structures as part of a systematic study of the loss and gain of various types of functional attribute potentially underlying pathogenic molecular alterations. We first propose a formal model to assess probabilistically function-impacting variants. We then develop an array of structure-based functional residue predictors, evaluate their performance, and use them to quantify the impact of disease-causing amino acid substitutions on catalytic activity, metal binding, macromolecular binding, ligand binding, allosteric regulation and post-translational modifications. We show that our methodology generates actionable biological hypotheses for up to 41% of disease-causing genetic variants mapped to protein structures suggesting that it can be reliably used to guide experimental validation. Our results suggest that a significant fraction of disease-causing human variants mapping to protein structures are function-altering both in the presence and absence of stability disruption

Path to Facilitate the Prediction of Functional Amino Acid Substitutions in Red Blood Cell Disorders – A Computational Approach

A major area of effort in current genomics is to distinguish mutations that are functionally neutral from those that contribute to disease. Single Nucleotide Polymorphisms (SNPs) are amino acid substitutions that currently account for approximately half of the known gene lesions responsible for human inherited diseases. As a result, the prediction of non-synonymous SNPs (nsSNPs) that affect protein functions and relate to disease is an important task.In this study, we performed a comprehensive analysis of deleterious SNPs at both functional and structural level in the respective genes associated with red blood cell metabolism disorders using bioinformatics tools. We analyzed the variants in Glucose-6-phosphate dehydrogenase (G6PD) and isoforms of Pyruvate Kinase (PKLR & PKM2) genes responsible for major red blood cell disorders. Deleterious nsSNPs were categorized based on empirical rule and support vector machine based methods to predict the impact on protein functions. Furthermore, we modeled mutant proteins and compared them with the native protein for evaluation of protein structure stability.We argue here that bioinformatics tools can play an important role in addressing the complexity of the underlying genetic basis of Red Blood Cell disorders. Based on our investigation, we report here the potential candidate SNPs, for future studies in human Red Blood Cell disorders. Current study also demonstrates the presence of other deleterious mutations and also endorses with in vivo experimental studies. Our approach will present the application of computational tools in understanding functional variation from the perspective of structure, expression, evolution and phenotype

Alpha-helical destabilization of the Bcl-2-BH4-domain peptide abolishes its ability to inhibit the IP3 receptor

The anti-apoptotic Bcl-2 protein is the founding member and namesake of the Bcl-2-protein family. It has recently been demonstrated that Bcl-2, apart from its anti-apoptotic role at mitochondrial membranes, can also directly interact with the inositol 1,4,5-trisphosphate receptor (IP3R), the primary Ca2+-release channel in the endoplasmic reticulum (ER). Bcl-2 can thereby reduce pro-apoptotic IP3R-mediated Ca2+ release from the ER. Moreover, the Bcl-2 homology domain 4 (Bcl-2-BH4) has been identified as essential and sufficient for this IP3R-mediated anti-apoptotic activity. In the present study, we investigated whether the reported inhibitory effect of a Bcl-2-BH4 peptide on the IP (3)R1 was related to the distinctive alpha-helical conformation of the BH4 domain peptide. We therefore designed a peptide with two glycine "hinges" replacing residues I14 and V15, of the wild-type Bcl-2-BH4 domain (Bcl-2-BH4-IV/GG). By comparing the structural and functional properties of the Bcl-2-BH4-IV/GG peptide with its native counterpart, we found that the variant contained reduced alpha-helicity, neither bound nor inhibited the IP (3)R1 channel, and in turn lost its anti-apoptotic effect. Similar results were obtained with other substitutions in Bcl-2-BH4 that destabilized the alpha-helix with concomitant loss of IP3R inhibition. These results provide new insights for the further development of Bcl-2-BH4-derived peptides as specific inhibitors of the IP3R with significant pharmacological implications

Physicochemical properties of pore residues predict activation gating of CaV1.2: A correlation mutation analysis

Single point mutations in pore-forming S6 segments of calcium channels may transform a high-voltage-activated into a low-voltage-activated channel, and resulting disturbances in calcium entry may cause channelopathies (Hemara-Wahanui et al., Proc Natl Acad Sci U S A 102(21):7553–7558, 16). Here we ask the question how physicochemical properties of amino acid residues in gating-sensitive positions on S6 segments determine the threshold of channel activation of CaV1.2. Leucine in segment IS6 (L434) and a newly identified activation determinant in segment IIIS6 (G1193) were mutated to a variety of amino acids. The induced leftward shifts of the activation curves and decelerated current activation and deactivation suggest a destabilization of the closed and a stabilisation of the open channel state by most mutations. A selection of 17 physicochemical parameters (descriptors) was calculated for these residues and examined for correlation with the shifts of the midpoints of the activation curve (ΔVact). ΔVact correlated with local side-chain flexibility in position L434 (IS6), with the polar accessible surface area of the side chain in position G1193 (IIIS6) and with hydrophobicity in position I781 (IIS6). Combined descriptor analysis for positions I781 and G1193 revealed that additional amino acid properties may contribute to conformational changes during the gating process. The identified physicochemical properties in the analysed gating-sensitive positions (accessible surface area, side-chain flexibility, and hydrophobicity) predict the shifts of the activation curves of CaV1.2