The viral oncoprotein E6 is related with the development of cervical cancer in HPV positive patients. E6 has been reported to bind to TNF-α, p53, p600 and p21 (Moody, C. et al., 2010), blocking the traditional pathways to apoptosis and leading to unregulated cell growth. This work aimed to use computational methods such as virtual screening (both target-based and ligand-based), molecular dynamics simulation, toxicity prediction and rational ligand design in order to search for potential inhibitors to Human Papillomavirus E6 oncoprotein. We also conducted a parallel study using other in silico methods such as binding site prediction and pharmacophore identification in order to thoroughly uncover the oncoprotein’s properties. The most common European variants of E6 oncoproteins, previously modeled by Tamarozzi & Giuliatti (2015), were submitted to target-based virtual screening by using GOLD and Autodock Vina software, and the molecules within the FDA-Approved database as ligand pool. The highest scoring protein/ligand complexes were submitted to molecular dynamics simulation using GROMACS, and then underwent visual analysis in both PyMOL and Chimera, in order to evaluate their molecular interactions with the active site of E6. Finally, we performed theoretical toxicity predictions and pharmacokinetics simulations through a series of open-source web servers, in order to select only the safest ligands, which were used as pivot molecules for the rational design of hybrids – novel molecules that combine the desired characteristics of the best ranked ligands, leading us to a considerable pool of potential inhibitors to Oncoprotein E6.