My passion for extreme conditions to solve the cytochrome P450 and NO synthase reaction mechanisms

Cytochrome P450, and especially its reaction mechanism, has always been a major research interest of Professor I.C. Gunsalus. The reaction cycle of this enzyme is complex, containing many elementary steps and intermediates, and constitutes a challenge for the scientific community. During his repeated stays in our laboratories in Paris and in Montpellier, he contributed decisively to our study of the P450 reaction mechanism under extreme conditions, i.e., at subzero temperatures and at high pressure. From this initial impulse, we continued the work with different forms of cytochrome P450 and later on with nitric oxide synthase. This paper gives an overview of the insights into these enzymes gained by the use of extreme conditions. These exciting achievements were initiated by numerous discussions with Professor Gunsalus and also reflect a long-lasting collaboration and friendship.     

Science with Gunny

Professor I.C. Gunsalus has not only performed outstanding research for more than 60 years but also is largely responsible for a much closer interaction between physicists and biologists. In this brief paper, I sketch the story of how this interaction came about.     

Global gene expression profiling in Escherichia coli K12: effects of oxygen availability and ArcA

The ArcAB two-component system of Escherichia coli regulates the aerobic/anaerobic expression of genes that encode respiratory proteins whose synthesis is coordinated during aerobic/anaerobic cell growth. A genomic study of E. coli was undertaken to identify other potential targets of oxygen and ArcA regulation. A group of 175 genes generated from this study and our previous study on oxygen regulation (Salmon, K., Hung, S. P., Mekjian, K., Baldi, P., Hatfield, G. W., and Gunsalus, R. P. (2003) J. Biol. Chem. 278, 29837-29855), called our gold standard gene set, have p values <0.00013 and a posterior probability of differential expression value of 0.99. These 175 genes clustered into eight expression patterns and represent genes involved in a large number of cell processes, including small molecule biosynthesis, macromolecular synthesis, and aerobic/anaerobic respiration and fermentation. In addition, 119 of these 175 genes were also identified in our previous study of the fnr allele. A MEME/weight matrix method was used to identify a new putative ArcA-binding site for all genes of the E. coli genome. 16 new sites were identified upstream of genes in our gold standard set. The strict statistical analyses that we have performed on our data allow us to predict that 1139 genes in the E. coli genome are regulated either directly or indirectly by the ArcA protein with a 99% confidence level.     

Nucleotide sequence of the Pseudomonas putida cytochrome P-450cam gene and its expression in Escherichia coli

Cytochrome P-450cam catalyzes the stereospecific methylene hydroxylation of camphor to form 5-exohydroxycamphor and is encoded by the camC gene on the CAM plasmid of Pseudomonas putida, ATCC 17453. The cytochrome P-450cam structural gene has been cloned by mutant complementation in P. putida (Koga, H., Rauchfuss, B., and Gunsalus, I. C. (1985) Biochem. Biophys. Res. Commun. 130, 412-417). We report the complete nucleotide sequence of the camC gene along with 155 base pairs of 5' and 175 base pairs of 3' flanking sequence. Upon comparison of the amino acid sequence derived from the gene sequence to the one obtained from the purified protein (Haniu, M., Armes, L. G., Yasunobu, K. T., Shastry, B. A., and Gunsalus, I. C. (1982) J. Biol. Chem. 257, 12664-12671), five differences were found. The most significant was the addition of a Trp and a Thr residue between Val-54 and Arg-55, thereby increasing the amino acid numbering scheme by 2 after Val-54, bringing the total number of amino acids to 414. Other differences were: Gln-274----Glu-276, Ser-359----His-361, and Asn-405----Asp-407. N-terminal amino acid sequence analysis of the cloned cytochrome P-450cam enzyme expressed in Escherichia coli under the lac promoter showed a faithful translation of the hemo-protein, with the N-terminal Met removed by processing as found in P. putida. Purification to homogeneity of the cloned protein was accomplished by the method used for the CAM plasmid-encoded enzyme of P. putida. The G + C content of the camC gene was found to be 59.0%, caused by a preferred usage of G and C terminated codons. The gene encoding putidaredoxin reductase, camA, was located 22 nucleotides downstream from the cytochrome P-450cam gene. The camA gene initiated with a novel GUG codon, the first such initiator documented in Pseudomonas.     

Book Reviews

Book reviewed in this article:
FOOD MICROSCOPY (1979). Edited by J. C. Vaughan.
A GUIDE TO IDENTIFYING AND CLASSIFYING YEASTS (1979) J. A. Barnett, R. W. Payne and D. Yarrow.
CURRENT PROCEDURES IN CLINICAL BACTERIOLOGY (1978). P. D. Ellner.
YEAR BOOK OF PATHOLOGY AND CLINICAL PATHOLOGY (1979). Edited by F. A. Carone and R. B. Conn.
METHODS IN MICROBIOLOGY Volume 12 (1979). Edited by T. Bergan and J. R. Norris.
THE BACTERIA—A TREATISE ON STRUCTURE AND FUNCTION, VOLUME VII (1979). MECHANISMS OF ADAPTATION. Editor-in-Chief, L. C. Gunsalus.
PROCEEDINGS IN LIFE SCIENCES—MICROBIAL ECOLOGY (1978). Proceedings of the 1st International Microbial Ecology Symposium, 1977. Edited by M. W. Loutit and J. A. R. Miles.     

Isocyanide binding kinetics to monomeric hemoproteins. A study on the ligand partition between solvent and heme pocket.

The kinetics of methyl-, ethyl-, iso-propyl-, and ter-butyl-isocyanide binding to Aplysia limacina myoglobin (distal His----Lys) and the isolated beta chains from hemoglobin Zurich (distal His----Arg) have been investigated by flash photolysis at various temperatures above 0 degrees C. Sperm whale (Physter catodon) myoglobin and the isolated beta chains from normal adult hemoglobin have been used as references. In most reaction systems investigated the apparent extent of photolysis increases with temperature. For sperm whale myoglobin and the normal beta chains the increase is of the same magnitude and not correlated to the type of ligand used. On the contrary, for the two proteins lacking the distal histidine, the phenomenon is dependent on the size of the alkyl side chain of the ligand. The results, analyzed on the basis of the multibarrier model (Austin, R.H., K.W. Beeson, L. Eisenstein, H. Frauenfelder, and I.C. Gunsalus, 1975, Biochemistry, 16:5355-5373), suggest that the partition of the ligand molecules between the solvent and the heme pocket, occurring during the photolysis process, is primarily determined by interactions between the ligand and residues in the heme cavity rather than by diffusion through the protein matrix.     

Nucleotide sequence analysis of the Pseudomonas putida PpG7 salicylate hydroxylase gene (nahG) and its 3'-flanking region

Gene nahG of naphthalene/salicylate catabolic plasmid NAH7 encodes a protein of molecular weight 45,000, salicylate hydroxylase. This enzyme catalyzes the formation of catechol from salicylate, a key intermediate in naphthalene catabolism. DNA sequence analysis of the 3.1-kilobase HindIII fragment containing the nahG locus reveals an open reading frame (ORF) of 1305 base pairs that corresponds to a protein of 434 amino acid residues. The predicted amino acid sequence of salicylate hydroxylase is in agreement with the molecular weight, NH2-terminal amino acid sequence, and total amino acid composition of the purified salicylate hydroxylase [You, I.-S., Murray, R. I., Jollie, D., & Gunsalus, I. C. (1990) Biochem. Biophys. Res. Commun. 169, 1049-1054]. The amino acid sequence between positions 8 and 37 of salicylate hydroxylase shows homology with known ADP binding sites of other FAD-containing oxidoreductases, thus confirming its biochemical function. The sequence of the Pseudomonas putida salicylate hydroxylase was compared with those of other similar flavoproteins. A small DNA segment (831 base pairs) disrupts the continuity of the known gene order nahG and nahH, the latter encoding catechol 2,3-dioxygenase. The complete nucleotide sequence of the intergenic region spanning genes nahG and nahH has been determined and its biological role proposed.     

Mass-weighted molecular dynamics simulation of the protein-ligand complex of rhizopuspepsin and inhibitor.

The mass-weighted molecular dynamics simulation method was developed previously for sampling the multidimensional conformational space of linear and cyclic polypeptides and studying their conformational flexibility. Herein results from molecular dynamics simulations of the protein-ligand complex of the aspartyl protease rhizopuspepsin and a polypeptide inhibitor are reported. The dihedral conformational space sampling for the linear peptide inhibitor in situ was found to be increased in the mass-weighted simulation as in other molecular systems previously studied. More significantly, the physical space of the enzyme binding pocket was also sampled efficiently in the simulations and multiple binding sites were identified for the inhibitor. These results suggest that it may be possible now to study, by computer simulations, the putative initial enzyme-inhibitor complex suggested experimentally from the time-dependent kinetics of enzyme inhibition by slow-binding inhibitors (Morrison, J. F., and C. T. Walsh. 1988. Adv. Enzymol. 61:201), and/or conformational substates in protein-ligand complexes suggested in the study of reassociation dynamics of myoglobin and carbon monoxide following photolysis (Austin, R. H., K. W. Beeson, L. Eisenstein, H. Frauenfelder, and I. C. Gunsalus. 1975. Biochemistry. 14:5355). Moreover, the intermediate binding steps and the molecular flexibility of the inhibitor shown in the MWMD simulation may have crucial roles in the ligand binding process.     

Genes for two herbicide-inducible cytochromes P-450 from Streptomyces griseolus.

Streptomyces griseolus ATCC 11796 contains two inducible, herbicide-metabolizing cytochromes P-450 previously designated P-450SU1 and P-450SU2 (P-450CVA1 and P-450CVB1, respectively, using nomenclature of Nebert et al. [D. W. Nebert, M. Adesnik, M. J. Coon, R. W. Estabrook, F. J. Gonzalez, F. P. Guengerich, I. C. Gunsalus, E. F. Johnson, B. Kemper, W. Levin, I. R. Phillips, R. Sato, and M. R. Waterman, DNA 6:1-11, 1987]). Using antibodies directed against cytochrome P-450SU1, its N-terminal amino acid sequence, and amino acid composition, we cloned the suaC gene encoding cytochrome P-450SU1. Similar information about the cytochrome P-450SU2 protein confirmed that a gene cloned by cross-hybridization to the suaC gene was the subC gene encoding cytochrome P-450SU2. The suaC and subC genes were expressed in Escherichia coli, DNA for both genes was sequenced, and the deduced amino acid sequences were compared with that of the well-characterized cytochrome P-450CAM from Pseudomonas putida. Both cytochromes P-450SU1 and P-450SU2 contain several regions of strong similarity with the amino acid sequence of P-450CAM, primarily in regions of the protein responsible for attachment and coordination of the heme prosthetic group.     

Nickel,a component of factor F 430 from Methanobacterium thermoautotrophicum

Methanobacterium thermoautotrophicum, growing on medium supplemented with 2 mol ⁶³NiCl₂/l, was found to take up 1.2 mol ⁶³Ni per g cells (dry weight). More than 70% of the radioisotope was incorporated into a compound, which dissociated from the protein fraction after heat treatment, was soluble in 70% acetone, and could be purified by chromatography on QAE-Sephadex A-25, Sephadex G-25, and DEAE cellulose. The purified ⁶³Ni labelled compound had an absorption spectrum and properties identical to those of factor F ₄₃₀ and is therefore considered to be identical with factor F ₄₃₀.Factor F ₄₃₀, a compound of molecular weight higher than 1000 with an absorbance maximum at 430 nm, has recently been purified from Methanobacterium thermoautotrophicum (Gunsalus and Wolfe, 1978). The structure and function of this compound are not yet known.     

Zinc-substituted cytochrome P450cam: characterization of protein conformers F420 and F450 by photoinduced electron transfer

Metal substitution of heme proteins is widely applied in the study of biologically relevant electron transfer (ET) reactions. It has been shown that many modified proteins remain in their native conformation and can provide useful insights into the molecular mechanism of electron transfer between the native protein and its substrates. We investigated ET reactions between zinc-substituted cytochrome P450(cam) and small organic compounds such as quinones and ferrocene, which are capable of accessing the protein's hydrophobic channel and binding close to the active site, like its native substrate, camphor. Following the substitution method developed by Gunsalus and co-workers [Wagner, G. C., et al. (1981) J. Biol. Chem. 256, 6262-6265], we have identified two dominant forms of the zinc-substituted protein, F450 and F420, that exhibit different photophysical and photochemical properties. The ET behavior of F420 suggests that hydrophobic redox-active ligands are able to penetrate the hydrophobic channel and place themselves in the direct vicinity of the Zn-porphyrin. In contrast, the slower ET quenching rates observed in the case of F450 indicate that the association is weak and occurs outside of the protein channel. Therefore, we conclude that F420 corresponds to the open structure of the native cytochrome P450(cam) while F450 has a closed or partially closed channel that is characteristic of the camphor-containing cytochrome P450(cam). The existence of two distinct conformers of Zn-bound P450(cam) is consistent with the findings of Goodin and co-workers [Lee, Y.-T., et al. (2010) Biochemistry 49, 3412-3419] and has significant consequences for future electron transfer studies on this popular metalloenzyme.     

Protein dynamics. Comparative investigation on heme-proteins with different physiological roles.

We report the low temperature carbon monoxide recombination kinetics after photolysis and the temperature dependence of the visible absorption spectra of the isolated alpha SH-CO and beta SH-CO subunits from human hemoglobin A in ethylene glycol/water and in glycerol/water mixtures. Kinetic measurements on sperm whale (Physeter catodon) myoglobin and previously published optical spectroscopy data on the latter protein and on human hemoglobin A, in both solvents, (Cordone, L., A. Cupane, M. Leone, E. Vitrano, and D. Bulone. 1988. J. Mol. Biol. 199:312-218) are taken as reference. Low temperature flash photolysis data are analyzed within the multiple substates model proposed by Frauenfelder and co-workers (Austin, R. H., K. W. Beeson, L. Eisenstein, H. Frauenfelder, and I. C. Gunsalus. 1975. Biochemistry. 14:5355-5373). Within this model a distribution of activation enthalpies for ligand binding accounts for the structural heterogeneity of the protein, while the preexponential factor, containing also the entropic contribution to the free energy of the process, is considered to be constant for all conformational substates. Optical spectra are deconvoluted in gaussian components and the temperature dependence of the moments of the resulting bands is analyzed, within the harmonic Frank-Condon approximation, to obtain information on the stereodynamic properties of the heme pocket. The kinetic and spectral parameters thus obtained are found to be protein dependent also with respect to their sensitivity to changes in the composition of the external medium. A close correlation between the kinetic and spectral features is observed for the proteins examined under all experimental conditions studied. The results reported are discussed in terms of differences in the heme pocket structure and in the conformational heterogeneity among the various proteins, as related to their different capability to accommodate constraints imposed by the external medium.     

They are not too excited: the possible role of adult-born neurons in epilepsy

In a new study in this issue of Neuron, Jakubs and colleagues report that adult-generated hippocampal granule cells develop particular functional properties when their birth is induced by epileptic seizures. The new neurons showed reduced excitatory synaptic input and decreased excitability. Their functional integration was thus adjusted to the prevailing functional state in the network. By this means, adult neurogenesis might contribute to network homeostasis in the epileptic temporal lobe.     

Development of polymer latex particles for selective cleavage of mismatched DNA and their application for DNA diagnosis

We developed functional polymer latex particles that can catch and cleave mismatched DNA selectively and propose a new mismatch detection system using the functional particles. The aimed particles possess two functional units composed of mismatch binding protein (MutS) and an anthraquinone derivative (AQ), a light-activated agent that photocleaves dsDNA. Use of the functional particles made it possible to discriminate complementary and mismatched DNAs and photocleave mismatched DNA selectively. The efficiency of photocleavage of mismatched DNA by the functional particles increased with UV irradiation time. It was also found that the functional particles were reusable and had dissociation constants (K(d)) of 1000 and 68.5 nM for G/C homoduplex and G/T heteroduplex, respectively. Using the functional particles and a dsDNA-binding fluorescent dye, SYBR-Gold, we could construct the system for detection of mismatched DNA that was 40 base pairs. The functional particles prepared in this study will be an absolutely new tool for mismatch detection in DNA diagnosis.     

miRDB: a microRNA target prediction and functional annotation database with a wiki interface

MicroRNAs (miRNAs) are short noncoding RNAs that are involved in the regulation of thousands of gene targets. Recent studies indicate that miRNAs are likely to be master regulators of many important biological processes. Due to their functional importance, miRNAs are under intense study at present, and many studies have been published in recent years on miRNA functional characterization. The rapid accumulation of miRNA knowledge makes it challenging to properly organize and present miRNA function data. Although several miRNA functional databases have been developed recently, this remains a major bioinformatics challenge to miRNA research community. Here, we describe a new online database system, miRDB, on miRNA target prediction and functional annotation. Flexible web search interface was developed for the retrieval of target prediction results, which were generated with a new bioinformatics algorithm we developed recently. Unlike most other miRNA databases, miRNA functional annotations in miRDB are presented with a primary focus on mature miRNAs, which are the functional carriers of miRNA-mediated gene expression regulation. In addition, a wiki editing interface was established to allow anyone with Internet access to make contributions on miRNA functional annotation. This is a new attempt to develop an interactive community-annotated miRNA functional catalog. All data stored in miRDB are freely accessible at http://mirdb.org.     

Remembering: functional organization of the declarative memory system

How do brain systems support our subjective experience of recollection and our senses of familiarity and novelty? A new functional imaging study concludes that each of these functions is accomplished by a distinct component of the medial temporal lobe, shedding new light on the functional organization of this memory system.