System for molecular packing calculations

A method of engineering a variant protein relative to a parent protein, the parent protein comprising a first residue comprising one or more residue atoms, the method comprising:

• (a) calculating a summation of one or more close contact potentials of the parent protein based on respective distances between the one or more residue atoms of the first residue of the parent protein and one or more environment atoms of the one or more residue atoms of the first residue of the parent protein, in accordance with a close contact potential formula of the form:
• U
• a
• ⁢
• b
• ⁡
• (
• r
• a
• ⁢
• b
• )
• ⁢
• =
• {
• 1
• if
• ⁢
• ⁢
• r
• a
• ⁢
• b
• <
• r
• c
• 1
• d
• ⁢
• (
• r
• c
• -
• r
• a
• ⁢
• b
• )
• ⁢
• if
• ⁢
• ⁢
• r
• c
• <
• r
• a
• ⁢
• b
• <
• r
• c
• +
• d
• 0
• if
• ⁢
• ⁢
• r
• a
• ⁢
• b
• <
• r
• c
• +
• d
• 0
• if
• ⁢
• ⁢
• β
• wherein,
• Uab(rab) is the close contact potential between a residue atom a of the first residue and an environmental atom b,
• d is a non-zero offset,
• rc, is equal to or greater than a first distance threshold,
• rab is the distance between a residue atom a of the first residue and an environmental atom b,
• βrefers to atoms a and b that are either members of the same residue or that are bonded either directly (1-2 bonded), through one intervening atom (1-3 bonded) or through two intervening atoms (1-4 bonded),
• (b) using a computer system to calculate a first contact area of the one or more residue atoms of the first residue of the parent protein that are exposed to the one or more environment atoms as a difference between (i) a total solvent accessible surface area of the one or more residue atoms of the first residue of the parent protein and (ii) a solvent accessible area of the one or more residue atoms of the first residue of the parent protein in the presence of the parent protein, wherein the method further comprises calculating the total solvent accessible surface area of the one or more residue atoms of the first residue of the parent protein and calculating the solvent accessible area of the one or more residue atoms of the first residue of the parent protein in the presence of the parent protein using the computer system,
• (c) calculating a first close contact surface density (CCSD) by dividing the summation of the one or more close contact potentials of the parent protein by the first contact area,
• (d) mutating in silico one or more residues of the first protein to provide the variant protein,
• (e) calculating a summation of one or more close contact potentials of the variant protein based on respective distances between the one or more residue atoms of the first residue of the variant protein and one or more environment atoms of the one or more residue atoms of the first residue of the variant protein in accordance with the close contact potential formula of the calculating (a),
• (f) calculating a second contact area of the one or more residue atoms of the first residue of the variant protein that are exposed to the one or more environment atoms as a difference between (i) a total solvent accessible surface area of the one or more residue atoms of the first residue of the variant protein and (ii) the solvent accessible area of the one or more residue atoms of the first residue of the variant protein in the presence of the variant protein,
• (g) calculating a second CCSD by dividing the summation of the one or more close contact potentials of the variant protein by the second contact area,
• (h) confirming the variant protein has improved stability relative to the parent protein by comparing the second CCSD to the first CCSD;
• and
• (i) making the variant protein, wherein the variant protein is other than a wild-type protein.