next up previous
Next: [3] read.expr Up: CHNOSZ examples Previous: [2] ionize.aa

more.aa

## Don't show: data(thermo)
thermo$obigt: 1809 aqueous, 3368 total species
## End(Don't show) # the first 13 names in UniProt for "aminotransferase ecoli" at.ecoli <- c("BIOA", "ARNB", "SERC", "AAT", "TYRB", "ARGD", "ILVE", "ALAA", "ALAC", "YBDL", "AVTA", "GLMS", "PUUE") # get the amino acid compositions # note that ALAA and ALAC are not matched at.aa <- more.aa(at.ecoli, "Eco")
more.aa: ALAA ALAC were not matched
# what are their formulas and oxidation states protein.formula(at.aa)
C H N O S BIOA_Eco 2103 3279 577 604 32 ARNB_Eco 1861 2946 518 539 14 SERC_Eco 1781 2762 478 527 15 AAT_Eco 1936 3028 536 584 13 TYRB_Eco 1950 3048 528 571 16 ARGD_Eco 1964 3043 535 569 16 ILVE_Eco 1523 2357 417 452 11 ALAA_Eco NA NA NA NA NA ALAC_Eco NA NA NA NA NA YBDL_Eco 1937 2990 524 560 12 AVTA_Eco 2114 3311 555 594 22 GLMS_Eco 2958 4746 832 897 18 PUUE_Eco 1990 3138 546 594 16
ZC(protein.formula(at.aa))
[1] -0.1312411 -0.1536808 -0.1301017 -0.1167286 -0.1487179 -0.1364562 -0.1181878 NA NA [10] -0.1414559 -0.1958373 -0.1419878 -0.1407035


next up previous
Next: [3] read.expr Up: CHNOSZ examples Previous: [2] ionize.aa