Novel strategies to overcome expression problems encountered with toxic proteins: application to the production of Lac repressor proteins for NMR studies

NMR studies of structural aspects of allosteric regulation by the Lac repressor requires overexpression and isotope labeling of the protein. The size of the repressor makes it a challenging target, putting constraints on both expression conditions and sample preparation methods to overcome problems associated with studies of larger proteins by NMR. We optimized protocols for the production of deuterated functionally active thermostable dimeric Lac repressor and its core domain mutants. The Lac repressor core domain has never been obtained as a recombinant protein, possibly due to the observed toxicity to the host cells. We overcame the core domain induced toxicity by co-expression of this domain with the full length Lac repressor, combined with a stringent control of culture conditions. Significant overexpression was only obtained if during all stages of pre-culturing the bacteria were kept in their exponential growth phase at low density. The sensitivity of NMR measurements is dramatically affected by buffer conditions; we therefore used a thermofluor buffer optimization screen to determine the optimal buffer conditions. The combined thermofluor and NMR screening method yielded thermostable fully functional Lac repressor domain samples suitable for high-resolution NMR studies. The optimized procedures to adapt Escherichia coli to growth in D2O, to overcome toxicity, and to optimize protein sample conditions provides a broad range of universally applicable techniques for production of larger proteins for NMR spectroscopy.

Specificity and affinity of Lac repressor for the auxiliary operators O2 and O3 are explained by the structures of their protein-DNA complexes

The structures of a dimeric mutant of the Lac repressor DNA-binding domain complexed with the auxiliary operators O2 and O3 have been determined using NMR spectroscopy and compared to the structures of the previously determined Lac-O1 and Lac-nonoperator complexes. Structural analysis of the Lac-O1 and Lac-O2 complexes shows highly similar structures with very similar numbers of specific and nonspecific contacts, in agreement with similar affinities for these two operators. The left monomer of the Lac repressor in the Lac-O3 complex retains most of these specific contacts. However, in the right half-site of the O3 operator, there is a significant loss of protein-DNA contacts, explaining the low affinity of the Lac repressor for the O3 operator. The binding mode in the right half-site resembles that of the nonspecific complex. In contrast to the Lac-nonoperator DNA complex where no hinge helices are formed, the stability of the hinge helices in the weak Lac-O3 complex is the same as in the Lac-O1 and Lac-O2 complexes, as judged from the results of hydrogen/deuterium experiments.

The solution structure of DNA-free Pax-8 paired box domain accounts for redox regulation of transcriptional activity in the pax protein family

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 redox-based 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.

E2-c-Cbl recognition is necessary but not sufficient for ubiquitination activity

The E2 ubiquitin-conjugating enzymes UbcH7 and UbcH5B both show specific binding to the RING (really interesting new gene) domain of the E3 ubiquitin-protein ligase c-Cbl, but UbcH7 hardly supports ubiquitination of c-Cbl and substrate in a reconstituted system. Here, we found that neither structural changes nor subtle differences in the E2-E3 interaction surface are possible explanations for the functional specificity of UbcH5B and UbcH7 in their interaction with c-Cbl. The quick transfer of ubiquitin from the UbcH5B-Ub thioester to c-Cbl or other ubiquitin acceptors suggests that UbcH5B might functionally be a relatively pliable E2 enzyme. In contrast, the UbcH7-Ub thioester is too stable to transfer ubiquitin under our assay conditions, indicating that UbcH7 might be a more specific E2 enzyme. Our results imply that the interaction specificity between c-Cbl and E2 is required but not sufficient for transfer of ubiquitin to potential targets.

Structural insight into the recognition of the H3K4me3 mark by the TFIID subunit TAF3

Trimethylation of lysine residue K4 of histone H3 (H3K4me3) strongly correlates with active promoters for RNA polymerase II-transcribed genes. Several reader proteins, including the basal transcription factor TFIID, for this nucleosomal mark have been identified. Its TAF3 subunit specifically binds the H3K4me3 mark via its conserved plant homeodomain (PHD) finger. Here, we report the solution structure of the TAF3-PHD finger and its complex with an H3K4me3 peptide. Using a combination of NMR, mutagenesis, and affinity measurements, we reveal the structural basis of binding affinity, methylation-state specificity, and crosstalk with asymmetric dimethylation of R2. A unique local structure rearrangement in the K4me3-binding pocket of TAF3 due to a conserved sequence insertion underscores the requirement for cation-pi interactions by two aromatic residues. Interference by asymmetric dimethylation of arginine 2 suggests that a H3R2/K4 "methyl-methyl" switch in the histone code dynamically regulates TFIID-promoter association.

The HhH domain of the human DNA repair protein XPF forms stable homodimers

The human XPF-ERCC1 protein complex plays an essential role in nucleotide excision repair by catalysing positioned nicking of a DNA strand at the 5' side of the damage. We have recently solved the structure of the heterodimeric complex of the C-terminal domains of XPF and ERCC1 (Tripsianes et al., Structure 2005;13:1849-1858). We found that this complex comprises a pseudo twofold symmetry axis and that the helix-hairpin-helix motif of ERCC1 is required for DNA binding, whereas the corresponding domain of XPF is functioning as a scaffold for complex formation with ERCC1. Despite the functional importance of heterodimerization, the C-terminal domain of XPF can also form homodimers in vitro. We here compare the stabilities of homodimeric and heterodimeric complexes of the C-terminal domains of XPF and ERCC1. The higher stability of the XPF HhH complexes under various experimental conditions, determined using CD and NMR spectroscopy and mass spectrometry, is well explained by the structural differences that exist between the HhH domains of the two complexes. The XPF HhH homodimer has a larger interaction interface, aromatic stacking interactions, and additional hydrogen bond contacts as compared to the XPF/ERCC1 HhH complex, which accounts for its higher stability.

The structural basis of the difference in sensitivity for PNGase F in the de-N-glycosylation of the native bovine pancreatic ribonucleases B and BS

In glycoanalysis protocols, N-glycans from glycoproteins are most frequently released with peptide- N (4)-( N-acetyl-beta-glucosaminyl)asparagine amidase F (PNGase F). As the enzyme is an amidase, it cleaves the NH-CO linkage between the Asn side chain and the Asn-bound GlcNAc residue. Usually, the enzyme has a low activity, or is not active at all, on native glycoproteins. A typical example is native bovine pancreatic ribonuclease B (RNase B) with oligomannose-type N-glycans at Asn-34. However, native RNase BS, generated by subtilisin digestion of native RNase B, which comprises amino acid residues 21-124 of RNase B, is sensitive to PNGase F digestion. The same holds for carboxymethylated RNase B (RNase B (cm)). In this study, NMR spectroscopy and molecular modeling have been used to explain the differences in PNGase F activity for native RNase B, native RNase BS, and RNase B (cm). NMR analysis combined with literature data clearly indicated that the N-glycan at Asn-34 is more mobile in RNase BS than in RNase B. MD simulations showed that the region around Asn-34 in RNase B is not very flexible, whereby the alpha-helix of the amino acid residues 1-20 has a stabilizing effect. In RNase BS, the alpha-helix formed by amino acid residues 23-32 is significantly more flexible. Using these data, the possibilities for complex formation of both RNase B and RNase BS with PNGase F were studied, and a model for the RNase BS-PNGase F complex is proposed.

High-level expression of biologically active glycoprotein hormones in Pichia pastoris strains--selection of strain GS115, and not X-33, for the production of biologically active N-glycosylated 15N-labeled phCG

The methylotrophic yeast Pichia pastoris is widely used for the production of recombinant glycoproteins. With the aim to generate biologically active ¹⁵N-labeled glycohormones for conformational studies focused on the unravelling of the NMR structures in solution, the P. pastoris strains GS115 and X-33 were explored for the expression of human chorionic gonadotropin (phCG) and human follicle-stimulating hormone (phFSH). In agreement with recent investigations on the N-glycosylation of phCG, produced in P. pastoris GS115, using ammonia/glycerol-methanol as nitrogen/carbon sources, the N-glycosylation pattern of phCG, synthesized using NH₄Cl/glucose–glycerol–methanol, comprised neutral and charged, phosphorylated high-mannose-type N-glycans (Man_8–15GlcNAc₂). However, the changed culturing protocol led to much higher amounts of glycoprotein material, which is of importance for an economical realistic approach of the aimed NMR research. In the context of these studies, attention was also paid to the site specific N-glycosylation in phCG produced in P. pastoris GS115. In contrast to the rather simple N-glycosylation pattern of phCG expressed in the GS115 strain, phCG and phFSH expressed in the X-33 strain revealed, besides neutral high-mannose-type N-glycans, also high concentrations of neutral hypermannose-type N-glycans (Man_up-to-30GlcNAc₂). The latter finding made the X-33 strain not very suitable for generating ¹⁵N-labeled material. Therefore, ¹⁵N-phCG was expressed in the GS115 strain using the new optimized protocol. The ¹⁵N-enrichment was evaluated by ¹⁵N-HSQC NMR spectroscopy and GLC-EI/MS. Circular dichroism studies indicated that ¹⁵N-phCG/GS115 had the same folding as urinary hCG. Furthermore, ¹⁵N-phCG/GS115 was found to be similar to the unlabeled protein in every respect as judged by radioimmunoassay, radioreceptor assays, and in vitro bioassays.

Structural dynamics in the activation of Epac

Epac1 is a cAMP-responsive exchange factor for the small G-protein Rap. It consists of a regulatory region containing a cyclic nucleotide binding (CNB) domain and a catalytic region that activates Rap. In the absence of cAMP, access of Rap to the catalytic site is blocked by the regulatory region. We analyzed the conformational states of the CNB domain in the absence and in the presence of cAMP and cAMP analogues by NMR spectroscopy, resulting in the first direct insights into the activation mechanism of Epac. We prove that the CNB domain exists in equilibrium between the inactive and the active conformation, which is shifted by binding of cAMP. cAMP binding results in conformational changes in both the ligand binding pocket and the outer helical segments. We used two different cAMP antagonists that block these successive changes to elucidate the steps of this process. Highlighting the role of dynamics, the superactivator 8-pCPT-2'-O-Me-cAMP induces similar conformational changes as cAMP but causes different internal mobility. The results reveal the critical elements of the CNB domain of Epac required for activation and highlight the role of dynamics in this process.

Analysis of the XPA and ssDNA-binding surfaces on the central domain of human ERCC1 reveals evidence for subfunctionalization

Human ERCC1/XPF is a structure-specific endonuclease involved in multiple DNA repair pathways. We present the solution structure of the non-catalytic ERCC1 central domain. Although this domain shows structural homology with the catalytically active XPF nuclease domain, functional investigation reveals a completely distinct function for the ERCC1 central domain by performing interactions with both XPA and single-stranded DNA. These interactions are non-competitive and can occur simultaneously through distinct interaction surfaces. Interestingly, the XPA binding by ERCC1 and the catalytic function of XPF are dependent on a structurally homologous region of the two proteins. Although these regions are strictly conserved in each protein family, amino acid composition and surface characteristics are distinct. We discuss the possibility that after XPF gene duplication, the redundant ERCC1 central domain acquired novel functions, thereby increasing the fidelity of eukaryotic DNA repair.

Interactions between the toxin Kid of the bacterial parD system and the antitoxins Kis and MazE

The proteins Kid and Kis are the toxin and antitoxin, respectively, encoded by the parD operon of Escherichia coli plasmid R1. Kis prevents the inhibition of E. coli cell growth caused by the RNA cleavage activity of Kid. Overproduction of MazE, the chromosome-encoded homologue of Kis, has been demonstrated to neutralize Kid toxicity to a certain extent in the absence of native Kis. Here, we show that a high structural similarity exists between these antitoxins, using NMR spectroscopy. We report about the interactions between Kid and Kis that are responsible for neutralization of Kid toxicity and enhance autoregulation of parD transcription. Native macromolecular mass spectrometry data demonstrate that Kid and Kis form multiple complexes. At Kis:Kid ratios equal to or exceeding 1:1, as found in vivo in a plasmid-containing cell, various complexes are present, ranging from Kid(2)-Kis(2) tetramer up to Kis(2)-Kid(2)-Kis(2)-Kid(2)-Kis(2) decamer. When Kid is in excess of Kis, corresponding to an in vivo situation immediately after loss of the plasmid, the Kid(2)-Kis(2)-Kid(2) heterohexamer is the most abundant species. NMR chemical shift and intensity perturbations in the (1)H (15)N HSQC spectra of Kid and Kis, observed when titrating the partner protein, show that the interaction sites of Kid and Kis resemble those within the previously reported MazF(2)-MazE(2)-MazF(2) complex. Furthermore, we demonstrate that Kid(2)-MazE(2) tetramers can be formed via weak interactions involving a limited part of the Kis-binding residues of Kid. The functional roles of the identified Kid-Kis and Kid-MazE interaction sites and complexes in toxin neutralization and repression of transcription are discussed.

On the role of aromatic side chains in the photoactivation of BLUF domains

BLUF (blue-light sensing using FAD) domain proteins are a novel group of blue-light sensing receptors found in many microorganisms. The role of the aromatic side chains Y21 and W104, which are in close vicinity to the FAD cofactor in the AppA BLUF domain from Rhodobacter sphaeroides, is investigated through the introduction of several amino acid substitutions at these positions. NMR spectroscopy indicated that in the W104F mutant, the local structure of the FAD binding pocket was not significantly perturbed as compared to that of the wild type. Time-resolved fluorescence and absorption spectroscopy was applied to explore the role of Y21 and W104 in AppA BLUF photochemistry. In the Y21 mutants, FADH*-W* radical pairs are transiently formed on a ps time scale and recombine to the ground state on a ns time scale. The W104F mutant shows a spectral evolution similar to that of wild type AppA but with an increased yield of signaling state formation. In the Y21F/W104F double mutant, all light-driven electron-transfer processes are abolished, and the FAD singlet excited-state evolves by intersystem crossing to the triplet state. Our results indicate that two competing light-driven electron-transfer pathways are available in BLUF domains: one productive pathway that involves electron transfer from the tyrosine, which leads to signaling state formation, and one nonproductive electron-transfer pathway from the tryptophan, which leads to deactivation and the effective lowering of the quantum yield of the signaling state formation. Our results are consistent with a photoactivation mechanism for BLUF domains where signaling state formation proceeds via light-driven electron and proton transfer from the conserved tyrosine to FAD, followed by a hydrogen-bond rearrangement and radical-pair recombination.

Modeling protein-protein complexes involved in the cytochrome C oxidase copper-delivery pathway

Proper assembly and function of cytochrome c oxidase, which catalyzes the reduction of O2 and generates the proton gradient driving ATP synthesis, depend on correct copper delivery and incorporation. Structural details about the protein-protein complexes involved in this process are still missing. We describe here models of four complexes along this pathway obtained by combining bioinformatics interface predictions with information-driven docking and discuss their relevance with respect to known and pathogenic mutations.

Interactions of Kid-Kis toxin-antitoxin complexes with the parD operator-promoter region of plasmid R1 are piloted by the Kis antitoxin and tuned by the stoichiometry of Kid-Kis oligomers

The parD operon of Escherichia coli plasmid R1 encodes a toxin–antitoxin system, which is involved in plasmid stabilization. The toxin Kid inhibits cell growth by RNA degradation and its action is neutralized by the formation of a tight complex with the antitoxin Kis. A fascinating but poorly understood aspect of the kid–kis system is its autoregulation at the transcriptional level. Using macromolecular (tandem) mass spectrometry and DNA binding assays, we here demonstrate that Kis pilots the interaction of the Kid–Kis complex in the parD regulatory region and that two discrete Kis-binding regions are present on parD. The data clearly show that only when the Kis concentration equals or exceeds the Kid concentration a strong cooperative effect exists between strong DNA binding and Kid₂–Kis₂–Kid₂–Kis₂ complex formation. We propose a model in which transcriptional repression of the parD operon is tuned by the relative molar ratio of the antitoxin and toxin proteins in solution. When the concentration of the toxin exceeds that of the antitoxin tight Kid₂–Kis₂–Kid₂ complexes are formed, which only neutralize the lethal activity of Kid. Upon increasing the Kis concentration, (Kid₂–Kis₂)_n complexes repress the kid–kis operon.

Structure and Function of Bacterial Kid-Kis and Related Toxin-Antitoxin Systems

Toxin-antitoxin systems were discovered as plasmid auxiliary maintenance cassettes. In recent years, an increasing amount of structural and functional information has become available about the proteins involved, allowing the understanding of bacterial cell growth inhibition by the toxins on a molecular level. A well-studied TA system is formed by the proteins Kid and Kis, encoded by the parD operon of the Escherichia coli plasmid R1. The toxicity of Kid has been related to its endoribonuclease activity, which is counteracted by binding of the antitoxin Kis at the proposed active site. In this review, the structural studies on the Kid-Kis system are compared to those of three related toxin-antitoxin systems: MazF-MazE, CcdB-CcdA and RelE-RelB.

Light-induced flipping of a conserved glutamine sidechain and its orientation in the AppA BLUF domain

The AppA BLUF domain is a blue light photoreceptor containing flavin. Conserved glutamine 63 is necessary for the photocycle of the protein, and its side chain has been proposed to flip in response to blue light illumination. Recently published crystal structures of AppA WT and the AppA mutant C20S describe contradictory conclusions regarding the orientation of the conserved glutamine 63 side chain in the dark. Here, we present evidence from NMR spectroscopy confirming light-induced flipping of the glutamine side chain to form a strong hydrogen bond between the glutamine 63 side chain carbonyl group and the tyrosine 21 side chain hydroxyl proton in the light-induced state. Our conclusions are consistent with published data from UV/vis absorbance and FTIR spectroscopy, as well as the crystal structure of AppA WT.

Solution structure of a chemosensory protein from the desert locust Schistocerca gregaria

Chemical stimuli, generally constituted by small volatile organic molecules, are extremely important for the survival of different insect species. In the course of evolution, insects have developed very sophisticated biochemical systems for the binding and the delivery of specific semiochemicals to their cognate membrane-bound receptors. Chemosensory proteins (CSPs) are a class of small soluble proteins present at high concentration in insect chemosensory organs; they are supposed to be involved in carrying the chemical messages from the environment to the chemosensory receptors. In this paper, we report on the solution structure of CSPsg4, a chemosensory protein from the desert locust Schistocerca gregaria, which is expressed in the antennae and other chemosensory organs. The 3D NMR structure revealed an overall fold consisting of six alpha-helices, spanning residues 13-18, 20-31, 40-54, 62-78, 80-90, and 97-103, connected by loops which in some cases show dihedral angles typical of beta-turns. As in the only other chemosensory protein whose structure has been solved so far, namely, CSP from the moth Mamestra brassicae, four helices are arranged to form a V-shaped motif; another helix runs across the two V's, and the last one is packed against the external face. Analysis of the tertiary structure evidenced multiple hydrophobic cavities which could be involved in ligand binding. In fact, incubation of the protein with a natural ligand, namely, oleamide, produced substantial changes to the NMR spectra, suggesting extensive conformational transitions upon ligand binding.

The nucleotide-binding site of bacterial translation initiation factor 2 (IF2) as a metabolic sensor

Translational initiation factor 2 (IF2) is a guanine nucleotide-binding protein that can bind guanosine 3',5'-(bis) diphosphate (ppGpp), an alarmone involved in stringent response in bacteria. In cells growing under optimal conditions, the GTP concentration is very high, and that of ppGpp very low. However, under stress conditions, the GTP concentration may decline by as much as 50%, and that of ppGpp can attain levels comparable to those of GTP. Here we show that IF2 binds ppGpp at the same nucleotide-binding site and with similar affinity as GTP. Thus, GTP and the alarmone ppGpp can be considered two alternative physiologically relevant IF2 ligands. ppGpp interferes with IF2-dependent initiation complex formation, severely inhibits initiation dipeptide formation, and blocks the initiation step of translation. Our data suggest that IF2 has the properties of a cellular metabolic sensor and regulator that oscillates between an active GTP-bound form under conditions allowing active protein syntheses and an inactive ppGpp-bound form when shortage of nutrients would be detrimental, if not accompanied by slackening of this synthesis.

Comparative NMR study on the impact of point mutations on protein stability of Pseudomonas mendocina lipase

In this work we compare the dynamics and conformational stability of Pseudomonas mendocina lipase enzyme and its F180P/S205G mutant that shows higher activity and stability for use in washing powders. Our NMR analyses indicate virtually identical structures but reveal remarkable differences in local dynamics, with striking correspondence between experimental data (i.e., (15)N relaxation and H/D exchange rates) and data from Molecular Dynamics simulations. While overall the cores of both proteins are very rigid on the pico- to nanosecond timescale and are largely protected from H/D exchange, the two point mutations stabilize helices alpha1, alpha4, and alpha5 and locally destabilize the H-bond network of the beta-sheet (beta7-beta9). In particular, it emerges that helix alpha5, undergoing some fast destabilizing motions (on the pico- to nanosecond timescale) in wild-type lipase, is substantially rigidified by the mutation of Phe180 for a proline at its N terminus. This observation could be explained by the release of some penalizing strain, as proline does not require any "N-capping" hydrogen bond acceptor in the i+3 position. The combined experimental and simulated data thus indicate that reduced molecular flexibility of the F180P/S205G mutant lipase underlies its increased stability, and thus reveals a correlation between microscopic dynamics and macroscopic thermodynamic properties. This could contribute to the observed altered enzyme activity, as may be inferred from recent studies linking enzyme kinetics to their local molecular dynamics.

Information-driven protein-DNA docking using HADDOCK: it is a matter of flexibility

Intrinsic flexibility of DNA has hampered the development of efficient protein−DNA docking methods. In this study we extend HADDOCK (High Ambiguity Driven DOCKing) [C. Dominguez, R. Boelens and A. M. J. J. Bonvin (2003) J. Am. Chem. Soc.125, 1731–1737] to explicitly deal with DNA flexibility. HADDOCK uses non-structural experimental data to drive the docking during a rigid-body energy minimization, and semi-flexible and water refinement stages. The latter allow for flexibility of all DNA nucleotides and the residues of the protein at the predicted interface. We evaluated our approach on the monomeric repressor−DNA complexes formed by bacteriophage 434 Cro, the Escherichia coli Lac headpiece and bacteriophage P22 Arc. Starting from unbound proteins and canonical B-DNA we correctly predict the correct spatial disposition of the complexes and the specific conformation of the DNA in the published complexes. This information is subsequently used to generate a library of pre-bent and twisted DNA structures that served as input for a second docking round. The resulting top ranking solutions exhibit high similarity to the published complexes in terms of root mean square deviations, intermolecular contacts and DNA conformation. Our two-stage docking method is thus able to successfully predict protein−DNA complexes from unbound constituents using non-structural experimental data to drive the docking.

Cooperative α-helix unfolding in a protein-DNA complex from hydrogen-deuterium exchange

We present experimental evidence for a cooperative unfolding transition of an alpha-helix in the lac repressor headpiece bound to a symmetric variant of the lac operator, as inferred from hydrogen-deuterium (H-D) exchange experiments monitored by NMR spectroscopy. In the EX1 limit, observed exchange rates become pH-independent and exclusively sensitive to local structure fluctuations that expose the amide proton HN to exchange. Close to this regime, we measured decay rates of individual backbone HN signals in D2O, and of their mutual HN-HN NOE by time-resolved two-dimensional (2D) NMR experiments. The data revealed correlated exchange at the center of the lac headpiece recognition helix, Val20-Val23, and suggested that the correlation breaks down at Val24, at the C terminus of the helix. A lower degree of correlation was observed for the exchange of Val9 and Ala10 at the center of helix 1, while no correlation was observed for Val38 and Glu39 at the center of helix 3. We conclude that HN exchange in the recognition helix and, to some extent, in helix 1 is a cooperative event involving the unfolding of these helices, whereas the HN exchange in helix 3 is dominated by random local structure fluctuations.

Direct Use of Unassigned Resonances in NMR Structure Calculations with Proxy Residues

We present a method that significantly enhances the robustness of (automated) NMR structure determination by allowing the NOE data corresponding to unassigned NMR resonances to be used directly in the calculations. The unassigned resonances are represented by additional atoms or groups of atoms that have no interaction with the regular protein atoms except through distance restraints. These so-called "proxy" residues can be used to generate NOE-based distance restraints in a similar fashion as for the assigned part of the protein. If sufficient NOE information is available, the restraints are expected to place the proxies at positions close to the correct atoms for the unassigned resonance, which can facilitate subsequent assignment. Convergence can be further improved by supplying additional information about the possible identities of the unassigned resonances. We have implemented this approach in the widely used automated assignment and structure calculation protocols ARIA and CANDID. We find that it significantly increases the robustness of structure calculations with regard to missing assignments and yields structures of higher quality. Our approach is still able to find correctly folded structures with up to 30% randomly missing resonance assignments, and even when only backbone and beta resonances are present! This should be of significant value to NMR-based structural proteomics initiatives.

The solution structure of the AppA BLUF domain: insight into the mechanism of light-induced signaling

The transcriptional antirepressor AppA from the photosynthetic bacterium Rhodobacter sphaeroides senses both the light climate and the intracellular redox state. Under aerobic conditions in the dark, AppA binds to and thereby blocks the function of PpsR, a transcriptional repressor. Absorption of a blue photon dissociates AppA from PpsR and allows the latter to repress photosynthesis gene expression. The N terminus of AppA contains sequence homology to flavin-containing photoreceptors that belong to the BLUF family. Structural and chemical aspects of signal transduction mediated by AppA are still largely unknown. Here we present NMR studies of the N-terminal flavin-binding BLUF domain of AppA. Its solution structure adopts an alpha/beta-sandwich fold with a beta alpha beta beta alpha beta beta topology, which represents a new flavin-binding fold. Considerable disorder is observed for residues near the chromophore due to conformational exchange. This disorder is observed both in the dark and in the light-induced signaling state of AppA. Furthermore, we detect light-induced structural changes in a patch of surface residues that provide a structural link between light absorption and signal-transduction events.

Combining NMR relaxation with chemical shift perturbation data to drive protein-protein docking

The modeling of biomolecular complexes by computational docking using the known structures of their constituents is developing rapidly to become a powerful tool in structural biology. It is especially useful in combination with even limited experimental information describing the interface. Here we demonstrate for the first time the use of diffusion anisotropy in combination with chemical shift perturbation data to drive protein-protein docking. For validation purposes we make use of simulated diffusion anisotropy data. Inclusion of this information, which can be derived from NMR relaxation rates and reports on the orientation of the components of a complex with respect to the rotational diffusion tensor, substantially improves the docking results.