Molecular Dynamics (MD) is a method used to calculate the movement of atoms and molecules broadly applied to several aspects of science. It involves computational simulation, which makes it, at first glance, not easily accessible. The rise of several automated tools to perform molecular simulations has allowed researchers to navigate through the various steps of MD. This enables to elucidate structural properties of proteins that could not be analyzed otherwise, such as the impact of glycosylation. Glycosylation dictates the physicochemical and biological properties of a protein modulating its solubility, stability, resistance to proteolysis, interaction partners, enzymatic activity, binding and recognition. Given the high conformational and compositional diversity of the glycan chains, assessing their influence on the protein structure is challenging using conventional analytical techniques. In this manuscript, we present a step-by-step workflow to build and perform MD analysis of glycoproteins focusing on the SPIKE glycoprotein of SARS-CoV-2 to appraise the impact of glycans in structure stabilization and antibody occlusion.

Carboxyl-modified substrates are the most common chemical moieties that are frequently used as protein defibrillators. We studied the stability of protein-benzoic acid complexes with bovine serum albumin (BSA), zein and lysozyme proteins using various computational methods. Structural model for zein was built using homology modelling technique and molecular docking was used to prepare complex structures of all three proteins with benzoic acid. Molecular dynamics calculations performed on these complex structures provided a strong support for the stability of protein-benzoic acid complexes. The results from various analyses including root-mean-square deviation (RMSD) and radius of gyration showed the stability and compactness of all proteins-benzoic acid complexes. Moreover, exploration of structural fluctuations in proteins revealed the stability of active site residues. Two potential binding modes of benzoic acid with all three proteins were identified via cluster analysis. The binding mode which was retrieved from top cluster containing 86-91% of total conformations displayed very strong binding interactions for zein, BSA and lysozyme proteins. In addition, the results of binding mode showed that various interactions, including hydrogen binding, hydrophobic and electrostatic interactions were important for the optimal binding of benzoic acid with the active sites of proteins. Exploration of solvent accessible surface area showed that lysozyme-binding cavity was more exposed to the surface as compared to the other two proteins. Free energy analysis of all protein systems showed the stability of protein-benzoic acid complexes with lysozyme and BSA relatively more stable than zein system. The results of our study provided important insights to the dynamic and structural information about protein-benzoic acid interactions with BSA, zein and lysozyme proteins. This work is important in enhancing the stability of therapeutic protein drugs loaded on carboxyl substrates.

Meprins are complex and highly glycosylated multi-domain enzymes that require post-translational modifications to reach full activity. Meprins are metalloproteases of the astacin family characterized by a conserved zinc-binding motif (HExxHxxGFxHExxRxDR). Human meprin-α and -β protease subunits are 55% identical at the amino acid level, however the substrate and peptide bond specificities vary markedly. Current work focuses on the critical amino acid residues in the non-primed subsites of human meprins-α and -β involved in inhibitor/ligand binding. To compare the molecular events underlying ligand affinity, homology modeling of the protease domain of humep-α and -β based on the astacin crystal structure followed by energy minimization and molecular dynamics simulation of fully solvated proteases with inhibitor Pro-Leu-Gly-hydroxamate in S subsites were performed. The solvent accessible surface area curve shows a decrease in solvent accessibility values at specific residues upon inhibitor binding. The potential energy, total energy, H-bond interactions, root mean square deviation and root mean square fluctuation plot reflect the subtle differences in the S subsite of the enzymes which interact with different residues at P site of the inhibitor.