With reference to electromagnetic interaction and Abdus Salam’s strong (nuclear) gravity, 1) Square root of ‘reciprocal’ of the strong coupling constant can be considered as the strength of nuclear elementary charge. 2) ‘Reciprocal’ of the strong coupling constant can be considered as the maximum strength of nuclear binding energy. 3) In deuteron, strength of nuclear binding energy is around unity and there exists no strong interaction in between neutron and proton. $G_s\cong 3.32688\times {10}^{28}{\mathrm{m}}^3{\mathrm{kg}}^{-1}{\sec }^{-2}$ being the nuclear gravitational constant, nuclear charge radius can be shown to be, $R_0\cong \frac{2G_s{m}_p}{c^2}\cong 1.24\mathrm{fm}.$ $e_s\cong \left(\frac{G_s{m}_p^2}{\hslash c}\right)e\cong 4.716785\times {10}^{-19}\mathrm{C}$ being the nuclear elementary charge, proton magnetic moment can be shown to be, ${\mu }_p\cong \frac{e_s\hslash }{2m_p}\cong \frac{e{G}_s{m}_p}{2c}\cong 1.48694\times {10}^{-26}\mathrm{J}{.\mathrm{T}}^{-1}.$ ${\alpha }_s\cong {\left(\frac{\hslash c}{G_s{m}_p^2}\right)}^2\cong 0.1153795$ being the strong coupling constant, strong interaction range can be shown to be proportional to $\exp \left(\frac{1}{{\alpha }_s^2}\right).$ Interesting points to be noted are: An increase in the value of ${\alpha }_s$ helps in decreasing the interaction range indicating a more strongly bound nuclear system. A decrease in the value of ${\alpha }_s$ helps in increasing the interaction range indicating a more weakly bound nuclear system. From $Z\cong 30$ onwards, close to stable mass numbers, nuclear binding energy can be addressed with, ${\left(B\right)}_{A_s}\cong Z\times \left\{\left(\frac{1}{{\alpha }_s}+1\right)+\sqrt{\sqrt{30\times 31}}\right\}\left(m_n-m_p\right)c^2\approx Z\times 19.66\mathrm{MeV}.$ With further study, magnitude of the Newtonian gravitational constant can be estimated with nuclear elementary physical constants. One sample relation is, $\left(\frac{G_N}{G_s}\right)\cong \frac{1}{2}{\left(\frac{m_e}{m_p}\right)}^{10}\left[\sqrt{\frac{G_F}{\hslash c}}/\left(\frac{\hslash }{m_{e}c}\right)\right]$ where $G_N$ represents the Newtonian gravitational constant and $G_F$ represents the Fermi’s weak coupling constant. Two interesting coincidences are, ${\left(m_p/m_e\right)}^{10}\cong \exp \left(1/{\alpha }_s^2\right)$ and $2G_s{m}_e/c^2\cong \sqrt{G_F/\hslash c}.$