Estimation of the number of alpha-helical and beta-strand segments in proteins using circular dichroism spectroscopy

A simple approach to estimate the number of alpha-helical and beta-strand segments from protein circular dichroism spectra is described. The alpha-helix and beta-sheet conformations in globular protein structures, assigned by DSSP and STRIDE algorithms, were divided into regular and distorted fractions by considering a certain number of terminal residues in a given alpha-helix or beta-strand segment to be distorted. The resulting secondary structure fractions for 29 reference proteins were used in the analyses of circular dichroism spectra by the SELCON method. From the performance indices of the analyses, we determined that, on an average, four residues per alpha-helix and two residues per beta-strand may be considered distorted in proteins. The number of alpha-helical and beta-strand segments and their average length in a given protein were estimated from the fraction of distorted alpha-helix and beta-strand conformations determined from the analysis of circular dichroism spectra. The statistical test for the reference protein set shows the high reliability of such a classification of protein secondary structure. The method was used to analyze the circular dichroism spectra of four additional proteins and the predicted structural characteristics agree with the crystal structure data.     

A model for the aggregation of the acylphosphatase from Sulfolobus solfataricus in its native-like state

Evidence is accumulating that normally folded proteins retain a significant tendency to form amyloid fibrils through a direct assembly of monomers in their native-like conformation. However, the factors promoting such processes are not yet well understood. The acylphosphatase from Sulfolobus solfataricus (Sso AcP) aggregates under conditions in which a native-like state is initially populated and forms, as a first step, aggregates in which the monomers maintain their native-like topology. An unstructured N-terminal segment and an edge beta-strand were previously shown to play a major role in the process. Using kinetic experiments on a set of Sso AcP variants we shall show that the major event of the first step is the establishment of an inter-molecular interaction between the unstructured segment of one Sso AcP molecule and the globular unit of another molecule. This interaction is determined by the primary sequence of the unstructured segment and not by its physico-chemical properties. Moreover, we shall show that the conversion of these initial aggregates into amyloid-like protofibrils is an intra-molecular process in which the Sso AcP molecules undergo conformational modifications. The obtained results allow the formulation of a model for the assembly of Sso AcP into amyloid-like aggregates at a molecular level.     

Side-chain entropy effects on protein secondary structure formation

Loss of conformational entropy is one of the primary factors opposing protein folding. Both the backbone and side-chain of each residue in a protein will have their freedom of motion restricted in the final folded structure. The type of secondary structure of which a residue is part will have a significant impact on how much side-chain entropy is lost. Side-chain conformational entropies have previously been determined for folded proteins, simple models of unfolded proteins, alpha-helices, and a dipeptide model for beta-strands, but not for polyproline II (PII) helices. In this work, we present side-chain conformational estimates for the three regular secondary structure types: alpha-helices, beta-strands, and PII helices. Entropies are estimated from Monte Carlo computer simulations. Beta-strands are modeled as two structures, parallel and antiparallel beta-strands. Our data indicate that restraining a residue to the PII helix or antiparallel beta-strand conformations results in side-chain entropies equal to or higher than those obtained by restraining residues to the parallel beta-strand conformation. Side-chains in the alpha-helix conformation have the lowest side-chain entropies. The observation that extended structures retain the most side-chain entropy suggests that such structures would be entropically favored in unfolded proteins under folding conditions. Our data indicate that the PII helix conformation would be somewhat favored over beta-strand conformations, with antiparallel beta-strand favored over parallel. Notably, our data imply that, under some circumstances, residues may gain side-chain entropy upon folding. Implications of our findings for protein folding and unfolded states are discussed.     

Physical-chemical determinants of turn conformations in globular proteins

Globular proteins are assemblies of alpha-helices and beta-strands, interconnected by reverse turns and longer loops. Most short turns can be classified readily into a limited repertoire of discrete backbone conformations, but the physical-chemical determinants of these distinct conformational basins remain an open question. We investigated this question by exhaustive analysis of all backbone conformations accessible to short chain segments bracketed by either an alpha-helix or a beta-strand (i.e., alpha-segment-alpha, beta-segment-beta, alpha-segment-beta, and beta-segment-alpha) in a nine-state model. We find that each of these four secondary structure environments imposes its own unique steric and hydrogen-bonding constraints on the intervening segment, resulting in a limited repertoire of conformations. In greater detail, an exhaustive set of conformations was generated for short backbone segments having reverse-turn chain topology and bracketed between elements of secondary structure. This set was filtered, and only clash-free, hydrogen-bond-satisfied conformers having reverse-turn topology were retained. The filtered set includes authentic turn conformations, observed in proteins of known structure, but little else. In particular, over 99% of the alternative conformations failed to satisfy at least one criterion and were excluded from the filtered set. Furthermore, almost all of the remaining alternative conformations have close tolerances that would be too tight to accommodate side chains longer than a single beta-carbon. These results provide a molecular explanation for the observation that reverse turns between elements of regular secondary can be classified into a small number of discrete conformations.     

Substrate-induced conformational changes of the periplasmic N-terminus of an outer-membrane transporter by site-directed spin labeling

The structure and dynamics of the N-terminal and core regions of BtuB, an outer membrane vitamin B(12) transporter from Escherichia coli, were investigated by site-directed spin labeling. Cysteine mutants were generated by site-directed mutagenesis to place spin labels in the N-terminal region (residues 1-17), the core region (residues 25-30), and double labels into the Ton box (residues 6-12). BtuB mutants were expressed, spin labeled, purified, and reconstituted into phosphatidylcholine. In the presence of substrate (vitamin B(12)), EPR spectroscopy demonstrates that there is a conformational change in the Ton box similar to that seen previously for BtuB in intact outer membranes. The Ton box is positioned within the beta-barrel of BtuB in the absence of substrate (docked configuration) but becomes unfolded and increases its aqueous exposure upon substrate binding (undocked configuration). This conformational change and the similarity in the EPR spectra between reconstituted and native membranes indicate that BtuB is correctly folded and functional in the reconstituted system. The protein segment on the N-terminal side of the Ton box is highly mobile, and it becomes more mobile in the presence of substrate. Side chains in the region C-terminal to the Ton box also show increases in mobility with substrate addition, but position 16 appears to define a hinge point for this conformation change. EPR line shapes and relaxation data indicate that residues 25-30 form a beta-strand structure, which is analogous to the first beta-strand in the cores of the homologous iron transporters. When substrate binds to BtuB, this first beta-strand remains folded. The EPR spectra of double-nitroxide labels within the Ton box are broadened because of dipolar and collisional exchange interactions. The broadening pattern indicates that the Ton box is not helical but is in an extended or beta-strand structure.     

Solution structure of Escherichia coli Par10: The prototypic member of the Parvulin family of peptidyl-prolyl cis/trans isomerases

E. coli Par10 is a peptidyl-prolyl cis/trans isomerase (PPIase) from Escherichia coli catalyzing the isomerization of Xaa-Pro bonds in oligopeptides with a broad substrate specificity. The structure of E. coli Par10 has been determined by multidimensional solution-state NMR spectroscopy based on 1207 conformational constraints (1067 NOE-derived distances, 42 vicinal coupling-constant restraints, 30 hydrogen-bond restraints, and 68 phi/psi restraints derived from the Chemical Shift Index). Simulated-annealing calculations with the program ARIA and subsequent refinement with XPLOR yielded a set of 18 convergent structures with an average backbone RMSD from mean atomic coordinates of 0.50 A within the well-defined secondary structure elements. E. coli Par10 is the smallest known PPIase so far, with a high catalytic efficiency comparable to that of FKBPs and cyclophilins. The secondary structure of E. coli Par10 consists of four helical regions and a four-stranded antiparallel beta-sheet. The N terminus forms a beta-strand, followed by a large stretch comprising three alpha-helices. A loop region containing a short beta-strand separates these helices from a fourth alpha-helix. The C terminus consists of two more beta-strands completing the four-stranded anti-parallel beta-sheet with strand order 2143. Interestingly, the third beta-strand includes a Gly-Pro cis peptide bond. The curved beta-strand forms a hydrophobic binding pocket together with alpha-helix 4, which also contains a number of highly conserved residues. The three-dimensional structure of Par10 closely resembles that of the human proteins hPin1 and hPar14 and the plant protein Pin1At, belonging to the same family of highly homologous proteins.     

A conformational analysis of Walker motif A [GXXXXGKT (S)] in nucleotide-binding and other proteins

The sequence GXXXXGKT/S, popularly known as Walker motif A, is widely believed to be the site for binding nucleotides in many proteins. Examination of the crystal structures in the Protein Data Bank showed that about half of the examples having these sequences do not bind or use nucleotides. Data analyses showed 92 different Walker sequences of the variable quartet (XXXX). Ramachandran angles in this segment revealed conformational similarity in the group of 45 proteins, known to bind or utilize nucleotides. The conformations of this segment in other proteins differ widely and it is not known whether they play any role in their functions. A flip of a peptide unit at different locations, with little change in the backbone conformation was noted in nine pairs of these proteins having same Walker sequence. An examination of the immediate neighborhood of the Walker sequence indicates that this region is preceded by a beta-strand and followed by an alpha-helix, resulting in the motif beta-W-alpha, an invariant feature amongst nucleotide-binding proteins.     

Three-dimensional structure of the rSly1 N-terminal domain reveals a conformational change induced by binding to syntaxin 5

Sec1/Mun18-like (SM) proteins and soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) play central roles in intracellular membrane fusion. Diverse modes of interaction between SM proteins and SNAREs from the syntaxin family have been described. However, the observation that the N-terminal domains of Sly1 and Vps45, the SM proteins involved in traffic at the endoplasmic reticulum, the Golgi, the trans-Golgi network and the endosomes, bind to similar N-terminal sequences of their cognate syntaxins suggested a unifying theme for SM protein/SNARE interactions in most internal membrane compartments. To further understand this mechanism of SM protein/SNARE coupling, we have elucidated the structure in solution of the isolated N-terminal domain of rat Sly1 (rSly1N) and analyzed its complex with an N-terminal peptide of rat syntaxin 5 by NMR spectroscopy. Comparison with the crystal structure of a complex between Sly1p and Sed5p, their yeast homologues, shows that syntaxin 5 binding requires a striking conformational change involving a two-residue shift in the register of the C-terminal beta-strand of rSly1N. This conformational change is likely to induce a significant alteration in the overall shape of full-length rSly1 and may be critical for its function. Sequence analyses indicate that this conformational change is conserved in the Sly1 family but not in other SM proteins, and that the four families represented by the four SM proteins found in yeast (Sec1p, Sly1p, Vps45p and Vps33p) diverged early in evolution. These results suggest that there are marked distinctions between the mechanisms of action of each of the four families of SM proteins, which may have arisen from different regulatory requirements of traffic in their corresponding membrane compartments.     

Protein folding pathways and kinetics: molecular dynamics simulations of beta-strand motifs

The folding pathways and the kinetic properties for three different types of off-lattice four-strand antiparallel beta-strand protein models interacting via a hybrid Go-type potential have been investigated using discontinuous molecular dynamics simulations. The kinetic study of protein folding was conducted by temperature quenching from a denatured or random coil state to a native state. The progress parameters used in the kinetic study include the squared radius of gyration R(2)(g), the fraction of native contacts within the protein as a whole Q, and between specific strands Q(ab). In the time series of folding, the denatured proteins undergo a conformational change toward the native state. The model proteins exhibit a variety of kinetic folding pathways that include a fast-track folding pathway without passing through an intermediate and multiple pathways with trapping into more than one intermediate. The kinetic folding behavior of the beta-strand proteins strongly depends on the native-state geometry of the model proteins and the size of the bias gap g, an artificial measure of a model protein's preference for its native state.     

Retrospective analysis of a secondary structure prediction: the catalytic domain of matrix metalloproteinases.

Secondary structure prediction of the catalytic domain of matrix metalloproteinases is evaluated in the light of recently published experimentally determined structures. The prediction was made by combining conformational propensity, surface probability, and residue conservation calculated for an alignment of 19 sequences. The position of each observed secondary structure element was correctly predicted with a high degree of accuracy, with a single beta-strand falsely predicted. The domain fold was also anticipated from the prediction by analogy with the structural elements found in the distantly related metalloproteinases thermolysin, astacin, and adamalysin.     

Conversion of a beta-strand to an alpha-helix induced by a single-site mutation observed in the crystal structure of Fis mutant Pro26Ala.

The conversion from an alpha-helix to a beta-strand has received extensive attention since this structural change may induce many amyloidogenic proteins to self-assemble into fibrils and cause fatal diseases. Here we report the conversion of a peptide segment from a beta-strand to an alpha-helix by a single-site mutation as observed in the crystal structure of Fis mutant Pro26Ala determined at 2.0 A resolution. Pro26 in Fis occurs at the point where a flexible extended beta-hairpin arm leaves the core structure. Thus it can be classified as a "hinge proline" located at the C-terminal end of the beta2-strand and the N-terminal cap of the A alpha-helix. The replacement of Pro26 to alanine extends the A alpha-helix for two additional turns in one of the dimeric subunits; therefore, the structure of the peptide from residues 22 to 26 is converted from a beta-strand to an alpha-helix. This result confirms the structural importance of the proline residue located at the hinge region and may explain the mutant''s reduced ability to activate Hin-catalyzed DNA inversion. The peptide (residues 20 to 26) in the second monomer subunit presumably retains its beta-strand conformation in the crystal; therefore, this peptide shows a "chameleon-like" character since it can adopt either an alpha-helix or a beta-strand structure in different environments. The structure of Pro26Ala provides an additional example where not only the protein sequence, but also non-local interactions determine the secondary structure of proteins.     

The importance of surface loops for stabilizing an eightfold beta alpha barrel protein.

An important step in understanding how a protein folds is to determine those regions of the sequence that are critical to both its stability and its folding pathway. We chose phosphoribosyl anthranilate isomerase from Escherichia coli, which is a monomeric representative of the (beta alpha)8 barrel family of proteins, to construct a variant that carries an internal tandem duplication of the fifth beta alpha module. This (beta alpha)9 variant was enzymically active and therefore must have a wild-type (beta alpha)8 core. It had a choice a priori to fold to three different folding frames, which are distinguished by carrying the duplicated segment as an insert into one out of three different loops. Steady-state kinetic constants, the fluorescence properties of a crucial tryptophan residue, and limited proteolysis showed that the stable (beta alpha)9 variant carries the insertion between beta-strand 5 and alpha-helix 5. This preference can be explained by the important role of loops between alpha helices and beta strands in stabilizing the structure of the enzyme.     

Prediction of polyelectrolyte polypeptide structures using Monte Carlo conformational search methods with implicit solvation modeling.

Many interesting proteins possess defined sequence stretches containing negatively charged amino acids. At present, experimental methods (X-ray crystallography, NMR) have failed to provide structural data for many of these sequence domains. We have applied the dihedral probability grid-Monte Carlo (DPG-MC) conformational search algorithm to a series of N- and C-capped polyelectrolyte peptides, (Glu)20, (Asp)20, (PSer)20, and (PSer-Asp)10, that represent polyanionic regions in a number of important proteins, such as parathymosin, calsequestrin, the sodium channel protein, and the acidic biomineralization proteins. The atomic charges were estimated from charge equilibration and the valence and van der Waals parameters are from DREIDING. Solvation of the carboxylate and phosphate groups was treated using sodium counterions for each charged side chain (one Na+ for COO-; two Na for CO(PO3)-2) plus a distance-dependent (shielded) dielectric constant, epsilon = epsilon 0 R, to simulate solvent water. The structures of these polyelectrolyte polypeptides were obtained by the DPG-MC conformational search with epsilon 0 = 10, followed by calculation of solvation energies for the lowest energy conformers using the protein dipole-Langevin dipole method of Warshel. These calculations predict a correlation between amino acid sequence and global folded conformational minima: 1. Poly-L-Glu20, our structural benchmark, exhibited a preference for right-handed alpha-helix (47% helicity), which approximates experimental observations of 55-60% helicity in solution. 2. For Asp- and PSer-containing sequences, all conformers exhibited a low preference for right-handed alpha-helix formation (< or = 10%), but a significant percentage (approximately 20% or greater) of beta-strand and beta-turn dihedrals were found in all three sequence cases: (1) Aspn forms supercoil conformers, with a 2:1:1 ratio of beta-turn:beta-strand:alpha-helix dihedral angles; (2) PSer20 features a nearly 1:1 ratio of beta-turn:beta-sheet dihedral preferences, with very little preference for alpha-helical structure, and possesses short regions of strand and turn combinations that give rise to a collapsed bend or hairpin structure; (3) (PSer-Asp)10 features a 3:2:1 ratio of beta-sheet:beta-turn:alpha-helix and gives rise to a superturn or C-shaped structure.     

Homology modeling in a dynamical world

A key concept in template‐based modeling (TBM) is the high correlation between sequence and structural divergence, with the practical consequence that homologous proteins that are similar at the sequence level will also be similar at the structural level. However, conformational diversity of the native state will reduce the correlation between structural and sequence divergence, because structural variation can appear without sequence diversity. In this work, we explore the impact that conformational diversity has on the relationship between structural and sequence divergence. We find that the extent of conformational diversity can be as high as the maximum structural divergence among families. Also, as expected, conformational diversity impairs the well‐established correlation between sequence and structural divergence, which is nosier than previously suggested. However, we found that this noise can be resolved using a priori information coming from the structure‐function relationship. We show that protein families with low conformational diversity show a well‐correlated relationship between sequence and structural divergence, which is severely reduced in proteins with larger conformational diversity. This lack of correlation could impair TBM results in highly dynamical proteins. Finally, we also find that the presence of order/disorder can provide useful beforehand information for better TBM performance.     

Construction and characterization of protein libraries composed of secondary structure modules

Only a minute fraction of all possible protein sequences can exist in the genomes of all life forms. To explore whether physicochemical constraints or a lack of need causes the paucity of different protein folds, we set out to construct protein libraries without any restriction of topology. We generated different libraries (all alpha-helix, all beta-strand, and alpha-helix plus beta-strand) with an average length of 100 amino acid residues, composed of designed secondary structure modules (alpha-helix, beta-strand, and beta-turn) in various proportions, based primarily on the patterning of polar and nonpolar residues. We wished to explore that part of sequence space that is rich in secondary structure. The analysis of randomly chosen clones from each of the libraries showed that, despite the low sequence homology to known protein sequences, a substantial proportion of the library members containing alpha-helix modules were indeed helical, possess a defined oligomerization state, and showed cooperative chemical unfolding behavior. On the other hand, proteins composed of mainly beta-strand modules tended to form amyloid-like fibrils and were among the least soluble proteins ever reported. We found that a large fraction of members in non-beta-strand-containing protein libraries that are distant from natural proteins in sequence space possess unexpectedly favorable properties. These results reinforce the efficacy of applying binary patterning to design proteins with native-like properties despite lack of restriction in topology. Because of the intrinsic tendency of beta-strand modules to aggregate, their presence requires precise topologic arrangement to prevent fibril formation.     

Improved beta-protein structure prediction by multilevel optimization of nonlocal strand pairings and local backbone conformation

Proteins with complex, nonlocal beta-sheets are challenging for de novo structure prediction, due in part to the difficulty of efficiently sampling long-range strand pairings. We present a new, multilevel approach to beta-sheet structure prediction that circumvents this difficulty by reformulating structure generation in terms of a folding tree. Nonlocal connections in this tree allow us to explicitly sample alternative beta-strand pairings while simultaneously exploring local conformational space using backbone torsion-space moves. An iterative, energy-biased resampling strategy is used to explore the space of beta-strand pairings; we expect that such a strategy will be generally useful for searching large conformational spaces with a high degree of combinatorial complexity.     

A simple and fast approach to prediction of protein secondary structure from multiply aligned sequences with accuracy above 70%.

To improve secondary structure predictions in protein sequences, the information residing in multiple sequence alignments of substituted but structurally related proteins is exploited. A database comprised of 70 protein families and a total of 2,500 sequences, some of which were aligned by tertiary structural superpositions, was used to calculate residue exchange weight matrices within alpha-helical, beta-strand, and coil substructures, respectively. Secondary structure predictions were made based on the observed residue substitutions in local regions of the multiple alignments and the largest possible associated exchange weights in each of the three matrix types. Comparison of the observed and predicted secondary structure on a per-residue basis yielded a mean accuracy of 72.2%. Individual alpha-helix, beta-strand, and coil states were respectively predicted at 66.7, and 75.8% correctness, representing a well-balanced three-state prediction. The accuracy level, verified by cross-validation through jack-knife tests on all protein families, dropped, on average, to only 70.9%, indicating the rigor of the prediction procedure. On the basis of robustness, conceptual clarity, accuracy, and executable efficiency, the method has considerable advantage, especially with its sole reliance on amino acid substitutions within structurally related proteins.     

Structures, basins, and energies: a deconstruction of the Protein Coil Library

Globular proteins adopt complex folds, composed of organized assemblies of alpha-helix and beta-sheet together with irregular regions that interconnect these scaffold elements. Here, we seek to parse the irregular regions into their structural constituents and to rationalize their formative energetics. Toward this end, we dissected the Protein Coil Library, a structural database of protein segments that are neither alpha-helix nor beta-strand, extracted from high-resolution protein structures. The backbone dihedral angles of residues from coil library segments are distributed indiscriminately across the phi,psi map, but when contoured, seven distinct basins emerge clearly. The structures and energetics associated with the two least-studied basins are the primary focus of this article. Specifically, the structural motifs associated with these basins were characterized in detail and then assessed in simple simulations designed to capture their energetic determinants. It is found that conformational constraints imposed by excluded volume and hydrogen bonding are sufficient to reproduce the observed ,psi distributions of these motifs; no additional energy terms are required. These three motifs in conjunction with alpha-helices, strands of beta-sheet, canonical beta-turns, and polyproline II conformers comprise approximately 90% of all protein structure.     

Modulation of Polyglutamine Conformations and Dimer Formation by the N-Terminus of Huntingtin

Polyglutamine expansions within different proteins are associated with nine different neurodegenerative diseases. There is growing interest in understanding the roles of flanking sequences from disease-relevant proteins in the intrinsic conformational and aggregation properties of polyglutamine. We report results from atomistic simulations and circular dichroism experiments that quantify the effect of the N-terminal 17-residue (Nt17) segment of the huntingtin protein on polyglutamine conformations and intermolecular interactions. We show that the Nt17 segment and polyglutamine domains become increasingly disordered as polyglutamine length (N) increases in Nt17-Q_N constructs. Hydrophobic groups within Nt17 become sequestered in intramolecular interdomain interfaces. We also show that the Nt17 segment suppresses the intrinsic propensity of polyglutamine aggregation. This inhibition arises from the incipient micellar structures adopted by monomeric forms of the peptides with Nt17 segments. The degree of intermolecular association increases with increasing polyglutamine length and is governed mainly by associations between polyglutamine domains. Comparative analysis of intermolecular associations for different polyglutamine-containing constructs leads to clearer interpretations of recently published experimental data. Our results suggest a framework for fibril formation and identify roles for flanking sequences in the modulation of polyglutamine aggregation.     

Variability of the canonical loop conformations in serine proteinases inhibitors and other proteins

Canonical loops of protein inhibitors of serine proteinases occur in proteins having completely different folds. In this article, conformations of the loops have been analyzed for inhibitors belonging to 10 structurally different families. Using deviation in C_α-C_α distances as a criterion for loop similarity, we found that the P3-P3′ segment defines most properly the length of the loop. When conformational differences among loops of individual inhibitors were compared using root mean square deviation (rmsd) in atomic coordinates for all main chain atoms (Δr method) and rmsd operating in main chain torsion angles (Δt method), differences of up to 2.1 Å and 72.3°, respectively, were observed. Such large values indicate significant conformational differences among individual loops. Nevertheless, the overall geometry of the inhibitor-proteinase interaction is very well preserved, as judged from the similarity of C_α-C_α distances between C_α of catalytic Ser and C_α of P3-P3′ residues in various enzyme-inhibitor complexes. The mode of interaction is very well preserved both in the chymotrypsin and subtilisin families, as distances calculated for subtilisin-inhibitor complexes are almost always within the range of those for chymotrypsin-inhibitor complexes. Complex formation leads to conformational changes of up to 160° for χ₁ angle. Side chains of residue P1 and P2′ adopt in different complexes a similar orientation (χ₁ angle = −60° and −180°, respectively). To check whether the canonical conformation can be found among non–proteinase-inhibitor Brookhaven Protein Data Bank structures, two selection criteria—the allowed main chain dihedral angles and C_α-C_α distances for the P3-P3′ segment—were applied to all these structures. This procedure detected 132 unique hexapeptide segments in 121 structurally and functionally unrelated proteins. Partial preferences for certain amino acids occurring at particular positions in these hexapeptides could be noted. Further restriction of this set to hexapeptides with a highly exposed P1 residue side chain resulted in 40 segments. The possibility of complexes formation between these segments and serine proteinases was ruled out in molecular modeling due to steric clashes. Several structural features that determine the canonical conformation of the loop both in inhibitors and in other proteins can be distinguished. They include main chain hydrogen bonds both within the P3-P3′ segment and with the scaffold region, P3-P4 and P3′-P4′ hydrophobic interactions, and finally either hydrophobic or polar interactions involving the P1′ residue. Proteins 32:459–474, 1998. © 1998 Wiley-Liss, Inc.