Screening for transient biological interactions as applied to albumin ligands: a new concept for drug discovery

The notion that many biological interactions are based on transient binding (dissociation constants (K(d)) in the range of 10-0.01 mM) is familiar, yet the implications for biological sciences have been realized only recently. An important area of biological sciences is drug design, where the traditional "lock and key" view of binding has prevailed and drug candidates are usually selected on their merits as being tight binders. However, the rationale that transient interactions are of importance for drug discovery is slowly gaining acceptance. These interactions may relate not only to the desired target interaction but also to unwanted interactions creating, for example, toxicity problems. Here we demonstrate, in a high-throughput screening format, affinity selection of weak binders to a model target of albumin by zonal retardation chromatography. It is perceived that this approach can define the "transient drug" as a complement to current drug discovery procedures.     

Not so different after all: a comparison of methods for detecting amino acid sites under selection

We consider three approaches for estimating the rates of nonsynonymous and synonymous changes at each site in a sequence alignment in order to identify sites under positive or negative selection: (1) a suite of fast likelihood-based "counting methods" that employ either a single most likely ancestral reconstruction, weighting across all possible ancestral reconstructions, or sampling from ancestral reconstructions; (2) a random effects likelihood (REL) approach, which models variation in nonsynonymous and synonymous rates across sites according to a predefined distribution, with the selection pressure at an individual site inferred using an empirical Bayes approach; and (3) a fixed effects likelihood (FEL) method that directly estimates nonsynonymous and synonymous substitution rates at each site. All three methods incorporate flexible models of nucleotide substitution bias and variation in both nonsynonymous and synonymous substitution rates across sites, facilitating the comparison between the methods. We demonstrate that the results obtained using these approaches show broad agreement in levels of Type I and Type II error and in estimates of substitution rates. Counting methods are well suited for large alignments, for which there is high power to detect positive and negative selection, but appear to underestimate the substitution rate. A REL approach, which is more computationally intensive than counting methods, has higher power than counting methods to detect selection in data sets of intermediate size but may suffer from higher rates of false positives for small data sets. A FEL approach appears to capture the pattern of rate variation better than counting methods or random effects models, does not suffer from as many false positives as random effects models for data sets comprising few sequences, and can be efficiently parallelized. Our results suggest that previously reported differences between results obtained by counting methods and random effects models arise due to a combination of the conservative nature of counting-based methods, the failure of current random effects models to allow for variation in synonymous substitution rates, and the naive application of random effects models to extremely sparse data sets. We demonstrate our methods on sequence data from the human immunodeficiency virus type 1 env and pol genes and simulated alignments.     

Block and inactivation of sodium channels in nerve by amino acid derivatives. II. Dependence on temperature and drug concentration.

The interaction of n-propylguanidinium (nPG) with sodium channels has been further characterized. From experiments at varying temperatures, the Q10 for the sodium current decay time constant in the two [Na+] gradients is 2.6-2.9 independent of drug. Testing several nPG concentrations we find that peak sodium current declines sharply with [nPG] at all levels, but the decay time constant approaches an asymptote above 4 mM. No "hooks" in sodium tail currents are seen. If the sodium current is allowed to decay completely before repolarization no tail current is observed. We have developed a kinetic model in which nPG acts at a single site within the sodium channel. Reaction of nPG with its receptor requires two steps. Fitting the temperature data shows that the first step involves diffusion of the drug to the site and close association with it. The second step may include molecular reorganization of the complex. The rate constants for the reaction are all simple exponential functions of voltage. Using them, the model successfully predicts decay time constants and peak currents, and their dependence on potential, [Na+] gradient, temperature, and nPG concentration. The results are consistent with the idea that an arginine residue may be closely associated with inactivation.     

Techniques for pharmacological and toxicological studies with isolated hepatocyte suspensions.

Since its introduction in 1969, the high-yield preparation of isolated hepatocytes has become a frequently used tool for the study of hepatic uptake, excretion, metabolism and toxicity of drugs and other xenobiotics. Basic preparative methods are now firmly established involving perfusion of the liver with a balanced-saline solution containing collagenase. Satisfactory procedures are available for determining cell yields, for expressing cellular activities and for establishing optimal incubation conditions. Gross cellular damage can be detected by means of trypan blue or by measuring enzyme leakage, and damaged cells can be removed from the preparation. Specialized techniques are available for preparing hepatocyte couplets and suspensions enriched with periportal or perivenous hepatocytes. The isolated hepatocyte preparation is particularly convenient for the study of the kinetics of hepatic drug uptake and excretion because the cells can be rapidly separated from the incubation medium. Isolated liver cells have also proved valuable for investigating drug metabolism since they show many of the features of the intact liver. However, they also show important differences such as losses of membrane specialization, some degree of cell polarity and the capacity to form bile. The many consequences of the hepatic toxicity of xenobiotics including lipid peroxidation, free radical formation, glutathione depletion, and covalent binding to macromolecules are also readily studied with the isolated liver cell preparation. A particular advantage is the ease with which morphological changes as a result of drug exposure can be observed in isolated hepatocytes. However, it must be remembered that the isolation procedure inevitably introduces changes that may make the cells more susceptible than the normal liver to damage by xenobiotic agents. Despite its limitations, the isolated hepatocyte preparation is now firmly established in the armamentarium of the investigator examining the interaction of the liver with xenobiotics.     

Interaction of a spin-labeled adenine-acridine conjugate with a DNA duplex containing an abasic site model

The abasic site is a common lesion in DNA that is also formed as an intermediate in the base excision repair of damaged bases. We have previously reported the adenine-acridine conjugate 1 that was designed to bind to the abasic site and interfere with the repair process. High-field NMR had shown that 1 forms specific complexes with a DNA duplex containing an apurinic abasic site model. We report here the dynamics of the interaction of the nitroxide-labeled analogue 3 of the conjugate 1 with the same apurinic oligonucleotide and with the parent unmodified duplex. Identical study of the labeled acridine subunit 5 used as a reference is also reported. In the presence of the apurinic duplex and depending on the concentrations and drug ratios, three species are observed: the radical "free in solution", the "intercalation" complex characterized by its similarity to that observed in the presence of the parent unmodified duplex, and the "abasic-site-specific" complex which is the sole species visible at low drug ratios. The experimental data reinforced by molecular modeling of the complex and theoretical calculation of correlation times suggest (i) the most immobilized form corresponds to that observed by NMR and (ii) complexation of the drug is little or not modified by the spin-label. We also show that the abasic site constitutes a binding site for the propylaminoacridine intercalator 5.     

Predicting Drug-Target Interactions for New Drug Compounds Using a Weighted Nearest Neighbor Profile

In silico discovery of interactions between drug compounds and target proteins is of core importance for improving the efficiency of the laborious and costly experimental determination of drug-target interaction. Drug-target interaction data are available for many classes of pharmaceutically useful target proteins including enzymes, ion channels, GPCRs and nuclear receptors. However, current drug-target interaction databases contain a small number of drug-target pairs which are experimentally validated interactions. In particular, for some drug compounds (or targets) there is no available interaction. This motivates the need for developing methods that predict interacting pairs with high accuracy also for these ''new'' drug compounds (or targets). We show that a simple weighted nearest neighbor procedure is highly effective for this task. We integrate this procedure into a recent machine learning method for drug-target interaction we developed in previous work. Results of experiments indicate that the resulting method predicts true interactions with high accuracy also for new drug compounds and achieves results comparable or better than those of recent state-of-the-art algorithms. Software is publicly available at http://cs.ru.nl/~tvanlaarhoven/drugtarget2013/.     

Chemical, biochemical and genetic endeavours characterizing the interaction of sparsomycin with the ribosome.

Sparsomycin interaction with the ribosome and characteristics of the drug binding site in the particle were studied using chemical modification of the drug, affinity labeling methods and isolation of drug resistant mutants. The structure-function relationship studies, performed with a large number of drug derivatives, indicate that the drug interacts with the ribosome by its western and eastern moieties. The uracil ring, in the western end of the drug molecule, probably forms hydrogen bonds with the rRNA, while the apolar CH3-S-CH3 group in the eastern end interacts with a hydrophobic ribosomal domain that affinity labeling results seem to indicate is formed by protein. An increase in lipophilicity in this part of the antibiotic results in a dramatic increase in the inhibitory activity of the drug. The sparsomycin binding site is not accessible in free ribosomes, but the presence of an N-blocked amino acyl-tRNA at the P-site turns the particles capable of reversible interaction with the drug. After failure using Escherichia coli, a sparsomycin-resistant mutant was obtained by direct mutagenesis on Halobacterium halobium, a species with a unique copy of rRNA genes, stressing the role of rRNA on the drug interaction site.     

Predicting new molecular targets for rhein using network pharmacology

Background

Drugs can influence the whole biological system by targeting interaction reactions. The existence of interactions between drugs and network reactions suggests a potential way to discover targets. The in silico prediction of potential interactions between drugs and target proteins is of core importance for the identification of new drugs or novel targets for existing drugs. However, only a tiny portion of drug-targets in current datasets are validated interactions. This motivates the need for developing computational methods that predict true interaction pairs with high accuracy. Currently, network pharmacology has used in identifying potential drug targets to predicting the spread of drug activity and greatly contributed toward the analysis of biological systems on a much larger scale than ever before.

Methods

In this article, we present a computational method to predict targets for rhein by exploring drug-reaction interactions. We have implemented a computational platform that integrates pathway, protein-protein interaction, differentially expressed genome and literature mining data to result in comprehensive networks for drug-target interaction. We used Cytoscape software for prediction rhein-target interactions, to facilitate the drug discovery pipeline.

Results

Results showed that 3 differentially expressed genes confirmed by Cytoscape as the central nodes of the complicated interaction network (99 nodes, 153 edges). Of note, we further observed that the identified targets were found to encompass a variety of biological processes related to immunity, cellular apoptosis, transport, signal transduction, cell growth and proliferation and metabolism.

Conclusions

Our findings demonstrate that network pharmacology can not only speed the wide identification of drug targets but also find new applications for the existing drugs. It also implies the significant contribution of network pharmacology to predict drug targets.     

Bayesian restoration of single-channel patch clamp recordings.

The technique of patch clamp recording makes possible the measurement of current flowing through a single ion channel in a cell membrane. Examination of such recordings suggests that the current is quantal in nature, alternating in a seemingly random manner between "on" and "off," but the recordings are corrupted by noise from a variety of sources. In this paper we propose and illustrate methods for restoring the underlying quantal signal from such noisy measurements. The methods use a Markov chain prior distribution for the underlying quantal process and base the restoration on the resulting posterior distribution.     

A null model for randomization tests of nestedness in species assemblages

Analysis of the degree of order in species assemblages in terms of nested subsets has received increased interest during the last decade. However, recently a series of papers have questioned the validity of methods employed for testing whether observed patterns deviate from random expectations. The current view seems to be that the randomization procedure should control for both number of species per site and species frequencies. The randomization procedures used also choose to keep the total number of observations constant in each resample. In this paper I question some of these assumptions when analyzing species-by-site matrices for detecting whether the biota is significantly nested or not. My basic assumption is that the observed species frequency is only an estimate of the probability of occurrence for the particular species. For a test of degree of nestedness all sites should be regarded as being equal. To what extent size, isolation or habitat quality may influence species distribution is a secondary question if nestedness can be statistically proven. This implies that generation of random matrices should only consider the frequency of the species (as an estimate of their probability of occurring in any patch). Such matrices are computationally simple and besides providing a test of nestedness also open the possibility of testing whether the range in species richness is smaller or larger than expected under random expectations. The choice of null model for the test should always be viewed in relation to the question asked. If nestedness is concerned the methods proposed here should be used. However, if other questions are at hand the restrictions of previous approaches may be valid. This is for instance the case if pairwise species co-occurrences are analyzed. In this case, the richness of each site should obviously be incorporated in the randomization to control for the higher probability of co-occurrence at species-rich sites.     

TrioMDR: Detecting SNP interactions in trio families with model-based multifactor dimensionality reduction

Single nucleotide polymorphism (SNP) interactions can explain the missing heritability of common complex diseases. Many interaction detection methods have been proposed in genome-wide association studies, and they can be divided into two types: population-based and family-based. Compared with population-based methods, family-based methods are robust vs. population stratification. Several family-based methods have been proposed, among which Multifactor Dimensionality Reduction (MDR)-based methods are popular and powerful. However, current MDR-based methods suffer from heavy computational burden. Furthermore, they do not allow for main effect adjustment. In this work we develop a two-stage model-based MDR approach (TrioMDR) to detect multi-locus interaction in trio families (i.e., two parents and one affected child). TrioMDR combines the MDR framework with logistic regression models to check interactions, so TrioMDR can adjust main effects. In addition, unlike consuming permutation procedures used in traditional MDR-based methods, TrioMDR utilizes a simple semi-parameter P-values correction procedure to control type I error rate, this procedure only uses a few permutations to achieve the significance of a multi-locus model and significantly speeds up TrioMDR. We performed extensive experiments on simulated data to compare the type I error and power of TrioMDR under different scenarios. The results demonstrate that TrioMDR is fast and more powerful in general than some recently proposed methods for interaction detection in trios. The R codes of TrioMDR are available at: https://github.com/TrioMDR/TrioMDR.     

Induced absorption band of holotransketolase and its interpretation

It has long been known that formation of a catalytically active holotransketolase from the apoenzyme and thiamine diphosphate (ThDP) is accompanied by appearance, in both the absorption and CD spectra, of a new band. Binding and subsequent conversion of transketolase substrates bring about changes in the intensity of this band. The observation of these changes allows the investigator to monitor the coenzyme-to-apoenzyme binding and the conversion of the substrates during the transketolase reaction and thus to kinetically characterize its individual steps. As regards the new absorption band induced by ThDP binding, its nature, until recently, remained unknown. The reason for its appearance was considered to be either the formation of a charge transfer complex between ThDP and tryptophan (phenylalanine) residue or stacking interaction between the residues of aromatic amino acids. They are thought to be brought together as a result of conformational changes of the apoenzyme during its interaction with the coenzyme. However none of these hypotheses had been substantiated experimentally. According to our hypothesis, the induced absorption band is that of the imino form of ThDP resulting from three contributing features of the ThDP binding site of transketolase: the relative hydrophobicity of this site, hydrogen bonding of the N1"-atom of the ThDP aminopyrimidine ring to Glu418, and base stacking interactions between the aminopyrimidine ring of ThDP and Phe445.     

Interactions of monovalent cations with sodium channels in squid axon. I. Modification of physiological inactivation gating

Inactivation of Na channels has been studied in voltage-clamped, internally perfused squid giant axons during changes in the ionic composition of the intracellular solution. Peak Na currents are reduced when tetramethylammonium ions (TMA+) are substituted for Cs ions internally. The reduction reflects a rapid, voltage-dependent block of a site in the channel by TMA+. The estimated fractional electrical distance for the site is 10% of the channel length from the internal surface. Na tail currents are slowed by TMA+ and exhibit kinetics similar to those seen during certain drug treatments. Steady state INa is simultaneously increased by TMA+, resulting in a "cross-over" of current traces with those in Cs+ and in greatly diminished inactivation at positive membrane potentials. Despite the effect on steady state inactivation, the time constants for entry into and exit from the inactivated state are not significantly different in TMA+ and Cs+. Increasing intracellular Na also reduces steady state inactivation in a dose-dependent manner. Ratios of steady state INa to peak INa vary from approximately 0.14 in Cs+- or K+-perfused axons to approximately 0.4 in TMA+- or Na+-perfused axons. These results are consistent with a scheme in which TMA+ or Na+ can interact with a binding site near the inner channel surface that may also be a binding or coordinating site for a natural inactivation particle. A simple competition between the ions and an inactivation particle is, however, not sufficient to account for the increase in steady state INa, and changes in the inactivation process itself must accompany the interaction of TMA+ and Na+ with the channel.     

PCR fragmentation of DNA

A method has been developed to prepare random DNA fragments using PCR. First, two cycles are carried out at 16 degrees C with the Klenow's fragment and oligonucleotides (random primers) with random 3'-sequences and the 5'-constant part containing the site for cloning with the site-specific endonuclease. The random primers can link to any DNA site, and random DNA fragments are formed during DNA synthesis. During the second cycle, after denaturation of the DNA and addition of the Klenow's fragment, the random primers can link to newly synthesized DNA strands, and after DNA synthesis single-stranded DNA fragments are produced which have a constant primer sequence at the 5'-end and a complementary to it sequence at the 3'-end. During the third cycle, the constant primer is added and double-stranded fragments with the constant primer sequences at both ends are formed during DNA synthesis. Incubation for 1 h at 37 degrees C degrades the oligonucleotides used at the first stage due to endonuclease activity of the Klenow's fragment. Then routine PCR amplification is carried out using the constant primer. This method is more advantageous than hydrodynamic methods of DNA fragmentation widely used for "shotgun" cloning.     

Methods for site-specific drug conjugation to antibodies

Antibody drug conjugates (ADCs) are an emerging class of targeted therapeutics with the potential to improve therapeutic index over traditional chemotherapy. Drugs and linkers have been the current focus of ADC development, in addition to antibody and target selection. Recently, however, the importance of conjugate homogeneity has been realized. The current methods for drug attachment lead to a heterogeneous mixture, and some populations of that mixture have poor in vivo performance. New methods for site-specific drug attachment lead to more homogeneous conjugates and allow control of the site of drug attachment. These subtle improvements can have profound effects on in vivo efficacy and therapeutic index. This review examines current methods for site-specific drug conjugation to antibodies, and compares in vivo results with their non-specifically conjugated counterparts. The apparent improvement in pharmacokinetics and the reduced off target toxicity warrant further development of this site-specific modification approach for future ADC development.     

Kinetics of the daunomycin--DNA interaction

The kinetics of the interaction of daunomycin with calf thymus DNA are described. Stopped-flow and temperature-jump relaxation methods, using absorption detection, were used to study the binding reaction. Three relaxation times were observed, all of which are concentration dependent, although the two slower relaxations approach constant values at high reactant concentrations. Relaxation times over a wide range of concentrations were gathered, and the data were fit by a minimal mechanism in which a rapid bimolecular association step is followed by two sequential isomerization steps. The six rate constants for this mechanism were extracted from our data by relaxation analysis. The values determined for the six rate constants may be combined to calculate an overall equilibrium constant that is in excellent agreement with that obtained by independent equilibrium measurements. Additional stopped-flow experiments, using first sodium dodecyl sulfate to dissociate bound drug and second pseudo-first-order conditions to study the fast bimolecular step, provide independent verification of three of the six rate constants. The temperature dependence of four of the six rate constants was measured, allowing estimates of the activation energy of some of the steps to be made. We speculate that the three steps in the proposed mechanism may correspond to a rapid "outside" binding of daunomycin to DNA, followed by intercalation of the drug, followed by either conformational adjustment of the drug or DNA binding site or redistribution of bound drug to preferred sites.     

Stimulation of P-glycoprotein ATPase by analogues of tetramethylrosamine: coupling of drug binding at the "R" site to the ATP hydrolysis transition state

The multidrug resistance efflux pump P-glycoprotein (Pgp) couples drug export to ATP binding and hydrolysis. Details regarding drug trajectory, as well as the molecular basis for coupling, remain unknown. Nearly all drugs exported by Pgp have been assayed for competitive behavior with rhodamine123 transport at a canonical "R" drug binding site. Tetramethylrosamine (TMR) displays a relatively high affinity for Pgp when compared to other rhodamines. Here, we present the construction and characterization of a library of compounds based upon the TMR scaffold and use this set to assess the determinants of drug binding to the "R" site of Pgp. This set contained modifications in (1) the number, location, and conformational mobility of hydrogen-bond acceptors; (2) the heteroatom in the xanthylium core; and (3) the size of the substituent in the 9-position of the xanthylium core. Relative specificity for coupling to the distal ATP catalytic site was assessed by ATPase stimulation. We found marked ( approximately 1000-fold) variation in the ATPase specificity constant within the library of TMR analogues. Using established methods involving ADP-Vi trapping by wild-type Pgp and ATP binding by catalytic carboxylate mutant Pgp, these effects can be extended to ATP hydrolysis transition-state stabilization and ATP occlusion at a single site. These data support the idea that drugs trigger the engagement of ATP catalytic site residues necessary for hydrolysis. Further, the nature of the drug binding site and coupling mechanism may be dissected by variation of a drug-like scaffold. These studies may facilitate development of novel competitive inhibitors at the "R" drug site.     

Methods for site-specific drug conjugation to antibodies

Antibody drug conjugates (ADCs) are an emerging class of targeted therapeutics with the potential to improve therapeutic index over traditional chemotherapy. Drugs and linkers have been the current focus of ADC development, in addition to antibody and target selection. Recently, however, the importance of conjugate homogeneity has been realized. The current methods for drug attachment lead to a heterogeneous mixture, and some populations of that mixture have poor in vivo performance. New methods for site-specific drug attachment lead to more homogeneous conjugates and allow control of the site of drug attachment. These subtle improvements can have profound effects on in vivo efficacy and therapeutic index. This review examines current methods for site-specific drug conjugation to antibodies, and compares in vivo results with their non-specifically conjugated counterparts. The apparent improvement in pharmacokinetics and the reduced off target toxicity warrant further development of this site-specific modification approach for future ADC development.     

Probabilistic secretion of quanta from nerve terminals at synaptic sites on muscle cells: non-uniformity, autoinhibition and the binomial hypothesis

A model of the secretion of a quantum at a release site is proposed in which, following the influx of calcium ions, synaptic vesicles are made available for release by the activation of kappa phosphorylation steps with rate alpha. At any time during this process the vesicles may become unavailable for secretion at rate gamma. On completion of the kappa phosphorylation steps the vesicles participate in the formation of a fusion pore with the terminal membrane to give exocytosis at rate delta. Changes in alpha, delta and kappa are shown to produce characteristic changes in the number and timecourse of quantal secretions following a nerve impulse, which are similar to those observed following drug treatments that are thought to act selectively on each of these processes. The number of quanta secreted from nerve terminals that consist of many release sites does not fluctuate much during a low frequency train of impulses: the variance is small compared with the mean level, so secretion follows binomial rather than Poisson statistics. A theory is derived that shows that variations in the probability of secretion amongst these release sites of any particular kind fails to reduce the variance of the total secretion from the terminal; Poisson rather than binomial statistics then still apply. The theory shows that an interaction between release sites is required to reduce this variance and such an effect is provided if secretion at a site inhibits secretion at nearby sites. Simulations show that incorporating this process of autoinhibition into the model reproduces the experimental observations on the effects of calcium ions on the binomial parameters p and n as well as on the relative constancy of p during facilitation and depression of quantal secretion. Methods for estimating the timecourse of changes in the probability of secretion at release sites following an impulse, by using either the time of occurrence of first, second, third or later quantal latencies, are given. These procedures show that current methods for estimating the time-dependent probability changes are inadequate for detecting interaction between release sites, such as autoinhibition, unless this is relatively large. Therefore, estimates from third quantal latencies are used.     

Run-Reversal Equilibrium for Clinical Trial Randomization

In this paper, we describe a new restricted randomization method called run-reversal equilibrium (RRE), which is a Nash equilibrium of a game where (1) the clinical trial statistician chooses a sequence of medical treatments, and (2) clinical investigators make treatment predictions. RRE randomization counteracts how each investigator could observe treatment histories in order to forecast upcoming treatments. Computation of a run-reversal equilibrium reflects how the treatment history at a particular site is imperfectly correlated with the treatment imbalance for the overall trial. An attractive feature of RRE randomization is that treatment imbalance follows a random walk at each site, while treatment balance is tightly constrained and regularly restored for the overall trial. Less predictable and therefore more scientifically valid experiments can be facilitated by run-reversal equilibrium for multi-site clinical trials.