Radical SAM, a novel protein superfamily linking unresolved steps in familiar biosynthetic pathways with radical mechanisms: functional characterization using new analysis and information visualization methods

A novel protein superfamily with over 600 members was discovered by iterative profile searches and analyzed with powerful bioinformatics and information visualization methods. Evidence exists that these proteins generate a radical species by reductive cleavage of S-adenosylmethionine (SAM) through an unusual Fe-S center. The superfamily (named here Radical SAM) provides evidence that radical-based catalysis is important in a number of previously well- studied but unresolved biochemical pathways and reflects an ancient conserved mechanistic approach to difficult chemistries. Radical SAM proteins catalyze diverse reactions, including unusual methylations, isomerization, sulfur insertion, ring formation, anaerobic oxidation and protein radical formation. They function in DNA precursor, vitamin, cofactor, antibiotic and herbicide biosynthesis and in biodegradation pathways. One eukaryotic member is interferon-inducible and is considered a candidate drug target for osteoporosis; another is observed to bind the neuronal Cdk5 activator protein. Five defining members not previously recognized as homologs are lysine 2,3-aminomutase, biotin synthase, lipoic acid synthase and the activating enzymes for pyruvate formate-lyase and anaerobic ribonucleotide reductase. Two functional predictions for unknown proteins are made based on integrating other data types such as motif, domain, operon and biochemical pathway into an organized view of similarity relationships.

Dynamic properties of the Ras switch I region and its importance for binding to effectors

We have investigated the dynamic properties of the switch I region of the GTP-binding protein Ras by using mutants of Thr-35, an invariant residue necessary for the switch function. Here we show that these mutants, previously used as partial loss-of-function mutations in cell-based assays, have a reduced affinity to Ras effector proteins without Thr-35 being involved in any interaction. The structure of Ras(T35S)⋅GppNHp was determined by x-ray crystallography. Whereas the overall structure is very similar to wildtype, residues from switch I are completely invisible, indicating that the effector loop region is highly mobile. 31P-NMR data had indicated an equilibrium between two rapidly interconverting conformations, one of which (state 2) corresponds to the structure found in the complex with the effectors. 31P-NMR spectra of Ras mutants (T35S) and (T35A) in the GppNHp form show that the equilibrium is shifted such that they occur predominantly in the nonbinding conformation (state 1). On addition of Ras effectors, Ras(T35S) but not Ras(T35A) shift to positions corresponding to the binding conformation. The structural data were correlated with kinetic experiments that show two-step binding reaction of wild-type and (T35S)Ras with effectors requires the existence of a rate-limiting isomerization step, which is not observed with T35A. The results indicate that minor changes in the switch region, such as removing the side chain methyl group of Thr-35, drastically affect dynamic behavior and, in turn, interaction with effectors. The dynamics of the switch I region appear to be responsible for the conservation of this threonine residue in GTP-binding proteins.

cis-Isomers of Cytokinins Predominate in Chickpea Seeds throughout Their Development1

Trans-isomers of cytokinins (CK) are thought to predominate and have greater biological activity than corresponding cis-isomers in higher plants. However, this study demonstrates a system within which the predominant CK are cis-isomers. CK were measured at four developmental stages in developing chickpea (Cicer arietinum L. cultivar Kaniva) seeds by gas chromatography-mass spectrometry. Concentrations were highest at an early endospermic fluid stage and fell considerably when the cotyledons expanded. The cis-isomers of zeatin nucleotide ([9R-MP]Z), zeatin riboside ([9R]Z), and zeatin (Z) were present in greater concentrations than those of corresponding trans-isomers: (trans)[9R-MP]Z, (trans)[9R]Z, (trans)Z, or dihydrozeatin riboside. Dihydrozeatin, dihydrozeatin nucleotide, and the isopentenyl-type CK concentrations were either low or not detectable. Root xylem exudates also contained predominantly cis-isomers of [9R-MP]Z and [9R]Z. Identities of (cis)[9R]Z and (cis)Z were confirmed by comparison of ion ratios and retention indices, and a full spectrum was obtained for (cis)[9R]Z. Tissues were extracted under conditions that minimized the possibility of RNase hydrolysis of tRNA following tissue disruption, being a significant source of the cis-CK. Since no isomerization of (trans)[2H]CK internal standards occurred, it is unlikely that the cis-CK resulted from enzymic or nonenzymic isomerization during extraction. Although quantities of total CK varied, similar CK profiles were found among three different chickpea cultivars and between adequately watered and water-stressed plants. Developing chickpea seeds will be a useful system for investigating the activity of cis-CK or determining the origin and metabolism of free CK.

Dynamical hierarchy in transition states: Why and how does a system climb over the mountain?

How a reacting system climbs through a transition state during the course of a reaction has been an intriguing subject for decades. Here we present and quantify a technique to identify and characterize local invariances about the transition state of an N-particle Hamiltonian system, using Lie canonical perturbation theory combined with microcanonical molecular dynamics simulation. We show that at least three distinct energy regimes of dynamical behavior occur in the region of the transition state, distinguished by the extent of their local dynamical invariance and regularity. Isomerization of a six-atom Lennard–Jones cluster illustrates this: up to energies high enough to make the system manifestly chaotic, approximate invariants of motion associated with a reaction coordinate in phase space imply a many-body dividing hypersurface in phase space that is free of recrossings even in a sea of chaos. The method makes it possible to visualize the stable and unstable invariant manifolds leading to and from the transition state, i.e., the reaction path in phase space, and how this regularity turns to chaos with increasing total energy of the system. This, in turn, illuminates a new type of phase space bottleneck in the region of a transition state that emerges as the total energy and mode coupling increase, which keeps a reacting system increasingly trapped in that region.

The relaxation dynamics of the excited electronic states of retinal in bacteriorhodopsin by two-pump-probe femtosecond studies

We present the results of two-pump and probe femtosecond experiments designed to follow the relaxation dynamics of the lowest excited state (S1) populated by different modes. In the first mode, a direct (S0 → S1) radiative excitation of the ground state is used. In the second mode, an indirect excitation is used where the S1 state is populated by the use of two femtosecond laser pulses with different colors and delay times between them. The first pulse excites the S0 → S1 transition whereas the second pulse excites the S1 → Sn transition. The nonradiative relaxation from the Sn state populates the lowest excited state. Our results suggest that the S1 state relaxes faster when populated nonradiatively from the Sn state than when pumped directly by the S0 → S1 excitation. Additionally, the Sn → S1 nonradiative relaxation time is found to change by varying the delay time between the two pump pulses. The observed dependence of the lowest excited state population as well as its dependence on the delay between the two pump pulses are found to fit a kinetic model in which the Sn state populates a different surface (called S′1) than the one being directly excited (S1). The possible involvement of the Ag type states, the J intermediate, and the conical intersection leading to the S0 or to the isomerization product (K intermediate) are discussed in the framework of the proposed model.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

The multiple phenotypes of allosteric receptor mutants.

Channel-linked neurotransmitter receptors are membrane-bound heterooligomers made up of distinct, although homologous, subunits. They mediate chemo-electrical signal transduction and its regulation via interconversion between multiple conformations that exhibit distinct pharmacological properties and biological activities. The large diversity of functional properties and the widely pleiotropic phenotypes, which arise from point mutations in their subunits (or from subunit substitutions), are interpreted in terms of an allosteric model that incorporates multiple discrete conformational states. The model predicts that three main categories of phenotypes may result from point mutations, altering selectively one (or more) of the following features: (i) the properties of individual binding sites (K phenotype), (ii) the biological activity of the ion channel (gamma phenotype) of individual conformations, or (iii) the isomerization constants between receptor conformations (L phenotype). Several nicotinic acetylcholine and glycine receptor mutants with complex phenotypes are quantitatively analyzed in terms of the model, and the analogies among phenotypes are discussed.

Transcriptional activation by protein-induced DNA bending: evidence for a DNA structural transmission model.

Integration host factor (IHF) is a DNA-bending protein that binds to an upstream activating sequence (UAS1) and, on a negatively supercoiled DNA template, activates transcription from the ilvPG promoter of the ilvG-MEDA operon of Escherichia coli. The transcriptional initiation site of the ilvGMEDA operon is located 92 bp downstream of UAS1. Activation is still observed when the orientation of the upstream IHF binding site is reversed. This manipulation places the IHF binding site on the opposite face of the DNA helix, directs the IHF-induced DNA bend in the opposite direction, and presents the opposite face of the nonsymmetrical, heterodimeric, IHF molecule to the downstream RNA polymerase. Lymphoid enhancer-binding factor, LEF-1, is a DNA-bending, lymphoid-specific, mammalian transcription factor that shares no amino acid sequence similarity with IHF. When the IHF site in UAS1 is replaced with a LEF-1 site, LEF-1 activates transcription from the downstream ilvPG promoter in E. coli as well as it is activated by its natural activator, IHF. These results suggest that specific interactions between IHF and RNA polymerase are not required for activation. The results of DNA structural studies show that IHF forms a protein-DNA complex in the UAS1 region that, in the absence of RNA polymerase, alters the structure of the DNA helix in the -10 hexanucleotide region of the downstream ilvPG promoter. The results of in vitro abortive transcription assays show that IIIF also increases the apparent rate of RNA polymerase isomerization from a closed to an open complex. We suggest, therefore, that IHF activates transcription by forming a higher-order protein-DNA complex in the UAS1 region that structurally alters the DNA helix in a way that facilitates open complex formation at the downstream ilvPG promoter site.

Phosphinate analogs of D-, D-dipeptides: slow-binding inhibition and proteolysis protection of VanX, a D-, D-dipeptidase required for vancomycin resistance in Enterococcus faecium.

VanX is a D-Ala-D-Ala dipeptidase that is essential for vancomycin resistance in Enterococcus faecium. Contrary to most proteases and peptidases, it prefers to hydrolyze the amino substrate but not the related kinetically and thermodynamically more favorable ester substrate D-Ala-D-lactate. The enzymatic activity of VanX was previously found to be inhibited by the phosphinate analogs of the proposed tetrahedral intermediate for hydrolysis of D-Ala-D-Ala. Here we report that such phosphinates are slow-binding inhibitors. D-3-[(1-Aminoethyl)phosphinyl]-D-2-methylpropionic acid I showed a time-dependent onset of inhibition of VanX and a time-dependent return to uninhibited steady-state rates upon dilution of the enzyme/inhibitor mixture. The initial inhibition constant Ki after immediate addition of VanX to phosphinate I to form the E-I complex is 1.5 microM but is then lowered by a relatively slow isomerization step to a second complex, E-I*, with a final K*i of 0.47 microM. This slow-binding inhibition reflects a Km/K*i ratio of 2900:1. The rate constant for the slow dissociation of complex E-I* is 0.24 min-1. A phosphinate analog with an ethyl group replacing what would be the side chain of the second D-alanyl residue in the normal tetrahedral adduct gives a K*i value of 90 nM. Partial proteolysis of VanX reveals two protease-sensitive loop regions that are protected by the intermediate analog phosphinate, indicating that they may be part of the VanX active site.

Molecular mechanism of protein-retinal coupling in bacteriorhodopsin.

Bacteriorhodopsin is a membrane protein that functions as a light-driven proton pump. Each cycle of proton transport is initiated by the light-induced isomerization of retinal from the all-trans to 13-cis configuration and is completed by the protein-driven reisomerization of retinal to the all-trans configuration. Previous studies have shown that replacement of Leu-93, a residue in close proximity to the 13-methyl group of retinal, by alanine, resulted in a 250-fold increase in the time required to complete each photocycle. Here, we show that the kinetic defect in the photocycle of the Leu-93-->Ala mutant occurs at a stage after the completion of proton transport and can be overcome in the presence of strong background illumination. Time-resolved retinal-extraction experiments demonstrate the continued presence of a 13-cis intermediate in the photocycle of the Leu-93-->Ala mutant well after the completion of proton release and uptake. These results indicate that retinal reisomerization is kinetically the rate-limiting step in the photocycle of this mutant and that the slow thermal reisomerization can be bypassed by the absorption of a second photon. The effects observed for the Leu-93-->Ala mutant are not observed upon replacement of any other residue in van der Waals contact with retinal or upon replacement of Leu-93 by valine. We conclude that the contact between Leu-93 and the 13-methyl group of retinal plays a key role in controlling the rate of protein conformational changes associated with retinal reisomerization and return of the protein to the initial state.

Protein crosslinking studies suggest that Rhizobium meliloti C4-dicarboxylic acid transport protein D, a sigma 54-dependent transcriptional activator, interacts with sigma 54 and the beta subunit of RNA polymerase.

Rhizobium meliloti C4-dicarboxylic acid transport protein D (DCTD) activates transcription by a form of RNA polymerase holoenzyme that has sigma 54 as its sigma factor (referred to as E sigma 54). DCTD catalyzes the ATP-dependent isomerization of closed complexes between E sigma 54 and the dctA promoter to transcriptionally productive open complexes. Transcriptional activation probably involves specific protein-protein interactions between DCTD and E sigma 54. Interactions between sigma 54-dependent activators and E sigma 54 are transient, and there has been no report of a biochemical assay for contact between E sigma 54 and any activator to date. Heterobifunctional crosslinking reagents were used to examine protein-protein interactions between the various subunits of E sigma 54 and DCTD. DCTD was crosslinked to Salmonella typhimurium sigma 54 with the crosslinking reagents succinimidyl 4-(N-maleimidomethyl)cyclohexane-1-carboxylate and N-hydroxysulfosuccinimidyl-4-azidobenzoate. Cys-307 of sigma 54 was identified by site-directed mutagenesis as the residue that was crosslinked to DCTD. DCTD was also crosslinked to the beta subunit of Escherichia coli core RNA polymerase with succinimidyl 4-(N-maleimidomethyl)cyclohexane-1-carboxylate, but not with N-hydroxysulfosuccinimidyl-4-azidobenzoate. These data suggest that interactions of DCTD with sigma 54 and the beta subunit may be important for transcriptional activation and offer evidence for interactions between a sigma 54-dependent activator and sigma 54, as well as the beta subunit of RNA polymerase.

Human immunodeficiency virus reverse transcriptase substrate-induced conformational changes and the mechanism of inhibition by nonnucleoside inhibitors.

A combination of transient kinetic and equilibrium titration methods has been used to show that both primer/template and nucleotide binding to human immunodeficiency virus type 1 (HIV-1) reverse transcriptase are two-step processes. In both cases, after initial formation of relatively weakly bound states, isomerization reactions lead to tightly bound states. In the case of deoxynucleotide binding to the reverse transcriptase-primer/template complex, the second step in the interaction is rate-limiting in the overall reaction during processive polymerization. Discrimination against incorrect nucleotides occurs both in the initial weak binding and in the second step but is purely kinetic in the second step (as opposed to thermodynamic in the first step). Nonnucleoside inhibitors have a relatively small effect on nucleotide-binding steps (overall affinity is reduced by a factor of ca. 10), while the affinity of the primer/template duplex is increased by at least a factor of 10. The major effect of nonnucleoside inhibitors is on the chemical step (nucleotide transfer).

Cyclophilin catalyzes protein folding in yeast mitochondria.

Cyclophilins are a family of ubiquitous proteins that are the intracellular target of the immunosuppressant drug cyclosporin A. Although cyclophilins catalyze peptidylprolyl cis-trans isomerization in vitro, it has remained open whether they also perform this function in vivo. Here we show that Cpr3p, a cyclophilin in the matrix of yeast mitochondria, accelerates the refolding of a fusion protein that was synthesized in a reticulocyte lysate and imported into the matrix of isolated yeast mitochondria. The fusion protein consisted of the matrix-targeting sequence of subunit 9 of F1F0-ATPase fused to mouse dihydrofolate reductase. Refolding of the dihydrofolate reductase moiety in the matrix was monitored by acquisition of resistance to proteinase K. The rate of refolding was reduced by a factor of 2-6 by 2.5 microM cyclosporin A. This reduced rate of folding was also observed with mitochondria lacking Cpr3p. In these mitochondria, protein folding was insensitive to cyclosporin A. The rate of protein import was not affected by cyclosporin A or by deletion of Cpr3p.

Evolutionary importance for the membrane enhancement of the production of vitamin D3 in the skin of poikilothermic animals.

The photoproduction of vitamin D in the skin was essential for the evolutionary development of terrestrial vertebrates. During exposure to sunlight, previtamin D3 formed in the skin is isomerized to vitamin D3 (calciol) by a temperature-dependent process. Since early land vertebrates were poikilothermic, the relatively slow conversion of previtamin D3 to vitamin D3 at ambient temperature put them at serious risk for developing vitamin D deficiency, thus leading to a poorly mineralized skeleton that could have ultimately halted further evolutionary development of vertebrates on land. We evaluated the rate of isomerization of previtamin D3 to vitamin D3 in the skin of iguanas and found the isomerization rate was enhanced by 1100% and 1700% at 25 degrees C and 5 degrees C, respectively. It is likely that the membrane entrapment of previtamin D3 in its s-cis,s-cis conformation is responsible for the markedly enhanced conversion of previtamin D3 to vitamin D3. The membrane-enhanced production of vitamin D3 ensures the critical supply of vitamin D3 to poikilothermic animals such as iguanas.

A cyclophilin from the polycentric anaerobic rumen fungus Orpinomyces sp. strain PC-2 is highly homologous to vertebrate cyclophilin B.

A cyclophilin (CyP) purified to homogeneity from the polycentric anaerobic rumen fungus Orpinomyces sp. strain PC-2 had a molecular mass of 20.5 kDa and a pI of 8.1. The protein catalyzed the isomerization of the prolyl peptide bond of N-succinyl-Ala-Ala-(cis,trans)-Pro-Phe p-nitroanilide with a kcat/Km value of 9.3 x 10(6) M-1.s-1 at 10 degrees C and pH 7.8. Cyclosporin A strongly inhibited this peptidylprolyl cis-trans isomerase activity with an IC50 of 19.6 nM. The sequence of the first 30 N-terminal amino acids of this CyP had high homology with the N-terminal sequences of other eukaryotic CyPs. By use of a DNA hybridization probe amplified by PCR with degenerate oligonucleotide primers designed based on the amino acid sequences of the N terminus of this CyP and highly conserved internal regions of other CyPs, a full-length cDNA clone was isolated. It possessed an open reading frame encoding a polypeptide of 203 amino acids with a calculated molecular weight of 21,969, containing a putative hydrophobic signal peptide sequence of 22 amino acids preceding the N terminus of the mature enzyme and a C-terminal sequence, Lys-Ala-Glu-Leu, characteristic of an endoplasmic reticulum retention signal. The Orpinomyces PC-2 CyP is a typical type B CyP. The amino acid sequence of the Orpinomyces CyP exhibits striking degrees of identity with the corresponding human (70%), bovine (69%), mouse (68%), chicken (66%), maize (61%), and yeast (54%) proteins. Phylogenetic analysis based on the CyP sequences indicated that the evolutionary origin of the Orpinomyces CyP was closely related with CyPs of animals.