N-glycans (oligosaccharides linked to proteins via the side chain of asparagine) play an important role in protein folding control, cell adhesion, organism immune response, etc. The functional versatility comes from structural variety and the inherit flexibility of glycosidic linkages. Each N-glycan has several co-existing conformers posing a challenge for experimental characterisation and computational modelling. In this work, we assess the enhanced sampling capabilities of the replica-exchange MD (REMD) and metadynamics (MTD). While REMD has been used in the previous studies of N-glycans, it is the first application of MTD, as we know. Our proposed approach is based on the multi-replica MTD, where some degrees of freedom are biased explicit by the MTD procedure and the others implicitly by using replicas at higher temperatures. In addition, we apply the replica state exchange MTD (RSE-MTD), our proposed enhancement to increase replica exchange efficiency. Finally, the methods have been tested on Bi9, a biantennary complex-type N-glycan. We demonstrate that REMD, following a protocol from the previous studies, fails to provide consistent results, in contrast, our MTD approach enhances sampling along the critical collective variables and ensures converged results. The estimated conformer populations of Bi9 are in agreement with the experimental results obtained from the NMR pseudo-contact shift (PCS) measurements.
Read full abstract