Rapid measurement of 3J(H N-H alpha) and 3J(N-H beta) coupling constants in polypeptides

1Hβ chemical shift-based phase modulated NMR methods for fast identification of amino acid types in proteins

Here we report two phase modulated NMR experiments: PM-2D HN(CACBHB) and PM-2D HN(HB), that use 1Hβ chemical shifts to rapidly identify amino acid type in proteins. The magnetization on the 1Hβ spins during the experiments is allowed to evolve for a fixed evolution period that results in phase modulation (positive or negative) of the cross peaks corresponding to various amino acid residues on their 2D HN projections, resembling a typical 2D [1H-15N]-HSQC spectrum. All amino acids except glycine can be categorized into three discernible groups based on their 1Hβ chemical shifts, resulting in unique phase patterns at different fixed evolution periods for 1Hβ, thus facilitating their identification. Remarkably, the PM-2D HN(HB) stands out among all amino acid type identification NMR techniques for its applicability with cost-effective and most routinely employed 15N-labeled protein samples for NMR studies. Furthermore, when combined effectively with the 13Cβ chemical shift-based phase modulated NMR method (PM-2D HN(CACB)), these methods resolved the identification of large groups of amino acids into relatively smaller groups. Moreover, these techniques can accelerate the sequence-specific sequential resonance assignment (SSRA) process and would help in fast tracking of assigned NMR signals exhibiting chemical shift perturbation on the 2D [1H-15N]-HSQC spectrum of proteins during various experiments (e.g., temperature change, pH change, and protein or ligand or cofactor binding) as well as in site-directed mutagenesis.

IGF-dependent dynamic modulation of a protease cleavage site in the intrinsically disordered linker domain of human IGFBP2

Functional regulation via conformational dynamics is well known in structured proteins but less well characterized in intrinsically disordered proteins and their complexes. Using NMR spectroscopy, we have identified a dynamic regulatory mechanism in the human insulin-like growth factor (IGF) system involving the central, intrinsically disordered linker domain of human IGF-binding protein-2 (hIGFBP2). The bioavailability of IGFs is regulated by the proteolysis of IGF-binding proteins. In the case of hIGFBP2, the linker domain (L-hIGFBP2) retains its intrinsic disorder upon binding IGF-1, but its dynamics are significantly altered, both in the IGF binding region and distantly located protease cleavage sites. The increase in flexibility of the linker domain upon IGF-1 binding may explain the IGF-dependent modulation of proteolysis of IGFBP2 in this domain. As IGF homeostasis is important for cell growth and function, and its dysregulation is a key contributor to several cancers, our findings open up new avenues for the design of IGFBP analogs inhibiting IGF-dependent tumors.

Predicting Genetic Variation Severity Using Machine Learning to Interpret Molecular Simulations

Distinct missense mutations in a specific gene have been associated with different diseases as well as differing severity of a disease. Current computational methods predict the potential pathogenicity of a missense variant but fail to differentiate between separate disease or severity phenotypes. We have developed a method to overcome this limitation by applying machine learning to features extracted from molecular dynamics simulations, creating a way to predict the effect of novel genetic variants in causing a disease, drug resistance, or another specific trait. As an example, we have applied this novel approach to variants in calmodulin associated with two distinct arrhythmias as well as two different neurodegenerative diseases caused by variants in amyloid-β peptide. The new method successfully predicts the specific disease caused by a gene variant and ranks its severity with more accuracy than existing methods. We call this method molecular dynamics phenotype prediction model.

Comparative structure, dynamics and evolution of acyl-carrier proteins from Borrelia burgdorferi, Brucella melitensis and Rickettsia prowazekii

Acyl carrier proteins (ACPs) are small helical proteins found in all kingdoms of life, primarily involved in fatty acid and polyketide biosynthesis. In eukaryotes, ACPs are part of the fatty acid synthase (FAS) complex, where they act as flexible tethers for the growing lipid chain, enabling access to the distinct active sites in FAS. In the type II synthesis systems found in bacteria and plastids, these proteins exist as monomers and perform various processes, from being a donor for synthesis of various products such as endotoxins, to supplying acyl chains for lipid A and lipoic acid FAS (quorum sensing), but also as signaling molecules, in bioluminescence and activation of toxins. The essential and diverse nature of their functions makes ACP an attractive target for antimicrobial drug discovery. Here, we report the structure, dynamics and evolution of ACPs from three human pathogens: Borrelia burgdorferi, Brucella melitensis and Rickettsia prowazekii, which could facilitate the discovery of new inhibitors of ACP function in pathogenic bacteria.

NMR structure and dynamics of inhibitory repeat domain variant 12, a plant protease inhibitor from Capsicum annuum, and its structural relationship to other plant protease inhibitors

Although several plant protease inhibitors have been structurally characterized using X-ray crystallography, very few have been studied using NMR techniques. Here, we report an NMR study of the solution structure and dynamics of an inhibitory repeat domain (IRD) variant 12 from the wound-inducible Pin-II type proteinase inhibitor from Capsicum annuum. IRD variant 12 (IRD12) showed strong anti-metabolic activity against the Lepidopteran insect pest, Helicoverpa armigera. The NMR-derived three-dimensional structure of IRD12 reveals a three-stranded anti-parallel β-sheet rigidly held together by four disulfide bridges and shows structural homology with known IRDs. It is interesting to note that the IRD12 structure containing ∼75% unstructured part still shows substantial amount of rigidity of N-H bond vectors with respect to its molecular motion.Communicated by Ramaswamy H. Sarma.

1H, 13C and 15N NMR assignments of two plant protease inhibitors (IRD7 and IRD12) from the plant Capsicum annuum

Helicoverpa species are polyphagous pests, with the larval stages causing major damage to economically valuable crops such as cotton, tomato, corn, sorghum, peas, sunflower, wheat and other pulses. Over the years, Helicoverpa armigera has developed resistance to most classes of chemical insecticides, and consequently it is now largely controlled on cotton plants via the use of Bt transgenic crops that express insecticidal Cry toxins which in-turn expedited resistance development in a number of pest species including H. armigera. In a hope to provide other eco-friendly alternatives solutions to counter the effect of the pest, people have identified a number of protease inhibitors (PIs) from the domesticated capsicum species Capsicum annuum, several of which potently inhibited H. armigera gut proteases and impeded growth of H. armigera larva. With a view to explore and enhance the specific nature or properties of these PIs on the mechanism of inhibition, structural and functional characterization of these PIs are inevitable. Towards this goal, we have carried out complete ¹H, ¹³C and ¹⁵N resonance assignments of two of these PIs, identified as IRD7 and IRD12, using a suite of 2D and 3D multi-dimensional and multi-nuclear NMR experiments.

Is the Conformational Ensemble of Alzheimer's Aβ10-40 Peptide Force Field Dependent?

By applying REMD simulations we have performed comparative analysis of the conformational ensembles of amino-truncated Aβ10-40 peptide produced with five force fields, which combine four protein parameterizations (CHARMM36, CHARMM22*, CHARMM22/cmap, and OPLS-AA) and two water models (standard and modified TIP3P). Aβ10-40 conformations were analyzed by computing secondary structure, backbone fluctuations, tertiary interactions, and radius of gyration. We have also calculated Aβ10-40 ³J_HNHα-coupling and RDC constants and compared them with their experimental counterparts obtained for the full-length Aβ1-40 peptide. Our study led us to several conclusions. First, all force fields predict that Aβ adopts unfolded structure dominated by turn and random coil conformations. Second, specific TIP3P water model does not dramatically affect secondary or tertiary Aβ10-40 structure, albeit standard TIP3P model favors slightly more compact states. Third, although the secondary structures observed in CHARMM36 and CHARMM22/cmap simulations are qualitatively similar, their tertiary interactions show little consistency. Fourth, two force fields, OPLS-AA and CHARMM22* have unique features setting them apart from CHARMM36 or CHARMM22/cmap. OPLS-AA reveals moderate β-structure propensity coupled with extensive, but weak long-range tertiary interactions leading to Aβ collapsed conformations. CHARMM22* exhibits moderate helix propensity and generates multiple exceptionally stable long- and short-range interactions. Our investigation suggests that among all force fields CHARMM22* differs the most from CHARMM36. Fifth, the analysis of ³J_HNHα-coupling and RDC constants based on CHARMM36 force field with standard TIP3P model led us to an unexpected finding that in silico Aβ10-40 and experimental Aβ1-40 constants are generally in better agreement than these quantities computed and measured for identical peptides, such as Aβ1-40 or Aβ1-42. This observation suggests that the differences in the conformational ensembles of Aβ10-40 and Aβ1-40 are small and the former can be used as proxy of the full-length peptide. Based on this argument, we concluded that CHARMM36 force field with standard TIP3P model produces the most accurate representation of Aβ10-40 conformational ensemble.

Dependence of protein conformational ensembles on force field parameterizations limits the predictive power of molecular dynamics simulations. To address this problem, we evaluated five all-atom force fields for their consistency in reproducing the conformational ensemble of Alzheimer’s Aβ10-40 peptide. To generate conformational ensembles, we have used replica exchange molecular dynamics and computed Aβ10-40 secondary and tertiary structures. We found that, although all force fields predict Aβ10-40 unfolded structure, they strongly disagree on helix and β propensities and tertiary structure distributions. We have also calculated Aβ10-40 J-coupling and residual dipolar coupling constants and compared them with the experimental data for the full-length Aβ1-40 peptide. Unexpectedly, we determined that in silico Aβ10-40 and experimental Aβ1-40 constants are in better agreement than these quantities computed and measured previously for identical peptides, such as Aβ1-40 or Aβ1-42. We then concluded that the conformational ensembles of Aβ10-40 and Aβ1-40 are similar and on this basis argue that CHARMM36 force field with standard TIP3P water model provides the most accurate description of Aβ10-40. Although our objective was not to evaluate the biomolecular force fields in general, our study is expected to facilitate their proper selection for the simulations of Alzheimer’s peptides.

1H, 13C and 15N NMR assignments of a calcium-binding protein from Entamoeba histolytica

We report almost complete sequence specific (1)H, (13)C and (15)N NMR assignments of a 150-residue long calmodulin-like calcium-binding protein from Entamoeba histolytica (EhCaBP6), as a prelude to its structural and functional characterization.

(1)H, (13)C and (15)N NMR assignments of Mg (2+) bound form of UV inducible transcript protein (UVI31+) from Chlamydomonas reinhardtii

Almost complete sequence specific (1)H, (13)C and (15)N resonance assignments of Mg(2+) bound form of UV inducible transcript protein (UVI31+) from Chlamydomonas reinhardtii are reported, as a prelude to its structural and functional characterization.

Functional manipulation of a calcium-binding protein from Entamoeba histolytica guided by paramagnetic NMR

EhCaBP1, one of the calcium-binding proteins from Entamoeba histolytica, is a two-domain EF-hand protein. The two domains of EhCaBP1 are structurally and functionally different from each other. However, both domains are required for structural stability and a full range of functional diversity. Analysis of sequence and structure of EhCaBP1 and other CaBPs indicates that the C-terminal domain of EhCaBP1 possesses a unique structure compared with other family members. This had been attributed to the absence of a Phe-Phe interaction between highly conserved Phe residues at the -4 position in EF-hand III (F[-4]; Tyr(81)) and at the 13th position in EF-hand IV (F[+13]; Phe(129)) of the C-terminal domain. Against this backdrop, we mutated the Tyr residue at the -4th position of EF III to the Phe residue (Y81F), to bring in the Phe-Phe interaction and understand the nature of structural and functional changes in the protein by NMR spectroscopy, molecular dynamics (MD) simulation, isothermal titration calorimetry (ITC), and biological assays, such as imaging and actin binding. The Y81F mutation in EhCaBP1 resulted in a more compact structure for the C-terminal domain of the mutant as in the case of calmodulin and troponin C. The compact structure is favored by the presence of a π-π interaction between Phe(81) and Phe(129) along with several hydrophobic interactions of Phe(81), which are not seen in the wild-type protein. Furthermore, the biological assays reveal preferential membrane localization of the mutant, loss of its colocalization with actin in the phagocytic cups, whereas retaining its ability to bind G- and F-actin.

Guanidine-HCl dependent structural unfolding of M-crystallin: fluctuating native state like topologies and intermolecular association

Numerous experimental techniques and computational studies, proposed in recent times, have revolutionized the understanding of protein-folding paradigm. The complete understanding of protein folding and intermediates are of medical relevance, as the aggregation of misfolding proteins underlies various diseases, including some neurodegenerative disorders. Here, we describe the unfolding of M-crystallin, a βγ-crystallin homologue protein from archaea, from its native state to its denatured state using multidimensional NMR and other biophysical techniques. The protein, which was earlier characterized to be a predominantly β-sheet protein in its native state, shows different structural propensities (α and β), under different denaturing conditions. In 2 M GdmCl, the protein starts showing two distinct sets of peaks, with one arising from a partially unfolded state and the other from a completely folded state. The native secondary structural elements start disappearing as the denaturant concentration approaches 4 M. Subsequently, the protein is completely unfolded when the denaturant concentration is 6 M. The (15)N relaxation data (T(1)/T(2)), heteronuclear (1)H-(15)N Overhauser effects (nOes), NOESY data, and other biophysical data taken together indicate that the protein shows a consistent, gradual change in its structural and motional preferences with increasing GdmCl concentration.

NMR structure of an acyl-carrier protein from Borrelia burgdorferi

Nearly complete resonance assignment and the high-resolution NMR structure of the acyl-carrier protein from Borrelia burgdorferi, a target of the Seattle Structural Genomics Center for Infectious Disease (SSGCID) structure-determination pipeline, are reported. This protein was chosen as a potential target for drug-discovery efforts because of its involvement in fatty-acid biosynthesis, an essential metabolic pathway, in bacteria. It was possible to assign >98% of backbone resonances and >92% of side-chain resonances using multidimensional NMR spectroscopy. The NMR structure was determined to a backbone r.m.s.d. of 0.4 Å and contained four α-helices and two 3(10)-helices. A structure-homology search revealed that this protein is highly similar to the acyl-carrier protein from Aquifex aeolicus.

Temperature-dependent oligomerization in M-crystallin: lead or lag toward cataract, an NMR perspective

The oligomerization and/or aggregation of proteins is of critical importance in a wide variety of biomedical situations, ranging from abnormal disease states like Alzheimer's and Parkinson's disease to the production of inclusion bodies, stability, and delivery of protein drugs. In the case of eye-lens proteins, oligomerization is implicated in cataract formation. In the present study, we have investigated the temperature driven oligomerization of M-crystallin, a close homologue of eye-lens proteins, using NMR spectroscopy and dynamic-light scattering (DLS). The NMR data primarily included R(1), R(2) relaxation rates and nOes of the backbone amide groups recorded at three different temperatures, 25, 20, and 15° C. The major outcome of the study is the two fold increase in the overall tumbling time (τ(c)) of M-crystallin on lowering the temperature from 25 to 15° C. An extrapolation of τ(c) to a further lower temperature (5° C) may lead to a τ(c) of ∼19 ns that would correspond to a τ(c) value of a tetrameric M-crystallin. These results also validate the observed changes in the hydrodynamic radius of M-crystallin, determined using DLS data. Further, the temperature-dependent protein dynamics of M-crystallin reveal considerable variation at/near the Ca(2+)-binding sites. A concerted analysis of the temperature dependent relaxation parameters and DLS data reveals that the self-association of the protein is not only a monomer-dimer equilibrium, but also goes to tetramers or other multimeric states. These higher states may co-exist in fast exchange with the monomeric and dimeric M-crystallin at milli-molar to sub-millimolar concentrations and at lower temperature.

Synthesis and reactivity of (benzoxazol-2-ylmethyl)phosphonic acid

An efficient three step synthesis of (benzoxazol-2-ylmethyl)phosphonic acid (6-H(2)) is described along with IR, mass spectrometry (MS), and (1)H, (13)C, and (31)P NMR spectroscopic characterization data, and a single crystal X-ray diffraction structure determination. 6-H(2) is unstable in acidic aqueous solutions (pH < 4) undergoing ring-opening to give [(2-hydroxyphenylcarbamoyl)methyl] phosphonic acid (7-H(2)) that is characterized by IR, MS, and NMR methods. The protonation constants (pK(a)) for 7-H(2) have been measured, and crystal structure determinations for (NH(4))(7-H) and K(7-H)·DMF are described. Reactions of NaOH and KOH with 6-H(2) in MeOH/H(2)O solutions led to isolation and crystal structure determinations of the salts [Na(6-H)·H(2)O](2), K(6-H), Na(3)(6)(6-H)·H(2)O, and [K(2)(6)](2)·3H(2)O. The complexation reactions of 7-H(2) with La(III), Nd(III), and Gd(III), as a function of pH, were also examined by titrametric methods, and a model for the 1:1 anion binding with Ln(III) cations is proposed.

Calmodulin-like protein from Entamoeba histolytica: solution structure and calcium-binding properties of a partially folded protein

The mechanism of Ca(2+)-signaling in the protozoan parasite Entamoeba histolytica is yet to be understood as many of the key regulators are still to be identified. E. histolytica encodes a number of multi-EF-hand Ca(2+)-binding proteins (EhCaBPs). Functionally only one of these molecules, EhCaBP1, has been characterized to date. The calmodulin-like protein from E. histolytica (abbreviated as EhCaM or EhCaBP3) is a 17.23 kDa monomeric protein that shows maximum sequence identity with heterologous calmodulins (CaMs). Though CaM activity has been biochemically shown in E. histolytica, there are no reports on the presence of a typical CaM. In an attempt to understand the structural and functional similarity of EhCaM with CaM, we have determined the three-dimensional (3D) solution structure of EhCaM using NMR. The EhCaM has a well-folded N-terminal domain and an unstructured C-terminal counterpart. Further, it sequentially binds only two calcium ions, an unusual mode of Ca(2+)-binding among the known CaBPs, notably both in the N-terminal domain of EhCaM. Further, EhCaM is present in the nucleus in addition to the cytoplasm as detected by immunofluorescence staining, unlike other EhCaBPs that are detected only in the cytoplasm. Therefore, this protein is likely to have a different function. The presence of unusual and a diverse set of CaBPs in E. histolytica suggests a distinct Ca(2+)-signaling process in E. histolytica. The results reported here help in understanding the structure-function relationship of CaBPs including their Ca(2+)-binding properties.

Root-mean-square-deviation-based rapid backbone resonance assignments in proteins

We have shown that the methodology based on the estimation of root-mean-square deviation (RMSD) between two sets of chemical shifts is very useful to rapidly assign the spectral signatures of (1)H(N), (13)C(α), (13)C(β), (13)C', (1)H(α) and (15)N spins of a given protein in one state from the knowledge of its resonance assignments in a different state, without resorting to routine established procedures (manual and automated). We demonstrate the utility of this methodology to rapidly assign the 3D spectra of a metal-binding protein in its holo-state from the knowledge of its assignments in apo-state, the spectra of a protein in its paramagnetic state from the knowledge of its assignments in diamagnetic state and, finally, the spectra of a mutant protein from the knowledge of the chemical shifts of the corresponding wild-type protein. The underlying assumption of this methodology is that, it is impossible for any two amino acid residues in a given protein to have all the six chemical shifts degenerate and that the protein under consideration does not undergo large conformational changes in going from one conformational state to another. The methodology has been tested using experimental data on three proteins, M-crystallin (8.5 kDa, predominantly β-sheet, for apo- to holo-state), Calbindin (7.5 kDa, predominantly α-helical, for diamagnetic to paramagnetic state and apo to holo) and EhCaBP1 (14.3 kDa, α-helical, the wild-type protein with one of its mutant). In all the cases, the extent of assignment is found to be greater than 85%.

Aliasing in reduced dimensionality NMR spectra: (3,2)D HNHA and (4,2)D HN(COCA)NH experiments as examples

Reduced dimensionality NMR spectra usually require very large spectral widths in the shared dimension. In this paper we show that aliasing can be introduced in reduced dimensionality NMR spectra either to decrease the acquisition time or increase the resolution of the experiments without losing information. The gains of introducing aliasing are illustrated with two examples, the (3,2)D HNHA and the (4,2)D HN(COCA)NH experiments. In both cases a reduction of the spectral width of more than 50% was possible.

CSSI-PRO: a method for secondary structure type editing, assignment and estimation in proteins using linear combination of backbone chemical shifts

Estimation of secondary structure in polypeptides is important for studying their structure, folding and dynamics. In NMR spectroscopy, such information is generally obtained after sequence specific resonance assignments are completed. We present here a new methodology for assignment of secondary structure type to spin systems in proteins directly from NMR spectra, without prior knowledge of resonance assignments. The methodology, named Combination of Shifts for Secondary Structure Identification in Proteins (CSSI-PRO), involves detection of specific linear combination of backbone (1)H(alpha) and (13)C' chemical shifts in a two-dimensional (2D) NMR experiment based on G-matrix Fourier transform (GFT) NMR spectroscopy. Such linear combinations of shifts facilitate editing of residues belonging to alpha-helical/beta-strand regions into distinct spectral regions nearly independent of the amino acid type, thereby allowing the estimation of overall secondary structure content of the protein. Comparison of the predicted secondary structure content with those estimated based on their respective 3D structures and/or the method of Chemical Shift Index for 237 proteins gives a correlation of more than 90% and an overall rmsd of 7.0%, which is comparable to other biophysical techniques used for structural characterization of proteins. Taken together, this methodology has a wide range of applications in NMR spectroscopy such as rapid protein structure determination, monitoring conformational changes in protein-folding/ligand-binding studies and automated resonance assignment.

Complete backbone assignment of a Ca2+-binding protein of the betagamma-crystallin superfamily from Methanosarcina acetivorans, at two denaturant concentrations

We report here almost complete backbone assignment of a Ca(2+)-binding protein of the betagamma-crystallin superfamily from Methanosarcina acetivorans, at two denaturant (GdmCl) concentrations, using double and triple resonance experiments. These NMR assignments will be useful to understand the unfolding path of this protein.

Chemical shift based editing of CH3 groups in fractionally 13C-labelled proteins using GFT (3, 2)D CT-HCCH-COSY: stereospecific assignments of CH3 groups of Val and Leu residues

We propose a (3, 2)D CT-HCCH-COSY experiment to rapidly collect the data and provide significant dispersion in the spectral region containing (13)C-(1)H cross peaks of CH(3) groups belonging to Ala, Ile, Leu, Met, Thr and Val residues. This enables one to carry out chemical shift based editing and grouping of all the (13)C-(1)H cross peaks of CH(3) groups belonging to Ala, Ile, Leu, Met, Thr and Val residues in fractionally (10%) (13)C-labelled proteins, which in turn aids in the sequence-specific resonance assignments in general and side-chain resonance assignments in particular, in any given protein. Further, we demonstrate the utility of this experiment for stereospecific assignments of the pro-R and pro-S methyl groups belonging to the Leu and Val residues in fractionally (10%) (13)C-labelled proteins. The proposed experiment opens up a wide range of applications in resonance assignment strategies and structure determination of proteins.

Sequence specific 1H, 13C and 15N resonance assignments of a calmodulin-like calcium-binding protein from the protozoan parasite Entamoeba histolytica (EhCaM)

We report almost complete sequence specific 1H, 13C and 15N NMR assignments of a 151-residue long calmodulin-like calcium-binding protein from Entamoeba histolytica (EhCaM).