# Option 1
G = max(nx.connected_component_subgraphs(G), key=len)

# Option 2
for component in list(nx.connected_components(G)):
    if len(component) < 16451: # Size of the largest component
        for node in component:
            G.remove_node(node)