Heteronuclear relaxation study of the PH domain of beta-spectrin: restriction of loop motions upon binding inositol trisphosphate

Non-inositol 1,4,5-trisphosphate (IP3) receptor IP3-binding proteins

Conventionally, myo-D-inositol 1, 4,5-trisphosphate (IP₃) is thought to exert its second messenger effects through the gating of IP₃R Ca²⁺ release channels, located in Ca²⁺-storage organelles like the endoplasmic reticulum. However, there is considerable indirect evidence to support the concept that IP₃ might interact with other, non-IP₃R proteins within cells. To explore this possibility further, the Protein Data Bank was searched using the term "IP3". This resulted in the retrieval of 203 protein structures, the majority of which were members of the IP₃R/ryanodine receptor superfamily of channels. Only 49 of these structures were complexed with IP₃. These were inspected for their ability to interact with the carbon-1 phosphate of IP₃, since this is the least accessible phosphate group of its precursor, phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂). This reduced the number of structures retrieved to 35, of which 9 were IP₃Rs. The remaining 26 structures represent a diverse range of proteins, including inositol-lipid metabolizing enzymes, signal transducers, PH domain containing proteins, cytoskeletal anchor proteins, the TRPV4 ion channel, a retroviral Gag protein and fibroblast growth factor 2. Such proteins may impact on IP₃ signalling and its effects on cell-biology. This represents an area open for exploration in the field of IP₃ signalling.

Structure and dynamics of the von Willebrand Factor C6 domain

Von Willebrand disease (VWD) is a bleeding disorder with different levels of severity. VWD-associated mutations are located in the von Willebrand factor (VWF) gene, coding for the large multidomain plasma protein VWF with essential roles in hemostasis and thrombosis. On the one hand, a variety of mutations in the C-domains of VWF are associated with increased bleeding upon vascular injury. On the other hand, VWF gain-of-function (GOF) mutations in the C4 domain have recently been identified, which induce an increased risk of myocardial infarction. Mechanistic insights into how these mutations affect the molecular behavior of VWF are scarce and holistic approaches are challenging due to the multidomain and multimeric character of this large protein. Here, we determine the structure and dynamics of the C6 domain and the single nucleotide polymorphism (SNP) variant G2705R in C6 by combining nuclear magnetic resonance spectroscopy, molecular dynamics simulations and aggregometry. Our findings indicate that this mutation mostly destabilizes VWF by leading to a more pronounced hinging between both subdomains of C6. Hemostatic parameters of variant G2705R are close to normal under static conditions, but the missense mutation results in a gain-of-function under flow conditions, due to decreased VWF stem stability. Together with the fact that two C4 variants also exhibit GOF characteristics, our data underline the importance of the VWF stem region in VWF's hemostatic activity and the risk of mutation-associated prothrombotic properties in VWF C-domain variants due to altered stem dynamics.

Paip2A inhibits translation by competitively binding to the RNA recognition motifs of PABPC1 and promoting its dissociation from the poly(A) tail

Eukaryotic mRNAs possess a poly(A) tail at their 3′-end, to which poly(A)-binding protein C1 (PABPC1) binds and recruits other proteins that regulate translation. Enhanced poly(A)-dependent translation, which is also PABPC1 dependent, promotes cellular and viral proliferation. PABP-interacting protein 2A (Paip2A) effectively represses poly(A)-dependent translation by causing the dissociation of PABPC1 from the poly(A) tail; however, the underlying mechanism remains unknown. This study was conducted to investigate the functional mechanisms of Paip2A action by characterizing the PABPC1–poly(A) and PABPC1–Paip2A interactions. Isothermal titration calorimetry and NMR analyses indicated that both interactions predominantly occurred at the RNA recognition motif (RRM)2–RRM3 regions of PABPC1, which have comparable affinities for poly(A) and Paip2A (dissociation constant, K_d = 1 nM). However, the K_d values of isolated RRM2 were 200 and 4 μM in their interactions with poly(A) and Paip2A, respectively; K_d values of 5 and 1 μM were observed for the interactions of isolated RRM3 with poly(A) and Paip2A, respectively. NMR analyses also revealed that Paip2A can bind to the poly(A)-binding interfaces of the RRM2 and RRM3 regions of PABPC1. Based on these results, we propose the following functional mechanism for Paip2A: Paip2A initially binds to the RRM2 region of poly(A)-bound PABPC1, and RRM2-anchored Paip2A effectively displaces the RRM3 region from poly(A), resulting in dissociation of the whole PABPC1 molecule. Together, our findings provide insight into the translation repression effect of Paip2A and may aid in the development of novel anticancer and/or antiviral drugs.

A Distinct, Sequence-Induced Conformation Is Required for Recognition of the Constitutive Decay Element RNA by Roquin

The constitutive decay element (CDE) of tumor necrosis factor α (TNF-α) mRNA (Tnf) represents the prototype of a class of RNA motifs that mediate rapid degradation of mRNAs encoding regulators of the immune response and development. CDE-type RNAs are hairpin structures featuring a tri-nucleotide loop. The protein Roquin recognizes CDE-type stem loops and recruits the Ccr4-Caf1-Not deadenylase complex to the mRNA, thereby inducing its decay. Stem recognition does not involve nucleotide bases; however, there is a strong stem sequence requirement for functional CDEs. Here, we present the solution structures of the natural Tnf CDE and of a CDE mutant with impaired Roquin binding. We find that the two CDEs adopt unique and distinct structures in both the loop and the stem, which explains the ability of Roquin to recognize stem loops in a sequence-specific manner. Our findings result in a relaxed consensus motif for prediction of new CDE stem loops.

Molecular basis of UG-rich RNA recognition by the human splicing factor TDP-43

TDP-43 encodes an alternative-splicing regulator with tandem RNA-recognition motifs (RRMs). The protein regulates cystic fibrosis transmembrane regulator (CFTR) exon 9 splicing through binding to long UG-rich RNA sequences and is found in cytoplasmic inclusions of several neurodegenerative diseases. We solved the solution structure of the TDP-43 RRMs in complex with UG-rich RNA. Ten nucleotides are bound by both RRMs, and six are recognized sequence specifically. Among these, a central G interacts with both RRMs and stabilizes a new tandem RRM arrangement. Mutations that eliminate recognition of this key nucleotide or crucial inter-RRM interactions disrupt RNA binding and TDP-43-dependent splicing regulation. In contrast, point mutations that affect base-specific recognition in either RRM have weaker effects. Our findings reveal not only how TDP-43 recognizes UG repeats but also how RNA binding-dependent inter-RRM interactions are crucial for TDP-43 function.

Rescue of the aggregation prone Itk Pleckstrin Homology domain by two mutations derived from the related kinases, Btk and Tec

IL-2 inducible T-cell kinase (Itk) is a Tec family non-receptor tyrosine kinase involved in signaling downstream of the T-cell receptor. Itk contains an amino-terminal Pleckstrin Homology (PH) domain that binds phosphatidylinositol (3,4,5)-trisphosphate, recruiting Itk to the plasma membrane upon T-cell receptor activation. In addition to phosphoinositide binding, accumulating data suggest that the Itk PH domain likely mediates additional interactions outside of the phosphoinositide ligand binding pocket. The structural basis for additional PH domain functions remains elusive because of the poor recombinant expression and in vitro solution behavior of the Itk PH domain. Here, we determine that the lone α-helix in the Itk PH domain is responsible for the poor solution properties and that mutation of just two residues in the Itk α-helix to the corresponding amino acids in Btk or Tec dramatically improves the soluble recombinant expression and solution behavior of the Itk PH domain. We present this double mutant as a valuable tool to characterize the structure and function of the Itk PH domain. It is also interesting to note that the precise sites of mutation identified in this study appear as somatic mutations associated with cancerous tissue. Collectively, the findings suggest that the two helical residues in the Itk PH domain may serve an important and unique structural role in wild-type Itk that differentiates this tyrosine kinase from its related family members.

Structural determinants of the specificity of a membrane binding domain of the scaffold protein Ste5 of budding yeast: implications in signaling by the scaffold protein in MAPK pathway

In the mitogen activated protein kinase (MAPK) cascades of budding yeast, the scaffold protein Ste5 is recruited to the plasma membrane to transmit pheromone induced signal. A region or domain of Ste5 i.e. residues P44-R67, referred here as Ste5PM24, has been known to be involved in direct interactions with the membrane. In order to gain structural insights into membrane interactions of Ste5, here, we have investigated structures and interactions of two synthetic peptide fragments of Ste5, Ste5PM24, and a hyperactive mutant, Ste5PM24LM, by NMR, ITC, and fluorescence spectroscopy, with lipid membranes. We observed that Ste5PM24 predominantly interacted only with the anionic lipid vesicles. By contrast, Ste5PM24LM exhibited binding with negatively charged as well as zwitterionic or mixed lipid vesicles. Binding of Ste5 peptides with the negatively charged lipid vesicles were primarily driven by hydrophobic interactions. NMR studies revealed that Ste5PM24 assumes dynamic or transient conformations in zwitterionic dodecylphosphocholine (DPC) micelles. By contrast, NMR structure, obtained in anionic sodium dodecyl sulphate (SDS), demonstrated amphipathic helical conformations for the central segment of Ste5PM24. The hydrophobic surface of the helix was found to be buried inside the micelles. Taken together, these results provide important insights toward the structure and specificity determinants of the scaffold protein interactions with the plasma membrane.

An NMR database for simulations of membrane dynamics

Computational methods are powerful in capturing the results of experimental studies in terms of force fields that both explain and predict biological structures. Validation of molecular simulations requires comparison with experimental data to test and confirm computational predictions. Here we report a comprehensive database of NMR results for membrane phospholipids with interpretations intended to be accessible by non-NMR specialists. Experimental ¹³C-¹H and ²H NMR segmental order parameters (S(CH) or S(CD)) and spin-lattice (Zeeman) relaxation times (T(1Z)) are summarized in convenient tabular form for various saturated, unsaturated, and biological membrane phospholipids. Segmental order parameters give direct information about bilayer structural properties, including the area per lipid and volumetric hydrocarbon thickness. In addition, relaxation rates provide complementary information about molecular dynamics. Particular attention is paid to the magnetic field dependence (frequency dispersion) of the NMR relaxation rates in terms of various simplified power laws. Model-free reduction of the T(1Z) studies in terms of a power-law formalism shows that the relaxation rates for saturated phosphatidylcholines follow a single frequency-dispersive trend within the MHz regime. We show how analytical models can guide the continued development of atomistic and coarse-grained force fields. Our interpretation suggests that lipid diffusion and collective order fluctuations are implicitly governed by the viscoelastic nature of the liquid-crystalline ensemble. Collective bilayer excitations are emergent over mesoscopic length scales that fall between the molecular and bilayer dimensions, and are important for lipid organization and lipid-protein interactions. Future conceptual advances and theoretical reductions will foster understanding of biomembrane structural dynamics through a synergy of NMR measurements and molecular simulations.

Periplasmic loop P2 of the MalF subunit of the maltose ATP binding cassette transporter is sufficient to bind the maltose binding protein MalE

The Escherichia coli maltose transporter belongs to the ATP binding cassette (ABC) transporter superfamily. Recently, the crystal structure of the full transporter MalFGK2 in complex with the maltose binding protein (MBP) was determined [Oldham, M. L., et al. (2007) Crystal structure of a catalytic intermediate of the maltose transporter. Nature 450, 515-522]. Using liquid-state NMR, we find that the periplasmic loop P2 of MalF (MalF-P2) folds independently in solution and adopts a well-defined tertiary structure which is similar to the one found in the crystal. MalF-P2 interacts with the maltose binding protein, independent of the transmembrane region of MalF and MalG with an affinity of 10-20 microM, in the presence and absence of substrate. Analysis of residual dipolar coupling (RDC) experiments shows that the conformation of the two individual domains of MalF-P2 is preserved in the absence of MalE and resembles the conformation in the X-ray structure. Upon titration of MalE to MalF-P2, the two domains of MalF-P2 change their relative orientation to accommodate the ligand. In particular, a conformational change of domain 2 of MalF-P2 is induced, which is distinct from the conformation found in the X-ray structure.

NMR structural studies of the Ste11 SAM domain in the dodecyl phosphocholine micelle

The sterile alpha-motif (SAM), a relatively small ( approximately 70 amino acids) interaction domain, is found in a variety of proteins involved in cell signaling, transcription regulation, and scaffolding. The Ste11 protein kinase from the mitogen activated protein kinase (MAPK) signaling cascades of the budding yeast is regulated by a SAM domain located at the N-terminus of full-length protein. In solution, the Ste11 SAM domain exists as a well-folded dimeric structure that is involved in interaction with the cognate SAM domain from an adaptor protein Ste50. In this work, we show that the Ste11 SAM domain has an intrinsic affinity towards the lipid membranes. The solution conformation of the Ste11 SAM determined in perdeuterated DPC micelle, using NMR spectroscopy, is defined by five helices of different lengths connected by a number of loops. In the micelle bound state, the non-polar and aromatic residues of the Ste11 SAM lack a native-like packing and are presumably engaged in interactions with the micelle. Using two different paramagnetic doxyl-lipids; we have mapped out localization of Ste11 SAM residues at the micelle surface. Most of the residues appear to localize at the interfacial region of the micelle. However, a number of non-polar residues from the central region of the domain are found to be located inside the core of the micelle including residues from the helix 4 and a loop between helix 2 and helix 3. Isothermal titration calorimetry studies demonstrate that a facile insertion of the Ste11 SAM into the DPC micelle is primarily driven by a large change in enthalpy, -50 kcal/mol with an apparent equilibrium association constant (Ka) of 7.86 x 10(6) M(-1). Interestingly, an interfacial mutant L60R of the Ste11 SAM lacking the dimeric structure does not show detectable interactions with the lipid micelle. The micelle-bound structure of the Ste11 SAM domain described in this work may have potential implications in the regulation of MAPK signaling whereby positioning of the Ste11 protein in close proximity to the membrane may facilitate efficient phosphorylation of the Ste11 kinase by the membrane attached upstream Ste20/pak kinase.

Interactions of the pleckstrin homology domain with phosphatidylinositol phosphate and membranes: characterization via molecular dynamics simulations

The mechanism of interaction of pleckstrin homology (PH) domains with phosphatidylinositol 4,5-bisphosphate (PIP 2)-containing lipid bilayers remains uncertain. While crystallographic studies have emphasized PH-inositol 1,4,5-trisphosphate (IP 3) interactions, biophysical studies indicate a degree of less specific protein-bilayer interactions. We have used molecular dynamics simulations to characterize the interactions of the PH domain from phospholipase C-delta1 with IP 3 and with PIP 2, the latter in lipid bilayers and in detergent micelles. Simulations of the PH domain in water reveal a reduction in protein flexibility when IP 3 is bound. Simulations of the PH domain bound to PIP 2 in lipid bilayers indicate a tightening of ligand-protein interactions relative to the PH-IP 3 complex, alongside formation of H-bonds between PH side chains and lipid (PC) headgroups, and a degree of penetration of hydrophobic side chains into the core of the bilayer. Comparison with simulations of the PH-bound domain to a PC bilayer in the absence of PIP 2 suggests that the presence of PIP 2 increases the extent of PH-membrane interactions. Thus, comparative molecular dynamics simulations reveal how a PI-binding domain undergoes changes in conformational dynamics on binding to a PIP 2-containing membrane and how interactions additional to those with the PI headgroup are formed.

Solution structures, dynamics, and lipid-binding of the sterile alpha-motif domain of the deleted in liver cancer 2

The deleted in liver cancer 2 (DLC2) is a tumor suppressor gene, frequently found to be underexpressed in hepatocellular carcinoma. DLC2 is a multidomain protein containing a sterile alpha-motif (SAM) domain, a GTPase-activating protein (GAP) domain, and a lipid-binding StAR-related lipid-transfer (START) domain. The SAM domain of DLC2, DLC2-SAM, exhibits a low level of sequence homology (15-30%) with other SAM domains, and appears to be the prototype of a new subfamily of SAM domains found in DLC2-related proteins. In the present study, we have determined the three-dimensional solution structure of DLC2-SAM using NMR methods together with molecular dynamics simulated annealing. In addition, we performed a backbone dynamics study. The DLC2-SAM packed as a unique four alpha-helical bundle stabilized by interhelix hydrophobic interactions. The arrangement of the four helices is distinct from all other known SAM domains. In contrast to some members of the SAM domain family which form either dimers or oligomers, both biochemical analyses and rotational correlation time (tau(c)) measured by backbone 15N relaxation experiments indicated that DLC2-SAM exists as a monomer in solution. The interaction of DLC2-SAM domain with sodium dodecyl sulfate (SDS) micelles and 1,2-dimyristoyl-sn-glycerol-3-phosphatidylglycerol (DMPG) phospholipids was examined by CD and NMR spectroscopic techniques. The DLC2-SAM exhibits membrane binding properties accompanied by minor loss of the secondary structure of the protein. Deletion studies showed that the self-association of DLC2 in vivo does not require SAM domain, instead, a protein domain consisting of residues 120-672 mediates the self-association of DLC2.

Backbone dynamics of SDF-1alpha determined by NMR: interpretation in the presence of monomer-dimer equilibrium

SDF-1alpha is a member of the chemokine family implicated in various reactions in the immune system. The interaction of SDF-1alpha with its receptor, CXCR4, is responsible for metastasis of a variety of cancers. SDF-1alpha is also known to play a role in HIV-1 pathogenesis. The structures of SDF-1alpha determined by NMR spectroscopy have been shown to be monomeric while X-ray structures are dimeric. Biochemical data and in vivo studies suggest that dimerization is likely to be important for the function of chemokines. We report here the dynamics of SDF-1alpha determined through measurement of main chain (15)N NMR relaxation data. The data were obtained at several concentrations of SDF-1alpha and used to determine a dimerization constant of approximately 5 mM for a monomer-dimer equilibrium. The dimerization constant was subsequently used to extrapolate values for the relaxation data corresponding to monomeric SDF-1alpha. The experimental relaxation data and the extrapolated data for monomeric SDF-1alpha were analyzed using the model free approach. The model free analysis indicated that SDF-1alpha is rigid on the nano- to picosecond timescale with flexible termini. Several residues involved in the dimer interface display slow micro- to millisecond timescale motions attributable to chemical exchange such as monomer-dimer equilibrium. NMR relaxation measurements are shown to be applicable for studying oligomerization processes such as the dimerization of SDF-1alpha.

Backbone dynamics of human parathyroid hormone (1-34): flexibility of the central region under different environmental conditions

The presence of a stable tertiary structure in the bioactive N-terminal portion of parathyroid hormone (PTH), a major hormone in the maintenance of extracellular calcium homeostasis, is still debated. In this work, 15N relaxation parameters of the 33 backbone amides of human PTH(1-34) were determined in phosphate-buffered saline solution (PBS) and in the presence of dodecylphosphocholine (DPC) micelles. The relaxation parameters were analyzed using both the model-free formalism (G. Lipari and A. Szabo, Journal of the American Chemical Society, 1982, Vol. 104, pp. 4546-4549) and the reduced spectral density functions approach (J.-F. Lefevre, K. T. Dayie, J. W. Peng, and G. Wagner, Biochemistry, 1996, Vol. 35, pp. 2674-2686). In PBS, the region around Gly12 possesses a high degree of flexibility and the C-terminal helix is less flexible than the N-terminal one. In the presence of DPC micelles, the mobility of the entire molecule is reduced, but the stability of the N-terminal helix increases relative to the C-terminal one. A point of relatively higher mobility at residue Gly12 is still present and a new site of local mobility at residues 16-17 is generated. These results justify the lack of experimental nuclear Overhauser effect (NOE) restraints with lack of tertiary structure and support the hypothesis that, in the absence of the receptor, the relative spatial orientation of the two N- and C-terminal helices is undefined. The flexibility in the midregion of PTH(1-34), maintained in the presence of the membrane-mimetic environment, may enable the correct relative disposition of the two helices, favoring a productive interaction with the receptor.

Structure and dynamics of the human pleckstrin DEP domain: distinct molecular features of a novel DEP domain subfamily

Pleckstrin1 is a major substrate for protein kinase C in platelets and leukocytes, and comprises a central DEP (disheveled, Egl-10, pleckstrin) domain, which is flanked by two PH (pleckstrin homology) domains. DEP domains display a unique alpha/beta fold and have been implicated in membrane binding utilizing different mechanisms. Using multiple sequence alignments and phylogenetic tree reconstructions, we find that 6 subfamilies of the DEP domain exist, of which pleckstrin represents a novel and distinct subfamily. To clarify structural determinants of the DEP fold and to gain further insight into the role of the DEP domain, we determined the three-dimensional structure of the pleckstrin DEP domain using heteronuclear NMR spectroscopy. Pleckstrin DEP shares main structural features with the DEP domains of disheveled and Epac, which belong to different DEP subfamilies. However, the pleckstrin DEP fold is distinct from these structures and contains an additional, short helix alpha4 inserted in the beta4-beta5 loop that exhibits increased backbone mobility as judged by NMR relaxation measurements. Based on sequence conservation, the helix alpha4 may also be present in the DEP domains of regulator of G-protein signaling (RGS) proteins, which are members of the same DEP subfamily. In pleckstrin, the DEP domain is surrounded by two PH domains. Structural analysis and charge complementarity suggest that the DEP domain may interact with the N-terminal PH domain in pleckstrin. Phosphorylation of the PH-DEP linker, which is required for pleckstrin function, could regulate such an intramolecular interaction. This suggests a role of the pleckstrin DEP domain in intramolecular domain interactions, which is distinct from the functions of other DEP domain subfamilies found so far.

Solution structure of the dimeric SAM domain of MAPKKK Ste11 and its interactions with the adaptor protein Ste50 from the budding yeast: implications for Ste11 activation and signal transmission through the Ste50-Ste11 complex

Ste11, a homologue of mammalian MAPKKKs, together with its binding partner Ste50 works in a number of MAPK signaling pathways of Saccharomyces cerevisiae. Ste11/Ste50 binding is mediated by their sterile alpha motifs or SAM domains, of which homologues are also found in many other intracellular signaling and regulatory proteins. Here, we present the solution structure of the SAM domain or residues D37-R104 of Ste11 and its interactions with the cognate SAM domain-containing region of Ste50, residues M27-Q131. NMR pulse-field-gradient (PFG) and rotational correlation time measurements (tauc) establish that the Ste11 SAM domain exists predominantly as a symmetric dimer in solution. The solution structure of the dimeric Ste11 SAM domain consists of five well-defined helices per monomer packed into a compact globular structure. The dimeric structure of the SAM domain is maintained by a novel dimer interface involving interactions between a number of hydrophobic residues situated on helix 4 and at the beginning of the C-terminal long helix (helix 5). The dimer structure may also be stabilized by potential salt bridge interactions across the interface. NMR H/2H exchange experiments showed that binding of the Ste50 SAM to the Ste11 SAM very likely involves the positively charged extreme C-terminal region as well as exposed hydrophobic patches of the dimeric Ste11 SAM domain. The dimeric structure of the Ste11 SAM and its interactions with the Ste50 SAM may have important roles in the regulation and activation of the Ste11 kinase and signal transmission and amplifications through the Ste50-Ste11 complex.

Mapping of the interaction interface of DNA polymerase beta with XRCC1

Residues of DNA polymerase beta (beta-Pol) that interact with the DNA repair protein XRCC1 have been determined by NMR chemical shift mapping (CSM) and mutagenesis. 15N/(13)C/(2)H/(1)H,(13)C-methyl(Leu,Ile,Val)-labeled beta-Pol palm-thumb domain was used for assignments of the 1H, 15N, and 13C resonances used for CSM of the palm-thumb on forming the 40 kDa complex with the XRCC1 N-terminal domain (NTD). Large chemical shift changes were observed in the thumb on complexation. 15N relaxation data indicate reduction in high-frequency motion for a thumb loop and three palm turn/loops, which showed concomitant chemical shift changes on complexation. A deltaV303-V306 deletion and an L301R/V303R/V306R triple mutation abolished complex formation due to loss in hydrophobicity. In an updated model, the thumb-loop of beta-Pol contacts an edge/face region of the beta sheet of the XRCC1 NTD, while the beta-Pol palm weakly contacts the alpha2 helix.

Simulated and NMR-derived backbone dynamics of a protein with significant flexibility: a comparison of spectral densities for the betaARK1 PH domain

A 7.6 ns molecular dynamics trajectory of the betaARK1 PH domain in explicit water with appropriate ions was calculated at 300 K. Spectral densities at omega = 0, omega(N), and 0.87omega(H) and the model-free parameters were evaluated from the experimental as well as the simulated data, taking the anisotropic overall motion of the protein into account. Experimental and simulated spectral densities are in reasonable general agreement for NH bond vectors, where the corresponding motions have converged within the simulation time. A sufficient sampling of the motions for NH bonds within flexible parts of the protein requires a longer simulation time. The simulated spectral densities J(0) and J(omega(N)) are, on average, 4.5% and 16% lower than the experimental data; the corresponding numbers for the core residues are about 6%; the high-frequency spectral densities J(0.87omega(H)) are lower by, on average, 16% (21% for the core). The simulated order parameters, S(2), are also lower, although the overall disagreement between the simulation and experiment is less pronounced: 1% for all residues and 6% for the core. The observed systematic decrease of simulated spectral density and the order parameters compared to the experimental data can be partially attributed to the ultrafast librational motion of the NH bonds with respect to their peptide plane, which was analyzed in detail. This systematic difference is most pronounced for J(0.87omega(H)), which appears to be most sensitive to the slow, subnanosecond time scale of internal motion, whereas J(0) and J(omega(N)) are dominated by the overall rotational tumbling of the protein. Similar discrepancies are observed between the experimentally measured (15)N relaxation parameters (R(1), R(2), NOE) and their values calculated from the simulated spectral densities. The analysis of spectral densities provides additional information regarding the comparison of the simulated and experimental data, not available from the model-free analysis.

Intramolecular dynamics and conformational transition in proteins studied by biophysical labelling methods. Common and specific features of proteins from thermophylic micro-organisms

A general survey is carried out on the theoretical grounds for methods of spin, luminescence and Mössbauer labels, as well as their application in the study of protein intramolecular dynamics. When combined, these methods allow the protein dynamics to be investigated within a wide range of correlation times (tau c = 10(2) - 10(-10) s) and amplitudes. The purposeful application of the methods to various proteins at different temperatures (30-330 K), water content, substrate addition, etc., revealed a number of dynamical processes and conformational transitions in proteins. The experiments indicated correlations between the local segmental mobility of protein globules in a nanosecond temporal scale and biochemical reactions, such as long-distance electron transfer, hydrolysis and photoreactions. The biophysical labelling methods results were analysed together with the data on dynamics obtained using complementary physico-chemical methods and theoretical calculations. Special emphasis is given to recent results on proteins from thermophylic micro-organisms. The mechanisms of protein intramolecular dynamics and their role in the stability and functions of proteins and enzymes are discussed.

Applications of NMR in drug discovery

In the half-century since its discovery, nuclear magnetic resonance (NMR) has become the single most powerful form of spectroscopy in both chemistry and structural biology. The dramatic technical advances over the past 10-15 years, which continue apace, have markedly increased the range of applications for NMR in the study of protein-ligand interactions. These form the basis for its most exciting uses in the drug discovery process, which range from the simple identification of whether a compound (or a component of a mixture) binds to a given protein, through to the determination of the full three-dimensional structure of the complex, with all the information this yields for structure-based drug design.

Efficient sampling in collective coordinate space

Collective motions in biological macromolecules have been shown to be important for function. The most important collective motions occur on slow time scales, which poses a sampling problem in dynamic simulation of biomolecules. We present a novel method for efficient conformational sampling. The method combines the simulation of an ensemble of concurrent trajectories with restraints acting on the ensemble of structures as a whole. Two properties of the ensemble may be restrained: (i) the variance of the ensemble and (ii) the average position of the ensemble. Both properties are defined in a subspace of collective coordinate space spanned by an arbitrary number of modes. We show that weak restraints on the ensemble variance suffice for an increase in sampling efficiency along soft modes by two orders of magnitudes. The resulting trajectories exhibit virtually the same structural quality as trajectories generated by restraint-free-molecular dynamics simulation, as judged by standard structure validation tools. The method is used to probe the resistance of a structure against conformational changes along collective modes and clearly distinguishes soft from stiff modes. Further applications are discussed. Proteins 2000;39:82-88.

Backbone dynamics and refined solution structure of the N-terminal domain of DNA polymerase beta. Correlation with DNA binding and dRP lyase activity

Mammalian DNA polymerase beta functions in the base excision DNA repair pathway filling in short patches (1-5 nt) in damaged DNA and removing deoxyribose 5'-phosphate from the 5'-side of damaged DNA. The backbone dynamics and the refined solution structure of the N-terminal domain of beta-Pol have been characterized in order to establish the potential contribution(s) of backbone motion to the DNA binding and deoxyribose 5'-phosphate lyase function of this domain. The N-terminal domain is formed from four helices packed as two antiparallel pairs with a 60 degrees crossing between the pairs. The RMSD of the NMR conformers (residues 13-80) is 0.37 A for the backbone heavy atoms and 0.78 A for all heavy atoms. NMR characterization of the binding site(s) for a ssDNA-5mer, ssDNA-8mer, ssDNA-9mer, and dsDNA-12mer shows a consensus surface for the binding of these various DNA oligomers, that surrounds and includes the deoxyribose 5'-phosphate lyase active site region. Connection segments between helices 1 and 2 and between helices 3 and 4 each contribute to DNA binding. Helix-3-turn-helix-4 forms a helix-hairpin-helix motif. The highly conserved hairpin sequence (LPGVG) displays a significant degree of picosecond time-scale motion within the backbone, that is possibly important for DNA binding at the phosphodiester backbone. An Omega-loop connecting helices 1 and 2 and helix-2 itself display significant exchange contributions (R(ex)) at the backbone amides due to apparent conformational type motion on a millisecond time-scale. This motion is likely important in allowing the Omega-loop and helix-2 to shift toward, and productively interact with, gapped DNA. The deoxyribose 5'-phosphate lyase catalytic residues that include K72 which forms the Schiff's base, Y39 which is postulated to promote proton transfer to the aldehyde, and K35 which assists in phosphate elimination, show highly restricted backbone motion. H34, which apparently participates in detection of the abasic site hole and assists in the opening of the hemiacetal, shows conformational exchange.

Pleckstrin homology domains of tec family protein kinases

Pleckstrin homology (PH) domains have been shown to be involved in different interactions, including binding to inositol compounds, protein kinase C isoforms, and heterotrimeric G proteins. In some cases, the most important function of PH domains is transient localisation of proteins to membranes, where they can interact with their partners. Tec family protein tyrosine kinases contain a PH domain. In Btk, also PH domain mutations lead into an immunodeficiency, X-linked agammaglobulinemia (XLA). A new disease-causing mutation was identified in the PH domain. The structures for the PH domains of Bmx, Itk, and Tec were modelled based on Btk structure. The domains seem to have similar scaffolding and electrostatic polarisation but to have some differences in the binding regions. The models provide new insight into the specificity, function, and regulation of Tec family kinases.