An accurate description of hydrogen bonds is essential to identify the determinants of protein stability and function as well as folding and misfolding behavior. We describe a method of using J couplings through hydrogen bonds as ensemble-averaged restraints in molecular dynamics simulations. Applications to the cases of ubiquitin and protein G show that these scalar couplings provide powerful structural information that, when used through the methodology that we present here, enables the description of the geometry and energetics of hydrogen bonds with an accuracy approaching that of high-resolution X-ray structures.