Homology model of the human tRNA splicing ligase RtcB

RtcB is an essential human tRNA ligase required for ligating the 2',3'-cyclic phosphate and 5'-hydroxyl termini of cleaved tRNA halves during tRNA splicing and XBP1 fragments during endoplasmic reticulum stress. Activation of XBP1 has been implicated in various human tumors including breast cancer. Here we present, for the first time, a homology model of human RtcB (hRtcB) in complex with manganese and covalently bound GMP built from the Pyrococcus horikoshii RtcB (bRtcB) crystal structure, PDB ID 4DWQA. The structure is analyzed in terms of stereochemical quality, folding reliability, secondary structure similarity with bRtcB, druggability of the active site binding pocket and its metal-binding microenvironment. In comparison with bRtcB, loss of a manganese-coordinating water and movement of Asn226 (Asn202 in 4DWQA) to form metal-ligand coordination, demonstrates the uniqueness of the hRtcB model. Rotation of GMP leads to the formation of an additional metal-ligand coordination (Mn-O). Umbrella sampling simulations of Mn binding in wild type and the catalytically inactive C122A mutant reveal a clear reduction of Mn binding ability in the mutant, thus explaining the loss of activity therein. Our results furthermore clearly show that the GTP binding site of the enzyme is a well-defined pocket that can be utilized as target site for in silico drug discovery.

Analysis of aquaporins from the euryhaline barnacle Balanus improvisus reveals differential expression in response to changes in salinity

Barnacles are sessile macro-invertebrates, found along rocky shores in coastal areas worldwide. The euryhaline bay barnacle Balanus improvisus (Darwin, 1854) (= Amphibalanus improvisus) can tolerate a wide range of salinities, but the molecular mechanisms underlying the osmoregulatory capacity of this truly brackish species are not well understood. Aquaporins are pore-forming integral membrane proteins that facilitate transport of water, small solutes and ions through cellular membranes, and that have been shown to be important for osmoregulation in many organisms. The knowledge of the function of aquaporins in crustaceans is, however, limited and nothing is known about them in barnacles. We here present the repertoire of aquaporins from a thecostracan crustacean, the barnacle B. improvisus, based on genome and transcriptome sequencing. Our analyses reveal that B. improvisus contains eight genes for aquaporins. Phylogenetic analysis showed that they represented members of the classical water aquaporins (Aqp1, Aqp2), the aquaglyceroporins (Glp1, Glp2), the unorthodox aquaporin (Aqp12) and the arthropod-specific big brain aquaporin (Bib). Interestingly, we also found two big brain-like proteins (BibL1 and BibL2) constituting a new group of aquaporins not yet described in arthropods. In addition, we found that the two water-specific aquaporins were expressed as C-terminal splice variants. Heterologous expression of some of the aquaporins followed by functional characterization showed that Aqp1 transported water and Glp2 water and glycerol, agreeing with the predictions of substrate specificity based on 3D modeling and phylogeny. To investigate a possible role for the B. improvisus aquaporins in osmoregulation, mRNA expression changes in adult barnacles were analysed after long-term acclimation to different salinities. The most pronounced expression difference was seen for AQP1 with a substantial (>100-fold) decrease in the mantle tissue in low salinity (3 PSU) compared to high salinity (33 PSU). Our study provides a base for future mechanistic studies on the role of aquaporins in osmoregulation.

Analysis of Biphenyl-Type Inhibitors Targeting the Eg5 α4/α6 Allosteric Pocket

Eg5 is a mitotic kinesin protein that plays an important role in the formation and maintenance of the bipolar spindle during the mitotic phase. Due to its potentially reduced side effects in cancer therapy, Eg5 is considered to be an attractive target for developing anticancer inhibitors. Herein, we report a computational modeling study involving biphenyl-type inhibitors known to interact with the α4/α6 allosteric pocket of Eg5. Compared to the well-known α2/L5/α3 allosteric inhibitors, biphenyl-type inhibitors show a unique activity profile. In the Eg5-PVZB1194 (a biphenyl-type inhibitor) crystal structure, loop L11, which is located in the entrance of the α4/α6 allosteric-binding pocket, is missing due to crystal-packing effects. To better understand the role of this flexible loop upon biphenyl-type inhibitor-binding, MD simulations were performed to observe the L11 conformations from different states. It was demonstrated that L11 was more stabilized and showed less fluctuation when PVZB1194 was bound to Eg5. Residue Asn287 from L11 forms hydrogen bonding to the sulfone group of PVZB1194, whereby L11 moves inward to the α4/α6 allosteric pocket and moves away from the pocket in absence of the inhibitor. Pharmacophore, three-dimensional (3D)-QSAR, and ADME studies of biphenyl-type inhibitors of Eg5 were also performed. A best pharmacophore model, DDRRH.6, was generated, having correlation coefficients in the 3D-QSAR study of R² = 0.81 and Q² = 0.64. Furthermore, docking studies were carried out to observe the interaction between the remaining biphenyl-type inhibitors with Eg5. In addition, on the basis of fragment docking, a structure-based pharmacophore was generated, which shares good overlap of the DHRR features of the pharmacophore model DDHRR.6. The structure-based pharmacophore also contains extra hydrogen-bond acceptors and hydrophobic groups, features which provide possibilities in developing new or improved series of compounds.

Ligand Selectivity between the ADP-Ribosylating Toxins: An Inverse-Docking Study for Multitarget Drug Discovery

Bacterial adenosine 5'-diphosphate-ribosylating toxins are encoded by several human pathogens, such as Pseudomonas aeruginosa (exotoxin A (ETA)), Corynebacterium diphtheriae (diphtheria toxin (DT)), and Vibrio cholerae (cholix toxin (CT)). The toxins modify eukaryotic elongation factor 2, an essential human enzyme in protein synthesis, thereby causing cell death. Targeting external virulence factors, such as the above toxins, is a promising alternative for developing new antibiotics, while at the same time avoiding drug resistance. This study aims to establish a reliable computational methodology to find a "silver bullet" able to target all three toxins. Herein, we have undertaken a detailed analysis of the active sites of ETA, DT, and CT, followed by the determination of the most appropriate selection of the size of the docking sphere. Thereafter, we tested two different approaches for normalizing the docking scores and used these to verify the best target (toxin) for each ligand. The results indicate that the methodology is suitable for identifying selective as well as multitoxin inhibitors, further validating the robustness of inverse docking for target-fishing experiments.

Defects in the calcium-binding region drastically affect the cadherin-like domains of RET tyrosine kinase

Mutations in the rearranged during transfection (RET) tyrosine kinase gene leading to gain or loss of function have been associated with the development of several human cancers and Hirschsprung's disease (HSCR). However, to what extent these mutations affect individual bio-molecular functions remains unclear. In this article, the functionally significant mutations in the RET CLD1-4 calcium-binding site which lead to HSCR, and depletion of calcium ions in the RET CLD1-4 calcium binding site, were investigated by molecular dynamics simulations--to understand the mechanistic action of the mutations or loss of calcium ions in altering the protein kinase structure, dynamics, and stability. The mutations or loss of calcium ions change the local conformation and change the free energy landscape. Specifically, the mutations and loss of calcium ions decrease the radius of gyration of the whole structure, leading to improper protein folding and GFL-GFRα contact site reduction. Furthermore, based on the most populated conformation in the wildtype MD simulations, a pharmacophore was generated by fragment docking to identify key features of the possible inhibitors targeting the calcium binding site. Overall, the findings may provide useful structural insights into the molecular mechanism underlying RET calcium-binding site mutations and assist in development of novel drugs targeting the extracellular ligand contact site of wildtype RET.

Estimation of Liposome Penetration Barriers of Drug Molecules with All-Atom and Coarse-Grained Models

Liposomes are common carriers of drug molecules, providing enhanced delivery and accumulation of hydrophilic agents or larger biomolecules. Molecular simulations can be used to estimate key features of the drug molecules upon interaction with the liposomes, such as penetration barriers and localization. Herein, we investigate several aspects of the computational estimation of penetration barriers, viz. the potential of mean force (PMFs) along a vector spanning the membrane. First, we provide an evaluation of the all-atom (AA) and coarse-grained (CG) parametrization of 5-aminolevulinic acid (5-ALA) and two of its alkyl esters by computing n-octanol/water partition coefficients. We find that the CG parametrization of the esters performs significantly better than the CG model of 5-ALA, highlighting the difficulty to coarse-grain small, polar molecules. However, the expected trend in partition coefficients is reproduced also with the CG models. Second, we compare PMFs in a small membrane slab described with either the AA or CG models. Here, we are able to reproduce the all-atom PMF of 5-ALA with CG. However, for the alkyl esters it is unfortunately not possible to correctly reproduce both the depth and the penetration barrier of the PMF seen in the AA simulations with any of the tested CG models. We argue that it is more important to choose a CG parametrization that reproduces the depth of the PMF. Third, we compare, using the CG model, PMFs in the membrane slab with PMFs in a large, realistic liposome. We find similar depths but slightly different penetration barriers most likely due to differences in the lipid density along the membrane axis. Finally, we compute PMFs in liposomes with three different lipid compositions. Unfortunately, differences in the PMFs could not be quantified, and it remains to be investigated to what extent liposome simulations can fully reproduce experimental findings.

Structures, Properties, and Dynamics of Intermediates in eEF2-Diphthamide Biosynthesis

The eukaryotic translation Elongation Factor 2 (eEF2) is an essential enzyme in protein synthesis. eEF2 contains a unique modification of a histidine (His699 in yeast; HIS) into diphthamide (DTA), obtained via 3-amino-3-carboxypropyl (ACP) and diphthine (DTI) intermediates in the biosynthetic pathway. This essential and unique modification is also vulnerable, in that it can be efficiently targeted by NAD(+)-dependent ADP-ribosylase toxins, such as diphtheria toxin (DT). However, none of the intermediates in the biosynthesis path is equally vulnerable against the toxins. This study aims to address the different susceptibility of DTA and its precursors against bacterial toxins. We have herein undertaken a detailed in silico study of the structural features and dynamic motion of different His699 intermediates along the diphthamide synthesis pathway (HIS, ACP, DTI, DTA). The study points out that DTA forms a strong hydrogen bond with an asparagine which might explain the ADP-ribosylation mechanism caused by the diphtheria toxin (DT). Finally, in silico mutagenesis studies were performed on the DTA modified protein, in order to hamper the formation of such a hydrogen bond. The results indicate that the mutant structure might in fact be less susceptible to attack by DT and thereby behave similarly to DTI in this respect.

Secretion of protein disulphide isomerase AGR2 confers tumorigenic properties

The extracellular matrix (ECM) plays an instrumental role in determining the spatial orientation of epithelial polarity and the formation of lumens in glandular tissues during morphogenesis. Here, we show that the Endoplasmic Reticulum (ER)-resident protein anterior gradient-2 (AGR2), a soluble protein-disulfide isomerase involved in ER protein folding and quality control, is secreted and interacts with the ECM. Extracellular AGR2 (eAGR2) is a microenvironmental regulator of epithelial tissue architecture, which plays a role in the preneoplastic phenotype and contributes to epithelial tumorigenicity. Indeed, eAGR2, is secreted as a functionally active protein independently of its thioredoxin-like domain (CXXS) and of its ER-retention domain (KTEL), and is sufficient, by itself, to promote the acquisition of invasive and metastatic features. Therefore, we conclude that eAGR2 plays an extracellular role independent of its ER function and we elucidate this gain-of-function as a novel and unexpected critical ECM microenvironmental pro-oncogenic regulator of epithelial morphogenesis and tumorigenesis.

DOI:http://dx.doi.org/10.7554/eLife.13887.001

Cancer cells multiply abnormally fast and therefore produce protein molecules faster than normal cells. To avoid becoming stressed by this overproduction, cancer cells make use of proteins that fold the new proteins inside the cell. One of these protein folders is called anterior gradient-2 (or AGR2 for short) and is produced at high levels in so-called epithelial cancers, such as breast and lung cancer.

Previous research has shown that AGR2 inside cancer cells can help them grow and survive and AGR2 can also be found outside cells, such as in the blood or the urine of cancer patients. Therefore some researchers have suggested that measuring the levels of AGR2 in bodily fluids may be a useful marker for detecting cancers. Fessart et al. hypothesized that – apart from becoming a promising diagnostic tool – the AGR2 protein itself, specifically when found outside cells, might make cancer cells more aggressive.

Fessart et al. used a range of techniques to test this hypothesis. For example, healthy lung cells and lung cancer cells were grown into miniature replicas of lung organs in the laboratory, and in a key experiment, AGR2 was added to the lung organoids grown from the healthy cells. The addition of AGR2 protein was enough to change the non-tumor organoids into tumor organoids and boosted their growth about ten-fold. Further experiments then revealed that AGR2 also makes cells more invasive and capable of moving, both important features of aggressive cancer cells.

Overall, Fessart et al. have proven that AGR2 is a signalling molecule found outside cancer cells that makes them more aggressive. In future, more research addressing how AGR2 achieves this may lead to new therapeutic strategies against some forms of cancer.

DOI:http://dx.doi.org/10.7554/eLife.13887.002

DNA Distortion Caused by Uracil-Containing Intrastrand Cross-Links

Four uracil-containing intrastrand cross-links have been detected in human cells upon UV irradiation of 5-bromouracil-containing DNA, namely 5'-G[8-5]U-3', 5'-U[5-8]G-3', 5'-A[8-5]U-3', and 5'-A[2-5]U-3'. These lesions feature unique composition and connectivity compared with other intrastrand cross-links reported in the literature. For the first time, structural information obtained using molecular dynamics (MD) simulations reveal that all four lesions distort the DNA helix, which can involve an extrahelical location of the cross-link, changes in the helical interactions of the complementary nucleotides, or disruption of hydrogen bonding in the flanking base pairs up to two positions from the cross-linked site; however, the degree of distortion varies between the cross-links, being affected by the sequence, nucleobase-nucleobase connectivity, and the purine involved. Most importantly, the relative distortion of the damaged DNA provides the first structural explanation for the observed abundances of the four uracil-containing cross-links. Furthermore, the highly distorted conformations suggest that these lesions will likely have severe implications for DNA replication and repair processes in cells.

Molecular Dynamics Studies of Liposomes as Carriers for Photosensitizing Drugs: Development, Validation, and Simulations with a Coarse-Grained Model

Liposomes are proposed as drug delivery systems and can in principle be designed so as to cohere with specific tissue types or local environments. However, little detail is known about the exact mechanisms for drug delivery and the distributions of drug molecules inside the lipid carrier. In the current work, a coarse-grained (CG) liposome model is developed, consisting of over 2500 lipids, with varying degrees of drug loading. For the drug molecule, we chose hypericin, a natural compound proposed for use in photodynamic therapy, for which a CG model was derived and benchmarked against corresponding atomistic membrane bilayer model simulations. Liposomes with 21-84 hypericin molecules were generated and subjected to 10 microsecond simulations. Distribution of the hypericins, their orientations within the lipid bilayer, and the potential of mean force for transferring a hypericin molecule from the interior aqueous "droplet" through the liposome bilayer are reported herein.

Sphingolipids contribute to acetic acid resistance in Zygosaccharomyces bailii

Lignocellulosic raw material plays a crucial role in the development of sustainable processes for the production of fuels and chemicals. Weak acids such as acetic acid and formic acid are troublesome inhibitors restricting efficient microbial conversion of the biomass to desired products. To improve our understanding of weak acid inhibition and to identify engineering strategies to reduce acetic acid toxicity, the highly acetic‐acid‐tolerant yeast Zygosaccharomyces bailii was studied. The impact of acetic acid membrane permeability on acetic acid tolerance in Z. bailii was investigated with particular focus on how the previously demonstrated high sphingolipid content in the plasma membrane influences acetic acid tolerance and membrane permeability. Through molecular dynamics simulations, we concluded that membranes with a high content of sphingolipids are thicker and more dense, increasing the free energy barrier for the permeation of acetic acid through the membrane. Z. bailii cultured with the drug myriocin, known to decrease cellular sphingo­lipid levels, exhibited significant growth inhibition in the presence of acetic acid, while growth in medium without acetic acid was unaffected by the myriocin addition. Furthermore, following an acetic acid pulse, the intracellular pH decreased more in myriocin‐treated cells than in control cells. This indicates a higher inflow rate of acetic acid and confirms that the reduction in growth of cells cultured with myriocin in the medium with acetic acid was due to an increase in membrane permeability, thereby demonstrating the importance of a high fraction of sphingolipids in the membrane of Z. bailii to facilitate acetic acid resistance; a property potentially transferable to desired production organisms suffering from weak acid stress. Biotechnol. Bioeng. 2016;113: 744–753. © 2015 The Authors. Biotechnology and Bioengineering Published by Wiley Periodicals, Inc.

Homology Models and Molecular Modeling of Human Retinoic Acid Metabolizing Enzymes Cytochrome P450 26A1 (CYP26A1) and P450 26B1 (CYP26B1)

Homology models of cytochrome P450 26A1 and cytochrome P450 26B1 were constructed using the crystal structures of human, CYP2C8, CYP2C9, and CYP3A4 as templates for the model building. The homology models generated were investigated for their docking capacities against the natural substrate all-trans-retinoic acid (atRA), five different tetralone-derived retinoic acid metabolizing blocking agents (RAMBAs), and R115866. Interaction energies (IE) and linear interaction energies (LIE) were calculated for all inhibitors in both homology models after molecular dynamics (MD) simulation of the enzyme-ligand complexes. The results revealed that the homologues had the capacity to distinguish between strong and weak inhibitors. Important residues in the active site were identified from the CYP26A1/B1-atRA complexes. Residues involved in hydrophobic interactions with atRA were Pro113, Phe222, Phe299, Val370, Pro371, and Phe374 in CYP26A1 and Leu88, Pro118, Phe222, Phe295, Ile368, and Tyr272 in CYP26B1. Hydrogen bonding interactions were observed between the atRA carboxylate group and Arg 90 in CYP26A1 and with Arg76, Arg95, and Ser369 in CYP26B1.

Characterization of interactions and pharmacophore development for DFG-out inhibitors to RET tyrosine kinase

RET (rearranged during transfection) tyrosine kinase is a promising target for several human cancers. Abt-348, Birb-796, Motesanib and Sorafenib are DFG-out multi-kinase inhibitors that have been reported to inhibit RET activity with good IC50 values. Although the DFG-out conformation has attracted great interest in the design of type II inhibitors, the structural requirements for binding to the RET DFG-out conformation remains unclear. Herein, the DFG-out conformation of RET was determined by homology modelling, the four inhibitors were docked, and the binding modes investigated by molecular dynamics simulation. Binding free energies were calculated using the molecular mechanics/Poisson-Bolzmann surface area (MM/PBSA) method. The trends in predicted binding free affinities correlated well with experimental data and were used to explain the activity difference of the studied inhibitors. Per-residue energy decomposition analyses provided further information on specific interaction properties. Finally, we also conducted a detailed e-pharmacophore modelling of the different RET-inhibitor complexes, explaining the common and specific pharmacophore features of the different complexes. The results reported herein will be useful in future rational design of novel DFG-out RET inhibitors.

De novo tertiary structure prediction using RNA123--benchmarking and application to Macugen

The present benchmarking study utilizes the RNA123 program for de novo prediction of tertiary structures of a set of 50 RNA molecules for which X-ray/NMR structures are available, based on the nucleic acid sequence only. All molecules contain a hairpin loop motif and a helical structure of canonical and non-canonical base pairs, interrupted by bulges and internal loops to various degrees. RNA molecules with double helices made up purely by canonical base pairing, and molecules containing symmetric internal loops of non-canonical base pairing are, overall, very well predicted. Structures containing bulges and asymmetric internal loops, and more complex structures containing multiple bulges and internal loops in the same molecule, result in larger deviations from their X-ray/NMR predicted structures due to higher degree of flexibility of the nucleotide bases in these regions. In a majority of the molecules included herein, the RNA123 program was, however, able to predict the tertiary structure with a heavy atom RMSD of less than 5 Å to the X-ray/NMR structure, and the models were in most cases structurally closer to the X-ray/NMR structures than models predicted by MC-Fold and MC-Sym. A set of RNA molecules containing pseudoknot tertiary structure motifs were included, but neither of the programs was able to predict the folding of the single-stranded stem onto the helix without additional structural input. The RNA123 program was then applied to predict the tertiary structure of the RNA segment of Macugen®, the first RNA aptamer approved for clinical use, and for which no tertiary structure has yet been solved. Four possible tertiary structures were predicted for this 27-nucleic-acid-long RNA molecule, which will be used in constructing a full model of the PEGylated aptamer and its interaction with the vascular endothelial growth factor target.

Rational design and validation of a Tip60 histone acetyltransferase inhibitor

Histone acetylation is required for many aspects of gene regulation, genome maintenance and metabolism and dysfunctional acetylation is implicated in numerous diseases, including cancer. Acetylation is regulated by histone acetyltransferases (HATs) and histone deacetylases and currently, few general HAT inhibitors have been described. We identified the HAT Tip60 as an excellent candidate for targeted drug development, as Tip60 is a key mediator of the DNA damage response and transcriptional co-activator. Our modeling of Tip60 indicated that the active binding pocket possesses opposite charges at each end, with the positive charges attributed to two specific side chains. We used structure based drug design to develop a novel Tip60 inhibitor, TH1834, to fit this specific pocket. We demonstrate that TH1834 significantly inhibits Tip60 activity in vitro and treating cells with TH1834 results in apoptosis and increased unrepaired DNA damage (following ionizing radiation treatment) in breast cancer but not control cell lines. Furthermore, TH1834 did not affect the activity of related HAT MOF, as indicated by H4K16Ac, demonstrating specificity. The modeling and validation of the small molecule inhibitor TH1834 represents a first step towards developing additional specific, targeted inhibitors of Tip60 that may lead to further improvements in the treatment of breast cancer.

Deamination features of 5-hydroxymethylcytosine, a radical and enzymatic DNA oxidation product

The 5-methylcytosine derivative 5-hydroxymethylcytosine (5hmCyt), which is generated via enzymatic oxidation, is sometimes referred to as the sixth nucleobase due to its widespread presence in the DNA of brain and embryonic stem cells. In this study, we used density functional based methods and reactivity indices from conceptual DFT to explore the mechanism and key features of the hydrolytic deamination of 5hmCyt. The data obtained are used to compare and contrast this deamination reaction with those of other cytosine derivatives. The deamination process for 5hmCyt is similar to the corresponding processes for other unsaturated derivatives in that the amino form is the reactive one and water addition is the rate-limiting step. However, several differences due to the rotameric asymmetry of the current system are also noted.

Hydrogen peroxide contributes to the ultraviolet-B (280-315 nm) induced oxidative stress of plant leaves through multiple pathways

Solar UV-B (280-315 nm) radiation is a developmental signal in plants but may also cause oxidative stress when combined with other environmental factors. Using computer modeling and in solution experiments we show that UV-B is capable of photosensitizing hydroxyl radical production from hydrogen peroxide. We present evidence that the oxidative effect of UV-B in leaves is at least twofold: (i) it increases cellular hydrogen peroxide concentrations, to a larger extent in pyridoxine antioxidant mutant pdx1.3-1 Arabidopsis and; (ii) is capable of a partial photo-conversion of both 'natural' and 'extra' hydrogen peroxide to hydroxyl radicals. As stress conditions other than UV can increase cellular hydrogen peroxide levels, synergistic deleterious effects of various stresses may be expected already under ambient solar UV-B.

Photochemical reaction mechanism of UV-B-induced monomerization of UVR8 dimers as the first signaling event in UV-B-regulated gene expression in plants

The Arabidopsis thaliana UV RESISTANCE LOCUS8 (UVR8) protein has been identified to specifically mediate photomorphogenic UV-B responses by acting as a UV-B photoreceptor. The dimeric structure of the UVR8 protein dissociates into signaling-active monomers upon UV-B exposure, and the monomers rapidly interact with downstream signaling components to regulate gene expression. UVR8 monomers revert to dimers in the absence of UV-B radiation, thereby reversing transcription activation. UVR8 amino acid residues W233 and W285 have been identified to play critical roles in the UVR8 dimer for the response to UV-B irradiation. In the present work, the photoreaction mechanism for UVR8 monomerization is explored with quantum chemical cluster calculations and evaluated by molecular dynamics simulations using the wild-type UVR8 dimer and novel force field parameters developed for intermediate radicals formed in the photochemical process. Three different models are investigated, which show that the preferred mechanism for UVR8 monomerization involves electron transfer from residue W233 to W285 and onward to R338 initiated by UV-B irradiation, coupled to simultaneous proton transfer from W233 to D129 leading to the formation of protonated D129, a deprotonated W233 radical, and a neutral R338 radical. Due to the formation of the neutral R338 radical, salt bridges involving this residue are disrupted together with the concomitant interruption of several other key salt bridges R286-D96, R286-D107, R338-D44, R354-E43, and R354-E53. The resulting large decrease in protein-protein interaction energy arising from this sequence of events leads to the monomerization of the UVR8 dimer. The mechanism presented is in accord with all experimental data available to date.

Modification of the anticancer drug tamoxifen to avoid CYP2D6 polymorphism

The prodrug tamoxifen (TAM) is the most widely used drug to treat breast cancer, and is metabolised to the active 4-hydroxy derivatives dominantly by hepatic CYP2D6. However, the application to patients with different polymorphic CYP2D6 has been under debate, because the efficacy of TAM is suspected to be suppressed in patients who have diminished CYP2D6 activity, resulting in inadequate active metabolites. We here propose modified structures, such as 4-methylTAM, which is highly possible to be activated by CYP3A, the most abundant CYP isoforms in the liver, whereby the genetic polymorphism of CYP2D6 is avoided. The diversity of CYP catalyzed metabolic paths for TAM and its derivatives are studied by quantum chemistry calculations on the reaction energies of the initial H atom abstraction steps. The ability of forming DNA adducts is compared through the formation enthalpy of the carbocation intermediate. The results suggest that the modified structures are safe with regard to forming DNA adducts and may be used as prodrugs in a wide range of patients, due to CYP3A, rather than CYP2D6, mediated activation.

Predicting enzyme-substrate specificity with QM/MM methods: a case study of the stereospecificity of (D)-glucarate dehydratase

The stereospecificity of d-glucarate dehydratase (GlucD) is explored by QM/MM calculations. Both the substrate binding and the chemical steps of GlucD contribute to substrate specificity. Although the identification of transition states remains computationally intensive, we suggest that QM/MM computations on ground states or intermediates can capture aspects of specificity that cannot be obtained using docking or molecular mechanics methods.

Catalysts or initiators? Beckmann rearrangement revisited

The catalytic mechanism of the organo-mediated Beckmann rearrangement has been modeled using DFT calculations. Five representative promoters were shown to be initiators rather than catalysts. A self-propagating mechanism is shown to be energetically much more favored than the previously proposed mechanisms involving a Meisenheimer complex.

Of mice and men: Dissecting the interaction between Listeria monocytogenes Internalin A and E-cadherin

We report a study of the interaction between internalin A (inlA) and human or murine E-cadherin (Ecad). inlA is used by Listeria monocytogenes to internalize itself into host cell, but the bacterium is unable to invade murine cells, which has been attributed to the difference in sequence between hEcad and mEcad. Using molecular dynamics simulations, MM/GBSA free energy calculations, hydrogen bond analysis, water characterization and umbrella sampling, we provide a complete atomistic picture of the binding between inlA and Ecad. We dissect key residues in the protein–protein interface and analyze the energetics using MM/GBSA. From this analysis it is clear that the binding of inlA–mEcad is weaker than inlA–hEcad, on par with the experimentally observed inability of inlA to bind to mEcad. However, extended MD simulations of 200 ns in length show no destabilization of the inlA–mEcad complex and the estimation of the potential of mean force (PMF) using umbrella sampling corroborates this conclusion. The binding strength computed from the PMFs show no significant difference between the two protein complexes. Hence, our study suggests that the inability of L. monocytogenes to invade murine cells cannot be explained by processes at the nanosecond to sub-microsecond time scale probed by the simulations performed here.

Catalytic mechanism of porphobilinogen synthase: the chemical step revisited by QM/MM calculations

Porphobilinogen synthase (PBGS) catalyzes the asymmetric condensation and cyclization of two 5-aminolevulinic acid (5-ALA) substrate molecules to give porphobilinogen (PBG). The chemical step of PBGS is herein revisited using QM/MM (ONIOM) calculations. Two different protonation states and several different mechanisms are considered. Previous mechanisms based on DFT-only calculations are shown unlikely to occur. According to these new calculations, the deprotonation step rather than ring closure is rate-limiting. Both the C-C bond formation first mechanism and the C-N bond formation first mechanism are possible, depending on how the A-site ALA binds to the enzyme. We furthermore propose that future work should focus on the substrate binding step rather than the enzymatic mechanism.

Identifying the sarco(endo)plasmic reticulum Ca2+ ATPase (SERCA) as a potential target for hypericin--a theoretical study

The exact cellular target for the potent anti-cancer agent hypericin has not yet been determined; this thus encourages the application of computational chemistry tools to be employed in order to provide insights that can be employed in further drug development studies. In the present study computational docking and molecular dynamics simulations are applied to investigate possible interactions between hypericin and the Ca(2+) pump SERCA as proposed in the literature. Hypericin was found to bind strongly both in pockets within the transmembrane region and in the cytosolic region of the protein, although the two studied isoforms of SERCA differ slightly in their preferred binding sites. The calculated binding energies for hypericin in the four investigated sites were of the same magnitude as for thapsigargin (TG), the most potent SERCA inhibitor, or in the range between TG and di-tert-butylhydroquinone (BHQ), which is also known to possess inhibitory activity. The hydrophobic character of hypericin indicates that the molecule initially binds in the ER membrane from which it diffuses into the transmembrane region of the protein and to binding pockets therein. The transmembrane TG and BHQ binding pockets provide suitable locations for hypericin as they allow for favourable interactions with the lipid tails that surround these. High binding energies were noted for hypericin in these pockets and are expected to constitute highly possible binding sites due to their accessibility from the ER membrane. Hypericin most likely binds to both isoforms of SERCA and acts as an inhibitor or, under light irradiation, as a singlet oxygen generator that in turn degrades the protein or induces lipid peroxidation.

Stability and iron coordination in DNA adducts of Anthracycline based anti-cancer drugs

There is evidence that the interaction of the α-ketol group of the Doxorubicin and Epirubicin anti-cancer drugs with Fe(III) generates hydroxyl radicals under aerobic conditions, causing cardiotoxicity in patients. Considering that the formation of DNA adducts is one of the main targets of Anthracycline drugs, we have in the present study characterized several [Anthracycline-DNA]Fe(III) complexes with respect to their stability and Fe(III) coordination, by means of MD simulations. Iron is found to coordinate well to the drugs containing an α-ketol group, this being the only group of the drug that binds to the metal. The complexes containing an α-ketol group, [Doxorubicin-DNA]Fe(III) and [Epirubicin-DNA]Fe(III), thus show greater stability than those not containing it, i.e., [Daunorubicin-DNA]Fe(III), [Idarubicin-DNA]Fe(III) and [5-Imino-Daunorubicin]Fe(III). Metal attachment to the α-ketol group is furthermore facilitated by the phosphate groups of DNA. The coordination to iron in the [Doxorubicin-DNA]Fe(III) system is smaller than that found for the [Epirubicin-DNA]Fe(III) system, and the corresponding number of coordinating waters in the former is larger than in the latter. This may in turn result in higher hydroxyl radical production, thus explaining the increased cardiotoxicity noted for Doxorubicin.

Cloning and functional studies of a splice variant of CYP26B1 expressed in vascular cells

Background

All-trans retinoic acid (atRA) plays an essential role in the regulation of gene expression, cell growth and differentiation and is also important for normal cardiovascular development but may in turn be involved in cardiovascular diseases, i.e. atherosclerosis and restenosis. The cellular atRA levels are under strict control involving several cytochromes P450 isoforms (CYPs). CYP26 may be the most important regulator of atRA catabolism in vascular cells. The present study describes the molecular cloning, characterization and function of atRA-induced expression of a spliced variant of the CYP26B1 gene.

Methodology/Principal Findings

The coding region of the spliced CYP26B1 lacking exon 2 was amplified from cDNA synthesized from atRA-treated human aortic smooth muscle cells and sequenced. Both the spliced variant and full length CYP26B1 was found to be expressed in cultured human endothelial and smooth muscle cells, and in normal and atherosclerotic vessel. atRA induced both variants of CYP26B1 in cultured vascular cells. Furthermore, the levels of spliced mRNA transcript were 4.5 times higher in the atherosclerotic lesion compared to normal arteries and the expression in the lesions was increased 20-fold upon atRA treatment. The spliced CYP26B1 still has the capability to degrade atRA, but at an initial rate one-third that of the corresponding full length enzyme. Transfection of COS-1 and THP-1 cells with the CYP26B1 spliced variant indicated either an increase or a decrease in the catabolism of atRA, probably depending on the expression of other atRA catabolizing enzymes in the cells.

Conclusions/Significance

Vascular cells express the spliced variant of CYP26B1 lacking exon 2 and it is also increased in atherosclerotic lesions. The spliced variant displays a slower and reduced degradation of atRA as compared to the full-length enzyme. Further studies are needed, however, to clarify the substrate specificity and role of the CYP26B1 splice variant in health and disease.

Characterization of agonist binding to His524 in the estrogen receptor α ligand binding domain

The bioactivities of the natural steroidal estrogen 17β-estradiol (E2), the synthetic estrogen diethylstilbestrol (DES), and the phytoestrogen genistein (GEN) are intimately associated with their binding to the estrogen receptor α ligand binding domain (ERα LBD) and accordingly allostery. Molecular modeling techniques have been performed on agonists in complex with the LBD, focusing on the pivotal role of His524 modeled as the ε-tautomer and the protonated form (depending on pH). It is found that E2 binds to the active LBD with the aid of Leu525, showing existing stable patterns of an H-binding network with Glu419 via His524 in all models. The main difference seen in the effect is that the full agonists E2 and DES have higher binding energies to the protonated His524 than the partial agonists GEN and Way-169916 (W), which is in line with noted experimental transcriptional activities. In conclusion, the study demonstrates that the phytoestrogen GEN interacts differently with the LBD than what E2 and DES do, which explains the observed signaling differences.

COMPARE analysis of the toxicity of an iminoquinone derivative of the imidazo[5,4-f]benzimidazoles with NAD(P)H:quinone oxidoreductase 1 (NQO1) activity and computational docking of quinones as NQO1 substrates

Synthesis and cytotoxicity of imidazo[5,4-f]benzimidazolequinones and iminoquinone derivatives is described, enabling structure-activity relationships to be obtained. The most promising compound (an iminoquinone derivative) has undergone National Cancer Institute (NCI) 60 cell line (single and five dose) screening, and using the NCI COMPARE program, has shown correlation to NQO1 activity and to other NQO1 substrates. Common structural features suggest that the iminoquinone moiety is significant with regard to NQO1 specificity. Computational docking into the active site of NQO1 was performed, and the first comprehensive mitomycin C (MMC)-NQO1 docking study is presented. Small distances for hydride reduction and high binding affinities are characteristic of MMC and of iminoquinones showing correlations with NQO1 via COMPARE analysis. Docking also indicated that the presence of a substituent capable of hydrogen bonding to the His194 residue is important in influencing the orientation of the substrate in the NQO1 active site, leading to more efficient reduction.

Formation mechanism and structure of a guanine-uracil DNA intrastrand cross-link

The formation and structure of the 5'-G[8-5]U-3' intrastrand cross-link are studied using density functional theory and molecular dynamics simulations due to the potential role of this lesion in the activity of 5-halouracils in antitumor therapies. Upon UV irradiation of 5-halouracil-containing DNA, a guanine radical cation reacts with the uracil radical to form the cross-link, which involves phosphorescence or an intersystem crossing and a rate-determining step of bond formation. Following ionizing radiation, guanine and the uracil radical react, with a rate-limiting step involving hydrogen atom removal. Although cross-link formation from UV radiation is favored, comparison of calculated reaction thermokinetics with that for related experimentally observed purine-pyrimidine cross-links suggests this lesion is also likely to form from ionizing radiation. For the first time, the structure of 5'-G[8-5]U-3' within DNA is identified by molecular dynamics simulations. Furthermore, three conformations of cross-linked DNA are revealed, which differ in the configuration of the complementary bases. Distortions, such as unwinding, are localized to the cross-linked dinucleotide and complementary nucleotides, with minimal changes to the flanking bases. Global changes to the helix, such as bending and groove alterations, parallel cisplatin-induced distortions, which indicate 5'-G[8-5]U-3', may contribute to the cytotoxicity of halouracils in tumor cell DNA using similar mechanisms.

Catalytic mechanism and roles of Arg197 and Thr183 in the Staphylococcus aureus sortase A enzyme

The sortase A enzyme, which catalyzes the peptidoglycan cell wall anchoring reaction of LPXTG surface proteins, has been proposed to be a universal target for therapeutic agents against Gram-positive bacteria. The catalytic mechanism of the Staphylococcus aureus sortase A enzyme has been systematically studied using molecular dynamics simulations, ONIOM(DFT:MM) calculations, and QM/MM charge deletion analysis. The catalytic roles of Arg197 and Thr183 were analyzed. Our calculations show that Arg197 has several important roles in the mechanism. It is crucial for substrate binding, and is capable of reversible shift of its hydrogen bonds between the LP and TG carbonyls of the LPXTG substrate motif, depending on the protonation state of the catalytic Cys184-His120 dyad. Arg197 stabilizes the catalytic dyad in the active ion pair form but at the same time raises the barrier to acylation by approximately 8 kcal/mol. Thr183 is also essential for the catalytic reaction in that it correspondingly lowers the barrier by the same amount via electrostatic interactions. The catalytic mechanism proceeds via proton transfer from His120, followed by nucleophilic attack from the thiolate anion of Cys184. The data thus supports the proposed reverse protonation mechanism, and disproves the hypothesis of the Arg197 generating an oxyanion hole to stabilize the tetrahedral intermediate of the reaction.

1,2-Bis(methylsulfonyl)-1-(2-chloroethyl)-2-[[1-(4-nitrophenyl)ethoxy]carbonyl]hydrazine (KS119): a cytotoxic prodrug with two stable conformations differing in biological and physical properties

The anticancer prodrug 1,2-bis(methylsulfonyl)-1-(2-chloroethyl)-2-[[1-(4-nitrophenyl)ethoxy]carbonyl]hydrazine (KS119) selectively releases a short-lived cytotoxin following enzymatic reduction in hypoxic environments found in solid tumors. KS119, in addition to two enantiomers, has two stable atropisomers (conformers differing in structure owing to hindered bond rotation) that interconvert at 37 °C in aqueous solution by first-order kinetics with t(1/2) values of ∼50 and ∼64 h. The atropisomers differ in physical properties such as partition coefficients that allow their chromatographic separation on non-chiral columns. A striking difference in the rate of metabolism of the two atropisomers occurs in intact EMT6 murine mammary carcinoma cells under oxygen-deficient conditions. A structurally related molecule, 1,2-bis(methylsulfonyl)-1-(2-chloroethyl)-2-[[1-(3-hydroxy-4-nitrophenyl)ethoxy]carbonyl]hydrazine (KS119WOH), was also found to exist in similar stable atropisomers. The ratio of the atropisomers of KS119 and structurally related agents has the potential to impact the bioavailability, activation, and therapeutic activity. Thus, thermally stable atropisomers/conformers in small molecules can result in chemically and enantiomerically pure compounds having differences in biological activities.

Computational evidence for the role of Arabidopsis thaliana UVR8 as UV-B photoreceptor and identification of its chromophore amino acids

A homology model of the Arabidopsis thaliana UV resistance locus 8 (UVR8) protein is presented herein, showing a seven-bladed β-propeller conformation similar to the globular structure of RCC1. The UVR8 amino acid sequence contains a very high amount of conserved tryptophans, and the homology model shows that seven of these tryptophans cluster at the 'top surface' of the UVR8 protein where they are intermixed with positive residues (mainly arginines) and a couple of tyrosines. Quantum chemical calculations of excitation spectra of both a large cluster model involving all twelve above-mentioned residues and smaller fragments thereof reveal that absorption maxima appearing in the 280-300 nm range for the full cluster result from interactions between the central tryptophans and surrounding arginines. This observation coincides with the published experimentally measured action spectrum for the UVR8-dependent UV-B stimulation of HY5 transcription in mature A. thaliana leaf tissue. In total these findings suggest that UVR8 has in fact in itself the ability to be an ultraviolet-B photoreceptor in plants.

The pea SAD short-chain dehydrogenase/reductase: quinone reduction, tissue distribution, and heterologous expression

The pea (Pisum sativum) tetrameric short-chain alcohol dehydrogenase-like protein (SAD) family consists of at least three highly similar members (SAD-A, -B, and -C). According to mRNA data, environmental stimuli induce SAD expression. The aim of this study was to characterize the SAD proteins by examining their catalytic function, distribution in pea, and induction in different tissues. In enzyme activity assays using a range of potential substrates, the SAD-C enzyme was shown to reduce one- or two-ring-membered quinones lacking long hydrophobic hydrocarbon tails. Immunological assays using a specific antiserum against the protein demonstrated that different tissues and cell types contain small amounts of SAD protein that was predominantly located within epidermal or subepidermal cells and around vascular tissue. Particularly high local concentrations were observed in the protoderm of the seed cotyledonary axis. Two bow-shaped rows of cells in the ovary and the placental surface facing the ovule also exhibited considerable SAD staining. Ultraviolet-B irradiation led to increased staining in epidermal and subepidermal cells of leaves and stems. The different localization patterns of SAD suggest functions both in development and in responses to environmental stimuli. Finally, the pea SAD-C promoter was shown to confer heterologous wound-induced expression in Arabidopsis (Arabidopsis thaliana), which confirmed that the inducibility of its expression is regulated at the transcriptional level.

Structural changes of Listeria monocytogenes Sortase A: a key to understanding the catalytic mechanism

Listeria monocytogenes is one of the most virulent foodborne pathogens. L. monocytogenes Sortase A (SrtA) enzyme, which catalyzes the cell wall anchoring reaction of the leucine, proline, X, threonine, and glycine proteins (LPXTG, where X is any amino acid), is a target for the development of antilisteriosis drugs. In this study, the structure of the L. monocytogenes SrtA enzyme-substrate complex was obtained using homology modeling, molecular docking and molecular dynamics simulations. Explicit enzyme-substrate interactions in the inactive and active forms of the enzyme were compared, based on 30 ns simulations on each system. The active site arginine (Arg 197) was found to be able change its hydrogen donor interactions from the LP backbone carbonyl groups of the LPXTG substrate in the inactive form, to the TG backbone carbonyls in the active form, which could be of importance for holding the substrate in position for the catalytic process.

The role of the pyridoxine (vitamin B6) biosynthesis enzyme PDX1 in ultraviolet-B radiation responses in plants

Ultraviolet-B radiation regulates plant growth and morphology at low and ambient fluence rates but can severely impact on plants at higher doses. Some plant UV-B responses are related to the formation of reactive oxygen species (ROS) and pyridoxine (vitamin B(6)) has been reported to be a quencher of ROS. UV-B irradiation of Arabidopsis Col-0 plants resulted in increased levels of PDX1 protein, compared with UV-A-exposed plants. This was shown by immunoblot analysis using specific polyclonal antibodies raised against the recombinant PDX1.3 protein and confirmed by mass spectrometry analysis of immunoprecipitated PDX1. The protein was located mainly in the cytosol but also to a small extent in the membrane fraction of plant leaves. Immunohistochemical analysis performed in pea revealed that PDX1 is present in UV-B-exposed leaf mesophyll and palisade parenchyma but not in epidermal cells. Pyridoxine production increased in Col-0 plants exposed to 3 days of UV-B, whereas in an Arabidopsis pdx1.3 mutant UV-B did not induce pyridoxine biosynthesis. In gene expression studies performed after UV-B exposure, the pdx1.3 mutant showed elevated transcript levels for the LHCB1*3 gene (encoding a chlorophyll a/b-binding protein of the photosystem II light-harvesting antenna complex) and the pathogenesis-related protein 5 (PR-5) gene, compared with wild type.