Fluctuations in the bending angles at internal irregularities of DNA and RNA (such as symmetric loops, bulges, and nicks/gaps) have been observed from various experiments. However, little effort has been made to computationally predict and explain the statistical behavior of semiflexible chains with internal defects. In this paper, we describe the general structure of these macromolecular chains as inextensible elastic chains with one or more internal joints which have limited ranges of rotation, and propose a method to compute the probability density functions of the end-to-end pose of these macromolecular chains. Our method takes advantage of the operational properties of the non-commutative Fourier transform for the group of rigid-body motions in three-dimensional space, SE(3). Two representative types of joints, the hinge for planar rotation and the ball joint for spatial rotation, are discussed in detail. The proposed method applies to various stiffness models of semi-flexible chain-like macromolecules. Examples are calculated using the Kratky-Porod model with specified stiffness, angular fluctuation, and joint locations. Entropic effects associated with internal angular fluctuations of semi-flexible macromolecular chains with internal joints can be computed using this formulation. Our method also provides a potential tool to detect the existence of internal irregularities.