In this work, we draw connections between the classical Shannon interpolation of bandlimited deterministic signals and the literature on estimating continuous-time random processes from their samples (known in various communities under different names, such as Wiener-Kolmogorov filtering, Gaussian process regression, and kriging). This leads to the realization that Shannon interpolation can be interpreted as implicitly expecting that the unknown signal has uniform spectral density within its bandwidth. However, in many practical applications, we often expect more energy at some frequencies than at others. This leads us to propose novel Shannon-like interpolators that are optimal with respect to appropriately-constructed weighted Hilbert spaces, where weighting enables us to accommodate prior information about nonuniform spectral density. Although our new interpolants are derived for data obtained with any sampling rate, we observe that they are particularly useful for interpolating sub-Nyquist data. In addition to theory, we also discuss aspects of practical implementation and show illustrative examples to demonstrate the potential benefits of the approach.