Molecular dynamics studies on mutants of Cu,Zn superoxide dismutase: The functional role of charged residues in the electrostatic loop VII

Molecular dynamics (MD) calculations have been performed on mutants of superoxide dismutase (SOD) on some residues present in the electrostatic loop. These calculations have provided the solution structures for the mutants Thr-137 → IIe and Arg; Lys-136 → Ala; Glu-132 → Gln; Glu-133 → Gln; Glu-132, Glu-133 → Gln-132, Gln-133 and → Gln-132, Lys-133. The structural and dynamic properties of these mutants have been correlated with the catalytic properties and available spectroscopic data. The water molecule present in the active site close to the copper ion in wild type (WT) SOD is missing in the MD average structure of the Thr-137 → IIe mutant, while this molecule is present in the MD average structures of all the other mutants and of WT SOD. This agrees with the experimental data. This is an important result that shows the validity of our calculations and their ability to reproduce even subtle structural features. Addition of one or more positive charges on the 132 and/or 133 positions does not sizably perturb the structure of the active site channel, while the introduction of a positively charged residue (Arg) on position 137 has a large effect on the structure of the electrostatic loop. Analysis of the MD average structures of these mutants has pointed out that the simple electrostatic effects of charged residues in the channel are not the only factor relevant for enzymatic behavior but that the structure of the electrostatic loop and the location of the charged residues also contribute to the catalytic properties of SOD. © 1994 John Wiley & Sons, Inc.     

Carbon monoxide binding to the heme group at the dimeric interface modulates structure and copper accessibility in the Cu,Zn superoxide dismutase from Haemophilus ducreyi: in silico and in vitro evidences

X-ray absorption near-edge structure (XANES) spectroscopy and molecular dynamics (MD) simulations have been jointly applied to the study of the Cu,Zn superoxide dismutase from Haemophilus ducreyi (HdSOD) in interaction with the carbon monoxide molecule. The configurational flexibility of the Fe(II)-heme group, intercalated between the two subunits, has been sampled by MD simulations and included in the XANES data analysis without optimization in the structural parameter space. Our results provide an interpretation of the observed discrepancy in the Fe-heme distances as detected by extended X-ray absorption fine structure (EXAFS) spectroscopy and the classical XANES analysis, in which the structural parameters are optimized in a unique structure. Moreover, binding of the CO molecule to the heme induces a long range effect on the Cu,Zn active site, as evidenced by both MD simulations and in vitro experiments. MD simulation of the CO bound system, in fact, highlighted a structural rearrangement of the protein-protein hydrogen bond network in the region of the Cu,Zn active site, correlated with an increase in water accessibility at short distance from the copper atom. In line, in vitro experiments evidenced an increase of copper accessibility to a chelating agent when the CO molecule binds to the heme group, as compared to a heme deprived HdSOD. Altogether, our results support the hypothesis that the HdSOD is a heme-sensor protein, in which binding to small gaseous molecules modulates the enzyme superoxide activity as an adaptive response to the bacterial environment.     

Simulation of the diffusion-controlled reaction between superoxide and superoxide dismutase. I. Simple models

A Brownian dynamics simulation method is used to study the diffusion-influenced bimolecular reaction between superoxide and superoxide dismutase (SOD). Using simple models, the details of which are based on the crystallographic structure of SOD, it is found that the electrostatic charge distribution of SOD serves to guide superoxide into the active site and enhance the diffusion-controlled rate constant by about 40%.     

Atomic crystal and molecular dynamics simulation structures of human carbonic anhydrase II: insights into the proton transfer mechanism

Human carbonic anhydrase II (HCA II) is a zinc-metalloenzyme that catalyzes the reversible interconversion of CO2 and HCO3-. The rate-limiting step of this catalysis is the transfer of a proton between the Zn-bound solvent molecule and residue His64. In order to fully characterize the active site structural features implicated in the proton transfer mechanism, the refined X-ray crystal structure of uncomplexed wild type HCA II to 1.05 A resolution with an Rcryst value of 12.0% and an Rfree value of 15.1% has been elucidated. This structure provides strong clues as to the pathway of the intramolecular proton transfer between the Zn-bound solvent and His64. The structure emphasizes the role of the solvent network, the unique positioning of solvent molecule W2, and the significance of the dual conformation of His64 in the active site. The structure is compared with molecular dynamics (MD) simulation calculations of the Zn-bound hydroxyl/His64+ (charged) and the Zn-bound water/His64 (uncharged) HCA II states. A comparison of the crystallographic anisotropic atomic thermal parameters and MD simulation root-mean-square fluctuation values show excellent agreement in the atomic motion observed between the two methods. It is also interesting that the observed active site solvent positions in the crystal structure are also the most probable positions of the solvent during the MD simulations. On the basis of the comparative study of the MD simulation results, the HCA II crystal structure observed is most likely in the Zn-bound water/His64 state. This conclusion is based on the following observations: His64 is mainly (80%) orientated in an inward conformation; electron density omit maps infer that His64 is not charged in an either inward or outward conformation; and the Zn-bound solvent is most likely a water molecule.     

Cluster analysis of hydration waters around the active sites of bacterial alanine racemase using a 2-ns MD simulation

Structural data produced by a 2-ns molecular dynamics (MD) simulation on Geobacillus alanine racemase (AlaR; PDB: 1SFT) was used to study hydration around the two AlaR active sites. AlaR is a crucial enzyme for bacterial cell wall biosynthesis. It has been shown previously that the potency of an inhibitor can be increased by incorporating a functional group or atom that displaces hydration sites close to the substrate binding pocket of its target enzyme. The complete linkage algorithm was used for cluster analysis of the active site water positions from 126 solvent configurations sampled at regular intervals from the 2-ns MD simulation. Crystal waters in the 1SFT X-ray structure occupy most of the tightly bound water sites that were discovered. We show here that tightly bound water sites can be identified by cluster analysis of MD-generated coordinates starting with data supplied by a single X-ray structure, and we predict a highly conserved hydration site close to the carboxyl oxygen of L-Ala substrate. This approach holds promise for accelerating the drug design process. We also discuss an analysis of the well-known notion of residence time and introduce a new measure called retention time.     

Structural Evidence for a Copper-Bound Carbonate Intermediate in the Peroxidase and Dismutase Activities of Superoxide Dismutase

Copper-zinc superoxide dismutase (SOD) is of fundamental importance to our understanding of oxidative damage. Its primary function is catalysing the dismutation of superoxide to O₂ and H₂O₂. SOD also reacts with H₂O₂, leading to the formation of a strong copper-bound oxidant species that can either inactivate the enzyme or oxidise other substrates. In the presence of bicarbonate (or CO₂) and H₂O₂, this peroxidase activity is enhanced and produces the carbonate radical. This freely diffusible reactive oxygen species is proposed as the agent for oxidation of large substrates that are too bulky to enter the active site. Here, we provide direct structural evidence, from a 2.15 Å resolution crystal structure, of (bi)carbonate captured at the active site of reduced SOD, consistent with the view that a bound carbonate intermediate could be formed, producing a diffusible carbonate radical upon reoxidation of copper. The bound carbonate blocks direct access of substrates to Cu(I), suggesting that an adjunct to the accepted mechanism of SOD catalysed dismutation of superoxide operates, with Cu(I) oxidation by superoxide being driven via a proton-coupled electron transfer mechanism involving the bound carbonate rather than the solvent. Carbonate is captured in a different site when SOD is oxidised, being located in the active site channel adjacent to the catalytically important Arg143. This is the probable route of diffusion from the active site following reoxidation of the copper. In this position, the carbonate is poised for re-entry into the active site and binding to the reduced copper.     

Probing the Structural Basis for Enzyme-Substrate Recognition In Cu,Zn Superoxide Dismutase

A full understanding of enzyme-substrate interactions requires a detailed knowledge of their structural basis at atomic resolution. Crystallographic and biochemical data have been analyzed with coupled computational and computer graphic approaches to characterize the molecular basis for recognition of the superoxide anion substrate by Cu. Zn superoxide dismutase (SOD). Detailed analysis of the bovine SOD structure aligned with SOD sequences from 15 species provides new results concerning the significance and molecular basis for sequence conservation. Specific roles have been assigned for all 23 invariant residues and additional residues exhibiting functional equivalence. Sequence invariance is dominated by 15 residues that form the active site stcreochemistry. supporting a primary biological function of superoxide dismutation. Using data from crystallographic structures and site-directed mutants, we are testing the role of individual residues in the active site channel, including (in human SOD) Glu132, Glu133, Lys136, Thr137, and Arg 143. Electrostatic calculations incorporating molecular flexibility suggest that the region of positive electrostatic potential in and over the active site channel above the Cu ion sweeps through space during molecular motion to enhance the facilitated diffusion responsible for the enzyme's rapid catalytic rate.     

Structural Analysis of Peroxide-Soaked MnSOD Crystals Reveals Side-On Binding of Peroxide to Active-Site Manganese

The superoxide dismutase (SOD) enzymes are important antioxidant agents that protect cells from reactive oxygen species. The SOD family is responsible for catalyzing the disproportionation of superoxide radical to oxygen and hydrogen peroxide. Manganese- and iron-containing SOD exhibit product inhibition whereas Cu/ZnSOD does not. Here, we report the crystal structure of Escherichia coli MnSOD with hydrogen peroxide cryotrapped in the active site. Crystallographic refinement to 1.55 Å and close inspection revealed electron density for hydrogen peroxide in three of the four active sites in the asymmetric unit. The hydrogen peroxide molecules are in the position opposite His26 that is normally assumed by water in the trigonal bipyramidal resting state of the enzyme. Hydrogen peroxide is present in active sites B, C, and D and is side-on coordinated to the active-site manganese. In chains B and D, the peroxide is oriented in the plane formed by manganese and ligands Asp167 and His26. In chain C, the peroxide is bound, making a 70° angle to the plane. Comparison of the peroxide-bound active site with the hydroxide-bound octahedral form shows a shifting of residue Tyr34 towards the active site when peroxide is bound. Comparison with peroxide-soaked Cu/ZnSOD indicates end-on binding of peroxide when the SOD does not exhibit inhibition by peroxide and side-on binding of peroxide in the product-inhibited state of MnSOD.     

Effect of local environment and protein on the mechanism of action of superoxide dismutase

Quantum mechanical simulations of the mechanism of action of superoxide dismutase (SOD) indicate that the presence of Arg-141 in the active site of the enzyme is responsible for the formation of an intermediate complex between superoxide and the enzyme in which the copper is not reduced. The analysis of the local environmental effects of Arg-141 shows that this residue prevents the reduction of copper by forming a hydrogen bond to superoxide and by generating an electric field in the active site that opposes the transfer of an electron from superoxide to copper. The protein enhances the effect of the opposing field generated by Arg-141. Local changes in the environment of the copper ion, simulated by stretching the Cu-NE2 (His-61) bond also do not induce an electron transfer from superoxide to copper. The protein increases the energies required for this stretch through the electric field it generates near the active site. These results are further support for the new proposed mechanism of action of SOD which is based on the inability of superoxide to reduce the cupric ion in the enzyme.     

Catalytic and structural effects of amino acid substitution at histidine 30 in human manganese superoxide dismutase: insertion of valine C gamma into the substrate access channel

Catalysis of the disproportionation of superoxide by human manganese superoxide dismutase (MnSOD) is characterized by an initial burst of catalysis followed by a much slower region that is zero order in superoxide and due to a product inhibition by peroxide anion. We have prepared site-specific mutants with replacements at His30, the side chain of which lies along the substrate access channel and is about 5.8 A from the metal. Using pulse radiolysis to generate superoxide, we have determined that kcat/K(m) was decreased and product inhibition increased for H30V MnSOD, both by 1-2 orders of magnitude, compared with wild type, H30N, and H30Q MnSOD. These effects are not attributed to the redox potentials, which are similar for all of these variants. An investigation of the crystal structure of H30V Mn(III)SOD compared with wild type, H30Q, and H30N Mn(III)SOD showed the positions of two gamma carbons of Val30 in the active site; Cgamma1 overlaps Cgamma of His30 in wild type, and Cgamma2 extends into the substrate access channel and occupies the approximate position of a water molecule in the wild type. The data suggest that Cgamma2 of the Val side chain has significantly interrupted catalysis by this overlap into the access channel with possible overlap with the substrate-product binding site. This is supported by comparison of the crystal structure of H30V MnSOD with that of azide bound to Mn(III)SOD from Thermus thermophilus and by visible absorption spectra showing that azide binding to the metal in H30V Mn(III)SOD is abolished. Moreover, the presence of Val30 caused a 100-fold decrease in the rate constant for dissociation of the product-inhibited complex compared with wild type.     

Computational study of IAG-nucleoside hydrolase: determination of the preferred ground state conformation and the role of active site residues

The mechanism of action of inosine-adenosine-guanosine nucleoside hydrolase (IAG-NH) has been investigated by long-term molecular dynamics (MD) simulation in TIP3P water using stochastic boundary conditions. Special attention has been given to the role of leaving group pocket residues and conformation of the bound substrate at the active site of IAG-NH. We also describe the positioning of the residues of an important flexible loop at the active site, which was previously unobservable by X-ray crystallography due to high B-factors. Five MD simulations have been performed with the Enzyme x Substrate complexes: Enzyme x anti-Adenosine with Asp40-COOH [E(40H) x Ade(a)], Enzyme x anti-Adenosine with Asp40-COO- [E(40) x Ade(a)], Enzyme x syn-Adenosine with Asp40-COOH [E(40H) x Ade(s)], Enzyme x syn-Adenosine with Asp40-COO- [E(40) x Ade(s)], and Enzyme x anti-Inosine with Asp40-COO- [E(40) x Ino(a)]. Overall, the structures generated from the MD simulation of E(40H) x Ade(s) preserve the catalytically important hydrogen bonds as well as electrostatic and hydrophobic interactions to provide a plausible catalytic structure. When deprotonated Asp40 (Asp4-COO-) is present, the active site is open to water solvent which interferes with the base stacking between Trp83 and nucleobase. A calculation using Poisson-Boltzmann equation module supports that Asp40 indeed has an elevated pK(app). Solvent accessible surface area (SASA) calculations on all the five MD structures shows that systems with protonated Asp40, namely, E(40H) x Ade(a) and E(40H) x Ade(s), have zero SASA. It is found that a water molecule is hydrogen-bonded to the N7 of the nucleobase and is probably the essential general acid to protonate the departing nucleobase anion. The N7-bonded water is in turn hydrogen-bonded to waters in a channel, held in place by the residues of the flexible loop, Tyr257, His247, and Cys245. Using normal-mode analysis with elastic network model, we find that the flexible loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.     

Alterations in local stability and dynamics of A4V SOD1 in the presence of trifluoroethanol

Alterations in the local dynamics of Cu/Zn Superoxide dismutase (SOD1) due to mutations affect the protein folding, stability, and function leading to misfolding and aggregation seen in amyotrophic lateral sclerosis (ALS). Here, we study the structure and dynamics of the most devastating ALS mutation, A4V SOD1 in aqueous trifluoroethanol (TFE) through experiments and simulation. Far‐UV circular dichroism (CD) studies shows that TFE at intermediate concentrations (～15% ‐ 30%) induce partially unfolded β‐sheet‐rich extended conformations in A4V SOD1 which subsequently aggregates. Molecular dynamics (MD) simulation results shows that A4V SOD1 increases local dynamics in the active site loops that leads to the destabilization of the β‐barrel and loss of hydrophobic contacts, thus stipulating a basis for aggregation. Free energy landscape (FEL) and essential dynamics (ED) analysis demonstrates the conformational heterogeneity in A4V SOD1. Our results thus shed light on the role of local unfolding and conformational dynamics in aggregation of SOD1.     

An alternative mechanism of bicarbonate-mediated peroxidation by copper-zinc superoxide dismutase: rates enhanced via proposed enzyme-associated peroxycarbonate intermediate

Hydrogen peroxide can interact with the active site of copper-zinc superoxide dismutase (SOD1) to generate a powerful oxidant. This oxidant can either damage amino acid residues at the active site, inactivating the enzyme (the self-oxidative pathway), or oxidize substrates exogenous to the active site, preventing inactivation (the external oxidative pathway). It is well established that the presence of bicarbonate anion dramatically enhances the rate of oxidation of exogenous substrates. Here, we show that bicarbonate also substantially enhances the rate of self-inactivation of human wild type SOD1. Together, these observations suggest that the strong oxidant formed by hydrogen peroxide and SOD1 in the presence of bicarbonate arises from a pathway mechanistically distinct from that producing the oxidant in its absence. Self-inactivation rates are further enhanced in a mutant SOD1 protein (L38V) linked to the fatal neurodegenerative disorder, familial amyotrophic lateral sclerosis. The 1.4 A resolution crystal structure of pathogenic SOD1 mutant D125H reveals the mode of oxyanion binding in the active site channel and implies that phosphate anion attenuates the bicarbonate effect by competing for binding to this site. The orientation of the enzyme-associated oxyanion suggests that both the self-oxidative and external oxidative pathways can proceed through an enzyme-associated peroxycarbonate intermediate.     

Molecular dynamic simulations of the N-terminal receiver domain of NtrC reveal intrinsic conformational flexibility in the inactive state

The N-terminal receiver domain of NtrC is the molecular switch in the two-component signal transduction. It is the first protein where structures of both the active (phosphyroylated) and inactive (unphosphyroylated) states are determined experimentally. Phosphorylation of the NtrC at the active site induces large structural change. NMR experiments suggested that the wild type unphosphorylated NtrC adopts both the active and the inactive conformations and the phosphorylation stabilizes the active conformations. We applied free (unconstrained) molecular dynamic (MD) simulation to examine the intrinsic flexibilities and stabilities of the NtrC receiver domain in both the active and inactive conformations. Molecular dynamic simulations showed that the inactive state of NtrC receiver domain is more flexible than the active state. There were large movements in helix 4 and loop beta3-alpha3 which coincide with major structural differences between the inactive and active states. We observed large root-mean-square deviations from the initial starting structure and the large root-mean-square fluctuations during MD simulation for the inactive state. We then investigated the activation pathway with Targeted MD simulation. We show that the intrinsic flexibility in the loop beta3-alpha3 plays an important role in triggering the conformational change. Phosphorylation at the active site may serve to stabilize the conformational change. These results together suggest that the unphosphorylated NtrC receiver domain could be involved in a conformational equilibrium between two different states.     

Pathways of H2 toward the active site of [NiFe]-hydrogenase

Hydrogenases catalyze the reversible oxidation of molecular hydrogen (H(2)), but little is known about the diffusion of H(2) toward the active site. Here we analyze pathways for H(2) permeation using molecular dynamics (MD) simulations in explicit solvent. Various MD simulation replicates were done, to improve the sampling of the system states. H(2) easily permeates hydrogenase in every simulation and it moves preferentially in channels. All H(2) molecules that reach the active site made their approach from the side of the Ni ion. H(2) is able to reach distances of <4 A from the active site, although after 6 A permeation is difficult. In this region we mutated Val-67 into alanine and perform new MD simulations. These simulations show an increase of H(2) inside the protein and at lower distances from the active site. This valine can be a control point in the H(2) access to the active center.     

Is the activity-linked electrostatic gradient of bovine Cu, Zn superoxide dismutases conserved in homologous enzymes irrespective of the number and distribution of charges?

Electrostatic potential calculations have been performed on three different Cu, Zn superoxide dismutases (superoxide: superoxide oxidoreductase, EC 1.15 1.1 ), in order to evaluate the degree of conservation of the pattern of the electrostatic interactions between O₂⁻ and the active site recently pointed out in bovine Cu Zn SOD. The three Cu, Zn SODs that have been selected for this study, namely the bovine, ovine, and porcine enzymes, are highly homologous as to reasonably assume identical three-dimensional structure but display large differences in their net charge, as shown by their pI's which span over a wide range pHrange: 8.0 (sheep), 6.5(pig), 5.2(ox). Despite such a large difference in the net protein charge and in the spatial arrangement of electrostatic charges, electrostatic potential calculations show that the electrostatic channel directing the negatively charged substrate toward the positive catalytic site is strictly preserved with the same features for the three proteins. This suggests that the electrostatic funnel for conducting small anions into the active site is a highly conservative property in the evolution of Cu, Zn SOD.     

Crystal structure of a nucleoside diphosphate kinase from Bacillus halodenitrificans: coexpression of its activity with a Mn-superoxide dismutase

We found that when grown under anaerobic conditions the moderate halophile, gram-positive bacterium Bacillus halodenitrificans (ATCC 49067) synthesizes large amounts of a polypeptide complex that contains a heme center capable of reversibly bind nitric oxide. This complex, when exposed to air, dissociates and reassociates into two active components, a Mn-containing superoxide dismutase (SOD) and a nucleoside diphosphate kinase (BhNDK). The crystal structure of this latter enzyme has been determined at 2.2A resolution using molecular replacement method, based on the crystal structure of Drosophila melanogaster NDK. The model contains 149 residues of a total 150 residues and 34 water molecules. BhNDK consists of a four-stranded antiparallel beta-sheet, whose surfaces are partially covered by six alpha-helices, and its overall and active site structures are similar to those of homologous enzymes. However, the hexameric packing of BhNDK shows that this enzyme is different from both eukaryotic and gram-negative bacteria. The need for the bacterium to presynthesize both SOD and NDK precursors which are activated during the anaerobic-aerobic transition is discussed.     

Molecular Dynamic Simulations of the N-Terminal Receiver Domain of NtrC Reveal Intrinsic Conformational Flexibility in the Inactive State

Abstract

The N-terminal receiver domain of NtrC is the molecular switch in the two-component signal transduction. It is the first protein where structures of both the active (phosphyroylated) and inactive (unphosphyroylated) states are determined experimentally. Phosphorylation of the NtrC at the active site induces large structural change. NMR experiments suggested that the wild type unphosphorylated NtrC adopts both the active and the inactive conformations and the phosphorylation stabilizes the active conformations. We applied free (unconstrained) molecular dynamic (MD) simulation to examine the intrinsic flexibilities and stabilities of the NtrC receiver domain in both the active and inactive conformations. Molecular dynamic simulations showed that the inactive state of NtrC receiver domain is more flexible than the active state. There were large movements in helix 4 and loop β3-α3 which coincide with major structural differences between the inactive and active states. We observed large root-mean-square deviations from the initial starting structure and the large root-mean-square fluctuations during MD simulation for the inactive state. We then investigated the activation pathway with Targeted MD simulation. We show that the intrinsic flexibility in the loop β3-α3 plays an important role in triggering the conformational change. Phosphorylation at the active site may serve to stabilize the conformational change. These results together suggest that the unphosphorylated NtrC receiver domain could be involved in a conformational equilibrium between two different states.