No Arabic abstract
Ligand-receptor binding and unbinding are fundamental biomolecular processes and particularly essential to drug efficacy. Environmental water fluctuations, however, impact the corresponding thermodynamics and kinetics and thereby challenge theoretical descriptions. Here, we devise a holistic, implicit-solvent, multi-method approach to predict the (un)binding kinetics for a generic ligand-pocket model. We use the variational implicit-solvent model (VISM) to calculate the solute-solvent interfacial structures and the corresponding free energies, and combine the VISM with the string method to obtain the minimum energy paths and transition states between the various metastable (dry and wet) hydration states. The resulting dry-wet transition rates are then used in a spatially-dependent multi-state continuous-time Markov chain Brownian dynamics simulations, and the related Fokker-Planck equation calculations, of the ligand stochastic motion, providing the mean first-passage times for binding and unbinding. We find the hydration transitions to significantly slow down the binding process, in semi-quantitative agreement with existing explicit-water simulations, but significantly accelerate the unbinding process. Moreover, our methods allow the characterization of non-equilibrium hydration states of pocket and ligand during the ligand movement, for which we find substantial memory and hysteresis effects for binding versus unbinding. Our study thus provides a significant step forward towards efficient, physics-based interpretation and predictions of the complex kinetics in realistic ligand-receptor systems.
Rebinding kinetics of molecular ligands plays a critical role in biomachinery, from regulatory networks to protein transcription, and is also a key factor for designing drugs and high-precision biosensors.In this study, we investigate initial release and rebinding of ligands to their binding sites grafted on a planar surface, a situation commonly observed in single molecule experiments and which occurs during exocytosis in vivo. Via scaling arguments and molecular dynamic simulations, we analyze the dependence of non-equilibrium rebinding kinetics on two intrinsic length scales: average separation distance between the binding sites and dimensions of diffusion volume (e.g., height of the experimental reservoir in which diffusion takes place or average distance between receptor-bearing surfaces). We obtain time-dependent scaling laws for on rates and for the cumulative number of rebinding events for various regimes. Our analyses reveal that, for diffusion-limited cases, the on rate decreases via multiple power law regimes prior to the terminal steady-state regime, in which the on rate becomes constant. At intermediate times, at which particle density has not yet become uniform throughout the reservoir, the number of rebindings exhibits a distinct plateau regime due to the three dimensional escape process of ligands from their binding sites. The duration of this regime depends on the average separation distance between binding sites. Following the three-dimensional diffusive escape process, a one-dimensional diffusive regime describes on rates. In the reaction-limited scenario, ligands with higher affinity to their binding sites delay the power laws. Our results can be useful for extracting hidden time scales in experiments where kinetic rates for ligand-receptor interactions are measured in microchannels, as well as for cell signaling via diffusing molecules.
We study unbinding of multivalent cationic ligands from oppositely charged polymeric binding sites sparsely grafted on a flat neutral substrate. Our molecular dynamics (MD) simulations are suggested by single-molecule studies of protein-DNA interactions. We consider univalent salt concentrations spanning roughly a thousandfold range, together with various concentrations of excess ligands in solution. To reveal the ionic effects on unbinding kinetics of spontaneous and facilitated dissociation mechanisms, we treat electrostatic interactions both at a Debye-H{u}ckel (DH, or `implicit ions, i.e., use of an electrostatic potential with a prescribed decay length) level, as well as by the more precise approach of considering all ionic species explicitly in the simulations. We find that the DH approach systematically overestimates unbinding rates, relative to the calculations where all ion pairs are present explicitly in solution, although many aspects of the two types of calculation are qualitatively similar. For facilitated dissociation (FD, acceleration of unbinding by free ligands in solution) explicit ion simulations lead to unbinding at lower free ligand concentrations. Our simulations predict a variety of FD regimes as a function of free ligand and ion concentrations; a particularly interesting regime is at intermediate concentrations of ligands where non-electrostatic binding strength controls FD. We conclude that explicit-ion electrostatic modeling is an essential component to quantitatively tackle problems in molecular ligand dissociation, including nucleic-acid-binding proteins.
Determining the complete set of ligands binding/unbinding pathways is important for drug discovery and to rationally interpret mutation data. Here we have developed a metadynamics-based technique that addressed this issue and allows estimating affinities in the presence of multiple escape pathways. Our approach is shown on a Lysozyme T4 variant in complex with the benzene molecule. The calculated binding free energy is in agreement with experimental data. Remarkably, not only we were able to find all the previously identified ligand binding pathways, but also we uncovered 3 new ones. This results were obtained at a small computational cost, making this approach valuable for practical applications, such as screening of small compounds libraries.
We investigated the yield stress and the apparent viscosity of sand with and without small amounts of liquid. By pushing the sand through a tube with an enforced Poiseuille like profile we minimize the effect of avalanches and shear localization. We find that the system starts to flow when a critical shear of the order of one particle diameter is exceeded. In contrast to common believe, we observe that the resistance against the flow of wet sand is much smaller than that of dry sand. For the dissipative flow we propose a non-equilibrium state equation for granular fluids.
Molecular interactions are key to many chemical and biological processes like protein function. In many signaling processes they occur in sub-cellular areas displaying nanoscale organizations and involving molecular assemblies. The nanometric dimensions and the dynamic nature of the interactions make their investigations complex in live cells. While super-resolution fluorescence microscopies offer live-cell molecular imaging with sub-wavelength resolutions, they lack specificity for distinguishing interacting molecule populations. Here we combine super-resolution microscopy and single-molecule Forster Resonance Energy Transfer (FRET) to identify dimers of receptors induced by ligand binding and provide super-resolved images of their membrane distribution in live cells. By developing a two-color universal-Point-Accumulation-In-the-Nanoscale-Topography (uPAINT) method, dimers of epidermal growth factor receptors (EGFR) activated by EGF are studied at ultra-high densities, revealing preferential cell-edge sub-localization. This methodology which is specifically devoted to the study of molecules in interaction, may find other applications in biological systems where understanding of molecular organization is crucial.