Use of chemical shifts in macromolecular structure determination

Conformational Dynamics of an Amyloidogenic Intermediate of Transthyretin: Implications for Structural Remodeling and Amyloid Formation

The aggregation pathway of transthyretin (TTR) proceeds through rate-limiting dissociation of the tetramer (a dimer of dimers) and partial misfolding of the resulting monomer, which assembles into amyloid structures through a downhill polymerization mechanism. The structural features of the aggregation-prone monomeric intermediate are poorly understood. NMR relaxation dispersion offers a unique opportunity to characterize amyloidogenic intermediates when they exchange on favorable timescales with NMR-visible ground states. Here we use NMR to characterize the structure and conformational dynamics of the monomeric F87E mutant of human TTR. Chemical shifts derived from analysis of multinuclear relaxation dispersion data provide insights into the structure of a low-lying excited state that exchanges with the ground state of the F87E monomer at a rate of 3800 s-1. Disruption of the subunit interfaces of the TTR tetramer leads to destabilization of edge strands in both β-sheets of the F87E monomer. Conformational fluctuations are propagated through the entire hydrogen bonding network of the DAGH β-sheet, from the inner β-strand H, which forms the strong dimer-dimer interface in the TTR tetramer, to outer strand D which is unfolded in TTR fibrils. Fluctuations are also propagated from the AB loop in the weak dimer-dimer interface to the EF helix, which undergoes structural remodeling in fibrils. The conformational fluctuations in both regions are enhanced at acidic pH where amyloid formation is most favorable. The relaxation dispersion data provide insights into the conformational dynamics of the amyloidogenic state of monomeric TTR that predispose it for structural remodeling and progression to amyloid fibrils.

NMR-Based Characterization of the Interaction between Yeast Oxa1-CTD and Ribosomes

In mitochondria, the major subunits of oxidative phosphorylation complexes are translated by the mitochondrial ribosome (mito-ribosome). The correct insertion and assembly of these subunits into the inner mitochondrial membrane (IMM) are facilitated by mitochondrial oxidase assembly protein 1 (Oxa1) during the translation process. This co-translational insertion process involves an association between the mito-ribosome and the C-terminus of Oxa1 (Oxa1-CTD) Nuclear magnetic resonance (NMR) methods were mainly used to investigate the structural characterization of yeast Oxa1-CTD and its mode of interaction with the E. coli 70S ribosome. Oxa1-CTD forms a transient α-helical structure within the residues P342-Q385, which were reported to form an α-helix when combining with the ribosome. Two conserved contact sites that could interact with the ribosome were further identified. The first site was located on the very end of the N-terminus (V321-I327), and the second one encompassed a stretch of amino acid residues I348-Q370. Based on our discoveries and previous reports, a model has been proposed in which Oxa1-CTD interacts with ribosomes, accompanied by transient-to-stable transitions at the second contact site. These observations may enhance our understanding of the potential role of Oxa1-CTD in facilitating the assembly of oxidative phosphorylation complexes and provide insight into the structural characteristics of Oxa1-CTD.

NMR-Chemical-Shift-Driven Protocol Reveals the Cofactor-Bound, Complete Structure of Dynamic Intermediates of the Catalytic Cycle of Oncogenic KRAS G12C Protein and the Significance of the Mg2+ Ion

In this work, catalytically significant states of the oncogenic G12C variant of KRAS, those of Mg2+-free and Mg2+-bound GDP-loaded forms, have been determined using CS-Rosetta software and NMR-data-driven molecular dynamics simulations. There are several Mg2+-bound G12C KRAS/GDP structures deposited in the Protein Data Bank (PDB), so this system was used as a reference, while the structure of the Mg2+-free but GDP-bound state of the RAS cycle has not been determined previously. Due to the high flexibility of the Switch-I and Switch-II regions, which also happen to be the catalytically most significant segments, only chemical shift information could be collected for the most important regions of both systems. CS-Rosetta was used to derive an "NMR ensemble" based on the measured chemical shifts, which, however, did not contain the nonprotein components of the complex. We developed a torsional restraint set for backbone torsions based on the CS-Rosetta ensembles for MD simulations, overriding the force-field-based parametrization in the presence of the reinserted cofactors. This protocol (csdMD) resulted in complete models for both systems that also retained the structural features and heterogeneity defined by the measured chemical shifts and allowed a detailed comparison of the Mg2+-bound and Mg2+-free states of G12C KRAS/GDP.

Atomic Resolution Insights into pH Shift Induced Deprotonation Events in LS-Shaped Aβ(1-42) Amyloid Fibrils

Alzheimer's disease is a neurodegenerative disorder associated with the deposition of misfolded aggregates of the amyloid-β protein (Aβ). Aβ(1-42) is one of the most aggregation-prone components in senile plaques of AD patients. We demonstrated that relatively homogeneous Aβ(1-42) fibrils with one predominant fold visible in solid-state NMR spectra can be obtained at acidic pH. The structure of these fibrils differs remarkably from some other polymorphs obtained at neutral pH. In particular, the entire N-terminal region is part of the rigid fibril core. Here, we investigate the effects of a pH shift on the stability and the fold of these fibrils at higher pH values. Fibril bundling at neutral pH values renders cryo-EM studies impractical, but solid-state NMR spectroscopy, molecular dynamics simulations, and biophysical methods provide residue-specific structural information under these conditions. The LS-fold of the Aβ(1-42) fibrils does not change over the complete pH range from pH 2 to pH 7; in particular, the N-terminus remains part of the fibril core. We observe changes in the protonation state of charged residues starting from pH 5 on a residue-specific level. The deprotonation of the C-terminal carboxyl group of A42 in the intermolecular salt bridge with D1 and K28 is slow on the NMR time scale, with a local pK_a of 5.4, and local conformations of the involved residues are affected by deprotonation of A42. Thus, we demonstrate that this fibril form is stable at physiological pH values.

Structural and dynamical determinants of a β-sheet-enriched intermediate involved in amyloid fibrillar assembly of human prion protein

The conformational conversion of the cellular prion protein (PrP^C) into a misfolded, aggregated and infectious scrapie isoform is associated with prion disease pathology and neurodegeneration. Despite the significant number of experimental and theoretical studies the molecular mechanism regulating this structural transition is still poorly understood. Here, via Nuclear Magnetic Resonance (NMR) methodologies we investigate at the atomic level the mechanism of the human HuPrP(90–231) thermal unfolding and characterize the conformational equilibrium between its native structure and a β-enriched intermediate state, named β-PrPI. By comparing the folding mechanisms of metal-free and Cu²⁺-bound HuPrP(23–231) and HuPrP(90–231) we show that the coupling between the N- and C-terminal domains, through transient electrostatic interactions, is the key molecular process in tuning long-range correlated μs–ms dynamics that in turn modulate the folding process. Moreover, via thioflavin T (ThT)-fluorescence fibrillization assays we show that β-PrPI is involved in the initial stages of PrP fibrillation, overall providing a clear molecular description of the initial phases of prion misfolding. Finally, we show by using Real-Time Quaking-Induced Conversion (RT-QuIC) that the β-PrPI acts as a seed for the formation of amyloid aggregates with a seeding activity comparable to that of human infectious prions.

The N-ter domain in HuPrP regulates the folding mechanism by tuning the long-range μs–ms dynamics. Removal of the N-ter domain triggers the formation of a stable β-enriched intermediate state inducing amyloid aggregates with HuPrP^Sc seeding activity.

The auto-inhibition mechanism of transcription factor Ets-1 induced by phosphorylation on the intrinsically disordered region

As the most abundant post-translation modifications (PTMs), the phosphorylation usually occurred on the intrinsically disordered regions (IDRs). The regulation on the structures and interactions of IDRs induced by phosphorylation is critical to the function performing. The eukaryotic transcription factor 1 (Ets-1) is a member of transcription factor family, which participates in many important biological processes. The DNA-binding ability of Ets-1 is auto-inhibited by a disordered serine-rich region (SRR) on the Ets-1. The inhibition ability of SRR is greatly enhanced by the phosphorylation of the serine on the SRR. Nevertheless, the molecular mechanisms of the phosphorylation regulation on the structure and activity of Ets-1 are still unclear and under debates. By using both of the molecular simulations and biochemical experiments, we studied the molecule mechanism of phosphorylation regulation on the auto-inhibition of the Ets-1. The reasons of stabilization of Ets-1 core by phosphorylation on SRR region were elucidated. More important, the free energy landscapes (FEL) show that both of the steric hindrance and allosteric regulation are responsible for the DNA-binding inhibitory induced by phosphorylation, but the steric effects contribute greater than the allosteric regulation. The phosphorylation not only enhances the electrostatic interactions to facilitate the steric impedance, but also promotes the formation of hydrophobic residue clusters, which provide major driven force for the allosteric regulation. The structural basis of auto-inhibition of Ets-1 induced by the phosphorylation revealed in this study would great help the developing of inhibitor for the cancer therapy.

An Improved Turn Structure for Inducing β-Hairpin Formation in Peptides

Although β-hairpins are widespread in proteins, there is no tool to coax any small peptide to adopt a β-hairpin conformation, regardless of sequence. Here, we report that δ-linked γ(R)-methyl-ornithine (^δ MeOrn) provides an improved β-turn template for inducing a β-hairpin conformation in peptides. We developed a synthesis of protected ^δ MeOrn as a building block suitable for use in Fmoc-based solid-phase peptide synthesis. The synthesis begins with l-leucine and affords gram quantities of the N^α -Boc-N^δ -Fmoc-γ(R)-methyl-ornithine building block. X-ray crystallography confirms that the ^δ MeOrn turn unit adopts a folded structure in a macrocyclic β-hairpin peptide. CD and NMR spectroscopy allow comparison of the ^δ MeOrn turn template to the δ-linked ornithine (^δ Orn) turn template that we previously introduced and to the popular d-Pro-Gly turn template. These studies show that the folding of the ^δ MeOrn turn template is substantially better than that of ^δ Orn and is comparable to d-Pro-Gly.

Design, Optimization, and Structural Characterization of an Apoptosis-Inducing Factor Peptide Targeting Human Cyclophilin A to Inhibit Apoptosis Inducing Factor-Mediated Cell Death

Blocking the interaction between the apoptosis-inducing factor (AIF) and cyclophilin A (CypA) by the AIF fragment AIF(370-394) is protective against glutamate-induced neuronal cell death and brain injury in mice. Starting from AIF(370-394), we report the generation of the disulfide-bridged and shorter variant AIF(381-389) and its structural characterization by nuclear magnetic resonance (NMR) in the free and CypA-bound state. AIF(381-389) in both the free and bound states assumes a β-hairpin conformation similar to that of the fragment in the AIF protein and shows a highly reduced conformational flexibility. This peptide displays a similar in vitro affinity for CypA, an improved antiapoptotic activity in cells and an enhanced proteolytic stability compared to the parent peptide. The NMR-based 3D model of the AIF(381-389)/CypA complex provides a better understanding of the binding hot spots on both the peptide and the protein and can be exploited to design AIF/CypA inhibitors with improved pharmacokinetic and pharmacodynamics features.

Informing NMR experiments with molecular dynamics simulations to characterize the dominant activated state of the KcsA ion channel

As the first potassium channel with an x-ray structure determined, and given its homology to eukaryotic channels, the pH-gated prokaryotic channel KcsA has been extensively studied. Nevertheless, questions related, in particular, to the allosteric coupling between its gates remain open. The many currently available x-ray crystallography structures appear to correspond to various stages of activation and inactivation, offering insights into the molecular basis of these mechanisms. Since these studies have required mutations, complexation with antibodies, and substitution of detergents in place of lipids, examining the channel under more native conditions is desirable. Solid-state nuclear magnetic resonance (SSNMR) can be used to study the wild-type protein under activating conditions (low pH), at room temperature, and in bacteriomimetic liposomes. In this work, we sought to structurally assign the activated state present in SSNMR experiments. We used a combination of molecular dynamics (MD) simulations, chemical shift prediction algorithms, and Bayesian inference techniques to determine which of the most plausible x-ray structures resolved to date best represents the activated state captured in SSNMR. We first identified specific nuclei with simulated NMR chemical shifts that differed significantly when comparing partially open vs fully open ensembles from MD simulations. The simulated NMR chemical shifts for those specific nuclei were then compared to experimental ones, revealing that the simulation of the partially open state was in good agreement with the SSNMR data. Nuclei that discriminate effectively between partially and fully open states belong to residues spread over the sequence and provide a molecular level description of the conformational change.

Thermodynamics and kinetics of the amyloid-β peptide revealed by Markov state models based on MD data in agreement with experiment

The amlyoid-β peptide (Aβ) is closely linked to the development of Alzheimer's disease. Molecular dynamics (MD) simulations have become an indispensable tool for studying the behavior of this peptide at the atomistic level. General key aspects of MD simulations are the force field used for modeling the peptide and its environment, which is important for accurate modeling of the system of interest, and the length of the simulations, which determines whether or not equilibrium is reached. In this study we address these points by analyzing 30-μs MD simulations acquired for Aβ40 using seven different force fields. We assess the convergence of these simulations based on the convergence of various structural properties and of NMR and fluorescence spectroscopic observables. Moreover, we calculate Markov state models for the different MD simulations, which provide an unprecedented view of the thermodynamics and kinetics of the amyloid-β peptide. This further allows us to provide answers for pertinent questions, like: which force fields are suitable for modeling Aβ? (a99SB-UCB and a99SB-ILDN/TIP4P-D); what does Aβ peptide really look like? (mostly extended and disordered) and; how long does it take MD simulations of Aβ to attain equilibrium? (at least 20-30 μs). We believe the analyses presented in this study will provide a useful reference guide for important questions relating to the structure and dynamics of Aβ in particular, and by extension other similar disordered proteins.

Multiple copies of a novel amphipathic α-helix forming segment in Physcomitrella patens dehydrin play a key role in abiotic stress mitigation

Plants use a diverse set of proteins to mitigate various abiotic stresses. The intrinsically disordered protein dehydrin is an important member of this repertoire of proteins, characterized by a canonical amphipathic K-segment. It can also contain other stress-mitigating noncanonical segments—a likely reflection of the extremely diverse nature of abiotic stress encountered by plants. Among plants, the poikilohydric mosses have no inbuilt mechanism to prevent desiccation and therefore are likely to contain unique noncanonical stress-responsive motifs in their dehydrins. Here we report the recurring occurrence of a novel amphipathic helix-forming segment (D-segment: EGφφD(R/K)AKDAφ, where φ represents a hydrophobic residue) in Physcomitrella patens dehydrin (PpDHNA), a poikilohydric moss. NMR and CD spectroscopic experiments demonstrated the helix-forming tendency of the D-segment, with the shuffled D-segment as control. PpDHNA activity was shown to be size as well as D-segment dependent from in vitro, in vivo, and in planta studies using PpDHNA and various deletion mutants. Bimolecular fluorescence complementation studies showed that D-segment-mediated PpDHNA self-association is a requirement for stress abatement. The D-segment was also found to occur in two rehydrin proteins from Syntrichia ruralis, another poikilohydric plant like P. patens. Multiple occurrences of the D-segment in poikilohydric plant dehydrins/rehydrins, along with the experimental demonstration of the role of D-segment in stress abatement, implies that the D-segment mediates unique resurrection strategies, which may be employed by plant dehydrins that are capable of mitigating extreme stress.

Chemical shift prediction of RNA imino groups: application toward characterizing RNA excited states

NH groups in proteins or nucleic acids are the most challenging target for chemical shift prediction. Here we show that the RNA base pair triplet motif dictates imino chemical shifts in its central base pair. A lookup table is established that links each type of base pair triplet to experimental chemical shifts of the central base pair, and can be used to predict imino chemical shifts of RNAs to remarkable accuracy. Strikingly, the semiempirical method can well interpret the variations of chemical shifts for different base pair triplets, and is even applicable to non-canonical motifs. This finding opens an avenue for predicting chemical shifts of more complicated RNA motifs. Furthermore, we combine the imino chemical shift prediction with NMR relaxation dispersion experiments targeting both ¹⁵N and ¹H^N of the imino group, and verify a previously characterized excited state of P5abc subdomain including an earlier speculated non-native G•G mismatch.

Prediction of chemical shifts is critical for extracting structural and dynamic information from biomolecular NMR data. Here the authors report an RNA imino group chemical shift predictor, showing that the imino chemical shifts of a residue are dictated by the surrounding base pair triplet.

The accuracy of force fields on the simulation of intrinsically disordered proteins: A benchmark test on the human p53 tumor suppressor

Intrinsically disordered proteins (IDPs) are a class of proteins without stable three-dimensional structures under physiological conditions. IDPs exhibit high dynamic nature and could be described by structural ensembles. As one of the most widely used tools, molecular dynamics (MD) simulation could provide the atomic descriptions of the structural ensemble of IDPs. However, the accuracy of the MD simulation largely depends on the accuracy of the force field. In this paper, we compared the structural ensembles of the activation domain 1 (AD1) in p53 tumor suppressor obtained from the widely used force fields, AMBER99SB-ILDN, CHARMM27, CHARMM36m with different water models. The results show that CHARMM36m generates more extended conformations than other force fields, while CHARMM27 prefers to sample the [Formula: see text]-helical structure. Moreover, the chemical shifts obtained by CHARMM36m are the closest to the experimental measurements. These results indicate that the CHARMM36m force field performs best in characterizing the structure properties of p53 AD1. Water models are also critical to describe the structural ensemble of IDPs. TIP4P water model can obtain more extended conformations and produce more local helical conformations than the TIP3P model in our simulation. In addition, we also compare the chemical shifts predicted by different chemical shift predicting programs with experimental measurements, the results show that SHIFTX2 obtains the best performance in the chemical shifts prediction.

Recent Development of Computational Predicting Bioluminescent Proteins

Bioluminescent Proteins (BLPs) are widely distributed in many living organisms that act as a key role of light emission in bioluminescence. Bioluminescence serves various functions in finding food and protecting the organisms from predators. With the routine biotechnological application of bioluminescence, it is recognized to be essential for many medical, commercial and other general technological advances. Therefore, the prediction and characterization of BLPs are significant and can help to explore more secrets about bioluminescence and promote the development of application of bioluminescence. Since the experimental methods are money and time-consuming for BLPs identification, bioinformatics tools have played important role in fast and accurate prediction of BLPs by combining their sequences information with machine learning methods. In this review, we summarized and compared the application of machine learning methods in the prediction of BLPs from different aspects. We wish that this review will provide insights and inspirations for researches on BLPs.

Impact of the Hereditary P301L Mutation on the Correlated Conformational Dynamics of Human Tau Protein Revealed by the Paramagnetic Relaxation Enhancement NMR Experiments

Tau forms intracellular insoluble aggregates as a neuropathological hallmark of Alzheimer’s disease. Tau is largely unstructured, which complicates the characterization of the tau aggregation process. Recent studies have demonstrated that tau samples two distinct conformational ensembles, each of which contains the soluble and aggregation-prone states of tau. A shift to populate the aggregation-prone ensemble may promote tau fibrillization. However, the mechanism of this ensemble transition remains elusive. In this study, we explored the conformational dynamics of a tau fragment by using paramagnetic relaxation enhancement (PRE) and interference (PRI) NMR experiments. The PRE correlation map showed that tau is composed of segments consisting of residues in correlated motions. Intriguingly, residues forming the β-structures in the heparin-induced tau filament coincide with residues in these segments, suggesting that each segment behaves as a structural unit in fibrillization. PRI data demonstrated that the P301L mutation exclusively alters the transiently formed tau structures by changing the short- and long-range correlated motions among residues. The transient conformations of P301L tau expose the amyloid motif PHF6 to promote tau self-aggregation. We propose the correlated motions among residues within tau determine the population sizes of the conformational ensembles, and perturbing the correlated motions populates the aggregation-prone form.

The dynein light chain 8 (LC8) binds predominantly "in-register" to a multivalent intrinsically disordered partner

Dynein light chain 8 (LC8) interacts with intrinsically disordered proteins (IDPs) and influences a wide range of biological processes. It is becoming apparent that among the numerous IDPs that interact with LC8, many contain multiple LC8-binding sites. Although it is established that LC8 forms parallel IDP duplexes with some partners, such as nucleoporin Nup159 and dynein intermediate chain, the molecular details of these interactions and LC8's interactions with other diverse partners remain largely uncharacterized. LC8 dimers could bind in either a paired "in-register" or a heterogeneous off-register manner to any of the available sites on a multivalent partner. Here, using NMR chemical shift perturbation, analytical ultracentrifugation, and native electrospray ionization MS, we show that LC8 forms a predominantly in-register complex when bound to an IDP domain of the multivalent regulatory protein ASCIZ. Using saturation transfer difference NMR, we demonstrate that at substoichiometric LC8 concentrations, the IDP domain preferentially binds to one of the three LC8 recognition motifs. Further, the differential dynamic behavior for the three sites and the size of the fully bound complex confirmed an in-register complex. Dynamics measurements also revealed that coupling between sites depends on the linker length separating these sites. These results identify linker length and motif specificity as drivers of in-register binding in the multivalent LC8-IDP complex assembly and the degree of compositional and conformational heterogeneity as a promising emerging mechanism for tuning of binding and regulation.

Intrinsic disorder and amino acid specificity modulate binding of the WW2 domain in kidney and brain protein (KIBRA) to synaptopodin

The second WW domain (WW2) of the kidney and brain scaffolding protein, KIBRA, has an isoleucine (Ile-81) rather than a second conserved tryptophan and is primarily unstructured. However, it adopts the canonical triple-stranded antiparallel β-sheet structure of WW domains when bound to a two-PPXY motif peptide of the synaptic protein Dendrin. Here, using a series of biophysical experiments, we demonstrate that the WW2 domain remains largely disordered when bound to a 69-residue two-PPXY motif polypeptide of the synaptic and podocyte protein synaptopodin (SYNPO). Isothermal titration calorimetry and CD experiments revealed that the interactions of the disordered WW2 domain with SYNPO are significantly weaker than SYNPO's interactions with the well-folded WW1 domain and that an I81W substitution in the WW2 domain neither enhances binding affinity nor induces substantial WW2 domain folding. In the tandem polypeptide, the two WW domains synergized, enhancing the overall binding affinity with the I81W variant tandem polypeptide 2-fold compared with the WT polypeptide. Solution NMR results showed that SYNPO binding induces small but definite chemical shift perturbations in the WW2 domain, confirming the disordered state of the WW2 domain in this complex. These analyses also disclosed that SYNPO binds the tandem WW domain polypeptide in an antiparallel manner, that is, the WW1 domain binds the second PPXY motif of SYNPO. We propose a binding model consisting of a bipartite interaction mode in which the largely disordered WW2 forms a "fuzzy" complex with SYNPO. This binding mode may be important for specific cellular functions.

Exploring the role of post-translational modifications in regulating α-synuclein interactions by studying the effects of phosphorylation on nanobody binding

Intracellular deposits of α-synuclein in the form of Lewy bodies are major hallmarks of Parkinson's disease (PD) and a range of related neurodegenerative disorders. Post-translational modifications (PTMs) of α-synuclein are increasingly thought to be major modulators of its structure, function, degradation and toxicity. Among these PTMs, phosphorylation near the C-terminus at S129 has emerged as a dominant pathogenic modification as it is consistently observed to occur within the brain and cerebrospinal fluid (CSF) of post-mortem PD patients, and its level appears to correlate with disease progression. Phosphorylation at the neighboring tyrosine residue Y125 has also been shown to protect against α-synuclein toxicity in a Drosophila model of PD. In the present study we address the potential roles of C-terminal phosphorylation in modulating the interaction of α-synuclein with other protein partners, using a single domain antibody fragment (NbSyn87) that binds to the C-terminal region of α-synuclein with nanomolar affinity. The results reveal that phosphorylation at S129 has negligible effect on the binding affinity of NbSyn87 to α-synuclein while phosphorylation at Y125, only four residues away, decreases the binding affinity by a factor of 400. These findings show that, despite the fact that α-synuclein is intrinsically disordered in solution, selective phosphorylation can modulate significantly its interactions with other molecules and suggest how this particular form of modification could play a key role in regulating the normal and aberrant function of α-synuclein.

Conformational studies of RGDechi peptide by natural-abundance NMR spectroscopy

Integrins are heterodimeric cell-surface proteins that play important roles during developmental and pathological processes. Diverse human pathologies involve integrin adhesion including thrombotic diseases, inflammation, tumour progression, fibrosis, and infectious diseases. Although in the past decade, novel integrin-inhibitor drugs have been developed for integrin-based medical applications, the structural determinants modulating integrin-ligands recognition mechanisms are still poorly understood, reducing the number of integrin subtype exclusive antagonists. In this scenario, we have very recently showed, by means of chemical and biological assays, that a chimeric peptide (named RGDechi), containing a cyclic RGD motif linked to an echistatin C-terminal fragment, is able to interact with the components of integrin family with variable affinities, the highest for α_v β3. Here, in order to understand the mechanistic details driving the molecular recognition mechanism of α_v β₃ by RGDechi, we have performed a detailed structural and dynamics characterization of the free peptide by natural abundance nuclear magnetic resonance (NMR) spectroscopy. Our data indicate that RGDechi presents in solution an heterogeneous conformational ensemble characterized by a more constrained and rigid pentacyclic ring and a largely unstructured acyclic region. Moreover, we propose that the molecular recognition of α_v β₃ integrin by RGDechi occurs by a combination of conformational selection and induced fit mechanisms. Finally, our study indicates that a detailed NMR characterization, by means of natural abundance ¹⁵ N and ¹³ C, of a mostly unstructured bioactive peptide may provide the molecular basis to get essential structural insights into the binding mechanism to the biological partner.

Concomitant disorder and high-affinity zinc binding in the human zinc- and iron-regulated transport protein 4 intracellular loop

The human zinc- and iron-regulated transport protein 4 (hZIP4) protein is the major plasma membrane protein responsible for the uptake of zinc in the body, and as such it plays a key role in cellular zinc homeostasis. hZIP4 plasma membrane levels are regulated through post-translational modification of its large, disordered, histidine-rich cytosolic loop (ICL2) in response to intracellular zinc concentrations. Here, structural characteristics of the isolated disordered loop region, both in the absence and presence of zinc, were investigated using nuclear magnetic resonance (NMR) spectroscopy. NMR chemical shifts, coupling constants and temperature coefficients of the apoprotein, are consistent with a random coil with minor propensities for transient polyproline Type II helices and β-strand in regions implicated in post-translational modifications. The ICL2 protein remains disordered upon zinc binding, which induces exchange broadening. Paramagnetic relaxation enhancement experiments reveal that the histidine-rich region in the apoprotein makes transient tertiary contacts with predicted post-translational modification sites. The residue-specific data presented here strengthen the relationship between hZIP4 post-translational modifications, which impact its role in cellular zinc homeostasis, and zinc sensing by the intracellular loop. Furthermore, the zinc sensing mechanism employed by the ICL2 protein demonstrates that high-affinity interactions can occur in the presence of conformational disorder.

Conformational Changes in the Cytoplasmic Region of KIR3DL1 upon Interaction with SHP-2

KIR3DL1 is an inhibitory killer cell immunoglobulin-like receptor (KIR) that negatively regulates natural killer cell cytotoxicity. The KIR3DL1 cytoplasmic region (3DL1-cyto) is disordered and can be dissected into three segments: (I) H340-V351; (II) M352-D371; and (III) P372-P423. NMR studies indicate that segment II can dynamically adopt a loop-like conformation, and segments I and III can form dynamic helices that may mediate binding to membranes, particularly in the region around the N-terminal (N) immunoreceptor tyrosine-based inhibitory motif (ITIM), consistent with its role in signaling. Furthermore, individual SH2 domains of SHP-2 strongly engage with the unphosphorylated N-ITIM of 3DL1-cyto, while binding of the tandem SHP-2 SH2 domains to the bis-phosphorylated ITIMs results in more extensive conformational changes in segments I and III. The findings enhance our understanding of KIR function and how ITIMs in a target receptor operate in concert to engage the tandem SH2 domains of SHP-2.

Enhanced dynamics of conformationally heterogeneous T7 bacteriophage lysozyme native state attenuates its stability and activity

Proteins are dynamic in nature and exist in a set of equilibrium conformations on various timescale motions. The flexibility of proteins governs various biological functions, and therefore elucidation of such functional dynamics is essential. In this context, we have studied the structure–dynamics–stability–activity relationship of bacteriophage T7 lysozyme/endolysin (T7L) native-state ensemble in the pH range of 6–8. Our studies established that T7L native state is conformationally heterogeneous, as several residues of its C-terminal half are present in two conformations (major and minor) in the slow exchange time scale of nuclear magnetic resonance (NMR). Structural and dynamic studies suggested that the residues belonging to minor conformations do exhibit native-like structural and dynamic features. Furthermore, the NMR relaxation experiments unraveled that the native state is highly dynamic and the dynamic behavior is regulated by the pH, as the pH 6 conformation exhibited enhanced dynamics compared with pH 7 and 8. The stability measurements and cell-based activity studies on T7L indicated that the native protein at pH 6 is ∼2 kcal less stable and is ∼50% less active than those of pH 7 and 8. A comprehensive analysis of the T7L active site, unfolding initiation sites and the residues with altered dynamics outlined that the attenuation of stability and activity is a resultant of its enhanced dynamic properties, which, in turn, can be attributed to the protonation/deprotonation of its partially buried His residues. Our study on T7L structure–dynamics–activity paradigm could assist in engineering novel amidase-based endolysins with enhanced activity and stability over a broad pH range.

A Slice-based 13C-detected NMR Spin System Forming and Resonance Assignment Method

Nuclear magnetic resonance (NMR) spectroscopy is attracting more attention in the field of computational structural biology. Till recently, ¹H-detected experiments are the dominant NMR technique used due to the high sensitivity of ¹H nuclei. However, the current availability of high magnetic fields and cryogenically cooled probe heads allow researchers to overcome the low sensitivity of ¹³C nuclei. Consequently, ¹³C-detected experiments have become a popular technique in different NMR applications especially resonance assignment and structure determination of large proteins. In this paper, we propose the first spin system forming method for ¹³C-detected NMR spectra. Our method is able to accurately form spin systems based on as few as two ¹³C-detected spectra, CBCACON, and CBCANCO. Our method picks slices from the more trusted spectrum and uses them as feedback to direct the slice picking in the less trusted one. This feedback leads to picking the accurate slices that consequently helps to form better spin systems. We tested our method on a real dataset of 'Ubiquitin' and a benchmark simulated dataset consisting of 12 proteins. We fed our spin systems as inputs to a genetic algorithm to generate the chemical shift assignment, and obtained 92 percent correct chemical shift assignment for Ubiquitin. For the simulated dataset, we obtained an average recall of 86 percent and an average precision of 88 percent. Finally, our chemical shift assignment of Ubiquitin was given as an input to CS-ROSETTA server that generated structures close to the experimentally determined structure.

Automatic 13C chemical shift reference correction for unassigned protein NMR spectra

Poor chemical shift referencing, especially for ¹³C in protein Nuclear Magnetic Resonance (NMR) experiments, fundamentally limits and even prevents effective study of biomacromolecules via NMR, including protein structure determination and analysis of protein dynamics. To solve this problem, we constructed a Bayesian probabilistic framework that circumvents the limitations of previous reference correction methods that required protein resonance assignment and/or three-dimensional protein structure. Our algorithm named Bayesian Model Optimized Reference Correction (BaMORC) can detect and correct ¹³C chemical shift referencing errors before the protein resonance assignment step of analysis and without three-dimensional structure. By combining the BaMORC methodology with a new intra-peaklist grouping algorithm, we created a combined method called Unassigned BaMORC that utilizes only unassigned experimental peak lists and the amino acid sequence. Unassigned BaMORC kept all experimental three-dimensional HN(CO)CACB-type peak lists tested within ± 0.4 ppm of the correct ¹³C reference value. On a much larger unassigned chemical shift test set, the base method kept ¹³C chemical shift referencing errors to within ± 0.45 ppm at a 90% confidence interval. With chemical shift assignments, Assigned BaMORC can detect and correct ¹³C chemical shift referencing errors to within ± 0.22 at a 90% confidence interval. Therefore, Unassigned BaMORC can correct ¹³C chemical shift referencing errors when it will have the most impact, right before protein resonance assignment and other downstream analyses are started. After assignment, chemical shift reference correction can be further refined with Assigned BaMORC. These new methods will allow non-NMR experts to detect and correct ¹³C referencing error at critical early data analysis steps, lowering the bar of NMR expertise required for effective protein NMR analysis.

NMR Measurements Reveal the Structural Basis of Transthyretin Destabilization by Pathogenic Mutations

Inherited mutations of transthyretin (TTR) destabilize its structure, leading to aggregation and familial amyloid disease. Although numerous crystal structures of wild-type (WT) and mutant TTRs have been determined, they have failed to yield a comprehensive structural explanation for destabilization by pathogenic mutations. To identify structural and dynamic variations that are not readily observed in the crystal structures, we used NMR to study WT TTR and three kinetically and/or thermodynamically destabilized pathogenic variants (V30M, L55P, and V122I). Sequence-corrected chemical shifts reveal important structural differences between WT and mutant TTR. The L55P mutation linked to aggressive early onset cardiomyopathy and polyneuropathy induces substantial structural perturbations in both the DAGH and CBEF β-sheets, whereas the V30M polyneuropathy-linked substitution perturbs primarily the CBEF sheet. In both variants, the structural perturbations propagate across the entire width of the β-sheets from the site of mutation. Structural changes caused by the V122I cardiomyopathy-associated mutation are restricted to the immediate vicinity of the mutation site, directly perturbing the subunit interfaces. NMR relaxation dispersion measurements show that WT TTR and the three pathogenic variants undergo millisecond time scale conformational fluctuations to populate a common excited state with an altered structure in the subunit interfaces. The excited state is most highly populated in L55P. The combined application of chemical shift analysis and relaxation dispersion to these pathogenic variants reveals differences in ground state structure and in the population of a transient excited state that potentially facilitates tetramer dissociation, providing new insights into the molecular mechanism by which mutations promote TTR amyloidosis.

Molecular dynamics ensemble refinement of the heterogeneous native state of NCBD using chemical shifts and NOEs

Many proteins display complex dynamical properties that are often intimately linked to their biological functions. As the native state of a protein is best described as an ensemble of conformations, it is important to be able to generate models of native state ensembles with high accuracy. Due to limitations in sampling efficiency and force field accuracy it is, however, challenging to obtain accurate ensembles of protein conformations by the use of molecular simulations alone. Here we show that dynamic ensemble refinement, which combines an accurate atomistic force field with commonly available nuclear magnetic resonance (NMR) chemical shifts and NOEs, can provide a detailed and accurate description of the conformational ensemble of the native state of a highly dynamic protein. As both NOEs and chemical shifts are averaged on timescales up to milliseconds, the resulting ensembles reflect the structural heterogeneity that goes beyond that probed, e.g., by NMR relaxation order parameters. We selected the small protein domain NCBD as object of our study since this protein, which has been characterized experimentally in substantial detail, displays a rich and complex dynamical behaviour. In particular, the protein has been described as having a molten-globule like structure, but with a relatively rigid core. Our approach allowed us to describe the conformational dynamics of NCBD in solution, and to probe the structural heterogeneity resulting from both short- and long-timescale dynamics by the calculation of order parameters on different time scales. These results illustrate the usefulness of our approach since they show that NCBD is rather rigid on the nanosecond timescale, but interconverts within a broader ensemble on longer timescales, thus enabling the derivation of a coherent set of conclusions from various NMR experiments on this protein, which could otherwise appear in contradiction with each other.

Structure and dynamics of the extended-helix state of alpha-synuclein: Intrinsic lability of the linker region

The Parkinson's protein alpha-synuclein binds to synaptic vesicles in vivo and adopts a highly extended helical conformation when binding to lipid vesicles in vitro. High-resolution structural analysis of alpha-synuclein bound to small lipid or detergent micelles revealed two helices connected by a non-helical linker, but corresponding studies of the vesicle-bound extended-helix state are hampered by the size and heterogeneity of the protein-vesicle complex. Here we employ fluorinated alcohols (FAs) to induce a highly helical aggregation-resistant state of alpha-synuclein in solution that resembles the vesicle-bound extended-helix state but is amenable to characterization using high-resolution solution-state NMR. Analysis of chemical shift, NOE, coupling constant, PRE and relaxation measurements shows that the lipid-binding domain of alpha-synuclein in FA solutions indeed adopts a single continuous helix and that the ends of this helix do not come into detectable proximity to each other. The helix is well ordered in the center, but features an increase in fast internal motions suggestive of helix fraying approaching the termini. The central region of the helix exhibits slower time scale motions that likely result from flexing of the highly anisotropic structure. Importantly, weak or missing short- and intermediate-range NOEs in the region corresponding to the non-helical linker of micelle-bound alpha-synuclein indicate that the helical structure in this region of the protein is intrinsically unstable. This suggests that conversion of alpha-synuclein from the extended-helix to the broken-helix state represents a functionally relevant structural transition.

Construction of Fluorescent Analogs to Follow the Uptake and Distribution of Cobalamin (Vitamin B12) in Bacteria, Worms, and Plants

Vitamin B₁₂ is made by only certain prokaryotes yet is required by a number of eukaryotes such as mammals, fish, birds, worms, and Protista, including algae. There is still much to learn about how this nutrient is trafficked across the domains of life. Herein, we describe ways to make a number of different corrin analogs with fluorescent groups attached to the main tetrapyrrole-derived ring. A further range of analogs were also constructed by attaching similar fluorescent groups to the ribose ring of cobalamin, thereby generating a range of complete and incomplete corrinoids to follow uptake in bacteria, worms, and plants. By using these fluorescent derivatives we were able to demonstrate that Mycobacterium tuberculosis is able to acquire both cobyric acid and cobalamin analogs, that Caenorhabditis elegans takes up only the complete corrinoid, and that seedlings of higher plants such as Lepidium sativum are also able to transport B₁₂.

A gratuitous β-Lactamase inducer uncovers hidden active site dynamics of the Staphylococcus aureus BlaR1 sensor domain

Increasing evidence shows that active sites of proteins have non-trivial conformational dynamics. These dynamics include active site residues sampling different local conformations that allow for multiple, and possibly novel, inhibitor binding poses. Yet, active site dynamics garner only marginal attention in most inhibitor design efforts and exert little influence on synthesis strategies. This is partly because synthesis requires a level of atomic structural detail that is frequently missing in current characterizations of conformational dynamics. In particular, while the identity of the mobile protein residues may be clear, the specific conformations they sample remain obscure. Here, we show how an appropriate choice of ligand can significantly sharpen our abilities to describe the interconverting binding poses (conformations) of protein active sites. Specifically, we show how 2-(2'-carboxyphenyl)-benzoyl-6-aminopenicillanic acid (CBAP) exposes otherwise hidden dynamics of a protein active site that binds β-lactam antibiotics. When CBAP acylates (binds) the active site serine of the β-lactam sensor domain of BlaR1 (BlaRS), it shifts the time scale of the active site dynamics to the slow exchange regime. Slow exchange enables direct characterization of inter-converting protein and bound ligand conformations using NMR methods. These methods include chemical shift analysis, 2-d exchange spectroscopy, off-resonance ROESY of the bound ligand, and reduced spectral density mapping. The active site architecture of BlaRS is shared by many β-lactamases of therapeutic interest, suggesting CBAP could expose functional motions in other β-lactam binding proteins. More broadly, CBAP highlights the utility of identifying chemical probes common to structurally homologous proteins to better expose functional motions of active sites.

Nonenzymatic acetylation of ubiquitin Lys side chains is modulated by their neighboring residues

Nonenzymatic acetylation of Lys side chains (Lys-SCs) by various in vivo reactive molecules has been suggested to play novel regulatory roles. Ubiquitin (UB) has seven Lys residues that are utilized for synthesis of specific poly-UB chains. To understand the nature of these Lys-SC modifications, the chemical acetylation rate and pK_a and Hill coefficient of each UB-Lys-SC were measured. Mutagenesis studies combined with the determination of activation energy indicated that specific neighboring residues of the Lys-SCs have a potential catalytic activity during nonenzymatic acetylation. Based on the shared chemistry between nonenzymatic Lys acetylation and ubiquitylation, the characterized chemical properties of the UB-Lys-SCs could be a reference for deciphering both mechanisms. Our NMR approaches could be useful for studying general nonenzymatic Lys acylations of various proteins.

Rapid and reliable protein structure determination via chemical shift threading

Protein structure determination using nuclear magnetic resonance (NMR) spectroscopy can be both time-consuming and labor intensive. Here we demonstrate how chemical shift threading can permit rapid, robust, and accurate protein structure determination using only chemical shift data. Threading is a relatively old bioinformatics technique that uses a combination of sequence information and predicted (or experimentally acquired) low-resolution structural data to generate high-resolution 3D protein structures. The key motivations behind using NMR chemical shifts for protein threading lie in the fact that they are easy to measure, they are available prior to 3D structure determination, and they contain vital structural information. The method we have developed uses not only sequence and chemical shift similarity but also chemical shift-derived secondary structure, shift-derived super-secondary structure, and shift-derived accessible surface area to generate a high quality protein structure regardless of the sequence similarity (or lack thereof) to a known structure already in the PDB. The method (called E-Thrifty) was found to be very fast (often < 10 min/structure) and to significantly outperform other shift-based or threading-based structure determination methods (in terms of top template model accuracy)-with an average TM-score performance of 0.68 (vs. 0.50-0.62 for other methods). Coupled with recent developments in chemical shift refinement, these results suggest that protein structure determination, using only NMR chemical shifts, is becoming increasingly practical and reliable. E-Thrifty is available as a web server at http://ethrifty.ca .

Defining the Structural Basis for Allosteric Product Release from E. coli Dihydrofolate Reductase Using NMR Relaxation Dispersion

The rate-determining step in the catalytic cycle of E. coli dihydrofolate reductase is tetrahydrofolate (THF) product release, which can occur via an allosteric or an intrinsic pathway. The allosteric pathway, which becomes accessible when the reduced cofactor NADPH is bound, involves transient sampling of a higher energy conformational state, greatly increasing the product dissociation rate as compared to the intrinsic pathway that obtains when NADPH is absent. Although the kinetics of this process are known, the enzyme structure and the THF product conformation in the transiently formed excited state remain elusive. Here, we use side-chain proton NMR relaxation dispersion measurements, X-ray crystallography, and structure-based chemical shift predictions to explore the structural basis of allosteric product release. In the excited state of the E:THF:NADPH product release complex, the reduced nicotinamide ring of the cofactor transiently enters the active site where it displaces the pterin ring of the THF product. The p-aminobenzoyl-l-glutamate tail of THF remains weakly bound in a widened binding cleft. Thus, through transient entry of the nicotinamide ring into the active site, the NADPH cofactor remodels the enzyme structure and the conformation of the THF to form a weakly populated excited state that is poised for rapid product release.

Unraveling the meaning of chemical shifts in protein NMR

Chemical shifts are among the most informative parameters in protein NMR. They provide wealth of information about protein secondary and tertiary structure, protein flexibility, and protein-ligand binding. In this report, we review the progress in interpreting and utilizing protein chemical shifts that has occurred over the past 25years, with a particular focus on the large body of work arising from our group and other Canadian NMR laboratories. More specifically, this review focuses on describing, assessing, and providing some historical context for various chemical shift-based methods to: (1) determine protein secondary and super-secondary structure; (2) derive protein torsion angles; (3) assess protein flexibility; (4) predict residue accessible surface area; (5) refine 3D protein structures; (6) determine 3D protein structures and (7) characterize intrinsically disordered proteins. This review also briefly covers some of the methods that we previously developed to predict chemical shifts from 3D protein structures and/or protein sequence data. It is hoped that this review will help to increase awareness of the considerable utility of NMR chemical shifts in structural biology and facilitate more widespread adoption of chemical-shift based methods by the NMR spectroscopists, structural biologists, protein biophysicists, and biochemists worldwide. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

A Combined NMR and SAXS Analysis of the Partially Folded Cataract-Associated V75D γD-Crystallin

A cataract is a pathological condition characterized by the clouding of the normally clear eye lens brought about by deposition of crystallin proteins in the lens fiber cells. These protein aggregates reduce visual acuity by scattering or blocking incoming light. Chemical damage to proteins of the crystallin family, accumulated over a lifetime, leads to age-related cataract, whereas inherited mutations are associated with congenital or early-onset cataract. The V75D mutant of γD-crystallin is associated with congenital cataract in mice and was previously shown to un/fold via a partially folded intermediate. Here, we structurally characterized the stable equilibrium urea unfolding intermediate of V75D at the ensemble level using solution NMR and small-angle x-ray scattering. Our data show that, in the intermediate, the C-terminal domain retains a folded conformation that is similar to the native wild-type protein, whereas the N-terminal domain is unfolded and comprises an ensemble of random conformers, without any detectable residual structural propensities.

Structural basis for ligand binding to an enzyme by a conformational selection pathway

Proteins can bind target molecules through either induced fit or conformational selection pathways. In the conformational selection model, a protein samples a scarcely populated high-energy state that resembles a target-bound conformation. In enzymatic catalysis, such high-energy states have been identified as crucial entities for activity and the dynamic interconversion between ground states and high-energy states can constitute the rate-limiting step for catalytic turnover. The transient nature of these states has precluded direct observation of their properties. Here, we present a molecular description of a high-energy enzyme state in a conformational selection pathway by an experimental strategy centered on NMR spectroscopy, protein engineering, and X-ray crystallography. Through the introduction of a disulfide bond, we succeeded in arresting the enzyme adenylate kinase in a closed high-energy conformation that is on-pathway for catalysis. A 1.9-Å X-ray structure of the arrested enzyme in complex with a transition state analog shows that catalytic sidechains are properly aligned for catalysis. We discovered that the structural sampling of the substrate free enzyme corresponds to the complete amplitude that is associated with formation of the closed and catalytically active state. In addition, we found that the trapped high-energy state displayed improved ligand binding affinity, compared with the wild-type enzyme, demonstrating that substrate binding to the high-energy state is not occluded by steric hindrance. Finally, we show that quenching of fast time scale motions observed upon ligand binding to adenylate kinase is dominated by enzyme-substrate interactions and not by intramolecular interactions resulting from the conformational change.

Multi-probe relaxation dispersion measurements increase sensitivity to protein dynamics

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion measurements are a valuable tool for the characterization of structural transitions on the micro-millisecond timescale. While the measurement of (15)N relaxation dispersion is now routine, the measurements with alternative nuclei remain limited. Here we report (15)N as well as (1)H R2 relaxation dispersion measurements of the N23PP/S148A "dynamic knockout" mutant of dihydrofolate reductase. The (1)H dispersion measurements are complementary to (15)N data as many additional residues are observed to have dispersive behavior for the (1)H nucleus. Simultaneous fitting of the dispersion profiles for the two nuclei increases the accuracy of exchange parameters determined for individual residues and clustered groups of residues. The different sensitivity of the two nuclei to changes in backbone torsional angles, ring currents, and hydrogen bonding effects provides important insights into the nature of the structural changes that take place during the exchange process. We observe clear evidence of direct and indirect hydrogen bond effects for the (15)N and (1)H chemical shift changes in the active-site, modulation of ring current shielding in the CD-loop and backbone torsional changes in a cluster of residues associated with the C-terminus. This work demonstrates the power of combined (1)H and (15)N probes for the study of backbone dynamics on the micro-millisecond timescale though the analysis of chemical shift changes.

Natural structural diversity within a conserved cyclic peptide scaffold

We recently isolated and described the evolutionary origin of a diverse class of small single-disulfide bonded peptides derived from Preproalbumin with SFTI-1 (PawS1) proteins in the seeds of flowering plants (Asteraceae). The founding member of the PawS derived peptide (PDP) family is the potent trypsin inhibitor SFTI-1 (sunflower trypsin inhibitor-1) from Helianthus annuus, the common sunflower. Here we provide additional structures and describe the structural diversity of this new class of small peptides, derived from solution NMR studies, in detail. We show that although most have a similar backbone framework with a single disulfide bond and in many cases a head-to-tail cyclized backbone, they all have their own characteristics in terms of projections of side-chains, flexibility and physiochemical properties, attributed to the variety of their sequences. Small cyclic and constrained peptides are popular as drug scaffolds in the pharmaceutical industry and our data highlight how amino acid side-chains can fine-tune conformations in these promising peptides.

The (unusual) aspartic acid in the metal coordination sphere of the prokaryotic zinc finger domain

The possibility of choices of protein ligands and coordination geometries leads to diverse Zn(II) binding sites in zinc-proteins, allowing a range of important biological roles. The prokaryotic Cys2His2 zinc finger domain (originally found in the Ros protein from Agrobacterium tumefaciens) tetrahedrally coordinates zinc through two cysteine and two histidine residues and it does not adopt a correct fold in the absence of the metal ion. Ros is the first structurally characterized member of a family of bacterial proteins that presents several amino acid changes in the positions occupied in Ros by the zinc coordinating residues. In particular, the second position is very often occupied by an aspartic acid although the coordination of structural zinc by an aspartate in eukaryotic zinc fingers is very unusual. Here, by appropriately mutating the protein Ros, we characterize the aspartate role within the coordination sphere of this family of proteins demonstrating how the presence of this residue only slightly perturbs the functional structure of the prokaryotic zinc finger domain while it greatly influences its thermodynamic properties.

Towards a protocol for solution structure determination of copper(II) proteins: the case of Cu(II)Zn(II) superoxide dismutase

We have developed an optimized protocol to solve the solution structure of copper(II) proteins. After assignment, proton-proton NOEs are used for the shell where 1H spectra are conveniently observed. In a shell closer to the metal ion, 13C NMR spectra with band-selective homonuclear decoupling provide the assignment of all nuclei except for those of the metal ligands. A convenient method for the measurement of 13C longitudinal-relaxation rates (R1) of carbonyls and carboxylate moieties is proposed. 1H NOEs and 1H and 13C R1 data are sufficient to produce a good/reasonable solution structure, as demonstrated for a monomeric species of superoxide dismutase, a 153-residue protein.

There is Diversity in Disorder-"In all Chaos there is a Cosmos, in all Disorder a Secret Order"

The protein universe consists of a continuum of structures ranging from full order to complete disorder. As the structured part of the proteome has been intensively studied, stably folded proteins are increasingly well documented and understood. However, proteins that are fully, or in large part, disordered are much less well characterized. Here we collected NMR chemical shifts in a small database for 117 protein sequences that are known to contain disorder. We demonstrate that NMR chemical shift data can be brought to bear as an exquisite judge of protein disorder at the residue level, and help in validation. With the help of secondary chemical shift analysis we demonstrate that the proteins in the database span the full spectrum of disorder, but still, largely segregate into two classes; disordered with small segments of order scattered along the sequence, and structured with small segments of disorder inserted between the different structured regions. A detailed analysis reveals that the distribution of order/disorder along the sequence shows a complex and asymmetric distribution, that is highly protein-dependent. Access to ratified training data further suggests an avenue to improving prediction of disorder from sequence.

Structure-Activity Relationship-based Optimization of Small Temporin-SHf Analogs with Potent Antibacterial Activity

Short antimicrobial peptides represent attractive compounds for the development of new antibiotic agents. Previously, we identified an ultrashort hydrophobic and phenylalanine-rich peptide, called temporin-SHf, representing the smallest natural amphibian antimicrobial peptide known to date. Here, we report on the first structure-activity relationship study of this peptide. A series of temporin-SHf derivatives containing insertion of a basic arginine residue as well as residues containing neutral hydrophilic (serine and α-hydroxymethylserine) and hydrophobic (α-methyl phenylalanine and p-(t)butyl phenylalanine) groups were designed to improve the antimicrobial activity, and their α-helical structure was investigated by circular dichroism and nuclear magnetic resonance spectroscopy. Three compounds were found to display higher antimicrobial activity with the ability to disrupt (permeabilization/depolarization) the bacterial membrane while retaining the nontoxic character of the parent peptide toward rat erythrocytes and human cells (THP-1 derived macrophages and HEK-293). Antimicrobial assays were carried out to explore the influence of serum and physiological salt concentration on peptide activity. Analogs containing d-amino acid residues were also tested. Our study revealed that [p-(t)BuF(2), R(5)]SHf is an attractive ultrashort candidate that is highly potent (bactericidal) against Gram-positive bacteria (including multidrug resistant S. aureus) and against a wider range of clinically interesting Gram-negative bacteria than temporin-SHf, and also active at physiological salt concentrations and in 30% serum.

Identify Beta-Hairpin Motifs with Quadratic Discriminant Algorithm Based on the Chemical Shifts

Successful prediction of the beta-hairpin motif will be helpful for understanding the of the fold recognition. Some algorithms have been proposed for the prediction of beta-hairpin motifs. However, the parameters used by these methods were primarily based on the amino acid sequences. Here, we proposed a novel model for predicting beta-hairpin structure based on the chemical shift. Firstly, we analyzed the statistical distribution of chemical shifts of six nuclei in not beta-hairpin and beta-hairpin motifs. Secondly, we used these chemical shifts as features combined with three algorithms to predict beta-hairpin structure. Finally, we achieved the best prediction, namely sensitivity of 92%, the specificity of 94% with 0.85 of Mathew’s correlation coefficient using quadratic discriminant analysis algorithm, which is clearly superior to the same method for the prediction of beta-hairpin structure from 20 amino acid compositions in the three-fold cross-validation. Our finding showed that the chemical shift is an effective parameter for beta-hairpin prediction, suggesting the quadratic discriminant analysis is a powerful algorithm for the prediction of beta-hairpin.

The molecular mechanism of nuclear transport revealed by atomic-scale measurements

Nuclear pore complexes (NPCs) form a selective filter that allows the rapid passage of transport factors (TFs) and their cargoes across the nuclear envelope, while blocking the passage of other macromolecules. Intrinsically disordered proteins (IDPs) containing phenylalanyl-glycyl (FG)-rich repeats line the pore and interact with TFs. However, the reason that transport can be both fast and specific remains undetermined, through lack of atomic-scale information on the behavior of FGs and their interaction with TFs. We used nuclear magnetic resonance spectroscopy to address these issues. We show that FG repeats are highly dynamic IDPs, stabilized by the cellular environment. Fast transport of TFs is supported because the rapid motion of FG motifs allows them to exchange on and off TFs extremely quickly through transient interactions. Because TFs uniquely carry multiple pockets for FG repeats, only they can form the many frequent interactions needed for specific passage between FG repeats to cross the NPC.

DOI:http://dx.doi.org/10.7554/eLife.10027.001

Eukaryotic cells have a nucleus that contains most of the organism's genetic material. Two layers of membrane form an envelope around the nucleus and protect its contents from the rest of the cell's interior. However, this protective barrier must also allow certain proteins and nucleic acids(collectively called ‘cargo’) to move in and out of the nucleus.

Cargo molecules can pass through channel-like structures called nuclear pore complexes, which are embedded in the nuclear envelope. However, transport across this barrier is highly selective. While small molecules can pass freely through nuclear pore complexes, larger cargo can only be transported when they are bound to so-called transport factors. The nuclear pore complex is a large structure made up of more than 30 different proteins called nucleoporins. Like all proteins, nucleoporins are built from amino acids. Many nucleoporins contain repeating units of two amino acids, namely phenylalanine (which is often referred to as ‘F’) and glycine (or ‘G’). These ‘FG nucleoporins’ are found on the inside of the nuclear pore complex and interact with transport factors to allow them to transit across the nuclear envelope.

Several models have been put forward to explain how FG nucleoporins block the passage of most molecules. But it was unclear from these models how these nucleoporins could do this while simultaneously allowing the selective and fast transport of nuclear transport receptors. There was also a major lack of experimental data that probed the behavior of FG nucleoporins in detail.

Hough, Dutta et al. have now used a technique called nuclear magnetic resonance spectroscopy (or NMR for short) to address this issue. NMR can be used to analyze the structure of proteins and how they interact with other molecules. This analysis revealed that FG nucleoporins never adopt an ordered three-dimensional shape, even briefly; instead they remain unfolded or disordered, moving constantly. Nevertheless, and unlike many other unfolded proteins, FG nucleoporins do not aggregate into clumps. This is because they are constantly changing and continuously interacting with other molecules present inside the cell, which prevents them from aggregating.

Hough, Dutta et al. also observed that the repeating units in the FG nucleoporins engaged briefly with a large number of sites or pockets present on the transport factors. These FG repeats can bind and then release the transport factors at unusually high speeds, which enables the transport factors to move quickly through the nuclear pore complex. This transit is specific because only transport factors have a high capacity for interacting with the FG repeats. These findings provide an explanation for how the nuclear pore complex achieves fast and selective transport. Further work is needed to see whether certain FG nucleoporins specifically interact with a particular type of transport factor, to provide preferred transport routes through the nuclear pore complex.

DOI:http://dx.doi.org/10.7554/eLife.10027.002