Preoperative diagnostic magnetic resonance (MR) breast images can provide good contrast between different tissues and 3-D information about suspicious tissues. Aligning preoperative diagnostic MR images with a patient in the theatre during breast conserving surgery could assist surgeons in achieving the complete excision of cancer with sufficient margins. Typically, preoperative diagnostic MR breast images of a patient are obtained in the prone position, while surgery is performed in the supine position. The significant shape change of breasts between these two positions due to gravity loading, external forces and related constraints makes the alignment task extremely difficult. Our previous studies have shown that either nonrigid intensity-based image registration or biomechanical modelling alone are limited in their ability to capture such a large deformation. To tackle this problem, we proposed in this paper a nonlinear biomechanical model-based image registration method with a simultaneous optimization procedure for both the material parameters of breast tissues and the direction of the gravitational force. First, finite element (FE) based biomechanical modelling is used to estimate a physically plausible deformation of the pectoral muscle and the major deformation of breast tissues due to gravity loading. Then, nonrigid intensity-based image registration is employed to recover the remaining deformation that FE analyses do not capture due to the simplifications and approximations of biomechanical models and the uncertainties of external forces and constraints. We assess the registration performance of the proposed method using the target registration error of skin fiducial markers and the Dice similarity coefficient (DSC) of fibroglandular tissues. The registration results on prone and supine MR image pairs are compared with those from two alternative nonrigid registration methods for five breasts. Overall, the proposed algorithm achieved the best registration p...