Synthesis and unexpected binding of monofluorinated N,N'-diacetylchitobiose and LacdiNAc to wheat germ agglutinin

Fluorination of carbohydrate ligands of lectins is a useful approach to examine their binding profile, improve their metabolic stability and lipophilicity, and convert them into 19F NMR-active probes. However, monofluorination of monovalent carbohydrate ligands often leads to a decreased or completely lost affinity. By chemical glycosylation, we synthesized the full series of methyl β-glycosides of N,N'-diacetylchitobiose (GlcNAcβ(1-4)GlcNAcβ1-OMe) and LacdiNAc (GalNAcβ(1-4)GlcNAcβ1-OMe) systematically monofluorinated at all hydroxyl positions. A competitive enzyme-linked lectin assay revealed that the fluorination at the 6'-position of chitobioside resulted in an unprecedented increase in affinity to wheat germ agglutinin (WGA) by one order of magnitude. For the first time, we have characterized the binding profile of a previously underexplored WGA ligand LacdiNAc. Surprisingly, 4'-fluoro-LacdiNAc bound WGA even stronger than unmodified LacdiNAc. These observations were interpreted using molecular dynamic calculations along with STD and transferred NOESY NMR techniques, which gave evidence for the strengthening of CH/π interactions after deoxyfluorination of the side chain of the non-reducing GlcNAc. These results highlight the potential of fluorinated glycomimetics as high-affinity ligands of lectins and 19F NMR-active probes.

Binding thermodynamics and interaction patterns of human purine nucleoside phosphorylase-inhibitor complexes from extensive free energy calculations

Human purine nucleoside phosphorylase (hPNP) plays a significant role in the catabolism of deoxyguanosine. The trimeric protein is an important target in the treatment of T-cell cancers and autoimmune disorders. Experimental studies on the inhibition of the hPNP observe that the first ligand bound to one of three subunits effectively inhibits the protein, while the binding of more ligands to the subsequent sites shows negative cooperativities. In this work, we performed extensive end-point and alchemical free energy calculations to determine the binding thermodynamics of the trimeric protein-ligand system. 13 Immucillin inhibitors with experimental results are under calculation. Two widely accepted charge schemes for small molecules including AM1-BCC and RESP are adopted for ligands. The results of RESP are in better agreement with the experimental reference. Further investigations of the interaction networks in the protein-ligand complexes reveal that several residues play significant roles in stabilizing the complex structure. The most commonly observed ones include PHE200, GLU201, MET219, and ASN243. The conformations of the protein in different protein-ligand complexes are observed to be similar. We expect these insights to aid the development of potent drugs targeting hPNP.

SAMPL7 TrimerTrip host-guest binding affinities from extensive alchemical and end-point free energy calculations

The prediction of host-guest binding affinities with computational modelling is still a challenging task. In the 7th statistical assessment of the modeling of proteins and ligands (SAMPL) challenge, a new host named TrimerTrip was synthesized and the thermodynamic parameters of 16 structurally diverse guests binding to the host were characterized. In the TrimerTrip-guest challenge, only structures of the host and the guests are provided, which indicates that the predictions of both the binding poses and the binding affinities are under assessment. In this work, starting from the binding poses obtained from our previous enhanced sampling simulations in the configurational space, we perform extensive alchemical and end-point free energy calculations to calculate the host-guest binding affinities retrospectively. The alchemical predictions with two widely accepted charge schemes (i.e. AM1-BCC and RESP) are in good agreement with the experimental reference, while the end-point estimates perform poorly in reproducing the experimental binding affinities. Aside from the absolute value of the binding affinity, the rank of binding free energies is also crucial in drug design. Surprisingly, the end-point MM/PBSA method seems very powerful in reproducing the experimental rank of binding affinities. Although the length of our simulations is long and the intermediate spacing is dense, the convergence behavior is not very good, which may arise from the flexibility of the host molecule. Enhanced sampling techniques in the configurational space may be required to obtain fully converged sampling. Further, as the length of sampling in alchemical free energy calculations already achieves several hundred ns, performing direct simulations of the binding/unbinding event in the physical space could be more useful and insightful. More details about the binding pathway and mechanism could be obtained in this way. The nonequilibrium method could also be a nice choice if one insists to use the alchemical method, as the intermediate sampling is avoided to some extent.