This study aimed to investigate whether running kinematics can be accurately estimated through an artificial neural network (ANN) model containing GPS-based accelerometer variables and anthropometric data. Thirteen male participants with extensive running experience completed treadmill running trials at several speeds. Participants wore a GPS device containing a triaxial accelerometer, and running kinematics were captured by an 18-camera motion capture system for each trial. Multiple multilayer perceptron neural network models were constructed to estimate participants’ 3D running kinematics. The models consisted of the following input variables: 3D peak accelerometer acceleration during foot stance (g), stance time (s), running speed (km/h), participant height (cm), leg length (cm), and mass (kg). Pearson’s correlation coefficient (r), root mean squared error (RMSE), and relative root mean squared error (rRMSE) showed that ANN models provide accurate estimations of joint/segment angles (mean rRMSE = 13.0 ± 4.3%) and peak segment velocities (mean rRMSE = 22.1 ± 14.7%) at key gait phases across foot stance. The highest accuracies were achieved for flexion/extension angles of the thorax, pelvis, and hip, and peak thigh flexion/extension and vertical velocities (rRMSE < 10%). The current findings offer sports science and medical practitioners working with this data a method of conducting field-based analyses of running kinematics using a single IMU.