We present a joint experimental and computational terahertz (THz) spectroscopy study of the most stable polymorph (form I) of an antihypertensive pharmaceutical solid, felodipine (FLD). The vibrational response has been analyzed at room temperature by combining optical (THz-TDS, FT-IR, THz-Raman) and neutron (INS) terahertz spectroscopy. With the challenging example of a large and flexible molecular solid, we illustrate the complementarity of the experimental techniques. We show how the results can be understood by employing ab initio modeling and discuss current progress in the field. To this end, we employ plane wave formulation of density functional theory (plane wave DFT) along with harmonic lattice dynamics calculations (HLD) and ab initio molecular dynamics (AIMD) simulations. Based on a comprehensive theoretical analysis, we discover an inconsistency in the commonly accepted structural model, which can be linked to a distinct librational dynamics of the side ester chains. As a result, only a moderate agreement with the experimental spectra can be achieved. We, therefore, propose an alternative structural model, effectively accounting for the influence of the large-amplitude librations and allowing for a comprehensive analysis of the vibrational resonances up to 4.5 THz. In that way, we illustrate the applicability of the computationally supported THz spectroscopy to detect subtle structural issues in molecular solids. While the provided structural model can be treated as a guess, the problem calls for further revision by means of high-resolution crystallography. The problem also draws a need of extending the THz experiments toward low-temperature conditions and single-crystal samples. On the other hand, the studied system emerges as a challenge for the DFT modeling, being extremely sensitive to the level of the theory used and the resulting description of the intermolecular forces. FLD form I can be, hence, considered as a testbed for the