The occurrence of large domain motions associated with the mechanism of action of many proteins is well established. We present a general method of predicting domain closure applicable to proteins containing domains separated by an apparent hinge. The method attempts to allow for natural directional bias within the closing protein by repeatedly applying a weak pulling force over a short distance between pairs of atoms chosen at random in the two domains in question. Appropriate parameters governing the pulling function were determined empirically. The method was applied to the bi-lobal protein PGK and a closed-form activated ternary complex generated for Bacillus stearothermophilus PGK. This model was compared with the recently determined crystal structure of closed-form Trypanosoma brucei PGK. The model predicts the correct hinge regions, although the magnitude of movement at one hinge point was overestimated, and provides a reasonable representation of the closed-form ternary complex.