Diproline Templates as Folding Nuclei in Designed Peptides. Conformational Analysis of Synthetic Peptide Helices Containing Amino Terminal Pro-Pro Segments

The effect of N-terminal diproline segments in nucleating helical folding in designed peptides has been studied in two model sequences Piv-Pro-Pro-Aib-Leu-Aib-Phe-OMe (1) and Boc-Aib-Pro-Pro-Aib-Val-Ala-Phe-OMe (2). The structure of 1 in crystals, determined by X-ray diffraction, reveals a helical (RR) conformation for the segment residues 2 to 5, stabilized by one 4 -> 1 hydrogen bond and two 5 -> 1 interactions. The N-terminus residue, Pro(1) adopts a polyproline II (P-II) conformation. NMR studies in three different solvent systems support a conformation similar to that observed in crystals. In the apolar solvent CDCl3, NOE data favor the population of both completely helical and partially unfolded structures. In the former, the Pro-Pro segment adopts an alpha(R)-alpha(R) conformation, whereas in the latter, a P-II-alpha(R) structure is established. The conformational equilibrium shifts in favor of the P-II-alpha(R) structure in solvents like methanol and DMSO. A significant population of the Pro(1)- Pro(2) cis conformer is also observed. The NMR results are consistent with the population of at least three conformational states about Pro- Pro segment: trans alpha(R)-alpha(R), trans P-II-alpha(R) and cis P-II-alpha(R). Of these, the two trans conformers are in rapid dynamic exchange on the NMR time scale, whereas the interconversion between cis and trans form is slow. Similar results are obtained with peptide 2. Analysis of 462 diproline segments in protein crystal structures reveals 25 examples of the alpha(R)-alpha(R) conformation followed by a helix. Modeling and energy minimization studies suggest that both P-II-alpha(R) and alpha(R)-alpha(R) conformations have very similar energies in the model hexapeptide 1

Structural Investigations on Poly(4-hydroxy-L-proline). 1. Theoretical Studies

Model building studies on poly(hydroxypro1ine) indicate that in addition to the well-known helical structure of form A, a left-handed helical structure with trans peptide units and with h = 2.86 A and n = 2.67 (i.e., 8 residues in 3 turns) is also possible. In this structure which is shown to be in agreement with X-ray data of the form B in the next paper, the y-hydroxyl group of an (i + 1)th Hyp residue is hydrogen bonded to the carbonyl oxygen of an (i - 1)th residue. The possibility of a structure with cis peptide units is ruled out. It is shown that both forms A and B are equally favorable from considerations of intramolecular energies. Since form B is further stabilized by intrachain hydrogen bonds, we believe that this is likely to be the ordered conformation for poly(hydroxypro1ine) in water.

Modes of binding of guanosine monophosphates to ribonuclease T1 - A computer modelling study

Computer-modelling studies on the modes of binding of the three guanosine monophosphate inhibitors 2'-GMP, 3'-GMP, and 5'-GMP to ribonuclease (RNase) T1 have been carried out by energy minimization in Cartesian-coordinate space. The inhibitory power was found to decrease in the order 2'-GMP > 3'-GMP > 5'-GMP in agreement with the experimental observations. The ribose moiety was found to form hydrogen bonds with the protein in all the enzyme-inhibitor complexes, indicating that it contributes to the binding energy and does not merely act as a spacer between the base and the phosphate moieties as suggested earlier. 2'-GMP and 5'-GMP bind to RNase T1 in either of the two ribose puckered forms (with C3'-endo more favoured over the C2'-endo) and 3'-GMP binds to RNase T1 predominantly in C3'-endo form. The catalytically important residue His-92 was found to form hydrogen bond with the phosphate moiety in all the enzyme-inhibitor complexes, indicating that this residue may serve as a general acid group during catalysis. Such an interaction was not found in either X-ray or two-dimensional NMR studies.

Identification and characterization of a library of microheterogeneous cyclohexa-depsipeptides from the fungus

Ten new cyclic hexadepsipeptides, six isariins and four isaridins, from the fungus Isaria have been identified and characterized by high-performance liquid chromatography, coupled to tandem electrospray ionization mass spectrometry (LC-ESIMS/MS). The isariins possess a beta-hydroxy acid residue and five alpha-amino acids, while isaridins contain a beta-amino acid, an alpha-hydroxy acid, and four alpha-amino acids. One- and two-dimensional NMR spectroscopy confirmed the chemical identity of some of the isariin fractions. Mass spectral fragmentation patterns of [M + H](+) ions reveal clear diagnostic fragment ions for the isariins and isaridins. Previously described cyclic depsipeptides, isarfelins from Isaria felina (Guo, Y. X.; Liu, Q. H.; Ng, T. B.; Wang H. X. Peptides 2005, 26, 2384), are now reassigned as members of the isaridin family. Examination of isaridin sequences revealed significant similarities with cyclic hexadepsipeptides such as destruxins and roseotoxins. The structure of an isariin (isariin A) investigated by NMR spectroscopy indicated the presence of a hybrid alpha beta C-11 turn, formed by the beta-hydroxy acid and glycine residues and a (D)Leu-(L)Ala type II' beta-turn. Additionally, the inhibitory effect of isariins and an isaridin on the intra-erythrocytic growth of Plasmodium falciparum is presented.

Structure of bovine prothrombin fragment 1 refined at 2.25 Å resolution

The structure of bovine prothrombin fragment 1 has been refined at 2.25 Å resolution using high resolution measurements made with the synchrotron beam at CHESS. The synchrotron data were collected photographically by oscillation methods (R-merge = 0.08). These were combined with lower order diffractometer data for refinement purposes. The structure was refined using restrained least-squares methods with the program PROLSQ to a crystallographic R-value of 0.175. The structure includes 105 water molecules with occupancies of >0·6. The first 35 residues (Ala1-Leu35) of the N-terminal ?-carboxy glutamic acid-domain (Ala1-Cys48) of fragment 1 are disordered as are two carbohydrate chains of Mr ? 5000; the latter two combine to render 40% of the structure disordered. The folding of the kringle of fragment 1 is related to the close intramolecular contact between the inner loop disulfide groups. Half of the conserved sequence of the kringle forms an inner core surrounding these disulfide groups. The remainder of the sequence conservation is associated with the many turns of the main chain. The Pro95 residue of the kringle has a cis conformation and Tyr74 is ordered in fragment 1, although nuclear magnetic resonance studies indicate that the comparable residue of plasminogen kringle 4 has two positions. Surface accessibility calculations indicate that none of the disulfide groups of fragment 1 is accessible to solvent.

A Deterministic Optimization Approach to Protein Sequence Design Using Continuous Models

Determining the sequence of amino acid residues in a heteropolymer chain of a protein with a given conformation is a discrete combinatorial problem that is not generally amenable for gradient-based continuous optimization algorithms. In this paper we present a new approach to this problem using continuous models. In this modeling, continuous "state functions" are proposed to designate the type of each residue in the chain. Such a continuous model helps define a continuous sequence space in which a chosen criterion is optimized to find the most appropriate sequence. Searching a continuous sequence space using a deterministic optimization algorithm makes it possible to find the optimal sequences with much less computation than many other approaches. The computational efficiency of this method is further improved by combining it with a graph spectral method, which explicitly takes into account the topology of the desired conformation and also helps make the combined method more robust. The continuous modeling used here appears to have additional advantages in mimicking the folding pathways and in creating the energy landscapes that help find sequences with high stability and kinetic accessibility. To illustrate the new approach, a widely used simplifying assumption is made by considering only two types of residues: hydrophobic (H) and polar (P). Self-avoiding compact lattice models are used to validate the method with known results in the literature and data that can be practically obtained by exhaustive enumeration on a desktop computer. We also present examples of sequence design for the HP models of some real proteins, which are solved in less than five minutes on a single-processor desktop computer Some open issues and future extensions are noted.

Functional Mimics of Glutathione Peroxidase: Bioinspired Synthetic Antioxidants

Oxidative stress is caused by an imbalance between the production of reactive oxygen species (ROS) and the biological system's ability to detoxify these reactive intermediates. Mammalian cells have elaborate antioxidant defense mechanisms to control the damaging effects of ROS. Glutathione peroxidase (GPx), a selenoenzyme, plays a key role in protecting the organism from oxidative damage by catalyzing the reduction of harmful hydroperoxides with thiol a ``catalytic triad'' with tryptophan and glutamine, which cofactors. The selenocysteine residue at the active site forms activates the selenium moiety for an efficient reduction of peroxides. After the discovery that ebselen, a synthetic organoselenium compound, mimics the catalytic activity of GPx both in vitro and in vivo, several research groups developed a number of small-molecule selenium compounds as functional mimics of GPx, either by modifying the basic structure of ebselen or by incorporating some structural features of the native enzyme. The synthetic mimics reported in the literature can be classified in three major categories: (i) cyclic selenenyl amides having a Se-N bond, (ii) diaryl diselenides, and (iii) aromatic or aliphatic monoselenides. Recent studies show that ebselen exhibits very poor GPx activity when aryl or benzylic thiols such as PhSH or BnSH are used as cosubstrates. Because the catalytic activity of each GPx mimic largely depends on the thiol cosubstrates used, the difference in the thiols causes the discrepancies observed in different studies. In this Account, we demonstrate the effect of amide and amine substituents on the GPx activity of various organoselenium compounds. The existence of strong Se ... O/N interactions in the selenenyl sulfide intermediates significantly reduces the GPx activity. These interactions facilitate an attack of thiol at selenium rather than at sulfur, leading to thiol exchange reactions that hamper the formation of catalytically active selenol. Therefore, any substituent capable of enhancing the nucleophilic attack of thiol at sulfur in the selenenyl sulfide state would enhance the antioxidant potency of organoselenium compounds. Interestingly, replacement of the sec-amide substituent by a tert-amide group leads to a weakening of Se ... 0 interactions in the selenenyl sulfide intermediates. This modification results in 10- to 20-fold enhancements in the catalytic activities. Another strategy involving the replacement of tert-amide moieties by tert-amino substituents further increases the activity by 3- to 4-fold. The most effective modification so far in benzylamine-based GPx mimics appears to be either the replacement of a tert-amino substituent by a sec-amino group or the introduction of an additional 6-methoxy group in the phenyl ring. These strategies can contribute to a remarkable enhancement in the GPx activity. In addition to enhancing catalytic activity, a change in the substituents near the selenium moiety alters the catalytic mechanisms. The mechanistic investigations of functional mimics are useful not only for understanding the complex chemistry at the active site of GPx but also for designing and synthesizing novel antioxidants and anti-inflammatory agents.

HELANAL: A program to characterise helix geometry in proteins,

A detailed analysis of structural and position dependent characteristic features of helices will give a better understanding of the secondary structure formation in globular proteins. Here we describe an algorithm that quantifies the geometry of helices in proteins on the basis of their C-alpha atoms alone. The Fortran program HELANAL can extract the helices from the PDB files and then characterises the overall geometry of each helix as being linear, curved or kinked, in terms of its local structural features, viz. local helical twist and rise, virtual torsion angle, local helix origins and bending angles between successive local helix axes. Even helices with large radius of curvature are unambiguously identified as being linear or curved. The program can also be used to differentiate a kinked helix and other motifs, such as helix-loop-helix or a helix-turn-helix (with a single residue linker) with the help of local bending angles. In addition to these, the program can also be used to characterise the helix start and end as well as other types of secondary structures.

Error-free matrix symmetrizers and equivalent symmetric matrices

A symmetrizer of the matrix A is a symmetric solution X that satisfies the matrix equation XA=AprimeX. An exact matrix symmetrizer is computed by obtaining a general algorithm and superimposing a modified multiple modulus residue arithmetic on this algorithm. A procedure based on computing a symmetrizer to obtain a symmetric matrix, called here an equivalent symmetric matrix, whose eigenvalues are the same as those of a given real nonsymmetric matrix is presented.

Conformational studies on beta-bend containing a cis peptide unit

Conformational studies have been carried out on the X-cis-Pro tripeptide system (a system of three linked peptide units, in the trans-cis-trans configuration) using energy minimization techniques. For X, residues Gly, L-Ala, D-Ala and L-Pro have been used. The energy minima have been classified into different groups based upon the conformational similarity. There are 15, 20, 18 and 6 minima that are possible for the four cases respectively and these fall into 11 different groups. A study of these minima shows that, (i) some minima contain hydrogen bonds - either 4-->1 or 1-->2 type, (ii) the low energy minima qualify themselves as bend conformations, (iii) cis' and trans' conformations are possible for the prolyl residue as also the C(gamma)-endo and C(gamma)-exo puckerings, and (iv) for Pro-cis-Pro, cis' at the first prolyl residue is ruled out, due to the high energy. The available crystal structure data on proteins and peptides, containing cis-Pro segment have been examined with a view to find the minima that occur in solid state. The data from protein show that they fall under two groups. The conformation at X in X-cis-Pro is near extended when it is a non-glycyl residue. In both peptides and proteins there exists a preference for trans' conformation at prolyl residue over cis' when X is a non-glycyl residue. The minima obtained can be useful in modelling studies.

(2E)-3-(4-Ethoxyphenyl)-1-(2-methyl-4-phenylquinolin-3-yl)prop-2-en-1-on e monohydrate

The title hydrate, C27H23NO2 center dot H2O, features an almost planar quinoline residue (r.m.s. deviation = 0.015 angstrom) with the benzene dihedral angle = 63.80 (7)degrees] and chalcone C-C-C-O torsion angle = -103.38 (18)degrees] substituents twisted significantly out of its plane. The configuration about the C=C bond 1.340 (2) angstrom] is E. In the crystal, molecules related by the 21 symmetry operation are linked along the b axis via water molecules that form O-H center dot center dot center dot O-c and O-H center dot center dot center dot N-q hydrogen bonds (c = carbonyl and q = quinoline). A C-H center dot center dot center dot O interaction also occurs.

Effect of constraints by threonine on proline containing αa-helix—A molecular dynamics approach

Proline plays an important role in the secondary structure of proteins. In the pursuit of understanding its structural role, Proline containing helices with constraints have been studied by employing molecular dynamics (MD) technique. In the present study, the constraint introduced is a threonine residue, whose sidechain has intramolecular hydrogen bond interaction with the backbone oxygen atom. The three systems that have been chosen for characterization are: (1) Ace-(Ala)12−Thr-Pro-(Ala)10−NHMe, (2) Ace-(Ala)13-Pro-Ala-Thr- (Ala)8-NHMe and (3) Ace-(Ala)13-Pro-(Ala)3-Thr-(Ala)6-NHMe. The equilibrium structures and structural transitions have been identified by monitoring the backbone dihedral angles, bend related parameters and the hydrogen bond interactions. The MD averages and root mean square (r.m.s.) fluctuations are compared and discussed. Energy minimization has been carried out on selected MD simulated points in order to analyze the characteristics of different conformations.

Identification of Critical Residues of the Mycobacterial Dephosphocoenzyme A Kinase by Site-Directed Mutagenesis

Dephosphocoenzyme A kinase performs the transfer of the c-phosphate of ATP to dephosphocoenzyme A, catalyzing the last step of coenzyme A biosynthesis. This enzyme belongs to the P-loop-containing NTP hydrolase superfamily, all members of which posses a three domain topology consisting of a CoA domain that binds the acceptor substrate, the nucleotide binding domain and the lid domain. Differences in the enzymatic organization and regulation between the human and mycobacterial counterparts, have pointed out the tubercular CoaE as a high confidence drug target (HAMAP database). Unfortunately the absence of a three-dimensional crystal structure of the enzyme, either alone or complexed with either of its substrates/regulators, leaves both the reaction mechanism unidentified and the chief players involved in substrate binding, stabilization and catalysis unknown. Based on homology modeling and sequence analysis, we chose residues in the three functional domains of the enzyme to assess their contributions to ligand binding and catalysis using site-directed mutagenesis. Systematically mutating the residues from the P-loop and the nucleotide-binding site identified Lys14 and Arg140 in ATP binding and the stabilization of the phosphoryl intermediate during the phosphotransfer reaction. Mutagenesis of Asp32 and Arg140 showed catalytic efficiencies less than 5-10% of the wild type, indicating the pivotal roles played by these residues in catalysis. Non-conservative substitution of the Leu114 residue identifies this leucine as the critical residue from the hydrophobic cleft involved in leading substrate, DCoA binding. We show that the mycobacterial enzyme requires the Mg2+ for its catalytic activity. The binding energetics of the interactions of the mutant enzymes with the substrates were characterized in terms of their enthalpic and entropic contributions by ITC, providing a complete picture of the effects of the mutations on activity. The properties of mutants defective in substrate recognition were consistent with the ordered sequential mechanism of substrate addition for CoaE.

Amino acid selective unlabeling for sequence specific resonance assignments in proteins

Sequence specific resonance assignment constitutes an important step towards high-resolution structure determination of proteins by NMR and is aided by selective identification and assignment of amino acid types. The traditional approach to selective labeling yields only the chemical shifts of the particular amino acid being selected and does not help in establishing a link between adjacent residues along the polypeptide chain, which is important for sequential assignments. An alternative approach is the method of amino acid selective `unlabeling' or reverse labeling, which involves selective unlabeling of specific amino acid types against a uniformly C-13/N-15 labeled background. Based on this method, we present a novel approach for sequential assignments in proteins. The method involves a new NMR experiment named, {(CO)-C-12 (i) -N-15 (i+1)}-filtered HSQC, which aids in linking the H-1(N)/N-15 resonances of the selectively unlabeled residue, i, and its C-terminal neighbor, i + 1, in HN-detected double and triple resonance spectra. This leads to the assignment of a tri-peptide segment from the knowledge of the amino acid types of residues: i - 1, i and i + 1, thereby speeding up the sequential assignment process. The method has the advantage of being relatively inexpensive, applicable to H-2 labeled protein and can be coupled with cell-free synthesis and/or automated assignment approaches. A detailed survey involving unlabeling of different amino acid types individually or in pairs reveals that the proposed approach is also robust to misincorporation of N-14 at undesired sites. Taken together, this study represents the first application of selective unlabeling for sequence specific resonance assignments and opens up new avenues to using this methodology in protein structural studies.

Analysis of the effects of amino-acid-sequence on the structure of proteins

The conformation of amino acid side chains as observed in well-determined structures of globular proteins has earlier been extensively investigated. In contrast, the structural features of the polypeptide backbone that result from the occurrence of specific amino acids along the polypeptide have not been analysed. In this article, we present the statistically significant features in the backbone geometry that appear to be a consequence of the occurrence of rotamers of different amino acid side chains by analysing 102 well-refined structures that form a random collection of proteins. It is found that the persistence of helical segments around each residue is influenced by the residue type. Several residues exert asymmetrical influence between the carboxyl and amino terminal polypeptide segments. The degree to which secondary structures depart from an average geometry also appears to depend on residue type. These departures are correlated to the corresponding Chou and Fasman parameters of amino acid residues. The frequency distribution of the side chain rotamers is influenced by polypeptide secondary structure. In turn, the rotamer conformation of side chain affects the extension of the secondary structure of the backbone. The strongest correlation is found between the occurrence of g+ conformation and helix propagation on the carboxyl side of many residues.