Simultaneous single-structure and bundle representation of protein NMR structures in torsion angle space

A method is introduced to represent an ensemble of conformers of a protein by a single structure in torsion angle space that lies closest to the averaged Cartesian coordinates while maintaining perfect covalent geometry and on average equal steric quality and an equally good fit to the experimental (e.g. NMR) data as the individual conformers of the ensemble. The single representative ‘regmean structure’ is obtained by simulated annealing in torsion angle space with the program CYANA using as input data the experimental restraints, restraints for the atom positions relative to the average Cartesian coordinates, and restraints for the torsion angles relative to the corresponding principal cluster average values of the ensemble. The method was applied to 11 proteins for which NMR structure ensembles are available, and compared to alternative, commonly used simple approaches for selecting a single representative structure, e.g. the structure from the ensemble that best fulfills the experimental and steric restraints, or the structure from the ensemble that has the lowest RMSD value to the average Cartesian coordinates. In all cases our method found a structure in torsion angle space that is significantly closer to the mean coordinates than the alternatives while maintaining the same quality as individual conformers. The method is thus suitable to generate representative single structure representations of protein structure ensembles in torsion angle space. Since in the case of NMR structure calculations with CYANA the single structure is calculated in the same way as the individual conformers except that weak positional and torsion angle restraints are added, we propose to represent new NMR structures by a ‘regmean bundle’ consisting of the single representative structure as the first conformer and all but one original individual conformers (the original conformer with the highest target function value is discarded in order to keep the number of conformers in the bundle constant). In this way, analyses that require a single structure can be carried out in the most meaningful way using the first model, while at the same time the additional information contained in the ensemble remains available.     

Conformational study of palindromic tripeptides (GPG, IPI and KPK) in HIV-1 protease--a density functional theory study

A comparative study has been carried out on three palindromic tripeptides Gly-Pro-Gly, Ile-Pro-Ile and Lys-Pro-Lys which were present in HIV protein along with their analogues applying density functional computation at B3LYP/6-31G* level of theory. Discrepancy from the structural analysis has been noted for all the systems and it was found to be more for amide capped structure at the C terminal of proline. The puckering amplitude A and Phase angle P of the pyrrolidine ring of proline in the chosen palindromic tripeptides and their analogues were calculated from the endocyclic torsion angles. The minimum energy conformers lying well within the prescribed region of proline were obtained for the derived compounds from potential energy surface scan mentioning that no role has been played by its terminal residues. This is further supported by the simulated amide bands identifying the helical structure for all three palindromic tripeptides signifying the importance of proline. The molecular properties such as stabilization energy, chemical hardness along with dipole moment were calculated and interpreted. The values of Calpha-H(s) and the peptide backbone N-Calpha-CO for all the selected conformers specify the three palindromic tripeptides to have a symmetrical achiral structure.     

Ab initio quantum mechanical models of peptide helices and their vibrational spectra

Structural parameters for standard peptide helices (alpha, 3(10), 3(1) left-handed) were fully ab initio optimized for Ac-(L-Ala)(9)-NHMe and for Ac-(L-Pro)(9)-NHMe (poly-L-proline-PLP I and PLP II-forms), in order to better understand the relative stability and minimum energy geometries of these conformers and the dependence of the ir absorption and vibrational CD (VCD) spectra on detailed variation in these conformations. Only the 3(10)-helical Ala-based conformation was stable in vacuum for this decaamide structure, but both Pro-based conformers minimized successfully. Inclusion of solvent effects, by use of the conductor-like screening solvent model (COSMO), enabled ab initio optimizations [at the DFT/B3LYP/SV(P) level] without any constraints for the alpha- and 3(10)-helical Ala-based peptides as well as the two Pro-based peptides. The geometries obtained compare well with peptide chain torsion angles and hydrogen-bond distances found for these secondary structure types in x-ray structures of peptides and proteins. For the simulation of VCD spectra, force field and intensity response tensors were obtained ab initio for the complete Ala-based peptides in vacuum, but constrained to the COSMO optimized torsional angles, due to limitations of the solvent model. Resultant spectral patterns reproduce well many aspects of the experimental spectra and capture the differences observed for these various helical types.     

Functional domains of Drosophila UNR in translational control

Translational repression of male-specific-lethal 2 (msl-2) mRNA by Sex-lethal (SXL) is an essential regulatory step of X chromosome dosage compensation in Drosophila. Translation inhibition requires that SXL recruits the protein upstream of N-ras (UNR) to the 3' UTR of msl-2 mRNA. UNR is a conserved, ubiquitous protein that contains five cold-shock domains (CSDs). Here, we dissect the domains of UNR required for translational repression and complex formation with SXL and msl-2 mRNA. Using gel-mobility shift assays, the domain involved in interactions with SXL and msl-2 was mapped specifically to the first CSD (CSD1). Indeed, excess of a peptide containing this domain derepressed msl-2 translation in vitro. The CSD1 of human UNR can also form a complex with SXL and msl-2. Comparative analyses of the CSDs of the Drosophila and human proteins together with site-directed mutagenesis experiments revealed that three exposed residues within CSD1 are required for complex formation. Tethering assays showed that CSD1 is not sufficient for translational repression, indicating that UNR binding to SXL and msl-2 can be distinguished from translation inhibition. Repression by tethered UNR requires residues from both the amino-terminal Q-rich stretch and the two first CSDs, indicating that the translational effector domain of UNR resides within the first 397 amino acids of the protein. Our results identify domains and residues required for UNR function in translational control.     

Structural insights for designed alanine-rich helices: comparing NMR helicity measures and conformational ensembles from molecular dynamics simulation

The temperature dependence of helical propensities for the peptides Ac-ZGG-(KAAAA)(3)X-NH(2) (Z = Y or G, X = A, K, and D-Arg) were studied both experimentally and by MD simulations. Good agreement is observed in both the absolute helical propensities as well as relative helical content along the sequence; the global minimum on the calculated free energy landscape corresponds to a single alpha-helical conformation running from K4 to A18 with some terminal fraying, particularly at the C-terminus. Energy component analysis shows that the single helix state has favorable intramolecular electrostatic energy due to hydrogen bonds, and that less-favorable two-helix globular states have favorable solvation energy. The central lysine residues do not appear to increase helicity; however, both experimental and simulation studies show increasing helicity in the series X = Ala --> Lys --> D-Arg. This C-capping preference was also experimentally confirmed in Ac-(KAAAA)(3)X-GY-NH(2) and (KAAAA)(3)X-GY-NH(2) sequences. The roles of the C-capping groups, and of lysines throughout the sequence, in the MD-derived ensembles are analyzed in detail.     

Identification of neutral and anionic 8 alpha-substituted flavin semiquinones in flavoproteins by electron spin resonance spectroscopy

Circular dichroic spectra of metmyoglobin and apomyoglobin were measured in neutral and acidic solution. Addition of sodium dodecyl sulfate (NaDodSO₄) slightly reduces the helicity (based on the circular dichroic magnitude) of both proteins probably because of the loss of long-range interactions among helical segments. Lowering the pH of the protein-surfactant solution to 3 slightly enhances the helical conformation of myoglobin due to the protonation of acidic side groups and thereby the reduction of coulombic repulsion among negative charges. For BrCN-digested fragments the COOH-terminal peptide (22 residues) loses its helicity which can be restored by addition of NaDodSO₄. The middle fragment (76 residues) retains a considerable amount of helicity in water alone, which further increases in the presence of NaDodSO₄. The NH₂-terminal fragment (55 residues) also has some helical conformation in water, which is enhanced by the addition of NaDodSO₄. The circular dichroic spectrum of an equimolar mixture of the three peptides in NaDodSO₄ solution is the same as that calculated from the spectra of isolated peptides under similar conditions.     

Conformational solution studies of neuropeptide gamma using CD and NMR spectroscopy.

Neuropeptide gamma is one of the largest members of the tachykinin family of peptides, exhibiting strong agonistic activity towards the NK-2 tachykinin receptor. This peptide was synthesized by the solid-phase method using the Fmoc chemistry. Circular-dichroism spectroscopy (CD) investigations of this peptide were performed in phosphate buffer, in the presence of sodium dodecylsulphate (SDS) micelles and trifluoroethanol (TFE) solutions and in DMSO-d6 using the 2D NMR technique in conjunction with two different theoretical approaches. The first assumes multiconformational equilibrium of the peptide studied characterized by the values of statistical weights of low-energy conformations. These calculations were performed using three different force fields ECEPP/3, AMBER4.1 and CHARMM (implemented in the X-PLOR program). The second method incorporates interproton distance and dihedral angle constraints into the starting conformation using the Simulated Annealing algorithm (X-PLOR program). The CD experiments revealed that although the peptide studied is flexible in polar solvents, a tendency to adopt a helical structure was observed in the hydrophobic environment. The NMR data (NOE effects) indicate a helical or reverse structure in the Ile7-His12 fragment of the peptide studied in DMSO-d6 solution. The results obtained cannot be interpreted in terms of a single conformation. Most of the conformations obtained with the ECEPP/3 force field possess a high content of a helical structure. None of the conformers, obtained with the AMBER4.1 and CHARMM force fields, can be considered as the dominant one. In all conformations several beta-turns were detected and in some cases gamma-turns were also found. But in fact, it is rather difficult to select the position of the secondary element(s) present in the structure of NPgamma in solution. All conformers calculated with the X-PLOR program (with using NMR derived distance and torsion angle constraints) are stabilized by several beta-turns. Common structural motives are a type IV beta-turn in the Gln6-His12 fragment. All conformations obtained using two approaches adopt very similar turn shapes in the middle region of molecule and a random structure on the N- and C-terminal fragments.     

Structure–function relationships in histidine-rich antimicrobial peptides from Atlantic cod

Gad-1 and Gad-2 are antimicrobial peptide (AMP) sequences encoded by paralogous genes. They are rich in histidine, which suggests that their activity might be pH-dependent. We examined their structure–function relationships with a view to learning how to improve AMP therapeutic ratios. Activity assays with Gram-negative bacteria and cancer cell lines demonstrate that Gad-2 is substantially more active at slightly acidic pH than it is at neutral pH. By contrast, the activity of Gad-1 at lower pH is similar to its activity at pH 7. Circular dichroism spectra indicate that the greater functional plasticity of Gad-2 correlates with a greater structural plasticity; Gad-2's percent helicity varies dramatically with altered pH and lipid environment. Interestingly, Gad-2's highest levels of helicity do not correspond to the conditions where it is most active. High resolution solution NMR structures were determined in SDS micelles at pH 5, conditions that induce an intermediate level of helicity in the peptides. Gad-1 is more helical than Gad-2, with both peptides exhibiting the greatest helical tendencies in their central region and lowest helicity in their N-termini. The high resolution structures suggest that maximum activity relies on the appropriate balance between an N-terminal region with mixed hydrophobic/hydrophilic structure features and an amphipathic central and C-terminal region. Taken together with previous studies, our results suggest that to improve the therapeutic ratio of AMPs, consideration should be given to including sequential histidine-pairs, keeping the overall charge of the peptide modest, and retaining a degree of structural plasticity and imperfect amphipathicity.     

Solution structure of human growth hormone-releasing factor fragment (1-29) by CD: characteristic conformational change on phospholipid membrane

The conformations of synthetic human growth hormone-releasing factor fragment (1-29) in the presence and the absence of 1,2-dimyristoyl-sn-glycero-3-phosphorylglycerol liposome as well as in aqueous 2,2,2-trifluoroethanol solution were investigated by CD spectroscopy. The secondary structure of the peptide in each solution was analyzed by two methods. Both results show that the peptide has an unordered structure in the aqueous solution, whereas it folds into helical structure in the aqueous alcohol and in the phospholipid solution. In addition, although the peptide exists as almost complete helix in the 50 vol% aqueous alcohol (80-90% helicity), it does not reach full helicity even in the solution containing excess amount of phospholipid liposome (maximum 65-70% helicity). The conformational difference is explained by the characteristic amphipathy of the peptide, i.e., the necessity to twist the separated amphipathic helical parts in the interaction with the phospholipid membrane probably makes the helicity of the peptide decrease.     

The unusual conformation of cross-strand disulfide bonds is critical to the stability of β-hairpin peptides

The cross-strand disulfides (CSDs) found in β-hairpin antimicrobial peptides (β-AMPs) show a unique disulfide geometry that is characterized by unusual torsion angles and a short Cα-Cα distance. While the sequence and disulfide bond connectivity of disulfide-rich peptides is well studied, much less is known about the disulfide geometry found in CSDs and their role in the stability of β-AMPs. To address this, we solved the nuclear magnetic resonance (NMR) structure of the β-AMP gomesin (Gm) at 278, 298, and 310 K, examined the disulfide bond geometry of over 800 disulfide-rich peptides, and carried out extensive molecular dynamics (MD) simulation of the peptides Gm and protegrin. The NMR data suggests Cα-Cα distances characteristic for CSDs are independent of temperature. Analysis of disulfide-rich peptides from the Protein Data Bank revealed that right-handed and left-handed rotamers are equally likely in CSDs. The previously reported preference for right-handed rotamers was likely biased by restricting the analysis to peptides and proteins solved using X-ray crystallography. Furthermore, data from MD simulations showed that the short Cα-Cα distance is critical for the stability of these peptides. The unique disulfide geometry of CSDs poses a challenge to biomolecular force fields and to retain the stability of β-hairpin fold over long simulation times, restraints on the torsion angles might be required.     

Conformational analysis of bacterial cell wall peptides indicates how particular conformations have influenced the evolution of penicillin-binding proteins,beta-lactam antibiotics and antibiotic resistance mechanisms

Our aim was to use a conformational analysis technique developed for peptides to identify structural relationships between bacterial cell wall peptides and beta-lactam antibiotics that might help to explain their different actions as substrates and inhibitors of penicillin binding proteins (PBPs). The conformational forms of the model cell wall peptide Ac-L-Lys(Ac)-D-Ala-D-Ala are described by just a few backbone torsion combinations: three C-terminal carboxylate regions, with Tor8 (psi(i+1)) ranges of D3 region (50 degrees to 70 degrees ), D6 region (140 degrees to 170 degrees ) and D9 region (-50 degrees to -70 degrees ) are combined with either of two Tor6 (phi(i))-Tor4 (psi(i)) combinations, C4 region (-50 degrees to -80 degrees ) with B8 region (-40 degrees to -70 degrees ) or C11 region (30 degrees to 50 degrees ) with B2 region (30 degrees to 70 degrees ). From these results, and comparisons with conformational analyses of various beta-lactams and Ac-L-Lys(Ac)-D-Ala-D-Lac, it is concluded that molecular recognition of cell wall peptide substrates by PBPs requires conformers with backbone torsion angles of D3C4B8. beta-Lactam antibiotics are constrained compounds with fewer conformational forms; these match well the backbone torsions of cell wall peptides at D3C4, allowing their recognition and acylation by PBPs, whereas their unique Tor4 produces differently orientated CO and N atoms that appear to prevent subsequent deacylation, leading to their action as suicide substrates. The results are also related to the selective pressures involved in evolution of beta-lactamases from PBPs. From analysis of conformers of Ac-L-Lys(Ac)-D-Ala-D-Ala and the vancomycin-resistant analogue Ac-L-Lys(Ac)-D-Ala-D-Lac, it is concluded that vancomycin may recognise D6C11B2 conformers, giving it complementary substrate specificity to PBPs. This approach could have applications in the rational design of antibiotics targeted against PBPs and their substrates.     

Conformational analysis of endothelin-1: Effects of solvation free energy

In order to investigate conformational preferences of the 21-residue peptide hormone endothelin-1 (ET-1), an extensive conformational search was carried out in vacuo using a combination of high temperature molecular dynamics / annealing and a Monte Carlo / minimization search in torsion angle space. Fully minimized conformations from the search were grouped into families using a clustering technique based on rms fitting over the Cartesian coordinates of the atoms of the peptide backbone of the ring region. A wide range of local energy minima were identified even though two disulfide bridges (Cys¹-Cys¹⁵ and Cys³-Cys¹¹) constrain the structure of the peptide. Low energy conformers of ET-1 as a nonionized species in vacuo arestabilized by intramolecular interaction of the ring region (residues 1-15) with the tail (residues 16–21). Strained conformations for individual residues are observed. Conformational similarity to protein loops is established by matching to protein crystal structures In order to assess the influence of aqueous environment on conformational preference, the electrostatic contribution to the solvation energy was calculated for ET-1 as a fully ionized species (Asp⁸, Lys⁹, Glu¹⁰, Asp¹⁸, N- and C-terminus) using a continuum electrostatics model (DelPhi) for each of the conformed generated in vacuo, and the total solvation free energy was estimated by adding a hydrophobic contribution proportional to solvent accessible surface area. Solvation dramatically alters the relative energetics of ET-1 conformers from that calculated in vacuo. Conformers of ET-1 favored by the electrostatic salvation energy in water include conformers with helical secondary structure in the region of residues 9–15. Perhaps of most importance, it was demonstrated that the contribution tosolvation by an individual charge depends not only on its solvent accessibility but on the proximity of other charges, i.e., it is a cooperative effect. This was shown by the calculation of electrostatic solvation energy as afunction of conformation with individual charges systematically turned “on” and “off”. The cooperative effect of multiple charges on solvation demonstrated in this manner calls into question models that relate solvation energysimply to solvent accessibility by atom or residue alone. © 1995 John Wiley & Sons, Inc.     

Folding of peptide fragments comprising the complete sequence of proteins. Models for initiation of protein folding. I. Myohemerythrin.

In an attempt to delineate potential folding initiation sites for different protein structural motifs, we have synthesized series of peptides that span the entire length of the polypeptide chain of two proteins, and examined their conformational preferences in aqueous solution using proton nuclear magnetic resonance and circular dichroism spectroscopy. We describe here the behavior of peptides derived from a simple four-helix bundle protein, myohemerythrin. The peptides correspond to the sequences of the four long helices (the A, B, C and D helices), the N- and C-terminal loops and the connecting sequences between the helices. The peptides corresponding to the helices of the folded protein all exhibit preferences for helix-like conformations in solution. The conformational ensembles of the A- and D-helix peptides contain ordered helical forms, as shown by extensive series of medium-range nuclear Overhauser effect connectivities, while the B- and C-helix peptides exhibit conformational preferences for nascent helix. All four peptides adopt ordered helical conformations in mixtures of trifluoroethanol and water. The terminal and interconnecting loop peptides also appear to contain appreciable populations of conformers with backbone phi and psi angles in the alpha-region and include highly populated hydrophobic cluster and/or turn conformations in some cases. Trifluoroethanol is unable to drive these peptides towards helical conformations. Overall, the peptide fragments of myohemerythrin have a marked preference towards secondary structure formation in aqueous solution. In contrast, peptide fragments derived from the beta-sandwich protein plastocyanin are relatively devoid of secondary structure in aqueous solution (see accompanying paper). These results suggest that the two different protein structural motifs may require different propensities for formation of local elements of secondary structure to initiate folding, and that there is a prepartitioning of conformational space determined by the local amino acid sequence that is different for the helical and beta-sandwich structural motifs.     

Conformational limitations of glycylsarcosine as a prototypic substrate for peptide transporters

Peptide transporters are present in all species to absorb the small peptides that occur ubiquitously as products of proteolysis. The broad substrate specificities of these systems allow them to be exploited therapeutically for delivery of peptidomimetic drugs in microbes and man. To this end, glycylsarcosine is currently used as a standard substrate for assaying peptidomimetic transport by peptide transporters. However, in this study we find it is unsuitable as a general substrate, based on assays of its transport by model bacterial peptide transporters and computer-based conformational analysis of its structure. Of the two generic transporters for di- and tripeptides, exemplified by Dpp and Tpp in Escherichia coli, only Dpp can transport glycylsarcosine. The explanation for this finding came from molecular modelling, which indicated that glycylsarcosine can adopt only a restricted range of conformers compared with typical dipeptides, and that of the conformers with a trans peptide bond, the majority have the specific psi and phi backbone torsion angles needed for molecular recognition and transport by Dpp but none possessed psi and phi torsions required for recognition by Tpp; moreover, 38% of its conformers have cis peptide bonds that are not substrates for any peptide transporter. Thus, using glycylsarcosine as substrate in competition assays with compounds that typically form conformers recognised by both types of peptide transporter will underestimate their transport. These findings have implications for assays of oral availability of peptidomimetic drugs such as beta-lactams, ACE inhibitors and anti-viral compounds, for which glycylsarcosine is routinely used.     

Distributed torsion angle grid search in high dimensions: A systematic approach to NMR structure determination

Two complementary approaches for systematic search in torsion angle space are described for the generation of all conformations of polypeptides which satisfy experimental NMR restraints, hard-sphere van der Waals radii, and rigid covalent geometry. The first procedure is based on a recursive, tree search algorithm for the examination of linear chains of torsion angles, and uses a novel treatment to propagate the search results to neighboring regions so that the structural consequences of the restraints are fully realized. The second procedure is based on a binary combination of torsion vector spaces for connected submolecules, and produces intermediate results in Cartesian space for a more robust restraint analysis. Restraints for NMR applications include bounds on torsion angles and internuclear distances, including relational and degenerate restraints involving equivalent and nonstereoassigned protons. To illustrate these methods, conformation search results are given for the tetrapeptide APGA restrained to an idealized -turn conformation, an alanine octapeptide restrained to a right-handed helical conformation, and the structured region of the peptide SYPFDV.     

The helical structure of surfactant peptide KL4 when bound to POPC: POPG lipid vesicles

KL 4 is a 21-residue peptide employed as a functional mimic of lung surfactant protein B, which successfully lowers surface tension in the alveoli. A mechanistic understanding of how KL 4 affects lipid properties has proven elusive as the secondary structure of KL 4 in lipid preparations has not been determined at high resolution. The sequence of KL 4 is based on the C-terminus of SP-B, a naturally occurring helical protein that binds to lipid interfaces. The spacing of the lysine residues in KL 4 precludes the formation of a canonical amphipathic alpha-helix; qualitative measurements using Raman, CD, and FTIR spectroscopies have given conflicting results as to the secondary structure of the peptide as well as its orientation in the lipid environment. Here, we present a structural model of KL 4 bound to lipid bilayers based on solid state NMR data. Double-quantum correlation experiments employing (13)C-enriched peptides were used to quantitatively determine the backbone torsion angles in KL 4 at several positions. These measurements, coupled with CD experiments, verify the helical nature of KL 4 when bound to lipids, with (phi, psi) angles that differ substantially from common values for alpha-helices of (-60, -45). The average torsion angles found for KL 4 bound to POPC:POPG lipid vesicles are (-105, -30); this deviation from ideal alpha-helical structure allows KL 4 to form an amphipathic helix at the lipid interface.     

PONDEROSA, an automated 3D-NOESY peak picking program, enables automated protein structure determination

SUMMARY: PONDEROSA (Peak-picking Of Noe Data Enabled by Restriction of Shift Assignments) accepts input information consisting of a protein sequence, backbone and sidechain NMR resonance assignments, and 3D-NOESY ((13)C-edited and/or (15)N-edited) spectra, and returns assignments of NOESY crosspeaks, distance and angle constraints, and a reliable NMR structure represented by a family of conformers. PONDEROSA incorporates and integrates external software packages (TALOS+, STRIDE and CYANA) to carry out different steps in the structure determination. PONDEROSA implements internal functions that identify and validate NOESY peak assignments and assess the quality of the calculated three-dimensional structure of the protein. The robustness of the analysis results from PONDEROSA's hierarchical processing steps that involve iterative interaction among the internal and external modules. PONDEROSA supports a variety of input formats: SPARKY assignment table (.shifts) and spectrum file formats (.ucsf), XEASY proton file format (.prot), and NMR-STAR format (.star). To demonstrate the utility of PONDEROSA, we used the package to determine 3D structures of two proteins: human ubiquitin and Escherichia coli iron-sulfur scaffold protein variant IscU(D39A). The automatically generated structural constraints and ensembles of conformers were as good as or better than those determined previously by much less automated means. AVAILABILITY: The program, in the form of binary code along with tutorials and reference manuals, is available at http://ponderosa.nmrfam.wisc.edu/.     

DNAase I hypersensitive sites may be correlated with genomic regions of large structural variation

Helical-twist, roll and torsion-angle variations calculated by the Calladine (1982)-Dickerson (1983) rules were scanned along several nucleotide sequences for which DNAase I cleavage data are available. It has been shown that for short synthetic oligomers DNAase I cuts preferentially at positions of high helical twist (Dickerson & Drew, 1981; Lomonossoff et al., 1981). Our calculations indicate that DNAase I sensitive and hypersensitive sites in chromatin are correlated with regions of successive, large, helical-twist angle variations from regular B-DNA. In many cases these regions exhibit large variations in base-pair roll and backbone torsion angles as well. It has been suggested that DNAase I cuts in the vicinity of cruciforms. However, it was recently demonstrated by Courey & Wang (1983) and Gellert et al. (1983) that such cruciform formation in a negatively supercoiled DNA is kinetically forbidden under physiological conditions. We thus propose that clustering of large twist-angle (and/or roll and backbone torsion angle) variations may be among the conformational features recognized by the enzyme. Specific cuts can then preferentially occur at base-pair steps with high helical twists.