Molecules with multiple sites of hydrogen bonding attached to suitable cores tend to crystallize as open networks. The resulting crystals can have the following unusual properties: They can include significant amounts of guest molecules; the guests are typically located in channels and can be exchanged without loss of crystallinity; and the geometry of the networks can change in response to new guests. We have found that DFT calculations can provide accurate simulations of the unusual structure and properties of such materials, represented by crystals of prototypic tetrapyridinone 1. These calculations have yielded three key insights that cannot be obtained directly from experiments. (1) The hypothetical porous network obtained by removing guests from crystals of compound 1 is highly flexible, and its deformations are inherently anisotropic, leading to lengthening or shortening of the channels along the c axis and no significant changes along the a and b axes. (2) Quantitative analysis of the total cohesive energy has revealed that hydrogen bonding within the network makes a dominant contribution, along with interactions of guests with the network. (3) Differences in the overall stability of crystals of compound 1 as the guests are varied do not arise primarily from significant changes in the cohesive energy of the network itself; instead, differences in guest-guest interactions play a key role, resulting from the nature of the guests and constraints imposed by the surrounding network. These insights, together with the results of ab initio molecular dynamics, help explain how hydrogen-bonded networks can be robust yet permit molecular movement that underlies the exchange of guests and adaptive porosity. These insights promise to be of general value to scientists studying ordered molecular materials in which strong directional interactions are prominent.