Aberrantly increased hydrophobicity shared by mutants of Cu,Zn-superoxide dismutase in familial amyotrophic lateral sclerosis

More than 100 different mutations in the gene encoding Cu,Zn-superoxide dismutase (SOD1) cause preferential motor neuron degeneration in familial amyotrophic lateral sclerosis (ALS). Although the cellular target(s) of mutant SOD1 toxicity have not been precisely specified, evidence to date supports the hypothesis that ALS-related mutations may increase the burden of partially unfolded SOD1 species. Influences that may destabilize SOD1 in vivo include impaired metal ion binding, reduction of the intrasubunit disulfide bond, or oxidative modification. In this study, we observed that metal-deficient as-isolated SOD1 mutants (H46R, G85R, D124V, D125H, and S134N) with disordered electrostatic and zinc-binding loops exhibited aberrant binding to hydrophobic beads in the absence of other destabilizing agents. Other purified ALS-related mutants that can biologically incorporate nearly normal amounts of stabilizing zinc ions (A4V, L38V, G41S, D90A, and G93A) exhibited maximal hydrophobic behavior after exposure to both a disulfide reducing agent and a metal chelator, while normal SOD1 was more resistant to these agents. Moreover, we detected hydrophobic SOD1 species in lysates from affected tissues in G85R and G93A mutant but not wildtype SOD1 transgenic mice. These findings suggest that a susceptibility to the cellular disulfide reducing environment and zinc loss may convert otherwise stable SOD1 mutants into metal-deficient forms with locally destabilized electrostatic and zinc-binding loops. These abnormally hydrophobic SOD1 species may promote aberrant interactions of the enzyme with itself or with other cellular constituents to produce toxicity in familial ALS.     

Evidence for a dual role of an active site histidine in α-amino-β-carboxymuconate-ε-semialdehyde decarboxylase

The previously reported crystal structures of α-amino-β-carboxymuconate-ε-semialdehyde decarboxylase (ACMSD) show a five-coordinate Zn(II)(His)(3)(Asp)(OH(2)) active site. The water ligand is H-bonded to a conserved His228 residue adjacent to the metal center in ACMSD from Pseudomonas fluorescens (PfACMSD). Site-directed mutagenesis of His228 to tyrosine and glycine in this study results in a complete or significant loss of activity. Metal analysis shows that H228Y and H228G contain iron rather than zinc, indicating that this residue plays a role in the metal selectivity of the protein. As-isolated H228Y displays a blue color, which is not seen in wild-type ACMSD. Quinone staining and resonance Raman analyses indicate that the blue color originates from Fe(III)-tyrosinate ligand-to-metal charge transfer. Co(II)-substituted H228Y ACMSD is brown in color and exhibits an electron paramagnetic resonance spectrum showing a high-spin Co(II) center with a well-resolved (59)Co (I = 7/2) eight-line hyperfine splitting pattern. The X-ray crystal structures of as-isolated Fe-H228Y (2.8 ?) and Co-substituted (2.4 ?) and Zn-substituted H228Y (2.0 ? resolution) support the spectroscopic assignment of metal ligation of the Tyr228 residue. The crystal structure of Zn-H228G (2.6 ?) was also determined. These four structures show that the water ligand present in WT Zn-ACMSD is either missing (Fe-H228Y, Co-H228Y, and Zn-H228G) or disrupted (Zn-H228Y) in response to the His228 mutation. Together, these results highlight the importance of His228 for PfACMSD's metal specificity as well as maintaining a water molecule as a ligand of the metal center. His228 is thus proposed to play a role in activating the metal-bound water ligand for subsequent nucleophilic attack on the substrate.     

Structural characterization of zinc-deficient human superoxide dismutase and implications for ALS

Over 130 mutations to copper, zinc superoxide dismutase (SOD) are implicated in the selective death of motor neurons found in 25% of patients with familial amyotrophic lateral sclerosis (ALS). Despite their widespread distribution, ALS mutations appear positioned to cause structural and misfolding defects. Such defects decrease SOD's affinity for zinc, and loss of zinc from SOD is sufficient to induce apoptosis in motor neurons in vitro. To examine the importance of the zinc site in the structure and pathogenesis of human SOD, we determined the 2.0-A-resolution crystal structure of a designed zinc-deficient human SOD, in which two zinc-binding ligands have been mutated to hydrogen-bonding serine residues. This structure revealed a 9 degrees twist of the subunits, which opens the SOD dimer interface and represents the largest intersubunit rotational shift observed for a human SOD variant. Furthermore, the electrostatic loop and zinc-binding subloop were partly disordered, the catalytically important Arg143 was rotated away from the active site, and the normally rigid intramolecular Cys57-Cys146 disulfide bridge assumed two conformations. Together, these changes allow small molecules greater access to the catalytic copper, consistent with the observed increased redox activity of zinc-deficient SOD. Moreover, the dimer interface is weakened and the Cys57-Cys146 disulfide is more labile, as demonstrated by the increased aggregation of zinc-deficient SOD in the presence of a thiol reductant. However, equimolar Cu,Zn SOD rapidly forms heterodimers with zinc-deficient SOD (t1/2 approximately 15 min) and prevents aggregation. The stabilization of zinc-deficient SOD as a heterodimer with Cu,Zn SOD may contribute to the dominant inheritance of ALS mutations. These results have general implications for the importance of framework stability on normal metalloenzyme function and specific implications for the role of zinc ion in the fatal neuropathology associated with SOD mutations.     

Solution structure of Apo Cu,Zn superoxide dismutase: role of metal ions in protein folding

The solution structure of the demetalated copper, zinc superoxide dismutase is obtained for the monomeric Glu133Gln/Phe50Glu/Gly51Glu mutant through NMR spectroscopy. The demetalated protein still has a well-defined tertiary structure; however, two beta-strands containing two copper ligands (His46 and His48, beta4) and one zinc ligand (Asp83, beta5) are shortened, and the sheet formed by these strands and strands beta7 and beta8 moves away from the other strands of the beta-barrel to form an open clam with respect to a closed conformation in the holoprotein. Furthermore, loop IV which contains three zinc ligands (His63, His71, and His80) and loop VII which contributes to the definition of the active cavity channel are severely disordered, and experience extensive mobility as it results from thorough (15)N relaxation measurements. These structural and mobility data, if compared with those of the copper-depleted protein and holoprotein, point out the role of each metal ion in the protein folding, leading to the final tertiary structure of the holoprotein, and provide hints for the mechanisms of metal delivery by metal chaperones.     

Protein-protein interactions in actin-myosin binding and structural effects of R405Q mutation: a molecular dynamics study

Detailed residue-wise interactions involved in the binding of myosin to actin in the rigor conformation without nucleotides have been examined using molecular dynamics simulations of the chicken skeletal myosin head complexed with two actin monomers, based on the cryo-microscopic model of Holmes et al. (Nature 2003;425:423-427). The overall interaction is largely electrostatic in nature, because of the charged residues in the four loops surrounding the central primary binding site. The 50k/20k loop, disordered in crystal structures and in simulations of free myosin in solution, was found to be in a conformation stabilized with 1 - 2 internal salt bridges. The cardiomyopathy loop forms 2 - 3 interprotein salt bridges with actin monomers upon binding, whereas its Arg405 residue, the mutation site associated with the hypertrophic cardiomyopathy, forms a strong salt bridge with Glu605 in the neighboring helix away from actin in the actin-bound myosin. The myopathy loop of the R405Q mutant maintains a high degree of two-strand beta-sheet character when bound to actin with the corresponding salt bridges broken.     

Copper and zinc metallation status of copper-zinc superoxide dismutase from amyotrophic lateral sclerosis transgenic mice

Mutations in the metalloenzyme copper-zinc superoxide dismutase (SOD1) cause one form of familial amyotrophic lateral sclerosis (ALS), and metals are suspected to play a pivotal role in ALS pathology. To learn more about metals in ALS, we determined the metallation states of human wild-type or mutant (G37R, G93A, and H46R/H48Q) SOD1 proteins from SOD1-ALS transgenic mice spinal cords. SOD1 was gently extracted from spinal cord and separated into insoluble (aggregated) and soluble (supernatant) fractions, and then metallation states were determined by HPLC inductively coupled plasma MS. Insoluble SOD1-rich fractions were not enriched in copper and zinc. However, the soluble mutant and WT SOD1s were highly metallated except for the metal-binding-region mutant H46R/H48Q, which did not bind any copper. Due to the stability conferred by high metallation of G37R and G93A, it is unlikely that these soluble SOD1s are prone to aggregation in vivo, supporting the hypothesis that immature nascent SOD1 is the substrate for aggregation. We also investigated the effect of SOD1 overexpression and disease on metal homeostasis in spinal cord cross-sections of SOD1-ALS mice using synchrotron-based x-ray fluorescence microscopy. In each mouse genotype, except for the H46R/H48Q mouse, we found a redistribution of copper between gray and white matters correlated to areas of high SOD1. Interestingly, a disease-specific increase of zinc was observed in the white matter for all mutant SOD1 mice. Together these data provide a picture of copper and zinc in the cell as well as highlight the importance of these metals in understanding SOD1-ALS pathology.     

Crystal structures of 26 kDa Clonorchis sinensis glutathione S-transferase reveal zinc binding and putative metal binding

The crystal structures of CsGST in two different space groups revealed that Asp26 and His79 coordinate a zinc ion. In one space group, His46 of an adjacent molecule participates in the coordination within 2.0 Å. In the other space group, Asp26, His79 and a water molecule coordinate a zinc ion. The CsGST–D26H structure showed that four histidine residues – His26 and His79 from one molecule and the same residues from a symmetry-related neighboring molecule – coordinate a zinc ion. The coordinated zinc ions are located between two molecules and mediate molecular contacts within the crystal.     

Crystal structure of ribonuclease T1 complexed with adenosine 2'-monophosphate at 1.8-A resolution.

Ribonuclease T1 was purified from an Escherichia coli overproducing strain and co-crystallized with adenosine 2'-monophosphate (2'-AMP) by microdialysis against 50% (v/v) 2-methyl-2,4-pentanediol in 20 mM sodium acetate, 2 mM calcium acetate, pH 4.2. The crystals have orthorhombic space group P2(1)2(1)2(1), with cell dimensions a = 48.93(1), b = 46.57(4), c = 41.04(2) A; Z = 4 and V = 93520 A3. The crystal structure was determined on the basis of the isomorphous structure of uncomplexed RNase T1 (Martinez-Oyanedel et al. (1991) submitted for publication) and refined by least squares methods using stereochemical restraints. The refinement was based on Fhkl of 7,445 reflections with Fo greater than or equal to 1 sigma (Fo) in the resolution range of 10-1.8 A, and converged at a crystallographic R factor of 0.149. The phosphate group of 2'-AMP is tightly hydrogen-bonded to the side chains of the active site residues Tyr38, His40, Glu58, Arg77, and His92, comparable with vanadate binding in the respective complex (Kostrewa, D., Choe, H.-W., Heinemann, U., and Saenger, W. (1989) Biochemistry 28, 7592-7600) and different from the complex with guanosine 2'-monophosphate (Arni, R., Heinemann, U., Tokuoka, R., and Saenger, W. (1988) J. Biol. Chem. 263, 15358-15368) where the phosphate does not interact with Arg77 and His92. The adenosine moiety is not located in the guanosine recognition site but stacked on Gly74 carbonyl and His92 imidazole, which serve as a subsite, as shown previously (Lenz, A., Cordes, F., Heinemann, U., and Saenger, W. (1991) J. Biol. Chem. 266, 7661-7667); in addition, there are hydrogen bonds adenine N6H . . . O Gly74 (minor component of three-center hydrogen bond) and adenosine O5' . . . O delta Asn36. These binding interactions readily explain why RNase T1 has some affinity for 2'-AMP. The molecular structure of RNase T1 is only marginally affected by 2'-AMP binding. Its "empty" guanosine-binding site features a flipped Asn43-Asn44 peptide bond and the side chains of Tyr45, Glu46 adopt conformations typical for RNase T1 not involved in guanosine binding. The side chains of amino acids Leu26, Ser35, Asp49, Val78 are disordered. The disorder of Val78 is of interest since this amino acid is located in a hydrophobic cavity, and the disorder appears to be correlated with an "empty" guanosine-binding site. The two Asp15 carboxylate oxygens and six water molecules coordinate a Ca2+ ion 8-fold in the form of a square antiprism.     

Molecular structure of dihydroorotase: a paradigm for catalysis through the use of a binuclear metal center

Dihydroorotase plays a key role in pyrimidine biosynthesis by catalyzing the reversible interconversion of carbamoyl aspartate to dihydroorotate. Here we describe the three-dimensional structure of dihydroorotase from Escherichia coli determined and refined to 1.7 A resolution. Each subunit of the homodimeric enzyme folds into a "TIM" barrel motif with eight strands of parallel beta-sheet flanked on the outer surface by alpha-helices. Unexpectedly, each subunit contains a binuclear zinc center with the metal ions separated by approximately 3.6 A. Lys 102, which is carboxylated, serves as a bridging ligand between the two cations. The more buried or alpha-metal ion in subunit I is surrounded by His 16, His 18, Lys 102, Asp 250, and a solvent molecule (most likely a hydroxide ion) in a trigonal bipyramidal arrangement. The beta-metal ion, which is closer to the solvent, is tetrahedrally ligated by Lys 102, His 139, His 177, and the bridging hydroxide. L-Dihydroorotate is observed bound to subunit I, with its carbonyl oxygen, O4, lying 2.9 A from the beta-metal ion. Important interactions for positioning dihydroorotate into the active site include a salt bridge with the guanidinium group of Arg 20 and various additional electrostatic interactions with both protein backbone and side chain atoms. Strikingly, in subunit II, carbamoyl L-aspartate is observed binding near the binuclear metal center with its carboxylate side chain ligating the two metals and thus displacing the bridging hydroxide ion. From the three-dimensional structures of the enzyme-bound substrate and product, it has been possible to propose a unique catalytic mechanism for dihydroorotase. In the direction of dihydroorotate hydrolysis, the bridging hydroxide attacks the re-face of dihydroorotate with general base assistance by Asp 250. The carbonyl group is polarized for nucleophilic attack by the bridging hydroxide through a direct interaction with the beta-metal ion. During the cyclization of carbamoyl aspartate, Asp 250 initiates the reaction by abstracting a proton from N3 of the substrate. The side chain carboxylate of carbamoyl aspartate is polarized through a direct electrostatic interaction with the binuclear metal center. The ensuing tetrahedral intermediate collapses with C-O bond cleavage and expulsion of the hydroxide which then bridges the binuclear metal center.     

The crystal structure of the monomeric human SOD mutant F50E/G51E/E133Q at atomic resolution. The enzyme mechanism revisited.

The crystal structure of the engineered monomeric human Cu,ZnSOD triple mutant F50E/G51E/E133Q (Q133M2SOD) is reported at atomic resolution (1.02 A). This derivative has about 20 % of the wild-type activity. Crystals of Q133M2SOD have been obtained in the presence of CdCl2. The metal binding site is disordered, with both cadmium and copper ions simultaneously binding to the copper site. The cadmium (II) ions occupy about 45 % of the copper sites by binding the four histidine residues which ligate copper in the native enzyme, and two further water molecules to complete octahedral coordination. The copper ion is tri-coordinate, and the fourth histidine (His63) is detached from copper and bridges cadmium and zinc. X-ray absorption spectroscopy performed on the crystals suggests that the copper ion has undergone partial photoreduction upon exposure to the synchrotron light. The structure is also disordered in the disulfide bridge region of loop IV that is located at the subunit/subunit interface in the native SOD dimer. As a consequence, the catalytically relevant Arg143 residue is disordered. The present structure has been compared to other X-ray structures on various isoenzymes and to the solution structure of the same monomeric form. The structural results suggest that the low activity of monomeric SOD is due to the disorder in the conformation of the side-chain of Arg143 as well as of loop IV. It is proposed that the subunit-subunit interactions in the multimeric forms of the enzyme are needed to stabilize the correct geometry of the cavity and the optimal orientation of the charged residues in the active channel. Furthermore, the different coordination of cadmium and copper ions, contemporaneously present in the same site, are taken as models for the oxidized and reduced copper species, respectively. These properties of the structure have allowed us to revisit the enzymatic mechanism.     

Amyloid-like filaments and water-filled nanotubes formed by SOD1 mutant proteins linked to familial ALS

Mutations in the SOD1 gene cause the autosomal dominant, neurodegenerative disorder familial amyotrophic lateral sclerosis (FALS). In spinal cord neurons of human FALS patients and in transgenic mice expressing these mutant proteins, aggregates containing FALS SOD1 are observed. Accumulation of SOD1 aggregates is believed to interfere with axonal transport, protein degradation and anti-apoptotic functions of the neuronal cellular machinery. Here we show that metal-deficient, pathogenic SOD1 mutant proteins crystallize in three different crystal forms, all of which reveal higher-order assemblies of aligned beta-sheets. Amyloid-like filaments and water-filled nanotubes arise through extensive interactions between loop and beta-barrel elements of neighboring mutant SOD1 molecules. In all cases, non-native conformational changes permit a gain of interaction between dimers that leads to higher-order arrays. Normal beta-sheet-containing proteins avoid such self-association by preventing their edge strands from making intermolecular interactions. Loss of this protection through conformational rearrangement in the metal-deficient enzyme could be a toxic property common to mutants of SOD1 linked to FALS.     

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.     

Design, characterization, and structure of a biologically active single-chain mutant of human IFN-gamma

A mutant form of human interferon-gamma (IFN-gamma SC1) that binds one IFN-gamma receptor alpha chain (IFN-gamma R alpha) has been designed and characterized. IFN-gamma SC1 was derived by linking the two peptide chains of the IFN-gamma dimer by a seven-residue linker and changing His111 in the first chain to an aspartic acid residue. Isothermal titration calorimetry shows that IFN-gamma SC1 forms a 1:1 complex with its high-affinity receptor (IFN-gamma R alpha) with an affinity of 27(+/- 9) nM. The crystal structure of IFN-gamma SC1 has been determined at 2.9 A resolution from crystals grown in 1.4 M citrate solutions at pH 7.6. Comparison of the wild-type receptor-binding domain and the Asp111-containing domain of IFN-gamma SC1 show that they are structurally equivalent but have very different electrostatic surface potentials. As a result, surface charge rather than structural changes is likely responsible for the inability of the His111-->Asp domain of to bind IFN-gamma R alpha. The AB loops of IFN-gamma SC1 adopt conformations similar to the ordered loops of IFN-gamma observed in the crystal structure of the IFN-gamma/IFN-gamma R alpha complex. Thus, IFN-gamma R alpha binding does not result in a large conformational change in the AB loop as previously suggested. The structure also reveals the final six C-terminal amino acid residues of IFN-gamma SC1 (residues 253-258) that have not been observed in any other reported IFN-gamma structures. Despite binding to only one IFN-gamma R alpha, IFN-gamma SC1 is biologically active in cell proliferation, MHC class I induction, and anti-viral assays. This suggests that one domain of IFN-gamma is sufficient to recruit IFN-gamma R alpha and IFN-gamma R beta into a complex competent for eliciting biological activity. The current data are consistent with the main role of the IFN-gamma dimer being to decrease the dissociation constant of IFN-gamma for its cellular receptors.     

Regulation of yeast actin behavior by interaction of charged residues across the interdomain cleft

His(73) participates in the regulation of the nucleotide binding cleft conformation in yeast actin. Earlier molecular dynamics studies suggested that Asp(184) interacts with His(73) thereby stabilizing a "closed-cleft" G-actin. However, beta-actin in the open-cleft state shows a closer interaction of His(73) with Asp(179) than with Asp(184). We have thus assessed the relative importance of Asp(184) and Asp(179) on yeast actin stability and function. Neutral substitutions at 184 or 179 alone had little adverse effect on the monomer and polymerization behavior of actin. Arg or His at 184 in H73E actin partially rescued the monomeric properties of H73E actin, as demonstrated by near-normal thermostability and wild-type (WT)-like protease digestion patterns. ATP exchange was still considerably faster than with WT-actin although slower than that of H73E alone. However, polymerization of H73E/D184R and H73E/D184H is worse than with H73E alone. Conversely, D179R rescued all monomeric properties of H73E to near WT values and largely restored polymerization rate and filament thermostability. These results and new simulations of G-actin in the "open" state underscore the importance of the His(73)-Asp(179) interaction and suggest that the open and not the closed state of yeast actin may be favored in the absence of the methyl group of His(73).     

Disulfide bridge engineering in the tachykinin NK1 receptor

As in most other seven-transmembrane receptors, the central disulfide bridge from the extracellular end of TM-III to the middle of the second extracellular loop was essential for ligand binding in the NK1 receptor. However, introduction of "extra", single Cys residues in the second extracellular loop, at positions where disease-associated Cys substitutions impair receptor function in the vasopressin V2 receptor and in rhodopsin, did not cause mispairing with the Cys residues involved in this central disulfide bridge. Cys residues were introduced in the N-terminal extension and in the third extracellular loop, respectively, in such a way that disulfide bridge formation could be monitored by loss of substance P binding and breakage of the bridge could be monitored by gain of ligand binding. This disulfide bridge formed spontaneously in the whole population of receptors and could be titrated with low concentrations of reducing agent, dithiothreitol. Another putative disulfide bridge "switch" was constructed at the extracellular ends of TM-V and -VI, i.e., at positions where a high-affinity zinc site previously had been constructed with His substitutions. Disulfide bridge formation at this position, monitored by loss of binding of the nonpeptide antagonist [3H]LY303.870, occurred spontaneously only in a small fraction of the receptors. It is concluded that disulfide bridges form readily between Cys residues introduced appropriately in the N-terminal extension and the third extracellular loop, whereas they form with more difficulty between Cys residues placed at the extracellular ends of the transmembrane segments even at positions where high-affinity metal ion sites can be constructed with His residues.     

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.     

Decreased metallation and activity in subsets of mutant superoxide dismutases associated with familial amyotrophic lateral sclerosis

Over 90 different mutations in the gene encoding copper/zinc superoxide dismutase (SOD1) cause approximately 2% of amyotrophic lateral sclerosis (ALS) cases by an unknown mechanism. We engineered 14 different human ALS-related SOD1 mutants and obtained high yields of biologically metallated proteins from an Sf21 insect cell expression system. Both the wild type and mutant "as isolated" SOD1 variants were deficient in copper and were heterogeneous by native gel electrophoresis. By contrast, although three mutant SOD1s with substitutions near the metal binding sites (H46R, G85R, and D124V) were severely deficient in both copper and zinc ions, zinc deficiency was not a consistent feature shared by the as isolated mutants. Eight mutants (A4V, L38V, G41S, G72S, D76Y, D90A, G93A, and E133 Delta) exhibited normal SOD activity over pH 5.5-10.5, per equivalent of copper, consistent with the presumption that bound copper was in the proper metal-binding site and was fully active. The H48Q variant contained a high copper content yet was 100-fold less active than the wild type enzyme and exhibited a blue shift in the visible absorbance peak of bound Cu(II), indicating rearrangement of the Cu(II) coordination geometry. Further characterization of these as-isolated SOD1 proteins may provide new insights regarding mutant SOD1 enzyme toxicity in ALS.     

Proton transfer dynamics of GART: the pH-dependent catalytic mechanism examined by electrostatic calculations

The enzyme glycinamide ribonucleotide transformylase (GART) catalyzes the transfer of a formyl group from formyl tetrahydrofolate (fTHF) to glycinamide ribonucleotide (GAR), a process that is pH-dependent with pK(a) of approximately 8. Experimental studies of pH-rate profiles of wild-type and site-directed mutants of GART have led to the proposal that His108, Asp144, and GAR are involved in catalysis, with His108 being an acid catalyst, while forming a salt bridge with Asp144, and GAR being a nucleophile to attack the formyl group of fTHF. This model implied a protonated histidine with pK(a) of 9.7 and a neutral GAR with pK(a) of 6.8. These proposed unusual pK(a)s have led us to investigate the electrostatic environment of the active site of GART. We have used Poisson-Boltzmann-based electrostatic methods to calculate the pK(a)s of all ionizable groups, using the crystallographic structure of a ternary complex of GART involving the pseudosubstrate 5-deaza-5,6,7,8-THF (5dTHF) and substrate GAR. Theoretical mutation and deletion analogs have been constructed to elucidate pairwise electrostatic interactions between key ionizable sites within the catalytic site. Also, a construct of a more realistic catalytic site including a reconstructed pseudocofactor with an attached formyl group, in an environment with optimal local van der Waals interactions (locally minimized) that imitates closely the catalytic reactants, has been used for pK(a) calculations. Strong electrostatic coupling among catalytic residues His108, Asp144, and substrate GAR was observed, which is extremely sensitive to the initial protonation and imidazole ring flip state of His108 and small structural changes. We show that a proton can be exchanged between GAR and His108, depending on their relative geometry and their distance to Asp144, and when the proton is attached on His108, catalysis could be possible. Using the formylated locally minimized construct of GART, a high pK(a) for His108 was calculated, indicating a protonated histidine, and a low pK(a) for GAR(NH(2)) was calculated, indicating that GAR is in neutral form. Our results are in qualitative agreement with the current mechanistic picture of the catalytic process of GART deduced from the experimental data, but they do not reproduce the absolute magnitude of the pK(a)s extracted from fits of k(cat)-pH profiles, possibly because the static time-averaged crystallographic structure does not describe adequately the dynamic nature of the catalytic site during binding and catalysis. In addition, a strong effect on the pK(a) of GAR(NH(2)) is produced by the theoretical mutations of His108Ala and Asp144Ala, which is not in agreement with the observed insensitivity of the pK(a) of GAR(NH(2)) modeled from the experimental data using similar mutations. Finally, we show that important three-way electrostatic interactions between highly conserved His137, with His108 and Asp144, are responsible for stabilizing the electrostatic microenvironment of the catalytic site. In conclusion, our data suggest that further detailed computational and experimental work is necessary.     

Nonnative interactions in coupled folding and binding processes of intrinsically disordered proteins

Proteins function by interacting with other molecules, where both native and nonnative interactions play important roles. Native interactions contribute to the stability and specificity of a complex, whereas nonnative interactions mainly perturb the binding kinetics. For intrinsically disordered proteins (IDPs), which do not adopt rigid structures when being free in solution, the role of nonnative interactions may be more prominent in binding processes due to their high flexibilities. In this work, we investigated the effect of nonnative hydrophobic interactions on the coupled folding and binding processes of IDPs and its interplay with chain flexibility by conducting molecular dynamics simulations. Our results showed that the free-energy profiles became rugged, and intermediate states occurred when nonnative hydrophobic interactions were introduced. The binding rate was initially accelerated and subsequently dramatically decreased as the strength of the nonnative hydrophobic interactions increased. Both thermodynamic and kinetic analysis showed that disordered systems were more readily affected by nonnative interactions than ordered systems. Furthermore, it was demonstrated that the kinetic advantage of IDPs ("fly-casting" mechanism) was enhanced by nonnative hydrophobic interactions. The relationship between chain flexibility and protein aggregation is also discussed.