Biomedical Informatics Training at the University of Wisconsin-Madison

The purpose of this paper is to describe biomedical informatics training at the University of Wisconsin-Madison (UW Madison).

We reviewed biomedical informatics training, research, and faculty/trainee participation at UW-Madison.

There are three primary approaches to training 1) The Computation & Informatics in Biology & Medicine Training Program, 2) formal biomedical informatics offered by various campus departments, and 3) individualized programs. Training at UW-Madison embodies the features of effective biomedical informatics training recommended by the American College of Medical Informatics that were delineated as: 1) curricula that integrate experiences among computational sciences and application domains, 2) individualized and interdisciplinary cross training among adiverse cadre of trainees to develop key competencies that he or she does not initially possess, 3) participation in research and development activities, and 4) exposure to a range of basic informational and computational sciences.

The three biomedical informatics training approaches immerse students in multidisciplinary training and education that is supported by faculty trainers who participate in collaborative research across departments. Training is provided across a range of disciplines and available at different training stages. Biomedical informatics training at UW-Madison illustrates how a large research University, with multiple departments across biological, computational and health fields, can provide effective and productive biomedical informatics training via multiple bioinformatics training approaches.

Ligand pathways in myoglobin: A review of trp cavity mutations

The pathways for ligand entry and exit in myoglobin have now been well established by a wide variety of experimental results, including pico- to nano- to microsecond transient absorbance measurements and time-resolved X-ray crystallographic measurements. Trp insertions have been used to block, one at a time, the three major cavities occupied by photodissociated ligands. In this work, we review the effects of the L29(B10)W mutation, which places a large indole ring in the initial 'docking site' for photodissociated ligands. Then, the effects of blocking the Xe4 site with I28W, V68W, and I107W mutations and the Xe1 cavity with L89W, L104W, and F138W mutations are described. The structures of four of these mutants are shown for the first time (Trp28, Trp68, Trp107, and Trp 138 sperm whale metMb). All available results support a 'side path' mechanism in which ligands move into and out of myoglobin by outward rotation of the HisE7 side chain, but after entry can migrate into internal cavities, including the distal Xe4 and proximal Xe1 binding sites. The distal cavities act like the pocket of a baseball glove, catching the ligand and holding it long enough for the histidine gate to close and facilitate internal coordination with the heme iron atom. The physiological role of the proximal Xe1 site is less clear because changes in the size of this cavity have minimal effects on overall O(2) binding parameters.

Biomedical informatics training at the University of Wisconsin-Madison

OBJECTIVES

The purpose of this paper is to describe biomedical informatics training at the University of Wisconsin-Madison (UW-Madison).

METHODS

We reviewed biomedical informatics training, research, and faculty/trainee participation at UW-Madison.

RESULTS

There are three primary approaches to training 1) The Computation & Informatics in Biology & Medicine Training Program, 2) formal biomedical informatics offered by various campus departments, and 3) individualized programs. Training at UW-Madison embodies the features of effective biomedical informatics training recommended by the American College of Medical Informatics that were delineated as: 1) curricula that integrate experiences among computational sciences and application domains, 2) individualized and interdisciplinary cross-training among a diverse cadre of trainees to develop key competencies that he or she does not initially possess, 3) participation in research and development activities, and 4) exposure to a range of basic informational and computational sciences.

CONCLUSIONS

The three biomedical informatics training approaches immerse students in multidisciplinary training and education that is supported by faculty trainers who participate in collaborative research across departments. Training is provided across a range of disciplines and available at different training stages. Biomedical informatics training at UW-Madison illustrates how a large research University, with multiple departments across biological, computational and health fields, can provide effective and productive biomedical informatics training via multiple bioinformatics training approaches.

Analysis of diffuse scattering from yeast initiator tRNA crystals

Yeast initiator tRNA crystals exhibit strong X-ray diffuse scattering. This scattering can be used to extract information about lattice-coupled and intramolecular motions in the crystals. The amplitudes and correlation distances of these motions can be estimated by calculating the diffuse scattering and comparing the results with the observed scattering. Results indicate that both anisotropic, lattice-coupled motions as well as short-range correlated local disorder in the anticodon arm contribute to the overall disorder in the crystals. These types of motions can be correlated with aspects of tRNA function. This additional information complements the results from analysis of crystallographic data and provides a more detailed picture of the structure and dynamics of the molecule. The degree to which the methodology presented here can account for the observed diffuse scattering from tRNA represents a significant step forward in the ability to use this conventionally discarded information, and encourages the ultimate extension of these ideas to a wide variety of macromolecular systems.

Role of CO2 in the formation of gold deposits

Much of global gold production has come from deposits with uneconomic concentrations of base metals, such as copper, lead and zinc. These 'gold-only' deposits are thought to have formed from hot, aqueous fluids rich in carbon dioxide, but only minor significance has been attached to the role of the CO2 in the process of gold transport. This is because chemical bonding between gold ions and CO2 species is not strong, and so it is unlikely that CO2 has a direct role in gold transport. An alternative indirect role for CO2 as a weak acid that buffers pH has also appeared unlikely, because previously inferred pH values for such gold-bearing fluids are variable. Here we show that such calculated pH values are unlikely to record conditions of gold transport, and propose that CO2 may play a critical role during gold transport by buffering the fluid in a pH range where elevated gold concentration can be maintained by complexation with reduced sulphur. Our conclusions, which are supported by geochemical modelling, may provide a platform for new gold exploration methods.

A fast Newton method for entropy maximization in statistical phase estimation

A fast Newton method is presented for solving the entropy maximization problem in the Bayesian statistical approach to phase estimation. The method requires only O(n log n) instead of standard O(n3) floating point operations per iteration, while converging in the same rate as the standard Newton method. The method is described and related computational issues are discussed. Numerical results on simple test cases are also presented to demonstrate the behavior of the method.

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.

How the CO in myoglobin acquired its bend: lessons in interpretation of crystallographic data

Contrary to the expectation of chemists, the first X-ray structures of carbon monoxide bound to myoglobin (Mb) showed a highly distorted Fe-C-O bond system. These results appeared to support the idea of a largely steric mechanism for discrimination by the protein against CO binding, a lethal act for the protein in terms of its physiological function. The most recent independently determined high-resolution structures of Mb-CO have allowed the 25 year old controversy concerning the mode of CO binding to be resolved. The CO is now seen to bind in a roughly linear fashion without substantial bending, consistent with chemical expectations and spectroscopic measurements. Access to deposited diffraction data prompted a reevaluation of the sources of the original misinterpretation. A series of careful refinements of models against the data at high (1.1 A) and modest resolutions (1.5 A) have been performed in anisotropic versus isotropic modes. The results suggest that the original artifact was a result of lower quality crystals combined with anisotropic motion and limited resolution of the diffraction data sets. This retrospective analysis should serve as a caution for all researchers using structural tools to draw far-reaching biochemical conclusions.

Waterproofing the heme pocket. Role of proximal amino acid side chains in preventing hemin loss from myoglobin

The ability of myoglobin to bind oxygen reversibly depends critically on retention of the heme prosthetic group. Globin side chains at the Leu(89)(F4), His(97)(FG3), Ile(99)(FG5), and Leu(104)(G5) positions on the proximal side of the heme pocket strongly influence heme affinity. The roles of these amino acids in preventing heme loss have been examined by determining high resolution structures of 14 different mutants at these positions using x-ray crystallography. Leu(89) and His(97) are important surface amino acids that interact either sterically or electrostatically with the edges of the porphyrin ring. Ile(99) and Leu(104) are located in the interior region of the proximal pocket beneath ring C of the heme prosthetic group. The apolar amino acids Leu(89), Ile(99), and Leu(104) "waterproof" the heme pocket by forming a barrier to solvent penetration, minimizing the size of the proximal cavity, and maintaining a hydrophobic environment. Substitutions with smaller or polar side chains at these positions result in exposure of the heme to solvent, the appearance of crystallographically defined water molecules in or near the proximal pocket, and large increases in the rate of hemin loss. Thus, the naturally occurring amino acid side chains at these positions serve to prevent hydration of the His(93)-Fe(III) bond and are highly conserved in all known myoglobins and hemoglobins.

Crystal structure of a nonsymbiotic plant hemoglobin

BACKGROUND

Nonsymbiotic hemoglobins (nsHbs) form a new class of plant proteins that is distinct genetically and structurally from leghemoglobins. They are found ubiquitously in plants and are expressed in low concentrations in a variety of tissues including roots and leaves. Their function involves a biochemical response to growth under limited O(2) conditions.

RESULTS

The first X-ray crystal structure of a member of this class of proteins, riceHb1, has been determined to 2.4 A resolution using a combination of phasing techniques. The active site of ferric riceHb1 differs significantly from those of traditional hemoglobins and myoglobins. The proximal and distal histidine sidechains coordinate directly to the heme iron, forming a hemichrome with spectral properties similar to those of cytochrome b(5). The crystal structure also shows that riceHb1 is a dimer with a novel interface formed by close contacts between the G helix and the region between the B and C helices of the partner subunit.

CONCLUSIONS

The bis-histidyl heme coordination found in riceHb1 is unusual for a protein that binds O(2) reversibly. However, the distal His73 is rapidly displaced by ferrous ligands, and the overall O(2) affinity is ultra-high (K(D) approximately 1 nM). Our crystallographic model suggests that ligand binding occurs by an upward and outward movement of the E helix, concomitant dissociation of the distal histidine, possible repacking of the CD corner and folding of the D helix. Although the functional relevance of quaternary structure in nsHbs is unclear, the role of two conserved residues in stabilizing the dimer interface has been identified.

Crystal structure of tropomyosin at 7 Angstroms resolution

Tropomyosin is a 400A-long coiled coil that polymerizes to form a continuous filament that associates with actin in muscle and numerous non-muscle cells. Tropomyosin and troponin together form a calcium-sensitive switch that is responsible for thin-filament regulation of striated muscle. Subtle structural features of the molecule, including non-canonical aspects of its coiled-coil motif, undoubtedly influence its association with f-actin and its role in thin filament regulation. Previously, careful inspection of native diffraction intensities was sufficient to construct a model of tropomyosin at 9A resolution in a spermine-induced crystal form that diffracts anisotropically to 4A resolution. Single isomorphous replacement (SIR) phasing has now provided an empirical determination of the structure at 7A resolution. A novel method of heavy-atom analysis was used to overcome difficulties in interpretation of extremely anisotropic diffraction. The packing arrangement of the molecules in the crystal, and important aspects of the tropomyosin geometry such as non-uniformities of the pitch and variable bending and radius of the coiled coil are evident.

Protein structure determination using a database of interatomic distance probabilities

The accelerated pace of genomic sequencing has increased the demand for structural models of gene products. Improved quantitative methods are needed to study the many systems (e.g., macromolecular assemblies) for which data are scarce. Here, we describe a new molecular dynamics method for protein structure determination and molecular modeling. An energy function, or database potential, is derived from distributions of interatomic distances obtained from a database of known structures. X-ray crystal structures are refined by molecular dynamics with the new energy function replacing the Van der Waals potential. Compared to standard methods, this method improved the atomic positions, interatomic distances, and side-chain dihedral angles of structures randomized to mimic the early stages of refinement. The greatest enhancement in side-chain placement was observed for groups that are characteristically buried. More accurate calculated model phases will follow from improved interatomic distances. Details usually seen only in high-resolution refinements were improved, as is shown by an R-factor analysis. The improvements were greatest when refinements were carried out using X-ray data truncated at 3.5 A. The database potential should therefore be a valuable tool for determining X-ray structures, especially when only low-resolution data are available.

Protein-assisted pericyclic reactions: an alternate hypothesis for the action of quantal receptors

The rules for allowable pericyclic reactions indicate that the photoisomerizations of retinals in rhodopsins can be formally analogous to thermally promoted Diels-Alder condensations of monoenes with retinols. With little change in the seven-transmembrane helical environment these latter reactions could mimic the retinal isomerization while providing highly sensitive chemical reception. In this way archaic progenitors of G-protein-coupled chemical quantal receptors such as those for pheromones might have been evolutionarily plagiarized from the photon quantal receptor, rhodopsin, or vice versa. We investigated whether the known structure of bacteriorhodopsin exhibited any similarity in its active site with those of the two known antibody catalysts of Diels-Alder reactions and that of the photoactive yellow protein. A remarkable three-dimensional motif of aromatic side chains emerged in all four proteins despite the drastic differences in backbone structure. Molecular orbital calculations supported the possibility of transient pericyclic reactions as part of the isomerization-signal transduction mechanisms in both bacteriorhodopsin and the photoactive yellow protein. It appears that reactions in all four of the proteins investigated may be biological analogs of the organic chemists' chiral auxiliary-aided Diels-Alder reactions. Thus the light receptor and the chemical receptor subfamilies of the heptahelical receptor family may have been unified at one time by underlying pericyclic chemistry.

Metal-ion affinity and specificity in EF-hand proteins: coordination geometry and domain plasticity in parvalbumin

BACKGROUND

The EF-hand family is a large set of Ca(2+)-binding proteins that contain characteristic helix-loop-helix binding motifs that are highly conserved in sequence. Members of this family include parvalbumin and many prominent regulatory proteins such as calmodulin and troponin C. EF-hand proteins are involved in a variety of physiological processes including cell-cycle regulation, second messenger production, muscle contraction, microtubule organization and vision.

RESULTS

We have determined the structures of parvalbumin mutants designed to explore the role of the last coordinating residue of the Ca(2+)-binding loop. An E101D substitution has been made in the parvalbumin EF site. The substitution decreases the Ca(2+)-binding affinity 100-fold and increases the Mg(2+)-binding affinity 10-fold. Both the Ca(2+)- and Mg(2+)-bound structures have been determined, and a structural basis has been proposed for the metal-ion-binding properties.

CONCLUSIONS

The E101D mutation does not affect the Mg(2+) coordination geometry of the binding loop, but it does pull the F helix 1.1 A towards the loop. The E101D-Ca(2+) structure reveals that this mutant cannot obtain the sevenfold coordination preferred by Ca(2+), presumably because of strain limits imposed by tertiary structure. Analysis of these results relative to previously reported structural information supports a model wherein the characteristics of the last coordinating residue and the plasticity of the Ca(2+)-binding loop delimit the allowable geometries for the coordinating sphere.

Effects of the location of distal histidine in the reaction of myoglobin with hydrogen peroxide

To clarify how the location of distal histidine affects the activation process of H2O2 by heme proteins, we have characterized reactions with H2O2 for the L29H/H64L and F43H/H64L mutants of sperm whale myoglobin (Mb), designed to locate the histidine farther from the heme iron. Whereas the L29H/H64L double substitution retarded the reaction with H2O2, an 11-fold rate increase versus wild-type Mb was observed for the F43H/H64L mutant. The Vmax values for 1-electron oxidations by the myoglobins correlate well with the varied reactivities with H2O2. The functions of the distal histidine as a general acid-base catalyst were examined based on the reactions with cumene hydroperoxide and cyanide, and only the histidine in F43H/H64L Mb was suggested to facilitate heterolysis of the peroxide bond. The x-ray crystal structures of the mutants confirmed that the distal histidines in F43H/H64L Mb and peroxidase are similar in distance from the heme iron, whereas the distal histidine in L29H/H64L Mb is located too far to enhance heterolysis. Our results indicate that the proper positioning of the distal histidine is essential for the activation of H2O2 by heme enzymes.

Crystal structures of Bacillus stearothermophilus adenylate kinase with bound Ap5A, Mg2+ Ap5A, and Mn2+ Ap5A reveal an intermediate lid position and six coordinate octahedral geometry for bound Mg2+ and Mn2+

Crystal structures of Bacillus stearothermophilus adenylate kinase with bound Ap5A, Mn2+ Ap5A, and Mg2+ Ap5A have been determined by X-ray crystallography to resolutions of 1.6 A, 1.85 A, and 1.96 A, respectively. The protein's lid domain is partially open, being both rotated and translated away from bound Ap5A. The flexibility of the lid domain in the ternary state and its ability to transfer force directly to the the active site is discussed in light of our proposed entropic mechanism for catalytic turnover. The bound Zn2+ atom is demonstrably structural in nature, with no contacts other than its ligating cysteine residues within 5 A. The B. stearothermophilus adenylate kinase lid appears to be a truncated zinc finger domain, lacking the DNA binding finger, which we have termed a zinc knuckle domain. In the Mg2+ Ap5A and Mn2+ Ap5A structures, Mg2+ and Mn2+ demonstrate six coordinate octahedral geometry. The interactions of the Mg2+-coordinated water molecules with the protein and Ap5A phosphate chain demonstrate their involvement in catalyzing phosphate transfer. The protein selects for beta-y (preferred by Mg2+) rather than alpha-gamma (preferred by Mn2+) metal ion coordination by forcing the ATP phosphate chain to have an extended conformation.

Characterizing global substates of myoglobin

BACKGROUND

The massive amount of information generated from current molecular dynamics simulations makes the data difficult to analyze efficiently. Principal component analysis has been used for almost a century to detect and characterize data relationships and to reduce the dimensionality for problems in many fields. Here, we present an adaptation of principal component analysis using a partial singular value decomposition (SVD) for investigating both the localized and global motions of macromolecules.

RESULTS

Configuration space projections from the SVD analysis of a variety of myoglobin simulations are used to characterize the dynamics of the protein. This technique reveals new dynamical motifs, which quantify proposed hierarchical structures of conformational substates for proteins and provide a means by which configuration space sampling efficiency may be probed. The SVD clearly shows that solvent effects facilitate transitions between global conformational substates for myoglobin molecular dynamics simulations. Lyapunov exponents calculated from the configuration space divergence of 15 trajectories agree with previous predictions for the chaotic behavior of complex protein systems.

CONCLUSIONS

Configuration space projections provide invaluable information about protein motions that would be extremely difficult to obtain otherwise. While the configuration space for myoglobin is quite large, it does have structure. Our analysis of this structure shows that the protein hops between a number of distinct global conformational states, much like the local behavior observed for an individual residue.

Solution and crystal structures of a sperm whale myoglobin triple mutant that mimics the sulfide-binding hemoglobin from Lucina pectinata

The bivalve mollusc Lucina pectinata harbors sulfide-oxidizing chemoautotrophic bacteria and expresses a monomeric hemoglobin I, HbI, with normal O2, but extraordinarily high sulfide affinity. The crystal structure of aquomet Lucina HbI has revealed an active site with three residues not commonly found in vertebrate globins: Phe(B10), Gln(E7), and Phe(E11) (Rizzi, M., Wittenberg, J. B., Coda, A., Fasano, M., Ascenzi, P., and Bolognesi, M. (1994) J. Mol. Biol. 244, 86-89). Engineering these three residues into sperm whale myoglobin results in a triple mutant with approximately 700-fold higher sulfide affinity than for wild-type. The single crystal x-ray structure of the aquomet derivative of the myoglobin triple mutant and the solution 1H NMR active site structures of the cyanomet derivatives of both the myoglobin mutant and Lucina HbI have been determined to examine further the structural origin of their unusually high sulfide affinities. The major differences in the distal pocket is that in the aquomet form the carbonyl of Gln64(E7) serves as a H-bond acceptor, whereas in the cyanomet form the amido group acts as H-bond donor to the bound ligand. Phe68(E11) is rotated approximately 90 degrees about chi2 and located approximately 1-2 A closer to the iron atom in the myoglobin triple mutant relative to its conformation in Lucina HbI. The change in orientation potentially eliminates the stabilizing interaction with sulfide and, together with the decrease in size of the distal pocket, accounts for the 7-fold lower sulfide affinity of the myoglobin mutant compared with that of Lucina HbI.

Nitric oxide myoglobin: crystal structure and analysis of ligand geometry

The structure of the ferrous nitric oxide form of native sperm whale myoglobin has been determined by X-ray crystallography to 1.7 angstroms resolution. The nitric oxide ligand is bent with respect to the heme plane: the Fe-N-O angle is 112 degrees. This angle is smaller than those observed in model compounds and in lupin leghemoglobin. The exact angle appears to be influenced by the strength of the proximal bond and hydrogen bonding interactions between the distal histidine and the bound ligand. Specifically, the N(epsilon) atom of histidine64 is located 2.8 angstroms away from the nitrogen atom of the bound ligand, implying electrostatic stabilization of the FeNO complex. This interpretation is supported by mutagenesis studies. When histidine64 is replaced with apolar amino acids, the rate of nitric oxide dissociation from myoglobin increases tenfold.

Motions of calmodulin characterized using both Bragg and diffuse X-ray scattering

BACKGROUND

Calmodulin is a calcium-activated regulatory protein which can bind to many different targets. The protein resembles a highly flexible dumbbell, and bends in the middle as it binds. This and other motions must be understood to formulate a realistic model of calmodulin function.

RESULTS

Using the Bragg reflections from X-ray crystallography, a multiple-conformer refinement of a calmodulin-peptide complex shows anisotropic displacements, with high variations of dihedral angles in several nonhelical domains: the flexible linker; three of the four calcium-binding sites (including both of the N-terminal sites); and a turn connecting the C-terminal EF-hand calcium-binding domains. Three-dimensional maps of the large scale diffuse X-ray scattering data show isotropic liquid-like motions with an unusually small correlation length. Three-dimensional maps of the small scale diffuse streaks show highly coupled, anisotropic motions along the head-to-tail molecular packing direction in the unit cell. There is also weak coupling perpendicular to the head-to-tail packing direction, particularly across a cavity occupied by the disordered linker domain of the molecule.

CONCLUSIONS

Together, the Bragg and diffuse scattering present a self-consistent description of the motions in the flexible linker of calmodulin. The other mobile regions of the protein are also of great interest. In particular, the high variations in the calcium-binding sites are likely to influence how strongly they bind ions. This is especially important in the N-terminal sites, which regulate the activity of the molecule.

Structure and dynamics of green fluorescent protein

Many marine organisms are luminescent. The proteins that produce the light include a primary light producer (aequorin or luciferase) and often a secondary photoprotein that red shifts the light for better penetration in the ocean. Green fluorescent protein is one such secondary protein. It is remarkable in that it autocatalyzes the formation of its own fluorophore and thus can be expressed in a variety of organisms in its fluorescent form. The recent determination of its 3D structure and other physical characterizations are revealing its molecular mechanism of action.

Characterization of recombinant soybean leghemoglobin a and apolar distal histidine mutants

The cDNA for soybean leghemoglobin a (Lba) was cloned from a root nodule cDNA library and expressed in Escherichia coli. The crystal structure of the ferric acetate complex of recombinant wild-type Lba was determined at a resolution of 2.2 A. Rate constants for O2, CO and NO binding to recombinant Lba are identical with those of native soybean Lba. Rate constants for hemin dissociation and auto-oxidation of wild-type Lba were compared with those of sperm whale myoglobin. At 37 degrees C and pH 7, soybean Lba is much less stable than sperm whale myoglobin due both to a fourfold higher rate of auto-oxidation and to a approximately 600-fold lower affinity for hemin. The role of His61(E7) in regulating oxygen binding was examined by site-directed mutagenesis. Replacement of His(E7) with Ala, Val or Leu causes little change in the equilibrium constant for O2 binding to soybean Lba, whereas the same mutations in sperm whale myoglobin cause 50 to 100-fold decreases in K(O2). These results show that, at neutral pH, hydrogen bonding with His(E7) is much less important in regulating O2 binding to the soybean protein. The His(E7) to Phe mutation does cause a significant decrease in K(O2) for Lba, apparently due to steric hindrance of the bound ligand. The rate constants for O2 dissociation from wild-type and native Lba decrease significantly with decreasing pH. In contrast, the O2 dissociation rate constants for mutants with apolar E7 residues are independent of pH, suggesting that hydrogen bonding to the distal histidine residue in the native protein is enhanced under acid conditions. All of these results support the hypothesis that the high affinity of Lba for oxygen and other ligands is determined primarily by enhanced accessibility and reactivity of the heme group.

High resolution crystal structures of the deoxy, oxy, and aquomet forms of cobalt myoglobin

The structures of the deoxy, oxy, and aquomet forms of native sperm whale myoglobin reconstituted with cobalt protoporphyrin IX have been determined by x-ray crystallography. As expected, cobalt myoglobin closely resembles native iron myoglobin in overall structure, especially in their respective aquomet forms. In the cobalt oxymyoglobin structure, the Nepsilon of distal histidine 64 lies within hydrogen bonding distance to both the oxygen atom directly bonded to the cobalt and the terminal oxygen atom, in agreement with previous EPR and resonance Raman studies. The metal atom in cobaltous myoglobin does show a small 0.06-A out-of-porphyrin plane displacement when moving from the oxy to deoxy state. In the case of the native iron-containing myoglobin, the oxy to deoxy transition results in a larger 0.16-A displacement of the metal farther out of the porphyrin plane, attributed to an increase in spin from S = 0 to S = 2. The small displacement in cobalt myoglobin is due to a change in coordination geometry, not spin state (S = 1/2 for both cobalt deoxy- and oxymyoglobin). The small out-of-porphyrin plane movement of cobalt which accompanies deoxygenation of myoglobin also occurs in cobalt hemoglobin and serves to explain why cooperativity, although reduced, is still preserved when iron is replaced by cobalt in human hemoglobin.

The molecular structure of green fluorescent protein

The crystal structure of recombinant wild-type green fluorescent protein (GFP) has been solved to a resolution of 1.9 A by multiwavelength anomalous dispersion phasing methods. The protein is in the shape of a cylinder, comprising 11 strands of beta-sheet with an alpha-helix inside and short helical segments on the ends of the cylinder. This motif, with beta-structure on the outside and alpha-helix on the inside, represents a new protein fold, which we have named the beta-can. Two protomers pack closely together to form a dimer in the crystal. The fluorophores are protected inside the cylinders, and their structures are consistent with the formation of aromatic systems made up of Tyr66 with reduction of its C alpha-C beta bond coupled with cyclization of the neighboring glycine and serine residues. The environment inside the cylinder explains the effects of many existing mutants of GFP and suggests specific side chains that could be modified to change the spectral properties of GFP. Furthermore, the identification of the dimer contacts may allow mutagenic control of the state of assembly of the protein.

Affinity and structure of complexes of tropomyosin and caldesmon domains

The interaction of caldesmon domains with tropomyosin has been studied using x-ray crystallography and an optical biosensor. Only whole caldesmon and the carboxyl-terminal domain of caldesmon (CaD-4, chicken gizzard residues 597-756) bound to tropomyosin with greater than millimolar affinity at 100 and 150 microM salt. Under these conditions the affinities of whole caldesmon and CaD-4 were both in the micromolar range. Data from the x-ray studies showed that whole caldesmon bound to tropomyosin in several places, with the region of tightest interaction being at tropomyosin residues 70-100 and/or 230-260. Studies with CaD-4 revealed that this region corresponded to the strong binding site seen with whole caldesmon. Weaker association of other regions of caldesmon to tropomyosin residues 180-210 and 5-50 was also observed. The results suggest that the carboxyl-terminus of caldesmon binds tightly to tropomyosin and that other regions of caldesmon may interact with tropomyosin tightly only when they are held close to tropomyosin by the carboxyl-terminal domain. Four models are presented to show the possible interactions of caldesmon with tropomyosin.

Mechanism of NO-induced oxidation of myoglobin and hemoglobin

Nitric oxide (NO) has been implicated as mediator in a variety of physiological functions, including neurotransmission, platelet aggregation, macrophage function, and vasodilation. The consumption of NO by extracellular hemoglobin and subsequent vasoconstriction have been suggested to be the cause of the mild hypertensive events reported during in vivo trials of hemoglobin-based O2 carriers. The depletion of NO from endothelial cells is most likely due to the oxidative reaction of NO with oxyhemoglobin in arterioles and surrounding tissue. In order to determine the mechanism of this key reaction, we have measured the kinetics of NO-induced oxidation of a variety of different recombinant sperm whale myoglobins (Mb) and human hemoglobins (Hb). The observed rates depend linearly on [NO] but show no dependence on [O2]. The bimolecular rate constants for NO-induced oxidation of MbO2 and HbO2 are large (k.ox,NO = 30-50 microM-1 s-1 for the wild-type proteins) and similar to those for simple nitric oxide binding to deoxygenated Mb and Hb. Both reversible NO binding and NO-induced oxidation occur in two steps: (1) bimolecular entry of nitric oxide into the distal portion of the heme pocket and (2) rapid reaction of noncovalently bound nitric oxide with the iron atom to produce Fe(2+)-N=O or with Fe(2+)-O-O delta- to produce Fe(3+)-OH2 and nitrate. Both the oxidation and binding rate constants for sperm whale Mb were increased when His(E7) was replaced by aliphatic residues. These mutants lack polar interactions in the distal pocket which normally hinder NO entry into the protein. Decreasing the volume of the distal pocket by replacing Leu(B10) and Val(E11) with aromatic amino acids markedly inhibits NO-induced oxidation of MbO2. The latter results provide a protein engineering strategy for reducing hypertensive events caused by extracellular hemoglobin-based O2 carriers. This approach has been explored by examining the effects of Phe(B10) and Phe(E11) substitutions on the rates of NO-induced oxidation of the alpha and beta subunits in recombinant human hemoglobin.

Study of global motions in proteins by weighted masses molecular dynamics: adenylate kinase as a test case

The weighted masses molecular dynamics (WMMD) technique is applied to the protein adenylate kinase. A novel set of restraints has been developed to allow the use of this technique with proteins. The WMMD simulation is successful in predicting the flexibility of the two mobile domains of the protein. The end product of the simulation is similar to the known open and AMP bound forms of the enzyme. The biological relevance of the restraints used and potential methods of improving the technique are discussed.

Crystal structures of CO-, deoxy- and met-myoglobins at various pH values

The distal histidine residue, His64(E7), and the proximal histidine residue, His93(F8), in myoglobin (Mb) are important for the function of the protein. For example, the increase in the association rate constant for CO binding at low pH has been suggested to be caused by the protonation of these histidine residues. In order to investigate the influence of protonation on the structure of myoglobin, we determined the crystal structures of sperm whale myoglobin to 2.0 A or better in different states of ligation (MbCO, deoxyMb and metMb) at pH values of 4, 5 and 6. The most dramatic change found at low pH is that His64 swings out of the distal pocket in the MbCO structure at pH 4, opening a direct channel from the solvent to the iron atom. This rotation seems to be facilitated by conformational changes in the CD corner. The benzyl side-chain of Phe46(CD4), which has been suggested to be a critical residue in controlling the rotation of His64, moves away from His64 at pH 4 in the deoxyMb structure, allowing more free rotation of His64. Arg45(CD3) is also important for the dynamics of myoglobin, since it influences the pK(a) of His64 and forms a hydrogen bond lattice that hinders the rotation of His64 at neutral pH. This hydrogen-bond lattice disappears at low pH. Although His64 rotates out of the distal pocket in the MbCO structure at pH 4, leaving more space for the CO ligand, the Fe-C-O angle refines to about 130 degrees, the same as those at pH 5 and 6. In the MbCO structure at pH 4, significant conformational changes appear in the EF corner. The peptide plane between Lys79(EF2) and Gly80(EF3) flips about 150 degrees. The occupancy of this conformation in the MbCO structures increases with decreases in pH. On the proximal side of the heme, the bond between the heme iron atom and N(epsilon) of His93 remains intact under the experimental conditions in the MbCO and deoxyMb structures, but appears elongated in the metMb structure at pH 4, representing either a weakened bond or the breakage of the bond in some fraction of the molecules in the crystal.

Mechanical properties of tropomyosin and implications for muscle regulation

Tropomyosin is a protein that controls the interactions of actin and myosin as a part of the regulation of muscle contraction. The 420 A long alpha-helical coiled-coil molecules form long filaments, both in muscle and in crystals. The x-ray diffraction data from tropomyosin crystals have indicated large scale motions of the filaments that can be related to the inherent mechanical properties of the molecule, and by extension, to the role of tropomyosin in the cooperative activation of the thin filaments of muscle. Diffuse scattering analysis has provided information about the amplitudes of the motions that has been used to calculate the intrinsic flexibility of the molecule. It can then be shown that each tropomyosin molecule by itself can only mediate interactions of the nearest-neighboring tropomyosin molecules along the filament. The repeating nature of the thin filament, however, allows the entire filament to activate cooperatively.

XRayView: a teaching aid for X-ray crystallography

A software package, XRayView, has been developed that uses interactive computer graphics to introduce basic concepts of x-ray diffraction by crystals, including the reciprocal lattice, the Ewald sphere construction, Laue cones, the wavelength dependence of the reciprocal lattice, primitive and centered lattices and systematic extinctions, rotation photography. Laue photography, space group determination and Laue group symmetry, and the alignment of crystals by examination of reciprocal space. XRayView is designed with "user-friendliness" in mind, using pull-down menus to control the program. Many of the experiences of using real x-ray diffraction equipment to examine crystalline diffraction can be simulated. Exercises are available on-line to guide the users through many typical x-ray diffraction experiments.

A double mutant of sperm whale myoglobin mimics the structure and function of elephant myoglobin

The functional, spectral, and structural properties of elephant myoglobin and the L29F/H64Q mutant of sperm whale myoglobin have been compared in detail by conventional kinetic techniques, infrared and resonance Raman spectroscopy, 1H NMR, and x-ray crystallography. There is a striking correspondence between the properties of the naturally occurring elephant protein and those of the sperm whale double mutant, both of which are quite distinct from those of native sperm whale myoglobin and the single H64Q mutant. These results and the recent crystal structure determination by Bisig et al. (Bisig, D. A., Di Iorio, E. E., Diederichs, K., Winterhalter, K. H., and Piontek, K. (1995) J. Biol. Chem. 270, 20754-20762) confirm that a Phe residue is present at position 29 (B10) in elephant myoglobin, and not a Leu residue as is reported in the published amino acid sequence. The single Gln64(E7) substitution lowers oxygen affinity approximately 5-fold and increases the rate of autooxidation 3-fold. These unfavorable effects are reversed by the Phe29(B10) replacement in both elephant myoglobin and the sperm whale double mutant. The latter, genetically engineered protein was originally constructed to be a blood substitute prototype with moderately low O2 affinity, large rate constants, and increased resistance to autooxidation. Thus, the same distal pocket combination that we designed rationally on the basis of proposed mechanisms for ligand binding and autooxidation is also found in nature.

Phe-46(CD4) orients the distal histidine for hydrogen bonding to bound ligands in sperm whale myoglobin

The role of Phe-46(CD4) in modulating the functional properties of sperm whale myoglobin was investigated by replacing this residue with Leu, Ile, Val, Ala, Trp, Tyr, and Glu. This highly conserved amino acid almost makes direct contact with the distal histidine and has been postulated to affect ligand binding. The overall association rate constants for CO, O2, and NO binding were little affected by decreasing the size of residue 46 step-wise from Phe to Leu to Val to Ala. In contrast, the rates of CO, O2, and NO dissociation increased 4-, 10-, and 25-fold, respectively, for the same series of mutants, causing large decreases in the affinity of myoglobin for all three diatomic gases. The rates of autooxidation at 37 degrees C, pH 7.0 increased dramatically from approximately 0.1-0.3 h-1 for wild-type, Tyr-46, and Trp-46 myoglobins to 1.5, 5.2, 4.9, and 5.0 h-1 for the Leu-46, Ile-46, Val-46 and Ala-46 mutants, respectively. Rates of NO and O2 geminate recombination were measured using 35 ps and 9 ns laser excitation pulses. Decreasing the size of residue 46 causes significant decreases in the extent of both picosecond and nanosecond rebinding processes. High resolution structures of Leu-46 and Val-46 metmyoglobins, Val-46 CO-myoglobin, and Val-46 deoxymyoglobin were determined by X-ray crystallography. When Phe-46 is replaced by Val, the loss of internal packing volume is compensated by (1) contraction of the CD corner toward the core of the protein, (2) movement of the E-helix toward the mutation site, (3) greater exposure of the distal pocket to intruding solvent molecules, and (4) large disorder in the position of the side chain of the distal histidine (His-64). In wild-type myoglobin, the van der Waals contact between C zeta of Phe-46 and C beta of His-64 appears to restrict rotation of the imidazole side chain. Insertion of Val at position 46 relieves this steric restriction, allowing the imidazole side chain to rotate about the C alpha - C beta bond toward the surface of the globin and about the C beta - C gamma bond toward the space previously occupied by the native Phe-46 side chain. This movement disrupts hydrogen bonding with bound ligands, causing significant decreases in affinity, and opens the distal pocket to solvent water molecules, causing marked increases in the rate of autooxidation.(ABSTRACT TRUNCATED AT 400 WORDS)

Automatic identification of discrete substates in proteins: singular value decomposition analysis of time-averaged crystallographic refinements

The singular value decomposition (SVD) provides a method for decomposing a molecular dynamics trajectory into fundamental modes of atomic motion. The right singular vectors are projections of the protein conformations onto these modes showing the protein motion in a generalized low-dimensional basis. Statistical analysis of the right singular vectors can be used to classify discrete configurational substates in the protein. The configuration space portraits formed from the right singular vectors can also be used to visualize complex high-dimensional motion and to examine the extent of configuration space sampling by the simulation.

The D-helix in myoglobin and in the beta subunit of hemoglobin is required for the retention of heme

All globins consist of eight helices and interconnecting loops except alpha hemoglobin subunits which lack the D-helix due to deletion of five consecutive residues. Previous site-directed mutagenesis work suggested that this deletion is a neutral modification in hemoglobin with respect to equilibrium O2 binding [Komiyama, N. H., Shih, T.-B., Looker, D., Tame, J., & Nagai, K. (1991) Nature 352, 349-351]. To examine the role of the D-helix in myoglobin, we have measured the O2 and CO binding and hemin dissociation properties of recombinant sperm whale myoglobin mutants in which residues 52-56 have been deleted, Mb(-D), replaced by five alanines, Mb(Ala52-56), and substituted with four alanines and a methionine, Mb(Ala52-55Met56). Crystal structures of aquometMb(-D) and aquometMb(Ala52-55Met56) were determined to 2.0 A resolution and show that the conformation of the distal pocket is little affected by removal of the D-helix or mutations in this region. As a result, these mutations have little effect on O2 and CO binding. Diffuse electron density is observed in the region between the C- and E-helices of Mb(-D), indicating a highly mobile or heterogeneous conformation in this portion of the tertiary structure. This flexibility provides an explanation for the 50-fold higher rate of hemin loss from Mb(-D) as compared to that from wild-type myoglobin. Hemin loss from Mb(Ala52-56) is also rapid. In contrast, Mb(Ala52-55Met56) shows a well-defined D-helix and has a rate of hemin loss identical to that of wild-type holoprotein [corrected].(ABSTRACT TRUNCATED AT 250 WORDS)

A sampling problem in molecular dynamics simulations of macromolecules

Correlations in low-frequency atomic displacements predicted by molecular dynamics simulations on the order of 1 ns are undersampled for the time scales currently accessible by the technique. This is shown with three different representations of the fluctuations in a macromolecule: the reciprocal space of crystallography using diffuse x-ray scattering data, real three-dimensional Cartesian space using covariance matrices of the atomic displacements, and the 3N-dimensional configuration space of the protein using dimensionally reduced projections to visualize the extent to which phase space is sampled.

Visualization of dynamic simulations of muscle thin filaments

Following a novel computational formalism, the thin filament of muscle can be modeled by a computational machine containing a large number of finite automata that have one-to-one correspondence with the constituent protein molecules. Computer graphics can be used to visualize the correspondence between the states of finite automata and the configurations of protein molecules according to the structural data. The dynamic simulation of the muscle filament that corresponds to the concurrent state transitions of finite automata can be represented as a sequence of video images. The kinetic and structural knowledge of individual protein molecules is, therefore, integrated into a coherent and functional system. This type of computation and visualization can also be useful for the investigation of molecular structure, function, and interaction in various complex biological systems.

Structure and dynamics of the water around myoglobin

The interplay between simulations at various levels of hydration and experimental observables has led to a picture of the role of solvent in thermodynamics and dynamics of protein systems. One of the most studied protein-solvent systems is myoglobin, which serves as a paradigm for the development of structure-function relationships in many biophysical studies. We review here some aspects of the solvation of myoglobin and the resulting implications. In particular, recent theoretical and simulation studies unify much of the diverse set of experimental results on water near proteins.

Structural and functional effects of apolar mutations of the distal valine in myoglobin

High-resolution structures of the aquomet, deoxy, and CO forms of Ala68, Ile68, Leu68, and Phe68 sperm whale myoglobins have been determined by X-ray crystallography. These 12 new structures, plus those of wild-type myoglobin, have been used to interpret the effects of mutations at position 68 and the effects of cobalt substitution on the kinetics of O2, CO, and NO binding. Molecular dynamics simulations based on crystal structures have provided information about the time-dependent behavior of photolyzed ligands for comparison with picosecond geminate recombination studies. The Val68-->Ala mutation has little effect on the structure and function of myoglobin. In Ala68 deoxymyoglobin, as in the wild-type protein, a water molecule hydrogen-bonded to the N epsilon atom of the distal histidine restricts ligand binding and appears to be more important in regulating the function of myoglobin than direct steric interactions between the ligand and the C gamma atoms of the native valine side-chain. This distal pocket water molecule is displaced by the larger side-chains at position 68 in the crystal structures of Leu68 and Ile68 deoxymyoglobins. The Leu68 side-chain can rotate about its C alpha-C beta and C beta-C gamma bonds to better accommodate bound ligands, resulting in net increases in overall association rate constants and affinities due to the absence of the distal pocket water molecule. However, the flexibility of Leu68 makes simulation of picosecond NO recombination difficult since multiple starting conformations are possible. In the case of Ile68, rotation of the substituted side-chain is restricted due to branching at the beta carbon, and as a result, the delta methyl group is located close to the iron atom in both the deoxy and liganded structures. The favorable effect of displacing the distal pocket water molecule is offset by direct steric hindrance between the bound ligand and the terminal carbon atom of the isoleucine side-chain, resulting in net decreases in affinity for all three ligands and inhibition of geminate recombination which is reproduced in the molecular dynamics simulations. In Phe68 myoglobin, the benzyl side-chain is pointed away from the ligand binding site, occupying a region in the back of the distal pocket. As in wild-type and Ala68 myoglobins, a well-defined water molecule is found hydrogen bonded to the distal histidine in Phe68 deoxymyoglobin. This water molecule, in combination with the large size of the benzyl side-chain, markedly reduces the speed and extent of ligand movement into the distal pocket. (ABSTRACT TRUNCATED AT 400 WORDS)

Crystal structure of photolysed carbonmonoxy-myoglobin

Myoglobin is a globular haem protein that reversibly binds ligands such as O2 and CO. Single photons of visible light can break the covalent bond between CO and the haem iron in carbon-monoxy-myoglobin (MbCO) and thus form an unstable intermediate, Mb*CO, with the CO inside the protein. The ensuing rebinding process has been extensively studied as a model for the interplay of dynamics, structure and function in protein reactions. We have used X-ray crystallography at liquid-helium temperatures to determine the structure of Mb*CO to a resolution of 1.5 A. The photodissociated CO lies on top of the haem pyrrole ring C. Comparison with the CO-bound and unligated myoglobin structures reveals that on photodissociation of the CO, the haem 'domes', the iron moves partially out of the haem plane, the iron-proximal histidine bonds is compressed, the F helix is strained and the distal histidine swings towards the outside of the ligand-binding pocket.

Nitric oxide recombination to double mutants of myoglobin: role of ligand diffusion in a fluctuating heme pocket

Picosecond recombination of nitric oxide to the double mutants of myoglobin, His64Gly-Val68Ala and His64Gly.Val68Ile, at E7 and E11, has been studied experimentally and by computation. It is shown that distal residues have a profound effect on NO recombination. Recombination in the mutants may be explained in terms of fluctuating free volume and structure of the heme pocket. The double mutants provide insight into the effects of free volume and steric hindrance on rates of ligand rebinding following photolysis. Water molecules of the first solvation shell replace surface residues deleted by mutation and can block apparent holes in the protein structure. Thus, water molecules extend the time required for ligands to escape significantly to a nanosecond time scale, which is much longer than would be expected for an open heme pocket. Both nearly exponential (G64A68) and markedly nonexponential (native and G64I68) kinetics are observed, a result at variance with expectation from the model of Petrich et al. [Petrich, J.W., Lambry, J.C., Kuczera, K., Karplus, M., Poyart, C., & Martin, J.L. (1991) Biochemistry 30, 3975-3987], which attributes nonexponential kinetics to proximal effects.