BACKGROUND
two extraordinary families of proteins have evolved to influence ice growth in opposite ways: antifreeze proteins  that irreversibly adsorb to the surface of ice crystals to prevent their further growth  <cit> ; and ice-nucleation proteins  that cause ice to form in solution at high sub-zero temperatures  <cit> .
whereas afps are small , and generally monomeric proteins, inps are large  and function as multimers  <cit> .
the tertiary structures of many afps are known, but none has been experimentally determined for an inp.
most inps contain three distinct domains, with the majority of their mass residing within a highly repetitive central domain  <cit> .
this domain consists of a variable number  of tandem 16-amino-acid  repeats, with each repeat following the general consensus sequence of gygstxtaxxxsxlxa  <cit> .
nmr and cd studies of synthetic bacterial inp peptides have not revealed a basic folding unit  <cit> .
molecular models of inp from the bacterium pseudomonas syringae have included a planar array of anti-parallel β-strands  <cit>  and a left-handed β-helix  <cit> .
the latter was modelled on udp-acetylglucosamine acyltransferease as an afp-like β-helix with the β-stranded txt motifs located within each 16-aa repeat aligning down one side of the protein and functioning as the site of ice nucleation.
interestingly, a 96-aa recombinantly-expressed fragment of psinp was shown to produce moderate levels of afp activity  <cit> , hinting that inps and afps share a similar mechanism of action.
the gram-negative bacterium p. borealis produces a 1244-aa inp  similar to all other known bacterial inps  <cit> .
here we have predicted the structure of a 128-aa segment of the protein using a combination of homology-based modelling and molecular dynamics  simulations.
the right-handed β-helical model is stabilized by internal serine and glutamine ladders, and is predicted to dimerize via a highly conserved solvent-exposed tyrosine ladder.
ach chain of the dimer contains two putative ice-nucleation sites, located opposite one another, and comprised of repetitive tqta and slta β-strands.
each surface is flat and relatively hydrophobic, but also replete with hydrogen bond donors and acceptors; hallmarks of a typical afp ice-binding site .
md simulations show each site is capable of ordering water molecules into an ice-like lattice, indicating that inps use the same anchored clathrate water  mechanism of action that was recently elucidated for all afps  <cit> .
indeed, acws align across the entire width of both sides of the inp dimer, dramatically increasing the 'active' surface area of the protein and further strengthening the idea that size, and not fold or mechanism of action, is the primary discriminating factor between the two families of proteins with diametrically opposite functions regarding ice growth.
RESULTS
pbinp contains two ice-nucleating motifs per 16-aa repeat
the highly repetitive central domain of pbinp consists of  <dig> tandem 16-aa repeats .
each repeat contains four distinct tetra-peptide motifs; two putative β-strands with consensus sequences of tqta and slta, and two glycine- and serine-rich turns with consensus sequences of xxxs and gygs .
each 16-aa repeat forms part of a higher order 48-aa repeat, consisting of three tandem 16-aa sequences .
the tetra-peptide tqta and lta β-strands  are highly similar to the txt motifs that define the ice-binding sites  of β-helical insect and grass afps  <cit> .
previous modelling of the homologous p. syringae inp  based on udp-acetylglucosamine acyltransferase  predicted each 16-aa repeat of the protein formed one loop of a left-handed β-helix  <cit> .
this aligned the tqta motifs down one side of he protein, forming a flat β-sheet that functioned as the site of ice nucleation.
since the 16-aa repeats of pbinp  contain two putative txt-like β-strands separated by glycine- and serine-rich tetra-peptide turns, each strand of pbinp was predicted to align on opposite sides of a β-helical structure, forming two parallel β-sheets.
for this reason, the β-roll from alkaline protease  was chosen here as a starting template from which to model pbinp.
located between its ca2+-binding turns are two parallel β-sheets packed against one another .
the model of pbinp was built by excising this central portion of the β-roll from alkaline protease and changing its residues to the appropriate tqta and slta motifs .
the xxxs and gygs loops were manually built to connect the strands , and the newly built 16-aa loop was duplicated, substituted with the appropriate pbinp sequence, and aligned to the ensuing β-strands of the structure.
this was performed until an eight-loop parallel β-helix was built that corresponds to residues 217- <dig> of binp .
pbinp folds as a parallel β-helix
pbinp was modelled as a right-handed β-helix due to the right-handedness of the modelling template.
mpafp_riv was previously modelled as a right-handed β-helix using the same β-roll template from alkaline protease and proved to be right-handed based on its crystal structure  <cit> .
however, pbinp could just as likely be left-handed, and the handedness of an ice-binding protein does not affect its ability to bind ice.
both l- and d-enantiomers of type i fish afp and snowflea afp display identical levels of antifreeze activity  <cit> , while the non-homologous β-helical insect afps produced by tenebrio molitor  <cit>  and choristoneura fumiferana  <cit>  are right- and left-handed respectively, yet display near identical ice-binding sites and activity levels.
therefore, while we modelled pbinp as a right-handed β-helix, we think that either a left- or right-handed version of the protein would be equally capable of efficiently nucleating ice.
to test the stability of the initial pbinp model, a solvated 5-ns md simulation was performed.
the model stabilized at the 2-ns mark of the simulation and remained remarkably constant throughout the remainder of the trajectory .
an energy-minimized average structure of the model was calculated from the simulation's final  <dig> ns.
the β-helix is roughly  <dig> Å in length ×  <dig> Å wide × 5- <dig> Å thick .
a ramachandran plot of the structure shows that all residues fall within allowable regions .
the protein contains a modest hydrophobic core, consisting of the alanine methyl groups at the end of each β-strand , the leucine side chain from the slta β-strands, and the side-chain methylene groups of the inward-projecting glutamine and serine residues from the tqta and xxxs motifs respectively .
the leucine residues interdigitate with the side-chain methylene groups of the glutamine residues, further increasing stabilization.
the polar groups of the internal glutamine and serine residues point toward the turns of the protein and away from the hydrophobic core, and main-chain inter-loop hydrogen bonds shield the hydrophobic core from the solvent.
the inward-pointing serine and glutamine residues each form a ladder in opposite corners of the structure.
the serines of the ladder reside in the final position of the xxxs motif that connects the tqta β-strand to the slta β-strand .
each ser hydroxyl group hydrogen bonds to main-chain carbonyl oxygens and amide nitrogens of the same and ensuing xxxs motifs of the structure, while the main-chain amide nitrogen of each serine hydrogen bonds with the main-chain carbonyl oxygen of the second x residue of the xxxs motif of the ensuing loop .
the three x residues of the xxxs motif adopt an xbl β-arcade as defined by hennetin et al.
<cit> .
β-arcades of this configuration are typically followed by residues with short polar side chains that point inwards, allowing them to form stablizing h-bonds with the protein's main chain.
this is clearly occurring in pbinp with the highly conserved inward-pointing row of ser residues.
also, the hyperactive β-helical afp from the beetle tenebrio molitor  uses an internal serine ladder to stabilize turns within its structure  <cit> .
the glutamine ladder is formed by the inward-pointing gln residue of each tqta motif .
glutamines typically repeat two out of every three loops at this position, with the third loop usually having a glycine substitution.
glutamine side chains form hydrogen bonds with main-chain carbonyl oxygens and amide nitrogens of the turn immediately preceding the tqta motif, as shown in stereo in figure 4b.
the void produced by the regularly-spaced glycine substitutions is filled by an internal water molecule that bridges adjacent glutamines and also forms hydrogen bonds to the main chain of the protein .
the hyperactive β-helical tmafp uses internal water molecules to stabilize its structure  <cit> .
here it should be noted that these proteins have evolved to function at sub-zero temperatures where hydrogen bonding is relatively strong.
while no structure in the protein data bank contains a glutamine ladder, several β-helical proteins use very similar asparagine ladders to stabilize their folds  <cit> , and tiba, a 104-kda bacterial glycoprotein with both adhesin and invasin properties, is also predicted to contain an internal glutamine ladder  <cit> .
pbinp dimerizes via a solvent-exposed tyrosine ladder
the tightly packed core of the protein left no room for the tyrosine residue of the gygs motifs, and as such, they were forced to adopt a solvent-exposed orientation during initial model building.
during the course of the md simulation, the tyrosine residues stacked on top of one another , and each tri-peptide ygs motif adopted a β-stranded conformation.
stacking of solvent-exposed aromatic residues has been observed within various proteins, and they typically reside in areas involved in receptor binding and/or dimerization.
for example, inij, a β-helical leucine-rich repeat produced by the bacterium listeria monocytogenes, contains multiple solvent-exposed and stacked aromatic residues within its predicted receptor recognition domain  <cit> .
most interestingly, engineered tyrosine ladders have been shown to promote the flatness  <cit>  and dimerization  <cit>  of ospa, a single-layer-β-sheet containing protein from the bacterium borrelia burgdorferi.
these results hinted at the possibility that pbinp's tyrosine ladder might play a role in dimerization.
to test this idea, a model dimer of pbinp was built by aligning two identical chains of pbinp in a parallel manner such that their tyrosine ladders were oriented towards one another and their n termini were at the same end of the structure.
this placed the slta surface of one chain on the same side of the structure as the tqta surface from the opposite chain.
a 10-ns md simulation was then performed where water molecules were allowed to equilibrate at the dimerization interface prior to the start of the trajectory .
immediately following the start of the simulation, the two chains moved towards one another, driven by the desolvation of the aromatic rings of the tyrosine residues .
the dimeric structure stabilized at the 5-ns mark of the trajectory and remained completely stable thereafter .
the energy-minimized average dimeric structure is shown in figure  <dig>  all waters were excluded from the dimerization interface as the aromatic rings of each tyrosine residue packed tightly against the β-stranded glycine residues of the opposite chain's gygs motif .
interestingly, the hydroxyl group of each tyrosine formed part of an intercalated hydrogen-bond network with the side-chain hydroxyl groups of the outward-projecting serine residues of the opposite chain's gygs motifs .
this is exemplified by the number of side-chain hydrogen bonds that formed during the simulation, increasing from  <dig> at the start to an average of  <dig> upon stabilization of the dimerization interface .
to test the stability of the dimer, 10-ns md simulations were performed on the energy-minimized average dimeric structure at the elevated temperatures of  <dig> k and  <dig> k .
in each case, the dimer remained stable throughout the trajectory, with no dissociation of individual chains, and complete exclusion of water molecules at the dimerization interface.
identical md simulations were performed on the pbinp monomer, and in each case, a partial unravelling of the c terminus of the protein occurred  .
these results demonstrate that dimerization increases the stability of pbinp.
we also attempted to model pbinp as an anti-parallel dimer .
the model remained intact following a 10-ns solvated md simulation, however, a twist developed in the structure that did not develop in the parallel dimer .
in an anti-parallel orientation, the hydroxyl groups of the solvent-exposed rank of serine residues from one chain are unable to form an intercalated hydrogen-bond network with the hydroxyl groups of the tyrosine residues of the opposite chain.
this lack of interaction allowed the twist to occur and therefore spoiled the flatness of the dimer.
flatness is a defining characteristic of ice-binding proteins, and as such, we think it is unlikely that pbinp would dimerize in an anti-parallel manner.
the ice nucleating surfaces of pbinp are highly conserved
as previously mentioned, when pbinp is modelled as a parallel dimer, the slta surface of one chain resides on the same side of the dimer as the tqta surface of the other chain .
this creates a flat surface on both sides of the dimer, each spanning the entire width of the structure.
the amino acid composition of both the tqta and slta surfaces is highly conserved, as is the area in between them created by the dimerization interface.
more specifically, the solvent-exposed first and third position of the tqta motifs have threonines present at these positions 87% and 92% of the time respectively.
the solvent-exposed first and third position of the slta motifs show greater variation, however, serines and threonines are present at the first position 71% of the time , while threonine and isoleucine are present at the third position 98% of the time .
the rare substitutions at these positions are aspartate, valine and leucine.
the dimerization interface is the most highly conserved portion of the whole protein, with the intercalated hydroxyl groups of the tyrosine and serine residues present 90% and 98% of the time, respectively.
pbinp orders water molecules via the anchored clathrate water mechanism
this high degree of amino acid conservation endows the pbinp dimer with two sides that are flat, and relatively hydrophobic, but also contain many hydrogen bond donors and acceptors.
these characteristics define the ibs's of all afps, and as such, raise the possibility that inps function by the acw mechanism of afp action  <cit> .
the acw mechanism states that the relative hydrophobicity of an afp's ibs orders water molecules into an ice-like lattice, and this lattice is then anchored to the surface of the protein via hydrogen bonds.
to investigate the water ordering potential of pbinp, a 2-ns md simulation was performed on the energy-minimized average dimeric structure at a temperature of  <dig> k. the tip5p water model was used as it accurately represents the behaviour of water during an md simulation  <cit> .
waters within  <dig> Å of pbinp's surface were extracted at 20-ps intervals throughout the course of the simulation and their crystallographic structure factors were calculated followed by a fourier transform, giving the electron density of all extracted waters.
strong electron density with ice-like spacing was present across the entire width of the protein's surface .
waters built into this density aligned down the troughs created by the outward projecting residues of the tqta and slta motifs, down the flat areas immediately preceding and following these motifs, and across the dimerization interface as well .
these waters closely match the  <dig>  Å ×  <dig>  Å spacing of waters on the primary prism plane of ice, and as such, have revealed a potential orientation of a nascent ice crystal on the surface of the inp .
as a positive control, the same simulation was performed on mpafp_riv, a hyperactive ca2+-dependent β-helical afp produced by the antarctic bacterium marinomonas primoryensis.
mpafp_riv was chosen as the positive control because the positions of the acws on its ibs are known  <cit> .
the simulation accurately predicted the position of the ice-like waters on the ibs of the protein .
indeed, aligning  <dig> manually-built waters to their crystallographic counterparts on the ibs of chain b of mpafp_riv produced an rmsd of only  <dig>  Å. further proof of the accuracy of the md simulation is evidenced by the identical coordination of waters present down the trough created by the tqta motifs of pbinp as compared to the waters located down the ibs txt troughs of the β-helical insect afps  as revealed by x-ray crystallography .
dimerization increases the active surface area of pbinp
the temperature at which an inp nucleates ice is dependent upon its oligomerization state  <cit> .
approximately one inp monomer is capable of nucleating ice at a temperature of -12°c, while an aggregate of at least  <dig> inp monomers is required for ice nucleation in the -2°c to -3°c temperature range.
however, due to the lack of an experimentally-determined inp structure, the mechanism by which an inp oligomerizes has remained a mystery.
previously, there have been only two attempts to predict how an inp might oligomerize.
wu et al.
<cit>  speculated that overlapping protein-protein interactions between the two flat ibs motifs could generate 'stairs' of inps, facilitating ice growth along the discontinuities.
prior to that, however, kajava and lindow  <cit>  modelled psinp as an array of interdigitating anti-parallel β-strands that formed a flat ice-nucleating array upon oligomerization.
while intriguing, the model was created prior to the structural determination of several β-helical afps  <cit> , and as such, it overlooked the potential water-organizing capabilities of the tqta and slta motifs when aligned as parallel β-strands in a flat β-sheet.
as previously mentioned, graether and jia  <cit>  modelled an inp as a β-helix, but the potential for oligomerization was not discussed.
in that model, the serine of the gygs motif pointed inwards and was not solvent exposed as predicted in this study.
this would prevent formation of the intercalated hydrogen-bond network between the serine and tyrosine residues of the gygs motifs from opposing chains in the dimer.
nevertheless, both β-helical models predict a solvent-exposed orientation for the tyrosine residue of the gygs motif, and this argues against inp peptide studies that suggest it points towards the interior of the structure  <cit> .
it is likely that short peptides adopt conformations in solution that are not representative of the full-length protein.
as a dimer with two bridged active sites on opposite sides of the structure, the potential ice-nucleating surface area of pbinp is dramatically increased .
a full-length pbinp dimer would have an active surface area of ca.
<dig>   <dig> Å <dig> , a value slightly larger than the  <dig>   <dig> Å <dig> minimum surface area required for an ice embryo to continue growing at -12°c  <cit>  .
by dimerizing in an offset manner, pbinp monomers could form extended oligomeric structures that would increase the active surface area of the complex without occluding it, therefore raising the temperature at which ice nucleation occurs .
in this scenario, the repetitive central domain of each pbinp chain runs parallel to the surface of the bacterium, with the non-repetitive n terminus of each chain serving as a membrane anchor.
multiple arrays of offset dimers could increase the active pbinp surface area on the outer membrane of the bacterium.
following initial construction of the monomeric model of pbinp, it was thought the solvent-exposed tyrosine ladder might play a role in anchoring the protein to the membrane.
however, several studies have shown the n-terminal domain alone is capable of performing this task  <cit> , and it therefore would seem redundant to have other segments of the protein anchored to the bacterium, which would restrict the protein's access to solvent.
further investigation into the tertiary and quaternary structures of bacterial inps, both in vitro and in vivo, are required and this oligomerization scenario is just one theoretical way of extending the ice nucleation surface.
even without oligomerization, the size of pbinp is significantly larger than any known afp .
this fact, combined with an inp's ability to order water molecules via the acw mechanism, further strengthens the idea that size, and not a particular protein fold or mechanism of action, is the primary discriminating factor between the two families of proteins with diametrically opposite functions.
CONCLUSIONS
this paper demonstrates through the use of homology modelling and md simulations that bacterial inps are able to fold as novel β-helical dimers.
de-solvation of solvent-exposed tyrosine ladders drives inp dimerization, and this allows the ice nucleation sites of the protein to extend as a continuum across the width of the dimer.
both sides of the dimer order water molecules into an ice-like lattice using the anchored clathrate water mechanism of action.
offset dimerization can allow inp oligomerization without active site occlusion.
this increases the active surface area of the inp, therefore raising the temperature at which ice nucleates.

 two extraordinary families of proteins have evolved to influence ice growth in opposite ways: antifreeze proteins that irreversibly adsorb to the surface of ice crystals to prevent their further growth; and ice-nucleation proteins that cause ice to form in solution at high sub-zero temperatures. whereas afps are small, and generally monomeric proteins, inps are large and function as multimers. the tertiary structures of many afps are known, but none has been experimentally determined for an inp. most inps contain three distinct domains, with the majority of their mass residing within a highly repetitive central domain. this domain consists of a variable number of tandem 16-amino-acid repeats, with each repeat following the general consensus sequence of gygstxtaxxxsxlxa. nmr and cd studies of synthetic bacterial inp peptides have not revealed a basic folding unit. molecular models of inp from the bacterium pseudomonas syringae have included a planar array of anti-parallel β-strands and a left-handed β-helix. the latter was modelled on udp-acetylglucosamine acyltransferease as an afp-like β-helix with the β-stranded txt motifs located within each 16-aa repeat aligning down one side of the protein and functioning as the site of ice nucleation. interestingly, a 96-aa recombinantly-expressed fragment of psinp was shown to produce moderate levels of afp activity, hinting that inps and afps share a similar mechanism of action. the gram-negative bacterium p. borealis produces a 1244-aa inp similar to all other known bacterial inps. here we have predicted the structure of a 128-aa segment of the protein using a combination of homology-based modelling and molecular dynamics simulations. the right-handed β-helical model is stabilized by internal serine and glutamine ladders, and is predicted to dimerize via a highly conserved solvent-exposed tyrosine ladder. ach chain of the dimer contains two putative ice-nucleation sites, located opposite one another, and comprised of repetitive tqta and slta β-strands. each surface is flat and relatively hydrophobic, but also replete with hydrogen bond donors and acceptors; hallmarks of a typical afp ice-binding site. md simulations show each site is capable of ordering water molecules into an ice-like lattice, indicating that inps use the same anchored clathrate water mechanism of action that was recently elucidated for all afps. indeed, acws align across the entire width of both sides of the inp dimer, dramatically increasing the'active' surface area of the protein and further strengthening the idea that size, and not fold or mechanism of action, is the primary discriminating factor between the two families of proteins with diametrically opposite functions regarding ice growth.  pbinp contains two ice-nucleating motifs per 16-aa repeat the highly repetitive central domain of pbinp consists of tandem 16-aa repeats. each repeat contains four distinct tetra-peptide motifs; two putative β-strands with consensus sequences of tqta and slta, and two glycine- and serine-rich turns with consensus sequences of xxxs and gygs. each 16-aa repeat forms part of a higher order 48-aa repeat, consisting of three tandem 16-aa sequences. the tetra-peptide tqta and lta β-strands are highly similar to the txt motifs that define the ice-binding sites of β-helical insect and grass afps. previous modelling of the homologous p. syringae inp based on udp-acetylglucosamine acyltransferase predicted each 16-aa repeat of the protein formed one loop of a left-handed β-helix. this aligned the tqta motifs down one side of he protein, forming a flat β-sheet that functioned as the site of ice nucleation. since the 16-aa repeats of pbinp contain two putative txt-like β-strands separated by glycine- and serine-rich tetra-peptide turns, each strand of pbinp was predicted to align on opposite sides of a β-helical structure, forming two parallel β-sheets. for this reason, the β-roll from alkaline protease was chosen here as a starting template from which to model pbinp. located between its ca2+-binding turns are two parallel β-sheets packed against one another. the model of pbinp was built by excising this central portion of the β-roll from alkaline protease and changing its residues to the appropriate tqta and slta motifs. the xxxs and gygs loops were manually built to connect the strands, and the newly built 16-aa loop was duplicated, substituted with the appropriate pbinp sequence, and aligned to the ensuing β-strands of the structure. this was performed until an eight-loop parallel β-helix was built that corresponds to residues 217- of binp. pbinp folds as a parallel β-helix pbinp was modelled as a right-handed β-helix due to the right-handedness of the modelling template. mpafp_riv was previously modelled as a right-handed β-helix using the same β-roll template from alkaline protease and proved to be right-handed based on its crystal structure. however, pbinp could just as likely be left-handed, and the handedness of an ice-binding protein does not affect its ability to bind ice. both l- and d-enantiomers of type i fish afp and snowflea afp display identical levels of antifreeze activity, while the non-homologous β-helical insect afps produced by tenebrio molitor and choristoneura fumiferana are right- and left-handed respectively, yet display near identical ice-binding sites and activity levels. therefore, while we modelled pbinp as a right-handed β-helix, we think that either a left- or right-handed version of the protein would be equally capable of efficiently nucleating ice. to test the stability of the initial pbinp model, a solvated 5-ns md simulation was performed. the model stabilized at the 2-ns mark of the simulation and remained remarkably constant throughout the remainder of the trajectory. an energy-minimized average structure of the model was calculated from the simulation's final ns. the β-helix is roughly Å in length × Å wide × 5- Å thick. a ramachandran plot of the structure shows that all residues fall within allowable regions. the protein contains a modest hydrophobic core, consisting of the alanine methyl groups at the end of each β-strand, the leucine side chain from the slta β-strands, and the side-chain methylene groups of the inward-projecting glutamine and serine residues from the tqta and xxxs motifs respectively. the leucine residues interdigitate with the side-chain methylene groups of the glutamine residues, further increasing stabilization. the polar groups of the internal glutamine and serine residues point toward the turns of the protein and away from the hydrophobic core, and main-chain inter-loop hydrogen bonds shield the hydrophobic core from the solvent. the inward-pointing serine and glutamine residues each form a ladder in opposite corners of the structure. the serines of the ladder reside in the final position of the xxxs motif that connects the tqta β-strand to the slta β-strand. each ser hydroxyl group hydrogen bonds to main-chain carbonyl oxygens and amide nitrogens of the same and ensuing xxxs motifs of the structure, while the main-chain amide nitrogen of each serine hydrogen bonds with the main-chain carbonyl oxygen of the second x residue of the xxxs motif of the ensuing loop. the three x residues of the xxxs motif adopt an xbl β-arcade as defined by hennetin et al.. β-arcades of this configuration are typically followed by residues with short polar side chains that point inwards, allowing them to form stablizing h-bonds with the protein's main chain. this is clearly occurring in pbinp with the highly conserved inward-pointing row of ser residues. also, the hyperactive β-helical afp from the beetle tenebrio molitor uses an internal serine ladder to stabilize turns within its structure. the glutamine ladder is formed by the inward-pointing gln residue of each tqta motif. glutamines typically repeat two out of every three loops at this position, with the third loop usually having a glycine substitution. glutamine side chains form hydrogen bonds with main-chain carbonyl oxygens and amide nitrogens of the turn immediately preceding the tqta motif, as shown in stereo in figure 4b. the void produced by the regularly-spaced glycine substitutions is filled by an internal water molecule that bridges adjacent glutamines and also forms hydrogen bonds to the main chain of the protein. the hyperactive β-helical tmafp uses internal water molecules to stabilize its structure. here it should be noted that these proteins have evolved to function at sub-zero temperatures where hydrogen bonding is relatively strong. while no structure in the protein data bank contains a glutamine ladder, several β-helical proteins use very similar asparagine ladders to stabilize their folds, and tiba, a 104-kda bacterial glycoprotein with both adhesin and invasin properties, is also predicted to contain an internal glutamine ladder. pbinp dimerizes via a solvent-exposed tyrosine ladder the tightly packed core of the protein left no room for the tyrosine residue of the gygs motifs, and as such, they were forced to adopt a solvent-exposed orientation during initial model building. during the course of the md simulation, the tyrosine residues stacked on top of one another, and each tri-peptide ygs motif adopted a β-stranded conformation. stacking of solvent-exposed aromatic residues has been observed within various proteins, and they typically reside in areas involved in receptor binding and/or dimerization. for example, inij, a β-helical leucine-rich repeat produced by the bacterium listeria monocytogenes, contains multiple solvent-exposed and stacked aromatic residues within its predicted receptor recognition domain. most interestingly, engineered tyrosine ladders have been shown to promote the flatness and dimerization of ospa, a single-layer-β-sheet containing protein from the bacterium borrelia burgdorferi. these results hinted at the possibility that pbinp's tyrosine ladder might play a role in dimerization. to test this idea, a model dimer of pbinp was built by aligning two identical chains of pbinp in a parallel manner such that their tyrosine ladders were oriented towards one another and their n termini were at the same end of the structure. this placed the slta surface of one chain on the same side of the structure as the tqta surface from the opposite chain. a 10-ns md simulation was then performed where water molecules were allowed to equilibrate at the dimerization interface prior to the start of the trajectory. immediately following the start of the simulation, the two chains moved towards one another, driven by the desolvation of the aromatic rings of the tyrosine residues. the dimeric structure stabilized at the 5-ns mark of the trajectory and remained completely stable thereafter. the energy-minimized average dimeric structure is shown in figure all waters were excluded from the dimerization interface as the aromatic rings of each tyrosine residue packed tightly against the β-stranded glycine residues of the opposite chain's gygs motif. interestingly, the hydroxyl group of each tyrosine formed part of an intercalated hydrogen-bond network with the side-chain hydroxyl groups of the outward-projecting serine residues of the opposite chain's gygs motifs. this is exemplified by the number of side-chain hydrogen bonds that formed during the simulation, increasing from at the start to an average of upon stabilization of the dimerization interface. to test the stability of the dimer, 10-ns md simulations were performed on the energy-minimized average dimeric structure at the elevated temperatures of k and k. in each case, the dimer remained stable throughout the trajectory, with no dissociation of individual chains, and complete exclusion of water molecules at the dimerization interface. identical md simulations were performed on the pbinp monomer, and in each case, a partial unravelling of the c terminus of the protein occurred. these results demonstrate that dimerization increases the stability of pbinp. we also attempted to model pbinp as an anti-parallel dimer. the model remained intact following a 10-ns solvated md simulation, however, a twist developed in the structure that did not develop in the parallel dimer. in an anti-parallel orientation, the hydroxyl groups of the solvent-exposed rank of serine residues from one chain are unable to form an intercalated hydrogen-bond network with the hydroxyl groups of the tyrosine residues of the opposite chain. this lack of interaction allowed the twist to occur and therefore spoiled the flatness of the dimer. flatness is a defining characteristic of ice-binding proteins, and as such, we think it is unlikely that pbinp would dimerize in an anti-parallel manner. the ice nucleating surfaces of pbinp are highly conserved as previously mentioned, when pbinp is modelled as a parallel dimer, the slta surface of one chain resides on the same side of the dimer as the tqta surface of the other chain. this creates a flat surface on both sides of the dimer, each spanning the entire width of the structure. the amino acid composition of both the tqta and slta surfaces is highly conserved, as is the area in between them created by the dimerization interface. more specifically, the solvent-exposed first and third position of the tqta motifs have threonines present at these positions 87% and 92% of the time respectively. the solvent-exposed first and third position of the slta motifs show greater variation, however, serines and threonines are present at the first position 71% of the time, while threonine and isoleucine are present at the third position 98% of the time. the rare substitutions at these positions are aspartate, valine and leucine. the dimerization interface is the most highly conserved portion of the whole protein, with the intercalated hydroxyl groups of the tyrosine and serine residues present 90% and 98% of the time, respectively. pbinp orders water molecules via the anchored clathrate water mechanism this high degree of amino acid conservation endows the pbinp dimer with two sides that are flat, and relatively hydrophobic, but also contain many hydrogen bond donors and acceptors. these characteristics define the ibs's of all afps, and as such, raise the possibility that inps function by the acw mechanism of afp action. the acw mechanism states that the relative hydrophobicity of an afp's ibs orders water molecules into an ice-like lattice, and this lattice is then anchored to the surface of the protein via hydrogen bonds. to investigate the water ordering potential of pbinp, a 2-ns md simulation was performed on the energy-minimized average dimeric structure at a temperature of k. the tip5p water model was used as it accurately represents the behaviour of water during an md simulation. waters within Å of pbinp's surface were extracted at 20-ps intervals throughout the course of the simulation and their crystallographic structure factors were calculated followed by a fourier transform, giving the electron density of all extracted waters. strong electron density with ice-like spacing was present across the entire width of the protein's surface. waters built into this density aligned down the troughs created by the outward projecting residues of the tqta and slta motifs, down the flat areas immediately preceding and following these motifs, and across the dimerization interface as well. these waters closely match the Å × Å spacing of waters on the primary prism plane of ice, and as such, have revealed a potential orientation of a nascent ice crystal on the surface of the inp. as a positive control, the same simulation was performed on mpafp_riv, a hyperactive ca2+-dependent β-helical afp produced by the antarctic bacterium marinomonas primoryensis. mpafp_riv was chosen as the positive control because the positions of the acws on its ibs are known. the simulation accurately predicted the position of the ice-like waters on the ibs of the protein. indeed, aligning manually-built waters to their crystallographic counterparts on the ibs of chain b of mpafp_riv produced an rmsd of only Å. further proof of the accuracy of the md simulation is evidenced by the identical coordination of waters present down the trough created by the tqta motifs of pbinp as compared to the waters located down the ibs txt troughs of the β-helical insect afps as revealed by x-ray crystallography. dimerization increases the active surface area of pbinp the temperature at which an inp nucleates ice is dependent upon its oligomerization state. approximately one inp monomer is capable of nucleating ice at a temperature of -12°c, while an aggregate of at least inp monomers is required for ice nucleation in the -2°c to -3°c temperature range. however, due to the lack of an experimentally-determined inp structure, the mechanism by which an inp oligomerizes has remained a mystery. previously, there have been only two attempts to predict how an inp might oligomerize. wu et al. speculated that overlapping protein-protein interactions between the two flat ibs motifs could generate'stairs' of inps, facilitating ice growth along the discontinuities. prior to that, however, kajava and lindow modelled psinp as an array of interdigitating anti-parallel β-strands that formed a flat ice-nucleating array upon oligomerization. while intriguing, the model was created prior to the structural determination of several β-helical afps, and as such, it overlooked the potential water-organizing capabilities of the tqta and slta motifs when aligned as parallel β-strands in a flat β-sheet. as previously mentioned, graether and jia modelled an inp as a β-helix, but the potential for oligomerization was not discussed. in that model, the serine of the gygs motif pointed inwards and was not solvent exposed as predicted in this study. this would prevent formation of the intercalated hydrogen-bond network between the serine and tyrosine residues of the gygs motifs from opposing chains in the dimer. nevertheless, both β-helical models predict a solvent-exposed orientation for the tyrosine residue of the gygs motif, and this argues against inp peptide studies that suggest it points towards the interior of the structure. it is likely that short peptides adopt conformations in solution that are not representative of the full-length protein. as a dimer with two bridged active sites on opposite sides of the structure, the potential ice-nucleating surface area of pbinp is dramatically increased. a full-length pbinp dimer would have an active surface area of ca. Å, a value slightly larger than the Å minimum surface area required for an ice embryo to continue growing at -12°c. by dimerizing in an offset manner, pbinp monomers could form extended oligomeric structures that would increase the active surface area of the complex without occluding it, therefore raising the temperature at which ice nucleation occurs. in this scenario, the repetitive central domain of each pbinp chain runs parallel to the surface of the bacterium, with the non-repetitive n terminus of each chain serving as a membrane anchor. multiple arrays of offset dimers could increase the active pbinp surface area on the outer membrane of the bacterium. following initial construction of the monomeric model of pbinp, it was thought the solvent-exposed tyrosine ladder might play a role in anchoring the protein to the membrane. however, several studies have shown the n-terminal domain alone is capable of performing this task, and it therefore would seem redundant to have other segments of the protein anchored to the bacterium, which would restrict the protein's access to solvent. further investigation into the tertiary and quaternary structures of bacterial inps, both in vitro and in vivo, are required and this oligomerization scenario is just one theoretical way of extending the ice nucleation surface. even without oligomerization, the size of pbinp is significantly larger than any known afp. this fact, combined with an inp's ability to order water molecules via the acw mechanism, further strengthens the idea that size, and not a particular protein fold or mechanism of action, is the primary discriminating factor between the two families of proteins with diametrically opposite functions.  this paper demonstrates through the use of homology modelling and md simulations that bacterial inps are able to fold as novel β-helical dimers. de-solvation of solvent-exposed tyrosine ladders drives inp dimerization, and this allows the ice nucleation sites of the protein to extend as a continuum across the width of the dimer. both sides of the dimer order water molecules into an ice-like lattice using the anchored clathrate water mechanism of action. offset dimerization can allow inp oligomerization without active site occlusion. this increases the active surface area of the inp, therefore raising the temperature at which ice nucleates. 
the protein uses internal serine and glutamine ladders for stabilization and is predicted to dimerize via the burying of a solvent-exposed tyrosine ladder to make an intimate hydrophobic contact along the dimerization interface. both sides of the pbinp structure have tandem arrays of amino acids that can organize waters into the ice-like clathrate structures seen on antifreeze proteins. as pbinp is highly similar to all known bacterial inps, we predict its fold and mechanism of action will apply to these other inps. the manner in which pbinp dimerizes also allows for its multimerization, which could explain the aggregation-dependence of inp activity. dimerization dramatically increases the'ice-active' surface area of the protein by doubling its width, increasing its length, and presenting identical ice-forming surfaces on both sides of the protein.
 ice nucleation proteins allow water to freeze at high subzero temperatures. due to their large size, membrane association, and tendency to aggregate, an experimentally-determined tertiary structure of an inp has yet to be reported. how they function at the molecular level therefore remains unknown.  here we have predicted a novel β-helical fold for the inp produced by the bacterium pseudomonas borealis. the protein uses internal serine and glutamine ladders for stabilization and is predicted to dimerize via the burying of a solvent-exposed tyrosine ladder to make an intimate hydrophobic contact along the dimerization interface. the manner in which pbinp dimerizes also allows for its multimerization, which could explain the aggregation-dependence of inp activity. both sides of the pbinp structure have tandem arrays of amino acids that can organize waters into the ice-like clathrate structures seen on antifreeze proteins.  dimerization dramatically increases the'ice-active' surface area of the protein by doubling its width, increasing its length, and presenting identical ice-forming surfaces on both sides of the protein. we suggest that this allows sufficient anchored clathrate waters to align on the inp surface to nucleate freezing. as pbinp is highly similar to all known bacterial inps, we predict its fold and mechanism of action will apply to these other inps. 