Extracellular matrix (ECM) remodeling is a key physiological process that occurs in a number of contexts, including cell migration, and is especially important for cellular form and function in three-dimensional (3D) matrices. However, there have been few attempts to computationally model how cells modify their environment in a manner that accounts for both cellular properties and the architecture of the surrounding ECM. To this end, we have developed and validated a novel model to simulate matrix remodeling that explicitly defines cells in a 3D collagenous matrix. In our simulation, cells can degrade, deposit, or pull on local fibers, depending on the fiber density around each cell. The cells can also move within the 3D matrix. Different cell phenotypes can be modeled by varying key cellular parameters. Using the model we have studied how two model cancer cell lines, of differing invasiveness, modify matrices with varying fiber density in their vicinity by tracking the metric of fraction of matrix occupied by fibers. Our results quantitatively demonstrate that in low density environments, cells deposit more collagen to uniformly increase fibril fraction. On the other hand, in higher density environments, the less invasive model cell line reduced the fibril fraction as compared to the highly invasive phenotype. These results show good qualitative and quantitative agreement with existing experimental literature. Our simulation is therefore able to function as a novel platform to provide new insights into the clinically relevant and physiologically critical process of matrix remodeling by helping identify critical parameters that dictate cellular behavior in complex native-like environments.