Inhibitor selectivity of CNTs and ENTs.

The concentrative nucleoside transporters (CNT; solute carrier family 28 (SLC28)) and the equilibrative nucleoside transporters (ENT; solute carrier family 29 (SLC29)) are important therapeutic targets but may also mediate toxicity or adverse events. To explore the relative role of the base and the monosaccharide moiety in inhibitor selectivity we selected compounds that either harbor an arabinose moiety or a cytosine moiety, as these groups had several commercially available drug members. The screening data showed that more compounds harboring a cytosine moiety displayed potent interactions with the CNTs than compounds harboring the arabinose moiety. In contrast, ENTs showed a preference for compounds with an arabinose moiety. The correlation between CNT1 and CNT3 was good as five of six compounds displayed IC50 values within the threefold threshold and one displayed a borderline 4-fold difference. For CNT1 and CNT2 as well as for CNT2 and CNT3 only two of six IC50 values correlated and one displayed a borderline 4-fold difference. Interestingly, of the six compounds that potently interacted with both ENT1 and ENT2 only nelarabine displayed selectivity. Our data show differences between inhibitor selectivities of CNTs and ENTs as well as differences within the CNT family members.

Three-dimensional structure of the human breast cancer resistance protein (BCRP/ABCG2) in an inward-facing conformation.

ABCG2 is an efflux drug transporter that plays an important role in drug resistance and drug disposition. In this study, the first three-dimensional structure of human full-length ABCG2 analysed by electron crystallography from two-dimensional crystals in the absence of nucleotides and transported substrates is reported at 2 nm resolution. In this state, ABCG2 forms a symmetric homodimer with a noncrystallographic twofold axis perpendicular to the two-dimensional crystal plane, as confirmed by subtomogram averaging. This configuration suggests an inward-facing configuration similar to murine ABCB1, with the nucleotide-binding domains (NBDs) widely separated from each other. In the three-dimensional map, densities representing the long cytoplasmic extensions from the transmembrane domains that connect the NBDs are clearly visible. The structural data have allowed the atomic model of ABCG2 to be refined, in which the two arms of the V-shaped ABCG2 homodimeric complex are in a more closed and narrower conformation. The structural data and the refined model of ABCG2 are compatible with the biochemical analysis of the previously published mutagenesis studies, providing novel insight into the structure and function of the transporter.

Complex activities of the human Bloom's syndrome helicase are encoded in a core region comprising the RecA and Zn-binding domains.

Bloom's syndrome DNA helicase (BLM), a member of the RecQ family, is a key player in homologous recombination (HR)-based error-free DNA repair processes. During HR, BLM exerts various biochemical activities including single-stranded (ss) DNA translocation, separation and annealing of complementary DNA strands, disruption of complex DNA structures (e.g. displacement loops) and contributes to quality control of HR via clearance of Rad51 nucleoprotein filaments. We performed a quantitative mechanistic analysis of truncated BLM constructs that are shorter than the previously identified minimal functional module. Surprisingly, we found that a BLM construct comprising only the two conserved RecA domains and the Zn(2+)-binding domain (residues 642-1077) can efficiently perform all mentioned HR-related activities. The results demonstrate that the Zn(2+)-binding domain is necessary for functional interaction with DNA. We show that the extensions of this core, including the winged-helix domain and the strand separation hairpin identified therein in other RecQ-family helicases, are not required for mechanochemical activity per se and may instead play modulatory roles and mediate protein-protein interactions.

Predicting substrates of the human breast cancer resistance protein using a support vector machine method.

BACKGROUND: Human breast cancer resistance protein (BCRP) is an ATP-binding cassette (ABC) efflux transporter that confers multidrug resistance in cancers and also plays an important role in the absorption, distribution and elimination of drugs. Prediction as to if drugs or new molecular entities are BCRP substrates should afford a cost-effective means that can help evaluate the pharmacokinetic properties, efficacy, and safety of these drugs or drug candidates. At present, limited studies have been done to develop in silico prediction models for BCRP substrates. In this study, we developed support vector machine (SVM) models to predict wild-type BCRP substrates based on a total of 263 known BCRP substrates and non-substrates collected from literature. The final SVM model was integrated to a free web server. RESULTS: We showed that the final SVM model had an overall prediction accuracy of ~73% for an independent external validation data set of 40 compounds. The prediction accuracy for wild-type BCRP substrates was ~76%, which is higher than that for non-substrates. The free web server (http://bcrp.althotas.com) allows the users to predict whether a query compound is a wild-type BCRP substrate and calculate its physicochemical properties such as molecular weight, logP value, and polarizability. CONCLUSIONS: We have developed an SVM prediction model for wild-type BCRP substrates based on a relatively large number of known wild-type BCRP substrates and non-substrates. This model may prove valuable for screening substrates and non-substrates of BCRP, a clinically important ABC efflux drug transporter.

Evaluation of drug-human serum albumin binding interactions with support vector machine aided online automated docking.

MOTIVATION: Human serum albumin (HSA), the most abundant plasma protein is well known for its extraordinary binding capacity for both endogenous and exogenous substances, including a wide range of drugs. Interaction with the two principal binding sites of HSA in subdomain IIA (site 1) and in subdomain IIIA (site 2) controls the free, active concentration of a drug, provides a reservoir for a long duration of action and ultimately affects the ADME (absorption, distribution, metabolism, and excretion) profile. Due to the continuous demand to investigate HSA binding properties of novel drugs, drug candidates and drug-like compounds, a support vector machine (SVM) model was developed that efficiently predicts albumin binding. Our SVM model was integrated to a free, web-based prediction platform (http://albumin.althotas.com). Automated molecular docking calculations for prediction of complex geometry are also integrated into the web service. The platform enables the users (i) to predict if albumin binds the query ligand, (ii) to determine the probable ligand binding site (site 1 or site 2), (iii) to select the albumin X-ray structure which is complexed with the most similar ligand and (iv) to calculate complex geometry using molecular docking calculations. Our SVM model and the potential offered by the combined use of in silico calculation methods and experimental binding data is illustrated.

Functional adaptation of the switch-2 nucleotide sensor enables rapid processive translocation by myosin-5.

Active site loops that are conserved across superfamilies of myosins, kinesins, and G proteins play key roles in allosteric coupling of NTP hydrolysis to interaction with track filaments or effector proteins. In this study, we investigated how the class-specific natural variation in the switch-2 active site loop contributes to the motor function of the intracellular transporter myosin-5. We used single-molecule, rapid kinetic and spectroscopic experiments and semiempirical quantum chemical simulations to show that the class-specific switch-2 structure including a tyrosine (Y439) in myosin-5 enables rapid processive translocation along actin filaments by facilitating Mg(2+)-dependent ADP release. Using wild-type control and Y439 point mutant myosin-5 proteins, we demonstrate that the translocation speed precisely correlates with the kinetics of nucleotide exchange. Switch-2 variants can thus be used to fine-tune translocation speed while maintaining high processivity. The class-specific variation of switch-2 in various NTPase superfamilies indicates its general role in the kinetic tuning of Mg(2+)-dependent nucleotide exchange.

Sulfobutylether-beta-cyclodextrin-enabled antiviral remdesivir: Characterization of electrospun- and lyophilized formulations.

Veklury™ by Gilead Sciences, Inc., containing antiviral drug, remdesivir (REM) has received emergency authorization in the USA and in Europe for COVID-19 therapy. Here, for the first time, we describe details of the non-covalent, host-guest type interaction between REM and the solubilizing excipient, sulfobutylether-beta-cyclodextrin (SBECD) that results in significant solubility enhancement. Complete amorphousness of the cyclodextrin-enabled REM formulation was demonstrated by X-ray diffraction, thermal analysis, Raman chemical mapping and electron microscopy/energy dispersive spectroscopy. The use of solubilizing carbohydrate resulted in a 300-fold improvement of the aqueous solubility of REM, and enhanced dissolution rate of the drug enabling the preparation of stable infusion solutions for therapy. 2D ROESY NMR spectroscopy provided information on the nature of REM-excipient interaction and indicated the presence of inclusion phenomenon and the electrostatic attraction between anionic SBECD and nitrogen-containing REM in aqueous solution.

Cyclodextrin KnowledgeBase a web-based service managing CD-ligand complexation data.

Cyclodextrins are cyclic oligosaccharides that are able to form water-soluble inclusion complexes with small molecules. Because of their complexing ability, they are widely applied in food, pharmaceutical and chemical industries. In this paper we describe the development of a free web-service, Cyclodextrin KnowledgeBase: ( http://www.cyclodextrin.net ). The database contains four modules: the Publication, Interaction, Chirality and Analysis Modules. In the Publication Module, almost 50,000 publication details are collected that can be retrieved by text search. In the Interaction and Chirality Modules relevant literature data on cyclodextrin complexation and chiral recognition are collected that can be retrieved by both text and structural searches. Moreover, in the Analysis Module, the geometries of small molecule-cyclodextrin complexes can be predicted using molecular docking tools in order to explore the structures and interaction energies of the inclusion complexes. Complex geometry prediction is made possible by the built-in database of 95 cyclodextrin derivatives, where the 3D structures as well as the partial charges are calculated and stored for further utilization. The use of the database is demonstrated by several examples.

In silico description of differential enantioselectivity in methoxychlor O-demethylation by CYP2C enzymes.

Methoxychlor undergoes metabolism by cytochrome P450 (CYP) enzymes forming a chiral mono-phenolic derivative (Mono-OH-M) as main metabolite. In the current study, members of the CYP2C family were examined for their chiral preference in Mono-OH-M formation. CYP2C9 and CYP2C19 possessed high enantioselectivity favoring the formation of S-Mono-OH-M; CYP2C3 showed no enantioselectivity, whereas CYP2C5 slightly favored the formation of R-Mono-OH-M. Molecular modeling calculations were utilized in order to explain the observed differences in chiral preference of CYP2C enzymes. Molecular docking calculations could describe neither the existence of chiral preference in metabolism, nor the enantiomer which is preferentially formed. Molecular dynamic calculations were also carried out and were found to be useful for accurate description of chiral preference in biotransformation of methoxychlor by CYP2C enzymes. An in silico model capable of predicting chiral preference in cytochrome P450 enzymes in general can be developed based on the analysis of the stability and rigidity parameters of interacting partners during molecular dynamic simulation.

Organogold complexes probe a large beta-barrel cavity for human serum alpha1-acid glycoprotein.

Human alpha(1)-acid glycoprotein (AAG) is an acute phase component of the plasma, binding numerous drugs and natural compounds with high-affinity. Using circular dichroism (CD) spectroscopy, strong AAG binding of organogold complexes was found, the molecular size and chemical structure of which differ from known AAG binding agents. The 16-membered Au(2)P(4)C(8)O(2) macrocycles interconvert rapidly between two helical forms and produce enantiomeric conformations which are in dynamic equilibrium in solution. AAG binds preferentially one of the chiral conformers as indicated by strong Cotton effects generated by intramolecular exciton coupling between the pairs of hetercyclic chromophores. Lipophilic nature of the guest molecules suggests the dominant contribution of hydrophobic interactions in the AAG binding. Comparison of the main genetic variants of AAG revealed that both the 'F1/S' and 'A' variants bind with high-affinity the gold(I) macrocycles (K(a) approximately 10(6) M(-1)). CD/fluorescence displacement, and fluorescence quenching experiments indicated inclusion of the compounds into the central beta-barrel cavity of AAG of which exact tertiary structure is yet unknown. Molecular dimensions of the gold(I) macrocycles (13 x 14 x 14 A) indicate that the principal ligand binding cavity of both the 'F1/S' and 'A' variants must be larger compared to the models published to date. Based on these findings, a novel homology model of AAG 'F1' variant was constructed using the human neutrophil gelatinase-associated lipocalin as a template. The organogold complexes were successfully docked into the central cavity of this model.

A comprehensive study on the 5-hydroxytryptamine(3A) receptor binding of agonists serotonin and m-chlorophenylbiguanidine.

Serotonin type 3 receptors (5-HT(3)R) are members of the ligand gated ion channel receptor family. In this study, the interactions of the agonists serotonin (5-HT) and m-chlorophenylbiguanidine (mCPBG) at the binding site of the 5-HT(3A)R were investigated at an atomic level. Site-directed mutagenesis studies in Loop B and E along with our earlier published results from mutations within Loops A, C, and D provide comprehensive data on the interaction of 5-HT and mCPBG with 5-HT(3A)Rs. Using this data we have constructed a refined homology model of the 5-HT(3A)R that considers all of the available experimental data. 5-HT and mCPBG were docked into the newly constructed homology model and the amino acid residues critical in binding of these agonists were compared and analyzed. Our docking results reveal many similar binding interactions for 5-HT and mCPBG. Namely, residues THR181, TRP183, PHE226, ILE228, TYR234 and GLU129 were all found to play key roles in binding of both 5-HT and mCPBG. However, the results also revealed two important differences that exist between the interactions of the two agonists. In our model, a hydrogen bond is formed between the indole hydrogen of 5-HT and the residue TYR153. This interaction is not present in the case of mCPBG. Conversely, a hydrogen bond exists between SER182 and a protonated nitrogen of mCPBG, which does not exist in 5-HT. Our modeling results were found to be in accordance with experimental data.

[DockingServer: molecular docking calculations online]. / DockingServer: molekuláris dokkolási számítások a böngészoben.

Over the last years, the use of bioinformatics tools such as molecular docking has become an essential part of research focused at prediction of the binding of small molecules to their target proteins. DockingServer offers a web-based, easy to use interface that handles all aspects of molecular docking from ligand and pro-tein set-up through results representation integrating a number of software frequently used in computational chemistry. While its user friendly interface enables docking calculation and results evaluation carried out by researchers coming from all fields of biochemistry, DockingServer also provides full control on the setting of specific parameters of ligand and protein set up and docking calculations for more advanced users. The application can be used for docking and analysis of single ligands as well as for high throughput docking of ligand libraries to target proteins. The use of "DockingServer" is illustrated by the formation of acetaminophene (paracetamol)-CYP2E1 complex.

Homology modeling of breast cancer resistance protein (ABCG2).

BCRP (also known as ABCG2, MXR, and ABC-P) is a member of the ABC family that transports a wide variety of substrates. BCRP is known to play a key role as a xenobiotic transporter. Since discovering its role in multidrug resistance, considerable efforts have been made in order to gain deeper understanding of BCRP structure and function. The recent study was aimed at predicting BCRP structure by creating a homology model. Based on sequence similarity with known structures of full-length, NB and TM domain of ABC transporters, TM, NB, and linker regions of BCRP were defined. The NB domain of BCRP was modeled using MalK as a template. Based on secondary structure prediction of BCRP and comparison of the transmembrane connecting regions of known structures of ABC transporters, the TM domain arrangement of BCRP was established and was found to resemble to that of the recently published crystal structure of Sav1866. Thus, an initial alignment of TM domain of BCRP was established using Sav1866 as a template. This alignment was subsequently refined using constrains derived from secondary structure and TM predictions and the final model was built. Finally, the complete homodimer ABCG2 model was generated using Sav1866 as template. Furthermore, known ligands of BCRP were docked to our model in order to define possible binding sites. The results of molecular dockings of known BCRP substrates to the BCRP model were in agreement with recently published experimental data indicating multiple binding sites in BCRP.

Cyclodextrin complexes: chiral recognition and complexation behaviour.

Cyclodextrins are well known in supramolecular chemistry as host molecules capable of engulfing molecules in their hydrophobic cavity via noncovalent interactions. Although cyclodextrins are frequently used for chiral separation of racemates, the mechanism of chiral recognition has not yet been fully characterised. The current investigation was aimed at examining chiral recognition mechanism in order to construct an in silico method for prediction of chiral recognition. Amino acids were selected as model guest, whereas alphaCD was used as model host. The results of molecular docking and molecular dynamic calculations were compared to results of stability constant determination and capillary electrophoresis measurements of enantioseparations. Positive correlation between binding strength and chiral separation ability was found. However, the small energy differences between interaction energy of each enantiomer with alphaCD fell into the range of standard error of molecular docking calculations limiting its applicability for in silico prediction. Examining the stability of complex geometry during molecular dynamics simulation revealed that stable complex geometry is likely to be a prerequisite for chiral recognition. This hypothesis was tested on methylderivatized tryptophan. Indeed, chiral separation of beta-methyl-tryptophans by alphaCD could be successfully predicted by examining the complex geometries during molecular dynamic simulation.

Aggregation of cyclodextrins as an important factor to determine their complexation behavior.

Here, we report a study on the complexation behavior of carotenoids with cyclodextrins (CDs) using solubility experiments and molecular-modelling methods. Carotenoids are an important group of naturally occurring dyes found in vegetables and fruits. Their antioxidant property has initiated investigations on their possible use as drugs. However, carotenoids are lipophilic molecules with very little inherent aqueous solubility. Cyclodextrin complexation has been widely used in order to increase the potential applications of hydrophobic compounds. Thus, the aim of our investigation was to design carotenoids with enhanced water solubility by cyclodextrin complexation. Molecular modelling of carotenoid-cyclodextrin complexes with a 1 : 1 stoichiometry successfully explained the experimentally observed capability of beta-cyclodextrins (beta-CDs) to form complexes with carotenoids as opposed to alpha-cyclodextrins (alpha-CDs) and gamma-cyclodextrins (gamma-CDs). Furthermore, molecular-dynamics calculations revealed that the aggregation properties of CD derivatives significantly influence their complexation behavior. Our docking calculations showed that RAMEB (random methylated beta-CD) is the beta-CD derivative that possesses the lowest tendency to aggregate. Solubility experiments yielded the same results, namely, RAMEB complexes possess the best water solubility. Our results showed that complexation of a ligand not buried inside of the CD cavity is dependent on two factors: i) the geometry of the inclusion part of the complex; ii) the self-aggregation property of the CD itself. The lower affinity the CDs possess for self-aggregation, the more likely are they involved in interactions with carotenoids. These results suggest that self-aggregation of CDs should be considered as an important parameter determining complexation in general.

Binding of the pepper alkaloid piperine to bovine beta-lactoglobulin: circular dichroism spectroscopy and molecular modeling study.

The pepper alkaloid piperine is a nontoxic, natural dietary compound with a broad range of physiological activity. The present work is the first demonstration of its interaction with a mammalian protein. Circular dichroism (CD) spectroscopy was used to reveal and analyze the binding of piperine to a lipocalin protein. Induced CD spectra measured in pH 7.7 phosphate buffer at 37 degrees C demonstrated reversible, non-covalent association of piperine with bovine beta-lactoglobulin (BLG), the major whey protein in milk. The binding parameters (K(a) approximately 8 x 10(4) M(-1), n = 0.8) determined from the CD titration data showed no significant differences between the piperine binding properties of the two main genetic variants of BLG (A and B). The vanishing extrinsic CD signal obtained upon acidification of the piperine-BLG sample solution (Tanford transition) suggested that the ligand binds in the central hydrophobic cavity of the beta-barrel. The cavity binding concept was further supported by a CD displacement experiment using palmitic acid, the well-known hydrophobic ligand of BLG. Molecular docking calculations showed that piperine can be efficiently accommodated within the calyx of BLG. Additional molecular modeling calculations indicated that the beta-barrel of human tear lipocalin, human serum retinol binding protein, and human neutrophil gelatinase associated lipocalin might also accommodate a piperine molecule.

Inhibition of cytochrome P450 enzymes participating in p-nitrophenol hydroxylation by drugs known as CYP2E1 inhibitors.

p-Nitrophenol hydroxylation is widely used as a probe for microsomal CYP2E1. Several drugs are known as CYP2E1 inhibitors because of their capability to inhibit p-nitrophenol hydroxylation. Our results suggest further participation of CYP2A6 and CYP2C19 enzymes in p-nitrophenol hydroxylation. Moreover, CYP2A6 and CYP2C19 may be considered as the primary catalysts, whereas CYP2E1 can also contribute to the hydroxylation of p-nitrophenol. Further aim of our study was to evaluate the selectivity of p-nitrophenol hydroxylase inhibitors towards cytochrome P450 enzymes. The effects of antifungals: bifonazole, econazole, clotrimazole, ketoconazole, miconazole; CNS-active drugs: chlorpromazine, desipramine, fluphenazine, thioridazine; and the non-steroidal anti-inflammatory drug: diclofenac were investigated on the enzyme activities selective for CYP2A6, CYP2C9, CYP2C19, CYP2E1 and CYP3A4. None of the drugs could be considered as a potent inhibitor of CYP2E1. Strong inhibition was observed for CYP3A4 by antifungals with IC(50) values in submicromolar range. However, ketoconazole was the only imidazole derivative that could be considered as a selective inhibitor of CYP3A4. The CNS-active drugs investigated were found to be weak inhibitors of CYP2A6, CYP2C9, CYP2C19, CYP2E1 and CYP3A4. Diclofenac efficiently inhibited CYP2C9 and to a less extent CYP3A4 enzyme.

A detailed study of antibacterial 3-acyltetramic acids and 3-acylpiperidine-2,4-diones.

Inspired by the core fragment of antibacterial natural products such as streptolydigin, 3-acyltetramic acids and 3-acylpiperidine-2,4-diones have been synthesised from the core heterocycle by direct acylation with the substituted carboxylic acids using a strategy which permits ready access to a structurally diverse compound library. The antibacterial activity of these systems has been established against a panel of Gram-positive and Gram-negative bacteria, with activity mostly against the former, which in some cases is very potent. Data consistent with modes of action against undecaprenylpyrophosphate synthase (UPPS) and/or RNA polymerase (RNAP) for a small subset of the library has been obtained. The most active compounds have been shown to exhibit binding at known binding sites of streptolydigin and myxopyronin at UPPS and RNAP. These systems offer potential for their antibacterial activity, and further demonstrate the use of natural products as biologically validated starting points for drug discovery.

Syntheses and antiproliferative effects of D-homo- and D-secoestrones.

Substituted and/or heterocyclic d-homoestrone derivatives were synthetized via the intramolecular cyclization of a δ-alkenyl-d-secoaldehyde, -d-secoalcohol or -d-secocarboxylic acid of estrone 3-benzyl ether. The d-secoalcohol was modified at three sites in the molecule. The in vitro antiproliferative activities of the new d-homo- and d-secoestrone derivatives were determined on HeLa, MCF-7, A431 and A2780 cells through use of MTT assay. d-Homoalcohols 3 and 5 displayed cell line-selective cytostatic effects against ovarian and cervical cell lines, respectively. Two d-secoestrones (6 and 12c) proved to be effective, with IC50 values comparable with those of the reference agent cisplatin. A selected compound (6) was tested by tubulin polymerization assay and its cancer specificity was additionally determined by using noncancerous human fibroblast cells.

Disordered proteins and network disorder in network descriptions of protein structure, dynamics and function: hypotheses and a comprehensive review.

During the last decade, network approaches became a powerful tool to describe protein structure and dynamics. Here we review the links between disordered proteins and the associated networks, and describe the consequences of local, mesoscopic and global network disorder on changes in protein structure and dynamics. We introduce a new classification of protein networks into 'cumulus-type', i.e., those similar to puffy (white) clouds, and 'stratus-type', i.e., those similar to flat, dense (dark) low-lying clouds, and relate these network types to protein disorder dynamics and to differences in energy transmission processes. In the first class, there is limited overlap between the modules, which implies higher rigidity of the individual units; there the conformational changes can be described by an 'energy transfer' mechanism. In the second class, the topology presents a compact structure with significant overlap between the modules; there the conformational changes can be described by 'multi-trajectories'; that is, multiple highly populated pathways. We further propose that disordered protein regions evolved to help other protein segments reach 'rarely visited' but functionally-related states. We also show the role of disorder in 'spatial games' of amino acids; highlight the effects of intrinsically disordered proteins (IDPs) on cellular networks and list some possible studies linking protein disorder and protein structure networks.