Analysis of a three-dimensional structure of human acidic mammalian chitinase obtained by homology modeling and ligand binding studies

The three-dimensional (3D) model of the human acidic mammalian chitinase (hAMCase) was constructed based on the crystal structure of the human chitotriosidase (EC 3.2.1.44, PDB code 1HKK) by using InsightII/Homology module. With the aid of molecular mechanics and molecular dynamics methods, the last refined model was obtained and further assessed by Profile-3D and Procheck, which confirms that the refined model is reliable. Furthermore, the docking results of the ligands (allosamidin and NAG₂) into the active site of hAMCase indicate that allosamidin is a more preferred ligand than NAG₂, and that Glu119 forms hydrogen bond with allosamidin, which is in good agreement with the experimental results. From the docking studies, we also suggest that Trp10, Glu49, Asp192, and Glu276 in hAMCase are four important determinant residues in binding as they have strong van-der-Waals and electrostatic interactions with the ligand, respectively.     

Homology modeling and docking study of cyclin-dependent kinase (CDK) 10

In order to understand the mechanisms of ligand binding and the interaction between the ligand and the cyclin-dependent kinase 10 (CDK10), a three-dimensional (3D) model of the CDK10 is generated based on the crystal structure of the cyclin-dependent kinase 2 (CDK2) (PDB code 1AQ1) by using InsightII/Homology module. With the aid of the molecular mechanics and molecular dynamics methods, the last refined model is obtained and is further assessed by PROFILE-3D and PROSTAT, which show that the refined model is reliable. With this model, a flexible docking study is performed and the results indicate that the Lys39 and Asp94 form hydrogen bonds and have strong nonbonding interaction with adenosine 5'-triphosphate (ATP). From the docking studies, we also suggest that the Leu141, Tyr21, and Val24 in CDK10 are three important determinant residues in binding as they have strong nonbonding interaction with ATP. The hydrogen bonding interactions also play an important role for the stability of the complex. Our results may be helpful for further experimental investigations.     

In silico modeling and hydrogen peroxide binding study of rice catalase

Homology modeling of the catalase, CatC cloned and sequenced from rice (Oryza sativa L., cv Ratna an Indica cultivar) has been performed based on the crystal structure of the catalase CatF (PDB code 1m7s) by using the software MODELLER. With the aid of molecular mechanics and molecular dynamics methods, the final model is obtained and is further assessed by PROCHECK and VERIFY - 3D graph, which show that the final refined model is reliable. With this model, a flexible docking study with the hydrogen peroxide, the substrate for catalase, is performed and the results indicate that Arg310, Asp343 and Arg346 in catalase are three important determinant residues in binding as they have strong hydrogen bonding contacts with the substrate. These hydrogen-bonding interactions play an important role for the stability of the complex. Our results may be helpful for further experimental investigations.     

Structural basis for different substrate specificities of two ADP-ribose pyrophosphatases from Thermus thermophilus HB8

ADP-ribose (ADPR) is one of the main substrates of Nudix proteins. Among the eight Nudix proteins of Thermus thermophilus HB8, we previously determined the crystal structure of Ndx4, an ADPR pyrophosphatase (ADPRase). In this study we show that Ndx2 of T. thermophilus also preferentially hydrolyzes ADPR and flavin adenine dinucleotide and have determined its crystal structure. We have determined the structures of Ndx2 alone and in complex with Mg2+, with Mg2+ and AMP, and with Mg2+ and a nonhydrolyzable ADPR analogue. Although Ndx2 recognizes the AMP moiety in a manner similar to those for other ADPRases, it recognizes the terminal ribose in a distinct manner. The residues responsible for the recognition of the substrate in Ndx2 are not conserved among ADPRases. This may reflect the diversity in substrate specificity among ADPRases. Based on these results, we propose the classification of ADPRases into two types: ADPRase-I enzymes, which exhibit high specificity for ADPR; and ADPRase-II enzymes, which exhibit low specificity for ADPR. In the active site of the ternary complexes, three Mg2+ ions are coordinated to the side chains of conserved glutamate residues and water molecules. Substitution of Glu90 and Glu94 with glutamine suggests that these residues are essential for catalysis. These results suggest that ADPRase-I and ADPRase-II enzymes have nearly identical catalytic mechanisms but different mechanisms of substrate recognition.     

Homology modeling and molecular dynamics study on N-acetylneuraminate lyase

With homology modeling techniques, molecular mechanics and molecular dynamics methods, a 3D structure model of N-acetylneuraminate lyase from human (hNAL, EC 4.1.3.3) was created and refined. This model was further assessed by Profile-3D and PROCHECK, which confirms that the refined model is reliable. Furthermore, the docking results of the substrates (sialic acid and KDO) into the active site of hNAL indicate that hNAL can cleave the sialic acid and KDO. Thr51 and Tyr143 may be the key amino acids residues as they have strong hydrogen bonding interactions with the substrates, which is in good agreement with the experimental results by Izard et al. (Structure 2:361–369. doi:10.1016/S0969-2126(00)00038-1 (1994)). From the docking studies, we also suggest that Asp176 and Ser218 only form hydrogen bonds with sialic acid, therefore, they may help sialic acid interact with hNAL steadly.     

Crystal structure of family GH-8 chitosanase with subclass II specificity from Bacillus sp. K17

Crystal structures of chitosanase from Bacillus sp. K17 (ChoK) have been determined at 1.5 A resolution in the active form and at 2.0 A resolution in the inactive form. This enzyme belongs to the family GH-8, out of 93 glycoside hydrolase families, and exhibits the substrate specificity of subclass II chitosanase. The catalytic site is constructed on the scaffold of a double-alpha(6)/alpha(6)-barrel, which is formed by six repeating helix-loop-helix motifs. This structure is quite different from those of the GH-46 chitosanases and of GH-5. Structural comparison with CelA (a cellulase belonging to the same family GH-8) suggests that the proton donor Glu122 is conserved, but the proton acceptor is the inserted Glu309 residue, and that the corresponding Asp278 residue in CelA is inactivated in ChoK. The four acidic residues, Asp179, Glu309, Asp183 and Glu107, can be involved in substrate recognition through interactions with the amino groups of the glucosamine residues bound in the -3, -2, -1 and +1 sites, respectively. The hydrophobic Trp235, Trp166, Phe413 and Tyr318 residues are highly conserved for binding of the hexose rings at the -3, -2, +1 and +2 sites, respectively. These structural features indicate that enzymes in GH-8 can be further divided into three subfamilies. Different types of chitosanases are discussed in terms of convergent evolution from different structural ancestors.     

The Nudix hydrolase Ndx1 from Thermus thermophilus HB8 is a diadenosine hexaphosphate hydrolase with a novel activity

The ndx1 gene, which encodes a Nudix protein, was cloned from the extremely thermophilic bacterium Thermus thermophilus HB8. This gene encodes a 126-amino acid protein that includes the characteristic Nudix motif conserved among Nudix proteins. Ndx1 was overexpressed in Escherichia coli and purified. Ndx1 was stable up to 95 degrees C and at extreme pH. Size exclusion chromatography indicated that Ndx1 was monomeric in solution. Ndx1 specifically hydrolyzed (di)adenosine polyphosphates but not ATP or diadenosine triphosphate, and it always generated ATP as the product. Diadenosine hexaphosphate (Ap(6)A), the most preferred substrate, was hydrolyzed to produce two ATP molecules, which is a novel hydrolysis mode for Ap(6)A, with a K(m) of 1.4 microm and a k(cat) of 4.1 s(-1). These results indicate that Ndx1 is a (di)adenosine polyphosphate hydrolase. Ndx1 activity required the presence of the divalent cations Mn(2+), Mg(2+), Zn(2+), and Co(2+), whereas Ca(2+), Ni(2+), and Cu(2+) were not able to activate Ndx1. Fluoride ion inhibited Ndx1 activity via a non-competitive mechanism. Optimal activity for Ap(6)A was observed at around pH 8.0 and about 70 degrees C. We found two important residues with pK(a) values of 6.1 and 9.6 in the free enzyme and pK(a) values of 7.9 and 10.0 in the substrate-enzyme complex. Kinetic studies of proteins with amino acid substitutions suggested that Glu-46 and Glu-50 were conserved residues in the Nudix motif and were involved in catalysis. Trp-26 was likely involved in enzyme-substrate interactions based on fluorescence measurements. Based on these results, the mechanism of substrate recognition and catalysis are discussed.     

Homology modelling and molecular dynamics study of human fatty acid amide hydrolase

The three-dimensional (3D) model of the human fatty acid amide hydrolase (hFAAH) was constructed based on the crystal structure of the rat FAAH (PDB code 1MT5) in complex with a substrate using Modeller9v2 program. With the aid of molecular mechanics and molecular dynamics method, the last model was obtained and further assessed by Profile-3D, Prosa2003 and Procheck, which confirms that the refined model is reliable. Furthermore, the docking results of propofol and its structural analogue into the active site of hFAAH indicate that 2,6-di-sec-butyl phenol is a more preferred ligand than others, which is in good agreement with the experimental results. From the docking studies, we also suggest that Phe192, Ile238, Thr377, Leu380, Phe381, Phe388 and Leu404 in the hFAAH are seven important determinant residues in binding as they have strong van der Waal interactions with the ligand.     

Homology modeling and molecular dynamics studies of a novel C3-like ADP-ribosyltransferase

The novel C3-like ADP-ribosyltransferase is produced by a Staphylococcus aureus strain that especially ADP-ribosylates RhoE/Rnd3 subtype proteins, and its three-dimensional (3D) structure has not known. In order to understand the catalytic mechanism, the 3D structure of the protein is built by using homology modeling based on the known crystal structure of exoenzyme C3 from Clostridium botulinum (1G24). Then the model structure is further refined by energy minimization and molecular dynamics methods. The putative nicotinamide adenine dinucleotide (NAD(+))-binding pocket of exoenzyme C3(Stau) is determined by Binding-Site Search module. The NAD(+)-enzyme complex is developed by molecular dynamics simulation and the key residues involved in the combination of enzyme binding to the ligand-NAD(+) are determined, which is helpful to guide the experimental realization and the new mutant designs as well. Our results indicated that the key binding-site residues of Arg48, Glu180, Ser138, Asn134, Arg85, and Gln179 play an important role in the catalysis of exoenzyme C3(Stau), which is in consistent with experimental observation.     

Molecular dynamics simulation of the P2Y14 receptor. Ligand docking and identification of a putative binding site of the distal hexose moiety

A rhodopsin-based homology model of the P2Y14 receptor was inserted into a phospholipid bilayer and refined by molecular dynamics (MD) simulation. The binding modes of several known agonists, namely UDP-glucose and its analogues, were proposed using automatic molecular docking combined with Monte Carlo Multiple Minimum calculations. Compared to other P2Y receptors, the P2Y14 receptor has an atypical binding mode of the nucleobase, ribose, and phosphate moieties. The diphosphate moiety interacts with only one cationic residue, namely Lys171 of EL2, while in other P2Y receptor subtypes three Arg or Lys residues interact with the phosphate chain. Two other conserved cationic residues, namely Arg253 (6.55) and Lys277 (7.35) of the P2Y14 receptor together with two anionic residues (Glu166 and Glu174, located in EL2), are likely involved in interactions with the distal hexose moiety.     

A computational H5N1 neuraminidase model and its binding to commercial drugs

In order to understand the mechanisms of ligand binding and interaction between two commercial drugs (ligands), zanamivir and oseltamivir and H5N1 Influenza Virus Neuraminidase subtype N1, a three-dimensional model of N1-ligand (GenBank accession no. AAS654617) was initially generated by homology modeling using the 13 high-resolution X-ray structures of neuraminidase N2 and N9 as the template. With the aid of the molecular mechanics and molecular dynamics methods, the final implicit solvent refined model was obtained. It was, then, assessed by PROCHECK, PROSA and VERIFY3D. With this model, a flexible docking study was performed. The results show strong hydrogen bond interactions between the glycerol side chains of zanamivir and Arg29 of the N1. Common hydrogen bonds between the carboxyl groups and Arg279 were found for both drugs. It was also found that the Glu30, Asp62, Arg63, Arg204, Trp310, Tyr313, Glu336, Ile338, Trp348, Ala349 were observed to facilitate the enzyme-ligand non-bonding interactions as they are located within the radius of 5 Å from all atoms of both drugs. Charge distribution was evaluated using the semi-empirical AM1 method. The results show that the total net charges of the –NH side chain of zanamivir is less negative than that of oseltamivir. This is in contrast to what is observed for the amide and alkyl (ether/glycerol) side chains. In comparison of the binding free energies between the X-ray N2-ligand and N9-ligand complexes, N1-ligand binding is found to be less potent than N2 and N9 subtypes, while N2-ligand and N9-ligand are roughly comparable. In addition, it is interesting to observe that the binding free energies for all three subtypes of the zanamivir complexes are lower than those of oseltamivir.     

Homology modeling and PAPS ligand (cofactor) binding study of bovine phenol sulfotransferase

In order to understand the mechanisms of ligand binding and the interaction between the ligand and the bovine phenol sulfotransferase, (bSULT1A1, EC 2.8.2.1) a three-dimensional (3D) model of the bSULT1A1 is generated based on the crystal structure of the estrogen sulfotransferase (PDB code 1AQU) by using the InsightII/Homology module. With the aid of the molecular mechanics and molecular dynamics methods, the final refined model is obtained and is further assessed by Profile-3D and ProStat, which show that the refined model is reliable. With this model, a flexible docking study is performed and the results indicate that 3-phosphoadenosine-5- phosphosulfate (PAPS) is a more preferred ligand than coenzyme A (CoA), and that His108 forms hydrogen bond with PAPS, which is in good agreement with the experimental results. From these docking studies, we also suggest that Phe255, Phe24 and Tyr169 in bSULT1A1 are three important determinant residues in binding as they have strong van-der-Waals contacts with the ligand. The hydrogen–bonding interactions also play an important role for the stability of the complex. Our results may be helpful for further experimental investigations.Figure The final 3D-structure of bSULT1A1. The structure is obtained by energy minimizing an average conformation over the last 100 ps of MD simulation. The -helix is represented in red and the -sheet in yellow.     

Steered molecular dynamics simulation of the binding of the β2 and β3 regions in domain-swapped human cystatin C dimer

The crystal structure of the human cystatin C (hCC) dimer revealed that a stable twofold-symmetric dimer was formed via 3D domain swapping. Domain swapping with the need for near-complete unfolding has been proposed as a possible route for amyloid fibril initiation. Thus, the interesting interactions that occur between the two molecules may be important for the further aggregation of the protein. In this work, we performed steered molecular dynamics (SMD) simulations to investigate the dissociation of the β2 and β3 strands in the hCC dimer. The energy changes observed during the SMD simulations showed that electrostatic interactions were the dominant interactions involved in stabilizing the two parts of the dimer during the early stages of SMD simulation, whereas van der Waals (VDW) interactions and electrostatic interactions were equally matched during the latter stages. Furthermore, our data indicated that the two parts of the dimer are stabilized by intermolecular hydrogen bonds among the residues Arg51 (β2), Gln48 (β2), Asp65 (β3), and Glu67 (β3), salt bridges among the residues Arg53 (β2), Arg51 (β2), and Asp65 (β3), and VDW interactions among the residues Gln48 (β2), Arg51 (β2), Glu67 (β3), Asp65 (β3), Phe63 (β3), and Asn61 (β3). The residues Gln48 (β2), Arg51 (β2), Asp65 (β3) and Glu67 (β3) appear to be crucial, as they play important roles in both electrostatic and VDW interactions. Thus, the present study determined the key residues involved in the stabilization of the domain-swapped dimer structure, and also provided molecular-level insights into the dissociation process of the hCC dimer.     

Bundles of amphipathic transmembrane alpha-helices as a structural motif for ion-conducting channel proteins: studies on sodium channels and acetylcholine receptors

Channel proteins are transmembrane symmetric (or pseudosymmetric) oligomers organized around a central ionic pore. We present here a molecular model of the pore forming structures of two channel proteins with different primary structures and oligomeric size: the voltage-sensitive sodium channel and the nicotinic cholinergic receptor. We report low-energy arrangements of alpha-helical bundles calculated by semiempiricial potential energy functions and optimization routines and further refined using molecular dynamics. The ion-conducting pore is considered to be a symmetric or pseudosymmetric homooligomer of 3-5 amphipathic alpha-helices arranged such that the polar residues line a central hydrophilic pathway and the apolar residues face the hydrophobic bilayer interior. The channel lining exposes either charged (Asp, Glu, Arg, Lys) or polar-neutral (Ser, Thr) residues. A bundle of four parallel helices constrained to C4 symmetry, the helix axis aligned with the symmetry axis, and the helices constrained to idealized dihedral angles, produces a structure with a pore of the size inferred for the sodium channel protein (area approximately 16 A2). Similarly, a pentameric array optimized with constraints to maintain C5 symmetry and backbone torsions characteristic of alpha-helices adopts a structure that appears well suited to form the lining of the nicotinic cholinergic receptor (pore area approximately 46 A2). Thus, bundles of amphipathic alpha-helices satisfy the structural, energetic, and dynamic requirements to be the molecular structural motif underlying the function of ionic channels.     

Refined X-ray structure of the low pH form of ribonuclease T1-2'-guanylic acid complex at 1.9 A resolution

The three-dimensional X-ray structure of the RNase T1[EC 3.1.27.3]-2'GMP complex crystallized at low pH value (4.0) was determined, and refined to 1.9 A resolution to give a final R value of 0.203. The refined model includes 781 protein atoms, 24 inhibitor atoms, and 43 solvent molecules. The imidazole rings of His27 and His40 interact with the carboxyl side chains of Glu82 and Glu58, respectively, whereas that of His92 is in contact with the main chain carbonyl oxygen of Ala75. In the complex, the ribose ring of the 2'GMP molecule adopts a C2'-endo puckering, and the exocyclic conformation is gauche(-)-gauche(+). The glycosyl torsion angle is in the syn range with an intramolecular hydrogen bond between N3 and O5', and the 2'-phosphate orientation is trans-gauche(-). The guanine base of the inhibitor is tightly bound to the base recognition site with five hydrogen bonds (N1--Glu46O epsilon 2, N2---Asn98O,O6---Asn44N, and N7 ---Asn43N delta 2/Asn43N) and is sandwiched between the phenolic ring portions of Tyr42 and Tyr45 by stacking interactions. The 2'-phosphate group interacts with Arg77N eta 2, Glu58O episilon 2, and Tyr 38O eta but not with any of the histidine residues. Arg77N eta 2 also interacts with Tyr38O eta. There is no interaction between the ribose moiety of the inhibitor and the enzyme.     

Structure prediction and R115866 binding study of human CYP26A1: homology modelling,fold recognition,molecular docking and MD simulations

Two three-dimensional (3D) models of human cytochrome P450 26A1 (CYP26A1) were constructed using the programs Modeller and Sybyl-GeneFold, respectively. After refinement by molecular mechanics and molecular dynamics (MD) simulations, the two models were validated by structure analysis-validation online server. Subsequently, a flexible docking study was performed on the model constructed by GeneFold with the potent and specific inhibitor R115866 to examine the enzyme–inhibitor interactions. From the docking results, we can see R115866 interacts with amino acid residues at the active site by multiple hydrophobic interactions including the side chains of His111, Trp112, Ser115, Val116, Leu125, Ser126, Leu221, Phe222, Glu296, Phe299, Gly300, Glu303, Thr304, Pro371 and the cofactor heme. Trp112 and Thr304 form hydrogen bonds with R115866 and play important roles in stabilising the complex. This constructed CYP26A1 model may provide an opportunity to understand the action mode of the enzyme and could be useful in designing novel retinoic acid metabolism blocking agents (RAMBAs).     

Structure-function map of the receptor site for β-scorpion toxins in domain II of voltage-gated sodium channels

Voltage-gated sodium (Na(v)) channels are the molecular targets of β-scorpion toxins, which shift the voltage dependence of activation to more negative membrane potentials by a voltage sensor-trapping mechanism. Molecular determinants of β-scorpion toxin (CssIV) binding and action on rat brain sodium channels are located in the S1-S2 (IIS1-S2) and S3-S4 (IIS3-S4) extracellular linkers of the voltage-sensing module in domain II. In IIS1-S2, mutations of two amino acid residues (Glu(779) and Pro(782)) significantly altered the toxin effect by reducing binding affinity. In IIS3-S4, six positions surrounding the key binding determinant, Gly(845), define a hot spot of high-impact residues. Two of these substitutions (A841N and L846A) reduced voltage sensor trapping. The other three substitutions (N842R, V843A, and E844N) increased voltage sensor trapping. These bidirectional effects suggest that the IIS3-S4 loop plays a primary role in determining both toxin affinity and efficacy. A high resolution molecular model constructed with the Rosetta-Membrane modeling system reveals interactions of amino acid residues in sodium channels that are crucial for toxin action with residues in CssIV that are required for its effects. In this model, the wedge-shaped CssIV inserts between the IIS1-S2 and IIS3-S4 loops of the voltage sensor, placing key amino acid residues in position to interact with binding partners in these extracellular loops. These results provide new molecular insights into the voltage sensor-trapping model of toxin action and further define the molecular requirements for the development of antagonists that can prevent or reverse toxicity of scorpion toxins.     

NMR determination of pKa values in α-synuclein

The intrinsically unfolded protein α-synuclein has an N-terminal domain with seven imperfect KTKEGV sequence repeats and a C-terminal domain with a large proportion of acidic residues. We characterized pK(a) values for all 26 sites in the protein that ionize below pH 7 using 2D (1) H-(15) N HSQC and 3D C(CO)NH NMR experiments. The N-terminal domain shows systematically lowered pK(a) values, suggesting weak electrostatic interactions between acidic and basic residues in the KTKEGV repeats. By contrast, the C-terminal domain shows elevated pK(a) values due to electrostatic repulsion between like charges. The effects are smaller but persist at physiological salt concentrations. For α-synuclein in the membrane-like environment of sodium dodecylsulfate (SDS) micelles, we characterized the pK(a) of His50, a residue of particular interest since it is flanked within one turn of the α-helix structure by the Parkinson's disease-linked mutants E46K and A53T. The pK(a) of His50 is raised by 1.4 pH units in the micelle-bound state. Titrations of His50 in the micelle-bound states of the E46K and A53T mutants show that the pK(a) shift is primarily due to interactions between the histidine and the sulfate groups of SDS, with electrostatic interactions between His50 and Glu46 playing a much smaller role. Our results indicate that the pK(a) values of uncomplexed α-synuclein differ significantly from random coil model peptides even though the protein is intrinsically unfolded. Due to the long-range nature of electrostatic interactions, charged residues in the α-synuclein sequence may help nucleate the folding of the protein into an α-helical structure and confer protection from misfolding.     

Rational proteomics V: structure-based mutagenesis has revealed key residues responsible for substrate recognition and catalysis by the dehydrogenase and isomerase activities in human 3beta-hydroxysteroid dehydrogenase/isomerase type 1

Mammalian 3beta-hydroxysteroid dehydrogenase/isomerase (3beta-HSD) is a member of the short chain dehydrogenase/reductase. It is a key steroidogenic enzyme that catalyzes the first step of the multienzyme pathway conversion of circulating dehydroepiandrosterone and pregnenolone to active steroid hormones. A three dimensional model of a ternary complex of human 3beta-HSD type 1 (3beta-HSD_1) with an NAD cofactor and androstenedione product has been developed based upon X-ray structures of the ternary complex of E. coli UDP-galactose 4-epimerase (UDPGE) with an NAD cofactor and substrate (PDB_AC: 1NAH) and the ternary complex of human type 1 17beta-hydroxysteroid dehydrogenase (17beta-HSD_1) with an NADP cofactor and androstenedione (PDB_AC: 1QYX). The dimeric structure of the enzyme was built from two monomer models of 3beta-HSD_1 by respective 3D superposition with A and B subunits of the dimeric structure of Streptococcus suis DTDP-D-glucose 4,6-dehydratase (PDB_AC: 1KEP). The 3D model structure of 3beta-HSD_1 has been successfully used for the rational design of mutagenic experiments to further elucidate the key substrate binding residues in the active site as well as the basis for dual function of the 3beta-HSD_1 enzyme. The structure based mutant enzymes, Asn100Ser, Asn100Ala, Glu126Leu, His232Ala, Ser322Ala and Asn323Leu, have been constructed and functionally characterized. The mutagenic experiments have confirmed the predicted roles of the His232 and Asn323 residues in recognition of the 17-keto group of the substrate and identified Asn100 and Glu126 residues as key residues that participate for the dehydrogenase and isomerization reactions, respectively.     

X-ray crystal structures of the oxidized and reduced forms of the rubredoxin from the marine hyperthermophilic archaebacterium Pyrococcus furiosus.

The structures of the oxidized and reduced forms of the rubredoxin from the archaebacterium, Pyrococcus furiosus, an organism that grows optimally at 100 degrees C, have been determined by X-ray crystallography to a resolution of 1.8 A. Crystals of this rubredoxin grow in space group P2(1)2(1)2(1) with room temperature cell dimensions a = 34.6 A, b = 35.5 A, and c = 44.4 A. Initial phases were determined by the method of molecular replacement using the oxidized form of the rubredoxin from the mesophilic eubacterium, Clostridium pasteurianum, as a starting model. The oxidized and reduced models of P. furiosus rubredoxin each contain 414 nonhydrogen protein atoms comprising 53 residues. The model of the oxidized form contains 61 solvent H2O oxygen atoms and has been refined with X-PLOR and TNT to a final R = 0.178 with root mean square (rms) deviations from ideality in bond distances and bond angles of 0.014 A and 2.06 degrees, respectively. The model of the reduced form contains 37 solvent H2O oxygen atoms and has been refined to R = 0.193 with rms deviations from ideality in bond lengths of 0.012 A and in bond angles of 1.95 degrees. The overall structure of P. furiosus rubredoxin is similar to the structures of mesophilic rubredoxins, with the exception of a more extensive hydrogen-bonding network in the beta-sheet region and multiple electrostatic interactions (salt bridge, hydrogen bonds) of the Glu 14 side chain with groups on three other residues (the amino-terminal nitrogen of Ala 1; the indole nitrogen of Trp 3; and the amide nitrogen group of Phe 29). The influence of these and other features upon the thermostability of the P. furiosus protein is discussed.