Molecular simulation of adsorption of water molecules in nanoporous amorphous biopolymers, e.g., cellulose, reveals nonlinear swelling and nonlinear mechanical response with the increase in fluid content. These nonlinearities result from hydrogen bond breakage by water molecules. Classical poroelastic models, employing porosity and pore pressure as basic variables for describing the "pore fluid," are not adequate for the description of these systems. There is neither a static geometric structure to which porosity can sensibly be assigned nor arrangements of water molecules that are adequately described by giving them a pressure. We employ molar concentration of water and chemical potential to describe the state of the "pore fluid" and stress-strain as mechanical variables. A thermodynamic description is developed using a model energy function having mechanical, fluid, and fluid-mechanical coupling contributions. The parameters in this model energy are fixed by the output of the initial simulation and validated with the results of further simulation. The poroelastic properties, e.g., swelling and mechanical response, are found to be functions both of the molar concentration of water and the stress. The basic fluid-mechanical coupling coefficient, the swelling coefficient, depends on the molar concentration of water and stress and is interpreted in terms of porosity change and solid matrix deformation. The difference between drained and undrained bulk stiffness is explained as is the dependence of these moduli on concentration and stress.