Human sweet taste receptor (hSTR) recognizes a wide array of sweeteners, resulting in sweet taste perception. Maltitol and lactitol have been extensively used in place of sucrose due to their capability to prevent dental caries. Herein, several molecular modeling approaches were applied to investigate the structural and energetic properties of these two polyols/hSTR complexes. Triplicate 500 ns molecular dynamics (MD) simulations and molecular mechanics/generalized Born surface area (MM/GBSA)-based free energy calculations revealed that the TAS1R2 monomer is the preferential binding site for maltitol and lactitol rather than the TAS1R3 region. Several polar residues (D142, S144, Y215, D278, E302, R383, and especially N143) were involved in polyols binding through electrostatic attractions and H-bond formations. The molecular complexation process not only induced the stable form of ligands but also stimulated the conformational adaptation of the TAS1R2 monomer to become a close-packed structure through an induced-fit mechanism. Notably, the binding affinity of the maltitol/TAS1R2 complex (ΔGbind of -17.93 ± 1.49 kcal/mol) was significantly higher than that of the lactitol/TAS1R2 system (-8.53 ± 1.78 kcal/mol), in line with the experimental relative sweetness. These findings provide an in-depth understanding of the differences in the sweetness response between maltitol and lactitol, which could be helpful to design novel polyol derivatives with higher sweet taste perception.
Global public health has been a critical problem by the sudden increase of the COVID-19 outbreak. The papain-like protease (PLpro) of SARS-CoV-2 is a key promising target for antiviral drug development since it plays a pivotal role in viral replication and innate immunity. Here, we employed the all-atom molecular dynamics (MD) simulations and binding free energy calculations based on MM-PB(GB)SA and SIE methods to elucidate and compare the binding behaviors of five inhibitors derived from peptidomimetic inhibitors (VIR250 and VIR251) and naphthalene-based inhibitors (GRL-0617, compound 3, and compound Y96) against SARS-CoV-2 PLpro. The obtained results revealed that all inhibitors interacting within the PLpro active site are mostly driven by vdW interactions, and the hydrogen bond formation in residues G163 and G271 with peptidomimetics and the Q269 residue with naphthalene-based inhibitors was essential for stabilizing the protein-ligand complexes. Among the five studied inhibitors, VIR250 exhibited the most binding efficiency with SARS-CoV-2 PLpro, and thus, it was chosen for the rational drug design. Based on the computationally designed ligand-protein complexes, the replacement of aromatic rings including heteroatoms (e.g., thiazolopyridine) at the P2 and P4 sites could help to improve the inhibitor-binding efficiency. Furthermore, the hydrophobic interactions with residues at P1-P3 sites can be increased by enlarging the nonpolar moieties (e.g., ethene) at the N-terminal of VIR250. We expect that the structural data obtained will contribute to the development of new PLpro inhibitors with more inhibitory potency for COVID-19.
The emergence outbreak caused by a novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has received significant attention on the global risks. Due to itscrucial role in viral replication, the main protease 3CLpro is an important target for drug discovery and development to combat COVID-19. In this work, the structural and dynamic behaviors as well as binding efficiency of the four peptidomimetic inhibitors (N3, 11a, 13b, and 14b) recently co-crystalized with SARS-CoV-2 3CLpro were studied and compared using all-atom molecular dynamics (MD) simulations and solvated interaction energy-based binding free energy calculations. The per-residue decomposition free energy results suggested that the key residues involved in inhibitors binding were H41, M49, L141-C145, H163-E166, P168, and Q189-T190 in the domains I and II. The van der Waals interaction yielded the main energy contribution stabilizing all the focused inhibitors. Besides, their hydrogen bond formations with F140, G143, C145, H164, E166, and Q189 residues in the substrate-binding pocket were also essential for strengthening the molecular complexation. The predicted binding affinity of the four peptidomimetic inhibitors agreed with the reported experimental data, and the 13b showed the most efficient binding to SARS-CoV-2 3CLpro. From rational drug design strategies based on 13b, the polar moieties (e.g., benzamide) and the bulky N-terminal protecting groups (e.g., thiazole) should be introduced to P1' and P4 sites in order to enhance H-bonds and hydrophobic interactions, respectively. We hope that the obtained structural and energetic information could be beneficial for developing novel SARS-CoV-2 3CLpro inhibitors with higher inhibitory potency to combat COVID-19.