Stereospecific assignments of protein NMR resonances based on the tertiary structure and 2D/3D NOE data

In many cases of protein structure determination by NMR a high-quality structure is required. An important contribution to structural precision is stereospecific assignment of magnetically nonequivalent prochiral methylene and methyl groups, eliminating the need for introducing pseudoatoms and pseudoatom corrections in distance restraint lists. Here, we introduce the stereospecific assignment program that uses the resonance assignment, a preliminary 3D structure and 2D and/or 3D nuclear Overhauser effect spectroscopy peak lists for stereospecific assignment. For each prochiral group the algorithm automatically calculates a score for the two different stereospecific assignment possibilities, taking into account the presence and intensity of the nuclear Overhauser effect (NOE) peaks that are expected from the local environment of each prochiral group (i.e., the close neighbors). The performance of the algorithm has been tested and used on NMR data of alpha-helical and beta-sheet proteins using homology models and/or X-ray structures. The program produced no erroneous stereospecific assignments provided the NOEs were carefully picked and the 3D model was sufficiently accurate. The set of NOE distance restraints produced by nmr2st using the results of the SSA module was superior in generating good-quality ensembles of NMR structures (low deviations from upper limits in conjunction with low root-mean-square-deviation values) in the first round of structure calculations. The program uses a novel approach that employs the entire 3D structure of the protein to obtain stereospecific assignment; it can be used to speed up the NMR structure refinement and to increase the quality of the final NMR ensemble even when no scalar or residual dipolar coupling information is available.     

New and old NMR experiments for the resonance assignment of complex oligosaccharides--application to a cyclodextrin derivative

The complete assignment of the (1)H and (13)C sugar resonances in mono-3,6-anhydro-heptakis(2,3-O-methyl)-hexakis(6-O-methyl)-β-cyclodextrin, an asymmetrically functionalized β-cyclodextrin, was carried out by means of 2D NMR experiments. The TOCSY and the homonuclear multiple relay COSY spectra provided most of the (1)H assignments. The multiplicity edited HSQC and a set of F(1) selective HSQC-TOCSY and multiple relay HSQC-COSY spectra gave access to most of the (13)C chemical shifts. The latter were fully and accurately determined by means of a pair of complementary, highly folded HSQC-TOCSY spectra. The TOCSY-ROESY and ROESY-TOCSY spectra yielded the sequential assignment of the sugar units. A high resolution F(1) selective F(1) decoupled version of the TOCSY-ROESY experiment was recorded.     

Influence of NMR Data Completeness on Structure Determinations of Homodimeric Proteins

The structure determination of homodimeric proteins by NMR using conventional NOESY experiments is still challenging due to the degeneracy of the chemical shifts in the identical monomers, which causes ambiguity in the NOE assignments. Residues involved in the interface between two monomers provide essential intermolecular NOEs for the structure determinations of homodimeric proteins. Hence NMR data, such as NOE peak lists and chemical shift assignments of these interface residues, play a crucial role for the successful structure determination of homodimeric proteins. This paper extends our previous report (Lin, Y.‐J.; Kirchner, D. K.; Güntert, P. J. Magn. Reson.­ 2012 , 222, 96) and investigates the influence of incomplete NOESY peak lists combined with incomplete ¹H chemical shift assignments of the interface residues on the structure determination of homodimeric proteins using the program CYANA. Data incompleteness was simulated by random omission of both NOESY cross peaks and interface ¹H chemical shifts. Our results for three proteins with different percentages of interface residues reveal that the algorithm can tolerate about 40–50% NOESY peak omission with complete interface chemical shift assignments, which indicates that partial NOESY peak omission does not cause severe problems when the interface chemical shifts are completely assigned. Combining NOESY peak omission with incomplete interface chemical shift assignments, the tolerance for interface chemical shift omission decreases with the extent of omitted NOESY peaks. The tolerance for unassigned interface side chain, methyl and aromatic chemical shifts is affected more strongly by NOESY peak omission than that for the omission of general interface ¹H chemical shifts including the backbone. In general about 10–30% peaks omission is tolerated in conjunction with missing chemical shift assignments. If more NOESY peaks are omitted calculations gradually become unstable and tend not to tolerate any missing interface chemical shifts. A large amount of omitted NOESY peaks, for instance 30% omission in our calculations, could decrease the tolerance for missing aromatic or methyl interface ¹H chemical shifts to as few as 2–4 missing chemical shifts, suggesting that complete aromatic and methyl ¹H chemical shift assignments are important when the NOESY peak data is significantly incomplete. Finally, for homodimeric proteins with a low percentage of interface residues, our results reveal that the omission of NOESY peaks, even at an extent of only 10%, can result in no tolerance against the omission of interface ¹H chemical shifts, suggesting that the completeness of both interface ¹H chemical shift assignments and NOESY peaks are important for the successful structure determination of proteins with a small homodimer interface.     

NMReDATA,a standard to report the NMR assignment and parameters of organic compounds

Even though NMR has found countless applications in the field of small molecule characterization, there is no standard file format available for the NMR data relevant to structure characterization of small molecules. A new format is therefore introduced to associate the NMR parameters extracted from 1D and 2D spectra of organic compounds to the proposed chemical structure. These NMR parameters, which we shall call NMReDATA (for nuclear magnetic resonance extracted data), include chemical shift values, signal integrals, intensities, multiplicities, scalar coupling constants, lists of 2D correlations, relaxation times, and diffusion rates. The file format is an extension of the existing Structure Data Format, which is compatible with the commonly used MOL format. The association of an NMReDATA file with the raw and spectral data from which it originates constitutes an NMR record. This format is easily readable by humans and computers and provides a simple and efficient way for disseminating results of structural chemistry investigations, allowing automatic verification of published results, and for assisting the constitution of highly needed open‐source structural databases.     

Prediction of the structures of helical membrane proteins based on a minimum unfavorable contacts approach

An understanding of structure–function relationships of membrane proteins continues to be a challenging problem, owing to the difficulty in obtaining their structures experimentally. This study suggests a method for modeling membrane protein structures that can be used to generate a reliable initial conformation prior to the use of other approaches for sampling conformations. It involves optimizing the orientation of hydrophilic residues so as to minimize unfavorable contacts with the hydrophobic tails of the lipid bilayer. Starting with the optimized initial conformation for three different proteins modeled based on this method, two independent approaches have been used for sampling the conformational space of the proteins. Both approaches are able to predict structures reasonably close to experimental structures, indicating that the initial structure enables the sampling of conformations that are close to the native structure. Possible improvements in the method for making it broadly applicable to helical membrane proteins are discussed. © 2015 Wiley Periodicals, Inc.     

hNCOcanH pulse sequence and a robust protocol for rapid and unambiguous assignment of backbone (1HN, 15N and 13C′) resonances in 15 N/13C‐labeled proteins

A three‐dimensional nuclear magnetic resonance (NMR) pulse sequence named as hNCOcanH has been described to aid rapid sequential assignment of backbone resonances in ¹⁵N/¹³C‐labeled proteins. The experiment has been derived by a simple modification of the previously described HN(C)N pulse sequence [Panchal et al., J. Biomol. NMR 20 (2001) 135–147]; t₂ evolution is used to frequency label ¹³C′ rather than ¹⁵N (similar trick has also been used in the design of hNCAnH pulse sequence from hNcaNH [Frueh et al., JACS, 131 (2009) 12880–12881]). The modification results in a spectrum equivalent to HNCO, but in addition to inter‐residue correlation peaks (i.e. H_i, C_i?1), the spectrum also contains additional intra‐residue correlation peaks (i.e. H_i?1, C_i?1) in the direct proton dimension which has maximum resolution. This is the main strength of the experiment and thus, even a small difference in amide ¹H chemical shifts (5–6 Hz) can be used for establishing a sequential connectivity. This experiment in combination with the HNN experiment described previously [Panchal et al., J. Biomol. NMR 20 (2001) 135–147] leads to a more robust assignment protocol for backbone resonances (¹H^N, ¹⁵N) than could be derived from the combination of HNN and HN(C)N experiments [Bhavesh et al., Biochemistry, 40 (2001) 14727–14735]. Further, this new protocol enables assignment of ¹³C′ resonances as well. We believe that the experiment and the protocol presented here will be of immense value for structural—and functional—proteomics research by NMR. Performance of this experiment has been demonstrated using ¹³C/¹⁵N labeled ubiquitin. Copyright © 2011 John Wiley & Sons, Ltd.     

Complete assignment of the NMR spectra of [D‐Leu1]‐microcystin‐LR and analysis of its solution structure

[D ‐Leu¹]‐microcystin‐LR is a recently discovered microcystin. We report the isolation of this microcystin analogue from a Microcystis aeruginosa strain isolated from the Lagoa de Iquipari, Rio de Janeiro, Brazil. The ¹H and ¹³C NMR spectra were completely assigned in both MeOH‐d₄ and DMSO‐d₆. Further, the solution structure of this compound was investigated with the use of two‐dimensional NMR and the amide proton temperature dependence, and was compared with those of its analogs, microcystin‐RR and microcystin‐LR. Copyright © 2002 John Wiley & Sons, Ltd.     

Quindolinocryptotackieine: the elucidation of a novel indoloquinoline alkaloid structure through the use of computer‐assisted structure elucidation and 2D NMR

Numerous indoloquinoline alkaloid structures have been identified from extracts of the West African plant Cryptolepis sanguinolenta. Recently, through the use of 2D NMR methods and cryogenic NMR probe technology in conjunction with computer‐assisted structure elucidation (CASE) methods, the structures of some chemical degradation products of this family of alkaloids have also been reported. We now report the characterization of a novel indoloquinoline dimeric alkaloid, quindolinocryptotackieine, through the extensive utilization of CASE methods. The NMR data presented here were collected over a decade earlier before the elucidation of the structure was possible, since manual analysis did not present a conclusive structure, whereas CASE produced a series of structures from which the structure could be verified. The original mass spectrometric (MS) data collected for the sample were problematic. Contemporary MS data were instead recollected from remaining small quantities of this alkaloid using modern instrumentation. The re‐collected data gave a usable molecular ion and several key fragment ions that were diagnostically useful. Copyright © 2003 John Wiley & Sons, Ltd.     

Automated structure verification based on a combination of 1D (1)H NMR and 2D (1)H - (13)C HSQC spectra

A method for structure validation based on the simultaneous analysis of a 1D (1)H NMR and 2D (1)H - (13)C single-bond correlation spectrum such as HSQC or HMQC is presented here. When compared with the validation of a structure by a 1D (1)H NMR spectrum alone, the advantage of including a 2D HSQC spectrum in structure validation is that it adds not only the information of (13)C shifts, but also which proton shifts they are directly coupled to, and an indication of which methylene protons are diastereotopic. The lack of corresponding peaks in the 2D spectrum that appear in the 1D (1)H spectrum, also gives a clear picture of which protons are attached to heteroatoms. For all these benefits, combined NMR verification was expected and found by all metrics to be superior to validation by 1D (1)H NMR alone. Using multiple real-life data sets of chemical structures and the corresponding 1D and 2D data, it was possible to unambiguously identify at least 90% of the correct structures. As part of this test, challenging incorrect structures, mostly regioisomers, were also matched with each spectrum set. For these incorrect structures, the false positive rate was observed as low as 6%.     

Hybrid Structure of the Type 1 Pilus of Uropathogenic Escherichia coli

Type 1 pili are filamentous protein assemblies on the surface of Gram‐negative bacteria that mediate adhesion to host cells during the infection process. The molecular structure of type 1 pili remains elusive on the atomic scale owing to their insolubility and noncrystallinity. Herein we describe an approach for hybrid‐structure determination that is based on data from solution‐state NMR spectroscopy on the soluble subunit and solid‐state NMR spectroscopy and STEM data on the assembled pilus. Our approach is based on iterative modeling driven by structural information extracted from different sources and provides a general tool to access pseudo atomic structures of protein assemblies with complex subunit folds. By using this methodology, we determined the local conformation of the FimA pilus subunit in the context of the assembled type 1 pilus, determined the exact helical pilus architecture, and elucidated the intermolecular interfaces contributing to pilus assembly and stability with atomic detail.     

Scaling law for the radius of gyration of proteins and its dependence on hydrophobicity

The scaling law between the radius of gyration and the length of a polymer chain has long been an interesting topic since the Flory theory. In this article, we seek to derive a unified formula for the scaling exponent of proteins under different solvent conditions. The formula is obtained by considering the balance between the excluded volume effect and elastic interactions among monomers. Our results show that the scaling exponent is closely related to the fractal dimension of a protein's structure at the equilibrium state. Applying this formula to natural proteins yields a 2/5 law with fractal dimension 2 at the native state, which is in good agreement with other studies based on Protein Data Bank analysis. We also study the dependence of the scaling exponent on the hydrophobicity of a protein chain through a simple two‐letters HP model. The results provides a way to estimate the globular structure of a protein, and could be helpful for the investigation of the mechanisms of protein folding. © 2008 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 47: 207–214, 2009     

Effect of Site-Specific O-Glycosylation on the Structural Behavior of NOTCH1 Receptor Extracellular EGF-like Domains 11 and 10

Human NOTCH1 receptor contains 36 epidermal growth factor (EGF)-like repeating domains, in which O-glycosylation status of EGF12 domain regulates the interaction with Notch ligands. Our interest is focused on the effect of specific O-glycosylation states on the structural behavior of EGF11 and EGF10, because they appeared to affect molecular mechanism in receptor–ligand interactions by inducing some conformational alterations in these domains and/or the regions connecting two domains. To understand the structural impact of various O-glycosylation patterns on the pivotal EGF-like repeats 10, 11, and 12, we performed chemical synthesis and NMR studies of site-specifically O-glycosylated EGF11 and EGF10. Our strategy enabled us to synthesize four EGF11 and five EGF10 modules. The specific O-glycosylation states affected in vitro folding of EGF10 more than EGF11, while calcium ion had a larger effect on EGF11 folding. Comprehensive NMR studies shed light on the new type “sugar bridges” crosslinking Thr-O-GlcNAc in the consensus sequence C5-X-X-G-X-(T/S)-G-X-X-C6 and an amino acid in the hinge region between the domains, 445Thr-O-GlcNAc—IIe451 in domain 11 and 405Thr-O-GlcNAc—Gln411 in domain 10, respectively.     

Stereospecificity of 1H, 13C and 15N shielding constants in the isomers of methylglyoxal bisdimethylhydrazone: problem with configurational assignment based on 1H chemical shifts

In the ¹³C NMR spectra of methylglyoxal bisdimethylhydrazone, the ¹³C‐5 signal is shifted to higher frequencies, while the ¹³C‐6 signal is shifted to lower frequencies on going from the EE to ZE isomer following the trend found previously. Surprisingly, the ¹H‐6 chemical shift and ¹J(C‐6,H‐6) coupling constant are noticeably larger in the ZE isomer than in the EE isomer, although the configuration around the –CH═N– bond does not change. This paradox can be rationalized by the C–H?N intramolecular hydrogen bond in the ZE isomer, which is found from the quantum‐chemical calculations including Bader's quantum theory of atoms in molecules analysis. This hydrogen bond results in the increase of δ(¹H‐6) and ¹J(C‐6,H‐6) parameters. The effect of the C–H?N hydrogen bond on the ¹H shielding and one‐bond ¹³C–¹H coupling complicates the configurational assignment of the considered compound because of these spectral parameters. The ¹H, ¹³C and ¹⁵N chemical shifts of the 2‐ and 8‐(CH₃)₂N groups attached to the –C(CH₃)═N– and –CH═N– moieties, respectively, reveal pronounced difference. The ab initio calculations show that the 8‐(CH₃)₂N group conjugate effectively with the π‐framework, and the 2‐(CH₃)₂N group twisted out from the plane of the backbone and loses conjugation. As a result, the degree of charge transfer from the N‐2– and N‐8– nitrogen lone pairs to the π‐framework varies, which affects the ¹H, ¹³C and ¹⁵N shieldings. Copyright © 2012 John Wiley & Sons, Ltd.     

Identification of Potential HCV Inhibitors Based on the Interaction of Epigallocatechin-3-Gallate with Viral Envelope Proteins

Hepatitis C is affecting millions of people around the globe annually, which leads to death in very high numbers. After many years of research, hepatitis C virus (HCV) remains a serious threat to the human population and needs proper management. The in silico approach in the drug discovery process is an efficient method in identifying inhibitors for various diseases. In our study, the interaction between Epigallocatechin-3-gallate, a component of green tea, and envelope glycoprotein E2 of HCV is evaluated. Epigallocatechin-3-gallate is the most promising polyphenol approved through cell culture analysis that can inhibit the entry of HCV. Therefore, various in silico techniques have been employed to find out other potential inhibitors that can behave as EGCG. Thus, the homology modelling of E2 protein was performed. The potential lead molecules were predicted using ligand-based as well as structure-based virtual screening methods. The compounds obtained were then screened through PyRx. The drugs obtained were ranked based on their binding affinities. Furthermore, the docking of the topmost drugs was performed by AutoDock Vina, while its 2D interactions were plotted in LigPlot+. The lead compound mms02387687 (2-[[5-[(4-ethylphenoxy) methyl]-4-prop-2-enyl-1,2,4-triazol-3-yl] sulfanyl]-N-[3(trifluoromethyl) phenyl] acetamide) was ranked on top, and we believe it can serve as a drug against HCV in the future, owing to experimental validation.     

Effect of Phosphorylation on the Structural Behaviour of Peptides Derived from the Intrinsically Disordered C-Terminal Domain of Histone H1.0

To investigate the structural impact of phosphorylation on the human histone H1.0 C-terminal domain, we performed NMR structural studies of model peptides containing a single phosphorylation site: T¹¹⁸-H1.0 (T¹¹⁸PKK motif) and T¹⁴⁰-H1.0 (T¹⁴⁰PVK motif). Both model peptides are mainly disordered in aqueous solution in their non-phosphorylated and phosphorylated forms, but become structured in the presence of trifluoroethanol. The peptides T¹¹⁸-H1.0 and pT¹¹⁸-H1.0 contain two helical regions, a long amphipathic α helix spanning residues 104–115 and a short α/3₁₀ helix (residues 119–123), that are almost perpendicular in T¹¹⁸-H1.0 but have a poorly defined orientation in pT¹¹⁸-H1.0. Peptides T¹⁴⁰-H1.0 and pT¹⁴⁰-H1.0 form very similar α helices between residues 141–147. The TPKK and TPVK motifs show the same backbone conformation, but differ in their side-chain contacts; the Thr and pThr side chains interact with the i+2 Lys side chain in the TPKK motif, and with the i+3 Lys side chain in the TPVK motif. The pT phosphate group in pT¹¹⁸-H1.0 and pT¹⁴⁰-H1.0 has pK_a values below the intrinsic values, which can be explained by non-specific charge–charge interactions with nearby Lys. The non-polar Val in the TPVK motif accounts for the pT¹⁴⁰ pK_a being closer to the intrinsic pK_a value than the pT¹¹⁸ pK_a. Altogether, these results validate that minimalist strategies using model peptides can provide structural details difficult to obtain in short-lived intrinsically disordered proteins and domains.     

Identification of β-barrel membrane proteins based on amino acid composition properties and predicted secondary structure

Unlike all-helices membrane proteins, β-barrel membrane proteins can not be successfully discriminated from other proteins, especially from all-β soluble proteins. This paper performs an analysis on the amino acid composition in membrane parts of 12 β-barrel membrane proteins versus β-strands of 79 all-β soluble proteins. The average and variance of the amino acid composition in these two classes are calculated. Amino acids such as Gly, Asn, Val that are most likely associated with classification are selected based on Fishers discriminant ratio. A linear classifier built with these selected amino acids composition in observed β-strands achieves 100% classification accuracy for 12 membrane proteins and 79 soluble proteins in a four-fold cross-validation experiment. Since at present the accuracy of secondary structure prediction is quite high, a promising method to identify β-barrel membrane proteins is presented based on the linear classifier coupled with predicted secondary structure. Applied to 241 β-barrel membrane proteins and 3855 soluble proteins with various structures, the method achieves 85.48% (206/241) sensitivity and 92.53% specificity (3567/3855).