Regulatory networks consist of interacting molecules with a high degree of
mutual chemical specificity. How can these molecules evolve when their function
depends on maintenance of interactions with cognate partners and simultaneous
avoidance of deleterious "crosstalk" with non-cognate molecules? Although
physical models of molecular interactions provide a framework in which
co-evolution of network components can be analyzed, most theoretical studies
have focused on the evolution of individual alleles, neglecting the network. In
contrast, we study the elementary step in the evolution of gene regulatory
networks: duplication of a transcription factor followed by selection for TFs
to specialize their inputs as well as the regulation of their downstream genes.
We show how to coarse grain the complete, biophysically realistic
genotype-phenotype map for this process into macroscopic functional outcomes
and quantify the probability of attaining each. We determine which evolutionary
and biophysical parameters bias evolutionary trajectories towards fast
emergence of new functions and show that this can be greatly facilitated by the
availability of "promiscuity-promoting" mutations that affect TF specificity