System-environment interactions are essential in determining charge-transfer (CT) rates and mechanisms. We developed a computationally accessible method, suitable to simulate CT in flexible molecules (i.e., DNA) with hundreds of sites, where the system-environment interactions are explicitly treated with numerical noise modeling of time-dependent site energies and couplings. The properties of the noise are tunable, providing us a flexible tool to investigate the detailed effects of correlated thermal fluctuations on CT mechanisms. The noise is parametrizable by molecular simulation and quantum calculation results of specific molecular systems, giving us better molecular resolution in simulating the system-environment interactions than sampling fluctuations from generic spectral density functions. The spatially correlated thermal fluctuations among different sites are naturally built-in in our method but are not readily incorporated using approximate spectral densities. Our method has quantitative accuracy in systems with small redox potential differences (