Summary: An analytic solution is presented for stresses induced in an infinite plate with two unequal circular holes by remote uniform loadings and arbitrary internal pressures in the holes. The solution is obtained by using the general expression for a biharmonic function in bipolar coordinates. The Airy stress function is decomposed in the sum of a fundamental stress function for an infinite plate remotely loaded, which gives non vanishing tractions on the circular boundaries, and an auxiliary stress function required to satisfy the boundary conditions on the pressures at the edges of the holes, which produces vanishing stresses at infinity. Correspondingly, the variations of the stress concentration factor are determined in terms of the holes geometry and loading conditions. The path independent \(J_{k}\)-(\(k=1,2\)), \(M\)- and \(L\)-integrals are analytically calculated on a closed contour encircling the two holes, under remote loading, in order to evaluate the energy release rates accompanying unit translation, self similar expansion and rotation of the holes, respectively. Results are then presented for varying loading orientation angle, biaxial loading ratio and holes geometry.