Dynamical principles in biological processes: A model of charge migration in proteins and DNA

The generalized master equations (GMEs) that contain multiple time scales have been derived quantum mechanically. The GME method has then been applied to a model of charge migration in proteins that invokes the hole hopping between local amino acid sites driven by the torsional motions of the floppy backbones. This model is then applied to analyze the experimental results for sequence-dependent long-range hole transport in DNA reported by Meggers et al. [Meggers, E., Michel-Beyerle, M. E., & Giese, B. (1998) J. Am. Chem. Soc. 120, 12950–12955]. The model has also been applied to analyze the experimental results of femtosecond dynamics of DNA-mediated electron transfer reported by Zewail and co-workers [Wan, C., Fiebig, T., Kelley, S. O., Treadway, C. R., Barton, J. K. & Zewail, A. H. (1999) Proc. Natl. Acad. Sci. USA 96, 6014–6019]. The initial events in the dynamics of protein folding have begun to attract attention. The GME obtained in this paper will be applicable to this problem.

Interactions between tRNA identity nucleotides and their recognition sites in glutaminyl-tRNA synthetase determine the cognate amino acid affinity of the enzyme.

Sequence-specific interactions between aminoacyl-tRNA synthetases and their cognate tRNAs both ensure accurate RNA recognition and prevent the binding of noncognate substrates. Here we show for Escherichia coli glutaminyl-tRNA synthetase (GlnRS; EC 6.1.1.18) that the accuracy of tRNA recognition also determines the efficiency of cognate amino acid recognition. Steady-state kinetics revealed that interactions between tRNA identity nucleotides and their recognition sites in the enzyme modulate the amino acid affinity of GlnRS. Perturbation of any of the protein-RNA interactions through mutation of either component led to considerable changes in glutamine affinity with the most marked effects seen at the discriminator base, the 10:25 base pair, and the anticodon. Reexamination of the identity set of tRNA(Gln) in the light of these results indicates that its constituents can be differentiated based upon biochemical function and their contribution to the apparent Gibbs' free energy of tRNA binding. Interactions with the acceptor stem act as strong determinants of tRNA specificity, with the discriminator base positioning the 3' end. The 10:25 base pair and U35 are apparently the major binding sites to GlnRS, with G36 contributing both to binding and recognition. Furthermore, we show that E. coli tryptophanyl-tRNA synthetase also displays tRNA-dependent changes in tryptophan affinity when charging a noncognate tRNA. The ability of tRNA to optimize amino acid recognition reveals a novel mechanism for maintaining translational fidelity and also provides a strong basis for the coevolution of tRNAs and their cognate synthetases.

Postsynaptic clustering of γ-aminobutyric acid type A receptors by the γ3 subunit in vivo

Synaptic localization of γ-aminobutyric acid type A (GABAA) receptors is a prerequisite for synaptic inhibitory function, but the mechanism by which different receptor subtypes are localized to postsynaptic sites is poorly understood. The γ2 subunit and the postsynaptic clustering protein gephyrin are required for synaptic localization and function of major GABAA receptor subtypes. We now show that transgenic overexpression of the γ3 subunit in γ2 subunit-deficient mice restores benzodiazepine binding sites, benzodiazepine-modulated whole cell currents, and postsynaptic miniature currents, suggesting the formation of functional, postsynaptic receptors. Moreover, the γ3 subunit can substitute for γ2 in the formation of GABAA receptors that are synaptically clustered and colocalized with gephyrin in vivo. These clusters were formed even in brain regions devoid of endogenous γ3 subunit, indicating that the factors present for clustering of γ2 subunit-containing receptors are sufficient to cluster γ3 subunit-containing receptors. The GABAA receptor and gephyrin-clustering properties of the ectopic γ3 subunit were also observed for the endogenous γ3 subunit, but only in the absence of the γ2 subunit, suggesting that the γ3 subunit is at a competitive disadvantage with the γ2 subunit for clustering of postsynaptic GABAA receptors in wild-type mice.

A point mutation in the γ2 subunit of γ-aminobutyric acid type A receptors results in altered benzodiazepine binding site specificity

Benzodiazepines allosterically modulate γ-aminobutyric acid (GABA) evoked chloride currents of γ-aminobutyric acid type A (GABAA) receptors. Coexpression of either rat γ2 or γ3, in combination with α1 and β2 subunits, results both in receptors displaying high [3H]Ro 15-1788 affinity. However, receptors containing a γ3 subunit display a 178-fold reduced affinity to zolpidem as compared with γ2-containing receptors. Eight chimeras between γ2 and γ3 were constructed followed by nine different point mutations in γ2, each to the homologous amino acid residue found in γ3. Chimeric or mutant γ subunits were coexpressed with α1 and β2 in human embryonic kidney 293 cells to localize amino acid residues responsible for the reduced zolpidem affinity. Substitution of a methionine-to-leucine at position 130 of γ2 (γ2M130L) resulted in a 51-fold reduction in zolpidem affinity whereas the affinity to [3H]Ro 15-1788 remained unchanged. The affinity for diazepam was only decreased by about 2-fold. The same mutation resulted in a 9-fold increase in Cl 218872 affinity. A second mutation (γ2M57I) was found to reduce zolpidem affinity by about 4-fold. Wild-type and γ2M130L-containing receptors were functionally expressed in Xenopus oocytes. Upon mutation allosteric coupling between agonist and modulatory sites is preserved. Dose–response curves for zolpidem and for diazepam showed that the zolpidem but not the diazepam apparent affinity is drastically reduced. The apparent GABA affinity is not significantly affected by the γ2M130L mutation. The identified amino acid residues may define part of the benzodiazepine binding pocket of GABAA receptors. As the modulatory site in the GABAA receptor is homologous to the GABA site, and to all agonist sites of related receptors, γ2M130 may either point to a homologous region important for agonist binding in all receptors or define a new region not underlying this principle.

Induction of minisatellite mutation in NIH 3T3 cells by treatment with the tumor promoter okadaic acid

Okadaic acid (OA) is a strong tumor promoter of mouse skin carcinogenesis and also a potent inhibitor of serine/threonine protein phosphatases. OA induces various genetic alterations in cultured cells, such as diphtheria-toxin-resistance mutations, sister chromatid exchange, exclusion of exogenous transforming oncogenes, and gene amplification. The present study revealed that it caused minisatellite mutation (MSM) at a high frequency in NIH 3T3 cells, although no microsatellite mutation was found. Nine of 31 clones (29%) exhibited MSM after 6 days of OA treatment, as opposed to only 1 of 30 clones (3%) without OA exposure. Moreover, NIH 3T3 cells treated with OA acquired tumorigenicity in nude mice, giving rise to 7 tumors within 25 weeks in 20 sites where 3 × 106 cells were injected. In contrast, the same numbers of untreated cells gave rise to only one tumor, and the tumor grew much slower. All of three OA-induced tumors examined manifested the MSM. The findings thus point to a molecular mechanism by which OA could function as a tumor promoter, and also the biological relevance of the induction of MSM in the tumorigenic process by OA.

Both N-terminal myosin-binding and C-terminal actin-binding sites on smooth muscle caldesmon are required for caldesmon-mediated inhibition of actin filament velocity

It has been suggested that the tethering caused by binding of the N-terminal region of smooth muscle caldesmon (CaD) to myosin and its C-terminal region to actin contributes to the inhibition of actin-filament movement over myosin heads in an in vitro motility assay. However, direct evidence for this assumption has been lacking. In this study, analysis of baculovirus-generated N-terminal and C-terminal deletion mutants of chicken-gizzard CaD revealed that the major myosin-binding site on the CaD molecule resides in a 30-amino acid stretch between residues 24 and 53, based on the very low level of binding of CaDΔ24–53 lacking the residues 24–53 to myosin compared with the level of binding of CaDΔ54–85 missing the adjacent residues 54–85 or of the full-length CaD. As expected, deletion of the region between residues 24 and 53 or between residues 54 and 85 had no effect on either actin-binding or inhibition of actomyosin ATPase activity. Deletion of residues 24–53 nearly abolished the ability of CaD to inhibit actin filament velocity in the in vitro motility experiments, whereas CaDΔ54–85 strongly inhibited actin filament velocity in a manner similar to that of full-length CaD. Moreover, CaD1–597, which lacks the major actin-binding site(s), did not inhibit actin-filament velocity despite the presence of the major myosin-binding site. These data provide direct evidence for the inhibition of actin filament velocity in the in vitro motility assay caused by the tethering of myosin to actin through binding of both the CaD N-terminal region to myosin and the C-terminal region to actin.

Identification of active sites in amidase: Evolutionary relationship between amide bond- and peptide bond-cleaving enzymes

Mainly based on various inhibitor studies previously performed, amidases came to be regarded as sulfhydryl enzymes. Not completely satisfied with this generally accepted interpretation, we performed a series of site-directed mutagenesis studies on one particular amidase of Rhodococcus rhodochrous J1 that was involved in its nitrile metabolism. For these experiments, the recombinant amidase was produced as the inclusion body in Escherichia coli to greatly facilitate its recovery and subsequent purification. With regard to the presumptive active site residue Cys203, a Cys203 → Ala mutant enzyme still retained 11.5% of the original specific activity. In sharp contrast, substitutions in certain other positions in the neighborhood of Cys203 had a far more dramatic effect on the amidase. Glutamic acid substitution of Asp191 reduced the specific activity of the mutant enzyme to 1.33% of the wild-type activity. Furthermore, Asp191 → Asn substitution as well as Ser195 → Ala substitution completely abolished the specific activity. It would thus appear that, among various conserved residues residing within the so-called signature sequence common to all amidases, the real active site residues are Asp191 and Ser195 rather than Cys203. Inasmuch as an amide bond (CO-NH2) in the amide substrate is not too far structurally removed from a peptide bond (CO-NH-), the signature sequences of various amidases were compared with the active site sequences of various types of proteases. It was found that aspartic acid and serine residues corresponding to Asp191 and Ser195 of the Rhodococcus amidase are present within the active site sequences of aspartic proteinases, thus suggesting the evolutionary relationship between the two.

Cloning of a gene (RIG-G) associated with retinoic acid-induced differentiation of acute promyelocytic leukemia cells and representing a new member of a family of interferon-stimulated genes

In a cell line (NB4) derived from a patient with acute promyelocytic leukemia, all-trans-retinoic acid (ATRA) and interferon (IFN) induce the expression of a novel gene we call RIG-G (for retinoic acid-induced gene G). This gene codes for a 58-kDa protein containing 490 amino acids with several potential sites for post-translational modification. In untreated NB4 cells, the expression of RIG-G is undetectable. ATRA treatment induces the transcriptional expression of RIG-G relatively late (12–24 hr) in a protein synthesis-dependent manner, whereas IFN-α induces its expression early (30 min to 3 hr). Database search has revealed a high-level homology between RIG-G and several IFN-stimulated genes in human (ISG54K, ISG56K, and IFN-inducible and retinoic acid-inducible 58K gene) and some other species, defining a well conserved gene family. The gene is composed of two exons and has been mapped by fluorescence in situ hybridization to chromosome 10q24, where two other human IFN-stimulated gene members are localized. A synergistic induction of RIG-G expression in NB4 cells by combined treatment with ATRA and IFNs suggests that a collaboration exists between their respective signaling pathways.

A single amino acid substitution converts a carboxylesterase to an organophosphorus hydrolase and confers insecticide resistance on a blowfly

Resistance to organophosphorus (OP) insecticides is associated with decreased carboxylesterase activity in several insect species. It has been proposed that the resistance may be the result of a mutation in a carboxylesterase that simultaneously reduces its carboxylesterase activity and confers an OP hydrolase activity (the “mutant ali-esterase hypothesis”). In the sheep blowfly, Lucilia cuprina, the association is due to a change in a specific esterase isozyme, E3, which, in resistant flies, has a null phenotype on gels stained using standard carboxylesterase substrates. Here we show that an OP-resistant allele of the gene that encodes E3 differs at five amino acid replacement sites from a previously described OP-susceptible allele. Knowledge of the structure of a related enzyme (acetylcholinesterase) suggests that one of these substitutions (Gly137 → Asp) lies within the active site of the enzyme. The occurrence of this substitution is completely correlated with resistance across 15 isogenic strains. In vitro expression of two natural and two synthetic chimeric alleles shows that the Asp137 substitution alone is responsible for both the loss of E3’s carboxylesterase activity and the acquisition of a novel OP hydrolase activity. Modeling of Asp137 in the homologous position in acetylcholinesterase suggests that Asp137 may act as a base to orientate a water molecule in the appropriate position for hydrolysis of the phosphorylated enzyme intermediate.

A bZIP factor, TRAB1, interacts with VP1 and mediates abscisic acid-induced transcription

The transcription factor VP1 regulates maturation and dormancy in plant seeds by activating genes responsive to the stress hormone abscisic acid (ABA). Although activation involves ABA-responsive elements (ABREs), VP1 itself does not specifically bind ABREs. Instead, we have identified and cloned a basic region leucine zipper (bZIP) factor, TRAB1, that interacts with both VP1 and ABREs. Transcription from a chimeric promoter with GAL4-binding sites was ABA-inducible if cells expressed a GAL4 DNA-binding domain∷TRAB1 fusion protein. Results indicate that TRAB1 is a true trans-acting factor involved in ABA-regulated transcription and reveal a molecular mechanism for the VP1-dependent, ABA-inducible transcription that controls maturation and dormancy in plant embryos.

Structure and function in rhodopsin: Rhodopsin mutants with a neutral amino acid at E134 have a partially activated conformation in the dark state*

The Glu-134–Arg-135 residues in rhodopsin, located near the cytoplasmic end of the C helix, are involved in G protein binding, or activation, or both. Furthermore, the charge-neutralizing mutation Glu-134 to Gln-134 produces hyperactivity in the activated state and produces constitutive activity in opsin. The Glu/Asp-Arg charge pair is highly conserved in equivalent positions in other G protein-coupled receptors. To investigate the structural consequences of charge-neutralizing mutations at Glu-134 and Arg-135 in rhodopsin, single spin-labeled side chains were introduced at sites in the cytoplasmic domains of helices C (140), E (227), F (250), or G (316) to serve as “molecular sensors” of the local helix bundle conformation. In each of the spin-labeled rhodopsins, a Gln substitution was introduced at either Glu-134 or Arg-135, and the electron paramagnetic resonance spectrum of the spin label was used to monitor the structural response of the helix bundle. The results indicate that a Gln substitution at Glu-134 induces a photoactivated conformation around helices C and G even in the dark state, an observation of potential relevance to the hyperactivity and constitutive activity of the mutant. In contrast, little change is induced in helix F, which has been shown to undergo a dominant motion upon photoactivation. This result implies that the multiple helix motions accompanying photoactivation are not strongly coupled and can be induced to take place independently. Gln substitution at Arg-135 produces only minor structural changes in the dark- or light-activated conformation, suggesting that this residue is not a determinant of structure in the regions investigated, although it may be functionally important.

Myosin Heavy Chain Phosphorylation Sites Regulate Myosin Localization during Cytokinesis in Live Cells

Conventional myosin II plays a fundamental role in the process of cytokinesis where, in the form of bipolar thick filaments, it is thought to be the molecular motor that generates the force necessary to divide the cell. In Dictyostelium, the formation of thick filaments is regulated by the phosphorylation of three threonine residues in the tail region of the myosin heavy chain. We report here on the effects of this regulation on the localization of myosin in live cells undergoing cytokinesis. We imaged fusion proteins of the green-fluorescent protein with wild-type myosin and with myosins where the three critical threonines had been changed to either alanine or aspartic acid. We provide evidence that thick filament formation is required for the accumulation of myosin in the cleavage furrow and that if thick filaments are overproduced, this accumulation is markedly enhanced. This suggests that myosin localization in dividing cells is regulated by myosin heavy chain phosphorylation.

Specific molecular recognition of mixed nucleic acid sequences: An aromatic dication that binds in the DNA minor groove as a dimer

Phenylamidine cationic groups linked by a furan ring (furamidine) and related compounds bind as monomers to AT sequences of DNA. An unsymmetric derivative (DB293) with one of the phenyl rings of furamidine replaced with a benzimidazole has been found by quantitative footprinting analyses to bind to GC-containing sites on DNA more strongly than to pure AT sequences. NMR structural analysis and surface plasmon resonance binding results clearly demonstrate that DB293 binds in the minor groove at specific GC-containing sequences of DNA in a highly cooperative manner as a stacked dimer. Neither the symmetric bisphenyl nor bisbenzimidazole analogs of DB293 bind significantly to the GC containing sequences. DB293 provides a paradigm for design of compounds for specific recognition of mixed DNA sequences and extends the boundaries for small molecule-DNA recognition.

Selenium redox biochemistry of zinc–sulfur coordination sites in proteins and enzymes

Selenium has been increasingly recognized as an essential element in biology and medicine. Its biochemistry resembles that of sulfur, yet differs from it by virtue of both redox potentials and stabilities of its oxidation states. Selenium can substitute for the more ubiquitous sulfur of cysteine and as such plays an important role in more than a dozen selenoproteins. We have chosen to examine zinc–sulfur centers as possible targets of selenium redox biochemistry. Selenium compounds release zinc from zinc/thiolate-coordination environments, thereby affecting the cellular thiol redox state and the distribution of zinc and likely of other metal ions. Aromatic selenium compounds are excellent spectroscopic probes of the otherwise relatively unstable functional selenium groups. Zinc-coordinated thiolates, e.g., metallothionein (MT), and uncoordinated thiolates, e.g., glutathione, react with benzeneseleninic acid (oxidation state +2), benzeneselenenyl chloride (oxidation state 0) and selenocystamine (oxidation state −1). Benzeneseleninic acid and benzeneselenenyl chloride react very rapidly with MT and titrate substoichiometrically and with a 1:1 stoichiometry, respectively. Selenium compounds also catalyze the release of zinc from MT in peroxidation and thiol/disulfide-interchange reactions. The selenoenzyme glutathione peroxidase catalytically oxidizes MT and releases zinc in the presence of t-butyl hydroperoxide, suggesting that this type of redox chemistry may be employed in biology for the control of metal metabolism. Moreover, selenium compounds are likely targets for zinc/thiolate coordination centers in vivo, because the reactions are only partially suppressed by excess glutathione. This specificity and the potential to undergo catalytic reactions at low concentrations suggests that zinc release is a significant aspect of the therapeutic antioxidant actions of selenium compounds in antiinflammatory and anticarcinogenic agents.

Calcium release from the endoplasmic reticulum of higher plants elicited by the NADP metabolite nicotinic acid adenine dinucleotide phosphate

Higher plants share with animals a responsiveness to the Ca2+ mobilizing agents inositol 1,4,5-trisphosphate (InsP3) and cyclic ADP-ribose (cADPR). In this study, by using a vesicular 45Ca2+ flux assay, we demonstrate that microsomal vesicles from red beet and cauliflower also respond to nicotinic acid adenine dinucleotide phosphate (NAADP), a Ca2+-releasing molecule recently described in marine invertebrates. NAADP potently mobilizes Ca2+ with a K1/2 = 96 nM from microsomes of nonvacuolar origin in red beet. Analysis of sucrose gradient-separated cauliflower microsomes revealed that the NAADP-sensitive Ca2+ pool was derived from the endoplasmic reticulum. This exclusively nonvacuolar location of the NAADP-sensitive Ca2+ pathway distinguishes it from the InsP3- and cADPR-gated pathways. Desensitization experiments revealed that homogenates derived from cauliflower tissue contained low levels of NAADP (125 pmol/mg) and were competent in NAADP synthesis when provided with the substrates NADP and nicotinic acid. NAADP-induced Ca2+ release is insensitive to heparin and 8-NH2-cADPR, specific inhibitors of the InsP3- and cADPR-controlled mechanisms, respectively. However, NAADP-induced Ca2+ release could be blocked by pretreatment with a subthreshold dose of NAADP, as previously observed in sea urchin eggs. Furthermore, the NAADP-gated Ca2+ release pathway is independent of cytosolic free Ca2+ and therefore incapable of operating Ca2+-induced Ca2+ release. In contrast to the sea urchin system, the NAADP-gated Ca2+ release pathway in plants is not blocked by L-type channel antagonists. The existence of multiple Ca2+ mobilization pathways and Ca2+ release sites might contribute to the generation of stimulus-specific Ca2+ signals in plant cells.