Conformational heterogeneity is essential for protein function, yet validating theoretical molecular dynamics (MD) ensembles remains a significant challenge. In this study, we present an approach that integrates free MD simulations, starting from an AlphaFold-generated structure, with refined experimental NMR-relaxation data to identify biologically relevant holistic time-resolved 4D conformational ensembles. Specifically, we select trajectory segments (RMSD plateaus) consistent with experimental observables. For the extracellular region of PsrSp, we found that only specific segments of the long MD trajectory aligned well with experimental data. The resulting ensembles revealed two regions with increased flexibility, both of which play important functional roles.