The paper presents modeling of steels strength characteristics in dependence from their alloying components quantities using neural networks as nonlinear approximation functions. Further, for optimization purpose the neural network models are used. The gradient descent algorithm based on utility function backpropagation through the models is applied. The approach is aimed at synthesis of steel alloys compositions with improved strength characteristics by solving multi-criteria optimization task. The obtained optimal alloying compositions fall into martenzite region of steels. They will be subject of further experimental testing in order to synthesize new steels with desired characteristics.