Structural alterations by five disease-causing mutations in the low-pH conformation of human dihydrolipoamide dehydrogenase (hLADH) analyzed by molecular dynamics – Implications in functional loss and modulation of reactive oxygen species generation by pathogenic hLADH forms

Human dihydrolipoamide dehydrogenase (hLADH) is a flavoenzyme component (E3) of the human alpha-ketoglutarate dehydrogenase complex (α-KGDHc) and few other dehydrogenase complexes. Pathogenic mutations of hLADH cause severe metabolic diseases (atypical forms of E3 deficiency) that often escalate to cardiological or neurological presentations and even premature death; the pathologies are generally accompanied by lactic acidosis. hLADH presents a distinct conformation under acidosis (pH 5.5–6.8) with lower physiological activity and the capacity of generating reactive oxygen species (ROS). It has been shown by our laboratory that selected pathogenic mutations, besides lowering the physiological activity of hLADH, significantly stimulate ROS generation by hLADH, especially at lower pH, which might play a role in the pathogenesis of E3-deficiency in respective cases. Previously, we generated by molecular dynamics (MD) simulation the low-pH hLADH structure and analyzed the structural changes induced in this structure by eight of the pathogenic mutations of hLADH. In the absence of high resolution mutant structures these pieces of information are crucial for the mechanistic investigation of the molecular pathogeneses of the hLADH protein. In the present work we analyzed by molecular dynamics simulation the structural changes induced in the low-pH conformation of hLADH by five pathogenic mutations of hLADH; the structures of these disease-causing mutants of hLADH have never been examined before.     

Inhibition of the alpha-ketoglutarate dehydrogenase-mediated reactive oxygen species generation by lipoic acid

Dihydrolipoamide dehydrogenase (LADH) is a flavo-enzyme that serves as a subunit of α-ketoglutarate dehydrogenase complex (α-KGDHC). Reactive oxygen species (ROS) generation by α-KGDHC has been assigned to LADH (E3 subunit) and explained by the diaphorase activity of E3. Dysfunctions of α-KGDHC and concurrent ROS production have been implicated in neurodegeneration, ischemia-reperfusion, and other pathological conditions. In this work we investigated the in-depth details of ROS generation by isolated LADH and α-KGDHC. We found a parallel generation of superoxide and hydrogen peroxide by the E3 subunit of α-KGDHC which could be blocked by lipoic acid (LA) acting on a site upstream of the E3 subunit. The pathologically relevant ROS generation (at high NADH/NAD⁺ ratio and low pH) in the reverse mode of α-KGDHC could also be inhibited by LA. Our results contradict the previously proposed mechanism for pH-dependent ROS generation by LADH, showing no disassembling of the E3 functional homodimer at acidic pH using a physiologically relevant method for the examination. It is also suggested that LA could be beneficial in reducing the cell damage related to excessive ROS generation under pathological conditions.     

Periplasmic cold expression and one-step purification of human dihydrolipoamide dehydrogenase

Dihydrolipoamide dehydrogenase (LADH) is a FAD-linked subunit of alpha-ketoglutarate, pyruvate and branched-chain amino acid dehydrogenases and the glycine cleavage system. As an oxidoreductase it transfers electrons from the dihydrolipoic acid prosthetic group to the NAD(+) cofactor via its FAD center. Besides its physiological function it is capable of generating harmful reactive oxygen species (ROS) in pathological settings therefore it is implicated in neurodegeneration, ischemia-reperfusion, cancer and several other disorders. Pathological mutants of the enzyme cause severe, sometimes lethal syndromes like hypotonia, metabolic acidosis or inefficiency in development. Recently it has been revealed that LADH is a moonlighting protease when specific mutations in the dimerization surface destabilize the functional homodimer and expose a serine-protease-like catalytic dyad. As the basis of versatile functions of LADH is far from elucidation, there is a constant need for a pure and functional enzyme product for investigations. Several studies used recombinant human LADH before, however, it was generated by more complicated and/or physiologically less compatible protocols than reported here; most papers on functional and structural studies do not even report detailed protocols and characteristics (most importantly the purity) of their protein products. Here we describe the details of an optimized, easy-to-use periplasmic expression and one-step purification protocol for obtaining a highly pure, active and authentic (tag-cleaved) enzyme with the characterization of the protein product. The purified LADH can be used in biophysical and structural studies while the published protocol is easily convertible to a protein labeling procedure.     

Evaluation of Influence of Single Nucleotide Polymorphisms in Cytochrome P450 2B6 on Substrate Recognition Using Computational Docking and Molecular Dynamics Simulation

In this study, we investigated the influence of single nucleotide polymorphisms on the conformation of mutated cytochrome P450 (CYP) 2B6 proteins using molecular dynamics (MD) simulation. Some of these mutations influence drug metabolism activities, leading to individual variations in drug efficacy and pharmacokinetics. Using computational docking, we predicted the structure of the complex between the antimalarial agent artemether and CYP2B6 whose conformations were obtained by MD simulation. The simulation demonstrated that the entire structure of the protein changes even when a single residue is mutated. Moreover, the structural flexibility is affected by the mutations and it may influence the enzyme activity. The results suggest that some of the inactive mutants cannot recognize artemether due to structural changes caused by the mutation.     

The effect of aqueous surfactant solutions on alcohol dehydrogenase (LADH)

Alcohol dehydrogenase (LADH) was studied in aqueous solutions of surfactants to determine its structural and catalytic characteristics. Fluorescence, circular dichroism (CD), and electron paramagnetic resonance (ERP) techniques were used to study structural changes to the enzyme. The activity of LADH in catalyzing the oxidation of ethanol was investigated. Short-chain alkyl sulfonates and sulfates did not deactivate LADH or alter its structure. Longer and branched alkyl sulfates and sulfonates, as well as a cationic surfactant (CTAB), affected both LADH activity and conformation. (c) 1993 John Wiley & Sons, Inc.     

Dissociation of outer-sphere water is rate-limiting for the binding of ligands in the active site of horse liver alcohol dehydrogenase

The association of imidazole and auramine O to native horse-liver alcohol dehydrogenase [Zn(II)LADH] and active-site specifically cobalt(II)-substituted horse-liver alcohol dehydrogenase [Co(II)LADH], respectively, has been investigated. In all cases [except imidazole binding to Zn(II)LADH in the presence of auramine O] the association rates approached an upper limit (kmax). The kmax values were compared for the metal ligands imidazole (monodentate), 1,10-phenanthroline and 2,2'-bipyridine (bidentate; see also the preceding paper), and for auramine O which does not coordinate to the catalytic metal ion. Independent of the large differences in their structure and metal-bonding capability, all these compounds exhibit common, maximum, limiting rate constants of about 60 s-1 and 200 s-1 for Co(II)LADH and Zn(II)LADH, respectively. These results demonstrate that kmax is strongly dependent on the catalytic metal ion but not on the ligand. The absence of spectral changes in the d-d transitions of the catalytic Co(II) ion upon auramine O binding to Co(II)LADH indicates that the rate-limiting step is not accompanied by a major conformational change. Finally, it is concluded that reactions in the inner coordination sphere of the catalytic metal ion (i.e. the metal-bound water molecule) are not responsible for the step characterized by kmax. We propose the rate-limiting step to consist of the dissociation of one or several water molecules from the second coordination sphere of the catalytic metal ion in the active site of LADH in its open conformation.     

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

The molecular basis of conformational diseases frequently resides in mutant proteins constituting a subset of the vast mutational space. While the subtleties of protein structure point to molecular dynamics (MD) techniques as promising tools for an efficient exploration of such a space, the average size of proteins and the time scale of unfolding events make this goal difficult with present computational capabilities. We show here, nevertheless, that an efficient approach is already feasible for modular proteins. Familial hypercholesterolemia (FH) is a conformational disease linked to mutations in the gene encoding the low density lipoprotein receptor. A high percentage of these mutations has been found in the seven small modular binding repeats of the receptor. Taking advantage of its small size, we have performed an in depth MD study of the fifth binding repeat. Fast unfolding dynamics have been observed in the absence of a structural bound calcium ion, which agrees with its reported essential role in the stability of the module. In addition, several mutations detected in FH patients have been analyzed, starting from the native conformation. Our results indicate that in contrast with the wild type protein and an innocuous control mutant, disease-related mutants experience, in short simulation times (2-8 ns), gross departures from the native state that lead to unfolded conformations and, in some cases, to binding site desorganization deriving in calcium release. Computational diagnosis of mutations leading to conformational diseases seems thus feasible, at least for small or modular pathogenic proteins.     

Rapid upregulation of endothelial P-selectin expression via reactive oxygen species generation

Endothelial cell ICAM-1 upregulation in response to TNF-alpha is mediated in part by reactive oxygen species (ROS) generated by the endothelial membrane-associated NADPH oxidase and occurs maximally after 4 h as the synthesis of new protein is required. However, thrombin-stimulated P-selectin upregulation is bimodal, the first peak occurring within minutes. We hypothesize that this early peak, which results from the release of preformed P-selectin from within Weibel-Palade bodies, is mediated in part by ROS generated from the endothelial membrane-associated xanthine oxidase. We found that this rapid expression of P-selectin on the surface of endothelial cells was accompanied by qualitatively parallel increases in ROS generation. Both P-selectin expression and ROS generation were inhibited, dose dependently, by the exogenous administration of disparate cell-permeable antioxidants and also by the inhibition of either of the known membrane-associated ROS-generating enzymes NADPH oxidase or xanthine oxidase. This rapid, posttranslational cell signaling response, mediated by ROS generated not only by the classical NADPH oxidase but also by xanthine oxidase, may well represent an important physiological trigger of the microvascular inflammatory response.     

Essential role of complex II of the respiratory chain in hypoxia-induced ROS generation in the pulmonary vasculature

In the pulmonary vasculature, the mechanisms responsible for oxygen sensing and the initiation of hypoxia-induced vasoconstriction and vascular remodeling are still unclear. Nitric oxide (NO) and reactive oxygen species (ROS) are discussed as early mediators of the hypoxic response. Here, we describe a quantitative analysis of NO- and ROS-producing cells within the vascular walls of murine lung sections cultured at normoxia or hypoxia. Whereas the number of NO-producing cells was not changed by hypoxia, the number of ROS-generating cells was significantly increased. Addition of specific inhibitors revealed that mitochondria were the source of ROS. The participation of the individual mitochondrial complexes differed in normoxic and hypoxic ROS generation. Whereas normoxic ROS production required complexes I and III, hypoxic ROS generation additionally demanded complex II. Histochemically demonstrable succinate dehydrogenase activity of complex II in the arterial wall decreased during hypoxia. Inhibition of the reversed enzymatic reaction, i.e., fumarate reductase, by application of succinate, specifically abolished hypoxic, but not normoxic, ROS generation. Thus complex II plays an essential role in hypoxic ROS production. Presumably, its catalytic activity switches from succinate dehydrogenase to fumarate reductase at reduced oxygen tension, thereby modulating the directionality of the electron flow.     

Molecular dynamic studies on the impact of mutations on the structure,stability, and N‐terminal orientation of annexin A1: Implications for membrane aggregation

Multiple MD simulations were performed for the full‐length wild‐type A1, the full length A1 mutations S27E and S27A, as well as the N‐terminal peptide (AMVSEFLKQAWFIDNEEQEYIKTVKG S ²⁷ KGGPGSAVSPYPTFN) of wild‐type A1 and mutations S27E and S27A. The MD simulation trajectories of about 350 ns were generated and analyzed to examine the changes of core domain calcium binding affinity, core domain and N‐terminal domain structures, and N‐terminal domain orientation. Our results indicated that S27A and S27E mutations caused little changes on the calcium‐binding affinity of the core domain of A1. However, the S27A mutation made the N‐terminal domain of A1 less helical, and made the N‐terminal domain migrate faster toward the core domain; these impacts on A1 are beneficial to the membrane aggregation process. On the contrary, the S27E mutation made the N‐terminal domain of A1 more stable, and hindered the migration to the core domain; these changes on A1 are antagonistic for the membrane aggregation process. Our results using MD simulations provide an atomistic explanation for experimental observations that the S27E mutant showed a higher calcium concentration requirement and lower maximal extent of aggregation, while the wild‐type and two mutants S27E and S27A required identical calcium concentrations for liposome binding. Proteins 2014; 82:3327–3334. © 2014 Wiley Periodicals, Inc.     

Copper increases the damage to DNA and proteins caused by reactive oxygen species

Copper [Cu(II)] is an ubiquitous transition and trace element in living organisms. It increases reactive oxygen species (ROS) and free-radical generation that might damage biomolecules like DNA, proteins, and lipids. Furthermore, ability of Cu(II) greatly increases in the presence of oxidants. ROS, like hydroxyl (·OH) and superoxide (·O₂) radicals, alter both the structure of the DNA double helix and the nitrogen bases, resulting in mutations like the AT→GC and GC→AT transitions. Proteins, on the other hand, suffer irreversible oxidations and loss in their biological role. Thus, the aim of this investigation is to characterize, in vitro, the structural effects caused by ROS and Cu(II) on bacteriophage λ DNA or proteins using either hydrogen peroxide (H₂O₂) or ascorbic acid with or without Cu(II). Exposure of DNA to ROS-generating mixtures results in electrophoretic (DNA breaks), spectrophotometric (band broadening, hypochromic, hyperchromic, and bathochromic effects), and calorimetric (denaturation temperature [T _d], denaturation enthalpy [ΔH], and heat capacity [C _p] values) changes. As for proteins, ROS increased their thermal stability. However, the extent of the observed changes in DNA and proteins were distinct, depending on the efficiency of the systems assayed to generate ROS. The resulting effects were most evident when Cu(II) was present. In summary, these results show that the ROS, ·O₂ and ·OH radicals, generated by the Cu(II) systems assayed deeply altered the chemical structure of both DNA and proteins. The physiological relevance of these structural effects should be further investigated.     

Proliferation and wound healing of vascular cells trigger the generation of extracellular reactive oxygen species and LDL oxidation

Cell proliferation of vascular cells is a key feature in vascular biology, wound healing, and pathophysiological processes such as atherosclerosis and restenosis. In atherosclerotic intima, cell proliferation colocalizes with oxidized LDL that indicate a local oxidative stress. This study aims to investigate whether cell proliferation is causally related with extracellular ROS generation and subsequent LDL oxidation. Sparse proliferating endothelial and smooth muscle cells generate higher levels of extracellular ROS (O₂^•− and H₂O₂) and LDL oxidation than confluent contact-inhibited cells. During wound healing of confluent cell layer, cell proliferation associated with healing also induced enhanced extracellular ROS generation and LDL oxidation. Proliferation-associated extracellular ROS generation is mediated through mitogenic signaling pathways, involving ERK1/2 and PKC, but is independent of de novo DNA synthesis, gene expression and protein synthesis. Data obtained with inhibitors of oxidases suggest that proliferation-associated extracellular ROS are not generated by a single ROS-generating system and are not essential for cell proliferation. In conclusion, our data show that proliferating vascular cells (in sparse culture or during wound healing) generate high levels of extracellular ROS and LDL oxidation through regulation of ROS-generating systems by mitogenic signaling. This constitutes a link between proliferative events and oxidative stress/LDL oxidation in atherosclerotic lesions and restenosis.     

High hydrophobic amino acid exposure is responsible of the neurotoxic effects induced by E200K or D202N disease-related mutations of the human prion protein

Mutations in prion protein are thought to be causative of inherited prion diseases favoring the spontaneous conversion of the normal prion protein into the scrapie-like pathological prion protein. We previously reported that, by controlled thermal denaturation, human prion protein fragment 90-231 acquires neurotoxic properties when transformed in a β-rich conformation, resembling the scrapie-like conformation. In this study we generated prion protein fragment 90-231 bearing mutations identified in familial prion diseases (D202N and E200K), to analyze their role in the induction of a neurotoxic conformation. Prion protein fragment 90-231(wild type) and the D202N mutant were not toxic in native conformation but induced cell death only after thermal denaturation. Conversely, prion protein fragment 90-231(E200K) was highly toxic in its native structure, suggesting that E200K mutation per se favors the acquisition of a peptide neurotoxic conformation. To identify the structural determinants of prion protein fragment 90-231 toxicity, we show that while the wild type peptide is structured in α-helix, hPrP90-231 E200K is spontaneously refolded in a β-structured conformer characterized by increased proteinase K resistance and propensity to generate fibrils. However, the most significant difference induced by E200K mutation in prion protein fragment 90-231 structure in native conformation we observed, was an increase in the exposure of hydrophobic amino-acids on protein surface that was detected in wild type and D202N proteins only after thermal denaturation. In conclusion, we propose that increased hydrophobicity is one of the main determinants of toxicity induced by different mutations in prion protein-derived peptides.     

Activity and flexibility of alcohol dehydrogenase in organic solvents

The oxidation of cinnamyl alcohol to cinnamaldehyde by horse liver alcohol dehydrogenase (LADH) was carried out in nearly anhydrous organic solvents and in solvents containing from 0.1 to 10% added water. In nearly anhydrous solvents containing less than 0.02% water, the oxidation rate increased as the water solubility in the solvent decreased, but the reaction did not require active LADH. Moreover, the highest activity in nearly anhydrous heptane was obtained by lyophilizing the enzyme from a solution of pH 2.0, even though LADH exhibits virtually no enzymatic activity in water at this pH. The catalytic activity of LADH was restored and increased dramatically as small amounts of water were added to each solvent. In conjunction with the activity measurements, electron paramagnetic resonance (EPR) spectroscopy and two active-site directed spin labels were used to examine solvent-dependent structural features of LADH. The EPR spectra indicated that LADH became more rigid as the dielectric constant of the solvent decreased. The degree of rigidity also depended on the pH from which the enzyme was lyophilized, indicating that the ionization state of the enzyme can have an important influence on its dynamics in organic solvents. Finally, adding 1% water to organic solvents had no apparent effect on the enzyme's conformation or flexibility near the spin label, even though enzyme activity was an order of magnitude higher when 1% water was present.     

Raman study of reduced nicotinamide adenine dinucleotide bound to liver alcohol dehydrogenase

We report the first Raman spectra of reduced nicotinamide adenine dinucleotide (NADH) when bound to an enzymatic active site, that of liver alcohol dehydrogenase (LADH). This was obtained by subtracting the Raman spectrum of LADH from that of the binary LADH/NADH complex. There are significant changes in the spectrum of bound NADH as compared to that in solution. The data indicate that both the nicotinamide moiety and the adenine moiety are involved in the binding. At least one of the two NH2 moieties of NADH also participates.     

Molecular Dynamics Simulations and Principal Component Analysis on Human Laforin Mutation W32G and W32G/K87A

Mutations in human laforin lead to an autosomal neurodegenerative disorder Lafora disease. In N-terminal carbohydrate binding domain of laforin, two mutations W32G and K87A are reported as highly disease causing laforin mutants. Experimental studies reported that mutations are responsible for the abolishment of glycogen binding which is a critical function of laforin. Our current computational study focused on the role of conformational changes in human laforin structure due to existing single mutation W32G and prepared double mutation W32G/K87A related to loss of glycogen binding. We performed 10 ns molecular dynamics (MD) simulation studies in the Gromacs package for both mutations and analyzed the trajectories. From the results, the global properties like root mean square deviation, root mean square fluctuation, radius of gyration, solvent accessible surface area and hydrogen bonds showed structural changes in atomic level observed in W32G and W32G/K87A laforin mutants. The conformational change induced by mutants influenced the loss of the overall stability of the native laforin. Moreover, the change in overall motion of protein was analyzed by principal component analysis and results showed protein clusters expanded more than native and also change in direction in case of double mutant in conformational space. Overall, our report provides theoretical information on loss of structure–function relationship due to flexible nature of laforin mutants. In conclusion, comparative MD simulation studies support the experimental data on W32G and W32G/K87A related to the lafora disease mechanism on glycogen binding.     

Exploring structural and thermodynamic stabilities of human prion protein pathogenic mutants D202N, E211Q and Q217R

The central event in the pathogenesis of prion protein (PrP) is a profound conformational change from its α-helical (PrP(C)) to its β-sheet-rich isoform (PrP(Sc)). Many single amino acid mutations of PrP are associated with familial prion diseases, such as D202N, E211Q, and Q217R mutations located at the third native α-helix of human PrP. In order to explore the underlying structural and dynamic effects of these mutations, we performed all-atom molecular dynamics (MD) simulations for the wild-type (WT) PrP and its mutants. The obtained results indicate that these amino acid substitutions have subtle effects on the protein structures, but show large changes of the overall electrostatic potential distributions. We can infer that the changes of PrP electrostatic surface due to the studied mutations may influence the intermolecular interactions during the aggregation process. In addition, the mutations also affect the thermodynamic stabilities of PrP.     

Computational compensatory mutation discovery approach: Predicting a PARP1 variant rescue mutation

The prediction of protein mutations that affect function may be exploited for multiple uses. In the context of disease variants, the prediction of compensatory mutations that reestablish functional phenotypes could aid in the development of genetic therapies. In this work, we present an integrated approach that combines coevolutionary analysis and molecular dynamics (MD) simulations to discover functional compensatory mutations. This approach is employed to investigate possible rescue mutations of a poly(ADP-ribose) polymerase 1 (PARP1) variant, PARP1 V762A, associated with lung cancer and follicular lymphoma. MD simulations show PARP1 V762A exhibits noticeable changes in structural and dynamical behavior compared with wild-type (WT) PARP1. Our integrated approach predicts A755E as a possible compensatory mutation based on coevolutionary information, and molecular simulations indicate that the PARP1 A755E/V762A double mutant exhibits similar structural and dynamical behavior to WT PARP1. Our methodology can be broadly applied to a large number of systems where single-nucleotide polymorphisms have been identified as connected to disease and can shed light on the biophysical effects of such changes as well as provide a way to discover potential mutants that could restore WT-like functionality. This can, in turn, be further utilized in the design of molecular therapeutics that aim to mimic such compensatory effect.