def test_result_vec_to_features(self): "Test lassopeptides.result_vec_to_features()" orig_feature = DummyCDS(0, 165, strand=1, locus_tag='FAKE0001') seq = "TAILTAILTAILTAILTAILTAILTAILTAILTAILCCTIP" vec = Lassopeptide(23, 42, 51, "HEADHEADHEAD", seq) vec.c_cut = "TIP" motif = result_vec_to_motif(orig_feature, vec) assert motif.locus_tag == "FAKE0001" assert motif.location == orig_feature.location assert motif.leader == "HEADHEADHEAD" assert motif.tail == "TIP" assert motif.core == seq[:-len(motif.tail)] assert motif.detailed_information.num_bridges == 1 assert motif.peptide_subclass == "Class III" assert motif.detailed_information.rodeo_score == 51 assert motif.score == 42 assert motif.detailed_information.macrolactam == "" self.assertAlmostEqual(motif.monoisotopic_mass, 4103.5, delta=1) self.assertAlmostEqual(motif.molecular_weight, 4106.1, delta=1) self.assertAlmostEqual(motif.detailed_information.cut_mass, 3792.3, delta=1) self.assertAlmostEqual(motif.detailed_information.cut_weight, 3794.7, delta=1)
def test_molecular_weight(self): "Test Lassopeptide.molecular_weight" lasso = Lassopeptide(23, 42, 17, 51, "lead", "MAGICHATTIP") lasso.c_cut = "TIP" analysis = ProteinAnalysis("MAGICHATTIP", monoisotopic=False) cut_analysis = ProteinAnalysis("MAGICHAT", monoisotopic=False) mw = analysis.molecular_weight() - 18.02 cut_mw = cut_analysis.molecular_weight() - 18.02 self.assertAlmostEqual(mw, lasso.molecular_weight) self.assertAlmostEqual(cut_mw, lasso.cut_weight)