The reported work aims to improve the performance of LSTM-based (Long Short-Term Memory) forecasting algorithms in cases of NARX (Nonlinear Autoregressive with eXogenous input) models by using evolutionary search. The proposed approach, ES-LSTM, combines a two-membered ES local search procedure (2MES) with an ADAM optimizer to train more accurate LSTMs. The accuracy is measured from both error and trend prediction points of view. The method first computes the learnable parameters of an LSTM, using a subset of the training data, and applies a modified version of 2MES optimization to tune them. In the second stage, all available training data are used to update the LSTM’s weight parameters. The performance of the resulting algorithm is assessed versus the accuracy of a standard trained LSTM in the case of multiple financial time series. The tests are conducted on both training and test data, respectively. The experimental results show a significant improvement in the forecasting of the direction of change without damaging the error measurements. All quality measures are better than in the case of the standard algorithm, while error measures are insignificantly higher or, in some cases, even better. Together with theoretical consideration, this proves that the new method outperforms the standard one.