Current models of brain networks may potentially be improved by integrating our knowledge of structural connections, within and between circuits, with metrics of functional interactions between network nodes. The former may be obtained from diffusion MRI of white matter (WM), while the latter may be derived by measuring correlations between resting state BOLD signals from pairs of gray matter (GM) regions. From inspection of diffusion MRI data, it is clear that each WM voxel within a 3D image array may be traversed by multiple WM structural tracts, each of which connects a pair of GM nodes. We hypothesized that by appropriately weighting and then integrating the functional connectivity of each such connected pair, the overall engagement of any WM voxel in brain functions could be evaluated. This model introduces a structural constraint to earlier studies of WM engagement and addresses some limitations of previous efforts to relate structure and function. Using concepts derived from graph theory, we obtained spatial maps of WM engagement which highlight WM regions critical for efficient communications across the brain. The distributions of WM engagement are highly reproducible across subjects and depict a notable interdependence between the distribution of GM activities and the detailed organization of WM. Additionally, we provide evidence that the engagement varies over time and shows significant differences between genders. These findings suggest the potential of WM engagement as a measure of the integrity of normal brain functions and as a biomarker for neurological and cognitive disorders.