We have determined the three-dimensional solution structure of the complex of Lactobacillus casei dihydrofolate reductase and the anticancer drug trimetrexate. Two thousand seventy distance, 345 dihedral angle, and 144 hydrogen bond restraints were obtained from analysis of multidimensional NMR spectra recorded for complexes containing 15N-labeled protein. Simulated annealing calculations produced a family of 22 structures fully consistent with the constraints. Several intermolecular protein-ligand NOEs were obtained by using a novel approach monitoring temperature effects of NOE signals resulting from dynamic processes in the bound ligand. At low temperature (5 °C) the trimethoxy ring of bound trimetrexate is flipping sufficiently slowly to give narrow signals in slow exchange, which give good NOE cross peaks. At higher temperature these broaden and their NOE cross peaks disappear thus allowing the signals in the lower-temperature spectrum to be identified as NOEs involving ligand protons. The binding site for trimetrexate is well defined and this was compared with the binding sites in related complexes formed with methotrexate and trimethoprim. No major conformational differences were detected between the different complexes. The 2,4-diaminopyrimidine-containing moieties in the three drugs bind essentially in the same binding pocket and the remaining parts of their molecules adapt their conformations such that they can make effective van der Waals interactions with essentially the same set of hydrophobic amino acids, the side-chain orientations and local conformations of which are not greatly changed in the different complexes (similar χ1 and χ2 values).

The high affinity interleukin-2 receptor is composed of three cell surface subunits, IL-2Rα, IL-2Rβ, and IL-2Rγ. Functional forms of the IL-2 receptor exist, however, that enlist only two of the three subunits. On activated T-cells, the α- and β-subunits combine as a preformed heterodimer (the pseudo-high affinity receptor) that serves to capture IL-2. On a subpopulation of natural killer cells, the β- and γ-subunits interact in a ligand-dependent manner to form the intermediate affinity receptor site. Previously, we have demonstrated the feasibility of employing coiled-coil molecular recognition for the solution assembly of a heteromeric IL-2 receptor complex. In that study, although the receptor was functional, the coiled-coil complex was a trimer rather than the desired heterodimer. We have now redesigned the hydrophobic heptad sequences of the coiled-coils to generate soluble forms of both the pseudo-high affinity and the intermediate affinity heterodimeric IL-2 receptors. The properties of these complexes were examined and their relevance to the physiological IL-2 receptor mechanism is discussed.

The heme-containing catalase HPII of Escherichia coli consists of a homotetramer in which each subunit contains a core region with the highly conserved catalase tertiary structure, to which are appended N- and C-terminal extensions making it the largest known catalase. HPII does not bind NADPH, a cofactor often found in catalases. In HPII, residues 585–590 of the C-terminal extension protrude into the pocket corresponding to the NADPH binding site in the bovine liver catalase. Despite this difference, residues that define the NADPH pocket in the bovine enzyme appear to be well preserved in HPII. Only two residues that interact ionically with NADPH in the bovine enzyme (Asp212 and His304) differ in HPII (Glu270 and Glu362), but their mutation to the bovine sequence did not promote nucleotide binding. The active-site heme groups are deeply buried inside the molecular structure requiring the movement of substrate and products through long channels. One potential channel is about 30 Å in length, approaches the heme active site laterally, and is structurally related to the branched channel associated with the NADPH binding pocket in catalases that bind the dinucleotide. In HPII, the upper branch of this channel is interrupted by the presence of Arg260 ionically bound to Glu270. When Arg260 is replaced by alanine, there is a threefold increase in the catalytic activity of the enzyme. Inhibitors of HPII, including azide, cyanide, various sulfhydryl reagents, and alkylhydroxylamine derivatives, are effective at lower concentration on the Ala260 mutant enzyme compared to the wild-type enzyme. The crystal structure of the Ala260 mutant variant of HPII, determined at 2.3 Å resolution, revealed a number of local structural changes resulting in the opening of a second branch in the lateral channel, which appears to be used by inhibitors for access to the active site, either as an inlet channel for substrate or an exhaust channel for reaction products.

This paper reports on an insulin analogue with 12.5-fold receptor affinity, the highest increase observed for a single replacement, and on its solution structure, determined by NMR spectroscopy.

The analogue is [d-AlaB26]des-(B27–B30)-tetrapeptide-insulin-B26-amide. C-terminal truncation of the B-chain by four (or five) residues is known not to affect the functional properties of insulin, provided the new carboxylate charge is neutralized. As opposed to the dramatic increase in receptor affinity caused by the substitution of d-Ala for the wild-type residue TyrB26 in the truncated molecule, this very substitution reduces it to only 18% of that of the wild-type hormone when the B-chain is present in full length.

The insulin molecule in solution is visualized as an ensemble of conformers interrelated by a dynamic equilibrium. The question is whether the “active” conformation of the hormone, sought after in innumerable structure/function studies, is or is not included in the accessible conformational space, so that it could be adopted also in the absence of the receptor. If there were any chance for the active conformation, or at least a predisposed state to be populated to a detectable extent, this chance should be best in the case of a superpotent analogue. This was the motivation for the determination of the three-dimensional structure of [d-AlaB26]des-(B27–B30)-tetrapeptide-insulin-B26-amide. However, neither the NMR data nor CD spectroscopic comparison of a number of related analogues provided a clue concerning structural features predisposing insulin to high receptor affinity. After the present study it seems more likely than before that insulin will adopt its active conformation only when exposed to the force field of the receptor surface.

The most common form of systemic amyloidosis originates from antibody light chains. The large number of amino acid variations that distinguish amyloidogenic from nonamyloidogenic light chain proteins has impeded our understanding of the structural basis of light-chain fibril formation. Moreover, even among the subset of human light chains that are amyloidogenic, many primary structure differences are found. We compared the thermodynamic stabilities of two recombinant κ4 light-chain variable domains (VLs) derived from amyloidogenic light chains with a VL from a benign light chain. The amyloidogenic VLs were significantly less stable than the benign VL. Furthermore, only the amyloidogenic VLs formed fibrils under native conditions in an in vitro fibril formation assay. We used site-directed mutagenesis to examine the consequences of individual amino acid substitutions found in the amyloidogenic VLs on stability and fibril formation capability. Both stabilizing and destabilizing mutations were found; however, only destabilizing mutations induced fibril formation in vitro. We found that fibril formation by the benign VL could be induced by low concentrations of a denaturant. This indicates that there are no structural or sequence-specific features of the benign VL that are incompatible with fibril formation, other than its greater stability. These studies demonstrate that the VL β-domain structure is vulnerable to destabilizing mutations at a number of sites, including complementarity determining regions (CDRs), and that loss of variable domain stability is a major driving force in fibril formation.

CyclicAMP receptor protein (CRP) regulates transcription of numerous genes in Escherichia coli. Both cAMP and cGMP bind CRP, but only cAMP induces conformational changes that dramatically increase the specific DNA binding activity of the protein. We have shown previously that our protein footprinting technique is sensitive enough to detect conformational changes in CRP by cAMP [Baichoo N, Heyduk T. 1997. Biochemistry 36:10830–10836]. In this work, conformational changes in CRP induced by cAMP and cGMP binding were mapped and quantitatively analyzed by protein footprinting using iron complexed to diethylenetriaminepentaacetic acid ([Fe-DTPA]2−), iron complexed to ethylenediaminediacetic acid ([Fe-EDDA]), iron complexed to desferrioxamine mesylate ([Fe-HDFO]+), and copper complexed to o-phenanthroline ([(OP)2Cu]+) as proteases. These chemical proteases differ in size, charge, and hydrophobicity. Binding of cAMP to CRP resulted in changes in susceptibility to cleavage by all four proteases. Cleavage by [Fe-EDDA] and [Fe-DTPA]2− of CRP-cAMP detected hypersensitivities in the DNA-binding F α-helix, the interdomain hinge, and the ends of the C α-helix, which is involved in intersubunit interactions. [Fe-EDDA] and [Fe-DTPA]2− also detected reductions in cleavage in the D and E α-helices, which are involved in DNA recognition. Cleavage by [Fe-HDFO]+ of CRP-cAMP detected hypersensitivities in β-strand 8, the B α-helix, as well as in parts of the F and C α-helices. [Fe-HDFO]+ also detected protections from cleavage in β-strands 4 to 5 and their intervening loop, β-strand 7, which is part of the nucleotide binding pocket, as well as in the D and E α-helices. Cleavage by [(OP)2Cu]+ of CRP-cAMP detected hypersensitivities in β-strands 9 and 11 as well as in the D and E α-helices. [(OP)2Cu]+ also detected protections in the C α-helix , the interdomain hinge, and β-strands 2–7. Binding of cGMP to CRP resulted in changes in susceptibility to cleavage only by [(OP)2Cu]+, which detected minor protections in β-strands 3–7, the interdomain hinge, and the C α-helix. These results show that binding of cAMP causes structural changes in CRP in the nucleotide binding domain, the interdomain hinge, the DNA binding domain, and regions involved in intersubunit interaction. Structural changes induced by binding of cGMP appear to be very minor and confined to the nucleotide binding domain, the interdomain hinge, and regions involved in intersubunit interaction. Use of different cleaving agents in protein footprinting seems to give a more detailed picture of structural changes than the use of a single protease alone.

δ-Crystallin, the major soluble protein component of the avian and reptilian eye lens, is homologous to the urea cycle enzyme argininosuccinate lyase (ASL). In duck lenses there are two δ crystallins, denoted δ1 and δ2. Duck δ2 is both a major structural protein of the lens and also the duck orthologue of ASL, an example of gene recruitment. Although 94% identical to δ2/ASL in the amino acid sequence, δ1 is enzymatically inactive. A series of hybrid proteins have been constructed to assess the role of each structural domain in the enzymatic mechanism. Five chimeras—221, 122, 121, 211, and 112, where the three numbers correspond to the three structural domains and the value of 1 or 2 represents the protein of origin, δ1 or δ2, respectively—were constructed and thermodynamically and kinetically analyzed. The kinetic analysis indicates that only domain 1 is crucial for restoring ASL activity to δ1 crystallin, and that amino acid substitutions in domain 2 may play a role in substrate binding. These results confirm the hypothesis that only one domain, domain 1, is responsible for the loss of catalytic activity in δ1. The thermodynamic characterization of human ASL (hASL) and duck δ1 and δ2 indicate that δ crystallins are slightly less stable than hASL, with the δ1 being the least stable. The ΔGs of unfolding are 57.25, 63.13, and 70.71 kcal mol−1 for δ1, δ2, and hASL, respectively. This result was unexpected, and we speculate that δ crystallins have adapted to their structural role by adopting a slightly less stable conformation that might allow for enhanced protein–protein and protein–solvent interactions.

Thymidylate synthase (TS) converts dUMP to dTMP by reductive methylation, where 5,10-methylenetetrahydrofolate is the source of both the methylene group and reducing equivalents. The mechanism of this reaction has been extensively studied, mainly using the enzyme from Escherichia coli. Bacillus subtilis contains two genes for TSs, ThyA and ThyB. The ThyB enzyme is very similar to other bacterial TSs, but the ThyA enzyme is quite different, both in sequence and activity. In ThyA TS, the active site histidine is replaced by valine. In addition, the B. subtilis enzyme has a 2.4-fold greater kcat than the E. coli enzyme. The structure of B. subtilis thymidylate synthase in a ternary complex with 5-fluoro-dUMP and 5,10-methylenetetrahydrofolate has been determined to 2.5 Å resolution. Overall, the structure of B. subtilis TS (ThyA) is similar to that of the E. coli enzyme. However, there are significant differences in the structures of two loops, the dimer interface and the details of the active site. The effects of the replacement of histidine by valine and a serine to glutamine substitution in the active site area, and the addition of a loop over the carboxy terminus may account for the differences in kcat found between the two enzymes.

The heat stable inhibitor of cAMP-dependent protein kinase (PKIα) contains both a nuclear export signal (NES) and a high affinity inhibitory region that is essential for inhibition of the catalytic subunit of the kinase. These functions are sequentially independent. Two-dimensional NMR spectroscopy was performed on uniformly [15N]-labeled PKIα to examine its structure free in solution. Seventy out of 75 residues were identified, and examination of the CαH chemical shifts revealed two regions of upfield chemical shifts characteristic of α-helices. When PKIα was fragmented into two functionally distinct peptides for study at higher concentrations, no significant alterations in chemical shifts or secondary structure were observed. The first ordered region, identified in PKIα (1–25), contains an α-helix from residues 1–13. This helix extends by one turn the helix observed in the crystal structure of a PKIα (5–24) peptide bound to the catalytic subunit. The second region of well-defined secondary structure, residues 35–47, overlaps with the nuclear export signal in the PKIα (26–75) fragment. This secondary structure consists of a helix with a hydrophobic face comprised of Leu37, Leu41, and Leu44, followed by a flexible turn containing Ile46. These four residues are critical for nuclear export function. The remainder of the protein in solution appears relatively unstructured, and this lack of structure surrounding a few essential and well-defined signaling elements may be characteristic of a growing family of small regulatory proteins that interact with protein kinases.

The energetics of structural changes in the holo and apo forms of α-lactalbumin and the transition between their native and denatured states induced by binding Ca2+ and Na+ have been studied by differential scanning and isothermal titration microcalorimetry and circular dichroism spectroscopy under various solvent conditions. Removal of Ca2+ from the protein enhances its sensitivity to pH and ionic conditions due to noncompensated negative charge–charge interactions at the cation binding site, which significantly reduces its overall stability. At neutral pH and low ionic strength, the native structure of apo-α-lactalbumin is stable below 14 °C and undergoes a conformational change to a native-like molten globule intermediate at temperatures above 25 °C. The denaturation of either holo- or apo-α-lactalbumin is a highly cooperative process that is characterized by an enthalpy of similar magnitude when calculated at the same temperature. Measured by direct calorimetric titration, the enthalpy of Ca2+-binding to apo-LA at pH 7.5 is −7.1 kJ mol−1 at 5.0 °C, which is essentially invariant to protonation effects. This small enthalpy effect infers that stabilization of α-lactalbumin by Ca2+ is primarily an entropy driven process, presumably arising from electrostatic interactions and the hydration effect. In contrast to the binding of calcium, the interaction of sodium with apo-LA does not produce a noticeable heat effect and is characterized by its ionic nature rather than specific binding to the metal-binding site. Characterization of the conformational stability and ligand binding energetics of α-lactalbumin as a function of solvent conditions furnishes significant insight regarding the molecular flexibility and regulatory mechanism mediated by this protein.

To explore the structure of the pore-forming fragment of colicin E1 in membranes, a series of 23 consecutive single cysteine substitution mutants was prepared in the sequence 402–424. Each mutant was reacted with a sulfhydryl-specific reagent to generate a nitroxide labeled side chain, and the mobility of the side chain and its accessibility to collision with paramagnetic reagents was determined from the electron paramagnetic resonance spectrum. Individual values of these quantities were used to identify tertiary contact sites and the nature of the surrounding solvent, while their periodic dependence on sequence position was used to identify secondary structure. In solution, the data revealed a regular helix of 11 residues in the region 406–416, consistent with helix IV of the crystal structure. Upon binding to negatively charged membranes at pH 4.0, helix IV apparently grows to a length of 19 residues, extending from 402–420. One face of the helix is solvated by the lipid bilayer, and the other by an environment of a polar nature. Surprisingly, a conserved charged pair, D408–R409, is located on the lipid-exposed face. Evidence is presented to suggest a transmembrane orientation of this new helix, although other topographies may exist in equilibrium.

An α-helix and a β-strand are said to be interactively packed if at least one residue in each of the secondary structural elements loses 10% of its solvent accessible contact area on association with the other secondary structural element. An analysis of all such 5,975 nonidentical α/β units in protein structures, defined at ≤2.5 Å resolution, shows that the interaxial distance between the α-helix and the β-strand is linearly correlated with the residue-dependent function, log[(V/nda)/n-int], where V is the volume of amino acid residues in the packing interface, nda is the normalized difference in solvent accessible contact area of the residues in packed and unpacked secondary structural elements, and n-int is the number of residues in the packing interface. The β-sheet unit (βu), defined as a pair of adjacent parallel or antiparallel hydrogen-bonded β-strands, packing with an α-helix shows a better correlation between the interaxial distance and log(V/nda) for the residues in the packing interface. This packing relationship is shown to be useful in the prediction of interaxial distances in α/β units using the interacting residue information of equivalent α/β units of homologous proteins. It is, therefore, of value in comparative modeling of protein structures.

Left-handed polyproline II helices (PPII) are contiguous elements of protein secondary structure in which the φ and ψ angles of constituent residues are restricted to around −75° and 145°, respectively. They are important in structural proteins, in unfolded states and as ligands for signaling proteins. Here, we present a survey of 274 nonhomologous polypeptide chains from proteins of known structure for regions that form these structures. Such regions are rare, but the majority of proteins contain at least one PPII helix. Most PPII helices are shorter than five residues, although the longest found contained 12 amino acids. Proline predominates in PPII, but Gln and positively charged residues are also favored. The basis of Gln's prevalence is its ability to form an i, i + 1 side-chain to main-chain hydrogen bond with the backbone carbonyl oxygen of the proceeding residue; this helps to fix the ψ angle of the Gln and the φ and ψ of the proceeding residue in PPII conformations and explains why Gln is favored at the first position in a PPII helix. PPII helices are highly solvent exposed, which explains why apolar amino acids are disfavored despite preferring this region of φ/ψ space when in isolation. PPII helices have perfect threefold rotational symmetry and within these structures we find significant correlation between the hydrophobicity of residues at i and i + 3; thus, PPII helices in globular proteins can be considered to be amphipathic.

Glucosamine 6-phosphate synthase converts fructose-6P into glucosamine-6P or glucose-6P depending on the presence or absence of glutamine. The isomerase activity is associated with a 40-kDa C-terminal domain, which has already been characterized crystallographically. Now the three-dimensional structures of the complexes with the reaction product glucose-6P and with the transition state analog 2-amino-2-deoxyglucitol-6P have been determined. Glucose-6P binds in a cyclic form whereas 2-amino-2-deoxyglucitol-6P is in an extended conformation. The information on ligand-protein interactions observed in the crystal structures together with the isotope exchange and site-directed mutagenesis data allow us to propose a mechanism of the isomerase activity of glucosamine-6P synthase. The sugar phosphate isomerization involves a ring opening step catalyzed by His504 and an enolization step with Glu488 catalyzing the hydrogen transfer from C1 to C2 of the substrate. The enediol intermediate is stabilized by a helix dipole and the ε-amino group of Lys603. Lys485 may play a role in deprotonating the hydroxyl O1 of the intermediate.

The protein–protein interaction energy of 12 nonhomologous serine protease-inhibitor and 15 antibody-antigen complexes is calculated using a molecular mechanics formalism and dissected in terms of the main-chain vs. side-chain contribution, nonrotameric side-chain contributions, and amino acid residue type involvement in the interface interaction.

There are major differences in the interactions of the two types of protein–protein complex. Protease-inhibitor complexes interact predominantly through a main-chain–main-chain mechanism while antibody-antigen complexes interact predominantly through a side-chain–side-chain or a side-chain–main-chain mechanism. However, there is no simple correlation between the main-chain–main-chain interaction energy and the percentage of main-chain surface area buried on binding. The interaction energy is equally effected by the presence of nonrotameric side-chain conformations, which constitute ∼20% of the interaction energy. The ability to reproduce the interface interaction energy of the crystal structure if original side-chain conformations are removed from the calculation is much greater in the protease-inhibitor complexes than the antibody-antigen complexes. The success of a rotameric model for protein–protein docking appears dependent on the extent of the main-chain–main-chain contribution to binding.

Analysis of (1) residue type and (2) residue pair interactions at the interface show that antibody-antigen interactions are very restricted with over 70% of the antibody energy attributable to just six residue types (Tyr > Asp > Asn > Ser > Glu > Trp) in agreement with previous studies on residue propensity. However, it is found here that 50% of the antigen energy is attributable to just four residue types (Arg = Lys > Asn > Asp). On average just 12 residue pair interactions (6%) contribute over 40% of the favorable interaction energy in the antibody-antigen complexes, with charge-charge and charge/polar-tyrosine interactions being prominent. In contrast protease inhibitors use a diverse set of residue types and residue pair interactions.

Fold assignments for proteins from the Escherichia coli genome are carried out using BASIC, a profile–profile alignment algorithm, recently tested on fold recognition benchmarks and on the Mycoplasma genitalium genome and PSI BLAST, the newest generation of the de facto standard in homology search algorithms. The fold assignments are followed by automated modeling and the resulting three-dimensional models are analyzed for possible function prediction.

Close to 30% of the proteins encoded in the E. coli genome can be recognized as homologous to a protein family with known structure. Most of these homologies (23% of the entire genome) can be recognized both by PSI BLAST and BASIC algorithms, but the latter recognizes an additional 260 homologies. Previous estimates suggested that only 10–15% of E. coli proteins can be characterized this way. This dramatic increase in the number of recognized homologies between E. coli proteins and structurally characterized protein families is partly due to the rapid increase of the database of known protein structures, but mostly it is due to the significant improvement in prediction algorithms.

Knowing protein structure adds a new dimension to our understanding of its function and the predictions presented here can be used to predict function for uncharacterized proteins. Several examples, analyzed in more detail in this paper, include the DPS protein protecting DNA from oxidative damage (predicted to be homologous to ferritin with iron ion acting as a reducing agent) and the ahpC/tsa family of proteins, which provides resistance to various oxidating agents (predicted to be homologous to glutathione peroxidase).

The trifluoroethanol (TFE)-induced structural changes of two proteins widely used in folding experiments, bovine α-lactalbumin, and bovine pancreatic ribonuclease A, have been investigated. The experiments were performed using circular dichroism spectroscopy in the far- and near-UV region to monitor changes in the secondary and tertiary structures, respectively, and dynamic light scattering to measure the hydrodynamic dimensions and the intermolecular interactions of the proteins in different conformational states. Both proteins behave rather differently under the influence of TFE: α-lactalbumin exhibits a molten globule state at low TFE concentrations before it reaches the so-called TFE state, whereas ribonuclease A is directly transformed into the TFE state at TFE concentrations above 40% (v/v). The properties of the TFE-induced states are compared with those of equilibrium and kinetic intermediate states known from previous work to rationalize the use of TFE in yielding information about the folding of proteins. Additionally, we report on the properties of TFE/water and TFE/buffer mixtures derived from dynamic light scattering investigations under conditions used in our experiments.

Human pancreatic α-amylase (HPA) was expressed in the methylotrophic yeast Pichia pastoris and two mutants (D197A and D197N) of a completely conserved active site carboxylic acid were generated. All recombinant proteins were shown by electrospray ionization mass spectrometry (ESI-MS) to be glycosylated and the site of attachment was shown to be Asn461 by peptide mapping in conjunction with ESI-MS. Treatment of these proteins with endoglycosidase F demonstrated that they contained a single N-linked oligosaccharide and yielded a protein product with a single N-acetyl glucosamine (GlcNAc), which could be crystallized. Solution of the crystal structure to a resolution of 2.0 Å confirmed the location of the glycosyl group as Asn461 and showed that the recombinant protein had essentially the same conformation as the native enzyme. The kinetic parameters of the glycosylated and deglycosylated wild-type proteins were the same while the kcat/Km values for D197A and D197N were 106–107 times lower than the wild-type enzyme. The decreased kcat/Km values for the mutants confirm that D197 plays a crucial role in the hydrolytic activity of HPA, presumably as the catalytic nucleophile.

The DNA repair protein RecA of Mycobacterium tuberculosis contains an intein, a self-splicing protein element. We have employed this Mtu recA intein to create a selection system for successful intein splicing by inserting it into a kanamycin-resistance gene so that functional antibiotic resistance can only be restored upon protein splicing. We then proceeded to develop an ORFTRAP, i.e., a selection system for the cloning of open reading frames (ORFs). The ORFTRAP exploits the self-splicing properties of inteins (which depend on full-length in-frame translation of a precursor protein) by allowing protein splicing to occur when DNA fragments encoding ORFs are inserted into the Mtu recA intein, whereas DNA fragments containing non-ORFs are selected against. Regions of the Mtu recA intein that tolerate the insertion of additional amino acids were identified by Bgl II linker scanning mutagenesis, and a respective construct was chosen as the ORFTRAP. To test the maximum insert size that could be cloned into ORFTRAP, DNA fragments of increasing length from the Listeria monocytogenes hly gene as well as a genomic library of Haemophilus influenzae were inserted and it was found that the longest permissive inserts were 425 bp and 251 bp, respectively. The H. influenzae ORFTRAP library also demonstrated the strength (strong selection power) and weakness (insertion of very small fragments) of the system. Further modifications should make the ORFTRAP useful for protein expression, epitope mapping, and antigen screening.

A protein structure comparison method is described that allows the generation of large populations of high-scoring alternate alignments. This was achieved by incorporating a random element into an iterative double dynamic programming algorithm. The maximum scores from repeated comparisons of a pair of structures converged on a value that was taken as the global maximum. This lay 15% over the score obtained from the single fixed (unrandomized) calculation. The effect of the gap penalty was observed through the shift of the alignment populations, characterized by their alignment length and root-mean-square deviation (RMSD). The best (lowest RMSD) values found in these populations provided a base-line against which other methods were compared.