Identification of amino acid residues essential for enzyme activity of sheep liver 5,10-methylenetetrahydrofolate reductase

Sheep liver 5,10-methylenetetrahydrofolate reductase was subjected to specific chemical modification with phenylglyoxal, diethyl pyrocarbonate and N-bromosuccinimide. The second-order rate constants for inactivation were calculated to be 54 M-1 X min-1, 103 M-1 X min-1 and 154 M-1 X min-1 respectively. This inactivation could be prevented by incubation with substrates or products, suggesting that the residues modified, namely arginine, histidine and tryptophan, are essential for enzyme activity.

Stereochemical restrictions on the occurrence of amino acid residues in the collagen structure

The primary structure of collagen is characterized by the repeating tripeptide sequence (Gly-R2-R3)n. The results of theoretical studies, carried out using contact criteria to compute the stereochemically allowed orientations for various side chains at locations 2 and 3, are reported here. It is found that side chains with only γ-atoms, as in valine, serine and threonine, or with only one δ-methyl group, as in isoleucine, can occur equally well at locations 2 and 3, as is actually the case in collagen. Side chains with two Cδ-atoms, as in leucine and phenyl-alanine, can also be accommodated at both positions. However, if they occur as R3 their freedom of orientation is severely restricted in the presence of a proline residue as R2 in a neighbouring chain. If water molecules bound to the chains of the triple helix are assumed to be present, then location 3 is virtually impossible for leucine and phenylalanine residues. Location 2 is, however, unaffected, and their presence as R2 can help to shield the water molecules from disturbance by the solvent medium. This may be the reason for the preferential occurrence of Leu and Phe residues in location 2 in the collagen triplets, although the polypeptides (Gly-Pro-Leu)n and (Gly-Pro-Phe)n form collagen-like structures.

Disallowed Ramachandran Conformations of Amino Acid Residues in Protein Structures

An analysis of the nature and distribution of disallowed Ramachandran conformations of amino acid residues observed in high resolution protein crystal structures has been carried out. A data set consisting of 110 high resolution, non-homologous, protein crystal structures from the Brookhaven Protein Data Bank was examined. The data set consisted of a total of 18,708 non-Gly residues, which were characterized on the basis of their backbone dihedral angles (φ, ψ). Residues falling outside the defined “broad allowed limits” on the Ramachandran map were chosen and the reportedB-factor value of the α-carbon atom was used to further select well defined disallowed conformations. The conformations of the selected 66 disallowed residues clustered in distinct regions of the Ramachandran map indicating that specific φ, ψ angle distortions are preferred under compulsions imposed by local constraints. The distribution of various amino acid residues in the disallowed residue data set showed a predominance of small polar/charged residues, with bulky hydrophobic residues being infrequent. As a further check, for all the 66 cases non-hydrogen van der Waals short contacts in the protein structures were evaluated and compared with the ideal “Ala-dipeptide” constructed using disallowed dihedral angle (φ, ψ) values. The analysis reveals that short contacts are eliminated in most cases by local distortions of bond angles. An analysis of the conformation of the identified disallowed residues in related protein structures reveals instances of conservation of unusual stereochemistry.

Empirical torsional potential functions from protein structure data. Phi- and psi-potentials for non-glycyl amino acid residues.

The torsional potential functions Vt(phi) and Vt(psi) around single bonds N--C alpha and C alpha--C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (phi, psi)-plane with the value of Vtot(phi, psi), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in L-configuration, are Vt(phi) = 1.0 cos (phi + 60 degrees); Vt(psi) = 0.5 cos (psi + 60 degrees) - 1.0 cos (2 psi + 30 degrees) - 0.5 cos (3 psi + 30 degrees). The dipeptide energy maps Vtot(phi, psi) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line psi = 0 degrees. These functions derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.

Two Novel Hexadepsipeptides with Several Modified Amino Acid Residues Isolated from the Fungus

Two new cyclohexadepsipeptides have been isolated from the fungus Isaria. Fungal growth in solid media yielded hyphal strands from which peptide fractions were readily isolable by organic-solvent extraction. Two novel cyclodepsipeptides, isaridin A and isaridin B, have been isolated by reverse-phase HPLC, and characterized by ESI-MS and 1H-NMR. Single crystals of both peptides have been obtained, and their 3D structures were elucidated by X-ray diffraction. The isaridins contain several unusual amino acid residues. The sequences are cyclo(β-Gly-HyLeu-Pro-Phe-NMeVal-NMePhe) and cyclo(β-Gly-HyLeu-β-MePro-Phe-NMeVal-NMePhe), where NMeVal is N-methylvaline, NMePhe N-methylphenylalanine, and HyLeu hydroxyleucine (=2-hydroxy-4-methylpentanoic acid). The two peptides differ from one another at residue 3, isaridin A having an (S)-proline at this position, while β-methyl-(S)-proline (=(2S,3S)-2,3,4,5-tetrahydro-3-methyl-1H-pyrrole-2-carboxylic acid) is found in isaridin B. The solid-state conformations of both cyclic depsipeptides are characterized by the presence of two cis peptide bonds at HyLeu(2)-Pro(3)/HyLeu(2)-β-MePro(3) and NMeVal(5)-NMePhe(6), respectively. In isaridin A, a strong intramolecular H-bond is observed between Phe(4)CO⋅⋅⋅HNβ-Gly(1), and a similar, but weaker, interaction is observed between β-Gly(1)CO⋅⋅⋅HNPhe(4). In contrast, in isaridin B, only a single intramolecular H-bond is observed between β-Gly(1)CO⋅⋅⋅HNPhe(4

Empirical torsional potential functions from protiesn structure DATA Phi- and Psi-Potentials for Non-glycyl Amino Acid Residues

The torsional potential functions Vt(φ) and Vt(ψ) around single bonds N–Cα and Cα-C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (φ, ψ)-plane with the value of Vtot(φ, ψ), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in l-configuration, are Vt(φ) = – 1.0 cos (φ + 60°); Vt(ψ) = – 0.5 cos (ψ + 60°) – 1.0 cos (2ψ + 30°) – 0.5 cos (3ψ + 30°). The dipeptide energy maps Vtot(φ, ψ) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line ψ = 0°. These functions, derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.

Identification of key amino acid residues in the catalytic mechanism of diaminopropionate ammonialyase from Salmonella typhimurium

Diaminopropionate ammonialyase (DAPAL), a fold-typeII pyridoxal 5-phosphate-dependent enzyme, catalyzes the ,-elimination of diaminopropionate (DAP) to pyruvate and ammonia. DAPAL was able to utilize both d- and l-DAP as substrates with almost equal efficiency. Mutational analysis of functionally important residues such as Thr385, Asp125 and Asp194 was carried out to understand the mechanism by which the isomers are hydrolyzed. Further, the putative residues involved in the formation of disulfide bond Cys271 and Cys299 were also mutated. T385S, T385D sDAPAL were as active with dl-DAP as substrate as sDAPAL, whereas the later exhibited a threefold increase in catalytic efficiency with d-Ser as substrate. Further analysis of these mutants suggested that DAPAL might follow an anti-E-2 mechanism of catalysis that does not involve the formation of a quinonoid intermediate. Of the two mutants of Asp125, D125E showed complete loss of activity with d-DAP as substrate, whereas the reaction with l-DAP was not affected significantly, demonstrating that Asp125 was essential for abstraction of protons from the d-isomer. By contrast, mutational analysis of Asp194 showed that the residue may not be directly involved in proton abstraction from l-DAP. sDAPAL does not form a disulfide bond in solution, although the position of Cys299 and Cys271 in the modeled structure of sDAPAL favored the formation of a disulfide bond. Further, unlike eDAPAL, sDAPAL could be activated by monovalent cations. Mutation of the cysteine residues showed that Cys271 may be involved in coordinating the monovalent cation, as observed in the case of other fold-typeII enzymes.

PROPENSITIES OF AMINO ACID RESIDUES IN PROTEINS FOR DIFFERENT REGIONS OF THE RAMACHANDRAN MAP

Using a dataset of 1164 crystal structures of largely non-homologous proteins defined at a resolution of 1.5 angstrom or better, we have investigated the (phi,psi) preferences of 20 residue types by considering the residues which occur in loops. Propensities of residue types to occur in the loops with (phi,psi) values in the aa region of the Ramachandran map has a poor correlation coefficient of 0.48 to the Chou-Fasman propensities of the residue types to occur in the a-helical segments. However the correlation coefficient between propensities of residues in loops to adopt beta conformations and those in beta-sheet is much higher (0.95). These observations suggest that a-helix formation is well influenced by the local amino acid sequence while intrinsic preference of residue types for beta-sheet plays a major role in the formation of beta-sheet. The main chain polar groups of residues in loops, that can affect the (phi,psi) values, can be involved in intra-molecular hydrogen bonding. Therefore we investigated further by considering subset of residues in loops with low (0 to 2) number of intra-molecular hydrogen bonds per residue involving main chain polar atoms. For this subset, the correlation coefficients between propensities for alpha-helix and alpha(R) region and between beta-sheet and beta-region are 0.26 and 0.64 respectively. This reiterates higher intrinsic tendency of beta-region favouring residues to adopt beta-sheet than alpha(R) region favouring residues to adopt alpha-helical structure.

Identification of key amino acid residues in the catalytic mechanism of diaminopropionate ammonialyase from Salmonella typhimurium

Diaminopropionate ammonialyase (DAPAL), a fold-typeII pyridoxal 5-phosphate-dependent enzyme, catalyzes the ,-elimination of diaminopropionate (DAP) to pyruvate and ammonia. DAPAL was able to utilize both d- and l-DAP as substrates with almost equal efficiency. Mutational analysis of functionally important residues such as Thr385, Asp125 and Asp194 was carried out to understand the mechanism by which the isomers are hydrolyzed. Further, the putative residues involved in the formation of disulfide bond Cys271 and Cys299 were also mutated. T385S, T385D sDAPAL were as active with dl-DAP as substrate as sDAPAL, whereas the later exhibited a threefold increase in catalytic efficiency with d-Ser as substrate. Further analysis of these mutants suggested that DAPAL might follow an anti-E-2 mechanism of catalysis that does not involve the formation of a quinonoid intermediate. Of the two mutants of Asp125, D125E showed complete loss of activity with d-DAP as substrate, whereas the reaction with l-DAP was not affected significantly, demonstrating that Asp125 was essential for abstraction of protons from the d-isomer. By contrast, mutational analysis of Asp194 showed that the residue may not be directly involved in proton abstraction from l-DAP. sDAPAL does not form a disulfide bond in solution, although the position of Cys299 and Cys271 in the modeled structure of sDAPAL favored the formation of a disulfide bond. Further, unlike eDAPAL, sDAPAL could be activated by monovalent cations. Mutation of the cysteine residues showed that Cys271 may be involved in coordinating the monovalent cation, as observed in the case of other fold-typeII enzymes.

Conformational Choices for the Stereochemically Constrained gamma-Amino Acid Residue Gabapentin: Theoretical Studies and Correlation With Experimental Results

Gabapentin (1-aminomethylcyclohexaneacetic acid, Gpn) is an achiral, conformationally constrained gamma amino acid residue. A survey of available crystal structures of Gpn peptides reveals that the torsion angles about the C-gamma-C-beta (theta(1)) and C-beta-C-alpha(theta(2)) bonds are overwhelmingly limited to gauche, gauche (g(+)g(+)/g(-)g(-)) conformations. The Gpn residue forms C-7 and C-9 hydrogen bonds in which the donor and acceptor atoms come from the flanking peptide units. In combination with alpha amino acid residues alpha gamma and gamma alpha segments can adopt C-12 hydrogen bonded structures. The conformational choices available to the Gpn residue have been probed using energy calculations, adopting a grid search strategy. Ramachandran phi-psi maps have been constructed for fixed values of theta(1) and theta(2), corresponding to the gauche and trans conformations. The sterically allowed and energetically favorable regions of conformational space have been defined and experimental observations compared. C-7 and C-9 hydrogen bonded conformational families have been identified using a grid search approach in which theta(1) and theta(2) values are varied over a range of +/- 10 degrees about ideal values at 1 degrees intervals. The theoretical analysis together with experimental observations for 59 Gpn residues from 35 crystal structures permits definition of the limited range of conformational possibilities at this gamma amino acid residue. .

Gabapentin: A Stereochemically Constrained gamma Amino Acid Residue in Hybrid Peptide Design

Nature has used the all-alpha-polypeptide backbone of proteins to create a remarkable diversity of folded structures. Sequential patterns of 20 distinct amino adds, which differ only in their side chains, determine the shape and form of proteins. Our understanding of these specific secondary structures is over half a century old and is based primarily on the fundamental elements: the Pauling alpha-helix and beta-sheet. Researchers can also generate structural diversity through the synthesis of polypeptide chains containing homologated (omega) amino acid residues, which contain a variable number of backbone atoms. However, incorporating amino adds with more atoms within the backbone introduces additional torsional freedom into the structure, which can complicate the structural analysis. Fortunately, gabapentin (Gpn), a readily available bulk drug, is an achiral beta,beta-disubstituted gamma amino add residue that contains a cyclohexyl ring at the C-beta carbon atom, which dramatically limits the range of torsion angles that can be obtained about the flanking C-C bonds. Limiting conformational flexibility also has the desirable effect of increasing peptide crystallinity, which permits unambiguous structural characterization by X-ray diffraction methods. This Account describes studies carried out in our laboratory that establish Gpn as a valuable residue in the design of specifically folded hybrid peptide structures. The insertion of additional atoms into polypeptide backbones facilitates the formation of intramolecular hydrogen bonds whose directionality is opposite to that observed in canonical alpha-peptide helices. If hybrid structures mimic proteins and biologically active peptides, the proteolytic stability conferred by unusual backbones can be a major advantage in the area of medicinal chemistry. We have demonstrated a variety of internally hydrogen-bonded structures in the solid state for Gpn-containing peptides, including the characterization of the C-7 and C-9 hydrogen bonds, which can lead to ribbons in homo-oligomeric sequences. In hybrid alpha gamma sequences, district C-12 hydrogen-bonded turn structures support formation of peptide helices and hairpins in longer sequences. Some peptides that include the Gpn residue have hydrogen-bond directionality that matches alpha-peptide helices, while others have the opposite directionality. We expect that expansion of the polypeptide backbone will lead to new classes of foldamer structures, which are thus far unknown to the world of alpha-polypeptides. The diversity of internally hydrogen-bonded structures observed in hybrid sequences containing Gpn shows promise for the rational design of novel peptide structures incorporating hybrid backbones.

Hybrid peptide design. Hydrogen bonded conformations in peptides containing the stereochemically constrained gamma-amino acid residue, gabapentin

he crystal structure of 12 peptides containing the conformationally constrained 1-(aminomethyl)cyclohexaneacetic acid, gabapentin (Gpn), are reported. In all the 39 Gpn residues conformationally characterized so far, the torsion angles about the C-alpha-C-beta and C-beta-C-gamma bonds are restricted to the gauche conformation (+/- 60 degrees). The Gpn residue is constrained to adopt folded conformations resulting in the formation of intramolecularly hydrogen-bonded structures even in short peptides. The peptides Boc-Ac(6)c-Gpn-OMe 1 and Boc-Gpn-Aib-Gpn-Aib-OMe 2 provide examples of C-7 conformation; peptides Boc-Gpn-Aib-OH 3, Boc-Ac(6)c-Gpn-OH 4, Boc-Val-Pro-Gpn-OH 5, Piv-Pro-Gpn-Val-OMe 6, and Boc-Gpn-Gpn-Leu-OMe 7 provide examples of C-9 conformation; peptide Boc-Ala-Aib-Gpn-Aib-Ala-OMe 8 provides an example of C-12 conformation and peptides Boc-beta Leu-Gpn-Val-OMe 9 and Boc-beta Phe-Gpn-Phe-OMe 10 provide examples of C-13 conformation. Gpn peptides provide examples of backbone expanded mimetics for canonical alpha-peptide turns like the gamma (C-7) and the beta (C-10) turns. The hybrid beta gamma sequences provide an example of a mimetic of the C-13 alpha-turn formed by three contiguous alpha-amino acid residues. Two examples of folded tripeptide structures, Boc-Gpn-beta Phe-Leu-OMe 11 and Boc-Aib-Gpn-beta Phg-NHMe 12, lacking internal hydrogen bonds are also presented. An analysis of available Gpn residue conformations provides the basis for future design of folded hybrid peptides.

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

C-H…O hydrogen bond mediated chain reversal in a peptide containing g-amino acid residue, determined directly from powder X-ray diffraction data

The finding that peptides containing -amino acid residues give rise to folding patterns hitherto unobserved in -amino acid peptides[1] has stimulated considerable interest in the conformational properties of peptides built from , and residues,[2] as the introduction of additional methylene (CH2) units into peptide chains provides further degrees of conformational freedom.