Protein ligand complex for aggregation study

Could anybody suggest how to perform MD simulations for protein ligand complex to study their aggregation??? Even suggest for Research articles are welcome