Minimal ensembles of side chain conformers for modeling protein–protein interactions

The goal of this article is to reduce the complexity of the side chain search within docking problems. We apply six methods of generating side chain conformers to unbound protein structures and determine their ability of obtaining the bound conformation in small ensembles of conformers. Methods are evaluated in terms of the positions of side chain end groups. Results for 68 protein complexes yield two important observations. First, the end‐group positions change less than 1 Å on association for over 60% of interface side chains. Thus, the unbound protein structure carries substantial information about the side chains in the bound state, and the inclusion of the unbound conformation into the ensemble of conformers is very beneficial. Second, considering each surface side chain separately in its protein environment, small ensembles of low‐energy states include the bound conformation for a large fraction of side chains. In particular, the ensemble consisting of the unbound conformation and the two highest probability predicted conformers includes the bound conformer with an accuracy of 1 Å for 78% of interface side chains. As more than 60% of the interface side chains have only one conformer and many others only a few, these ensembles of low‐energy states substantially reduce the complexity of side chain search in docking problems. This approach was already used for finding pockets in protein–protein interfaces that can bind small molecules to potentially disrupt protein–protein interactions. Side‐chain search with the reduced search space will also be incorporated into protein docking algorithms. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

IRAA: A statistical tool for investigating a protein–protein interaction interface from multiple structures

Understanding protein–protein interactions (PPIs) is fundamental to infer how different molecular systems work. A major component to model molecular recognition is the buried surface area (BSA), that is, the area that becomes inaccessible to solvent upon complex formation. To date, many attempts tried to connect BSA to molecular recognition principles, and in particular, to the underlying binding affinity. However, the most popular approach to calculate BSA is to use a single (or in some cases few) bound structures, consequently neglecting a wealth of structural information of the interacting proteins derived from ensembles corresponding to their unbound and bound states. Moreover, the most popular method inherently assumes the component proteins to bind as rigid entities. To address the above shortcomings, we developed a Monte Carlo method-based Interface Residue Assessment Algorithm (IRAA), to calculate a combined distribution of BSA for a given complex. Further, we apply our algorithm to human ACE2 and SARS-CoV-2 Spike protein complex, a system of prime importance. Results show a much broader distribution of BSA compared to that obtained from only the bound structure or structures and extended residue members of the interface with implications to the underlying biomolecular recognition. We derive that specific interface residues of ACE2 and of S-protein are consistently highly flexible, whereas other residues systematically show minor conformational variations. In effect, IRAA facilitates the use of all available structural data for any biomolecular complex of interest, extracting quantitative parameters with statistical significance, thereby providing a deeper biophysical understanding of the molecular system under investigation.     

Shape complementarity of protein-protein complexes at multiple resolutions

Biological complexes typically exhibit intermolecular interfaces of high shape complementarity. Many computational docking approaches use this surface complementarity as a guide in the search for predicting the structures of protein-protein complexes. Proteins often undergo conformational changes to create a highly complementary interface when associating. These conformational changes are a major cause of failure for automated docking procedures when predicting binding modes between proteins using their unbound conformations. Low resolution surfaces in which high frequency geometric details are omitted have been used to address this problem. These smoothed, or blurred, surfaces are expected to minimize the differences between free and bound structures, especially those that are due to side chain conformations or small backbone deviations. Despite the fact that this approach has been used in many docking protocols, there has yet to be a systematic study of the effects of such surface smoothing on the shape complementarity of the resulting interfaces. Here we investigate this question by computing shape complementarity of a set of 66 protein-protein complexes represented by multiresolution blurred surfaces. Complexed and unbound structures are available for these protein-protein complexes. They are a subset of complexes from a nonredundant docking benchmark selected for rigidity (i.e. the proteins undergo limited conformational changes between their bound and unbound states). In this work, we construct the surfaces by isocontouring a density map obtained by accumulating the densities of Gaussian functions placed at all atom centers of the molecule. The smoothness or resolution is specified by a Gaussian fall-off coefficient, termed "blobbyness." Shape complementarity is quantified using a histogram of the shortest distances between two proteins' surface mesh vertices for both the crystallographic complexes and the complexes built using the protein structures in their unbound conformation. The histograms calculated for the bound complex structures demonstrate that medium resolution smoothing (blobbyness = -0.9) can reproduce about 88% of the shape complementarity of atomic resolution surfaces. Complexes formed from the free component structures show a partial loss of shape complementarity (more overlaps and gaps) with the atomic resolution surfaces. For surfaces smoothed to low resolution (blobbyness = -0.3), we find more consistency of shape complementarity between the complexed and free cases. To further reduce bad contacts without significantly impacting the good contacts we introduce another blurred surface, in which the Gaussian densities of flexible atoms are reduced. From these results we discuss the use of shape complementarity in protein-protein docking.     

The relationship between the flexibility of proteins and their conformational states on forming protein-protein complexes with an application to protein-protein docking

We investigate the extent to which the conformational fluctuations of proteins in solution reflect the conformational changes that they undergo when they form binary protein-protein complexes. To do this, we study a set of 41 proteins that form such complexes and whose three-dimensional structures are known, both bound in the complex and unbound. We carry out molecular dynamics simulations of each protein, starting from the unbound structure, and analyze the resulting conformational fluctuations in trajectories of 5 ns in length, comparing with the structure in the complex. It is found that fluctuations take some parts of the molecules into regions of conformational space close to the bound state (or give information about it), but at no point in the simulation does each protein as whole sample the complete bound state. Subsequent use of conformations from a clustered MD ensemble in rigid-body docking is nevertheless partially successful when compared to docking the unbound conformations, as long as the unbound conformations are themselves included with the MD conformations and the whole globally rescored. For one key example where sub-domain motion is present, a ribonuclease inhibitor, principal components analysis of the MD was applied and was also able to produce conformations for docking that gave enhanced results compared to the unbound. The most significant finding is that core interface residues show a tendency to be less mobile (by size of fluctuation or entropy) than the rest of the surface even when the other binding partner is absent, and conversely the peripheral interface residues are more mobile. This surprising result, consistent across up to 40 of the 41 proteins, suggests different roles for these regions in protein recognition and binding, and suggests ways that docking algorithms could be improved by treating these regions differently in the docking process.     

Relative solvent accessible surface area predicts protein conformational changes upon binding

Protein interactions are often accompanied by significant changes in conformation. We have analyzed the relationships between protein structures and the conformational changes they undergo upon binding. Based upon this, we introduce a simple measure, the relative solvent accessible surface area, which can be used to predict the magnitude of binding-induced conformational changes from the structures of either monomeric proteins or bound subunits. Applying this to a large set of protein complexes suggests that large conformational changes upon binding are common. In addition, we observe considerable enrichment of intrinsically disordered sequences in proteins predicted to undergo large conformational changes. Finally, we demonstrate that the relative solvent accessible surface area of monomeric proteins can be used as a simple proxy for protein flexibility. This reveals a powerful connection between the flexibility of unbound proteins and their binding-induced conformational changes, consistent with the conformational selection model of molecular recognition.     

Sampling the conformation of protein surface residues for flexible protein docking

Background  

The problem of determining the physical conformation of a protein dimer, given the structures of the two interacting proteins in their unbound state, is a difficult one. The location of the docking interface is determined largely by geometric complementarity, but finding complementary geometry is complicated by the flexibility of the backbone and side-chains of both proteins. We seek to generate candidates for docking that approximate the bound state well, even in cases where there is backbone and/or side-chain difference from unbound to bound states.     

Docking without docking: ISEARCH--prediction of interactions using known interfaces

The increasing number of solved protein structures provides a solid number of interfaces, if protein-protein interactions, domain-domain contacts, and contacts between biological units are taken into account. An interface library gives us the opportunity to identify surface regions on a target molecule that are similar by local structure and residue composition. If both unbound components of a possible protein complex exhibit structural similarities to a known interface, the unbound structures can be superposed onto the known interfaces. The approach is accompanied by two mathematical problems. Protein surfaces have to be quickly screened by thousands of patches, and similarity has to be evaluated by a suitable scoring scheme. The used algorithm (NeedleHaystack) identifies similar patches within minutes. Structurally related sites are recognized even if only parts of the template patches are structurally related to the interface region. A successful prediction of the protein complex depends on a suitable template of the library. However, the performed tests indicate that interaction sites are identified even if the similarity is very low. The approach complements existing ab initio methods and provides valuable results on standard benchmark sets.     

A preference-based free-energy parameterization of enzyme-inhibitor binding. Applications to HIV-1-protease inhibitor design.

The interface between protein receptor-ligand complexes has been studied with respect to their binary interatomic interactions. Crystal structure data have been used to catalogue surfaces buried by atoms from each member of a bound complex and determine a statistical preference for pairs of amino-acid atoms. A simple free energy model of the receptor-ligand system is constructed from these atom-atom preferences and used to assess the energetic importance of interfacial interactions. The free energy approximation of binding strength in this model has a reliability of about +/- 1.5 kcal/mol, despite limited knowledge of the unbound states. The main utility of such a scheme lies in the identification of important stabilizing atomic interactions across the receptor-ligand interface. Thus, apart from an overall hydrophobic attraction (Young L, Jernigan RL, Covell DG, 1994, Protein Sci 3:717-729), a rich variety of specific interactions is observed. An analysis of 10 HIV-1 protease inhibitor complexes is presented that reveals a common binding motif comprised of energetically important contacts with a rather limited set of atoms. Design improvements to existing HIV-1 protease inhibitors are explored based on a detailed analysis of this binding motif.     

PepBind: A Comprehensive Database and Computational Tool for Analysis of Protein-peptide Interactions

Protein–peptide interactions, where one partner is a globular protein (domain) and the other is a flexible linear peptide, are key components of cellular processes predominantly in signaling and regulatory networks, hence are prime targets for drug design. To derive the details of the protein–peptide interaction mechanism is often a cumbersome task, though it can be made easier with the availability of specific databases and tools. The Peptide Binding Protein Database (PepBind) is a curated and searchable repository of the structures, sequences and experimental observations of 3100 protein–peptide complexes. The web interface contains a computational tool, protein inter-chain interaction (PICI), for computing several types of weak or strong interactions at the protein–peptide interaction interface and visualizing the identified interactions between residues in Jmol viewer. This initial database release focuses on providing protein–peptide interface information along with structure and sequence information for protein–peptide complexes deposited in the Protein Data Bank (PDB). Structures in PepBind are classified based on their cellular activity. More than 40% of the structures in the database are found to be involved in different regulatory pathways and nearly 20% in the immune system. These data indicate the importance of protein–peptide complexes in the regulation of cellular processes. PepBind is freely accessible at http://pepbind.bicpu.edu.in/.     

Preorientation of protein and RNA just before contacting

Protein and RNA molecules interact and form complexes in many biological processes. However, it is still unclear how they can find the correct docking direction before forming complex. In this paper, we study preorientation of RNA and protein separated at a distance of 5–7?Å just before they form contacts and interact with each other only through pure electrostatic interaction when neglecting the influence of other molecules and complicated environment. Since geometric complementary has no meaning at such a distance, this is not a docking problem and so the conventional docking methods, like FTDock, are inapplicable. However, like the usual docking problem, we need to sample all the positions and orientations of RNA surrounding the protein to find the lowest energy orientations between RNA and protein. Therefore, we propose a long-range electrostatic docking-like method using Fast Fourier Transform-based sampling, LEDock, to study this problem. Our results show that the electrostatically induced orientations between RNA and protein at a distance of 5–7?Å are very different from the random ones and are much closer to those in their native complexes. Meanwhile, electrostatic funnels are found around the RNA-binding sites of the proteins in 62 out of 78 bound protein–RNA complexes. We also tried to use LEDock to find RNA-binding residues and it seems to perform slightly better than BindN Server for 23 unbound protein–RNA complexes.     

Statistical analysis of structural determinants for protein–DNA‐binding specificity

DNA‐binding proteins play critical roles in biological processes including gene expression, DNA packaging and DNA repair. They bind to DNA target sequences with different degrees of binding specificity, ranging from highly specific (HS) to nonspecific (NS). Alterations of DNA‐binding specificity, due to either genetic variation or somatic mutations, can lead to various diseases. In this study, a comparative analysis of protein–DNA complex structures was carried out to investigate the structural features that contribute to binding specificity. Protein–DNA complexes were grouped into three general classes based on degrees of binding specificity: HS, multispecific (MS), and NS. Our results show a clear trend of structural features among the three classes, including amino acid binding propensities, simple and complex hydrogen bonds, major/minor groove and base contacts, and DNA shape. We found that aspartate is enriched in HS DNA binding proteins and predominately binds to a cytosine through a single hydrogen bond or two consecutive cytosines through bidentate hydrogen bonds. Aromatic residues, histidine and tyrosine, are highly enriched in the HS and MS groups and may contribute to specific binding through different mechanisms. To further investigate the role of protein flexibility in specific protein–DNA recognition, we analyzed the conformational changes between the bound and unbound states of DNA‐binding proteins and structural variations. The results indicate that HS and MS DNA‐binding domains have larger conformational changes upon DNA‐binding and larger degree of flexibility in both bound and unbound states. Proteins 2016; 84:1147–1161. © 2016 Wiley Periodicals, Inc.     

Species-dependent stereoselective drug binding to albumin: a circular dichroism study

Drug binding to albumins from different mammalian species was investigated to disclose evidence of species-dependent stereoselectivity in drug-binding processes and affinities. This aspect is important for evaluating the reliability of extrapolating distribution data among species. The circular dichroism (CD) signal induced by drug binding to the albumins [human serum albumin (HSA), bovine serum albumin (BSA), rat serum albumin (RSA), and dog serum albumin (DSA)] were measured and analyzed. The binding of selected drugs and metabolites to HSA significantly differed from the binding to the other albumins in terms of affinity and conformation of the bound ligands. In particular, phenylbutazone, a marker of site one on HSA, showed a higher affinity for binding to BSA with respect to RSA, HSA, and DSA, respectively. In the case of diazepam, a marker of site two on HSA, the affinity decreased in order from HSA to DSA, RSA, and BSA. The induced CD spectra were similar in terms of energy and band signs, suggesting almost the same conformation for the bound drug to the different albumins. Stereoselectivity was high for the binding of ketoprofen to HSA and RSA. A different sign was observed for the CD spectra induced by the drug to the two albumins because of the prevalence of a different conformation of the bound drug. Interestingly, the same induced CD spectra were obtained using either the racemic form or the (S)-enantiomer. Finally, significant differences were observed in the affinity of bilirubin, being highest for BSA, then decreasing for RSA, HSA, and DSA. A more complex conformational equilibrium was observed for bound bilirubin.     

Locating the binding sites of Pb(II) ion with human and bovine serum albumins

Lead is a potent environmental toxin that has accumulated above its natural level as a result of human activity. Pb cation shows major affinity towards protein complexation and it has been used as modulator of protein-membrane interactions. We located the binding sites of Pb(II) with human serum (HSA) and bovine serum albumins (BSA) at physiological conditions, using constant protein concentration and various Pb contents. FTIR, UV-visible, CD, fluorescence and X-ray photoelectron spectroscopic (XPS) methods were used to analyse Pb binding sites, the binding constant and the effect of metal ion complexation on HSA and BSA stability and conformations. Structural analysis showed that Pb binds strongly to HSA and BSA via hydrophilic contacts with overall binding constants of K(Pb-HSA)?=?8.2 (±0.8)×10(4) M(-1) and K(Pb-BSA)?=?7.5 (±0.7)×10(4) M(-1). The number of bound Pb cation per protein is 0.7 per HSA and BSA complexes. XPS located the binding sites of Pb cation with protein N and O atoms. Pb complexation alters protein conformation by a major reduction of α-helix from 57% (free HSA) to 48% (metal-complex) and 63% (free BSA) to 52% (metal-complex) inducing a partial protein destabilization.     

A non‐redundant protein–RNA docking benchmark version 2.0

We present an updated version of the protein–RNA docking benchmark, which we first published four years back. The non‐redundant protein–RNA docking benchmark version 2.0 consists of 126 test cases, a threefold increase in number compared to its previous version. The present version consists of 21 unbound–unbound cases, of which, in 12 cases, the unbound RNAs are taken from another complex. It also consists of 95 unbound–bound cases where only the protein is available in the unbound state. Besides, we introduce 10 new bound–unbound cases where only the RNA is found in the unbound state. Based on the degree of conformational change of the interface residues upon complex formation the benchmark is classified into 72 rigid‐body cases, 25 semiflexible cases and 19 full flexible cases. It also covers a wide range of conformational flexibility including small side chain movement to large domain swapping in protein structures as well as flipping and restacking in RNA bases. This benchmark should provide the docking community with more test cases for evaluating rigid‐body as well as flexible docking algorithms. Besides, it will also facilitate the development of new algorithms that require large number of training set. The protein–RNA docking benchmark version 2.0 can be freely downloaded from http://www.csb.iitkgp.ernet.in/applications/PRDBv2 . Proteins 2017; 85:256–267. © 2016 Wiley Periodicals, Inc.     

Carbon 13 NMR studies of saturated fatty acids bound to bovine serum albumin. I. The filling of individual fatty acid binding sites

13C NMR chemical shift and intensity results for a series of carboxyl 13C-enriched saturated fatty acids (8-18 carbons) bound to bovine serum albumin (BSA) are presented as a function of increasing fatty acid (FA)/BSA mole ratio. Spectra for long-chain (greater than or equal to 12 carbons) FA X BSA complexes exhibited up to five FA carboxyl resonances, designated a, b, b', c, and d. Only three resonances (peaks b, b', and d) were observed below 3:1 FA X BSA mole ratio, and at greater than or equal to 3:1 mole ratio, two additional resonances were observed (peaks c and a). In a spectrum of 5:1 stearic acid X BSA complexes, peaks b, b', and d each represented approximately one-fifth, and peak c approximately two-fifths, of the total FA carboxyl intensity. Plots of total carboxyl/carbonyl intensity ratio as a function of FA X BSA mole ratio were linear up to 7-9 mole ratio. Deviation from linearity at mole ratios greater than or equal to 7 was accompanied by the detection of crystalline unbound FA (as 1:1 acid/soap) by X-ray diffraction. In contrast to long-chain FA X BSA complexes, 13C NMR spectra of octanoic acid X BSA complexes yielded only one FA carboxyl resonance (peak c) at FA X BSA mole ratios between 1 and 20. We conclude: peaks b, b', and d represent FA bound to three individual high affinity (primary) long-chain FA binding sites on BSA; peak c represents FA bound to several secondary long-chain (or primary short-chain) FA binding sites on BSA; peak a represents long-chain FA bound to an additional lower affinity binding site. We present a model that correlates the observed 13C NMR resonances with individual binding site locations predicted by a recent three-dimensional model of BSA.     

A dissection of the protein-protein interfaces in icosahedral virus capsids

We selected 49 icosahedral virus capsids whose crystal structures are reported in the Protein Data Bank. They belong to the T=1, T=3, pseudo T=3 and other lattice types. We identified in them 779 unique interfaces between pairs of subunits, all repeated by icosahedral symmetry. We analyzed the geometric and physical chemical properties of these interfaces and compared with interfaces in protein-protein complexes and homodimeric proteins, and with crystal packing contacts. The capsids contain one to 16 subunits implicated in three to 66 unique interfaces. Each subunit loses 40-60% of its accessible surface in contacts with an average of 8.5 neighbors. Many of the interfaces are very large with a buried surface area (BSA) that can exceed 10,000 A(2), yet 39% are small with a BSA<800 A(2) comparable to crystal packing contacts. Pairwise capsid interfaces overlap, so that one-third of the residues are part of more than one interface. Those with a BSA>800 A(2) resemble homodimer interfaces in their chemical composition. Relative to the protein surface, they are non-polar, enriched in aliphatic residues and depleted of charged residues, but not of neutral polar residues. They contain one H-bond per about 200 A(2) BSA. Small capsid interfaces (BSA<800 A(2)) are only slightly more polar. They have a similar amino acid composition, but they bury fewer atoms and contain fewer H-bonds for their size. Geometric parameters that estimate the quality of the atomic packing suggest that the small capsid interfaces are loosely packed like crystal packing contacts, whereas the larger interfaces are close-packed as in protein-protein complexes and homodimers. We discuss implications of these findings on the mechanism of capsid assembly, assuming that the larger interfaces form first to yield stable oligomeric species (capsomeres), and that medium-size interfaces allow the stepwise addition of capsomeres to build larger intermediates.     

Solvated protein–DNA docking using HADDOCK

Interfacial water molecules play an important role in many aspects of protein–DNA specificity and recognition. Yet they have been mostly neglected in the computational modeling of these complexes. We present here a solvated docking protocol that allows explicit inclusion of water molecules in the docking of protein–DNA complexes and demonstrate its feasibility on a benchmark of 30 high-resolution protein–DNA complexes containing crystallographically-determined water molecules at their interfaces. Our protocol is capable of reproducing the solvation pattern at the interface and recovers hydrogen-bonded water-mediated contacts in many of the benchmark cases. Solvated docking leads to an overall improvement in the quality of the generated protein–DNA models for cases with limited conformational change of the partners upon complex formation. The applicability of this approach is demonstrated on real cases by docking a representative set of 6 complexes using unbound protein coordinates, model-built DNA and knowledge-based restraints. As HADDOCK supports the inclusion of a variety of NMR restraints, solvated docking is also applicable for NMR-based structure calculations of protein–DNA complexes.     

Energy landscape and transition state of protein-protein association

Formation of a stereospecific protein complex is favored by specific interactions between two proteins but disfavored by the loss of translational and rotational freedom. Echoing the protein folding process, we have previously proposed a transition state for protein-protein association. Here we clarify the specification of the transition state by working with two types of toy models for protein association. A “hemisphere” model consists of two matching hemispheres as associating proteins, and a “crater” model consists of a spherical protein with a crater to which another spherical protein fits snugly. Short-range pairwise interactions between sites across the interface hold together the bound complex. Small relative translation and rotation between the subunits quickly destroy the interactions, leading to a sharp transition between the bound state with numerous short-range interactions but restricted translation and rotational freedom and the unbound state with, at most, a small number of interactions but expanded configurational freedom. This transition sets the outer boundary of the bound state as well as the transition state for association. The energy landscape is funnel-like, with the deep well of the bound state surrounded by a broad shallow basin. Calculations with the toy models suggest that mutational change in the interaction energy in the x-ray structure of a protein-protein complex, commonly used to approximate the effect on the association constant, overestimates the effect of mutation by 10–20%. With an eye toward specifying the transition states of actual protein complexes, we find that the total number of contacts between the subunits serves as a good surrogate of the interaction energy. This formalism of protein-protein association is applied to the barnase-barstar complex, reproducing the experimental results for the association rate over a wide range of ionic strength.     

Molecular Dynamics on Complexes of ras-p21 and Its Inhibitor Protein, rap-1A, Bound to the ras-Binding Domain of the raf-p74 Protein: Identification of Effector Domains in the raf Protein

We have computed the average structures for the ras-p21 protein and its strongly homologous inhibitor protein, rap-1A, bound to the ras-binding domain (RBD) of the raf protein, using molecular dynamics. Our purpose is to determine the differences in structure between these complexes that would result in no mitogenic activity of rap-1A-RBD but full activity of p21-RBD. We find that despite the similarities of the starting structures for both complexes, the average structures differ considerably, indicating that these two proteins do not interact in the same way with this vital target protein. p21 does not undergo major changes in conformation when bound to the RBD, while rap-1 A undergoes significant changes in structure on binding to the RBD, especially in the critical region around residue 61. The p21 and rap-1A make substantially different contacts with the RBD. For example, the loop region from residues 55–71 of rap-la makes extensive hydrogen-bond contacts with the RBD, while the same residues of p21 do not. Comparison of the structures of the RBD in both complexes reveals that it undergoes considerable changes in structure when its structure bound to p21 is compared with that bound to rap-1A. These changes in structure are due to displacements of regular structure (e.g., α-helices and β-sheets) rather than to changes in the specific conformations of the segments themselves. Three regions of the RBD have been found to differ significantly from one another in the two complexes: the binding interface between the two proteins at residues 60 and 70, the region around residues 105–106, and 118–120. These regions may constitute effector domains of the RBD whose conformations determine whether or not mitogenic signal transduction will occur.     

Protein-DNA interactions: A structural analysis

A detailed analysis of the DNA-binding sites of 26 proteins is presented using data from the Nucleic Acid Database (NDB) and the Protein Data Bank (PDB). Chemical and physical properties of the protein-DNA interface, such as polarity, size, shape, and packing, were analysed. The DNA-binding sites shared common features, comprising many discontinuous sequence segments forming hydrophilic surfaces capable of direct and water-mediated hydrogen bonds. These interface sites were compared to those of protein-protein binding sites, revealing them to be more polar, with many more intermolecular hydrogen bonds and buried water molecules than the protein-protein interface sites. By looking at the number and positioning of protein residue-DNA base interactions in a series of interaction footprints, three modes of DNA binding were identified (single-headed, double-headed and enveloping). Six of the eight enzymes in the data set bound in the enveloping mode, with the protein presenting a large interface area effectively wrapped around the DNA.A comparison of structural parameters of the DNA revealed that some values for the bound DNA (including twist, slide and roll) were intermediate of those observed for the unbound B-DNA and A-DNA. The distortion of bound DNA was evaluated by calculating a root-mean-square deviation on fitting to a canonical B-DNA structure. Major distortions were commonly caused by specific kinks in the DNA sequence, some resulting in the overall bending of the helix. The helix bending affected the dimensions of the grooves in the DNA, allowing the binding of protein elements that would otherwise be unable to make contact. From this structural analysis a preliminary set of rules that govern the bending of the DNA in protein-DNA complexes, are proposed.