A systematic computational study is carried out to investigate hydrogen bond (HB) interactions in the real crystalline structures of L-cysteine at 30 and 298 K by density functional theory (DFT) calculations of electric field gradient (EFG) tensors at the sites of O-17, N-14, and H-2 nuclei. One-molecule (monomer) and nine-molecule (cluster) models of L-cysteine are created by available crystal coordinates at both temperatures and the EFG tensors are calculated for both models to indicate the effect of HB interactions on the tensors. The calculated EFG tensors at the level of B3LYP and B3PW91 DFT methods and 6-311++G** and cc-pVTZ basis sets are converted to those experimentally measurable nuclear quadrupole resonance (NQR) parameters i.e. quadrupole coupling constants (qcc) and asymmetry parameters (eta(Q)). The evaluated NQR parameters reveal that the EFG tensors of (17)O, (14)N, and (2)H are influenced and show particular trends from monomer to the target molecule in the cluster due to the contribution of target molecule to classic N-H...O, and non-classic S-H...O and S-H...S types of HB interactions. On the other hand, atoms in molecules (AIM) analyses confirm the presence of HB interactions and rationalize the observed EFG trends. The results indicate different contribution of various nuclei to HB interactions in the cluster where O2 and N1 have major contributions. The EFG tensors as well as AIM analysis at the H6 site show that the N1-H6...O2 HB undergoes a significant change from 30 to 298 K where changes in other N-H...O interactions are almost negligible. There is a good agreement between the calculated (14)N NQR parameters and reported experimental data.