Structural compromise of disallowed conformations in peptide and protein structures

Using a data set of 454 crystal structures of peptides and 80 crystal structures of non-homologous proteins solved at ultra high resolution of 1.2 A or better we have analyzed the occurrence of disallowed Ramachandran (phi, psi) angles. Out of 1492 and 13508 non-glycyl residues in peptides and proteins respectively 12 and 76 residues in the two datasets adopt clearly disallowed combinations of Ramachandran angles. These examples include a number of conformational points which are far away from any of the allowed regions in the Ramachandran map. According to the Ramachandran map a given (phi, psi) combination is considered disallowed when two non-bonded atoms in a system of two-linked peptide units with ideal geometry are prohibitively proximal in space. However, analysis of the disallowed conformations in peptide and protein structures reveals that none of the observations of disallowed conformations in the crystal structures correspond to a short contact between non-bonded atoms. A further analysis of deviations of bond lengths and angles, from the ideal peptide geometry, at the residue positions of disallowed conformations in the crystal structures suggest that individual bond lengths and angles are all within acceptable limits. Thus, it appears that the rare tolerance of disallowed conformations is possible by gentle and acceptable deviations in a number of bond lengths and angles, from ideal geometry, over a series of bonds resulting in a net gross effect of acceptable non-bonded inter-atomic distances.     

Sparsely populated residue conformations in protein structures: Revisiting “experimental” Ramachandran maps

The Ramachandran map clearly delineates the regions of accessible conformational (φ–ψ) space for amino acid residues in proteins. Experimental distributions of φ, ψ values in high‐resolution protein structures, reveal sparsely populated zones within fully allowed regions and distinct clusters in apparently disallowed regions. Conformational space has been divided into 14 distinct bins. Residues adopting these relatively rare conformations are presented and amino acid propensities for these regions are estimated. Inspection of specific examples in a completely “arid”, fully allowed region in the top left quadrant establishes that side‐chain and backbone interactions may provide the energetic compensation necessary for populating this region of φ–ψ space. Asn, Asp, and His residues showed the highest propensities in this region. The two distinct clusters in the bottom right quadrant which are formally disallowed on strict steric considerations correspond to the gamma turn (C7 axial) conformation (Bin 12 ) and the i + 1 position of Type II′ β turns (Bin 13) . Of the 516 non‐Gly residues in Bin 13 , 384 occupied the i + 1 position of Type II′ β turns. Further examination of these turn segments revealed a high propensity to occur at the N‐terminus of helices and as a tight turn in β hairpins. The β strand–helix motif with the Type II′ β turn as a connecting element was also found in as many as 57 examples. Proteins 2014; 82:1101–1112. © 2013 Wiley Periodicals, Inc.     

Descriptive statistics of disallowed regions and various protein secondary structures in the context of studying twisted β-hairpins

A detailed analysis of polypeptide-chain backbone conformations was carried out for polypeptide-chain segments adjacent to β-turn regions, including the sites of disallowed conformations. A cross comparison of conformations was performed for disallowed regions of the Ramachandran plot and main types of β-turns and adjacent secondary structures. Based on the results, disallowed region 2 (II, II') in the Ramachandran plot was shown to coincide mainly with β-hairpins and, more exactly, twisted β-hairpins. The frequency of residues with angles ?_i, ψ_i that fall in region 2 (II, II') in the latter is 140 times higher than in common β-hairpins.     

Structural assessment of glycyl mutations in invariantly conserved motifs

Motifs that are evolutionarily conserved in proteins are crucial to their structure and function. In one of our earlier studies, we demonstrated that the conserved motifs occurring invariantly across several organisms could act as structural determinants of the proteins. We observed the abundance of glycyl residues in these invariantly conserved motifs. The role of glycyl residues in highly conserved motifs has not been studied extensively. Thus, it would be interesting to examine the structural perturbations induced by mutation in these conserved glycyl sites. In this work, we selected a representative set of invariant signature (IS) peptides for which both the PDB structure and mutation information was available. We thoroughly analyzed the conformational features of the glycyl sites and their local interactions with the surrounding residues. Using Ramachandran angles, we showed that the glycyl residues occurring in these IS peptides, which have undergone mutation, occurred more often in the L-disallowed as compared with the L-allowed region of the Ramachandran plot. Short range contacts around the mutation site were analyzed to study the steric effects. With the results obtained from our analysis, we hypothesize that any change of activity arising because of such mutations must be attributed to the long-range interaction(s) of the new residue if the glycyl residue in the IS peptide occurred in the L-allowed region of the Ramachandran plot. However, the mutation of those conserved glycyl residues that occurred in the L-disallowed region of the Ramachandran plot might lead to an altered activity of the protein as a result of an altered conformation of the backbone in the immediate vicinity of the glycyl residue, in addition to long range effects arising from the long side chains of the new residue. Thus, the loss of activity because of mutation in the conserved glycyl site might either relate to long range interactions or to local perturbations around the site depending upon the conformational preference of the glycyl residue.     

Beyond basins: φ,ψ preferences of a residue depend heavily on the φ,ψ values of its neighbors

The Ramachandran plot distributions of nonglycine residues from experimentally determined structures are routinely described as grouping into one of six major basins: β, P_II, α, α_L, ξ and γ'. Recent work describing the most common conformations adopted by pairs of residues in folded proteins [i.e., (φ,ψ)₂‐motifs] showed that commonly described major basins are not true single thermodynamic basins, but are composed of distinct subregions that are associated with various conformations of either the preceding or following neighbor residue. Here, as documentation of the extent to which the conformational preferences of a central residue are influenced by the conformations of its two neighbors, we present a set of φ,ψ‐plots that are delimited simultaneously by the φ,ψ‐angles of its neighboring residues on both sides. The level of influence seen here is typically greater than the influence associated with considering the identities of neighboring residues, implying that the use of this heretofore untapped information can improve the accuracy of structure prediction algorithms and low resolution protein structure refinement.     

Crystallographic characterization of the α,γ C12 helix in hybrid peptide sequences

The solid‐state conformations of two αγ hybrid peptides Boc‐[Aib‐γ⁴(R)Ile]₄‐OMe 1 and Boc‐[Aib‐γ⁴(R)Ile]₅‐OMe 2 are described. Peptides 1 and 2 adopt C₁₂‐helical conformations in crystals. The structure of octapeptide 1 is stabilized by six intramolecular 4 → 1 hydrogen bonds, forming 12 atom C₁₂ motifs. The structure of peptide 2 reveals the formation of eight successive C₁₂ hydrogen‐bonded turns. Average backbone dihedral angles for αγ C₁₂ helices are peptide 1 , Aib; φ (°) = ?57.2 ± 0.8, ψ (°) = ?44.5 ± 4.7; γ⁴(R)Ile; φ (°) = ?127.3 ± 7.3, θ₁ (°) = 58.5 ± 12.1, θ₂ (°) = 67.6 ± 10.1, ψ (°) = ?126.2 ± 16.1; peptide 2 , Aib; φ (°) = ?58.8 ± 5.1, ψ (°) = ?40.3 ± 5.5; ψ⁴(R)Ile; φ (°) = ?123.9 ± 2.7, θ₁ (°) = 53.3 θ 4.9, θ ₂ (°) = 61.2 ± 1.6, ψ (°) = ?121.8 ± 5.1. The tendency of γ⁴‐substituted residues to adopt gauche–gauche conformations about the C^α–C^β and C^β–C^γ bonds facilitates helical folding. The αγ C₁₂ helix is a backbone expanded analog of α peptide 3₁₀ helix. The hydrogen bond parameters for α peptide 3₁₀ and α‐helices are compared with those for αγ hybrid C₁₂ helix. Copyright © 2016 European Peptide Society and John Wiley & Sons.     

On residues in the disallowed region of the Ramachandran map.

An analysis of the occurrence of nonglycyl residues in conformations disallowed in the Ramachandran plot is presented. Ser, Asn, Thr, and Cys have the highest propensities to exhibit such conformations, and the branched aliphatic residues the lowest. Residues cluster in five regions and there are some trends in the types of residues and their side-chain conformations (chi(1)) occupying these. Majority of the residues are found at the edge of helices and strands and in short loops, and are involved in different types of weak, stabilizing interactions. A structural motif has been identified where a residue in disallowed conformation occurs as the first residue of a short 3(10)-helix. On the basis of the types of neighboring residues, the location in the three-dimensional structure and accessibility, there are similarities with the occurrence of cis peptide bonds in protein structures.     

Toward engineering efficient peptidomimetics. Screening conformational landscape of two modified dehydroaminoacids

Effective peptidomimetics should posses structural rigidity and appropriate interaction pattern leading to potential spatial and electronic matching to the target receptor site. Rational design of such small bioactive molecules could push chemical synthesis and molecular modeling toward faster progress in medicinal chemistry. Conformational properties of N‐t‐butoxycarbonyl‐glycine‐(E/Z)‐dehydrophenylalanine N′,N′‐dimethylamides (Boc‐Gly‐(E/Z)‐ΔPhe‐NMe₂) in chloroform were studied by NMR and IR spectroscopy. The experimental findings were supported by extensive calculations at DFT(B3LYP, M06‐2X) and MP2 levels of theory and the β‐turn tendency for both isomers of the studied dipeptide were determined in vacuum and in solution. The theoretical data and experimental IR results were used as an additional information for the NMR‐based determination of the detailed solution conformations of the peptides. The obtained results reveal that N‐methylation of C‐terminal amide group changes dramatically the conformational properties of studied dehydropeptides. Theoretical conformational analysis reveals that the tendency to adopt β‐turn conformations is much weaker for the N‐methylated Z isomer (Boc‐Gly‐(Z)‐ΔPhe‐NMe₂), both in vacuum and in polar environment. On the contrary, N‐methylated E isomer (Boc‐Gly‐(E)‐ΔPhe‐NMe₂) can easily adopt β‐turn conformation, but the backbone torsion angles (φ₁, ψ₁, φ₂, ψ₂) are off the limits for common β‐turn types. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 28–40, 2014.     

Influence of solvent molecules on the stereochemical code of glycyl residues in proteins

The Ramachandran steric map and energy diagrams of the glycyl residue are symmetric. A plot of (phi,psi) angles of glycyl residues in 250 nonhomologous and high-resolution protein structures is also largely symmetric. However, there is a clear aberration in the symmetry. Although there is a cluster of points corresponding to the right-handed alpha-helical region, the "equivalent" cluster is clearly shifted to in and around the (phi,psi) values of (90 degrees, 0 degrees ) instead of being centered at the left-handed alpha-helical region of (60 degrees, 40 degrees ). This lack of symmetry exists even in the (phi,psi) distribution of residues from non-alpha-helical regions in proteins. Here we provide an explanation for this observation. An analysis of glycyl conformations in small peptide structures and in "coil" proteins, which are largely devoid of helical and sheet regions, shows that glycyl residues prefer to adopt conformations around (+/-90 degrees, 0 degrees ) instead of right- and left-handed alpha-helical regions. By using theoretical calculations, such conformations are shown to have highest solvent accessibility in a system of two-linked peptide units with glycyl residue at the central C(alpha) atom. This finding is consistent with the observations from 250 nonhomologous protein structures where glycyl residues with conformations close to (+/-90 degrees, 0 degrees ) are seen to have high solvent accessibility. Analysis of a subset of nonhomologous structures with very high resolution (1.5 A or better) shows that water molecules are indeed present at distances suitable for hydrogen bond interaction with glycyl residues possessing conformations close to (+/-90 degrees, 0 degrees ). It is suggested that water molecules play a key role in determining and stabilizing these conformations of glycyl residues and explain the aberration in the symmetry of glycyl conformations in proteins.     

Conformational characteristics of polypeptides containing isolated L-proline residues with cis peptide bonds

The conformational characteristics of the peptide sequence X-l-Pro, where X  Gly or l-Ala and the peptide bond joining X and l-Pro is cis, are evaluated. Semi-empirical potential functions are used to estimate the contributions to the conformational energy made by the non-bonded van der Waals' and electrostatic interactions and the intrinsic torsional potentials about the NC^a and C^aC′ bonds. Rotations φ₁ and ψ₁ about the NC^a and C^aC′ bonds in residue X and rotation ψ₂ about the C^aC′ bond in l-Pro are permitted, while the angle of rotation φ₂ about the NC^a bond in l-Pro is fixed at 120 ° by the pyrrolidine ring. The presence of the cis peptide bond connecting X and l-Pro renders the backbone rotations φ₁, ψ₁ in X dependent upon the rotation ψ₂ about the C^aC′ bond in l-Pro. (Interdependence of rotations in neighboring residues joined by a cis peptide bond was previously observed in l-alanine oligomers.) The number of energetically allowed conformations for the Gly and l-Ala residues preceding a cis peptide bond l-Pro residue are found to be substantially reduced from those permitted when the peptide bond is trans or when l-Pro is replaced by an amino acid residue. On the other hand, ψ₂ = 100 to 160 ° (cis′) and 300 to 0 ° (trans′) are found to be the lowest energy conformations of the l-Pro residue irrespective of the cis or trans conformation of the X-l-Pro peptide bond.     

A high-resolution 13C-nmr study of collagenlike polypeptides and collagen fibrils in solid state studied by the cross-polarization–magic angle-spinning method. Manifestation of conformation-dependent 13C chemical shifts and application to conformational characterization

We have recorded high-resolution ¹³C-nmr spectra of collagen fibrils in the solid state by the cross-polarization–magic-angle-spinning(CP–MAS)method and analyzed the spectra with reference to those of collagenlike polypeptides. We used two kinds of model polypeptides to obtain reference ¹³C chemical shifts of major amino acid residues of collagen (Gly, Pro, Ala, and Hyp): the 3₁-helical polypeptides [(Gly)_nII, (Pro)_nII, (Hyp)_n, and (Ala? Gly? Gly)_nII], and the triple-helical polypeptides [(Pro? Gly? Pro)_n and (Pro? Ala? Gly)_n]. Examination of the ¹³C chemical shifts of these polypeptides, together with our previous data, showed that the ¹³C chemical shifts of individual amino acid residues are the same, within experimental error (±0.5 ppm), among different polypeptides with different primary sequences, if the conformations are the same. We found that the ¹³C chemical shifts of Ala residues of the 3₁-helical (Ala? Gly? Gly)_n and triple-helical (Pro? Ala? Gly)_n are significantly displaced, compared with those of the α-helix, β-sheet, and silk I form, and can be utilized as excellent probes to examine conformational features of collagen-like polypeptides. Further, the ¹³C chemical shifts of Gly and Pro residues in the triple-helical polypeptides are substantially displaced from those found in (Gly)_nII and (Pro)_nII of the 3₁-helix, reflecting further conformational change from the 3₁-helix to the supercoiled triple helix. In particular, the ¹³C chemical shifts of Gly C ? O carbons of the triple-helical polypeptides are substantially displaced upfield (4.1–5.1 ppm), with respect to those of the 3₁-helical polypeptides. These displacements are interpreted by that Gly C ? O of the former is not involved in NH …? O ? C hydrogen bonds, while this carbon of the latter is linked by these kinds of hydrogen bonds. On the basis of these ¹³C chemical shifts, as reference data for the collagenlike structure, we were able to assign the ¹³C-nmr peaks of Gly, Ala, Pro, and Hyp residues of collagen fibrils, which are in good agreement with the values expected from the model polypeptides mentioned above. We also discuss a plausible conformational change of collagen fibrils during denaturation.     

On the presence of short-range periodicities in protein structures that are not related to established secondary structure elements

Standard secondary structure elements such as α-helices or β-sheets, are characterized by repeating backbone torsion angles (φ,ψ) at the single residue level. Two-residue motifs of the type (φ,ψ)₂ are also observed in nonlinear conformations, mainly turns. Taking these observations a step further, it can be argued that there is no a priori reason why the presence of higher order periodicities can not be envisioned in protein structures, such as, for example, periodic transitions between successive residues of the type (…-α-β-α-β-α-…), or (…-β-α_L-β-α_L-β-…), or (…-α-β-α_L-α-β-α_L-…), and so forth, where the symbols (α,β,α_L) refer to the established Ramachandran-based residue conformations. From all such possible higher order periodicities, here we examine the deposited (with the PDB) protein structures for the presence of short-range periodical conformations comprising five consecutive residues alternating between two (and only two) distinct Ramachandran regions, for example, conformations of the type (α-β-α-β-α) or (β-α_L-β-α_L-β), and so forth. Using a probabilistic approach, we have located several thousands of such peptapeptides, and these were clustered and analyzed in terms of their structural characteristics, their sequences, and their putative functional correlations using a gene ontology-based approach. We show that such nonstandard short-range periodicities are present in a large and functionally diverse sample of proteins, and can be grouped into two structurally conserved major types. Examination of the structural context in which these peptapeptides are observed gave no conclusive evidence for the presence of a persistent structural or functional role of these higher order periodic conformations.     

De novo design,synthesis and solution conformational study of two didehydroundecapeptides: effect of nature and number of amino acids interspersed between Phe residues

De novo design of peptides and proteins has recently surfaced as an approach for investigating protein structure and function. This approach vitally tests our knowledge of protein folding and function, while also laying the groundwork for the fabrication of proteins with properties not precedented in nature. The success relies heavily on the ability to design relatively short peptides that can espouse stable secondary structures. To this end, substitution with α,β‐didehydroamino acids, especially α,β‐didehydrophenylalanine (Δ^zPhe), comes in use for spawning well‐defined structural motifs. Introduction of ΔPhe induces β‐bends in small and 3₁₀‐helices in longer peptide sequences. The present work aims to investigate the effect of nature and the number of amino acids interspersed between two ΔPhe residues in two model undecapeptides, Ac‐Gly‐Ala‐ΔPhe‐Ile‐Val‐ΔPhe‐Ile‐Val‐ΔPhe‐Ala‐Gly‐NH₂ (I) and Boc‐Val‐ΔPhe‐Phe‐Ala‐Phe‐ΔPhe‐Phe‐Leu‐Ala‐ΔPhe‐Gly‐OMe (II). Peptide I was synthesized using solid‐phase chemistry and characterized using circular dichroism spectroscopy. Peptide II was synthesized using solution‐phase chemistry and characterized using circular dichroism and nuclear magnetic resonance spectroscopy. Peptide I was designed to examine the effect of incorporating β‐strand‐favoring residues like valine and isoleucine as spacers between two ΔPhe residues on the final conformation of the resulting peptide. Circular dichroism studies on this peptide have shown the existence of a 3₁₀‐helical conformation. Peptide II possesses three amino acids as spacers between ΔPhe residues and has been reported to adopt a mixed 3₁₀/α‐helical conformation using circular dichroism and nuclear magnetic resonance spectroscopy studies. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Identification of functionally important conserved trans‐membrane residues of bacterial PIB‐type ATPases

Powered by ATP hydrolysis, P_IB‐ATPases drive the energetically uphill transport of transition metals. These high affinity pumps are essential for heavy metal detoxification and delivery of metal cofactors to specific cellular compartments. Amino acid sequence alignment of the trans‐membrane (TM) helices of P_IB‐ATPases reveals a high degree of conservation, with ～60–70 fully conserved positions. Of these conserved positions, 6–7 were previously identified to be important for transport. However, the functional importance of the majority of the conserved TM residues remains unclear. To investigate the role of conserved TM residues of P_IB‐ATPases we conducted an extensive mutagenesis study of a Zn²⁺/Cd²⁺ P_IB‐ATPase from Rhizobium radiobacter (rrZntA) and seven other P_IB‐ATPases. Of the 38 conserved positions tested, 24 had small effects on metal tolerance. Fourteen mutations compromised in vivo metal tolerance and in vitro metal‐stimulated ATPase activity. Based on structural modelling, the functionally important residues line a constricted ‘channel’, tightly surrounded by the residues that were found to be inconsequential for function. We tentatively propose that the distribution of the mutable and immutable residues marks a possible trans‐membrane metal translocation pathway. In addition, by substituting six trans‐membrane amino acids of rrZntA we changed the in vivo metal specificity of this pump from Zn²⁺/Cd²⁺ to Ag⁺.     

Site-directed mutagenesis of the conserved Ala348 and Gly350 residues at the putative active site of Bacillus kaustophilus leucine aminopeptidase

Leucine aminopeptidase (LAP) is an exopeptidase that catalyzes the hydrolysis of amino acid residues from the amino terminus of proteins and peptides. Sequence alignment shows that the conserved Ala348 and Gly350 residues of Bacillus kaustophilus LAP (BkLAP) are located right next to a coordinated ligand. We further investigated the roles of these two residues by performing computer modeling and site-directed mutagenesis. Based on the modeling, the carbonyl group of Ala348 interacts with Asn345 and Asn435, and that of Gly350 with Ile353 and Leu354, where these interactions might maintain the zinc-coordinated residues at their correct positions. Replacement of Ala348 with Arg resulted in a dramatic reduction in LAP activity. A complete loss of the activity was also observed in A348E, A348V, and the Gly350 variants. Measurement of intrinsic tryptophan fluorescence revealed alteration of the microenvironment of aromatic amino acid residues, while circular dichroism spectra were nearly identical for wild-type and all mutant enzymes. Protein modeling and site-directed mutagenesis suggest that residues Ala348 and Gly350 are essential for BkLAP in maintaining a stable active-site environment for the catalytic reaction.     

The impact of β‐azido(or 1‐piperidinyl)methylamino acids in position 2 or 3 on biological activity and conformation of dermorphin analogues

The synthesis of new dermorphin analogues is described. The (R)‐alanine or phenylalanine residues of natural dermorphin were substituted by the corresponding α‐methyl‐β‐azidoalanine or α‐benzyl‐β‐azido(1‐piperidinyl)alanine residues. The potency and selectivity of the new analogues were evaluated by a competitive receptor binding assay in rat brain using [³H]DAMGO (a μ ligand) and [³H]DELT (a δ ligand). The most active analogue in this series, Tyr‐(R)‐Ala‐(R)‐α‐benzyl‐β‐azidoAla‐Gly‐Tyr‐Pro‐Ser‐NH₂ and its epimer were analysed by ¹H and ¹³C NMR spectroscopy and restrained molecular dynamics simulations. The dominant conformation of the investigated peptides depended on the absolute configuration around C^α in the α‐benzyl‐β‐azidoAla residue in position 3. The (R) configuration led to the formation of a type I β‐turn, whilst switching to the (S) configuration gave rise to an inverse β‐turn of type I′, followed by the formation of a very short β‐sheet. The selectivity of Tyr‐(R)‐Ala‐(R) and (S)‐α‐benzyl‐β‐azidoAla‐Gly‐Tyr‐Pro‐Ser‐NH₂ was shown to be very similar; nevertheless, the two analogues exhibited different conformational preferences. Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.     

NMR study of silk I structure of Bombyx mori silk fibroin with 15N- and 13C-NMR chemical shift contour plots

The metastable state silk I structures of Bombyx mori silk fibroin in the solid state were studied on the basis of ¹⁵N- and ¹³C-nmr chemical shifts of Ala, Ser, and Gly residues. The ¹⁵N cross-polarization magic angle spinning (CP/MAS) nmr spectra of the precipitated fraction after chymotrypsin hydrolysis of B. mori silk fibroin with the silk I and silk II forms were measured to determine the ¹⁵N chemical shifts of Gly, Ala, and Ser residues. For comparison, ¹⁵N CP/MAS nmr chemical shifts of Ala were measured for [¹⁵N] Ala Philosamia cynthia ricini silk fibroin with antiparallel β-sheet and α-helix forms. The ¹³C CP/MAS nmr chemical shifts of Ala, Ser, and Gly residues of B. mori silk fibroin with the silk I and silk II forms, as well as ¹³C CP/MAS nmr chemical shifts of Ala residue of P. c. ricini silk fibroin with β-sheet and α-helix forms, are used for the examination of the silk I structure. Both silk I and α-helix peaks are shifted to a lower field than silk II (β-sheet) for the Cα carbons of the Ala residues, while both Cβ carbon peaks are shifted to higher field. However, the silk I peak of the ¹⁵N nucleus of the Ala residue is shifted to lower field than the silk II peak, but the α-helix peak is shifted to high field. Thus, the difference in the structure between the silk I and α-helix is reflected in a different manner between the ¹³C and ¹⁵N chemical shifts. The Cα and Cβ chemical shift contour plots for Ala and Ser residues, and the Cα plot for the Gly residue, were prepared from the Protein Data Bank data obtained for 12 proteins and used for discussing the silk I structure quantitatively from the conformation-dependent chemical shifts. The plots reported by Le and Oldfield for ¹⁵N chemical shifts were also used for the purpose. All these chemical shift data support Fossey's model (Ala: ϕ = −80°, φ = 150°, Gly: ϕ = −150°, φ = 80°) and do not support Lotz and Keith's model (Ala: ϕ = −104.6°, φ = 112.2°, Gly: ϕ = 79.8°, φ = 49.7° or Ala: ϕ = −124.5°, φ = 88.2°, Gly: ϕ = −49.8°, φ = −76.1°) as the silk I structure. © 1997 John Wiley & Sons, Inc.     

Conformational preferences of amino acid side chains in collagen

Conformational energy computations were carried out on collagenlike triple-stranded conformations of several poly(tripeptide)s with the general structure CH₃CO? (Gly? X? Y)₃? NHCH₃. The sequences considered had various amino acid residues in position X or Y of the central tripeptide, with either Pro or Ala as a neighbor, i.e., Gly-X-Pro, Gly-X-Ala, Gly-Pro-Y, and Gly-Ala-Y. Minimum-energy conformations were computed for the side chains, and their distributions were compared for the four sequences. The residues used were Abu (= α-aminobutyric acid), Leu, Phe, Ser, Asp, Asn, Val, Ile, and Thr. The conformational energy of a ? Ch₂? CH₃ side chain in Abu was mapped as a function of the dihedral angle χ¹. Intrastrand interactions with neighboring residues do not affect the conformations of a side chain in position Y, and they have a minor effect on it in the X-Ala sequence, but they strongly restrict the conformational freedom of the side chain in the X-Pro sequence. Conversely, interstrand interactions do not affect side chains in position X, but they strongly restrict the conformational freedom of a side chain in position Y if there is a nearby Pro residue in a neighboring strand. Hydrogen bonds with the backbone can be formed in some conformations of long polar side chains, such as Asp, Asn, or Gln. All amino acid residues can be accommodated in collagen. Because of the interactions mentioned above, steric and energetic constraints can be correlated with observed preferences of certain amino acids for positions X or Y in collagen. Hence, these preferences may be explained, in part, in terms of differences in the conformational freedom of the side chains in the triple-stranded structure.     

Insect chymotrypsins: chloromethyl ketone inactivation and substrate specificity relative to possible coevolutional adaptation of insects and plants

Insect digestive chymotrypsins are present in a large variety of insect orders but their substrate specificity still remains unclear. Four insect chymotrypsins from 3 different insect orders (Dictyoptera, Coleoptera, and two Lepidoptera) were isolated using affinity chromatography. Enzymes presented molecular masses in the range of 20 to 31 kDa and pH optima in the range of 7.5 to 10.0. Kinetic characterization using different colorimetric and fluorescent substrates indicated that insect chymotrypsins differ from bovine chymotrypsin in their primary specificity toward small substrates (like N‐benzoyl‐L‐Tyr p‐nitroanilide) rather than on their preference for large substrates (exemplified by Succynil‐Ala‐Ala‐Pro‐Phe p‐nitroanilide). Chloromethyl ketones (TPCK, N‐ α‐tosyl‐L‐Phe chloromethyl ketone and Z‐GGF‐CK, N‐ carbobenzoxy‐Gly‐Gly‐Phe‐CK) inactivated all chymotrypsins tested. Inactivation rates follow apparent first‐order kinetics with variable second order rates (TPCK, 42 to 130 M^?1 s^?1; Z‐GGF‐CK, 150 to 450 M^?1 s^?1) that may be remarkably low for S. frugiperda chymotrypsin (TPCK, 6 M^?1 s^?1; Z‐GGF‐CK, 6.1 M^?1 s^?1). Homology modelling and sequence alignment showed that in lepidopteran chymotrypsins, differences in the amino acid residues in the neighborhood of the catalytic His 57 may affect its pKa value. This is proposed as the cause of the decrease in His 57 reactivity toward chloromethyl ketones. Such amino acid replacement in the active site is proposed to be an adaptation to the presence of dietary ketones. © 2009 Wiley Periodicals, Inc.     

Implication of an unfavorable residue (Thr346) in intrinsic flexibility of firefly luciferase

In order to better understand the functional role of an unusual residue (Thr346) of firefly luciferase mutagenesis at this residue was performed. Firefly luciferase, catalyzes the bioluminescence reaction and is an excellent tool as a reporter in nano-system biology studies. Nonetheless, the enzyme rapidly loses its activity at temperatures above 30°C and this leads to reduced sensitivity and precision in analytical applications. Residue Thr346 in a connecting loop (341-348) of firefly luciferase is located in a disallowed region of Ramachandran plot. In this study, we have substituted this residue (T346) with anomalous dihedral angles with Val, Gly and Pro to clarify the role of this residue in structure and function of the enzyme using site-directed mutagenesis. Substitution of this unfavorable residue (T346) with atypical dihedral angles (ψ, φ) with other residues brought about an increase of thermostability and decrease of specific activity. Structural and functional properties of the mutants were analyzed using different spectroscopic methods. It seems that this residue is a critically conserved residue to support the functional flexibility for a fast kinetic bioluminescence reaction at the expense of lower stability.