BPPred: a Web-based computational tool for predicting biophysical parameters of proteins

Bioactivity of dihydropyrimidinone derivatives as inhibitors of cyclooxygenase-2 (COX-2): an in silico approach

Cyclooxygenase-2 (COX-2) is an enzyme involved in inflammation. The overexpression of COX-2 causes chronic inflammation, which can be prevented by COX-2 inhibitors. Generally, COX-2 inhibitors possess a carboxyl group and an aromatic ring in their molecular structure. These moieties are involved in the interaction with the active site of COX-2, thus playing a pivotal role in the inhibitory activity. Regarding the requisite molecular structure of COX-2 inhibitors, derivatives of dihydropyrimidinone (DHPM) are ideal candidates to be explored as COX-2 inhibitors, due to the ease of synthesis and their versatility to be transformed chemically. In this study, we prepared a novel small library consisting of 288 designed DHPM derivatives by varying the constituent components. The selection criteria of potential candidates for the COX-2 inhibitor of the data bank involve in silico studies via molecular docking investigations, prediction of ADMET and druglikeness, as well as molecular dynamics (MD) simulations. Molecular docking served as the initial step of selection, based on the comparison of grid score, docking pose, and interactions with those of lumiracoxib (LUR) as the original ligand of COX-2. The next criteria of selection were scores obtained from the ADMET and druglikeness by comparing the designed candidates with COX-2 inhibitors that were already marketed. Compound RDUE2 and SDT29 were the most potential candidates, which were further analyzed using the MD simulation. The results of the MD simulation indicated that RDUE2 and SDT29 interacted stably with amino acid residues on the active site of COX-2. The estimation of binding free energy indicated that SDT29 exhibited an inhibitory activity comparable to that of LUR, whereas RDUE2 showed a lower inhibitory activity than that of SDT29 and LUR.

In silico design of novel CDK2 inhibitors through QSAR, ADMET, molecular docking and molecular dynamics simulation studies

The present study aims to investigate about the quantitative structure-activity relationship (QSAR) of a series of Thiazole derivatives reported as anticancer agents (hepatocellular carcinoma), using principally the electronic descriptors calculated by the DFT method and by applying the multiple linear regression method. The developed model showed good statistical parameters (R2 = 0.725, R2adj = 0.653, MSE = 0.060, R2test = 0.827, Q2cv = 0.536). The energy EHOMO orbital, electronic energy (TE), shape coefficient (I), number of rotatable bonds (NROT), and index of refraction (n) were revealed to be the main descriptors influencing the anti-cancer activity. Further, new Thiazole derivatives have been designed and their activities and pharmacokinetic properties have been predicted using the validated QSAR model. The designed molecules were then assessed to molecular docking (MD), and molecular dynamic (MDs) simulation accompanied by the calculation of the binding affinity using MMPBSA script according to 100 ns a simulation trajectory, to study both their affinity and their stability towards CDK2 as a target protein for the cancer disease treatment. This research concluded with the identification of four new CDK2 inhibitors which are A1, A3, A5, and A6 showing good pharmacokinetic properties. The MDs results revealed that the newly designed compound A5 remained stable in the active center of the discovered CDK2 protein, indicating its potential as a novel inhibitor for the treatment of hepatocellular carcinoma. The current findings may eventually contribute to the development of robust CDK2 inhibitors in the future.Communicated by Ramaswamy H. Sarma.

Mechanism exploration and prognosis study of Astragali Radix-Spreading hedyotis herb for the treatment of lung adenocarcinoma based on bioinformatics approaches and molecular dynamics simulation

Background: Herb pair of Astragali Radix (AR) and Spreading Hedyotis Herb (SH) has been frequently prescribed in clinical for the treatment of lung cancer owing to its favorable efficacy. Yet, the mechanism under the therapeutic effects remained unveiled, which has limited its clinical applications, and new drug development for lung cancer.

Methods: The bioactive ingredients of AR and SH were retrieved from the Traditional Chinese Medicine System Pharmacology Database, with the targets of obtained components predicted by Swiss Target Prediction. Genes related to lung adenocarcinoma (LUAD) were acquired from GeneCards, OMIM and CTD databases, with the hub genes of LUAD screened by CTD database. The intersected targets of LUAD and AR-SH were obtained by Venn, with David Database employed to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Survival analysis of the hub genes of LUAD was carried out using TCGA-LUAD dataset. Molecular docking of core proteins and active ingredients was performed by Auto-Dock Vina software, followed by molecular dynamics simulations of protein-ligand complexes with well-docked conformations.

Results: 29 active ingredients were screened out with 422 corresponding targets predicted. It is revealed that AR-SH can act on various targets such as EGFR, MAPK1, and KARS by ursolic acid (UA), Astragaloside IV(ASIV), and Isomucronulatol 7,2′-di-O-glucoside (IDOG) to alleviate the symptoms of LUAD. Biological processes involved are protein phosphorylation, negative regulation of apoptotic process, and pathways involved are endocrine resistance, EGFR tyrosine kinase inhibitor resistance, PI3K-Akt, and HIF-1 pathway. Molecular docking analysis indicated that the binding energy of most of the screened active ingredients to proteins encoded by core genes was less than −5.6 kcal/mol, with some active ingredients showing even lower binding energy to EGFR than Gefitinib. Three ligand-receptor complexes including EGFR-UA, MAPK1-ASIV, and KRAS-IDOG were found to bind relatively stable by molecular dynamics simulation, which was consistent with the results of molecule docking.

Conclusion: We suggested that the herb pair of AR-SH can act on targets like EGFR, MAPK1 and KRAS by UA, ASIV and IDOG, to play a vital role in the treatment and the enhancement of prognosis of LUAD.

Exploring the Folding Mechanism of Dimeric Superoxide Dismutase

The Cu/Zn Human Superoxide Dismutase (SOD1) is a dimeric metalloenzyme whose genetic mutations are directly related to amyotrophic lateral sclerosis (ALS), so understanding its folding mechanism is of fundamental importance. Currently, the SOD1 dimer formation is studied via molecular dynamics simulations using a simplified structure-based model and an all-atom model. Results from the simplified model reveal a mechanism dependent on distances between monomers, which are limited by constraints to mimic concentration dependence. The stability of intermediates (during the int state) is significantly affected by this distance, as well as by the presence of two folded monomers prior to dimer formation. The kinetics of interface formation are also highly dependent on the separation distance. The folding temperature of the dimer is about 4.2% higher than that of the monomer, a value not too different from experimental data. All-atom simulations on the apo dimer give binding free energy between monomers similar to experimental values. An intermediate state is evident for the apo form at a separation distance between monomers slightly larger than the native distance which has little formed interface between monomers. We have shown that this intermediate is stabilized by non-native intra- and intercontacts.

Cryptic binding properties of a transient folding intermediate in a PDZ tandem repeat

PDZ domains are the most diffused protein-protein interaction modules of the human proteome and are often present in tandem repeats. An example is PDZD2, a protein characterized by the presence of six PDZ domains that undergoes a proteolytic cleavage producing sPDZD2, comprising a tandem of two PDZ domains, namely PDZ5 and PDZ6. Albeit the physiopathological importance of sPDZD2 is well-established, the interaction with endogenous ligands has been poorly characterized. To understand the determinants of the stability and function of sPDZD2, we investigated its folding pathway. Our data highlights the presence of a complex scenario involving a transiently populated folding intermediate that may be accumulated from the concurrent denaturation of both PDZ5 and PDZ6 domains. Importantly, double jump kinetic experiments allowed us to pinpoint the ability of this transient intermediate to bind the physiological ligand of sPDZD2 with increased affinity compared to the native state. In summary, our results provide an interesting example of a functionally competent misfolded intermediate, which may exert a cryptic function that is not captured from the analysis of the native state only.

Natural flavonoids effectively block the CD81 receptor of hepatocytes and inhibit HCV infection: a computational drug development approach

Hepatitis C virus (HCV) infection is a major public health concern, and almost two million people are infected per year globally. This is occurred by the diverse spectrum of viral genotypes, which are directly associated with chronic liver disease (fibrosis, and cirrhosis). Indeed, the viral genome encodes three principal proteins as sequentially core, E1, and E2. Both E1 and E2 proteins play a crucial role in the attachment of the host system, but E2 plays a more fundamental role in attachment. The researchers have found the "E2-CD81 complex" at the entry site, and therefore, CD81 is the key receptor for HCV entrance in both humans, and chimpanzees. So, the researchers are trying to block the host CD81 receptor and halt the virus entry within the cellular system via plant-derived compounds. Perhaps that is why the current research protocol is designed to perform an in silico analysis of the flavonoid compounds for targeting the tetraspanin CD81 receptor of hepatocytes. To find out the best flavonoid compounds from our library, web-based tools (Swiss ADME, pKCSM), as well as computerized tools like the PyRx, PyMOL, BIOVIA Discovery Studio Visualizer, Ligplot+ V2.2, and YASARA were employed. For molecular docking studies, the flavonoid compounds docked with the targeted CD81 protein, and herein, the best-outperformed compounds are Taxifolin, Myricetin, Puerarin, Quercetin, and (-)-Epicatechin, and outstanding binding affinities are sequentially - 7.5, - 7.9, - 8.2, - 8.4, and - 8.5 kcal/mol, respectively. These compounds have possessed more interactions with the targeted protein. To validate the post docking data, we analyzed both 100 ns molecular dynamic simulation, and MM-PBSA via the YASARA simulator, and finally finds the more significant outcomes. It is concluded that in the future, these compounds may become one of the most important alternative antiviral agents in the fight against HCV infection. It is suggested that further in vivo, and in vitro research studies should be done to support the conclusions of this in silico research workflow.

Target specificity of selective bioactive compounds in blocking α-dystroglycan receptor to suppress Lassa virus infection: an in silico approach

Lassa hemorrhagic fever, caused by Lassa mammarenavirus (LASV) infection, accumulates up to 5000 deaths every year. Currently, there is no vaccine available to combat this disease. In this study, a library of 200 bioactive compounds was virtually screened to study their drug-likeness with the capacity to block the α-dystroglycan (α-DG) receptor and prevent LASV influx. Following rigorous absorption, distribution, metabolism, and excretion (ADME) and quantitative structure-activity relationship (QSAR) profiling, molecular docking was conducted with the top ligands against the α-DG receptor. The compounds chrysin, reticuline, and 3-caffeoylshikimic acid emerged as the top three ligands in terms of binding affinity. Post-docking analysis revealed that interactions with Arg76, Asn224, Ser259, and Lys302 amino acid residues of the receptor protein were important for the optimum binding affinity of ligands. Molecular dynamics simulation was performed comprehensively to study the stability of the protein-ligand complexes. In-depth assessment of root-mean-square deviation (RMSD), root mean square fluctuation (RMSF), polar surface area (PSA), B-Factor, radius of gyration (Rg), solvent accessible surface area (SASA), and molecular surface area (MolSA) values of the protein-ligand complexes affirmed that the candidates with the best binding affinity formed the most stable protein-ligand complexes. To authenticate the potentialities of the ligands as target-specific drugs, an in vivo study is underway in real time as the continuation of the research.

Small Neutral Crowding Solute Effects on Protein Folding Thermodynamic Stability and Kinetics

Molecular crowding is a ubiquitous phenomenon in biological systems, with significant consequences on protein folding and stability. Small compounds, such as the osmolyte trimethylamine N-oxide (TMAO), can also present similar effects. To analyze the effects arising from these solute-like molecules, we performed a series of crowded coarse-grained folding simulations. Two well-known proteins were chosen, CI2 and SH3, modeled by the alpha-carbon-structure-based model. In the simulations, the crowding molecules were represented by low-sized neutral atom beads in different concentrations. The results show that a low level of the volume fraction occupied by neutral agents can change protein stability and folding kinetics for the two systems. However, the kinetics were shown to be unaffected in their respective folding temperatures. The faster kinetics correlates with changes in the folding route for systems at the same temperature, where non-native contacts were shown to be relevant. The transition states of the two systems with and without crowders are similar. In the case of SH3, there are differences in the structuring of two strands, which may be associated with the increase in its folding rate, in addition to the destabilization of the denatured ensemble. The present study also detected a crossover in the thermodynamic stability behavior, previously observed experimentally and theoretically. As the temperature increases, crowders change from destabilizing to stabilizing agents.

Computational and Experimental Assessment of Backbone Templates for Computational Redesign of the Thioredoxin Fold

Computational protein design has taken big strides in recent years; however, the tools available are still not at a state where a sequence can be designed to fold into a given protein structure at will and with high probability. We have applied here a recent release of Rosetta Design to redesign a set of structurally very similar proteins belonging to the thioredoxin fold. We used a genetic screening tool to estimate solubility/folding of the designed proteins in E. coli and to select the best hits from this for further biochemical characterization. We have previously used this set of template proteins for redesign and found that success was highly dependent on template structure, a trait which was also found in this study. Nevertheless, state-of-the-art design software is now able to predict the best template, most likely due to the introduction of an energy term that reports on stress in covalent bond lengths and angles. The template that led to the greatest fraction of successful designs was the same (a thioredoxin from spinach) as that identified in our previous study. Our previously described redesign of thioredoxin, which also used the spinach protein as a template, however also performed well as a template. In the present study, both of these templates yielded proteins with compact folded structures and enforced the conclusion that any design project must carefully consider different design templates. Fortunately, selecting designs based on energies appears to correctly identify such templates.

Dissecting Inter-domain Cooperativity in the Folding of a Multi Domain Protein

Correct protein folding underlies all cellular functions. While there are detailed descriptions and a good understanding of protein folding pathways for single globular domains there is a paucity of quantitative data regarding folding of multidomain proteins. We have here investigated the folding of a three-domain supramodule from the protein PSD-95, consisting of one PDZ domain, one SH3 domain and one guanylate kinase-like (GK) domain. This supramodule has previously been shown to work as one functional unit with regard to ligand binding. We used equilibrium and kinetic folding experiments to demonstrate that the PDZ domain folds faster and independently from the SH3-GK tandem, which folds as one cooperative unit. However, concurrent folding of the PDZ domain slows down folding of SH3-GK by non-native interactions, resulting in an off-pathway folding intermediate. Our data contribute to an emerging description of multidomain protein folding in which individual domains cannot a priori be viewed as separate folding units.

Linking thermodynamics and measurements of protein stability

We review the background, theory and general equations for the analysis of equilibrium protein unfolding experiments, focusing on denaturant and heat-induced unfolding. The primary focus is on the thermodynamics of reversible folding/unfolding transitions and the experimental methods that are available for extracting thermodynamic parameters. We highlight the importance of modelling both how the folding equilibrium depends on a perturbing variable such as temperature or denaturant concentration, and the importance of modelling the baselines in the experimental observables.

Global analysis of protein stability by temperature and chemical denaturation

The stability of a protein is a fundamental property that determines under which conditions, the protein is functional. Equilibrium unfolding with denaturants requires preparation of several samples and only provides the free energy of folding when performed at a single temperature. The typical sample requirement is around 0.5-1 mg of protein. If the stability of many proteins or protein variants needs to be determined, substantial protein production may be needed. Here we have determined the stability of acyl-coenzyme A binding protein at pH 5.3 and chymotrypsin inhibitor 2 at pH 3 and pH 6.25 by combined temperature and denaturant unfolding. We used a setup where tryptophan fluorescence is measured in quartz capillaries where only 10 μl is needed. Temperature unfolding of a series of 15 samples at increasing denaturant concentrations provided accurate and precise thermodynamic parameters. We find that the number of samples may be further reduced and less than 10 μg of protein in total are needed for reliable stability measurements. For assessment of stability of protein purified in small scale e.g. in micro plate format, our method will be highly applicable. The routine for fitting the experimental data is made available as a python notebook.

Dissecting the statistical properties of the linear extrapolation method of determining protein stability

The linear extrapolation method to determine protein stability from denaturant-induced unfolding experiments is based on the observation that the free energy of unfolding is often a linear function of the denaturant concentration. The value in the absence of denaturant is then estimated by extrapolation from this linear relationship. Parameters and their confidence intervals are typically estimated by nonlinear least-squares regression. We have compared different methods for calculating confidence intervals and found that a simple method based on linear theory gives accurate results. We have also compared three different parameterizations of the linear extrapolation method and show that the most commonly used form is problematic since the stability and m-value are correlated in the nonlinear least-squares analysis. Parameter correlation can in some cases causes problems in the estimation of confidence intervals and regions and should be avoided when possible. Two alternative parameterizations show much less correlation between parameters.

Consensus sequence design as a general strategy to create hyperstable, biologically active proteins

Consensus sequence design offers a promising strategy for designing proteins of high stability while retaining biological activity since it draws upon an evolutionary history in which residues important for both stability and function are likely to be conserved. Although there have been several reports of successful consensus design of individual targets, it is unclear from these anecdotal studies how often this approach succeeds and how often it fails. Here, we attempt to assess generality by designing consensus sequences for a set of six protein families with a range of chain lengths, structures, and activities. We characterize the resulting consensus proteins for stability, structure, and biological activities in an unbiased way. We find that all six consensus proteins adopt cooperatively folded structures in solution. Strikingly, four of six of these consensus proteins show increased thermodynamic stability over naturally occurring homologs. Each consensus protein tested for function maintained at least partial biological activity. Although peptide binding affinity by a consensus-designed SH3 is rather low, K _m values for consensus enzymes are similar to values from extant homologs. Although consensus enzymes are slower than extant homologs at low temperature, they are faster than some thermophilic enzymes at high temperature. An analysis of sequence properties shows consensus proteins to be enriched in charged residues, and rarified in uncharged polar residues. Sequence differences between consensus and extant homologs are predominantly located at weakly conserved surface residues, highlighting the importance of these residues in the success of the consensus strategy.

The phosphoglycerate kinase 1 variants found in carcinoma cells display different catalytic activity and conformational stability compared to the native enzyme

Cancer cells are able to survive in difficult conditions, reprogramming their metabolism according to their requirements. Under hypoxic conditions they shift from oxidative phosphorylation to aerobic glycolysis, a behavior known as Warburg effect. In the last years, glycolytic enzymes have been identified as potential targets for alternative anticancer therapies. Recently, phosphoglycerate kinase 1 (PGK1), an ubiquitous enzyme expressed in all somatic cells that catalyzes the seventh step of glycolysis which consists of the reversible phosphotransfer reaction from 1,3-bisphosphoglycerate to ADP, has been discovered to be overexpressed in many cancer types. Moreover, several somatic variants of PGK1 have been identified in tumors. In this study we analyzed the effect of the single nucleotide variants found in cancer tissues on the PGK1 structure and function. Our results clearly show that the variants display a decreased catalytic efficiency and/or thermodynamic stability and an altered local tertiary structure, as shown by the solved X-ray structures. The changes in the catalytic properties and in the stability of the PGK1 variants, mainly due to the local changes evidenced by the X-ray structures, suggest also changes in the functional role of PGK to support the biosynthetic need of the growing and proliferating tumour cells.

Insights on the conformational dynamics of human frataxin through modifications of loop-1

Human frataxin (FXN) is a highly conserved mitochondrial protein involved in iron homeostasis and activation of the iron-sulfur cluster assembly. FXN deficiency causes the neurodegenerative disease Friedreich's Ataxia. Here, we investigated the effect of alterations in loop-1, a stretch presumably essential for FXN function, on the conformational stability and dynamics of the native state. We generated four loop-1 variants, carrying substitutions, insertions and deletions. All of them were stable and well-folded proteins. Fast local motions (ps-ns) and slower long-range conformational dynamics (μs-ms) were altered in some mutants as judged by NMR. Particularly, loop-1 modifications impact on the dynamics of a distant region that includes residues from the β-sheet, helix α1 and the C-terminal. Remarkably, all the mutants retain the ability to activate cysteine desulfurase, even when two of them exhibit a strong decrease in iron binding, revealing a differential sensitivity of these functional features to loop-1 perturbation. Consequently, we found that even for a small and relatively rigid protein, engineering a loop segment enables to alter conformational dynamics through a long-range effect, preserving the native-state structure and important aspects of function.

Predicting the impact of Lynch syndrome-causing missense mutations from structural calculations

Accurate methods to assess the pathogenicity of mutations are needed to fully leverage the possibilities of genome sequencing in diagnosis. Current data-driven and bioinformatics approaches are, however, limited by the large number of new variations found in each newly sequenced genome, and often do not provide direct mechanistic insight. Here we demonstrate, for the first time, that saturation mutagenesis, biophysical modeling and co-variation analysis, performed in silico, can predict the abundance, metabolic stability, and function of proteins inside living cells. As a model system, we selected the human mismatch repair protein, MSH2, where missense variants are known to cause the hereditary cancer predisposition disease, known as Lynch syndrome. We show that the majority of disease-causing MSH2 mutations give rise to folding defects and proteasome-dependent degradation rather than inherent loss of function, and accordingly our in silico modeling data accurately identifies disease-causing mutations and outperforms the traditionally used genetic disease predictors. Thus, in conclusion, in silico biophysical modeling should be considered for making genotype-phenotype predictions and for diagnosis of Lynch syndrome, and perhaps other hereditary diseases.

The protein quality control system targets misfolded proteins for degradation. So far it has not been possible from sequence or structural data to predict the biological stability of a misfolded protein, or the effect of mutations on intracellular protein levels. Here we demonstrate that in silico saturation mutagenesis and biophysical calculations of the structural stability of the human mismatch repair protein MSH2 correlate with cellular protein levels, turnover and function. Of 24 different MSH2 variants, some of which are linked to Lynch syndrome, a destabilization of as little as 3 kcal/mol is sufficient to cause rapid degradation via the ubiquitin-proteasome pathway. Thus, biophysical modeling can, to a large extent, predict the metabolic stability of proteins. We also show that the same biophysical calculations can be used to distinguish with high accuracy neutral sequence variation from pathogenic variants, and that the calculations outperform several traditionally used disease predictors. We therefore suggest the method to be of potential value for patient stratification in Lynch syndrome, and perhaps other hereditary diseases.

Single-Nucleotide Polymorphism of PPARγ, a Protein at the Crossroads of Physiological and Pathological Processes

Genome polymorphisms are responsible for phenotypic differences between humans and for individual susceptibility to genetic diseases and therapeutic responses. Non-synonymous single-nucleotide polymorphisms (nsSNPs) lead to protein variants with a change in the amino acid sequence that may affect the structure and/or function of the protein and may be utilized as efficient structural and functional markers of association to complex diseases. This study is focused on nsSNP variants of the ligand binding domain of PPARγ a nuclear receptor in the superfamily of ligand inducible transcription factors that play an important role in regulating lipid metabolism and in several processes ranging from cellular differentiation and development to carcinogenesis. Here we selected nine nsSNPs variants of the PPARγ ligand binding domain, V290M, R357A, R397C, F360L, P467L, Q286P, R288H, E324K, and E460K, expressed in cancer tissues and/or associated with partial lipodystrophy and insulin resistance. The effects of a single amino acid change on the thermodynamic stability of PPARγ, its spectral properties, and molecular dynamics have been investigated. The nsSNPs PPARγ variants show alteration of dynamics and tertiary contacts that impair the correct reciprocal positioning of helices 3 and 12, crucially important for PPARγ functioning.

Principles for designing proteins with cavities formed by curved β sheets

Active sites and ligand-binding cavities in native proteins are often formed by curved β sheets, and the ability to control β-sheet curvature would allow design of binding proteins with cavities customized to specific ligands. Toward this end, we investigated the mechanisms controlling β-sheet curvature by studying the geometry of β sheets in naturally occurring protein structures and folding simulations. The principles emerging from this analysis were used to design, de novo, a series of proteins with curved β sheets topped with α helices. Nuclear magnetic resonance and crystal structures of the designs closely match the computational models, showing that β-sheet curvature can be controlled with atomic-level accuracy. Our approach enables the design of proteins with cavities and provides a route to custom design ligand-binding and catalytic sites.

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.

A Soluble, Folded Protein without Charged Amino Acid Residues

Charges are considered an integral part of protein structure and function, enhancing solubility and providing specificity in molecular interactions. We wished to investigate whether charged amino acids are indeed required for protein biogenesis and whether a protein completely free of titratable side chains can maintain solubility, stability, and function. As a model, we used a cellulose-binding domain from Cellulomonas fimi, which, among proteins of more than 100 amino acids, presently is the least charged in the Protein Data Bank, with a total of only four titratable residues. We find that the protein shows a surprising resilience toward extremes of pH, demonstrating stability and function (cellulose binding) in the pH range from 2 to 11. To ask whether the four charged residues present were required for these properties of this protein, we altered them to nontitratable ones. Remarkably, this chargeless protein is produced reasonably well in Escherichia coli, retains its stable three-dimensional structure, and is still capable of strong cellulose binding. To further deprive this protein of charges, we removed the N-terminal charge by acetylation and studied the protein at pH 2, where the C-terminus is effectively protonated. Under these conditions, the protein retains its function and proved to be both soluble and have a reversible folding-unfolding transition. To the best of our knowledge, this is the first time a soluble, functional protein with no titratable side chains has been produced.

Understanding the effect of alternative splicing in the folding and function of the second PDZ from protein tyrosine phosphatase-BL

PDZ domains are the most prominent biological structural domains involved in protein-protein interactions in the human cell. The second PDZ domain of the protein tyrosine phosphatase BL (PDZ2) interacts and binds the C-termini of the tumour suppressor protein APC and of the LIM domain-containing protein RIL. One isoform of PDZ2 (PDZ2as) involves an alternative spliced form that exhibits an insertion of 5 residues in a loop. PDZ2as abrogates binding to its partners, even if the insertion is directly located in its binding pocket. Here, we investigate the folding and function of PDZ2as, in comparison to the previously characterized PDZ2 domain. Data reveal that, whilst the thermodynamic stability of PDZ2as appears as nearly identical to that of PDZ2, the insertion of 5 amino acids induces formation of some weak transient non-native interactions in the folding transition state, as mirrored by a concomitant increase of both the folding and unfolding rate constants. From a functional perspective, we show that the decrease in affinity is caused by a pronounced decrease of the association rate constants (by nearly ten fold), with no effect on the microscopic dissociation rate constants. The results are briefly discussed in the context of previous work on PDZ domains.

Interplay between partner and ligand facilitates the folding and binding of an intrinsically disordered protein

Protein-protein interactions are at the heart of regulatory and signaling processes in the cell. In many interactions, one or both proteins are disordered before association. However, this disorder in the unbound state does not prevent many of these proteins folding to a well-defined, ordered structure in the bound state. Here we examine a typical system, where a small disordered protein (PUMA, p53 upregulated modulator of apoptosis) folds to an α-helix when bound to a groove on the surface of a folded protein (MCL-1, induced myeloid leukemia cell differentiation protein). We follow the association of these proteins using rapid-mixing stopped flow, and examine how the kinetic behavior is perturbed by denaturant and carefully chosen mutations. We demonstrate the utility of methods developed for the study of monomeric protein folding, including β-Tanford values, Leffler α, Φ-value analysis, and coarse-grained simulations, and propose a self-consistent mechanism for binding. Folding of the disordered protein before binding does not appear to be required and few, if any, specific interactions are required to commit to association. The majority of PUMA folding occurs after the transition state, in the presence of MCL-1. We also examine the role of the side chains of folded MCL-1 that make up the binding groove and find that many favor equilibrium binding but, surprisingly, inhibit the association process.

High thermodynamic stability of parametrically designed helical bundles

We describe a procedure for designing proteins with backbones produced by varying the parameters in the Crick coiled coil-generating equations. Combinatorial design calculations identify low-energy sequences for alternative helix supercoil arrangements, and the helices in the lowest-energy arrangements are connected by loop building. We design an antiparallel monomeric untwisted three-helix bundle with 80-residue helices, an antiparallel monomeric right-handed four-helix bundle, and a pentameric parallel left-handed five-helix bundle. The designed proteins are extremely stable (extrapolated ΔGfold > 60 kilocalories per mole), and their crystal structures are close to those of the design models with nearly identical core packing between the helices. The approach enables the custom design of hyperstable proteins with fine-tuned geometries for a wide range of applications.

Evolution and design of protein structure by folding nucleus symmetric expansion

Models of symmetric protein evolution typically invoke gene duplication and fusion events, in which repetition of a structural motif generates foldable, stable symmetric protein architecture. Success of such evolutionary processes suggests that the duplicated structural motif must be capable of nucleating protein folding. If correct, symmetric expansion of a folding nucleus sequence derived from an extant symmetric fold may be an elegant and computationally tractable solution to de novo protein design. We report the efficient de novo design of a β-trefoil protein by symmetric expansion of a β-trefoil folding nucleus, previously identified by ɸ-value analysis. The resulting protein, having exact sequence symmetry, exhibits superior folding properties compared to its naturally evolved progenitor-with the potential for redundant folding nuclei. In principle, folding nucleus symmetric expansion can be applied to any given symmetric protein fold (that is, nearly one-third of the known proteome) provided information of the folding nucleus is available.

Decreased amyloidogenicity caused by mutational modulation of surface properties of the immunoglobulin light chain BRE variable domain

Amyloid formation by immunoglobulin light chain (LC) proteins is associated with amyloid light chain (AL) amyloidosis. Destabilization of the native state of the variable domain of the LC (VL) is known to be one of the critical factors in promoting the formation of amyloid fibrils. However, determining the key residues involved in this destabilization remains challenging, because of the existence of a number of intrinsic sequence variations within VL. In this study, we identified the key residues for destabilization of the native state of amyloidogenic VL in the LC of BRE by analyzing the stability of chimeric mutants of BRE and REI VL; the latter immunoglobulin is not associated with AL amyloidosis. The results suggest that the surface-exposed residues N45 and D50 are the key residues in the destabilization of the native state of BRE VL. Point mutations at the corresponding residues in REI VL (K45N, E50D, and K45N/E50D) destabilized the native state and increased amyloidogenicity. However, the reverse mutations in BRE VL (N45K, D50E, and N45K/D50E) re-established the native state and decreased amyloidogenicity. Thus, analyses using chimeras and point mutants successfully elucidated the key residues involved in BRE VL destabilization and increased amyloidogenic propensity. These results also suggest that the modulation of surface properties of wild-type VL may improve their stability and prevent the formation of amyloid fibrils.

Random single amino acid deletion sampling unveils structural tolerance and the benefits of helical registry shift on GFP folding and structure

Altering a protein’s backbone through amino acid deletion is a common evolutionary mutational mechanism, but is generally ignored during protein engineering primarily because its effect on the folding-structure-function relationship is difficult to predict. Using directed evolution, enhanced green fluorescent protein (EGFP) was observed to tolerate residue deletion across the breadth of the protein, particularly within short and long loops, helical elements, and at the termini of strands. A variant with G4 removed from a helix (EGFP^G4Δ) conferred significantly higher cellular fluorescence. Folding analysis revealed that EGFP^G4Δ retained more structure upon unfolding and refolded with almost 100% efficiency but at the expense of thermodynamic stability. The EGFP^G4Δ structure revealed that G4 deletion caused a beneficial helical registry shift resulting in a new polar interaction network, which potentially stabilizes a cis proline peptide bond and links secondary structure elements. Thus, deletion mutations and registry shifts can enhance proteins through structural rearrangements not possible by substitution mutations alone.

•

Using directed evolution, the impact of amino acid deletion on EGFP is explored

•

Loops, helices, and strand termini are especially tolerant to amino acid deletion

•

A deletion mutant that enhances cellular production and fluorescence is identified

•

Structure reveals that a helical registry shift creates a new polar network

Structural basis of α-catenin recognition by EspB from enterohaemorrhagic E. coli based on hybrid strategy using low-resolution structural and protein dissection

Enterohaemorrhagic E. coli (EHEC) induces actin reorganization of host cells by injecting various effectors into host cytosol through type III secretion systems. EspB is the natively partially folded EHEC effector which binds to host α-catenin to promote the actin bundling. However, its structural basis is poorly understood. Here, we characterize the overall structural properties of EspB based on low-resolution structural data in conjunction with protein dissection strategy. EspB showed a unique thermal response involving cold denaturation in the presence of denaturant according to far-UV circular dichroism (CD). Small angle X-ray scattering revealed the formation of a highly extended structure of EspB comparable to the ideal random coil. Various disorder predictions as well as CD spectra of EspB fragments identified the presence of α-helical structures around G41 to Q70. The fragment corresponding to this region indicated the thermal response similar to EspB. Moreover, this fragment showed a high affinity to C-terminal vinculin homology domain of α-catenin. The results clarified the importance of preformed α-helix of EspB for recognition of α-catenin.

Effect of single amino acid substitution observed in cancer on Pim-1 kinase thermodynamic stability and structure

Pim-1 kinase, a serine/threonine protein kinase encoded by the pim proto-oncogene, is involved in several signalling pathways such as the regulation of cell cycle progression and apoptosis. Many cancer types show high expression levels of Pim kinases and particularly Pim-1 has been linked to the initiation and progression of the malignant phenotype. In several cancer tissues somatic Pim-1 mutants have been identified. These natural variants are nonsynonymous single nucleotide polymorphisms, variations of a single nucleotide occurring in the coding region and leading to amino acid substitutions. In this study we investigated the effect of amino acid substitution on the structural stability and on the activity of Pim-1 kinase. We expressed and purified some of the mutants of Pim-1 kinase that are expressed in cancer tissues and reported in the single nucleotide polymorphisms database. The point mutations in the variants significantly affect the conformation of the native state of Pim-1. All the mutants, expressed as soluble recombinant proteins, show a decreased thermal and thermodynamic stability and a lower activation energy values for kinase activity. The decreased stability accompanied by an increased flexibility suggests that Pim-1 variants may be involved in a wider network of protein interactions. All mutants bound ATP and ATP mimetic inhibitors with comparable IC50 values suggesting that the studied Pim-1 kinase mutants can be efficiently targeted with inhibitors developed for the wild type protein.

Slow, reversible, coupled folding and binding of the spectrin tetramerization domain

Many intrinsically disordered proteins (IDPs) are significantly unstructured under physiological conditions. A number of these IDPs have been shown to undergo coupled folding and binding reactions whereby they can gain structure upon association with an appropriate partner protein. In general, these systems display weaker binding affinities than do systems with association between completely structured domains, with micromolar K(d) values appearing typical. One such system is the association between α- and β-spectrin, where two partially structured, incomplete domains associate to form a fully structured, three-helix bundle, the spectrin tetramerization domain. Here, we use this model system to demonstrate a method for fitting association and dissociation kinetic traces where, using typical biophysical concentrations, the association reactions are expected to be highly reversible. We elucidate the unusually slow, two-state kinetics of spectrin assembly in solution. The advantages of studying kinetics in this regime include the potential for gaining equilibrium constants as well as rate constants, and for performing experiments with low protein concentrations. We suggest that this approach would be particularly appropriate for high-throughput mutational analysis of two-state reversible binding processes.

Structural stability of human protein tyrosine phosphatase ρ catalytic domain: effect of point mutations

Protein tyrosine phosphatase ρ (PTPρ) belongs to the classical receptor type IIB family of protein tyrosine phosphatase, the most frequently mutated tyrosine phosphatase in human cancer. There are evidences to suggest that PTPρ may act as a tumor suppressor gene and dysregulation of Tyr phosphorylation can be observed in diverse diseases, such as diabetes, immune deficiencies and cancer. PTPρ variants in the catalytic domain have been identified in cancer tissues. These natural variants are nonsynonymous single nucleotide polymorphisms, variations of a single nucleotide occurring in the coding region and leading to amino acid substitutions. In this study we investigated the effect of amino acid substitution on the structural stability and on the activity of the membrane-proximal catalytic domain of PTPρ. We expressed and purified as soluble recombinant proteins some of the mutants of the membrane-proximal catalytic domain of PTPρ identified in colorectal cancer and in the single nucleotide polymorphisms database. The mutants show a decreased thermal and thermodynamic stability and decreased activation energy relative to phosphatase activity, when compared to wild- type. All the variants show three-state equilibrium unfolding transitions similar to that of the wild- type, with the accumulation of a folding intermediate populated at ~4.0 M urea.

Secretion from bacterial versus mammalian cells yields a recombinant scFv with variable folding properties

Escherichia coli (E. coli) is the most commonly used organism for expressing antibody fragments such as single chain antibody Fvs (scFvs). Previously, we have utilized E. coli to express well-folded scFvs for characterization and engineering purposes with the goal of using these engineered proteins as building blocks for generating IgG-like bispecific antibodies (BsAbs). In the study, described here, we observed a significant difference in the secondary structure of an scFv produced in E. coli and the same scFv expressed and secreted from chinese hamster ovary (CHO) cells as part of a BsAb. We devised a proteolytic procedure to separate the CHO-derived scFv from its antibody-fusion partner and compared its properties with those of the E. coli-derived scFv. In comparison to the CHO-derived scFv, the E. coli-derived scFv was found trapped in a misfolded, but monomeric state that was stable for months at 4 °C. The misfolded state bound antigen in a heterogeneous fashion that included non-specific binding, which made functional characterization challenging. This odd incidence of obtaining a misfolded scFv from bacteria suggests careful characterization of the folded properties of bacterially expressed scFvs is warranted if anomalous issues with antigen-binding or non-specificity occur during an engineering campaign. Additionally, our proteolytic methodology for obtaining significant levels of intact scFvs from highly expressed IgG-like antibody proteins serves as a robust method for producing scFvs in CHO without the use of designed cleavage motifs.

Probing protein surface with a solvent mimetic carbene coupled to detection by mass spectrometry

Much knowledge into protein folding, ligand binding, and complex formation can be derived from the examination of the nature and size of the accessible surface area (SASA) of the polypeptide chain, a key parameter in protein science not directly measurable in an experimental fashion. To this end, an ideal chemical approach should aim at exerting solvent mimicry and achieving minimal selectivity to probe the protein surface regardless of its chemical nature. The choice of the photoreagent diazirine to fulfill these goals arises from its size comparable to water and from being a convenient source of the extremely reactive methylene carbene (:CH(2)). The ensuing methylation depends primarily on the solvent accessibility of the polypeptide chain, turning it into a valuable signal to address experimentally the measurement of SASA in proteins. The superb sensitivity and high resolution of modern mass spectrometry techniques allows us to derive a quantitative signal proportional to the extent of modification (EM) of the sample. Thus, diazirine labeling coupled to electrospray mass spectrometry (ESI-MS) detection can shed light on conformational features of the native as well as non-native states, not easily addressable by other methods. Enzymatic fragmentation of the polypeptide chain at the level of small peptides allows us to locate the covalent tag along the amino acid sequence, therefore enabling the construction of a map of solvent accessibility. Moreover, by subsequent MS/MS analysis of peptides, we demonstrate here the feasibility of attaining amino acid resolution in defining the target sites.

In silico Evaluation of Crosslinking Effects on Denaturant m(eq) values and ΔCp upon Protein Unfolding

Important thermodynamic parameters including denaturant equilibrium m values (m(eq)) and heat capacity changes (ΔCp) can be predicted based on changes in Solvent Accessible Surface Area (SASA) upon unfolding. Crosslinks such as disulfide bonds influence the stability of the proteins by decreasing the entropy gain as well as reduction of SASA of unfolded state. The aim of the study was to develop mathematical models to predict the effect of crosslinks on ΔSASA and ultimately on m(eq) and ΔCp based on in silico methods. Changes of SASA upon computationally simulated unfolding were calculated for a set of 45 proteins with known m(eq) and ΔCp values and the effect of crosslinks on ΔSASA of unfolding was investigated. The results were used to predict the m(eq) of denaturation for guanidine hydrochloride and urea, as well as ΔCp for the studied proteins with overall error of 20%, 31% and 17%, respectively. The results of the current study were in close agreement with those obtained from the previous studies.

Fine details of IGF-1R activation, inhibition, and asymmetry determined by associated hydrogen /deuterium-exchange and peptide mass mapping

The structural features of the asymmetric activated states of the insulin receptor family are still poorly understood. We investigated hydrogen/deuterium (H/D)-exchange within the extracellular domain of the type-I insulin-like growth factor receptor (IGF-1R) in the absence and presence of IGF-1 (active state) and in the presence of antibody inhibitors (inactive state). Near complete coverage of the 210 kDa receptor sequence was obtained by mass mapping of proteolytically derived peptides at all H/D-exchange time points. The data provide details regarding solvent exposure and dynamics across the extracellular region as well as conformational changes induced by activation or inactivation. Multiple peptides, distant in structure, individually demonstrated two distinct H/D-exchange rates, suggesting that each of these peptides exists in two separate environments in IGF-1R. The dual-exchange behavior of these peptides was enhanced on ligand binding and eliminated on inhibitor binding, clearly associating these regions with active state asymmetry and enabling them to serve as reporters of receptor activity.

Calorimetric study of a series of designed repeat proteins: modular structure and modular folding

Repeat proteins comprise tandem arrays of a small structural motif. Their structure is defined and stabilized by interactions between residues that are close in the primary sequence. Several studies have investigated whether their structural modularity translates into modular thermodynamic properties. Tetratricopeptide repeat proteins (TPRs) are a class in which the repeated unit is a 34 amino acid helix-turn-helix motif. In this work, we use differential scanning calorimetry (DSC) to study the equilibrium stability of a series of TPR proteins with different numbers of an identical consensus repeat, from 2 to 20, CTPRa2 to CTPRa20. The DSC data provides direct evidence that the folding/unfolding transition of CTPR proteins does not fit a two-state folding model. Our results confirm and expand earlier studies on TPR proteins, which showed that apparent two-state unfolding curves are better fit by linear statistical mechanics models: 1D Ising models in which each repeat is treated as an independent folding unit.

The denatured state dictates the topology of two proteins with almost identical sequence but different native structure and function

The protein folding problem is often studied by comparing the mechanisms of proteins sharing the same structure but different sequence. The recent design of the two proteins G(A)88 and G(B)88, displaying different structures and functions while sharing 88% sequence identity (49 out of 56 amino acids), allows the unique opportunity for a complementary approach. At which stage of its folding pathway does a protein commit to a given topology? Which residues are crucial in directing folding mechanisms to a given structure? By using a combination of biophysical and computational techniques, we have characterized the folding of both G(A)88 and G(B)88. We show that, contrary to expectation, G(B)88, characterized by a native α+β fold, displays in the denatured state a content of native-like helical structure greater than G(A)88, which is all-α in its native state. Both experiments and simulations indicate that such residual structure may be tuned by changing pH. Thus, despite the high sequence identity, the folding pathways for these two proteins appear to diverge as early as in the denatured state. Our results suggest a mechanism whereby protein topology is committed very early along the folding pathway, being imprinted in the residual structure of the denatured state.

Decreased stability and increased formation of soluble aggregates by immature superoxide dismutase do not account for disease severity in ALS

Protein aggregation is a hallmark of many diseases, including amyotrophic lateral sclerosis (ALS), where aggregation of Cu/Zn superoxide dismutase (SOD1) is implicated in causing neurodegeneration. Recent studies have suggested that destabilization and aggregation of the most immature form of SOD1, the disulfide-reduced, unmetallated (apo) protein is particularly important in causing ALS. We report herein in depth analyses of the effects of chemically and structurally diverse ALS-associated mutations on the stability and aggregation of reduced apo SOD1. In contrast with previous studies, we find that various reduced apo SOD1 mutants undergo highly reversible thermal denaturation with little aggregation, enabling quantitative thermodynamic stability analyses. In the absence of ALS-associated mutations, reduced apo SOD1 is marginally stable but predominantly folded. Mutations generally result in slight decreases to substantial increases in the fraction of unfolded protein. Calorimetry, ultracentrifugation, and light scattering show that all mutations enhance aggregation propensity, with the effects varying widely, from subtle increases in most cases, to pronounced formation of 40-100 nm soluble aggregates by A4V, a mutation that is associated with particularly short disease duration. Interestingly, although there is a correlation between observed aggregation and stability, there is minimal to no correlation between observed aggregation, predicted aggregation propensity, and disease characteristics. These findings suggest that reduced apo SOD1 does not play a dominant role in modulating disease. Rather, additional and/or multiple forms of SOD1 and additional biophysical and biological factors are needed to account for the toxicity of mutant SOD1 in ALS.

Stability curves of laboratory evolved thermostable mutants of a Bacillus subtilis lipase

Shape of the protein stability curves changes to achieve higher melting temperature. Broadly, these changes have been classified as upward shift (increased G(s)), rightward shift (increase in T(s)) and flattening of the stability curves (decrease in C(p)). Comparative studies on homologous mesophilic-thermophilic protein pairs highlighted the differential contribution of these three strategies amongst proteins. But unambiguous way of identification of the strategies, which will be preferred for a protein, is still not achieved. We have performed comparative thermodynamic studies using differential scanning calorimeter (DSC) on thermostable variants of a lipase from Bacillus subtilis. These variants are products of 1, 2, 3 and 4 rounds of directed evolution and harbor mutations having definite contribution in thermostability unlike natural thermophilic proteins. We have shown that upward and rightward shift in stability curves are prime strategies in this lipase. Our results along with that from the other study on laboratory evolved xylanase A suggest that optimization of suboptimal thermodynamic parameters is having a dominant influence in selection of thermodynamic strategies for higher thermostability.

Folding and association of thermophilic dimeric and trimeric DsrEFH proteins: Tm0979 and Mth1491

Although the majority of natural proteins exist as protein-protein complexes, the molecular basis for the formation and regulation of such interactions and the evolution of protein interfaces remain poorly understood. We have investigated these phenomena by characterizing the thermal and chemical denaturation of thermophilic DsrEFH proteins that have a common subunit fold but distinct quaternary structures: homodimeric Tm0979 and homotrimeric Mth1491. Tm0979 forms a moderate affinity dimer, and a monomeric intermediate is readily populated at equilibrium and during folding kinetics. In contrast, the Mth1491 trimer has extremely high stability, so that a monomeric form is not measurably populated at equilibrium, although it may be during folding kinetics. A common mechanism for evolution of quaternary structures may be facile formation of a relatively stable monomeric species, with stabilizing intermolecular interactions centering on alternative environments for a beta-strand at the edge of the monomer, augmented by malleable hydrophobic interactions. The exceptional trimer stability arises from a remarkably slow unfolding rate constant, 6.5 x 10(-13) s(-1), which is a common characteristic of highly stable thermophilic and/or oligomeric proteins. The folding characteristics of Tm0979 and Mth1491 have interesting implications for assembly and regulation of homo- and heterooligomeric proteins in vivo.

Molecular origin of constant m-values, denatured state collapse, and residue-dependent transition midpoints in globular proteins

Experiments show that for many two-state folders the free energy of the native state, DeltaG(ND)([C]), changes linearly as the denaturant concentration, [C], is varied. The slope {m = [dDeltaG(ND)([C])]/(d[C])}, is nearly constant. According to the transfer model, the m-value is associated with the difference in the surface area between the native (N) and denatured (D) state, which should be a function of DeltaR(g)(2), the difference in the square of the radius of gyration between the D and N states. Single-molecule experiments show that the R(g) of the structurally heterogeneous denatured state undergoes an equilibrium collapse transition as [C] decreases, which implies m also should be [C]-dependent. We resolve the conundrum between constant m-values and [C]-dependent changes in R(g) using molecular simulations of a coarse-grained representation of protein L, and the molecular transfer model, for which the equilibrium folding can be accurately calculated as a function of denaturant (urea) concentration. In agreement with experiment, we find that over a large range of denaturant concentration (>3 M) the m-value is a constant, whereas under strongly renaturing conditions (<3 M), it depends on [C]. The m-value is a constant above [C] > 3 M because the [C]-dependent changes in the surface area of the backbone groups, which make the largest contribution to m, are relatively narrow in the denatured state. The burial of the backbone and hydrophobic side chains gives rise to substantial surface area changes below [C] < 3 M, leading to collapse in the denatured state of protein L. Dissection of the contribution of various amino acids to the total surface area change with [C] shows that both the sequence context and residual structure are important. There are [C]-dependent variations in the surface area for chemically identical groups such as the backbone or Ala. Consequently, the midpoints of transition of individual residues vary significantly (which we call the Holtzer effect) even though global folding can be described as an all-or-none transition. The collapse is specific in nature, resulting in the formation of compact structures with appreciable populations of nativelike secondary structural elements. The collapse transition is driven by the loss of favorable residue-solvent interactions and a concomitant increase in the strength of intrapeptide interactions with a decreasing [C]. The strength of these interactions is nonuniformly distributed throughout the structure of protein L. Certain secondary structure elements have stronger [C]-dependent interactions than others in the denatured state.

ProtSA: a web application for calculating sequence specific protein solvent accessibilities in the unfolded ensemble

Background

The stability of proteins is governed by the heat capacity, enthalpy and entropy changes of folding, which are strongly correlated to the change in solvent accessible surface area experienced by the polypeptide. While the surface exposed in the folded state can be easily determined, accessibilities for the unfolded state at the atomic level cannot be obtained experimentally and are typically estimated using simplistic models of the unfolded ensemble. A web application providing realistic accessibilities of the unfolded ensemble of a given protein at the atomic level will prove useful.

Results

ProtSA, a web application that calculates sequence-specific solvent accessibilities of the unfolded state ensembles of proteins has been developed and made freely available to the scientific community. The input is the amino acid sequence of the protein of interest. ProtSA follows a previously published calculation protocol which uses the Flexible-Meccano algorithm to generate unfolded conformations representative of the unfolded ensemble of the protein, and uses the exact analytical software ALPHASURF to calculate atom solvent accessibilities, which are averaged on the ensemble.

Conclusion

ProtSA is a novel tool for the researcher investigating protein folding energetics. The sequence specific atom accessibilities provided by ProtSA will allow obtaining better estimates of the contribution of the hydrophobic effect to the free energy of folding, will help to refine existing parameterizations of protein folding energetics, and will be useful to understand the influence of point mutations on protein stability.

Folding mechanism of the C-terminal domain of nucleophosmin: residual structure in the denatured state and its pathophysiological significance

Nucleophosmin (NPM1) is a ubiquitously expressed protein and is one of the most abundant proteins found in the nucleolus. Naturally occurring mutations in the C-terminal domain of nucleophosmin (Cter-NPM1) are found in approximately 30% of patients with acute myeloid leukemia (AML). These mutations cause changes at the C terminus of NPM1 that lead to denaturation of the protein, a critical factor in determining aberrant translocation of NPM1 to the cytosol. Hence, this protein system represents an ideal candidate to investigate the relations between folding and unfolding and disease. Here we report the characterization of the folding and unfolding kinetics of Cter-NPM1. Data reveal that this small helical domain folds via a compact denatured state, displaying a malleable residual structure. Moreover, analysis of folding rate constants measured under different experimental conditions suggests that the existence of a preorganized structure in the denatured state accelerates folding, implying a native-like residual structure. Because a major feature of Cter-NPM1 mutants responsible for AML is a reduction in stability of the protein and thus prevalence of a denatured state even under physiological conditions, our findings may pave the way to further studies with the aim of designing chemicals capable of interacting with the "pathological" mutants to stabilize the native conformation.

Studying the unfolding kinetics of proteins under pressure using long molecular dynamic simulation runs

The usefulness of computational methods such as molecular dynamics simulation has been extensively established for studying systems in equilibrium. Nevertheless, its application to complex non-equilibrium biological processes such as protein unfolding has been generally regarded as producing results which cannot be interpreted straightforwardly. In the present study, we present results for the kinetics of unfolding of apomyoglobin, based on the analysis of long simulation runs of this protein in solution at 3 kbar (1 atm = 1.01325, bar = 101,325 Pa). We hereby demonstrate that the analysis of the data collected within a simulated time span of 0.18 mus suffices for producing results, which coincide remarkably with the available unfolding kinetics experimental data. This not only validates molecular dynamics simulation as a valuable alternative for studying non-equilibrium processes, but also enables a detailed analysis of the actual structural mechanism which underlies the unfolding process of proteins under elusive denaturing conditions such as high pressure.

Thermodynamic analysis of denaturant-induced unfolding of HodC69S protein supports a three-state mechanism

Thermodynamic stability parameters and the equilibrium unfolding mechanism of His 6HodC69S, a mutant of 1 H-3-hydroxy-4-oxoquinaldine 2,4-dioxygenase (Hod) having a Cys to Ser exchange at position 69 and an N-terminal hexahistidine tag (His 6HodC69S), have been derived from isothermal unfolding studies using guanidine hydrochloride (GdnHCl) or urea as denaturants. The conformational changes were monitored by following changes in circular dichroism (CD), fluorescence, and dynamic light scattering (DLS), and the resulting transition curves were analyzed on the basis of a sequential three-state model N = I = D. The structural changes have been correlated to catalytic activity, and the contribution to stability of the disulfide bond between residues C37 and C184 in the native protein has been established. A prominent result of the present study is the finding that, independent of the method used for denaturing the protein, the unfolding mechanism always comprises three states which can be characterized by, within error limits, identical sets of thermodynamic parameters. Apparent deviations from three-state unfolding can be rationalized by the inability of a spectroscopic probe to discriminate clearly between native, intermediate, and unfolded ensembles. This was the case for the CD-monitored urea unfolding curve.