code: def rmsd_eval(rmsd_params):
      
      specification, sequence, parsed_ind, reference_pdb = rmsd_params
      model = specification(*parsed_ind)
      model.pack_new_sequences(sequence)
      ca, bb, aa = run_profit(model.pdb, reference_pdb, path1=False, path2=False)
      return bb
docs: Builds a model and runs profit against a reference model.
  
      Parameters
      ----------
      rmsd_params
  
      Returns
      -------
      rmsd: float
          rmsd against reference model as calculated by profit.