To extract topological data of the global electric field, several parameters need to be determined for every system: the volume to calculate the 3-D 𝐸𝐸 "⃗ (𝑟𝑟), the number of streamlines, and the step size. First, the volume and location where the streamlines are computed should be based on a chemically relevant space, and is reaction specific. We define this space as one in which r undergoes significant deformation for a reaction to proceed . For example, in a bond breaking reaction, this volume constitutes the bond itself. However, if a conjugated system is involved, the reaction may involve perturbation of r over the entire region of the conjugation. For the sake of understanding the isolated effect of 𝐸𝐸 "⃗ (𝑟𝑟) that the protein scaffold exerts, we zero charges in the chemically relevant region. For Cytochrome C, we zeroed the charges on the heme and coordinating residues. In Alcohol Dehydrogenase, we zeroed only the substrate, as we wished to compare to an experimental study that probed the effect of the entire protein scaffold and coordinating metal on the substrate . For HIV-1 Protease, we zeroed the catalytic aspartate residues as well as the peptide substrate.