Search

Recent technological advances in NMR spectroscopy have alleviated the size limitations for the determination of biomolecular structures in solution. At the same time, novel NMR parameters such as residual dipolar couplings are providing greater accuracy. As this review shows, the structures of protein-protein and protein-nucleic acid complexes up to 50 kDa can now be accurately determined. Although de novo structure determination still requires considerable effort, information on interaction surfaces from chemical shift perturbations is much easier to obtain. Advances in modelling and data-driven docking procedures allow this information to be used for determining approximate structures of biomolecular complexes. As a result, a wealth of information has become available on the way in which proteins interact with other biomolecules. Of particular interest is the fact that these NMR-based methods can be applied to weak and transient protein-protein complexes that are difficult to study by other structural methods.

We present two new databases of NMR-derived distance and dihedral angle restraints: the Database Of Converted Restraints (DOCR) and the Filtered Restraints Database (FRED). These databases currently correspond to 545 proteins with NMR structures deposited in the Protein Databank (PDB). The criteria for inclusion were that these should be unique, monomeric proteins with author-provided experimental NMR data and coordinates available from the PDB capable of being parsed and prepared in a consistent manner. The Wattos program was used to parse the files, and the CcpNmr FormatConverter program was used to prepare them semi-automatically. New modules, including a new implementation of Aqua in the BioMagResBank (BMRB) software Wattos were used to analyze the sets of distance restraints (DRs) for inconsistencies, redundancies, NOE completeness, classification and violations with respect to the original coordinates. Restraints that could not be associated with a known nomenclature were flagged. The coordinates of hydrogen atoms were recalculated from the positions of heavy atoms to allow for a full restraint analysis. The DOCR database contains restraint and coordinate data that is made consistent with each other and with IUPAC conventions. The FRED database is based on the DOCR data but is filtered for use by test calculation protocols and longitudinal analyses and validations. These two databases are available from websites of the BMRB and the Macromolecular Structure Database (MSD) in various formats: NMR-STAR, CCPN XML, and in formats suitable for direct use in the software packages CNS and CYANA.

We present a CPU efficient protocol for refinement of protein structures in a thin layer of explicit solvent and energy parameters with completely revised dihedral angle terms. Our approach is suitable for protein structures determined by theoretical (e.g., homology modeling or threading) or experimental methods (e.g., NMR). In contrast to other recently proposed refinement protocols, we put a strong emphasis on consistency with widely accepted covalent parameters and computational efficiency. We illustrate the method for NMR structure calculations of three proteins: interleukin- 4, ubiquitin, and crambin. We show a comparison of their structure ensembles before and after refinement in water with and without a force field energy term for the dihedral angles; crambin was also refined in DMSO. Our results demonstrate the significant improvement of structure quality by a short refinement in a thin layer of solvent. Further, they show that a dihedral angle energy term in the force field is beneficial for structure calculation and refinement. We discuss the optimal weight for the energy constant for the backbone angle omega and include an extensive discussion of meaning and relevance of the calculated validation criteria, in particular root mean square Z scores for covalent parameters such as bond lengths. Proteins 2003;50: 496–506.

Nisin is an example of type-A lantibiotics that contain cyclic lanthionine rings and unusual dehydrated amino acids. Among the numerous pore-forming antimicrobial peptides, type-A lantibiotics form an unique family of post-translationally modified peptides. Via the recognition of cell wall precursor lipid II, nisin has the capacity to form pores against Gram-positive bacteria with an extremely high activity in the nanomolar (nM) range. Here we report a high-resolution NMR spectroscopy study of nisin/lipid II interactions in SDS micelles as a model membrane system in order to elucidate the mechanism of molecular recognition at residue level. The binding to lipid II was studied through (15)N-(1)H HSQC titration, backbone amide proton temperature coefficient analysis, and heteronuclear (15)N[(1)H]-NOE relaxation dynamics experiments. Upon the addition of lipid II, significant changes were monitored in the N-terminal part of nisin. An extremely low amide proton temperature coefficient (Delta delta/Delta T) was found for the amide proton of Ala3 (> -0.1 ppb/K) in the complex form. This suggests tight hydrogen bonding and/or isolation from the bulk solvent for this residue. Large chemical shift perturbations were also observed in the first two rings. In contrast, the C-terminal part of nisin was almost unaffected. This part of the molecule remains flexible and solvent-exposed. On the basis of our results, a multistep pore-forming mechanism is proposed. The N-terminal part of nisin first binds to lipid II, and a subsequent structural rearrangement takes place. The C-terminal part of nisin is possibly responsible for the activation of the pore formation. In light of the emerging antibiotic resistance problems, an understanding of the specific recognition mechanism of nisin with lipid II at the residue specific level may therefore aid in the development of novel antibiotics.

Proteins involved in signal transduction can usually be present in two states: an inactive and an active (signaling) state. In the case of photoreceptors such as PYP, it has been shown that the signaling state has a large degree of structural and dynamic disorder. Conventional structural NMR approaches present difficulties in describing such partially unfolded states. Owing to the disordered dynamical and transient nature of such states classical NOE-based information, when present, is sparse. Chemical shift changes upon partial unfolding can, however, be easily monitored from HSQC spectra. We show here that such states can be modeled by defining native-like inter-residue contacts for those residues that do not shift significantly upon partial unfolding. The feasibility of this approach is demonstrated using lysozyme as a test case and applied to model the partially unfolded signaling state (pB) of a truncated form of the photoactive yellow protein for which a ‘‘classical’’ NOE-based structure is available for validation. This approach should be generally applicable to systems in which part of the structure remains in a well-defined native-like conformation.

Contains fulltext : 57416.pdf (Publisher’s version ) (Closed access) Several studies have shown that biomolecular NMR structures are often of lower quality when compared to crystal structures, and consequently they are often excluded from structural analyses. We present a publicly available database of re-refined NMR structures, exhibiting significantly improved quality. This database (available at http://www.cmbi.kun.nl/dress/) presents a uniformly refined and validated set of structural models that improves the value of these NMR structures as input for experimental and theoretical studies in many fields of research.

The retinoid X receptor (RXR) is a prominent member of the nuclear receptor family of ligand-inducible transcription factors. Many proteins of this family exert their function as heterodimers with RXR as a common upstream partner. Studies of the DNA-binding domains of several nuclear receptors reveal differences in structure and dynamics, both between the different proteins and between the free- and DNA-bound receptor DBDs. We investigated the differences in dynamics between RXR free in solution and in complex with a 14 base-pair oligonucleotide, using (1)H and (15)N relaxation studies. Nano- to picosecond dynamics were probed on (15)N, employing Lipari-Szabo analysis with an axially symmetric tumbling model to estimate the exchange contributions to the transverse relaxation rates. Furthermore, milli- to microsecond dynamics were estimated qualitatively for (1)H and (15)N, using CPMG-HSQC and CPMG-T(2) measurements with differential pulse spacing. RXR shows hardly any nano- to picosecond time-scale internal motion. Upon DNA binding, the order parameters show a tiny increase. Dynamics in the milli- to microsecond time scale is more prevalent. It is localized in the first and second zinc fingers of the free RXR. Upon DNA-binding, exchange associated with specific/aspecific DNA-binding of RXR is observed throughout the sequence, whereas conformational flexibility of the D-box and the second zinc finger of RXR is greatly reduced. Since this DNA-binding induced folding transition occurs remote from the DNA in a region which is involved in protein-protein interactions, it may very well be related to the cooperativity of dimeric DNA binding.

Experimental autoimmune encephalomyelitis can be induced in susceptible animals by immunodominant determinants of myelin basic protein (MBP). To characterize the molecular features of antigenic sites important for designing experimental autoimmune encephalomyelitis suppressing molecules, we report structural studies, based on NMR experimental data in conjunction with molecular dynamic simulations, of the potent linear dodecapeptide epitope of guinea pig MBP, Gln74-Lys75-Ser76-Gln77-Arg78-Ser79-Gln80-Asp81-Glu82-Asn83-Pro84-Val85 [MBP(74-85)], and its antagonist analogue Ala81MBP(74-85). The two peptides were studied in both water and Me(2)SO in order to mimic solvent-dependent structural changes in MBP. The agonist MBP(74-85) adopts a compact conformation because of electrostatic interactions of Arg78 with the side chains of Asp81 and Glu82. Arg78 is 'locked' in a well-defined conformation, perpendicular to the peptide backbone which is practically solvent independent. These electrostatic interactions are, however, absent from the antagonist Ala81MBP(74-85), resulting in great flexibility of the side chain of Arg78. Sequence alignment of the two analogues with several species of MBP suggests a critical role for the positively charged residue Arg78, firstly, in the stabilization of the local microdomains (epitopes) of the integral protein, and secondly, in a number of post-translational modifications relevant to multiple sclerosis, such as the conversion of charged arginine residues to uncharged citrullines. Flexible docking calculations on the binding of the MBP(74-85) antigen to the MHC class II receptor site I-A(u) using haddock indicate that Gln74, Ser76 and Ser79 are MHC II anchor residues. Lys75, Arg78 and Asp81 are prominent, solvent-exposed residues and, thus, may be of importance in the formation of the trimolecular T-cell receptor-MBP(74-85)-MHC II complex.

Specific protein–protein interactions play crucial roles in the regulation of any biological process. Recently, a new protein–protein interaction domain termed PB1 (Phox and Bem1) was identified, which is conserved throughout evolution and present in diverse proteins functioning in signal transduction, cell polarity and survival. Here, we investigated the structure and molecular interactions of the PB1 heterodimer complex composed of the PB1 domains of the yeast proteins Bem1 and Cdc24. A structural model of the Cdc24 PB1 was built by homology modeling and molecular dynamics simulations, and experimentally validated by 15N nuclear Overhauser effect spectroscopy (NOESY)–heteronuclear single quantum coherence (HSQC) analysis. Residues at the interface of the complex for both proteins were identified by NMR titration experiments. A model of the heterodimer was obtained by docking of the two PB1 domains with HADDOCK, which applies ambiguous interaction restraints on residues at the interface to drive the docking procedure. The refined model was validated by site-directed mutagenesis on both Bem1 and Cdc24. Finally, the docking was repeated from the newly published NMR structure of Cdc24, allowing us to assess the performance of homology-based docking. Our results provide insight into the molecular structure of the Bem1–Cdc24 PB1-mediated heterodimer, which allowed identification of critical residues at the binding interface.

The high sensitivity of modern NMR instrumentation, in combination with full deuteration, enabled the measurement of long-range NOEs between amide protons in a fully deuterated protein corresponding to distances up to 8 Å. These are beyond the limit normally observed in protonated samples. Such long-distance NOEs could be observed using long mixing times, which became possible due to reduced spin diffusion and T1 relaxation of the amide protons in the fully deuterated sample. This information was used in combination with secondary structure restraints derived from secondary chemical shifts for structure calculations. With these backbone amide proton NOEs only, a unique fold could be obtained with positional root mean square deviations from the average of 1.30 and 2.25 Å for backbone and heavy atoms, respectively. Despite the low density of restraints, no mirror image problems were observed. Addition of sidechain NOE information increased the precision of the ensemble and in particular of the core packing. The structures obtained in this way were close to the published crystal structure. NOE completeness analysis revealed that the cumulative completeness is still more than 80% for an 8.0 Å cut-off distance.

Interaction of regulatory DNA binding proteins with their target sites is usually preceded by binding to nonspecific DNA. This speeds up the search for the target site by several orders of magnitude. We report the solution structure and dynamics of the complex of a dimeric lac repressor DNA binding domain with nonspecific DNA. The same set of residues can switch roles from a purely electrostatic interaction with the DNA backbone in the nonspecific complex to a highly specific binding mode with the base pairs of the cognate operator sequence. The protein-DNA interface of the nonspecific complex is flexible on biologically relevant time scales that may assist in the rapid and efficient finding of the target site.

The glucocorticoid receptor (GR) is a hormone dependent nuclear receptor that regulates gene transcription when bound to the glucocorticoid response element (GRE). The GRE acts as an allosteric effector, inducing a structural change in the GR DNA-binding domain (GR DBD) upon binding, thereby switching the GR to an active conformation. A similar conformational change can be induced by two single point mutations: Ser459Ala and Pro493Arg . Structural and dynamical aspects of the conformational switch have been investigated by molecular dynamics simulations. Our results indicate that these two mutants, which share a similar phenotype, exert their action at a structural level through different mechanisms. In the Arg 493 mutant, the D-loop and the second helix are stabilized in a conformation that preforms the protein–protein dimer interface. In the Ala 459 mutant, the structurally important hydrogen bond between Arg 496 and Ser 459 is missing, which leads to a core reorganization and a reorientation of the second helical region. Although remote, both in sequence and three dimensional structure, these two mutations induce structural changes that are ultimately reflected in similar regions of the GR DBD structure, namely the D-loop and the short second helical region. These correspond to hot area of the GR DBD that are important both for DNA-binding and for the proper formation of the protein–protein interface. The conformational rearrangements in these area are proposed to decrease unfavorable protein–DNA and protein–protein contacts and allow unspecific DNA-binding leading to the squelching phenotype of the mutants. The GR DBD can thus exist in two states, a transcriptionally active and an inactive state. Switching between these states can be accomplished either by GRE binding or by the described mutations.

Herpes simplex virion protein 16 (VP16) contains two strong activation regions that can independently and cooperatively activate transcription in vivo. We have identified the regions and residues involved in the interaction with the human transcriptional coactivator positive cofactor 4 (PC4) and the general transcription factor TFIIB. NMR and biochemical experiments revealed that both VP16 activation regions are required for the interaction and undergo a conformational transition from random coil to alpha-helix upon binding to its target PC4. The interaction is strongly electrostatically driven and the binding to PC4 is enhanced by the presence of its amino-terminal domain. We propose models for binding of VP16 to the core domains of PC4 and TFIIB that are based on two independent docking approaches using NMR chemical shift changes observed in titration experiments. The models are consistent with results from site-directed mutagenesis and provide an explanation for the contribution of both acidic and hydrophobic residues for transcriptional activation by VP16. Both intrinsically unstructured activation domains are attracted to their interaction partner by electrostatic interactions, and adopt an alpha-helical conformation around the important hydrophobic residues. The models showed multiple distinct binding surfaces upon interaction with various partners, providing an explanation for the promiscuous properties, cooperativity, and the high activity of this activation domain.

Background: Lactose repressor protein (Lac) controls the expression of the lactose metabolic genes in Escherichia coli by binding to an operator sequence in the promoter of the lac operon. Binding of inducer molecules to the Lac core domain induces changes in tertiary structure that are propagated to the DNA-binding domain through the connecting hinge region, thereby reducing the affinity for the operator. Protein–protein and protein–DNA interactions involving the hinge region play a crucial role in the allosteric changes occurring upon induction, but have not, as yet, been analyzed in atomic detail. Results: We have used nuclear magnetic resonance (NMR) spectroscopy and restrained molecular dynamics (rMD) to determine the structure of the Lac repressor DNA-binding domain (headpeice 62; HP62) in complex with a symmetrized lac operator. Analysis of the structures reveals specific interactions between Lac repressor and DNA that were not found in previously investigated Lac repressor–DNA complexes. Important differences with the previously reported structures of the HP56–DNA complex were found in the loop following the helix-turn-helix (HTH) motif. The protein–protein and protein–DNA interactions involving the hinge region and the deformations in the DNA structure could be delineated in atomic detail. The structures were also used for comparison with the available crystallographic data on the Lac and Pur repressor–DNA complexes. Conclusions: The structures of the HP62–DNA complex provide the basis for a better understanding of the specific recognition in the Lac repressor–operator complex. In addition, the structural features of the hinge region provide detailed insight into the protein–protein and protein–DNA interactions responsible for the high affinity of the repressor for operator DNA.

Pax-8 is a transcription factor belonging to the PAX genes superfamily and its crucial role has been proven both in embryo and in the adult organism. Pax-8 activity is regulated via a redoxbased mechanism centered on the glutathionylation of specific cysteines in the N-terminal region (Cys45 and Cys57). These residues belong to a highly evolutionary conserved DNA binding site: the Paired Box (Prd) domain. Crystallographic protein-DNA complexes of the homologues Pax-6 and Pax-5 showed a bipartite Prd domain consisting of two helix-turn-helix (HTH) motifs separated by an extended linker region. Here, by means of nuclear magnetic resonance, we show for the first time that the HTH motifs are largely defined in the unbound Pax-8 Prd domain. Our findings contrast with previous induced fit models, in which Pax-8 is supposed to largely fold upon DNA binding. Importantly, our data provide the structural basis for the enhanced chemical reactivity of residues Cys45 and Cys57 and explain clinical missense mutations that are not obviously related to the DNA binding interface of the paired box domain. Finally, sequence conservation suggests that our findings could be a general feature of the Pax family transcription factors.

Sequence-specific protein-DNA interactions are responsible for the regulation of key biological functions such as replication of the genome, initiation of transcription, and repair of damaged DNA. All of these regulatory pathways are built on the foundation that proteins are able to bind selectively to a particular DNA site in the genome. The most challenging issue for specific protein-DNA recognition is that the target sequence is immersed in a huge molar excess of nonspecific DNA sequences, which also compete for the same protein. The distribution of protein between regulatory and nonspecific sequences determines the occupancy of the target site in the cell and, hence, the transcriptional activity of the corresponding promoter. The problem is much more acute than it is for other protein-ligand interactions, because obviously the protein cannot differentiate between these sites on the basis of their size or shape.1 For repressor molecules, such as the lactose (lac) repressor, the problem of finding the correct site is even more severe, as only a few copies ( 10) of the protein exist in the cell. Thus, to compete with the large excess of nonspecific genomic DNA, these proteins have evolved to bind with a very high specificity ratio ( 107).

Structural genomics, the determination of protein structures on a genome-wide scale, is still in its infancy for eukaryotes due to the number and size of their genes. Low protein expression and solubility of eukaryotic geneproducts are the major bottlenecks in high-throughput (HTP) recombinant protein production with the E. coli expression systems. To circumvent this problem we decided to focus on separate protein domains. We describe here a fast microtiterplate based, expression and solubility screening procedure, using a combination of in vitro and in vivo expression, and purification with nickel-NTA magnetic beads. All steps are optimized for automatic HTP processing using a liquid handling station. Furthermore, large-scale expression and protein purification conditions are optimized, permitting the purification of 24 protein samples per week. We further show that results obtained from the expression screening can be extrapolated to the production of protein samples for NMR. Starting with 81 cloned human protein domains, in vivo expression was detected in 54 cases, and from 28 of those milligrams of protein were purified. An informative HSQC spectrum was recorded for 18 proteins (22%), half of which were indicative of a folded protein. The success rate and quality of the HSQC spectra suggest that the domain approach holds promise for human proteins.