The identification and optimization of mutations in nanobodies are crucial for enhancing their therapeutic potential in disease prevention and control. However, this process is often complex and time-consuming, which limit its widespread application in practice. In this study, we developed a workflow, named Evolutionary-Nanobody (EvoNB), to predict key mutation sites of nanobodies by combining protein language models (PLMs) and molecular dynamic (MD) simulations. By fine-tuning the ESM2 model on a large-scale nanobody dataset, the ability of EvoNB to capture specific sequence features of nanobodies was significantly enhanced. The fine-tuned EvoNB model demonstrated higher predictive accuracy in the conserved framework and highly variable complementarity-determining regions of nanobodies. Additionally, we selected four widely representative nanobody-antigen complexes to verify the predicted effects of mutations. MD simulations analyzed the energy changes caused by these mutations to predict their impact on binding affinity to the targets. The results showed that multiple mutations screened by EvoNB significantly enhanced the binding affinity between nanobody and its target, further validating the potential of this workflow for designing and optimizing nanobody mutations. Additionally, sequence-based predictions are generally less dependent on structural absence, allowing them to be more easily integrated with tools for structural predictions, such as AlphaFold 3. Through mutation prediction and systematic analysis of key sites, we can quickly predict the most promising variants for experimental validation without relying on traditional evolutionary or selection processes. The EvoNB workflow provides an effective tool for the rapid optimization of nanobodies and facilitates the application of PLMs in the biomedical field.