The conformational properties of α,β‐dehydroamino acids with a C‐terminal ester group

α,β‐Dehydroamino acid esters occur in nature. To investigate their conformational properties, a systematic theoretical analysis was performed on the model molecules Ac‐ΔXaa‐OMe [ΔXaa = ΔAla, (E)‐ΔAbu, (Z)‐ΔAbu, ΔVal] at the B3LYP/6‐311+ + G(d,p) level in the gas phase as well as in chloroform and water solutions with the self‐consistent reaction field‐polarisable continuum model method. The Fourier transform IR spectra in CCl₄ and CHCl₃ have been analysed as well as the analogous solid state conformations drawn from The Cambridge Structural Database. The ΔAla residue has a considerable tendency to adopt planar conformations C5 (?, ψ ≈ ? 180°, 180°) and β2 (?, ψ ≈ ? 180°, 0°), regardless of the environment. The ΔVal residue prefers the conformation β2 (?, ψ ≈ ? 120°, 0°) in a low polar environment, but the conformations α (?, ψ ≈ ? 55°, 35°) and β (?, ψ ≈ ? 55°, 145°) when the polarity increases. The ΔAbu residues reveal intermediate properties, but their conformational dispositions depend on configuration of the side chain of residue: (E)‐ΔAbu is similar to ΔAla, whereas (Z)‐ΔAbu to ΔVal. Results indicate that the low‐energy conformation β2 is the characteristic feature of dehydroamino acid esters. The studied molecules constitute conformational patterns for dehydroamino acid esters with various side chain substituents in either or both Z and E positions. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Conformation analysis of the N-glycosylation site Asn-X-Thr/Ser in glycoproteins

Theoretical conformational analysis of oligopeptides CH3CO-Asn-X-Thr-NHCH3 (X = Gly, Ala, Pro), modelling N-glycosylation site, and their glycosylated derivatives CH3CO-(GlcNAc beta 1-4GlcNAc beta 1) Asn-X-Thr-NHCH3 has been carried out. Active conformations of the site are found, corresponding to structural prerequisities of N-glycosylation: Asn residue's position in beta-turn and hydrogen bond formation between side chains of Asn and Thr/Ser residues. In this case the L conformation of the central residue X is most probable. Since Pro residue does not possess this conformation, sequences with X = Pro are not glycosylated. It is shown that glycosylation of the above-mentioned sites is accompanied by reorientation of the Asn residue's side chains.     

Sparse networks of directly coupled,polymorphic, and functional side chains in allosteric proteins

Recent studies have highlighted the role of coupled side‐chain fluctuations alone in the allosteric behavior of proteins. Moreover, examination of X‐ray crystallography data has recently revealed new information about the prevalence of alternate side‐chain conformations (conformational polymorphism), and attempts have been made to uncover the hidden alternate conformations from X‐ray data. Hence, new computational approaches are required that consider the polymorphic nature of the side chains, and incorporate the effects of this phenomenon in the study of information transmission and functional interactions of residues in a molecule. These studies can provide a more accurate understanding of the allosteric behavior. In this article, we first present a novel approach to generate an ensemble of conformations and an efficient computational method to extract direct couplings of side chains in allosteric proteins, and provide sparse network representations of the couplings. We take the side‐chain conformational polymorphism into account, and show that by studying the intrinsic dynamics of an inactive structure, we are able to construct a network of functionally crucial residues. Second, we show that the proposed method is capable of providing a magnified view of the coupled and conformationally polymorphic residues. This model reveals couplings between the alternate conformations of a coupled residue pair. To the best of our knowledge, this is the first computational method for extracting networks of side chains' alternate conformations. Such networks help in providing a detailed image of side‐chain dynamics in functionally important and conformationally polymorphic sites, such as binding and/or allosteric sites. Proteins 2015; 83:497–516. © 2014 Wiley Periodicals, Inc.     

Conformation of the central sequence of angiotensin II and analogs

We describe the synthesis and the conformational analysis by ir, CD, and proton-nmr spectroscopy of four model peptides of the type N-Ac-Tyr-X-His-NH₂ with X = Val, Leu, Ala, Gly. These peptides represent the central sequence of the hormone angiotensin II and its position-5 analogs. We studied their conformational behavior in aqueous solution during pH titration and in organic solvents. For specific purposes of spectral analysis (ir band assignment, proton-nmr signal assignment, heteronuclear vicinal coupling constants), we synthesized three isotopically enriched homologs of the mother sequence, i.e., N-Ac-(¹⁵N-Tyr)-Val-His-NH₂, N-Ac-(¹³C, ²H, Tyr)-Val-His-NH₂, and N-Ac-Tyr-(¹³C, ²H, Val)-His-NH₂. Results are summarized as follows: the tyrosine and the histidine side chains influence each other through space; this mutual influence is modulated by the nature of the side chain in position X and decreases in going from X?Val to X?Gly as a consequence of two simultaneous events, changes in the side-chain rotamer distribution and changes in the φ and ψ angles of residue X. The decrease in the bulkiness of the side-chain X (Val → Gly) leads to increased flexibility of the peptide backbone at this site, which is also reflected in the apparent ratio of C₅, C₇, and intermediate conformations present in equilibrium. The three spectroscopic techniques, in addition to the results of chymotryptic degradation experiments, show a high level of agreement, and all reflect the dynamic conformation of these peptides in a different manner.     

Asparagine and glutamine differ in their propensities to form specific side chain-backbone hydrogen bonded motifs in proteins

Short range side chain‐backbone hydrogen bonded motifs involving Asn and Gln residues have been identified from a data set of 1370 protein crystal structures (resolution ≤ 1.5 Å). Hydrogen bonds involving residues i ? 5 to i + 5 have been considered. Out of 12,901 Asn residues, 3403 residues (26.4%) participate in such interactions, while out of 10,934 Gln residues, 1780 Gln residues (16.3%) are involved in these motifs. Hydrogen bonded ring sizes (C_n, where n is the number of atoms involved), directionality and internal torsion angles are used to classify motifs. The occurrence of the various motifs in the contexts of protein structure is illustrated. Distinct differences are established between the nature of motifs formed by Asn and Gln residues. For Asn, the most highly populated motifs are the C₁₀ (CO^δ_i …NH_{i + 2}), C₁₃ (CO^δ_i …NH_{i + 3}) and C₁₇ (N^δH_i …CO_{i ? 4}) structures. In contrast, Gln predominantly forms C₁₆ (CO^ε_i …NH_{i ? 3}), C₁₂ (N^εH_i …CO_{i ? 2}), C₁₅ (N^εH_i …CO_{i ? 3}) and C₁₈ (N^εH_i …CO_{i ? 4}) motifs, with only the C₁₈motif being analogous to the Asn C₁₇structure. Specific conformational types are established for the Asn containing motifs, which mimic backbone β‐turns and α‐turns. Histidine residues are shown to serve as a mimic for Asn residues in side chain‐backbone hydrogen bonded ring motifs. Illustrative examples from protein structures are considered. Proteins 2012; © 2011 Wiley Periodicals, Inc.     

Stabilization of collagen fibrils by hydroxyproline

The substitution of hydroxyproline for proline in position Y of the repeating Gly-X-Y tripeptide sequence of collagen-like poly(tripeptide)s (i.e., in the position in which Hyp occurs naturally) is predicted to enhance the stability of aggregates of triple helices, while the substitution of Hyp in position X (where no Hyp occurs naturally) is predicted to decrease the stability of aggregates. Earlier conformational energy computations have indicated that two triple helices composed of poly(Gly-Pro-Pro) polypeptide chains pack preferentially with a nearly parallel orientation of the helix axes [Nemethy, G., & Scheraga, H.A. (1984) Biopolymers 23, 2781-2799]. Conformational energy computations reported here indicate that the same packing arrangement is preferred for the packing of two poly(Gly-Pro-Hyp) triple helices. The OH groups of the Hyp residues can be accommodated in the space between the two packed triple helices without any steric hindrance. They actually contribute about 1.9 kcal/mol per Gly-Pro-Hyp tripeptide to the packing energy, as a result of the formation of weak hydrogen bonds and other favorable noncovalent interatomic interactions. On the other hand, the substitution of Hyp in position X weakens the packing by about 1.7 kcal/mol per Gly-Hyp-Pro tripeptide. Numerous published experimental studies have established that Hyp in position Y stabilizes an isolated triple helix relative to dissociated random coils, while Hyp in position X has the opposite effect. We propose that Hyp in position Y also enhances the stability of the assembly of collagen into microfibrils while, in position X, it decreases this stability.     

Side chain interactions in aromatic dipeptides

A study of the 100-MHZ nuclear magnetic resonace spectra in D₂O solution was made of a series of linear dipeptides of the types L -phenylalanine-L -and-D -X, and L -phenylalanine-L -and-D -Y, where X comprised a group of amino acid residues with polar side chains (X = glutamine, glutamic acid, arginine, and N^ε-acetyllysine) and Y comprised amino acid residues with purely aliphatic side chains (Y = α-aminobutyric acid and norvaline). It was found that regardless of the side chain length, resonances due to the α-methylene protons in the X and Y side chains of the L -Phe-D -Y series consistently exhibited upfield shifts greater than any other protons in these side chains, when compared to the corresponding side chain resonances of the nonaromatic dipeptide series L -Ala-L -X and L -Ala-L -Y. The magnitudes of these shielding effects were consistently and considerably greater for the L -Phe-D -X series than for the L -Phe-D -Y series. An intramolecular complex–formed by association of armatic π-electrons with the positive end of the dipole in the polar side chains—was proposed as one plausible interpretation of the enhanced shielding effects. An increase in temperature from 32 to 70–80° was sufficient to overcome the enhanced shielding attributable to the suggested complex.     

Monte Carlo study of the effect of β2-microglobulin on the binding cleft of the HLA-A2 complex

Peptide recognition by class I products of the major histocompatibility complex requires association of the class I heavy chain with β₂-microglobulin. We present results of Monte Carlo simulations of the β-pleated sheet floor of the human class I MHC molecule, HLA-A2, with and without β₂-microglobulin. We find a significant effect of β₂-microglobulin on the side chains of residues near a region that would accommodate the C-terminus of a bound peptide. By modeling simultaneously each loop and its neighboring strand at either end of the class I cleft, we find that β₂-microglobulin restricts the conformational space of residues that are central to binding peptides. The effect is most pronounced for R97 and H114 and somewhat less important for Y99 and Y116, the latter forming strong hydrogen bonds with neighboring residues in the heavy chain itself.     

Conformational constraints of amino acid side chains in alpha-helices

The conformational freedom of amino acid side chains is strongly reduced when the side chains occur on an α-helix. A quantitative evaluation of this freedom has been carried out by means of conformational energy computations for all naturally occurring amino acids and for α-aminobutyric acid when they are placed in the middle of a right-handed poly(L-alanine) α-helix. One of the three possible rotameric states for rotation around the C^α ? C^β bond (viz. g⁺) is excluded completely on the helix because of steric hindrance, and the relative populations of the other two rotamers (t and g^?) are altered because of steric interactions and the reduction of hydrogen-bonding possibilities. The computed tendencies of the changes in distributions of rotamers, on going from an ensemble of all backbone conformations to the α-helix, agree with the observed tendencies in proteins. Minimum-energy side-chain conformations in an α-helix have been tabulated for use in conformational energy computations on polypeptides.     

Conformational characteristics of polyglutamic acid esters: Persistence of orientational correlation along the side chain flanking the α-helical backbone

The side chain conformations of α-helical poly(L -glutamic acid) esters $ \rlap{--}[NHCH(CH_2 CH_2 COOR)CO\rlap{--}]_x $, carrying a homologous series of ester residues such as R = ? (CH₂)_n? with n = 1–3, have been studied in the lyotropic liquid crystalline state (chloroform 20 v/v%) by the deuterium nmr method. In order to study the surface chirality of the molecule, the phenyl groups situated at the terminal of the side chain have been deuterated. From the observed deuterium quadrupolar splittings, the average inclination θ_p of the para-axis of the phenyl group with respect to the α-helical backbone was elucidated. A distinct odd–even oscillation in the quantity such as 〈 cos² θ_p〉 was observed with the number of methylene units n. A rotational isomeric state analysis has indicated that the observed orientational correlation arises from the interdependence of the neighboring bond rotation along the side chain. Preference of the “extended” conformations is also enhanced by the mutual conformational exclusion of neighboring side chains.     

The conformational properties of dehydrobutyrine and dehydrovaline: theoretical and solid‐state conformational studies

Dehydrobutyrine is the most naturally occurring dehydroamino acid. It is also the simplest dehydroamino acid having the geometrical isomers E/Z. To investigate its conformational properties, a theoretical analysis was performed on N‐acetyl‐α,β‐dehydrobutyrine N′‐methylamides, Ac‐(E)‐ΔAbu‐NHMe and Ac‐(Z)‐ΔAbu‐NHMe, as well as the dehydrovaline derivative Ac‐ΔVal‐NHMe. The ?, ψ potential energy surfaces and the localised conformers were calculated at the B3LYP/6‐311 + + G(d,p) level of theory both in vacuo and with inclusion of the solvent (chloroform, water) effect (SCRF method). The X‐ray crystal structures of Ac‐(Z)‐ΔAbu‐NHMe and Ac‐ΔVal‐NHMe were determined at 85 and 100 K, respectively. The solid‐state conformational preferences for the studied residues have been analysed and compared with the other related structures. Despite the limitations imposed by the C^α = C^β double bond on the topography of the side chains, the main chains of the studied dehydroamino acids are more flexible than in standard alanine. The studied dehydroamino acids differ in their conformational preferences, which depend on the polarity of the environment. This might be a reason why the nature quite precisely differentiates between ΔVal and each of the ΔAbu isomers, and why, particularly so with the latter, they are used as a conformational tool to influence the biological action of usually small, cyclic dehydropeptides. Copyright © 2010 European Peptide Society and John Wiley & Sons, Ltd.     

Influence of local interactions on protein structure. I. Conformational energy studies of N-acetyl-N′-methylamides of pro-X and X-pro dipeptides

Conformational energy calculations using an Empirical Conformational Energy Program for Peptides (ECEPP) were carried out on the N-acetyl-N′-methylamides of Pro-X, where X = Ala, Asn, Asp, Gly, Leu, Phe, Ser, and Val, and of X-Pro, where X = Ala, Asn, Gly, and Pro. The conformational energy was minimized from starting conformations which included all combinations of low-energy single-residue minima and several standard bend structures. It was found that almost all resulting minima are combinations of low-energy single-residue minima, suggesting that intra residue interactions predominate in determining conformation. The calculations also indicate, however, that inter residue interactions can be important. In addition, librational entropy was found to influence the relative stabilities of some minima. Because of the existence of 10–100 low-energy minima for each dipeptide, the normalized statistical weight of an individual minimum rarely exceeds 0.3, suggesting that these dipeptides have considerable conformational flexibility and exist as statistical ensembles of low-energy structures. The propensity of each dipeptide to form bend conformations was calculated, and the results were compared with available experimental data. It was found that bends are favored in Pro-X dipeptides because ?_Pro is fixed by the pyrrolidine ring in a conformation which is frequently found in bends, but that bends are not favored in X-Pro dipeptides because interactions between the X residue and the pyrrolidine ring restrict the X residue to conformations which are not usually found in bends.     

Structural properties of amyloid β(1‐40) dimer explored by replica exchange molecular dynamics simulations

Replica exchange molecular dynamics simulations (300 ns) were used to study the dimerization of amyloid β(1‐40) (Aβ(1‐40)) polypeptide. Configurational entropy calculations revealed that at physiological temperature (310 K, 37°C) dynamic dimers are formed by randomly docked monomers. Free energy of binding of the two chains to each other was ?93.56 ± 6.341 kJ mol^?1. Prevalence of random coil conformations was found for both chains with the exceptions of increased β‐sheet content from residues 16‐21 and 29‐32 of chain A and residues 15‐21 and 30‐33 of chain B with β‐turn/β‐bend conformations in both chains from residues 1‐16, 21‐29 of chain A, 1‐16, and 21‐29 of chain B. There is a mixed β‐turn/β‐sheet region from residues 33‐38 of both chains. Analysis of intra‐ and interchain residue distances shows that, although the individual chains are highly flexible, the dimer system stays in a loosely packed antiparallel β‐sheet configuration with contacts between residues 17‐21 of chain A with residues 17‐21 and 31‐36 of chain B as well as residues 31‐36 of chain A with residues 17‐21 and 31‐36 of chain B. Based on dihedral principal component analysis, the antiparallel β‐sheet‐loop‐β‐sheet conformational motif is favored for many low energy sampled conformations. Our results show that Aβ(1‐40) can form dynamic dimers in aqueous solution that have significant conformational flexibility and are stabilized by collapse of the central and C‐terminal hydrophobic cores with the expected β‐sheet‐loop‐β‐sheet conformational motif. Proteins 2017; 85:1024–1045. © 2017 Wiley Periodicals, Inc.     

Disorder and order in unfolded and disordered peptides and proteins: A view derived from tripeptide conformational analysis. II. Tripeptides with short side chains populating asx and β‐type like turn conformations

In the preceding paper, we found that ensembles of tripeptides with long or bulky chains can include up to 20% of various turns. Here, we determine the structural and thermodynamic characteristics of GxG peptides with short polar and/or ionizable central residues (D, N, C), whose conformational distributions exhibit higher than average percentage (>20%) of turn conformations. To probe the side‐chain conformations of these peptides, we determined the ³J(H^α,H^β) coupling constants and derived the population of three rotamers with χ₁‐angles of ?60°, 180° and 60°, which were correlated with residue propensities by DFT‐calculations. For protonated GDG, the rotamer distribution provides additional evidence for asx‐turns. A comparison of vibrational spectra and NMR coupling constants of protonated GDG, ionized GDG, and the protonated aspartic acid dipeptide revealed that side chain protonation increases the pPII content at the expense of turn populations. The charged terminal groups, however, have negligible influence on the conformational properties of the central residue. Like protonated GDG, cationic GCG samples asx‐turns to a significant extent. The temperature dependence of the UVCD spectra and ³J(H^NH^α) constants suggest that the turn populations of GDG and GNG are practically temperature‐independent, indicating enthalpic and entropic stabilization. The temperature‐independent J‐coupling and UVCD spectra of GNG require a three‐state model. Our results indicate that short side chains with hydrogen bonding capability in GxG segments of proteins may serve as hinge regions for establishing compact structures of unfolded proteins and peptides. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

The effects of isolated N-methylated residues on the conformational characteristics of polypeptides

Conformational energies have been estimated for the tripeptide fragments L -Ala-N-methyl-L -Ala-L -Ala, L -Ala-L -Ala-N-methyl-L -Ala, L -Ala-Sar-L -Ala, and L -Ala-Gly-N-methyl-L -Ala. The peptide bonds connecting L -Ala and Gly with N-methyl-L -Ala and L -Ala with Sar were permitted to adopt the planar cis as well as the usual trans conformation. Contour maps of the conformational energies of the central residue in these tripeptide fragments are presented and compared to the conformational energy maps previously calculated for unmethylated L -Ala and Gly surrounded by residues which are also unmethylated. In generl it is observed that L -Ala and Gly residues that are either N-methylated in their conformational freedom relative to the same residues in an unmethylated polypeptide chain.     

Exploring the conformational space of protein side chains using dead-end elimination and the A* algorithm

We describe an algorithm which enables us to search the conformational space of the side chains of a protein to identify the global minimum energy combination of side chain conformations as well as all other conformations within a specified energy cutoff of the global energy minimum. The program is used to explore the side chain conformational energy surface of a number of proteins, to investigate how this surface varies with the energy model used to describe the interactions within the system and the rotamer library. Enumeration of the rotamer combinations enables us to directly evaluate the partition function, and thus calculate the side chain contribution to the conformational entropy of the folded protein. An investigation of these conformations and the relationships between them shows that most of the conformations near to the global energy minimum arise from changes in side chain conformations that are essentially independent; very few result from a concerted change in conformation of two or more residues. Some of the limitations of the approach are discussed. Proteins 33:227–239, 1998. © 1998 Wiley-Liss, Inc.     

Atypical gly-X-Y sequences surround interruptions in the repeating tripeptide pattern of basement membrane collagen

The triple-helical domains of type IV collagen chains have more than 20 sites at which the repeating (Gly-X-Y)_n pattern is interrupted. Analysis of α1 (IV) and α2 (IV) chains indicates the residues in the three Gly-X-Y triplets preceding or following interruptions differ statistically from the rest of the chain. Unusually high frequencies of charged residues are seen at a number of X and Y sites, with the charge density being particularly high C-terminal to the interruption site. Analyses were carried out on individual categories of interruptions, classified as insertions or deletions in the Y position. All of the residues in the X and Y positions of the triplets flanking insertion sites are atypical, with a high concentration of charged residues. Triplets flanking sites where there has been a deletion in the Y position show unusually high frequencies of charged residues at some sites, hydrophobic residues at other sites, and an invariant imino acid N-terminal to the interruption. The presence of atypical sequences surrounding interruptions could be important at a molecular level, related to triple-helix stability, or at a supramolecular level, related to the association of molecules to form networks in basement membranes. © 1995 John Wiley & Sons, Inc.     

A Hydrogen Bonded Chain in Bactereorhodopsin by Computer Modelling Approach

Abstract

The seven α-helical segments of Bacteriorhodopsin (BR) passing through the membrane are investigated for a continuous Hydrogen Bonded Chain (HBC). The study is carried out by computer modelling approach. It is assumed that the seven helices are placed as (AGFEDCB), which has been accepted as the best model by several groups. Helices A, D, E and G are considered to be present in right handed α-helical conformation. The inter-orientation of these helices are represented by Eulerian angles α, β and γ. For the helices B, C and F which contain Proline in the middle, several conformational possibilities were considered. In these cases apart from the Eulerian angles α, β and γ, the dihedral angles φ_{p_1} and ψ_{p_1} of the residues that are succeeded by Proline residue in the helical regions were also used in fixing the position of the helices with respect to each other. All these parameters were varied to fit with the top, middle and bottom distances reported by electron diffraction studies. Good fit was obtained for all right handed α-helical conformations and also for helices B, C and F with a left handed turn at the residue preceeding proline. Hence two structures were analysed for continuous HBC, Structure I which contained all the seven helices in right handed α-helical conformation and Structure II, which had the helices A, D, E and G in right handed conformation and the helices B, C and F in right handed α-helical conformation with a left handed turn at the residue preceeding proline. All possible staggered conformations were considered for the side chains and the inter atomic distances were analysed for Hydrogen bonds. It was possible to obtain a continuous chain in both the structures which includes most of the residues found to be important by the experiments. However Lys-216 has to be considered in two different conformations to connect the cytoplasmic side with the extra cellular side. The overall height spanned by HBC is about 25Å. The chains obtained by both the structures I and II are analysed in terms of the conformational parameters. It has also been possible to place the retinal in the region as predicted by the experiments. The Tryptophan residues which affect the spectral characterestics can be aligned on either side of the retinal.     

Flexibility of aromatic residues in the active-site gorge of acetylcholinesterase: X-ray versus molecular dynamics

The high aromatic content of the deep and narrow active-site gorge of acetylcholinesterase (AChE) is a remarkable feature of this enzyme. Here, we analyze conformational flexibility of the side chains of the 14 conserved aromatic residues in the active-site gorge of Torpedo californica AChE based on the 47 three-dimensional crystal structures available for the native enzyme, and for its complexes and conjugates, and on a 20-ns molecular dynamics (MD) trajectory of the native enzyme. The degree of flexibility of these 14 aromatic side chains is diverse. Although the side-chain conformations of F330 and W279 are both very flexible, the side-chain conformations of F120, W233, W432, Y70, Y121, F288, F290 and F331 appear to be fixed. Residues located on, or adjacent to, the Ω-loop (C67-C94), namely W84, Y130, Y442, and Y334, display different flexibilities in the MD simulations and in the crystal structures. An important outcome of our study is that the majority of the side-chain conformations observed in the 47 Torpedo californica AChE crystal structures are faithfully reproduced by the MD simulation on the native enzyme. Thus, the protein can assume these conformations even in the absence of the ligand that permitted their experimental detection. These observations are pertinent to structure-based drug design.