For hepatitis B virus (HBV)-related chronic infection under treatment by nucleos(t)ide analogues (NUCs), HBsAg clearance is the ultimate therapeutic goal but very infrequent. We investigated how HBV envelope protein variability could lead to differential HBsAg clearance on NUCs. For 12 HBV genotype D patients receiving NUCs, six resolvers (HBsAg clearance) were compared to six matched nonresolvers (HBsAg persistence). PreS/S amino acid (aa) sequences were analysed with bioinformatics to predict HBV envelope antigenicity and aa covariance. To enrich our analyses on very rare resolvers, these were compared with other HBV genotype D strains in three characterized clinical cohorts including common chronically infected patients. The sT125M+sP127T combination was observed in four nonresolvers of six, corroborated by aa covariance analysis, associated with a lower predicted antigenicity than sT125T+sP127P. Concordant features within this HBV key functional domain, at positions 125 and 127, were reported from two of the three comparative cohorts. In our hands, a lower ELISA reactivity of HBV-vaccinated mice sera was observed against the sT125M mutant. In the S gene, 56 aa changes in minor variants were detected in non-resolvers, mainly in the major hydrophilic region, vs 28 aa changes in resolvers. Molecular features in patients showing HBsAg persistence on NUCs argue in favour of a different aa pattern in the HBV S gene compared to those showing HBsAg clearance. In nonresolvers, a decrease in HBs ‘a’ determinant antigenicity and more frequent mutations in the S gene suggest a role for the HBV envelope characteristics in HBsAg persistence.