In order to perform a systematic study of the interaction between grain boundaries (GBs) and dislocations using molecular dynamics (MD), several tools need to be available. A combination of computational geometry and MD was used to build the foundations of what we call a virtual laboratory. First, an algorithm to generate GBs on face-centered cubic bicrystals was developed. Two crystals with different orientations are placed together. Then, by applying “microscopic” rigid body translations along the GB plane to one of the crystals and removing overlapping atoms, a set of initial configurations is sampled and a minimum energy configuration is found. Second, to classify the geometry of the GBs a local symmetry type (LST) describing the angular environment of each atom is calculated. It is found that for a given relaxed GB the number of atoms with different LSTs is not very large and that it is possible to find unique geometrical patterns in each GB. For instance, the LSTs of two GBs having the same “macroscopic” configuration but different “microscopic” degrees of freedom can be dissimilar: the configurations with higher GB energy tend to have a higher number of atoms with different LSTs. Third, edge dislocations are introduced into the bicrystals. We see that full edge dislocations split into Shockley partials. Finally, by loading the bicrystals with tensile stresses the edge dislocations are put into motion. Various examples of dislocation-GB interactions in Cu are presented.