The domain architecture of the protozoan protein J-DNA-binding protein 1 suggests synergy between base J DNA binding and thymidine hydroxylase activity.

J-DNA-binding protein 1 (JBP1) contributes to the biosynthesis and maintenance of base J (ß-d-glucosyl-hydroxymethyluracil), an epigenetic modification of thymidine (T) confined to pathogenic protozoa such as Trypanosoma and Leishmania JBP1 has two known functional domains: an N-terminal T hydroxylase (TH) homologous to the 5-methylcytosine hydroxylase domain in TET proteins and a J-DNA-binding domain (JDBD) that resides in the middle of JBP1. Here, we show that removing JDBD from JBP1 results in a soluble protein (Δ-JDBD) with the N- and C-terminal regions tightly associated together in a well-ordered structure. We found that this Δ-JDBD domain retains TH activity in vitro but displays a 15-fold lower apparent rate of hydroxylation compared with JBP1. Small-angle X-ray scattering (SAXS) experiments on JBP1 and JDBD in the presence or absence of J-DNA and on Δ-JDBD enabled us to generate low-resolution three-dimensional models. We conclude that Δ-JDBD, and not the N-terminal region of JBP1 alone, is a distinct folding unit. Our SAXS-based model supports the notion that binding of JDBD specifically to J-DNA can facilitate T hydroxylation 12-14 bp downstream on the complementary strand of the J-recognition site. We postulate that insertion of the JDBD module into the Δ-JDBD scaffold during evolution provided a mechanism that synergized J recognition and T hydroxylation, ensuring inheritance of base J in specific sequence patterns following DNA replication in kinetoplastid parasites.

An allosteric binding site of the α7 nicotinic acetylcholine receptor revealed in a humanized acetylcholine-binding protein.

Nicotinic acetylcholine receptors (nAChRs) belong to the family of pentameric ligand-gated ion channels and mediate fast excitatory transmission in the central and peripheral nervous systems. Among the different existing receptor subtypes, the homomeric α7 nAChR has attracted considerable attention because of its possible implication in several neurological and psychiatric disorders, including cognitive decline associated with Alzheimer's disease or schizophrenia. Allosteric modulators of ligand-gated ion channels are of particular interest as therapeutic agents, as they modulate receptor activity without affecting normal fluctuations of synaptic neurotransmitter release. Here, we used X-ray crystallography and surface plasmon resonance spectroscopy of α7-acetylcholine-binding protein (AChBP), a humanized chimera of a snail AChBP, which has 71% sequence similarity with the extracellular ligand-binding domain of the human α7 nAChR, to investigate the structural determinants of allosteric modulation. We extended previous observations that an allosteric site located in the vestibule of the receptor offers an attractive target for receptor modulation. We introduced seven additional humanizing mutations in the vestibule-located binding site of AChBP to improve its suitability as a model for studying allosteric binding. Using a fragment-based screening approach, we uncovered an allosteric binding site located near the ß8-ß9 loop, which critically contributes to coupling ligand binding to channel opening in human α7 nAChR. This work expands our understanding of the topology of allosteric binding sites in AChBP and, by extrapolation, in the human α7 nAChR as determined by electrophysiology measurements. Our insights pave the way for drug design strategies targeting nAChRs involved in ion channel-mediated disorders.

A series of PDB-related databanks for everyday needs.

We present a series of databanks (http://swift.cmbi.ru.nl/gv/facilities/) that hold information that is computationally derived from Protein Data Bank (PDB) entries and that might augment macromolecular structure studies. These derived databanks run parallel to the PDB, i.e. they have one entry per PDB entry. Several of the well-established databanks such as HSSP, PDBREPORT and PDB_REDO have been updated and/or improved. The software that creates the DSSP databank, for example, has been rewritten to better cope with π-helices. A large number of databanks have been added to aid computational structural biology; some examples are lists of residues that make crystal contacts, lists of contacting residues using a series of contact definitions or lists of residue accessibilities. PDB files are not the optimal presentation of the underlying data for many studies. We therefore made a series of databanks that hold PDB files in an easier to use or more consistent representation. The BDB databank holds X-ray PDB files with consistently represented B-factors. We also added several visualization tools to aid the users of our databanks.

Detection of trans-cis flips and peptide-plane flips in protein structures.

A coordinate-based method is presented to detect peptide bonds that need correction either by a peptide-plane flip or by a trans-cis inversion of the peptide bond. When applied to the whole Protein Data Bank, the method predicts 4617 trans-cis flips and many thousands of hitherto unknown peptide-plane flips. A few examples are highlighted for which a correction of the peptide-plane geometry leads to a correction of the understanding of the structure-function relation. All data, including 1088 manually validated cases, are freely available and the method is available from a web server, a web-service interface and through WHAT_CHECK.

Data mining in the Life Sciences with Random Forest: a walk in the park or lost in the jungle?

In the Life Sciences 'omics' data is increasingly generated by different high-throughput technologies. Often only the integration of these data allows uncovering biological insights that can be experimentally validated or mechanistically modelled, i.e. sophisticated computational approaches are required to extract the complex non-linear trends present in omics data. Classification techniques allow training a model based on variables (e.g. SNPs in genetic association studies) to separate different classes (e.g. healthy subjects versus patients). Random Forest (RF) is a versatile classification algorithm suited for the analysis of these large data sets. In the Life Sciences, RF is popular because RF classification models have a high-prediction accuracy and provide information on importance of variables for classification. For omics data, variables or conditional relations between variables are typically important for a subset of samples of the same class. For example: within a class of cancer patients certain SNP combinations may be important for a subset of patients that have a specific subtype of cancer, but not important for a different subset of patients. These conditional relationships can in principle be uncovered from the data with RF as these are implicitly taken into account by the algorithm during the creation of the classification model. This review details some of the to the best of our knowledge rarely or never used RF properties that allow maximizing the biological insights that can be extracted from complex omics data sets using RF.

Structural and dynamic aspects of Ca2+ and Mg2+ binding of the regulatory domains of the Na+/Ca2+ exchanger.

Intracellular Ca2+ regulates the activity of the NCX (Na+/Ca2+ exchanger) through binding to the cytosolic CBD (Ca2+-binding domain) 1 and CBD2. In vitro studies of the structure and dynamics of CBD1 and CBD2, as well as studies of their kinetics and thermodynamics of Ca2+ binding, greatly enhanced our understanding of NCX regulation. We describe the fold of the CBDs in relation to other known structures and review Ca2+ binding of the different CBD variants from a structural perspective. We also report on new findings concerning Mg2+ binding to the CBDs and finally we discuss recent results on CBD1-CBD2 interdomain interactions.

The second Ca(2+)-binding domain of NCX1 binds Mg2+ with high affinity.

We report the effects of binding of Mg(2+) to the second Ca(2+)-binding domain (CBD2) of the sodium-calcium exchanger. CBD2 is known to bind two Ca(2+) ions using its Ca(2+)-binding sites I and II. Here, we show by nuclear magnetic resonance (NMR), circular dichroism, isothermal titration calorimetry, and mutagenesis that CBD2 also binds Mg(2+) at both sites, but with significantly different affinities. The results from Mg(2+)-Ca(2+) competition experiments show that Ca(2+) can replace Mg(2+) from site I, but not site II, and that Mg(2+) binding affects the affinity for Ca(2+). Furthermore, thermal unfolding circular dichroism data demonstrate that Mg(2+) binding stabilizes the domain. NMR chemical shift perturbations and (15)N relaxation data reveal that Mg(2+)-bound CBD2 adopts a state intermediate between the apo and fully Ca(2+)-loaded forms. Together, the data show that at physiological Mg(2+) concentrations CBD2 is loaded with Mg(2+) preferentially at site II, thereby stabilizing and structuring the domain and altering its affinity for Ca(2+).

On the complexity of Engh and Huber refinement restraints: the angle τ as example.

The Engh and Huber parameters for bond lengths and bond angles have been used uncontested in macromolecular structure refinement from 1991 until very recently, despite critical discussion of their ubiquitous validity by many authors. An extensive analysis of the backbone angle τ (N-C(α)-C) illustrates that the Engh and Huber parameters can indeed be improved and a recent study [Tronrud et al. (2010), Acta Cryst. D66, 834-842] confirms these ideas. However, the present study of τ shows that improving the Engh and Huber parameters will be considerably more complex than simply making the parameters a function of the backbone ϕ, ψ angles. Many other aspects, such as the cooperativity of hydrogen bonds, the bending of secondary-structure elements and a series of biophysical aspects of the 20 amino-acid types, will also need to be taken into account. Different sets of Engh and Huber parameters will be needed for conceptually different refinement programs.

Homology-based hydrogen bond information improves crystallographic structures in the PDB.

The Protein Data Bank (PDB) is the global archive for structural information on macromolecules, and a popular resource for researchers, teachers, and students, amassing more than one million unique users each year. Crystallographic structure models in the PDB (more than 100,000 entries) are optimized against the crystal diffraction data and geometrical restraints. This process of crystallographic refinement typically ignored hydrogen bond (H-bond) distances as a source of information. However, H-bond restraints can improve structures at low resolution where diffraction data are limited. To improve low-resolution structure refinement, we present methods for deriving H-bond information either globally from well-refined high-resolution structures from the PDB-REDO databank, or specifically from on-the-fly constructed sets of homologous high-resolution structures. Refinement incorporating HOmology DErived Restraints (HODER), improves geometrical quality and the fit to the diffraction data for many low-resolution structures. To make these improvements readily available to the general public, we applied our new algorithms to all crystallographic structures in the PDB: using massively parallel computing, we constructed a new instance of the PDB-REDO databank (https://pdb-redo.eu). This resource is useful for researchers to gain insight on individual structures, on specific protein families (as we demonstrate with examples), and on general features of protein structure using data mining approaches on a uniformly treated dataset.

New Biological Insights from Better Structure Models.

Structure validation is a key component of all steps in the structure determination process, from structure building, refinement, deposition, and evaluation all the way to post-deposition optimisation of structures in the Protein Data Bank (PDB) by re-refinement and re-building. Today, many aspects of protein structures are understood better than 10years ago, and combined with improved software and more computing power, the automated PDB_REDO procedure can significantly improve about 85% of all X-ray structures ever deposited in the PDB. We review structure validation, structure improvement, and a series of validation resources and facilities that give access to improved PDB files and to reports on the quality of the original and the improved structures. Post-deposition optimisation generally leads to improved protein structures and a series of examples will illustrate how that, in turn, leads to improved or even novel biological insights.

Validation and correction of Zn-CysxHisy complexes.

Many crystal structures in the Protein Data Bank contain zinc ions in a geometrically distorted tetrahedral complex with four Cys and/or His ligands. A method is presented to automatically validate and correct these zinc complexes. Analysis of the corrected zinc complexes shows that the average Zn-Cys distances and Cys-Zn-Cys angles are a function of the number of cysteines and histidines involved. The observed trends can be used to develop more context-sensitive targets for model validation and refinement.

BDB: databank of PDB files with consistent B-factors.

Protein structures available from the PDB contain for each atom the coordinates, the occupancy and the B-factor that indicates the mobility of the atom. The values that should represent B-factors can relate to atomic motions in different ways. We present here a databank in which all B-factors have been converted to the one, homogeneous representation that is most useful for protein engineering applications. The Databank of PDB files with consistent B-factors (BDB) is freely available through http://www.cmbi.umcn.nl/bdb/.