Predicting the 3D structure of protein interactions remains a challenge in the field of computational structural biology. This is in part due to difficulties in sampling the complex energy landscape of multiple interacting flexible polypeptide chains. Coarse-graining approaches, which reduce the number of degrees of freedom of the system, help address this limitation by smoothing the energy landscape, allowing an easier identification of the global energy minimum. They also accelerate the calculations, allowing to model larger assemblies. Here, we present the implementation of the MARTINI coarse-grained force field for proteins into HADDOCK, our integrative modelling platform. Docking and refinement are performed at the coarse-grained level and the resulting models are then converted back to atomistic resolution through a distance restraints-guided morphing procedure. Our protocol, tested on the largest complexes of the protein docking benchmark 5, shows an overall ~7-fold speed increase compared to standard all-atom calculations, while maintaining a similar accuracy and yielding substantially more near-native solutions. To showcase the potential of our method, we performed simultaneous 7 body docking to model the 1:6 KaiC-KaiB complex, integrating mutagenesis and hydrogen/deuterium exchange data from mass spectrometry with symmetry restraints, and validated the resulting models against a recently published cryo-EM structure.between the number of estimated protein-protein interactions and those deposited in the Protein Data Bank 4 can be overcome based solely on experimental methods 5 .Computational docking has come of age as a complement to experimental methods in order to generate 3D models of protein assemblies. In particular, data-or information-driven docking and other integrative approaches are particularly appealing 1,6-8 . While docking performs sufficiently well for small-and medium-sized proteins, applications to large biological systems, either containing large individual molecules or a large number of interactors, are limited by the significant computational cost of thoroughly sampling complex conformational landscapes.Coarse-grained (CG) models mitigate this limitation by grouping atoms into larger "pseudoatoms" or beads 9-11 , thus reducing the number of particles to consider in the computations. These models were used in the very first energy minimization of a protein in 1969 12 and again in the first docking simulation 13 .Since then, several CG models have been developed and applied to study different aspects of protein structural biology 14 . For protein docking in particular, of the CG models developed over the years, three standouts for their performance and/or success in community assessment experiments: Those implemented in ATTRACT, CABS-dock, and RosettaDock. The ATTRACT model 15,16 , developed by Zacharias and coworkers for flexible protein docking, represents the protein backbone by two pseudo-atoms and the side chains by an additional particle (or two in the case of larger amino acids). ...