Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

Design, synthesis and cellular and molecular pharmacology of 2$\sp\prime$-bromo-4$\sp\prime$-epidaunorubicin (WP401): A novel non-cross-resistant anthracycline antibiotic with the intrinsic ability to circumvent P -glycoprotein-mediated drug resistance

We designed and synthesized a novel daunorubicin (DNR) analogue that effectively circumvents P-glycoprotein (P-gp)-mediated drug resistance. The fully protected carbohydrate intermediate 1,2-dibromoacosamine was prepared from acosamine and effectively coupled to daunomycinone in high yield. Deprotection under alkaline conditions yielded 2$\sp\prime$-bromo-4$\sp\prime$-epidaunorubicin (WP401). The in vitro cytotoxicity and cellular and molecular pharmacology of WP401 were compared with those of DNR in a panel of wild-type cell lines (KB-3-1, P388S, and HL60S) and their multidrug-resistant (MDR) counterparts (KB-V1, P388/DOX, and HL60/DOX). Fluorescent spectrophotometry, flow cytometry, and confocal laser scanning microscopy were used to measure intracellular accumulation, retention, and subcellular distribution of these agents. All MDR cell lines exhibited reduced DNR uptake that was restored, upon incubation with either verapamil (VER) or cyclosporin A (CSA), to the level found in sensitive cell lines. In contrast, the uptake of WP401 was essentially the same in the absence or presence of VER or CSA in all tested cell lines. The in vitro cytotoxicity of WP401 was similar to that of DNR in the sensitive cell lines but significantly higher in resistant cell lines (resistance index (RI) of 2-6 for WP401 vs 75-85 for DNR). To ascertain whether drug-mediated cytotoxicity and retention were accompanied by DNA strand breaks, DNA single- and double-strand breaks were assessed by alkaline elution. High levels of such breaks were obtained using 0.1-2 $\mu$g/mL of WP401 in both sensitive and resistant cells. In contrast, DNR caused strand breaks only in sensitive cells and not much in resistant cells. We also compared drug-induced DNA fragmentation similar to that induced by DNR. However, in P-gp-positive cells, WP401 induced 2- to 5-fold more DNA fragmentation than DNR. This increased DNA strand breakage by WP401 was correlated with its increased uptake and cytotoxicity in these cell lines. Overall these results indicate that WP401 is more cytotoxic than DNR in MDR cells and that this phenomenon might be related to the reduced basicity of the amino group and increased lipophilicity of WP401. ^

Immobilisation dose rates and photogrammetry based mass calculations for southern elephant seals during expedition JUB2010

Immobilization and anaesthesia of adult male southern elephant seals (Mirounga leonina) is potentially risky for animals and scientists. A tiletamine/zolazepam injection is considered the most appropriate drug combination for field application in this species. Since appropriate dosages are difficult to assess due to uncertainties in weight estimation, we used photogrammetry-derived weight estimates to ensure precise post hoc calculations of dosages. We report on 15 intramuscular tiletamine/zolazepam immobilizations of post-moult males of the upper weight class at King George Island/Isla 25 de Mayo, in April 2010. Initial injections were made using blowpipe syringes. Mean tiletamine/zolazepam combined dosages of 0.71 mg/kg (SD ± 0.16) ranged between 0.46 and 1.01 mg/kg. In four cases, ketamine was added in dosages between 0.96 and 2.61 mg/kg. Mean induction period was 23 min (± 15), and the mean duration of the procedures from first injection to release of the animals required 96 min (± 51). Four seals exhibited periods of apnoea, and one case of an extended, repetitive, and potentially critical apnoea (> 25 and 8 min) required intervention in order to successfully re-initiate spontaneous respiration. All procedures resulted in proper immobilizations allowing for the deployment of the satellite tags on the seals' heads. The fact that even substantial deviations between the initial weight estimates and the photogrammetry-derived weight estimates had no apparent effect on the course of the immobilization underlines the drugs' wide safety margin in this species.

Equivalent dose, dose rate and age determination of cobble surfaces in the Antarctic

Most current methods of reconstructing past sea levels within Antarctica rely on radiocarbon dating. However, radiocarbon dating is limited by the availability of material for dating and problems inherent with radiocarbon reservoirs in Antarctic marine systems. Here we report on the success of a new approach to dating raised beach deposits in Antarctica for the purpose of reconstructing past sea levels. This new approach is the use of optically stimulated luminescence (OSL) on quartz-grains obtained from the underside of cobbles within raised beaches and boulder pavements. We obtained eight OSL dates from three sites along the shores of Maxwell Bay in the South Shetland Islands of the Antarctic Peninsula. These dates are internally consistent and fit well with previously published radiocarbon ages obtained from the same deposits. In addition, when the technique was applied to a modern beach, it resulted in an age of zero. Our results suggest that this method will provide a valuable tool in the reconstruction of past sea levels in Antarctica and other coarse-grained beach deposits across the globe.

Analysis of the operator action and the single failure criteria in a SGTR sequence using best estimate assumptions with TRACE 5.0

Steam Generator Tube Rupture (SGTR) sequences in Pressurized Water Reactors are known to be one of the most demanding transients for the operating crew. SGTR are a special kind of transient as they could lead to radiological releases without core damage or containment failure, as they can constitute a direct path from the reactor coolant system to the environment. The first methodology used to perform the Deterministic Safety Analysis (DSA) of a SGTR did not credit the operator action for the first 30 min of the transient, assuming that the operating crew was able to stop the primary to secondary leakage within that period of time. However, the different real SGTR accident cases happened in the USA and over the world demonstrated that the operators usually take more than 30 min to stop the leakage in actual sequences. Some methodologies were raised to overcome that fact, considering operator actions from the beginning of the transient, as it is done in Probabilistic Safety Analysis. This paper presents the results of comparing different assumptions regarding the single failure criteria and the operator action taken from the most common methodologies included in the different Deterministic Safety Analysis. One single failure criteria that has not been analysed previously in the literature is proposed and analysed in this paper too. The comparison is done with a PWR Westinghouse three loop model in TRACE code (Almaraz NPP) with best estimate assumptions but including deterministic hypothesis such as single failure criteria or loss of offsite power. The behaviour of the reactor is quite diverse depending on the different assumptions made regarding the operator actions. On the other hand, although there are high conservatisms included in the hypothesis, as the single failure criteria, all the results are quite far from the regulatory limits. In addition, some improvements to the Emergency Operating Procedures to minimize the offsite release from the damaged SG in case of a SGTR are outlined taking into account the offsite dose sensitivity results.

Choroid plexus epithelial expression of MDR1 P glycoprotein and multidrug resistance-associated protein contribute to the blood–cerebrospinal-fluid drug-permeability barrier

The blood–brain barrier and a blood–cerebrospinal-fluid (CSF) barrier function together to isolate the brain from circulating drugs, toxins, and xenobiotics. The blood–CSF drug-permeability barrier is localized to the epithelium of the choroid plexus (CP). However, the molecular mechanisms regulating drug permeability across the CP epithelium are defined poorly. Herein, we describe a drug-permeability barrier in human and rodent CP mediated by epithelial-specific expression of the MDR1 (multidrug resistance) P glycoprotein (Pgp) and the multidrug resistance-associated protein (MRP). Noninvasive single-photon-emission computed tomography with 99mTc-sestamibi, a membrane-permeant radiopharmaceutical whose transport is mediated by both Pgp and MRP, shows a large blood-to-CSF concentration gradient across intact CP epithelium in humans in vivo. In rats, pharmacokinetic analysis with 99mTc-sestamibi determined the concentration gradient to be greater than 100-fold. In membrane fractions of isolated native CP from rat, mouse, and human, the 170-kDa Pgp and 190-kDa MRP are identified readily. Furthermore, the murine proteins are absent in CP isolated from their respective mdr1a/1b(−/−) and mrp(−/−) gene knockout littermates. As determined by immunohistochemical and drug-transport analysis of native CP and polarized epithelial cell cultures derived from neonatal rat CP, Pgp localizes subapically, conferring an apical-to-basal transepithelial permeation barrier to radiolabeled drugs. Conversely, MRP localizes basolaterally, conferring an opposing basal-to-apical drug-permeation barrier. Together, these transporters may coordinate secretion and reabsorption of natural product substrates and therapeutic drugs, including chemotherapeutic agents, antipsychotics, and HIV protease inhibitors, into and out of the central nervous system.

Inter- and intraspecies polymorphisms in the cholecystokinin-B/gastrin receptor alter drug efficacy

The brain cholecystokinin-B/gastrin receptor (CCK-BR) is a major target for drug development because of its postulated role in modulating anxiety, memory, and the perception of pain. Drug discovery efforts have resulted in the identification of small synthetic molecules that can selectively activate this receptor subtype. These drugs include the peptide-derived compound PD135,158 as well as the nonpeptide benzodiazepine-based ligand, L-740,093 (S enantiomer). We now report that the maximal level of receptor-mediated second messenger signaling that can be achieved by these compounds (drug efficacy) markedly differs among species homologs of the CCK-BR. Further analysis reveals that the observed differences in drug efficacy are in large part explained by single or double aliphatic amino acid substitutions between respective species homologs. This interspecies variability in ligand efficacy introduces the possibility of species differences in receptor-mediated function, an important consideration when selecting animal models for preclinical drug testing. The finding that even single amino acid substitutions can significantly affect drug efficacy prompted us to examine ligand-induced signaling by a known naturally occurring human CCK-BR variant (glutamic acid replaced by lysine in position 288; 288E → K). When examined using the 288E → K receptor, the efficacies of both PD135,158 and L-740,093 (S) were markedly increased compared with values obtained with the wild-type human protein. These observations suggest that functional variability resulting from human receptor polymorphisms may contribute to interindividual differences in drug effects.

Graded transcriptional response to different concentrations of a single transactivator

Threshold mechanisms of transcriptional activation are thought to be critical for translating continuous gradients of extracellular signals into discrete all-or-none cellular responses, such as mitogenesis and differentiation. Indeed, unequivocal evidence for a graded transcriptional response in which the concentration of inducer directly correlates with the level of gene expression in individual eukaryotic cells is lacking. By using a novel binary tetracycline regulatable retroviral vector system, we observed a graded rather than a threshold mechanism of transcriptional activation in two different model systems. When polyclonal populations of cells were analyzed at the single cell level, a dose-dependent, stepwise increase in expression of the reporter gene, green fluorescent protein (GFP), was observed by fluorescence-activated cell sorting. These data provide evidence that, in addition to the generally observed all-or-none switch, the basal transcription machinery also can respond proportionally to changes in concentration of extracellular inducers and trancriptional activators.

Inverse radiation dose-rate effects on somatic and germ-line mutations and DNA damage rates

The mutagenic effect of low linear energy transfer ionizing radiation is reduced for a given dose as the dose rate (DR) is reduced to a low level, a phenomenon known as the direct DR effect. Our reanalysis of published data shows that for both somatic and germ-line mutations there is an opposite, inverse DR effect, with reduction from low to very low DR, the overall dependence of induced mutations being parabolically related to DR, with a minimum in the range of 0.1 to 1.0 cGy/min (rule 1). This general pattern can be attributed to an optimal induction of error-free DNA repair in a DR region of minimal mutability (MMDR region). The diminished activation of repair at very low DRs may reflect a low ratio of induced (“signal”) to spontaneous background DNA damage (“noise”). Because two common DNA lesions, 8-oxoguanine and thymine glycol, were already known to activate repair in irradiated mammalian cells, we estimated how their rates of production are altered upon radiation exposure in the MMDR region. For these and other abundant lesions (abasic sites and single-strand breaks), the DNA damage rate increment in the MMDR region is in the range of 10% to 100% (rule 2). These estimates suggest a genetically programmed optimatization of response to radiation in the MMDR region.

Nuclear foci of mammalian recombination proteins are located at single-stranded DNA regions formed after DNA damage

A sensitive and rapid in situ method was developed to visualize sites of single-stranded (ss) DNA in cultured cells and in experimental test animals. Anti-bromodeoxyuridine antibody recognizes the halogenated base analog incorporated into chromosomal DNA only when substituted DNA is in the single strand form. After treatment of cells with DNA-damaging agents or γ irradiation, ssDNA molecules form nuclear foci in a dose-dependent manner within 60 min. The mammalian recombination protein Rad51 and the replication protein A then accumulate at sites of ssDNA and form foci, suggesting that these are sites of recombinational DNA repair.

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

Relationship between psychostimulant-induced "high" and dopamine transporter occupancy.

The ability of cocaine to inhibit the dopamine transporter (DAT) appears to be crucial for its reinforcing properties. The potential use of drugs that produce long-lasting inhibition of the DAT as a mean of preventing the "high" and reducing drug-seeking behavior has become a major strategy in medication development. However, neither the relation between the high and DAT inhibition nor the ability to block the high by prior DAT blockade have ever been demonstrated. To evaluate if DAT could prevent the high induced by methylphenidate (MP), a drug which like cocaine inhibits the DAT, we compared the responses in eight non-drug-abusing subjects between the first and the second of two MP doses (0.375 mg/kg, i.v.) given 60 min apart. At 60 min the high from MP has returned to baseline, but 75-80% of the drug remains in brain. Positron-emission tomography and [11C]d-threo-MP were used to estimate DAT occupancies at different times after MP. DAT inhibition by MP did not block or attenuate the high from a second dose of MP given 60 min later, despite a 80% residual transporter occupancy from the first dose. Furthermore some subjects did not perceive a high after single or repeated administration despite significant DAT blockade. These results indicate that DAT occupancy is not sufficient to account for the high, and that for DAT inhibitors to be therapeutically effective, occupancies > 80% may be required.

Host-parasite dynamics and outgrowth of virus containing a single K70R amino acid change in reverse transcriptase are responsible for the loss of human immunodeficiency virus type 1 RNA load suppression by zidovudine.

The association between human immunodeficiency virus type I (HIV-1) RNA load changes and the emergence of resistant virus variants was investigated in 24 HIV-1-infected asymptomatic persons during 2 years of treatment with zidovudine by sequentially measuring serum HIV-1 RNA load and the relative amounts of HIV-1 RNA containing mutations at reverse transcriptase (RT) codons 70 (K-->R), 41 (M-->L), and 215 (T-->Y/F). A mean maximum decline in RNA load occurred during the first month, followed by a resurgence between 1 and 3 months, which appeared independent of drug-resistance. Mathematical modeling suggests that this resurgence is caused by host-parasite dynamics, and thus reflects infection of the transiently increased numbers of CD4+ lymphocytes. Between 3 and 6 months of treatment, the RNA load returned to baseline values, which was associated with the emergence of virus containing a single lysine to arginine amino acid change at RT codon 70, only conferring an 8-fold reduction in susceptibility. Despite the relative loss of RNA load suppression, selection toward mutations at RT codons 215 and 41 continued. Identical patterns were observed in the mathematical model. While host-parasite dynamics and outgrowth of low-level resistant virus thus appear responsible for the loss of HIV-1 RNA load suppression, zidovudine continues to select for alternative mutations, conferring increasing levels of resistance.

Molecular determinants of drug access to the receptor site for antiarrhythmic drugs in the cardiac Na+ channel.

The clinical efficacy of local anesthetic and antiarrhythmic drugs is due to their voltage- and frequency-dependent block of Na+ channels. Quaternary local anesthetic analogs such as QX-314, which are permanently charged and membrane-impermeant, effectively block cardiac Na+ channels when applied from either side of the membrane but block neuronal Na+ channels only from the intracellular side. This difference in extracellular access to QX-314 is retained when rat brain rIIA Na+ channel alpha subunits and rat heart rH1 Na+ channel alpha subunits are expressed transiently in tsA-201 cells. Amino acid residues in transmembrane segment S6 of homologous domain IV (IVS6) of Na+ channel alpha subunits have important effects on block by local anesthetic drugs. Although five amino acid residues in IVS6 differ between brain rIIA and cardiac rH1, exchange of these amino acid residues by site-directed mutagenesis showed that only conversion of Thr-1755 in rH1 to Val as in rIIA was sufficient to reduce the rate and extent of block by extracellular QX-314 and slow the escape of drug from closed channels after use-dependent block. Tetrodotoxin also reduced the rate of block by extracellular QX-314 and slowed escape of bound QX-314 via the extracellular pathway in rH1, indicating that QX-314 must move through the pore to escape. QX-314 binding was inhibited by mutation of Phe-1762 in the local anesthetic receptor site of rH1 to Ala whether the drug was applied extracellularly or intracellularly. Thus, QX-314 binds to a single site in the rH1 Na+ channel alpha subunit that contains Phe-1762, whether it is applied from the extracellular or intracellular side of the membrane. Access to that site from the extracellular side of the pore is determined by the amino acid at position 1755 in the rH1 cardiac Na+ channel.