Though significant efforts are in progress for developing drugs and vaccines against COVID-19, limited therapeutic agents are available currently. Thus, it is essential to undertake COVID-19 research and to identify therapeutic interventions in which computational modeling and virtual screening of lead molecules provide significant insights. The present study aimed to predict the interaction potential of natural lead molecules against prospective protein targets of SARS-CoV-2 by molecular modeling, docking, and dynamic simulation. Based on the literature survey and database search, fourteen molecular targets were selected and the three targets which lack the native structures were computationally modeled. The drug-likeliness and pharmacokinetic features of ninety-two natural molecules were predicted. Four lead molecules with ideal drug-likeliness and pharmacokinetic properties were selected and docked against fourteen targets, and their binding energies were compared with the binding energy of the interaction between Chloroquine and Hydroxychloroquine to their usual targets. The stabilities of selected docked complexes were confirmed by MD simulation and energy calculations. Four natural molecules demonstrated profound binding to most of the prioritized targets, especially, Hyoscyamine and Tamaridone to spike glycoprotein and Rotiorinol-C and Scutifoliamide-A to replicase polyprotein-1ab main protease of SARS-CoV-2 showed better binding energy, conformational and dynamic stabilities compared to the binding energy of Chloroquine and its usual target glutathione-S-transferase. The aforementioned lead molecules can be used to develop novel therapeutic agents towards the protein targets of SARS-CoV-2, and the study provides significant insight for structure-based drug development against COVID-19.