Unveiling the Unfolding Pathway of F5F8D Disorder-Associated D81H/V100D Mutant of MCFD2 via Multiple Molecular Dynamics Simulations

Abstract

Combined factor deficiency (F5F8D) is a rare autosomal recessive disorder caused by mutations in the LMAN1 or MCFD2 genes. It has been proposed that this pathogenic process occurs via a multi-step pathway involving metal loss, EF-hand-Ca²⁺ dissociation and assembly of misfolded MCFD2-LMAN1 complex. Here, we have investigated the solution conformations of the MCFD2^(D81H,V100D) protein mutant through extensive molecular dynamics (MD) simulations. The V100D, one of the many MCFD2 mutations known to be associated to F5F8D, is difficult to be reconciled with the pathway model because it is located far from the metal sites and the MCFD2/LMAN1 interface. Consequently, an inspection of all the steps involved in D81H/V100D MCFD2 misfolding is expected to provide hints in the understanding of the molecular basis of the disease. A comparison with parallel studies carried out for the Wild-Type (WT) MCFD2 pointed out that the mutation decreases the affinity of the protein for the Ca²⁺ ion. Multiple explicit solvents MD simulations (50 ns) performed on the two proteins revealed that in the WT protein, stable H-bond network and compact hydrophobic core region are created thus confirming a pivotal role of this region in driving the biophysical properties of the entire protein. In fact it is shown that the V100D mutation, although located far away the EF-hand domain, may induce subtle modification in the structural core of MCFD2 leading to the loosening of metal binding and to the formation of metastable intermediate states along the unfolding pathway. The native-like hydrophobic cluster formed near the V100 residue in the wild-type protein is disrupted by the negatively charged Asparagine residue. Furthermore, the presence of the D81H mutation in the EF-1 hand domain may also increase the protein unfolding rate and consequently prevent the formation of the MCFD2-LMAN1 complex. The detailed structural insights obtained from our large-scale simulations complement the clinical features and offer useful insights into the mechanism behind MCFD2 protein misfolding.     

Crystal structure of the LMAN1-CRD/MCFD2 transport receptor complex provides insight into combined deficiency of factor V and factor VIII

LMAN1 is a glycoprotein receptor, mediating transfer from the ER to the ER-Golgi intermediate compartment. Together with the co-receptor MCFD2, it transports coagulation factors V and VIII. Mutations in LMAN1 and MCFD2 can cause combined deficiency of factors V and VIII (F5F8D). We present the crystal structure of the LMAN1/MCFD2 complex and relate it to patient mutations. Circular dichroism data show that the majority of the substitution mutations give rise to a disordered or severely destabilized MCFD2 protein. The few stable mutation variants are found in the binding surface of the complex leading to impaired LMAN1 binding and F5F8D.

Structured summary

MINT-7557086: lman1 (uniprotkb:P49257) and mcfd2 (uniprotkb:Q8NI22) bind (MI:0407) by X-ray crystallography (MI:0114)     

Structural Characterization of Carbohydrate Binding by LMAN1 Protein Provides New Insight into the Endoplasmic Reticulum Export of Factors V (FV) and VIII (FVIII)

LMAN1 (ERGIC-53) is a key mammalian cargo receptor responsible for the export of a subset of glycoproteins from the endoplasmic reticulum. Together with its soluble coreceptor MCFD2, LMAN1 transports coagulation factors V (FV) and VIII (FVIII). Mutations in LMAN1 or MCFD2 cause the genetic bleeding disorder combined deficiency of FV and FVIII (F5F8D). The LMAN1 carbohydrate recognition domain (CRD) binds to both glycoprotein cargo and MCFD2 in a Ca²⁺-dependent manner. To understand the biochemical basis and regulation of LMAN1 binding to glycoprotein cargo, we solved crystal structures of the LMAN1-CRD bound to Man-α-1,2-Man, the terminal carbohydrate moiety of high mannose glycans. Our structural data, combined with mutagenesis and in vitro binding assays, define the central mannose-binding site on LMAN1 and pinpoint histidine 178 and glycines 251/252 as critical residues for FV/FVIII binding. We also show that mannobiose binding is relatively independent of pH in the range relevant for endoplasmic reticulum-to-Golgi traffic, but is sensitive to lowered Ca²⁺ concentrations. The distinct LMAN1/MCFD2 interaction is maintained at these lowered Ca²⁺ concentrations. Our results suggest that compartmental changes in Ca²⁺ concentration regulate glycoprotein cargo binding and release from the LMAN1·MCFD2 complex in the early secretory pathway.     

LMAN1 and MCFD2 form a cargo receptor complex and interact with coagulation factor VIII in the early secretory pathway

Mutations in LMAN1 (ERGIC-53) and MCFD2 are the causes of a human genetic disorder, combined deficiency of coagulation factor V and factor VIII. LMAN1 is a type 1 transmembrane protein with homology to mannose-binding lectins. MCFD2 is a soluble EF-hand-containing protein that is retained in the endoplasmic reticulum through its interaction with LMAN1. We showed that endogenous LMAN1 and MCFD2 are present primarily in complex with each other with a 1:1 stoichiometry, although MCFD2 is not required for oligomerization of LMAN1. Using a cross-linking-immunoprecipitation assay, we detected a specific interaction of both LMAN1 and MCFD2 with factor VIII, with the B domain as the most likely site of interaction. We also present evidence that this interaction is independent of the glycosylation state of factor VIII but requires native calcium concentration in the endoplasmic reticulum. The interaction of MCFD2 with factor VIII appeared to be independent of LMAN1-MCFD2 complex formation. These results suggest that LMAN1 and MCFD2 form a cargo receptor complex and that the primary sorting signals residing in the B domain direct the binding of factor VIII to LMAN1-MCFD2 through calcium-dependent protein-protein interactions. MCFD2 may function to specifically recruit factor V and factor VIII to sites of transport vesicle budding within the endoplasmic reticulum lumen.     

Role of hydrophobic core on the thermal stability of proteins - molecular dynamics simulations on a single point mutant of Sso7d abstract

The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 s long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Role of Hydrophobic Core on the Thermal Stability of Proteins—Molecular Dynamics Simulations on a Single Point Mutant of Sso7d

Abstract

The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 μs long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

New insights into multiple coagulation factor deficiency from the solution structure of human MCFD2

Human MCFD2 (multiple coagulation factor deficiency 2) is a 16-kDa protein known to participate in transport of the glycosylated human coagulation factors V and VIII along the secretory pathway. Mutations in MCFD2 or in its binding partner, the membrane-bound transporter ERGIC (endoplasmic reticulum-Golgi intermediate compartment)-53, cause a mild form of inherited hemophilia known as combined deficiency of factors V and VIII (F5F8D). While ERGIC-53 is known to be a lectin-type mannose binding protein, the role of MCFD2 in the secretory pathway is comparatively unclear. MCFD2 has been shown to bind both ERGIC-53 and the blood coagulation factors, but little is known about the binding sites or the true function of the protein. In order to facilitate understanding of the function of MCFD2 and the mechanism by which mutations in the protein cause F5F8D, we have determined the structure of human MCFD2 in solution by NMR. Our results show the folding of MCFD2 to be dependent on availability of calcium ions. The protein, which is disordered in the apo state, folds upon binding of Ca²⁺ to the two EF-hand motifs of its C-terminus, while retaining some localized disorder in the N-terminus. NMR studies on two disease-causing mutant variants of MCFD2 show both to be predominantly disordered, even in the presence of calcium ions. These results provide an explanation for the previously observed calcium dependence of the MCFD2-ERGIC-53 interaction and, furthermore, clarify the means by which mutations in this protein result in inefficient secretion of blood coagulation factors V and VIII.     

Kinetic stability of Cu/Zn superoxide dismutase is dependent on its metal ligands: implications for ALS

Over 100 mutants of the enzyme Cu/Zn superoxide dismutase (SOD) have been implicated in the neurodegenerative disease familial amyotrophic lateral sclerosis (FALS). Growing evidence suggests that the aggregation of SOD mutants may play a causative role in FALS and that aberrant copper chemistry, decreased thermodynamic stability, and decreased affinity for metals may contribute independently or synergistically to this process. Since the loss of the copper and zinc ions significantly decreases the thermodynamic stability of SOD, it is expected that this would also decrease its kinetic stability, thereby facilitating partial or global unfolding transitions that may lead to misfolding and aggregation. Here we used wild-type (WT) SOD and five FALS-related mutants (G37R, H46R, G85R, D90A, and L144F) to show that the metals contribute significantly to the kinetic stability of the protein, with demetalated (apo) SOD showing acid-induced unfolding rates about 60-fold greater than the metalated (holo) protein. However, the unfolding rates of SOD WT and mutants were similar to each other in both the holo and apo states, indicating that regardless of the effect of mutation on thermodynamic stability, the kinetic barrier toward SOD unfolding is dependent on the presence of metals. Thus, these results suggest that pathogenic SOD mutations that do not significantly alter the stability of the protein may still lead to SOD aggregation by compromising its ability to bind or retain its metals and thereby decrease its kinetic stability. Furthermore, the mutant-like decrease in the kinetic stability of apo WT SOD raises the possibility that the loss of metals in WT SOD may be involved in nonfamilial forms of ALS.     

Role of hydrophobic core on the thermal stability of proteins-molecular dynamics simulations on a single point mutant of sso7d

Abstract The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 μs long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Effects of helix and fingertip mutations on the thermostability of xyn11A investigated by molecular dynamics simulations and enzyme activity assays

Local conformational changes and global unfolding pathways of wildtype xyn11A recombinant and its mutated structures were studied through a series of atomistic molecular dynamics (MD) simulations, along with enzyme activity assays at three incubation temperatures to investigate the effects of mutations at three different sites to the thermostability. The first mutation was to replace an unstable negatively charged residue at a surface beta turn near the active site (D32G) by a hydrophobic residue. The second mutation was to create a disulphide bond (S100C/N147C) establishing a strong connection between an alpha helix and a distal beta hairpin associated with the thermally sensitive Thumb loop, and the third mutation add an extra hydrogen bond (A155S) to the same alpha helix. From the MD simulations performed, MM/PBSA energy calculations of the unfolding energy were in a good agreement with the enzyme activities measured from the experiment, as all mutated structures demonstrated the improved thermostability, especially the S100C/N147C proved to be the most stable mutant both by the simulations and the experiment. Local conformational analysis at the catalytic sites and the xylan access region also suggested that mutated xyn11A structures could accommodate xylan binding. However, the analysis of global unfolding pathways showed that structural disruptions at the beta sheet regions near the N-terminal were still imminent. These findings could provide the insight on the molecular mechanisms underlying the enhanced thermostability due to mutagenesis and changes in the protein unfolding pathways for further protein engineering of the GH11 family xylanase enzymes.     

A hidden aggregation‐prone structure in the heart of hypoxia inducible factor prolyl hydroxylase

Prolyl hydroxylase domain‐containing protein 2 (PHD2), as one of the most important regulators of angiogenesis and metastasis of cancer cells, is a promising target for cancer therapy drug design. Progressive studies imply that abnormality in PHD2 function may be due to misfolding. Therefore, study of the PHD2 unfolding pathway paves the way for a better understanding of the influence of PHD2 mutations and cancer cell metabolites on the protein folding pathway. We study the unfolding of the PHD2 catalytic domain using differential scanning calorimetry (DSC), fluorescence spectroscopy, and discrete molecular dynamics simulations (DMD). Using computational and experimental techniques, we find that PHD2 undergoes four transitions along the thermal unfolding pathway. To illustrate PHD2 unfolding events in atomic detail, we utilize DMD simulations. Analysis of computational results indicates an intermediate species in the PHD2 unfolding pathway that may enhance aggregation propensity, explaining mutation‐independent PHD2 malfunction. Proteins 2016; 84:611–623. © 2016 Wiley Periodicals, Inc.     

Familial amyotrophic lateral sclerosis-associated mutations decrease the thermal stability of distinctly metallated species of human copper/zinc superoxide dismutase

We report the thermal stability of wild type (WT) and 14 different variants of human copper/zinc superoxide dismutase (SOD1) associated with familial amyotrophic lateral sclerosis (FALS). Multiple endothermic unfolding transitions were observed by differential scanning calorimetry for partially metallated SOD1 enzymes isolated from a baculovirus system. We correlated the metal ion contents of SOD1 variants with the occurrence of distinct melting transitions. Altered thermal stability upon reduction of copper with dithionite identified transitions resulting from the unfolding of copper-containing SOD1 species. We demonstrated that copper or zinc binding to a subset of "WT-like" FALS mutants (A4V, L38V, G41S, G72S, D76Y, D90A, G93A, and E133Delta) conferred a similar degree of incremental stabilization as did metal ion binding to WT SOD1. However, these mutants were all destabilized by approximately 1-6 degrees C compared with the corresponding WT SOD1 species. Most of the "metal binding region" FALS mutants (H46R, G85R, D124V, D125H, and S134N) exhibited transitions that probably resulted from unfolding of metal-free species at approximately 4-12 degrees C below the observed melting of the least stable WT species. We conclude that decreased conformational stability shared by all of these mutant SOD1s may contribute to SOD1 toxicity in FALS.     

A shift in the salt bridge interaction of residues D620 and E621 mediates the constitutive activation of Jak2-H538Q/K539L

Jak2 mutations in the exon 14 and exon 12 regions that cause constitutive activation have been associated with myeloproliferative neoplasms. We have previously shown that a pi stacking interaction between F617 and F595 is important for the constitutive activation of Jak2-V617F (Gnanasambandan et al., Biochemistry 49:9972-9984, 2010). Here, using a combination of molecular dynamics (MD) simulations and in vitro mutagenesis, we studied the molecular mechanism for the constitutive activation of the Jak2 exon 12 mutation, H538Q/K539L. The activation levels of Jak2-H538Q/K539L were found to be similar to that of Jak2-V617F, and Jak2-H538Q/K539L/V617F. Data from MD simulations indicated a shift in the salt bridge interactions of D620 and E621 with K539 in Jak2-WT to R541 in Jak2-H538Q/K539L. When compared to Jak2-WT, K539A mutation resulted in increased activation, while K539D or K539E mutations diminished Jak2 activation by 50 %. In the context of Jak2-H538Q/K539L, R541A mutation reduced its activation by 50 %, while R541D and R541E mutations returned its activation levels to that of Jak2-WT. Collectively, these results indicate that a shift in the salt bridge interaction of D620 and E621 with K539 in Jak2-WT to R541 in Jak2-H538Q/K539L is critical for constitutive activation of this Jak2 exon 12 mutant.     

HIV-1整合酶核心区野生型和F185K突变型活性和溶解性的比较及分子动力学模拟分析

表达纯化了野生型(WT)及F185K突变型HIV-1整合酶核心区蛋白(IN_C),并对二者的溶解性和活性进行了比较．实验结果表明：F185K 突变后IN_C溶解性显著提高,活性有一定程度降低．对WT和F185K IN_C体系进行了1800 ps的分子动力学模拟．模拟结果表明：F185K IN_C功能loop区柔性和蛋白质整体运动性降低,使蛋白质活性降低,F185K突变后盐桥网络的变化驱动了IN_C局部构象改变,引起IN_C表面的部分疏水残基被包埋,亲水残基暴露,相对亲水溶剂可接近面积增大,同时,突变后IN_C与水之间形成氢键的数量增加,与水之间作用加强,以上变化使IN_C溶解性提高．分子动力学模拟与实验结果相吻合．为理解蛋白质溶解性和对蛋白质进行可溶性改造提供了一定的理论依据．     

Different misfolding mechanisms converge on common conformational changes

Prion diseases are caused by misfolding and aggregation of the prion protein (PrP). Pathogenic mutations such as Y218N and E196K are known to cause Gerstmann-Sträussler-Scheinker syndrome and Creutzfeldt-Jakob disease, respectively. Here we describe molecular dynamics simulations of these mutant proteins to better characterize the detailed conformational effects of these sequence substitutions. Our results indicate that the mutations disrupt the wild-type native PrP^C structure and cause misfolding. Y218N reduced hydrophobic packing around the X-loop (residues 165–171), and E196K abolished an important wild-type salt bridge. While differences in the mutation site led PrP mutants to misfold along different pathways, we observed multiple traits of misfolding that were common to both mutants. Common traits of misfolding included: 1) detachment of the short helix (HA) from the PrP core; 2) exposure of side chain F198; and 3) formation of a nonnative strand at the N-terminus. The effect of the E196K mutation directly abolished the wild-type salt bridge E196-R156, which further destabilized the F198 hydrophobic pocket and HA. The Y218N mutation propagated its effect by increasing the HB-HC interhelical angle, which in turn disrupted the packing around F198. Furthermore, a nonnative contact formed between E221 and S132 on the S1-HA loop, which offered a direct mechanism for disrupting the hydrophobic packing between the S1-HA loop and HC. While there were common misfolding features shared between Y218N and E196K, the differences in the orientation of HB and HC and the X-loop conformation might provide a structural basis for identifying different prion strains.     

Hydrophobic Core Mutations Associated with Cataract Development in Mice Destabilize Human ��D-Crystallin

The human eye lens is composed of fiber cells packed with crystallins up to 450 mg/ml. Human γD-crystallin (HγD-Crys) is a monomeric, two-domain protein of the lens central nucleus. Both domains of this long lived protein have double Greek key β-sheet folds with well packed hydrophobic cores. Three mutations resulting in amino acid substitutions in the γ-crystallin buried cores (two in the N-terminal domain (N-td) and one in the C-terminal domain (C-td)) cause early onset cataract in mice, presumably an aggregated state of the mutant crystallins. It has not been possible to identify the aggregating precursor within lens tissues. To compare in vivo cataract-forming phenotypes with in vitro unfolding and aggregation of γ-crystallins, mouse mutant substitutions were introduced into HγD-Crys. The mutant proteins L5S, V75D, and I90F were expressed and purified from Escherichia coli. WT HγD-Crys unfolds in vitro through a three-state pathway, exhibiting an intermediate with the N-td unfolded and the C-td native-like. L5S and V75D in the N-td also displayed three-state unfolding transitions, with the first transition, unfolding of the N-td, shifted to significantly lower denaturant concentrations. I90F destabilized the C-td, shifting the overall unfolding transition to lower denaturant concentrations. During thermal denaturation, the mutant proteins exhibited lowered thermal stability compared with WT. Kinetic unfolding experiments showed that the N-tds of L5S and V75D unfolded faster than WT. I90F was globally destabilized and unfolded more rapidly. These results support models of cataract formation in which generation of partially unfolded species are precursors to the aggregated cataractous states responsible for light scattering.     

Free energies of binding of <Emphasis Type="Italic">R</Emphasis>- and <Emphasis Type="Italic">S</Emphasis>-propranolol to wild-type and F483A mutant cytochrome P450 2D6 from molecular dynamics simulations

Detailed molecular dynamics (MD) simulations have been performed to reproduce and rationalize the experimental finding that the F483A mutant of CYP2D6 has lower affinity for R-propranolol than for S-propranolol. Wild-type (WT) CYP2D6 does not show this stereospecificity. Four different approaches to calculate the free energy differences have been investigated and were compared to the experimental binding data. From the differences between calculations based on forward and backward processes and the closure of thermodynamic cycles, it was clear that not all simulations converged sufficiently. The approach that calculates the free energies of exchanging R-propranolol with S-propranolol in the F483A mutant relative to the exchange free energy in WT CYP2D6 accurately reproduced the experimental binding data. Careful inspection of the end-points of the MD simulations involved in this approach, allowed for a molecular interpretation of the observed differences. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation

The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.