def read_pdbatoms(pdbfile):
    """Read pdbfile and return list of ATOM entries"""
    pdblines = []
    for line in open(pdbfile, 'r'):
        if line[:4] == "ATOM":
            pdblines.append(line)
    return pdblines