Cooperative interactions link the behavior of different amino acid residues within a protein molecule. As a result, the effects of chemical or physical perturbations to any given residue are propagated to other residues by an intricate network of interactions. Very often, amino acids "sense" the effects of perturbations occurring at very distant locations in the protein molecule. In these studies, we have investigated by computer simulation the structural distribution of those interactions. We show here that cooperative interactions are not intrinsically bi-directional and that different residues play different roles within the intricate network of interactions existing in a protein. The effect of a perturbation to residue j on residue k is not necessarily equal to the effect of the same perturbation to residue k on residue j. In this paper, we introduce a computer algorithm aimed at mapping the network of cooperative interactions within a protein. This algorithm exhaustively performs single site thermodynamic mutations to each residue in the protein and examines the effects of those mutations on the distribution of conformational states. The algorithm has been applied to three different proteins (lambda repressor fragment 6-85, chymotrypsin inhibitor 2, and barnase). This algorithm accounts well for the observed behavior of these proteins.