A complete Hermitian operator basis set for any spin quantum number.

A new Hermitian operator basis set for spins of any quantum number is presented for use in simulations of NMR experiments. The advantage with a Hermitian operator basis is that the Liouville-von Neumann equation, including relaxation with dynamic frequency shifts, is real. Real algebra makes numerical calculations faster and simplifies physical interpretation of the equation system as compared to complex algebra. The unity operator is included in the Hermitian operator basis, which makes it easy to rewrite the inhomogeneous Liouville-von Neumann equation into a homogeneous form. The unity operator also simplifies physical interpretation of the equation system for coupled spin systems.

Both proline-rich sequences in the TH region of Bruton's tyrosine kinase stabilize intermolecular interactions with the SH3 domain.

The Tec homology (TH) region located N-terminal to the Src homology 3 (SH3) domain of Bruton's tyrosine kinase (Btk) contains two proline-rich SH3-binding sequences (PRRs). We have previously demonstrated that the TH region acts to stabilize intermolecular interactions in N-terminally extended SH3 (PRR-SH3) fragments. Here, we analyze six PRR-SH3 fragments with different proline-to-alanine substitutions in the two PRRs. Gel permeation chromatography and nuclear magnetic resonance spectroscopy show that both PRRs can stabilize self-association. This observation provides an explanation to why the TH region of Btk makes intermolecular interactions, whereas the corresponding interaction in the related Itk kinase with only one PRR, is intramolecular.

The N-terminal regions of estrogen receptor alpha and beta are unstructured in vitro and show different TBP binding properties.

The N-terminal regions of the estrogen receptor alpha (ER alpha-N) and beta (ER beta-N) were expressed and purified to homogeneity. Using NMR and circular dichroism spectroscopy, we conclude that both ER alpha-N and ER beta-N are unstructured in solution. The TATA box-binding protein (TBP) has been shown previously to interact with ER alpha-N in vitro and to potentiate ER-activated transcription. We used surface plasmon resonance and circular dichroism spectroscopy to confirm and further characterize the ER-N-TBP interaction. Our results show that the intrinsically unstructured ER alpha-N interacts with TBP, and suggest that structural changes are induced in ER alpha-N upon TBP interaction. Conformational changes upon target factor interaction have not previously been demonstrated for any N-terminal region of nuclear receptors. In addition, no binding of ER beta-N to TBP was detected. This difference in TBP binding could imply differential recruitment of target proteins by ER alpha-N and ER beta-N. The affinity of the ER alpha-N-TBP interaction was determined to be in the micromolar range (K(D) = 10(-6) to 10(-5) m). Our results support models of TBP as a target protein for the N-terminal activation domain of ER alpha. Further, our results suggest that target proteins can induce and/or stabilize ordered structure in N-terminal regions of nuclear receptors upon interaction.

Structure and backbone dynamics of Apo-CBFbeta in solution.

Runx proteins constitute a family of mammalian transcription factors that interact with DNA through their evolutionarily conserved Runt domain. CBFbeta, alternatively denoted PEBP2beta, is the non-DNA-binding heterodimer partner and acts to enhance the DNA binding affinity of Runx proteins. Runx proteins and CBFbeta are associated with a variety of biological functions and human diseases; they are, for example, together the most frequent targets for chromosomal rearrangements in acute human leukemias. We have determined the solution structure and characterized the backbone dynamics of C-terminally truncated fragments containing residues 1-141 of CBFbeta. The present apo-CBFbeta structure is very similar to that seen in a Runt-CBFbeta complex. An evaluation of backbone (15)N NMR relaxation parameters shows that CBFbeta is a rigid molecule with high order parameters throughout the backbone; the only regions displaying significant dynamics are a long loop and the C-terminal alpha-helix. A few residues display relaxation behavior indicating conformational exchange on microsecond to millisecond time scales, but only one of these is located at the Runt binding surface. Our structure and dynamics analysis of CBFbeta therefore suggests that the protein binds to Runt without large conformational changes or induced folding ("lock-and-key" interaction). The apo-CBFbeta structure presented here exhibits several significant differences with two other published NMR ensembles of very similar protein fragments. The differences are located in four regions outside of the central beta-barrel, whereas the beta-barrel itself is almost identical in the three NMR structures. The comparison illustrates that independently determined NMR structures may display rather large differences in backbone conformation in regions that appear to be well-defined in each of the calculated NMR ensembles.

Intermolecular interactions between the SH3 domain and the proline-rich TH region of Bruton's tyrosine kinase.

The SH3 domain of Bruton's tyrosine kinase (Btk) is preceded by the Tec homology (TH) region containing proline-rich sequences. We have studied a protein fragment containing both the Btk SH3 domain and the proline-rich sequences of the TH region (PRR-SH3). Intermolecular NMR cross-relaxation measurements, gel permeation chromatography profiles, titrations with proline-rich peptides, and (15)N NMR relaxation measurements are all consistent with a monomer-dimer equilibrium with a dissociation constant on the order of 60 microM. The intermolecular interactions do, at least in part, involve proline-rich sequences in the TH region. This behavior of Btk PRR-SH3 may have implications for the functional action of Btk.

Chloride binding by the AML1/Runx1 transcription factor studied by NMR.

It is known that the DNA binding Runt domain of the AML1/Runx1 transcription factor coordinates Cl(-) ions. In this paper we have determined Cl(-) binding affinities of AML1 by (35)Cl nuclear magnetic resonance (NMR) linewidth analysis. The Runt domain binds Cl(-) with a dissociation constant (K(d,Cl)) of 34 mM. If CBFbeta is added to form a 1:1 complex, the K(d,Cl) value increases to 56 mM. Homology modeling suggests that a high occupancy Cl(-) binding site overlaps with the DNA binding surface. NMR data show that DNA displaces this Cl(-) ion. Possible biological roles of Cl(-) binding are discussed based on these findings.

Simulations of NMR pulse sequences during equilibrium and non-equilibrium chemical exchange.

The McConnell equations combine the differential equations for a simple two-state chemical exchange process with the Bloch differential equations for a classical description of the behavior of nuclear spins in a magnetic field. This equation system provides a useful starting point for the analysis of slow, intermediate and fast chemical exchange studied using a variety of NMR experiments. The McConnell equations are in the mathematical form of an inhomogeneous system of first-order differential equations. Here we rewrite the McConnell equations in a homogeneous form in order to facilitate fast and simple numerical calculation of the solution to the equation system. The McConnell equations can only treat equilibrium chemical exchange. We therefore also present a homogeneous equation system that can handle both equilibrium and non-equilibrium chemical processes correctly, as long as the kinetics is of first-order. Finally, the same method of rewriting the inhomogeneous form of the McConnell equations into a homogeneous form is applied to a quantum mechanical treatment of a spin system in chemical exchange. In order to illustrate the homogeneous McConnell equations, we have simulated pulse sequences useful for measuring exchange rates in slow, intermediate and fast chemical exchange processes. A stopped-flow NMR experiment was simulated using the equations for non-equilibrium chemical exchange. The quantum mechanical treatment was tested by the simulation of a sensitivity enhanced 15N-HSQC with pulsed field gradients during slow chemical exchange and by the simulation of the transfer efficiency of a two-dimensional heteronuclear cross-polarization based experiment as a function of both chemical shift difference and exchange rate constants.

Another piece of the ribosome: solution structure of S16 and its location in the 30S subunit.

BACKGROUND: X-ray crystallography has recently yielded much-improved electron-density maps of the bacterial ribosome and its two subunits and many structural details of bacterial ribosome subunits are now being resolved. One approach to complement the structures and elucidate the details of rRNA and protein packing is to determine structures of individual protein components and model these into existing intermediate resolution electron density. RESULTS: We have determined the solution structure of the ribosomal protein S16 from Thermus thermophilus. S16 is a mixed alpha/beta protein with a novel folding scaffold based on a five-stranded antiparallel/parallel beta sheet. Three large loops, which are partially disordered, extend from the sheet and two alpha helices are packed against its concave surface. Calculations of surface electrostatic potentials show a large continuous area of positive electrostatic potential and smaller areas of negative potential. S16 was modeled into a 5.5 A electron-density map of the T. thermophilus 30S ribosomal subunit. CONCLUSIONS: The location and orientation of S16 in a narrow crevice formed by helix 21 and several other unassigned rRNA helices is consistent with electron density corresponding to the shape of S16, hydroxyl radical protection data, and the electrostatic surface potential of S16. Two protein neighbors to S16 are S4 and S20, which facilitate binding of S16 to the 30S subunit. Overall, this work exemplifies the benefits of combining high-resolution nuclear magnetic resonance (NMR) structures of individual components with low-resolution X-ray maps to elucidate structures of large complexes.

Sequence-specific DNA binding by the glucocorticoid receptor DNA-binding domain is linked to a salt-dependent histidine protonation.

We used isothermal titration calorimetry in the temperature range 21-25 degrees C to investigate the effect of pH on the calorimetric enthalpy (delta H(cal)) for sequence specific DNA-binding of the glucocorticoid receptor DNA-binding domain (GR DBD). Titrations were carried out in solutions containing 100 mM NaCl, 1 mM dithiothreitol, 5% glycerol by volume, and 20 mM Tris, Hepes, Mops, or sodium phosphate buffers at pH 7.5. A strong dependence of delta H(cal) on the buffer ionization enthalpy is observed, demonstrating that the DNA binding of the GR DBD is linked to proton uptake at these conditions. The apparent increase in the pK(a) for an amino acid side chain upon DNA binding is supported by the results of complementary titrations, where delta H(cal) shows a characteristic dependence on the solution pH. delta H(cal) is also a function of the NaCl concentration, with opposite dependencies in Tris and Hepes buffers, respectively, such that a similar delta H(cal) value is approached at 300 mM NaCl. This behavior shows that the DNA-binding induced protonation is inhibited by increased concentrations of NaCl. A comparison with structural data suggests that the protonation involves a histidine (His451) in the GR DBD, because in the complex this residue is located close to a DNA phosphate at an orientation that is consistent with a charged-charged hydrogen bond in the protonated state. NMR spectra show that His451 is not protonated in the unbound protein at pH 7.5. The pH dependence in delta H(cal) can be quantitatively described by a shift of the pK(a) of His451 from approximately 6 in the unbound state to close to 8 when bound to DNA at low salt concentration conditions. A simple model involving a binding competition between a proton and a Na(+) counterion to the GR DBD-DNA complex reproduces the qualitative features of the salt dependence.

The solution structure of ribosomal protein L36 from Thermus thermophilus reveals a zinc-ribbon-like fold.

We have determined the solution NMR structure of the ribosomal protein L36 from Thermus thermophilus. L36 is the smallest protein in the large subunit of the prokaryotic ribosome. The sequence contains three completely conserved cysteine residues and one conserved histidine residue in a C-X(2)-C-X(12)-C-X(4)-H motif. Extended X-ray absorption fine structure spectroscopy was used to confirm that a purified L36 sample contains an equimolar amount of zinc. The structure of L36 was determined using simulated annealing based on NOE distance restraints, dihedral angle restraints and hydrogen bond distance restraints derived from NMR spectra of (15)N-labeled and non-labeled L36 samples at pH 7 and 12 degrees C, and by imposing tetrahedral zinc ion coordination geometry. The L36 fold is characterized by a triple-stranded antiparallel beta-sheet with the zinc-binding site at one end. The structure of the zinc site is well-determined and shows that the three cysteine sulphur atoms are supported by hydrogen bonds to backbone amide protons. The conserved histidine residue is located in a short 3(10)-helix and coordinates zinc by the N(delta1) atom. The electrostatic surface potential and location of conserved Arg, Lys and His side-chains suggest a large continuous L36-rRNA interaction interface. The folding topology as well as position and conformation of many conserved side-chains in L36 are very similar to those of zinc-ribbon domains found in the archaeal transcription factor TFIIB N terminus and the eukaryal transcription elongation factor hTFIIS C terminus. Given the relative antiquity of the ribosome it is possible that L36 reflects the parent of transcription-related zinc ribbons.

Ansig for Windows: an interactive computer program for semiautomatic assignment of protein NMR spectra.

Assignment of NMR spectra is a prerequisite for structure determination of proteins using NMR. The time spent on the assignment is comparatively long compared to that spent on other parts in the protein structure determination process, but it can be shortened by using either interactive or fully automated computer programs. To benefit from the advantages of both types of program we have developed a version of the interactive assignment program ANSIG to include automatized, yet user-supervised, routines. The new program includes tools for (i) semiautomatic sequential assignment, (ii) plotting of distances from PDB structure files directly in NMR spectra and (iii) statistical analysis of distance restraint violations with the possibility to directly zoom to violated NOEs in NOESY spectra.

Solution structure of the ribosomal protein S19 from Thermus thermophilus.

Ribosomal protein S19 is a 10.6 kDa protein in the small subunit of the prokaryotic ribosome. We have determined a high-resolution solution structure of S19 from Thermus thermophilus. Structures were calculated using 1160 distance and dihedral angle restraints derived from (1)H, (15)N and (13)C NMR spectra. The structures show that S19 is a mixed alpha/beta protein with long disordered tails. The folding topology is not homologous to that of any other known protein structure. Potential rRNA and protein binding sites have been identified on the S19 surface.

Conformational dynamics and molecular recognition: backbone dynamics of the estrogen receptor DNA-binding domain.

We examined the internal mobility of the estrogen receptor DNA-binding domain (ER DBD) using NMR15N relaxation measurements and compared it to that of the glucocorticoid receptor DNA-binding domain (GR DBD). The studied protein fragments consist of residues Arg183-His267 of the human ER and residues Lys438-Gln520 of the rat GR. The15N longitudinal (R1) and transverse (R2) relaxation rates and steady state {1H}-15N nuclear Overhauser enhancements (NOEs) were measured at 30 degrees C at1H NMR frequencies of 500 and 600 MHz. The NOE versus sequence profile and calculated order parameters for ER DBD backbone motions indicate enhanced internal dynamics on pico- to nanosecond time-scales in two regions of the core DBD. These are the extended strand which links the DNA recognition helix to the second zinc domain and the larger loop region of the second zinc domain. The mobility of the corresponding regions of the GR DBD, in particular that of the second zinc domain, is more limited. In addition, we find large differences between the ER and GR DBDs in the extent of conformational exchange mobility on micro- to millisecond time-scales. Based on measurements of R2as a function of the15N refocusing (CPMG) delay and quantitative (Lipari-Szabo-type) analysis, we conclude that conformational exchange occurs in the loop of the first zinc domain and throughout most of the second zinc domain of the ER DBD. The conformational exchange dynamics in GR DBD is less extensive and localized to two sites in the second zinc domain. The different dynamical features seen in the two proteins is consistent with previous studies of the free state structures in which the second zinc domain in the ER DBD was concluded to be disordered whereas the corresponding region of the GR DBD adopts a stable fold. Moreover, the regions of the ER DBD that undergo conformational dynamics on the micro- to millisecond time-scales in the free state are involved in intermolecular protein-DNA and protein-protein interactions in the dimeric bound state. Based on the present data and the previously published dynamical and DNA binding properties of a GR DBD triple mutant which recognize an ER binding site on DNA, we argue that the free state dynamical properties of the nuclear receptor DBDs is an important element in molecular recognition upon DNA binding.

Solution properties of the free and DNA-bound Runt domain of AML1.

The Runt domain is responsible for specific DNA and protein-protein interactions in a family of transcription factors which includes human AML1. Structural data on the Runt domain has not yet become available, possibly due to solubility and stability problems with expressed protein fragments. Here we describe the optimization and characterization of a 140-residue fragment, containing the Runt domain of AML1, which is suitable for structural studies. The fragment of AML1 including amino acids 46-185 [AML1 Dm(46-185)] contains a double cysteine-->serine mutation which does not affect Runt domain structure or DNA-binding affinity. Purified AML1 Dm(46-185) is soluble and optimally stable in a buffer containing 200 mm MgSO4 and 20 mm sodium phosphate at pH 6.0. Nuclear magnetic resonance and circular dichroism spectroscopy indicate that the Runt domain contains beta-sheet, but little or no alpha-helical secondary structure elements. The 45 N-terminal residues of AML1 are unstructured and removal of the N-terminal enhances sequence-specific DNA binding. The NMR spectrum of AML1 Dm(46-185) displays a favorable chemical shift dispersion and resolved NOE connectivities are readily identified, suggesting that a structure determination of this Runt domain fragment is feasible. A titration of 15N-labelled AML1 Dm(46-185) with a 14-bp cognate DNA duplex results in changes in the 15N NMR heteronuclear single quantum coherence spectrum which indicate the formation of a specific complex and structural changes in the Runt domain upon DNA binding.

The complete homogeneous master equation for a heteronuclear two-spin system in the basis of cartesian product operators

The complete homogeneous form of the quantum mechanical master equation for a heteronuclear two-spin system is presented in the basis of Cartesian product operators. The homogeneous master equation is useful since it allows fast, single-step computation of the density operator during pulse sequences, without neglecting relaxation effects. The homogeneous master equation is also a prerequisite for an expansion of the average Hamiltonian theory to include relaxation, thus forming average Liouvillian theory. The coherences of the two-spin system are assumed to be relaxed both by mutual dipole-dipole interaction and by chemical shift anisotropy interaction with the static magnetic field. The cross-correlation between dipole-dipole and chemical shift anisotropy relaxation mechanisms is also considered. To illustrate the applicability of the developed formalism we simulate the overall transfer efficiency of three different inverse detection 1H-15N correlation experiments with parameters corresponding to a large protein. Copyright 1998 Academic Press.

Architecture of nonspecific protein-DNA interactions in the Sso7d-DNA complex.

Many biochemical processes, including DNA packing, maintenance and control, rely on non-sequence specific protein-DNA interactions. Nonspecific DNA-binding proteins have evolved to tolerate a wide range of DNA sequences, yet bind with a respectable affinity. The nonspecific binding requirement is in contrast to that imposed on, for example, transcription factors and implies a different structural basis for the biomolecular recognition process. To address this issue, and the mechanism for archaeal DNA packing, we determined the structure of the Sso7d protein from Sulfolobus solfataricus in complex with DNA. Sso7d binds DNA by placing a triple-stranded beta-sheet across the DNA minor groove. The protein is anchored in this position by the insertion of hydrogen bond-donating side chains into the groove and additionally stabilized by electrostatic and non-polar interactions with the DNA backbone. This structure explains how strong binding can be achieved independent of DNA sequence. Sso7d binding also distorts the DNA conformation and introduces significant unwinding of the helix. This effect suggests a mechanism for DNA packing in Sulfolobus based on negative DNA supercoiling.

Thermodynamic characterization of non-sequence-specific DNA-binding by the Sso7d protein from Sulfolobus solfataricus.

We used isothermal titration calorimetry and fluorescence spectroscopy to investigate the thermodynamics of non-sequence-specific DNA-binding by the Sso7d protein from the archaeon Sulfolobus solfataricus. We report the Sso7d-poly(dGdC) binding thermodynamics as a function of buffer composition (Tris-HCl or phosphate), temperature (15 to 45 degrees C), pH (7.1 to 8.0), osmotic stress and solvent (H2O/2H2O), and compare it to poly (dAdT) binding; and we have previously also reported the salt concentration dependence. Binding isotherms can be represented by the McGhee-von Hippel model for non-cooperative binding, with a binding site size of four to five DNA base-pairs and binding free energies in the range DeltaG degrees approximately -7 to DeltaG degrees approximately -10 kcal mol-1, depending on experimental conditions. The non-specific nature of the binding is reflected in similar thermodynamics for binding to poly(dAdT) and poly(dGdC). The native lysine methylation of Sso7d has only minor effects on the binding thermodynamics. Sso7d binding to poly(dGdC) is endothermic at 25 degrees C with a binding enthalpy DeltaH degrees approximately 10 kcal mol-1 in both phosphate and Tris-HCl buffers at pH 7.6, indicating that DeltaH degrees does not include large contributions from coupled buffer ionization equilibria at this pH. The binding enthalpy is temperature dependent with a measured heat capacity change DeltaCp degrees=-0.25(+/-0.01) kcal mol-1 K-1 and extrapolations of thermodynamic data indicate that the complex is heat stable with exothermic binding close to the growth temperature (75 to 80 degreesC) of S. solfataricus. Addition of neutral solutes (osmotic stress) has minor effects on DeltaG degrees and the exchange of H2O for 2H2O has only a small effect on DeltaH degrees, consistent with the inference that complex formation is not accompanied by net changes in surface hydration. Thus, other mechanisms for the heat capacity change must be found. The observed thermodynamics is discussed in relation to the nature of non-sequence-specific DNA-binding by proteins.

Solution structure of the SH3 domain from Bruton's tyrosine kinase.

X-linked agammaglobulinemia (XLA) is a heritable immunodeficiency caused by mutations in the gene coding for Bruton's tyrosine kinase (Btk). Btk belongs to the Tec family of tyrosine kinases. Each member of the family contains five regions and mutations causing XLA have been isolated in all five regions. We have determined the solution structure of the Src homology 3 (SH3) domain of Btk using two- and three-dimensional nuclear magnetic resonance (NMR) spectroscopy on natural abundance and 15N-labeled protein material. The structure determination is complemented by investigation of backbone dynamics based on 15N NMR relaxation. The Btk SH3 forms a well-defined structure and shows the typical SH3 topology of two short antiparallel beta-sheets packed almost perpendicular to each other in a sandwich-like fold. The N- and C-termini are more flexible as are peptide fragments in the RT and n-Src loops. The studied Btk SH3 fragment adopts two slowly interconverting conformations with a relative concentration ratio of 7:1. The overall fold of the minor form is similar to that of the major form, as judged on the basis of observed NOE connectivities and small chemical shift differences. A tryptophan (W251) ring flip is the favored mechanism for interconversion, although other possibilities cannot be excluded. The side chain of Y223, which becomes autophosphorylated upon activation of Btk, is exposed within the potential SH3 ligand binding site. Finally, we compare the present Btk SH3 structure with other SH3 structures.