En los últimos años, la necesidad de reducir la dependencia del petróleo, ha promovido el uso de la biomasa vegetal no solo como fuente de energía, sino también de abundantes compuestos bioactivos, que poseen diversas propiedades organolépticas, físicoquímicas y biológicas de interés en el control de insectos como el mosquito Aedes aegypti, principal transmisor del virus Zika, la fiebre Chikungunya y el dengue. La posibilidad de contar con una metodología teórica simple para predecir la actividad larvicida de las moléculas, resulta una ventaja en el diseño racional de fitosanitarios basados en materiales renovables. Se emplea la relación cuantitativa estructura-actividad (QSAR), para predecir la concentración letal media (CL50) de 60 compuestos naturales, cuyos datos experimentales se obtienen de la literatura. La representación estructural se realiza a través del programa Chemsketch. Posteriormente, se calculan 15.254 descriptores moleculares mediante los programas de acceso libre PaDEL, Mold2 y EPI SuiteTM. La relación descriptores-propiedad se establece a través del ̳método de reemplazo‘, un novedoso algoritmo programado por nuestro grupo, que se basa en regresión lineal multivariable (MLR). La capacidad predictiva del modelo se define con un conjunto de moléculas de validación; además, se utilizan las técnicas de validación cruzada ̳dejar-uno-afuera‘ (loo) y aleatorización-Y (rand). Todos los algoritmos usados en este trabajo son programados en Matlab y se encuentran a disposición del público.