ObjectiveThe coronavirus disease 2019 (COVID-19) pandemic has caused an exponential rise in death rates and hospitalizations. The aim of this study was to characterize the D614G substitution in the severe acute respiratory syndome coronavirus 2 (SARS-CoV-2) spike glycoprotein (S protein), which may affect viral infectivity.MethodsThe effect of D614G substitution on the structure and thermodynamic stability of the S protein was analyzed with use of DynaMut and SCooP. HDOCK and PRODIGY were used to model furin protease binding to the S protein RRAR cleavage site and calculate binding affinities. Molecular dynamics simulations were used to predict the S protein apo structure, the S protein–furin complex structure, and the free binding energy of the complex.ResultsThe D614G substitution in the G clade of SARS-CoV-2 strains introduced structural mobility and decreased the thermal stability of the S protein (ΔΔG = −0.086 kcal mol−1). The substitution resulted in stronger binding affinity (Kd = 1.6 × 10−8) for furin, which may enhance S protein cleavage. The results were corroborated by molecular dynamics simulations demonstrating higher binding energy of furin and the S protein D614G mutant (−61.9 kcal mol−1 compared with −56.78 kcal mol−1 for wild-type S protein).ConclusionsThe D614G substitution in the G clade induced flexibility of the S protein, resulting in increased furin binding, which may enhance S protein cleavage and infiltration of host cells. Therefore, the SARS-CoV-2 D614G substitution may result in a more virulent strain.