A concept of dual-responsive prodrugs based on oligomerization-controlled reactivity of ester groups: an improvement of cancer cells versus neutrophils selectivity of camptothecin.

Many known chemotherapeutic anticancer agents exhibit neutropenia as a dose-limiting side effect. In this paper we suggest a prodrug concept solving this problem for camptothecin (HO-cpt). The prodrug is programmed according to Boolean "AND" logic. In the absence of H2O2 (trigger T1), e.g. in the majority of normal cells, it exists as an inactive oligomer. In cancer cells and in primed neutrophils (high H2O2), the oligomer is disrupted forming intermediate (inactive) lipophilic cationic species. These are accumulated in mitochondria (Mit) of cancer cells, where they are activated by hydrolysis at mitochondrial pH 8 (trigger T2) with formation of camptothecin. In contrast, the intermediates remain stable in neutrophils lacking Mit and therefore a source of T2. In this paper we demonstrated a proof-of-concept. Our prodrug exhibits antitumor activity both in vitro and in vivo, but is not toxic to normal cell and neutrophils in contrast to known single trigger prodrugs and the parent drug HO-cpt.

In Silico Study of Camptothecin-Based Pro-Drugs Binding to Human Carboxylesterase 2.

Pro-drugs, which ideally release their active compound only at the site of action, i.e., in a cancer cell, are a promising approach towards an increased specificity and hence reduced side effects in chemotherapy. A popular form of pro-drugs is esters, which are activated upon their hydrolysis. Since carboxylesterases that catalyse such a hydrolysis reaction are also abundant in normal tissue, it is of great interest whether a putative pro-drug is a probable substrate of such an enzyme and hence bears the danger of being activated not just in the target environment, i.e., in cancer cells. In this work, we study the binding mode of carboxylesters of the drug molecule camptothecin, which is an inhibitor of topoisomerase I, of varying size to human carboxylesterase 2 (HCE2) by molecular docking and molecular dynamics simulations. A comparison to irinotecan, known to be a substrate of HCE2, shows that all three pro-drugs analysed in this work can bind to the HCE2 protein, but not in a pose that is well suited for subsequent hydrolysis. Our data suggest, moreover, that for the irinotecan substrate, a reactant-competent pose is stabilised once the initial proton transfer from the putative nucleophile Ser202 to the His431 of the catalytic triad has already occurred. Our simulation work also shows that it is important to go beyond the static models obtained from molecular docking and include the flexibility of enzyme-ligand complexes in solvents and at a finite temperature. Under such conditions, the pro-drugs studied in this work are unlikely to be hydrolysed by the HCE2 enzyme, indicating a low risk of undesired drug release in normal tissue.

Protonation-State Dependence of Hydration and Interactions in the Two Proton-Conducting Channels of Cytochrome c Oxidase.

Cytochrome c Oxidase (CcO), a membrane protein of the respiratory chain, pumps protons against an electrochemical gradient by using the energy of oxygen reduction to water. The ("chemical") protons required for this reaction and those pumped are taken up via two distinct channels, named D-channel and K-channel, in a step-wise and highly regulated fashion. In the reductive phase of the catalytic cycle, both channels transport protons so that the pumped proton passes the D-channel before the "chemical" proton has crossed the K-channel. By performing molecular dynamics simulations of CcO in the O→E redox state (after the arrival of the first reducing electron) with various combinations of protonation states of the D- and K-channels, we analysed the effect of protonation on the two channels. In agreement with previous work, the amount of water observed in the D-channel was significantly higher when the terminal residue E286 was not (yet) protonated than when the proton arrived at this end of the D-channel and E286 was neutral. Since a sufficient number of water molecules in the channel is necessary for proton transport, this can be understood as E286 facilitating its own protonation. K-channel hydration shows an even higher dependence on the location of the excess proton in the K-channel. Also in agreement with previous work, the K-channel exhibits a very low hydration level that likely hinders proton transfer when the excess proton is located in the lower part of the K-channel, that is, on the N-side of S365. Once the proton has passed S365 (towards the reaction site, the bi-nuclear centre (BNC)), the amount of water in the K-channel provides hydrogen-bond connectivity that renders proton transfer up to Y288 at the BNC feasible. No significant direct effect of the protonation state of one channel on the hydration level, hydrogen-bond connectivity, or interactions between protein residues in the other channel could be observed, rendering proton conductivity in the two channels independent of each other. Regulation of the order of proton uptake and proton passage in the two channels such that the "chemical" proton leaves its channel last must, therefore, be achieved by other means of communication, such as the location of the reducing electron.

Interplay of Hydration and Protonation Dynamics in the K-Channel of Cytochrome c Oxidase.

Cytochrome c oxidase is a membrane protein of the respiratory chain that consumes protons and molecular oxygen to produce water and uses the resulting energy to pump protons across the membrane. Our molecular dynamics simulations with an excess proton located at different positions in one of the proton-conducting channels, the K-channel, show a clear dependence of the number of water molecules inside the channel on the proton position. A higher hydration level facilitates the formation of hydrogen-bonded chains along which proton transfer can occur. However, a sufficiently high hydration level for such proton transport is observed only when the excess proton is located above S365, i.e., the lower third of the channel. From the channel entrance up to this point, proton transport is via water molecules as proton carriers. These hydronium ions move with their surrounding water molecules, up to K362, filling and widening the channel. The conformation of K362 depends on its own protonation state and on the hydration level, suggesting its role to be proton transport from a hydronium ion at the height of K362 to the upper part of the channel via a conformational change. The protonation-dependent conformational dynamics of E101 at the bottom of the channel renders proton transfer via E101 unlikely. Instead, its role is rather that of an amplifier of H96's proton affinity, suggesting H96 as the initial proton acceptor.

Effect of a U:G mispair on the water around DNA.

DNA repair proteins are able to discriminate DNA lesions among an abundance of intact DNA with high selectivity. To investigate detectable characteristics of one specific lesion, we compare statistical results from molecular dynamics simulations of two different DNA in water, one with an intact C:G pair and one that contains a U:G mispair, and perform a comparative analysis of the water dynamics around the two. Our data show that in addition to the local DNA conformation, also the surrounding water shell exhibits significant differences that may help mispair discrimination. The chemical groups which account for a U:G mispair to exhibit a wobble conformation instead of the 'proper' Watson-Crick pairing of a C:G pair, that is an oxygen atom (in uracil) instead of an amino group (in cytosine), also order the water molecules around the bases such that they act predominantly as hydrogen-bond donor or acceptor to the uracil or cytosine base, respectively. These changes in water conformation stretch into the second solvation shell, which may be exploited by repair enzymes to achieve lesion detection with high efficiency.

Oxidation Enhances Binding of Extrahelical 5-Methyl-Cytosines by Thymine DNA Glycosylase.

The DNA repair protein thymine DNA glycosylase (TDG) removes mispaired or damaged bases, such as oxidized methyl-cytosine, from DNA by cleavage of the glycosidic bond between the sugar and the target base flipped into the enzyme's active site. The enzyme is active against formyl-cytosine and carboxyl-cytosine, whereas the lower oxidized hydroxymethyl-cytosine and methyl-cytosine itself are not processed by the enzyme. Molecular dynamics simulations with thermodynamic integration of TDG complexed to DNA carrying one of four different (oxidized) methyl-cytosine bases in extrahelcial conformation, methyl-cytosine (mC), hydroxymethyl-cytosine (hmC), formyl-cytosine (fC), or carboxyl-cytosine (caC), show a more favorable binding affinity of the higher oxidized forms, fC and caC, than the nonsubstrate bases hmC and mC. Despite rather comparable, reaction-competent conformations of the flipped bases in the active site of the enzyme, more and stronger interactions with active site residues account for the preferred binding of the higher oxidized bases. Binding of the negatively charged caC and the neutral fC are strengthened by interactions with positively charged His151. Our calculated proton affinities find this protonation state of His151 the preferred one in the presence of caC and conceivable in the presence of fC as well as increasing the binding affinity toward the two bases. Discrimination of the substrate bases is further achieved by the backbone of Tyr152 that forms a strong hydrogen bond to the carboxyl and formyl oxygen atoms of caC and fC, respectively, a contact that is completely lacking in mC and much weaker in hmC. Overall, our computational results indicate that the enzyme discriminates the different oxidation forms of methyl-cytosine already at the formation of the extrahelical complexes.

Structural, spectroscopic, in silico, in vitro and DNA binding evaluations of tyrosyl-lysyl-threonine.

The main objective of the present study is to investigate the molecular structure and DNA binding interaction of the tyrosyl-lysyl-threonine (YKT) tripeptide, which has anticancer, antioxidant and analgesic properties, using various in silico (MD, QM, molecular docking), spectroscopic (UV, FT-IR, FTIR-ATR, Raman, gel electrophoresis) and in vitro (MCF-7 and HeLa cancer cell lines and BEAS-2B cell line) methods. The optimized geometry, vibrational wavenumbers, molecular electrostatic potential (MEP), natural bond orbital (NBO) and HOMO-LUMO (highest occupied molecular orbital- lowest unoccupied molecular orbital) calculations were carried out with Density Functional Theory (DFT) using B3LYP/6-311++G(d,p) basis set to indicate conformational, vibrational and intramolecular charge transfer characteristics. The assignment of all fundamental theoretical vibration wavenumbers was performed using potential energy distribution analysis (PED). DNA is a significant pharmacological target of drugs in several diseases such as cancer. For this reason, molecular docking calculation was used to elucidate the binding and interaction between YKT tripeptide and DNA at the atomic level. Also, the dynamic behaviors of YKT and DNA was examined using MD simulations. Besides, the interaction of YKT with DNA was experimentally examined by UV titration method and agarose gel electrophoresis method. Experimental results showed that YKT was intercalatively and electrostatically bound to CT-DNA (Calf thymus DNA) and cleavage pBR322 DNA in the presence of H2O2. The pharmacokinetic profile of YKT was also obtained. Cytotoxic effect of YKT was evaluated on MCF-7, HeLa and BEAS-2B cell lines. Hence, these studies about YKT tripeptide may pave the way for the development of various cancer drugs. Communicated by Ramaswamy H. Sarma.

Interaction of Thymine DNA Glycosylase with Oxidised 5-Methyl-cytosines in Their Amino- and Imino-Forms.

Thymine DNA Glycosylase (TDG) is an enzyme of the base excision repair mechanism and removes damaged or mispaired bases from DNA via hydrolysis of the glycosidic bond. Specificity is of high importance for such a glycosylase, so as to avoid the damage of intact DNA. Among the substrates reported for TDG are mispaired uracil and thymine but also formyl-cytosine and carboxyl-cytosine. Methyl-cytosine and hydroxylmethyl-cytosine are, in contrast, not processed by the TDG enzyme. We have in this work employed molecular dynamics simulations to explore the conformational dynamics of DNA carrying a formyl-cytosine or carboxyl-cytosine and compared those to DNA with the non-cognate bases methyl-cytosine and hydroxylmethyl-cytosine, as amino and imino tautomers. Whereas for the mispairs a wobble conformation is likely decisive for recognition, all amino tautomers of formyl-cytosine and carboxyl-cytosine exhibit the same Watson-Crick conformation, but all imino tautomers indeed form wobble pairs. The conformational dynamics of the amino tautomers in free DNA do not exhibit differences that could be exploited for recognition, and also complexation to the TDG enzyme does not induce any alteration that would indicate preferable binding to one or the other oxidised methyl-cytosine. The imino tautomers, in contrast, undergo a shift in the equilibrium between a closed and a more open, partially flipped state, towards the more open form upon complexation to the TDG enzyme. This stabilisation of the more open conformation is most pronounced for the non-cognate bases methyl-cytosine and hydroxyl-cytosine and is thus not a likely mode for recognition. Moreover, calculated binding affinities for the different forms indicate the imino forms to be less likely in the complexed DNA. These findings, together with the low probability of imino tautomers in free DNA and the indifference of the complexed amino tautomers, suggest that discrimination of the oxidised methyl-cytosines does not take place in the initial complex formation.

Long-range allostery mediates cooperative adenine nucleotide binding by the Ski2-like RNA helicase Brr2.

Brr2 is an essential Ski2-like RNA helicase that exhibits a unique structure among the spliceosomal helicases. Brr2 harbors a catalytically active N-terminal helicase cassette and a structurally similar but enzymatically inactive C-terminal helicase cassette connected by a linker region. Both cassettes contain a nucleotide-binding pocket, but it is unclear whether nucleotide binding in these two pockets is related. Here we use biophysical and computational methods to delineate the functional connectivity between the cassettes and determine whether occupancy of one nucleotide-binding site may influence nucleotide binding at the other cassette. Our results show that Brr2 exhibits high specificity for adenine nucleotides, with both cassettes binding ADP tighter than ATP. Adenine nucleotide affinity for the inactive C-terminal cassette is more than two orders of magnitude higher than that of the active N-terminal cassette, as determined by slow nucleotide release. Mutations at the intercassette surfaces and in the connecting linker diminish the affinity of adenine nucleotides for both cassettes. Moreover, we found that abrogation of nucleotide binding at the C-terminal cassette reduces nucleotide binding at the N-terminal cassette 70 Å away. Molecular dynamics simulations identified structural communication lines that likely mediate these long-range allosteric effects, predominantly across the intercassette interface. Together, our results reveal intricate networks of intramolecular interactions in the complex Brr2 RNA helicase, which fine-tune its nucleotide affinities and which could be exploited to regulate enzymatic activity during splicing.

Effect of the Hydration Shell on the Carbonyl Vibration in the Ala-Leu-Ala-Leu Peptide.

The vibrational spectrum of the Ala-Leu-Ala-Leu peptide in solution, computed from first-principles simulations, shows a prominent band in the amide I region that is assigned to stretching of carbonyl groups. Close inspection reveals combined but slightly different contributions by the three carbonyl groups of the peptide. The shift in their exact vibrational signature is in agreement with the different probabilities of these groups to form hydrogen bonds with the solvent. The central carbonyl group has a hydrogen bond probability intermediate to the other two groups due to interchanges between different hydrogen-bonded states. Analysis of the interaction energies of individual water molecules with that group shows that shifts in its frequency are directly related to the interactions with the water molecules in the first hydration shell. The interaction strength is well correlated with the hydrogen bond distance and hydrogen bond angle, though there is no perfect match, allowing geometrical criteria for hydrogen bonds to be used as long as the sampling is sufficient to consider averages. The hydrogen bond state of a carbonyl group can therefore serve as an indicator of the solvent's effect on the vibrational frequency.

Tacrine-sugar mimetic conjugates as enhanced cholinesterase inhibitors.

We have used the Cu(i)-catalyzed azide-alkyne Huisgen cycloaddition reaction to obtain two families of bivalent heterodimers where tacrine is connected to an azasugar or iminosugar, respectively, via linkers of variable length. The heterodimers were investigated as cholinesterase inhibitors and it was found that their activity increased with the length of the linker. Two of the heterodimers were significantly stronger acetylcholinesterase inhibitors than the monomeric tacrine. Molecular modelling indicated that the longer heterodimers fitted better into the active gorge of acetylcholinesterase than the shorter counterparts and the former provided more efficient simultaneous interaction with the tryptophan residues in the catalytic anionic binding site (CAS) and the peripheral anionic binding site (PAS).

Quantifying how step-wise fluorination tunes local solute hydrophobicity, hydration shell thermodynamics and the quantum mechanical contributions of solute-water interactions.

The ability to locally tune solute-water interactions and thus control the hydrophilic/hydrophobic character of a solute is key to control molecular self-assembly and to develop new drugs and biocatalysts; it has been a holy grail in synthetic chemistry and biology. To date, the connection between (i) the hydrophobicity of a functional group; (ii) the local structure and thermodynamics of its hydration shell; and (iii) the relative influence of van der Waals (dispersion) and electrostatic interactions on hydration remains unclear. We investigate this connection using spectroscopic, classical simulation and ab initio methods by following the transition from hydrophile to hydrophobe induced by the step-wise fluorination of methyl groups. Along the transition, we find that water-solute hydrogen bonds are progressively transformed into dangling hydroxy groups. Each structure has a distinct thermodynamic, spectroscopic and quantum-mechanical signature connected to the associated local solute hydrophobicity and correlating with the relative contribution of electrostatics and dispersion to the solute-water interactions.

Proton transfer in the D-channel of cytochrome c oxidase modeled by a transition network approach.

BACKGROUND: Determination of proton uptake pathways in Cytochrome c Oxidase is difficult due to the complexity of the system. The transition networks approach allows sampling of proton transfer pathways without predefined reaction coordinate. METHODS: Computation of the proton transfer pathways in a model of the D-channel of cytochrome c oxidase has been performed by a transition network approach that combines discrete, optimisation based and molecular dynamics based sampling. RESULTS: The optimal pathway involves an opening of the so-called asparagine gate, hydration of the asparagine region, the formation of a hydrogen-bonded chain, and finally concerted proton hole transport along this chain. The optimal pathway finds the protonation of residue H26 close to the channel entrance favourable for lowering the transition energies of subsequent steps, in particular, opening of the Asn gate and formation of a hydrogen-bonded chain. Residue Y33 plays an important role in shuttling the transferred proton hole. CONCLUSIONS: The optimal pathway found by the transition network approach shows the same important characteristics as pathways determined earlier by other methods. The computed barrier and reaction energies are also in good agreement with previous studies. The transition network approach provides an alternative to explore pathways in complex systems. GENERAL SIGNIFICANCE: The correct function of the enzyme as oxidase and proton pump depends on the interplay of several redox and proton transport steps. Understanding the proton transport mechanism is therefore key to understanding the protein's function. The complex nature of long- distances proton transfer through a protein requires a non-trivial simulation strategy.

Molecular Dynamics Simulations of a Chimeric Androgen Receptor Protein (SPARKI) Confirm the Importance of the Dimerization Domain on DNA Binding Specificity.

The DNA binding domains of Androgen/Glucocorticoid receptors (AR/GR), members of class I steroid receptors, bind as a homo-dimer to a cis-regulatory element. These response elements are arranged as inverted repeat (IR) of hexamer "AGAACA", separated with a 3 base pairs spacer. DNA binding domains of the Androgen receptor, AR-DBDs, in addition, selectively recognize a direct-like repeat (DR) arrangement of this hexamer. A chimeric AR protein, termed SPARKI, in which the second zinc-binding motif of AR is swapped with that of GR, however, fails to recognize DR-like elements. By molecular dynamic simulations, we identify how the DNA binding domains of the wild type AR/GR, and also the chimeric SPARKI model, distinctly interact with both IR and DR response elements. AR binds more strongly to DR than GR binds to IR elements. A SPARKI model built from the structure of the AR (SPARKI-AR) shows significantly fewer hydrogen bond interactions in complex with a DR sequence than with an IR sequence. Moreover, a SPARKI model based on the structure of the GR (SPARKI-GR) shows a considerable distortion in its dimerization domain when complexed to a DR-DNA whereas it remains in a stable conformation in a complex with an IR-DNA. The diminished interaction of SPARKI-AR with and the instability of SPARKI-GR on DR response elements agree with SPARKI's lack of affinity for these sequences. The more GR-like binding specificity of the chimeric SPARKI protein is further emphasized by both SPARKI models binding even more strongly to IR elements than observed for the DNA binding domain of the GR.

A Conserved Kinase-Based Body-Temperature Sensor Globally Controls Alternative Splicing and Gene Expression.

Homeothermic organisms maintain their core body temperature in a narrow, tightly controlled range. Whether and how subtle circadian oscillations or disease-associated changes in core body temperature are sensed and integrated in gene expression programs remain elusive. Furthermore, a thermo-sensor capable of sensing the small temperature differentials leading to temperature-dependent sex determination (TSD) in poikilothermic reptiles has not been identified. Here, we show that the activity of CDC-like kinases (CLKs) is highly responsive to physiological temperature changes, which is conferred by structural rearrangements within the kinase activation segment. Lower body temperature activates CLKs resulting in strongly increased phosphorylation of SR proteins in vitro and in vivo. This globally controls temperature-dependent alternative splicing and gene expression, with wide implications in circadian, tissue-specific, and disease-associated settings. This temperature sensor is conserved across evolution and adapted to growth temperatures of diverse poikilotherms. The dynamic temperature range of reptilian CLK homologs suggests a role in TSD.

Hydrophobic but Water-Friendly: Favorable Water-Perfluoromethyl Interactions Promote Hydration Shell Defects.

Although perfluorination is known to enhance hydrophobicity and change protein activity, its influence on hydration-shell structure and thermodynamics remains an open question. Here we address that question by combining experimental Raman multivariate curve resolution spectroscopy with theoretical classical simulations and quantum mechanical calculations. Perfluorination of the terminal methyl group of ethanol is found to enhance the disruption of its hydration-shell hydrogen bond network. Our results reveal that this disruption is not due to the associated volume change but rather to the electrostatic stabilization of the water dangling OH···F interaction. Thus, the hydration shell structure of fluorinated methyl groups results from a delicate balance of solute-water interactions that is intrinsically different from that associated with a methyl group.

Quantum Mechanical/Molecular Mechanical Analysis of the Catalytic Mechanism of Phosphoserine Phosphatase.

Phosphoserine phosphatase (PSP), a member of the haloacid dehalogenase (HAD) superfamily that comprises the vast majority of phosphotransferases, is likely a steady-state regulator of the level of d-serine in the brain. The proposed catalytic cycle of PSP consists of a two-step mechanism: formation of a phospho-enzyme intermediate by phosphate transfer to Asp11 and its subsequent hydrolysis. Our combined quantum mechanical/molecular mechanical (QM/MM) calculations of the reaction pathways favour a dissociative mechanism of nucleophilic substitution via a trigonal-planar metaphosphate-like configuration for both steps, associated with proton transfer to the leaving group or from the nucleophile. This proton transfer is facilitated by active site residue Asp13 that acts as both a general base and a general acid. Free energy calculation on the reaction pathways further support the structural role of the enzymatic environment and the active site architecture. The choice of a proper reaction coordinate along which to bias the free energy calculations can be guided by a projection of the canonical reaction coordinate obtained from a chain-of-state optimisation onto important internal coordinates.

Prediction of perturbed proton transfer networks.

The transfer of protons through proton translocating channels is a complex process, for which direct samplings of different protonation states and side chain conformations in a transition network calculation provide an efficient, bias-free description. In principle, a new transition network calculation is required for every unsampled change in the system of interest, e.g. an unsampled protonation state change, which is associated with significant computational costs. Transition networks void of or including an unsampled change are termed unperturbed or perturbed, respectively. Here, we present a prediction method, which is based on an extensive coarse-graining of the underlying transition networks to speed up the calculations. It uses the minimum spanning tree and a corresponding sensitivity analysis of an unperturbed transition network as initial guess and refinement parameter for the determination of an unknown, perturbed transition network. Thereby, the minimum spanning tree defines a sub-network connecting all nodes without cycles and minimal edge weight sum, while the sensitivity analysis analyzes the stability of the minimum spanning tree towards individual edge weight reductions. Using the prediction method, we are able to reduce the calculation costs in a model system by up to 80%, while important network properties are maintained in most predictions.

Interactions of the DNA Repair Enzyme Human Thymine DNA Glycosylase with Cognate and Noncognate DNA.

Glycosylases specifically recognize and flip their target base out of the DNA helix into the enzyme's active site. Our simulations show that a partially flipped state, already present in free DNA carrying a T:G mispair, becomes the more probable state compared to the closed state after binding of thymine DNA glycosylase (TDG). Paired thymine (T:A) or methyl-cytosine (mC:G) does not exhibit a partially flipped state in free or complexed DNA. Important enzyme-DNA interactions exhibit significant strength in the intrahelical and extrahelical TDG-DNA complexes. The computed binding free energy differences suggest these interactions account for the stabilization of the partially flipped state, thereby driving the T:G mispair toward base flip. In the fully flipped state, the cognate base thymine is significantly better accommodated in the enzyme's active site than noncognate bases are, suggesting the hydrolysis step as the last of several stages at which base recognition can be achieved.