The Coulomb interaction between the two protons is fully included in the calculation of proton-deuteron breakup with realistic interactions for the first time. The hadron dynamics is based on the purely nucleonic charge-dependent (CD) Bonn potential and its realistic extension CD Bonn +Delta to a coupled-channel two-baryon potential, allowing for single virtual Delta-isobar excitation. Calculations are done using integral equations in momentum space. The screening and renormalization approach is employed for including the Coulomb interaction. The Coulomb effect on breakup observables is seen at all energies in particular kinematic regimes.