Automated protein structure calculation from NMR data

Current software is almost at the stage to permit completely automatic structure determination of small proteins of < 15 kDa, from NMR spectra to structure validation with minimal user interaction. This goal is welcome, as it makes structure calculation more objective and therefore more easily validated, without any loss in the quality of the structures generated. Moreover, it releases expert spectroscopists to carry out research that cannot be automated. It should not take much further effort to extend automation to ca 20 kDa. However, there are technological barriers to further automation, of which the biggest are identified as: routines for peak picking; adoption and sharing of a common framework for structure calculation, including the assembly of an automated and trusted package for structure validation; and sample preparation, particularly for larger proteins. These barriers should be the main target for development of methodology for protein structure determination, particularly by structural genomics consortia

Pressure-dependent 13C chemical shifts in proteins: Origins and applications

Pressure-dependent (13)C chemical shifts have been measured for aliphatic carbons in barnase and Protein G. Up to 200 MPa (2 kbar), most shift changes are linear, demonstrating pressure-independent compressibilities. CH(3), CH(2) and CH carbon shifts change on average by +0.23, -0.09 and -0.18 ppm, respectively, due to a combination of bond shortening and changes in bond angles, the latter matching one explanation for the gamma-gauche effect. In addition, there is a residue-specific component, arising from both local compression and conformational change. To assess the relative magnitudes of these effects, residue-specific shift changes for protein G were converted into structural restraints and used to calculate the change in structure with pressure, using a genetic algorithm to convert shift changes into dihedral angle restraints. The results demonstrate that residual (13)C alpha shifts are dominated by dihedral angle changes and can be used to calculate structural change, whereas (13)C beta shifts retain significant dependence on local compression, making them less useful as structural restraints

The interaction of amyloid A beta(1-40) with lipid bilayers and ganglioside as studied by P-31 solid-state NMR

Amyloid P-peptide (A beta) is a major component of plaques in Alzheimer's disease, and formation of senile plaques has been suggested to originate fro m regions of neuronal membrane rich in gangliosides. We analyzed the mode of interaction of A beta with lipid bilayers by multinuclear NMR using P-31 nuclei. We found that A beta (1-40) strongly perturbed the bilayer structure of dimyristoylphosphatidylcholine (DMPQ, to form a non-lamellar phase (most likely micellar). The ganglioside GM1 potentiated the effect of A beta (1-40), as viewed from P-31 NMR. The difference of the isotropic peak intensity between DMPC/A beta and DMPC/GM1/A beta suggests a specific interaction between A beta and GM1. We show that in the DMPC/GM1/A beta system there are three lipid phases, namely a lamellar phase, a hexagonal phase and non-oriented lipids. The latter two phases are induced by the presence of the A beta peptide, and facilitated by GM1. 9) 2008 Elsevier Ireland Ltd. All rights reserved

Drawing single NMR spins and understanding relaxation

How should we draw vectors to represent individual nuclear spins? Vectors in 3D space are merely a way of understanding the mathematics of quantum mechanics, which provides the “true” description of a single spin-½ nucleus. They are a useful aid to understanding, but there is no single “correct” vector representation, and the different vector models that are used have advantages and disadvantages. Here, we discuss the 2 standard vector models for a nuclear magnetic resonance spin: the up/down or alignment model and the 2-cone model, and we show how they relate to quantum mechanics. We show why both of these models are limited and discuss a third model, the uniform model, in which individual spins can be in any orientation. We demonstrate how the uniform model presents a clear and logically coherent description for spins: at equilibrium; following a 90° pulse; and during the subsequent relaxation back to equilibrium. The uniform model is fully consistent with quantum mechanics and leads to an understanding of coherence and relaxation that cannot be obtained from the other 2 models. We suggest that the uniform model is more helpful than the other 2 for most purposes

Bioinformatic analysis of WxL domain proteins

The WxL domain is found on the cell surface of many bacteria, most of which are commensal gut bacteria. Its functions are generally identified as being related to virulence and/or peptidoglycan attachment, but there is so far no clear function or structure for this domain. Here, a range of bioinformatics tools were used to clarify the structure and function. These indicate that WxL domains occur in cell surface-associated gene clusters that always contain a small WxL, large WxL and DUF916 domain; and that the small and large WxL proteins have distinct structure despite sharing two conserved WxL motifs. The two WxL motifs form a hydrophobic surface buried inside the protein. The likely function of the WxL domain is to attach to bacterial peptidoglycan, forming a platform to allow associated domains in the cluster to interact with host proteins

A review on the structure of Bombyx mori silk fibroin fiber studied using solid-state NMR: an antipolar lamella with an 8-residue repeat

Silk fibroin (SF) fiber from the silkworm Bombyx mori in the Silk II form has been used as an excellent textile fiber for over 5000 years. Recently it has been developed for a range of biomedical applications. Further expansion of these uses builds on the excellent mechanical strength of SF fiber, which derives from its structure. This relationship between strength and SF structure has been studied for over 50 years, but it is still not well understood. In this review, we report the use of solid-state NMR to study stable-isotope labeled SF fiber and stable-isotope labeled peptides including (Ala-Gly)15 and (Ala-Gly-Ser-Gly-Ala-Gly)5 as models of the crystalline fraction. We show that the crystalline fraction is a lamellar structure with a repetitive folding using β-turns every eighth amino acid, and that the sidechains adopt an antipolar arrangement rather than the more well-known polar structure described by Marsh, Corey and Pauling (that is, the Ala methyls in each layer point in opposite directions in alternate strands). The amino acids Ser, Tyr and Val are the next most common in B. mori SF after Gly and Ala, and occur in the crystalline and semi-crystalline regions, probably defining the edges of the crystalline region. Thus, we now have an understanding of the main features of Silk II but there is still a long way to go

The accuracy of protein structures in solution determined by AlphaFold and NMR

In the recent Critical Assessment of Structure Prediction (CASP) competition, AlphaFold2 performed outstandingly. Its worst predictions were for nuclear magnetic resonance (NMR) structures, which has two alternative explanations: either the NMR structures were poor, implying that Alpha-Fold may be more accurate than NMR, or there is a genuine difference between crystal and solution structures. Here, we use the program Accuracy of NMR Structures Using RCI and Rigidity (ANSURR), which measures the accuracy of solution structures, and show that one of the NMR structures was indeed poor. We then compare Alpha-Fold predictions to NMR structures and show that Alpha-Fold tends to be more accurate than NMR ensembles. There are, however, some cases where the NMR ensembles are more accurate. These tend to be dynamic structures, where Alpha-Fold had low confidence. We suggest that Alpha-Fold could be used as the model for NMR-structure refinements and that Alpha-Fold structures validated by ANSURR may require no further refinement

Molecular Mechanism for the Hofmeister Effect Derived from NMR and DSC Measurements on Barnase

The effects of sodium thiocyanate, sodium chloride, and sodium sulfate on the ribonuclease barnase were studied using differential scanning calorimetry (DSC) and NMR. Both measurements reveal specific and saturable binding at low anion concentrations (up to 250 mM), which produces localized conformational and energetic effects that are unrelated to the Hofmeister series. The binding of sulfate slows intramolecular motions, as revealed by peak broadening in 13 C heteronuclear single quantum coherence spectroscopy. None of the anions shows significant binding to hydrophobic groups. Above 250 mM, the DSC results are consistent with the expected Hofmeister effects in that the chaotropic anion thiocyanate destabilizes barnase. In this higher concentration range, the anions have approximately linear effects on protein NMR chemical shifts, with no evidence for direct interaction of the anions with the protein surface. We conclude that the effects of the anions on barnase are mediated by solvent interactions. The results are not consistent with the predictions of the preferential interaction, preferential hydration, and excluded volume models commonly used to describe Hofmeister effects. Instead, they suggest that the Hofmeister anion effects on both stability and solubility of barnase are due to the way in which the protein interacts with water molecules, and in particular with water dipoles, which are more ordered around sulfate anions and less ordered around thiocyanate anions

Structural investigation into the threading intercalation of a chiral dinuclear ruthenium(II) polypyridyl complex through a B-DNA oligonucleotide

Herein we report the separation of the three stereoisomers of the DNA light-switch compound [{Ru(bpy)2}2(tpphz)]4+ (tpphz = tetrapyrido[3,2-a:2',3'-c:3″,2″-h:2‴,3‴-j]phenazine) by column chromatography and the characterization of each stereoisomer by X-ray crystallography. The interaction of these compounds with a DNA octanucleotide d(GCATATCG).d(CGATATGC) has been studied using NMR techniques. Selective deuteration of the bipyridyl rings was needed to provide sufficient spectral resolution to characterize structures. NMR-derived structures for these complexes show a threading intercalation binding mode with slow and chirality-dependent rates. This represents the first solution structure of an intercalated bis-ruthenium ligand. Intriguingly, we find that the binding site selectivity is dependent on the nature of the stereoisomer employed, with Λ RuII centers showing a better intercalation fit

DUF916 and DUF3324 in the WxL protein cluster bind to WxL and link bacterial and host surfaces

Bacterial WxL proteins contain peptidoglycan-binding WxL domains, which have a dual Trp-x-Leu motif and are involved in virulence. It was recently shown that WxL proteins occur in gene clusters, containing typically a small WxL protein (which in the mature protein consists only of a WxL domain), a large WxL protein (which contains a C-terminal WxL domain with N-terminal host-binding domains), and a conserved protein annotated as a Domain of Unknown Function (DUF). Here we analyze this DUF and show that it contains two tandem domains—DUF916 and DUF3324—which both have an IgG-like fold and together form a single functional unit, connected to a C-terminal transmembrane helix. DUF3324 is a stable domain, while DUF916 is less stable and is likely to require a stabilizing interaction with WxL. The protein is suggested to have an important role to bind and stabilize WxL on the peptidoglycan surface, via the DUF916 domain, and to bind to host cells via the DUF3324 domain. AlphaFold2 predicts that a β-hairpin strand from DUF916 inserts into WxL adjacent to its N-terminus. We therefore propose to rename the DUF916-DUF3324 pair as WxL Interacting Protein (WxLIP), with DUF916, DUF3324 and the transmembrane helix forming the first, second and third domains of WxLIP, which we characterize as peptidoglycan binding domain (PGBD), host binding domain (HBD), and transmembrane helix (TMH) respectively