Ligand Binding Pathway Elucidation for Cryptophane Host-Guest Complexes

Systematic Dissociation Pathway Searches Guided by Principal Component Modes

We introduce a novel method, Pathway Search guided by Internal Motions (PSIM), that efficiently finds molecular dissociation pathways of a ligand-receptor system with guidance from principal component (PC) modes obtained from molecular dynamics (MD) simulations. Modeling ligand-receptor dissociation pathways can provide insights into molecular recognition and has practical applications, including understanding kinetic mechanisms and barriers to binding/unbinding as well as design of drugs with desired kinetic properties. PSIM uses PC modes in multilayer internal coordinates to identify natural molecular motions that guide the search for conformational switches and unbinding pathways. The new multilayer internal coordinates overcome problems with Cartesian and classical internal coordinates that fail to smoothly present dihedral rotation or generate nonphysical distortions. We used HIV-1 protease, which has large-scale flap motions, as an example protein to demonstrate use of the multilayer internal coordinates. We provide examples of algorithms and implementation of PSIM with alanine dipeptide and chemical host-guest systems, 2-naphthyl ethanol-β-cyclodextrin and tetramethylammonium-cryptophane complexes. Tetramethylammonium-cryptophane has slow binding/unbinding kinetics. Its residence time, the length to dissociate tetramethylammonium from the host, is ∼14 s from experiments, and PSIM revealed 4 dissociation pathways in approximately 150 CPU h. We also searched the releasing pathways for the product glyceraldehyde-3-phosphate from tryptophan synthase, and one complete dissociation pathway was constructed after running multiple search iterations in approximately 300 CPU h. With guidance by internal PC modes from MD simulations, the PSIM method has advantages over simulation-based methods to search for dissociation pathways of molecular systems with slow noncovalent kinetic behavior.

Xe affinities of water-soluble cryptophanes and the role of confined water

Simulations provide molecular insight on the aqueous binding of Xe to cryptophanes.

Given their relevance to drug design and chemical sensing, host–guest interactions are of broad interest in molecular science. Natural and synthetic host molecules provide vehicles for understanding selective molecular recognition in aqueous solution. Here, cryptophane–Xe host–guest systems are considered in aqueous media as a model molecular system that also has important applications. ¹²⁹Xe–cryptophane systems can be used in the creation of biosensors and powerful contrast agents for magnetic resonance imaging applications. Detailed molecular information on the determinants of Xe affinity is difficult to obtain experimentally. Thus, molecular simulation and free energy perturbation methods were applied to estimate the affinities of Xe for six water-soluble cryptophanes. The calculated affinities correlated well with the previously measured experimental values. The simulations provided molecular insight on the differences in affinities and the roles of conformational fluctuations, solvent, and counter ions on Xe binding to these host molecules. Displacement of confined water from the host interior cavity is a key component of the binding equilibrium, and the average number of water molecules within the host cavity is correlated with the free energy of Xe binding to the different cryptophanes. The findings highlight roles for molecular simulation and design in modulating the relative strengths of host–guest and host–solvent interactions.