Skin aging is the most intuitive manifestation of aging. Skin aging inevitably leads to cosmetic and psychological problems, and even diseases. The present study aims to research the pathological and molecular mechanisms underlying skin aging and identify the therapeutic agents for reversing skin aging. Two available gene expression datasets (GSE55118 and GSE72264) for skin aging were downloaded from Gene Expression Omnibus, followed by bioinformatic analyses performed on the datasets. Firstly, 169 crucial mRNAs, 27 crucial miRNAs and 50 crucial lncRNAs closely related to skin aging were identified by weighted gene co-expression network analysis. Then, function Enrichment Analysis performed by Metascape database showed that skin aging involves a variety of biological functions, such as detection of stimulus, response to steroid hormone and water channel activity, regulation of muscle contraction. Next, ten hub genes including AQP4, TRPM8, TBR1, NTSR2, MPPED1, BARHL2, PAX9, CPN1, CES3, and CHGB were screened out by the protein-protein interaction analysis. Next, the "lncRNA-miRNA-mRNA" network and the "lncRNA-miRNA-hub mRNA" network were constructed to explore the competing endogenous RNAs mechanism of skin aging. Finally, ten significant potential small molecules mitigating skin aging were screened using CMAP platform, including tretinoin, pifithrin, selamectin, entinostat, bretazenil, syringic-acid, BRD-K96475865, emedastine, abacavir, and rotenone, and their reliability was verified by molecular docking experiments. The present study provided basis for revealing the molecular mechanism of skin aging and identified the potential candidate drugs for mitigating skin aging.