Clustered protocadherins are a large family of paralogous proteins that play important roles in neuronal development. The more than 50 clustered protocadherin isoforms have remarkable homophilic specificity for interactions between cellular surfaces that is controlled by a large antiparallel dimer interface formed by the first four extracellular cadherin (EC) domains. To understand how specificity is achieved between the numerous paralogs, we used a combination of structural and computational approaches. Molecular dynamics simulations revealed that individual EC interactions are weak and go through binding and unbinding events, but together they form a stable complex through polyvalency. Using sequence coevolution, we generated a statistical model of interaction energy for the clustered protocadherin family that measures the contributions of all amino acid pairs in the interface. Our interaction energy model assesses specificity for all possible pairs of isoforms, recapitulating known pairings and predicting the effects of experimental changes in isoform specificity that are consistent with literature results. Our results show that sequence coevolution can be used to understand specificity determinants in a protein family and prioritize interface amino acid substitutions to reprogram specific protein-protein interactions. Clustered protocadherins (Pcdhs) are a large protein family (53 isoforms in humans) that play roles in vertebrate nervous system development, including neuronal survival, axon targeting, neuronal arborization, and dendritic self-avoidance (1-9). Dendritic self-avoidance is mediated by formation of a clustered Pcdh assembly between two dendrites (10). This assembly relies on individual recognition units formed in trans across two cellular membranes that consist of homodimers of the first four extracellular cadherin-repeat (EC) domains in an antiparallel arrangement ( Figure 1A). These homodimers are highly specific such that no cross interactions are observed in even the most similar isoforms (11-13). The trans EC1-4 interaction is also found in other non-clustered Pcdhs (14-16), which have roles in nervous system development and maintenance (17), indicating the importance of this recognition unit in cognitive function. Given that there are many isoforms per vertebrate genome, we sought to understand how specificity is achieved in this large interface.Structures of these recognition units (14,(18)(19)(20) have revealed idiosyncratic characteristics of individual dimer structures, such as the lack of the EC1/EC4 interaction in the structure of PcdhγA1 and PcdhγA8 (20) and the small interface of the EC2/EC3 interaction in PcdhγB3 (14). More subtle structural differences can be observed broadly between the three clustered Pcdh subfamilies, α, β, and γ (14,19,20). Based on the variety of interfaces found in the existing crystal structures (14,(18)(19)(20), it is possible that every isoform achieves specificity by adopting a different static interface conformation, or that isoforms sample a distributi...