Ritonavir and xk263 Binding-Unbinding with HIV-1 Protease: Pathways, Energy and Comparison

Understanding non-covalent biomolecular recognition, which includes drug–protein bound states and their binding/unbinding processes, is of fundamental importance in chemistry, biology, and medicine. Fully revealing the factors that govern the binding/unbinding processes can further assist in designing drugs with desired binding kinetics. HIV protease (HIVp) plays an integral role in the HIV life cycle, so it is a prime target for drug therapy. HIVp has flexible flaps, and the binding pocket can be accessible by a ligand via various pathways. Comparing ligand association and dissociation pathways can help elucidate the ligand–protein interactions such as key residues directly involved in the interaction or specific protein conformations that determine the binding of a ligand under certain pathway(s). Here, we investigated the ligand unbinding process for a slow binder, ritonavir, and a fast binder, xk263, by using unbiased all-atom accelerated molecular dynamics (aMD) simulation with a re-seeding approach and an explicit solvent model. Using ritonavir-HIVp and xk263-HIVp ligand–protein systems as cases, we sampled multiple unbinding pathways for each ligand and observed that the two ligands preferred the same unbinding route. However, ritonavir required a greater HIVp motion to dissociate as compared with xk263, which can leave the binding pocket with little conformational change of HIVp. We also observed that ritonavir unbinding pathways involved residues which are associated with drug resistance and are distal from catalytic site. Analyzing HIVp conformations sampled during both ligand–protein binding and unbinding processes revealed significantly more overlapping HIVp conformations for ritonavir-HIVp rather than xk263-HIVp. However, many HIVp conformations are unique in xk263-HIVp unbinding processes. The findings are consistent with previous findings that xk263 prefers an induced-fit model for binding and unbinding, whereas ritonavir favors a conformation selection model. This study deepens our understanding of the dynamic process of ligand unbinding and provides insights into ligand–protein recognition mechanisms and drug discovery.

GENERIC-compliant simulations of Brownian multi-particle systems: modeling stochastic lubrication

A stochastic Lagrangian model for simulating the dynamics and rheology of a Brownian multi-particle system interacting with a simple liquid medium is presented. The discrete particle model is formulated within the GENERIC framework for Non-Equilibrium Thermodynamics and therefore it satisfies discretely the First/Second Laws of Thermodynamics and the Fluctuation Dissipation Theorem (FDT). Long-range fluctuating hydrodynamics interactions between suspended particles are described by an explicit solvent model. To this purpose, the Smoothed Dissipative Particle Dynamics method is adopted, which is a GENERIC-compliant Lagrangian meshless discretization of the fluctuating Navier–Stokes equations. In dense multi-particle systems, the average inter-particle distance is typically small compared to the particle size and short-range hydrodynamics interactions play a major role. In order to bypass an explicit—computationally costly—solution for these forces, a lubrication correction is introduced based on semi-analytical expressions for spheres under Stokes flow conditions. We generalize here the lubrication formalism to Brownian conditions, where an additional thermal-lubrication contribution needs to be taken into account in a way that discretely satisfies FDT. The coupled lubrication dynamics is integrated in time using a generalized semi-implicit splitting scheme for stochastic differential equations. The model is finally validated for a single particle diffusion as well as for a Brownian multi-particle system under homogeneous shear flow. Results for the diffusional properties as well as the rheological behavior of the whole suspension are presented and discussed.

Hydrogen bonding interactions between arsenious acid and dithiothreitol/dithioerythritol at different pH values: A computational study under explicit solvent model

Sulfhydryl compounds are regarded as potential functional monomers for arsenious acid imprinted polymers due to their high affinity to arsenious acid in aqueous solution. However, the recognition and binding mechanism...

Transport of polymer-modified nanoparticles in nanochannels coated with polymers

Using molecular dynamics simulations based on explicit-solvent model, we study migration of polymer-modified nanoparticles through nanochannels coated with polymers.

Unexpected solvent effects on the UV/Vis absorption spectra of o -cresol in toluene and benzene: in contrast with non-aromatic solvents

Cresol is a prototype molecule in understanding intermolecular interactions in material and biological systems, because it offers different binding sites with various solvents and protonation states under different pH values. It is found that the UV/Vis absorption spectra of o -cresol in aromatic solvents (benzene, toluene) are characterized by a sharp peak, unlike the broad double-peaks in 11 non-aromatic solvents. Both molecular dynamics simulations and electronic structure calculations revealed the formation of intermolecular π-complexation between o -cresol and aromatic solvents. The thermal movements of solvent and solute molecules render the conformations of o -cresol changing between trans and cis isomers. The π-interaction makes the cis configuration a dominant isomer, hence leading to the single keen-edged UV/Vis absorption peak at approximately 283 nm. The free conformation changes between trans and cis in aqueous solution rationalize the broader absorption peaks in the range of 260–280 nm. The pH dependence of the UV/Vis absorption spectra in aqueous solutions is also rationalized by different protonation states of o -cresol. The explicit solvent model with long-ranged interactions is vital to describe the effects of π-complexation and electrostatic interaction on the UV/Vis absorption spectra of o -cresol in toluene and alkaline aqueous (pH > 10.3) solutions, respectively.