Before revision:  We apply the method to simulate multivalent ligand-receptor interaction systems . Simulation results reveal insights into ligand-receptor binding kinetics that are not available from previously developed equilibrium models .
After revision:  We apply the method to simulate multivalent ligand-receptor interaction systems . Simulation results show that the rejection-free method is more efficient over wide parameter ranges than a related method that relies on rejection sampling .