Formation of transient protein complexes is an important process in cells. Details of the association process as well as the energy landscapes of association are not well understood. In particular, the nature, height and position of the energy barriers during complexation are debated. Computational studies are well suited for atomistically investigating protein association processes. The Barnase-Barstar complex constitutes a well-studied target for computational studies as a small system with fast association rates. Here, we performed constraint biased Molecular Dynamics simulations along the reaction coordinate reaching from the diffusion regime to the bound complex. We simulated the wild-type and different mutants at different salt concentrations. A structural analysis of our simulation trajectories revealed not a single, but two distinct association patterns dominated by an interplay between two charged contact points near the binding site. Electrostatics and/or mutations influence the relative population of these patterns. Further, we computed the energy landscape of association as PMF (Potential of Mean Force) profiles within a reasonable agreement to experiment. We find a single energy barrier at a distance of similar to 0.3 nm, which corresponds to the final desolvation transition. Electrostatics has a profound influence on the height of this energy barrier, but not on its position. (C) 2010 Elsevier Inc. All rights reserved.