Water Accessibility Refinement of the Extended Structure of KirBac1.1 in the Closed State

NMR structures of membrane proteins are often hampered by poor chemical shift dispersion and internal dynamics which limit resolved distance restraints. However, the ordering and topology of these systems can be defined with site-specific water or lipid proximity. Membrane protein water accessibility surface area is often investigated as a topological function via solid-state NMR. Here we leverage water-edited solid-state NMR measurements in simulated annealing calculations to refine a membrane protein structure. This is demonstrated on the inward rectifier K+ channel KirBac1.1 found in Burkholderia pseudomallei. KirBac1.1 is homologous to human Kir channels, sharing a nearly identical fold. Like many existing Kir channel crystal structures, the 1p7b crystal structure is incomplete, missing 85 out of 333 residues, including the N-terminus and C-terminus. We measure solid-state NMR water proximity information and use this for refinement of KirBac1.1 using the Xplor-NIH structure determination program. Along with predicted dihedral angles and sparse intra- and inter-subunit distances, we refined the residues 1–300 to atomic resolution. All structural quality metrics indicate these restraints are a powerful way forward to solve high quality structures of membrane proteins using NMR.

Mapping the electrostatic potential of the nucleosome acidic patch

The nucleosome surface contains an area with negative electrostatic potential known as the acidic patch, which functions as a binding platform for various proteins to regulate chromatin biology. The dense clustering of acidic residues may impact their effective pKa and thus the electronegativity of the acidic patch, which in turn could influence nucleosome-protein interactions. We here set out to determine the pKa values of residues in and around the acidic patch in the free H2A-H2B dimer using NMR spectroscopy. We present a refined solution structure of the H2A-H2B dimer based on intermolecular distance restraints, displaying a well-defined histone-fold core. We show that the conserved histidines H2B H46 and H106 that line the acidic patch have pKa of 5.9 and 6.5, respectively, and that most acidic patch carboxyl groups have pKa values well below 5.0. For H2A D89 we find strong evidence for an elevated pKa of 5.3. Our data establish that the acidic patch is highly negatively charged at physiological pH, while protonation of H2B H106 and H2B H46 at slightly acidic pH will reduce electronegativity. These results will be valuable to understand the impact of pH changes on nucleosome-protein interactions in vitro, in silico or in vivo.

Paramagnetic spin labeling of a bacterial DnaB helicase for solid-state NMR

Labeling of biomolecules with a paramagnetic probe for nuclear magnetic resonance (NMR) spectroscopy enables determining long-range distance restraints, which are otherwise not accessible by classically used dipolar coupling-based NMR approaches. Distance restraints derived from paramagnetic relaxation enhancements (PREs) can facilitate the structure determination of large proteins and protein complexes. We herein present the site-directed labeling of the large oligomeric bacterial DnaB helicase from Helicobacter pylori with cysteine-reactive maleimide tags carrying either a nitroxide radical or a lanthanide ion. The success of the labeling reaction was followed by quantitative continuous-wave electron paramagnetic resonance (EPR) experiments performed on the nitroxide-labeled protein. PREs were extracted site-specifically from 2D and 3D solid-state NMR spectra. A good agreement with predicted PRE values, derived by computational modeling of nitroxide and Gd3+ tags in the low-resolution DnaB crystal structure, was found. Comparison of experimental PREs and model-predicted spin label-nucleus distances indicated that the size of the "blind sphere" around the paramagnetic center, in which NMR resonances are not detected, is slightly larger for Gd3+ (~14 Å) than for nitroxide (~11 Å) in 13C-detected 2D spectra of DnaB. We also present Gd3+-Gd3+ dipolar electron-electron resonance EPR experiments on DnaB supporting the conclusion that DnaB was present as a hexameric assembly.

Rotating-Frame Overhauser Transfer via Long-Lived Coherences

Solution-state distance restraints for protein structure determination with Ångström-level resolution rely on through-space transfer of magnetization between nuclear spins. Such magnetization transfers, named Overhauser effects, occur via dipolar magnetic couplings. We demonstrate improvements in magnetization transfer using long-lived coherences (LLCs)—singlet-triplet superpositions that are antisymmetric with respect to spin-permutation within pairs of coupled magnetic nuclei—as the magnetization source. Magnetization transfers in the presence of radio-frequency irradiation, known as ‘rotating-frame’ Overhauser effects (ROEs), are predicted by theory to improve by the use of LLCs; calculations are matched by preliminary experiments herein. The LLC-ROE transfers were compared to the transmission of magnetization via classical transverse routes. Long-lived coherences accumulate magnetization on an external third proton, K, with transfer rates that depended on the tumbling regime. I,S →K transfers in the LLC configuration for (I,S) are anticipated to match, and then overcome, the same transfer rates in the classical configuration as the molecular rotational correlation times increase. Experimentally, we measured the LLC-ROE transfer in dipeptide AlaGly between aliphatic protons in different residues K = Ala − Hα and (I,S) = Gly − Hα1,2 over a distance dK,I,S = 2.3 Å. Based on spin dynamics calculations, we anticipate that, for such distances, a superior transfer of magnetization occurs using LLC-ROE compared to classical ROE at correlation times above τC=10 ns. The LLC-ROE effect shows potential for improving structural studies of large proteins and offering constraints of increased precision for high-affinity protein-ligand complexes in slow tumbling in the liquid state.

Distance-based Reconstruction of Protein Quaternary Structures from Inter-Chain Contacts

Predicting the quaternary structure of protein complex is an important problem. Inter-chain residue-residue contact prediction can provide useful information to guide the ab initio reconstruction of quaternary structures. However, few methods have been developed to build quaternary structures from predicted inter-chain contacts. Here, we introduce a gradient descent optimization algorithm (GD) to build quaternary structures of protein dimers utilizing inter-chain contacts as distance restraints. We evaluate GD on several datasets of homodimers and heterodimers using true or predicted contacts. GD consistently performs better than a simulated annealing method and a Markov Chain Monte Carlo simulation method. Using true inter-chain contacts as input, GD can reconstruct high-quality structural models for homodimers and heterodimers with average TM-score ranging from 0.92 to 0.99 and average interface root mean square distance (I-RMSD) from 0.72 Å to 1.64 Å. On a dataset of 115 homodimers, using predicted inter-chain contacts as input, the average TM-score of the structural models built by GD is 0.76. For 46% of the homodimers, high-quality structural models with TM-score >= 0.9 are reconstructed from predicted contacts. There is a strong correlation between the quality of the reconstructed models and the precision and recall of predicted contacts. If the precision or recall of predicted contacts is >20%, GD can reconstruct good models for most homodimers, indicating only a moderate precision or recall of inter-chain contact prediction is needed to build good structural models for most homodimers. Moreover, the accuracy of reconstructed models positively correlates with the contact density in dimers.

Methodology for rigorous modeling of protein conformational changes by Rosetta using DEER Distance Restraints

We describe an approach for integrating distance restraints from Double Electron-Electron Resonance (DEER) spectroscopy into Rosetta with the purpose of modeling alternative protein conformations from an initial experimental structure. Fundamental to this approach is a multilateration algorithm that harnesses sets of interconnected spin label pairs to identify optimal rotamer ensembles at each residue that fit the DEER decay in the time domain. Benchmarked relative to data analysis packages, the algorithm yields comparable distance distributions with the advantage that fitting the DEER decay and rotamer ensemble optimization are coupled. We demonstrate this approach by modeling the protonation-dependent transition of the multidrug transporter PfMATE to an inward facing conformation with a deviation to the experimental structure of less than 2Å Cα RMSD. By decreasing spin label rotamer entropy, this approach engenders more accurate Rosetta models that are also more closely clustered, thus setting the stage for more robust modeling of protein conformational changes.

NMR refinement and peptide folding using the GROMACS software

Nuclear magnetic resonance spectroscopy is used routinely for studying the three-dimensional structures and dynamics of proteins and nucleic acids. Structure determination is usually done by adding restraints based upon NMR data to a classical energy function and performing restrained molecular simulations. Here we report on the implementation of a script to extract NMR restraints from a NMR-STAR file and export it to the GROMACS software. With this package it is possible to model distance restraints, dihedral restraints and orientation restraints. The output from the script is validated by performing simulations with and without restraints, including the ab initio refinement of one peptide.

NMR Refinement and Peptide Folding Using the GROMACS Software

<div><div><div><p>Nuclear magnetic resonance spectroscopy is used routinely for studying the three-dimensional structures and dynamics of proteins. Structure determination is usually done by adding restraints based upon NMR data to a classical energy function and performing restrained molecular simulations. Here we report on the implementation of a script to extract NMR restraints from a NMR-STAR file and export it to the GROMACS software. With this package it is possible to model distance restraints, dihedral restraints and orientation restraints. The output from the script is validated by performing simulations with and without restraints, including the ab initio refinement of one peptide.</p></div></div></div>

NMR Refinement and Peptide Folding Using the GROMACS Software

<div><div><div><p>Nuclear magnetic resonance spectroscopy is used routinely for studying the three-dimensional structures and dynamics of proteins. Structure determination is usually done by adding restraints based upon NMR data to a classical energy function and performing restrained molecular simulations. Here we report on the implementation of a script to extract NMR restraints from a NMR-STAR file and export it to the GROMACS software. With this package it is possible to model distance restraints, dihedral restraints and orientation restraints. The output from the script is validated by performing simulations with and without restraints, including the ab initio refinement of one peptide.</p></div></div></div>