Structure of porcine pancreatic elastase complexed with FR901277, a novel macrocyclic inhibitor of elastases, at 1.6 A resolution

Human leukocyte elastase (HLE) is a serine protease that contributes to tissue destruction in various disease states-for example, in emphysema. FR901277 is a natural product isolated from the culture filtrate of Streptomyces resistomicificus and is a potent inhibitor of both HLE and porcine pancreatic elastase (PPE). FR901277 consists of four normal amino acids and three unusual amino acids, and is a unique bicyclic peptide compound. The crystal structure of PPE complexed with FR901277 has been determined at 1.6 A resolution. The Ogamma atom of Ser-195 in PPE did not form a covalent bond with FR901277, but formed a hydrogen bond with the Nvarepsilon atom of His-57. On the other hand, the portion from L-Orn(1) through dehydroxyThr(3) in FR901277 formed an antiparallel beta-sheet structure with the backbone of the active site in PPE. The S4 through S2' binding subsites in PPE were all occupied by the hydrophobic side chains of the inhibitor molecule. Especially, the ethylidene moiety of FR901277 occupied the S1 specific pocket, indicating a CH/pi interaction. In addition, the isopropyl side chain of L-Val(7) was located at the enzyme surface between the S2 and S1' pockets with several van der Waals contacts. However, the amino acid (4) residue was not involved in a significant interaction with PPE. Comparison of inhibitor structures in different environments showed that FR901277 has a highly rigid bicyclic framework; however, it can slightly change its conformation according to the circumstances. The binding mode of FR901277 at the active site of PPE was directly applicable to that in HLE, after consideration of induced fit. The structure of the PPE-FR901277 complex provided much information regarding potential sites for modification of the physicochemical properties of FR901277.     

The solution structure of echistatin: evidence for disulphide bond rearrangement in homologous snake toxins.

The solution structure of the fibrinogen antagonist, echistatin, has been determined by a combination of NMR and simulated annealing methods. While the structure of the disulphide-linked core is well-defined by the NMR data, the N- and C-termini and the loop bearing the RGD sequence (which is responsible for the fibrinogen antagonist properties) are poorly defined. The pattern of disulphide bridges, which could not be determined by classical methods, was predicted by a statistical analysis of the simulated annealing structures. This pattern is distinct from that for the homologous protein kistrin, leading to the novel suggestion that homologous proteins possess non-conserved patterns of disulphide bridges.     

Three-dimensional structure of a complement control protein module in solution.

The complement control protein (CCP) modules (also known as short consensus repeats) are defined by a consensus sequence within a stretch of about 60 amino acid residues. These modules have been identified more than 140 times in over 20 proteins, including 12 proteins of the complement system. The solution structure of the 16th CCP module from human complement factor H has been determined by a combination of 2-dimensional nuclear magnetic resonance spectroscopy and restrained simulated annealing. In all, 548 structurally important nuclear Overhauser enhancement cross-peaks were quantified as distance restraints and, together with 41 experimentally measured angle restraints, were incorporated into a simulated annealing protocol to determine a family of closely related structures that satisfied the experimental observations. The CCP structure is shown to be based on a beta-sandwich arrangement; one face made up of three beta-strands hydrogen-bonded to form a triple-stranded region at its centre and the other face formed from two separate beta-strands. Both faces of the molecule contribute highly conserved hydrophobic side-chains to a compact core. The regions between the beta-strands are composed of both well-defined turns and less well-defined loops. Analysis of CCP sequence alignments, in light of the determined structure, reveals a high degree of conservation amongst residues of obvious structural importance, while almost all insertions, deletions or replacements observed in the known sequences are found in the less well-defined loop regions. On the basis of these observations it is postulated that models of other CCP modules that are based on the structure presented here will be accurate. Certain families of CCP modules differ from the consensus in that they contain extra cysteine residues. As a test of structural consensus, the extra disulphide bridges are shown to be easily accommodated within the determined CCP model.     

Solution structure of alpha t alpha, a helical hairpin peptide of de novo design.

alpha t alpha is a 38-residue peptide designed to adopt a helical hairpin conformation in solution (Fezoui Y, Weaver DL Osterhout JJ, 1995, Protein Sci 4:286-295). A previous study of the carboxylate form of alpha t alpha by CD and two-dimensional NMR indicated that the peptide was highly helical and that the helices associated in approximately the intended orientation (Fezoui Y, Weaver DL, Osterhout JJ, 1994, Proc Natl Acad Sci USA 91:3675-3679). Here, the solution structure of alpha t alpha as determined by two-dimensional NMR is reported. A total of 266 experimentally derived distance restraints and 20 dihedral angle restraints derived from J-couplings were used. One-hundred initial structures were generated by distance geometry and refined by dynamical simulated annealing. Twenty-three of the lowest-energy structures consistent with the experimental restraints were analyzed. The results presented here show that alpha t alpha is comprised of two associating helices connected by a turn region.     

Three-dimensional solution structure of the E3-binding domain of the dihydrolipoamide succinyltransferase core from the 2-oxoglutarate dehydrogenase multienzyme complex of Escherichia coli.

The three-dimensional solution structure of a 51-residue synthetic peptide comprising the dihydrolipoamide dehydrogenase (E3)-binding domain of the dihydrolipoamide succinyltransferase (E2) core of the 2-oxoglutarate dehydrogenase multienzyme complex of Escherichia coli has been determined by nuclear magnetic resonance spectroscopy and hybrid distance geometry-dynamical simulated annealing calculations. The structure is based on 630 approximate interproton distance and 101 torsion angle (phi, psi, chi 1) restraints. A total of 56 simulated annealing structures were calculated, and the atomic rms distribution about the mean coordinate positions for residues 12-48 of the synthetic peptide is 1.24 A for the backbone atoms, 1.68 A for all atoms, and 1.33 A for all atoms excluding the six side chains which are disordered at chi 1 and the seven which are disordered at chi 2; when the irregular partially disordered loop from residues 31 to 39 is excluded, the rms distribution drops to 0.77 A for the backbone atoms, 1.55 A for all atoms, and 0.89 A for ordered side chains. Although proton resonance assignments for the N-terminal 11 residues and the C-terminal 3 residues were obtained, these two segments of the polypeptide are disordered in solution as evidenced by the absence of nonsequential nuclear Overhauser effects. The solution structure of the E3-binding domain consists of two parallel helices (residues 14-23 and 40-48), a short extended strand (24-26), a five-residue helical-like turn, and an irregular (and more disordered) loop (residues 31-39). This report presents the first structure of an E3-binding domain from a 2-oxo acid dehydrogenase complex.(ABSTRACT TRUNCATED AT 250 WORDS)     

Design of lambda Cro fold: solution structure of a monomeric variant of the de novo protein

One of the classical DNA-binding proteins, bacteriophage lambda Cro, forms a homodimer with a unique fold of alpha-helices and beta-sheets. We have computationally designed an artificial sequence of 60 amino acid residues to stabilize the backbone tertiary structure of the lambda Cro dimer by simulated annealing using knowledge-based structure-sequence compatibility functions. The designed amino acid sequence has 25% identity with that of natural lambda Cro and preserves Phe58, which is important for formation of the stably folded structure of lambda Cro. The designed dimer protein and its monomeric variant, which was redesigned by the insertion of a beta-hairpin sequence at the C-terminal region to prevent dimerization, were synthesized and biochemically characterized to be well folded. The designed protein was monomeric under a wide range of protein concentrations and its solution structure was determined by NMR spectroscopy. The solved structure is similar to that of a monomeric variant of natural lambda Cro with a root-mean-square deviation of the polypeptide backbones at 2.1A and has a well-packed protein core. Thus, our knowledge-based functions provide approximate but essential relationships between amino acid sequences and protein structures, and are useful for finding novel sequences that are foldable into a given target structure.     

Optimal molecular structures of prion AGAAAAGA amyloid fibrils formatted by simulated annealing

To date, there is little structural data available on the AGAAAAGA palindrome in the hydrophobic region (113–120) of prion proteins, although many experimental studies have shown that this region has amyloid fibril forming properties. This region belongs to the N-terminal unstructured region (1–123) of prions, the structure of which has proved hard to determine using NMR or X-ray crystallography. This paper reports the successful construction of three amyloid fibril models for this region. The models were formatted by standard simulated annealing using suitable templates from the Protein Data Bank, and were refined using several traditional optimization methods within AMBER. Because the NMR or X-ray structure of the hydrophobic region AGAAAAGA of prion proteins has not yet been determined, these models can be used as a reference for experimental studies on this region. The results presented here confirm standard simulated annealing as an effective tool in molecular modeling. The three constructed models for amyloid fibrils may be useful in furthering the goals of medicinal chemistry in this field.     

NMR studies on turn mimetic analogs derived from melanocyte-stimulating hormones

Oligomers with alpha-aminooxy acids are reported to form very stable turn and helix structures, and they are supposed to be useful peptidomimetics for drug design. A recent report suggested that homochiral oxa-peptides form a strong eight-member-ring structure by a hydrogen bond between adjacent aminooxy-acid residues in a CDCl3 solution. In order to design an alpha-MSH analog with a stable turn conformation, we synthesized four tetramers and one pentamer, based on alpha-MSH sequence, and determined the solution structures of the molecules by two-dimensional NMR spectroscopy and simulated annealing calculations. The solution conformations of the three peptidomimetic molecules (TLV, TDV, and TLL) in DMSO-d6 contain a stable 7-membered-ring structure that is similar to a gamma-turn in normal peptides. Newly-designed tetramer TDF and pentamer PDF have a ball-type rigid structure that is induced by strong hydrogen bonds between adjacent amide protons and carbonyl oxygens. In conclusion, the aminooxy acids, easily prepared from natural or unnatural amino acids, can be employed to prepare peptidomimetic analogues with well-defined turn structures for pharmaceutical interest.     

Trajectory analysis of NMR structure calculations

Summary NMR as well as X-ray crystallography are used to determine the three-dimensional structures of macromolecules at atomic resolution. Structure calculation generates coordinates that are compatible with NMR data from randomly generated initial structures. We analyzed the trajectory taken by structures during NMR structure calculation in conformational space, assuming that the distance between two structures in conformational space is the root-mean-square deviation between the two structures. The coordinates of a structure in conformational space were obtained by applying the metric multidimensional scaling method. As an example, we used a 22-residue peptide, -Conotoxin GIIIA, and a simulated annealing protocol of XPLOR. We found that the three-dimensional solution of the multidimensional scaling analysis is sufficient to describe the overall configuration of the trajectories in conformational space. By comparing the trajectories of the entire calculation with those of the converged calculation, random sampling of conformational space is readily discernible. Trajectory analysis can also be used for optimization of protocols of NMR structure calculation, by examining individual trajectories.Abbreviations MD molecular dynamics - MDS multidimensional scaling - rmsd root-mean-square deviation - armsd angular rmsd - R multiple correlation coefficient - YASAP yet another simulated annealing protocol - PCA principal component analysis     

Crystal and molecular structure of RNase Rh, a new class of microbial ribonuclease from Rhizopus niveus.

The crystal structure of RNase Rh, a new class of microbial ribonuclease from Rhizopus niveus, has been determined at 2.5 A resolution by the multiple isomorphous replacement method. The crystal structure was refined by simulated annealing with molecular dynamics. The current crystallographic R-factor is 0.200 in the 10-2.5 A resolution range. The molecular structure which is completely different from the known structures of RNase A and RNase T1 consists of six alpha-helices and seven beta-strands, belonging to the alpha+beta type structure. Two histidine and one glutamic acid residues which were predicted as the most probably functional residues by chemical modification studies are found to be clustered. The steric nature of the active site taken together with the relevant site-directed mutagenesis experiments (Irie et al.) indicates that: (i) the two histidine residues are the general acid and base; and (ii) an aspartic acid residue plays a role of recognizing adenine moiety of the substrate.     

Determination of the complete three-dimensional structure of the trypsin inhibitor from squash seeds in aqueous solution by nuclear magnetic resonance and a combination of distance geometry and dynamical simulated annealing

The complete three-dimensional structure of the trypsin inhibitor from seeds of the squash Cucurbita maxima in aqueous solution was determined on the basis of 324 interproton distance constraints, 80 non-nuclear Overhauser effect distances, and 22 hydrogen-bonding constraints, supplemented by 27 phi backbone angle constraints derived from nuclear magnetic resonance measurements. The nuclear magnetic resonance input data were converted to the distance constraints in a semiquantitative manner after a sequence specific assignment of 1H spectra was obtained using two-dimensional nuclear magnetic resonance techniques. Stereospecific assignments were obtained for 17 of the 48 prochiral centers of the squash trypsin inhibitor using the floating chirality assignment introduced at the dynamical simulated annealing stage of the calculations. A total of 34 structures calculated by a hybrid distance geometry-dynamical simulated annealing method exhibit well-defined positions for both backbone and side-chain atoms. The average atomic root-mean-square difference between the individual structures and the minimized mean structure is 0.35(+/- 0.08) A for the backbone atoms and 0.89(+/- 0.17) A for all heavy atoms. The precision of the structure determination is discussed and correlated to the experimental input data.     

Solid phase synthesis and NMR conformational studies on cyclic decapeptide template molecule

Selectively addressable topological templates represent a key feature in the de novo design of proteins using the TASP concept (Template Assembled Synthetic Proteins). The regioselectively addressable (orthogonally protected) lysine-containing cyclic templates are especially interesting for combinatorial chemistry. We report the synthesis and structural analysis of a series of cyclic and bicyclic decapeptide templates (model of TASP molecules). The peptides were synthesized via solid phase synthesis and followed by solution cyclization. The conformation of the peptides was studied by proton nmr spectroscopy in dimethylsulfoxide and in trifluoroethanol/water. The structure of the peptide template was calculated with the program DIANA and followed by simulated annealing from the nmr experimental constraints. The peptides adopts a fold comprising two beta-strands and two type II beta-turns predicted for conformationally well-defined templates. The design of such a restained cyclic decapeptide template will be discussed along with Regioselectively Addressable Functionalized Template (RAFT) molecular recognition and template for combinatorial synthesis.     

Automated method for modeling seven-helix transmembrane receptors from experimental data.

A rule-based automated method is presented for modeling the structures of the seven transmembrane helices of G-protein-coupled receptors. The structures are generated by using a simulated annealing Monte Carlo procedure that positions and orients rigid helices to satisfy structural restraints. The restraints are derived from analysis of experimental information from biophysical studies on native and mutant proteins, from analysis of the sequences of related proteins, and from theoretical considerations of protein structure. Calculations are presented for two systems. The method was validated through calculations using appropriate experimental information for bacteriorhodopsin, which produced a model structure with a root mean square (rms) deviation of 1.87 A from the structure determined by electron microscopy. Calculations are also presented using experimental and theoretical information available for bovine rhodopsin to assign the helices to a projection density map and to produce a model of bovine rhodopsin that can be used as a template for modeling other G-protein-coupled receptors.     

Low resolution structure of interleukin-1 beta in solution derived from 1H-15N heteronuclear three-dimensional nuclear magnetic resonance spectroscopy

A low resolution solution structure of the cytokine interleukin-1 beta, a 153 residue protein of molecular weight 17,400, has been determined on the basis of 446 nuclear Overhauser effect (NOE) derived approximate interproton distance restraints involving solely NH, C alpha H and C beta H protons, supplemented by 90 distance restraints for 45 hydrogen bonds, and 79 phi torsion angle restraints. With the exception of 27 C alpha H-C alpha H NOEs, all the NOEs were assigned from a three-dimensional 1H-1H NOE 15N-1H heteronuclear multiple quantum coherence (HMQC) spectrum. The torsion angle restraints were obtained from accurate 3JHN alpha coupling constants measured from a HMQC-J spectrum, while the hydrogen bonds were derived from a qualitative analysis of the NOE, coupling constant and amide exchange data. A total of 20 simulated annealing (SA) structures was computed using the hybrid distance geometry-dynamical simulated annealing method. The solution structure of IL-1 beta comprises 12 beta-strands arranged in three pseudo-symmetrical topological units (each consisting of 5 anti-parallel beta-strands), joined by turns, short loops and long loops. The core of the structure, which is made up of the 12 beta-strands, together with the turns joining strands I and II, strands VIII and IX and strands X and XI, is well determined with a backbone atomic root-mean-square (r.m.s.) distribution about the mean co-ordinate positions of 1.2(+/- 0.1) A. The loop conformations, on the other hand, are poorly determined by the current data. A comparison of the core of the low resolution solution structure of IL-1 beta with that of the X-ray structure indicates that they are similar, with a backbone atomic r.m.s. difference of only 1.5 A between the co-ordinates of the restrained minimized mean of the SA structures and the X-ray structure.     

Solution structure of endothelin-3 determined using NMR spectroscopy.

The aqueous solution structure of the 21-residue vasoactive peptide hormone endothelin-3 has been determined using high-resolution NMR spectroscopy. A total of 177 proton-proton distance measurements and 5 chi 1 dihedral angle constraints derived from NMR spectra were used to calculate the structure using a combination of distance geometry and dynamical simulated annealing calculations. The calculations reveal a highly ordered, compact conformation in which a helical region extending from K9 to C15 lies in close apposition with the C-terminal hexapeptide; this interaction seems to be largely driven by hydrophobic interactions. Structure-activity studies are interpreted in terms of the conformational features of the calculated endothelin-3 structure.     

NMR three-dimensional solution structure of the serine protease inhibitor cyclotheonamide A

The nmr solution conformation of cyclotheonamide A (CtA) was determined in aqueous media. The data produced 15 distance and 10 torsional constraints which were used to generate conformations using restrained simulated annealing (SA) and distance geometry/simulated annealing (DG/SA) calculations. Two different calculation protocols were performed to ensure proper sampling of conformational space and even though the torsional restraints were input differently, both calculation methods yielded the same conformation of CtA. In the structure calculations, all solutions of the Karplus equation were sampled simultaneously using the restrained SA protocol and large ranges were used for the dihedral restraints in the DG/SA protocol because all solutions to the Karplus equation could not be sampled simultaneously. The solution conformation was also compared to the solid state x-ray conformations of CtA bound to thrombin and trypsin. The conformation of the residues important for active site binding (d -Phe, h-Arg, and Pro) are nearly identical in aqueous solution and solid state with largest differences at the a-Ala and v-Tyr residues. CtA appears to be preordered in structure and does not undergo a significant conformational change upon binding to the enzyme active site. © 1997 John Wiley & Sons, Inc.     

Three-dimensional structure of soybean trypsin/chymotrypsin Bowman-Birk inhibitor in solution.

The three-dimensional structure of soybean trypsin/chymotrypsin Bowman-Birk inhibitor in solution has been determined by two-dimensional 1H nuclear magnetic resonance spectroscopy and dynamical simulated annealing using the program XPLOR. The structure was defined by 907 NOEs involving intra- and interresidue contacts which served as distance constraints for a protocol of dynamical simulated annealing. In addition, 48 phi angle constraints involving non-proline amino acids, 29 chi angle constraints, six omega angle constraints for the X-Pro peptide bond, and 35 stereoassignments for prochiral centers were incorporated during the course of the calculation. The protein is characterized by two distinct binding domains for serine protease. Each domain is comprised of a beta-hairpin (antiparallel beta-sheet and a cis-proline-containing type VIb reverse turn) with a short segment making a third strand of antiparallel beta-sheet. The structure determination and refinement are described, and the structure is compared to other structures of Bowman-Birk inhibitors as well as other families of serine protease inhibitors.     

Solution structure of a sweet protein single-chain monellin determined by nuclear magnetic resonance and dynamical simulated annealing calculations

Single-chain monellin (SCM), which is an engineered 94-residue polypeptide, has proven to be as sweet as native two-chain monellin. SCM is more stable than the native monellin for both heat and acidic environments. Data from gel filtration HPLC and NMR indicate that the SCM exists as a monomer in aqueous solution. The solution structure of SCM has been determined by nuclear magnetic resonance (NMR) spectroscopy and dynamical simulated annealing calculations. A stable alpha-helix spanning residues Phe11-Ile26 and an antiparallel beta-sheet formed by residues 2-5, 36-38, 41-47, 54-64, 69-75, and 83-88 have been identified. The sheet was well defined by backbone-backbone NOEs, and the corresponding beta-strands were further confirmed by hydrogen bond networks based on amide hydrogen exchange data. Strands beta2 and beta3 are connected by a small bulge comprising residues Ile38-Cys41. A total of 993 distance and 56 dihedral angle restraints were used for simulated annealing calculations. The final simulated annealing structures (k) converged well with a root-mean-square deviation (rmsd) between backbone atoms of 0.49 A for secondary structural regions and 0.70 A for backbone atoms excluding two loop regions. The average restraint energy-minimized (REM) structure exhibited root-mean-square deviations of 1.19 A for backbone atoms and 0.85 A for backbone atoms excluding two loop regions with respect to 20 k structures. The solution structure of SCM revealed that the long alpha-helix was folded into the concave side of a six-stranded antiparallel beta-sheet. The side chains of Tyr63 and Asp66 which are common to all sweet peptides showed an opposite orientation relative to H1 helix, and they were all solvent-exposed. Residues at the proposed dimeric interface in the X-ray structure were observed to be mostly solvent-exposed and demonstrated high degrees of flexibility.