Structural Interpretation of Paramagnetic Relaxation Enhancement-Derived Distances for Disordered Protein States

Paramagnetic relaxation enhancement (PRE) is a powerful technique for studying transient tertiary organizations of unfolded and partially folded proteins. The heterogeneous and dynamic nature of disordered protein states, together with the r⁻⁶ dependence of PRE, presents significant challenges for reliable structural interpretation of PRE-derived distances. Without additional knowledge of accessible conformational substates, ensemble-simulation-based protocols have been used to calculate structure ensembles that appear to be consistent with the PRE distance restraints imposed on the ensemble level with the proper r⁻⁶ weighting. However, rigorous assessment of the reliability of such protocols has been difficult without intimate knowledge of the true nature of disordered protein states. Here we utilize sets of theoretical PRE distances derived from simulated structure ensembles that represent the folded, partially folded and unfolded states of a small protein to investigate the efficacy of ensemble-simulation-based structural interpretation of PRE distances. The results confirm a critical limitation that, due to r⁻⁶ weighting, only one or a few members need to satisfy the distance restraints and the rest of the ensemble are essentially unrestrained. Consequently, calculated structure ensembles will appear artificially heterogeneous no matter whether the PRE distances are derived from the folded, partially unfolded or unfolded state. Furthermore, the nature of the heterogeneous ensembles is largely determined by the protein model employed in structure calculation and reflects little on the true nature of the underlying disordered state. These findings suggest that PRE measurements on disordered protein states alone generally do not contain enough information for a reliable structural interpretation and that the latter will require additional knowledge of accessible conformational substates. Interestingly, when a very large number of PRE measurements is available, faithful structural interpretation might be possible with intermediate ensemble sizes under ideal conditions.     

Ionic strength‐dependent conformations of a ubiquitin‐like small archaeal modifier protein (SAMP1) from Haloferax volcanii

Eukaryotic ubiquitin and ubiquitin‐like systems play crucial roles in various cellular biological processes. In this work, we determined the solution structure of SAMP1 from Haloferax volcanii by NMR spectroscopy. Under low ionic conditions, SAMP1 presented two distinct conformations, one folded β‐grasp and the other disordered. Interestingly, SAMP1 underwent a conformational conversion from disorder to order with ion concentration increasing, indicating that the ordered conformation is the functional form of SAMP1 under the physiological condition of H. volcanii. Furthermore, SAMP1 could interact with proteasome‐activating nucleotidase B, supposing a potential role of SAMP1 in the protein degradation pathway mediated by proteasome.     

Dynamic aspects of pressure and temperature-stabilized intermediates of outer surface protein A

Structural characterization of alternatively folded and partially disordered protein conformations remains challenging. Outer surface protein A (OspA) is a pivotal protein in Borrelia infection, which is the etiological agent of Lyme disease. OspA exists in equilibrium with intermediate conformations, in which the central and the C-terminal regions of the protein have lower stabilities than the N-terminal. Here, we characterize pressure- and temperature-stabilized intermediates of OspA by nuclear magnetic resonance spectroscopy combined with paramagnetic relaxation enhancement (PRE). We found that although the C-terminal region of the intermediate was partially disordered, it retains weak specific contact with the N-terminal region, owing to a twist of the central β-sheet and increased flexibility in the polypeptide chain. The disordered C-terminal region of the pressure-stabilized intermediate was more compact than that of the temperature-stabilized form. Further, molecular dynamics simulation demonstrated that temperature-induced disordering of the β-sheet was initiated at the C-terminal region and continued through to the central region. An ensemble of simulation snapshots qualitatively described the PRE data from the intermediate and indicated that the intermediate structures of OspA may expose tick receptor-binding sites more readily than does the basic folded conformation.     

VCD Robustness of the Amide‐I and Amide‐II Vibrational Modes of Small Peptide Models

The rotational strengths and the robustness values of amide‐I and amide‐II vibrational modes of For(AA)_nNHMe (where AA is Val, Asn, Asp, or Cys, n = 1–5 for Val and Asn; n = 1 for Asp and Cys) model peptides with α‐helix and β‐sheet backbone conformations were computed by density functional methods. The robustness results verify empirical rules drawn from experiments and from computed rotational strengths linking amide‐I and amide‐II patterns in the vibrational circular dichroism (VCD) spectra of peptides with their backbone structures. For peptides with at least three residues (n ≥ 3) these characteristic patterns from coupled amide vibrational modes have robust signatures. For shorter peptide models many vibrational modes are nonrobust, and the robust modes can be dependent on the residues or on their side chain conformations in addition to backbone conformations. These robust VCD bands, however, provide information for the detailed structural analysis of these smaller systems. Chirality 27:625–634, 2015 © 2015 Wiley Periodicals, Inc.     

Molecular Conformation of Alkaline Proteinase of Aspergillus sojae in Aqueous Solution

The molecular conformation of the alkaline proteinase of Aspergillus sojae in aqueous solution was investigated by the optical rotatory dispersion, Cotton effects, infra-red absorption spectra (amide I and V bands), ultraviolet difference spectra, etc. It is concluded that; (1) there are about 10 to 15% of the α-helix and a small amount of the β-structure in the enzyme molecule, but most parts of the peptide chain are present as the disordered structure; (2) the enzyme molecule is compactly folded even in the disordered parts; and (3) the two tryptophan residues involved in the peptide chain are burried in the interior of the molecule.     

Solvent effects on the conformational preferences of model peptoids. MP2 study

The influence of aqueous environment on the main‐chain conformation (ω₀, ?, and ψ dihedral angles) of two model peptoids: N‐acetyl‐N‐methylglycine N’‐methylamide (Ac‐N(Me)‐Gly‐NHMe) ( 1 ) and N‐acetyl‐N‐methylglycine N’,N’‐dimethylamide (Ac‐N(Me)‐Gly‐NMe₂) ( 2 ) was investigated by MP2/6‐311++G(d,p) method. The Ramachandran maps of both studied molecules with cis and trans configuration of the N‐terminal amide bond in the gas phase and in water environment were obtained and all energy minima localized. The polarizable continuum model was applied to estimate the solvation effect on conformation. Energy minima of the Ac‐N(Me)‐Gly‐NHMe and Ac‐N(Me)‐Gly‐NMe₂ have been analyzed in terms of the possible hydrogen bonds and C = O dipole attraction. To validate the theoretical results obtained, conformations of the similar structures gathered in the Cambridge Crystallographic Data Centre were analyzed. Obtained results indicate that aqueous environment in model peptoids 1 and 2 favors the conformation F (? and ψ = ?70º, 180º), and additionally significantly increases the percentage of structures with cis configuration of N‐terminal amide bond in studied compounds. Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

A conformational study of the tetrapeptide CH3CO-Ala-Asp-Gly-Lys-NHCH3 corresponding to a β-bend in staphylococcal nuclease

The tetrapeptide sequence Ala-Asp-Gly-Lys occurs as a type I′ β-bend at residues 94–97 in staphylococcal nuclease. We have synthesized theN-acetyl,N′-methylamide derivative of this tetrapeptide and studied its conformation in solution, using nuclear magnetic resonance (NMR) and circular dichroism (CD) spectroscopy. In the synthesis, special attention was paid to the possibility of cyclic aspartimide formation giving rise to mixtures of α- and β-Asp-Gly products. The presence of such a mixture was excluded by infrared, NMR, and other analytical procedures applied to the products and to models for α- and β-linked aspartyl residues. The CD spectra of the protected tetrapeptide in water, methanol, and trifluoroethanol show no evidence of preferred chain conformations. In dimethylsulfoxide-d ₆ , however, the NMR spectra are consistent with the presence of a population of conformers in which the Lys and C-terminal NHCH₃ amide protons are shielded from solvent. Taken together with the observed³J_NH-C ^α _H coupling constants for all residues, this permitted the construction and energetic evaluation of possible conformations in solution. Only one such conformation was fully compatible with the NMR data; this is a type II β-bend in which the Lys and C-terminal NHCH₃ amide protons are close to the Ala C=O group and may form bifurcated hydrogen bonds with it. This conformation can be converted into the conformation existing in staphylococcal nuclease by rotating the plane of the Ala-Asp peptide group by about 120° around a line connecting the Ala and Asp C_α atoms and by making small shifts in dihedral angles elsewhere in the peptide.     

Assessing the native state conformational distribution of ubiquitin by peptide acidity

At equilibrium, every energetically feasible conformation of a protein occurs with a non-zero probability. Quantitative analysis of protein flexibility is thus synonymous with determining the proper Boltzmann-weighting of this conformational distribution. The exchange reactivity of solvent-exposed amide hydrogens greatly varies with conformation, while the short-lived peptide anion intermediate implies an insensitivity to the dynamics of conformational motion. Amides that are well-exposed in model conformational ensembles of ubiquitin vary a million-fold in exchange rates which continuum dielectric methods can predict with an rmsd of 3. However, the exchange rates for many of the more rarely exposed amides are markedly overestimated in the PDB-deposited 2K39 and 2KN5 ubiquitin ensembles, while the 2NR2 ensemble predictions are largely consistent with those of the Boltzmann-weighted conformational distribution sampled at the level of 1%. The correlation between the fraction of solvent-accessible conformations for a given amide hydrogen and the exchange rate constant for that residue provides a useful monitor of the degree of completeness with which a given ensemble has sampled the energetically accessible conformational space. These exchange predictions correlate with the degree to which each ensemble deviates from a set of 46 ubiquitin X-ray structures. Kolmogorov-Smirnov analysis for the distribution of intra- and inter-ensemble pairwise structural rmsd values assisted the identification of a subensemble of 2K39 that eliminates the overestimations of hydrogen exchange rates observed for the full ensemble. The relative merits of incorporating experimental restraints into the conformational sampling process are compared to using these restraints as filters to select subpopulations consistent with the experimental data.     

Simultaneous detection of intra- and inter-molecular paramagnetic relaxation enhancements in protein complexes

Paramagnetic relaxation enhancement (PRE) measurements constitute a powerful approach for detecting both permanent and transient protein–protein interactions. Typical PRE experiments require an intrinsic or engineered paramagnetic site on one of the two interacting partners; while a second, diamagnetic binding partner is labeled with stable isotopes (¹⁵N or ¹³C). Multiple paramagnetic labeled centers or reversed labeling schemes are often necessary to obtain sufficient distance restraints to model protein–protein complexes, making this approach time consuming and expensive. Here, we show a new strategy that combines a modified pulse sequence (¹H_N-Γ₂-CCLS) with an asymmetric labeling scheme to enable the detection of both intra- and inter-molecular PREs simultaneously using only one sample preparation. We applied this strategy to the non-covalent dimer of ubiquitin. Our method confirmed the previously identified binding interface for the transient di-ubiquitin complex, and at the same time, unveiled the internal structural dynamics rearrangements of ubiquitin upon interaction. In addition to reducing the cost of sample preparation and speed up PRE measurements, by detecting the intra-molecular PRE this new strategy will make it possible to measure and calibrate inter-molecular distances more accurately for both symmetric and asymmetric protein–protein complexes.     

Solution structure of human calcitonin gene-related peptide by 1H NMR and distance geometry with restrained molecular dynamics

The structure of human calcitonin gene-related peptide 1 (hCGRP-1) has been determined by 1H NMR in a mixed-solvent system of 50% trifluoroethanol/50% H2O at pH 3.7 and 27 degrees C. Complete resonance assignment was achieved by using two-dimensional methods. Distance restraints for structure calculations were obtained by semiquantitative analysis of intra- and interresidue nuclear Overhauser effects; in addition, stereospecific or X1 rotamer assignments were obtained for certain side chains. Structures were generated from the distance restraints by distance geometry, followed by refinement using molecular dynamics, and were compared with experimental NH-C alpha H coupling constants and amide hydrogen exchange data. The structure of hCGRP-1 in this solvent comprises an amino-terminal disulfide-bonded loop (residues 2-7) leading into a well-defined alpha-helix between residues 8 and 18; thereafter, the structure is predominantly disordered, although there are indications of a preference for a turn-type conformation between residues 19 and 21. Comparison of spectra for the homologous hCGRP-2 with those of hCGRP-1 indicates that the conformations of these two forms are essentially identical.     

Conformational energy calculations on glycosylated turns in glycoproteins

Analysis of the amino acid sequence of glycoproteins has suggested the β-turn as a likely site of glycosylation in glycoproteins. According to this model, the peptide chain traverses the interior of a globular protein, reversing its direction at the protein surface, a likely point for the attachment of hydrophilic carbohydrate residues. In order to search for plausible conformations of glycosylated β-turns in asparagine-linked glycoproteins, we have adapted the conformational energy calculation method of Scheraga and coworkers for use in carbohydrates. The parameters for nonbonded and hydrogen-bonded interactions have been published, and electrostatic parameters are derived from a CNDO calculation on a model glycopeptide. Our results indicate that the orientation of the glycosyl amide bond having the amide proton nearly trans to the anomeric proton of the sugar has the lowest energy. Although CD and nmr experiments in our laboratory have consistently found this conformation, our calculations show the conformation having these two protons in a cis relationship to lie very close in energy. Calculations on the glycopeptide linkage model, α-N-acetyl, δ-N(2-acetamido-1,2-dideoxy-β-D -glucopyranosyl)-N′-methyl-L -asparaginyl amide show that several distinct geometries are allowed for glycosylated β-turns. For a type I β-turn, three conformations of the glycosylated side chain are found within 4 kcal of the minimum, while two conformations of the glycosylated side chain are allowed for a type II turn. The hydrogen-bonded C7 conformation is also allowed. Stereoviews of the low-energy conformations reveal no major hydrogen-bonding interaction between the peptide and sugar.     

Conformational interconversions in peptide β-turns: Discrimination between enantiomeric conformations by chiral perturbation

The crystal structure of the model tripeptide Boc-Aib-Gly-Leu-OMe ( 1 ) reveals two independent molecules in the asymmetric unit that adopt “enantiomeric” type I and type I′ β-turn conformations with the Aib and Gly residues occupying the corner (i + 1 and i + 2) positions. ¹³C cross polarization and magic angle sample spinning spectra in the solid state also support the coexistence of two conformational species. ¹³C-nmr in CDCl₃ establishes the presence of a single species or rapid exchange between conformations. 400 MHz ¹H-nmr provides evidence for conformational exchange involving a major and minor species, with β-turn conformations supported by the low solvent exposure of Leu(3) NH and the observation of N_iH ↔ N_i+1H nuclear Overhauser effects. CD bands in the region 190–230 nm are positive, supporting a major population of type I′ β-turns. The isomeric peptide, Boc-Gly-Leu-Aib-OMe ( 2 ), adopts an “open” type II′ β-turn conformation in crystals. Solid state and solution nmr support population of a single conformational species. Chiral perturbation introduced outside the folded region of peptides may provide a means of modulating screw sense in achiral sequences. © 1998 John Wiley & Sons, Inc. Biopoly 45: 191–202, 1998     

Energetically unfavorable amide conformations for N6‐acetyllysine side chains in refined protein structures

The reversible acetylation of lysine to form N6‐acetyllysine in the regulation of protein function is a hallmark of epigenetics. Acetylation of the positively charged amino group of the lysine side chain generates a neutral N‐alkylacetamide moiety that serves as a molecular “switch” for the modulation of protein function and protein–protein interactions. We now report the analysis of 381 N6‐acetyllysine side chain amide conformations as found in 79 protein crystal structures and 11 protein NMR structures deposited in the Protein Data Bank (PDB) of the Research Collaboratory for Structural Bioinformatics. We find that only 74.3% of N6‐acetyllysine residues in protein crystal structures and 46.5% in protein NMR structures contain amide groups with energetically preferred trans or generously trans conformations. Surprisingly, 17.6% of N6‐acetyllysine residues in protein crystal structures and 5.3% in protein NMR structures contain amide groups with energetically unfavorable cis or generously cis conformations. Even more surprisingly, 8.1% of N6‐acetyllysine residues in protein crystal structures and 48.2% in NMR structures contain amide groups with energetically prohibitive twisted conformations that approach the transition state structure for cis‐trans isomerization. In contrast, 109 unique N‐alkylacetamide groups contained in 84 highly accurate small molecule crystal structures retrieved from the Cambridge Structural Database exclusively adopt energetically preferred trans conformations. Therefore, we conclude that cis and twisted N6‐acetyllysine amides in protein structures deposited in the PDB are erroneously modeled due to their energetically unfavorable or prohibitive conformations. Proteins 2013; © 2012 Wiley Periodicals, Inc.     

Omega amino acids in peptide design: incorporation into helices

Incorporation of easily available achiral ω-amino acid residues into an oligopeptide results in substitution of amide bonds by polymethylene units of an aliphatic chain, thereby providing a convenient strategy for constructing a peptidomimetic. The central Gly-Gly segment of the helical octapeptide Boc-Leu-Aib-Val-Gly-Gly-Leu-Aib-Val-Ome(1) has been replaced by δ-amino-valeric acid (δ-Ava) residue in the newly designed peptide Boc-Leu-Aib-Val-δ-Ava-Leu-Aib-Val-OMe(2). ¹H-nmr results clearly suggest that in the apolar solvent CDCl₃, the δ-Ava residue is accommodated into a folded helical conformation, stabilized by successive hydrogen bonds involving the NH groups of Val(3), δ-Ava(4), and Leu(5). The δ-Ava residue must adopt a gauche-gauche-trans-gauche-gauche conformation along the central polymethylene unit of the aliphatic segment, a feature seen in an energy-minimized model conformation based on nmr parameters. The absence of hydrogen bonding functionalities, however, limits the elongation of the helix. In fact, in CDCl₃, the folded conformation consists of an N-terminal helix spanning residues 1–4, followed by a Type II β-turn at residues 5 and 6, whereas in strongly solvating media like (CD₃)₂SO, the unfolding of the N-terminal helix results in β-turn conformations at Leu(1)-Aib(2). The Type II β-turn at the Leu(5)-Aib(6) segment remains intact even in (CD₃)₂SO. CD comparisons of peptides 1 and 2 reveal a “nonhelical” spectrum for 2 in 2,2,2-trifluoroethanol. © 1996 John Wiley & Sons, Inc.     

Laser Raman spectroscopic analysis of angiotensin peptides' conformation

In this study we examined the conformation and side chain environments of angiotensins I, II, III, and [Sar¹-Ile⁵-Ala⁸]angiotensin II using laser Raman spectroscopy. The positions of the amide I bands for all four peptides were found between 1664 and 1673 cm^?1. D₂O exchange studies confirmed the positions of the amide I and amide III bands. The positions of the amide I bands for all the angiotensins were found at approximately 1665 cm^?1 and the amide III bands were all located between 1265 and 1278 cm^?1. From the positions and intensities of the amide I and III bands we concluded that all peptides share the same overall conformation consisting of β-turn structure. Spectral analysis indicated that although the spectra for all the peptides were qualitatively identical there was evidence that the angiotensin conformations were more flexible in the aqueous phase than the solid phase. Examination of the $\text{850}\text{830}$ cm^?1 tyrosine doublet suggested that the tyrosine residue in the peptides is exposed to the solvent environment and becomes more exposed as the peptide length is decreased. Therefore, there are some localized conformational differences among the angiotensins. The conformational data yielded by this study leads us to conclude that the various biological properties ascribed to the angiotensins are not due to different conformations of the peptides. The biological differences could perhaps be attributed to localized interactions of the individual amino acid residues with themselves and with the hormone receptors.     

The amide III vibrational circular dichroism band as a probe to detect conformational preferences of alanine dipeptide in water

The conformational preferences of blocked alanine dipeptide (ADP), Ac‐Ala‐NHMe, in aqueous solution were studied using vibrational circular dichroism (VCD) together with density functional theory (DFT) calculations. DFT calculations of three most representative conformations of ADP surrounded by six explicit water molecules immersed in a dielectric continuum have proven high sensitivity of amide III VCD band shape that is characteristic for each conformation of the peptide backbone. The polyproline II (P_II) and α_R conformation of ADP are associated with a positive VCD band while β conformation has a negative VCD band in amide III region. Knowing this spectral characteristic of each conformation allows us to assign the experimental amide III VCD spectrum of ADP. Moreover, the amide III region of the VCD spectrum was used to determine the relative populations of conformations of ADP in water. Based on the interpretation of the amide III region of VCD spectrum we have shown that dominant conformation of ADP in water is P_II which is stabilized by hydrogen bonded water molecules between CO and NH groups on the peptide backbone. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 814–818, 2014.     

Raman and infrared spectroscopic study of gramicidin A conformations

Raman scattering and infrared spectroscopic techniques were used to study the vibrational spectrum and conformation of the membrane channel protein gramicidin A in the solid state, in organic solutions and, using Raman scattering only, in a phospholipid environment. The investigation also includes measurements on head- and tail-group-modifled gramicidin A and a potassium thiocyanate-gramicidin A complex. Tentative identification of the molecular vibrations is proposed on the basis of the data on model compounds. The existence of four distinct conformations of the gramicidin A chain is established: conformation I present in the solid state, and CH₃OH and CD₃OD solutions; conformation II present in films cast from CHCl₃ solution; conformation III present in (CH₃)₂SO and (CD₃)₂SO solutions at concentrations below 0.5 m gramicidin A; and conformation IV present in the potassium thiocyanate-gramicidin A complex. The data obtainable on a gramicidin A-phospholipid suspension indicate a gramicidin A conformation in this environment corresponding either to the conformation I or II. The details of the spectra in the amide I region are shown to be consistent with a β-parallel hydrogen-bonded π_LD helix for conformational I, in terms of the polypeptide vibrational calculations of Nevskaya and co-workers. Conformation II is found to be consistent with an antiparallel double-stranded π_LD helix, while conformations III and IV probably have π-helical structures with larger channel diameters. The data on head- and tail-modified gramicidin A molecules indicate that their conformations are only slightly different from that of gramicidin A in conformation I.     

β-turn tendency in N-methylated peptides with dehydrophenylalanine residue: DFT study

The tendency to adopt β‐turn conformation by model dipeptides with α,β‐dehydrophenylalanine (ΔPhe) residue in the gas phase and in solution is investigated by theoretical methods. We pay special attention to a dependence of conformational properties on the side‐chain configuration of dehydro residue and the influence of N‐methylation on β‐turn stability. An extensive computational study of the conformational preferences of Z and E isomers of dipeptides Ac‐Gly‐(E/Z)‐ΔPhe‐NHMe ( 1a / 1b ) and Ac‐Gly‐(E/Z)‐ΔPhe‐NMe₂ ( 2a / 2b ) by B3LYP/6‐311++G(d,p) and MP2/6‐311++G(d,p) methods is reported. It is shown that, in agreement with experimental data, Ac‐Gly‐(Z)‐ΔPhe‐NHMe has a great tendency to adopt β‐turn conformation. In the gas phase the type II β‐turn is preferred, whereas in the polar environment, the type I. On the other hand, dehydro residue in Ac‐Gly‐(E)‐ΔPhe‐NHMe has a preference to adopt extended conformations in all environments. N‐methylation of C‐terminal amide group, which prevents the formation of 1←4 intramolecular hydrogen bond, change dramatically the conformational properties of studied dehydropeptides. Especially, the tendency to adopt β‐turn conformations is much weaker for the N‐methylated Z isomer (Ac‐Gly‐(Z)‐ΔPhe‐NMe₂), both in vacuo and in the polar environment. On the contrary, N‐methylated E isomer (Ac‐Gly‐(E)‐ΔPhe‐NMe₂) can easier adopt β‐turn conformation, but the backbone torsion angles (?₁, ψ₁, ?₂, ψ₂) are off the limits for common β‐turn types. © 2012 Wiley Periodicals, Inc. Biopolymers 97:518–528, 2012.     

Raman study of crystalline peptides containing beta turns

The Raman spectra of crystalline H-ProLeuGlyNH₂ which has a type II β turn, crystalline S-benzylCysProLeuGlyNH₂ which has a type I β-turn, and crystalline gramicidin S which has two β turns and β-sheet structure in its conformation, were investigated. The amide I and amide III bands of the peptides with β turns were generally different from those which are diagnostic for α-helix and β-sheet conformations. The patterns of the amide I and amide III bands, when examined together, indicate that Raman spectra can provide diagnostic evidence for β-turn structure in peptides.