monomer_binding_rule = Rule(
    name = 'monomer_binding_rule',
    reactants = {'ligand_reactant': free_ligand, 'free_receptor': free_receptor},
    actions = ['ligand_reactant.fc1.add_bond(free_receptor.alpha)'],
    rate_prefix = 'association_constant',