Diffusion-influenced ligand binding to buried sites in macromolecules and transmembrane channels

Solute translocation probability, lifetime, and "rectification" in membrane channels with localized constriction

We study the translocation probability and lifetime of a solute molecule in a cylindrical membrane channel that contains a localized constriction at an arbitrary location. Using a one-dimensional continuous diffusion description of solute dynamics in the channel, we explore two models. The first one describes a molecule's interaction with the constriction in terms of a narrow rectangular barrier in the potential of mean force. The second novel model proposed here represents this interaction by introducing an infinitely thin permeable partition. It is shown that when the parameters of the two models are chosen to warrant the same translocation probability, both models predict the same mean lifetime of the molecule in the channel. While the translocation probability is independent of the constriction location, the mean lifetime is a function of the location. The benefit of the thin partition model is that it allows one to lump together the height and length of the potential barrier into a single parameter, which is the partition's permeability. It is shown that in the case of an asymmetric location of the localized constriction and strong repulsion between the solutes, the solute flux through the channel is a function of the direction in which it goes, analogous to the phenomenon known in ion channel electrophysiology as rectification.

A multiscale predictive digital twin for neurocardiac modulation

Cardiac function is tightly regulated by the autonomic nervous system (ANS). Activation of the sympathetic nervous system increases cardiac output by increasing heart rate and stroke volume, while parasympathetic nerve stimulation instantly slows heart rate. Importantly, imbalance in autonomic control of the heart has been implicated in the development of arrhythmias and heart failure. Understanding of the mechanisms and effects of autonomic stimulation is a major challenge because synapses in different regions of the heart result in multiple changes to heart function. For example, nerve synapses on the sinoatrial node (SAN) impact pacemaking, while synapses on contractile cells alter contraction and arrhythmia vulnerability. Here, we present a multiscale neurocardiac modelling and simulator tool that predicts the effect of efferent stimulation of the sympathetic and parasympathetic branches of the ANS on the cardiac SAN and ventricular myocardium. The model includes a layered representation of the ANS and reproduces firing properties measured experimentally. Model parameters are derived from experiments and atomistic simulations. The model is a first prototype of a digital twin that is applied to make predictions across all system scales, from subcellular signalling to pacemaker frequency to tissue level responses. We predict conditions under which autonomic imbalance induces proarrhythmia and can be modified to prevent or inhibit arrhythmia. In summary, the multiscale model constitutes a predictive digital twin framework to test and guide high-throughput prediction of novel neuromodulatory therapy. KEY POINTS: A multi-layered model representation of the autonomic nervous system that includes sympathetic and parasympathetic branches, each with sparse random intralayer connectivity, synaptic dynamics and conductance based integrate-and-fire neurons generates firing patterns in close agreement with experiment. A key feature of the neurocardiac computational model is the connection between the autonomic nervous system and both pacemaker and contractile cells, where modification to pacemaker frequency drives initiation of electrical signals in the contractile cells. We utilized atomic-scale molecular dynamics simulations to predict the association and dissociation rates of noradrenaline with the β-adrenergic receptor. Multiscale predictions demonstrate how autonomic imbalance may increase proclivity to arrhythmias or be used to terminate arrhythmias. The model serves as a first step towards a digital twin for predicting neuromodulation to prevent or reduce disease.

Diffusion Resistance of Segmented Channels

The geometry of ion and metabolite channels in membranes of biological cells and organelles is usually far from that of a regular right cylinder. Rather, the channels have complex shapes that are characterized by the so-called vestibules and constriction zones which play roles of molecular filters determining the channel selectivity. In the present paper we discuss several channel structures with varying radius that approximate most of the cases found in nature, specifically, channels of smoothly varying radius and channels composed of multiple cylindrical sections of different lengths and radii including channels containing very thin circular constrictions. We consider diffusive transport of electrically neutral molecules driven by the concentration gradient and derive analytical expressions for the diffusion resistance - the integral parameter that describes steady-state transport properties of membrane channels.

Permeability and diffusion resistance of porous membranes: Analytical theory and its numerical test

This study is devoted to the transport of neutral solutes through porous flat membranes, driven by the solute concentration difference in the reservoirs separated by the membrane. Transport occurs through membrane channels, which are assumed to be non-overlapping, identical, straight cylindrical pores connecting the reservoirs. The key quantities characterizing transport are membrane permeability and its diffusion resistance. Such transport problems arising in very different contexts, ranging from plant physiology and cell biology to chemical engineering, have been studied for more than a century. Nevertheless, an expression giving the permeability for a membrane of arbitrary thickness at arbitrary surface densities of the channel openings is still unknown. Here, we fill in the gap and derive such an expression. Since this expression is approximate, we compare its predictions with the permeability obtained from Brownian dynamics simulations and find good agreement between the two.

Blocker Effect on Diffusion Resistance of a Membrane Channel: Dependence on the Blocker Geometry

Being motivated by recent progress in nanopore sensing, we develop a theory of the effect of large analytes, or blockers, trapped within the nanopore confines, on diffusion flow of small solutes. The focus is on the nanopore diffusion resistance which is the ratio of the solute concentration difference in the reservoirs connected by the nanopore to the solute flux driven by this difference. Analytical expressions for the diffusion resistance are derived for a cylindrically symmetric blocker whose axis coincides with the axis of a cylindrical nanopore in two limiting cases where the blocker radius changes either smoothly or abruptly. Comparison of our theoretical predictions with the results obtained from Brownian dynamics simulations shows good agreement between the two.

Effective diffusivity of a Brownian particle in a two-dimensional periodic channel of abruptly alternating width

We study diffusion of a Brownian particle in a two-dimensional periodic channel of abruptly alternating width. Our main result is a simple approximate analytical expression for the particle effective diffusivity, which shows how the diffusivity depends on the geometric parameters of the channel: lengths and widths of its wide and narrow segments. The result is obtained in two steps: first, we introduce an approximate one-dimensional description of particle diffusion in the channel, and second, we use this description to derive the expression for the effective diffusivity. While the reduction to the effective one-dimensional description is standard for systems of smoothly varying geometry, such a reduction in the case of abruptly changing geometry requires a new methodology used here, which is based on the boundary homogenization approach to the trapping problem. To test the accuracy of our analytical expression and thus establish the range of its applicability, we compare analytical predictions with the results obtained from Brownian dynamics simulations. The comparison shows excellent agreement between the two, on condition that the length of the wide segment of the channel is equal to or larger than its width.

Localized potential well vs binding site: Mapping solute dynamics in a membrane channel onto one-dimensional description

In the one-dimensional description, the interaction of a solute molecule with the channel wall is characterized by the potential of mean force U(x), where the x-coordinate is measured along the channel axis. When the molecule can reversibly bind to certain amino acid(s) of the protein forming the channel, this results in a localized well in the potential U(x). Alternatively, this binding can be modeled by introducing a discrete localized site, in addition to the continuum of states along x. Although both models may predict identical equilibrium distributions of the coordinate x, there is a fundamental difference between the two: in the first model, the molecule passing through the channel unavoidably visits the potential well, while in the latter, it may traverse the channel without being trapped at the discrete site. Here, we show that when the two models are parameterized to have the same thermodynamic properties, they automatically yield identical translocation probabilities and mean translocation times, yet they predict qualitatively different shapes of the translocation time distribution. Specifically, the potential well model yields a narrower distribution than the model with a discrete site, a difference that can be quantified by the distribution's coefficient of variation. This coefficient turns out to be always smaller than unity in the potential well model, whereas it may exceed unity when a discrete trapping site is present. Analysis of the translocation time distribution beyond its mean thus offers a way to differentiate between distinct translocation mechanisms.

Trapping of particles diffusing in two dimensions by a hidden binding site

We study trapping of particles diffusing in a two-dimensional rectangular chamber by a binding site located at the end of a rectangular sleeve. To reach the site a particle first has to enter the sleeve. After that it has two options: to come back to the chamber or to diffuse to the site where it is trapped. The main result of the present work is a simple expression for the mean particle lifetime as a function of its starting position and geometric parameters of the system. This expression is obtained by an approximate reduction of the initial two-dimensional problem to the effective one-dimensional one which can be solved with relative ease. Our analytical predictions are tested against the results obtained from Brownian dynamics simulations. The test shows excellent agreement between the two for a wide range of the geometric parameters of the system.

Detailed balance for diffusion in a potential with trapping and forward-backward symmetry of trapping time distributions

For particles diffusing in a potential, detailed balance guarantees the absence of net fluxes at equilibrium. Here, we show that the conventional detailed balance condition is a special case of a more general relation that works when the diffusion occurs in the presence of a distributed sink that eventually traps the particle. We use this relation to study the lifetime distribution of particles that start and are trapped at specified initial and final points. It turns out that when the sink strength at the initial point is nonzero, the initial and final points are interchangeable, i.e., the distribution is independent of which of the two points is initial and which is final. In other words, this conditional trapping time distribution possesses forward–backward symmetry.

Blocker escape kinetics from a membrane channel analyzed by mapping blocker diffusive dynamics onto a two-site model

When a large solute molecule enters a membrane channel from the membrane-bathing electrolyte solution, it blocks the small-ion current flowing through the channel. If the molecule spends in the channel sufficiently long time, individual blockades can be resolved in single-channel experiments. In this paper, we develop an analytical theory of the blocker escape kinetics from the channel, assuming that a charged blocking molecule cannot pass through a constriction region (bottleneck). We focus on the effect of the external voltage bias on the blocker survival probability in the channel. The bias creates a potential well for the charged blocker in the channel with the minimum located near the bottleneck. When the bias is strong, the well is deep, and escape from the channel is a slow process that allows for time-resolved observation of individual blocking events. Our analysis is performed in the framework of a two-site model of the blocker dynamics in the channel. Importantly, the rate constants, fully determining this model, are derived from a more realistic continuum diffusion model. This is done by mapping the latter onto its two-site counterpart which, while being much simpler, captures the main features of the blocker escape kinetics at high biases.

Mapping Intrachannel Diffusive Dynamics of Interacting Molecules onto a Two-Site Model: Crossover in Flux Concentration Dependence

This study focuses on how interactions of solute molecules affect the concentration dependence of their flux through narrow membrane channels. It is assumed that the molecules cannot bypass each other because of their hard-core repulsion. In addition, other short- and long-range solute-solute interactions are included into consideration. These interactions make it impossible to develop an analytical theory for the flux in the framework of a diffusion model of solute dynamics in the channel. To overcome this difficulty, we course-grain the diffusion model by mapping it onto a two-site one, where the rate constants describing the solute dynamics are expressed in terms of the parameters of the initial diffusion model. This allows us (i) to find an analytical solution for the flux as a function of the solute concentration and (ii) to characterize the solute-solute interactions by two dimensionless parameters. Such a characterization proves to be very informative as it results in a clear classification of the effects of the solute-solute interactions on the concentration dependence of the flux. Unexpectedly, this dependence can be nonmonotonic, exhibiting a sharp maximum in a certain parameter range. We hypothesize that such a behavior may constitute an element of a regulatory mechanism, wherein maximal flux reports on the optimal solute concentration in the bulk near the channel entrance.

Theory of Diffusion-Influenced Reaction Networks

A formalism is developed to describe how diffusion alters the kinetics of coupled reversible association-dissociation reactions in the presence of conformational changes that can modify the reactivity. The major difficulty in constructing a general theory is that, even to the lowest order, diffusion can change the structure of the rate equations of chemical kinetics by introducing new reaction channels (i.e., modifies the kinetic scheme). Therefore, the right formalism must be found that allows the influence of diffusion to be described in a concise and elegant way for networks of arbitrary complexity. Our key result is a set of non-Markovian rate equations involving stoichiometric matrices and net reaction rates (fluxes), in which these rates are coupled by a time-dependent pair association flux matrix, whose elements have a simple physical interpretation. Specifically, each element is the probability density that an isolated pair of reactants irreversibly associates at time t via one reaction channel on the condition that it started out with the dissociation products of another (or the same) channel. In the Markovian limit, the coupling of the chemical rates is described by committors (or splitting/capture probabilities). The committor is the probability that an isolated pair of reactants formed by dissociation at one site will irreversibly associate at another site rather than diffuse apart. We illustrate the use of our formalism by considering three reversible reaction schemes: (1) binding to a single site, (2) binding to two inequivalent sites, and (3) binding to a site whose reactivity fluctuates. In the first example, we recover the results published earlier, while in the second one we show that a new reaction channel appears, which directly connects the two bound states. The third example is particularly interesting because all species become coupled and an exchange-type bimolecular reaction appears. In the Markovian limit, some of the diffusion-modified rate constants that describe new transitions become negative, indicating that memory effects cannot be ignored.

Structural and dynamic basis of substrate permissiveness in hydroxycinnamoyltransferase (HCT)

Substrate permissiveness has long been regarded as the raw materials for the evolution of new enzymatic functions. In land plants, hydroxycinnamoyltransferase (HCT) is an essential enzyme of the phenylpropanoid metabolism. Although essential enzymes are normally associated with high substrate specificity, HCT can utilize a variety of non-native substrates. To examine the structural and dynamic basis of substrate permissiveness in this enzyme, we report the crystal structure of HCT from Selaginella moellendorffii and molecular dynamics (MD) simulations performed on five orthologous HCTs from several major lineages of land plants. Through altogether 17-μs MD simulations, we demonstrate the prevalent swing motion of an arginine handle on a submicrosecond timescale across all five HCTs, which plays a key role in native substrate recognition by these intrinsically promiscuous enzymes. Our simulations further reveal how a non-native substrate of HCT engages a binding site different from that of the native substrate and diffuses to reach the catalytic center and its co-substrate. By numerically solving the Smoluchowski equation, we show that the presence of such an alternative binding site, even when it is distant from the catalytic center, always increases the reaction rate of a given substrate. However, this increase is only significant for enzyme-substrate reactions heavily influenced by diffusion. In these cases, binding non-native substrates ‘off-center’ provides an effective rationale to develop substrate permissiveness while maintaining the native functions of promiscuous enzymes.

Examples abound of enzymes that can process substrates other than their native ones. However, the structural and dynamic basis of this promiscuity remains to be fully understood. In this work, we examine HCT, an intrinsically promiscuous acyltransferase with conserved function in all land plants. We uncover the sub-microsecond swing motion of a key arginine residue facilitating the recognition of both native and non-native substrates of HCT. We also quantify the impact of an off-center binding site on the non-native reaction rate. Although our calculations were inspired by HCT, the results apply in general, i.e., for enzymes heavily influenced by diffusion, binding non-native substrates ‘off-center’, even with rather weak affinity, can accelerate non-native reactions to appreciable levels.

Effect of surface curvature on diffusion-limited reactions on a curved surface

To investigate how the curvature of a reactive surface can affect reaction kinetics, we use a simple model in which a diffusion-limited bimolecular reaction occurs on a curved surface that is hollowed inward, flat, or extended outward while keeping the reactive area on the surface constant. By numerically solving the diffusion equation for this model using the finite element method, we find that the rate constant is a non-linear function of the surface curvature and that there is an optimal curvature providing the maximum value of the rate constant, which indicates that a spherical reactant whose entire surface is reactive (a uniformly reactive sphere) is not the most reactive species for a given reactive surface area. We discuss how this result arises from the interplay between two opposing effects: the exposedness of the reactive area to its partner reactants, which causes the rate constant to increase as the curvature increases, and the competition occurring on the reactive surface, which decreases the rate constant. This study helps us to understand the role of curvature in surface reactions and allows us to rationally design reactants that provide a high reaction rate.

Effect of stochastic gating on channel-facilitated transport of non-interacting and strongly repelling solutes

Ligand- or voltage-driven stochastic gating-the structural rearrangements by which the channel switches between its open and closed states-is a fundamental property of biological membrane channels. Gating underlies the channel's ability to respond to different stimuli and, therefore, to be functionally regulated by the changing environment. The accepted understanding of the gating effect on the solute flux through the channel is that the mean flux is the product of the flux through the open channel and the probability of finding the channel in the open state. Here, using a diffusion model of channel-facilitated transport, we show that this is true only when the gating is much slower than the dynamics of solute translocation through the channel. If this condition breaks, the mean flux could differ from this simple estimate by orders of magnitude.

Kinetics for Drug Discovery: an industry-driven effort to target drug residence time

A considerable number of approved drugs show non-equilibrium binding characteristics, emphasizing the potential role of drug residence times for in vivo efficacy. Therefore, a detailed understanding of the kinetics of association and dissociation of a target-ligand complex might provide crucial insight into the molecular mechanism-of-action of a compound. This deeper understanding will help to improve decision making in drug discovery, thus leading to a better selection of interesting compounds to be profiled further. In this review, we highlight the contributions of the Kinetics for Drug Discovery (K4DD) Consortium, which targets major open questions related to binding kinetics in an industry-driven public-private partnership.

Reversible Stochastically Gated Diffusion-Influenced Reactions

An approximate but accurate theory is developed for the kinetics of reversible binding of a ligand to a macromolecule when either can stochastically fluctuate between reactive and unreactive conformations. The theory is based on a set of reaction-diffusion equations for the deviations of the pair distributions from their bulk values. The concentrations are shown to satisfy non-Markovian rate equations with memory kernels that are obtained by solving an irreversible geminate (i.e., two-particle) problem. The relaxation to equilibrium is not exponential but rather a power law. In the Markovian limit, the theory reduces to a set of ordinary rate equations with renormalized rate constants.

Diffusive flux in a model of stochastically gated oxygen transport in insect respiration

Oxygen delivery to insect tissues is controlled by transport through a branched tubular network that is connected to the atmosphere by valve-like gates, known as spiracles. In certain physiological regimes, the spiracles appear to be randomly switching between open and closed states. Quantitative analysis of this regime leads a reaction-diffusion problem with stochastically switching boundary condition. We derive an expression for the diffusive flux at long times in this problem. Our approach starts with the derivation of the passage probability for a single particle that diffuses between a stochastically gated boundary, which models the opening and closing spiracle, and the perfectly absorbing boundary, which models oxygen absorption by the tissue. This passage probability is then used to derive an expression giving the diffusive flux as a function of the geometric parameters of the tube and characteristic time scales of diffusion and gate dynamics.

Predicting the influence of long-range molecular interactions on macroscopic-scale diffusion by homogenization of the Smoluchowski equation

The macroscopic diffusion constant for a charged diffuser is in part dependent on (1) the volume excluded by solute "obstacles" and (2) long-range interactions between those obstacles and the diffuser. Increasing excluded volume reduces transport of the diffuser, while long-range interactions can either increase or decrease diffusivity, depending on the nature of the potential. We previously demonstrated [P. M. Kekenes-Huskey et al., Biophys. J. 105, 2130 (2013)] using homogenization theory that the configuration of molecular-scale obstacles can both hinder diffusion and induce diffusional anisotropy for small ions. As the density of molecular obstacles increases, van der Waals (vdW) and electrostatic interactions between obstacle and a diffuser become significant and can strongly influence the latter's diffusivity, which was neglected in our original model. Here, we extend this methodology to include a fixed (time-independent) potential of mean force, through homogenization of the Smoluchowski equation. We consider the diffusion of ions in crowded, hydrophilic environments at physiological ionic strengths and find that electrostatic and vdW interactions can enhance or depress effective diffusion rates for attractive or repulsive forces, respectively. Additionally, we show that the observed diffusion rate may be reduced independent of non-specific electrostatic and vdW interactions by treating obstacles that exhibit specific binding interactions as "buffers" that absorb free diffusers. Finally, we demonstrate that effective diffusion rates are sensitive to distribution of surface charge on a globular protein, Troponin C, suggesting that the use of molecular structures with atomistic-scale resolution can account for electrostatic influences on substrate transport. This approach offers new insight into the influence of molecular-scale, long-range interactions on transport of charged species, particularly for diffusion-influenced signaling events occurring in crowded cellular environments.

Biased diffusion in three-dimensional comb-like structures

In this paper, we study biased diffusion of point Brownian particles in a three-dimensional comb-like structure formed by a main cylindrical tube with identical periodic cylindrical dead ends. It is assumed that the dead ends are thin cylinders whose radius is much smaller than both the radius of the main tube and the distance between neighboring dead ends. It is also assumed that in the main tube, the particle, in addition to its regular diffusion, moves with a uniform constant drift velocity. For such a system, we develop a formalism that allows us to derive analytical expressions for the Laplace transforms of the first two moments of the particle displacement along the main tube axis. Inverting these Laplace transforms numerically, one can find the time dependences of the two moments for arbitrary values of both the drift velocity and the dead-end length, including the limiting case of infinitely long dead ends, where the unbiased diffusion becomes anomalous at sufficiently long times. The expressions for the Laplace transforms are used to find the effective drift velocity and diffusivity of the particle as functions of its drift velocity in the main tube and the tube geometric parameters. As might be expected from common-sense arguments, the effective drift velocity monotonically decreases from the initial drift velocity to zero as the dead-end length increases from zero to infinity. The effective diffusivity is a more complex, non-monotonic function of the dead-end length. As this length increases from zero to infinity, the effective diffusivity first decreases, reaches a minimum, and then increases approaching a plateau value which is proportional to the square of the particle drift velocity in the main tube.

Diffusion-influenced reactions in a hollow nano-reactor with a circular hole

The problem of diffusion to a reactive sphere inside a cavity through a hole is solved for arbitrary aperture values.

From normal to anomalous diffusion in comb-like structures in three dimensions

Diffusion in a comb-like structure, formed by a main cylindrical tube with identical periodic dead ends of cylindrical shape, occurs slower than that in the same system without dead ends. The reason is that the particle, entering a dead end, interrupts its propagation along the tube axis. The slowdown becomes stronger and stronger as the dead end length increases, since the particle spends more and more time in the dead ends. In the limiting case of infinitely long dead ends, diffusion becomes anomalous with the exponent equal to 1/2. We develop a formalism which allows us to study the mean square displacement of the particle along the tube axis in such systems. The formalism is applicable for an arbitrary dead end length, including the case of anomalous diffusion in a tube with infinitely long dead ends. In particular, we demonstrate how intermediate anomalous diffusion arises when the dead ends are long enough.

Channel-facilitated diffusion boosted by particle binding at the channel entrance

We investigate single-file diffusion of Brownian particles in arrays of closely confining microchannels permeated by a variety of attractive optical potentials and connecting two baths with equal particle concentration. We simultaneously test free diffusion in the channel, diffusion in optical traps coupled in the center of the channel, and diffusion in traps extending into the baths. We found that both classes of attractive optical potentials enhance the translocation rate through the channel with respect to free diffusion. Surprisingly, for the latter class of potentials we measure a 40-fold enhancement in the translocation rate with respect to free diffusion and find a sublinear power law dependence of the translocation rate on the average number of particles in the channel. Our results reveal the function of particle binding at the channel entrances for diffusive transport and open the way to a better understanding of membrane transport and design of synthetic membranes with enhanced diffusion rate.

Microsecond folding experiments and simulations: a match is made

For the past two decades, protein folding experiments have been speeding up from the second or millisecond time scale to the microsecond time scale, and full-atom simulations have been extended from the nanosecond to the microsecond and even millisecond time scale. Where the two meet, it is now possible to compare results directly, allowing force fields to be validated and refined, and allowing experimental data to be interpreted in atomistic detail. In this perspective we compare recent experiments and simulations on the microsecond time scale, pointing out the progress that has been made in determining native structures from physics-based simulations, refining experiments and simulations to provide more quantitative underlying mechanisms, and tackling the problems of multiple reaction coordinates, downhill folding, and complex underlying structure of unfolded or misfolded states.

BDflex: a method for efficient treatment of molecular flexibility in calculating protein-ligand binding rate constants from brownian dynamics simulations

A method developed by Northrup et al. [J. Chem. Phys. 80, 1517 (1984)] for calculating protein-ligand binding rate constants (k(a)) from brownian dynamics (BD) simulations has been widely used for rigid molecules. Application to flexible molecules is limited by the formidable computational cost to treat conformational fluctuations during the long BD simulations necessary for k(a) calculation. Here, we propose a new method called BDflex for k(a) calculation that circumvents this problem. The basic idea is to separate the whole space into an outer region and an inner region, and formulate k(a) as the product of k(E) and η(d), which are obtained by separately solving exterior and interior problems. k(E) is the diffusion-controlled rate constant for the ligand in the outer region to reach the dividing surface between the outer and inner regions; in this exterior problem conformational fluctuations can be neglected. η(d) is the probability that the ligand, starting from the dividing surface, will react at the binding site rather than escape to infinity. The crucial step in reducing the determination of η(d) to a problem confined to the inner region is a radiation boundary condition imposed on the dividing surface; the reactivity on this boundary is proportional to k(E). By confining the ligand to the inner region and imposing the radiation boundary condition, we avoid multiple-crossing of the dividing surface before reaction at the binding site and hence dramatically cut down the total simulation time, making the treatment of conformational fluctuations affordable. BDflex is expected to have wide applications in problems where conformational fluctuations of the molecules are crucial for productive ligand binding, such as in cases where transient widening of a bottleneck allows the ligand to access the binding pocket, or the binding site is properly formed only after ligand entrance induces the closure of a lid.

Survival probability of an immobile target in a sea of evanescent diffusive or subdiffusive traps: a fractional equation approach

We calculate the survival probability of an immobile target surrounded by a sea of uncorrelated diffusive or subdiffusive evanescent traps (i.e., traps that disappear in the course of their motion). Our calculation is based on a fractional reaction-subdiffusion equation derived from a continuous time random walk model of the system. Contrary to an earlier method valid only in one dimension (d=1), the equation is applicable in any Euclidean dimension d and elucidates the interplay between anomalous subdiffusive transport, the irreversible evanescence reaction, and the dimension in which both the traps and the target are embedded. Explicit results for the survival probability of the target are obtained for a density ρ(t) of traps which decays (i) exponentially and (ii) as a power law. In the former case, the target has a finite asymptotic survival probability in all integer dimensions, whereas in the latter case there are several regimes where the values of the decay exponent for ρ(t) and the anomalous diffusion exponent of the traps determine whether or not the target has a chance of eternal survival in one, two, and three dimensions.

Finite Element Estimation of Protein-Ligand Association Rates with Post-Encounter Effects: Applications to Calcium binding in Troponin C and SERCA

We introduce a computational pipeline and suite of software tools for the approximation of diffusion-limited binding based on a recently developed theoretical framework. Our approach handles molecular geometries generated from high-resolution structural data and can account for active sites buried within the protein or behind gating mechanisms. Using tools from the FEniCS library and the APBS solver, we implement a numerical code for our method and study two Ca(2+)-binding proteins: Troponin C and the Sarcoplasmic Reticulum Ca(2+) ATPase (SERCA). We find that a combination of diffusional encounter and internal 'buried channel' descriptions provide superior descriptions of association rates, improving estimates by orders of magnitude.

Role of diffusion in the kinetics of reversible enzyme-catalyzed reactions

The accurate expression for the steady-state velocity of an irreversible enzyme-catalyzed reaction obtained by Shin and co-workers is generalized to allow for the rebinding of the product. The amplitude of the power-law (t(-1/2)) relaxation of the free- and bound-enzyme concentrations to steady-state values is expressed in terms of the steady-state velocity and the intrinsic (chemical) rate constants. This result is conjectured to be exact, even though our expression for the steady-state velocity in terms of microscopic parameters is only approximate.

Theory and simulation of diffusion-influenced, stochastically gated ligand binding to buried sites

We consider the diffusion-influenced rate coefficient of ligand binding to a site located in a deep pocket on a protein; the binding pocket is flexible and can reorganize in response to ligand entrance. We extend to this flexible protein-ligand system a formalism developed previously [A. M. Berezhkovskii, A, Szabo, and H.-X. Zhou, J. Chem. Phys. 135, 075103 (2011)] for breaking the ligand-binding problem into an exterior problem and an interior problem. Conformational fluctuations of a bottleneck or a lid and the binding site are modeled as stochastic gating. We present analytical and Brownian dynamics simulation results for the case of a cylindrical pocket containing a binding site at the bottom. Induced switch, whereby the conformation of the protein adapts to the incoming ligand, leads to considerable rate enhancement.