From binding data to pharmacokinetics: a novel approach to active drug absorbtion prediction

Oral administered drugs are mainly absorbed in the small intestine. Here, depending on drug composition and size, absorption can happen through a variety of processes . Through the epithelial cells and the lamina pro- pria the drug passes from the lumen into the blood stream in the capillaries. On its way it might be metabolised, transported away from the tract where absorption is possible or accumulate in organs other than those of treatment. Besides a fundamental interest in understanding the basic mechanisms by which a drug is assimilated by the human body, the kinetics of drug absorption is also a topic of much practical interest. A detailed knowledge of this process, resulting in the prediction of the drug absorption profile, can be of much help in the drug development stage .

To this end, several kinetic models for drug absorption within the body have been introduced (see e.g. ). They necessarily introduce some simplifications belonging to the category of the so-called three-compartment models where the substances (such as drugs or nutrients) move between three volumes (e.g. the human organs). In fact the models require two kinds of molecular properties. First are purely physical characteristics, such as solubility, differential solubility, LogP etc. These quantities are easy to measure or to calcualte, have direct physical meaning and sufficient to predict absorbtion profile of passively absorbed drugs. Actively transported molecules interact with protein transporters and therefore prediction for actively transporting compounds require a lot of separate knowledge of binding to and kinetics of the transporting proteins.





The major objective of this investigation was to develop a drug absorbtion prediction approach based on entirely different paradigm, thus avoiding difficulties of both knowledge-based and QSAR-based models, and therefore capable of better predictions. Recently it was observed that experimental values of molecular activities against a large proteins set can be used for predicting broad biological effects . In this investigation we take advantage of this concept and develop a novel quntitative method for identification of actively transported drugs. To do that we performed a docking study of a few hundreds small molecules (mostly drugs) against a diversified 510 proteins set representing human proteom. Using available absorbtion data for each of the molecules we obtained a support vector classifier capable to identify proteins which affinity for drugs correlates well with the active absorption of these drugs in 81% cases. The observation helped us improve our passive absorbtion model by adding non-liner fluxes associated with the transporting protein to obtain also a quantitative model of the passively absorbed drugs.

Ref: arXiv:0810.2617 [ps, pdf, other]

Title: From protein binding to pharmacokinetics: a novel approach to active drug absorption prediction

Authors: P.O. Fedichev, T.V. Kolesnikova, A.A. Vinnik

Comments: 9 pages, 5 eps figures

Subjects: Quantitative Methods (q-bio.QM); Biomolecules (q-bio.BM)

The nature of percolation phase transition in films of hydration water around immersed bodies.

In a set of molecular dynamics calculations (MD) the percolation phase transition in water layer absorbed on a body surface was revealed at definite temperature. Below this temperature the infinite network of unbroken hydrogen bonds exists. Above it this network decays on islands. This conclusion corresponds also with measurements of conduction of moisture contained disperse materials as quartz, for example: the conductivity drops almost to zero value while heating the specimens up to definite temperature. It is known that the water conductance dominates by the “estafette” mechanism in which protons are transferred over the hydrogen bonds. The breakdown of network means the conductivity drop. These phenomena are explained in the paper in frames of early published continuous vector model of polar liquids. It is shown that the immersed bodies are surrounded by the ferroelectric film, in which the dipole moments of water molecules are ordered, arranged in one direction parallel to the interface. It is the physics behind above mentioned MD results. In addition of our previous papers the stability of this ferroelectric order is proved. The character of phase transition to the paraelectric phase is discussed and its temperature is estimated that is in agreement with MD results. Below the critical temperature the polarization vector field contains the structures as “vortex-antivortex pairs”. These pairs dissociate above this temperature that means the order breaking. The boundary conditions for the polarization vector field of molecular dipole moments are derived that is necessary to enclose the vector model equations.

Reference: accepted for publication to Journal of Structual Chemistry (Russian Journal of), 2008

Spontaneous polarization of a polar liquid next to nano-scale impurities

Numerous properties of water are determined by the hydrogen bonds between its molecules. Water does not form hydrogen bonds with hydrophobic materials, henceforth, dipole moments of its molecules are arranged mainly parallel to the interfaces with such substances. According to molecular dynamics calculations (MD) at such orientation molecules save the maximal number of hydrogen bonds: three of fourth. It is shown in this Letter that in the layer of water or ice next to surface the long-range order spontaneously forms: remaining parallel to the surface dipole moment vectors arrange in one direction. Some fraction of dipole moments form the vortex structures on the surface. At low temperatures the ordered state has small admixture of vortex-antivortex pairs. The interaction energy of vortexes in this pairs arises proportional to the distance between them. A definite temperature the phase transition takes place: pairs suffer the dissociation, the molecular dipole moments order disappears. This conclusion agrees with he results of MD calculations, in which the percolation phase transition was revealed in the hydrogen bond network of water molecules absorbed on a surface.

The spontaneous polarization of liquid induced by the immersed in it nano-size bodies (proteins, peptides, …) results in the additional long-range interaction between them that depends on their relative orientation. Polarization of liquid in this case looks like that presented in Fig.1 in agreement with MD. All mentioned MD results can not be explained in frames of standard continuous scalar theory of water. These phenomena were analyzed here in frames of continuous vector model of polar liquids applications of which looks like promising to speed the simulations of macromolecular complexes.

Reference: arXiv:cond-mat/0601129 [ps, pdf, other]

Title: Long-Range Order and Interactions of Macroscopic Objects in Polar Liquids

Authors: P.O. Fedichev, L.I. Men'shikov

Comments: 11 pages, 6 figures

Subjects: Soft Condensed Matter (cond-mat.soft); Chemical Physics (physics.chem-ph); Biomolecules (q-bio.BM)

Accepted for publication in Journal of Physical Chemistry A (Russian Journal of), 2009

What's an ultimate value of reversible drug binding constant?

Traditional opinion is that a good drug must have a high value of the absolute meaning of the binding energy with target protein in order to prevent the thermal dissociation of the drug-protein complex. In this case an essential deformation of protein arises, which has to be taken into account in developing different models of protein-small molecule and protein-protein interaction, and computing affinity constants in drug discovery in-silico methods. The effect of essential perturbation of protein molecule is ignored in standard computational methods of drug design that can contribute a large mistake to results of calculation, to binding energy, for example.
To demonstrate the existence of the ultimate value of the binding energy two models are considered: macroscopic and microscopic, both giving the same conclusions: the critical value of absolute meaning of binding energy is 50-100kJ/M. If the binding energy exceeds this value, then drug essentially perturbs protein configuration. In a microscopic picture this perturbation is a sequence of irreversible conformational transitions in protein body. In a macroscopic one it is an inelastic deformation of a protein substance. Our estimation agrees with the experimental value (50 kJ /M) of the ultimate energy that can be stored in a protein molecule without its destruction.
The existence of the critical value of binding energy should be accounted in structure based drug design methods where protein molecule is considered in an elastic deformation approximation.

Reference: accepted in Russian Journal of Biophysics, 2008

Ferro-electric phase transition in a polar liquid and the nature of lambda-transition in supercooled water

Water is a major and all-important example of a strongly interacting polar liquid. Dielectric properties of water surrounding nano-scale objects pose a fundamentally important problem in physics, chemistry, structural biology and in silica drug design. The issue of temperature dependence of dielectric constant, the role of fluctuations and a possibility of a ferro-electric phase transition in a polar liquid is fairly old . It attracted a new attention when a new phase transition (so called lambda-transition) was observed in supercooled water at critical temperatures between T_{c}=228K and T_{c}=231.4K . Isothermal compressibility, density, diffusion coefficient, viscosity and static dielectric constant \epsilon and other quantities diverge as T_{c} is approached, which is signature of a second order phase transition. The singularity of \epsilon is a feature of a ferro-electric transition . However, given a complexity of interactions between water molecules, the physical picture behind this phenomenon is not entirely clear . In the phase transition is explained as a formation of a rigid network of hydrogen bonds. On the other hand a ferro-electric hypothesis was also proposed and supported by molecular-dynamics simulations (MD). For example, a ferro-electric liquid phase was observed in a model of the so called ``soft spheres'' with static dipole moments . In fact, the existence of a ferro-electric phase appears to be model independent: domains where formed both in MD calculations with hard spheres with point dipoles and with soft spheres with extended dipoles .

In the our latest publication, Ferro-electric phase transition in a polar liquid and the nature of lambda-transition in supercooled water, we develop two related approaches to calculate free energy of a polar liquid. We show that long range nature of dipole interactions between the molecules leads to para-electric state instability at sufficiently low temperatures and to a second-order phase transition. We establish the transition temperature, T_{c}, both within mean field and ring diagrams approximation and demonstrate that the ferro-electric transition is a sound physical explanation behind the experimentally observed \lambda-transition in supercooled water. Finally we discuss dielectric properties, the role of fluctuations and establish connections with earlier phenomenological models of polar liquids.

Reference: arXiv:0808.0991 [ps, pdf, other]

Title: Ferro-electric phase transition in a polar liquid and the nature of \lambda-transition in supercooled water

Authors: P.O. Fedichev, L.I. Menshikov

Comments: 4 pages, 1 eps figure

Subjects: Statistical Mechanics (cond-mat.stat-mech); Soft Condensed Matter (cond-mat.soft)