Dynamic properties of the G93A mutant of copper-zinc superoxide dismutase as detected by NMR spectroscopy: implications for the pathology of familial amyotrophic lateral sclerosis

The backbone assignment of the copper-zinc superoxide dismutase amyotrophic lateral sclerosis G93A mutant was performed on an (15)N-enriched protein sample. (15)N R(1), R(2), and R(1)(rho) and (15)N-(1)H NOE experiments were then carried out at 600 MHz on G93A Cu(2)Zn(2)SOD and the values compared to the dynamics data for the "wild-type" protein. In addition, (15)N and (1)H chemical shift comparisons between wild-type Cu(2)Zn(2)SOD and its G93A mutant were also made. G93A exhibits a higher mobility than wild-type Cu(2)Zn(2)SOD, particularly in loops III and V, on a time scale faster than the rate of protein tumbling. There are also distinct chemical shift and NOE differences in residues 35-42 and 92-95, which comprise these loops. These two regions of Cu(2)Zn(2)SOD form the end of the beta-barrel termed the "beta-barrel plug" [Tainer, J. A., Getzoff, E. D., Beem, K. M., Richardson, J. S., and Richardson, D. C. (1982) J. Mol. Biol. 160, 181-217]. The increased mobility and reduction of the number of observed NOEs in this region indicate an opening of the beta-barrel that may lead to amyloid fibrillogenesis. Alternatively, a motor neuron-specific substrate may bind this region of the protein, leading to deleterious modifications and/or reactions.     

Characterizing a partially folded intermediate of the villin headpiece domain under non-denaturing conditions: contribution of His41 to the pH-dependent stability of the N-terminal subdomain

The contribution of interactions involving the imidazole ring of His41 to the pH-dependent stability of the villin headpiece (HP67) N-terminal subdomain has been investigated by nuclear magnetic resonance (NMR) spin relaxation. NMR-derived backbone N-H order parameters (S2) for wild-type (WT) HP67 and H41Y HP67 indicate that reduced conformational flexibility of the N-terminal subdomain in WT HP67 is due to intramolecular interactions with the His41 imidazole ring. These interactions, together with desolvation effects, contribute to significantly depress the pKa of the buried imidazole ring in the native state. 15N R1rho relaxation dispersion data indicate that WT HP67 populates a partially folded intermediate state that is 10.9 kJ mol(-1) higher in free energy than the native state under non-denaturing conditions at neutral pH. The partially folded intermediate is characterized as having an unfolded N-terminal subdomain while the C-terminal subdomain retains a native-like fold. Although the majority of the residues in the N-terminal subdomain sample a random-coil distribution of conformations, deviations of backbone amide 1H and 15N chemical shifts from canonical random-coil values for residues within 5A of the His41 imidazole ring indicate that a significant degree of residual structure is maintained in the partially folded ensemble. The pH-dependence of exchange broadening is consistent with a linear three-state exchange model whereby unfolding of the N-terminal subdomain is coupled to titration of His41 in the partially folded intermediate with a pKa,I=5.69+/-0.07. Although maintenance of residual interactions with the imidazole ring in the unfolded N-terminal subdomain appears to reduce pKa,I compared to model histidine compounds, protonation of His41 disrupts these interactions and reduces the difference in free energy between the native state and partially folded intermediate under acidic conditions. In addition, chemical shift changes for residues Lys70-Phe76 in the C-terminal subdomain suggest that the HP67 actin binding site is disrupted upon unfolding of the N-terminal subdomain, providing a potential mechanism for regulating the villin-dependent bundling of actin filaments.     

Surface interactions in the complex between cytochrome f and the E43Q/D44N and E59K/E60Q plastocyanin double mutants as determined by 1H-NMR chemical shift analysis

A combination of site-directed mutagenesis and NMR chemical shift perturbation analysis of backbone and side-chain protons has been used to characterize the transient complex of the photosynthetic redox proteins plastocyanin and cytochrome f. To elucidate the importance of charged residues on complex formation, the complex of cytochrome f and E43Q/D44N or E59K/E60Q spinach plastocyanin double mutants was studied by full analysis of the (1)H chemical shifts by use of two-dimensional homonuclear NMR spectra. Both mutants show a significant overall decrease in chemical shift perturbations compared with wild-type plastocyanin, in agreement with a large decrease in binding affinity. Qualitatively, the E43Q/D44N mutant showed a similar interaction surface as wild-type plastocyanin. The interaction surface in the E59K/E60Q mutant was distinctly different from wild type. It is concluded that all four charged residues contribute to the affinity and that residues E59 and E60 have an additional role in fine tuning the orientation of the proteins in the complex.     

15N-labeled tRNA. Identification of 4-thiouridine in Escherichia coli tRNASer1 and tRNATyr2 by 1H-15N two-dimensional NMR spectroscopy

Uridine is uniquely conserved at position 8 in elongator tRNAs and binds to A14 to form a reversed Hoogsteen base pair which folds the dihydrouridine loop back into the core of the L-shaped molecule. On the basis of 1H NMR studies, Hurd and co-workers (Hurd, R. E., Robillard, G. T., and Reid, B. R. (1977) Biochemistry 16, 2095-2100) concluded that the interaction between positions 8 and 14 is absent in Escherichia coli tRNAs with only 3 base pairs in the dihydrouridine stem. We have taken advantage of the unique 15N chemical shift of N3 in thiouridine to identify 1H and 15N resonances for the imino units of S4U8 and s4U9 in E. coli tRNASer1 and tRNATyr2. Model studies with chloroform-soluble derivatives of uridine and 4-thiouridine show that the chemical shifts of the protons in the imino moieties move downfield from 7.9 to 14.4 ppm and from 9.1 to 15.7 ppm, respectively; whereas, the corresponding 15N chemical shifts move downfield from 157.5 to 162.5 ppm and from 175.5 to 180.1 ppm upon hydrogen bonding to 5'-O-acetyl-2',3'-isopropylidene adenosine. The large difference in 15N chemical shifts for U and s4U allows one to unambiguously identify s4U imino resonances by 15N NMR spectroscopy. E. coli tRNASer1 and tRNATyr2 were selectively enriched with 15N at N3 of all uridines and modified uridines. Two-dimensional 1H-15N chemical shift correlation NMR spectroscopy revealed that both tRNAs have resonances with 1H and 15N chemical shifts characteristic of s4UA pairs. The 1H shift is approximately 1 ppm upfield from the typical s4U8 resonance at 14.8 ppm, presumably as a result of local diamagnetic anisotropies. An additional s4U resonance with 1H and 15N shifts typical of interaction of a bound water or a sugar hydroxyl group with s4U9 was discovered in the spectrum of tRNATyr2. Our NMR results for tRNAs with 3-base pair dihydrouridine stems suggest that these molecules have an U8A14 tertiary interaction similar to that found in tRNAs with 4-base pair dihydrouridine stems.     

Protein-protein interaction revealed by NMR T(2) relaxation experiments: the lipoyl domain and E1 component of the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus

T(2) relaxation experiments in combination with chemical shift and site-directed mutagenesis data were used to identify sites involved in weak but specific protein-protein interactions in the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus. The pyruvate decarboxylase component, a heterotetramer E1(alpha(2)beta(2)), is responsible for the first committed and irreversible catalytic step. The accompanying reductive acetylation of the lipoyl group attached to the dihydrolipoyl acetyltransferase (E2) component involves weak, transient but specific interactions between E1 and the lipoyl domain of the E2 polypeptide chain. The interactions between the free lipoyl domain (9 kDa) and free E1alpha (41 kDa), E1beta (35 kDa) and intact E1alpha(2)beta(2) (152 kDa) components, all the products of genes or sub-genes over-expressed in Escherichia coli, were investigated using heteronuclear 2D NMR spectroscopy. The experiments were conducted with uniformly (15)N-labeled lipoyl domain and unlabeled E1 components. Major contact points on the lipoyl domain were identified from changes in the backbone (15)N spin-spin relaxation time in the presence and absence of E1(alpha(2)beta(2)) or its individual E1alpha or E1beta components. Although the E1alpha subunit houses the sequence motif associated with the essential cofactor, thiamin diphosphate, recognition of the lipoyl domain was distributed over sites in both E1alpha and E1beta. A single point mutation (N40A) on the lipoyl domain significantly reduces its ability to be reductively acetylated by the cognate E1. None the less, the N40A mutant domain appears to interact with E1 similarly to the wild-type domain. This suggests that the lipoyl group of the N40A lipoyl domain is not being presented to E1 in the correct orientation, owing perhaps to slight perturbations in the lipoyl domain structure, especially in the lipoyl-lysine beta-turn region, as indicated by chemical shift data. Interaction with E1 and subsequent reductive acetylation are not necessarily coupled.     

Dynamic conformations compared for IgE and IgG1 in solution and bound to receptors.

Dynamic conformations of two distinct immunoglobulin (Ig) isotypes, murine IgE and human IgG1, were examined with fluorescence resonance energy transfer measurements. The IgE mutant epsilon/C gamma 3* and the IgG1 mutant gamma/C gamma 3* each bind [5-(dimethylamino)naphthalen-1-yl]sulfonyl (DNS) in two identical antigen binding sites at the amino (N)-terminal ends of the Ig in the Fab segments. Eosin-DNS bound in these Fab sites served as the acceptor probe in these studies. Both Ig have a carboxy (C)-terminal domain (C gamma 3*) which contains genetically introduced cysteine residues. Modification of these cysteine sulfhydryls with fluorescein maleimide provided donor probes near the C-terminal ends of the Ig in the Fc segment. Energy transfer between the C-terminal and N-terminal ends was compared for these two Ig in solution and when they were found to their respective high-affinity receptors on plasma membranes: IgE-Fc epsilon RI on RBL cell membranes and IgG1-Fc gamma RI on U937 cell membranes. Previous energy-transfer measurements with these probes yielded an average end-to-end distance of 71 A for IgE in solution and 69 A for IgE bound to Fc epsilon RI, indicating that in both situations IgE is bent such that the axes of the Fab segments and the axis of the Fc segment do not form a planar Y-shape [Zheng, Shopes, Holowka, & Baird (1991) Biochemistry 30, 9125]. In the current study we found the average end-to-end distance for IgG1 in solution is 75 A and greater than or equal to 85 A for IgG1 bound to Fc gamma RI, suggesting an average bend conformation for IgG1 as well. The contributions of segmental flexibility to the average distances were assessed directly by measuring the efficiency of energy transfer as a function of variations in donor quantum yield caused by a collisional quencher and using these data to extract a Gaussian distribution of end-to-end distances. The distribution average (rho) and half-width (hw) were determined to be as follows: rho = 75 A, hw = 24 A for IgE in solution; rho = 71 A, hw = 12 A for IgE bound to Fc epsilon RI; and rho = 100 A, hw = 88 A for IgG in solution.(ABSTRACT TRUNCATED AT 400 WORDS)     

Kinetic control of thiamin diphosphate activation in enzymes studied by proton-nitrogen correlated NMR spectroscopy

Proton-nitrogen correlated NMR studies were performed on thiamin diphosphate, which has been specifically labeled with (15)N at the 4'-amino group. After reconstitution of the labeled coenzyme with the apoenzymes of both wild-type pyruvate decarboxylase from Zymomonas mobilis and the E50Q variant, a high-field shift of the (15)N signal of approximately 4 ppm is observed at pH 5.9 when compared to that of the free coenzyme, indicating a higher electron density at the 4'-amino nitrogen in the enzyme-bound state. The pH dependence of the chemical shift of the (15)N signals in the (1)H-(15)N heteronuclear single-quantum coherence NMR spectra reveals typical titration curves for the free as well as the reconstituted coenzyme with nearly identical chemical shift end points. The midpoints of the transitions are at pH 5.3 and 5.0 for the free and enzyme-bound coenzyme, respectively. We conclude that the tremendous rate acceleration of C2-H deprotonation in ThDP enzymes is mainly the result of the enforced V conformation of the cofactor in the active site being perfectly suited to allowing intramolecular acid-base catalysis.     

Influence of N-H...O and O-H...O hydrogen bonds on the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen: a density functional theory study

A computational investigation was carried out to characterize the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen. We found that N-H...O and O-H...O hydrogen bonds around the acetaminophen molecule in the crystal lattice have different influences on the calculated (17)O, (15)N and (13)C chemical shielding eigenvalues and their orientations in the molecular frame of axes. The calculations were performed with the B3LYP method and 6-311++G(d, p) and 6-311+G(d) standard basis sets using the Gaussian 98 suite of programs. Calculated chemical shielding tensors were used to evaluate the (17)O, (15)N, and (13)C NMR chemical shift tensors in crystalline acetaminophen, which are in reasonable agreement with available experimental data. The difference between the calculated NMR parameters of the monomer and molecular clusters shows how much hydrogen-bonding interactions affect the chemical shielding tensors of each nucleus. The computed (17)O chemical shielding tensor on O(1), which is involved in two intermolecular hydrogen bonds, shows remarkable sensitivity toward the choice of the cluster model, whereas the (17)O chemical shielding tensor on O(2) involved in one N-H...O hydrogen bond, shows smaller improvement toward the hydrogen-bonding interactions. Also, a reasonably good agreement between the experimentally obtained solid-state (15)N and (13)C NMR chemical shifts and B3LYP/6-311++G(d, p) calculations is achievable only in molecular cluster model where a complete hydrogen-bonding network is considered. Moreover, at the B3LYP/6-311++G(d, p) level of theory, the calculated (17)O, (15)N and (13)C chemical shielding tensor orientations are able to reproduce the experimental values to a reasonably good degree of accuracy.     

The ATP binding site on rho protein. Affinity labeling of Lys181 by pyridoxal 5'-diphospho-5'-adenosine

We have labeled the nucleoside triphosphate-binding domain of Escherichia coli rho factor with the ATP affinity analog [3H]pyridoxal 5'-diphospho-5'-adenosine (PLP-AMP). PLP-AMP completely inactivates the RNA-dependent ATPase activity of rho upon incorporation of 3 mol of reagent/mol of hexameric rho protein. Although the potency of PLP-AMP is enhanced when an RNA substrate such as poly(C) is present, the stoichiometry for inhibition remains the same as in the absence of poly(C). The nucleotide substrate ATP competes very effectively for the binding site and protects against PLP-AMP inactivation. A domain of rho called N2, which comprises the distal two-thirds of the molecule (residues 152-419) and encompasses the region proposed to bind ATP, is labeled specifically in the presence of poly(C). Amino acid sequence analysis of the single [3H]PLP-AMP labeled proteolytic fragment showed Lys181 to be the site of modification, suggesting that this residue normally interacts with the gamma-phosphoryl of bound ATP. These results agree with our proposed tertiary structure for the ATP-binding domain of rho that places this lysine residue in a flexible loop above a hydrophobic nucleotide-binding pocket comprised of several parallel beta-strands, similar to adenylate kinase, F1-ATPase, and related ATP-binding proteins. Parallel studies of rho structure and function by site-directed mutagenesis and chemical modification support this interpretation.     

Kinetic properties of GABA rho1 homomeric receptors expressed in HEK293 cells

The rho1 subunit of the ionotropic GABA receptors is thought to contribute to the formation of the GABA(C) receptors with pharmacological and physiological properties distinct from those of GABA(A) receptors. Previous characterization of this subunit expressed in the Xenopus oocytes revealed an ion channel with slow activation and deactivation and no desensitization, quite different from the properties of GABA(C) receptors observed in native cells. We expressed the human rho1 subunit in human embryonic kidney (HEK) 293 cells and quantitatively characterized the kinetic properties of these receptors using a rapid drug application device. The rho1 subunit expressed in HEK293 cells exhibited pharmacological and kinetic properties qualitatively identical to those described when rho1 was expressed in the oocytes. An apparent desensitizing current observed during a constant GABA application was determined to be secondary to an E(Cl) shift. Detailed kinetic analyses and parameter estimation for a five-state kinetic model revealed that the channel is best described by a set of rate constants with a notably faster GABA unbinding K(off) rate compared to the parameters proposed for the same subunit expressed in the oocytes. The same subunit expressed in hippocampal neurons showed activation and deactivation kinetics identical to the current characterized in HEK293 cells. The kinetic properties of rho1 subunit expressed in a nonoocyte model system may be better described quantitatively by the rate constants presented here.     

Recognition of the lipoyl domain is the ultimate determinant of substrate channelling in the pyruvate dehydrogenase multienzyme complex

Reductive acetylation of the lipoyl domain (E2plip) of the dihydrolipoyl acetyltransferase component of the pyruvate dehydrogenase multienzyme complex of Escherichia coli is catalysed specifically by its partner pyruvate decarboxylase (E1p), and no productive interaction occurs with the analogous 2-oxoglutarate decarboxylase (E1o) of the 2-oxoglutarate dehydrogenase complex. Residues in the lipoyl-lysine beta-turn region of the unlipoylated E2plip domain (E2plip(apo)) undergo significant changes in both chemical shift and transverse relaxation time (T(2)) in the presence of E1p but not E1o. Residue Gly11, in a prominent surface loop between beta-strands 1 and 2 in the E2plip domain, was also observed to undergo a significant change in chemical shift. Addition of pyruvate to the mixture of E2plip(apo) and E1p caused larger changes in chemical shift and the appearance of multiple cross-peaks for certain residues, suggesting that the domain was experiencing more than one type of interaction. Residues in both beta-strands 4 and 5, together with those in the prominent surface loop and the following beta-strand 2, appeared to be interacting with E1p, as did a small patch of residues centred around Glu31. The values of T(2) across the polypeptide chain backbone were also lower than in the presence of E1p alone, suggesting that E2plip(apo) binds more tightly after the addition of pyruvate. The lipoylated domain (E2plip(holo)) also exhibited significant changes in chemical shift and decreases in the overall T(2) relaxation times in the presence of E1p, the residues principally affected being restricted to the half of the domain that contains the lipoyl-lysine (Lys41) residue. In addition, small chemical shift changes and a general drop in T(2) times in the presence of E1o were observed, indicating that E2plip(holo) can interact, weakly but non-productively, with E1o. It is evident that recognition of the protein domain is the ultimate determinant of whether reductive acetylation of the lipoyl group occurs, and that this is ensured by a mosaic of interactions with the Elp.     

Nitrogen-15-labeled yeast tRNAPhe: double and two-dimensional heteronuclear NMR of guanosine and uracil ring NH groups

5N1-Labeled hypoxanthine and 1,3-15N-labeled uracil were synthesized chemically and used to prepare labeled yeast tRNAPhe biosynthetically. Maps (500 MHz) of 15N chemical shift vs. proton chemical shift were obtained, for each ring NH group, by means of INDOR (difference heterodecoupling) and also by means of a proton-observe two-dimensional method involving coherences of forbidden resonances of the NH system. Resonances of GC11, T54-m1A58, GU4, and A psi 31 were confirmed, assigned, or reassigned. psi 39 was found to be in anti conformation, not syn as previously stated. Almost all the uracil NH group resonances could be separated, but most of the GC resonances are too close even in two dimensions to be separately resolved with the observed 20-Hz 15N line width.     

Interactions of SARS-CoV-2 envelope protein with amilorides correlate with antiviral activity

SARS-CoV-2 is the novel coronavirus that is the causative agent of COVID-19, a sometimes-lethal respiratory infection responsible for a world-wide pandemic. The envelope (E) protein, one of four structural proteins encoded in the viral genome, is a 75-residue integral membrane protein whose transmembrane domain exhibits ion channel activity and whose cytoplasmic domain participates in protein-protein interactions. These activities contribute to several aspects of the viral replication-cycle, including virion assembly, budding, release, and pathogenesis. Here, we describe the structure and dynamics of full-length SARS-CoV-2 E protein in hexadecylphosphocholine micelles by NMR spectroscopy. We also characterized its interactions with four putative ion channel inhibitors. The chemical shift index and dipolar wave plots establish that E protein consists of a long transmembrane helix (residues 8–43) and a short cytoplasmic helix (residues 53–60) connected by a complex linker that exhibits some internal mobility. The conformations of the N-terminal transmembrane domain and the C-terminal cytoplasmic domain are unaffected by truncation from the intact protein. The chemical shift perturbations of E protein spectra induced by the addition of the inhibitors demonstrate that the N-terminal region (residues 6–18) is the principal binding site. The binding affinity of the inhibitors to E protein in micelles correlates with their antiviral potency in Vero E6 cells: HMA ≈ EIPA > DMA >> Amiloride, suggesting that bulky hydrophobic groups in the 5’ position of the amiloride pyrazine ring play essential roles in binding to E protein and in antiviral activity. An N15A mutation increased the production of virus-like particles, induced significant chemical shift changes from residues in the inhibitor binding site, and abolished HMA binding, suggesting that Asn15 plays a key role in maintaining the protein conformation near the binding site. These studies provide the foundation for complete structure determination of E protein and for structure-based drug discovery targeting this protein.     

NMR identification of left-handed polyproline type II helices

NMR characteristics of a model left-handed 3(1)-helical peptide are reported in this study. With temperature and sequence corrections on the predicted random coil (15)N chemical shifts, a significant (15)N chemical shift deviation is observed for the model 3(1) peptide. The (15)N chemical shift differences also correlate well with the molar ellipticities (at 220 nm) of the CD spectra at different temperatures, indicating that the (15)N chemical shift is a sensitive probe for 3(1)-helices. The average (3)J(HNalpha) and (1)J(CalphaHalpha) values of the model peptide are determined to be 6.5 and 142.6 Hz, respectively, which are consistent with the values calculated from the geometry of 3(1)-helices. With careful measurements of amide (15)N chemical shifts and incorporating temperature and sequence effect corrections, the (15)N chemical shifts can be used together with (3)J(HNalpha) and (1)J(CalphaHalpha) to differentiate 3(1)-helices from random coils with high confidence. Based on the observed NMR characteristics, a strategy is developed for probing left-handed 3(1)-helical structures from other secondary structures.     

Carbon-13 NMR in conformational analysis of nucleic acid fragments. Heteronuclear chemical shift correlation spectroscopy of RNA constituents.

The assignment of the non-quaternary 13C resonances by means of two-dimensional heteronuclear chemical shift correlation spectroscopy is presented for several oligoribonucleotides: The dimers m6(2)AU, m6(2)Am6(2)A and mpUm6(2)A and the trimers m6(2)AUm6(2)A and m4(2)Cm4(2)Cm6(2)A. The temperature and concentration dependency of the 13C chemical shifts are studied with emphasis on the behaviour of the dimer m6(2)AU. The present study shows that in the 5-50 mM range the concentration-dependent chemical shift changes of the ribose carbons are negligible compared to chemical shift changes due to intramolecular events. All compounds studied show a surprising correlation between the chemical shifts of the carbon atoms of the ribose ring and the sugar conformational equilibrium as expressed by the percentage N or S conformer. Thus the chemical shift data can be used to obtain the thermodynamical parameters of the two-state N/S equilibrium. Parameters deduced for m6(2)AU are Tm = 306 K and delta S = -25 cal mol-1 K-1, which values are in satisfactory agreement with results obtained earlier from 1H NMR and from Circular Dichroism.