Recently, two-dimensional (2D) THz-Raman spectroscopy has been used to investigate the intermolecular modes of liquid water. We examine such 2D spectroscopy signals by means of full molecular dynamics (MD) simulations. In this way, we carry out a detailed analysis of intermolecular interactions that play an essential role in many important chemical processes. We calculate 2D Raman-THz-THz (RTT), THz-Raman-THz (TRT), and 2D Raman signals for liquid water, methanol, formamide, acetonitrile, formaldehyde, and dimethyl sulfoxide using an equilibrium-non-equilibrium hybrid MD simulation algorithm originally developed for 2D Raman spectroscopy. These signals are briefly analyzed in terms of anharmonicity and nonlinear polarizability of vibrational modes on the basis of the 2D Raman signals calculated from a Brownian oscillator model with a nonlinear system-bath interaction. We find that the anharmonic contribution is dominant in the RTT case, while the nonlinear polarizability contribution is dominant in the TRT case. For water and methanol, we observed vibrational echo peaks of librational motion in the 2D TRT signals. The predicted signal profiles and intensities that we obtained provide valuable information that can be applied to 2D spectroscopy experiments, allowing them to be carried out more efficiently.