New methods of structure refinement for macromolecular structure determination by NMR

Recent advances in multidimensional NMR methodology have permitted solution structures of proteins in excess of 250 residues to be solved. In this paper, we discuss several methods of structure refinement that promise to increase the accuracy of macromolecular structures determined by NMR. These methods include the use of a conformational database potential and direct refinement against three-bond coupling constants, secondary 13C shifts, 1H shifts, T1/T2 ratios, and residual dipolar couplings. The latter two measurements provide long range restraints that are not accessible by other solution NMR parameters.

Endogenous intracellular glutathionyl radicals are generated in neuroblastoma cells under hydrogen peroxide oxidative stress.

We report the detection of endogenous intracellular glutathionyl (GS.) radicals in the intact neuroblastoma cell line NCB-20 under oxidative stress. Spin-trapping and electron paramagnetic resonance (EPR) spectroscopic methods were used for monitoring the radicals. The cells incubated with the spin trap 5,5-dimethyl-1-pyrroline 1-oxide (DMPO) were challenged with H2O2 generated by the enzymic reaction of glucose/glucose oxidase. These cells exhibit the EPR spectrum of the GS. radical adduct of DMPO (DMPO-.SG) without exogenous reduced glutathione (GSH). The identity of this radical adduct was confirmed by observing hyperfine coupling constants identical to previously reported values in in vitro studies, which utilized known enzymic reactions, such as horseradish peroxidase and Cu/Zn superoxide dismutase, with GSH and H2O2 as substrates. The formation of the GS. radicals required viable cells and continuous biosynthesis of GSH. No significant effect on the resonance amplitude by the addition of a membrane-impermeable paramagnetic broadening agent indicated that these radicals were located inside the intact cell. N-Acetyl-L-cysteine (NAC)-treated cells produced NAC-derived free radicals (NAC.) in place of GS. radicals. The time course studies showed that DMPO-.SG formation exhibited a large increase in its concentration after a lag period, whereas DMPO-NAC. formation from NAC-treated cells did not show this sudden increase. These results were discussed in terms of the limit of antioxidant enzyme defenses in cells and the potential role of the GS. radical burst in activation of the transcription nuclear factor NF-kappa B in response to oxidative stress.

Equilibrium constants and free energies in unfolding of proteins in urea solutions

A novel thermodynamic approach to the reversible unfolding of proteins in aqueous urea solutions has been developed based on the premise that urea ligands are bound cooperatively to the macromolecule. When successive stoichiometric binding constants have values larger than expected from statistical effects, an equation for moles of bound urea can be derived that contains imaginary terms. For a very steep unfolding curve, one can then show that the fraction of protein unfolded, B̄, depends on the square of the urea concentration, U, and is given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\bar {B}=\frac{{\mathit{A}}^{{\mathit{2}}}_{{\mathit{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}n\bar {B}}{\mathit{U}}^{{\mathit{2}}}}{{\mathrm{1\hspace{.167em}+\hspace{.167em}}}{\mathit{A}}^{{\mathrm{2}}}_{{\mathrm{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}\bar {B}}{\mathit{U}}^{{\mathrm{2}}}}{\mathrm{.}}\end{equation*}\end{document} A12 is the binding constant as B̄→ 0, and λ is a parameter that reflects the augmentation in affinities of protein for urea as the moles bound increases to the saturation number, n. This equation provides an analytic expression that reproduces the unfolding curve with good precision, suggests a simple linear graphical procedure for evaluating A12 and λ, and leads to the appropriate standard free energy changes. The calculated ΔG° values reflect the coupling of urea binding with unfolding of the protein. Some possible implications of this analysis to protein folding in vivo are described.

Kinetic coupling between electron and proton transfer in cytochrome c oxidase: simultaneous measurements of conductance and absorbance changes.

Bovine heart cytochrome c oxidase is an electron-current driven proton pump. To investigate the mechanism by which this pump operates it is important to study individual electron- and proton-transfer reactions in the enzyme, and key reactions in which they are kinetically and thermodynamically coupled. In this work, we have simultaneously measured absorbance changes associated with electron-transfer reactions and conductance changes associated with protonation reactions following pulsed illumination of the photolabile complex of partly reduced bovine cytochrome c oxidase and carbon monoxide. Following CO dissociation, several kinetic phases in the absorbance changes were observed with time constants ranging from approximately 3 microseconds to several milliseconds, reflecting internal electron-transfer reactions within the enzyme. The data show that the rate of one of these electron-transfer reactions, from cytochrome a3 to a on a millisecond time scale, is controlled by a proton-transfer reaction. These results are discussed in terms of a model in which cytochrome a3 interacts electrostatically with a protonatable group, L, in the vicinity of the binuclear center, in equilibrium with the bulk through a proton-conducting pathway, which determines the rate of proton transfer (and indirectly also of electron transfer). The interaction energy of cytochrome a3 with L was determined independently from the pH dependence of the extent of the millisecond-electron transfer and the number of protons released, as determined from the conductance measurements. The magnitude of the interaction energy, 70 meV (1 eV = 1.602 x 10(-19) J), is consistent with a distance of 5-10 A between cytochrome a3 and L. Based on the recently determined high-resolution x-ray structures of bovine and a bacterial cytochrome c oxidase, possible candidates for L and a physiological role for L are discussed.

Simultaneous determination of helical unwinding angles and intrinsic association constants in ligand–DNA complexes: The interaction between DNA and calichearubicin B

We present a helical unwinding assay for reversibly binding DNA ligands that uses closed circular DNA, topoisomerase I (Topo I), and two-dimensional agarose gel electrophoresis. Serially diluted Topo I relaxation reactions at constant DNA/ligand ratio are performed, and the resulting apparent unwinding of the closed circular DNA is used to calculate both ligand unwinding angle (φ) and intrinsic association constant (Ka). Mathematical treatment of apparent unwinding is formally analogous to that of apparent extinction coefficient data for optical binding titrations. Extrapolation to infinite DNA concentration yields the true unwinding angle of a given ligand and its association constant under Topo I relaxation conditions. Thus this assay delivers simultaneous structural and thermodynamic information describing the ligand–DNA complex. The utility of this assay has been demonstrated by using calichearubicin B (CRB), a synthetic hybrid molecule containing the anthraquinone chromophore of (DA) and the carbohydrate domain of calicheamicin γ1I. The unwinding angle for CRB calculated by this method is −5.3 ± 0.5°. Its Ka value is 0.20 × 106 M−1. For comparison, the unwinding angles of ethidium bromide and DA have been independently calculated, and the results are in agreement with canonical values for these compounds. Although a stronger binder to selected sites, CRB is a less potent unwinder than its parent compound DA. The assay requires only small amounts of ligand and offers an attractive option for analysis of DNA binding by synthetic and natural compounds.

Versatile 5′ phosphoryl coupling of small and large molecules to an RNA

A Ca2+-requiring catalytic RNA is shown to create 5′ phosphate–phosphate linkages with all nucleotides and coenzymes including CoA, nicotinamide adenine dinucleotide phosphate, thiamine phosphate, thiamine pyrophosphate, and flavin mononucleotide. In addition to these small molecules, macromolecules such as RNAs with 5′-diphosphates, and nonnucleotide molecules like Nɛ-phosphate arginine and 6-phosphate gluconic acid also react. That is, the self-capping RNA isolate 6 is an apparently universal 5′ phosphate-linker, reacting with any nucleophile containing an unblocked phosphate. These RNA reactions demonstrate a unique RNA catalytic capability and imply versatile and specific posttranscriptional RNA modification by RNA catalysis.

Translational coupling by modulation of feedback repression in the IF3 operon of Escherichia coli

A pseudoknot formed by a long-range interaction in the mRNA of the initiation factor 3 (IF3) operon is involved in the translational repression of the gene encoding ribosomal protein L35 by another ribosomal protein, L20. The nucleotides forming the 5′ strand of the key stem of the pseudoknot are located within the gene for IF3, whereas those forming the 3′ strand are located 280 nt downstream, immediately upstream of the Shine–Dalgarno sequence of the gene for L35. Here we show that premature termination of IF3 translation at a nonsense codon introduced upstream of the pseudoknot results in a substantial enhancement of L20-mediated repression of L35 expression. Conversely, an increase of IF3 translation decreases repression. These results, in addition to an analysis of the effect of mutations in sequences forming the pseudoknot, indicate that IF3 translation decreases L20-mediated repression of L35 expression. We propose that ribosomes translating IF3 disrupt the pseudoknot and thereby attenuate repression. The result is a novel type of translational coupling, where unfolding of the pseudoknot by ribosomes translating IF3 does not increase expression of L35 directly, but alleviates its repression by L20.

Analysis of exocytosis mutants indicates close coupling between regulated secretion and transcription activation in Tetrahymena

Stimulation of regulated secretory cells promotes protein release via the fusion of cytoplasmic storage vesicles with the plasma membrane. In Tetrahymena thermophila, brief exposure to secretagogue results in synchronous fusion of the entire set of docked dense-core granules with the plasma membrane. We show that stimulation is followed by rapid new dense-core granule synthesis involving gene induction. Two genes encoding granule matrix proteins, GRL1 and GRL4, are shown to undergo induction following stimulation, resulting in ≈10-fold message accumulation within 1 h. The mechanism of induction involves transcriptional regulation, and the upstream region of GRL1 functions in vivo as an inducible promoter in a heterologous reporter construct using the gene encoding green fluorescent protein. Taking advantage of the characterized exocytosis (exo−) mutants available in this system, we asked whether the signals for regranulation were generated directly by the initial stimulation, or whether downstream events were required for transcription activation. Three mutants, with defects at three distinct stages in the regulated secretory pathway, failed to show induction of GRL1 and GRL4 after exposure to secretagogue. These results argue that regranulation depends upon signals generated by the final steps in exocytosis.

Transcriptional coupling between the divergent promoters of a prototypic LysR-type regulatory system, the ilvYC operon of Escherichia coli

The twin-domain model [Liu, L. F. & Wang, J. C. (1987) Proc. Natl. Acad. Sci. USA 84, 7024–7027] suggests that closely spaced, divergent, superhelically sensitive promoters can affect the transcriptional activity of one another by transcriptionally induced negative DNA supercoiling generated in the divergent promoter region. This gene arrangement is observed for many LysR-type-regulated operons in bacteria. We have examined the effects of divergent transcription in the prototypic LysR-type system, the ilvYC operon of Escherichia coli. Double-reporter constructs with the lacZ gene under transcriptional control of the ilvC promoter and the galK gene under control of the divergent ilvY promoter were used to demonstrate that a down-promoter mutation in the ilvY promoter severely decreases in vivo transcription from the ilvC promoter. However, a down-promoter mutation in the ilvC promoter only slightly affects transcription from the ilvY promoter. In vitro transcription assays with DNA topoisomers showed that transcription from the ilvC promoter increases over the entire range of physiological superhelical densities, whereas transcription initiation from the ilvY promoter exhibits a broad optimum at a midphysiological superhelical density. Evidence that this promoter coupling is DNA supercoiling-dependent is provided by the observation that a novobiocin-induced decrease in global negative superhelicity results in an increase in ilvY promoter activity and a decrease in ilvC promoter activity predicted by the in vitro data. We suggest that this transcriptional coupling is important for coordinating basal level expression of the ilvYC operon with the nutritional and environmental conditions of cell growth.

Reducing the flexibility of retinal restores a wild-type-like photocycle in bacteriorhodopsin mutants defective in protein–retinal coupling

The thermal re-isomerization of retinal from the 13-cis to the all-trans state is a key step in the final stages of the photocycle of the light-driven proton pump, bacteriorhodopsin. This step is greatly slowed upon replacement of Leu-93, a residue in van der Waals contact with retinal. The most likely role of this key interaction is that it restricts the flexibility of retinal. To test this hypothesis, we have exchanged native retinal in Leu-93 mutants with bridged retinal analogs that render retinal less flexible by restricting free rotation around either the C10—C11 (9,11-bridged retinal) or C12—C13 (11,13-bridged retinal) single bonds. The effect of the analogs on the photocycle was then determined spectroscopically by taking advantage of the previous finding that the decay of the O intermediate in the Leu-93 mutants provides a convenient marker for retinal re-isomerization. Time-resolved spectroscopic studies showed that both retinal analogs resulted in a dramatic acceleration of the photocycling time by increasing the rate of decay of the O intermediate. In particular, exchange of native retinal in the Leu-93 → Ala mutant with the 9,11-bridged retinal resulted in an acceleration of the decay of the O intermediate to a rate similar to that seen in wild-type bacteriorhodopsin. We conclude that the protein-induced restriction of conformational flexibility in retinal is a key structural requirement for efficient protein–retinal coupling in the bacteriorhodopsin photocycle.

Ab initio calculation of electronic coupling in the photosynthetic reaction center

We have carried out an ab initio electronic structure calculations of electron transfer couplings between chromophores in the bacterial photosynthetic reaction center. The couplings agree remarkably well with parameters obtained from recent quantum dynamical modeling of experimental data assuming an explicit intermediate mechanism. We also have computed couplings on the M-side of the reaction center and have found that the interaction of the primary donor to the M-side intermediate bacteriochlorophyll is quite small because of destructive interference of the two localized coupling matrix elements. This may explain the slow rate of electron transfer down the M-side of the reaction center.

Effective intercellular communication distances are determined by the relative time constants for cyto/chemokine secretion and diffusion

A cell’s ability to effectively communicate with a neighboring cell is essential for tissue function and ultimately for the organism to which it belongs. One important mode of intercellular communication is the release of soluble cyto- and chemokines. Once secreted, these signaling molecules diffuse through the surrounding medium and eventually bind to neighboring cell’s receptors whereby the signal is received. This mode of communication is governed both by physicochemical transport processes and cellular secretion rates, which in turn are determined by genetic and biochemical processes. The characteristics of transport processes have been known for some time, and information on the genetic and biochemical determinants of cellular function is rapidly growing. Simultaneous quantitative analysis of the two is required to systematically evaluate the nature and limitations of intercellular signaling. The present study uses a solitary cell model to estimate effective communication distances over which a single cell can meaningfully propagate a soluble signal. The analysis reveals that: (i) this process is governed by a single, key, dimensionless group that is a ratio of biological parameters and physicochemical determinants; (ii) this ratio has a maximal value; (iii) for realistic values of the parameters contained in this dimensionless group, it is estimated that the domain that a single cell can effectively communicate in is ≈250 μm in size; and (iv) the communication within this domain takes place in 10–30 minutes. These results have fundamental implications for interpretation of organ physiology and for engineering tissue function ex vivo.

Attenuated T2 relaxation by mutual cancellation of dipole–dipole coupling and chemical shift anisotropy indicates an avenue to NMR structures of very large biological macromolecules in solution

Fast transverse relaxation of 1H, 15N, and 13C by dipole-dipole coupling (DD) and chemical shift anisotropy (CSA) modulated by rotational molecular motions has a dominant impact on the size limit for biomacromolecular structures that can be studied by NMR spectroscopy in solution. Transverse relaxation-optimized spectroscopy (TROSY) is an approach for suppression of transverse relaxation in multidimensional NMR experiments, which is based on constructive use of interference between DD coupling and CSA. For example, a TROSY-type two-dimensional 1H,15N-correlation experiment with a uniformly 15N-labeled protein in a DNA complex of molecular mass 17 kDa at a 1H frequency of 750 MHz showed that 15N relaxation during 15N chemical shift evolution and 1HN relaxation during signal acquisition both are significantly reduced by mutual compensation of the DD and CSA interactions. The reduction of the linewidths when compared with a conventional two-dimensional 1H,15N-correlation experiment was 60% and 40%, respectively, and the residual linewidths were 5 Hz for 15N and 15 Hz for 1HN at 4°C. Because the ratio of the DD and CSA relaxation rates is nearly independent of the molecular size, a similar percentagewise reduction of the overall transverse relaxation rates is expected for larger proteins. For a 15N-labeled protein of 150 kDa at 750 MHz and 20°C one predicts residual linewidths of 10 Hz for 15N and 45 Hz for 1HN, and for the corresponding uniformly 15N,2H-labeled protein the residual linewidths are predicted to be smaller than 5 Hz and 15 Hz, respectively. The TROSY principle should benefit a variety of multidimensional solution NMR experiments, especially with future use of yet somewhat higher polarizing magnetic fields than are presently available, and thus largely eliminate one of the key factors that limit work with larger molecules.

The role of members of the pertussis toxin-sensitive family of G proteins in coupling receptors to the activation of the G protein-gated inwardly rectifying potassium channel

Inwardly rectifying potassium (K+) channels gated by G proteins (Kir3.x family) are widely distributed in neuronal, atrial, and endocrine tissues and play key roles in generating late inhibitory postsynaptic potentials, slowing the heart rate and modulating hormone release. They are directly activated by Gβγ subunits released from G protein heterotrimers of the Gi/o family upon appropriate receptor stimulation. Here we examine the role of isoforms of pertussis toxin (PTx)-sensitive G protein α subunits (Giα1–3 and GoαA) in mediating coupling between various receptor systems (A1, α2A, D2S, M4, GABAB1a+2, and GABAB1b+2) and the cloned counterpart of the neuronal channel (Kir3.1+3.2A). The expression of mutant PTx-resistant Gi/oα subunits in PTx-treated HEK293 cells stably expressing Kir3.1+3.2A allows us to selectively investigate that coupling. We find that, for those receptors (A1, α2A) known to interact with all isoforms, Giα1–3 and GoαA can all support a significant degree of coupling to Kir3.1+3.2A. The M4 receptor appears to preferentially couple to Giα2 while another group of receptors (D2S, GABAB1a+2, GABAB1b+2) activates the channel predominantly through Gβγ liberated from GoA heterotrimers. Interestingly, we have also found a distinct difference in G protein coupling between the two splice variants of GABAB1. Our data reveal selective pathways of receptor activation through different Gi/oα isoforms for stimulation of the G protein-gated inwardly rectifying K+ channel.

Proton long-range migration along protein monolayers and its consequences on membrane coupling

It has been shown with lipid layers and more recently with purple membranes that protons have slow surface-to-bulk transfer. This results in long-range proton lateral conduction along membranes. We report here that such lateral transfer can take place along a pure protein film. It is strongly controlled by the packing. Subtle reorganizations of the protein–protein contact can be biological switches between interfacial and delocalized proton pathways between sources and sinks.