How to connect proteins via disulfide bonds computationally?

I want to connect proteins together to form a dimer. As seen in the picture, the monomers come close to each other along the edge of the higher-order structure (forget about the sulfate there). How can I find out, where I need to introduce an Xaa->Cys mutation so that the most stable disulfide bond is formed between the proteins? It is easy to get a visual clue where about a mutation might be good, but is there a good computational method to calculate the energies and scan through some possibilities?


(I would also go with hydrophobic or van de Waals interactions. Basically, I just want to glue those proteins together.)


[![Edge of two proteins (modelled in Chimera)](https://i.stack.imgur.com/JfVAD.png)](https://i.stack.imgur.com/JfVAD.png)




If you have a structural model for the interaction, i.e. atomic coordinates, then you could use our software called Disulfide by Design. See the following links.


<http://cptweb.cpt.wayne.edu/DbD2/>


<http://www.biomedcentral.com/1471-2105/14/346>


 1