[1] We have recently extended the global magnetohydrodynamic (MHD) model BATS-R-US to account for pressure anisotropy. Since the inner magnetosphere dynamics cannot be fully described even by anisotropic MHD, we coupled our anisotropic MHD model with two inner magnetospheric models: the Rice Convection Model (RCM) and the Comprehensive Ring Current Model (CRCM). The coupled models provide better representations of the near-Earth plasma, especially during geomagnetic storms. In this paper, we present the two-way coupling algorithms with both ring current models. The major difference between these two couplings is that the RCM assumes isotropic and constant pressures along closed field lines, while the CRCM resolves pitch angle anisotropy. For model validation, we report global magnetosphere simulations performed by the coupled models. The simulation results are compared to the results given by the coupled isotropic MHD and ring current models. We find that in the global MHD simulations coupled with ring current models, pressure anisotropy results in a thinner magnetosheath, a shorter tail, a much smaller Earthward plasma jet from the tail reconnection site, and is also important in controlling the magnetic field configuration. The comparisons with satellite data for the magnetospheric event simulations show improvements on reproducing the measured tail magnetic field and inner magnetospheric flow velocity when including pressure anisotropy in the ring current model coupled global MHD model.