Structure and dynamics of Helicobacter pylori nickel-chaperone HypA: an integrated approach using NMR spectroscopy, functional assays and computational tools.

Helicobacter pylori HypA (HpHypA) is a metallochaperone necessary for maturation of [Ni,Fe]-hydrogenase and urease, the enzymes required for colonization and survival of H. pylori in the gastric mucosa. HpHypA contains a structural Zn(II) site and a unique Ni(II) binding site at the N-terminus. X-ray absorption spectra suggested that the Zn(II) coordination depends on pH and on the presence of Ni(II). This study was performed to investigate the structural properties of HpHypA as a function of pH and Ni(II) binding, using NMR spectroscopy combined with DFT and molecular dynamics calculations. The solution structure of apo,Zn-HpHypA, containing Zn(II) but devoid of Ni(II), was determined using 2D, 3D and 4D NMR spectroscopy. The structure suggests that a Ni-binding and a Zn-binding domain, joined through a short linker, could undergo mutual reorientation. This flexibility has no physiological effect on acid viability or urease maturation in H. pylori. Atomistic molecular dynamics simulations suggest that Ni(II) binding is important for the conformational stability of the N-terminal helix. NMR chemical shift perturbation analysis indicates that no structural changes occur in the Zn-binding domain upon addition of Ni(II) in the pH 6.3-7.2 range. The structure of the Ni(II) binding site was probed using 1H NMR spectroscopy experiments tailored to reveal hyperfine-shifted signals around the paramagnetic metal ion. On this basis, two possible models were derived using quantum-mechanical DFT calculations. The results provide a comprehensive picture of the Ni(II) mode to HpHypA, important to rationalize, at the molecular level, the functional interactions of this chaperone with its protein partners.

The Structural Basis of Calcium-Dependent Inactivation of the Transient Receptor Potential Vanilloid 5 Channel.

The transient receptor potential vanilloid channel subfamily member 5 (TRPV5) is a highly selective calcium ion channel predominately expressed in the kidney epithelium that plays an essential role in calcium reabsorption from renal infiltrate. In order to maintain Ca2+ homeostasis, TRPV5 possesses a tightly regulated negative feedback mechanism, where the ubiquitous Ca2+ binding protein calmodulin (CaM) directly binds to the intracellular TRPV5 C-terminus, thus regulating TRPV5. Here we report on the characterization of the TRPV5 C-terminal CaM binding site and its interaction with CaM at an atomistic level. We have solved the de novo solution structure of the TRPV5 C-terminus in complex with a CaM mutant, creating conditions that mimic the cellular basal Ca2+ state. We demonstrate that under these conditions the TRPV5 C-terminus is exclusively bound to the CaM C-lobe only, while it confers conformational freedom to the CaM N-lobe. We also show that at elevated calcium levels, additional interactions between the TRPV5 C-terminus and CaM N-lobe occur, resulting in formation of a tight 1:1 complex, effectively making the N-lobe the calcium sensor. Together, these data are consistent with and support the novel model for Ca2+/CaM-dependent inactivation of TRPV channels as proposed by Bate and co-workers [ Bate , N. , et al. ( 2018 ) Biochemistry , ( 57), DOI: 10.1021/acs.biochem.7b01286 ].

The solution structure of the MANEC-type domain from hepatocyte growth factor activator inhibitor-1 reveals an unexpected PAN/apple domain-type fold.

A decade ago, motif at N-terminus with eight-cysteines (MANEC) was defined as a new protein domain family. This domain is found exclusively at the N-terminus of >400 multi-domain type-1 transmembrane proteins from animals. Despite the large number of MANEC-containing proteins, only one has been characterized at the protein level: hepatocyte growth factor activator inhibitor-1 (HAI-1). HAI-1 is an essential protein, as knockout mice die in utero due to placental defects. HAI-1 is an inhibitor of matriptase, hepsin and hepatocyte growth factor (HGF) activator, all serine proteases with important roles in epithelial development, cell growth and homoeostasis. Dysregulation of these proteases has been causatively implicated in pathological conditions such as skin diseases and cancer. Detailed functional understanding of HAI-1 and other MANEC-containing proteins is hampered by the lack of structural information on MANEC. Although many MANEC sequences exist, sequence-based database searches fail to predict structural homology. In the present paper, we present the NMR solution structure of the MANEC domain from HAI-1, the first three-dimensional (3D) structure from the MANEC domain family. Unexpectedly, MANEC is a new subclass of the PAN/apple domain family, with its own unifying features, such as two additional disulfide bonds, two extended loop regions and additional α-helical elements. As shown for other PAN/apple domain-containing proteins, we propose a similar active role of the MANEC domain in intramolecular and intermolecular interactions. The structure provides a tool for the further elucidation of HAI-1 function as well as a reference for the study of other MANEC-containing proteins.

Overlapping transport and chaperone-binding functions within a bacterial twin-arginine signal peptide.

The twin-arginine translocation (Tat) pathway is a protein targeting system present in many prokaryotes. The physiological role of the Tat pathway is the transmembrane translocation of fully-folded proteins, which are targeted by N-terminal signal peptides bearing conserved SRRxFLK 'twin-arginine' amino acid motifs. In Escherichia coli the majority of Tat targeted proteins bind redox cofactors and it is important that only mature, cofactor-loaded precursors are presented for export. Cellular processes have been unearthed that sequence these events, for example the signal peptide of the periplasmic nitrate reductase (NapA) is bound by a cytoplasmic chaperone (NapD) that is thought to regulate assembly and export of the enzyme. In this work, genetic, biophysical and structural approaches were taken to dissect the interaction between NapD and the NapA signal peptide. A NapD binding epitope was identified towards the N-terminus of the signal peptide, which overlapped significantly with the twin-arginine targeting motif. NMR spectroscopy revealed that the signal peptide adopted a α-helical conformation when bound by NapD, and substitution of single residues within the NapA signal peptide was sufficient to disrupt the interaction. This work provides an increased level of understanding of signal peptide function on the bacterial Tat pathway.

Structural diversity in twin-arginine signal peptide-binding proteins.

The twin-arginine transport (Tat) system is dedicated to the translocation of folded proteins across the bacterial cytoplasmic membrane. Proteins are targeted to the Tat system by signal peptides containing a twin-arginine motif. In Escherichia coli, many Tat substrates bind redox-active cofactors in the cytoplasm before transport. Coordination of cofactor insertion with protein export involves a "Tat proofreading" process in which chaperones bind twin-arginine signal peptides, thus preventing premature export. The initial Tat signal-binding proteins described belonged to the TorD family, which are required for assembly of N- and S-oxide reductases. Here, we report that E. coli NapD is a Tat signal peptide-binding chaperone involved in biosynthesis of the Tat-dependent nitrate reductase NapA. NapD binds tightly and specifically to the NapA twin-arginine signal peptide and suppresses signal peptide translocation activity such that transport via the Tat pathway is retarded. High-resolution, heteronuclear, multidimensional NMR spectroscopy reveals the 3D solution structure of NapD. The chaperone adopts a ferredoxin-type fold, which is completely distinct from the TorD family. Thus, NapD represents a new family of twin-arginine signal-peptide-binding proteins.

Fast empirical pKa prediction by Ewald summation.

pK(a) calculations for macromolecules are normally performed by solving the Poisson-Boltzmann equation, accounting for the different dielectric constants of solvent and solute, as well as the ionic strength. Despite the large number of successful applications, there are some situations where the current algorithms are not suitable: (1) large scale, high-throughput analysis which requires calculations to be completed within a fraction of a second, e.g. when permanently monitoring pK(a) shifts during a molecular dynamics simulation; (2) prediction of pK(a)s in periodic boundaries, e.g. when reconstructing entire protein crystal unit cells from PDB files, including the correct protonation patterns at experimental pH. Such in silico crystals are needed by 'self-parameterizing' molecular dynamics force fields like YASARA YAMBER, that optimize their parameters while energy-minimizing high-resolution protein crystals. To address both problems, we define an empirical equation that expresses the pK(a) as a function of electrostatic potential, hydrogen bonds and accessible surface area. The electrostatic potential is evaluated by Ewald summation, which captures periodic crystal environments and the uncertainty in atom positions using Gaussian charge densities. The empirical proportionality constants are derived from 217 experimentally determined pK(a)s, and despite its simplicity, this pK(a) calculation method reaches a high overall jack-knifed accuracy, and is fast enough to be used during a molecular dynamics simulation. A reliable null-model to judge pK(a) prediction accuracies is also presented.

Traditional biomolecular structure determination by NMR spectroscopy allows for major errors.

One of the major goals of structural genomics projects is to determine the three-dimensional structure of representative members of as many different fold families as possible. Comparative modeling is expected to fill the remaining gaps by providing structural models of homologs of the experimentally determined proteins. However, for such an approach to be successful it is essential that the quality of the experimentally determined structures is adequate. In an attempt to build a homology model for the protein dynein light chain 2A (DLC2A) we found two potential templates, both experimentally determined nuclear magnetic resonance (NMR) structures originating from structural genomics efforts. Despite their high sequence identity (96%), the folds of the two structures are markedly different. This urged us to perform in-depth analyses of both structure ensembles and the deposited experimental data, the results of which clearly identify one of the two models as largely incorrect. Next, we analyzed the quality of a large set of recent NMR-derived structure ensembles originating from both structural genomics projects and individual structure determination groups. Unfortunately, a visual inspection of structures exhibiting lower quality scores than DLC2A reveals that the seriously flawed DLC2A structure is not an isolated incident. Overall, our results illustrate that the quality of NMR structures cannot be reliably evaluated using only traditional experimental input data and overall quality indicators as a reference and clearly demonstrate the urgent need for a tight integration of more sophisticated structure validation tools in NMR structure determination projects. In contrast to common methodologies where structures are typically evaluated as a whole, such tools should preferentially operate on a per-residue basis.

Definition of a new information-based per-residue quality parameter.

For biomolecular NMR structures typically only a poor correspondence is observed between statistics derived from the experimental input data and structural quality indicators obtained from the structure ensembles. Here, we investigate the relationship between the amount of available NMR data and structure quality. By generating datasets with a predetermined information content and evaluating the quality of the resulting structure ensembles we show that there is, in contrast to previous findings, a linear relation between the information contained in experimental data and structural quality. From this relation, a new quality parameter is derived that provides direct insight, on a per-residue basis, into the extent to which structural quality is governed by the experimental input data.

Solution structure of the second PDZ domain of the neuronal adaptor X11alpha and its interaction with the C-terminal peptide of the human copper chaperone for superoxide dismutase.

Protection against reactive oxygen species is provided by the copper containing enzyme superoxide dismutase 1 (SOD1). The copper chaperone CCS is responsible for copper insertion into apo-SOD1. This role is impaired by an interaction between the second PDZ domain (PDZ2alpha) of the neuronal adaptor protein X11alpha and the third domain of CCS (McLoughlin et al. (2001) J. Biol. Chem., 276, 9303-9307). The solution structure of the PDZ2alpha domain has been determined and the interaction with peptides derived from CCS has been explored. PDZ2alpha binds to the last four amino acids of the CCS protein (PAHL) with a dissociation constant of 91 +/- 2 microM. Peptide variants have been used to map the interaction areas on PDZ2alpha for each amino acid, showing an important role for the C-terminal leucine, in line with canonical PDZ-peptide interactions.

RECOORD: a recalculated coordinate database of 500+ proteins from the PDB using restraints from the BioMagResBank.

State-of-the-art methods based on CNS and CYANA were used to recalculate the nuclear magnetic resonance (NMR) solution structures of 500+ proteins for which coordinates and NMR restraints are available from the Protein Data Bank. Curated restraints were obtained from the BioMagResBank FRED database. Although the original NMR structures were determined by various methods, they all were recalculated by CNS and CYANA and refined subsequently by restrained molecular dynamics (CNS) in a hydrated environment. We present an extensive analysis of the results, in terms of various quality indicators generated by PROCHECK and WHAT_CHECK. On average, the quality indicators for packing and Ramachandran appearance moved one standard deviation closer to the mean of the reference database. The structural quality of the recalculated structures is discussed in relation to various parameters, including number of restraints per residue, NOE completeness and positional root mean square deviation (RMSD). Correlations between pairs of these quality indicators were generally low; for example, there is a weak correlation between the number of restraints per residue and the Ramachandran appearance according to WHAT_CHECK (r = 0.31). The set of recalculated coordinates constitutes a unified database of protein structures in which potential user- and software-dependent biases have been kept as small as possible. The database can be used by the structural biology community for further development of calculation protocols, validation tools, structure-based statistical approaches and modeling. The RECOORD database of recalculated structures is publicly available from http://www.ebi.ac.uk/msd/recoord.

Combining data from genomes, Y2H and 3D structure indicates that BolA is a reductase interacting with a glutaredoxin.

Genomes, functional genomics data and 3D structure reflect different aspects of protein function. Here, we combine these data to predict that BolA, a widely distributed protein family with unknown function, is a reductase that interacts with a glutaredoxin. Comparisons at the 3D structure level as well as at the sequence profile level indicate homology between BolA and OsmC, an enzyme that reduces organic peroxides. Complementary to this, comparative analyses of genomes and genomics data provide strong evidence of an interaction between BolA and the mono-thiol glutaredoxin family. The interaction between BolA and a mono-thiol glutaredoxin is of particular interest because BolA does not, in contrast to its homolog OsmC, have evolutionarily conserved cysteines to provide it with reducing equivalents. We propose that BolA uses the mono-thiol glutaredoxin as the source for these.

Quantitative evaluation of experimental NMR restraints.

Nuclear Overhauser effect (NOE) data are an indispensable source of structural information in biomolecular structure determination by NMR spectroscopy. The number and type of experimental restraints used in the structure calculation and the RMS deviation of the restraints are usually reported. We present a new method for quantifying the information contained in the experimental NMR restraints. The method is based on a description of the structure in distance space and concepts derived from information theory. It allows for an objective description of the amount of available experimental information, which we show to be related to the positional uncertainty of the NMR ensemble. The measure of information presented is not affected by redundancy in the experimental restraints. Using various examples, we show that the method successfully identifies the crucial restraints in a structure determination: those restraints that are both important and unique. Finally, we demonstrate that the method can detect a wider range of redundancy in experimental datasets when compared to currently available methods. Because our method describes the quantitative evaluation of experimental NMR restraints, we propose the acronym QUEEN.

The precision of NMR structure ensembles revisited.

Biomolecular structures provide the basis for many studies in research areas such as structure-based drug design and homology modeling. In order to use molecular coordinates it is important that they are reliable in terms of accurate description of the experimental data and in terms of the overall and local geometry. Besides these primary quality criteria an indication is needed for the uncertainty in the atomic coordinates that may arise from the dynamic behavior of the considered molecules as well as from experimental- and computational procedures. In contrast to the crystallographic B-factor, a good measure for the uncertainty in NMR-derived atomic coordinates is still not available. It has become clear in recent years that the widely used atomic Root Mean Square Deviation (RMSD), which is a measure for the precision of the data, overestimates the accuracy of NMR structure ensembles and therefore is a problematic measure for the uncertainty in the atomic coordinates. In this study we report a method that yields a more realistic estimate of the uncertainty in the atomic coordinates by maximizing the RMSD of an ensemble of structures, while maintaining the accordance with the experimentally derived data. The results indicate that the RMSD of most NMR structure ensembles can be significantly increased compromising neither geometric quality nor NMR data. This maximized RMSD therefore seems a better estimate of the true uncertainty in the atomic coordinates.