We have combined molecular dynamics simulations with first principles calculations to study electron 4 transport in a single molecule of perylene tetracarboxylic diimide (PTCDI) sandwiched between two gold electrodes with an aqueous electrolyte. This combination has for the first time allowed one to reveal statistical behavior of molecular conductance in solution at different temperatures and to produce conductance histograms that can be directly compared with experiments. Our calculations show that experimentally observed temperature-dependent conductance ran be attributed to the thermal effect on the hydrogen bonding network around the molecule and can be described by the radial distribution of water molecules surrounding the oxygen atom in the PTCDI molecule.