Computational analysis of folding and mutation properties of C5 domain of myosin binding protein C

Thermal folding molecular dynamics simulations of the domain C5 of Myosin binding protein C were performed using a native-centric model to study the role of three mutations related to Familial Hypertrophic Cardiomyopathy. Mutation of Asn755 causes the largest shift of the folding temperature, and the residue is located in the CFGA' beta-sheet featuring the highest phi-values. The mutation thus appears to reduce the thermodynamic stability in agreement with experimental data. The mutations on Arg654 and Arg668, conversely, cause little change in the folding temperature and they reside in the low phi-value BDE beta-sheet, so that their pathological role cannot be related to impairment of the folding process but possibly to the binding with target molecules. As the typical signature of Domain C5 is the presence of a longer and destibilizing CD-loop with respect to the other Ig-like domains, we completed the work with a bioinformatic analysis of this loop showing a high density of negative charge and low hydrophobicity. This indicates the CD-loop as a natively unfolded sequence with a likely coupling between folding and ligand binding.     

Analyzing pathogenic mutations of C5 domain from cardiac myosin binding protein C through MD simulations

The folding properties of wild type and mutants of domain C5 from cardiac myosin binding protein C have been investigated via molecular dynamics simulations within the framework of a native-centric and coarse-grained model. The relevance of a mutation has been assessed through the shift in the unfolding temperature, the change in the unfolding rate it determines and Phi-values analysis. In a previous paper (Guardiani et al. Biophys J 94:1403-1411, 2008), we performed Kinetic simulations on native contact formation revealing an entropy-driven folding pathway originating near the FG and DE loops. This folding mechanism allowed also a possible interpretation of the molecular impact of the three mutations, Arg14His, Arg28His and Asn115Lys involved in the Familial Hypertrophic Cardiomyopathy. Here we extend that analysis by enriching the mutant pool and we identify a correlation between unfolding rates and the number of native contacts retained in the transition state.     

In silico Investigation of the Disease-Associated Retinoschisin C110Y and C219G Mutants

The juvenile X-linked retinoschisis (XLRS) is a retinal disease caused by mutations in the secretory protein, retinoschisin (RS1). Majority of the disease is resulted from single point mutations on the RS1 discoidin domain with cysteine mutations being related to some of the more severe cases of XLRS. Previous studies have indicated that two mutations (C110Y and C219G), which involve cysteines that form intramolecular disulfide bonds in the native discoidin domain, resulted in different oligomerization states of the proteins and did not correlate with the degree of protein stability as calculated by the change in folding free energy. Through homology modeling, bioinformatics predictions, molecular dynamics (MD) and docking simulations, we attempt to investigate the effects of these two mutations on the structure of the RS1 discoidin domain in relevance to the discrepancy found between structural stability and aggregation propensity. Based on our findings, this discrepancy can be explained by the ability of C110Y mutant to establish suitable modules for initiating amorphous aggregation and to expand the aggregating mass through predominantly hydrophobic interactions. The low capability of C219G mutant to oligomerize, on the other hand, may be due to its greater structural instability and lesser hydrophobic tendency, two properties that may be unsupportive of aggregation. The results, altogether, indicate that aggregation propensity in the RS1 C110Y mutant is dependent upon the formation of suitable aggregating substrates for propagation of aggregation and not directly related to or determined by overall structural instability. As for the wildtype protein, the binding specificity of the spikes for biological function and the formation of octameric structure are contributed by important loop interactions, as well as evolved structural and sequence-based properties that prevent aggregation.     

The Saccharomyces cerevisiae YFR041C/ERJ5 gene encoding a type I membrane protein with a J domain is required to preserve the folding capacity of the endoplasmic reticulum

YFR041C/ERJ5 was identified in Saccharomyces cerevisiae as a gene regulated by the unfolded protein response pathway (UPR). The open reading frame of the gene has a J domain characteristic of the DnaJ chaperone family of proteins that regulate the activity of Hsp70 chaperones. We determined the expression and topology of Erj5p, a type I membrane protein with a J domain in the lumen of the endoplasmic reticulum (ER) that colocalizes with Kar2p, the major Hsp70 in the yeast ER. We identified synthetic interactions of Deltaerj5 with mutations in genes involved in protein folding in the ER (kar2-159, Deltascj1Deltajem1) and in the induction of the unfolded protein response (Deltaire1). Loss of Erj5p in yeast cells with impaired ER protein folding capacity increased sensitivity to agents that cause ER stress. We identified the ERJ5 mRNA and confirmed that agents that promote accumulation of misfolded proteins in the ER regulate its abundance. We found that loss of the non-essential ERJ5 gene leads to a constitutively induced UPR, indicating that ERJ5 is required for maintenance of an optimal folding environment in the yeast ER.     

NMR structure of the LCCL domain and implications for DFNA9 deafness disorder

The LCCL domain is a recently discovered, conserved protein module named after its presence in Limulus factor C, cochlear protein Coch-5b2 and late gestation lung protein Lgl1. The LCCL domain plays a key role in the autosomal dominant human deafness disorder DFNA9. Here we report the nuclear magnetic resonance (NMR) structure of the LCCL domain from human Coch-5b2, where dominant mutations leading to DFNA9 deafness disorder have been identified. The fold is novel. Four of the five known DFNA9 mutations are shown to involve at least partially solvent-exposed residues. Except for the Trp91Arg mutant, expression of these four LCCL mutants resulted in misfolded proteins. These results suggest that Trp91 participates in the interaction with a binding partner. The unexpected sensitivity of the fold with respect to mutations of solvent-accessible residues might be attributed to interference with the folding pathway of this disulfide-containing domain.     

Mutations of the anti-mullerian hormone gene in patients with persistent mullerian duct syndrome: biosynthesis, secretion, and processing of the abnormal proteins and analysis using a three-dimensional model

Anti-Müllerian hormone (AMH), a TGF-beta family member, determines whether an individual develops a uterus and Fallopian tubes. Mutations in the AMH gene lead to persistent Müllerian duct syndrome in males. The wild-type human AMH protein is synthesized as a disulfide-linked dimer of two identical 70-kDa polypeptides, which undergoes proteolytic processing to generate a 110-kDa N-terminal dimer and a bioactive 25-kDa TGF-beta-like C-terminal dimer. We have studied the biosynthesis and secretion of wild-type AMH and of seven persistent Müllerian duct syndrome proteins, containing mutations in either the N- or C-terminal domain. Mutant proteins lacking the C-terminal domain are secreted more rapidly than full-length AMH, whereas single amino acid changes in both domains can have profound effects on protein stability and folding. The addition of a cysteine in an N-terminal domain mutant, R194C, prevents proper folding, whereas the elimination of the cysteine involved in forming the interchain disulfide bond, in a C-terminal domain mutant, C525Y, leads to a truncation at the C terminus. A molecular model of the AMH C-terminal domain provides insights into how some mutations could affect biosynthesis and function.     

Phage P22 tailspike protein: removal of head-binding domain unmasks effects of folding mutations on native-state thermal stability.

A shortened, recombinant protein comprising residues 109-666 of the tailspike endorhamnosidase of Salmonella phage P22 was purified from Escherichia coli and crystallized. Like the full-length tailspike, the protein lacking the amino-terminal head-binding domain is an SDS-resistant, thermostable trimer. Its fluorescence and circular dichroism spectra indicate native structure. Oligosaccharide binding and endoglycosidase activities of both proteins are identical. A number of tailspike folding mutants have been obtained previously in a genetic approach to protein folding. Two temperature-sensitive-folding (tsf) mutations and the four known global second-site suppressor (su) mutations were introduced into the shortened protein and found to reduce or increase folding yields at high temperature. The mutational effects on folding yields and subunit folding kinetics parallel those observed with the full-length protein. They mirror the in vivo phenotypes and are consistent with the substitutions altering the stability of thermolabile folding intermediates. Because full-length and shortened tailspikes aggregate upon thermal denaturation, and their denaturant-induced unfolding displays hysteresis, kinetics of thermal unfolding were measured to assess the stability of the native proteins. Unfolding of the shortened wild-type protein in the presence of 2% SDS at 71 degrees C occurs at a rate of 9.2 x 10(-4) s(-1). It reflects the second kinetic phase of unfolding of the full-length protein. All six mutations were found to affect the thermal stability of the native protein. Both tsf mutations accelerate thermal unfolding about 10-fold. Two of the su mutations retard thermal unfolding up to 5-fold, while the remaining two mutations accelerate unfolding up to 5-fold. The mutational effects can be rationalized on the background of the recently determined crystal structure of the protein.     

Cofactors‐loaded quaternary structure of lysine‐specific demethylase 5C (KDM5C) protein: Computational model

The KDM5C gene (also known as JARID1C and SMCX) is located on the X chromosome and encodes a ubiquitously expressed 1560‐aa protein, which plays an important role in lysine methylation (specifically reverses tri‐ and di‐methylation of Lys4 of histone H3). Currently, 13 missense mutations in KDM5C have been linked to X‐linked mental retardation. However, the molecular mechanism of disease is currently unknown due to the experimental difficulties in expressing such large protein and the lack of experimental 3D structure. In this work, we utilize homology modeling, docking, and experimental data to predict 3D structures of KDM5C domains and their mutual arrangement. The resulting quaternary structure includes KDM5C JmjN, ARID, PHD1, JmjC, ZF domains, substrate histone peptide, enzymatic cofactors, and DNA. The predicted quaternary structure was investigated with molecular dynamic simulation for its stability, and further analysis was carried out to identify features measured experimentally. The predicted structure of KDM5C was used to investigate the effects of disease‐causing mutations and it was shown that the mutations alter domain stability and inter‐domain interactions. The structural model reported in this work could prompt experimental investigations of KDM5C domain‐domain interaction and exploration of undiscovered functionalities. Proteins 2016; 84:1797–1809. © 2016 Wiley Periodicals, Inc.     

Plasticity and steric strain in a parallel beta-helix: rational mutations in the P22 tailspike protein

By means of genetic screens, a great number of mutations that affect the folding and stability of the tailspike protein from Salmonella phage P22 have been identified. Temperature-sensitive folding (tsf) mutations decrease folding yields at high temperature, but hardly affect thermal stability of the native trimeric structure when assembled at low temperature. Global suppressor (su) mutations mitigate this phenotype. Virtually all of these mutations are located in the central domain of tailspike, a large parallel beta-helix. We modified tailspike by rational single amino acid replacements at three sites in order to investigate the influence of mutations of two types: (1) mutations expected to cause a tsf phenotype by increasing the side-chain volume of a core residue, and (2) mutations in a similar structural context as two of the four known su mutations, which have been suggested to stabilize folding intermediates and the native structure by the release of backbone strain, an effect well known for residues that are primarily evolved for function and not for stability or folding of the protein. Analysis of folding yields, refolding kinetics and thermal denaturation kinetics in vitro show that the tsf phenotype can indeed be produced rationally by increasing the volume of side chains in the beta-helix core. The high-resolution crystal structure of mutant T326F proves that structural rearrangements only take place in the remarkably plastic lumen of the beta-helix, leaving the arrangement of the hydrogen-bonded backbone and thus the surface of the protein unaffected. This supports the notion that changes in the stability of an intermediate, in which the beta-helix domain is largely formed, are the essential mechanism by which tsf mutations affect tailspike folding. A rational design of su mutants, on the other hand, appears to be more difficult. The exchange of two residues in the active site expected to lead to a drastic release of steric strain neither enhanced the folding properties nor the stability of tailspike. Apparently, side-chain interactions in these cases overcompensate for backbone strain, illustrating the extreme optimization of the tailspike protein for conformational stability. The result exemplifies the view arising from the statistical analysis of the distribution of backbone dihedral angles in known three-dimensional protein structures that the adoption of straight phi/psi angles other than the most favorable ones is often caused by side-chain interactions. Proteins 2000;39:89-101.     

Experiment and theory highlight role of native state topology in SH3 folding.

We use a combination of experiments, computer simulations and simple model calculations to characterize, first, the folding transition state ensemble of the src SH3 domain, and second, the features of the protein that determine its folding mechanism. Kinetic analysis of mutations at 52 of the 57 residues in the src SH3 domain revealed that the transition state ensemble is even more polarized than suspected earlier: no single alanine substitution in the N-terminal 15 residues or the C-terminal 9 residues has more than a two-fold effect on the folding rate, while such substitutions at 15 sites in the central three-stranded beta-sheet cause significant decreases in the folding rate. Molecular dynamics (MD) unfolding simulations and ab initio folding simulations on the src SH3 domain exhibit a hierarchy of folding similar to that observed in the experiments. The similarity in folding mechanism of different SH3 domains and the similar hierarchy of structure formation observed in the experiments and the simulations can be largely accounted for by a simple native state topology-based model of protein folding energy landscapes.     

Dissecting a bacterial collagen domain from Streptococcus pyogenes: sequence and length-dependent variations in triple helix stability and folding

To better investigate the relationship between sequence, stability, and folding, the Streptococcus pyogenes collagenous domain CL (Gly-Xaa-Yaa)(79) was divided to create three recombinant triple helix subdomains A, B, and C of almost equal size with distinctive amino acid features: an A domain high in polar residues, a B domain containing the highest concentration of Pro residues, and a very highly charged C domain. Each segment was expressed as a monomer, a linear dimer, and a linear trimer fused with the trimerization domain (V domain) in Escherichia coli. All recombinant proteins studied formed stable triple helical structures, but the stability varied depending on the amino acid sequence in the A, B, and C segments and increased as the triple helix got longer. V-AAA was found to melt at a much lower temperature (31.0 °C) than V-ABC (V-CL), whereas V-BBB melted at almost the same temperature (～36-37 °C). When heat-denatured, the V domain enhanced refolding for all of the constructs; however, the folding rate was affected by their amino acid sequences and became reduced for longer constructs. The folding rates of all the other constructs were lower than that of the natural V-ABC protein. Amino acid substitution mutations at all Pro residues in the C fragment dramatically decreased stability but increased the folding rate. These results indicate that the thermostability of the bacterial collagen is dominated by the most stable domain in the same manner as found with eukaryotic collagens.     

Folding and oxidation of the antibody domain C(H)3

The non-covalent homodimer formed by the C-terminal domains of the IgG1 heavy chains (C(H)3) is the simplest naturally occurring model system for studying immunoglobulin folding and assembly. In the native state, the intrachain disulfide bridge, which connects a three-stranded and a four-stranded beta-sheet is buried in the hydrophobic core of the protein. Here, we show that the disulfide bridge is not required for folding and association, since the reduced C(H)3 domain folds to a dimer with defined secondary and tertiary structure. However, the thermodynamic stability of the reduced C(H)3 dimer is much lower than that of the oxidized state. This allows the formation of disulfide bonds either concomitant with folding (starting from the reduced, denatured state) or after folding (starting from the reduced dimer). The analysis of the two processes revealed that, under all conditions investigated, one of the cysteine residues, Cys 86, reacts preferentially with oxidized glutathione to a mixed disulfide that subsequently interacts with the less-reactive second thiol group of the intra-molecular disulfide bond. For folded C(H)3, the second step in the oxidation process is slow. In contrast, starting from the unfolded and reduced protein, the oxidation reaction is faster. However, the overall folding reaction of C(H)3 during oxidative folding is a slow process. Especially, dimerization is slow, compared to the association starting from the denatured oxidized state. This deceleration may be due to misfolded conformations trapped by the disulfide bridge.     

Structural basis for type I and type II deficiencies of antithrombotic plasma protein C: Patterns revealed by three-dimensional molecular modelling of mutations of the protease domain

Familial deficiency of protein C is associated with inherited thrombophilia. To explore how specific missense mutations might cause observed clinical phenotypes, known protein C missense mutations were mapped onto three-dimensional homology models of the protein C protease domain, and the implications for domain folding and structure were evaluated. Most Type I missense mutations either replaced internal hydrophobic residues (I201T, L223F, A259V, A267T, A346T, A346V, G376D) or nearby interacting residues (I403M, T298M, Q184H), thus disrupting the packing of internal hydrophobic side chains, or changed hydrophilic residues, thus disrupting ion pairs (N256D, R178W). Mutations (P168L, R169W) at the activation site destabilized the region containing the activation peptide structure. Most Type II mutations involved solvent-exposed residues and were clustered either in a positively charged region (R147W, R157Q, R229Q, R352W) or were located in or near the active site region (S252N, D359N, G381S, G391S, H211Q). The cluster of arginines 147, 157, 229, and 352 may identify a functionally important exosite. Identification of the spatial relationships of natural mutations in the protein C model is helpful for understanding manifestations of protein C deficiency and for identification of novel, functionally important molecular features and exosites. © 1994 John Wiley & Sons, Inc.     

Folding of the multidomain ribosomal protein L9: the two domains fold independently with remarkably different rates

The folding and unfolding behavior of the multidomain ribosomal protein L9 from Bacillus stearothermophilus was studied by a novel combination of stopped-flow fluorescence and nuclear magnetic resonance (NMR) spectroscopy. One-dimensional 1H spectra acquired at various temperatures show that the C-terminal domain unfolds at a lower temperature than the N-terminal domain (Tm = 67 degrees C for the C-terminal domain, 80 degrees C for the N-terminal domain). NMR line-shape analysis was used to determine the folding and unfolding rates for the N-terminal domain. At 72 degrees C, the folding rate constant equals 2980 s-1 and the unfolding rate constant equals 640 s-1. For the C-terminal domain, saturation transfer experiments performed at 69 degrees C were used to determine the folding rate constant, 3.3 s-1, and the unfolding rate constant, 9.0 s-1. Stopped-flow fluorescence experiments detected two resolved phases: a fast phase for the N-terminal domain and a slow phase for the C-terminal domain. The folding and unfolding rate constants determined by stopped-flow fluorescence are 760 s-1 and 0.36 s-1, respectively, for the N-terminal domain at 25 degrees C and 3.0 s-1 and 0.0025 s-1 for the C-terminal domain. The Chevron plots for both domains show a V-shaped curve that is indicative of two-state folding. The measured folding rate constants for the N-terminal domain in the intact protein are very similar to the values determined for the isolated N-terminal domain, demonstrating that the folding kinetics of this domain is not affected by the rest of the protein. The remarkably different rate constants between the N- and C-terminal domains suggest that the two domains can fold and unfold independently. The folding behavior of L9 argues that extremely rapid folding is not necessarily functionally important.     

The R1441C mutation alters the folding properties of the ROC domain of LRRK2

LRRK2 is a 250 kDa multidomain protein, mutations in which cause familial Parkinson's disease. Previously, we have demonstrated that the R1441C mutation in the ROC domain decreases GTPase activity. Here we show that the R1441C alters the folding properties of the ROC domain, lowering its thermodynamic stability. Similar to small GTPases, binding of different guanosine nucleotides alters the stability of the ROC domain, suggesting that there is an alteration in conformation dependent on GDP or GTP occupying the active site. GTP/GDP bound state also alters the self-interaction of the ROC domain, accentuating the impact of the R1441C mutation on this property. These data suggest a mechanism whereby the R1441C mutation can reduce the GTPase activity of LRRK2, and highlights the possibility of targeting the stability of the ROC domain as a therapeutic avenue in LRRK2 disease.     

Formation of the VHL-elongin BC tumor suppressor complex is mediated by the chaperonin TRiC

von Hippel-Lindau (VHL) disease is caused by loss of function of the VHL tumor suppressor protein. Here, we demonstrate that the folding and assembly of VHL into a complex with its partner proteins, elongin B and elongin C (herein, elongin BC), is directly mediated by the chaperonin TRiC/CCT. Association of VHL with TRiC is required for formation of the VHL-elongin BC complex. A 55-amino acid domain of VHL is both necessary and sufficient for binding to TRiC. Importantly, mutation or deletion of this domain is associated with VHL disease. We identified two mutations that disrupt the normal interaction with TRiC and impair VHL folding. Our results define a novel role for TRiC in mediating oligomerization and suggest that inactivating mutations can impair polypeptide function by interfering with chaperone-mediated folding.     

Isolation and characterization of temperature-sensitive and thermostable mutants of the human receptor-like protein tyrosine phosphatase LAR.

Human LAR is a transmembrane receptor-like protein whose cytoplasmic region contains two tandemly duplicated domains homologous to protein tyrosine phosphatases (PTPases). Whereas the membrane-proximal domain I has enzymatic activity, the membrane-distal domain II has no apparent catalytic activity but seems to have a regulatory function. In order to study structure-function relationships of the LAR PTPase, LAR domain I was expressed in Escherichia coli, and mutants that have reduced catalytic activity or reduced thermostability were isolated and characterized. We isolated 18 unique hydroxylamine-induced missense mutations in the LAR domain I segment, of which three were temperature-sensitive. Five additional temperature-sensitive mutations were isolated using N-methyl-N'-nitro-N-nitrosoguanidine. All eight temperature-sensitive mutations are confined within a short segment of the LAR domain I sequence between amino acid positions 1329 and 1407. To examine whether this region is particularly prone to temperature-sensitive mutations, tyrosine at amino acid position 1379 was changed to a phenylalanine by oligonucleotide-directed mutagenesis. This mutant, Y1379-F, was indeed temperature-sensitive. We also isolated a revertant of a temperature-sensitive mutant. The revertant contained a second-site mutation (C1446-Y) that suppresses several temperature-sensitive mutations and also enhances the folding of LAR protein produced in E. coli.     

Molecular recognition via coupled folding and binding in a TPR domain

The majority of known tetratricopeptide repeat (TPR) domains consist of three copies of the helix-turn-helix TPR motif, together with a seventh C-terminal helix. TPR domains function as protein-protein recognition modules in intracellular signalling. This function is exemplified by the TPR domain of protein phosphatase 5 (PP5), which binds to the C terminus of the chaperone protein Hsp90. Here, we report NMR and CD spectroscopic studies that reveal that this domain is largely unfolded at physiological temperatures, and that interaction with an MEEVD pentapeptide derived from Hsp90 stabilises a folded structure. This complex, coupled folding-binding mechanism is characterised further by its observed enthalpy change on binding (determined by isothermal titration calorimetry), which displays a markedly non-linear relationship with temperature. A nested Gibbs-Helmholtz model is used in a novel combined analysis of the CD and ITC data to determine separately the thermodynamic contributions of the intrinsic folding and binding events to the overall coupled process. The analysis shows that, despite the expected large entropic opposition to the folding process, a nearly equal favourable folding enthalpy means the net effect of coupled folding on the observed affinity is small across a broad range of temperature. We hypothesise that a coupled folding-binding mechanism is common in this class of domains.     

Mapping the folding pathway of an immunoglobulin domain: structural detail from Phi value analysis and movement of the transition state

BACKGROUND: Do proteins that have the same structure fold by the same pathway even when they are unrelated in sequence? To address this question, we are comparing the folding of a number of different immunoglobulin-like proteins. Here, we present a detailed protein engineering phi value analysis of the folding pathway of TI I27, an immunoglobulin domain from human cardiac titin. RESULTS: TI I27 folds rapidly via a kinetic intermediate that is destabilized by most mutations. The transition state for folding is remarkably native-like in terms of solvent accessibility. We use phi value analysis to map this transition state and show that it is highly structured; only a few residues close to the N-terminal region of the protein remain completely unfolded. Interestingly, most mutations cause the transition state to become less native-like. This anti-Hammond behavior can be used as a novel means of obtaining additional structural information about the transition state. CONCLUSIONS: The residues that are involved in nucleating the folding of TI I27 are structurally equivalent to the residues that form the folding nucleus in an evolutionary unrelated fibronectin type III protein. These residues form part of the common structural core of Ig-like domains. The data support the hypothesis that interactions essential for defining the structure of these beta sandwich proteins are also important in nucleation of folding.