Neuroglobin (Ngb), a recently discovered member of the globin family, is overexpressed in the brain tissues over oxygen deprivation. Unlike more classical globins, such as myoglobin and hemoglobin, it is characterized by a hexacoordinated heme, and its physiological role is still unknown, despite the numerous investigations made oil the protein in recent years. Another important specific feature of human Ngb is the presence of two cysteine residues (Cys46 and Cys55), which are known to form an intramolecular disulfide bridge. Since previous work on human Ngb reported that its ligand binding properties could be controlled by the coordination state of the Fe2+ atom (in the heme moiety) and the redox state of the thiol groups, we choose to develop a Simulation approach combining coarse-grain Brownian dynamics and all-atom Molecular dynamics and metadynamics. We have Studied the diffusion of small ligands (CO, NO, and O-2) in the globin internal cavity network for Various states of human Ngb. Our results show how the structural and mechanical properties of the protein can be related to the ligand migration pathway, which can be extensively modified when changing the thiol's redox state and the iron's coordination state. We suggest that ligand binding is favored in the pentacoordinated species bearing an internal disulfide bridge.