Structural, dynamical, and vibrational properties, of complexes made of metal cobalt(III) coordinated to different amounts of cysteine molecules were investigated with DFT-based Car-Parrinello molecular dynamics (CPMD) simulations in liquid water solution. The systems are composed of Co(III):3Cys and Co(III):2Cys immersed in liquid water which are. modeled by about 110 explicit water molecules, thus one of the biggest molecular systems studied with ab initio molecular simulations so far. In such a way, we were able to investigate structural and dynamical properties of a model of a typical metal binding site used by several proteins. Cobalt, mainly a toxicological agent, can replace the natural binding metal and thus modify the biochemical activity. The structure of the surrounding solvent around the metal-ligands complexes is reported in detail, as well as the metal-ligands coordination bonds, using radial distribution functions and electronic analyses with Mayer bond orders. Structures of the Cocysteine complexes are found in very good agreement with EXAFS experimental data, stressing the importance of considering the surrounding solvent in the modeling. A vibrational. analysis is also conducted and compared. to experiment, which strengthens the reliability of the solvent interactions with the Cocysteine complexes from our molecular dynamics simulations, as well as the dynamics of the systems. From this preliminary analysis, we could suggest a vibrational fingerprint able to distinguish Co(III):2Cys from Co(III):3Cys. Our simulations also show the importance of considering a quantum explicit solvent, as solute-to-sol vent proton transfer events have been observed.