High-Resolution Hydrogen-Deuterium Protection Factors from Sparse Mass Spectrometry Data Validated by Nuclear Magnetic Resonance Measurements.

Experimental measurement of time-dependent spontaneous exchange of amide protons with deuterium of the solvent provides information on the structure and dynamical structural variation in proteins. Two experimental techniques are used to probe the exchange: NMR, which relies on different magnetic properties of hydrogen and deuterium, and MS, which exploits the change in mass due to deuteration. NMR provides residue-specific information, that is, the rate of exchange or, analogously, the protection factor (i.e., the unitless ratio between the rate of exchange for a completely unstructured state and the observed rate). MS provides information that is specific to peptides obtained by proteolytic digestion. The spatial resolution of HDX-MS measurements depends on the proteolytic pattern of the protein, the fragmentation method used, and the overlap between peptides. Different computational approaches have been proposed to extract residue-specific information from peptide-level HDX-MS measurements. Here, we demonstrate the advantages of a method recently proposed that exploits self-consistency and classifies the possible sets of protection factors into a finite number of alternative solutions compatible with experimental data. The degeneracy of the solutions can be reduced (or completely removed) by exploiting the additional information encoded in the shape of the isotopic envelopes. We show how sparse and noisy MS data can provide high-resolution protection factors that correlate with NMR measurements probing the same protein under the same conditions.

Partial Opening of Cytochrome P450cam (CYP101A1) Is Driven by Allostery and Putidaredoxin Binding.

Cytochrome P450cam (CYP101A1) catalyzes the regio- and stereo-specific 5-exo-hydroxylation of camphor via a multistep catalytic cycle that involves two-electron transfer steps, with an absolute requirement that the second electron be donated by the ferrodoxin, putidaredoxin (Pdx). Whether P450cam, once camphor has bound to the active site and the substrate entry channel has closed, opens up upon Pdx binding, during the second electron transfer step, or it remains closed is still a matter of debate. A potential allosteric site for camphor binding has been identified and postulated to play a role in the binding of Pdx. Here, we have revisited paramagnetic NMR spectroscopy data and determined a heterogeneous ensemble of structures that explains the data, provides a complete representation of the P450cam/Pdx complex in solution, and reconciles alternative hypotheses. The allosteric camphor binding site is always present, and the conformational changes induced by camphor binding to this site facilitates Pdx binding. We also determined that the state to which Pdx binds comprises an ensemble of structures that have features of both the open and closed state. These results demonstrate that there is a finely balanced interaction between allosteric camphor binding and the binding of Pdx at high camphor concentrations.

Estimating Constraints for Protection Factors from HDX-MS Data.

Hydrogen/deuterium exchange monitored by mass spectrometry is a promising technique for rapidly fingerprinting structural and dynamical properties of proteins. The time-dependent change in the mass of any fragment of the polypeptide chain depends uniquely on the rate of exchange of its amide hydrogens, but determining the latter from the former is generally not possible. Here, we show that, if time-resolved measurements are available for a number of overlapping peptides that cover the whole sequence, rate constants for each amide hydrogen exchange (or equivalently, their protection factors) may be extracted and the uniqueness of the solutions obtained depending on the degree of peptide overlap. However, in most cases, the solution is not unique, and multiple alternatives must be considered. We provide a statistical method that clusters the solutions to further reduce their number. Such analysis always provides meaningful constraints on protection factors and can be used in situations in which obtaining more refined experimental data is impractical. It also provides a systematic way to improve data collection strategies to obtain unambiguous information at single-residue level (e.g., for assessing protein structure predictions at atomistic level).

Farseer-NMR: automatic treatment, analysis and plotting of large, multi-variable NMR data.

We present Farseer-NMR ( https://git.io/vAueU ), a software package to treat, evaluate and combine NMR spectroscopic data from sets of protein-derived peaklists covering a range of experimental conditions. The combined advances in NMR and molecular biology enable the study of complex biomolecular systems such as flexible proteins or large multibody complexes, which display a strong and functionally relevant response to their environmental conditions, e.g. the presence of ligands, site-directed mutations, post translational modifications, molecular crowders or the chemical composition of the solution. These advances have created a growing need to analyse those systems' responses to multiple variables. The combined analysis of NMR peaklists from large and multivariable datasets has become a new bottleneck in the NMR analysis pipeline, whereby information-rich NMR-derived parameters have to be manually generated, which can be tedious, repetitive and prone to human error, or even unfeasible for very large datasets. There is a persistent gap in the development and distribution of software focused on peaklist treatment, analysis and representation, and specifically able to handle large multivariable datasets, which are becoming more commonplace. In this regard, Farseer-NMR aims to close this longstanding gap in the automated NMR user pipeline and, altogether, reduce the time burden of analysis of large sets of peaklists from days/weeks to seconds/minutes. We have implemented some of the most common, as well as new, routines for calculation of NMR parameters and several publication-quality plotting templates to improve NMR data representation. Farseer-NMR has been written entirely in Python and its modular code base enables facile extension.

A Novel Mechanism for Calmodulin-Dependent Inactivation of Transient Receptor Potential Vanilloid 6.

The paralogues TRPV5 and TRPV6 belong to the vanilloid subfamily of the transient receptor potential (TRP) superfamily of ion channels, and both play an important role in overall Ca2+ homeostasis. The functioning of the channels centers on a tightly controlled Ca2+-dependent feedback mechanism in which the direct binding of the universal Ca2+-binding protein calmodulin (CaM) to the channel's C-terminal tail is required for channel inactivation. We have investigated this interaction at the atomic level and propose that under basal cellular Ca2+ concentrations CaM is constitutively bound to the channel's C-tail via CaM C-lobe only contacts. When the cytosolic Ca2+ concentration increases charging the apo CaM N-lobe with Ca2+, the CaM:TRPV6 complex rearranges and the TRPV6 C-tail further engages the CaM N-lobe via a crucial interaction involving L707. In a cellular context, mutation of L707 significantly increased the rate of channel inactivation. Finally, we present a model for TRPV6 CaM-dependent inactivation, which involves a novel so-called "two-tail" mechanism whereby CaM bridges two TRPV6 monomers resulting in closure of the channel pore.

Methyl group assignment using pseudocontact shifts with PARAssign.

A new version of the program PARAssign has been evaluated for assignment of NMR resonances of the 76 methyl groups in leucines, isoleucines and valines in a 25 kDa protein, using only the structure of the protein and pseudocontact shifts (PCS) generated with a lanthanoid tag at up to three attachment sites. The number of reliable assignments depends strongly on two factors. The principle axes of the magnetic susceptibility tensors of the paramagnetic centers should not be parallel so as to avoid correlated PCS. Second, the fraction of resonances in the spectrum of a paramagnetic sample that can be paired with the diamagnetic counterparts is critical for the assignment. With the data from two tag positions a reliable assignment could be obtained for 60% of the methyl groups and for many of the remaining resonances the number of possible assignments is limited to two or three. With a single tag, reliable assignments can be obtained for methyl groups with large PCS near the tag. It is concluded that assignment of methyl group resonances by paramagnetic tagging can be particularly useful in combination with some additional data, such as from mutagenesis or NOE-based experiments. Approaches to yield the best assignment results with PCS generating tags are discussed.

CcpNmr AnalysisAssign: a flexible platform for integrated NMR analysis.

NMR spectroscopy is an indispensably powerful technique for the analysis of biomolecules under ambient conditions, both for structural- and functional studies. However, in practice the complexity of the technique has often frustrated its application by non-specialists. In this paper, we present CcpNmr version-3, the latest software release from the Collaborative Computational Project for NMR, for all aspects of NMR data analysis, including liquid- and solid-state NMR data. This software has been designed to be simple, functional and flexible, and aims to ensure that routine tasks can be performed in a straightforward manner. We have designed the software according to modern software engineering principles and leveraged the capabilities of modern graphics libraries to simplify a variety of data analysis tasks. We describe the process of backbone assignment as an example of the flexibility and simplicity of implementing workflows, as well as the toolkit used to create the necessary graphics for this workflow. The package can be downloaded from www.ccpn.ac.uk/v3-software/downloads and is freely available to all non-profit organisations.

Delicate conformational balance of the redox enzyme cytochrome P450cam.

The energy landscapes of proteins are highly complex and can be influenced by changes in physical and chemical conditions under which the protein is studied. The redox enzyme cytochrome P450cam undergoes a multistep catalytic cycle wherein two electrons are transferred to the heme group and the enzyme visits several conformational states. Using paramagnetic NMR spectroscopy with a lanthanoid tag, we show that the enzyme bound to its redox partner, putidaredoxin, is in a closed state at ambient temperature in solution. This result contrasts with recent crystal structures of the complex, which suggest that the enzyme opens up when bound to its partner. The closed state supports a model of catalysis in which the substrate is locked in the active site pocket and the enzyme acts as an insulator for the reactive intermediates of the reaction.

Structure calculation, refinement and validation using CcpNmr Analysis.

CcpNmr Analysis provides a streamlined pipeline for both NMR chemical shift assignment and structure determination of biological macromolecules. In addition, it encompasses tools to analyse the many additional experiments that make NMR such a pivotal technique for research into complex biological questions. This report describes how CcpNmr Analysis can seamlessly link together all of the tasks in the NMR structure-determination process. It details each of the stages from generating NMR restraints [distance, dihedral, hydrogen bonds and residual dipolar couplings (RDCs)], exporting these to and subsequently re-importing them from structure-calculation software (such as the programs CYANA or ARIA) and analysing and validating the results obtained from the structure calculation to, ultimately, the streamlined deposition of the completed assignments and the refined ensemble of structures into the PDBe repository. Until recently, such solution-structure determination by NMR has been quite a laborious task, requiring multiple stages and programs. However, with the new enhancements to CcpNmr Analysis described here, this process is now much more intuitive and efficient and less error-prone.

A two-armed lanthanoid-chelating paramagnetic NMR probe linked to proteins via thioether linkages.

Paramagnetic NMR probes provide valuable long-range structural information on proteins and protein complexes. A new, stable, two-armed lanthanoid probe is reported that can be attached to a protein site-specifically via chemically inert thioether linkages.

The structure of the cytochrome p450cam-putidaredoxin complex determined by paramagnetic NMR spectroscopy and crystallography.

Cytochrome P450cam catalyzes the hydroxylation of camphor in a complex process involving two electron transfers (ETs) from the iron-sulfur protein putidaredoxin. The enzymatic control of the successive steps of catalysis is critical for a highly efficient reaction. The injection of the successive electrons is part of the control system. To understand the molecular interactions between putidaredoxin and cytochrome P450cam, we determined the structure of the complex both in solution and in the crystal state. Paramagnetic NMR spectroscopy using lanthanide tags yielded 446 structural restraints that were used to determine the solution structure. An ensemble of 10 structures with an RMSD of 1.3Å was obtained. The crystal structure of the complex was solved, showing a position of putidaredoxin that is identical with the one in the solution structure. The NMR data further demonstrate the presence of a minor state or set of states of the complex in solution, which is attributed to the presence of an encounter complex. The structure of the major state shows a small binding interface and a metal-to-metal distance of 16Å, with two pathways that provide strong electronic coupling of the redox centers. The interpretation of these results is discussed in the context of ET. The structure indicates that the ET rate can be much faster than the reported value, suggesting that the process may be gated.

PARAssign--paramagnetic NMR assignments of protein nuclei on the basis of pseudocontact shifts.

The use of paramagnetic NMR data for the refinement of structures of proteins and protein complexes is widespread. However, the power of paramagnetism for protein assignment has not yet been fully exploited. PARAssign is software that uses pseudocontact shift data derived from several paramagnetic centers attached to the protein to obtain amide and methyl assignments. The ability of PARAssign to perform assignment when the positions of the paramagnetic centers are known and unknown is demonstrated. PARAssign has been tested using synthetic data for methyl assignment of a 47 kDa protein, and using both synthetic and experimental data for amide assignment of a 14 kDa protein. The complex fitting space involved in such an assignment procedure necessitates that good starting conditions are found, both regarding placement and strength of paramagnetic centers. These starting conditions are obtained through automated tensor placement and user-defined tensor parameters. The results presented herein demonstrate that PARAssign is able to successfully perform resonance assignment in large systems with a high degree of reliability. This software provides a method for obtaining the assignments of large systems, which may previously have been unassignable, by using 2D NMR spectral data and a known protein structure.

Small-molecule binding sites on proteins established by paramagnetic NMR spectroscopy.

Determining the three-dimensional structure of a small molecule-protein complex with weak affinity can be a significant challenge. We present a paramagnetic NMR method to determine intermolecular structure restraints based on pseudocontact shifts (PCSs). Since the ligand must be in fast exchange between free and bound states and the fraction bound can be as low as a few percent, the method is ideal for ligands with high micromolar to millimolar dissociation constants. Paramagnetic tags are attached, one at a time, in a well-defined way via two arms at several sites on the protein surface. The ligand PCSs were measured from simple 1D (1)H spectra and used as docking restraints. An independent confirmation of the complex structure was carried out using intermolecular NOEs. The results show that structures derived from these two approaches are similar. The best results are obtained if the magnetic susceptibility tensors of the tags are known, but it is demonstrated that with two-armed probes, the magnetic susceptibility tensor can be predicted with sufficient accuracy to provide a low-resolution model of the ligand orientation and the location of the binding site in the absence of isotope-labeled protein. This approach can facilitate fragment-based drug discovery in obtaining structural information on the initial fragment hits.