Computational diagnosis of protein conformational diseases: Short molecular dynamics simulations reveal a fast unfolding of r-LDL mutants that cause familial hypercholesterolemia

Interpreting the Mechanism of APOE (p.Leu167del) Mutation in the Incidence of Familial Hypercholesterolemia; An In-silico Approach

Background:

Apolipoprotein E (APOE) gene is a ligand protein in humans which mediates the metabolism of cholesterol by binding to the low-density lipoprotein receptor (LDLR). P.Leu167del mutation in APOE gene was recently connected with Familial Hypercholesterolemia, a condition associated with premature cardiovascular disease. The consequences of this mutation on the protein structure and its receptor binding capacity remain largely unknown.

Objective:

The current study aims to further decipher the underlying mechanism of this mutation using advanced software-based algorithms. The consequences of disrupting the leucine zipper by this mutation was studied at the structural and functional level of the APOE protein.

Methods:

3D protein modeling for both APOE and LDLR (wild types), along with APOE (p.Leu167del) mutant type were generated using homology modeling template-based alignment. Structural deviation analysis was performed to evaluate the spatial orientation and the stability of the mutant APOE structure. Molecular docking analysis simulating APOE-LDLR protein interaction was carried out, in order to evaluate the impact of the mutation on the binding affinity.

Result:

Structural deviation analysis for APOE mutated model showed low degree of deviance scoring root-mean-square deviation, (RMSD) = 0.322 Å. Whereas Docking simulation revealed an enhanced molecular interaction towards the LDLR with an estimation of +171.03 kJ/mol difference in binding free energy.

Conclusion:

This in-silico study suggests that p.Leu167del is causing the protein APOE to associate strongly with its receptor, LDLR. This gain-of-function is likely hindering the ability of LDLR to be effectively recycled back to the surface of the hepatocytes to clear cholesterol from the circulation therefore leading to FH.

Exploring the complete mutational space of the LDL receptor LA5 domain using molecular dynamics: linking SNPs with disease phenotypes in familial hypercholesterolemia

Familial hypercholesterolemia (FH), a genetic disorder with a prevalence of 0.2%, represents a high-risk factor to develop cardiovascular and cerebrovascular diseases. The majority and most severe FH cases are associated to mutations in the receptor for low-density lipoproteins receptor (LDL-r), but the molecular basis explaining the connection between mutation and phenotype is often unknown, which hinders early diagnosis and treatment of the disease. We have used atomistic simulations to explore the complete SNP mutational space (227 mutants) of the LA5 repeat, the key domain for interacting with LDL that is coded in the exon concentrating the highest number of mutations. Four clusters of mutants of different stability have been identified. The majority of the 50 FH known mutations (33) appear distributed in the unstable clusters, i.e. loss of conformational stability explains two-third of FH phenotypes. However, one-third of FH phenotypes (17 mutations) do not destabilize the LR5 repeat. Combining our simulations with available structural data from different laboratories, we have defined a consensus-binding site for the interaction of the LA5 repeat with LDL-r partner proteins and have found that most (16) of the 17 stable FH mutations occur at binding site residues. Thus, LA5-associated FH arises from mutations that cause either the loss of stability or a decrease in domain's-binding affinity. Based on this finding, we propose the likely phenotype of each possible SNP in the LA5 repeat and outline a procedure to make a full computational diagnosis for FH.

Sorting an LDL receptor with bound PCSK9 to intracellular degradation

OBJECTIVE

This article reviews the mechanism by which the low density lipoprotein receptor (LDLR) that has bound proprotein convertase subtilisin/kexin type 9 (PCSK9), is rerouted to intracellular degradation instead of being recycled.

METHODS

A search of relevant published literature has been conducted.

RESULTS

PCSK9 binds to the LDLR at the cell surface. It is the catalytic domain of PCSK9 that binds to the epidermal growth factor repeat A of the LDLR. The LDLR:PCSK9 complex is internalized through clathrin-mediated endocytosis. Due to an additional electrostatic interaction at acidic pH between the C-terminal domain of PCSK9 and the ligand-binding domain of the LDLR, PCSK9 remains bound to the LDLR in the sorting endosome. As a consequence, the LDLR fails to adopt a closed conformation and is degraded instead of being recycled. The mechanism for the failure of the LDLR to recycle appears to involve ectodomain cleavage of the extended LDLR by a cysteine cathepsin in the sorting endosome. The cleaved LDLR ectodomain will be confined to the vesicular part of the sorting endosome for degradation in the endosomal/lysosomal tract.

CONCLUSION

Ectodomain cleavage of an LDLR with bound PCSK9 in the sorting endosome disrupts the normal recycling of the LDLR.

PCSK9-mediated degradation of the LDL receptor generates a 17 kDa C-terminal LDL receptor fragment

Proprotein convertase subtilisin/kexin type 9 (PCSK9) binds to the LDL receptor (LDLR) at the cell surface and reroutes the internalized LDLR to intracellular degradation. In this study, we have shown that PCSK9-mediated degradation of the full-length 160 kDa LDLR generates a 17 kDa C-terminal LDLR fragment. This fragment was not generated from mutant LDLRs resistant to PCSK9-mediated degradation or when degradation was prevented by chemicals such as ammonium chloride or the cysteine cathepsin inhibitor E64d. The observation that the 17 kDa fragment was only detected when the cells were cultured in the presence of the γ-secretase inhibitor DAPT indicates that this 17 kDa fragment undergoes γ-secretase cleavage within the transmembrane domain. The failure to detect the complementary 143 kDa ectodomain fragment is likely to be due to its rapid degradation in the endosomal lumen. The 17 kDa C-terminal LDLR fragment was also generated from a Class 5 mutant LDLR undergoing intracellular degradation. Thus, one may speculate that an LDLR with bound PCSK9 and a Class 5 LDLR with bound LDL are degraded by a similar mechanism that could involve ectodomain cleavage in the endosome.

Quantitative fluorescence imaging reveals point of release for lipoproteins during LDLR-dependent uptake

The LDL receptor (LDLR) supports efficient uptake of both LDL and VLDL remnants by binding lipoprotein at the cell surface, internalizing lipoprotein through coated pits, and releasing lipoprotein in endocytic compartments before returning to the surface for further rounds of uptake. While many aspects of lipoprotein binding and receptor entry are well understood, it is less clear where, when, and how the LDLR releases lipoprotein. To address these questions, the current study employed quantitative fluorescence imaging to visualize the uptake and endosomal processing of LDL and the VLDL remnant β-VLDL. We find that lipoprotein release is rapid, with most release occurring prior to entry of lipoprotein into early endosomes. Published biochemical studies have identified two mechanisms of lipoprotein release: one that involves the β-propeller module of the LDLR and a second that is independent of this module. Quantitative imaging comparing uptake supported by the normal LDLR or by an LDLR variant incapable of β-propeller-dependent release shows that the β-propeller-independent process is sufficient for release for both lipoproteins but that the β-propeller process accelerates both LDL and β-VLDL release. Together these findings define where, when, and how lipoprotein release occurs and provide a generalizable methodology for visualizing endocytic handling in situ.

Structure and folding of disulfide-rich miniproteins: insights from molecular dynamics simulations and MM-PBSA free energy calculations

The fold of small disulfide-rich proteins largely relies on two or more disulfide bridges that are main components of the hydrophobic core. Because of the small size of these proteins and their high cystine content, the cysteine connectivity has been difficult to ascertain in some cases, leading to uncertainties and debates in the literature. Here, we use molecular dynamics simulations and MM-PBSA free energy calculations to compare similar folds with different disulfide pairings in two disulfide-rich miniprotein families, namely the knottins and the short-chain scorpion toxins, for which the connectivity has been discussed. We first show that the MM-PBSA approach is able to discriminate the correct knotted topology of knottins from the laddered one. Interestingly, a comparison of the free energy components for kalata B1 and MCoTI-II suggests that cyclotides and squash inhibitors, although sharing the same scaffold, are stabilized through different interactions. Application to short-chain scorpion toxins suggests that the conventional cysteine pairing found in many homologous toxins is significantly more stable than the unconventional pairing reported for maurotoxin and for spinoxin. This would mean that native maurotoxin and spinoxin are not at the lowest free energy minimum and might result from kinetically rather than thermodynamically driven oxidative folding processes. For both knottins and toxins, the correct or conventional disulfide connectivities provide lower flexibilities and smaller deviations from the initial conformations. Overall, our work suggests that molecular dynamics simulations and the MM-PBSA approach to estimate free energies are useful tools to analyze and compare disulfide bridge connectivities in miniproteins.

A novel 9-bp insertion detected in steroid 21-hydroxylase gene (CYP21A2): prediction of its structural and functional implications by computational methods

BACKGROUND

Steroid 21-hydroxylase deficiency is the most common cause of congenital adrenal hyperplasia (CAH). Detection of underlying mutations in CYP21A2 gene encoding steroid 21-hydroxylase enzyme is helpful both for confirmation of diagnosis and management of CAH patients. Here we report a novel 9-bp insertion in CYP21A2 gene and its structural and functional consequences on P450c21 protein by molecular modeling and molecular dynamics simulations methods.

METHODS

A 30-day-old child was referred to our laboratory for molecular diagnosis of CAH. Sequencing of the entire CYP21A2 gene revealed a novel insertion (duplication) of 9-bp in exon 2 of one allele and a well-known mutation I172N in exon 4 of other allele. Molecular modeling and simulation studies were carried out to understand the plausible structural and functional implications caused by the novel mutation.

RESULTS

Insertion of the nine bases in exon 2 resulted in addition of three valine residues at codon 71 of the P450c21 protein. Molecular dynamics simulations revealed that the mutant exhibits a faster unfolding kinetics and an overall destabilization of the structure due to the triple valine insertion was also observed.

CONCLUSION

The novel 9-bp insertion in exon 2 of CYP21A2 genesignificantly lowers the structural stability of P450c21 thereby leading to the probable loss of its function.

Mechanism of low density lipoprotein (LDL) release in the endosome: implications of the stability and Ca2+ affinity of the fifth binding module of the LDL receptor

Uptake of low density lipoproteins (LDL) by their receptor, LDLR, is the primary mechanism by which cells incorporate cholesterol from plasma. Mutations in LDLR lead to familial hypercholesterolemia, a common disease affecting 1 in 500 of the human population. LDLR is a modular protein that uses several small repeats to bind LDL. The repeats contain around 40 residues, including three disulfide bonds and a calcium ion. Repeat 5 (LR5) is critical for LDL and beta-migrating very low density lipoprotein binding. Based on the crystal structure of LDLR at endosomal pH (but close to extracellular calcium concentration), LR5 has been proposed to bind to the epidermal growth factor (EGF) precursor domain of LDLR in the endosome, thus releasing the LDL particles previously bound in extracellular conditions. We report here the conformational stability of LR5 as a function of temperature and calcium concentration under both extracellular and endosomal pH conditions. The repeat was very stable when it bore a bound calcium ion but was severely destabilized in the absence of calcium and even further destabilized at acidic versus neutral pH. The temperature and calcium concentration dependence of LR5 stability clearly indicate that under endosomal conditions the unfolded conformation of the repeat is largely dominant. We thus propose a new mechanism for LDL release in the endosome in which calcium depletion and decreased stability at acidic pH drives LR5 unfolding, which triggers LDL release from the receptor. Subsequent binding of LR5 to the EGF precursor domain, if it takes place at low calcium concentrations, would contribute to a further shifting of the equilibrium toward dissociation.

Scrambled isomers as key intermediates in the oxidative folding of ligand binding module 5 of the low density lipoprotein receptor

The ligand binding module five (LA5) of the low density lipoprotein receptor is a small, single-domain protein of 40 residues and three disulfide bonds with a calcium binding motif that is essential for its structure and function. Several mutations in LA5 have been reported to cause familial hypercholesterolemia by impairing a proper folding of the module. The current study reports the oxidative folding and reductive unfolding pathways of wild type and mutant LA5 modules through kinetic and structural analysis of the trapped intermediates. Wild type LA5 folding involves an initial phase of nonspecific packing where the sequential oxidation of its cysteines gives rise to complex equilibrated populations of intermediates. In the presence of calcium, the attainment of a coordination-competent conformation becomes the rate-limiting step of folding while binding of the ion funnels both thermodynamically and kinetically the folding reaction toward the native state. In the absence of calcium, a scrambled isomer (termed Xa) constitutes the global free energy minimum of the folding process. Xa and the native form are stable, inter-convertible species whose relative populations at equilibrium appear displaced in disease-linked mutants toward the scrambled form. Because stable scrambled isomers such as Xa avoid the exposition of reactive cysteines in misfolded modules, they might constitute a strategy to prevent wrong interactions with other domains during folding of the receptor. Comparison of the folding pathways of wild type and mutant LA5 provides the molecular basis to understand how LA modules fold into a functional conformation or upon mutation misfold and lead to disease.