Likelihood estimates of local rates of evolution within proteins reveal that selective constraints on structure and function are quantitatively stable over billions of years of divergence. The stability of constraints produces an intramolecular clock that gives each protein a characteristic pattern of evolutionary rates along its sequence. This pattern allows the identification of constrained regions and, because the rate of evolution is a quantitative measure of the strength of the constraint, of their functional importance. We show that results from such analyses, which require only sequence alignments, are consistent with experimental and mutational data. The methodology has significant predictive power and may be used to guide structure-function studies for any protein represented by a modest number of homologs in sequence databases.T he principle that the rate of molecular evolution is inversely correlated with the strength of selective constraints has long been known (1, 2). The average evolutionary rate of a protein reflects the overall importance of the protein for organismal functions, whereas rate variation within the protein reflects intramolecular differences in structural and functional constraints. Intramolecular rate variation has been the subject of many studies focused on devising more realistic models of sequence evolution that do not assume rate constancy among sites (e.g., refs. 3-6). A more recent application of estimating rate variation within proteins has been the inference of structural and functional constraints (7-9).To identify evolutionarily constrained regions (ECRs) we devised a general approach to inferring rate variation within proteins. We construct a multiple sequence alignment of orthologs and͞or closely related paralogs and build the maximum likelihood tree. Holding the branching structure of the tree fixed, we then calculate the number of substitutions in each window of a fixed width over the entire alignment. The ''relative rate'' in the window is obtained by dividing the number of substitutions per site in the window by the average of all windows. Plotting the windows' relative rates as a function of their position generates a rate profile (RP), and a heuristic algorithm automatically identifies ECRs and ranks them by their rate of evolution. This approach allows us to infer both the existence of constrained regions in a protein and, because the rate of evolution is a quantitative measure of the strength of the constraint, the relative importance of the identified region.Our method requires only a multiple sequence alignment and is sufficiently powerful to allow analyses involving a relatively small number of fairly closely related sequences. It enables us (i) to use sequences for which the quality of the alignment over most its length is indisputably robust, and (ii) to use alignments of orthologs and closely related paralogs for which conservation of structural and functional constraints can be reasonably assumed. We show below that the method identifies known domains wit...