Zinc ions are found at mM concentration in amyloid plaques of Alzheimer's disease and the role of zinc in protein oligomerization is the object of intense investigations. As an in vitro model for studying interactions between Zn(2+) and the Abeta peptide, that is the main component of plaques, the N- and C-termini protected Abeta(1-16) fragment has been chosen because reliable spectroscopic studies in water solution are possible due to the low propensity for oligomerization at pH approximately 6.5, and because all the Zn binding sites of Abeta have been identified in the 1-16 region. In this work we present the results of first principle simulations of several initial models of Zn-Abeta(1-16) complexes. The NMR results about the same system, where His 6, 13, 14 and Glu 11 side-chains coordinate the Zn ion, are strongly supported by these models. Coordination of Asp 1 to Zn drives the complex towards the expulsion of one of initially bonded His side-chains. Coordination of Tyr 10 to Zn is possible only when Tyr 10 is deprotonated. The interplay between physico-chemical properties of the Abeta ligand and the Zn coordination is discussed.