BACKGROUND
motions play an essential role in catalysis and protein-ligand interactions.
in particular, hinge bending motions account for 45% of motions in a representative set from the database of macromolecular motions <cit>  comprising domain hinge motions  and fragment hinge motions   <cit> .
thus understanding fundamental aspects of hinge bending mechanisms may lead to an improved understanding of the relationship between structure and function.
there are three levels of hinge prediction.
the easiest case occurs when the atomic coordinates are available for two or more conformations of a given protein.
in this case it is possible to visually inspect the motion to determine the hinge location, as we have done here.
the process can also be automated with various available packages, including flexprot <cit> , hingefind <cit>  and dyndom <cit> .
a much more difficult problem is that of predicting hinges when only one set of structural atomic coordinates is available.
several algorithms have been developed for this purpose  <cit> .
the very hardest case occurs when the sequence is known but no atomic coordinates are available at all.
the problem of finding flexible hinges between rigid regions based on sequence is in some ways similar to the problem of finding domain boundaries, which can be flexible or inflexible.
although little work has been done on the former problem, several algorithms exist to address the latter.
in one significant contribution, nagarajan and yona <cit>  analyzed multiple sequence alignments and were able to identify domains with some accuracy.
marsden et al <cit>  focused on the case of proteins with no significant sequence homology to well characterized proteins and found that predicted secondary structure contained information about domain boundaries.
jones et al.
combined puu <cit> , detective <cit> , and domak <cit>  to make a consensus-based domain boundary predictor <cit> .
heger et al.
<cit>  created the automatic domain decomposition algorithm  <cit>  and associated online database <cit> .
murzin et al.
created the scop  database <cit> .
for purposes of the above algorithms and classifications, however, domains are defined as proteins or regions of proteins having a common evolutionary origin <cit> ; flexibility is not a consideration.
indeed most small and medium sized proteins, such as those prevalent in the hinge atlas, consist of a single domain <cit> .
therefore the problem of finding flexible hinges is not solved by finding domain boundaries as defined for these methods.
schlessinger et al <cit>  developed a method to predict b-factors from sequence, but it is not clear that b-factors obtained in this way would yield accurate flexibility predictions <cit> .
in light of the limitations of existing methods, the prediction of domain hinges from sequence is considered an open problem <cit> .
in this article we focus on the characterization of these hinges based on sequence.
to that end, we compiled the hinge atlas, a manually annotated dataset of hinge bending motions, as well as a separate computer annotated dataset, both available for further studies.
the hinge atlas has several applications.
first, the statistical properties of hinges can be studied .
second, it can be used to benchmark hinge prediction programs.
third, by homology hinge annotations could potentially be transferred to proteins where the existence and location of a hinge are unknown.
fourth, the annotations could conceivably be used in future protein motion prediction programs.
the first application was of most interest to us in the current work.
our molecular motions database serves a wide variety of purposes, helping investigators understand the motion characteristics of individual proteins, as well as statistical properties of large groups of motions.
it is the ideal platform for the current study, since it contains over  <dig> morphs.
a morph is a set of atomic coordinates for two homologous protein structures , plus several structures which our morph server generates as interpolations between the two <cit> .
our server displays these structures in succession as a "movie" which suggests a possible trajectory of motion between the two conformations.
in this study we compiled two representative sets of morphs with hinge annotation: a computer annotated set and a manually annotated set, the hinge atlas.
using mostly the latter, we addressed the following questions:
<dig>  are certain residue types differentially represented in hinges?
<dig>  do certain pairs of amino acids coincide with hinges?
<dig>  can sequence be used to predict hinges?
<dig>  do hinges coincide with active sites?
<dig>  do hinges prefer certain secondary structural elements?
<dig>  do hinge residues share physicochemical or steric properties?
<dig>  are hinge residues conserved in evolution?
as our first task, we computed the rate of occurrence of each residue type in the hinge atlas.
certain amino acids were found to be differentially represented in hinges in a statistically significant fashion.
we also investigated whether certain consecutive pairs of residues were differentially represented in hinges.
in the course of the above, we observed that one of the overrepresented residues  is potentially catalytic; this was the original motivation for question  <dig> above.
to answer that question, we searched the catalytic sites atlas  <cit>  for close homologs to the proteins in our dataset, and extracted the active site residue numbers from those proteins for comparison to the hinge atlas annotation.
our next task was to investigate hinge coincidence with secondary structure.
hinges are generally believed to occur in disordered regions, but this belief has never been tested or quantified rigorously to our knowledge.
following up on our finding that hinges coincide with active site residues, we went on to the question, are hinge residues more likely to be conserved than other residues, as active sites are?
we ranked the residues by relative conservation and examined the differences between hinge and non-hinge residues.
significant correlations between sequence features and hinges were found in the above analyses.
we computed hinge indices for each of these which may be used to relate sequence features to flexibility.
we then sought to determine what predictive value sequence might have on its own and whether various sequence features collectively could be used for prediction.
we first made a simple gor   <cit> -like predictor.
we computed the log-odds rate of occurrence for residues located at the - <dig> to + <dig> positions along the sequence in the training set.
we used this table to make predictions on the test set and examined their predictive power.
as a second approach, we made a composite hinge index, which we call hingeseq, from the hinge indices of each of the sequence features found to be the strongest indicators of flexibility.
the statistical significance of this measure was computed much as for the individual sequence features.
to show that the measure is predictive, we again divided the hinge atlas into training and test sets and recomputed the relevant hinge indices to include only training set data.
we used the regenerated hingeseq to predict hinges in the test set and generated a receiver operating characteristic  curve.
as a final step, we examined molmovdb as a whole to determine whether any particular database bias was in evidence.
we also used resampling <cit>  to check for sampling artifacts in the hinge atlas.
lastly, we compared the hinge atlas to our computer annotated dataset.
the resulting work provides insight into the composition, physicochemical properties, geometry, and evolution of hinge regions in proteins.
methods
preparation of computer annotated hinge dataset
prior to generating the manually annotated hinge atlas, we used computational methods to generate a dataset of hinge residues for our statistical studies.
we began by running flexprot <cit> , a leading hinge identification tool, on all morphs  in the database of macromolecular motions <cit>  flexprot works by matching and structurally aligning fragments in one structure with corresponding fragments in the other.
the goal is to find fragment pairs which  have minimal rmsd and  are maximal in size.
the hinges are then reported as the boundaries separating those fragments.
goal  is equivalent to minimizing the number of these hinges.
since domains are never completely rigid, rmsd tends to grow with fragment size and therefore goal  is in conflict with goal .
this conflict is dealt with by providing the user with a series of adjustable parameters, and further by reporting not one but several alternative hinge locations from which the user can choose.
we used a combination of computer and manual culling to select those morphs for which the identified hinges met the following criteria:
<dig>  motion was domain wise, i.e.
two or more domains could be observed moving approximately as rigid bodies with respect to each other.
<dig>  the identified hinge was located in the flexible region connecting two rigid domains, rather than in the domains themselves.
<dig>  the morph trajectory was sterically reasonable, i.e.
chains were not broken in the attempt to interpolate motion.
we found that flexprot's maximal <cit>  rmsd  parameter had a strong effect on the results.
therefore when flexprot gave visibly incorrect results for a given morph, we reran the program, systematically varying this parameter.
if one of these runs gave sufficiently accurate results, the annotation for that morph was entered into the database.
we discarded immediately those morphs that did not exhibit clear hinge bending motion.
lastly, we removed redundant morphs using nrdb <dig> <cit> .
note that the definition of a hinge given in the introduction allows for a hinge of zero length.
flexprot indeed often returned such hinges.
to deal with this, in all cases one residue on each side of the hinge, was taken to also belong to the hinge.
thus most hinges are two residues long.
at the end of this process, the computer annotated set contained  <dig> morphs.
as described, the computer annotation of hinges requires significant human intervention and the results were often debatable.
many of the hinge annotations differed slightly but visibly from the boundary between rigid domains, such that the backbone flexions that could account for the domain motion were not seen in the predicted hinge region.
in other cases hinges were missed, and some annotations appeared where no hinge existed.
the more flagrantly misannotated hinges were removed from the dataset, but making the manual culling too stringent would simply have resulted in a dataset too small to be statistically meaningful.
for these reasons, the computer annotated dataset was not used in most of this work.
nonetheless, the computer annotated dataset is arguably more objective then the manually annotated set described below, and so is made available to the community.
to address the accuracy issues, we decided to generate a manually annotated set of hinges – the hinge atlas.
to generate this set we first created the hinge annotation tool which can also be used by the public as we will now explain.
the hinge annotation tool
the creation of publicly accessible tools for manual annotation of hinges involved significant changes to the morph page.
the morph page is the primary point on molmovdb <cit>  for analyzing single morphs.
it is accessible from the "movies" page or through our search tool, both linked to or visible on our front page.
our server also provides a link to this page in an email sent to the submitter of each morph request.
we added all of the new tools to the "hinge analysis" tab on this page.
the first of these is the hinge annotation tool.
each of three rows of "arrow" buttons on this tool move a highlighted window of two residues along the protein chain, allowing the user to highlight up to three hinges in a protein.
the "show all" button then highlights all selected residues in the jmol viewer window.
once the user is satisfied with the hinge selection, clicking "submit" records this selection in the database.
once the morph page is regenerated, a "show public hinge" button will be visible which, when clicked, highlights the selected residues.
lastly, the user can use a pointing device to reorient the protein in the jmol window to his/her liking.
a gif image based on that view can be generated by clicking on the "color by domain" link.
the animation will be rendered using vmd's <cit>  "new cartoon" style, with the identified hinge region and two rigid domains each colored distinctively.
the hinge annotations made in this way persist in our database for visualization and use by others, until overwritten.
with minor modification, these tools were used to generate the hinge atlas dataset of manually annotated hinges.
the criteria we used for selection are described in the following section.
highlighting the hinge atlas hinges  on the animated morph movie is a matter of going to the morph page and clicking on the "hinge analysis" tab as above and clicking the "show hinge atlas hinge" button.
the annotated hinge location will be rendered in green spacefill style, which contrasts with the white trace used elsewhere in the protein.
construction of the hinge atlas
the tools described above answer only the technical question of how we annotated hinges.
in this section we clarify the motivation for the hinge atlas and its applications and answer the scientific question of how we decided on the precise location of the hinge for each morph.
for each morph in the hinge atlas, we used the hinge annotation tool as described to select the hinge location.
motivated in part by our long term goal of providing a resource that could be used in motion prediction work, and in part by a desire to deepen basic understanding of protein motion, we asked ourselves the following question:
would it be possible to approximately reproduce the observed motion by allowing flexure at the hinge points but keeping the regions between hinges rigid?
in order for this question to be answered in the affirmative, the hinge selection should be the one to best meet the following criteria:
<dig>  the φ, ψ, and α  <cit>  torsion angles of hinge residues may often  be larger than those of their neighbors.
<dig>  amino acids on either side of the hinge residues must be co-moving with their respective rigid regions.
<dig>  rotations of one of the rigid regions about the hinge region must not result in significant and irreconcilable steric clashes.
in order to use  as a useful guide to selecting the hinge location, we made use of the torsion angle charts and graphs in the structure analysis tools section on the morph page.
however often large rotations of the main chain are induced by multiple cooperative torsions in the hinge, and these may be individually small, particularly in α-helices <cit> .
the usefulness of this flexibility measure is further limited by the frequent occurrence of large torsion angles which do not coincide with hinges <cit> .
nonetheless, when the precise location of the hinge was otherwise unclear, torsion angles were often examined to help adjust the selection.
criterion  is a definition of a hinge.
sometimes the hinge was slightly longer than others, and in those cases we added more residues to the hinge, up to a limit of about five residues in total.
if the hinge was distributed over too many residues such that no one short stretch could be said to constitute the entire hinge, then the morph was discarded from the hinge atlas, since the motion was not hingelike.
criterion  is a practical requirement of a working hinge.
if substantial flexure at points outside the hinge is required to avoid domain interpenetration, then the choice of hinge location is incorrect, or the motion is not hinge but rather shear or unclassifiable <cit> .
the next question was, how to select the morphs which would be annotated and included in the hinge atlas.
the entire database of macromolecular motions  <cit>  with  over  <dig> morphs, could clearly not all be annotated given limited manpower.
further, only a minority of morphs  exhibited hinge bending motion, and even within this group much redundancy existed.
to address these issues and make the annotation work manageable, we first selected a nonredundant subset of the morphs in molmovdb by aligning all sequences to nrdb <dig> <cit> .
this reduced the dataset to  <dig> morphs.
this was more manageable, but still the set contained many proteins which did not exhibit hinge bending motions.
fortunately we found that the score output by flexprot, normalized by dividing by the number of residues, provided an accurate measure of the degree to which a protein exhibited hinge bending.
high scores, close to unity, indicated proteins more likely to exhibit hinge bending motion.
lower scores, below  <dig>  or so, were very unlikely to do so.
we sorted the  <dig> nonredundant morphs by descending normalized flexprot score  and annotated them in that order.
those proteins for which we could find hinges allowing a positive answer to the question above were annotated and added to the hinge atlas.
those proteins which did not exhibit hinge bending motion or for which no suitable hinge could be found were discarded.
at the end of this culling and annotation effort, the hinge atlas contained  <dig> nonredundant annotated morphs.
we also manually annotated a small set of specifically fragment  hinge bending motions which may be useful for some studies, described below.
availability of datasets
in the course of this study we compiled a number of sets of morphs which can be viewed on our online galleries listed and linked to on our sets page <cit> .
the hinge atlas and computer annotated sets are compared more rigorously in the "statistical comparison of datasets" section.
the galleries provide easy browsing and visual inspection of morph movies sharing certain characteristics.
the sets offered include:
nonredundant
no two morphs in this set have more than 90% sequence homology.
this set was compiled by alignment to proteins in nrdb <dig>
catalytic sites atlas
all morphs in this set have annotated active sites which can be highlighted in the jmol viewer.
catalytic sites atlas
same as above, but with redundant morphs removed by comparison to nrdb <dig>
flexprot hinges
computer annotated set used in parts of this study and described above.
we consider it to be less useful than the hinge atlas, but the data is nonetheless made available.
fragment hinge motions
a small set of hinge bending motions involving fragments smaller than domains, as alluded to in the previous section.
hinge atlas
contains the manually annotated protein pairs used in this study.
a link on the sets page permits the download of the sequence data  in mysql format.
the same data is available in tab-delimited text format which is human readable and importable into ms excel and other packages.
another link on the same page facilitates the download of the interpolated structure files associated with each morph in the hinge atlas set.
clicking on the thumbnail image leads to the "movies" page, where users can browse through the  <dig> proteins in the hinge atlas.
clicking on any of the protein thumbnail images, in turn, leads to the corresponding morph page, where the hinge annotation can be viewed as described in the "hinge annotation tool" section
method for analyzing relative frequency of occurrence
throughout this study, we will be comparing how often a particular entity  occurs in hinges versus everywhere in the hinge atlas or another of the datasets described above.
the statistical analysis will be the same regardless of the particulars, so we will here present the general approach and later only mention adjustments particular to the specific question addressed.
first we defined the following variables:
d = total number of residues in the dataset
h = total number of residues in hinges in the dataset
c = classification scheme used to create groups of residue positions.
for example, c could be secondary structure, degree of conservation, etc.
c = a particular grouping of residues, where c ∈ c. for instance, if c = secondary structure, then c = helix is the class of all residues in helices, c = strand is the class of all residues in strands, etc.
another example might be c = evolutionary conservation, with c = cons <dig> = top 20% most conserved residues, c = cons <dig> = second 20% most conserved, etc.
ac = set of all residues of class c in the dataset.
dc = number of times residues of class c occurred anywhere in the dataset.
hc = number of times residues of a particular class c occurred in hinges.
these can be used to estimate various probabilities as follows:
p = dc/d is the prior probability of c – in other words, the probability that residues of class c occur anywhere in the dataset.
p = hc/h is the conditional probability that a residue belongs to class c, given it is a hinge.
a quantity that is of interest in hinge prediction is the posterior probability p, the probability that a residue is a hinge given it is in ac.
we obtain this from bayes' rule:
p=p⋅pp=hcdc
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgwbaccqggoaakcqwgobaacqgg8bafcqwghbqydawgaawcbagaem4yamgabeaakiabcmcapiabg2da9maalaaabagaemicaanaeiikagiaemyyae2aasbaasqaaiabdogajbqabagccqgg8bafcqwgobaacqggpaqkcqghfly1cqwgwbaccqggoaakcqwgobaacqggpaqkaeaacqwgwbaccqggoaakcqwghbqydawgaawcbagaem4yamgabeaakiabcmcapaaacqgh9aqpdawcaaqaaiabdigaonaabaaaleaacqwgjbwyaeqaaagcbagaemizaq2aasbaasqaaiabdogajbqabaaaaaaa@52b3@
equation 1
where the prior probability that a residue is a hinge is given by p=
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgwbaccqggoaakcqwgobaacqggpaqkcqgh9aqpcqggoaakdawcaaqaaiabdieaibqaaiabdseaebaacqggpaqkaaa@3612@.
we further define the hinge index hi, similar to the domain linker index used in armadillo <cit> :
hi=log⁡10pp
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgibascqwgjbqscqggoaakcqwghbqydawgaawcbagaem4yamgabeaakiabcmcapiabg2da9igbcygasjabc+gavjabcegannaabaaaleaacqaixaqmcqaiwaamaeqaaowaasaaaeaacqwgwbaccqggoaakcqwghbqydawgaawcbagaem4yamgabeaakiabcyha8jabdigaojabcmcapaqaaiabdchawjabcicaoiabdggahnaabaaaleaacqwgjbwyaeqaaogaeiykakcaaaaa@4959@
equation 2
the argument of the log is the ratio of the observed frequency of occurrence of classes of amino acids ac in hinges, over the expected.
note that this argument is close to the likelihood ratio pp
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaadawcaaqaaiabdchawjabcicaoiabdggahnaabaaaleaacqwgjbwyaeqaaogaeiifawnaemiaagmaeiykakcabagaemicaanaeiikagiaemyyae2aasbaasqaaiabdogajbqabagccqgg8bafcqgg+bgfcqwgobaacqggpaqkaaaaaa@3fc8@ used in bayesian statistics <cit>  because h is so small compared to d. the quantity hi yields an intuitive measure of the enrichment of certain classes of residues in hinges, with positive numbers indicating enrichment and negative numbers indicating scarcity.
just because the hi is nonzero, however, does not mean that the differential representation has statistical significance.
to establish the latter, we considered two statistical hypotheses:
h0: the null hypothesis.
assume hc is a randomly distributed random variable with mean μh.
the null hypothesis states that:
μh=dcd⋅h
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaaiigacqwf8oqbdawgaawcbagaemiaaggabeaakiabg2da9maalaaabagaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwgebaraagaeyyxictaemisageaaa@3858@
if this is true, then the hinge set is chosen without replacement in an unbiased fashion from the dataset, and p is given by the hypergeometric distribution .
h <dig>  the alternate hypothesis.
this states that:
μh≠dcd⋅h.
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaaiigacqwf8oqbdawgaawcbagaemiaaggabeaakiabgcmi5oaalaaabagaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwgebaraagaeyyxictaemisagkaeiola4caaa@39fd@
it is equivalent to saying that p is not p, and therefore the null hypothesis can be rejected.
we test this as follows.
if it is the case that
hch>dcd,
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaadawcaaqaaiabdigaonaabaaaleaacqwgjbwyaeqaaagcbagaemisageaaiabg6da+maalaaabagaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwgebaraagaeiilawcaaa@3692@
and if we choose a significance threshold of  <dig> , we can reject h <dig> iff our p-value ∑x=hh∞hyp< <dig>
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaadaaewbqaaiabdieaijabdmfazjabdcfaqjabcicaoiabdieaijabcycasiabdseaejabcycasiabdiha4jabcycasiabdsgaknaabaaaleaacqwgjbwyaeqaaogaeiykakiaeyipawjaegimaajaeiola4iaegimaajaegynaudaleaacqwg4baecqgh9aqpcqwgobaacqggoaakcqwghbqycqggpaqkdawgaaadbagaemiaaggabeaaasqaaiabg6hilcqdcqghris5aaaa@4bc9@.
the left hand side of the latter inequality can be interpreted as the probability that hc or more residues of class ac could be found in hinges, assuming h <dig> and given h, d, and dc.
the argument of the sum is the hypergeometric function, which gives the probability that dc residues taken without replacement from a set of d residues of which h are hinges, would contain exactly x hinges:
hyp
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgibascqwgzbqwcqwgqbaucqggoaakcqwgibascqggsaalcqwgebarcqggsaalcqwg4baecqggsaalcqwgkbazdawgaawcbagaem4yamgabeaakiabcmcapmaalaaabawaaewaaeaafaqaaegabaaabagaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwg4baeaaaacagloagaayzkaawaaewaaeaafaqaaegabaaabagaemiraqkaeyoei0iaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwgibascqghsislcqwg4baeaaaacagloagaayzkaaaabawaaewaaeaafaqaaegabaaabagaemiraqeabagaemisageaaagaayjkaiaawmcaaaaaaaa@4e8f@
equation 3
otherwise, if it is the case that
hch<dcd,
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaadawcaaqaaiabdigaonaabaaaleaacqwgjbwyaeqaaagcbagaemisageaaiabgyda8maalaaabagaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwgebaraagaeiilawcaaa@368e@
then we reject h <dig> iff our p-value ∑x=−∞hhyp)< <dig>
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaadaaewbqaaiabdieaijabdmfazjabdcfaqjabcicaoiabdieaijabcycasiabdseaejabcycasiabdiha4jabcycasiabdsgakjabcicaoiabdggahnaabaaaleaacqwgpbqaaeqaaogaeiykakiaeiykakiaeyipawjaegimaajaeiola4iaegimaajaegynaudaleaacqwg4baecqgh9aqpcqghsislcqgheispaeaacqwgobaacqggoaakcqwghbqycqggpaqka0gaeyyeiuoaaaa@4e2e@.
RESULTS
are certain amino acids more likely to occur in hinges?
we applied the described statistical formalism to the problem of amino acid frequency of occurrence in hinges by taking c = amino acid type, and c to designate each of the  <dig> canonical amino acids.
hi scores and p-values were thus calculated for each of  <dig> identifications of c corresponding to the  <dig> canonical amino acids.
we found that glycine and serine are overrepresented in a highly significant fashion.
we also found phenylalanine, valine, alanine, and leucine to be underrepresented, albeit with lower significance .
we also investigated the frequency of occurrence of sequential pairs of amino acids in hinges, but since  <dig> sequential pairs are possible the significance of the results was much lower and no conclusion could be drawn.
are residues within a certain distance of an active site more likely to be hinge residues?
as mentioned earlier, the fact that one of the overrepresented residues is potentially catalytic led us to suspect that hinge residues are more likely to occur in active sites, or within a few residues of an active site, than would be expected by chance.
this would make sense from a biochemical and mechanical perspective.
hinge motions are often opening and closing motions of domains intended to expose the active site, which often would be located at the center of the motion, i.e.
the hinge.
prior work <cit>  shows that active sites are more likely to occur at regions of low first normal mode displacement.
such regions have been shown to coincide with hinges <cit> .
here we close the loop, comparing active sites directly with the hinge atlas annotation and quantifying the correspondence.
in order to annotate the active site locations, we blasted <cit>  the morph sequences in the computer annotated dataset against the sequences in the catalytic sites atlas and considered a morph in the hinge dataset to match a protein in the csa if they had sequence identity ≥ 99%.
this high threshold was chosen to minimize the possibility of incorrectly labeling a residue in the hinge atlas and thereby diminishing the significance of the results.
for each such pair, we transferred the catalytic site annotation to the morph.
we described earlier how to browse the csa morphs online.
of the  <dig> proteins in the hinge atlas,  <dig> were annotated with active site information from the csa; the rest had no close csa homologs.
the  <dig> proteins comprised the dataset for this calculation.
we analyzed this set using the statistical formalism described earlier, with the following variable definitions:
c = distance from the nearest active site, in residues.
c = successively: active site residues, amino acids  <dig> residue away from the nearest active site residue,  <dig> residues away, etc.
d =  <dig> residues in the dataset of  <dig> proteins
h =  <dig> hinge residues in the dataset
dc = residues of class c in the dataset
hc = residues of class c in hinges.
the results are shown in figure  <dig> and table  <dig>  at short distances from the active sites, hinge residues were overrepresented.
the active site residues and residues as much as four residues away from the nearest active site were significantly overrepresented in hinges.
are hinges segregated by secondary structure?
it is generally accepted that hinges tend to avoid secondary structure.
however this belief has, to our knowledge, never been tested on a quantitative basis, and indeed numerous counterexamples can be found.
for instance, the hinge in calmodulin <cit>  and troponin c <cit>  occurs in an α-helix, and in glutamine binding protein it occurs in two parallel beta strands <cit> .
thus we do not know which particular types of secondary structure are avoided or preferred, or to what degree.
to obtain this information, we tabulated the number of hinge residues occurring in the various types of secondary structural elements, and compared this with the distribution of all residues, proceeding as follows.
stride <cit>  recognizes secondary structural elements from atomic coordinates.
we used this program to assign secondary structural classes to all residues in the hinge atlas.
we then tabulated the number of residues assigned to each class, both in hinges and elsewhere in the dataset.
lastly, we calculated the hi scores and the p-values as before, letting c = secondary structural element type and c designate e.g.
helix, coil, etc.
we found that three types of secondary structure were differentially represented in hinges with extremely high significance.
we conclude that hinges are less likely to occur in α-helices, and are more likely to occur in turns or random coils .
for the user's convenience, secondary structure assignments for individual morphs can be obtained in the 'hinge analysis tools' box on the molmovdb morph page <cit>  mentioned earlier.
are certain physicochemical properties preferred in hinge residues?
it is intuitive that certain physicochemical classes of residues  would occur more frequently in hinges, and this would help explain the amino acid propensities reported earlier in this work.
to check and quantify this, we grouped amino acids into several non-exlusive categories <cit> .
following again our statistical treatment, we calculated hi scores and p-values, letting c = physicochemical grouping, and c = aliphatic, polar, charged, etc.
we discovered that aliphatic and hydrophobic residues were very significantly underrepresented.
overrepresented were small and tiny residues .
are hinge residues conserved in evolution?
we next investigated whether hinge residues are conserved.
since certain residue classes are preferred in hinges, one might suspect that hinge residues would be conserved.
first, we blasted <cit>  each of the hinge atlas sequences against nrdb <dig>  a non-redundant sequence database in which protein sequences have no more than 90% sequence identity with each other  <cit>  next we extracted up to  <dig> top-aligned sequences to a given morph to generate a multiple sequence alignment using clustal w <cit> .
for each position in the multiple sequence alignment, we used the formalism developed by schneider et al <cit>  to compute the information content associated with a column in the multiple sequence alignment at this position <cit> .
we sorted the residues in hinge atlas morphs according to the magnitude of the information content scores.
we then divided the residues into five bins of equal size.
if hinge residues are conserved, then there should be an enrichment of hinge residues in the top bins, which correspond to the most conserved residues.
on the other hand, if hinge residues are hypermutable, there should be more of them in the bottom bins, corresponding to the least conserved residues.
because it is widely agreed that active sites should be conserved, we used the conservation of active sites as a control.
to quantify the enrichment, we calculated the hi scores as described previously.
here, c is a label applied to residues that ranked in a given percentile bin, e.g.
the top 20% most conserved.
for that bin p = hc/h is thus the ratio of the number of hinge residues in the bin divided by the total number of hinge residues.
similarly, p = dc/d is the ratio of the number of residues in the dataset in the bin divided by the grand total of residues in the dataset.
to determine the statistical significance of hi scores, we calculated the p-values using the hypergeometric distribution with the dc, hc, d, h defined above.
for the control set, we performed the same calculation but made the following changes to the variable definitions:
<dig>  our dataset was no longer the hinge atlas, but rather the "catalytic sites atlas " set described earlier.
d is the total number of residues in this set.
<dig>  ac still represents residues in the dataset belonging to a given conservation rank bin.
dc is the total number of residues in that bin.
<dig>  hc now represents the number of active site residues in a given bin corresponding to c. similarly, h represents the total number of active site residues in the dataset.
we found that hinge residues distribute evenly in the top 80%, and have a slight but statistically very significant enrichment in the bottom 20% bin .
thus hinge residues are hypermutable.
we observe a highly significant enrichment of active site residues in the top 20% bin, as expected.
among the  <dig> active site residues,  <dig> of them  are in the highest bin, and the numbers progressively decrease in lower bins.
the hinge atlas pools enzymes together with non-catalytic proteins.
we reasoned therefore that perhaps only hinges in non-catalytic proteins are hypermutable, and that if we analyzed a set consisting only of enzymes, then the propensity of active sites to occur in hinges would lead to conservation, rather than hypermutability of hinge residues for that set.
to test this idea, we decided to calculate the propensity of hinges to occur in specific bins of conservation score, for the  <dig> proteins in the hinge atlas with csa annotation, rather than for the larger set of  <dig>  for this set we also found that hinge residues occur more frequently among the 20% least conserved residues for each protein .
at a p-value of  <dig> , the confidence in this result is high.
even this test, however pools together hinges that are near the active site  with hinges that occur at some distance from it.
so we selected from the  <dig> proteins a small set that had at least one active site residue in the hinge, and removed the active site residues themselves.
we then calculated the propensity of hinge residues to occur in the five conservation bins.
this set was found to be too small, however, and statistical significance was too low to draw a conclusion .
a study using the set of fragment hinge motions described earlier was similarly inconclusive.
the hypermutability of hinge residues that we found is reasonable because hinge residues tend to be on the surface of proteins  rather than in the more highly conserved core.
hinges are less likely to be buried inside domains because they would then be highly coordinated with near neighbors and hence less flexible.
the apparent contradiction of hypermutability on the one hand and enrichment of active sites on the other is dealt with in the discussion section.
are hinge residues more likely to occur on the surface?
to support our argument that hinge residues are hypermutable partly because they occur on the surface, we quantified the degree to which the latter is the case.
to do this, we used a solvent accessible surface area  calculation program  <cit>  with a probe radius of  <dig> Å. the asa of each of the backbone heavy atoms  was calculated and summed for each residue in each protein in the hinge atlas.
we then binned the residues by this quantity.
lastly, we counted the number of hinge residues in each bin and calculated hi and p-value as before.
as expected, bin # <dig>  was significantly enriched with hinges .
bin # <dig> was also highly enriched, while bins # <dig> and # <dig> had fewer hinges, all with extremely high significance.
bins 1– <dig> contain  <dig> total residues and bin  <dig> contains  <dig>
how many hinge sites appear in each protein?
perhaps the simplest hinge consists of a single point on the chain separating two rigid regions.
however it is also possible for the chain to pass multiple times through the same region, or to have multiple independent hinge regions.
this leads to the question, how many proteins had single hinge points, versus a larger number of hinge points?
we answer this question in table  <dig>  most morphs had three or fewer hinge points.
can hinges be predicted by a simple gor-like method?
gor <cit>  method is useful for predicting secondary structure from sequence with fair accuracy.
we implemented a gor-like method to determine whether sequence contained enough information for hinge prediction.
we divided the dataset into a training set and a test set for this study.
the log-odds frequency of occurrence of amino acids in the training set were tabulated not only at a given hinge residue, but also at positions ranging from - <dig> to + <dig> from the given residue in sequence space.
for simplicity, hinge residues at positions less than eight residues from either end of the chain were not included.
once the table was generated, it was used on the test set.
the score for a given residue was taken to be the sum of the scores for the residues in positions - <dig> to + <dig> from that residue.
the scores were computed for all residues in the test set, except those less than eight residues from either end of the chain.
the idea is that a threshold score can be chosen and residues scoring higher than this threshold are considered more likely to be hinges.
note that where robson and suzuki used a different fitting parameter for each type of secondary structure, we used no fitting parameter, since we were interested in only one "secondary structure": the hinges.
the rates of true and false positives and negatives were calculated for each choice of score threshold over a range.
our training set numbered  <dig> proteins from the computer annotated set.
we tested the method on a test set of  <dig> proteins from the same set and obtained a roc curve <cit>  .
the area under this curve was nearly  <dig> , indicating negligible predictive value.
hinge prediction by combining sequence features
as the gor-like method did not work well, we sought to measure the predictive power of the various sequence features studied above.
the hi scores we have reported provide an intuitive means of weighing the relative predictive value of each sequence feature.
we show how to combine the hi scores for several features in order to make a more powerful predictor, which we call hingeseq.
we define this predictor as follows:
hs=log⁡10ppppp)=hiamino⋅acid+hisecondary⋅structure+hiactive⋅site
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgibascqwgtbwucqggoaakcqwgpbqacqggpaqkcqgh9aqpcyggsbabcqggvbwbcqggnbwzdawgaawcbagaegymaejaegimaadabeaakmaabmaabawaasaaaeaacqwgwbaccqggoaakcqwghbqydawgaawcbagaemoaaogabeaakiabcyha8jabdigaojabcmcapiabdchawjabcicaoiabdggahnaabaaaleaacqwgrbwaaeqaaogaeiifawnaemiaagmaeiykakiaemicaanaeiikagiaemyyae2aasbaasqaaiabdygasbqabagccqgg8bafcqwgobaacqggpaqkaeaacqwgwbaccqggoaakcqwghbqydawgaawcbagaemoaaogabeaakiabcmcapiabdchawjabcicaoiabdggahnaabaaaleaacqwgrbwaaeqaaogaeiykakiaemicaanaeiikagiaemyyae2aasbaasqaaiabdygasbqabagccqggpaqkaaaacagloagaayzkaagaeyypa0jaemisagkaemysak0aasbaasqaagqaciab=fgahjab=1gatjab=lgapjab=5gaujab=9gavjabgwsixlabdggahjabdogajjabdmgapjabdsgakbqabagccqggoaakcqwgpbqacqggpaqkcqghrawkcqwgibascqwgjbqsdawgaawcbagae83camnae8xzaumae83yammae83ba8mae8nba4mae8hzaqmae8xyaemae8ncainae8xeaknaeyyxictaem4camnaemidaqnaemocainaemydaunaem4yammaemidaqnaemydaunaemocainaemyzaugabeaakiabcicaoiabdmgapjabcmcapiabgucariabdieaijabdmeajnaabaaaleaacqwghbqycqwgjbwycqwg0badcqwgpbqacqwg2bgdcqwglbqzcqghfly1cqwgzbwccqwgpbqacqwg0badcqwglbqzaeqaaogaeiikagiaemyaakmaeiykakcaaa@b40a@
equation 4
for simplicity, statistical independence of the various features was assumed in creating this definition.
here the i's correspond to individual amino acids in the protein sequence.
for each i, j designates one of the  <dig> amino acid types, k designates the secondary structural classification, and l designates active site versus non-active site classification.
thus hiamino·acid is assigned according to residue type by looking up the corresponding value in table  <dig>  similarly, hisecondary·structure isobtained according to secondary structure type from table  <dig>  following table  <dig> approximately, we assign hiactive·site as  <dig>  for residues four or fewer amino acid positions away from the nearest active site residue, and  <dig>  elsewhere.
the highest values of hs correspond to residues most likely to occur in hinges.
clearly, extending this method is only a matter of obtaining amino acid propensities to occur in hinges according to additional classifications.
the resulting index can then simply be included as an additional term in the above formula, with no need for adjustable weighting factors.
we evaluated the statistical significance of this measure much as for the individual sequence features.
we counted the number of residues in the hinge atlas with a hingeseq score above  <dig> , and within that set the number of hinge residues.
we compared this to the total number of hinges and the population size of the hinge atlas .
using the cumulative hypergeometric distribution as before, we computed a p-value of order 10- <dig>  thus the measure shows high statistical significance.
however since only about 5% of the residues scoring over  <dig>  were annotated hinges, hingeseq is not likely to be sensitive enough to be used alone for hinge prediction.
the low p-value indicates that the predictor results have high statistical significance.
however the low sensitivity limits its potential predictive value.
we nonetheless wished to show that hingeseq is predictive, rather simply reflectling peculiarities of the dataset.
to this end, we divided the  <dig> proteins of the hinge atlas into a training set numbering  <dig> proteins, and a test set numbering  <dig>  of the  <dig> hinge atlas proteins, the  <dig> proteins with annotation from the csa were apportioned such that  <dig> were included in the training set and  <dig> in the test set.
we tested the performance of the predictor by means of roc  curves.
we need to define a few terms in order to use these:
test positives: residues with hsgreater than or equal to a certain threshold.
test negatives: residues with hsless than a certain threshold.
gold standard positives: residues annotated as hinges in the hinge atlas.
gold standard negatives: residues which are not in hinges according to the hinge atlas annotation.
true positives : those residues that are both test positives and gold standard positives.
true negatives : residues that are both test negatives and gold standard negatives.
false positives : residues that are test positives and gold standard negatives.
false negatives : residues that are test negatives and gold standard positives.
sensitivity=tptp+fn
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgzbwccqwglbqzcqwgubgbcqwgzbwccqwgpbqacqwg0badcqwgpbqacqwg2bgdcqwgpbqacqwg0badcqwg5bqecqgh9aqpdawcaaqaaiabdsfaujabdcfaqbqaaiabdsfaujabdcfaqjabgucariabdaeagjabd6eaobaaaaa@450b@
specificity=tnfp+tn
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgzbwccqwgwbaccqwglbqzcqwgjbwycqwgpbqacqwgmbgzcqwgpbqacqwgjbwycqwgpbqacqwg0badcqwg5bqecqgh9aqpdawcaaqaaiabdsfaujabd6eaobqaaiabdaeagjabdcfaqjabgucariabdsfaujabd6eaobaaaaa@44a9@
1−specificity=fpfp+tn
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqaixaqmcqghsislcqwgzbwccqwgwbaccqwglbqzcqwgjbwycqwgpbqacqwgmbgzcqwgpbqacqwgjbwycqwgpbqacqwg0badcqwg5bqecqgh9aqpdawcaaqaaiabdaeagjabdcfaqbqaaiabdaeagjabdcfaqjabgucariabdsfaujabd6eaobaaaaa@466e@
the roc curve is simply a plot of the true positive rate  vs. false positive rate , for each value of the threshold, as the threshold is varied from + <dig> to - <dig>  a range which included all possible values of hs.
for a good predictor, the true positive rate will increase faster than the false positive rate as the threshold is lowered, and the area under the curve will be significantly greater than  <dig> .
the roc curve is shown in figure  <dig>  although work remains to be done before sequence-based hinge prediction can be relied upon exclusively, hingeseq displays significant ability to detect potential for flexibility.
checking for dataset bias
these findings assume that the dataset used does not contain significant bias or artifacts, either in the composition of the entire dataset or of the hinges within it.
to substantiate this, we performed various studies as follows.
bias in amino acid composition and functional classification
in order to find out whether the molmovdb database contained any bias in amino acid composition, we extracted the sequences of all the morphs in molmovdb and counted the total occurrence of each residue type.
suspecting that redundancies might bias the result, we clustered the sequences and recounted the amino acid residues in the same way.
we compared these numbers to publicly available amino acid frequencies of occurrence for the pdb  <cit>  .
the amino acid frequency of occurrence for the clustered molmovdb morphs was found to be essentially that of the pdb, from which it was created, therefore no particular database bias is in evidence.
we also sought to determine whether there existed a bias towards particular protein classes, in either the hinge atlas or the nonredundant set of molmovdb morphs from which it was compiled.
to do this, we first counted the number of times each top-level gene ontology  term under the "molecular function" ontology was associated with a protein in the hinge atlas.
where the annotation was given for deeper levels, we traced up the hierarchical tree to retrieve the corresponding top level term in the ontology.
thus we found, for example, that  <dig> proteins in the hinge atlas were associated with the term "nucleic acid binding."
we repeated this procedure for the pdb as a whole as well as for the non-redundant set of  <dig> morphs in molmovdb from which the hinge atlas was compiled.
the results for the  <dig> most frequently encountered go terms are shown in table  <dig>
to compare the hinge atlas counts to the pdb counts in an overall fashion, we used the chi-square distribution with  <dig> degrees of freedom  and obtained a chi-square value of  <dig> .
this corresponds to a p-value of  <dig> , so there is no statistically significant difference in the distribution of these terms in the hinge atlas vs. the entire protein data bank.
statistical comparison of datasets
the hinge atlas and computer annotated sets were compiled differently, therefore one might suspect that the hinges from one set might comprise a statistically different population from the hinges of the other set.
if this were the case, then one of the two sets would be preferable to the other, otherwise if the populations were essentially the same then the two sets could potentially be used interchangeably.
it is therefore necessary to quantitatively compare these two populations.
it is also necessary to confirm that within one set, the hinge residues are a statistically distinct population from the rest of the set; if this were not true then the amino acid propensity data reported earlier would not be meaningful.
although the hinge atlas and the computer annotated set share a total of  <dig> residues, only  <dig>   of these are hinge residues.
this is another reason to suspect that the hinge population of the hinge atlas is statistically different from the hinge population of the computer annotated set.
to test this, we computed the chi-square value for hinge atlas hinges vs. computer annotated hinges, and obtained a p-value of  <dig> .
therefore, the hinge atlas hinges are different from the computer annotated hinges.
the chi-square value describing the difference between amino acid frequency of occurrence in the hinge vs. non-hinge subsets of the hinge atlas was  <dig> .
with  <dig> degrees of freedom  this corresponds to a p-value below 10- <dig> .
therefore, the hinge residues are shown with high confidence to be different from non-hinge residues in the hinge atlas.
a similar calculation yielded a p-value of  <dig>  for the computer annotated set.
therefore we conclude that the hinge and non-hinge populations are different for the computer annotated set, as well.
the same is true for the hinge atlas set.
hinges in the hinge atlas are likely to be a distinct population from hinges in the computer annotated set ; therefore for future studies one or the other should be used, and we recommend that it be the hinge atlas.
we conclude from this calculation for both the hinge atlas and computer annotated set, the hinge population is different from the non-hinge population, therefore statistically significant information can be extracted from both.
however the hinge population of the hinge atlas is different from that of the computer annotated set, albeit with much lower significance.
we argue that one of the two sets should therefore be preferred for statistical studies.
the preferred set should be the hinge atlas since the computer annotated set contains numerous annotations which are slightly different from the correct and visually verifiable hinge location.
similarity within morph pairs
we next asked the question, do the morphs in the hinge atlas reflect intrinsic flexibility of the protein, or is the apparent conformational change the result of sequence differences between the two structures in the pair?
that is, do the morphs display motions observable in a single protein, or do they instead represent evolutionary change?
to answer this we counted the number of times both structures in the morph came from the same vs different organisms.
of the  <dig> morphs,  <dig> had structures downloaded directly from the pdb rather than uploaded by users, and also had valid source organism data.
for  <dig> of the  <dig>  both proteins in the pair came from the same species, while for  <dig> the two proteins came from different species.
of the  <dig>   <dig> pairs were of proteins that were somewhat related to each other , while only three pairs were comprised of two proteins from different kingdoms.
thus the conformational changes are likely to reflect experimentally observable motions rather than evolutionary effects.
resampling of hinge residues
as a further test of confidence in the hinge atlas, we decided to look for sampling artifacts in the hinge set.
resampling <cit>  or bootstrapping <cit>  is a technique suited for this purpose.
we bootstrapped the frequency of occurrence of amino acid types.
the method consists of drawing random samples and computing the frequency of occurrence of a given amino acid type in that sample.
we present the results for glycine, the residue type most overrepresented in hinges.
we randomly chose 1/ <dig> of the  <dig> proteins in the hinge atlas.
the sample was labelled with an index j. within that sample we counted the following:
hj*
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgibasdaqhaawcbagaemoaaogabagamaivcqcaqaaaaaa@316f@ : the number of hinge residues of all amino acid types in sample j,
h˜j*
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacuwgibasgaacamaadaaaleaacqwgqbgaaeaacwag8kokaocaaaaa@317e@ : the number of non-hinge residues of all amino acid types in sample j,
hj*
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgobaadaqhaawcbagaemoaaogabagamaivcqcaqaaaaaa@31af@  : the number of glycine residues in hinges in sample j,
h˜j*
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacuwgobaagaacamaadaaaleaacqwgqbgaaeaacwag8kokaocaaaaa@31be@  : the number of glycines in non-hinge residues in sample j,
fj*=hj*hj*
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgmbgzdaqhaawcbagaemoaaogabagamaipcqcaqaaakiabcicaoiabdggahnaabaaaleaacqwghbwrcqwgmbatcqwgzbqwaeqaaogaeiykakiaeyypa0zaasaaaeaacqwgobaadaqhaawcbagaemoaaogabagamaivcqcaqaaakiabcicaoiabdggahnaabaaaleaacqwghbwrcqwgmbatcqwgzbqwaeqaaogaeiykakcabagaemisag0aa0baasqaaiabdqgaqbqaaiadasqgqagkaaaaaaaa@49b7@ : the sample frequency of occurrence of glycines within hinges within sample j, and
f˜j*=h˜j*h˜j*
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacuwgmbgzgaacamaadaaaleaacqwgqbgaaeaacwagkkokaocaaogaeiikagiaemyyae2aasbaasqaaiabdeeahjabdyeamjabdmfazbqabagccqggpaqkcqgh9aqpdawcaaqaaiqbdigaozaaiawaa0baasqaaiabdqgaqbqaaiadayrgqagkaagccqggoaakcqwghbqydawgaawcbagaem4rackaemitawkaemywakfabeaakiabcmcapaqaaiqbdieaizaaiawaa0baasqaaiabdqgaqbqaaiadasqgqagkaaaaaaaa@49e4@ : the sample frequency of occurrence of glycines among non-hinge residues in sample j.
we repeated the above for for j =  <dig> to  <dig>  randomizing the sample each time.
for the case of agly = glycines, we generated bins  <dig>  wide and counted the number of times values of fj*
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgmbgzdaqhaawcbagaemoaaogabagamaiscqcaqaaaaaa@319f@ and f˜j*
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacuwgmbgzgaacamaadaaaleaacqwgqbgaaeaacwagwkokaocaaaaa@31ae@ occurred in each interval.
the results for glycine are shown in figure  <dig>  if sampling artifacts were present in the dataset, these might manifest themselves as departures from the gaussian distribution.
this was not observed; the distribution is approximately normal.
as a further measure of our confidence in the overrepresentation of this particular amino acid, we can obtain the significance by an alternate  test as follows.
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacuwgmbgzgaacamaadaaaleaacqwgqbgaaeaacwagwkokaocaaaaa@31ae@  and fj*
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgmbgzdaqhaawcbagaemoaaogabagamaiscqcaqaaaaaa@319f@ .
the sample frequency of glycine residues among non-hinge residues in bins containing 1/8th of all hinge atlas proteins was found to average  <dig> .
the sample frequency of glycine among hinge residues in bins containing 1/8th of all hinge atlas proteins was found to average  <dig> .
the standard deviation was considerably larger for the hinge set, since this is a small subset of the hinge atlas.
since the distribution is approximately gaussian the standard deviation of the difference between means should be obtainable by summing the standard deviations of the two sample frequencies in quadrature <cit> .
the z-score of the difference is obtained by dividing the difference between the average sample frequencies by the the thus-obtained standard deviation:
z=f¯gly*−f˜¯gly*2+2= <dig>
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwg6bgecqgh9aqpdawcaaqaaiqbdagamzaarawaa0baasqaaiabdeeahjabdyeamjabdmfazbqaaiadasqgqagkaagccqghsislcuwgmbgzgaacgaqeamaadaaaleaacqwghbwrcqwgmbatcqwgzbqwaeaacwagkkokaocaaagcbawaaoaaaeaacqggoaakiigacqwfdpwcdawgaawcbagaem4rackaemitawkaemywakfabeaakiabcmcapmaacaaaleqabagaegomaidaaogaey4kasiaeiikagiaf83wdmnbagaadawgaawcbagaem4rackaemitawkaemywakfabeaakiabcmcapmaacaaaleqabagaegomaidaaaqabaaaaogaeyypa0jaegymaejaeiola4iaeginaqjaegomaidaaa@54cb@
from the cumulative gaussian distribution <cit> , events  <dig>  or more standard deviations from the mean have a probability of occurrence of  <dig> , giving us an additional measure of confidence that the distribution of glycines is different in the hinge vs. non-hinge sets.
note that we would expect this to be a conservative estimate of the significance .
the main point of this analysis is that the sample is not biased by particular anomalous proteins.
discussion
correlations were found between hinges and several sequence features.
we found that some amino acid types are overrepresented in hinges, and much of this can be explained on the basis of physicochemical properties.
small residues appear to be preferred, especially the "tiny ser, gly, and ala. aliphatic and hydrophobic residues tend not to be in hinges.
we found that residues within four amino acid positions of an active site are significantly more likely to be hinges.
this is most likely related to the fact that hinge bending motion is often related to the catalytic mechanism of the enzyme.
active site residues most logically occur inside the binding cleft and therefore are likely to be in the hinge or close by.
some of these results are intuitive, but are nonetheless useful in buttressing the less expected results.
further, even the intuitive results have in many cases never been rigorously tested or put on a quantitative footing.
surprisingly, hypermutable residues are more likely than conserved residues to occur in hinges.
this was found to be true not only for the hinge atlas set of  <dig> proteins , but also for the subset of  <dig> enzymes with csa annotation .
this may appear to contradict our earlier result that active site residues and their near neighbors are enriched in hinges.
however although the catalytic residue enrichment has very high statistical significance, the number of active site residues in hinges is still small compared to the total number of residues in hinges.
thus their presence is insufficient to counter the wider tendency of hinge residues to be hypermutable.
also, the near neighbors of active site residues have no particular reason to be conserved and thus their enrichment in hinges seems unlikely to counter the tendency toward hypermutability.
this raises the question, why would residues that are functionally important not be conserved?
the answer may be that it is the intricate network of interactions within the hydrophobic core of rigid regions on either side of the hinge that needs to be conserved <cit> , and not the hinges themselves.
the importance of the stability of these domains rather than of any detailed properties of the hinges themselves is underscored by the significant success of structure-based hinge predictors which analyze the interactions within the domains and between the domains and the solvent, but which pay no particular attention to the hinge region itself , or which implicitly <cit>  or explicitly <cit>  find highly interconnected regions of the protein.
one might also ask, is it possible that co-evolution  occurs in hinge residues even in the absence of independent  conservation?
repeatedly investigators have found that co-evolving residue pairs tend to be proximal in space <cit>  and stabilize proteins, for instance by periodically bridging consecutive turns of α-helices or by interacting across the contact interface between two such helices <cit> .
this is an active area of research with possible future implications on hinge finding.
sequence in the immediate neighborhood of a hinge was not found to be sufficient for substantive hinge prediction by a gor-like method, although the latter is successful at predicting secondary structure.
similiarly, no particular sequential pairs of amino acid types were found to be overrepresented in hinges.
however, we did find that combining amino acid propensity data with hinge propensities of active sites and secondary structure yielded some predictive information.
the prediction method we present can easily be extended as additional hinge propensity data is reported.
indeed the publicly available hinge atlas can be used not only to obtain such data but also to test the resulting predictors.
as an additional application, the hinge atlas can potentially be used to help find hinges by homology.
we note, for instance, that a hinge occurring  in the helix connecting the two ef hands of calmodulin has also been found in the evolutionarily related troponin c.
CONCLUSIONS
we found that the amino acids glycine and serine are more likely to occur in hinges, whereas phenylalanine, alanine, valine, and leucine are less likely to occur.
no evidence was found for sequence bias in hinges by a gor-like method, nor for propensity towards sequential pairs of residues.
hinges tend to be small, but not hydrophobic or aliphatic.
they are found less often in α-helices, and more often in turns or random coils.
active site residues were found to coincide significantly with hinges.
interestingly, however, the latter were not conserved.
lastly, hinges are also more likely to occur on the protein surface than in the core.
a consistent picture of hinge residues is suggested.
in this view, hinges often occur near the active site, probably to participate in the bending motion needed for catalysis.
they avoid regions of secondary structure.
they are hypermutable, possibly due to the fact that they occur more often on the surface than in the core.
these correlations yield insights into protein flexibility and the structure-function relationship.
strong sequence-based hinge prediction, however, remains a goal for future work.
authors' contributions
sf annotated hinge locations, performed the statistical studies and wrote the manuscript, web tools, and most of the algorithms.
ll computed the evolutionary conservation of hinge and active site residues.
nc ran flexprot and generated graphics for all morphs in molmovdb, and in other ways provided high performance supercomputing support to the hinge prediction project.
mg supervised the project and edited the paper.
all authors read and approved the final manuscript.

 motions play an essential role in catalysis and protein-ligand interactions. in particular, hinge bending motions account for 45% of motions in a representative set from the database of macromolecular motions comprising domain hinge motions and fragment hinge motions. thus understanding fundamental aspects of hinge bending mechanisms may lead to an improved understanding of the relationship between structure and function. there are three levels of hinge prediction. the easiest case occurs when the atomic coordinates are available for two or more conformations of a given protein. in this case it is possible to visually inspect the motion to determine the hinge location, as we have done here. the process can also be automated with various available packages, including flexprot, hingefind and dyndom. a much more difficult problem is that of predicting hinges when only one set of structural atomic coordinates is available. several algorithms have been developed for this purpose. the very hardest case occurs when the sequence is known but no atomic coordinates are available at all. the problem of finding flexible hinges between rigid regions based on sequence is in some ways similar to the problem of finding domain boundaries, which can be flexible or inflexible. although little work has been done on the former problem, several algorithms exist to address the latter. in one significant contribution, nagarajan and yona analyzed multiple sequence alignments and were able to identify domains with some accuracy. marsden et al focused on the case of proteins with no significant sequence homology to well characterized proteins and found that predicted secondary structure contained information about domain boundaries. jones et al. combined puu, detective, and domak to make a consensus-based domain boundary predictor. heger et al. created the automatic domain decomposition algorithm and associated online database. murzin et al. created the scop database. for purposes of the above algorithms and classifications, however, domains are defined as proteins or regions of proteins having a common evolutionary origin; flexibility is not a consideration. indeed most small and medium sized proteins, such as those prevalent in the hinge atlas, consist of a single domain. therefore the problem of finding flexible hinges is not solved by finding domain boundaries as defined for these methods. schlessinger et al developed a method to predict b-factors from sequence, but it is not clear that b-factors obtained in this way would yield accurate flexibility predictions. in light of the limitations of existing methods, the prediction of domain hinges from sequence is considered an open problem. in this article we focus on the characterization of these hinges based on sequence. to that end, we compiled the hinge atlas, a manually annotated dataset of hinge bending motions, as well as a separate computer annotated dataset, both available for further studies. the hinge atlas has several applications. first, the statistical properties of hinges can be studied. second, it can be used to benchmark hinge prediction programs. third, by homology hinge annotations could potentially be transferred to proteins where the existence and location of a hinge are unknown. fourth, the annotations could conceivably be used in future protein motion prediction programs. the first application was of most interest to us in the current work. our molecular motions database serves a wide variety of purposes, helping investigators understand the motion characteristics of individual proteins, as well as statistical properties of large groups of motions. it is the ideal platform for the current study, since it contains over morphs. a morph is a set of atomic coordinates for two homologous protein structures, plus several structures which our morph server generates as interpolations between the two. our server displays these structures in succession as a "movie" which suggests a possible trajectory of motion between the two conformations. in this study we compiled two representative sets of morphs with hinge annotation: a computer annotated set and a manually annotated set, the hinge atlas. using mostly the latter, we addressed the following questions: are certain residue types differentially represented in hinges? do certain pairs of amino acids coincide with hinges? can sequence be used to predict hinges? do hinges coincide with active sites? do hinges prefer certain secondary structural elements? do hinge residues share physicochemical or steric properties? are hinge residues conserved in evolution? as our first task, we computed the rate of occurrence of each residue type in the hinge atlas. certain amino acids were found to be differentially represented in hinges in a statistically significant fashion. we also investigated whether certain consecutive pairs of residues were differentially represented in hinges. in the course of the above, we observed that one of the overrepresented residues is potentially catalytic; this was the original motivation for question above. to answer that question, we searched the catalytic sites atlas for close homologs to the proteins in our dataset, and extracted the active site residue numbers from those proteins for comparison to the hinge atlas annotation. our next task was to investigate hinge coincidence with secondary structure. hinges are generally believed to occur in disordered regions, but this belief has never been tested or quantified rigorously to our knowledge. following up on our finding that hinges coincide with active site residues, we went on to the question, are hinge residues more likely to be conserved than other residues, as active sites are? we ranked the residues by relative conservation and examined the differences between hinge and non-hinge residues. significant correlations between sequence features and hinges were found in the above analyses. we computed hinge indices for each of these which may be used to relate sequence features to flexibility. we then sought to determine what predictive value sequence might have on its own and whether various sequence features collectively could be used for prediction. we first made a simple gor -like predictor. we computed the log-odds rate of occurrence for residues located at the - to + positions along the sequence in the training set. we used this table to make predictions on the test set and examined their predictive power. as a second approach, we made a composite hinge index, which we call hingeseq, from the hinge indices of each of the sequence features found to be the strongest indicators of flexibility. the statistical significance of this measure was computed much as for the individual sequence features. to show that the measure is predictive, we again divided the hinge atlas into training and test sets and recomputed the relevant hinge indices to include only training set data. we used the regenerated hingeseq to predict hinges in the test set and generated a receiver operating characteristic curve. as a final step, we examined molmovdb as a whole to determine whether any particular database bias was in evidence. we also used resampling to check for sampling artifacts in the hinge atlas. lastly, we compared the hinge atlas to our computer annotated dataset. the resulting work provides insight into the composition, physicochemical properties, geometry, and evolution of hinge regions in proteins. methods preparation of computer annotated hinge dataset prior to generating the manually annotated hinge atlas, we used computational methods to generate a dataset of hinge residues for our statistical studies. we began by running flexprot, a leading hinge identification tool, on all morphs in the database of macromolecular motions flexprot works by matching and structurally aligning fragments in one structure with corresponding fragments in the other. the goal is to find fragment pairs which have minimal rmsd and are maximal in size. the hinges are then reported as the boundaries separating those fragments. goal is equivalent to minimizing the number of these hinges. since domains are never completely rigid, rmsd tends to grow with fragment size and therefore goal is in conflict with goal. this conflict is dealt with by providing the user with a series of adjustable parameters, and further by reporting not one but several alternative hinge locations from which the user can choose. we used a combination of computer and manual culling to select those morphs for which the identified hinges met the following criteria: motion was domain wise, i.e. two or more domains could be observed moving approximately as rigid bodies with respect to each other. the identified hinge was located in the flexible region connecting two rigid domains, rather than in the domains themselves. the morph trajectory was sterically reasonable, i.e. chains were not broken in the attempt to interpolate motion. we found that flexprot's maximal rmsd parameter had a strong effect on the results. therefore when flexprot gave visibly incorrect results for a given morph, we reran the program, systematically varying this parameter. if one of these runs gave sufficiently accurate results, the annotation for that morph was entered into the database. we discarded immediately those morphs that did not exhibit clear hinge bending motion. lastly, we removed redundant morphs using nrdb. note that the definition of a hinge given in the introduction allows for a hinge of zero length. flexprot indeed often returned such hinges. to deal with this, in all cases one residue on each side of the hinge, was taken to also belong to the hinge. thus most hinges are two residues long. at the end of this process, the computer annotated set contained morphs. as described, the computer annotation of hinges requires significant human intervention and the results were often debatable. many of the hinge annotations differed slightly but visibly from the boundary between rigid domains, such that the backbone flexions that could account for the domain motion were not seen in the predicted hinge region. in other cases hinges were missed, and some annotations appeared where no hinge existed. the more flagrantly misannotated hinges were removed from the dataset, but making the manual culling too stringent would simply have resulted in a dataset too small to be statistically meaningful. for these reasons, the computer annotated dataset was not used in most of this work. nonetheless, the computer annotated dataset is arguably more objective then the manually annotated set described below, and so is made available to the community. to address the accuracy issues, we decided to generate a manually annotated set of hinges – the hinge atlas. to generate this set we first created the hinge annotation tool which can also be used by the public as we will now explain. the hinge annotation tool the creation of publicly accessible tools for manual annotation of hinges involved significant changes to the morph page. the morph page is the primary point on molmovdb for analyzing single morphs. it is accessible from the "movies" page or through our search tool, both linked to or visible on our front page. our server also provides a link to this page in an email sent to the submitter of each morph request. we added all of the new tools to the "hinge analysis" tab on this page. the first of these is the hinge annotation tool. each of three rows of "arrow" buttons on this tool move a highlighted window of two residues along the protein chain, allowing the user to highlight up to three hinges in a protein. the "show all" button then highlights all selected residues in the jmol viewer window. once the user is satisfied with the hinge selection, clicking "submit" records this selection in the database. once the morph page is regenerated, a "show public hinge" button will be visible which, when clicked, highlights the selected residues. lastly, the user can use a pointing device to reorient the protein in the jmol window to his/her liking. a gif image based on that view can be generated by clicking on the "color by domain" link. the animation will be rendered using vmd's "new cartoon" style, with the identified hinge region and two rigid domains each colored distinctively. the hinge annotations made in this way persist in our database for visualization and use by others, until overwritten. with minor modification, these tools were used to generate the hinge atlas dataset of manually annotated hinges. the criteria we used for selection are described in the following section. highlighting the hinge atlas hinges on the animated morph movie is a matter of going to the morph page and clicking on the "hinge analysis" tab as above and clicking the "show hinge atlas hinge" button. the annotated hinge location will be rendered in green spacefill style, which contrasts with the white trace used elsewhere in the protein. construction of the hinge atlas the tools described above answer only the technical question of how we annotated hinges. in this section we clarify the motivation for the hinge atlas and its applications and answer the scientific question of how we decided on the precise location of the hinge for each morph. for each morph in the hinge atlas, we used the hinge annotation tool as described to select the hinge location. motivated in part by our long term goal of providing a resource that could be used in motion prediction work, and in part by a desire to deepen basic understanding of protein motion, we asked ourselves the following question: would it be possible to approximately reproduce the observed motion by allowing flexure at the hinge points but keeping the regions between hinges rigid? in order for this question to be answered in the affirmative, the hinge selection should be the one to best meet the following criteria: the φ, ψ, and α torsion angles of hinge residues may often be larger than those of their neighbors. amino acids on either side of the hinge residues must be co-moving with their respective rigid regions. rotations of one of the rigid regions about the hinge region must not result in significant and irreconcilable steric clashes. in order to use as a useful guide to selecting the hinge location, we made use of the torsion angle charts and graphs in the structure analysis tools section on the morph page. however often large rotations of the main chain are induced by multiple cooperative torsions in the hinge, and these may be individually small, particularly in α-helices. the usefulness of this flexibility measure is further limited by the frequent occurrence of large torsion angles which do not coincide with hinges. nonetheless, when the precise location of the hinge was otherwise unclear, torsion angles were often examined to help adjust the selection. criterion is a definition of a hinge. sometimes the hinge was slightly longer than others, and in those cases we added more residues to the hinge, up to a limit of about five residues in total. if the hinge was distributed over too many residues such that no one short stretch could be said to constitute the entire hinge, then the morph was discarded from the hinge atlas, since the motion was not hingelike. criterion is a practical requirement of a working hinge. if substantial flexure at points outside the hinge is required to avoid domain interpenetration, then the choice of hinge location is incorrect, or the motion is not hinge but rather shear or unclassifiable. the next question was, how to select the morphs which would be annotated and included in the hinge atlas. the entire database of macromolecular motions with over morphs, could clearly not all be annotated given limited manpower. further, only a minority of morphs exhibited hinge bending motion, and even within this group much redundancy existed. to address these issues and make the annotation work manageable, we first selected a nonredundant subset of the morphs in molmovdb by aligning all sequences to nrdb. this reduced the dataset to morphs. this was more manageable, but still the set contained many proteins which did not exhibit hinge bending motions. fortunately we found that the score output by flexprot, normalized by dividing by the number of residues, provided an accurate measure of the degree to which a protein exhibited hinge bending. high scores, close to unity, indicated proteins more likely to exhibit hinge bending motion. lower scores, below or so, were very unlikely to do so. we sorted the nonredundant morphs by descending normalized flexprot score and annotated them in that order. those proteins for which we could find hinges allowing a positive answer to the question above were annotated and added to the hinge atlas. those proteins which did not exhibit hinge bending motion or for which no suitable hinge could be found were discarded. at the end of this culling and annotation effort, the hinge atlas contained nonredundant annotated morphs. we also manually annotated a small set of specifically fragment hinge bending motions which may be useful for some studies, described below. availability of datasets in the course of this study we compiled a number of sets of morphs which can be viewed on our online galleries listed and linked to on our sets page. the hinge atlas and computer annotated sets are compared more rigorously in the "statistical comparison of datasets" section. the galleries provide easy browsing and visual inspection of morph movies sharing certain characteristics. the sets offered include: nonredundant no two morphs in this set have more than 90% sequence homology. this set was compiled by alignment to proteins in nrdb catalytic sites atlas all morphs in this set have annotated active sites which can be highlighted in the jmol viewer. catalytic sites atlas same as above, but with redundant morphs removed by comparison to nrdb flexprot hinges computer annotated set used in parts of this study and described above. we consider it to be less useful than the hinge atlas, but the data is nonetheless made available. fragment hinge motions a small set of hinge bending motions involving fragments smaller than domains, as alluded to in the previous section. hinge atlas contains the manually annotated protein pairs used in this study. a link on the sets page permits the download of the sequence data in mysql format. the same data is available in tab-delimited text format which is human readable and importable into ms excel and other packages. another link on the same page facilitates the download of the interpolated structure files associated with each morph in the hinge atlas set. clicking on the thumbnail image leads to the "movies" page, where users can browse through the proteins in the hinge atlas. clicking on any of the protein thumbnail images, in turn, leads to the corresponding morph page, where the hinge annotation can be viewed as described in the "hinge annotation tool" section method for analyzing relative frequency of occurrence throughout this study, we will be comparing how often a particular entity occurs in hinges versus everywhere in the hinge atlas or another of the datasets described above. the statistical analysis will be the same regardless of the particulars, so we will here present the general approach and later only mention adjustments particular to the specific question addressed. first we defined the following variables: d = total number of residues in the dataset h = total number of residues in hinges in the dataset c = classification scheme used to create groups of residue positions. for example, c could be secondary structure, degree of conservation, etc. c = a particular grouping of residues, where c ∈ c. for instance, if c = secondary structure, then c = helix is the class of all residues in helices, c = strand is the class of all residues in strands, etc. another example might be c = evolutionary conservation, with c = cons = top 20% most conserved residues, c = cons = second 20% most conserved, etc. ac = set of all residues of class c in the dataset. dc = number of times residues of class c occurred anywhere in the dataset. hc = number of times residues of a particular class c occurred in hinges. these can be used to estimate various probabilities as follows: p = dc/d is the prior probability of c – in other words, the probability that residues of class c occur anywhere in the dataset. p = hc/h is the conditional probability that a residue belongs to class c, given it is a hinge. a quantity that is of interest in hinge prediction is the posterior probability p, the probability that a residue is a hinge given it is in ac. we obtain this from bayes' rule: p=p⋅pp=hcdc mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgwbaccqggoaakcqwgobaacqgg8bafcqwghbqydawgaawcbagaem4yamgabeaakiabcmcapiabg2da9maalaaabagaemicaanaeiikagiaemyyae2aasbaasqaaiabdogajbqabagccqgg8bafcqwgobaacqggpaqkcqghfly1cqwgwbaccqggoaakcqwgobaacqggpaqkaeaacqwgwbaccqggoaakcqwghbqydawgaawcbagaem4yamgabeaakiabcmcapaaacqgh9aqpdawcaaqaaiabdigaonaabaaaleaacqwgjbwyaeqaaagcbagaemizaq2aasbaasqaaiabdogajbqabaaaaaaa@52b3@ equation 1 where the prior probability that a residue is a hinge is given by p= mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgwbaccqggoaakcqwgobaacqggpaqkcqgh9aqpcqggoaakdawcaaqaaiabdieaibqaaiabdseaebaacqggpaqkaaa@3612@. we further define the hinge index hi, similar to the domain linker index used in armadillo: hi=log⁡10pp mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgibascqwgjbqscqggoaakcqwghbqydawgaawcbagaem4yamgabeaakiabcmcapiabg2da9igbcygasjabc+gavjabcegannaabaaaleaacqaixaqmcqaiwaamaeqaaowaasaaaeaacqwgwbaccqggoaakcqwghbqydawgaawcbagaem4yamgabeaakiabcyha8jabdigaojabcmcapaqaaiabdchawjabcicaoiabdggahnaabaaaleaacqwgjbwyaeqaaogaeiykakcaaaaa@4959@ equation 2 the argument of the log is the ratio of the observed frequency of occurrence of classes of amino acids ac in hinges, over the expected. note that this argument is close to the likelihood ratio pp mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaadawcaaqaaiabdchawjabcicaoiabdggahnaabaaaleaacqwgjbwyaeqaaogaeiifawnaemiaagmaeiykakcabagaemicaanaeiikagiaemyyae2aasbaasqaaiabdogajbqabagccqgg8bafcqgg+bgfcqwgobaacqggpaqkaaaaaa@3fc8@ used in bayesian statistics because h is so small compared to d. the quantity hi yields an intuitive measure of the enrichment of certain classes of residues in hinges, with positive numbers indicating enrichment and negative numbers indicating scarcity. just because the hi is nonzero, however, does not mean that the differential representation has statistical significance. to establish the latter, we considered two statistical hypotheses: h0: the null hypothesis. assume hc is a randomly distributed random variable with mean μh. the null hypothesis states that: μh=dcd⋅h mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaaiigacqwf8oqbdawgaawcbagaemiaaggabeaakiabg2da9maalaaabagaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwgebaraagaeyyxictaemisageaaa@3858@ if this is true, then the hinge set is chosen without replacement in an unbiased fashion from the dataset, and p is given by the hypergeometric distribution. h the alternate hypothesis. this states that: μh≠dcd⋅h. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaaiigacqwf8oqbdawgaawcbagaemiaaggabeaakiabgcmi5oaalaaabagaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwgebaraagaeyyxictaemisagkaeiola4caaa@39fd@ it is equivalent to saying that p is not p, and therefore the null hypothesis can be rejected. we test this as follows. if it is the case that hch>dcd, mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaadawcaaqaaiabdigaonaabaaaleaacqwgjbwyaeqaaagcbagaemisageaaiabg6da+maalaaabagaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwgebaraagaeiilawcaaa@3692@ and if we choose a significance threshold of, we can reject h iff our p-value ∑x=hh∞hyp< mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaadaaewbqaaiabdieaijabdmfazjabdcfaqjabcicaoiabdieaijabcycasiabdseaejabcycasiabdiha4jabcycasiabdsgaknaabaaaleaacqwgjbwyaeqaaogaeiykakiaeyipawjaegimaajaeiola4iaegimaajaegynaudaleaacqwg4baecqgh9aqpcqwgobaacqggoaakcqwghbqycqggpaqkdawgaaadbagaemiaaggabeaaasqaaiabg6hilcqdcqghris5aaaa@4bc9@. the left hand side of the latter inequality can be interpreted as the probability that hc or more residues of class ac could be found in hinges, assuming h and given h, d, and dc. the argument of the sum is the hypergeometric function, which gives the probability that dc residues taken without replacement from a set of d residues of which h are hinges, would contain exactly x hinges: hyp mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgibascqwgzbqwcqwgqbaucqggoaakcqwgibascqggsaalcqwgebarcqggsaalcqwg4baecqggsaalcqwgkbazdawgaawcbagaem4yamgabeaakiabcmcapmaalaaabawaaewaaeaafaqaaegabaaabagaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwg4baeaaaacagloagaayzkaawaaewaaeaafaqaaegabaaabagaemiraqkaeyoei0iaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwgibascqghsislcqwg4baeaaaacagloagaayzkaaaabawaaewaaeaafaqaaegabaaabagaemiraqeabagaemisageaaagaayjkaiaawmcaaaaaaaa@4e8f@ equation 3 otherwise, if it is the case that hch<dcd, mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaadawcaaqaaiabdigaonaabaaaleaacqwgjbwyaeqaaagcbagaemisageaaiabgyda8maalaaabagaemizaq2aasbaasqaaiabdogajbqabaaakeaacqwgebaraagaeiilawcaaa@368e@ then we reject h iff our p-value ∑x=−∞hhyp)< mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaadaaewbqaaiabdieaijabdmfazjabdcfaqjabcicaoiabdieaijabcycasiabdseaejabcycasiabdiha4jabcycasiabdsgakjabcicaoiabdggahnaabaaaleaacqwgpbqaaeqaaogaeiykakiaeiykakiaeyipawjaegimaajaeiola4iaegimaajaegynaudaleaacqwg4baecqgh9aqpcqghsislcqgheispaeaacqwgobaacqggoaakcqwghbqycqggpaqka0gaeyyeiuoaaaa@4e2e@.  are certain amino acids more likely to occur in hinges? we applied the described statistical formalism to the problem of amino acid frequency of occurrence in hinges by taking c = amino acid type, and c to designate each of the canonical amino acids. hi scores and p-values were thus calculated for each of identifications of c corresponding to the canonical amino acids. we found that glycine and serine are overrepresented in a highly significant fashion. we also found phenylalanine, valine, alanine, and leucine to be underrepresented, albeit with lower significance. we also investigated the frequency of occurrence of sequential pairs of amino acids in hinges, but since sequential pairs are possible the significance of the results was much lower and no conclusion could be drawn. are residues within a certain distance of an active site more likely to be hinge residues? as mentioned earlier, the fact that one of the overrepresented residues is potentially catalytic led us to suspect that hinge residues are more likely to occur in active sites, or within a few residues of an active site, than would be expected by chance. this would make sense from a biochemical and mechanical perspective. hinge motions are often opening and closing motions of domains intended to expose the active site, which often would be located at the center of the motion, i.e. the hinge. prior work shows that active sites are more likely to occur at regions of low first normal mode displacement. such regions have been shown to coincide with hinges. here we close the loop, comparing active sites directly with the hinge atlas annotation and quantifying the correspondence. in order to annotate the active site locations, we blasted the morph sequences in the computer annotated dataset against the sequences in the catalytic sites atlas and considered a morph in the hinge dataset to match a protein in the csa if they had sequence identity ≥ 99%. this high threshold was chosen to minimize the possibility of incorrectly labeling a residue in the hinge atlas and thereby diminishing the significance of the results. for each such pair, we transferred the catalytic site annotation to the morph. we described earlier how to browse the csa morphs online. of the proteins in the hinge atlas, were annotated with active site information from the csa; the rest had no close csa homologs. the proteins comprised the dataset for this calculation. we analyzed this set using the statistical formalism described earlier, with the following variable definitions: c = distance from the nearest active site, in residues. c = successively: active site residues, amino acids residue away from the nearest active site residue, residues away, etc. d = residues in the dataset of proteins h = hinge residues in the dataset dc = residues of class c in the dataset hc = residues of class c in hinges. the results are shown in figure and table at short distances from the active sites, hinge residues were overrepresented. the active site residues and residues as much as four residues away from the nearest active site were significantly overrepresented in hinges. are hinges segregated by secondary structure? it is generally accepted that hinges tend to avoid secondary structure. however this belief has, to our knowledge, never been tested on a quantitative basis, and indeed numerous counterexamples can be found. for instance, the hinge in calmodulin and troponin c occurs in an α-helix, and in glutamine binding protein it occurs in two parallel beta strands. thus we do not know which particular types of secondary structure are avoided or preferred, or to what degree. to obtain this information, we tabulated the number of hinge residues occurring in the various types of secondary structural elements, and compared this with the distribution of all residues, proceeding as follows. stride recognizes secondary structural elements from atomic coordinates. we used this program to assign secondary structural classes to all residues in the hinge atlas. we then tabulated the number of residues assigned to each class, both in hinges and elsewhere in the dataset. lastly, we calculated the hi scores and the p-values as before, letting c = secondary structural element type and c designate e.g. helix, coil, etc. we found that three types of secondary structure were differentially represented in hinges with extremely high significance. we conclude that hinges are less likely to occur in α-helices, and are more likely to occur in turns or random coils. for the user's convenience, secondary structure assignments for individual morphs can be obtained in the'hinge analysis tools' box on the molmovdb morph page mentioned earlier. are certain physicochemical properties preferred in hinge residues? it is intuitive that certain physicochemical classes of residues would occur more frequently in hinges, and this would help explain the amino acid propensities reported earlier in this work. to check and quantify this, we grouped amino acids into several non-exlusive categories. following again our statistical treatment, we calculated hi scores and p-values, letting c = physicochemical grouping, and c = aliphatic, polar, charged, etc. we discovered that aliphatic and hydrophobic residues were very significantly underrepresented. overrepresented were small and tiny residues. are hinge residues conserved in evolution? we next investigated whether hinge residues are conserved. since certain residue classes are preferred in hinges, one might suspect that hinge residues would be conserved. first, we blasted each of the hinge atlas sequences against nrdb a non-redundant sequence database in which protein sequences have no more than 90% sequence identity with each other next we extracted up to top-aligned sequences to a given morph to generate a multiple sequence alignment using clustal w. for each position in the multiple sequence alignment, we used the formalism developed by schneider et al to compute the information content associated with a column in the multiple sequence alignment at this position. we sorted the residues in hinge atlas morphs according to the magnitude of the information content scores. we then divided the residues into five bins of equal size. if hinge residues are conserved, then there should be an enrichment of hinge residues in the top bins, which correspond to the most conserved residues. on the other hand, if hinge residues are hypermutable, there should be more of them in the bottom bins, corresponding to the least conserved residues. because it is widely agreed that active sites should be conserved, we used the conservation of active sites as a control. to quantify the enrichment, we calculated the hi scores as described previously. here, c is a label applied to residues that ranked in a given percentile bin, e.g. the top 20% most conserved. for that bin p = hc/h is thus the ratio of the number of hinge residues in the bin divided by the total number of hinge residues. similarly, p = dc/d is the ratio of the number of residues in the dataset in the bin divided by the grand total of residues in the dataset. to determine the statistical significance of hi scores, we calculated the p-values using the hypergeometric distribution with the dc, hc, d, h defined above. for the control set, we performed the same calculation but made the following changes to the variable definitions: our dataset was no longer the hinge atlas, but rather the "catalytic sites atlas " set described earlier. d is the total number of residues in this set. ac still represents residues in the dataset belonging to a given conservation rank bin. dc is the total number of residues in that bin. hc now represents the number of active site residues in a given bin corresponding to c. similarly, h represents the total number of active site residues in the dataset. we found that hinge residues distribute evenly in the top 80%, and have a slight but statistically very significant enrichment in the bottom 20% bin. thus hinge residues are hypermutable. we observe a highly significant enrichment of active site residues in the top 20% bin, as expected. among the active site residues, of them are in the highest bin, and the numbers progressively decrease in lower bins. the hinge atlas pools enzymes together with non-catalytic proteins. we reasoned therefore that perhaps only hinges in non-catalytic proteins are hypermutable, and that if we analyzed a set consisting only of enzymes, then the propensity of active sites to occur in hinges would lead to conservation, rather than hypermutability of hinge residues for that set. to test this idea, we decided to calculate the propensity of hinges to occur in specific bins of conservation score, for the proteins in the hinge atlas with csa annotation, rather than for the larger set of for this set we also found that hinge residues occur more frequently among the 20% least conserved residues for each protein. at a p-value of, the confidence in this result is high. even this test, however pools together hinges that are near the active site with hinges that occur at some distance from it. so we selected from the proteins a small set that had at least one active site residue in the hinge, and removed the active site residues themselves. we then calculated the propensity of hinge residues to occur in the five conservation bins. this set was found to be too small, however, and statistical significance was too low to draw a conclusion. a study using the set of fragment hinge motions described earlier was similarly inconclusive. the hypermutability of hinge residues that we found is reasonable because hinge residues tend to be on the surface of proteins rather than in the more highly conserved core. hinges are less likely to be buried inside domains because they would then be highly coordinated with near neighbors and hence less flexible. the apparent contradiction of hypermutability on the one hand and enrichment of active sites on the other is dealt with in the discussion section. are hinge residues more likely to occur on the surface? to support our argument that hinge residues are hypermutable partly because they occur on the surface, we quantified the degree to which the latter is the case. to do this, we used a solvent accessible surface area calculation program with a probe radius of Å. the asa of each of the backbone heavy atoms was calculated and summed for each residue in each protein in the hinge atlas. we then binned the residues by this quantity. lastly, we counted the number of hinge residues in each bin and calculated hi and p-value as before. as expected, bin # was significantly enriched with hinges. bin # was also highly enriched, while bins # and # had fewer hinges, all with extremely high significance. bins 1– contain total residues and bin contains how many hinge sites appear in each protein? perhaps the simplest hinge consists of a single point on the chain separating two rigid regions. however it is also possible for the chain to pass multiple times through the same region, or to have multiple independent hinge regions. this leads to the question, how many proteins had single hinge points, versus a larger number of hinge points? we answer this question in table most morphs had three or fewer hinge points. can hinges be predicted by a simple gor-like method? gor method is useful for predicting secondary structure from sequence with fair accuracy. we implemented a gor-like method to determine whether sequence contained enough information for hinge prediction. we divided the dataset into a training set and a test set for this study. the log-odds frequency of occurrence of amino acids in the training set were tabulated not only at a given hinge residue, but also at positions ranging from - to + from the given residue in sequence space. for simplicity, hinge residues at positions less than eight residues from either end of the chain were not included. once the table was generated, it was used on the test set. the score for a given residue was taken to be the sum of the scores for the residues in positions - to + from that residue. the scores were computed for all residues in the test set, except those less than eight residues from either end of the chain. the idea is that a threshold score can be chosen and residues scoring higher than this threshold are considered more likely to be hinges. note that where robson and suzuki used a different fitting parameter for each type of secondary structure, we used no fitting parameter, since we were interested in only one "secondary structure": the hinges. the rates of true and false positives and negatives were calculated for each choice of score threshold over a range. our training set numbered proteins from the computer annotated set. we tested the method on a test set of proteins from the same set and obtained a roc curve. the area under this curve was nearly, indicating negligible predictive value. hinge prediction by combining sequence features as the gor-like method did not work well, we sought to measure the predictive power of the various sequence features studied above. the hi scores we have reported provide an intuitive means of weighing the relative predictive value of each sequence feature. we show how to combine the hi scores for several features in order to make a more powerful predictor, which we call hingeseq. we define this predictor as follows: hs=log⁡10ppppp)=hiamino⋅acid+hisecondary⋅structure+hiactive⋅site mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgibascqwgtbwucqggoaakcqwgpbqacqggpaqkcqgh9aqpcyggsbabcqggvbwbcqggnbwzdawgaawcbagaegymaejaegimaadabeaakmaabmaabawaasaaaeaacqwgwbaccqggoaakcqwghbqydawgaawcbagaemoaaogabeaakiabcyha8jabdigaojabcmcapiabdchawjabcicaoiabdggahnaabaaaleaacqwgrbwaaeqaaogaeiifawnaemiaagmaeiykakiaemicaanaeiikagiaemyyae2aasbaasqaaiabdygasbqabagccqgg8bafcqwgobaacqggpaqkaeaacqwgwbaccqggoaakcqwghbqydawgaawcbagaemoaaogabeaakiabcmcapiabdchawjabcicaoiabdggahnaabaaaleaacqwgrbwaaeqaaogaeiykakiaemicaanaeiikagiaemyyae2aasbaasqaaiabdygasbqabagccqggpaqkaaaacagloagaayzkaagaeyypa0jaemisagkaemysak0aasbaasqaagqaciab=fgahjab=1gatjab=lgapjab=5gaujab=9gavjabgwsixlabdggahjabdogajjabdmgapjabdsgakbqabagccqggoaakcqwgpbqacqggpaqkcqghrawkcqwgibascqwgjbqsdawgaawcbagae83camnae8xzaumae83yammae83ba8mae8nba4mae8hzaqmae8xyaemae8ncainae8xeaknaeyyxictaem4camnaemidaqnaemocainaemydaunaem4yammaemidaqnaemydaunaemocainaemyzaugabeaakiabcicaoiabdmgapjabcmcapiabgucariabdieaijabdmeajnaabaaaleaacqwghbqycqwgjbwycqwg0badcqwgpbqacqwg2bgdcqwglbqzcqghfly1cqwgzbwccqwgpbqacqwg0badcqwglbqzaeqaaogaeiikagiaemyaakmaeiykakcaaa@b40a@ equation 4 for simplicity, statistical independence of the various features was assumed in creating this definition. here the i's correspond to individual amino acids in the protein sequence. for each i, j designates one of the amino acid types, k designates the secondary structural classification, and l designates active site versus non-active site classification. thus hiamino·acid is assigned according to residue type by looking up the corresponding value in table similarly, hisecondary·structure isobtained according to secondary structure type from table following table approximately, we assign hiactive·site as for residues four or fewer amino acid positions away from the nearest active site residue, and elsewhere. the highest values of hs correspond to residues most likely to occur in hinges. clearly, extending this method is only a matter of obtaining amino acid propensities to occur in hinges according to additional classifications. the resulting index can then simply be included as an additional term in the above formula, with no need for adjustable weighting factors. we evaluated the statistical significance of this measure much as for the individual sequence features. we counted the number of residues in the hinge atlas with a hingeseq score above, and within that set the number of hinge residues. we compared this to the total number of hinges and the population size of the hinge atlas. using the cumulative hypergeometric distribution as before, we computed a p-value of order 10- thus the measure shows high statistical significance. however since only about 5% of the residues scoring over were annotated hinges, hingeseq is not likely to be sensitive enough to be used alone for hinge prediction. the low p-value indicates that the predictor results have high statistical significance. however the low sensitivity limits its potential predictive value. we nonetheless wished to show that hingeseq is predictive, rather simply reflectling peculiarities of the dataset. to this end, we divided the proteins of the hinge atlas into a training set numbering proteins, and a test set numbering of the hinge atlas proteins, the proteins with annotation from the csa were apportioned such that were included in the training set and in the test set. we tested the performance of the predictor by means of roc curves. we need to define a few terms in order to use these: test positives: residues with hsgreater than or equal to a certain threshold. test negatives: residues with hsless than a certain threshold. gold standard positives: residues annotated as hinges in the hinge atlas. gold standard negatives: residues which are not in hinges according to the hinge atlas annotation. true positives: those residues that are both test positives and gold standard positives. true negatives: residues that are both test negatives and gold standard negatives. false positives: residues that are test positives and gold standard negatives. false negatives: residues that are test negatives and gold standard positives. sensitivity=tptp+fn mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgzbwccqwglbqzcqwgubgbcqwgzbwccqwgpbqacqwg0badcqwgpbqacqwg2bgdcqwgpbqacqwg0badcqwg5bqecqgh9aqpdawcaaqaaiabdsfaujabdcfaqbqaaiabdsfaujabdcfaqjabgucariabdaeagjabd6eaobaaaaa@450b@ specificity=tnfp+tn mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgzbwccqwgwbaccqwglbqzcqwgjbwycqwgpbqacqwgmbgzcqwgpbqacqwgjbwycqwgpbqacqwg0badcqwg5bqecqgh9aqpdawcaaqaaiabdsfaujabd6eaobqaaiabdaeagjabdcfaqjabgucariabdsfaujabd6eaobaaaaa@44a9@ 1−specificity=fpfp+tn mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqaixaqmcqghsislcqwgzbwccqwgwbaccqwglbqzcqwgjbwycqwgpbqacqwgmbgzcqwgpbqacqwgjbwycqwgpbqacqwg0badcqwg5bqecqgh9aqpdawcaaqaaiabdaeagjabdcfaqbqaaiabdaeagjabdcfaqjabgucariabdsfaujabd6eaobaaaaa@466e@ the roc curve is simply a plot of the true positive rate vs. false positive rate, for each value of the threshold, as the threshold is varied from + to - a range which included all possible values of hs. for a good predictor, the true positive rate will increase faster than the false positive rate as the threshold is lowered, and the area under the curve will be significantly greater than. the roc curve is shown in figure although work remains to be done before sequence-based hinge prediction can be relied upon exclusively, hingeseq displays significant ability to detect potential for flexibility. checking for dataset bias these findings assume that the dataset used does not contain significant bias or artifacts, either in the composition of the entire dataset or of the hinges within it. to substantiate this, we performed various studies as follows. bias in amino acid composition and functional classification in order to find out whether the molmovdb database contained any bias in amino acid composition, we extracted the sequences of all the morphs in molmovdb and counted the total occurrence of each residue type. suspecting that redundancies might bias the result, we clustered the sequences and recounted the amino acid residues in the same way. we compared these numbers to publicly available amino acid frequencies of occurrence for the pdb. the amino acid frequency of occurrence for the clustered molmovdb morphs was found to be essentially that of the pdb, from which it was created, therefore no particular database bias is in evidence. we also sought to determine whether there existed a bias towards particular protein classes, in either the hinge atlas or the nonredundant set of molmovdb morphs from which it was compiled. to do this, we first counted the number of times each top-level gene ontology term under the "molecular function" ontology was associated with a protein in the hinge atlas. where the annotation was given for deeper levels, we traced up the hierarchical tree to retrieve the corresponding top level term in the ontology. thus we found, for example, that proteins in the hinge atlas were associated with the term "nucleic acid binding." we repeated this procedure for the pdb as a whole as well as for the non-redundant set of morphs in molmovdb from which the hinge atlas was compiled. the results for the most frequently encountered go terms are shown in table to compare the hinge atlas counts to the pdb counts in an overall fashion, we used the chi-square distribution with degrees of freedom and obtained a chi-square value of. this corresponds to a p-value of, so there is no statistically significant difference in the distribution of these terms in the hinge atlas vs. the entire protein data bank. statistical comparison of datasets the hinge atlas and computer annotated sets were compiled differently, therefore one might suspect that the hinges from one set might comprise a statistically different population from the hinges of the other set. if this were the case, then one of the two sets would be preferable to the other, otherwise if the populations were essentially the same then the two sets could potentially be used interchangeably. it is therefore necessary to quantitatively compare these two populations. it is also necessary to confirm that within one set, the hinge residues are a statistically distinct population from the rest of the set; if this were not true then the amino acid propensity data reported earlier would not be meaningful. although the hinge atlas and the computer annotated set share a total of residues, only of these are hinge residues. this is another reason to suspect that the hinge population of the hinge atlas is statistically different from the hinge population of the computer annotated set. to test this, we computed the chi-square value for hinge atlas hinges vs. computer annotated hinges, and obtained a p-value of. therefore, the hinge atlas hinges are different from the computer annotated hinges. the chi-square value describing the difference between amino acid frequency of occurrence in the hinge vs. non-hinge subsets of the hinge atlas was. with degrees of freedom this corresponds to a p-value below 10-. therefore, the hinge residues are shown with high confidence to be different from non-hinge residues in the hinge atlas. a similar calculation yielded a p-value of for the computer annotated set. therefore we conclude that the hinge and non-hinge populations are different for the computer annotated set, as well. the same is true for the hinge atlas set. hinges in the hinge atlas are likely to be a distinct population from hinges in the computer annotated set; therefore for future studies one or the other should be used, and we recommend that it be the hinge atlas. we conclude from this calculation for both the hinge atlas and computer annotated set, the hinge population is different from the non-hinge population, therefore statistically significant information can be extracted from both. however the hinge population of the hinge atlas is different from that of the computer annotated set, albeit with much lower significance. we argue that one of the two sets should therefore be preferred for statistical studies. the preferred set should be the hinge atlas since the computer annotated set contains numerous annotations which are slightly different from the correct and visually verifiable hinge location. similarity within morph pairs we next asked the question, do the morphs in the hinge atlas reflect intrinsic flexibility of the protein, or is the apparent conformational change the result of sequence differences between the two structures in the pair? that is, do the morphs display motions observable in a single protein, or do they instead represent evolutionary change? to answer this we counted the number of times both structures in the morph came from the same vs different organisms. of the morphs, had structures downloaded directly from the pdb rather than uploaded by users, and also had valid source organism data. for of the both proteins in the pair came from the same species, while for the two proteins came from different species. of the pairs were of proteins that were somewhat related to each other, while only three pairs were comprised of two proteins from different kingdoms. thus the conformational changes are likely to reflect experimentally observable motions rather than evolutionary effects. resampling of hinge residues as a further test of confidence in the hinge atlas, we decided to look for sampling artifacts in the hinge set. resampling or bootstrapping is a technique suited for this purpose. we bootstrapped the frequency of occurrence of amino acid types. the method consists of drawing random samples and computing the frequency of occurrence of a given amino acid type in that sample. we present the results for glycine, the residue type most overrepresented in hinges. we randomly chose 1/ of the proteins in the hinge atlas. the sample was labelled with an index j. within that sample we counted the following: hj* mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgibasdaqhaawcbagaemoaaogabagamaivcqcaqaaaaaa@316f@: the number of hinge residues of all amino acid types in sample j, h˜j* mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacuwgibasgaacamaadaaaleaacqwgqbgaaeaacwag8kokaocaaaaa@317e@: the number of non-hinge residues of all amino acid types in sample j, hj* mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgobaadaqhaawcbagaemoaaogabagamaivcqcaqaaaaaa@31af@: the number of glycine residues in hinges in sample j, h˜j* mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacuwgobaagaacamaadaaaleaacqwgqbgaaeaacwag8kokaocaaaaa@31be@: the number of glycines in non-hinge residues in sample j, fj*=hj*hj* mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgmbgzdaqhaawcbagaemoaaogabagamaipcqcaqaaakiabcicaoiabdggahnaabaaaleaacqwghbwrcqwgmbatcqwgzbqwaeqaaogaeiykakiaeyypa0zaasaaaeaacqwgobaadaqhaawcbagaemoaaogabagamaivcqcaqaaakiabcicaoiabdggahnaabaaaleaacqwghbwrcqwgmbatcqwgzbqwaeqaaogaeiykakcabagaemisag0aa0baasqaaiabdqgaqbqaaiadasqgqagkaaaaaaaa@49b7@: the sample frequency of occurrence of glycines within hinges within sample j, and f˜j*=h˜j*h˜j* mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacuwgmbgzgaacamaadaaaleaacqwgqbgaaeaacwagkkokaocaaogaeiikagiaemyyae2aasbaasqaaiabdeeahjabdyeamjabdmfazbqabagccqggpaqkcqgh9aqpdawcaaqaaiqbdigaozaaiawaa0baasqaaiabdqgaqbqaaiadayrgqagkaagccqggoaakcqwghbqydawgaawcbagaem4rackaemitawkaemywakfabeaakiabcmcapaqaaiqbdieaizaaiawaa0baasqaaiabdqgaqbqaaiadasqgqagkaaaaaaaa@49e4@: the sample frequency of occurrence of glycines among non-hinge residues in sample j. we repeated the above for for j = to randomizing the sample each time. for the case of agly = glycines, we generated bins wide and counted the number of times values of fj* mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgmbgzdaqhaawcbagaemoaaogabagamaiscqcaqaaaaaa@319f@ and f˜j* mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacuwgmbgzgaacamaadaaaleaacqwgqbgaaeaacwagwkokaocaaaaa@31ae@ occurred in each interval. the results for glycine are shown in figure if sampling artifacts were present in the dataset, these might manifest themselves as departures from the gaussian distribution. this was not observed; the distribution is approximately normal. as a further measure of our confidence in the overrepresentation of this particular amino acid, we can obtain the significance by an alternate test as follows. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacuwgmbgzgaacamaadaaaleaacqwgqbgaaeaacwagwkokaocaaaaa@31ae@ and fj* mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwgmbgzdaqhaawcbagaemoaaogabagamaiscqcaqaaaaaa@319f@. the sample frequency of glycine residues among non-hinge residues in bins containing 1/8th of all hinge atlas proteins was found to average. the sample frequency of glycine among hinge residues in bins containing 1/8th of all hinge atlas proteins was found to average. the standard deviation was considerably larger for the hinge set, since this is a small subset of the hinge atlas. since the distribution is approximately gaussian the standard deviation of the difference between means should be obtainable by summing the standard deviations of the two sample frequencies in quadrature. the z-score of the difference is obtained by dividing the difference between the average sample frequencies by the the thus-obtained standard deviation: z=f¯gly*−f˜¯gly*2+2= mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaach8aky=wiffydh8gipec8eeeu0xxdbba9frfj0=oqffea0dxdd9vqai=hguq8kuc9pgc9s8qqaq=dirpe0xb9q8qilsfr0=vr0=vr0dc8meaabaqaciaacagaaeqabaqabegadaaakeaacqwg6bgecqgh9aqpdawcaaqaaiqbdagamzaarawaa0baasqaaiabdeeahjabdyeamjabdmfazbqaaiadasqgqagkaagccqghsislcuwgmbgzgaacgaqeamaadaaaleaacqwghbwrcqwgmbatcqwgzbqwaeaacwagkkokaocaaagcbawaaoaaaeaacqggoaakiigacqwfdpwcdawgaawcbagaem4rackaemitawkaemywakfabeaakiabcmcapmaacaaaleqabagaegomaidaaogaey4kasiaeiikagiaf83wdmnbagaadawgaawcbagaem4rackaemitawkaemywakfabeaakiabcmcapmaacaaaleqabagaegomaidaaaqabaaaaogaeyypa0jaegymaejaeiola4iaeginaqjaegomaidaaa@54cb@ from the cumulative gaussian distribution, events or more standard deviations from the mean have a probability of occurrence of, giving us an additional measure of confidence that the distribution of glycines is different in the hinge vs. non-hinge sets. note that we would expect this to be a conservative estimate of the significance. the main point of this analysis is that the sample is not biased by particular anomalous proteins. discussion correlations were found between hinges and several sequence features. we found that some amino acid types are overrepresented in hinges, and much of this can be explained on the basis of physicochemical properties. small residues appear to be preferred, especially the "tiny ser, gly, and ala. aliphatic and hydrophobic residues tend not to be in hinges. we found that residues within four amino acid positions of an active site are significantly more likely to be hinges. this is most likely related to the fact that hinge bending motion is often related to the catalytic mechanism of the enzyme. active site residues most logically occur inside the binding cleft and therefore are likely to be in the hinge or close by. some of these results are intuitive, but are nonetheless useful in buttressing the less expected results. further, even the intuitive results have in many cases never been rigorously tested or put on a quantitative footing. surprisingly, hypermutable residues are more likely than conserved residues to occur in hinges. this was found to be true not only for the hinge atlas set of proteins, but also for the subset of enzymes with csa annotation. this may appear to contradict our earlier result that active site residues and their near neighbors are enriched in hinges. however although the catalytic residue enrichment has very high statistical significance, the number of active site residues in hinges is still small compared to the total number of residues in hinges. thus their presence is insufficient to counter the wider tendency of hinge residues to be hypermutable. also, the near neighbors of active site residues have no particular reason to be conserved and thus their enrichment in hinges seems unlikely to counter the tendency toward hypermutability. this raises the question, why would residues that are functionally important not be conserved? the answer may be that it is the intricate network of interactions within the hydrophobic core of rigid regions on either side of the hinge that needs to be conserved, and not the hinges themselves. the importance of the stability of these domains rather than of any detailed properties of the hinges themselves is underscored by the significant success of structure-based hinge predictors which analyze the interactions within the domains and between the domains and the solvent, but which pay no particular attention to the hinge region itself, or which implicitly or explicitly find highly interconnected regions of the protein. one might also ask, is it possible that co-evolution occurs in hinge residues even in the absence of independent conservation? repeatedly investigators have found that co-evolving residue pairs tend to be proximal in space and stabilize proteins, for instance by periodically bridging consecutive turns of α-helices or by interacting across the contact interface between two such helices. this is an active area of research with possible future implications on hinge finding. sequence in the immediate neighborhood of a hinge was not found to be sufficient for substantive hinge prediction by a gor-like method, although the latter is successful at predicting secondary structure. similiarly, no particular sequential pairs of amino acid types were found to be overrepresented in hinges. however, we did find that combining amino acid propensity data with hinge propensities of active sites and secondary structure yielded some predictive information. the prediction method we present can easily be extended as additional hinge propensity data is reported. indeed the publicly available hinge atlas can be used not only to obtain such data but also to test the resulting predictors. as an additional application, the hinge atlas can potentially be used to help find hinges by homology. we note, for instance, that a hinge occurring in the helix connecting the two ef hands of calmodulin has also been found in the evolutionarily related troponin c.  we found that the amino acids glycine and serine are more likely to occur in hinges, whereas phenylalanine, alanine, valine, and leucine are less likely to occur. no evidence was found for sequence bias in hinges by a gor-like method, nor for propensity towards sequential pairs of residues. hinges tend to be small, but not hydrophobic or aliphatic. they are found less often in α-helices, and more often in turns or random coils. active site residues were found to coincide significantly with hinges. interestingly, however, the latter were not conserved. lastly, hinges are also more likely to occur on the protein surface than in the core. a consistent picture of hinge residues is suggested. in this view, hinges often occur near the active site, probably to participate in the bending motion needed for catalysis. they avoid regions of secondary structure. they are hypermutable, possibly due to the fact that they occur more often on the surface than in the core. these correlations yield insights into protein flexibility and the structure-function relationship. strong sequence-based hinge prediction, however, remains a goal for future work. authors' contributions sf annotated hinge locations, performed the statistical studies and wrote the manuscript, web tools, and most of the algorithms. ll computed the evolutionary conservation of hinge and active site residues. nc ran flexprot and generated graphics for all morphs in molmovdb, and in other ways provided high performance supercomputing support to the hinge prediction project. mg supervised the project and edited the paper. all authors read and approved the final manuscript. 
therefore it is important to relate the hinge location to sequence features such as residue type, physicochemical class, secondary structure, solvent exposure, evolutionary conservation, and proximity to active sites. to do this, we first generated the hinge atlas, a set of protein motions with the hinge locations manually annotated, and then studied the coincidence of these features with the hinge location. motions play an important role in catalysis and protein-ligand interactions. we evaluate the potential for hinge prediction based on sequence. using the molecular motions database we have created a hinge atlas of manually annotated hinges and a statistical formalism for calculating the enrichment of various types of residues in these hinges. less surprisingly, we found that hinge residues tend to be small, not hydrophobic or aliphatic, and occur in turns and random coils on the surface. in particular, we found that hinges tend to coincide with active sites, but unlike the latter they are not at all conserved in evolution. most interestingly, we found that hinges tend to occur at or near active sites and yet unlike the latter are not conserved. the hinge atlas is made available to the community for further flexibility studies. we found that all of the features have bearing on the hinge location. relating features of protein sequences to structural hinges is important for identifying domain boundaries, understanding structure-function relationships, and designing flexibility into proteins.
 relating features of protein sequences to structural hinges is important for identifying domain boundaries, understanding structure-function relationships, and designing flexibility into proteins. efforts in this field have been hampered by the lack of a proper dataset for studying characteristics of hinges.  using the molecular motions database we have created a hinge atlas of manually annotated hinges and a statistical formalism for calculating the enrichment of various types of residues in these hinges.  we found various correlations between hinges and sequence features. some of these are expected; for instance, we found that hinges tend to occur on the surface and in coils and turns and to be enriched with small and hydrophilic residues. others are less obvious and intuitive. in particular, we found that hinges tend to coincide with active sites, but unlike the latter they are not at all conserved in evolution. we evaluate the potential for hinge prediction based on sequence. motions play an important role in catalysis and protein-ligand interactions. hinge bending motions comprise the largest class of known motions. therefore it is important to relate the hinge location to sequence features such as residue type, physicochemical class, secondary structure, solvent exposure, evolutionary conservation, and proximity to active sites. to do this, we first generated the hinge atlas, a set of protein motions with the hinge locations manually annotated, and then studied the coincidence of these features with the hinge location. we found that all of the features have bearing on the hinge location. most interestingly, we found that hinges tend to occur at or near active sites and yet unlike the latter are not conserved. less surprisingly, we found that hinge residues tend to be small, not hydrophobic or aliphatic, and occur in turns and random coils on the surface. a functional sequence based hinge predictor was made which uses some of the data generated in this study. the hinge atlas is made available to the community for further flexibility studies. 